A comprehensive analysis of recently integrated human Ta L1 elements

The Ta (transcribed, subset a) subfamily of L1 LINEs (long interspersed elements) is characterized by a 3-bp ACA sequence in the 3' untranslated region and contains approximately 520 members in the human genome. Here, we have extracted 468 Ta L1Hs (L1 human specific) elements from the draft human genomic sequence and screened individual elements using polymerase-chain-reaction (PCR) assays to determine their phylogenetic origin and levels of human genomic diversity. One hundred twenty-four of the elements amenable to complete sequence analysis were full length ( approximately 6 kb) and have apparently escaped any 5' truncation. Forty-four of these full-length elements have two intact open reading frames and may be capable of retrotransposition. Sequence analysis of the Ta L1 elements showed a low level of nucleotide divergence with an estimated age of 1.99 million years, suggesting that expansion of the L1 Ta subfamily occurred after the divergence of humans and African apes. A total of 262 Ta L1 elements were screened with PCR-based assays to determine their phylogenetic origin and the level of human genomic variation associated with each element. All of the Ta L1 elements analyzed by PCR were absent from the orthologous positions in nonhuman primate genomes, except for a single element (L1HS72) that was also present in the common (Pan troglodytes) and pygmy (P. paniscus) chimpanzee genomes. Sequence analysis revealed that this single exception is the product of a gene conversion event involving an older preexisting L1 element. One hundred fifteen (45%) of the Ta L1 elements were polymorphic with respect to insertion presence or absence and will serve as identical-by-descent markers for the study of human evolution.     

Straightening out the LINEs: LINE-1 orthologous loci

The L1Hs preTa subfamily of long interspersed elements (LINEs) originated after the divergence of human and chimpanzee and is therefore found only in the human genome. Thirty-three of the 254 L1Hs preTa elements are polymorphic for the absence/presence of the insertion, making them useful markers for studying human population genetics. The problem of homoplasy, however, can diminish the value of LINEs as phylogenetic and population genetic markers. We examined anomalous orthologous sites in a range of nonhuman primates. Only two cases of other mobile elements inserting near the preintegration sites of L1Hs preTa elements were observed: an AluY insertion in Chlorocebus and an L1PA8 insertion in Aotus. Sequence analysis showed that both elements were clearly distinguishable from their human counterparts. We conclude that L1 elements can continue to be regarded as essentially homoplasy-free genetic characters.     

Non-traditional Alu evolution and primate genomic diversity

Alu elements belonging to the previously identified "young" subfamilies are thought to have inserted in the human genome after the divergence of humans from non-human primates and therefore should not be present in non-human primate genomes. Polymerase chain reaction (PCR) based screening of over 500 Alu insertion loci resulted in the recovery of a few "young" Alu elements that also resided at orthologous positions in non-human primate genomes. Sequence analysis demonstrated these "young" Alu insertions represented gene conversion events of pre-existing ancient Alu elements or independent parallel insertions of older Alu elements in the same genomic region. The level of gene conversion between Alu elements suggests that it may have a significant influence on the single nucleotide diversity within the genome. All the instances of multiple independent Alu insertions within the same small genomic regions were recovered from the owl monkey genome, indicating a higher Alu amplification rate in owl monkeys relative to many other primates. This study suggests that the majority of Alu insertions in primate genomes are the products of unique evolutionary events.     

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.     

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.     

Analysis of the human Alu Ya-lineage

