Genome-wide targeted search for human specific and polymorphic L1 integrations

Retroelements (REs) occupy up to 40% of the human genome. Newly integrated REs can change the pattern of expression of pre-existing host genes and therefore might play a significant role in evolution. In particular, human- and primate-specific REs could affect the divergence of the Hominoidea superfamily. A comparative genome-wide analysis of RE sites of integration, neighboring genes, and their regulatory interplay in human and ape genomes would be of help in understanding the impact of REs on evolution and genome regulation. We have developed a technique for the genome-wide comparison of the integrations of transposable elements in genomic DNAs of closely related species. The technique called targeted genome differences analysis (TGDA) is also useful for the detection of deletion/insertion polymorphisms of REs. The technique is based on an enhanced version of subtractive hybridization and does not require preliminary knowledge of the genome sequences under comparison. In this report, we describe its application to the detection and analysis of human specific L1 integrations and their polymorphisms. We obtained a library highly enriched in human-specific L1 insertions and identified 24 such new insertions. Many of these insertions are polymorphic in human populations. The total number of human-specific L1 inserts was estimated to be approximately 4000. The results suggest that TGDA is a universal method that can be successfully used for the detection of evolutionary and polymorphic markers in any closely related genomes.     

Genome-wide comparison of differences in the integration sites of interspersed repeats between closely related genomes

A technique for genome-wide detection of differences in the integration site positions of interspersed repeats in related genomes (DiffIR) is described. The technique is based on a whole- genome selective PCR amplification of the repeats’ flanking regions followed by a differential hybridization screening of the arrayed library of the selected amplicons. The technique was successfully applied to the comparison of the integration sites in the human and chimpanzee genomes, allowing us to discover 11 new human-specific integrations of human endogenous retrovirus, K family (HML-2) long terminal repeats.     

Endogenous retroviruses and human evolution

Humans share about 99% of their genomic DNA with chimpanzees and bonobos; thus, the differences between these species are unlikely to be in gene content but could be caused by inherited changes in regulatory systems. Endogenous retroviruses (ERVs) comprise approximately 5% of the human genome. The LTRs of ERVs contain many regulatory sequences, such as promoters, enhancers, polyadenylation signals and factor-binding sites. Thus, they can influence the expression of nearby human genes. All known human-specific LTRs belong to the HERV-K (human ERV) family, the most active family in the human genome. It is likely that some of these ERVs could have integrated into regulatory regions of the human genome, and therefore could have had an impact on the expression of adjacent genes, which have consequently contributed to human evolution. This review discusses possible functional consequences of ERV integration in active coding regions.     

Identification of large-scale human-specific copy number differences by inter-species array comparative genomic hybridization

Copy number differences (CNDs), and the concomitant differences in gene number, have contributed significantly to the genomic divergence between humans and other primates. To assess its relative importance, the genomes of human, common chimpanzee, bonobo, gorilla, orangutan and macaque were compared by comparative genomic hybridization using a high-resolution human BAC array (aCGH). In an attempt to avoid potential interference from frequent intra-species polymorphism, pooled DNA samples were used from each species. A total of 322 sites of large-scale inter-species CND were identified. Most CNDs were lineage-specific but frequencies differed considerably between the lineages; the highest CND frequency among hominoids was observed in gorilla. The conserved nature of the orangutan genome has already been noted by karyotypic studies and our findings suggest that this degree of conservation may extend to the sub-microscopic level. Of the 322 CND sites identified, 14 human lineage-specific gains were observed. Most of these human-specific copy number gains span regions previously identified as segmental duplications (SDs) and our study demonstrates that SDs are major sites of CND between the genomes of humans and other primates. Four of the human-specific CNDs detected by aCGH map close to the breakpoints of human-specific karyotypic changes [e.g., the human-specific inversion of chromosome 1 and the polymorphic inversion inv(2)(p11.2q13)], suggesting that human-specific duplications may have predisposed to chromosomal rearrangement. The association of human-specific copy number gains with chromosomal breakpoints emphasizes their potential importance in mediating karyotypic evolution as well as in promoting human genomic diversity. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Following the LINEs: an analysis of primate genomic variation at human-specific LINE-1 insertion sites

