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

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 utilize an experimental evolution framework, combined with novel population genetic method development for inference from time-sampled data, in order to evaluate the effectiveness of favipiravir against IAV. Evaluating whole genome polymorphism data across fifteen time points under multiple drug concentrations and in controls, we present the first evidence for the ability of viral 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

Positive Selection Drives Preferred Segment Combinations during Influenza Virus Reassortment

Influenza A virus (IAV) has a segmented genome that allows for the exchange of genome segments between different strains. This reassortment accelerates evolution by breaking linkage, helping IAV cross species barriers to potentially create highly virulent strains. Challenges associated with monitoring the process of reassortment in molecular detail have limited our understanding of its evolutionary implications. We applied a novel deep sequencing approach with quantitative analysis to assess the in vitro temporal evolution of genomic reassortment in IAV. The combination of H1N1 and H3N2 strains reproducibly generated a new H1N2 strain with the hemagglutinin and nucleoprotein segments originating from H1N1 and the remaining six segments from H3N2. By deep sequencing the entire viral genome, we monitored the evolution of reassortment, quantifying the relative abundance of all IAV genome segments from the two parent strains over time and measuring the selection coefficients of the reassorting segments. Additionally, we observed several mutations coemerging with reassortment that were not found during passaging of pure parental IAV strains. Our results demonstrate how reassortment of the segmented genome can accelerate viral evolution in IAV, potentially enabled by the emergence of a small number of individual mutation

Positive Selection Drives Preferred Segment Combinations during Influenza Virus Reassortment

Influenza A virus (IAV) has a segmented genome that allows for the exchange of genome segments between different strains. This reassortment accelerates evolution by breaking linkage, helping IAV cross species barriers to potentially create highly virulent strains. Challenges associated with monitoring the process of reassortment in molecular detail have limited our understanding of its evolutionary implications. We applied a novel deep sequencing approach with quantitative analysis to assess the in vitro temporal evolution of genomic reassortment in IAV. The combination of H1N1 and H3N2 strains reproducibly generated a new H1N2 strain with the hemagglutinin and nucleoprotein segments originating from H1N1 and the remaining six segments from H3N2. By deep sequencing the entire viral genome, we monitored the evolution of reassortment, quantifying the relative abundance of all IAV genome segments from the two parent strains over time and measuring the selection coefficients of the reassorting segments. Additionally, we observed several mutations coemerging with reassortment that were not found during passaging of pure parental IAV strains. Our results demonstrate how reassortment of the segmented genome can accelerate viral evolution in IAV, potentially enabled by the emergence of a small number of individual mutation

Evolution of the influenza A virus genome during development of oseltamivir resistance in vitro

Influenza A virus (IAV) is a major cause of morbidity and mortality throughout the world. Current antiviral therapies include oseltamivir, a neuraminidase inhibitor that prevents the release of nascent viral particles from infected cells. However, the IAV genome can evolve rapidly, and oseltamivir resistance mutations have been detected in numerous clinical samples. Using an in vitro evolution platform and whole-genome population sequencing, we investigated the population genomics of IAV during the development of oseltamivir resistance. Strain A/Brisbane/59/2007 (H1N1) was grown in Madin-Darby canine kidney cells with or without escalating concentrations of oseltamivir over serial passages. Following drug treatment, the H274Y resistance mutation fixed reproducibly within the population. The presence of the H274Y mutation in the viral population, at either a low or a high frequency, led to measurable changes in the neuraminidase inhibition assay. Surprisingly, fixation of the resistance mutation was not accompanied by alterations of viral population diversity or differentiation, and oseltamivir did not alter the selective environment. While the neighboring K248E mutation was also a target of positive selection prior to H274Y fixation, H274Y was the primary beneficial mutation in the population. In addition, once evolved, the H274Y mutation persisted after the withdrawal of the drug, even when not fixed in viral populations. We conclude that only selection of H274Y is required for oseltamivir resistance and that H274Y is not deleterious in the absence of the drug. These collective results could offer an explanation for the recent reproducible rise in oseltamivir resistance in seasonal H1N1 IAV strains in humans

Mutations in influenza A virus neuraminidase and hemagglutinin confer resistance against a broadly neutralizing hemagglutinin stem antibody

