Fixation of mutations in the viral genome during an outbreak of foot-and-mouth disease: heterogeneity and rate variations

Rates of fixation of mutations during the evolution of the foot-and-mouth disease virus (FMDV) C₁ in nature have been estimated by hybridization of viral RNA to cloned cDNAs representing defined FMDV genome segments, and comparison of the selected RNAs by T₁ RNase oligonucleotide fingerprinting. Values ranged from <0.04 × 10⁻² to 4.5 × 10⁻² substitutions per nucleotide per year (s/nt/yr), depending on the time period and the genomic segment considered. Rates for viral structural protein genes were up to sixfold higher than for nonstructural protein genes. Values in excess of 10⁻² s/nt/yr have been measured for the RNA region that encodes VP1–VP3. The nucleotide sequences of the major immunogenic region of capsid protein VP1 have been determined for six new FMDV C₁ isolates, and they are compared with the two previously known sequences of FMDV C₁ (C-S8 and C₁-O). Both oligonucleotide fingerprinting of selected RNA fragments and direct nucleotide sequencing demonstrate that genetic heterogeneity exists among three viruses isolated on the same day, introducing a significant indetermination in the evaluation of fixation rates of mutations. During the FMDV C₁ outbreak, amino acid substitutions did occur that are known to affect the immunological properties of the virus. The proportion of mutations between two viral RNAs does not increase significantly with the time elapsed between the two isolations, suggesting a cocirculation of multiple, related, nonidentical FMDVs (‘evolving quasispecies’) as the mode of evolution of this agent.     

Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing

Accurate identification of the transmitted virus and sequences evolving from it could be instrumental in elucidating the transmission of human immunodeficiency virus type 1 (HIV-1) and in developing vaccines, drugs, or microbicides to prevent infection. Here we describe an experimental approach to analyze HIV-1 env genes as intact genetic units amplified from plasma virion RNA by single-genome amplification (SGA), followed by direct sequencing of uncloned DNA amplicons. We show that this strategy precludes in vitro artifacts caused by Taq-induced nucleotide substitutions and template switching, provides an accurate representation of the env quasispecies in vivo, and has an overall error rate (including nucleotide misincorporation, insertion, and deletion) of less than 8 x 10(-5). Applying this method to the analysis of virus in plasma from 12 Zambian subjects from whom samples were obtained within 3 months of seroconversion, we show that transmitted or early founder viruses can be identified and that molecular pathways and rates of early env diversification can be defined. Specifically, we show that 8 of the 12 subjects were each infected by a single virus, while 4 others acquired more than one virus; that the rate of virus evolution in one subject during an 80-day period spanning seroconversion was 1.7 x 10(-5) substitutions per site per day; and that evidence of strong immunologic selection can be seen in Env and overlapping Rev sequences based on nonrandom accumulation of nonsynonymous mutations. We also compared the results of the SGA approach with those of more-conventional bulk PCR amplification methods performed on the same patient samples and found that the latter is associated with excessive rates of Taq-induced recombination, nucleotide misincorporation, template resampling, and cloning bias. These findings indicate that HIV-1 env genes, other viral genes, and even full-length viral genomes responsible for productive clinical infection can be identified by SGA analysis of plasma virus sampled at intervals typical in large-scale vaccine trials and that pathways of viral diversification and immune escape can be determined accurately.     

Evaluation of accumulation of hepatitis C virus mutations in a chronically infected chimpanzee: comparison of the core, E1, HVR1, and NS5b regions

