Sequence amplification via cell passaging creates spurious signals of positive adaptation in influenza virus H3N2 hemagglutinin

Epistasis mediates the evolution of the receptor binding mode in recent human H3N2 hemagglutinin

The receptor-binding site of influenza A virus hemagglutinin partially overlaps with major antigenic sites and constantly evolves. In this study, we observe that mutations G186D and D190N in the hemagglutinin receptor-binding site have coevolved in two recent human H3N2 clades. X-ray crystallography results show that these mutations coordinately drive the evolution of the hemagglutinin receptor binding mode. Epistasis between G186D and D190N is further demonstrated by glycan binding and thermostability analyses. Immunization and neutralization experiments using mouse and human samples indicate that the evolution of receptor binding mode is accompanied by a change in antigenicity. Besides, combinatorial mutagenesis reveals that G186D and D190N, along with other natural mutations in recent H3N2 strains, alter the compatibility with a common egg-adaptive mutation in seasonal influenza vaccines. Overall, our findings elucidate the role of epistasis in shaping the recent evolution of human H3N2 hemagglutinin and substantiate the high evolvability of its receptor-binding mode.

Glyco-engineered MDCK cells display preferred receptors of H3N2 influenza absent in eggs used for vaccines

Evolution of human H3N2 influenza viruses driven by immune selection has narrowed the receptor specificity of the hemagglutinin (HA) to a restricted subset of human-type (Neu5Acα2-6 Gal) glycan receptors that have extended poly-LacNAc (Galβ1-4GlcNAc) repeats. This altered specificity has presented challenges for hemagglutination assays, growth in laboratory hosts, and vaccine production in eggs. To assess the impact of extended glycan receptors on virus binding, infection, and growth, we have engineered N-glycan extended (NExt) cell lines by overexpressing β3-Ν-acetylglucosaminyltransferase 2 in MDCK, SIAT, and hCK cell lines. Of these, SIAT-NExt cells exhibit markedly increased binding of H3 HAs and susceptibility to infection by recent H3N2 virus strains, but without impacting final virus titers. Glycome analysis of these cell lines and allantoic and amniotic egg membranes provide insights into the importance of extended glycan receptors for growth of recent H3N2 viruses and relevance to their production for cell- and egg-based vaccines.

Egg-adaptive mutations of human influenza H3N2 virus are contingent on natural evolution

Egg-adaptive mutations in influenza hemagglutinin (HA) often emerge during the production of egg-based seasonal influenza vaccines, which contribute to the largest share in the global influenza vaccine market. While some egg-adaptive mutations have minimal impact on the HA antigenicity (e.g. G186V), others can alter it (e.g. L194P). Here, we show that the preference of egg-adaptive mutation in human H3N2 HA is strain-dependent. In particular, Thr160 and Asn190, which are found in many recent H3N2 strains, restrict the emergence of L194P but not G186V. Our results further suggest that natural amino acid variants at other HA residues also play a role in determining the preference of egg-adaptive mutation. Consistently, recent human H3N2 strains from different clades acquire different mutations during egg passaging. Overall, these results demonstrate that natural mutations in human H3N2 HA can influence the preference of egg-adaptation mutation, which has important implications in seed strain selection for egg-based influenza vaccine.

MADE: A Computational Tool for Predicting Vaccine Effectiveness for the Influenza A(H3N2) Virus Adapted to Embryonated Eggs

Seasonal Influenza H3N2 virus poses a great threat to public health, but its vaccine efficacy remains suboptimal. One critical step in influenza vaccine production is the viral passage in embryonated eggs. Recently, the strength of egg passage adaptation was found to be rapidly increasing with time driven by convergent evolution at a set of functionally important codons in the hemagglutinin (HA1). In this study, we aim to take advantage of the negative correlation between egg passage adaptation and vaccine effectiveness (VE) and develop a computational tool for selecting the best candidate vaccine virus (CVV) for vaccine production. Using a probabilistic approach known as mutational mapping, we characterized the pattern of sequence evolution driven by egg passage adaptation and developed a new metric known as the adaptive distance (AD) which measures the overall strength of egg passage adaptation. We found that AD is negatively correlated with the influenza H3N2 vaccine effectiveness (VE) and ~75% of the variability in VE can be explained by AD. Based on these findings, we developed a computational package that can Measure the Adaptive Distance and predict vaccine Effectiveness (MADE). MADE provides a powerful tool for the community to calibrate the effect of egg passage adaptation and select more reliable strains with minimum egg-passaged changes as the seasonal A/H3N2 influenza vaccine.

