The E2 glycoprotein holds key residues for Mayaro virus adaptation to the urban Aedes aegypti mosquito.

Adaptation to mosquito vectors suited for transmission in urban settings is a major driver in the emergence of arboviruses. To better anticipate future emergence events, it is crucial to assess their potential to adapt to new vector hosts. In this work, we used two different experimental evolution approaches to study the adaptation process of an emerging alphavirus, Mayaro virus (MAYV), to Ae. aegypti, an urban mosquito vector of many other arboviruses. We identified E2-T179N as a key mutation increasing MAYV replication in insect cells and enhancing transmission after escaping the midgut of live Ae. aegypti. In contrast, this mutation decreased viral replication and binding in human fibroblasts, a primary cellular target of MAYV in humans. We also showed that MAYV E2-T179N generates reduced viremia and displays less severe tissue pathology in vivo in a mouse model. We found evidence in mouse fibroblasts that MAYV E2-T179N is less dependent on the Mxra8 receptor for replication than WT MAYV. Similarly, exogenous expression of human apolipoprotein receptor 2 and Mxra8 enhanced WT MAYV replication compared to MAYV E2-T179N. When this mutation was introduced in the closely related chikungunya virus, which has caused major outbreaks globally in the past two decades, we observed increased replication in both human and insect cells, suggesting E2 position 179 is an important determinant of alphavirus host-adaptation, although in a virus-specific manner. Collectively, these results indicate that adaptation at the T179 residue in MAYV E2 may result in increased vector competence-but coming at the cost of optimal replication in humans-and may represent a first step towards a future emergence event.

RNA-mediated interference and reverse transcription control the persistence of RNA viruses in the insect model Drosophila.

How persistent viral infections are established and maintained is widely debated and remains poorly understood. We found here that the persistence of RNA viruses in Drosophila melanogaster was achieved through the combined action of cellular reverse-transcriptase activity and the RNA-mediated interference (RNAi) pathway. Fragments of diverse RNA viruses were reverse-transcribed early during infection, which resulted in DNA forms embedded in retrotransposon sequences. Those virus-retrotransposon DNA chimeras produced transcripts processed by the RNAi machinery, which in turn inhibited viral replication. Conversely, inhibition of reverse transcription hindered the appearance of chimeric DNA and prevented persistence. Our results identify a cooperative function for retrotransposons and antiviral RNAi in the control of lethal acute infection for the establishment of viral persistence.

APOBEC3C S188I Polymorphism Enhances Context-Specific Editing of Hepatitis B Virus Genome.

Single-nucleotide polymorphism in APOBEC3C (resulting in a serine to isoleucine in position 188) is present in approximately 10% of African populations and greatly enhances restriction against human immunodeficiency virus-1 and simian immunodeficiency virus by improving dimerization and DNA processivity of the enzyme. In this study, we demonstrated in culture and in infected patients that hepatitis B virus (HBV) could be edited by APOBEC3CS188I. Using next-generation sequencing, we demonstrated that APOBEC3CS188I led to enhanced editing activity in 5'TpCpA→5'TpTpA context. This constitutes a new hallmark of this enzyme, which could be used to determine its impact on HBV or nuclear DNA.

Host nutritional status affects alphavirus virulence, transmission, and evolution.

Malnourishment, specifically overweight/obesity and undernourishment, affects more than 2.5 billion people worldwide, with the number affected ever-increasing. Concurrently, emerging viral diseases, particularly those that are mosquito-borne, have spread dramatically in the past several decades, culminating in outbreaks of several viruses worldwide. Both forms of malnourishment are known to lead to an aberrant immune response, which can worsen disease outcomes and reduce vaccination efficacy for viral pathogens such as influenza and measles. Given the increasing rates of malnutrition and spread of arthropod-borne viruses (arboviruses), there is an urgent need to understand the role of host nutrition on the infection, virulence, and transmission of these viruses. To address this gap in knowledge, we infected lean, obese, and undernourished mice with arthritogenic arboviruses from the genus Alphavirus and assessed morbidity, virus replication, transmission, and evolution. Obesity and undernourishment did not consistently influence virus replication in the blood of infected animals except for reductions in virus in obese mice late in infection. However, morbidity was increased in obese mice under all conditions. Using Mayaro virus (MAYV) as a model arthritogenic alphavirus, we determined that both obese and undernourished mice transmit virus less efficiently to mosquitoes than control (lean) mice. In addition, viral genetic diversity and replicative fitness were reduced in virus isolated from obese compared to lean controls. Taken together, nutrition appears to alter the course of alphavirus infection and should be considered as a critical environmental factor during outbreaks.

