Isolation, identification and complete genome sequence analysis of a strain of foot-and-mouth disease virus serotype Asia1 from pigs in southwest of China
© Yang et al; licensee BioMed Central Ltd. 2011
Received: 18 February 2011
Accepted: 16 April 2011
Published: 16 April 2011
Foot-and-mouth disease virus (FMDV) serotype Asia1 generally infects cattle and sheep, while its infection of pigs is rarely reported. In 2005-2007, FMD outbreaks caused by Asia1 type occurred in many regions of China, as well as some parts of East Asia countries. During the outbreaks, there was not any report that pigs were found to be clinically infected.
In this study, a strain of FMDV that isolated from pigs was identified as serotype Asia1, and designated as "Asia1/WHN/CHA/06". To investigate the genomic feature of the strain, complete genome of Asia1/WHN/CHA/06 was sequenced and compared with sequences of other FMDVs by phylogenetic and recombination analysis. The complete genome of Asia1/WHN/CHA/06 was 8161 nucleotides (nt) in length, and was closer to JS/CHA/05 than to all other strains. Potential recombination events associated with Asia1/WHN/CHA/06 were found between JS/CHA/05 and HNK/CHA/05 strains with partial 3B and 3C fragments.
This is the first report of the isolation and identification of a strain of FMDV type Asia1 from naturally infected pigs. The Asia1/WHN/CHA/06 strain may evolve from the recombination of JS/CHA/05 and HNK/CHA/05 strains.
Foot-and-mouth Disease is a highly contagious and economically important disease of cloven-hoofed animals, predominantly for cattle, sheep, and pigs. The aetiological agent, foot-and-mouth disease virus, is classified as small icosahedral virus of the Aphthovirus group within the Picornaviridae family. There are seven immunologically distinct serotypes of the virus, namely types O, A, C, SAT1, SAT2, SAT3 and Asia1, and subtypes have also been found within some serotypes .
FMD serotype Asia1 has been epidemic in China for more than 50 years. This serotype usually infects cattle and sheep, and its infection of pigs is seldom reported . In 2005-2007, FMD outbreaks caused by Asia1 type occurred in many regions of China, as well as some parts of East Asia countries. During the outbreaks, there was not any report that pigs were found to be clinically infected .
FMDV has a single-stranded positive sense RNA genome of approximately 8.5 kb in length, including the 5' untranslated region (5'UTR), a large singe open reading frame (ORF), and the 3' untranslated region (3'UTR) [3, 4]. The 5' UTR consists of a short (S) fragment, a poly (C) tract, and a long fragment (5'LF-UTR), which contains three or four tandemly repeated pseudoknots (PKs) and an internal ribosome entry site (IRES) . The ORF encodes a polyprotein that can be cleaved to form four structural proteins (VP4, VP2, VP3 and VP1) and 8 non-structural proteins (L, 2A, 2B, 2C, 3A, 3B, 3C and 3D) . The VP1 protein plays an important role in virus attachment, protective immunity and serotype specificity, and nucleotide sequencing of this region has been extensively used for molecular epidemiology studies on FMD [6, 7]. The G-H loop of the VP1 protein of FMDV spanning residues 134-158 contains conserved Arg-Gly-Asp (RGD) tripeptide, which is considered to be a ligand for cell-surface attachment . In addition, FMDV 3A region has been implicated in virus virulence and host range, similar to the 3A proteins of other picornaviruses . The 3'UTR of about 90 residues with a poly (A) tail (35-100nt) at 3'-end is likely to be a site of interaction with viral and host proteins for RNA replication [10, 11].
This study, for the first time to our knowledge, described the isolation and identification of a strain of FMDV type Asia1 from pigs in China. To investigate the genomic feature of the strain and further understand its role in epidemiology of FMDV, the complete genome of the strain was sequenced and compared with sequences of other strains of FMDV Asia1 by phylogenetic and recombination analysis.
Materials and methods
Sample collection and clinical samples treatment
Three samples of ruptured vesicular fluids were collected from FMD-suspected pigs in a pig farm in southwest of China in 2006. The samples were transported from the collection site to diagnostic laboratory in 0.04 M phosphate buffer (pH 7.2) with 50% glycerol at 4°C and stored at -20°C until tested.
Virus isolation and identification
Established cell layer of BHK-21 cells were inoculated with 0.2 ml the three samples of vesicular fluids, respectively. The cell cultures would be examined for Cytopathic effects (CPE) for 48 hours. If CPE was formed, the cells would be collected for subsequent experiments. If no CPE was detected, the cells would be frozen and thawed, and used to inoculate fresh cultures of 0.2 ml and examined for CPE for another 48 hours. The infected BHK-21 monolayer cells were subjected to three freeze-thaw cycles to release the viral particles. The viral suspension was clarified from the cell debris by centrifugation at 800 × g for 10 min and stored at -70°C for the following experiments.
The FMDV O, Asia1 and A type positive serums were chosen as antiserum in complement fixation test (CFT). The CFT was performed by adding 0.2 ml virus sample, each of the 3 type of antiserum and complement in tubes, respectively. After incubating the mixtures at 37°C for 1.5 h, sensitized sheep erythrocytes was added to each tube, and the mixtures were incubated for 1 h at 37°C. Degrees of hemolysis were determined by visual reading. The tests were set up in triplicate. Also, virus sample, antisera, complement, and erythrocytes controls were set up with each test.
RNA extraction and RT-PCR
Synthetic oligonucleotides used for amplification of the FMDV genomea
Used for 5'-RACE
Used for 3'-RACE
The target fragment of PCR products were purified using Gel Extraction Kit (OMEGA bio-tek). The purified products were ligated to pMD18-T vector (TaKaRa, Japan) and transfected into E. coli JM109. The positive clones were selected to provide for sequencing by Shanghai Sangon Biological Engineering Technology and Services Co., Ltd. Each nucleotide was determined from three identical results. Sequences were assembled into complete genome sequence using SeqMan || program of DNAstar software package (DNAstar Inc., USA)
Sequence comparison and phylogenetic analysis
Reference FMDV strains selected from GenBank
Genbank Accession no.
Species of origin
Country and year of isolation
Genbank Accession no.
Species of origin
Country and year of isolation
Hong Kong, 2005
Hong Kong, 2005
Viet Nam, 2005
Viet Nam, 2005
Hong Kong, 2005
Viet Nam, 2005
Hong Kong, 2002
Putative recombination analysis
The complete genome sequences based on the multiple alignment result were used for recombination analysis with the SimPlot program version 3.5.1. Similarity plot were built using Asia1/WHN/CHA/06 as the query sequence against the complete genomic sequences of 8 strains of serotype Asia1. To further investigate the potential recombination event, the Recombination Detection Program (RDP) , GENECONV , BOOTSCAN , and MaxChi  methods were carried to employ Recombination Detection Program 3 (RDP3)  software. Statistical significance was set at the P = 0.05 level. Each analysis was conducted twice to ensure repeatability of results.
Virus isolation and identification
Virus from one sample was adapted to BHK-21 cells and CPE occurred when inoculated fresh cultures first time after 24 hours, but the other two didn't show any CPE by blind passage 5 times. Therefore, one strain of virus was isolated, and the first passage of the strain in BHK cells was used for following experiments. For complement fixation test, hemolysis occurred in the mixtures containing FMDV O and A type positive serums, whereas the mixtures containing FMDV Asia1 type positive serums didn't show any virtual change. Thus, the isolate was identified as FMDV of serotype Asia1, and was designated as Asia1/WHN/CHA/06.
RT-PCR and complete genome of Asia1/WHN/CHA/06
In order to determine the full-length genome sequence of Asia1/WHN/CHA/06 strain, the seven RT-PCR cDNA amplicons covering the entire RNA genome were cloned and each nucleotide was determined from three identical results. The genome sequence was submitted to GenBank with the accession number of FJ906802. The complete genome sequence of Asia1/WHN/CHA/06, excluding poly(C) and poly(A) tail, is 8161 nt in length, containing an ORF, a 5' UTR, and a 3' UTR. ORF is 6990 nt (encoding 2329 aa)in length, which consists of L (603 nt), P1 (2199 nt), P2 (1464 nt), and P3 (2724 nt) genes. P1 protein is predicted to be cleaved into four structural proteins, including VP4 (85 aa), VP2 (219 aa), VP3 (218 aa), and VP1 (211 aa). There are three non-structural proteins in P2, including 2A (16 aa), 2B (154 aa), and 2C (318 aa). P3 protein contains 3A (153 aa), 3B (71 aa), 3C (213 aa), and 3D (471 aa) four proteins. The 5' UTR consists of S, FUR, and IRES structures, which is 366, 258, and 454 nt in length, respectively. The 3' UTR of Asia1/WHN/CHA/06 was 93 nt in length, followed by a 66 bp poly (A) tail at least.
Similarity comparison of nucleotide and amino acid sequences of Asia1/WHN/CHA/06 to other FMDV Asia1 type isolates
Nucleotide (amino acid) sequence similarity (%) with Asia1/WHN/CHA/06
The ORF of Asia1/WHN/CHA/06 was 89.6 - 98.5% identical with the six reference sequences at the nucleotide level with no deletion or insertion. Compared with 40 sequences of VP1 genes, Asia1/WHN/CHA/06 exhibited identities ranging from 81.0 to 99.2% (data not shown). By comparing the 40 amino acid sequences of G-H loop region in VP1 protein, Asia1/WHN/CHA/06 got highly consistency with the stains isolated after the year of 2004, and no mutation was observed in RGD tripeptide. Interestingly, RGD mutants (RGD→HGD) were examined in Asia1/3kimron/61 and KRG/2/04. For bordering region of RGD of Asia1/WHN/CHA/06, there was no amino acids' mutation at RGD +1 and +4 sides, but a mutation at the RGD +9 (N→D). The 3A gene of Asia1/WHN/CHA/06 strain was 459 nt without insertions and deletions at this region. However, three amino acids' mutations occurred at the residues 86 (D→N), 93 (I→N), and 131 (H→R). Interestingly, Asia1/WHN/CHA/06 displayed the highest nucleotide similarity with JS/CHA/05 at all the regions except 3B gene (93.4%), whereas it showed the highest similarity (97.7%) with HNK/CHA/05 at 3B gene.
The 5' UTR of Asia1/WHN/CHA/06 shared 96.8% identity with the JS/CHA/05 strain, whereas it showed only 84.4-86.7% identity with other 5 reference strains. For the S fragment of 5'UTR, Asia1/WHN/CHA/06 displayed 3 nucleotide deletions located at nt 72, 74, and 184. The 712-nt of Asia1/WHN/CHA/06 of 5' LF-UTRs was observed, whereas within other reference strains, the length of 5' LF-UTRs ranged from 690 nts to 714 nts. Asia1/WHN/CHA/06 strain showed four PKs structure in 5' LF-UTR, which was the same as other reference isolates, and no nucleotide deletion was observed. For the IRES fragment of 5' UTR, Asia1/WHN/CHA/06 was more genetic conservative than that of S fragment, displaying 86.8-99.3% homology to other Asia1 type isolates, and only found to be 2 nucleotide deletions. The sequence comparison of the 3' UTR regions showed that Asia1/WHN/CHA/06 exhibited 74.2-87.1% homology to other six Asia1 type isolates, indicating a relatively higher level of genetic variation within the 3' UTR than the 5'UTR.
Asia1/WHN/CHA/06 was isolated from pigs in a pig farm in 2006. According to the investigation of the pig farm during the outbreak, about 20% of sows presented typical clinical symptoms of FMD, whereas none of the piglets and boars showed any symptoms and no pigs died. The FMD in the pig farm lasted for about one month. Moreover, there were not any reports about outbreak of FMD in other pig farms around the area at that time. Therefore, we speculated that Asia1/WHN/CHA/06 had low pathogenicity or infectivity to pigs, though high lethality of the strain to sulking mouth was observed in this study. But further study should be carried out to test the virulence of the strain to pigs to approve the inference.
The complete genome sequence of Asia1/WHN/CHA/06 was determined, except Poly(C) tract. Poly(C) tract of FMDV is usually ranging from 100 to 420 cytosine residues interspersed with an occasional uracil residue , so it is difficult to amplify the true sequence of the tract by RT-PCR from the viral RNA genome . In this study, we didn't manage to get the poly(C) tract sequence of Asia1/WHN/CHA/06 either. Pseudoknots (PKs) have been found to be 3-4 tandem repeats of RNA secondary structure in FMDV genome, which probably had some role during viral replication . Additionally, it is reported that partial deletions of FMDV strains in PKs' region correlated with the host range . However, no deletions occurred in PKs' region of Asia1/WHN/CHA/06. In sequence of VP1 protein, the RGD motif is the molecular basis of binding of FMDV to its cellular receptor . Although the RGD bordering region sequence of Asia1/WHN/CHA/06 is highly consistent with the reference strains isolated after the year of 2004, an amino acids' position at RGD +9 mutated (N→D). As several lines of evidence show that the vicinity of the RGD motif, including RGD +1 (leucine) and +4 (leucine) positions, are necessary for the FMDV RGD-mediated receptor recognition [24, 25], it is unclear whether the variation at RGD +9 site influence the receptor recognition of Asia1/WHN/CHA/06. Although non-structural proteins which play important roles for viral propagation are generally more conserved than the structural proteins, 3A is one of the most variable proteins of FMDV . In addition, deletions in 3A have been shown to associate with FMDV attenuation in cattle and high virulence for pigs . Meanwhile, a single amino acid mutation at position 44 in the 3A protein is present in guinea pig-adapted isolates examined . In our study, we found no deletion but 3 amino acid substitutions in the 3A region of Asia1/WHN/CHA/06. Considering that the isolate isolated from a pig, while other FMDV Asia1 type isolates came from cattle, it would be interesting to examine by reverse genetics if the 3 amino acid substitutions may contribute to adaptation of host range.
From a practical point of view, it is important to analyze the nucleotide sequence of the virus from new hosts for subsequent comparison with other isolates to trace the origin and route of the viral spread. To date, there are over 40 complete genome sequences for FMDV Asia1 isolated from cattle available in Genbank, but no data are available on pigs. The present study reported that the complete nucleotide sequence of Asia1/WHN/CHA/06 was isolated from pigs in southwest of China in 2006. In general, the complete genome sequence of Asia1/WHN/CHA/06 exhibited a high level of similarity with JS/CHA/05 strain, which indicated that Asia1/WHN/CHA/06 had a close genetic relationship with JS/CHA/05 strain. Phylogenetic analysis of the virus VP1 region of FMDV has been used extensively to investigate the molecular epidemiology of the disease worldwide . In this study, phylogenetic tree based on VP1 region of Asia1 isolates identified the existence of three groups. Group1a mainly consisted of 11 isolates, which collected in East Asia (China, East of Russia, East of Mongolia, and North Korea) during 2005-2007. However, India/80 and India/81, which were both collected from India in 1980s, also clustered with Group1a. These data suggest that the East Asia strains, including Asia1/WHN/CHA/06, have a close genetic relationship with the two India strains emerged two decades ago.
Recombination plays an important role in FMDV evolution [28, 29]. In this study, supported by Simplot and RDP program, Asia1/WHN/CHA/06 was found to be a mosaic between JS/CHA/05 and HNK/CHA/05, whose breakpoint encompassed partial 3B and 3C genes. Protein 3B, which is present in three similar but nonidentical copies (3B1, 3B2 and 3B3) in FMDV , functions in priming picornavirus RNA synthesis , and the 3C protease is responsible for most of the cleavages during the FMDV polyprotein processing . Recombination in FMDV occurs mainly in non-structural genes , but the recombination over 3B/3C gene boundaries appears to be less frequent . Thus, the role of the mosaicism observed in 3B/3C of FMDV should be targeted for further study. Moreover, JS/CHA/05 and HNK/CHA/05 strains were both isolated from the east of China in 2005, so there was a geographical opportunity for the recombination of the two strains. However, as Asia1/WHN/CHA/06 was isolated from west of China, it was possible that Asia1/WHN/CHA/06 was transmitted from east of China by introduction of pigs or other transmission routes.
Here we reported, for the first time to our knowledge, the isolation and identification of a strain of FMDV type Asia1 from naturally infected pigs, and described the features of the genome of Asia1/WHN/CHA/06. Sequence analysis showed that it was belonged to East Asian strains, and had a close genetic relationship to JS/CHA/05. Besides, potential recombination events associated with Asia1/WHN/CHA/06 were found between JS/CHA/05 and HNK/CHA/05 strains with partial 3B and 3C fragments, which may trace the origin and evolution of Asia1/WHN/CHA/06 strain.
This research was supported by China National '11th Five-Year Plan' Science and Technology Supporting Program: Research and Demonstration for Healthy Breeding of new model in intensive pig farms (2006BAD14B05).
- Domingo E, Escarmis C, Baranowski E, Ruiz-Jarabo CM, Carrillo E, Nunez JI, Sobrino F: Evolution of foot-and-mouth disease virus. Virus Res 2003, 91:47–63.PubMedView Article
- Zhang Q, Li D, Liu X, Liu Z, Cai X, Wu G, Qi S, Yang S, Yan X, Shang Y, He J, Ma J, Li J, Ma W, Han R, Liu X, Zhang J, Xie Q, Zhang Z: Experimental studies with foot-and-mouth disease virus type Asia-1, responsible for the 2005 epidemic in China. Res Vet Sci 2008, 85:368–371.PubMedView Article
- Lau LT, Reid SM, King DP, Lau AMF, Shaw AE, Ferris NP, Yu ACH: Detection of foot-and-mouth disease virus by nucleic acid sequence-based amplification. Vet Microb 2008, 126:101–110.View Article
- Grubman MJ, Baxt B: Foot-and-mouth disease. Clin Microb 2004, 17:465–493.View Article
- Belsham GJ: Translation and replication of FMDV RNA. Curr T Micr 2005, 288:43–70.View Article
- Du J, Chang H, Cong G, Shao J, Lin T, Shang Y, Liu Z, Liu X, Cai X, Xie Q: Complete nucleotide sequence of a Chinese serotype Asia1 vaccine strain of foot-and-mouth disease virus. Virus Genes 2007, 35:635–642.PubMedView Article
- Malirat V, Barros JJF de, Bergmann IE, Campos RM, Neitzert E, Eliane VC, Da Silva EE, Falczuk AJ, Pinheiro DSB, Natalia V, Jose LQC, Eduardo M, Rosa DL: Phylogenetic analysis of foot-and-mouth disease virus type O re-emerging in free areas of South America. Virus Res 2007, 124:22–28.PubMedView Article
- Brown F, Benkirane N, Limal D, Halimi H, Newman JFE, Regenmortel MHV, Briand JP, Muller S: Delineation of a neutralizing subregion within the immunodominant epitope (GH loop) of foot-and-mouth disease virus VP1 which does not contain the RGD motif. Vaccine 2008, 18:50–56.View Article
- Lama J, Sanz MA, Carrasco L: Genetic analysis of poliovirus protein 3A: characterization of a non-cytopathic mutant virus defective in killing Vero cells. J Gen Virol 1998, 79:1911–1921.PubMed
- Domingo E, Baranowski E, Escarmis C, Sobrino F: Foot-and-mouth disease virus. Comp Immu M I 2002, 25:297–308.
- Carrillo C, Tulman ER, Delhon G, Lu Z, Carreno A, Vagnozzi A, Kutish GF, Rock DL: Comparative genomics of foot-and-mouth disease virus. J Virol 2005, 79:6487–6504.PubMedView Article
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Ev 2007, 24:1596–1599.View Article
- Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol 1999,73(1):152–160.PubMed
- Martin DP, Rybicki E: RDP: detection of recombination amongst aligned sequences. Bioinformatics 2000, 16:562–563.PubMedView Article
- Padidam M, Sawyer S, Fauquet CM: Possible emergence of new geminiviruses by frequent recombination. Virology 1999, 265:218–225.PubMedView Article
- Martin DP, Posada D, Crandall KA, Williamson C: A modified BOOTSCAN algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res Hum Retroviruses 2005, 21:98–102.PubMedView Article
- Smith JM: Analyzing the mosaic structure of genes. J Mol Evol 1992, 34:126–129.PubMed
- Heath L, Walt E, Varsani A, Martin DP: Recombination patterns in aphthoviruses mirror those found in other picornaviruses. J Virol 2006, 80:11827–11832.PubMedView Article
- Escarmis C, Toja M, Medina M, Domingo E: Modification of the 5' untranslated region of foot-and-mouth disease virus after prolonged persistence in cell culture. Virus Res 1992, 26:113–125.PubMedView Article
- Rieder E, Bunch T, Brown F, Mason PW: Genetically engineered foot-and-mouth disease viruses with poly(C) tracts of two nucleotides are virulent in mice. J Virol 1993, 67:5139–5145.PubMed
- Biswas S, Sanyal A, Hemadri D, Tosh C, Mohapatra JK, Manoj KR, Bandyopadhyay SK: Sequence analysis of the non-structural 3A and 3C protein-coding regions of foot-and-mouth disease virus serotype Asia1 field isolates from an endemic country. Vet Microb 2006, 116:187–193.View Article
- Dunn CS, Donaldson AI: Natural adaption to pigs of a Taiwanese isolate of foot-and-mouth disease virus. Vet Rec 1997, 141:174–175.PubMedView Article
- Storey P, Theron J, Maree FF, O'Neill HG: A second RGD motif in the 1D capsid protein of a SAT1 type foot-and-mouth disease virus field isolate is not essential for attachment to target cells. Virus Res 2007, 124:184–192.PubMedView Article
- Jackson T, Blakemore W, Newman JWI, Knowles NJ, Mould AP, Humphries MJ, King AMQ: Foot-and-mouth disease virus is a ligand for the high-affinity binding conformation of integrinα5β1: influence of the leucine residue within the RGDL motif on selectivity of integrin binding. J Gen Virol 2000, 81:1383–1391.PubMed
- Jackson T, Mould AP, Sheppard D, King AMQ: Integrin αvβ1 is a receptor for foot-and-mouth disease virus. J Virol 2002, 76:935–941.PubMedView Article
- Nunez JI, Baranowski E, Molina N, Ruiz-Jarabo CM, Sanchez C, Domingo E, Sobrino F: A single amino acid substitution in non-structural protein 3A can mediate adaptation of foot-and-mouth disease virus to the guinea pig. J Virol 2001, 75:3977–3983.PubMedView Article
- Ayelet G, Mahapatra M, Gelaye E, Egziabher BG, Rufeal T, Sahle M, Ferris NP, Wadsworth J, Hutchings GH, Knowles NJ: Genetic Characterization of Foot-and-Mouth Disease Viruses, Ethiopia, 1981–2007. Em Infect D 2009, 15:1409–1417.View Article
- Haydon DT, Bastos AD, Awadalla P: Low linkage disequilibrium indicative of recombination in foot-and-mouth disease virus gene sequence alignments. J Gen Virol 2004, 85:1095–1100.PubMedView Article
- Rensburg HG, Haydon D, Joubert F, Bastos AD, Heath L, Nel L: Genetic heterogeneity in the foot-and-mouth disease virus Leader and 3C proteinases. Gene 2002, 289:19–29.PubMedView Article
- Forss S, Schaller H: A tandem repeat gene in a picornavirus. Nucleic Acids Res 1982, 10:6441–6450.PubMedView Article
- Wimmer E: Genome-linked proteins of viruses. Cell 1982, 28:199–201.PubMedView Article
- Vakharia VN, Devaney MA, Moore DM, Dunn JJ, Grubman MJ: Proteolytic processing of foot-and-mouth disease virus polyproteins expressed in a cell-free system from clone-derived transcripts. J Virol 1987, 61:3199–3207.PubMed
- Klein J, Hussain M, Ahmad M, Normann P, Afzal M, Alexandersen S: Genetic characterisation of the recent foot-and-mouth disease virus subtype A/IRN/2005. Virol J 2007, 4:122.PubMedView Article
- Jackson AL, O'Neill H, Maree F, Blignaut B, Carrillo C, Rodriguez L, Haydon DT: Mosaic structure of foot-and-mouth disease virus genomes. J Gen Virol 2007, 88:487–492.PubMedView Article
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.