Evolutionary dynamics of influenza A/H1N1 virus circulating in India from 2011 to 2021

Influenza A viruses (IAV) are fast-evolving pathogens with a very high mutation rate (2.0 × 10^-6 to 2.0 × 10^-4) compared to the influenza B (IBV) and influenza C (ICV) viruses. Generally, the tropical regions are considered as the reservoir for the IAV's genetic and antigenic evolutionary modification to be reintroduced into the temperate region. Therefore, in connection to the above facts, the present study emphasized on the evolutionary dynamic of the pandemic-2009 H1N1 (pdmH1N1) influenza virus in India. A total of Ninety-two whole genome sequences of pdmH1N1 viruses circulating in India during the 2009 post-pandemic era were analysed. The temporal signal of the study, indicating strict molecular clock evolutionary process and the overall substitution rate is 2.21 × 10^-3/site/year. We are using the nonparametric Bayesian Skygrid coalescent model to estimates the effective past population dynamic or size over time. The study exhibits a strong relation between the genetic distances and collection dates of the Indian pdmH1N1 strain. The skygrid plot represents the highest exponential growth of IAV in rainy and winter seasons. All the genes of Indian pdmH1N1 were under purifying selective pressure. The Bayesian time-imprinted phylogenetic tree represents the following clade distributions in the country within the last 10 years; I) clade 6, 6C, and 7 were co-circulating between the 2011 to 2012 flu season; II) the clade 6B was introduced into circulation in the late seasons of 2012; III) lastly, the clade 6B remain existing in the circulation and segregated into subclade 6B.1 with five different subgroup (6B.1A, 6B.1A.1, 6B.1A.5a, 6B.1A.5a.2, 6B.1A.7). The recent circulating strain of Indian H1N1 strain represent the insertion of basic-amino acid arginine (R) in the cleavage site (325/K-R) of HA protein and amino acid mutation (314/I-M) on the lateral head surface domain of NA protein. Moreover, the study indicates the sporadic presence of the oseltamivir-resistant (275/H-Y) H1N1 variant in circulation. The present study suggests the purifying selective pressure and stochastic ecological factors for the existence and adaptation of a certain clade 6B in the host populations and additional information on the emergence of mutated strains in the circulation.