- Research article
- Open access
- Published:
Genomic analysis and replication kinetics of the closely related EHV-1 neuropathogenic 21P40 and abortigenic 97P70 strains
Veterinary Research volume 56, Article number: 12 (2025)
Abstract
Varicellovirus equidalpha 1, formerly known as Equid alphaherpesvirus 1 (EHV-1), is highly prevalent and can lead to various problems, such as respiratory problems, abortion, neonatal foal death, and neurological disorders. The latter is known as equine herpes myeloencephalopathy (EHM). Cases of EHM have significantly increased since the beginning of the twenty-first century. The genomic sequences of five isolates associated with the fatal neurological outbreak in Valencia, Spain, in 2021 were analyzed and documented. The genome and replication kinetics of the Belgian EHM isolate 21P40, associated with the Valencia outbreak, and the well-characterized abortigenic strain 97P70 were compared. Both strains exhibited a nucleotide identity of 99.96%, with only seven genetic mutations in ORFs 13, 24, 30, 32, 40, 65, and 71. Isoleucine and asparagine at loci 291 and 207 of ORF30 (DNA polymerase) and ORF65 (ICP22), respectively, were unique to isolates from the Valencia outbreak. The replication kinetics of these two genetically closely related strains were determined in rabbit kidney (RK-13), equine respiratory, and vaginal mucosal explant cells, as well as equine blood monocytes (CD172a+). Both strains replicated equally well in RK-13 cells. The neuropathogenic isolate 21P40 exhibited a more extensive infection in respiratory explants and blood monocytes, as demonstrated by more plaques and single infected leukocytes, and a higher percentage of infected monocytes. In contrast, vaginal explants infected with the abortigenic strain 97P70 demonstrated more plaques and single infected leukocytes. In conclusion, 21P40 replicated significantly different compared to 97P70 but shared similarities with the Belgian well-studied neuropathogenic EHV-1 strain 03P37.
Introduction
Varicellovirus equidalpha 1, formerly known as Equid alphaherpesvirus 1 (EHV-1), belongs to the subfamily Alphaherpesvirinae causing primarily infections of the upper respiratory tract of horses. The genome of EHV-1 consists of a single linear double-stranded DNA of approximately 150 kbp and comprises unique long (UL) and short (US) regions flanked by identical internal (IR) and terminal (TR) repeats. It carries at least 76 open reading frames (ORFs), of which ORF30 encodes for the DNA polymerase, ORF64 encodes for the immediate early regulatory protein ICP4, ORFs 5, 63, 65, and 77 encode for the early regulatory proteins ICP27, ICP0, ICP22, and IR2, respectively, and ORF12 encodes for the late regulatory protein VP16 [1,2,3,4,5,6,7,8].
The most important impact of EHV-1 infection is the manifestation of respiratory problems in young horses. However, it can also result in abortions, premature neonatal mortality in foals, and a neurological condition referred to as equine herpesvirus myeloencephalopathy (EHM) [9, 10]. The emergence of EHM outbreaks can result in severe economic losses and have negative consequences for the equine industry, as evidenced by recent worldwide outbreaks [11,12,13,14].
During and after the International Valencia Spring Tour held in Spain in February 2021, an EHM outbreak of significant severity was recorded, making it one of the most severe outbreaks in Europe so far. The event in question showcased a total of 752 horses. A total of 18 dead horses were reported and confirmed to be associated with cases in ten countries, including Belgium, Denmark, France, Germany, Italy, Qatar, Spain, Slovakia, Sweden, and Switzerland, in which quarantine regulations were immediately implemented to prevent further dissemination [15].
The occurrence of EHM has been suggested to be associated with a single nucleotide polymorphism (SNP) located at position 2254 of ORF30, encoding the DNA polymerase. This SNP results in the substitution of the original asparagine (N) with aspartic acid (D) at position 752 within the DNA polymerase. Still, previous studies, conducted in several countries, showed a big variability (52 to 97%), linking the N752 (A2254) marker to abortion cases [13, 16,17,18,19,20,21].
Five EHV-1 isolates from horses in Belgium and France, linked to the neurological outbreak during the Valencia tour in 2021, were sequenced and reported in May 2021. All strains exhibited the H250/N752/Y753/K990 genotype of the DNA polymerase and were closely related to each other and phylogenetically grouped into clade 10. In addition to those neuropathogenic isolates, clade 10 contains four sequenced abortigenic isolates, including the well-characterized Belgian abortigenic 97P70 strain [13].
This study aimed to conduct a thorough comparative analysis of the available EHV-1 genomes and replication kinetics of the Belgian neuropathogenic 21P40 strain from the Valencia outbreak and the abortigenic 97P70 strain.
Materials and methods
Viruses
Two Belgian EHV-1 strains were grown in RK-13 cells. The neuropathogenic strain 21P40 was originally isolated in 2021 from a nasal swab of a Belgian horse that had attended the Valencia Tour in February 2021 [13]. The abortigenic isolate 97P70 was first isolated from the lungs of an aborted fetus in 1997 [22].
Genomic analysis
Both strains were sequenced at PathoSense BV using Oxford Nanopore Technology (ONT) on a MinION R.9.4 flow cell [13]. The downstream analysis involved pairwise identity determination and multiple-sequence alignment using the Basic Local Alignment Search Tool (BLAST v2.12.0 +) and Multiple Alignment using Fast Fourier Transform (MAFFT v7.471). The tools were executed with their default configurations. The comparative analysis of genomes encompassed all strains of EHV-1 (n = 127) that were accessible in June 2021 in the GenBank repository and previous studies.
Tissue samples and cells
Respiratory and vaginal explants
Samples of tissues were collected from three healthy horses post-slaughter, specifically from the deep intranasal septum, the proximal segment of the trachea, and the middle portion of the vagina. The tissues were submerged in transport medium consisting of 500 mL of phosphate-buffered saline (PBS) containing calcium and magnesium, supplemented with 100 U/mL penicillin (ThermoFisher Scientific, Paisley, UK), 0.1 mg/mL streptomycin (ThermoFisher Scientific), 0.1 mg/mL gentamicin (ThermoFisher Scientific), 0.1 mg/mL kanamycin (Merck, Darmstadt, Germany), and 0.25 μg/mL amphotericin B (ThermoFisher Scientific) and transported to the laboratory on ice.
Upon arrival, the mucosae were separated from the underlying tissues and divided into uniform explants measuring 25 mm2 using sterile tweezers and surgical blades No. 24 (Swann-Morton, Sheffield, England). To evaluate the viability of the mucosal explants, a single explant per tissue type was embedded in methylcellulose 2% (Methocel®; Sigma-Aldrich, St. Louis, USA) and snap-frozen using dry ice and ethanol 100%. The explants were subjected to an air–liquid interface cultivation method for 18 h at 37 °C with 5% CO2, as previously described by Vandekerckhove et al. in 2009 for the respiratory tissues [23] and Negussie et al. in 2016 for the vagina [24]. Briefly, explants were positioned onto sterilized fine-meshed gauzes with the epithelium facing upwards. The gauzes were then placed in 6-well plates containing serum-free medium, composed of a 50% mixture of Dulbecco’s modified eagle medium including Glutamax (DMEM; ThermoFisher Scientific) and Roswell Park Memorial Institute medium including Glutamax and HEPES (RPMI; ThermoFisher Scientific), supplemented with 100 U/mL penicillin, 0.1 mg/mL streptomycin, 0.1 mg/mL gentamicin and 0.25 μg/mL amphotericin B. To mimic the air–liquid interface, the medium covered the sides of explants, leaving the epithelium exposed to the air.
The ciliary beating of the epithelial cells of the nasal and tracheal explants was assessed by an Olympus IX50 light microscope, along with an in-situ cell death detection kit (Roche Diagnostics Corporation, Basel, Switzerland). This kit, that relies on the TUNEL technology, allowed to determine the number of apoptotic cells in the epithelium and lamina propria of all explants, which is related to the tissue viability.
Blood monocytes (CD172a+)
Twenty milliliters of blood were obtained from three horses through jugular venipuncture and collected into heparin 15 U/mL (Leo, Zaventem, Belgium). The samples were diluted with an equivalent volume of Dulbecco’s phosphate-buffered saline (DPBS) without calcium and magnesium (Gibco, Invitrogen, Paisley, UK). The isolation of peripheral blood mononuclear cells (PBMCs) was carried out through density centrifugation on Ficoll-Paque (density = 1.077 g mL−1) (GE Healthcare, Life Sciences) at 800 × g for 30 min at 18 °C. The PBMCs were obtained from the interphase and subjected to three rounds of washing with DPBS. The cells were suspended in leukocyte medium (LM), composed of RPMI, supplemented with 5% fetal calf serum (FCS, Gibco), 0.1 mM non-essential amino acids (NEAA, Gibco), 1 mM sodium pyruvate (Gibco), 100 U/mL penicillin, 0.1 mg/mL streptomycin, 0.1 mg/mL gentamicin, and 0.25 μg/mL amphotericin B. Subsequently, the cells were cultured on inserts pre-coated with 500 µL of 0.47 mM Poly-L-Lysine (Sigma Aldrich) and placed in 24-well plates (Nunc A/S, Roskilde, Denmark) at a density of one million cells per mL and maintained at 37 °C with 5% CO2. After 12 h, adherent monocytic cells were washed three times with RPMI to eliminate nonadherent lymphocytes.
Rabbit kidney (RK-13) cells
Rabbit kidney (RK-13) cells on inserts (0.47 mM Poly-L-Lysine) in 24-well plates were used as a control culture in this study. They were maintained in RK-13 medium, composed of Modified Eagle's medium (MEM, ThermoFisher), supplemented with 5% FCS, 100 U/mL penicillin, 0.1 mg/mL streptomycin, 0.1 mg/mL gentamicin, and 0.25 μg/mL amphotericin B.
Virus inoculation
To study the replication kinetics, three biological replicates per tissue and cell type were used. The virus stocks used for virus inoculation were at the third and sixth passage for EHV-1 21P40 and EHV-1 97P70, respectively.
Respiratory and vaginal explants
Following overnight cultivation, the explants were inoculated in a 24-well plate by submerging them in 1 mL of inoculum containing 106.5 TCID50 of EHV-1 21P40 or 97P70 and maintained at 37 °C and 5% CO2 for 1 h. Following two rounds of washing with a warm serum-free medium, the explants were repositioned onto their respective gauzes. At 1, 24, and 48 h post-inoculation (hpi), the medium of the explants was collected for extracellular virus titration. Subsequently, the explants were embedded in Methocel®, snap-frozen, and stored at -70 °C. Mock-inoculated explants with only serum-free medium were conducted concurrently.
Blood monocytes (CD172a+) and RK-13 cells
Adherent blood monocytes and RK-13 cells were exposed to 200 µL of EHV-1 21P40 or EHV-1 97P70 in a 24-well plate; two wells per time point, at a multiplicity of infection (MOI) of 5 for 1 h at 37 °C and 5% CO2. The blood monocytes and RK-13 cells were washed twice with a warm respective medium to remove the viral inoculum. At 1, 12, 24, and 48 hpi, cells from the first well were fixed in 100% methanol at -20 °C for 20 min to assess the quantity and percentage of EHV-1 infected cells through immunofluorescent staining. Moreover, the second well’s supernatant and cells were collected for extracellular and intracellular virus titration, respectively. Mock-infected cells with only the respective medium were processed in parallel.
Replication kinetics
Respiratory and vaginal explants
Plaque analysis and penetration through the basement membrane
At corresponding time points, 50 cryosections, each measuring 16 µm, were cut from explants and subsequently fixed in 100% methanol at -20 °C for 20 min. Firstly, endogenous biotin was blocked using the Avidin/Biotin Blocking Kit (Invitrogen). The basement membrane was subsequently stained using monoclonal antibodies (mAb) against mouse anti-collagen VII 1/300 (Sigma-Aldrich), followed by incubation in secondary Texas Red®-labelled goat anti-mouse IgG antibodies 1/100 (Molecular Probes, Eugene). Subsequently, the viral proteins were stained with biotinylated equine polyclonal anti-EHV-1 IgG 0.15 mg/mL, previously prepared in our laboratory [25] (1 h at 37 °C) and streptavidin FITC® 1/200 (Molecular Probes; 1 h at 37 °C). Between the two staining steps and at the end, the cryosections were washed three times with PBS. Finally, the cryosections were mounted with glycerin DABCO (Janssen Pharmaceutica, Beerse, Belgium). The stained cryosections were analyzed using a Leica TCS SP2 laser scanning spectral confocal system (Leica Microsystems, GmbH, Wetzlar, Germany), using an argon 488-nm laser line and a Gre/Ne 543-nm laser line to excite FITC and Texas red, respectively.
Quantification and identification of single EHV-1-infected cells
The number and identity of single EHV-1-infected epithelial cells and leukocytes in the epithelium and the lamina propria were assessed by a double immunofluorescence staining. At corresponding time points, 40 cryosections of 16 µm were prepared and fixed in 100% methanol for 20 min at -20 °C. Twenty cryosections were stained for each cell-surface marker separately. In the first step, the mAb DH59B IgG1 1/100 (VMRD Inc., Pullman, USA) and UC F6G-3 IgG1 1/50 (University of California, Davis, USA) were used as markers for CD172a+ cells from the monocytic lineage or CD3+ cells (pan T-lymphocytes), respectively. Subsequently, cryosections were incubated with secondary Texas Red®-labelled goat anti-mouse antibodies 1/100. In the second step, all cryosections were stained for viral proteins by incubation of biotinylated equine polyclonal anti-EHV-1 IgG 0.15 mg/mL, followed by streptavidin FITC® 1/200. Isotype-corresponding control (IgG1) mouse monoclonal anti-PRV gD antibody 13D12 was included [26]. In general, antibodies were incubated for 1 h at 37 °C. Finally, cryosections were washed three times in PBS and mounted with glycerin-DABCO. The number of individual infected cells and their identity were determined by confocal laser-scanning microscopy. EHV-1 infected cells were counted in five randomly chosen fields of 100 cells.
Virus titration
At each time point, serum-free medium from each well was collected and centrifuged at 400 × g for 10 min at 4 °C and the supernatant was stored at -70 °C. Titration was performed on RK-13 cells and incubated at 37 °C and 5% CO2 and monitored for 7 days. The final virus titer was calculated as 50% tissue culture infective dose (TCID50) according to the Reed and Muench formula (1938) [27].
Blood monocytes (CD172a+) and RK-13 cells
Percentage of EHV-1-infected cells
Inserts with adherent blood monocytes and RK-13 cells were fixed in 100% methanol at -20 °C for 20 min. In the first step, mAb DH59B 1/100 was used as a marker for CD172a+ cells. Subsequently, inserts were incubated with secondary Texas Red®-labelled goat anti-mouse IgG antibodies 1/100. In a second step, inserts were stained for viral proteins by incubation with biotinylated equine polyclonal anti-EHV-1 IgG 0.15 mg/mL, followed by streptavidin-FITC® 1/200. In all steps, antibodies were incubated for 1 h at 37 °C. Finally, inserts were washed three times in PBS and mounted with glycerine-DABCO. The EHV-1-infected blood monocytes and RK-13 were detected, and their percentage was assessed by confocal laser scanning microscopy.
Extracellular and intracellular virus titration
At corresponding sample collection times (1, 12, 24, 48 hpi), both extracellular and intracellular virus titers were determined as described above. The supernatant, containing extracellular virus, was collected, and centrifuged at 400 × g for 10 min at 4 °C prior to storage at -70 °C. Cells containing the intracellular virus were harvested by flushing and scraping the cells in fresh 1 mL of respective cell medium. Afterwards, they were centrifuged at 400 × g for 10 min at 4 °C. The cells were lysed through three cycles of freeze-thawing and then stored at -70 °C. Viral titrations were conducted on RK-13 cells as described previously.
Statistical analysis
For statistical significance, the data were evaluated by ordinary two-way analysis of variance (ANOVA) with a post hoc Tukey’s multiple comparisons test. When p-values were ≤ 0.05, the differences in the results were considered significant. The statistical analysis of the data was done using Prism 9 for macOS, version 9.3.0.
Results
Genomic analysis
Pairwise identity determination revealed a high nucleotide identity between the abortigenic EHV-1 97P70 and neurovirulent EHV-1 21P40 strain (99.96%) with 98% query coverage. The sequence alignment showed differences in seven ORFs, including five SNPs in ORFs 13, 30, 32, 40, 65, one deletion in ORF24 and one insertion in ORF71 (Table 1).
Upon the inclusion of other EHV-1 isolates available in the GenBank repository (n = 127 of which 28 are complete genomes) in the analysis, an amino acid substitution of the isoleucine at position 291 of ORF30 (DNA polymerase) and the asparagine at position 207 of ORF65 (regulator protein ICP22) were unique to the isolates from the Valencia neurological outbreak in 2021, including the 21P40 isolate. On the contrary, the isoleucine at position 180 of ORF32 (DNA packaging terminase subunit 2 ICP18.5) and tyrosine at position 258 of ORF40 (tegument protein) were unique for the abortigenic EHV-1 97P70 strain. The nucleotide substitution at position 927 of ORF13 (Tegument protein VP13/14) represents a silent mutation at the amino acid level at position 309. For ORF24 and 71, also structural variations (insertions/deletions) within their coding sequence were identified among isolates. In Additional file 1, a detailed analysis of the four single-point amino acid substitutions reported in EHV-1 97P70 and 21P40, in addition to locus 752 of the DNA polymerase, compared to other EHV-1 strains in the GenBank repository can be found.
Replication kinetics in respiratory and vaginal explants
The replication of both strains progressed in a plaque-wise manner, spreading laterally. Plaques were detectable from 24 hpi. The plaques induced by both strains increased in number and size but did not traverse the basement membrane at all observed time points (Figure 1A).
Replication kinetics in the respiratory and vaginal mucosa. Representative confocal microscopy images of nasal septum, trachea and vagina infected with EHV-1 abortigenic 97P70 and neuropathogenic 21P40 at 1, 24 and 48 hpi (viral proteins in green and basement membrane in red) A Plaque number and average plaque diameter. B Number of single infected cells C Virus titer. Data are presented as means of three experiments. Asterisk represents statistically significant difference. Scale bar 50 µm. Dashed line represents the limit of detection (LOD).
In the nasal septum mucosa, plaques induced by both strains appeared from 24 hpi and increased at 48 hpi (Figure 1A left). The neuropathogenic EHV-1 21P40 exhibited more plaques at both time points with a significant difference at 48 hpi compared to the abortigenic EHV-1 97P70 (6 ± 3 and 15 ± 5 plaques for EHV-1 21P40 and 2 ± 1 and 9 ± 2 plaques for EHV-1 97P70 at 24 and 48 hpi, respectively). The plaque diameter increased over time for both strains. The abortigenic EHV-1 97P70 strain showed slightly but not significantly bigger (p = 0.77) plaques at 24 hpi (41.00 ± 8.89 µm and 32.73 ± 7.92 µm for EHV-1 97P70 and 21P40, respectively). However, at 48 hpi, the neuropathogenic EHV-1 21P40-induced plaques were significantly bigger (112.10 ± 23.29 µm and 65.20 ± 9.18 µm for EHV-1 21P40 and 97P70, respectively). Single infected cells were visible in the epithelium and lamina propria for both strains from 24 hpi and increased over time (Figure 1B left and Table 2). Both strains infected more CD172a+ monocytic cells in the epithelium than CD3+ cells. In the lamina propria, the neurovirulent EHV-1 21P40 strain infected more CD172a+ monocytic cells than CD3+ cells, while the abortigenic EHV-1 97P70 infected more CD3+ cells than CD172a+ monocytic cells. The titer of both strains showed a steady increase over time with no statistically significant differences between the two strains (Figure 1C left).
In the tracheal mucosa, plaques were observed at 24 hpi for both strains and increased over time (Figure 1A middle). Interestingly, at 24 hpi, the number of EHV-1 97P70 plaques was over two-fold the number of plaques induced by EHV-1 21P40 (29 ± 6 and 13 ± 5 plaques for EHV-1 97P70 and 21P40, respectively). However, at 48 hpi, the EHV-1 21P40 plaques number increased twice to nearly the same number of EHV-1 97P70 plaques, which remained constant (29 ± 7 and 27 ± 4 plaques for EHV-1 97P70 and 21P40). The average plaque diameter increased over time for both strains. EHV-1 97P70 plaques were significantly bigger at 48 hpi (76.20 ± 12.03 µm and 159.17 ± 16.00 µm for EHV-1 97P70 and 71.93 ± 9.25 µm and 127.40 ± 16.75 µm for EHV-1 21P40 at 24 and 48 hpi, respectively). Single infected cells were present in the epithelium and lamina propria for both strains from 24 hpi and increased at 48 hpi (Figure 1B middle and Table 2). The behavior of both strains in the epithelium and the lamina propria was similar to their behavior in the nasal septum. The titer of both strains increased over time, but no statistical significance was observed between the two strains (Figure 1C middle).
In the vaginal mucosa (Figure 1A right), both strains induced plaques at 24 hpi which increased over time. The abortigenic EHV-1 97P70 showed a better replication with over two-fold difference in the number and size of induced plaques (6 ± 3 and 12 ± 6 plaques of 46.90 ± 12.64 µm and 94.77 ± 25.13 µm at 24 and 48 hpi, respectively), with a significant difference at 48 hpi compared to EHV-1 21P40 (2 ± 1 and 5 ± 2 plaques of 21.80 ± 4.56 µm and 40.36 ± 19.96 µm at 24 and 48 hpi, respectively). Both strains infected single cells in the epithelium and lamina propria at 24 hpi and increased at 48 hpi (Figure 1B right and Table 2). More infected CD172a+ monocytic cells than CD3+ cells in the epithelium and lamina propria were observed in both strains. Like the respiratory mucosa, the titer of both strains increased over time with no significant difference (Figure 1C right).
Replication kinetics in blood monocytes (CD172a+) and RK-13 cells
The replication of both strains in blood monocytes was restricted compared to RK-13 cells (Figure 2). The proportion of infected monocytes (Figure 2A) exhibited a temporal pattern whereby it rose steadily, reaching its maximum at 24 hpi, and subsequently declined at 48 hpi. At 12 hpi, both strains exhibited a similar level of infection (2 ± 1% for EHV-1 21P40 and 2 ± 1% for EHV-1 97P70). However, at 24 and 48 hpi, the neuropathogenic 21P40 isolate replicated better in comparison to the abortigenic 97P70 strain (6 ± 2% and 4 ± 1% for EHV-1 21P40 and 3 ± 1% and 2 ± 1% for EHV-1 97P70 at 24 and 48 hpi, respectively) with a statistically significant difference at 24 hpi. On the other hand, the percentage of EHV-1 97P70-infected RK-13 cells showed a steady increase till it reached 98 ± 2% at 24 hpi, while the percentage of EHV-1 21P40-infected RK-13 cells surged from 23 ± 6 to 97 ± 2% at 12 and 24 hpi, respectively. The significant difference between both strains was observed at 12 hpi (51 ± 11% for EHV-1 97P70 and 23 ± 6% for EHV-1 21P40). Throughout the experiment, the restricted replication behavior of both isolates in blood monocytes was consistent (Figure 2B), as manifested by the nearly constant intracellular and extracellular virus titer at different time points. In contrast, in RK-13 cells, the intracellular and extracellular virus titer steadily increased through the experiment time points.
Replication kinetics in blood monocytes and RK-13 cells. Representative confocal microscopy images of blood monocytes and RK-13 cells infected with EHV-1 abortigenic 97P70 and neuropathogenic 21P40 at 1, 12, 24 and 48 hpi (viral proteins in green and CD172a marker in red). A Percentage of infected cells. B Extracellular and intracellular virus titers. Data are presented as means of three experiments. Values with an asterisk are statistically different. At 48 hpi, RK-13 cells were detached. Scale bar 25 μm. Dashed line represents the limit of detection (LOD).
Discussion
Equine herpesvirus-1 (EHV-1) has drawn more attention in the past decade due to recurrent outbreaks of abortions and EHM in China, Europe, New Zealand, and the USA that have caused significant financial losses in the global equine industry. Although experimental and field studies have been conducted to understand the pathogenesis of these abortions and EHM, the mechanism by which the viral factors impact the clinical picture of EHV-1 infection is still poorly understood [11, 13, 28,29,30,31,32,33,34,35,36]. This work sheds new light on the viral factors that are presumably linked to the higher susceptibility to developing EHM.
Previously, a link between a unique SNP in the catalytic subunit of the viral DNA polymerase (ORF30) and EHV-1 outcome had been reported. It was suggested that strains encoding D752 were associated with neurological outbreaks, while strains carrying N752 were more likely linked with non-neurological cases [16, 37, 38]. Surprisingly, this single amino acid substitution (N752D), located in the palm region of the DNA polymerase (ORF30), was not reported in the five strains isolated from the Valencia outbreak and documented in the study of Vereecke et al. [13]. Here, our Valencia strains (including the Belgian 21P40 isolate from the Valencia outbreak) showed a genotype with an isoleucine at position 291 and not the serine within the N-terminal domain of the DNA polymerase (ORF30), which guides the single-stranded template to the polymerase active site [39]. Similar findings were documented in the examination of ORF30 among isolates associated with neurological outbreaks in the United Kingdom between 2012 and 2013, where solely the H250R mutation in the N-terminal domain was identified [38]. Furthermore, it is noteworthy that the EHV-1 isolates identified in an outbreak characterized by a combination of EHM and abortion cases in Germany in 2012 exhibited the G2254/D752 DNA polymerase genotype [33]. A novel EHV-1 DNA polymerase genotype C2254/H752 was identified by Sutton et al. during an outbreak that was documented in France in 2018 and Pusterla et al. during an investigation of EHV-1 isolates from 2019 and 2022 in the USA [39, 40]. These studies confirm that the N752D polymorphism in ORF30 does not exclusively serve as a neuropathogenic marker. Other mutations in the DNA polymerase, other ORFs, or intergenic regions may also contribute to the pathogenesis and outcome of EHV-1. Hence, it is imperative to conduct additional research on ORF30 mutations apart from N752D.
Another unique substitution in Valencia strains was found in regulatory early infected cell protein 22 (ICP22) (ORF65). This mutation involves a change from aspartic acid (D) in the 97P70 strain to asparagine (N) in 21P40 at position 207. This position is located in the acidic-rich area of region 3, which mediates the interaction with immediate early protein 4 (ICP4) and its co-transactivation [41]. Unfortunately, the sequence of ICP22 is only available in a limited number of EHV-1 strains on GenBank, which limits our interpretation of the significance of this mutation in EHV-1 interaction with the host cell and its clinical picture. This is an additional confirmation that the lack of complete EHV-1 genomes minimizes our understanding of the impact of viral factors on the mechanism of EHV-1 infections and the development of EHM [13]. Moreover, to understand the possible impact of mutations in ORF30 and ORF65, as observed in our study, we plan to induce SNP at those two positions in EHV-1 97P70 and determine whether its phenotype changes.
The respiratory mucosal explant model, developed by Vandekerckhove et al. [23], provides a valuable tool as an in vitro model that mimics the in vivo horse situation to have better insights into the pathogenesis of horse respiratory pathogens. This model has been used in previous studies to assess differences in replication between EHV-1 strains, especially between neuropathogenic and abortigenic isolates. In addition, the literature has shown that monocytic cells are the primary carrier cells of EHV-1, which facilitates its dissemination to target organs and has been correlated with the virus's virulence and neuropathogenicity [23, 42,43,44,45,46,47,48,49,50,51]. Therefore, in this study, we used the more complex respiratory and vaginal mucosal explant model in parallel to blood monocytes and an RK-13 control cell line to determine the replication behavior of the neuropathogenic 21P40 and abortigenic 97P70 strains.
Upon conducting a comprehensive assessment of the replication and invasion mechanisms of EHV-1 isolates, differences were noted based on the phenotype of the isolate and the organs utilized. The EHV-1 strains, namely neuropathogenic 21P40 and abortigenic 97P70, exhibited varying levels of infection efficiency and kinetics in different mucosal regions and/or cell types. Specifically, 21P40 demonstrated greater efficacy in the upper respiratory mucosa, whereas 97P70 was more efficient at replicating in the vaginal mucosa. This is in agreement with previous studies to compare the neuropathogenic 03P37 and the abortigenic 97P70 strains in nasal [42, 43, 46] and vaginal mucosa [24].
The CD172a+ monocytic cells were the predominant infected cells in the epithelium, which confirms the role of monocytic cells as the main target cell and carrier for EHV-1. In our study, some EHV-1 infected cells were neither identified as CD172a+ monocytic cells nor CD3+ T-lymphocytes. In the epithelium, the explanation could be that most infected cells are epithelial cells. In the lamina propria, unidentified EHV-1 infected cells may be natural killer (NK) or endothelial cells. It is also possible that CD172a and CD3 might be present but at a low level that is not detectable by immunofluorescence.
Cell-associated viremia is a prerequisite for EHV-1 to disseminate to specific organs, including the central nervous system (CNS) and the pregnant uterus. The manifestation of neurologic disease in horses infected with EHV-1 is associated with a long duration of cell-associated viremia at a sustained and elevated level [17, 18, 37]. Past research by our research group indicated that the neuropathogenic 03P37 variant replicated better in blood monocytes than in the abortigenic 97P70 isolate [51, 52]. Our study revealed a comparable finding wherein more infected monocytes were observed with the neurovirulent 21P40 isolate in contrast to the abortigenic 97P70 strain.
The clade 10 neuropathogenic EHV-1 21P40 strain exhibits distinct replication and invasion patterns in the respiratory and vaginal epithelium, despite having minimal genetic alterations (99.96% nucleotide identity and seven SNPs) when compared to the clade 10 abortigenic EHV-1 97P70. Furthermore, the ability of the 21P40 isolate to infect blood monocytes suggests a higher likelihood of inducing neurological complications. The presence of the ORF30 N752D polymorphism cannot be used as a sole marker for the neuropathogenicity of an isolate. Therefore, it is imperative to investigate other potential genetic changes, such as mutations S291I and D207N in ORFs 30 and 65, respectively. Also, we encourage the scientific community to focus on delivering complete high-quality EHV-1 genomic data, including the under-sampled Unique short (Us) region, flanked by repeat regions IR and UR, as they may be associated with alterations in tropism and replication characteristics.
Availability of data and materials
The accession numbers of all strains used in the study are available in the Additional file 1.
Change history
07 February 2025
A Correction to this paper has been published: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13567-025-01464-5
References
Telford EAR, Watson MS, McBride K, Davison AJ (1992) The DNA sequence of equine herpesvirus-1. Virology 189:304–316
Oladunni FS, Horohov DW, Chambers TM (2019) EHV-1: a constant threat to the horse industry. Front Microbiol 10:2668
Caughman GB, Staczek J, O’Callaghan DJ (1985) Equine herpesvirus type 1 infected cell polypeptides: evidence for immediate early/early/late regulation of viral gene expression. Virology 145:49–61
Gray WL, Baumann RP, Robertson AT, Caughman GB, O’Callaghan DJ, Staczek J (1987) Regulation of equine herpesvirus type 1 gene expression: characterization of immediate early, early, and late transcription. Virology 158:79–87
Bowles DE, Holden VR, Zhao Y, O’Callaghan DJ (1997) The ICP0 protein of equine herpesvirus 1 is an early protein that independently transactivates expression of all classes of viral promoters. J Virol 71:4904–4914
Bowles DE, Kim SK, O’Callaghan DJ (2000) Characterization of the trans-activation properties of equine herpesvirus 1 EICP0 protein. J Virol 74:1200–1208
Buczynski KA, Kim SK, O’Callaghan DJ (1999) Characterization of the transactivation domain of the equine herpesvirus type 1 immediate-early protein. Virus Res 65:131–140
Derbigny WA, Kim SK, Jang HK, O’Callaghan DJ (2002) EHV-1 EICP22 protein sequences that mediate its physical interaction with the immediate-early protein are not sufficient to enhance the trans-activation activity of the IE protein. Virus Res 84:1–15
Ostlund EN (1993) The equine herpesviruses. Vet Clin North Am Equine Pract 9:283–294
Pusterla N, Hussey GS, Goehring LS (2022) Equine herpesvirus-1 myeloencephalopathy. Vet Clin N Am Equine Pract 38:339–362
McFadden AMJ, Hanlon D, McKenzie RK, Gibson I, Bueno IM, Pulford DJ, Orr D, Dunowska M, Stanislawek WL, Spence RP, McDonald WL, Munro G, Mayhew IG (2016) The first reported outbreak of equine herpesvirus myeloencephalopathy in New Zealand. N Z Vet J 64:125–134
Pusterla N, Barnum S, Miller J, Varnell S, Dallap-Schaer B, Aceto H, Simeone A (2021) Investigation of an EHV-1 outbreak in the United States caused by a new H752 genotype. Pathogens 10:747
Vereecke N, Carnet F, Pronost S, Vanschandevijl K, Theuns S, Nauwynck H (2021) Genome sequences of equine herpesvirus 1 strains from a European outbreak of neurological disorders linked to a horse gathering in Valencia, Spain, in 2021. Microbiol Resour Announc 10:e00333-e421
Tong P, Duan R, Palidan N, Deng H, Duan L, Ren M, Song X, Jia C, Tian S, Yang E, Kuang L, Xie J (2022) Outbreak of neuropathogenic equid herpesvirus 1 causing abortions in Yili horses of Zhaosu, North Xinjiang, China. BMC Vet Res 18:83
Termine C, Akerstrom G, Paixão G (2021) Management of an EHV-1 outbreak at FEI events and its international impact. Vet Rec 189:e905
Nugent J, Birch-Machin I, Smith KC, Mumford JA, Swann Z, Newton JR, Bowden RJ, Allen GP, Davis-Poynter N (2006) Analysis of equid herpesvirus 1 strain variation reveals a point mutation of the DNA polymerase strongly associated with neuropathogenic versus nonneuropathogenic disease outbreaks. J Virol 80:4047–4060
Goodman LB, Wagner B, Flaminio MJBF, Sussman KH, Metzger SM, Holland R, Osterrieder N (2006) Comparison of the efficacy of inactivated combination and modified-live virus vaccines against challenge infection with neuropathogenic equine herpesvirus type 1 (EHV-1). Vaccine 24:3636–3645
Van de Walle GR, Goupil R, Wishon C, Damiani A, Perkins GA, Osterrieder N (2009) A single-nucleotide polymorphism in a herpesvirus DNA polymerase is sufficient to cause lethal neurological disease. J Infect Dis 200:20–25
Sutton G, Garvey M, Cullinane A, Jourdan M, Fortier C, Moreau P, Foursin M, Gryspeerdt A, Maisonnier V, Marcillaud-Pitel C, Legrand L, Paillot R, Pronost S (2019) Molecular surveillance of EHV-1 strains circulating in France during and after the major 2009 outbreak in Normandy involving respiratory infection, neurological disorder, and abortion. Viruses 11:916
Smith KL, Allen GP, Branscum AJ, Frank Cook R, Vickers ML, Timoney PJ, Balasuriya UBR (2010) The increased prevalence of neuropathogenic strains of EHV-1 in equine abortions. Vet Microbiol 141:5–11
Perkins GA, Goodman LB, Tsujimura K, Van De Walle GR, Kim SG, Dubovi EJ, Osterrieder N (2009) Investigation of the prevalence of neurologic equine herpes virus type 1 (EHV-1) in a 23-year retrospective analysis (1984–2007). Vet Microbiol 139:375–378
van der Meulen KM, Nauwynck HJ, Buddaert W, Pensaert MB (2000) Replication of equine herpesvirus type 1 in freshly isolated equine peripheral blood mononuclear cells and changes in susceptibility following mitogen stimulation. J Gen Virol 81:21–25
Vandekerckhove A, Glorieux S, Broeck WVD, Gryspeerdt A, van der Meulen KM, Nauwynck HJ (2009) In vitro culture of equine respiratory mucosa explants. Vet J 181:280–287
Negussie H, Li Y, Tessema TS, Nauwynck HJ (2016) Replication characteristics of equine herpesvirus 1 and equine herpesvirus 3: comparative analysis using ex vivo tissue cultures. Vet Res 47:19
van der Meulen KM, Nauwynck HJ, Pensaert MB (2003) Absence of viral antigens on the surface of equine herpesvirus-1-infected peripheral blood mononuclear cells: a strategy to avoid complement-mediated lysis. J Gen Virol 84:93–97
Nauwynck HJ, Pensaert MB (1995) Effect of specific antibodies on the cell-associated spread of pseudorabies virus in monolayers of different cell types. Arch Virol 140:1137–1146
Reed LJ, Muench H (1938) A simple method of estimating fifty per cent endpoints. Am J Hyg 27:493–497
Henninger RW, Reed SM, Saville WJ, Allen GP, Hass GF, Kohn CW, Sofaly C (2007) Outbreak of neurologic disease caused by equine herpesvirus-1 at a University Equestrian Center. J Vet Intern Med 21:157–165
Tsujimura K, Oyama T, Katayama Y, Muranaka M, Bannai H, Nemoto M, Yamanaka T, Kondo T, Kato M, Matsumura T (2011) Prevalence of equine herpesvirus type 1 strains of neuropathogenic genotype in a major breeding area of Japan. J Vet Med Sci 73:1663–1667
Burgess BA, Tokateloff N, Manning S, Lohmann K, Lunn DP, Hussey SB, Morley PS (2012) Nasal shedding of equine herpesvirus-1 from horses in an outbreak of equine herpes myeloencephalopathy in Western Canada. J Vet Intern Med 26:384–392
Pronost S, Legrand L, Pitel P-H, Wegge B, Lissens J, Freymuth F, Richard E, Fortier G (2012) Outbreak of equine herpesvirus myeloencephalopathy in France: a clinical and molecular investigation. Transbound Emerg Dis 59:256–263
Traub-Dargatz JL, Pelzel-McCluskey AM, Creekmore LH, Geiser-Novotny S, Kasari TR, Wiedenheft AM, Bush EJ, Bjork KE (2013) Case–control study of a multistate equine herpesvirus myeloencephalopathy outbreak. J Vet Intern Med 27:339–346
Walter J, Seeh C, Fey K, Bleul U, Osterrieder N (2013) Clinical observations and management of a severe equine herpesvirus type 1 outbreak with abortion and encephalomyelitis. Acta Vet Scand 55:19
Estell KE, Dawson DR, Magdesian KG, Swain E, Laing ST, Siso S, Mapes S, Pusterla N (2015) Quantitative molecular viral loads in 7 horses with naturally occurring equine herpesvirus-1 infection. Equine Vet J 47:689–693
van Galen G, Leblond A, Tritz P, Martinelle L, Pronost S, Saegerman C (2015) A retrospective study on equine herpesvirus type-1 associated myeloencephalopathy in France (2008–2011). Vet Microbiol 179:304–309
Negussie H, Gizaw D, Tessema TS, Nauwynck HJ (2017) Equine herpesvirus-1 myeloencephalopathy, an emerging threat of working equids in Ethiopia. Transbound Emerg Dis 64:389–397
Allen GP, Breathnach CC (2006) Quantification by real-time PCR of the magnitude and duration of leucocyte-associated viraemia in horses infected with neuropathogenic vs. non-neuropathogenic strains of EHV- 1. Equine Vet J 38:252–257
Bryant NA, Wilkie GS, Russell CA, Compston L, Grafham D, Clissold L, McLay K, Medcalf L, Newton R, Davison AJ, Elton DM (2018) Genetic diversity of equine herpesvirus 1 isolated from neurological, abortigenic and respiratory disease outbreaks. Transbound Emerg Dis 65:817–832
Sutton G, Thieulent C, Fortier C, Hue ES, Marcillaud-Pitel C, Pléau A, Deslis A, Guitton E, Paillot R, Pronost S (2020) Identification of a new equid herpesvirus 1 DNA polymerase (ORF30) genotype with the isolation of a C2254/H752 strain in French horses showing no major impact on the strain behaviour. Viruses 12:1160
Pusterla N, Barnum S, Lawton K, Wademan C, Corbin R, Hodzic E (2023) Investigation of the EHV-1 genotype (N752, D752, and H752) in swabs collected from equids with respiratory and neurological disease and abortion from the United States (2019–2022). J Equine Vet Sci 123:104244
He Q, Wu Y, Wang M, Chen S, Jia R, Yang Q, Zhu D, Liu M, Zhao X, Zhang S, Huang J, Ou X, Mao S, Gao Q, Sun D, Tian B, Cheng A (2021) ICP22/IE63 mediated transcriptional regulation and immune evasion: two important survival strategies for alphaherpesviruses. Front Immunol 12:743466
Vandekerckhove AP, Glorieux S, Gryspeerdt AC, Steukers L, Duchateau L, Osterrieder N, Van de Walle GR, Nauwynck HJ (2010) Replication kinetics of neurovirulent versus non-neurovirulent equine herpesvirus type 1 strains in equine nasal mucosal explants. J Gen Virol 91:2019–2028
Gryspeerdt AC, Vandekerckhove AP, Garré B, Barbé F, Van de Walle GR, Nauwynck HJ (2010) Differences in replication kinetics and cell tropism between neurovirulent and non-neurovirulent EHV1 strains during the acute phase of infection in horses. Vet Microbiol 142:242–253
Van Cleemput J, Poelaert KCK, Laval K, Maes R, Hussey GS, Van den Broeck W, Nauwynck HJ (2017) Access to a main alphaherpesvirus receptor, located basolaterally in the respiratory epithelium, is masked by intercellular junctions. Sci Rep 7:16656
Laval K (2016) Equine CD172a+ monocytic cells, the “Trojan horse” for equine herpesvirus type 1 (EHV-1) dissemination in the horse. PhD Thesis, Ghent University, Faculty of Veterinary Medicine, Belgium
Zhao J, Poelaert KCK, Van Cleemput J, Nauwynck HJ (2017) CCL2 and CCL5 driven attraction of CD172a+ monocytic cells during an equine herpesvirus type 1 (EHV-1) infection in equine nasal mucosa and the impact of two migration inhibitors, rosiglitazone (RSG) and quinacrine (QC). Vet Res 48:14
Laval K, Poelaert KCK, Van Cleemput J, Zhao J, Vandekerckhove AP, Gryspeerdt AC, Garré B, van der Muelen K, Baghi HB, Dubale HN, Zarak I, Van Crombrugge E, Nauwynck HJ (2021) The pathogenesis and immune evasive mechanisms of equine herpesvirus type 1. Front Microbiol 12:662686
Laval K, Favoreel HW, Nauwynck HJ (2015) Equine herpesvirus type 1 replication is delayed in CD172a+ monocytic cells and controlled by histone deacetylases. J Gen Virol 96:118–130
Laval K, Favoreel HW, Poelaert KCK, Van Cleemput J, Nauwynck HJ (2015) Equine Herpesvirus type 1 enhances viral replication in CD172a+ monocytic cells upon adhesion to endothelial cells. J Virol 89:10912–10923
Laval K, Favoreel HW, Van Cleemput J, Poelaert KCK, Brown IK, Verhasselt B, Nauwynck HJ (2016) Entry of equid herpesvirus 1 into CD172a+ monocytic cells. J Gen Virol 97:733–746
Laval K, Van Cleemput J, Poelaert KC, Brown IK, Nauwynck HJ (2017) Replication of neurovirulent equine herpesvirus type 1 (EHV-1) in CD172a+ monocytic cells. Comp Immunol Microbiol Infect Dis 50:58–62
Van Cleemput J, Poelaert KCK, Laval K, Nauwynck HJ (2019) Unravelling the first key steps in equine herpesvirus type 5 (EHV5) pathogenesis using ex vivo and in vitro equine models. Vet Res 50:13
Acknowledgements
The authors acknowledge Carine Boone, Charlotte Sonck, and Joridan Rombaut for their excellent technical support and The Egyptian Ministry of Higher Education and Scientific Research for the financial support of Eslam Mohamed.
Funding
This research was funded by the Special Research Fund from Ghent University (BOF, GOA project, grant number GOA013-17), and the Research Foundation—Flanders (FWO, senior project, grant number G035920N). EM was funded by the Egyptian Ministry of Higher Education and Scientific Research. NV was funded by the Flemish Agency for Innovation and Entrepreneurship (Baekeland mandate HBC.2020.2889).
Author information
Authors and Affiliations
Contributions
Conceptualization, EM, KL, and HN; data curation, EM, IZ, and NV; writing—original draft, EM; writing—review and editing, EM, IZ, NV, ST, KL and HN. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Additional information
Handling editor: Marie Galloux
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original version of this article was revised: the authors identified that there is an overlap in Figure 2.
Supplementary Information
Additional file 1. Genomic analysis of EHV-1 isolates available in the GenBank repository
. The four single-point amino acid mutations in ORFs 30, 32, 40 and 65 of EHV-1 97P70 and 21P40, in addition to locus 752 of the DNA polymerase (ORF30), were used in the analysis.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Mohamed, E., Zarak, I., Vereecke, N. et al. Genomic analysis and replication kinetics of the closely related EHV-1 neuropathogenic 21P40 and abortigenic 97P70 strains. Vet Res 56, 12 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13567-024-01434-3
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13567-024-01434-3