Complete DNA sequence of the rat cytomegalovirus genome

We have determined the complete genome sequence of the Maastricht strain of rat cytomegalovirus (RCMV). The RCMV genome has a length of 229,896 bp and is arranged as a single unique sequence flanked by 504-bp terminal direct repeats. RCMV was found to have counterparts of all but one of the open reading frames (ORFs) that are conserved between murine CMV (MCMV) and human CMV (HCMV). Like HCMV, RCMV lacks homologs of the genes belonging to the MCMV m02 glycoprotein gene family. However, RCMV contains 15 ORFs with homology to members of the MCMV m145 glycoprotein gene family. Four ORFs are predicted to encode homologs of host proteins; R33 and R78 both putatively encode G protein-coupled receptors, whereas r144 and r131 encode homologs of major histocompatibility class I heavy chains and CC chemokines, respectively. An intriguing feature of the RCMV genome is the presence of an ORF, r127, with similarity to the rep gene of parvoviruses as well as ORF U94 of human herpesvirus 6A (HHV-6A) and HHV-6B. Counterparts of these ORFs have not been found in the other sequenced herpesviruses.     

Analysis and characterization of the complete genome of tupaia (tree shrew) herpesvirus

The tupaia herpesvirus (THV) was isolated from spontaneously degenerating tissue cultures of malignant lymphoma, lung, and spleen cell cultures of tree shrews (Tupaia spp.). The determination of the complete nucleotide sequence of the THV strain 2 genome resulted in a 195,857-bp-long, linear DNA molecule with a G+C content of 66.5%. The terminal regions of the THV genome and the loci of conserved viral genes were found to be G+C richer. Furthermore, no large repetitive DNA sequences could be identified. This is in agreement with the previous classification of THV as the prototype species of herpesvirus genome group F. The search for potential coding regions resulted in the identification of 158 open reading frames (ORFs) regularly distributed on both DNA strands. Seventy-six out of the 158 ORFs code for proteins that are significantly homologous to known herpesvirus proteins. The highest homologies found were to primate and rodent cytomegaloviruses. Biological properties, protein homologies, the arrangement of conserved viral genes, and phylogenetic analysis revealed that THV is a member of the subfamily Betaherpesvirinae. The evolutionary lineages of THV and the cytomegaloviruses seem to have branched off from a common ancestor. In addition, it was found that the arrangements of conserved genes of THV and murine cytomegalovirus strain Smith, both of which are not able to form genomic isomers, are colinear with two different human cytomegalovirus (HCMV) strain AD169 genomic isomers that differ from each other in the orientation of the long unique region. The biological properties and the high degree of relatedness of THV to the mammalian cytomegaloviruses allow the consideration of THV as a model system for investigation of HCMV pathogenicity.     

In silico pattern-based analysis of the human cytomegalovirus genome

