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.     

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.     

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.     

Concerted evolution of members of the multisequence family chAB4 located on various nonhomologous chromosomes

During the last years it became obvious that a lot of families of long-range repetitive DNA elements are located within the genomes of mammals. The principles underlying the evolution of such families, therefore, may have a greater impact than anticipated on the evolution of the mammalian genome as a whole. One of these families, called chAB4, is represented with about 50 copies within the human and the chimpanzee genomes and with only a few copies in the genomes of gorilla, orang-utan, and gibbon. Members of chAB4 are located on 10 different human chromosomes. FISH of chAB4-specific probes to chromosome preparations of the great apes showed that chAB4 is located, with only one exception, at orthologous places in the human and the chimpanzee genome. About half the copies in the human genome belong to two species-specific subfamilies that evolved after the divergence of the human and the chimpanzee lineages. The analysis of chAB4-specific PCR-products derived from DNA of rodent/human cell hybrids showed that members of the two human-specific subfamilies can be found on 9 of the 10 chAB4-carrying chromosomes. Taken together, these results demonstrate that the members of DNA sequence families can evolve as a unit despite their location at multiple sites on different chromosomes. The concerted evolution of the family members is a result of frequent exchanges of DNA sequences between copies located on different chromosomes. Interchromosomal exchanges apparently take place without greater alterations in chromosome structure. Received: 20 March 1997 / Accepted: 13 September 1997     

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.     

Characterization of a Novel Class of Interspersed LTR Elements in Primate Genomes: Structure, Genomic Distribution, and Evolution

Retrovirus-like sequences and their solitary (solo) long terminal repeats (LTRs) are common repetitive elements in eukaryotic genomes. We reported previously that the tandemly arrayed genes encoding U2 snRNA (the RNU2 locus) in humans and apes contain a solo LTR (U2-LTR) which was presumably generated by homologous recombination between the two LTRs of an ancestral provirus that is retained in the orthologous baboon RNU2 locus. We have now sequenced the orthologous U2-LTRs in human, chimpanzee, gorilla, orangutan, and baboon and examined numerous homologs of the U2-LTR that are dispersed throughout the human genome. Although these U2-LTR homologs have been collectively referred to as LTR13 in the literature, they do not display sequence similarity to any known retroviral LTRs; however, the structure of LTR13 closely resembles that of other retroviral LTRs with a putative promoter, polyadenylation signal, and a tandemly repeated 53-bp enhancer-like element. Genomic blotting indicates that LTR13 is primate-specific; based on sequence analysis, we estimate there are about 2,500 LTR13 elements in the human genome. Comparison of the primate U2-LTR sequences suggests that the homologous recombination event that gave rise to the solo U2-LTR occurred soon after insertion of the ancestral provirus into the ancestral U2 tandem array. Phylogenetic analysis of the LTR13 family confirms that it is diverse, but the orthologous U2-LTRs form a coherent group in which chimpanzee is closest to the humans; orangutan is a clear outgroup of human, chimpanzee, and gorilla; and baboon is a distant relative of human, chimpanzee, gorilla, and orangutan. We compare the LTR13 family with other known LTRs and consider whether these LTRs might play a role in concerted evolution of the primate RNU2 locus. Received: 29 September 1997 / Accepted: 16 January 1998     

The evolutionary analysis of the Tnt1 retrotransposon in Nicotiana species reveals the high variability of its regulatory sequences

We studied the evolution of the tobacco Tnt1 retrotransposon by analyzing Tnt1 partial sequences containing both coding domains and U3 regulatory sequences obtained from a number of Nicotiana species. We detected three different subfamilies of Tnt1 elements, Tnt1A, Tnt1B, and Tnt1C, that differ completely in their U3 regions but share conserved flanking coding and LTR regions. U3 divergence between the three subfamilies is found in the region that contains the regulatory sequences that control the expression of the well-characterized Tnt1-94 element. This suggests that expression of the three Tnt1 subfamilies might be differently regulated. The three Tnt1 subfamilies were present in the Nicotiana genome at the time of species divergence, but have evolved independently since then in the different genomes. Each Tnt1 subfamily seems to have conserved its ability to transpose in a limited and different number of Nicotiana species. Our results illustrate the high variability of Tnt1 regulatory sequences. We propose that this high sequence variability could allow these elements to evolve regulatory mechanisms in order to optimize their coexistence with their host genome.      

