Reduction of synonymous substitutions in the core protein gene of hepatitis C virus

Molecular evolutionary analyses were carried out to elucidate the phylogenetic relationships, the evolutionary rate, and the divergence times of hepatitis C viruses. Using the nucleotide sequences of the viruses isolated from various locations in the world, we constructed phylogenetic trees. The trees showed that strains isolated from a single location were not necessarily clustered as a group. This suggests that the viruses may be transferred with blood on a worldwide scale. We estimated the evolutionary rates at synonymous and nonsynonymous sites for all genes in the viral genome. We then found that the rate (1.35 × 10^–3 per site per year) at synonymous sites for the C gene was much smaller than those for the other genes (e.g., 6.29 × 10^–3 per site per year for the E gene). This indicates that a special type of functional constraint on synonymous substitutions may exist in the C gene. Because we found an open reading frame (ORF) with the C gene region, the possibility exists that synonymous substitutions for the C gene are constrained by the overlapping ORF whose reading frame is different from that of the C gene. Applying the evolutionary rates to the trees, we also suggest that major groups of hepatitis C viruses diverged from their common ancestor several hundred years ago. Correspondence to: T. Gojobori     

Avian influenza virus exhibits rapid evolutionary dynamics

Influenza A viruses from wild aquatic birds, their natural reservoir species, are thought to have reached a form of stasis, characterized by low rates of evolutionary change. We tested this hypothesis by estimating rates of nucleotide substitution in a diverse array of avian influenza viruses (AIV) and allowing for rate variation among lineages. The rates observed were extremely high, at >10(-3) substitutions per site, per year, with little difference among wild and domestic host species or viral subtypes and were similar to those seen in mammalian influenza A viruses. Influenza A virus therefore exhibits rapid evolutionary dynamics across its host range, consistent with a high background mutation rate and rapid replication. Using the same approach, we also estimated that the common ancestors of the hemagglutinin and neuraminidase sequences of AIV arose within the last 3,000 years, with most intrasubtype diversity emerging within the last 100 years and suggestive of a continual selective turnover.     

Extremely rapid spread of human immunodeficiency virus type 1 BF recombinants in Argentina

The epidemic of human immunodeficiency virus type 1 (HIV-1) in Argentina is distinctive in that many infections are caused by subtype BF recombinant viruses. To determine their demographic history, we estimated the evolutionary rate, mode of population growth, and age of genetic diversity among 40 BF vpu sequences. This revealed one of the highest substitution rates reported for HIV-1, at 10.793 x 10(-3) substitutions per site per year, and a very rapid rate of population growth, with an initial mean epidemic doubling time of 3.72 months. This rapid population growth is compatible with an elevated fitness for subtype BF compared to that for "pure" B and F viruses.     

Phylogenetic evidence for rapid rates of molecular evolution in the single-stranded DNA begomovirus tomato yellow leaf curl virus

Geminiviruses are devastating viruses of plants that possess single-stranded DNA (ssDNA) DNA genomes. Despite the importance of this class of phytopathogen, there have been no estimates of the rate of nucleotide substitution in the geminiviruses. We report here the evolutionary rate of the tomato yellow leaf curl disease-causing viruses, an intensively studied group of monopartite begomoviruses. Sequences from GenBank, isolated from diseased plants between 1988 and 2006, were analyzed using Bayesian coalescent methods. The mean genomic substitution rate was estimated to be 2.88 x 10(-4) nucleotide substitutions per site per year (subs/site/year), although this rate could be confounded by frequent recombination within Tomato yellow leaf curl virus genomes. A recombinant-free data set comprising the coat protein (V1) gene in isolation yielded a similar mean rate (4.63 x 10(-4) subs/site/year), validating the order of magnitude of genomic substitution rate for protein-coding regions. The intergenic region, which is known to be more variable, was found to evolve even more rapidly, with a mean substitution rate of approximately 1.56 x 10(-3) subs/site/year. Notably, these substitution rates, the first reported for a plant DNA virus, are in line with those estimated previously for mammalian ssDNA viruses and RNA viruses. Our results therefore suggest that the high evolutionary rate of the geminiviruses is not primarily due to frequent recombination and may explain their ability to emerge in novel hosts.     

Phylogenetic evidence for the rapid evolution of human B19 erythrovirus

