Heterogeneity of the mutation rates of influenza A viruses: isolation of mutator mutants

Investigating the Genetic Diversity of H5 Avian Influenza Viruses in the United Kingdom from 2020-2022

Since 2020, the United Kingdom and Europe have experienced annual epizootics of high-pathogenicity avian influenza virus (HPAIV). The first epizootic, during the autumn/winter of 2020-2021, involved six H5Nx subtypes, although H5N8 HPAIV dominated in the United Kingdom. While genetic assessments of the H5N8 HPAIVs within the United Kingdom demonstrated relative homogeneity, there was a background of other genotypes circulating at a lower degree with different neuraminidase and internal genes.  Following a small number of detections of H5N1 in wild birds over the summer of 2021, the autumn/winter of 2021-2022 saw another European H5 HPAIV epizootic that dwarfed the prior epizootic. This second epizootic was dominated almost exclusively by H5N1 HPAIV, although six distinct genotypes were defined. We have used genetic analysis to evaluate the emergence of different genotypes and proposed reassortment events that have been observed. The existing data suggest that the H5N1 viruses circulating in Europe during late 2020 continued to circulate in wild birds throughout 2021, with minimal adaptation, but then went on to reassort with AIVs in the wild bird population. We have undertaken an in-depth genetic assessment of H5 HPAIVs detected in the United Kingdom over two winter seasons and demonstrate the utility of in-depth genetic analyses in defining the diversity of H5 HPAIVs circulating in avian species, the potential for zoonotic risk, and whether incidents of lateral spread can be defined over independent incursions of infections from wild birds. This provides key supporting data for mitigation activities. IMPORTANCE High-pathogenicity avian influenza virus (HPAIV) outbreaks devastate avian species across all sectors, having both economic and ecological impacts through mortalities in poultry and wild birds, respectively. These viruses can also represent a significant zoonotic risk. Since 2020, the United Kingdom has experienced two successive outbreaks of H5 HPAIV. While H5N8 HPAIV was predominant during the 2020-2021 outbreak, other H5 subtypes were also detected. The following year, there was a shift in the subtype dominance to H5N1 HPAIV, but multiple H5N1 genotypes were detected. Through the thorough utilization of whole-genome sequencing, it was possible to track and characterize the genetic evolution of these H5 HPAIVs in United Kingdom poultry and wild birds. This enabled us to assess the risk posed by these viruses at the poultry-wild bird and the avian-human interfaces and to investigate the potential lateral spread between infected premises, a key factor in understanding the threat to the commercial sector.

Predicting influenza with pandemic-awareness via Dynamic Virtual Graph Significance Networks

Every year, influenza spreads worldwide and burdens people's health substantially. We need a reliable model to help hospitals, pharmaceutical companies, and governments better prepare for influenza outbreaks in a timely manner. However, the domain knowledge for such public health events, such as the variable influenza seasonality and occasional pandemics, poses significant challenges in predicting influenza outbreaks. The existing methods use current and historical values in a user-defined time window as input to predict future values but lack considering the situations outside the window. To address these limitations, we proposed Dynamic Virtual Graph Significance Networks (DVGSN). The graph-based algorithm can supervisedly and dynamically learn the implied knowledge from similar "infection situations" in all the historical timepoints without the limitation of time window. Furthermore, representation learning on the dynamic virtual graph can tackle the varied seasonality with pandemic-awareness without requiring domain knowledge input. The extensive experiments on real-world influenza data demonstrate that DVGSN significantly outperforms the state-of-the-art methods. To the best of our knowledge, this is the first attempt to supervisedly learn a dynamic virtual graph for time-series prediction tasks. Moreover, the proposed method has rich interpretabilities, which makes the method more acceptable in the fields of public health, life sciences, and so on. Our source code and dataset are available at https://github.com/aI-area/DVGSN.

Generation of Human Lung Organoid Cultures from Healthy and Tumor Tissue to Study Infectious Diseases

Three-dimensional (3D) human lung organoids reflect the native cell composition of the lung as well as its physiological properties. Human 3D lung organoids offer ideal conditions, such as timely availability in large quantities and high physiological relevance for reassessment and prediction of disease outbreaks of respiratory pathogens and pathogens that use the lung as a primary entry portal.

Mutations in PB2 and HA are crucial for the increased virulence and transmissibility of H1N1 swine influenza virus in mammalian models

Genetic analyses indicated that the pandemic H1N1/2009 influenza virus originated from a swine influenza virus (SIV). However, SIVs bearing the same constellation of genetic features as H1N1/2009 have not been isolated. Understanding the adaptation of SIVs with such genotypes in a new host may provide clues regarding the emergence of pandemic strains such as H1N1/2009. In this study, an artificial SIV with the H1N1/2009 genotype (rH1N1) was sequentially passaged in mice through two independent series, yielding multiple mouse-adapted mutants with high genetic diversity and increased virulence. These experiments were meant to mimic genetic bottlenecks during adaptation of wild viruses with rH1N1 genotypes in a new host. Molecular substitutions in the mouse-adapted variants mainly occurred in genes encoding surface proteins (hemagglutinin [HA] and neuraminidase [NA]) and polymerase proteins (polymerase basic 2 [PB2], polymerase basic 1 [PB1], polymerase acid [PA] proteins and nucleoprotein [NP]). The PB2D309N and HAL425M substitutions were detected at high frequencies in both passage lines and enhanced the replication and pathogenicity of rH1N1 in mice. Moreover, these substitutions also enabled direct transmission of rH1N1 in other mammals such as guinea pigs. PB2D309N showed enhanced polymerase activity and HAL425M showed increased stability compared with the wild-type proteins. Our findings indicate that if SIVs with H1N1/2009 genotypes emerge in pigs, they could undergo rapid adaptive changes during infection of a new host, especially in the PB2 and HA genes. These changes may facilitate the emergence of pandemic strains such as H1N1/2009.

Mutation Rates, Mutation Frequencies, and Proofreading-Repair Activities in RNA Virus Genetics

The error rate displayed during template copying to produce viral RNA progeny is a biologically relevant parameter of the replication complexes of viruses. It has consequences for virus-host interactions, and it represents the first step in the diversification of viruses in nature. Measurements during infections and with purified viral polymerases indicate that mutation rates for RNA viruses are in the range of 10-3 to 10-6 copying errors per nucleotide incorporated into the nascent RNA product. Although viruses are thought to exploit high error rates for adaptation to changing environments, some of them possess misincorporation correcting activities. One of them is a proofreading-repair 3' to 5' exonuclease present in coronaviruses that may decrease the error rate during replication. Here we review experimental evidence and models of information maintenance that explain why elevated mutation rates have been preserved during the evolution of RNA (and some DNA) viruses. The models also offer an interpretation of why error correction mechanisms have evolved to maintain the stability of genetic information carried out by large viral RNA genomes such as the coronaviruses.

Quantifying within-host diversity of H5N1 influenza viruses in humans and poultry in Cambodia

Avian influenza viruses (AIVs) periodically cross species barriers and infect humans. The likelihood that an AIV will evolve mammalian transmissibility depends on acquiring and selecting mutations during spillover, but data from natural infection is limited. We analyze deep sequencing data from infected humans and domestic ducks in Cambodia to examine how H5N1 viruses evolve during spillover. Overall, viral populations in both species are predominated by low-frequency (<10%) variation shaped by purifying selection and genetic drift, and half of the variants detected within-host are never detected on the H5N1 virus phylogeny. However, we do detect a subset of mutations linked to human receptor binding and replication (PB2 E627K, HA A150V, and HA Q238L) that arose in multiple, independent humans. PB2 E627K and HA A150V were also enriched along phylogenetic branches leading to human infections, suggesting that they are likely human-adaptive. Our data show that H5N1 viruses generate putative human-adapting mutations during natural spillover infection, many of which are detected at >5% frequency within-host. However, short infection times, genetic drift, and purifying selection likely restrict their ability to evolve extensively during a single infection. Applying evolutionary methods to sequence data, we reveal a detailed view of H5N1 virus adaptive potential, and develop a foundation for studying host-adaptation in other zoonotic viruses.

