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.     

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.     

Molecular epidemiology of the foot-and-mouth disease virus outbreak in the United Kingdom in 2001

The objective of this study was to quantify the extent to which the genetic diversity of foot-and-mouth disease virus (FMDV) arising over the course of infection of an individual animal becomes fixed, is transmitted to other animals, and thereby accumulates over the course of an outbreak. Complete consensus sequences of 23 genomes (each of 8,200 nucleotides) of FMDV were recovered directly from epithelium tissue acquired from 21 farms infected over a nearly 7-month period during the 2001 FMDV outbreak in the United Kingdom. An analysis of these consensus sequences revealed very few apparently ambiguous sites but clear evidence of 197 nucleotide substitutions at 191 different sites. We estimated the rate of nucleotide substitution to be 2.26 x 10(-5) per site per day (95% confidence interval [CI], 1.75 x 10(-5) to 2.80 x 10(-5)) and nucleotide substitutions to accrue in the consensus sequence at an average rate of 1.5 substitutions per farm infection. This is a sufficiently high rate showing that detailed histories of the transmission pathways can be reliably reconstructed. Coalescent methods indicated that the date at which FMDV first infected livestock in the United Kingdom was 7 February 2001 (95% CI, 20 January to 19 February 2001), which was identical to estimates obtained on the basis of purely clinical evidence. Nucleotide changes appeared to have occurred evenly across the genome, and within the open reading frame, the ratio of nonsynonymous-to-synonymous change was 0.09. The ability to recover particular transmission pathways of acutely acting RNA pathogens from genetic data will help resolve uncertainties about how virus is spread and could help in the control of future epidemics.     

Frequency of spontaneous mutations in an avian hepadnavirus infection

In this study, we measured the frequency of revertants of a cytopathic strain of the duck hepatitis B virus that bears a single nucleotide substitution in the pre-S envelope protein open reading frame, resulting in the amino acid substitution G133E. Cytopathic virus mixed with known amounts of a genetically marked wild-type virus was injected into ducklings. Virus outgrowth was accompanied by a coselection of wild-type and spontaneous revertants during recovery of the ducklings from the acute liver injury caused by death of the G133E-infected cells. The frequency of individual revertants in the selected noncytopathic virus population was estimated by determining the ratio of each revertant to the wild-type virus. Spontaneous revertants were found to be present at frequencies of 1 x 10(-5) to 6 x 10(-5) per G133E genome inoculated. A mathematical model was used to estimate that the mutation rate was 0.8 x 10(-5) to 4.5 x 10(-5) per nucleotide per generation.     

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.     

Founder virus population related to route of virus transmission: a determinant of intrahost human immunodeficiency virus type 1 evolution?

We and others have shown that in individual human immunodeficiency virus type 1 (HIV-1) infection, the adaptive evolution of HIV-1 is influenced by host immune competence. In this study, we tested the hypothesis that in addition to selective forces operating within the host, transmission bottlenecks have an impact on HIV-1 intrahost evolution. Therefore, we studied the intrahost evolution of the V3 region of the external glycoprotein gp120 of HIV-1 during the 3- and 5-year periods following seroconversion after parenteral versus sexual (male-to-male) transmission in 41 participants of the Amsterdam prospective cohorts of homosexual men (n = 31) and intravenous drug users (IVDUs; n = 10) who were AIDS free and had comparable numbers of CD4+ cells. We observed that HIV-1 strains in homosexual men accumulated over 5 years more nonsynonymous substitutions within the V3 loop than HIV-1 strains in IVDUs as a result of lower rates of nonsynonymous evolution in both the initial 3-year period from seroconversion and the following 2-year period as well as a larger proportion of nonsynonymous back substitutions in IVDUs. The mean numbers of synonymous substitutions did not differ between the two risk groups. Since HIV-1 strains in IVDUs could be distinguished from the viruses of homosexual men based on several nucleotide substitutions of which the most conserved is a synonymous substitution at the tip of the V3 loop (GGC pattern), we studied whether the founder virus population itself has an impact on the intrahost evolution of HIV-1. The mean number of nonsynonymous substitutions accumulated over 5 years within the V3 loop was lower in 10 IVDUs infected by the HIV-1 strains with the GGC signature than in 4 IVDUs infected by HIV-1 strains lacking this pattern, while the mean numbers of synonymous substitutions were similar in the two groups.     

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.     

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.     

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.     

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.     

Genotypes and phylogeographical relationships of infectious hematopoietic necrosis virus in California, USA

