DATA REPORT article

Front. Vet. Sci., 29 July 2025

Sec. Veterinary Epidemiology and Economics

Volume 12 - 2025 | https://doi.org/10.3389/fvets.2025.1638067

First detection of clade 2.3.4.4b H5N1 highly pathogenic avian influenza virus in a wild leopard cat (Prionailurus bengalensis) in South Korea

  • 1. Wildlife Disease Research Team, National Institute of Wildlife Disease Control and Prevention, Gwangju, Republic of Korea

  • 2. Wildlife Health Laboratory, College of Veterinary Medicine, Konkuk University, Seoul, Republic of Korea

  • 3. Avian Disease Laboratory, College of Veterinary Medicine, Konkuk University, Seoul, Republic of Korea

1 Introduction

The A/goose/Guangdong/1/1996 (Gs/GD) lineage of highly pathogenic avian influenza (HPAI) viruses initially identified in China in 1996 and have evolved over subsequent decades, characterized by accumulation of point mutations and multiple reassortment events with low pathogenic avian influenza (LPAI) viruses (, ). Particularly those of the H5N1 subtype within clade 2.3.4.4b have emerged as a significant threat to poultry, wild birds, and mammals worldwide (). Since their widespread dissemination in wild bird population in the early 2020s, these viruses have caused multiple sporadic infections in mammals, including small carnivores, marine mammals, cattle, and humans (). These cross-species transmissions, often linked to the consumption of infected birds or exposure to contaminated environments, have raised significant concerns about the zoonotic potential of clade 2.3.4.4b H5N1 and its capacity to evolve into a strain with pandemic potential.

Since its initial detection in 2014 with the H5N8 subtype in wild birds, multiple subtypes of clade 2.3.4.4b highly pathogenic avian influenza viruses (HPAIVs) have caused recurrent outbreaks in South Korea including H5N1, H5N6, and H5N8 (). Notable epidemics occurred during 2021–2022, with 44 cases primarily of H5N1 identified in wild birds (), and 2022–2023, with 174 wild bird cases, highlighting its persistent circulation wild waterfowl population (). These outbreaks underscore the role of wild birds in introduction and spread of the virus across South Korea, posing ongoing challenges to poultry industry and public health. In addition, although HPAIV infection in domestic cats were reported in 2016 and 2023 (), no instances of wild mammalian infection were documented in South Korea during these epidemics.

Here, we report the first documented case of H5N1 HPAI in a wild mammal in South Korea, identified in a leopard cat (Prionailurus bengalensis). On March 18, 2025, a wild leopard cat discovered moribund near a freshwater reservoir in Hwasun County of Jeollanam-do Province and submitted to the National Institute of Wildlife Disease Control and Prevention (NIWDC) of Korea. We isolated the H5N1 virus from this leopard cat, sequenced, and assessed its evolutionary history and molecular markers indicative of mammalian adaptation.

2 Materials and methods

2.1 Sample collection and virus isolation

On March 18, 2025, a wild leopard cat was found moribund near the reservoir in Hwasun County, Jeollanam-do Province, South Korea (GPS coordinates ≈ 35°03′N, 126°59′E). The wild leopard cat transported to the Jeollanamdo wild animal rescue center and died within a few hours. The carcass was submitted to the biosafety level 3 facility of NIWDC and organs including brain, trachea, and lung were collected. Samples were placed in phosphate-buffered saline with 400 mg/ml gentamicin, homogenized by vortexing, and filtered using a 0.45-μm Minisart Syringe Filter (Sartorius, Göttingen, Germany) following centrifugation at 3,000 rpm for 10 min. Filtered supernatant was inoculated into 10-day-old specific-pathogen-free embryonated chicken eggs and incubated at 37°C for 72 h. Allantoic fluids were harvested and tested for hemagglutination activity (HA) using 0.5% chicken red blood cells. RNA was extracted from tissue samples and HA-positive allantoic fluids using the Maxwell RSC simply RNA Tissue Kit (Promega, Madison, WI, USA) and screened for influenza A matrix and H5 genes via real-time reverse transcription-PCR (rRT-PCR), following established protocols ().

2.2 Whole genome sequencing and sequence analysis

