Next Article in Journal
MapReduce-Based Parallel Genetic Algorithm for CpG-Site Selection in Age Prediction
Next Article in Special Issue
Coding-Gene Coevolution Analysis of Rotavirus Proteins: A Bioinformatics and Statistical Approach
Previous Article in Journal
Novel TNXB Variants in Two Italian Patients with Classical-Like Ehlers-Danlos Syndrome
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 12
10.3390/genes10120968
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

RNA-Dependent RNA Polymerase Speed and Fidelity are not the Only Determinants of the Mechanism or Efficiency of Recombination

by
Hyejeong Kim
†,
Victor D. Ellis III
,
Andrew Woodman
,
Yan Zhao
,
Jamie J. Arnold
† and
Craig E. Cameron
*,†
Department of Biochemistry and Molecular Biology, The Pennsylvania State University, 201 Althouse Laboratory, University Park, PA 16802, USA
*
Author to whom correspondence should be addressed.
†
Present address: Department of Microbiology and Immunology, School of Medicine, University of North Carolina at Chapel Hill, 125 Mason Farm Rd., Chapel Hill, NC 27599-7290, USA.
Genes 2019, 10(12), 968; https://doi.org/10.3390/genes10120968
Submission received: 15 October 2019 / Revised: 14 November 2019 / Accepted: 21 November 2019 / Published: 25 November 2019
(This article belongs to the Special Issue Virus Evolution and Mutagenesis)
Browse Figures
Versions Notes

Abstract

:
Using the RNA-dependent RNA polymerase (RdRp) from poliovirus (PV) as our model system, we have shown that Lys-359 in motif-D functions as a general acid in the mechanism of nucleotidyl transfer. A K359H (KH) RdRp derivative is slow and faithful relative to wild-type enzyme. In the context of the KH virus, RdRp-coding sequence evolves, selecting for the following substitutions: I331F (IF, motif-C) and P356S (PS, motif-D). We have evaluated IF-KH, PS-KH, and IF-PS-KH viruses and enzymes. The speed and fidelity of each double mutant are equivalent. Each exhibits a unique recombination phenotype, with IF-KH being competent for copy-choice recombination and PS-KH being competent for forced-copy-choice recombination. Although the IF-PS-KH RdRp exhibits biochemical properties within twofold of wild type, the virus is impaired substantially for recombination in cells. We conclude that there are biochemical properties of the RdRp in addition to speed and fidelity that determine the mechanism and efficiency of recombination. The interwoven nature of speed, fidelity, the undefined property suggested here, and recombination makes it impossible to attribute a single property of the RdRp to fitness. However, the derivatives described here may permit elucidation of the importance of recombination on the fitness of the viral population in a background of constant polymerase speed and fidelity.

1. Introduction

In spite of the substantial resources that have been allocated by the World Health Organization, the National Institute of Allergy and Infectious Diseases, the European Centre for Disease Prevention and Control, and the Centers for Disease Control and Prevention to support prediction of emerging viral pathogens [1], viral outbreaks over the past few decades have been caused by viruses for which surveillance was not considered a priority. Rapid response to an outbreak caused by an unexpected viral pathogen requires, minimally, the existence of broad-spectrum, antiviral therapeutics. Prevention requires the availability of vaccines, development of which could take years. Indeed, approved vaccines still do not exist to prevent infections by West Nile virus or severe acute respiratory syndrome (SARS) coronavirus, and these outbreaks occurred more than one decade ago.
Because all RNA viruses encode an RNA-dependent RNA polymerase (RdRp) with conserved function and mechanism, this enzyme has emerged as an attractive target for development of broad-spectrum therapeutics [2,3,4] and a target for function/mechanism-based strategies for viral attenuation [5,6,7,8,9,10]. One function/mechanism of the RdRp that has been targeted most is that required for faithful incorporation of nucleotides [5,6,7,8,9,10]. Changing RdRp fidelity decreases or increases the genetic variation of the viral population, which, in turn, decreases fitness and virulence of the viral population [11,12,13,14,15,16].
Known RdRp variants exhibiting a high-fidelity phenotype can also exhibit a reduced speed of nucleotide incorporation, at least at the biochemical level [17]. Replication speed is also a determinant of viral fitness and virulence [18]. So, is it the increased fidelity or decreased speed of nucleotide addition that gives rise to the attenuated phenotype? Further complicating the fidelity-versus-speed question are the recent observations that changes to fidelity also have consequence for the efficiency of recombination [19,20,21,22,23,24]. Increased RdRp fidelity decreases recombination efficiency and vice versa [19,20,21,22,23,24]. It will likely be impossible to attribute a single biochemical property of the RdRp to biological outcome.
The most extensively studied PV fidelity mutants encode RdRps with amino acid substitutions at sites remote from the active site. Constructing equivalent substitutions conferring equivalent phenotypes in RdRps other than PV is difficult if not impossible. Our laboratory has therefore pursued active-site-based strategies to manipulate the fidelity, speed, and/or recombination efficiency of the RdRp [25]. We have shown that a lysine (Lys-359 in PV) present in conserved structural motif D of the RdRp contributes to the efficiency of nucleotidyl transfer [25]. A PV mutant encoding a K359R RdRp is attenuated but elicits a protective immune response in mice that is at least as robust as the immune response elicited by the type 1 Sabin vaccine strain [7].
In this study, we characterize a second motif-D mutant of PV (K359H). Unlike the K359R RdRp-encoding virus [7], K359H PV is genetically unstable and acquires mutations encoding two second-site amino acid substitutions after a few passages in cell culture. Together, the two substitutions restore all biochemical properties of the derivative to near wild-type levels. Individually, however, we observe differences in the mechanism (copy-choice vs. forced-copy-choice) and efficiency of recombination by each derivative, although each derivative exhibits equivalent speed and fidelity. We conclude that biochemical properties in addition to speed and fidelity must exist and contribute to both the mechanism and efficiency of recombination. The desire to attribute single, biochemical properties of the viral RdRp to fitness, virulence, and pathogenesis may be futile.

2. Materials and Methods

2.1. Cells and Viruses

Adherent monolayers of HeLa and L929 fibroblasts were grown in DMEM/F-12 media. Media was supplemented with 100 U/mL penicillin, 100 µg/mL streptomycin, and 10% heat inactivated (HI)-FBS. All cells were passaged in the presence of trypsin-EDTA. Wild-type and recombinant PV viruses were recovered after transfection of RNA produced in vitro (see below) from full-length cDNA or from the CRE-REP assay parental partners [23]. PV type 1 (Mahoney) was used throughout this study (GenBank accession number: V01149.1).

2.2. Plasmids, In Vitro Transcription, Cell Transfection, and Virus Quantification