Influenza A virus (IAV), a major cause of human morbidity and mortality, continuously evolves in response to selective pressures. Stem-directed, broadly neutralizing antibodies (sBnAbs) targeting influenza hemagglutinin (HA) are a promising therapeutic strategy, but neutralization escape mutants can develop. We used an integrated approach combining viral passaging, deep sequencing, and protein structural analyses to define escape mutations and mechanisms of neutralization escape in vitro for the F10 sBnAb. IAV was propagated with escalating concentrations of F10 over serial passages in cultured cells to select for escape mutations. Viral sequence analysis revealed three mutations in HA and one in neuraminidase (NA). Introduction of these specific mutations into IAV through reverse genetics confirmed their roles in resistance to F10. Structural analyses revealed that the selected HA mutations (S123G, N460S, and N203V) are away from the F10 epitope but may indirectly impact influenza receptor binding, endosomal fusion, or budding. The NA mutation E329K, which was previously identified to be associated with antibody escape, affects the active site of NA, highlighting the importance of the balance between HA and NA function for viral survival. Thus, whole genome population sequencing enables the identification of viral resistance mutations responding to antibody-induced selective pressure.IMPORTANCE Influenza A virus is a public health threat for which currently available vaccines are not always effective. Broadly neutralizing antibodies that bind to the highly-conserved stem region of influenza hemagglutinin (HA) can neutralize many influenza strains. To understand how influenza virus can become resistant or escape such antibodies, we propagated influenza A virus in vitro with escalating concentrations of antibody and analyzed viral populations with whole genome sequencing. We identified HA mutations near and distal to the antibody binding epitope that conferred resistance to antibody neutralization. Additionally, we identified a neuraminidase (NA) mutation that allowed the virus to grow in the presence of high concentrations of the antibody. Virus carrying dual mutations in HA and NA also grew under high antibody concentrations. We show that NA mutations mediate the escape of neutralization by antibodies against HA, highlighting the importance of a balance between HA and NA for optimal virus function

Differential Requirements of Two recA Mutants for Constitutive SOS Expression in Escherichia coli K-12

Background Repairing DNA damage begins with its detection and is often followed by elicitation of a cellular response. In E. coli, RecA polymerizes on ssDNA produced after DNA damage and induces the SOS Response. The RecA-DNA filament is an allosteric effector of LexA auto-proteolysis. LexA is the repressor of the SOS Response. Not all RecA-DNA filaments, however, lead to an SOS Response. Certain recA mutants express the SOS Response (recAC) in the absence of external DNA damage in log phase cells. Methodology/Principal Findings Genetic analysis of two recAC mutants was used to determine the mechanism of constitutive SOS (SOSC) expression in a population of log phase cells using fluorescence of single cells carrying an SOS reporter system (sulAp-gfp). SOSC expression in recA4142 mutants was dependent on its initial level of transcription, recBCD, recFOR, recX, dinI, xthA and the type of medium in which the cells were grown. SOSC expression in recA730 mutants was affected by none of the mutations or conditions tested above. Conclusions/Significance It is concluded that not all recAC alleles cause SOSC expression by the same mechanism. It is hypothesized that RecA4142 is loaded on to a double-strand end of DNA and that the RecA filament is stabilized by the presence of DinI and destabilized by RecX. RecFOR regulate the activity of RecX to destabilize the RecA filament. RecA730 causes SOSC expression by binding to ssDNA in a mechanism yet to be determined

A novel role for RecA under non-stress: promotion of swarming motility in Escherichia coli K-12

BACKGROUND: Bacterial motility is a crucial factor in the colonization of natural environments. Escherichia coli has two flagella-driven motility types: swimming and swarming. Swimming motility consists of individual cell movement in liquid medium or soft semisolid agar, whereas swarming is a coordinated cellular behaviour leading to a collective movement on semisolid surfaces. It is known that swimming motility can be influenced by several types of environmental stress. In nature, environmentally induced DNA damage (e.g. UV irradiation) is one of the most common types of stress. One of the key proteins involved in the response to DNA damage is RecA, a multifunctional protein required for maintaining genome integrity and the generation of genetic variation. RESULTS: The ability of E. coli cells to develop swarming migration on semisolid surfaces was suppressed in the absence of RecA. However, swimming motility was not affected. The swarming defect of a ΔrecA strain was fully complemented by a plasmid-borne recA gene. Although the ΔrecA cells grown on semisolidsurfaces exhibited flagellar production, they also presented impaired individual movement as well as a fully inactive collective swarming migration. Both the comparative analysis of gene expression profiles in wild-type and ΔrecA cells grown on a semisolid surface and the motility of lexA1 [Ind-] mutant cells demonstrated that the RecA effect on swarming does not require induction of the SOS response. By using a RecA-GFP fusion protein we were able to segregate the effect of RecA on swarming from its other functions. This protein fusion failed to regulate the induction of the SOS response, the recombinational DNA repair of UV-treated cells and the genetic recombination, however, it was efficient in rescuing the swarming motility defect of the ΔrecA mutant. The RecA-GFP protein retains a residual ssDNA-dependent ATPase activity but does not perform DNA strand exchange. CONCLUSION: The experimental evidence presented in this work supports a novel role for RecA: the promotion of swarming motility. The defective swarming migration of ΔrecA cells does not appear to be associated with defective flagellar production; rather, it seems to be associated with an abnormal flagellar propulsion function. Our results strongly suggest that the RecA effect on swarming motility does not require an extensive canonical RecA nucleofilament formation. RecA is the first reported cellular factor specifically affecting swarming but not swimming motility in E. coli. The integration of two apparently disconnected biologically important processes, such as the maintenance of genome integrity and motility in a unique protein, may have important evolutive consequences

