Nucleotide sequence of an infectious molecularly cloned genome of ground squirrel hepatitis virus.

We have determined the complete nucleotide sequence of an infectious cloned genome of ground squirrel hepatitis virus (GSHV), a nonpathogenic member of the hepadnavirus group. The genome is 3,311 base pairs long and contains the major open reading frames described for the related human and woodchuck hepatitis B viruses (HBV and WHV, respectively). These reading frames include genes for the major structural proteins (the surface and core antigens), unassigned open reading frames (A and B), the longer of which is presumed to encode the viral DNA polymerase, and an open reading frame preceding and continuous with the surface antigen gene. The arrangement of these open reading frames is similar to that encountered in the genomes of HBV and WHV: all of the reading frames are encoded on the same strand, they are positioned in the same fashion with respect to each other, and a large portion (at least 51%) of the genome can be translated in two reading frames. Comparisons of the predicted translational products of the three mammalian hepadnaviruses reveal 78% amino acid homology between the proteins of GSHV and WHV and 43% homology between those of GSHV and HBV. In addition, a perfect direct repeat of 10 to 11 base pairs, separated by ca. 46 to 223 base pairs, is present in the three mammalian viruses and in duck hepatitis B virus; the position of the repeats near the 5' termini of the two strands of virion DNA suggests a role in viral replication.     

Subtype ayw variant of hepatitis B virus. DNA primary structure analysis

The entire genome of human hepatitis B virus (HBV) occurring in Latvia was sequenced. This sequence, which is 3182 nucleotides long, was compared with the other previously published HBV genomes and was shown to share maximum homology with HBV subtype ayw DNA. The coordinates of 4 main open reading frames as well as hairpin structures are very well conserved in the two genomes. The distribution of nucleotide substitutions among different HBV genomes suggest that the open reading frames P and X can fulfil a coding function. On the basis of primary structure comparison for hepadnaviral DNAs several evolutionary conclusions can be drawn.     

A revised evolutionary history of hepatitis B virus (HBV)

Previous studies of the evolutionary history of hepatitis B virus (HBV) have been compromised by intergenotype recombination and complex patterns of nucleotide substitution, perhaps caused by differential selection pressures. We examined the phylogenetic distribution of recombination events among human HBV genotypes and found that genotypes A plus D, and genotypes B plus C, had distinct patterns of recombination suggesting differing epidemiological relationships among them. By analyzing the nonoverlapping regions of the viral genome we found strong bootstrap support for some intergenotypic groupings, with evidence of a division between human genotypes A–E from the viruses sampled from apes and human genotype F. However, the earliest events in the divergence of HBV remain uncertain. These uncertainties could not be explained by differential selection pressures, as the ratio of nonsynonymous-to-synonymous substitutions (d _N/d _S) did not vary extensively among lineages and there is no strong evidence for positive selection across the whole tree. Finally, we provide a new estimate of the mean substitution rate in HBV, 4.2 × 10⁻⁵, which suggests that divergence of HBV in humans and apes has occurred only in the last 6000 years.     

Annotation of selection strengths in viral genomes

MOTIVATION: Viral genomes tend to code in overlapping reading frames to maximize informational content. This may result in atypical codon bias and particular evolutionary constraints. Due to the fast mutation rate of viruses, there is additional strong evidence for varying selection between intra- and intergenomic regions. The presence of multiple coding regions complicates the concept of K(a)/K(s) ratio, and thus begs for an alternative approach when investigating selection strengths. Building on the paper by McCauley and Hein, we develop a method for annotating a viral genome coding in overlapping reading frames. We introduce an evolutionary model capable of accounting for varying levels of selection along the genome, and incorporate it into our prior single sequence HMM methodology, extending it now to a phylogenetic HMM. Given an alignment of several homologous viruses to a reference sequence, we may thus achieve an annotation both of coding regions as well as selection strengths, allowing us to investigate different selection patterns and hypotheses. RESULTS: We illustrate our method by applying it to a multiple alignment of four HIV2 sequences, as well as of three Hepatitis B sequences. We obtain an annotation of the coding regions, as well as a posterior probability for each site of the strength of selection acting on it. From this we may deduce the average posterior selection acting on the different genes. Whilst we are encouraged to see in HIV2, that the known to be conserved genes gag and pol are indeed annotated as such, we also discover several sites of less stringent negative selection within the env gene. To the best of our knowledge, we are the first to subsequently provide a full selection annotation of the Hepatitis B genome by explicitly modelling the evolution within overlapping reading frames, and not relying on simple K(a)/K(s) ratios.     

Human endogenous retrovirus K106 (HERV-K106) was infectious after the emergence of anatomically modern humans

