Positive selection on hemagglutinin and neuraminidase genes of H1N1 influenza viruses

Enhanced Diversifying Selection on Polymerase Genes in H5N1 Clade 2.3.4.4b: A Key Driver of Altered Species Tropism and Host Range Expansion

Highly pathogenic avian influenza H5N1 clade 2.3.4.4b viruses have shown unprecedented host range and pathogenicity, including infections in cattle, previously not susceptible to H5N1. We investigated whether selection pressures on clade 2.3.4.4b viral genes could shed light on their unique epidemiological features. Our analysis revealed that while the gene products of clade 2.3.4.4b H5N1 primarily undergo purifying selection, there are notable instances of episodic diversifying selection. Specifically, the polymerase genes PB2, PB1, and PA exhibit significantly greater selection pressures in clade 2.3.4.4b than all earlier H5N1 virus clades. Polymerases play critical roles in influenza virus adaptation, including viral fitness, interspecies transmission, and virulence. Our findings provide evidence that significant selection pressures have shaped the evolution of the H5N1 clade 2.3.4.4b viruses, facilitating their expanded host tropism and the potential for further adaptation to mammalian hosts. We discuss how exogenous factors, such as altered bird migration patterns and increased host susceptibility, may have contributed to the expanded host range. As H5N1 viruses continue to infect new hosts, there is a greater risk of emergent novel variants with increased pathogenicity in humans and animals. Thus, comprehensive One Health surveillance is critical to monitor transmission among avian and mammalian hosts.

Genetic diversification patterns in swine influenza A virus (H1N2) in vaccinated and nonvaccinated animals

Influenza A viruses (IAVs) are characterized by having a segmented genome, low proofreading polymerases, and a wide host range. Consequently, IAVs are constantly evolving in nature causing a threat to animal and human health. In 2009 a new human pandemic IAV strain arose in Mexico because of a reassortment between two strains previously circulating in pigs; Eurasian "avian-like" (EA) swine H1N1 and "human-like" H1N2, highlighting the importance of swine as adaptation host of avian to human IAVs. Nowadays, although of limited use, a trivalent vaccine, which include in its formulation H1N1, H3N2, and, H1N2 swine IAV (SIAV) subtypes, is one of the most applied strategies to reduce SIAV circulation in farms. Protection provided by vaccines is not complete, allowing virus circulation, potentially favoring viral evolution. The evolutionary dynamics of SIAV quasispecies were studied in samples collected at different times from 8 vaccinated and 8 nonvaccinated pigs, challenged with H1N2 SIAV. In total, 32 SIAV genomes were sequenced by next-generation sequencing, and subsequent variant-calling genomic analysis was carried out. Herein, a total of 364 de novo single nucleotide variants (SNV) were found along all genetic segments in both experimental groups. The nonsynonymous substitutions proportion found was greater in vaccinated animals suggesting that H1N2 SIAV was under positive selection in this scenario. The impact of each substitution with an allele frequency greater than 5% was hypothesized according to previous literature, particularly in the surface glycoproteins hemagglutinin and neuraminidase. The H1N2 SIAV quasispecies evolution capacity was evidenced, observing different evolutionary trends in vaccinated and nonvaccinated animals.

Vaccination against swine influenza in pigs causes different drift evolutionary patterns upon swine influenza virus experimental infection and reduces the likelihood of genomic reassortments

Influenza A viruses (IAVs) can infect a wide variety of bird and mammal species. Their genome is characterized by 8 RNA single stranded segments. The low proofreading activity of their polymerases and the genomic reassortment between different IAVs subtypes allow them to continuously evolve, constituting a constant threat to human and animal health. In 2009, a pandemic of an IAV highlighted the importance of the swine host in IAVs adaptation between humans and birds. The swine population and the incidence of swine IAV is constantly growing. In previous studies, despite vaccination, swine IAV growth and evolution were proven in vaccinated and challenged animals. However, how vaccination can drive the evolutionary dynamics of swine IAV after coinfection with two subtypes is poorly studied. In the present study, vaccinated and nonvaccinated pigs were challenged by direct contact with H1N1 and H3N2 independent swine IAVs seeder pigs. Nasal swab samples were daily recovered and broncho-alveolar lavage fluid (BALF) was also collected at necropsy day from each pig for swine IAV detection and whole genome sequencing. In total, 39 swine IAV whole genome sequences were obtained by next generation sequencing from samples collected from both experimental groups. Subsequently, genomic, and evolutionary analyses were carried out to detect both, genomic reassortments and single nucleotide variants (SNV). Regarding the segments found per sample, the simultaneous presence of segments from both subtypes was much lower in vaccinated animals, indicating that the vaccine reduced the likelihood of genomic reassortment events. In relation to swine IAV intra-host diversity, a total of 239 and 74 SNV were detected within H1N1 and H3N2 subtypes, respectively. Different proportions of synonymous and nonsynonymous substitutions were found, indicating that vaccine may be influencing the main mechanism that shape swine IAV evolution, detecting natural, neutral, and purifying selection in the different analyzed scenarios. SNV were detected along the whole swine IAV genome with important nonsynonymous substitutions on polymerases, surface glycoproteins and nonstructural proteins, which may have an impact on virus replication, immune system escaping and virulence of virus, respectively. The present study further emphasized the vast evolutionary capacity of swine IAV, under natural infection and vaccination pressure scenarios.

Vaccination decreases the risk of influenza A virus reassortment but not genetic variation in pigs

Although vaccination is broadly used in North American swine breeding herds, managing swine influenza is challenging primarily due to the continuous evolution of influenza A virus (IAV) and the ability of the virus to transmit among vaccinated pigs. Studies that have simultaneously assessed the impact of vaccination on the emergence of IAV reassortment and genetic variation in pigs are limited. Here, we directly sequenced 28 bronchoalveolar lavage fluid (BALF) samples collected from vaccinated and unvaccinated pigs co-infected with H1N1 and H3N2 IAV strains, and characterized 202 individual viral plaques recovered from 13 BALF samples. We identified 54 reassortant viruses that were grouped in 17 single and 16 mixed genotypes. Notably, we found that prime-boost vaccinated pigs had less reassortant viruses than nonvaccinated pigs, likely due to a reduction in the number of days pigs were co-infected with both challenge viruses. However, direct sequencing from BALF samples revealed limited impact of vaccination on viral variant frequency, evolutionary rates, and nucleotide diversity in any IAV coding regions. Overall, our results highlight the value of IAV vaccination not only at limiting virus replication in pigs but also at protecting public health by restricting the generation of novel reassortants with zoonotic and/or pandemic potential.

Swine influenza A viruses cause severe illness among pigs and financial losses on pig farms worldwide. These viruses can also infect humans and have caused deadly human pandemics in the past. Influenza A viruses are dangerous because viruses can be transferred between humans, birds and pigs. These co-infections can allow the viruses to swap genetic material. Viral genetic exchanges can result in new virus strains that are more dangerous or that can infect other types of animals more easily.

Farmers vaccinate their pigs to control the swine influenza A virus. The vaccines are regularly updated to match circulating virus strains. But the virus evolves rapidly to escape vaccine-induced immunity, and infections are common even in vaccinated pigs. Learning about how vaccination affects the evolution of influenza A viruses in pigs could help scientists prevent outbreaks on pig farms and avoid spillover pandemics in humans.

