Phylogeny of SINE-R retroposons in Asian apes.

The SINE-R retroposon family was derived from the long terminal repeats (LTRs) of human endogenous retrovirus K (HERV-K) that had been active during the hominoid evolution. The retroposons and HERV-K LTR elements have potential relevance to structural change and genetic variation of the hominoid genome. In our previous study, we found that the SINE-R retroposons were hominoid specific. Here we identified seventeen new SINE-R retroposons (14 from orangutan and 3 from gibbon) from Asian apes and phylogenetically analysed them in comparison with those of the humans and African great apes. None of the retroposons from Asian apes were closely related to SINE-R.C2 that is human specific, and originally identified in the gene for the C2 component of complement, whereas some retroposons (Ch-M10, Ch-M16, Gor-M, Gor-F1, Gor-M6, and Gor-F9) from African great apes showed very close relationship with that of the SINE-R.C2 retroposon. The phylogenetic tree based on the SINE-R retroposons revealed wide overlap of the retroposons across species, suggesting that the SINE-R retroposons have been evolved parallel pattern in the course of hominoid evolution.     

Phylogenetic analysis of retroposon family as exemplified on human chromosome 13: further evidence for recent proliferation

The retroposon SINE-R.C2 was first identified as a human-specific insertion in the complement C2 gene. In our previous study, SINE-R type retroposons, derived from the endogenous retrovirus HERV-K family, have been found to be hominoid specific. In this report on human chromosome 13, we identified eighteen new SINE-R retroposons resembling those we have previously reported on the sex chromosomes and on chromosomes 7 and 17. Phylogenetic analysis using the neighbor-joining method revealed that four SINE-R retroposons (13-16, 21, 23, 25) on chromosome 13 were closely related to the human-specific retroposon SINE-R.C2, with a high degree of sequence homology (95-97%). Such elements differ from the HERV-K10. LTR sequence from which they are derived in being deleted for the promoter region. Therefore while the evidence adds to the case that some classes of SINE-R element have continued to proliferate in hominid and hominoid evolution and may, as in the case of Fukuyama type muscular dystrophy, be a cause of insertional mutagenesis, they are less likely than the HERV-K10 LTR to have a positive effect on host gene activity.     

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     

Phylogenetic analysis of a retroposon family as represented on the human X chromosome

SINE-R elements constitute a class of retroposons derived from the long terminal repeat (LTR) of the human endogenous retrovirus HERV-K family that are present in hominoid primates and active in the human genome. In an investigation of the X chromosome, we identified twenty-five SINE-R elements with between 89.6 and 97.7% homology with the SINE-R.C2 element that is human specific, originally identified in the gene for the C2 component of complement. SINE-R.C2 and a sequence HS307 that we previously identified in a region of Xq21.3 that has a recently created homology with a 4 Mb block in Yp11.2 are amongst the group of elements that have diverged furthest from the parent HERV-K10 sequence. The sequence on the X chromosome resemble those that we previously described on chromosomes 7 and 17 and the Y chromosome, with a similar range of variation. Phylogenetic analysis from the retroposon family including those of African great apes using the neighbor-joining method suggests that the SINE-R retroposon family have evolved independently during primate evolution. Further investigation of SINE-R elements on the sex chromosomes, particularly in recently created regions of X-Y homology, may cast light on the timing of the retroposition process and its possible relevance to recent evolutionary change.     

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.     

Evolutionary dynamics of the human endogenous retrovirus family HERV-K inferred from full-length proviral genomes

Several distinct families of endogenous retroviruses exist in the genomes of primates. Most of them are remnants of ancient germ-line infections. The human endogenous retrovirus family HERV-K represents the unique known case of endogenous retrovirus that amplified in the human genome after the divergence of human and chimpanzee lineages. There are two types of HERV-K proviral genomes differing by the presence or absence of 292 bp in the pol-env boundary. Human-specific insertions exist for both types. The analyses shown in the present work reveal that several lineages of type 1 and type 2 HERV-K proviruses remained transpositionally active after the human/chimpanzee split. The data also reflect the important role of mosaic evolution (either by recombination or gene conversion) during the evolutionary history of HERV-K. Received: 5 February 2001 / Accepted: 22 March 2001     

Comparison Between Two Human Endogenous Retrovirus (HERV)-Rich Regions Within the Major Histocompatibility Complex

