key: cord-0963703-qcbskifq
authors: Smith, Everett Clinton; Blanc, Hervé; Vignuzzi, Marco; Denison, Mark R.
title: Coronaviruses Lacking Exoribonuclease Activity Are Susceptible to Lethal Mutagenesis: Evidence for Proofreading and Potential Therapeutics
date: 2013-08-15
journal: PLoS Pathog
DOI: 10.1371/journal.ppat.1003565
sha: 68422d2849ca1b297e561f73160f9a74a7077797
doc_id: 963703
cord_uid: qcbskifq

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 potential for CoVs to cause significant human disease is well demonstrated, with six known HCoVs-HKU1, OC43, NL63, 229E, SARS-CoV and MERS-CoV-causing colds, pneumonia, systemic infection, and severe or lethal disease [1] [2] [3] [4] [5] . Four of these viruses have been identified in just the last 10 years, with two, SARS-CoV and MERS-CoV, causing lethal respiratory and systemic infection [1, [3] [4] [5] [6] . Studies over the past 10 years have expanded the known phylogenetic, geographic, and species diversity of CoVs, and support multiple emergence events of CoVs into humans from bats and other zoonotic pools [7] [8] [9] [10] . The most recent evidence for CoV trans-species movement comes from the emergence of the novel MERS-CoV [1, 11, 12] . From April 2012 to June 2013 MERS-CoV has caused 72 laboratory confirmed cases and up to 50% mortality from severe respiratory and systemic disease in at least 8 countries, with evidence for human-to-human transmission [13] . MERS-CoV is most closely related to the bat CoVs HKU4 and HKU5 [11] , and the recently identified receptor dipeptidyl peptidase 4 (DPP4) is present on both human and bat cells [14] , providing a compelling argument that zoonotic CoV infections resulting in severe human disease may be more frequent events than previously thought. Because of the lack of epidemiological data, it remains unknown whether multiple introductions from a zoonotic source or human transmission of a mild or asymptomatic disease is responsible for these continuing cases of sporadic severe infections. However, based on the high mortality rates associated with SARS-CoV and those reported for MERS-CoV [13] , this novel virus potentially represents a serious threat to global health for which no vaccines or therapeutics currently exist.

CoVs contain the largest known RNA genomes (27-32 kb) and encode an array of 16 viral replicase proteins, including a 39-to-59 exoribonuclease (ExoN) domain within nonstructural protein 14 (nsp14) [2, [15] [16] [17] . Similar to the proofreading subunit (e) of E. coli DNA polymerase III, CoV nsp14-ExoN is a member of the DEDD superfamily of DNA and RNA exonucleases [15, 18] . This superfamily contains four conserved D-E-D-D acidic residues that are required for enzymatic activity, and mutation of these critical residues within CoV ExoN ablates or significantly reduces ExoN activity [15] . Studies from our group have demonstrated that ExoN activity is essential for high-fidelity replication in both the model CoV murine hepatitis virus (MHV) and SARS-CoV [19, 20] . Inactivation of ExoN activity due to alanine substitution of the first two active site residues results in 15-to 20-fold reduced replication fidelity in cell culture [19, 20] and a 12-fold reduction during SARS-CoV infection in vivo [21] , associated with profound and stable attenuation of SARS-CoV virulence and replication. A recent study has shown that bacterially-expressed SARS-CoV nsp14-ExoN can remove mismatched nucleotides in vitro, and that ExoN activity is stimulated in vitro through interactions with the non-enzymatic CoV protein nsp10 [22] . Thus all bioinformatic, genetic and biochemical studies to date support the hypothesis that nsp14-ExoN is the first identified proofreading enzyme for an RNA virus and functions together with other CoV replicase proteins to perform the crucial role of maintaining CoV replication fidelity.

Retrospective clinical studies during the SARS epidemic ultimately concluded that treatment with ribavirin (RBV), an antiviral drug shown to be mutagenic for some RNA viruses [23, 24] , was ineffective against SARS-CoV [25] [26] [27] [28] . Because ExoN activity is required for CoV high-fidelity replication [19] [20] [21] , we sought to determine if ExoN was responsible for CoV resistance to RNA mutagens. Using the nucleoside analog RBV and the base analog 5-fluorouracil (5-FU; [29] ) we show that CoVs lacking ExoN activity (ExoN2) are up to 300-fold more sensitive to inhibition than wild-type CoVs (ExoN+). Additionally, using full-genome nextgeneration sequencing we show that ExoN2 viruses accumulate 15to 20-fold more A:G and U:C transitions, consistent with 5-FU incorporation during RNA synthesis. Ultimately our results suggest the exciting possibility that small-molecule inhibitors of ExoN activity could be potential pan-CoV therapeutics, especially when used in combination with RBV or RNA mutagens.

Murine astrocytoma delayed brain tumor cells (DBT cells) were grown at 37uC and maintained in DMEM (Invitrogen) containing 10% FBS, supplemented with penicillin, streptomycin, HEPES, and amphotericin B. VeroE6 (Vero) cells were grown at 37uC and maintained in MEM (Invitrogen) containing 10% FBS supplemented with penicillin, streptomycin, and amphotericin B. All work with MHV was performed using the reverse genetics infectious clone based on strain MHV-A59 [30] , and work with SARS-CoV was performed using the reverse genetics infectious clone based on the Urbani strain [31] . Viral studies using SARS-CoV were performed in Select Agent certified BSL-3 laboratories using protocols reviewed and approved by the Institutional Biosafety Committee of Vanderbilt University and the Centers for Disease Control for the safe study and maintenance of SARS-CoV.

Compounds and cell viability studies 5-fluorouracil (5-FU), ribavirin (RBV), guanosine (GUA) and mycophenolic acid (MPA) were obtained from Sigma. 5-FU and RBV were made as 200 mM stock solutions, and were prepared in DMSO and sterile water, respectively. GUA and MPA were prepared in DMSO as 40 mM or 100 mM stocks, respectively. Low concentration (mM) working stocks were prepared as needed in sterile water prior to dilution in DMEM. Viability of DBT and Vero cells was assessed using CellTiter-Glo (Promega) in 96-well plate format according to manufacturer's instructions. DBT and Vero cells were seeded into opaque tissue culture grade 96-well plates, and DMEM containing RBV or 5-FU was added to each well to achieve the concentrations indicated. Water or DMSO vehicle controls were performed, in addition to a 20% ethanol control for cell death. The cells were then incubated at 37uC for either 12 or 24 h, and cell viability was determined using a Veritas Microplate Luminometer (Promega). The resultant values were then normalized to untreated cells.

Drug sensitivity studies and plaque assays Subconfluent monolayers of DBT cells in 6-well plates were pretreated for 30 min at 37uC with 1 mL of DMEM containing vehicle or the indicated concentration of RBV, 5-FU, MPA, or GUA. The drug was then removed and cells were infected with MHV-ExoN+ or ExoN2 viruses at an MOI of 1 plaque forming units (PFU)/cell (single-cycle) or 0.01 (multi-cycle) for 30 min at 37uC. Virus was then removed and 1 mL of DMEM containing vehicle, RBV, 5-FU, MPA, or GUA was added to each well. Cells were then incubated at 37uC for either 12 (single-cycle) or 24 (multi-cycle) h. The supernatant was harvested and virus titer was determined by plaque assay on DBT cells. For SARS-CoV studies, subconfluent monolayers of Vero cells in T25 flasks were pretreated for 30 min at 37uC with DMEM containing vehicle, RBV, or 5-FU. The drug was removed and cells were infected with either SARS-ExoN+ or ExoN2 viruses at an MOI of 0.1 PFU/cell (single-cycle) for 30 min. The virus was removed and DMEM containing vehicle, RBV, or 5-FU was added back. Cells were then incubated for 24 h, at which point the supernatant was harvested and virus titer was determined by plaque assay on Vero cells. All treated samples were normalized to the untreated vehicle control, and values were expressed as fold change from untreated virus titers.

RNA viruses have high mutation rates (10 23 to 10 25 mutations/nucleotide/round of replication), allowing for rapid viral adaptation in response to selective pressure. While RNA viruses have long been considered unable to correct mistakes during replication, CoVs such as SARS-CoV and the recently emerged MERS-CoV are important exceptions to this paradigm. All CoVs encode an exoribonuclease activity in nonstructural protein 14 (nsp14-ExoN) that is proposed to prevent and/or remove misincorporated nucleotides. Because of the demonstrated resistance of SARS-CoV to the antiviral drug ribavirin (RBV), we hypothesized that ExoN is responsible for CoV resistance to RNA mutagens. Using RBV and the RNA mutagen 5fluorouracil (5-FU), we show that CoVs lacking ExoN activity (ExoN2) are highly susceptible to RBV and 5-FU, in contrast to wild-type (ExoN+) CoVs. The inhibitory activity of 5-FU against ExoN2 viruses resulted specifically from 5-FU incorporation during viral RNA synthesis that lead to extensive mutagenesis within the viral population, and was associated with a profound decrease in virus specific infectivity. These results demonstrate the proofreading activity of ExoN during virus replication and suggest that inhibitors of ExoN activity could be broadly useful inhibitors of CoV replication in combination with RBV or RNA mutagens.

Real-time quantitative reverse transcription PCR (realtime qRT-PCR) of viral genomic RNA Viral RNA was harvested from infected cell monolayers using TRIzol reagent (Invitrogen), and was reverse transcribed (RT) using SuperScript III (Invitrogen). Random hexamers (1 mL of 50 mM stock) and 1 mg of total RNA were incubated for 5 min at 70uC. The remaining reagents were then added according to the manufacturer's protocol, and the mixture was incubated at 50uC for 1 h and then at 85uC for 5 min. All RT reactions were performed in a final volume of 20 mL. Real-time qRT-PCR was performed on the RT product using the Applied Biosciences 7500 Real-Time PCR System with Power SYBR Green PCR Master Mix (Life Technologies). Each reaction was performed in a total volume of 25 mL containing 12.5 mL of the Power SYBR Green PCR Master Mix, 125 ng each of the forward and reverse primers and 1 mL of the RT product which was diluted 1:1000. Viral genomic RNA was detected using primers (forward: ACAGGGTGGAGTTCCCGTTA and reverse: ACGGAAG-CACCACCATAAGA) optimized to generate a ,120 nt portion of ORF1a. These values were normalized using the 2 2DDCt method [32] to endogenous expression of the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) using primers (forward: GGGTGTGAACCACGAGAAAT and reverse: CCTTCCACAATGCCAAAGTT) optimized to yield a ,120 nt portion of GAPDH [33, 34] . Triplicate wells of each sample were analyzed, and averaged into one value representing a single replicate to minimize well-to-well variation. The cycle parameters were as follows: Stage 1, (1 rep) at 50uC for 2 min; Stage 2, (1 rep) 95uC for 10 min; Stage 3, (40 reps) at 95uC for 15 sec and 57uC for 1 min. One representative product from each treatment was verified by melting curve analysis and agarose gel electrophoresis.

Viral RNA from SARS-ExoN+ or ExoN2 infected Vero monolayers was harvested using TRIzol reagent, and was reverse transcribed (RT) using SuperScript III as described above except with 5 mL of random hexamers (50 mM stock), 5 mg of total RNA, and in a final volume of 100 mL for each reaction. Four microliters of RT product was then used to generate 12 overlapping ,3 kb amplicons for each virus treated with either 0 or 400 mM 5-FU by PCR. The high-fidelity polymerase Easy A (Agilent) was used to ensure that errors were minimal during PCR. All primer sets generated single bands which were then purified using the Wizard SV Gel and PCR Clean-Up System (Promega).

Prior to sequencing, cDNA amplicons were fragmented (Fragmentase, NEB), clustered, and sequenced with Illumina cBot and GAIIX technology as previously described [35] . Between 1.4610 8 and 4.5610 8 bases, comprised of ,69-nt reads, were obtained per virus, and CASAVA 1.8.2 was used to demultiplex and create the fastq files. Low quality bases from the ends of each sequence read were then trimmed, using Phred scores as the guiding metric (error probabilities higher than 0.001), and sequences with less than 16 bases after trimming were discarded to reduce false alignment and subsequent false variant calls. The program fastq-clipper (http://hannonlab.cshl.edu/fastx_toolkit/ index. html) was used for this quality filtering. The Burrows-Wheeler Alignment tool was then used to align reads to the SARS-CoV ExoN+ or ExoN2 reference genomes with a maximum of two mismatches per read [36] . Base calling at each position was determined using SAMTOOLS [37] . After the pileup, an in-house script collected the data per-position. For each position throughout the viral genome, the bases and their qualities were gathered, each variant allele's rate was initially modified according to its covering read qualities based on a maximum likelihood estimation and test for significance using Wilks' theorem. Additionally, an allele confidence interval was calculated and output for each allele. Only alleles with statistically significant p,0.05 values were retained and considered to be true variants. Above 0.01% all variants were found to be statistically significant, while below 0.01% many variants could not be distinguished from background error. Thus, the background noise caused by sequencing error was determined to be 0.01% or less.

Statistical tests were applied where noted within the figure legends and were determined using GraphPad Prism (La Jolla, CA) software. Statistical significance is denoted (*P,0.05, **P,0.01, ***P,0.0001) and was determined using an unpaired, two-tailed Student's t test compared to either untreated samples or to the corresponding ExoN+ sample. For the cell viability studies, treated samples were compared to the DMEM sample containing DMSO.

Because RBV has been shown to be incorporated as ribavirin monophosphate (RMP) into viral RNA during replication [23, 24, [38] [39] [40] [41] [42] , the presence of a proofreading enzyme would be predicted to exclude and/or remove nucleotide misincorporation [43] [44] [45] [46] [47] . If ExoN is responsible for the resistance phenotype, viruses lacking ExoN activity (ExoN2) should demonstrate increased titer reduction following RBV treatment as compared to wild-type viruses containing ExoN activity (ExoN+). To test this hypothesis, we examined the sensitivity of MHV-ExoN+ and ExoN2 viruses to RBV during single-cycle (MOI = 1 PFU/cell) replication in murine astrocytoma delayed brain tumor cells (DBT cells). No toxicity was observed in DBT cells following treatment with up to 400 mM RBV ( Figure 1A ). MHV-ExoN+ viruses were resistant to 10 mM RBV ( Figure 1B) , while MHV-ExoN2 virus titers decreased by ,200-fold following treatment with 10 mM RBV. The capacity of 10 mM RBV to inhibit MHV-ExoN2 replication is surprising because at least 10-fold higher concentrations of RBV are required to inhibit poliovirus and chikungunya viruses [48] [49] [50] . This observation could be due to the longer genomes of CoVs or to the mechanism(s) by which RBV inhibits CoV replication.

If RBV is exerting antiviral activity primarily through mutagenesis following incorporation of RMP, MHV-ExoN2 viruses should exhibit increased sensitivity during multi-cycle replication. To test this, we determined the sensitivity of MHV-ExoN+ and ExoN2 viruses to RBV at a low multiplicity of infection (MOI = 0.01 PFU/cell). Unexpectedly, multi-cycle replication of MHV-ExoN2 viruses in the presence of RBV ( Figure 1B ) was indistinguishable from single-cycle replication.

RBV has been reported to exert antiviral activity through numerous mechanisms [38] including disruption of viral RNA synthesis and inhibition of the cellular enzyme inosine monophosphate dehydrogenase (IMPDH). To determine if RBV treatment was affecting CoV RNA synthesis, we performed two-step realtime quantitative reverse transcription PCR (real-time qRT-PCR) MHV-ExoN2 (green) viruses following treatment with RBV is shown. DBT cells were infected with MHV-ExoN+ or MHV-ExoN2 in the presence or absence of RBV, and virus titer was determined by plaque assay. Genomic RNA levels were determined using two-step real-time qRT-PCR and primers to determine viral genomic RNA levels in the presence or absence of RBV. Similar to Figure 1B , MHV-ExoN+ titers were unaffected, whereas there was a dose-dependent reduction in MHV-ExoN2 titers following RBV treatment ( Figure 1C , filled bars). Corresponding dose-dependent reductions in MHV-ExoN2 genomic RNA were observed ( Figure 1C , hatched bars) following RBV treatment, demonstrating that treatment with 10 mM RBV decreased MHV-ExoN2 RNA synthesis by nearly 100-fold during replication. Because RBV caused decreased RNA synthesis in MHV-ExoN2 viruses, we calculated the relative specific infectivities of both viruses at each RBV concentration ( Table 1 ). The relative specific infectivity of MHV-ExoN2 viruses was decreased by 6-to 9-fold following treatment with RBV, while MHV-ExoN+ viruses were unaffected.

In addition to decreasing viral RNA synthesis, RBV could be exerting antiviral activity against MHV-ExoN2 through competitive inhibition of IMPDH by RMP [51] . To test this possible mechanism, we treated MHV-ExoN+ and MHV-ExoN2 viruses with the specific IMPDH inhibitor mycophenolic acid (MPA; [52] [53] [54] ) during both single-and multi-cycle replication. A concentration-dependent decrease in MHV-ExoN2 virus titer was observed following MPA treatment during single-cycle replication ( Figure 1D ). MHV-ExoN+ titers were reduced by less than 10-fold, consistent with what was observed following RBV treatment ( Figure 1B) . Similar to RBV, increased sensitivity of MHV-ExoN2 viruses to MPA was not observed during multicycle replication. If RBV is acting via IMDPH inhibition, addition of extracellular guanosine (GUA) should restore virus titers, as has been demonstrated previously for Dengue virus [55] . Addition of 100 mM GUA following RBV or MPA pretreatment and viral infection had no effect on MHV-ExoN+ viruses ( Figure 1E ), but completely restored MHV-ExoN2 titer even in the continued presence of 10 mM RBV or 1 mM MPA ( Figure 1F ). These data indicate that the antiviral activity of RBV against MHV-ExoN2 viruses is occurring at least in part through decreasing viral RNA synthesis and inhibition of IMPDH. Because our primary goal was to test the role of nsp14-ExoN in the prevention and/or removal of nucleotide misincorporation we did not further investigate how RBV was specifically inhibiting ExoN2 viruses. However, these results do show that the presence of ExoN activity is capable of preventing RBV inhibition of CoV replication.

The increased sensitivity of MHV-ExoN2 viruses to 5-FU treatment is consistent with mutagenesis We next examined the sensitivity of MHV-ExoN+ and ExoN2 viruses to the pyrimidine base analog 5-FU, which has been shown to be mutagenic for many RNA viruses [29, 56] . Treatment of DBT cells with up to 400 mM 5-FU did not result in any detectable cellular toxicity (Figure 2A ). Following treatment with up to 200 mM 5-FU ( Figure 2B ) during single-cycle infections, MHV-ExoN+ titers were inhibited less than 3-fold, while titers of MHV-ExoN2 decreased ,900 fold, representing a ,300-fold increase in sensitivity as compared to MHV-ExoN+. During multicycle replication, MHV-ExoN+ virus titers were reduced by less than 10-fold following 5-FU treatment, while MHV-ExoN2 showed a ,50,000-fold reduction in titer ( Figure 2B ). Virus was undetectable by plaque assay at 5-FU concentrations above 80 mM. Analysis of viral RNA synthesis by two-step real-time qRT-PCR demonstrated that MHV-ExoN+ RNA levels were not reduced following 5-FU treatment, while 5-FU treatment resulted in minimal two-to-five fold decreases in MHV-ExoN2 RNA ( Figure 2C ). The specific infectivity of MHV-ExoN2 was decreased by 14-and 128-fold following treatment with 100 mM and 200 mM 5-FU, respectively (Table 1) . These results demonstrate that ExoN activity confers resistance to 5-FU, and support the hypothesis that 5-FU is driving increased genomic mutagenesis in MHV-ExoN2 virus populations, leading to lethal mutagenesis and extinction.

To determine whether SARS-CoV viruses lacking ExoN activity (SARS-ExoN2) also were inhibited by RBV and 5-FU, we infected Vero cells with either SARS-ExoN+ or ExoN2 viruses in the presence or absence of RBV or 5-FU. Treatment of Vero cells with up to 400 mM RBV or 5-FU did not decrease cell viability by more than 20% ( Figure 3A ). Recent reports have optimized to amplify a ,120 nt region of ORF1a [33] . The change in genomic RNA levels (2 2DDCt ) is shown relative to endogenous GAPDH expression and was normalized to RNA levels from described the lack of RBV uptake by Vero cells due to the absence of specific equilibrative nucleoside transporters [57, 58] . Additionally, previous studies have shown that RBV failed to inhibit SARS-CoV replication in Vero cells [59] . Consistent with those reports, in our experiments both SARS-ExoN+ and ExoN2 viruses were unaffected by treatment with up to 400 mM RBV ( Figure 3B ). We therefore performed subsequent experiments with 5-FU. SARS-ExoN+ titers were reduced 3-and 10-fold following treatment with 200 or 400 mM 5-FU, respectively ( Figure 3C ). In contrast, SARS-ExoN2 titers were reduced ,300-fold by 200 mM 5-FU ( Figure 3C ), similar to MHV-ExoN2 viruses. At 400 mM 5-FU, SARS-ExoN2 virus was inhibited 2,000-fold during a single replication cycle, representing a ,160-fold increase in 5-FU sensitivity compared to SARS-ExoN+ viruses. Thus, our data indicate that increased sensitivity of CoVs to RNA mutagens in the absence of ExoN activity is conserved across diverse members of the CoV family. Of interest, our studies with SARS-ExoN+ also indicate that ExoN-mediated protection from nucleotide misincorporation can be overcome at higher concentrations of mutagen.

Studies with the RNA viruses lymphocytic choriomeningitis virus (LCMV), foot-and-mouth disease virus (FMDV) and vesicular stomatitis virus (VSV) have demonstrated that 5-FU is incorporated as 5-fluorouridine monophosphate (FUMP) into replicating viral RNA, thus increasing genomic mutations [60] [61] [62] .

To determine whether 5-FU was causing increased mutagenesis in SARS-CoV populations, we performed full-genome NGS analysis of both virus populations replicating in the presence or absence of 5-FU. To analyze the entire spectrum of mutations arising during replication, we extracted total intracellular RNA from Vero cells infected with either SARS-ExoN+ or ExoN2 viruses following treatment with either 0 mM or 400 mM 5-FU. We then generated 12 overlapping cDNA amplicons of approximately 3 kb in length for each sample. For each of the four samples, 1.4610 8 to 4.5610 8 bases were sequenced, corresponding to an average coverage depth of between 4,600 and 15,000 at each nucleotide position. We compared the statistically significant minority variants, defined as having a p-value of #0.05 following a multiple-testing correction (Benjamini-Hochberg), between the untreated and 5-FU-treated SARS-ExoN+ and ExoN2 populations. Following treatment with 400 mM 5-FU ( Figure 3D) , there was an increase in mutations within the SARS-ExoN+ population from 11 to 259 (24-fold). In contrast, for SARS-ExoN2 there were 3648 mutations present within the 5-FU-treated SARS-ExoN2 population compared to the 99 mutations in the untreated population (40-fold increase). Most remarkably, this represented a 16-fold increase in the number of statistically significant minority variants between 5-FU treated ExoN+ and ExoN2 SARS-CoV. Thus, these data support our hypothesis that 5-FU was increasing genomic mutations through incorporation of FUMP into viral genomes in the absence of ExoN activity. The change in titer (filled bars) and genomic RNA levels (hatched bars) of MHV-ExoN+ (blue) and MHV-ExoN2 (green) viruses following treatment with 5-FU is shown. DBT cells were infected with MHV-ExoN+ or MHV-ExoN2 in the presence or absence of 5-FU, and virus titer was determined by plaque assay. Genomic RNA levels were determined using two-step real-time qRT-PCR and primers optimized to amplify a ,120 nt region of ORF1a [33] . The change in genomic RNA levels (2 2DDCt ) is shown relative to endogenous GAPDH expression and was normalized to RNA levels from untreated samples. Mean values 6 S.E.M. are shown, n = 6. For all parts, statistical significance was determined using an unpaired, two-tailed Student's t test (*P,0.05, **P,0.01, ***P,0.0001). doi:10.1371/journal.ppat.1003565.g002 5-FU-associated A-to-G and U-to-C transitions are highly represented and distributed across the genome Incorporation of FUMP instead of uracil into replicating RNA allows FUMP to base pair with both guanosine and adenine [61, 63] . This decreased specificity in base pairing has been shown in studies with LCMV and primarily results in A-to-G (A:G) and U-to-C (U:C) transitions [29, 61, 63] . To determine if FUMP was being incorporated at higher levels in the absence of ExoNmediated proofreading, we analyzed the numbers and types of transitions and transversions occurring in each virus population (Figure 4 ). Transitions are indicated in grey boxes and transversions in white boxes, with the number for each shown. Transversions comprised the majority of variants for both untreated ExoN2 and ExoN+ viruses. Treatment with 5-FU caused the number of U:C and A:G transitions to increase in both ExoN+ and ExoN2 populations, from 2 to 197 for SARS-ExoN+ and from 16 to 3304 for SARS-ExoN2 (Figures 4A and B) . This increase and bias toward U:C and A:G transitions is consistent with FUMP being incorporated into both minus-and plus-strand RNA [63] during both ExoN+ and ExoN2 replication; however the absolute numbers were dramatically increased (16-fold) during ExoN2 replication compared to ExoN+. In untreated cells, A:G and U:C transitions accounted for less than 25% of the total minority variants within each population ( Figure 4C ). Following 5-FU treatment, A:G and U:C transitions accounted for 70-95% of the total minority variants within each population.

To further examine the genomic distribution of these two transitions, we plotted the total number of A:G and U:C transitions occurring at a frequency of between 0.1% and 1% ( Figure 5 ). Approximately 75% and 90% of the total minority variants occurring at a frequency between 0.1 and 1% following 5-FU treatment were due to A:G or U:C transitions ( Figure 5 ), for the SARS-ExoN+ and ExoN2 populations, respectively. In both populations, these mutations were distributed across the entire genome following treatment with 400 mM 5-FU. Thus our data provide direct evidence indicating that 5-FU drives increased genomic mutations within SARS-CoV in the absence of ExoN proofreading activity. For panels A-C statistical significance was determined using an unpaired, two-tailed Student's t test (*P,0.05, **P,0.01, ***P,0.0001). doi:10.1371/journal.ppat.1003565.g003

Discussion Viral sensitivity to RNA mutagens is determined by several factors including polymerase selectivity [39, 40, [64] [65] [66] [67] , mutational robustness [68] , and the acquisition of mutations that increase or decrease replication fidelity. Increased and decreased fidelity mutants have been described for picornaviruses and arboviruses [35, 48, 50, 69] , all of which have occurred in the viral RdRp. The CoV nsp14-ExoN is the first identified RNA virus protein distinct from the RdRp that affects replication fidelity [19] [20] [21] 70] . While the G641D mutation within the chikungunya (CHIKV) nonstructural protein 2 (nsP2) has been implicated in CHIKV resistance to RBV, a direct role for this protein in fidelity regulation has not been described [48] . A Sindbis virus variant containing mutations within nsP1, a viral guanylyl-and methyltransferase [71] , has been shown to be resistant to both RBV and MPA [72] . However, this phenotype is related to viral RNA capping and not replication fidelity [72] . In this report, we identify CoV ExoN activity as a critical determinant of viral sensitivity to RNA mutagens. Using two phylogenetically distant b-CoVs we demonstrate that this phenotype is well conserved across CoVs. Clearly, there is a profound increase both in overall mutations and in specific 5-FUassociated mutations within the ExoN2 population as compared to the ExoN+ wild-type population. Furthermore, the vast majority of statistically significant mutations were distributed genome-wide at frequencies between 0.2 and 1%, providing strong evidence supporting ExoN-mediated proofreading during CoV replication. Of interest, our experiments also revealed that ExoNmediated prevention and/or removal of misincorporated nucleotides is not absolute, especially in the setting of higher concentrations of mutagen. ExoN+ SARS-CoV populations demonstrated 24-fold more mutations following 5-FU treatment, suggesting that ExoN proofreading can be overwhelmed by higher concentrations of mutagens and likely by other nucleoside or base analogs. This raises the further possibility that ExoN may be less efficient at recognizing or removing some types of nucleoside or base analogs than others, and that such approaches to virus inhibition might be viable, particularly in combination with inhibitors that target ExoN activity.

The antiviral nucleoside analog RBV is currently used to treat hepatitis C virus (HCV; [73] [74] [75] ), Lassa virus [76] and respiratory syncytial virus (RSV) infections [77, 78] . The potential clinical use of RBV for CoV infections is complicated by the multiple mechanisms of action that have been reported [38] , and by the potential for disease exacerbation, as reported during the SARS-CoV epidemic [25] [26] [27] [28] . Our data suggest that RBV primarily inhibits MHV-ExoN2 virus replication through decreasing viral RNA synthesis and inhibition of IMPDH ( Figure 1 ). Inhibition of IMPDH by RMP has been shown to decrease intracellular GTP pools [51] , thus altering the balance of nucleoside triphosphates (NTPs) within the cell. Decreased GTP levels could result in forced misincorporations due to NTP imbalances in the absence of ExoN activity [72] . However, the moderate 6-to 9-fold decreases in relative specific infectivity observed for MHV-ExoN2 following RBV treatment (Table 1) suggests that mutagenesis is not the primary mechanism by which RBV is exerting an antiviral effect. An additional possibility is that the antiviral activity of RBV against ExoN2 viruses is unrelated to the putative proofreading function of this enzyme. Both biochemical and cell culture studies have demonstrated that loss of ExoN activity leads to impaired RNA synthesis [15, 19, 20] . Furthermore, in addition to ExoN activity, nsp14 contains N7-methyltransferase (N7-MTase) activity, a critical step in RNA capping [79, 80] . A recent report has demonstrated that the ExoN and N7-MTase domains are structurally inseparable, and that residues within the ExoN domain are important for N7-MTase activity [81] . Thus, the increased sensitivity of MHV-ExoN2 to RBV could result from the impairment of undefined functions of ExoN during CoV replication, particularly during RNA synthesis. The parallel use of ExoN+ and ExoN2 viruses with RBV may allow us to define how RBV is exerting an antiviral effect against CoVs and the potentially novel mechanisms by which ExoN may act to counter that inhibition.

Since the identification of nsp14-ExoN activity [15] and studies demonstrating the requirement for ExoN in high-fidelity replication [19] [20] [21] , mounting evidence points to a role for nsp14-ExoN in proofreading activity during RNA virus replication [22] . Here we used NGS to determine the number of mutations present in SARS-ExoN+ and ExoN2 populations. The characteristic 5-FUmediated transitions U:C and A:G comprised 90% of the total statistically significant minority variants within SARS-ExoN2 population, and were present at levels 15-and 20-fold higher than those same transitions within the ExoN+ population ( Figure 4) . Overall, our data represent the first direct test of ExoN proofreading during SARS-CoV replication in the absence of ExoN. Furthermore, the sequencing depth attained using NGS shows that ExoN inactivation likely skews the spectrum of spontaneous mutations present within the untreated population ( Figure 4 ). Such overrepresentation of specific mutations in the context of ExoN inactivation is similar to studies of S. cerevisiae DNA polymerases e and d containing mutations within their respective 39-to-59 DEDD exonucleases [82] [83] [84] [85] [86] . This altered distribution due to ExoN inactivation could have profound implications for CoV adaptation and evolution.

Lethal mutagenesis occurs through the accumulation of mutations within the viral genome during replication, and ultimately results in virus extinction (reviewed in [56, 87] ). While lethal mutagenesis has been studied extensively [87] , our work is the first to identify an RNA virus protein distinct from the RdRp that directly regulates the sensitivity of RNA viruses to genomic mutations resulting from mutagen incorporation. Currently, RBV is the only FDA-approved antiviral with demonstrated mutagenic activity. The first demonstration of RBV acting as a mutagen was performed using poliovirus [23, 24] almost 30 years after the antiviral activity of RBV was described [88] . The nucleoside analog T-705 (Favipiravir; [89] ) is currently in clinical development, and has been shown recently to drive lethal mutagenesis of influenza virus [90] . We have shown that ExoN+ viruses replicate well in the presence of RBV or 5-FU. However, we also have shown that ExoN2 mutants of SARS-CoV and MHV have 15-to-20-fold decreased fidelity [19, 20] , are attenuated, are subject to rapid loss of replication and clearance in vivo [21] , and are highly susceptible to low concentrations of RNA mutagens. An exciting possibility is that this conserved CoV proofreading enzyme could be targeted for inhibition, thus leading to the development of broadly useful CoV therapeutics. While ExoN inhibitors alone might be efficacious, combining an inhibitor of CoV fidelity with an RNA mutagen would magnify the intrinsic fidelity defect of ExoN inhibition and drive high-level mutagenesis. A potential advantage of such an approach would be to rapidly drive the virus to extinction, while limiting or blocking the capacity of the virus to overcome inhibition by reversion. ExoN2 mutants of both MHV and SARS-CoV have shown no reversion over multiple passages in culture or during persistent infections in vivo [19] [20] [21] . Furthermore, we did not observe any primary reversions within the ExoN DEDD motif following 5-FU treatment. While mutations within the CoV RdRp could emerge during acute treatment, mutations within other RNA virus RdRps have demonstrated that the maximum tolerance for increased or decreased fidelity without loss of virus viability is between ,3to 6-fold [35, 48, 69, 91] . In addition, our data demonstrate that ExoN2 viruses are profoundly sensitive to inhibition by lower concentrations of mutagen, providing a possible improved therapeutic index and margin of safety for use.

In summary, this study provides the most direct evidence to date that CoV ExoN provides a proofreading function during virus replication, and identifies ExoN as the critical determinant of CoV sensitivity to RNA mutagens. Because CoV replication fidelity is likely determined by the concerted effort of multiple virus proteins [19, 20, 22] , our data suggest the exciting possibility that significant attenuation of CoV fitness and pathogenesis could be achieved by targeting the conserved process of CoV replication fidelity. Ultimately, uncovering the mechanism of fidelity regulation and methodologies to disrupt this critical process will be vital to responding to both endemic and future emerging CoVs such as SARS-CoV and MERS-CoV.

Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia

Coronaviruses post-SARS: update on replication and pathogenesis

Identification of a novel coronavirus in patients with severe acute respiratory syndrome

A novel coronavirus associated with severe acute respiratory syndrome

Coronavirus as a possible cause of severe acute respiratory syndrome

Novel coronavirus -Saudi Arabia: human isolate

Evolutionary insights into the ecology of coronaviruses

Distant relatives of severe acute respiratory syndrome coronavirus and close relatives of human coronavirus 229E in bats

Detection and prevalence patterns of group I coronaviruses in bats

Evidence supporting a zoonotic origin of human coronavirus strain NL63

Genomic characterization of a newly discovered coronavirus associated with acute respiratory distress syndrome in humans

Severe respiratory illness caused by a novel coronavirus

Dipeptidyl peptidase 4 is a functional receptor for the emerging human coronavirus-EMC

Discovery of an RNA virus 39R59 exoribonuclease that is critically involved in coronavirus RNA synthesis

Nidovirales: evolving the largest RNA virus genome

Implications of altered replication fidelity on the evolution and pathogenesis of coronaviruses

Exoribonuclease superfamilies: structural analysis and phylogenetic distribution

Infidelity of SARS-CoV Nsp14-exonuclease mutant virus replication is revealed by complete genome sequencing

High fidelity of murine hepatitis virus replication is decreased in nsp14 exoribonuclease mutants

A live, impaired-fidelity coronavirus vaccine protects in an aged, immunocompromised mouse model of lethal disease

RNA 39-end mismatch excision by the severe acute respiratory syndrome coronavirus nonstructural protein nsp10/nsp14 exoribonuclease complex

RNA virus error catastrophe: direct molecular test by using ribavirin

The broadspectrum antiviral ribonucleoside ribavirin is an RNA virus mutagen

Adverse effects of ribavirin and outcome in severe acute respiratory syndrome: experience in two medical centers

Adverse events associated with high-dose ribavirin: evidence from the Toronto outbreak of severe acute respiratory syndrome

Enhancement of the infectivity of SARS-CoV in BALB/c mice by IMP dehydrogenase inhibitors, including ribavirin

SARS: systematic review of treatment effects

Molecular indetermination in the transition to error catastrophe: systematic elimination of lymphocytic choriomeningitis virus through mutagenesis does not correlate linearly with large increases in mutant spectrum complexity

Systematic assembly of a fulllength infectious cDNA of mouse hepatitis virus strain A59

Reverse genetics with a full-length infectious cDNA of severe acute respiratory syndrome coronavirus

Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method

Murine hepatitis virus replicase protein nsp10 is a critical regulator of viral RNA synthesis

Mutational analysis of MHV-A59 replicase protein-nsp10

Coxsackievirus B3 mutator strains are attenuated in vivo

Fast and accurate short read alignment with Burrows-Wheeler transform

The Sequence Alignment/Map format and SAMtools

Ribavirin's antiviral mechanism of action: lethal mutagenesis?

A multistep process of viral adaptation to a mutagenic nucleoside analogue by modulation of transition types leads to extinction-escape

Determinants of RNA-dependent RNA polymerase (in)fidelity revealed by kinetic analysis of the polymerase encoded by a foot-and-mouth disease virus mutant with reduced sensitivity to ribavirin

Hepatitis C virus RNAdependent RNA polymerase (NS5B) as a mediator of the antiviral activity of ribavirin

Foot-and-mouth disease virus mutant with decreased sensitivity to ribavirin: implications for error catastrophe

Inhibition of the 39R59 exonuclease of human DNA polymerase epsilon by fludarabine-terminated DNA

Differential incorporation and removal of antiviral deoxynucleotides by human DNA polymerase gamma

Toxicity of antiviral nucleoside analogs and the human mitochondrial DNA polymerase

Toxicity of nucleoside analogues used to treat AIDS and the selectivity of the mitochondrial DNA polymerase

Sensitivity of mitochondrial transcription and resistance of RNA polymerase II dependent nuclear transcription to antiviral ribonucleosides

Arbovirus high fidelity variant loses fitness in mosquitoes and mice

Engineering attenuated virus vaccines by controlling replication fidelity

A single mutation in poliovirus RNAdependent RNA polymerase confers resistance to mutagenic nucleotide analogs via increased fidelity

Mechanism of action of 1--D-ribofuranosyl-1,2,4-triazole-3-carboxamide (Virazole), a new broad-spectrum antiviral agent

Mycophenolic acid: an anti-cancer compound with unusual properties

Characterization of human type I and type II IMP dehydrogenases

The inhibition of nucleic acid synthesis by mycophenolic acid

Inhibition of dengue virus replication by mycophenolic acid and ribavirin

Viral quasispecies evolution

Reduced ribavirin antiviral efficacy via nucleoside transporter-mediated drug resistance

Cell type mediated resistance of vesicular stomatitis virus and Sendai virus to ribavirin

Glycyrrhizin, an active component of liquorice roots, and replication of SARSassociated coronavirus

Mutagenesis-mediated virus extinction: virus-dependent effect of viral load on sensitivity to lethal defection

Molecular characterization of a dual inhibitory and mutagenic activity of 5-fluorouridine triphosphate on viral RNA synthesis. Implications for lethal mutagenesis

Arenavirus extinction through lethal mutagenesis

5-fluorouracil in lethal mutagenesis of foot-and-mouth disease virus

Remote site control of an active site fidelity checkpoint in a viral RNA-dependent RNA polymerase

Structural basis for active site closure by the poliovirus RNA-dependent RNA polymerase

Stabilization of poliovirus polymerase by NTP binding and fingers-thumb interactions

A substitution in the fingers domain of DNA polymerase delta reduces fidelity by altering nucleotide discrimination in the catalytic site

Mutational robustness of an RNA virus influences sensitivity to lethal mutagenesis

Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population

Coronaviruses: an RNA proofreading machine regulates replication fidelity and diversity

Association of the Sindbis virus RNA methyltransferase activity with the nonstructural protein nsP1

Mutations that confer resistance to mycophenolic acid and ribavirin on Sindbis virus map to the nonstructural protein nsP1

The way forward in HCV treatment-finding the right path

Interferon and ribavirin vs interferon alone in the re-treatment of chronic hepatitis C previously nonresponsive to interferon: A meta-analysis of randomized trials

Interferon alfa-2b alone or in combination with ribavirin for the treatment of relapse of chronic hepatitis C. International Hepatitis Interventional Therapy Group

Lassa fever. Effective therapy with ribavirin

Ribavirin treatment of experimental respiratory syncytial viral infection. A controlled double-blind study in young adults

Respiratory syncytial virus (RSV) disease and prospects for its control

In vitro reconstitution of SARS-coronavirus mRNA cap methylation

Functional screen reveals SARS coronavirus nonstructural protein nsp14 as a novel cap N7 methyltransferase

Structure-function Analysis of SARS Coronavirus RNA Cap Guanine-N7 Methyltransferase

Saccharomyces cerevisiae DNA polymerase delta: high fidelity for base substitutions but lower fidelity for single-and multi-base deletions

Evidence from mutational specificity studies that yeast DNA polymerases delta and epsilon replicate different DNA strands at an intracellular replication fork

DNA polymerization in the absence of exonucleolytic proofreading: in vivo and in vitro studies

The 39R59 exonucleases of both DNA polymerases delta and epsilon participate in correcting errors of DNA replication in Saccharomyces cerevisiae

Unique error signature of the four-subunit yeast DNA polymerase epsilon

Lethal mutagenesis of viruses

Broad-spectrum antiviral activity of Virazole: 1-beta-D-ribofuranosyl-1,2,4-triazole-3-carboxamide

vitro and in vivo activities of anti-influenza virus compound T-705

T-705 (favipiravir) induces lethal mutagenesis in influenza A H1N1 viruses in vitro

Increased fidelity reduces poliovirus fitness and virulence under selective pressure in mice

We thank Ofer Isakov and Noam Shomron at Tel Aviv University for their help in NGS bioinformatics analysis and Michelle Becker at Vanderbilt University Medical Center for critical reading of the manuscript.