Four hepatitis C virus genome regions (the core, E1, HVR1, and NS5b) were amplified and sequenced from yearly samples obtained from a chronically infected chimpanzee over a 12-year span. Nucleotide substitutions were found to accumulate in the core, E1, and HVR1 regions during the course of chronic infection; substitutions within the NS5b region were not detected for the first 8 years and were found to be minimal during the last 4 years. The rate of accumulation of mutations in the core and E1 regions, based on a direct comparison between the first 1979 sequence and the last 1990 sequence, was 1.120 x 10(-3), while phylogenetic ancestral comparison using the 12 yearly sequences showed a rate of 0.816 x 10(-3) bases per site per year. Temporal evaluation of the sequences revealed that there appeared to be periods in which substitutions accumulated and became fixed, followed by periods with relative stasis or random substitutions that did not persist. Synonymous and nonsynonymous substitutions within the core, E1, and HVR1 regions were also analyzed. In the core and E1 regions, synonymous substitutions predominated and gradually increased over time. However, within the HVR1 region, nonsynonymous substitutions predominated but gradually decreased over time.     

Concurrent evolution of human immunodeficiency virus type 1 in patients infected from the same source: rate of sequence change and low frequency of inactivating mutations.

Direct sequencing of segments of the envelope gene of human immunodeficiency virus type 1 proviruses in peripheral blood mononuclear cells has revealed that a cohort of hemophiliacs who were infected after exposure to a single common batch of factor VIII share closely related virus strains. Seventy-four sequences extending from hypervariable regions V4 through V5 from nine patients yielded a mean intrapatient nucleotide distance of 5.5%, while a mean of 4.2% was observed in 39 sequences of the V3 loop (six patients). Phylogenetic analysis revealed that sequences of six Edinburgh patients were particularly closely related and those from a patient infected in the United States were very distinct. The mean nucleotide distance among these six was 8.3%, while the mean distance from the U.S.-derived sequences was 25.5% in the V4-V5 region. The rate of sequence change across this patient group has been estimated to be 0.4% per year in the V4-V5 region and 0.5% per year in the V3 region, with at least a twofold range across patients. Only two inactivating nucleotide substitutions have been observed in a total of 42 kb of sequence obtained from the env and gag genes during this study.     

Antigenic and genetic variation in influenza A (H1N1) virus isolates recovered from a persistently infected immunodeficient child.

Antigenic and genetic variations have been analyzed in eight consecutive isolates recovered from a child with severe combined immunodeficiency syndrome persistently infected with naturally acquired type A (H1N1) influenza virus over a 10-month period. Hemagglutination inhibition reactions and T1 oligonucleotide fingerprinting demonstrated that these viruses were related to strains causing outbreaks in the United States at that time (1983 to 1984) but that antigenic and genetic differences between consecutive isolates could be detected. This variation between isolates was examined further by sequencing the RNAs encoding the HA1 region of the hemagglutinin (HA) and the nucleoprotein (NP) in five of the consecutive isolates. Multiple point mutations were detected in both genes, and a deletion of one amino acid was detected in the HA. Depending on the isolates compared, 5.8 x 10(-3) to 17 x 10(-3) substitutions per nucleotide site per year were detected in the RNAs encoding the HA1, and 3.5 x 10(-3) to 24 x 10(-3) substitutions per nucleotide site per year were detected in the NP gene. Fifty-four percent of the base changes in the HA1 and 73% in the NP led to amino acid substitutions. A progressive accumulation of mutations over time was not observed, suggesting that the genetic diversity of these viruses may best be interpreted as the result of shifts in the population equilibrium (quasi-species) of replicating variant genomes.     

Rapid selection of genetic and antigenic variants of foot-and-mouth disease virus during persistence in cattle.

