Dynamics of widespread foot-and-mouth disease virus serotypes A, O and Asia-1 in southern Asia: A Bayesian phylogenetic perspective

Foot-and-mouth disease (FMD) is, arguably, the animal disease with the most devastating global economic impact owing in part, to the severe trade restrictions imposed upon affected countries and regions. South Asia is one of the regions where widespread lineages of the FMDV virus (FMDV) have emerged. Here, we performed an integrative phylogenetic analysis of all FMDV serotypes (A, O and Asia-1) circulating in southern Asia, including viral sequences collected until 2013. Our results describe the occurrence of FMD caused by different serotypes and lineages, focusing in the cycles where a specific lineage predominates within a region for a protracted period and then are rapidly or progressively replaced by an emergent or re-emergent strain that is introduced from an adjacent region. Transmission between the two main regions in southern Asia (the Indian subcontinent and the region comprised by Afghanistan, Iran and Pakistan) has been limited. Results of time divergence estimation of lineages that currently circulate in this region indicate that the most recent common ancestor of endemic lineages are: 1992 [1989-1995] for lineage O/PanAsia; 1997 [1995-1999] for PanAsia2; 2001 [1998-2004] for O/Ind2001; 2001 [2000-2002] for A/Iran-05; 1990 [1988-1991] for A/G-18 (G-VII); 2003 [2000-2006] for Asia-1 Sindh08 and 2002 [1999-2004] for Asia-1 G-VIII. We estimated the mean of the overall substitution rate of the VP1 coding region (substitution/site/year) for serotype O (5.95 × 10^-3 ), serotype A (1.19 × 10^-2 ) and serotype Asia-1 (3.08 × 10^-3 ). The potential factors driving the lineage turnover are discussed. Our results provide insights into the ecological and evolutionary factors driving the emergence of FMDV.

Foot-and-mouth disease virus transmission dynamics and persistence in a herd of vaccinated dairy cattle in India

Foot-and-mouth disease (FMD) is an important transboundary disease with substantial economic impacts. Although between-herd transmission of the disease has been well studied, studies focusing on within-herd transmission using farm-level outbreak data are rare. The aim of this study was to estimate parameters associated with within-herd transmission, host physiological factors and FMD virus (FMDV) persistence using data collected from an outbreak that occurred at a large, organized dairy farm in India. Of 1,836 regularly vaccinated, adult dairy cattle, 222 had clinical signs of FMD over a 39-day period. Assuming homogenous mixing, a frequency-dependent compartmental model of disease transmission was built. The transmission coefficient and basic reproductive number were estimated to be between 16.2-18.4 and 67-88, respectively. Non-pregnant animals were more likely to manifest clinical signs of FMD as compared to pregnant cattle. Based on oropharyngeal fluid (probang) sampling and FMDV-specific RT-PCR, four of 36 longitudinally sampled animals (14%) were persistently infected carriers 10.5 months post-outbreak. There was no statistical difference between subclinical and clinically infected animals in the duration of the carrier state. However, prevalence of NSP-ELISA antibodies differed significantly between subclinical and clinically infected animals 12 months after the outbreak with 83% seroprevalence amongst clinically infected cattle compared to 69% of subclinical animals. This study further elucidates within-herd FMD transmission dynamics during the acute-phase and characterizes duration of FMDV persistence and seroprevalence of FMD under natural conditions in an endemic setting.

Quantitative characteristics of the foot-and-mouth disease carrier state under natural conditions in India

The goal of this study was to characterize the properties and duration of the foot-and-mouth disease (FMD) carrier state and associated serological responses subsequent to vaccination and naturally occurring infection at two farms in northern India. Despite previous vaccination of cattle in these herds, clinical signs of FMD occurred in October 2013 within a subset of animals at the farms containing juvenile-yearling heifers and steers (Farm A) and adult dairy cattle (Farm B). Subsequent to the outbreak, FMD virus (FMDV) asymptomatic carriers were identified in both herds by seroreactivity to FMDV non-structural proteins and detection of FMDV genomic RNA in oropharyngeal fluid. Carriers' seroreactivity and FMDV genome detection status were subsequently monitored monthly for 23 months. The mean extinction time of the carrier state was 13.1 ± 0.2 months, with extinction having occurred significantly faster amongst adult dairy cattle at Farm B compared to younger animals at Farm A. The rate of decrease in the proportion of carrier animals was calculated to be 0.07 per month. Seroprevalence against FMDV non-structural proteins decreased over the course of the study period, but was found to increase transiently following repeated vaccinations. These data provide novel insights into viral and host factors associated with the FMDV carrier state under natural conditions. The findings reported herein may be relevant to field veterinarians and governmental regulatory entities engaged in FMD response and control measures.