H5N1 avian influenza viruses can cross species barriers and cause severe disease in humans. H5N1 viruses currently cannot replicate and transmit efficiently among humans, but animal infection studies and modeling experiments have suggested that human adaptation may require only a few mutations. However, data from natural spillover infection has been limited, posing a challenge for risk assessment. Here, we analyze a unique dataset of deep sequence data from H5N1 virus-infected humans and domestic ducks in Cambodia. We find that well-known markers of human receptor binding and replication arise in multiple, independent humans. We also find that 3 mutations detected within-host are enriched along phylogenetic branches leading to human infections, suggesting that they are likely human-adapting. However, we also show that within-host evolution in both humans and ducks are shaped heavily by purifying selection and genetic drift, and that a large fraction of within-host variation is never detected on the H5N1 phylogeny. Taken together, our data show that H5N1 viruses do generate human-adapting mutations during natural infection. However, short infection times, purifying selection, and genetic drift may severely limit how much H5N1 viruses can evolve during the course of a single infection.

Molecular basis of genetic variation of viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation; mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physical-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate. A distinction is made between mechanistically unavoidable and evolutionarily relevant mutation and recombination.

Selection on mutators is not frequency-dependent

The evolutionary fate of mutator mutations – genetic variants that raise the genome-wide mutation rate – in asexual populations is often described as being frequency (or number) dependent. Mutators can invade a population by hitchhiking with a sweeping beneficial mutation, but motivated by earlier experiments results, it has been repeatedly suggested that mutators must be sufficiently frequent to produce such a driver mutation before non-mutators do. Here, we use stochastic, agent-based simulations to show that neither the strength nor the sign of selection on mutators depend on their initial frequency, and while the overall probability of hitchhiking increases predictably with frequency, the per-capita probability of fixation remains unchanged.

Tyr82 Amino Acid Mutation in PB1 Polymerase Induces an Influenza Virus Mutator Phenotype

In various positive-sense single-stranded RNA viruses, a low-fidelity viral RNA-dependent RNA polymerase (RdRp) confers attenuated phenotypes by increasing the mutation frequency. We report a negative-sense single-stranded RNA virus RdRp mutant strain with a mutator phenotype. Based on structural data of RdRp, rational targeting of key residues, and screening of fidelity variants, we isolated a novel low-fidelity mutator strain of influenza virus that harbors a Tyr82-to-Cys (Y82C) single-amino-acid substitution in the PB1 polymerase subunit. The purified PB1-Y82C polymerase indeed showed an increased frequency of misincorporation compared with the wild-type PB1 in an in vitro biochemical assay. To further investigate the effects of position 82 on PB1 polymerase fidelity, we substituted various amino acids at this position. As a result, we isolated various novel mutators other than PB1-Y82C with higher mutation frequencies. The structural model of influenza virus polymerase complex suggested that the Tyr82 residue, which is located at the nucleoside triphosphate entrance tunnel, may influence a fidelity checkpoint. Interestingly, although the PB1-Y82C variant replicated with wild-type PB1-like kinetics in tissue culture, the 50% lethal dose of the PB1-Y82C mutant was 10 times lower than that of wild-type PB1 in embryonated chicken eggs. In conclusion, our data indicate that the Tyr82 residue of PB1 has a crucial role in regulating polymerase fidelity of influenza virus and is closely related to attenuated pathogenic phenotypes in vivo IMPORTANCE Influenza A virus rapidly acquires antigenic changes and antiviral drug resistance, which limit the effectiveness of vaccines and drug treatments, primarily owing to its high rate of evolution. Virus populations formed by quasispecies can contain resistance mutations even before a selective pressure is applied. To study the effects of the viral mutation spectrum and quasispecies, high- and low-fidelity variants have been isolated for several RNA viruses. Here, we report the discovery of a low-fidelity RdRp variant of influenza A virus that contains a substitution at Tyr82 in PB1. Viruses containing the PB1-Y82C substitution showed growth kinetics and viral RNA synthesis levels similar to those of the wild-type virus in cell culture; however, they had significantly attenuated phenotypes in a chicken egg infection experiment. These data demonstrated that decreased RdRp fidelity attenuates influenza A virus in vivo, which is a desirable feature for the development of safer live attenuated vaccine candidates.

Host susceptibility to severe influenza A virus infection

Most people exposed to a new flu virus do not notice any symptoms. A small minority develops critical illness. Some of this extremely broad variation in susceptibility is explained by the size of the initial inoculum or the influenza exposure history of the individual; some is explained by generic host factors, such as frailty, that decrease resilience following any systemic insult. Some demographic factors (pregnancy, obesity, and advanced age) appear to confer a more specific susceptibility to severe illness following infection with influenza viruses. As with other infectious diseases, a substantial component of susceptibility is determined by host genetics. Several genetic susceptibility variants have now been reported with varying levels of evidence. Susceptible hosts may have impaired intracellular controls of viral replication (e.g. IFITM3, TMPRS22 variants), defective interferon responses (e.g. GLDC, IRF7/9 variants), or defects in cell-mediated immunity with increased baseline levels of systemic inflammation (obesity, pregnancy, advanced age). These mechanisms may explain the prolonged viral replication reported in critically ill patients with influenza: patients with life-threatening disease are, by definition, abnormal hosts. Understanding these molecular mechanisms of susceptibility may in the future enable the design of host-directed therapies to promote resilience.

Influenza A Virus Subpopulations and Their Implication in Pathogenesis and Vaccine Development

The concept of influenza A virus (IAV) subpopulations emerged approximately 75 years ago, when Preben von Magnus described "incomplete" virus particles that interfere with the replication of infectious virus. It is now widely accepted that infectious particles constitute only a minor portion of biologically active IAV subpopulations. The IAV quasispecies is an extremely diverse swarm of biologically and genetically heterogeneous particle subpopulations that collectively influence the evolutionary fitness of the virus. This review summarizes the current knowledge of IAV subpopulations, focusing on their biologic and genomic diversity. It also discusses the potential roles IAV subpopulations play in virus pathogenesis and live attenuated influenza vaccine development.

Single-Cell Virus Sequencing of Influenza Infections That Trigger Innate Immunity

Influenza virus-infected cells vary widely in their expression of viral genes and only occasionally activate innate immunity. Here, we develop a new method to assess how the genetic variation in viral populations contributes to this heterogeneity. We do this by determining the transcriptome and full-length sequences of all viral genes in single cells infected with a nominally "pure" stock of influenza virus. Most cells are infected by virions with defects, some of which increase the frequency of innate-immune activation. These immunostimulatory defects are diverse and include mutations that perturb the function of the viral polymerase protein PB1, large internal deletions in viral genes, and failure to express the virus's interferon antagonist NS1. However, immune activation remains stochastic in cells infected by virions with these defects and occasionally is triggered even by virions that express unmutated copies of all genes. Our work shows that the diverse spectrum of defects in influenza virus populations contributes to-but does not completely explain-the heterogeneity in viral gene expression and immune activation in single infected cells.IMPORTANCE Because influenza virus has a high mutation rate, many cells are infected by mutated virions. But so far, it has been impossible to fully characterize the sequence of the virion infecting any given cell, since conventional techniques such as flow cytometry and single-cell transcriptome sequencing (scRNA-seq) only detect if a protein or transcript is present, not its sequence. Here we develop a new approach that uses long-read PacBio sequencing to determine the sequences of virions infecting single cells. We show that viral genetic variation explains some but not all of the cell-to-cell variability in viral gene expression and innate immune induction. Overall, our study provides the first complete picture of how viral mutations affect the course of infection in single cells.

Migration promotes mutator alleles in subdivided populations