Rapid evolution of foot-and-mouth disease virus (FMDV) is documented during persistent infections of cattle. The carrier state was established experimentally with plaque-purified FMDV of serotype C3. Virus was recovered from the esophageal pharyngeal area of the animals up to 539 days postinfection. Analysis of capsid proteins by electrofocusing and by electrophoretic mobility of the genomic poly(C)-rich tract suggested heterogeneity in several isolates and sequential dominance of viral subpopulations. Nucleotide sequences of the VP1-coding region of the parental FMDV C3 clones and of seven isolates from the carrier cattle showed point mutations that represented rates of fixation of mutations of 0.9 X 10(-2) to 7.4 X 10(-2) substitutions per nucleotide per year; 59% of the base changes led to amino acid substitutions, some of which were located within residues 135 to 151, a region involved in neutralization of FMDV. In the esophageal pharyngeal fluid samples, FMDV C3-neutralizing activity was present. Antigenic variation was demonstrated with monoclonal antibodies raised against FMDV C3. Two isolates from carrier cattle differed from the parental virus by 10(2)- or 10(3)-fold decreased reactivity with neutralizing monoclonal antibodies. We suggest that persistent, inapparent infections of ruminants, in addition to being a reservoir of virus, may promote the rapid selection of antigenically variant FMDVs.     

Sequencing and analysis for the full-length genome RNA of foot-and-mouth disease virus China/99

The complete nucleotide sequence of genomic RNA of foot and mouth disease virus (FMDV) strain China/99 from infected bovine tongue epithelium is presented. The nucleotide sequence extending from the 5' end of the genomic RNA to the 5' end of poly (A) tail contains 8173 nucleotides (nt). Its open reading frame, which encodes a single polypeptide of 2332 amino acids, encompasses 6999 nt starting from the initiation codon AUG and terminating at the UAA codon 93 bases upstream from the 5' end of poly (A) tract. The 5' untranslated region (UTR) is composed of 1081 nt. The consensus of the 1d gene of FMDV strain China/99 compared with that of UKG/6/2001, UKG/12/2001, China/99HN4 and China/3/Tibet is over 97%. The result showed the stains belong to the members of the Pan-Asia family. There is a remarkable differentiation in the function-unknown (FUR), p2 and p3 regions between FMDV isolates from infected cattle and swine, especially in 3a gene. No deletion was found in genes /, 1a, 1b, 2a, 2c, 3b, and 3d. Thes     

Dynamics of mutation and recombination in a replicating population of complementing, defective viral genomes

In a previous study, we documented that serial passage of a biological clone of foot-and-mouth disease virus (FMDV) at high multiplicity of infection (moi) in cell culture resulted in viral populations dominated by defective genomes that included internal in-frame deletions, affecting the L and capsid-coding regions, and were infectious by complementation. In the present study, analyses of the defective genomes present in individual viral plaques, and of consensus nucleotide sequences determined for the entire genomes of sequential samples, have revealed a continuous dynamics of mutation and recombination. At some points of high genetic instability, multiple minority genomes with different internal deletions co-existed in the population. At later passages, a new defective RNA arose and displaced a related, previously dominant RNA. Nucleotide sequences of the different genomic forms found in sequential isolates have revealed an accumulation of mutations at an average rate of 0.12 substitutions per genome per passage. At the regions around the deletion sites, substantial, minor or no nucleotide sequence identity is found, suggesting relaxed sequence requirements for the occurrence of internal deletions. Competition experiments indicate a selective advantage of late phase defective genomes over their precursor forms. The defective genome-based FMDV retained an expansion of host cell tropism, undergone by the standard virus at a previous stage of the same evolutionary lineage. Thus, despite a complex dynamics of mutation and recombination, and phases of high genetic instability, a biologically relevant phenotypic trait was stably maintained after the evolutionary transition towards a primitive genome segmentation. The results extend the concept of a complex spectrum of mutant genomes to a complex spectrum of defective genomes in some evolutionary transitions of RNA viruses.     

Beyond the consensus: dissecting within-host viral population diversity of foot-and-mouth disease virus by using next-generation genome sequencing

