Thomas S. Mehoke

Flexibility in surface-exposed loops in a virus capsid mediates escape from antibody neutralization.

ABSTRACT New human norovirus strains emerge every 2 to 3 years, partly due to mutations in the viral capsid that allow escape from antibody neutralization and herd immunity. To understand how noroviruses evolve antibody resistance, we investigated the structural basis for the escape of murine norovirus (MNV) from antibody neutralization. To identify specific residues in the MNV-1 protruding (P) domain of the capsid that play a role in escape from the neutralizing monoclonal antibody (MAb) A6.2, 22 recombinant MNVs were generated with amino acid substitutions in the A′B′ and E′F′ loops. Six mutations in the E′F′ loop (V378F, A382K, A382P, A382R, D385G, and L386F) mediated escape from MAb A6.2 neutralization. To elucidate underlying structural mechanisms for these results, the atomic structure of the A6.2 Fab was determined and fitted into the previously generated pseudoatomic model of the A6.2 Fab/MNV-1 virion complex. Previously, two distinct conformations, A and B, of the atomic structures of the MNV-1 P domain were identified due to flexibility in the two P domain loops. A superior stereochemical fit of the A6.2 Fab to the A conformation of the MNV P domain was observed. Structural analysis of our observed escape mutants indicates changes toward the less-preferred B conformation of the P domain. The shift in the structural equilibrium of the P domain toward the conformation with poor structural complementarity to the antibody strongly supports a unique mechanism for antibody escape that occurs via antigen flexibility instead of direct antibody-antigen binding. IMPORTANCE Human noroviruses cause the majority of all nonbacterial gastroenteritis worldwide. New epidemic strains arise in part by mutations in the viral capsid leading to escape from antibody neutralization. Herein, we identify a series of point mutations in a norovirus capsid that mediate escape from antibody neutralization and determine the structure of a neutralizing antibody. Fitting of the antibody structure into the virion/antibody complex identifies two conformations of the antibody binding domain of the viral capsid: one with a superior fit and the other with an inferior fit to the antibody. These data suggest a unique mode of antibody neutralization. In contrast to other viruses that largely escape antibody neutralization through direct disruption of the antibody-virus interface, we identify mutations that acted indirectly by limiting the conformation of the antibody binding loop in the viral capsid and drive the antibody binding domain into the conformation unable to be bound by the antibody.

Isolation and Analysis of Rare Norovirus Recombinants from Coinfected Mice Using Drop-Based Microfluidics

ABSTRACT Human noroviruses (HuNoVs) are positive-sense RNA viruses that can cause severe, highly infectious gastroenteritis. HuNoV outbreaks are frequently associated with recombination between circulating strains. Strain genotyping and phylogenetic analyses show that noroviruses often recombine in a highly conserved region near the junction of the viral polyprotein (open reading frame 1 [ORF1]) and capsid (ORF2) genes and occasionally within the RNA-dependent RNA polymerase (RdRP) gene. Although genotyping methods are useful for tracking changes in circulating viral populations, they report only the dominant recombinant strains and do not elucidate the frequency or range of recombination events. Furthermore, the relatively low frequency of recombination in RNA viruses has limited studies to cell culture or in vitro systems, which do not reflect the complexities and selective pressures present in an infected organism. Using two murine norovirus (MNV) strains to model coinfection, we developed a microfluidic platform to amplify, detect, and recover individual recombinants following in vitro and in vivo coinfection. One-step reverse transcriptase PCR (RT-PCR) was performed in picoliter drops with primers that identified the wild-type and recombinant progenies and scanned for recombination breakpoints at ∼1-kb intervals. We detected recombination between MNV strains at multiple loci spanning the viral protease, RdRP, and capsid ORFs and isolated individual recombinant RNA genomes that were present at a frequency of 1/300,000 or higher. This study is the first to examine norovirus recombination following coinfection of an animal and suggests that the exchange of RNA among viral genomes in an infected host occurs in multiple locations and is an important driver of genetic diversity. IMPORTANCE RNA viruses increase diversity and escape host immune barriers by genomic recombination. Studies using a number of viral systems indicate that recombination occurs via template switching by the virus-encoded RNA-dependent RNA polymerase (RdRP). However, factors that govern the frequency and positions of recombination in an infected organism remain largely unknown. This work leverages advances in the applied physics of drop-based microfluidics to isolate and sequence rare recombinants arising from the coinfection of mice with two distinct strains of murine norovirus. This study is the first to detect and analyze norovirus recombination in an animal model.

A high-throughput drop microfluidic system for virus culture and analysis

High mutation rates and short replication times lead to rapid evolution in RNA viruses. New tools for high-throughput culture and analysis of viral phenotypes will enable more effective studies of viral evolutionary processes. A water-in-oil drop microfluidic system to study virus-cell interactions at the single event level on a massively parallel scale is described here. Murine norovirus (MNV-1) particles were co-encapsulated with individual RAW 264.7 cells in 65 pL aqueous drops formed by flow focusing in 50 μm microchannels. At low multiplicity of infection (MOI), viral titers increased greatly, reaching a maximum 18 h post-encapsulation. This system was employed to evaluate MNV-1 escape from a neutralizing monoclonal antibody (clone A6.2). Further, the system was validated as a means for testing escape from antibody neutralization using a series of viral point mutants. Finally, the replicative capacity of single viral particles in drops under antibody stress was tested. Under standard conditions, many RNA virus stocks harbor minority populations of genotypic and phenotypic variants, resulting in quasispecies. These data show that when single cells are encapsulated with single viral particles under antibody stress without competition from other virions, the number of resulting infectious particles is nearly equivalent to the number of viral genomes present. These findings suggest that lower fitness virions can infect cells successfully and replicate, indicating that the microfluidics system may serve as an effective tool for isolating mutants that escape evolutionary stressors.

