The influenza virus resource at the National Center for Biotechnology Information

High performance screening, structural and molecular dynamics analysis to identify H1 inhibitors from TCM Database@Taiwan

New-type oseltamivir-resistant H1N1 influenza viruses have been a major threat to human health since the 2009 flu pandemic. To resolve the drug resistance issue, we aimed to identify a new type of inhibitors against H1 from traditional Chinese medicine (TCM) by employing the world's largest TCM database () for virtual screening and molecular dynamics (MD). From the virtual screening results, sodium (+)-isolaricireinol-2 alpha-sulfate, sodium 3,4-dihydroxy-5-methoxybenzoic acid methyl ester-4-sulfate, sodium (E)-7-hydroxy-1,7-bis(4-hydroxyphenyl)hept-5-ene-3S-sulfonate, and 3-methoxytyramine-betaxanthin were identified as potential drug-like compounds. MD simulation of the binding poses with the key residues Asp103 and Glu83, as well as other binding site residues, identified higher numbers of hydrogen bonds than N-Acetyl-D-Glucosamine (NAG), the natural ligand of the esterase domain in H1. Ionic bonds, salt bridges, and electrostatic energy also contribute to binding stability. Key binding residues include Lys71, Glu83, Asp103, and Arg238. Structural moieties promoting H-bond or salt bridge formations at these locations greatly contribute to a stable ligand-protein complex. An available sodium atom for ionic interactions with Asp103 can further stabilize the ligands. Based on virtual screening, MD simulation, and interaction energy evaluation, TCM candidates demonstrate good potential as novel H1 inhibitors. In addition, the identified stabilizing features can provide insights for designing highly stable H1 inhibitors.

Migratory flyway and geographical distance are barriers to the gene flow of influenza virus among North American birds

Despite the importance of migratory birds in the ecology and evolution of avian influenza virus (AIV), there is a lack of information on the patterns of AIV spread at the intra-continental scale. We applied a variety of statistical phylogeographic techniques to a plethora of viral genome sequence data to determine the strength, pattern and determinants of gene flow in AIV sampled from wild birds in North America. These analyses revealed a clear isolation-by-distance of AIV among sampling localities. In addition, we show that phylogeographic models incorporating information on the avian flyway of sampling proved a better fit to the observed sequence data than those specifying homogeneous or random rates of gene flow among localities. In sum, these data strongly suggest that the intra-continental spread of AIV by migratory birds is subject to major ecological barriers, including spatial distance and avian flyway.

Gnarled-trunk evolutionary model of influenza A virus hemagglutinin

Human influenza A viruses undergo antigenic changes with gradual accumulation of amino acid substitutions on the hemagglutinin (HA) molecule. A strong antigenic mismatch between vaccine and epidemic strains often requires the replacement of influenza vaccines worldwide. To establish a practical model enabling us to predict the future direction of the influenza virus evolution, relative distances of amino acid sequences among past epidemic strains were analyzed by multidimensional scaling (MDS). We found that human influenza viruses have evolved along a gnarled evolutionary pathway with an approximately constant curvature in the MDS-constructed 3D space. The gnarled pathway indicated that evolution on the trunk favored multiple substitutions at the same amino acid positions on HA. The constant curvature was reasonably explained by assuming that the rate of amino acid substitutions varied from one position to another according to a gamma distribution. Furthermore, we utilized the estimated parameters of the gamma distribution to predict the amino acid substitutions on HA in subsequent years. Retrospective prediction tests for 12 years from 1997 to 2009 showed that 70% of actual amino acid substitutions were correctly predicted, and that 45% of predicted amino acid substitutions have been actually observed. Although it remains unsolved how to predict the exact timing of antigenic changes, the present results suggest that our model may have the potential to recognize emerging epidemic strains.

[Bioinformatics technologies for the analysis of antigenic evolution of influenza viruses]

Human influenza viruses mutate from time to time, causing annual epidemics worldwide. The strong immune pressure in the human population selects a new variant every year, and the antigenic change is one of the primary reasons why vaccination is not a perfect measure to control seasonal influenza. Thus prediction of antigenic change of influenza A virus has been one of the major public health goals. In this review bioinformatics technologies that have been developed to achieve this goal were summarized.

Reassortment events among swine influenza A viruses in China: implications for the origin of the 2009 influenza pandemic

That pigs may play a pivotal role in the emergence of pandemic influenza was indicated by the recent H1N1/2009 human pandemic, likely caused by a reassortant between viruses of the American triple-reassortant (TR) and Eurasian avian-like (EA) swine influenza lineages. As China has the largest human and pig populations in the world and is the only place where both TR and EA viruses have been reported to cocirculate, it is potentially the source of the H1N1/2009 pandemic virus. To examine this, the genome sequences of 405 swine influenza viruses from China were analyzed. Thirty-six TR and EA reassortant viruses were identified before and after the occurrence of the pandemic. Several of these TR-EA reassortant viruses had genotypes with most segments having the same lineage origin as the segments of the H1N1/2009 pandemic virus. However, these viruses were generated from independent reassortment events throughout our survey period and were not associated with the current pandemic. One TR-EA reassortant, which is least similar to the pandemic virus, has persisted since 2007, while all the other variants appear to be transient. Despite frequent reassortment events between TR and EA lineage viruses in China, evidence for the genesis of the 2009 pandemic virus in pigs in this region is still absent.

BFluenza: A Proteomic Database on Bird Flu

Influenza A virus subtype H5N1, also known as "bird flu" has been documented to cause an outbreak of respiratory diseases in humans. The unprecedented spread of highly pathogenic avian influenza type A is a threat to veterinary and human health. The BFluenza is a relational database which is solely devoted to proteomic information of H5N1 subtype. Bfluenza has novel features including computed physico-chemical properties data of H5N1 viral proteins, modeled structures of viral proteins, data of protein coordinates, experimental details, molecular description and bibliographic reference. The database also contains nucleotide and their decoded protein sequences data. The database can be searched in various modes by setting search options. The structure of viral protein could be visualized by JMol viewer or by Discovery Studio.

AVAILABILITY

The database is available for free at http://www.bfluenza.info.

Environmental levels of the antiviral oseltamivir induce development of resistance mutation H274Y in influenza A/H1N1 virus in mallards

Oseltamivir (Tamiflu®) is the most widely used drug against influenza infections and is extensively stockpiled worldwide as part of pandemic preparedness plans. However, resistance is a growing problem and in 2008–2009, seasonal human influenza A/H1N1 virus strains in most parts of the world carried the mutation H274Y in the neuraminidase gene which causes resistance to the drug. The active metabolite of oseltamivir, oseltamivir carboxylate (OC), is poorly degraded in sewage treatment plants and surface water and has been detected in aquatic environments where the natural influenza reservoir, dabbling ducks, can be exposed to the substance. To assess if resistance can develop under these circumstances, we infected mallards with influenza A/H1N1 virus and exposed the birds to 80 ng/L, 1 µg/L and 80 µg/L of OC through their sole water source. By sequencing the neuraminidase gene from fecal samples, we found that H274Y occurred at 1 µg/L of OC and rapidly dominated the viral population at 80 µg/L. IC₅₀ for OC was increased from 2–4 nM in wild-type viruses to 400–700 nM in H274Y mutants as measured by a neuraminidase inhibition assay. This is consistent with the decrease in sensitivity to OC that has been noted among human clinical isolates carrying H274Y. Environmental OC levels have been measured to 58–293 ng/L during seasonal outbreaks and are expected to reach µg/L-levels during pandemics. Thus, resistance could be induced in influenza viruses circulating among wild ducks. As influenza viruses can cross species barriers, oseltamivir resistance could spread to human-adapted strains with pandemic potential disabling oseltamivir, a cornerstone in pandemic preparedness planning. We propose surveillance in wild birds as a measure to understand the resistance situation in nature and to monitor it over time. Strategies to lower environmental levels of OC include improved sewage treatment and, more importantly, a prudent use of antivirals.

