Molecular epidemiology of foot-and-mouth disease virus serotypes A and O with emphasis on Korean isolates: temporal and spatial dynamics

Phylogenetic and evolutionary analysis of foot-and-mouth disease virus A/ASIA/Sea-97 lineage

Foot-and-mouth disease virus (FMDV) A/ASIA/Sea-97 is a predominant lineage in Southeast Asia and East Asia. However, Sea-97 lineage has not been well studied since its first outbreak in Thailand in 1997. Thus, we conducted phylogenetic and evolutionary analysis of Sea-97 using 224 VP1 sequences of FMDV A/ASIA during 1960 and 2018. Phylogenetic analysis revealed that Sea-97 lineage can be classified into five groups (G1–G5). After the emergence of G2 from G1, the genetic diversity of Sea-97 increased sharply, causing divergence into G3, G4 and G5. During this evolutionary process, Sea-97 lineage, which was initially found only in some countries in Southeast Asia, gradually spread to East Asia. The evolution rate of this lineage was estimated to be 1.2 × 10^–2 substitutions/site/year and there were many differences in amino acid residues compared to vaccine strain. Substitutions at antigenically important sites may affect the efficacy of the vaccine, suggesting the need for appropriate vaccine strains. Our results could provide meaningful information to understand comprehensive characteristic of Sea-97 lineage.

Picornaviruses: A View from 3A

Picornaviruses are comprised of a positive-sense RNA genome surrounded by a protein shell (or capsid). They are ubiquitous in vertebrates and cause a wide range of important human and animal diseases. The genome encodes a single large polyprotein that is processed to structural (capsid) and non-structural proteins. The non-structural proteins have key functions within the viral replication complex. Some, such as 3D^pol (the RNA dependent RNA polymerase) have conserved functions and participate directly in replicating the viral genome, whereas others, such as 3A, have accessory roles. The 3A proteins are highly divergent across the Picornaviridae and have specific roles both within and outside of the replication complex, which differ between the different genera. These roles include subverting host proteins to generate replication organelles and inhibition of cellular functions (such as protein secretion) to influence virus replication efficiency and the host response to infection. In addition, 3A proteins are associated with the determination of host range. However, recent observations have challenged some of the roles assigned to 3A and suggest that other viral proteins may carry them out. In this review, we revisit the roles of 3A in the picornavirus life cycle. The 3AB precursor and mature 3A have distinct functions during viral replication and, therefore, we have also included discussion of some of the roles assigned to 3AB.

Antigenic properties of a novel vaccine strain for type Asia1 foot-and-mouth disease in pigs

Newly developed vaccine strains to prevent foot-and-mouth disease caused by the emerging serotype Asia1 virus were evaluated. To protect against the group (G)-VIII strain, which occurred recently, we produced an infectious cDNA clone of Asia1 Shamir cDNA (Asia1 Shamir-R). In addition, by adding a site 1 epitope of VP1 of the G-VIII lineage virus to this virus, we produced a new virus (Sham GVIII- EPI), and another virus(Sham GVIII-VP1) was replaced with that of G-VIII lineage in the VP1 region of Shamir. Test vaccines were produced using these three types of vaccine virus, and their immunogenicity and protection capabilities were evaluated in mice. Immunized mice were challenged with the Asia1 Shamir or G-VIII virus, and the results show that all the vaccines have similar protective effects. As they showed similar antigenicity, we chose the Shamir-R vaccine. Pigs maintained relatively high neutralizing antibody levels against homologous viruses of the Shamir and G-VII or G-VIII lineage three to four weeks after immunization. However, they formed relatively low levels of antibodies to G-IV and G-V viruses. In conclusion, we produced a vaccine candidate capable of protection against the G-VIII virus in the vaccine experiment for the type Asia1 serotype vaccine. This Shamir-R vaccine virus was found to protect against the viruses of the Asia1 genotype G-VII and G-VIII lineages, which occurred recently in Asia.

Evidence for multiple recombination events within foot-and-mouth disease viruses circulating in West Eurasia

Phylogenetic studies on foot-and-mouth disease viruses (FMDVs) circulating in the West Eurasian region have largely focused on the genomic sequences encoding the structural proteins that determine the serotype. The present study has compared near-complete genome sequences of FMDVs representative of the viruses that circulate in this region. The near-complete genome sequences (ca. 7,600 nt) were generated from multiple overlapping RT-PCR products. These amplicons were from FMDVs belonging to serotypes O, A and Asia-1, including members of the O-PanAsia-II and the A-Iran05 lineages, and of Group-II and Group-VII (Sindh-08) within serotype Asia-1, which are currently predominant and widespread in West Eurasia. These new sequences were analysed together with other sequences obtained from GenBank. Comparison of different regions of the FMDVs genomes revealed evidence for multiple, inter-serotypic, recombination events between FMDVs belonging to the serotypes O, A and Asia-1. It is concluded from the present study that dramatic changes in virus sequences can occur in the field through recombination between different FMDV genomes. These analyses provide information about the ancestry of the serotype O, A and Asia-1 FMDVs that are currently circulating within the West Eurasian region.