Comprehensive analysis of two Alu Yd subfamilies

Alu elements have inserted in the human genome throughout primate evolution. A small number of Alu insertions have occurred after the divergence of humans from nonhuman primates and therefore should not be present in nonhuman primate genomes. Most of these recently integrated Alu elements are contained with a series of discrete Alu subfamilies that are related to each other based upon diagnostic nucleotide substitutions. We have extracted members of the Alu Yd subfamily that are derivatives of the Alu Y subfamily that share a common 12-bp deletion that defines the Yd lineage from the draft sequence of the human genome. Analysis of the Yd Alu elements resulted in the recovery of two new Alu subfamilies, Yd3 and Yd6, which contain a total of 295 members (198 Yd3 and 97 Yd6). DNA sequence analysis of each of the Alu Yd subfamilies yielded age estimates of 8.02 and 1.20 million years old for the Alu Yd3 and Yd6 subfamilies, respectively. Two hundred Alu Yd3 and Yd6 loci were screened using polymerase chain reaction (PCR) assays to determine their phylogenetic origin and associated levels of human genomic diversity. The Alu Yd3 subfamily appears to have started amplifying relatively early in primate evolution and continued propagating albeit at a low level as many of its members are found in a variety of hominoid (humans, greater and lesser ape) genomes. Only two of the elements are polymorphic in the human genome and absent from the genomes of nonhuman primates. By contrast all of the members of the Alu Yd6 subfamily are restricted to the human genome, with 12% of the elements representing insertion polymorphisms in human populations. A single Alu Yd6 locus contained an independent parallel forward insertion of a paralogous Alu Sq sequence in the owl monkey. These Alu subfamilies are a source of genomic fossil relics for the study of primate phylogenetics and human population genetics.     

Evolutionary history of the human endogenous retrovirus family ERV9

Several distinct families of endogenous retrovirus-like elements (ERVs) exist in the genomes of primates. Despite the important evolutionary consequences that carrying these intragenomic parasites may have for their hosts, our knowledge about their evolution is still scarce. A matter of particular interest is whether evolution of ERVs occurs via a master lineage or through several lineages coexisting over long periods of time. In this work, the paleogenomic approach has been applied to the study of the evolution of ERV9, one of the human endogenous retrovirus families mobilized during primate evolution. By searching the GenBank database with the first 676 bp of the ERV9 long terminal repeat, we identified 156 different element insertions into the human genome. These elements were grouped into 14 subfamilies based on several characteristic nucleotide differences. The age of each subfamily was roughly estimated based on the average sequence divergence of its members from the subfamily consensus sequence. Determination of the sequential order of diagnostic substitutions led to the identification of four distinct lineages, which retained their capacity of transposition over extended periods of evolution. Strong evidence for mosaic evolution of some of these lineages is presented. Taken altogether, the available data indicate that the possibility of ERV9 still being active in the human lineage can not be discarded.     

SVA elements: a hominid-specific retroposon family

SVA is a composite repetitive element named after its main components, SINE, VNTR and Alu. We have identified 2762 SVA elements from the human genome draft sequence. Genomic distribution analysis indicates that the SVA elements are enriched in G+C-rich regions but have no preferences for inter- or intragenic regions. A phylogenetic analysis of the elements resulted in the recovery of six subfamilies that were named SVA_A to SVA_F. The composition, age and genomic distribution of the subfamilies have been examined. Subfamily age estimates based upon nucleotide divergence indicate that the expansion of four SVA subfamilies (SVA_A, SVA_B, SVA_C and SVA_D) began before the divergence of human, chimpanzee and gorilla, while subfamilies SVA_E and SVA_F are restricted to the human lineage. A survey of human genomic diversity associated with SVA_E and SVA_F subfamily members showed insertion polymorphism frequencies of 37.5% and 27.6%, respectively. In addition, we examined the amplification dynamics of SVA elements throughout the primate order and traced their origin back to the beginnings of hominid primate evolution, approximately 18 to 25 million years ago. This makes SVA elements the youngest family of retroposons in the primate order.     

Evolution of the master Alu gene(s)

