Comparative Genomic Sequencing and Pathogenic Properties of Equine Herpesvirus 1 KyA and RacL11

Equine herpesvirus 1 (EHV-1) is a major pathogen affecting equines worldwide. The virus causes respiratory disease, abortion, and, in some cases, neurological disease. EHV-1 Kentucky A (KyA) is attenuated in the mouse and equine, whereas wild-type pathogenic strain RacL11 induces severe inflammatory infiltration of the lung, causing infected mice to succumb. The complete DNA sequencing of the KyA genome revealed that genes UL17 (ORF17), US6 (ORF73; gI), US7 (ORF74; gE), and US8 (ORF75; 10 K) are deleted as compared to the RacL11 and Ab4 genomes. In-frame deletions in the US1 (ORF68), US4 (ORF71; gp2), and UL63 (ORF63; EICP0) genes and point mutations in 14 different open reading frames (ORFs) were detected in the KyA genome. Interestingly, UL1 (ORF1) and UL2 (ORF2) were deleted in both KyA and RacL11. Our previous studies showed that EHV-1 glycoproteins gI, gE, and full-length gp2 contribute to the pathogenesis of the RacL11 strain. The confirmation of these gene deletions in KyA suggests their contribution to the attenuation of this virus. The growth kinetics results revealed that KyA replicates to high titers in cell culture as compared to RacL11 and Ab4, indicating that the above genomic deletions and mutations in KyA do not have an inhibitory effect on KyA replication in cells of mouse, rabbit, equine, or human origin. Studies of EHV-1 pathogenesis in CBA mice showed that KyA is attenuated whereas mice infected with RacL11 succumbed by 3–6 days post-infection, which is consistent with our previous results.

Introduction

Equine herpesvirus 1 (EHV-1) is a common pathogen of horses worldwide. EHV-1 is a member of the Varicellovirus genus within the subfamily Alphaherpesvirinae (1, 2) Infection causes upper respiratory tract disease, abortion in pregnant mares, neonatal death, and fatal myeloencephalopathy (3–5). The role of viral proteins ORF30 and glycoprotein D in the pathogenicity of myeloencephalopathy was reported recently as these two viral proteins are responsible for modulating host immune response (6). One of the main causes of abortion and neurological disease in horses is thrombosis in placental and spinal cord vessels. The findings of Stokol et al. (7) suggested that activation of platelets is a feature of EHV-1 pathogenicity and could be a factor leading to thrombosis. Important in EHV-1 pathogenicity is its ability to infect lymphoid tissue of the upper respiratory tract and spread to monocytes that mediate viremia (8, 9). By this means, the virus can reach the vasculature of the uterus and the central nervous system, resulting in severe clinical outcome (9, 10).

Several vaccines are in current use to prevent the most severe threat of EHV-1 infection. However, the currently available vaccines do not confer long-term protection against the severe manifestations of EHV-1 disease, such as fatal myeloencephalopathy. Therefore, there is an urgent need to develop vaccines that provide long-term protection. A better understanding of EHV-1 pathogenesis in relation to virus genomic sequences and specific alterations in the genome will help in the design of effective vaccines to reduce mortality and morbidity.

Equine herpesvirus 1 RacL11, isolated from an aborted foal (11), causes a severe infection of mice, resulting in significant weight loss, lethargy, ruffled fur, inflammation of the lungs, and a high mortality rate (12–14). The Ab4 strain of EHV-1 was isolated from a quadriplegic gelding (15), and it has been reported that deletion of gene US4 (ORF71) in the Ab4 genome resulted in attenuation. The gene 71-deleted virus mutant was evaluated as an experimental vaccine and was found to confer protection against pulmonary disease in mice after challenge with wild-type (wt) EHV-1 (16). The EHV-1 strain Kentucky A (KyA), a candidate vaccine strain, is non-pathogenic for both mice and horses after serial passage in mouse fibroblast L-M cells (17–20). The repeated passage in cells other than those of the natural host has resulted in the accumulation of genomic alterations of the KyA chromosome, including deletion of several genes or portions of open reading frames (ORFs).

The structure of the EHV-1 genome was first elucidated by detailed restriction enzyme mapping (21) and was confirmed by electron microscopic analyses (22). The EHV-1 genome comprises a unique long (U_L) region and a short (S) region that consists of a unique segment (Us) bracketed by identical internal and terminal repeat, designated as IR and TR, respectively. The complete genomic sequence of the EHV-1 Ab4 was reported (23), and found to be 150,223 bp in length and to contain 76 distinct ORFs (23, 24). Six identical ORFs are present in the internal and terminal repeats of the EHV-1 genome as regulatory genes IR2 (ORF77) (25, 26) and IR3 (ORF78) (27–29) were identified in the KyA genome in addition to genes 64, 65, 66, and 67 as reported for Ab4 by Telford et al. (23).

Our previous DNA sequencing analyses of portions of the KyA genome revealed several gene deletions UL1 [ORF1], UL2 [ORF2], UL17 [ORF17], US6 [ORF73; gI], US7 [ORF74; gE], and US8 [ORF75] and in-frame deletions in UL63 (ORF63; EICP0) and US4 (ORF71; gp2) relative to the sequence of the wt Ab4 strain of EHV-1 (30–35). EHV-1 glycoproteins gI, gE, and gp2 contribute to the pathogenesis of the wt RacL11 strain (36–38). The deletion of ORF1/2 was found to impact clinical disease as ponies infected with the UL1/UL2 deleted virus had significantly shorter primary pyrexia, reduced nasal shedding, and a decrease in PBMC IL-8 as compared to ponies infected with wild-type virus (39). A single nucleotide polymorphism, N752 replaced by D752, in the gene encoding the viral DNA polymerase (ORF30) played a significant role in the magnitude and longer duration of viremia and enhancement of neurovirulence of naturally occurring EHV-1 strains (40).

