key: cord-0764156-t8grxsxb
authors: Steinbach, F.; Westcott, D. G.; McGowan, S. L.; Grierson, S. S.; Frossard, J. P.; Choudhury, B.
title: Re-emergence of a genetic outlier strain of equine arteritis virus: Impact on phylogeny
date: 2015-04-16
journal: Virus Research
DOI: 10.1016/j.virusres.2014.12.009
sha: a29da3ff89510405fc9ebc914f6614236e439fbc
doc_id: 764156
cord_uid: t8grxsxb

Abstract Equine arteritis virus (EAV) is the causative agent of equine viral arteritis (EVA), a respiratory and reproductive disease of equids, which is notifiable in some countries including the Great Britain (GB) and to the OIE. Herein, we present the case of a persistently infected stallion and the phylogenetic tracing of the virus strain isolated. Discussing EAV occurrence and phylogenetic analysis we review features, which may aid to harmonise and enhance the classification of EAV.

EAV is an enveloped, linear, positive sense, single stranded RNA virus representing the prototype virus of the Arteriviridae family, in the order Nidovirales. This family also includes simian haemorrhagic fever virus, porcine reproductive and respiratory syndrome virus (PRRSV) and lactate dehydrogenase-elevating virus (Cavanagh, 1997) .

The genome of EAV ( Fig. 1) is approximately 12.7 kb in length and consists of 10 overlapping open reading frames (ORFs), which are flanked by 5 and 3 untranslated regions (UTRs) (Snijder and Meulenberg, 1998; Firth et al., 2011; Balasuriya et al., 2014) . ORFs 1a and b occupy approximately 75% of the genome and both are translated to produce polyproteins, which are then extensively processed into 13 non-structural proteins (Table 1) (Snijder and Meulenberg, 1998) . At the 3 end, occupying approximately 25% of the genome, the remaining eight ORFs, 2-7 are nested and encode seven envelope proteins and the viral nucleocapsid (Table 1) (Snijder and Meulenberg, 1998; Firth et al., 2011) .

Although there is only one serotype of EAV, variation exists between field strains of the virus and their neutralising properties, i.e. neutralising epitopes are not shared across all strains as evident from the use of monoclonal antibodies (Glaser et al., 1995; Balasuriya et al., 1997) . ORF5 encodes the major envelope protein, GP5, which contains the best characterised neutralisation determinants. This is also a most variable region within the genome and hence is commonly used for phylogenetic analysis (Balasuriya et al., 1995; Stadejek et al., 1999; Hornyak et al., 2005; Mittelholzer et al., 2006; Metz et al., 2011) . Phylogenetic analysis has also been conducted using ORFs 3, 6 and 7 (Lepage et al., 1996; Hedges et al., 2001; Chirnside et al., 1994) .

Despite variations in gene, phylogenetic analysis models and number of sequences used the general consensus of such phylogenetic studies show a segregation of sequences into European and North American lineages of which the European lineage can be further divided into two clusters. Unsurprisingly, some European lineage viruses have been documented in North America and vice versa, highlighting the transport of viruses between the two continents via trade of horses and semen.

Excluding patent applications and molecular clones there are 800 plus EAV sequences available in GenBank (as of June 2014). Thirty-two of these are full genomes: two (accession numbers: NC002532 and DQ846750) are different iterations of Bucyrus, sharing 99.6% homology at the nucleotide level; five are from the French outbreak in 2007 and the other 25 are from the USA.

Of the remaining sequences, around 25% of sequences deposited are of North American origin, another 25% were reported in Western Europe and a further 25% from Eastern Europe. There is only one Asian sequence, from China (JX868590) and three from the Middle East (Qatar; accession no. AY453329-32). The African continent is also poorly represented with only 16 sequences, all reported from South Africa. There are around sixty sequences for which the location was not stated in the GenBank record or associated manuscripts.