Single-cell analysis reveals the relevance of foot-and-mouth disease virus persistence to emopamil-binding protein gene expression in host cells

Foot-and-mouth disease virus (FMDV) infects host cells in either an acute or persistent manner. In this study, we examined the relevance of the establishment of FMDV persistence to the expression of the emopamil-binding protein (EBP) gene in 231 individual persistently infected baby hamster kidney (BHK-21) cells after passages 28, 38, and 68 (PI28, PI38, and PI68). At PI28, the stage at which persistent infection of FDMV becomes unstable, the percentage of cells carrying FMDV was 66.7%, while 80.2% of cells were EBP positive. Additionally, in 55.6% of the EBP-positive cells at PI28, EBP expression was upregulated approximately 149.9% compared to uninfected BHK-21 cells. This was the highest expression level among all cell passages measured. Interestingly, in a parallel experiment, the average EBP expression level in the whole cell population at PI28 was only slightly higher (108.2%) than that in uninfected BHK-21 cells. At PI38, 98.7% of the cells were positive for FMDV 3D (an RNA-dependent RNA polymerase enzyme gene), and its maximum expression level observed at this passage. The expression level of EBP in 78.2% of the total cells, however, was reduced significantly. At PI68, 95.8% of the cells were 3D positive, and the expression of both the EBP and 3D genes were at the lowest levels of all the passages. Our studies using single cells yielded data that are otherwise inaccessible a using whole cell population. These results suggest that the establishment of persistent infection by FMDV is a dynamic process that results from the continuous adaptation and coevolution of viruses and cells to reach an equilibrium.

Phylodynamics of foot-and-mouth disease virus O/PanAsia in Vietnam 2010-2014

Foot-and-mouth disease virus (FMDV) is endemic in Vietnam, a country that plays an important role in livestock trade within Southeast Asia. The large populations of FMDV-susceptible species in Vietnam are important components of food production and of the national livelihood. In this study, we investigated the phylogeny of FMDV O/PanAsia in Vietnam, reconstructing the virus' ancestral host species (pig, cattle or buffalo), clinical stage (subclinical carrier or clinically affected) and geographical location. Phylogenetic divergence time estimation and character state reconstruction analyses suggest that movement of viruses between species differ. While inferred transmissions from cattle to buffalo and pigs and from pigs to cattle are well supported, transmission from buffalo to other species, and from pigs to buffalo may be less frequent. Geographical movements of FMDV O/PanAsia virus appears to occur in all directions within the country, with the South Central Coast and the Northeast regions playing a more important role in FMDV O/PanAsia spread. Genetic selection of variants with changes at specific sites within FMDV VP1 coding region was different depending on host groups analyzed. The overall ratio of non-synonymous to synonymous nucleotide changes was greater in pigs compared to cattle and buffalo, whereas a higher number of individual amino acid sites under positive selection were detected in persistently infected, subclinical animals compared to viruses collected from clinically diseased animals. These results provide novel insights to understand FMDV evolution and its association with viral spread within endemic countries. These findings may support animal health organizations in their endeavor to design animal disease control strategies in response to outbreaks.

Evolutionary dynamics of foot-and-mouth disease virus O/ME-SA/Ind2001 lineage

Foot-and-mouth disease (FMD) virus serotype O Ind2001 lineage within the Middle East-South Asia topotype is the major cause of recent FMD incidences in India. A sub-lineage of Ind2001 caused severe outbreaks in the southern region of the country during 2013 and also reported for the first time from Libya. In this study, we conducted a detailed evolutionary analysis of Ind2001 lineage. Phylogenetic analysis of Ind2001 lineage based on maximum likelihood method revealed two major splits and three sub-lineages. The mean nucleotide substitution rate for this lineage was calculated to be 6.338×10(-3)substitutions/site/year (s/s/y), which is similar to those of PanAsian sub-lineages. Evolutionary time scale analysis indicated that the Ind2001 lineage might have originated in 1989. The sub-lineage Ind2001d that caused 2013 outbreaks seems to be relatively more divergent genetically from other Ind2001 sub-lineages. Seven codons in the VP1 region of Ind2001 were found to be under positive selection. Four out of 24 recent Ind2001 strains tested in 2D-MNT had antigenic relationship value of <0.3 with the serotype O vaccine strain indicating intra-epidemic antigenic diversity. Amino acid substitutions found in these minor variants with reference to antigenic diversity have been discussed. The dominance of antigenically homologous strains indicates absence of vaccine immunity in the majority of the affected hosts. Taken together, the evolution of Ind2001 lineage deviates from the strict molecular clock and a typical lineage evolutionary dynamics characterized by periodic emergence and re-emergence of Ind2001 and PanAsia lineage have been observed in respect of serotype O.

