Human endogenous retrovirus K106 (HERV-K106) was infectious after the emergence of anatomically modern humans

HERV-K113 and HERV-K115 have been considered to be among the youngest HERVs because they are the only known full-length proviruses that are insertionally polymorphic and maintain the open reading frames of their coding genes. However, recent data suggest that HERV-K113 is at least 800,000 years old, and HERV-K115 even older. A systematic study of HERV-K HML2 members to identify HERVs that may have infected the human genome in the more recent evolutionary past is lacking. Therefore, we sought to determine how recently HERVs were exogenous and infectious by examining sequence variation in the long terminal repeat (LTR) regions of all full-length HERV-K loci. We used the traditional method of inter-LTR comparison to analyze all full length HERV-Ks and determined that two insertions, HERV-K106 and HERV-K116 have no differences between their 5' and 3' LTR sequences, suggesting that these insertions were endogenized in the recent evolutionary past. Among these insertions with no sequence differences between their LTR regions, HERV-K106 had the most intact viral sequence structure. Coalescent analysis of HERV-K106 3' LTR sequences representing 51 ethnically diverse individuals suggests that HERV-K106 integrated into the human germ line approximately 150,000 years ago, after the emergence of anatomically modern humans.     

Molecular characterization of the HERV-W env gene in humans and primates: expression, FISH, phylogeny, and evolution

We characterized the human endogenous retrovirus (HERV-W) family in humans and primates. In silico expression data indicated that 22 complete HERV-W families from human chromosomes 1-3, 5-8, 10-12, 15, 19, and X are randomly expressed in various tissues. Quantitative real-time RT-PCR analysis of the HERV-W env gene derived from human chromosome 7q21.2 indicated predominant expression in the human placenta. Several copies of repeat sequences (SINE, LINE, LTR, simple repeat) were detected within the complete or processed pseudo HERV-W of the human, chimpanzee, and rhesus monkey. Compared to other regions (5'LTR, Gag, Gag-Pol, Env, 3'LTR), the repeat family has been mainly integrated into the region spanning the 5'LTRs of Gag (1398 bp) and Pol (3242 bp). FISH detected the HERV-W probe (fosWE1) derived from a gorilla fosmid library in the metaphase chromosomes of all primates (five hominoids, three Old World monkeys, two New World monkeys, and one prosimian), but not in Tupaia. This finding was supported by molecular clock and phylogeny data using the divergence values of the complete HERV-W LTR elements. The data suggested that the HERV-W family was integrated into the primate genome approximately 63 million years (Myr) ago, and evolved independently during the course of primate radiation.     

Negative regulatory element associated with potentially functional promoter and enhancer elements in the long terminal repeats of endogenous murine leukemia virus-related proviral sequences.

Three series of recombinant DNA clones were constructed, with the bacterial chloramphenicol acetyltransferase (CAT) gene as a quantitative indicator, to examine the activities of promoter and enhancer sequence elements in the 5' long terminal repeat (LTR) of murine leukemia virus (MuLV)-related proviral sequences isolated from the mouse genome. Transient CAT expression was determined in mouse NIH 3T3, human HT1080, and mink CCL64 cultured cells transfected with the LTR-CAT constructs. The 700-base-pair (bp) LTRs of three polytropic MuLV-related proviral clones and the 750-bp LTRs of four modified polytropic proviral clones, in complete structures either with or without the adjacent downstream sequences, all showed very little or negligible activities for CAT expression, while ecotropic MuLV LTRs were highly active. The MuLV-related LTRs were divided into three portions and examined separately. The 3' portion of the MuLV-related LTRs that contains the CCAAC and TATAA boxes was found to be a functional promoter, being about one-half to one-third as active as the corresponding portion of ecotropic MuLV LTRs. A MboI-Bg/II fragment, representing the distinct 190- to 200-bp inserted segment in the middle, was found to be a potential enhancer, especially when examined in combination with the simian virus 40 promoter in CCL64 cells. A PstI-MboI fragment of the 5' portion, which contains the protein-binding motifs of the enhancer segment as well as the upstream LTR sequences, showed moderate enhancer activities in CCL6 cells but was virtually inactive in NIH 3T3 cells and HT1080 cells; addition of this fragment to the ecotropic LTR-CAT constructs depressed CAT expression. Further analyses using chimeric LTR constructs located the presence of a strong negative regulatory element within the region containing the 5' portion of the enhancer and the immediate upstream sequences in the MuLV-related LTRs.     

Human-specific subfamilies of HERV-K (HML-2) long terminal repeats: three master genes were active simultaneously during branching of hominoid lineages

