Whole genome characterization, phylogenetic and genome signature analysis of human pandemic H1N1 virus in Thailand, 2009-2012

Evidence of avian and human influenza A virus infection in farmed Siamese crocodiles (Crocodylus siamensis) in Thailand

Crocodilians are susceptible to a range of virus infection including influenza A virus (IAV). However, little is known about the ecology and epidemiology of IAV in crocodile species. This study aimed to investigate IAV infection in farmed Siamese crocodiles in central Thailand. We collected plasma samples and pharyngeal swab samples from Siamese crocodiles residing in 13 crocodile farms in 9 provinces of central Thailand during 2019. Additional archival plasma samples of Siamese crocodiles collected in 2012 and 2018 were also included in the study. Plasma samples were screened for influenza A antibodies by a hemagglutination inhibition (HI) assay and positive were evaluated by a cytopathic effect/hemagglutination based-microneutralization (MN) assay. Swab samples were tested for influenza viral RNA by a real-time RT-PCR assay targeting the influenza matrix (M) gene. Among 246 tested plasma samples, the overall seroprevalence of antibodies against IAV in farmed Siamese crocodiles was 17.5% (43/246). The most common hemagglutinin (HA) subtype was H2 (46.5%, 20/43) followed by H9 (39.5%, 17/43), human H1 (14%, 6/43) and H1 (7%, 3/43). Multiple HA subtypes were also detected in 7% (3/43) of infected crocodiles with combination of H1 and H2 subtypes. All 126 tested swab samples were negative for influenza viral RNA. In addition, we demonstrated the ability of wild-type IAV subtypes (H1, H2, H9 and human H1) to infect primary Siamese crocodile fibroblast cells. To our knowledge, this is the first report of serological evidences of avian and human IAV infection in Siamese crocodiles. Our findings highlighted the role of crocodile species in the ecology of IAV particularly the potential to serve as the reservoir or mixing vessel for the viruses that significantly threaten both human and animal health.

Global Prevalence and Hemagglutinin Evolution of H7N9 Avian Influenza Viruses from 2013 to 2022

H7N9 avian influenza viruses have caused severe harm to the global aquaculture industry and human health. For further understanding of the characteristics of prevalence and hemagglutinin evolution of H7N9 avian influenza viruses, we generated the global epidemic map of H7N9 viruses from 2013 to 2022, constructed a phylogenetic tree, predicted the glycosylation sites and compared the selection pressure of the hemagglutinin. The results showed that although H7N9 avian influenza appeared sporadically in other regions worldwide, China had concentrated outbreaks from 2013 to 2017. The hemagglutinin genes were classified into six distinct lineages: A, B, C, D, E and F. After 2019, H7N9 viruses from the lineages B, E and F persisted, with the lineage B being the dominant. The hemagglutinin of highly pathogenic viruses in the B lineage has an additional predicted glycosylation site, which may account for their persistent pandemic, and is under more positive selection pressure. The most recent ancestor of the H7N9 avian influenza viruses originated in September 1991. The continuous evolution of hemagglutinin has led to an increase in virus pathogenicity in both poultry and humans, and sustained human-to-human transmission. This study provides a theoretical basis for better prediction and control of H7N9 avian influenza.

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.

Genomic Signature in Evolutionary Biology: A Review

Organisms are unique physical entities in which information is stored and continuously processed. The digital nature of DNA sequences enables the construction of a dynamic information reservoir. However, the distinction between the hardware and software components in the information flow is crucial to identify the mechanisms generating specific genomic signatures. In this work, we perform a bibliometric analysis to identify the different purposes of looking for particular patterns in DNA sequences associated with a given phenotype. This study has enabled us to make a conceptual breakdown of the genomic signature and differentiate the leading applications. On the one hand, it refers to gene expression profiling associated with a biological function, which may be shared across taxa. This signature is the focus of study in precision medicine. On the other hand, it also refers to characteristic patterns in species-specific DNA sequences. This interpretation plays a key role in comparative genomics, identifying evolutionary relationships. Looking at the relevant studies in our bibliographic database, we highlight the main factors causing heterogeneities in genome composition and how they can be quantified. All these findings lead us to reformulate some questions relevant to evolutionary biology.

