The role of genomics in tracking the evolution of influenza A virus

Wastewater-based epidemiology of influenza viruses: a systematic review

INTRODUCTION

Wastewater-based epidemiology (WBE) has emerged as a valuable public health tool for monitoring the circulation of many pathogens, including influenza viruses (IVs). The general aim of this study is to systematically retrieve and summarize evidence on the use of WBE for supporting influenza surveillance. Specific objectives are: (i) to map influenza monitoring activities using WBE; (ii) to assess the performance of viral recovery methods; (iii) to explore association with clinical data; (iv) to evaluate the feasibility of typing/subtyping IVs directly from wastewater.

METHODS

We conducted a systematic review following the PRISMA guidelines, focusing on original data from peer-reviewed studies identified through PubMed/Medline, Scopus, and Web of Science.

RESULTS

Of 882 identified citations, 42 studies were included in the review. IVs detection was reported in all but one study, although typically at lower concentration than SARS-CoV-2. Thirteen studies (38.09 %) performed comparative analysis of different protocols, with mostly inconclusive results. Detection of IVs in the solid fraction of wastewater samples generally outperformed detection in the supernatant/liquid. Additionally, we describe the findings from 22 studies (52.38 %) that examined the link between environmental viral concentrations and clinical data, and 14 studies (33.33 %) that described IVs subtyping in wastewater.

CONCLUSION

WBE has the potential to monitor influenza circulation in humans and animals, offering insights into outbreak size and circulating IVs subtypes. However, several key areas remain unexplored. Further research is needed to refine experimental techniques and standardize protocols, and to understand how to successfully integrate WBE data into public health strategies for influenza control.

VirusWarn: A mutation-based early warning system to prioritize concerning SARS-CoV-2 and influenza virus variants from sequencing data

The rapid evolution of respiratory viruses is characterized by the emergence of variants with concerning phenotypes that are efficient in antibody escape or show high transmissibility. This necessitates timely identification of such variants by surveillance networks to assist public health interventions. Here, we introduce VirusWarn, a comprehensive system designed for detecting, prioritizing, and warning of emerging virus variants from large genomic datasets. VirusWarn uses both manually-curated rules and machine-learning (ML) classifiers to generate and rank pathogen sequences based on mutations of concern and regions of interest. Validation results for SARS-CoV-2 showed that VirusWarn successfully identifies variants of concern in both assessments, with manual- and ML-derived criteria from positive selection analyses. Although initially developed for SARS-CoV-2, VirusWarn was adapted to Influenza viruses and their dynamics, and provides a robust performance, integrating a scheme that accounts for fixed mutations from past seasons. HTML reports provide detailed results with searchable tables and visualizations, including mutation plots and heatmaps. Because VirusWarn is written in Nextflow, it can be easily adapted to other pathogens, demonstrating its flexibility and scalability for genomic surveillance efforts.

Genomic diversity and evolutionary dynamics of Influenza A viruses in Colombian swine: implications for one health surveillance and control

Influenza A viruses (IAV) impose significant respiratory disease burdens in both swine and humans worldwide, with frequent human-to-swine transmission driving viral evolution in pigs and highlighting the risk at the animal-human interface. Therefore, a comprehensive One Health approach (interconnection among human, animal, and environmental health) is needed for IAV prevention, control, and response. Animal influenza genomic surveillance remains limited in many Latin American countries, including Colombia. To address this gap, we genetically characterized 170 swine specimens from Colombia (2011-2017). Whole genome sequencing revealed a predominance of pandemic-like H1N1 lineage, with a minority belonging to H3N2 and H1N2 human seasonal-like lineage and H1N1 early classical swine lineages. Significantly, we have identified reassortant and recombinant viruses (H3N2, H1N1) not previously reported in Colombia. This suggests a broad genotypic viral diversity, likely resulting from reassortment between classical endemic viruses and new introductions established in Colombia's swine population (e.g. the 2009 H1N1 pandemic). Our study highlights the importance of a One Health approach in disease control, particularly in an ecosystem where humans are a main source of IAV to swine populations, and emphasizes the need for continued surveillance and enhanced biosecurity measures. The co-circulation of multiple subtypes in regions with high swine density facilitates viral exchange, underscoring the importance of monitoring viral evolution to inform vaccine selection and public health policies locally and globally.

Whole genome sequence analysis of shallot virus X from India reveals it to be a natural recombinant with positive selection pressure

BACKGROUND

Shallots are infected by various viruses like Onion yellow dwarf virus (OYDV), Leek yellow stripe virus (LYSV), Shallot latent virus (SLV) and Shallot virus X (ShVX). In India, they have been found to be persistently infected by ShVX. ShVX also infects onion and garlic in combination with other carlaviruses and potyviruses. ShVX is a member of genus Allexivirus of family Alphaflexiviridae. ShVX has a monopartite genome, which is represented by positive sense single-stranded RNA. Globally, only six complete and 3 nearly complete genome sequences of ShV X are reported to date. This number is insufficient to measure a taxon's true molecular diversity. Moreover, the complete genome sequence of ShVX from Asia has not been reported as yet. Therefore, this study was undertaken to generate a complete genome sequence of ShVX from India.

RESULTS

Shallot virus X (ShVX) is one of the significant threats to Allium crop production. In this study, we report the first complete genome sequence of the ShVX from India through Next-generation sequencing (NGS). The complete genome of the ShVX (Accession No. OK104171), from this study comprised 8911 nucleotides. In-silico analysis of the sequence revealed variability between this isolate and isolates from other countries. The dissimilarities are spread all over the genome specifically some non-coding intergenic regions. Statistical analysis of individual genes for site-specific selection indicates a positive selection in NABP region. The presence of a recombination event was detected in coat protein region. The sequence similarity percentage and phylogenetic analysis indicate ShVX Indian isolate is a distinctly different isolate. Recombination and site-specific selection may have a function in the evolution of this isolate. This is the first detailed study of the ShVX complete genome sequence from Southeast Asia.

CONCLUSION

This study presents the first report of the entire genome sequence of an Indian isolate of ShVX along with an in-depth exploration of its evolutionary traits. The findings highlight the Indian variant as a naturally occurring recombinant, emphasizing the substantial role of recombination in the evolution of this viral species. This insight into the molecular diversity of strains within a specific geographical region holds immense significance for comprehending and forecasting potential epidemics. Consequently, the insights garnered from this research hold practical value for shaping ShVX management strategies and providing a foundation for forthcoming studies delving into its evolutionary trajectory.

Prediction of Antigenic Distance in Influenza A Using Attribute Network Embedding

Owing to the rapid changes in the antigenicity of influenza viruses, it is difficult for humans to obtain lasting immunity through antiviral therapy. Hence, tracking the dynamic changes in the antigenicity of influenza viruses can provide a basis for vaccines and drug treatments to cope with the spread of influenza viruses. In this paper, we developed a novel quantitative prediction method to predict the antigenic distance between virus strains using attribute network embedding techniques. An antigenic network is built to model and combine the genetic and antigenic characteristics of the influenza A virus H3N2, using the continuous distributed representation of the virus strain protein sequence (ProtVec) as a node attribute and the antigenic distance between virus strains as an edge weight. The results show a strong positive correlation between supplementing genetic features and antigenic distance prediction accuracy. Further analysis indicates that our prediction model can comprehensively and accurately track the differences in antigenic distances between vaccines and influenza virus strains, and it outperforms existing methods in predicting antigenic distances between strains.

Retrospective immunogenicity analysis of seasonal flu H3N2 vaccines recommended in the past ten years using immunized animal sera