Separation of Recombination and SOS Response in Escherichia coli RecA Suggests LexA Interaction Sites

RecA plays a key role in homologous recombination, the induction of the DNA damage response through LexA cleavage and the activity of error-prone polymerase in Escherichia coli. RecA interacts with multiple partners to achieve this pleiotropic role, but the structural location and sequence determinants involved in these multiple interactions remain mostly unknown. Here, in a first application to prokaryotes, Evolutionary Trace (ET) analysis identifies clusters of evolutionarily important surface amino acids involved in RecA functions. Some of these clusters match the known ATP binding, DNA binding, and RecA-RecA homo-dimerization sites, but others are novel. Mutation analysis at these sites disrupted either recombination or LexA cleavage. This highlights distinct functional sites specific for recombination and DNA damage response induction. Finally, our analysis reveals a composite site for LexA binding and cleavage, which is formed only on the active RecA filament. These new sites can provide new drug targets to modulate one or more RecA functions, with the potential to address the problem of evolution of antibiotic resistance at its root

A Wide Extent of Inter-Strain Diversity in Virulent and Vaccine Strains of Alphaherpesviruses

Alphaherpesviruses are widespread in the human population, and include herpes simplex virus 1 (HSV-1) and 2, and varicella zoster virus (VZV). These viral pathogens cause epithelial lesions, and then infect the nervous system to cause lifelong latency, reactivation, and spread. A related veterinary herpesvirus, pseudorabies (PRV), causes similar disease in livestock that result in significant economic losses. Vaccines developed for VZV and PRV serve as useful models for the development of an HSV-1 vaccine. We present full genome sequence comparisons of the PRV vaccine strain Bartha, and two virulent PRV isolates, Kaplan and Becker. These genome sequences were determined by high-throughput sequencing and assembly, and present new insights into the attenuation of a mammalian alphaherpesvirus vaccine strain. We find many previously unknown coding differences between PRV Bartha and the virulent strains, including changes to the fusion proteins gH and gB, and over forty other viral proteins. Inter-strain variation in PRV protein sequences is much closer to levels previously observed for HSV-1 than for the highly stable VZV proteome. Almost 20% of the PRV genome contains tandem short sequence repeats (SSRs), a class of nucleic acids motifs whose length-variation has been associated with changes in DNA binding site efficiency, transcriptional regulation, and protein interactions. We find SSRs throughout the herpesvirus family, and provide the first global characterization of SSRs in viruses, both within and between strains. We find SSR length variation between different isolates of PRV and HSV-1, which may provide a new mechanism for phenotypic variation between strains. Finally, we detected a small number of polymorphic bases within each plaque-purified PRV strain, and we characterize the effect of passage and plaque-purification on these polymorphisms. These data add to growing evidence that even plaque-purified stocks of stable DNA viruses exhibit limited sequence heterogeneity, which likely seeds future strain evolution

Specific capture and whole-genome sequencing of viruses from clinical samples.

Whole genome sequencing of viruses directly from clinical samples is integral for understanding the genetics of host-virus interactions. Here, we report the use of sample sparing target enrichment (by hybridisation) for viral nucleic acid separation and deep-sequencing of herpesvirus genomes directly from a range of clinical samples including saliva, blood, virus vesicles, cerebrospinal fluid, and tumour cell lines. We demonstrate the effectiveness of the method by deep-sequencing 13 highly cell-associated human herpesvirus genomes and generating full length genome alignments at high read depth. Moreover, we show the specificity of the method enables the study of viral population structures and their diversity within a range of clinical samples types