The main objective of the present study was to determine the complete nucleotide sequence of the genomes of the attenuated EHV-1 KyA strain and the wt RacL11 strain and to compare these genomes in order to identify mutations associated with viral growth, cell tropism, and virulence. Comparison of the KyA and RacL11 genomes and additional comparison of these two genomes to that of wt Ab4 will help to identify viral gene(s) that are required for virulence and may trigger host cell gene expression, including innate immune responses. In addition, we determined the growth kinetics of the EHV-1 KyA, RacL11, and Ab4 strains in four different cell types, information that may eventually be helpful in the design of efficacious vaccines.

Materials and Methods

Cells and Viruses

Equine dermis (NBL6), rabbit kidney (RK-13), human epithelial kidney (HEK293), and mouse lung epithelial (MLE12) (ATCC: CRL-2110) cells were used for one-step growth experiments. NBL6, RK-13, and HEK293 cells were maintained in Eagle’s minimal medium supplemented with 100 units of penicillin/ml, 100 µg of streptomycin/ml, 1× non-essential amino acids, and 5% (or 10%) fetal bovine serum (FBS). MLE12 cells were maintained in HITES medium (ATCC) containing 2% FBS. Three EHV-1 strains KyA, RacL11, and Ab4 were used in the study. All three strains were grown in NBL6 cells.

Mice

Three- to four-week-old female CBA mice were obtained from Harlan Sprague-Dawley (Indianapolis, IN, USA) and maintained in filter-topped cages at Animal Resource Facility of the Louisiana State University Health Science Center (LSUHSC), Shreveport. All mice were rested for 5 to 7 days prior to use.

Infection and Assessment of Pathogenicity of EHV-1 Strains

To assess viral virulence, CBA mice (n = 5 per group) were anesthetized by isoflurane (Sigma Chemical Co., St. Louis, MO, USA) inhalation and were inoculated intranasally either with of 1.5 × 10⁶ plaque forming unit (PFU)/mouse of KyA, wild-type (wt) RacL11, or wt Ab4. Mice were monitored daily for observable clinical signs, such as abnormal respiration, crouching, ruffled fur, and huddling. All mice were also monitored for the loss of body weight. Animals exhibiting profound deterioration or a progressive decline in mobility were euthanized to prevent suffering. The significance of the percent survival to judge virulence was determined by GraphPad Prism software (GraphPad Software, Inc., La Jolla, CA, USA). Mice were sacrificed by CO₂ inhalation, the lungs were removed at 2, 3, and 4 days post-infection (dpi), and virus titers were determined by plaque assay on NBL6 cells.

Growth Kinetics in Cell Culture

NBL6, RK-13, HEK293, and MLE12 cells (8 × 10⁵) were seeded in six-well plates and infected at a multiplicity of infection (MOI) of 0.05. After 1 h of attachment at 37°C in 5% CO₂, cells were washed with medium three times and overlaid with Dulbecco’s modification of Eagle’s medium (DMEM; Corning Inc., Manassas, VA, USA) containing 10% FBS. The plates were incubated in 5% CO₂ at 37°C and harvested at 2, 12, 24, 36, and 48 h post-infection. Virus quantification was done by plaque assay on equine NBL6 cells. Briefly, freshly grown confluent NBL6 cells in 24 well plates (Corning, Inc., Manassas, VA, USA) were infected with a 10-fold serial dilution of each virus. After 1 h adsorption at 37°C, virus-infected cells were covered with medium containing 1.25% methylcellulose and incubated in 5% CO₂ at 37°C. After 4 dpi, cells were fixed with 10% formalin and stained with 0.5% crystal violet to visualize the plaques. Growth kinetics results of KyA, RacL11, and Ab4 were statistically evaluated by analysis of variance (ANOVA) using GraphPad PrismTM version 5.0 software (GraphPad Software Inc., San Diego, CA, USA).

Viral DNA Extraction and High-Throughput Sequencing

The genomes of KyA and RacL11 strains of EHV-1 were sequenced. EHV-1 KyA is attenuated in the mouse and equine (19, 20), whereas wild-type pathogenic strain RacL11 induces severe inflammatory infiltration of the lung causing infected mice to succumb (12). Total viral DNA was extracted from the supernatant of NBL6 cells infected with EHV-1 KyA or RacL11 using the QiAampMinElute virus spin kit (Qiagen, Valencia, CA, USA), according to the manufacturer’s instructions. Next-generation DNA sequencing was done on an Ion-torrent platform at PrimBio Research Institute, Philadelphia, PA, USA. The libraries of the EHV-1 strains were prepared using Ion Xpress Plus Fragment Library Kit (Thermo Fisher, South San Francisco, CA, USA). Briefly, 100 ng of virus DNA was enzymatically fragmented. Ion P1 adapter and barcodes were ligated to the fragmented DNA, and size selection was performed to obtain the fragments in the range of 200–350 bp. The libraries were further amplified with the Library Amplification Primer Mix. The final libraries were analyzed using the Agilent High Sensitivity DNA Analysis Kit (Agilent, Santa Clara, CA, USA). The library (26 pM) was templated and enriched using Ion PGM Hi-Q OT2 kit (Thermo Fisher). Sequencing was performed using an Ion Torrent PGM Sequencer and the Ion PGM Hi-Q Sequencing Kit (Thermo Fisher) with 850 sequencing flows. Point mutations were confirmed by Sanger sequencing at the DNA Core Facility of Iowa State University.

Sequencing Data Analysis

DNA sequence assembly was done by AssemblerSPAdes (v5.0.0.0) and Mira Assembly. The complete genome sequences were assembled by reference sequence mapping with Mira Assembly using the genomic sequence of EHV-1 Ab4 strain (GenBank accession number AY665713). Alignment and comparison of DNA sequences were done by ClustlW2 and BioEdit sequence alignment editor. The phylogenetic tree was constructed to compare the genomes of KyA and RacL11 to those of 27 other equine herpesvirus strains of EHV-1, EHV-4, and EHV-9. The nucleotide sequences of the 27 other equine herpesvirus strains were obtained from the GenBank database. The phylogenetic tree was constructed by the neighbor-joining method using the Molecular Evolutionary Genetic Analysis (MEGA7) software (41) based on the nucleotide sequences. Node support was evaluated using 1,000 bootstrap replicates.