Genetic and Antigenic Evolution of European Swine Influenza A Viruses of HA-1C (Avian-Like) and HA-1B (Human-Like) Lineages in France from 2000 to 2018

This study evaluated the genetic and antigenic evolution of swine influenza A viruses (swIAV) of the two main enzootic H1 lineages, i.e., HA-1C (H1_av) and -1B (H1_hu), circulating in France between 2000 and 2018. SwIAV RNAs extracted from 1220 swine nasal swabs were hemagglutinin/neuraminidase (HA/NA) subtyped by RT-qPCRs, and 293 virus isolates were sequenced. In addition, 146 H1_avNy and 105 H1_huNy strains were submitted to hemagglutination inhibition tests. H1_avN1 (66.5%) and H1_huN2 (25.4%) subtypes were predominant. Most H1 strains belonged to HA-1C.2.1 or -1B.1.2.3 clades, but HA-1C.2, -1C.2.2, -1C.2.3, -1B.1.1, and -1B.1.2.1 clades were also detected sporadically. Within HA-1B.1.2.3 clade, a group of strains named "Δ146-147" harbored several amino acid mutations and a double deletion in HA, that led to a marked antigenic drift. Phylogenetic analyses revealed that internal segments belonged mainly to the "Eurasian avian-like lineage", with two distinct genogroups for the M segment. In total, 17 distinct genotypes were identified within the study period. Reassortments of H1_av/H1_hu strains with H1N1pdm virus were rarely evidenced until 2018. Analysis of amino acid sequences predicted a variability in length of PB1-F2 and PA-X proteins and identified the appearance of several mutations in PB1, PB1-F2, PA, NP and NS1 proteins that could be linked to virulence, while markers for antiviral resistance were identified in N1 and N2. Altogether, diversity and evolution of swIAV recall the importance of disrupting the spreading of swIAV within and between pig herds, as well as IAV inter-species transmissions.

Evolution of the neuraminidase gene of seasonal influenza A and B viruses in Thailand between 2010 and 2015

The neuraminidase inhibitors (NAIs) oseltamivir and zanamivir are commonly used for the treatment and control of influenza A and B virus infection. However, the emergence of new influenza virus strains with reduced susceptibility to NAIs may appear with the use of these antivirals or even naturally. We therefore screened the neuraminidase (NA) sequences of seasonal influenza virus A(H1N1), A(H1N1)pdm09, A(H3N2), and influenza B virus strains identified in Thailand for the presence of substitutions previously reported to reduce susceptibility to NAIs. We initially examined oseltamivir resistance (characterized by the H275Y mutation in the NA gene) in 485 A(H1N1)pdm09 strains circulating in Thailand and found that 0.82% (4/485) had this substitution. To further evaluate the evolution of the NA gene, we also randomly selected 98 A(H1N1)pdm09, 158 A(H3N2), and 69 influenza B virus strains for NA gene amplification and sequencing, which revealed various amino acid mutations in the active site of the NA protein previously shown to be associated with reduced susceptibility to NAIs. Phylogenetic analysis of the influenza virus strains from this study and elsewhere around the world, together with the estimations of nucleotide substitution rates and selection pressure, and the predictions of B-cell epitopes and N-linked glycosylation sites all provided evidence for the ongoing evolution of NA. The overall rates of NA evolution for influenza A viruses were higher than for influenza B virus at the nucleotide level, although influenza B virus possessed more genealogical diversity than that of influenza A viruses. The continual surveillance of the antigenic changes associated with the NA protein will not only contribute to the influenza virus database but may also provide a better understanding of selection pressure exerted by antiviral use.

Age-specific genetic and antigenic variations of influenza A viruses in Hong Kong, 2013-2014

Age-specific genetic and antigenic variations of influenza viruses have not been documented in tropical and subtropical regions. We implemented a systematic surveillance program in two tertiary hospitals in Hong Kong Island, to collect 112 A(H1N1)pdm09 and 254 A(H3N2) positive specimens from 2013 to 2014. Of these, 56 and 72 were identified as genetic variants of the WHO recommended vaccine composition strains, respectively. A subset of these genetic variants was selected for hemagglutination-inhibition (HI) tests, but none appeared to be antigenic variants of the vaccine composition strains. We also found that genetic and antigenicity variations were similar across sex and age groups of ≤18 yrs, 18 to 65 yrs, and ≥65 yrs. Our findings suggest that none of the age groups led other age groups in genetic evolution of influenza virus A strains. Future studies from different regions and longer study periods are needed to further investigate the age and sex heterogeneity of influenza viruses.