HERV-K113 and HERV-K115 have been considered to be among the youngest HERVs because they are the only known full-length proviruses that are insertionally polymorphic and maintain the open reading frames of their coding genes. However, recent data suggest that HERV-K113 is at least 800,000 years old, and HERV-K115 even older. A systematic study of HERV-K HML2 members to identify HERVs that may have infected the human genome in the more recent evolutionary past is lacking. Therefore, we sought to determine how recently HERVs were exogenous and infectious by examining sequence variation in the long terminal repeat (LTR) regions of all full-length HERV-K loci. We used the traditional method of inter-LTR comparison to analyze all full length HERV-Ks and determined that two insertions, HERV-K106 and HERV-K116 have no differences between their 5' and 3' LTR sequences, suggesting that these insertions were endogenized in the recent evolutionary past. Among these insertions with no sequence differences between their LTR regions, HERV-K106 had the most intact viral sequence structure. Coalescent analysis of HERV-K106 3' LTR sequences representing 51 ethnically diverse individuals suggests that HERV-K106 integrated into the human germ line approximately 150,000 years ago, after the emergence of anatomically modern humans.     

Immune-induced evolutionary selection focused on a single reading frame in overlapping hepatitis B virus proteins

Viruses employ various means to evade immune detection. Reduction of CD8(+) T cell epitopes is one of the common strategies used for this purpose. Hepatitis B virus (HBV), a member of the Hepadnaviridae family, has four open reading frames, with about 50% overlap between the genes they encode. We computed the CD8(+) T cell epitope density within HBV proteins and the mutations within the epitopes. Our results suggest that HBV accumulates escape mutations that reduce the number of epitopes. These mutations are not equally distributed among genes and reading frames. While the highly expressed core and X proteins are selected to have low epitope density, polymerase, which is expressed at low levels, does not undergo the same selection. In overlapping regions, mutations in one protein-coding sequence also affect the other protein-coding sequence. We show that mutations lead to the removal of epitopes in X and surface proteins even at the expense of the addition of epitopes in polymerase. The total escape mutation rate for overlapping regions is lower than that for nonoverlapping regions. The lower epitope replacement rate for overlapping regions slows the evolutionary escape rate of these regions but leads to the accumulation of mutations more robust in the transfer between hosts, such as mutations preventing proteasomal cleavage into epitopes.     

Molecular evolution of the hepatitis B virus genome

The hepatitis B virus (HBV) has a circular DNA genome of about 3,200 base pairs. Economical use of the genome with overlapping reading frames may have led to severe constraints on nucleotide substitutions along the genome and to highly variable rates of substitution among nucleotide sites. Nucleotide sequences from 13 complete HBV genomes were compared to examine such variability of substitution rates among sites and to examine the phylogenetic relationships among the HBV variants. The maximum likelihood method was employed to fit models of DNA sequence evolution that can account for the complexity of the pattern of nucleotide substitution. Comparison of the models suggests that the rates of substitution are different in different genes and codon positions; for example, the third codon position changes at a rate over ten times higher than the second position. Furthermore, substantial variation of substitution rates was detected even after the effects of genes and codon positions were corrected; that is, rates are different at different sites of the same gene or at the same codon position. Such rates after the correction were also found to be positively correlated at adjacent sites, which indicated the existence of conserved and variable domains in the proteins encoded by the viral genome. A multiparameter model validates the earlier finding that the variation in nucleotide conservation is not random around the HBV genome. The test for the existence of a molecular clock suggests that substitution rates are more or less constant among lineages. The phylogenetic relationships among the viral variants were examined. Although the data do not seem to contain sufficient information to resolve the details of the phylogeny, it appears quite certain that the serotypes of the viral variants do not reflect their genetic relatedness. Correspondence to: Z. Yang     

Complete genome sequences of novel canine noroviruses in Hong Kong

We report the complete genome sequences of two novel isolates of norovirus isolated from the fecal swab specimens of dogs in Hong Kong. The complete viral genome is approximately 7.6 kb in length and consists of 3 overlapping open reading frames encoding the ORF1 polyprotein, VP1, and VP2, respectively. Analysis of the VP1 sequence suggested that these noroviruses are divergent from known noroviruses and may represent a novel phylogenetic clade within the genus.     

The genetic organization of integrated hepatitis B virus DNA in the human hepatoma cell line PLC/PRF/5.

Hepatitis B virus (HBV) DNA is often found integrated in the genome of infected human liver cells and is supposed to be related to the development of primary liver carcinoma (PLC). Four clones of HBV DNA-containing sequences derived from DNA of the human PLC-derived cell line PLC/PRF/5 are discussed. The viral sequences show no intricate rearrangements excepting for a duplication and an inversion in one case, and a deletion in another. In all cases integration of the viral DNA was seen to be in a region which is single-stranded in the unintegrated HBV DNA. Sequence homologies between human and viral DNA flanking the integration sites have been detected. That may have a functional role in integration. Nucleotide sequence analyses of regions encompassing the viral-human junctions reveal open reading frames which consist of viral and/or human information. The possible expression of chimeric or cellular proteins may play a role in tumour development, and offers directions for further investigations.     