Sequence dynamics of three influenza A virus strains grown in different MDCK cell lines, including those expressing different sialic acid receptors

Viruses are often cultured in cell lines for research and vaccine development, and those often differ from the natural hosts or tissues. Cell lines can also differ in the presence of virus receptors, such as the sialic acid (Sia) receptors used by influenza A viruses (IAV), which can vary in linkage (α2,3- or α2,6-linkage) and form (N-glycolylneuraminic acid [Neu5Gc] or N-acetylneuraminic acid [Neu5Ac]). The selective pressures resulting from passaging viruses in cell types with host-specific variations in viral receptors are still only partially understood. IAV are commonly cultured in MDCK cells which are both derived from canine kidney tubule epithelium and inherently heterogeneous. MDCK cells naturally present Neu5Ac and α2,3-linked Sia forms. Here, we examine natural MDCK variant lineages, as well as engineered variants that synthesize Neu5Gc and/or α2,6-linkages. We determined how viral genetic variation occurred within human H3N2, H1N1 pandemic and canine H3N2 IAV populations when serially passaged in MDCK cell lines that vary in cell type (MDCK-Type I or MDCK-Type II clones) and in Sia display. Deep sequencing of viral genomes showed small numbers of consensus-level mutations, mostly within the hemagglutinin (HA) gene. Both human IAV showed variants in the HA stem and the HA receptor-binding site of populations passaged in cells displaying Neu5Gc. Canine H3N2 showed variants near the receptor-binding site when passaged in cells displaying Neu5Gc or α2,6-linkages. Viruses replicated to low titres in MDCK-Type II cells, suggesting that not all cell types in heterogeneous MDCK cell populations are equally permissive to infection.

Passage Adaptation Correlates With the Reduced Efficacy of the Influenza Vaccine

BACKGROUND

As a dominant seasonal influenza virus, H3N2 virus rapidly evolves in humans and is a constant threat to public health. Despite sustained research efforts, the efficacy of H3N2 vaccine has decreased rapidly. Even though antigenic drift and passage adaptation (substitutions accumulated during vaccine production in embryonated eggs) have been implicated in reduced vaccine efficacy (VE), their respective contributions to the phenomenon remain controversial.

METHODS

We utilized mutational mapping, a powerful probabilistic method for studying sequence evolution, to analyze patterns of substitutions in different passage conditions for an unprecedented amount of H3N2 hemagglutinin sequences (n = 32 278).

RESULTS

We found that passage adaptation in embryonated eggs is driven by repeated convergent evolution over 12 codons. Based on substitution patterns at these sites, we developed a metric, adaptive distance (AD), to quantify the strength of passage adaptation and subsequently identified a strong negative correlation between AD and VE.

CONCLUSIONS

The high correlation between AD and VE implies that passage adaptation in embryonated eggs may be a strong contributor to the recent reduction in H3N2 VE. We developed a computational package called MADE (Measuring Adaptive Distance and vaccine Efficacy based on allelic barcodes) to measure the strength of passage adaptation and predict the efficacy of a candidate vaccine strain. Our findings shed light on strategies for reducing Darwinian evolution within the passaging medium in order to potentially restore an effective vaccine program in the future.

Different genetic barriers for resistance to HA stem antibodies in influenza H3 and H1 viruses

The discovery and characterization of broadly neutralizing human antibodies (bnAbs) to the highly conserved stem region of influenza hemagglutinin (HA) have contributed to considerations of a universal influenza vaccine. However, the potential for resistance to stem bnAbs also needs to be more thoroughly evaluated. Using deep mutational scanning, with a focus on epitope residues, we found that the genetic barrier to resistance to stem bnAbs is low for the H3 subtype but substantially higher for the H1 subtype owing to structural differences in the HA stem. Several strong resistance mutations in H3 can be observed in naturally circulating strains and do not reduce in vitro viral fitness and in vivo pathogenicity. This study highlights a potential challenge for development of a truly universal influenza vaccine.

Within-Host Viral Diversity: A Window into Viral Evolution

