Next Article in Journal
Identification of Human Norovirus GII.3 Blockade Antibody Epitopes
Previous Article in Journal
RETRACTED: SARS–CoV–2 Spike Impairs DNA Damage Repair and Inhibits V(D)J Recombination In Vitro
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 10
10.3390/v13102057
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 and Pathogenicity of H7N7 and H7N9 Avian Influenza Viruses Isolated from South Korea

by
Eun-Jee Na
,
Young-Sik Kim
,
Yoon-Ji Kim
,
Jun-Soo Park
and
Jae-Ku Oem
*
Laboratory of Veterinary Infectious Diseases, College of Veterinary Medicine, Jeonbuk National University, Iksan 54596, Korea
*
Author to whom correspondence should be addressed.
Viruses 2021, 13(10), 2057; https://doi.org/10.3390/v13102057
Submission received: 25 September 2021 / Revised: 8 October 2021 / Accepted: 9 October 2021 / Published: 13 October 2021
(This article belongs to the Topic Veterinary Infectious Diseases)
Browse Figures
Review Reports Versions Notes

Abstract

:
H7 low pathogenic avian influenza viruses (LPAIVs) can mutate into highly pathogenic avian influenza viruses (HPAIVs). In addition to avian species, H7 avian influenza viruses (AIVs) also infect humans. In this study, two AIVs, H7N9 (20X-20) and H7N7 (34X-2), isolated from the feces of wild birds in South Korea in 2021, were genetically analyzed. The HA cleavage site of the two H7 Korean viruses was confirmed to be ELPKGR/GLF, indicating they are LPAIVs. There were no amino acid substitutions at the receptor-binding site of the HA gene of two H7 Korean viruses compared to that of A/Anhui/1/2013 (H7N9), which prefer human receptors. In the phylogenetic tree analysis, the HA gene of the two H7 Korean viruses shared the highest nucleotide similarity with the Korean H7 subtype AIVs. In addition, the HA gene of the two H7 Korean viruses showed high nucleotide similarity to that of the A/Jiangsu/1/2018(H7N4) virus, which is a human influenza virus originating from avian influenza virus. Most internal genes (PB2, PB1, PA, NP, NA, M, and NS) of the two H7 Korean viruses belonged to the Eurasian lineage, except for the M gene of 34X-2. This result suggests that active reassortment occurred among AIVs. In pathogenicity studies of mice, the two H7 Korean viruses replicated in the lungs of mice. In addition, the body weight of mice infected with 34X-2 decreased 7 days post-infection (dpi) and inflammation was observed in the peribronchiolar and perivascular regions of the lungs of mice. These results suggest that mammals can be infected with the two H7 Korean AIVs. Our data showed that even low pathogenic H7 AIVs may infect mammals, including humans, as confirmed by the A/Jiangsu/1/2018(H7N4) virus. Therefore, continuous monitoring and pathogenicity assessment of AIVs, even of LPAIVs, are required.

1. Introduction

Avian influenza viruses (AIVs) are eight segmented, single-stranded negative-sense RNA viruses [1]. Based on viral surface proteins, hemagglutinin (HA) and neuraminidase (NA), influenza viruses are classified into 18 HA and 11 NA subtypes. H1-16 and N1-9 have been detected in avian species, but H17-18 and N10-11 have been discovered only in bats [2,3,4]. In particular, H5 and H7 are important subtypes because they have the potential to mutate into highly pathogenic avian influenza viruses (HPAIVs), causing severe clinical signs in poultry [5,6,7]. According to previous research, the H7 low pathogenic avian influenza virus (LPAIV) is a precursor to the H7 HPAIV [8,9,10,11,12,13,14,15].
H7 LPAIVs have been detected in poultry farms worldwide, and wild birds are regarded as the origin of H7 LPAIVs [16]. In wild birds, all nine neuraminidase (N1–N9) are found in H7 AIVs, and H7N7 viruses are the subtypes reported in the largest number of countries geographically [17]. In South Korea, two LPAIVs, H7N3 and H7N8, were detected in domestic ducks in 2007 [18]. Since then, through active national surveillance, H7 LPAIVs have been isolated continuously from poultry farms and wild birds, and H7N7 viruses are the most common subtypes in wild birds [19,20]. Viruses isolated from poultry are closely related to wild bird isolates [19]. Moreover, some H7 LPAIVs isolated from wild birds showed pathogenicity in chickens in laboratory experiments.
AIVs can cause cross-species transmission from birds to mammals [21]. Avian-like H7N7 viruses have been detected in equines and pinnipeds [22]. In addition, there are many cases of H7 AIV infection in humans. The first case of H7N7 LPAIV directly transmitted from avian to human was reported in 1996 [23]. The infected woman had contact with a duck; she suffered from conjunctivitis and recovered naturally. In 2003, 89 cases of H7N7 HPAIV human infections occurred in the Netherlands [24]. One of the 89 patients died of acute respiratory distress syndrome, and three were confirmed to have contracted the virus through transmission between humans. A recent human infection of H7N7 AIVs was observed in three poultry workers in Italy in 2013, and they showed only conjunctivitis but no respiratory syndromes [25].
In 2013, the first case of human infection with H7N9 AIVs was reported in Shanghai, and since then, the H7N9 virus has shown five epidemic patterns in China [26,27]. The H7N9 viruses that occurred in 2013 were low pathogenic viruses but gradually mutated, resulting in high pathogenic strains in 2017 [6,28]. By 2019, 1568 cases of H7N9 human infections were confirmed worldwide, of which 616 deaths occurred [29]. However, after the introduction of H7N9 vaccine in poultry in 2017 in China, human infection with H7N9 virus has decreased dramatically [30,31]. Since September 2019, no human infection has been reported according to the FAO H7N9 situation update. Moreover, human infections related to other subtypes of H7 LPAIVs, such as H7N2 and H7N3, have often been reported since the 2000s [32]. Severe human infection with H7N4 originating from backyard poultry was reported in Jiangsu in 2018 [33].
In 2021, two AIVs, H7N7 and H7N9, were isolated from wild birds in South Korea. The molecular and phylogenetic characterizations of the two H7 Korean AIVs were analyzed in this study. In addition, pathogenicity in mammals was evaluated using a murine model.

2. Materials and Methods

2.1. Sampling

Fresh fecal samples (2499) of wild birds were collected between October 2020 and March 2021 by active monitoring of wild bird habitats. The collected samples were immediately shipped to the laboratory and stored between 2 and 8 °C until further analysis.

2.2. Virus Isolation

Each sample was suspended in phosphate-buffered saline (PBS, pH 7.4) supplemented with 100 U/μL of penicillin and 100 U/μL of streptomycin and centrifuged at 2800× g for 10 min. Supernatants were filtered with a 0.45-μm pore size syringe filter (GVS, Sanford, ME, USA). The filtered samples were inoculated into the allantoic cavities of 9–11-day-old specific-pathogen-free embryonated chicken eggs (Seong-Min Inc., Hwaseong, Korea) and incubated at 37 °C for 72 h. After chilling at 4 °C, the allantoic fluids were harvested and checked for AIV using hemagglutinin assay (HA assay) according to the World Organization for Animal Health (OIE) standard [34]. The AIV was verified in the HA-positive samples using the M gene quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) kit (iNtRON Biotechnology, Seongnam, Korea).

2.3. Subtyping and Species Identification

Viral RNA was extracted using a Miracle-AutoXT Automated Nucleic Acid Extraction System (iNtRON Biotechnology, Seongnam, Korea). To identify subtypes, qRT-PCR was conducted using the TOPreal™ One-step RT qPCR Kit (Enzynomics, Daejeon, Korea) with influenza-specific primers as previously described [35]. Species identification was performed using mitochondrial cytochrome C oxidase I mitochondrial gene in the DNA of fecal samples, as previously described [36].

2.4. Sequencing and Phylogenetic Analysis