Sixteen human endogenous retrovirus (HERV) sequences were detected within 656 kb of genomic sequence obtained from the alpha- and beta-block of the class I region of the major histocompatibility complex (MHC). The HERVs were identified and characterized as family members of HERV-16 (11 copies), HERV-L (1 copy), HERV-I (2 copies), HERV-K91 (1 copy), and HARLEQUIN (1 copy) by sequence comparison using CENSOR or Repeat Masker, BLAST searches, and dot plots. The 11 copies of HERV-16 arose as products of duplication of genomic segments containing HLA class I (HLAcI) and PERB11 (MIC) genes inter alia, whereas the other five HERVs arose after duplication probably as a consequence of single insertion events or translocations. HERV-L and HERV-I are located between the duplicated genes PERB11.2 (MICB) and PERB11.1 (MICA), and HLA-B and HLA-C, respectively, whereas HERV-K91 and HARLEQUIN are located telomeric of HLA-C. A highly fragmented copy of HERV-I was also found telomeric of PERB11.4. Structural analysis of open reading frames (ORFs) revealed the absence of intact coding sequence within the putative gag, pol, and env gene regions of all the HERVs with the exception of HERV-K91, which had two large ORFs within the region of the putative protease and pol genes. In addition, the 5′-LTR of HERV-L contained a 2.5-kb element that was AT-rich and large ORFs with putative amino acid sequences rich in tyrosines and isoleucines. HERV-I, HARLEQUIN, and at least four copies of HERV-16 appear to have been receptors for the insertion of other retrotransposons including Alu elements and fragments of L1 and THE1. Examination of flanking sequences suggests that HERV-I and HERV-L had occurred by insertion into ancient L1 fragments. This study has revealed that the alpha- and beta-block region within the MHC is rich in HERV sequences occurring at a much higher ratio (10 to 1) than normally observed in the human genome. These HERV sequences will therefore enhance further studies on disease associations and differences between human haplotypes and primates and their role in the evolution of class I genes in the MHC. Received: 17 September 1998 / Accepted: 8 January 1999     

A novel human nonviral retroposon derived from an endogenous retrovirus.

In a human genome, we found dispersed repetitive sequences homologous to part of a human endogenous retrovirus termed HERV-K which resembled mouse mammary tumor virus. For elucidation of their structure and organization, we cloned some of these sequences from a human gene library. The sequence common to the cloned DNA was ca. 630 base-pairs (bp) in length with an A-rich tail at the 3' end and was found to be a SINE (short interspersed repeated sequence) type nonviral retroposon. In this retroposon, the 5' end had multiple copies of a 40 bp direct repeat very rich in GC content and about the next 510 nucleotides were homologous to the 3' long terminal repeat and its upstream flanking region of the HERV-K genome. This retroposon was thus given the name, SINE-R element since most of it derived from a retrovirus. SINE-R elements were present at 4,000 to 5,000 copies per haploid human genome. The nucleotide sequence was ca. 90% homologous among the cloned elements.     

Structure and Origin of a Novel Dimeric Retroposon B1-dID

Here we describe a new short retroposon family of rodents. Like the primate Alu element consisting of two similar monomers, it is dimeric, but the left and right monomers are different and descend from B1 and ID short retroposons, respectively. Such elements (B1-dID) were found in the genomes of Gliridae, Sciuridae, Castoridae, Caviidae, and Hystricidae. Nucleotide sequences of this retroposon can be assigned to several structural variants. Phylogenetic analysis of B1-dID and related sequences suggests a possible scenario of B1-dID evolution in the context of rodent evolution. Received: 30 August 1999 / Accepted: 20 March 2000     

Ancestry of SINE-R.C2 a human-specific retroposon

We have reported previously that a retroposon, containing a variable number of tandemly repeated nucleotide sequences, is present in the third intron of the human C2 gene. This element, termed SINE-R.C2, is a member of a large retroposon family derived from the endogenous retrovirus HERV-K10 and estimated to include a few thousand copies per haploid human genome. In the present study we analyzed genomic DNA from 175 humans from several ethnic groups including Americans of European and African descent, Chinese, Africans, Australians, Pacific Islanders, Japanese, and Koreans. They all contained SINE-R.C2, as indicated by Southern blotting. However, SINE-R.C2 was absent from the genome of nonhuman primates, although SINE-R-type elements were present in chimpanzees and gorillas and the HERV-K10 genome was apparently present in all primates except for New World monkeys. These results indicate that HERV-K10 was inserted into the genome after the divergence of New World monkeys; the prototype SINE-R element, after divergence of orangutans; and SINE-R.C2, after the split between humans and chimpanzees.     

