Antigenic drift in H5N1 avian influenza virus in poultry is driven by mutations in major antigenic sites of the hemagglutinin molecule analogous to those for human influenza virus

H5N1 highly pathogenic avian influenza virus has been endemic in poultry in Egypt since 2008, notwithstanding the implementation of mass vaccination and culling of infected birds. Extensive circulation of the virus has resulted in a progressive genetic evolution and an antigenic drift. In poultry, the occurrence of antigenic drift in avian influenza viruses is less well documented and the mechanisms remain to be clarified. To test the hypothesis that H5N1 antigenic drift is driven by mechanisms similar to type A influenza viruses in humans, we generated reassortant viruses, by reverse genetics, that harbored molecular changes identified in genetically divergent viruses circulating in the vaccinated population. Parental and reassortant phenotype viruses were antigenically analyzed by hemagglutination inhibition (HI) test and microneutralization (MN) assay. The results of the study indicate that the antigenic drift of H5N1 in poultry is driven by multiple mutations primarily occurring in major antigenic sites at the receptor binding subdomain, similarly to what has been described for human influenza H1 and H3 subtype viruses.     

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

ABSTRACT: BACKGROUND: Influenza virus undergoes rapid evolution by both antigenic shift and antigenic drift. Antibodies, particularly those binding near the receptor-binding site of hemagglutinin (HA) or the neuraminidase (NA) active site, are thought to be the primary defense against influenza infection, and mutations in antibody binding sites can reduce or eliminate antibody binding. The binding of antibodies to their cognate antigens is governed by such biophysical properties of the interacting surfaces as shape, non-polar and polar surface area, and charge. Methods: To understand forces shaping evolution of influenza virus, we have examined HA sequences of human influenza A and B viruses, assigning each amino acid values reflecting total accessible surface area, non-polar and polar surface area, and net charge due to the side chain. Changes in each of these values between neighboring sequences were calculated for each residue and mapped onto the crystal structures. Results: Areas of HA showing the highest frequency of changes agreed well with previously identified antigenic sites in H3 and H1 HAs, and allowed us to propose more detailed antigenic maps and novel antigenic sites for H1 and influenza B HA. Changes in biophysical properties differed between HAs of different subtypes, and between different antigenic sites of the same HA. For H1, statistically significant differences in several biophysical quantities compared to residues lying outside antigenic sites were seen for some antigenic sites but not others. Influenza B antigenic sites all show statistically significant differences in biophysical quantities for all antigenic sites, whereas no statistically significant differences in biophysical quantities were seen for any antigenic site is seen for H3. In many cases, residues previously shown to be under positive selection at the genetic level also undergo rapid change in biophysical properties. Conclusions: The biophysical consequences of amino acid changes introduced by antigenic drift vary from subtype to subtype, and between different antigenic sites. This suggests that the significance of antibody binding in selecting new variants may also be variable for different antigenic sites and influenza subtypes.     

Phylodynamic Inference and Model Assessment with Approximate Bayesian Computation: Influenza as a Case Study

A key priority in infectious disease research is to understand the ecological and evolutionary drivers of viral diseases from data on disease incidence as well as viral genetic and antigenic variation. We propose using a simulation-based, Bayesian method known as Approximate Bayesian Computation (ABC) to fit and assess phylodynamic models that simulate pathogen evolution and ecology against summaries of these data. We illustrate the versatility of the method by analyzing two spatial models describing the phylodynamics of interpandemic human influenza virus subtype A(H3N2). The first model captures antigenic drift phenomenologically with continuously waning immunity, and the second epochal evolution model describes the replacement of major, relatively long-lived antigenic clusters. Combining features of long-term surveillance data from the Netherlands with features of influenza A (H3N2) hemagglutinin gene sequences sampled in northern Europe, key phylodynamic parameters can be estimated with ABC. Goodness-of-fit analyses reveal that the irregularity in interannual incidence and H3N2''s ladder-like hemagglutinin phylogeny are quantitatively only reproduced under the epochal evolution model within a spatial context. However, the concomitant incidence dynamics result in a very large reproductive number and are not consistent with empirical estimates of H3N2''s population level attack rate. These results demonstrate that the interactions between the evolutionary and ecological processes impose multiple quantitative constraints on the phylodynamic trajectories of influenza A(H3N2), so that sequence and surveillance data can be used synergistically. ABC, one of several data synthesis approaches, can easily interface a broad class of phylodynamic models with various types of data but requires careful calibration of the summaries and tolerance parameters.     