Human B19 erythrovirus is a ubiquitous viral pathogen, commonly infecting individuals before adulthood. As with all autonomous parvoviruses, its small single-stranded DNA genome is replicated with host cell machinery. While the mechanism of parvovirus genome replication has been studied in detail, the rate at which B19 virus evolves is unknown. By inferring the phylogenetic history and evolutionary dynamics of temporally sampled B19 sequences, we observed a surprisingly high rate of evolutionary change, at approximately 10(-4) nucleotide substitutions per site per year. This rate is more typical of RNA viruses and suggests that high mutation rates are characteristic of the Parvoviridae.     

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.     

Sequence Heterogeneity in Closed Simian Virus 40 Deoxyribonucleic Acid

The heteroduplex molecules formed by self-annealing of denatured, singly nicked simian virus 40 (SV40) deoxyribonucleic acid (DNA) prepared from closed viral DNA were examined by formamide-protein film electron microscopy to test the DNA for sequence homogeneity. Sequence inhomogeneity appears in the heteroduplexes as single-strand loops. These result from sequence deletion or from sequence substitution, if regions greater than 50 nucleotides are involved. The undenatured DNA from viruses passaged twice at multiplicities of infection much less than 1 plaque-forming unit (PFU) per cell appeared to be homogeneous in size. The heteroduplexes formed by this DNA indicated that approximately 2% of the molecules carried deletions, but that substitutions were below the level of detection. In contrast, undenatured DNA from viruses grown by passaging undiluted lysates seven times or by infection with stock virus at a multiplicity of infection of 5 PFU per cell contained a large frequency of molecules shorter than the full length. The heteroduplex samples indicated that 12 and 7% of the undenatured material contained base substitutions, and 13 and 11% contained deletions. The deletions and substitutions appear to occur in separate molecules. Length measurements on heteroduplexes displaying the loop characteristic of substitutions have established that these molecules are from true sequence substitutions, and not from adjacent or overlapping deletions. More than 80% of the molecules carrying substitutions are shorter than the native SV40 length. On the average, the substituted sequence is about 20% of the length of SV40, but it replaces a sequence about 30% of the native length. The substituted sequences may be host cell nuclear DNA, possibly arising from integration of SV40 into the chromosome followed by excision of the SV40 DNA together with chromosomal DNA.     

Impact of human immunodeficiency virus type 1 gp41 amino acid substitutions selected during enfuvirtide treatment on gp41 binding and antiviral potency of enfuvirtide in vitro

Enfuvirtide (ENF), a novel human immunodeficiency virus type 1 (HIV-1) fusion inhibitor, has potent antiviral activity against HIV-1 both in vitro and in vivo. Resistance to ENF observed after in vitro passaging was associated with changes in a three-amino-acid (aa) motif, GIV, at positions 36 to 38 of gp41. Patients with ongoing viral replication while receiving ENF during clinical trials acquired substitutions within gp41 aa 36 to 45 in the first heptad repeat (HR-1) of gp41 in both population-based plasma virus sequences and proviral DNA sequences from isolates showing reduced susceptibilities to ENF. To investigate their impact on ENF susceptibility, substitutions were introduced into a modified pNL4-3 strain by site-directed mutagenesis, and the susceptibilities of mutant viruses and patient-derived isolates to ENF were tested. In general, susceptibility decreases for single substitutions were lower than those for double substitutions, and the levels of ENF resistance seen for clinical isolates were higher than those observed for the site-directed mutant viruses. The mechanism of ENF resistance was explored for a subset of the substitutions by expressing them in the context of a maltose binding protein chimera containing a portion of the gp41 ectodomain and measuring their binding affinity to fluorescein-labeled ENF. Changes in binding affinity for the mutant gp41 fusion proteins correlated with the ENF susceptibilities of viruses containing the same substitutions. The combined results support the key role of gp41 aa 36 to 45 in the development of resistance to ENF and illustrate that additional envelope regions contribute to the ENF susceptibility of fusion inhibitor-na?ve viruses and resistance to ENF.     

Comparison of an avian osteopetrosis virus with an avian lymphomatosis virus by RNA-DNA hybridization.

Myeloblastosis-associated virus (MAV)-2(0), a virus which was derived from avian myeloblastosis virus and induced a high incidence of osteopetrosis, was compared with avian lymphomatosis virus 5938, a recent field isolate which induced a high incidence of lymphomatosis. The following information was obtained. (i) MAV-2(0) induced osteopetrosis, nephroblastoma, and a very low incidence of hepatocellular carcinoma. No difference was seen in the oncogenic spectrum of end point and plaque-purified MAV-2(0). (ii) 125I-labeled RNA sequences from MAV-2(0) formed hybrids with DNA extracted from osteopetrotic bone at a rate suggesting five proviral copies per haploid cell genome. The extent of hybridization of MAV-2(0) RNA with DNA from osteopetrotic tissue was more extensive (87%) than was observed in reactions with DNA from uninfected chicken embryos (52%). (iii) Competition of unlabeled viral RNA in hybridization reactions between the radioactive RNA from the two viruses and their respective proviral sequences present in tumor tissues showed that 15 to 20% of the viral sequences detected in these reactions were unshared. In contrast, no differences were detected in competition analyses of RNA sequences from the two viruses detected in DNA of normal chicken cells. (iv) MAV-2(0) 35S RNA was indistinguishable in size from avian lymphomatosis virus 5938 35S RNA by polyacrylamide gel electrophoresis.     

