Contribution of the gag and pol sequences to the leukemogenicity of Friend murine leukemia virus.

Friend murine leukemia virus (F-MuLV) is a highly leukemogenic replication-competent murine retrovirus. Both the F-MuLV envelope gene and the long terminal repeat (LTR) contribute to its pathogenic phenotype (A. Oliff, K. Signorelli, and L. Collins, J. Virol. 51:788-794, 1984). To determine whether the F-MuLV gag and pol genes also possess sequences that affect leukemogenicity, we generated recombinant viruses between the F-MuLV gag and pol genes and two other murine retroviruses, amphotrophic clone 4070 (Ampho) and Friend mink cell focus-inducing virus (Fr-MCF). The F-MuLV gag and pol genes were molecularly cloned on a 5.8-kilobase-pair DNA fragment. This 5.8-kilobase-pair F-MuLV DNA was joined to the Ampho envelope gene and LTR creating a hybrid viral DNA, F/A E+L. A second hybrid viral DNA, F/Fr ENV, was made by joining the 5.8-kilobase-pair F-MuLV DNA to the Fr-MCF envelope gene plus the F-MuLV LTR. F/A E+L and F/Fr ENV DNAs generated recombinant viruses upon transfection into NIH 3T3 cells. F/A E+L virus (F-MuLV gag and pol, Ampho env and LTR) induced leukemia in 20% of NIH Swiss mice after 6 months. Ampho-infected mice did not develop leukemia. F/Fr ENV virus (F-MuLV gag and pol, Fr-MCV env, F-MuLV LTR) induced leukemia in 46% of mice after 3 months. Recombinant viruses containing the Ampho gag and pol, Fr-MCF env, and F-MuLV LTR caused leukemia in 38% of mice after 6 months. We conclude that the F-MuLV gag and pol genes contain sequences that contribute to the pathogenicity of murine retroviruses. These sequences can convert a nonpathogenic virus into a leukemia-causing virus or increase the pathogenicity of viruses that are already leukemogenic.     

Evolutionary insights into the ecology of coronaviruses

Although many novel members of the Coronaviridae have recently been recognized in different species, the ecology of coronaviruses has not been established. Our study indicates that bats harbor a much wider diversity of coronaviruses than any other animal species. Dating of different coronavirus lineages suggests that bat coronaviruses are older than those recognized in other animals and that the human severe acute respiratory syndrome (SARS) coronavirus was directly derived from viruses from wild animals in wet markets of southern China. Furthermore, the most closely related bat and SARS coronaviruses diverged in 1986, an estimated divergence time of 17 years prior to the outbreak, suggesting that there may have been transmission via an unknown intermediate host. Analysis of lineage-specific selection pressure also indicated that only SARS coronaviruses in civets and humans were under significant positive selection, also demonstrating a recent interspecies transmission. Analysis of population dynamics revealed that coronavirus populations in bats have constant population growth, while viruses from all other hosts show epidemic-like increases in population. These results indicate that diverse coronaviruses are endemic in different bat species, with repeated introductions to other animals and occasional establishment in other species. Our findings suggest that bats are likely the natural hosts for all presently known coronavirus lineages and that all coronaviruses recognized in other species were derived from viruses residing in bats. Further surveillance of bat and other animal populations is needed to fully describe the ecology and evolution of this virus family.     

Novel Endogenous Type C Retrovirus in Baboons: Complete Sequence, Providing Evidence for Baboon Endogenous Virus gag-pol Ancestry

A complete endogenous type C viral genome has been isolated from a baboon genomic library. The provirus, Papio cynocephalus endogenous retrovirus (PcEV), is 8,572 nucleotides long, and 38 to 59 proviral copies per baboon genome are found. The PcEV provirus possesses the typical simple retroviral gene organization, including two long terminal repeats and genes encoding gag, pol, and env proteins. The open reading frames for gag-pol and env are complete but have premature stop codons or frameshift mutations. The primer binding site of PcEV is complementary to tRNAGly. The gag and pol genes of PcEV are closely related to those of the baboon endogenous virus (BaEV). The env coding region of PcEV is related to the env genes of type C retroviruses. This suggests that PcEV is one of the ancestors of BaEV contributing the type C gag-pol genome fragment to the type C/D recombinant virus BaEV. Earlier it was shown that another endogenous type D virus (simian endogenous retrovirus) provided the env gene for BaEV (A. C. van der Kuyl et al., J. Virol. 71:3666-3676, 1997).     

