Rates of evolution of the retroviral oncogene of Moloney murine sarcoma virus and of its cellular homologues

A bioactive phlebovirus-like envelope protein in a hookworm endogenous virus

Endogenous viral elements (EVEs), accounting for 15% of our genome, serve as a genetic reservoir from which new genes can emerge. Nematode EVEs are particularly diverse and informative of virus evolution. We identify Atlas virus-an intact retrovirus-like EVE in the human hookworm Ancylostoma ceylanicum, with an envelope protein genetically related to GN-GC glycoproteins from the family Phenuiviridae. A cryo-EM structure of Atlas GC reveals a class II viral membrane fusion protein fold not previously seen in retroviruses. Atlas GC has the structural hallmarks of an active fusogen. Atlas GC trimers insert into membranes with endosomal lipid compositions and low pH. When expressed on the plasma membrane, Atlas GC has cell-cell fusion activity. With its preserved biological activities, Atlas GC has the potential to acquire a cellular function. Our work reveals structural plasticity in reverse-transcribing RNA viruses.

An open chat with…Takashi Gojobori

Takashi Gojobori was one of the original founding members of the Editorial Board of FEBS Open Bio and is easily one of our most hard‐working editors, having carefully evaluated hundreds of manuscripts on bioinformatics for us over the last decade. Takashi Gojobori received his doctorate at Kyushu University, Japan, before joining the University of Texas Health Science Center as a Research Associate and then Research Assistant Professor. Takashi was formerly the Vice‐Director and Professor of National Institute of Genetics (NIG), Mishima, Japan, and at present, he is Distinguished Professor of Bioscience and Bioscience Acting Director of the Computational Bioscience Research Center at King Abdullah University of Science and Technology (KAUST), Saudi Arabia. Takashi was formerly the Editor‐in‐Chief of the journal Gene and is also a long‐term member of the Editorial Board of FEBS Letters, among other editorial appointments. In honour of his 10 years (and counting!) of service to the journal, we present here this interview with him on his research and experiences.

Long-term virus evolution in nature

Viruses spread to give rise to epidemics and pandemics, and some key parameters that include virus and host population numbers determine virus persistence or extinction in nature. Viruses evolve at different rates depending on the polymerase copying fidelity during genome replication and a number of environmental influences. Calculated rates of evolution in nature vary depending on the time interval between virus isolations. In particular, intrahost evolution is generally more rapid that interhost evolution, and several possible mechanisms for this difference are considered. The mechanisms by which the error-prone viruses evolve are very unlikely to render the operation of a molecular clock (constant rate of incorporation of mutations in the evolving genomes), although a clock is assumed in many calculations. Several computational tools permit the alignment of viral sequences and the establishment of phylogenetic relationships among viruses. The evolution of the virus in the form of dynamic mutant clouds in each infected individual, together with multiple environmental parameters renders the emergence and reemergence of viral pathogens an unpredictable event, another facet of biological complexity.

Possible Roles of New Mutations Shared by Asian and American Zika Viruses

Originating in Africa, the Zika virus (ZIKV) has spread to Asia, Pacific Islands and now to the Americas and beyond. Since the first isolation in 1947, ZIKV strains have been sampled at various times in the last 69 years, but this history has not been reflected in studying the patterns of mutation accumulation in their genomes. Implementing the viral history, we show that the ZIKV ancestor appeared sometime in 1930–1945 and, at that point, its mutation rate was probably less than 0.2 × 10⁻³/nucleotide site/year and subsequently increased significantly in most of its descendants. Sustaining a high mutation rate of 4 × 10⁻³/site/year throughout its evolution, the Ancestral Asian strain, which was sampled from a mosquito in Malaysia, accumulated 13 mutations in the 3′-untranslated region of RNA stem-loops prior to 1963, seven of which generate more stable stem-loop structures and are likely to inhibit cellular antiviral activities, including immune and RNA interference (RNAi) pathways. The seven mutations have been maintained in all Asian and American strains and may be responsible for serious medical problems we are facing today and offer testable hypotheses to examine their roles in molecular interactions during brain development.

Long-Term Virus Evolution in Nature

Viruses spread to give rise to epidemics and pandemics, and some key parameters that include virus and host population numbers determine virus persistence or extinction in nature. Viruses evolve at different rates of evolution depending on the polymerase copying fidelity during genome replication. Calculated rates of evolution in nature vary depending on the time interval between virus isolations. In particular, intra-host evolution is generally more rapid that inter-host evolution and several possible mechanisms for this difference are considered. The mechanisms by which the error-prone viruses evolve render very unlikely the operation of a molecular clock (constant rate of incorporation of mutations in the evolving genomes). Several computational methods are reviewed that permit the alignment of viral sequences and the establishment of phylogenetic relationships among viruses. The evolution of virus in the form of dynamic mutant clouds in each infected individual, together with multiple environmental influences, render the emergence and reemergence of viral pathogens an unpredictable event, another example of biological complexity.

Clonality and intracellular polyploidy in virus evolution and pathogenesis

In the present article we examine clonality in virus evolution. Most viruses retain an active recombination machinery as a potential means to initiate new levels of genetic exploration that go beyond those attainable solely by point mutations. However, despite abundant recombination that may be linked to molecular events essential for genome replication, herein we provide evidence that generation of recombinants with altered biological properties is not essential for the completion of the replication cycles of viruses, and that viral lineages (near-clades) can be defined. We distinguish mechanistically active but inconsequential recombination from evolutionarily relevant recombination, illustrated by episodes in the field and during experimental evolution. In the field, recombination has been at the origin of new viral pathogens, and has conferred fitness advantages to some viruses once the parental viruses have attained a sufficient degree of diversification by point mutations. In the laboratory, recombination mediated a salient genome segmentation of foot-and-mouth disease virus, an important animal pathogen whose genome in nature has always been characterized as unsegmented. We propose a model of continuous mutation and recombination, with punctuated, biologically relevant recombination events for the survival of viruses, both as disease agents and as promoters of cellular evolution. Thus, clonality is the standard evolutionary mode for viruses because recombination is largely inconsequential, since the decisive events for virus replication and survival are not dependent on the exchange of genetic material and formation of recombinant (mosaic) genomes.

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.

