INTRODUCTION

Foot-and-mouth disease (FMD) is a highly contagious viral disease affecting cloven-hoofed livestock, including cattle, goats, sheep, and pigs, and it also infects wild ungulates, particularly species at risk of extinction [1]. Caused by an aphthovirus, FMD spreads rapidly through direct contact with infected animals or via contaminated fomites. Clinically, the disease is characterized by fever, vesicular lesions in the oral cavity and feet, lameness, and reduced milk yield, leading to substantial economic losses for livestock producers. Vaccination programs combined with strict biosecurity measures are therefore essential for effective disease prevention and control. Although FMD has been eradicated in several regions, including Australia, Canada, Greenland, Mexico, and New Zealand, it remains a major transboundary animal disease threat in many parts of Africa, Asia, and South America.

Foot-and-mouth disease virus (FMDV) is a positive-sense single-stranded RNA virus with a genome of approximately 8400 nucleotides. Viral particles are spherical, with diameters ranging from 25 to 30 nm, and consist of an RNA genome enclosed within a protein capsid. The capsid is composed of 60 capsomeres, each containing four structural polypeptides: VP1, VP2, VP3, and VP4. The proteins VP1, VP2, and VP3 are exposed on the viral surface, whereas VP4 is located internally [2]. The surface-exposed proteins (VP1, VP2, and VP3) constitute the major antigenic components of the virus [3]. The VP1 protein contains two key immunogenic regions, namely the G–H loop (amino acid positions 141–160) and the C-terminus (positions 200–213). The G–H loop harbors the arginine–glycine–aspartic acid (RGD) motif, which is essential for viral attachment to host cells through integrin receptors [4]. Owing to its critical role in host cell attachment, immune recognition, and serotype specificity, the VP1 coding region is widely used for genetic characterization of FMDV strains. Sequences of VP1 are also routinely employed in phylogenetic analyses to investigate epidemiological relationships and to trace the origin and spread of virus strains during outbreaks [5].

In 2022, FMD outbreaks continued to affect several Asian countries, including Cambodia, China, peninsular Malaysia, Indonesia, Thailand, and Vietnam, whereas no new outbreaks were reported in Lao PDR, Mongolia, or Myanmar. The importation of cattle and meat products markedly increases the risk of FMD transmission, as the virus can disseminate rapidly through livestock movement between regions. After more than 30 years of disease freedom, Indonesia experienced a major re-emergence of FMD in 2022. The outbreak initially reappeared in Gresik, East Java, and subsequently spread swiftly across multiple provinces, including West Java, Bali, West Kalimantan, Central Kalimantan, East Kalimantan, and West Nusa Tenggara, ultimately affecting almost all provinces on Sumatra.

Phylogenetic analysis reported in a previous study [6] indicated that the virus responsible for the 2022 outbreak originated from Gresik and Banjarnegara in East Java Province. Clinical samples analyzed in the present study were obtained from FMD outbreaks in Indonesia during 2022–2023. Indonesia is now considered endemic for FMD, and the introduction of novel viral lineages into Southeast Asia has been increasingly reported [6].

Despite extensive molecular surveillance of FMDV in Southeast Asia, most studies conducted after the 2022 re-emergence in Indonesia have primarily focused on partial genomic characterization and lineage assignment based on VP1 phylogeny. While these analyses have clarified the geographic origin and spread of circulating strains, they provide limited insight into the functional immunological consequences of emerging genetic variation. In particular, there is a lack of integrated evidence linking sequence divergence in VP1, VP2, and VP3 with changes in T- and B-cell epitope profiles, antigenicity, and innate immune recognition. Moreover, structural features such as ligand-binding pocket variability and host receptor interactions, especially involving TLR7, remain largely unexplored for Indonesian FMDV isolates. This gap hampers the ability to assess potential vaccine mismatch, host-specific immune responsiveness, and the emergence of micro-lineages that may compromise current vaccination strategies. Consequently, a comprehensive, multi-layered analysis combining molecular, immunoinformatic, and structural approaches is needed to better understand the evolving antigenic landscape of FMDV in Indonesia.

The present study aimed to comprehensively characterize FMDV serotype O isolates collected during the 2022–2023 outbreaks in Indonesia using an integrated molecular, immunoinformatic, and structural framework. Specifically, this study sought to (i) determine the genetic diversity and phylogenetic relationships of circulating isolates based on VP1, VP2, and VP3 sequences, (ii) predict and compare VP1-derived T-cell and B-cell epitopes to identify potential antigenic shifts relative to recent reference strains, (iii) evaluate the antigenicity, allergenicity, and toxicity of predicted epitopes using in silico approaches, (iv) analyze variation in VP1 ligand-binding pockets that may influence antigen–antibody interactions, and (v) assess the interaction potential between VP1 and TLR7 through molecular docking. Through this integrated approach, the study aimed to generate functional immunogenomic evidence relevant to vaccine matching and to support adaptive FMD control strategies in Indonesia.

