Skip to main content
  • Research Article
  • Open access
  • Published:

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

Abstract

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.

Introduction

Foot-and-mouth disease (FMD) is a highly transmissible viral disease of cloven hooved animals and is considered one of the most important diseases of livestock. Countries in Southeast and East Asia have varying levels of FMD endemicity, with Cambodia, Thailand, Laos, China and Vietnam having relatively high FMD incidence throughout the year [1]. Vietnam is the largest pig trader in Southeast Asia, so FMD control in this country is critical for the entire region [2]. The economic burden of FMD is substantial for large and small-scale pig producers and pig owners [3]. Cattle and water buffalo (Bubalus bubalis) have similar and important roles in agricultural practice in Vietnam [2, 4]. Both species are kept in varying degrees of intensity for dairy and meat production and are additionally used for draught purposes. Both species are often allowed to range freely for variable periods of time, and are frequently moved across country borders [5].

A study carried out in Vietnam, targeting areas with recent history of FMDV, found 22.3% seropositivity (to non-structural proteins) amongst subclinical buffalo and cattle sampled [6]. FMDV serotypes A and O currently circulate in Vietnam, while serotype Asia 1 has not been reported in Vietnam since 2008 [79]. FMDV serotype O is the most prevalent in the country; a recent study sequenced 71 serotype O viruses from samples collected between 2009 and 2013 in Vietnam, 65 of these viruses belonged to O/ME-SA/PanAsia lineage, while only 6 were classified as O/SEA/Mya-98 [6, 8]. However, in Vietnam from 2014–2016 it is believed that there has been a resurgence of O/SEA/Mya-98 (Dung and Long, unpublished data). Additionally, an incursion of O/ME-SA/Ind2001d lineage was reported for the first time in the country in 2016 [10].

Ruminants infected with FMDV may either clear the virus within 1–2 weeks after initial infection or develop a subclinical, persistent infection [1113]. The World Organization for Animal Health (OIE) defines persistent FMDV infection as the recovery of FMDV from oropharyngeal fluids >28 days post-infection (dpi). However, recent work has demonstrated that, under experimental conditions, cattle that clear the infection can be differentiated from those that develop persistent infection at 14 dpi for vaccinated animals and 21 dpi amongst non-vaccinated animals [14, 15]. These persistently infected animals are referred to as FMDV carriers [13, 16] and earlier studies have estimated that the proportion of infected cattle that become carriers range from 50 to 65% [14, 17]. The role of carriers in disease transmission amongst cattle has been extensively debated [1719]. Researchers have conducted numerous experimental studies, but have failed to detect transmission from carrier cattle to susceptible animals, while others have concluded that transmission from persistently infected cattle may occur to a very limited extent [20, 21].

There are intrinsic challenges to investigating FMDV transmission between herds and between different host species under field conditions in areas where animal movements and husbandry practices are variable and inconsistently regulated. Molecular techniques used to study the evolution of a pathogen through phylogenetic reconstruction can help understand the geographical spread, transmission between species and transmission from carrier and acutely infected animals, respectively when records about animal movements and contacts are not available [22, 23].

The objective of this study was to reconstruct the phylogeny, inferred inter-species transmission, and geographic spread of FMDV using VP1 coding sequences obtained from different host species, clinical stages, and locations in Vietnam. These analyses can help identify characteristics of the host species, geographic location, and disease status that are associated with specific changes in the viral genome and selection pressures.

Materials and methods

Data source

This investigation of FMDV phylodynamics in Vietnam was based on 125 FMDV VP1 coding region sequences [639 nucleotides (nt) total length]. The VP1 capsid protein comprises approximately 7.6% of the FMDV genome, and is commonly used for first-line phylogenetic analyses because it is known to be the most variable region of FMDV due to selective pressure on the immunogenic epitopes contained therein [24]. These sequences were obtained from three sources: (1) previous studies conducted by our group wherein 77 FMDVs collected in Vietnam were sequenced at the Plum Island Animal Disease Center, United States and identified as FMDV O/ME-SA topotype, PanAsia lineage as described in previous publications [6, 25], (2) eleven sequences generated by the OIE/FAO World Reference Laboratory for Foot-and-mouth disease (WRLFMD, The Pirbright Institute, UK), delivered directly to the Vietnam Department of Animal Health and deposited in GenBank for this study, and (3) 37 recently described, genetically related sequences that were retrieved from GenBank, 34 of which are from Vietnam, two from Kazakhstan and one from China. All sequences used in the study are currently available (see Additional file 1; Table 1). No animal experimentation or euthanasia was performed for the sake of completing this study. Because several VP1 sequences were shorter than complete length of the protein-coding segment (639 nt), we trimmed the alignment to the first 621 nt, to have consistent data for phylogenetic reconstruction.

Table 1 Bayesian stochastic search variable selection analysis results

Phylogenetic reconstruction using divergence time estimation

We reconstructed the phylogeny of FMDV and estimated divergence times. Sequences were aligned using the MUSCLE algorithm [26]. To estimate the best codon partition and best substitution model, we analyzed the sequences using Partition Finder [27] and selected the best partition scheme based on the Bayesian Information Criterion (BIC). The purpose of this selection is to identify the best codon partition scheme and the best nt substitution model for each partition.