Precise identification of endogenous proviruses of NFS/N mice participating in recombination with moloney ecotropic murine leukemia virus (MuLV) to generate polytropic MuLVs

Polytropic murine leukemia viruses (MuLVs) are generated by recombination of ecotropic MuLVs with env genes of a family of endogenous proviruses in mice, resulting in viruses with an expanded host range and greater virulence. Inbred mouse strains contain numerous endogenous proviruses that are potential donors of the env gene sequences of polytropic MuLVs; however, the precise identification of those proviruses that participate in recombination has been elusive. Three different structural groups of proviruses in NFS/N mice have been described and different ecotropic MuLVs preferentially recombine with different groups of proviruses. In contrast to other ecotropic MuLVs such as Friend MuLV or Akv that recombine predominantly with a single group of proviruses, Moloney MuLV (M-MuLV) recombines with at least two distinct groups. In this study, we determined that only three endogenous proviruses, two of one group and one of another group, are major participants in recombination with M-MuLV. Furthermore, the distinction between the polytropic MuLVs generated by M-MuLV and other ecotropic MuLVs is the result of recombination with a single endogenous provirus. This provirus exhibits a frameshift mutation in the 3' region of the surface glycoprotein-encoding sequences that is excluded in recombinants with M-MuLV. The sites of recombination between the env genes of M-MuLV and endogenous proviruses were confined to a short region exhibiting maximum homology between the ecotropic and polytropic env sequences and maximum stability of predicted RNA secondary structure. These observations suggest a possible mechanism for the specificity of recombination observed for different ecotropic MuLVs.

Lethal mutagenesis of HIV

HIV-1 and other retroviruses exhibit mutation rates that are 1,000,000-fold greater than their host organisms. Error-prone viral replication may place retroviruses and other RNA viruses near the threshold of "error catastrophe" or extinction due to an intolerable load of deleterious mutations. Strategies designed to drive viruses to error catastrophe have been applied to HIV-1 and a number of RNA viruses. Here, we review the concept of extinguishing HIV infection by "lethal mutagenesis" and consider the utility of this new approach in combination with conventional antiretroviral strategies.

HIV, 'an evolving species'. Roles of cellular activation and co-infections

Each small variation of the genome of a species can be preserved if it is useful for the survival of the species in a given environment. Within this framework, the finality of the biological cycle of HIV consists in a search for harmony (biological coherence) with man, which is to say a stable condition. Cellular activation appears to be the strategy developed by HIV in order to achieve this coherence. The price of this strategy is the AIDS. The first contact between HIV and immune system appears to determine the subsequent clinical outcome and the future of HIV. Lymphocytic activation varies during the course of the vital cycle of HIV. For each individual, this lymphocytic activation depends on both the HLA repertoire acquired during thymic ontogenesis and the antigenic experience before and after HIV infection. Thus intercurrent infections alter the immune condition of the organism and influence the outcome of HIV. We described a synthetic analysis of the effects of HIV on the surface protein expression and the cellular activation pathways which should provide insights in the evolutionary relationship between HIV and man and should permit to do a more physiological therapeutic approach.

Cloning, expression and evolution of the gene encoding the elongation factor 1alpha from a low thermophilic Sulfolobus solfataricus strain

The gene encoding the elongation factor 1alpha (EF-1alpha) from the archaeon Sulfolobus solfataricus strain MT3 (optimum growth temperature 75 degrees C) was cloned, sequenced and expressed in Escherichia coli. The structural and biochemical properties of the purified enzyme were compared to those of EF-1alpha isolated from S. solfataricus strain MT4 (optimum growth temperature 87 degrees C). Only one amino acid change (Val15-->Ile) was found. Interestingly, the difference was in the first guanine nucleotide binding consensus sequence G(13)HIDHGK and was responsible for a reduced efficiency in protein synthesis, which was accompanied by an increased affinity for both guanosine diphosphate (GDP) and guanosine triphosphate (GTP), and an increased efficiency in the intrinsic GTPase activity. Despite the different thermophilicities of the two microorganisms, only very marginal effects on the thermal properties of the enzyme were observed. Molecular evolution among EF-1alpha genes from Sulfolobus species showed that the average rate of nucleotide substitution per site per year (0.0312x10(-9)) is lower than that reported for other functional genes.

Repeat-induced point mutation (RIP) in Magnaporthe grisea: implications for its sexual cycle in the natural field context

Repeat-induced point mutation (RIP) is a process that detects DNA duplications and peppers their sequences with C:G to T:A transitions in the sexual phase of the life cycle. So far, this unique mechanism has been identified as a currently active process in only two fungal species, Neurospora crassa and Podospora anserina. To determine whether a RIP-like process operates in the plant pathogenic fungus Magnaporthe grisea, the retrotransposon MAGGY and the hygromycin B phosphotransferase gene were introduced into the fungus as multiple transgenes and examined for sequence alterations after a cross. Frequent C:G to T:A transitions in the transgenes were found in the descendants, preferentially in (A/Tp)Cp(A/T)contexts, suggesting that a process similar to RIP functions in M.grisea. We also examined the sequence of another retrotransposon Pyret in six field isolates of M. grisea. Even though no perfect stage has been known in M. grisea under field conditions to date, RIP-like transitions were found in all the field isolates tested. Interestingly, the frequency of the transitions mostly correlated with the fertility of the isolates examined under laboratory conditions. These results imply that the sexual cycle of this fungus exists or existed in the natural field context.

Frequency of missense mutations in the coding region of a eukaryotic gene transferred by retroviral vectors

