Permissive secondary mutations enable the evolution of influenza oseltamivir resistance

Antiviral susceptibility profile of influenza A viruses; keep an eye on immunocompromised patients under prolonged treatment

There was an increase in severe and fatal influenza cases in Greece during the 2011-2015 post-pandemic period. To investigate causality, we determined neuraminidase (NA) inhibitor susceptibility and resistance-conferring NA and hemagglutinin (HA) mutations in circulating influenza type A viruses during the pandemic (2009-2010) and post-pandemic periods in Greece. One hundred thirty-four influenza A(H1N1)pdm09 and 95 influenza A(H3N2) viruses submitted to the National Influenza Reference Laboratory of Southern Greece were tested for susceptibility to oseltamivir and zanamivir. Antiviral resistance was assessed by neuraminidase sequence analysis, as well as the fluorescence-based 50 % inhibitory concentration (IC₅₀) method. Five influenza A(H1N1)pdm09 viruses (2.2 %) showed significantly reduced inhibition by oseltamivir (average IC₅₀ 300.60nM vs. 1.19nM) by Gaussian kernel density plot analysis. These viruses were isolated from immunocompromised patients and harbored the H275Y oseltamivir resistance-conferring NA substitution. All A(H1N1)pdm09 viruses were zanamivir-susceptible, and all A(H3N2) viruses were susceptible to both drugs. Oseltamivir-resistant viruses did not form a distinct cluster by phylogenetic analysis. Permissive mutations were detected in immunogenic and non immunogenic NA regions of both oseltamivir- resistant and susceptible viruses in the post-pandemic seasons. Several amino acid substitutions in the HA1 domain of the HA gene of post-pandemic viruses were identified. This study indicated low resistance to NAIs among tested influenza viruses. Antiviral resistance emerged only in immunocompromised patients under long-term oseltamivir treatment. Sequential sample testing in this vulnerable group of patients is recommended to characterise resistance or reinfection and viral evolution.

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.

An experimental evaluation of drug-induced mutational meltdown as an antiviral treatment strategy

The rapid evolution of drug resistance remains a critical public health concern. The treatment of influenza A virus (IAV) has proven particularly challenging, due to the ability of the virus to develop resistance against current antivirals and vaccines. Here, we evaluate a novel antiviral drug therapy, favipiravir, for which the mechanism of action in IAV involves an interaction with the viral RNA-dependent RNA polymerase resulting in an effective increase in the viral mutation rate. We used an experimental evolution framework, combined with novel population genetic method development for inference from time-sampled data, to evaluate the effectiveness of favipiravir against IAV. Evaluating whole genome polymorphism data across 15 time points under multiple drug concentrations and in controls, we present the first evidence for the ability of IAV populations to effectively adapt to low concentrations of favipiravir. In contrast, under high concentrations, we observe population extinction, indicative of mutational meltdown. We discuss the observed dynamics with respect to the evolutionary forces at play and emphasize the utility of evolutionary theory to inform drug development.

Risk of resistant avian influenza A virus in wild waterfowl as a result of environmental release of oseltamivir

Oseltamivir is the best available anti-influenza drug and has therefore been stockpiled worldwide in large quantities as part of influenza pandemic preparedness planning. The active metabolite oseltamivir carboxylate (OC) is stable and is not removed by conventional sewage treatment. Active OC has been detected in river water at concentrations up to 0.86 µg/L. Although the natural reservoir hosts of influenza A virus (IAV) are wild waterfowl that reside in aquatic environments, the ecologic risks associated with environmental OC release and its potential to generate resistant viral variants among wild birds has largely been unknown. However, in recent years a number of in vivo mallard (Anas platyrhynchos) studies have been conducted regarding the potential of avian IAVs to become resistant to OC in natural reservoir birds if these are drug exposed. Development of resistance to OC was observed both in Group 1 (N1) and Group 2 (N2, N9) neuraminidase subtypes, when infected ducks were exposed to OC at concentrations between 0.95 and 12 µg/L in their water. All resistant variants maintained replication and transmission between ducks during drug exposure. In an A(H1N1)/H274Y virus, the OC resistance mutation persisted without selective drug pressure, demonstrating the potential of an IAV with a permissive genetic background to acquire and maintain OC resistance, potentially allowing circulation of the resistant variant among wild birds. The experimental studies have improved the appreciation of the risks associated with the environmental release of OC related to resistance development of avian IAVs among wild birds. Combined with knowledge of efficient methods for improved sewage treatment, the observations warrant implementation of novel efficient wastewater treatment methods, rational use of anti-influenza drugs, and improved surveillance of IAV resistance in wild birds.

Platform for determining the inhibition profile of neuraminidase inhibitors in an influenza virus N1 background

Efforts to develop novel neuraminidase inhibitors (NAIs) for the treatment of influenza are ongoing. Novel NAIs should in particular be also effective against seasonal and/or pandemic N1 that carry a H274Y or N294S substitution (N2 numbering), which are most commonly linked to oseltamivir resistance. Here we report a platform for profiling the efficacy of novel NAIs in the N1 genetic background of influenza A virus. Employing reverse genetics, a set of influenza virus variants containing an amino acid substitution associated with oseltamivir resistance in N1 isolates (H274Y, N294S, Y155H or Q136L) was generated. In parallel, so far unreported mutations of I427 (I427Q and I427M) were investigated. These possibly interfere with the side chain orientation of R371 and alter the binding affinity of most relevant NAIs. The profiling platform was validated with both oseltamivir and zanamivir and exemplarily applied to three analogs with differing decorations at positions 4 and 5. Besides confirming the inhibition profile of zanamivir and oseltamivir, the distinct effect of I427Q/M on the activity of both NAIs was shown. For 5-amidino and 5-guanidino analogs of oseltamivir a significantly stronger inhibition of virus variants carrying a NA-H274Y was confirmed, and additionally shown for NA-N294S and NA-Y155H substitutions as compared to the parent compound. Hence, the herein presented profiling platform is a valid tool for defining the inhibition profile of novel NAIs in the N1 background.

Mutations at highly conserved residues in influenza A(H1N1)pdm09 virus affect neuraminidase activity

Influenza virus neuraminidase (NA) plays a pivotal role during viral growth since its sialidase activity allows the efficient release of nascent virions from infected cells. Consequently, mutations in the NA catalytic site affecting sialic acid (SA) cleavage may influence the biological properties of influenza viruses. This study reports two amino acid substitutions (N386K and P431S) in the NA of the influenza A(H1N1)pdm09 virus that emerged in 2009 in Mexico. The NA sialidase activity to cleave SA-like substrates, and viral growth were examined and the mutant viruses had various deficiencies. NA mutations N386K and P431S together or separately, and in the presence or absence of H275Y were further evaluated using recombinant influenza A/California/04/2009 (pH1N1) viruses containing single, double, or triple mutations. Viral growth was reduced in the presence of mutation P431S alone or combined with N386K and/or H275Y. Substrates hydrolysis was reduced when recombinant pH1N1 viruses were analyzed by NA inhibitory assays. Moreover, elution assays with guinea pig red blood cells indicated an unbalanced hemagglutinin (HA):NA functionality. Altogether, our data underline the functional significance of mutations at highly conserved sites in influenza virus NA glycoprotein and the occurrence of permissive mutations to compensate virus viability in vitro.

Emergence and spread of drug resistant influenza: A two-population game theoretical model

Background

The potential for emergence of antiviral drug resistance during influenza pandemics has raised great concern for public health. Widespread use of antiviral drugs is a significant factor in producing resistant strains. Recent studies show that some influenza viruses may gain antiviral drug resistance without a fitness penalty. This creates the possibility of strategic interaction between populations considering antiviral drug use strategies.

Methods