Time scale of Poxvirus evolution

Unlike in vertebrates and RNA viruses, the molecular clock has not been estimated so far for DNA viruses. The extended conserved central region (102 kb) of the orthopoxvirus genome and the DNA polymerase gene (3 kb) were analyzed in viruses representing several genera of the family Poxviridae. Analysis was based on the known dating of the variola virus (VARV) transfer from Western Africa to South America and previous data on the phylogenetic relatedness of modern West African and South American isolates of VARV. The mutation accumulation rate was for the first time estimated for these DNA viruses at (0.9–1.2) × 10^?6 substitutions per site per year. It was assumed that poxviruses diverged from an ancestor approximately 500,000 years ago to form the recent species and that the ancestor of the genus Orthopoxvirus emerged approximately 300,000 years ago and gave origin to the modern species approximately 14,000 years ago.     

Rate of divergence of cellular sequences homologous to segments of Moloney sarcoma virus.

The RNA genome of the Moloney isolate of murine sarcoma virus (M-MSV) consists of two parts--a sarcoma-specific region with no homology to known leukemia viral RNAs, and a shared region present also in Moloney murine leukemia virus RNA. Complementary DNA was isolated which was specific for each part of the M-MSV genome. The DNA of a number of mammalian species was examined for the presence of nucleotide sequences homologous with the two M-MSV regions. Both sets of viral sequences had homologous nucleotide sequences present in normal mouse cellular DNA. MSV-specific sequences found in mouse cellular DNA closely matched those nucleotide sequences found in M-MSV as seen by comparisons of thermal denaturation profiles. In all normal mouse cells tested, the cellular set of M-MSV-specific nucleotide sequences was present in DNA as one to a few copies per cell. The rate of base substitution of M-MSV nucleotide sequences was compared with the rate of evolution of both unique sequences and the hemoglobin gene of various species. Conservation of MSV-specific nucleotide sequences among species was similar to that of mouse globin gene(s) and greater than that of average unique cellular sequences. In contrast, cellular nucleotide sequences that are homologous to the M-MSV-murine leukemia virus "common" nucleotide region were present in multiple copies in mouse cells and were less well matched, as seen by reduced melting profiles of the hybrids. The cellular common nucleotide sequences diverged very rapidly during evolution, with a base substitution rate similar to that reported for some primate and avian endogenous virogenes. The observation that two sets of covalently linked viral sequences evolved at very different rates suggests that the origin of M-MSV may be different from endogenous helper viruses and that cellular sequences homologous to MSV-specific nucleotide sequences may be important to survival.     

Gnarled-trunk evolutionary model of influenza A virus hemagglutinin

Human influenza A viruses undergo antigenic changes with gradual accumulation of amino acid substitutions on the hemagglutinin (HA) molecule. A strong antigenic mismatch between vaccine and epidemic strains often requires the replacement of influenza vaccines worldwide. To establish a practical model enabling us to predict the future direction of the influenza virus evolution, relative distances of amino acid sequences among past epidemic strains were analyzed by multidimensional scaling (MDS). We found that human influenza viruses have evolved along a gnarled evolutionary pathway with an approximately constant curvature in the MDS-constructed 3D space. The gnarled pathway indicated that evolution on the trunk favored multiple substitutions at the same amino acid positions on HA. The constant curvature was reasonably explained by assuming that the rate of amino acid substitutions varied from one position to another according to a gamma distribution. Furthermore, we utilized the estimated parameters of the gamma distribution to predict the amino acid substitutions on HA in subsequent years. Retrospective prediction tests for 12 years from 1997 to 2009 showed that 70% of actual amino acid substitutions were correctly predicted, and that 45% of predicted amino acid substitutions have been actually observed. Although it remains unsolved how to predict the exact timing of antigenic changes, the present results suggest that our model may have the potential to recognize emerging epidemic strains.     

Evolutionary pattern of the hemagglutinin gene of influenza B viruses isolated in Japan: cocirculating lineages in the same epidemic season.