Isolation and characterization of foot-and-mouth disease virus from Odisha, India

Foot-and-mouth disease (FMD) is a highly contagious and rapidly transmissible disease of cloven footed animals. Emergence of genetically divergent strains of FMD virus (FMDV) is a major concern globally. FMD is endemic in India and three serotypes (O, A and Asia 1) prevail. The study was undertaken to characterize the isolates from the state of Odisha, India both genetically and antigenically. FMDV was detected in 7 of the 17 clinical samples collected from FMD affected/suspected animals, in which serotype O and A were found in three and four samples, respectively. Serotype O field isolates clustered in an unnamed group (designated here as Eastern cluster) circulating mostly in the Eastern region of the country and had 10-12.7% divergence from the Ind2001 lineage circulating predominantly throughout the country. The serotype A isolates sequenced in this study was grouped within VP3⁵⁹-deletion group of genotype 18, precisely in clade 18c, having high genetic homology to the virus circulating in the neighboring states, suggesting interstate movement. Both the serotype O and A isolates showed good antigenic relationship value with the respective vaccine strains currently used in the country.

Spectrum of VP1 region genetic variants in the foot-and-mouth disease virus serotype O populations derived from infected cattle tongue epithelium

RNA virus population exists as a complex distribution of non-identical but closely related sequences known as viral quasispecies. Variant strains are selected from this quasispecies population in response to changing environment. The quasispecies dynamics of a virus existing within an infected host differs from that in a cell culture-adapted population. This study was carried out to explore the genetic variations present in the VP1 coding region of the foot-and-mouth disease (FMD) virus serotype O derived directly from infected cattle tongue epithelium. Molecular clonal populations of two serotype O strains belonging to lineages Ind2001 (IND 30/2011) and PanAsia2 (IND 5/2011) were sequenced at VP1 coding region. For IND 30/2011, 19 clones were sequenced and analysis showed variations at 12 nucleotide positions (nt) resulting in 8 amino acid (aa) replacements. Similarly, for IND 5/2011 virus, 18 clones were sequenced, of which six showed nt variations leading to 3 aa replacements. Most of the variable positions mapped to the surface-exposed loops and some of them were found in the neutralizing antigenic sites (position 81, 149, 169, 186 and 202 of IND 30/2011 and 141 of IND 5/2011), which potentially could be beneficial in rapid adaptive evolution of the virus by giving rise to antigenic variants to overcome neutralizing antibodies. These findings encourage further research into the landscape of the viral quasispecies population in vivo and its implication for viral ecology.

Sequence analysis of capsid coding region of foot-and-mouth disease virus type A vaccine strain during serial passages in BHK-21 adherent and suspension cells

Sequence variability within the capsid coding region of the foot-and-mouth disease virus type A vaccine strain during serial in vitro passage was investigated. Specifically, two methods of virus propagation were utilized, a monolayer and suspension culture of BHK-21 cells. At three positions (VP2(131) E-K in both monolayer and suspension passages, VP3(85) H-R in late monolayer passages and VP3(139) K-E in only suspension passages), all mapped to surface exposed loops, amino acid substitutions were apparently fixed without reversion till the end of the passage regime. Interestingly, VP2(131, 121) and VP3(85) which form part of the heparan sulphate binding pocket, showed a tendency to acquire positively charged amino acids in either monolayer or suspension environment probably to better interact with the negatively charged cell surface glycosaminoglycans. At three identified antigenically critical positions (VP2(79), VP3(139) and VP1(154)), amino acids substitutions even in the absence of immune pressure were noticed. Hence both random drift and adaptive mutations attributable to the strong selective pressure exerted by the proposed cell surface alternate receptors could play a role in modifying the capsid sequence of cell culture propagated FMDV vaccine virus, which in turn may alter the desired potency of the vaccine formulations.

Experimental evidence for competitive growth advantage of genotype VII over VI: implications for foot-and-mouth disease virus serotype A genotype turnover in nature

In India, systematic genotype replacement has been observed for serotype A foot-and-mouth disease virus. After a decade of co-circulation of genotypes VI and VII, genotype VII emerged as the single dominant genotype since 2001. To derive possible explanations for such epochal evolution dynamics, in vitro intergenotype growth competition experiments involving both co- and superinfection regimes were conducted. Coinfection of BHK-21 cells demonstrated abrupt loss in the genotype VI viral load with commensurate increase in the load of genotype VII as measured by the genotype differentiating ELISA, RT-PCR and real-time RT-PCR. The superinfection dynamics was shaped by temporal spacing of infection, where the invading genotype VII took more number of passages than coinfection to eventually overtake the resident genotype VI. It was speculated that such superior replicative fitness of genotype VII could have been a possible factor for the ultimate dominance of genotype VII in nature.