Molecular fossils illuminate the evolution of retroviruses following a macroevolutionary transition from land to water

The ancestor of cetaceans underwent a macroevolutionary transition from land to water early in the Eocene Period >50 million years ago. However, little is known about how diverse retroviruses evolved during this shift from terrestrial to aquatic environments. Did retroviruses transition into water accompanying their hosts? Did retroviruses infect cetaceans through cross-species transmission after cetaceans invaded the aquatic environments? Endogenous retroviruses (ERVs) provide important molecular fossils for tracing the evolution of retroviruses during this macroevolutionary transition. Here, we use a phylogenomic approach to study the origin and evolution of ERVs in cetaceans. We identify a total of 8,724 ERVs within the genomes of 25 cetaceans, and phylogenetic analyses suggest these ERVs cluster into 315 independent lineages, each of which represents one or more independent endogenization events. We find that cetacean ERVs originated through two possible routes. 298 ERV lineages may derive from retrovirus endogenization that occurred before or during the transition from land to water of cetaceans, and most of these cetacean ERVs were reaching evolutionary dead-ends. 17 ERV lineages are likely to arise from independent retrovirus endogenization events that occurred after the split of mysticetes and odontocetes, indicating that diverse retroviruses infected cetaceans through cross-species transmission from non-cetacean mammals after the transition to aquatic life of cetaceans. Both integration time and synteny analyses support the recent or ongoing activity of multiple retroviral lineages in cetaceans, some of which proliferated into hundreds of copies within the host genomes. Although ERVs only recorded a proportion of past retroviral infections, our findings illuminate the complex evolution of retroviruses during one of the most marked macroevolutionary transitions in vertebrate history.     

Molecular epidemiology of simian immunodeficiency virus SIVsm in U.S. primate centers unravels the origin of SIVmac and SIVstm

Retrospective molecular epidemiology was performed on samples from four sooty mangabey (SM) colonies in the United States to characterize simian immunodeficiency virus SIVsm diversity in SMs and to trace virus circulation among different primate centers (PCs) over the past 30 years. The following SIVsm sequences were collected from different monkeys: 55 SIVsm isolates from the Tulane PC sampled between 1984 and 2004, 10 SIVsm isolates from the Yerkes PC sampled in 2002, 7 SIVsm isolates from the New Iberia PC sampled between 1979 and 1986, and 8 SIVsm isolates from the California PC sampled between 1975 and 1977. PCR and sequencing were done to characterize the gag, pol, and env gp36 genes. Phylogenetic analyses were correlated with the epidemiological data. Our analysis identified nine different divergent phylogenetic lineages that cocirculated in these four SM colonies in the Unites States in the past 30 years. Lineages 1 to 5 have been identified previously. Two of the newly identified SIVsm lineages found in SMs are ancestral to SIVmac251/SIVmac239/SIVmne and SIVstm. We further identified the origin of these two macaque viruses in SMs from the California National Primate Research Center. The diversity of SIVsm isolates in PCs in the United States mirrors that of human immunodeficiency virus type 1 (HIV-1) group M subtypes and offers a model for the molecular epidemiology of HIV and a new approach to vaccine testing. The cocirculation of divergent SIVsm strains in PCs resulted in founder effects, superinfections, and recombinations. This large array of SIVsm strains showing the same magnitude of diversity as HIV-1 group M subtypes should be extremely useful for modeling the efficacy of vaccination strategies under the real-world conditions of HIV-1 diversity. The genetic variability of SIVsm strains among PCs may influence the diagnosis and monitoring of SIVsm infection and, consequently, may bias the results of pathogenesis studies.     

