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.     

Conservation of gene organization in the lymphotropic herpesviruses herpesvirus Saimiri and Epstein-Barr virus.

By analyses of short DNA sequences, we have deduced the overall arrangement of genes in the (A + T)-rich coding sequences of herpesvirus saimiri (HVS) relative to the arrangements of homologous genes in the (G + C)-rich coding sequences of the Epstein-Barr virus (EBV) genome and the (A + T)-rich sequences of the varicella-zoster virus (VZV) genome. Fragments of HVS DNA from 13 separate sites within the 111 kilobase pairs of the light DNA coding sequences of the genome were subcloned into M13 vectors, and sequences of up to 350 bases were determined from each of these sites. Amino acid sequences predicted for fragments of open reading frames defined by these sequences were compared with a library of the protein sequences of major open reading frames predicted from the complete DNA sequences of VZV and EBV. Of the 13 short amino acid sequences obtained from HVS, only 3 were recognizably homologous to proteins encoded by VZV, but all 13 HVS sequences were unambiguously homologous to gene products encoded by EBV. The HVS reading frames identified by this method included homologs of the major capsid polypeptides, glycoprotein H, the major nonstructural DNA-binding protein, thymidine kinase, and the homolog of the regulatory gene product of the BMLF1 reading frame of EBV. Locally as well as globally, the order and relative orientation of these genes resembled that of their homologs on the EBV genome. Despite the major differences in their nucleotide compositions and in the nature and arrangements of reiterated DNA sequences, the genomes of the lymphotropic herpesviruses HVS and EBV encode closely related proteins, and they share a common organization of these coding sequences which differs from that of the neurotropic herpesviruses, VZV and herpes simplex virus.     

Identification of the Rhesus Macaque Rhadinovirus Lytic Origin of DNA Replication

We have identified a lytic origin of DNA replication (oriLyt) for rhesus macaque rhadinovirus (RRV), the rhesus macaque homolog of human herpesvirus 8 (HHV-8), also known as Kaposi's sarcoma-associated herpesvirus. RRV oriLyt maps to the region of the genome between open reading frame 69 (ORF69) and ORF71 (vFLIP) and is composed of an upstream A+T-rich region followed by a short (300-bp) downstream G+C-rich DNA sequence. A set of overlapping cosmids corresponding to the entire genome of RRV was capable of complementing oriLyt-dependent DNA replication only when additional ORF50 was supplied as an expression plasmid in the transfection mixture, suggesting that the level of ORF50 protein originating from input cosmid DNA was insufficient. The requirement of RRV ORF50 in the cotransfection replication assay may also suggest a direct role for this protein in DNA replication. RRV oriLyt shares a high degree of nucleotide sequence and G+C base distribution with the corresponding loci in HHV-8.     

Kaposi's sarcoma-associated herpesvirus contains G protein-coupled receptor and cyclin D homologs which are expressed in Kaposi's sarcoma and malignant lymphoma.

A new human herpesvirus was recently identified in all forms of Kaposi's sarcoma (Kaposi's sarcoma-associated herpesvirus [KSHV] or human herpesvirus 8), as well as in primary effusion (body cavity-based) lymphomas (PELs). A 12.3-kb-long KSHV clone was obtained from a PEL genomic library. Sequencing of this clone revealed extensive homology and colinearity with the right end of the herpesvirus saimiri (HVS) genome and more limited homology to the left end of the Epstein-Barr virus genome. Four open reading frames (ORFs) were sequenced and characterized; these are homologous to the following viral and/or cellular genes: (i) Epstein-Barr virus membrane antigen p140 and HVS p160, (ii) HVS and cellular type D cyclins, (iii) HVS and cellular G protein-coupled receptors, and (iv) HVS. Since there is considerable evidence that cyclin D1 and some G protein-coupled receptors contribute to the development of specific cancers, the presence of KSHV homologs of these genes provides support for a role for KSHV in malignant transformation. All ORFs identified are transcribed in PELs and Kaposi's sarcoma tissues, further suggesting an active role for KSHV in these diseases.     

Primary structure of the herpesvirus saimiri genome.