Rapid PCR-based molecular pathotyping of H5 and H7 avian influenza viruses

While the majority of avian influenza virus (AIV) subtypes are classified as low-pathogenicity avian influenza viruses (LPAIV), the H5 and H7 subtypes have the ability to mutate to highly pathogenic avian influenza viruses (HPAIV) in poultry and therefore are the etiological agents of notifiable AIV (NAIV). It is of great importance to distinguish HPAIV from LPAIV variants during H5/H7 outbreaks and surveillance. To this end, a novel and fast strategy for the molecular pathotyping of H5/H7 AIVs is presented. The differentiation of the characteristic hemagglutinin (HA) protein cleavage sites (CSs) of HPAIVs and LPAIVs is achieved by a novel PCR method where the samples are interrogated for all existing CSs with a 484-plex primer mixture directly targeting the CS region. CSs characteristic for HP or LP H5/H7 viruses are distinguished in a seminested duplex real-time PCR format using plexor fluorogenic primers. Eighty-six laboratory isolates and 60 characterized NAIV-positive clinical specimens from poultry infected with H5/H7 both experimentally and in the field were successfully pathotyped in the validation. The method has the potential to substitute CS sequencing in the HA gene for the determination of the molecular pathotype, thereby providing a rapid means to acquire additional information concerning NAIV outbreaks, which may be critical to their management. The new assay may be extended to the LP/HP differentiation of previously unknown H5/H7 isolates. It may be considered for integration into surveillance and control programs in both domestic and wild bird populations.

The quality preserving database: a computational framework for encouraging collaboration, enhancing power and controlling false discovery

The common scenario in computational biology in which a community of researchers conduct multiple statistical tests on one shared database gives rise to the multiple hypothesis testing problem. Conventional procedures for solving this problem control the probability of false discovery by sacrificing some of the power of the tests. We suggest a scheme for controlling false discovery without any power loss by adding new samples for each use of the database and charging the user with the expenses. The crux of the scheme is a carefully crafted pricing system that fairly prices different user requests based on their demands while keeping the probability of false discovery bounded. We demonstrate this idea in the context of HIV treatment research, where multiple researchers conduct tests on a repository of HIV samples.

Evolutionary pathways of the pandemic influenza A (H1N1) 2009 in the UK

The emergence of the influenza (H1N1) 2009 virus provided a unique opportunity to study the evolution of a pandemic virus following its introduction into the human population. Virological and clinical surveillance in the UK were comprehensive during the first and second waves of the pandemic in 2009, with extensive laboratory confirmation of infection allowing a detailed sampling of representative circulating viruses. We sequenced the complete coding region of the haemagglutinin (HA) segment of 685 H1N1 pandemic viruses selected without bias during two waves of pandemic in the UK (April-December 2009). Phylogenetic analysis showed that although temporal accumulation of amino acid changes was observed in the HA sequences, the overall diversity was less than that typically seen for seasonal influenza A H1N1 or H3N2. There was co-circulation of multiple variants as characterised by signature amino acid changes in the HA. A specific substitution (S203T) became predominant both in UK and global isolates. No antigenic drift occurred during 2009 as viruses with greater than four-fold reduction in their haemagglutination inhibition (HI) titre ("low reactors") were detected in a low proportion (3%) and occurred sporadically. Although some limited antigenic divergence in viruses with four-fold reduction in HI titre might be related to the presence of 203T, additional studies are needed to test this hypothesis.

Putative amino acid determinants of the emergence of the 2009 influenza A (H1N1) virus in the human population

The emergence of the unique H1N1 influenza A virus in 2009 resulted in a pandemic that has spread to over 200 countries. The constellation of molecular factors leading to the emergence of this strain is still unclear. Using a computational approach, we identified molecular determinants that may discriminate the hemagglutinin protein of the 2009 human pandemic H1N1 (pH1N1) strain from that of other H1N1 strains. As expected, positions discriminating the pH1N1 from seasonal human strains were located in or near known H1N1 antigenic sites, thus camouflaging the pH1N1 strain from immune recognition. For example, the alteration S145K (an antigenic position) was found as a characteristic of the pH1N1 strain. We also detected positions in the hemagglutinin protein differentiating classical swine viruses from pH1N1. These positions were mostly located in and around the receptor-binding pocket, possibly influencing binding affinity to the human cell. Such alterations may be liable in part for the virus's efficient infection and adaptation to humans. For instance, 133(A) and 149 were identified as discriminative positions. Significantly, we showed that the substitutions R133(A)K and R149K, predicted to be pH1N1 characteristics, each altered virus binding to erythrocytes and conferred virulence to A/swine/NC/18161/02 in mice, reinforcing the computational findings. Our findings provide a structural explanation for the deficient immunity of humans to the pH1N1 strain. Moreover, our analysis points to unique molecular factors that may have facilitated the emergence of this swine variant in humans, in contrast to other swine variants that failed.

Glycosylation site alteration in the evolution of influenza A (H1N1) viruses

Influenza virus typically alters protein glycosylation in order to escape immune pressure from hosts and hence to facilitate survival in different host environments. In this study, the patterns and conservation of glycosylation sites on HA and NA of influenza A/H1N1 viruses isolated from various hosts at different time periods were systematically analyzed, by employing a new strategy combining genome-based glycosylation site prediction and 3D modeling of glycoprotein structures, for elucidation of the modes and laws of glycosylation site alteration in the evolution of influenza A/H1N1 viruses. The results showed that influenza H1N1 viruses underwent different alterations of protein glycosylation in different hosts. Two alternative modes of glycosylation site alteration were involved in the evolution of human influenza virus: One was an increase in glycosylation site numbers, which mainly occurred with high frequency in the early stages of evolution. The other was a change in the positional conversion of the glycosylation sites, which was the dominating mode with relatively low frequency in the later evolutionary stages. The mechanisms and possibly biological functions of glycosylation site alteration for the evolution of influenza A/H1N1 viruses were also discussed. Importantly, the significant role of positional alteration of glycosylation sites in the host adaptation of influenza virus was elucidated. Although the results still need to be supported by experimental data, the information here may provide some constructive suggestions for research into the glycosylation of influenza viruses as well as even the design of surveillance and the production of viral vaccines.

Molecular characteristics of outbreaks of nosocomial infection with influenza A/H3N2 virus variants

OBJECTIVE

To describe outbreaks of nosocomial influenza infection with molecular methods and to elucidate the viral linkages among outbreak case patients including both inpatients and healthcare workers (HCWs).

SETTING

A 180-bed acute and long-term care hospital in Japan.

METHODS

Retrospective observational study of nosocomial outbreaks of infection with influenza A/H3N2. Together with information about onset dates and vaccination history, we obtained nasopharyngeal swab samples from individuals with cases of influenza or influenza-like illness (ILI). The hemagglutinin genes of the recovered viruses were sequenced and compared, along with those of community-circulating strains, for similarity by phylogenetic tree analysis.

RESULTS

The outbreaks occurred from February 26 through April 3, 2007, during the 2006-2007 epidemic season, and they involved 11 patients and 13 HCWs. The 2 outbreaks involved 2 different genotypes of influenza A/H3N2 viruses. These virus variants were closely related to the influenza strains that were circulating in the community during the same epidemic season.

CONCLUSION

This study showed the dissemination of highly homologous influenza virus variants among inpatients and HCWs within a short period, as a result of nosocomial transmission. These strains were also similar to influenza strains that were circulating in the community.

Strength and tempo of selection revealed in viral gene genealogies

Background

