Next Article in Journal
Serum Lactate Dehydrogenase Level One Week after Admission Is the Strongest Predictor of Prognosis of COVID-19: A Large Observational Study Using the COVID-19 Registry Japan
Next Article in Special Issue
Deoptimization of FMDV P1 Region Results in Robust Serotype-Independent Viral Attenuation
Previous Article in Journal
Multiplex RT-qPCR Application in Early Detection of Bovine Respiratory Disease in Healthy Calves
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 15
Issue 3
10.3390/v15030670
viruses-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

Evaluation of Potential In Vitro Recombination Events in Codon Deoptimized FMDV Strains

by
Edward Spinard
1,2,
Ian Fish
1,2,
Paul A. Azzinaro
1,
Monica Rodriguez-Calzada
1,3,
Ethan J. Hartwig
1,
George R. Smoliga
1,
Aishwarya Mogulothu
1,4,
Jonathan Arzt
1,
Teresa de los Santos
1,* and
Gisselle N. Medina
1,5,*
1
Plum Island Animal Disease Center (PIADC), ARS, USDA, P.O. Box 848, Greenport, NY 11944, USA
2
College of Veterinary Medicine, Kansas State University, Manhattan, KS 66506, USA
3
ORISE-PIADC Research Participation Program, Oak Ridge, TN 37831, USA
4
Department of Pathobiology and Veterinary Science, University of Connecticut, Storrs, CT 06269, USA
5
National Bio and Agro-Defense Facility (NBAF), ARS, USDA, Manhattan, KS 66502, USA
*
Authors to whom correspondence should be addressed.
Viruses 2023, 15(3), 670; https://doi.org/10.3390/v15030670
Submission received: 3 January 2023 / Revised: 10 February 2023 / Accepted: 25 February 2023 / Published: 2 March 2023
(This article belongs to the Special Issue Impact of Synonymous Mutations on the Evolution, Fitness and Pathogenesis of RNA Viruses)
Browse Figures
Versions Notes

Abstract

:
Codon deoptimization (CD) has been recently used as a possible strategy to derive foot-and-mouth disease (FMD) live-attenuated vaccine (LAV) candidates containing DIVA markers. However, reversion to virulence, or loss of DIVA, from possible recombination with wild-type (WT) strains has yet to be analyzed. An in vitro assay was developed to quantitate the levels of recombination between WT and a prospective A24-P2P3 partially deoptimized LAV candidate. By using two genetically engineered non-infectious RNA templates, we demonstrate that recombination can occur within non-deoptimized viral genomic regions (i.e., 3′end of P3 region). The sequencing of single plaque recombinants revealed a variety of genome compositions, including full-length WT sequences at the consensus level and deoptimized sequences at the sub-consensus/consensus level within the 3′end of the P3 region. Notably, after further passage, two recombinants that contained deoptimized sequences evolved to WT. Overall, recombinants featuring large stretches of CD or DIVA markers were less fit than WT viruses. Our results indicate that the developed assay is a powerful tool to evaluate the recombination of FMDV genomes in vitro and should contribute to the improved design of FMDV codon deoptimized LAV candidates.

1. Introduction

Foot-and-mouth disease (FMD) is a highly contagious disease of cloven-hoofed animals, including cows, pigs, sheep, and goats that is characterized by fever, lameness, and the formation of blisters on the toes, tongue, and utters of infected animals [1,2]. The causative agent of FMD, FMD virus (FMDV), belongs to the genus Aphthovirus of the Picornaviridae family.
The genome of FMDV is an ~8.2 kb single-stranded, positive sense RNA molecule that encodes a single open-reading frame (ORF) [3,4]. The ORF is composed of four coding regions: Leader (Lpro), P1 capsid coding region (VP4, VP2, VP3, and VP1), P2 non-structural coding region (2A, 2B, and 2C), and P3 non-structural coding region (3A, 3B1-3, 3C, and 3D) [3,4]. The ORF is flanked by the highly structured 5′ and 3′ untranslated regions [3,4]. Lpro, 2A, and 3C process the polyprotein into mature proteins [4]. The genome is replicated by the RNA-dependent RNA polymerase (RdRp) encoded by 3D [4]. In short, a complex composed of RdRp and other viral proteins synthesize a negative-strand intermediate that serves as a template for positive-strand synthesis [4,5].
Although mortality is low, outbreaks of FMD may cause severe economic losses in endemic regions. FMD is estimated to account for a yearly loss of between US $6.5 and $21 billion worldwide [6]. Control strategies in FMD endemic areas include the use of vaccination with an adjuvanted chemically inactivated FMDV preparation; however, use of these vaccines is marked by several difficulties: (1) production requires a specialized high containment biosafety level III facility to limit risk of virus release, (2) a chain of cold storage is required to prevent spoilage, and (3) multiple injections are required due to the short-lived vaccine-induced immunity (reviewed in [7]). Current vaccine improvement efforts focus on addressing these challenges and, among others, include the development of novel vaccine platforms such as live attenuated vaccines (LAVs). In general, LAVs induce more rapid and long-term protection [8]. Optimally, vaccines should also contain DIVA (Differentiating Infected from Vaccinated Animals) markers to allow for the immunological differentiation between naturally infected and previously vaccinated animals [7,9].
Potential LAV strains of FMDV have been created by deleting regions or incorporating mutations that directly target the coding and non-coding regions of the genome. For instance, the deletion of Lpro (leaderless, LLV) or specific mutations in the Lpro SAP domain rendered attenuated viruses in vitro and in vivo [10,11,12]. However, these strategies have proved insufficient, as inoculation with the live LLV vaccine did not induce a protective neutralizing antibody response in swine or cattle, which was presumably due to excessive attenuation. Additionally, the attenuation of the FMDV Lpro SAP mutant only relied on two amino acid substitutions, thus posing a significant risk of reversion given the known high error rate of the FMDV RdRp [13]. Therefore, improved FMDV LAV candidates must contain the right level of attenuation to offer protection and an increased number of attenuating mutations to reduce the chances of reversion to virulence.
Recently, the creation of attenuated FMDV strains has shifted from the direct targeting of virulence factors to altering the codon bias (CB) or codon-pair bias (CPB) of the genome [14,15]. In short, all organisms exhibit a preference toward the use of certain codons or separately, codon pairs. Modifying the CB or CPB can lead to mutants harboring hundreds of synonymous mutations. Complete or partial genome codon deoptimization (CD) or codon-pair deoptimization (CPD) has been leveraged to create a multitude of attenuated viruses [16,17,18,19,20,21,22]. Indeed, CPD of the FMDV P1 capsid coding region of strain A12 and CD of the P2 and/or P3 non-structural (NS) coding region of FMDV A24 has led to attenuated strains both in vitro and in vivo [14,15].
Still, historically, a major criticism against the use of LAVs is the possibility of recombination with circulating WT viruses [23]. Furthermore, DIVA markers would be rendered useless if they were lost or passed to a WT strain via recombination. In endemic areas, there is evidence of FMDV transmission between wild animals and domesticated livestock and, since the late 2000s, next-generation sequencing (NGS) and phylogenic analysis have demonstrated extensive recombination within circulating FMDV strains [24,25,26]. The mechanism of Picornavirus recombination is believed to take place by way of template switching during genome replication; this involves the movement of the RdRp from one parental RNA genome (donor) to another genome (acceptor) midway through negative strand synthesis, thereby yielding a hybrid product [27,28]. Recombination hotspots have been identified throughout the genome with a substantial number in the conserved regions encoding NS genes [26,29,30,31,32,33]. More recently, both intra- and inter-serotypic recombination was observed after experimental infection of natural hosts [34,35,36]. Furthermore, it was demonstrated that recombination occurred during the carrier phase of superinfected cattle but not during acute co-infection [36].
The recombination of FMDV has been studied in cell culture since the 1960s using temperature-sensitive and replication-inhibitor resistant mutants [37,38,39,40,41,42]. More recently, a molecular technique was developed to quantify the levels of recombination in vitro for two positive-sense RNA viruses within the Picornaviridae family, poliovirus and senecavirus A [43,44,45]. This assay involves the co-transfection of cells with equimolar amounts of two non-infectious in vitro synthesized RNA templates. One template includes a deletion of three amino acids in RdRp’s active site (ΔGDD), rendering it non-replicative. The other contains deletions in the capsid-encoding region (ΔP1), without which the virus cannot be encapsulated and transmitted. Thus, viable virus can only be generated through genomic recombination events, i.e., the RdRp undergoes template switching from the replication competent donor template (ΔP1) to the replication incompetent acceptor template (ΔGDD), successfully complementing the deleted regions without the generation of additional lethal mutations during negative-strand replication.
In this study, the above methodology was developed for use with A24-FMDV to quantitate the levels of recombination in cell culture using two in vitro synthesized constructs: the donor template A24-ΔP1 (ΔP1) and the acceptor template A24-ΔGDD (ΔGDD). After verifying the validity and reproducibility of the assay, the effects of deoptimization on recombination within the conserved P2P3 region were tested by replacing the donor template ΔP1 with A24-ΔP1-P2P3Deopt (ΔP1-P2P3Deopt). A24-ΔP1-P2P3Deopt is based on the previously described LAV DIVA candidate strain A24-P2P3Deopt [14]. Following transfection of the two genetically engineered non-infectious RNAs, viable recombinant viruses from various time points were plaque purified and analyzed by NGS to identify the general sites of recombination and the stability of the P2P3 deoptimized regions post-recombination. The results highlighted that our in vitro assay was functional. Additionally, codon deoptimization of most of the P2P3 coding region did not completely prevent recombination, although recombinants that contained large regions of deoptimized sequences, including DIVA markers, were less fit than WT in vitro.