BACKGROUND

Influenza A (H3N2) virus (A/H3N2) has complex antigenic evolution, resulting in frequent vaccine strain updates. We aimed to evaluate the protective effect of the vaccine strains on the circulating strains from past ten years and provide a basis for finding a broader and more efficient A/H3N2 vaccine strain.

METHODS

Eighty-four representative circulating A/H3N2 strains were selected from 65,791 deposited sequences in 2011-2020 and pseudotyped viruses were constructed with the VSV vector. We immunized guinea pigs with DNA vaccine containing the A/H3N2 components of the vaccine strains from 2011 to 2021 and tested neutralizing antibody against the pseudotyped viruses. We used a hierarchical clustering method to classify the eighty-four representative strains into different antigenic clusters. We also immunized animals with monovalent vaccine stock of the vaccine strains for the 2020-2021 and 2021-2022 seasons and tested neutralizing antibody against the pseudotyped viruses.

FINDINGS

The vaccine strains PE/09, VI/11 and TE/12 induced higher levels of neutralizing antibody against representative strains circulating in recommended year and the year immediately prior whereas vaccine strains HK/14, HK/19 and CA/20 induced poor neutralization against all representative strains. The representative strains were divided into five antigenic clusters (AgV), which were not identical to gene clades. The AgV5 strains were most difficult to be protected among the five clusters. Compared with single-dose immunization, three doses of monovalent vaccine stock (HK/19 or CA/20) could induce stronger and broader neutralizing antibodies against strains in each of the antigenic clusters.

INTERPRETATION

The protective effect of vaccine strains indicated that the accurate selection of A/H3N2 vaccine strains must remain a top priority. By increasing the frequency of immunization, stronger and broader neutralizing antibodies against strains in all antigenic clusters were induced, which provides direction for a new immunization strategy.

FUNDING

This work was supported by a grant from National Key R&D Program of China (No. 2021YFC2301700).

Influenza A Virus Exacerbates Group A Streptococcus Infection and Thwarts Anti-Bacterial Inflammatory Responses in Murine Macrophages

Seasonal influenza epidemics pose a considerable hazard for global health. In the past decades, accumulating evidence revealed that influenza A virus (IAV) renders the host vulnerable to bacterial superinfections which in turn are a major cause for morbidity and mortality. However, whether the impact of influenza on anti-bacterial innate immunity is restricted to the vicinity of the lung or systemically extends to remote sites is underexplored. We therefore sought to investigate intranasal infection of adult C57BL/6J mice with IAV H1N1 in combination with bacteremia elicited by intravenous application of Group A Streptococcus (GAS). Co-infection in vivo was supplemented in vitro by challenging murine bone marrow derived macrophages and exploring gene expression and cytokine secretion. Our results show that viral infection of mice caused mild disease and induced the depletion of CCL2 in the periphery. Influenza preceding GAS infection promoted the occurrence of paw edemas and was accompanied by exacerbated disease scores. In vitro co-infection of macrophages led to significantly elevated expression of TLR2 and CD80 compared to bacterial mono-infection, whereas CD163 and CD206 were downregulated. The GAS-inducible upregulation of inflammatory genes, such as Nos2, as well as the secretion of TNFα and IL-1β were notably reduced or even abrogated following co-infection. Our results indicate that IAV primes an innate immune layout that is inadequately equipped for bacterial clearance.

Characterization of SARS-CoV-2 ORF6 deletion variants detected in a nosocomial cluster during routine genomic surveillance, Lyon, France

During routine molecular surveillance of SARS-CoV-2 performed at the National Reference Center of Respiratory Viruses (Lyon, France) (n = 229 sequences collected February-April 2020), two frameshifting deletions were detected in the open reading frame 6, at the same position (27267). While a 26-nucleotide deletion variant (D26) was only found in one nasopharyngeal sample in March 2020, the 34-nucleotide deletion (D34) was found within a single geriatric hospital unit in 5/9 patients and one health care worker in April 2020. Phylogeny analysis strongly suggested a nosocomial transmission of D34, with potential fecal transmission, as also identified in a stool sample. No difference in disease severity was observed between patients hospitalized in the geriatric unit infected with WT or D34. In vitro D26 and D34 characterization revealed comparable replication kinetics with the wild-type (WT), but differential host immune responses. While interferon-stimulated genes were similarly upregulated after infection with WT and ORF6 variants, the latter specifically induced overexpression of 9 genes coding for inflammatory cytokines in the NF-kB pathway, including CCL2/MCP1, PTX3, and TNFα, for which high plasma levels have been associated with severe COVID-19. Our findings emphasize the need to monitor the occurrence of ORF6 deletions and assess their impact on the host immune response.

Reassortment with dominant chicken H9N2 influenza virus contributed to the fifth H7N9 virus human epidemic

H9N2 Avian influenza virus (AIV) is regarded as a principal donor of viral genes through reassortment to co-circulating influenza viruses that can result in zoonotic reassortants. Whether H9N2 virus can maintain sustained evolutionary impact on such reassortants is unclear. Since 2013, avian H7N9 virus had caused five sequential human epidemics in China; the fifth wave in 2016-2017 was by far the largest but the mechanistic explanation behind the scale of infection is not clear. Here, we found that, just prior to the fifth H7N9 virus epidemic, H9N2 viruses had phylogenetically mutated into new sub-clades, changed antigenicity and increased its prevalence in chickens vaccinated with existing H9N2 vaccines. In turn, the new H9N2 virus sub-clades of PB2 and PA genes, housing mammalian adaptive mutations, were reassorted into co-circulating H7N9 virus to create a novel dominant H7N9 virus genotype that was responsible for the fifth H7N9 virus epidemic. H9N2-derived PB2 and PA genes in H7N9 virus conferred enhanced polymerase activity in human cells at 33°C and 37°C, and increased viral replication in the upper and lower respiratory tracts of infected mice which could account for the sharp increase in human cases of H7N9 virus infection in the 2016-2017 epidemic. The role of H9N2 virus in the continual mutation of H7N9 virus highlights the public health significance of H9N2 virus in the generation of variant reassortants of increasing zoonotic potential.IMPORTANCEAvian H9N2 influenza virus, although primarily restricted to chicken populations, is a major threat to human public health by acting as a donor of variant viral genes through reassortment to co-circulating influenza viruses. We established that the high prevalence of evolving H9N2 virus in vaccinated flocks played a key role, as donor of new sub-clade PB2 and PA genes in the generation of a dominant H7N9 virus genotype (G72) with enhanced infectivity in humans during the 2016-2017 N7N9 virus epidemic. Our findings emphasize that the ongoing evolution of prevalent H9N2 virus in chickens is an important source, via reassortment, of mammalian adaptive genes for other influenza virus subtypes. Thus, close monitoring of prevalence and variants of H9N2 virus in chicken flocks is necessary in the detection of zoonotic mutations.

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.

Synthetic biology for bioengineering virus-like particle vaccines

