Next Article in Journal
Autoimmunity and Infection in Glomerular Disease
Next Article in Special Issue
Prevalence and Molecular Characteristics of FAdV-4 from Indigenous Chicken Breeds in Yunnan Province, Southwestern China
Previous Article in Journal
Discovery of Antimicrobial Peptides That Can Accelerate Culture Diagnostics of Slow-Growing Mycobacteria Including Mycobacterium tuberculosis
Previous Article in Special Issue
Environmental Sampling Methods for Detection of Salmonella Infections in Laying Hens: A Systematic Review and Meta-Analysis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 11
Issue 9
10.3390/microorganisms11092226
microorganisms-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

Emergence and Persistent Circulation of Highly Pathogenic Avian Influenza Virus A (H5N8) in Kosovo, May 2021–May 2022

by
Armend Cana
1,2,*,
Bianca Zecchin
3,
Xhavit Merovci
1,
Alice Fusaro
3,
Edoardo Giussani
3,
Sadik Heta
1,
Kiril Krstevski
4,
Dafina Mehmetukaj
1,
Izedin Goga
5,
Beqe Hulaj
1,
Bafti Murati
1,
Calogero Terregino
3 and
Aleksandar Dodovski
4
1
Kosovo Food and Veterinary Agency, Industrial Zone, 10 000 Prishtina, Kosovo
2
UBT—Higher Education Institution, Lagjja Kalabria, 10 000 Prishtina, Kosovo
3
Istituto Zooprofilattico Sperimentale Delle Venezie (IZSVe), 35020 Legnaro, Italy
4
Veterinary Institute, Faculty of Veterinary Medicine in Skopje, Ss. Cyril and Methodius University in Skopje, Lazar Pop Trajkov 5-7, MK-1000 Skopje, North Macedonia
5
Agricultural and Veterinary Faculty, University of Prishtina, Bulevardi Bill Clinton, 10 000 Prishtina, Kosovo
*
Author to whom correspondence should be addressed.
Microorganisms 2023, 11(9), 2226; https://doi.org/10.3390/microorganisms11092226
Submission received: 4 August 2023 / Revised: 28 August 2023 / Accepted: 29 August 2023 / Published: 2 September 2023
(This article belongs to the Special Issue Poultry Pathogens and Poultry Diseases)
Browse Figures
Review Reports Versions Notes

Abstract

:
In this study, we report the first outbreak of highly pathogenic avian influenza (HPAI) A H5N8, clade 2.3.4.4b in Kosovo on 19 May 2021. The outbreak consisted of three phases: May–June 2021, September–November 2021, and January–May 2022. In total, 32 backyards and 10 commercial holdings tested positive for the virus. Interestingly, the third and last phase of the outbreak coincided with the massive H5N1 clade 2.3.4.4b epidemic in Europe. Phylogenetic analyses of 28 viral strains from Kosovo revealed that they were closely related to the H5N8 clade 2.3.4.4.b viruses that had been circulating in Albania, Bulgaria, Croatia, Hungary, and Russia in early 2021. Whole genome sequencing of the 25 and partial sequencing of three H5N8 viruses from Kosovo showed high nucleotide identity, forming a distinctive cluster and suggesting a single introduction. The results of the network analysis were in accordance with the three epidemic waves and suggested that the viral diffusion could have been caused by secondary spreads among farms and/or different introductions of the same virus from wild birds. The persistent circulation of the same virus over a one-year period highlights the potential risk of the virus becoming endemic, especially in settings with non-adequate biosecurity.

1. Introduction

Avian Influenza (AI) A viruses are a diverse group of viruses causing mild to severe disease in poultry, pigeons, wild birds, and aquatic birds. The natural reservoirs for most low pathogenicity avian influenza (LPAI) viruses are wild birds, particularly birds of wetlands and aquatic environments such as the Anseriformes (particularly ducks, geese, and swans) and Charadriiformes (particularly gulls, terns, and waders) [1]. Highly pathogenic avian influenza (HPAI) viruses, with the potential for severe clinical signs and high mortality rates among poultry, originate from the changes in the hemagglutinin proteolytic cleavage site of H5 or H7 LPAI viruses. HPAI had not been present in wild bird host reservoirs before 2005 [2]. In 1996, a new H5 HPAI strain A/goose/Guangdong/1/1996 (Gs/GD) appeared for the first time in poultry in China, where it circulated for several years and spread to the surrounding countries. However, in 2005, this virus jumped from domestic to wild migratory birds. Since then, the H5 virus of the Gs/GD lineage has spread worldwide, causing mass mortality events in wild birds, great losses in the poultry industry and intermittent human cases [3,4,5]. While co-circulating with LPAI strains, it evolved into multiple genetic clades and variants with significant antigenic diversity [6]. A novel clade 2.3.4.4 reassortant H5N8, known as group A, first appeared in China in 2010, spread through Eurasia, and reached Europe in autumn 2014 [7,8]. In May 2016, an HPAI virus of subtype H5N8, clade 2.3.4.4 of group B (2.3.4.4b) was discovered in waterfowl and shorebirds in Siberia. As of March 2017, the virus had spread across Europe, the Middle East, and Africa [9,10,11,12], causing the most devastating epizootic ever recorded in domestic poultry [13]. Since then, multiple incursions of clade 2.3.4.4b have been reported in Europe in 2019–2020, 2020–2021, and 2021–2022 [14]. Increased deaths among a wide range of wild bird species, which were observed during the 2016–2017 HPAI H5N8 outbreaks in Europe, have shifted the paradigm of wild birds as unaffected agents of HPAI viruses, with increasing concerns about the potential effects on their populations as well as on the biodiversity of the ecosystems [15]. With a total of more than 2300 outbreaks in poultry, 46 million birds culled at the affected premises, and 2700 HPAI occurrences in wild birds throughout 36 European countries, the 2021–2022 HPAI epidemic season has been the largest pandemic to date to be documented in Europe [16]. The virus soon reached the Balkans [13,17,18] through wild bird migratory routes, the Black Sea–Mediterranean Flyway and the Adriatic Flyway. The first outbreak in Kosovo was identified on 19th May, 2021 [19]. Kosovo is a landlocked country in Southeast Europe with an area of 10,887 km2. It has a high density of backyard poultry farms as well as a high population of Anseriformes (geese, and domestic ducks), thus representing a favorable ecosystem for the introduction and spread of the infection throughout the country. Kosovo has around 120 commercial poultry farms distributed throughout the country, which host around 1.3 million birds. There are no data on the number of backyard farms. In this study, we report the first outbreak of HPAI A H5N8, clade 2.3.4.4b in Kosovo and the persistent circulation of the virus during the period May 2021–May 2022. Here, we also provide the most common necropsy findings observed, and to investigate the possible transmission dynamics among the affected farms, the complete genome sequences of 25 viruses were subjected to a more in-depth study by using a genetic network analysis.