Li et al. show that influenza A viruses are less likely to swap genetic material in vaccinated and boosted pigs than in unvaccinated animals. In the experiments, Li et al. collected swine influenza A samples from the lungs of pigs that had received different vaccination protocols. Next, Li et al. used next-generation sequencing to identify new mutations in the virus or genetic swaps among different strains. In pigs infected with both the H1N1 and H3N2 strains of influenza, the two viruses began trading genes within a week. But less genetic mixing occurred in vaccinated and boosted pigs because they spent less time infected with both viruses than in unvaccinated pigs. The vaccination status of the pig did not have much effect on how many new mutations occurred in the viruses.

The experiments show that vaccinating and boosting pigs against influenza A viruses may protect against genetic swapping among influenza viruses. If future studies on pig farms confirm the results, the information gleaned from the study could help scientists improve farm vaccine protocols to further reduce influenza risks to animals and people.

Molecular Epidemiology and Evolutionary Analysis of Avian Influenza A(H5) Viruses Circulating in Egypt, 2019–2021

The highly pathogenic avian influenza (HPAI) H5N8 virus was first detected in Egypt in late 2016. Since then, the virus has spread rapidly among different poultry sectors, becoming the dominant HPAI H5 subtype reported in Egypt. Different genotypes of the HPAI H5N8 virus were reported in Egypt; however, the geographic patterns and molecular evolution of the Egyptian HPAI H5N8 viruses are still unclear. Here, extensive epidemiological surveillance was conducted, including more than half a million samples collected from different poultry sectors (farms/backyards/live bird markets) from all governorates in Egypt during 2019–2021. In addition, genetic characterization and evolutionary analyses were performed using 47 selected positive H5N8 isolates obtained during the same period. The result of the conducted surveillance showed that HPAI H5N8 viruses of clade 2.3.4.4b continue to circulate in different locations in Egypt, with an obvious seasonal pattern, and no further detection of the HPAI H5N1 virus of clade 2.2.1.2 was observed in the poultry population during 2019–2021. In addition, phylogenetic and Bayesian analyses revealed that two major genotypes (G5 and G6) of HPAI H5N8 viruses were continually expanding among the poultry sectors in Egypt. Notably, molecular dating analysis suggested that the Egyptian HPAI H5N8 virus is the potential ancestral viruses of the European H5N8 viruses of 2020–2021. In summary, the data of this study highlight the current epidemiology, diversity, and evolution of HPAI H5N8 viruses in Egypt and call for continuous monitoring of the genetic features of the avian influenza viruses in Egypt.

A retrospective comparison of COVID-19 and seasonal influenza mortality and outcomes in the ICUs of a French university hospital

BACKGROUND

SARS-Cov-2 (COVID-19) has become a major worldwide health concern since its appearance in China at the end of 2019.

OBJECTIVE

To evaluate the intrinsic mortality and burden of COVID-19 and seasonal influenza pneumonia in ICUs in the city of Lyon, France.

DESIGN

A retrospective study.

SETTING

Six ICUs in a single institution in Lyon, France.

PATIENTS

Consecutive patients admitted to an ICU with SARS-CoV-2 pneumonia from 27 February to 4 April 2020 (COVID-19 group) and seasonal influenza pneumonia from 1 November 2015 to 30 April 2019 (influenza group). A total of 350 patients were included in the COVID-19 group (18 refused to consent) and 325 in the influenza group (one refused to consent). Diagnosis was confirmed by RT-PCR. Follow-up was completed on 1 April 2021.

MAIN OUTCOMES AND MEASURES

Differences in 90-day adjusted-mortality between the COVID-19 and influenza groups were evaluated using a multivariable Cox proportional hazards model.

RESULTS

COVID-19 patients were younger, mostly men and had a higher median BMI, and comorbidities, including immunosuppressive condition or respiratory history were less frequent. In univariate analysis, no significant differences were observed between the two groups regarding in-ICU mortality, 30, 60 and 90-day mortality. After Cox modelling adjusted on age, sex, BMI, cancer, sepsis-related organ failure assessment (SOFA) score, simplified acute physiology score SAPS II score, chronic obstructive pulmonary disease and myocardial infarction, the probability of death associated with COVID-19 was significantly higher in comparison to seasonal influenza [hazard ratio 1.57, 95% CI (1.14 to 2.17); P = 0.006]. The clinical course and morbidity profile of both groups was markedly different; COVID-19 patients had less severe illness at admission (SAPS II score, 37 [28 to 48] vs. 48 [39 to 61], P < 0.001 and SOFA score, 4 [2 to 8] vs. 8 [5 to 11], P < 0.001), but the disease was more severe considering ICU length of stay, duration of mechanical ventilation, PEEP level and prone positioning requirement.

CONCLUSION

After ICU admission, COVID-19 was associated with an increased risk of death compared with seasonal influenza. Patient characteristics, clinical course and morbidity profile of these diseases is markedly different.

Evolutionary conservation and positive selection of Influenza A Nucleoprotein CTL epitopes for universal vaccination: a proof-of-concept

Influenza (flu) infection is a leading cause of respiratory disease and death worldwide. While seasonal flu vaccines are effective at reducing morbidity and mortality, such effects rely on the odds of successful prediction of the upcoming viral strains. Additional threats from emerging flu viruses that we cannot predict and avian flu viruses that can be directly transmitted to humans, urge the strategic development of universal vaccinations that can protect against flu viruses of different subtypes and across species. Annual flu vaccines elicit mainly humoral responses. Under circumstances when antibodies induced by vaccination fail to recognize and neutralize the emerging virus adequately, virus-specific cytotoxic T lymphocytes (CTLs) are the major contributors to the control of viral replication and elimination of infected cells. Our studies exploited the evolutionary conservation of influenza A nucleoprotein (NP) and the fact that NP-specific CTL responses pose a constant selecting pressure on functional CTL epitopes, to screen for NP epitopes that are highly conserved among heterosubtypes but are subjected to positive selection historically. We identified a region on NP that is evolutionarily conserved and historically positively selected (NP_137-182 ) and validated that it contains an epitope that is functional in eliciting NP-specific CTL responses and immunity that can partially protect immunized mice against lethal dose infection of a heterosubtypic influenza A virus. Our proof-of-concept study supports the hypothesis that evolutionary conservation and positive selection of influenza nucleoprotein can be exploited to identify functional CTL epitope to elicit cross protection against different heterosubtypes, therefore, to help develop strategies to modify flu vaccine formula for a broader and more durable protective immunity. This article is protected by copyright. All rights reserved.

Identification and Characterization of Swine Influenza Virus H1N1 Variants Generated in Vaccinated and Nonvaccinated, Challenged Pigs