Murine norovirus (MNV-1) exposure in vitro to the purine nucleoside analog Ribavirin increases quasispecies diversity

Ribavirin is a pharmaceutical antiviral used for the treatment of RNA virus infections including norovirus, hepatitis C virus, hepatitis E virus, Lassa virus, respiratory syncytial virus, and rhinovirus. Despite the drugs history and documented efficacy, the antiviral mechanism of Ribavirin remains unclear. Mechanisms proposed include depletion of the intracellular GTP pool, immunomodulatory effects, induction of error catastrophe, inhibition of viral polymerase activity, and/or inhibition of viral capping. In the present study, we leveraged deep sequencing data to demonstrate that Ribavirin increases murine norovirus (MNV-1) viral diversity. By serial passaging MNV-1 in RAW 264.7 cells for twenty generations in the presence of Ribavirin, we demonstrated statistically significant increases in both the number of unique haplotypes and the average pairwise difference (APD). Based on statistically significant differences in the probability of nucleotide mutations based on Roche 454 sequencing, we also demonstrated that single nucleotide substitutions are increased in the presence of Ribavirin. Finally, we demonstrated Ribavirins impact on statistically significantly reducing the relative proportion of the dominant sequence within the quasispecies.

Tuning the course of evolution on the biophysical fitness landscape of an RNA virus

Predicting viral evolution remains a major challenge with profound implications for public health. Viral evolutionary pathways are determined by the fitness landscape, which maps viral genotype to fitness. However, a quantitative description of the landscape and the evolutionary forces on it remain elusive. Here, we apply a biophysical fitness model based on capsid folding stability and antibody binding affinity to predict the evolutionary pathway of norovirus escaping a neutralizing antibody. The model is validated by experimental evolution in bulk culture and in a drop-based microfluidics device, the “Evolution Chip”, which propagates millions of independent viral sub-populations. We demonstrate that along the axis of binding affinity, selection for escape variants and drift due to random mutations have the same direction. However, along folding stability, selection and drift are opposing forces whose balance is tuned by viral population size. Our results demonstrate that predictable epistatic tradeoffs shape viral evolution.

Statistical analysis of co-occurrence patterns in microbial presence-absence datasets

Drawing on a long history in macroecology, correlation analysis of microbiome datasets is becoming a common practice for identifying relationships or shared ecological niches among bacterial taxa. However, many of the statistical issues that plague such analyses in macroscale communities remain unresolved for microbial communities. Here, we discuss problems in the analysis of microbial species correlations based on presence-absence data. We focus on presence-absence data because this information is more readily obtainable from sequencing studies, especially for whole-genome sequencing, where abundance estimation is still in its infancy. First, we show how Pearson’s correlation coefficient (r) and Jaccard’s index (J)–two of the most common metrics for correlation analysis of presence-absence data–can contradict each other when applied to a typical microbiome dataset. In our dataset, for example, 14% of species-pairs predicted to be significantly correlated by r were not predicted to be significantly correlated using J, while 37.4% of species-pairs predicted to be significantly correlated by J were not predicted to be significantly correlated using r. Mismatch was particularly common among species-pairs with at least one rare species (<10% prevalence), explaining why r and J might differ more strongly in microbiome datasets, where there are large numbers of rare taxa. Indeed 74% of all species-pairs in our study had at least one rare species. Next, we show how Pearson’s correlation coefficient can result in artificial inflation of positive taxon relationships and how this is a particular problem for microbiome studies. We then illustrate how Jaccard’s index of similarity (J) can yield improvements over Pearson’s correlation coefficient. However, the standard null model for Jaccard’s index is flawed, and thus introduces its own set of spurious conclusions. We thus identify a better null model based on a hypergeometric distribution, which appropriately corrects for species prevalence. This model is available from recent statistics literature, and can be used for evaluating the significance of any value of an empirically observed Jaccard’s index. The resulting simple, yet effective method for handling correlation analysis of microbial presence-absence datasets provides a robust means of testing and finding relationships and/or shared environmental responses among microbial taxa.

Evolution on the Biophysical Fitness Landscape of an RNA Virus

Abstract Viral evolutionary pathways are determined by the fitness landscape, which maps viral genotype to fitness. However, a quantitative description of the landscape and the evolutionary forces on it remain elusive. Here, we apply a biophysical fitness model based on capsid folding stability and antibody binding affinity to predict the evolutionary pathway of norovirus escaping a neutralizing antibody. The model is validated by experimental evolution in bulk culture and in a drop-based microfluidics that propagates millions of independent small viral subpopulations. We demonstrate that along the axis of binding affinity, selection for escape variants and drift due to random mutations have the same direction, an atypical case in evolution. However, along folding stability, selection and drift are opposing forces whose balance is tuned by viral population size. Our results demonstrate that predictable epistatic tradeoffs between molecular traits of viral proteins shape viral evolution.