RNA viruses evolve extremely quickly, allowing them to rapidly adapt to new environmental conditions. Viral pathogens, such as influenza virus, exploit this capacity for evolutionary change to persist within the human population despite substantial immune pressure. Understanding the process of adaptation in these viral systems is essential to our efforts to combat infectious disease.

Results

Through analysis of simulated populations and sequence data from influenza A (H3N2) and measles virus, we show how phylogenetic and population genetic techniques can be used to assess the strength and temporal pattern of adaptive evolution. The action of natural selection affects the shape of the genealogical tree connecting members of an evolving population, causing deviations from the neutral expectation. The magnitude and distribution of these deviations lends insight into the historical pattern of evolution and adaptation in the viral population. We quantify the degree of ongoing adaptation in influenza and measles virus through comparison of census population size and effective population size inferred from genealogical patterns, finding a 60-fold greater deviation in influenza than in measles. We also examine the tempo of adaptation in influenza, finding evidence for both continuous and episodic change.

Conclusions

Our results have important consequences for understanding the epidemiological and evolutionary dynamics of the influenza virus. Additionally, these general techniques may prove useful to assess the strength and pattern of adaptive evolution in a variety of evolving systems. They are especially powerful when assessing selection in fast-evolving populations, where temporal patterns become highly visible.

A computational-experimental approach identifies mutations that enhance surface expression of an oseltamivir-resistant influenza neuraminidase

The His274→Tyr (H274Y) oseltamivir (Tamiflu) resistance mutation causes a substantial decrease in the total levels of surface-expressed neuraminidase protein and activity in early isolates of human seasonal H1N1 influenza, and in the swine-origin pandemic H1N1. In seasonal H1N1, H274Y only became widespread after the occurrence of secondary mutations that counteracted this decrease. H274Y is currently rare in pandemic H1N1, and it remains unclear whether secondary mutations exist that might similarly counteract the decreased neuraminidase surface expression associated with this resistance mutation in pandemic H1N1. Here we investigate the possibility of predicting such secondary mutations. We first test the ability of several computational approaches to retrospectively identify the secondary mutations that enhanced levels of surface-expressed neuraminidase protein and activity in seasonal H1N1 shortly before the emergence of oseltamivir resistance. We then use the most successful computational approach to predict a set of candidate secondary mutations to the pandemic H1N1 neuraminidase. We experimentally screen these mutations, and find that several of them do indeed partially counteract the decrease in neuraminidase surface expression caused by H274Y. Two of the secondary mutations together restore surface-expressed neuraminidase activity to wildtype levels, and also eliminate the very slight decrease in viral growth in tissue-culture caused by H274Y. Our work therefore demonstrates a combined computational-experimental approach for identifying mutations that enhance neuraminidase surface expression, and describes several specific mutations with the potential to be of relevance to the spread of oseltamivir resistance in pandemic H1N1.

A highly conserved neutralizing epitope on group 2 influenza A viruses

Current flu vaccines provide only limited coverage against seasonal strains of influenza viruses. The identification of V(H)1-69 antibodies that broadly neutralize almost all influenza A group 1 viruses constituted a breakthrough in the influenza field. Here, we report the isolation and characterization of a human monoclonal antibody CR8020 with broad neutralizing activity against most group 2 viruses, including H3N2 and H7N7, which cause severe human infection. The crystal structure of Fab CR8020 with the 1968 pandemic H3 hemagglutinin (HA) reveals a highly conserved epitope in the HA stalk distinct from the epitope recognized by the V(H)1-69 group 1 antibodies. Thus, a cocktail of two antibodies may be sufficient to neutralize most influenza A subtypes and, hence, enable development of a universal flu vaccine and broad-spectrum antibody therapies.

Evolutionary dynamics of influenza A nucleoprotein (NP) lineages revealed by large-scale sequence analyses

Influenza A viral nucleoprotein (NP) plays a critical role in virus replication and host adaptation, however, the underlying molecular evolutionary dynamics of NP lineages are less well-understood. In this study, large-scale analyses of 5094 NP nucleotide sequences revealed eight distinct evolutionary lineages, including three host-specific lineages (human, classical swine and equine), two cross-host lineages (Eurasian avian-like swine and swine-origin human pandemic H1N1 2009) and three geographically isolated avian lineages (Eurasian, North American and Oceanian). The average nucleotide substitution rate of the NP lineages was estimated to be 2.4 × 10(-3) substitutions per site per year, with the highest value observed in pandemic H1N1 2009 (3.4 × 10(-3)) and the lowest in equine (0.9 × 10(-3)). The estimated time of most recent common ancestor (TMRCA) for each lineage demonstrated that the earliest human lineage was derived around 1906, and the latest pandemic H1N1 2009 lineage dated back to December 17, 2008. A marked time gap was found between the times when the viruses emerged and were first sampled, suggesting the crucial role for long-term surveillance of newly emerging viruses. The selection analyses showed that human lineage had six positive selection sites, whereas pandemic H1N1 2009, classical swine, Eurasian avian and Eurasian swine had only one or two sites. Protein structure analyses revealed several positive selection sites located in epitope regions or host adaptation regions, indicating strong adaptation to host immune system pressures in influenza viruses. Along with previous studies, this study provides new insights into the evolutionary dynamics of influenza A NP lineages. Further lineage analyses of other gene segments will allow better understanding of influenza A virus evolution and assist in the improvement of global influenza surveillance.

Structural basis and sequence co-evolution analysis of the hemagglutinin protein of pandemic influenza A/H1N1 (2009) virus

Severe pandemic influenza A H1N1 (2009) infection, especially in the lower respiratory tract, is often associated with the virus carrying a D222G substitution in the hemagglutinin (HA) protein of the virus. The mechanism for this association has not been fully explored. In the in vitro binding assay, it was found that clinical isolates carrying D222G substitution exhibit higher binding avidity to 2,3-linked sialic acids than the wild-type virus. The receptor binding pocket of the pandemic influenza (H1N1) HA was found to be smaller than those of other influenza A strains, allowing tighter binding of the virus with the receptor, yet also inducing steric stress for the binding. Our homology modeling and molecular docking calculations implicated that residue 222 may affect the positioning of the conserved Q223 residue, hence modulating flexibility of the binding pocket and steric hindrance during receptor binding. The molecular property of residue 222 can also directly influence the 'lysine fence' via the polarity of the amino acid residue where D222G substitution will enhance the electrostatic interactions between the receptor and the protein. The potential importance of residue 222 was illustrated by evolutionary analysis, which showed that this site is under intense selection pressure during adaptation of the virus to human host. Our findings provide a useful reference for follow-up studies in monitoring the ongoing evolution of the pandemic influenza A H1N1 (2009) virus.

Pandemic (H1N1) 2009 virus revisited: an evolutionary retrospective

The pandemic (H1N1) 2009 virus is unique in many aspects, especially in its genetics and evolution. In this paper, we examine the molecular mechanisms underlying the evolution of this novel virus through a comprehensive bioinformatics analysis, and present results in the context of a review of the literature. The pandemic virus was found to arise from a reassortment of two swine viruses, each of which ultimately arose from interspecies transmission. It experienced fast evolutionary rates and strong selection pressures, diverging into two different clusters at the early pandemic stage. Cluster I became extinct at the end of 2009 whereas Cluster II continued to circulate at much lower rates in 2010. Therefore, on August 10 of 2010 the WHO declared the end of the pandemic. Important mutations associated with host specificity, virulence, and drug resistance were detected in the pandemic virus, indicating effective transmission and increased severity in humans. Much has been learned about the evolutionary dynamics of this pandemic virus; however, it is still impossible to predict when the next pandemic will occur and which virus will be responsible. Improved surveillance at different levels (both national and international) and in different hosts (especially in swine) appears to be crucial for early detection and prevention of future influenza pandemics.

Extensive geographical mixing of 2009 human H1N1 influenza A virus in a single university community