Evolution of the Hemagglutinin Protein of the New Pandemic H1N1 Influenza Virus: Maintaining Optimal Receptor Binding by Compensatory Substitutions

Pandemic influenza A H1N1 (pH1N1) virus emerged in 2009. In the subsequent 4 years, it acquired several genetic changes in its hemagglutinin (HA). Mutations may be expected while virus is adapting to the human host or upon evasion from adaptive immune responses. However, pH1N1 has not displayed any major antigenic changes so far. We examined the effect of the amino acid substitutions found to be most frequently occurring in the pH1N1 HA protein before 1 April 2012 on the receptor-binding properties of the virus by using recombinant soluble HA trimers. Two changes (S186P and S188T) were shown to increase the receptor-binding avidity of HA, whereas two others (A137T and A200T) decreased binding avidity. Construction of an HA protein tree revealed the worldwide emergence of several HA variants during the past few influenza seasons. Strikingly, two major variants harbor combinations of substitutions (S186P/A137T and S188T/A200T, respectively) with opposite individual effects on binding. Stepwise reconstruction of the HA proteins of these variants demonstrated that the mutations that increase receptor-binding avidity are compensated for by the acquisition of subsequent mutations. The combination of these substitutions restored the receptor-binding properties (avidity and specificity) of these HA variants to those of the parental virus. The results strongly suggest that the HA of pH1N1 was already optimally adapted to the human host upon its emergence in April 2009. Moreover, these results are in agreement with a recent model for antigenic drift, in which influenza A virus mutants with high and low receptor-binding avidity alternate.     

Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus

The number of N-linked glycosylation sites in the globular head of hemagglutinin (HA) has increased during evolution of H3N2 human influenza A virus. Here natural selection operating on the gains of N-linked glycosylation sites was examined by using the single-site analysis and the single-substitution analysis. In the single-site analysis, positive selection was not inferred at the amino acid sites where the substitutions generating N-linked glycosylation sites were observed, but was detected at antigenic sites. In contrast, in the single-substitution analysis, positive selection was detected for the amino acid substitutions generating N-linked glycosylation sites. The single-site analysis and the single-substitution analysis appeared to be suitable for detecting recurrent and episodic natural selection, respectively. The gains of N-linked glycosylation sites were likely to be positively selected for the function of shielding antigenic sites from immune responses. At the antigenic sites, positive selection appeared to have operated not only on the radical substitution but also on the conservative substitution in terms of the charge of amino acids, suggesting that the antigenic drift is not a by-product of the evolution of receptor binding avidity in HA of human H3N2 virus.     

Influenza A Gradual and Epochal Evolution: Insights from Simple Models

The recurrence of influenza A epidemics has originally been explained by a “continuous antigenic drift” scenario. Recently, it has been shown that if genetic drift is gradual, the evolution of influenza A main antigen, the haemagglutinin, is punctuated. As a consequence, it has been suggested that influenza A dynamics at the population level should be approximated by a serial model. Here, simple models are used to test whether a serial model requires gradual antigenic drift within groups of strains with the same antigenic properties (antigenic clusters). We compare the effect of status based and history based frameworks and the influence of reduced susceptibility and infectivity assumptions on the transient dynamics of antigenic clusters. Our results reveal that the replacement of a resident antigenic cluster by a mutant cluster, as observed in data, is reproduced only by the status based model integrating the reduced infectivity assumption. This combination of assumptions is useful to overcome the otherwise extremely high model dimensionality of models incorporating many strains, but relies on a biological hypothesis not obviously satisfied. Our findings finally suggest the dynamical importance of gradual antigenic drift even in the presence of punctuated immune escape. A more regular renewal of susceptible pool than the one implemented in a serial model should be part of a minimal theory for influenza at the population level.     