Influenza viruses represent a continuous threat to both animal and human health. The 2009 H1N1 A influenza pandemic highlighted the importance of a swine host in the adaptation of influenza viruses to humans. Nowadays, one of the most extended strategies used to control swine influenza viruses (SIVs) is the trivalent vaccine application, whose formulation contains the most frequently circulating SIV subtypes H1N1, H1N2, and H3N2. These vaccines do not provide full protection against the virus, allowing its replication, evolution, and adaptation. To better understand the main mechanisms that shape viral evolution, here, the SIV intra-host diversity was analyzed in samples collected from both vaccinated and nonvaccinated animals challenged with the H1N1 influenza A virus. Twenty-eight whole SIV genomes were obtained by next-generation sequencing, and differences in nucleotide variants between groups were established. Substitutions were allocated along all influenza genetic segments, while the most relevant nonsynonymous substitutions were allocated in the NS1 protein on samples collected from vaccinated animals, suggesting that SIV is continuously evolving despite vaccine application. Moreover, new viral variants were found in both vaccinated and nonvaccinated pigs, showing relevant substitutions in the HA, NA, and NP proteins, which may increase viral fitness under field conditions.

Substantial Antigenic Drift in the Hemagglutinin Protein of Swine Influenza A Viruses

The degree of antigenic drift in swine influenza A viruses (swIAV) has historically been regarded as minimal compared to that of human influenza A virus strains. However, as surveillance activities on swIAV have increased, more isolates have been characterized, revealing a high level of genetic and antigenic differences even within the same swIAV lineage. The objective of this study was to investigate the level of genetic drift in one enzootically infected swine herd over one year. Nasal swabs were collected monthly from sows (n = 4) and piglets (n = 40) in the farrowing unit, and from weaners (n = 20) in the nursery. Virus from 1-4 animals were sequenced per month. Analyses of the sequences revealed that the hemagglutinin (HA) gene was the main target for genetic drift with a substitution rate of 7.6 × 10-3 substitutions/site/year and evidence of positive selection. The majority of the mutations occurred in the globular head of the HA protein and in antigenic sites. The phylogenetic tree of the HA sequences displayed a pectinate typology, where only a single lineage persists and forms the ancestor for subsequent lineages. This was most likely caused by repeated selection of a single immune-escape variant, which subsequently became the founder of the next wave of infections.

A prime-boost concept using a T-cell epitope-driven DNA vaccine followed by a whole virus vaccine effectively protected pigs in the pandemic H1N1 pig challenge model

Influenza A virus (IAV) vaccines in pigs generally provide homosubtypic protection but fail to prevent heterologous infections. In this pilot study, the efficacy of an intradermal pDNA vaccine composed of conserved SLA class I and class II T cell epitopes (EPITOPE) against a homosubtypic challenge was compared to an intramuscular commercial inactivated whole virus vaccine (INACT) and a heterologous prime boost approach using both vaccines. Thirty-nine IAV-free, 3-week-old pigs were randomly assigned to one of five groups including NEG-CONTROL (unvaccinated, sham-challenged), INACT-INACT-IAV (vaccinated with FluSure XP® at 4 and 7 weeks, pH1N1 challenged), EPITOPE-INACT-IAV (vaccinated with PigMatrix EDV at 4 and FluSure XP® at 7 weeks, pH1N1 challenged), EPITOPE-EPITOPE-IAV (vaccinated with PigMatrix EDV at 4 and 7 weeks, pH1N1 challenged), and a POS-CONTROL group (unvaccinated, pH1N1 challenged). The challenge was done at 9 weeks of age and pigs were necropsied at day post challenge (dpc) 5. At the time of challenge, all INACT-INACT-IAV pigs, and by dpc 5 all EPITOPE-INACT-IAV pigs were IAV seropositive. IFNγ secreting cells, recognizing vaccine epitope-specific peptides and pH1N1 challenge virus were highest in the EPITOPE-INACT-IAV pigs at challenge. Macroscopic lung lesion scores were reduced in all EPITOPE-INACT-IAV pigs while INACT-INACT-IAV pigs exhibited a bimodal distribution of low and high scores akin to naïve challenged animals. No IAV antigen in lung tissues was detected at necropsy in the EPITOPE-INACT-IAV group, which was similar to naïve unchallenged pigs and different from all other challenged groups. Results suggest that the heterologous prime boost approach using an epitope-driven DNA vaccine followed by an inactivated vaccine was effective against a homosubtypic challenge, and further exploration of this vaccine approach as a practical control measure against heterosubtypic IAV infections is warranted.

Molecular footprints of selective pressure in the neuraminidase gene of currently circulating human influenza subtypes and lineages

Influenza neuraminidase (NA) is under selective pressure (SP) of both host immune system and drug use. Here, we assembled large datasets of NA sequences of worldwide circulating viruses to estimate the global and site-specific SP acting on all current subtypes/lineages of human influenza NA. An overall negative SP of similar magnitude and a prevalence of negatively selected sites were observed for all subtypes/lineages. Positively selected sites varied according to the subtype/lineage, including N1-NA sites 247 and 275, N2-NA sites 148 and 151, and B/Victoria-NA site 395 associated with drug-resistance or reduced susceptibility. These results evidenced a potential role of positive selection in the low-level spread of A(H1N1)pdm09-H275Y drug-resistant viruses, and alerted for a potential higher risk of spread of a synergistic A(H1N1)pdm09 drug-resistant variant (H275Y/S247N). The positive selection detected at N2-NA sites 148 and 151 was probably an artefact from cell-culture. Overall mapping revealed six potential new druggable regions.

Noninvasive vaccination against infectious diseases

The development of a successful vaccine, which should elicit a combination of humoral and cellular responses to control or prevent infections, is the first step in protecting against infectious diseases. A vaccine may protect against bacterial, fungal, parasitic, or viral infections in animal models, but to be effective in humans there are some issues that should be considered, such as the adjuvant, the route of vaccination, and the antigen-carrier system. While almost all licensed vaccines are injected such that inoculation is by far the most commonly used method, injection has several potential disadvantages, including pain, cross contamination, needlestick injury, under- or overdosing, and increased cost. It is also problematic for patients from rural areas of developing countries, who must travel to a hospital for vaccine administration. Noninvasive immunizations, including oral, intranasal, and transcutaneous administration of vaccines, can reduce or eliminate pain, reduce the cost of vaccinations, and increase their safety. Several preclinical and clinical studies as well as experience with licensed vaccines have demonstrated that noninvasive vaccine immunization activates cellular and humoral immunity, which protect against pathogen infections. Here we review the development of noninvasive immunization with vaccines based on live attenuated virus, recombinant adenovirus, inactivated virus, viral subunits, virus-like particles, DNA, RNA, and antigen expression in rice in preclinical and clinical studies. We predict that noninvasive vaccine administration will be more widely applied in the clinic in the near future.

Neuraminidase inhibitor susceptibility and neuraminidase enzyme kinetics of human influenza A and B viruses circulating in Thailand in 2010-2015