MATERIALS AND METHODS

Ethical approval

All animal procedures were conducted in accordance with Indonesia’s national guidelines for animal welfare and biosafety for transboundary animal diseases. The study protocol, including field sampling, animal handling, and laboratory processing of infectious materials, was reviewed and approved by the Ethical Clearance Committee of the National Research and Innovation Agency (BRIN), Indonesia (approval no. 065/KE.02/SK/10/2022, October 6, 2022). Only naturally infected animals identified during the FMD outbreak were included; no experimental infection, disease induction, or invasive procedures were performed. Oropharyngeal epithelium and lesion scrapings were collected by licensed veterinarians and trained animal health officers following World Organisation for Animal Health (WOAH) guidelines for safe FMD sampling. Animals were handled using minimal-restraint techniques without sedation to avoid pain or distress. Farmers or animal owners provided prior informed consent, and all activities were supervised by the District Livestock Services in West Java and South Sumatra. Sample collection was conducted using appropriate personal protective equipment, and transport complied with national regulations for infectious substances. Laboratory work, including RNA extraction, cDNA synthesis, polymerase chain reaction (PCR) amplification, and sample storage, was performed in a BSL-2 facility with Class II biological safety cabinets. Waste materials were decontaminated by autoclaving or chemical disinfection. The study adhered to Animal Research: Reporting of In Vivo Experiments 2.0, the WOAH Terrestrial Manual (Chapter 3.1.8), and the Three Rs principles.

Study period and location

Samples were obtained from suspected FMD cases based on clinical signs (hypersalivation and vesicular lesions of the oral cavity and feet) in West Java and South Sumatra, Indonesia, during 2022–2023 (Figure 1). Sampling was conducted on smallholder cattle farms raising local breeds (Bos indicus × Bos taurus; Peranakan Ongole, crossbred, and Bali), with herd sizes of approximately 50–100 animals.

Sample collection and transportation

The methodology used in this study is presented in Figure 2. Animals with clinical scores of 2 and 3 (moderate and severe, respectively) presenting vesicles and ulceration of the oral cavity and interdigital spaces of the feet, hypersalivation, and lameness were selected. A total of 28 clinical samples (oral and foot epithelium) were collected and preserved in RNA Shield (Zymo Research, Irvine, CA, USA) (catalog no. R1100-50). The samples were transported in insulated containers with ice packs at approximately 4°C, with temperature validation using temperature logs, over 48–72 h. RNA Shield (Zymo Research) was used to protect nucleic acids at ambient temperature for up to 7 days. All samples were stored at −20°C in the laboratory for further analysis.

RNA extraction and cDNA synthesis

Total RNA extraction was performed according to the manufacturer’s instructions using the Quick-RNA Miniprep Kit (Zymo Research) (catalog no. R1054), with a total RNA elution volume of 30 µL. RNA quality and concentration were evaluated using a NanoDrop spectrophotometer, and RNA with an A260/280 ratio of 1.8–2.0 was selected for further processing. cDNA synthesis was performed using the RevertAid cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA) (catalog no. K1622) to generate complementary DNA. Negative controls (RNase/DNase-free water) were included to exclude contamination. The reaction setup involved incubation at 65°C for 5 min, followed by cDNA synthesis at 25°C for 5 min, 42°C for 60 min, and enzyme inactivation at 70°C for 5 min.

PCR amplification and capsid gene analysis

The cDNA was amplified by conventional PCR using DreamTaq PCR Master Mix (Thermo Fisher Scientific) (catalog no. K1081) with 3D primers for FMD diagnosis [7]. PCR conditions included initial denaturation at 95°C for 5 min, followed by 35 cycles of denaturation at 95°C for 10 s, annealing at 60°C for 20 s, and extension at 72°C for 30 s, with a final extension at 72°C for 5 min. The capsid genes in positive samples were analyzed using primers designed based on the FMDV serotype O reference sequence (MN095357.1) (Table 1). Gene-specific annealing temperatures were applied for VP1, VP2, and VP3. PCR amplification was performed in triplicate, and amplicons were pooled. RNase/DNase-free water was used as a negative control. Amplicons were confirmed on 1.5% agarose gels and purified using a PCR Clean-Up Kit (Geneaid, New Taipei City, Taiwan) (catalog no. DFC100). Purified products were sequenced using Sanger sequencing (BigDye Terminator v3.1, 1st BASE), and sequence identity was confirmed using BLAST against the NCBI GenBank database.

Sequence quality assessment and data submission

Sequence data were analyzed using MEGA-X v10.2.6. Low-quality regions were trimmed based on chromatogram inspection, Phred quality scores VP1 sequences were deposited in GenBank on October 2, 2025 (accession nos. PX406137–PX406155).

Phylogenetic analysis of VP1