Vaccination is the most effective method of disease prevention and control. Many viruses and bacteria that once caused catastrophic pandemics (e.g., smallpox, poliomyelitis, measles, and diphtheria) are either eradicated or effectively controlled through routine vaccination programs. Nonetheless, vaccine manufacturing remains incredibly challenging. Viruses exhibiting high antigenic diversity and high mutation rates cannot be fairly contested using traditional vaccine production methods and complexities surrounding the manufacturing processes, which impose significant limitations. Virus-like particles (VLPs) are recombinantly produced viral structures that exhibit immunoprotective traits of native viruses but are noninfectious. Several VLPs that compositionally match a given natural virus have been developed and licensed as vaccines. Expansively, a plethora of studies now confirms that VLPs can be designed to safely present heterologous antigens from a variety of pathogens unrelated to the chosen carrier VLPs. Owing to this design versatility, VLPs offer technological opportunities to modernize vaccine supply and disease response through rational bioengineering. These opportunities are greatly enhanced with the application of synthetic biology, the redesign and construction of novel biological entities. This review outlines how synthetic biology is currently applied to engineer VLP functions and manufacturing process. Current and developing technologies for the identification of novel target-specific antigens and their usefulness for rational engineering of VLP functions (e.g., presentation of structurally diverse antigens, enhanced antigen immunogenicity, and improved vaccine stability) are described. When applied to manufacturing processes, synthetic biology approaches can also overcome specific challenges in VLP vaccine production. Finally, we address several challenges and benefits associated with the translation of VLP vaccine development into the industry.

Influenza virus RNA polymerase may be activated inside the virion

Influenza A and B virions are packaged with their polymerases to catalyse RNA-dependent RNA polymerase activity. Since there is no evidence to rule in or out the permissiveness of influenza virions to triphosphate ribonucleotides, we functionally evaluated this. We found the means to stimulate influenza A and B RNA polymerase activity inside the virion, called natural endogenous RNA polymerase (NERP) activity. Stimulation of NERP activity increased up to 3 log10 viral RNA content, allowing the detection of influenza virus in otherwise undetectable clinical samples. NERP activation also improved our capacity to sequence misidentified regions of the influenza genome from clinical samples. By treating the samples with the ribavirin triphosphate we inhibited NERP activity, which confirms our hypothesis and highlights that this assay could be used to screen antiviral drugs. Altogether, our data show that NERP activity could be explored to increase molecular diagnostic sensitivity and/or to develop antiviral screening assays.

Variant antigen repertoires in Trypanosoma congolense populations and experimental infections can be profiled from deep sequence data using universal protein motifs

African trypanosomes are vector-borne hemoparasites of humans and animals. In the mammal, parasites evade the immune response through antigenic variation. Periodic switching of the variant surface glycoprotein (VSG) coat covering their cell surface allows sequential expansion of serologically distinct parasite clones. Trypanosome genomes contain many hundreds of VSG genes, subject to rapid changes in nucleotide sequence, copy number, and chromosomal position. Thus, analyzing, or even quantifying, VSG diversity over space and time presents an enormous challenge to conventional techniques. Indeed, previous population genomic studies have overlooked this vital aspect of pathogen biology for lack of analytical tools. Here we present a method for analyzing population-scale VSG diversity in Trypanosoma congolense from deep sequencing data. Previously, we suggested that T. congolense VSGs segregate into defined “phylotypes” that do not recombine. In our data set comprising 41 T. congolense genome sequences from across Africa, these phylotypes are universal and exhaustive. Screening sequence contigs with diagnostic protein motifs accurately quantifies relative phylotype frequencies, providing a metric of VSG diversity, called the “variant antigen profile.” We applied our metric to VSG expression in the tsetse fly, showing that certain, rare VSG phylotypes may be preferentially expressed in infective, metacyclic-stage parasites. Hence, variant antigen profiling accurately and rapidly determines the T. congolense VSG gene and transcript repertoire from sequence data, without need for manual curation or highly contiguous sequences. It offers a tractable approach to measuring VSG diversity across strains and during infections, which is imperative to understanding the host–parasite interaction at population and individual scales.

Evolutionary Pattern and Large-Scale Architecture of Mutation Networks of 2009 A (H1N1) Influenza A Virus

The adaptive evolution of influenza virus is an important question, but predicting its evolutionary future will be more challenging. Here, we investigated the mutation characteristic of influenza virus based on the complete genome data of 2009 (H1N1) influenza A virus. By assuming that evolution proceeds via the accumulation of mutations, we analyzed the mutation networks at four different time stages and found that the network structure follows the characteristics of a scale-free network. These results will be important for epidemiology and the future control of influenza viruses. Furthermore, we predicted the predominant mutation virus strain by using the early mutation network of influenza viruses, and this result was consistent with the WHO recommendation for the candidate vaccine of influenza virus. The key contribution of this study is that we explained the biological significance of this scale-free network for influenza pandemic and provided a potential method for predicting the candidate vaccine by using the early-stage network.

Sweep Dynamics (SD) plots: Computational identification of selective sweeps to monitor the adaptation of influenza A viruses

Monitoring changes in influenza A virus genomes is crucial to understand its rapid evolution and adaptation to changing conditions e.g. establishment within novel host species. Selective sweeps represent a rapid mode of adaptation and are typically observed in human influenza A viruses. We describe Sweep Dynamics (SD) plots, a computational method combining phylogenetic algorithms with statistical techniques to characterize the molecular adaptation of rapidly evolving viruses from longitudinal sequence data. SD plots facilitate the identification of selective sweeps, the time periods in which these occurred and associated changes providing a selective advantage to the virus. We studied the past genome-wide adaptation of the 2009 pandemic H1N1 influenza A (pH1N1) and seasonal H3N2 influenza A (sH3N2) viruses. The pH1N1 influenza virus showed simultaneous amino acid changes in various proteins, particularly in seasons of high pH1N1 activity. Partially, these changes resulted in functional alterations facilitating sustained human-to-human transmission. In the evolution of sH3N2 influenza viruses, we detected changes characterizing vaccine strains, which were occasionally revealed in selective sweeps one season prior to the WHO recommendation. Taken together, SD plots allow monitoring and characterizing the adaptive evolution of influenza A viruses by identifying selective sweeps and their associated signatures.

Pathogen genomic surveillance elucidates the origins, transmission and evolution of emerging viral agents in China

In the past twenty years, numerous novel zoonotic viral agents with pandemic potential have emerged in China, such as the severe acute respiratory syndrome (SARS) coronavirus and, more recently, the avian-origin influenza A/H7N9 virus, which have caused outbreaks among humans with high morbidity and mortality. In addition, several emerging and re-emerging viral pathogens have also been imported into China from travelers, e.g. the Middle East respiratory syndrome (MERS) coronavirus and Zika virus (ZIKV). Herein, we review these emerging viral pathogens in China and focus on how surveillance by pathogen genomics has been employed to discover and annotate novel pathogenic agents, identify natural reservoirs, monitor the transmission events and delineate their evolution and adaption to the human host. We also highlight the application of genomic sequencing in the recent Ebola epidemics in Western Africa. In summary, genomic sequencing has become a standard research tool in the field of emerging infectious diseases which has been proven invaluable in containing these viral infections and reducing burden of disease in humans and animals. Genomic surveillance of pathogenic agents will serve as a key epidemiological and research tool in the modern era of precision infectious diseases and in the future studies of virosphere.

Genetic variation and co-evolutionary relationship of RNA polymerase complex segments in influenza A viruses

The RNA polymerase complex (RNApc) in influenza A viruses (IVs) is composed of the PB2, PB1 and PA subunits, which are encoded by the three longest genome segments (Seg1-3) and are responsible for the replication of vRNAs and transcription of viral mRNAs. However, the co-evolutionary relationships of the three segments from the known 126 subtypes IVs are unclear. In this study, we performed a detailed analysis based on a total number of 121,191 nucleotide sequences. Three segment sequences were aligned before the repeated, incomplete and mixed sequences were removed for homologous and phylogenetic analyses. Subsequently, the estimated substitution rates and TMRCAs (Times for Most Recent Common Ancestor) were calculated by 175 representative IVs. Tracing the cladistic distribution of three segments from these IVs, co-evolutionary patterns and trajectories could be inferred. The further correlation analysis of six internal protein coding segments reflect the RNApc segments have the closer correlation than others during continuous reassortments. This global approach facilitates the establishment of a fast antiviral strategy and monitoring of viral variation.