Mutator alleles that elevate the genomic mutation rate may invade nonrecombining populations by hitchhiking with beneficial mutations. Mutators have been repeatedly observed to take over adapting laboratory populations and have been found at high frequencies in both microbial pathogen and cancer populations in nature. Recently, we have shown that mutators are only favored by selection in sufficiently large populations and transition to being disfavored as population size decreases. This population size-dependent sign inversion in selective effect suggests that population structure may also be an important determinant of mutation rate evolution. Although large populations may favor mutators, subdividing such populations into sufficiently small subpopulations (demes) might effectively inhibit them. On the other hand, migration between small demes that otherwise inhibit hitchhiking may promote mutator fixation in the whole metapopulation. Here, we use stochastic, agent-based simulations and evolution experiments with the yeast Saccharomyces cerevisiae to show that mutators can, indeed, be favored by selection in subdivided metapopulations composed of small demes connected by sufficient migration. In fact, we show that population structure plays a previously unsuspected role in promoting mutator success in subdivided metapopulations when migration is rare.

The parts are greater than the whole: the role of semi-infectious particles in influenza A virus biology

The influenza A virus (IAV) genome is incorporated into newly produced virions through a tightly orchestrated process that is one of the best studied examples of genome packaging by a segmented virus. Despite the remarkable selectivity and efficiency of this process, it is clear that the vast majority of IAV virions are unable to express the full set of essential viral gene products and are thus incapable of productive replication in the absence of complementation. Here, we attempt to reconcile the widespread production of these semi-infectious particles (SIPs) with the high efficiency and selectivity of IAV genome packaging. We also cover what is known and what remains unknown about the consequences of SIP production for the replication and evolution of viral populations.

Multi-step prediction for influenza outbreak by an adjusted long short-term memory

Influenza results in approximately 3–5 million annual cases of severe illness and 250 000–500 000 deaths. We urgently need an accurate multi-step-ahead time-series forecasting model to help hospitals to perform dynamical assignments of beds to influenza patients for the annually varied influenza season, and aid pharmaceutical companies to formulate a flexible plan of manufacturing vaccine for the yearly different influenza vaccine. In this study, we utilised four different multi-step prediction algorithms in the long short-term memory (LSTM). The result showed that implementing multiple single-output prediction in a six-layer LSTM structure achieved the best accuracy. The mean absolute percentage errors from two- to 13-step-ahead prediction for the US influenza-like illness rates were all <15%, averagely 12.930%. To the best of our knowledge, it is the first time that LSTM has been applied and refined to perform multi-step-ahead prediction for influenza outbreaks. Hopefully, this modelling methodology can be applied in other countries and therefore help prevent and control influenza worldwide.

Sign of selection on mutation rate modifiers depends on population size

The influence of population size (N) on natural selection acting on alleles that affect fitness has been understood for almost a century. As N declines, genetic drift overwhelms selection and alleles with direct fitness effects are rendered neutral. Often, however, alleles experience so-called indirect selection, meaning they affect not the fitness of an individual but the fitness distribution of its offspring. Some of the best-studied examples of indirect selection include alleles that modify aspects of the genetic system such as recombination and mutation rates. Here, we use analytics, simulations, and experimental populations of Saccharomyces cerevisiae to examine the influence of N on indirect selection acting on alleles that increase the genomic mutation rate (mutators). Mutators experience indirect selection via genomic associations with beneficial and deleterious mutations they generate. We show that, as N declines, indirect selection driven by linked beneficial mutations is overpowered by drift before drift can neutralize the cost of the deleterious load. As a result, mutators transition from being favored by indirect selection in large populations to being disfavored as N declines. This surprising phenomenon of sign inversion in selective effect demonstrates that indirect selection on mutators exhibits a profound and qualitatively distinct dependence on N.

Extreme heterogeneity of influenza virus infection in single cells

Viral infection can dramatically alter a cell's transcriptome. However, these changes have mostly been studied by bulk measurements on many cells. Here we use single-cell mRNA sequencing to examine the transcriptional consequences of influenza virus infection. We find extremely wide cell-to-cell variation in the productivity of viral transcription - viral transcripts comprise less than a percent of total mRNA in many infected cells, but a few cells derive over half their mRNA from virus. Some infected cells fail to express at least one viral gene, but this gene absence only partially explains variation in viral transcriptional load. Despite variation in viral load, the relative abundances of viral mRNAs are fairly consistent across infected cells. Activation of innate immune pathways is rare, but some cellular genes co-vary in abundance with the amount of viral mRNA. Overall, our results highlight the complexity of viral infection at the level of single cells.

Population Diversity and Collective Interactions during Influenza Virus Infection

Influenza A virus (IAV) continues to pose an enormous and unpredictable global public health threat, largely due to the continual evolution of escape from preexisting immunity and the potential for zoonotic emergence. Understanding how the unique genetic makeup and structure of IAV populations influences their transmission and evolution is essential for developing more-effective vaccines, therapeutics, and surveillance capabilities. Owing to their mutation-prone replicase and unique genome organization, IAV populations exhibit enormous amounts of diversity both in terms of sequence and functional gene content. Here, I review what is currently known about the genetic and genomic diversity present within IAV populations and how this diversity may shape the replicative and evolutionary dynamics of these viruses.

Epistatic Interactions within the Influenza A Virus Polymerase Complex Mediate Mutagen Resistance and Replication Fidelity

RNA viruses exist as genetically diverse populations. This standing genetic diversity gives them the potential to adapt rapidly, evolve resistance to antiviral therapeutics, and evade immune responses. Viral mutants with altered mutation rates or mutational tolerance have provided insights into how genetic diversity arises and how it affects the behavior of RNA viruses. To this end, we identified variants within the polymerase complex of influenza virus that are able tolerate drug-mediated increases in viral mutation rates. We find that drug resistance is highly dependent on interactions among mutations in the polymerase complex. In contrast to other viruses, influenza virus counters the effect of higher mutation rates primarily by maintaining high levels of genome replication. These findings suggest the importance of maintaining large population sizes for viruses with high mutation rates and show that multiple proteins can affect both mutation rate and genome synthesis.

Lethal mutagenesis is a broad-spectrum antiviral strategy that employs mutagenic nucleoside analogs to exploit the high mutation rate and low mutational tolerance of many RNA viruses. Studies of mutagen-resistant viruses have identified determinants of replicative fidelity and the importance of mutation rate to viral population dynamics. We have previously demonstrated the effective lethal mutagenesis of influenza A virus using three nucleoside analogs as well as the virus’s high genetic barrier to mutagen resistance. Here, we investigate the mutagen-resistant phenotypes of mutations that were enriched in drug-treated populations. We find that PB1 T123A has higher replicative fitness than the wild type, PR8, and maintains its level of genome production during 5-fluorouracil (2,4-dihydroxy-5-fluoropyrimidine) treatment. Surprisingly, this mutagen-resistant variant also has an increased baseline rate of C-to-U and G-to-A mutations. A second drug-selected mutation, PA T97I, interacts epistatically with PB1 T123A to mediate high-level mutagen resistance, predominantly by limiting the inhibitory effect of nucleosides on polymerase activity. Consistent with the importance of epistatic interactions in the influenza virus polymerase, our data suggest that nucleoside analog resistance and replication fidelity are strain dependent. Two previously identified ribavirin {1-[(2R,3R,4S,5R)-3,4-dihydroxy-5-(hydroxymethyl)oxolan-2-yl]-1H-1,2,4-triazole-3-carboxamide} resistance mutations, PB1 V43I and PB1 D27N, do not confer drug resistance in the PR8 background, and the PR8-PB1 V43I polymerase exhibits a normal baseline mutation rate. Our results highlight the genetic complexity of the influenza A virus polymerase and demonstrate that increased replicative capacity is a mechanism by which an RNA virus can counter the negative effects of elevated mutation rates.

IMPORTANCE RNA viruses exist as genetically diverse populations. This standing genetic diversity gives them the potential to adapt rapidly, evolve resistance to antiviral therapeutics, and evade immune responses. Viral mutants with altered mutation rates or mutational tolerance have provided insights into how genetic diversity arises and how it affects the behavior of RNA viruses. To this end, we identified variants within the polymerase complex of influenza virus that are able to tolerate drug-mediated increases in viral mutation rates. We find that drug resistance is highly dependent on interactions among mutations in the polymerase complex. In contrast to other viruses, influenza virus counters the effect of higher mutation rates primarily by maintaining high levels of genome replication. These findings suggest the importance of maintaining large population sizes for viruses with high mutation rates and show that multiple proteins can affect both mutation rate and genome synthesis.