Aligned VP1 nucleotide and amino acid sequences were analyzed using MEGA-X v10.2.6 and compared with reference FMDV isolates [8]. Phylogenetic trees were constructed using the maximum likelihood method under the generalized time-reversible model with the Kimura 2-parameter substitution model. Bootstrap analysis was conducted with 1000 replicates, and gaps were treated using partial deletion.

Epitope prediction for MHC class I and B cells

Translated VP1 amino acid sequences were analyzed using IEDB tools to predict MHC class I–binding epitopes employing NetMHCpan 4.1 EL [9]. BoLA alleles were prioritized. Linear B-cell epitope prediction was performed using IEDB B-cell tools (threshold 0.5) [10]. Experimentally validated high- and low-immunogenic peptides were used as positive and negative controls, respectively.

Antigenicity, allergenicity, and toxicity assessment

Antigenicity was evaluated using VaxiJen v2.0 [11]. Allergenicity was assessed using AllerTOP v2.1 and AllergenFP v1.1 [12, 13]. Toxicity prediction was conducted using ToxinPred [14]. Predicted epitopes were benchmarked against experimentally validated epitopes.

Prediction of VP1 ligand-binding sites

Ligand-binding site prediction and 3D visualization were performed using PrankWeb [15]. VP1 structural models were generated using Swiss-Model and assessed using standard quality parameters. Binding pockets were ranked based on calibrated probability scores.

Molecular docking between VP1 and TLR7

Molecular docking between VP1 and TLR7 was performed using the HDOCK server [16]. Docking scores and confidence scores were used to evaluate binding probability, and docking complexes were visualized using Swiss-PdbViewer v4.1.0.

RESULTS

Detection and sequencing of FMDV serotype O capsid genes

Local FMDV serotype O strains were obtained from naturally infected cattle in South Sumatra and West Java during the 2022–2023 outbreak period (Table 2). PCR amplification was performed using serotype O–specific primers targeting the capsid gene coding regions VP1, VP2, and VP3. Sequencing quality was evaluated using Phred quality scores, and only reads meeting a Q20 threshold were included for downstream analysis. Comparative sequence analysis demonstrated that the VP1, VP2, and VP3 genes showed reduced nucleotide identity relative to the most recent Indonesian reference strain, FMDV serotype O isolate ISA/1/2022 (Table 3). For VP1, samples Depok 2, Palembang 8, Sukabumi 1, and Sukabumi 2 exhibited the lowest sequence identity (98%). A similar pattern was observed for VP2, with Sukabumi 1 and Sukabumi 2 also showing 98% identity. In contrast, VP3 identity decreased to 98% in most samples, except for Depok 1 and Bogor 1, which retained higher similarity (99%).

+ detected, − undetected

FMDV = Foot-and-mouth disease virus, ISA = Indonesian isolate.

Phylogenetic analysis of VP1

Phylogenetic reconstruction based on nucleotide substitution models classified the recent FMDV serotype O isolates into three principal clusters. Cluster I comprised Depok 2, Sukabumi 1, Sukabumi 2, and reference FMDV serotype O strains from the database (O/ME-SA/Ind-2001d, Myanmar, India, Oman, Kuwait, Iran, and Bhutan). Cluster II included Bandung 9 and Palembang 1, 2, 3, 4, 6, 7, and 8. Cluster III consisted of Bogor 1, Bandung 1, 3, 4, 5, 8, and Depok 1, which showed close genetic relatedness to Indonesian reference isolates ISA/Gresik/2022 and ISA/Banjarnegara/2022. In contrast, phylogenetic analysis based on amino acid substitutions resolved only two main clusters: Cluster I included Depok 2, Sukabumi 1, Sukabumi 2, and related database isolates, whereas Cluster II comprised the remaining samples together with ISA/Gresik/2022 and ISA/Banjarnegara/2022 (Figure 3).

Comparative analysis of VP1 sequences identified multiple nucleotide variations relative to Indonesian and other serotype O reference sequences, particularly at positions 10, 37, 84, 102, 103, 171, 219, 225, 386, 408, 451, 459, 462, 465, 515, 526, 582, 589, 615, and 627 (Table 4). These substitutions resulted in amino acid changes at positions 34, 113, 129, 172, and 197 (Table 5).

ISA = Indonesian isolate.

ISA = Indonesian isolate

Epitope prediction and immunological profiling of VP1

Based on MHC class I–binding predictions, the VP1 protein yielded eight peptides capable of binding T-cell epitopes associated with BoLA-1, BoLA-2, BoLA-3, BoLA-4, BoLA-5, BoLA-6, BoLA-T5, and BoLA-T7 alleles. Compared with Ind/2022, ISA/Banjarnegara/Gresik, and other FMDV serotype O isolates, samples from Bandung 1, 3, 4, 5, 8, 9, and Palembang 8 exhibited distinct T-cell epitope binding profiles. Samples Depok 2 and Sukabumi 1 and 2 differed from Ind/2022 and ISA/Banjarnegara/Gresik only at amino acid positions 4 (alanine) and 13 (tyrosine), localized within residues 2–15 of VP1 corresponding to the BoLA-T7 allele. No epitope differences were observed for other alleles. All predicted peptides were antigenic and non-toxic, although several were classified as probable allergens (Table 6).