Amino acid substitutions within or near the active site of the viral neuraminidase (NA) may affect influenza virus fitness. In influenza A(H3N2) and B viruses circulating in Thailand between 2010 and 2015, we identified several NA substitutions that were previously reported to be associated with reduced inhibition by NA inhibitors (NAIs). To study the effect of these substitutions on the enzymatic properties of NA and on virus characteristics, we generated recombinant influenza viruses possessing either a wild type (WT) NA or an NA with a single I222V, S331G, or S331R substitution [in influenza A(H3N2) viruses] or a single D342S, A395T, A395V, or A395D NA substitution (in influenza B viruses). We generated recombinant (7:1) influenza A and B viruses on the genetic background of A/Puerto Rico/8/1934 (A/PR/8, H1N1) or B/Yamanashi/166/1998 (B/YAM) viruses, respectively. In contrast to the expected phenotypes, all the recombinant influenza A(H3N2) and B viruses carrying putative NA resistance substitutions were susceptible to NAIs. The Km and Vmax for the NAs of A/PR8-S331G and A/PR8-S331R viruses were higher than for the NA of WT virus, and the corresponding values for the B/YAM-D342S virus were lower than for the NA of WT virus. Although there was initial variation in the kinetics of influenza A and B viruses' replication in MDCK cells, their titers were comparable to each other and to WT viruses at later time points. All introduced substitutions were stable except for B/YAM-D342S and B/YAM-A395V which reverted to WT sequences after three passages. Our data suggest that inferring susceptibility to NAIs based on sequence information alone should be cautioned. The impact of NA substitution on NAI resistance, viral growth, and enzymatic properties is viral context dependent and should be empirically determined.

Adaptive molecular evolution of MC1R gene reveals the evidence for positive diversifying selection in indigenous goat populations

Detecting signatures of selection can provide a new insight into the mechanism of contemporary breeding and artificial selection and further reveal the causal genes associated to the phenotypic variation. However, the signatures of selection on genes entailing for profitable traits between Chinese commercial and indigenous goats have been poorly interpreted. We noticed footprints of positive selection at MC1R gene containing SNPs genotyped in five Chinese native goat breeds. An experimental distribution of F_ST was built based on approximations of F_ST for each SNP across five breeds. We identified selection using the high F_ST outlier method and found that MC1R candidate gene show evidence of positive selection. Furthermore, adaptive selection pressure on specific codons was determined using different codon based on maximum‐likelihood methods; signature of positive selection in mammalian MC1R was explored in individual codons. Evolutionary analyses were inferred under maximum likelihood models, the HyPhy package implemented in the DATAMONKEY Web Server. The results of codon selection displayed positive diversifying selection at the sites were mainly involved in development of genetic variations in coat color in various mammalian species. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat MC1R provides new insights that the gene evolution may have been modulated by domestication events in goats.

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.

Molecular characterization and phylogenetic analysis of human influenza A viruses isolated in Iran during the 2014-2015 season

Influenza A viruses are an important cause of severe infectious diseases in humans and are characterized by their fast evolution rate. Global monitoring of these viruses is critical to detect newly emerging variants during annual epidemics. Here, we sought to genetically characterize influenza A/H1N1pdm09 and A/H3N2 viruses collected in Iran during the 2014-2015 influenza season. A total of 200 nasopharyngeal swabs were collected from patients with influenza-like illnesses. Swabs were screened for influenza A and B using real-time PCR. Furthermore, positive specimens with high virus load underwent virus isolation and genetic characterization of their hemagglutinin (HA) and M genes. Of the 200 specimens, 80 were influenza A-positive, including 44 A/H1N1pdm09 and 36 A/H3N2, while 18 were influenza B-positive. Phylogenetic analysis of the HA genes of the A/H1N1pdm09 viruses revealed the circulation of clade 6C, characterized by amino acid substitutions D97N, V234I and K283E. Analysis of the A/H3N2 viruses showed a genetic drift from the vaccine strain A/Texas/50/2012 with 5 mutations (T128A, R142G, N145S, P198S and S219F) belonging to the antigenic sites A, B, and D of the HA protein. The A/H3N2 viruses belonged to phylogenetic clades 3C.2 and 3C.3. The M gene trees of the Iranian A/H1N1pdm09 and A/H3N2 mirrored the clustering patterns of their corresponding HA trees. Our results reveal co-circulation of several influenza A virus strains in Iran during the 2014-2015 influenza season.

Evolutionary and functional implications of hypervariable loci within the skin virome

Localized genomic variability is crucial for the ongoing conflicts between infectious microbes and their hosts. An understanding of evolutionary and adaptive patterns associated with genomic variability will help guide development of vaccines and antimicrobial agents. While most analyses of the human microbiome have focused on taxonomic classification and gene annotation, we investigated genomic variation of skin-associated viral communities. We evaluated patterns of viral genomic variation across 16 healthy human volunteers. Human papillomavirus (HPV) and Staphylococcus phages contained 106 and 465 regions of diversification, or hypervariable loci, respectively. Propionibacterium phage genomes were minimally divergent and contained no hypervariable loci. Genes containing hypervariable loci were involved in functions including host tropism and immune evasion. HPV and Staphylococcus phage hypervariable loci were associated with purifying selection. Amino acid substitution patterns were virus dependent, as were predictions of their phenotypic effects. We identified diversity generating retroelements as one likely mechanism driving hypervariability. We validated these findings in an independently collected skin metagenomic sequence dataset, suggesting that these features of skin virome genomic variability are widespread. Our results highlight the genomic variation landscape of the skin virome and provide a foundation for better understanding community viral evolution and the functional implications of genomic diversification of skin viruses.

Genome Hotspots for Nucleotide Substitutions and the Evolution of Influenza A (H1N1) Human Strains

In recent years a number of studies have brought attention to the role of positive selection during the evolution of antigenic escape by influenza strains. Particularly, the identification of positively selected sites within antigenic domains of viral surface proteins has been used to suggest that the evolution of viral–host receptor binding specificity is driven by selection. Here we show that, following the 1918 outbreak, the antigenic sites of the hemagglutinin (HA) viral surface protein and the stalk region of neuraminidase became substitution hotspots. The hotspots show similar patterns of nucleotide substitution bias at synonymous and nonsynonymous sites. Such bias imposes directionality in amino acid replacements that can influence signals of selection at antigenic sites. Our results suggest that the high accumulation of substitutions within the antigenic sites of HA can explain not only cases of antigenic escape by antigenic drift but also lead to occasional episodes of viral extinction.

Dynamic reassortments and genetic heterogeneity of the human-infecting influenza A (H7N9) virus

Influenza A (H7N9) virus has been causing human infections in China since February 2013, raising serious concerns of potential pandemics. Previous studies demonstrate that human infection is directly linked to live animal markets, and that the internal genes of the virus are derived from H9N2 viruses circulating in the Yangtze River Delta area in Eastern China. Here following analysis of 109 viruses, we show a much higher genetic heterogeneity of the H7N9 viruses than previously reported, with a total of 27 newly designated genotypes. Phylogenetic and genealogical inferences reveal that genotypes G0 and G2.6 dominantly co-circulate within poultry, with most human isolates belonging to the genotype G0. G0 viruses are also responsible for the inter- and intra-province transmissions, leading to the genesis of novel genotypes. These observations suggest the province-specific H9N2 virus gene pools increase the genetic diversity of H7N9 via dynamic reassortments and also imply that G0 has not gained overwhelming fitness and the virus continues to undergo reassortment.

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.

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.

Phylogenetic analysis of the neuraminidase gene of pandemic H1N1 influenza A virus circulating in the South American region

