Background: Influenza A viruses pose a significant public health threat globally and are characterized by rapid evolution of the hemagglutinin (HA) gene causing seasonal epidemics. The aim of this study was to investigate the evolutionary dynamics of A(H3N2) circulating in Bhutan during 2022 and 2023.

Methods: We analysed 166 whole-genome sequences of influenza A(H3N2) from Bhutan, obtained from the GISAID database. We employed a Bayesian Markov Chain Monte Carlo (MCMC) framework, with a curated global dataset of HA sequences from regions with significant migration links to Bhutan. Phylogenetic, temporal, and phylogeographic analyses were conducted to elucidate the evolutionary dynamics and spatial dissemination of the viruses.

Results: Our phylogenetic analysis identified the circulation of influenza A(H3N2) Clade 3C.2a1b.2a.2 in Bhutan during 2022 and 2023, with viruses further classified into three subclades: 2a.3 (39/166), 2a.3a.1 (58/166) and 2a.3b (69/166). The TMRCA estimates suggest that these viral lineages originated approximately 1.93 years prior to their detection. Phylogeographic analysis indicates introductions from the United States in 2022 and Australia in 2023. The mean evolutionary rate across all gene segments was calculated to be 4.42 × 10-3 substitutions per site per year (95% HPD: 3.19 × 10-3 to 5.84 × 10-3), with evidence of purifying selection and limited genetic diversity. Furthermore, reassortment events were rare, with an estimated rate of 0.045 events per lineage per year.

Conclusion: Our findings show that primary forces shaping the local evolution of the influenza A(H3N2) in Bhutan are largely stochastic, with only sporadic instances of adaptive change, and thus underscore the importance of continuous surveillance to mitigate the impact of evolving strains.