The diverse sequences of viral populations within individual hosts are the starting material for selection and subsequent evolution of RNA viruses such as foot-and-mouth disease virus (FMDV). Using next-generation sequencing (NGS) performed on a Genome Analyzer platform (Illumina), this study compared the viral populations within two bovine epithelial samples (foot lesions) from a single animal with the inoculum used to initiate experimental infection. Genomic sequences were determined in duplicate sequencing runs, and the consensus sequence of the inoculum determined by NGS was identical to that previously determined using the Sanger method. However, NGS revealed the fine polymorphic substructure of the viral population, from nucleotide variants present at just below 50% frequency to those present at fractions of 1%. Some of the higher-frequency polymorphisms identified encoded changes within codons associated with heparan sulfate binding and were present in both foot lesions, revealing intermediate stages in the evolution of a tissue culture-adapted virus replicating within a mammalian host. We identified 2,622, 1,434, and 1,703 polymorphisms in the inoculum and in the two foot lesions, respectively: most of the substitutions occurred in only a small fraction of the population and represented the progeny from recent cellular replication prior to onset of any selective pressures. We estimated the upper limit for the genome-wide mutation rate of the virus within a cell to be 7.8 × 10(-4) per nucleotide. The greater depth of detection achieved by NGS demonstrates that this method is a powerful and valuable tool for the dissection of FMDV populations within hosts.     

Molecular evolution of the major capsid protein VP1 of enterovirus 70.

Nucleotide sequences of the genome RNA encoding capsid protein VP1 (918 nucleotides) of 18 enterovirus 70 (EV70) isolates collected from various parts of the world in 1971 to 1981 were determined, and nucleotide substitutions among them were studied. The genetic distances between isolates were calculated by the pairwise comparison of nucleotide difference. Regression analysis of the genetic distances against time of isolation of the strains showed that the synonymous substitution rate was very high at 21.53 x 10(-3) substitution per nucleotide per year, while the nonsynonymous rate was extremely low at 0.32 x 10(-3) substitution per nucleotide per year. The rate estimated by the average value of synonymous and nonsynonymous substitutions (W.-H. Li, C.-C. Wu, and C.-C. Luo, Mol. Biol. Evol. 2:150-174, 1985) was 5.00 x 10(-3) substitution per nucleotide per year. Taking the average value of synonymous and nonsynonymous substitutions as genetic distances between isolates, the phylogenetic tree was inferred by the unweighted pairwise grouping method of arithmetic average and by the neighbor-joining method. The tree indicated that the virus had evolved from one focal place, and the time of emergence was estimated to be August 1967 +/- 15 months, 2 years before first recognition of the pandemic of acute hemorrhagic conjunctivitis. By superimposing every nucleotide substitution on the branches of the phylogenetic tree, we analyzed nucleotide substitution patterns of EV70 genome RNA. In synonymous substitutions, the proportion of transitions, i.e., C<==>U and G<==>A, was found to be extremely frequent in comparison with that reported on other viruses or pseudogenes. In addition, parallel substitutions (independent substitutions at the same nucleotide position on different branches, i.e., different isolates, of the tree) were frequently found in both synonymous and nonsynonymous substitutions. These frequent parallel substitutions and the low nonsynonymous substitution rate despite the very high synonymous substitution rate described above imply a strong restriction on nonsynonymous substitution sites of VP1, probably due to the requirement for maintaining the rigid icosahedral conformation of the virus.     

The low evolutionary rate of human T-cell lymphotropic virus type-1 confirmed by analysis of vertical transmission chains

The evolutionary rate of the human T-cell lymphotropic virus type-1 (HTLV-1) is considered to be very low, in strong contrast to the related human retrovirus HIV. However, current estimates of the HTLV-1 rate rely on the anthropological calibration of phylogenies using assumed dates of human migration events. To obtain an independent rate estimate, we analyzed two variable regions of the HTLV-1 genome (LTR and env) from eight infected families. Remarkable genetic stability was observed, as only two mutations in LTR (756 bp) and three mutations in env (522 bp) occurred within the 16 vertical transmission chains, including one ambiguous position in each region. The evolutionary rate in HTLV-1 was then calculated using a maximum-likelihood approach that used the highest and lowest possible times of HTLV-1 shared ancestry, given the known transmission histories. The rates for the LTR and env regions were 9.58 x 10(-8)-1.25 x 10(-5) and 7.84 x 10(-7) -2.33 x 10(-5)nucleotide substitutions per site per year, respectively. A more precise estimate was obtained for the combined LTR-env data set, which was 7.06 x 10(-7)-1.38 x 10(-5)substitutions per site per year. We also note an interesting correlation between the occurrence of mutations in HTLV-1 and the age of the individual infected.     