Full-Length L1 Elements Have Arisen Recently in the Same 1-kb Region of the Gorilla and Human Genomes

New copies of the mammalian retrotransposon L1 arise in the germline at an undetermined rate. Each new L1 copy appears at a specific evolutionary time point that can be estimated by phylogenetic analysis. In humans, the active L1 sequence L1.2 resides at the genomic locus LRE1. Here we analyzed the region surrounding the LRE1 locus in humans and gorillas to determine the evolutionary history of the region and to estimate the age of L1.2. We found that the region was composed of an ancient L1, L1Hs-Lrg, which was significantly divergent from all other L1 sequences available in the databases. We also determined that L1.2 was absent from the gorilla genome and arose in humans after the divergence of gorilla and human lineages. In the gorilla LRE1 region, we discovered a different full-length L1 element, L1Gg-1, which was allelic and present at a high gene frequency in gorillas but absent from other primates. We determined the nucleotide sequence of L1Gg-1 and found that it was 98% identical to L1.2, suggesting a close relationship between active L1s in gorillas and humans. Received: 28 December 1997 / Accepted: 20 March 1998     

Evolution of the glycophorin gene family in the hominoid primates

Analysis of nucleotide sequences of the human glycophorin A (GPA) and glycophorin B (GPB) genes has indicated that the GPA gene most closely resembles the ancestral gene, whereas the GPB gene likely arose from the GPA gene by homologous recombination. To study the evolution of the glycophorin gene family in the hominoid primates, restricted DNA on Southern blots from man, pygmy chimpanzee, common chimpanzee, gorilla, orangutan, and gibbon was probed with cDNA fragments encoding the human GPA and GPB coding and 3-untranslated regions. This showed the presence in all of the hominoid primates of at least one GPA-like gene. In addition, at least one GPB-like gene was detected in man, both chimpanzee species, and gorilla, strongly suggesting that the event that produced the GPB gene occurred in the common ancestor of man-chimpanzee-gorilla. An unexpected finding in this study was the conservation ofEcoRI restriction sites relative to those of the other four enzymes used; the significance of this observation is unclear, but raises the question of nonrandomness ofEcoRI restriction sites in noncoding regions. Further analysis of the evolution of this multigene family, including nucleotide sequence analysis, will be useful in clarification of the evolutionary relationships of the hominoid primates, in correlation with the structure and function of the glycophorin molecules, and in assessment of the role of evolution in the autogenicity of glycophorin determinants.This work was supported in part by National Institutes of Health Grants AM33463 and CA33000.     

Ancient Origin of the Null Allele se 428 of the Human ABO-Secretor Locus (FUT2)

In human populations, a null allele having several nucleotide differences from the wild-type allele is segregating at the FUT2 locus (the ABO-Secretor locus) encoding α(1,2)fucosyltransferase. To estimate the age of the most recent common ancestor (MRCA) of these two alleles, we sequenced FUT2 homologues from chimpanzee, gorilla, orangutan, and green monkey. Since we did not detect acceleration or any heterogeneity in the substitution rate at this locus among these species, the age of the MRCA was estimated to be around 3 MYA, assuming the divergence time of human and chimpanzee to be 5 MYA. We developed a simple test to examine whether or not the old age of the MRCA of the FUT2 is consistent with that expected for two divergent neutral alleles sampled from a random mating population. An application of the test to the data at FUT2 indicated that the age of the MRCA is too old to be explained by the simple neutral assumptions, although our test depends on accurate estimation of the divergence time of human and chimpanzee in units of twice the human population size. Various possibilities including balancing selection are discussed to explain this old age of the MRCA. Received: 9 May 1999 / Accepted: 20 September 1999     

Using Alu J Elements as Molecular Clocks to Trace the Evolutionary Relationships Between Duplicated HLA Class I Genomic Segments