To explain why, we develop and analyze a classical 2-player game theoretical model where each player chooses from a range of possible rates of antiviral drug use, and payoffs are derived as a function of final size of epidemic with the regular and mutant strain. Final sizes are derived from a stochastic compartmental epidemic model that captures transmission within each population and between populations, and the stochastic emergence of antiviral drug resistance. High treatment levels not only increase the spread of the resistant strain in the subject population but also affect the other population by increasing the density of the resistant strain infectious individuals due to travel between populations.

Results

We found two Nash equilibria where both populations treat at a high rate, or both treat at a low rate. Hence the game theoretical analysis predicts that populations will not choose different treatment strategies than other populations, under these assumptions. The populations may choose to cooperate by maintaining a low treatment rate that does not increase the incidence of mutant strain infections or cause case importations to the other population. Alternatively, if one population is treating at a high rate, this will generate a large number of mutant infections that spread to the other population, in turn incentivizing that population to also treat at a high rate. The prediction of two separate Nash equilibria is robust to the mutation rate and the effectiveness of the drug in preventing transmission, but it is sensitive to the volume of travel between the two populations.

Conclusions

Model-based evaluations of antiviral influenza drug use during a pandemic usually consider populations in isolation from one another, but our results show that strategic interactions could strongly influence a population's choice of antiviral drug use policy. Furthermore, the high treatment rate Nash equilibrium has the potential to become socially suboptimal (i.e. non-Pareto optimal) under model assumptions that might apply under other conditions. Because of the need for players to coordinate their actions, we conclude that communication and coordination between jurisdictions during influenza pandemics is a priority, especially for influenza strains that do not evolve a fitness penalty under antiviral drug resistance.

Flow Cytometric and Cytokine ELISpot Approaches To Characterize the Cell-Mediated Immune Response in Ferrets following Influenza Virus Infection

Influenza virus infections represent a significant socioeconomic and public health burden worldwide. Although ferrets are considered by many to be ideal for modeling human responses to influenza infection and vaccination, efforts to understand the cellular immune response have been severely hampered by a paucity of standardized procedures and reagents. In this study, we developed flow cytometric and T cell enzyme-linked immunosorbent spot (ELISpot) approaches to characterize the leukocyte composition and antigen-specific T cell response within key lymphoid tissues following influenza virus infection in ferrets. Through a newly designed and implemented set of serological reagents, we used multiparameter flow cytometry to directly quantify the frequency of CD4(+) and CD8(+) T cells, Ig(+) B cells, CD11b(+) myeloid-derived cells, and major histocompatibility complex (MHC) class II-positive antigen-presenting cells (APCs) both prior to and after intranasal infection with A/California/04/09 (H1N1). We found that the leukocyte composition was altered at 10 days postinfection, with notable gains in the frequency of T cells and myeloid cells within the draining lymph node. Furthermore, these studies revealed that the antigen specificity of influenza virus-reactive CD4 and CD8 T cells was very broad, with recognition of the viral HA, NA, M1, NS1, and NP proteins, and that total reactivity to influenza virus postinfection represented approximately 0.1% of the circulating peripheral blood mononuclear cells (PBMC). Finally, we observed distinct patterns of reactivity between individual animals, suggesting heterogeneity at the MHC locus in ferrets within commercial populations, a finding of considerable interest in efforts to move the ferret model forward for influenza vaccine and challenge studies.

IMPORTANCE

Ferrets are an ideal animal model to study transmission, diseases, and vaccine efficacies of respiratory viruses because of their close anatomical and physiological resemblances to humans. However, a lack of reagents has limited our understanding of the cell-mediated immune response following infection and vaccination. In this study, we used cross-reactive and ferret-specific antibodies to study the leukocyte composition and antigen-specific CD4 and CD8 T cell responses following influenza A/California/04/09 (H1N1) virus infection. These studies revealed strikingly distinct patterns of reactivity between CD4 and CD8 T cells, which were overlaid with differences in protein-specific responses between individual animals. Our results provide a first, in-depth look at the T cell repertoire in response to influenza infection and suggest that there is considerable heterogeneity at the MHC locus, which is akin to that in humans and an area of intense research interest.

A Novel Endonuclease Inhibitor Exhibits Broad-Spectrum Anti-Influenza Virus Activity In Vitro

Antiviral drugs are important in preventing and controlling influenza, particularly when vaccines are ineffective or unavailable. A single class of antiviral drugs, the neuraminidase inhibitors (NAIs), is recommended for treating influenza. The limited therapeutic options and the potential risk of antiviral resistance are driving the search for additional small-molecule inhibitors that act on influenza virus proteins. The acid polymerase (PA) of influenza viruses is a promising target for new antivirals because of its essential role in initiating virus transcription. Here, we characterized a novel compound, RO-7, identified as a putative PA endonuclease inhibitor. RO-7 was effective when added before the cessation of genome replication, reduced polymerase activity in cell-free systems, and decreased relative amounts of viral mRNA and genomic RNA during influenza virus infection. RO-7 specifically inhibited the ability of the PA endonuclease domain to cleave a nucleic acid substrate. RO-7 also inhibited influenza A viruses (seasonal and 2009 pandemic H1N1 and seasonal H3N2) and B viruses (Yamagata and Victoria lineages), zoonotic viruses (H5N1, H7N9, and H9N2), and NAI-resistant variants in plaque reduction, yield reduction, and cell viability assays in Madin-Darby canine kidney (MDCK) cells with nanomolar to submicromolar 50% effective concentrations (EC50s), low toxicity, and favorable selective indices. RO-7 also inhibited influenza virus replication in primary normal human bronchial epithelial cells. Overall, RO-7 exhibits broad-spectrum activity against influenza A and B viruses in multiple in vitro assays, supporting its further characterization and development as a potential antiviral agent for treating influenza.

The Mutational Robustness of Influenza A Virus

A virus’ mutational robustness is described in terms of the strength and distribution of the mutational fitness effects, or MFE. The distribution of MFE is central to many questions in evolutionary theory and is a key parameter in models of molecular evolution. Here we define the mutational fitness effects in influenza A virus by generating 128 viruses, each with a single nucleotide mutation. In contrast to mutational scanning approaches, this strategy allowed us to unambiguously assign fitness values to individual mutations. The presence of each desired mutation and the absence of additional mutations were verified by next generation sequencing of each stock. A mutation was considered lethal only after we failed to rescue virus in three independent transfections. We measured the fitness of each viable mutant relative to the wild type by quantitative RT-PCR following direct competition on A549 cells. We found that 31.6% of the mutations in the genome-wide dataset were lethal and that the lethal fraction did not differ appreciably between the HA- and NA-encoding segments and the rest of the genome. Of the viable mutants, the fitness mean and standard deviation were 0.80 and 0.22 in the genome-wide dataset and best modeled as a beta distribution. The fitness impact of mutation was marginally lower in the segments coding for HA and NA (0.88 ± 0.16) than in the other 6 segments (0.78 ± 0.24), and their respective beta distributions had slightly different shape parameters. The results for influenza A virus are remarkably similar to our own analysis of CirSeq-derived fitness values from poliovirus and previously published data from other small, single stranded DNA and RNA viruses. These data suggest that genome size, and not nucleic acid type or mode of replication, is the main determinant of viral mutational fitness effects.

Like other RNA viruses, influenza virus has a very high mutation rate. While high mutation rates may increase the rate at which influenza virus will adapt to a new host, acquire a new route of transmission, or escape from host immune surveillance, data from model systems suggest that most new viral mutations are either lethal or highly detrimental. Mutational robustness refers to the ability of a virus to tolerate, or buffer, these mutations. The mutational robustness of a virus will determine which mutations are maintained in a population and may have a greater impact on viral evolution than mutation rate. We defined the mutational robustness of influenza A virus by measuring the fitness of a large number of viruses, each with a single point mutation. We found that the overall robustness of influenza was similar to that of poliovirus and other viruses of similar size. Interestingly, mutations appeared to be more easily accommodated in hemagglutinin and neuraminidase than elsewhere in the genome. This work will inform models of influenza evolution at the global and molecular scale.

