Introduction
The foot-and-mouth disease virus (FMDV), which causes severe vesicular disease in livestock, is a member of the Aphthovirus genus in the Picornaviridae family [2,3]. The FMDV has high potential for antigenic and genetic variation; based on their induction of cross-protection in host animals, seven serotypes (A, O, C, Asia1, SAT 1, SAT 2, and SAT3) of FMDV have been identified [5,20]. Additionally, advances in DNA sequencing have dramatically increased the rate at which genotypic and phenotypic variants of FMDV are identified [4].
Replication and translation of FMDV RNA occur in association with the cell membrane in the cytoplasm of infected cells. The most critical step of FMDV replication is RNA-dependent RNA synthesis by 3D polymerase, which requires a regulatory network involving viral-encoded proteins (3B and 3D), various host factors, and a non-coding structural RNA element. The 5’ UTR contains two highly structured RNA sequences; the cloverleaf, required for genome replication, and the internal ribosome entry site, which directs translation initiation[8]. RNA replication is carried out on membranous structures by viral RNA-dependent RNA polymerase in conjunction with other viral and cellular proteins and cis-acting replication element (CRE) [10,19]. The viral RNA structure is critical for several essential functions, including replication, translation, and encapsidation [14]. Determining the structure of viral RNA has broadened our understanding of its involvement in the viral infection cycle[9,24].
Replication of the FMDV is initiated by the 3B protein, which acts as a primer [7]. Uridylylation by FMDV VPg occurs in a template-dependent manner and requires a small stem loop structure in the CRE as a natural template [17, 18, 22]. This reflects a common theme among many eukaryote-infecting viruses, which have evolved a variety of mechanisms to manipulate cellular transcription and translation machinery, in many cases via elegant RNA-centered strategies.
FMDV can spread through direct or indirect contact with infected animals and related products or by long-distance airborne transmission. Such spreading can occur in an extremely rapid manner for a variety of reasons, including the small amount of virus required to initiate infection, the large amount of virus excreted by affected animals, and the multiple routes of infection [23]. Additionally, the rate by which a population of viruses, including the FMDV, evolves can be influenced by genomic mutation rates, genomic architecture, and the speed of replication and recombination [1].
The biological significance of the CRE secondary structure in the 5’ UTR of FMDV is unclear. RNA functionality arises from its ability to fold into complex 3D structures and, often, its ability to change conformations to enable different functions, such as binding other types of RNA molecules or proteins. Interestingly, an extensive pan-flavivirus sequence analysis proposed that duplicated or repeated RNA motifs are associated with the acquisition of multiple hosts during viral evolution [9]. Experimental data from a dengue virus model indicated that RNA replication allows viruses to accumulate mutations that are beneficial in one host but deleterious in another, conferring robustness during host switching [24]. Additional studies have revealed the existence of different dengue virus RNA structures in two types of hosts[6, 15, 25]. FMDV infection is usually easier in cows than in pigs, but FMDV replication and transmission are more rapid in infected pigs than in infected cows. Since 2010, FMDV epidemics in Korea seem to be sensitive and specific to one type of host. In this study, we determined that the secondary structure of the FMDV CRE plays a role in determining the host specificity of an infection. These results, thus, assist in shedding light on FMDV evolution and host adaptation.
Materials and Methods
FMDV stock production
Viruses were isolated from LFBK (porcine kidney), ZZR127 (goat fetal tongue epithelium) and BHK-21 (baby hamster kidney) cell lines obtained from the ATCC (LGC Standard). Viurs isolation was performed according to the OIE manual (http://www.oie.int/fileadmin/Home/eng/ Health_standards/tahm/2.01.05_FMD.pdf). To produce an amplified stock of inoculum for viral serotypes A and O, two cell batches per serotype were inoculated with plaque forming units (PFUs) of field-collected vesicular fluid. For amplification of FMDV, cells were harvested after 24 hr infection and subsequently repeated freezing-thawing three times. The supernatants were collected in tubes and an equal volume of MEM with 25 mM HEPES was added to each tube before freezing at -70°C.
RT-PCR and FMDV RNA sequencing
FMDV RNA was extracted using an automatic RNA extraction machine (MagNA Pure 96, Roche), according to the manufacturer’s instructions. RNA was stored at -70℃ until use. cDNA was synthesized using a PrimeScript™ II 1st strand cDNA Synthesis Kit (TAKARA). Briefly, a 10 μl reaction mixture was prepared containing 10 mM dNTP, 1 μl oligo dT primer (50 μM), 3 μl RNase Free dH2O, and 5 μl viral RNA. The mixture was incubated for 5 min at 65°C then cooled immediately on ice. Next, the reaction was mixed with 10 μl of a second reaction mixture containing 5X PrimeScript II buffer, 0.5 μl 40 U/μl RNase inhibitor, 1μl enzyme, and RNase-free dH2O. The mixture was then incubated at 42°C for 45 min, then 70°C for 15 min.
The entire genome was amplified using AccuPower® ProFi Taq PCR PreMix (BIONEER, Korea), according to the manufacturer’s instructions, with nine overlapping pairs of FMDV-specific primers. RT-PCR products were analyzed by QIAxcel (Qiagen).
Purified PCR products were either sequenced directly or after cloning into the pGEM-T easy vector (Promega, USA). DNA sequencing was performed using an automatic DNA sequencer (ABI 3730) using the BigDye Terminator v3.1 cycle sequencing kit (ABI, USA). Analyses of sequence identity and divergence were carried out using BioEdit software (version 7.2.5.). PCR product sequences were assembled with SeqMan Pro software (DNASTAR, Inc., Madison, WI, USA) using default parameters.
Analysis of sequence arrangement
The FMDV genomic sequence was confirmed based on our sequencing results and the NCBI database. Viral gene sequences were arranged using a ClustalW multiple sequence alignment of full FMDV genome sequences. Amino acid sequence alignments of FMDV genes were also performed.
Phylogenetic analysis
Aligned sequences were used to construct a phylogenetic tree via the neighbor-joining method using MEGA 2.1 software. Evolutionary distances were calculated using Kimura’s two parameter model. The reliability of branching orders was estimated by bootstrapping (1000 replicates). The tree was rooted using the FMDV serotype O strains [O1/ Kaufbeuren (GenBank accession no: X00871), O1/Campos (AJ320488) and O2/Brescia (M55287)] as outgroups. Percentage similarities/differences were estimated using the Meg Align program from the Lasergene package (DNASTAR Inc., USA).
Modelling the RNA secondary structure
We predicted the RNA secondary structure of the FMDV CRE region using the Fold Web Server (http://rna.urmc. rochester. edu/RNAstructureWeb/Servers/Fold/Fold.html) with the lowest free energy calculation.
Results
Sample collection and sequence variation analysis of FMDV non-coding regions
Clinical isolates were collected from suspected FMDV outbreaks in Korea from 2010 to 2017 (Table 1). These isolates were classified as either A or O. Strain names were generated and recorded as serotype/outbreak place/SKR (South Korea)/outbreak year. Serotype O strains were amplified in cows and labelled O/GH/SKR/2010, O/AD/ SKR/2010, O/BE/SKR/2017, and O/JE/SKR/2017. Serotype O strains were amplified in pigs and labelled O/AD/ SKR/2010, O/YJ/SKR/2010, O/JC/SKR/2014, O/US/SKR/2014, O/GJ/SKR/2016, and O/GC/SKR/2016. Serotype A strains amplified in cows were labelled A/PC/SKR/2010 and A/YC/SKR/2017 and produced yields above 100 PFU/ mL in BHK-21 cells. Viral RNA was extracted using an RNeasy mini kit (Qiagen), and the full genome was PCR-amplified using a Qiagen one-step reverse transcriptase PCR (RT-PCR) kit and 19 overlapping pairs of FMDV-specific primers based on the sequence of O/SKR/JC/2014 (GenBank accession no. MG257782). PCR products were directly sequenced using an ABI 3730XL with the BigDye Terminator v3.1 cycle sequencing kit. PCR product sequences were assembled using ClustalW multiple sequence alignments with default parameters using BioEdit software (version 7.2.5.). A single open reading frame (ORF; protein-coding region) was predicted by comparing the O/SKR/JC/2014 sequence with the FMDV reference sequence (GenBank accession no. MG257782). Sequence homology was determined using BLAST via the NCBI website.
Table 1. Foot-and-Mouth Disease Viruses used in this study
* Unpublished sequence by APQA
After obtaining the full sequences of 11 different FMDV strains, we focused on the genetic variation of CRE sequences. Fig. 1A shows 55 nucleotide-long CRE sequences of 11 FMDV strains isolated during epidemics from 2010 to 2017 based on the O/AD/SKR/2010/C reference sequence. Genetic variation patterns of the CRE differed based on identification as either serotype O or A. All nucleotide variations at individual positions of these 55 nucleotide CRE sequences are shown in Fig. 1B. The five nucleotides, 23A24A25A26C27A, were well conserved in all virus strains.
Fig. 1. Sequence alignment and nucleotide variation of FMDV CRE regions. (A) Nucleotide variations in the 55 nucleotides of the FMDV CRE. The FMDV strains were isolated in Korea from 2010 to 2017. (B) Summary of the variations at each 55 nucleotide position.
Classification of CRE nucleotide variations based on genotype and host species
All FMDV strains isolated from 2010 to 2017 in Korea included serotypes O and A and were obtained from either cow or pig host species. To define the relationship between CRE variations, genotypes, and host species, we classified CRE nucleotide variations as shown in Fig. 2A and 2B present information regarding CRE variants classified as either serotype O or A. The CRE nucleotide sequences appeared to vary based on serotypes. For serotype O, the 13 nucleotides of the CRE sequence were either G or A, with A appearing in samples from 2014 onward. Serotype O isolated from 2010 to 2016 showed similar CRE variation patterns as the 2017 endemic FMDV isolate (Fig. 2A and 2B). Comparisons of the CRE variation between host species indicated that CRE sequences from pig isolates of FMDVs were extremely similar (Fig. 2D), while CREs from cow isolates of FMDVs showed large variations in divergence (Fig. 2C). Interestingly, comparisons of nucleotide sequences of FMDV isolates from cows, A/PC/SKR/2010/C, and two O serotypes from 2017 showed similar CRE variation patterns regardless of serotype. These results demonstrate that specific CRE variation patterns are associated with either different genotypes or host species.
Fig. 2. Sequence alignment and nucleotide variation of the FMDV CRE region classified by O/A serotypes and host species. (A) Nine virus strains had expected nucleotide sequences of the O serotype. (B) Two strains (A/PC/SKR/2010/C and A/YC/ SKR/2017/C) had expected nucleotide sequences of the A serotype. (C) Seven FMDV strains were isolated from cows. (D) Four FMDV strains from 2014 and 2016 were isolated from pigs.
Adaptive mutations of the full and CRE sequences of FMDV show grouped phylogenetic tree patterns
Phylogenetic trees were constructed based on the complete sequences (Fig. 3A) and CRE regions (Fig. 3B) of FMDV isolates. In the full sequence-based tree (Fig. 3A), SKR/2010, SKR/2014, and SKR/2016 isolates clustered closely and did not share a branch with SKR/2017. In addition, phylogenetic analysis based on CRE regions yielded similar results. Interestingly, in the phylogenetic tree constructed using CRE sequences, the strain O/GH/SKR/2010/C was more closely related to the O/GC/SKR/2016 strain compared to the phylogenetic tree constructed using total FMDV sequences. In CRE and full sequence-based trees, SKR/2016 strains were more proximal to SKR/2014 rather than SKR/2017. Results from phylogenetic analyses of serotype A isolates from 2010 to 2017 indicate high genetic homogeneity across cow, but not pig, isolates (Fig. 2), which is similar to what has been observed from published sequences from neighboring countries, confirming previous findings. These results, thus, indicate that estimations of viral evolution using CRE sequences correlate with results generated by comparisons of entire FMDV sequences.
Fig. 3. Maximum Likelihood tree showing phylogenetic relationships of FMDV isolates based on the complete genomic sequence (A) and CRE region sequence (B). Three trees were rooted using the FMDV type O strain (O/AD/SKR/2010/C). Bootstrap support values above 70% out of 1,000 replicates are shown near the major nodes. Horizontal branch lengths are drawn to scale.
Comparisons of secondary structures of FMDV CREs
Minimum free energy and thermodynamically stable secondary structures of CREs were predicted using the m-fold web server (http://bioweb.pasteur.fr/seqanal/interfaces/ mfold-simple.html) and revealed a single stem loop structure, similar to other lineages. Each 55 nucleotide-long CRE structure included the conserved sequence motif AAACA located within a loop at the end of a stable stem. The A23A24 A25C26A27 motif in the CRE (Domain I) is essential for virus replication, and VPg uridylylation [17] was fully conserved amongst all examined strains. Since the CRE region is essential for FMDV replication, we focused on sequence variations within the CRE region and identified multiple nucleotide differences between FMDV isolates. We converted these RNA sequences into secondary RNA structures using the Fold Web Server (http://rna.urmc.rochester.edu/RNAstructure Web/Servers/Fold/Fold.html) program and generated RNA structure predictions with the lowest free energy. The FMDV CRE region forms a stem-loop structure composed of 55 nucleotides. The 15 or 17 nucleotide-long circle loop structure consistently contained the five conserved nucleotides (AAACA) for all CRE sequences. In addition, we detected other structures, including pentagon, hexagon, and heptagon structures, within the CREs of various strains (Fig. 4). In summary, the secondary structure of the CRE is composed of a 15~17 nucleotide-long loop, a 4~6 nucleotide-long stem, two hexagons, and one additional structure of either a pentagon or a heptagon.
Fig. 4. RNA secondary structures of individual FMDV strains depend on the host species infected. The secondary structure of the CRE RNA of FMDVs was determined using the Fold Web Server with the lowest free energy calculation. The blue-colored box indicates relatively variable regions of RNA secondary structures of FMDV CRE.
Host species-associated differences in the RNA secondary structure of the FMDV CRE
Among FMDV epidemics in Korea over the last decade, those in 2014 and 2016 were primarily detected in pigs, while those in 2010, 2011, and 2017 were primarily detected in cows. For efficient FMDV replication in infected host cells, the RNA replication complex requires non-coding regulatory CREs, FMDV 3B/3D proteins, and host cell-specific replication factors. In addition, species-specific replication factors in cows or pigs may interact with the RNA replication complex of the FMDV while the CRE motif acts as a platform for the RNA replication complex. Therefore, we analyzed whether there are host species-dependent differences in the CRE secondary structure. The CRE secondary structures showed four different RNA structure patterns (Fig. 4). One CRE structural pattern was associated with pig-specific FMDV strains, including O/US/SKR/2014/Pig, O/JC/SKR/2014/Pig, O/GC/SKR/2016/Pig, and O/GJ/SKR/2016/Pig. The CRE associated with pig infections showed a 17 nucleotide-long loop, a 5 nucleotide-long stem, a heptagon, and two hexagons. CREs from cow infections showed three different RNA secondary structure patterns with different nucleotide lengths of loops and stems (Fig. 4). Additionally, heptagon structures were unique to the CREs of pig FMDV isolates and did not appear in cow isolates.
Discussion
The functional cooperation of viral replication factors and host cell proteins plays a critical role in the viral replication of several picornaviruses. Based on these observations, FMDVs may acquire a distinct mechanism for efficient viral replication of their genomic RNA dependent on regulatory RNA elements in addition to host factors and FMDV nonstructural proteins [15]. FMDV 3B protein shows weak association with the RNA replication complex of the FMDV genome, suggesting the assistance of another host cellular factor for the establishment of a strong viral replication complex. The formation of a complete functional complex is a rate-limiting step for FMDV replication [21]. For initial FMDV replication procedures, FMDV 3B and 3D proteins should recognize a cognate CRE RNA site [16]. These RNA secondary structures of FMDV CREs are crucial for acting as docking sites of host replication proteins as well as 3B and 3D viral proteins.
As the FMDV is an RNA virus with a single positive strand, it has shown high genetic variation in viral replication. For initial FMDV infection, only low levels of the virus may be required, and several rounds of viral replication with higher amounts of FMDV can transmit into a host animal with multiple virus variants. The rapid rate of FMDV replication induces an immune response in host animals within a short time, resulting in the decrease of other viral infections by related FMDV variants [12]. For the comparison of the RNA substitution of the FMDV, no equal distribution was shown within the FMDV full genome. The highest gene variation was shown in the VP1 region, resulting from the development of viral escape by vaccine usage and the host immune response. Even though there is only a small portion (less than 8%) of the VP1 gene in the entire FMDV genome, it has been used in calculations of the phylogenetic relationship of FMDV variants, since the VP1 site is important for host cell attachment and entry. In addition, VP1 plays a critical role in the induction of the immune response and the determination of serotype specificity [5,13]. The frequency of genetic variation of the non-coding region of the FMDV is lower than that of the coding region. CRE nucleotide substitutions show low rates of variation, but only one nucleotide variation within the CRE contributes to different RNA secondary structures for the recognition of CRE-binding proteins.
Here, we provide novel results regarding the epidemiological trends of FMDV outbreaks in Korea over a recent 10-year period, and provide viral determinants for host susceptibility. The endemic viruses from clinically infected animals showed different genetic variations dependent on the number of endemic years and the host species. Early identification of the host species of FMDV susceptibility can contribute toward efficient FMDV control programs, including vaccine application and quarantine countermeasures. The inter species transmission of FMDV between pig and cow might provide the genetic background of host specific susceptibility. However, in most of countries including Korea, since pig and cow farms were constituted of separated places, there were not research reports about FMDV inter-transmission of cow and pig. Generally, while cow might be highly infected by FMDV, pig could transmit FMDV stronger than cow.
The sequence-associated RNA secondary structure analysis of the FMDV CRE obtained from infected animal tissues may elucidate relationships between an outbreak and host susceptibility. Despite the presence of the same functional proteins in different host animals, the differential amino compositions of host cellular proteins may provide different CRE RNA recognition and discriminate in favor of specific host species for viral replication and transmission.
The pattern of genetic variation of the FMDV requires host adaptability for viral entry, intracellular replication, and establishment of the cytoplasmic viral assembly environment [11]. Accumulated information regarding the genetic variation of the FMDV may lead to the early determination of susceptible host species and provide molecular characteristics of emerging FMDV variants.
Acknowledgements
This work was supported by the Korean Animal and Plant Quarantine Agency (Z-1543082-2018-19-0102) and Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, Forestry (IPET) through Animal Disease Management Technology Development Program, funded by Ministry of Agriculture, Food and Rural Affairs (MAFRA) (320058-02).
The Conflict of Interest Statement
The authors declare that they have no conflicts of interest with the contents of this article.
References
- Alexandersen, S., Zhang, Z., Donaldson, A. I. and Garland, A. 2003. The pathogenesis and diagnosis of foot-and-mouth disease. J. Comp. Pathol. 129, 1-36. https://doi.org/10.1016/S0021-9975(03)00041-0
- Arzt, J., Juleff, N., Zhang, Z. and Rodriguez, L. 2011. The pathogenesis of foot and mouth disease I: viral pathways in cattle. Transbound. Emerg. Dis. 58, 291-304. https://doi.org/10.1111/j.1865-1682.2011.01204.x
- Arzt, J., Baxt, B., Grubman, M., Jackson, T., Juleff, N., Rhyan, J., Rieder, E., Waters, R. and Rodriguez, L. 2011. The pathogenesis of foot and mouth disease II: Viral pathways in swine, small ruminants, and wildlife; myotropism, chronic syndromes, and molecular virus-host interactions. Transbound. Emerg. Dis. 58, 305-326. https://doi.org/10.1111/j.1865-1682.2011.01236.x
- Carrillo, C., Lu, Z., Borca, M. V., Vagnozzi, A., Kutish, G. F. and Rock, D. L. 2007. Genetic and phenotypic variation of foot-and-mouth disease virus during serial passages in a natural host. J. Virol. 81, 11341-11351. https://doi.org/10.1128/JVI.00930-07
- Carrillo, C., Tulman, E. R., Delhon, G., Lu, Z., Carreno, A., Vagnozzi, A., Kutish, G. F. and Rock, D. L. 2005. Comparative genomics of foot-and-mouth disease virus. J. Virol. 79, 6487-6504. https://doi.org/10.1128/JVI.79.10.6487-6504.2005
- de Borba, L., Villordo, S. M., Iglesias, N. G., Filomatori, C. V., Gebhard, L. G. and Gamarnik, A. V. 2015. Overlapping local and long-range RNA-RNA interactions modulate dengue virus genome cyclization and replication. J. Virol. 89, 3430-3437. https://doi.org/10.1128/jvi.02677-14
- Ferrer-Orta, C., Agudo, R., Domingo, E. and Verdaguer, N. 2009. Structural insights into replication initiation and elongation processes by the FMDV RNA-dependent RNA polymerase. Curr. Opin. Struct. Biol. 19, 752-758. https://doi.org/10.1016/j.sbi.2009.10.016
- Gao, Y., Sun, S. and Guo, H. 2016. Biological function of Foot-and-mouth disease virus non-structural proteins and non-coding elements. Virol. J. 13, 1-17. https://doi.org/10.1186/s12985-015-0456-4
- Gritsun, T. and Gould, E. 2006. Origin and evolution of 3' UTR of flaviviruses: long direct repeats as a basis for the formation of secondary structures and their significance for virus transmission. Adv. Virus Res. 69, 203-248. https://doi.org/10.1016/S0065-3527(06)69005-2
- Herod, M. R., Ferrer-Orta, C., Loundras, E. A., Ward, J. C., Verdaguer, N., Rowlands, D. J. and Stonehouse, N. J. 2016. Both cis and trans activities of foot-and-mouth disease virus 3D polymerase are essential for viral RNA replication. J. Virol. 90, 6864-6883. https://doi.org/10.1128/JVI.00469-16
- Horsington, J. and Zhang, Z. 2007. Consistent change in the B-C loop of VP2 observed in foot-and-mouth disease virus from persistently infected cattle: Implications for association with persistence. Virus Res. 125, 114-118. https://doi.org/10.1016/j.virusres.2006.12.008
- Juleff, N., Windsor, M., Lefevre, E. A., Gubbins, S., Hamblin, P., Reid, E., McLaughlin, K., Beverley, P. C., Morrison, I. W. and Charleston, B. 2009. Foot-and-mouth disease virus can induce a specific and rapid CD4+ T-cell-independent neutralizing and isotype class-switched antibody response in naive cattle. J. Virol. 83, 3626-3636. https://doi.org/10.1128/JVI.02613-08
- Knowles, N. and Samuel, A. 2003. Molecular epidemiology of foot-and-mouth disease virus. Virus Res. 91, 65-80. https://doi.org/10.1016/S0168-1702(02)00260-5
- Liu, Y., Wimmer, E. and Paul, A. V. 2009. Cis-acting RNA elements in human and animal plus-strand RNA viruses. Biochim. Biophys. Acta 1789, 495-517. https://doi.org/10.1016/j.bbagrm.2009.09.007
- Markoff, L., Pang, X., Houng Hs, H. S., Falgout, B., Olsen, R., Jones, E. and Polo, S. 2002. Derivation and characterization of a dengue type 1 host range-restricted mutant virus that is attenuated and highly immunogenic in monkeys. J. Virol. 76, 3318-3328. https://doi.org/10.1128/JVI.76.7.3318-3328.2002
- Mason, P. W., Bezborodova, S. V. and Henry, T. M. 2002. Identification and characterization of a cis-acting replication element (cre) adjacent to the internal ribosome entry site of foot-and-mouth disease virus. J. Virol. 76, 9686-9694. https://doi.org/10.1128/JVI.76.19.9686-9694.2002
- Nayak, A., Goodfellow, I. G., Woolaway, K. E., Birtley, J., Curry, S. and Belsham, G. J. 2006. Role of RNA structure and RNA binding activity of foot-and-mouth disease virus 3C protein in VPg uridylylation and virus replication. J. Virol. 80, 9865-9875. https://doi.org/10.1128/JVI.00561-06
- Paul, A. V. and Wimmer, E. 2015. Initiation of proteinprimed picornavirus RNA synthesis. Virus Res. 206, 12-26. https://doi.org/10.1016/j.virusres.2014.12.028
- Salonen, A., Ahola, T. and Kaariainen, L. 2004, "Viral RNA replication in association with cellular membranes" in Membrane trafficking in viral replication Springer, pp. 139-173.
- Samuel, A. and Knowles, N. 2001. Foot-and-mouth disease type O viruses exhibit genetically and geographically distinct evolutionary lineages (topotypes). J. Gen. Virol. 82, 609-621. https://doi.org/10.1099/0022-1317-82-3-609
- Spear, A., Sharma, N. and Flanegan, J. B. 2008. Protein-RNA tethering: The role of poly (C) binding protein 2 in poliovirus RNA replication. Virology 374, 280-291. https://doi.org/10.1016/j.virol.2007.12.039
- Steil, B. P. and Barton, D. J. 2009. Conversion of VPg into VPgpUpUOH before and during poliovirus negative-strand RNA synthesis. J. Virol. 83, 12660-12670. https://doi.org/10.1128/JVI.01676-08
- Stenfeldt, C., Pacheco, J. M., Brito, B. P., Moreno-Torres, K. I., Branan, M. A., Delgado, A. H., Rodriguez, L. L. and Arzt, J. 2016. Transmission of foot-and-mouth disease virus during the incubation period in pigs. Front. Vet. Sci. 3, 105.
- Villordo, S. M., Filomatori, C. V., Sanchez-Vargas, I., Blair, C. D. and Gamarnik, A. V. 2015. Dengue virus RNA structure specialization facilitates host adaptation. PLoS Pathog. 11, e1004604. https://doi.org/10.1371/journal.ppat.1004604
- Villordo, S. M. and Gamarnik, A. V. 2013. Differential RNA sequence requirement for dengue virus replication in mosquito and mammalian cells. J. Virol. 87, 9365-9372. https://doi.org/10.1128/JVI.00567-13