As mentioned the majority of phylogenetic analysis so far has been conducted on ORF5. Of the other segments used for phylogenetic analysis, ORF3 offers over 250 sequences, ORF6 over 210 sequences and ORF7 over 280 sequences, compared to nearly 600 Snijder and Kikkert, 2013 and GenBank's genome section on EAV). The ORF1ab (blue) is expressed as a fusion protein, from the genomic RNA, despite the −1 programmed ribosomal frameshift required for ORF1b (light blue). The genes encoding the structural proteins (red) are translated from subgenomic RNAs. Polyprotein cleavage sites are depicted above ORF1ab. Blue arrowheads: sites cleaved by PLP domains of nsp1/nsp2; red arrowheads: sites cleaved by the nsp4. The graph displays the relative sizes and locations of the resulting proteins. Proteins expressed in frame to each other are displayed on the same axis. NSP: non-structural protein; GP: glycoprotein; M: membrane protein; E: envelope protein; N: nucleocapsid protein. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) sequences which will allow use of ORF5 e.g. including those which cover ORFs 2-7. However, where the focus is singularly on ORF5 only, the majority (188/278) sequences deposited in GenBank are partial sequences which are 518 nucleotides in length and cover genome positions 11296-11813 (Bucyrus) representing the most variable region within EAV ORF5.

EAV was first isolated in Bucyrus, OH, USA in 1953, from foetal lung tissue during an outbreak of respiratory disease and abortions in horses (Doll et al., 1957; reviewed in Balasuriya et al., 2014) . The virus is thought to have an almost global distribution: via serological surveys EAV has been noted in equine populations in Europe, Africa, Asia, the Americas and Australia (Timoney and McCollum, 1993) . The lack of unified control measures, poor biosecurity and the increased international trade in horses has assisted in the spread of EAV. The only countries considered to be disease free are Iceland and Japan. Recently, McFadden et al. (2013) demonstrated a lack of EAV prevalence in the general horse population in New Zealand.

A hallmark of EAV is its persistence in and spread through stallions. Larger clinical outbreaks occur only periodically and genetic markers for pathogenicity have been suggested, but lack further characterisation (Balasuriya et al., 2014) . The first clinical appearance of EAV in the United Kingdom was in 1993 (Wood et al., 1995) . The source of the outbreak was a stallion imported from Poland. The outbreak began on a non-thoroughbred stud and spread to five other premises via chilled semen used for artificial insemination. Around 100 animals became infected and the outbreak was contained by means of voluntary movement restrictions. Subsequently EAV was made notifiable in the UK (EAV order 1995). Since then, there have been sporadic cases up to and including the case of the imported Spanish stallion in 2012, which is described in this manuscript.

During the last decade clinical outbreaks of EAV occurred in several countries on both sides of the Atlantic. In the summer of 2007, the largest EVA outbreak in France occurred causing mortality and economic loss due to disruption of race schedules (Pronost et al., 2010) . The outbreak affected 18 premises in 5 counties in western France, which represented the index, 8 primary and 9 secondary premises. Artificial insemination in draught horses was responsible for the index case.

In a similar timeframe (2006-2007), a multistate US outbreak occurred and for the first time disease was observed in Quarter Horses (Powell and Timoney, 2006) . Due to the nature of EAV dissemination through infectious semen, EVA prevalence can show disparity between horse breeds e.g. prior to that outbreak seroprevalance in Quarter Horses was recorded at 0.6% (Anonymous, 2000) .

There is limited information regarding EAV prevalence in South America. Echeverría et al. (2003) reported the first isolation of EAV in Argentina, and perhaps South America, in 2001. The virus was isolated from the semen of an imported seropositive stallion held in isolation at a breeding farm in the Buenos Aires Province.

Three additional isolates were obtained in 2007 with phylogenetic analysis suggesting an European lineage (Echeverría et al., 2007) . Shortly thereafter in 2010 another significant outbreak occurred in Argentina, which was linked to a semen import from the Netherlands (Barrandeguy, 2010) . Table 1 EAV genes, proteins and their functions.

Protein Function 1a nsp1 Accessory protease: important role in virion biogenesis and crucial in processing of the replicase polyprotein and the production of subgenomic RNA (Tijms et al., 2007) nsp2 Accessory-protease: auto-protease activity between nsp2 and 3. Co-factor for the successful processing of nsp4/5 (Snijder et al., 1994) . Deubiquitinating enzyme activity (van Kasteren et al., 2013) . nsp3

Transmembrane non-structural protein thought to be involved in double-membrane vesicle formation (Posthuma et al., 2008) . nsp4

Main protease: responsible for cleaving nsp3-8 and nsp3-12 polypeptides (Snijder et al., 1994) . nsp5

Suggested role in the membrane association of the EAV replication complex. (Snijder et al., 1994) . nsp6

Function not determined. nsp7

Further cleaved into nsp7␣ and nsp7␤, function not fully determined, critical role in RNA synthesis (Manolaridis et al., 2011) . nsp8