Molecular characterization of circulating influenza A viruses (IAV) in all regions of the world is essential to detect mutations potentially involved in increased virulence, anti-viral resistance and immune escape. In order to gain insight into these matters, a phylogenetic analysis of the neuraminidase (NA) gene of 146 pandemic H1N1 (H1N1pdm) influenza A virus strains isolated in Argentina, Brazil, Chile, Paraguay, Peru and Uruguay from 2009 to 2013 was performed. Comparison of vaccine strain A/California/7/2009 included in the influenza vaccine recommended for the Southern hemisphere from 2010 through 2013 influenza seasons and strains isolated in South America revealed several amino acid substitutions. Mapping of these substitutions revealed that most of them are located at the surface of the protein and do not interfere with the active site. 3.4% of the strains enrolled in these studies carried the H275Y substitution that confers resistance to oseltamivir. Strains isolated in South America differ from vaccine in two predicted B-cell epitope regions present at positions 102-103 and 351-352 of the NA protein. Moreover, vaccine and strains isolated in Paraguay differ also in an epitope present at position 229. These differences among strains isolated in South America and vaccine strain suggests that these epitopes may not be present in strains isolated in this region. A potential new N-linked glycosylation site was observed in the NA protein of an H1N1pdm IAV strain isolated in Brazil. The results of these studies revealed several genetic and antigenic differences in the NA of H1N1pdm IAV among vaccine and strains circulating in South America. All these findings contribute to our understanding of the course of genetic and antigenic evolution of H1N1pdm IAV populations circulating in the South American region and, consequently, contribute to the study and selection of future and more appropriate vaccines and anti-viral drugs.

Adaptive evolution of a novel avian-origin influenza A/H7N9 virus

In China, the recent outbreak of novel influenza A/H7N9 virus has been assumed to be severe, and it may possibly turn brutal in the near future. In order to develop highly protective vaccines and drugs for the A/H7N9 virus, it is critical to find out the selection pressure of each amino acid site. In the present study, six different statistical methods consisting of four independent codon-based maximum likelihood (CML) methods, one hierarchical Bayesian (HB) method and one branch-site (BS) method, were employed to determine if each amino acid site of A/H7N9 virus is under natural selection pressure. Functions for both positively and negatively selected sites were inferred by annotating these sites with experimentally verified amino acid sites. Comprehensively, the single amino acid site 627 of PB2 protein was inferred as positively selected and it function was identified as a T-cell epitope (TCE). Among the 26 negatively selected amino acid sites of PB2, PB1, PA, HA, NP, NA, M1 and NS2 proteins, only 16 amino acid sites were identified to be involved in TCEs. In addition, 7 amino acid sites including, 608 and 609 of PA, 480 of NP, and 24, 25, 109 and 205 of M1, were identified to be involved in both B-cell epitopes (BCEs) and TCEs. Conversely, the function of positions 62 of PA, and, 43 and 113 of HA was unknown. In conclusion, the seven amino acid sites engaged in both BCEs and TCEs were identified as highly suitable targets, as these sites will be predicted to play a principal role in inducing strong humoral and cellular immune responses against A/H7N9 virus.

Mutations associated with severity of the pandemic influenza A(H1N1)pdm09 in humans: a systematic review and meta-analysis of epidemiological evidence

Mutations in the haemagglutinin (HA), non-structural protein 1 (NS1) and polymerase basic protein 2 (PB2) of influenza viruses have been associated with virulence. This study investigated the association between mutations in these genes in influenza A(H1N1)pdm09 virus and the risk of severe or fatal disease. Searches were conducted on the MEDLINE, EMBASE and Web of Science electronic databases and the reference lists of published studies. The PRISMA and STROBE guidelines were followed in assessing the quality of studies and writing-up. Eighteen (18) studies, from all continents, were included in the systematic review (recruiting patients 0 - 77 years old). The mutation D222G was associated with a significant increase in severe disease (pooled RD: 11 %, 95 % CI: 3.0 % - 18.0 %, p = 0.004) and the risk of fatality (RD: 23 %, 95 % CI: 14.0 %-31.0 %, p = < 0.0001). No association was observed between the mutations HA-D222N, D222E, PB2-E627K and NS1-T123V and severe/fatal disease. The results suggest that no virus quasispecies bearing virulence-conferring mutations in the HA, PB2 and NS1 predominated. However issues of sampling bias, and bias due to uncontrolled confounders such as comorbidities, and viral and bacterial coinfection, should be born in mind. Influenza A viruses should continue to be monitored for the occurrence of virulence-conferring mutations in HA, PB2 and NS1. There are suggestions that respiratory virus coinfections also affect virus virulence. Studies investigating the role of genetic mutations on disease outcome should make efforts to also investigate the role of respiratory virus coinfections.

Genetic dynamic analysis of the influenza A H5N1 NS1 gene in China

The direct precursors of the A/Goose/Guangdong/1/1996 (GS/GD) virus lineage and its reassortants have been established geographically and ecologically. To investigate the variation and evolutionary dynamics of H5N1 viruses, whole-genome viral sequences (n = 164) were retrieved from the NCBI Influenza Virus Resource. Here, we present phylogenetic evidence for intrasubtype reassortments among H5N1 viruses isolated from China during 1996–2012. On the basis of phylogenetic analysis, we identified four major groups and further classified the reassortant viruses into three subgroups. Putative mosaic structures were mostly found in the viral ribonucleoprotein (vRNP) complexes and 91.0% (10/11) mosaics were obtained from terrestrial birds. Sequence variability and selection pressure analyses revealed that both surface glycoproteins (HA and NA) and nonstructural protein 1 (NS1) have higher dN/dS ratio and variability than other internal proteins. Furthermore, we detected 47 positively selected sites in genomic segments with the exception of PB2 and M1 genes. Hemagglutinin (HA) and neuraminidase (NA) are considered highly variable due to host immune pressure, however, it is not known what drives NS1 variability. Therefore, we performed a thorough analysis of the genetic variation and selective pressure of NS1 protein (462 available NS1 sequences). We found that most of positively selected sites and variable amino acids were located in the C-terminal effector domain (ED) of NS1. In addition, we focused on the NS1–RNA and NS1–protein interactions that were involved in viral replication mechanisms and host immune response. Transcriptomic analysis of H5N1-infected monkey lungs showed that certain PI3K-related genes were up-regulated.

Virus-specific factors associated with zoonotic and pandemic potential

Influenza A is a highly contagious respiratory virus in constant evolution and represents a threat to both veterinary and human public health. IA viruses (IAVs) originate in avian reservoirs but may adapt to humans, either directly or through the spillover to another mammalian species, to the point of becoming pandemic. IAVs must successfully be able to (i) transmit from animal to human, (ii) interact with host cells, and (iii) transmit from human to human. The mechanisms by which viruses evolve, cause zoonotic infections, and adapt to a new host species are indeed complex and appear to be a heterogeneous collection of viral evolutionary events rather than a single phenomenon. Progress has been made in identifying some of the genetic markers mainly associated with virulence and transmission; this achievement has improved our knowledge of how to manage a pandemic event and of how to identify IAVs with pandemic potential. Early evidence of emerging viruses and surveillance of animal IAVs is made possible only by strengthening the collaboration between the public and veterinary health sectors.

Kinetics of pulmonary immune cells, antibody responses and their correlations with the viral clearance of influenza A fatal infection in mice