A novel twelve class fluctuation test reveals higher than expected mutation rates for influenza A viruses

Influenza virus’ low replicative fidelity contributes to its capacity for rapid evolution. Clonal sequencing and fluctuation tests have suggested that the influenza virus mutation rate is 2.7 × 10^–6 - 3.0 × 10^–5 substitutions per nucleotide per strand copied (s/n/r). However, sequencing assays are biased toward mutations with minimal fitness impacts and fluctuation tests typically investigate only a subset of all possible single nucleotide mutations. We developed a fluctuation test based on reversion to fluorescence in a set of virally encoded mutant green fluorescent proteins, which allowed us to measure the rates of selectively neutral mutations representative of the twelve different mutation types. We measured an overall mutation rate of 1.8 × 10^–4 s/n/r for PR8 (H1N1) and 2.5 × 10^–4 s/n/r for Hong Kong 2014 (H3N2) and a transitional bias of 2.7–3.6. Our data suggest that each replicated genome will have an average of 2–3 mutations and highlight the importance of mutational load in influenza virus evolution.

DOI:http://dx.doi.org/10.7554/eLife.26437.001

Known mutator alleles do not markedly increase mutation rate in clinical Saccharomyces cerevisiae strains

Natural selection has the potential to act on all phenotypes, including genomic mutation rate. Classic evolutionary theory predicts that in asexual populations, mutator alleles, which cause high mutation rates, can fix due to linkage with beneficial mutations. This phenomenon has been demonstrated experimentally and may explain the frequency of mutators found in bacterial pathogens. By contrast, in sexual populations, recombination decouples mutator alleles from beneficial mutations, preventing mutator fixation. In the facultatively sexual yeast Saccharomyces cerevisiae, segregating alleles of MLH1 and PMS1 have been shown to be incompatible, causing a high mutation rate when combined. These alleles had never been found together naturally, but were recently discovered in a cluster of clinical isolates. Here we report that the incompatible mutator allele combination only marginally elevates mutation rate in these clinical strains. Genomic and phylogenetic analyses provide no evidence of a historically elevated mutation rate. We conclude that the effect of the mutator alleles is dampened by background genetic modifiers. Thus, the relationship between mutation rate and microbial pathogenicity may be more complex than once thought. Our findings provide rare observational evidence that supports evolutionary theory suggesting that sexual organisms are unlikely to harbour alleles that increase their genomic mutation rate.

Population Heterogeneity in Mutation Rate Increases the Frequency of Higher-Order Mutants and Reduces Long-Term Mutational Load

Mutation rate is a crucial evolutionary parameter that has typically been treated as a constant in population genetic analyses. However, the propensity to mutate is likely to vary among co-existing individuals within a population, due to genetic polymorphisms, heterogeneous environmental influences, and random physiological fluctuations. We review the evidence for mutation rate heterogeneity and explore its consequences by extending classic population genetic models to allow an arbitrary distribution of mutation rate among individuals, either with or without inheritance. With this general new framework, we rigorously establish the effects of heterogeneity at various evolutionary timescales. In a single generation, variation of mutation rate about the mean increases the probability of producing zero or many simultaneous mutations on a genome. Over multiple generations of mutation and selection, heterogeneity accelerates the appearance of both deleterious and beneficial multi-point mutants. At mutation-selection balance, higher-order mutant frequencies are likewise boosted, while lower-order mutants exhibit subtler effects; nonetheless, population mean fitness is always enhanced. We quantify the dependencies on moments of the mutation rate distribution and selection coefficients, and clarify the role of mutation rate inheritance. While typical methods of estimating mutation rate will recover only the population mean, analyses assuming mutation rate is fixed to this mean could underestimate the potential for multi-locus adaptation, including medically relevant evolution in pathogenic and cancerous populations. We discuss the potential to empirically parameterize mutation rate distributions, which have to date hardly been quantified.

Development of an Influenza A Master Virus for Generating High-Growth Reassortants for A/Anhui/1/2013(H7N9) Vaccine Production in Qualified MDCK Cells

In 2013, the first case of human infection with an avian influenza A virus (H7N9) was reported in China, and the human infection with this virus has continued as of 2016. At the request of the WHO, we have successfully developed candidate reassortant vaccine virus using A/Anhui/1/2013 and the high egg-growth master virus A/PR/8/1934. Recent plans regarding influenza vaccine production include using cell-cultured systems in Japan and several other countries. However, egg-based vaccine viruses are not always suitable for cell-cultured vaccine production due to potential issues with growth, protein yield and antigenic stability. Therefore, in this study, we have developed a high-growth master virus (hg-PR8) adapted to qualified NIID-MDCK cells that are competent for vaccine production. The virus hg-PR8 was obtained after 20 serial passages of A/Puerto Rico/8/1934 (PR8) in NIID-MDCK cells. The viral titer of hg-PR8 was 10^8.6 plaque-forming units per milliliter (PFU/mL). Seven amino acid substitutions were identified in the amino acid sequences of PB2, PB1, PA, NA, M and NS of hg-PR8 compared to the sequence of the original PR8 (org-PR8) strain. The growth capacities of the reassortant viruses, which possess heterologous internal genes from hg-PR8 or org-PR8, indicated that the amino acid changes in PB2 and NS2 similarly affected growth capacity in NIID-MDCK cells. To assess the suitability of hg-PR8 as a master virus, we generated 6:2 reassortant viruses possessing the HA and NA segments from A/Anhui/1/2013 (H7N9) and the remaining segments from hg-PR8. The virus titers of the reassortant strains were 10⁷−10⁸ PFU/mL. The antigenicity of the viruses was stable during ten passages of the viruses in NIID-MDCK cells. In comparison with the egg-based reassortant vaccine viruses with identical HA and NA segments, the hg-PR8-based viruses showed 1.5- to 2-fold higher protein yields in NIID-MDCK cells.

Fidelity Variants and RNA Quasispecies

By now, it is well established that the error rate of the RNA-dependent RNA polymerase (RdRp) that replicates RNA virus genomes is a primary driver of the mutation frequencies observed in RNA virus populations-the basis for the RNA quasispecies. Over the last 10 years, a considerable amount of work has uncovered the molecular determinants of replication fidelity in this enzyme. The isolation of high- and low-fidelity variants for several RNA viruses, in an expanding number of viral families, provides evidence that nature has optimized the fidelity to facilitate genetic diversity and adaptation, while maintaining genetic integrity and infectivity. This chapter will provide an overview of what fidelity variants tell us about RNA virus biology and how they may be used in antiviral approaches.

Full Genome Characterization of Human Influenza A/H3N2 Isolates from Asian Countries Reveals a Rare Amantadine Resistance-Conferring Mutation and Novel PB1-F2 Polymorphisms

Influenza A viruses evolve at a high rate requiring continuous monitoring to maintain the efficacy of vaccines and antiviral drugs. We performed next generation sequencing analysis of 100 influenza A/H3N2 isolates collected in four Asian countries (Japan, Lebanon, Myanmar, and Vietnam) during 2012-2015. Phylogenetic analysis revealed several reassortment events leading to the circulation of multiple clades within the same season. This was particularly evident during the 2013 and 2013/2014 seasons. Importantly, our data showed that certain lineages appeared to be fitter and were able to persist into the following season. The majority of A/H3N2 viruses continued to harbor the M2-S31N mutation conferring amantadine-resistance. In addition, an S31D mutation in the M2-protein, conferring a similar level of resistance as the S31N mutation, was detected in three isolates obtained in Japan during the 2014/2015 season. None of the isolates possessed the NA-H274Y mutation conferring oseltamivir-resistance, though a few isolates were found to contain mutations at the catalytic residue 151 (D151A/G/N or V) of the NA protein. These variations did not alter the susceptibility to neuraminidase inhibitors and were not detected in the original clinical specimens, suggesting that they had been acquired during their passage in MDCK cells. Novel polymorphisms were detected in the PB1-F2 open-reading frame resulting in truncations in the protein of 24-34 aminoacids in length. Thus, this study has demonstrated the utility of monitoring the full genome of influenza viruses to allow the detection of the potentially fittest lineages. This enhances our ability to predict the strain(s) most likely to persist into the following seasons and predict the potential degree of vaccine match or mismatch with the seasonal influenza season for that year. This will enable the public health and clinical teams to prepare for any related healthcare burden, depending on whether the vaccine match is predicted to be good or poor for that season.

