- Open Access
Phylodynamic reconstruction of O CATHAY topotype foot-and-mouth disease virus epidemics in the Philippines
© Di Nardo et al.; licensee BioMed Central Ltd. 2014
Received: 21 March 2014
Accepted: 8 August 2014
Published: 24 August 2014
Reconstructing the evolutionary history, demographic signal and dispersal processes from viral genome sequences contributes to our understanding of the epidemiological dynamics underlying epizootic events. In this study, a Bayesian phylogenetic framework was used to explore the phylodynamics and spatio-temporal dispersion of the O CATHAY topotype of foot-and-mouth disease virus (FMDV) that caused epidemics in the Philippines between 1994 and 2005. Sequences of the FMDV genome encoding the VP1 showed that the O CATHAY FMD epizootic in the Philippines resulted from a single introduction and was characterised by three main transmission hubs in Rizal, Bulacan and Manila Provinces. From a wider regional perspective, phylogenetic reconstruction of all available O CATHAY VP1 nucleotide sequences identified three distinct sub-lineages associated with country-based clusters originating in Hong Kong Special Administrative Region (SAR), the Philippines and Taiwan. The root of this phylogenetic tree was located in Hong Kong SAR, representing the most likely source for the introduction of this lineage into the Philippines and Taiwan. The reconstructed O CATHAY phylodynamics revealed three chronologically distinct evolutionary phases, culminating in a reduction in viral diversity over the final 10 years. The analysis suggests that viruses from the O CATHAY topotype have been continually maintained within swine industries close to Hong Kong SAR, following the extinction of virus lineages from the Philippines and the reduced number of FMD cases in Taiwan.
Foot-and-mouth disease (FMD) is an economically devastating transboundary disease of cloven-hoofed domestic and wild ruminants, causing an acute and highly contagious vesicular disease which can develop into a persistent infection. The aetiological agent is FMD virus (FMDV), a single-stranded RNA virus belonging to the Aphthovirus genus, family Picornaviridae. FMDV is characterised by high genetic variability and exists as seven different serotypes named as O, A, C, Asia 1, Southern African Territories (SAT) 1, SAT 2, and SAT 3 . As a consequence of their high mutation rate, FMDV lineages quickly diverge as they replicate and spread into new areas. Therefore, transmission of the virus through space and time directly defines the evolutionary patterns observed between related FMDV strains . In addition to the accumulation of nucleotide substitutions through errors, large block of sequence changes can be mediated via recombination between different FMDV genomes, further expanding its evolutionary repertoire. In this context, FMDV populations often exhibit extensive genetic and antigenic heterogeneity at both the molecular and geographical level, driven by co-circulation of multiple lineages, heterogenic mixed host populations, extensive animal movements and trade patterns . FMDV serotypes have evolved independently in different geographical regions to give rise to distinct genetic lineages, designated topotypes. Eleven topotypes have been defined for serotype O, based on phylogenetic relationships between available sequence data and a value of ~15% of nucleotide (nt) sequence difference in the VP1 coding region ,.
The O CATHAY FMDV topotype
The first FMDV strain belonging to the O CATHAY topotype was isolated from Hong Kong SAR from pig samples collected during 1970 (HKN/21/70, GenBank accession no. AJ294911) and was characterised by a 93-102 nt deletion within the 3A coding region that is associated with the atypical porcinophilic phenotype of this FMDV lineage . Subsequently, O CATHAY isolates have been confirmed in several Southeast and East Asian countries (including Malaysia, the Philippines, Taiwan, Thailand and Vietnam), although since 1970, the majority of field cases due to this topotype have been reported in Hong Kong SAR and China -. The O CATHAY FMD outbreak in Taiwan which began during 1997 resulted in the stamping-out of more than 4 million pigs and generated economic losses of over 6 billion US dollars . Outside of Asia, viruses belonging to the O CATHAY topotype have been responsible for isolated FMD outbreaks that occurred in Europe in 1981 (Thalheim, Austria), 1982 (Wuppertal, Germany) and 1995 (Moscow, Russia). In the last ten years, O CATHAY FMDV strains causing epizootics have been collected in Hong Kong SAR on a yearly basis, where the last reported outbreak occurred during March 2014. However, FMD viruses belonging to the type O CATHAY topotype are sampled on a more sporadic basis from countries in Southeast Asia, and it is currently unclear where this topotype is maintained and/or how it is dispersed.
FMDV in the Philippines
The introduction of FMD into the Philippines can be dated back to 1902 as a result of the importation of infected cattle from Hong Kong SAR to Manila. Following large epidemics reported in Sorsogon and Bukidnon Provinces in 1920, FMD became widespread in the entire Philippines. FMDV lineages belonging to serotypes A, O and C were identified in samples collected from outbreaks occurring in the Philippines during the period between 1954 and 2005. Major epidemics were caused by type O (from 1972 to 1991), type A (from 1975 to 1983) and type C (from 1976 to 1995) strains . The O CATHAY topotype was first detected in August 1994 in a backyard piggery located in Rizal Province. More recently, this FMDV topotype has been the sole lineage responsible for epidemics in the Philippines until December 2005, when the last detected case was confirmed in Quezon Province. The majority of the cases due to O CATHAY were located on Luzon Island, from where FMD spread to 27 provinces. It has been estimated that wholesale market prices of both pork and even chicken in Central Luzon dropped significantly following the start of the epidemic in 1995, highlighting the economic impact of FMD across the entire supply chain . Since June 2011, the Philippines have been officially declared as FMD-free (without vaccination).
This study explored the phylodynamics of these O CATHAY outbreaks reconstructed through molecular epidemiological analyses of VP1 coding sequences (n = 112) collected between 1994 and 2005. In addition, a wider picture of the O CATHAY topotype phylogenetics was determined from a larger database of currently available VP1 coding sequences (n = 322) to enable the characterisation of geographical movements of this FMDV lineage across historically affected countries of Southeast and East Asia.
Materials and methods
This study accessed archived vesicular fluid and/or epithelium samples (n = 112) from the FAO World Reference Laboratory for FMD (WRLFMD) at The Pirbright Institute, United Kingdom, which had been stored at −20°C in 0.04 M phosphate buffer (M25; disodium hydrogen phosphate, potassium dihydrogen phosphate, pH 7.5) and 50% (vol/vol) glycerol. This dataset represented clinical samples collected in the Philippines from 22 provinces in the period between 1994 and 2005 (Additional file 1). In addition, a further 210 VP1 coding region sequences and representing isolates collected from Austria, China, Germany, Hong Kong SAR, Malaysia, Russia, Taiwan, Thailand and Vietnam ,- were retrieved from both GenBank at NCBI  and the WRLFMD sequence archive and, then, integrated with the Philippines collection to comprise a total dataset of 322 VP1 coding sequences (Additional file 2) These VP1 coding region sequences have been submitted to GenBank as have been assigned the following accession numbers: KM243030-KM243172.
Viral RNA detection and sequencing
Oligonucleotide primers used for either RT-PCR or cycle sequencing of the VP1 region from the FMDV isolates
Primer sequence (5' to 3')
Start - end
3558 – 3578
3624 – 3649
3372 – 3390
3181 – 3199
2469 – 2494
2497 – 2519
2508 – 2530
3181 – 3199
2724 – 2740
2830 – 2855
2830 – 2855
Before performing the phylogenetic reconstruction, jModelTest 2.1.4 analysis , was undertaken to determine the best fitting nucleotide substitution model using the Bayesian Information Criterion (BIC) . Statistical parsimony  was used for reconstructing the genealogical networks as implemented in the TCS 1.21 program . The network generated was then edited and plotted in yEd Graph Editor 3.12.
A Bayesian analysis framework was employed for phylogenetic and demographic inferences using a Markov chain Monte Carlo (MCMC) method implemented in the BEAST 1.8.0 package . The analysis was performed using the Hasegawa-Kishino-Yano substitution model plus gamma-distributed rates (HKY85 + Γ4), and the relaxed uncorrelated lognormal molecular clock model ,. Demographic reconstruction was employed using the Bayesian skyline model . Spatial patterns of FMDV dispersal were estimated through a probabilistic discrete asymmetric diffusion model using a continuous-time Markov chain process, adopting a Bayesian stochastic search variable selection (BSSVS) procedure to select among all possible migration pathways . Nonzero rates of virus movement between countries were judged to be supported when the associated Bayes factor (BF) exceeded 3. The MCMCs were run for 150 million iterations, sub-sampling every 15 000 states. Convergence of the chain was assessed using Tracer 1.5 removing the initial 10% of the chain as burn-in. The maximum clade credibility (MCC) tree was summarised using TreeAnnotator 1.8.0 and constructed using FigTree 1.4.0. Phylogeographic maps were constructed using ArcGIS 10.2.1 (Environmental Systems Research Institute, Inc.).
The epidemic curve was constructed using the Handistatus II data for the Philippines retrieved from the OIE website . Statistical computations were performed in R 3.0.3  and graphs were plotted using the ggplot2 package for R , whereas complex vector images were rendered using Inkscape 0.48.4. To determine the potential extent of recombination in the genetic structuring of the virus population, ratios of per-site recombination rate to the per-site mutation rate (r) were estimated using LAMARC 2.1.9 .
O CATHAY FMDV country based phylodynamics: the Philippines
A FASTA search  of all publically available VP1 coding sequences was completed to identify a candidate for the most likely common ancestor for the Philippines lineage: the closest match was identified as a sequence from Hong Kong SAR with 99.2% nt identity (HKN/12/91, GenBank accession no. AJ294921).
Genetic, time and geographical pairwise distances (with corresponding standard deviation values) calculated for the within-year Philippines O CATHAY FMDV isolates groups and for each of the country based data from the earliest samples collected within the specific group
No of samples
0.010 ± 0.002
0.23 ± 0.05
0.013 ± 0.006
0.63 ± 0.77
124.72 ± 95.81
0.020 ± 0.004
0.10 ± 0.03
235.66 ± 252.21
0.016 ± 0.013
0.38 ± 0.38
119.65 ± 117.12
0.011 ± 0.009
51.62 ± 42.08
0.011 ± 0.002
146.83 ± 222.79
0.054 ± 0.021
0.16 ± 0.13
296.55 ± 130.74
0.057 ± 0.008
14.22 ± 4.28
0.051 ± 0.013
0.19 ± 0.08
344.45 ± 26.98
0.010 ± 0.004
0.41 ± 0.23
322.32 ± 38.83
0.005 ± 0.005
0.03 ± 0.01
12.48 ± 17.65
0.148 ± 0.081
31.93 ± 19.15
0.157 ± 0.022
32.95 ± 7.24
0.047 ± 0.022
4.59 ± 3.50
0.015 ± 0.025
1.38 ± 3.01
0.104 ± 0.016
7.84 ± 1.65
The molecular clock for the O CATHAY Philippines lineage was estimated to be 1.25 × 10−2 nt/site/yr (95%HPD 9.47 × 10−3 to 1.57 × 10−2) with a standard deviation of 0.70 (95%HPD 0.49 to 0.91). No evidence of autocorrelation of rates in the reconstructed phylogeny was provided by the covariance value of 2.65 × 10−3. The introduction date, the time of the MRCA (TMRCA), of the type O CATHAY topotype FMDV lineage into the Philippines was calculated to be the 30th of March 1994 (95%HPD 07/08/1993 to 08/08/1994, a time interval which included the date of the first officially reported case).
The reconstructed FMDV population dynamics from the skyline plot (Figure 3) describes a demographic history characterised by three phases. In the first phase, after an initial exponential increase from mid-1994 until late 1997 at a rate that decreased from late 1996, a sudden and short period of decline was observed, resulting in a population bottleneck. Since genetic bottlenecks correspond to significant reductions in population size, these changes in the O CATHAY population dynamics in the Philippines probably link to the launch of an extensive control plan in 1996 that was successful in limiting the further spread of FMD and thereby reducing the number of outbreaks . However, during 1999 a new FMD outbreak occurred within an already declared FMD-free zone, the Panay region. Therefore in the second phase, the skyline trajectory recorded a second rapidly increasing viral population size starting in mid-1998 and lasting up to the first months of 1999, which resulted in a diversification of viral lineages. In the third phase, the viral population size reached a plateau until late 2002, when further control policies resulted in a steady decline in FMD prevalence until eradication.
The epidemic curve drawn from the field epidemiological data from the OIE for the period 1995-2005  described an oscillatory trend in the number of FMD outbreaks reported in the Philippines, with times of high epidemic peaks interleaved by low-level FMD circulation. The frequency of these oscillations was higher between 1997 and the beginning of 2000 (a monthly average of 37.9 FMD outbreaks), after which the number of FMD outbreaks started to decline following periods of low reporting (with a monthly average of 19.8 FMD outbreaks). However, the reported epidemic trend did not overlap with the skyline plot trajectory, although the epidemic window from mid-2000 to 2005 characterised by a reduced number of outbreaks could be evinced by the plateauing and subsequent decrease in the genetic diversity of the skyline plot. It should be noted that although more than 300 outbreaks were officially reported through OIE during 2002, no clinical samples (and thus genetic information) were collected within that time window.
According to the results obtained by the discrete states phylogeography analysis, the root of the Philippines taxon was found to be from Rizal Province, consistent with the location of the first officially reported cases of O CATHAY topotype in the Philippines during August 1994 (Figure 3). Three main epidemic hubs could be identified from the analysis: the first from the beginning of the epizootic up to mid-1996, where outbreaks were found to be seeded from Rizal Province; the second lasting until 2001, where Bulacan Province was estimated to be the main source of FMD spread; and lastly, Manila Province as the last epidemic hub. The movement transitions between the three main epidemic hubs were supported by Bayes factor values of > 24 [posterior probability (pk) = 1.0] for movements from Rizal to Bulacan and from Bulacan to Manila, respectively.
Global and regional phylodynamics of O CATHAY topotype FMDV
The molecular clock rate for all the O CATHAY topotype VP1 data was estimated to be 1.06 × 10−2 nt/site/yr (95%HPD 8.99 × 10−3 to 1.23 × 10−2), with a standard deviation of 0.81 (95% HPD 0.67 to 0.94). This value was comparable with the molecular clock rate reported for the Philippine isolates only. The MRCA for the O CATHAY topotype was estimated to have been present between 1955 and 1960. The r recombination parameter returned a value of 8.3 × 10−9 site/generation indicating a very low influence of recombination relative to mutation.
The MRCA shared between the Philippines and Taiwanese phylogenetic clusters was estimated to have been circulating in 1993 (95%HPD 1992 to 1994), whereas the origin of the MRCA for the more recent O CATHAY FMD epidemics in the Southeast and East Asia regions was dated 1991 (95%HPD 1990 to 1992). No other virus introduction or escape was ascribed to the Philippines O CATHAY FMD epidemic history, suggesting the Philippines sub-lineage to be monophyletic. In contrast, Hui and Leung  described two different FMDV introductions into the Philippines, assigning the PHI/5/95 isolates within the phylogenetic cluster which includes the Taiwanese isolates. However, the tree node that governed this inclusion had a bootstrap value of < 70, suggesting uncertainty in the assignment of these descendants.
The evolutionary dynamics of the O CATHAY topotype of FMDV have been analysed allowing the transmission dynamics to be reconstructed across countries in Southeast Asia that have been impacted by this lineage. The O CATHAY FMDV strains isolated from outbreaks reported in Hong Kong, Taiwan and Philippines were defined as belonging to three different sublineages, which were related by a shared common ancestry to an unsampled FMDV strain sourced from Hong Kong SAR. The O CATHAY FMD epizootic in the Philippines resulted from a single introduction and was characterised by three main transmission hubs in Rizal, Bulacan and Manila. Although the evolutionary dynamics of the O CATHAY FMDV lineage were described by three phases from the skyline reconstruction, this was not entirely consistent with the monthly epidemic curve (Figure 3). This could be either due to a spatio-temporal bias in the genetic information analysed or in the incompleteness of the outbreak reporting database used, or both.
The phylodynamics of FMDV reconstructed from the FMDV type O CATHAY VP1 coding sequences indicates a marked reduction in viral diversity in the last 10 years, corresponding to the eradication of FMD in the Philippines and the more limited disease events experienced in Taiwan. Furthermore, the introduction of the FMDV type O Southeast Asia (SEA) topotype Mya-98 lineage into Hong Kong SAR during 2010 could have reduced the genetic diversity within O CATHAY lineages through direct competition with available hosts, as well as the presence of cross-protective antibodies in convalescent animals. These findings indicate that the O CATHAY topotype is maintained in the Hong Kong SAR ecosystem and sporadically spread from there to other Southeast Asian countries, as would be the case for the Philippines in 1994 and Vietnam in 1997. However, few O CATHAY FMDV strains have been reported from mainland China, which has the largest swine production industry in the world (representing over 51% of the world's pig population). These few isolates were collected in 1986, 2000, 2001 and 2003, therefore sampling bias or underreporting of epidemic events occurring in China would likely have an impact on assessing the geographical movements of the FMDV type O CATHAY topotype. It is, nevertheless, clear from the analysis that a transmission link exists between China and Hong Kong SAR, thus indicating a historically southward movement of the O CATHAY FMDV lineage.
The molecular clock estimated here for the O CATHAY topotype is at the high end of evolutionary rate estimates for FMDV. Previously estimates reported an average evolutionary rate across all FMDV serotype of 2.48 × 10−3 nt/site/yr , while rates of 3.14 × 10−3, 1.3 × 10−3 and 4.8 × 10−3 nt/site/yr were reported for serotype O -. In addition, lineage-based FMDV molecular clock rates of 2.8 × 10−3, 6.65 × 10−3, 7.81 × 10−3 and 2.7 × 10−3 nt/site/yr were previously estimated for the O-PanAsia lineage in India, O-PanAsia-2 sublineage in Pakistan and Afghanistan, and type O in East Africa, respectively -. The higher rate of FMDV evolution reported for the A-Iran 05 FMDV lineage in Afghanistan and Pakistan (1.2 ± 10−2 nt/site/yr)  was similar to the molecular clock for the O CATHAY topotype estimated by this study. Therefore, genotypically and regionally variable evolutionary rates may in fact reflect real differences in the epidemiological dynamics and host-interaction of FMDV.
Although using a large database of FMDV isolates and generating a comprehensive picture of the O CATHAY topotype evolutionary history, this study has some limitations largely derived from the nature of the genetic data used for the analysis. The VP1 coding region, although defining only ~8% (639 nt of length) of the complete FMDV genome, is the most variable section of the FMDV genome and is historically used for tracing the movement and spread of FMD globally , and, furthermore, provides the basis for FMDV genotype definition . Analysing a larger part of the FMDV genome, such as the whole capsid region or the full-length genome, would produce results with a higher resolution ,. However, it should be noted that recombination events seem to be more widespread in other part of the genome ,,, thus representing a limitation in interpreting results based on full-length genome analysis of large scale FMDV evolutionary studies. The ratio of per-site recombination to mutation rate here estimated from the full currently available FMDV type O CATHAY topotype VP1 coding sequences database is very low indicating that these results are not influenced by the process of recombination.
This study was supported by the UK Biotechnology and Biological Science Research Council (BBSRC) Doctoral Training Award (DTA) granted to ADN. DPK, NJK and JW are partially funded by a grant from the UK Department for Environment, Food and Rural Affairs (DEFRA SE2940) that supports molecular characterisation work at The Pirbright Institute under the auspices of the WRLFMD, led by DPK. The authors would like to thank the field teams in the Philippines for their effort in collecting the clinical samples and the laboratory staff at the Pirbright Institute for all their help and maintenance of the WRLFMD sample archive.
- Knowles NJ, Samuel AR: Molecular epidemiology of foot-and-mouth disease virus. Virus Res. 2003, 91: 65-80. 10.1016/S0168-1702(02)00260-5.View ArticlePubMedGoogle Scholar
- Knowles NJ, Wadsworth J, Parlak U, Ozyoruk F, Nazem Shirazi MH, Ferris NP, Hutchings GH, Stirling JM, Hammond JM, King DP: Recent events in the evolution of foot-and-mouth disease in the Middle East. Open Session of the European Commission for the Control of Foot-and-Mouth Disease Standing Technical Committee: Vienna, 28 September - 1 October 2010. 2010, Food and Agriculture Organization of the United Nations, Rome, ItalyGoogle Scholar
- Di Nardo A, Knowles NJ, Paton DJ: Combining livestock trade patterns with phylogenetics to help understand the spread of foot and mouth disease in sub-Saharan Africa, the Middle East and Southeast Asia. Rev Sci Tech. 2011, 30: 63-85.PubMedGoogle Scholar
- Samuel AR, Knowles NJ: Foot-and-mouth disease type O viruses exhibit genetically and geographically distinct evolutionary lineages (topotypes). The J Gen Virol. 2001, 82: 609-621.View ArticlePubMedGoogle Scholar
- Knowles NJ, Wadsworth J, Hammond JM, King DP: Foot-and-mouth disease virus genotype definitions and nomenclature. Open Session of the European Commission for the Control of Foot-and-Mouth Disease Standing Technical Committee: Vienna, 28 September - 1 October 2010. 2010, Food and Agriculture Organization of the United Nations, Rome, ItalyGoogle Scholar
- Knowles NJ, Davies PR, Henry T, O'Donnell V, Pacheco JM, Mason PW: Emergence in Asia of foot-and-mouth disease viruses with altered host range: characterization of alterations in the 3A protein. J Virol. 2001, 75: 1551-1556. 10.1128/JVI.75.3.1551-1556.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Gleeson LJ: A review of the status of foot and mouth disease in South-East Asia and approaches to control and eradication. Rev Sci Tech. 2002, 21: 465-475.PubMedGoogle Scholar
- Knowles NJ, Samuel AR, Davies PR, Midgley RJ, Valarcher JF: Pandemic strain of foot-and-mouth disease virus serotype O. Emerg Infect Dis. 2005, 11: 1887-1893. 10.3201/eid1112.050908.PubMed CentralView ArticlePubMedGoogle Scholar
- Cao Y, Lu Z, Li D, Fan P, Sun P, Bao H, Fu Y, Li P, Bai X, Chen Y, Xie B, Liu Z: Evaluation of cross-protection against three topotypes of serotype O foot-and-mouth disease virus in pigs vaccinated with multi-epitope protein vaccine incorporated with poly(I:C). Vet Microbiol. 2014, 168: 294-301. 10.1016/j.vetmic.2013.11.023.View ArticlePubMedGoogle Scholar
- Yang PC, Chu RM, Chung WB, Sung HT: Epidemiological characteristics and financial costs of the 1997 foot-and-mouth disease epidemic in Taiwan. Vet Rec. 1999, 145: 731-734.PubMedGoogle Scholar
- Randolph TF, Perry BD, Benigno CC, Santos IJ, Agbayani AL, Coleman P, Webb R, Gleeson LJ: The economic impact of foot and mouth disease control and eradication in the Philippines. Rev Sci Tech. 2002, 21: 645-661.PubMedGoogle Scholar
- Abao LNB, Kono H, Gunarathne A, Promentilla RR, Gaerlan MZ: Impact of foot-and-mouth disease on pork and chicken prices in Central Luzon, Philippines. Prev Vet Med. 2014, 113: 398-406. 10.1016/j.prevetmed.2013.12.005.View ArticlePubMedGoogle Scholar
- Beard CW, Mason PW: Genetic determinants of altered virulence of Taiwanese foot-and-mouth disease virus. J Virol. 2000, 74: 987-991. 10.1128/JVI.74.2.987-991.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsai CP, Pan CH, Liu MY, Lin YL, Chen CM, Huang TS, Cheng IC, Jong MH, Yang PC: Molecular epidemiological studies on foot-and-mouth disease type O Taiwan viruses from the 1997 epidemic. Vet Microbiol. 2000, 74: 207-216. 10.1016/S0378-1135(00)00182-6.View ArticlePubMedGoogle Scholar
- Knowles NJ, Samuel AR, Davies PR, Kitching RP, Donaldson AI: Outbreak of foot-and-mouth disease virus serotype O in the UK caused by a pandemic strain. Vet Rec. 2001, 148: 258-259.PubMedGoogle Scholar
- 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. 10.1128/JVI.79.10.6487-6504.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Abdul-Hamid NF, Hussein NM, Wadsworth J, Radford AD, Knowles NJ, King DP: Phylogeography of foot-and-mouth disease virus types O and A in Malaysia and surrounding countries. Infect Genet Evol. 2011, 11: 320-328. 10.1016/j.meegid.2010.11.003.View ArticlePubMedGoogle Scholar
- Hui RK, Leung FC: Evolutionary trend of foot-and-mouth disease virus in Hong Kong. Vet Microbiol. 2012, 159: 221-229. 10.1016/j.vetmic.2012.03.016.View ArticlePubMedGoogle Scholar
- Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW: GenBank. Nucleic Acids Res. 2013, 41: D36-D42. 10.1093/nar/gks1195. [http://www.ncbi.nlm.nih.gov/genbank/], [http://www.ncbi.nlm.nih.gov/genbank/]PubMed CentralView ArticlePubMedGoogle Scholar
- Knowles NJ, Nazem Shirazi MH, Wadsworth J, Swabey KG, Stirling JM, Statham RJ, Li Y, Hutchings GH, Ferris NP, Parlak U, Ozyörük F, Sumption KJ, King DP, Paton DJ: Recent spread of a new strain (A-Iran-05) of foot-and-mouth disease virus type A in the Middle East. Transbound Emerg Dis. 2009, 56: 157-169. 10.1111/j.1865-1682.2009.01074.x.View ArticlePubMedGoogle Scholar
- Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Soding J, Thompson JD, Higgins DG: Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011, 7: 539-10.1038/msb.2011.75.PubMed CentralView ArticlePubMedGoogle Scholar
- Forss S, Strebel K, Beck E, Schaller H: Nucleotide sequence and genome organization of foot-and-mouth disease virus. Nucleic Acids Res. 1984, 12: 6587-6601. 10.1093/nar/12.16.6587.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D: jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012, 9: 772-10.1038/nmeth.2109.PubMed CentralView ArticlePubMedGoogle Scholar
- Posada D, Buckley TR: Model selection and model averaging in phylogenetics: advantages of akaike information criterion and bayesian approaches over likelihood ratio tests. Syst Biol. 2004, 53: 793-808. 10.1080/10635150490522304.View ArticlePubMedGoogle Scholar
- Templeton AR, Crandall KA, Sing CF: A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data: III: cladogram estimation. Genetics. 1992, 132: 619-633.PubMed CentralPubMedGoogle Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.View ArticlePubMedGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.PubMed CentralView ArticlePubMedGoogle Scholar
- Hasegawa M, Kishino H, Yano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985, 22: 160-174. 10.1007/BF02101694.View ArticlePubMedGoogle Scholar
- Drummond AJ, Ho SY, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: e88-10.1371/journal.pbio.0040088.PubMed CentralView ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005, 22: 1185-1192. 10.1093/molbev/msi103.View ArticlePubMedGoogle Scholar
- Lemey P, Rambaut A, Drummond AJ, Suchard MA: Bayesian phylogeography finds its roots. PLoS Comput Biol. 2009, 5: e1000520-10.1371/journal.pcbi.1000520.PubMed CentralView ArticlePubMedGoogle Scholar
- World Organisation for Animal Health: Animal health data (Handistatus II). 2014 , [http://www.oie.int/hs2/report.asp]
- R: A Language and Environment for Statistical Computing. 2013, R Foundation for Statistical Computing, Vienna, AustriaGoogle Scholar
- Wickham H: Ggplot2: Elegant Graphics for Data Analysis. 2009, Springer, New YorkView ArticleGoogle Scholar
- Kuhner MK: LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics. 2006, 22: 768-770. 10.1093/bioinformatics/btk051.View ArticlePubMedGoogle Scholar
- McWilliam H, Li W, Uludag M, Squizzato S, Park YM, Buso N, Cowley AP, Lopez R: Analysis tool web services from the EMBL-EBI. Nucleic Acids Res. 2013, 41: W597-W600. 10.1093/nar/gkt376. [http://www.ebi.ac.uk/Tools/sss/fasta/nucleotide.html], [http://www.ebi.ac.uk/Tools/sss/fasta/nucleotide.html]PubMed CentralView ArticlePubMedGoogle Scholar
- Sinnott RW: Virtues of the haversine. Sky Telescope. 1984, 68: 159-Google Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59: 307-321. 10.1093/sysbio/syq010.View ArticlePubMedGoogle Scholar
- Tully DC, Fares MA: The tale of a modern animal plague: tracing the evolutionary history and determining the time-scale for foot-and-mouth disease virus. Virology. 2008, 382: 250-256. 10.1016/j.virol.2008.09.011.View ArticlePubMedGoogle Scholar
- Jenkins GM, Rambaut A, Pybus OG, Holmes EC: Rates of molecular evolution in RNA viruses: a quantitative phylogenetic analysis. J Mol Evol. 2002, 54: 156-165. 10.1007/s00239-001-0064-3.View ArticlePubMedGoogle Scholar
- Yoon SH, Lee KN, Park JH, Kim H: Molecular epidemiology of foot-and-mouth disease virus serotypes A and O with emphasis on Korean isolates: temporal and spatial dynamics. Arch Virol. 2011, 156: 817-826. 10.1007/s00705-011-0921-3.View ArticlePubMedGoogle Scholar
- Hemadri D, Tosh C, Sanyal A, Venkataramanan R: Emergence of a new strain of type O foot-and-mouth disease virus: its phylogenetic and evolutionary relationship with the PanAsia pandemic strain. Virus Genes. 2002, 25: 23-34. 10.1023/A:1020165923805.View ArticlePubMedGoogle Scholar
- Jamal SM, Ferrari G, Ahmed S, Normann P, Belsham GJ: Genetic diversity of foot-and-mouth disease virus serotype O in Pakistan and Afghanistan, 1997-2009. Infect Genet Evol. 2011, 11: 1229-1238. 10.1016/j.meegid.2011.03.006.View ArticlePubMedGoogle Scholar
- Balinda SN, Sangula AK, Heller R, Muwanika VB, Belsham GJ, Masembe C, Siegismund HR: Diversity and transboundary mobility of serotype O foot-and-mouth disease virus in East Africa: implications for vaccination policies. Infect Genet Evol. 2010, 10: 1058-1065. 10.1016/j.meegid.2010.06.017.View ArticlePubMedGoogle Scholar
- Jamal SM, Ferrari G, Ahmed S, Normann P, Curry S, Belsham GJ: Evolutionary analysis of serotype A foot-and-mouth disease viruses circulating in Pakistan and Afghanistan during 2002-2009. J Gen Virol. 2011, 92: 2849-2864. 10.1099/vir.0.035626-0.View ArticlePubMedGoogle Scholar
- Cottam EM, Haydon DT, Paton DJ, Gloster J, Wilesmith JW, Ferris NP, Hutchings GH, King DP: Molecular epidemiology of the foot-and-mouth disease virus outbreak in United Kingdom in 2001. J Virol. 2006, 80: 11274-11282. 10.1128/JVI.01236-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Cottam EM, Wadsworth J, Shaw AE, Rowlands RJ, Goatley L, Maan S, Maan NS, Mertens PP, Ebert K, Li Y, Ryan ED, Juleff N, Ferris NP, Wilesmith JW, Haydon DT, King DP, Paton DJ, Knowles NJ: Transmission pathways of foot-and-mouth disease virus in the United Kingdom in 2007. PLoS Pathog. 2008, 4: e1000050-10.1371/journal.ppat.1000050.PubMed CentralView ArticlePubMedGoogle Scholar
- Jackson AL, O’Neill H, Maree F, Blignaut B, Carrillo C, Rodriguez L, Haydon DT: Mosaic structure of foot-and-mouth disease virus. J Gen Virol. 2007, 88: 487-492. 10.1099/vir.0.82555-0.View ArticlePubMedGoogle Scholar
- Wright CF, Knowles NJ, Di Nardo A, Paton DJ, Haydon DT, King DP: Reconstructing the origin and transmission dynamics of the 1967-68 foot-and-mouth disease epidemic in the United Kingdom. Infect Genet Evol. 2013, 20: 230-238. 10.1016/j.meegid.2013.09.009.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin YL, Jong MH, Huang CC, Shieh HK, Chang PC: Genetic and antigenic characterization of foot-and-mouth disease viruses isolated in Taiwan between 1998 and 2009. Vet Microbiol. 2010, 145: 34-40. 10.1016/j.vetmic.2010.03.003.View ArticlePubMedGoogle Scholar
- Feng Q, Yu H, Liu Y, He C, Hu J, Sang H, Ding N, Ding M, Fung YW, Lau LT, Yu AC, Chen J: Genome comparison of a novel foot-and-mouth disease virus with other FMDV strains. Biochem Biophys Res Commun. 2004, 323: 254-263. 10.1016/j.bbrc.2004.08.086.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.