Molecular evolution of kobuviruses in cats

Aichi virus, a causative agent of human gastroenteritis, is one of a number of animal viruses belonging to the genus Kobuvirus within the family Picornaviridae. The kobuvirus genome encodes several structural and nonstructural proteins; the capsid proteins encoded by the VP1 gene are key immunogenic factors. Here, we used the VP1 region to determine substitution rates and the time to the most recent common ancestor (TMRCA) by comparing feline kobuvirus (FKoVs) sequences with kobuvirus sequences isolated from members of other species. The substitution rate for FKoVs was 1.29 × 10(-2 )substitutions/site/year (s/s/y) and the TMRCA was 5.3 years.

Phylodynamic reconstruction of O CATHAY topotype foot-and-mouth disease virus epidemics in the Philippines

Reconstructing the evolutionary history, demographic signal and dispersal processes from viral genome sequences contributes to our understanding of the epidemiological dynamics underlying epizootic events. In this study, a Bayesian phylogenetic framework was used to explore the phylodynamics and spatio-temporal dispersion of the O CATHAY topotype of foot-and-mouth disease virus (FMDV) that caused epidemics in the Philippines between 1994 and 2005. Sequences of the FMDV genome encoding the VP1 showed that the O CATHAY FMD epizootic in the Philippines resulted from a single introduction and was characterised by three main transmission hubs in Rizal, Bulacan and Manila Provinces. From a wider regional perspective, phylogenetic reconstruction of all available O CATHAY VP1 nucleotide sequences identified three distinct sub-lineages associated with country-based clusters originating in Hong Kong Special Administrative Region (SAR), the Philippines and Taiwan. The root of this phylogenetic tree was located in Hong Kong SAR, representing the most likely source for the introduction of this lineage into the Philippines and Taiwan. The reconstructed O CATHAY phylodynamics revealed three chronologically distinct evolutionary phases, culminating in a reduction in viral diversity over the final 10 years. The analysis suggests that viruses from the O CATHAY topotype have been continually maintained within swine industries close to Hong Kong SAR, following the extinction of virus lineages from the Philippines and the reduced number of FMD cases in Taiwan.

Foot-and-mouth disease: past, present and future

Foot-and-mouth disease (FMD) is a highly contagious disease of cloven-hoofed animals including cattle, pigs, sheep and many wildlife species. It can cause enormous economic losses when incursions occur into countries which are normally disease free. In addition, it has long-term effects within countries where the disease is endemic due to reduced animal productivity and the restrictions on international trade in animal products. The disease is caused by infection with foot-and-mouth disease virus (FMDV), a picornavirus. Seven different serotypes (and numerous variants) of FMDV have been identified. Some serotypes have a restricted geographical distribution, e.g. Asia-1, whereas others, notably serotype O, occur in many different regions. There is no cross-protection between serotypes and sometimes protection conferred by vaccines even of the same serotype can be limited. Thus it is important to characterize the viruses that are circulating if vaccination is being used for disease control. This review describes current methods for the detection and characterization of FMDVs. Sequence information is increasingly being used for identifying the source of outbreaks. In addition such information can be used to understand antigenic change within virus strains. The challenges and opportunities for improving the control of the disease within endemic settings, with a focus on Eurasia, are discussed, including the role of the FAO/EuFMD/OIE Progressive Control Pathway. Better control of the disease in endemic areas reduces the risk of incursions into disease-free regions.

Genetic basis of antigenic variation in foot-and-mouth disease serotype A viruses from the Middle East

