Comparative complete genome analysis of Indian type A foot-and-mouth disease virus field isolates.

Comparative complete genome analysis of 17 serotype A Indian field isolates representing different genotypes and sub-lineages is presented in this report. Overall 79% of amino acids were invariant in the coding region. Chunk deletion of nucleotide was observed in S and L fragment of 5'-UTR. More variability which is comparable to that of capsid coding region was found in L and 3A region. Functional motifs and residues critical for virus biology were conserved most. Polyprotein cleavage sites accepted few changes. Many sites were detected to be under positive selection in L, P1, 2C, 3A, 3C, and 3D region and of which some are functionally important and antigenically critical. Genotype/lineage specific signature residues could be identified which implies evolution under different selection pressure. Transmembrane domain could be predicted in 2B, 2C, 3A, and 3C proteins in agreement with their membrane binding properties. Phylogenetic analysis at complete coding region placed the isolates in genotype IV, VI, and VII and two broad clusters comprising VP3(59)-deletion and non-deletion group within genotypes VII. The VP3(59)-deletion group has diversified genetically with time giving rise to three lineages. Incongruence in tree topology observed for different non structural protein coding region and UTRs-based phylogeny indicate suspected recombination.

Multiple origins of foot-and-mouth disease virus serotype Asia 1 outbreaks, 2003-2007.

We investigated the molecular epidemiology of foot-and-mouth disease virus (FMDV) serotype Asia 1, which caused outbreaks of disease in Asia during 2003-2007. Since 2004, the region affected by outbreaks of this serotype has increased from disease-endemic countries in southern Asia (Afghanistan, India, Iran, Nepal, Pakistan) northward to encompass Kyrgyzstan, Tajikistan, Uzbekistan, several regions of the People's Republic of China, Mongolia, Eastern Russia, and North Korea. Phylogenetic analysis of complete virus capsid protein 1 (VP1) gene sequences demonstrated that the FMDV isolates responsible for these outbreaks belonged to 6 groups within the Asia 1 serotype. Some contemporary strains were genetically closely related to isolates collected historically from the region as far back as 25 years ago. Our analyses also indicated that some viruses have spread large distances between countries in Asia within a short time.

Analysis of the leader proteinase (L(pro)) region of type A foot-and-mouth disease virus with due emphasis on phylogeny and evolution of the emerging VP3(59)-deletion lineage from India.

Genotype inclusive grouping of Indian type A isolates as observed in 1D region based phylogeny was distorted at complete L(pro) region, where the VP3(59)-deletion group lineages of genotype VII clustered away from both genotypes VII and VI, confirming its uniqueness and independent evolution of L(pro) and 1D region. Akin to the 1D region, this deletion group is gradually diverging genetically even at L region forming more number of lineages and inter-lineage distance at L region is considerably more than that for 1D region. The deletion group is restricted to India only as none of the exotic sequences clustered within this group. Notably, L protein exhibited variability comparable to external capsid proteins as evident from its high d(N)/d(S) ratio (0.105), number of variable amino acid positions (41%), low Ts/Tv ratio (3.47) and alignment revealed N-terminal region, beta2 sheet and C-terminal extension to be extremely variable. Basic residues at P1, P3 and only leucine at P2 were predicted to provide an optimum autocatalytic cleavage site at L/P1 junction. All of the eight sites identified to be under positive selection revealed aa substitutions of varied physicochemical properties and at two positions lineage specific signatures were observed, which supports the contention that lineages are evolving under differential selection pressure to adapt to the varied ecological environment.

Comparative genomics of serotype Asia 1 foot-and-mouth disease virus isolates from India sampled over the last two decades.