Patterns of predicted T-cell epitopes associated with antigenic drift in influenza H3N2 hemagglutinin

Antigenic drift allowing escape from neutralizing antibodies is an important feature of transmission and survival of influenza viruses in host populations. Antigenic drift has been studied in particular detail for influenza A H3N2 and well defined antigenic clusters of this virus documented. We examine how host immunogenetics contributes to determination of the antibody spectrum, and hence the immune pressure bringing about antigenic drift. Using uTOPE™ bioinformatics analysis of predicted MHC binding, based on amino acid physical property principal components, we examined the binding affinity of all 9-mer and 15-mer peptides within the hemagglutinin 1 (HA1) of 447 H3N2 virus isolates to 35 MHC-I and 14 MHC-II alleles. We provide a comprehensive map of predicted MHC-I and MHC-II binding affinity for a broad array of HLA alleles for the H3N2 influenza HA1 protein. Each HLA allele exhibited a characteristic predicted binding pattern. Cluster analysis for each HLA allele shows that patterns based on predicted MHC binding mirror those described based on antibody binding. A single amino acid mutation or position displacement can result in a marked difference in MHC binding and hence potential T-helper function. We assessed the impact of individual amino acid changes in HA1 sequences between 10 virus isolates from 1968–2002, representative of antigenic clusters, to understand the changes in MHC binding over time. Gain and loss of predicted high affinity MHC-II binding sites with cluster transitions were documented. Predicted high affinity MHC-II binding sites were adjacent to antibody binding sites. We conclude that host MHC diversity may have a major determinant role in the antigenic drift of influenza A H3N2.     

Paradox of peramorphic paedomorphosis: Heterochrony and human evolution

This paper reviews Gould's clock model for heterochronic processes and uses that model to develop simple matrix representations of growth and shape change. Matrix representations of growth and development provide a common formulation for all heterochronic processes. In particular, we show how neoteny can be diagnosed using such a matrix approach. The literature is rife with contradictory representations of how neoteny affects growth allometries and the timing of developmental events, and therefore of the role of neoteny in human evolution. Through the use of multivariate models, we explore these relationships and the internal consistency of opposing views. Gould's neoteny hypothesis for human evolution has been criticized for a number of reasons. Humans do not grow slowly. The slopes of our growth allometries show no common pattern of change vis-a-vis those of our closest relatives. Humans prolong rather than reduce rates of growth and development of body parts; the brain, for example, ceases growing later in humans than in apes, but during this prolonged period of early ontogeny, it grows at a rapid pace. This paper evaluates Gould's hypothesis and its critiques by focusing on particular questions. Does neoteny imply slow growth? Does it imply a unidirectional change in the rates of growth of traits? Under neoteny, should the brain cease growing in ancestor and descendant at the same age? Does prolongation of phases of growth and development confute neoteny? On the other hand, is paedomorphosis an inevitable consequence of prolonged growth and development? We show that, for all of these questions, the answer is no. © 1996 Wiley-Liss, Inc.     

Single Hemagglutinin Mutations That Alter both Antigenicity and Receptor Binding Avidity Influence Influenza Virus Antigenic Clustering