In Silico Vaccine Strain Prediction for Human Influenza Viruses

Vaccines preventing seasonal influenza infections save many lives every year; however, due to rapid viral evolution, they have to be updated frequently to remain effective. To identify appropriate vaccine strains, the World Health Organization (WHO) operates a global program that continually generates and interprets surveillance data. Over the past decade, sophisticated computational techniques, drawing from multiple theoretical disciplines, have been developed that predict viral lineages rising to predominance, assess their suitability as vaccine strains, link genetic to antigenic alterations, as well as integrate and visualize genetic, epidemiological, structural, and antigenic data. These could form the basis of an objective and reproducible vaccine strain-selection procedure utilizing the complex, large-scale data types from surveillance. To this end, computational techniques should already be incorporated into the vaccine-selection process in an independent, parallel track, and their performance continuously evaluated.

Public health concerns of highly pathogenic avian influenza H5N1 endemicity in Africa

Highly pathogenic avian influenza virus (HPAIV) H5N1 was first officially reported in Africa in 2006; thereafter this virus has spread rapidly from Nigeria to 11 other African countries. This study was aimed at utilizing data from confirmed laboratory reports to carry out a qualitative evaluation of the factors responsible for HPAI H5N1 persistence in Africa and the public health implications; and to suggest appropriate control measures. Relevant publications were sought from data banks and repositories of FAO, OIE, WHO, and Google scholars. Substantiated data on HPAI H5N1 outbreaks in poultry in Africa and in humans across the world were mined. HPAI H5N1 affects poultry and human populations, with Egypt having highest human cases (346) globally. Nigeria had a reinfection from 2014 to 2015, with outbreaks in Côte d'Ivoire, Ghana, Niger, Nigeria, and Burkina Faso throughout 2016 unabated. The persistence of this virus in Africa is attributed to the survivability of HPAIV, ability to evolve other subtypes through genetic reassortment, poor biosecurity compliance at the live bird markets and poultry farms, husbandry methods and multispecies livestock farming, poultry vaccinations, and continuous shedding of HPAIV, transboundary transmission of HPAIV through poultry trades; and transcontinental migratory birds. There is, therefore, the need for African nations to realistically reassess their status, through regular surveillance and be transparent with HPAI H5N1 outbreak data. Also, it is important to have an understanding of HPAIV migration dynamics which will be helpful in epidemiological modeling, disease prevention, control and eradication measures.

N-glycan profiles in H9N2 avian influenza viruses from chicken eggs and human embryonic lung fibroblast cells

N-glycosylation can affect the host specificity, virulence and infectivity of influenza A viruses (IAVs). In this study, the distribution and evolution of N-glycosylation sites in the hemagglutinin (HA) and neuraminidase (NA) of H9N2 virus were explored using phylogenetic analysis. Then, one strain of the H9N2 subtypes was proliferated in the embryonated chicken eggs (ECE) and human embryonic lung fibroblast cells (MRC-5) system. The proliferated viral N-glycan profiles were analyzed by a glycomic method that combined the lectin microarray and MALDI-TOF/TOF-MS. As a result, HA and NA of H9N2 viruses prossess six and five highly conserved N-glycosylation sites, respectively. Sixteen lectins (e.g., MAL-II, SNA and UEA-I) had increased expression levels of the glycan structures in the MRC-5 compared with the ECE system; however, 6 lectins (e.g., PHA-E, PSA and DSA) had contrasting results. Eleven glycans from the ECE system and 13 glycans from the MRC-5 system were identified. Our results showed that the Fucα-1,6GlcNAc(core fucose) structure was increased, and pentaantennary N-glycans were only observed in the ECE system. The SAα2-3/6Gal structures were highly expressed and Fucα1-2Galβ1-4GlcNAc structures were only observed in the MRC-5 system. We conclude that the existing SAα2-3/6Gal sialoglycans make the offspring of the H9N2 virus prefer entially attach to each other, which decreases the virulence. Alterations in the glycosylation sites for the evolution and role of IAVs have been widely described; however, little is known about the exact glycan structures for the same influenza strain from different hosts. Our findings may provide a novel way for further discussing the molecular mechanism of the viral transmission and virulence associated with viral glycosylation in avian and human hosts as well as vital information for designing a vaccine against influenza and other human viruses.

Semi-quantitative Influenza A population averages from a multiplex respiratory viral panel (RVP): potential for reflecting target sequence changes affecting the assay

Background

Yearly influenza virus mutations potentially affect the performance of molecular assays, if nucleic acid changes involve the sequences in the assay. Because individual patient viral loads depend on variables such as duration of illness, specimen type, age, and immunosuppression, we examined seasonal population averages of positive tests to smooth inherent variability.

Methods

We studied the population seasonal averages of the semi-quantitative nAMPs for the influenza matrix and hemagglutinin genes in the GenMark (Carlsbad, CA) Respiratory Viral Panel assay between 3 institutions over 3 Influenza seasons.

Results

Population average nAMPs were strikingly consistent between separate institutions, but differed substantially between H3N2 and H1N1 seasons. In the 2012–2013 and 2014–2015 influenza seasons, matrix gene H3N2 nAMP averages were 50–70% less than those of the same assay in the 2013–2014 H1N1 season. Influenza strains representative of these seasons were grown in tissue culture and when the supernatant virus was adjusted to the same copy number using a TaqMan assay, the same relative differences were reproduced in the RVP assay. Because the sequences for the PCR and PCR product detection in the GenMark assay are proprietary, the manufacturer provided single stranded DNA matching the capture probe for the representative H3N2 (3 mismatches) and H1N1 strains (2 different mismatches). Equimolar concentrations of these synthetic DNA sequences gave average nAMP values that closely correlated with the average nAMPS of the representative strains and their respective seasonal averages.

Conclusions

Seasonal averages of semi-quantitative data may provide a means to follow assay performance as a reflection of the effects of molecular drift.

Asymmetric structure of the influenza A virus and novel function of the matrix protein M1

Influenza virus is an enveloped virus. It comprises two major modules: external lipoprotein envelope and internal ribonucleoprotein (RNP) containing the genomic negative-strand RNA. Lipoprotein envelope contains four vital proteins: hemagglutinin (HA), neuraminidase (NA), transmembrane ionic channel M2, and minor amounts of nuclear export protein NEP. RNP contains RNA and four polypeptides: major nucleocapsid protein NP and three polymerase subunits PB1, PB2, PA. Both modules are linked with each other by matrix M1 maintaining the virus integrity. According to the structural function, NP and M1 are predominant in virus particle in the amounts of 1000 and 3000 molecules, respectively. In addition to the structural function, M1 plays a role in regulation of intracellular and nuclear migration of viral RNP and virus assembly, referred as budding process, at the plasma membrane in infected cells. The bipolar structure of the influenza virus characterized by asymmetric location of RNP and nonregular distribution of M1 and M2 inside the virion is reviewed. The role of M1 in maintaining the asymmetric structure of the virus particle and regulation of RNP transport inside virus particle is considered. First experimental data confirming (i) intravirion RNP transport and its outside exit directed by the M1 and (ii) the importance of this process in virus uncoating and initiation of infection in target cell are discussed. A novel class of antiviral agents activating ATP-ase of the early endosome compartment in the target cell is discussed.