Fatal influenza A virus infection is a major threat to public health throughout the world. Lung macrophages and neutrophils have critical roles for both the pathogenesis and viral clearance of fatal viral infections. These are complicated by the interaction of innate immunity and adaptive immunity against viral infection. In this study, we investigated the overall kinetics of lung macrophages, neutrophils, CD4⁺T cells, CD8⁺T cells, CD38⁺ cells, and CD138⁺ cells, the levels of antibody and cytokine responses, both in the early and late phases of fatal infection with A/PR/8/34 (H1N1) virus in mice. The changes in lung viral load were also evaluated. We found that pulmonary macrophages and neutrophils both accumulated in the early and late phases of fatal infections and they positively correlated with the lung and serum antibody titers, and negatively correlated with the viral load locally. The secretion of IL-6 might relate to high numbers of macrophages and neutrophils in the early infection. The work implies that pulmonary macrophages, neutrophils and the antibody response all have an essential role in virus elimination of fatal influenza A viral infection. These findings may have implications for the development of prophylactic and therapeutic strategies in fatal influenza A viral infection. Further evaluation of the cooperation among macrophages, neutrophils and antibody responses in eliminating the virus with fatal infection is needed.

Swine influenza: clinical, serological, pathological, and virological cross-sectional studies in nine farms in Argentina

BACKGROUND

Influenza A viruses (IAV) are important pathogens responsible for economic losses in the swine industry and represent a threat to public health. In Argentina, clinical, pathological, and virological findings suggest that IAV infection is widespread among pig farms. In addition, several subtypes of IAV, such as pH1N1, H3N2, δ1H1N1, and δ2H1N2, have been reported.

OBJECTIVES

To evaluate the infection patterns of influenza virus in nine pig farms in Argentina.

METHODS

Clinical, serological, pathological, and virological cross-sectional studies were conducted.

RESULTS

Clinical and pathological results were characteristic of endemic influenza infection in eight of the nine farms studied. By rRT-PCR, six of the nine farms were positive to influenza. Five IAV were obtained. Genome analysis determined that four of the isolations were pH1N1 and that the remaining one was a reassortant human origin H3N2 virus containing pandemic internal genes. Serological results showed that all farms were positive to influenza A antibodies. Moreover, the hemagglutination inhibition test showed that infection with viruses containing HA's from different subtypes (pH1, δ1H1, δ2H1, and H3) is present among the farms studied and that coinfections with two or more subtypes were present in 80.5% of positive pigs.

CONCLUSIONS

Because vaccines against IAV are not licensed in Argentina, these results reflect the situation of IAV infection in non-vaccinated herds. This study provides more information about the circulation and characteristics of IAV in a poorly surveyed region. This study provides more data that will be used to evaluate the tools necessary to control this disease.

Different evolutionary trends of swine H1N2 influenza viruses in Italy compared to European viruses

European H1N2 swine influenza viruses (EU H1N2SIVs) arose from multiple reassortment events among human H1N1, human H3N2, and avian influenza viruses. We investigated the evolutionary dynamics of 53 Italian H1N2 strains by comparing them with EU H1N2 SIVs. Hemagglutinin (HA) phylogeny revealed Italian strains fell into four groups: Group A and B (41 strains) had a human H1 similar to EU H1N2SIVs, which probably originated in 1986. However Group B (38 strains) formed a subgroup that had a two-amino acid deletion at positions 146/147 in HA. Group C (11 strains) contained an avian H1 that probably originated in 1996, and Group D (1 strain) had an H1 characteristic of the 2009 pandemic strain. Neuraminidase (NA) phylogeny suggested a series of genomic reassortments had occurred. Group A had an N2 that originated from human H3N2 in the late 1970s. Group B had different human N2 that most likely arose from a reassortment with the more recent human H3N2 virus, which probably occurred in 2000. Group C had an avian-like H1 combined with an N2 gene from one of EU H1N2SIVs, EU H3N2SIVs or Human H3N2. Group D was part of the EU H3N2SIVs clade. Although selection pressure for HA and NA was low, several positively selected sites were identified in both proteins, some of which were antigenic, suggesting selection influenced the evolution of SIV. The data highlight different evolutionary trends between European viruses and currently circulating Italian B strains and show the establishment of reassortant strains involving human viruses in Italian pigs.

Genomic analysis of pandemic and post-pandemic influenza A pH1N1 viruses isolated in Rio Grande do Sul, Brazil

During the 2009 influenza A pH1N1 pandemics in Brazil, the state that was most affected was Rio Grande do Sul (RS), with over 3,000 confirmed cases, including 298 deaths. While no cases were confirmed in 2010, 103 infections with 14 deaths by pH1N1 were reported in 2011. Genomic analysis of the circulating viruses is fundamental for understanding viral evolution and supporting vaccine development against these pathogens. This study investigated whole genomes of six pH1N1 virus isolates from pandemic and post-pandemic periods in RS, Brazil. Phylogenetic analysis using the concatenated genome segments demonstrated that at least two lineages of the virus co-circulated in RS during the 2009 pandemic period. Moreover, our analysis showed that the post-pandemic pH1N1 virus from 2011 constitutes a distinct clade whose ancestor belongs to clade 7. All six isolates contained amino acid substitutions in their proteins when compared to the archetype strains California/04/2009 and California/07/2009. The 2011 isolates contained more amino acid substitutions, and most of their genes were under purifying selection. Based on the amino acid substitutions in HA epitopes from strains isolated in RS, Brazil, in silico analysis predicted a decrease in vaccine efficacy against post-pandemic strains (median 31.562 %) in relation to pandemic ones (median 39.735 %).

Evidence of selection pressures of neuraminidase gene (NA) of influenza A virus subtype H5N1 on different hosts in Guangxi Province of China

In the present study, the possible evidence of positive selection was analyzed for the neuraminidase (NA) sequences of Guangxi H5N1 strains of China. Based on an overall site-specific positive selection analysis, it was found that NA gene of H5N1 Guangxi strains underwent purifying selection and no significant positively selected sites were identified. For the branch-specific positive selection analysis, there was no positive selection evidence for the branches leading to different poultry hosts (chicken, duck and goose). Conclusively, positive selection seems not possible (if not rare) for the NA gene in influenza H5N1 subtype, at least for the samples found in Guangxi Province of China.

Molecular evolution of the hemagglutinin and neuraminidase genes of pandemic (H1N1) 2009 influenza viruses in Sendai, Japan, during 2009-2011

Analyzing the evolutionary pattern of the influenza A(H1N1)pdm09 strain in different regions is important for understanding its diversification. We therefore conducted this study to elucidate the genetic variability and molecular evolution of the influenza A(H1N1)pdm09 strains that circulated during the 2009-2010 and 2010-2011 influenza seasons in Sendai, Japan. Nasopharyngeal swab specimens were collected from patients with influenza-like illnesses who visited outpatient clinics in Sendai City, Japan, from September 2009 to April 2011. A total of 75 isolates were selected from September 2009 to April 2011 to analyze the genetic changes in the entire hemagglutinin 1 (HA1) segment of the HA gene and the neuraminidase (NA) gene based on sequence analysis. Bayesian coalescent Markov chain Monte Carlo analyses of HA1 and NA gene sequences were performed for further analysis. High sequence identities were observed for HA1 and NA in influenza A(H1N1)pdm09, displaying 99.06 and 99.33 % nucleotide identities, respectively, with the A(H1N1)pdm09 vaccine strain A/California/07/2009. The substitution rates of nucleotides for HA1 in the 2009-2010 and 2010-2011 were 1.5 × 10-3 and 1.6 × 10-3 substitutions per site per year, respectively. Phylogenetic tree analysis demonstrated that Sendai isolates were clustered into global clade 7, which is characterized by an S203T mutation in the HA1 gene. Moreover, two distinct circulation clusters were present in the 2010-2011 season. Mutations were present in antigenic or receptor-binding domains of the HA1 segment, including A141V, S143G, S183P, S185T, and S203T. The Bayesian skyline plot model illustrated a steady rate for the maintenance of genetic diversity, followed by a slight increase in the later part of the 2010-2011 season. Selection analysis revealed that the HA1 (position 197) and NA (position 46) sites were under positive selection; however, no known mutation conferring resistance to NA inhibitors such as H275Y was observed. The effect on control of the influenza A(H1N1)pdm09 virus, including vaccine strain selection, requires continuous monitoring of the strain by genetic surveillance.