More than 200 open reading frames (ORFs) from the human cytomegalovirus genome have been reported as potentially coding for proteins. We have used two pattern-based in silico approaches to analyze this set of putative viral genes. With the help of an objective annotation method that is based on the Bio-Dictionary, a comprehensive collection of amino acid patterns that describes the currently known natural sequence space of proteins, we have reannotated all of the previously reported putative genes of the human cytomegalovirus. Also, with the help of MUSCA, a pattern-based multiple sequence alignment algorithm, we have reexamined the original human cytomegalovirus gene family definitions. Our analysis of the genome shows that many of the coded proteins comprise amino acid combinations that are unique to either the human cytomegalovirus or the larger group of herpesviruses. We have confirmed that a surprisingly large portion of the analyzed ORFs encode membrane proteins, and we have discovered a significant number of previously uncharacterized proteins that are predicted to be G-protein-coupled receptor homologues. The analysis also indicates that many of the encoded proteins undergo posttranslational modifications such as hydroxylation, phosphorylation, and glycosylation. ORFs encoding proteins with similar functional behavior appear in neighboring regions of the human cytomegalovirus genome. All of the results of the present study can be found and interactively explored online (http://cbcsrv.watson.ibm.com/virus/).     

Complete sequence and genomic analysis of murine gammaherpesvirus 68.

Murine gammaherpesvirus 68 (gammaHV68) infects mice, thus providing a tractable small-animal model for analysis of the acute and chronic pathogenesis of gammaherpesviruses. To facilitate molecular analysis of gammaHV68 pathogenesis, we have sequenced the gammaHV68 genome. The genome contains 118,237 bp of unique sequence flanked by multiple copies of a 1,213-bp terminal repeat. The GC content of the unique portion of the genome is 46%, while the GC content of the terminal repeat is 78%. The unique portion of the genome is estimated to encode at least 80 genes and is largely colinear with the genomes of Kaposi's sarcoma herpesvirus (KSHV; also known as human herpesvirus 8), herpesvirus saimiri (HVS), and Epstein-Barr virus (EBV). We detected 63 open reading frames (ORFs) homologous to HVS and KSHV ORFs and used the HVS/KSHV numbering system to designate these ORFs. gammaHV68 shares with HVS and KSHV ORFs homologous to a complement regulatory protein (ORF 4), a D-type cyclin (ORF 72), and a G-protein-coupled receptor with close homology to the interleukin-8 receptor (ORF 74). One ORF (K3) was identified in gammaHV68 as homologous to both ORFs K3 and K5 of KSHV and contains a domain found in a bovine herpesvirus 4 major immediate-early protein. We also detected 16 methionine-initiated ORFs predicted to encode proteins at least 100 amino acids in length that are unique to gammaHV68 (ORFs M1 to 14). ORF M1 has striking homology to poxvirus serpins, while ORF M11 encodes a potential homolog of Bcl-2-like molecules encoded by other gammaherpesviruses (gene 16 of HVS and KSHV and the BHRF1 gene of EBV). In addition, clustered at the left end of the unique region are eight sequences with significant homology to bacterial tRNAs. The unique region of the genome contains two internal repeats: a 40-bp repeat located between bp 26778 and 28191 in the genome and a 100-bp repeat located between bp 98981 and 101170. Analysis of the gammaHV68, HVS, EBV, and KSHV genomes demonstrated that each of these viruses have large colinear gene blocks interspersed by regions containing virus-specific ORFs. Interestingly, genes associated with EBV cell tropism, latency, and transformation are all contained within these regions encoding virus-specific genes. This finding suggests that pathogenesis-associated genes of gammaherpesviruses, including gammaHV68, may be contained in similarly positioned genome regions. The availability of the gammaHV68 genomic sequence will facilitate analysis of critical issues in gammaherpesvirus biology via integration of molecular and pathogenetic studies in a small-animal model.     

Sequence analysis of the 4.7-kb BamHI-EcoRI fragment of the equine herpesvirus type-1 short unique region.

To localize gene that may encode immunogens potentially important for recombinant vaccine design, we have analysed a region of the equine herpesvirus type-1 (EHV-1) genome where a glycoprotein-encoding gene had previously been mapped. The 4707-bp BamHI-EcoRI fragment from the short unique region of the EHV-1 genome was sequenced. This sequence contains three entire open reading frames (ORFs), and portions of two more. ORF1 codes for 161 amino acids (aa), and represents the C terminus of a possible membrane-bound protein. ORF2 (424 aa) and ORF3 (550 aa) are potential glycoprotein-encoding genes; the predicted aa sequences contain possible signal sequences, N-linked glycosylation sites and transmembrane domains; they also show homology to the glycoproteins gI and gE of herpes simplex virus type-1 (HSV-1), and the related proteins of pseudorabies virus and varicella-zoster virus. The predicted aa sequence of ORF4 shares no homology with other known herpesvirus proteins, but the nucleotide sequence shows a high level of homology with the corresponding region of the EHV-4 genome. ORF5 may be related to US9 of HSV-1.     

Psittacid herpesvirus 1 and infectious laryngotracheitis virus: Comparative genome sequence analysis of two avian alphaherpesviruses

Psittacid herpesvirus 1 (PsHV-1) is the causative agent of Pacheco's disease, an acute, highly contagious, and potentially lethal respiratory herpesvirus infection in psittacine birds, while infectious laryngotracheitis virus (ILTV) is a highly contagious and economically significant avian herpesvirus which is responsible for an acute respiratory disease limited to galliform birds. The complete genome sequence of PsHV-1 has been determined and compared to the ILTV sequence, assembled from published data. The PsHV-1 and ILTV genomes exhibit similar structural characteristics and are 163,025 bp and 148,665 bp in length, respectively. The PsHV-1 genome contains 73 predicted open reading frames (ORFs), while the ILTV genome contains 77 predicted ORFs. Both genomes contain an inversion in the unique long region similar to that observed in pseudorabies virus. PsHV-1 is closely related to ILTV, and it is proposed that it be assigned to the Iltovirus genus. These two avian herpesviruses represent a phylogenetically unique clade of alphaherpesviruses that are distinct from the Marek's disease-like viruses (Mardivirus). The determination of the complete genomic nucleotide sequences of PsHV-1 and ILTV provides a tool for further comparative and functional analysis of this unique class of avian alphaherpesviruses.     

Structure of the rat cytomegalovirus genome termini.

The lytic replication cycle of herpesviruses can be divided into the following three steps: (i) circularization, in which, after infection, the termini of the linear double-stranded viral genome are fused; (ii) replication, in which the circular DNA serves as template for DNA replication, which generates large DNA concatemers; and (iii) maturation, in which the concatemeric viral DNA is processed into unit-length genomes, which are packaged into capsids. Sequences at the termini of the linear virion DNA are thought to play a key role in both genome circularization and maturation. To investigate the mechanism of these processes in the replication of rat cytomegalovirus (RCMV), we cloned, sequenced, and characterized the genomic termini of this betaherpesvirus. Both RCMV genomic termini were found to contain a single copy of a direct terminal repeat (TR). The TR sequence is 504 bp in length, has a high GC content (76%), and is not repeated at internal sites within the RCMV genome. The TR comprises several small internal direct repeats as well as two sequences which are homologous to herpesvirus pac-1 and pac-2 sites, respectively. The organization of the RCMV TR is unique among cytomegaloviruses with respect to the position of the pac sequences: pac-1 is located near the left end of the TR, whereas pac-2 is present near the right end. Both RCMV DNA termini carry an extension of a single nucleotide at the 3' end. Since these nucleotides are complementary, circularization of the viral genome is likely to occur via a simple ligation reaction.     

The genome of S-PM2, a "photosynthetic" T4-type bacteriophage that infects marine Synechococcus strains

Bacteriophage S-PM2 infects several strains of the abundant and ecologically important marine cyanobacterium Synechococcus. A large lytic phage with an isometric icosahedral head, S-PM2 has a contractile tail and by this criterion is classified as a myovirus (1). The linear, circularly permuted, 196,280-bp double-stranded DNA genome of S-PM2 contains 37.8% G+C residues. It encodes 239 open reading frames (ORFs) and 25 tRNAs. Of these ORFs, 19 appear to encode proteins associated with the cell envelope, including a putative S-layer-associated protein. Twenty additional S-PM2 ORFs have homologues in the genomes of their cyanobacterial hosts. There is a group I self-splicing intron within the gene encoding the D1 protein. A total of 40 ORFs, organized into discrete clusters, encode homologues of T4 proteins involved in virion morphogenesis, nucleotide metabolism, gene regulation, and DNA replication and repair. The S-PM2 genome encodes a few surprisingly large (e.g., 3,779 amino acids) ORFs of unknown function. Our analysis of the S-PM2 genome suggests that many of the unknown S-PM2 functions may be involved in the adaptation of the metabolism of the host cell to the requirements of phage infection. This hypothesis originates from the identification of multiple phage-mediated modifications of the host's photosynthetic apparatus that appear to be essential for maintaining energy production during the lytic cycle.     

The primary sequence of rhesus monkey rhadinovirus isolate 26-95: sequence similarities to Kaposi's sarcoma-associated herpesvirus and rhesus monkey rhadinovirus isolate 17577

The primary sequence of the long unique region L-DNA (L for low GC) of rhesus monkey rhadinovirus (RRV) isolate 26-95 was determined. The L-DNA consists of 130,733 bp that contain 84 open reading frames (ORFs). The overall organization of the RRV26-95 genome was found to be very similar to that of human Kaposi sarcoma-associated herpesvirus (KSHV). BLAST search analysis revealed that in almost all cases RRV26-95 coding sequences have a greater degree of similarity to corresponding KSHV sequences than to other herpesviruses. All of the ORFs present in KSHV have at least one homologue in RRV26-95 except K3 and K5 (bovine herpesvirus-4 immediate-early protein homologues), K7 (nut-1), and K12 (Kaposin). RRV26-95 contains one MIP-1 and eight interferon regulatory factor (vIRF) homologues compared to three MIP-1 and four vIRF homologues in KSHV. All homologues are correspondingly located in KSHV and RRV with the exception of dihydrofolate reductase (DHFR). DHFR is correspondingly located near the left end of the genome in RRV26-95 and herpesvirus saimiri (HVS), but in KSHV the DHFR gene is displaced 16,069 nucleotides in a rightward direction in the genome. DHFR is also unusual in that the RRV26-95 DHFR more closely resembles HVS DHFR (74% similarity) than KSHV DHFR (55% similarity). Of the 84 ORFs in RRV26-95, 83 contain sequences similar to the recently determined sequences of the independent RRV isolate 17577. RRV26-95 and RRV17577 sequences differ in that ORF 67.5 sequences contained in RRV26-95 were not found in RRV17577. In addition, ORF 4 is significantly shorter in RRV26-95 than was reported for RRV17577 (395 versus 645 amino acids). Only four of the corresponding ORFs between RRV26-95 and RRV17577 exhibited less than 95% sequence identity: glycoproteins H and L, uracil DNA glucosidase, and a tegument protein (ORF 67). Both RRV26-95 and RRV17577 have unique ORFs between positions 21444 to 21752 and 110910 to 114899 in a rightward direction and from positions 116524 to 111082 in a leftward direction that are not found in KSHV. Our analysis indicates that RRV26-95 and RRV17577 are clearly independent isolates of the same virus species and that both are closely related in structural organization and overall sequence to KSHV. The availability of detailed sequence information, the ability to grow RRV lytically in cell culture, and the ability to infect monkeys experimentally with RRV will facilitate the construction of mutant strains of virus for evaluating the contribution of individual genes to biological properties.     

Identification and sequencing of a novel rodent gammaherpesvirus that establishes acute and latent infection in laboratory mice

Gammaherpesviruses encode numerous immunomodulatory molecules that contribute to their ability to evade the host immune response and establish persistent, lifelong infections. As the human gammaherpesviruses are strictly species specific, small animal models of gammaherpesvirus infection, such as murine gammaherpesvirus 68 (γHV68) infection, are important for studying the roles of gammaherpesvirus immune evasion genes in in vivo infection and pathogenesis. We report here the genome sequence and characterization of a novel rodent gammaherpesvirus, designated rodent herpesvirus Peru (RHVP), that shares conserved genes and genome organization with γHV68 and the primate gammaherpesviruses but is phylogenetically distinct from γHV68. RHVP establishes acute and latent infection in laboratory mice. Additionally, RHVP contains multiple open reading frames (ORFs) not present in γHV68 that have sequence similarity to primate gammaherpesvirus immunomodulatory genes or cellular genes. These include ORFs with similarity to major histocompatibility complex class I (MHC-I), C-type lectins, and the mouse mammary tumor virus and herpesvirus saimiri superantigens. As these ORFs may function as immunomodulatory or virulence factors, RHVP presents new opportunities for the study of mechanisms of immune evasion by gammaherpesviruses.     

Sequences, annotation and single nucleotide polymorphism of the major histocompatibility complex in the domestic cat

Two sequences of major histocompatibility complex (MHC) regions in the domestic cat, 2.976 and 0.362 Mbps, which were separated by an ancient chromosome break (55-80 MYA) and followed by a chromosomal inversion were annotated in detail. Gene annotation of this MHC was completed and identified 183 possible coding regions, 147 human homologues, possible functional genes and 36 pseudo/unidentified genes) by GENSCAN and BLASTN, BLASTP RepeatMasker programs. The first region spans 2.976 Mbp sequence, which encodes six classical class II antigens (three DRA and three DRB antigens) lacking the functional DP, DQ regions, nine antigen processing molecules (DOA/DOB, DMA/DMB, TAPASIN, and LMP2/LMP7,TAP1/TAP2), 52 class III genes, nineteen class I genes/gene fragments (FLAI-A to FLAI-S). Three class I genes (FLAI-H, I-K, I-E) may encode functional classical class I antigens based on deduced amino acid sequence and promoter structure. The second region spans 0.362 Mbp sequence encoding no class I genes and 18 cross-species conserved genes, excluding class I, II and their functionally related/associated genes, namely framework genes, including three olfactory receptor genes. One previously identified feline endogenous retrovirus, a baboon retrovirus derived sequence (ECE1) and two new endogenous retrovirus sequences, similar to brown bat endogenous retrovirus (FERVmlu1, FERVmlu2) were found within a 140 Kbp interval in the middle of class I region. MHC SNPs were examined based on comparisons of this BAC sequence and MHC homozygous 1.9x WGS sequences and found that 11,654 SNPs in 2.84 Mbp (0.00411 SNP per bp), which is 2.4 times higher rate than average heterozygous region in the WGS (0.0017 SNP per bp genome), and slightly higher than the SNP rate observed in human MHC (0.00337 SNP per bp).     