Chikungunya Virus Vaccine Candidates with Decreased Mutational Robustness Are Attenuated In Vivo and Have Compromised Transmissibility.

Chikungunya virus (CHIKV) is a reemerged arbovirus, a member of the Togaviridae family. It circulates through mosquito vectors mainly of the Aedes family and a mammalian host. CHIKV causes chikungunya fever, a mild to severe disease characterized by arthralgia, with some fatal outcomes described. In the past years, several outbreaks mainly caused by enhanced adaptation of the virus to the vector and ineffective control of the contacts between infected mosquito populations and the human host have been reported. Vaccines represent the best solution for the control of insect-borne viruses, including CHIKV, but are often unavailable. We designed live attenuated CHIKVs by applying a rational genomic design based on multiple replacements of synonymous codons. In doing so, the virus mutational robustness (capacity to maintain phenotype despite introduction of mutations to genotype) is decreased, driving the viral population toward deleterious evolutionary trajectories. When the candidate viruses were tested in the insect and mammalian hosts, we observed overall strong attenuation in both and greatly diminished signs of disease. Moreover, we found that the vaccine candidates elicited protective immunity related to the production of neutralizing antibodies after a single dose. During an experimental transmission cycle between mosquitoes and naive mice, vaccine candidates could be transmitted by mosquito bite, leading to asymptomatic infection in mice with compromised dissemination. Using deep-sequencing technology, we observed an increase in detrimental (stop) codons, which confirmed the effectiveness of this genomic design. Because the approach involves hundreds of synonymous modifications to the genome, the reversion risk is significantly reduced, rendering the viruses promising vaccine candidates.IMPORTANCE Chikungunya fever is a debilitating disease that causes severe pain to the joints, which can compromise the patient's lifestyle for several months and even in some grave cases lead to death. The etiological agent is chikungunya virus, an alphavirus transmitted by mosquito bite. Currently, there are no approved vaccines or treatments against the disease. In our research, we developed novel live attenuated vaccine candidates against chikungunya virus by applying an innovative genomic design. When tested in the insect and mammalian host, the vaccine candidates did not cause disease, elicited strong protection against further infection, and had low risk of reversion to pathogenic phenotypes.

piRNA pathway is not required for antiviral defense in Drosophila melanogaster.

Since its discovery, RNA interference has been identified as involved in many different cellular processes, and as a natural antiviral response in plants, nematodes, and insects. In insects, the small interfering RNA (siRNA) pathway is the major antiviral response. In recent years, the Piwi-interacting RNA (piRNA) pathway also has been implicated in antiviral defense in mosquitoes infected with arboviruses. Using Drosophila melanogaster and an array of viruses that infect the fruit fly acutely or persistently or are vertically transmitted through the germ line, we investigated in detail the extent to which the piRNA pathway contributes to antiviral defense in adult flies. Following virus infection, the survival and viral titers of Piwi, Aubergine, Argonaute-3, and Zucchini mutant flies were similar to those of wild type flies. Using next-generation sequencing of small RNAs from wild type and siRNA mutant flies, we showed that no viral-derived piRNAs were produced in fruit flies during different types of viral infection. Our study provides the first evidence, to our knowledge, that the piRNA pathway does not play a major role in antiviral defense in adult Drosophila and demonstrates that viral-derived piRNA production depends on the biology of the host-virus combination rather than being part of a general antiviral process in insects.

Uncovering the Repertoire of Endogenous Flaviviral Elements in Aedes Mosquito Genomes.