A relatively high mutation rate is probably a major factor in the evolutionary success of retroviruses, because it generates the genetic diversity that helps them to cope with changes in the environment. When using recombinant retroviruses as vectors for gene transfer and gene therapy, it is important to consider the implications of this biological characteristic. Until now, the mutation rate has been studied by using noneukaryotic genes as reporters. Here we report point mutations in the human glucose-6-phosphate dehydrogenase (hG6PD) gene transferred by Moloney murine leukemia virus-based vectors into murine bone marrow cells and NIH 3T3 murine fibroblasts. After bone marrow transplantation, we observed an hG6PD with abnormal electrophoretic mobility for 2 out of 34 mice. Next, we studied this phenomenon quantitatively and found 1 electrophoretically abnormal hG6PD variant among 93 independently isolated NIH 3T3 clones, from which we estimate a mutation rate of 1.4 x 10(-5) per base pair per replication cycle. Mutations in the transferred gene can thus contribute to the impairment of the effectiveness of retrovirus-mediated gene transfer.

Drug discovery and development of antiviral agents for the treatment of chronic hepatitis B virus infection

A safe and effective vaccine for hepatitis B virus (HBV) has been available for nearly twenty years and currently campaigns to provide universal vaccination in developing countries are underway. Nevertheless, chronic HBV infection remains a leading cause of chronic hepatitis worldwide and there is a strong need for safe and effective antiviral therapies. Attempts to identify and develop antiviral agents to treat chronic HBV infection remains focused on nucleoside analogs such as 3TC (lamivudine), adefovir dipivoxil, (bis-POMPMEA), and others. However, advances in our understanding of the molecular biology of HBV and the development of new assays for HBV polymerase activity, such as the reconstitution of active HBV polymerase in vitro, should facilitate large screening efforts for non-nucleoside reverse transcriptase inhibitors. Recent advances have furthered our understanding of clinical resistance to lamivudine, have provided new approaches to treatment, and have offered new perspectives on the major challenges to the identification and development of antiviral agents for chronic HBV infection. Here, in an update to our previous review article that appeared in this series [59a], we focus on recent advances that have occurred in the areas of virus structure and replication, in vitro viral polymerase assays, cell culture systems, and animal models.

Genomic stability of murine leukemia viruses containing insertions at the Env-3' untranslated region boundary

Retroviruses containing inserts of exogenous sequences frequently eliminate the inserted sequences upon spread in susceptible cells. We have constructed replication-competent murine leukemia virus (MLV) vectors containing internal ribosome entry site (IRES)-transgene cassettes at the env-3' untranslated region boundary in order to examine the effects of insert sequence and size on the loss of inserts during viral replication. A virus containing an insertion of 1.6 kb replicated with greatly attenuated kinetics relative to wild-type virus and lost the inserted sequences in a single infection cycle. In contrast, MLVs containing inserts of 1.15 to 1.30 kb replicated with kinetics only slightly attenuated compared to wild-type MLV and exhibited much greater stability, maintaining their genomic integrity over multiple serial infection cycles. Eventually, multiple species of deletion mutants were detected simultaneously in later infection cycles; once detected, these variants rapidly dominated the population and thereafter appeared to be maintained at a relative equilibrium. Sequence analysis of these variants identified preferred sites of recombination in the parental viruses, including both short direct repeats and inverted repeats. One instance of insert deletion through recombination with an endogenous retrovirus was also observed. When specific sequences involved in these recombination events were eliminated, deletion variants still arose with the same kinetics upon virus passage and by apparently similar mechanisms, although at different locations in the vectors. Our results suggest that while lengthened, insert-containing genomes can be maintained over multiple replication cycles, preferential deletions resulting in loss of the inserted sequences confer a strong selective advantage.

Mutation patterns for two flaviviruses: hepatitis C virus and GB virus C/hepatitis G virus

We studied the mutation patterns of hepatitis C virus (HCV) and GB virus C/hepatitis G virus (HGV). Although the mutation patterns of the two viruses were similar to each other, they were quite different from that of HIV. In particular, the similarity of the patterns between HCV or HGV and human nuclear pseudogenes was statistically significant whereas there was no similarity between HIV and human nuclear pseudogenes. This finding suggests that the mutation patterns of HCV and HGV are similar to the patterns of spontaneous substitution mutations of human genes, implying that nucleotide analogues which are effective against HCV and HGV may have a side effect on the normal cells of humans.

Genetic control and dynamics of the cellular immune response to the human T-cell leukaemia virus, HTLV-I

About 1% of people infected with the human T-cell leukaemia virus, type 1 (HTLV-I) develop a disabling chronic inflammatory disease of the central nervous system known as HTLV-I-associated myelopathy/tropical spastic paraparesis (HAM/TSP). Patients with HAM/TSP have a vigorous immune response to HTLV-I, and it has been widely suggested that this immune response, particularly the HTLV-I-specific cytotoxic T-lymphocyte (CTL) response, causes the tissue damage that is seen in HAM/TSP. In this paper we summarize recent evidence that a strong CTL response to HTLV-I does in fact protect against HAM/TSP by reducing the proviral load of HTLV-I. We conclude that HTLV-I is persistently replicating at a high level, despite the relative constancy of its genome sequence. These results imply that antiretroviral drugs could reduce the risk of HAM/TSP by reducing the viral load, and that an effective anti-HTLV-I vaccine should elicit a strong CTL response to the virus. The dynamic nature of the infection also has implications for the epidemiology and the evolution of HTLV-I.

Antigenic variation of varicella zoster virus Fc receptor gE: loss of a major B cell epitope in the ectodomain

Varicella zoster virus (VZV) is considered to possess a genetically stable genome; only one serotype is recognized around the world. The 125-kbp genome contains approximately 70 open reading frames. One that has received particular attention is open reading frame 68, which codes for glycoprotein gE, the predominant 623-residue viral envelope product that harbors both B and T cell epitopes. This report describes the initial characterization of a community-acquired VZV isolate that was a distinguishable second serotype (i.e., it had lost a major B cell epitope defined on the gE ectodomain by a murine monoclonal antibody called mAb 3B3). The mAb 3B3 epitope was found not only on the prototype sequenced Dumas strain from Holland and all previously tested North American isolates but also on the varicella vaccine Oka strain originally attenuated in Japan. Sequencing of the mutated gE ectodomain demonstrated that codon 150 exhibited a single base change that led to an amino acid change (aspartic acid to asparagine). Observation of the monolayers infected with the mutant VZV strain also led to the surprising discovery that the topography of egress was altered. Wild-type VZV emerges along distinctive viral highways, whereas the mutant strain virions were nearly uniformly distributed over the cell surface in a pattern more closely resembling egress of herpes simplex virus 1. The mutant VZV strain was designated VZV-MSP because it was isolated in Minnesota.