The Alu Ya-lineage is a group of related, short interspersed elements (SINEs) found in primates. This lineage includes subfamilies Ya1-Ya5, Ya5a2 and others. Some of these subfamilies are still actively mobilizing in the human genome. We have analyzed 2482 elements that reside in the human genome draft sequence and focused our analyses on the 2318 human autosomal Ya Alu elements. A total of 1470 autosomal loci were subjected to polymerase chain reaction (PCR)-based assays that allow analysis of individual Ya-lineage Alu elements. About 22% (313/1452) of the Ya-lineage Alu elements were polymorphic for the insertion presence on human autosomes. Less than 0.01% (5/1452) of the Ya-lineage loci analyzed displayed insertions in orthologous loci in non-human primate genomes. DNA sequence analysis of the orthologous inserts showed that the orthologous loci contained older pre-existing Y, Sc or Sq Alu subfamily elements that were the result of parallel forward insertions or involved in gene conversion events in the human lineage. This study is the largest analysis of a group of "young", evolutionarily related human subfamilies. The size, evolutionary age and variable allele insertion frequencies of several of these subfamilies makes members of the Ya-lineage useful tools for human population studies and primate phylogenetics.     

Recently integrated human Alu repeats: finding needles in the haystack

Alu elements undergo amplification through retroposition and integration into new locations throughout primate genomes. Over 500,000 Alu elements reside in the human genome, making the identification of newly inserted Alu repeats the genomic equivalent of finding needles in the haystack. Here, we present two complementary methods for rapid detection of newly integrated Alu elements. In the first approach we employ computational biology to mine the human genomic DNA sequence databases in order to identify recently integrated Alu elements. The second method is based on an anchor-PCR technique which we term Allele-Specific Alu PCR (ASAP). In this approach, Alu elements are selectively amplified from anchored DNA generating a display or 'fingerprint' of recently integrated Alu elements. Alu insertion polymorphisms are then detected by comparison of the DNA fingerprints generated from different samples. Here, we explore the utility of these methods by applying them to the identification of members of the smallest previously identified subfamily of Alu repeats in the human genome termed Ya8. This subfamily of Alu repeats is composed of about 50 elements within the human genome. Approximately 50% of the Ya8 Alu family members have inserted in the human genome so recently that they are polymorphic, making them useful markers for the study of human evolution. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

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.     

Recently integrated Alu elements and human genomic diversity

A comprehensive analysis of two Alu Y lineage subfamilies was undertaken to assess Alu-associated genomic diversity and identify new Alu insertion polymorphisms for the study of human population genetics. Recently integrated Alu elements (283) from the Yg6 and Yi6 subfamilies were analyzed by polymerase chain reaction (PCR), and 25 of the loci analyzed were polymorphic for insertion presence/absence within the genomes of a diverse array of human populations. These newly identified Alu insertion polymorphisms will be useful tools for the study of human genomic diversity. Our screening of the Alu insertion loci also resulted in the recovery of several "young" Alu elements that resided at orthologous positions in nonhuman primate genomes. Sequence analysis demonstrated these "young" Alu insertions were the products of gene conversion events of older, preexisting Alu elements or independent parallel forward insertions of older Alu elements in the same short genomic region. The level of gene conversion between Alu elements suggests that it may have an influence on the single nucleotide polymorphism within Alu elements in the genome. We have also identified two genomic deletions associated with the retroposition and insertion of Alu Y lineage elements into the human genome. This type of Alu retroposition-mediated genomic deletion is a novel source of lineage-specific evolution within primate genomes.     

Analysis of newly identified low copy AluYj subfamily

Human specific AluY elements were investigated by comparative analysis between human chromosome 21 and chimpanzee chromosome 22. Human specific AluY element was identified on human chromosome 21q22 (accession no. AL163282), and then that was a new member of AluYj subfamily. From the bioinformatic analysis, AluYj subfamily was investigated in human whole genome using AluYj4 consensus sequence (accession no. AL163282). Thirteen members of the AluYj4 elements (4 diagnostic mutations) and eight members of the AluYj3 elements (3 diagnostic mutations) were identified with distinct diagnostic mutation from AluY consensus sequence. The results of the molecular clock calculation of non-CpG region substitution indicated that, AluYj4 elements (2.1 million years old) may be proliferated more recent time than AluYj3 elements (14.1 million years old). For the verification of recent insertion time, four of AluYj4 elements (ch2-AC017101, ch10-AC044786, ch12-AC007656 and ch21-AL163282) from human chromosomes 2, 10, 12, 21 were analyzed by PCR amplification using various human and primate DNA samples. Though, no polymorphism was detected in human population, we identified the new AluYj4 subfamily as the human specific elements.     

