- Open Access
Full-length genome sequence analysis of enzootic nasal tumor virus isolated from goats in China
© The Author(s). 2017
- Received: 23 November 2016
- Accepted: 3 July 2017
- Published: 26 July 2017
Enzootic nasal tumor virus (ENTV) is a betaretrovirus of sheep (ENTV-1) and goats (ENTV-2) associated with neoplastic transformation of epithelial cells of the ethmoid turbinate. Confirmation of the role of ENTV in the pathogenesis of enzootic nasal adenocarcinoma (ENA) has yet to be resolved due to the inability to culture the virus. Very little is known about the prevalence of this disease, particularly in China.
To evaluate the genetic diversity of ENTV-2 from Shaanxi province of China, the complete genome sequence of four isolates from Shaanxi province was determined by RT-PCR. These sequences were analyzed to evaluate their genetic relatedness with other small ruminant betaretroviruses. Phylogenetic analyses based on the gag gene and env gene were performed.
The ENTV-2-Shaanxi1 genome shared 97.0% sequence identity with ENTV-2-SC (accession number HM104174.1), and 89.6% sequence identity with the ENTV-2 sequences (accession number AY197548.1). ENTV-2 is closely related to the ENTV-1 and jaagsiekte retrovirus (JSRV). The main sequence differences between these viruses reside in LTR, two small regions of Gag, Orf-x, and the transmembrane (TM) region of Env. A stretch of 6 consecutive proline residues exists in VR1 of the ENTV-2-Shaanxi1 ~ 4 isolates. All the ENTV-2-Shaanxi isolates have the YXXM motif in the cytoplasmic tail of the Env. Phylogenetic analysis by nucleotide sequences showed that ENTV-2-Shaanxi1 ~ 4 isolates were closest related to two ENTV-2 isolates published in NCBI, especially with ENTV-2-SC strain.
This finding indicates that ENA most likely was introduced to Shaanxi province by the movement of contaminated goats from other areas in China. This study adds to understand the circulation, variation and distribution of ENTV-2, and may prove beneficial in future control or eradication programmes.
- Enzootic nasal tumor virus
- Genome sequence
- Phylogenetic analysis
Enzootic nasal adenocarcinoma (ENA) of sheep and goats and ovine pulmonary adenocarcinoma (OPA or jaagsiekte) are contagious diseases characterized by neoplastic transformation of secretory epithelial cells in the respiratory tract [1, 2]. OPA is a tumour derived from type II pneumocytes and Clara cells in the lung, whereas ENA arises from secretory epithelial cells of the ethmoid turbinate . The gross pathology and histology of ENA in sheep and goats appear to be identical [1, 4, 5]. Jaagsiekte sheep retrovirus (JSRV), an ovine betaretrovirus, is the causative agent of OPA . Enzootic nasal adenocarcinoma (ENA) is an economically important contagious tumour of the nasal mucosa of sheep and goats and is associated with enzootic nasal tumour virus (ENTV) [3–5]. ENTV belongs to the genus Betaretrovirus in the family Retroviridae. ENTV display characteristics of both D- and B-type viruses. The ENTV genome consists of a single, positive stranded RNA of about 7.5 kb and has a structure analogous to that of cellular mRNA, with 5′ and 3′ untranslated regions (UTR), a 7-methylguanosine cap and a polyadenylated 3′ end. Like other retroviral genomes, ENTV genome has a basic canonical structure, 5′-U5-Gag-Pro-Pol-Env-U3-3′, with four open reading frames (Orf) and flanking untranslated regions and terminal repeats (LTR) on either end. The gag (group antigen) gene encodes the structural proteins that make up the capsid (CA) and matrix (MA) layer or shell, as well as the nucleocapsid (NC) protein, which interacts directly with the genomic RNA. The pro gene encodes the viral protease (PR) responsible for proteolytic processing of viral proteins. The pol gene encodes reverse transcriptase (RT) and integrase (IN), the replicative enzymes required for reverse transcription and integration of the genome. The env gene encodes both SU (surface) and TM (transmembrane) subunits of the envelope glycoprotein . With the exception of Australia and New Zealand, ENA has been recorded worldwide wherever sheep and goats are farmed, with a prevalence of up to 10% in some areas [4, 6]. Very little is known about the prevalence of this disease due to the lack of an infectious molecular clone and the inability to culture the virus, and only two full-length sequences are available for ENTV-2 (ENTV-2 European strain: accession number AY197548.1, and Chinese ENTV-2-SC strain: accession number HM104174.1) [1, 7, 8].
Goat enzootic nasal tumors appeared enzootically in China in recent years. Clinically, the affected goats showed copious serous nasal discharge, then snuffle and progressively developed dyspnea, ocular protrusion and skull deformations, eventually death from suffocation. The epidemiology, clinical and pathohistological pattern of goat intranasal adenoma and adenocarcinoma in Shaanxi province were similar to those of goat enzootic intranasal tumors that had in Spain, France and other provinces of China . In order to understand the molecular evolution of ENTV-2, the full length genome sequence of four ENTV-2 derived from nasal fluid of ENA isolated from conventionally reared goats in China was determined and analyzed with other retroviruses. The results showed that the main sequence differences between these viruses reside in LTR, two small regions of Gag, Orf-x, and the transmembrane (TM) region of Env. Phylogenetic analysis revealed that these Shaanxi isolates were closely related to all the known ENTV-2 isolates, especially with ENTV-2-SC strain.
The information of the 4 ENTV-2
2015, Yangling, Shaanxi
2016, Baoji, Shaanxi
2016, Yangling, Shaanxi
2016, Lantian, Shaanxi
RNA extraction and RT-PCR
Viral RNA were extracted from the purified virus using the TIANamp Virus RNA Kit (Tiangen, Beijing, China) according to the manufacture’s protocol, and cDNA was synthesized by oligo(dT)-primed or random-primed reverse transcription, using the Moloney murine leukaemia virus M-MLV (H-) riboclone cDNA synthesis system (Promega, Madison, USA). The cDNA was used as a template for subsequent PCR analysis.
Primers for amplifying the complete genome of ENTV-2
Sequence and phylogenetic analysis
Profile of betaretroviruses isolates used for sequence analyses
Nucleotide sequence accession numbers
The information of goats with nasal tumors was in Table 1. Goat 1 and 3 were collected at different time points in the same farm in Yangling city, thereby making it possible to assess the level of ENTV-2 nucleotide divergence within a given farm. All of the sequences described in this manuscript are available in the GenBank Nucleotide Sequence Database under the accession numbers found in Table 1.
Clinical findings and gross pathology
Comparison of percentage nucleotide (Ntd) and amino acid (Aa) identity of the ENTV-2-Shaanxi1 consensus sequence with ENTV-2, ENTV-2-SC, ENTV-1-1NA1 and JSRV-USA
The pro open reading frame of ENTV-2 encodes a bifunctional protease/dUTPase of 308 amino acids (predicted molecular mass, 33 kDa) . As with other retroviruses, Pro is expressed as a polyprotein with Gag by a mechanism involving ribosomal frameshifting-the exact site of which and relative efficiency have yet to be determined . ENTV-2-Shaanxi1 Pro was greater than 97.7% identical at the amino acid level to ENTV-2-SC, and 94.8% identical at the amino acid level to ENTV-2 (accession number AY197548.1). Overall, the Pro region of ENTV-2-Shaanxi1 ~ 4 were highly homologous to other ENTV-2 and showed very little nucleotide or amino acid variability.
The pol protein of ENTV-2 encodes a 869-amino-acid (aa) polypeptide (predicted molecular mass, 99 kDa) . The Pol protein is synthesized as a polyprotein with Gag-Pro and is cleaved by Pro to produce reverse transcriptase (RT) and integrase (IN) after virus assembly. As with pro, the pol genes of ENTV-2-Shaanxi isolates were highly homologous to that of the ENTV-2-SC with 98.7% amino acid identity. Nearly every amino acid difference was the result of a synonymous substitution.
Phylogenetic analysis of small ruminant betaretroviruses
These features of ENA in Shaanxi province are similar to bighorn sheep (Ovis canadensis) sinus tumors [24, 25]. Clinically, domestic sheep and goats affected by ENA have abundant seromucinous nasal discharge. Similarities between ENA and bighorn sheep sinus tumors include the presence of seromucinous nasal discharge clinically, the gross finding of a soft white mass in the sinus cavity, and classification of some masses as adenocarcinoma [24, 26]. Additionally, the inflammatory nasal polyps often associated with ENA share characteristics with the hyperplastic masses described here for bighorn sheep, although in bighorn sheep the sinus tumors is a diffuse thickening of the sinus lining and not a discrete polypoid mass. Prominent differences between ENA and the sinus tumors described here are location (sinus tumors in bighorn sheep and nasal in domestic sheep and goats) and malignancy (predominantly, hyperplastic masses in bighorn sheep and neoplastic masses in domestic sheep and goats). Other prominent differences between the 2 entities include the papillary appearance and often grey-red color of ENA tumors not characteristic of bighorn sheep sinus tumors, as well as the prominent mesenchymal population histologically present in bighorn sheep sinus tumors but infrequently described for ENA .
This study represented four full-length ENTV-2 sequences from goat flocks in Shaanxi province of China, and made genomic analysis with other betaretrovirus genus. Understanding the genetic heterogeneity of ENTV-2 is important both as a tool for epidemiological studies and as means to clarify the origin and future evolution of ENTV-2 . The molecular sequence of this virus is closely related to sequences of JSRV and ENTV-1 [11, 27]. The ENTV-2-Shaanxi1 ~ 4 genome shared greater than 97.0% with ENTV-2-SC (accession number HM104174.1), and 89.6% sequence identity with the ENTV-2 sequence (accession number AY197548.1). ENTV-2-Shaanxi1 ~ 4 and ENTV-2-SC differ significantly to the ENTV-2 sequence (accession number AY197548.1). The isolation between different continents probably leads to the relatively large difference. ENTV-2 is closely related to ENTV-1 which is associated with enzootic adenocarcinoma of sheep, and to jaagsiekte retrovirus. The main sequence differences between these viruses resided in LTR, Orf-x, two small regions in Gag and the transmembrane (TM) region of Env. The LTR of ENTV-2-Shaanxi1 ~ 4 were similar to ENTV-2. The role of ENTV LTR in pathogenesis has yet to be determined. There was very little amino acid diversity in the Gag polyprotein except in the two domains defined as variable (VR1 and VR2) , suggesting that this part of the genomemay be able to withstand change. A stretch of 6 consecutive proline residues(Aa 122 ~ 127) exists in all ENTV-2. JSRV-USA contains a stretch of 7 consecutive proline residues at the region. To JSRV, the PPPPPPPS motif of the exogenous VR1 is neither necessary nor sufficient for particle release [12, 14]. Maybe the fuction of the proline-rich region in ENTV-2 is similar with JSRV. Orf-x is the most genetically diverse protein coding sequence of ENTV-2. Most of the amino acid differences were found in Orf-x, which in the corresponding ENTV-2-Shaanxi1 ~ 4 and ENTV-2-SC genome were 108 amino acids, but ENTV-2 sequence (accession number AY197548.1) is 166 amino acids. ENTV-2 in China potentially encode a truncated Orf-x protein compared to ENTV-2 (accession number AY197548.1). The same situation also appears in ENTV-1. Orf-x is the most genetically diverse protein coding sequence of ENTV-1 . One study showed that a complete Orf-x is not required for pathogenesis of JSRV . Truncation of Orf-x in JSRV did not alter the pathogenesis of the virus compared to wild type JSRV in experimental infection. So the large genetic changes of Orf-x in ENTV-2-Shaanxi may indicate that Orf-x does not play a significant role in in the pathogenesis of ENA. The sequence differences of Orf-x between different strains may be due to the regional differences. ENTV-2 is also more like JSRV than ENTV-1 in the C-terminal region of the Env TM cytoplasmic tail, so it is possible that some aspect of this region is also important in dissemination. Apart from Orf-x, two small regions in Gag, the transmembrane (TM) region, especially the cytoplasmic tail (CT) of Env, were the main differences between ENTV-2-Shaanxi1 ~ 4 and other ENTV-2. ENTV-2-Shaanxi1 ~ 4 were found to be highly conserved with less than 3.4% amino acid differences and less than 9.3% nucleotide differences in the gag coding region to ENTV-2 and ENTV-2-SC, and less than 0.3% amino acid differences and less than 5.1% nucleotide differences in the env coding region to ENTV-2-SC. The three tyrosine residues (Y590, Y592, and Y596), especially the Y590 in the ENTV-1NA1 (accession number GU292317.1) CT are known to be essential for Env mediated transformation [17, 18, 29], but only two tyrosine residues (Y598 and Y602) were found among ENTV-2-Shaanxi2 ~ 4 isolates, and one tyrosine residues (Y598) were found in ENTV-2-Shaanxi. Y590 residue is essential for tumorigenesis and/or replication of JSRV . Maeda and others highlighted the importance of a single tyrosine residue in the cytoplasmic tail of JSRV TM, Y590. The Y590 is conserved in the ENTV-1 cytoplasmic tail, and mutation of this residue greatly reduces transformation; in contrast, mutation of two other tyrosine residues in the ENTV-1 TM (Y592 and Y596) had relatively less effect [17, 18, 29]. Similarly, the Y598 of ENTV-2 is corresponds to Y590 of ENTV-1 or JSRV. The Y598 is conserved in ENTV-2 cytoplasmic tail. Except ENTV-2-Shaanxi, the other five ENTV-2 sequences have another tyrosine residue Y602. Y590 residue is essential for tumorigenesis and/or replication of JSRV. It likely that only Y598 is important to signal transduction for transformation of ENTV-2. The cytoplasmic tail of the envelope transmembrane (TM) protein is necessary for transformation, and in particular a consensus binding motif (YXXM) for phosphatidylinositol 3-kinase (PI3K) is important. YXXM motif is a reliable molecular marker for the infectious exogenous virus . All the ENTV-2-Shaanxi isolates have the YXXM motif. Maybe the YXXM motif is an essential part for transformation, rather than single amino acids.
The importance of dissemination in the virus life-cycle/pathogenesis remains to be proven, particularly as ENTV-1 and ENTV-2 induce very similar diseases even though dissemination in the host appears to be more limited for ENTV-1 . Future efforts will involve application of our RT-PCR protocol for preclinical diagnosis of ENTV-2 from nasal swabs and construction of an infectious molecular clone of ENTV-2 for subsequent pathogenesis studies. The next studies of molecular epidemiological characterization of ENTV-2 isolates will increase our understanding and will provide knowledge on how to achieve a better understanding of pathways that have led to later transmission of ENTV-2 and how control management should be implemented to prevent the spread of disease.
In this study, the complete sequences were determined from four isolates of Shaanxi province. The ENTV-2-Shaanxi genomes shared 97.0% sequence identity with ENTV-2-SC (accession number HM104174.1), and 89.6% sequence identity with the ENTV-2 sequence (accession number AY197548.1). ENTV-2 is closely related to the ENTV-1 and JSRV. The main sequence differences between these viruses reside in LTR, VR1 and VR2 of Gag, Orf-x, and the transmembrane (TM) region of Env. Phylogenetic analysis by nucleotide sequences showed that four ENTV-2 isolates of shaanxi province were closest related to three ENTV-2 isolates published in NCBI, especially with ENTV-2-SC strain. This finding indicates that ENA most likely was introduced to Shaanxi province by the movement of contaminated goats from other areas in China. This study adds to understand the circulation, variation and distribution of ENTV-2, and may prove beneficial in future control or eradication programmes.
We would like to thank Dr. Kangkang Guo (College of Veterinary Medicine, Northwest A&F University) for technical assistance. We thank Dr. Sanke Yu (College of Veterinary Medicine, Northwest A&F University) for the advice and help in English editing of the manuscript.
This work was supported by grants from the Major Industrial Innovation Chain Projects of Shaanxi Province (No. 2016KTZDNY02-06) and the innovation project for agro-technology of Shaanxi Province (No. 2016NY-092), China.
Availability of data and materials
All relevant information is provided in this current manuscript.
YH and QZ performed the majority of experiments and involved in manuscript preparation, JW participated in editing of the manuscript. MZ participated part of the experiments. MF and XX conceived of the study, participate in its design and coordination, and revised the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
The study was approved by the Ethics Committee of Northwest A&F University. All the animals were treated in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the Ministry of Science and Technology of the People’s Republic of China.
Consent for publication
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ortín A, Cousens C, Minguijón E, Pascual Z, Villarreal MP, Sharp JM, De IHM. Characterization of enzootic nasal tumour virus of goats: complete sequence and tissue distribution. J Gen Virol. 2003;84:2245–52.View ArticlePubMedGoogle Scholar
- Palmarini M, Sharp J, de las Heras M, Fan H. Jaagsiekte sheep retrovirus is necessary and sufficient to induce a contagious lung cancer in sheep. J Virol. 1999;73:6964–72.PubMedPubMed CentralGoogle Scholar
- Walsh SR, Linnerth-Petrik NM, Yu DL, Foster RA, Menzies PI, Diaz-Méndez A, Chalmers HJ, Wootton SK. Experimental transmission of enzootic nasal adenocarcinoma in sheep. Vet Res. 2013;44:66–78.View ArticlePubMedPubMed CentralGoogle Scholar
- De las Heras M, Ortín A, Cousens C, Minguijón E, Sharp J. Enzootic nasal adenocarcinoma of sheep and goats. Curr Top Microbiol Immunol. 2003;275:201–23.PubMedGoogle Scholar
- Walsh SR, Stinson KJ, Wootton SK. Seroconversion of sheep experimentally infected with enzootic nasal tumor virus. BMC Res Notes. 2016;9:1–4.View ArticleGoogle Scholar
- Walsh SR. Pathogenesis of enzootic nasal tumor virus. Guelph: University of Guelph; 2014. https://atrium.lib.uoguelph.ca/xmlui/handle/10214/7808.
- Walsh SR, Linnerth-Petrik NM, Laporte AN, Menzies PI, Foster RA, Wootton SK. Full-length genome sequence analysis of enzootic nasal tumor virus reveals an unusually high degree of genetic stability. Virus Res. 2010;151:74–87.View ArticlePubMedGoogle Scholar
- Feng YC, Yan QG, Guo WZ, Wang XY, Shu L. Construction and bioinformatics analysis of cDNA library of goat enzootic nasal tumor virus SC strain. Chin Vet Sci. 2011;41:126–30.Google Scholar
- Herring AJ, Sharp JM, Scott F, Angus KW. Further evidence for a retrovirus as the aetiological agent of sheep pulmonary adenomatosis (jaagsiekte). Vet Microbiol. 1983;8(3):237–49.View ArticlePubMedGoogle Scholar
- York DF, Vigne R, Verwoerd DW, Querat G. Isolation, identification, and partial cdna cloning of genomic rna of jaagsiekte retrovirus, the etiological agent of sheep pulmonary adenomatosis. J Virol. 1991;65(9):5061–7.PubMedPubMed CentralGoogle Scholar
- Cousens C, Minguijon E, Dalziel R, Ortin A, Garcia M, Park J, Gonzalez L, Sharp J, de las Heras M. Complete sequence of enzootic nasal tumor virus, a retrovirus associated with transmissible intranasal tumors of sheep. J Virol. 1999;73:3986–93.PubMedPubMed CentralGoogle Scholar
- Hallwirth C, Maeda N, York D, Fan H. Variable regions 1 and 2 (VR1 and VR2) in JSRV gag are not responsible for the endogenous JSRV particle release defect. Virus Genes. 2005;30:59–68.View ArticlePubMedGoogle Scholar
- Palmarini M, Hallwirth C, York D, Murgia C, de Oliveira T, Spencer T, Fan H. Molecular cloning and functional analysis of three type D endogenous retroviruses of sheep reveal a different cell tropism from that of the highly related exogenous jaagsiekte sheep retrovirus. J Virol. 2000;74:8065–76.View ArticlePubMedPubMed CentralGoogle Scholar
- Walsh SR, Gerpe MCR, Wootton SK. Construction of a molecular clone of ovine enzootic nasal tumor virus. Virol J. 2016;13:209–118.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosati S, Pittau M, Alberti A, Pozzi S, York D, Sharp J, Palmarini M. An accessory open reading frame (Orf-x) of jaagsiekte sheep retrovirus is conserved between different virus isolates. Virus Res. 2000;66:109–16.View ArticlePubMedGoogle Scholar
- Yang SF, Si JQ, Zhen QW, Liang T, Xiao YY, Chen XK, Yang X, Sheng JL. Epidemiological Survey of Ovine Pulmonary Adenomatosis in Some Regions of Xinjiang. J Shihezi Univ Nat Sci. 2015;33:568–73.Google Scholar
- Maeda N, Inoshima Y, Fruman DA, Brachmann SM, Fan H. Transformation of mouse fibroblasts by jaagsiekte sheep retrovirus envelope does not require phosphatidylinositol 3-kinase. J Virol. 2003;77(18):9951–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu S, Lerman M, Miller A. Putative phosphatidylinositol 3-kinase (PI3K) binding motifs in ovine betaretrovirus Env proteins are not essential for rodent fibroblast transformation and PI3K/Akt activation. J Virol. 2003;77:7924–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Walsh SR, Stinson KJ, Menzies PI, Wootton SK. Development of an ante-mortem diagnostic test for enzootic nasal tumor virus and detection of neutralizing antibodies in host serum. J Gen Virol. 2014;95:1843–54.View ArticlePubMedGoogle Scholar
- Liu FX, Lin X, FYH H, Hao XP, Zhao ZH, Yu WS, Wang HZ. An Investigation and pathologicmorphological Study on Tumors of Domestic Animals in Inner Mongolia. Vet Orientation. 1984;3:15–25.Google Scholar
- Lei HY, Su JM, Ning LZ, Kang Y, Chen KY, Zeng DN. Investigation and diagnosis of an enzootic nasal tumor of goat. Prog Vet Med. 2006;27:112–4.Google Scholar
- Liu F, Feng YC, Yan QG, Han GQ. Diagnosis of an enzootic nasal tumor of goat in Sichuan. Anim Husb Vet Med. 2011;43:83–6.Google Scholar
- Yu YD, Wei LF, Huang XJ, Zhang B, Yu XH, Tang C. Diagnosis of Four Cases of Enzootic Nasal Adenocarcinoma of Goats. Prog Vet Med. 2014;35:129–31.Google Scholar
- Fox KA, Wootton SK, Quackenbush SL, Wolfe LL, Levan IK, Miller MW, et al. Paranasal sinus masses of rocky mountain bighorn sheep (Ovis canadensis canadensis). Vet Pathol. 2011;48:706–12.View ArticlePubMedGoogle Scholar
- Fox KA, Rouse NM, Huyvaert KP, Griffin KA, Killion HJ, Jennings-Gaines J, et al. Bighorn sheep (Ovis canadensis) sinus tumors are associated with coinfections by potentially pathogenic bacteria in the upper respiratory tract. J Wildl Dis. 2014;51(1):19–27.View ArticleGoogle Scholar
- Fox KA, Wootton S, Marolf A, Rouse N, Levan I, Spraker T, et al. Experimental transmission of bighorn sheep sinus tumors to bighorn sheep (Ovis canadensis canadensis) and domestic sheep. Vet Pathol. 2016;53(6):1164–71.View ArticlePubMedGoogle Scholar
- York DF, Vigne R, Verwoerd DW, Querat G. Nucleotide sequence of the jaagsiekte retrovirus, an exogenous and endogenous type D and B retrovirus of sheep and goats. J Virol. 1992;66:4930–9.PubMedPubMed CentralGoogle Scholar
- Cousens C, Maeda N, Murgia C, Dagleish MP, Palmarini M, Fan H. In vivo tumorigenesis by jaagsiekte sheep retrovirus (jsrv) requires y590 in env tm, but not full-length orfx open reading frame. Virologie. 2007;367:413–21.View ArticleGoogle Scholar
- Maeda N, Fan H. Signal transduction pathways utilized by enzootic nasal tumor virus (entv-1) envelope protein in transformation of rat epithelial cells resemble those used by jaagsiekte sheep retrovirus. Virus Genes. 2008;36:147–55.View ArticlePubMedGoogle Scholar