Despite growing interest in the molecular epidemiology of influenza virus, the pattern of viral spread within individual communities remains poorly understood. To determine the phylogeography of influenza virus in a single population, we examined the spatial diffusion of H1N1/09 influenza A virus within the student body of the University of California, San Diego (UCSD), sampling for a 1-month period between October and November 2009. Despite the highly focused nature of our study, an analysis of complete viral genome sequences revealed between 24 and 33 independent introductions of H1N1/09 into the UCSD community, comprising much of the global genetic diversity in this virus. These data were also characterized by a relatively low level of on-campus transmission as well as extensive spatial mixing, such that there was little geographical clustering by either student residence or city ZIP code. Most notably, students experiencing illness on the same day and residing in the same dorm possessed phylogenetically distinct lineages. H1N1/09 influenza A virus is therefore characterized by a remarkable spatial fluidity, which is likely to impede community-based methods for its control, including class cancellations, quarantine, and chemoprophylaxis.

Host shifts and molecular evolution of H7 avian influenza virus hemagglutinin

Evolutionary consequences of host shifts represent a challenge to identify the mechanisms involved in the emergence of influenza A (IA) viruses. In this study we focused on the evolutionary history of H7 IA virus in wild and domestic birds, with a particular emphasis on host shifts consequences on the molecular evolution of the hemagglutinin (HA) gene. Based on a dataset of 414 HA nucleotide sequences, we performed an extensive phylogeographic analysis in order to identify the overall genetic structure of H7 IA viruses. We then identified host shift events and investigated viral population dynamics in wild and domestic birds, independently. Finally, we estimated changes in nucleotide substitution rates and tested for positive selection in the HA gene. A strong association between the geographic origin and the genetic structure was observed, with four main clades including viruses isolated in North America, South America, Australia and Eurasia-Africa. We identified ten potential events of virus introduction from wild to domestic birds, but little evidence for spillover of viruses from poultry to wild waterbirds. Several sites involved in host specificity (addition of a glycosylation site in the receptor binding domain) and virulence (insertion of amino acids in the cleavage site) were found to be positively selected in HA nucleotide sequences, in genetically unrelated lineages, suggesting parallel evolution for the HA gene of IA viruses in domestic birds. These results highlight that evolutionary consequences of bird host shifts would need to be further studied to understand the ecological and molecular mechanisms involved in the emergence of domestic bird-adapted viruses.

PhyloMap: an algorithm for visualizing relationships of large sequence data sets and its application to the influenza A virus genome

Background

Results of phylogenetic analysis are often visualized as phylogenetic trees. Such a tree can typically only include up to a few hundred sequences. When more than a few thousand sequences are to be included, analyzing the phylogenetic relationships among them becomes a challenging task. The recent frequent outbreaks of influenza A viruses have resulted in the rapid accumulation of corresponding genome sequences. Currently, there are more than 7500 influenza A virus genomes in the database. There are no efficient ways of representing this huge data set as a whole, thus preventing a further understanding of the diversity of the influenza A virus genome.

Results

Here we present a new algorithm, "PhyloMap", which combines ordination, vector quantization, and phylogenetic tree construction to give an elegant representation of a large sequence data set. The use of PhyloMap on influenza A virus genome sequences reveals the phylogenetic relationships of the internal genes that cannot be seen when only a subset of sequences are analyzed.

Conclusions

The application of PhyloMap to influenza A virus genome data shows that it is a robust algorithm for analyzing large sequence data sets. It utilizes the entire data set, minimizes bias, and provides intuitive visualization. PhyloMap is implemented in JAVA, and the source code is freely available at http://www.biochem.uni-luebeck.de/public/software/phylomap.html

Extensive geographic mosaicism in avian influenza viruses from gulls in the northern hemisphere

Due to limited interaction of migratory birds between Eurasia and America, two independent avian influenza virus (AIV) gene pools have evolved. There is evidence of low frequency reassortment between these regions, which has major implications in global AIV dynamics. Indeed, all currently circulating lineages of the PB1 and PA segments in North America are of Eurasian origin. Large-scale analyses of intercontinental reassortment have shown that viruses isolated from Charadriiformes (gulls, terns, and shorebirds) are the major contributor of these outsider events. To clarify the role of gulls in AIV dynamics, specifically in movement of genes between geographic regions, we have sequenced six gull AIV isolated in Alaska and analyzed these along with 142 other available gull virus sequences. Basic investigations of host species and the locations and times of isolation reveal biases in the available sequence information. Despite these biases, our analyses reveal a high frequency of geographic reassortment in gull viruses isolated in America. This intercontinental gene mixing is not found in the viruses isolated from gulls in Eurasia. This study demonstrates that gulls are important as vectors for geographically reassorted viruses, particularly in America, and that more surveillance effort should be placed on this group of birds.

Analysis of antigenically important residues in human influenza A virus in terms of B-cell epitopes

In this paper we undertake an analysis of the antigenicity of influenza A virus hemagglutinin. We developed a novel computational approach to the identification of antigenically active regions and showed that the amino acid substitutions between successive predominant seasonal strains form clusters that are consistent, in terms of both their location and their size, with the properties of B-cell epitopes in general and with those epitopes that have been identified experimentally in influenza A virus hemagglutinin to date. Such an interpretation provides a biologically plausible framework for an understanding of the location of antigenically important substitutions that is more specific than the canonical "antigenic site" model and provides an effective basis for deriving models that predict antigenic escape in the H3N2 subtype. Our results support recent indications that antibodies binding to the "stalk" region of hemagglutinin are found in the human population and exert evolutionary pressure on the virus. Our computational approach provides a possible method for identifying antigenic escape through evolution in this region, which in some cases will not be identified by the hemagglutinin inhibition assay.

Spatial dynamics of human-origin H1 influenza A virus in North American swine

The emergence and rapid global spread of the swine-origin H1N1/09 pandemic influenza A virus in humans underscores the importance of swine populations as reservoirs for genetically diverse influenza viruses with the potential to infect humans. However, despite their significance for animal and human health, relatively little is known about the phylogeography of swine influenza viruses in the United States. This study utilizes an expansive data set of hemagglutinin (HA1) sequences (n = 1516) from swine influenza viruses collected in North America during the period 2003–2010. With these data we investigate the spatial dissemination of a novel influenza virus of the H1 subtype that was introduced into the North American swine population via two separate human-to-swine transmission events around 2003. Bayesian phylogeographic analysis reveals that the spatial dissemination of this influenza virus in the US swine population follows long-distance swine movements from the Southern US to the Midwest, a corn-rich commercial center that imports millions of swine annually. Hence, multiple genetically diverse influenza viruses are introduced and co-circulate in the Midwest, providing the opportunity for genomic reassortment. Overall, the Midwest serves primarily as an ecological sink for swine influenza in the US, with sources of virus genetic diversity instead located in the Southeast (mainly North Carolina) and South-central (mainly Oklahoma) regions. Understanding the importance of long-distance pig transportation in the evolution and spatial dissemination of the influenza virus in swine may inform future strategies for the surveillance and control of influenza, and perhaps other swine pathogens.

Since 1998, genetically and antigenically diverse influenza A viruses have circulated in North American swine due to continuous cross-species transmission and reassortment with avian and human influenza viruses, presenting a pandemic threat to humans. Millions of swine are transported year-round from the southern United States into the corn-rich Midwest, but the importance of these movements in the spatial dissemination and evolution of the influenza virus in swine is unknown. Using a large data set of influenza virus sequences collected in North American swine during 2003–2010, we investigated the spatial dynamics of two influenza viruses of the H1 subtype that were introduced into swine from humans around 2003. Employing recently developed Bayesian phylogeography methods, we find that the spread of this influenza virus follows the large-scale transport of swine from the South to the Midwest. Based on this pattern of viral migration, we suggest that the genetic diversity of swine influenza viruses in the Midwest is continually augmented by the importation of viruses from source populations located in the South. Understanding the importance of long-distance pig movements in the evolution and spatial dissemination of influenza virus in swine may inform future strategies for the surveillance and control of influenza, and perhaps other swine pathogens.