Extracted viral RNAs were reverse transcribed to complementary DNAs (cDNAs) using a cDNA Synthesis kit (Wizbiosolutions, Seongnam, Korea) according to the manufacturer’s instructions. The synthesized cDNA was amplified using PCR and universal primers [37]. The amplicons were sequenced using an ABI 3730xl DNA analyzer (Applied Biosystems, Foster City, CA, USA). Nucleotide sequences were aligned using the ClustalW algorithm in BioEdit v.7.0.9.0. Phylogenetic trees were constructed using the maximum-likelihood (ML) method with 1000 bootstrap replications in MEGA (version 7.0) software [38]. The best substitution model for each segment was estimated using jModelTest (version 2.1.10) software [39,40]. The best models were the general time reverse model with invariable sites (GTR + I) for HA; GTR with gamma distribution (GTR + G) for PB2, PA, and NP; GTR + G + I for PB1 and M; Hasegawa Kishino Yano model with invariable sites (HKY + I) for NA7; and HKY + G for NA9 and NS. The reference sequences were retrieved from the National Center for Biotechnology Information Influenza Virus Resource (http://www.ncbi.nlm.nih.gov/genomes/FLU/FLU.html, accessed on 15 March 2021) and the Global Initiative on Sharing Avian Influenza Data (GISAID; http://www.gisaid.org, accessed on 15 March 2021). The nucleotide sequences of two H7 Korean AIVs were compared with sequences of NCBI-Influenza Virus Resource and GISAID database and reference sequences with high nucleotide identity were downloaded.
Maximum clade credibility (MCC) trees of the HA gene were constructed using Bayesian Evolution Analysis Sampling Trees (BEAST) version 1.10.4, to estimate the time to the most recent common ancestor [41]. We constructed a phylogenetic tree using H7 gene sequences obtained from GISAID between 2016 and 2021 in Asia. We used the HKY nucleotide substitution model with a gamma distribution and uncorrelated lognormal distribution relaxed clock method. All chains were run in 100,000,000 generations and the effective sample size values were greater than 200 in Tracer 1.7.1 (http://tree.bio.ed.ac.uk/software/tracer/, accessed on 27 April 2021). The MCC tree was constructed using TreeAnnotator v1.8.1 (http://beast.bio.ed.ac.uk/TreeAnnotator/, accessed on 6 May 2021) with 10% burn-in cutoffs and visualized with FigTree 1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/, accessed on 6 May 2021).

2.5. Pathogenicity in Mice

Forty female 6-week-old BALB/c mice (Samtaco, Osan, Korea) were used in this study. The mice were randomly divided into three groups, with 15 mice in the two experimental groups and 10 mice in the control group. After being anesthetized with avertin (Sigma-Aldrich, Milan, Italy), each group of 15 mice was inoculated intranasally with 107 egg infectivity dose (EID50) in a volume of 50 μL of 20X-20 (H7N9) and 34X-2 (H7N7) strains. Ten BALB/c mice were inoculated intranasally with 50 μL of PBS as the negative control group. The mice were monitored daily for body weight, body temperature, and clinical signs for 14 days post-infection (dpi). Three mice from the experimental group and two from the control group were sacrificed at 1, 3, 5, 7, and 14 dpi. The right lungs of mice were collected and were analyzed for virus detection using the M gene qRT-PCR kit (iNtRON Biotechnology, Seongnam, Korea), and the positive lung tissue as detected by qRT-PCR was used for virus titration using EID50 via the Reed and Muench method [42]. The left lobe of the mice was collected for pathology study and fixed in formalin for hematoxylin and eosin staining. The paraffin-embedded lung tissue was sectioned at 3 μm. For seroconversion and antigenic analysis, serum samples were collected from three mice in the experimental group and two mice from the control group at 14 dpi. The collected serum was treated with a receptor-destroying enzyme (Denka Seiken, Tokyo, Japan) and incubated for 18 h in a 37 °C water bath. Serum was used for hemagglutination inhibition (HI) tests following the OIE standard [34].

2.6. Statistical Analysis

The data and graphs shown in this study were analyzed using analysis of variance (ANOVA) in GraphPad Prism (version 6.0; GraphPad Software Inc., San Diego, CA, USA). The values shown are the mean with standard deviation (SD) or standard error of the mean (SEM). Statistical significance was set at p < 0.05.

2.7. Ethical Approval

All experiments were conducted at the animal biosafety level 2+ (ABSL + 2) facility. This study was approved by the Jeonbuk National University Laboratory Animal Research Center (JBNU 2021-0129).

3. Results

3.1. Identification of H7NX AIV Isolates

In January 2021, H7N9 and H7N7 AIVs were isolated from JeonJucheon (35°52′9.00′′ N, 127°6′20.28′′ E) and Dongjingang (35°41′22.53′′ N, 126°50′50.19′′ E), respectively. Cytochrome C oxidase I gene sequencing revealed that fecal samples came from Anas platyrhynchos (mallards). The H7N9 and N7N7 AIVs were designated as A/mallard/South Korea/JB20X-20/2021(H7N9; abbreviated as 20X-20) and A/mallard/South Korea/JB34X-2/2021(H7N7; abbreviated as 34X-2), respectively.

3.2. Molecular Characterization of H7 AIV Isolates

The whole genomes of 20X-20 and 34X-2 viruses were obtained and deposited in GenBank with the accession numbers MZ803114-MZ803121 and MZ803125-MZ803132, respectively. The obtained full gene sequences of 20X-20 and 34X-2 were compared with four selected reference H7 viruses: A/chicken/Jiangsu/1/2018 (H7N4; chicken/1), A/Jiangsu/1/2018 (H7N4; Jiangsu/1), A/Anhui/1/2013 (H7N9; Anhui/1), and A/Italy/3/2013 (H7N7; Italy/3) (Table 1) (14, 15, 19). The Jiangsu/1, Anhui/1, and Italy/3 viruses originate from the H7 AIV and cause human infections. The chicken/1 virus was isolated from chickens in the backyard of humans infected with the Jiangsu/1 virus. The amino acid sequences of the HA cleavage sites of 20X-20 and 34X-2 were ELPKGR/GLF, indicating they are LPAIVs. The amino acids at positions 186, 190, 225, 226, 227, and 228 of the HA gene (numbering based on H3) are related to host receptor binding efficiency [43,44,45]. Except for 186V and 226L in Anhui/1, all H7 viruses had 186G, 190E, 225G, 226Q, 227S, and 228G, indicating that they prefer avian receptors over human receptors. NA stalk deletions (69–73) were not detected in any of the H7 viruses except for the Anhui/1 virus. However, drug resistance-associated mutations were identified at the 117 position in the NA genes of 20X-20 and 34X-2.
The amino acid mutations of the two H7 Korean AIVs were mostly similar, except for K142R in the PA gene of 20X-20 and M105V in the NP gene of 34X-2. The PA gene of the two H7 Korean AIVs was clustered in the same subgroup in the phylogenetic tree but showed a different amino acid mutation at the 142 position. The amino acid substitutions of E627K in the PB2 gene, known to increase virulence in mice and human adaptation markers, were observed in Jiangsu/1 but not in the two H7 Korean AIVs [46,47]. Additionally, the amino acid substitutions of I292V of the PB2 gene, which increase the polymerase activity in both and avian cells were not identified in the two H7 Korean AIVs. The amino acid substitutions of N66S in the PB1-F2 gene, which increased virulence in mice, were observed in the two H7 Korean AIVs compared to the chicken/1, which most recently caused human infections [48]. The other mutations related to the virulence of the virus and drug resistance are shown in Table 1.

3.3. Phylogenetic Analysis of H7 AIVs Isolates

An ML tree analysis of the HA gene showed that the HA genes of 20X-20 and 34X-2 belonged to the Eurasian lineage (Figure 1). The HA genes of 20X-20 and 34X-2 showed 98.5% nucleotide identities with each other and were most closely related to A/wild duck/South Korea/KNU18-106/2018(H7N7), with identities of 97.8% and 98.5%, respectively (Table 2). Notably, the two H7 Korean AIVs showed high nucleotide similarity to the A/Jiangsu/1/2018(H7N4) virus (97.5% and 97%, respectively), which are human influenza viruses originating from AIVs. However, the HA genes of 20X-20 and 34X-2 viruses showed low nucleotide similarity to the HA gene of A/Italy/3/2013(H7N7) (90.2% each) and A/Anhui/1/2013(H7N9) (88.6% and 88.9%, respectively), which cause human infections.
To estimate the times of origins of the HA gene of the two H7 Korean AIVs, we reconstructed the MCC tree. As the HA gene of the two H7 Korean AIVs belonged to the Asian lineage, we reconstructed the MCC tree using the H7 gene of LPAIVs isolated between 2016 and 2021 in Asia. In the MCC tree analysis, HA genes formed two distinct genetic subgroups (G1 and G2) (Figure 2). The H7 LPAIV in the G1 subgroup was identified as the one circulating in South Korea in 2017, and the H7 LPAIVs in the G2 subgroup were closely related to the H7 LPAIV isolated in China between 2016 and 2017. The HA genes of 20X-20 and 34X-2 belonged to the G1 subgroup.
The results of the phylogenetic tree of the NA gene showed that the NA gene of the two H7 Korean AIVs belonged to the Eurasian lineage (Supplementary Materials Figure S1). The NA genes of 20X-20 and 34X-2 showed the highest nucleotide identities with the NA gene of A/Anas platyrhynchos/South Korea/JB31/96/2019(H11N9) (98.2%) and A/mallard/Korea/A15/2016(H7N7) (97.4%), respectively. However, the NA gene of the two H7 Korean AIVs clustered distinct subgroups with human influenza viruses.
Most of the internal genes (PB2, PB1, PA, NP, M, and NS) of the two H7 Korean AIVs belonged to the Eurasian lineage, except for the M gene of 34X-2 (Supplementary Materials Figure S2). The PB2, PB1, PA, NP, M, and NS genes of 20X-20 and 34X-2 shared 95.1%. 95%, 96.7%, 94.1%, 93.9%, and 99.5% nucleotide identities, respectively. The NP and M genes, which showed less than 95% nucleotide similarity, were clustered into different subgroups. The PA gene of 20X-20 was clustered into the same subgroup as the PA gene of A/Mandarin duck/Korea/WB246/2016(H5N6), which is an internal gene of HPAIV.

3.4. Pathogenicity in Mice

To determine the pathogenicity of 20X-20 and 34X-2 viruses in a murine model, 6-week-old BALB/c mice were intranasally inoculated with each virus. All mice were monitored for 14 days post-infection (dpi). Although mice infected with 20X-20 showed no clinical symptoms, mice infected with the 34X-2 showed ruffle fur and significantly lower body weight than mice infected with 20X-20 and PBS on 7 dpi (p < 0.01) (Figure 3). The mortality rate was 0% in all mice.
To observe the viral load in mouse lungs, the mice were sacrificed on 1, 3, 5, 7, and 14 dpi. Until 7 dpi, the virus antigens were detected in the lungs of mice, but there was no viral detection in the lungs of all groups on 14 dpi (Table 3). Further, 20X-20 replicated in the lungs with the highest titer on 5 dpi (2.80 ± 2.43 EID50/mL), which decreased dramatically on 7 dpi (0.48 ± 0.83 EID50/mL). However, 34X-2 showed a higher titer than 20X-20 in the lungs on 1, 3, 5, and 7 dpi. Moreover, 34X-2 replicated in the lung with the highest titer on 3 dpi (4.28 ± 0.14 EID50/mL), which decreased slightly until 7 dpi (3.86 ± 0.80 EID50/mL) but maintained relatively high viral titer.
Histopathological analysis of the lung tissues showed different degrees of inflammation between the groups challenged with 20X-20 and 34X-2 viruses. The inflammatory cells infiltrated into the peribronchiolar and perivascular regions of the lungs of mice infected with 20X-20 (5 and 7 dpi) and 34X-2 (3 and 7 dpi) are shown in Figure 4. Severe interstitial pneumonia and degeneration of the bronchial epithelium were observed in the lungs of mice infected with 34X-2 on 7 dpi.
Seroconversion and antigenic analyses were performed using mouse serum on 14 dpi. Seroconversion was observed in all mice and two of three mice challenged with 20X-20 and 34X-2, respectively (Table 4). The antigenic result of the cross-HI test revealed no difference in antigenic patterns between the two H7 Korean viruses. The 20X-20 virus showed identical HI titers against its homologous antiserum (8, 32, 32) and 34X-2 antiserum (8, 32, 32). The 34X-2 virus also showed identical HI titers against its homologous antiserum (8, 32) and 20X-2 antiserum (8, 32).

4. Discussion

LPAIVs cause mild respiratory syndromes and decreased feed intake and egg production in poultry, leading to significant economic losses [73,74]. Although H7 HPAIVs are yet to cause an outbreak in South Korea, H7 LPAIVs have caused an outbreak, resulting in economic losses. In 2010 and 2011, H7N2, H7N6, H7N7, and H7N8 were isolated from domestic ducks and over 197,000 poultry were culled to eradicate the disease [17]. Moreover, H7 LPAIVs can infect humans and cause pathogenicity [75]. Therefore, in this study, we analyzed the genetic characterization and pathogenicity of two LPAIV viruses, H7N7 and H7N9, isolated from wild birds in South Korea in 2021.
We isolated H7N9 (20X-20) and H7N7 (34X-2) AIVs through active influenza surveillance from October 2020 to March 2021. Both viruses have ELPKGR/GLF at the HA cleavage site, indicating they are LPAIVs, and the same cleavage site has been reported in human influenza virus [33]. In particular, the receptor binding sites involved in interspecies transmission are located in the HA gene [76]. There were no amino acid substitutions of G186V, E190D, Q/G225D, Q226L, S227N, and G228S in the HA gene with known effect on receptor binding to mammalian tissue, demonstrating that the two Korean H7 AIVs would likely have an avian receptor binding specificity. However, since the characterization was analyzed only at the sequence level, an actual receptor binding experiment would need to be performed to demonstrate avian receptor binding specificity.
The amino acid substitutions of E627K and D701N in PB2, which is a well-known mammalian virulence factor, were not observed in the two H7 Korean AIVs. However, several amino acid substitutions related to enhanced virulence in mice were identified in the two H7 Korean AIVs: L89V, G309D, and T339R in the PB2 gene; D622G in the PB1 gene; N66S in the PB1-F2 gene; V286A and M437T in the NP gene; N30D, I43M, and T215A in the M1 gene; and P42S, L103F, and C138F in the NS1 gene. Although the two H7 Korean AIVs have avian receptor binding preference in the HA gene, they showed pathogenicity in the murine model. This result suggests that the presence and combination of other amino acid mutations may affect virulence in mice. Therefore, continuous studies are needed on how amino acid mutations play an important role in LPAIVs replication in the mammalian cell. The two H7 Korean AIVs carried I117T amino acid mutations in the NA gene, which reduced the antiviral drug susceptibility. A/Jiangsu/1/2018 (H7N4; Jiangsu/1), which is closely related to chicken/1, was susceptible to antivirals (oseltamivir) and they had no I117T mutations in the NA gene [66]. This result suggests that the two H7 Korean AIVs may be resistant to antiviral drugs, but further studies are required to confirm this finding.
The eight segments of the two H7 Korean AIVs belonged to the Eurasian lineage. The MCMC tree of the HA gene revealed that the HA gene of the two H7 Korean AIVs clustered in the G1 group, suggesting that the HA gene of these two strains likely originated from the AIV circulating in Korea in 2017. Moreover, the PB2, PB1, PA, HA, and NA genes of 20X-20 and the PB2, PB1, HA, and NA genes of 34X-2 showed high nucleotide similarity to AIVs of Korean origin. The NP, M, and NS genes of 20X-20 and the PA, NP, and NS genes of 34X-2 showed high nucleotide identity with Asian AIV strains, suggesting that the segments were introduced to Korea by intracontinental reassortment. However, the M gene of the 34X-2 virus belonged to the American lineage, and this result demonstrates intercontinental reassortment of AIVs.
In the ML tree analysis of the HA gene, two H7 Korean AIVs showed low nucleotide similarity to the HA gene of A/Italy/3/2013(H7N7) and A/Anhui/1/2013(H7N9), which originate from AIVs and cause human infections. The ML tree analysis of the NA gene of the two H7 Korean AIVs showed the same results as the HA gene. However, the HA gene of the two H7 Korean AIVs showed a high similarity to the A/Jiangsu/1/2018(H7N4) virus, which recently caused a severe human infection in China [33]. China is adjacent to South Korea, and there is a possibility of novel viruses being introduced into South Korea. Therefore, continuous surveillance of AIVs is required.
Previous studies have shown that H7 LPAIVs can replicate in the lungs of murine models [77,78,79]. The two H7 Korean AIVs were isolated from the lungs of infected mice without prior host adaptation. All infected mice showed a 0% mortality rate. However, mice infected with 34X-2 showed body weight loss and clinical symptoms of ruffled fur on 7 dpi. In the histopathology analysis, inflammation in the peribronchiolar and perivascular regions was observed in both lungs of mice infected with the two H7 Korean AIVs. Inflammation was most severe in 34X-2-infected mice on 7 dpi. Further, 20X-20 and 34X-2 showed the highest viral titer of 2.80 log10EID50/mL on 5 dpi and 4.28 log10EID50/mL on 3 dpi, respectively. However, the two H7 Korean AIVs were not detected in the lungs of mice on 14 dpi. Moreover, the two H7 Korean AIVs indicated a low seropositive rate on 14 dpi, and even one of three mice infected with 34X-2 viruses showed no seroconversion. This result suggests that the two H7 Korean AIVs can replicate in mice and that the viruses can recover naturally. However, the lungs of the mice showed severe histopathological lesions, suggesting that AIVs can also show pathogenicity in other mammals. The antigenic analysis using the cross-HI test demonstrated that there were no significant antigenic differences between 20X-20 and 34X-2. However, 20X-20 and 34X-2 showed differences in pathogenicity in the murine model. AIVs can increase virulence by reassortment with other AIVs and the introduction of another segment may facilitate transmission of AIVs [80,81]. This result suggests that internal genes of 34X-2, especially the M, NP, and NA genes, which clustered into distinct subgroups with 20X-20, may enhance the pathogenicity of 34X-2.
Our study provides information about genetic characterization and pathogenicity in murine models of H7H7 and H7H9 AIVs isolated from wild birds in South Korea in 2021. AIVs have been a concern because of the potential risks of human infection. Human infection with H7N9 AIVs has been reported to be more prevalent than other subtypes [29]. Although the H7N7 and H7N9 AIVs are LPAIVs, they exhibited pathogenicity in a murine model. Our findings suggest that the two H7 Korean AIVs have the potential to infect mammals, including humans. Moreover, H7 LPAIVs have the potential to mutate to HPAIV. Therefore, continuous surveillance of H7 AIVs is required.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13102057/s1, Figure S1: Phylogenetic tree of the NA9 (a) and NA7 (b) gene based on nucleotides. The tree was analyzed using the maximum-likelihood (ML) method with 1000 bootstrap replication and only bootstrap values more than 50% are shown. The two H7 Korean avian influenza viruses (AIVs) isolated in this study are shown in bold. The human influenza viruses are indicated by red circles. Table S1: The list of reference sequences used in Figure 1 and Figure 2.

Author Contributions

Conceptualization: J.-K.O. Methodology: Y.-S.K. and E.-J.N. Data curation: Y.-J.K. and J.-S.P. Writing and editing: E.-J.N. Supervision: J.-K.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Education (2019R1A6A1A03033084). This subject is supported by the National Institute of Wildlife Disease Control and Prevention as “Specialized Graduate School Support Project for Wildlife Disease Specialists”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

References

  1. Lamb, R.A.; Krug, R.M. Orthomyxoviridae: The viruses and their replication. Fields Virol. 2001, 1, 1487–1531. [Google Scholar]
  2. Alexander, D.J. An overview of the epidemiology of avian influenza. Vaccine 2007, 25, 5637–5644. [Google Scholar] [CrossRef] [PubMed]
  3. Tong, S.; Zhu, X.; Li, Y.; Shi, M.; Zhang, J.; Bourgeois, M.; Yang, H.; Chen, X.; Recuenco-Cabrera, S.; Gómez, J.; et al. New world bats harbor diverse influenza a viruses. PLoS Pathog. 2013, 9, e1003657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tong, S.; Li, Y.; Rivailler, P.; Conrardy, C.; Castillo, D.A.A.; Chen, L.-M.; Recuenco-Cabrera, S.; Ellison, J.; Davis, C.T.; York, I.; et al. A distinct lineage of influenza A virus from bats. Proc. Natl. Acad. Sci. USA 2012, 109, 4269–4274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Gonzales, J.; Elbers, A.; Bouma, A.; Koch, G.; de Wit, J.; Stegeman, J. Transmission characteristics of low pathogenic avian influenza virus of H7N7 and H5N7 subtypes in layer chickens. Vet. Microbiol. 2012, 155, 207–213. [Google Scholar] [CrossRef]
  6. Shi, J.; Deng, G.; Kong, H.; Gu, C.; Ma, S.; Yin, X.; Zeng, X.; Cui, P.; Chen, Y.; Yang, H.; et al. H7N9 virulent mutants detected in chickens in China pose an increased threat to humans. Cell Res. 2017, 27, 1409–1421. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Horimoto, T.; Rivera, E.; Pearson, J.; Senne, D.; Krauss, S.; Kawaoka, Y.; Webster, R. Origin and molecular changes associated with emergence of a highly pathogenic H5N2 influenza virus in Mexico. Virology 1995, 213, 223–230. [Google Scholar] [CrossRef] [Green Version]
  8. Banks, J.; Speidel, E.S.; Moore, E.; Plowright, L.; Piccirillo, A.; Capua, I.; Cordioli, P.; Fioretti, A.; Alexander, D.J. Changes in the haemagglutinin and the neuraminidase genes prior to the emergence of highly pathogenic H7N1 avian influenza viruses in Italy. Arch. Virol. 2001, 146, 963–973. [Google Scholar] [CrossRef]
  9. Berhane, Y.; Hisanaga, T.; Kehler, H.; Neufeld, J.; Manning, L.; Argue, C.; Handel, K.; Hooper-McGrevy, K.; Jonas, M.; Robinson, J.; et al. Highly pathogenic avian influenza virus A (H7N3) in Domestic Poultry, Saskatchewan, Canada, 2007. Emerg. Infect. Dis. 2009, 15, 1492–1495. [Google Scholar] [CrossRef] [PubMed]
  10. Suarez, D.L.; Senne, D.A.; Banks, J.; Brown, I.H.; Essen, S.C.; Lee, C.-W.; Manvell, R.J.; Mathieu-Benson, C.; Moreno, V.; Pedersen, J.C.; et al. Recombination resulting in virulence shift in Avian Influenza Outbreak, Chile. Emerg. Infect. Dis. 2004, 10, 693–699. [Google Scholar] [CrossRef]
  11. Beerens, N.; Heutink, R.; Harders, F.; Bossers, A.; Koch, G.; Peeters, B. Emergence and selection of a highly pathogenic avian influenza H7N3 virus. J. Virol. 2020, 94, e01818-19. [Google Scholar] [CrossRef] [PubMed]
  12. Lee, D.-H.; Torchetti, M.K.; Killian, M.L.; Berhane, Y.; Swayne, D.E. Highly pathogenic avian influenza A (H7N9) virus, Tennessee, USA, March 2017. Emerg. Infect. Dis. 2017, 23, 1860. [Google Scholar] [CrossRef] [PubMed]
  13. Killian, M.L.; Kim-Torchetti, M.; Hines, N.; Yingst, S.; DeLiberto, T.; Lee, D.-H. Outbreak of h7n8 low pathogenic avian influenza in commercial turkeys with spontaneous mutation to highly pathogenic avian influenza. Genome Announc. 2016, 4, e00457-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Youk, S.; Lee, D.-H.; Killian, M.L.; Pantin-Jackwood, M.J.; Swayne, D.E.; Torchetti, M.K. Highly pathogenic avian influenza A (H7N3) virus in Poultry, United States, 2020. Emerg. Infect. Dis. 2020, 26, 2966. [Google Scholar] [CrossRef]
  15. Byrne, A.; Reid, S.; Seekings, A.; Núñez, A.; Prieto, A.O.; Ridout, S.; Warren, C.; Puranik, A.; Ceeraz, V.; Essen, S.; et al. H7N7 avian influenza virus mutation from low to high pathogenicity on a layer chicken farm in the UK. Viruses 2021, 13, 259. [Google Scholar] [CrossRef]
  16. Naguib, M.M.; Verhagen, J.H.; Mostafa, A.; Wille, M.; Li, R.; Graaf, A.; Järhult, J.D.; Ellström, P.; Zohari, S.; Lundkvist, Å.; et al. Global patterns of avian influenza A (H7): Virus evolution and zoonotic threats. FEMS Microbiol. Rev. 2019, 43, 608–621. [Google Scholar] [CrossRef] [Green Version]
  17. Abdelwhab, E.M.; Veits, J.; Mettenleiter, T.C. Prevalence and control of H7 avian influenza viruses in birds and humans. Epidemiol. Infect. 2014, 142, 896–920. [Google Scholar] [CrossRef]
  18. Kim, M.-C.; Jeong, O.-M.; Kang, H.-M.; Paek, M.-R.; Kwon, J.-S.; Song, C.-S.; Kwon, Y.-K.; Lee, J.-G.; Kwon, J.-H.; Lee, Y.-J. Pathogenicity and transmission studies of H7N7 avian influenza virus isolated from feces of magpie origin in chickens and magpie. Vet. Microbiol. 2010, 141, 268–274. [Google Scholar] [CrossRef]
  19. Kim, H.-R.; Park, C.-K.; Lee, Y.-J.; Oem, J.-K.; Kang, H.-M.; Choi, J.-G.; Lee, O.-S.; Bae, Y.-C. Low pathogenic H7 subtype avian influenza viruses isolated from domestic ducks in South Korea and the close association with isolates of wild birds. J. Gen. Virol. 2012, 93, 1278–1287. [Google Scholar] [CrossRef]
  20. Lee, E.-K.; Lee, Y.-N.; Kye, S.-J.; Lewis, N.S.; Brown, I.H.; Sagong, M.; Heo, G.-B.; Kang, Y.-M.; Cho, H.-K.; Kang, H.-M.; et al. Characterization of a novel reassortant H5N6 highly pathogenic avian influenza virus clade 2.3.4.4 in Korea, 2017. Emerg. Microbes Infect. 2018, 7, 103. [Google Scholar] [CrossRef] [Green Version]
  21. Herfst, S.; Imai, M.; Kawaoka, Y.; Fouchier, R.A.M. Avian influenza virus transmission to mammals. In Influenza Pathogenesis and Control; Current Topics in Microbiology and Immunology 385; Springer: Berlin, Germany, 2014; Volume I, pp. 137–155. [Google Scholar] [CrossRef]
  22. Reperant, L.; Rimmelzwaan, G.; Kuiken, T. Avian influenza viruses in mammals. Rev. Sci. Tech. L’oie 2009, 28, 137–159. [Google Scholar] [CrossRef]
  23. Kurtz, J.; Manvell, R.J.; Banks, J. Avian influenza virus isolated from a woman with conjunctivitis. Lancet 1996, 348, 901–902. [Google Scholar] [CrossRef]
  24. Fouchier, R.A.M.; Schneeberger, P.M.; Rozendaal, F.W.; Broekman, J.M.; Kemink, S.A.G.; Munster, V.; Kuiken, T.; Rimmelzwaan, G.F.; Schutten, M.; van Doornum, G.J.J.; et al. Avian influenza A virus (H7N7) associated with human conjunctivitis and a fatal case of acute respiratory distress syndrome. Proc. Natl. Acad. Sci. USA 2004, 101, 1356–1361. [Google Scholar] [CrossRef] [Green Version]
  25. Puzelli, S.; Rossini, G.; Facchini, M.; Vaccari, G.; Di Trani, L.; Di Martino, A.; Gaibani, P.; Vocale, C.; Cattoli, G.; Bennett, M.; et al. Human infection with highly pathogenic A (H7N7) avian influenza virus, Italy, 2013. Emerg. Infect. Dis. 2014, 20, 1741–1745. [Google Scholar] [CrossRef] [Green Version]
  26. Gao, R.; Cao, B.; Hu, Y.; Feng, Z.; Wang, D.; Hu, W.; Chen, J.; Jie, Z.; Qiu, H.; Xu, K.; et al. Human infection with a novel avian-origin influenza A (H7N9) virus. N. Engl. J. Med. 2013, 368, 1888–1897. [Google Scholar] [CrossRef] [Green Version]
  27. Wu, X.; Xiao, L.; Li, L. Research progress on human infection with avian influenza H7N9. Front. Med. 2020, 14, 8–20. [Google Scholar] [CrossRef] [Green Version]
  28. Shi, J.; Deng, G.; Ma, S.; Zeng, X.; Yin, X.; Li, M.; Zhang, B.; Cui, P.; Chen, Y.; Yang, H.; et al. Rapid evolution of H7N9 highly pathogenic viruses that emerged in China in 2017. Cell Host Microbe 2018, 24, 558–568.e7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Philippon, D.A.M.; Wu, P.; Cowling, B.J.; Lau, E.H.Y. Avian influenza human infections at the human-animal interface. J. Infect. Dis. 2020, 222, 528–537. [Google Scholar] [CrossRef] [PubMed]
  30. Yin, X.; Deng, G.; Zeng, X.; Cui, P.; Hou, Y.; Liu, Y.; Fang, J.; Pan, S.; Wang, D.; Chen, X.; et al. Genetic and biological properties of H7N9 avian influenza viruses detected after application of the H7N9 poultry vaccine in China. PLoS Pathog. 2021, 17, e1009561. [Google Scholar] [CrossRef] [PubMed]
  31. Zeng, X.; Tian, G.; Shi, J.; Deng, G.; Li, C.; Chen, H. Vaccination of poultry successfully eliminated human infection with H7N9 virus in China. Sci. China Life Sci. 2018, 61, 1465–1473. [Google Scholar] [CrossRef]
  32. Belser, J.A.; Bridges, C.B.; Katz, J.M.; Tumpey, T.M. Past, present, and possible future human infection with influenza virus a subtype H7. Emerg. Infect. Dis. 2009, 15, 859–865. [Google Scholar] [CrossRef]
  33. Huo, X.; Cui, L.-B.; Chen, C.; Wang, D.; Qi, X.; Zhou, M.-H.; Guo, X.; Wang, F.; Liu, W.J.; Kong, W.; et al. Severe human infection with a novel avian-origin influenza A(H7N4) virus. Sci. Bull. 2018, 63, 1043–1050. [Google Scholar] [CrossRef] [Green Version]
  34. Stear, M.J. OIE manual of diagnostic tests and vaccines for terrestrial animals (mammals, birds and bees) 5th Edn. Volumes 1 & 2. World organization for animal health 2004. ISBN 92 9044 622 6. €140. Parasitology 2005, 130, 727. [Google Scholar] [CrossRef]
  35. Hoffmann, B.; Hoffmann, D.; Henritzi, D.; Beer, M.; Harder, T.C. Riems influenza a typing array (RITA): An RT-qPCR-based low density array for subtyping avian and mammalian influenza a viruses. Sci. Rep. 2016, 6, 27211. [Google Scholar] [CrossRef] [Green Version]
  36. Lee, D.-H.; Lee, H.-J.; Kang, H.-M.; Jeong, O.-M.; Kim, M.-C.; Kwon, J.-S.; Kwon, J.-H.; Kim, C.-B.; Lee, J.-B.; Park, S.-Y.; et al. DNA barcoding techniques for avian influenza virus surveillance in migratory bird habitats. J. Wildl. Dis. 2010, 46, 649–654. [Google Scholar] [CrossRef] [PubMed]
  37. Hoffmann, E.; Stech, J.; Guan, Y.; Webster, R.G.; Perez, D. Universal primer set for the full-length amplification of all influenza A viruses. Arch. Virol. 2001, 146, 2275–2289. [Google Scholar] [CrossRef] [PubMed]
  38. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  39. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Guindon, S.; Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 2003, 52, 696–704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Drummond, A.J.; Rambaut, A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 2007, 7, 214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Reed, L.; Muench, H. A simple method of estimating fifty percent Endpoints12. Am. J. Epidemiol. 1938, 27, 493–497. [Google Scholar] [CrossRef]
  43. Long, J.S.; Mistry, B.; Haslam, S.M.; Barclay, W.S. Host and viral determinants of influenza a virus species specificity. Nat. Rev. Microbiol. 2019, 17, 67–81. [Google Scholar] [CrossRef] [PubMed]
  44. Chen, L.M.; Blixt, O.; Stevens, J.; Lipatov, A.S.; Davis, C.T.; Collins, B.E.; Cox, N.J.; Paulson, J.C.; Donis, R.O. In vitro evolution of H5N1 avian influenza virus toward human-type receptor specificity. Virology 2012, 422, 105–113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Shi, Y.; Zhang, W.; Wang, F.; Qi, J.; Wu, Y.; Song, H.; Gao, F.; Bi, Y.; Zhang, Y.; Fan, Z.; et al. Structures and receptor binding of hemagglutinins from human-infecting H7N9 influenza viruses. Science 2013, 342, 243–247. [Google Scholar] [CrossRef] [PubMed]
  46. 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]
  47. Jonges, M.; Welkers, M.; Jeeninga, R.; Meijer, A.; Schneeberger, P.; Fouchier, R.; De Jong, M.; Koopmans, M. Emergence of the virulence-associated PB2 E627K substitution in a fatal human case of highly pathogenic avian influenza virus A(H7N7) infection as determined by illumina ultra-deep sequencing. J. Virol. 2013, 88, 1694–1702. [Google Scholar] [CrossRef] [Green Version]
  48. Conenello, G.M.; Zamarin, D.; Perrone, L.A.; Tumpey, T.; Palese, P. A single mutation in the PB1-F2 of H5N1 (HK/97) and 1918 influenza a viruses contributes to increased virulence. PLoS Pathog. 2007, 3, e141. [Google Scholar] [CrossRef]
  49. Senne, D.A.; Panigrahy, B.; Kawaoka, Y.; Pearson, J.E.; Süss, J.; Lipkind, M.; Kida, H.; Webster, R.G. Survey of the hemagglutinin (HA) cleavage site sequence of H5 and H7 avian influenza viruses: Amino acid sequence at the HA cleavage site as a marker of pathogenicity potential. Avian Dis. 1996, 40, 425–437. [Google Scholar] [CrossRef]
  50. Fan, S.; Hatta, M.; Kim, J.H.; Halfmann, P.; Imai, M.; Macken, C.A.; Le, M.Q.; Nguyen, T.; Neumann, G.; Kawaoka, Y. Novel residues in avian influenza virus PB2 protein affect virulence in mammalian hosts. Nat. Commun. 2014, 5, 5021. [Google Scholar] [CrossRef] [Green Version]
  51. Kong, H.; Ma, S.; Wang, J.; Gu, C.; Wang, Z.; Shi, J.; Deng, G.; Guan, Y.; Chen, H. Identification of key amino acids in the PB2 and M1 proteins of H7N9 influenza virus that affect its transmission in guinea pigs. J. Virol. 2019, 94, e01180-19. [Google Scholar] [CrossRef] [Green Version]
  52. Hu, M.; Yuan, S.; Zhang, K.; Singh, K.; Ma, Q.; Zhou, J.; Chu, H.; Zheng, B.-J. PB2 substitutions V598T/I increase the virulence of H7N9 influenza A virus in mammals. Virology 2017, 501, 92–101. [Google Scholar] [CrossRef] [PubMed]
  53. Subbarao, E.K.; London, W.; Murphy, B.R. A single amino acid in the PB2 gene of influenza A virus is a determinant of host range. J. Virol. 1993, 67, 1761–1764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Gao, Y.; Zhang, Y.; Shinya, K.; Deng, G.; Jiang, Y.; Li, Z.; Guan, Y.; Tian, G.; Li, Y.; Shi, J.; et al. Identification of amino acids in HA and PB2 critical for the transmission of H5N1 avian influenza viruses in a mammalian host. PLoS Pathog. 2009, 5, e1000709. [Google Scholar] [CrossRef] [Green Version]
  55. Liang, L.; Jiang, L.; Li, J.; Zhao, Q.; Wang, J.; He, X.; Huang, S.; Wang, Q.; Zhao, Y.; Wang, G.; et al. Low polymerase activity attributed to pa drives the acquisition of the PB2 E627K mutation of H7N9 avian influenza virus in mammals. Mbio 2019, 10, e01162-19. [Google Scholar] [CrossRef] [Green Version]
  56. Suzuki, Y.; Uchida, Y.; Tanikawa, T.; Maeda, N.; Takemae, N.; Saito, T. Amino acid substitutions in PB1 of avian influenza viruses influence pathogenicity and transmissibility in chickens. J. Virol. 2014, 88, 11130–11139. [Google Scholar] [CrossRef] [Green Version]
  57. Feng, X.; Wang, Z.; Shi, J.; Deng, G.; Kong, H.; Tao, S.; Li, C.; Liu, L.; Guan, Y.; Chen, H. Glycine at position 622 in PB1 contributes to the virulence of H5N1 avian influenza virus in mice. J. Virol. 2016, 90, 1872–1879. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. 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] [Green Version]
  59. Song, J.; Feng, H.; Xu, J.; Zhao, D.; Shi, J.; Li, Y.; Deng, G.; Jiang, Y.; Li, X.; Zhu, P.; et al. The PA protein directly contributes to the virulence of H5N1 avian influenza viruses in domestic ducks. J. Virol. 2011, 85, 2180–2188. [Google Scholar] [CrossRef] [Green Version]
  60. 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] [Green Version]
  61. Tada, T.; Suzuki, K.; Sakurai, Y.; Kubo, M.; Okada, H.; Itoh, T.; Tsukamoto, K. NP body domain and PB2 contribute to increased virulence of H5N1 highly pathogenic avian influenza viruses in chickens. J. Virol. 2010, 85, 1834–1846. [Google Scholar] [CrossRef] [Green Version]
  62. Wasilenko, J.L.; Sarmento, L.; Pantin-Jackwood, M.J. A single substitution in amino acid 184 of the NP protein alters the replication and pathogenicity of H5N1 avian influenza viruses in chickens. Arch. Virol. 2009, 154, 969–979. [Google Scholar] [CrossRef]
  63. 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. 2016, 19, e12692. [Google Scholar] [CrossRef] [Green Version]
  64. Ma, S.; Zhang, B.; Shi, J.; Yin, X.; Wang, G.; Cui, P.; Liu, L.; Deng, G.; Jiang, Y.; Li, C.; et al. Amino acid mutations A286V and T437M in the nucleoprotein attenuate H7N9 viruses in mice. J. Virol. 2020, 94, e01530-19. [Google Scholar] [CrossRef] [Green Version]
  65. Bi, Y.; Xiao, H.; Chen, Q.; Wu, Y.; Fu, L.; Quan, C.; Wong, G.; Liu, J.; Haywood, J.; Liu, Y.; et al. Changes in the length of the neuraminidase stalk region impact H7N9 virulence in mice. J. Virol. 2016, 90, 2142–2149. [Google Scholar] [CrossRef] [Green Version]
  66. Kode, S.S.; Pawar, S.D.; Tare, D.S.; Keng, S.S.; Hurt, A.C.; Mullick, J. A novel I117T substitution in neuraminidase of highly pathogenic avian influenza H5N1 virus conferring reduced susceptibility to oseltamivir and zanamivir. Vet. Microbiol. 2019, 235, 21–24. [Google Scholar] [CrossRef] [PubMed]
  67. Nao, N.; Kajihara, M.; Manzoor, R.; Maruyama, J.; Yoshida, R.; Muramatsu, M.; Miyamoto, H.; Igarashi, M.; Eguchi, N.; Sato, M.; et al. A single amino acid in the M1 protein responsible for the different pathogenic potentials of H5N1 highly pathogenic avian influenza virus strains. PLoS ONE 2015, 10, e0137989. [Google Scholar] [CrossRef] [Green Version]
  68. Cheung, C.; Rayner, J.M.; Smith, G.J.; Wang, P.; Naipospos, T.S.P.; Zhang, J.; Yuen, K.-Y.; Webster, R.G.; Peiris, J.S.M.; Guan, Y.; et al. Distribution of amantadine-resistant H5N1 avian influenza variants in Asia. J. Infect. Dis. 2006, 193, 1626–1629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. 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] [PubMed] [Green Version]
  70. Li, W.; Wang, G.; Zhang, H.; Xin, G.; Zhang, D.; Zeng, J.; Chen, X.; Xu, Y.; Cui, Y.; Li, K. Effects of NS1 variants of H5N1 influenza virus on interferon induction, TNFα response and p53 activity. Cell. Mol. Immunol. 2010, 7, 235–242. [Google Scholar] [CrossRef] [Green Version]
  71. Dankar, S.K.; Wang, S.; Ping, J.; Forbes, N.E.; Keleta, L.; Li, Y.; Brown, E.G. Influenza a virus NS1 gene mutations F103L and M106I increase replication and virulence. Virol. J. 2011, 8, 13. [Google Scholar] [CrossRef] [Green Version]
  72. Li, J.; Zhang, K.; Chen, Q.; Zhang, X.; Sun, Y.; Bi, Y.; Zhang, S.; Gu, J.; Li, J.; Liu, D.; et al. Three amino acid substitutions in the NS1 protein change the virus replication of H5N1 influenza virus in human cells. Virology 2018, 519, 64–73. [Google Scholar] [CrossRef]
  73. Rutten, N.; Gonzales, J.L.; Elbers, A.R.W.; Velthuis, A. Cost analysis of various low pathogenic avian influenza surveillance systems in the Dutch egg layer sector. PLoS ONE 2012, 7, e33930. [Google Scholar] [CrossRef]
  74. Thompson, J.; Seitzinger, A.H. Economic evaluation of low pathogenic avian influenza in northeastern US live bird markets. J. Appl. Poult. Res. 2019, 28, 78–84. [Google Scholar] [CrossRef]
  75. Liu, J.; Xiao, H.; Wu, Y.; Liu, D.; Qi, X.; Shi, Y.; Gao, G.F. H7N9: A low pathogenic avian influenza A virus infecting humans. Curr. Opin. Virol. 2014, 5, 91–97. [Google Scholar] [CrossRef]
  76. Shi, Y.; Wu, Y.; Zhang, W.; Qi, J.; Gao, G.F. Enabling the ‘host jump’: Structural determinants of receptor-binding specificity in influenza A viruses. Nat. Rev. Genet. 2014, 12, 822–831. [Google Scholar] [CrossRef] [PubMed]
  77. Zanin, M.; Koçer, Z.A.; Poulson, R.L.; Gabbard, J.D.; Howerth, E.W.; Jones, C.A.; Friedman, K.; Seiler, J.; Danner, A.; Kercher, L.; et al. Potential for low-pathogenic avian H7 influenza a viruses to replicate and cause disease in a mammalian model. J. Virol. 2017, 91, e01934-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Tang, W.; Li, X.; Tang, L.; Wang, T.; He, G. Characterization of the low-pathogenic H7N7 avian influenza virus in Shanghai, China. Poult. Sci. 2020, 100, 565–574. [Google Scholar] [CrossRef] [PubMed]
  79. Belser, J.A.; Sun, X.; Brock, N.; Pulit-Penaloza, J.A.; Jones, J.; Zanders, N.; Davis, C.T.; Tumpey, T.M.; Maines, T.R. Mammalian pathogenicity and transmissibility of low pathogenic avian influenza H7N1 and H7N3 viruses isolated from North America in 2018. Emerg. Microbes Infect. 2020, 9, 1037–1045. [Google Scholar] [CrossRef]
  80. Zhu, W.; Dong, J.; Zhang, Y.; Yang, L.; Li, X.; Chen, T.; Zhao, X.; Wei, H.; Bo, H.; Zeng, X.; et al. A gene constellation in avian influenza a (H7N9) viruses may have facilitated the fifth wave outbreak in China. Cell Rep. 2018, 23, 909–917. [Google Scholar] [CrossRef] [Green Version]
  81. Tscherne, D.M.; García-Sastre, A. Virulence determinants of pandemic influenza viruses. J. Clin. Investig. 2011, 121, 6–13. [Google Scholar] [CrossRef]