Using the best partition scheme and corresponding substitution model, we reconstructed the phylogeny of O/ME-SA/PanAsia FMD virus using a Bayesian statistical approach (Bayesian Evolutionary Analysis by Sampling Trees), implemented in BEAST 1.8.2 [28]. The sampling dates were specified to estimate times of divergence. We used the lognormal uncorrelated relaxed clock model, and the Coalescent Bayesian Skyline tree prior. The analysis was run for 2 × 108 iterations within a web-based platform with access to computational resources available in CIPRES [29]. Convergence of the chain was assessed using Tracer 1.6, by visualizing traces of parameters of trees sampled and confirming that mixing of the chain had been achieved so that the effective sample size of all parameters was >200 [30]. From all trees sampled, the maximum clade credibility tree (MCC) was annotated and depicted using FigTree 1.4.2 [31]. The initial 10% sampled trees were discarded as burn-in. Time to most recent common ancestor (tMRCA) of all nodes and 95% highest posterior densities (95% HPD) were obtained from the MCC tree.

Phylogeographic analysis and ancestral character state reconstruction

We used discrete ancestral character state reconstruction to estimate the viral history, specifying traits according to 3 different characteristics: the host species, the clinical status, and the location. To reconstruct the host species, we defined three discrete characters: cattle, pig or buffalo. For the clinical status, we used two characters: sequences described as “Clinical” were obtained from samples derived from vesicular lesions of animals during outbreaks of clinical FMD in Vietnam, whereas “Carrier” sequences were obtained from oropharyngeal fluid (probang) from subclinical cattle and buffalo identified through active surveillance as previously described [6]. To infer ancestral states with respect to geographical location we categorized each viral sequence assigning one of the 8 different geographical regions defined within Vietnam: Northeast, Northwest, Red River Delta, North Central Coast, South Central Coast, Central Highlands, Southeast, and the Mekong River Delta.

For each host species, clinical disease and geographic location traits, we reconstructed the inferred character state for each node within the phylogenetic tree. We analyzed the inferred transmission rates between character states using an asymmetric model for the discrete traits and estimated the significance of the network with Bayesian stochastic search variable selection (BSSVS), which tests the hypothesis of non-zero transmission rates between discrete characters [22]. Statistical support was assessed using Bayes Factor (BF) for discrete traits implemented in SPREAD1.0.6 [32], we considered BF > 3 as significant non-zero transmission. The analysis was carried out in BEAST 1.8.2 and the number of iterations and assessment of chains were performed as described above. Additionally, we obtained the count of character transitions (“jumps”) from all trees (excluding initial burning). The 95% high-density interval (95% HDI) of the values collected for each character state change that had a non-zero rate BF > 3.0 was computed using HDInterval package in R [33].

Evolutionary selection of nucleotide sites across different species and between clinically versus persistently infected animals

The VP1 coding region of the 122 FMDV O/ME-SA/PanAsia sequences described above and the additional related sequences (n = 3) from Kazakhstan and China were analyzed for positive and negative selection. Sequence groups from different species (cattle, pigs and buffalo-excluding samples from persistently infected animals) as well as clinical or persistent status were analyzed independently to study differences in viral selection by nt site and overall selection. The mean ratio of non-synonymous (dN) and synonymous changes (dS) (global ω = dN/dS ratio) was computed using the single-likelihood ancestor counting (SLAC) method, and 95% confidence intervals estimated from the data and the likelihood profile [34]. Individual site selection was also computed by the fixed effects likelihood (FEL) method and by the random effects likelihood (REL; which was used for the species groups but not for the clinical/carrier groups due to the size restriction -number of sequences- of the analysis) [34]. Only unique sequences were used in this analysis (i.e. identical sequences were removed). The HKY85 nt substitution model was used. A site was considered positively or negatively selected if identified by at least one of the methods described (dN-dS > 0 with a p value cutoff 0.1 for SLAC, 0.1 for FEL and 50 BF for REL). We displayed the dN-dS values computed for each codon position and indicated the statistically significant positively selected sites. The analyses were performed in the HyPhy2.2.1 software package [34, 35].

Results

Divergence time estimation and ancestral character reconstruction

The partition scheme and substitution models selected were the K80 + I for codon positions 1 + 2, and the HKY + G substitution model, for codon position 3. The mean substitution/site/year for VP1 coding segment the O/ME-SA/PanAsia phylogeny reconstructed was 1.66 × 10−2 (95% HPD 1.21–2.11 × 10−2).

Ancestral state character reconstruction

Host species

The MCC tree depicting the reconstruction based on FMDV host species is presented in Figure 1. The current FMD O/ME-SA/PanAsia viruses circulating in Vietnam diverged into two different clades in June 2010 (95% HPD March 2010 to September 2010). One of these lineages (Figure 1, branch A), initially found in all three hosts species, diverged into a cluster of viruses that subsequently were found predominantly in pigs. Transitions between host species character state inferred from the MCC tree (and 95% HDI) occurred from cattle to pigs (7 events 95% HDI 5–10), pigs to cattle (2 events 95% HDI 1–3), and from cattle to buffalo (8 events 95% HDI 7–11). Results of the analysis to estimate significant transmission rates between species, BSSVS, are shown in Table 1. With respect to host species, non-zero inferred transmission rates were detected from cattle to pig, from cattle to buffalo and from pig to cattle.