B-cell epitope prediction indicated that Sukabumi 1, Sukabumi 2, and Depok 2 differed from ISA/2022 and ISA/Banjarnegara/Gresik by an alanine-to-threonine substitution at position 13 within the 6–30 epitope region. These profiles were otherwise consistent with other FMDV serotype O database isolates. Notably, Palembang 8 exhibited a glutamic acid–to–lysine substitution at position 196 within the 192–205 epitope region, distinguishing it from Ind/2022, ISA/Banjarnegara/Gresik, and other reference isolates. Peptides spanning positions 93–107 were predicted to be non-antigenic across all samples, whereas other regions demonstrated antigenic potential. Several peptides showed predicted allergenicity, while all remained non-toxic (Table 7).

BoLA = Bovine leukocyte antigen.

Prediction of VP1 protein–ligand binding sites

Protein–ligand binding site analysis of VP1 revealed notable variation among isolates (Table 8). Palembang 8 exhibited the highest number of predicted binding pockets, enriched in threonine, tyrosine, glutamine, and asparagine residues capable of forming hydrogen bonds. Differences in pocket composition between ISA/Gresik/2022, ISA/Banjarnegara/2022, and the field samples suggested altered ligand-binding specificity. Structural visualization showed five high-probability binding pockets in Palembang 8, whereas other isolates displayed fewer than five pockets (Figure 4). Binding pockets ranked 1 or 2 were prioritized as the most probable functional sites.

Molecular docking analysis between VP1 and TLR7

Molecular docking analysis demonstrated variable interaction strength between VP1 and TLR7. Palembang 8 exhibited the lowest docking score and highest confidence score, indicating the strongest predicted interaction with TLR7 (Table 9). A more negative docking score reflected thermodynamically favorable binding, while confidence scores >0.7 indicated a high-probability of interaction; scores between 0.5 and 0.7 suggested possible binding, and values <0.5 indicated unlikely interactions. Palembang 8 formed a dominant docking cluster comprising 65 structures, indicating a stable and reliable binding mode. The associated Z-score (−1.2) confirmed significantly higher stability compared with other clusters. Energy decomposition analysis revealed strong electrostatic and van der Waals contributions, consistent with dense hydrogen-bonding networks and optimal shape complementarity. Residue-level interaction analysis further showed that Palembang 8 had the highest proportion of hydrogen bonds with TLR7, supporting its enhanced binding specificity (Table 10). Visualization of the VP1–TLR7 complex confirmed that key binding epitopes were surface-exposed, facilitating effective receptor engagement (Figure 5).

RMSD = Root mean square deviation.

A = O/ME-SA/Ind, B = MYA/13/2017 (Myanmar), C = PD287/2015 (India), D = PD118/2017 (India), E = PD73/2017 (India), F = IRN/72/2009 (Iran), BHU/3/2009 (Bhutan), G = O/KUW/3/97 (Kuwait), H = INDIA 2013, I = ISA/1/2022, ISA/Gresik/2022, ISA/Banjarnegara/2022, J = Bandung 1, 3, 4, 5, 9, Depok 1, Bogor 1, Palembang 1, 2, 3, 4, 6, 7, 9, K = Depok 2, Sukabumi 1, Sukabumi 2, L = Palembang 8, M = Bandung 8, TLR7 = Toll-like receptor 7, VP1 = Viral protein 1, aa = Amino acid.

DISCUSSION

Persistence of FMD despite vaccination in Indonesia

FMD remains the most damaging viral disease affecting cattle production in Indonesia. Although a national FMD control program has been implemented using an inactivated vaccine containing serotype O strains O/Mya98/XJ/2010, O1 Campos, O1 Manisa, and IND/O/R2/75 (Decree of the Minister of Agriculture of the Republic of Indonesia No. 738/KPTS/PK.300/M/10/2022), FMD continues to circulate within the country. The vaccination campaign has reduced the frequency of clinical disease; however, unregulated animal movement facilitates the ongoing circulation of FMDV strains. In addition, Indonesia’s geographic position increases the risk of incursion of circulating strains from other Asian countries. Consequently, serotype O remains endemic and continues to circulate in Indonesia [6].

Molecular diversity of capsid genes during the 2022–2023 outbreaks