Endogenous viral elements derived from nonretroviral RNA viruses have been described in various animal genomes. Whether they have a biological function, such as host immune protection against related viruses, is a field of intense study. Here, we investigated the repertoire of endogenous flaviviral elements (EFVEs) in Aedes mosquitoes, the vectors of arboviruses such as dengue and chikungunya viruses. Previous studies identified three EFVEs from Aedes albopictus cell lines and one from Aedes aegypti cell lines. However, an in-depth characterization of EFVEs in wild-type mosquito populations and individual mosquitoes in vivo has not been performed. We detected the full-length DNA sequence of the previously described EFVEs and their respective transcripts in several A. albopictus and A. aegypti populations from geographically distinct areas. However, EFVE-derived proteins were not detected by mass spectrometry. Using deep sequencing, we detected the production of PIWI-interacting RNA-like small RNAs, in an antisense orientation, targeting the EFVEs and their flanking regions in vivo The EFVEs were integrated in repetitive regions of the mosquito genomes, and their flanking sequences varied among mosquito populations. We bioinformatically predicted several new EFVEs from a Vietnamese A. albopictus population and observed variation in the occurrence of those elements among mosquitoes. Phylogenetic analysis of an A. aegypti EFVE suggested that it integrated prior to the global expansion of the species and subsequently diverged among and within populations. The findings of this study together reveal the substantial structural and nucleotide diversity of flaviviral integrations in Aedes genomes. Unraveling this diversity will help to elucidate the potential biological function of these EFVEs.IMPORTANCE Endogenous viral elements (EVEs) are whole or partial viral sequences integrated in host genomes. Interestingly, some EVEs have important functions for host fitness and antiviral defense. Because mosquitoes also have EVEs in their genomes, characterizing these EVEs is a prerequisite for their potential use to manipulate the mosquito antiviral response. In the study described here, we focused on EVEs related to the Flavivirus genus, to which dengue and Zika viruses belong, in individual Aedes mosquitoes from geographically distinct areas. We show the existence in vivo of flaviviral EVEs previously identified in mosquito cell lines, and we detected new ones. We show that EVEs have evolved differently in each mosquito population. They produce transcripts and small RNAs but not proteins, suggesting a function at the RNA level. Our study uncovers the diverse repertoire of flaviviral EVEs in Aedes mosquito populations and contributes to an understanding of their role in the host antiviral system.

Understanding the Mechanism of the Broad-Spectrum Antiviral Activity of Favipiravir (T-705): Key Role of the F1 Motif of the Viral Polymerase.

Favipiravir (T-705) is a broad-spectrum antiviral agent that has been approved in Japan for the treatment of influenza virus infections. T-705 also inhibits the replication of various RNA viruses, including chikungunya virus (CHIKV). We demonstrated earlier that the K291R mutation in the F1 motif of the RNA-dependent RNA polymerase (RdRp) of CHIKV is responsible for low-level resistance to T-705. Interestingly, this lysine is highly conserved in the RdRp of positive-sense single-stranded RNA (+ssRNA) viruses. To obtain insights into the unique broad-spectrum antiviral activity of T-705, we explored the role of this lysine using another +ssRNA virus, namely, coxsackievirus B3 (CVB3). Introduction of the corresponding K-to-R substitution in the CVB3 RdRp (K159R) resulted in a nonviable virus. Replication competence of the K159R variant was restored by spontaneous acquisition of an A239G substitution in the RdRp. A mutagenesis analysis at position K159 identified the K159M variant as the only other viable variant which had also acquired the A239G substitution. The K159 substitutions markedly decreased the processivity of the purified viral RdRp, which was restored by the introduction of the A239G mutation. The K159R A239G and K159M A239G variants proved, surprisingly, more susceptible than the wild-type virus to T-705 and exhibited lower fidelity in polymerase assays. Furthermore, the K159R A239G variant was found to be highly attenuated in mice. We thus demonstrate that the conserved lysine in the F1 motif of the RdRp of +ssRNA viruses is involved in the broad-spectrum antiviral activity of T-705 and that it is a key amino acid for the proper functioning of the enzyme.IMPORTANCE In this study, we report the key role of a highly conserved lysine residue of the viral polymerase in the broad-spectrum antiviral activity of favipiravir (T-705) against positive-sense single-stranded RNA viruses. Substitutions of this conserved lysine have a major negative impact on the functionality of the RdRp. Furthermore, we show that this lysine is involved in the fidelity of the RdRp and that the RdRp fidelity influences the sensitivity of the virus for the antiviral efficacy of T-705. Consequently, these results provide insights into the mechanism of the antiviral activity of T-705 and may lay the basis for the design of novel chemical scaffolds that may be endowed with a more potent broad-spectrum antiviral activity than that of T-705.