Intrahost human immunodeficiency virus type 1 evolution is related to length of the immunocompetent period.

The antigenic diversity threshold theory predicts that antigenic sites of human immunodeficiency virus type 1, such as the V3 region of the external glycoprotein gp120, evolve more rapidly during the symptom-free period in individuals progressing to AIDS than in those who remain asymptomatic for a long time. To test this hypothesis, genomic RNA sequences were obtained from the sera of 44 individuals at seroconversion and 5 years later. The mean number of nonsynonymous nucleotide substitutions in the V3 region of the viruses circulating in 31 nonprogressors (1.1 x 10(-2) +/- 0.1 x 10(-2) per site per year) was higher than the corresponding value for 13 progressors (0.66 x 10(-2) +/- 0.1 x 10(-2) per site per year) (P < 0.01), while no difference between the mean numbers of synonymous substitutions in the two groups was seen (0.37 x 10(-2) +/- 0.1 x 10(-2) and 0.51 x 10(-2) +/- 0.2 x 10(-2) per site per year for nonprogressors and progressors, respectively; P > 0.1). The mean ratios of synonymous nucleotide p distance to nonsynonymous p distance were 0.35 for nonprogressors and 0.62 for progressors. The number of nonsynonymous substitutions was not associated with virus load or virus phenotype, which are established predictors of disease progression, but correlated strongly with the duration of the immunocompetent period (r2 = 0.41; P = 0.001). This indicates that there is no causative relationship between intrahost evolution and CD4+ cell decline. Our data suggest that intrahost evolution in human immunodeficiency virus type 1 infection is driven by selective forces, the strength of which is related to the duration of the immunocompetent period.     

Fixation of mutations at the VP1 gene of foot-and-mouth disease virus. Can quasispecies define a transient molecular clock?

The number of nucleotide (nt) substitutions found in the VP1 gene (encoding viral capsid protein) between any two of 16 closely related isolates of foot-and-mouth disease virus (FMDV) has been quantified as a function of the time interval between isolations [Villaverde et al., J. Mol. Biol. 204 (1988) 771-776]. One of them (isolate C-S12) includes some replacements found in isolates that preceded it and other replacements found in later isolates. The study has revealed alternating periods of rapid evolution and of relative genetic stability of VP1. During a defined period of acute disease, the rate of fixation of replacements at the VP1 coding segment was 6 x 10(-3) substitutions per nt per year. Only small differences in the rate of evolution were observed between subsegments within the VP1 gene. The observation of a relatively constant rate of evolution during a disease episode was unexpected. We propose that such constancy may be a consequence of random sampling of mutants from the FMDV quasispecies, followed by their amplification in susceptible hosts (to generate a new quasispecies). Successive sampling and amplification events may result in a steady accumulation of mutations.     

Human glyceraldehyde-3-phosphate dehydrogenase pseudogenes: Molecular evolution and a possible mechanism for amplification