This study presents the first molecular analysis of VP1, VP2, and VP3 sequences from the 2022–2023 outbreaks in South Sumatra and West Java. The nucleotide identity of VP1 decreased to 98.26–99.05% in several samples, while VP3 identity declined to 98.48% in multiple isolates. These findings suggest ongoing antigenic drift that may influence vaccine matching. Phylogenetic analysis based on VP1 nucleotide and amino acid sequences demonstrated that the collected samples formed clusters distinct from recent Indonesian reference isolates ISA/Banjarnegara/2022 and ISA/Gresik/2022. Specifically, three clusters were identified using nucleotide variation and two clusters using amino acid sequences. This pattern contrasts with previous reports in which 2022 Indonesian isolates clustered tightly within Ind-2001d or Mya-98 lineages, highlighting novel intra-country antigenic diversification within a single outbreak period. Multiple nucleotide and amino acid substitutions were detected relative to ISA/Banjarnegara/2022 and ISA/Gresik/2022, including changes at amino acid positions 4, 13, 129, 172, and 197 (Table 6). These substitutions are located near known immunogenic regions, including the G–H loop, B–C loop, and N-terminus of VP1, and may influence receptor binding and antigen recognition.

Evolutionary mechanisms shaping FMDV diversity

The observed genetic heterogeneity of serotype O FMDV across different regions is likely driven by unrestricted livestock movement. Endemic co-circulation of multiple genotypes may facilitate recombination when closely related strains co-infect a host [17]. As an RNA virus lacking proofreading mechanisms, FMDV evolves rapidly, with mutation rates ranging from 10^-3 to 10^-5 substitutions per nucleotide site per year [18]. Each viral population exists as a quasispecies comprising genetically and antigenically diverse variants [19]. Because the principal immunogenic region of FMDV is located within the 1D region, VP1 remains the most informative gene for investigating viral variation and evolutionary relationships among isolates [20].

Implications for vaccine effectiveness and disease control

FMDV evolution predominantly occurs through point mutations, and the extent of genetic change influences viral adaptability, transmissibility, and pathogenicity [21, 22]. The emergence of new strains may compromise disease control in endemic regions when vaccine strains fail to provide adequate cross-protection [23]. Moreover, current vaccines do not prevent the establishment of a carrier state, as infected cattle and buffaloes may excrete the virus for extended periods regardless of vaccination status. Although mortality is uncommon, FMD causes severe disease, and vaccination primarily reduces clinical severity rather than preventing infection or halting transmission.

Immunoinformatic insights into VP1 epitope variation

Epitope prediction analysis of VP1 revealed variation in peptide binding to host T- and B-cell receptors among the collected samples compared with database sequences. All predicted peptides were antigenic and non-toxic, although some showed potential allergenicity based on in silico analysis. The BLV envelope protein (UniProt: P03380) served as a positive control and was predicted to be antigenic, non-allergenic, and non-toxic, supporting the robustness of the analytical pipeline. Certain isolates exhibited unique T-cell binding shifts, such as BoLA-T7–associated changes linked to A4T and A13T substitutions. Given the genetic diversity of Indonesian cattle populations, allele-specific immune responses may influence vaccine performance at regional levels. These findings indicate host-allele-dependent immunogenic divergence, which has not been previously reported for Indonesian FMD outbreaks. Key B-cell antigenic sites of FMDV are located within surface-exposed structural motifs of the capsid, including the N-terminus, B–C, E–F, F–G, G–H loops, and the C-terminus [24]. The conserved RGD motif within the G–H loop of VP1 is essential for integrin-mediated viral entry [25], although linear peptides spanning this region provide limited protection in natural hosts [26]. While neutralizing antibodies are central to protective immunity, T-cell responses are also required for robust immune protection [27], mediated through peptide presentation by MHC molecules to T-cell receptors [28].

Structural and functional relevance of VP1 binding sites

Previous epitope studies of serotype O have primarily focused on neutralizing domains of VP1 at amino acid positions 133, 143–144, 147–149, 151–153, 161, and 164–167 [29], as well as positions 50 and 52 [30] and 199–202 [31]. Computational analyses have predicted multiple B- and T-cell epitopes for serotypes A and O [32], while secondary structure and antigenicity profiling of VP1 has further identified dominant B-cell epitopes [33, 34]. In the present study, ligand-binding pocket analysis revealed notable differences in binding specificity among Indonesian isolates. Palembang 8 exhibited the highest number of binding pockets, enriched with threonine, tyrosine, glutamine, and asparagine residues, distinguishing it from ISA/2022 and other reference isolates. These findings represent the first structural binding-pocket mapping of Indonesian FMDV serotype O isolates.

Molecular docking of VP1 with TLR7 and study limitations

Protein–protein docking analysis demonstrated that Palembang 8 showed the lowest docking and confidence scores for interaction between VP1 and TLR7, indicating a higher probability of binding [39]. This isolate exhibited extensive hydrogen-bonding, strong electrostatic and van der Waals interactions, and high structural complementarity, suggesting enhanced engagement with innate immune receptors. To our knowledge, this is the first report describing the VP1–TLR7 binding potential of FMDV isolates, providing new insights into innate immune recognition mechanisms. Nevertheless, this study is limited by its reliance on in silico predictions without experimental neutralization assays or empirical antigenicity testing. Future work should incorporate wet-lab validation, refined structural modeling, and epitope conservancy monitoring to strengthen these findings.