Coronavirus Host Range Expansion and Middle East Respiratory Syndrome Coronavirus Emergence: Biochemical Mechanisms and Evolutionary Perspectives

Coronaviruses have frequently expanded their host range in recent history, with two events resulting in severe disease outbreaks in human populations. Severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in 2003 in Southeast Asia and rapidly spread around the world before it was controlled by public health intervention strategies. The 2012 Middle East respiratory syndrome coronavirus (MERS-CoV) outbreak represents another prime example of virus emergence from a zoonotic reservoir. Here, we review the current knowledge of coronavirus cross-species transmission, with particular focus on MERS-CoV. MERS-CoV is still circulating in the human population, and the mechanisms governing its cross-species transmission have been only partially elucidated, highlighting a need for further investigation. We discuss biochemical determinants mediating MERS-CoV host cell permissivity, including virus spike interactions with the MERS-CoV cell surface receptor dipeptidyl peptidase 4 (DPP4), and evolutionary mechanisms that may facilitate host range expansion, including recombination, mutator alleles, and mutational robustness. Understanding these mechanisms can help us better recognize the threat of emergence for currently circulating zoonotic strains.

Expected Effect of Deleterious Mutations on Within-Host Adaptation of Pathogens

Adaptation is a common theme in both pathogen emergence, for example, in zoonotic cross-species transmission, and pathogen control, where adaptation might limit the effect of the immune response and antiviral treatment. When such evolution requires deleterious intermediate mutations, fitness ridges and valleys arise in the pathogen's fitness landscape. The effect of deleterious intermediate mutations on within-host pathogen adaptation is examined with deterministic calculations, appropriate for pathogens replicating in large populations with high error rates. The effect of deleterious intermediate mutations on pathogen adaptation is smaller than their name might suggest: when two mutations are required and each individual single mutation is fully deleterious, the pathogen can jump across the fitness valley by obtaining two mutations at once, leading to a proportion of adapted mutants that is 20-fold lower than that in the situation where the fitness of all mutants is neutral. The negative effects of deleterious intermediates are typically substantially smaller and outweighed by the fitness advantages of the adapted mutant. Moreover, requiring a specific mutation order has a substantially smaller effect on pathogen adaptation than the effect of all intermediates being deleterious. These results can be rationalized when the number of routes of mutation available to the pathogen is calculated, providing a simple approach to estimate the effect of deleterious mutations. The calculations discussed here are applicable when the effect of deleterious mutations on the within-host adaptation of pathogens is assessed, for example, in the context of zoonotic emergence, antigenic escape, and drug resistance.

IMPORTANCE

Adaptation is critical for pathogens after zoonotic transmission into a new host species or to achieve antigenic immune escape and drug resistance. Using a deterministic approach, the effects of deleterious intermediate mutations on pathogen adaptation were calculated while avoiding commonly made simplifications that do not apply to large pathogen populations replicating with high mutation rates. Perhaps unexpectedly, pathogen adaptation does not halt when the intermediate mutations are fully deleterious. The negative effects of deleterious mutations are substantially outweighed by the fitness gains of adaptation. To gain an understanding of the effect of deleterious mutations on pathogen adaptation, a simple approach that counts the number of routes available to the pathogen with and without deleterious intermediate mutations is introduced. This methodology enables a straightforward calculation of the proportion of the pathogen population that will cross a fitness valley or traverse a fitness ridge, without reverting to more complicated mathematical models.

Experimental evolution and the dynamics of genomic mutation rate modifiers

Because genes that affect mutation rates are themselves subject to mutation, mutation rates can be influenced by natural selection and other evolutionary forces. The population genetics of mutation rate modifier alleles has been a subject of theoretical interest for many decades. Here, we review experimental contributions to our understanding of mutation rate modifier dynamics. Numerous evolution experiments have shown that mutator alleles (modifiers that elevate the genomic mutation rate) can readily rise to high frequencies via genetic hitchhiking in non-recombining microbial populations. Whereas these results certainly provide an explanatory framework for observations of sporadically high mutation rates in pathogenic microbes and in cancer lineages, it is nonetheless true that most natural populations have very low mutation rates. This raises the interesting question of how mutator hitchhiking is suppressed or its phenotypic effect reversed in natural populations. Very little experimental work has addressed this question; with this in mind, we identify some promising areas for future experimental investigation.

Variation in RNA virus mutation rates across host cells

It is well established that RNA viruses exhibit higher rates of spontaneous mutation than DNA viruses and microorganisms. However, their mutation rates vary amply, from 10⁻⁶ to 10⁻⁴ substitutions per nucleotide per round of copying (s/n/r) and the causes of this variability remain poorly understood. In addition to differences in intrinsic fidelity or error correction capability, viral mutation rates may be dependent on host factors. Here, we assessed the effect of the cellular environment on the rate of spontaneous mutation of the vesicular stomatitis virus (VSV), which has a broad host range and cell tropism. Luria-Delbrück fluctuation tests and sequencing showed that VSV mutated similarly in baby hamster kidney, murine embryonic fibroblasts, colon cancer, and neuroblastoma cells (approx. 10⁻⁵ s/n/r). Cell immortalization through p53 inactivation and oxygen levels (1–21%) did not have a significant impact on viral replication fidelity. This shows that previously published mutation rates can be considered reliable despite being based on a narrow and artificial set of laboratory conditions. Interestingly, we also found that VSV mutated approximately four times more slowly in various insect cells compared with mammalian cells. This may contribute to explaining the relatively slow evolution of VSV and other arthropod-borne viruses in nature.

RNA viruses show high rates of spontaneous mutation, a feature that profoundly influences viral evolution, disease emergence, the appearance of drug resistances, and vaccine efficacy. However, RNA virus mutation rates vary substantially and the factors determining this variability remain poorly understood. Here, we investigated the effects of host factors on viral replication fidelity by measuring the viral mutation rate in different cell types and under various culturing conditions. To carry out these experiments we chose the vesicular stomatitis virus (VSV), an insect-transmitted mammalian RNA virus with an extremely wide cellular and host tropism. We found that the VSV replication machinery was robust to changes in cellular physiology driven by cell immortalization or shifts in temperature and oxygen levels. In contrast, VSV mutated significantly more slowly in insect cells than in mammalian cells, a finding may help us to understand why arthropod-borne viruses tend to evolve more slowly than directly transmitted viruses in nature.

Alphavirus mutator variants present host-specific defects and attenuation in mammalian and insect models

Arboviruses cycle through both vertebrates and invertebrates, which requires them to adapt to disparate hosts while maintaining genetic integrity during genome replication. To study the genetic mechanisms and determinants of these processes, we use chikungunya virus (CHIKV), a re-emerging human pathogen transmitted by the Aedes mosquito. We previously isolated a high fidelity (or antimutator) polymerase variant, C483Y, which had decreased fitness in both mammalian and mosquito hosts, suggesting this residue may be a key molecular determinant. To further investigate effects of position 483 on RNA-dependent RNA-polymerase (RdRp) fidelity, we substituted every amino acid at this position. We isolated novel mutators with decreased replication fidelity and higher mutation frequencies, allowing us to examine the fitness of error-prone arbovirus variants. Although CHIKV mutators displayed no major replication defects in mammalian cell culture, they had reduced specific infectivity and were attenuated in vivo. Unexpectedly, mutator phenotypes were suppressed in mosquito cells and the variants exhibited significant defects in RNA synthesis. Consequently, these replication defects resulted in strong selection for reversion during infection of mosquitoes. Since residue 483 is conserved among alphaviruses, we examined the analogous mutations in Sindbis virus (SINV), which also reduced polymerase fidelity and generated replication defects in mosquito cells. However, replication defects were mosquito cell-specific and were not observed in Drosophila S2 cells, allowing us to evaluate the potential attenuation of mutators in insect models where pressure for reversion was absent. Indeed, the SINV mutator variant was attenuated in fruit flies. These findings confirm that residue 483 is a determinant regulating alphavirus polymerase fidelity and demonstrate proof of principle that arboviruses can be attenuated in mammalian and insect hosts by reducing fidelity.