Uncovering the Potential Pan Proteomes Encoded by Genomic Strand RNAs of Influenza A Viruses

Influenza A virus genomes are composed of eight negative sense RNAs. In total, 16 proteins encoded by eight positive sense RNAs were identified. One putative protein coding sequence (PCS) encoded by genomic strand RNA of segment 8 has been previously proposed. In this study, 95,608, 123,965 and 35,699 genomic strand RNA sequences from influenza A viruses from avian, human and mammalian hosts, respectively, were used to identify PCSs encoded by the genomic strand RNAs. In total, 326,069 PCSs with lengths equal to or longer than 80 amino acids were identified and clustered into 270 PCS groups. Twenty of the 270 PCS groups which have greater than 10% proportion in influenza A viruses from avian, human or mammalian hosts were selected for detailed study. Maps of the 20 PCSGs in the influenza A virus genomes were constructed. The proportions of the 20 PCSGs in influenza A viruses from different hosts and serotypes were analyzed. One secretory and five membrane proteins predicted from the PCS groups encoded by genomic strand RNAs of segments 1, 2, 4, 6, 7 and 8 were identified. These results suggest the possibility of the ambisense nature of the influenza A virus genomic RNAs and a potential coding sequence reservoir encoding potential pan proteomes of influenza A viruses.

New genetic variants of influenza A(H1N1)pdm09 detected in Cuba during 2011-2013

Influenza A(H1N1)pdm09 virus has evolved continually since its emergence in 2009. For influenza virus strains, genetic changes occurring in HA1 domain of the hemagglutinin cause the emergence of new variants. The aim of our study is to establish genetic associations between 35 A(H1N1)pdm09 viruses circulating in Cuba in 2011-2012 and 2012-2013 seasons, and A/California/07/2009 strain recommended by WHO as the H1N1 component of the influenza vaccine. The phylogenetic analysis revealed the circulation of clades 3, 6A, 6B, 6C and 7. Mutations were detected in the antigenic site or in the receptor-binding domains of HA1 segment, including S174P, S179N, K180Q, S202T, S220T and R222K. Substitutions S174P, S179N, K180Q and R222K were detected in Cuban strains for the first time.

Diversifying Selection Analysis Predicts Antigenic Evolution of 2009 Pandemic H1N1 Influenza A Virus in Humans