CONCLUSION

This study provides comprehensive molecular, immunoinformatic, and structural evidence that FMDV serotype O circulating in Indonesia during the 2022–2023 outbreaks is undergoing active genetic and antigenic diversification. Sequence analysis demonstrated reduced nucleotide identity in VP1 (98.26–99.05%) and VP3 (down to 98.48%) compared with recent Indonesian reference isolates. Phylogenetic reconstruction revealed three nucleotide-based clusters and two amino acid–based clusters, indicating intra-country diversification distinct from previously reported tight clustering within Ind-2001d or Mya-98 lineages. Multiple novel amino acid substitutions in VP1 (positions 4, 13, 129, 172, and 197) were identified near key immunogenic regions, including the G–H loop, B–C loop, and N-terminus. Immunoinformatic analyses showed isolate-specific variation in T- and B-cell epitope binding, including BoLA-T7–associated shifts, while structural analyses identified isolate-dependent variation in VP1 ligand-binding pockets and enhanced VP1–TLR7 interaction potential, particularly in Palembang 8.

These findings have direct implications for FMD control in Indonesia. The observed antigenic drift in VP1 suggests that current vaccine strains may gradually lose matching efficiency against emerging local variants. Host-allele-dependent epitope variation further indicates that vaccine responsiveness may differ among cattle populations. The identification of isolates with enhanced TLR7 interaction potential highlights the relevance of innate immune engagement in shaping host responses and may inform future vaccine adjuvant or immunogen design. Continuous molecular surveillance integrating phylogenetics and immunoinformatics is therefore essential to support adaptive vaccine updating and region-specific control strategies.

The major strength of this work lies in its integrated approach, combining molecular sequencing, phylogenetic analysis, epitope prediction, antigenicity/allergenicity/toxicity screening, structural pocket mapping, and protein–protein docking. This multi-layered framework enabled functional interpretation of genetic variation beyond descriptive phylogeny and represents the first report linking VP1 sequence diversity with predicted TLR7 binding potential in Indonesian FMDV isolates.

This study is limited by its reliance on in silico analyses without experimental validation. Neutralization assays, empirical antigenicity testing, and functional immune assays were not performed. In addition, the number of isolates analyzed, although representative, may not capture the full national diversity of circulating FMDV serotype O.

Future studies should include wet-lab validation of predicted epitopes, vaccine cross-neutralization assays, and longitudinal surveillance to monitor epitope conservancy. Structural refinement using higher-resolution models and experimental confirmation of VP1–TLR7 interactions will further strengthen mechanistic understanding. Expanding analyses to additional provinces and integrating host genetic data may also clarify host–virus co-evolution dynamics.

Overall, this study demonstrates that FMDV serotype O in Indonesia is evolving through measurable genetic, antigenic, and structural changes with potential consequences for vaccine effectiveness and immune recognition. The findings underscore the need for continuous, integrative surveillance and evidence-based vaccine adaptation to sustain effective FMD control in endemic settings.

DATA AVAILABILITY

The data supporting this study are included in the manuscript. All full sequence VP1 data were deposited in GeneBank with accession numbers PX406137-PX406155.

AUTHORS’ CONTRIBUTIONS

RIA1, ATW, RIA2, NLPIM, and TPP: Drafted the manuscript. RIA1, SS, TPP, and HR: Collected the sample and performed data analysis. RIA1: Extracted RNA, primer design, and gene amplification. RIA1 and TP: Performed epitope and docking analysis. ATW, RIA2, NLPIM, and TPP: Methodology and software validation. All authors have read and approved the final version of the manuscript.

COMPETING INTERESTS

The authors declare that they have no competing interests.

PUBLISHER’S NOTE

Veterinary World remains neutral with regard to jurisdictional claims in the published institutional affiliations.

ACKNOWLEDGMENTS

This research was supported by RIIM Kompetisi Funding from The Indonesia Endowment Fund for Education Agency (LPDP) No: B-1744/II.7.5/FR/11/2022, B-14145/III.11/KS.00.00/11/2022 and Beasiswa Pendidikan Indonesia (BPI) LPDP

REFERENCES

