Phylogenomic assessment of 23 equid alphaherpesvirus 1 isolates obtained from USA-based equids

Abstract

Abstract Background Equid alphaherpesvirus 1 (EHV-1) is a global viral pathogen of domestic equids which causes reproductive, respiratory and neurological disease. Few isolates acquired from naturally infected USA-based hosts have been fully sequenced and analyzed to date. An ORF 30 (DNA polymerase) variant (A2254G) has previously been associated with neurological disease in host animals. The purpose of this study was to perform phylogenomic analysis of EHV-1 isolates acquired from USA-based hosts and compare these isolates to previously sequenced global isolates. Methods EHV-1 was isolated from 23 naturally infected USA-based equids (6 different states, 15 disease outbreaks) with reproductive (22/23) or neurological disease (1/23). Following virus isolation, EHV-1 DNA was extracted for sequencing using Illumina MiSeq. Following reference-based assembly, whole viral genomes were annotated and assessed. Previously sequenced EHV-1 isolates (n = 114) obtained from global host equids were included in phylogenomic analyses. Results The overall average genomic distance was 0.0828% (SE 0.004%) for the 23 newly sequenced USA isolates and 0.0705% (SE 0.003%) when all 137 isolates were included. Clade structure was predominantly based on geographic origin. Numerous nucleotide substitutions (mean [range], 179 [114–297] synonymous and 81 [38–120] non-synonymous substitutions per isolate) were identified throughout the genome of the newly sequenced USA isolates. The previously described ORF 30 A2254G substitution (associated with neurological disease) was found in only one isolate obtained from a host with non-neurological clinical signs (reproductive disease), six additional, unique, non-synonymous ORF 30 substitutions were detected in 22/23 USA isolates. Evidence of recombination was present in most (22/23) of the newly sequenced USA isolates. Conclusions Overall, the genomes of the 23 newly sequenced EHV-1 isolates obtained from USA-based hosts were broadly similar to global isolates. The previously described ORF 30 A2254G neurological substitution was infrequently detected in the newly sequenced USA isolates, most of which were obtained from host animals with reproductive disease. Recombination was likely to be partially responsible for genomic diversity in the newly sequenced USA isolates.