Figure 1
figure 1

Maximum clade credibility tree of FMDVs in Vietnam and related viruses from China and Kazakhstan between 2010 and 2014. The color of tree branches and nodes indicates the ancestral host species for the reconstructed phylogeny. Clades A and B represent the two main O/ME-SA/PanAsia sublineages that have diverged recently in Vietnam. Characteristics of the sampled viruses (clinical stage: C., host species: H., and location: L.) are indicated in colored columns aligned to the right of the tree, color coding is indicated in the metadata legend.

Clinical status

The MCC tree showing the ancestral reconstruction of clinically affected animals (“Clinical”) and subclinically infected animals (“Carriers”) is shown in Figure 2. In the upper main branch of the tree (Figure 2, branch A), there are two viruses from persistently infected animals (KT153098-O/VIT/12/2012pro and KT153128-O/VIT/25/2012pro) likely originating from previous outbreaks. However the tMRCA with the most closely genetically related ‘clinical’ viruses was relatively long (compared to tMRCA of related viruses), suggesting that they had already diverged approximately one year (1.32 and 0.53 years for viruses O/VIT/12/2012pro and O/VIT/25/2012pro respectively) from their ancestral “clinical” sample sequences. The lower branch of the tree (Figure 2, branch B) includes closely related viruses from 9 different persistently infected animals. The ancestral reconstruction of the set of viruses analyzed suggests that these viruses may have initiated an outbreak with clinical disease, which includes several closely related viruses (Figure 2, *).

Figure 2
figure 2

Maximum clade credibility FMDV O/ME-SA/PanAsia viruses collected in Vietnam (and additional 3 sequences from China and Kazakhstan) between 2010 and 2014. Nodes and branches of the trees are colored according to the clinical stage reconstructed in the phylogeny. The ancestral reconstruction of the viruses analyzed suggests 1 instance where outbreak viruses may have originated from carriers (*). Clades A and B represent the two main O/ME-SA/PanAsia sublineages that have diverged recently in Vietnam. Characteristics of the sampled viruses (clinical stage: C., host species: H., and location: L.) are indicated in colored columns aligned to the right of the tree (color coding legend for the columns is indicated in Figure 1).

Inferred non-zero transmission rates estimated by BSSVS analysis (Table 1), which is used to detect significant non-zero transmission from carrier to clinical and from clinical to carrier categories.

Phylogeographic analysis

The phylogeographic reconstruction of the viruses is shown in Figure 3. The common ancestor of all O/ME-SA/PanAsia sequences included in this study was located in the northeast region. Putative movements of this lineage from the Northeast of Vietnam into China and from the Northwest Region of Vietnam into Kazakhstan can be observed in Figure 3—branch A and may have occurred in December of 2010. The phylogeographic analysis demonstrates that the upper part of Figure 3—branch A, which contains a cluster of pig-derived viruses (Figure 2, branch A), may have initially originated from the Northeast and then spread into the Northwest, and the Red River Delta Regions.

Figure 3
figure 3

Maximum clade credibility FMDV O/ME-SA/PanAsia viruses collected in Vietnam (and additional 3 sequences from China and Kazakhstan) between 2010 and 2014. Nodes and branches of the trees are colored according to the location with the higher posterior probability in the reconstructed phylogeny. Ancestors of the two divergent clusters are inferred to exist in the Northeast. Clades A and B represent the two main O/ME-SA/PanAsia sublineages that have diverged recently in Vietnam. Characteristics of the sampled viruses (clinical stage: C., host species: H., and location: L.) are indicated in colored columns aligned to the right of the tree, color coding is indicated in the metadata legend.

To visualize the significant inferred transmission rates (BF > 3) between the geographical regions we overlaid these results on a map of Vietnam (Figure 4). The most significant inferred transmissions were detected from Northeast and South Central Coast to their corresponding adjacent regions, as well as transmission from Mekong River Delta into Southeast region. Inferred disease spread occurred in northern and southern directions, however there was a general trend of higher BF for southbound transfer of viruses. BF-inferred transmission is also shown in the heatmap displayed in Figure 4, and it further evidences that the South Central Coast and the Northeast are the regions from where the viruses are more frequently spread into other regions.

Figure 4
figure 4

Results from the Bayesian stochastic search variable selection of the Phylogeographic reconstruction of O/ME-SA/PanAsia in Vietnam regions, and related sequences from Kazakhstan and China. Only Bayes factor > 3 are represented as arrows as significant non-zero transmission of O/ME-SA/PanAsia. Most significant transmissions were inferred for some adjacent regions, although inferred transmission between some distant regions was also statistically supported. The heatmap in the lower area of the figure depicts the magnitude of the statistical support (Bayes factor) for transmission rate between geographic regions in Vietnam. This heatmap allows visualizing that more transmission occurred from South Central Coast and the Northeast regions into other areas, whereas the Red River Delta and Mekong Delta were the ones with more incoming transmission from other regions.

Evolutionary selection of sites in different species and clinically versus persistently infected animals

Results of the global dN/dS ratio (ratio of non-synonymous to synonymous changes) estimated by category are shown in Table 2. Considering only viruses from clinical outbreaks, pig-derived viruses had a higher overall positive selection ratio (dN/dS) compared to buffalo and cattle; however, only the difference with cattle was statistically significant at p < 0.05 (based on non-overlapping 95% confidence intervals). The extent of positive selection was similar in buffalo and cattle. Similarly, viruses collected from persistently infected animals and from clinically affected ones had almost the same dN/dS values (Table 2).

Table 2 Results of the global dN/dS ratio estimated for each of the categories and corresponding 95% confidence interval

Results of the specific site selection for every codon (dN-dS) in all categories are shown in Figures 5A and B. No specific sites under positive selection were found to be statistically significant when analyzing viral sequences collected from cattle and buffalo, whereas few sites were detected when analyzing viral sequences from pig (site numbers: 1, 152, 153 and 172). Several statistically significant negatively selected sites were found for cattle (23 sites), pigs (28 sites), and buffalo (7 sites). Overall positive selection was mainly found within the known antigenic sites (GH loop, BC loop). Two statistically significant positively selected sites were found in the GH loop in pigs, compared to none in buffalo and cattle groups.

Figure 5
figure 5

Site selection (dN-dS) by VP1 coding region in alignments of FMDV sequences. A Site selection (dN-dS) results per VP1 coding region in alignments of cattle viral sequences (excluding persistently infected), buffalo (excluding persistently infected) and pigs. Values >0 represent positive selection. Bars colored in red indicate sites where selection is statistically significant. Grey shaded areas correspond to known antigenic sites within VP1: BC loop (sites 43–45), GH loop (sites 238–254), and C terminus (sites 200–207). Pink-shaded sites correspond to the RGD integrin-binding motif. B Site selection (dN-dS) results per VP1 coding region in alignments including of all viral sequences, “clinical” (cattle and buffalo) viral sequences and “carrier” (cattle and buffalo) infected animals.

Although the global positive selection ratio was the same in carrier and clinically diseased animals, when dissecting the specific positively selected amino acid sites in the antigenic domains, carrier animals had a higher number of statistically significant positively selected sites compared to selection in viruses from clinical samples: Four statistically significant positively selected sites in the GH loop (near the RGD motif) and one in the BC loop were found in carriers, whereas two positively sites in the GH loop and two sites in the BC loop without statistical significance were found in clinical samples from cattle and buffalo (Figure 5B). Three additional sites of significant positive selection were also identified in viruses from carriers, but in non-antigenic regions. When analyzing all viral sequences together, several positively selected (non significant) sites were found within the GH-loop and the BC loop, whereas numerous significantly negatively selected sites were found throughout VP1. Site 73 was consistently negatively selected (higher expected number of synonymous changes) in all group categories. Only one site (154) was detected as both positive and negatively selected depending on the category; while positive selection was detected within the carrier animals category, this site was identified as a negatively selected for pig viral sequences.

Discussion

In the current study, we analyzed phylodynamics of FMDV O/ME-SA/PanAsia viruses recovered from livestock in Vietnam between 2010 and 2014. The results presented herein provide a novel overview of comparative viral evolution of field samples collected from different host species, locations, and different clinical stages of infection. These data suggest differences in viral evolution within distinct animal groups, which may contribute to understanding of the mechanisms of maintenance, emergence and spread of viruses across Vietnam.

Using viruses sampled between 2010 and 2014 from different species, stages of infection and provinces of Vietnam, we found two main sublineages of FMDV O/ME-SA/PanAsia. One of these sublineages has been most frequently recovered from pigs since 2011 while the other main sublineage has been most frequently found in cattle. Some FMDV strains may have a predilection to certain host species; this phenomenon has been extensively described with porcinophilic FMDV serotype O of the Cathay topotype [36, 37]. However, the data presented herein is not sufficient to determine limited host range of the mostly pig-derived FMDV subclade found in this study. This is because several other factors, including sampling bias, could have determined the apparent species-specificity. The hypothesis of limited host range could be confirmed by follow up studies in tissue culture or in vivo to determine if these viruses have evolved to a specific pig host predilection.

The substitution rate calculated herein for VP1 of FMDV O/ME-SA/PanAsia circulating in Vietnam between 2009 and 2014 was similar to a previous estimate of serotype O Cathay topotype (1.06 × 10−2) [38]. However, the rate estimated for this study was higher than those computed across other serotype O viruses estimated at 6.65 × 10−3 [39], 6.34 × 10−3 [40], or 4.81 × 10−3 (serotype O collected globally between 1939 and 2010) [41] substitutions/site/year. The high mutation rate found in our analyses may be a result of capturing mutations that do not become fixed in a population, as a consequence of sampling a high number of outbreaks in a relatively short period of time. In contrast, isolated sequences from different regions sampled over longer periods may reflect fixed mutations, which may result in computing a different substitution rate. However, the specific underlying mechanisms driving these differences across studies using different geotemporal conditions remain undetermined.

Inferred transmission as deduced by ancestral character reconstruction suggests that cattle have a relevant role in inter-species dissemination. Various factors may contribute to this phenomenon including subclinical infection in vaccinated cattle and a longer duration of FMDV infectiousness in cattle compared to other species [12, 13]. Similarly, our study also provides evidence of inferred transmission of viruses from pigs to cattle. This may be explained by the large pig population in the country, and husbandry practices including the highly common practice of co-mingling pigs with other species in rural households and villages. In contrast, based on the sequences obtained, we did not identify putative transmission from buffalo to either cattle or pigs. Although transmission from cattle to buffalo, and buffalo to cattle is known to occur [42, 43], to our knowledge, comparative quantitative estimates of transmission rates of these two species have not been estimated. We detected several instances of transfer of viruses from cattle to buffalo, but none from buffalo to cattle. The lack of detection may either represent lesser ability of buffalo to transmit the disease to cattle, differential ranging and transport of buffalo and cattle within and between different areas in Vietnam, or it may be an effect of sampling bias.

Phylogenetic analyses in this study, suggest that there are instances in which viruses found in carrier animals may have been the ancestors of viruses that later caused outbreaks. However, it is also possible that undersampling of outbreak sequences around the period of the apparent carrier-to-clinical transfer may have created sampling bias which influenced these results (Figure 2). Transmission from carriers to susceptible animals has been either low or inexistent in controlled experiments. However, the potential epidemiological role of the FMDV carrier state in maintaining FMDV and being the source of new outbreaks is still highly controversial [13, 17, 44]. Quantitatively, it has been estimated that the transmission rate from acutely (clinically) infected animals is more than 500 times that of the carrier [20]. However even though the probability is exceedingly low, it would only require one successful transmission event within the thousands of contacts to start an outbreak and thus have a substantial impact.

Overall, the number of outbreaks and potential infected animals included herein is low compared to the actual number of outbreaks and infections that occurred in the field. Thus, such findings must be interpreted conservatively. Because of the inability to obtain every relevant virus from the field, analyses derived from studies such as this, contain intrinsic sampling bias and must be interpreted as generating, rather than confirming hypotheses. Similarly, these phylogenetic analyses were carried out using only VP1 coding sequence of FMDV, which is 639 nt long of the ~7000 nt ORF length. This genetic segment is the most variable within FMDV genome and contains the most relevant antigenic sites. Full genome sequences in sufficient quantities to enable these studies of phylodynamics of Vietnam, were not available at the time of this work.

Geographic spread of FMD within Vietnam inferred by viral sequence suggests that transmission occurs between several different regions of the country. However transmission seems to be more frequent from South Central Coast and Northeast into other parts of the country. As expected, the greatest extent of transmission was to immediately adjacent regions, but also to distant parts of the country. These findings are partially consistent with previous studies conducted by the OIE, which assessed the livestock movement in Vietnam and neighboring countries [5]. That report describes that substantial ruminant movement occurs from central areas of the country into northern and southern regions. Within the northern regions, this flow tends to be towards northeastern region, whereas in the south, animal movement converges into northern areas of the Mekong River Delta region. In contrast, pig movement occurs mostly from northern into southern areas of Vietnam and from northern areas into China. Transmission occurring in different geographic directions inferred by our analysis reflects the combination of these reported pig and ruminant species movement. The movement of large ruminants and pigs into China, however, is supported by our current inferred viral transmission data. Additionally, animal movement from Thailand, Cambodia and Laos into Vietnam occurs frequently [5].

Selection pressure may be an important driver of viral evolution [45]. It is biologically plausible that selection pressures differ in different host species, or through different phases of infection within the same host (e.g. acute clinical infection in contrast to persistent infection). Viral evolution comparing different groups of infected animals has not been extensively explored previously. Potential specific molecular changes in the FMDV genome that establish persistent infection have been reported in few controlled studies [17, 21, 46]. Here, we found that the global positive selection of VP1 tended to be similar in carrier animals compared to viruses recovered from clinical outbreaks. Although the value of the global positive selection rate was almost equal in carrier and outbreak groups, the quantity of individual (statistically significant) positive selected codon sites was higher in carrier animals compared with those with acute clinical infection, especially in the antigenic regions.

Because animal movement and mixing animals from different origins is a common practice in Vietnam, persistently infected animals may play an important role in disease spread. Understanding specific molecular (genomic and/or antigenic) changes can eventually help to understand the mechanisms through which FMDV persistence is established and maintained. Although analyses of our data demonstrate potentially relevant trends, results should be interpreted carefully with mindfulness of the relatively low quantity of samples from carrier animals (n = 12).

Additionally, we found that statistically significant positive selection was greatest in pigs compared to buffalo and cattle, and that these sites were found mostly in regions that code for the known antigenic domains. This suggests that the viral evolution in the pig population may provide an important contribution to antigenic diversity and strain emergence compared to other species. Previous work has shown that most of the genomic variation occurs early in the course of infection, when there is greater viral replication [47].

In conclusion, the current study suggests that inferred virus transmission patterns in Vietnam may differ depending on the host species and clinical status of infected hosts. This may be related to differences in viral selection between species and in the persistently infected animals compared to clinically affected individuals. These differences can help to elucidate viral evolution within-host, across host species, and within populations. Application and combination of the methodologies described herein to study specific aspects of FMDV evolution may help to gain new knowledge in our understanding of FMDV, ultimately contributing to disease control and eradication in endemic countries.

Abbreviations

BF:

Bayes factor

BSSVS:

Bayesian stochastic search variable selection

dN:

substitution rates at non-synonymous sites

dpi:

days post-infection

dS:

substitution rates at synonymous sites

FMD:

foot-and-mouth disease

FMDV:

foot-and-mouth disease virus

MCC:

maximum clade credibility

Nt:

nucleotide

OIE:

World Organization for Animal Health

tMRCA:

time to most recent common ancestor

References

  1. Sumption K, Rweyemamu M, Wint W (2008) Incidence and distribution of foot-and-mouth disease in Asia, Africa and South America; combining expert opinion, official disease information and livestock populations to assist risk assessment. Transbound Emerg Dis 55:5–13

    Article  CAS  PubMed  Google Scholar 

  2. Di Nardo A, Knowles NJ, Paton DJ (2011) Combining livestock trade patterns with phylogenetics to help understand the spread of foot and mouth disease in sub-Saharan Africa, the Middle East and Southeast Asia. Rev Sci Tech 30:63–85

    PubMed  Google Scholar 

  3. Pham HT, Antoine-Moussiaux N, Grosbois V, Moula N, Truong BD, Phan TD, Vu TD, Trinh TQ, Vu CC, Rukkwamsuk T, Peyre M (2016) Financial impacts of priority swine diseases to pig farmers in Red River and Mekong River Delta, Vietnam. Transbound Emerg Dis. doi:10.1111/tbed.12482

    PubMed  Google Scholar 

  4. Gleeson LJ (2002) A review of the status of foot and mouth disease in South-East Asia and approaches to control and eradication. Rev Sci Tech 21:465–475

    CAS  PubMed  Google Scholar 

  5. Smith P, Bourgeois Lüthi N, Huachun L, Naing Oo K, Phonvisay A, Premashthira S, Abila R, Widders P, Kukreja K, Miller C (2015) Movement pathways and market chains of large ruminants in the Greater Mekong Sub-region. Report of the World Organisation for Animal Health (OIE). Sub-regional representation for South-East Asia, Thailand. http://www.rr-asia.oie.int/fileadmin/SRR_Activities/documents/movement.pdf. Accessed Jan 2017

  6. de Carvalho Ferreira HC, Pauszek SJ, Ludi A, Huston CL, Pacheco JM, Le VT, Nguyen PT, Bui HH, Nguyen TD, Nguyen T, Nguyen TT, Ngo LT, Do DH, Rodriguez L, Arzt J (2017) An integrative analysis of foot-and-mouth disease virus carriers in Vietnam achieved through targeted surveillance and molecular epidemiology. Transbound Emerg Dis 64:547–563

    Article  PubMed  Google Scholar 

  7. Brito BP, Rodriguez LL, Hammond JM, Pinto J, Perez AM (2015) Review of the global distribution of foot-and-mouth disease virus from 2007 to 2014. Transbound Emerg Dis 64:316–332

    Article  PubMed  Google Scholar 

  8. Le VP, Nguyen T, Park JH, Kim SM, Ko YJ, Lee HS, Nguyen VC, Mai TD, Do TH, Cho IS, Lee KN (2010) Heterogeneity and genetic variations of serotypes O and Asia 1 foot-and-mouth disease viruses isolated in Vietnam. Vet Microbiol 145:220–229

    Article  PubMed  Google Scholar 

  9. Le PV, Vu TTH, Duong HQ, Than VT, Song D (2016) Evolutionary phylodynamics of foot-and-mouth disease virus serotypes O and A circulating in Vietnam. BMC Vet Res 12:269

    Article  PubMed  PubMed Central  Google Scholar 

  10. WRLFMD (2015) Reference laboratory reports. In: The Fao World Reference Laboratory for Foot-and-mouth Disease. http://www.wrlfmd.org/ref_labs/fmd_ref_lab_reports.htm. Accessed 15 Sept 2016

  11. Mardones F, Perez A, Sanchez J, Alkhamis M, Carpenter T (2010) Parameterization of the duration of infection stages of serotype O foot-and-mouth disease virus: an analytical review and meta-analysis with application to simulation models. Vet Res 41:45

    Article  PubMed  PubMed Central  Google Scholar 

  12. Alexandersen S, Zhang Z, Donaldson AI, Garland AJ (2003) The pathogenesis and diagnosis of foot-and-mouth disease. J Comp Pathol 129:1–36

    Article  CAS  PubMed  Google Scholar 

  13. Arzt J, Juleff N, Zhang Z, Rodriguez LL (2011) The pathogenesis of foot-and-mouth disease I: viral pathways in cattle. Transbound Emerg Dis 58:291–304

    Article  CAS  PubMed  Google Scholar 

  14. Stenfeldt C, Eschbaumer M, Rekant SI, Pacheco JM, Smoliga GR, Hartwig EJ, Rodriguez LL, Arzt J (2016) The foot-and-mouth disease carrier state divergence in cattle. J Virol 90:6344–6364

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Stenfeldt C, Belsham GJ (2012) Detection of foot-and-mouth disease virus RNA in pharyngeal epithelium biopsy samples obtained from infected cattle: investigation of possible sites of virus replication and persistence. Vet Microbiol 154:230–239

    Article  CAS  PubMed  Google Scholar 

  16. Sutmoller P, McVicar JW, Cottral GE (1968) The epizootiological importance of foot-and-mouth disease carriers. I. Experimentally produced foot-and-mouth disease carriers in susceptible and immune cattle. Arch Gesamte Virusforsch 23:227–235

    Article  CAS  PubMed  Google Scholar 

  17. Salt J (2004) Persistence of Foot-and-mouth Disease Virus. In: Sobrino F, Domingo E (eds) Foot-and-mouth disease current perspectives. Horizon Bioscience, Wymondham, pp 103–144

    Google Scholar 

  18. Kitching RP (2002) Identification of foot and mouth disease virus carrier and subclinically infected animals and differentiation from vaccinated animals. Rev Sci Tech 21:531–538

    CAS  PubMed  Google Scholar 

  19. Garland AJ, de Clercq K (2011) Cattle, sheep and pigs vaccinated against foot and mouth disease: does trade in these animals and their products present a risk of transmitting the disease? Rev Sci Tech 30:189–206

    CAS  PubMed  Google Scholar 

  20. Tenzin Dekker A, Vernooij H, Bouma A, Stegeman A (2008) Rate of foot-and-mouth disease virus transmission by carriers quantified from experimental data. Risk Anal 28:303–309

    Article  CAS  PubMed  Google Scholar 

  21. Parthiban AB, Mahapatra M, Gubbins S, Parida S (2015) Virus excretion from foot-and-mouth disease virus carrier cattle and their potential role in causing new outbreaks. PLoS One 10:e0128815

    Article  PubMed  PubMed Central  Google Scholar 

  22. Lemey P, Rambaut A, Drummond AJ, Suchard MA (2009) Bayesian phylogeography finds its roots. PLoS Comput Biol 5:e1000520

    Article  PubMed  PubMed Central  Google Scholar 

  23. Baele G, Suchard MA, Rambaut A, Lemey P (2017) Emerging concepts of data integration in pathogen phylodynamics. Syst Biol 66:e47–e65

    PubMed  Google Scholar 

  24. Jackson T, King AM, Stuart DI, Fry E (2003) Structure and receptor binding. Virus Res 91:33–46

    Article  CAS  PubMed  Google Scholar 

  25. Pauszek SJ, Eschbaumer M, Brito B, de Carvalho Ferreira HC, Vu LT, Phuong NT, Hoang BH, Tho ND, Dong PV, Tung N, Long NT, Dung DH, Rodriguez LL, Arzt J (2016) Site-specific substitution (Q172R) in the VP1 protein of FMDV isolates collected from subclinical field cases in Vietnam. Virol Rep 6:90–96

    Article  Google Scholar 

  26. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32:1792–1797

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Lanfear R, Calcott B, Ho SY, Guindon S (2012) Partitionfinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol 29:1695–1701

    Article  CAS  PubMed  Google Scholar 

  28. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7:214

    Article  PubMed  PubMed Central  Google Scholar 

  29. Miller MA, Pfeiffer W, Schwartz T (2010) Creating the CIPRES Science Gateway for inference of large phylogenetic trees. Gateway Computing Environments Workshop (GCE), New Orleans

    Book  Google Scholar 

  30. Rambaut A, Drummond A (2007) Tracer version 1.4, Computer program and documentation distributed by the author. http://tree.bio.ed.ac.uk/software/tracer/. Accessed 15 Sept 2016

  31. Rambault A (2006–2009) Fig Tree Tree Figure Drawing Tool, Version 13.1. http://tree.bio.ed.ac.uk/software/figtree/. Accessed 15 Sept 2016

  32. Bielejec F, Rambaut A, Suchard MA, Lemey P (2011) SPREAD: spatial phylogenetic reconstruction of evolutionary dynamics. Bioinformatics 27:2910–2912

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Meredith M, Kruschke J (2016) HDInterval: highest (posterior) density intervals. R package version 0.1.3. https://CRAN.R-project.org/package=HDInterval. Accessed 3 Mar 2017

  34. Kosakovsky Pond SL, Frost SDW (2005) Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol 22:1208–1222

    Article  PubMed  Google Scholar 

  35. Kosakovsky Pond SL, Frost SDW, Muse SV (2005) HyPhy: hypothesis testing using phylogenies. Bioinformatics 21:676–679

    Article  Google Scholar 

  36. Oem JK, Yeh MT, McKenna TS, Hayes JR, Rieder E, Giuffre AC, Robida JM, Lee KN, Cho IS, Fang X, Joo YS, Park JH (2008) Pathogenic characteristics of the Korean 2002 isolate of foot-and-mouth disease virus serotype O in pigs and cattle. J Comp Pathol 138:204–214

    Article  CAS  PubMed  Google Scholar 

  37. Beard CW, Mason PW (2000) Genetic determinants of altered virulence of Taiwanese foot-and-mouth disease virus. J Virol 74:987–991

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Di Nardo A, Knowles NJ, Wadsworth J, Haydon DT, King DP (2014) Phylodynamic reconstruction of O CATHAY topotype foot-and-mouth disease virus epidemics in the Philippines. Vet Res 45:90

    Article  PubMed  PubMed Central  Google Scholar 

  39. Jamal SM, Ferrari G, Ahmed S, Normann P, Belsham G (2011) Genetic diversity of foot-and-mouth disease virus serotype O in Pakistan and Afghanistan, 1997-2009. Infect Genet Evol 11:1229–1238

    Article  PubMed  Google Scholar 

  40. Subramaniam S, Mohapatra JK, Sharma GK, Biswal JK, Ranjan R, Rout M, Das B, Dash BB, Sanyal A, Pattnaik B (2015) Evolutionary dynamics of foot-and-mouth disease virus O/ME-SA/Ind2001 lineage. Vet Microbiol 178:181–189

    Article  PubMed  Google Scholar 

  41. Yoon SH, Lee KN, Park JH, Kim H (2011) Molecular epidemiology of foot-and-mouth disease virus serotypes A and O with emphasis on Korean isolates: temporal and spatial dynamics. Arch Virol 156:817–826

    Article  CAS  PubMed  Google Scholar 

  42. Madhanmohan M, Yuvaraj S, Nagendrakumar SB, Srinivasan VA, Gubbins S, Paton DJ, Parida S (2014) Transmission of foot-and-mouth disease virus from experimentally infected Indian buffalo (Bubalus bubalis) to in-contact naïve and vaccinated Indian buffalo and cattle. Vaccine 32:5125–5130

    Article  CAS  PubMed  Google Scholar 

  43. Gomes I, Ramalho AK, de Mello PA (1997) Infectivity assays of foot-and-mouth disease virus: contact transmission between cattle and buffalo (Bubalus bubalis) in the early stages of infection. Vet Rec 140:43–47

    Article  CAS  PubMed  Google Scholar 

  44. Alexandersen S, Mowat N (2005) Foot-and-mouth disease: host range and pathogenesis. Curr Top Microbiol Immunol 288:9–42

    CAS  PubMed  Google Scholar 

  45. Domingo E, Sheldon J, Perales C (2012) Viral quasispecies evolution. Microbiol Mol Biol Rev 76:159–216

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Kopliku L, Relmy A, Romey A, Gorna K, Zientara S, Bakkali-Kassimi L, Blaise-Boisseau S (2015) Establishment of persistent foot-and-mouth disease virus (FMDV) infection in MDBK cells. Arch Virol 160:2503–2516

    Article  CAS  PubMed  Google Scholar 

  47. Malirat V, De Mello PA, Tiraboschi B, Beck E, Gomes I, Bergmann IE (1994) Genetic variation of foot-and-mouth disease virus during persistent infection in cattle. Virus Res 34:31–48

    Article  CAS  PubMed  Google Scholar 

Download references

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JA, LR, BB, NTL, DHD, PVD, and SP conceived the study. SP, HF, PQM, LTV, NTP, BHH, NDT, ME, DK and NK performed the laboratory analysis, BB performed the sequence and statistical analyses, LTV, NDT, BHH, NTP, DK and NK coordinated and performed the sample collection. BB, JA, CS and SP wrote the manuscript. All authors read and approved the manuscript.

Acknowledgements

This research was funded in part by ARS-CRIS Project 1940-32000-057-00D and through an interagency agreement with the Science and Technology Directorate of the U.S. Department of Homeland Security under Award Number HSHQDC-12-X-0060. Additional funding was provided by the U.S. Department of State, Biosecurity Engagement Program through the USDA, ARS Office of International Research Programs and the Cooperative Biological Engagement Program of the U.S. Department of Defense, Defense Threat Reduction Agency. We acknowledge Ethan J. Hartwig and George R. Smoliga for expert technical support. Barbara Brito, Carolina Stenfeldt, Michael Eschbaumer, and Helena C. de Carvalho Ferreira are the recipients of a Plum Island Animal Disease Center Research Participation Program fellowship, administered by the Oak Ridge Institute for Science and Education (ORISE) through an interagency agreement between the U.S. Department of Energy (DOE) and the U.S. Department of Agriculture (USDA). All opinions expressed in this paper are the author’s and do not necessarily reflect the policies and views of the USDA-ARS FADRU, Vietnam DAH or ORAU/ORISE. The Vietnam DAH and USDA-ARS FADRU are members of the Global Foot-and-Mouth Disease Research Alliance (GFRA). The work contributed by the WRLFMD was supported by the Department for Environment, Food and Rural Affairs (Project SE2943: Defra, UK), and funding provided to the EuFMD from the European Union.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jonathan Arzt.

Additional file

13567_2017_424_MOESM1_ESM.docx

Additional file 1. Foot-and-mouth disease viruses from Vietnam, Kazakhstan and China included in the study analyses. Genbank accession numbers, virus name, geographical location, collection date, host species and clinical stage (outbreak or persistently infected animals) are tabulated.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Brito, B., Pauszek, S.J., Eschbaumer, M. et al. Phylodynamics of foot-and-mouth disease virus O/PanAsia in Vietnam 2010–2014. Vet Res 48, 24 (2017). https://0-doi-org.brum.beds.ac.uk/10.1186/s13567-017-0424-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13567-017-0424-7

Keywords