Genetic and phenotypic changes to Venezuelan equine encephalitis virus following treatment with β-D-N4-hydroxycytidine, an RNA mutagen

The high mutation rate of RNA viruses provides viral populations with the ability to adapt to new environments but also makes them vulnerable to extinction due to the deleterious effects of mutations, which is the conceptual basis for the antiviral activity of RNA mutagens. However, there are still gaps in the quantitative understanding of the dynamics between the mutations induced by an RNA mutagen and its effects on viral fitness. To address this, we used Venezuelan Equine Encephalitis Virus (VEEV) and the potent RNA mutagen β-d-N4-hydroxycytidine (NHC) as a model to analyze virus replication competency and mutation frequency following treatment in the total and replication-competent viral populations separately. We found that NHC induced transition mutations in a concentration dependent manner in the total population, while the replication-competent population maintained itself within an increased, yet narrow, mutation spectrum. The incorporation of NHC mainly happened during the positive sense RNA synthesis of VEEV. A growth kinetic analysis of VEEV population treated with NHC pointed to a lower but more diverse distribution in mutational fitness, demonstrating that NHC-induced mutations negatively and broadly affect the fitness of the virus. Together, our study provides mechanistic insight into how RNA mutagens affect viral population landscape and the potential of RNA mutagens as an antiviral strategy for alphaviruses. Each sample is a PCR-amplified region of TC-83 genome encoding the non-structural proteins 2, 3, and 4. The samples were derived from viruses that were either (1) untreated, (2) treated with an RNA mutatgen, or (3) treated both with a direct antiviral and RNA mutatgen. The majority of the samples were from individually isolated virus clones (Nanopore sequencing), while the samples sequenced by PacBio were from TC-83 virus stock (heterogeneous population).