Glycan receptor binding of the influenza A virus H7N9 hemagglutinin

The advent of H7N9 in early 2013 is of concern for a number of reasons, including its capability to infect humans, the lack of clarity in the etiology of infection, and because the human population does not have pre-existing immunity to the H7 subtype. Earlier sequence analyses of H7N9 hemagglutinin (HA) point to amino acid changes that predicted human receptor binding and impinge on the antigenic characteristics of the HA. Here, we report that the H7N9 HA shows limited binding to human receptors; however, should a single amino acid mutation occur, this would result in structural changes within the receptor binding site that allow for extensive binding to human receptors present in the upper respiratory tract. Furthermore, a subset of the H7N9 HA sequences demarcating coevolving amino acids appears to be in the antigenic regions of H7, which, in turn, could impact effectiveness of the current WHO-recommended prepandemic H7 vaccines.

Monitoring the antigenic evolution of human influenza A viruses to understand how and when viruses escape from existing immunity

Background

The World Health Organization (WHO) organizes consultations in February and September of each year, spearheaded by an advisory group of experts to analyze influenza surveillance data generated by the WHO Global Influenza Surveillance and Response System (GISRS). The purpose of these consultations is to recommend the composition on influenza virus vaccines for the northern and southern hemispheres, respectively. The latest news of influenza viruses is made available to the public and updated on the WHO website. Although WHO discloses the manner in which it has made the recommendation, usually by considering epidemiological and clinical information to analyze the antigenic and genetic characteristics of seasonal influenza viruses, most individuals do not possess an understanding of antigenic drift and when it occurs.

Findings

We have constructed a web server, named Fluctrl, and implemented a pipeline whereby HA sequence data is downloaded from the Influenza Virus Resource at NCBI along with their isolation information including isolation year and location, which are parsed and managed in MySQL database. By analyzing the frequency of each amino acid residue of the HA1 domain expressed by the viruses on annual basis, users are able to obtain evolutionary dynamics of human influenza viruses corresponding with epidemics. Users are able to upload and analyze their HA1 sequences for generating evolutionary dynamics. In addition, a distribution of amino acid residues at a particular site is represented geographically to trace the location where antigenic variants are seeded.

Conclusions

Fluctrl is constructed for monitoring the antigenic evolution of human influenza A viruses. This tool is intended to inform the general public how and when influenza viruses evade the human body's immunity. Furthermore, leveraging the geographic information, the original locations of emerging influenza viruses can be traced. Fluctrl is freely accessible at http://sb.nhri.org.tw/fluctrl.

Molecular epidemiology and complete genome characterization of H1N1pdm virus from India

Background

Influenza A virus is one of world’s major uncontrolled pathogen, causing seasonal epidemic as well as global pandemic. This was evidenced by recent emergence and continued prevalent 2009 swine origin pandemic H1N1 Influenza A virus, provoking first true pandemic in the past 40 years. In the course of its evolution, the virus acquired many mutations and multiple unidentified molecular determinants are likely responsible for the ability of the 2009 H1N1 virus to cause increased disease severity in humans. Availability of limited data on complete genome hampers the continuous monitoring of this type of events. Outbreaks with considerable morbidity and mortality have been reported from all parts of the country.

Methods/Results

Considering a large number of clinical cases of infection complete genome based sequence characterization of Indian H1N1pdm virus and their phylogenetic analysis with respect to circulating global viruses was undertaken, to reveal the phylodynamic pattern of H1N1pdm virus in India from 2009–2011. The Clade VII was observed as a major circulating clade in phylogenetic analysis. Selection pressure analysis revealed 18 positively selected sites in major surface proteins of H1N1pdm virus.

Conclusions

This study clearly revealed that clade VII has been identified as recent circulating clade in India as well globally. Few clade VII specific well identified markers undergone positive selection during virus evolution. Continuous monitoring of the H1N1pdm virus is warranted to track of the virus evolution and further transmission. This study will serve as a baseline data for future surveillance and also for development of suitable therapeutics.

Reconstruction of the evolutionary dynamics of the A(H1N1)pdm09 influenza virus in Italy during the pandemic and post-pandemic phases

The aim of this study was to reconstruct the evolutionary dynamics of the A(H1N1)pdm09 influenza virus in Italy during two epidemic seasons (2009/2010 and 2010/2011) in the light of the forces driving the evolution of the virus. Nearly six thousands respiratory specimens were collected from patients with influenza-like illness within the framework of the Italian Influenza Surveillance Network, and the A(H1N1)pdm09 hemagglutinin (HA) gene was amplified and directly sequenced from 227 of these. Phylodynamic and phylogeographical analyses were made using a Bayesian Markov Chain Monte Carlo method, and codon-specific positive selection acting on the HA coding sequence was evaluated. The global and local phylogenetic analyses showed that all of the Italian sequences sampled in the post-pandemic (2010/2011) season grouped into at least four highly significant Italian clades, whereas those of the pandemic season (2009/2010) were interspersed with isolates from other countries at the tree root. The time of the most recent common ancestor of the strains circulating in the pandemic season in Italy was estimated to be between the spring and summer of 2009, whereas the Italian clades of the post-pandemic season originated in the spring of 2010 and showed radiation in the summer/autumn of the same year; this was confirmed by a Bayesian skyline plot showing the biphasic growth of the effective number of infections. The local phylogeography analysis showed that the first season of infection originated in Northern Italian localities with high density populations, whereas the second involved less densely populated localities, in line with a gravity-like model of geographical dispersion. Two HA sites, codons 97 and 222, were under positive selection. In conclusion, the A(H1N1)pdm09 virus was introduced into Italy in the spring of 2009 by means of multiple importations. This was followed by repeated founder effects in the post-pandemic period that originated specific Italian clades.

Bayesian coalescent analysis of pandemic H1N1 influenza A virus circulating in the South American region