This study deals with a comparative analysis of complete genome sequences of twenty-one serotype Asia 1 foot-and-mouth disease (FMD) field viruses isolated over a period of two decades from India, two vaccine strains and seven exotic sequences. The Indian viruses could be grouped in to three distinct lineages at the entire coding region, evolving independently probably under differential selection pressure as evident from the lineage-specific signatures identified. This comparison revealed 80% of amino acids at the polyprotein region to be invariant. Twenty-one residues in L, 3A and P1 region were identified to be under positive selection of which some are antigenically critical. Analysis at functionally crucial motifs, receptor contact residues, polyprotein cleavage sites and at putative T-cell epitopes expands the knowledge beyond other serotypes. Antigenic site II in betaB-betaC loop of VP2 was highly unstable suggesting its exposure to extreme immune pressure. A single cross-over at the L proteinase region in an isolate from buffalo, also featuring an extensive deletion at the 5' untranslated region (UTR), reflects the role of intraserotypic genetic recombination in natural evolution. The likely biological relevance of deletions/insertions observed at UTRs, VP1 and 3A could not be deduced. Altogether, a substantial amount of data raised on full length genomes of type Asia 1 virus adds value to the FMD virus genomics.

Assessment of suitability of two serotype A candidate vaccine strains for inclusion in FMD vaccine in India.

The recent type A foot and mouth disease virus field isolates recovered in India are shown to be antigenically quite divergent from the in-use vaccine strain (IND 17/82), warranting the selection of a suitable vaccine strain which can cover this diversity in antigenic spectrum. In earlier studies employing neutralization test with anti-146S rabbit sera raised against eight candidate vaccine strains, IND 81/00 and IND 40/00 belonging to genotype VII were found to offer the best antigenic coverage. In order to assess the credibility of IND 81/00 and IND 40/00 as vaccine strains, 17 recent isolates received during 2005-2006 and representative isolates from older genotypes were subjected to two-dimensional micro-neutralization assay using bovine convalescent serum (against IND 81/00 and IND 40/00) and bovine vaccinate serum (against IND 40/00). From the results it is evident that both the isolates IND 81/00 (antigenic relationship 'r-value' >0.40 with 86% of isolates) and IND 40/00 ('r-value' >0.40 with 78% of isolates) show nearly equal antigenic relatedness with the recent field viruses and hence both of these are effective vaccine candidates in present context. Though very limited in its extent, these useful data obtained with antisera raised in homologous host system are logical extension of the on going quest for the appropriate vaccine strain and circumvents species disparities in the immune recognition of epitopes.

Genotype differentiating RT-PCR and sandwich ELISA: handy tools in epidemiological investigation of foot and mouth disease.

Years of molecular epidemiological surveillance has revealed co-circulation of two antigenically divergent genotypes of foot and mouth disease virus serotype A in India. Genotype differentiating RT-PCR and sandwich ELISA were developed as fast, cost-effective and user-friendly alternatives to 1D region based phylogeny for detection and differentiation of genotype VI and VII. The RT-PCR assay targeting 1D region was found to be more sensitive and authentic in distinguishing genotypes than sandwich ELISA. These assays promise to be reliable tools in the epidemiological investigation of foot and mouth disease in the country.

Development and comparison of genome detection assays for the diagnosis of foot-and-mouth disease suspected clinical samples.

Detection of foot-and-mouth disease virus (FMDV) from clinical specimens by conventional sandwich enzyme-linked immunosorbent assay (ELISA) and virus isolation in cell culture is often compromised owing to limited sensitivity and inactivation during transit, respectively. A RT-PCR (oligoprobing) ELISA in both solid and aqueous phase hybridization formats targeting an across serotype conserved site at 3C-3D region was developed and its effectiveness was compared with that of the known targets at the IRES region. A non-isotopic RNA dot hybridization assay with colorimetric detection targeting both the IRES and the 3D region were also validated, which is capable of handling high throughput samples with ease. RT-PCR (oligoprobing) ELISA and dot hybridization assay showed 1000- and 10-fold greater sensitivity than the sandwich ELISA, respectively. Robustness of these diagnostic methods was explored by examining on sandwich ELISA-negative clinical samples. Both the assays developed in the present study were able to detect viral genomes in samples undetectable by conventional ELISA, thereby demonstrating 'proof of sensitivity'. Although the potential of these assays for providing definitive diagnosis in carrier hosts and in species where clinical disease is inapparent remains to be examined, nevertheless these assays can be adapted for comprehensive surveillance of foot-and-mouth disease in India.