Large-Scale Phylogenomic Analysis Reveals the Complex Evolutionary History of Rabies Virus in Multiple Carnivore Hosts.

The natural evolution of rabies virus (RABV) provides a potent example of multiple host shifts and an important opportunity to determine the mechanisms that underpin viral emergence. Using 321 genome sequences spanning an unprecedented diversity of RABV, we compared evolutionary rates and selection pressures in viruses sampled from multiple primary host shifts that occurred on various continents. Two major phylogenetic groups, bat-related RABV and dog-related RABV, experiencing markedly different evolutionary dynamics were identified. While no correlation between time and genetic divergence was found in bat-related RABV, the evolution of dog-related RABV followed a generally clock-like structure, although with a relatively low evolutionary rate. Subsequent molecular clock dating indicated that dog-related RABV likely underwent a rapid global spread following the intensification of intercontinental trade starting in the 15th century. Strikingly, although dog RABV has jumped to various wildlife species from the order Carnivora, we found no clear evidence that these host-jumping events involved adaptive evolution, with RABV instead characterized by strong purifying selection, suggesting that ecological processes also play an important role in shaping patterns of emergence. However, specific amino acid changes were associated with the parallel emergence of RABV in ferret-badgers in Asia, and some host shifts were associated with increases in evolutionary rate, particularly in the ferret-badger and mongoose, implying that changes in host species can have important impacts on evolutionary dynamics.

Homology-Based Identification of a Mutation in the Coronavirus RNA-Dependent RNA Polymerase That Confers Resistance to Multiple Mutagens.

UNLABELLED: Positive-sense RNA viruses encode RNA-dependent RNA polymerases (RdRps) essential for genomic replication. With the exception of the large nidoviruses, such as coronaviruses (CoVs), RNA viruses lack proofreading and thus are dependent on RdRps to control nucleotide selectivity and fidelity. CoVs encode a proofreading exonuclease in nonstructural protein 14 (nsp14-ExoN), which confers a greater-than-10-fold increase in fidelity compared to other RNA viruses. It is unknown to what extent the CoV polymerase (nsp12-RdRp) participates in replication fidelity. We sought to determine whether homology modeling could identify putative determinants of nucleotide selectivity and fidelity in CoV RdRps. We modeled the CoV murine hepatitis virus (MHV) nsp12-RdRp structure and superimposed it on solved picornaviral RdRp structures. Fidelity-altering mutations previously identified in coxsackie virus B3 (CVB3) were mapped onto the nsp12-RdRp model structure and then engineered into the MHV genome with [nsp14-ExoN(+)] or without [nsp14-ExoN(-)] ExoN activity. Using this method, we identified two mutations conferring resistance to the mutagen 5-fluorouracil (5-FU): nsp12-M611F and nsp12-V553I. For nsp12-V553I, we also demonstrate resistance to the mutagen 5-azacytidine (5-AZC) and decreased accumulation of mutations. Resistance to 5-FU, and a decreased number of genomic mutations, was effectively masked by nsp14-ExoN proofreading activity. These results indicate that nsp12-RdRp likely functions in fidelity regulation and that, despite low sequence conservation, some determinants of RdRp nucleotide selectivity are conserved across RNA viruses. The results also indicate that, with regard to nucleotide selectivity, nsp14-ExoN is epistatic to nsp12-RdRp, consistent with its proposed role in a multiprotein replicase-proofreading complex. IMPORTANCE: RNA viruses have evolutionarily fine-tuned replication fidelity to balance requirements for genetic stability and diversity. Responsibility for replication fidelity in RNA viruses has been attributed to the RNA-dependent RNA polymerases, with mutations in RdRps for multiple RNA viruses shown to alter fidelity and attenuate virus replication and virulence. Coronaviruses (CoVs) are the only known RNA viruses to encode a proofreading exonuclease (nsp14-ExoN), as well as other replicase proteins involved in regulation of fidelity. This report shows that the CoV RdRp (nsp12) likely functions in replication fidelity; that residue determinants of CoV RdRp nucleotide selectivity map to similar structural regions of other, unrelated RNA viral polymerases; and that for CoVs, the proofreading activity of the nsp14-ExoN is epistatic to the function of the RdRp in fidelity.