The evolutionary dynamics of a virus can differ within hosts and across populations. Studies of within-host evolution provide an important link between experimental studies of virus evolution and large-scale phylodynamic analyses. They can determine the extent to which global processes are recapitulated on local scales and how accurately experimental infections model natural ones. They may also inform epidemiologic models of disease spread and reveal how host-level dynamics contribute to a virus's evolution at a larger scale. Over the last decade, advances in viral sequencing have enabled detailed studies of viral genetic diversity within hosts. I review how within-host diversity is sampled, measured, and expressed, and how comparative studies of viral diversity can be leveraged to elucidate a virus's evolutionary dynamics. These concepts are illustrated with detailed reviews of recent research on the within-host evolution of influenza virus, dengue virus, and cytomegalovirus.

Linking influenza virus evolution within and between human hosts

Influenza viruses rapidly diversify within individual human infections. Several recent studies have deep-sequenced clinical influenza infections to identify viral variation within hosts, but it remains unclear how within-host mutations fare at the between-host scale. Here, we compare the genetic variation of H3N2 influenza within and between hosts to link viral evolutionary dynamics across scales. Synonymous sites evolve at similar rates at both scales, indicating that global evolution at these putatively neutral sites results from the accumulation of within-host variation. However, nonsynonymous mutations are depleted between hosts compared to within hosts, suggesting that selection purges many of the protein-altering changes that arise within hosts. The exception is at antigenic sites, where selection detectably favors nonsynonymous mutations at the global scale, but not within hosts. These results suggest that selection against deleterious mutations and selection for antigenic change are the main forces that act on within-host variants of influenza virus as they transmit and circulate between hosts.

Evolution and rapid spread of a reassortant A(H3N2) virus that predominated the 2017-2018 influenza season

The 2017-2018 North American influenza season caused more hospitalizations and deaths than any year since the 2009 H1N1 pandemic. The majority of recorded influenza infections were caused by A(H3N2) viruses, with most of the virus's North American diversity falling into the A2 clade. Within A2, we observe a subclade which we call A2/re that rose to comprise almost 70 per cent of A(H3N2) viruses circulating in North America by early 2018. Unlike most fast-growing clades, however, A2/re contains no amino acid substitutions in the hemagglutinin (HA) segment. Moreover, hemagglutination inhibition assays did not suggest substantial antigenic differences between A2/re viruses and viruses sampled during the 2016-2017 season. Rather, we observe that the A2/re clade was the result of a reassortment event that occurred in late 2016 or early 2017 and involved the combination of the HA and PB1 segments of an A2 virus with neuraminidase (NA) and other segments a virus from the clade A1b. The success of this clade shows the need for antigenic analysis that targets NA in addition to HA. Our results illustrate the potential for non-HA drivers of viral success and necessitate the need for more thorough tracking of full viral genomes to better understand the dynamics of influenza epidemics.

Cell culture propagation of foot-and-mouth disease virus: adaptive amino acid substitutions in structural proteins and their functional implications

Foot-and-mouth disease is endemic in livestock in large parts of Africa and Asia, where it is an important driver of food insecurity and a major obstacle to agricultural development and the international trade in animal products. Virtually all commercially available vaccines are inactivated whole-virus vaccines produced in cell culture, but the adaptation of a field isolate of the virus to growth in culture is laborious and time-consuming. This is of particular concern for the development of vaccines to newly emerging virus lineages, where long lead times from virus isolate to vaccine can delay the implementation of effective control programs. High antigen yields in production cells are also necessary to make vaccines affordable for less developed countries in endemic areas. Therefore, a rational approach to cell culture adaptation that combines prior knowledge of common adaptive mutations and reverse genetics techniques is urgently required. This review provides an overview of amino acid exchanges in the viral capsid proteins in the context of adaptation to cell culture.

Allele-specific nonstationarity in evolution of influenza A virus surface proteins

Influenza A virus (IAV) is a major public health problem and a pandemic threat. Its evolution is largely driven by diversifying positive selection so that relative fitness of different amino acid variants changes with time due to changes in herd immunity or genomic context, and novel amino acid variants attain fitness advantage. Here, we hypothesize that diversifying selection also has another manifestation: the fitness associated with a particular amino acid variant should decline with time since its origin, as the herd immunity adapts to it. By tracing the evolution of antigenic sites at IAV surface proteins, we show that an amino acid variant becomes progressively more likely to become replaced by another variant with time since its origin-a phenomenon we call "senescence." Senescence is particularly pronounced at experimentally validated antigenic sites, implying that it is largely driven by host immunity. By contrast, at internal sites, existing variants become more favorable with time, probably due to arising contingent mutations at other epistatically interacting sites. Our findings reveal a previously undescribed facet of adaptive evolution and suggest approaches for prediction of evolutionary dynamics of pathogens.