Within-population variation in prevalence and lineage distribution of avian malaria in blue tits, Cyanistes caeruleus

The development of molecular genetic screening techniques for avian blood parasites has revealed many novel aspects of their ecology, including greatly elevated diversity and complex host-parasite relationships. Many previous studies of malaria in birds have treated single study populations as spatially homogeneous with respect to the likelihood of transmission of malaria to hosts, and we have very little idea whether any spatial heterogeneity influences different malaria lineages similarly. Here, we report an analysis of variation in the prevalence and cytochrome b lineage distribution of avian malaria infection with respect to environmental and host factors, and their interactions, in a single blue tit (Cyanistes caeruleus) population. Of 11 Plasmodium and Haemoproteus cytochrome b lineages found in 997 breeding individuals, the three most numerous (pSGS1, pTURDUS1 and pBT7) were considered separately, in addition to analyses of all avian malaria lineages pooled. Our analyses revealed marked spatial differences in the prevalence and distribution of these lineages, with local prevalence of malaria within the population ranging from over 60% to less than 10%. In addition, we found several more complex patterns of prevalence with respect to local landscape features, host state, parasite genotype, and their interactions. We discuss the implications of such heterogeneity in parasite infection at a local scale for the study of the ecology and evolution of infectious diseases in natural populations. The increased resolution afforded by the combination of molecular genetic and geographical information systems (GIS) tools has the potential to provide many insights into the epidemiology, evolution and ecology of these parasites in the future.     

Determinants of high titer in recombinant porcine endogenous retroviruses

Porcine endogenous retroviruses (PERVs) pose a potential stumbling block for therapeutic xenotransplantation, with the greatest threat coming from viruses generated by recombination between members of the PERV subgroup A (PERV-A) and PERV-C families (PERV-A/C recombinants). PERV-A and PERV-B have been shown to infect human cells in culture, albeit with low titers. PERV-C has a more restricted host range and cannot infect human cells. A recombinant PERV-A/C virus (PERV-A14/220) contains the PERV-A sequence between the end of pol and the middle of the SU region in env. The remaining sequence is derived from PERV-C. PERV-A14/220 is approximately 500-fold more infectious than PERV-A. To determine the molecular basis for the increased infectivity of PERV-A14/220, we have made a series of vector constructs. The primary determinant for the enhanced replicative potential of the recombinant virus appeared to be the env gene. Using a series of chimeric env genes, we could identify two determinants of high infectivity; one was an isoleucine to valine substitution at position 140 between variable regions A and B, and the other lies within the proline rich region. Taken together, these results show that the novel juxtaposition of env gene sequences enhanced the infectivity of PERV-A14/220 for human cells, perhaps by stabilization of the envelope glycoprotein or increased receptor binding.     

Heteroduplex analysis of molecular clones of the pathogenic Friend virus complex: Friend murine leukemia virus, Friend mink cell focus-forming virus, and the polycythemia- and anemia-inducing strains of Friend spleen focus-forming virus.

The pathogenic Friend virus complex is of considerable interest in that, although members of this group are genetically related, they differ markedly in biochemical and biological properties. Heteroduplex mapping of molecular clones of the Friend virus complex, which includes the replication-competent ecotropic Friend murine leukemia virus (F-MuLV) and mink cell focus-forming virus (F-MCF) and replication-defective polycythemia- and anemia-inducing strains of spleen focus-forming virus (SFFVp and SFFVa, respectively), was employed to provide insight into the molecular basis of their relationships. In heteroduplexes of F-MuLV X F-MCF, a major substitution of 0.89 kilobases in the env gene of F-MCF was discerned. Heteroduplexes of SFFVp X F-MuLV or F-MCF and SFFVa X F-MuLV or F-MCF showed several major deletions in the pol gene region and a single major deletion in the 3' half of the env gene region of SFFVp and SFFVa. A major substitution of 0.89 kilobases was mapped to the 5' end of the env deletion of SFFVp and SFFVa in heteroduplexes with F-MuLV, similar to that seen in F-MuLV X F-MCF heteroduplexes. In contrast, this env gene region was totally homologous in F-MCF X SFFVp or SFFVa and SFFVp X SFFVa heteroduplexes. Our results suggest that (i) both SFFVp and SFFVa lack part of the env gene at its 3' end, corresponding to the p15(E) coding region, (ii) major deletions occur in the pol and env genes which account for the replication defectiveness of SFFVp and SFFVa, (iii) minor substitutions occur in the gag gene region of SFFVa that are not present in SFFVp, F-MuLV, or F-MCF, (iv) a major substitution exists in the gp70 region of the env gene between F-MuLV and F-MCF that probably accounts for the differences in their host range specificities, (v) this substitution in F-MCF is identical to the gp70 part of the gp52 coding region of SFFVp and SFFVa, and (vi) heteroduplexes to F-MCF show unambiguously that no additional large substitutions are present in SFFVp or SFFVa that could account for differences in their leukemogenicity.     

