|Year : 2021 | Volume
| Issue : 7 | Page : 309-315
Spatio-temporal history of H9N2 viruses in Iran and neighbor countries by Bayesian analysis and molecular characterization
Nima Ghalekhani1, Saied Bokaie2, Sana Eybpoosh3, Hesameddin Akbarein2
1 Divisions of Epidemiology & Zoonoses, Department of Food Hygiene and Quality Control, Faculty of Veterinary Medicine, University of Tehran; HIV/STI Surveillance Research Center, and WHO Collaborating Center for HIV Surveillance, Institute for Futures Studies in Health, Kerman University of Medical Sciences, Kerman, Iran
2 Divisions of Epidemiology & Zoonoses, Department of Food Hygiene and Quality Control, Faculty of Veterinary Medicine, University of Tehran, Iran
3 Department of Epidemiology and Biostatistics, Research Centre for Emerging and Reemerging Infectious Diseases, Pasteur Institute of Iran, Tehran, Iran
|Date of Submission||01-Apr-2021|
|Date of Decision||26-Jun-2021|
|Date of Acceptance||29-Jun-2021|
|Date of Web Publication||12-Jul-2021|
Divisions of Epidemiology & Zoonoses, Department of Food Hygiene and Quality Control, Faculty of Veterinary Medicine, University of Tehran
Source of Support: None, Conflict of Interest: None
Objective: To delineate the H9N2 influenza virus circulation within Iran and its neighboring countries, the potential source of the epidemic in these countries, and its date of origin.
Methods: We obtained all hemagglutinin (HA) and neuraminidase (NA) nucleotide sequences of influenza H9N2 available up to December 25, 2020 from Iran and its neighboring countries (i.e., Pakistan, Afghanistan, Turkmenistan, Armenia, Azerbaijan, Turkey, and Iraq). We also performed a Bayesian Markov chain Monte Carlo method to infer the evolutionary dynamic and the most recent common ancestor for the HA and NA sequences.
Results: H9N2 epidemic may have started in Iran and Pakistan much earlier than the other investigated countries in the region, and an ongoing bidirectional dispersion of the virus between the investigated countries was also observed. The mean time of the most recent common ancestor of H9N2 viruses was 1988 for HA, and 1992 for NA.
Conclusions: Strains from investigated countries rooted in Pakistan and Iran. Regular surveillance of H9N2 viruses, especially in the live bird markets, enhancing the biosecurity of poultry industry and screening newly arriving immigrants and tourists from neighboring countries at border should be considered to control spread of the virus. Furthermore, surveillance of viral molecular evolution should be initiated for effective prevention of epidemic and pandemic spreads.
Keywords: Influenza virus; H9N2; Bayesian phylogeographic; Iran
|How to cite this article:|
Ghalekhani N, Bokaie S, Eybpoosh S, Akbarein H. Spatio-temporal history of H9N2 viruses in Iran and neighbor countries by Bayesian analysis and molecular characterization. Asian Pac J Trop Med 2021;14:309-15
|How to cite this URL:|
Ghalekhani N, Bokaie S, Eybpoosh S, Akbarein H. Spatio-temporal history of H9N2 viruses in Iran and neighbor countries by Bayesian analysis and molecular characterization. Asian Pac J Trop Med [serial online] 2021 [cited 2021 Jul 24];14:309-15. Available from: https://www.apjtm.org/text.asp?2021/14/7/309/320521
| 1. Introduction|| |
Avian influenza is an important public health threat worldwide. Influenza A viruses are classified based on their hemagglutinin (HA) and neuraminidase (NA) surface proteins. So far, 18 HA and 11 NA subsets have been identified. Based on their pathogenicity in birds, these strains are classified into low and highly pathogenic avian influenza. The compatibility and survival of these viruses in waterfowl made these birds natural reservoirs of the viruses. Birds can transmit the virus to a wide range of mammals, such as pigs, horses, and humans. Among these strains, highly acute ones with H5 and H7 surface antigens and low pathogenicity strains with H9 have received more attention due to their potentials for causing pandemics through infected humans. In 1999 and following in 2003, H9N2 (A/HK/2108/03) and (99/A/HK/1073) viruses were reported to infect humans, and since then, the H9N2 virus has been considered as a potential candidate for the development of new pandemics. The H9N2 avian influenza virus is one of the major serotypes of influenza virus that has recently been distributed in different regions of the world.
Studies have shown that H9N2 is genetically similar to the H5N1 isolated from humans in 1997. Six genetic segments of H5N1 are shared with the H9N2 virus, which was found in a large number of Asian birds, and that virus had the potential to infect 18 people and kill six cases. H9N2 viruses are classified into two major lineages, i. e., the “North American” and “Eurasian” lineages, circulating in poultry and wild birds. The Eurasian lineage is subdivided into three major lineages: a) Korean lineage [A/duck/Hong Kong/Y439/97 (Y439-like) and A/chicken/Korea/38349-p96323/96 (Korean-like)]; b) Y280 lineage [A/duck/Hong Kong/Y280/97-like (Y280-like), A/Chicken/Hong Kong/G9/97 (G9-like) and A/Chicken/Beijing/1/94 (BJ94-like)]; and c) 3-G1 lineage [A/quail/Hong Kong/G1/97-like (G1-like)].
The H9N2 serotype was first isolated and identified in 1998 in Iran known as A/chicken/Iran/259/1998 in the Poultry Diseases Research and Diagnosis Department of Razi Institute and Faculty of Veterinary Medicine, University of Tehran. Unfortunately, despite the initial plan to remove infected mother flocks and eggs, hatching and vaccination with inactivated vaccines, the disease outbreak spread rapidly through the country and has become endemic. It has caused a lot of damage to the country’s poultry industry.
Simultaneously with the spread of viral diseases such as influenza, the Bayesian phylogeographic method in 2009 was initiated and developed. For example, Lemey et al. employed this method to investigate the global spread of influenza A (H5N1) and the rabies virus. Several studies have also estimated the Spatio-temporal spread of some HIV-1 clades in different countries,,.
Despite the H9N2 avian influenza virus was mostly circulated in Iran and most of its neighboring countries, such as Pakistan, Afghanistan, Turkmenistan, Azerbaijan, Turkey, and Iraq, over the past decade,,, little is known about the virus’ circulation within and between these countries, the potential source of the epidemic in these countries, or its date of origin. To shed more light on the above knowledge gaps, we performed this Bayesian phylogeographic study on H9N2 genetic sequences available in the region. The results may be helpful in designing evidence-based control measures and prevention of onward virus transmissions.
| 2. Materials and methods|| |
We downloaded all HA and NA nucleotide sequences of influenza H9N2 available up to December 25, 2020 from Iran and its neighboring countries (i.e., Pakistan, Afghanistan, Armenia, Turkmenistan, Azerbaijan, Turkey and Iraq) from the NCBI Influenza Virus Resource. Using information deposited in the NCBI Influenza Virus Resource, the data was annotated with sampling date and location, as well as host types (chickens, ducks, wild birds, etc.,. After the exclusion of identical sequences and isolates with unknown sampling time and location, a total of 268 HA and 211 NA sequences remained in our dataset. Multiple sequence alignment was performed using the ClustalW algorithm in MEGA v7.0 and manually inspected and trimmed for phylogenetic analysis. The initial maximum-likelihood phylogenetic analysis was conducted using Hasegawa-Kishino-Yano for HA and General Time Reversible for N2 nucleotide substitution models, as identified by an online execution on the ATGC bioinformatics platform (http://www.atgc-montpelier. fr/). Analysis was performed with 1 000 bootstraps in MEGA v7.0.
We also performed a Bayesian Markov chain Monte Carlo (MCMC) method to infer the evolutionary dynamic and the most recent common ancestor for the HA and NA sequences. Bayesian MCMC analysis was carried out using the BEAST v.2.5.241 package, with a Hasegawa-Kishino-Yano+G (HA sequences) and General Time Reversible+G+I (NA sequences) model as the substitution models along with an uncorrelated lognormal relaxed clock model and Bayesian skyline coalescent tree priors. Two independent MCMC chains were run for 50 million generations (sampling every 5 000 steps) and were combined using the LogCombiner program v.1.54. Convergence was assessed based on the effective sampling size (ESS) after a 10% burn-in using Tracer software version 1.5. ESSs of 200 and more were accepted and a minimum ESS of 200 was required. The maximum clade credibility tree was generated in the TreeAnnotator program, while the initial 10% of trees were discarded as burn-in,. The maximum clade credibility tree was visualized in the FigTree program version 1.2.3.
Inferences about spatial dynamic and potential viral migration patterns of H9N2 within Iran and between Iran and its neighboring countries were made based on a discrete Bayesian phylogeographic model developed in BEAST software v.2.5.241. The symmetric substitution model with the Bayesian Stochastic Search Variable Selection approach was performed in the SPREAD program v.1.0.7.
| 3. Results|| |
3.1. Phylogenetic characteristics of H9N2 virus
The topology of the phylogenetic tree for the HA region suggested five genetic lineages. In the NA phylogenetic tree, we identified two lineages [Supplementary Figure 1] [Additional file 1] and [Supplementary Figure 2] [Additional file 2]. The biggest lineage in the HA tree is labeled as lineage 1 and consists of strains with Pakistan, Afghanistan, Iran, and Iraq origins. Chicken was the major host of the isolates within this lineage. The Lineage 4 is comprised only of strains that belong to Iran, but other lineages are formed by Pakistan and Iran’s strains [Figure 1] and [Figure 2].
|Figure 1: Maximum likelihood phylogenetic tree of the hemagglutinin gene (Lineage 1).|
Click here to view
|Figure 2: Maximum likelihood phylogenetic tree of the hemagglutinin gene (Lineage 4).|
Click here to view
The NA tree showed that the main locations of viruses in lineage 1 and 2 were Pakistan and Iran and for the NA tree chicken was the major host of the strains.
3.2. Phylogeography reconstruction of the H9N2 viruses in Iran and its neighboring countries
The basal branches of the HA tree were mainly occupied by Iran, with the root of the tree dating back to 1989. The virus later spread to Pakistan, most likely through a series of transmission events that began around 1992. The epidemic in Pakistan seems to be widely spread. The few available sequences from Iraq are dispersed independently throughout the tree, probably suggesting multiple independent clusters of infection between 1997 and 2008 and in 2013. Many sequences may not be available during this period [Figure 3].
|Figure 3: Maximum clade credibility phylogenies for the hemagglutinin gene of avian influenza A H9N2. The branches are colored according to the location of their nodes. The scale bar at the bottom indicates the years before the most recent sampling time (2019).|
Click here to view
A new cluster was also discovered in Pakistan, which diverged from Iranian lineages or lineages from intermediate locations in 1995, but whose data was not included in the current study. Between 2004-2008 and around 2012, we observed multiple transmission events from Pakistan to Iran [Figure 3]. The introduced lineages have experienced successful evolution and spread in Iran. The few sequences available from Afghanistan were nested within the Iranian lineages and dated back to 2014. Due to the small number of sequences available from Afghanistan and Iraq, it is not possible to provide a robust estimate of the date of origin as well as the magnitude of the outbreak in these countries [Figure 3].
Due to the high age and divergence of sequences in Iran and Pakistan in comparison with other countries in this sequence dataset, it seems that these two countries may have played a more highlighted role in the transmission pattern, but due to the insufficient sequences from other countries, conclusions about the role of the other countries may not be reliable.
Like the HA tree, the basal branches of the NA tree were also occupied by the sequences from Iran and Pakistan, respectively. The root height of the tree, which shows the estimated date of the epidemic’s onset in the region, was around 1992. The virus was first detected in Pakistan in 1996, and there appeared to be multiple back-transmission events from Pakistan to Iran in 2004, 2006, 2007, and 2010. The few available sequences from Afghanistan are nested within the Iranian branches, with a divergence date of 2008 and 2012 [Figure 4].
|Figure 4: Maximum clade credibility phylogenies for the neuraminidase gene of avian influenza A H9N2. The branches are colored according to the location of their nodes. The scale bar at the bottom indicates the years before the most recent sampling time (2019).|
Click here to view
The evolutionary rate for the HA and NA genes of the H9N2 viruses were estimated at 8.40 × 10-3 [95% HPD: (6.34-10.00) × 10-3] and 4.22 × 10-3 [95% HPD: (3.07-5.52) × 10-3] substitutions/site/year, respectively. The mean time of the most recent common ancestor of H9N2 viruses was 1988 (95% HPD: 1980-1995) for HA, and 1992 (95% HPD: 1986-1996) for NA.
| 4. Discussion|| |
To the best of our knowledge, this is the first Bayesian analysis of H9N2 viruses that attempts to fill the knowledge gap about the date of origin and circulation pathways of the influenza H9N2 virus in Iran and its neighboring countries. The results showed that the H9N2 epidemic may have started in Iran and Pakistan much earlier than the other investigated countries in the region, including Pakistan, Afghanistan, Armenia, Turkmenistan, Azerbaijan, Turkey, and Iraq. The onset of the H9N2 epidemic in the region was estimated to be between 1989 and 1992. An ongoing bidirectional dispersion of the virus between the investigated countries was also observed.
Our results showed that strains from the investigated countries were rooted in Pakistan and Iran. One justification for this may be the role of Pakistan and Iran in initiating the epidemic in the region. However, the abundance of sequence data available from Pakistan and Iran, and scarcity of data from other investigated countries, may have biased the topology of the phylogenetic tree. Lemey et al. discussed that sampling bias presents a critical challenge for popular discrete trait ancestral reconstruction procedures. As further sequence data become available from currently under-represented countries, more rigorous conclusions could be made.
Molecular clock analysis estimates the onset date of the epidemic in Pakistan and Iran was around 1992 and 1989, respectively. This finding suggests that the case detection and reporting along with the influenza surveillance program in the initial years had some failures and the laboratory detection of these serotypes occurred several years after their introduction into these countries.
There was evidence of the virus moving in both directions, or bidirectionally, between Iran and Pakistan. A cluster in Pakistan showed divergence from Iranian lineages or intermediate locations whose data was not included in this analysis. The large genetic distance between Iran and Pakistan suggests that this transmission is probably an indirect transmission that caused an outbreak in Pakistan in 1995. Between 2004-2008 and around 2012, we observed multiple introduction events from Pakistan to Iran. The introduced lineages have experienced successful evolution and spread in Iran.
As small number of sequences were available from Afghanistan and Iraq, so the epidemic onset date could not be robustly estimated in these countries. The few available sequences from Iraq are dispersed independently throughout the tree, suggesting multiple independent introductions of the virus and multiple epidemic clusters probably occurred between 1997 and 2008 and in 2013. The few sequences available from Afghanistan were also nested within the Iranian lineages and dated back to 2014.
This circulation may be associated with the fact that Iran is one of the destinations for religious and recreational tourists and refugees from neighboring countries, especially from Iraq and Afghanistan, due to the low socioeconomic level and the civil war. This could have also contributed to between-country virus transmission where neighboring countries’ tourists and refugees might have exchanged the virus with the Iranians.
A phylogenetic study showed that H9N2 viruses isolated from birds in the UAE, Iran, and Pakistan were closely related to the H9N2 viruses that infected humans in Hong Kong. In a broad phylogenetic study, all Iranian strains were located in one large cluster that also contained the sequences from neighboring countries, suggesting the epidemic in the region had originated from a common source. A study by Toroghi et al. found a high homology between Iranian and German H9N2 viruses and suggested that the epidemic might have originated in Germany or in intermediate locations in Europe. Detection of virus isolation from Pakistani parrots that were exported to Japan in 1998 was before the detection of the virus in Iran and with all these considerations it could be inferred that maybe Pakistan being infected before Iran,. The importance of ornamental birds in the spread of influenza viruses was then given more attention. Furthermore, given the high similarity of the NA and HA genes of the first Iranian influenza viruses isolated with the parakeet/chibal/1/97 virus, as well as the timing of virus isolation in Pakistan (1997, prior to the outbreak of Iranian avian influenza), the German origin of the epidemic in Iran may not be the case. Instead, the source of the epidemic in Iran may probably be the wildlife of its eastern countries, such as Pakistan. The wide common land border between Iran and Pakistan and the intensive trade and transportation of poultry products, as well as the trade of ornamental birds, may support the latter hypothesis.
Furthermore, some of the region’s people’s traditions must be considered for their effects on virus transmission dynamics. Some people living near the common border between Iran, Afghanistan, and Pakistan engage in the live bird trade as a source of food. It was recently discovered that live bird markets with a lack of biosecurity measures may be linked to the accelerated spread of the H7N9 virus in China.
Therefore, regular surveillance of H9N2 viruses in the country, especially in the live bird markets, enhancing the biosecurity of the poultry industry and screening newly arriving immigrants and tourists from neighboring countries at the border should be considered to control the spread of the virus, which is a notable threat to public health and subsequent human infections have been reported. Also, it is important that countries have a regional program in addition to domestic programs, and it will not be possible to control this disease without such a program. Furthermore, given the growing potential threat of the H9N2 virus to public health, surveillance of viral molecular evolution should be initiated. This is also critical for developing effective epidemic and pandemic prevention and control measures.
This study has some limitations. We could not deeply investigate the epidemiological parameters of the H9N2 virus in the neighboring countries like Iraq, Turkey, Armenia, Afghanistan and the Persian Gulf countries, due to the limited available sequence data. In addition, we could not directly estimate the evolutionary rates from the sequence alignments as the timespan of these datasets was not wide enough. Our model and dating estimates should be investigated further in the future as new sequences from additional countries and timespans become available, allowing for a more comprehensive evaluation of the H9N2 virus’s evolutionary history.
In summary, our results suggest that the H9N2 epidemic may have started in Iran and Pakistan much earlier than the other investigated countries in the region and was estimated to be between 1989 and 1992. Therefore, regular surveillance of H9N2 viruses in the country, especially in the live bird markets, enhancing the biosecurity of the poultry industry and screening newly arriving immigrants and tourists from neighboring countries at border, should be considered to control the spread of the virus. Also, it is important that countries have a regional program in addition to domestic programs and also surveillance of viral molecular evolution should be initiated. In the future, our model and dating estimates should be further investigated as new sequences from additional countries and timespans become available.
Conflict of interest statement
We declare that we have no conflict of interest.
We would like to express our gratitude to Institute for Futures Studies in Health, Kerman University of Medical Sciences, Kerman, Iran, especially Reza Sheikhzadeh (IT Manager) for his kind collaboration to access server computer in order to conduct computational analysis.
N.G, S.B and S.E had made substantial contributions to design and perform the research. N.G and S.E had been involved in acquisition of data, or analysis and interpretation of data. NG and S.E had been involved in drafting the manuscript. N.G, S.B, S.E and H.A had been revised initial draft critically for important intellectual content and given final approval of the version to be published.
| References|| |
Talat S, Abouelmaatti RR, Almeer R, Abdel-Daim MM, Elfeil WK. Comparison of the effectiveness of two different vaccination regimes for avian influenza H9N2 in broiler chicken. Animals
Shirvani E, Varghese BP, Paldurai A, Samal SK. A recombinant avian paramyxovirus serotype 3 expressing the hemagglutinin protein protects chickens against H5N1 highly pathogenic avian influenza virus challenge. Sci Rep
Philippon DA, Wu P, Cowling BJ, Lau EH. Avian influenza human infections at the human-animal interface. J Infect Dis
Khan M, Chaudhry M, Fatima Z, Khan RU, Ahmad B, Ullah R, et al. Effect of avian influenza H9N2 subtype virus infection on backyard poultry production. Sci Lett
Lee YJ, Shin JY, Song MS, Lee YM, Choi JG, Lee EK, et al. Continuing evolution of H9 influenza viruses in Korean poultry. Virol
Vasfi MM, Bozorg MM, editors. An outbreak of non-highly pathogenic avian influenza in chickens in Iran 1999
. Proceedings meeting of World Veterinary Association. Cairo, 2002.
Fallah Mehrabadi MH, Shoushtari A, Tehrani F, Motamed N, Haerian B, Ghalyanchilangeroudi A, et al. Serological and molecular survey of avian influenza H9N2 subtype in live birds markets-2016. J Vet Res
Cabello M, Mendoza Y, Bello G. Spatiotemporal dynamics of dissemination of non-pandemic HIV-1 subtype B clades in the Caribbean region. PloS One
(8): e106045. doi: 10.1371/journal. pone.0106045.
Elgendy EM, Arai Y, Kawashita N, Isobe A, Daidoji T, Ibrahim MS, et al. Double mutations in the H9N2 avian influenza virus PB2
gene act cooperatively to increase viral host adaptation and replication for human infections. J General Virol
(6): 001612. doi: 10.1099/ jgv.0.001612.
Paraskevis D, Pybus O, Magiorkinis G, Hatzakis A, Wensing AM, van de Vijver DA, et al. Tracing the HIV-1 subtype B mobility in Europe: A phylogeographic approach. Retrovirol
Song Xf, Han P, Chen YPP. Genetic variation of the hemagglutinin of avian influenza virus H9N2. J Med Virol
Fusaro A, Monne I, Salviato A, Valastro V, Schivo A, Amarin NM, et al. Phylogeography and evolutionary history of reassortant H9N2 viruses with potential human health implications. J Virol
Lee DH, Swayne DE, Sharma P, Rehmani SF, Wajid A, Suarez DL, et al. H9N2 low pathogenic avian influenza in Pakistan (2012–2015). Vet Rec Open
(1): e000171. doi: 10.1136/vetreco-2016-000171.
Rich SN, Richards VL, Mavian CN, Switzer WM, Rife Magalis B, Poschman K, et al. Employing molecular phylodynamic methods to identify and forecast HIV transmission clusters in public health settings: A qualitative study. Viruses
Yang J, Müller NF, Bouckaert R, Xu B, Drummond AJ. Bayesian phylodynamics of avian influenza A virus H9N2 in Asia with time-dependent predictors of migration. PLoS Comput Biol
(8): e1007189. doi: https://doi.org/10.1371/journal.pcbi.1007189
Kumar S, Stecher G, Tamura K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol
Drummond A, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol
Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol
Jin Y, Yu D, Ren H, Yin Z, Huang Z, Hu M, et al. Phylogeography of Avian influenza A H9N2 in China. BMC Genom
Suchard A, Lemey P, Baele G, Ayres D, Drummond A, Rambaut A. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol
(1): vey016. doi: 10.1093/ve/vey016.
Bielejec F, Rambaut A, Suchard MA, Lemey P. SPREAD: Spatial phylogenetic reconstruction of evolutionary dynamics. Bioinform
Lemey P, Hong SL, Hill V, Baele G, Poletto C, Colizza V, et al. Accommodating individual travel history and unsampled diversity in Bayesian phylogeographic inference of SARS-CoV-2. Nat Commun
Aamir U, Wernery U, Ilyushina N, Webster R. Characterization of avian H9N2 influenza viruses from United Arab Emirates 2000 to 2003. Virology
Toroghi R, Momayez R. Biological and molecular characterization of Avian influenza virus (H9N2) isolates from Iran. Acta Virol
Karimi V, Bozorgmehri Fard MH, Shahbazzadeh D, Esmaelizad M, Pourbakhsh SA. Sequence analysis and phylogenetic study of hemagglutinin gene of H9N2 subtype of avian influenza virus isolated during 1998-2002 in Iran. Iran Biomed J
Naeem K. An outbreak of avian influenza in poultry in Pakistan. Vet Rec
Soltanialvar M, Shoushtari H, Bozorgmehrifard M, Charkhkar S, Eshratabadi F. Molecular characterization of hemagglutinin and neuraminidase genes of H9N2 avian influenza viruses isolated from commercial broiler chicken in Iran. J Biol Sci
Bao CJ, Cui LB, Zhou MH, Hong L, Gao GF, Wang H. Live-animal markets and influenza A (H7N9) virus infection. N Engl J Med
Naguib MM, Arafa ASA, El-Kady MF, Selim AA, Gunalan V, Maurer-Stroh S, et al. Evolutionary trajectories and diagnostic challenges of potentially zoonotic avian influenza viruses H5N1 and H9N2 co-circulating in Egypt. Infect Gen Evol
[Figure 1], [Figure 2], [Figure 3], [Figure 4]