Identification of potential conserved RNA secondary structure throughout influenza A coding regions

Influenza A is a negative sense RNA virus of significant public health concern. While much is understood about the life cycle of the virus, knowledge of RNA secondary structure in influenza A virus is sparse. Predictions of RNA secondary structure can focus experimental efforts. The present study analyzes coding regions of the eight viral genome segments in both the (+) and (-) sense RNA for conserved secondary structure. The predictions are based on identifying regions of unusual thermodynamic stabilities and are correlated with studies of suppression of synonymous codon usage (SSCU). The results indicate that secondary structure is favored in the (+) sense influenza RNA. Twenty regions with putative conserved RNA structure have been identified, including two previously described structured regions. Of these predictions, eight have high thermodynamic stability and SSCU, with five of these corresponding to current annotations (e.g., splice sites), while the remaining 12 are predicted by the thermodynamics alone. Secondary structures with high conservation of base-pairing are proposed within the five regions having known function. A combination of thermodynamics, amino acid and nucleotide sequence comparisons along with SSCU was essential for revealing potential secondary structures.

Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG)

In April 2008, a nucleotide-sequence-based, complete genome classification system was developed for group A rotaviruses (RVs). This system assigns a specific genotype to each of the 11 genome segments of a particular RV strain according to established nucleotide percent cutoff values. Using this approach, the genome of individual RV strains are given the complete descriptor of Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx. The Rotavirus Classification Working Group (RCWG) was formed by scientists in the field to maintain, evaluate and develop the RV genotype classification system, in particular to aid in the designation of new genotypes. Since its conception, the group has ratified 51 new genotypes: as of April 2011, new genotypes for VP7 (G20-G27), VP4 (P[28]-P[35]), VP6 (I12-I16), VP1 (R5-R9), VP2 (C6-C9), VP3 (M7-M8), NSP1 (A15-A16), NSP2 (N6-N9), NSP3 (T8-T12), NSP4 (E12-E14) and NSP5/6 (H7-H11) have been defined for RV strains recovered from humans, cows, pigs, horses, mice, South American camelids (guanaco), chickens, turkeys, pheasants, bats and a sugar glider. With increasing numbers of complete RV genome sequences becoming available, a standardized RV strain nomenclature system is needed, and the RCWG proposes that individual RV strains are named as follows: RV group/species of origin/country of identification/common name/year of identification/G- and P-type. In collaboration with the National Center for Biotechnology Information (NCBI), the RCWG is also working on developing a RV-specific resource for the deposition of nucleotide sequences. This resource will provide useful information regarding RV strains, including, but not limited to, the individual gene genotypes and epidemiological and clinical information. Together, the proposed nomenclature system and the NCBI RV resource will offer highly useful tools for investigators to search for, retrieve, and analyze the ever-growing volume of RV genomic data.

Adamantane- and oseltamivir-resistant seasonal A (H1N1) and pandemic A (H1N1) 2009 influenza viruses in Guangdong, China, during 2008 and 2009

Adamantane and oseltamivir resistance among influenza viruses is a major concern to public health officials. To determine the prevalence of antiviral-resistant influenza viruses in Guangdong, China, 244 seasonal A (H1N1) and 222 pandemic A (H1N1) 2009 viruses were screened for oseltamivir resistance by a fluorescence-based neuraminidase (NA) inhibition assay along with NA gene sequencing. Also, 147 seasonal A (H1N1) viruses were sequenced to detect adamantane resistance markers in M2. Adamantane-resistant seasonal A (H1N1) viruses clustering to clade 2C were dominant in 2008, followed by oseltamivir-resistant seasonal A (H1N1) viruses, clustering to clade 2B during January and May 2009. In June 2009, a lineage of double-resistant seasonal A (H1N1) viruses emerged, until it was replaced by the pandemic A (H1N1) 2009 viruses. The lineage most likely resulted from reassortment under the pressure of the overuse of adamantanes. As all viruses were resistant to at least one of the two types of antiviral agents, the need for close monitoring of the prevalence of antiviral resistance is stressed.

Genetic composition of contemporary swine influenza viruses in the West Central region of the United States of America

BACKGROUND

Because of continuous circulation in different animal species and humans, influenza viruses have host-specific phenotypic and genetic features. Reassortment of the genome segments can significantly change virus phenotype, potentially generating virus with pandemic potential. In 2009, a new pandemic influenza virus emerged.

OBJECTIVES

In this study, we attempted to find precursor viruses or genes of pandemic H1N1 influenza 2009 among 25 swine influenza viruses, isolated in the West Central region of the United States of America (USA), between 2007 and 2009. The Phylogenetically Similar Triple-Reassortant Internal Genes (PSTRIG) cassette of all the viruses studied here as well as the PSTRIG cassette of pandemic H1N1 viruses have close but equidistant phylogenetic relationships to the early triple-reassortant swine H3N2 influenza A isolated in the USA in 1998.

METHODS

Samples (nasal swabs and lung tissue lavage) were taken from swine with or without clinical signs of respiratory disease via farmer-funded syndromic surveillance. All studied viruses were isolated in Madin-Darby Canine Kidney cell cultures from the above-mentioned samples according to standard protocols recommended for influenza virus isolation. Sequences were obtained using BigDye Terminator v3.1 Cycle Sequencing kit. Phylogenetic trees were built with MEGA 4.0 software using maximum composite likelihood algorithm and neighbor-joining method for tree topology reconstruction.

RESULTS

Among the 25 viruses studied, we have not found any gene segments of Eurasian origin. Our results suggest that pandemic H1N1 viruses diverged and are not directly descended from swine viruses that have been circulating in USA since 1998.

Mixed influenza A and B infections complicate the detection of influenza viruses with altered sensitivities to neuraminidase inhibitors

Previously, three influenza A(H3N2) isolates with a reduced susceptibility to the neuraminidase inhibitors (NAIs) zanamivir and oseltamivir were identified during screening by the Neuraminidase Inhibitor Susceptibility Network (NISN). The isolates were from untreated patients from the first three years post licensure of the NAIs. We plaque-purified progeny from each of these isolates and determined the NAI sensitivity of each plaqued population. Sequencing and serology for each population revealed that the isolates contained a mix of wild type influenza A(H3N2) and influenza B. The NAI susceptibility reductions that had originally been reported were a consequence of influenza B neuraminidases that have lower relative NAI sensitivities, rather than being due to resistant influenza A(H3N2) viruses. Our study highlights the need to check for mixed influenza infections when isolates with potentially lower sensitivities to NAIs are identified.

Evolution of the influenza A virus untranslated regions

In case of the influenza virus, untranslated regions are considered to possess plural functions. The specific packaging model suggests that each of eight segments contains a unique "packaging signal", situated within the untranslated region and coding regions at the 5' and 3' ends of the genomic RNA. In addition to packaging signal, UTRs play key roles in polyadenylation and signals for genomic replication. In the present study, we report the evolutionary characteristics of untranslated regions of influenza virus. We analyzed 574 sequences of the human influenza A (H3N2) virus. The sequence of the untranslated region is highly conserved. Our analysis produced several observations regarding the untranslated region: (1) it has an extremely low divergence; (2) its evolutionary speed is not associated with the protein(s) it encodes; (3) its evolutionary rate is smaller than that of the open reading frame. We revealed the evolutionary characteristics of the untranslated region; this is the first study on the non-coding region of the influenza virus from the perspective of evolutionary molecular biology. In the future, clarifying the detailed mechanisms of packaging as well as the function of the RNA sequence in the untranslated region will increase the understanding of viral life cycles and evolution.