Nucleotide sequence of Mason-Pfizer monkey virus: an immunosuppressive D-type retrovirus

The genetic structure of Mason-Pfizer monkey virus (MPMV), a D-type retrovirus, has been determined. In addition to the viral gag, pol, and env genes is an ORF overlapping both gag and pol and that encodes the viral protease. Surprisingly, the MPMV env protein is highly homologous to that of the avian C-type virus, reticuloendotheliosis associated virus REV-A. The env sequence encodes an immunosuppressive peptide, which suggests that MPMV, like REV-A, may transiently induce a T-suppressor cell population. The different phylogenies of the MPMV pol and env genes indicate a recombinatorial origin for the D-type viruses. Sequence comparisons show that SRV-1, an MPMV-like virus etiologically linked to simian AIDS (SAIDS), is in fact a variant of MPMV. While MPMV-like viruses cannot be used as direct models for the AIDS/SAIDS associated with lentiviruses, they provide an important system for studying the molecular basis of immunosuppressive diseases in primates.     

Genetic diversity of human immunodeficiency virus type 2: evidence for distinct sequence subtypes with differences in virus biology.

The virulence properties of human immunodeficiency virus type 2 (HIV-2) are known to vary significantly and to range from relative attenuation in certain individuals to high-level pathogenicity in others. These differences in clinical manifestations may, at least in part, be determined by genetic differences among infecting virus strains. Evaluation of the full spectrum of HIV-2 genetic diversity is thus a necessary first step towards understanding its molecular epidemiology, natural history of infection, and biological diversity. In this study, we have used nested PCR techniques to amplify viral sequences from the DNA of uncultured peripheral blood mononuclear cells from 12 patients with HIV-2 seroreactivity. Sequence analysis of four nonoverlapping genomic regions allowed a comprehensive analysis of HIV-2 phylogeny. The results revealed (i) the existence of five distinct and roughly equidistant evolutionary lineages of HIV-2 which, by analogy with HIV-1, have been termed sequence subtypes A to E; (ii) evidence for a mosaic HIV-2 genome, indicating that coinfection with genetically divergent strains and recombination can occur in HIV-2-infected individuals; and (iii) evidence supporting the conclusion that some of the HIV-2 subtypes may have arisen from independent introductions of genetically diverse sooty mangabey viruses into the human population. Importantly, only a subset of HIV-2 strains replicated in culture: all subtype A viruses grew to high titers, but attempts to isolate representatives of subtypes C, D, and E, as well as the majority of subtype B viruses, remained unsuccessful. Infection with all five viral subtypes was detectable by commercially available serological (Western immunoblot) assays, despite intersubtype sequence differences of up to 25% in the gag, pol, and env regions. These results indicate that the genetic and biological diversity of HIV-2 is far greater than previously appreciated and suggest that there may be subtype-specific differences in virus biology. Systematic natural history studies are needed to determine whether this heterogeneity has clinical relevance and whether the various HIV-2 subtypes differ in their in vivo pathogenicity.     

Nucleotide sequence of human endogenous retrovirus genome related to the mouse mammary tumor virus genome.