Positive Selection Drives Evolution at the Host-Filovirus Interaction Surface

Filovirus infection is mediated by engagement of the surface-exposed glycoprotein (GP) by its cellular receptor, NPC1 (Niemann-Pick C1). Two loops in the C domain of NPC1 (NPC1-C) bind filovirus GP. Herein, we show that filovirus GP and NPC1-C evolve under mutual selective pressure. Analysis of a large mammalian phylogeny indicated that strong functional/structural constraints limit the NPC1 sequence space available for adaptive change and most sites at the contact interface with GP are under negative selection. These constraints notwithstanding, we detected positive selection at NPC1-C in all mammalian orders, from Primates to Xenarthra. Different codons evolved adaptively in distinct mammals, and most selected sites are located within the two NPC1-C loops that engage GP, or at their anchor points. In Homininae, NPC1-C was a preferential selection target, and the T419I variant possibly represents a human-specific adaptation to filovirus infection. On the other side of the arms-race, GP evolved adaptively during filovirus speciation. One of the selected sites (S142Q) establishes several atom-to-atom contacts with NPC1-C. Additional selected sites are located within epitopes recognized by neutralizing antibodies, including the 14G7 epitope, where sites selected during the recent EBOV epidemic also map. Finally, pairs of co-evolving sites in Marburgviruses and Ebolaviruses were found to involve antigenic determinants. These findings suggest that the host humoral immune response was a major selective pressure during filovirus speciation. The S142Q variant may contribute to determine Ebolavirus host range in the wild. If this were the case, EBOV/BDBV (S142) and SUDV (Q142) may not share the same reservoir(s).

Competitive Fitness of Influenza B Viruses Possessing E119A and H274Y Neuraminidase Inhibitor Resistance-Associated Substitutions in Ferrets

Neuraminidase (NA) inhibitors (NAIs) are the only antiviral drugs recommended for influenza treatment and prophylaxis. Although NAI-resistant influenza B viruses that could pose a threat to public health have been reported in the field, their fitness is poorly understood. We evaluated in ferrets the pathogenicity and relative fitness of reverse genetics (rg)-generated influenza B/Yamanashi/166/1998-like viruses containing E119A or H274Y NA substitutions (N2 numbering). Ferrets inoculated with NAI-susceptible rg-wild-type (rg-WT) or NAI-resistant (rg-E119A or rg-H274Y) viruses developed mild infections. Growth of rg-E119A virus in the nasal cavities was delayed, but the high titers at 3 days post-inoculation (dpi) were comparable to those of the rg-WT and rg-H274Y viruses (3.6-4.1 log10TCID50/mL). No virus persisted beyond 5 dpi and replication did not extend to the trachea or lungs. Positive virus antigen-staining of the nasal turbinate epithelium was intermittent with the rg-WT and rg-H274Y viruses; whereas antigen-staining for the rg-E119A virus was more diffuse. Virus populations in ferrets coinoculated with NAI-susceptible and -resistant viruses (1:1 mixture) remained heterogeneous at 5 dpi but were predominantly rg-WT (>70%). Although the E119A substitution was associated with delayed replication in ferrets, the H274Y substitution did not measurably affect viral growth properties. These data suggest that rg-H274Y has undiminished fitness in single virus inoculations, but neither rg-E119A nor rg-H274Y gained a fitness advantage over rg-WT in direct competition experiments without antiviral drug pressure. Taken together, our data suggest the following order of relative fitness in a ferret animal model: rg-WT > rg-H274Y > rg-E119A.

Answering evolutionary questions: A guide for mechanistic biologists

The questions and methods of molecular biology and evolutionary biology are clearly distinct, yet a unified approach can lead to deep insights. Unfortunately, attempts to unify these approaches are fraught with pitfalls. In this informal series of questions and answers, we offer the mechanistically oriented biologist a set of steps to come up with evolutionarily reasonable and meaningful hypotheses. We emphasize the critical power and importance of carefully constructed null hypotheses, and we illustrate our ideas with examples representing a range of topics, from the biology of aging, to protein structure, to speciation, and more. We also stress the importance of mathematics as the lingua franca for biologists of all stripes, and encourage mechanistic biologists to seek out quantitative collaborators to build explicit mathematical models, making their assumptions explicit, and their logic clear and testable. Biologists in all realms of inquiry stand to gain from strong bridges between our disciplines.

Epistasis in protein evolution

The structure, function, and evolution of proteins depend on physical and genetic interactions among amino acids. Recent studies have used new strategies to explore the prevalence, biochemical mechanisms, and evolutionary implications of these interactions-called epistasis-within proteins. Here we describe an emerging picture of pervasive epistasis in which the physical and biological effects of mutations change over the course of evolution in a lineage-specific fashion. Epistasis can restrict the trajectories available to an evolving protein or open new paths to sequences and functions that would otherwise have been inaccessible. We describe two broad classes of epistatic interactions, which arise from different physical mechanisms and have different effects on evolutionary processes. Specific epistasis-in which one mutation influences the phenotypic effect of few other mutations-is caused by direct and indirect physical interactions between mutations, which nonadditively change the protein's physical properties, such as conformation, stability, or affinity for ligands. In contrast, nonspecific epistasis describes mutations that modify the effect of many others; these typically behave additively with respect to the physical properties of a protein but exhibit epistasis because of a nonlinear relationship between the physical properties and their biological effects, such as function or fitness. Both types of interaction are rampant, but specific epistasis has stronger effects on the rate and outcomes of evolution, because it imposes stricter constraints and modulates evolutionary potential more dramatically; it therefore makes evolution more contingent on low-probability historical events and leaves stronger marks on the sequences, structures, and functions of protein families.

Antiviral drug resistance as an adaptive process

Antiviral drug resistance is a matter of great clinical importance that, historically, has been investigated mostly from a virological perspective. Although the proximate mechanisms of resistance can be readily uncovered using these methods, larger evolutionary trends often remain elusive. Recent interest by population geneticists in studies of antiviral resistance has spurred new metrics for evaluating mutation and recombination rates, demographic histories of transmission and compartmentalization, and selective forces incurred during viral adaptation to antiviral drug treatment. We present up-to-date summaries on antiviral resistance for a range of drugs and viral types, and review recent advances for studying their evolutionary histories. We conclude that information imparted by demographic and selective histories, as revealed through population genomic inference, is integral to assessing the evolution of antiviral resistance as it pertains to human health.

High-Efficiency Capture of Drug Resistant-Influenza Virus by Live Imaging of Sialidase Activity

Influenza A and B viruses possess a neuraminidase protein that shows sialidase activity. Influenza virus-specific neuraminidase inhibitors (NAIs) are commonly used for clinical treatment of influenza. However, some influenza A and B viruses that are resistant to NAIs have emerged in nature. NAI-resistant viruses have been monitored in public hygiene surveys and the mechanism underlying the resistance has been studied. Here, we describe a new assay for selective detection and isolation of an NAI-resistant virus in a speedy and easy manner by live fluorescence imaging of viral sialidase activity, which we previously developed, in order to achieve high-efficiency capture of an NAI-resistant virus. An NAI-resistant virus maintains sialidase activity even at a concentration of NAI that leads to complete deactivation of the virus. Infected cells and focuses (infected cell populations) of an oseltamivir-resistant virus were selectively visualized by live fluorescence sialidase imaging in the presence of oseltamivir, resulting in high-efficiency isolation of the resistant viruses. The use of a combination of other NAIs (zanamivir, peramivir, and laninamivir) in the imaging showed that the oseltamivir-resistant virus isolated in 2008 was sensitive to zanamivir and laninamivir but resistant to peramivir. Fluorescence imaging in the presence of zanamivir also succeeded in selective live-cell visualization of cells that expressed zanamivir-resistant NA. Fluorescence imaging of NAI-resistant sialidase activity will be a powerful method for study of the NAI resistance mechanism, for public monitoring of NAI-resistant viruses, and for development of a new NAI that shows an effect on various NAI-resistant mutations.

