Background: Ongoing outbreaks of H5N1 highly pathogenic avian influenza (HPAI) viruses and the emergence of the genetic-related hemagglutinin (HA) gene of reassortant H5Nx viruses currently circulating in wild birds and poultries pose a great global public health concern. In this study, we comprehensively analyzed the genetic evolution of Thai H5N1 HA and neuraminidase (NA) genes between 2003 and 2010. The H5N1 Thailand virus clade 2.3.4 was also genetically compared to the currently circulating clade 2.3.4.4 of H5Nx viruses.

Methods: Full-length nucleotide sequences of 178 HA and 143 NA genes of H5N1 viruses circulating between 2003 and 2010 were phylogenetically analyzed using maximum likelihood (ML) phylogenetic construction. Bayesian phylogenetic trees were reconstructed using BEAST analysis with a Bayesian Markov chain Monte Carlo (MCMC) approach. The maximum clade credibility (MCC) tree was determined, and the time of the most recent common ancestor (tMRCA) was estimated. The H5N1 HA nucleotide sequences of clade 2.3.4 Thailand viruses were phylogenetically analyzed using ML phylogenetic tree construction and analyzed for nucleotide similarities with various subtypes of reassortant H5Nx HA clade 2.3.4.4.

Results: ML phylogenetic analysis revealed two distinct HA clades, clade 1 and clade 2.3.4, and two distinct NA groups within the corresponding H5 clade 1 viruses. Bayesian phylogenetic reconstruction for molecular clock suggested that the Thai H5N1 HA and NA emerged in 2001.87 (95% HPD: 2001.34-2002.49) and 2002.38 (95% HPD: 2001.99-2002.82), respectively, suggesting that the virus existed before it was first reported in 2004. The Thai H5N1 HA clade 2.3.4 was grouped into corresponding clades 2.3.4, 2.3.4.1, 2.3.4.2, and 2.3.4.3, and shared nucleotide similarities to reassortant H5Nx clade 2.3.4.4 ranged from 92.4-96.8%. Phylogenetic analysis revealed monophyletic H5Nx clade 2.3.4.4 evolved from H5N1 clade 2.3.4.

Conclusion: H5N1 viruses existed, and were presumably introduced and circulated in avian species in Thailand, before they were officially reported in 2004. HA and NA genes continuously evolved during circulation between 2004 and 2010. This study provides a better understanding of genetic evolution with respect to molecular epidemiology. Monitoring and surveillance of emerging variants/reassortants should be continued.