The L1 Ta subfamily of long interspersed elements (LINEs) consists exclusively of human-specific L1 elements. Polymerase chain reaction-based screening in nonhuman primate genomes of the orthologous sites for 249 human L1 Ta elements resulted in the recovery of various types of sequence variants for approximately 12% of these loci. Sequence analysis was employed to capture the nature of the observed variation and to determine the levels of gene conversion and insertion site homoplasy associated with LINE elements. Half of the orthologous loci differed from the predicted sizes due to localized sequence variants that occurred as a result of common mutational processes in ancestral sequences, often including regions containing simple sequence repeats. Additional sequence variation included genomic deletions that occurred upon L1 insertion, as well as successive mobile element insertions that accumulated within a single locus over evolutionary time. Parallel independent mobile element insertions at orthologous loci in distinct species may introduce homoplasy into retroelement-based phylogenetic and population genetic data. We estimate the overall frequency of parallel independent insertion events at L1 insertion sites in seven different primate species to be very low (0.52%). In addition, no cases of insertion site homoplasy involved the integration of a second L1 element at any of the loci, but rather largely involved secondary insertions of Alu elements. No independent mobile element insertion events were found at orthologous loci in the human and chimpanzee genomes. Therefore, L1 insertion polymorphisms appear to be essentially homoplasy free characters well suited for the study of population genetics and phylogenetic relationships within closely related species.     

Molecular cloning and phylogenetic analysis of the human endogenous retrovirus HERV-K long terminal repeat elements in various cancer cells

Long terminal repeats (LTRs) of the human endogenous retroviruses K family (HERV-K) have been found to affect expression of genes located nearby. It has been suggested that the HERV-K LTR elements contributed to the structural change in the genome and genetic variation connected to various diseases. We examined the HERV-K LTR elements in human cancer cells. Using genomic DNA from various cancer cells, we performed PCR amplification and identified forty-nine HERV-K LTR elements. Those LTR elements showed a high degree of sequence similarity with human-specific HERV-K LTR elements. A phylogenetic tree, obtained by the neighbor-joining method, revealed that twelve HERV-K LTR elements were closely related to human-specific HERV-K LTR elements. These elements proliferated recently and were detectable in many human cancer cell lines. These results suggest that HERV-K LTR could be implicated in a pathogenic role, although this phenomenon may not directly lead to human cancers. Further studies on the biological function and expression of HERV-K LTR elements in cancer cells are indicated.     

Full-sized HERV-K (HML-2) human endogenous retroviral LTR sequences on human chromosome 21: map locations and evolutionary history

One of the evolutionary mechanisms for acquisition of novel functional sequences can be domestication of exogenous retroviruses that have been integrated into the germ line. The whole genome mapping of such elements in various species could reveal differences in positions of the retroviral integration and suggest possible roles of these differences in speciation. Here, we describe the number, locations and sequence features of the human endogenous retrovirus HERV-K (HML-2) long terminal repeat (LTR) sequences on human chromosome 21. We show that their distribution along the chromosome is not only non-random but also roughly correlated with the gene density. Amplification of orthologous LTR sites from a number of primate genomes produced patterns of presence and absence for each LTR sequence and allowed determination of the phylogenetic ages and evolutionary order of appearance of individual LTRs. The identity level and phylogenetic age of the LTRs did not correlate with their map locations. Thus, despite the non-random distribution of LTRs, they have apparently been inserted randomly into the chromosome relative to each other. As evidenced in previous studies of chromosomes 19 and 22, this is a characteristic of HERV-K integration.     

Differences in HERV-K LTR insertions in orthologous loci of humans and great apes

The classification of the long terminal repeats (LTRs) of the human endogenous retrovirus HERV-K (HML-2) family was refined according to diagnostic differences between the LTR sequences. The mutation rate was estimated to be approximately equal for LTRs belonging to different families and branches of human endogenous retroviruses (HERVs). An average mutation rate value was calculated based on differences between LTRs of the same HERV and was found to be 0.13% per million years (Myr). Using this value, the ages of different LTR groups belonging to the LTR HML-2 subfamily were found to vary from 3 to 50Myr. Orthologous potential LTR-containing loci from different primate species were PCR amplified using primers corresponding to the genomic sequences flanking LTR integration sites. This allowed us to calculate the phylogenetic times of LTR integrations in primate lineages in the course of the evolution and to demonstrate that they are in good agreement with the LTR ages calculated from the mutation rates. Human-specific integrations for some very young LTRs were demonstrated. The possibility of LTRs and HERVs involvement in the evolution of primates is discussed.     

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.     