Both Epistasis and Diversifying Selection Drive the Structural Evolution of the Ebola Virus Glycoprotein Mucin-Like Domain

Throughout the last 3 decades, Ebola virus (EBOV) outbreaks have been confined to isolated areas within Central Africa; however, the 2014 variant reached unprecedented transmission and mortality rates. While the outbreak was still under way, it was reported that the variant leading up to this outbreak evolved faster than previous EBOV variants, but evidence for diversifying selection was undetermined. Here, we test this selection hypothesis and show that while previous EBOV outbreaks were preceded by bursts of diversification, evidence for site-specific diversifying selection during the emergence of the 2014 EBOV clade is weak. However, we show strong evidence supporting an interplay between selection and correlated evolution (epistasis), particularly in the mucin-like domain (MLD) of the EBOV glycoprotein. By reconstructing ancestral structures of the MLD, we further propose a structural mechanism explaining how the substitutions that accumulated between 1918 and 1969 distorted the MLD, while more recent epistatic substitutions restored part of the structure, with the most recent substitution being adaptive. We suggest that it is this complex interplay between weak selection, epistasis, and structural constraints that has shaped the evolution of the 2014 EBOV variant.

IMPORTANCE

The role that selection plays in the emergence of viral epidemics remains debated, particularly in the context of the 2014 EBOV outbreak. Most critically, should such evidence exist, it is generally unclear how this relates to function and increased virulence. Here, we show that the viral lineage leading up to the 2014 outbreak underwent a complex interplay between selection and correlated evolution (epistasis) in a protein region that is critical for immune evasion. We then reconstructed the three-dimensional structure of this domain and showed that the initial mutations in this lineage deformed the structure, while subsequent mutations restored part of the structure. Along this mutational path, the first and last mutations were adaptive, while the intervening ones were epistatic. Altogether, we provide a mechanistic model that explains how selection and epistasis acted on the structural constraints that materialized during the 2014 EBOV outbreak.

Development and applications of single-cycle infectious influenza A virus (sciIAV)

The diverse host range, high transmissibility, and rapid evolution of influenza A viruses justify the importance of containing pathogenic viruses studied in the laboratory. Other than physically or mechanically changing influenza A virus containment procedures, modifying the virus to only replicate for a single round of infection similarly ensures safety and consequently decreases the level of biosafety containment required to study highly pathogenic members in the virus family. This biological containment is more ideal because it is less apt to computer, machine, or human error. With many necessary proteins that can be deleted, generation of single-cycle infectious influenza A viruses (sciIAV) can be achieved using a variety of approaches. Here, we review the recent burst in sciIAV generation and summarize the applications and findings on this important human pathogen using biocontained viral mimics.

Functional Sites Induce Long-Range Evolutionary Constraints in Enzymes

Functional residues in proteins tend to be highly conserved over evolutionary time. However, to what extent functional sites impose evolutionary constraints on nearby or even more distant residues is not known. Here, we report pervasive conservation gradients toward catalytic residues in a dataset of 524 distinct enzymes: evolutionary conservation decreases approximately linearly with increasing distance to the nearest catalytic residue in the protein structure. This trend encompasses, on average, 80% of the residues in any enzyme, and it is independent of known structural constraints on protein evolution such as residue packing or solvent accessibility. Further, the trend exists in both monomeric and multimeric enzymes and irrespective of enzyme size and/or location of the active site in the enzyme structure. By contrast, sites in protein-protein interfaces, unlike catalytic residues, are only weakly conserved and induce only minor rate gradients. In aggregate, these observations show that functional sites, and in particular catalytic residues, induce long-range evolutionary constraints in enzymes.

Zanamivir-resistant influenza viruses with Q136K or Q136R neuraminidase residue mutations can arise during MDCK cell culture creating challenges for antiviral susceptibility monitoring

Surveillance of circulating influenza strains for antiviral susceptibility is important to ensure patient treatment guidelines remain appropriate. Influenza A(H3N2) and A(H1N1)pdm09 virus isolates containing mutations at the Q136 residue of the neuraminidase (NA) that conferred reduced susceptibility to the NA inhibitor (NAI) zanamivir were detected during antiviral susceptibility monitoring. Interestingly, the mutations were not detectable in the viruses from respective clinical specimens, only in the cultured isolates. We showed that variant viruses containing the Q136K and Q136R NA mutations were preferentially selected in Madin-Darby canine kidney epithelial (MDCK) cells, but were less well supported in MDCK-SIAT1 cells and embryonated eggs. The effect of Q136K, Q136R, Q136H and Q136L substitutions in NA subtypes N1 and N2 on NAI susceptibility and in vitro viral fitness was assessed. This study highlights the challenges that cell culture derived mutations can pose to the NAI susceptibility analysis and interpretation and reaffirms the need to sequence viruses from respective clinical specimens to avoid misdiagnosis. However, we also demonstrate that NA mutations at residue Q136 can confer reduced zanamivir, peramivir or laninamivir susceptibility, and therefore close monitoring of viruses for mutations at this site from patients being treated with these antivirals is important.

CD4 Receptor is a Key Determinant of Divergent HIV-1 Sensing by Plasmacytoid Dendritic Cells

Plasmacytoid dendritic cells (pDC) are innate immune cells that sense viral nucleic acids through endosomal Toll-like receptor (TLR) 7/9 to produce type I interferon (IFN) and to differentiate into potent antigen presenting cells (APC). Engagement of TLR7/9 in early endosomes appears to trigger the IRF7 pathway for IFN production whereas engagement in lysosomes seems to trigger the NF-κB pathway for maturation into APC. We showed previously that HIV-1 (HIV) localizes predominantly to early endosomes, not lysosomes, and mainly stimulate IRF7 rather than NF-κB signaling pathways in pDC. This divergent signaling may contribute to disease progression through production of pro-apoptotic and pro-inflammatory IFN and inadequate maturation of pDCs. We now demonstrate that HIV virions may be re-directed to lysosomes for NF-κB signaling by either pseudotyping HIV with influenza hemagglutinin envelope or modification of CD4 mediated-intracellular trafficking. These data suggest that HIV envelope-CD4 receptor interactions drive pDC activation toward an immature IFN producing phenotype rather than differentiation into a mature dendritic cell phenotype.

Plasmacytoid dendritic cells (pDC) are innate immune cells that are specialized to produce type I interferon (IFN) and to activate adaptive immune responses. Although IFN is an anti-viral cytokine, it may contribute more to pathogenesis than to protection during chronic viral infections, including chronic HIV infection. pDC sense HIV to produce abundant IFN but minimal NF- κB–dependent production of TNFα and minimal up-regulation of co-stimulatory molecules, suggesting that HIV promotes pDC to become interferon producing cells (IPC) rather than antigen presenting cells (APC). Here, we use florescent HIV virions pseudotyped with influenza hemagglutinin (HA) envelope and a cell system expressing CD4 molecules with modified intracellular trafficking. We found that HIV virions pseudotyped with HA stimulate pDC to mature, similar to influenza-stimulated pDC, and traffic intracellularly similarly to influenza. We also find that CD4-mediated intracellular trafficking guides HIV trafficking and downstream signaling. Our study presents new and important findings which demonstrate that divergent HIV sensing by pDC to produce IFN, rather than to become mature antigen presenting cells, is mediated specifically by CD4-HIV envelope interactions.

Neuraminidase Activity and Resistance of 2009 Pandemic H1N1 Influenza Virus to Antiviral Activity in Bronchoalveolar Fluid