Identification and characterization of novel polymorphic LINE-1 insertions through comparison of two human genome sequence assemblies

Mobile elements represent a relatively new class of markers for the study of human evolution. Long interspersed elements (LINEs) belong to a group of retrotransposons comprising approximately 21% of the human genome. Young LINE-1 (L1) elements that have integrated recently into the human genome can be polymorphic for insertion presence/absence in different human populations at particular chromosomal locations. To identify putative novel L1 insertion polymorphisms, we computationally compared two draft assemblies of the whole human genome (Public and Celera Human Genome assemblies). We identified a total of 148 potential polymorphic L1 insertion loci, among which 73 were candidates for novel polymorphic loci. Based on additional analyses we selected 34 loci for further experimental studies. PCR-based assays and DNA sequence analysis were performed for these 34 loci in 80 unrelated individuals from four diverse human populations: African-American, Asian, Caucasian, and South American. All but two of the selected loci were confirmed as polymorphic in our human population panel. Approximately 47% of the analyzed loci integrated into other repetitive elements, most commonly older L1s. One of the insertions was accompanied by a BC200 sequence. Collectively, these mobile elements represent a valuable source of genomic polymorphism for the study of human population genetics. Our results also suggest that the exhaustive identification of L1 insertion polymorphisms is far from complete, and new whole genome sequences are valuable sources for finding novel retrotransposon insertion polymorphisms.     

The Alu Yc1 subfamily: sorting the wheat from the chaff

Members of the Alu Yc1 subfamily are distinguished from the older Alu Y subfamily by a signature G-->A substitution at base 148 of their 281-bp consensus sequence. Members of the much older and larger Alu Y subfamily could have by chance accumulated this signature G-->A substitution and be misclassified as belonging to the Alu Yc1 subfamily. Using a Mahanalobis classification method, it was estimated that the "authentic" Alu Yc1 subfamily consists of approximately 262 members in the human genome. PCR amplification and further analysis was successfully completed on 225 of the Yc1 Alu family members. One hundred and seventy-seven Yc1 Alu elements were determined to be monomorphic (fixed for presence) in a panel of diverse human genomes. Forty-eight of the Yc1 Alu elements were polymorphic for insertion presence/absence in diverse human genomes. The insertion polymorphism rate of 21% in the human genome is similar to rates reported previously for other "young" Alu subfamilies. The polymorphic Yc1 Alu elements will be useful genetic loci for the study of human population genetics.     

A new 5' sequence associated with mouse L1 elements is representative of a major class of L1 termini.