Summary A comparison of Alu sequences that comprise more recently amplified Alu subfamilies was made. There are 18 individual diagnostic mutations associated with the different subfamilies. This analysis confirmed that the formation of each subfamily can be explained by the sequential accumulation of mutations relative to the previous subfamily. Polymerase chain reaction amplification of orthologous loci in several primate species allowed us to determine the time of insertion of Alu sequences in individual loci. These data suggest that the vast majority of Alu elements amplified at any given time comprised a single Alu subfamily. We find that, although the individual divergence relative to a consensus sequence correlate reasonably well with sequence age, the diagnostic mutations are a more accurate measure of the age of any individual Alu family member. Our data are consistent with a model in which all Alu family members have been made from a single master gene or from a series of sequential master genes. This master gene(s) accumulated diagnostic base changes, resulting in the amplification of different subfamilies from the master gene at different times in primate evolution. The changes in the master gene(s) probably occurred individually, but their appearance is clearly punctuated. Ten of them have occurred within an 15-million-year time span, 40–25 million years ago, and 8 changes have occurred within the last 5 million years. Surprisingly, no changes appeared in the 20 milion years separating these periods.     

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.     

HERV-K(OLD): ancestor sequences of the human endogenous retrovirus family HERV-K(HML-2)

Sequences homologous to the human endogenous retrovirus (HERV) family HERV-K(HML-2) are present in all Old World primate species. A previous study showed that a central region of the HERV-K(HML-2) gag genes in Hominoidea species displays a 96-bp deletion compared to the gag genes in lower Old World primates. The more ancient HERV-K(HML-2) sequences present in lower Old World primates were apparently not conserved during hominoid evolution, as opposed to the deletion variants. To further clarify the evolutionary origin of the HERV-K(HML-2) family, we screened GenBank with the 96-bp gag-sequence characteristic of lower Old World primates and identified, to date, 10 human sequence entries harboring either full-length or partially deleted proviral structures, probably representing remnants of a more ancient HERV-K(HML-2) variant. The high degree of mutations demonstrates the long-time presence of these HERV-K(OLD) proviruses in the genome. Nevertheless, they still belong to the HML-2 family as deduced from dot matrix and phylogenetic analyses. We estimate, based on the family ages of integrated Alu elements and on long terminal repeat (LTR) divergence data, that the average age of HERV-K(OLD) proviruses is ca. 28 million years, supporting an integration time before the evolutionary split of Hominoidea from lower Old World primates. Analysis of HERV-K(OLD) LTR sequences led to the distinction of two subgroups, both of which cluster with LTRs belonging to an evolutionarily older cluster. Taken together, our data give further insight into the evolutionary history of the HERV-K(HML-2) family during primate evolution.     

Human endogenous retrovirus family HERV-K(HML-5): status, evolution, and reconstruction of an ancient betaretrovirus in the human genome

The human genome harbors numerous distinct families of so-called human endogenous retroviruses (HERV) which are remnants of exogenous retroviruses that entered the germ line millions of years ago. We describe here the hitherto little-characterized betaretrovirus HERV-K(HML-5) family (named HERVK22 in Repbase) in greater detail. Out of 139 proviruses, only a few loci represent full-length proviruses, and many lack gag protease and/or env gene regions. We generated a consensus sequence from multiple alignment of 62 HML-5 loci that displays open reading frames for the four major retroviral proteins. Four HML-5 long terminal repeat (LTR) subfamilies were identified that are associated with monophyletic proviral bodies, implying different evolution of HML-5 LTRs and genes. Sequence analysis indicated that the proviruses formed approximately 55 million years ago. Accordingly, HML-5 proviral sequences were detected in Old World and New World primates but not in prosimians. No recent activity is associated with this HERV family. We also conclude that the HML-5 consensus sequence primer binding site is identical to methionine tRNA. Therefore, the family should be designated HERV-M. Our study provides important insights into the structure and evolution of the oldest betaretrovirus in the primate genome known to date.     

The Cassandra retrotransposon landscape in sugar beet (Beta vulgaris) and related Amaranthaceae: recombination and re-shuffling lead to a high structural variability

