| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Journal of Virology, May 2007, p. 5058-5065, Vol. 81, No. 10
0022-538X/07/$08.00+0 doi:10.1128/JVI.00146-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Laboratory of Genome Science, Tokyo University of Marine Science and Technology, Konan 4-5-7 Minato, 108-8477 Tokyo, Japan,1 Laboratory of Comparative Genomics, Nara Institute of Science and Technology, Keihanna Science City, 8916-5, Takayama-cho, Ikoma, Nara 630-0192, Japan,2 Institute of Microbiology, Department of Clinical Microbiology, The Hebrew University-Hadassah Medical School, Ein Karen, Jerusalem, Israel,3 Kimron Veterinary Institute, Bet Dagan, Israel,4 MRC Virology Unit, Church St., Glasgow G11 5JR, United Kingdom,5 Department of Medicine and Epidemiology, School of Veterinary Medicine, University of California, Davis, California 956166
Received 22 January 2007/ Accepted 20 February 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Koi herpesvirus (KHV) was first identified in 1998 as the cause of mass mortality among juvenile and adult koi and among common carp cultured in Israel, the United States, and Germany (5, 24, 25). Considerably more information has since been obtained on the properties of this virus, including its host range, the effects of water temperature on disease outcome, development of detection methods, and novel attempts at control (1, 4, 15-17, 20, 26, 27, 34, 37-39). The virus has now been found associated with mass mortality events on most continents, including countries throughout Europe, the United States, Japan, Indonesia, South Africa, Thailand, Taiwan, China, and Malaysia (22, 29, 42, 45). Outbreaks continue to occur among cultured populations of koi in most of these areas and, alarmingly, epidemics have been observed among wild common carp populations in the United Kingdom, Japan and, most recently, several states in the United States (42; P. Bowser and J. Grizzle, unpublished data; K. Way, unpublished data; R. Hedrick, unpublished data).
Taxonomically, members of the family Herpesviridae consist of a linear, double-stranded DNA genome packaged within an icosahedral capsid of characteristic architecture that is surrounded by a proteinaceous tegument layer and finally by a host-derived envelope containing virus glycoproteins (11). Phylogenetically, herpesviruses (HVs) fall into three major clades that are barely related at the sequence level: those that infect mammals, birds, and reptiles (the mammalian HV clade), those that infect fish and amphibians (the fish HV clade), and a single virus that infects an invertebrate, the oyster (the bivalve HV clade). A substantial number of viruses in the mammalian HV clade have been sequenced, but only three in the fish HV clade (ictalurid HV 1, IcHV-1, also known as channel catfish virus, and two ranid HVs [7, 10]) and one in the bivalve HV clade (ostreid HV 1, OsHV-1, also known as oyster HV [12]). It has been proposed that these clades should be classified as three families encompassed by an order, the Herpesvirales (8, 12). The families would be the revised Herpesviridae, the Alloherpesviridae, and the Malacoherpesviridae, respectively.
Early observations demonstrated the clear developmental and morphological affinities of KHV to the HVs (24), but subsequent findings prompted others to suggest that the virus should not be considered a member of the Herpesviridae (26, 38, 39). Key factors in this argument included the estimated genome size of KHV (277 kbp), which is greater than that observed previously among HVs (125 to 245 kbp), and the failure to demonstrate convincing genetic relationships between KHV and recognized HVs from limited sequence analyses. However, comparisons conducted on a more appropriate basis made it clear that KHV is related to fish HVs (46). KHV is related closely to cyprinid HVs 1 and 2 (CyHV-1 and CyHV-2), the agents associated with carp pox and hematopoietic necrosis in goldfish, respectively (28, 41), and distantly to IcHV-1 and ranid HV 1 (RaHV-1, also known as Lucké tumor herpesvirus of frog) (46). Therefore, KHV has been proposed formally as a member of the Alloherpesviridae under the species name Cyprinid herpesvirus 3.
In the present study, we describe the complete genome sequences of three KHV strains isolated from Japan (strain J), the United States (strain U), and Israel (strain I). We confirm that KHV shares significant similarities to fish HVs and belongs as a new species in the family Herpesviridae. The comparative information gained from the three strains also provides insights into epidemiologic features of an emerging disease threatening koi and carp worldwide.
| MATERIALS AND METHODS |
|---|
|
|
|---|
DNA sequencing. The complete genome sequences were determined commercially (Hitachi Instruments Service Co., Ltd., Tokyo, Japan) by standard shotgun sequencing. Briefly, fragments in the 2- to 5-kbp size range generated from purified viral DNA by passage through a syringe needle were cloned into pUC118 to construct a plasmid library. The KHV DNA inserts were sequenced from both ends of the plasmids by using the M13 forward and reverse primers with a BigDye Terminator v3.1 cycle sequencing kit and ABI 3700 and 3730xl DNA analyzers.
DNA sequence analysis. The DNA sequences were assembled by using Phrap (14), employing the quality files and default settings to produce a consensus sequence that was edited by using Consed (18). The final consensus sequence in each case represented an average eightfold redundancy at each base position. Gap closure was achieved by primer walking of gap-spanning clones and sequencing of PCR products. The sequence traces were checked carefully at positions where coding regions were fragmented in order to rule out errors. The genome termini in U were located precisely by PCR amplification of virion DNA that had been flush ended and ligated to a partially double-stranded adaptor oligonucleotide as described previously (13), using a primer matching part of the adaptor plus 5'-CGCAGTAGGCCTTGACCAGCA-3' for the left terminus or 5'-AGGGTCTGAACATGGCTTAGG-3' for the right terminus. The locations of the termini determined from sequencing of 12 plasmid clones of each product were assumed to be identical in J and I.
The GCG suite (Accelrys) and other programs (19, 33) were used for standard analyses of the sequences. ATG-initiated open reading frames (ORFs) of >50 codons were evaluated initially for protein coding potential by using GAMBLER, which semiautomates analysis of the output from genome assembler software, assigns ORFs automatically, and carries out homology searching (40). Additional homology searches were carried out by using the FastA (36) and BLAST (2) tools, and alignments were made by using CLUSTAL W (44). The preliminary gene map was refined by using GCG Codonpreference, which measures codon preference and third position codon G+C bias on the basis of statistics generated from standard genes. Since KHV happens to be strongly biased in these measures, the analysis proved useful in revising the preliminary map, as a result of which most ORFs were confirmed as likely to be protein coding and several ambiguities were clarified. Four large KHV ORFs with a high probability of encoding proteins were used as the standards in Codonpreference: ORF79, ORF92, ORF107, and ORF62. Additional criteria routinely used in HV genome analysis were used, in particular the lack of extensive overlap between coding regions and the presence of polyadenylation signals downstream from ORFs or families of similarly oriented ORFs. Splicing is relatively rare in HVs, and predictions for KHV based on the presence of appropriately located splice donor and acceptor signals were kept to a minimum.
The availability of an incomplete sequence database for the majority of the related CyHV-1 genome (T. Waltzek and R. Hedrick, unpublished data) enabled most predicted protein-coding regions in KHV (including splice sites) to be confirmed by homology and helped clarify many remaining difficulties. A few aspects of the interpretation were also strengthened by comparisons with homologous ORFs in the more distantly related IcHV-1 genome (7).
Analysis of mutated regions. Regions containing frameshift mutations in one or another strain were amplified by standard PCR from U at passages 2 and 12 (the former from infected cell DNA and latter from infected cell DNA and purified virion DNA). In addition, data were obtained for a fourth KHV strain (CO5-53-3) at passage 0 (i.e., directly from infected fish tissue) and 2. The CO5-53-3 strain was obtained in 2005 from a moribund wild common carp during a mortality episode in the lower San Joacquin River, California. Standard PCRs were carried out, and the amplicons were sequenced directly using the PCR primers.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
The genomes are highly similar to each other at the sequence level, with U and I more closely related to each other than either is to J. For example, in respect of single nucleotide substitutions (not counting duplicates in the terminal repeat), J differs from U and I at 181 of the 217 loci that are not conserved in all three strains; that is, there is a nucleotide difference every 1.5 kbp on average. Of the remaining 36 nonconserved loci, I differs from J and U at 32 loci, and U differs from J and I at 4 loci. These relationships imply a history in which an ancestral KHV strain gave rise initially to two branches: the J lineage that led eventually to J and the U/I lineage, which subsequently split into the branches leading to U and I. This conclusion is also consistent with the pattern of differences due to insertions and deletions.
Genetic content. The close relationship among the three strains resulted in the same predicted gene map for each, with the exception of a number of fragmented ORFs discussed below. Figure 1 shows the predicted layout of protein-coding genes in an ancestral (wild-type) KHV as envisaged prior to gene fragmentation events. The genetic complement is 156, with 8 duplicated in the terminal repeat, yielding a total of 164 in the genome. Table 1 lists the features of KHV genes for which information was derived from bioinformatic analysis.
|
|
The KHV genome contains 15 genes (colored red in Fig. 1) that have clear homologs in IcHV-1. The data for IcHV-1 and other HVs (7-9) indicate that these genes encode proteins involved in capsid morphogenesis (ORF92, ORF72, and ORF78, encoding two structural components of the capsid shell and a candidate protease involved in capsid maturation), nucleotide metabolism (ORF19 and ORF123 encoding deoxynucleoside kinase and deoxyuridine triphosphatase, respectively), DNA replication (ORF79, ORF71, and ORF46, encoding DNA polymerase, helicase, and a candidate primase), and DNA packaging (ORF33, which encodes the putative ATPase subunit of terminase; the three exons that comprise this coding region are shown in Fig. 1 as connected by introns). In addition, the conserved genes include one encoding a large membrane glycoprotein (ORF99) and five encoding proteins with unknown functions (ORF47, ORF61, ORF80, ORF90, and ORF56). This modest number of convincingly conserved genes highlights the substantial evolutionary distance between KHV and IcHV-1 (26, 46). When the ranid HV sequences (10) are taken into account, only 13 genes are convincingly conserved: the 15 genes conserved between KHV and IcHV-1 as described above, less ORF19 and ORF123. These findings indicate that the fish HV clade is considerably more divergent overall than the mammalian HV clade, in which 43 genes have been inherited from a common ancestor (8). Nonetheless, KHV clearly represents a new HV species, as concluded previously (23, 46). It is likely that the number of conserved genes in the fish HV clade will increase as data for other species facilitate more sensitive comparisons, but not to the number observed in the mammalian HV clade. As is the case with distantly related viruses in the mammalian HV clade (8), conserved genes are located centrally in the KHV and IcHV-1 genomes, and their arrangement is not conserved.
The remaining 141 KHV genes are in the nonconserved set. These are colored pink in Fig. 1, except for five families of related genes. The RING family encodes four proteins containing a zinc ion-coordinating motif known as the RING finger. This motif is ubiquitous among the HVs and is encoded by representatives of the three clades. The tumor necrosis factor receptor (TNFR) family encodes two versions of a secreted form of TNFR, which presumably have roles in immune evasion, as demonstrated with the poxviruses (43). The ORF2 and ORF22 families encode two and three proteins, respectively. The ORF25 family encodes six membrane glycoproteins, of which two are fragmented, one (ORF27) at different locations in U and I and one (ORF26) in all three strains.
The putative products of several KHV genes are convincingly related to enzymes in addition to those described above. These include the large and small subunits of ribonucleotide reductase (ORF141 and ORF23, respectively), thymidine kinase (ORF55), thymidylate kinase (ORF140), uracil-DNA glycosylase (ORF98), serine protease (ORF94), and serine-threonine protein kinase (ORF104). KHV also encodes proteins related to G protein-coupled receptors (ORF16), eukaryotic DUF614 proteins (ORF31), a family of iridovirus proteins (ORF32), nucleoside transporters (ORF64), a cellular protein of unknown function (ORF114), and a poxvirus B22R protein (ORF139) that is likely to be involved in immune evasion (30). Like some members of the mammalian HV clade (6), KHV encodes a protein that is clearly related to interleukin-10 (ORF134), which may modulate host immune responses. Other KHV proteins bear limited similarities to bacterial NAD-dependent epimerase/dehydratases (ORF28), protein kinases (ORF48), zinc-binding proteins (ORF54), cysteine proteases (ORF62), myosin (ORF68), and nucleic acid-binding proteins (ORF112). It is apparent from the conserved and nonconserved gene sets that KHV evolution has been characterized by gene capture from the cell or other viruses.
Strain evolution. As mentioned above, the three KHV strains differ from each other at only a small number of loci. A few of the insertions and deletions represent mutations (usually frameshifts) in one or more strains that disrupt coding regions. It should be noted that the identification of affected ORFs depends upon the accuracy of the predicted gene set and, moreover, that genes mutated in all strains may only be identified in certain circumstances. This is possible when such genes are members of a family or are related to genes in other organisms or where the encoded proteins contain characteristic features, such as those of a membrane protein with a signal sequence and transmembrane region.
Despite these analytical limitations, 15 ORFs (10% of the complement) appear to be mutated in one or more strains (Table 1). Of these, the unmutated forms of 11 ORFs encode proteins with features implying that their assignment as protein-coding regions at least is correct. It is notable that the majority of mutated ORFs encode membrane glycoproteins and that none of the conserved genes is affected. Nine genes (listed as "broken" in Table 1) are probably rendered nonfunctional by frameshifts located centrally in the coding regions, and the wild-type gene is identified readily. Six genes (listed as "frameshifted") might retain function since the mutations are close (sometimes very close) to the 3' ends of the ORFs, making it more difficult to identify the wild-type gene (see the assumptions explained in the footnotes to Tables 1 and 2). Nonetheless, even though identification of the precise number of mutated genes in each genome is problematic, all three sequenced strains are evidently multiple mutants derived from a wild-type ancestor.
|
The sequence comparisons indicate that the three sequenced KHV strains arose via the loss of genetic functions, as evidenced by frameshifted coding regions, with at least some of the cognate mutations having occurred in vivo. Elsewhere among large DNA viruses, fragmented genes have been documented extensively in the Poxviridae, particularly the Orthopoxvirus genus, where up to 16% of genes in a single virus may be fragmented (21). These genes are invariably not required for virus growth in cell culture and are therefore presumed to contribute to some aspect of growth in the host. The simplest interpretation is that some orthopoxviruses have become associated with their hosts relatively recently and have diverged from an ancestral virus by losing certain functions either because these functions are not required or because they reduce fitness. Gene fragmentation was also observed in the bivalve HV, OsHV-1, which had not been passaged in cell culture, suggesting that this might have contributed to the striking pathogenicity (and perhaps increased host range) of this emerging virus in farmed shellfish (12). The apparent loss of KHV gene functions, particular among those encoding membrane glycoproteins that may be associated with host specificity and critical to virulence, presents a provocative parallel. Nonetheless, a role for gene loss in the emergence of KHV is currently speculative.
Origins of KHV. The epidemiologic and pathogenic features of KHV-associated disease are new, and we doubt that they were overlooked previously (23, 46). The earliest known archival evidence indicates the presence of KHV among wild common carp in the United Kingdom as early as 1996 (K. Way, unpublished data), preceding observations of the disease in koi in Germany that were first recorded in 1997 (5). The active and often unregulated movements of large numbers of koi have contributed to a rapid spread of the virus presumably from these origins (22). The first cases of the disease in the eastern United States (strain U) occurred in 1998 following a koi show in New York that involved fish from Israel (24), a finding consistent with the high degree of similarity between U and I. Causal links for the origin of KHV among common carp in Japan are less defined (42).
Despite these observations, the origins of KHV remain obscure. It is possible that it has derived from an innocuous virus of C. carpio or another of the more than 2,000 species of cyprinid fishes, escaping from a long-standing equilibrium with its host and adapting via increased virulence to the conditions of aquaculture, where susceptible hosts are in abundant supply, constantly renewed, and transported worldwide. In this respect, findings with KHV have relevance beyond carp aquaculture, since similar phenomena have been observed with herpesviruses of other intensively cultured animals, including Marek's disease virus in chickens, where the evolution of increasingly virulent strains is exacerbated by vaccination (32), and perhaps OsHV-1 in oysters (12).
The sequence comparisons prompt the hypothesis that intensive culture of common carp and koi, combined with large-scale movements of live koi, may have favored transmission of genetically deficient KHV strains of enhanced virulence. Testing of this hypothesis will include further extensive comparisons of KHV isolates that have not been passaged in cell culture in order to understand the extent and timescale of gene loss and specific mutagenesis studies in order to assess any contribution of gene loss to pathogenicity. Whether or not further adaptation might occur, more active control measures, including listing of the KHV disease by international organizations charged with disease control, are probably needed to reduce future economic and ecologic impacts of this important viral pathogen.
Nucleotide sequence accession numbers. The genome sequences for KHV strains J, U, and I were submitted to GenBank/DDJB under accession numbers AP008984, DQ657948, and DQ177346, respectively.
| ACKNOWLEDGMENTS |
|---|
We are grateful to Charles Cunningham for technical support in preparing and sequencing PCR products and to Duncan McGeoch for helpful comments on the manuscript.
| FOOTNOTES |
|---|
Published ahead of print on 28 February 2007. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Mol. Cell. Biol. | Microbiol. Mol. Biol. Rev. |
|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|