The unexpectedly low efficacy of influenza vaccine during school outbreaks of influenza B virus in the spring of 1987 in Japan was probably attributable to a poor antibody response of vaccinees to the epidemic viruses. An antigenic analysis of the causative B viruses isolated in 1987 and 1988 showed much variation in hemagglutination inhibition patterns. The nucleotide sequences that code for the HA1 domain of B/Fukuoka/c-27/81, B/Ibaraki/2/85, B/Nagasaki/1/87, and B/Yamagata/16/88 viruses were determined and compared with those of the previously reported hemagglutinin genes. The nucleotide sequences of the hemagglutinin gene of a new variant, B/Yamagata/16/88, had only 93.4% homology with those of two other viruses from the same epidemic. An analysis of nucleotide and amino acid substitutions of the hemagglutinin genes of influenza B viruses revealed that new and some old variants could cocirculate in the same epidemic. A phylogenetic tree constructed by the neighbor-joining method allowed estimation of an evolutionary rate of 2.3 x 10(-3) synonymous (silent) substitutions per nucleotide site per year in the hemagglutinin gene.     

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.     

Analysis of the evolution and variation of the human influenza A virus nucleoprotein gene from 1933 to 1990.

This study examined the evolution and variation of the human influenza virus nucleoprotein gene from the earliest isolates to the present. Phylogenetic reconstruction of the most parsimonious evolutionary path connecting 49 nucleoprotein sequences yielded a single lineage. The average calculated rate of mutation was 3.6 nucleotide substitutions per year (2.3 x 10(-3) substitutions per site per year). Thirty-two percent of these mutations resulted in amino acid substitutions, and the remainder were silent mutations. Analysis of virus isolates from China and elsewhere showed no significant differences in their rate of evolution, genetic diversity, or mean survival time. The nearly constant rate of change was maintained through the two antigenic shifts, and there were no obvious changes in the number or types of mutations associated with the changes in the surface proteins. A detailed comparison of the changes that have occurred on the main evolutionary path with those that have occurred on the side branches of the phylogenetic tree was made. This showed that while 35% of the mutations on the side branches resulted in amino acid changes, only 21% of those on the main path affected the protein sequence. These results suggest that although the rate of change of the human influenza virus nucleoprotein is much higher than that previously described for avian influenza viruses, there are measurable constraints on the evolution of the surviving virus lineage. Comparison of the nucleoproteins of virus isolates adapted to chicken embryos with the nucleoproteins of those grown only in MDCK cells revealed no consistent differences between the virus pairs. Thus, although the nucleoprotein is known to be critical for host specificity, its adaptation to growth in eggs apparently involves no immediate selective pressures, such as are found with hemagglutinin.     

3D gene of foot-and-mouth disease virus. Conservation by convergence of average sequences

The nucleotide sequence of the 3D (polymerase) gene of eight epidemiologically related isolates of foot-and-mouth disease virus of serotype C1 is reported. The genetic heterogeneity of 3D RNA is compared with that of the VP1-coding RNA of the same viruses. Regression lines of substitutions per nucleotide that distinguish any pair of viruses as a function of the time interval between the corresponding isolations show: (1) the slope (substitutions/nucleotide per month) is 2.1 times larger for the VP1 RNA than for the 3D RNA region; (2) the intercept with the ordinate (substitutions/nucleotide) for VP1 RNA is indistinguishable from that for 3D RNA. Thus, the average heterogeneity of the VP1-coding region is very similar to that of the 3D-coding region only among co-circulating viruses. Nine mutations and points of heterogeneity occurred within nucleotide residues 883 to 1026, which encode an amino acid segment, extremely conserved among many different RNA viruses. The results suggest that, rather than due to inherently lower mutability, the conservation of 3D genes is caused by a limitation in the fixation of substitutions in viable genomes.     

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 origin and evolution of Ebola and Marburg viruses

Molecular evolutionary analyses for Ebola and Marburg viruses were conducted with the aim of elucidating evolutionary features of these viruses. In particular, the rate of nonsynonymous substitutions for the glycoprotein gene of Ebola virus was estimated to be, on the average, 3.6 x 10(-5) per site per year. Marburg virus was also suggested to be evolving at a similar rate. Those rates were a hundred times slower than those of retroviruses and human influenza A virus, but were of the same order of magnitude as that of the hepatitis B virus. When these rates were applied to the degree of sequence divergence, the divergence time between Ebola and Marburg viruses was estimated to be more than several thousand years ago. Moreover, most of the nucleotide substitutions were transitions and synonymous for Marburg virus. This suggests that purifying selection has operated on Marburg virus during evolution.