The molecular epidemiology of Foot-and-Mouth Disease virus serotypes A and O from 1998 to 2004 in Turkey
© Klein et al; licensee BioMed Central Ltd. 2006
Received: 07 July 2006
Accepted: 04 December 2006
Published: 04 December 2006
Foot-and-Mouth Disease (FMD) causes significant economic losses in Turkish livestock.
We have analysed the genetic diversity of the 1D sequences, encoding the hypervariable surface protein VP1, of Turkish isolates of serotype A and O collected from 1998 to 2004 in order to obtain epidemiological and immunological information.
The 1D coding region of 33 serotype O and 20 serotype A isolates, obtained from outbreaks of FMD between 1998 and 2004, was sequenced.
For serotype A, we confirmed the occurrence of the two subtypes IRN99 and IRN96. These subtypes are most divergent within the region encoding the immuno-dominant GH-loop. Also a close relationship to Foot-and-Mouth Disease virus (FMDV) serotype A isolates obtained from outbreaks in Iraq and Iran were detected and a clustering of isolates collected during the same period of time were found.
The analysis of the deduced amino-acid sequences of these subtypes revealed evidence of positive selection in one site and one deletion, both within the GH-loop region. By inferring the ancestral history of the positively selected codon, two potential precursors were found. Furthermore, the structural alignment of IRN99 and IRN96 revealed differences between the tertiary structures of these subtypes.
The similarity plot of the serotype O isolates suggested a more homogeneous group than the serotype A isolates. However, phylogenetic analysis revealed two major groups, each further divided in subgroups, of which some only consisted of Turkish isolates.
Positively selected sites and structural differences of the Turkish isolates analysed, were not found.
The sequence and structural analysis of the IRN99 strains is indicative of positive selection suggesting an immunological advantage compared to IRN96. However, results of antigenic comparison reported elsewhere do not substantiate such a conclusion. There is evidence that IRN99 was introduced to Turkey, in all probability from Iran.
Since, a member of the IRN96 lineage was included as a component of the FMDV vaccine produced since 2000, the outbreaks caused by IRN96 strains in 2004 could be due to incomplete vaccine coverage.
The Turkish type O strains, all with a VP1 structure similar to the O1/Manisa/69 vaccine, appear in several sublineages. Whether these sublineages reflect multiple samplings from a limited number of outbreaks, or if they reflect cross-boundary introductions is not clear.
Foot-and-Mouth Disease (FMD) is an acute, highly communicable and economically important disease caused by Foot-and-Mouth Disease virus (FMDV). Animals that can be affected include cattle, swine, sheep, goats, wild pigs, wild ruminants and buffaloes . FMDV is a positive sense single-stranded RNA virus (genus Aphthovirus, family Picornaviridae) occurring in seven serotypes, O, A, Asia 1, C, SAT 1, SAT 2 and SAT 3, each with a wide spectrum of antigenic and epidemiological different subtypes. The wide diversity is considered a consequence of the high mutation rate and quasi-species dynamic .
FMD has been one of the most important diseases causing significant economic losses in the Turkish livestock sector. Together with production losses, export restrictions on several agricultural products cause additional losses to the Turkish economy .
Turkey is of special interest for the European countries, because it is the natural connection to Asia, where FMD is enzootic and thus serves as buffer zone against FMD .
The two prevalent serotypes in Turkey are A and O , and especially the pandemic FMDV serotype O is a major threat to Europe, because it is the most aggressive serotype . However, this type O is still covered by the O1/Manisa/69 vaccine .
Serotype A displays a great antigenic diversity of subtypes and there is also often no cross-protection between them . Since 1999 two subtypes of serotype A have been circulating in Turkey, designated by the World reference Laboratory IAH Pirbright as IRN96 and IRN99 .
Turkey practices vaccination against FMD, using a self produced trivalent vaccine consisting of O1/Manisa/69, Asia1 and since September 2000 a member of the IRN96 lineage as a replacement for A22, previously included [personal communication].
Due to the introduction of a twice-yearly national vaccination campaign in 1996, the number of FMD outbreaks has decreased significantly, but not homogeneously over the whole country, which may be explained by vaccine failure or incomplete vaccine coverage .
A common and major epitope of FMDV is located within the surface protein VP1, containing the immuno-dominant GH-loop and the RGD-integrin binding motif, essential for cell attachment . Changes in this protein may cause vaccine failure and changes in host specificity .
We describe the genetic diversity of the 1D nucleotide sequence which encodes for the VP1 protein, in isolates obtained from outbreaks in Turkish cattle herds between 1998 and 2004 and analyse the data with regard to potential epidemiological information and to establish possibly whether FMD outbreaks are caused by viruses persistently circulating and evolving in Turkey or by strains introduced to Turkey.
Serotype O, with O1/Manisa/69 as reference sequence, displays one homogeneous group, with a diversity in the range 9%.
The IRN96 cluster is further subdivided into two temporal divergent groups, I and II. The isolates of group I all originate from outbreaks in 1998, it was two years before this lineage was included in the Turkish FMD vaccine.
The IRN99 related viruses display a single, not further subdivided, phylogenetic structure; with the oldest isolate DQ296550 Igdir/06.99, obtained from an outbreak in June 1999 in the province Igdir, close to the Iranian border. Another member of IRN99 lineage is the isolate A/IRN/22/99# (sequence information kindly provided from IAH Pirbright) obtained from an outbreak in Iran in 1999.
Cluster A is split into the two groups I and II. Group I consists of sequences derived from FMD outbreaks in Pakistan and Iran in addition to a single outbreak in January/February 2004 in the Mersin region of Turkey. The differential success, expressed as correlation to the number of substitutions, with regard to the vaccine lineage, is lowest in the group members obtained from outbreaks in Pakistan in 2002 and 2003, intermediate in the Iranian isolates from 2003 and 2004 and most advanced in the Turkish isolates from 2004.
Group II of cluster A consists of Turkish isolates obtained from outbreaks in the provinces Sivas in 2001 and Igdir in 1998. The latter province is close to the Iranian border.
The majority of Turkish isolates belong to cluster B, which can also be subdivided into two groups I and II.
Cluster I of group B was further subdivided into two lineages. One group consists of isolates, obtained only from outbreaks in Turkey and the other of isolates known as the Panasia lineage. This partition is supported by a clade credibility value of 0.45. The Turkish lineage displays an increasing differential success rate, with regard to the vaccine strain O1/Manisa/69, starting from the DQ296516 Ankara/12.99 isolate, obtained in the Ankara province in 1999, towards the DQ296528 Kutahya/05.04 isolate, obtained in the Kutahya province in 2004.
Cluster II of group B consists of viruses obtained from outbreaks in Turkey, Israel and Lebanon, with the earliest obtained from an outbreak in the Turkish Ankara region in 2001.
Positive selection analysis
By investigating the inferred ancestor history of site 151 (additional file 2 and 3), it can be seen that all but two of the IRN96 isolates, DQ296532 and DQ296533, each collected from the same outbreak, have the common ancestor codon 'AGC' (Ser), whereas IRN99 has most likely the 'GGG' (Gly) codon as precursor.
From the hydrophobicity plot, DQ296505 Sivas/02.01 seems to represent the most distinct lineage of the Turkish serotype O isolates (Figure 7a). However, the structural alignment between DQ296505 Sivas/02.01 and the vaccine strain O1/Manisa/69, displays a high similarity of 99% (Figure 7b).
The similarity plot (Figure 1), as well as the phylogenetic analysis of serotype A (Figures 2), confirms the presence of two independent lineages within our collection of isolates. A member of the IRN96 cluster was not included in vaccines before September 2000, which means that the vaccines before that did not provide sufficient protection against IRN96 strains. However, the IRN96 isolates (Figure 2, IRN96 II) collected from outbreaks in 2004, and the clustering of these isolates with the IRN96 isolates from Iran and Iraq alternatively suggests multiple introductions from one or both of these countries.
The IRN99 lineage is suggested to be a newly introduced lineage rather than an escape mutant lineage evolved from IRN96. This is supported by the lack of evolutionary intermediates between the two clusters, the consistent differences between isolates of the two clusters and the low probability of a common ancestor of residue 151 potentially subject to positive selection. In addition that the first IRN99 isolate reported in Turkey originate from an outbreak in the province of Igdir, close to the Iranian border. The structural alignment of IRN99 and IRN96 (Figure 6) revealed significant conformational differences within this important immuno-dominant protein suggesting an immunological advantage of IRN 99 strains given that a member of the IRN96 is used as a vaccine component. However, antigenic comparisons based on virus neutralisation analysis and epidemiological data previously reported  do not substantiate such a conclusion.
Recalling the fact that structures are more conserved in evolution than sequences, these observed differences support the idea of an independently evolved IRN99 lineage. Recombination events between the two serotype A types IRN99 and IRN96 within the 1D region were not detected by the similarity scan. This may indicate that both serotypes A types are in the phase of high fitness and have found their particular ecological niche.
In contrast to the serotype A strains, the serotype O isolates included in this study constitute a more homogeneous group, as shown by the similarity plot (Figure 1).
The phylogenetic analysis of serotype O isolates demonstrate the presence of two major subtypes, cluster A and B (Figure 3), each further subdivided in distinct lineages. Cluster A I contained isolates from one single outbreak in the Turkish Mersin province, as well as isolates obtained from earlier outbreaks in Pakistan and Iran. The temporal distribution and the differential success of these isolates, i.e. the branch lengths, suggest that the virus was introduced from Pakistan to Iran and then to Turkey. Cluster A II, consisting of the Turkish isolates DQ296505 Sivas/02.01 and DQ296504 Igdir/08.98, is well separated from A I and represent a sublineage only observed in Turkey. The same conclusion may be valid for the cluster B I, which also exclusively comprised Turkish isolates. However, the intermediate clade credibility value of 0.45 at the node between cluster B I and the Panasia lineage may indicate a closer relationship to the pandemic Panasia lineage than displayed in the phylogenetic tree. An explanation for these purely Turkish sublineages may be that the virus is circulating in Turkish cattle herds, either through subclinically infected cattle or by transmission to cattle from other reservoirs such as sheep and goats.
Cluster B II is clearly separated from the Panasia lineage, consisting of isolates from newer outbreaks in Turkey, Israel and Lebanon, indicating the spread of this sublineage.
Nevertheless, the structural comparison between the field isolates and the vaccine strain O1/Manisa/69 (Figure 8) displayed no distinct differences suggesting that the vaccine strain, O1/Manissa/69, in current use should provide full protection. This suggestion is concordant with other studies [11, 12].
In conclusion, this study confirmed the presence of two major serotype A lineages, namely IRN 99 and IRN96, in Turkey. There was no evidence that the A22 strain was present in Turkey during the examined period of time, however this may be related to the small sample size used in this study. The sequence and structural analysis of the IRN99 related viruses revealed potential conformational differences of the IRN99 VP1 possibly affecting the antigenic properties of particular regions. There is evidence that IRN99 is introduced from outside to Turkey, in all probability from Iran.
Since, the IRN96 lineage is a component of the FMDV vaccine used since 2000, the outbreaks in 2004 occurred most probably due to incomplete vaccine coverage.
The endemic, slowly evolving Turkish O lineages, with a predicted VP1 structure similar to the O1 Manisa vaccine, occurs in several sublineages in Turkey. Some of these sublineages are most probably introduced from neighbouring countries, whereas others are circulating within Turkish cattle or are frequently reintroduced to cattle from persistently infected small ruminants, i.e. sheep and goats, which share pastures and transport vehicles.
Thirty-three serotype O (Accession numbers: DQ296497 – DQ296531) and twenty serotype A (Accession numbers: DQ296532 – DQ296552) field isolates were collected from 1998 to 2004 from bovine herds widely distributed in Turkey. All isolates were epidermal tissue samples representing vesicular lesions and had been stored in TRIZOL©-reagent. Province and time (month/year) of collection are indicated with the accession numbers.
RNA extraction, reverse transcriptase – PCR and cycle sequencing
Tissue (100–150 mg) was homogenized in 1 ml RNApro™ Solution (Qbiogene, USA) in a Lysing Matrix D tube (Qbiogene, Inc., USA) using a FP 120 Fast Prep™ Cell Disruptor (Qbiogene, USA). Total RNA was extracted using RNeasy-Mini Kit™ (Qiagen, Germany) according to manufacturer's instructions. cDNA synthesis was done using
Ready-To-Go™ You-Prime First-Strand Beads (GE Healthcare Life Sciences, Sweden), employing the primers NV27T and random hexamers pdN6.
Primers used in this study
Sequence (5' to 3')
The resulting PCR product was examined by electrophoresis, using a 1,2% agarose gel, containing 0,005% ethidium bromide, with a separation time of 1,5 hours at 6,5 V/cm.
Amplicons were extracted and purified from agarosegel with QIAquick Gel Extraction kit (Qiagen, Germany) and cycle-sequencing was then performed by Agowa GmbH, Germany.
Sequence assembling was performed with ContigExpress (VectorNTI©-software) and multiple alignment was performed by log-expectation comparison, using the MUSCLE (v.3.6) software . For the serotype A alignment manual editing at codon position 153 was done.
Similarity plots were performed using Simplot software  with the following parameters: For serotype A a windows size of 100 bp (step: 20 bp) with gap-stripping and Jukes-cantor correction was chosen, whereas a window of 200 bp (step:20 bp) with gap-stripping and Kimura (2-parameter) correction was chosen for serotype O.
For serotype A the HKY+G was used and Bayesian analysis was performed using MrBayes (v3.2)  with the following settings. The maximum likelihood model employed 2 substitution types ("nst = 2"), with base frequencies set to fixed values ("statefreqpr = fixed"). Rate variation across sites was modelled using a gamma distribution (rates = "gamma"). The Markov chain Monte Carlo search was run with 4 chains for 500000 generations, with trees begin sampled every 100 generations (the first 1000 trees were discarded as "burnin").
For serotype O the GTR+G was used and Bayesian analysis was performed using MrBayes (v3.2)  with the following settings. The maximum likelihood model employed 6 substitution types ("nst = 6"), with base frequencies set to variable values ("statefreqpr = dirichlet(1,1,1,1)"). Rate variation across sites was modelled using a gamma distribution (rates = "invgamma"). The Markov chain Monte Carlo search was run with 4 chains for 500000 generations, with trees begin sampled every 100 generations (the first 1000 trees were discarded as "burnin").
Selection pressure analysis
Positively selected sites (Codon-specific analyses of dN/dS) was identified using the Single Likelihood Ancestor Counting (SLAC) analysis , a modification of the Suzuki-Gojobori method , available at the Datamonkey web side . For both serotypes, the analysis was done with the HKY85 substitution model on phylogenetic trees inferred using the Neighbour-Joining method with a cut-off p value of 0.5.
Three-dimensional structure of the VP1 protein was obtained using comparative protein modelling , employing the CPHmodels 2.0 Server  and pairwise structure comparison was done using the DaliLite workbench . Molecule Visualization was done using RASMOL 2.7 software .
The structural analysis of the A serotype is based on the comparison of the translated amino acid consensus sequences from the IRN96 and IRN99 lineage. Both, IRN99 consensus and IRN96 consensus sequences, used the structure of serotype A-1061 (1ZBA:1; Resolution: 2 Å)  as template.
For serotype O viruses the structural analysis was based on the translated amino acid sequences of DQ296505 Sivas/02.01 and the vaccine strain O1/Manisa/69. Both used the structure of strain O-1860 (1FOD:1; Resolution: 2,6 Å)  as template.
We thank Gert Olsen and Anita Bagge Hemmingsen for excellent technical assistance, Søren Alexandersen, Jens Nielsen and Thomas Bruun Rasmussen for useful discussions and Nick Knowles for supplying sequence information.
This work was funded by the FMD improcon project of the European Union Improvement of Foot and Mouth disease control by ethically acceptable methods based on scientifically validated assays and new knowledge on FMD vaccines, including the impact of vaccination (grant no. SSPE-CT-2003-503603).
- Rohrer H, Olechnowitz AF: Maul-und Klauenseuche. Jena, Germany: VEB Gustav Fischer Verlag: 1980.Google Scholar
- Domingo E, Pariente N, Airaksinen A, Gonzalez-Lopez C, Sierra S, Herrera M, Grande-Perez A, Lowenstein PR, Manrubia SC, Lazaro E, Escarmis C: Foot-and-mouth disease virus evolution: exploring pathways towards virus extinction. Curr Top Microbiol Immunol. 2005, 288: 149-173.PubMedGoogle Scholar
- Askaroglu H: Report on the FMD situation and control programme in Turkey. European commission for the control of Foot-and-Mouth Disease Report of the Sixty-eighth Session of the Executive Committee Vilnius, Lithuania 7 and 8 November 2002. 2002, FAO RomeGoogle Scholar
- EUFMD-REPORTS 63rd Session – Appendix 3. [http://www.fao.org/AG/AGAInfo/commissions/en/eufmd/app/63/app03.html]
- Kitching R: Global epidemiology and prospects for control of foot-and-mouth disease. Curr Top Microbiol Immunol. 2005, 288: 133-148.PubMedGoogle Scholar
- FMD Global Update 2003–2004. [http://www.iah.bbsrc.ac.uk/virus/Picornaviridae/Aphthovirus/WRLFMD/EUFMD-Crete/Valarcher.ppt]
- Gilbert M, Aktas S, Mohammed H, Roeder P, Sumption K, Tufan M, Slingenbergh J: Patterns of spread and persistence of foot-and-mouth disease types A, O and Asia-1 in Turkey: a meta-population approach. Epidemiol Infect. 2005, 133: 537-545. 10.1017/S0950268804003516.PubMed CentralView ArticlePubMedGoogle Scholar
- Fox G, Parry NR, Barnett PV, McGinn B, Rowlands DJ, Brown F: The cell attachment site on foot-and-mouth disease virus includes the amino acid sequence RGD (arginine-glycine-aspartic acid). J Gen Virol. 1989, 70: 625-637.View ArticlePubMedGoogle Scholar
- Hernandez J, Valero ML, Andreu D, Domingo E, Mateu MG: Antibody and host cell recognition of foot-and-mouth disease virus (serotype C) cleaved at the Arg-Gly-Asp (RGD) motif: a structural interpretation. J Gen Virol. 1996, 77: 257-264.View ArticlePubMedGoogle Scholar
- Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982, 157: 105-132. 10.1016/0022-2836(82)90515-0.View ArticlePubMedGoogle Scholar
- Aktas S: Molecular epidemiology of foot and mouth disease types O and A in Turkey. Reading University Library 1998.Google Scholar
- Parlak Ü, Özyörük F, Knowles N, Armstrong R, Aktas S, Feray A, Cokcaliskan C, Christensen L: Characterisation of foot-and-mouth disease virus strains circulating in Turkey during 1996–2004. Submitted to Archives of Virology.
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.PubMed CentralView ArticlePubMedGoogle Scholar
- 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: 152-160.PubMed CentralPubMedGoogle Scholar
- Swofford D: PAUP*: Phylogenetic Analysis Using Parsimony (*And Other Methods), version b4. Sinauer Associates, Sunderland, MA: 1998.Google Scholar
- Nylander J: MrModeltest v2. Program distributed by the author. Evolutionary Biology Centre, Uppsala University, Sweden: 2004.Google Scholar
- Huelsenbeck J, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Pond SL, Frost SD: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics. 2005, 21: 2531-2533. 10.1093/bioinformatics/bti320.View ArticlePubMedGoogle Scholar
- Suzuki Y, Gojobori T: A method for detecting positive selection at single amino acid sites. Mol Biol Evol. 1999, 16: 1315-1328.View ArticlePubMedGoogle Scholar
- Marti-Renom MA, Madhusudhan MS, Fiser A, Rost B, Sali A: Reliability of assessment of protein structure prediction methods. Structure. 2002, 10: 435-440. 10.1016/S0969-2126(02)00731-1.View ArticlePubMedGoogle Scholar
- Lund O, Nielsen M, Lundegaard C, Worning P: CPHmodels 2.0: X3M a Computer Program to Extract 3D Models. Abstract at the CASP5 conference A102. 2002Google Scholar
- Holm L, Park J: DaliLite workbench for protein structure comparison. Bioinformatics. 2000, 16: 566-567. 10.1093/bioinformatics/16.6.566.View ArticlePubMedGoogle Scholar
- Roger A, Sayle E, James Milner-White E: RASMOL: biomolecular graphics for all. Trends Biochem Sci. 1995, 20: 374-10.1016/S0968-0004(00)89080-5.View ArticleGoogle Scholar
- Fry E, Newman J, Curry S, Najjam S, Jackson T, Blakemore W, Lea S, Miller L, Burman A, King A, Stuart D: Structure of Foot-and-mouth disease virus serotype A1061 alone and complexed with oligosaccharide receptor: receptor conservation in the face of antigenic variation. J Gen Virol. 2005, 86: 1909-1920. 10.1099/vir.0.80730-0.View ArticlePubMedGoogle Scholar
- Logan D, Abu-Ghazaleh R, Blakemore W, Curry S, Jackson T, King A, Lea S, Lewis R, Newman J, Parry N, Rowlands D, Stuart D, Fry E: Structure of a major immunogenic site on foot-and-mouth disease virus. Nature. 1993, 362: 566-568. 10.1038/362566a0.View ArticlePubMedGoogle Scholar
- Knowles NJ, Davies PR, Midgley RJ, Valarcher JF: Identification of a Ninth Foot-and-Mouth Disease Virus Type O Topotype and Evidence for a Recombination Event in its Evolution. 1993, Available online (10.05.2006):, [http://www.fao.org/AG/AGAInfo/commissions/docs/greece04/App24.pdf]Google Scholar
- Locher F, Suryanarayana VV, Tratschin JD: Rapid detection and characterization of foot-and-mouth disease virus by restriction enzyme and nucleotide sequence analysis of PCR products. J Clin Microbiol. 1995, 33: 440-444.PubMed CentralPubMedGoogle Scholar
- Mason PW, Pacheco JM, Zhao QZ, Knowles NJ: Comparisons of the complete genomes of Asian, African and European isolates of a recent foot-and-mouth disease virus type O pandemic strain (PanAsia). J Gen Virol. 2003, 84: 1583-1593. 10.1099/vir.0.18669-0.View ArticlePubMedGoogle Scholar