The identification and development of antiviral agents for the treatment of chronic hepatitis B virus infection

Hepatitis B virus (HBV) is the leading cause of chronic hepatitis throughout the world. Notwithstanding the availability of a safe and effective vaccine, the world prevalence of HBV has not declined significantly, thus resulting in the need for a selective antiviral agent. HBV is a small, partially double-stranded DNA virus which replicates through an RNA intermediate. Most efforts to develop anti-HBV agents have been targeted to the viral DNA polymerase which possesses reverse transcriptase activity. Currently, the most promising anti-HBV agents are nucleoside analogs which interfere with viral DNA replication. Although earlier nucleoside analogs such as vidarabine (ara-A) and fialuridine (FIAU) have displayed unacceptable toxicities, newer analogs such as lamivudine (3TC), bis-POM PMEA (GS-840), lobucavir, and BMS-200,475 have demonstrated clinical utility. In particular, the use of lamivudine has generated considerable interest in the development of other L-enantiomeric nucleoside analogs for use against HBV. Here, we provide an overview of HBV structure and replication strategy and discuss the use of cell culture systems, in vitro viral polymerase systems, and animal models to identify and evaluate anti-HBV agents. We also discuss the various classes of nucleoside analogs in terms of structure, mechanism of action, status in clinical development, ability to select for resistant HBV variants, and use in combination therapies. Finally, we present a discussion of novel antiviral approaches, including antisense and gene therapy, and address the various challenges to successful anti-HBV chemotherapeutic intervention.

The origin and evolution of human T-cell lymphotropic virus types I and II

Studies on human T-cell lymphotropic virus types I (HTLV-I) and II (HTLV-II) are briefly reviewed from the viewpoint of molecular evolution, with special reference to the evolutionary rate and evolutionary relationships among these viruses. In particular, it appears that, in contrast to the low level of variability of HTLV-I among different isolates, individual isolates form quasispecies structures. Elucidating the mechanisms connecting these two phenomena will be one of the future problems in the study of the molecular evolution of HTLV-I and HTLV-II.

Evolutionary mechanisms and population dynamics of the third variable envelope region of HIV within single hosts

Clonal diversifications of HIV virus were monitored by periodic samplings on each of the six patients with regard to 183- to 335-bp segments of the env gene, which invariably included the functionally critical V3 region. Subsequently, six individual phylogenetic trees of viral variants were constructed. It was found that at one time or another during the course of disease progression, viral variants were inexplicably released from a strong negative selection against nonsynonymous base substitutions, possibly indicating positive selection. This resulted in concentrated amino acid substitutions at five specific sites within the V3 region. It was noted that these sites were often involved as antigenic determinants that provoked the host immune response and that these sites were also involved in the determination of viral phenotypes as to their cell tropism, syncytium formation capability, and replication rates.

Molecular phylogeny of human T-cell leukemia virus type I and II of Amerindians in Colombia and Chile

Six human T-cell leukemia virus type I (HTLV-I) and eight human T-cell leukemia virus type II (HTLV-II) cases newly isolated from the South American countries of Colombia and Chile were analyzed together with the two Amerindian HTLV-I isolates previously reported. All of the HTLV-I isolates belonged to the transcontinental subgroup of the "cosmopolitan" group, and Colombian isolates, including those from native Amerindians and Negroes, formed a single tight cluster within this subgroup. The transcontinental subgroup consisted of isolates from various regions such as the Caribbean basin, India, Iran, South Africa, Sakhalin, and Japan, and included isolates from the "Ainu" and "Okinawa" people, regarded as relatively pure Japanese descended from the prehistoric "Jomon" period which began more than 10,000 years ago. This implied a dissemination of the subgroup associated with the movement of human beings in ancient times. On the other hand, all of the HTLV-II isolates from native Amerindians in Colombia and Chile belonged to the HTLV-IIb subtype which has previously been reported to be mainly endemic in certain populations of native Amerindians. The southernmost isolate from Chile, showing wide distribution of the IIb subtype in native South Amerindians and largest heterogeneity of the subtype in Colombian isolates, supported the idea that the HTLV-IIb subtype has been endemic for a long time in native Indians of South America.

Constrained evolution with respect to gene overlap of hepatitis B virus

With the aim of elucidating the evolution of a hepadnavirus family, we constructed molecular phylogenetic trees for 27 strains of hepatitis B virus (HBV) using both the unweighted pair-grouping and neighbor-joining methods. All five gene regions, P, C, S, X, and preS, were used to construct the phylogenetic trees. Using the phylogenetic trees obtained, we classified these strains into five major groups in which the strains were closely related to each other. Our classification reinforced our previous view that genetic classification is not always compatible with conventional classification determined by serological subtypes. Moreover, constraints on the evolutionary process of HBV were analyzed for amino-acid-altering (nonsynonymous) and silent (synonymous) substitutions, because two-thirds of the open reading frame (ORF), P, contains alternating overlapping ORFs. In our unique analysis of this interesting gene structure of HBV, the most frequent synonymous substitutions were observed in the nonoverlapped parts of the P and C genes. On the other hand, the number of synonymous substitutions per nucleotide site for the S gene was quite low and appeared a strongly constrained evolution. Because the P gene overlaps the S gene in a different frame, the low rate of synonymous substitution for the S gene can be explained by the evolutionary constraints which are imposed on the overlapping gene region. In other words, synonymous substitutions in the S gene can cause amino acid changes in its overlapping region in a different frame. Thus, the evolution of HBV is constrained evolutionarily by the overlapping genes. We propose calling this mode of viral evolution "constrained evolution." The evolution of HBV represents a typical constrained evolution.

Introduction of a new hepatitis agent in retrospect: genetic studies of Swedish hepatitis D virus strains