Increasing Clinical Severity during a Dengue Virus Type 3 Cuban Epidemic: Deep Sequencing of Evolving Viral Populations.

UNLABELLED: During the dengue virus type 3 (DENV-3) epidemic that occurred in Havana in 2001 to 2002, severe disease was associated with the infection sequence DENV-1 followed by DENV-3 (DENV-1/DENV-3), while the sequence DENV-2/DENV-3 was associated with mild/asymptomatic infections. To determine the role of the virus in the increasing severity demonstrated during the epidemic, serum samples collected at different time points were studied. A total of 22 full-length sequences were obtained using a deep-sequencing approach. Bayesian phylogenetic analysis of consensus sequences revealed that two DENV-3 lineages were circulating in Havana at that time, both grouped within genotype III. The predominant lineage is closely related to Peruvian and Ecuadorian strains, while the minor lineage is related to Venezuelan strains. According to consensus sequences, relatively few nonsynonymous mutations were observed; only one was fixed during the epidemic at position 4380 in the NS2B gene. Intrahost genetic analysis indicated that a significant minor population was selected and became predominant toward the end of the epidemic. In conclusion, greater variability was detected during the epidemic's progression in terms of significant minority variants, particularly in the nonstructural genes. An increasing trend of genetic diversity toward the end of the epidemic was observed only for synonymous variant allele rates, with higher variability in secondary cases. Remarkably, significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in the structural proteins premembrane (PrM) and envelope (E). Therefore, the dynamic of evolving viral populations in the context of heterotypic antibodies could be related to the increasing clinical severity observed during the epidemic. IMPORTANCE: Based on the evidence that DENV fitness is context dependent, our research has focused on the study of viral factors associated with intraepidemic increasing severity in a unique epidemiological setting. Here, we investigated the intrahost genetic diversity in acute human samples collected at different time points during the DENV-3 epidemic that occurred in Cuba in 2001 to 2002 using a deep-sequencing approach. We concluded that greater variability in significant minor populations occurred as the epidemic progressed, particularly in the nonstructural genes, with higher variability observed in secondary infection cases. Remarkably, for the first time significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in structural proteins. These findings indicate that high-resolution approaches are needed to unravel molecular mechanisms involved in dengue pathogenesis.

Group Selection and Contribution of Minority Variants during Virus Adaptation Determines Virus Fitness and Phenotype.

Understanding how a pathogen colonizes and adapts to a new host environment is a primary aim in studying emerging infectious diseases. Adaptive mutations arise among the thousands of variants generated during RNA virus infection, and identifying these variants will shed light onto how changes in tropism and species jumps can occur. Here, we adapted Coxsackie virus B3 to a highly permissive and less permissive environment. Using deep sequencing and bioinformatics, we identified a multi-step adaptive process to adaptation involving residues in the receptor footprints that correlated with receptor availability and with increase in virus fitness in an environment-specific manner. We show that adaptation occurs by selection of a dominant mutation followed by group selection of minority variants that together, confer the fitness increase observed in the population, rather than selection of a single dominant genotype.

Mutations in coronavirus nonstructural protein 10 decrease virus replication fidelity.