2. Materials and Methods

Cells. BHK-21 (Baby hamster kidney [C-13] CCL-10, ATCC) cells were grown in minimum essential media (MEM; Gibco Waltham, MA, USA) supplemented with 10% Fetal Bovine Serum Characterized (HyClone, Logan, UT, USA), 1% Anti-Anti (Gibco, Waltham, MA, USA), 1% MEM Non-Essential Amino Acids (Gibco, Waltham, MA, USA), and 10% Tryptose Phosphate Broth (Teknova, Hollister, CA, USA) at 37 °C in a 5% CO2 atmosphere. LFBK-αvβ6 (Fetal porcine kidney expressing the bovine αVβ6 integrin (canonical FMDV receptor)) cells were grown in Dulbecco’s modified Eagle’s medium (DMEM; Gibco, Waltham, MA) supplemented with 10% Fetal Bovine Serum Characterized, 1% Anti-Anti, and 1% MEM Non-Essential Amino Acids at 37 °C in 5% CO2 atmosphere [46].
Plasmid Construction and RNA Synthesis. pA24Cru and pA24Cru-NheI P2P3Deopt have been previously described [14]. pA24-ΔP1 was constructed by replacing the P1 coding region of pA24Cru-(containing unique NheI and FseI) with a 60 bp de novo synthesized fragment (GenScript, Piscataway, NJ, United States) that contained 15 bp of VP4 directly followed by 27 bp of VP1 via subcloning (Figure S1) [47]. pA24-ΔP1-P2P3Deopt was subsequently constructed from this vector by replacing the P2P3 fragment with the equivalent fragment from the previously described pA24Cru-NheI P2P3Deopt [14]. pA24-ΔGDD was constructed by replacing a BamHI/SwaI fragment of pA24Cru with a 925 nt de novo synthesized fragment (GenScript, Piscataway, NJ, United States) that contained the 9-nucleotide deletion of the GDD active site of 3D (Figure S1). Modified plasmids were digested with SwaI for linearization and viral RNAs were synthesized by in vitro transcription with T7 polymerase using a MEGAscript T7 kit (Thermo Fisher, Waltham, MA, USA) digested with the included TurboDNase and purified with a RNeasy kit (Qiagen, Hilden, Germany) following the manufacturer’s directions. RNA quantification was performed on a Nanodrop Spectrophotometer ND 1000, and integrity was analyzed using an Agilent Bioanalyzer (Santa Clara, CA, USA).
Cell-Based Recombination Assay. BHK-21 cells were seeded at 1 × 106 cells/well in a 6-well tissue culture dish (Costar or Corning, Corning, NY, USA). Following an overnight incubation, cell media was replaced by 500 µL of complete cell culture media. Cellular monolayers were then transfected with 1.7 × 1011 copy number of the indicated RNA (equivalent to 580–781 ng) in 15 µL of Lipofectamine 2000 reagent (Invitrogen, Waltham, MA, USA) diluted in 250 μL of Opti-MEM (Gibco, Waltham, MA) for 2–3 biological repeats per reaction. Complete media (1.5 mL) was added to the wells after a 3 h incubation at 37 °C in 5% CO2. Plates were incubated at 37 °C in 5% CO2 until frozen at 24 or 48 hpt.
Plaque Assay and Plaque Isolation. Virus from thawed cell lysates was titered in duplicate in LFBK-αvβ6 cells. Plaques were visualized by staining cell cultures with crystal violet, which was followed by quantitation at 48 hpt. Plaque isolation was performed for each biological repeat of the ΔP1-P2P3Deopt × ΔGDD recombinant product as previously described [35]. Briefly, 200 µL of diluted cell lysate was incubated for 1 h at 37 °C in 5% CO2, which was followed by the addition of 2 mL of agarose overlay media (1:1 ratio of 2% UltraPure Low Melting Point Agarose (Invitrogen; Thermo Fisher, Waltham, MA, USA) and 2× MEM (Gibco, Waltham, MA, USA)), supplemented with 1% Anti-Anti and 1% MEM Non-Essential Amino Acids. At 48 hpt, 200 µL of clearly detected isolated plaques were randomly chosen, quickly collected before agar solidification, diluted in 100 µL of PBS++ (Gibco, Waltham, MA, USA) and frozen. Specific plaques were chosen, which was followed by passing in monolayers of LFBK-αvβ6 in T25 flasks, and they were frozen after detection of CPE (24–48 hpt).
Viral RNA Extraction, Sequencing, and NGS analysis. FMDV full-length genome sequencing was performed as previously described [36]. In brief, FMDV RNA was extracted from either isolated plaques or total cell culture supernatants using the MagMAX Viral RNA Isolation Kit (Thermo Fisher, Waltham, MA, USA), reverse-transcribed using Superscript II (Invitrogen; Thermo Fisher, Waltham, MA, USA) using random hexameric primers, followed by Nextera XT (Illumina, San Diego, CA, USA) library construction. Reads were filtered for quality and length and competitively mapped to either A24 Cruzeiro (GenBank #AY593768) or A24-P2P3Deopt DIVA in CLC Genomics Workbench v. 21 (Qiagen; Hilden, Germany) [14].
Statistical Analysis. Data handling, analysis, and graphic representations were performed using Prism 9 (GraphPad Software, San Diego, CA, United States), Microsoft Excel (Microsoft, Redmond, WA, United States) and RStudio 4.1.2 (RStudio Team (2020). RStudio: Integrated Development for R. RStudio, PBC, Boston, MA, USA). Statistical significance was determined using one-tailed Student’s t test; p values of <0.05 were considered significant.

3. Results

Recombination Rates of ΔP1 and ΔGDD. To quantify recombination, a cell culture-based assay previously developed for two members of the Picornaviridae family, poliovirus and senecavirus A, was implemented for FMDV [43,44,45]. The FMDV serotype A24Cru infectious clone (pA24-Cruz) [48] was modified to create two vectors (Figure 1 and Table 1): (a) pA24-ΔP1 has a 2172 in-frame nucleotide deletion of the P1 coding region that leaves the first five amino acids of VP4 and the last nine amino acids of VP1 intact; (b) pA24-ΔGDD contains a nine-nucleotide deletion that removes the GDD motif from the active site of 3D, rendering RdRp unable to replicate the viral genome. Viral RNA was synthesized from both vectors to create the genomic templates ΔP1 and ΔGDD and from an unmodified pA24-Cru to create the A24-WT genome. To produce viable virus, recombination must follow the diagram detailed in Figure 2A. In short, the 3D RdRp coding region must be translated from the ΔP1 template. During negative-strand synthesis, RdRp must switch templates from the ΔP1 donor to the ΔGDD recipient between the two indicated black arrows (~3800 nt window, Figure 2A) to successfully synthesize a negative-strand template without any deletions or additional lethal mutations (Figure 2A). As expected, the transfection of A24-WT RNA produced viable infectious particles at 24 hpt (1.2 × 107 plaque-forming units (PFU)/mL ± 1.3 × 106) and 48 hpt (7.0 × 106 PFU/mL ± 7.5 × 105) and, as expected, single transfections of either non-infectious ΔP1 and ΔGDD RNAs failed to produce infectious particles (Figure 2B). In contrast, infectious particles were successfully generated following the co-transfection of both non-infectious RNAs (donor) ΔP1 and (recipient) ΔGDD. The appearance of infectious viral particles indicated successful recombination, although lower viral titers were derived by using these modified RNA molecules compared to the unmodified WT control RNAs at 24 hpt (1.7 × 104 ± 5.1 × 103 PFU/mL) and 48 hpt (5.5 × 105 ± 2.1 × 105 PFU/mL). Viral plaque morphologies were similar for WT and ΔP1 × ΔGDD recombinants examined at 24 hpt; however, by 48 h, the WT plaques were larger than those observed in the recombinants, confirming adequate fitness for viruses derived from the WT template (Figure 2C) [49]. The above assays indicate that recombination between non-infectious FMDV RNAs generated infectious virus that allowed for quantification in PFU/mL in vitro.
Recombination efficiency of ΔP1-P2P3Deopt and ΔGDD. Previously, it was demonstrated that codon deoptimization of the conserved non-structural coding P2 and P3 regions of FMDV A24Cru gave rise to attenuated viruses in cell culture, laboratory animal models, and natural hosts including swine [14] and cattle (unpublished). To determine if CD of the P2-P3 genomic regions would prevent the likelihood of recombination with respect to its parental non-deoptimized clone, the pA24-ΔP1 template mentioned above was modified to contain P2-P3 deoptimized coding sequences. Of note, although highly codon deoptimized, this construct was designed to contain stretches of non-deoptimized (WT) sequences, which were purposely preserved to maintain the predicted RNA secondary structure contained within those regions (Figure S2). RNA transfections followed the diagram detailed in Step 1 of Figure 3; the generation of a negative-strand template via successful template switching that is suitable to generate infectious particles occurs as detailed in Figure 4A. The transfection of A24-WT RNA produced viable infectious particles at 24 hpt (5.8 × 106 ± 8.5 × 104 PFU/mL) and at 48 hpt (4.7 × 106 ± 1.3 × 105 PFU/mL). No infectious particles were generated when the non-infectious RNA templates (ΔP1, ΔP1-P2P3Deopt, or ΔGDD) were transfected alone (Figure 4B). Unexpectedly, the transfection of full-length control A24-P2P3Deopt RNA did not produce infectious virus at 24 or 48 hpt; granted, the original A24-P2P3Deopt virus previously characterized in cell culture and in animals was derived after the transfection of 10,000 ng of RNA, utilized different transfection conditions (electroporation), and required several rounds of blind passes in BHK-21 cells before viable virus was observed [14]. No significant differences (p = 0.26) were detected between the number of ΔP1 × ΔGDD (1.3 × 103 ± 5.4 × 102 PFU/mL) and ΔP1-P2P3Deopt × ΔGDD (7.2 × 102 ± 2.9×102 PFU/mL) infectious recombinants generated at 24 hpt (Figure 4B). Additionally, plaque morphologies were similar between WT and the ΔP1 × ΔGDD recombinants and ΔP1-P2P3Deopt × ΔGDD recombinants at 24 hpt in BHK-21 (Figure 4C) and in porcine LFBK-αvβ6 cells (Figure 4D). Increased viral production at a later time point was detected in recombinants derived from non-deoptimized donor sequences. Specifically, at 48 hpt, there was a statistically significant 10-fold difference between the number of infectious particles generated by ΔP1 × ΔGDD (3.3 × 105 ± 1.2 × 105 PFU/mL) and ΔP1-P2P3Deopt × ΔGDD (3.3 × 104 ± 1.3 × 104 PFU/mL) (p = 0.03). At 48 hpt, plaque morphologies were again similar between the ΔP1 × ΔGDD recombinants and ΔP1-P2P3Deopt × ΔGDD recombinants in both BHK-21 and LFBK-αvβ6 cells; however, WT had larger plaques than either recombinant. These results indicated that in our system, the presence of deoptimized sequences in donor RNA P2 and P3 regions fails to abolish recombination when examined at 24 or 48 hpt.
NGS analysis of recombinants. As detailed in Figure 3 (Steps 4 and 6), RNA extractions followed by Illumina sequencing was performed to assess the sites of recombination in ΔP1-P2P3Deopt × ΔGDD recombinant plaque isolates originating from the 24 hpt (Figure 5 and Figure S3), and the 48 hpt cell lysates (Figure 6), and also from selected 48 hpt plaques passed in LFBK-αvβ6 cells (Figure 7). Non-infectious ΔP1-P2P3Deopt and ΔGDD transcribed templates were included as controls. Following full-length genome sequencing, using CLC Genomics Workbench, reads were competitively mapped to three reference genomic regions that represented either the WT or deoptimized reference genomes: a 256 nt region of 2A and 2B, a 1413 nt region including 2C and 3A, and a 1059 nt region of 3C and 3D (graphically illustrated in Figure S3). The depth of coverage for reads that preferentially mapped to the WT (pink) or deoptimized (green) reference at each position are indicated in the graphs (Figure 5, Figure 6, Figure 7 and Figure S3). As expected, reads derived from supernatant transfected with the RNA templates ΔP1-P2P3Deopt or ΔGDD mapped exclusively to the transfected template (Figure 5). NGS revealed that all plaques isolated from the 24 hpt sample had recombinant genomes with at least sub-consensus levels of deoptimized sequence; however, most of the genomes within twelve plaques (ES06…08, ES10…12, and ES-14…ES-19) were composed of a WT sequence at the consensus level at all locations within the genome (Figure 5 and Figure S3). Only three plaques (ES10, ES17, and ES19) isolated from the 24 hpt transfection had an extremely low amount (maximum depth of coverage of 2 reads/position) of sub-consensus deoptimized sequence present in the 2AB region (Figure S3). Comparatively, the average read depth mapped to the corresponding WT sequence within each sample ranged between 38 and 253 reads/nt. Similarly, within the 2C-3A region, fifteen plaques isolated from the 24 hpt transfections had low levels of sub-consensus deoptimized sequence (Figure S3). Within the 3CD coding region, evidence of recombinants composed of genomes with mixed WT and deoptimized sequence was much more apparent with all but one plaque (ES06) containing deoptimized sequence. Indeed, multiple plaque isolates from the 24 hpt transfection (ES01…05, ES09, and ES20) contained a deoptimized sequence at the consensus level at the 3′ end of the 3CD region (Figure 5).
Plaques isolated from the 48 hpt transfection did not contain any deoptimized sequence in the 2AB or 2C-3A region (data not shown). Like the plaques isolated from samples acquired at 24 hpt, multiple sites of obvious recombination were recorded in the 3CD region and, within the 3′ end of the indicated area, plaques ES08, ES10, ES11, ES14, ES16, ES17, and ES20 were composed of a deoptimized sequence at the consensus level (Figure 6). In contrast, as previously seen with plaques isolated at 24 hpt, the genomes within plaques ES01, ES02, ES03, ES04, ES05, ES06, ES07, ES09, ES12, ES13, ES15, ES18, and ES19 were composed of consensus-level WT sequences in the entire 1059 nt region (Figure 6).
To determine the stability of the deoptimized sequences, plaques ES08, ES10, ES11, ES14, ES15, ES16, ES17, and ES20 isolated from the 48 hpt infection were passed in cell culture. Plaques ES01 and ES02 were also included as a control. Total cell culture supernatants were collected, and viral RNA was sequenced and analyzed as detailed above. As expected, after passage, samples ES01 and ES02 remained almost exclusively WT (Figure 7). Interestingly, plaques ES15 and ES17 were composed of consensus-level WT sequence losing all but an extremely low level of deoptimized sequences (relative to WT 0.01% for ES15 and 0.16% for ES17). In contrast, the other samples (ES08, ES10, ES11, ES14, ES16, and ES20) retained the deoptimized codons. Furthermore, the breakpoint nt position in these samples (ES08, ES10, ES11, ES14, ES16, and ES20) remained within 3–80 nt of the previously observed boundaries (Table S1), and the relative depth of coverage between WT and deoptimized reads did not change, indicating that these samples did not gain WT or deoptimized sequences.

4. Discussion

The work presented herein demonstrates that the in vitro recombination assay previously developed for poliovirus and senecavirus A can be adapted to quantify and evaluate recombination events for FMDV [43,44,45].
The co-transfection of ΔP1-P2P3Deopt × ΔGDD produced approximately the same number of viable particles as co-transfection of ΔP1 × ΔGDD RNAs at 24 hpt, indicating that recombination could still occur despite the presence of CD sequences. Furthermore, our data suggest that even if recombination occurred between a WT strain and the A24-P2P3Deopt LAV candidate in vitro, the resulting recombinant would not display a gain of fitness. Illumina sequencing of the ΔP1-P2P3Deopt × ΔGDD recombinants revealed that at 24 or 48 hpt, most recombinants were fully composed of WT sequences at the consensus level and only very few recombinants contained deoptimized sequences at the consensus level. Furthermore, even at the sub-consensus level, recombinants were much more likely to contain deoptimized sequences in the 3CD coding region, as only a few recombinants isolated at 24 hpt contained deoptimized sequences in the 2AB and 2C-3A coding regions (Figure S3). Thus, our assay demonstrates that the DIVA markers were not passed and maintained; all recombinants that contained consensus level CD contained the deoptimized sequence downstream of the DIVA markers. Of note, only one recombinant (ES06) isolated at 48 hpt contained the 3D DIVA markers at an appreciable coverage, although it was still at a sub-consensus level. There are conflicting reports regarding the influence of sequence homology in viral recombination [50]. Thus, the absence of plaque isolated recombinants containing consensus levels of deoptimized sequences in 2AB and 2C-3A coding regions may not be simply linked to a decrease in sequence homology compared to WT. More likely, recombinants did not contain deoptimized sequence at the consensus level in the 2AB and 2C-3A coding regions because large regions of deoptimization can cause an attenuation in FMDV [14,15]. The assay was specifically performed in a cell line that lacked an intact innate immune system to limit the number of host-associated variables. Still, A24-P2P3Deopt virus has been shown to have a slower growth rate in BHK-21 cells compared to the parental WT [14]. How codon deoptimization affects a virus is a very complex issue and has not been thoroughly examined in FMDV. In other viruses, codon deoptimization has been shown to cause a decrease in the rate of translation, aberrant folding of the nascent amino acid chain during translation, a decreased replication rate, an increase in dinucleotide frequencies CpG and UpA, and destabilization of conserved RNA structures [17]. Perhaps a combination of these factors affected the recovery of viable A24-P2P3 deoptimized, but not A24 WT viruses, after the transfection of full-length viral RNA in our in vitro assay. Furthermore, the ability to generate a less deoptimized genome via recombination of the two non-infectious RNA molecules (ΔP1-P2P3Deopt × ΔGDD RNA) may have facilitated template switching and the overall rate of replication and translation required for the generation of viable virus from non-defective full genomes [51]. Hence, the absence of plaque isolated recombinants containing large regions of deoptimized sequence at the consensus level could have been caused by a number of different and difficult to dissect molecular mechanisms, the elucidation of which, was beyond the scope of this study. Future studies targeted to understand those mechanisms will hopefully answer these questions.
Due to the molecular biology of FMDV, certain limitations were inherent to the study. In the recoding design of the original A24-P2P3deopt virus, codon deoptimization of FMDV P2 and P3 was excluded from regions that were predicted to form highly conserved secondary RNA structures. These regions have been confirmed to be critical for viral genome replication, irrespective of FMDV serotype [52]. Within the 3D coding region, only the first 272 of the 470 total codons were deoptimized. Accordingly, ΔP1-P2P3Deopt and ΔGDD contain a 185 nt homologous sequence between the boundary of deoptimized 3D and the ΔGDD site. If recombination were to occur in this region, a genome containing only WT sequence would be generated. We could further conclude that the emergence of WT sequence via the re-optimization of all deoptimized codons is extremely unlikely. In fact, we have previously demonstrated that the genome of the A24-P2P3Deopt virus was stable after several passages in vitro [14].
Alternatively, instead of the formation of recombinants with full-length genomes, the generation of infectious particles could be explained by the presence of defective interfering (DI) particles, which is a phenomena that has been previously described for FMDV [53]. Noteworthy, these DI particles only became apparent after multiple passages (over 200) in cell culture without innate immune pressure. In our study, the assay determines the products of recombination after only one blind passage post-RNA transfection in cell culture. Furthermore, the absence of reads containing either ΔP1 or ΔGDD sequences suggests that no DI particles were present in our samples.
In our assay, we cannot assume that all infectious particles titered at 24 or 48 hpt arose from a unique transfection event. Given the fast replication kinetics of FMDV, it is possible that infectious particles generated early post-transfection underwent multiple rounds of replication, thus inflating the total observed number of viable viruses. While plaque sequencing results indicate that all ΔP1-P2P3Deopt × ΔGDD RNA recombinants contained unique genome sequences, i.e., different breakpoints, a much larger sampling size would be needed to definitively support this hypothesis. An additional complication is that multiple viruses can contribute to a single viral plaque formation. This could result in the observation of multiple sub-consensus genomes within the purified plaques [54]. Lastly, recombination has been described as a frequent and continuous process during viral replication for picornaviruses and other positive-sense RNA viruses [50,55]. Although FMDV has not been investigated as deeply as other picornavirus family members, a similar situation could also explain the detection of multiple sub-consensus genomes in the isolated FMDV recombinants.
Additional time points and NGS that produces long reads could further refine our assay. As demonstrated in poliovirus, when imprecise recombination occurs, duplications within the genome are generated. Additional recombination events can later refine and correct the introduced errors, thus generating fitter viruses [56]. For FMDV, the NGS data indicate that even when a large portion of the viral genome sequence was deoptimized, recombination could still occur within relatively short unchanged sequence stretches, deriving fitter WT genomes. It is unclear from our limited sampling if the full-length WT virus emerged from one or multiple recombination events within the homologous region to eventually outcompete the recombinants containing deoptimized sequence or if, alternatively, the recombinants containing deoptimized sequence gained more WT sequence through multiple recombination events.
Additionally, in its current iteration, our assay could only quantify and sequence recombination events that should have resulted in viable virus. However, based on the reported frequency of recombination during replication of poliovirus [56], it is likely that a higher number of recombination events may occur within FMDV sequences, as compared to those detected in this study.
It is possible that the recently described RNA secondary structures within 3D contributed to the higher levels of recombination in the non-deoptimized 3D region [52]. Indeed, in the 1980s, early software predicted the formation of RNA stem loop structures in recombination hotspots observed during FMDV recombination in cell culture; however, other more recent studies have demonstrated conflicting data in other viruses [50,57]. Thus, it is unclear if the essential RNA structures in 3D contributed to a global pattern of higher recombination frequency in that area or if the potential decrease in the stability of RNA structures within the P2P3 region may have contributed to a lower level of recombination in those areas. Additional studies are needed to determine the role of RNA secondary structures during FMDV recombination.
Although informative, care must be taken when interpreting the results from the assay presented in this manuscript and linking it to what could occur during a co-infection of two strains of FMDV in a host. Deoptimization of the P2 and P3-coding regions resulted in attenuated phenotypes in at least two FMDV susceptible hosts of importance, swine, and cattle [14] and unpublished]. Of course, the in vitro assay described in this manuscript does not consider any recombination that would allow for the newly generated virus to escape the innate or adaptive immune system of the host (e.g., swapping capsids, gaining mutations that improve fitness, etc.) nor has attenuation through the codon deoptimization of the P2-P3 coding region in different serotypes been demonstrated. Notwithstanding, a backbone of A24 P2P3 with DIVA markers has been shown to be a reliable platform and vaccine candidate, allowing for the insertion of different serotype specific capsids [47].