We determined the complete nucleotide sequence of the human endogenous retrovirus genome HERV-K10 isolated as the sequence homologous to the Syrian hamster intracisternal A-particle (type A retrovirus) genome. HERV-K10 is 9,179 base pairs long with long terminal repeats of 968 base pairs at both ends; a sequence 290 base pairs long, however, was found to be deleted. It was concluded that a composite genome having the 290-base-pair fragment is the prototype HERV-K provirus gag (666 codons), protease (334 codons), pol (937 codons), and env (618 codons) genes. The size of the protease gene product of HERV-K is essentially the same as that of A- and D-type oncoviruses but nearly twice that of other retroviruses. A comparison of the deduced amino acid sequences encoded by the pol region showed HERV-K to be closely related to types A and D retroviruses and even more so to type B retrovirus. It was noted that the env gene product of HERV-K structurally resembles the mouse mammary tumor virus (type B retrovirus) env protein, and the possible expression of the HERV-K env gene in human breast cancer cells is discussed.     

Different models, different trees: the geographic origin of PTLV-I

Using nucleotide sequences from three genomic regions of the human and simian T-cell lymphotropic virus type I (HTLV-I/STLV-I)-consisting of 69 sequences from a 140-bp segment of the pol region, 98 sequences from a 503-bp segment of the LTR, and 154 sequences from a 386-bp segment of the env region-we tested two hypotheses concerning the geographic origin and evolution of STLV-I and HTLV-I. First, we tested the assumption of equal rates of evolution along STLV-I and HTLV-I lineages using a likelihood ratio test to ascertain whether current levels of genomic diversity can be used to determine ancestry. We demonstrated that unequal rates of evolution along HTLV-I and STLV-I lineages have occurred throughout evolutionary time, thus calling into question the use of pairwise distances to assign ancestry. Second, we constructed phylogenetic trees using multiple phylogenetic techniques to test for the geographic origin of STLV-I and HTLV-I. Using the principle of likelihood, we chose a statistically justified model of evolution for each data set. We demonstrated the utility of the likelihood ratio test to determine which model of evolution should be chosen for phylogenetic analyses, revealing that using different models of evolution produces conflicting results, and neither the hypothesis of an African origin nor the hypothesis of an Asian origin can be rejected statistically. Our best estimates of phylogenetic relationships, however, support an African origin of PTLV for each gene region.     

Discovery of a New Endogenous Type C Retrovirus (FcEV) in Cats: Evidence for RD-114 Being an FcEVGag-Pol/Baboon Endogenous Virus BaEVEnv Recombinant

Analysis of a cat genomic DNA library showed that cats harbor a previously unrecognized endogenous type C retrovirus, whose env gene has homology to the murine Fv-4 resistance gene. This unique retrovirus, designated FcEV (Felis catus endogenous retrovirus), has a type C pol gene, closely related to the primate Papio cynocephalus endogenous virus (PcEV) pol, not overlapping the env gene, unlike in other type C retroviruses, and is presumably present in a higher copy number than RD-114. Phylogenetic analysis of FcEV and RD-114 fragments amplified from cat species and comparison with baboon endogenous virus (BaEV) fragments from monkeys suggested that RD-114 does not represent the cat strain of BaEV but is actually a new recombinant between FcEV type C genes and the env gene of BaEV. Although BaEV did appear to have infected an ancestor of the domestic cat lineage, it was a de novo recombinant that made its way into the cat germ line.     

Protease gene structure and env gene variability of the AIDS virus

The protease gene structure and the env gene variability have been precisely compared between the AIDS virus and members of the HTLV/BLV family. The conserved amino acid sequence (LVDT) which is repeated in the proteases of the HTLV/BLV family is not repeated in AIDS virus. Comparative analysis of the env gene sequences reveals the striking fact that the env gene of AIDS virus is 8-12-times more variable than those of the HTLV/BLV family. Within the AIDS virus env gene, the surface glycoprotein region is more liable to vary than is the transmembrane region; unexpectedly, however, this liability is not a characteristic feature of the AIDS virus because it is more prominent in other retroviruses including members of the HTLV/BLV family.