1. Li P, Huang S, Zha J, Sun P, Li D, Bao H. Evaluation of immunogenicity and cross-reactive responses of vaccines prepared from two chimeric serotype O foot-and-mouth disease viruses in pigs and cattle. Vet. Res 2022;53(56):1-12. [Google Scholar] | [Crossref]
2. Mason P. W, Grubman M. J, Baxt B. Molecular basis of pathogenesis of FMDV. Virus. Res 2003;91:9-32. [Google Scholar] | [Crossref]
3. Liu W, Yang B, Wang M, Liang W, Wang H, Yang D. Identification of a conserved conformational epitope in the VP2 protein of foot-and-mouth disease virus. Arch. Virol 2017;162(7):1877-1885. [Google Scholar] | [Crossref]
4. Zhu Z, Yang F, Cao W, Liu H, Zhang K, Tian H. The pseudoknot region of the 5′untranslated region is a determinant of viral tropism and virulence of foot-and-mouth disease virus. J. Virol 2019;93(8):e02039-18. [Google Scholar] | [Crossref]
5. Abeyratne S. A. E, Amarasekera S. S. C, Ranaweera L. T, Salpadoru T. B, Thilakarathne S. M. N. K, Knowles N. J. Phylogenetic analysis of the VP1 genomic region in foot-and-mouth disease virus serotype O isolates in Sri Lanka reveals the existence of “Srl-97”, a newly named endemic lineage. PLoS One 2018;13(4):e0196491. [Google Scholar] | [Crossref]
6. Zainuddin N, Susila E. B, Wibawa H, Daulay R. S. D, Wijayanti P. E, Fitriani D. Genome sequence of a foot-and-mouth disease virus detected in Indonesia in 2022. Microbiol. Resour. Announc 2023;12(2):e0108122. [Google Scholar] | [Crossref]
7. Manual of Diagnostic Tests and Vaccines for Terrestrial Animals. Paris: WOAH; 2025. [Google Scholar]
8. Zhang X, Ma W, Liu B, Shen C, Yang F, Yang Y. Phylogenetic analyses and antigenic characterization of foot-and-mouth disease virus PanAsia lineage circulating in China between 1999 and 2023. Virol. Sin 2024;39(5):747-754. [Google Scholar] | [Crossref]
9. Reynisson B, Alvarez B, Paul S, Peters B, Nielsen M. NetMHCpan-4.1 and NetMHCIIpan-4.0:Improved predictions of MHC antigen presentation by concurrent motif deconvolution and integration of MS MHC eluted ligand data. Nucleic. Acids. Res 2020;48(1):449-454. [Google Scholar] | [Crossref]
10. Jespersen M. C, Peters B, Nielsen M, Marcatili P. BepiPred-2.0:Improving sequence-based B-cell epitope prediction using conformational epitopes. Nucleic. Acids. Res 2017;45(1):24-29. [Google Scholar] | [Crossref]
11. Doytchinova I. A, Flower D. R. VaxiJen:A server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinform 2007;8((4)). [Google Scholar] | [Crossref]
12. Dimitrov I, Bangov I, Flower D. R. Doytchinova I. AllerTOP v.2 –A server for in silico prediction of allergens. J. Mol. Model 2014;20(6):2278. [Google Scholar] | [Crossref]
13. Dimitrov I, Naneva L, Doytchinova I, Bangov I. AllergenFP:Allergenicity prediction by descriptor fingerprints. Bioinformatics 2014;30(6):846-851. [Google Scholar] | [Crossref]
14. Gupta S, Kapoor P, Chaudhary K, Gautam A, Kumar R, Raghava G. P. In silico approach for predicting toxicity of peptides and proteins. PLoS One 2013;8(9):e73957. [Google Scholar] | [Crossref]
15. Jendele L, Krivak R, Skoda P, Novotny M, Hoksza D. PrankWeb:A web server for ligand-binding site prediction and visualization. Nucleic. Acids. Res 2019;47(2):345-349. [Google Scholar] | [Crossref]
16. Li H, Huang E, Zhang Y, Huang S. Y, Xiao Y. HDOCK update for modeling protein–RNA/DNA complex structures. Protein. Sci 2022;31(11):e4441. [Google Scholar] | [Crossref]
17. Subramaniam S, Mohapatra J. K, Sahoo N. R, Sahoo A. P, Dahiya S. S, Rout M. Foot-and-mouth disease status in India during the second decade of the twenty-first century (2011–2020). Vet. Res. Commun 2022;46:1011-1022. [Google Scholar] | [Crossref]
18. Domingo E, García-Crespo C, Lobo-Vega R, Perales C. Mutation rates, mutation frequencies, and proofreading-repair activities in RNA virus genetics. Viruses 2021;13(9):1882. [Google Scholar] | [Crossref]
19. Fish I, Stenfeldt C, Palinski R. M, Pauszek S. J, Arzt J. Into the deep (sequence) of the foot-and-mouth disease virus gene pool:Bottlenecks and adaptation during infection in naïve and vaccinated cattle. Pathogens 2020;9(3):208. [Google Scholar] | [Crossref]
20. Hassan A. M, Zaher M. R, Hassanien R. T, Abd-El-Moniem M. I, Habashi A. R, Ibraheem E. M. Molecular detection, phylogenetic analysis and genetic diversity of recently isolated foot-and-mouth disease virus serotype A African topotype, Genotype IV. Virol. J 2022;19(1):1. [Google Scholar] | [Crossref]
21. Arzt J, Fish I, Pauszek S. J, Johnson S. L, Chain P. S, Rai D. K. The evolution of a super-swarm of foot-and-mouth disease virus in cattle. PLoS One 2019;14(4):e0210847. [Google Scholar] | [Crossref]
22. Dahiya S. S, Subramaniam S, Mohapatra J. K, Rout M, Biswal J. K, Giri P. Foot-and-mouth disease virus serotype O exhibits phenomenal genetic lineage diversity in India during 2018–2022. Viruses 2023;15(7):1529. [Google Scholar] | [Crossref]
23. Paton D. J, Sumption K. J, Charleston B. Options for control of foot-and-mouth disease:Knowledge, capability and policy. Philos. Trans. R Soc Lond. B Biol. Sci 2009;364(1530):2657-2667. [Google Scholar] | [Crossref]
24. Woldemariyam F, Paeshuyse J. Viral protein 1 (VP1) sequence-based genetic diversity of SAT 2 FMDV circulating in Ethiopia from 1990 to 2015. Vet. Med 2023;14:91-101. [Google Scholar] | [Crossref]
25. Elrashedy A, Nayel M, Salama A, Zaghawa A, El-Shabasy R. M, Hasan M. E. Foot-and-mouth disease:Genomic and proteomic structure, antigenic sites, serotype relationships, immune evasion, recent vaccine development strategies, and future perspectives. Vet. Res 2025;56(78):1-29. [Google Scholar] | [Crossref]
26. Hussein H. A. M, Walker L. R, Abdel-Raouf U. M, Desouky S. A, Montasser A. K. M, Akula S. M. Beyond RGD:Virus interactions with integrins. Arch. Virol 2015;160(11):2669-2681. [Google Scholar] | [Crossref]
27. Medina G. N, Diaz San Segundo F. Virulence and immune evasion strategies of FMDV:Implications for vaccine design. Vaccines 2024;12(9):1071. [Google Scholar] | [Crossref]
28. Babiuk S, Horseman B, Zhang C, Bickis M, Kusalik A, Schook L. B. BoLA class I allele diversity and polymorphism in a herd of cattle. Immunogenetics 2007;59(2):167-176. [Google Scholar] | [Crossref]
29. Zhao Q, Li F, Huang S, Xing X, Sun Y, Li P. Remapping the spatial distribution of neutralizing sites and their immunodominance on the capsid of different topotypes of FMDV serotype O by site-directed competitive ELISA for detection of neutralizing antibodies. Microbiol. Spectr 2025;13(16):e03344-24. [Google Scholar] | [Crossref]
30. He Y, Li K, Cao Y, Sun Z, Li P, Bao H. Structures of foot-and-mouth disease virus with neutralizing antibodies derived from recovered natural host reveal a mechanism for cross-serotype neutralization. PLoS. Pathog 2021;17(4):1009507. [Google Scholar] | [Crossref]
31. Dong H, Liu P, Bai M, Wang K, Feng R, Zhu D. Structural and molecular basis for foot-and-mouth disease virus neutralization by two potent protective antibodies. Protein. Cell 2022;13:446-453. [Google Scholar] | [Crossref]
32. Momtaz S, Rahman A, Sultana M, Hossain M. A. Evolutionary analysis and prediction of peptide vaccine candidates for foot-and-mouth disease virus types A and O in Bangladesh. Evol. Bioinform. Online 2014;10:187-196. [Google Scholar] | [Crossref]
33. Yun T, Liu G. Q, Leng Q. W, Guo H. C, Liu Z. X, Zhang J. N, Xie Q. G. Prediction of secondary structure and B-cell epitope for capsid protein of FMDV. Chin. J. Vet. Sci. Technol 2004;34:23-29. [Google Scholar] | [Crossref]
34. Aggarwal N, Zhang Z, Cox S, Statham R, Alexandersen S, Kitching R. P, Barnett P. V. Experimental studies with foot-and-mouth disease virus, strain O, responsible for the 2001 epidemic in the United Kingdom. Vaccine 2002;20((19–20)):2508-2515. [Google Scholar] | [Crossref]
35. Zhang Y, Liu T, Wang J, Zou B, Li L, Yao L. Cellinker:A platform of ligand–receptor interactions for intercellular communication analysis. Bioinformatics 2021;37:1-8. [Google Scholar] | [Crossref]
36. Wang H. Prediction of protein–ligand binding affinity via deep learning models. Brief. Bioinform 2024;25(2):1-17. [Google Scholar] | [Crossref]
37. Macalino S. J, Gosu V, Hong S, Choi S. Role of computer-aided drug design in modern drug discovery. Arch. Pharm. Res 2015;38(9):1686-1701. [Google Scholar] | [Crossref]
38. Anderson A. C. The process of structure-based drug design. Chem Biol 2003;10(9):787-797. [Google Scholar] | [Crossref]
39. Dominguez C, Boelens R, Bonvin A. M. HADDOCK:A protein–protein docking approach based on biochemical or biophysical information. J. Am. Chem. Soc 2003;125(7):1731-1737. [Google Scholar] | [Crossref]