The hemagglutination inhibition (HAI) assay is the primary measurement used for identifying antigenically novel influenza virus strains. HAI assays measure the amount of reference sera required to prevent virus binding to red blood cells. Receptor binding avidities of viral strains are not usually taken into account when interpreting these assays. Here, we created antigenic maps of human H3N2 viruses that computationally account for variation in viral receptor binding avidities. These new antigenic maps differ qualitatively from conventional antigenic maps based on HAI measurements alone. We experimentally focused on an antigenic cluster associated with a single N145K hemagglutinin (HA) substitution that occurred between 1992 and 1995. Reverse-genetics experiments demonstrated that the N145K HA mutation increases viral receptor binding avidity. Enzyme-linked immunosorbent assays (ELISA) revealed that the N145K HA mutation does not prevent antibody binding; rather, viruses possessing this mutation escape antisera in HAI assays simply by attaching to cells more efficiently. Unexpectedly, we found an asymmetric antigenic effect of the N145K HA mutation. Once H3N2 viruses acquired K145, an epitope involving amino acid 145 became antigenically dominant. Antisera raised against an H3N2 strain possessing K145 had reduced reactivity to H3N2 strains possessing N145. Thus, individual mutations in HA can influence antigenic groupings of strains by altering receptor binding avidity and by changing the dominance of antibody responses. Our results indicate that it will be important to account for variation in viral receptor binding avidity when performing antigenic analyses in order to identify genuine antigenic differences among influenza virus variants.     

Antigenic and genetic evolution of swine influenza A (H3N2) viruses in Europe

In the early 1970s, a human influenza A/Port Chalmers/1/73 (H3N2)-like virus colonized the European swine population. Analyses of swine influenza A (H3N2) viruses isolated in The Netherlands and Belgium revealed that in the early 1990s, antigenic drift had occurred, away from A/Port Chalmers/1/73, the strain commonly used in influenza vaccines for pigs. Here we show that Italian swine influenza A (H3N2) viruses displayed antigenic and genetic changes similar to those observed in Northern European viruses in the same period. We used antigenic cartography methods for quantitative analyses of the antigenic evolution of European swine H3N2 viruses and observed a clustered virus evolution as seen for human viruses. Although the antigenic drift of swine and human H3N2 viruses has followed distinct evolutionary paths, potential cluster-differentiating amino acid substitutions in the influenza virus surface protein hemagglutinin (HA) were in part the same. The antigenic evolution of swine viruses occurred at a rate approximately six times slower than the rate in human viruses, even though the rates of genetic evolution of the HA at the nucleotide and amino acid level were similar for human and swine H3N2 viruses. Continuous monitoring of antigenic changes is recommended to give a first indication as to whether vaccine strains may need updating. Our data suggest that humoral immunity in the population plays a smaller role in the evolutionary selection processes of swine H3N2 viruses than in human H3N2 viruses.     

Dynamics of glycoprotein charge in the evolutionary history of human influenza

Background

Influenza viruses show a significant capacity to evade host immunity; this is manifest both as large occasional jumps in the antigenic phenotype of viral surface molecules and in gradual antigenic changes leading to annual influenza epidemics in humans. Recent mouse studies show that avidity for host cells can play an important role in polyclonal antibody escape, and further that electrostatic charge of the hemagglutinin glycoprotein can contribute to such avidity.

Methodology/Principal Findings

We test the role of glycoprotein charge on sequence data from the three major subtypes of influenza A in humans, using a simple method of calculating net glycoprotein charge. Of all subtypes, H3N2 in humans shows a striking pattern of increasing positive charge since its introduction in 1968. Notably, this trend applies to both hemagglutinin and neuraminidase glycoproteins. In the late 1980s hemagglutinin charge reached a plateau, while neuraminidase charge started to decline. We identify key groups of amino acid sites involved in this charge trend.

Conclusions/Significance

To our knowledge these are the first indications that, for human H3N2, net glycoprotein charge covaries strongly with antigenic drift on a global scale. Further work is needed to elucidate how such charge interacts with other immune escape mechanisms, such as glycosylation, and we discuss important questions arising for future study.     

Long-term excretion of vaccine-derived poliovirus by a healthy child