Homologous recombination between different genotypes of hepatitis B virus

Phylogenetic analysis was used to examine the evolutionary relationships among 99 complete HBV sequences. Analysis revealed nine viral genomes clustered with different genotypes depending on genome region analyzed. This discordance indicated that recombination events occurred during HBV history. The putative breakpoints between genomes of different genotypes have been mapped. Six mosaic genomes representing B/C hybrids were isolated in East Asia and three A/D hybrids in Italy. At least some recombinant strains appear to be fully viable and possess high evolutionary potential. As a result, B/C recombinants overspread through the East Asia region. They were found among the isolates from Japan, China and Indonesia. Our results suggest that recombination is a significant and relatively frequent event in the evolution of HBV genome. A possible mechanism and the implications of recombination for the natural history of HBV, clinically important properties, and phylogenetic reconstruction are discussed.     

Recombination between sequences of hepatitis B virus from different genotypes

A comparison of 25 hepatitis B virus (HBV) isolates for which complete genome sequences are available revealed two that occupied different positions in phylogenetic trees reconstructed from different open reading frames. Further analysis indicated that this incongruence was the result of recombination between viruses of different genomic and antigenic types. Both putative recombinants originated from geographic regions where multiple genotypes are known to cocirculate. A search of the sequence databases showed evidence of similar intergenotypic recombinants. These observations indicate that recombination between divergent strains may represent an important source of genetic variation in HBV.     

The tendency of lentiviral open reading frames to become A-rich: constraints imposed by viral genome organization and cellular tRNA availability

Human immunodeficiency virus type 1 (HIV-1) and other lentiviridae demonstrate a strong preference for the A-nucleotide, which can account for up to 40% of the viral RNA genome. The biological mechanism responsible for this nucleotide bias is currently unknown. The increased A-content of these viral genomes corresponds to the typical use of synonymous codons by all members of the lentiviral family (HIV, SIV, BIV, FIV, CAEV, EIAV, visna) and the human spuma retrovirus, but not by other retroviruses like the human T-cell leukemia viruses HTLV-I and HTLV-II. In this article, we analyzed A-bias for all codon groups in all open reading frames of several lentiviruses. The extent of lentiviral codon bias could be related to host cellular translation. By calculating codon bias indices (CBIs), we were able to demonstrate an inverse correlation between the extent of codon bias and the rate of translation of individual reading frames in these viruses. Specifically, the shift toward A-rich codons is more pronounced in pol than in gag lentiviral genes. Since it is known that Gag synthesis exceeds Pol synthesis by a factor of 20 due to infrequent ribosomal frame-shifting during translation of the gap-pol mRNA molecule, we propose that the aminoacyl-tRNA availability in the host cell restricts the lentiviral preference for A-rich codons. In addition, less A-nucleotides were found in regions of the viral genome encoding multiple functions; e.g., overlapping reading frames (tat-rev-env) or in genes that overlap regulatory sequences (nef-LTR region). Finally, the characteristics of lentiviral codon usage are presented as a phylogenetic tree without the need for prior sequence alignment.Correspondence to: B. Berkhout     

Origin and evolution of avian microchromosomes

The origin of avian microchromosomes has long been the subject of much speculation and debate. Microchromosomes are a universal characteristic of all avian species and many reptilian karyotypes. The typical avian karyotype contains about 40 pairs of chromosomes and usually 30 pairs of small to tiny microchromosomes. This characteristic karyotype probably evolved 100-250 million years ago. Once the microchromosomes were thought to be a non-essential component of the avian genome. Recent work has shown that even though these chromosomes represent only 25% of the genome; they encode 50% of the genes. Contrary to popular belief, microchromosomes are present in a wide range of vertebrate classes, spanning 400-450 million years of evolutionary history. In this paper, comparative gene mapping between the genomes of chicken, human, mouse and zebrafish, has been used to investigate the origin and evolution of avian microchromosomes during this period. This analysis reveals evidence for four ancient syntenies conserved in fish, birds and mammals for over 400 million years. More than half, if not all, microchromosomes may represent ancestral syntenies and at least ten avian microchromosomes are the product of chromosome fission. Birds have one of the smallest genomes of any terrestrial vertebrate. This is likely to be the product of an evolutionary process that minimizes the DNA content (mostly through the number of repeats) and maximizes the recombination rate of microchromosomes. Through this process the properties (GC content, DNA and repeat content, gene density and recombination rate) of microchromosomes and macrochromosomes have diverged to create distinct chromosome types. An ancestral genome for birds likely had a small genome, low in repeats and a karyotype with microchromosomes. A "Fission-Fusion Model" of microchromosome evolution based on chromosome rearrangement and minimization of repeat content is discussed.     

