Emergence and Magnitude of ML336 Resistance in Venezuelan Equine Encephalitis Virus Depend on the Microenvironment

RNA viruses, including Venezuelan equine encephalitis virus (VEEV), have high mutation rates that allow for rapid adaptation to selective pressures in their environment. Antiviral compounds exert one such pressure on virus populations during infections. Next-generation sequencing allows for examination of viruses at the population level, which enables tracking of low levels of single-nucleotide polymorphisms in the population over time. Therefore, the timing and extent of the emergence of resistance to antivirals can be tracked and assessed. We show here that in VEEV, the trajectory and penetration of antiviral resistance reflected the microenvironment in which the virus population replicates. In summary, we show the diversity of VEEV within a single population under antiviral pressure and two distinct cell types, and we show that population dynamics in these viruses can be examined to better understand how they evolve over time. ABSTRACT Venezuelan equine encephalitis virus (VEEV) is a New World Alphavirus that can cause neurological disease and death in humans and equines following transmission from infected mosquitoes. Despite the continued epidemic threat of VEEV, and its potential use as a bioterrorism agent, there are no FDA-approved antivirals or vaccines for treatment or prevention. Previously, we reported the discovery of a small molecule, ML336, with potent antiviral activity against VEEV. To further explore the population-level resistance profiles of ML336, we developed a whole-genome next-generation sequencing (NGS) approach to examine single nucleotide polymorphisms (SNPs) from virus passaged in dose escalation studies in a nonhuman primate kidney epithelial and a human astrocyte cell line, Vero 76 and SVGA, respectively. We passaged VEEV TC-83 in these two cell lines over seven concentrations of ML336, starting at 50 nM. NGS revealed several prominent mutations in the nonstructural protein (nsP) 3 and nsP4 genes that emerged consistently in these two distinct in vitro environments—notably, a mutation at Q210 in nsP4. Several of these mutations were stable following passaging in the absence of ML336 in Vero 76 cells. Network analyses showed that the trajectory of resistance differed between Vero and SVGA. Moreover, the penetration of SNPs was lower in SVGA. In conclusion, we show that the microenvironment influenced the SNP profile of VEEV TC-83. Understanding the dynamics of resistance in VEEV against newly developed antiviral compounds will guide the design of optimal drug candidates and dosing regimens for minimizing the emergence of resistant viruses. IMPORTANCE RNA viruses, including Venezuelan equine encephalitis virus (VEEV), have high mutation rates that allow for rapid adaptation to selective pressures in their environment. Antiviral compounds exert one such pressure on virus populations during infections. Next-generation sequencing allows for examination of viruses at the population level, which enables tracking of low levels of single-nucleotide polymorphisms in the population over time. Therefore, the timing and extent of the emergence of resistance to antivirals can be tracked and assessed. We show here that in VEEV, the trajectory and penetration of antiviral resistance reflected the microenvironment in which the virus population replicates. In summary, we show the diversity of VEEV within a single population under antiviral pressure and two distinct cell types, and we show that population dynamics in these viruses can be examined to better understand how they evolve over time.