Re-emergence of H3N2 strains carrying potential neutralizing mutations at the N-linked glycosylation site at the hemagglutinin head, post the 2009 H1N1 pandemic

Background

Seasonally prevalent H1N1 and H3N2 influenza A viruses have evolved by antigenic drift; this evolution has resulted in the acquisition of asparagine (N)-linked glycosylation sites (NGSs) in the globular head of hemagglutinin (HA), thereby affecting the antigenic and receptor-binding properties, as well as virulence. An epidemiological survey indicated that although the traditional seasonal H1N1 strain had disappeared, H3N2 became predominant again in the seasons (2010–11 and 2011–12) immediately following the H1N1 pandemic of 2009. Interestingly, although the 2009 pandemic H1N1 strain (H1N1pdm09) lacks additional NGSs, clinically isolated H3N2 strains obtained during these seasons gained N (Asn) residues at positions 45 and 144 of HA that forms additional NGSs.

Methods

To investigate whether these NGSs are associated with re-emergence of H3N2 within the subtype, we tested the effect of amino acid substitutions on neutralizing activity by using the antisera raised against H3N2 strains with or without additional NGSs. Furthermore, because the N residue at position 144 of HA was identified as the site of mismatch between the vaccine and epidemic strains of 2011–2012, we generated mutant viruses by reverse genetics and tested the functional importance of this particular NGS for antibody-mediated neutralization by intranasal inoculation of mice.

Results

The results indicated that amino acid substitution at residue 144 significantly affected neutralization activity, acting as an escape mutation.

Conclusions

Our data suggest that the newly acquired NGSs in the HA globular head may play an important role in the re-emergence of endemic seasonal H3N2 strain by aiding the escape from humoral immunity.

Gene overlapping and size constraints in the viral world

Background

Viruses are the simplest replicating units, characterized by a limited number of coding genes and an exceptionally high rate of overlapping genes. We sought a unified evolutionary explanation that accounts for their genome sizes, gene overlapping and capsid properties.

Results

We performed an unbiased statistical analysis of ~100 families within ~400 genera that comprise the currently known viral world. We found that the volume utilization of capsids is often low, and greatly varies among viral families. Furthermore, although viruses span three orders of magnitude in genome length, they almost never have over 1500 overlapping nucleotides, or over four significantly overlapping genes per virus.

Conclusions

Our findings undermine the generality of the compression theory, which emphasizes optimal packing and length dependency to explain overlapping genes and capsid size in viral genomes. Instead, we propose that gene novelty and evolution exploration offer better explanations to size constraints and gene overlapping in all viruses.

Reviewers

This article was reviewed by Arne Elofsson and David Kreil.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-016-0128-3) contains supplementary material, which is available to authorized users.

Determination of antigenicity-altering patches on the major surface protein of human influenza A/H3N2 viruses

Human influenza viruses are rapidly evolving RNA viruses that cause short-term respiratory infections with substantial morbidity and mortality in annual epidemics. Uncovering the general principles of viral coevolution with human hosts is important for pathogen surveillance and vaccine design. Protein regions are an appropriate model for the interactions between two macromolecules, but the currently used epitope definition for the major antigen of influenza viruses, namely hemagglutinin, is very broad. Here, we combined genetic, evolutionary, antigenic, and structural information to determine the most relevant regions of the hemagglutinin of human influenza A/H3N2 viruses for interaction with human immunoglobulins. We estimated the antigenic weights of amino acid changes at individual sites from hemagglutination inhibition data using antigenic tree inference followed by spatial clustering of antigenicity-altering protein sites on the protein structure. This approach determined six relevant areas (patches) for antigenic variation that had a key role in the past antigenic evolution of the viruses. Previous transitions between successive predominating antigenic types of H3N2 viruses always included amino acid changes in either the first or second antigenic patch. Interestingly, there was only partial overlap between the antigenic patches and the patches under strong positive selection. Therefore, besides alterations of antigenicity, other interactions with the host may shape the evolution of human influenza A/H3N2 viruses.

Antigenic evolution of H9N2 chicken influenza viruses isolated in China during 2009-2013 and selection of a candidate vaccine strain with broad cross-reactivity

We previously demonstrated that H9N2 subtype avian influenza viruses (AIVs) isolated from 1994 to 2008 evolved into distinct antigenic groups (C, D, and E) and then underwent antigenic drift from commercial vaccines, causing a country-wide outbreak during 2010-2013. In this study, H9N2 AIVs isolated from chickens during 2009-2013 were antigenically analyzed by performing hemagglutination inhibition and neutralization assays using a panel of polyclonal antibodies. Our findings confirmed the antigenic drift of recent H9N2 viruses from the commercial vaccine and showed that most of these antigenic variants form a novel HI antigenic group, F, with a few belonging to groups D and E. Slight antigenic variation was observed in group F viruses. Genetic analysis of amino acid sequences deduced from hemagglutinin (HA) gene sequences indicated that 9 of 15 mutations predominant in the 2009-2013 viruses can be mapped to known antigenic sites, which might be responsible for the novel antigenicity of group F. These antigenic changes make it necessary to modify the influenza vaccine to ensure efficient protection. A vaccine candidate, Ck/HeB/YT/10, was selected and provided significant protection against viruses from different antigenic groups in terms of reduction in virus shedding, suggesting broad cross-reactivity. Taken together, our results indicate that the H9N2 chicken influenza viruses in China have evolved from distinct antigenic groups into a novel group F that became dominant during the country-wide outbreak and now seems to be undergoing new antigenic divergence. Systematic surveillance and timely updating of vaccine strains are important for viral prevention and control in the future.

Evolution of the H9N2 influenza genotype that facilitated the genesis of the novel H7N9 virus

The emergence of human infection with a novel H7N9 influenza virus in China raises a pandemic concern. Chicken H9N2 viruses provided all six of the novel reassortant's internal genes. However, it is not fully understood how the prevalence and evolution of these H9N2 chicken viruses facilitated the genesis of the novel H7N9 viruses. Here we show that over more than 10 y of cocirculation of multiple H9N2 genotypes, a genotype (G57) emerged that had changed antigenicity and improved adaptability in chickens. It became predominant in vaccinated farm chickens in China, caused widespread outbreaks in 2010-2013 before the H7N9 viruses emerged in humans, and finally provided all of their internal genes to the novel H7N9 viruses. The prevalence and variation of H9N2 influenza virus in farmed poultry could provide an important early warning of the emergence of novel reassortants with pandemic potential.

Reassortment patterns of avian influenza virus internal segments among different subtypes

Background

The segmented RNA genome of avian Influenza viruses (AIV) allows genetic reassortment between co-infecting viruses, providing an evolutionary pathway to generate genetic innovation. The genetic diversity (16 haemagglutinin and 9 neuraminidase subtypes) of AIV indicates an extensive reservoir of influenza viruses exists in bird populations, but how frequently subtypes reassort with each other is still unknown. Here we quantify the reassortment patterns among subtypes in the Eurasian avian viral pool by reconstructing the ancestral states of the subtypes as discrete states on time-scaled phylogenies with respect to the internal protein coding segments. We further analyzed how host species, the inferred evolutionary rates and the d_N/d_S ratio varied among segments and between discrete subtypes, and whether these factors may be associated with inter-subtype reassortment rate.

Results