The supposedly first outbreak of hepatitis D virus (HDV) infection in Sweden occurred among intravenous drug addicts in the Malmö area in the mid-1970s. Stored sera from this outbreak were used for viral RNA extraction and analysis. By sequence comparisons, the HDV genomes from those Swedish patients fell into two separate clusters, within which the RNA sequences were closely related. These two HDV groups genetically resembled the French and US-1 isolates of genotype I, respectively, indicating that there had been at least two separate sources of HDV infection. The genetic alterations of the HDV RNA were investigated by sequence analysis of nine annually drawn serum samples from one patient and paired samples collected between 2 and more than 10 years apart from six patients with chronic HDV infection. Only mutational changes were observed, and no insertion or deletion appeared throughout the periods observed. It was found that the Swedish HDV isolates mutationally evolved at an average rate of 1.1 x 10(-3) substitutions per nucleotide per year over a long time course of chronic HDV infection, which is of the same magnitude as that of other RNA viruses.

Quasispecies and naturally occurring superinfection in feline immunodeficiency virus infection

Analysis of individual clones containing the V1 and V2 domains of the segment of the FIV env gene present in a naturally infected cat (T) was carried out. The polymerase chain reaction (PCR) was used to amplify proviral FIV DNA extracted from peripheral blood mononuclear cells (PBMCs) obtained in October 1994 from this cat. The PCR products were cloned and the DNA sequences determined for 11 clones. Sequences obtained were aligned with sequences corresponding to FIV isolates (T90, T91, T92) previously obtained from the same cat in 1990, 1991 and 1992. Phylogenetic analysis was performed which included consensus sequences of another Australian isolate, N91, as well as UK, US, Swiss and Japanese isolates of FIV. All clones varied from each other, and none of these clones was identical to the consensus sequences of the isolates obtained previously from the same cat (the T-series). However, most of these clones appeared to have originated from the ancestor of the most recent isolate (T92). In addition, 2 of the clones (7&11) are closely related to another Australian isolate N91, obtained from a different cat (N) in 1991. Because these two cats (T and N) were housed together for at least 3 years (1990-1993) it is suggested that the first cat (T) has become superinfected with an isolate from a second cat (N) under natural conditions. The identification of clones of differing sequences, which were not identical to each other nor to their ancestors, emphasises the rapid mutation of lentiviruses within the env region, and the difficulty of developing an effective FIV vaccine. More importantly, the possibility of natural superinfection with FIV in cats has implications for the development of a successful lentiviral vaccine.

Evolution of multiple families of non-LTR retrotransposons in phlebotomine sandflies

In this paper we report on the diversity and distribution of a set of non-LTR retrotransposon (RTP) reverse transcriptase (RT) sequences isolated from phlebotomine sandflies, and their potential for investigating the evolutionary histories of members of this subfamily of flies (Diptera:Psychodidae, Phlebotominae). The phlebotomine RT sequence families derived from one species were as different from each other as they were from RT sequences derived from other species. When each was used to probe Southern blots of sandfly genomic DNA they hybridized only to the species of source and, usually, to others of the same subgenus, but not to DNA from other subgenera-a hybridization pattern consistent with vertical evolution. There was considerable intraspecific variation in hybridization pattern, suggesting the RTs were part of non-LTR RTPs that are (or were recently) subject to flux in genomic position and copy number. Most of the RT families detected in phlebotomines are monophyletic with respect to previously described RTs, and all are monophyletic with RTs of the F/Jockey (Drosophila melanogaster) type of RTP. Orthologous sequences were isolated from the closely related species Phlebotomus perniciosus and P. tobbi (subgenus Larroussius), and different populations of P. perniciosus. The level of sequence divergence among these orthologous RTs, the subgeneric distribution of each RT family, and the intraspecific variation in hybridization pattern of many of them, indicate this class of sequence will provide genetic markers at the sub-generic level.

Substitution rates in hepatitis delta virus

Substitution rates were estimated for the coding and noncoding regions of the hepatitis delta virus (HDV). The estimated rates of synonymous substitution in HDV were lower than the rates of substitution at non-synonymous sites and in the noncoding region. HDV has lower synonymous substitution rates than the hepatitis C virus, though both are RNA viruses. The relatively low rate of synonymous substitution in HDV may be due to a strong preference of G and C nucleotides at third codon positions. Variation in substitution rate among HDV lineages may be correlated with the clinical development of the HDV-induced hepatitis. The phylogenetic tree inferred for 24 HDV strains reveals similarities between lineages isolated from the same geographic region.

Analysis and chromosomal localization of retrotransposons in sugar beet (Beta vulgaris L.): LINEs and Ty1-copia-like elements as major components of the genome

DNA sequences of the reverse transcriptase gene of long terminal repeat (LTR) and non-LTR (non-viral) retrotransposons have been isolated and cloned from the genome of sugar beet (Beta vulgaris). Both retrotransposon types are highly amplified in sugar beet and may account for 2-5% of the genome. The BNR1 family, representing the first non-viral retrotransposon reported from a dicotyledonous species, shows homology to the mammalian L1 family of long interspersed repeated sequences (LINEs) and to retrotransposable elements from maize and lily. Sequences of the Tbv family are homologous to the Ty1-copia class of LTR retrotransposons. The BNR1 and Tbv retrotransposon families are characterized by sequence heterogeneity and are probably defective. The deduced peptide sequences were used to investigate the relation to other retroelements from plants, insects and mammals. Fluorescence in situ hybridization was used to investigate the physical distribution and revealed that both retrotransposon families are present on all sugar beet chromosomes and largely excluded from chromosomal regions harbouring the 18S-5.8S-25S rRNA genes. The BNR1 family is organized in discrete clusters, while the Tbv family of Ty1-copia-like retrotransposons shows a more uniform distribution along chromosome arms and is absent from some chromosomal regions. These contrasting distributions emphasize the differences in evolutionary amplification and dispersion mechanisms between the two types of retrotransposons. The in situ results of both elements reflect significant features of a higher order structure of the genome, as it is known for both short interspersed repeated sequences (SINEs) and LINEs in human.

