- Research article
- Open Access
Low genetic diversity among historical and contemporary clinical isolates of felid herpesvirus 1
BMC Genomics volume 17, Article number: 704 (2016)
Felid herpesvirus 1 (FHV-1) causes upper respiratory tract diseases in cats worldwide, including nasal and ocular discharge, conjunctivitis and oral ulceration. The nature and severity of disease can vary between clinical cases. Genetic determinants of virulence are likely to contribute to differences in the in vivo phenotype of FHV-1 isolates, but to date there have been limited studies investigating FHV-1 genetic diversity. This study used next generation sequencing to compare the genomes of contemporary Australian clinical isolates of FHV-1, vaccine isolates and historical clinical isolates, including isolates that predated the introduction of live attenuated vaccines into Australia. Analysis of the genome sequences aimed to assess the level of genetic diversity, identify potential genetic markers that could influence the in vivo phenotype of the isolates and examine the sequences for evidence of recombination.
The full genome sequences of 26 isolates of FHV-1 were determined, including two vaccine isolates and 24 clinical isolates that were collected over a period of approximately 40 years. Analysis of the genome sequences revealed a remarkably low level of diversity (0.0–0.01 %) between the isolates. No potential genetic determinants of virulence were identified, but unique single nucleotide polymorphisms (SNPs) in the UL28 and UL44 genes were detected in the vaccine isolates that were not present in the clinical isolates. No evidence of FHV-1 recombination was detected using multiple methods of recombination detection, even though many of the isolates originated from cats housed in a shelter environment where high infective pressures were likely to exist. Evidence of displacement of dominant FHV-1 isolates with other (genetically distinct) FHV-1 isolates over time was observed amongst the isolates obtained from the shelter-housed animals.
The results show that FHV-1 genomes are highly conserved. The lack of recombination detected in the FHV-1 genomes suggests that the risk of attenuated vaccines recombining to generate virulent field viruses is lower than has been suggested for some other herpesviruses. The SNPs detected only in the vaccine isolates offer the potential to develop PCR-based methods of differentiating vaccine and clinical isolates of FHV-1 in order to facilitate future epidemiological studies.
Felid herpesvirus 1 (FHV-1), an alphaherpesvirus, causes upper respiratory tract disease in cats that is characterised by pyrexia, severe nasal and ocular discharge, conjunctivitis and oral ulceration. It is a common cause of necrotising upper respiratory tract infections (URI) in cats housed in shelters, with more than 50 % of cats shedding virus only one week after admission . The virus typically replicates in epithelial cells of the conjunctiva and the upper respiratory tract and infects local neurons, establishing latency in the trigeminal, pterygopalatine and cranial cervical ganglia . Although virus replication is generally limited to the upper respiratory tract and conjunctiva, viraemia has been detected during the acute phase of infection [3, 4] and a recent report has described a case of non-suppurative meningoencephalitis in cats, suggesting that FHV-1 has potential to cause more invasive disease . Vaccination has been used to help control disease caused by FHV-1 since the mid-1970s. Both attenuated and inactivated vaccines are in use currently, and while they can reduce the severity of disease, they do not prevent infection [6–9].
Early studies demonstrated that isolates of FHV-1 were antigenically indistinguishable using conventional serological techniques. Additionally, a high level of genomic homogeneity was identified using restriction endonuclease cleavage analyses [10, 11]. A subsequent restriction endonuclease cleavage analysis study of 78 Japanese isolates detected three FHV-1 genotypes, although these genotypes still appeared to be very similar . Since these early studies there has been only limited investigation of the extent of variation between vaccine and wild type strains of FHV-1.
This study aimed to use next generation sequencing to compare the genomes of contemporary Australian clinical isolates of FHV-1 to those of historical clinical isolates, including isolates that predate the introduction of live attenuated vaccines into Australia. Analysis of the genome sequences aimed to assess the level of genetic diversity, identify potential genetic markers that could influence the in vivo phenotype of the isolates and examine the sequences for evidence of recombination.
Genome sequence analysis of 26 FHV-1 isolates
The full genome sequences of a panel of 26 FHV-1 isolates were assembled and analysed. This included 24 Australian clinical field isolates spanning nearly 40 years, as well as isolates from two currently available vaccines. The genome sequence of the prototype C-27 FHV-1 strain (GenBank:NC_013590), a historical isolate from North America, was also included in the analysis . The genomes of the Feligen vaccine isolate of FHV-1 assembled by two different methods (de novo assembly and map to reference assembly) were identical, with the exception of the regions containing short sequence repeats. This provided a high level of confidence in genome assembly by mapping to a reference sequence and this method was used for all other FHV-1 isolates in this study. Genome sequence data have been deposited into GenBank under the accession numbers as listed in Table 1 [GenBank:KR296657 and GenBank:KR381779 - KR381803]. The size of the FHV-1 full genomes ranged from 133.5 kbp (isolate 3227–05) to 135.1 kbp (isolate 3225–05). Genome length variation occurred primarily as a result of differences in the large sequence repeat regions (IRS/TRS) and in the unique short region as a result of tandem repeat iterations.
Despite the long time span over which the isolates were obtained, alignment of the FHV-1 genomes (Fig. 1a) revealed that diversity between them only ranged from 0.0 to 0.01 %. While the limited geographical range of many of the Australian isolates could have influenced this, three of the genomes analysed in this study were from isolates of North American origin. The most genetically diverse isolates (135/68 from 3225/05 and 3230/05), differed genome-wide by only 83 single nucleotide polymorphisms (SNPs), while some isolates shared 100 % nucleotide identity (excluding large regions of tandem repeat reiterations). The identical viruses included 356/75 with 571/79, while the eleven 2004/2005 field isolates differed by 0 to 4 SNPs, and the three 2006 isolates differed by 1 or 2 SNPs. The two vaccines were identical (excluding large reiterative repeat regions) and so are both likely to have been derived from the widely used F2 strain. The C-27 strain, a field strain isolated contemporaneously with the F2 vaccine strain from the USA, varied from the F2-derived vaccine strains by only 5 SNPs genome-wide (excluding large reiterative repeat regions) [14, 15]. The non-synonymous nucleotide differences (SNPs and insertions/deletions, or indels), as indicated by predicted amino acid residue differences between isolates, are provided in Table 2.
Phylogenetic and recombination analyses of FHV-1 isolates
Phylogenetic analyses revealed that FHV-1 isolates fell into two main groups (Fig. 1b). Whilst some isolates with known geographical or temporal associations did cluster together (for example all 2004 and 2005 isolates), overall there was no clear relationship between viral genome sequences, the geographical location where the virus was isolated, the year of virus isolation or whether they predated the date of introduction of the attenuated FHV-1 vaccine into Australia in 1977. No clear association was detected between the genome sequence and the presence or nature of disease in the cat from which the virus was isolated.
The genome sequences of Australian isolates that predated the introduction of the FHV-1 vaccine into Australia (in 1977) had a high sequence identity with isolates obtained since then, and also with isolates originating from the USA. Overall, our data demonstrated that FHV-1 genomes are highly conserved, more so than any of the other herpesviruses analysed to date, with a similar study of equine herpesvirus 1 (EHV-1) isolates finding variation of between 5 and 173 SNPs between EHV-1 isolates . Previous studies have suggested greater differences between FHV-1 isolates, including re-arrangements in the UL44 (gC) gene , restriction endonuclease cleavage site differences in the UL17 gene , and differences in MluI cleavage patterns across the genome . However, all the isolates examined in our study had conserved gene arrangements, with only one SNP detected in UL44 (a non-synonymous substitution), and the predicted restriction endonuclease cleavage patterns generated from the sequence data varied mainly as a result of variations in short tandem repeat iterations.
This study included analyses of fourteen FHV-1 isolates from cats housed in an animal shelter between 2004 and 2006. The cluster of viruses from this shelter in 2004 and 2005 is of particular interest because of the apparent transience of a dominant virus in a shelter environment. The SNP profile of the viruses from 2004 to 2005 was distinct from those of all the other isolates sequenced, but this pattern was not detected in the 2006 isolates from the same shelter. The close phylogenetic relationship between the viruses isolated in 2004/2005, and between those collected in 2006, suggests endemic circulation in each of these periods, with clear evidence of displacement of the 2004/2005 viruses by those detected in 2006. This displacement may have occurred as a result of increased transmissibility of the newer viruses, or management changes, such as seasonal crowding , depopulation and repopulation of the shelter, or sampling bias (only samples free of feline calicivirus were selected for use in this study).
In vitro studies have reported a high rate of recombination in FHV-1 . However, no evidence of recombination was detected between the 27 FHV-1 isolates, irrespective of genomic region analysed or the program used (P > 0.05). It is possible that recombination events may not have been detected because of the low levels of diversity among the viruses in our study. It is also possible that previous in vitro studies may have overestimated the likelihood of FHV-1 recombination because of the use of high multiplicities of infection and simultaneous infection, as both these factors have been shown in previous studies to increase the rate of recombination in vitro in other herpesviruses . Recombination would be predicted in circumstances where the opportunity for co-infection is high and where the recombinant progeny have a survival advantage over the parent viruses and dominant field strains. While infective pressure in a shelter environment would be expected to allow concurrent infections with different FHV-1 isolates, there may be insufficient selective pressure on recombinants for them to emerge and persist in vivo in a natural (shelter) environment, particularly as the low levels of diversity among FHV-1 viruses may result in insufficient selective distinction for specific viruses to emerge in a population.
Attenuated herpesvirus vaccines are used not only against FHV-1 in cats, but also against Marek’s disease virus and infectious laryngotracheitis virus (ILTV) in poultry, bovine herpesvirus 1 in cattle, pseudorabies virus in swine, and EHV-1 and equine herpes virus 4 (EHV-4) in horses. Previous studies of ILTV in poultry have demonstrated that attenuated vaccine strains can participate in natural recombination events that lead to the generation of more virulent progeny . The risk of recombination involving other attenuated veterinary herpesvirus vaccines has not been evaluated comprehensively, but the risk appears to higher in some viruses (such as EHV-4) and lower in others (such as EHV-1) . The results from this study suggest that the risk for FHV-1 vaccine recombination may also be low. Although the FHV-1 genome sequencing data from this study confirms the findings of earlier studies that FHV-1 s are relatively homogeneous  this study did detect SNPs unique to F2-derived vaccine strains, in genes UL28 and UL44. A search of all FHV-1 sequences available on GenBank identified only nine sequences covering these regions of the FHV-1 genomes, all relating to the UL44 gene. None of the sequences contained the vaccine-specific SNP, however the sequences were mostly of unknown origin apart from two; a BAC of the C-27 strain [GenBank:GM036023] and strain C7301 [GenBank:D86616]. There is the potential to use the F2-derived vaccine unique SNPs to differentiate vaccine FHV-1 isolates from other circulating FHV-1 s, and thus facilitate future epidemiological studies in FHV-1, especially to investigate the epidemiology of attenuated FHV-1 vaccines. Sequencing of more geographically diverse isolates will help confirm whether these SNPs are suitable markers for such studies.
Viruses and genome sequencing
Ten historical and fourteen recent Australian FHV-1 isolates, as well as two live attenuated vaccines currently available in Australia (Table 1) were selected from our laboratory archive. Historical isolates were obtained from clinical samples submitted to our laboratory as a component of past diagnostic services. Recent samples were collected from shelter-housed cats in 2004 to 2006 with approval from the Faculty of Veterinary Science, Animal Ethics Committee at The University of Melbourne (reference number 0004055.1) and with the consent from the shelter owners. Viruses were propagated in cultures of Crandell Rees feline kidney (CRFK) cells  for six passages, during which three plaque purifications were performed. CRFK cells are derived from a normal, female, 12 weeks old cat and are commercially available (ATCC CCL-94). Isolates selected for sequencing dating from 2004 to 2006 were all from animals that were negative for feline calicivirus, as determined by RT-PCR . Viral nucleocapsid genomic DNA was purified and sequenced as described previously . Libraries were prepared using 50 ng of viral genomic DNA with the Illumina Nextera DNA library preparation kits according to manufacturer’s instructions, and loaded onto an Illumina MiSeq. Sequencing was carried out using a 300 cycle V2 SBS kit (Illumina) in paired-end format. Reads were trimmed to an error probability limit of 0.5 % and mapped against the prototype FHV-1 sequence, strain C-27 [GenBank:NC_013590], using medium-low sensitivity options in the bioinformatics package Geneious V6.1.7 . For comparison, a de novo assembly of the genome of the F2 strain of FHV-1 from the Feligen vaccine (Virbac) was also performed using the medium-low sensitivity settings in Geneious. The details of the sequencing metrics are provided in Additional file 1: Table S1.
Phylogenetic and recombination analyses
Alignments of the complete genome sequences of the FHV-1 isolates, excluding the terminal repeat regions (TRS), were prepared using the Multiple Alignment with Fast Fourier Transformation (MAFFT) version 7 plugin in Geneious . The F2 vaccine strain of FHV-1 was used as the reference sequence. The historical C-27 strain of FHV-1 was also included throughout the analyses. Phylogenetic analyses were performed on the complete genome sequences, excluding the TRS and large sequence gaps, using the neighbor-joining method in Geneious with the Jukes Cantor model of nucleotide substitution . One thousand bootstrap replicates were used to assess the significance of the phylogenetic tree topology.
To detect evidence of historical recombination events between the isolates, two programs were used, SplitsTree4  and RDP4 , as previously described . Short tandem repeat regions were identified in the aligned genomes prior to analysis using the Phobos plugin in Geneious V6.1.7. Default settings and score constraints for satellites with a maximum repeat unit length of 50 nucleotides were used. The identified repeat regions were removed from all viruses in the alignment as the methods used in this study do not allow the length of these regions to be accurately determined, and differences in the length of these regions could give rise to false positive recombination events being detected. The resultant whole genome alignments, as well as the separate UL, US and IRS genomic region alignments, were analysed for evidence of recombination. Splits network trees were generated with an uncorrected P characters transformation model, ignoring gapped sites. Other models were tested but yielded no significant topological differences. Statistical analyses of the recombination networks were performed using the Phi test as implemented in SplitsTree4 . In RDP4 six different methods were used to assess the sequences for recombination breakpoints, MaxChi, Bootscan, SiScan, Chimaera, GENECONV and RDP, with default RDP4 settings used throughout.
Felid herpesvirus 1
Infectious laryngotracheitis virus
Internal/external repeat short
Single nucleotide polymorphism
Crandell Rees feline kidney cells
Polymerase chain reaction
Recombination detection program version 4
Pedersen NC, Sato R, Foley JE, Poland AM. Common virus infections in cats, before and after being placed in shelters, with emphasis on feline enteric coronavirus. J Feline Med Surg. 2004;6(2):83–8.
Townsend WM, Jacobi S, Tai SH, Kiupel M, Wise AG, Maes RK. Ocular and neural distribution of feline herpesvirus-1 during active and latent experimental infection in cats. BMC Vet Res. 2013;9:185.
Swenson CL, Gardner K, Arnoczky SP. Infectious feline herpesvirus detected in distant bone and tendon following mucosal inoculation of specific pathogen-free cats. Vet Microbiol. 2012;160(3–4):484–7.
Westermeyer HD, Thomasy SM, Kado-Fong H, Maggs DJ. Assessment of viremia associated with experimental primary feline herpesvirus infection or presumed herpetic recrudescence in cats. Am J Vet Res. 2009;70(1):99–104.
Hora AS, Tonietti PO, Guerra JM, Leme MC, Pena HF, Maiorka PC, Brandão PE. Felid herpesvirus 1 as a causative agent of severe nonsuppurative meningoencephalitis in a domestic cat. J Clin Microbiol. 2013;51(2):676.
Jas D, Aeberle C, Lacombe V, Guiot AL, Poulet H. Onset of immunity in kittens after vaccination with a non-adjuvanted vaccine against feline panleucopenia, feline calicivirus and feline herpesvirus. Vet J. 2009;182(1):86–93.
Lappin MR, Sebring RW, Porter M, Radecki SJ, Veir J. Effects of a single dose of an intranasal feline herpesvirus 1, calicivirus, and panleukopenia vaccine on clinical signs and virus shedding after challenge with virulent feline herpesvirus 1. J Feline Med Surg. 2006;8(3):158–63.
Slater E, York C. Comparative studies on parenteral and intranasal inoculation of an attenuated feline herpes virus. Dev Biol Stand. 1976;33:410–6.
Gaskell R, Dawson S, Radford A, Thiry E. Feline herpesvirus. Vet Res. 2007;38(2):337–54.
Grail A, Harbour DA. Restriction endonuclease analysis of DNA from isolates of feline herpesvirus type 1. Nihon Juigaku Zasshi. 1990;52(5):1007–13.
Horimoto T, Limcumpao JA, Xuan X, Ono M, Maeda K, Kawaguchi Y, Kai C, Takahashi E, Mikami T. Heterogeneity of feline herpesvirus type 1 strains. Arch Virol. 1992;126(1–4):283–92.
Maeda K, Kawaguchi Y, Ono M, Tajima T, Mikami T. Restriction endonuclease analysis of field isolates of feline herpesvirus type 1 and identification of heterogeneous regions. J Clin Microbiol. 1995;33(1):217–21.
Tai SH, Niikura M, Cheng HH, Kruger JM, Wise AG, Maes RK. Complete genomic sequence and an infectious BAC clone of feline herpesvirus-1 (FHV-1). Virology. 2010;401(2):215–27.
Crandell RA, Maurer FD. Isolation of a feline virus associated with intranuclear inclusion bodies. Proc Soc Exp Biol Med. 1958;97(3):487–90.
Bittle JL. Serologic relationship of new feline cytopathogenic viruses. Am J Vet Res. 1960;21:547.
Vaz PK, Horsington J, Hartley CA, Browning GF, Ficorilli NP, Studdert MJ, Gilkerson JR, Devlin JM. Evidence of widespread natural recombination among field isolates of equine herpesvirus 4 but not among field isolates of equine herpesvirus 1. J Gen Virol. 2015;97(3):747-55.
Hamano M, Maeda K, Mizukoshi F, Mochizuki M, Tohya Y, Akashi H, Kai K. Genetic rearrangements in the gC gene of the feline herpesvirus type 1. Virus Genes. 2004;28(1):55–60.
Hamano M, Maeda K, Kai K, Mochizuki M, Tohya Y, Akashi H. A novel genetic marker to differentiate feline herpesvirus type 1 field isolates. Vet Microbiol. 2005;106(3–4):195–200.
Sykes JE, Anderson GA, Studdert VP, Browning GF. Prevalence of feline Chlamydia psittaci and feline herpesvirus 1 in cats with upper respiratory tract disease. J Vet Intern Med. 1999;13(3):153–62.
Fujita K, Maeda K, Yokoyama N, Miyazawa T, Kai C, Mikami T. In vitro recombination of feline herpesvirus type 1. Arch Virol. 1998;143(1):25–34.
Thiry E, Meurens F, Muylkens B, McVoy M, Gogev S, Thiry J, Vanderplasschen A, Epstein A, Keil G, Schynts F. Recombination in alphaherpesviruses. Rev Med Virol. 2005;15(2):89–103.
Lee SW, Markham PF, Coppo MJ, Legione AR, Markham JF, Noormohammadi AH, Browning GF, Ficorilli N, Hartley CA, Devlin JM. Attenuated vaccines can recombine to form virulent field viruses. Science. 2012;337(6091):188.
Crandell RA, Fabricant CG, Nelson-Rees WA. Development, characterization, and viral susceptibility of a feline (Felis catus) renal cell line (CRFK). In Vitro. 1973;9(3):176–85.
Sykes JE, Allen JL, Studdert VP, Browning GF. Detection of feline calicivirus, feline herpesvirus 1 and Chlamydia psittaci mucosal swabs by multiplex RT-PCR/PCR. Vet Microbiol. 2001;81(2):95–108.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9.
Katoh K, Standley DM. MAFFT Multiple Sequence Alignment Software Version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.
Jukes TH. DDT: Maxwell’s Demon. Science. 1969;166(3901):44.
Huson DH. SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics. 1998;14(1):68–73.
Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evolution 2015;1(1). doi:10.1093/ve/vev003.
Bruen TC, Philippe H, Bryant D. A simple and robust statistical test for detecting the presence of recombination. Genetics. 2006;172(4):2665–81.
Studdert MJ, Martin MC. Virus diseases of the respiratory tract of cats. 1. Isolation of feline rhinotracheitis virus. Aust Vet J. 1970;46(3):99–104.
The authors thank Cynthia Brown for her excellent technical assistance.
The Australian Research Council (ARC DP130103991) and the Feline Health Research Fund supported this work. JMD is supported by ARC FT140101287.
Availability of data and materials
The genome sequence data of the 26 FHV-1 isolates have been submitted to GenBank and the accession numbers are KR296657, KR381779 - KR381803. Figure 1 is available from Figshare [DOI:10.4225/49/57BA69ECC410A].
PKV, NJ, JH, NF, MJS and JRG carried out the experiments, including sample collection, sample preparation, genome assembly and annotation and the molecular genetic studies discussed. CAH, GFB, JRG and JMD conceived of the study, designed the study and participated in the analysis of the data. All authors participated in the drafting of the manuscript, and read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Samples from Shelter B animals (2004–2006) were collected with approval from the Faculty of Veterinary Science, Animal Ethics Committee at The University of Melbourne (reference number 0004055.1) and with the consent from the shelter owners. Historical archived samples (1968–1983) were submitted to the laboratory as a component of institutional diagnostic services.
About this article
Cite this article
Vaz, P.K., Job, N., Horsington, J. et al. Low genetic diversity among historical and contemporary clinical isolates of felid herpesvirus 1. BMC Genomics 17, 704 (2016). https://doi.org/10.1186/s12864-016-3050-2
- Felid herpesvirus