UNLABELLED: Coronaviruses (CoVs) are unique in encoding a 3'→5' exoribonuclease within nonstructural protein 14 (nsp14-ExoN) that is required for high-fidelity replication, likely via proofreading. nsp14 associates with the CoV RNA-dependent RNA polymerase (nsp12-RdRp), and nsp14-ExoN activity is enhanced by binding nsp10, a small nonenzymatic protein. However, it is not known whether nsp10 functions in the regulation of CoV replication fidelity. To test this, we engineered single and double alanine substitution mutations into the genome of murine hepatitis virus (MHV-A59) containing ExoN activity [ExoN(+)] at positions within nsp10 known to disrupt the nsp10-nsp14 interaction in vitro. We show that an nsp10 mutant, R80A/E82A-ExoN(+), was five to ten times more sensitive to treatment with the RNA mutagen 5-fluorouracil (5-FU) than wild-type (WT)-ExoN(+), suggestive of decreased replication fidelity. This decreased-fidelity phenotype was confirmed using two additional nucleoside analogs, 5-azacytidine and ribavirin. R80A/E82A-ExoN(+) reached a peak titer similar to and demonstrated RNA synthesis kinetics comparable to those seen with WT-ExoN(+). No change in 5-FU sensitivity was observed for R80A/E82A-ExoN(-) relative to MHV-ExoN(-), indicating that the decreased-fidelity phenotype of R80A/E82A-ExoN(-) is linked to the presence of ExoN activity. Our results demonstrate that nsp10 is important for CoV replication fidelity and support the hypothesis that nsp10 functions to regulate nsp14-ExoN activity during virus replication. IMPORTANCE: The adaptive capacity of CoVs, as well as all other RNA viruses, is partially attributed to the presence of extensive population genetic diversity. However, decreased fidelity is detrimental to CoV replication and virulence; mutant CoVs with decreased replication fidelity are attenuated and more sensitive to inhibition by RNA mutagens. Thus, identifying the viral protein determinants of CoV fidelity is important for understanding CoV replication, pathogenesis, and virulence. In this report, we show that nsp10, a small, nonenzymatic viral protein, contributes to CoV replication fidelity. Our data support the hypothesis that CoVs have evolved multiple proteins, in addition to nsp14-ExoN, that are responsible for maintaining the integrity of the largest known RNA genomes.

Viral Polymerase-Helicase Complexes Regulate Replication Fidelity To Overcome Intracellular Nucleotide Depletion.

UNLABELLED: To date, the majority of work on RNA virus replication fidelity has focused on the viral RNA polymerase, while the potential role of other viral replicase proteins in this process is poorly understood. Previous studies used resistance to broad-spectrum RNA mutagens, such as ribavirin, to identify polymerases with increased fidelity that avoid misincorporation of such base analogues. We identified a novel variant in the alphavirus viral helicase/protease, nonstructural protein 2 (nsP2) that operates in concert with the viral polymerase nsP4 to further alter replication complex fidelity, a functional linkage that was conserved among the alphavirus genus. Purified chikungunya virus nsP2 presented delayed helicase activity of the high-fidelity enzyme, and yet purified replication complexes manifested stronger RNA polymerization kinetics. Because mutagenic nucleoside analogs such as ribavirin also affect intracellular nucleotide pools, we addressed the link between nucleotide depletion and replication fidelity by using purine and pyrimidine biosynthesis inhibitors. High-fidelity viruses were more resistant to these conditions, and viral growth could be rescued by the addition of exogenous nucleosides, suggesting that mutagenesis by base analogues requires nucleotide pool depletion. This study describes a novel function for nsP2, highlighting the role of other components of the replication complex in regulating viral replication fidelity, and suggests that viruses can alter their replication complex fidelity to overcome intracellular nucleotide-depleting conditions. IMPORTANCE: Previous studies using the RNA mutagen ribavirin to select for drug-resistant variants have highlighted the essential role of the viral RNA-dependent RNA polymerase in regulating replication fidelity. However, the role of other viral replicase components in replication fidelity has not been studied in detail. We identified here an RNA mutagen-resistant variant of the nsP2 helicase/protease that conferred increased fidelity and yet could not operate in the same manner as high-fidelity polymerases. We show that the alphavirus helicase is a key component of the fidelity-regulating machinery. Our data show that the RNA mutagenic activity of compounds such as ribavirin is coupled to and potentiated by nucleotide depletion and that RNA viruses can fine-tune their replication fidelity when faced with an intracellular environment depleted of nucleotides.

Deep sequencing analysis of viral infection and evolution allows rapid and detailed characterization of viral mutant spectrum.