A child was found to be excreting type 1 vaccine-derived poliovirus (VDPV) with a 1.1% sequence drift from Sabin type 1 vaccine strain in the VP1 coding region 6 months after he was immunized with oral live polio vaccine. Seventeen type 1 poliovirus isolates were recovered from stools taken from this child during the following 4 months. Contrary to expectation, the child was not deficient in humoral immunity and showed high levels of serum neutralization against poliovirus. Selected virus isolates were characterized in terms of their antigenic properties, virulence in transgenic mice, sensitivity for growth at high temperatures, and differences in nucleotide sequence from the Sabin type 1 strain. The VDPV isolates showed mutations at key nucleotide positions that correlated with the observed reversion to biological properties typical of wild polioviruses. A number of capsid mutations mapped at known antigenic sites leading to changes in the viral antigenic structure. Estimates of sequence evolution based on the accumulation of nucleotide changes in the VP1 coding region detected a "defective" molecular clock running at an apparent faster speed of 2.05% nucleotide changes per year versus 1% shown in previous studies. Remarkably, when compared to several type 1 VDPV strains of different origins, isolates from this child showed a much higher proportion of nonsynonymous versus synonymous nucleotide changes in the capsid coding region. This anomaly could explain the high VP1 sequence drift found and the ability of these virus strains to replicate in the gut for a longer period than expected.     

Simulating the Dynamics of Scale-Free Networks via Optimization

We deal here with the issue of complex network evolution. The analysis of topological evolution of complex networks plays a crucial role in predicting their future. While an impressive amount of work has been done on the issue, very little attention has been so far devoted to the investigation of how information theory quantifiers can be applied to characterize networks evolution. With the objective of dynamically capture the topological changes of a network''s evolution, we propose a model able to quantify and reproduce several characteristics of a given network, by using the square root of the Jensen-Shannon divergence in combination with the mean degree and the clustering coefficient. To support our hypothesis, we test the model by copying the evolution of well-known models and real systems. The results show that the methodology was able to mimic the test-networks. By using this copycat model, the user is able to analyze the networks behavior over time, and also to conjecture about the main drivers of its evolution, also providing a framework to predict its evolution.     

Evolution of cooperation through indirect reciprocity

How can cooperation through indirect reciprocity evolve and what would it be like? This problem has previously been studied by simulating evolution in a small group of interacting individuals, assuming no gene flow between groups. In these simulations, certain ''image scoring'' strategies were found to be the most successful. However, analytical arguments show that it would not be in an individual''s interest to use these strategies. Starting with this puzzle, we investigate indirect reciprocity in simulations based on an island model. This has an advantage in that the role of genetic drift can be examined. Our results show that the image scoring strategies depend on very strong drift or a very small cost of giving help. As soon as these factors are absent, selection eliminates image scoring. We also consider other possibilities for the evolution of indirect reciprocity. In particular, we find that the strategy of aiming for ''good standing'' has superior properties. It can be an evolutionarily stable strategy and, even if not, it usually beats image scoring. Furthermore, by introducing quality variation among individuals into the model, we show that the standing strategy can be quality revealing, adding a new dimension to indirect reciprocity. Finally, we discuss general problems with currently popular modelling styles.     

Towards a more biologically realistic use of Droop's equations to model growth under multiple nutrient limitation

Droop's model was originally designed to describe the growth of unicellular phytoplankton species in chemostats but it is now commonly used for a variety of organisms in models of trophic interactions, ecosystem functioning, and evolution. Despite its ubiquitous use, Droop's model is still limited by several simplifying assumptions. For example, the assumption of equal theoretical maximum growth rates for all nutrients is commonly used to describe growth limited by multiple nutrients. This assumption, however, is both biologically unrealistic and potentially misleading. We propose the alternative hypothesis of equal realized maximum growth rates for all nutrients. We support our hypothesis with empirical and theoretical arguments and discuss how it may improve our understanding of the biology of growth, while avoiding some of the pitfalls of the previous assumption.     

High avidity antigen-specific CTL identified by CD8-independent tetramer staining