Reverse transcription loop-mediated isothermal amplification assays for rapid identification of eastern and western strains of bluetongue virus in India.

Bluetongue virus (BTV) infects all ruminants, including cattle, goats and camelids, causing bluetongue disease (BT) that is often severe in naïve deer and sheep. Reverse-transcription-loop-mediated-isothermal-amplification (RT-LAMP) assays were developed to detect eastern or western topotype of BTV strains circulating in India. Each assay uses four primers recognizing six distinct sequences of BTV genome-segment 1 (Seg-1). The eastern (e)RT-LAMP and western (w)RT-LAMP assay detected BTV RNA in all positive isolates that were tested (n=52, including Indian BTV-1, -2, -3, -5, -9, -10, -16, -21 -23, and -24 strains) with high specificity and efficiency. The analytical sensitivity of the RT-LAMP assays is comparable to real-time RT-PCR, but higher than conventional RT-PCR. The accelerated eRT-LAMP and wRT-LAMP assays generated detectable levels of amplified DNA, down to 0.216 fg of BTV RNA template or 108 fg of BTV RNA template within 60-90min respectively. The assays gave negative results with RNA from foot-and-mouth-disease virus (FMDV), peste des petits ruminants virus (PPRV), or DNA from Capripox viruses and Orf virus (n=10), all of which can cause clinical signs similar to BT. Both RT-LAMP assays did not show any cross-reaction among themselves. The assays are rapid, easy to perform, could be adapted as a 'penside' test making them suitable for 'front-line' diagnosis, helping to identify and contain field outbreaks of BTV.

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.

Recent reports indicated presence of two antigenic and genetic groups (genotypes VI and VII) of foot-and-mouth disease virus (FMDV) type A in India and are divergent from the vaccine strains. In order to choose suitable field isolate as candidate vaccine strain, anti-sera against representative isolates from both the genotypes and two in-use vaccine strains are tested in neutralization assay. Two candidate vaccine strains from both the genotypes are selected with close antigenic match to the field isolates. From the result it is evident that IND 81/00 (genotypes VII), gave a better antigenic coverage (antigenic relationship (r)-value>0.40 with 79% of isolates of 2002--2003) than IND 258/99 (genotype VI; r-value>0.40 with 42% of 2002--2003 isolates). Phylogenetic analysis based on P1 genomic region placed all the recent isolates (2001--2003) into genotype VII, with emergence of a new variant virus (VIIb--VP3(59)del) having amino acid deletion at an antigenically critical residue (VP3(59)), indicating a major evolutionary jump probably due to immune selection. Though very limited in its extent, this data indicates an apparent dominance of genotype VII over genotype VI and underscores the need to continue further molecular epidemiological investigations to substantiate this finding.

Development and evaluation of a multiplex PCR for differentiation of foot-and-mouth disease virus strains native to India.

A multiplex PCR (mPCR) for the differentiation of Indian FMDV serotypes, O, A, Asia 1 and C was developed and evaluated on 142 clinical and 39 cell culture samples. On the latter samples both the tests worked well with 100% efficiency, whereas on clinical samples mPCR had better efficiency than ELISA. The test was found to be specific for FMDV. The detection limit of the assay was varied among the serotypes; it was most sensitive on types A and Asia 1 and least sensitive on type C. The mPCR clearly identified the serotype and in some cases detected dual infections. The test is sensitive and reliable and can be used for serotyping of ELISA negative samples.

Sequence and phylogenetic analysis of the L and VP1 genes of foot-and-mouth disease virus serotype Asia1.