MOTIVATION: The study of RNA virus populations is a challenging task. Each population of RNA virus is composed of a collection of different, yet related genomes often referred to as mutant spectra or quasispecies. Virologists using deep sequencing technologies face major obstacles when studying virus population dynamics, both experimentally and in natural settings due to the relatively high error rates of these technologies and the lack of high performance pipelines. In order to overcome these hurdles we developed a computational pipeline, termed ViVan (Viral Variance Analysis). ViVan is a complete pipeline facilitating the identification, characterization and comparison of sequence variance in deep sequenced virus populations. RESULTS: Applying ViVan on deep sequenced data obtained from samples that were previously characterized by more classical approaches, we uncovered novel and potentially crucial aspects of virus populations. With our experimental work, we illustrate how ViVan can be used for studies ranging from the more practical, detection of resistant mutations and effects of antiviral treatments, to the more theoretical temporal characterization of the population in evolutionary studies. AVAILABILITY AND IMPLEMENTATION: Freely available on the web at http://www.vivanbioinfo.org CONTACT: : nshomron@post.tau.ac.il SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Alphavirus mutator variants present host-specific defects and attenuation in mammalian and insect models.

Arboviruses cycle through both vertebrates and invertebrates, which requires them to adapt to disparate hosts while maintaining genetic integrity during genome replication. To study the genetic mechanisms and determinants of these processes, we use chikungunya virus (CHIKV), a re-emerging human pathogen transmitted by the Aedes mosquito. We previously isolated a high fidelity (or antimutator) polymerase variant, C483Y, which had decreased fitness in both mammalian and mosquito hosts, suggesting this residue may be a key molecular determinant. To further investigate effects of position 483 on RNA-dependent RNA-polymerase (RdRp) fidelity, we substituted every amino acid at this position. We isolated novel mutators with decreased replication fidelity and higher mutation frequencies, allowing us to examine the fitness of error-prone arbovirus variants. Although CHIKV mutators displayed no major replication defects in mammalian cell culture, they had reduced specific infectivity and were attenuated in vivo. Unexpectedly, mutator phenotypes were suppressed in mosquito cells and the variants exhibited significant defects in RNA synthesis. Consequently, these replication defects resulted in strong selection for reversion during infection of mosquitoes. Since residue 483 is conserved among alphaviruses, we examined the analogous mutations in Sindbis virus (SINV), which also reduced polymerase fidelity and generated replication defects in mosquito cells. However, replication defects were mosquito cell-specific and were not observed in Drosophila S2 cells, allowing us to evaluate the potential attenuation of mutators in insect models where pressure for reversion was absent. Indeed, the SINV mutator variant was attenuated in fruit flies. These findings confirm that residue 483 is a determinant regulating alphavirus polymerase fidelity and demonstrate proof of principle that arboviruses can be attenuated in mammalian and insect hosts by reducing fidelity.

Coronaviruses lacking exoribonuclease activity are susceptible to lethal mutagenesis: evidence for proofreading and potential therapeutics.

No therapeutics or vaccines currently exist for human coronaviruses (HCoVs). The Severe Acute Respiratory Syndrome-associated coronavirus (SARS-CoV) epidemic in 2002-2003, and the recent emergence of Middle East Respiratory Syndrome coronavirus (MERS-CoV) in April 2012, emphasize the high probability of future zoonotic HCoV emergence causing severe and lethal human disease. Additionally, the resistance of SARS-CoV to ribavirin (RBV) demonstrates the need to define new targets for inhibition of CoV replication. CoVs express a 3'-to-5' exoribonuclease in nonstructural protein 14 (nsp14-ExoN) that is required for high-fidelity replication and is conserved across the CoV family. All genetic and biochemical data support the hypothesis that nsp14-ExoN has an RNA proofreading function. Thus, we hypothesized that ExoN is responsible for CoV resistance to RNA mutagens. We demonstrate that while wild-type (ExoN+) CoVs were resistant to RBV and 5-fluorouracil (5-FU), CoVs lacking ExoN activity (ExoN-) were up to 300-fold more sensitive. While the primary antiviral activity of RBV against CoVs was not mutagenesis, ExoN- CoVs treated with 5-FU demonstrated both enhanced sensitivity during multi-cycle replication, as well as decreased specific infectivity, consistent with 5-FU functioning as a mutagen. Comparison of full-genome next-generation sequencing of 5-FU treated SARS-CoV populations revealed a 16-fold increase in the number of mutations within the ExoN- population as compared to ExoN+. Ninety percent of these mutations represented A:G and U:C transitions, consistent with 5-FU incorporation during RNA synthesis. Together our results constitute direct evidence that CoV ExoN activity provides a critical proofreading function during virus replication. Furthermore, these studies identify ExoN as the first viral protein distinct from the RdRp that determines the sensitivity of RNA viruses to mutagens. Finally, our results show the importance of ExoN as a target for inhibition, and suggest that small-molecule inhibitors of ExoN activity could be potential pan-CoV therapeutics in combination with RBV or RNA mutagens.