We screened two human genomic libraries and isolated 14 different clones, designated λG1 and EG1-EG13, homologous to human glyceraldehyde-3-phosphate dehydrogenase (GAPD) cDNA. Subcloning and sequencing these recombinant phages led us to classify them as five different pseudogenes (ψG1–ψG5). All these sequences show such features typical of processed pseudogenes as numerous mutations, insertions, and deletions. The identity of numerous mutated sites among these pseudogenes and the presence of two Alu sequences flanking both ends of ψG1 suggest that GAPD pseudogenes originated from a unique reverse transcribed mRNA followed by gene duplication. The rate of nucleotide substitutions per site per year for known GAPD functional genes is low both for the synonymous substitutions (1.87×10⁻⁹) and for the nonsynonymous substitutions (0.12￠10⁻⁹) and indicates that the GAPD cDNA sequence is well conserved not only at the amino acid level, but also at the nucleotide level. The rate of nucleotide substitutions per site per year for GAPD pseudogenes shows a higher value (5.9×10⁻⁹) and suggests that these pseudogenes do not have any functional role. This work was supported by grants from the Consiglio Nazionale delle Ricerche and the Ministero Pubblica Istruzione (Rome, Italy). Special acknowledgment is given to the “Progetto Finalizzato Ingegneria Genetica e Basi Molecolari delle Malattie Ereditarie.”     

Molecular evolution of the tick-borne encephalitis and Powassan viruses

The problem of emerging viruses, their genetic diversity and viral evolution in nature are attracting more attention. The phylogenetic analysis and evaluationary rate estimation were made for pathogenic flaviviruses such as tick-borne encephalitis virus (TBEV) and Powassan (PV) circulated in natural foci in Russia. 47 nucleotide sequences of encoded protein E of the TBEV and 17 sequences of NS5 genome region of the PV have been used. It was found that the rate of accumulation of nucleotide substitutions for E genome region of TBEV was approximately 1.4 x 10(-4) and 5.4 x 10(-5) substitutions per site per year for NS5 genome region of PV. The ratio of non-synonymous nucleotide substitutions to synonymous substitution (dN/dS) for viral sequences were estimated of 0.049 for TBEV and 0.098 for PV. Maximum value dN/dS was 0.201-0.220 for sub-cluster of Russian and Canadian strains of PV and the minimum - 0.024 for cluster of Russian and Chinese strains of Far Eastern genotype TBEV. Evaluation of time intervals of evolutionary events associated with these viruses showed that European subtype TBEV are diverged from all-TBEV ancestor within approximately 2750 years and the Siberian and Far Eastern subtypes are emerged about 2250 years ago. The PV was introduced into natural foci of the Primorsky Krai of Russia only about 70 years ago and PV is a very close to Canadian strains of PV. Evolutionary picture for PV in North America is similar to evolution of Siberian and Far Eastern subtypes TBEV in Asia. The divergence time for main genetic groups of TBEV and PV are correlated with historical periods of warming and cooling. These allow to propose a hypothesis that climate changes were essential to the evolution of the flaviviruses in the past millenniums.     

Sequencing and analysis for the full-length genome RNA of foot-and-mouth disease virus China/99

The complete nucleotide sequence of genomic RNA of foot and mouth disease virus (FMDV) strain China/99 from infected bovine tongue epithelium is presented. The nucleotide sequence extending from the 5’ end of the genomic RNA to the 5’ end of poly (A) tail contains 8173 nucleotides (nt). Its open reading frame, which encodes a single polypeptide of 2332 amino acids, encompasses 6999 nt starting from the initiation codon AUG and terminating at the UAA codon 93 bases upstream from the 5’ end of poly (A) tract. The 5’ untranslated region (UTR) is composed of 1081nt. The consensus of the 1d gene of FMDV strain China/99 compared with that of UKG/6/2001, UKG/12/2001, China/99HN4 and China/3/Tibet is over 97%. The result showed the stains belong to the members of the Pan-Asia family. There is a remarkable differentiation in the function-unknown (FUR), p2 and p3 regions between FMDV isolates from infected cattle and swine, especially in 3a gene. No deletion was found in genes /, 1a, 1b, 2a, 2c, 3b, and 3d. These genes might be indispensable to the surviving of FMDV. The secondary structures of small (S) fragments, FUR and an internal ribosome entry site can be classified into three types, and the S fragment and 3’ UTR of the positive-sense RNA fold into stem-loop structures similar to the shape of clover.     