Most of the molecular epidemiological studies of foot-and-mouth disease virus (FMDV) are based on comparison of VP1 gene sequence. In this report, we determine the nucleotide (nt) sequence of the L (603 nt) and VP1 (633 nt) genes of 27 FMDV serotype Asia 1 isolates recovered from different outbreaks in India, and compared with each other and the vaccine strain, IND 63/72, used in the country. Independent phylogenetic analyses on both the aligned gene sequences identified two major lineages (designated A & B) in the Asia 1 isolates. Both L- and VP1-based trees were congruent with respect to the major branching pattern of the isolates. The lineage A is represented by the isolates of 1986-2000 including the vaccine strain IND 63/72, whereas, lineage B appeared to be dominant and responsible for most of the recent outbreaks. A correlation was observed between the clustering of the isolates in the phylogenetic tree and the amino acid changes at many of the positions in VP1 as well as in L protein. The annual rate of evolution in L and VP1 genes was found similar and estimated to be 4.0 x 10(-3) and 3.8 x 10(-3) substitutions per nucleotide, respectively. Our result, largely from the congruence in phylogenetic trees and the rate of evolution in both the genes, suggests the possibility for the use of L gene sequence in phylogenetic comparison of FMDV.

Phylogenetic analysis of Indian serotype Asia1 foot-and-mouth-disease virus isolates revealed emergence and reemergence of different genetic lineages.

Foot-and-mouth-disease virus (FMDV) serotype Asia1 causes significant number of disease outbreaks in India. Indian Asia1 virus isolates were shown to be genetically heterogeneous and of the two lineages (lineage B and lineage C) described in India, lineage C caused majority of the outbreaks. Emergence of a novel divergent lineage (lineage D) within lineage C has been described in 2001. In the present report, the complete VP1 genomic region of 41 FMDV Asia1 field isolates collected between 2003 and 2008 was sequenced. Phylogenetic analysis revealed reemergence of lineage C since 2005 following exclusive dominance of lineage D in the period between 2002 and 2004. At many positions lineage specific signature residues were identified. The antigenic relationship of the field isolates with the currently used vaccine strain IND63/72 was also determined, which reflects antigenic stability of serotype Asia1 in-spite of genetic heterogeneity.

Genetic characterization of type A foot-and-mouth disease virus 3A region in context of the reemergence of VP359-deletion lineage in India.

The 3A region of foot-and-mouth disease virus has been implicated in host range and virulence. Here we analyzed the 3A region of serotype A virus in view of the emergence of a variant group in India with an amino acid deletion at an antigenically critical position of capsid protein, VP3. The 3A region exhibited extreme variability with 38% of the amino acid positions showing substitutions and the C-terminal third (127-151) region was most flexible. Genotype inclusive grouping of type A foot-and-mouth disease virus as observed in 1D region based phylogeny was much less apparent at 3A region possibly due to independent evolution of nonstructural and structural protein coding regions. Akin to the 1D region, the VP3(59)-deletion group maintained its phylogenetic distinctness even at the 3A region and was found to be diverging with time. Twelve lineage specific signature amino acid residues, of which four were identified to be experiencing positive selection, indicates fixation of advantageous mutations in a lineage specific manner. Six positions, all located in the hypervariable C-terminal third, were identified to be under positive selection and were presumed to be imparting the virus certain advantage accounting for its adaptability to wide host spectrum and rapid dissemination. A significant change of Q(44)H was noted only in the older lineage (VIIb) of the deletion group at a position where Q(44)R mutation is associated with guinea pig adaptation. As this site has been detected to be under positive selection, such a lineage specific substitution is thought to have imparted certain temporary advantage to the virus during its possible adaptation in wild or some understudied domestic hosts and must not have seriously compromised fitness upon readaptation in bovines. A conserved hydrophobic transmembrane domain from position 59 to 76 could be predicted which possibly anchors 3A to intracellular membranes for successful interaction with RNA replication complex.

Evaluation of in vitro inhibitory potential of small interfering RNAs directed against various regions of foot-and-mouth disease virus genome.