Full-length L1 elements have been shown to possess, at their 5' end, tandem repeats called "A" or "F" types. By sequencing the 5' region of two large L1 copies that did not hybridize to A or F probes, we have identified a new sequence that is found at the 5' end of many L1 elements and that we call "V." The element characterized has no 200-bp tandem repetitive structure, and the new 5' sequence is not similar to the A or F sequences. The study of the relationships between the V and L1 sequences has shown that only half of the V (i.e., V-specific 5') sequences in the genome are linked to the 5' end of L1 copies. In related rodent species, a comparative study by Southern blot and PCR analysis of the V sequence suggests that this L1 subfamily has an ancient origin and that V sequence isolated from the remainder of the L1 element has been amplified during the evolution of the mouse genome.     

Genomic rearrangements by LINE-1 insertion-mediated deletion in the human and chimpanzee lineages

Long INterspersed Elements (LINE-1s or L1s) are abundant non-LTR retrotransposons in mammalian genomes that are capable of insertional mutagenesis. They have been associated with target site deletions upon insertion in cell culture studies of retrotransposition. Here, we report 50 deletion events in the human and chimpanzee genomes directly linked to the insertion of L1 elements, resulting in the loss of ~18 kb of sequence from the human genome and ~15 kb from the chimpanzee genome. Our data suggest that during the primate radiation, L1 insertions may have deleted up to 7.5 Mb of target genomic sequences. While the results of our in vivo analysis differ from those of previous cell culture assays of L1 insertion-mediated deletions in terms of the size and rate of sequence deletion, evolutionary factors can reconcile the differences. We report a pattern of genomic deletion sizes similar to those created during the retrotransposition of Alu elements. Our study provides support for the existence of different mechanisms for small and large L1-mediated deletions, and we present a model for the correlation of L1 element size and the corresponding deletion size. In addition, we show that internal rearrangements can modify L1 structure during retrotransposition events associated with large deletions.     

An effort to use human-based exome capture methods to analyze chimpanzee and macaque exomes

Non-human primates have emerged as an important resource for the study of human disease and evolution. The characterization of genomic variation between and within non-human primate species could advance the development of genetically defined non-human primate disease models. However, non-human primate specific reagents that would expedite such research, such as exon-capture tools, are lacking. We evaluated the efficiency of using a human exome capture design for the selective enrichment of exonic regions of non-human primates. We compared the exon sequence recovery in nine chimpanzees, two crab-eating macaques and eight Japanese macaques. Over 91% of the target regions were captured in the non-human primate samples, although the specificity of the capture decreased as evolutionary divergence from humans increased. Both intra-specific and inter-specific DNA variants were identified; Sanger-based resequencing validated 85.4% of 41 randomly selected SNPs. Among the short indels identified, a majority (54.6%-77.3%) of the variants resulted in a change of 3 base pairs, consistent with expectations for a selection against frame shift mutations. Taken together, these findings indicate that use of a human design exon-capture array can provide efficient enrichment of non-human primate gene regions. Accordingly, use of the human exon-capture methods provides an attractive, cost-effective approach for the comparative analysis of non-human primate genomes, including gene-based DNA variant discovery.     

L1 repeat elements in the human epsilon-G gamma-globin gene intergenic region: sequence analysis and concerted evolution within this family

We have deduced the sequence of a composite long interspersed repeated DNA in primates and herein describe its relationship to a complex repeat element (L1Heg) located in the interval linking the human epsilon- and G gamma-globin genes. The main element of L1Heg is 3' truncated and interrupted by the insertion of the 3' end of a second L1 element. Transposition of L1Heg into this intergenic locus generated a 62-bp duplication of flanking sequences. In contrast, insertion of the second repeat may have been mediated by homology between donor and target sequences. The main repeat represents a novel class of abundant elements whose sequences have diverged from other rodent and primate LINES approximately 1.3 kb downstream from the 5' terminus of L1Heg. Comparison of L1Heg with the sequences of two other related L1 members revealed a complex set of rearrangements confined within a region that resembles the long terminal repeats of other types of retroposons. The boundaries of conversion-like events were defined on the basis of the clustering of nucleotide sequence variants common to two or more nonallelic 3' L1H elements. Several of these events are apparently initiated or resolved within a common 150-bp region that coincides with the 3' terminus of a pan-mammalian open reading frame. This analysis showed that concerted genetic interactions and random drift both contribute appreciably to sequence variation within this set of L1H members.     

Phylogenetic relationships among transposon-like elements in human and primate DNA

A THE-1 sequence in intron 7 of the human dystrophin gene has been found to represent a new subfamily of THE-1 elements. The sequence is closely related to the MstII family of repetitive sequences and is more like single-copy sequences found in the galago genome than any other THE-1 sequence previously reported. This new THE-1 sequence has been compared with two other complete THE-1 sequences and three related long-terminal repeat elements that we have previously found in intron 7 of the dystrophin gene, and with members of the same family from elsewhere in the primate genome. Parsimony and deletion analysis show that the cluster of THE-1 sequences in intron 7 of the dystrophin gene has arisen from at least three individual insertion events, rather than from the insertion and duplication of a single progenitor sequence. Correspondence to: G.B. Petersen