This report describes the complete nucleotide sequence of the genome of herpesvirus saimiri, the prototype of gammaherpesvirus subgroup 2 (rhadinoviruses). The unique low-G + C-content DNA region has 112,930 bp with an average base composition of 34.5% G + C and is flanked by about 35 noncoding high-G + C-content DNA repeats of 1,444 bp (70.8% G + C) in tandem orientation. We identified 76 major open reading frames and a set of seven U-RNA genes for a total of 83 potential genes. The genes are closely arranged, with only a few regions of sizable noncoding sequences. For 60 of the predicted proteins, homologous sequences are found in other herpesviruses. Genes conserved between herpesvirus saimiri and Epstein-Barr virus (gammaherpesvirus subgroup 1) show that their genomes are generally collinear, although conserved gene blocks are separated by unique genes that appear to determine the particular phenotype of these viruses. Several deduced protein sequences of herpesvirus saimiri without counterparts in most of the other sequenced herpesviruses exhibited significant homology with cellular proteins of known function. These include thymidylate synthase, dihydrofolate reductase, complement control proteins, the cell surface antigen CD59, cyclins, and G protein-coupled receptors. Searching for functional protein motifs revealed that the virus may encode a cytosine-specific methylase and a tyrosine-specific protein kinase. Several herpesvirus saimiri genes are potential candidates to cooperate with the gene for saimiri transformation-associated protein of subgroup A (STP-A) in T-lymphocyte growth stimulation.     

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.     

Diversity in G + C content at the third position of codons in vertebrate genes and its cause.

Correlation was positive between the G + C content at the codon third position in genes of vertebrates and the G + C content of the genome portion surrounding each gene. Exons of genes with a high G + C% at the codon 3rd position are surrounded by G + C-rich introns and G + C-rich flanking sequences, and those with a low G + C% at the position by A + T-rich introns and flanking sequences. Analysis of G + C content distribution along DNA sequences using a DNA Sequence Data Bank supported the view that the vertebrate genome is a mosaic of regions with clear differences in their G + C content. The biological significance of the variation in G + C content throughout the vertebrate genome is discussed in connection with chromosomal banding.     

Genome organization of herpesvirus aotus type 2.

Herpesvirus aotus type 2, a virus commonly found in owl monkeys without overt disease, has a similar genome structure to the oncogenic herpesviruses of nonhuman primates (herpesvirus saimiri, herpesvirus ateles). Virion DNA of herpesvirus aotus type 2 (M-DNA) has an unique 110-kilobase-pair region of low G + C content (40.2%, L-DNA), inserted between stretches of repetitive H-DNA (68.7% G + C, about 41 kilobase pairs per molecule) that are variable in length. A minority of virions contain defective genomes that consist of repetitive H-DNA only. The H-DNA is composed of various types of repeat units that are related in sequence with each other. The two dominant types of repeats (2.3 and 2.7 kilobase pairs) were cloned and compared by restriction enzyme cleavages and partial nucleotide sequencing. They are homologous in at least 1.3 kilobase pairs. The two forms of repeat units are randomly arranged and oriented in tandem. Reassociation kinetics did not allow detection of sequence homologies between H- and L-DNA of herpesvirus aotus type 2 and the respective sequences of oncogenic primate herpesviruses.     

The complete nucleotide sequence of the maize chlorotic mottle virus genome.

The complete nucleotide sequence of the maize chlorotic mottle virus (MCMV) genome has been determined to be 4437 nucleotides. The viral genome has four long open reading frames (ORFs) which could encode polypeptides of 31.6, 50, 8.9 and 25.1 kd. If the termination codons, for the polypeptides encoded by the 50 and 8.9 kd ORFs are suppressed, readthrough products of 111 and 32.7 kd result. The 31.6 and 50 kd ORFs overlap for nearly the entire length of the 31.6 kd ORF. Striking amino acid homology has been observed between two potential polypeptides encoded by MCMV and polypeptides encoded by carnation mottle virus (CarMV) and turnip crinkle virus (TCV). The 25.1 kd ORF most likely encodes the capsid protein. The similar genome organization and amino acid sequence homology of MCMV with CarMV and TCV suggest an evolutionary relationship with these members of the carmovirus group.     