Phylogenetic analysis of 3C protease (3C(pro)) coding region of Foot-and-mouth disease virus type A

Nucleotide sequence analysis of the 3C protease (3C(pro)) region of Foot-and-mouth disease virus type A (FMDV-A) isolates from India has revealed incongruous phylogenetic grouping between 3C(pro) and VP1 region possibly due to the genetic recombination or independent evolution of non-structural and structural protein coding regions. Similar to the VP1 region, the emerging VP3(59)-deletion group maintained its genetic distinctiveness at 3C(pro) region and was found to be diverging with time. Two lineage specific signature aa residues were detected for the deletion group in proof of lineage specific drift or selection events. 3C(pro) region exhibited high degree of conservation as evident from low dN/dS ratio (0.036) and percentage of variable aa positions (20%). A transmembrane domain from aa 27 to 44 could be predicted that possibly anchors 3C to intracellular membranes for better interaction with RNA replication complex. On the basis of sequence conservation, the likelihood that the region aa 121-150 was carrying a vaccine exploitable T-cell epitope was very high.

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

A total of 18 foot-and-mouth disease virus (FMDV) serotype Asia1 field isolates belonging to two different lineages (including the divergent group) as delineated earlier in VP1-based phylogeny were sequenced in the non-structural 3A and 3C protein-coding regions. The phylogenetic trees representing the regions coding for the non-structural proteins were very similar to that of the structural VP1 protein-coding region. Phylogenetic comparison at 3C region revealed clustering of Asia1 viruses with the isolates of serotypes O, A and C in the previously identified clade. Comparison of amino acid sequences identified lineage-specific signature residues in both the non-structural proteins. Overall analysis of the amino acid substitutions revealed that the 3A coding region was more prone to amino acid alterations than 3C region.

A novel genetic lineage differentiating RT-PCR as a useful tool in molecular epidemiology of foot-and-mouth disease in India

Comparison of nucleotide sequences at the VP1 coding region of foot-and-mouth disease serotype Asia1 viruses from India has revealed two genetic lineages with emergence of a genetically divergent group in recent years. In this study a simple, fast, relatively costeffective multi-primer RT-PCR assay to differentiate genetic lineages of type Asia1 viruses was developed. Efforts were made in the design of novel lineage-specific primers and in optimization of the multi-primer assay protocol in conjunction with the use of the serotype specific primer for confirmation of serotype Asia1 virus. This assay promises to be an effective tool in molecular epidemiological investigation of FMD in the country.

Genetic comparison of large fragment of the 5'untranslated region among foot-and-mouth disease viruses with special reference to serotype Asia1

Foot-and-mouth disease (FMD), the most economically important disease of cloven-hoofed animals, is endemic in India. Sequence analysis revealed that phylogenetic grouping of type Asia1 field isolates on the basis of the large fragment of the 5'untranslated region (5'LF-UTR) was quite similar to that based on the sequences of the capsid-coding (VP1) region of the same viruses. The existence of two distinct lineages of type Asia1 suggested by the study on the VP1 region was further supported by the detection of a difference in length and predicted secondary structure of the 5'LF-UTR between the two lineages. Sequence variability between the isolates of the two lineages was also observed within the different domains of the internal ribosome entry site (IRES) around conserved motifs like the GNRA,- RAAA,- and the polypyrimidine tract. Certain group and lineage-specific signature nucleotides pertaining to FMDV type Asia1 in the 5'LF-UTR have been identified. The present study shows that the 5'LF-UTR of FMDV serotype Asia1 field isolates are variable in relation to the length and probable secondary structure of the IRES.

Sequence variability in the structural protein-encoding region of foot-and-mouth disease virus serotype Asia1 field isolates

A total of 30 field isolates of foot-and-mouth disease virus (FMDV) serotype Asia1 belonging to two different lineages and five isolates belonging to a divergent group as delineated earlier in 1D (encodingVP1 protein) gene-based phylogeny were sequenced in the structural protein (P1) coding region. Phylogenetic comparison of these isolates along with some of the published exotic sequences revealed the presence of five different lineages around the world. Similar grouping pattern was observed for the P1 region and 1D gene-based phylogeny, where the Indian isolates were clustered in two genetic lineages. The recently identified divergent group of virus falls into a separate sub-cluster. Similar grouping was also observed in L gene-based phylogeny. Comparison of amino acid sequences identified lineage-specific signature residues in all the structural proteins. Comparison of Asia1 field isolates at the identified key residues of other FMD viruses involved in the formation of the heparan sulfate-binding ligand confirmed many of them to be conserved and the presence of VP3(56) Arg suggested their cell culture adaptation. Although a considerable genetic variation was observed among the isolates of present study, all of them tested in micro-neutralization test were serologically related to the vaccine strain.

