Large-scale analysis of the Alu Ya5 and Yb8 subfamilies and their contribution to human genomic diversity

We have utilized computational biology to screen GenBank for the presence of recently integrated Ya5 and Yb8 Alu family members. Our analysis identified 2640 Ya5 Alu family members and 1852 Yb8 Alu family members from the draft sequence of the human genome. We selected a set of 475 of these elements for detailed analyses. Analysis of the DNA sequences from the individual Alu elements revealed a low level of random mutations within both subfamilies consistent with the recent origin of these elements within the human genome. Polymerase chain reaction assays were used to determine the phylogenetic distribution and human genomic variation associated with each Alu repeat. Over 99 % of the Ya5 and Yb8 Alu family members were restricted to the human genome and absent from orthologous positions within the genomes of several non-human primates, confirming the recent origin of these Alu subfamilies in the human genome. Approximately 1 % of the analyzed Ya5 and Yb8 Alu family members had integrated into previously undefined repeated regions of the human genome. Analysis of mosaic Yb8 elements suggests gene conversion played an important role in generating sequence diversity among these elements. Of the 475 evaluated elements, a total of 106 of the Ya5 and Yb8 Alu family members were polymorphic for insertion presence/absence within the genomes of a diverse array of human populations. The newly identified Alu insertion polymorphisms will be useful tools for the study of human genomic diversity.     

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.     

Mobility of short interspersed repeats within the chimpanzee lineage

The PV subfamily of Alu repeats in human DNA is largely composed of recently inserted members. Here we document additional members of the PV subfamily that are found in chimpanzee but not in the orthologous loci of human and gorilla, confirming the relatively recent and independent expansion of this Alu subfamily in the chimpanzee lineage. As further evidence for the youth of this Alu subfamily, one PV Alu repeat is specific to Pan troglodytes, whereas others are present in Pan paniscus as well. The A-rich tails of these Alu repeats have different lengths in Pan paniscus and Pan troglodytes. The dimorphisms caused by the presence and absence of PV Alu repeats and the length polymorphisms attributed to their A-rich tails should provide valuable genetic markers for molecular-based studies of chimpanzee relationships. The existence of lineage-specific Alu repeats is a major sequence difference between human and chimpanzee DNAs. Correspondence to: C.W. Schmid     

Alu Evolution in Human Populations: Using the Coalescent to Estimate Effective Population Size

There are estimated to be ~1000 members of the Ya5 Alu subfamily of retroposons in humans. This subfamily has a distribution restricted to humans, with a few copies in gorillas and chimpanzees. Fifty-seven Ya5 elements were previously cloned from a HeLa-derived randomly sheared total genomic library, sequenced, and screened for polymorphism in a panel of 120 unrelated humans. Forty-four of the 57 cloned Alu repeats were monomorphic in the sample and 13 Alu repeats were dimorphic for insertion presence/absence. The observed distribution of sample frequencies of the 13 dimorphic elements is consistent with the theoretical expectation for elements ascertained in a single diploid cell line. Coalescence theory is used to compute expected total pedigree branch lengths for monomorphic and dimorphic elements, leading to an estimate of human effective population size of ~18,000 during the last one to two million years.     

A human-specific subfamily of Alu sequences.

Of a total of 500,000 Alu family members, approximately 500 are present as a human-specific (HS) subfamily. Each of the HS subfamily members shares a high degree of nucleotide identity and is not present at orthologous positions in other primate genomes, suggesting that HS subfamily members have recently inserted within the human genome. This confirms the hypothesis that the majority of Alu family members are amplified copies of a "master" gene(s). This master gene appears to be amplifying at a rate much slower than that seen earlier in primate evolution. Some of the HS Alu subfamily members have amplified so recently that they are dimorphic in the human population, making them a potentially powerful tool for studies of human populations.     

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.     

Retrotransposition of Alu elements: how many sources?

It is generally thought that only a few Alu elements are capable of retrotransposition and that these 'master' sources produce inactive copies. Here, we use a network phylogenetic approach to demonstrate that recently integrated human-specific Alu subfamilies typically contain 10-20% of secondary source elements that contributed 20-40% of all subfamily members. This multiplicity of source elements provides new insight into the remarkably successful amplification strategy of the Alu family.     

In search of polymorphic Alu insertions with restricted geographic distributions

Alu elements are transposable elements that have reached over one million copies in the human genome. Some Alu elements inserted in the genome so recently that they are still polymorphic for insertion presence or absence in human populations. Recently, there has been an increasing interest in using Alu variation for studies of human population genetic structure and inference of individual geographic origin. Currently, this requires a high number of Alu loci. Here, we used a linker-mediated polymerase chain reaction method to preferentially identify low-frequency Alu elements in various human DNA samples with different geographic origins. The candidate Alu loci were subsequently genotyped in 18 worldwide human populations (approximately 370 individuals), resulting in the identification of two new Alu insertions restricted to populations of African ancestry. Our results suggest that it may ultimately become possible to correctly infer the geographic affiliation of unknown samples with high levels of confidence without having to genotype as many as 100 Alu loci. This is desirable if Alu insertion polymorphisms are to be used for human evolution studies or forensic applications.     

Alu repeats in the human genome

Highly repetitive DNA sequences account for more than 50% of the human genome. The L1 and Alu families harbor the most common mammalian long (LINEs) and short (SINEs) interspersed elements. Alu elements are each a dimer of similar, but not identical, fragments of total size about 300 bp, and originate from the 7SL RNA gene. Each element contains a bipartite promoter for RNA polymerase III, a poly(A) tract located between the monomers, a 3'-terminal poly(A) tract, and numerous CpG islands, and is flanked by short direct repeats. Alu repeats comprise more than 10% of the human genome and are capable of retroposition. Possibly, these elements played an important part in genome evolution. Insertion of an Alu element into a functionally important genome region or other Alu-dependent alterations of gene functions cause various hereditary disorders and are probably associated with carcinogenesis. In total, 14 Alu families differing in diagnostic mutations are known. Some of these, which are present in the human genome, are polymorphic and relatively recently inserted into new loci. Alu copies transposed during ethnic divergence of the human population are useful markers for evolutionary genetic studies.     