Chikungunya (CHIKV) is a re-emerging mosquito-borne virus that constitutes a major and growing human health burden. Like all RNA viruses, during viral replication CHIKV copies its genome using a polymerase that makes an average of one mistake per replication cycle. Therefore, a single virus generates millions of viral progeny that carry a multitude of distinct mutations in their genomes. In this study, we isolated CHIKV mutators (strains that make more errors than the wildtype virus), to study how higher mutation rates affect fitness in arthropod-borne viruses (arboviruses). CHIKV mutators have reduced virulence in mice and severe replication defects in Aedes mosquito cells. However, these replication defects result in selective pressure for reversion of mutators to a wildtype polymerase in mosquito hosts. To examine how mutators would behave in an insect model in absence of this genetic instability, we isolated mutators of a related virus, Sindbis virus (SINV). SINV mutators had no replication defect in fruit fly (Drosophila) cells, and a SINV mutator strain was stable and attenuated in fruit flies. This work shows proof of principle that arbovirus mutators can exhibit attenuation in both mammalian and insect hosts, and may remain a viable vaccine strategy.

Cell tropism predicts long-term nucleotide substitution rates of mammalian RNA viruses

The high rates of RNA virus evolution are generally attributed to replication with error-prone RNA-dependent RNA polymerases. However, these long-term nucleotide substitution rates span three orders of magnitude and do not correlate well with mutation rates or selection pressures. This substitution rate variation may be explained by differences in virus ecology or intrinsic genomic properties. We generated nucleotide substitution rate estimates for mammalian RNA viruses and compiled comparable published rates, yielding a dataset of 118 substitution rates of structural genes from 51 different species, as well as 40 rates of non-structural genes from 28 species. Through ANCOVA analyses, we evaluated the relationships between these rates and four ecological factors: target cell, transmission route, host range, infection duration; and three genomic properties: genome length, genome sense, genome segmentation. Of these seven factors, we found target cells to be the only significant predictors of viral substitution rates, with tropisms for epithelial cells or neurons (P<0.0001) as the most significant predictors. Further, one-tailed t-tests showed that viruses primarily infecting epithelial cells evolve significantly faster than neurotropic viruses (P<0.0001 and P<0.001 for the structural genes and non-structural genes, respectively). These results provide strong evidence that the fastest evolving mammalian RNA viruses infect cells with the highest turnover rates: the highly proliferative epithelial cells. Estimated viral generation times suggest that epithelial-infecting viruses replicate more quickly than viruses with different cell tropisms. Our results indicate that cell tropism is a key factor in viral evolvability.

RNA viruses are the fastest evolving human pathogens, making their treatment and control difficult. Compared to DNA viruses, RNA viruses replicate with much lower fidelity, which can explain why RNA viruses evolve significantly faster than most DNA viruses. However, there is tremendous variation among the evolutionary rates of different RNA viruses, which is not explained by variation in mutation rates. Here we present a survey of mammalian RNA virus rates of evolution, and a comprehensive comparison of these rates to different properties of virus genomic architecture and ecology. We found that cell tropism is the most significant predictor of long-term rates of mammalian RNA virus evolution. For instance, viruses targeting epithelial cells evolve significantly faster than viruses that target neurons. Our results provide mechanistic insight into why viruses that infect respiratory and gastrointestinal epithelia have been difficult to control.

Biological activities of 'noninfectious' influenza A virus particles

Only a small fraction of influenza A virus (IAV) particles within a viral population register as infectious by traditional infectivity assays. Despite constituting the most abundant product of influenza infection, the role that the 'noninfectious' particle fraction plays in the biology of the virus has largely been ignored. This review shines a light on this oft-ignored population by highlighting studies, both old and new, that describe the unique biological activities of these particles, and discussing what this population can tell us about the biology of IAV evolution and disease.

The effect of population bottlenecks on mutation rate evolution in asexual populations

In the absence of recombination, a mutator allele can spread through a population by hitchhiking with beneficial mutations that appear in its genetic background. Theoretical studies over the past decade have shown that the survival and fixation probability of beneficial mutations can be severely reduced by population size bottlenecks. Here, we use computational modelling and evolution experiments with the yeast S. cerevisiae to examine whether population bottlenecks can affect mutator dynamics in adapting asexual populations. In simulation, we show that population bottlenecks can inhibit mutator hitchhiking with beneficial mutations and are most effective at lower beneficial mutation supply rates. We then subjected experimental populations of yeast propagated at the same effective population size to three different bottleneck regimes and observed that the speed of mutator hitchhiking was significantly slower at smaller bottlenecks, consistent with our theoretical expectations. Our results, thus, suggest that bottlenecks can be an important factor in mutation rate evolution and can in certain circumstances act to stabilize or, at least, delay the progressive elevation of mutation rates in asexual populations. Additionally, our findings provide the first experimental support for the theoretically postulated effect of population bottlenecks on beneficial mutations and demonstrate the usefulness of studying mutator frequency dynamics for understanding the underlying dynamics of fitness-affecting mutations.

Correlation between mutation rate and genome size in riboviruses: mutation rate of bacteriophage Qβ

Genome sizes and mutation rates covary across all domains of life. In unicellular organisms and DNA viruses, they show an inverse relationship known as Drake’s rule. However, it is still unclear whether a similar relationship exists between genome sizes and mutation rates in RNA genomes. Coronaviruses, the RNA viruses with the largest genomes (∼30 kb), encode a proofreading 3′ exonuclease that allows them to increase replication fidelity. However, it is unknown whether, conversely, the RNA viruses with the smallest genomes tend to show particularly high mutation rates. To test this, we measured the mutation rate of bacteriophage Qβ, a 4.2-kb levivirus. Amber reversion-based Luria–Delbrück fluctuation tests combined with mutant sequencing gave an estimate of 1.4 × 10⁻⁴ substitutions per nucleotide per round of copying, the highest mutation rate reported for any virus using this method. This estimate was confirmed using a direct plaque sequencing approach and after reanalysis of previously published estimates for this phage. Comparison with other riboviruses (all RNA viruses except retroviruses) provided statistical support for a negative correlation between mutation rates and genome sizes. We suggest that the mutation rates of RNA viruses might be optimized for maximal adaptability and that the value of this optimum may in turn depend inversely on genome size.

Reassortment complements spontaneous mutation in influenza A virus NP and M1 genes to accelerate adaptation to a new host

Influenza A virus (IAV) infects a remarkably wide variety of avian and mammalian hosts. Evolution finely hones IAV genes to optimally infect and be transmitted in a particular host species. Sporadically, IAV manages to jump between species, introducing novel antigenic strains into the new host population that wreak havoc until herd immunity develops. IAV adaptation to new hosts typically involves reassortment of IAV gene segments from coinfecting virus strains adapted to different hosts in conjunction with multiple adaptive mutations in the various IAV genes. To better understand host adaptation between mammalian species in real time, we passaged mouse-adapted A/PR8/34 (PR8) in guinea pigs. Guinea pigs, unlike mice, support spontaneous and robust IAV transmission. For some IAV strains, including PR8, adaptation is required for a virus to attain transmissibility, providing an opportunity to understand the evolution of transmissibility in guinea pigs. Multiple guinea pig-adapted PR8 mutants generated by serial nasal wash passaging in independent lines replicated more efficiently and were transmitted by cocaging. All transmissible variants possessed one of two nonsynonymous mutations in M1, either alone or in combination with mutations in PB2, HA, NP, or NA. Rapid reassortment between independently selected variants combined beneficial mutations in NP and M1 to form the fittest virus capable of being transmitted. These findings provide further insight into genetic determinants in NP and M1 involved in PR8 IAV adaptation to be transmitted in a new host and clearly show the benefit of a segmented genome in rapidly generating optimal combinations of mutations in IAV evolution.