Preventing an Antigenically Disruptive Mutation in Egg-Based H3N2 Seasonal Influenza Vaccines by Mutational Incompatibility

Egg-based seasonal influenza vaccines are the major preventive countermeasure against influenza virus. However, their effectiveness can be compromised when antigenic changes arise from egg-adaptive mutations on influenza hemagglutinin (HA). The L194P mutation is commonly observed in egg-based H3N2 vaccine seed strains and significantly alters HA antigenicity. An approach to prevent L194P would therefore be beneficial. We show that emergence of L194P during egg passaging can be impeded by preexistence of a G186V mutation, revealing strong incompatibility between these mutations. X-ray structures illustrate that individual G186V and L194P mutations have opposing effects on the HA receptor-binding site (RBS), and when both G186V and L194P are present, the RBS is severely disrupted. Importantly, wild-type HA antigenicity is maintained with G186V, but not L194P. Our results demonstrate that these epistatic interactions can be used to prevent the emergence of mutations that adversely alter antigenicity during egg adaptation.

•

Most H3N2 egg isolates carry hemagglutinin mutation G186V or L194P, but not both

•

Hemagglutinin double mutation G186V/L194P is highly deleterious to the virus

•

Hemagglutinin double mutation G186V/L194P disrupts the receptor-binding site

•

Wild-type hemagglutinin antigenicity is maintained in G186V, but not in L194P

Structures and functions linked to genome-wide adaptation of human influenza A viruses

Human influenza A viruses elicit short-term respiratory infections with considerable mortality and morbidity. While H3N2 viruses circulate for more than 50 years, the recent introduction of pH1N1 viruses presents an excellent opportunity for a comparative analysis of the genome-wide evolutionary forces acting on both subtypes. Here, we inferred patches of sites relevant for adaptation, i.e. being under positive selection, on eleven viral protein structures, from all available data since 1968 and correlated these with known functional properties. Overall, pH1N1 have more patches than H3N2 viruses, especially in the viral polymerase complex, while antigenic evolution is more apparent for H3N2 viruses. In both subtypes, NS1 has the highest patch and patch site frequency, indicating that NS1-mediated viral attenuation of host inflammatory responses is a continuously intensifying process, elevated even in the longtime-circulating subtype H3N2. We confirmed the resistance-causing effects of two pH1N1 changes against oseltamivir in NA activity assays, demonstrating the value of the resource for discovering functionally relevant changes. Our results represent an atlas of protein regions and sites with links to host adaptation, antiviral drug resistance and immune evasion for both subtypes for further study.

Insights into the antigenic advancement of influenza A(H3N2) viruses, 2011-2018

Influenza A(H3N2) viruses evade human immunity primarily by acquiring antigenic changes in the haemagglutinin (HA). HA receptor-binding features of contemporary A(H3N2) viruses hinder traditional antigenic characterization using haemagglutination inhibition and promote selection of HA mutants. Thus, alternative approaches are needed to reliably assess antigenic relatedness between circulating viruses and vaccines. We developed a high content imaging-based neutralization test (HINT) to reduce antigenic mischaracterization resulting from virus adaptation to cell culture. Ferret reference antisera were raised using clinical specimens containing viruses representing recent vaccine strains. Analysis of viruses circulating during 2011–2018 showed that gain of an N158-linked glycosylation in HA was a molecular determinant of antigenic distancing between A/Hong Kong/4801/2014-like (clade 3C.2a) and A/Texas/50/2012-like viruses (clade 3C.1), while multiple evolutionary HA F193S substitution were linked to antigenic distancing from A/Switzerland/97152963/2013-like (clade 3C.3a) and further antigenic distancing from A/Texas/50/2012-like viruses. Additionally, a few viruses carrying HA T135K and/or I192T showed reduced neutralization by A/Hong Kong/4801/2014-like antiserum. Notably, this technique elucidated the antigenic characteristics of clinical specimens, enabling direct characterization of viruses produced in vivo, and eliminating in vitro culture, which rapidly alters the genotype/phenotype. HINT is a valuable new antigenic analysis tool for vaccine strain selection.

Influenza passaging annotations: what they tell us and why we should listen