Human herpesvirus 6 is closely related to human cytomegalovirus.

A sequence of 21,858 base pairs from the genome of human herpesvirus 6 (HHV-6) strain U1102 is presented. The sequence has a mean composition of 41% G + C, and the observed frequency of CpG dinucleotides is close to that predicted from this mononucleotide composition. The sequence contains 17 complete open reading frames (ORFs) and part of another at the 5' end of the sequence. The predicted protein products of two of these ORFs have no recognizable homologs in the genomes of other sequenced human herpesviruses (i.e., Epstein-Barr virus [EBV], human cytomegalovirus [HCMV], herpes simplex virus [HSV], and varicella-zoster virus [VZV]). However, the products of nine other ORFs are clearly homologous to a set of genes that is conserved in all other sequenced herpesviruses, including homologs of the alkaline exonuclease, the phosphotransferase, the spliced ORF, and the major capsid protein genes. Measurements of similarity between these homologous sequences showed that HHV-6 is clearly most closely related to HCMV. The degree of relatedness between HHV-6 and HCMV was commensurate with that observed in comparisons between HSV and VZV or EBV and herpesvirus saimiri and significantly greater than its relatedness to EBV, HSV, or VZV. In addition, the gene for the major capsid protein and its 5' neighbor are reoriented with respect to the spliced ORFs in the genomes of both HHV-6 and HCMV relative to the organization observed in EBV, HSV, and VZV. Three ORFs in HHV-6 have recognizable homologs only in the genome of HCMV. Despite differences in gross composition and size, we conclude that the genomes of HHV-6 and HCMV are closely related.     

The Complete Nucleotide Sequence of Odontoglossum Ringspot Virus (Cy-1 Strain) Genomic RNA

The complete nucleotide sequence of the genomic RNA of odontoglossum ringspot virus Cy-1 strain (ORSV Cy-1) was determined using cloned cDNA. This sequence is 6611 nucleotides long containing four open reading frames, which correspond to 126 K, 183 K, 31 K, and 18 K proteins. Its genomic organization is similar to other tobamoviruses, TMV-V(vulgare), TMV-L (tomato strain), tobacco mild green mosaic virus (TMGMV) and cucumber green mottle mosaic virus (CGMMV). The 5′ non-coding regions of ORSV Cy-1 is 62 nucleotides. The ORFs encoded a 126 K polypeptide and a 183 K read-through product in which helicase-sequence and polymerase-sequence motifs are found. The ORFs encoding the 126 K and 183 K proteins have 61% and 63% identities with those of TMV-V. The third ORF encoded a 31 K protein homologous to TMV cell-to-cell movement protein. It has 63% identities with that of TMV-V. The fourth ORF encoded an 18 K coat protein. The 5′ non-coding region, which extends from base 1 to 62 has 2 G residues and a ribosome binding site (AUU). The 3′ non-coding region, 414 nucleotides in length, is entirely different from that of other tobamoviruses.     

ORF73 of herpesvirus Saimiri strain C488 tethers the viral genome to metaphase chromosomes and binds to cis-acting DNA sequences in the terminal repeats

Kaposi's sarcoma (KS)-associated herpesvirus (KSHV), a human oncogenic gamma-2-herpesvirus, transforms human endothelial cells and establishes latent infection at a low efficiency in vitro. During latent infection, only a limited number of genes are expressed, and the circularized viral genome is maintained as a multicopy episome. Latency-associated nuclear antigen (LANA), exclusively expressed during latency, has been shown to have a multifunctional role in KS pathogenesis. LANA tethers the viral episome to the host chromosome, thus ensuring efficient persistence of the viral genome during successive rounds of cell division. Besides episome maintenance, LANA modulates the expression of genes of various cellular and viral pathways, including those of retinoblastoma protein and p53. Herpesvirus saimiri (HVS), another gamma-2-herpesvirus, primarily infects New World primates. Orf73, encoding the nuclear antigen of HVS, is the positional homolog of the LANA gene, and the ORF73 protein has some sequence homology to KSHV LANA. However, the function of ORF73 of HVS has not been thoroughly investigated. In this report, we show that HVS ORF73 may be important for episome persistence and colocalizes with the HVS genomic DNA on metaphase chromosomes. Furthermore, HVS terminal repeats (TRs) contain a cis-acting sequence similar to that in KSHV TRs, suggesting that the LANA binding sequence is conserved between these two viruses. This cis-acting element is sufficient to bind HVS ORF73 from strains C488 and A11, and plasmids containing the HVS C488 TR element are maintained and replicate in HVS C488 ORF73-expressing cells.     