Complete genomic sequence of bacteriophage B3, a Mu-like phage of Pseudomonas aeruginosa

Bacteriophage B3 is a transposable phage of Pseudomonas aeruginosa. In this report, we present the complete DNA sequence and annotation of the B3 genome. DNA sequence analysis revealed that the B3 genome is 38,439 bp long with a G+C content of 63.3%. The genome contains 59 proposed open reading frames (ORFs) organized into at least three operons. Of these ORFs, the predicted proteins from 41 ORFs (68%) display significant similarity to other phage or bacterial proteins. Many of the predicted B3 proteins are homologous to those encoded by the early genes and head genes of Mu and Mu-like prophages found in sequenced bacterial genomes. Only two of the predicted B3 tail proteins are homologous to other well-characterized phage tail proteins; however, several Mu-like prophages and transposable phage D3112 encode approximately 10 highly similar proteins in their predicted tail gene regions. Comparison of the B3 genomic organization with that of Mu revealed evidence of multiple genetic rearrangements, the most notable being the inversion of the proposed B3 immunity/early gene region, the loss of Mu-like tail genes, and an extreme leftward shift of the B3 DNA modification gene cluster. These differences illustrate and support the widely held view that tailed phages are genetic mosaics arising by the exchange of functional modules within a diverse genetic pool.     