Background and AimsPlant genomes contain many retrotransposons and their derivatives, which are subject to rapid sequence turnover. As non-autonomous retrotransposons do not encode any proteins, they experience reduced selective constraints leading to their diversification into multiple families, usually limited to a few closely related species. In contrast, the non-coding Cassandra terminal repeat retrotransposons in miniature (TRIMs) are widespread in many plants. Their hallmark is a conserved 5S rDNA-derived promoter in their long terminal repeats (LTRs). As sugar beet (Beta vulgaris) has a well-described LTR retrotransposon landscape, we aim to characterize TRIMs in beet and related genomes.MethodsWe identified Cassandra retrotransposons in the sugar beet reference genome and characterized their structural relationships. Genomic organization, chromosomal localization, and distribution of Cassandra-TRIMs across the Amaranthaceae were verified by Southern and fluorescent in situ hybridization.Key resultsAll 638 Cassandra sequences in the sugar beet genome contain conserved LTRs and thus constitute a single family. Nevertheless, variable internal regions required a subdivision into two Cassandra subfamilies within B. vulgaris. The related Chenopodium quinoa harbours a third subfamily. These subfamilies vary in their distribution within Amaranthaceae genomes, their insertion times and the degree of silencing by small RNAs. Cassandra retrotransposons gave rise to many structural variants, such as solo LTRs or tandemly arranged Cassandra retrotransposons. These Cassandra derivatives point to an interplay of template switch and recombination processes – mechanisms that likely caused Cassandra’s subfamily formation and diversification.ConclusionsWe traced the evolution of Cassandra in the Amaranthaceae and detected a considerable variability within the short internal regions, whereas the LTRs are strongly conserved in sequence and length. Presumably these hallmarks make Cassandra a prime target for unequal recombination, resulting in the observed structural diversity, an example of the impact of LTR-mediated evolutionary mechanisms on the host genome.     

Characterization of the intragenomic spread of the human endogenous retrovirus family HERV-W

This study examines the intragenomic spread of the human endogenous retrovirus family HERV-W from insertions present within the draft sequence of the human genome. Identification of shared diagnostic differences and phylogenetic analyses revealed the existence of three main subfamilies. The average divergence between sequences for each of the subfamilies suggests that most of the HERV-W elements were inserted within the genome during a short period of evolutionary time. Each one of the subfamilies consists of two types of insertions, the expected proviral sequences and other sequences resembling the structure of processed retrogenes. These HERV-W retrosequences extend from the R region of the 5' long-terminal repeat (LTR) to the R region of the 3' LTR (as viral genomic RNAs), end in poly(A) 3' tails, and are flanked by direct repeats longer than the proviral integrations. Furthermore, several of the HERV-W retrosequences are 5'-truncated at different sites. I suggest the involvement of the L1 machinery in these integrations and discuss the characteristic features of the evolutionary history of HERV-W, with emphasis on the putative impact of HERV-W retrosequence integrations on the mammalian genome.     

Human-Specific Integrations of the HERV-K Endogenous Retrovirus Family

Several distinct families of endogenous retrovirus-like sequences (HERVs) exist in the genomes of humans and other primates. One of these families, the HERV-K group, contains members that encode functional proteins and that have been implicated in the etiology of insulin-dependent diabetes mellitus (IDDM). Because of potential functional and disease relevance, it is important to determine if there are HERV-K-associated genetic differences between individuals. In this study, we have investigated the divergence and evolutionary age of HERV-K long terminal repeats (LTRs). Thirty-seven LTRs, taken primarily from random human clones in GenBank, were aligned and grouped into nine clusters with decreasing sequence divergence. Cluster 1 sequences are 8.6% divergent, on average, whereas cluster 9 LTRs, represented by the LTRs of the fully sequenced HERV-K10 clone, show an average of only 1.1% divergence from each other. The evolutionary age of 18 LTRs from different clusters was then investigated by genomic PCR to determine presence or absence of the retroviral element in different primate species. LTRs from clusters of higher divergence were detected in monkeys and apes, whereas LTRs in clusters with lower divergence were acquired later in evolution. Notably, LTRs of cluster 9 were found only in humans at all nine loci examined. Genomic Southern analysis with an oligonucleotide probe specific for cluster 9 LTRs suggests that HERV-K elements with this type of LTR expanded independently in the genomes of humans and the great apes. This is the first report of endogenous retroviral integrations that are specific to humans and indicates that some HERVs have amplified much later than previously thought. These elements may still be actively transposing and may therefore represent a source of genetic variation linked to disease development.