Influenza databases now contain over 100,000 worldwide sequence records for strains influenza A(H3N2) and A(H1N1). Although these data facilitate global research efforts and vaccine development practices, they also represent a stumbling block for researchers because of their confusing and heterogeneous annotation. Unclear passaging annotations are particularly concerning given the recent work highlighting the presence and risk of false adaptation signals introduced by cell passaging of viral isolates. With this in mind, we aim to provide a concise outline of why viruses are passaged, a clear overview of passaging annotation nomenclature currently in use, and suggestions for a standardized nomenclature going forward. Our hope is that this summary will empower researchers and clinicians alike to more easily understand a virus sample's passage history when analyzing influenza sequences.

Deep mutational scanning of hemagglutinin helps predict evolutionary fates of human H3N2 influenza variants

A key goal in the study of influenza virus evolution is to forecast which viral strains will persist and which ones will die out. Here we experimentally measure the effects of all amino acid mutations to the hemagglutinin protein from a human H3N2 influenza strain on viral growth in cell culture. We show that these measurements have utility for distinguishing among viral strains that do and do not succeed in nature. Overall, our work suggests that new high-throughput experimental approaches may be useful for understanding virus evolution in nature.

Human influenza virus rapidly accumulates mutations in its major surface protein hemagglutinin (HA). The evolutionary success of influenza virus lineages depends on how these mutations affect HA’s functionality and antigenicity. Here we experimentally measure the effects on viral growth in cell culture of all single amino acid mutations to the HA from a recent human H3N2 influenza virus strain. We show that mutations that are measured to be more favorable for viral growth are enriched in evolutionarily successful H3N2 viral lineages relative to mutations that are measured to be less favorable for viral growth. Therefore, despite the well-known caveats about cell-culture measurements of viral fitness, such measurements can still be informative for understanding evolution in nature. We also compare our measurements for H3 HA to similar data previously generated for a distantly related H1 HA and find substantial differences in which amino acids are preferred at many sites. For instance, the H3 HA has less disparity in mutational tolerance between the head and stalk domains than the H1 HA. Overall, our work suggests that experimental measurements of mutational effects can be leveraged to help understand the evolutionary fates of viral lineages in nature—but only when the measurements are made on a viral strain similar to the ones being studied in nature.

The influenza virus hemagglutinin head evolves faster than the stalk domain

The limited ability of current influenza virus vaccines to protect from antigenically drifted or shifted viruses creates a public health problem that has led to the need to develop effective, broadly protective vaccines. While current influenza virus vaccines mostly induce an immune response against the immunodominant and variable head domain of the hemagglutinin, the major surface glycoprotein of the virus, the hemagglutinin stalk domain has been identified to harbor neutralizing B-cell epitopes that are conserved among and even between influenza A virus subtypes. A complete understanding of the differences in evolution between the main target of current vaccines and this more conserved stalk region are missing. Here, we performed an evolutionary analysis of the stalk domains of the hemagglutinin of pre-pandemic seasonal H1N1, pandemic H1N1, seasonal H3N2, and influenza B viruses and show quantitatively for the first time that the stalk domain is evolving at a rate that is significantly slower than that of the head domain. Additionally, we found that the cross-reactive epitopes in the stalk domain targeted by broadly neutralizing monoclonal antibodies are evolving at an even slower rate compared to the full head and stalk regions of the protein. Finally, a fixed-effects likelihood selection analysis was performed for these virus groups in both the head and stalk domains. While several positive selection sites were found in the head domain, only a single site in the stalk domain of pre-pandemic seasonal H1 hemagglutinin was identified at amino acid position 468 (H1 numbering from methionine). This site is not located in or close to the epitopes of cross-reactive anti-stalk monoclonal antibodies. Furthermore, we found that changes in this site do not significantly impact virus binding or neutralization by human anti-stalk antibodies, suggesting that some positive selection in the stalk domain is independent of immune pressures. We conclude that, while the stalk domain does evolve over time, this evolution is slow and, historically, is not directed to aid in evading neutralizing antibody responses.

Hsp90 shapes protein and RNA evolution to balance trade-offs between protein stability and aggregation