Biased Usages of Arginines and Lysines in Proteins Are Correlated with Local-Scale Fluctuations of the G + C Content of DNA Sequences

Amino acid residues arginine (R) and lysine (K) have similar physicochemical characteristics and are often mutually substituted during evolution without affecting protein function. Statistical examinations on human proteins show that more R than K residues are used in the proximity of R residues, whereas more K than R are used near K residues. This biased use occurs on both a global and a local scale (shorter than ∼100 residues). Even within a given exon, G + C-rich and A + T-rich short DNA segments preferentially encode R and K, respectively. The biased use of R and K on a local scale is also seen in Saccharomyces cerevisiae and Caenorhabdidtis elegans, which lack global-scale mosaic structures with varying GC%, or isochores. Besides R and K, several amino acids are also used with a positive or negative correlation with the local GC% of third codon bases. The local-, or ``within-gene'-, scale heterogeneity of the DNA sequence may influence the sequence of the encoded protein segment. Received: 2 March 1998 / Accepted: 23 April 1998     

Negative correlation of G+C content at silent substitution sites between orthologous human and mouse protein-coding sequences.

We conducted a genome-wide analysis of variations in guanine plus cytosine (G+C) content at the third codon position at silent substitution sites of orthologous human and mouse protein-coding nucleotide sequences. Alignments of 3776 human protein-coding DNA sequences with mouse orthologs having >50 synonymous codons were analyzed, and nucleotide substitutions were counted by comparing sequences in the alignments extracted from gap-free regions. The G+C content at silent sites in these pairs of genes showed a strong negative correlation (r = -0.93). Some gene pairs showed significant differences in G+C content at the third codon position at silent substitution sites. For example, human thymine-DNA glycosylase was A+T-rich at the silent substitution sites, while the orthologous mouse sequence was G+C-rich at the corresponding sites. In contrast, human matrix metalloproteinase 23B was G+C-rich at silent substitution sites, while the mouse ortholog was A+T-rich. We discuss possible implications of this significant negative correlation of G+C content at silent sites.     

Mitochondrial Genes Collectively Suggest the Paraphyly of Crustacea with Respect to Insecta

Complete sequences of seven protein coding genes from Penaeus notialis mitochondrial DNA were compared in base composition and codon usage with homologous genes from Artemia franciscana and four insects. The crustacean genes are significantly less A + T-rich than their counterpart in insects and the pattern of codon usage (ratio of G + C-rich versus A + T-rich codon) is less biased. A phylogenetic analysis using amino acid sequences of the seven corresponding polypeptides supports a sister-taxon status for mollusks–annelid and arthropods. Furthermore, a distance matrix-based tree and two most-parsimonious trees both suggest that crustaceans are paraphyletic with respect to insects. This is also supported by the inclusion of Panulirus argus COII (complete) and COI and COIII (partial) sequence data. From analysis of single and combined genes to infer phylogenies, it is observed that obtained from single genes are not well supported in most topologies cases and notably differ from that of the tree based on all seven genes. Received: 25 August 1998 / Accepted: 8 March 1999     

The mitochondrial genome of Ruspolia dubia (Orthoptera: Conocephalidae) contains a short A+T-rich region of 70 bp in length.

The complete sequence (14 971 bp) of the Ruspolia dubia mitochondrial genome was determined and annotated. The genome contains the gene content, base composition, and codon usage typical of metazoan mitochondrial genomes. All 37 genes are conserved in the positions observed most frequently in insect mitochondrial genome structures. The secondary structures of both small subunit and large subunit rRNA were predicted. The most unusual features found were the initiation codon (TTA) of COI and a short A+T-rich region of 70 bp in length. In addition, a short, highly conserved polythymidine stretch that was previously described in Orthoptera and Diptera was also present in the A+T-rich region.