Genetic evidence for containment of viruses in the first outbreak of influenza A pandemic (H1N1) 2009 in Kobe, Japan

Please cite this paper as: Ninomiya‐Mori et al. (2010) Genetic evidence for containment of viruses in the first outbreak of influenza A pandemic (H1N1) 2009 in Kobe, Japan. Influenza and Other Respiratory Viruses Doi: 10.1111/j.1750‐2659.2010.00188.x.

Background  On 16 May 2009, a high school student in Kobe with no history of overseas travel was reported as the first case of influenza A pandemic (H1N1) 2009 virus infection in Japan. Subsequently, it was revealed that the infection had spread to some cities in the Kansai region and most patients were high school students. The number of patients decreased rapidly within a week; however, it began to increase in the middle of July.

Methods  We phylogenetically analyzed viral characteristics using 27 viruses isolated from patients living in Kobe.

Results and conclusions  We demonstrated that viruses isolated in the early phase of the outbreak were distinguishable from those after the reappearance of patients. These findings provide genetic evidence for the effectiveness of public health containment measures in the Kansai region in preventing the progression of the outbreak.

The nucleoprotein and matrix protein segments of H5N1 influenza viruses are responsible for dominance in embryonated eggs

Since their emergence in 1996 in southern China, highly pathogenic H5N1 avian influenza viruses have spread widely and continue to circulate in some countries. Genetic reassortment has created multiple H5N1 virus lineages, some of which are dominant in nature. However, the mechanism by which certain H5N1 influenza virus lineages (or genotypes) become dominant in avian species remains unknown. Here, we used competitive inoculation and genetic analysis of the resultant viruses to show that the nucleoprotein (NP) and matrix protein (M) segments of Fujian-like viruses (clade 2.3.4), which became predominant in southern China in mid-2006, are responsible for viral dominance in embryonated eggs. We further found that specific residues in the NP and M proteins play key roles in conferring this viral dominance; specifically, a glutamic acid at position 66 in M2 was conserved among the Fujian-like viruses. These results suggest roles for these viral proteins in influenza virus dominance.

Reassortment and mutations associated with emergence and spread of oseltamivir-resistant seasonal influenza A/H1N1 viruses in 2005-2009

A dramatic increase in the frequency of the H275Y mutation in the neuraminidase (NA), conferring resistance to oseltamivir, has been detected in human seasonal influenza A/H1N1 viruses since the influenza season of 2007-2008. The resistant viruses emerged in the ratio of 14.3% and quickly reached 100% in Taiwan from September to December 2008. To explore the mechanisms responsible for emergence and spread of the resistant viruses, we analyzed the complete genome sequences of 25 viruses collected during 2005-2009 in Taiwan, which were chosen from various clade viruses, 1, 2A, 2B-1, 2B-2, 2C-1 and 2C-2 by the classification of hemagglutinin (HA) sequences. Our data revealed that the dominant variant, clade 2B-1, in the 2007-2008 influenza emerged through an intra-subtype 4+4 reassortment between clade 1 and 2 viruses. The dominant variant acquired additional substitutions, including A206T in HA, H275Y and D354G in NA, L30R and H41P in PB1-F2, and V411I and P453S in basic polymerase 2 (PB2) proteins and subsequently caused the 2008-2009 influenza epidemic in Taiwan, accompanying the widespread oseltamivir-resistant viruses. We also characterized another 3+5 reassortant virus which became double resistant to oseltamivir and amantadine. Comparison of oseltamivir-resistant influenza A/H1N1 viruses belonging to various clades in our study highlighted that both reassortment and mutations were associated with emergence and spread of these viruses and the specific mutation, H275Y, conferring to antiviral resistance, was acquired in a hitch-hiking mechanism during the viral evolutionary processes.

A comprehensive benchmark study of multiple sequence alignment methods: current challenges and future perspectives

Multiple comparison or alignmentof protein sequences has become a fundamental tool in many different domains in modern molecular biology, from evolutionary studies to prediction of 2D/3D structure, molecular function and inter-molecular interactions etc. By placing the sequence in the framework of the overall family, multiple alignments can be used to identify conserved features and to highlight differences or specificities. In this paper, we describe a comprehensive evaluation of many of the most popular methods for multiple sequence alignment (MSA), based on a new benchmark test set. The benchmark is designed to represent typical problems encountered when aligning the large protein sequence sets that result from today's high throughput biotechnologies. We show that alignmentmethods have significantly progressed and can now identify most of the shared sequence features that determine the broad molecular function(s) of a protein family, even for divergent sequences. However,we have identified a number of important challenges. First, the locally conserved regions, that reflect functional specificities or that modulate a protein's function in a given cellular context,are less well aligned. Second, motifs in natively disordered regions are often misaligned. Third, the badly predicted or fragmentary protein sequences, which make up a large proportion of today's databases, lead to a significant number of alignment errors. Based on this study, we demonstrate that the existing MSA methods can be exploited in combination to improve alignment accuracy, although novel approaches will still be needed to fully explore the most difficult regions. We then propose knowledge-enabled, dynamic solutions that will hopefully pave the way to enhanced alignment construction and exploitation in future evolutionary systems biology studies.

Prediction of directional changes of influenza A virus genome sequences with emphasis on pandemic H1N1/09 as a model case

Influenza virus poses a significant threat to public health, as exemplified by the recent introduction of the new pandemic strain H1N1/09 into human populations. Pandemics have been initiated by the occurrence of novel changes in animal sources that eventually adapt to human. One important issue in studies of viral genomes, particularly those of influenza virus, is to predict possible changes in genomic sequence that will become hazardous. We previously established a clustering method termed ‘BLSOM’ (batch-learning self-organizing map) that does not depend on sequence alignment and can characterize and compare even 1 million genomic sequences in one run. Strategies for comparing a vast number of genomic sequences simultaneously become increasingly important in genome studies because of remarkable progresses in nucleotide sequencing. In this study, we have constructed BLSOMs based on the oligonucleotide and codon composition of all influenza A viral strains available. Without prior information with regard to their hosts, sequences derived from strains isolated from avian or human sources were successfully clustered according to the hosts. Notably, the pandemic H1N1/09 strains have oligonucleotide and codon compositions that are clearly different from those of human seasonal influenza A strains. This enables us to infer future directional changes in the influenza A viral genome.

Genome-scale evolution and phylodynamics of equine H3N8 influenza A virus

Equine influenza viruses (EIVs) of the H3N8 and H7N7 subtypes are the causative agents of an important disease of horses. While EIV H7N7 apparently is extinct, H3N8 viruses have circulated for more than 50 years. Like human influenza viruses, EIV H3N8 caused a transcontinental pandemic followed by further outbreaks and epidemics, even in populations with high vaccination coverage. Recently, EIV H3N8 jumped the species barrier to infect dogs. Despite its importance as an agent of infectious disease, the mechanisms that underpin the evolutionary and epidemiological dynamics of EIV are poorly understood, particularly at a genomic scale. To determine the evolutionary history and phylodynamics of EIV H3N8, we conducted an extensive analysis of 82 complete viral genomes sampled during a 45-year span. We show that both intra- and intersubtype reassortment have played a major role in the evolution of EIV, and we suggest that intrasubtype reassortment resulted in enhanced virulence while heterosubtypic reassortment contributed to the extinction of EIV H7N7. We also show that EIV evolves at a slower rate than other influenza viruses, even though it seems to be subject to similar immune selection pressures. However, a relatively high rate of amino acid replacement is observed in the polymerase acidic (PA) segment, with some evidence for adaptive evolution. Most notably, an analysis of viral population dynamics provided evidence for a major population bottleneck of EIV H3N8 during the 1980s, which we suggest resulted from changes in herd immunity due to an increase in vaccination coverage.