India is endemic for foot-and-mouth disease and it continues to be a major threat to the livestock industry despite vaccination programmes. In the present study, the ability of specific small interfering (si)RNAs directed against different genomic regions of foot-and-mouth disease virus (FMDV) to inhibit virus replication in BHK-21 cells was examined. For preliminary evaluation of possible siRNA-mediated FMDV inhibition, a cocktail of several unique populations of 12-30bp siRNAs were successfully produced corresponding to three target regions located at structural (VP3-VP1), non-structural (2A-2C), and non-structural-untranslated (3D-3'UTR) region of serotype Asia1. Once the populations of siRNAs generated were found to reduce the virus titre significantly, two highly conserved 21bp siRNA duplexes were designed by analysing all FMDV sequence entries available in public-domain databases. In virus titration assay, more than 99% inhibition of virus yield for all the four serotypes (type Asia1, O, A, and C) could be demonstrated in cells transfected with each of the FMDV-specific siRNAs at 24h post-infection, compared to control cells transfected with scrambled siRNA. This was well supported by reduction in OD values in FMDV-specific sandwich ELISA. Although 100-fold reduction in virus titre with siRNA1 is substantial considering the transfection efficiency and fixed level of input siRNA, siRNA2 emerged to be a better choice as target where more than 300-fold reduction was observed and its inhibitory effect extended up to 48 h post-infection against all the serotypes. Interestingly, in the present study type A virus (IND 17/77) had a single mismatch at position 2 in the siRNA2 target region but it did not abrogate the inhibitory effect.

Evidence of recombination in the capsid-coding region of type A foot-and-mouth disease virus.

Recombination is one of the factors that contribute to genetic diversity in foot-and-mouth disease virus (FMDV). Similarity and bootscan analyses have provided evidence of recombination in the capsid-coding (P1) region of the virus. In the present study, of the 14 subtype A22 field isolates that were distributed in three previously described genotypes (IV, VI and VII) based on the 1D (VP1-encoding) gene sequence (Tosh et al., 2002 ), one isolate (IND 170/88) was found to be a hybrid of genotypes VI and VII in the P1 region. VP1, VP4, the 5' region of VP2 and the 3' region of VP3 of this virus were characteristic of genotype VI, whereas the remaining 3' region of VP2 and the 5' region of VP3 were characteristic of genotype VII. No insertion or deletion was observed in the recombinant virus. Recombination in the P1 region may provide an escape mechanism for the virus.

Emergence of a new strain of type O foot-and-mouth disease virus: its phylogenetic and evolutionary relationship with the PanAsia pandemic strain.

In India, Foot-and-mouth disease virus (FMDV) serotype O has been associated with more than 75% of the outbreaks. Previous studies with this serotype have indicated that the viruses circulating in India belong to a single genotype. Recent (February 2001) FMD epidemics in Europe have focussed global attention on the source of the virus and have been traced to a strain, PanAsia (serotype O), which is present in India since 1990. In this study, to further characterize the isolates belonging to the PanAsian strain, we sequenced the complete VP1-encoding (1 D) gene for 71 FMDV serotype O isolates from India recovered from the field outbreaks during the last 4 decades (1962-2001). All the isolates in the tree were distributed in to three major branches (designated as A, B and C); the branch C is further divided into four groups (I-IV), of which the group IV belongs to the PanAsia strain. Furthermore, we show that the PanAsia strain has been circulating endemically since 1982 (not 1990 as reported earlier) and has been the most dominant outbreak strain in the recent years and distributed at least in 17 states of the country. During the year 2001, another new group (group III) of virus with genetic divergence of 5.4-11.1% at nucleotide level from the PanAsia strain is found to co-circulate endemically, and is slowly replacing it. At amino acid level this strain differed from PanAsia strain at five amino acid positions in the VP1. Although these strains are divergent at nucleotide level, they maintained a good antigenic relationship with one of the vaccine strains (IND R2/75) widely used in the country. Given the ability of the PanAsia virus to persist, spread and to outcompete other strains, the present trend could be of serious concern as the newly emerging virus is replacing it. If this is true, then there is another equally divergent strain as PanAsia that may pose a serious threat to the global dairy and meat industries.