Influenza virus subpopulations: exchange of lethal H5N1 virus NS for H1N1 virus NS triggers de novo generation of defective-interfering particles and enhances interferon-inducing particle efficiency

Reassortment of influenza A viruses is known to affect viability, replication efficiency, antigenicity, host range, and virulence, and can generate pandemic strains. In this study, we demonstrated that the specific exchange of the NS gene segment from highly pathogenic A/HK/156/97 (H5N1) [E92 or E92D NS1] virus for the cognate NS gene segment of A/PR/834(H1N1) [D92 NS1] virus did not cause a significant change in the sizes of infectious particle subpopulations. However, it resulted in 2 new phenotypic changes: (1) de novo generation of large subpopulations of defective-interfering particles (DIPs); and (2) enhancement of interferon (IFN)-inducing particle efficiency leading to an order of magnitude or higher quantum (peak) yield of IFN in both avian and mammalian cells. These changes were attributed to loss of function of the H5N1-NS gene products. Most notably, the NS exchange obliterated the usual IFN-induction-suppressing capacity associated with expression of full-size NS1 proteins, and hence functionally mimicked deletions in the NS1 gene. The loss of NS1-mediated suppression of IFN induction, de novo generation of DIPs, and the concomitant enhancement of IFN-inducing particle efficiency suggest that in an attenuated background, the H5N1-NS could be used to formulate a self-adjuvanting live attenuated influenza vaccine similar to viruses with deletions in the NS1 gene.

Selection of variant viruses during replication and transmission of H7N1 viruses in chickens and turkeys

The influence of different glycosylation patterns of the haemagglutinin glycoprotein of H7N1 avian influenza viruses on virus replication in vivo was examined. Experimental infection of chickens and turkeys was carried out with H7N1 avian influenza viruses with alternative sites of glycosylation in the haemagglutinin and infected birds were sampled daily by swabbing the buccal and cloacal cavities. cDNAs of the HA1 coding region of the HA gene were prepared from the swabs and cloned into plasmids. Sequencing multiple plasmids made from individual swabs taken over the period of virus shedding showed that viruses with specific patterns of glycosylation near the receptor binding site were stable when birds were infected with a single variant, but when presented with a mixed population of viruses encoding differing patterns of glycosylation a specific variant was rapidly selected in the infected host.

Interrelationship between HIV-1 fitness and mutation rate

Differences in replication fidelity, as well as mutator and antimutator strains, suggest that virus mutation rates are heritable and prone to natural selection. Human immunodeficiency virus type 1 (HIV-1) has many distinct advantages for the study of mutation rate optimization given the wealth of structural and biochemical data on HIV-1 reverse transcriptase (RT) and mutants. In this study, we conducted parallel analyses of mutation rate and viral fitness. In particular, a panel of 10 RT mutants-most having drug resistance phenotypes-was analyzed for their effects on viral fidelity and fitness. Fidelity differences were measured using single-cycle vector assays, while fitness differences were identified using ex vivo head-to-head competition assays. As anticipated, virus mutants possessing either higher or lower fidelity had a corresponding loss in fitness. While the virus panel was not chosen randomly, it is interesting that it included more viruses possessing a mutator phenotype rather than viruses possessing an antimutator phenotype. These observations provide the first description of an interrelationship between HIV-1 fitness and mutation rate and support the conclusion that mutator and antimutator phenotypes correlate with reduced viral fitness. In addition, the findings here help support a model in which fidelity comes at a cost of replication kinetics and may help explain why retroviruses like HIV-1 and RNA viruses maintain replication fidelity near the extinction threshold.

Genomic and protein structural maps of adaptive evolution of human influenza A virus to increased virulence in the mouse

Adaptive evolution is characterized by positive and parallel, or repeated selection of mutations. Mouse adaptation of influenza A virus (IAV) produces virulent mutants that demonstrate positive and parallel evolution of mutations in the hemagglutinin (HA) receptor and non-structural protein 1 (NS1) interferon antagonist genes. We now present a genomic analysis of all 11 genes of 39 mouse adapted IAV variants from 10 replicate adaptation experiments. Mutations were mapped on the primary and structural maps of each protein and specific mutations were validated with respect to virulence, replication, and RNA polymerase activity. Mouse adapted (MA) variants obtained after 12 or 20–21 serial infections acquired on average 5.8 and 7.9 nonsynonymous mutations per genome of 11 genes, respectively. Among a total of 115 nonsynonymous mutations, 51 demonstrated properties of natural selection including 27 parallel mutations. The greatest degree of parallel evolution occurred in the HA receptor and ribonucleocapsid components, polymerase subunits (PB1, PB2, PA) and NP. Mutations occurred in host nuclear trafficking factor binding sites as well as sites of virus-virus protein subunit interaction for NP, NS1, HA and NA proteins. Adaptive regions included cap binding and endonuclease domains in the PB2 and PA polymerase subunits. Four mutations in NS1 resulted in loss of binding to the host cleavage and polyadenylation specificity factor (CPSF30) suggesting that a reduction in inhibition of host gene expression was being selected. The most prevalent mutations in PB2 and NP were shown to increase virulence but differed in their ability to enhance replication and demonstrated epistatic effects. Several positively selected RNA polymerase mutations demonstrated increased virulence associated with >300% enhanced polymerase activity. Adaptive mutations that control host range and virulence were identified by their repeated selection to comprise a defined model for studying IAV evolution to increased virulence in the mouse.

"Boom" and "Bust" cycles in virus growth suggest multiple selective forces in influenza a evolution

Background

Influenza A virus evolution in humans is driven at least in part by mutations allowing the virus to escape antibody neutralization. Little is known about the evolution of influenza in birds, a major reservoir of influenza A.

Methods

Neutralizing polyclonal antiserum was raised in chicken against reassortant influenza virus, CalX, bearing the hemagglutinin (HA) and neuraminidase (NA) of A/California/7/2004 [H3N2]. CalX was serially passaged in the presence of anti-CalX polyclonal IgY to derive viruses capable of growth in the presence of antibody.

Results

Polyclonal chicken antibody neutralized both HA activity and infection by CalX, but had no effect on a strain bearing an earlier human H3 and an irrelevant neuraminidase (A/Memphis/71-Bellamy/42 [H3N1]). Surprisingly, most of the antibody-resistant viruses were still at least partially sensitive to neutralization of HA activity and viral infection. Although mutant HA genes bearing changes that might affect antibody neutralization were identified, the vast majority of HA sequences obtained were identical to wild type, and no individual mutant sequence was found in more than one passage, suggesting that those mutations that were observed did not confer sufficient selective advantage to come to dominate the population. Different passages yielded infectious foci of varying size and plaques of varying size and morphology. Yields of infectious virus and relative frequency of different morphologies changed markedly from passage to passage. Sequences of bulk, uncloned PCR products from antibody-resistant passages indicated changes in the PB2 and PA proteins with respect to the wild type virus.

Conclusions

Each antibody-selected passage consisted of a variety of different cocirculating populations, rather than pure populations of virus able to escape antibody by changes in antibody epitopes. The ability to escape antibody is apparently due to changes in genes encoding the viral polymerase complex, probably resulting in more robust viral replication, allowing the few virus particles not completely neutralized by antibody to rapidly produce large numbers of progeny. Our data suggest that the relative success of an individual variant may depend on both its own gain and loss of fitness, as well as that of its cocirculating variants.

Role of evolved host breadth in the initial emergence of an RNA virus

Understanding how evolution promotes pathogen emergence would aid disease management, and prediction of future host shifts. Increased pathogen infectiousness of different hosts may occur through direct selection, or fortuitously via indirect selection. However, it is unclear which type of selection tends to produce host breadth promoting pathogen emergence. We predicted that direct selection for host breadth should foster emergence by causing higher population growth on new hosts, lower among-population variance in growth on new hosts, and lower population variance in growth across new hosts. We tested the predictions using experimentally evolved vesicular stomatitis virus populations, containing groups of host-use specialists, directly selected generalists, and indirectly selected generalists. In novel-host challenges, viruses directly selected for generalism showed relatively higher or equivalent host growth, lower among-population variance in host growth, and lower population variance in growth across hosts. Thus, two of three outcomes supported our prediction that directly selected host breadth should favor host colonization. Also, we observed that indirectly selected generalists were advantaged over specialist viruses, indicating that fortuitous changes in host breadth may also promote emergence. We discuss evolution of phenotypic plasticity versus environmental robustness in viruses, virus avoidance of extinction, and surveillance of pathogen niche breadth to predict future likelihood of emergence.