Long terminal repeat enhancer core sequences in proviruses adjacent to c-myc in T-cell lymphomas induced by a murine retrovirus

The transcriptional enhancer in the long terminal repeat (LTR) of the T-lymphomagenic retrovirus SL3-3 differs from that of the nonleukemogenic virus Akv at several sites, including a single base pair difference in an element termed the enhancer core. Mutation of this T-A base pair to the C-G C-G sequence found in Akv significantly attenuated the leukemogenicity of SL3-3. Thus, this difference is important for viral leukemogenicity. Since Akv is an endogenous virus, this suggests that the C-G in its core is an adaptation to being minimally pathogenic. Most tumors that occurred in mice inoculated with the mutant virus, called SAA, contained proviruses with reversion or potential suppressor mutations in the enhancer core. We also found that the 72-bp tandem repeats constituting the viral enhancer could vary in number. Most tumors contained mixtures of proviruses with various numbers of 72-bp units, usually between one and four. Variation in repeat number was most likely due to recombination events involving template misalignment during viral replication. Thus, two processes during viral replication, misincorporation and recombination, combined to alter LTR enhancer structure and generate more pathogenic variants from the mutant virus. In SAA-induced tumors, enhancers of proviruses adjacent to c-myc had the largest number of core reversion or suppressor mutations of all of the viral enhancers in those tumors. This observation was consistent with the hypothesis that one function of the LTR enhancers in leukemogenesis is to activate proto-oncogenes such as c-myc.

Evolution of pathogenic viruses with special reference to the rates of synonymous and nonsynonymous substitutions

For pathogenic viruses such as human immunodeficiency virus (HIV), hepatitis B virus (HBV), hepatitis C virus (HCV), human influenza A virus, and human T-cell leukemia virus type I (HTLV-I), the evolutionary features were briefly reviewed with special reference to the rates of synonymous and nonsynonymous substitutions. In particular, these rates were discussed in connection with the neutral theory of molecular evolution. It was common to all the five pathogenic viruses that the rate of synonymous substitution was higher than that of nonsynonymous substitution particularly when the entire gene regions were compared between different isolates. This suggests that the viral proteins are quite conservative to functional and structural changes even though most of these viral genomes are evolving at a speed extraordinarily higher than their host genomes. Thus, this feature is consistent with the neutral theory. However, it is also pointed out that positive selection may be operating on some specific sites such as antigenic sites in order for the pathogenic viruses to escape from the host immune system.

Molecular evolution: codes, clocks, genes and genomes

The discoveries, advancements and continuing controversies in the field of molecular evolution are reviewed. Topics summarized are (1) the evolution of the genetic code, (2) gene evolution including the demonstration of homology, estimation of sequence divergence, phylogenetic trees, the molecular clock and the origin of genes and gene families by various genetic mechanisms, and (3) eukaryotic genome evolution, including the highly repeated satellite sequences, the interspersed and potentially mobile repeated sequences and the unique sequence fraction of the genome.

Subclonal components of consensus fitness in an RNA virus clone

Most RNA virus populations exhibit extremely high mutation frequencies which generate complex, genetically heterogeneous populations referred to as quasi-species. Previous work has shown that when a large spectrum of the quasi-species is transferred, natural selection operates, leading to elimination of noncompetitive (inferior) genomes and rapid gains in fitness. However, whenever the population is repeatedly reduced to a single virion, variable declines in fitness occur as predicted by the Muller's ratchet hypothesis. Here, we quantitated the fitness of 98 subclones isolated from an RNA virus clonal population. We found a normal distribution around a lower fitness, with the average subclone being less fit than the parental clonal population. This finding demonstrates the phenotypic diversity in RNA virus populations and shows that, as expected, a large fraction of mutations generated during virus replication is deleterious. This clarifies the operation of Muller's ratchet and illustrates why a large number of virions must be transferred for rapid fitness gains to occur. We also found that repeated genetic bottleneck passages can cause irregular stochastic declines in fitness, emphasizing again the phenotypic heterogeneity present in RNA virus populations. Finally, we found that following only 60 h of selection (15 passages in which virus yields were harvested after 4 h), RNA virus populations can undergo a 250% average increase in fitness, even on a host cell type to which they were already well adapted. This is a remarkable ability; in population biology, even a much lower fitness gain (e.g., 1 to 2%) can represent a highly significant reproductive advantage. We discuss the biological implications of these findings for the natural transmission and pathogenesis of RNA viruses.

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 x 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 x 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.

A model for the mechanism of initial generation of short interspersed elements (SINEs)

Most animal genomes contain a large number of short interspersed elements (SINEs) that have a composite structure and contain a region that is homologous to a tRNA. The majority of SINEs have been found to be derived from a tRNA(Lys), being categorized as members of a superfamily of tRNA(Lys)-related SINEs. The consensus sequences of five SINEs that belong to this superfamily were aligned. It was found that, in the tRNA-unrelated region, there are two sequence motifs that are almost identical among these five SINEs and are at a distance of 10-11 nucleotides from each other. This observation suggests a common evolutionary origin of these SINEs and/or some function(s) for these motifs. Similar sequences were unexpectedly found to be present in the sequences complementary to the U5 regions of several mammalian retroviruses whose primer is a tRNA(Lys). On the basis of these findings, we propose a possible model for the generation of SINEs whereby they are derived from a "strong stop DNA" with a primer tRNA that is an intermediate in the process of reverse transcription of certain retroviruses.

BARE-1, a copia-like retroelement in barley (Hordeum vulgare L.)