All mutations were introduced into the pET26Ub-PV 3D [26] bacterial expression plasmids using overlap extension PCR. The presence of the desired mutations and the absence of additional mutations were verified by DNA sequencing. The expression plasmid encoding the various mutations was digested and inserted into an intermediate plasmid, pUC18-BglII-EcoRI-3CD (referred as pUC-3CD in [27], and then the fragment between BglII and ApaI was cloned into the subgenomic replicon, pRLucRA, or viral cDNA, pMovRA, or PV1ΔCRE. PV1ΔCRE is a full-length PV type 1 (Mahoney) cDNA bearing 8 synonymous substitutions in the oriI cis-acting replication element (CRE) located in the 2C-coding region and was described previously [23]. In order to create a fluorescence-based assay for recombination in cells, we modified the PV1ΔCRE acceptor RNA to include the coding sequence for UnaG green fluorescent protein. Briefly, UnaG-encoding sequence [28] carrying a 3C protease cleavage site at its carboxyl terminus was inserted between the IRES and the P1 region of the PV sequence. Translation occurred from the natural poliovirus initiation codon. Proteolytic cleavage and release of the fluorescent protein occurred by normal 3C protease activity. Plasmids encoding PV genomes (full length or subgenomic) were linearized with ApaI. All linearized cDNAs were transcribed in vitro using T7 RNA Polymerase and treated with 2U DNAse Turbo (Thermo Fisher, Waltham, MA, USA) to remove residual DNA template. The RNA transcripts were purified using RNeasy Mini Kit (Qiagen, Germantown, MD, USA) before spectrophotometric quantification. Purified RNA in RNase-free H2O was transfected into either HeLa or L929 fibroblasts using TransMessenger (Qiagen, Germantown, MD, USA). Virus yield was quantified by plaque assay. Briefly, cells and media were harvested at 2–3 days post-transfection, subjected to three freeze-thaw cycles, and clarified. Supernatant was then used on fresh HeLa cells in 6-well plates; virus infection was allowed to continue for 30 min. Media was then removed, and cells were washed twice with PBS (pH 7.4) washes before a 0.8% (w/v) agarose-media overlay was added. Cells were incubated for 2–3 days and then fixed and stained with crystal violet for virus quantification.

2.3. Luciferase Assays

Subgenomic luciferase assays were performed as described previously [29].

2.4. Virus Sequencing

Viral RNA was extracted from clarified culture supernatant using a Qiagen RNAeasy Mini kit, reverse transcribed using Superscript II reverse transcriptase (Invitrogen, Waltham, MA, USA) using an oligo-dT primer according to the manufacturer’s protocol. PCR amplification of the 3Dpol region in PV used template cDNA and appropriate oligonucleotides as listed in Table 1 by using Phusion high-fidelity DNA polymerase (NEB, Ipswich, MA, USA) according to the manufacturer’s protocol. PCR products were gel purified and sequenced by the Genomics Core Facility of the Pennsylvania State University.

2.5. One-Step Growth Analysis

HeLa cells in 12-well plates were infected by each virus at a MOI of 10. Following a 30-minute incubation, cells were washed twice with PBS and media was replaced. Virus was harvested at different time-points post infection and the virus yield was quantified by plaque assay.

2.6. Quantitative RT-PCR

Viral RNA was purified from virus stocks by using QiaAmp viral RNA purification kit (Qiagen, Germantown, MD, USA) and used for RT-qPCR to determine genome copies. This analysis was performed by the Genomics Core Facility of the Pennsylvania State University. DNAse- treated RNA was reverse-transcribed using the High Capacity cDNA Reverse Transcription kit (Applied Biosystems, Foster City, CA, USA) following the protocol provided with the kit. Quantification by real-time qPCR was done with 2X TaqMan Universal PCR Master Mix (Applied Biosystems, Foster City, CA, USA) in a volume of 20 μL, with primers 5′- ACCCCTGGTAGCAATCAATATCTTAC-3′ (forward) and 5′- TTCTTTACTTCACCGGGTATGTCA-3′ (reverse), and probe 5′-[6-Fam] TGTGCGCTGCCTGAATTTGATGTGA-3′ in a 7300 Real-Time qPCR System (Foster City, CA, USA) machine. A standard curve was generated using in vitro transcribed RNA.

2.7. Ribavirin-Sensitivity Assay

This assay was described previously [30]. HeLa cells were treated with 600 µM ribavirin for 1 h before infection. Ribavirin-treated cells were then infected at MOI 0.1 with each virus variant. Following infection, the cells were washed with PBS and media was replaced with ribavirin. Infection was allowed for 24 h. Cells and supernatant were subjected to three freeze–thaw cycles. Media was clarified and used for plaque assays. All yields were normalized to an untreated control.

2.8. Cell-Based Recombination Assay

This assay was developed in the laboratory of David Evans and was used as described previously [23,31]. L929 fibroblasts were transfected with the PV donor and accepter RNAs, both carrying the same mutations in the RdRp gene. For the fluorescence-based detection of recombinants, the UnaG green fluorescence reporter was encoded by the acceptor RNA. Supernatant at 2–3 days post-transfection was used to infect HeLa cells. Recombinant viruses were either quantified by plaque assay or by fluorescence imaging (488 nm (ex) and 509 nm (em)).

2.9. Polymerase Expression and Purification

Purified polymerases used for biochemical analysis were prepared as described previously [26,27].

2.10. Poly(rU) Polymerase Activity Assay

Reactions contained 50 mM HEPES, pH 7.5, 10 mM 2-mercaptoethanol, 5 mM MgCl2, 60 μM ZnCl2, dT15 (2 μM), poly(rA) (100 μM AMP), UTP (500 μM), α-32P-UTP (0.2 μCi/μL), and PV RdRp (0.2 μM). Reactions were initiated by addition of PV RdRp and incubated at 30 °C for 5 min at which time the reactions were quenched by addition of EDTA to a final concentration of 50 mM. Reaction volumes were 50 μL. Products were analyzed by DE81 filter binding, where 10 μL of the quenched reaction was spotted onto DE81 filter paper discs and dried completely. The discs were washed three times for 10 min in 250 mL of 5% dibasic sodium phosphate and rinsed in absolute ethanol. Bound radioactivity was quantitated by liquid scintillation counting in 5 mL of Ecoscint scintillation fluid (National Diagnostics, Atlanta, GA, USA).

2.11. Sym/Sub-Based Template Switching Assay

The sym/sub assay has been described previously [23,32]. Elongation complexes were assembled by incubating 5 μM PV RdRp with 1 μM sym/sub RNA primer-template and 500 μM ATP for 5 min (Mix 1). Template-switching reactions were initiated by addition of 60 μM RNA acceptor template and 500 μM CTP, GTP, and UTP (Mix 2) and then quenched at various times by addition of 50 mM EDTA. All reactions were performed at 30 °C in 50 mM HEPES, pH 7.5, 10 mM 2-mercaptoethanol, 60 μM ZnCl2, and 5 mM MgCl2. Products were analyzed by denaturing polyacrylamide gel electrophoresis, visualized using a PhosphorImager, and the transfer products quantified using ImageQuant TL software (GE Healthcare, Chicago, IL, USA).

3. Results

3.1. K359H PV Requires Two Second-Site Substitutions to Restore a “Wild-Type” Growth Phenotype

Most studies of RdRp fidelity have benefited from the selection of derivatives that were either more or less sensitive to a mutagenic nucleoside [33,34,35,36,37,38]. Almost invariably, the derivatives changed residues remote from the catalytic site, thus using an allosteric mechanism to perturb fidelity. Among the most famous of these is the G64S substitution in PV RdRp [16,17,34]. Many years ago, our laboratory showed that the active site of all polymerases contain a general acid, Lys-359 in the case of PV RdRp, that protonated the pyrophosphate leaving group during nucleotidyl transfer, thereby increasing the catalytic efficiency of RdRp [25]. Substitutions of Lys-359 in PV RdRp give rise to changes in fidelity [7,18]. The arginine substitution of Lys-359 (K359R RdRp) and the histidine substitution (K359H RdRp) catalyze nucleotidyl transfer at a rate 10-fold lower than wild type [25]. The impact of the arginine substitution of Lys-359 on PV fidelity and its potential application to vaccine development have been described [7].
Biological studies of K359H PV were not reported, because this virus was not genetically stable. When PV is rescued from in vitro transcribed RNA, four passages are required for the genetic diversity of the viral population to come to equilibrium [11]. During this time for K359H PV, mutations were observed in RdRp-coding sequence that changed the plaque phenotype of the virus. One change was in motif C, I331F, and the other change was in motif D, P356S (Figure 1A,B). We have adopted the following nomenclature to refer to the various mutant PVs and their corresponding RdRps: I331F-K359H, IF-KH; P356S-K359H, PS-KH; and I331F-P356S-K359H, IF-PS-KH. In order to compare the impact of these substitutions on PV genome replication without complication of reversion, we engineered the various substitutions into a subgenomic replicon, where the capsid coding P1 region was replaced with the gene for firefly luciferase, producing luciferase as an indirect measure of viral RNA produced [39,40,41]. Experiments in the presence of guanidine hydrochloride (GuHCl) report on translation of transfected RNA in the absence of replication (Figure 1C). Relative to WT PV, K359H was the most debilitated, with replication failing to reach an end point at 10 h post-transfection (Figure 1C). Each double mutant was markedly better than K359H alone, but the addition of I331F conferred a greater replication advantage than P356S (Figure 1C). The triple mutant replicated even faster but still exhibited a significant reduction in the rate of replication relative to WT (Figure 1C).

3.2. Characterization of the Double- and Triple-Mutant Viruses and Their Polymerases

At multiplicities of infection of one or higher, the double- and triple-mutant viruses were stable for multiple passages, thus permitting us to characterize the biological properties of these viruses. Interestingly, the growth properties of the viruses were not as expected based on the experiments with the replicon. Each mutant virus exhibited the same delay relative to WT virus prior to the first detection of infectious virus after four hours post-infection (hpi) (Figure 2A). Thereafter, growth of each virus was far more distinct than observed for the replicon, with PS-KH PV much slower than IF-KH PV and IF-PS-KH PV was substantially faster than both double mutants (Figure 2A). The differences in outcome could reflect the differences in replication efficiency and/or a direct consequence of the substitutions in 3D-coding sequence on virus assembly or spread caused by changes to 3CD protein or 3D-containing precursor protein [29,42]. After a productive infection, the exponential phase of genome replication by PV typically occurs between 2 and 6 h post-infection, with virus assembly occurring between 4 and 8 h post-infection. Release of infectious virus by lytic and/or non-lytic mechanisms can occur once virus is assembled and can persist for longer than 24 hours.
Our studies of RdRp fidelity mutants in cell culture have highlighted the fact that the use of the plaque-forming unit (PFU) as a measure of virus concentration will mask changes in the specific infectivity of the viral RNA [11]. As a result, our experiments generally use genomes instead of PFU as the unit of measure [11]. In doing so, we can use genomes/PFU as a measure of the specific infectivity. The instability of KH PV precluded rigorous, quantitative analysis of this virus, but at least 10-fold more KH PV than WT PV was required to observe a comparable number of plaques. Both substitutions were able to increase the efficiency of plaque formation relative to KH PV (Figure 2B). Indeed, the specific infectivity of IF-KH PV was equivalent to WT. The addition of PS to IF-KH reduced the efficiency of plaque formation; the specific infectivity of IF-PS-KH PV was reduced by twofold relative to IF-KH PV (Figure 2B). Therefore, the behavior of PS-KH PV in the PFU-based growth assay relative to the replicon assay likely reflects an additional defect to virus assembly and/or spread.
A change to the specific infectivity of the viral RNA caused by increased mutational load—that is, reduced fidelity of the PS-KH RdRp, would also explain the observed reduction in the specific infectivity of PS-KH PV relative to the other mutant PVs. To test this possibility, we evaluated the sensitivity of each virus population to growth in the presence of ribavirin. We reasoned that the mutagenic activity of ribavirin would exhibit the greatest negative impact on viruses with a mutator phenotype as described previously [11]. Relative to WT, K359H PV exhibited the highest fidelity based on the 2-log difference in sensitivity to ribavirin (Figure 2C). IF-KH and PS-KH PVs also exhibited a higher fidelity than WT PV (Figure 2C). Importantly, both mutants exhibited essentially equivalent fidelity phenotypes (Figure 2C), consistent with the suggestion above that the reduced efficiency of plaque formation observed for PS-KH PV is likely related to impairment of virus assembly and/or spread. Finally, evaluation of IF-PS-KH PV revealed equivalent sensitivity of this mutant to ribavirin as observed for WT PV (Figure 2C), suggesting that fidelity had been returned to normal.
We have previously reported biochemical properties of the RdRps for the repertoire of mutant viruses described above [18]. Consistent with these reports, each substitution increases the speed and reduces the fidelity to create a biochemical phenotype on par with that observed for WT (Figure 2D).

3.3. Biochemical Properties of the RdRp Other Than Speed and Fidelity Contribute to the Efficiency of Recombination in Cell Culture

Based on the myriad RdRp fidelity mutants that have been reported to date, there appears to be a direct correlation between the rate of nucleotide addition (speed) and the fidelity of nucleotide addition [7,11,18,23,43]. This correlation also extends to recombination efficiency [23]. If this is the case, then the ongoing debate of speed versus fidelity as the key determinant of viral fitness will become even more complicated to resolve [18].
The current state of the art for evaluation of recombination in cell culture is based on the co-transfection of two viral (sub) genomic RNAs incapable of producing infectious virus [31]. The donor RNA is a replication-competent, subgenomic RNA that encodes a luciferase reporter instead of the viral capsid (Figure 3A). The acceptor RNA is a replication-incompetent, genomic RNA that has a defective cis-acting replication element, termed oriI (Figure 3A) [31]. Initiation of replication on the donor followed by a switch to the acceptor at a site after the oriI locus will yield an infectious genome that can be scored by plaque assay (Figure 3A) [31].
We evaluated our panel of KH-containing PV mutants using this assay, and the outcomes were, in most cases, quite unexpected (Figure 3B). KH PV was unable to produce viable recombinants, as expected for a high-fidelity RdRp [23]. The first surprise was that IF-KH and PS-KH PVs did not exhibit the same phenotypes (Figure 3B). While IF-KH PV produced viable recombinant virus, PS-KH PV did not (Figure 3B). Based on the experiments performed above, these viruses and their polymerases replicate with comparable efficiency and fidelity (Figure 1A and Figure 2C,D). Even more surprising, however, was the observation that IF-PS-KH PV was impaired more than one log relative to WT PV in its ability to produce viable recombinants (Figure 3B). The replication efficiency and fidelity of this virus and its polymerase were always within twofold of that observed for WT (Figure 1A and Figure 2C,D). Together, these results suggest that there are properties of the PV RdRp important for recombination that are not revealed by our existing biological and biochemical assays.
The one caveat of the recombination assay is that the recombinant viruses produced must be able to form plaques, which may mean that the recombinant virus must spread by a lytic mechanism. It was possible that PS-KH PV was defective for virus assembly and/or spread, at least lytic spread (Figure 2B). In order to score for recombinant viruses that spread by either a lytic or non-lytic mechanism, we engineered the acceptor template to encode the UnaG [28] green fluorescent protein (Figure 3C). The polyprotein was designed such that UnaG protein is released from capsid precursor protein by 3C protease activity. Observation of green cells were visible from a donor-acceptor pair producing WT RdRp as early as 7 hpi, with increases continuing over a 48-h period (Figure 3D). With the sensitivity of this assay and the absence of a requirement for lytic spread, it is clear that PS-KH PV exhibits a recombination defect (Figure 3D).

3.4. Poly(rU) Polymerase Activity as a Predictor of the Efficiency of Copy-Choice Recombination in Cell Culture

Twenty years ago, we showed that template switching was the primary mechanism of product formation when oligo (dT) or oligo(rU) were used to prime poly(rU) RNA synthesis on oligo(rA) or poly(rA) RNA templates [32]. Briefly, polymerase engages a primed template, elongates that primer, disengages from the first template at some point during the elongation process to engage a second (acceptor) template and continue RNA synthesis (Figure 4A) [32]. This process occurs reiteratively, yielding products that are much greater than the average length of template used in the reaction [32].
The observation that only a subset of the KH mutants were competent for recombination, presumably by a copy-choice mechanism [44], provided the opportunity to determine the extent to which the poly (rU) polymerase activity can predict biological phenotypes. The congruence between the two experiments was remarkable (Figure 4B). Both KH and PS-KH RdRps were impaired for poly (rU) polymerase activity, both IF-KH and IF-PS-KH RdRps exhibited near-WT levels of poly (rU) polymerase activity (Figure 4B). These results further confirm template switching as the primary mechanism of poly (rU) polymerase activity and validate this assay as a screen for identification of RdRps with deficits in template switching.

3.5. Polymerase Determinants Supporting Efficient Copy-Choice Recombination do not Overlap Completely With Determinants Supporting Efficient Forced-Copy-Choice Recombination

Recently, our laboratory developed an assay for forced-copy-choice recombination [23,45], inspired by an analogous assay developed to study template switching by the reverse transcriptase from human immunodeficiency virus [46,47,48]. The concept of the assay and templates used are presented in Figure 5A. The experimental design is diagrammed in Figure 5B. Polymerase assembles on the primed template (sym/sub-U). In the presence of only the first nucleotide (ATP), the assembled complexes can be identified as the one-nucleotide-extended product (n + 1). Addition of the remaining nucleotides in the absence or presence of the acceptor template will yield a strong-stop RNA product. In the presence of a complementary acceptor RNA, the strong-strop RNA product (donor) will be extended, creating a transfer product.
We evaluated each KH-containing RdRp derivative in this assay. Substrate and product analysis are presented in Figure 5C, with the quantitation of the transfer product relative to WT presented in Figure 5D. KH RdRp was defective in this assay as well, as transfer product was not detected (Figure 5C). Importantly, the failure to transfer was not a reflection of the inability to assemble or produce strong-stop donor RNA (Figure 5C). In this assay, PS-KH RdRp outperformed IF-KH RdRp, and PS-KH and IF-PS-KH RdRp were essentially identical in activity.
We conclude that copy-choice and forced-copy-choice recombination use distinct mechanisms, requiring unique biochemical properties of the RdRp.

4. Discussion

For more than a decade now, our laboratory and others have published studies asserting connections between viral RdRp fidelity and viral fitness [11,12,13,14,15,16]. However, in the few instances in which RdRp derivatives exhibiting perturbed fidelity have been characterized, these derivatives also exhibit changes to the rate of nucleotide addition [7,11,18,23,43]. Relative to wild-type polymerase, a higher fidelity polymerase is a slower enzyme, and a lower fidelity polymerase is a faster enzyme [7,11,18,23,43]. Polymerase speed will clearly determine replication kinetics and therefore can contribute to viral fitness, as has been suggested recently [18]. The most extensively characterized fidelity mutants of PV harbor an RdRp with a substitution located at a remote site that likely causes substantial collateral damage, further confounding the fidelity-versus-speed debate [5,17]. This study was motivated by the need to better understand the relationship between the biochemical properties of the viral RdRp and viral fitness, virulence, and pathogenesis. We have had in hand for a long time a PV mutant whose speed and fidelity were perturbed by changing an RdRp active-site residue, K359H [25]. We did not publish this PV mutant until now because of its genetic instability (Figure 1A). We realized that this genetic instability might actually represent an opportunity, as the biochemical phenotypes reverted by the second-site suppressors might highlight the biochemical properties driving viral fitness.
Lys-359 protonates the pyrophosphate leaving group during nucleotidyl transfer [25,49]. The K359H substitution will change the efficiency of that protonation event for two reasons. First, the pKa value of histidine is substantially lower than lysine and will therefore not be protonated to the same extent. Second, the distance between the nucleotide phosphates and the imidazole ring of the histidine will be greater than that of the amino group of lysine. It is likely that the distance is more of an issue, because K359R PV is stable and K359R RdRp exhibits the same reduction in catalytic efficiency as K359H RdRp [7]. The changes of I331F and P356S in motifs C and D, respectively, may hypothetically serve to move His-359 closer to the nucleotide phosphates.
Each substitution (IF and PS) individually increased the replication efficiency of K359H PV, with the combination of the two substitutions producing an additive outcome (Figure 1C). Each substitution also decreased the fidelity of nucleotide addition in cells, as assessed by ribavirin sensitivity, to the same extent; the combination mutant yielded a virus with properties on par with wild type (Figure 2C). The observed changes to replication speed and fidelity for the mutant viruses were explained by the biochemical properties of the corresponding RdRps (Figure 2D). As observed in the past, both speed and fidelity appear to be correlated.
The substitutions at the active site of the polymerase do not only affect replication efficiency and population diversity but also impact other stages of the lifecycle. The kinetics of virus production for the different mutants did not correspond directly to the kinetics of replication (compare Figure 2A to Figure 1C). Virus production was significantly delayed for PS-KH PV (Figure 2A). This issue with virus production manifests as a substantial decrease in the specific infectivity of this virus (Figure 2B). One possible explanation for this outcome is that the RdRp is encoded by the 3D region of the viral genome and this region is also a component of the PV 3CD protein. The 3CD protein has well established roles in aspects of the lifecycle before, during, and after replication, including virion morphogenesis [29,42,50,51]. Observations such as these highlight the additional level of complexity associated with establishing a cause-and-effect relationship between biochemical properties of the RdRp and fitness, virulence, and/or pathogenesis.
We and others have observed correlations between RdRp fidelity and recombination efficiency, with higher fidelity suppressing recombination and vice versa [19,20,21,22,23,24]. An unexpected outcome of this study was the observation that IF-KH and PS-KH PVs exhibit substantially different propensities for recombination in cells (Figure 3B), in spite of very similar speed and fidelity phenotypes in cells (Figure 1C and Figure 2C) and in vitro (Figure 2D). PS-KH PV appeared to be completely incapable of supporting recombination in cells (Figure 3B). Because the cell-based recombination assay requires the ability of virus to plaque, it was conceivable that the defect of PS-KH PV was more a reflection of virus spread than impaired recombination. We established a more sensitive fluorescence-based assay (Figure 3C), which also indicated impaired recombination (Figure 3D). The impact of the PS substitution on recombination was also quite evident in the context of IF-PS-KH PV. The biological and biochemical properties of this triple mutant are within twofold of wild type (Figure 1 and Figure 2), but the recombination efficiency of this mutant is down by 30-fold relative to WT PV (Figure 3B). Together, these observations support the existence of biochemical properties of the RdRp other than speed and fidelity that are essential for recombination and impaired by changing Pro-356 to Ser. Further characterization of these mutants will be required to identify this undefined biochemical property.
Our laboratory has had a longstanding interest in the mechanism of recombination [19,20,23,32,45]. We have shown that the high poly (rU) polymerase activity of PV RdRp derives from template switching during elongation, thus mimicking copy-choice recombination (Figure 4A) [32]. A second assay that we developed is based on template switching from the end of template, thus mimicking forced-copy-choice recombination (Figure 5A) [23,45]. Based on the mutants and corresponding RdRp derivatives reported here, high fidelity impairs both types of recombination (see K359H in Figure 4B and Figure 5D). Interestingly, unique RdRp determinants exist for each mode of recombination. PS-KH RdRp is impaired for copy-choice recombination (Figure 4B); IF-KH is impaired for forced-copy-choice recombination (Figure 5D). Impairment of recombination by PS-KH PV in cells is consistent with template switching as the primary mechanism of copy-choice recombination in cells, as demonstrated by Kirkegaard and Baltimore [44]. The observation that IF-KH PV fails to exhibit a wild-type recombination phenotype (Figure 3B) may suggest that forced-copy-choice recombination contributes to recombination in cells as well.
The current narrative of the characterized PV mutants impaired for recombination is that recombination acts as a mechanism to purge deleterious mutations in the viral population [52,53,54]. Our study reveals several caveats for the interpretation of these previous studies. There are known and unknown biochemical properties of the RdRp that determine recombination efficiency, and there are at least two distinct mechanisms of recombination. Without understanding the biochemical properties of the mutants under investigation, it is difficult to compare one study to another that may explain disparate observations made between laboratories [30,55]. Further analysis of the panel of mutants reported here may help to resolve the controversy but should illuminate how changes to the mechanism and efficiency of recombination impact viral evolution in a background of constant speed and fidelity.
Viral fitness, virulence, and pathogenesis are determined collectively by the structure, dynamics, and activity of all virus-encoded functions. The conserved structure, dynamics, and mechanism of the viral RdRp makes this enzyme an attractive target for development of attenuated viruses by changing conserved residues capable of perturbing conserved, biochemical function [7,9,11]. As a field, we have focused on RdRp fidelity [7] but have come to the realization that RdRp speed should not be ignored [18]. This study shows that speed, fidelity, and undefined biochemical properties of the RdRp exist and contribute to both the mechanism and efficiency of recombination (Figure 6). This inextricable connection of the myriad biochemical properties of the RdRp precludes attribution of a single biochemical property to viral fitness, virulence, and pathogenesis (Figure 6). The field will be served best by continued emphasis on discovery of manipulatable functions of the RdRp, and other viral enzymes, instead of debating the importance of individual properties.

Author Contributions

Conceptualization, C.E.C. and J.J.A.; Methodology, A.W. and J.J.A.; Validation, H.K., A.W., V.D.E., Y.Z., and J.J.A.; Formal analysis, H.K., A.W., V.D.E., Y.Z., J.J.A., and C.E.C.; Investigation, H.K., A.W., V.D.E., Y.Z., and J.J.A.; Data curation, H.K., and Y.Z.; Writing—original draft preparation, H.K. and A.W.; Writing—review and editing, H.K., A.W., V.D.E., Y.Z., J.J.A., and C.E.C.; Visualization, H.K., A.W., V.D.E., and J.J.A.; Supervision, C.E.C., J.J.A., and A.W.; Project administration, C.E.C. and J.J.A.; Funding acquisition, C.E.C. and J.J.A.

Funding

This study was supported by a grant (AI045818) from NIAID, NIH to C.E.C. and A.W. is the recipient of a fellowship (18POST33960071) from the American Heart Association (AHA).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. NIAID Emerging Infectious Diseases/Pathogens. Available online: https://www.niaid.nih.gov/research/emerging-infectious-diseases-pathogens (accessed on 12 September 2019).
  2. Graci, J.D.; Cameron, C.E. Therapeutically targeting RNA viruses via lethal mutagenesis. Future Virol. 2008, 3, 553–566. [Google Scholar] [CrossRef] [PubMed]
  3. Perales, C.; Domingo, E. Antiviral Strategies Based on Lethal Mutagenesis and Error Threshold. Curr. Top. Microbiol. Immunol. 2016, 392, 323–339. [Google Scholar] [CrossRef] [PubMed]
  4. Furuta, Y.; Gowen, B.B.; Takahashi, K.; Shiraki, K.; Smee, D.F.; Barnard, D.L. Favipiravir (T-705), a novel viral RNA polymerase inhibitor. Antivir. Res. 2013, 100, 446–454. [Google Scholar] [CrossRef] [PubMed]
  5. Cameron, C.E.; Moustafa, I.M.; Arnold, J.J. Fidelity of Nucleotide Incorporation by the RNA-Dependent RNA Polymerase from Poliovirus. Enzymes 2016, 39, 293–323. [Google Scholar] [CrossRef] [PubMed]
  6. Lee, C.A.; August, A.; Arnold, J.J.; Cameron, C.E. Polymerase Mechanism-Based Method of Viral Attenuation. Methods Mol. Biol. 2016, 1349, 83–104. [Google Scholar] [CrossRef]
  7. Weeks, S.A.; Lee, C.A.; Zhao, Y.; Smidansky, E.D.; August, A.; Arnold, J.J.; Cameron, C.E. A Polymerase mechanism-based strategy for viral attenuation and vaccine development. J. Biol. Chem. 2012, 287, 31618–31622. [Google Scholar] [CrossRef]
  8. Graham, R.L.; Becker, M.M.; Eckerle, L.D.; Bolles, M.; Denison, M.R.; Baric, R.S. A live, impaired-fidelity coronavirus vaccine protects in an aged, immunocompromised mouse model of lethal disease. Nat. Med. 2012, 18, 1820–1826. [Google Scholar] [CrossRef]
  9. Vignuzzi, M.; Wendt, E.; Andino, R. Engineering attenuated virus vaccines by controlling replication fidelity. Nat. Med. 2008, 14, 154–161. [Google Scholar] [CrossRef]
  10. Li, C.; Wang, H.; Yuan, T.; Woodman, A.; Yang, D.; Zhou, G.; Cameron, C.E.; Yu, L. Foot-and-mouth disease virus type O specific mutations determine RNA-dependent RNA polymerase fidelity and virus attenuation. Virology 2018, 518, 87–94. [Google Scholar] [CrossRef]
  11. Korboukh, V.K.; Lee, C.A.; Acevedo, A.; Vignuzzi, M.; Xiao, Y.; Arnold, J.J.; Hemperly, S.; Graci, J.D.; August, A.; Andino, R.; et al. RNA virus population diversity, an optimum for maximal fitness and virulence. J. Biol. Chem. 2014, 289, 29531–29544. [Google Scholar] [CrossRef]
  12. Rozen-Gagnon, K.; Stapleford, K.A.; Mongelli, V.; Blanc, H.; Failloux, A.B.; Saleh, M.C.; Vignuzzi, M. Alphavirus mutator variants present host-specific defects and attenuation in mammalian and insect models. PLoS Pathog. 2014, 10, e1003877. [Google Scholar] [CrossRef] [PubMed]
  13. Zeng, J.; Wang, H.; Xie, X.; Li, C.; Zhou, G.; Yang, D.; Yu, L. Ribavirin-resistant variants of foot-and-mouth disease virus: The effect of restricted quasispecies diversity on viral virulence. J. Virol. 2014, 88, 4008–4020. [Google Scholar] [CrossRef] [PubMed]
  14. Gnadig, N.F.; Beaucourt, S.; Campagnola, G.; Borderia, A.V.; Sanz-Ramos, M.; Gong, P.; Blanc, H.; Peersen, O.B.; Vignuzzi, M. Coxsackievirus B3 mutator strains are attenuated in vivo. Proc. Natl. Acad. Sci. USA 2012, 109, E2294–E2303. [Google Scholar] [CrossRef]
  15. Coffey, L.L.; Beeharry, Y.; Borderia, A.V.; Blanc, H.; Vignuzzi, M. Arbovirus high fidelity variant loses fitness in mosquitoes and mice. Proc. Natl. Acad. Sci. USA 2011, 108, 16038–16043. [Google Scholar] [CrossRef] [PubMed]
  16. Vignuzzi, M.; Stone, J.K.; Arnold, J.J.; Cameron, C.E.; Andino, R. Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature 2006, 439, 344–348. [Google Scholar] [CrossRef] [PubMed]
  17. Arnold, J.J.; Vignuzzi, M.; Stone, J.K.; Andino, R.; Cameron, C.E. Remote site control of an active site fidelity checkpoint in a viral RNA-dependent RNA polymerase. J. Biol. Chem. 2005, 280, 25706–25716. [Google Scholar] [CrossRef]
  18. Fitzsimmons, W.J.; Woods, R.J.; McCrone, J.T.; Woodman, A.; Arnold, J.J.; Yennawar, M.; Evans, R.; Cameron, C.E.; Lauring, A.S. A speed-fidelity trade-off determines the mutation rate and virulence of an RNA virus. PLoS Biol. 2018, 16, e2006459. [Google Scholar] [CrossRef] [PubMed]
  19. Li, C.; Wang, H.; Shi, J.; Yang, D.; Zhou, G.; Chang, J.; Cameron, C.E.; Woodman, A.; Yu, L. Senecavirus-Specific Recombination Assays Reveal the Intimate Link between Polymerase Fidelity and RNA Recombination. J. Virol. 2019. [Google Scholar] [CrossRef]
  20. Woodman, A.; Lee, K.M.; Janissen, R.; Gong, Y.N.; Dekker, N.H.; Shih, S.R.; Cameron, C.E. Predicting Intraserotypic Recombination in Enterovirus 71. J. Virol. 2019, 93. [Google Scholar] [CrossRef]
  21. Te Velthuis, A.J.W.; Long, J.C.; Bauer, D.L.V.; Fan, R.L.Y.; Yen, H.L.; Sharps, J.; Siegers, J.Y.; Killip, M.J.; French, H.; Oliva-Martin, M.J.; et al. Mini viral RNAs act as innate immune agonists during influenza virus infection. Nat. Microbiol. 2018, 3, 1234–1242. [Google Scholar] [CrossRef]
  22. Rawson, J.M.O.; Nikolaitchik, O.A.; Keele, B.F.; Pathak, V.K.; Hu, W.S. Recombination is required for efficient HIV-1 replication and the maintenance of viral genome integrity. Nucleic Acids Res. 2018, 46, 10535–10545. [Google Scholar] [CrossRef] [PubMed]
  23. Woodman, A.; Arnold, J.J.; Cameron, C.E.; Evans, D.J. Biochemical and genetic analysis of the role of the viral polymerase in enterovirus recombination. Nucleic Acids Res. 2016, 44, 6883–6895. [Google Scholar] [CrossRef] [PubMed]
  24. Poirier, E.Z.; Mounce, B.C.; Rozen-Gagnon, K.; Hooikaas, P.J.; Stapleford, K.A.; Moratorio, G.; Vignuzzi, M. Low-Fidelity Polymerases of Alphaviruses Recombine at Higher Rates to Overproduce Defective Interfering Particles. J. Virol. 2015, 90, 2446–2454. [Google Scholar] [CrossRef] [PubMed]
  25. Castro, C.; Smidansky, E.D.; Arnold, J.J.; Maksimchuk, K.R.; Moustafa, I.; Uchida, A.; Gotte, M.; Konigsberg, W.; Cameron, C.E. Nucleic acid polymerases use a general acid for nucleotidyl transfer. Nat. Struct. Mol. Biol. 2009, 16, 212–218. [Google Scholar] [CrossRef] [PubMed]
  26. Gohara, D.W.; Ha, C.S.; Kumar, S.; Ghosh, B.; Arnold, J.J.; Wisniewski, T.J.; Cameron, C.E. Production of “authentic” poliovirus RNA-dependent RNA polymerase (3D(pol)) by ubiquitin-protease-mediated cleavage in Escherichia coli. Protein Expr. Purif. 1999, 17, 128–138. [Google Scholar] [CrossRef]
  27. Gohara, D.W.; Crotty, S.; Arnold, J.J.; Yoder, J.D.; Andino, R.; Cameron, C.E. Poliovirus RNA-dependent RNA polymerase (3Dpol): Structural, biochemical, and biological analysis of conserved structural motifs A and B. J. Biol. Chem. 2000, 275, 25523–25532. [Google Scholar] [CrossRef]
  28. Kumagai, A.; Ando, R.; Miyatake, H.; Greimel, P.; Kobayashi, T.; Hirabayashi, Y.; Shimogori, T.; Miyawaki, A. A bilirubin-inducible fluorescent protein from eel muscle. Cell 2013, 153, 1602–1611. [Google Scholar] [CrossRef]
  29. Oh, H.S.; Banerjee, S.; Aponte-Diaz, D.; Sharma, S.D.; Aligo, J.; Lodeiro, M.F.; Ning, G.; Sharma, R.; Arnold, J.J.; Cameron, C.E. Multiple poliovirus-induced organelles suggested by comparison of spatiotemporal dynamics of membranous structures and phosphoinositides. PLoS Pathog. 2018, 14, e1007036. [Google Scholar] [CrossRef]
  30. Kempf, B.J.; Peersen, O.B.; Barton, D.J. Poliovirus Polymerase Leu420 Facilitates RNA Recombination and Ribavirin Resistance. J. Virol. 2016, 90, 8410–8421. [Google Scholar] [CrossRef]
  31. Lowry, K.; Woodman, A.; Cook, J.; Evans, D.J. Recombination in enteroviruses is a biphasic replicative process involving the generation of greater-than genome length ‘imprecise’ intermediates. PLoS Pathog. 2014, 10, e1004191. [Google Scholar] [CrossRef]
  32. Arnold, J.J.; Cameron, C.E. Poliovirus RNA-dependent RNA polymerase (3Dpol) is sufficient for template switching in vitro. J. Biol. Chem. 1999, 274, 2706–2716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Crotty, S.; Maag, D.; Arnold, J.J.; Zhong, W.; Lau, J.Y.; Hong, Z.; Andino, R.; Cameron, C.E. The broad-spectrum antiviral ribonucleoside ribavirin is an RNA virus mutagen. Nat. Med. 2000, 6, 1375–1379. [Google Scholar] [CrossRef] [PubMed]
  34. Pfeiffer, J.K.; Kirkegaard, K. A single mutation in poliovirus RNA-dependent RNA polymerase confers resistance to mutagenic nucleotide analogs via increased fidelity. Proc. Natl. Acad. Sci. USA 2003, 100, 7289–7294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Severson, W.E.; Schmaljohn, C.S.; Javadian, A.; Jonsson, C.B. Ribavirin causes error catastrophe during Hantaan virus replication. J. Virol. 2003, 77, 481–488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Vignuzzi, M.; Stone, J.K.; Andino, R. Ribavirin and lethal mutagenesis of poliovirus: Molecular mechanisms, resistance and biological implications. Virus Res. 2005, 107, 173–181. [Google Scholar] [CrossRef]
  37. Sierra, M.; Airaksinen, A.; Gonzalez-Lopez, C.; Agudo, R.; Arias, A.; Domingo, E. Foot-and-mouth disease virus mutant with decreased sensitivity to ribavirin: Implications for error catastrophe. J. Virol. 2007, 81, 2012–2024. [Google Scholar] [CrossRef] [Green Version]
  38. Arias, A.; Arnold, J.J.; Sierra, M.; Smidansky, E.D.; Domingo, E.; Cameron, C.E. 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. J. Virol. 2008, 82, 12346–12355. [Google Scholar] [CrossRef] [Green Version]
  39. Pathak, H.B.; Ghosh, S.K.; Roberts, A.W.; Sharma, S.D.; Yoder, J.D.; Arnold, J.J.; Gohara, D.W.; Barton, D.J.; Paul, A.V.; Cameron, C.E. Structure-function relationships of the RNA-dependent RNA polymerase from poliovirus (3Dpol). A surface of the primary oligomerization domain functions in capsid precursor processing and VPg uridylylation. J. Biol. Chem. 2002, 277, 31551–31562. [Google Scholar] [CrossRef] [Green Version]
  40. Herold, J.; Andino, R. Poliovirus requires a precise 5′ end for efficient positive-strand RNA synthesis. J. Virol. 2000, 74, 6394–6400. [Google Scholar] [CrossRef] [Green Version]
  41. Andino, R.; Rieckhof, G.E.; Achacoso, P.L.; Baltimore, D. Poliovirus RNA synthesis utilizes an RNP complex formed around the 5′-end of viral RNA. EMBO J. 1993, 12, 3587–3598. [Google Scholar] [CrossRef]
  42. Banerjee, S.; Aponte-Diaz, D.; Yeager, C.; Sharma, S.D.; Ning, G.; Oh, H.S.; Han, Q.; Umeda, M.; Hara, Y.; Wang, R.Y.L.; et al. Hijacking of multiple phospholipid biosynthetic pathways and induction of membrane biogenesis by a picornaviral 3CD protein. PLoS Pathog. 2018, 14, e1007086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Dulin, D.; Vilfan, I.D.; Berghuis, B.A.; Hage, S.; Bamford, D.H.; Poranen, M.M.; Depken, M.; Dekker, N.H. Elongation-Competent Pauses Govern the Fidelity of a Viral RNA-Dependent RNA Polymerase. Cell Rep. 2015, 10, 983–992. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kirkegaard, K.; Baltimore, D. The mechanism of RNA recombination in poliovirus. Cell 1986, 47, 433–443. [Google Scholar] [CrossRef]
  45. Arnold, J.J.; Cameron, C.E. Poliovirus RNA-dependent RNA polymerase (3D(pol)). Assembly of stable, elongation-competent complexes by using a symmetrical primer-template substrate (sym/sub). J. Biol. Chem. 2000, 275, 5329–5336. [Google Scholar] [CrossRef] [Green Version]
  46. Peliska, J.A.; Benkovic, S.J. Fidelity of in vitro DNA strand transfer reactions catalyzed by HIV-1 reverse transcriptase. Biochemistry 1994, 33, 3890–3895. [Google Scholar] [CrossRef]
  47. Peliska, J.A.; Benkovic, S.J. Mechanism of DNA strand transfer reactions catalyzed by HIV-1 reverse transcriptase. Science 1992, 258, 1112–1118. [Google Scholar] [CrossRef]
  48. Peliska, J.A.; Balasubramanian, S.; Giedroc, D.P.; Benkovic, S.J. Recombinant HIV-1 nucleocapsid protein accelerates HIV-1 reverse transcriptase catalyzed DNA strand transfer reactions and modulates RNase H activity. Biochemistry 1994, 33, 13817–13823. [Google Scholar] [CrossRef]
  49. Castro, C.; Smidansky, E.; Maksimchuk, K.R.; Arnold, J.J.; Korneeva, V.S.; Gotte, M.; Konigsberg, W.; Cameron, C.E. Two proton transfers in the transition state for nucleotidyl transfer catalyzed by RNA- and DNA-dependent RNA and DNA polymerases. Proc. Natl. Acad. Sci. USA 2007, 104, 4267–4272. [Google Scholar] [CrossRef] [Green Version]
  50. Parsley, T.B.; Cornell, C.T.; Semler, B.L. Modulation of the RNA binding and protein processing activities of poliovirus polypeptide 3CD by the viral RNA polymerase domain. J. Biol. Chem. 1999, 274, 12867–12876. [Google Scholar] [CrossRef] [Green Version]
  51. Harris, K.S.; Reddigari, S.R.; Nicklin, M.J.; Hammerle, T.; Wimmer, E. Purification and characterization of poliovirus polypeptide 3CD, a proteinase and a precursor for RNA polymerase. J. Virol. 1992, 66, 7481–7489. [Google Scholar]
  52. Muller, H.J. The Relation of Recombination to Mutational Advance. Mutat. Res. 1964, 106, 2–9. [Google Scholar] [CrossRef]
  53. Felsenstein, J. The evolutionary advantage of recombination. Genetics 1974, 78, 737–756. [Google Scholar] [PubMed]
  54. Chao, L. Fitness of RNA virus decreased by Muller’s ratchet. Nature 1990, 348, 454–455. [Google Scholar] [CrossRef] [PubMed]
  55. Kempf, B.J.; Watkins, C.L.; Peersen, O.B.; Barton, D.J. Picornavirus RNA Recombination Counteracts Error Catastrophe. J. Virol. 2019, 93. [Google Scholar] [CrossRef] [Green Version]
  56. Furio, V.; Moya, A.; Sanjuan, R. The cost of replication fidelity in an RNA virus. Proc. Natl. Acad. Sci. USA 2005, 102, 10233–10237. [Google Scholar] [CrossRef] [Green Version]
  57. Elena, S.F.; Sanjuan, R. Adaptive value of high mutation rates of RNA viruses: Separating causes from consequences. J. Virol. 2005, 79, 11555–11558. [Google Scholar] [CrossRef] [Green Version]
  58. Furio, V.; Moya, A.; Sanjuan, R. The cost of replication fidelity in human immunodeficiency virus type 1. Proc. Biol. Sci. 2007, 274, 225–230. [Google Scholar] [CrossRef] [Green Version]
Genes 10 00968 g001 550
Figure 1. Serial passage of K359H PV in cell culture leads to additional changes to RdRp-coding sequence. (A) Changes in RdRp-coding sequence after serial passage of K359H PV. RNA was isolated from passage 2, converted to cDNA, and amplified by PCR. The electropherogram of the sequenced PCR product is shown, revealing two amino acid changes: I331F (motif C) and P356S (motif D). (B) Location of I331, P356, and K359 in structure of PV RdRp. Palm, fingers, and thumb subdomains are shown. Conserved structural motifs are colored: A, red; B, green; C, yellow; D, blue; E, purple; F, orange; G, black. (C) Replication phenotypes observed for K359H PV and KH-containing PV mutants by using a replicon assay. Replication was monitored using a subgenomic replicon firefly luciferase. Luciferase specific activity is reported in relative light units (RLU) per microgram of total protein in the extract as a function of time post-transfection. Shown is one representative data set. This experiment has been repeated more than three times with similar results.
Figure 1. Serial passage of K359H PV in cell culture leads to additional changes to RdRp-coding sequence. (A) Changes in RdRp-coding sequence after serial passage of K359H PV. RNA was isolated from passage 2, converted to cDNA, and amplified by PCR. The electropherogram of the sequenced PCR product is shown, revealing two amino acid changes: I331F (motif C) and P356S (motif D). (B) Location of I331, P356, and K359 in structure of PV RdRp. Palm, fingers, and thumb subdomains are shown. Conserved structural motifs are colored: A, red; B, green; C, yellow; D, blue; E, purple; F, orange; G, black. (C) Replication phenotypes observed for K359H PV and KH-containing PV mutants by using a replicon assay. Replication was monitored using a subgenomic replicon firefly luciferase. Luciferase specific activity is reported in relative light units (RLU) per microgram of total protein in the extract as a function of time post-transfection. Shown is one representative data set. This experiment has been repeated more than three times with similar results.
Genes 10 00968 g001
Genes 10 00968 g002 550
Figure 2. Biological and biochemical characterization of variants of K359H PV and the RdRp reveal that both substitutions are required to achieve a near wild-type phenotype. (A) One-step growth analysis. Cells were infected at MOI 10 with the following PVs: WT, IF-KH, PS-KH, and IF-PS-KH. Viral titer (pfu/mL) was plotted as a function of time post-infection. Duplicate infected samples were used for plaque assays. Error bars indicate SEM (n = 2), biological replicates. (B) Specific infectivity. Virus was isolated 24 h post-infection and used for qRT-PCR to determine genomes/mL or plaque assay to determine pfu/mL, with the quotient yielding specific infectivity, genomes/pfu. Two biological replicates performed in triplicate for qPCR analysis. Statistical analyses were performed using unpaired, two-tailed t-test (* significance level p ≤ 0.05, n.s. = non-significant). (C) Ribavirin sensitivity. HeLa cells were infected at a MOI 0.1 with each PV in the presence of 600 µM ribavirin. After a 24-h incubation at 37 °C, virus was isolated and used for plaque assay. Indicated is the titer of virus recovered in the presence of ribavirin normalized to that recovered in the absence of ribavirin. Solid bar indicates the mean of each virus yield. Error bars indicate SEM (n = 3). Statistical analyses were performed using unpaired, two-tailed t-test (** significance level p ≤ 0.005, n.s. = non-significant). (D) Biochemical analysis. The rate constant for incorporation of a single correct nucleotide, ATP, (kpol,corr) and a single incorrect nucleotide, GTP, (kpol,incorr) by each PV RdRp was performed as previously described [18]. Data are reported using one significant figure. The error reported is the standard error from the fit of the data to a single exponential [18].
Figure 2. Biological and biochemical characterization of variants of K359H PV and the RdRp reveal that both substitutions are required to achieve a near wild-type phenotype. (A) One-step growth analysis. Cells were infected at MOI 10 with the following PVs: WT, IF-KH, PS-KH, and IF-PS-KH. Viral titer (pfu/mL) was plotted as a function of time post-infection. Duplicate infected samples were used for plaque assays. Error bars indicate SEM (n = 2), biological replicates. (B) Specific infectivity. Virus was isolated 24 h post-infection and used for qRT-PCR to determine genomes/mL or plaque assay to determine pfu/mL, with the quotient yielding specific infectivity, genomes/pfu. Two biological replicates performed in triplicate for qPCR analysis. Statistical analyses were performed using unpaired, two-tailed t-test (* significance level p ≤ 0.05, n.s. = non-significant). (C) Ribavirin sensitivity. HeLa cells were infected at a MOI 0.1 with each PV in the presence of 600 µM ribavirin. After a 24-h incubation at 37 °C, virus was isolated and used for plaque assay. Indicated is the titer of virus recovered in the presence of ribavirin normalized to that recovered in the absence of ribavirin. Solid bar indicates the mean of each virus yield. Error bars indicate SEM (n = 3). Statistical analyses were performed using unpaired, two-tailed t-test (** significance level p ≤ 0.005, n.s. = non-significant). (D) Biochemical analysis. The rate constant for incorporation of a single correct nucleotide, ATP, (kpol,corr) and a single incorrect nucleotide, GTP, (kpol,incorr) by each PV RdRp was performed as previously described [18]. Data are reported using one significant figure. The error reported is the standard error from the fit of the data to a single exponential [18].
Genes 10 00968 g002
Genes 10 00968 g003 550
Figure 3. Parameters in addition to RdRp speed and fidelity contribute to recombination efficiency. (A) Schematic of the PV recombination assay used [31]. Two RNAs are used: a replication-competent subgenomic RNA lacking capsid-coding sequence (donor RNA); and a replication-incompetent full-length genomic RNA with a defective cis-acting replication element (CRE, indicated by hairpin with defective version indicated by X) in 2C-coding sequence (acceptor RNA). Co-transfection of these RNAs produces infectious virus if recombination occurs. (B) The indicated PV RdRp was engineered into both donor and accepter RNA and co-transfected into a L929 mouse fibroblast cell line. Infectious virus produced by recombination in L929 cells was determined by plaque assay using HeLa cells. Each point shown is an independent experiment reflecting the average of three replicates. When plaques could not be detected, n.d. is indicated. Mean and SEM (n = 5) are indicated. Statistical analyses were performed using unpaired, two-tailed t-test (* significance level p ≤ 0.05). (C) Schematic of a modified PV recombination assay using an acceptor RNA expressing the green fluorescent protein, UnaG [28]. Recombinants are scored by expression of green fluorescence instead of plaques, thereby increasing the sensitivity. (D) Infectious virus produced by recombination in L929 was scored for infectious virus in HeLa cells 7, 24, and 48 h post-infection for WT and at 48 h only for K359H and PS-KH. Shown is one representative image from three independent experiments.
Figure 3. Parameters in addition to RdRp speed and fidelity contribute to recombination efficiency. (A) Schematic of the PV recombination assay used [31]. Two RNAs are used: a replication-competent subgenomic RNA lacking capsid-coding sequence (donor RNA); and a replication-incompetent full-length genomic RNA with a defective cis-acting replication element (CRE, indicated by hairpin with defective version indicated by X) in 2C-coding sequence (acceptor RNA). Co-transfection of these RNAs produces infectious virus if recombination occurs. (B) The indicated PV RdRp was engineered into both donor and accepter RNA and co-transfected into a L929 mouse fibroblast cell line. Infectious virus produced by recombination in L929 cells was determined by plaque assay using HeLa cells. Each point shown is an independent experiment reflecting the average of three replicates. When plaques could not be detected, n.d. is indicated. Mean and SEM (n = 5) are indicated. Statistical analyses were performed using unpaired, two-tailed t-test (* significance level p ≤ 0.05). (C) Schematic of a modified PV recombination assay using an acceptor RNA expressing the green fluorescent protein, UnaG [28]. Recombinants are scored by expression of green fluorescence instead of plaques, thereby increasing the sensitivity. (D) Infectious virus produced by recombination in L929 was scored for infectious virus in HeLa cells 7, 24, and 48 h post-infection for WT and at 48 h only for K359H and PS-KH. Shown is one representative image from three independent experiments.
Genes 10 00968 g003
Genes 10 00968 g004 550
Figure 4. RdRp derivatives incapable of recombination in cells exhibit defects in copy-choice recombination (template switching) in vitro. (A) Products of reactions using oligo/poly(rA) as template reflect copy-choice recombination, which requires the template-switching activity of the RdRp [32]. RdRp engages primed template and initiates RNA synthesis. The RdRp-nascent RNA complex moves to new template from internal positions or as a consequence of reaching the end of template, thus creating greater-than-template-length products. (B) Specific poly (rU) polymerase activity (pmol/min/µg) of each RdRp. WT and the indicated PV RdRp variants assessed by using a poly (rU) polymerase activity assay. Error bars represent SEM (n = 3). Statistical analyses were performed using unpaired, two-tailed t-test (** significance level p ≤ 0.01, n.s. = non-significant).
Figure 4. RdRp derivatives incapable of recombination in cells exhibit defects in copy-choice recombination (template switching) in vitro. (A) Products of reactions using oligo/poly(rA) as template reflect copy-choice recombination, which requires the template-switching activity of the RdRp [32]. RdRp engages primed template and initiates RNA synthesis. The RdRp-nascent RNA complex moves to new template from internal positions or as a consequence of reaching the end of template, thus creating greater-than-template-length products. (B) Specific poly (rU) polymerase activity (pmol/min/µg) of each RdRp. WT and the indicated PV RdRp variants assessed by using a poly (rU) polymerase activity assay. Error bars represent SEM (n = 3). Statistical analyses were performed using unpaired, two-tailed t-test (** significance level p ≤ 0.01, n.s. = non-significant).
Genes 10 00968 g004
Genes 10 00968 g005 550
Figure 5. Properties of the RdRp required for copy-choice recombination in vitro are separable from those required for forced-copy-choice recombination in vitro. (A) The heteropolymeric, symmetrical, primed-template substrate (sym/sub) has been used to establish an assay for forced-copy-choice recombination in vitro [23,45]. RdRp assembles on sym/sub. ATP is added and incorporated to yield a stable elongation complex. It is this elongation complex that is monitored for extension and transfer. Transfer is strictly dependent on the presence of an acceptor RNA with complementarity to the 3′-end of donor RNA. (B) Schematic of the experimental design is indicated. Products in boxes are those observed and monitored by denaturing polyacrylamide gel electrophoresis. (C) Reaction products were resolved by electrophoresis and detected by phosphorimaging. The only regions of the gel with bands are shown; these correspond to the sym/sub primer, one-nucleotide-extended primer (n+1), four-nucleotides-extended product (strong stop) and non-templated addition of nucleotides to that product, and the transfer product. (D) Transfer products were quantified and are expressed as a percentage relative to the value observed for WT (The concentrations (µM) of n+1 formed for WT, K359H, IF-KH, PS-KH and IF-PS-KH are 0.22, 0.22, 0.32, 0.36, and 0.37, respectively. The concentrations (µM) of transfer product for WT, K359H, IF-KH, PS-KH, and IF-PS-KH are 0.011, 0.00011, 0.0016, 0.0072, and 0.01, respectively. Overall, transfer efficiency ranges from 1% to 5%. Error bars represent SEM (n = 3). Statistical analyses were performed using unpaired, two-tailed t-test (** significance level p = 0.0001).
Figure 5. Properties of the RdRp required for copy-choice recombination in vitro are separable from those required for forced-copy-choice recombination in vitro. (A) The heteropolymeric, symmetrical, primed-template substrate (sym/sub) has been used to establish an assay for forced-copy-choice recombination in vitro [23,45]. RdRp assembles on sym/sub. ATP is added and incorporated to yield a stable elongation complex. It is this elongation complex that is monitored for extension and transfer. Transfer is strictly dependent on the presence of an acceptor RNA with complementarity to the 3′-end of donor RNA. (B) Schematic of the experimental design is indicated. Products in boxes are those observed and monitored by denaturing polyacrylamide gel electrophoresis. (C) Reaction products were resolved by electrophoresis and detected by phosphorimaging. The only regions of the gel with bands are shown; these correspond to the sym/sub primer, one-nucleotide-extended primer (n+1), four-nucleotides-extended product (strong stop) and non-templated addition of nucleotides to that product, and the transfer product. (D) Transfer products were quantified and are expressed as a percentage relative to the value observed for WT (The concentrations (µM) of n+1 formed for WT, K359H, IF-KH, PS-KH and IF-PS-KH are 0.22, 0.22, 0.32, 0.36, and 0.37, respectively. The concentrations (µM) of transfer product for WT, K359H, IF-KH, PS-KH, and IF-PS-KH are 0.011, 0.00011, 0.0016, 0.0072, and 0.01, respectively. Overall, transfer efficiency ranges from 1% to 5%. Error bars represent SEM (n = 3). Statistical analyses were performed using unpaired, two-tailed t-test (** significance level p = 0.0001).
Genes 10 00968 g005
Genes 10 00968 g006 550
Figure 6. The biochemical properties of the RdRp are inextricably linked and collectively determine biological outcomes. PV RdRp is among the most extensively studied enzyme both in cells and in test tubes [5]. A vast majority of these studies have emphasized elaboration of mechanisms governing efficient incorporation of nucleotides (speed) with high specificity (fidelity). Whether or not speed is a consequence of fidelity, vice versa, or completely separable is not clear [18,43,56,57,58]. However, these two parameters are linked to recombination efficiency. Moreover, the studies reported herein point to the existence of a biochemical property other than speed and fidelity, referred to here as undefined, that contributes to recombination efficiency. Together, the observations reported herein demonstrate that the current state of the art precludes association of a single biochemical property to viral fitness, virulence, and/or pathogenesis.
Figure 6. The biochemical properties of the RdRp are inextricably linked and collectively determine biological outcomes. PV RdRp is among the most extensively studied enzyme both in cells and in test tubes [5]. A vast majority of these studies have emphasized elaboration of mechanisms governing efficient incorporation of nucleotides (speed) with high specificity (fidelity). Whether or not speed is a consequence of fidelity, vice versa, or completely separable is not clear [18,43,56,57,58]. However, these two parameters are linked to recombination efficiency. Moreover, the studies reported herein point to the existence of a biochemical property other than speed and fidelity, referred to here as undefined, that contributes to recombination efficiency. Together, the observations reported herein demonstrate that the current state of the art precludes association of a single biochemical property to viral fitness, virulence, and/or pathogenesis.
Genes 10 00968 g006
Table
Table 1. Oligonucleotides used in the study.
Table 1. Oligonucleotides used in the study.
OligonucleotideSequence (5′-3′)
3D EcoRI-RAAT TAA TTC ATC GAT GAA TTC GGG CC
3D BglII-FGGC AAA GAA GTG GAG ATC TTG GAT GC

Share and Cite

MDPI and ACS Style

Kim, H.; Ellis, V.D., III; Woodman, A.; Zhao, Y.; Arnold, J.J.; Cameron, C.E. RNA-Dependent RNA Polymerase Speed and Fidelity are not the Only Determinants of the Mechanism or Efficiency of Recombination. Genes 2019, 10, 968. https://doi.org/10.3390/genes10120968

AMA Style

Kim H, Ellis VD III, Woodman A, Zhao Y, Arnold JJ, Cameron CE. RNA-Dependent RNA Polymerase Speed and Fidelity are not the Only Determinants of the Mechanism or Efficiency of Recombination. Genes. 2019; 10(12):968. https://doi.org/10.3390/genes10120968

Chicago/Turabian Style

Kim, Hyejeong, Victor D. Ellis, III, Andrew Woodman, Yan Zhao, Jamie J. Arnold, and Craig E. Cameron. 2019. "RNA-Dependent RNA Polymerase Speed and Fidelity are not the Only Determinants of the Mechanism or Efficiency of Recombination" Genes 10, no. 12: 968. https://doi.org/10.3390/genes10120968

APA Style

Kim, H., Ellis, V. D., III, Woodman, A., Zhao, Y., Arnold, J. J., & Cameron, C. E. (2019). RNA-Dependent RNA Polymerase Speed and Fidelity are not the Only Determinants of the Mechanism or Efficiency of Recombination. Genes, 10(12), 968. https://doi.org/10.3390/genes10120968

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop