Low linkage disequilibrium indicative of recombination in foot-and-mouth disease virus gene sequence alignments

Demonstration of Co-Infection and Trans-Encapsidation of Viral RNA In Vitro Using Epitope-Tagged Foot-and-Mouth Disease Viruses

Foot-and-mouth disease, caused by foot-and-mouth disease virus (FMDV), is an economically devastating disease affecting several important livestock species. FMDV is antigenically diverse and exists as seven serotypes comprised of many strains which are poorly cross-neutralised by antibodies induced by infection or vaccination. Co-infection and recombination are important drivers of antigenic diversity, especially in regions where several serotypes co-circulate at high prevalence, and therefore experimental systems to study these events in vitro would be beneficial. Here we have utilised recombinant FMDVs containing an HA or a FLAG epitope tag within the VP1 capsid protein to investigate the products of co-infection in vitro. Co-infection with viruses from the same and from different serotypes was demonstrated by immunofluorescence microscopy and flow cytometry using anti-tag antibodies. FLAG-tagged VP1 and HA-tagged VP1 could be co-immunoprecipitated from co-infected cells, suggesting that newly synthesised capsids may contain VP1 proteins from both co-infecting viruses. Furthermore, we provide the first demonstration of trans-encapsidation of an FMDV genome into capsids comprised of proteins encoded by a co-infecting heterologous virus. This system provides a useful tool for investigating co-infection dynamics in vitro, particularly between closely related strains, and has the advantage that it does not depend upon the availability of strain-specific FMDV antibodies.

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.

Co-evolution of sites under immune selection shapes Epstein-Barr Virus population structure

Epstein–Barr virus (EBV) is one of the most common viral infections in humans and persists within its host for life. EBV therefore represents an extremely successful virus that has evolved complex strategies to evade the host’s innate and adaptive immune response during both initial and persistent stages of infection. Here, we conducted a comparative genomics analysis on 223 whole genome sequences of worldwide EBV strains. We recover extensive genome-wide linkage disequilibrium (LD) despite pervasive genetic recombination. This pattern is explained by the global EBV population being subdivided into three main subpopulations, one primarily found in East Asia, one in Southeast Asia and Oceania, and the third including most of the other globally distributed genomes we analyzed. Additionally, sites in LD were overrepresented in immunogenic genes. Taken together, our results suggest that host immune selection and local adaptation to different human host populations has shaped the genome-wide patterns of genetic diversity in EBV.

Within-Host Recombination in the Foot-and-Mouth Disease Virus Genome

Recombination is one of the determinants of genetic diversity in the foot-and-mouth disease virus (FMDV). FMDV sequences have a mosaic structure caused by extensive intra- and inter-serotype recombination, with the exception of the capsid-encoding region. While these genome-wide patterns of broad-scale recombination are well studied, not much is known about the patterns of recombination that may exist within infected hosts. In addition, detection of recombination among viruses evolving at the within-host level is challenging due to the similarity of the sequences and the limitations in differentiating recombination from point mutations. Here, we present the first analysis of recombination events between closely related FMDV sequences occurring within buffalo hosts. The detection of these events was made possible by the occurrence of co-infection of two viral swarms with about 1% nucleotide divergence. We found more than 15 recombination events, unequally distributed across eight samples from different animals. The distribution of these events along the FMDV genome was neither uniform nor related to the phylogenetic distribution of recombination breakpoints, suggesting a mismatch between within-host evolutionary pressures and long-term selection for infectivity and transmissibility.

Phylogeographic analysis of the 2000-2002 foot-and-mouth disease epidemic in Argentina

Foot-and-mouth disease (FMD) is a highly transmissible disease of hooved livestock. Although FMD has been eradicated from many countries, economic and social consequences of FMD reintroductions are devastating. After achieving disease eradication, Argentina was affected by a major epidemic in 2000-2002, and within few months, FMD virus spread throughout most of the country and affected >2500 herds. Available records and viral strains allowed us to assess the origins, spread and progression of this FMD epidemic, which remained uncertain. We used whole genome viral sequences and a continuous phylogeographic diffusion approach, which revealed that the viruses that caused the outbreaks spread fast in different directions from a central area in Argentina. The analysis also suggests that the virus that caused the outbreaks in the year 2000 was different from those found during the 2001 epidemic. To estimate if the approximate overall genetic diversity of the virus was related to disease transmission, we reconstructed the viral demographic variation in time using Bayesian Skygrid approach and compared it with the epidemic curve and the within-herd transmission rate and showed that the genetic temporal diversity of the virus was associated with the increasing number of outbreaks in the exponential phase of the epidemic. Results here provide new evidence of how the disease entered and spread throughout the country. We further demonstrate that genetic data collected during a FMD epidemic can be informative indicators of the progression of an ongoing epidemic.

The impact of within-herd genetic variation upon inferred transmission trees for foot-and-mouth disease virus

•

45 full FMDV genomes were analysed from farms infected in the UK during 2007.

•

Analyses revealed intra-herd clustering consistent with virus transmission events.

•

Inter-herd sequence differences were estimated to be 4.6 nucleotides.

•

Random selection of 1 sequence from each herd generated similar transmission trees.

•

These results will help to design cost-effective approaches for the control of FMD.

Full-genome sequences have been used to monitor the fine-scale dynamics of epidemics caused by RNA viruses. However, the ability of this approach to confidently reconstruct transmission trees is limited by the knowledge of the genetic diversity of viruses that exist within different epidemiological units. In order to address this question, this study investigated the variability of 45 foot-and-mouth disease virus (FMDV) genome sequences (from 33 animals) that were collected during 2007 from eight premises (10 different herds) in the United Kingdom. Bayesian and statistical parsimony analysis demonstrated that these sequences exhibited clustering which was consistent with a transmission scenario describing herd-to-herd spread of the virus. As an alternative to analysing all of the available samples in future epidemics, the impact of randomly selecting one sequence from each of these herds was used to assess cost-effective methods that might be used to infer transmission trees during FMD outbreaks. Using these approaches, 85% and 91% of the resulting topologies were either identical or differed by only one edge from a reference tree comprising all of the sequences generated within the outbreak. The sequence distances that accrued during sequential transmission events between epidemiological units was estimated to be 4.6 nucleotides, although the genetic variability between viruses recovered from chronic carrier animals was higher than between viruses from animals with acute-stage infection: an observation which poses challenges for the use of simple approaches to infer transmission trees. This study helps to develop strategies for sampling during FMD outbreaks, and provides data that will guide the development of further models to support control policies in the event of virus incursions into FMD free countries.

Recombination and evolutionary dynamics of human echovirus 6

Enterovirus (EV) infections are associated with a wide array of often severe disease presentations including aseptic meningitis, encephalitis, and acute flaccid paralysis. Surveillance for polioviruses and other EVs is therefore important as a public health measure both for patient management and epidemiological studies. From 1988 to 2008, echovirus (E) 30 was the predominant genotype in Spain (33.7% of the total typed EVs). E6 was also endemic throughout this period although isolated less frequently (12.5%). In 2009, however, a substantial increase in the incidence of E6 was detected (60%), displacing E30 type (2%). To investigate the evolution and recombination in the epidemiology and transmission of E6 in Spain, a genetic analysis in VP1 and 3Dpol regions of 67 Spanish strains collected during the period 2004-2010 was performed. All VP1 sequences clustered monophyletically in the assigned genogroup C, subgroup 9, currently the predominant circulating strains identified in Europe and elsewhere in the last 10 years. 3Dpol sequences were interspersed with other species B EVs resulting from several recombination events that generated at least 12 different recombinant forms (RFs) among study samples. These showed typically minimal divergence in VP1. The co-circulation of different lineages of E6 in the same geographical area associated with its mainly endemic pattern of transmission may have contributed to the extremely short estimated half-life of E6 RFs (0.87 years). This pattern contrasts markedly with other species B EVs and EV71 where VP1 lineage expansion and extinction occurred in step with defined recombination events and periodic changes in incidence.

Isolation, identification and complete genome sequence analysis of a strain of foot-and-mouth disease virus serotype Asia1 from pigs in southwest of China

Backgroud

Foot-and-mouth disease virus (FMDV) serotype Asia1 generally infects cattle and sheep, while its infection of pigs is rarely reported. In 2005-2007, FMD outbreaks caused by Asia1 type occurred in many regions of China, as well as some parts of East Asia countries. During the outbreaks, there was not any report that pigs were found to be clinically infected.

Results

In this study, a strain of FMDV that isolated from pigs was identified as serotype Asia1, and designated as "Asia1/WHN/CHA/06". To investigate the genomic feature of the strain, complete genome of Asia1/WHN/CHA/06 was sequenced and compared with sequences of other FMDVs by phylogenetic and recombination analysis. The complete genome of Asia1/WHN/CHA/06 was 8161 nucleotides (nt) in length, and was closer to JS/CHA/05 than to all other strains. Potential recombination events associated with Asia1/WHN/CHA/06 were found between JS/CHA/05 and HNK/CHA/05 strains with partial 3B and 3C fragments.

Conclusion

This is the first report of the isolation and identification of a strain of FMDV type Asia1 from naturally infected pigs. The Asia1/WHN/CHA/06 strain may evolve from the recombination of JS/CHA/05 and HNK/CHA/05 strains.

Loss of linkage disequilibrium and accelerated protein divergence in duplicated cytomegalovirus chemokine genes

Human CMV (hCMV) encodes several captured chemokine ligand and chemokine receptor genes that may play a role in immune evasion. The adjacent viral alpha-chemokine genes UL146 and UL147 appear to have duplicated subsequent to a recent gene capture event. Sequence data from multiple hCMV isolates suggest accelerated protein evolution in one of the paralogues, UL146. Extensive sequence variation was noted throughout the more rapidly evolving paralogue, although significant variation was also observed within the more slowly evolving gene, especially within a region corresponding to a possible signal peptide. In contrast to the haplotype structure observed for other hCMV genes, the distribution of nucleotide variants indicates a marked loss of linkage disequilibrium within UL146 and to a lesser extent UL147. Despite evidence of accelerated protein evolution, the rate of nonsynonymous to synonymous substitutions (d(N)/d(S)) in the more rapidly evolving paralogue was not indicative of neutral evolution, but of moderate purifying selection. The data presented here provides a unique opportunity to study the mechanisms by which a recently duplicated pair of genes has diverged and suggests a role for recombination.

Understanding the molecular epidemiology of foot-and-mouth-disease virus

The use of molecular epidemiology is an important tool in understanding and consequently controlling FMDV. In this review I will present basic information about the disease, needed to perform molecular epidemiology. I will give a short introduction to the history and impact of foot-and-mouth disease, clinical picture, infection route, subclinical and persistent infections, general aspects of the transmission of FMDV, serotype-specific epidemiological characteristics, field epidemiology of FMDV, evolution and molecular epidemiology of FMDV. This is followed by two chapters describing the molecular epidemiology of foot-and-mouth disease in global surveillance and molecular epidemiology of foot-and-mouth disease in outbreak investigation.

Multiplexed molecular assay for rapid exclusion of foot-and-mouth disease

A nucleic acid-based multiplexed assay was developed that combines detection of foot-and-mouth disease virus (FMDV) with rule-out assays for two other foreign animal diseases and four domestic animal diseases that cause vesicular or ulcerative lesions indistinguishable from FMDV infection in cattle, sheep and swine. The FMDV "look-alike" diagnostic assay panel contains 5 PCR and 12 reverse transcriptase PCR (RT-PCR) signatures for a total of 17 simultaneous PCR amplifications for 7 diseases plus incorporating 4 internal assay controls. It was developed and optimized to amplify both DNA and RNA viruses simultaneously in a single tube and employs Luminex liquid array technology. Assay development including selection of appropriate controls, a comparison of signature performance in single and multiplex testing against target nucleic acids, as well of limits of detection for each of the individual signatures is presented. While this assay is a prototype and by no means a comprehensive test for FMDV "look-alike" viruses, an assay of this type is envisioned to have benefit to a laboratory network in routine surveillance and possibly for post-outbreak proof of freedom from foot-and-mouth disease.

Mosaic structure of foot-and-mouth disease virus genomes

The results of a simple pairwise-scanning analysis designed to identify inter-serotype recombination fragments, applied to genome data from 156 isolates of Foot-and-mouth disease virus (FMDV) representing all seven serotypes, are reported. Large numbers of candidate recombinant fragments were identified from all parts of the FMDV genome, with the exception of the capsid genes, within which such fragments are infrequent. As expected, intertypic fragment exchange is most common between geographically sympatric FMDV serotypes. After accounting for the likelihood of intertypic convergence in highly conserved parts of the FMDV genome, it is concluded that intertypic recombination is probably widespread throughout the non-structural genes, but that recombination over the 2B/C and 3B/C gene boundaries appears to be less frequent than expected, given the large numbers of recombinant gene fragments arising in these genes.

Recombination patterns in aphthoviruses mirror those found in other picornaviruses

Foot-and-mouth disease virus (FMDV) is thought to evolve largely through genetic drift driven by the inherently error-prone nature of its RNA polymerase. There is, however, increasing evidence that recombination is an important mechanism in the evolution of these and other related picornoviruses. Here, we use an extensive set of recombination detection methods to identify 86 unique potential recombination events among 125 publicly available FMDV complete genome sequences. The large number of events detected between members of different serotypes suggests that horizontal flow of sequences among the serotypes is relatively common and does not incur severe fitness costs. Interestingly, the distribution of recombination breakpoints was found to be largely nonrandom. Whereas there are clear breakpoint cold spots within the structural genes, two statistically significant hot spots precisely separate these from the nonstructural genes. Very similar breakpoint distributions were found for other picornovirus species in the genera Enterovirus and Teschovirus. Our results suggest that genome regions encoding the structural proteins of both FMDV and other picornaviruses are functionally interchangeable modules, supporting recent proposals that the structural and nonstructural coding regions of the picornaviruses are evolving largely independently of one another.

Recombination and selection in the evolution of picornaviruses and other Mammalian positive-stranded RNA viruses

Picornaviridae are a large virus family causing widespread, often pathogenic infections in humans and other mammals. Picornaviruses are genetically and antigenically highly diverse, with evidence for complex evolutionary histories in which recombination plays a major part. To investigate the nature of recombination and selection processes underlying the evolution of serotypes within different picornavirus genera, large-scale analysis of recombination frequencies and sites, segregation by serotype within each genus, and sequence selection and composition was performed, and results were compared with those for other nonenveloped positive-stranded viruses (astroviruses and human noroviruses) and with flavivirus and alphavirus control groups. Enteroviruses, aphthoviruses, and teschoviruses showed phylogenetic segregation by serotype only in the structural region; lack of segregation elsewhere was attributable to extensive interserotype recombination. Nonsegregating viruses also showed several characteristic sequence divergence and composition differences between genome regions that were absent from segregating virus control groups, such as much greater amino acid sequence divergence in the structural region, markedly elevated ratios of nonsynonymous-to-synonymous substitutions, and differences in codon usage. These properties were shared with other picornavirus genera, such as the parechoviruses and erboviruses. The nonenveloped astroviruses and noroviruses similarly showed high frequencies of recombination, evidence for positive selection, and differential codon use in the capsid region, implying similar underlying evolutionary mechanisms and pressures driving serotype differentiation. This process was distinct from more-recent sequence evolution generating diversity within picornavirus serotypes, in which neutral or purifying selection was prominent. Overall, this study identifies common themes in the diversification process generating picornavirus serotypes that contribute to understanding of their evolution and pathogenicity.

A simple and robust statistical test for detecting the presence of recombination

Recombination is a powerful evolutionary force that merges historically distinct genotypes. But the extent of recombination within many organisms is unknown, and even determining its presence within a set of homologous sequences is a difficult question. Here we develop a new statistic, phi(w), that can be used to test for recombination. We show through simulation that our test can discriminate effectively between the presence and absence of recombination, even in diverse situations such as exponential growth (star-like topologies) and patterns of substitution rate correlation. A number of other tests, Max chi2, NSS, a coalescent-based likelihood permutation test (from LDHat), and correlation of linkage disequilibrium (both r2 and /D'/) with distance, all tend to underestimate the presence of recombination under strong population growth. Moreover, both Max chi2 and NSS falsely infer the presence of recombination under a simple model of mutation rate correlation. Results on empirical data show that our test can be used to detect recombination between closely as well as distantly related samples, regardless of the suspected rate of recombination. The results suggest that phi(w) is one of the best approaches to distinguish recurrent mutation from recombination in a wide variety of circumstances.

Comparative genomics of foot-and-mouth disease virus

Here we present complete genome sequences, including a comparative analysis, of 103 isolates of foot-and-mouth disease virus (FMDV) representing all seven serotypes and including the first complete sequences of the SAT1 and SAT3 genomes. The data reveal novel highly conserved genomic regions, indicating functional constraints for variability as well as novel viral genomic motifs with likely biological relevance. Previously undescribed invariant motifs were identified in the 5' and 3' untranslated regions (UTR), as was tolerance for insertions/deletions in the 5' UTR. Fifty-eight percent of the amino acids encoded by FMDV isolates are invariant, suggesting that these residues are critical for virus biology. Novel, conserved sequence motifs with likely functional significance were identified within proteins L(pro), 1B, 1D, and 3C. An analysis of the complete FMDV genomes indicated phylogenetic incongruities between different genomic regions which were suggestive of interserotypic recombination. Additionally, a novel SAT virus lineage containing nonstructural protein-encoding regions distinct from other SAT and Euroasiatic lineages was identified. Insights into viral RNA sequence conservation and variability and genetic diversity in nature will likely impact our understanding of FMDV infections, host range, and transmission.

Evolutionary timescale of rabies virus adaptation to North American bats inferred from the substitution rate of the nucleoprotein gene

Throughout North America, rabies virus (RV) is endemic in bats. Distinct RV variants exist that are closely associated with infection of individual host species, such that there is little or no sustained spillover infection away from the primary host. Using Bayesian methodology, nucleotide substitution rates were estimated from alignments of partial nucleoprotein (N) gene sequences of nine distinct bat RV variants from North America. Substitution rates ranged from 2.32 x 10(-4) to 1.38 x 10(-3) substitutions per site per year. A maximum-likelihood (ML) molecular clock model was rejected for only two of the nine datasets. In addition, using sequences from bat RV variants across the Americas, the evolutionary rate for the complete N gene was estimated to be 2.32 x 10(-4). This rate was used to scale trees using Bayesian and ML methods, and the time of the most recent common ancestor for current bat RV variant diversity in the Americas was estimated to be 1660 (range 1267-1782) and 1651 (range 1254-1773), respectively. Our reconstructions suggest that RV variants currently associated with infection of bats from Latin America (Desmodus and Tadarida) share the earliest common ancestor with the progenitor RV. In addition, from the ML tree, times were estimated for the emergence of the three major lineages responsible for bat rabies cases in North America. Adaptation to infection of the colonial bat species analysed (Eptesicus fuscus, Myotis spp.) appears to have occurred much quicker than for the solitary species analysed (Lasionycteris noctivagans, Pipistrellus subflavus, Lasiurus borealis, Lasiurus cinereus), suggesting that the process of virus adaptation may be dependent on host biology.