Tetrameric MHC/peptide complexes are important tools for enumerating, phenotyping, and rapidly cloning Ag-specific T cells. It remains however unclear whether they can reliably distinguish between high and low avidity T cell clones. In this report, tetramers with mutated CD8 binding site selectively stain higher avidity human and murine CTL capable of recognizing physiological levels of Ag. Furthermore, we demonstrate that CD8 binding significantly enhances the avidity as well as the stability of interactions between CTL and cognate tetramers. The use of CD8-null tetramers to identify high avidity CTL provides a tool to compare vaccination strategies for their ability to enhance the frequency of high avidity CTL. Using this technique, we show that DNA priming and vaccinia boosting of HHD A2 transgenic mice fail to selectively expand large numbers of high avidity NY-ESO-1(157-165)-specific CTL, possibly due to the large amounts of antigenic peptide delivered by the vaccinia virus. Furthermore, development of a protocol for rapid identification of high avidity human and murine T cells using tetramers with impaired CD8 binding provides an opportunity not only to monitor expansion of high avidity T cell responses ex vivo, but also to sort high avidity CTL clones for adoptive T cell transfer therapy.     

The roles of wing color pattern and geography in the evolution of Neotropical Preponini butterflies

Diversification rates and evolutionary trajectories are known to be influenced by phenotypic traits and the geographic history of the landscapes that organisms inhabit. One of the most conspicuous traits in butterflies is their wing color pattern, which has been shown to be important in speciation. The evolution of many taxa in the Neotropics has also been influenced by major geological events. Using a dated, species‐level molecular phylogenetic hypothesis for Preponini, a colorful Neotropical butterfly tribe, we evaluated whether diversification rates were constant or varied through time, and how they were influenced by color pattern evolution and biogeographical events. We found that Preponini originated approximately 28 million years ago and that diversification has increased through time consistent with major periods of Andean uplift. Even though some clades show evolutionarily rapid transitions in coloration, contrary to our expectations, these shifts were not correlated with shifts in diversification. Involvement in mimicry with other butterfly groups might explain the rapid changes in dorsal color patterns in this tribe, but such changes have not increased species diversification in this group. However, we found evidence for an influence of major Miocene and Pliocene geological events on the tribe''s evolution. Preponini apparently originated within South America, and range evolution has since been dynamic, congruent with Andean geologic activity, closure of the Panama Isthmus, and Miocene climate variability.     

Rapid Evolution of Pandemic Noroviruses of the GII.4 Lineage

Over the last fifteen years there have been five pandemics of norovirus (NoV) associated gastroenteritis, and the period of stasis between each pandemic has been progressively shortening. NoV is classified into five genogroups, which can be further classified into 25 or more different human NoV genotypes; however, only one, genogroup II genotype 4 (GII.4), is associated with pandemics. Hence, GII.4 viruses have both a higher frequency in the host population and greater epidemiological fitness. The aim of this study was to investigate if the accuracy and rate of replication are contributing to the increased epidemiological fitness of the GII.4 strains. The replication and mutation rates were determined using in vitro RNA dependent RNA polymerase (RdRp) assays, and rates of evolution were determined by bioinformatics. GII.4 strains were compared to the second most reported genotype, recombinant GII.b/GII.3, the rarely detected GII.3 and GII.7 and as a control, hepatitis C virus (HCV). The predominant GII.4 strains had a higher mutation rate and rate of evolution compared to the less frequently detected GII.b, GII.3 and GII.7 strains. Furthermore, the GII.4 lineage had on average a 1.7-fold higher rate of evolution within the capsid sequence and a greater number of non-synonymous changes compared to other NoVs, supporting the theory that it is undergoing antigenic drift at a faster rate. Interestingly, the non-synonymous mutations for all three NoV genotypes were localised to common structural residues in the capsid, indicating that these sites are likely to be under immune selection. This study supports the hypothesis that the ability of the virus to generate genetic diversity is vital for viral fitness.