JVI Figure table search 04
Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Carrillo, C.
Right arrow Articles by Rock, D. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Carrillo, C.
Right arrow Articles by Rock, D. L.

 Previous Article  |  Next Article 

Journal of Virology, May 2005, p. 6487-6504, Vol. 79, No. 10
0022-538X/05/$08.00+0     doi:10.1128/JVI.79.10.6487-6504.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.

Comparative Genomics of Foot-and-Mouth Disease Virus{dagger}

C. Carrillo,1 E. R. Tulman,1,3,4 G. Delhon,1,2 Z. Lu,1 A. Carreno,1 A. Vagnozzi,1,3,4 G. F. Kutish,1,3,4 and D. L. Rock1,3,4*

Plum Island Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Greenport, New York 11944,1 Area of Virology, School of Veterinary Sciences, University of Buenos Aires, 1427 Buenos Aires, Argentina,2 Department of Pathobiology and Veterinary Science,3 Center of Excellence for Vaccine Research, University of Connecticut, Storrs, Connecticut 062694

Received 4 August 2004/ Accepted 18 January 2005


   ABSTRACT
Top
 Abstract
Introduction
 Materials and Methods
 Results and Discussion
 References
 
Here we present complete genome sequences, including a comparative analysis, of 103 isolates of foot-and-mouth disease virus (FMDV) representing all seven serotypes and including the first complete sequences of the SAT1 and SAT3 genomes. The data reveal novel highly conserved genomic regions, indicating functional constraints for variability as well as novel viral genomic motifs with likely biological relevance. Previously undescribed invariant motifs were identified in the 5' and 3' untranslated regions (UTR), as was tolerance for insertions/deletions in the 5' UTR. Fifty-eight percent of the amino acids encoded by FMDV isolates are invariant, suggesting that these residues are critical for virus biology. Novel, conserved sequence motifs with likely functional significance were identified within proteins Lpro, 1B, 1D, and 3C. An analysis of the complete FMDV genomes indicated phylogenetic incongruities between different genomic regions which were suggestive of interserotypic recombination. Additionally, a novel SAT virus lineage containing nonstructural protein-encoding regions distinct from other SAT and Euroasiatic lineages was identified. Insights into viral RNA sequence conservation and variability and genetic diversity in nature will likely impact our understanding of FMDV infections, host range, and transmission.


   INTRODUCTION
Top
 Abstract
 Introduction
Materials and Methods
 Results and Discussion
 References
 
Foot-and-mouth disease (FMD) is an acute, systemic disease of domestic and wild cloven-hooved animal species and is caused by Foot-and-mouth disease virus (FMDV), the prototype member of the genus Aphthovirus of the family Picornaviridae. FMDV occurs as seven distinct serotypes (Euroasiatic serotypes A, O, C, and Asia1 and South African Territories [SAT] serotypes SAT1, SAT2, and SAT3) and multiple subtypes reflecting significant genetic variability. The highly contagious nature of FMDV and the associated productivity losses make it a primary animal health concern worldwide. Effective vaccines and stringent control measures have enabled FMD eradication in most developed countries, which maintain unvaccinated, seronegative herds in compliance with strict international trade policies. However, the disease remains enzootic in many regions of the world, posing a serious problem for commercial trade with FMD-free countries.

FMDV is a small nonenveloped virus with a pseudo T=3 icosahedral capsid made up of 60 copies each of four structural proteins. The capsid surrounds an 8.4-kilobase, positive-sense, single-stranded RNA genome which is covalently bound at its 5' end to the small viral protein 3B (or VPg) and is polyadenylated at its 3' end. Upon virus entry into a cell, the viral genome is rapidly translated into a polyprotein which is co- and posttranslationally cleaved by viral proteinases into several partially cleaved, likely functional, intermediates and ultimately into 12 mature proteins (87, 97).

The FMDV genome organization is similar to that of other picornaviruses, including a large single open reading frame (ORF) flanked by highly structured 5' and 3' untranslated regions (5' UTR and 3' UTR, respectively) (Fig. 1A). The 5' UTR consists of, from the 5' end, a 350- to 380-nucleotide (nt) "short" (S) fragment, a 100- to 420-nt poly(C) tract (90% C), and the approximately 700-nt 5' terminus of the genomic "long" (L) fragment, which contains three or four tandemly repeated pseudoknots, a stem-loop cis-acting replication element (cre), and a type II internal ribosome entry site (IRES) (69). The FMDV 5' UTR plays important roles in cap-independent translation initiation of the viral polyprotein and in viral genome replication (64). The 3' UTR is about 90 nt long and is thought to contain cis-acting elements required for efficient genome replication (4).




View larger version (64K):
[in this window]
[in a new window]
  		FIG. 1. FMDV UTR variability. (A) FMDV S and L genomic fragments showing the locations and predicted secondary structures of UTRs. Cn, poly(C) tract; AUG, start codons; pAn, poly(A); I, II, and III, S fragment, IRES, and 3' UTR, respectively. (B) Rates of substitution per nt between positions 1 and 1218 and positions 8277 and 8441 in a genomic ClustalW alignment. Arrows with numbers indicate positions with high rates of nt substitution. (C) 5' UTR nt conservation between positions 1 and 1153 for the 103 FMDV isolates examined in this study. Motifs previously defined as conserved between picornavirus IRESs are underlined. Boxes indicate primers used for amplification. Asterisks indicate nt positions tolerant to insertions/deletions in two or more isolates. (D) Nucleotide conservation within S fragment (I), IRES (II), and 3' UTR (III) predicted secondary structures. Arrows indicate regions used for PCR amplification. UTR regions from the isolate C5Argentina/69 were used to generate secondary structure predictions and plots with mfold and Squiggles. Black circles indicate invariant nt among the 103 isolates presented here; gray circles represent nt that are conserved in at least 100 of 103 isolates (96%); gray shading indicates residues comprising previously described motifs. Nucleotide position numbers correspond to the consensus sequence generated by a ClustalW alignment of the 103 isolates.

 
The four capsid proteins, 1A, 1B, 1C, and 1D (also known as VP4, VP2, VP3, and VP1, respectively), are encoded by the N-terminal half of the ORF, and with the exception of 1A, which is excluded from the virion surface, are involved in antigenicity and binding to a subset of RGD-dependent integrins and heparan sulfate proteoglycan receptors on the cell surface (reviewed in reference 45). Nonstructural proteins represent about two-thirds of the ORF-encoded proteins and include Lpro, 2A, 2B, 2C, 3A, 3B, 3Cpro, and 3Dpol (93). FMDV polyprotein processing is mediated by Lpro, 3Cpro, and 2A. Lpro is a papain-like protease that, in addition to excising itself from the polyprotein, cleaves the cellular translation initiation factor eIF4G, resulting in a shutoff of host cap-dependent translation (65). 3Cpro, a member of the trypsin family of serine proteinases, performs all but three of the cleavages leading to mature viral proteins and also cleaves host cell proteins (111). FMDV 2A mediates autocleavage at its C terminus, apparently by inducing a ribosomal skip during polyprotein synthesis (19, 20). Although the functions of the FMDV 2B and 2C proteins are unknown, preliminary work suggests that, similar to those of other picornaviruses, they localize to endoplasmic reticulum (ER)-derived vesicles, the sites of viral genome replication (106). 3A is thought to be a multifunctional integral membrane protein that enhances viral RNA synthesis by 3Dpol and stimulates cleavage of the 3CD precursor (91). FMDV encodes three nonidentical copies of genome-linked 3B, a protein required for viral RNA replication (26). The 3D gene encodes the viral RNA-dependent RNA polymerase, and it and 3A colocalize with ER membrane-associated replication complexes (80, 92).

Despite a basic understanding of many aspects of picornaviral biology, much information regarding FMDV UTR, protein, and protein precursor functions and roles in virulence, host range, and virus transmission remains poorly understood. Comparative genomic analysis of a large number of FMDV genomes may allow the identification of highly conserved regulatory or coding regions which are critical for aspects of virus biology. To date, the few comparative analyses of full-length FMDV genomes have relied on intratypic (serotype O) or intertypic (serotypes O, A, and C) comparisons of a small number of isolates and serotypes (70, 86). In addition, complete genomic sequences of some serotypes are not available (SAT1 and SAT3) or, for others, only represent highly cell culture-adapted isolates (serotypes A and C). Studies involving large numbers of virus isolates have largely focused on genomic regions encoding structural proteins associated with serotype specificity for phylogenetic purposes (reviewed in reference 55).