Interspecific exchange of avian influenza virus genes in Alaska: the influence of trans-hemispheric migratory tendency and breeding ground sympatry

The movement and transmission of avian influenza viral strains via wild migratory birds may vary by host species as a result of migratory tendency and sympatry with other infected individuals. To examine the roles of host migratory tendency and species sympatry on the movement of Eurasian low-pathogenic avian influenza (LPAI) genes into North America, we characterized migratory patterns and LPAI viral genomic variation in mallards (Anas platyrhynchos) of Alaska in comparison with LPAI diversity of northern pintails (Anas acuta). A 50-year band-recovery data set suggests that unlike northern pintails, mallards rarely make trans-hemispheric migrations between Alaska and Eurasia. Concordantly, fewer (14.5%) of 62 LPAI isolates from mallards contained Eurasian gene segments compared to those from 97 northern pintails (35%), a species with greater inter-continental migratory tendency. Aerial survey and banding data suggest that mallards and northern pintails are largely sympatric throughout Alaska during the breeding season, promoting opportunities for interspecific transmission. Comparisons of full-genome isolates confirmed near-complete genetic homology (>99.5%) of seven viruses between mallards and northern pintails. This study found viral segments of Eurasian lineage at a higher frequency in mallards than previous studies, suggesting transmission from other avian species migrating inter-hemispherically or the common occurrence of endemic Alaskan viruses containing segments of Eurasian origin. We conclude that mallards are unlikely to transfer Asian-origin viruses directly to North America via Alaska but that they are likely infected with Asian-origin viruses via interspecific transfer from species with regular migrations to the Eastern Hemisphere.

Prevalence of epistasis in the evolution of influenza A surface proteins

The surface proteins of human influenza A viruses experience positive selection to escape both human immunity and, more recently, antiviral drug treatments. In bacteria and viruses, immune-escape and drug-resistant phenotypes often appear through a combination of several mutations that have epistatic effects on pathogen fitness. However, the extent and structure of epistasis in influenza viral proteins have not been systematically investigated. Here, we develop a novel statistical method to detect positive epistasis between pairs of sites in a protein, based on the observed temporal patterns of sequence evolution. The method rests on the simple idea that a substitution at one site should rapidly follow a substitution at another site if the sites are positively epistatic. We apply this method to the surface proteins hemagglutinin and neuraminidase of influenza A virus subtypes H3N2 and H1N1. Compared to a non-epistatic null distribution, we detect substantial amounts of epistasis and determine the identities of putatively epistatic pairs of sites. In particular, using sequence data alone, our method identifies epistatic interactions between specific sites in neuraminidase that have recently been demonstrated, in vitro, to confer resistance to the drug oseltamivir; these epistatic interactions are responsible for widespread drug resistance among H1N1 viruses circulating today. This experimental validation demonstrates the predictive power of our method to identify epistatic sites of importance for viral adaptation and public health. We conclude that epistasis plays a large role in shaping the molecular evolution of influenza viruses. In particular, sites with , which would normally not be identified as positively selected, can facilitate viral adaptation through epistatic interactions with their partner sites. The knowledge of specific interactions among sites in influenza proteins may help us to predict the course of antigenic evolution and, consequently, to select more appropriate vaccines and drugs.

Epistasis describes non-additive interactions among genetic sites: the consequence of a mutation at one site may depend on the status of the genome at other sites. In an extreme case, a mutation may have no effect if it arises on one genetic background, but a strong effect on another background. Epistatic mutations in viruses and bacteria that live under severe conditions, such as antibiotic treatments or immune pressure, often allow pathogens to develop drug resistance or escape the immune system. In this paper we develop a new phylogenetic method for detecting epistasis, and we apply this method to the surface proteins of the influenza A virus, which are important targets of the immune system and drug treatments. The authors identify and characterize hundreds of epistatic mutations in these proteins. Among those identified, we find the specific epistatic mutations that were recently shown, experimentally, to confer resistance to the drug Tamiflu. The results of this study may help to predict the course of influenza's antigenic evolution and to select more appropriate vaccines and drugs.

Changed epitopes drive the antigenic drift for influenza A (H3N2) viruses

Background

In circulating influenza viruses, gradually accumulated mutations on the glycoprotein hemagglutinin (HA), which interacts with infectivity-neutralizing antibodies, lead to the escape of immune system (called antigenic drift). The antibody recognition is highly correlated to the conformation change on the antigenic sites (epitopes), which locate on HA surface. To quantify a changed epitope for escaping from neutralizing antibodies is the basis for the antigenic drift and vaccine development.

Results

We have developed an epitope-based method to identify the antigenic drift of influenza A utilizing the conformation changes on epitopes. A changed epitope, an antigenic site on HA with an accumulated conformation change to escape from neutralizing antibody, can be considered as a "key feature" for representing the antigenic drift. According to hemagglutination inhibition (HI) assays and HA/antibody complex structures, we statistically measured the conformation change of an epitope by considering the number of critical position mutations with high genetic diversity and antigenic scores. Experimental results show that two critical position mutations can induce the conformation change of an epitope to escape from the antibody recognition. Among five epitopes of HA, epitopes A and B, which are near to the receptor binding site, play a key role for neutralizing antibodies. In addition, two changed epitopes often drive the antigenic drift and can explain the selections of 24 WHO vaccine strains.

Conclusions

Our method is able to quantify the changed epitopes on HA for predicting the antigenic variants and providing biological insights to the vaccine updates. We believe that our method is robust and useful for studying influenza virus evolution and vaccine development.

A computationally optimized broadly reactive antigen (COBRA) based H5N1 VLP vaccine elicits broadly reactive antibodies in mice and ferrets

Pandemic outbreaks of influenza are caused by the emergence of a pathogenic and transmissible virus to which the human population is immunologically naïve. Recent outbreaks of highly pathogenic avian influenza (HPAI) of the H5N1 subtype are of particular concern because of the high mortality rate (60% case fatality rate) and novel subtype. In order to develop a vaccine that elicits broadly reactive antibody responses against emerging H5N1 isolates, we utilized a novel antigen design technique termed computationally optimized broadly reactive antigen (COBRA). The COBRA HA sequence was based upon HA amino acid sequences from clade 2 H5N1 human infections and the expressed protein retained the ability to bind the receptor, as well as mediate particle fusion. Non-infectious recombinant VLP vaccines using the COBRA HA were purified from a mammalian expression system. Mice and ferrets vaccinated with COBRA HA H5N1 VLPs had protective levels of HAI antibodies to a representative isolates from each subclade of clade 2. Furthermore, VLP vaccinated animals were completely protected from a lethal challenge of the clade 2.2 H5N1 virus A/Whooper Swan/Mongolia/244/2005. This is the first report describing the use of COBRA-based antigen design. The COBRA HA H5N1 VLP vaccine elicited broadly reactive antibodies and is an effective influenza vaccine against HPAI virus.

Proteotyping to establish gene origin within reassortant influenza viruses

The application of a rapid and direct proteotyping approach with which to identify the gene origin of viral antigens in a reassortant influenza strain is demonstrated. The reassortant strain, constructed for a vaccine against type A 2009 H1N1 pandemic influenza, contains genes derived from a wild-type pandemic strain (A/California/7/2009) and an egg adapted high-growth strain (denoted NYMC X-157) derived from an earlier A/Puerto Rico/8/34 strain. The proteotyping approach employs modern proteomics methods and high resolution mass spectrometry to correctly establish that the genes of the surface antigens, hemagglutinin and neuraminidase, are derived from the A/California/7/2009 strain while those for nucleoprotein and matrix protein M1 antigens are derived from the NYMC X-157 strain. This is achieved for both gel-separated antigens and those from a whole vaccine digest. Furthermore, signature peptides detected in the mass spectra of the digested antigens enable the engineered reassortant strain to be identified as a type A virus of the H1N1 subtype in accord with earlier studies. The results demonstrate that proteotyping approach provides a more direct and rapid approach over RT-PCR with which to characterize reassortant strains of the influenza virus at the molecular protein level. Given that these strains pose the greatest risk to human and animal health and have been responsible for all human pandemics of the 20th and 21st centuries, there is a vital need for the origins and evolutionary history of these strains to be rapidly established.