Evolution of alphaviruses in the eastern equine encephalomyelitis complex.

Evolution of viruses in the eastern equine encephalomyelitis (EEE) complex was studied by analyzing RNA sequences and oligonucleotide fingerprints from isolates representing the North and South American antigenic varieties. By using homologous sequences of Venezuelan equine encephalomyelitis virus as an outgroup, phylogenetic trees revealed three main EEE virus monophyletic groups. A North American variety group included all isolates from North America and the Caribbean. One South American variety group included isolates from the Amazon basin in Brazil and Peru, while the other included strains from Argentina, Guyana, Ecuador, Panama, Trinidad, and Venezuela. No evidence of heterologous recombination was obtained when three separate regions of the EEE virus genome were analyzed independently. Estimates of the overall rate of EEE virus evolution (nucleotide substitution) were 1.6 x 10(-4) substitution per nucleotide per year for the North American group and 4.3 x 10(-4) for the Argentina-Panama South American group. Evolutionary rate estimates for the North American group increased over 10-fold (from about 2 x 10(-5) to 4 x 10(-4)) concurrent with divergence of two monophyletic groups during the early 1970s. The North and South American antigenic varieties diverged roughly 1,000 years ago, while the two main South American groups diverged about 450 years ago. Analysis of multiple strains isolated from an upstate New York transmission focus during the same years suggested that, in certain locations, EEE virus may be relatively isolated for short time periods.     

Selection of genetic variants of simian immunodeficiency virus in persistently infected rhesus monkeys.

Genetic and antigenic variation may be one means by which lentiviruses that cause AIDS avoid elimination by host immune responses. Genetic variation in the envelope gene (env) was studied by comparing the nucleotide sequences of 27 clones obtained from two rhesus monkeys infected with molecularly cloned simian immunodeficiency virus. All 27 clones differed from each other and differed from the input clone in the gp120 (SU) portion of the envelope gene. Nucleotide substitutions were shown to accumulate with time at an average rate of 8.5 per 1,000 per year in SU. Surprisingly, the majority of nucleotide substitutions (81%) resulted in amino acid changes. Variation in SU was not random but occurred predominantly in five discrete regions. Within these variable regions, a remarkable 98% of the nucleotide substitutions changed the amino acid. These results demonstrate that extensive sequence variability accumulates in vivo after infection with molecularly cloned virus and that selection occurs in vivo for changes in distinct variable regions in env.     

Molecular evolution of hepatitis B virus over 25 years

Determining the longitudinal molecular evolution of hepatitis B virus (HBV) is difficult due to HBV's genomic complexity and the need to study paired samples collected over long periods of time. In this study, serial samples were collected from eight hepatitis B virus e antigen-negative asymptomatic carriers of HBV genotype B in 1979 and 2004, thus providing a 25-year period to document the long-term molecular evolution of HBV. The rate, nature, and distribution of mutations that emerged over 25 years were determined by phylogenetic and linear regression analysis of full-length HBV genome sequences. Nucleotide hypervariability was observed within the polymerase and pre-S/S overlap region and within the core gene. The calculated mean number of nucleotide substitutions/site/year (7.9 x 10(-5)) was slightly higher than published estimates (1.5 x 10(-5) to 5 x 10(-5)). Nucleotide changes in the quasispecies population did not significantly alter the molecular evolutionary rate, based on linear regression analysis of evolutionary distances among serial clone pre-S region sequences. Therefore, the directly amplified or dominant sequence was sufficient to estimate the putative molecular evolutionary rate for these long-term serial samples. On average, the ratio of synonymous (dS) to nonsynonymous (dN) substitutions was highest for the polymerase-coding region and lowest for the core-coding region. The low dS/dN ratios observed within the core suggest that selection occurs within this gene region, possibly as an immune evasion strategy. The results of this study suggest that HBV sequence divergence may occur more rapidly than previously estimated, in a host immune phase-dependent manner.