Acquisition of mutations is central to evolution; however, the detrimental effects of most mutations on protein folding and stability limit protein evolvability. Molecular chaperones, which suppress aggregation and facilitate polypeptide folding, may alleviate the effects of destabilizing mutations thus promoting sequence diversification. To illuminate how chaperones can influence protein evolution, we examined the effect of reduced activity of the chaperone Hsp90 on poliovirus evolution. We find that Hsp90 offsets evolutionary trade-offs between protein stability and aggregation. Lower chaperone levels favor variants of reduced hydrophobicity and protein aggregation propensity but at a cost to protein stability. Notably, reducing Hsp90 activity also promotes clusters of codon-deoptimized synonymous mutations at inter-domain boundaries, likely to facilitate cotranslational domain folding. Our results reveal how a chaperone can shape the sequence landscape at both the protein and RNA levels to harmonize competing constraints posed by protein stability, aggregation propensity, and translation rate on successful protein biogenesis.

Predictive Modeling of Influenza Shows the Promise of Applied Evolutionary Biology

Seasonal influenza is controlled through vaccination campaigns. Evolution of influenza virus antigens means that vaccines must be updated to match novel strains, and vaccine effectiveness depends on the ability of scientists to predict nearly a year in advance which influenza variants will dominate in upcoming seasons. In this review, we highlight a promising new surveillance tool: predictive models. Based on data-sharing and close collaboration between the World Health Organization and academic scientists, these models use surveillance data to make quantitative predictions regarding influenza evolution. Predictive models demonstrate the potential of applied evolutionary biology to improve public health and disease control. We review the state of influenza predictive modeling and discuss next steps and recommendations to ensure that these models deliver upon their considerable biomedical promise.

Cooperating H3N2 Influenza Virus Variants Are Not Detectable in Primary Clinical Samples

Viruses mutate rapidly, and recent studies of RNA viruses have shown that related viral variants can sometimes cooperate to improve each other’s growth. We previously described two variants of H3N2 influenza virus that cooperate in cell culture. The mutation responsible for cooperation is often observed when human samples of influenza virus are grown in the lab before sequencing, but it is unclear whether the mutation also exists in human infections or is exclusively the result of lab passage. We identified nine human isolates of influenza virus that had developed the cooperating mutation after being grown in the lab and performed highly sensitive deep sequencing of the unpassaged clinical samples to determine whether the mutation existed in the original human infections. We found no evidence of the cooperating mutation in the unpassaged samples, suggesting that the cooperation arises primarily under laboratory conditions.

The high mutation rates of RNA viruses lead to rapid genetic diversification, which can enable cooperative interactions between variants in a viral population. We previously described two distinct variants of H3N2 influenza virus that cooperate in cell culture. These variants differ by a single mutation, D151G, in the neuraminidase protein. The D151G mutation reaches a stable frequency of about 50% when virus is passaged in cell culture. However, it is unclear whether selection for the cooperative benefits of D151G is a cell culture phenomenon or whether the mutation is also sometimes present at appreciable frequency in virus populations sampled directly from infected humans. Prior work has not detected D151G in unpassaged clinical samples, but those studies have used methods like Sanger sequencing and pyrosequencing, which are relatively insensitive to low-frequency variation. We identified nine samples of human H3N2 influenza virus collected between 2013 and 2015 in which Sanger sequencing had detected a high frequency of the D151G mutation following one to three passages in cell culture. We deep sequenced the unpassaged clinical samples to identify low-frequency viral variants. The frequency of D151G did not exceed the frequency of library preparation and sequencing errors in any of the sequenced samples. We conclude that passage in cell culture is primarily responsible for the frequent observations of D151G in recent H3N2 influenza virus strains.

IMPORTANCE Viruses mutate rapidly, and recent studies of RNA viruses have shown that related viral variants can sometimes cooperate to improve each other’s growth. We previously described two variants of H3N2 influenza virus that cooperate in cell culture. The mutation responsible for cooperation is often observed when human samples of influenza virus are grown in the lab before sequencing, but it is unclear whether the mutation also exists in human infections or is exclusively the result of lab passage. We identified nine human isolates of influenza virus that had developed the cooperating mutation after being grown in the lab and performed highly sensitive deep sequencing of the unpassaged clinical samples to determine whether the mutation existed in the original human infections. We found no evidence of the cooperating mutation in the unpassaged samples, suggesting that the cooperation arises primarily under laboratory conditions.

Integrating patient and whole-genome sequencing data to provide insights into the epidemiology of seasonal influenza A(H3N2) viruses