Human bronchoalveolar fluid is known to have anti-influenza activity. It is believed to be a frontline innate defense against the virus. Several antiviral factors, including surfactant protein D, are believed to contribute to the activity. The 2009 pandemic H1N1 influenza virus was previously shown to be less sensitive to surfactant protein D. Nevertheless, whether different influenza virus strains have different sensitivities to the overall anti-influenza activity of human bronchoalveolar fluid was not known. We compared the sensitivities of 2009 pandemic H1N1, seasonal H1N1, and seasonal H3N2 influenza virus strains to inhibition by human bronchoalveolar lavage (BAL) fluid. The pandemic and seasonal H1N1 strains showed lower sensitivity to human BAL fluid than the H3N2 strains. The BAL fluid anti-influenza activity could be enhanced by oseltamivir, indicating that the viral neuraminidase (NA) activity could provide resistance to the antiviral defense. In accordance with this finding, the BAL fluid anti-influenza activity was found to be sensitive to sialidase. The oseltamivir resistance mutation H275Y rendered the pandemic H1N1 virus but not the seasonal H1N1 virus more sensitive to BAL fluid. Since only the seasonal H1N1 but not the pandemic H1N1 had compensatory mutations that allowed oseltamivir-resistant strains to maintain NA enzymatic activity and transmission fitness, the resistance to BAL fluid of the drug-resistant seasonal H1N1 virus might play a role in viral fitness.

IMPORTANCE

Human airway secretion contains anti-influenza activity. Different influenza strains may vary in their susceptibilities to this antiviral activity. Here we show that the 2009 pandemic and seasonal H1N1 influenza viruses were less sensitive to human bronchoalveolar lavage (BAL) fluid than H3N2 seasonal influenza virus. The resistance to the pulmonary innate antiviral activity of the pandemic virus was determined by its neuraminidase (NA) gene, and it was shown that the NA inhibitor resistance mutation H275Y abolished this resistance of the pandemic H1N1 but not the seasonal H1N1 virus, which had compensatory mutations that maintained the fitness of drug-resistant strains. Therefore, the innate respiratory tract defense may be a barrier against NA inhibitor-resistant mutants, and evasion of this defense may play a role in the emergence and spread of drug-resistant strains.

Targeted disruption of influenza A virus hemagglutinin in genetically modified mice reduces viral replication and improves disease outcome

Influenza A virus can cause acute respiratory infection in animals and humans around the globe, and is still a major threat to animal husbandry and public health. Due to antigenic drift and antigenic shift of the virus, development of novel anti-influenza strategies has become an urgent task. Here we generated transgenic (TG) mice stably expressing a short-hairpin RNA specifically targeting hemagglutinin (HA) of influenza A virus, and investigated the susceptibility of the mice to influenza virus infection. We found that HA expression was dramatically disrupted in TG mice infected with WSN or PR8 virus. Importantly, the animals showed reduced virus production in lungs, slower weight loss, attenuated acute organ injury and consequently increased survival rates as compared to wild type (WT) mice after the viral infection. Moreover, TG mice exhibited a normal level of white blood cells following the virus infection, whereas the number of these cells was significantly decreased in WT mice with same challenge. Together, these experiments demonstrate that the TG mice are less permissive for influenza virus replication, and suggest that shRNA-based efficient disruption of viral gene expression in animals may be a useful strategy for prevention and control of a viral zoonosis.

Causes of molecular convergence and parallelism in protein evolution

To what extent is the convergent evolution of protein function attributable to convergent or parallel changes at the amino acid level? The mutations that contribute to adaptive protein evolution may represent a biased subset of all possible beneficial mutations owing to mutation bias and/or variation in the magnitude of deleterious pleiotropy. A key finding is that the fitness effects of amino acid mutations are often conditional on genetic background. This context dependence (epistasis) can reduce the probability of convergence and parallelism because it reduces the number of possible mutations that are unconditionally acceptable in divergent genetic backgrounds. Here, I review factors that influence the probability of replicated evolution at the molecular level.

Human Sentinel Surveillance of Influenza and Other Respiratory Viral Pathogens in Border Areas of Western Cambodia

Little is known about circulation of influenza and other respiratory viruses in remote populations along the Thai-Cambodia border in western Cambodia. We screened 586 outpatients (median age 5, range 1–77) presenting with influenza-like-illness (ILI) at 4 sentinel sites in western Cambodia between May 2010 and December 2012. Real-time reverse transcriptase (rRT) PCR for influenza was performed on combined nasal and throat specimens followed by viral culture, antigenic analysis, antiviral susceptibility testing and full genome sequencing for phylogenetic analysis. ILI-specimens negative for influenza were cultured, followed by rRT-PCR for enterovirus and rhinovirus (EV/RV) and EV71. Influenza was found in 168 cases (29%) and occurred almost exclusively in the rainy season from June to November. Isolated influenza strains had close antigenic and phylogenetic relationships, matching vaccine and circulating strains found elsewhere in Cambodia. Influenza vaccination coverage was low (<20%). Western Cambodian H1N1(2009) isolate genomes were more closely related to 10 earlier Cambodia isolates (94.4% genome conservation) than to 13 Thai isolates (75.9% genome conservation), despite sharing the majority of the amino acid changes with the Thai references. Most genes showed signatures of purifying selection. Viral culture detected only adenovirus (5.7%) and parainfluenza virus (3.8%), while non-polio enteroviruses (10.3%) were detected among 164 culture-negative samples including coxsackievirus A4, A6, A8, A9, A12, B3, B4 and echovirus E6 and E9 using nested RT-PCR methods. A single specimen of EV71 was found. Despite proximity to Thailand, influenza epidemiology of these western Cambodian isolates followed patterns observed elsewhere in Cambodia, continuing to support current vaccine and treatment recommendations from the Cambodian National Influenza Center. Amino acid mutations at non-epitope sites, particularly hemagglutinin genes, require further investigation in light of an increasingly important role of permissive mutations in influenza virus evolution. Further research about the burden of adenovirus and non-polio enteroviruses as etiologic agents in acute respiratory infections in Cambodia is also needed.

Cooperation between distinct viral variants promotes growth of H3N2 influenza in cell culture

RNA viruses rapidly diversify into quasispecies of related genotypes. This genetic diversity has long been known to facilitate adaptation, but recent studies have suggested that cooperation between variants might also increase population fitness. Here, we demonstrate strong cooperation between two H3N2 influenza variants that differ by a single mutation at residue 151 in neuraminidase, which normally mediates viral exit from host cells. Residue 151 is often annotated as an ambiguous amino acid in sequenced isolates, indicating mixed viral populations. We show that mixed populations grow better than either variant alone in cell culture. Pure populations of either variant generate the other through mutation and then stably maintain a mix of the two genotypes. We suggest that cooperation arises because mixed populations combine one variant's proficiency at cell entry with the other's proficiency at cell exit. Our work demonstrates a specific cooperative interaction between defined variants in a viral quasispecies.

Engineering and Evolution of Molecular Chaperones and Protein Disaggregases with Enhanced Activity

Cells have evolved a sophisticated proteostasis network to ensure that proteins acquire and retain their native structure and function. Critical components of this network include molecular chaperones and protein disaggregases, which function to prevent and reverse deleterious protein misfolding. Nevertheless, proteostasis networks have limits, which when exceeded can have fatal consequences as in various neurodegenerative disorders, including Parkinson's disease and amyotrophic lateral sclerosis. A promising strategy is to engineer proteostasis networks to counter challenges presented by specific diseases or specific proteins. Here, we review efforts to enhance the activity of individual molecular chaperones or protein disaggregases via engineering and directed evolution. Remarkably, enhanced global activity or altered substrate specificity of various molecular chaperones, including GroEL, Hsp70, ClpX, and Spy, can be achieved by minor changes in primary sequence and often a single missense mutation. Likewise, small changes in the primary sequence of Hsp104 yield potentiated protein disaggregases that reverse the aggregation and buffer toxicity of various neurodegenerative disease proteins, including α-synuclein, TDP-43, and FUS. Collectively, these advances have revealed key mechanistic and functional insights into chaperone and disaggregase biology. They also suggest that enhanced chaperones and disaggregases could have important applications in treating human disease as well as in the purification of valuable proteins in the pharmaceutical sector.