Foot-and-mouth disease viruses (FMDV) from serotype A exhibit high antigenic diversity. Within the Middle East, a strain called A-Iran-05 emerged in 2003, and subsequently replaced the A-Iran-96 and A-Iran-99 strains that were previously circulating in the region. Viruses from this strain did not serologically match with the established A/Iran/96 vaccine, although most early samples matched with the older A22/Iraq vaccine. However, many viruses from this strain collected after 2006 had poor serological match with the A22/Iraq vaccine necessitating the development of a new vaccine strain (A/TUR/2006). More recently, viruses from the region now exhibit lower cross-reactivity with the A/TUR/2006 antisera highlighting the inadequacy of the serotype A vaccines used in the region. In order to understand the genetic basis of these antigenic phenotypes, we have determined the full capsid sequence for 57 Middle Eastern viruses isolated between 1996 and 2011 and analysed these data in context of antigenic relationship (r1) values that were generated using antisera to A22/Iraq and A/TUR/2006. Comparisons of capsid sequences identified substitutions in neutralising antigenic sites (1, 2 and 4), which either individually or together underpin these observed antigenic phenotypes.

Immunity of foot-and-mouth disease serotype Asia 1 by sublingual vaccination

Foot-and-mouth disease virus (FMDV) causes vesicular disease of cloven-hoofed animals, with severe agricultural and economic losses. Here we present study using a sublingual (SL) route with the killed serotype Asia 1 FMDV vaccine. Guinea pigs were vaccinated using a commercially available vaccine formulation at the manufacturer’s recommended full, 1/4, and 1/16 antigen doses. Animals were challenged with homologous FMDV Asia1 strain at various times following vaccination. All control guinea pigs exhibited clinical disease, including fever, viremia, and lesions, specifically vesicle formation in feet. Animals vaccinated with the 1/16 and 1/4 doses were protected after challenge at days 7, 28, and 35 post vaccination. These data suggest that effective protection against foot-and-mouth disease can be achieved with 1/16 of the recommended vaccine dose using SL vaccination, indicating that the sublingual route is an attractive alternative for the administration of the FMDV vaccine.

Comparative study of codon substitution patterns in foot-and-mouth disease virus (serotype O)

We compared genetic variations in the VP1 gene of foot-and-mouth disease viruses (FMDVs) isolated since 2000 from various region of the world. We analyzed relative synonymous codon usage (RSCU) and phylogenetic relationship between geographical regions, and calculated the genetic substitution patterns between Korean isolate and those from other countries. We calculated the ratios of synonymously substituted codons (SSC) to all observed substitutions and developed a new analytical parameter, EMC (the ratio of exact matching codons within each synonymous substitution group) to investigate more detailed substitution patterns within each synonymous codon group. We observed that FMDVs showed distinct RSCU patterns according to phylogenetic relationships in the same serotype (serotype O). Moreover, while the SSC and EMC values of FMDVs decreased according to phylogenetic distance, G + C composition at the third codon position was strictly conserved. Although there was little variation among the SSC values of 18 amino acids, more dynamic differences were observed in EMC values. The EMC values of 4- and 6-fold degenerate amino acids showed significantly lower values while most 2-fold degenerate amino acids showed no significant difference. Our findings suggest that different EMC patterns among the 18 amino acids might be an important factor in determining the direction of evolution in FMDV.

Evolutionary analysis of serotype A foot-and-mouth disease viruses circulating in Pakistan and Afghanistan during 2002-2009

Foot-and-mouth disease (FMD) is endemic in Pakistan and Afghanistan. Three different serotypes of the virus, namely O, A and Asia-1, are responsible for the outbreaks of this disease in these countries. In the present study, the nucleotide-coding sequences for the VP1 capsid protein (69 samples) or for all four capsid proteins (P1, seven representative samples) of the serotype A FMD viruses circulating in Pakistan and Afghanistan were determined. Phylogenetic analysis of the foot-and-mouth disease virus (FMDV) VP1-coding sequences from these countries collected between 2002 and 2009 revealed the presence of at least four lineages within two distinct genotypes, all belonging to the Asia topotype, within serotype A. The predominant lineage observed was A-Iran05 but three other lineages (a new one is named here A-Pak09) were also identified. The A-Iran05 lineage is still evolving as revealed by the presence of seven distinct variants, the dominant being the A-Iran05AFG-07 and A-Iran05BAR-08 sublineages. The rate of evolution of the A-Iran05 lineage was found to be about 1.2×10(-2) substitutions per nucleotide per year. This high rate of change is consistent with the rapid appearance of new variants of FMDV serotype A in the region. The A22/Iraq FMDV vaccine is antigenically distinct from the A-Iran05BAR-08 viruses. Mapping of the amino acid changes between the capsid proteins of the A22/Iraq vaccine strain and the A-Iran05BAR-08 viruses onto the A22/Iraq capsid structure identified candidate amino acid substitutions, exposed on the virus surface, which may explain this antigenic difference.