5. Conclusions

Taken together, our data indicate that: (1) an in vitro assay to measure recombination for FMDV genome sequences was successfully established; (2) the majority of observed recombinants were generated by template switching that presumably occurred positionally very early during the reading of the template strand, either within the 185 nt of non-deoptimized sequence or within the first 750 nt of the 3507 nt deoptimized region in the negative strand; (3) in our system, codon deoptimization did not completely prevent recombination; (4) however, recombinants containing large regions of deoptimized regions at the consensus level were not observed and presumably displayed lower fitness than WT. Thus, our results suggest that deoptimization of the P2 and P3 coding regions of FMDV is an effective strategy to cause viral attenuation and—at least in our in vitro system—prevent the passing of DIVA markers between intra-serotype genomes. This information may be considered in the development of subsequent generations of LAVs for FMDV.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15030670/s1, Figure S1: ΔP1 and ΔGDD Cloning Fragments; Figure S2: Alignment of the P2P3 coding region originating from A24WT and A24-P2P3Deopt; Figure S3: Sequence Analysis of 24-hpt ΔP1-P2P3Deopt x ΔGDD Recombinant Plaque Isolates; Table S1: Recombinant Breakpoint Location (nt) within the 3CD Coding Region.

Author Contributions

E.S. designed the studies, performed the experiments, and analyzed the data. I.F. performed the experiments and analyzed the data. P.A.A., M.R.-C., E.J.H., G.R.S., A.M. and G.N.M. performed the experiments. J.A. contributed to data analysis and critical reading of the manuscript. E.S., G.N.M. and T.d.l.S. conceived, analyzed and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the U.S. Department of Agriculture (USDA), Current Research Information System (CRIS) project number 1940-32000-061-00D, ARS, USDA; by the Plum Island Animal Disease Research Participation Program administered by the Oak Ridge Institute for Science and Education (ORISE) through an interagency agreement between USDA and the U.S. Department of Energy; and by Non Assistance Cooperative Agreements (NACAs) number 58-8064-8-010 between ARS and the Department of Diagnostic Medicine/Pathobiology, of the College of Veterinary Medicine at Kansas State University, Manhattan, Kansas.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Acknowledgments