Using 40 known human-specific LTR sequences, we have derived a consensus sequence for an evolutionary young HERV-K (HML-2) LTR family, which was named the HS family. In the human genome the HS family is represented by approximately 150-160 LTR sequences, 90% of them being human-specific (hs). The family can be subdivided into two subfamilies differing in five linked nucleotide substitutions: HS-a and HS-b of 5.8 and 10.3 Myr evolutionary ages, respectively. The HS-b subfamily members were transpositionally active both before the divergence of the human and chimpanzee ancestor lineages and after it in both lineages. The HS-a subfamily comprises only hs LTRs. These and other data strongly suggest that at least three "master genes" of HERV-K (HML-2) LTRs were active in the human ancestor lineage after the human-chimpanzee divergence. We also found hs HERV-K (HML-2) LTRs integrations in introns of 12 human genes and identified 13 new hs HERV-K (HML-2) LTRs.     

Molecular phylogeny and evolution of the human endogenous retrovirus HERV-W LTR family in hominoid primates

Long terminal repeats (LTRs) of human endogenous retrovirus (HERV) have contributed to the structural change or genetic variation of primate genome that are connected to speciation and evolution. Using genomic DNAs that were derived from hominoid primates (chimpanzee, gorilla, orangutan, and gibbon), we performed PCR amplification and identified thirty HERV-W LTR elements. These LTR elements showed a 82-98% sequence similarity with HERV-W LTR (AF072500). Specifically, additional sequences (GCCACCACCACTGTTT in the gorilla and TGCTGCTGACTCCCATCC in the gibbon) were noticed. Clone OR3 from the orangutan and clone GI2 from the gibbon showed a 100% sequence similarity, although they are different species. This indicates that both LTR elements were proliferated during the last 2 to 5 million years from the integration of the original LTR element. A phylogenetic tree that was obtained by the neighbor-joining method revealed a wide overlap of the LTR elements across species, suggesting that the HERV-W LTR family evolved independently during the hominoid evolution.     

The R-U5-5'' leader sequence of neurovirulent wild mouse retrovirus contains an element controlling the incubation period of neurodegenerative disease.

The wild mouse ecotropic retrovirus CasBrE causes a spongiform neurodegenerative disease after neonatal inoculation, with an incubation period ranging from 2 to 12 months. We previously showed that introduction of long terminal repeat (LTR) and gag-pol sequences from a strain of Friend murine leukemia virus (FB29) resulted in a dramatic acceleration of the onset of the disease. The chimeric virus FrCasE, which consisted of the FB29 genome containing 3' pol and env sequences from the wild mouse virus, induced a highly predictable, lethal neurodegenerative disease with an incubation period of only 16 days. Here we report that the sequences which are primary determinants of the length of the incubation period are located in the 5' end of the viral genome between a KpnI site in the R region of the LTR and a PstI site immediately 5' of the start codon for pr65gag (R-U5-5' leader). This region contains the tRNA primer binding site, splice donor site for the subgenomic env mRNA, and the packaging sequence. Computer-assisted sequence analysis failed to find evidence of a consensus sequence for a DNA enhancer in this region. In addition, sequences within a region of the genome between a ClaI site at the 3' end of env to the KpnI site in the R region of the LTR (inclusive of U3) also influenced the incubation period of the disease, but the effect was distinctly weaker than that of the R-U5-5' leader sequence. This U3 effect, however, appeared to be independent of the number of direct repeats, since deletion of one of two duplicated 42-base repeats containing consensus sequences of nuclear-factor binding domains had no effect on the incubation period of the disease. On the basis of Southern blot analysis of total viral DNA in the tissues, the effect of these sequences on the incubation period appeared to be related to the level of virus replication in the central nervous system. All of the chimeric viruses analyzed, irrespective of neurovirulence, replicated to comparable levels in the spleen and induced comparable levels of viremia.     

Alterations in the U3 region of the long terminal repeat of an infectious thymotropic type B retrovirus.

We isolated and characterized a type B thymotropic retrovirus (DMBA-LV) which is highly related to mouse mammary tumor virus (MMTV) isolates and which induces T-cell thymomas with a high incidence and a very short latent period. Regions of nonhomology between the DMBA-LV genome and the MMTV genome were identified by heteroduplex mapping and nucleotide sequence studies. In the electron microscope heteroduplex mapping studies the EcoRI-generated 5' and 3' fragments of the DMBA-LV genome were compared with the corresponding fragments of the MMTV (C3H and GR) genome isolated from mammary tumors. The results indicated that DMBA-LV contained a region of nonhomologous nucleotide sequences in the 3' half of the U3 region of the long terminal repeat (LTR). Nucleotide sequence studies confirmed these results and showed that in this region 440 nucleotides of the MMTV (C3H) sequences were deleted and substituted with a segment of 122 nucleotides. This substituted segment in the form of a tandem repeat structure contained nucleotide sequences derived exclusively from sequences which flanked the substitution loop. The distal glucocorticoid regulatory element was unaltered, and two additional copies of the distal glucocorticoid regulatory element-binding site were present in the substituted region. The restriction endonuclease map of the reconstructed molecular clone of DMBA-LV was identical to that corresponding to unintegrated linear DMBA-LV DNA present in DMBA-LV-induced tumor cell lines. Since the nucleotide sequences of the LTRs present in four different DMBA-LV proviral copies isolated from a single thymoma were identical, we concluded that they were derived from the same parental virus and that this type B retrovirus containing an alteration in the U3 region of its LTR could induce thymic lymphomas. Thus, DMBA-LV represents the first example of a productively replicating type B retrovirus that contains an LTR modified in the U3 region and that has target cell and disease specificity for T cells.     