Although a large number of immune epitopes have been identified in the influenza A virus (IAV) hemagglutinin (HA) protein using various experimental systems, it is unclear which are involved in protective immunity to natural infection in humans. We developed a data mining approach analyzing natural H1N1 human isolates to identify HA protein regions that may be targeted by the human immune system and can predict the evolution of IAV. We identified 16 amino acid sites experiencing diversifying selection during the evolution of prepandemic seasonal H1N1 strains and found that 11 sites were located in experimentally determined B-cell/antibody (Ab) epitopes, including three distinct neutralizing Caton epitopes: Sa, Sb, and Ca2 [A. J. Caton, G. G. Brownlee, J. W. Yewdell, and W. Gerhard, Cell 31:417-427, 1982, http://dx.doi.org/10.1016/0092-8674(82)90135-0]. We predicted that these diversified epitope regions would be the targets of mutation as the 2009 H1N1 pandemic (pH1N1) lineage evolves in response to the development of population-level protective immunity in humans. Using a chi-squared goodness-of-fit test, we identified 10 amino acid sites that significantly differed between the pH1N1 isolates and isolates from the recent 2012-2013 and 2013-2014 influenza seasons. Three of these sites were located in the same diversified B-cell/Ab epitope regions as identified in the analysis of prepandemic sequences, including Sa and Sb. As predicted, hemagglutination inhibition (HI) assays using human sera from subjects vaccinated with the initial pH1N1 isolate demonstrated reduced reactivity against 2013-2014 isolates. Taken together, these results suggest that diversifying selection analysis can identify key immune epitopes responsible for protective immunity to influenza virus in humans and thereby predict virus evolution.

IMPORTANCE

The WHO estimates that approximately 5 to 10% of adults and 20 to 30% of children in the world are infected by influenza virus each year. While an adaptive immune response helps eliminate the virus following acute infection, the virus rapidly evolves to evade the established protective memory immune response, thus allowing for the regular seasonal cycles of influenza virus infection. The analytical approach described here, which combines an analysis of diversifying selection with an integration of immune epitope data, has allowed us to identify antigenic regions that contribute to protective immunity and are therefore the key targets of immune evasion by the virus. This information can be used to determine when sequence variations in seasonal influenza virus strains have affected regions responsible for protective immunity in order to decide when new vaccine formulations are warranted.

Amino acid substitutions in hemagglutinin of the 2009 pandemic influenza A(H1N1) viruses that might affect the viral antigenicity

BACKGROUND

During 2009 to 2012, Thailand had encountered 4 distinctive waves of the 2009 pandemic influenza A(H1N1) (H1N1pdm) outbreaks. Considering the RNA nature of the influenza viral genome, a mutation in hemagglutinin (HA) gene which led to change in antigenicity of the strains circulating during those epidemic periods is anticipated. It is also uncertain whether the A/California/07/2009 (H1N1) (CA/07) vaccine strain still confers protective immunity against those evolved viruses, the causative agents of the later epidemic waves.

METHODS

HA gene segments of 10 H1N1pdm isolates obtained during 2009 to 2012 were sequenced and phylogenetically analysed using ClustalW and MEGA5 programs. A total of 124 convalescent serum samples collected from patients naturally infected during 3 epidemic waves were employed as tools to investigate for antigenic change in HA of these 10 circulating H1N1pdm viruses by hemagglutination inhibition (HI) assay.

RESULTS

A phylogenetic analysis showed that the 10 virus isolates were grouped into 4 clusters corresponding to the time of 4 consecutive outbreaks. An accumulation of amino acid substitutions in HA was observed in viruses derived from the late epidemic waves. Significantly lower antibody titers were observed when CA/07 was tested against convalescent sera collected from the 3 waves (p<0.05) compared to most of Thai isolates; and significantly lower antibody titers were also obtained when virus isolates, retrieved from the third epidemic wave were tested against convalescent sera collected during the first and second wave. These results were suggestive of change in antigenicity of the evolved viruses. Our results also showed some mutation position residing outside the previously reported antigenic site that may involve in an alteration of the viral antigenicity.

CONCLUSIONS

Our study demonstrated that convalescent sera collected from individuals naturally infected with H1N1pdm virus were successfully used to reveal a statistically significant change in antibody titers against the currently evolved H1N1pdm viruses as determined by HI assay. Nevertheless, the antibody titers of individual serum against various viruses were less than 4-folded difference as compared to that against the CA/07 vaccine strain. Therefore, CA/07 is still a potent vaccine strain for those evolved H1N1pdm viruses.

Molecular characterization of avian influenza H5N1 virus in Egypt and the emergence of a novel endemic subclade

Clade 2.2 highly pathogenic H5N1 viruses have been in continuous circulation in Egyptian poultry since 2006. Their persistence caused significant genetic drift that led to the reclassification of these viruses into subclades 2.2.1 and 2.2.1.1. Here, we conducted full-genome sequence and phylogenetic analyses of 45 H5N1 isolated during 2006-2013 through systematic surveillance in Egypt, and 53 viruses that were sequenced previously and available in the public domain. Results indicated that H5N1 viruses in Egypt continue to evolve and a new distinct cluster has emerged. Mutations affecting viral virulence, pathogenicity, transmission, receptor-binding preference and drug resistance were studied. In light of our findings that H5N1 in Egypt continues to evolve, surveillance and molecular studies need to be sustained.

Influenza A H1N1 pandemic strain evolution--divergence and the potential for antigenic drift variants

The emergence of a novel A(H1N1) strain in 2009 was the first influenza pandemic of the genomic age, and unprecedented surveillance of the virus provides the opportunity to better understand the evolution of influenza. We examined changes in the nucleotide coding regions and the amino acid sequences of the hemagglutinin (HA), neuraminidase (NA), and nucleoprotein (NP) segments of the A(H1N1)pdm09 strain using publicly available data. We calculated the nucleotide and amino acid hamming distance from the vaccine strain A/California/07/2009 for each sequence. We also estimated Pepitope-a measure of antigenic diversity based on changes in the epitope regions-for each isolate. Finally, we compared our results to A(H3N2) strains collected over the same period. Our analysis found that the mean hamming distance for the HA protein of the A(H1N1)pdm09 strain increased from 3.6 (standard deviation [SD]: 1.3) in 2009 to 11.7 (SD: 1.0) in 2013, while the mean hamming distance in the coding region increased from 7.4 (SD: 2.2) in 2009 to 28.3 (SD: 2.1) in 2013. These trends are broadly similar to the rate of mutation in H3N2 over the same time period. However, in contrast to H3N2 strains, the rate of mutation accumulation has slowed in recent years. Our results are notable because, over the course of the study, mutation rates in H3N2 similar to that seen with A(H1N1)pdm09 led to the emergence of two antigenic drift variants. However, while there has been an H1N1 epidemic in North America this season, evidence to date indicates the vaccine is still effective, suggesting the epidemic is not due to the emergence of an antigenic drift variant. Our results suggest that more research is needed to understand how viral mutations are related to vaccine effectiveness so that future vaccine choices and development can be more predictive.

Pathogenicity of pandemic H1N1 influenza A virus in immunocompromised cynomolgus macaques

Pandemic (H1N1) 2009 influenza virus spread throughout the world since most people did not have immunity against the virus. In the post pandemic phase when many humans might possess immunity against the pandemic virus, one of the concerns is infection in immunocompromised people. Therefore, we used an immunosuppressed macaque model to examine pathogenicity of the pandemic (H1N1) 2009 virus under an immunocompromised condition. The virus in nasal samples of immunosuppressed macaques infected with the pandemic (H1N1) 2009 virus was detected longer after infection than was the virus in nasal samples of immunocompetent macaques. As expected, not only virus amounts but also virus propagation sites in the immunosuppressed macaques were larger than those in lungs of the immunocompetent macaques when they were infected with the pandemic virus. Immunosuppressed macaques possessed low levels of immune cells producing cytokines and chemokines, but levels of inflammatory cytokines/chemokine interleukin (IL)-6, IL-18, and monocyte chemotactic protein (MCP)-1 in lungs of the immunosuppressed macaques were higher than those in lungs of the immunocompetent macaques, though the differences were not statistically significant. Therefore, under an immunosuppressive condition, the pandemic influenza (H1N1) 2009 virus might cause more severe morbidity with high cytokine/chemokine production by the host innate immune system than that seen in macaques under the immunocompetent condition.

Evaluation of neutralizing efficacy of monoclonal antibodies specific for 2009 pandemic H1N1 influenza A virus in vitro and in vivo

Pandemic influenza A virus (H1N1) 2009 poses a serious public-health challenge worldwide. To characterize the neutralizing epitopes of this virus, we generated a panel of eight monoclonal antibodies (mAbs) against the HA of the A/California/07/2009 virus. The antibodies were specific for the 2009 pdm H1N1 HA, as the antibodies displayed HA-specific ELISA, hemagglutination inhibition (HAI) and neutralization activity. One mAb (mAb12) showed significantly higher HAI and neutralizing titers than the other mAbs. We mapped the antigenic epitopes of the HA by characterizing escape mutants of a 2009 H1N1 vaccine strain (NYMC X-179A). The amino acid changes suggested that these eight mAbs recognized HA antigenic epitopes located in the Sa, Sb, Ca1 and Ca2 sites. Passive immunization with mAbs showed that mAb12 displayed more efficient neutralizing activity in vivo than the other mAbs. mAb12 was also found to be protective, both prophylactically and therapeutically, against a lethal viral challenge in mice. In addition, a single injection of 10 mg/kg mAb12 outperformed a 5-day course of treatment with oseltamivir (10 mg/kg/day by gavage) with respect to both prophylaxis and treatment of lethal viral infection. Taken together, our results showed that mouse-origin mAbs displayed neutralizing effectiveness in vitro and in vivo. One mAb in particular (mAb12) recognized an epitope within the Sb site and demonstrated outstanding neutralizing effectiveness.