We thank Elizabeth Rieder for supplying the infectious clone A24 cruzeiro to derive the constructs for the in vitro recombination assays.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Animal and Plant Health Inspection Service. Foot-and-Mouth Disease. 2021. Available online: https://www.aphis.usda.gov/publications/animal_health/fs-fmd-general.pdf (accessed on 1 March 2022).
  2. Arzt, J.; Juleff, N.; Zhang, Z.; Rodriguez, L.L. The Pathogenesis of Foot-and-Mouth Disease I: Viral Pathways in Cattle. Transbound. Emerg. Dis. 2011, 58, 291–304. [Google Scholar] [CrossRef]
  3. Carrillo, C.; Tulman, E.R.; Delhon, G.; Lu, Z.; Carreno, A.; Vagnozzi, A.; Kutish, G.F.; Rock, D.L. Comparative Genomics of Foot-and-Mouth Disease Virus. J. Virol. 2005, 79, 6487–6504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Grubman, M.J.; Baxt, B. Foot-and-Mouth Disease. Clin. Microbiol. Rev. 2004, 17, 465–493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ferrer-Orta, C.; Arias, A.; Agudo, R.; Pérez-Luque, R.; Escarmís, C.; Domingo, E.; Verdaguer, N. The structure of a protein primer–polymerase complex in the initiation of genome replication. EMBO J. 2006, 25, 880–888. [Google Scholar] [CrossRef] [PubMed]
  6. Knight-Jones, T.J.D.; Rushton, J. The economic impacts of foot and mouth disease—What are they, how big are they and where do they occur? Prev. Vet. Med. 2013, 112, 161–173. [Google Scholar] [CrossRef] [Green Version]
  7. Kamel, M.; El-Sayed, A.; Vazquez, H.C. Foot-and-mouth disease vaccines: Recent updates and future perspectives. Arch. Virol. 2019, 164, 1501–1513. [Google Scholar] [CrossRef]
  8. Minor, P.D. Live attenuated vaccines: Historical successes and current challenges. Virology 2015, 479, 379–392. [Google Scholar] [CrossRef] [Green Version]
  9. Animal and Plant Health Inspection Service. NAHEMS GUIDELINES. 2015. Available online: https://www.aphis.usda.gov/animal_health/emergency_management/downloads/nahems_guidelines/nahems_guidelines_appa_vacfmdv2.pdf (accessed on 1 March 2022).
  10. Segundo, F.D.-S.; Weiss, M.; Pérez-Martín, E.; Dias, C.C.; Grubman, M.J.; de los Santos, T. Inoculation of Swine with Foot-and-Mouth Disease SAP-Mutant Virus Induces Early Protection against Disease. J. Virol. 2012, 86, 1316–1327. [Google Scholar] [CrossRef] [Green Version]
  11. Mason, P.; Piccone, M.; McKenna, T.-C.; Chinsangaram, J.; Grubman, M. Evaluation of a Live-Attenuated Foot-and-Mouth Disease Virus as a Vaccine Candidate. Virology 1997, 227, 96–102. [Google Scholar] [CrossRef]
  12. Piccone, M.E.; Rieder, E.; Mason, P.W.; Grubman, M.J. The foot-and-mouth disease virus leader proteinase gene is not required for viral replication. J. Virol. 1995, 69, 5376–5382. [Google Scholar] [CrossRef] [Green Version]
  13. Rai, D.K.; Segundo, F.D.-S.; Campagnola, G.; Keith, A.; Schafer, E.A.; Kloc, A.; de los Santos, T.; Peersen, O.; Rieder, E. Attenuation of Foot-and-Mouth Disease Virus by Engineered Viral Polymerase Fidelity. J. Virol. 2017, 91, e00081-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Diaz-San Segundo, F.; Medina, G.N.; Spinard, E.; Kloc, A.; Ramirez-Medina, E.; Azzinaro, P.; Mueller, S.; Rieder, E.; de los Santos, T. Use of Synonymous Deoptimization to Derive Modified Live Attenuated Strains of Foot and Mouth Disease Virus. Front. Microbiol. 2021, 11, 610286. [Google Scholar] [CrossRef] [PubMed]
  15. Diaz-San Segundo, F.; Medina, G.N.; Ramirez-Medina, E.; Velazquez-Salinas, L.; Koster, M.; Grubman, M.J.; de los Santos, T. Synonymous Deoptimization of Foot-and-Mouth Disease Virus Causes Attenuation in Vivo while Inducing a Strong Neutralizing Antibody Response. J. Virol. 2016, 90, 1298–1310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Osterrieder, N.; Kunec, D. Attenuation of Viruses by Large-Scale Recoding of their Genomes: The Selection Is Always Biased. Curr. Clin. Microbiol. Rep. 2018, 5, 66–72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Gonçalves-Carneiro, D.; Bieniasz, P.D. Mechanisms of Attenuation by Genetic Recoding of Viruses. Mbio 2021, 12, e02238-20. [Google Scholar] [CrossRef] [PubMed]
  18. Le Nouën, C.; Collins, P.L.; Buchholz, U.J. Attenuation of Human Respiratory Viruses by Synonymous Genome Recoding. Front. Immunol. 2019, 10, 1250. [Google Scholar] [CrossRef] [Green Version]
  19. Pintó, R.M.; Pérez-Rodríguez, F.-J.; Andrea, L.D.; de Castellarnau, M.; Guix, S.; Bosch, A. Hepatitis A Virus Codon Usage: Implications for Translation Kinetics and Capsid Folding. Cold Spring Harb. Perspect. Med. 2018, 8, a031781. [Google Scholar] [CrossRef] [PubMed]
  20. Pintó, R.M.; Aragonès, L.; Costafreda, M.I.; Ribes, E.; Bosch, A. Codon usage and replicative strategies of hepatitis A virus. Virus Res. 2007, 127, 158–163. [Google Scholar] [CrossRef]
  21. Mueller, S.; Coleman, J.R.; Papamichail, D.; Ward, C.B.; Nimnual, A.; Futcher, B.; Skiena, S.; Wimmer, E. Live attenuated influenza virus vaccines by computer-aided rational design. Nat. Biotechnol. 2010, 28, 723–726. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, Y.; Yang, C.; Song, Y.T.; Coleman, J.R.; Stawowczyk, M.; Tafrova, J.; Tasker, S.; Boltz, D.; Baker, R.; Garcia, L.; et al. Scalable live-attenuated SARS-CoV-2 vaccine candidate demonstrates preclinical safety and efficacy. Proc. Natl. Acad. Sci. USA 2021, 118, e2102775118. [Google Scholar] [CrossRef]
  23. Bull, J.J. Evolutionary reversion of live viral vaccines: Can genetic engineering subdue it? Virus Evol. 2015, 1, vev005. [Google Scholar] [CrossRef]
  24. Miguel, E.; Grosbois, V.; Caron, A.; Boulinier, T.; Fritz, H.; Cornélis, D.; Foggin, C.; Makaya, P.V.; Tshabalala, P.T.; de Garine-Wichatitsky, M. Contacts and foot and mouth disease transmission from wild to domestic bovines in Africa. Ecosphere 2013, 4, 1–32. [Google Scholar] [CrossRef] [Green Version]
  25. Omondi, G.; Alkhamis, M.A.; Obanda, V.; Gakuya, F.; Sangula, A.; Pauszek, S.; Perez, A.; Ngulu, S.; van Aardt, R.; Arzt, J.; et al. Phylogeographical and cross-species transmission dynamics of SAT1 and SAT2 foot-and-mouth disease virus in Eastern Africa. Mol. Ecol. 2019, 28, 2903–2916. [Google Scholar] [CrossRef] [PubMed]
  26. Lewis-Rogers, N.; McClellan, D.A.; Crandall, K.A. The evolution of foot-and-mouth disease virus: Impacts of recombination and selection. Infect. Genet. Evol. 2008, 8, 786–798. [Google Scholar] [CrossRef] [PubMed]
  27. Simon-Loriere, E.; Holmes, E.C. Why do RNA viruses recombine? Nat. Rev. Microbiol. 2011, 9, 617–626. [Google Scholar] [CrossRef] [PubMed]
  28. Kirkegaard, K.; Baltimore, D. The mechanism of RNA recombination in poliovirus. Cell 1986, 47, 433–443. [Google Scholar] [CrossRef] [PubMed]
  29. Brito, B.; Pauszek, S.J.; Hartwig, E.J.; Smoliga, G.R.; Vu, L.T.; Dong, P.V.; Stenfeldt, C.; Rodriguez, L.L.; King, D.P.; Knowles, N.J.; et al. A traditional evolutionary history of foot-and-mouth disease viruses in Southeast Asia challenged by analyses of non-structural protein coding sequences. Sci. Rep. 2018, 8, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Ferretti, L.; Di Nardo, A.; Singer, B.; Lasecka-Dykes, L.; Logan, G.; Wright, C.F.; Pérez-Martín, E.; King, D.P.; Tuthill, T.J.; Ribeca, P. Within-Host Recombination in the Foot-and-Mouth Disease Virus Genome. Viruses 2018, 10, 221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Lasecka-Dykes, L.; Wright, C.F.; Di Nardo, A.; Logan, G.; Mioulet, V.; Jackson, T.; Tuthill, T.J.; Knowles, N.J.; King, D.P. Full Genome Sequencing Reveals New Southern African Territories Genotypes Bringing Us Closer to Understanding True Variability of Foot-and-Mouth Disease Virus in Africa. Viruses 2018, 10, 192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bertram, M.R.; Brito, B.; Palinski, R.M.; Fish, I.H.; Pauszek, S.J.; Hartwig, E.J.; Smoliga, G.R.; Vu, L.T.; Hoang, B.H.; Phuong, N.T.; et al. Novel Recombinant Foot-and-Mouth Disease Virus Circulating in Vietnam. Microbiol. Resour. Announc. 2021, 10, e01263-20. [Google Scholar] [CrossRef]
  33. Palinski, R.M.; Brito, B.; Jaya, F.R.; Sangula, A.; Gakuya, F.; Bertram, M.R.; Pauszek, S.J.; Hartwig, E.J.; Smoliga, G.R.; Obanda, V.; et al. Viral Population Diversity during Co-Infection of Foot-and-Mouth Disease Virus Serotypes SAT1 and SAT2 in African Buffalo in Kenya. Viruses 2022, 14, 897. [Google Scholar] [CrossRef]
  34. Ferretti, L.; Perez-Martin, E.; Zhang, F.Q.; Maree, F.; de Klerk-Lorist, L.M.; van Schalkwykc, L.; Juleff, N.D.; Charleston, B.; Ribeca, P. Pervasive within-host recombination and epistasis as major determinants of the molecular evolution of the foot-and-mouth disease virus capsid. PLoS Pathog 2020, 16, e1009050. [Google Scholar] [CrossRef]
  35. Fish, I.; Stenfeldt, C.; Spinard, E.; Medina, G.N.; Azzinaro, P.A.; Bertram, M.R.; Holinka, L.; Smoliga, G.R.; Hartwig, E.J.; de Los Santos, T.; et al. Foot-and-Mouth Disease Virus Interserotypic Recombination in Superinfected Carrier Cattle. Pathogens 2022, 11, 644. [Google Scholar] [CrossRef]
  36. Arzt, J.; Fish, I.H.; Bertram, M.R.; Smoliga, G.R.; Hartwig, E.J.; Pauszek, S.J.; Holinka-Patterson, L.; Diaz-San Segundo, F.C.; Sitt, T.; Rieder, E.; et al. Simultaneous and Staggered Foot-and-Mouth Disease Virus Coinfection of Cattle. J. Virol. 2021, 95, e0165021. [Google Scholar] [CrossRef]
  37. Pringle, C. Evidence of genetic recombination in foot-and-mouth disease virus. Virology 1965, 25, 48–54. [Google Scholar] [CrossRef]
  38. Pringle, C.R. Recombination between Conditional Lethal Mutants within a Strain of Foot-and-mouth Disease Virus. J. Gen. Virol. 1968, 2, 199–202. [Google Scholar] [CrossRef]
  39. Pringle, C.R.; Slade, W.R. The Origin of Hybrid Variants Derived from Subtype Strains of Foot-and-Mouth Disease Virus. J. Gen. Virol. 1968, 2, 319–329. [Google Scholar] [CrossRef] [PubMed]
  40. MacKenzie, J.S. Virulence of temperature-sensitive mutants of foot-and-mouth disease virus. Arch. Virol. 1975, 48, 1–8. [Google Scholar] [CrossRef] [PubMed]
  41. MacKenzie, J.S.; Slade, W.R. Evidence for recombination between two different immunological types of foot-and-mouth disease virus. Aust. J. Exp. Biol. Med. Sci. 1975, 53, 251–256. [Google Scholar] [CrossRef] [PubMed]
  42. MacKenzie, J.S.; Slade, W.R.; Lake, J.; Priston, R.A.J.; Bisby, J.; Laing, S.; Newman, J. Temperature-sensitive Mutants of Foot-and-Mouth Disease Virus: The Isolation of Mutants and Observations on their Properties and Genetic Recombination. J. Gen. Virol. 1975, 27, 61–70. [Google Scholar] [CrossRef]
  43. 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] [PubMed] [Green Version]
  44. Li, C.; Wang, H.W.; Shi, J.B.; Yang, D.C.; Zhou, G.H.; Chang, J.T.; Cameron, C.E.; Woodman, A.; Yu, L. Senecavirus-Specific Recombination Assays Reveal the Intimate Link between Polymerase Fidelity and RNA Recombination. J. Virol. 2019, 93, e00576-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kempf, B.J.; Watkins, C.L.; Peersen, O.B.; Barton, D.J. Picornavirus RNA Recombination Counteracts Error Catastrophe. J. Virol. 2019, 93, e00652-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. LaRocco, M.; Krug, P.W.; Kramer, E.; Ahmed, Z.; Pacheco, J.M.; Duque, H.; Baxt, B.; Rodriguez, L.L. A Continuous Bovine Kidney Cell Line Constitutively Expressing Bovine alpha(V)beta(6) Integrin Has Increased Susceptibility to Foot-and-Mouth Disease Virus. J. Clin. Microbiol. 2013, 51, 1714–1720. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Uddowla, S.; Hollister, J.; Pacheco, J.M.; Rodriguez, L.L.; Rieder, E. A Safe Foot-and-Mouth Disease Vaccine Platform with Two Negative Markers for Differentiating Infected from Vaccinated Animals. J. Virol. 2012, 86, 11675–11685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Rieder, E.; Henry, T.; Duque, H.; Baxt, B. Analysis of a Foot-and-Mouth Disease Virus Type A 24 Isolate Containing an SGD Receptor Recognition Site in Vitro and Its Pathogenesis in Cattle. J. Virol. 2005, 79, 12989–12998. [Google Scholar] [CrossRef] [Green Version]
  49. Cottral, G.E.; Patty, R.E.; Gailiunas, P.; Scott, F.W. Relationship of foot-and-mouth disease virus plaque size on cell cultures to infectivity for cattle by intramuscular inoculation. Arch. Ges. Virusforsch 1966, 18, 276–293. [Google Scholar] [CrossRef]
  50. Bentley, K.; Evans, D.J. Mechanisms and consequences of positive-strand RNA virus recombination. J. Gen. Virol. 2018, 99, 1345–1356. [Google Scholar] [CrossRef]
  51. Groenke, N.; Trimpert, J.; Merz, S.; Conradie, A.M.; Wyler, E.; Zhang, H.W.; Hazapis, O.G.; Rausch, S.; Landthaler, M.; Osterrieder, N.; et al. Mechanism of Virus Attenuation by Codon Pair Deoptimization. Cell Rep. 2020, 31, 107586. [Google Scholar] [CrossRef]
  52. Lasecka-Dykes, L.; Tulloch, F.; Simmonds, P.; Luke, G.A.; Ribeca, P.; Gold, S.; Knowles, N.J.; Wright, C.F.; Wadsworth, J.; Azhar, M.; et al. Mutagenesis Mapping of RNA Structures within the Foot-and-Mouth Disease Virus Genome Reveals Functional Elements Localized in the Polymerase (3D(pol))-Encoding Region. Msphere 2021, 6, e0001521. [Google Scholar] [CrossRef]
  53. García-Arriaza, J.; Domingo, E.; Escarmís, C. A segmented form of foot-and-mouth disease virus interferes with standard virus: A link between interference and competitive fitness. Virology 2005, 335, 155–164. [Google Scholar] [CrossRef] [Green Version]
  54. Aguilera, E.R.; Erickson, A.K.; Jesudhasan, P.R.; Robinson, C.M.; Pfeiffer, J.K. Plaques Formed by Mutagenized Viral Populations Have Elevated Coinfection Frequencies. Mbio 2017, 8, e02020-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Dolan, P.T.; Whitfield, Z.J.; Andino, R. Mechanisms and Concepts in RNA Virus Population Dynamics and Evolution. Annu. Rev. Virol. 2018, 5, 69–92. [Google Scholar] [CrossRef]
  56. Bentley, K.; Alnaji, F.G.; Woodford, L.; Jones, S.; Woodman, A.; Evans, D.J. Imprecise recombinant viruses evolve via a fitness-driven, iterative process of polymerase template-switching events. PLOS Pathog. 2021, 17, e1009676. [Google Scholar] [CrossRef] [PubMed]
  57. Wilson, V.; Taylor, P.; Desselberger, U. Crossover regions in foot-and-mouth disease virus (FMDV) recombinants correspond to regions of high local secondary structure. Arch. Virol. 1988, 102, 131–139. [Google Scholar] [CrossRef] [PubMed]