The significance of naturally occurring neuraminidase quasispecies of H5N1 avian influenza virus on resistance to oseltamivir: a point of concern

Viral adaptability and survival arise due to the presence of quasispecies populations that are able to escape the immune response or produce drug-resistant variants. However, the presence of H5N1 virus with natural mutations acquired without any drug selection pressure poses a great threat. Cloacal samples collected from the 2004-2005 epidemics in Thailand from Asian open-billed storks revealed one major and several minor quasispecies populations with mutations on the oseltamivir (OTV)-binding site of the neuraminidase gene (NA) without prior exposure to a drug. Therefore, this study investigated the binding between the NA-containing novel mutations and OTV drug using molecular dynamic simulations and plaque inhibition assay. The results revealed that the mutant populations, S236F mutant, S236F/C278Y mutant, A250V/V266A/P271H/G285S mutant and C278Y mutant, had a lower binding affinity with OTV as compared with the WT virus due to rearrangement of amino acid residues and increased flexibility in the 150-loop. This result was further emphasized through the IC50 values obtained for the major population and WT virus, 104.74 nM and 18.30 nM, respectively. Taken together, these data suggest that H5N1 viruses isolated from wild birds have already acquired OTV-resistant point mutations without any exposure to a drug.

Application of virus-like particles (VLP) to NMR characterization of viral membrane protein interactions

The membrane proteins of viruses play critical roles in the virus life cycle and are attractive targets for therapeutic intervention. Virus-like particles (VLP) present the possibility to study the biochemical and biophysical properties of viral membrane proteins in their native environment. Specifically, the VLP constructs contain the entire protein sequence and are comprised of native membrane components including lipids, cholesterol, carbohydrates and cellular proteins. In this study we prepare VLP containing full-length hemagglutinin (HA) or neuraminidase (NA) from influenza and characterize their interactions with small molecule inhibitors. Using HA-VLP, we first show that VLP samples prepared using the standard sucrose gradient purification scheme contain significant amounts of serum proteins, which exhibit high potential for non-specific interactions, thereby complicating NMR studies of ligand-target interactions. We then show that the serum contaminants may be largely removed with the addition of a gel filtration chromatography step. Next, using HA-VLP we demonstrate that WaterLOGSY NMR is significantly more sensitive than Saturation Transfer Difference (STD) NMR for the study of ligand interactions with membrane bound targets. In addition, we compare the ligand orientation to HA embedded in VLP with that of recombinant HA by STD NMR. In a subsequent step, using NA-VLP we characterize the kinetic and binding properties of substrate analogs and inhibitors of NA, including study of the H274Y-NA mutant, which leads to wide spread resistance to current influenza antivirals. In summary, our work suggests that VLP have high potential to become standard tools in biochemical and biophysical studies of viral membrane proteins, particularly when VLP are highly purified and combined with control VLP containing native membrane proteins.

Coupling high-throughput genetics with phylogenetic information reveals an epistatic interaction on the influenza A virus M segment

Background

Epistasis is one of the central themes in viral evolution due to its importance in drug resistance, immune escape, and interspecies transmission. However, there is a lack of experimental approach to systematically probe for epistatic residues.

Results

By utilizing the information from natural occurring sequences and high-throughput genetics, this study established a novel strategy to identify epistatic residues. The rationale is that a substitution that is deleterious in one strain may be prevalent in nature due to the presence of a naturally occurring compensatory substitution. Here, high-throughput genetics was applied to influenza A virus M segment to systematically identify deleterious substitutions. Comparison with natural sequence variation showed that a deleterious substitution M1 Q214H was prevalent in circulating strains. A coevolution analysis was then performed and indicated that M1 residues 121, 207, 209, and 214 naturally coevolved as a group. Subsequently, we experimentally validated that M1 A209T was a compensatory substitution for M1 Q214H.

Conclusions

This work provided a proof-of-concept to identify epistatic residues by coupling high-throughput genetics with phylogenetic information. In particular, we were able to identify an epistatic interaction between M1 substitutions A209T and Q214H. This analytic strategy can potentially be adapted to study any protein of interest, provided that the information on natural sequence variants is available.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2358-7) contains supplementary material, which is available to authorized users.

The population genetics of drug resistance evolution in natural populations of viral, bacterial and eukaryotic pathogens

Drug resistance is a costly consequence of pathogen evolution and a major concern in public health. In this review, we show how population genetics can be used to study the evolution of drug resistance and also how drug resistance evolution is informative as an evolutionary model system. We highlight five examples from diverse organisms with particular focus on: (i) identifying drug resistance loci in the malaria parasite Plasmodium falciparum using the genomic signatures of selective sweeps, (ii) determining the role of epistasis in drug resistance evolution in influenza, (iii) quantifying the role of standing genetic variation in the evolution of drug resistance in HIV, (iv) using drug resistance mutations to study clonal interference dynamics in tuberculosis and (v) analysing the population structure of the core and accessory genome of Staphylococcus aureus to understand the spread of methicillin resistance. Throughout this review, we discuss the uses of sequence data and population genetic theory in studying the evolution of drug resistance.

Distal substitutions drive divergent DNA specificity among paralogous transcription factors through subdivision of conformational space

Many genomes contain families of paralogs--proteins with divergent function that evolved from a common ancestral gene after a duplication event. To understand how paralogous transcription factors evolve divergent DNA specificities, we examined how the glucocorticoid receptor and its paralogs evolved to bind activating response elements [(+)GREs] and negative glucocorticoid response elements (nGREs). We show that binding to nGREs is a property of the glucocorticoid receptor (GR) DNA-binding domain (DBD) not shared by other members of the steroid receptor family. Using phylogenetic, structural, biochemical, and molecular dynamics techniques, we show that the ancestral DBD from which GR and its paralogs evolved was capable of binding both nGRE and (+)GRE sequences because of the ancestral DBD's ability to assume multiple DNA-bound conformations. Subsequent amino acid substitutions in duplicated daughter genes selectively restricted protein conformational space, causing this dual DNA-binding specificity to be selectively enhanced in the GR lineage and lost in all others. Key substitutions that determined the receptors' response element-binding specificity were far from the proteins' DNA-binding interface and interacted epistatically to change the DBD's function through DNA-induced allosteric mechanisms. These amino acid substitutions subdivided both the conformational and functional space of the ancestral DBD among the present-day receptors, allowing a paralogous family of transcription factors to control disparate transcriptional programs despite high sequence identity.

Residue-based design of small molecule inhibitor for H1N1, H5N1 and H7N1 mutants

Point mutations H274Y and N294S can lead to resistance of influenza virus strains to some drug molecules. Recently, a large number of experiments has focused on the many frameworks and catalytic residues thought to prevent the efficacy of anti-flu drugs. In the past, most research has considered the role of drugs in rigid proteins rather than in flexible proteins. In this study, we used molecular dynamics simulation (MD) combined with structure- and ligand-based drug design (SBDD and LBDD) methods to study dynamic interaction and protein dynamics correlation statistics between compounds and both the framework and catalytic residues in influenza virus N1 strains. Drug candidates were screened using the IC50 of the docking result predicted by support vector machine, multiple linear regression, and genetic function approximation (P < 0.001). As shown by MD, saussureamine C and diiodotyrosine have a protein dynamics correlation similar to that of sialic acid, and both can participate in hydrogen bond formation with loop, framework, and catalytic residues. Our in silico findings suggest that saussureamine C can inhibit H274Y and N294S mutants, and that diiodotyrosine can also inhibit N294S mutants. Therefore, the drugs saussureamine C and diiodotyrosine have the potential to produce inhibitory effects on wild-type influenza virus and some N1 mutants.

A Balance between Inhibitor Binding and Substrate Processing Confers Influenza Drug Resistance