For this work, we obtained and analyzed 103 complete genome sequences representing all FMDV serotypes, including the previously unavailable SAT1 and SAT3 genomes. Our analyses identified novel highly conserved genomic regions indicating functional constraints for variability, novel viral genomic motifs with likely biological relevance, and previously undescribed virus lineages.


   MATERIALS AND METHODS
Top
 Abstract
 Introduction
 Materials and Methods
Results and Discussion
 References
 
FMDV isolates. One hundred three FMDV isolates were obtained from the Animal Plant Health Inspection Service, United States Department of Agriculture (Table 1). Except for recent isolates from Argentina, the United Kingdom, and South Korea (2000-2001), most isolates belonged to the Animal Plant Health Inspection Service archive collection and represented a broad temporal and geographical distribution (Table 1). Although certain archive isolates have been adapted to guinea pigs or passaged in cell cultures, isolates obtained from their natural hosts were well represented. For some analyses, complete genome or whole polyprotein FMDV sequences currently available in GenBank release 141 were considered and included (Table 1).


View this table:
[in this window]
[in a new window]
  		TABLE 1. FMDV sequences

 
FMDV RNA isolation, RT-PCR amplification, and sequencing. Full-length FMDV genome sequences were obtained by reverse transcription (RT) of the complete viral genomic RNA followed by redundant amplification and sequencing of overlapping cDNA fragments spanning the entire viral genome. The total RNA was extracted directly from 140 µl of vesicular fluid, the supernatants of infected cell cultures, or 0.1 g of epithelial tissue from vesicular lesions by the use of RNeasy Mini kits (QIAGEN) and was then resuspended in 40 µl of double-distilled H2O.

Complete viral genomes were reverse transcribed to cDNAs by use of a Super SMART PCR cDNA synthesis kit (Clontech), using the supplied modified oligo(dT) according to the manufacturer's instructions. After spin column purification, the full-length double-stranded cDNA from each isolate was amplified in eight identical long-distance PCRs (2 µl cDNA in a 50-µl reaction volume; 20 cycles of 95°C for 15 s, 65°C for 15 s, and 68°C for 2 min) using primers supplied in the kit. Next, 21 overlapping 1.2-kb PCR amplicons spanning the whole FMDV genome were generated in 96-well plates (2.5 µl of template from each of the eight double-stranded cDNA replicates in 25-µl reactions using the primers described in Table S1 in the supplemental material). Adjacent amplicons overlapped for >50% of their lengths. To prevent truncation by subsequent rounds of PCR, we artificially created poly(C) sequences by substituting flanking sequences with specific primers (F1R and F2F) containing six C's or G's.

Amplified products of the correct size were purified (DNA filtration system; Eppendorf) and resuspended in 40 µl of sterile double-distilled H2O per well. All amplicon replicates (eight per FMDV subgenomic fragment per isolate) were subjected to direct sequencing (2.5 µl of template per reaction) using position-specific forward and reverse primers (see Table S2 in the supplemental material), Big Dye Terminator cycle sequencing kits (Applied Biosystems), and a PRISM 3730xl automated DNA sequencer (Applied Biosystems). To overcome virus variability and PCR/sequencing errors, we sequenced each amplicon in multiple reactions with both forward and reverse primers, specifically selected by serotype (see Table S2 in the supplemental material). Since each amplicon overlapped its neighbor >50%, the multiple direct sequencing reactions resulted in a redundancy rate ranging from 13 to 48 and averaged 28 sequence events per nucleotide.

Direct DNA sequencing of amplicons derived from a given FMDV isolate yielded a master sequence representing the most probable nt for each position of the sequence. This approach prevented analyses of minor sequence variants, polymerase misincorporation errors, and sequencing ambiguities through multiple independent cDNA synthesis, PCR amplification, and direct sequencing events. Due to the quasispecies nature of FMDV populations, polymorphisms were detected at some nt positions. Nevertheless, all positions could be unambiguously assigned to a single nt due to the high degree of redundancy generated by the sequencing strategy. The following sequence ambiguity code was used: K (T/G), M (A/C), R (A/G), S (C/G), W (A/T), Y (C/T), B (C/T/G), D (A/T/G), H (A/C/T), V (A/C/G), and N (A/C/G/T).