Viruses 15 00670 g001 550
Figure 1. Schematic diagram of the FMDV RNA templates used during transfection. DIVA markers are indicated by open (3B1 deletion, RQKP9-12→PVKV) and closed (3D H27N31→YR) triangles. Deoptimized regions are indicated by stripes. ΔGDD within the active site of RdRp encoded by 3D is indicated by an X.
Figure 1. Schematic diagram of the FMDV RNA templates used during transfection. DIVA markers are indicated by open (3B1 deletion, RQKP9-12→PVKV) and closed (3D H27N31→YR) triangles. Deoptimized regions are indicated by stripes. ΔGDD within the active site of RdRp encoded by 3D is indicated by an X.
Viruses 15 00670 g001
Viruses 15 00670 g002 550
Figure 2. (A). Schematic model of FMDV recombination during negative-strand synthesis with the given templates. ΔGDD within the active site of RdRp encoded by 3D is indicated by an X. Template switching must occur between the two arrows. (B). Number of infectious particles generated 24 and 48 h post-transfection in BHK-21 cells. Error bars represent 1 standard deviation. ND—not detected. (C). Plaque morphologies of WT and recombinants in LFBK-αVβ6 cells stained at 48 hpt.
Figure 2. (A). Schematic model of FMDV recombination during negative-strand synthesis with the given templates. ΔGDD within the active site of RdRp encoded by 3D is indicated by an X. Template switching must occur between the two arrows. (B). Number of infectious particles generated 24 and 48 h post-transfection in BHK-21 cells. Error bars represent 1 standard deviation. ND—not detected. (C). Plaque morphologies of WT and recombinants in LFBK-αVβ6 cells stained at 48 hpt.
Viruses 15 00670 g002
Viruses 15 00670 g003 550
Figure 3. Experimental design schematic of the ΔP1-P2P3Deopt × ΔGDD recombination assay and the proceeding viral quantification, plaque isolation, plaque passing, and sample types that were sequenced by NGS.
Figure 3. Experimental design schematic of the ΔP1-P2P3Deopt × ΔGDD recombination assay and the proceeding viral quantification, plaque isolation, plaque passing, and sample types that were sequenced by NGS.
Viruses 15 00670 g003
Viruses 15 00670 g004 550
Figure 4. (A). Schematic model of FMDV recombination during negative-strand synthesis with the given templates. DIVA markers are indicated by open (3B1 deletion, RQKP9-12→PVKV) and closed (3D H27N31→YR) DIVA markers. Deoptimized regions are indicated by stripes. The ΔGDD within the active site of RdRp encoded by 3D is indicated by an X. Template switching must occur between the two arrows. (B). Number of infectious particles generated 24 and 48 h post-transfection in BHK-21 cells. * Statistically significant difference between the number of infectious virions produced from the different donor RNA templates. Error bars represent 1 standard deviation. ND—not detected. (C,D). Plaque morphologies of WT and recombinants in BHK-21 or LFBK-αVβ6 cells respectively, stained at 48 hpt.
Figure 4. (A). Schematic model of FMDV recombination during negative-strand synthesis with the given templates. DIVA markers are indicated by open (3B1 deletion, RQKP9-12→PVKV) and closed (3D H27N31→YR) DIVA markers. Deoptimized regions are indicated by stripes. The ΔGDD within the active site of RdRp encoded by 3D is indicated by an X. Template switching must occur between the two arrows. (B). Number of infectious particles generated 24 and 48 h post-transfection in BHK-21 cells. * Statistically significant difference between the number of infectious virions produced from the different donor RNA templates. Error bars represent 1 standard deviation. ND—not detected. (C,D). Plaque morphologies of WT and recombinants in BHK-21 or LFBK-αVβ6 cells respectively, stained at 48 hpt.
Viruses 15 00670 g004
Viruses 15 00670 g005 550
Figure 5. Sequence analysis of 24-hpt ΔP1-P2P3Deopt × ΔGDD recombinant plaque isolates and ΔP1-P2P3Deopt or ΔGDD transfection total cell supernatant (shaded blue box). Locations (x-axis) and counts (y-axis) of reads that preferentially map to either a A24-WT or A24-P2P3Deopt reference within a 1059 nt section of the 3CD coding region as indicated by the pink shaded box. Triangles indicate the location of the 3B DIVA marker (white), 3D DIVA marker (black), the boundary between deoptimized and non-deoptimized regions (red), and the GDD deletion (blue).
Figure 5. Sequence analysis of 24-hpt ΔP1-P2P3Deopt × ΔGDD recombinant plaque isolates and ΔP1-P2P3Deopt or ΔGDD transfection total cell supernatant (shaded blue box). Locations (x-axis) and counts (y-axis) of reads that preferentially map to either a A24-WT or A24-P2P3Deopt reference within a 1059 nt section of the 3CD coding region as indicated by the pink shaded box. Triangles indicate the location of the 3B DIVA marker (white), 3D DIVA marker (black), the boundary between deoptimized and non-deoptimized regions (red), and the GDD deletion (blue).
Viruses 15 00670 g005
Viruses 15 00670 g006 550
Figure 6. Sequence analysis of 48-hpt ΔP1-P2P3Deopt × ΔGDD recombinant plaque isolates. Locations (x-axis) and counts (y-axis) of reads that preferentially map to either a A24-WT or A24-P2P3Deopt reference within a 1059 nt section of the 3CD coding region as indicated by the pink shaded box. Triangles indicate the location of the 3B DIVA marker (white), 3D DIVA marker (black), the boundary between deoptimized and non-deoptimized regions (red), and the GDD deletion (blue).
Figure 6. Sequence analysis of 48-hpt ΔP1-P2P3Deopt × ΔGDD recombinant plaque isolates. Locations (x-axis) and counts (y-axis) of reads that preferentially map to either a A24-WT or A24-P2P3Deopt reference within a 1059 nt section of the 3CD coding region as indicated by the pink shaded box. Triangles indicate the location of the 3B DIVA marker (white), 3D DIVA marker (black), the boundary between deoptimized and non-deoptimized regions (red), and the GDD deletion (blue).
Viruses 15 00670 g006
Viruses 15 00670 g007 550
Figure 7. Sequence analysis of total cell culture supernatant of 48-hpt ΔP1-P2P3Deopt × ΔGDD recombinants plaque isolates passed in LFBK-αvβ6 cells. Locations (x-axis) and counts (y-axis) of reads that preferentially map to either a A24-WT or A24-P2P3Deopt reference within a 1059 nt section of the 3CD coding region as indicated by the pink shaded box. Triangles indicate the location of the 3B DIVA marker (white), 3D DIVA marker (black), the boundary between deoptimized and non-deoptimized regions (red), and the GDD deletion (blue).
Figure 7. Sequence analysis of total cell culture supernatant of 48-hpt ΔP1-P2P3Deopt × ΔGDD recombinants plaque isolates passed in LFBK-αvβ6 cells. Locations (x-axis) and counts (y-axis) of reads that preferentially map to either a A24-WT or A24-P2P3Deopt reference within a 1059 nt section of the 3CD coding region as indicated by the pink shaded box. Triangles indicate the location of the 3B DIVA marker (white), 3D DIVA marker (black), the boundary between deoptimized and non-deoptimized regions (red), and the GDD deletion (blue).
Viruses 15 00670 g007
Table
Table 1. Genotypic and phenotypic characteristics of viral RNAs.
Table 1. Genotypic and phenotypic characteristics of viral RNAs.
TemplateInfectiousReplication Competent
A24WT✓✓
A24-P2P3Deopt✓✓
ΔP1×✓
ΔP1-P2P3Deopt×✓
ΔGDD××
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Spinard, E.; Fish, I.; Azzinaro, P.A.; Rodriguez-Calzada, M.; Hartwig, E.J.; Smoliga, G.R.; Mogulothu, A.; Arzt, J.; de los Santos, T.; Medina, G.N. Evaluation of Potential In Vitro Recombination Events in Codon Deoptimized FMDV Strains. Viruses 2023, 15, 670. https://doi.org/10.3390/v15030670