2. Materials and Methods

The sampling frame was part of the passive surveillance annual program of the Kosovo Food and Veterinary Agency for AI and Newcastle disease (ND). Official veterinary inspectors visited all the suspected farms. An epidemiological investigation was carried out. The number and composition of poultry species on the farms and the number of dead poultry at the time of the first visit were recorded. Selected fresh poultry carcasses, as well as tracheal swabs, were submitted to the Kosovo Food and Veterinary Laboratory (KFVL) in the shortest time possible. From the first detected positive H5 Avian influenza case on 19 May 2021 and right up to May 2022, 502 poultry carcasses from 77 backyard farms and 12 commercial farms were submitted to the KFVL. The breakdown by species was as follows: 459 chickens/broilers, 8 turkeys, 22 common quails, 8 domestic ducks, 2 geese, and 3 pigeons. On average, 4.9 carcasses from backyard farms and 9.5 carcasses from commercial farms were collected. In addition, a variable number of oropharyngeal swabs from live poultry were collected from commercial farms only. All animals underwent gross necropsy under appropriate biosafety measures. For further analyses, oropharyngeal and cloacal swabs and tissues such as the lung, trachea, proventriculus, intestine, spleen, and brain were sampled.

2.1. Laboratory Analysis

Oropharyngeal swabs from each carcass were processed individually, in 1 X PBS containing 10,000 IU/mL penicillin, 10 mg/mL streptomycin, 0.25 mg/mL gentamicin, and 5000 IU/mL nystatin. Viral RNA from swabs was extracted manually using the QIAamp Viral RNA Mini Kit (Qiagen®, Hilden, Germany), according to the manufacturer’s recommendations, whereas extraction from tissues was done manually, starting with the rapid partial disruption of the tissue with PBS supplemented with antibiotics, and continuing with the IndiSpin® Pathogen Kit (Indical Bioscience, Leipzig, Germany). One to two individual tissues/organs representative from each farm were tested. Lungs were the preferred choice, whereas the second choice was proventricle. They were selected based on the presence of extensive lesions. The rest of the organs were kept frozen and stored for further analyses. All samples were tested by real-time RT-PCR for the AI matrix gene [20] and for H5 [21], H7 [22], N1 [23], and N8 [24] genes. Probes and primers were synthesized by Eurofins Genomics (Ebersberg, Germany).
Amplification of target sequences was performed using AgPath-ID One-Step RT-PCR kit (Applied Biosystems™, Waltham, MA, USA) by adding 12.5 μL 2× RT-PCR master mix reagent, 1 μL 25× RT-PCR enzyme mix; the primer and probe concentration for each assay is listed in Supplementary Table S1. Five µL of RNA template for AI type A and N1/N8 subtyping or 2 µL of RNA template for H5 and H7 subtyping were used. The 25 µL final reaction volume was adjusted with nuclease-free water up to the final volume for each assay. All reactions were performed on the QuantStudio 5 Real-Time PCR system (Applied Biosystems™, USA). The following thermal profile was used: reverse transcription at 45 °C for 10 min, activation step at 95 °C for 10 min, and 45 amplification cycles with denaturation at 95 °C for 15 s and annealing at 60 °C for AIV type A, and 56 °C for 45 s for H5, H7, N1, and N8. Samples were tested in the same run by using the Veryflex option.

2.1.1. H5 Sanger Sequencing of the Cleavage Site and HPAI Identification

For case definition, selected H5N8-positive samples were submitted to the Faculty of Veterinary Medicine—Skopje, Ss Cyril and Methodius University in Skopje for further characterization of H5 AI virus by One-step RT-PCR and for Sanger sequencing of the hemagglutinin cleavage site according to the protocol [25]. To increase the probability of obtaining sequences of good quality, RNA samples with a high (RNA) viral load were selected. The viral RNA yield in the purified RNA samples was assessed based on the Ct values of the samples obtained with the AIV M-gene RT qPCR.

2.1.2. Whole Genome Sequencing, Phylogenetic and Network Analyses