Sequence analysis. Bases were called from chromatogram traces produced with Phred (25), which also produced a quality file containing a predicted probability of error at each base position. Viral sequences were assembled with the Phrap (http://www.phrap.org) and CAP3 assemblers (43). Gap closure was performed as described previously (3). Multiple sequence alignments were performed with the ClustalW (1.7) and Dialign (2.2) computer programs (73, 107). The positions of 5' UTR sequences presented in Results refer to the consensus sequence generated by ClustalW alignment of this region for the 103 isolates, with some manual editing. Nucleotide substitution analysis was carried out by use of the DISTREE (1.2) (103), DNArates (1.1) (105), ALISTAT (16), and PRETTY (GCG 10 software package [16]) programs. Analyses of codons and synonymous/nonsynonymous (syn/nonsyn) substitution ratios were calculated by using SNAP (77), CodonW (http://www.molbiol.ox.ac.uk/cu/), and codeml (PAML3.14 package), which was also used for statistical evaluations of heterogeneous selection pressures at amino acid sites (116). For protein analysis, the PRETTY program was used. Searches for motifs and/or signal sequences were performed with the MOTIFS (GCG package), HMM (58), pfscan (http://www.isrec.isb-sib.ch/ftp-server/pftools/pft2.2/pft2.2.tar.Z), and Blocks (36) programs of the PROSITE, Pfam, and Blocks databases. Transmembrane protein segment predictions were performed by using Memsat (50), Tmpred (41), Toppred (114), Psort (76), and Saps (11). Secondary structure predictions were performed for proteins by using GOR secondary structure prediction (29) and Pratt (48, 49) and for RNAs by using mfold (46, 117), rnafold (40), and Squiggles (83). Pseudoknot analysis was conducted with pknots (96).

Phylogenetic analysis was performed on aligned genomic and subgenomic regions of FMDV by utilizing the neighbor-joining and split decomposition methods as implemented in the Phylo_win (28) and/or SplitsTree 4.0 (http://www-ab.informatik.uni-tuebingen.de/software/jsplits/welcome_en.html) (44) software package and by utilizing maximum likelihood as implemented in Puzzle (104) and dnaml in the Phylip package (http://evolution.genetics.washington.edu/phylip.html) (27). Individual protein-coding regions were used to screen for incongruent tree topologies suggestive of genomic recombination, and split decomposition analysis (7) of genomic and subgenomic regions was utilized to graphically screen for reticulated branching patterns which were also suggestive of recombination. Similarity plots and bootscanning analysis (100), which compare a given query sequence to several reference sequences via incremental sliding sequence windows to yield corrected similarity values and bootstrap resampling frequencies, respectively, were performed as implemented in the SimPlot, v. 2.5, package (63), utilizing default settings and a window size of 400 nt. As a further test for potential recombination, Sawyer's run test (102), as implemented in GENECONV, v. 1.81 (http://www.math.wustl.edu/~sawyer), was used on Dialign-aligned FMDV polyprotein-encoding regions, using default settings.

Nucleotide sequence accession numbers. The GenBank accession numbers for the genome sequences of FMDV isolates sequenced for this study are listed in Table 1.


   RESULTS AND DISCUSSION
Top
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
References
 
With the exception of the poly(C) tract, full-length genomic sequences from 103 FMDV isolates representing all serotypes, including isolates from serotypes not previously fully sequenced (SAT1 and SAT3), were obtained and analyzed. Given that the structure and dynamics of FMDV populations in nature are complex and the relatively small number of SAT2 and SAT3 genomes presented here, it is possible that additional sequence data from field isolates will alter the status of residues and nucleotides presented below as invariant.

FMDV genomes ranged from 8,046 to 8,215 nt, consistent with previously reported genome lengths (86, 108). Although small insertions/deletions were observed in the coding region affecting Lpro, 1B, 1C, 1D, and 3A, most variability in genome size was due to insertions/deletions in the UTRs.

5' UTR. The approximately 1,300-base 5' UTR plays important roles in FMDV replication. However, the specific contribution of S fragment, poly(C) tract, and L fragment pseudoknots to FMDV biology is unknown. The cardiovirus poly(C) tract has been shown to affect virulence (23), while shortening of the FMDV poly(C) tract affected virus growth in vitro but not virulence in a suckling mouse model (95). The cre is essential for picornaviral replication and contains a conserved AAACA motif which in poliovirus functions as a template for 3Dpol-mediated uridylylation of 3B (75). Mutation of the FMDV cre stem region reduced virus RNA replication (68). The approximately 500-nt FMDV IRES, which is responsible for cap-independent polyprotein translation (60), is predicted to contain four structural domains (Fig. 1D, panel II, domains 2 to 5) that interact with cellular factors involved in host translation initiation (94). Domain 2 binds the polypyrimidine tract-binding protein at a UUUC motif. Domain 3 contains at least three loop structures that are conserved in all picornaviruses, including a GNRA motif essential for the maintenance of the IRES tertiary structure and potentially involved in long-range RNA-RNA interactions. The Y-shaped domain 4 and domain 5 interact with the host translational factors eIF4B and eIF4G (65, 99). A 22-nt spacer for ribosomal recognition separates the IRES from the first of two polyprotein start codons, and an additional 84 nt separates the first start codon from the second. Picornavirus utilization of the second start codon seems to be affected by its adjacent sequence and by sequences in the IRES (82). Unlike poliovirus, FMDV utilizes the second start codon twice as often as the first (66).

The analysis in this study of 5' UTRs from 103 FMDV isolates showed that only 12% and 33% of nt were invariant in the 5' UTR S and L fragments, respectively (Fig. 1B and C; Table 2). In pairwise comparisons, however, FMDV isolates averaged 80% and 85% nt identity for the S and L fragments, respectively, indicating a high degree of conservation between isolates. This suggests that although the FMDV 5' UTR can tolerate changes at most positions, selective pressures clearly favor overall conservation, likely for the maintenance of a functional secondary structure. Relatively high rates of substitution (>10) in the 5' UTR L fragment were observed at positions 513, 558, 615, 684, 696, and 1144 (Fig. 1B).


View this table:
[in this window]
[in a new window]
  		TABLE 2. FMDV genome variabilitya

 
5' UTR S fragment. S fragment lengths averaged 322 and 373 nt for SAT and Euroasiatic isolates, respectively, and spanned positions 1 to 393 in the 5' UTR alignment (Fig. 1B). SAT viruses contained insertions/deletions of 1 to 3 nt, with most occurring at positions 120 to 160 and 200 to 300, including SAT-specific nt present at positions 38, 39, 67, 95, 250, 350, 366, and 375. Euroasiatic isolates contained insertions/deletions of 1 to 5 nt in length (with the exception of isolates C Waldman strain 149, A Canefa 1/61, and A25 Argentina/59, which contained a 76-nt deletion between positions 153 and 228) and included Euroasiatic-specific nt present at positions 86, 123, 124, 144, 145, 155, 156, 157, 174, 181 to 191, 197 to 219, 240, 256, 277, 292 to 294, 300 to 306, and 324 (Table 3). Notably, SAT1/7 Isrl 4/62 and SAT2-3 Kenya 11/60 appeared to have the most complete S fragment sequences of those sampled, lacking both SAT and Euroasiatic-specific deletions, and contained nt at positions 142, 143, 239, 290, 291, and 307 which were absent from other viruses.


View this table:
[in this window]
[in a new window]
  		TABLE 3. Selected FMDV 5' UTR insertions/deletions

 
Although the overall nt substitution rates in the S fragment were low, specific regions of variability were observed, and only 19 nt were invariant, excluding the regions used to prime amplification (Fig. 1C and D). S fragment sequences generally grouped into SAT and Euroasiatic lineages, which averaged only 50% nt identity with each other; however, little or no correlation with other serotypic groupings was observed (data not shown). Notably, several isolates of disparate serotypes contained very similar S fragment sequences, including the C Waldman strain 149 and A12 Valle strain 119 S fragments, which shared 98% nt identity.

In most cases, secondary structure analyses of S fragments predicted a single stem-loop structure similar to that previously proposed (12, 24). However, the S fragment sequences of 16% of the isolates folded into alternative structures within similar free energy levels (e.g., strains A22 Turkey/65, A24 Argentina 65, and Asia 1 Leb 83) (data not shown).

5' UTR L fragment. The 5' UTR region of the L fragment ranged from 604 to 751 nt long. Remarkably, apparently unrelated virus isolates contained identical deletions in the L fragment (Table 3). SAT and Euroasiatic-specific L fragment insertions/deletions were identified, as were isolate-specific insertions/deletions, including an 18-nt insertion located 28 nt downstream of the poly(C) tract in strain A5 Westerwald/51.

The pseudoknot region between positions 403 and 600 was highly tolerant to changes, with no invariant nt located within the first 200 positions and insertions/deletions tolerated at each position (Fig. 1C). Previously undescribed invariant nt and motifs located between the pseudoknot region and the IRES included AGAAWYGGGACGU (positions 617 to 629), GCRCACGWAACGCGC (positions 632 to 646), and ACAAAC (positions 668 to 673) within the cre.

FMDV IRES sequences (positions 640 to 1151) showed 70 to 100% nt identity in pairwise comparisons, 47% invariant nt, and numerous invariant motifs in the predicted secondary structure domains (Fig. 1C and D). Domain 2 contained the polypyrimidine tract-binding protein motif UUUC and three previously undescribed motifs at the base, bulge (invariant GGUCUWGAG motif), and apex regions of the stem-loop. Domain 3 contained GYRA (corresponding to the conserved picornaviral GNRA motif) and CRAAA (except in isolates O1Argentina/65 and OAkesu/58) motifs in loops A and B, respectively. While the picornavirus domain 3 loop D motif ACCC was present in 99/103 FMDV isolates, the loop 3C motif ACAC was not conserved in FMDV. Domain 3 also contained novel highly conserved motifs, including the UCGUMGCGGAGCA (positions 823 to 835) motif and GRUACUGGUA and GRGACUGGUA motifs (positions 965 to 974), specific for Euroasiatic and SAT serotypes, respectively, and CUGGWGRCAGGCUAAGGAUGCCCU (positions 983 to 1006), predicted to form a prominent bulge at the base of the stem (Fig. 1C and D). Domain 4 was highly conserved, with two novel invariant motifs (GAUCUGAG [positions 1039 to 1046] and UUAAAAG [positions 1080 to 1087]) predicted to form prominent bulges in the secondary structure. In domain 5, 19 of the 21 nt were invariant (Fig. 1C and D). The conservation of these motifs suggests biological significance. Notably, the 22-nt region between the IRES and the first AUG was highly variable, possibly underlying FMDV's preferential use of the second start codon (9, 66). In summary, the FMDV 5' UTR, although tolerant for nt substitutions and insertions/deletions, showed significant conservation, especially in the IRES, where novel sequence motifs were identified.

Polyprotein region. (i) Nucleotide sequence and codon analysis. The nucleotide variability for the 103 FMDV polyprotein (ORF) and individual protein regions is summarized in Table 2. Overall, 46% of all nt were invariant, with 73% nt identity between the least similar pair of sequences. An average transition/transversion (Ts/Tv) rate of 2.4 and syn/nonsyn substitution ratio of 2.1 were observed. As expected, the region encoding 1D was the most variable, exhibiting the lowest percentage of invariant nt (21%), Ts/Tv rate (1.55), and syn/nonsyn ratio (1.03) of all regions. In fact, excluding the highly conserved 1A (VP4) protein, regions encoding structural proteins had significantly lower Ts/Tv rates, nonsyn substitutions/site, and syn/nonsyn ratios than the rest of the genome. Although the substitution rates given here are averages for each protein-coding region, specific regions or residues within each protein were observed to contain higher or lower substitution rates (Fig. 2B and C).



View larger version (38K):
[in this window]
[in a new window]
  		FIG. 2. FMDV coding region variability. (A) Schematic diagram of FMDV polyprotein coding region showing the positions of, from top to bottom, protein-encoding regions (open boxes), cleavage intermediates, and mature proteins (lines). (B) Graphic representation of rates of nt substitution per site as calculated with DNArates software. (C) Graphic representation of nonsynonymous substitutions per site as calculated with SNAP software.

 
In general, regions encoding FMDV nonstructural proteins were highly conserved between isolates, exhibiting higher percentages of invariant residues, higher Ts/Tv rates and syn/nonsyn substitution ratios, and lower overall numbers of nonsyn changes/site (Table 2; Fig. 2C). The most conserved FMDV regions were those encoding 2B and 3C, which exhibited the highest percentage of invariant nt (61% and 59%, respectively) and amino acids (76%) in the genome. Notably, the nonstructural proteins Lpro, 3A, and 3B exhibited variability comparable to that of the structural proteins 1B, 1C, and 1D (Table 4), suggesting that these proteins are subjected to selective pressures distinct from those on the other nonstructural proteins. A comparison of CODEML {omega} rates between different substitution models further indicated that four of these proteins (Lpro, 1D, 3A, and 3B) may indeed undergo diversifying selection (Table 4) (116). Notably, this was observed for relatively few codons (1.5% to 7.6%), with fewer highly significant ones (P > 0.9) tending to cluster within each protein (positions 19, 20, 22, 23, and 82 in Lpro; positions 45, 48, 142, 143, 144, and 146 in 1D; positions 44, 132, 135, 136, and 144 in 3A; positions 4 and 11 in 3B1; and positions 17, 18, and 19 in 3B2). Interestingly, the variable capsid proteins 1B and 1C did not appear in this analysis to undergo diversifying selection (Table 4).


View this table:
[in this window]
[in a new window]
  		TABLE 4. FMDV genome variability in encoded amino acid sequencesa

 
(ii) Amino acid sequence analysis. Fifty-eight percent of all amino acids (aa) were invariant among the FMDV isolates analyzed here, with the percentage of amino acid conservation for each viral protein correlating with nt conservation in the corresponding genomic region (Fig. 3; Tables 2 and 4). Positions at putative protein cleavage sites showed various degrees of sequence conservation, with Euroasiatic lineages being more conserved than SATs. Residues at the 1A/1B, 2A/2B, 2B/2C, 2C/3A, and 3B/3C cleavage sites exhibited high degrees of conservation (Table 5).



View larger version (48K):
[in this window]
[in a new window]
  		FIG. 3. FMDV invariant amino acids (single letter code). The consensus sequence is based on the comparison of 103 FMDV isolates conducted for this study, with the exception of the 1D region, which also included all 1D sequences available in GenBank (release 141).

 

View this table:
[in this window]
[in a new window]
  		TABLE 5. FMDV polyprotein cleavage site variabilitya

 
Leader proteinase. FMDV Lpro is expressed as long (Lab) and short (Lb) isoforms that result from alternative usage of either the first or second start codon, respectively, with the sequence between the two AUGs predicted to form a stable hairpin structure crucial to IRES activity (39, 66). Lpro isoforms have indistinguishable activities and specificities and have been proposed to play a role in virulence that may involve the regulation of host interferon responses (89).

Lpro lacked insertions/deletions, with the exception of two additional aa in Euroasiatic lineages (positions 22 and 23). Each methionine start codon was invariant, indicating that both Lpro isoforms are significant for aspects of FMDV biology. Consistent with the previously described hypervariability in the region flanked by the two AUGs, only one residue (C6) was invariant in this region (113). Only 44% of the residues in Lb were invariant, making it much less conserved than was previously reported (30). Although substitutions were concentrated in terminal regions of Lpro, the predicted secondary structure in these regions remained relatively unaltered (data not shown).

A detailed analysis of the Lpro/1A junction resulted in a more ambiguous cleavage motif than was previously described (Table 5) (113). Only GAGXS at the 1A N terminus and the previously described basic residue (K/R) required for cleavage at the Lpro C terminus were invariant. Residues required for Lpro catalytic activity (C52, H149, and D165) (32, 59), suggested to be involved in Lpro autocatalysis (E77) (89), and important for eIF4G cleavage (H110 and H139) (90) were invariant, except for an H110D substitution in all SAT viruses, an E77Q substitution in O6 Pirbright/65, and an E77K substitution in A Phillipines/75 (Table 6). The present FMDV analysis revealed that only 43 of the 65 residues previously identified as conserved between FMDV O1 Kaufbeuren/66 Lpro and the aphthovirus equine rhinitis A virus Lpro were invariant, with 35 of the 43 residues concentrated in three distinct regions (positions 44 to 63, 95 to 110, and 133 to 185) (38) (Table 6). Additional, previously undescribed invariant FMDV Lpro motifs are shown in Table 6.


View this table:
[in this window]
[in a new window]
  		TABLE 6. Amino acid conservation in FMDV and ERAV Lpro

 
These data indicate that despite a relatively low overall Lpro aa sequence conservation among the FMDV isolates examined here, functional elements such as the catalytic and cleavage sites and additional, previously undescribed motifs are invariant. The variability observed at the Lpro N and C termini may affect specific virus-host interactions, including ribosomal recognition of alternative start codons.

Structural proteins. FMDV structural proteins are involved in capsid assembly and stability, virus binding to target cells, and antigenic specificity, influencing significant aspects of virus infection and immunity (45). The high level of variability in FMDV external capsid proteins observed here likely reflects the selective pressures on them.

(i) 1A (VP4). 1A was the most conserved FMDV protein, with 81% of the aa being invariant (Fig. 3; Table 4), including the N-terminal myristylation site and a previously identified heterotypic swine and bovine T-cell epitope (1A positions 20 to 35) (10). Interestingly, the Q73 residue in the SAT viruses distinguishes them from Euroasiatic lineages. Residues potentially specific for SAT2 and SAT3 (I76) and for SAT1 (V80) were also identified.

(ii) 1B (VP2). 1B, a protein of 218 or 219 aa, plays a critical role in virion structural stability and maturation (14). N-terminal regions of 1B contain previously undescribed motifs. For example, DKKTEETTLLEDRILTTRNGHT(T/I)STTQSSVG and DKKTEETT(L/H)LEDRI(L/M/V)TT(S/R)H(G/N)TTTSTTQSSVG are conserved among Euroasiatic and SAT virus isolates, respectively. Three previously identified T-cell epitopes (positions 48 to 68, 114 to 132, and 179 to 187) were conserved (88). Notably, only the N-terminal half of each epitope was invariant, while their C-terminal regions were highly variable.

(iii) 1C (VP3). 1C, a protein of 219 to 221 aa, contains important conformational neutralizing epitopes and makes significant contributions to capsid stability (62). Among the FMDV isolates examined here, 1C was highly variable, with only 39% of the aa being invariant. Most amino acid substitutions were concentrated in four regions (positions 55 to 88, 130 to 140, 176 to 186, and 196 to 208), with insertions/deletions present in the first two regions and a previously described T-cell epitope present in the second region (88). The 1B/1C cleavage site was nearly invariant in Euroasiatic lineages but much less conserved among SAT isolates (Table 5).

(iv) 1D (VP1). 1D is the most studied FMDV protein due to its significance for virus attachment and entry, protective immunity, and serotype specificity (45). 1D ranges in size from 217 to 221 aa, with insertions/deletions contained in regions 140-150 and 166-170. Overall, only 26% of the 1D residues are invariant. The invariant residues and motifs are shown in Fig. 3. Given the available data and the proposed significance of the conserved RGD motif for virus reception and pathogenesis (71, 72), it is notable that complete integrin-binding RGD motifs were lacking in 9 of the 103 FMDV isolates examined here (GGD in C Waldman strain 149 and C4 Tierra del Fuego/66, TGD in C5 Argentina/69, RGE in Asia1/2 Isrl 3/63 and A21 Kenya/64, RDD in A25 Argentina/59 and A Canefa 1/61, HGD in Asia1/3 Kimrom, and PGD in A27 Colombia/67) and in additional sequences present in GenBank (RGN in Asia1 Pak/1/54 and O/Syria/1/87, RSG in A A/IND/110/99, RGE in O KEN/5/95, KGD in O PAK/1/94 and O/IRQ/26/2000, SGD in O Akesu/58, and IGD in O3/Venezuela/51), with substitutions tolerable at any one of the three residues.

Exposure of the 1D N terminus, a region which includes a 10-aa motif which is conserved in several other picornaviruses, is critical for poliovirus entry into the cell (13). The FMDV 1D sequences studied here lack this N-terminal motif, suggesting that some aspects of viral entry may differ from those of other picornaviruses.

Structural implications of capsid protein amino acid conservation. Picornaviral 1A is cleaved from the 1AB (VP0) precursor by unknown mechanisms, and this cleavage is required for virion maturation and infectivity (6, 74). During poliovirus entry into host cells, 1A is involved in receptor-mediated capsid conformational changes resulting in membrane ion channel formation and the release of viral RNA into the cytoplasm (reviewed in reference 42). Several FMDV residues which are known or suspected to affect 1AB cleavage based on their homology to other picornaviral 1ABs were identified in the FMDV isolates examined here. Notably, a number of these residues were located in previously unidentified invariant sequence motifs, suggesting that the critical function of these residues may be contextual. Specifically, these included 1B residues H145, P144, and L83 and 1C residues G39, F41 (replaced by Y in three isolates), and A50, which is included in the conserved motif LDVAEACPT (positions 45 to 53) (8). Additional residues contributing to the 1AB cleavage pocket include 1D P204, contained in the motif RMKRAE(T/L)YCPR (positions 195 to 205), and 1B V32 (I in SAT2 viruses), T33, and Y36, included in the motif TTSTTQSSVG(V/I)T(Y/F)GY (positions 22 to 36). C-terminal to this last motif lies a putative serotype sequence signature (positions 36 to 47) confirmed for 173 available 1B sequences (YSTXEDHXXGPN in A viruses, YXTXEDFVXGPN in O viruses, YATXEDXXGPN in C viruses, YXVXEDAVSGPN in Asia1 viruses, YAXXDXFLPGPN in SAT1 viruses, YADXDSFRXGPN in SAT2 viruses, and YXSADRFLPGPN in SAT3 viruses) and the motif GPNT(S/N)GLEXRVxQAER(F/Y)(F/Y)K (positions 45 to 63), which is present in all FMDV isolates examined here.

An H-rich region at the 1B/1C interphase is proposed to mediate 1B-1C hydrogen bonding and is likely involved in virion sensitivity to acid, a characteristic with implications for virus stability and transmission (2). 1B residues H21 (position 19 in SAT viruses), H145, H157, and H174 were invariant in all FMDV isolates studied here, while H87P and H168Y substitutions were present in serotype C and Euroasiatic viruses, respectively. In addition, 1C contained five invariant H residues, at positions 86 (84 in SAT viruses), 109, 146, 149, and 198 (196 in SAT viruses).

The N termini of five 1C molecules make up the ß annulus at the axis of virion fivefold symmetry, thus contributing to capsid stability (1). The 53 N-terminal residues of 1C, including P4, which is possibly involved in binding to the 3Cpro substrate-binding site, and C7, a residue that is invariant in Euroasiatic isolates and involved in disulfide bonding between 1C N termini, were invariant within serotypes (1). However, all SAT viruses demonstrated an unexpected nonconservative C7V substitution, suggesting that SAT virus capsids may exhibit distinct physical properties.

Nonstructural proteins. (i) 2A. FMDV 2A is an 18-aa peptide similar to the C terminus of cardiovirus 2A, a protein which induces a modification of the cellular translation apparatus resulting in 2A release (19, 20). The FMDV 2A proteins examined here averaged an 89% amino acid identity. Fourteen residues were identified as ≥98% invariant, including the DVEXNPG motif, which is essential for encephalomyocarditis virus 2A activity (21). The high level of 2A conservation likely reflects structural and functional constraints associated with the small size of the protein.

(ii) 2B. FMDV 2B localizes to sites of virus genome replication in ER-derived vesicles (106). The 2B protein in other picornaviruses enhances membrane permeability, blocks protein secretory pathways, suppresses apoptotic responses by affecting intracellular Ca2+ homeostasis, and is implicated in virus-induced cytopathic effects (17, 47, 112).

FMDV 2B, a protein of 154 aa, contains 117 invariant residues (Table 4), with amino acid substitutions limited to only one or two alternate residues per site. A previously undescribed, conserved transmembrane domain was predicted between positions 120 and 140, suggesting that 2B is an integral membrane protein and consistent with the proposed localization of 2B to ER-derived vesicles (106).

(iii) 2C. FMDV 2C is homologous to poliovirus 2C, an ATPase affecting the initiation of minus-strand RNA synthesis and whose precursor, 2BC, induces vesicle formation in the cytoplasm (53). FMDV 2C localizes to membrane-associated virus-replicating complexes (106).

FMDV 2C is 318 aa long and contains 72% invariant residues, including those of the putative ATP/GTP binding domain (positions 110 to 116, 160 to 163, and 243 to 246) (15). Amino acid substitutions unique to individual SAT serotypes were identified between positions 33 and 92.

(iv) 3A. FMDV 3A has been implicated in virus virulence and host range, similar to the 3A proteins of other picornaviruses (33, 61). Deletions in 3A have been associated with FMDV attenuation in cattle and with the porcinophilic phenotype of O Taiwan/97 (54). However, evidence suggests that other viral genetic determinants are necessary for these phenotypes (54, 81).

FMDV 3A ranges from 143 to 153 aa, with insertions/deletions preferentially occurring at positions 70 to 110 and 130 to 150, regions identified here as containing previously undescribed variability. In fact, 3A is one of the most variable proteins encoded by FMDV, with fewer invariant aa (37%) than the variable capsid proteins 1B and 1C and the highly variable Lpro and contrasting with the high degree of conservation previously reported for a limited number of PanAsia 3A sequences (54). As described above, 3A contains residues predicted to undergo positive selection, including the Q44 residue, at which a Q44R mutation was previously associated with a pathogenic phenotype in guinea pigs (78). Although the Q44R mutation was present in several guinea pig-adapted isolates examined here, this mutation was absent from other isolates adapted to guinea pigs and present in still others that were not adapted to guinea pigs, suggesting that alternative mutations may underlie this particular phenotype. A previously described transmembrane domain (positions 60 through 76) (54) which tolerated conservative amino acid substitutions at all positions except invariant residues L64, L68, A70, and I72 was also variable in this study. Our data revealed limited variability in the potential 3A T-cell epitope that was previously identified between positions 21 and 35 (10). The variable nature of 3A suggests that it may be highly informative for epidemiologic or forensic purposes. Additionally, the likely role for 3A in virulence and host range suggests that interactions with host factors underlie 3A's variability and the diversifying selection predicted to act upon it.

(v) 3B. FMDV is the only picornavirus to encode multiple 3B proteins (3B1, 23 aa; 3B2, 24 aa; and 3B3, 24 aa), whose homologue in poliovirus primes genomic RNA synthesis during virus replication (26, 85). In our study, 3B1, 3B2, and 3B3 were present in all FMDV isolates. The motif GPYXGP (except for the substitution GPYXRP in Sat 1/20 Rv 11/37), which contains a Y residue homologous to the poliovirus Y3 residue involved in phosphodiester linkage of 3B to the 5' end of the viral genome, was invariant in the N terminus of each protein (5) (Table 7). The 3B C-terminal regions contained more amino acid variability, including the majority of observed nonconservative substitutions and the fewest invariant residues. Notably, 3B3 was highly conserved in all isolates examined, supporting previous experimental evidence indicating that only this 3B isoform is essential for virus viability (26, 84). 3B1 and 3B2 were more variable, and in fact, contained residues predicted here to undergo diversifying selection (3B1 residues 4 and 11 and 3B2 residues 17, 18, and 19). These data, in conjunction with experimental data indicating that the deletion of 3B1 and 3B2 may affect FMDV virulence and host range (84), suggest that, similar to that of 3A, 3B1 and 3B2 protein variability reflects host range-specific functions.


View this table:
[in this window]
[in a new window]
  		TABLE 7. FMDV 3B (VPg) variability

 
(vi) 3C. FMDV 3Cpro is a 213-amino-acid proteinase involved in the cleavage of most viral precursor proteins and of host factors associated with translation (eIF4A and eIF4G) and transcription (histone H3). A previous mutational analysis of 3Cpro identified the catalytic triad residues C163, H46, and D84 and residues critical for proper protein folding, i.e., D84 and Y136, all of which were invariant in the isolates examined here (34). A previously undescribed invariant motif (VKGQDMLSDAALMVLH) and a predicted transmembrane domain were identified at positions 76 to 91 and 27 to 44, respectively. A high degree of overall 3Cpro conservation indicates a limited tolerance for mutations, likely due to significant structural/functional constraints. Despite the 3Cpro conservation, 36 SAT-specific residues were identified, suggesting that distinct selective pressures exist for the two virus groups.

(vii) 3D. The 469-amino-acid viral RNA-dependent RNA polymerase 3Dpol, responsible for generating minus- and plus-sense genomic RNA, is one of the most conserved FMDV proteins. Our analysis indicated that although it is conserved, 3Dpol is more tolerant of substitutions than was previously reported (56), as our results extended the proportion of variable residues from 8.6% to 26%. D245, N307, and G295, which are essential for maintaining the functional integrity of the picornaviral polymerases (56), were invariant in FMDV 3Dpol, along with other residues described as being highly conserved among all RNA-dependent RNA polymerases (31), including the NTP-binding residues G337, D338, and D339 (115). The previously described picornaviral polymerase peptide motifs KDELR, PSG, YGGD, and FLKR were conserved among all FMDV isolates analyzed (18, 51). Finally, the three previously described hypervariable, hydrophobic antigenic regions of 3Dpol (aa 1 to 12, 64 to 76, and 143 to 153) were also variable in all virus isolates examined here (31).

3' UTR. The picornaviral 3' UTR binds several viral and host proteins and is believed to contain structural cis-acting elements required for negative-strand RNA synthesis (4). Removal of the terminal poly(A) tract or mutagenesis of structural elements abrogates the infectivity of FMDV infectious clones (98). FMDV 3' UTR sequences stimulate virus-specific, IRES-dependent translation (67) and likely affect other aspects of viral replication, including genome circularization (37).

The FMDV 3' UTRs were highly variable in length among the isolates, ranging from 85 to 101 nt. A secondary structure analysis of the FMDV 3' UTR confirmed the Y shape which is also predicted for other picornaviral 3' UTRs, suggesting that the structure plays an important role in the 3' UTR function (Fig. 1D and data not shown). A previously unidentified motif was located at the vertex of the Y structure between positions 37 and 61 (Fig. 1C and D). In some cases, however, structural features which could conceivably affect the efficiency of ribosome/RNA dissociation and translation were observed (for example, the exceptionally long stem of A10 Holland/42 and the short stem of O2 Brescia/47) (data not shown).

FMDV phylogenetic and recombination analysis. To date, phylogenetic analyses have been performed largely on FMDV sequences from the 1D coding region. These analyses have permitted the discrimination among serotypically related FMDV isolates (55). Split decomposition, a clustering technique used previously on 1D sequences to suggest the occurrence of quasispecies evolution in Euroasiatic FMDV (22), was used to examine the complete genome sequences of all 103 isolates described here (Fig. 4A). The results indicated that complex phylogenetic relationships exist between members of different serotypes, including relationships between the A24 Argentina/65 and European/South American O1 viruses; between A12 Valle strain 119 and C Waldman strain 149; between the O1 M11, O2 Brescia/47, and O3 Venezuela/71 viruses and European serotype A viruses; between 05India/62 and Asia1 viruses from Israel; among six SAT1 and three SAT3 viruses; and most notably, between SAT1/7 Isrl 4/62 and SAT2/3 Kenya 11/60 and other SAT1 and SAT2 viruses (Fig. 4A and data not shown).



View larger version (54K):
[in this window]
[in a new window]
  		FIG. 4. FMDV phylogenetic and recombination analysis. (A) One hundred three FMDV genomes were aligned with ClustalW and visually screened for reticulated branching patterns by split decomposition cluster analysis, using SplitsTree with Hamming distances. Numbers at terminal nodes correspond to the isolates presented in Table 1, arcs and labels indicate serotype-specific groups, similar colors indicate isolates exhibiting obvious interserotypic relationships, red stars indicate SAT viruses with unique nonstructural protein sequences, and bars represent estimated distances. Similar results were obtained with alternative corrections for multiple substitution (data not shown). Subgenomic split decomposition analysis was conducted on structural protein (1A to 1D) (B) and nonstructural protein (Lpro and 2A to 3D) (C) coding regions aligned with Dialign. (D to G) Similarity (upper box) and bootscanning (lower box) plots were generated with a query isolate against reference isolates aligned with Dialign, indicating corrected percent similarities and bootstrap frequencies as percentages of permuted trees, respectively. Polyprotein schematics are shown above the similarity plots. Query isolates are indicated at the left within similarity plots and included the following: (D) SAT 1/7 Isrl 4/62 (no. 96), (E) A24 Argentina/65 (no. 3), (F) O1 M11 (no. 74), and (G) A10 Holland/42 (no. 25). Reference isolates are indicated individually or as groups in color, as labeled at the right within similarity plots, and included the following: (D) red, SAT 2/3 Kenya 11/60 (no. 99); green, SAT1 viruses (no. 88 to 95); blue, SAT2 and SAT3 viruses (no. 97 and 98 and 100 to 103, respectively); black, A viruses (no. 30, 38, and 46), O viruses (no. 80 and 86), C viruses (no. 57), and Asia1 viruses (no. 51); (E) green, A14 Spain/59 (no. 36); orange, A29 Peru/69 (no. 22); blue, A viruses (no. 2, 10, 18, 19, 22, 25, 27, 30, 38, 39, 42, 46, and 47); red, O1 Campos/58 (no. 70); gray, O viruses (no. 74, 78, 80, 83, and 87); (F) red, A10 Holland/42 (no. 25); pink, A1 Bayern/Bavaria (no. 27); green, A4 WG/72 (no. 30); blue, A viruses (no. 2, 10, 18, 19, 22, 38, 39, 42, 46, 47); black O1 Campos/58 (no. 70); orange, O1 Vallee/39 (no. 77); gray, O viruses (no. 63, 75, 77, 78, 80, 83, 87); (G) red, A12 Valle Strain 119 (no. 26); pink, A1 Bayern/Bavaria (no. 27); orange, A3 Mecklenburg/68 (no. 28); green, A4 WG/72 (no. 30); blue, A viruses (no. 2, 3, 5, 10, 13, 17, 19, 21, 22, 29, 32, 33, 36 to 39, 41, 46, and 47). Aligned polyprotein-encoding nt positions are indicated below the bootscan plots. Colored bars within similarity plot boxes indicate fragments identified by GENECONV as having significant global permutation values (P < 0.05) shared between the query sequence and particular reference sequences.

 
The data from analyses of individual protein-encoding regions indicated that incongruent tree topologies exist for different genomic regions (data not shown). An analysis of 1B, 1C, 1D, or all of the structural protein-encoding regions together indicated few obvious phylogenetic incongruities, placing viruses into serotypic, temporal, and/or geographic groups (Fig. 4B and data not shown). In general, analyses of the entire structural protein-encoding region improved bootstrap/confidence values relative to 1D analysis alone (data not shown). Conversely, analyses of the 5' UTR and regions encoding individual or all nonstructural proteins failed to produce the serotype-inclusive groups observed with capsid protein sequences, often yielding unresolved, star-like topologies among diverse serotypes (Fig. 4C and data not shown). Consistent with the complex relationships illustrated in Fig. 4A, analyses of nonstructural protein-encoding regions indicated groupings of serotypically disparate viruses, which extended previous but limited data (54, 110, 113) (Fig. 4C). Similarity plots and bootscanning analysis further indicated that several of these interserotypically related viruses contained nonstructural protein-encoding regions most similar to those of viruses with divergent capsid protein-encoding sequences, suggesting the occurrence of intertypic recombination during the history of these lineages (Fig. 4D to F and data not shown). Examples of potential intratypic recombination were also noted, more obviously in nonstructural protein-encoding regions but also in structural protein-encoding regions, as has been suggested previously (Fig. 4F and G and data not shown) (35, 109).

Taken together, these data suggest that FMDV capsid protein sequences may undergo intertypic recombination infrequently relative to that in other genomic regions, which conceivably undergo complex recombination events and fail to display serotype-specific phylogenetic relationships. These observations are consistent with and extend previous reports of FMDV recombination, which has been inferred from incongruent Lpro, 3Cpro, and 1D topologies and observed (and suggested to predominate) in C-terminal genomic regions (52, 57, 101, 113). Similar observations have recently been made for enterovirus genomes, which were suggested to undergo relatively extensive recombination in nonstructural gene regions while generally maintaining serospecific monophyly (79). These observations, including the inability to reliably define viral relationships based on nonstructural protein-encoding sequences and data suggestive of recombination, raise interesting questions about FMDV genome evolution in nature and the relative contribution of recombination to the generation of FMDV genetic and population diversity.

Notably, the analysis described here identified a novel SAT lineage represented by SAT1-7 Isrl 4/62 and SAT2-3 Kenya 11/60, which contain nonstructural protein-encoding regions that are divergent from each other but are also clearly distinct from those of the other SAT and Euroasiatic lineages presented here (Fig. 4C). The distinct nature of nonstructural proteins from these structurally/serotypically SAT-like viruses was supported by a previously available genome sequence similar to that of SAT2-3 Kenya (SAT2 Ken/3/57; 97% overall nucleotide identity) (Table 1) and by previous data indicating that the SAT2 Ken/3/57 Lpro- and 3C-encoding regions may be distinct from those of other SATs (102). Taken together, these results suggest that more FMDV genome diversity may exist in nature than is currently indicated by serology or 1D sequence analysis.


   ACKNOWLEDGMENTS
 
We thank J. Lubroth for providing FMDV isolates and A. Lakowitz and A. Zsak for providing excellent technical assistance.


   FOOTNOTES
 
* Corresponding author. Mailing address: Department of Pathobiology and Veterinary Science, University of Connecticut, 61 N. Eagleville Road, Unit-3089, Storrs, CT 06269-3089. Phone: (860) 486-2845. Fax: (860) 486-2794. E-mail: drock{at}cshore.com. Back

{dagger} Supplemental material for this article may be found at http://jvi.asm.org/. Back


   REFERENCES
Top
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 References
 

  1. Acharya, R., E. Fry, D. Stuart, G. Fox, D. Rowlands, and F. Brown. 1990. The structure of foot-and-mouth disease virus: implications for its physical and biological properties. Vet. Microbiol. 23:21-34.[CrossRef][Medline]
  2. Acharya, R., E. Fry, D. Stuart, G. Fox, D. Rowlands, and F. Brown. 1989. The three-dimensional structure of foot-and-mouth disease virus at 2.9 Å resolution. Nature (London) 337:709-716.[CrossRef][Medline]
  3. Afonso, C. L., E. R. Tulman, Z. Lu, E. Oma, G. F. Kutish, and D. L. Rock. 1999. The genome of Melanoplus sanguinipes entomopoxvirus. J. Virol. 73:533-552.[Abstract/Free Full Text]
  4. Agol, V. I., A. V. Paul, and E. Wimmer. 1999. Paradoxes of the replication of picornaviral genomes. Virus Res. 62:129-147.[CrossRef][Medline]
  5. Ambros, V., and D. Baltimore. 1978. Protein is linked to the 5' end of poliovirus RNA by a phosphodiester linkage to tyrosine. J. Biol. Chem. 253:5263-5266.[Abstract/Free Full Text]
  6. Ansardi, D. C., D. C. Porter, and C. D. Morrow. 1992. Myristylation of poliovirus capsid precursor P1 is required for assembly of subviral particles. J. Virol. 66:4556-4563.[Abstract/Free Full Text]
  7. Bandelt, H. J., and A. W. Dress. 1992. Split decomposition: a new and useful approach to phylogenetic analysis of distance data. Mol. Phylogenet. Evol. 1:242-252.[CrossRef][Medline]
  8. Basavappa, R., R. Syed, O. Flore, J. P. Icenogle, and D. J. Filman. 1994. Role and mechanism of the maturation cleavage of VP0 in poliovirus assembly: structure of the empty capsid assembly intermediate at 2.9 A resolution. Protein Sci. 3:1651-1669.[Abstract]
  9. Belsham, G. J. 1992. Dual initiation sites of protein synthesis on foot-and-mouth disease virus RNA are selected following internal site and scanning of ribosomes in vivo. EMBO J. 11:1105-1110.[Medline]
  10. Blanco, E., M. Garcia-Briones, A. Sanz-Parra, P. Gomes, E. De Oliveira, M. L. Valero, D. Andreu, V. Ley, and F. Sobrino. 2001. Identification of T-cell epitopes in nonstructural proteins of foot-and-mouth disease virus. J. Virol. 75:3164-3174.[Abstract/Free Full Text]
  11. Brendel, V., P. Bucher, I. Nourbakhsh, B. E. Blaisdell, and S. Karlin. 1992. Methods and algorithms for statistical analysis of protein sequences. Proc. Natl. Acad. Sci. USA 89:2002-2006.[Abstract/Free Full Text]
  12. Clarke, B. E., A. L. Brown, K. M. Currey, S. E. Newton, D. J. Rowlands, and A. R. Carroll. 1987. Potential secondary and tertiary structure in the genomic RNA of foot-and-mouth disease virus. Nucleic Acids Res. 15:7067-7079.[Abstract/Free Full Text]
  13. Curry, S., M. Chow, and J. M. Hogle. 1996. The poliovirus 135S particle is infectious. J. Virol. 70:7125-7131.[Abstract/Free Full Text]
  14. Curry, S., E. Fry, W. Blakemore, R. Abu-Ghazaleh, T. Jackson, A. M. King, S. Lea, J. Newman, and D. Stuart. 1997. Dissecting the roles of VP0 cleavage and RNA packaging in picornavirus capsid stabilization: the structure of empty capsids of foot-and-mouth disease virus. J. Virol. 71:9743-9752.[Abstract]
  15. Dever, T. E., M. J. Glynias, and W. C. Merrick. 1987. GTP-binding domain: three consensus sequence elements with distinct spacing. Proc. Natl. Acad. Sci. USA 84:1814-1818.[Abstract/Free Full Text]
  16. Devereux, J., P. Haeberli, and O. Smithies. 1984. A comprehensive set of sequence analysis programs for the VAX. Nucleic Acids Res. 12:387-395.
  17. Doedens, J. R., and K. Kirkegaard. 1995. Inhibition of cellular protein secretion by poliovirus proteins 2B and 3A. EMBO J. 14:894-907.[Medline]
  18. Doherty, M., D. Todd, N. McFerran, and E. M. Hoey. 1999. Sequence analysis of a porcine enterovirus serotype 1 isolate: relationships with other picornaviruses. J. Gen. Virol. 80:1929-1941.[Abstract/Free Full Text]
  19. Donnelly, M., L. Hughes, G. Luke, H. Mendoza, E. ten Dam, D. Gani, and M. Ryan. 2001. The "cleavage" activities of foot-and-mouth disease virus 2A site-directed mutants and naturally occurring "2A-like" sequences. J. Gen. Virol. 82:1027-1041.[Abstract/Free Full Text]
  20. Donnelly, M., G. Luke, A. Mehrotra, X. Li, L. Hughes, D. Gani, and M. Ryan. 2001. Analysis of the aphthovirus 2A/2B polyprotein "cleavage" mechanism indicates not a proteolytic reaction, but a novel translational effect: a putative ribosomal "skip." J. Gen. Virol. 82:1013-1025.[Abstract/Free Full Text]
  21. Donnelly, M. L., D. Gani, M. Flint, S. Monaghan, and M. D. Ryan. 1997. The cleavage activities of aphthovirus and cardiovirus 2A proteins. J. Gen. Virol. 78:13-21.[Abstract]
  22. Dopazo, J., A. Dress, and A. V. Haeseler. 1993. Split decomposition: a technique to analyze viral evolution. Proc. Natl. Acad. Sci. USA 90:10320-10324.[Abstract/Free Full Text]
  23. Duke, G. M., J. E. Osorio, and A. C. Palmenberg. 1990. Attenuation of Mengo virus through genetic engineering of the 5'-noncoding poly(C) tract. Nature 343:474-476.[CrossRef][Medline]
  24. Escarmis, C., M. Toja, M. Medina, and E. Domingo. 1992. Modification of the 5' untranslated region of foot-and-mouth disease virus after prolonged persistence in cell culture. Virus Res. 26:113-125.[CrossRef][Medline]
  25. Ewing, B., L. Hillier, M. C. Wendl, and P. Green. 1998. Base-calling of automated sequencer traces using Phred. I. Accuracy assessment. Genome Res. 8:175-185.[Abstract/Free Full Text]
  26. Falk, M., F. Sobrino, and E. Beck. 1992. VPg gene amplification correlates with infective particle formation in foot-and-mouth disease virus. J. Virol. 66:2251-2260.[Abstract/Free Full Text]
  27. Felsenstein, J. 1993. PHYLIP (Phylogeny Inference Package), version 3.5c. Distributed by the author. Department of Genetics, University of Washington, Seattle.
  28. Galtier, N., M. Gouy, and C. Gautier. 1996. SeaView and Phylo_win, two graphic tools for sequence alignment and molecular phylogeny. Comput. Applic. Biosci. 12:543-548.[Abstract/Free Full Text]
  29. Garnier, J., D. J. Osguthorpe, and B. Robson. 1978. Analysis of the accuracy and implications of simple methods for predicting the secondary structure of globular proteins. J. Mol. Biol. 120:97-120.[CrossRef][Medline]
  30. George, M., R. Venkataramanan, C. B. Gurumurthy, and D. Hemadri. 2001. The non-structural leader protein gene of foot-and-mouth disease virus is highly variable between serotypes. Virus Genes 22:271-278.[CrossRef][Medline]
  31. George, M., R. Venkataramanan, B. Pattnaik, A. Sanyal, C. B. Gurumurthy, D. Hemadri, and C. Tosh. 2001. Sequence analysis of the RNA polymerase gene of foot-and-mouth disease virus serotype Asia1. Virus Genes 22:21-26.[CrossRef][Medline]
  32. Gorbalenya, A. E., E. V. Koonin, and M. M. Lai. 1991. Putative papain-related thiol proteases of positive strand RNA viruses. Identification of rubi- and aphthovirus proteases and delineation of a novel conserved domain associated with proteases of rubi-, alpha-, and coronaviruses. FEBS Lett. 288:201-205.[CrossRef][Medline]
  33. Graff, J., C. Kasang, A. Normann, M. Pfisterer-Hunt, S. M. Feinstone, and B. Flehmig. 1994. Mutational events in consecutive passages of hepatitis A virus strain GBM during cell culture adaptation. Virology 204:60-68.[CrossRef][Medline]
  34. Grubman, M. J., M. Zellner, G. Bablanian, P. W. Mason, and M. E. Piccone. 1995. Identification of the active-site residues of the 3C proteinase of foot-and-mouth disease virus. Virology 213:581-589.[CrossRef][Medline]
  35. Haydon, D. T., A. D. S. Bastos, and P. Awadalla. 2004. Low linkage disequilibrium indicative of recombination in foot-and-mouth disease virus gene sequence alignments. J. Gen. Virol. 85:1095-1100.[Abstract/Free Full Text]
  36. Henikoff, S., and J. G. Henikoff. 1994. Protein family classification based on searching a database of blocks. Genomics 19:97-107.[CrossRef][Medline]
  37. Herold, J., and R. Andino. 2001. Poliovirus RNA replication requires genome circularization through a protein-protein bridge. Mol. Cell 7:581-591.[CrossRef][Medline]
  38. Hinton, T., N. Ross-Smith, S. Warner, G. Belsham, and B. Crabb. 2002. Conservation of L and 3C proteinase activities across distantly related aphthoviruses. J. Gen. Virol. 83:3111-3121.[Abstract/Free Full Text]
  39. Hinton, T. M., F. Li, and B. S. Crabb. 2000. Internal ribosomal entry site-mediated translation initiation in equine rhinitis A virus: similarities to and differences from that of foot-and-mouth disease virus. J. Virol. 74:11708-11716.[Abstract/Free Full Text]
  40. Hofacker, I. L. 2003. Vienna RNA secondary structure server. Nucleic Acids Res. 31:3429-3431.[Abstract/Free Full Text]
  41. Hofmann, K., and W. Stoffel. 1993. TMBase, a database of membrane spanning protein segments. Biol. Chem. Hoppe-Seyler 347:166.
  42. Hogle, J. M. 2002. Poliovirus cell entry: common structural themes in viral cell entry pathways. Annu. Rev. Microbiol. 56:677-702.[CrossRef][Medline]
  43. Huang, X., and A. Madan. 1999. CAP3: a DNA sequence assembly program. Genome Res. 9:868-877.[Abstract/Free Full Text]
  44. Huson, D. 1998. Splitstree: a program for analyzing and visualizing evolutionary data. Bioinformatics 14:68-73.[Abstract/Free Full Text]
  45. Jackson, T., A. M. Q. King, D. I. Stuart, and E. Fry. 2003. Structure and receptor binding. Virus Res. 91:33-46.[CrossRef][Medline]
  46. Jaeger, J. A., D. H. Turner, and M. Zuker. 1990. Predicting optimal and suboptimal secondary structure for RNA. Methods Enzymol. 183:281-306.[Medline]
  47. Jecht, M., C. Probst, and V. Gauss-Muller. 1998. Membrane permeability induced by hepatitis A virus proteins 2B and 2BC and proteolytic processing of HAV 2BC. Virology 252:218-227.[CrossRef][Medline]
  48. Jonassen, I. 1997. Efficient discovery of conserved patterns using a pattern graph. Comput. Appl. Biol. Sci. 13:509-522.[Abstract/Free Full Text]
  49. Jonassen, I., J. Collins, and D. Higgins. 1995. Finding flexible patterns in unaligned protein sequences. Protein Sci. 4:1587-1595.[Abstract]
  50. Jones, D. T., W. R. Taylor, and J. M. Thornton. 1994. A model recognition approach to the prediction of all-helical membrane protein structure and topology. Biochemistry 33:3038-3049.[CrossRef][Medline]
  51. Kaku, Y., S. Yamada, and Y. Murakami. 1999. Sequence determination and phylogenetic analysis of RNA-dependent RNA polymerase (RdRp) of the porcine enterovirus 1 (PEV-1) Talfan strain. Arch. Virol. 144:1845-1852.[CrossRef][Medline]
  52. King, A. M. Q., D. McCahon, K. Saunders, J. W. I. Newman, and W. R. Slade. 1985. Multiple sites of recombination within the RNA genome of foot-and-mouth disease virus. Virus Res. 3:373-384.[CrossRef][Medline]
  53. Klein, M., H. J. Eggers, and B. Nelsen-Salz. 1999. Echovirus 9 strain barty non-structural protein 2C has NTPase activity. Virus Res. 65:155-160.[CrossRef][Medline]
  54. Knowles, N. J., P. R. Davies, T. Henry, L. V. O'Donnel, J. M. Pacheco, and P. W. Mason. 2001. Emergence in Asia of foot-and-mouth disease viruses with altered host range: characterization of alterations in the 3A protein. J. Virol. 75:1551-1556.[Abstract/Free Full Text]
  55. Knowles, N. J., and A. R. Samuel. 2003. Molecular epidemiology of foot-and-mouth disease virus. Virus Res. 91:65-80.[CrossRef][Medline]
  56. Koonin, E. V. 1991. The phylogeny of RNA-dependent RNA polymerases of positive-strand RNA viruses. J. Gen. Virol. 72:2197-2206.[Abstract/Free Full Text]
  57. Krebs, O., and O. Marquardt. 1992. Identification and characterization of foot-and-mouth disease virus O1 Burgwedel/1987 as an intertypic recombinant. J. Gen. Virol. 73:613-619.[Abstract/Free Full Text]
  58. Krogh, A., M. Brown, I. S. Mian, K. Sjolander, and D. Haussler. 1994. Hidden Markov models in computational biology: applications to protein modeling. J. Mol. Biol. 235:1501-1531.[CrossRef][Medline]
  59. Kronovetr, J., and T. Skern. 2002. Foot-and-mouth disease virus leader proteinase: a papain-like enzyme requiring an acidic environment in the active site. FEBS Lett. 528:58-62.[CrossRef][Medline]
  60. Kuhn, R., N. Luz, and E. Beck. 1990. Functional analysis of the internal translation initiation site of foot-and-mouth disease virus. J. Virol. 64:4625-4631.[Abstract/Free Full Text]
  61. Lama, J., M. A. Sanz, and L. Carrasco. 1998. Genetic analysis of poliovirus protein 3A: characterization of a non-cytopathic mutant virus defective in killing Vero cells. J. Gen. Virol. 79:1911-1921.[Abstract]
  62. Logan, D., R. Abu-Ghazaleh, W. Blakemore, S. Curry, T. Jackson, A. King, S. Lea, R. Lewis, J. Newman, N. Parry, et al. 1993. Structure of a major immunogenic site on foot-and-mouth disease virus. Nature 362:566-568.[CrossRef][Medline]