Characterization of foot-and-mouth disease serotype Asia1 viruses grown in the presence of polyclonal antisera in serology and nucleotide sequence analysis

Foot-and-mouth disease viruses (FMDV) have a high rate of mutation and spontaneous mutants can be readily. isolated in the laboratory. In this study, plaque purified FMDV Asial vaccine strains (IND 63/72 and IND 491/97) were passaged in-vitro in Baby Hamster Kidney-21 cell monolayers in the presence of sub-neutralizing levels of antiviral polyclonal sera (APS), raised in guinea pigs against the purified and inactivated whole virus particles of IND 63/72, IND 491/97 and IND 13/01. After serial passages under selective immune pressure, the viruses starts growing in the presence of undiluted sera and showed certain characteristics like an increased resistance to neutralization by APS and reduction in plaque counts on titration in plaque assay. Cross-neutralization of these viruses with above-mentioned APS revealed selection of three complete and one partial polyclonal antibody resistant (PAR) viruses based on the 'r' value in micro neutralization test. Alterations were detected at several amino acid residues in the structural protein-coding P1 region. Many of the residues inferred to be positively selected sites in other serotypes of this virus were also prone to substitution under immune selection pressure in Asia1 virus. The present work extends the finding that selection exerted by host antibody also plays a major role in the rapid evolution of FMDV Asia1, as observed in other serotypes.

Complete nucleotide sequence analysis of a vaccine strain and a field isolate of foot-and-mouth disease virus serotype Asia1 with an insertion in VP1 genomic region

Complete nucleotide sequences except the poly (C) tract and poly (A) tail of a vaccine strain (IND 491/97) and an atypical field isolate (IND 321/01) of Foot-and-mouth disease virus (FMDV) serotype Asia1 are described. Amino acid (aa) sequence analysis of the VP1 protein of the field isolate revealed that the latter has 212 instead of 210 or 211 aa found in the so far available sequences of other FMDV isolates of Asia1 serotype. The insertion was localized in the hypervariable region of aa 130-160 of VP1 protein. Nucleotide sequencing of the entire genome was therefore carried out to detect changes in other parts of the genome, if any, besides VP1, which could contribute to its fitness. An 8.16 kb sequence of IND 491/97 and an 8.162 kb sequence of IND 321/01 were compared with each other and also with the known sequence of IND 63/72, another vaccine strain of serotype Asia1. Comparison of the entire polyprotein coding (L to 3D) region of IND 321/01 with those of the two Asia1 vaccine strains (IND 63/72 and IND 491/97) revealed no significant differences. A similar comparison of IND 491/97 with IND 63/72 revealed variability across the entire length of the genome. In addition to the capsid-coding region, sequence variability was also observed in non-structural proteins albeit to different extent. This study shows that in the gene pool of serotype Asia1 at least three groups of isolates/strains are present with respect to the length of VP1 protein.

Mutation in the 1D gene (VP1) of Foot-and-mouth disease virus serotype Asia1 during serial cytolytic infections in cell culture

Changes in the nucleotide sequence of the 1D gene of two vaccine strains (IND 63/72 and IND 491/97) of Foot-and-mouth disease virus (FMDV) serotype Asia1 during serial cytolytic infections in cell culture have been analyzed. Sequence comparisons revealed a majority of transition mutations in IND 491/97. The mutation frequency of the 1D gene of IND 491/97 was about 4.5 to 6.0 fold higher than that of IND 63/72. At the amino acids 40-60 and 140-160 regions the mutation frequency was higher compared to the whole VP1. Both viruses showed a constant change at certain residues of the G-H loop region with an accumulation of amino acid replacements during serial cytolytic passages in cell culture. The critical residues (145 and 153) identified previously using mAbs recognizing trypsin-sensitive epitopes were not substituted in the absence of immune selection but changes were observed at positions 142 and 148. Non-reactivity of IND 63/72 after 50(th) passage level onwards with a panel of mAbs indicated an alteration in the antigenic specificity of the virus. Comparison of amino acid sequences in the entire capsid coding region of the naturally occurring field isolates with that of the 50(th) and 100(th) passage level viruses of IND 63/72 revealed that the residues 56 and 74 of VP2 could be involved in mAb binding. The results suggest that fixation of amino acid replacements occurs in VP1 of Asia1 virus, which could play an important role in antigenic variation by modulating different antigenic epitopes located on the surface of the virus.