Infectious hematopoietic necrosis virus (IHNV) contains 3 major genogroups in North America with discreet geographic ranges designated as upper (U), middle (M), and lower (L). A comprehensive genotyping of 237 IHNV isolates from hatchery and wild salmonids in California revealed 25 different sequence types (a to y) all in the L genogroup; specifically, the genogroup contained 14 sequence types that were unique to individual isolates as well as 11 sequence types representing 2 or more identical isolates. The most evident trend was the phylogenetic and geographical division of the L genogroup into 2 distinct subgroups designated as LI and LII. Isolates within Subgroup LI were primarily found within waterways linked to southern Oregon and northern California coastal rivers. Isolates in Subgroup LII were concentrated within inland valley watersheds that included the Sacramento River, San Joaquin River, and their tributaries. The temporal and spatial patterns of virus occurrence suggested that infections among adult Chinook salmon in the hatchery or that spawn in the river are a major source of virus potentially infecting other migrating or resident salmonids in California. Serum neutralization results of the California isolates of IHNV corroborated a temporal trend of sequence divergence; specifically, 2 progressive shifts in which more recent virus isolates represent new serotypes. A comparison of the estimates of divergence rates for Subgroup LI (1 x 10(-5) mutations per nucleotide site per year) indicated stasis similar to that observed in the U genogroup, while the Subgroup LII rate (1 x 10(-3) mutations per nucleotide site per year) suggested a more active evolution similar to that of the M genogroup.     

Molecular evolution of a type 1 wild-vaccine poliovirus recombinant during widespread circulation in China

Type 1 wild-vaccine recombinant polioviruses were isolated from poliomyelitis patients in China from 1991 to 1993. We compared the sequences of 34 recombinant isolates over the 1,353-nucleotide (nt) genomic interval (nt 2480 to 3832) encoding the major capsid protein, VP1, and the protease, 2A. All recombinants had a 367-nt block of sequence (nt 3271 to 3637) derived from the Sabin 1 oral poliovirus vaccine strain spanning the 3'-terminal sequences of VP1 (115 nt) and the 5' half of 2A (252 nt). The remaining VP1 sequences were closely (up to 99.5%) related to those of a major genotype of wild type 1 poliovirus endemic to China up to 1994. In contrast, the non-vaccine-derived sequences at the 3' half of 2A were more distantly related (<90% nucleotide sequence match) to those of other contemporary wild polioviruses from China. The vaccine-derived sequences of the earliest (April 1991) isolates completely matched those of Sabin 1. Later isolates diverged from the early isolates primarily by accumulation of synonymous base substitutions (at a rate of approximately 3.7 x 10(-2) substitutions per synonymous site per year) over the entire VP1-2A interval. Distinct evolutionary lineages were found in different Chinese provinces. From the combined epidemiologic and evolutionary analyses, we propose that the recombinant virus arose during mixed infection of a single individual in northern China in early 1991 and that its progeny spread by multiple independent chains of transmission into some of the most populous areas of China within a year of the initiating infection.     

Molecular Evolution of Puumala Hantavirus

Puumala virus (PUUV) is a negative-stranded RNA virus in the genus Hantavirus, family Bunyaviridae. In this study, detailed phylogenetic analysis was performed on 42 complete S segment sequences of PUUV originated from several European countries, Russia, and Japan, the largest set available thus far for hantaviruses. The results show that PUUV sequences form seven distinct and well-supported genetic lineages; within these lineages, geographical clustering of genetic variants is observed. The overall phylogeny of PUUV is star-like, suggesting an early split of genetic lineages. The individual PUUV lineages appear to be independent, with the only exception to this being the Finnish and the Russian lineages that are closely connected to each other. Two strains of PUUV-like virus from Japan form the most ancestral lineage diverging from PUUV. Recombination points within the S segment were searched for and evidence for intralineage recombination events was seen in the Finnish, Russian, Danish, and Belgian lineages of PUUV. Molecular clock analysis showed that PUUV is a stable virus, evolving slowly at a rate of 0.7 x 10(-7) to 2.2 x 10(-6) nt substitutions per site per year.     

Rice yellow mottle virus, an RNA plant virus, evolves as rapidly as most RNA animal viruses

The rate of evolution of an RNA plant virus has never been estimated using temporally spaced sequence data, by contrast to the information available on an increasing range of animal viruses. Accordingly, the evolution rate of Rice yellow mottle virus (RYMV) was calculated from sequences of the coat protein gene of isolates collected from rice over a 40-year period in different parts of Africa. The evolution rate of RYMV was estimated by pairwise distance linear regression on five phylogeographically defined groups comprising a total of 135 isolates. It was further assessed from 253 isolates collected all over Africa by Bayesian coalescent methods under strict and relaxed molecular clock models and under constant size and skyline population genetic models. Consistent estimates of the evolution rate between 4 x 10(-4) and 8 x 10(-4) nucleotides (nt)/site/year were obtained whatever method and model were applied. The synonymous evolution rate was between 8 x 10(-4) and 11 x 10(-4) nt/site/year. The overall and synonymous evolution rates of RYMV were within the range of the rates of 50 RNA animal viruses, below the average but above the distribution median. Experimentally, in host change studies, substitutions accumulated at an even higher rate. The results show that an RNA plant virus such as RYMV evolves as rapidly as most RNA animal viruses. Knowledge of the molecular clock of plant viruses provides methods for testing a wide range of biological hypotheses.     