Viruses 13 02057 g001 550
Figure 1. Phylogenetic tree of the HA gene of the H7 avian influenza virus (AIV) nucleotides. The tree was analyzed using the maximum-likelihood (ML) method with 1000 bootstrap replication and only bootstrap values more than 50% are shown. The two H7 Korea AIVs isolated in this study are shown in bold. The human influenza viruses are indicated by red circles.
Figure 1. Phylogenetic tree of the HA gene of the H7 avian influenza virus (AIV) nucleotides. The tree was analyzed using the maximum-likelihood (ML) method with 1000 bootstrap replication and only bootstrap values more than 50% are shown. The two H7 Korea AIVs isolated in this study are shown in bold. The human influenza viruses are indicated by red circles.
Viruses 13 02057 g001
Viruses 13 02057 g002 550
Figure 2. Maximum clade credibility (MCC) trees of the HA gene of H7 LPAIVs isolated between 2016 and 2021 in Asia. The MCC tree were constructed using the uncorrelated lognormal distribution relaxed clock method in BEAST v1.10.4. The ESS values were greater than 200. Posterior probabilities >0.8 are provided in the tree. The horizontal axis indicates the time scale, and the unit is 10 years. The two H7 Korean AIVs are shown in red. The H7 LPAIVs circulating in South Korea in 2017 are shown in orange.
Figure 2. Maximum clade credibility (MCC) trees of the HA gene of H7 LPAIVs isolated between 2016 and 2021 in Asia. The MCC tree were constructed using the uncorrelated lognormal distribution relaxed clock method in BEAST v1.10.4. The ESS values were greater than 200. Posterior probabilities >0.8 are provided in the tree. The horizontal axis indicates the time scale, and the unit is 10 years. The two H7 Korean AIVs are shown in red. The H7 LPAIVs circulating in South Korea in 2017 are shown in orange.
Viruses 13 02057 g002
Viruses 13 02057 g003 550
Figure 3. Body weight changes in mice infected with 20X-20 and 34X-2 viruses. BALB/c were intranasally inoculated with 107EID50 of each virus. The mice were observed for 14 days. The data indicate the mean and standard error of the mean (SEM). The data were analyzed using two-way ANOVA with Tukey’s multiple comparisons test. The asterisks show that the body weight of 34X-2 is significantly different from that of 20X-20 and the PBS group (** p < 0.01).
Figure 3. Body weight changes in mice infected with 20X-20 and 34X-2 viruses. BALB/c were intranasally inoculated with 107EID50 of each virus. The mice were observed for 14 days. The data indicate the mean and standard error of the mean (SEM). The data were analyzed using two-way ANOVA with Tukey’s multiple comparisons test. The asterisks show that the body weight of 34X-2 is significantly different from that of 20X-20 and the PBS group (** p < 0.01).
Viruses 13 02057 g003
Viruses 13 02057 g004 550
Figure 4. Histopathology analysis of the lungs of mice inoculated with 20X-20 and 34X-2. Lung tissue sections from mice were stained with hematoxylin and eosin. Inflammatory cells infiltrated into the peribronchiolar (black arrow) and perivascular regions (red arrow) of the lungs of mice infected with 20X-20 (H7N9) and 34X-2 (H7N7) viruses can be visualized. Interstitial pneumonia and degeneration of bronchial epithelium (asterisk) were observed in 34X-2 (H7N7)-infected mice on 7 dpi. The original magnification is ×100 and the scale bar indicates 200 µm. Dpi, days post-infection.
Figure 4. Histopathology analysis of the lungs of mice inoculated with 20X-20 and 34X-2. Lung tissue sections from mice were stained with hematoxylin and eosin. Inflammatory cells infiltrated into the peribronchiolar (black arrow) and perivascular regions (red arrow) of the lungs of mice infected with 20X-20 (H7N9) and 34X-2 (H7N7) viruses can be visualized. Interstitial pneumonia and degeneration of bronchial epithelium (asterisk) were observed in 34X-2 (H7N7)-infected mice on 7 dpi. The original magnification is ×100 and the scale bar indicates 200 µm. Dpi, days post-infection.
Viruses 13 02057 g004
Table
Table 1. Molecular characterization of H7 avian influenza virus (AIV) isolates.
Table 1. Molecular characterization of H7 avian influenza virus (AIV) isolates.
Viral ProteinAmino Acid ResidueVirus StrainsCommentsReference
20X-20 (2021)34X-2
(2021)		Chicken/1
(2018)		Jiangsu/1
(2018)		Anhui/1
(2013)		Italy/3
(2013)
HACleavage siteELPKGR/GLFELPKGR/GLFELPKGR/GLFELPKGR/GLFELPKGR/GLFETPKRRERR/GLFLPAIV
-monobasic		[49]
G186VGGGGVGIncreased α2-6 binding[45]
E190DEEEEEE[43]
Q/G225DGGGGGG[43]
Q226LQQQQLQ[43]
S227NSSSSSS[44]
G228SGGGGGG[43]
PB2L89VVVVVVVIncreased polymerase activity in mammalian cell lines and mice[46]
I147TIIIIIT[50]
I292VIIIIVI[51]
G309DDDDDDD[46]
T339KKKKKKK[46]
K389RRRRRKR[52]
E627KEEEKKEIncreased virulence in mice[53]
D701NDDDDDD[54,55]
V598T/ITTTTVTIncreased polymerase activity in mammalian cell line[52]
PB1C38YYYYYYYIncreased polymerase activity in mammalian cell line and pathogenicity in chicken[56]
D622GGGGGGGIncreased polymerase activity and virulence in mice[57]
PB1-F2N66SSSNNNNIncreased virulence in mice[48]
PAS37AAAAASAIncreased polymerase activity in mammalian cell line[58]
K142RRKKKKK[55]
N383DDDDDDDIncreased pathogenicity in ducks[59]
N409SSSSSNSIncreased polymerase activity in mammalian cell line[58]
NPI41VIIIIIIIncreased polymerase activity in mammalian cell line[60]
M105VMVMMVVIncreased virulence in chicken[61]
A184KKKKKKK[62]
F253IIIIIIIIncreased virulence in mice[63]
V286AAAAAAA[64]
M437TTTTTTT
NA69–73 deletion
(QISNT)		NoNoNoNoYesNoIncreased virulence in mice[65]
I117TTTIITTIncreased resistance to antiviral drugs (oseltamivir and zanamivir)[66]
M1N30DDDDDDDIncreased virulence in mice[50]
I43MMMMMMM[67]
T215AAAAAAA[50]
M2L26FLLLLLLIncreased resistance to antiviral drugs (amantadine and rimantadine)[68]
S31NSSSSNS[68]
NS1P42SSSSSSSIncreased virulence in mice[69]
D92EDDDDDDIncreased virulence in mice[70]
L103FFFFFLFIncreased replication and virulence in mice[71]
C138FFFFFFFIncreased replication in mammalian cell line[72]
Table
Table 2. Sequence identities of the A/mallard/South Korea/JB20X-20/2021 (20X-20) and A/mallard/South Korea/JB34X-2/2021 (34X-2) genome.
Table 2. Sequence identities of the A/mallard/South Korea/JB20X-20/2021 (20X-20) and A/mallard/South Korea/JB34X-2/2021 (34X-2) genome.
Gene20X-20Genetic Identity34X-2Genetic Identity
PB2A/wild_duck/South_Korea/KNU18-106/2018(H7N7)98.3%A/red-crowned crane/South Korea/H1026/2017(H7N7)96.6%
PB1A/red-crowned crane/South Korea/H1026/2017(H7N7)95%A/wild_duck/South_Korea/KNU18-104/2018(H7N7)93.6%
PAA/wild_duck/South_Korea/KNU18-106/2018(H7N7)99.1%A/wild_bird/Eastern_China/1758/2017(H5N3)98.6%
HAA/wild_duck/South_Korea/KNU18-106/2018(H7N7)97.8%A/wild_duck/South_Korea/KNU18-106/2018(H7N7)97.3%
NPA/common teal/Shanghai/NH110923/2019(H1N1)98.8%A/duck/Mongolia/926/2019(H5N3)99.3%
NAA/Anas platyrhynchos/South Korea/JB31-96/2019(H11N9)98.2%A/mallard/Korea/A15/2016(H7N7)97.4%
MPA/duck/Mongolia/916/2018(H3N8)99.7%A/northern pintail/Alaska/362/2013(H3N8)99.4%
NSA/duck/Bangladesh/37509/2019(H8N4)99.6%A/duck/Bangladesh/37509/2019(H8N4)99.6%
Table
Table 3. Replication of H7 isolates in the lungs of 6-week-old BALB/c mouse. Values shown are the number of infected mice per number of inoculated mice. Values in parentheses are viral titers (log10EID50/mL). The viral titers are expressed as the mean ± standard deviation (SD) of the samples. The viral titers from the lungs of mice were calculated using the Reed and Muench method.
Table 3. Replication of H7 isolates in the lungs of 6-week-old BALB/c mouse. Values shown are the number of infected mice per number of inoculated mice. Values in parentheses are viral titers (log10EID50/mL). The viral titers are expressed as the mean ± standard deviation (SD) of the samples. The viral titers from the lungs of mice were calculated using the Reed and Muench method.
IsolateDays Post Infection (Mean ± SD, log10EID50/mL)
135714
20X-20
(H7N9)		3/3
(2.36 ± 1.01)		3/3
(1.31 ± 2.28)		3/3
(2.80 ± 2.43)		3/3
(0.48 ± 0.83)		0/3
34X-2
(H7N7)		3/3
(2.71 ± 2.36)		3/3
(4.28 ± 0.14)		3/3
(4.14 ± 0.10)		3/3
(3.86 ± 0.80)		0/3
Dpi, days post-infection; EID50, 50% Egg-infective dose.
Table
Table 4. Seroconversion in 6-week-old BALB/c mouse infected with H7 isolates and antigenic analysis of H7 isolates. Serum was collected on 14 dpi and treated for 18 h with a receptor-destroying enzyme. Seroconversion was confirmed using the hemagglutination inhibition (HI) assay.
Table 4. Seroconversion in 6-week-old BALB/c mouse infected with H7 isolates and antigenic analysis of H7 isolates. Serum was collected on 14 dpi and treated for 18 h with a receptor-destroying enzyme. Seroconversion was confirmed using the hemagglutination inhibition (HI) assay.
Seroconversion: Positive/Total (HI Titers)
SerumVirus20X-2034X-2
Antigen
20X-203/3 (8, 32, 32)2/3 (8, 32)
34X-23/3 (8, 32, 32)2/3 (8, 32)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Na, E.-J.; Kim, Y.-S.; Kim, Y.-J.; Park, J.-S.; Oem, J.-K. Genetic Characterization and Pathogenicity of H7N7 and H7N9 Avian Influenza Viruses Isolated from South Korea. Viruses 2021, 13, 2057. https://doi.org/10.3390/v13102057

AMA Style

Na E-J, Kim Y-S, Kim Y-J, Park J-S, Oem J-K. Genetic Characterization and Pathogenicity of H7N7 and H7N9 Avian Influenza Viruses Isolated from South Korea. Viruses. 2021; 13(10):2057. https://doi.org/10.3390/v13102057

Chicago/Turabian Style

Na, Eun-Jee, Young-Sik Kim, Yoon-Ji Kim, Jun-Soo Park, and Jae-Ku Oem. 2021. "Genetic Characterization and Pathogenicity of H7N7 and H7N9 Avian Influenza Viruses Isolated from South Korea" Viruses 13, no. 10: 2057. https://doi.org/10.3390/v13102057

APA Style

Na, E. -J., Kim, Y. -S., Kim, Y. -J., Park, J. -S., & Oem, J. -K. (2021). Genetic Characterization and Pathogenicity of H7N7 and H7N9 Avian Influenza Viruses Isolated from South Korea. Viruses, 13(10), 2057. https://doi.org/10.3390/v13102057

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