The general patterns of reassortment are similar among five internal segments with the exception of segment 8, encoding the Non-Structural genes, which has a more divergent phylogeny. However, significant variation in rates between subtypes was observed. In particular, hemagglutinin-encoding segments of subtypes H5 to H9 reassort at a lower rate compared to those of H1 to H4, and Neuraminidase-encoding segments of subtypes N1 and N2 reassort less frequently than N3 to N9. Both host species and d_N/d_S ratio were significantly associated with reassortment rate, while evolutionary rate was not associated. The d_N/d_S ratio was negatively correlated with reassortment rate, as was the number of negatively selected sites for all segments.

Conclusions

These results indicate that overall selective constraint and host species are both associated with reassortment rate. These results together identify the wild bird population as the major source of new reassortants, rather than domestic poultry. The lower reassortment rates observed for H5N1 and H9N2 may be explained by the large proportion of strains derived from domestic poultry populations. In contrast, the higher rates observed in the H1N1, H3N8 and H4N6 subtypes could be due to their primary origin as infections of wild birds with multiple low pathogenicity strains in the large avian reservoir.

Role of genomic and proteomic tools in the study of host-virus interactions and virus evolution

Viruses have short replication cycles and produce genomic variants within a host, a process that seems to adapt to their specific host and also enable them to infect new hosts. The recent emergence of viral genomic variants from the circulating pool within the host population and re-emergence of the old ones are posing serious threat to agriculture, animal husbandry and humanity as a whole. This review assesses the potential role of genomic and proteomic tools that can monitor not only the course of infection and pathogenesis, but also predict the pandemic or zoonotic epidemic potential of a virus in a previously exposed or immunologically naive biological population.

Loss and gain of N-linked glycosylation sites in globular head and stem of HA found in A/H3N2 flu fatal and severe cases during 2013 Tunisia flu seasonal survey

Glycosylation on the globular head of the hemagglutinin (HA) protein of influenza virus acts as an important target for recognition and destruction of virus by innate immune proteins of the collectin family. In the current study, we have characterized the dynamic amino acid changes at N-linked glycosylation sites of full length sequences of HA genes of 5 A/H3N2 Tunisian strains isolates from mild, severe, and fatal cases. Compared to the reference strain, A/Perth/16/2009 substitutions in potential N-glycosylation sites were observed in 5 HA genes at five different positions (45, 124, 128, 144, and 145) generating the losses and gains of N-linked glycosylation sites. Also the mutation N145S was presented in the receptor-binding site of all segments analyzed. Point mutations in several positions in the gene encoding the H3 of Tunisian strains were shown to ablate a glycan attachment site and also loss of a potential glycosylation site. The relation between these mutations and virulence of influenza A/H3N2 virus needed to be verified in the further experiments.

Evolving gene targets and technology in influenza detection

Influenza viruses cause recurring epidemic outbreaks every year associated with high morbidity and mortality. Despite extensive research and surveillance efforts to control influenza outbreaks, the primary mitigation treatment for influenza is the development of yearly vaccine mixes targeted for the most prevalent virus strains. Consequently, the focus of many detection technologies has evolved toward accurate identification of subtype and understanding the evolution and molecular determinants of novel and pathogenic forms of influenza. The recent availability of potential antiviral treatments are only effective if rapid and accurate diagnostic tests for influenza epidemic management are available; thus, early detection of influenza infection is still important for prevention, containment, patient management, and infection control. This review discusses the current and emerging technologies for detection and strain identification of influenza virus and their specific gene targets, as well as their implications in patient management.

PTree: pattern-based, stochastic search for maximum parsimony phylogenies

Phylogenetic reconstruction is vital to analyzing the evolutionary relationship of genes within and across populations of different species. Nowadays, with next generation sequencing technologies producing sets comprising thousands of sequences, robust identification of the tree topology, which is optimal according to standard criteria such as maximum parsimony, maximum likelihood or posterior probability, with phylogenetic inference methods is a computationally very demanding task. Here, we describe a stochastic search method for a maximum parsimony tree, implemented in a software package we named PTree. Our method is based on a new pattern-based technique that enables us to infer intermediate sequences efficiently where the incorporation of these sequences in the current tree topology yields a phylogenetic tree with a lower cost. Evaluation across multiple datasets showed that our method is comparable to the algorithms implemented in PAUP* or TNT, which are widely used by the bioinformatics community, in terms of topological accuracy and runtime. We show that our method can process large-scale datasets of 1,000-8,000 sequences. We believe that our novel pattern-based method enriches the current set of tools and methods for phylogenetic tree inference. The software is available under: http://algbio.cs.uni-duesseldorf.de/webapps/wa-download/.

Simplified large-scale Sanger genome sequencing for influenza A/H3N2 virus

Background

The advent of next-generation sequencing technologies and the resultant lower costs of sequencing have enabled production of massive amounts of data, including the generation of full genome sequences of pathogens. However, the small genome size of the influenza virus arguably justifies the use of the more conventional Sanger sequencing technology which is still currently more readily available in most diagnostic laboratories.

Results

We present a simplified Sanger-based genome sequencing method for sequencing the influenza A/H3N2 virus in a large-scale format. The entire genome sequencing was completed with 19 reverse transcription-polymerase chain reactions (RT-PCRs) and 39 sequencing reactions. This method was tested on 15 native clinical samples and 15 culture isolates, respectively, collected between 2009 and 2011. The 15 native clinical samples registered quantification cycle values ranging from 21.0 to 30.56, which were equivalent to 2.4×10³–1.4×10⁶ viral copies/µL of RNA extract. All the PCR-amplified products were sequenced directly without PCR product purification. Notably, high quality sequencing data up to 700 bp were generated for all the samples tested. The completed sequence covered 408,810 nucleotides in total, with 13,627 nucleotides per genome, attaining 100% coding completeness. Of all the bases produced, an average of 89.49% were Phred quality value 40 (QV40) bases (representing an accuracy of circa one miscall for every 10,000 bases) or higher, and an average of 93.46% were QV30 bases (one miscall every 1000 bases) or higher.

Conclusions

This sequencing protocol has been shown to be cost-effective and less labor-intensive in obtaining full influenza genomes. The constant high quality of sequences generated imparts confidence in extending the application of this non-purified amplicon sequencing approach to other gene sequencing assays, with appropriate use of suitably designed primers.

Design of a set of probes with high potential for influenza virus epidemiological surveillance

An Influenza Probe Set (IPS) consisting in 1,249 9-mer probes for genomic fingerprinting of closely and distantly related Influenza Virus strains was designed and tested in silico. The IPS was derived from alignments of Influenza genomes. The RNA segments of 5,133 influenza strains having diverse degree of relatedness were concatenated and aligned. After alignment, 9-mer sites having high Shannon entropy were searched. Additional criteria such as: G+C content between 35 to 65%, absence of dimer or trimer consecutive repeats, a minimum of 2 differences between 9mers and selecting only sequences with Tm values between 34.5 and 36.5oC were applied for selecting probes with high sequential entropy. Virtual Hybridization was used to predict Genomic Fingerprints to assess the capability of the IPS to discriminate between influenza and related strains. Distance scores between pairs of Influenza Genomic Fingerprints were calculated, and used for estimating Taxonomic Trees. Visual examination of both Genomic Fingerprints and Taxonomic Trees suggest that the IPS is able to discriminate between distant and closely related Influenza strains. It is proposed that the IPS can be used to investigate, by virtual or experimental hybridization, any new, and potentially virulent, strain.

Multiple infections and the evolution of virulence