Human cytomegalovirus clinical isolates carry at least 19 genes not found in laboratory strains.

Nucleotide sequence comparisons were performed on a highly heterogeneous region of three human cytomegalovirus strains, Toledo, Towne, and AD169. The low-passage, virulent Toledo genome contained a DNA segment of approximately 13 kbp that was not found in the Towne genome and a segment of approximately 15 kbp that was not found in the AD169 genome. The Towne strain contained approximately 4.7 kbp of DNA that was absent from the AD169 genome, and only about half of this segment was present, arranged in an inverted orientation, in the Toledo genome. These additional sequences were located at the unique long (UL)/b' (IRL) boundary within the L component of the viral genome. A region representing nucleotides 175082 to 178221 of the AD169 genome was conserved in all three strains; however, substantial reduction in the size of the adjacent b' sequence was found. The additional DNA segment within the Toledo genome contained 19 open reading frames not present in the AD169 genome. The additional DNA segment within the Towne genome contained four new open reading frames, only one of which shared homology with the Toledo genome. This comparison was extended to five additional clinical isolates, and the additional Toledo sequence was conserved in all. These findings reveal a dramatic level of genome sequence complexity that may explain the differences that these strains exhibit in virulence and tissue tropism. Although the additional sequences have not altered the predicted size of the viral genome (230 to 235 kbp), a total of 22 new open reading frames (denoted UL133 to UL154), many of which have sequence characteristics of glycoproteins, are now defined as cytomegalovirus specific. Our work suggests that wild-type virus carries more than 220 genes, some of which are lost by large-scale deletion and rearrangement of the UL/b' region during laboratory passage.     

Molecular evolution of the Paramyxoviridae and Rhabdoviridae multiple-protein-encoding P gene

Presented here is an analysis of the molecular evolutionary dynamics of the P gene among 76 representative sequences of the Paramyxoviridae and Rhabdoviridae RNA virus families. In a number of Paramyxoviridae taxa, as well as in vesicular stomatitis viruses of the Rhabdoviridae, the P gene encodes multiple proteins from a single genomic RNA sequence. These products include the phosphoprotein (P), as well as the C and V proteins. The complexity of the P gene makes it an intriguing locus to study from an evolutionary perspective. Amino acid sequence alignments of the proteins encoded at the P and N loci were used in independent phylogenetic reconstructions of the Paramyxoviridae and Rhabdoviridae families. P-gene-coding capacities were mapped onto the Paramyxoviridae phylogeny, and the most parsimonious path of multiple-coding-capacity evolution was determined. Levels of amino acid variation for Paramyxoviridae and Rhabdoviridae P-gene-encoded products were also analyzed. Proteins encoded in overlapping reading frames from the same nucleotides have different levels of amino acid variation. The nucleotide architecture that underlies the amino acid variation was determined in order to evaluate the role of selection in the evolution of the P gene overlapping reading frames. In every case, the evolution of one of the proteins encoded in the overlapping reading frames has been constrained by negative selection while the other has evolved more rapidly. The integrity of the overlapping reading frame that represents a derived state is generally maintained at the expense of the ancestral reading frame encoded by the same nucleotides. The evolution of such multicoding sequences is likely a response by RNA viruses to selective pressure to maximize genomic information content while maintaining small genome size. The ability to evolve such a complex genomic strategy is intimately related to the dynamics of the viral quasispecies, which allow enhanced exploration of the adaptive landscape.     

State of hepatitis B viral DNA in a human hepatoma cell line.

PLC/PRF/5, a tissue culture cell line isolated from a human hepatocellular carcinoma and producing hepatitis B surface antigen, was studied for the presence of hepatitis B virus (HBV)-specific DNA and RNA. PLC/PRF/5 cell DNA accelerated the rate of reassociation of HBV [32P]DNA, and quantitative experiments indicated that the cells contained approximately four copies of viral DNA per haploid, mammalian cell DNA equivalent. PLC/PRF/5 DNA accelerated the rate of reassociation of all individual restriction endonucleases HincII and HaeIII fragments of HBV [32P]DNA, indicating that DNA from all regions of the viral genome is present in the cells. This suggests that these cells contain at least most, and possibly all, of the viral genome. Digestion of PLC/PRF/5 cell DNA with restriction endonuclease HindIII (an enzyme found not to cleave the DNA of any HBV isolate so far examined) yielded only three fragments, all larger than virion DNA, which contained HBV DNA base sequences, suggesting that HBV DNA is integrated in high-molecular-weight DNA at three different sites in these cells and that there is no viral DNA in an episomal form. PLC/PRF/5 cell [32P]RNA was found to hybridize with all restriction fragments of HBV DNA adequately tested, indicating that at least most, and possibly all, of the viral DNA in these cells is transcribed.