Complete genome sequence of the shrimp white spot bacilliform virus.

We report the first complete genome sequence of a marine invertebrate virus. White spot bacilliform virus (WSBV; or white spot syndrome virus) is a major shrimp pathogen with a high mortality rate and a wide host range. Its double-stranded circular DNA genome of 305,107 bp contains 181 open reading frames (ORFs). Nine homologous regions containing 47 repeated minifragments that include direct repeats, atypical inverted repeat sequences, and imperfect palindromes were identified. This is the largest animal virus that has been completely sequenced. Although WSBV is morphologically similar to insect baculovirus, the two viruses are not detectably related at the amino acid level. Rather, some WSBV genes are more homologous to eukaryotic genes than viral genes. In fact, sequence analysis indicates that WSBV differs from all known viruses, although a few genes display a weak homology to herpesvirus genes. Most of the ORFs encode proteins that bear no homology to any known proteins, either suggesting that WSBV represents a novel class of viruses or perhaps implying a significant evolutionary distance between marine and terrestrial viruses. The most unique feature of WSBV is the presence of an intact collagen gene, a gene encoding an extracellular matrix protein of animal cells that has never been found in any viruses. Determination of the genome of WSBV will facilitate a better understanding of the molecular mechanism underlying the pathogenesis of the WSBV virus and will also provide useful information concerning the evolution and divergence of marine and terrestrial animal viruses at the molecular level.     