Infections that consist of multiple parasite strains or species are common in the wild and are a major public health concern. Theory suggests that these infections have a key influence on the evolution of infectious diseases and, more specifically, on virulence evolution. However, we still lack an overall vision of the empirical support for these predictions. We argue that within-host interactions between parasites largely determine how virulence evolves and that experimental data support model predictions. Then, we explore the main limitation of the experimental study of such 'mixed infections', which is that it draws conclusions on evolutionary outcomes from studies conducted at the individual level. We also discuss differences between coinfections caused by different strains of the same species or by different species. Overall, we argue that it is possible to make sense out of the complexity inherent to multiple infections and that experimental evolution settings may provide the best opportunity to further our understanding of virulence evolution.

Development of neuraminidase subtype-specific reference antisera by recombinant protein expressed in baculovirus

Outbreaks of avian influenza A virus infection, particularly the H5N1 strains that have affected birds and some humans for the past 15 years, have highlighted the need for increased surveillance and disease control. Such measures require diagnostic tests to detect and characterize the different subtypes of influenza virus. In the current study, a simple method for producing reference avian influenza virus antisera to be used in diagnostic tests was developed. Antisera of nine avian influenza A virus neuraminidases (NA) used for NA subtyping were produced using a recombinant baculovirus. The recombinant NA (rNA) proteins were expressed in Sf9 insect cells and inoculated intramuscularly into specific-pathogen-free chickens with the ISA70 adjuvant. The NA inhibition antibody titers of the rNA antiserum were in the ranges of 5 to 8 and 6 to 9 log(2) units after the primary and boost immunizations, respectively. The antisera were subtype specific, showing low cross-reactivity against every other NA subtype using the conventional thiobarbituric acid NA inhibition assay. These results suggest that this simple method for producing reference NA antisera without purification may be useful for the diagnosis and surveillance of influenza virus.

Reproductive clonality of pathogens: a perspective on pathogenic viruses, bacteria, fungi, and parasitic protozoa

We propose that clonal evolution in micropathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure, a definition already widely used for all kinds of pathogens, although not clearly formulated by many scientists and rejected by others. The two main manifestations of clonal evolution are strong linkage disequilibrium (LD) and widespread genetic clustering ("near-clading"). We hypothesize that this pattern is not mainly due to natural selection, but originates chiefly from in-built genetic properties of pathogens, which could be ancestral and could function as alternative allelic systems to recombination genes ("clonality/sexuality machinery") to escape recombinational load. The clonal framework of species of pathogens should be ascertained before any analysis of biomedical phenotypes (phylogenetic character mapping). In our opinion, this model provides a conceptual framework for the population genetics of any micropathogen.

Inferring the inter-host transmission of influenza A virus using patterns of intra-host genetic variation

Influenza A viruses (IAVs) cause acute, highly transmissible infections in a wide range of animal species. Understanding how these viruses are transmitted within and between susceptible host populations is critical to the development of effective control strategies. While viral gene sequences have been used to make inferences about IAV transmission dynamics at the epidemiological scale, their utility in accurately determining patterns of inter-host transmission in the short-term--i.e. who infected whom--has not been strongly established. Herein, we use intra-host sequence data from the viral HA1 (hemagglutinin) gene domain from two transmission studies employing different IAV subtypes in their natural hosts--H3N8 in horses and H1N1 in pigs-to determine how well these data recapitulate the known pattern of inter-host transmission. Although no mutations were fixed over the course of either experimental transmission chain, we show that some minor, transient alleles can provide evidence of host-to-host transmission and, importantly, can be distinguished from those that cannot.

Evasion of influenza A viruses from innate and adaptive immune responses

The influenza A virus is one of the leading causes of respiratory tract infections in humans. Upon infection with an influenza A virus, both innate and adaptive immune responses are induced. Here we discuss various strategies used by influenza A viruses to evade innate immune responses and recognition by components of the humoral and cellular immune response, which consequently may result in reduced clearing of the virus and virus-infected cells. Finally, we discuss how the current knowledge about immune evasion can be used to improve influenza A vaccination strategies.

Field monitoring of avian influenza viruses: whole-genome sequencing and tracking of neuraminidase evolution using 454 pyrosequencing

Adaptation of avian influenza viruses (AIVs) from waterfowl to domestic poultry with a deletion in the neuraminidase (NA) stalk has already been reported. The way the virus undergoes this evolution, however, is thus far unclear. We address this question using pyrosequencing of duck and turkey low-pathogenicity AIVs. Ducks and turkeys were sampled at the very beginning of an H6N1 outbreak, and turkeys were swabbed again 8 days later. NA stalk deletions were evidenced in turkeys by Sanger sequencing. To further investigate viral evolution, 454 pyrosequencing was performed: for each set of samples, up to 41,500 reads of ca. 400 bp were generated and aligned. Genetic polymorphisms between duck and turkey viruses were tracked on the whole genome. NA deletion was detected in less than 2% of reads in duck feces but in 100% of reads in turkey tracheal specimens collected at the same time. Further variations in length were observed in NA from turkeys 8 days later. Similarly, minority mutants emerged on the hemagglutinin (HA) gene, with substitutions mostly in the receptor binding site on the globular head. These critical changes suggest a strong evolutionary pressure in turkeys. The increasing performances of next-generation sequencing technologies should enable us to monitor the genomic diversity of avian influenza viruses and early emergence of potentially pathogenic variants within bird flocks. The present study, based on 454 pyrosequencing, suggests that NA deletion, an example of AIV adaptation from waterfowl to domestic poultry, occurs by selection rather than de novo emergence of viral mutants.

Inference of genotype-phenotype relationships in the antigenic evolution of human influenza A (H3N2) viruses

Distinguishing mutations that determine an organism's phenotype from (near-) neutral ‘hitchhikers’ is a fundamental challenge in genome research, and is relevant for numerous medical and biotechnological applications. For human influenza viruses, recognizing changes in the antigenic phenotype and a strains' capability to evade pre-existing host immunity is important for the production of efficient vaccines. We have developed a method for inferring ‘antigenic trees’ for the major viral surface protein hemagglutinin. In the antigenic tree, antigenic weights are assigned to all tree branches, which allows us to resolve the antigenic impact of the associated amino acid changes. Our technique predicted antigenic distances with comparable accuracy to antigenic cartography. Additionally, it identified both known and novel sites, and amino acid changes with antigenic impact in the evolution of influenza A (H3N2) viruses from 1968 to 2003. The technique can also be applied for inference of ‘phenotype trees’ and genotype–phenotype relationships from other types of pairwise phenotype distances.

The molecular evolution of any organism is described by changes in the genotype resulting from genetic drift or selection to maintain or establish fitness under the given environmental conditions. Identification of phenotype-defining changes and their distinction from (near-) neutral (‘hitchhikers’) ones is a fundamental challenge in genome research. The standard approach involves time- and cost-intensive mutation experiments, which are typically low throughput, due to their experimental nature. We have developed a computational method for the inference of phenotypic impact of genotypic changes that is applicable to any system, within or across species, where homologous genetic sequences and associated pairwise phenotype distances are available. We demonstrate the accuracy of our method by application to the human influenza A (H3N2) virus. This exemplary system is of particular interest, as recognizing changes in the antigenic phenotype and a viral strains' capability to evade pre-existing host immunity is important for the production of efficient vaccines. We accurately identified known sites and amino acid changes with antigenic impact over 35 years of evolution, and provide further details on individual antigenically relevant changes in the evolution of influenza A (H3N2) viruses.

Detecting patches of protein sites of influenza A viruses under positive selection