Nucleotide Sequence Accession Numbers

The genomic sequences of KyA and RacL11 strains described in this work were deposited in GenBank under accession numbers MF975655 and MF975656, respectively.

Results

To elucidate the genetic relationships between the attenuated EHV-1 KyA and wild-type (wt) RacL11 and identify genomic mutations and alterations possibly associated with attenuation and virulence, whole genome sequencing was done. The EHV-1 Ab4 sequence was used as a reference for mapping and annotation. The sequencing results showed that the KyA genome is 141,350 bp and the RacL11 genome is 147,469 bp as compared to the wt EHV-1 Ab4 genome of 150,223 bp (Figure 1A). The comparison of sequences with Ab4 showed that the KyA and RacL11 genomes are 8,873 bp and 2,754 bp shorter, respectively. The G+C content of the KyA and RacL11 genomes is 57 and 56%, respectively. The G+C content of the Ab4 genome is 56.7% (23). A high G+C content is one of the features of the herpes simples virus (HSV) genomes as HSV-1 and HSV-2 have a G+C content of 68 and 70%, respectively, while the human genome has a G+C content 41% (42). Table 1 presents a comparison of the gene products of EHV-1 Ab4, RacL11, and KyA. The genomic sequencing data of this study and the results of a large number of publications on EHV-1 gene functions indicate that the wild-type EHV-1 chromosome harbors 63 ORFs in the UL region, 6 ORFs in each inverted repeat (IR and TR), and 9 ORFs in the Us segment. Although the function of many of the 78 distinct EHV-1 ORFs remains to be elucidated, studies indicate that most EHV-1 gene products characterized to data have functions similar to those of their homologs in the genome of HSV-1.

FIGURE 1

Figure 1. Schematic illustration of deletions in the genomes of equine herpesvirus 1 (EHV-1) Kentucky A (KyA), RacL11, and Ab4 strains. EHV-1 Ab4 sequences were used as the reference for mapping and annotation. (A) The two unique genomic segments, unique long (U_L) and unique short (U_S), and the internal and terminal repeat sequences (IR and TR, respectively) are shown. The sequencing results revealed that the genomes of KyA and RacL11 are 141,350 bp and 147,469 bp, respectively. nt, nucleotide. (B) Sequence analysis of the U_L region showed that the UL1 (ORF1) and UL2 (ORF2) genes are deleted in both KyA and RacL11. Deletion of UL17 was also observed in KyA. An identical in-frame deletion in UL63 (ORF63; EICP0) of both RacL11 and KyA was detected. The sequence analyses of the U_S segment showed that US6 (ORF73; gI), US7 (ORF74; gE), and US8 (ORF75; 10 K) genes were deleted in KyA. An in-frame deletion in the US1 (ORF68) and US4 (ORF71; gp2) genes was also detected in the KyA genome.

TABLE 1

Table 1. Comparison of gene products of equine herpesvirus 1 (EHV-1) Ab4, RacL11, and Kentucky A (KyA).

To assess the point mutations, the amino acid (aa) sequences of all ORFs were compared between the KyA and RacL11 strains (Table 2). Thirty six point mutations were detected in 14 different ORFs of the KyA genome. Fifty nine ORFs were completely conserved among the sequences of KyA and RacL11 examined. Interestingly, both the IR2 (ORF77) gene (25, 26) and the IR3 (ORF78) ORF (27–29) that were discovered to be regulatory genes in the KyA genome are conserved in both the RacL11 and Ab4 genomes.

TABLE 2

Table 2. Amino acid substitutions in the Kentucky A (KyA) and RacL11 proteins compared to Ab4.

DNA Sequence Analysis of the U_L Region

The EHV-1 genome is composed of a unique long (U_L) region and a unique short segment (U_S) that is flanked by two inverted repeats (IR and TR (21);). The U_L region of EHV-1 Ab4 is 112,935 bp that includes the 32 bp repeats. The sequencing results revealed that the U_L region of KyA totals to 109,697 bp that includes an identical 32 bp repeat at each end. The U_L region of RacL11 totals to 111,148 bp that also includes the identical 32 bp repeat at each terminus. In the KyA chromosome, the UL1 (ORF1), UL2 (ORF2), and UL17 (ORF17) genes are deleted as compared with the Ab4 sequence (Figures 1 and 2). These three major deletions in the U_L region of KyA have been reported (30, 31, 34) and were confirmed in this study. With regard to the RacL11 strain, the sequence of the ORF17 gene was identical to that of the Ab4 virus, and the ORF would encode a protein of 401 residues (Figure 2). Interestingly, the UL1 and UL2 genes were also deleted in the RacL11 strain. The deletion of these two ORFs in the RacL11 genome was surprising as the UL1 and UL2 genes were reported to be important for disease outcome and the modulation of cytokine responses (39).

FIGURE 2