The DNA virus Invertebrate iridescent virus 6 is a target of the Drosophila RNAi machinery.

RNA viruses in insects are targets of an RNA interference (RNAi)-based antiviral immune response, in which viral replication intermediates or viral dsRNA genomes are processed by Dicer-2 (Dcr-2) into viral small interfering RNAs (vsiRNAs). Whether dsDNA virus infections are controlled by the RNAi pathway remains to be determined. Here, we analyzed the role of RNAi in DNA virus infection using Drosophila melanogaster infected with Invertebrate iridescent virus 6 (IIV-6) as a model. We show that Dcr-2 and Argonaute-2 mutant flies are more sensitive to virus infection, suggesting that vsiRNAs contribute to the control of DNA virus infection. Indeed, small RNA sequencing of IIV-6-infected WT and RNAi mutant flies identified abundant vsiRNAs that were produced in a Dcr-2-dependent manner. We observed a highly uneven distribution with strong clustering of vsiRNAs to small defined regions (hotspots) and modest coverage at other regions (coldspots). vsiRNAs mapped in similar proportions to both strands of the viral genome, suggesting that long dsRNA derived from convergent overlapping transcripts serves as a substrate for Dcr-2. In agreement, strand-specific RT-PCR and Northern blot analyses indicated that antisense transcripts are produced during infection. Moreover, we show that vsiRNAs are functional in silencing reporter constructs carrying fragments of the IIV-6 genome. Together, our data indicate that RNAi provides antiviral defense against dsDNA viruses in animals. Thus, RNAi is the predominant antiviral defense mechanism in insects that provides protection against all major classes of viruses.

Coxsackievirus B3 mutator strains are attenuated in vivo.

Based on structural data of the RNA-dependent RNA polymerase, rational targeting of key residues, and screens for Coxsackievirus B3 fidelity variants, we isolated nine polymerase variants with mutator phenotypes, which allowed us to probe the effects of lowering fidelity on virus replication, mutability, and in vivo fitness. These mutator strains generate higher mutation frequencies than WT virus and are more sensitive to mutagenic treatments, and their purified polymerases present lower-fidelity profiles in an in vitro incorporation assay. Whereas these strains replicate with WT-like kinetics in tissue culture, in vivo infections reveal a strong correlation between mutation frequency and fitness. Variants with the highest mutation frequencies are less fit in vivo and fail to productively infect important target organs, such as the heart or pancreas. Furthermore, whereas WT virus is readily detectable in target organs 30 d after infection, some variants fail to successfully establish persistent infections. Our results show that, although mutator strains are sufficiently fit when grown in large population size, their fitness is greatly impacted when subjected to severe bottlenecking, which would occur during in vivo infection. The data indicate that, although RNA viruses have extreme mutation frequencies to maximize adaptability, nature has fine-tuned replication fidelity. Our work forges ground in showing that the mutability of RNA viruses does have an upper limit, where larger than natural genetic diversity is deleterious to virus survival.

Mutations in GLIS3 are responsible for a rare syndrome with neonatal diabetes mellitus and congenital hypothyroidism.

We recently described a new neonatal diabetes syndrome associated with congenital hypothyroidism, congenital glaucoma, hepatic fibrosis and polycystic kidneys. Here, we show that this syndrome results from mutations in GLIS3, encoding GLI similar 3, a recently identified transcription factor. In the original family, we identified a frameshift mutation predicted to result in a truncated protein. In two other families with an incomplete syndrome, we found that affected individuals harbor deletions affecting the 11 or 12 5'-most exons of the gene. The absence of a major transcript in the pancreas and thyroid (deletions from both families) and an eye-specific transcript (deletion from one family), together with residual expression of some GLIS3 transcripts, seems to explain the incomplete clinical manifestations in these individuals. GLIS3 is expressed in the pancreas from early developmental stages, with greater expression in beta cells than in other pancreatic tissues. These results demonstrate a major role for GLIS3 in the development of pancreatic beta cells and the thyroid, eye, liver and kidney.