Genetic surveillance of seasonal influenza is largely focused on sequencing of the haemagglutinin gene. Consequently, our understanding of the contribution of the remaining seven gene segments to the evolution and epidemiological dynamics of seasonal influenza is relatively limited. The increased availability of next-generation sequencing technologies allows rapid and economic whole-genome sequencing (WGS) of influenza virus. Here, 150 influenza A(H3N2) positive clinical specimens with linked epidemiological data, from the 2014/15 season in Scotland, were sequenced directly using both Sanger sequencing of the HA1 region and WGS using the Illumina MiSeq platform. Sequences generated by the two methods were highly correlated, and WGS provided on average >90 % whole genome coverage. As reported in other European countries during 2014/15, all strains belonged to genetic group 3C, with subgroup 3C.2a predominating. Multiple inter-subgroup reassortants were identified, including three 3C.3 viruses descended from a single reassortment event, which had persisted in the population. Cases of severe acute respiratory illness were significantly clustered on phylogenies of multiple gene segments indicating potential genetic factors warranting further investigation. Severe cases were also more likely to be associated with reassortant viruses and to occur later in the season. These results suggest that WGS provides an opportunity to develop our understanding of the relationship between the influenza genome and disease severity and the epidemiological consequences of within-subtype reassortment. Therefore, increased levels of WGS, linked to clinical and epidemiological data, could improve influenza surveillance.

Fitness cost of reassortment in human influenza

Reassortment, which is the exchange of genome sequence between viruses co-infecting a host cell, plays an important role in the evolution of segmented viruses. In the human influenza virus, reassortment happens most frequently between co-existing variants within the same lineage. This process breaks genetic linkage and fitness correlations between viral genome segments, but the resulting net effect on viral fitness has remained unclear. In this paper, we determine rate and average selective effect of reassortment processes in the human influenza lineage A/H3N2. For the surface proteins hemagglutinin and neuraminidase, reassortant variants with a mean distance of at least 3 nucleotides to their parent strains get established at a rate of about 10⁻² in units of the neutral point mutation rate. Our inference is based on a new method to map reassortment events from joint genealogies of multiple genome segments, which is tested by extensive simulations. We show that intra-lineage reassortment processes are, on average, under substantial negative selection that increases in strength with increasing sequence distance between the parent strains. The deleterious effects of reassortment manifest themselves in two ways: there are fewer reassortment events than expected from a null model of neutral reassortment, and reassortant strains have fewer descendants than their non-reassortant counterparts. Our results suggest that influenza evolves under ubiquitous epistasis across proteins, which produces fitness barriers against reassortment even between co-circulating strains within one lineage.

The genome of the human influenza virus consists of 8 disjoint RNA polymer segments. These segments can undergo reassortment: when two viruses co-infect a host cell, they can produce viral offspring with a new combination of segments. In this paper, we show that reassortment within a given influenza lineage induces a fitness cost that increases in strength with increasing genetic distance of the parent viruses. Our finding suggests that evolution continuously produces viral proteins whose fitness depends on each other; reassortment reduces fitness by breaking up successful combinations of proteins. Thus, selection across proteins constrains viral evolution within a given lineage, and it may be an important factor in defining a viral species.

Full-genome dengue virus sequencing in mosquito saliva shows lack of convergent positive selection during transmission by Aedes aegypti

Like other pathogens with high mutation and replication rates, within-host dengue virus (DENV) populations evolve during infection of their main mosquito vector, Aedes aegypti. Within-host DENV evolution during transmission provides opportunities for adaptation and emergence of novel virus variants. Recent studies of DENV genetic diversity failed to detect convergent evolution of adaptive mutations in mosquito tissues such as midgut and salivary glands, suggesting that convergent positive selection is not a major driver of within-host DENV evolution in the vector. However, it is unknown whether this conclusion extends to the transmitted viral subpopulation because it is technically difficult to sequence DENV genomes in mosquito saliva. Here, we achieved DENV full-genome sequencing by pooling saliva samples collected non-sacrificially from 49 to 163 individual Ae. aegypti mosquitoes previously infected with one of two DENV-1 genotypes. We compared the transmitted viral subpopulations found in the pooled saliva samples collected in time series with the input viral population present in the infectious blood meal. In all pooled saliva samples examined, the full-genome consensus sequence of the input viral population was unchanged. Although the pooling strategy prevents analysis of individual saliva samples, our results demonstrate the lack of strong convergent positive selection during a single round of DENV transmission by Ae. aegypti. This finding reinforces the idea that genetic drift and purifying selection are the dominant evolutionary forces shaping within-host DENV genetic diversity during transmission by mosquitoes.