Cytomegalovirus US2 destroys two components of the MHC class II pathway, preventing recognition by CD4+ T cells.

Human cytomegalovirus (HCMV) is a ubiquitous herpesvirus that causes life-threatening disease in patients who are immunosuppressed for bone marrow or tissue transplantation or who have AIDS (ref. 1). HCMV establishes lifelong latent infections and, after periodic reactivation from latency, uses a panel of immune evasion proteins to survive and replicate in the face of robust, fully primed host immunity. Monocyte/macrophages are important host cells for HCMV, serving as a latent reservoir and as a means of dissemination throughout the body. Macrophages and other HCMV-permissive cells, such as endothelial and glial cells, can express MHC class II proteins and present antigens to CD4+ T lymphocytes. Here, we show that the HCMV protein US2 causes degradation of two essential proteins in the MHC class II antigen presentation pathway: HLA-DR-alpha and DM-alpha. This was unexpected, as US2 has been shown to cause degradation of MHC class I (refs. 5,6), which has only limited homology with class II proteins. Expression of US2 in cells reduced or abolished their ability to present antigen to CD4+ T lymphocytes. Thus, US2 may allow HCMV-infected macrophages to remain relatively 'invisible' to CD4+ T cells, a property that would be important after virus reactivation.     

Sequence and genomic analysis of a Rhesus macaque rhadinovirus with similarity to Kaposi's sarcoma-associated herpesvirus/human herpesvirus 8