Competition between high- and higher-mutating strains of Escherichia coli

Experimental studies have shown that a mutator allele can readily hitchhike to fixation with beneficial mutations in an asexual population having a low, wild-type mutation rate. Here, we show that a genotype bearing two mutator alleles can supplant a population already fixed for one mutator allele. Our results provide experimental support for recent theory predicting that mutator alleles will tend to accumulate in asexual populations by hitchhiking with beneficial mutations, causing an ever-higher genomic mutation rate.

Pathways to extinction: beyond the error threshold

Since the introduction of the quasispecies and the error catastrophe concepts for molecular evolution by Eigen and their subsequent application to viral populations, increased mutagenesis has become a common strategy to cause the extinction of viral infectivity. Nevertheless, the high complexity of virus populations has shown that viral extinction can occur through several other pathways apart from crossing an error threshold. Increases in the mutation rate enhance the appearance of defective forms and promote the selection of mechanisms that are able to counteract the accelerated appearance of mutations. Current models of viral evolution take into account more realistic scenarios that consider compensatory and lethal mutations, a highly redundant genotype-to-phenotype map, rough fitness landscapes relating phenotype and fitness, and where phenotype is described as a set of interdependent traits. Further, viral populations cannot be understood without specifying the characteristics of the environment where they evolve and adapt. Altogether, it turns out that the pathways through which viral quasispecies go extinct are multiple and diverse.

Viral mutation rates

Accurate estimates of virus mutation rates are important to understand the evolution of the viruses and to combat them. However, methods of estimation are varied and often complex. Here, we critically review over 40 original studies and establish criteria to facilitate comparative analyses. The mutation rates of 23 viruses are presented as substitutions per nucleotide per cell infection (s/n/c) and corrected for selection bias where necessary, using a new statistical method. The resulting rates range from 10(-8) to 10(-6) s/n/c for DNA viruses and from 10(-6) to 10(-4) s/n/c for RNA viruses. Similar to what has been shown previously for DNA viruses, there appears to be a negative correlation between mutation rate and genome size among RNA viruses, but this result requires further experimental testing. Contrary to some suggestions, the mutation rate of retroviruses is not lower than that of other RNA viruses. We also show that nucleotide substitutions are on average four times more common than insertions/deletions (indels). Finally, we provide estimates of the mutation rate per nucleotide per strand copying, which tends to be lower than that per cell infection because some viruses undergo several rounds of copying per cell, particularly double-stranded DNA viruses. A regularly updated virus mutation rate data set will be available at www.uv.es/rsanjuan/virmut.

Genetic and pathobiologic characterization of pandemic H1N1 2009 influenza viruses from a naturally infected swine herd

Since its initial identification in Mexico and the United States, concerns have been raised that the novel H1N1 influenza virus might cause a pandemic of severity comparable to that of the 1918 pandemic. In late April 2009, viruses phylogenetically related to pandemic H1N1 influenza virus were isolated from an outbreak on a Canadian pig farm. This outbreak also had epidemiological links to a suspected human case. Experimental infections carried out in pigs using one of the swine isolates from this outbreak and the human isolate A/Mexico/InDRE4487/2009 showed differences in virus recovery from the lower respiratory tract. Virus was consistently isolated from the lungs of pigs infected with A/Mexico/InDRE4487/2009, while only one pig infected with A/swine/Alberta/OTH-33-8/2008 yielded live virus from the lung, despite comparable amounts of viral RNA and antigen in both groups of pigs. Clinical disease resembled other influenza virus infections in swine, albeit with somewhat prolonged virus antigen detection and delayed viral-RNA clearance from the lungs. There was also a noteworthy amount of genotypic variability among the viruses isolated from the pigs on the farm. This, along with the somewhat irregular pathobiological characteristics observed in experimentally infected animals, suggests that although the virus may be of swine origin, significant viral evolution may still be ongoing.

Genome packaging in influenza A virus

The negative-sense RNA genome of influenza A virus is composed of eight segments, which encode 12 proteins between them. At the final stage of viral assembly, these genomic virion (v)RNAs are incorporated into the virion as it buds from the apical plasma membrane of the cell. Genome segmentation confers evolutionary advantages on the virus, but also poses a problem during virion assembly as at least one copy of each of the eight segments is required to produce a fully infectious virus particle. Historically, arguments have been presented in favour of a specific packaging mechanism that ensures incorporation of a full genome complement, as well as for an alternative model in which segments are chosen at random but packaged in sufficient numbers to ensure that a reasonable proportion of virions are viable. The question has seen a resurgence of interest in recent years leading to a consensus that the vast majority of virions contain no more than eight segments and that a specific mechanism does indeed function to select one copy of each vRNA. This review summarizes work leading to this conclusion. In addition, we describe recent progress in identifying the specific packaging signals and discuss likely mechanisms by which these RNA elements might operate.

Host restriction of avian influenza viruses at the level of the ribonucleoproteins

Although transmission of avian influenza viruses to mammals, particularly humans, has been repeatedly documented, adaptation and sustained transmission in the new host is a rare event that in the case of humans may result in pandemics. Host restriction involves multiple genetic determinants. Among the known determinants of host range, key determinants have been identified on the genes coding for the nucleoprotein and polymerase proteins that, together with the viral RNA segments, form the ribonucleoproteins (RNPs). The RNP genes form host-specific lineages and harbor host-associated genetic signatures. The functional significance of these determinants has been studied by reassortment and reverse genetics experiments, underlining the influence of the global genetic context. In some instances the molecular mechanisms have been approached, pointing to the importance of the polymerase activity and interaction with cellular host factors. Better knowledge of determinants of host restriction will allow monitoring of the pandemic potential of avian influenza viruses.

Theory of lethal mutagenesis for viruses

Mutation is the basis of adaptation. Yet, most mutations are detrimental, and elevating mutation rates will impair a population's fitness in the short term. The latter realization has led to the concept of lethal mutagenesis for curing viral infections, and work with drugs such as ribavirin has supported this perspective. As yet, there is no formal theory of lethal mutagenesis, although reference is commonly made to Eigen's error catastrophe theory. Here, we propose a theory of lethal mutagenesis. With an obvious parallel to the epidemiological threshold for eradication of a disease, a sufficient condition for lethal mutagenesis is that each viral genotype produces, on average, less than one progeny virus that goes on to infect a new cell. The extinction threshold involves an evolutionary component based on the mutation rate, but it also includes an ecological component, so the threshold cannot be calculated from the mutation rate alone. The genetic evolution of a large population undergoing mutagenesis is independent of whether the population is declining or stable, so there is no runaway accumulation of mutations or genetic signature for lethal mutagenesis that distinguishes it from a level of mutagenesis under which the population is maintained. To detect lethal mutagenesis, accurate measurements of the genome-wide mutation rate and the number of progeny per infected cell that go on to infect new cells are needed. We discuss three methods for estimating the former. Estimating the latter is more challenging, but broad limits to this estimate may be feasible.

Structure and Function of RNA Replication

Contrary to their host cells, many viruses contain RNA as genetic material and hence encode an RNA-dependent RNA polymerase to replicate their genomes. This review discusses the present status of our knowledge on the structure of these enzymes and the mechanisms of RNA replication. The simplest viruses encode only the catalytic subunit of the replication complex, but other viruses also contribute a variable number of ancillary factors. These and other factors provided by the host cell play roles in the specificity and affinity of template recognition and the assembly of the replication complex. Usually, these host factors are involved in protein synthesis or RNA modification in the host cell, but they play roles in remodeling RNA-RNA, RNA-protein, and protein-protein interactions during virus RNA replication. Furthermore, viruses take advantage of and modify previous cell structural elements, frequently membrane vesicles, for the formation of RNA replication complexes.