Whole Genome Characterization,Phylogenetic and Genome Signature Analysis of Human Pandemic H1N1 Virus in Thailand, 2009–2012

Background

Three waves of human pandemic influenza occurred in Thailand in 2009–2012. The genome signature features and evolution of pH1N1 need to be characterized to elucidate the aspects responsible for the multiple waves of pandemic.

Methodology/Findings

Forty whole genome sequences and 584 partial sequences of pH1N1 circulating in Thailand, divided into 1^st, 2^nd and 3^rd wave and post-pandemic were characterized and 77 genome signatures were analyzed. Phylogenetic trees of concatenated whole genome and HA gene sequences were constructed calculating substitution rate and d_N/d_S of each gene. Phylogenetic analysis showed a distinct pattern of pH1N1 circulation in Thailand, with the first two isolates from May, 2009 belonging to clade 5 while clades 5, 6 and 7 co-circulated during the first wave of pH1N1 pandemic in Thailand. Clade 8 predominated during the second wave and different proportions of the pH1N1 viruses circulating during the third wave and post pandemic period belonged to clades 8, 11.1 and 11.2. The mutation analysis of pH1N1 revealed many adaptive mutations which have become the signature of each clade and may be responsible for the multiple pandemic waves in Thailand, especially with regard to clades 11.1 and 11.2 as evidenced with V731I, G154D of PB1 gene, PA I330V, HA A214T S160G and S202T. The substitution rate of pH1N1 in Thailand ranged from 2.53×10⁻³±0.02 (M2 genes) to 5.27×10⁻³±0.03 per site per year (NA gene).

Conclusions

All results suggested that this virus is still adaptive, maybe to evade the host''s immune response and tends to remain in the human host although the d_N/d_S were under purifying selection in all 8 genes. Due to the gradual evolution of pH1N1 in Thailand, continuous monitoring is essential for evaluation and surveillance to be prepared for and able to control future influenza activities.     

Assessing Antigenic Drift of Seasonal Influenza A(H3N2) and A(H1N1)pdm09 Viruses

Under selective pressure from the host immune system, antigenic epitopes of influenza virus hemagglutinin (HA) have continually evolved to escape antibody recognition, termed antigenic drift. We analyzed the genomes of influenza A(H3N2) and A(H1N1)pdm09 virus strains circulating in Thailand between 2010 and 2014 and assessed how well the yearly vaccine strains recommended for the southern hemisphere matched them. We amplified and sequenced the HA gene of 120 A(H3N2) and 81 A(H1N1)pdm09 influenza virus samples obtained from respiratory specimens and calculated the perfect-match vaccine efficacy using the p _epitope model, which quantitated the antigenic drift in the dominant epitope of HA. Phylogenetic analysis of the A(H3N2) HA1 genes classified most strains into genetic clades 1, 3A, 3B, and 3C. The A(H3N2) strains from the 2013 and 2014 seasons showed very low to moderate vaccine efficacy and demonstrated antigenic drift from epitopes C and A to epitope B. Meanwhile, most A(H1N1)pdm09 strains from the 2012–2014 seasons belonged to genetic clades 6A, 6B, and 6C and displayed the dominant epitope mutations at epitopes B and E. Finally, the vaccine efficacy for A(H1N1)pdm09 (79.6–93.4%) was generally higher than that of A(H3N2). These findings further confirmed the accelerating antigenic drift of the circulating influenza A(H3N2) in recent years.     

Molecular Epidemiology and Phylogenetic Analyses of Influenza B Virus in Thailand during 2010 to 2014

Influenza B virus remains a major contributor to the seasonal influenza outbreak and its prevalence has increased worldwide. We investigated the epidemiology and analyzed the full genome sequences of influenza B virus strains in Thailand between 2010 and 2014. Samples from the upper respiratory tract were collected from patients diagnosed with influenza like-illness. All samples were screened for influenza A/B viruses by one-step multiplex real-time RT-PCR. The whole genome of 53 influenza B isolates were amplified, sequenced, and analyzed. From 14,418 respiratory samples collected during 2010 to 2014, a total of 3,050 tested positive for influenza virus. Approximately 3.27% (471/14,418) were influenza B virus samples. Fifty three isolates of influenza B virus were randomly chosen for detailed whole genome analysis. Phylogenetic analysis of the HA gene showed clusters in Victoria clades 1A, 1B, 3, 5 and Yamagata clades 2 and 3. Both B/Victoria and B/Yamagata lineages were found to co-circulate during this time. The NA sequences of all isolates belonged to lineage II and consisted of viruses from both HA Victoria and Yamagata lineages, reflecting possible reassortment of the HA and NA genes. No significant changes were seen in the NA protein. The phylogenetic trees generated through the analysis of the PB1 and PB2 genes closely resembled that of the HA gene, while trees generated from the analysis of the PA, NP, and M genes showed similar topology. The NS gene exhibited the pattern of genetic reassortment distinct from those of the PA, NP or M genes. Thus, antigenic drift and genetic reassortment among the influenza B virus strains were observed in the isolates examined. Our findings indicate that the co-circulation of two distinct lineages of influenza B viruses and the limitation of cross-protection of the current vaccine formulation provide support for quadrivalent influenza vaccine in this region.