The therapeutic benefits of the neuraminidase (NA) inhibitor oseltamivir are dampened by the emergence of drug resistance mutations in influenza A virus (IAV). To investigate the mechanistic features that underlie resistance, we developed an approach to quantify the effects of all possible single-nucleotide substitutions introduced into important regions of NA. We determined the experimental fitness effects of 450 nucleotide mutations encoding positions both surrounding the active site and at more distant sites in an N1 strain of IAV in the presence and absence of oseltamivir. NA mutations previously known to confer oseltamivir resistance in N1 strains, including H275Y and N295S, were adaptive in the presence of drug, indicating that our experimental system captured salient features of real-world selection pressures acting on NA. We identified mutations, including several at position 223, that reduce the apparent affinity for oseltamivir in vitro. Position 223 of NA is located adjacent to a hydrophobic portion of oseltamivir that is chemically distinct from the substrate, making it a hotspot for substitutions that preferentially impact drug binding relative to substrate processing. Furthermore, two NA mutations, K221N and Y276F, each reduce susceptibility to oseltamivir by increasing NA activity without altering drug binding. These results indicate that competitive expansion of IAV in the face of drug pressure is mediated by a balance between inhibitor binding and substrate processing.

The effects of a deleterious mutation load on patterns of influenza A/H3N2's antigenic evolution in humans

Recent phylogenetic analyses indicate that RNA virus populations carry a significant deleterious mutation load. This mutation load has the potential to shape patterns of adaptive evolution via genetic linkage to beneficial mutations. Here, we examine the effect of deleterious mutations on patterns of influenza A subtype H3N2's antigenic evolution in humans. By first analyzing simple models of influenza that incorporate a mutation load, we show that deleterious mutations, as expected, act to slow the virus's rate of antigenic evolution, while making it more punctuated in nature. These models further predict three distinct molecular pathways by which antigenic cluster transitions occur, and we find phylogenetic patterns consistent with each of these pathways in influenza virus sequences. Simulations of a more complex phylodynamic model further indicate that antigenic mutations act in concert with deleterious mutations to reproduce influenza's spindly hemagglutinin phylogeny, co-circulation of antigenic variants, and high annual attack rates.

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

Each year, up to 15% of the world's population experience symptoms of an influenza infection, also commonly known as flu. The most common culprit is a strain of the virus called influenza type A subtype H3N2. One reason that so many people become infected each year is that this virus evolves rapidly. Within a few years, proteins on the surface of the virus known as antigens become less recognizable to the immune system of a person who has been previously infected. This means that the person can become ill with the virus again because their immune system cannot mount an effective response to the evolved virus strain.

Influenza virus strains evolve rapidly because their genetic material accumulates mutations quickly. Although some of these mutations are beneficial to the virus, other mutations are harmful and reduce the ability of the virus to spread. Sometimes beneficial mutations may occur alongside harmful ones, but it is not known how the harmful mutations affect the evolution of the virus.

Here, Koelle and Rasmussen used computer models of H3N2 influenza to examine the effect of harmful mutations on the evolution of this virus population. The models show that harmful mutations limit how quickly the antigens can evolve. Also, the presence of these harmful mutations effectively acts as a sieve: they allow only large changes in the antigens to establish in the virus population.

The models suggest that there are three routes by which large changes in the antigens on H3N2 viruses may occur. The first is by a single mutation that has a big effect on the antigens in viruses that only carry a few harmful mutations, but these large mutations would not happen very often. Another route may be through more common mutations that have only a small or moderate benefit, which would allow the virus to become more common in the population before it acquires a beneficial mutation with a much greater effect. The third possibility is that a large beneficial mutation may arise in viruses that have many harmful mutations. These harmful mutations may initially limit the ability of the virus to spread, but over time, some of these harmful mutations may then be lost.

Koelle and Rasmussen found that the computer models could recreate the patterns of virus evolution that have been observed in real strains of H3N2. Researchers use predictions of influenza evolution to help them decide which virus strains should be included in flu vaccines each year. Koelle and Rasmussen findings indicate that harmful mutations should be considered when making these predictions.

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

Oseltamivir-resistant influenza A(H1N1)pdm09 viruses, United States, 2013-14

We report characteristics of oseltamivir-resistant influenza A(H1N1)pdm09 viruses and patients infected with these viruses in the United States. During 2013-14, fifty-nine (1.2%) of 4,968 analyzed US influenza A(H1N1)pdm09 viruses had the H275Y oseltamivir resistance-conferring neuraminidase substitution. Our results emphasize the need for local surveillance for neuraminidase inhibitor susceptibility among circulating influenza viruses.

New small-molecule drug design strategies for fighting resistant influenza A

Influenza A virus is the major cause of seasonal or pandemic flu worldwide. Two main treatment strategies-vaccination and small molecule anti-influenza drugs are currently available. As an effective vaccine usually takes at least 6 months to develop, anti-influenza small molecule drugs are more effective for the first line of protection against the virus during an epidemic outbreak, especially in the early stage. Two major classes of anti-influenza drugs currently available are admantane-based M2 protein blockers (amantadine and rimantadine) and neuraminidase (NA) inhibitors (oseltamivir, zanamivir, and peramivir). However, the continuous evolvement of influenza A virus and the rapid emergence of resistance to current drugs, particularly to amantadine, rimantadine, and oseltamivir, have raised an urgent need for developing new anti-influenza drugs against resistant forms of influenza A virus. In this review, we first give a brief introduction of the molecular mechanisms behind resistance, and then discuss new strategies in small-molecule drug development to overcome influenza A virus resistance targeting mutant M2 proteins and neuraminidases, and other viral proteins not associated with current drugs.

Detecting epistasis from an ensemble of adapting populations

The role that epistasis plays during adaptation remains an outstanding problem, which has received considerable attention in recent years. Most of the recent empirical studies are based on ensembles of replicate populations that adapt in a fixed, laboratory controlled condition. Researchers often seek to infer the presence and form of epistasis in the fitness landscape from the time evolution of various statistics averaged across the ensemble of populations. Here, we provide a rigorous analysis of what quantities, drawn from time series of such ensembles, can be used to infer epistasis for populations evolving under weak mutation on finite-site fitness landscapes. First, we analyze the mean fitness trajectory-that is, the time course of the ensemble average fitness. We show that for any epistatic fitness landscape and starting genotype, there always exists a non-epistatic fitness landscape that produces the exact same mean fitness trajectory. Thus, the presence of epistasis is not identifiable from the mean fitness trajectory. By contrast, we show that two other ensemble statistics-the time evolution of the fitness variance across populations, and the time evolution of the mean number of substitutions-can detect certain forms of epistasis in the underlying fitness landscape.

The virtue of innovation: innovation through the lenses of biological evolution

We rehearse the processes of innovation and discovery in general terms, using as our main metaphor the biological concept of an evolutionary fitness landscape. Incremental and disruptive innovations are seen, respectively, as successful searches carried out locally or more widely. They may also be understood as reflecting evolution by mutation (incremental) versus recombination (disruptive). We also bring a platonic view, focusing on virtue and memory. We use 'virtue' as a measure of efforts, including the knowledge required to come up with disruptive and incremental innovations, and 'memory' as a measure of their lifespan, i.e. how long they are remembered. Fostering innovation, in the evolutionary metaphor, means providing the wherewithal to promote novelty, good objective functions that one is trying to optimize, and means to improve one's knowledge of, and ability to navigate, the landscape one is searching. Recombination necessarily implies multi- or inter-disciplinarity. These principles are generic to all kinds of creativity, novel ideas formation and the development of new products and technologies.

Multiplexed detection of viral infections using rapid in situ RNA analysis on a chip

Viral infections are a major cause of human disease, but many require molecular assays for conclusive diagnosis. Current assays typically rely on RT-PCR or ELISA; however, these tests often have limited speed, sensitivity or specificity. Here, we demonstrate that rapid RNA FISH is a viable alternative method that could improve upon these limitations. We describe a platform beginning with software to generate RNA FISH probes both for distinguishing related strains of virus (even those different by a single base) and for capturing large numbers of strains simultaneously. Next, we present a simple fluidic device for reliably performing RNA FISH assays in an automated fashion. Finally, we describe an automated image processing pipeline to robustly identify uninfected and infected samples. Together, our results establish RNA FISH as a methodology with potential for viral point-of-care diagnostics.