Function not determined. 1b nsp9 RNA-dependent RNA polymerase (Snijder et al., 1994) . nsp10

Helicase (van Dinten et al., 1997) . nsp11

Multifunctional including endoribonuclease activity, critical role in RNA synthesis (Posthuma et al., 2006) . nsp12

Function not determined. 2a E Predicted to be a membrane protein which is essential for the production of infectious progeny (Snijder et al., 1999) . 2b GP2 Minor envelope protein: membrane protein which is essential for the production of infectious progeny (Snijder et al., 1999) . 3 GP3 Membrane associated protein: function not fully determined but thought essential for viral replication (Wieringa et al., 2002) . 4 GP4 Membrane glycoprotein: function not fully determined but thought essential for viral replication (Wieringa et al., 2002) . 5a

ORF5a

Membrane protein: function not fully determined but thought important for viral replication (Firth et al., 2011 (Table 2) with GB Glos 2012. The Maximumlikelihood phylogeny was constructed using the GTR +G+ I model. Sequences are listed by their accession numbers followed by strain and country of isolation. Percentage of a 1000 bootstrap replicates are cited at nodes for values ≥50%.

Scant information is only available regarding the disease situation in Africa. Paweska et al. (1996) reported high seroprevalence of an asinine strain of EAV in South African donkeys. Guthrie et al. (2003) reported lateral transmission among Lipizzaner stallions caused by an imported stallion from Yugoslavia in 1981. A few years later, Stadejek et al. (2006) reported the isolation of a highly divergent EAV strain from a South African donkey.

No internationally published information appears to be available for remainder of Africa, Asia and the Middle East; however the fact there are sequences from these regions available on GenBank shows the disease to be prevalent there.

In September 2012 we analysed a semen sample from a 9-yearold Andalusian stallion which was imported to Gloucestershire, UK from Spain. Virus was isolated following the protocol described in chapter 2.5.10 of the OIE Terrestrial Manual 2013 (Anonymous, 2013) . The resultant isolate was named GB Glos 2012, reflecting the country, province and year of isolation.

RNA was extracted from the semen using the QIAamp Viral RNA Mini Kit (Qiagen Ltd., Manchester, UK) according to the manufacturer's instructions. A 564 nucleotide portion of the ORF5 gene was amplified using the CR2 and EAV32 primers described by Stadejek et al. (1999) . In parallel, a 395 nucleotide portion of the ORF7 gene was amplified using the 14a and 15 primers described by Belak et al. (1994) . Following purification, the PCR products were Sanger sequenced using the same primer pair on ABI PRISM 3130xl Genetic Analyser machine (Applied Biosystems Inc.).

For phylogenetic analysis sequences which represent all described EAV genetic groups were downloaded from GenBank and compared to the sequence (GB Glos 2012) generated in this study. The sequences were aligned and maximum likelihood and neighbour joining phylogenies constructed using the MEGA5 programme (Tamura et al., 2011) . Table 2 The percentage sequence homology between GB Glos 2012 and the other selected full genomes.

The GB Glos 2012 isolate was also used as the starting material for next generation sequencing via the Roche 454 FLX platform following protocols described by Rasmussen et al. (2008 Rasmussen et al. ( , 2010 .

The semen sample tested positive for EAV both by virus isolation and PCRs targeting ORF5 and ORF7.

It was possible to recover the full genome of GB Glos 2012 (Gen-Bank: LC000003) via NGS. All of the full genome sequences available on GenBank (n = 32, excluding molecular clones) were compared phylogenetically. As expected comparing these genomes did not elucidate any details regarding origins as the limited geographic distribution of the available data hinders any meaningful discovery. However, GB Glos 2012 fell outside all existing full genome sequences and by using the existing clusters it was possible to take a reference sequence from each to compare genome features against GB Glos 2012 (Fig. 2) . The GB Glos 2012 isolate is 12,702 nucleotides in length, slightly shorter than Bucyrus (12,704 nucleotides), CW01 (12,708 nucleotides), F27 (12,710 nucleotides) and S4445 (12,731 nucleotides). The isolate shares greatest homology, 85.5% at the nucleotide level with Bucyrus and 84.4%, 84.7% and 84.6% and homology with CW01, F27 and S4445 respectively ( Table 2) . Analysing homology at the amino acid level for specific ORFs least similarity was noted across ORF3 at 78.5% with F27 and greatest similarity was across ORF7 at 99.1% with Bucyrus, CW01 and S4445. GB Glos 2012 had a three nucleotide deletion in nsp2 (codon 459 of protein 1ab) apart from which no other indels were noted.

After a thorough analysis of GenBank entries, phylogenies were initially constructed using 211 sequences. Sequences were selected to cover greatest geographic and time differences, thus molecular clones, Bucyrus derivatives and laboratory strains were not selected for the analysis. Reference sequences from each of the manuscripts cited in Fig. 3 were used to inform clusters. Phylogenies were constructed in both neighbour-joining and maximum likelihood methods, which made little difference to topology (data not shown). For ease of viewing Fig. 3 contains an abridged version of the phylogenetic tree containing 101 sequences on display.

Analysis of the GB Glos 2012 ORF5 PCR product sequence (Fig. 3 ) revealed the isolate to be an outlier, only clustering with isolates KY63 (GenBank: U81014), H10 (GenBank: AY453306) and H21 (GenBank: AY453306) supported by bootstrap values of 70% and 99% respectively. GB Glos 2012 shared 87%, 94.6% and 95% nucleotide homology across ORF5 GL segment with KY63, H10 and H21 respectively. The next closest neighbours were strains isolated in South Africa (GenBank: AY453337 and AY4533379) which were thought to be of Yugoslavian origin and Hungarian strains (Gen-Bank: JN314862 and JN314866).

Following tree construction 7 clusters alphabetically labelled A-G, became apparent (Fig. 3) . In our hands three new clusters A, B and G were visible, which were supported by bootstrap values of 96%, 70% and 100% respectively. Clusters A and B have previously been included in group EU1 and cluster G within group EU2 (Zhang et al., 2007) . However taking the bootstrap values and the within between group variation (shown in Table 3 ) into account separate clusters appear warranted.

There were three clusters, C, D, and F in which the majority of isolates segregated. The existence of one cluster, C, however, was not supported via bootstrap analysis, and regardless of the phylogenetic scheme used it has the greatest within group variation with 9.9%/10.10%/9.9% and 13.41% respectively (Table 3) . This, along with the lack of bootstrap support suggests this cluster could be further broken down. Comparing the previously published phylogenies the IA and IB groupings cited in Stadejek et al. (1999) were not sustainable with current information and thus merged into one group as already suggested by Mittelholzer et al. (2006) and Zhang et al. (2007) , here labelled group F in Fig. 3 . Regardless of the phylogenetic scheme used there are isolates which did not fit into any of the defined clusters (Fig. 3) : the use of full genome data and/or further sampling may aid to resolve this, either by generating further clusters or including the viruses into existing ones.

Reviewing the relation of GB Glos 2012 with other EAV strains, we focused our phylogeny on ORF5, as did previous authors (Balasuriya et al., 1995; Stadejek et al., 1999; Hornyak et al., 2005; Mittelholzer et al., 2006; Metz et al., 2011) . Whilst it would have been interesting to compare further gene segments appropriate geographic and time point distribution were observed across ORF5 sequences only. As seen from Fig. 2 limited information could be gleaned from phylogeny of full genome data as these data are limited in number and location. It is anticipated that with the development of full genome sequencing into a routine technology, more viruses from all currently existing clades will be added soon to allow for a wider comparison. With regards to GB Glos 2012 it would be particularly interesting to compare it against the South African and Hungarian strains with which the isolate clustered as here we were limited to the 518 nucleotides of the ORF5 available for the other strains. As with PRRSV, variability of nsp2 is also noted for EAV. We reported a three nucleotide deletion at codon 459 of protein 1ab and comparison with other full genome sequences highlighted numerous alterations in this region e.g. strain S4445 contains a 5 amino acid insertion between codons 448 and 489 of protein 1ab. As with the phylogeny, we were again limited in our ability to collate such observations as the publicly available sequence data for the nsp2 region is restricted to data from the USA or France with the majority of these sequences being linked to a single outbreak in each of these countries.

The re-appearance of a rare genotype of EAV provided us with reason to re-visit the relationship of EAV strains and their nomenclature. GB Glos 2012 demonstrates the re-emergence of a genotype across decades and around the globe. KY63 is a field strain isolated in Kentucky, USA in 1963; this isolate has always clustered separately from the other US isolates of the same time period suggesting a different origin. Strains H10 and H21 appear to have been isolated in the early 2000s in Hungary and the "South African" strains were isolated in the mid-1990s with a history of Eastern European imports, whereas the isolated characterised here (GB Glos 2012) is from a Spanish stallion. Even considering a potential lack of sampling the case demonstrates that established Fig. 3 . Phylogenetic analysis of ORF5. Maximum-likelihood tree using GTR +G+ I model with 101 EAV sequences shown. Phylogenetic grouping as suggested by Stadejek et al. (1999 and Zhang et al. (2007) are labelled on the right hand side of the tree. Phylogenetic scheme suggested from this work is displayed on the left of the tree. Sequences are listed by their accession numbers followed by strain name, country of isolation and where available year of isolation. Percentage of a 1000 bootstrap replicates are cited at nodes for values ≥50%.

The average percentage pairwise nucleotide distances within clades. Stadejek et al. (1999 ) Zhang et al. (2007 strains do not become extinct and through global trade may reappear unexpectedly. This is surprising since Arteriviruses, such as EAV or PRRSV are known for their rapid evolution, linked to the lack of a proofreading enzyme such as the exoribonuclease ExoN, which Coronaviruses for example contain and which is vital for their genetic integrity (Minskaia et al., 2006; Nga et al., 2011) . With viruses being subjected to host pressure the general notion was that they are constantly changing, perhaps as is reflected in the phylogenetic trees for EAV and PRRSV getting ever more complicated and seemingly highlighting new clades and genotypes. However, this case as well as the recent analysis of PRRSV and Bovine viral diarrhoea virus (BVDV; another persistent RNA virus of ruminants) in the United Kingdom demonstrates the persistence of strain variations in animal populations over decades (Frossard et al., 2013; Strong et al., 2013) . This significantly broadens our understanding of the evolution of such viruses: in the absence of appropriate control measure to eradicate certain viruses (strains), host pressure will lead to a continuous divergence of viruses, but with the parallel persistence of "old" variants.

If phylogenies are to be used for epidemiological tracing further data points to fill in sampling gaps need to be included. That will require a more active sampling and control policy than currently established in most countries. The fact that this case (as some recent cases) was detected in a horse of some economic value, which was supposed to be used for breeding purposes, emphasizes the need for further vigilance, in particular in relation to the health horse concept advocated by the OIE.

Additionally, while ORF5 represents a good starting point with its gene product GP5 being directly under host-pressure from the antibody response other genes are important for the pathogenicity too and recombination has been described for Arteriviruses, thus full genome data certainly represents the state of the art for the description of index cases in outbreaks. The available data needs to be labelled appropriately with the sequence record containing the location and date of isolation: as mentioned above there are multiple sequences at present in GenBank without any such information. Similarly manuscripts citing sequences should state accession numbers and particularly if presenting novel sequences provide details in the manuscript.

Following the emergence of the Asian lineage of highly pathogenic H5N1 influenza, the WHO along with OIE and FAO created a working group to implement a unified nomenclature system (Anonymous, 2012) . The system set out the minimum requirements to define a new clade including minimum bootstrap value after 1000 bootstrap replicates (≥60) and the average percentage pairwise nucleotide distances between and within clades (>1.5% and <1.5%, respectively). A similar approach could be undertaken for EAV and thus aid in resolving the differences observed between phylogenies. Such a system would allow a standard set of reference sequences to be available representing each clade, which would in turn allow future users to accurately define their data, harmonising the approach for all users.

From the analysis conducted within, we propose a novel nomenclature system for EAV and identify 3 novel clades labelled A, B and G, which are supported via bootstrap analysis. As all such systems the approach is not static: sequences that currently remaining outside the clades may form new clusters or existing clusters expand to include such current outliers. Existing clusters may divide into two or more new clusters (not least depending on criteria to be defined further) or even that clusters collapse (cluster C does not have strong bootstrap support at present for example). Importantly, we think that previous approaches by Mittelholzer et al. (2006) and Zhang et al. (2007) made reference to the nomenclature for PRRSV, which is confusing. The numerical nomenclature system (e.g. EAV-1) conveys the impression of subtypes, sub-species or alike. For EAV no biologically criteria have been identified to define such and serologically there is only one group. However, the strains identified in South Africa (Fig. 3 , cluster G) are the furthest away from all other viruses; this highlights that a second serogroup or (sub)-species may arise in the future. Similarly reference to and EU or NA origin is not warranted for EAV. More than with PRRSV, EAV strains are spreading globally and the historic origins are not traceable: since clinical outbreaks of EAV appear only periodically, it is assumed that EAV circulated already long before its first isolation.

USDA-APHIS VS CEAH, National Health Monitoring System, Fort Collins, CO. Anonymous, 2012. WHO/OIE/FAO H5N1 Evolution Working Group: continued evolution of highly pathogenic avian influenza A (H5N1): updated nomenclature

Equine viral arteritis

Phylogenetic analysis of open reading frame 5 of field isolates of equine arteritis virus and identification of conserved and nonconserved regions in the GL envelope glycoprotein

Neutralisation determinants of laboratory strains and field isolates of equine arteritis virus: identification of four neutralisation sites in the amino-terminal ectodomain of the G(L) envelope glycoprotein

Experiences with infectious cDNA clones of equine arteritis virus: lessons learned and insights gained

EVA outbreak in Argentina

Evaluation of a nested PCR assay for the detection of equine arteritis virus

Nidovirales: a new order comprising Coronaviridae and Arteriviridae

Comparison of M and N gene sequences distinguishes variation amongst equine arteritis virus isolates

The two major envelope proteins of equine arteritis virus associate into disulfide-linked heterodimers

Isolation of a filterable agent causing arteritis of horses and abortions in mares. Its differentiation from equine abortion (influenza) virus

The first isolation of equine arteritis virus in Argentina

Genetic typing of equine arteritis virus isolates from Argentina

Discovery of a small arterivirus gene that overlaps the GP5 coding sequence and is important for virus production

Porcine reproductive and respiratory syndrome virus: genetic diversity of recent British isolates

Comparison of equine arteritis virus isolates using neutralizing monoclonal antibodies and identification of sequence changes in GL associated with neutralization resistance

Lateral transmission of equine arteritis virus among Lipizzaner stallions in South Africa

Genetic variation of ORFs 3 and 4 of equine arteritis virus

A novel subgroup among genotypes of equine arteritis virus: genetic comparison of 40 strains

Comparison of nucleic and amino acid sequences and phylogenetic analysis of the Gs protein of various equine arteritis virus isolates

Structure and genetic analysis of the arterivirus non-structural protein 7alpha

Evidence for absence of equine arteritis virus in the horse population of New Zealand

Extended phylogeny of the equine arteritis virus sequences including South American sequences

Discovery of an RNA virus 3 → 5 exoribonuclease that is critically involved in coronavirus RNA synthesis

Extended phylogeny of equine arteritis virus: division into new subgroups

Discovery of the first insect nidovirus, a missing evolutionary link in the emergence of the largest RNA virus genomes

Transmission of the South African asinine strain of equine arteritis virus (EAV) among horses and between donkeys and horses

Site-directed mutagenesis of the Nidovirus replicative endoribonuclease NendoU exerts pleiotropic effects on the arterivirus life cycle

Formation of the arterivirus replication/transcription complex: a key role for non-structural protein 3 in the remodelling of intracellular membranes

EVA outbreak in New Mexico

Description of the first recorded major occurrence of equine viral arteritis in France

Direct recovery of infectious pestivirus from a full-length RT-PCR amplicon

Generation of recombinant pestiviruses using a full-genome amplification strategy

Arteriviruses

The molecular biology or arteriviruses

Proteolytic processing of the replicase ORF1a protein of equine arteritis virus

Identification of a novel structural protein of arteriviruses

Genetic diversity of equine arteritis virus

Highly diverse type of equine arteritis virus (EAV) from the semen of a South African donkey: short communication

Increased phylogenetic diversity of bovine viral diarrhoea virus type 1 isolates in England and Wales since

MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods

Nuclear localization of non-structural protein 1 and nucleocapsid protein of equine arteritis virus

Arterivirus subgenomic mRNA synthesis and virion biogenesis depend on the multifunctional nsp1 autoprotease

Equine viral arteritis

An infectious arterivirus cDNA clone: identification of a replicase point mutation that abolishes discontinuous mRNA transcription

Deubiquitinase function of arterivirus papain-like protease 2 suppresses the innate immune response in infected host cells

Characterisation of two new structural glycoproteins GP(3) and GP(4), of equine arteritis virus

First recorded outbreak of equine viral arteritis in the United Kingdom

Genetic variation and phylogenetic analysis of 22 French isolates of equine arteritis virus

The work was funded by the Department of Environment Food and Rural Affairs (Defra) (SV3300).