A structural explanation for the low effectiveness of the seasonal influenza H3N2 vaccine

The effectiveness of the annual influenza vaccine has declined in recent years, especially for the H3N2 component, and is a concern for global public health. A major cause for this lack in effectiveness has been attributed to the egg-based vaccine production process. Substitutions on the hemagglutinin glycoprotein (HA) often arise during virus passaging that change its antigenicity and hence vaccine effectiveness. Here, we characterize the effect of a prevalent substitution, L194P, in egg-passaged H3N2 viruses. X-ray structural analysis reveals that this substitution surprisingly increases the mobility of the 190-helix and neighboring regions in antigenic site B, which forms one side of the receptor binding site (RBS) and is immunodominant in recent human H3N2 viruses. Importantly, the L194P substitution decreases binding and neutralization by an RBS-targeted broadly neutralizing antibody by three orders of magnitude and significantly changes the HA antigenicity as measured by binding of human serum antibodies. The receptor binding mode and specificity are also altered to adapt to avian receptors during egg passaging. Overall, these findings help explain the low effectiveness of the seasonal vaccine against H3N2 viruses, and suggest that alternative approaches should be accelerated for producing influenza vaccines as well as isolating clinical isolates.

Seasonal influenza vaccine does not always confer protection in vaccinated individuals. Vaccine candidates are selected from clinical isolates based on their antigenic properties. It is common to use chicken eggs for culturing clinical isolates and for large-scale production of vaccines. However, influenza virus often mutates to adapt to being grown in chicken eggs, which can influence antigenicity and hence vaccine effectiveness. Here, we structurally characterize an egg-adaptive substitution, namely L194P, in H3N2 virus hemagglutinin. Our results reveal that the L194P substitution substantially increases the flexibility of an epitope region that is commonly targeted by antibodies. Based on the binding affinity of a broadly neutralizing antibody and a panel of human serum antibodies, we further show that the L194P substitution dramatically changes the HA antigenicity. The change of the receptor-binding mode associated with the L194P substitution provides an explanation for its ability to successfully grow in eggs. Our study describes a mechanism for the low influenza vaccine effectiveness and reaffirms the urgency for replacing the egg-based production of influenza vaccines.

In Silico Vaccine Strain Prediction for Human Influenza Viruses

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

Identification of positive selection in genes is greatly improved by using experimentally informed site-specific models

Background

Sites of positive selection are identified by comparing observed evolutionary patterns to those expected under a null model for evolution in the absence of such selection. For protein-coding genes, the most common null model is that nonsynonymous and synonymous mutations fix at equal rates; this unrealistic model has limited power to detect many interesting forms of selection.

Results

I describe a new approach that uses a null model based on experimental measurements of a gene’s site-specific amino-acid preferences generated by deep mutational scanning in the lab. This null model makes it possible to identify both diversifying selection for repeated amino-acid change and differential selection for mutations to amino acids that are unexpected given the measurements made in the lab. I show that this approach identifies sites of adaptive substitutions in four genes (lactamase, Gal4, influenza nucleoprotein, and influenza hemagglutinin) far better than a comparable method that simply compares the rates of nonsynonymous and synonymous substitutions.

Conclusions

As rapid increases in biological data enable increasingly nuanced descriptions of the constraints on individual protein sites, approaches like the one here can improve our ability to identify many interesting forms of selection in natural sequences.

Reviewers

This article was reviewed by Sebastian Maurer-Stroh, Olivier Tenaillon, and Tal Pupko. All three reviewers are members of the Biology Direct editorial board.

Electronic supplementary material

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

Viral factors in influenza pandemic risk assessment

The threat of an influenza A virus pandemic stems from continual virus spillovers from reservoir species, a tiny fraction of which spark sustained transmission in humans. To date, no pandemic emergence of a new influenza strain has been preceded by detection of a closely related precursor in an animal or human. Nonetheless, influenza surveillance efforts are expanding, prompting a need for tools to assess the pandemic risk posed by a detected virus. The goal would be to use genetic sequence and/or biological assays of viral traits to identify those non-human influenza viruses with the greatest risk of evolving into pandemic threats, and/or to understand drivers of such evolution, to prioritize pandemic prevention or response measures. We describe such efforts, identify progress and ongoing challenges, and discuss three specific traits of influenza viruses (hemagglutinin receptor binding specificity, hemagglutinin pH of activation, and polymerase complex efficiency) that contribute to pandemic risk.