The class I region of the major histocompatibility complex contains two subgenomic blocks (250–350 kb each), known as the alpha and beta blocks. These blocks contain members of multicopy gene families including HLA class I, HERV-16 (previously called P5 sequences), and PERB11 (MIC). We have previously shown that each block consists of imperfect duplicated segments (duplicons) containing linked members of different gene families, retroelements and transposons that have coevolved as part of two separate evolutionary events. Another region provisionally designated here as the kappa block is located between the alpha and the beta blocks and contains HLA-E, -30, and -92, HERV-16 (P5.3), and PERB11.3 (MICC) within about 250 kb of sequence. Using Alu elements to trace the evolutionary relationships between different class I duplicons, we have found that (a) the kappa block contains paralogous (duplicated) Alu J sequences and other retroelement patterns more in common with the beta than the alpha block; (b) the retroelement pattern associated with the HLA-E duplicon is different from all other HLA class I duplicons, indicating a more complex evolution; (c) the HLA-92 duplicon, although substantially shorter, is closely related in sequence to the HLA-B and -C duplicons; (d) two of the six paralogous Alu J elements within the HLA-B and -C duplicons are associated with the HLA-X duplicon, confirming their evolutionary relationships within the beta block; and (e) the paralogous Alu J elements within the alpha block are distinctly different from those identified within the beta and kappa blocks. The sequence conservation and location of duplicated (paralogous) Alu J elements in the MHC class I region show that the beta and kappa blocks have evolved separately from the alpha block beginning at a time before or during the evolution of Alu J elements in primates. Received: 22 September 1999 / Accepted: 24 January 2000     

Concerted evolution of the primate immunoglobulin alpha-gene through gene conversion.

We determined four nucleotide sequences of the hominoid immunoglobulin alpha (C alpha) genes (chimpanzee C alpha 2, gorilla C alpha 2, and gibbon C alpha 1 and C alpha 2 genes), which made possible the examination of gene conversions in all hominoid C alpha genes. The following three methods were used to detect gene conversions: 1) phenetic tree construction; 2) detection of a DNA segment with extremely low variability between duplicated C alpha genes; and 3) a site by site search of shared nucleotide changes between duplicated C alpha genes. Results obtained from method 1 indicated a concerted evolution of the duplicated C alpha genes in the human, chimpanzee, gorilla, and gibbon lineages, while results obtained from method 2 suggested gene conversions in the human, gorilla, and gibbon C alpha genes. With method 3 we identified clusters of shared nucleotide changes between duplicated C alpha genes in human, chimpanzee, gorilla, and gibbon lineages, and in their hypothetical ancestors. In the present study converted regions were identified over the entire C alpha gene region excluding a few sites in the coding region which have escaped from gene conversion. This indicates that gene conversion is a general phenomenon in evolution, that can be clearly observed in non-functional regions.     

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.     

Genomic and Phylogenetic Analysis of the Human CD1 and HLA Class I Multicopy Genes

The human CD1 proteins belong to a lipid-glycolipid antigen-presenting gene family and are related in structure and function to the MHC class I molecules. Previous mapping and DNA hybridization studies have shown that five linked genes located within a cluster on human chromosome 1q22-23 encode the CD1 protein family. We have analyzed the complete genomic sequence of the human CD1 gene cluster and found that the five active genes are distributed over 175,600 nucleotides and separated by four expanded intervening genomic regions (IGRs) ranging in length between 20 and 68 kb. The IGRs are composed mostly of retroelements including five full-length L1 PA sequences and various pseudogenes. Some L1 sequences have acted as receptors for other subtypes or families of retroelements. Alu molecular clocks that have evolved during primate history are found distributed within the HLA class I duplicated segments (duplicons) but not within the duplicons of CD1. Phylogeny of the alpha3 domain of the class I-like superfamily of proteins shows that the CD1 cluster is well separated from HLA class I by a number of superfamily members including MIC (PERB11), HFE, Zn-alpha2-GP, FcRn, and MR1. Phylogenetically, the human CD1 sequences are interspersed by CD1 sequences from other mammalian species, whereas the human HLA class I sequences cluster together and are separated from the other mammalian sequences. Genomic and phylogenetic analyses support the view that the human CD1 gene copies were duplicated prior to the evolution of primates and the bulk of the HLA class I genes found in humans. In contrast to the HLA class I genomic structure, the human CD1 duplicons are smaller in size, they lack Alu clocks, and they are interrupted by IGRs at least 4 to 14 times longer than the CD1 genes themselves. The IGRs seem to have been created as "buffer zones" to protect the CD1 genes from disruption by transposable elements.