Coordinated Evolution of Influenza A Surface Proteins

The surface proteins hemagglutinin (HA) and neuraminidase (NA) of human influenza A virus evolve under selection pressures to escape adaptive immune responses and antiviral drug treatments. In addition to these external selection pressures, some mutations in HA are known to affect the adaptive landscape of NA, and vice versa, because these two proteins are physiologically interlinked. However, the extent to which evolution of one protein affects the evolution of the other one is unknown. Here we develop a novel phylogenetic method for detecting the signatures of such genetic interactions between mutations in different genes – that is, inter-gene epistasis. Using this method, we show that influenza surface proteins evolve in a coordinated way, with mutations in HA affecting subsequent spread of mutations in NA and vice versa, at many sites. Of particular interest is our finding that the oseltamivir-resistance mutations in NA in subtype H1N1 were likely facilitated by prior mutations in HA. Our results illustrate that the adaptive landscape of a viral protein is remarkably sensitive to its genomic context and, more generally, that the evolution of any single protein must be understood within the context of the entire evolving genome.

The fitness of an organism depends on the coordinated function of many genes. Thus, how a mutation in one gene affects fitness often depends on what mutations are present in other genes. This dependence is called “genetic interaction” or “epistasis”. The prevalence and type of such interactions are not well understood. Epistasis can be inferred from time-series sequencing data when a mutation in one gene is observed to facilitate the spread of a mutation in another gene. However, the situation is much more complicated when new combinations of genes are formed by processes such as recombination or reassortment. In such cases, deducing the time and order of genetic changes is difficult. Here, we devise a method to infer pairs of mutations in different genes which closely follow one another in the presence of reassortment. We apply it to evolution of two surface proteins of influenza A virus, hemagglutinin and neuraminidase, which are important targets for the human immune system and drugs. We show that mutations in one of these proteins are often facilitated by prior mutations, or compensated by subsequent mutations, in the other protein. In particular, drug-resistance mutations in neuraminidase were likely made possible by prior mutation in hemagglutinin. Knowledge of such interactions is necessary to fully understand and predict evolution.

Sequence-Specific Modifications Enhance the Broad-Spectrum Antiviral Response Activated by RIG-I Agonists

The cytosolic RIG-I (retinoic acid-inducible gene I) receptor plays a pivotal role in the initiation of the immune response against RNA virus infection by recognizing short 5'-triphosphate (5'ppp)-containing viral RNA and activating the host antiviral innate response. In the present study, we generated novel 5'ppp RIG-I agonists of varieous lengths, structures, and sequences and evaluated the generation of the antiviral and inflammatory responses in human epithelial A549 cells, human innate immune primary cells, and murine models of influenza and chikungunya viral pathogenesis. A 99-nucleotide, uridine-rich hairpin 5'pppRNA termed M8 stimulated an extensive and robust interferon response compared to other modified 5'pppRNA structures, RIG-I aptamers, or poly(I·C). Interestingly, manipulation of the primary RNA sequence alone was sufficient to modulate antiviral activity and inflammatory response, in a manner dependent exclusively on RIG-I and independent of MDA5 and TLR3. Both prophylactic and therapeutic administration of M8 effectively inhibited influenza virus and dengue virus replication in vitro. Furthermore, multiple strains of influenza virus that were resistant to oseltamivir, an FDA-approved therapeutic treatment for influenza, were highly sensitive to inhibition by M8. Finally, prophylactic M8 treatment in vivo prolonged survival and reduced lung viral titers of mice challenged with influenza virus, as well as reducing chikungunya virus-associated foot swelling and viral load. Altogether, these results demonstrate that 5'pppRNA can be rationally designed to achieve a maximal RIG-I-mediated protective antiviral response against human-pathogenic RNA viruses.

IMPORTANCE

The development of novel therapeutics to treat human-pathogenic RNA viral infections is an important goal to reduce spread of infection and to improve human health and safety. This study investigated the design of an RNA agonist with enhanced antiviral and inflammatory properties against influenza, dengue, and chikungunya viruses. A novel, sequence-dependent, uridine-rich RIG-I agonist generated a protective antiviral response in vitro and in vivo and was effective at concentrations 100-fold lower than prototype sequences or other RNA agonists, highlighting the robust activity and potential clinical use of the 5'pppRNA against RNA virus infection. Altogether, the results identify a novel, sequence-specific RIG-I agonist as an attractive therapeutic candidate for the treatment of a broad range of RNA viruses, a pressing issue in which a need for new and more effective options persists.

A genomic and evolutionary approach reveals non-genetic drug resistance in malaria

BACKGROUND

Drug resistance remains a major public health challenge for malaria treatment and eradication. Individual loci associated with drug resistance to many antimalarials have been identified, but their epistasis with other resistance mechanisms has not yet been elucidated.

RESULTS

We previously described two mutations in the cytoplasmic prolyl-tRNA synthetase (cPRS) gene that confer resistance to halofuginone. We describe here the evolutionary trajectory of halofuginone resistance of two independent drug resistance selections in Plasmodium falciparum. Using this novel methodology, we discover an unexpected non-genetic drug resistance mechanism that P. falciparum utilizes before genetic modification of the cPRS. P. falciparum first upregulates its proline amino acid homeostasis in response to halofuginone pressure. We show that this non-genetic adaptation to halofuginone is not likely mediated by differential RNA expression and precedes mutation or amplification of the cPRS gene. By tracking the evolution of the two drug resistance selections with whole genome sequencing, we further demonstrate that the cPRS locus accounts for the majority of genetic adaptation to halofuginone in P. falciparum. We further validate that copy-number variations at the cPRS locus also contribute to halofuginone resistance.

CONCLUSIONS

We provide a three-step model for multi-locus evolution of halofuginone drug resistance in P. falciparum. Informed by genomic approaches, our results provide the first comprehensive view of the evolutionary trajectory malaria parasites take to achieve drug resistance. Our understanding of the multiple genetic and non-genetic mechanisms of drug resistance informs how we will design and pair future anti-malarials for clinical use.

A genomic and evolutionary approach reveals non-genetic drug resistance in malaria

Background

Drug resistance remains a major public health challenge for malaria treatment and eradication. Individual loci associated with drug resistance to many antimalarials have been identified, but their epistasis with other resistance mechanisms has not yet been elucidated.

Results

We previously described two mutations in the cytoplasmic prolyl-tRNA synthetase (cPRS) gene that confer resistance to halofuginone. We describe here the evolutionary trajectory of halofuginone resistance of two independent drug resistance selections in Plasmodium falciparum. Using this novel methodology, we discover an unexpected non-genetic drug resistance mechanism that P. falciparum utilizes before genetic modification of the cPRS. P. falciparum first upregulates its proline amino acid homeostasis in response to halofuginone pressure. We show that this non-genetic adaptation to halofuginone is not likely mediated by differential RNA expression and precedes mutation or amplification of the cPRS gene. By tracking the evolution of the two drug resistance selections with whole genome sequencing, we further demonstrate that the cPRS locus accounts for the majority of genetic adaptation to halofuginone in P. falciparum. We further validate that copy-number variations at the cPRS locus also contribute to halofuginone resistance.

Conclusions

We provide a three-step model for multi-locus evolution of halofuginone drug resistance in P. falciparum. Informed by genomic approaches, our results provide the first comprehensive view of the evolutionary trajectory malaria parasites take to achieve drug resistance. Our understanding of the multiple genetic and non-genetic mechanisms of drug resistance informs how we will design and pair future anti-malarials for clinical use.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0511-2) contains supplementary material, which is available to authorized users.