Retroviruses and retrotransposons make up the broad class of retroelements replicating and transposing via reverse transcriptase. Retroelements have recently been found to be ubiquitous in the plants. We report here the isolation, sequence and analysis of a retroelement from barley (Hordeum vulgare L.) with all the features of a copia-like retrotransposon. This is named BARE-1 (for BArley RetroElement 1), the first such element described for barley. BARE-1 is 12,088 bp, with long terminal repeats (LTRs) of 1829 bp containing perfect 6 bp inverted repeats at their ends and flanked by 4 bp direct repeats in the host DNA. Between the long terminal repeats is an internal domain with a derived amino acid sequence of 1285 residues, bearing homology to the gag, pro, int and rt domains of retroviruses and both plant and non-plant copia-like retrotransposons. Cultivated barley contains about 5000 elements in the genome similar to the BARE-1 putative gag domain, but ten-fold more hybridizing to rt or LTR probes. The particular BARE-1 element reported here appears to be inactive, as the putative protein-coding domain is interrupted by four stop codons and a frameshift. In addition, the 3' LTR is 4% divergent from the 5' LTR and contains a 3135 bp insertion. Nevertheless, we have recently detected transcripts hybridizing to BARE-1 on northern blots, presumably from active copies. Analysis of BARE-1 expression and function in barley is currently underway.

Several short interspersed repetitive elements (SINEs) in distant species may have originated from a common ancestral retrovirus: characterization of a squid SINE and a possible mechanism for generation of tRNA-derived retroposons

Using labeled transcripts generated in vitro from squid total genomic DNA as a probe, we isolated and characterized a SINE that is present in the squid genome. The squid SINE appears to be derived from a tRNA(Lys). When the consensus sequences of five different SINEs with a tRNA(Lys)-like structure from distantly related species, including squid, were aligned, we found in the tRNA-unrelated region two sequence motifs that were almost identical among these five SINEs. This observation suggests a common evolutionary origin for these SINEs and/or some function(s) for these motifs. Similar sequences were unexpectedly found to be present in sequences complementary to the U5 regions of several mammalian retroviruses whose primer is a tRNA(Lys). On the basis of these findings, we present a model for the generation of SINEs. We propose that they are derived from a "strong-stop DNA" with a primer tRNA(Lys) that is an intermediate in the reverse transcription of certain retroviruses. Our model suggests that a certain group of SINEs may have been generated by horizontal transmission, although it is not clear whether information was transmitted via a similar retrovirus or via an RNA or DNA of a SINE.

Nucleotide sequence and molecular evolution of the gene coding for glyceraldehyde-3-phosphate dehydrogenase in the thermoacidophilic archaebacterium Sulfolobus solfataricus

A Sulfolobus solfataricus genomic library cloned in the EMBL3 phage was screened using as probes synthetic oligonucleotides designed from the known amino acid sequence of a peptide obtained from the purified glyceraldehyde-3-phosphate dehydrogenase (aGAPD) protein. The screening led to the isolation of six recombinant phages (lambda G1-lambda G6) and one of them (lambda G4) contained the entire GAPD gene. The deduced amino acid sequence accounts for a protein made of 341 amino acids and the initial methionine is encoded by a GTG triplet. Alignment of the S. solfataricus aGAPD sequence versus GAPD from archaea, eukarya, and bacteria showed that aGAPD is very similar to other archaebacterial but not to eukaryotic or eubacterial GAPD. For known archaebacterial GAPD sequences, the rate of nucleotide substitutions per site per year showed that these sequences are homologous not only at the amino acid but also at the nucleotide level. The evolutionary rates are nearly similar to those reported for other eukaryotic genes.

Extensive sequence variation of feline immunodeficiency virus env genes in isolates from naturally infected cats

In an investigation of the evolution of feline immunodeficiency virus (FIV) in vivo, sequential isolates from a persistently infected cat were examined by direct sequencing following amplification of selected subgenomic regions by polymerase chain reaction (PCR). Three isolates, T90, T91, and T92, obtained over a three-year period revealed no changes to regions known to be conserved within gag and pol genes. Additionally, no change occurred within gag and pol in an isolate recovered from a second cat which was experimentally infected with T90. Changes were detected within an N-terminal region of the envelope glycoprotein gp 120 (env). These consisted of point mutations, some of which would result in amino acid substitutions and the predicted amino acid changes tended to cluster within variable domains. Inoculation of T90 into a second cat resulted in a different pattern of mutations than that observed for the three isolates from the first cat. In all cases, virus isolates derived from the same cat were much more highly related to each other (extent of env variation was 0.5-1.5%) than to isolates from other cats (10-12% env variation). The rate of change of FIV was estimated to be 3.4 x 10(-3) nucleotide substitutions per site per year for the env gene and less than 10(-4) nucleotide substitutions per site per year for the gag and pol genes, values concordant with that found for human immunodeficiency virus 1. Both nucleotide and amino acid changes in the gp 120 region were found to be directional, suggesting that selective pressures influence FIV envelope gene sequences.

Six strains of human immunodeficiency virus type 1 isolated in Japan and their molecular phylogeny

Five strains of human immunodeficiency virus type 1 (HIV-1) were isolated from five Japanese hemophilia patients. Two isolates, HIV-1[GUN-1] and HIV-1[GUN-2], were from brother patients with hemophilia B and the other three isolates, HIV-1[GUN-3], HIV-1[GUN-4], and HIV-1[GUN-5], were from hemophilia A patients. Another HIV-1 strain, HIV-1[GUN-6], was isolated from a Canadian male homosexual with AIDS. The restriction endonuclease cleavage maps of the proviral genomes of these six HIV-1 strains revealed that they were apparently different from each other. The phylogenetic trees constructed using restriction maps and nucleotide sequences were quite similar, indicating that phylogenetic analyses of Japanese HIV-1 isolates can be done using restriction maps of the proviruses. Phylogenetic analyses showed that they were more closely related to HIV-1s which had been reported to be isolated from homosexual patients in the United States than those isolated from African patients. In particular, GUN-1 and GUN-2 isolates were on the branch of a San Francisco isolate, ARV2, while GUN-5 and GUN-6 isolates were on the branch of HTLV-IIIB-related isolates.

Genetic mechanisms of antigenic variation in infectious bursal disease virus: analysis of a naturally occurring variant virus