We have sequenced the long unique region (LUR) and characterized the terminal repeats of the genome of a rhesus rhadinovirus (RRV), strain 17577. The LUR as sequenced is 131,364 bp in length, with a G+C content of 52.2% and a CpG ratio of 1.11. The genome codes for 79 open reading frames (ORFs), with 67 of these ORFs similar to genes found in both Kaposi's sarcoma-associated herpesvirus (KSHV) (formal name, human herpesvirus 8) and herpesvirus saimiri. Eight of the 12 unique genes show similarity to genes found in KSHV, including genes for viral interleukin-6, viral macrophage inflammatory protein, and a family of viral interferon regulatory factors (vIRFs). Genomic organization is essentially colinear with KSHV, the primary differences being the number of cytokine and IRF genes and the location of the gene for dihydrofolate reductase. Highly repetitive sequences are located in positions corresponding to repetitive sequences found in KSHV. Phylogenetic analysis of several ORFs supports the similarity between RRV and KSHV. Overall, the sequence, structural, and phylogenetic data combine to provide strong evidence that RRV 17577 is the rhesus macaque homolog of KSHV.     

Sequence of the short unique region, short repeats, and part of the long repeats of human cytomegalovirus

We have determined the DNA sequence (46 kilobases) of the short unique region, the short repeat, and part of the long repeat of human cytomegalovirus strain AD169. Analysis of the sequence has revealed at least 38 possible regions that may code for protein. Many of these open reading frames show homology to each other, and five groups of homologous reading frames are identified. Half of the predicted translation products appear to be membrane proteins, and fall into two distinct classes; those that have potential signal and anchor sequences, and those that have seven potential membrane-spanning regions and appear to be integral membrane proteins. A number of the former class contain sites for N-linked glycosylation and may therefore be glycoproteins. None of the 38 open reading frames shows homology to other known herpesvirus proteins.