Characterization of proviruses cloned from mink cell focus-forming virus-infected cellular DNA

Two proviruses were cloned from EcoRI-digested DNA extracted from mink cells chronically infected with AKR mink cell focus-forming (MCF) 247 murine leukemia virus (MuLV), using a lambda phage host vector system. One cloned MuLV DNA fragment (designated MCF 1) contained sequences extending 6.8 kilobases from an EcoRI restriction site in the 5' long terminal repeat (LTR) to an EcoRI site located in the envelope (env) region and was indistinguishable by restriction endonuclease mapping for 5.1 kilobases (except for the EcoRI site in the LTR) from the 5' end of AKR ecotropic proviral DNA. The DNA segment extending from 5.1 to 6.8 kilobases contained several restriction sites that were not present in the AKR ecotropic provirus. A 0.5-kilobase DNA segment located at the 3' end of MCF 1 DNA contained sequences which hybridized to a xenotropic env-specific DNA probe but not to labeled ecotropic env-specific DNA. This dual character of MCF 1 proviral DNA was also confirmed by analyzing heteroduplex molecules by electron microscopy. The second cloned proviral DNA (designated MCF 2) was a 6.9-kilobase EcoRI DNA fragment which contained LTR sequences at each end and a 2.0-kilobase deletion encompassing most of the env region. The MCF 2 proviral DNA proved to be a useful reagent for detecting LTRs electron microscopically due to the presence of nonoverlapping, terminally located LTR sequences which effected its circularization with DNAs containing homologous LTR sequences. Nucleotide sequence analysis demonstrated the presence of a 104-base-pair direct repeat in the LTR of MCF 2 DNA. In contrast, only a single copy of the reiterated component of the direct repeat was present in MCF 1 DNA.     

A rare event of insertion polymorphism of a HERV-K LTR in the human genome

Human endogenous retroviruses (HERVs), which constitute a significant part of the human genome, might have a serious impact on primate evolution. Over a hundred insertions of HERV-K(HML-2) family members distinguish the human genome from other primate genomes. However, only three cases of insertion polymorphisms have been reported so far, all for endogenous HERV-K proviruses. This suggests that some retroviral integrations occurred rather recently in human genome evolution. In this report, we describe a very rare case of true insertion polymorphism of a solitary HERV-K LTR in the human genome. Distribution of the LTR-containing allele was tested in 5 Africans and 83 individuals from three Russian populations. The allele frequency appeared to be relatively high in populations of both European and Asian origin. The detected polymorphic LTR could be a useful molecular genetic marker of the corresponding genomic region.     

Identification of paralogous HERV-K LTRs on human chromosomes 3, 4, 7 and 11 in regions containing clusters of olfactory receptor genes

A locus harboring a human endogenous retroviral LTR (long terminal repeat) was mapped on the short arm of human chromosome 7 (7p22), and its evolutionary history was investigated. Sequences of two human genome fragments that were homologous to the LTR-flanking sequences were found in human genome databases: (1) an LTR-containing DNA fragment from region 3p13 of the human genome, which includes clusters of olfactory receptor genes and pseudogenes; and (2) a fragment of region 21q22.1 lacking LTR sequences. PCR analysis demonstrated that LTRs with highly homologous flanking sequences could be found in the genomes of human, chimp, gorilla, and orangutan, but were absent from the genomes of gibbon and New World monkeys. A PCR assay with a primer set corresponding to the sequence from human Chr 3 allowed us to detect LTR-containing paralogous sequences on human chromosomes 3, 4, 7, and 11. The divergence times for the LTR-flanking sequences on chromosomes 3 and 7, and the paralogous sequence on chromosome 21, were evaluated and used to reconstruct the order of duplication events and retroviral insertions. (1) An initial duplication event that occurred 14-17 Mya and before LTR insertion - produced two loci, one corresponding to that located on Chr 21, while the second was the ancestor of the loci on chromosomes 3 and 7. (2) Insertion of the LTR (most probably as a provirus) into this ancestral locus took place 13 Mya. (3) Duplication of the LTR-containing ancestral locus occurred 11 Mya, forming the paralogous modern loci on Chr 3 and 7.