To confirm the results and further characterize viral strains, HPAI H5N8-positive representative samples from 10 farms sampled in May–June 2021, as well as from 17 farms sampled from September 2021 to April 2022, were sent to the EU Reference Laboratory for Avian Influenza and Newcastle Disease, Istituto Zooprofilattico Sperimentale delle Venezie (IZSVe, Legnaro, Italy). These samples consisted of oropharyngeal swabs and pools of organs such as the lung, trachea, proventriculus, intestine, and spleen from animals based on the low Ct values of the samples obtained previously with the AIV M-gene RT-qPCR. Total RNA was purified by using the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The SuperScript™ III One-Step RT-PCR System with Plati-num™ Taq High Fidelity DNA Polymerase (Invitrogen, Carlsbard, CA, USA) was used to obtain complete genomes as previously described by Zhou et al., 2009 [26]. The Agencourt AMPure XP (Beckman Coulter Inc., Brea, CA, USA) and Qubit™ DNA HS Assay (Thermo Fisher Scientific, Waltham, MA, USA) were, respectively, used to purify and quantify the amplicons that were then mixed in equimolar proportion. Illumina Nextera XT DNA Sample Preparation Kit (Illumina, San Diego, CA, USA) was used to prepare sequencing libraries. The Illumina MiSeq platform (2 × 250 bp Paired-End; Illumina, San Diego, CA, USA) was employed to perform the sequencing. FastQC v0.11.2 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 14 April 2022)) was used to assess the read quality and raw data were filtered by removing reads with more than 100 bases and with a Q score lower than 7, reads with more than 10% of undetermined bases, and duplicated paired-end reads. Illumina Nextera XT adaptors sequences (Illumina, San Diego, CA, USA) were clipped from reads with scythe v0.991 (https://github.com/vsbuffalo/scythe (accessed on 14 April 2022)) and trimmed with sickle v1.33 (https://github.com/najoshi/sickle (accessed on 14 April 2022)). Complete genomes were generated with BWA v0.7.12 (https://github.com/lh3/bwa (accessed on 14 April 2022)) [27] through a reference-based approach and the alignments were processed with Picard-tools v2.1.0 (http://picard.sourceforge.net) and GATK v3.5 (https://github.com/moka-guys/gatk_v3.5 (accessed on 14 April 2022)) [27,28,29]. LoFreq v2.1.2 (https://github.com/CSB5/lofreq (accessed on 14 April 2022)) [30] was used to call Single Nucleotide Polymorphisms (SNPs). Consensus sequences were submitted to the GISAID EpiFlu™ database (http://www.gisaid.org (accessed on 15 April 2022)) under the accession numbers reported in Table S2. Sequences of each gene segment were aligned in MAFFT v7 [31] and compared to the most related sequences available in GISAID obtained from a BLAST search. Maximum likelihood phylogenetic trees were generated in IQTREE v1.6.6 (https://github.com/iqtree/iqtree1 (accessed on 15 April 2022)) performing an ultrafast bootstrap resampling analysis (with 1000 replications). Phylogenetic trees were visualized in FigTree v1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 15 April 2022)).
The genetic network was obtained using the Median Joining (MJ) method implemented in NETWORK 10.2.0.0 [32] for the concatenated gene segments of 25 H5N8 non-reassortant viruses from Kosovo (the complete genome was available for 25 out of 28 sequenced viruses, Table 1). The MJ network uses a maximum parsimony approach to reconstruct the relationships between highly similar sequences, displayed by nodes and links connecting the nodes.

2.2. Ethical Statement

During the sample collection from live birds, tracheal and cloacal swabs were collected by expert veterinarians in accordance with ethical standards and animal welfare requirements on the Directive of the European Council on the protection of animals used for scientific purposes (2010/63/EU).

3. Results

3.1. Distribution of Positive Farms

Out of 502 tested carcasses, 177 out of 459 chickens or broilers, 4/8 ducks, 2/2 geese, 2/8 turkeys, and 16/22 common quails were positive for the virus. Tissue samples with positive results on AIV M-gene RT-qPCR showed a lower mean Ct value, with a mean Ct value of 19.78 (n = 58, standard deviation (SD) = 3.84), compared to individual oropharyngeal swabs from carcasses with positive results, with a mean Ct value of 23.80 (n = 201, SD = 4.88).
In the period 19 May 2021–12 May 2022, a total of 32 backyards and 10 commercial farms tested positive for HPAI H5N8 (Figure 1). The outbreak consisted of three phases, i.e., May–June 2021, September–November 2021, and January–May 2022, with the periods in between free of reported cases in poultry. During the first wave (May–June 2021), 20 farms tested positive, 19 of which were backyard farms, and in 11 of them, waterfowl such as geese and domestic ducks were present. The period between the 14 June and 29 September was free of reported cases. The recurrence (September–November 2021) was characterized by the involvement of 7 out of 11 commercial chicken farms. The third wave (January–May 2022) consisted of 11 outbreaks, with 2 out of 11 registered in commercial farms. The overall number of poultry in all the affected farms amounted to 179,198. Breakdown by species was of 175,058 chickens, 492 turkeys, 553 geese, 378 domestic ducks, and 2717 common quails. The number of affected poultry in commercial farms only was of 165,899, representing roughly 12.7% of the entire commercial poultry industry in the country.

3.2. Necropsy Finding

The most notable gross lesions at necropsy included: cyanotic comb and legs, multifocal petechiae to ecchymoses of the periventricular mucosa and intestine, and inflammations of the trachea, whereas multifocal necrosis in the pancreas was one of the most striking and consistent features observed during necropsy (Figure 2).

3.3. Virus Identification and Pathotyping

Virus isolation in SPF embryonated eggs (according to the Diagnostic Manual Directive 2006/437/CE) and viral intravenous inoculum in SPF chickens (in compliance with the Diagnostic Manual Directive 2006/437/CE)(*) were performed at the EURL—IZSVe on a chicken sample from the first outbreak case. The Intravenous Pathogenicity Index was 3.00.

3.4. Phylogenetic and Network Analyses

We characterized either the complete (n = 25) or partial (n = 3) genome of representative HPAI H5N8 viruses collected from backyards (n = 19) and commercial farms (n = 9) during the three epidemic waves. All the viruses from Kosovo belonged to clade 2.3.4.4b and clustered with the viruses that had been circulating in Russia and Europe since the end of 2020; in particular, they showed the highest similarity to the HPAI H5N8 viruses detected in Albania, Bulgaria, Croatia, Hungary, and Russia in early 2021 (Figure 3). The phylogenetic analyses of the eight gene segments (Figure 3, Supplementary Figures S1–S7) showed that all the HPAI H5N8 viruses from Kosovo cluster together, which suggests one single virus introduction into the country.
To investigate the possible transmission dynamics among the affected farms, the complete genome sequences of 25 viruses were subjected to a more in-depth study by using a genetic network analysis (Figure 4). The viruses showed a clear clustering by epidemic wave. The low number of nucleotide differences between the viruses collected from some of the outbreaks in commercial or backyard farms might be indicative of viral spread among them. In particular, viruses 3124-03, 3124-04, 3124-05, and 3124-07 collected from commercial farms cluster together, showing from 0 to 8 nucleotide differences along the entire genome, and separately from the backyard cases. On the other hand, some of the viruses collected from the backyard farms are closely related to each other, such as viruses 5162-15, 5162-19, 5162-07, 5162-04 (2 to 7 nucleotide differences) or viruses 5162-01, 5162-02, 5162-16 (1 to 5 nucleotide differences), indicating a potential virus spread among them.

4. Discussion

In 2020–2021, HPAI H5N8 was the main strain circulating in Europe, which was largely replaced by the HPAI H5N1 strain in the period 2021–2022 [16,19]. Currently, H5N1 has become the most widespread throughout Europe, although it has not been detected in Kosovo so far. Furthermore, Kosovo was one of the few countries where the circulation of HPAI H5N8 persisted for two consecutive epidemics [16,33]. Of particular interest is that the period between the first and second waves, lasting from 14 June to 29 September, despite a very active awareness campaign, was free of reported cases. Part of this period was covered by the strict housing order, which included the closing of all live bird markets and the prohibition of the poultry movement. After the abolishment of the housing order on 13 July 2021, 25 days after the last case during the first phase, the housing order was re-established again on 30 September 2021, after the reoccurrence of AIV H5N8 in a commercial poultry farm on 29 September 2021. The movement of poultry and all markets involving live poultry were prohibited again. All actions were based on Administrative Instruction No. 2005/24 Against Zoonotic Disease Avian Influenza and the National Contingency Plan for AI. The second phase during September raised concern regarding the involvement of a high number of commercial farms. Based on phylogenetic analyses, there is a clear indication of the single introduction of the virus in the country, and it is not possible to exclude a persistent circulation of the virus within the country in the time elapsing between the three epidemic waves in un-sampled hosts (e.g., wild birds). Due to the lack of sequences from wild birds, it is not clear whether the viral diffusion during the outbreak phases was caused by different introductions from wild birds or by secondary spread among farms. The virus could have circulated undetected between phases in resident wild birds or in the backyard poultry sector. However, farm-to-farm transmission of the virus, as suggested by the small number of nucleotide differences among some of the viruses recovered from various outbreaks during the same phase, is also supported by epidemiological investigation. For instance, the farmer of the outbreak on the commercial farm identified as 3124-03 had sold chickens or shared cages with farms identified as 3124-04, 3124-05, and 3124-07, confirming the hypothesis of secondary spread.
To date, no proper surveillance of the wild bird population has been carried out in Kosovo. At this step, it is of crucial importance for the authorities to start implementing active surveillance of bird flu in wild bird populations, as this is important to better evaluate the number of introductions in the poultry farms, as we cannot exclude the circulation in the wild population of highly related viruses. Moreover, surveillance in wild birds can be crucial to determine the possible source of infection for commercial poultry (backyard farms or wild birds, and if wild birds, which species?). This would help to improve surveillance and monitoring strategies. This publication covers the devastating impact that HPAI H5N8 had on the domestic backyard and commercial farms in Kosovo during the three outbreak phases (in 2020–2021 and 2021–2022). Twelve out of 120 chicken commercial farms in the country were affected in the period May 2021–May 2022. They represent around 12.7% of the entire commercial poultry population of the country. The number of backyard farms that may have been affected could be considerably higher than the one officially reported; in particular, those with a small number of poultry. Retrospectively, in identified HPAI H5N8 backyard farm outbreaks, epidemiological investigation revealed that there were sometimes unreported suspected cases within the same village. In fact, backyard farms could have played a significant role in the uncontrolled dissemination of the virus, and are more difficult to control. The persistent circulation of the same virus strain in the period under investigation highlights the potential risk of the virus to become endemic in the future, especially in settings with non-adequate biosecurity [34]. Considering the current pattern of the spread of H5Nx through the Black Sea-Mediterranean Flyway, it is likely that in the future, Kosovo might be exposed to different novel strains of AI. The great number of affected commercial farms highlights the urgent need to improve biosecurity measures and to prevent as much as possible contacts between wild birds and poultry.
The increasing number of reports of mammalian species infected by clade 2.3.4.4b as well as the sporadic human cases reported in recent years highlight the need to constantly monitor the evolution of this virus to promptly detect the emergence of new variants with possible increased zoonotic potential [35,36,37,38,39,40,41,42].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms11092226/s1, Supplementary Table S1. Primer/probe sequences used in the initial screening of the suspected outbreaks in the Kosovo Food and Veterinary Laboratory; Supplementary Table S2. Acknowledgment table of the authors, originating and submitting laboratories of the sequences from GISAID’s EpiFlu™ Database on which this work is based in part; all submitters of data may be contacted directly via www.gisaid.org (accessed on 15 April 2022); Supplementary Figure S1. Maximum likelihood phylogenetic tree of the PB2 gene segment, obtained in IQTREE V1.6.6. Viruses from Kosovo are colored in red; ultrafast bootstrap values higher than 80 are shown next to the nodes; Supplementary Figure S2. Maximum likelihood phylogenetic tree of the PB1 gene segment, obtained in IQTREE V1.6.6; Supplementary Figure S3. Maximum likelihood phylogenetic tree of the PA gene segment, obtained in IQTREE V1.6.6; Supplementary Figure S4. Maximum likelihood phylogenetic tree of the NP gene segment, obtained in IQTREE V1.6.6; Supplementary Figure S5. Maximum likelihood phylogenetic tree of the NA gene segment (N8 subtype), obtained in IQTREE V1.6.6; Supplementary Figure S6. Maximum likelihood phylogenetic tree of the M gene segment, obtained in IQTREE V1.6.6; Supplementary Figure S7. Maximum likelihood phylogenetic tree of the NS gene segment, obtained in IQTREE V1.6.6.

Author Contributions

Conceptualization, A.C., I.G., S.H. and B.M.; formal analysis, X.M., D.M. and B.H.; investigation, B.Z., A.C., E.G., X.M. and S.H., data curation, B.Z. and E.G.; writing—original draft preparation, A.C.; writing—review and editing, A.F., K.K., A.D. and C.T.; visualization, B.Z., A.C. and I.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The consensus sequences of the viruses analysed in this study were submitted to the GISAID EpiFlu™ database under the accession numbers reported in Table 1.

Acknowledgments

We are grateful to the Official Veterinarians, who gave great support in the epidemiological investigation activities and data collection, Laboratory staff of the Kosovo Food and Veterinary Laboratory, Istituto Zooprofilattico Sperimentale delle Venezie, and Veterinary Institute—Faculty of Veterinary Medicine in Skopje. We thank Francesca Ellero, Tony Wilsmore, and Malbora Shandro Zeqiri for their excellent assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fouchier, R.A.; Munster, V.J. Epidemiology of low pathogenic avian influenza viruses in wild birds. Rev. Sci. Tech. Off. Int. Epiz. 2009, 28, 49–58. [Google Scholar] [CrossRef] [PubMed]
  2. Lee, D.H.; Bertran, K.; Kwon, J.H.; Swayne, D.E. Evolution, global spread, and pathogenicity of highly pathogenic avian influenza H5Nx clade 2.3.4.4. J. Vet. Sci. 2019, 18, 269–280. [Google Scholar] [CrossRef] [PubMed]
  3. Lee, D.H.; Criado, M.F.; David, F.S. Pathobiological Origins and Evolutionary History of Highly Pathogenic Avian Influenza Viruses. Cold Spring Harb. Perspect. Med. 2021, 11, a038679. [Google Scholar] [CrossRef]
  4. WHO (World Health Organization). Cumulative Number of Confirmed Human Cases for Avian Influenza A(H5N1) Reported to WHO, 2003–2021, 15 April 2021. Available online: https://www.who.int/publications/m/item/cumulative-number-of-confirmed-human-cases-for-avian-influenza-a(h5n1)-reported-to-who-2003-2021-15-april-2021 (accessed on 15 April 2021).
  5. The Global Consortium for H5N8 and Related Influenza Viruses. Role for migratory wild birds in the global spread of avian influenza H5N8. Science 2016, 354, 213–217. [Google Scholar] [CrossRef]
  6. Baek, Y.G.; Lee, Y.N.; Lee, D.H.; Shin, J.I.; Lee, J.H.; Chung, D.H.; Lee, E.-K.; Heo, G.-B.; Sagong, M.; Kye, S.J.; et al. Multiple Reassortants of H5N8 Clade 2.3.4.4b Highly Pathogenic Avian Influenza Viruses Detected in South Korea during the Winter of 2020–2021. Viruses 2021, 13, 490. [Google Scholar] [CrossRef] [PubMed]
  7. Kwon, J.H.; Lee, D.H.; Swayne, D.E.; Noh, J.-Y.; Yuk, S.S.; Erdene-Ochir, T.O.; Hong, W.T.; Jeong, J.H.; Jeong, S.; Gwon, G.B.; et al. Highly Pathogenic Avian Influenza A(H5N8) Viruses Reintroduced into South Korea by Migratory Waterfowl, 2014–2015. Emerg. Infect. Dis. 2016, 22, 507–510. [Google Scholar] [CrossRef]
  8. Fan, S.; Zhou, L.; Wu, D.; Gao, X.; Pei, E.; Wang, T.; Gao, Y.; Xia, X. A novel highly pathogenic H5N8 avian influenza virus isolated from a wild duck in China. Influenza Other Respir. Viruses 2014, 8, 646–653. [Google Scholar] [CrossRef]
  9. Lee, D.H.; Sharshov, K.; Swayne, D.E.; Kurskaya, O.; Sobolev, I.; Kabilov, M.; Alekseev, A.; Irza, V.; Shestopalov, A. Novel Reassortant Clade 2.3.4.4 Avian Influenza A(H5N8) Virus in Wild Aquatic Birds, Russia, 2016. Emerg. Infect. Dis. 2017, 23, 359–360. [Google Scholar] [CrossRef]
  10. Koethe, S.; Ulrich, L.; Ulrich, R.; Amler, S.; Graaf, A.; Harder, T.C.; Grund, C.; Mettenleiter, T.C.; Conraths, F.J.; Beer, M.; et al. Modulation of lethal HPAIV H5N8 clade 2.3.4.4B infection in AIV pre-exposed mallards. Emerg. Microbes Infect. 2020, 9, 180–193. [Google Scholar] [CrossRef]
  11. Świętoń, E.; Fusaro, A.; Shittu, I.; Niemczuk, K.; Zecchin, B.; Joannis, T.; Bonfante, F.; Śmietanka, K.; Terregino, C. Sub-Saharan Africa and Eurasia Ancestry of Reassortant Highly Pathogenic Avian Influenza A(H5N8) Virus, Europe, December 2019. Emerg. Infect. Dis. 2020, 26, 1557–1561. [Google Scholar] [CrossRef]
  12. Brown, I.; Mulatti, P.; Smietanka, K.; Staubach, C.; Willeberg, P.; Adlhoch, C.; Candiani, D.; Fabris, C.; Zancanaro, G.; Morgado, J.; et al. Avian influenza overview October 2016–August 2017. EFSA J. 2017, 15, 5018. [Google Scholar]
  13. Adlhoch, C.; Fusaro, A.; Kuiken, T.; Monne, I.; Smietanka, K.; Staubach, C.; Guajardo, M.I.; Baldinelli, F. Avian influenza overview February–August 2019. EFSA J. 2019, 17, 5843. [Google Scholar]
  14. Adlhoch, C.; Fusaro, A.; Gonzales, J.L.; Kuiken, T.; Marangon, S.; Mirinaviciute, G.; Niqueux, É.; Stahl, K.; Staubach, C.; Terregino, C.; et al. Avian influenza overview December 2022–March 2023. EFSA J. 2023, 21, 07917. [Google Scholar]
  15. Kleyheeg, E.; Slaterus, R.; Bodewes, R.; Rijks, J.M.; Spierenburg, M.A.; Beerens, N.; Kelder, L.; Poen, M.P.; Stegeman, J.A.; Fouchier, R.A.; et al. Deaths among Wild Birds during Highly Pathogenic Avian Influenza A(H5N8) Virus Outbreak, the Netherlands. Emerg. Infect. Dis. 2017, 23, 2050–2054. [Google Scholar] [CrossRef]
  16. Adlhoch, C.; Fusaro, A.; Gonzales, J.L.; Kuiken, T.; Marangon, S.; Niqueux, É.; Staubach, C.; Terregino, C.; Aznar, I.; Guajardo, I.M.; et al. Avian influenza overview March–June 2022. EFSA J. 2022, 20, 7415. [Google Scholar]
  17. Zecchin, B.; Goujgoulova, G.; Monne, I.; Salviato, A.; Schivo, A.; Slavcheva, L.; Pastori, A.; Brown, I.H.; Lewis, N.S.; Terregino, C.; et al. Evolutionary Dynamics of H5 Highly Pathogenic Avian Influenza Viruses (Clade 2.3.4.4B) Circulating in Bulgaria in 2019–2021. Viruses 2021, 13, 2086. [Google Scholar] [CrossRef]
  18. Stoimenov, G.M.; Goujgoulova, G.V.; Nikolov, B.; Hristov, K.; Teneva, A. Pathological changes in natural infection of pheasants with highly pathogenic avian influenza A (H5N8) in Bulgaria. J. Vet. Res. 2019, 63, 497–502. [Google Scholar] [CrossRef]
  19. Adlhoch, C.; Fusaro, A.; Gonzales, J.L.; Kuiken, T.; Marangon, S.; Niqueux, É.; Staubach, C.; Terregino, C.; Aznar, I.; Guajardo, I.M.; et al. Avian influenza overview May–September 2021. EFSA J. 2021, 20, e07122. [Google Scholar]
  20. Heine, H.G.; Foord, A.J.; Wang, J.; Valdeter, S.; Walker, S.; Morrissy, C.; Wong, F.Y.; Meehan, B. Detection of highly pathogenic zoonotic influenza virus H5N6 by reverse-transcriptase quantitative polymerase chain reaction. Virol. J. 2015, 8, 12–18. [Google Scholar] [CrossRef]
  21. Slomka, M.J.; Pavlidis, T.; Banks, J.; Shell, W.; McNally, A.; Essen, S.; Brown, I.H. Validated H5 Eurasian real-time reverse transcriptase-polymerase chain reaction and its application in H5N1 outbreaks in 2005–2006. Avian Dis. 2007, 51, 373–377. [Google Scholar] [CrossRef]
  22. Van Borm, S.; Suarez, D.L.; Boschmans, M.; Ozhelvaci, O.; Marché, S.; van den Berg, T.P. Rapid detection of Eurasian and American H7 subtype influenza A viruses using a single TaqManMGB real-time RT-PCR. Avian Dis. 2010, 54, 632–638. [Google Scholar] [CrossRef] [PubMed]
  23. Payungporn, S.; Chutinimitkul, S.; Chaisingh, A.; Damrongwantanapokin, S.; Buranathai, C.; Amonsin, A.; Theamboonlers, A.; Poovorawan, Y. Single step multiplex real-time RT-PCR for H5N1 influenza A virus detection. J. Virol. Methods 2006, 131, 143–147. [Google Scholar] [CrossRef] [PubMed]
  24. James, J.; Slomka, M.J.; Reid, S.M.; Thomas, S.S.; Mahmood, S.; Byrne, A.M.P.; Cooper, J.; Russell, C.; Mollett, B.C.; Agyeman-Dua, E.; et al. Development and Application of Real-Time PCR Assays for Specific Detection of Contemporary Avian Influenza Virus Subtypes N5, N6, N7, N8, and N9. Avian Dis. 2019, 6, 209–218. [Google Scholar] [CrossRef] [PubMed]
  25. Slomka, M.J.; Coward, V.J.; Banks, J.; Löndt, B.Z.; Brown, I.H.; Voermans, J.; Koch, G.; Handberg, K.J.; JØrgensen, P.H.; Cherbonnel-Pansart, M.; et al. Identification of sensitive and specific avian influenza polymerase chain reaction methods through blind ring trials organized in the European Union. Avian Dis. 2007, 51, 227–234. [Google Scholar] [CrossRef]
  26. Zhou, B.; Donnelly, M.E.; Scholes, D.T.; George, K.S.; Hatta, M.; Kawaoka, Y.; Wentworth, D.E.; St. George, K.; Hatta, M.; Kawaoka, Y.; et al. Single-Reaction Genomic Amplification Accelerates Sequencing and Vaccine Production for Classical and Swine Origin Human Influenza A Viruses. J. Virol. 2009, 83, 10309–10313. [Google Scholar] [CrossRef]
  27. Van der Auwera, G.; Carneiro, M.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ data to high confidence variant calls: The Genome Analysis Toolkit best practices pipeline. Curr. Protoc. Bioinform. 2013, 43, 11.10.1–11.10.33. [Google Scholar] [CrossRef]
  28. Depristo, M.A.; Banks, E.; Poplin, R.E.; Garimella, K.V.; Maguire, J.R.; Hartl, C.; Philippakis, A.A.; Del Angel, G.; Rivas, M.A.; Hanna, M.; et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 2011, 43, 491–498. [Google Scholar] [CrossRef]
  29. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef]
  30. Wilm, A.; Aw, P.P.K.; Bertrand, D.; Yeo, G.H.T.; Ong, S.H.; Wong, C.H.; Khor, C.C.; Petric, R.; Hibberd, M.L.; Nagarajan, N. LoFreq: A sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic. Acids Res. 2012, 40, 11189–11201. [Google Scholar] [CrossRef]
  31. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef]
  32. Bandelt, H.J.; Forster, P.; Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef] [PubMed]
  33. Adlhoch, C.; Fusaro, A.; Gonzales, J.L.; Kuiken, T.; Marangon, S.; Niqueux, É.; Staubach, C.; Terregino, C.; Aznar, I.; Guajardo, I.M.; et al. Avian influenza overview September–December 2021. EFSA J. 2021, 19, e7108. [Google Scholar]
  34. El-Shesheny, R.; Kandeil, A.; Mostafa, A.; Ali, M.A.; Webby, R.J. H5 Influenza Viruses in Egypt. Cold Spring Harb. Perspect. Med. 2021, 11, a038745. [Google Scholar] [CrossRef] [PubMed]
  35. Pyankova, O.G.; Susloparov, I.M.; Moiseeva, A.A.; Kolosova, N.P.; Onkhonova, G.S.; Danilenko, A.V.V.; Vakalova, E.V.; Shendo, G.L.; Nekeshina, N.N.; Noskova, L.N.; et al. Isolation of clade 2.3.4.4b A(H5N8), a highly pathogenic avian influenza virus, from a worker during an outbreak on a poultry farm, Russia, December 2020. Eurosurveillance 2021, 26, 2100439. [Google Scholar] [CrossRef]
  36. WHO (World Health Organization). Avian Influenza A(H5N8) Infects Humans in Russian Federation; WHO: Copenhagen, Denmark, 2021. [Google Scholar]
  37. Floyd, T.; Banyard, A.C.; Lean, F.Z.X.; Byrne, A.M.P.; Fullick, E.; Whittard, E.; Mollett, B.C.; Bexton, S.; Swinson, V.; Macrelli, M.; et al. Encephalitis and Death in Wild Mammals at a Rehabilitation Center after Infection with Highly Pathogenic Avian Influenza A(H5N8) Virus, United Kingdom. Emerg. Infect. Dis. 2021, 27, 2856–2863. [Google Scholar] [CrossRef] [PubMed]
  38. Herve, S.; Schmitz, A.; Briand, F.X.; Gorin, S.; Quéguiner, S.; Niqueux, É.; Paboeuf, F.; Scoizec, A.; Le Bouquin-Leneveu, S.; Eterradossi, N.; et al. Serological Evidence of Backyard Pig Exposure to Highly Pathogenic Avian Influenza H5N8 Virus during 2016–2017 Epizootic in France. Pathogens 2021, 10, 621. [Google Scholar] [CrossRef]
  39. Schülein, A.; Ritzmann, M.; Christian, J.; Schneider, K.; Neubauer-Juric, A. Exposure of wild boar to Influenza A viruses in Bavaria: Analysis of seroprevalences and antibody subtype specificity before and after the panzootic of highly pathogenic avian influenza viruses A (H5N8). Zoonoses Public Health 2021, 68, 503–515. [Google Scholar] [CrossRef]
  40. Adlhoch, C.; Fusaro, A.; Gonzales, J.L.; Kuiken, T.; Meliodu, A.; Mirinaviciute, G.; Niqueux, É.; Stahl, K.; Staubach, C.; Terregino, C.; et al. Avian influenza overview April–June 2023. EFSA J. 2023, 21, e08191. [Google Scholar]
  41. Fusaro, A.; Monne, I.; Mulatti, P.; Zecchin, B.; Bonfanti, L.; Ormelli, S.; Milani, A.; Cecchettin, K.; Lemey, P.; Moreno, A.; et al. Genetic Diversity of Highly Pathogenic Avian Influenza A(H5N8/H5N5) Viruses in Italy, 2016–2017. Emerg. Infect. Dis. 2017, 23, 1543–1547. [Google Scholar] [CrossRef]
  42. Shi, W.; Gao, G.F. Emerging H5N8 avian influenza viruses. Science 2021, 372, 784–786. [Google Scholar] [CrossRef]