Figure 2. ORF17 and open reading frame (ORF) 68 in the Kentucky A (KyA) genome. (A) Predicted amino acid sequence of the ORF17 protein of RacL11 compared to that of Ab4 reveal perfect identity. Sequence analysis of the KyA genome revealed that ORF17 is deleted in KyA genome. (B) An identical in-frame deletion was detected in the US1 (OFR68) gene of the KyA and RacL11 genomes. Clustl W multiple sequence alignment software (http://www.ebi.ac.uk/Tools/msa/clustalw2/) was used for the alignment.

ORF17 is non-essential for EHV-1 replication in cell culture (34). An identical in-frame deletion in the UL63 (ORF63; EICP0) gene of KyA and RacL11 genomes was observed (Figure 1 and Table 3). The results revealed that sequences encoding 113aa were deleted in the EICP0 gene of both strains. Our previous results showed that the KyA EICP0 protein functions as a potent transcriptional activator of all classes of EHV-1 promoters and, thus, the 113 residues (aa319–432) of the Ab4 EICP0 are not essential for its transactivation activities (35, 36). The sequencing results also showed that the addition of eight nucleotides AGCCGCCC at position 1,847 was detected in ORF14 of both KyA and RacL11 as compared with ORF14 of Ab4. Due to the addition of these nucleotides, the ORF14 of KyA and RacL11 was shifted such that the final length of the ORF14 protein is 626aa in contrast to the Ab4 ORF14 protein of 747aa.

TABLE 3

Table 3. Deletions in the genomes and proteins of EHV-1 Kentucky A (KyA) and RacL11 as compared to equine herpesvirus 1 (EHV-1) Ab4.

DNA Sequence Analysis of the U_S Segment

The U_S segment of EHV-1 Ab4 is 11,861 bp long and is flanked with the IR and TR of 12,714 bp each (23). The DNA sequencing results showed that the KyA U_S segment is 6,759 bp flanked with IR and TR sequences of 12,446 bp each, while the U_S segment of RacL11 is 11,872 bp flanked with IR/TR of 12,223 bp each. EHV-1 KyA has a total of 268 bp deletions in each inverted repeat, located between ORF63 and ORF64. RacL11 has a total of 491 bp deletions in each inverted repeat, located upstream of ORF64 and upstream of ORF65. The genome of EHV-1 Ab4 does not have these deletions. DNA sequence analysis of the U_S segment of KyA revealed that genes US6 (ORF73; gI), US7 (ORF74; gE), and US8 (ORF75; 10 K) were deleted as compared to the sequences of Ab4 and RacL11 (Figure 1 and Table 3). An in-frame deletion in the KyA US4 (gp2) gene was detected. An identical in-frame deletion in the US1 (OFR68) gene of the KyA and RacL11 genomes was also detected, and this deletion would result in the ORF68 protein of KyA and RacL11 being 246aa in contrast to the 418aa ORF68 protein of Ab4 (Figure 2B). Whether the difference in the amino acid sequence of the carboxyl terminus of the KyA and RacL11 ORF68 protein affects the function of this virion protein remains to be determined. EHV-1 glycoproteins gI, gE, and full-length gp2 are involved in cell-to-cell spread and contribute to the pathogenesis of EHV-1 (37, 38, 50). These results suggest that the deletion of gI, gE, and gp2 in KyA has important roles in the attenuation of this EHV-1 strain.

Growth Kinetics of KyA, RacL11, and Ab4 in Cell Culture

One-step growth experiments of the three EHV-1 strains, KyA, RacL11, and Ab4, were conducted to examine whether mutations detected in KyA affect viral growth in mouse MLE12 rabbit RK-13, equine NBL6, and human HEK293 cells. The monolayers were infected with each of the three EHV-1 strains at a MOI of 0.05 and were harvested at 0, 12, 24, 36, and 48 h post-infection (hpi). Culture medium was also collected at the same time points. Extracellular and intracellular virus titers were determined by plaque assay on NBL6 cells. KyA and RacL11 exhibited identical growth kinetics in NBL6 cells (Figures 3A,B). Whereas the titer of Ab4 was 2 logs lower than those of KyA and RacL11 (Figure 3A). In RK-13 cells, the growth patterns of the RacL11 and Ab4 strains were similar, although the intracellular titer of KyA was higher than those of the RacL11 and Ab4 strains (Figures 3C,D). In HEK293 cells, the RacL11 and Ab4 showed similar growth patterns (Figures 3E,F). By contrast, the intracellular titer of KyA was significantly higher at all time points (Figure 3E). In MLE12 cells, the extracellular and intracellular titers of KyA significantly increased rapidly with time while Ab4 was not able to grow in these mouse cells (Figures 3G,H). Overall, these results demonstrated that deletions and point mutations in the KyA genome presented in this study did not impair KyA growth or tropism for any of the four cell types.

FIGURE 3

Figure 3. Growth kinetics of equine herpesvirus 1 (EHV-1) strains in four cell types. Intracellular and extracellular growth kinetics of EHV-1 Kentucky A (KyA), RacL11, and Ab4 in equine NBL6 cells (A,B), rabbit RK-13 cells (C,D), human HEK293 cells (E,F), and mouse MLE12 cells (G,H). Cell monolayers were infected at a multiplicity of infection of 0.05. Culture medium and cells were harvested at 2, 12, 24, 36, and 48 h post-infection. The intracellular virus was released by three freeze–thaw cycles. Virus titer was determined by plaque assay on NBL6 cells. Data are averages and are representative of three independent experiments. Error bars indicate SD. Statistical differences between KyA and RacL11 (or Ab4) viruses were labeled (*P < 0.05, **P < 0.01, and ***P < 0.001).

Pathogenicity of EHV-1 KyA, RacL11, and Ab4 strains in CBA Mice

CBA mice were infected intranasally with KyA, RacL11, or Ab4 (1.5 × 10⁶ PFU/mouse) and were monitored daily for weight gain or loss and morbidity. Mouse lungs were harvested and homogenized following EHV-1 infection, and the amount of infectious virus was determined. RacL11-infected mice lost more than 20% of their pre-infection body weight by day 4 pi and succumbed by 3 to 6 dpi (Figures 4A,B). By contrast, KyA- and Ab4-infected mice lost less than 10% of their pre-infection body weight and rapidly regained their weight. None of the mice succumbed to KyA infection (Figures 4A,B). Virus titers in the lungs of KyA-infected mice were 100-fold lower at 4 dpi than those of the lungs of RacL11-infected mice (Figure 4C). These results demonstrate that EHV-1 KyA and Ab4 are apathogenic in the CBA mouse model and are consistent with previous studies (14, 37, 38, 50, 51).

FIGURE 4

Figure 4. Pathogenicity of equine herpesvirus 1 (EHV-1) Kentucky A (KyA), RacL11, and Ab4 in CBA mice. Mice were infected intranasally with 1.5 × 10⁶ plaque forming unit of KyA, RacL11, or Ab4. Cell culture medium was used for mock infection. (A) Mice (n = 5 per group) were weighed daily at days 1–8. The body weights were expressed in average percentage gain or loss relative to the initial body weight. (B) Percentage survival at days post-infection in mock- or virus-infected mice. (C) Virus titers of lungs of EHV-1-infected CBA mice. At days 2, 3, and 4 pi, lungs were removed and homogenized, and the amounts of KyA, RacL11, and Ab4 were determined by standard plaque titration on NBL6 monolayers.

Phylogenetic Analysis

A phylogenetic tree was constructed to compare the genomes of KyA and RacL11 with the genomes of 27 other equine herpesviruses, including EHV-1, EHV-4, and EHV-9 (Figure 5). In the phylogenetic analysis, the 20 EHV-1 strains fell into two genetic groups. Group I consisted of 16 strains, while group II consisted of 4 strains. The KyA and RacL11 strains belonged to group I, and the analyses revealed that strains KyA and RacL11 are closely related to strains Ab4, 00c19, 89c105, 01c1, 90c16, VA02, NY05, NMKT04, and FL06. The eight strains of equine herpesvirus 4 were closely related and separated into a single group, group IV. Lastly, the sole strain of equine herpesvirus 9 fell into a separate group, designated as group III. Figure 5B presents the genetic relationships among EHV-1 strain of group 1 and reveals that KyA has a greater number of genetic changes as compared to RacL11 and Ab4.

FIGURE 5

Figure 5. Phylogenetic analyses of the complete genomes of strains of equine herpesvirus 1 (EHV-1), EHV-4, and EHV-9. (A) The phylogenetic tree indicating evolutionary relationships was generated from the genomic sequence alignments of the Kentucky A (KyA) (GenBank: MF975655), RacL11 (GenBank: MF975656), and 27 other EHV strains using Molecular Evolutionary Genetic Analysis software (41). (B) Subtree of group I reveals EHV-1 KyA and RacL11 genomes are closely related but the KyA genome exhibits a greater evolutionary distance. Bootstrap values out of 1,000 replicates are indicated as a percentage to the left of each branch of the tree. The bar at the bottom of the figure provides a scale that represents the amount of genetic change. The units of branch length are nucleotide substitutions per site divided by the length of the genomic sequence.

Discussion

It is noteworthy that the present study reveals the whole-genome sequence of attenuated EHV-1 KyA and wild-type (wt) RacL11 strains, both of which show several gene alterations as compared to the genome of wt EHV-1 Ab4 (Tables 2 and 3, Figures 1 and 2; KyA and RacL11 genomes GenBank MF975655 and MF975656, respectively). The sequencing results of the U_L region of KyA revealed genomic deletions of UL1 (ORF1), UL2 (ORF2), and UL17 (ORF17). Our previous sequencing results of portions of the KyA genome showed genomic deletions of UL1, UL2, UL17, US6 (ORF73; gI), US7 (ORF74; gE), and US8 (ORF75) and in-frame deletions in genes US1 (ORF68), US4 (ORF71; gp2), and UL63 (ORF63; EICP0) as compared to the sequences of Ab4 (30–35). The roles of the UL1 and UL2 in EHV-1 have been investigated (39), and the data showed that deletion of these two genes from the EHV-1 strain Ab4 attenuated disease despite the fact that both wt and mutant viruses showed similar in vitro growth kinetics. These results suggest that the deletion of the UL1 and UL2 genes is involved in the attenuation of the KyA strain.

We confirmed an identical in-frame deletion in the EICP0 gene in the genomes of KyA and RacL11 such that the EICP0 regulatory protein of these two strains would be 419 amino acids in contrast to the 532-amino acid protein of Ab4. The demonstration that the EICP0 protein of KyA has potent transactivation activity of all classes of EHV-1 promoters (35, 36) indicates that residues 319–432 of the Ab4 EICP0 protein are not essential for its transactivation activity. Our sequence analyses of the U_S segment showed genomic deletions of gI, gE, 10 K, and gp2 in KyA as compared to the Us of RacL11 and Ab4. The role of gI and gE in viral pathogenesis is well established. Previously, it was demonstrated that the restoration of gI and gE to attenuated KyA generated the KgI/gE virus that exhibited a virulence phenotype in the CBA mouse (38). The KgI/gE virus was able to spread to the brains of CBA mice after intranasal infection, resulting in meningoencephalitis (52), which is consistent with the role of gI and gE observed in EHV-1-infected horses exhibiting neurological signs (9). However, the deletion of glycoproteins gI and gE from wt RacL11 did not significantly reduce clinical signs of respiratory disease in mice (38), suggesting that multiple viral genes are involved in the virulence phenotype. The deletion of gene US4 (gp2) in EHV-1 Ab4 resulted in its attenuation, and the generated virus conferred protection against pulmonary disease in mice after challenge with wt virus (16). The insertion and expression of the full-length RacL11 gp2 (791aa) in KyA generated a recombinant virus that elicited fatal immunopathological responses in the lower respiratory tract of CBA mice, indicating a role of US4 in EHV-1 pathogenesis (50). Overall, we conclude that the six genomic deletions UL1, UL2, UL17, gI, gE, and US8 and two in-frame deletions in ECIP0 and gp2 may contribute to the attenuation of EHV-1 KyA.

The data presented in this study revealed that the KyA strain replicated to titers that exceeded those observed for strains RacL11 and Ab4 in all four cell types examined. The growth kinetics results revealed that mutations in KyA have no inhibitory effect on its replication, and the rapid replication of KyA to high titers makes this strain an attractive vaccine candidate. Collaborative studies with Japanese researchers demonstrated that the ΔgEΔgI-lacZ mutant of a pathogenic strain of EHV-1 replicated to high titers in cell culture, indicating that gE and gI are not required for viral growth in vitro (37). In another study, it was shown that deletion of UL1 and UL2 in EHV-1 Ab4 did not alter in vitro growth kinetics (39).

Our findings showed that all mice infected with RacL11 succumbed by 6 dpi while the KyA- or Ab4-infected mice survived. KyA and RacL11 were able to replicate in the lungs of CBA mice while Ab4 did not (Figure 3C). EHV-1 Ab4 is pathogenic in foals (39, 53) and is not pathogenic in mice because this strain is not mouse-adapted and cannot replicate in the lungs of mice. However, if the infectious dose of Ab4 was increased 10-fold to 1.5 × 10⁷ PFU, the CBA mice lost more body weight and 20% of the animals succumbed to infection by 6 dpi (data not shown). When infected with high dose (4 × 10⁷) of KyA, all the mice survived and showed weight loss/gain patterns similar to those of mock-infected mice (19). The observation that KyA and Ab4 are apathogenic in mice was also shown by van Woensel et al. (51). These results are consistent with our previously published findings that KyA is non-pathogenic in CBA mice and horses, whereas the wt RacL11 induces severe inflammatory cell infiltration in the lung, such that infected mice succumb at 4–6 dpi (14, 20, 38, 50).

In the phylogenetic study, KyA, RacL11, and Ab4 strains belong to group I (Figure 5). The group I consists of 16 EHV-1 strains including 11 EHV-1 strains sequenced in Japan and two strains T953 and T953-P15 sequenced in USA. The phylogenetic tree showed that both EHV-1 KyA and RacL11 are closely related to the Ab4 and Japanese strains.

Conclusion

The findings of the present study reveal several genomic deletions and in-frame mutations in the KyA genome as compared to those of RacL11 and Ab4. In vivo studies demonstrated that KyA-infected mice lost less than 10% of their body weight, rapidly regained weight, and all mice survived. By contrast, the RacL11 strain induced severe infection, as evident by significant weight loss, ruffled fur, severe inflammation of the lungs, and an increased mortality rate. Our previous results showed that the immunization of CBA mice with KyA protected animals from wt RacL11 challenge (14, 37, 38, 50). In vitro results demonstrated that KyA was able to grow rapidly in a variety of cell types, including mouse, rabbit, equine, and human cells. The attenuation of KyA in the mouse model and the ability of KyA to replicate to high titers make it an attractive vaccine candidate. We are planning to study EHV-1 KyA in horses. Our future study will give more insight into the suitability of KyA as a vaccine to combat EHV-1-mediated disease in the horse. Understanding the role in attenuation of genomic deletions and mutations may help in the development of a safe and effective EHV-1 vaccine.

Ethics Statement

All the animal experiments were done in accordance to the guidelines of the Louisiana State University Health Sciences Center at Shreveport Institutional Animal Care and Use Committee (IACUC). The LSUHSC-Shreveport OLAW assurance number is A3095-01. Ethical clearance for all animal studies was taken from the LSUHSC-Shreveport IACUC under protocols P-16-030 and P-15-038.

Author Contributions

AS performed most of the experiments, analyzed data, and drafted the manuscript. DO conceived of the study, analyzed data, and drafted the manuscript. SK conceived of the study, performed some of the experiments, analyzed data, and drafted the manuscript. All authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Mrs. Suzanne Zavecz for her excellent technical assistance. We are grateful to Elsevier for permission to use material from the paper of Telford et al. (23).

Funding

This work, including the efforts of Seong K. Kim and Dennis J. O’Callaghan, was funded by NIH grants P30GM110703 and AI-22001. This work, including the efforts of Seong K. Kim, was funded by USDA National Institute of Food and Agriculture (NIFA) (2013-67015-21311).

References

1. O’Callaghan DJ, Harty RN. In: Webster RG, Granoff A, editors. Equine Herpesvirus in Encyclopedia of Virology (Vol. 1), San Diego, CA: Academic Press, Harcourt Brace & Company (1994).

Google Scholar

2. O’Callaghan DJ, Osterrieder N. 2nd ed. In: Webster RG, Granoff A, editors. Equine Herpesvirus in Encyclopedia of Virology. London: Academic Press Ltd. Harcourt Brace & Company (1999).

Google Scholar

3. Carroll CL, Westbury HA. Isolation of equine herpesvirus 1 from the brain of a horse affected with paresis. Aust Vet J (1985) 62(10):345–6. doi:10.1111/j.1751-0813.1985.tb07660.x

CrossRef Full Text | Google Scholar

4. Allen GP, Bryans JT. Molecular epizootiology, pathogenesis, and prophylaxis of equine herpesvirus-1 infections. Prog Vet Microbiol Immunol (1986) 2:78–144.

Google Scholar

5. Patel JR, Heldens J. Equine herpesviruses 1 (EHV-1) and 4 (EHV-4) – epidemiology, disease and immunoprophylaxis: a brief review. Vet J (2005) 170(1):14–23. doi:10.1016/j.tvjl.2004.04.018

CrossRef Full Text | Google Scholar

6. Holz CL, Nelli RK, Wilson ME, Zarski LM, Azab W, Baumgardner R, et al. Viral genes and cellular markers associated with neurological complications during herpesvirus infections. J Gen Virol (2017) 98(6):1439–54. doi:10.1099/jgv.0.000773

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Stokol T, Yeo WM, Burnett D, DeAngelis N, Huang T, Osterrieder N, et al. Equid herpesvirus type 1 activates platelets. PLoS One (2015) 10(4):e0122640. doi:10.1371/journal.pone.0122640PONE-D-14-56878

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kydd JH, Smith KC, Hannant D, Livesay GJ, Mumford JA. Distribution of equid herpesvirus-1 (EHV-1) in respiratory tract associated lymphoid tissue: implications for cellular immunity. Equine Vet J (1994) 26(6):470–3. doi:10.1111/j.2042-3306.1994.tb04051.x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Wilson WD. Equine herpesvirus 1 myeloencephalopathy. Vet Clin North Am Equine Pract (1997) 13(1):53–72. doi:10.1016/S0749-0739(17)30255-9

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Smith KC, Mumford JA, Lakhani K. A comparison of equid herpesvirus-1 (EHV-1) vascular lesions in the early versus late pregnant equine uterus. J Comp Pathol (1996) 114(3):231–47. doi:10.1016/S0021-9975(96)80045-4

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Reczko E, Mayr A. On the fine structure of a virus of the herpes group isolated from horses (short report). Arch Gesamte Virusforsch (1963) 13:591–3. doi:10.1007/BF01267802

CrossRef Full Text | Google Scholar

12. Smith PM, Zhang Y, Grafton WD, Jennings SR, O’Callaghan DJ. Severe murine lung immunopathology elicited by the pathogenic equine herpesvirus 1 strain RacL11 correlates with early production of macrophage inflammatory proteins 1alpha, 1beta, and 2 and tumor necrosis factor alpha. J Virol (2000) 74(21):10034–40. doi:10.1128/JVI.74.21.10034-10040.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Zhang Y, Smith PM, Jennings SR, O’Callaghan DJ. Quantitation of virus-specific classes of antibodies following immunization of mice with attenuated equine herpesvirus 1 and viral glycoprotein D. Virology (2000) 268(2):482–92. doi:10.1006/viro.2000.0197S0042-6822(00)90197-1

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Kim SK, Shakya AK, O’Callaghan DJ. Immunization with attenuated equine herpesvirus 1 strain KyA induces innate immune responses that protect mice from lethal challenge. J Virol (2016) 90(18):8090–104. doi:10.1128/JVI.00986-16JVI.00986-16

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Edington N, Bridges CG, Patel JR. Endothelial cell infection and thrombosis in paralysis caused by equid herpesvirus-1: equine stroke. Arch Virol (1986) 90(1–2):111–24. doi:10.1007/BF01314149

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Marshall KR, Sun Y, Brown SM, Field HJ. An equine herpesvirus-1 gene 71 deletant is attenuated and elicits a protective immune response in mice. Virology (1997) 231(1):20–7. doi:10.1006/viro.1997.8483

PubMed Abstract | CrossRef Full Text | Google Scholar

17. O’Callaghan DJ, Cheevers WP, Gentry GA, Randall CC. Kinetics of cellular and viral DNA synthesis in equine abortion (herpes) virus infection of L-M cells. Virology (1968) 36(1):104–14. doi:10.1016/0042-6822(68)90120-7

CrossRef Full Text | Google Scholar

18. O’Callaghan DJ, Hyde JM, Gentry GA, Randall CC. Kinetics of viral deoxyribonucleic acid, protein, and infectious particle production and alterations in host macromolecular syntheses in equine abortion (herpes) virus-infected cells. J Virol (1968) 2(8):793–804.

Google Scholar

19. Colle CF III, Tarbet EB, Grafton WD, Jennings SR, O’Callaghan DJ. Equine herpesvirus-1 strain KyA, a candidate vaccine strain, reduces viral titers in mice challenged with a pathogenic strain, RacL. Virus Res (1996) 43(2):111–24. doi:10.1016/0168-1702(96)01324-X

CrossRef Full Text | Google Scholar

20. Matsumura T, O’Callaghan D, Kondo T, Kamada M. Lack of virulence of the murine fibroblast adapted strain, Kentucky A (KyA), of equine herpesvirus type 1 (EHV-1) in young horses. Vet Microbiol (1996) 48(3–4):353–65. doi:10.1016/0378-1135(09)59999-3

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Henry BE, Robinson RA, Dauenhauer SA, Atherton SS, Hayward GS, O’Callaghan DJ. Structure of the genome of equine herpesvirus type 1. Virology (1981) 115(1):97–114. doi:10.1016/0042-6822(81)90092-1

CrossRef Full Text | Google Scholar

22. Ruyechan WT, Dauenhauer SA, O’Callaghan DJ. Electron microscopic study of equine herpesvirus type 1 DNA. J Virol (1982) 42(1):297–300.

PubMed Abstract | Google Scholar

23. Telford EA, Watson MS, McBride K, Davison AJ. The DNA sequence of equine herpesvirus-1. Virology (1992) 189(1):304–16. doi:10.1016/0042-6822(92)90706-U

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Davison AJ, Eberle R, Ehlers B, Hayward GS, McGeoch DJ, Minson AC, et al. The order Herpesvirales. Arch Virol (2009) 154(1):171–7. doi:10.1007/s00705-008-0278-4

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Harty RN, O’Callaghan DJ. An early gene maps within and is 3’ coterminal with the immediate-early gene of equine herpesvirus 1. J Virol (1991) 65(7):3829–38.

PubMed Abstract | Google Scholar

26. Kim SK, Ahn BC, Albrecht RA, O’Callaghan DJ. The unique IR2 protein of equine herpesvirus 1 negatively regulates viral gene expression. J Virol (2006) 80(10):5041–9. doi:10.1128/JVI.80.10.5041-5049.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Holden VR, Harty RN, Yalamanchili RR, O’Callaghan DJ. The IR3 gene of equine herpesvirus type 1: a unique gene regulated by sequences within the intron of the immediate-early gene. DNA Seq (1992) 3(3):143–52. doi:10.3109/10425179209034010

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Ahn BC, Breitenbach JE, Kim SK, O’Callaghan DJ. The equine herpesvirus-1 IR3 gene that lies antisense to the sole immediate-early (IE) gene is trans-activated by the IE protein, and is poorly expressed to a protein. Virology (2007) 363(1):15–25. doi:10.1016/j.virol.2007.01.024

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Ahn BC, Zhang Y, O’Callaghan DJ. The equine herpesvirus-1 (EHV-1) IR3 transcript downregulates expression of the IE gene and the absence of IR3 gene expression alters EHV-1 biological properties and virulence. Virology (2010) 402(2):327–37. doi:10.1016/j.virol.2010.03.051S0042-6822(10)00234-5

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Yalamanchili RR, O’Callaghan DJ. Sequence and organization of the genomic termini of equine herpesvirus type 1. Virus Res (1990) 15(2):149–61. doi:10.1016/0168-1702(90)90005-V

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Yalamanchili RR, Raengsakulrach B, O’Callaghan DJ. Equine herpesvirus 1 sequence near the left terminus codes for two open reading frames. Virus Res (1991) 18(2–3):109–16. doi:10.1016/0168-1702(91)90012-K

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Colle CF III, Flowers CC, O’Callaghan DJ. Open reading frames encoding a protein kinase, homolog of glycoprotein gX of pseudorabies virus, and a novel glycoprotein map within the unique short segment of equine herpesvirus type 1. Virology (1992) 188(2):545–57. doi:10.1016/0042-6822(92)90509-N

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Flowers CC, O’Callaghan DJ. The equine herpesvirus type 1 (EHV-1) homolog of herpes simplex virus type 1 US9 and the nature of a major deletion within the unique short segment of the EHV-1 KyA strain genome. Virology (1992) 190(1):307–15. doi:10.1016/0042-6822(92)91217-I

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Matsumura T, Smith RH, O’Callaghan DJ. DNA sequence and transcriptional analyses of the region of the equine herpesvirus type 1 Kentucky A strain genome encoding glycoprotein C. Virology (1993) 193(2):910–23. doi:10.1006/viro.1993.1200

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Bowles DE, Holden VR, Zhao Y, O’Callaghan DJ. The ICP0 protein of equine herpesvirus 1 is an early protein that independently transactivates expression of all classes of viral promoters. J Virol (1997) 71(7):4904–14.

PubMed Abstract | Google Scholar

36. Bowles DE, Kim SK, O’Callaghan DJ. Characterization of the trans-activation properties of equine herpesvirus 1 EICP0 protein. J Virol (2000) 74(3):1200–8. doi:10.1128/JVI.74.3.1200-1208.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Matsumura T, Kondo T, Sugita S, Damiani AM, O’Callaghan DJ, Imagawa H. An equine herpesvirus type 1 recombinant with a deletion in the gE and gI genes is avirulent in young horses. Virology (1998) 242(1):68–79. doi:10.1006/viro.1997.8984

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Frampton AR Jr, Smith PM, Zhang Y, Matsumura T, Osterrieder N, O’Callaghan DJ. Contribution of gene products encoded within the unique short segment of equine herpesvirus 1 to virulence in a murine model. Virus Res (2002) 90(1–2):287–301. doi:10.1016/S0168-1702(02)00245-9

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Soboll Hussey G, Hussey SB, Wagner B, Horohov DW, Van de Walle GR, Osterrieder N, et al. Evaluation of immune responses following infection of ponies with an EHV-1 ORF1/2 deletion mutant. Vet Res (2011) 42:23. doi:10.1186/1297-9716-42-231297-9716-42-23

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Nugent J, Birch-Machin I, Smith KC, Mumford JA, Swann Z, Newton JR, et al. 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 (2006) 80(8):4047–60. doi:10.1128/JVI.80.8.4047-4060.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger aatasets. Mol Biol Evol (2016) 33(7):1870–4. doi:10.1093/molbev/msw054msw054

CrossRef Full Text | Google Scholar

42. Brown JC. High G+C content of herpes simplex virus DNA: proposed role in protection against retrotransposon insertion. Open Biochem J (2007) 1:33–42. doi:10.2174/1874091X00701010033

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Roizman B. The function of herpes simplex virus genes: a primer for genetic engineering of novel vectors. Proc Natl Acad Sci U S A (1996) 93:11307–12. doi:10.1073/pnas.93.21.11307

CrossRef Full Text | Google Scholar

44. Roizman B, Knipe D, Whitley R. Herpes simplex viruses. In: Knipe D, Howley P, editors. Fields Virology. Philadelphia: Lippincott, Williams, and Wilkins (2013). p. 1823–97.

Google Scholar

45. Ahn BC, Kim S, Zhang Y, Charvat RA, O’Callaghan DJ. The early UL3 gene of equine herpesvirus-1 encodes a tegument protein not essential for replication or virulence in the mouse. Virology (2011) 420(1):20–31. doi:10.1016/j.virol.2011.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Schimmer C, Neubauer A. The equine herpesvirus 1 UL11 gene product localizes to the trans-Golgi network and is involved in cell-to-cell spread. Virology (2003) 308(1):23–36. doi:10.1016/S0042-6822(02)00060-0

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Guggemoos S, Just FT, Neubauer A. The equine herpesvirus 1 UL20 product interacts with glycoprotein K and promotes egress of mature particles. J Virol (2006) 80(1):95–107. doi:10.1128/JVI.80.1.95-107.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Kim S, Ahn BC, O’Callaghan DJ, Kim SK. The early UL31 gene of equine herpesvirus 1 encodes a single-stranded DNA-binding protein that has a nuclear localization signal sequence at the C-terminus. Virology (2012) 432(2):306–15. doi:10.1016/j.virol.2012.05.031S0042-6822(12)00282-6

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Neubauer A, Rudolph J, Brandmuller C, Just FT, Osterrieder N. The equine herpesvirus 1 UL34 gene product is involved in an early step in virus egress and can be efficiently replaced by a UL34-GFP fusion protein. Virology (2002) 300(2):189–204. doi:10.1006/viro.2002.1488

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Smith PM, Kahan SM, Rorex CB, von Einem J, Osterrieder N, O’Callaghan DJ. Expression of the full-length form of gp2 of equine herpesvirus 1 (EHV-1) completely restores respiratory virulence to the attenuated EHV-1 strain KyA in CBA mice. J Virol (2005) 79(8):5105–15. doi:10.1128/JVI.79.8.5105-5115.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

51. van Woensel PA, Goovaerts D, Markx D, Visser N. A mouse model for testing the pathogenicity of equine herpes virus-1 strains. J Virol Methods (1995) 54(1):39–49. doi:10.1016/0166-0934(95)00024-O

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Frampton AR, Smith PM, Zhang Y, Grafton WD, Matusmura T, Osterrieder N, et al. Meningoencephalitis in mice infected with an equine herpesvirus 1strain KyA recombinant expressing glycoprotein I and glycoprotein E. Virus Genes (2004) 29:9–17. doi:10.1023/B:VIRU.0000032785.19420.14

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Gibson JS, Slater JD, Field HJ. The pathogenicity of Ab4p, the sequenced strain of equine herpesvirus-1, in specific pathogen-free foals. Virology (1992) 189(1):317–9. doi:10.1016/0042-6822(92)90707-V

PubMed Abstract | CrossRef Full Text | Google Scholar