Development of an adult mouse model for studies on protection against rotavirus.

Although mice have been used as an animal model for studies on rotavirus disease, these studies have been limited by the short time period after birth during which mice are susceptible to rotavirus illness (i.e., approximately 15 days). To overcome this limitation, an adult mouse model was developed in which the endpoint was infection rather than illness. The model developed utilized a strain of mouse rotavirus (EDIM) adapted to grow in culture by multiple passages in MA104 cells. The second cell culture passage of EDIM caused severe diarrhea in neonatal BALB/c mice, and little or no amelioration of disease was observed after nine cell culture passages, even when this preparation was plaque purified. Oral administration of 2 x 10(3) PFU of passage 9 also consistently caused infection of mice 4, 10, 15, 30, 60, 120, and 180 days of age as determined by viral shedding and seroconversion. Reinoculation of these mice with the same virus preparation at 2, 3, or 4 months after the first inoculation produced no evidence of reinfection. In contrast, infection of neonatal mice with the heterotypic WC3 bovine rotavirus did not prevent reinfection with culture-adapted EDIM. Thus, this strain of EDIM caused consistent infection of previously uninoculated neonatal and adult BALB/c mice and produced homotypic but not heterotypic protection against reinfection.     

Evolution of Envelope Sequences from the Genital Tract and Peripheral Blood of Women Infected with Clade A Human Immunodeficiency Virus Type 1

The development of viral diversity during the course of human immunodeficiency virus type 1 (HIV-1) infection may significantly influence viral pathogenesis. The paradigm for HIV-1 evolution is based primarily on studies of male cohorts in which individuals were presumably infected with a single virus variant of subtype B HIV-1. In this study, we evaluated virus evolution based on sequence information of the V1, V2, and V3 portions of HIV-1 clade A envelope genes obtained from peripheral blood and cervical secretions of three women with genetically heterogeneous viral populations near seroconversion. At the first sample following seroconversion, the number of nonsynonymous substitutions per potential nonsynonymous site (dn) significantly exceeded substitutions at potential synonymous sites (ds) in plasma viral sequences from all individuals. Generally, values of dn remained higher than values of ds as sequences from blood or mucosa evolved. Mutations affected each of the three variable regions of the envelope gene differently; insertions and deletions dominated changes in V1, substitutions involving charged amino acids occurred in V2, and sequential replacement of amino acids over time at a small subset of positions distinguished V3. The relationship among envelope nucleotide sequences obtained from peripheral blood mononuclear cells, plasma, and cervical secretions was evaluated for each individual by both phylogenetic and phenetic analyses. In all subjects, sequences from within each tissue compartment were more closely related to each other than to sequences from other tissues (phylogenetic tissue compartmentalization). At time points after seroconversion in two individuals, there was also greater genetic identity among sequences from the same tissue compartment than among sequences from different tissue compartments (phenetic tissue compartmentalization). Over time, temporal phylogenetic and phenetic structure was detectable in mucosal and plasma viral samples from all three women, suggesting a continual process of migration of one or a few infected cells into each compartment followed by localized expansion and evolution of that population.     

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.     

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.     

The visna virus genome: evidence for a hypervariable site in the env gene and sequence homology among lentivirus envelope proteins.

The complete nucleotide sequence of the visna virus 1514 genome was determined. Our sequence confirms the relationship of visna virus and other lentiviruses to human immunodeficiency virus (HIV) both at the level of sequence homology and of genomic organization. Sequence homology is shown to extend to the transmembrane proteins of lentivirus env genes; this homology is strongest in the extracellular domain, suggesting that close structural and functional similarities may also exist among these envelope proteins. Comparison of our data with the sequence of visna virus LV1-1, an antigenic variant derived from strain 1514, demonstrates that the rate of divergence has been about 1.7 x 10(-3) substitutions per nucleotide per year in vivo. This rate is orders of magnitude higher than that for most DNA genomes, but agrees well with estimates of the rate for HIV. A statistically significant cluster of mutations in the env gene appears to represent a hypervariable site and may correspond to the epitope responsible for the antigenic differences between 1514 and LV1-1. Analysis of the potential RNA folding pattern of the visna virus env gene shows that this hypervariable site falls within a region with little potential for intramolecular base pairing. This correlation of hypervariability with lack of RNA secondary structure is strengthened by the fact that it also holds for a hypervariable site in the env gene of HIV.     

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.