AMA Style

Spinard E, Fish I, Azzinaro PA, Rodriguez-Calzada M, Hartwig EJ, Smoliga GR, Mogulothu A, Arzt J, de los Santos T, Medina GN. Evaluation of Potential In Vitro Recombination Events in Codon Deoptimized FMDV Strains. Viruses. 2023; 15(3):670. https://doi.org/10.3390/v15030670

Chicago/Turabian Style

Spinard, Edward, Ian Fish, Paul A. Azzinaro, Monica Rodriguez-Calzada, Ethan J. Hartwig, George R. Smoliga, Aishwarya Mogulothu, Jonathan Arzt, Teresa de los Santos, and Gisselle N. Medina. 2023. "Evaluation of Potential In Vitro Recombination Events in Codon Deoptimized FMDV Strains" Viruses 15, no. 3: 670. https://doi.org/10.3390/v15030670

APA Style

Spinard, E., Fish, I., Azzinaro, P. A., Rodriguez-Calzada, M., Hartwig, E. J., Smoliga, G. R., Mogulothu, A., Arzt, J., de los Santos, T., & Medina, G. N. (2023). Evaluation of Potential In Vitro Recombination Events in Codon Deoptimized FMDV Strains. Viruses, 15(3), 670. https://doi.org/10.3390/v15030670

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