Influenza genome diversity and evolution

The influenza viruses contain highly variable genomes and are able to infect a wide range of host species. Large-scale sequencing projects have collected abundant influenza sequence data for assessing influenza genome diversity and evolution. This work reviews current influenza sequence databases characteristics and statistics, as well as recent studies utilizing these databases to unravel influenza virus diversity and evolution. Also discussed are the newest deep sequencing methods and their applications to influenza virus research.

Identification of novel conserved functional motifs across most Influenza A viral strains

Background

Influenza A virus poses a continuous threat to global public health. Design of novel universal drugs and vaccine requires a careful analysis of different strains of Influenza A viral genome from diverse hosts and subtypes. We performed a systematic in silico analysis of Influenza A viral segments of all available Influenza A viral strains and subtypes and grouped them based on host, subtype, and years isolated, and through multiple sequence alignments we extrapolated conserved regions, motifs, and accessible regions for functional mapping and annotation.

Results

Across all species and strains 87 highly conserved regions (conservation percentage > = 90%) and 19 functional motifs (conservation percentage = 100%) were found in PB2, PB1, PA, NP, M, and NS segments. The conservation percentage of these segments ranged between 94 - 98% in human strains (the most conserved), 85 - 93% in swine strains (the most variable), and 91 - 94% in avian strains. The most conserved segment was different in each host (PB1 for human strains, NS for avian strains, and M for swine strains). Target accessibility prediction yielded 324 accessible regions, with a single stranded probability > 0.5, of which 78 coincided with conserved regions. Some of the interesting annotations in these regions included sites for protein-protein interactions, the RNA binding groove, and the proton ion channel.

Conclusions

The influenza virus has evolved to adapt to its host through variations in the GC content and conservation percentage of the conserved regions. Nineteen universal conserved functional motifs were discovered, of which some were accessible regions with interesting biological functions. These regions will serve as a foundation for universal drug targets as well as universal vaccine design.

BiologicalNetworks--tools enabling the integration of multi-scale data for the host-pathogen studies

BACKGROUND

Understanding of immune response mechanisms of pathogen-infected host requires multi-scale analysis of genome-wide data. Data integration methods have proved useful to the study of biological processes in model organisms, but their systematic application to the study of host immune system response to a pathogen and human disease is still in the initial stage.

RESULTS

To study host-pathogen interaction on the systems biology level, an extension to the previously described BiologicalNetworks system is proposed. The developed methods and data integration and querying tools allow simplifying and streamlining the process of integration of diverse experimental data types, including molecular interactions and phylogenetic classifications, genomic sequences and protein structure information, gene expression and virulence data for pathogen-related studies. The data can be integrated from the databases and user's files for both public and private use.

CONCLUSIONS

The developed system can be used for the systems-level analysis of host-pathogen interactions, including host molecular pathways that are induced/repressed during the infections, co-expressed genes, and conserved transcription factor binding sites. Previously unknown to be associated with the influenza infection genes were identified and suggested for further investigation as potential drug targets. Developed methods and data are available through the Java application (from BiologicalNetworks program at http://www.biologicalnetworks.org) and web interface (at http://flu.sdsc.edu).

Effects of windowing and zero-padding on Complex Resonant Recognition Model for protein sequence analysis

Signal processing techniques such as Fourier Transform have widely been studied and successfully applied in many different areas. Techniques such as zero-padding and windowing have been developed and found very useful to improve the outcome of the signal processing methods. Resonant Recognition Model (RRM) and Complex Resonant Recognition Model (CRRM) that are based on the discrete Fourier Transform and widely used for the analysis of protein sequences do not consider such methods, which can however improve or alter the features extracted from the protein sequences. Therefore, in this paper, an extensive analysis was carried out to investigate into the influence of the zero-padding and windowing on the features extracted from the Complex Resonant Recognition Model. In order to present such effects, five different classes of influenza A virus Neuraminidase genes, which include H1N1, H1N2, H2N2, H3N2 and H5N1 genes, were used as a case study. For each of the Influenza A subtypes, two sets of Common Frequency Peaks (CFP) were extracted, one where windowing is applied and the other one where windowing is suppressed, for each signal length set for the analysis. In order to make all the signals (protein sequence) the same length, zero-padding was used. The signal lengths used in this study are set to 470, which is the maximum protein length, and also 512, 1024, 2048, 4096, 8192 and 16384 for further analysis. The results suggest that the windowing and zero-padding have key impact on CFP extracted from the Influenza A subtypes as the best match with CFP extracted from influenza A subtypes using CRRM is when the signal length of 4096 and windowing were both applied. Therefore, the outcome of this study should be taken into consideration for more accurate and reliable analysis of the protein sequences.

Infectious Disease Genomics

The history and development of infectious disease genomics are discussed in this chapter. HGP must not be restricted to the human genome and should include model organisms including mouse, bacteria, yeast, fruit fly, and worm. The completed or ongoing genome projects will provide enormous opportunities for the discovery of novel vaccines and drug targets against human pathogens as well as the improvement of diagnosis and discovery of infectious agents and the development of new strategies for invertebrate vector control. The polysaccharide capsule is important for meningococci to escape from complement-mediated killing. With the completion of the genome sequence of a virulent MenB strain, a “reverse vaccinology” approach was applied for the development of a universal MenB vaccine by Novartis. The indispensable fatty acid synthase (FAS) pathway in bacteria has been regarded as a promising target for the development of antimicrobial agents. Through a systematic screening of 250,000 natural product extracts, a Merck team identified a potent and broad-spectrum antibiotic, platensimycin, which is derived from Streptomyces platensis. Vector Biology Network was formed to achieve three goals (1) to develop basic tools for the stable transformation of anopheline mosquitoes by the year 2000; (2) to engineer a mosquito incapable of carrying the malaria parasite by 2005; and (3) to run controlled experiments to test how to drive the engineered genotype into wild mosquito populations by 2010. The most immediate impact of a completely sequenced pathogen genome is for infectious disease diagnosis.

Phylogeography of the spring and fall waves of the H1N1/09 pandemic influenza virus in the United States

Spatial variation in the epidemiological patterns of successive waves of pandemic influenza virus in humans has been documented throughout the 20th century but never understood at a molecular level. However, the unprecedented intensity of sampling and whole-genome sequencing of the H1N1/09 pandemic virus now makes such an approach possible. To determine whether the spring and fall waves of the H1N1/09 influenza pandemic were associated with different epidemiological patterns, we undertook a large-scale phylogeographic analysis of viruses sampled from three localities in the United States. Analysis of genomic and epidemiological data reveals distinct spatial heterogeneities associated with the first pandemic wave, March to July 2009, in Houston, TX, Milwaukee, WI, and New York State. In Houston, no specific H1N1/09 viral lineage dominated during the spring of 2009, a period when little epidemiological activity was observed in Texas. In contrast, major pandemic outbreaks occurred at this time in Milwaukee and New York State, each dominated by a different viral lineage and resulting from strong founder effects. During the second pandemic wave, beginning in August 2009, all three U.S. localities were dominated by a single viral lineage, that which had been dominant in New York during wave 1. Hence, during this second phase of the pandemic, extensive viral migration and mixing diffused the spatially defined population structure that had characterized wave 1, amplifying the one viral lineage that had dominated early on in one of the world's largest international travel centers.