The first influenza pandemic of this century was declared in April of 2009, with the emergence of a novel H1N1 influenza A virus strain (H1N1pdm). Understanding the evolution of H1N1pdm populations within the South American region is essential for studying global diversification, emergence, resistance and vaccine efficacy. In order to gain insight into these matters, we have performed a Bayesian coalescent Markov Chain Monte Carlo analysis of hemagglutinin (HA) and neuraminidase (NA) gene sequences of all available and comparable HA and NA sequences obtained from H1N1pdm IAV circulating in the South American region. High evolutionary rates and fast population growths characterize the population dynamics of H1N1pdm strains in this region of the world. A significant contribution of first codon position to the mean evolutionary rate was found for both genes studied, revealing a high contribution of non-synonymous substitutions to the mean substitution rate. In the 178days period covered by these studies, substitutions in all HA epitope regions can be observed. HA substitutions D239G/N and Q310H have been observed only in Brazilian patients. While substitution D239G/N is not particularly associated to a specific genetic lineage, all strains bearing substitution Q310H were assigned to clade 6, suggesting a founder effect. None of the substitutions found in the NA proteins of H1N1pdm strains isolated in South America appears sufficiently close to affect the drug binding pocket for the three NA inhibitor antivirals tested. A more detailed analysis of NA proteins revealed epitope differences among 2010 vaccine and H1N1pdm IAV strains circulating in the South American region.

Antigenic sites of H1N1 influenza virus hemagglutinin revealed by natural isolates and inhibition assays

The antigenic sites of hemagglutinin (HA) are crucial for understanding antigenic drift and vaccine strain selection for influenza viruses. In 1982, 32 epitope residues (called laboratory epitope residues) were proposed for antigenic sites of H1N1 HA based on the monoclonal antibody-selected variants. Interestingly, these laboratory epitope residues only cover 28% (23/83) mutation positions for 9 H1N1 vaccine strain comparisons (from 1977 to 2009). Here, we propose the entropy and likelihood ratio to model amino acid diversity and antigenic variant score for inferring 41 H1N1 HA epitope residues (called natural epitope residues) with statistically significant scores according to 1572 HA sequences and 197 pairs of HA sequences with hemagglutination inhibition (HI) assays of natural isolates. By combining both natural and laboratory epitope residues, we identified 62 (11 overlapped) residues clustered into five antigenic sites (i.e., A-E) which are highly correlated to the antigenic sites of H3N2 HA. Our method recognizes sites A, B and C as critical sites for escaping from neutralizing antibodies in H1N1 virus. Experimental results show that the accuracies of our models are 81.2% and 82.2% using 41 and 62 epitope residues, respectively, for predicting antigenic variants on 197 paring HA sequences. In addition, our model can detect the emergence of epidemic strains and reflect the genetic diversity and antigenic variant between the vaccine and circulating strains. Finally, our model is theoretically consistent with the evolution rates of H3N2 and H1N1 viruses and is often consistent to WHO vaccine strain selections. We believe that our models and the inferred antigenic sites of HA are useful for understanding the antigenic drift and evolution of influenza A H1N1 virus.

Subtype- and antigenic site-specific differences in biophysical influences on evolution of influenza virus hemagglutinin

Background

Influenza virus undergoes rapid evolution by both antigenic shift and antigenic drift. Antibodies, particularly those binding near the receptor-binding site of hemagglutinin (HA) or the neuraminidase (NA) active site, are thought to be the primary defense against influenza infection, and mutations in antibody binding sites can reduce or eliminate antibody binding. The binding of antibodies to their cognate antigens is governed by such biophysical properties of the interacting surfaces as shape, non-polar and polar surface area, and charge.

Methods

To understand forces shaping evolution of influenza virus, we have examined HA sequences of human influenza A and B viruses, assigning each amino acid values reflecting total accessible surface area, non-polar and polar surface area, and net charge due to the side chain. Changes in each of these values between neighboring sequences were calculated for each residue and mapped onto the crystal structures.

Results

Areas of HA showing the highest frequency of pairwise changes agreed well with previously identified antigenic sites in H3 and H1 HAs, and allowed us to propose more detailed antigenic maps and novel antigenic sites for H1 and influenza B HA. Changes in biophysical properties differed between HAs of different subtypes, and between different antigenic sites of the same HA. For H1, statistically significant differences in several biophysical quantities compared to residues lying outside antigenic sites were seen for some antigenic sites but not others. Influenza B antigenic sites all show statistically significant differences in biophysical quantities for all antigenic sites, whereas no statistically significant differences in biophysical quantities were seen for any antigenic site is seen for H3. In many cases, residues previously shown to be under positive selection at the genetic level also undergo rapid change in biophysical properties.

Conclusions

The biophysical consequences of amino acid changes introduced by antigenic drift vary from subtype to subtype, and between different antigenic sites. This suggests that the significance of antibody binding in selecting new variants may also be variable for different antigenic sites and influenza subtypes.

Palindromes drive the re-assortment in Influenza A

Different subtypes of Influenza A virus are associated with species specific, zoonotic or pandemic Influenza. The cause of its severity underlies in complicated evolution of its segmented RNA genome. Although genetic shift and genetic drift are well known in the evolution of this virus, we reported the significant role of unique RNA palindromes in its evolution. Our computational approach identified the existence of unique palindromes in each subtype of Influenza A virus with its absence in Influenza B relating the fact of virulence and vigorous genetic hitchhiking in Influenza A. The current study focused on the re-assortment event responsible for the emergence of pandemic-2009 H1N1 virus, which is associated with outgrow of new palindrome and in turn, changing its RNA structure. We hypothesize that the change in RNA structure due to the presence of palindrome facilitates the event of re-assortment in Influenza A. Thus the evolutionary process of Influenza A is much more complicated as previously known, and that has been demonstrated in this study.

Genetic structure of human A/H1N1 and A/H3N2 influenza virus on Corsica Island: phylogenetic analysis and vaccine strain match, 2006-2010

BACKGROUND

The aim of this study was to analyse the genetic patterns of Hemagglutinin (HA) genes of influenza A strains circulating on Corsica Island during the 2006-2009 epidemic seasons and the 2009-2010 pandemic season.

METHODS

Nasopharyngeal samples from 371 patients with influenza-like illness (ILI) were collected by General Practitioners (GPs) of the Sentinelles Network through a randomised selection routine.

RESULTS

Phylogenetic analysis of HA revealed that A/H3N2 strains circulating on Corsica were closely related to the WHO recommended vaccine strains in each analyzed season (2006-2007 to 2008-2009). Seasonal Corsican influenza A/H1N1 isolated during the 2007-2008 season had drifted towards the A/Brisbane/59/2007 lineage, the A/H1N1 vaccine strain for the 2008-2009 season. The A/H1N1 2009 (A/H1N1pdm) strains isolated on Corsica Island were characterized by the S220T mutation specific to clade 7 isolates. It should be noted that Corsican isolates formed a separate sub-clade of clade 7 as a consequence of the presence of the fixed substitution D222E. The percentages of the perfect match vaccine efficacy, estimated by using the p(epitope) model, against influenza viruses circulating on Corsica Island varied substantially across the four seasons analyzed, and tend to be highest for A/H1N1 compared with A/H3N2 vaccines, suggesting that cross-immunity seems to be stronger for the H1 HA gene.

CONCLUSION

The molecular analysis of the HA gene of influenza viruses that circulated on Corsica Island between 2006-2010 showed for each season the presence of a dominant lineage characterized by at least one fixed mutation. The A/H3N2 and A/H1N1pdm isolates were characterized by multiples fixation at antigenic sites. The fixation of specific mutations at each outbreak could be explained by the combination of a neutral phenomenon and a founder effect, favoring the presence of a dominant lineage in a closed environment such as Corsica Island.