Influenza A viruses are single-stranded RNA viruses capable of evolving rapidly to adapt to environmental conditions. Examples include the establishment of a virus in a novel host or an adaptation to increasing immunity within the host population due to prior infection or vaccination against a circulating strain. Knowledge of the viral protein regions under positive selection is therefore crucial for surveillance. We have developed a method for detecting positively selected patches of sites on the surface of viral proteins, which we assume to be relevant for adaptive evolution. We measure positive selection based on dN/dS ratios of genetic changes inferred by considering the phylogenetic structure of the data and suggest a graph-cut algorithm to identify such regions. Our algorithm searches for dense and spatially distinct clusters of sites under positive selection on the protein surface. For the hemagglutinin protein of human influenza A viruses of the subtypes H3N2 and H1N1, our predicted sites significantly overlap with known antigenic and receptor-binding sites. From the structure and sequence data of the 2009 swine-origin influenza A/H1N1 hemagglutinin and PB2 protein, we identified regions that provide evidence of evolution under positive selection since introduction of the virus into the human population. The changes in PB2 overlap with sites reported to be associated with mammalian adaptation of the influenza A virus. Application of our technique to the protein structures of viruses of yet unknown adaptive behavior could identify further candidate regions that are important for host–virus interaction.

L-ficolin binds to the glycoproteins hemagglutinin and neuraminidase and inhibits influenza A virus infection both in vitro and in vivo

L-ficolin, one of the complement lectins found in human serum, is a novel pattern recognition molecule that can specifically bind to microbial carbohydrates, thereby activating the lectin complement pathway and mounting a protective innate immune response. However, little is known about the role of L-ficolin during viral infections in vivo. In the present study, we used a mouse model of influenza A virus infection to demonstrate that the administration of exogenous L-ficolin or ficolin A (FCNA - an L-ficolin-like molecule in the mouse) is protective against the virus. Furthermore, FCNA-null mice have a greatly increased susceptibility to infection with the influenza A virus. Moreover, we found recombinant human L-ficolin inhibited influenza A virus entry into Madin-Darby canine kidney cells. More importantly, L-ficolin can recognize and bind hemagglutinin (HA) and neuraminidase (NA) glycoproteins and different subtypes of influenza A virus, and these interactions can be competitively inhibited by N-acetyl-D-glucosamine. In addition, the binding of L-ficolin and FCNA may lead to the activation of the lectin complement pathway. To our knowledge, this is the first report demonstrating that L-ficolin can block influenza virus infections both in vitro and in vivo using FCNA-knockout mice, possibly by interacting with the carbohydrates of HA and NA. Therefore, these data may provide new immunotherapeutic strategies based on the innate immune molecule L-ficolin against the influenza A virus.

Alignment free characterization of the influenza-A hemagglutinin genes by the ISSCOR method

Analyses and visualizations by the ISSCOR method of the influenza virus hemagglutinin genes of three different A-subtypes revealed some rather striking temporal (for A/H3N3), and spatial relationships (for A/H5N1) between groups of individual gene subsets. The application to the A/H1N1 set revealed also relationships between the seasonal H1, and the swine-like novel 2009 H1v variants in a quick and unambiguous manner. Based on these examples we consider the application of the ISSCOR method for analysis of large sets of homologous genes as a worthwhile addition to a toolbox of genomics-it allows a rapid diagnostics of trends, and possibly can even aid an early warning of newly emerging epidemiological threats.

Evidence for N-glycan shielding of antigenic sites during evolution of human influenza A virus hemagglutinin

After the emergence of influenza A viruses in the human population, the number of N-glycosylation sites (NGS) in the globular head region of hemagglutinin (HA) has increased continuously for several decades. It has been speculated that the addition of NGS to the globular head region of HA has conferred selective advantages to the virus by preventing the binding of antibodies (Ab) to antigenic sites (AS). Here, the effect of N-glycosylation on the binding of Ab to AS in human influenza A virus subtype H3N2 (A/H3N2) was examined by inferring natural selection at AS and other sites (NAS) that are located close to and distantly from the NGS in the three-dimensional structure of HA through a comparison of the rates of synonymous (d(S)) and nonsynonymous (d(N)) substitutions. When positions 63, 122, 126, 133, 144, and 246 in the globular head region of HA were non-NGS, the d(N)/d(S) was >1 and positive selection was detected at the AS located near these positions. However, the d(N)/d(S) value decreased and the evidence of positive selection disappeared when these positions became NGS. In contrast, d(N)/d(S) at the AS distantly located from the positions mentioned above and at the NAS of any location were generally <1 and did not decrease when these positions changed from non-NGS to NGS. These results suggest that the attachment of N-glycans to the NGS in the globular head region of HA prevented the binding of Ab to AS in the evolutionary history of human A/H3N2 virus.

Molecular epidemiology of influenza A/H3N2 viruses circulating in Uganda

The increasing availability of complete influenza virus genomes is deepening our understanding of influenza evolutionary dynamics and facilitating the selection of vaccine strains. However, only one complete African influenza virus sequence is available in the public domain. Here we present a complete genome analysis of 59 influenza A/H3N2 viruses isolated from humans in Uganda during the 2008 and 2009 season. Isolates were recovered from hospital-based sentinel surveillance for influenza-like illnesses and their whole genome sequenced. The viruses circulating during these two seasons clearly differed from each other phylogenetically. They showed a slow evolution away from the 2009/10 recommended vaccine strain (A/Brisbane/10/07), instead clustering with the 2010/11 recommended vaccine strain (A/Perth/16/09) in the A/Victoria/208/09 clade, as observed in other global regions. All of the isolates carried the adamantane resistance marker S31N in the M2 gene and carried several markers of enhanced transmission; as expected, none carried any marker of neuraminidase inhibitor resistance. The hemagglutinin gene of the 2009 isolates differed from that of the 2008 isolates in antigenic sites A, B, D, and to a lesser extent, C and E indicating evidence of an early phylogenetic shift from the 2008 to 2009 viruses. The internal genes of the 2009 isolates were similar to those of one 2008 isolate, A/Uganda/MUWRP-050/2008. Another 2008 isolate had a truncated PB1-F2 protein. Whole genome sequencing can enhance surveillance of future seasonal changes in the viral genome which is crucial to ensure that selected vaccine strains are protective against the strains circulating in Eastern Africa. This data provides an important baseline for this surveillance. Overall the influenza virus activity in Uganda appears to mirror that observed in other regions of the southern hemisphere.

Prospective of Genomics in Revealing Transmission, Reassortment and Evolution of Wildlife-Borne Avian Influenza A (H5N1) Viruses

The outbreak of highly pathogenic avian influenza (HPAI) H5N1 disease has led to significant loss of poultry and wild life and case fatality rates in humans of 60%. Wild birds are natural hosts for all avian influenza virus subtypes and over120 bird species have been reported with evidence of H5N1 infection. Influenza A viruses possess a segmented RNA genome and are characterized by frequently occurring genetic reassortment events, which play a very important role in virus evolution and the spread of novel gene constellations in immunologically naïve human and animal populations. Phylogenetic analysis of whole genome or sub-genomic sequences is a standard means for delineating genetic variation, novel reassortment events, and surveillance to trace the global transmission pathways. In this paper, special emphasis is given to the transmission and circulation of H5N1 among wild life populations, and to the reassortment events that are associated with inter-host transmission of the H5N1 viruses when they infect different hosts, such as birds, pigs and humans. In addition, we review the inter-subtype reassortment of the viral segments encoding inner proteins between the H5N1 viruses and viruses of other subtypes, such as H9N2 and H6N1. Finally, we highlight the usefulness of genomic sequences in molecular epidemiological analysis of HPAI H5N1 and the technical limitations in existing analytical methods that hinder them from playing a greater role in virological research.