Evolutionary Dynamics of H5 Highly Pathogenic Avian Influenza Viruses (Clade 2.3.4.4B) Circulating in Bulgaria in 2019–2021
by
, ,
Bianca Zecchin
1,*
,
Gabriela Goujgoulova
2,
Isabella Monne
1,
Annalisa Salviato
1,
Alessia Schivo
1,
Iskra Slavcheva
2,
Ambra Pastori
1,
Ian H. Brown
3,
Nicola S. Lewis
3,4,
Calogero Terregino
1 and
Alice Fusaro
1,* 1
EU/OIE/National Reference Laboratory for Avian Influenza and Newcastle Disease, FAO Reference Centre for Animal Influenza and Newcastle Disease, Istituto Zooprofilattico Sperimentale delle Venezie, 35020 Legnaro, Italy
2
National Reference Laboratory of Avian Influenza and Newcastle Disease, National Diagnostic and Research Veterinary Medical Institute, 1231 Sofia, Bulgaria
3
OIE/FAO International Reference Laboratory for Avian Influenza, Swine Influenza and Newcastle Disease Virus, Animal and Plant Health Agency-Weybridge, Addlestone, Surrey KT15 3NB, UK
4
Department of Pathobiology and Population Sciences, Royal Veterinary College, Hatfield, Hertfordshire AL9 7TA, UK
*
Authors to whom correspondence should be addressed.
Viruses 2021, 13(10), 2086; https://0-doi-org.brum.beds.ac.uk/10.3390/v13102086
Submission received: 20 August 2021
/
Revised: 22 September 2021
/
Accepted: 11 October 2021
/
Published: 16 October 2021
(This article belongs to the Special Issue Drivers of Evolution of Animal RNA Viruses)
Abstract
:The first detection of a Highly Pathogenic Avian Influenza (HPAI) H5N8 virus in Bulgaria dates back to December 2016. Since then, many outbreaks caused by HPAI H5 viruses from clade 2.3.4.4B have been reported in both domestic and wild birds in different regions of the country. In this study, we characterized the complete genome of sixteen H5 viruses collected in Bulgaria between 2019 and 2021. Phylogenetic analyses revealed a persistent circulation of the H5N8 strain for four consecutive years (December 2016–June 2020) and the emergence in 2020 of a novel reassortant H5N2 subtype, likely in a duck farm. Estimation of the time to the most recent common ancestor indicates that this reassortment event may have occurred between May 2019 and January 2020. At the beginning of 2021, Bulgaria experienced a new virus introduction in the poultry sector, namely a HPAI H5N8 that had been circulating in Europe since October 2020. The periodical identification in domestic birds of H5 viruses related to the 2016 epidemic as well as a reassortant strain might indicate undetected circulation of the virus in resident wild birds or in the poultry sector. To avoid the concealed circulation and evolution of viruses, and the risk of emergence of strains with pandemic potential, the implementation of control measures is of utmost importance, particularly in duck farms where birds display no clinical signs.
1. Introduction
In 2016, HPAI H5 viruses of clade 2.3.4.4B descending from the Gs/GD lineage caused one of the biggest epidemics in Europe by the number of countries affected and the species diversity of infected wild birds [1]. Different HPAI H5 subtypes were detected after the occurrence of reassortment events with LPAI viruses circulating in wild birds. Within the European territory, alongside the HPAI H5N8, which represented the main circulating subtype, HPAI H5N2, H5N5 and H5N6 were also identified [1]. The latter was subsequently responsible for a new epidemic wave in north-eastern Europe between December 2017 and January 2019 [2,3].
In December 2019 a new incursion of a HPAI H5N8 subtype was reported in Poland [4]. Since then, the virus has spread to 8 central-eastern European countries (Bulgaria, Czech Republic, Germany, Hungary, Poland, Romania, Slovakia, Ukraine), affecting mainly the poultry population, with more than 300 outbreaks notified in Hungary and Poland [5,6,7,8]. The last outbreak was reported in June 2020.
In October 2020, with the beginning of the migration of wild birds to their wintering areas, a new HPAI H5 incursion occurred in northern Europe. Similar to the 2016 epidemic wave, the virus had first been detected further east and arrived across a wide area causing variable mortality in the affected wild bird populations and rapidly spread to most of the European countries [9,10,11,12].
In Bulgaria, the first detection of a HPAI H5N8 virus dates back to December 2016. Since then, many outbreaks have been reported in domestic and wild birds in different administrative regions [13]. A recent work by Venkatesh et al. (2020) describes the 2017/2018 HPAI H5N8 epidemic in Bulgaria. The study reports the identification of two separate viral introductions, one in the north-eastern region of Dobrich and another in central and eastern Bulgaria, and suggests that the duck sector played a critical role in HPAI H5N8 virus spread and maintenance in the country [9]. Here, we describe the genetic diversity of HPAI H5 viruses of clade 2.3.4.4B in Bulgaria, analyzing sixteen newly sequenced H5 viruses identified in the country between 2019 and 2021 and the emergence of a new H5N2 subtype, likely from a local reassortment event.
2. Materials and Methods
2.1. Sample Collection
Between March 2019 and February 2021, sixteen outbreaks caused by the HPAI H5N8 and H5N2 subtypes were registered in chicken and duck farms in Plovdiv, Lovech and Pleven regions (Table 1).
The National Diagnostic and Research Veterinary Medical Institute of Sofia (Bulgaria) sent clinical samples (intestines, tracheal and cloacal swabs collected from H5Nx infected mule ducks, layer hens and domestic ducks) as well as infectious allantoic fluids (harvested after inoculating 9–11 day old embryonated specific pathogen free (SPF) hen’s eggs for virus isolation) to the EU/OIE/National Reference Laboratory for Avian Influenza and Newcastle Disease at the Istituto Zooprofilattico Sperimentale delle Venezie in Legnaro (Italy) for genetic characterization.
2.2. Whole Genome Sequencing
Total RNA was extracted using QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Complete genomes were obtained as previously described using the SuperScript™ III One-Step RT-PCR System with Platinum™ Taq High Fidelity DNA Polymerase (Invitrogen, Carlsbard, CA, USA) [14]. Amplicons were purified using Agencourt AMPure XP (Beckman Coulter Inc., Brea, CA, USA), quantified with Qubit™ DNA HS Assay (Thermo Fisher Scientific, Waltham, MA, USA), and mixed in equimolar proportion. Illumina Nextera XT DNA Sample Preparation Kit (Illumina, San Diego, CA, USA) was used to prepare sequencing libraries. Sequencing was performed on the Illumina MiSeq platform (2x250 bp Paired-End; Illumina, San Diego, CA, USA). The read quality was assessed by using FastQC v0.11.2 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/, accessed on 4 March 2021) and raw data were filtered by removing reads with more than 100 bases with Q score below seven, reads with more than 10% of undetermined (“N”) 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 4 March 2021) and trimmed with sickle v1.33 (https://github.com/najoshi/sickle, accessed on 4 March 2021). Complete genomes were generated through a reference-based approach using BWA v0.7.12 (https://github.com/lh3/bwa, accessed on 4 March 2021) [15] and processing the alignments with Picard-tools v2.1.0 (http://picard.sourceforge.net, accessed on 4 March 2021) and GATK v3.5 (https://github.com/moka-guys/gatk_v3.5, accessed on 4 March 2021) [15,16,17]. LoFreq v2.1.2 (https://github.com/CSB5/lofreq, accessed on 4 March 2021) [18] was used to call Single Nucleotide Polymorphisms (SNPs). The consensus sequences were submitted to the GISAID EpiFlu™ database (http://www.gisaid.org, accessed on 26 April 2021) under the accession numbers EPI1780066-EPI1780081, EPI1807268-EPI1807356 (Table 1).
2.3. Phylogenetic and Evolutionary Analyses
Sequences were aligned in MAFFT v7 [19]. Maximum likelihood phylogenetic trees were generated in IQTREE v1.6 (https://github.com/iqtree/iqtree1, accessed on 2 September 2021) [20,21], performing an ultrafast bootstrap resampling analysis (1000 replications). Phylogenetic trees were visualized in FigTree v1.4.2 (http://tree.bio.ed.ac.uk/software/figtree/, accessed on 2 September 2021).
The evolutionary rate and the time to the most common recent ancestor (tMRCA) for all the eight gene segments were estimated using Bayesian inference. Markov chain Monte Carlo (MCMC) sampling analyses were performed using BEAST v1.10.4 [22]. We employed an uncorrelated lognormal relaxed molecular clock that allows for rate variation across lineages. We used the SRD06 substitution model (HKY85 + Γ4 with two partitions - 1st + 2nd positions vs. 3rd position -, base frequencies and Γ-rate heterogeneity unlinked across all codon positions) [23]. Maximum clade credibility (MCC) trees were summarized using TreeAnnotator v1.10.4 (http://beast.bio.ed.ac.uk/TreeAnnotator/, accessed on 2 September 2021) and visualized in FigTree v1.4.2 (http://tree.bio.ed.ac.uk/software/figtree/, accessed on 2 September 2021).
3. Results
In this study, we analyzed the genetic characteristics of the HPAI H5N8 and H5N2 viruses which caused sixteen outbreaks notified in Bulgaria since the beginning of 2019, and specifically in the periods March-April 2019, February-June 2020 and February 2021. The phylogenetic analysis of the hemagglutinin (HA) gene showed that the HPAI H5N8 viruses collected in March-April 2019 (outbreaks 19/1-19/5) and the H5N8 (outbreaks 20/2, 20/5, 20/7) and H5N2 (outbreaks 20/3, 20/6, 20/8) viruses collected in the period February-June 2020 are related to HPAI H5N8 viruses of clade 2.3.4.4B, which have been identified in Bulgaria since the end of 2016. Despite the introduction in December 2019 in east-central Europe of a novel HPAI H5N8, which extensively circulated in Poland, Germany, Hungary, Czech Republic and Slovakia during the first half of 2020, the outbreaks reported in Bulgaria in the same time period were not caused by the strain newly introduced in Europe through wintering migrations of wild birds (Figure 1).
The first outbreak of 2020 (20/1) occurred in a mule duck fattening farm. The clinical sample and the respective isolate were both confirmed to be co-infected with H5N2 and H5N8 viral subtypes by Real Time PCR (using specific primers both for the N2 and for the N8 subtypes), as well as sequencing. Specifically, next generation sequencing (Illumina) of the clinical sample only determined an H5N2 subtype, which was the prevalent one, while a partial sequence of about 300 nucleotides of the N8 segment was obtained using Sanger sequencing. Differently, next generation sequencing of the isolate clearly showed a mixed H5N2/H5N8 virus population and the complete sequence of both N2 and N8 genes was obtained.
The analyses of the complete genomes confirm that the 2019–2020 H5N8 belong to the genotype responsible of the 2016 European epidemic (Figure 2), clustering with HPAI H5N8 viruses circulating in Bulgaria in 2017–2018, for all the eight gene segments. Specifically, they can be grouped into three different sub-clusters: the first one includes the HPAI H5N8 viruses collected in 2019 (outbreaks 19/1 to 19/5); the second one is composed of HPAI H5N8 viruses from outbreaks 20/2 (ducks) and 20/5 (mule ducks) from Plovdiv region; the third one contains HPAI H5N8 viruses from outbreaks 20/7 (hens) and 20/9 (hens), respectively from Kurdzhali and Plovdiv regions.
The H5N2 viruses, all identified in Plovdiv region from outbreaks 20/1 (H5N2/N8 co-infected mule ducks), 20/3 (hens), 20/6 (hens) and 20/8 (hens), derive from reassortment events between the HPAI H5N8 and LPAI viruses. Specifically, the HA and M genes (Figure 1 and Figure S7) cluster with HPAI H5N8 viruses circulating in Bulgaria in 2017–2020, showing the highest identity (99.1%–99.3% for the HA gene) with viruses of sub-cluster 3 (Figure 1 and Figures S1–S8); while the N2 gene clusters with a LPAI H4N2 virus identified in ducks in Bulgaria in 2018 (95.5%–95.6% identity) (Supplementary Figure S5). The internal gene segments of the H5N2 subtype show the highest similarity with different LPAI viruses identified in wild and domestic birds in Eurasia (Supplementary Figures S1–S4 and S8). In particular, the NS gene segment of the Bulgarian H5N2 viruses belongs to allele B, while the H5N8 viruses possess the allele A (Supplementary Figure S8). Of note, a sample co-infected with H5N8 and H5N2 viral subtypes had already been identified in ducks in Plovdiv region in 2018 (accession numbers EPI1743963-EPI1743964, GISAID EpiFlu database) but its N2 gene is very different from that of the 2020 H5N2 viruses (Supplementary Figure S5), indicating an independent reassortment event.
To estimate the time of emergence of the HPAI H5N8 and H5N2 viruses in Bulgaria we performed a Bayesian evolutionary analysis of all the eight gene segments, and we estimated the time to the most common recent ancestor (tMRCA). The HPAI H5N8 viruses belonging to the cluster of the 2018–2020 Bulgarian viruses seem to have emerged in the first half of 2017 (Dec. 2016–Aug. 2017) (Figure 3 and Figures S9–S18). Conversely, the reassortant H5N2 viruses collected in Bulgaria in 2020 appear to have been generated between May 2019 and January 2020, with the exception of NS gene segment, for which the limited number of allele B viruses widens the tMRCA estimate to between October 2017 and July 2019 (Figure 3 and Figures S9–S18). The 2020 reassortant H5N2 strain was identified only in duck and chicken farms located in Plovdiv region, suggesting that it likely emerged in this area, which was hardly affected by the H5N8 subtype outbreaks.
In February 2021 a new introduction of HPAI H5N8 viruses of clade 2.3.4.4B was recorded in the northern region of Pleven. The 2021 HPAI H5N8 viruses (outbreaks 21/1-21/3, Table 1) cluster separately from the viruses previously detected in the country and show the highest identity with the HPAI H5N8 viruses, which have been circulating in Europe since October 2020 [12]; in particular the three Bulgarian viruses are almost identical and cluster with HPAI H5N8 viruses identified in Poland, England and Russia (Krasnodar Region) in November 2020 and January 2021 (Figure 1 and Figure 2).
Molecular analysis of the HA protein revealed that the Bulgarian H5N2 viruses obtained from outbreaks 20/1 and 20/3 possess a potential additional glycosylation site at position 70–72 (D70N) (numbering starting from the first methionine), which likely indicates an adaptation of the virus to domestic birds [24]. For the NA protein, potential additional glycosylation sites were identified in A/chicken/Bulgaria/217_20VIR1724-1/2020 (H5N8) from outbreak 20/7 at position 411-413 (E411N). No evidence of known mutations associated with mammalian adaptation were recorded in any of the analysed viruses.
4. Discussion
In Europe, the last cases of HPAI H5N8 viruses related to the 2016–2017 epidemic were detected in March and July 2018 in Italy and Russia, respectively [25,26]. After the detection of the HPAI H5N8 in 2016, Bulgaria has continued to identify up to June 2020 H5 viruses related to the 2016–2017 European epidemic in domestic birds. This may indicate an undetected circulation of the virus in resident wild birds, or more likely in poultry sector, in particular in duck farms where the HPAI H5 viruses have shown to be able to circulate causing no or few clinical signs [27,28]. Extremely low biosecurity measures in foie gras farms and persistent circulation of a single strain in this poultry sector have been previously reported in the country [29]. Several studies based on epidemiological and phylogenetic analyses have recognized domestic Anseriformes to play a critical role in HPAI virus persistence, spread and evolution as well as being an important source for virus infection in chickens [30,31,32,33,34].
After France, Bulgaria is one of the main European foie gras producers, with a highly developed commercial poultry sector which includes mule ducks and hybrid ducks raised for foie gras and meat production. Following recovery from the 2016–2017 avian influenza crisis, production underwent a sharp increase in 2018. In 2019, the Bulgarian industry produced about 2600 tonnes of foie gras intended for export mainly to France (71%), Belgium (14%) and Spain (14%) (https://www.itavi.asso.fr/download/10930, accessed on 10 October 2021, Situation de la production et du marché du foie gras en avril 2020, accessed on 20 September 2021). This rapid expansion, not always accompanied by the implementation of risk-appropriate biosecurity measures with only very basic sanitation and hygiene, could represent the ideal environment for avian influenza viruses to be introduced, maintained undetected and reassort, as demonstrated by the identification of two independently originated H5N2 reassortant strains in 2018 and 2020 in one of the regions with the highest duck population density [35].
Of note, the 2020 H5N2 was firstly identified in a farm rearing mule duck, where the H5N8 was also detected. Additional testing would have been useful to elucidate the extent of circulation of the two subtypes in the farm. The subsequent identification of a genetically related H5N2 strain in hens (outbreaks 20/3, 20/6 and 20/8) may indicate a virus transmission between the two poultry sectors, which may have been fostered by the use of vehicles interchangeably employed to transport eggs, fodder and animals from/to different production holdings. However, we cannot exclude that the virus might have circulated undetected also in other poultry sectors (i.e., backyards). The estimation of tMRCA indicated that this subtype may have been generated between the end of 2019 and the beginning of 2020, which clearly indicates that this genotype was detected soon after its emergence, about two years after the introduction in Bulgaria of the H5N8 progenitor strain (first half of 2017 based on tMRCA estimation). However, given the lack of identification of the LPAI progenitor virus, it is difficult to determine whether the reassortment event took place within the mule ducks farm or in another setting (i.e., in another farm or in wild birds). Increasing active surveillance in wild birds may help to shed light on the ecological context which may favor the emergence of reassortant viruses.
The new introduction in February 2021 of a novel H5N8 strain, which has been showing a high propensity to reassort and the ability to acquire mutations associated with mammalian adaptation, is of concern [36]. The implementation of adequate control measures as well as a regular and capillary surveillance in the poultry industry, in particular in the duck farms, also in absence of evident clinical signs is of utmost importance to stop the circulation of the virus in the country. This could help preventing the emergence of novel reassortant viruses, or of variants with unexpected biological properties, which may represent a serious risk for the operators in the poultry sectors, but also for the general public health.
Supplementary Materials
The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/v13102086/s1, Figure S1: Maximum Likelihood phylogenetic tree of the PB2 gene, Figure S2: Maximum Likelihood phylogenetic tree of the PB1 gene, Figure S3: Maximum Likelihood phylogenetic tree of the PA gene, Figure S4: Maximum Likelihood phylogenetic tree of the NP gene, Figure S5: Maximum Likelihood phylogenetic tree of the NA (N2 subtype) gene, Figure S6: Maximum Likelihood phylogenetic tree of the NA (N8 subtype) gene, Figure S7: Maximum Likelihood phylogenetic tree of the M gene, Figure S8: Maximum Likelihood phylogenetic tree of the NS gene, Figure S9: Maximum Clade Credibility Tree of the PB2 gene, Figure S10: Maximum Clade Credibility Tree of the PB1 gene, Figure S11: Maximum Clade Credibility Tree of the PA gene, Figure S12: Maximum Clade Credibility Tree of the HA gene, Figure S13: Maximum Clade Credibility Tree of the NP gene, Figure S14: Maximum Clade Credibility Tree of the NA-N2 subtype gene, Figure S15: Maximum Clade Credibility Tree of the NA-N8 subtype gene, Figure S16: Maximum Clade Credibility Tree of the M gene, Figure S17: Maximum Clade Credibility Tree of the NS (Allele A) gene, Figure S18: Maximum Clade Credibility Tree of the NS (Allele B) gene.
Author Contributions
Conceptualization, A.F., B.Z., I.M.; methodology, A.F., B.Z.; formal analysis, A.P., A.S. (Annalisa Salviato), B.Z., A.S. (Alessia Schivo); resources, G.G., I.S.; writing—original draft preparation, B.Z., A.F.; writing—review and editing, A.F., I.M., B.Z., C.T., G.G., N.S.L., I.H.B.; supervision, C.T., A.F. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
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 wish to thank Francesca Ellero for her excellent technical assistance.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Alarcon, P.; Brouwer, A.; Venkatesh, D.; Duncan, D.; Dovas, C.I.; Georgiades, G.; Monne, I.; Fusaro, A.; Dan, A.; Śmietanka, K.; et al. Comparison of 2016–17 and Previous Epizootics of Highly Pathogenic Avian Influenza H5 Guangdong Lineage in Europe. Emerg. Infect. Dis. 2018, 24, 2270–2283. [Google Scholar] [CrossRef] [PubMed]
- Baldinelli, F.; Papanikolaou, A.; Stoicescu, A.; Van der Stede, Y.; Aznar, I. Annual Report on surveillance for Avian Influenza in poultry and wild birds in Member States of the European Union in 2019. EFSA J. 2020, 18, 6349. [Google Scholar] [CrossRef]
- Brouwer, A.; Gonzales, J.; Huneau, A.; Mulatti, P.; Kuiken, T.; Staubach, C.; Stegeman, A.; Antoniou, S.E.; Baldinelli, F.; Van der Stede, Y.; et al. Annual Report on surveillance for avian influenza in poultry and wild birds in Member States of the European Union in 2018. EFSA J. 2019, 17. [Google Scholar] [CrossRef]
- Śmietanka, K.; Świętoń, E.; Kozak, E.; Wyrostek, K.; Tarasiuk, K.; Tomczyk, G.; Konopka, B.; Welz, M.; Domańska-Blicharz, K.; Niemczuk, K. Highly pathogenic avian influenza H5N8 in Poland in 2019–2020. J. Vet. Res. 2020, 64, 469–476. [Google Scholar] [CrossRef]
- Adlhoch, C.; Fusaro, A.; Kuiken, T.; Niqueux, E.; Staubach, C.; Terregino, C.; Guajardo, I.M.; Baldinelli, F. Avian influenza overview November 2019–February2020. EFSA J. 2020, 18. [Google Scholar] [CrossRef]
- Adlhoch, C.; Fusaro, A.; Kuiken, T.; Niqueux, E.; Staubach, C.; Terregino, C.; Guajardo, I.M.; Baldinelli, F. Avian influenza overview February–May 2020. EFSA J. 2020, 18. [Google Scholar] [CrossRef]
- Adlhoch, C.; Fusaro, A.; Kuiken, T.; Niqueux, É.; Staubach, C.; Terregino, C.; Muñoz Guajardo, I.; Baldinelli, F. Avian influenza overview May–August 2020. EFSA J. 2020, 18. [Google Scholar] [CrossRef]
- Lewis, N.S.; Banyard, A.C.; Whittard, E.; Karibayev, T.; Al Kafagi, T.; Chvala, I.; Byrne, A.; Meruyert Akberovna, S.; King, J.; Harder, T.; et al. Emergence and spread of novel H5N8, H5N5 and H5N1 clade 2.3.4.4 highly pathogenic avian influenza in 2020. Emerg. Microbes. Infect. 2021, 10, 148–151. [Google Scholar] [CrossRef]
- Venkatesh, D.; Brouwer, A.; Ellis, R.; Goujgoulova, G.; Seekings, J.; Brown, I.; Lewis, N. Regional transmission and reassortment of 2.3.4.4b highly pathogenic avian influenza (HPAI) viruses in Bulgarian poultry 2017/18. Viruses 2020, 12, 605. [Google Scholar] [CrossRef]
- Adlhoch, C.; Fusaro, A.; Gonzales, J.L.; Kuiken, T.; Marangon, S.; Niqueux, É.; Staubach, C.; Smietanka, K.; Terregino, C.; Van der Stede, Y.; et al. Avian influenza overview—update on 19 November 2020, EU/EEA and the UK. EFSA J. 2020, 18, e06341. [Google Scholar] [CrossRef]
- Adlhoch, C.; Fusaro, A.; Gonzales, J.L.; Kuiken, T.; Marangon, S.; Niqueux, É.; Staubach, C.; Terregino, C.; Baldinelli, F. Avian influenza overview August–December 2020. EFSA J. 2020, 18, 6379. [Google Scholar] [CrossRef]
- Adlhoch, C.; Fusaro, A.; Gonzales, J.L.; Kuiken, T.; Marangon, S.; Niqueux, É.; Staubach, C.; Terregino, C.; Muñoz Guajardo, I.; Lima, E.; et al. Avian influenza overview December 2020–February 2021. EFSA J. 2021, 19, e06497. [Google Scholar] [CrossRef]
- 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] [Green Version]
- Zhou, B.; Donnelly, M.E.; Scholes, D.T.; St. George, K.; Hatta, M.; Kawaoka, Y.; Wentworth, D.E. 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] [PubMed] [Green Version]
- 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 Bioinforma. 2013, 43, 10.1–10.33. [Google Scholar] [CrossRef]
- 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]
- 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] [Green Version]
- 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] [Green Version]
- 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] [Green Version]
- Nguyen, L.-T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
- Hoang, D.T.; Chernomor, O.; von Haeseler, A.; Minh, B.Q.; Vinh, L.S. UFBoot2: Improving the Ultrafast Bootstrap Approximation. Mol. Biol. Evol. 2018, 35, 518–522. [Google Scholar] [CrossRef]
- Suchard, M.A.; Lemey, P.; Baele, G.; Ayres, D.L.; Drummond, A.J.; Rambaut, A. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 2018, 4, vey016. [Google Scholar] [CrossRef] [Green Version]
- Shapiro, B.; Rambaut, A.; Drummond, A.J. Choosing Appropriate Substitution Models for the Phylogenetic Analysis of Protein-Coding Sequences. Mol. Biol. Evol. 2006, 23, 7–9. [Google Scholar] [CrossRef] [Green Version]
- Matrosovich, M.; Zhou, N.; Kawaoka, Y.; Webster, R. The Surface Glycoproteins of H5 Influenza Viruses Isolated from Humans, Chickens, and Wild Aquatic Birds Have Distinguishable Properties. J. Virol. 1999, 73, 1146. [Google Scholar] [CrossRef] [Green Version]
- Adlhoch, C.; Brouwer, A.; Kuiken, T.; Mulatti, P.; Smietanka, K.; Staubach, C.; Muñoz Guajardo, I.; Verdonck, F.; Amato, L.; Baldinelli, F. Avian influenza overview February–May 2018. EFSA J. 2018, 16, e05358. [Google Scholar] [CrossRef]
- World Organisation for Animal Health (OIE). OIE Situation Report for Highly Pathogenic Avian Influenza; 2018; Available online: https://www.oie.int/fileadmin/Home/eng/Animal_Health_in_the_World/docs/pdf/OIE_AI_situation_report/OIE_SituationReport_AI_August2018.pdf (accessed on 10 October 2021).
- Kim, J.-K.; Negovetich, N.J.; Forrest, H.L.; Webster, R.G. Ducks: The “Trojan Horses” of H5N1 influenza. Influenza Other Respi. Viruses 2009, 3, 121–128. [Google Scholar] [CrossRef]
- Webster, R.G.; Webby, R.J.; Hoffmann, E.; Rodenberg, J.; Kumar, M.; Chu, H.J.; Seiler, P.; Krauss, S.; Songserm, T. The immunogenicity and efficacy against H5N1 challenge of reverse genetics-derived H5N3 influenza vaccine in ducks and chickens. Virology 2006, 351, 303–311. [Google Scholar] [CrossRef] [Green Version]
- Marinova-Petkova, A.; Georgiev, G.; Petkov, T.; Darnell, D.; Franks, J.; Kayali, G.; Walker, D.; Seiler, P.; Danner, A.; Graham, A.; et al. Influenza surveillance on “foie gras” duck farms in Bulgaria, 2008–2012. Influenza Other Respi. Viruses 2016, 10, 98–108. [Google Scholar] [CrossRef] [PubMed]
- Henning, J.; Wibawa, H.; Morton, J.; Usman, T.B.; Junaidi, A.; Meers, J. Scavenging Ducks and Transmission of Highly Pathogenic Avian Influenza, Java, Indonesia. Emerg. Infect. Dis. 2010, 16, 1244. [Google Scholar] [CrossRef] [PubMed]
- Gilbert, M.; Xiao, X.; Pfeiffer, D.U.; Epprecht, M.; Boles, S.; Czarnecki, C.; Chaitaweesub, P.; Kalpravidh, W.; Minh, P.Q.; Otte, M.J.; et al. Mapping H5N1 highly pathogenic avian influenza risk in Southeast Asia. Proc. Natl. Acad. Sci. USA 2008, 105, 4769–4774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Gilbert, M.; Chaitaweesub, P.; Parakamawongsa, T.; Premashthira, S.; Tiensin, T.; Kalpravidh, W.; Wagner, H.; Slingenbergh, J. Free-grazing Ducks and Highly Pathogenic Avian Influenza, Thailand. Emerg. Infect. Dis. 2006, 12, 227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Kwon, J.-H.; Lee, D.-H.; Criado, M.F.; Killmaster, L.; Ali, M.Z.; Giasuddin, M.; Samad, M.A.; Karim, M.R.; Hasan, M.; Brum, E.; et al. Genetic evolution and transmission dynamics of clade 2.3.2.1a highly pathogenic avian influenza A/H5N1 viruses in Bangladesh. Virus Evol. 2020, 6, veaa046. [Google Scholar] [CrossRef] [PubMed]
- Barman, S.; Marinova-Petkova, A.; Hasan, M.K.; Akhtar, S.; El-Shesheny, R.; Turner, J.C.; Franks, J.; Walker, D.; Seiler, J.; Friedman, K.; et al. Role of domestic ducks in the emergence of a new genotype of highly pathogenic H5N1 avian influenza A viruses in Bangladesh. Emerg. Microbes Infect. 2017, 6, e72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Stoimenov, G.M. Review Highly Pathogenic Avian Influenza in Bulgaria A Review. Bulg. J. Vet. Med. 2020. [Google Scholar] [CrossRef]
- 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 February–May 2021. EFSA J. 2021. Available online: https://www.ecdc.europa.eu/en/publications-data/avian-influenza-overview-february-may-2021#no-link (accessed on 10 October 2021).
Figure 1.
Maximum Likelihood phylogenetic tree of the HA gene (clade 2.3.4.4B). The tree was inferred using IQTREE v1.6.6. The grey box shows the cluster of Bulgarian H5 (N8 and N2) viruses identified in 2019–2020 together with some previously reported 2018 H5N8-H5N2 detections; the light green box shows the cluster of 2020 Bulgarian H5N2 viruses. The viruses are colored according to the years of collection. Ultrafast bootstrap supports higher than 80 are indicated next to the nodes.
Figure 1.
Maximum Likelihood phylogenetic tree of the HA gene (clade 2.3.4.4B). The tree was inferred using IQTREE v1.6.6. The grey box shows the cluster of Bulgarian H5 (N8 and N2) viruses identified in 2019–2020 together with some previously reported 2018 H5N8-H5N2 detections; the light green box shows the cluster of 2020 Bulgarian H5N2 viruses. The viruses are colored according to the years of collection. Ultrafast bootstrap supports higher than 80 are indicated next to the nodes.
Figure 2.
HPAI A(H5) genotypes identified in Bulgaria in 2019–2021. Blue bars: gene segments originating from Eurasian/African HPAI H5N8 viruses of the European epidemic started in 2016; violet bars: gene segments originating from Eurasian/African HPAI H5N8 viruses of the European epidemic of 2020–2021; pink bar: gene segment closely related to the LPAI virus A/mallard/Novosibirsk region/999k/2018; light-blue bars: gene segments closely related to the LPAI virus A/duck/Hungary/1147/2019; orange bar: gene segment closely related to the LPAI virus A/common teal/Chany/892/2018; yellow bar: gene segment closely related to the LPAI virus A/duck/Bulgaria_Plovdiv/75/2018; dark violet bar: gene segment closely related to viruses of which the NS gene belongs to Allele B.
Figure 2.
HPAI A(H5) genotypes identified in Bulgaria in 2019–2021. Blue bars: gene segments originating from Eurasian/African HPAI H5N8 viruses of the European epidemic started in 2016; violet bars: gene segments originating from Eurasian/African HPAI H5N8 viruses of the European epidemic of 2020–2021; pink bar: gene segment closely related to the LPAI virus A/mallard/Novosibirsk region/999k/2018; light-blue bars: gene segments closely related to the LPAI virus A/duck/Hungary/1147/2019; orange bar: gene segment closely related to the LPAI virus A/common teal/Chany/892/2018; yellow bar: gene segment closely related to the LPAI virus A/duck/Bulgaria_Plovdiv/75/2018; dark violet bar: gene segment closely related to viruses of which the NS gene belongs to Allele B.
Figure 3.
Time to the most common recent ancestor (tMRCA) calculated for each gene segment for the 2019–2020 Bulgarian H5N8 (grey box in Figure 1 and Figures S1–S8) and the H5N2 (light green box in Figure 1 and Figures S1–S8) subgroups.
Figure 3.
Time to the most common recent ancestor (tMRCA) calculated for each gene segment for the 2019–2020 Bulgarian H5N8 (grey box in Figure 1 and Figures S1–S8) and the H5N2 (light green box in Figure 1 and Figures S1–S8) subgroups.
Table 1.
Analyzed viruses, epidemiological information and GISAID EpiFlu™ database accession numbers.
Table 1.
Analyzed viruses, epidemiological information and GISAID EpiFlu™ database accession numbers.
| Outbreak | Virus | Subtype | Region | District | Collection Date (dd/mm/yyyy) | Species | Production Category | Accession Number (GISAID EpiFlu™) |
|---|---|---|---|---|---|---|---|---|
| 19/1 | A/duck/Bulgaria/130-1_19VIR3314-1/2019 | H5N8 | Lovech | Yoglav | 27/03/2019 | Mule duck | Fattening | EPI1807349-EPI1807356 |
| 19/2 | A/chicken/Bulgaria/136_19VIR3315-1/2019 | H5N8 | Plovdiv | Krumovo | 01/04/2019 | Chicken | Backyard | EPI1807268-EPI1807275 |
| 19/3 | A/chicken/Bulgaria/143_19VIR3315-2/2019 | H5N8 | Plovdiv | Asenograd | 04/04/2019 | Chicken | Laying hens | EPI1807276-EPI1807283 |
| 19/4 | A/chicken/Bulgaria/148-01_19VIR3315-3/2019 | H5N8 | Plovdiv | Asenograd | 05/04/2019 | Chicken | Laying hens | EPI1807284-EPI1807291 |
| 19/5 | A/chicken/Bulgaria/150-01_19VIR3315-5/2019 | H5N8 | Plovdiv | Asenograd | 05/04/2019 | Chicken | Laying hens | EPI1807292-EPI1807299 |
| 20/1 | A/mule_duck/Bulgaria/50-506_20VIR1414-1/2020 | H5N2/N8 | Plovdiv | Rakovski | 11/02/2020 | Mule duck | Fattening | EPI1807332-EPI1807339, EPI1807348 |
| 20/2 | A/duck/Bulgaria/78-4t_20VIR1416-3/2020 | H5N8 | Plovdiv | Padarsko | 21/02/2020 | Duck | Foie gras | EPI1807300-EPI1807307 |
| 20/3 | A/chicken/Bulgaria/77_20VIR1727/2020 | H5N2 | Plovdiv | Trilistnik | 21/02/2020 | Chicken | Layer | EPI1780067, EPI1780069, EPI1780071, EPI1780073, EPI1780075, EPI1780077, EPI1780079, EPI1780081 |
| 20/5 | A/mule_duck/Bulgaria/147_20VIR1721-1/2020 | H5N8 | Plovdiv | Bolyarino | 21/02/2020 | Mule duck | Foie gras | EPI1807340-EPI1807347 |
| 20/6 | A/chicken/Bulgaria/201_20VIR1723-1/2020 | H5N2 | Plovdiv | Trilistnik | 02/03/2020 | Chicken | Layer | EPI1807308-EPI1807315 |
| 20/7 | A/chicken/Bulgaria/217_20VIR1724-1/2020 | H5N8 | Kurdzhali | Perperek | 09/03/2020 | Chicken | Laying hens | EPI1780066, EPI1780068, EPI1780070, EPI1780072, EPI1780074, EPI1780076, EPI1780078, EPI1780080 |
| 20/8 | A/chicken/Bulgaria/221_20VIR1725-1/2020 | H5N2 | Plovdiv | Manole | 11/03/2020 | Chicken | Laying hens | EPI1807316-EPI1807323 |
| 20/9 | A/chicken/Bulgaria/380-1_20VIR3542-1/2020 | H5N8 | Plovdiv | Asenograd | 03/06/2020 | Chicken | Laying hens | EPI1807324-EPI1807331 |
| 21/1 | A/chicken/Bulgaria/39_21VIR1454-3/2021 | H5N8 | Pleven | Slavyanovo | 01/02/2021 | Chicken | Laying hens | EPI1858599-EPI1858606 |
| 21/2 | A/pekin_duck/Bulgaria/48-3_21VIR1454-7/2021 | H5N8 | Pleven | Slavyanovo | 04/02/2021 | Pekin duck | Fattening | EPI1858607-EPI1858614 |
| 21/3 | A/chicken/Bulgaria/50-1_21VIR1454-9/2021 | H5N8 | Pleven | Slavyanovo | 08/02/2021 | Chicken | Breeding hens | EPI1858615- EPI1858622 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Zecchin, B.; Goujgoulova, G.; Monne, I.; Salviato, A.; Schivo, A.; Slavcheva, I.; 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. https://0-doi-org.brum.beds.ac.uk/10.3390/v13102086
AMA Style
Zecchin B, Goujgoulova G, Monne I, Salviato A, Schivo A, Slavcheva I, Pastori A, Brown IH, Lewis NS, 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(10):2086. https://0-doi-org.brum.beds.ac.uk/10.3390/v13102086
Chicago/Turabian StyleZecchin, Bianca, Gabriela Goujgoulova, Isabella Monne, Annalisa Salviato, Alessia Schivo, Iskra Slavcheva, Ambra Pastori, Ian H. Brown, Nicola S. Lewis, Calogero Terregino, and et al. 2021. "Evolutionary Dynamics of H5 Highly Pathogenic Avian Influenza Viruses (Clade 2.3.4.4B) Circulating in Bulgaria in 2019–2021" Viruses 13, no. 10: 2086. https://0-doi-org.brum.beds.ac.uk/10.3390/v13102086
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