The major immunogenic protein VP2 from a pathogenic field isolate (variant A virus) of infectious bursal disease virus (IBDV) was cloned and sequenced to examine antigenic variations. The VP2 open reading frame consists of 1509 nucleotides and codes for a 503 amino acid protein. Overall, the VP2 amino acid sequence of the variant A virus shares 98.6% identity with VP2 genes from other published IBDV strains. However, within the central region of VP2 (amino acids 222-334) lies a highly divergent area that we have termed the variable domain. Relative to five other IBDV isolates, a total of six amino acid changes occur within the variable domain of the variant A virus. At positions 284-288, a substitution of isoleucine to threonine, a decrease in the number of Chou and Fasman beta turns, and a switch from a hydrophilic to a hydrophobic region are found only in the variant A virus. Together these changes predict a decrease in antigenicity as determined by calculation of potential antigenic sites. This suggests that only minor changes within VP2 contributed to the emergence of a variant virus that can cause disease in immunized birds.

Determination of the poliovirus RNA polymerase error frequency at eight sites in the viral genome

The poliovirus RNA polymerase error frequency was measured in vivo at eight sites in the poliovirus genome. The frequency at which specific G residues in poliovirion RNA changed to another base during one round of viral RNA replication was determined. Poliovirion RNA uniformly labeled with 32Pi was hybridized to a synthetic DNA oligonucleotide that was complementary to a sequence in the viral genome that contained a single internal G residue. The nonhybridized viral RNA was digested with RNase T1, and the protected RNA oligonucleotide was purified by gel electrophoresis. The base substitution frequency at the internal G residue was measured by finding the fraction of this RNA oligonucleotide that was resistant to RNase T1 digestion. A mean value of 2.0 x 10(-3) +/- 1.2 x 10(-3) was obtained at two sites. A modification of the above procedure involved the use of 5'-end-labeled RNA oligonucleotides. The mean value of the error frequency determined at eight sites in the viral genome by using this technique was 4.1 x 10(-3) +/- 0.6 x 10(-3). Sequencing two of the RNase T1-resistant RNA oligonucleotides confirmed that the internal G was changed to a C, A, or U residue in most of these oligonucleotides. Thus, our results indicated that the polymerase had a high error frequency in vivo and that there was no significant variation in the values determined at the specific sites examined in this study.

Direct determination of the point mutation rate of a murine retrovirus

The point mutation rate of a murine leukemia virus (MuLV) genome (AKV) was determined under conditions in which the number of replicative cycles was carefully controlled and the point mutation rate was determined by direct examination of the RNA genomes of progeny viruses. A clonal cell line infected at a low multiplicity of infection (2 x 10(-3)) was derived to provide a source of virus with high genetic homogeneity. Virus stocks from this cell line were used to infect cells at a low multiplicity of infection, and the cells were seeded soon after infection to obtain secondary clonal cell lines. RNase T1-oligonucleotide fingerprinting analyses of virion RNAs from 93 secondary lines revealed only 3 base changes in nearly 130,000 bases analyzed. To obtain an independent assessment of the mutation rate, we directly sequenced virion RNAs by using a series of DNA oligonucleotide primers distributed across the genome. RNA sequencing detected no mutations in over 21,000 bases analyzed. The combined fingerprinting and sequencing analyses yielded a mutation rate for infectious progeny viruses of one base change per 50,000 (2 x 10(-5)) bases per replication cycle. Our results suggest that over 80% of infectious progeny MuLVs may be replicated with complete fidelity and that only a low percentage undergo more than one point mutation during a replication cycle. Previous estimates of retroviral mutation rates suggest that the majority of infectious progeny viruses have undergone one or more point mutations. Recent studies of the mutation rates of marker genes in spleen necrosis virus-based vectors estimate a base substitution rate lower than estimates for infectious avian retroviruses and nearly identical to our determinations with AKV. The differences between mutation rates observed in studies of retroviruses may reflect the imposition of different selective conditions.

Ty1-copia group retrotransposons and the evolution of retroelements in the eukaryotes

Ty1-copia group retrotransposons are among the best studied transposable elements in the eukaryotes. This review discusses the extent of these transposons in the eukaryote kingdoms and compares models for the evolution of these genetic elements in the light of recent phylogenetic data. These data show that the Ty1-copia group is widespread among invertebrate eukaryotes, especially in the higher plant kingdom, where these genetic elements are unusually common and heterogeneous in their sequence. The phylogenetic data also suggest that the present day spectrum of Ty1-copia group retrotransposons has been influenced both by divergence during vertical transmission down evolving lineages and by horizontal transmission between distantly related species. Lastly, the factors affecting Ty1-copia group retrotransposon copy number and sequence heterogeneity in eukaryotic genomes and the effects of transpositional quiescence and defective retrotransposons upon evolution of Ty1-copia group retrotransposons are discussed.

Nucleotide sequence analysis of isolates of human T-lymphotropic virus type 1 of diverse geographical origins

Nucleotide sequences for long terminal repeat (LTR), gag, the protease gene, and pol of a human T-lymphotropic virus type 1 (HTLV-1) isolate of probable Caribbean origin (HTLV-1CH) and a Zairian isolate (HTLV-1EL) were determined providing complete proviral sequences for these isolates. These sequences were compared with those available from previously analyzed isolates. Nucleotide sequence differences of 1.2-3.3% were identified among isolates for which complete genetic information was available. Nucleotide sequence diversity was distributed relatively evenly over the genome with 1.3-5.2% differences in the LTR, 1.1-2.9% differences in gag, 0.7-2.1% differences in the protease gene, 0.9-2.5% differences in pol, 0.9-2.4% differences in env, 0.0-1.4% differences in rex, and 0.1-2.6% differences in tax. There were 1.2-2.3% amino acid differences overall, with 0.8-1.6% nonconservative amino acid alterations. Nucleotide differences were not found in regions of the LTR which are important for transcriptional activity or Tax response. Within the Rex-response element, nucleotide differences were found predominantly in loop rather than stem structures, thus, maintaining the overall secondary structure necessary for Rex activity. Evolutionary tree analysis of the sequence differences suggests a predominant clustering of different HTLV1 strains according to geographical origin. An open reading frame was also identified on the minus DNA strand situated between the env and rex/tax genes which exhibits 0.1-6.9% nucleotide sequence variation among HTLV1 strains. The limited sequence variation among HTLV-1 strains is in striking contrast to the extensive heterogeneity seen among human immunodeficiency virus (HIV) strains.