Human DNA quantitation using Alu element-based polymerase chain reaction

Human forensic casework requires sensitive quantitation of human nuclear DNA from complex sources. Widely used commercially available systems detect both nonhuman and human primate DNA, often require special equipment, and have a detection limit of approximately 0.1ng. Multicopy Alu elements include recently integrated subfamilies that are present in the human genome but are largely absent from nonhuman primates. Here, we present two Alu element-based alternative methods for the rapid identification and quantitation of human DNA, inter-Alu PCR and intra-Alu PCR. Using SYBR green-based detection, the effective minimum threshold level for human DNA quantitation was 0.01ng using inter-Alu- and 0.001ng using intra-Alu-based PCR. Background cross-amplification with nonhuman DNA templates was detected at low levels using inter-Alu-based PCR, but was negligible using intra-Alu-based PCR. These Alu-based methods have several advantages over currently available systems. First, the assays are PCR based and no additional unique equipment is required. Second, the high copy number of subfamily-specific Alu repeats in the human genome makes these assays human specific within a very sensitive linear range. The introduction of these assays to forensic laboratories will undoubtedly increase the sensitivity and specificity of human DNA detection and quantitation from complex sources.     

Phylogenetic affinities of tarsier in the context of primate Alu repeats

Related genomes tend to be colonized by the same or similar repetitive sequence elements. Analysis of these elements provides useful taxonomic information. We have sequenced Alu repeats from tarsier and compared them with those from strepsirhine prosimians (lemurs, sifaka, and galago) and the human genome. Tarsier elements cluster with Alu subfamilies from the human lineage. The oldest subfamily in tarsier and the most abundant human subfamilies share an RNA secondary structure motif which is absent both in the earliest dimeric Alu Jo and in the strepsirhine elements. These findings are consistent with the view that tarsiers form a sister clade with anthropoides rather than with other prosimians. Alu repeats in tarsier genome are relatively old, which indicates a dramatic slowdown or even an arrest of these elements' amplification about 20 Myr ago.     

Alu insertion profiling: array-based methods to detect Alu insertions in the human genome

The analysis of the genetic variability associated to Alu sequences was hampered by the absence of genome-wide methodologies able to efficiently detect new polymorphisms/mutations among these repetitive elements. Here we describe two Alu insertion profiling (AIP) methods based on the hybridization of Alu-flanking genomic fragments on tiling microarrays. Protocols are designed to preferentially detect active Alu subfamilies. We tested AIP methods by analyzing chromosomes 1 and 6 in two genomic samples. In genomic regions covered by array-features, with a sensitivity of 2% (AIP1) -4% (AIP2) and 5% (AIP1) -8% (AIP2) for the old J and S Alu lineages respectively, we obtained a sensitivity of 67% (AIP1) -90% (AIP2) for the young Ya subfamily. Among the loci showing sample-to-sample differences, 5 (AIP1) -8 (AIP2) were associated to known Alu polymorphisms. Moreover, we were able to confirm by PCR and DNA sequencing 4 new intragenic Alu elements, polymorphic in 10 additional individuals.     

Association of hY4 pseudogenes with Alu repeats and abundance of hY RNA-like sequences in the human genome.

Three loci having homology with the small human cytoplasmic RNA, hY4, were isolated from human genomic DNA libraries and sequenced. Each sequence contains dispersed mismatches as compared with hY4 RNA, is followed by an A-rich or A + T-rich sequence, and is bordered by direct repeats. Each of these loci, therefore, appears to constitute a small RNA class-III pseudogene. Surprisingly, two of the three loci are associated with Alu repeats. In the hY4.B7 locus, the hY4 sequence has integrated into the tail of an Alu element and in the hY4.F2 locus, an Alu sequence has inserted into the hY4 tail, confirming that A-rich tracts are preferential targets for retroposition. In addition, Southern blots with probes for each of the four hY RNAs indicate that hY RNA-like sequences are abundant in the human genome.     

Widespread occurrence of power-law distributions in inter-repeat distances shaped by genome dynamics

Repetitive DNA sequences derived from transposable elements (TE) are distributed in a non-random way, co-clustering with other classes of repeat elements, genes and other genomic components. In a previous work we reported power-law-like size distributions (linearity in log-log scale) in the spatial arrangement of Alu and LINE1 elements in the human genome. Here we investigate the large-scale features of the spatial arrangement of all principal classes of TEs in 14 genomes from phylogenetically distant organisms by studying the size distribution of inter-repeat distances. Power-law-like size distributions are found to be widespread, extending up to several orders of magnitude. In order to understand the emergence of this distributional pattern, we introduce an evolutionary scenario, which includes (i) Insertions of DNA segments (e.g., more recent repeats) into the considered sequence and (ii) Eliminations of members of the studied TE family. In the proposed model we also incorporate the potential for transposition events (characteristic of the DNA transposons' life-cycle) and segmental duplications. Simulations reproduce the main features of the observed size distributions. Furthermore, we investigate the effects of various genomic features on the presence and extent of power-law size distributions including TE class and age, mode of parental TE transmission, GC content, deletion and recombination rates in the studied genomic region, etc. Our observations corroborate the hypothesis that insertions of genomic material and eliminations of repeats are at the basis of power-laws in inter-repeat distances. The existence of these power-laws could facilitate the formation of the recently proposed "fractal globule" for the confined chromatin organization.