Volume 24, Number 6—June 2018
Foot-and-Mouth Disease in the Middle East Caused by an A/ASIA/G-VII Virus Lineage, 2015–2016
Foot-and-mouth disease (FMD) can decrease productivity in the cloven-hooved livestock industry. As this disease spreads rapidly over large distances, it is regarded to be one of the most economically devastating diseases of livestock. FMD is caused by FMD virus (FMDV; family Picornaviridae, genus Aphthovirus), which has 7 immunologically distinct serotypes, O, A, C, Asia 1, SAT 1, SAT 2, and SAT 3. Worldwide, ecologic niches of FMDV circulation have been defined as 7 virus pools. Pools 1–3 are present in Asia, where only serotypes O, A, and Asia 1 are present.
Serotype A viruses are considered to be the most variable (genetically and antigenically) of EurAsian serotypes. Three topotypes (ASIA, AFRICA, and EURO-SA [Europe–South America]) and multiple diverse lineages and sublineages have been identified (1). The ASIA topotype is widespread and is found in most countries in Asia; there have been sporadic incursions into North Africa. Although the G-VII lineage (also known as genotype 18) (2,3) usually circulates in countries containing virus pool 2 (commonly in Bangladesh and India, rarely in Bhutan and Nepal, but until now not in Sri Lanka), this lineage has also been reported in Saudi Arabia in 1995, Albania and the former Yugoslav Republic of Macedonia in 1996, and Myanmar in 2010.
In India, viruses of the A/ASIA/G-VII lineage were isolated in 1983 (4) and until 2001 were co-circulating with the A/ASIA/G-VI lineage (2). After 2001, only the G-VII lineage has been reported (5). Overall incidence of FMD outbreaks caused by serotype A in India during 2011–2016 was low (3.1% of the total reported outbreaks). During the same time, 86.8% of outbreaks were caused by serotype O, and 10.1% by serotype Asia 1 (6).
Despite the low number of outbreaks investigated, emergence of a capsid viral protein (VP3) 3 deletion variant, VP59, was reported in the A/ASIA/G-VII lineage during 2002 (4). Currently, a group of viruses described as clade C, a subgroup within the VP359-deletion variant, is speculated to be the dominant group of the A/ASIA/G-VII lineage prevalent in India (2). We report foot-and-mouth disease in the Middle East during 2015–2016 caused by an A/ASIA/G-VII virus lineage.
During outbreak investigations of FMD in Saudi Arabia in 2015, a virus of the A/ASIA/G-VII lineage (VP359-deletion variant), was identified in cattle (7). The outbreak spread quickly to several strictly monitored dairy farms that had high rates of vaccination, as well as to nomadic herds. Concurrently, related viruses were found in Armenia, Iran, and Turkey in 2015 and continued to circulate in Saudi Arabia, Turkey, and Iran in 2016. Sequences of VP1-coding regions from samples submitted to the Food and Agriculture Organization of the United Nations World Reference Laboratory for FMD (Pirbright, UK) were determined by using described methods (8). VP1 sequences from outbreaks in Armenia, Bangladesh, and Turkey were determined at the Federal Centre for Animal Health (Vladimir, Russia); the University of Dhaka (Dhaka, Bangladesh); and the Food-and-Mouth Disease Institute (Ankara, Turkey), respectively (Table).
We performed maximum-likelihood analyses to compare VP1 coding sequences with other contemporary sequences of the A/ASIA/G-VII lineage and grouped them within the VP359-deleted C clade (2) (Figure 1, panel A). We estimated time-resolved phylogenetic trees for 101 serotype A FMDV G-VII VP1 sequences by using BEAST version 1.8.4 (9) and incorporated the general time reversible model with gamma-distributed rate variation among sites and 0.5 prior proportion of invariant sites, the Bayesian Skyline tree before accounting for demographic uncertainty, and a log-normal uncorrelated relaxed clock across branches (10). We ran Markov Chain Monte Carlo analysis for 200 million steps and sampled trees every 20,000 steps after a burn-in of 20 million steps. We assessed convergence and good mixing of the Markov Chain Monte Carlo C chain by using Tracer version 1.6 (http://beast.community/tracer).
VP1 coding region–based Bayesian analyses identified >2 independent introductions of G-VII virus into the study region, >1 to Saudi Arabia and 1 to Iran (Figure 1, panel B). However, lack of sequences for recent viruses circulating in the Indian subcontinent makes it difficult to resolve more precisely the number of introductions. Closely related viruses might be circulating in a wider geographic area, thus being a source of the outbreaks. However, this speculation is not supported by available epidemiologic information.
Outbreaks in Armenia, Iran, and Turkey were closely related and most likely originated from the same source. The most recent common ancestor (MRCA) was dated to March 2015 (95% high posterior density [HPD] October 2014–July 2015). The MRCA of the C clade was dated to January 2006 (95% HPD June 2005–April 2006), which is consistent with the first isolate obtained in India during 2007 (2). The MRCA of the phylogenetic cluster grouping Middle East (Saudi Arabia, Iran, and Turkey) isolates was dated to July 2015 (95% HPD April–October 2015). Potential movement of G-VII FMDV lineages from the Indian subcontinent might be dated to October 2014 (95% HPD January 2014–February 2015). Evolution of the G-VII C clade lineage was estimated as having a mutation rate of 1.1 × 10−2 nt/site/y (95% HPD 8.0 × 10−3 to 1.4 × 10−2 nt/site/y).
Amino acid substitutions on the FMDV surface, particularly in the G-H loop of VP1 (antigenic site 1), have been implicated in antigenic variation of the virus in vitro and in vivo (11). Predicted amino acid sequences obtained from samples collected during outbreaks in 2015–2016 were compared with the A/BAN/CH/Sa-304/2016 sequence (most closely related virus from the Indian subcontinent), and 8 nonconservative substitutions were identified at 6 positions (138, 143, 145, 147, 148, and 151) within the VP1 antigenic site 1 (12,13) (Figure 2; Table).
Changes from hydrophobic (alanine and leucine) to hydrophilic (threonine, glutamine, and serine) amino acid residues were most common, found at 4 positions; changes from hydrophobic (glycine) to acidic (aspartic acid and glutamic acid) amino acids were found at 2 positions. Antigenic variation at site 1 of type C viruses is often based on alternate switching between alanine and threonine residues without accumulation of amino acid substitutions (14). In addition, 2 independent changes that did not alter the amino acid characteristics were identified at positions 142 and 147. We also showed that changes at antigenic site 1 were conserved mainly within but differed between farms, supporting independent selection pressures. Although the same vaccine was used, the intensive and frequent vaccination regimen routinely used on the affected farms in Saudi Arabia might have led to an independent antigenic evolution on an individual farm level from chance substitutions. Nevertheless, occasional substitutions within antigenic site 1 were also observed in Armenia, Turkey, and Iran (Figure 2).
As reported for the O/ME-SA/Ind-2001d virus lineage (15), A/ASIA/G-VII is the second FMDV lineage believed to have originated in Indian subcontinent since 2013 and resulted in extensive outbreaks outside its usual area of distribution. Similar to the Ind-2001d outbreaks, current outbreaks caused by the G-VII lineage appear to be linked to multiple introductions of the virus from the Indian subcontinent; the virus then spread among susceptible ruminant populations in Saudi Arabia, Iran, Turkey, and Armenia.
It is a concern that in vitro vaccine matching data (by virus neutralization) provides poor confidence that commercially available vaccines would offer effective protection against the G-VII lineage (A. Ludi, pers. comm., May 2016). To improve control programs, it is crucial to identify expected routes of FMDV escape (e.g., international trade in animals and animal products) outside historically defined geographic distribution, and to establish transmission pathways within affected areas. To reconstruct likely transmission pathways at greater resolution, genome sequencing of viruses described in this report is currently in progress.
Dr. Bachanek-Bankowska is a senior research scientist at the Food and Agriculture Organization of the United Nations World Reference Laboratory for FMD at The Pirbright Institute, United Kingdom. Her primary research interests are molecular epidemiology and development of molecular diagnostic methods for FMDV.
This study was supported by the Department for Environment, Food and Rural Affairs, United Kingdom (research grant SE2943). The World Reference Laboratory for Foot-and-Mouth Disease was supported by the European Union (contract from the European Commission for the Control of Foot-and-Mouth Disease, Rome, Italy). The Pirbright Institute was supported by a grant from the Biotechnology and Biological Sciences Research Council of the United Kingdom.
- Knowles NJ, Samuel AR. Molecular epidemiology of foot-and-mouth disease virus. Virus Res. 2003;91:65–80.
- Das B, Mohapatra JK, Pande V, Subramaniam S, Sanyal A. Evolution of foot-and-mouth disease virus serotype A capsid coding (P1) region on a timescale of three decades in an endemic context. Infect Genet Evol. 2016;41:36–46.
- Nandi SP, Rahman MZ, Momtaz S, Sultana M, Hossain MA. Emergence and distribution of foot-and-mouth disease virus serotype A and O in Bangladesh. Transbound Emerg Dis. 2015;62:328–31.
- Jangra RK, Tosh C, Sanyal A, Hemadri D, Bandyopadhyay SK. Antigenic and genetic analyses of foot-and-mouth disease virus type A isolates for selection of candidate vaccine strain reveals emergence of a variant virus that is responsible for most recent outbreaks in India. Virus Res. 2005;112:52–9.
- Mohapatra JK, Subramaniam S, Pandey LK, Pawar SS, De A, Das B, et al. Phylogenetic structure of serotype A foot-and-mouth disease virus: global diversity and the Indian perspective. J Gen Virol. 2011;92:873–9.
- Indian Council of Agricultural Research–Directorate on Foot and Mouth Disease. Annual report 2015–2016 [cited 2018 Feb 5]. http://www.pdfmd.ernet.in/index_files/Content/Reports/PDFMD_AR_2015_16_English.pdf
- Bachanek-Bankowska K, Wadsworth J, Thapa B, King DP, Knowles NJ. Complete genome sequence of a serotype A foot-and-mouth disease virus from an outbreak in Saudi Arabia during 2015. Genome Announc. 2016;4:e01591–15.
- Knowles NJ, Wadsworth J, Bachanek-Bankowska K, King DP. VP1 sequencing protocol for foot and mouth disease virus molecular epidemiology. Rev Sci Tech. 2016;35:741–55.
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
- Drummond AJ, Ho SY, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4:e88.
- Grubman MJ, Baxt B. Foot-and-mouth disease. Clin Microbiol Rev. 2004;17:465–93.
- Thomas AA, Woortmeijer RJ, Puijk W, Barteling SJ. Antigenic sites on foot-and-mouth disease virus type A10. J Virol. 1988;62:2782–9.
- Bolwell C, Clarke BE, Parry NR, Ouldridge EJ, Brown F, Rowlands DJ. Epitope mapping of foot-and-mouth disease virus with neutralizing monoclonal antibodies. J Gen Virol. 1989;70:59–68.
- Martínez MA, Dopazo J, Hernández J, Mateu MG, Sobrino F, Domingo E, et al. Evolution of the capsid protein genes of foot-and-mouth disease virus: antigenic variation without accumulation of amino acid substitutions over six decades. J Virol. 1992;66:3557–65.
- Knowles NJ, Bachanek-Bankowska K, Wadsworth J, Mioulet V, Valdazo-González B, Eldaghayes IM, et al. Outbreaks of foot-and-mouth disease in Libya and Saudi Arabia during 2013 due to an exotic O/ME-SA/Ind-2001 lineage virus. Transbound Emerg Dis. 2016;63:e431–5.
Suggested citation for this article: Bachanek-Bankowska K, Di Nardo A, Wadsworth J, Henry EKM, Parlak Ü, Timina A, et al. Foot-and-mouth disease in the Middle East caused by an A/ASIA/G-VII virus lineage, 2015–2016. Emerg Infect Dis. 2018 Jun [date cited]. https://doi.org/10.3201/eid2406.170715