Complementary DNA was synthesized using the SuperScript III First-Strand Synthesis System (Invitrogen, Carlsbad, CA, USA), and the eight gene segments were amplified with AccuPrime Pfx DNA Polymerase (Invitrogen, Carlsbad, CA, USA), as previously described (). DNA libraries were prepared using the Illumina DNA Prep Kit (Illumina, San Diego, CA, USA) and sequenced on the Illumina MiSeq platform (paired-end 150 bp). Raw reads were trimmed using BBDuk (v38.84) with a minimum quality threshold of 30 (), assembled de novo with SPAdes (v3.15.5). Trimmed reads were mapped to the top BLAST result from the GISAID EpiFlu database using Minimap2 (v2.24). Consensus sequences were generated using Geneious Prime software and deposited in GISAID Epiflu database (EPI_ISL_20051149).

2.3 Phylogenetic and mutation analysis

The top 250 hits for each segment query were retrieved and sequences with high identity (ranging from 99.5 to 99.9%, depending on the segment) were removed using CD-hit (). We also included genome sequences of four H5N1 HPAIVs [A/Wild_Duck/Korea/24WF364-8P/2024, A/Eurasian_wigeon/Korea/24WF382-7P/2024, A/Vulture/Korea/24WC103/2024, and A/Bean_goose/Korea/24WC196/2025] which were isolated from wild birds during the winter season of 2024-2025, all of which have been deposited in GISAID with their respective accession numbers (EPI_ISL_20051150, 19832581-19832583). Phylogenetic trees were constructed for each gene segment using RAxML v8.0 with the general time reversible model and 1,000 bootstrap replicates (). Interactive Tree of Life (iTOL) was used to visualize the tree of each gene (). A Bayesian relaxed-clock phylogeny of the hemagglutinin (HA) gene was reconstructed using BEAST version 1.10.4 (), employing the Hasegawa-Kishino-Yano substitution model with an uncorrelated log-normal distribution and a Gaussian Markov Random Field (GMRF) Bayesian skyride coalescent prior (). The Markov Chain Monte Carlo (MCMC) process was executed in parallel across three chains, each comprising 50 million iterations, with results combined after a 10% burn-in. All parameters achieved effective sample sizes (ESS) >200 and were assessed using TRACER v1.5 (http://tree.bio.ed.ac.uk/software/tracer/) (). A maximum clade credibility (MCC) tree was generated using TreeAnnotator and visualized with FigTree v1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/).

Molecular markers of mammalian adaptation, pathogenicity, and drug resistance were identified using the FluMut tool (). In the FluMut analysis, in addition to the sequence isolated from the leopard cat, seven of clade 2.3.4.4b HPAI H5N1 viruses were additionally analyzed including viruses reported in infected mammals in the United States (), viruses isolated from domestic cats in Korea in 2023 (), and viruses isolated from wild birds during the 2024–2025 winter season ().

3 Descriptive results

3.1 Isolation and genome sequencing of the virus

The brain and trachea from the submitted leopard cat tested positive for influenza A virus via chicken embryonated egg inoculation and rRT-PCR. The isolated virus, designated A/Leopard Cat/Korea/24WM130/2025(H5N1) (hereafter 24WM130) yielded 288,814 NGS reads, enabling assembly of complete coding genome sequences across all eight influenza virus segments.

3.2 Genome analysis

The virus was identified as HPAIV based on the presence of multiple basic amino acids at the HA proteolytic cleavage site (PLREKRRKR/G) (). All gene segments of the 24WM130 virus clustered with clade 2.3.4.4b H5N1 HPAIVs isolated from wild birds in South Korea during the 2024–2025 winter season, showing close genetic relatedness and a likely origin from infected wild birds through predation or scavenging (Supplementary Figure 1). We previously reported two genotypes of H5N1 clade 2.3.4.4b viruses in October 2024, the genotype 1 and 2, represented by A/Northern pintail/Korea/24WC025/2024 virus and A/Mandarin duck/Korea/24WS005-2/2024 virus, respectively (). Genotype 1 possessed a G2d-lineage HA gene and a genome constellation identical to strains circulating in Japan during 2023–2024. Genotype 2 carried a G2c-lineage HA gene, neuraminidase (NA) and M genes from H5Nx clade 2.3.4.4b viruses circulating in 2022–2024, and internal genes from LPAIVs in the East Asian–Australasian flyway. The 24WM130 virus had a G2d-lineage HA gene, while its remaining segments closely matched those of genotype 2, indicating a reassortant virus derived from early HPAI outbreaks in October 2024 () (Figure 1B). These reassortment events are likely driven by the high density and mobility of migratory birds, which promote co-infection and gene exchange. In the Bayesian phylogenetic analysis of the HA gene, the 24WM130 virus clustered with clade 2.3.4.4b H5N1 HPAIVs from wild birds in South Korea during December 2024–February 2025 and supported by a high posterior probability (0.99; Figure 1A). Their tMRCA was estimated to be August 2024 (95% BCI: April 24–November 20, 2024), suggesting that these viruses originated most likely during the late breeding season of waterfowl in Eurasia.

Figure 1

Leopard cat is an endangered, solitary, and opportunistic mesocarnivore species native to South Korea and other parts of Asia. Their solitary lifestyle, coupled with the absence of additional HPAI detections in other mammals, suggests that this outbreak was most likely a sporadic event. The leopard cat's ecological overlap with wild birds and free-ranging mammals raises concerns about potential disease spillovers and their role as an intermediate species in future pandemics ().

Among the 36 mammalian adaptation markers identified in 24WM130, three mutations (Table 1), I292V in PB2, D154N in HA, and D74N in NS, were unique compared to other H5N1 viruses identified in mammals. The I292V in PB2 and D74N in NS substitutions have been associated with increased polymerase activity in mammalian hosts and enhanced virulence in mice. The D154N in HA has been reported to increase binding affinity to α2,6-linked sialic acid receptors (, ). Notably, the 24WM130 virus did not possess the E627K or D701N in PB2 which were previously identified in clade 2.3.4.4b HPAI H5N1 viruses isolated from domestic cats in South Korea and M631L in PB2 which was a unique mutation found in dairy cows in the U.S. The HPAI H5N1 viruses from wild birds in Korea during the same wintering season also had 35 mammalian adaptation markers out of 36 detected in the 24WM130, highlighting the high potential for clade 2.3.4.4b H5N1 viruses in wild birds to spillover to mammalian hosts. The D74N substitution, detected exclusively in the 24WM130 virus and absent from closely related avian viruses in Korea during 2024–2025, suggests it may have emerged in the leopard cat.

Table 1

GeneMutationaA/Leopard Cat/ Korea/24WM 130/2025A/feline/South Korea/SNU-01/2023A/feline/ Korea/ M305-3/ 2023A/Mandarin duck/ Korea/24WS 005-2/2024A/Northern pintail/Korea/ 24WC25/2024A/Spot-billed duck/Korea/ 24WF364-8PA/dairy cow/USA/ 007549-005/2025A/Louisiana/ 12/2024Effecta
PB2L89VVVVVVVV• Increased polymerase activity in mammalian, virulence in mice
I292VIIVIVII
G309DDDDDDDD
T339KKKKKKKK
K389RRRRRRRR
R447GGGGGGGG
K482RKKKKRKKK
I495VVVVVVIV
V598TTTTTTTT
M631LMMMMMMLM
A676TTTTTTAT
E627KEEKEEEEE• Increased polymerase activity, replication in mammalian
• Increased virulence in mice, ferrets
• Contributes to airborne pathogenicity in ferrets
D701NDNDDDDDD• Increased polymerase activity, replication in mammalian
• Increased virulence in mice, guinea pigs
PB1D3VVVVVVVV• Increased polymerase activity in mammalian
F2:N66SSSSSSNS
R207KKKKKKKK
D622GGGGGGGG
PAS37AAAAAAAA• Increased polymerase activity in mammalian
N383DDDDDDDD
N409SSSSSSSS
K497RKKKKKRK
HAS133AAAAAAAA• Increase virus binding to α2–6
D154NDDDNNNN
K218QQQQQQQQ
S223RRRRRRRR
S107RRRRRRRR• Increase virulence in mice
T108IIIIIIII
T156AAAAAAAA• Increase virus binding to α2–6
• Increase transmission in guinea pigs
N189DNNNDNNNN• Increase virus binding to α2–6
• Transmissible among ferrets
NPN319KNNNNKNNN• Increased polymerase activity, replication in mammalian
MP
(M1)
N30DDDDDDDD• Increased virulence in mice
I43MMMMMMMM
T215AAAAAAAA
NS
(NS-1)
P3SSSSSSSS• Increase replication in mammalian
R41KKKKKKKK
K55EEEEEEEE
K66EEEEEEEE
C138FFFFFFFF
P42SSSSSSSS• Increase virulence in mice
D74NDDDDDDD• Increase virulence in mice
I106MMMMMMMM• Increase replication in mammalian
T48AAAAAAAA• Increase virulence in ferret

Comparison of mammalian adaptation markers among A/Leopard Cat/Korea/WM130/2025, other mammalian-derived strains, and clade 2.3.4.4b highly pathogenic H5N1 viruses isolated from Korean wild birds in the same season.

a

Mammalian adaptation mutations and their effects were identified using the FluMut tool.

b

Identified substitutions associated with mammalian adaptation are indicated in bold. Substitutions detected in A/Leopard Cat/Korea/WM130/2025 are highlighted in red.

4 Conclusion

The identification of clade 2.3.4.4b H5N1 in a wild leopard cat in South Korea highlights the evolving epidemiology of HPAI and the need for expanded surveillance in wild mammals. Although infections in domestic cats were reported in South Korea in 2016 and 2023 (), this case represents the first confirmed detection of HPAI H5N1 in a wild mammalian species in the country. Given the role of wild birds in the spread and maintenance of HPAIVs, the significant presence of mammalian adaptation markers in viruses circulating in wild bird populations is of great concern, particularly for potential widespread dissemination and interspecies transmissions. Genomic data from this study, shared via GISAID, will support ongoing efforts to track viral spread and assess zoonotic risks. Enhanced monitoring of predator-prey interactions and mammalian populations near avian habitats is recommended to mitigate future spillover events. Further in-depth genomic and phenotypic analyses are warranted to better understand the pathobiological characteristics and zoonotic potential of this virus.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Ethics statement

Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because it involved a naturally deceased wild animal (carcass), and no live animal experimentation was conducted.

Author contributions

Y-JS: Writing – review & editing, Resources, Formal analysis. S-HL: Investigation, Writing – review & editing, Writing – original draft, Formal analysis, Data curation, Methodology. D-JK: Resources, Formal analysis, Writing – review & editing. KL: Writing – review & editing, Resources. M-aL: Resources, Writing – review & editing. D-YL: Investigation, Writing – original draft, Formal analysis. Y-RS: Data curation, Formal analysis, Investigation, Writing – original draft. HJ: Funding acquisition, Resources, Writing – review & editing. SL: Funding acquisition, Writing – review & editing, Resources. D-HL: Conceptualization, Supervision, Writing – original draft, Investigation, Funding acquisition, Writing – review & editing, Data curation.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was financially supported by a grant from the National Institute of Wildlife Disease Control and Prevention (NIWDC), Ministry of Environment, Republic of Korea (grant number 2025-006).

Acknowledgments

We thank the GISAID Initiative for providing access to reference sequences. This paper was also supported by the KU Research Professor Program of Konkuk University.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2025.1638067/full#supplementary-material

Supplementary Figure 1

Maximum-likelihood tree constructed using the complete coding nucleotide sequences of (A) polymerase basic protein 2, (B) polymerase basic protein 1, (C) polymerase acidic protein, (D) hemagglutinin protein, (E) nucleoprotein, (F) neuraminidase protein, (G) matrix protein, and (H) non-structural protein. Red taxa marked with circles represent the H5N1 isolate from the leopard cat (Prionailurus bengalensis) in South Korea. Purple taxa indicate the H5N1 HPAI viruses isolated from wild birds during the 2024–2025 winter season. Squares denote HPAI viruses detected in wild birds in South Korea in October 2024. Background shading indicates HPAI isolates in blue. Numerical values at the nodes represent bootstrap support values (%) based on 1,000 replicates. Bootstrap values below 70% are not shown.

References

Summary

Keywords

HPAI H5N1, wild leopard cat, South Korea, clade 2.3.4.4b, mammalian adaptation, zoonotic potential

Citation

Si Y-J, Lee S-H, Kim D-J, Lee K, Lee M, Lee D-Y, Seo Y-R, Jeong H, Lee S and Lee D-H (2025) First detection of clade 2.3.4.4b H5N1 highly pathogenic avian influenza virus in a wild leopard cat (Prionailurus bengalensis) in South Korea. Front. Vet. Sci. 12:1638067. doi: 10.3389/fvets.2025.1638067

Received

30 May 2025

Accepted

30 June 2025

Published

29 July 2025

Volume

12 - 2025

Edited by

Iryna Goraichuk, Agricultural Research Service (USDA), United States

Reviewed by

Victor C. Huber, University of South Dakota, United States

Irene Iglesias, National Institute for Agricultural and Food Research and Technology, Spain

Judith Oguzie, University of Texas Medical Branch, United States

Updates

Copyright

*Correspondence: Dong-Hun Lee

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics