foot and mouth disease ( fmd ) is an economically important disease of cloven - hoofed livestock such as cattle , sheep , and other domestic animals in addition to several wild - life species .
the causative agent of the fmd is the foot and mouth disease virus ( fmdv ) .
the disease is considered one of the most important barriers to the worldwide trade of livestock and animal products .
fmdv isolates sampled from the world have been categorized into seven different serotypes , a , o , c , asia1 , and south african territories 1 ( sat1 ) , sat2 , and sat3 , according to antigenic levels ; and a large number of variants have appeared within each serotype . out of the three serotypes ( o , a , and asia1 ) found to be prevalent in iran in recent years , serotype a has been found to be more antigenically and genetically diverse than the others .
the fmdv possesses a single - stranded rna molecule consisting of about 8,200 nucleotides within an icosahedral capsid composed of structural proteins .
the open reading frame encodes a single polyprotein which can be cleaved into four structural proteins ( vp4 , vp2 , vp3 , and vp1 ) and eight non - structural proteins ( l , 2a , 2b , 2c , 3a , 3b , 3c , and 3d polymerase [ 3d ] ) . in general ,
mutations or deletions in structural proteins may help the fmdv to evade immune responses of the host whereas mutations or deletions in non - structural proteins can inhibit viral replication and protein processing .
the fmdv usually infects cells by binding to integrin receptors via a long flexible loop ( g - h loop ) of vp1 ( 1d ) .
the sequence of this loop contains a conserved arginineglycine - aspartic acid ( rgd ) tripeptide motif which is characteristic of ligands that bind to integrin receptors .
the vp1 protein is encoded by the vp1 coding region of viral rna ; this coding region is 627~639 bp long and produces a protein containing 209~213 amino acid residues depending on the serotype .
this gene contains a taa stop codon and is responsible for encoding a protein containing 470 amino acid residues .
similar to other picornaviruses , the fmdv 3d protein is a viral - encoded rna polymerase . among the different fmdv serotypes and subtypes , both the nucleotide and amino acid sequences of 3d
uninterrupted co - circulation of the virus in the environment and the absence of viral polymerase proofreading activity during the replication process lead to the appearance of various genetic isolates .
analysis of the viral genome sequence is a crucially important approach for monitoring field isolates in areas where the disease is endemic .
previous studies have shown that iran has one of the highest reported rates of fmd cases per year . in iran ,
fmd is largely controlled by vaccinating cattle and sheep with vaccines prepared against isolates that are likely to be encountered in the region .
spread of the disease is promoted by the presence of large populations of susceptible animals , low vaccination rates , prevalence of multiple serotypes ( including serotypes a , asia1 , and o ) , and unlimited movement of susceptible animals in the country .
furthermore , there is endemic co - circulation of multiple genotypes of type a virus .
the high incidence of fmd has allowed the identification of new variants of the virus over the last 7 years .
the general objective of this study was to determine the nucleotide sequences of genes encoding vp1 and 3d proteins of a type a fmdv ( a / iran87 ) highly adapted to cell culturing , and to compare them to other corresponding sequences available in the genbank database .
the razi vaccine and serum research institute ( iran ) generously provided samples of the fmdv used for our study
. clinical specimens of fmdv which included samples in tongue epithelium tissue were collected in 1987 from infected calves displaying clinical symptoms of fmd in an iranian field located in tehran . the a / iran87 strain was initially detected and isolated nearly more than two decades ago .
thus , the highly passaged virus investigated in this study was not virus isolated from the field .
the tissue samples were used to infect baby hamster kidney 21 ( bhk-21 ) cells .
the virulent isolate was passaged in ba ( cell line derived from pig kidney ) and bhk-21 cell monolayers .
propagation steps included six passages in ba cells and then four passages in bhk-21 cells .
thereafter , it has been used for the development of fmd vaccines during the past years as a vaccine seed .
the infected cell culture supernatant from the ~150th passage was clarified and stored at -70 before use .
total rna was extracted from fmdv - infected bhk-21 cells using an rnx plus kit ( cinnagen , iran ) according to the manufacturer 's instructions .
briefly , the cell culture was centrifuged at 85,000 g at 4 for 2 h. the pellet was resuspended in 200 l of pbs , mixed well with 400 l of rnx reagent , and incubated at room temperature for 5 min .
the pellet was washed with 75% ethanol and dissolved in 20 l of rnase - free water .
rt - pcr was carried out in 50 l of reaction mixture containing 10 l of 5 reaction buffer , 4 l of mixed dntps ( 2.5 mm each ) , 1 l of amv enzyme ( titan one tube rt - pcr system kit ; roche diagnostic , germany ) , 1 l of each primer ( 10 pmol each ) , 4 l of rna template , 2.5 l dtt , 3 l 25 mm mgcl2 , and 23.5 l of h2o . the following conditions were used for amplification : 42 for 30 min , 94 for 3 min ; 30 cycles of 94 for 30 sec , 52 for 30 sec , 72 for 40 sec ; then followed by 72 for 5 min .
the vp1 and 3d coding regions of fmdv ( 621 and 690 bp , respectively ) were amplified using standard methods with a one - step rt - pcr system and specific primer combinations .
5'-gacatgtcctggtgcatctg-3 ' for 3d , and forward : fmg13 5'-accaggatgatgattggcag-3 ' and reverse : fm15 5'-tttcactcctacggtgtcgc-3 ' for vp1 .
amplified pcr products of the expected length were separated by electrophoresis in a 1% agarose gel , stained with ethidium bromide , and visualized under a uv transilluminator .
after successful amplification of the target dna sequences , fragments were purified using a gel extraction kit ( roche , germany ) following the recommendations of the supplier .
ligation was performed with plasmid vector ptz57r / t ( fermentas , germany ) in 0.165 g , 0.18 pmol ends and 0.54 pmol ends purified pcr fragment in 1 ligation buffer , polyethylene glycol ( peg 4000 ) , 5 units of t4 dna ligase , and up to 30 l of deionized water at 22 for 16 h. the ligated products were used to transform chemically competent cells ( xl1-blue cells ) and white colonies grown on lb plates were randomly selected .
the cloned pcr products were purified using plasmid purification kit ( roche diagnostics , germany ) according to the manufacturer 's instructions and sequencing was carried out in both directions using a t7 promoter primer ( mwg biotech , germany ) .
the published sequences of 66 fmdv type a isolates recovered from different parts of the world were included in this analysis and compared to the corresponding sequence of the a / iran87 isolate .
reference fmdv sequences were obtained from the national center for biotechnology information ( ncbi , usa ) .
the sequences used for each gene were first examined to exclude ambiguous sequences which were incomplete or frame - shifted . to determine the degree to which genetic diversity
was observed in the vp1 and 3d proteins , multiple alignments and comparisons of the predicted amino acid sequences of isolates were carried out ( fig .
nucleotide sequence homology / divergence was calculated using the megalign project of the dnastar software package ( ver 5.1 ; dnastar , usa ; data not shown ) .
the histories of the type a fmdv field isolates , including year of isolation , accession number , and geographical distribution , are presented in table 1 .
the nucleotide sequences of the vp1 and 3d coding regions of a / iran87 have been submitted to genbank ( accession no .
total rna was extracted from fmdv - infected bhk-21 cells using an rnx plus kit ( cinnagen , iran ) according to the manufacturer 's instructions .
briefly , the cell culture was centrifuged at 85,000 g at 4 for 2 h. the pellet was resuspended in 200 l of pbs , mixed well with 400 l of rnx reagent , and incubated at room temperature for 5 min .
the pellet was washed with 75% ethanol and dissolved in 20 l of rnase - free water .
rt - pcr was carried out in 50 l of reaction mixture containing 10 l of 5 reaction buffer , 4 l of mixed dntps ( 2.5 mm each ) , 1 l of amv enzyme ( titan one tube rt - pcr system kit ; roche diagnostic , germany ) , 1 l of each primer ( 10 pmol each ) , 4 l of rna template , 2.5 l dtt , 3 l 25 mm mgcl2 , and 23.5 l of h2o .
the following conditions were used for amplification : 42 for 30 min , 94 for 3 min ; 30 cycles of 94 for 30 sec , 52 for 30 sec , 72 for 40 sec ; then followed by 72 for 5 min .
the vp1 and 3d coding regions of fmdv ( 621 and 690 bp , respectively ) were amplified using standard methods with a one - step rt - pcr system and specific primer combinations .
amplified pcr products of the expected length were separated by electrophoresis in a 1% agarose gel , stained with ethidium bromide , and visualized under a uv transilluminator .
after successful amplification of the target dna sequences , fragments were purified using a gel extraction kit ( roche , germany ) following the recommendations of the supplier .
ligation was performed with plasmid vector ptz57r / t ( fermentas , germany ) in 0.165 g , 0.18 pmol ends and 0.54 pmol ends purified pcr fragment in 1 ligation buffer , polyethylene glycol ( peg 4000 ) , 5 units of t4 dna ligase , and up to 30 l of deionized water at 22 for 16 h. the ligated products were used to transform chemically competent cells ( xl1-blue cells ) and white colonies grown on lb plates were randomly selected .
the cloned pcr products were purified using plasmid purification kit ( roche diagnostics , germany ) according to the manufacturer 's instructions and sequencing was carried out in both directions using a t7 promoter primer ( mwg biotech , germany ) .
the published sequences of 66 fmdv type a isolates recovered from different parts of the world were included in this analysis and compared to the corresponding sequence of the a / iran87 isolate .
reference fmdv sequences were obtained from the national center for biotechnology information ( ncbi , usa ) .
the sequences used for each gene were first examined to exclude ambiguous sequences which were incomplete or frame - shifted . to determine the degree to which genetic diversity
was observed in the vp1 and 3d proteins , multiple alignments and comparisons of the predicted amino acid sequences of isolates were carried out ( fig .
nucleotide sequence homology / divergence was calculated using the megalign project of the dnastar software package ( ver 5.1 ; dnastar , usa ; data not shown ) .
the histories of the type a fmdv field isolates , including year of isolation , accession number , and geographical distribution , are presented in table 1 .
the nucleotide sequences of the vp1 and 3d coding regions of a / iran87 have been submitted to genbank ( accession no .
genetic comparison of the field isolates to the vaccine strain is of significant importance for testing the usefulness of the existing vaccine strain as well as selecting new vaccine strain(s ) . during virus circulation ,
amino acid substitutions accumulate at the different positions of the protein sequence . in the present study ,
a highly passaged field isolate of fmdv ( a / iran87 ) named mardabad or a87 isolate was selected for genetic analysis .
sequences of the protein coding regions ( vp1 and 3d ) from the high - passage cell - adapted vaccine strain were determined by pcr amplification and sequencing .
comparison of the amino acid sequences revealed amino acid changes in both structural and non - structural proteins .
the number of sequence differences exhibited by each of the isolates showed that a / iran87 contains four amino acid substitutions at positions 17 , 26 , 50 , and 57 in the 3d coding region ( fig .
1 ) . there is clear evidence of novel amino acid substitutions ( alaser , aspglu , gluasn , and alagly ) in two domains of the 3d protein .
comparison of the vp1 protein among different variants showed that a - iran - vaccinal contained a change at position 179 .
1a shows that amino acid changes among the vp1 genes of the field isolates consist of numerous nucleotide substitutions and deletions compared to the consensus sequence .
comparison of the predicted amino acid sequences of vp1 region in the field isolates revealed amino acid deletions in three iranian vaccine isolates ( a - iran - vaccinal , a iran04 , and a / iran87 ) . as shown in fig .
1a , deletions occurred within the 13 amino acid positions ( 168 to 180 ) of the vp1 region .
three - dimensional analysis of 3d protein showed that the aspglu substitution occurred in a beta sheet located in a small groove of the protein .
the variations were not distributed uniformly along the genes ; there were areas of high and low incidence of nucleotide sequence variation ( data not shown ) .
the region of the vp1 gene between amino acids 157 and 183 contained a high degree of variation and encoded an important immunogenic site on the viral surface .
consequently , nucleotide changes in this region are most likely involved in the appearance of new antigenic variants .
the vp1 and 3d nucleotide sequences of the selected fmdv - a subtypes isolated from outbreaks in iran and other countries were used to construct a fmdv - a vp1 and 3d - based sequence similarity tree .
2 shows a phylogenetic tree that was constructed based on the sequence alignment of 35 genomes of the vp1 region and 42 genomes in the 3d region which are distinctly divided into different lineages . as depicted in fig .
2a , the a / iran87 isolate clustered with a saudi arabian and five iranian isolates into a branch separate from other type a isolates .
2b demonstrates that a / iran87 and all isolates examined in the present study originated from different geographical areas and did not cluster in relatively similar lineages based on 3d sequences .
all the field isolates shared comparatively lower homology ( 90~93% ) with the a / iran87 isolate in the 3d coding region ( data not shown ) .
the topology of the phylogenetic tree indicated that there was no remarkable similarity between a / iran87 and all other isolates in the 3d coding region .
space - filling and worm styles of the three - dimensional structures of the 3d and vp1 proteins from the a / iran87 isolate were observed .
analysis of the 3d protein structure showed that the aspglu substitution occurred in a beta sheet located within a small groove of the protein . in order to examine the long amino acid deletion in the three - dimensional structure of the vp1 protein , homology - based modeling of the protein domains
it indicated that , 13 amino acid deletions in the vp1 protein of a / iran87 isolate caused a change in three - dimensional structure of protein in the g - h loop region .
detailed knowledge of the molecular characteristics of the major fmdv immunogenic components would be useful for monitoring various processes like evolution , genetic diversity , and virus origin , and for developing protective vaccines . among the fmdv serotypes ,
fmdv variants also appear during continuous infection of animals on farms or proliferation in cell cultures .
furthermore , genetic diversification among fmdv subtypes can be caused by events occurring during recombination . along with
above mechanisms , there is continuous co - circulation of multiple genotypes of fmdv type a which may lead to recombination ; this progresses gradually between closely related isolates when multiple viral genotypes
exchange of different genomic regions between fmdvs by recombination is directly linked to fmdv diversification .
recombination , particularly in the structural protein - coding region , may offer selective advantages to the virus
. this can be great concern to areas where multiple fmdv genotypes co - circulate .
the presence of extensive antigenic variation among different fmdv isolates can hinder vaccination against fmdv .
genetic polymorphisms amongst different type a subtypes demonstrated that four nucleotides in the 3d coding region of isolate a / iran87 ( 49 , 78 , 148 , and 170 ) were changed from g to t , c to a , g to a , and c to g , respectively , resulting in amino acid substitutions ala17ser , asp26glu , glu50asn , and ala57gly , respectively .
in contrast to the a / iran87 isolate , the remaining 41 isolates showed no corresponding changes at these positions .
analysis of the sequence data revealed that the a / iran87 isolate also contained 30 nucleotide deletions in region 501~540 of the vp1 nucleotide sequence , which led to long amino acid deletions ( rgdlgslaarvaa ) in the g - h loop .
it should be noted that the rgd sequence was deleted from the g - h loop .
the three - dimensional and antigenic structures of many different serotypes of the fmdv have been examined . in spite of losing the rgd sequence and acquiring
a asp26glu substitution in a beta sheet located in a small groove of the 3d protein , the virus grew in a bhk-21 suspension cell culture and exerted cytopathic effects after 16~18 h. since this strain is used as a vaccine strain , it may be inferred that the rgd does not have a key role in the virus binding to cells during infection .
was found to be deleted , such a natural deletion in vp1 gene of fmdv is a novel phenomenon and has not been previously reported in compared isolates .
the vp1 capsid protein contains a mobile loop between the g and h strands on the virus surface ; this protein not only contains the major immunodominant epitopes of the virion but also a highly conserved rgd amino acid sequence motif .
several studies using different approaches have indicated that naturally developed fmdv isolates attach to cells via the highly conserved rgd motif .
fmdvs have been shown to use multiple rgd - dependent integrins of the v subgroup to initiate infection ; these include v3 , v6 , v1 , and v8 .
on the other hand , these viruses are capable of entering cells via non - integrin pathways .
tissue culture - adapted fmdvs are able to utilize heparan sulfate ( a glycosaminoglycan found on the cell surface ) as a receptor for entering the cells .
it has recently been found that fmdv can utilize unknown receptors except for integrin and heparan sulfate .
this indicates the presence of other receptors and possible alternative mechanisms for viral particle entrance into the cells .
this conclusion is also supported by the finding that destruction of the g-h loop rgd motif by site - directed mutagenesis created viruses which were able to replicate and assemble , but were not capable of initiating further infection cycles because of their inability to interact with cell receptors .
previous findings from studies of synthetic peptides suggested that the carboxy - terminus of the vp1 protein situated in the vicinity of the rgd motif is required for rgd - mediated cell binding .
the g - h loop of the vp1 protein on the viral particle surface is considered to be the dominant epitope for several fmdv types and the most variable part of the particle .
antibody - based responses against the loop have the ability to neutralize fmdvs in vitro and are protective in vivo .
for example , replacing the whole or partial g - h loop leads to the doubling of antibody titers against type a virus .
new findings in this field are in complete agreement with a previous report documenting that the only mutation which did not increase heterologous responses was deletion of the rgd - motif ( 3a ) . in another study , the rgd sequence was deleted from a genetically engineered fmdv ( type a12 ) , thereby rendering the virus particles incapable of binding to cell surfaces .
protective immunity in cattle has been obtained by utilizing particles from which the receptor binding site was deleted as a vaccine , revealing that this deletion has no significant effect on the antigenic properties of viruses with the rgd sequence .
many reports   have suggested that fmdv subtypes can be obtained in the field by immune selection during outbreaks of fmdv which have not been effectively neutralized in animal populations with incomplete immunization .
thus , key changes in protein sequences that give rise to differences in immunogenicity and virulence can be seen among viruses isolated from partially immunized animals .
in addition , it is possible that other significant changes take place when the virus moves from one host cell type to another .
for instance , this may occur when fmdv spreads from one species to another ( from cattle to swine and then back to cattle ) or when fmdv cultured in bhk cells infects animal herds . for a number of years
, it has been observed that fmdv easily adapts to many different tissue culture cell types .
a recent report on a single passage of subtype a12 ( avr1 ) in bhk cells showed that some viral variants selected from this passage contained amino acid changes at residues 148 and 153 of the vp1 , including ones that have a limited capability for eliciting cross - neutralizing antibodies .
the rapid generation of immunogenically distinguishable variants among fmdv subtype a viruses may affect the large - scale production of vaccines in bhk cells .
adaptation for growth in tissue culture can be used to select viruses that have antigenic and cell attachment site(s ) different from those of the parental virus .
amino acid residues 134~158 form a g-h loop structure which contains the main fmdv immunogenic epitope . in fmdv
type c , mutation of residues 138~140 and 148~150 in the vp1 region has been shown to affect fmdv antigenicity .
mutation frequency in the vp1 gene sequence of the fmdv has been estimated to be 1.6  10~6.4  10 substitutions for each nucleotide per year .
g - h loop - specific antibody responses are known to play a major role in immunity induced by the current fmdv vaccines . considering the high sequence variability of this loop within and among the serotypes , it will be valuable to direct immune responses against other protective epitopes that are not as variable .
residues 140~160 of vp1 have been previously shown to induce neutralizing antibodies against fmdv types o and a . due to a novel deletion of the 13 amino acid region including the rgd motif in the vpi capsid protein of a / iran87 , the effect of this motif in receptor binding and cell infection must be determined .
this protein was initially called the fmdv infection - associated antigen since antibodies to this antigen are be detected in serum from recovery animals .
the 3d position is located at the c - terminus of the polypeptide next to the 3c gene product .
the 3d coding region produces a 470-amino acid protein with a molecular weight of ~55 kda .
the pivotal role of 3d in the viral replication process is attributed to the fact that the 3d gene is highly conserved , particularly in the functional motifs .
since rna polymerase is necessary for viral replication , the degree of sequence identity among the 3d genes is much greater than that of vp1 gene sequences .
a wealth of three - dimensional structural information is currently available for a large number of rna - dependent rna polymerases from different families of positive- and double - stranded rna viruses , including three different members of the picornaviridae family .
considering the fact that this strain can infect bhk-21 cells and has been used as a vaccine strain for nearly two decades , it may be concluded that deletion of rgd has no critical role in virus binding to cells during the initiation of infection .
the findings of the current study demonstrated that the cell - adapted vaccine strain a / iran87 containing 13-amino acid deletions in the vp1 protein can be considered as a novel variant of fmdv type - a .
 the nucleotide sequence of the vp1 ( 1d ) and partial 3d polymerase ( 3dpol ) coding regions of the foot and mouth disease virus ( fmdv ) vaccine strain a / iran87 , a highly passaged isolate ( ~150 passages ) , was determined and aligned with previously published fmdv serotype a sequences . 
 overall analysis of the amino acid substitutions revealed that the partial 3dpol coding region contained four amino acid alterations . 
 amino acid sequence comparison of the vp1 coding region of the field isolates revealed deletions in the highly passaged iranian isolate ( a / iran87 ) . 
 the prominent g - h loop of the fmdv vp1 protein contains the conserved arginine - glycine - aspartic acid ( rgd ) tripeptide , which is a well - known ligand for a specific cell surface integrin . despite losing the rgd sequence of the vp1 protein and an asp26glu substitution in a beta sheet located within a small groove of the 3dpol protein , the virus grew in bhk 21 suspension cell cultures . 
 since this strain has been used as a vaccine strain , it may be inferred that the rgd deletion has no critical role in virus attachment to the cell during the initiation of infection . 
 it is probable that this fmdv subtype can utilize other pathways for cell attachment . 