Microorganisms 11 02226 g001
Figure 1. Spatial distribution and incidence timeline of HPAI H5N8-positive farms in Kosovo in the period May 2021–May 2022 (red color), and timeline distribution of suspected AI farms that resulted negative (light blue).
Figure 1. Spatial distribution and incidence timeline of HPAI H5N8-positive farms in Kosovo in the period May 2021–May 2022 (red color), and timeline distribution of suspected AI farms that resulted negative (light blue).
Microorganisms 11 02226 g001
Microorganisms 11 02226 g002
Figure 2. Gross lesions of chicken: (A,B) Cyanotic comb and legs, (C,D) Petechiae to ecchymoses of the proventricular mucosa and intestine, (E) Inflammations of the trachea and (F) multifocal necrosis in the pancreas.
Figure 2. Gross lesions of chicken: (A,B) Cyanotic comb and legs, (C,D) Petechiae to ecchymoses of the proventricular mucosa and intestine, (E) Inflammations of the trachea and (F) multifocal necrosis in the pancreas.
Microorganisms 11 02226 g002
Microorganisms 11 02226 g003
Figure 3. Maximum likelihood phylogenetic tree of the HA gene segment, obtained in IQTREE V1.6.6. Viruses from Kosovo are colored in red; ultrafast bootstrap values higher than 80 are shown next to the nodes.
Figure 3. Maximum likelihood phylogenetic tree of the HA gene segment, obtained in IQTREE V1.6.6. Viruses from Kosovo are colored in red; ultrafast bootstrap values higher than 80 are shown next to the nodes.
Microorganisms 11 02226 g003
Microorganisms 11 02226 g004
Figure 4. Genetic network generated using the Median Joining method implemented in NETWORK 10.2.0, for the eight concatenated gene segments of 25 non-reassortant H5N8 viruses from Kosovo. Circles represent viral types; the size is proportional to the number of viruses sharing the same type, and the length of the branches is proportional to the number of nucleotide differences. Samples are colored according to the month of collection. Commercial farms are identified by underlined IZSVe-IDs.
Figure 4. Genetic network generated using the Median Joining method implemented in NETWORK 10.2.0, for the eight concatenated gene segments of 25 non-reassortant H5N8 viruses from Kosovo. Circles represent viral types; the size is proportional to the number of viruses sharing the same type, and the length of the branches is proportional to the number of nucleotide differences. Samples are colored according to the month of collection. Commercial farms are identified by underlined IZSVe-IDs.
Microorganisms 11 02226 g004
Table 1. HPAI H5N8 viruses from Kosovo sequenced within the framework of the present study.
Table 1. HPAI H5N8 viruses from Kosovo sequenced within the framework of the present study.
VirusDistrictCategoryIDIZSVe-IDDateGenomeAccession Number
A/partridge/Kosovo/96_21VIR5162-14/2021PrizrenBackyard965162-1431 May 2021CompleteEPI_ISL_3128524
A/chicken/Kosovo/97_21VIR5162-15/2021GjilanBackyard975162-1531 May 2021CompleteEPI_ISL_3128525
A/chicken/Kosovo/98_21VIR5162-16/2021MitrovicaBackyard985162-1631 May 2021CompleteEPI_ISL_3128526
A/duck/Kosovo/107_21VIR5162-17/2021FerizajBackyard1075162-1702 June 2021CompleteEPI_ISL_3128527
A/chicken/Kosovo/113_21VIR5162-19/2021FerizajBackyard1135162-1903 June 2021CompleteEPI_ISL_3128528
A/chicken/Kosovo/82_21VIR5162-1/2021MitrovicaBackyard825162-119 May 2021CompleteEPI_ISL_3128529
A/chicken/Kosovo/84_21VIR5162-2/2021FerizajBackyard845162-225 May 2021CompleteEPI_ISL_3128530
A/chicken/Kosovo/90_21VIR5162-4/2021PrizrenBackyard905162-425 May 2021CompleteEPI_ISL_3128531
A/chicken/Kosovo/104_21VIR5162-7/2021PrizrenBackyard1045162-701 June 2021CompleteEPI_ISL_3128532
A/chicken/Kosovo/279_22VIR3124-2/2021PrizrenCommercial2793124-204 October 2021CompleteEPI_ISL_12176839
A/chicken/Kosovo/280_22VIR3124-3/2021GjakovaCommercial2803124-305 October 2021CompleteEPI_ISL_12176840
A/chicken/Kosovo/283_22VIR3124-4/2021FerizajCommercial2833124-407 October 2021CompleteEPI_ISL_12176841
A/chicken/Kosovo/284_22VIR3124-5/2021GjakovaCommercial2843124-511 October 2021CompleteEPI_ISL_12176843
A/common_quail/Kosovo/288_22VIR3124-6/2021PrizrenBackyard2883124-618 October 2021CompleteEPI_ISL_12176844
A/chicken/Kosovo/290_22VIR3124-7/2021PristinaCommercial2903124-721 October 2021CompleteEPI_ISL_12176845
A/chicken/Kosovo/299_22VIR3124-8/2021FerizajCommercial2993124-803 November 2021CompleteEPI_ISL_12176846
A/chicken/Kosovo/302_22VIR3124-9/2021PristinaBackyard3023124-905 November 2021CompleteEPI_ISL_12176847
A/chicken/Kosovo/303_22VIR3124-10/2021PristinaBackyard3033124-1005 November 2021CompleteEPI_ISL_12176848
A/chicken/Kosovo/22-2_22VIR3124-13/2022PristinaBackyard22-23124-1306 January 2022CompleteEPI_ISL_12176850
A/chicken/Kosovo/22-8_22VIR3124-14/2022PristinaBackyard22-83124-1414 January 2022CompleteEPI_ISL_12176851
A/chicken/Kosovo/22-9_22VIR3124-15/2022GjakovaCommercial22-93124-1519 January 2022CompleteEPI_ISL_12176852
A/chicken/Kosovo/126_22VIR3124-18/2021MitrovicaCommercial1263124-1814 June 2021CompleteEPI_ISL_12176853
A/chicken/Kosovo/22-50_22VIR3124-19/2022PrizrenBackyard22-503124-1915 February 2022CompleteEPI_ISL_12176854
A/chicken/Kosovo/22-59_22VIR3124-20/2022PristinaBackyard22-593124-2001 March 2022CompleteEPI_ISL_12176855
A/turkey/Kosovo/13-2_22VIR3124-31/2022PrizrenBackyard13-23124-3119 January 2022CompleteEPI_ISL_12176931
A/turkey/Kosovo/90_21VIR5162-5/2021PrizrenBackyard905162-525 May 2021PartialEPI_ISL_3142407
A/chicken/Kosovo/278_22VIR3124-1/2021PristinaCommercial2783124-129 September 2021PartialEPI_ISL_12176838
A/common_quail/Kosovo/304-2_22VIR3124-12/2021PrizrenBackyard304-23124-1208 November 2021PartialEPI_ISL_12176849
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cana, A.; Zecchin, B.; Merovci, X.; Fusaro, A.; Giussani, E.; Heta, S.; Krstevski, K.; Mehmetukaj, D.; Goga, I.; Hulaj, B.; et al. Emergence and Persistent Circulation of Highly Pathogenic Avian Influenza Virus A (H5N8) in Kosovo, May 2021–May 2022. Microorganisms 2023, 11, 2226. https://doi.org/10.3390/microorganisms11092226

AMA Style

Cana A, Zecchin B, Merovci X, Fusaro A, Giussani E, Heta S, Krstevski K, Mehmetukaj D, Goga I, Hulaj B, et al. Emergence and Persistent Circulation of Highly Pathogenic Avian Influenza Virus A (H5N8) in Kosovo, May 2021–May 2022. Microorganisms. 2023; 11(9):2226. https://doi.org/10.3390/microorganisms11092226

Chicago/Turabian Style

Cana, Armend, Bianca Zecchin, Xhavit Merovci, Alice Fusaro, Edoardo Giussani, Sadik Heta, Kiril Krstevski, Dafina Mehmetukaj, Izedin Goga, Beqe Hulaj, and et al. 2023. "Emergence and Persistent Circulation of Highly Pathogenic Avian Influenza Virus A (H5N8) in Kosovo, May 2021–May 2022" Microorganisms 11, no. 9: 2226. https://doi.org/10.3390/microorganisms11092226

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