Nucleotide sequence and phylogenetic analysis of long terminal repeats of human endogenous retrovirus K family (HERV-K) on human chromosomes

It has been suggested that human endogenous retroviruses K family (HERV-K) has a role in disease, and solitary long terminal repeats (LTRs) of HERV-K have been potentially capable of affecting the expression of closely located genes. Using the human monochromosomes 8, 9, 17, and 18, with specific PCR primers, we identified thirty-four sequences of new HERV-K LTRs. Those LTR elements were analyzed phylogenetically with the human-specific HERV-K LTRs using neighbor-joining and maximum parsimony methods. Clones HKL8-5, HKL9-5, and HKL9-8 are related by more than 99% homology with the human-specific HERV-K LTRs. The HKL9-5 clone on chromosome 9 was 100% identical with the sequences of human-specific LTR, AC002400, on chromosome 16. The findings suggest that there has been recent proliferation, transposition, or chromosomal translocation of HERV-K LTR elements on human chromosomes.     

Human-specific sequences: isolation of species-specific DNA regions by genome subtraction

Uniqueness is fundamental to the individuality of species, and this in turn is based on the uniqueness of their genomes. For the purpose of resolving the genetic basis of human uniqueness, we describe here the isolation of human-specific sequences using the technique of genome subtraction, i.e., competitive reassociation of genomic DNAs between two very closely related species. One such sequence, HS5, was found to be present only in the human genome and absent in the genomes of non-human primates including chimpanzees, the species most closely related to humans.     

High polymorphism level of genomic sequences flanking insertion sites of human endogenous retroviral long terminal repeats

The polymorphism at the multitude of loci adjacent to human endogenous retrovirus long terminal repeats (LTRs) was analyzed by a technique for whole genome differential display based on the PCR suppression effect that provides selective amplification and display of genomic sequences flanking interspersed repeated elements. This strategy is simple, target-specific, requires a small amount of DNA and provides reproducible and highly informative data. The average frequency of polymorphism observed in the vicinity of the LTR insertion sites was found to be about 12%. The high incidence of polymorphism within the LTR flanks together with the frequent location of LTRs near genes makes the LTR loci a useful source of polymorphic markers for gene mapping.     

Formation of a new solo-LTR of the human endogenous retrovirus H family in human chromosome 21

Human endogenous retroviruses (HERVs) contribute to various kinds of genomic instability via rearrangement and retrotransposition events. In the present study the formation of a new human-specific solo-LTR belonging to the HERV-H family (AP001667; chromosome 21q21) was detected by a comparative analysis of human chromosome 21 and chimpanzee chromosome 22. The solo-LTR was formed as a result of an equal homologous recombination excision event. Several evolutionary processes have occurred at this locus during primate evolution, indicating that mammalian-wide interspersed repeat (MIR) and full-length HERV-H elements integrated into hominoid genomes after the divergence of Old World monkeys and hominoids, and that the solo-LTR element was created by recombination excision of the HERV-H only in the human genome.     

Human-Specific HERV-K Insertion Causes Genomic Variations in the Human Genome

Human endogenous retroviruses (HERV) sequences account for about 8% of the human genome. Through comparative genomics and literature mining, we identified a total of 29 human-specific HERV-K insertions. We characterized them focusing on their structure and flanking sequence. The results showed that four of the human-specific HERV-K insertions deleted human genomic sequences via non-classical insertion mechanisms. Interestingly, two of the human-specific HERV-K insertion loci contained two HERV-K internals and three LTR elements, a pattern which could be explained by LTR-LTR ectopic recombination or template switching. In addition, we conducted a polymorphic test and observed that twelve out of the 29 elements are polymorphic in the human population. In conclusion, human-specific HERV-K elements have inserted into human genome since the divergence of human and chimpanzee, causing human genomic changes. Thus, we believe that human-specific HERV-K activity has contributed to the genomic divergence between humans and chimpanzees, as well as within the human population.