Alu recombination-mediated structural deletions in the chimpanzee genome

With more than 1.2 million copies, Alu elements are one of the most important sources of structural variation in primate genomes. Here, we compare the chimpanzee and human genomes to determine the extent of Alu recombination-mediated deletion (ARMD) in the chimpanzee genome since the divergence of the chimpanzee and human lineages (~6 million y ago). Combining computational data analysis and experimental verification, we have identified 663 chimpanzee lineage-specific deletions (involving a total of ~771 kb of genomic sequence) attributable to this process. The ARMD events essentially counteract the genomic expansion caused by chimpanzee-specific Alu inserts. The RefSeq databases indicate that 13 exons in six genes, annotated as either demonstrably or putatively functional in the human genome, and 299 intronic regions have been deleted through ARMDs in the chimpanzee lineage. Therefore, our data suggest that this process may contribute to the genomic and phenotypic diversity between chimpanzees and humans. In addition, we found four independent ARMD events at orthologous loci in the gorilla or orangutan genomes. This suggests that human orthologs of loci at which ARMD events have already occurred in other nonhuman primate genomes may be “at-risk” motifs for future deletions, which may subsequently contribute to human lineage-specific genetic rearrangements and disorders.     

High-resolution copy-number variation map reflects human olfactory receptor diversity and evolution

Olfactory receptors (ORs), which are involved in odorant recognition, form the largest mammalian protein superfamily. The genomic content of OR genes is considerably reduced in humans, as reflected by the relatively small repertoire size and the high fraction ( approximately 55%) of human pseudogenes. Since several recent low-resolution surveys suggested that OR genomic loci are frequently affected by copy-number variants (CNVs), we hypothesized that CNVs may play an important role in the evolution of the human olfactory repertoire. We used high-resolution oligonucleotide tiling microarrays to detect CNVs across 851 OR gene and pseudogene loci. Examining genomic DNA from 25 individuals with ancestry from three populations, we identified 93 OR gene loci and 151 pseudogene loci affected by CNVs, generating a mosaic of OR dosages across persons. Our data suggest that approximately 50% of the CNVs involve more than one OR, with the largest CNV spanning 11 loci. In contrast to earlier reports, we observe that CNVs are more frequent among OR pseudogenes than among intact genes, presumably due to both selective constraints and CNV formation biases. Furthermore, our results show an enrichment of CNVs among ORs with a close human paralog or lacking a one-to-one ortholog in chimpanzee. Interestingly, among the latter we observed an enrichment in CNV losses over gains, a finding potentially related to the known diminution of the human OR repertoire. Quantitative PCR experiments performed for 122 sampled ORs agreed well with the microarray results and uncovered 23 additional CNVs. Importantly, these experiments allowed us to uncover nine common deletion alleles that affect 15 OR genes and five pseudogenes. Comparison to the chimpanzee reference genome revealed that all of the deletion alleles are human derived, therefore indicating a profound effect of human-specific deletions on the individual OR gene content. Furthermore, these deletion alleles may be used in future genetic association studies of olfactory inter-individual differences.     

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.     

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.     

Rhesus macaque class I duplicon structures, organization, and evolution within the alpha block of the major histocompatibility complex

The alpha block of the human and chimpanzee major histocompatibility complex (MHC) class I genomic region contains 10 to 11 duplicated MHC class I genes, including the HLA/Patr-A, -G, and -F genes. In comparison, the alpha block of the rhesus macaque (Macaca mulatta, Mamu) has an additional 20 MHC class I genes within this orthologous region. The present study describes the identification and analysis of the duplicated segmental genomic structures (duplicons) and genomic markers within the alpha block of the rhesus macaque and their use to reconstruct the duplication history of the genes within this region. A variety of MHC class I genes, pseudogenes, transposons, and retrotransposons, such as Alu and ERV16, were used to categorize the 28 duplicons into four distinct structural categories. The phylogenetic relationship of MHC class I genes, Alu, and LTR16B sequences within the duplicons was examined by use of the Neighbor-Joining (NJ) method. Two single-duplicon tandem duplications, two polyduplicon tandem duplications with an accompanying inversion product per duplication, eight polyduplicon tandem duplications steps, 12 deletions, and at least two recombinations were reconstructed to explain the highly complex organization and evolution of the 28 duplicons (nine inversions) within the Mamu alpha block. On the basis of the phylogenetic evidence and the reconstructed tandem duplication history of the 28 duplicons, the Mamu/Patr/HLA-F ortholog was the first MHC class I gene to have been fixed without further duplication within the alpha block of primates. Assuming that the rhesus macaque and the chimpanzee/human lineages had started with the same number of MHC class I duplicons at the time of their divergence approximately 24 to 31 MYA, then the number of genes within the alpha block have been duplicated at an approximately three times greater rate in the rhesus macaque than in either the human or chimpanzee.     

The Contribution of Alu Elements to Mutagenic DNA Double-Strand Break Repair

Alu elements make up the largest family of human mobile elements, numbering 1.1 million copies and comprising 11% of the human genome. As a consequence of evolution and genetic drift, Alu elements of various sequence divergence exist throughout the human genome. Alu/Alu recombination has been shown to cause approximately 0.5% of new human genetic diseases and contribute to extensive genomic structural variation. To begin understanding the molecular mechanisms leading to these rearrangements in mammalian cells, we constructed Alu/Alu recombination reporter cell lines containing Alu elements ranging in sequence divergence from 0%-30% that allow detection of both Alu/Alu recombination and large non-homologous end joining (NHEJ) deletions that range from 1.0 to 1.9 kb in size. Introduction of as little as 0.7% sequence divergence between Alu elements resulted in a significant reduction in recombination, which indicates even small degrees of sequence divergence reduce the efficiency of homology-directed DNA double-strand break (DSB) repair. Further reduction in recombination was observed in a sequence divergence-dependent manner for diverged Alu/Alu recombination constructs with up to 10% sequence divergence. With greater levels of sequence divergence (15%-30%), we observed a significant increase in DSB repair due to a shift from Alu/Alu recombination to variable-length NHEJ which removes sequence between the two Alu elements. This increase in NHEJ deletions depends on the presence of Alu sequence homeology (similar but not identical sequences). Analysis of recombination products revealed that Alu/Alu recombination junctions occur more frequently in the first 100 bp of the Alu element within our reporter assay, just as they do in genomic Alu/Alu recombination events. This is the first extensive study characterizing the influence of Alu element sequence divergence on DNA repair, which will inform predictions regarding the effect of Alu element sequence divergence on both the rate and nature of DNA repair events.     

A young Alu subfamily amplified independently in human and African great apes lineages.

A variety of Alu subfamilies amplified in primate genomes at different evolutionary time periods. Alu Sb2 belongs to a group of young subfamilies with a characteristic two-nucleotide deletion at positions 65/66. It consists of repeats having a 7-nucleotide duplication of a sequence segment involving positions 246 through 252. The presence of Sb2 inserts was examined in five genomic loci in 120 human DNA samples as well as in DNAs of higher primates. The lack of the insertional polymorphism seen at four human loci and the absence of orthologous inserts in apes indicated that the examined repeats retroposed early in the human lineage, but following the divergence of great apes. On the other hand, similar analysis of the fifth locus (butyrylcholinesterase gene) suggested contemporary retropositional activity of this subfamily. By a semi-quantitative PCR, using a primer pair specific for Sb2 repeats, we estimated their copy number at about 1500 per human haploid genome; the corresponding numbers in chimpanzee and gorilla were two orders of magnitude lower, while in orangutan and gibbon the presence of Sb2 Alu was hardly detectable. Sequence analysis of PCR-amplified Sb2 repeats from human and African great apes is consistent with the model in which the founding of Sb2 subfamily variants occurred independently in chimpanzee, gorilla and human lineages.     

Evolutionary impact of human Alu repetitive elements

Early studies of human Alu retrotransposons focused on their origin, evolution and biological properties, but current focus is shifting toward the effect of Alu elements on evolution of the human genome. Recent analyses indicate that numerous factors have affected the chromosomal distribution of Alu elements over time, including male-driven insertions, deletions and rapid CpG mutations after their retrotransposition. Unequal crossing over between Alu elements can lead to local mutations or to large segmental duplications responsible for genetic diseases and long-term evolutionary changes. Alu elements can also affect human (primate) evolution by introducing alternative splice sites in existing genes. Studying the Alu family in a human genomic context is likely to have general significance for our understanding of the evolutionary impact of other repetitive elements in diverse eukaryotic genomes.     

ERVK9, transposons and the evolution of MHC class I duplicons within the alpha-block of the human and chimpanzee

The genomic sequences within the alpha-block (approximately 288-310 kb) of the human and chimpanzee MHC class I region contains ten MHC class I genes and three MIC gene fragments grouped together within alternating duplicated genomic segments or duplicons. In this study, the chimpanzee and human genomic sequences were analyzed in order to determine whether the remnants of the ERVK9 and other retrotransposon sequences are useful genomic markers for reconstructing the evolutionary history of the duplicated MHC gene families within the alpha-block. A variety of genes, pseudogenes, autologous DNA transposons and retrotransposons such as Alu and ERVK9 were used to categorize the ten duplicons into four distinct structural groups. The phylogenetic relationship of the ten duplicons was examined by using the neighbour joining method to analyze transposon sequence topologies of selected Alu members, LTR16B and Charlie9. On the basis of these structural groups and the phylogeny of the duplicated transposon sequences, a duplication model was reconstructed involving four multipartite tandem duplication steps to explain the organization and evolution of the ten duplicons within the alpha-block of the chimpanzee and human. The phylogenetic analysis and inferred duplication history suggests that the Patr/HLA-F was the first MHC class I gene to have been fixed and not required as a precursor for further duplication within the alpha-block of the ancestral species.     

Gross deletions involving IGHM, BTK, or Artemis: a model for genomic lesions mediated by transposable elements

Most genetic disruptions underlying human disease are microlesions, whereas gross lesions are rare with gross deletions being most frequently found (6%). Similar observations have been made in primary immunodeficiency genes, such as BTK, but for unknown reasons the IGHM and DCLRE1C (Artemis) gene defects frequently represent gross deletions ( approximately 60%). We characterized the gross deletion breakpoints in IGHM-, BTK-, and Artemis-deficient patients. The IGHM deletion breakpoints did not show involvement of recombination signal sequences or immunoglobulin switch regions. Instead, five IGHM, eight BTK, and five unique Artemis breakpoints were located in or near sequences derived from transposable elements (TE). The breakpoints of four out of five disrupted Artemis alleles were located in highly homologous regions, similar to Ig subclass deficiencies and Vh deletion polymorphisms. Nevertheless, these observations suggest a role for TEs in mediating gross deletions. The identified gross deletion breakpoints were mostly located in TE subclasses that were specifically overrepresented in the involved gene as compared to the average in the human genome. This concerned both long (LINE1) and short (Alu, MIR) interspersed elements, as well as LTR retrotransposons (ERV). Furthermore, a high total TE content (>40%) was associated with an increased frequency of gross deletions. Both findings were further investigated and confirmed in a total set of 20 genes disrupted in human disease. Thus, to our knowledge for the first time, we provide evidence that a high TE content, irrespective of the type of element, results in the increased incidence of gross deletions as gene disruption underlying human disease.     

Small deletion variants have stable breakpoints commonly associated with alu elements

Copy number variants (CNVs) contribute significantly to human genomic variation, with over 5000 loci reported, covering more than 18% of the euchromatic human genome. Little is known, however, about the origin and stability of variants of different size and complexity. We investigated the breakpoints of 20 small, common deletions, representing a subset of those originally identified by array CGH, using Agilent microarrays, in 50 healthy French Caucasian subjects. By sequencing PCR products amplified using primers designed to span the deleted regions, we determined the exact size and genomic position of the deletions in all affected samples. For each deletion studied, all individuals carrying the deletion share identical upstream and downstream breakpoints at the sequence level, suggesting that the deletion event occurred just once and later became common in the population. This is supported by linkage disequilibrium (LD) analysis, which has revealed that most of the deletions studied are in moderate to strong LD with surrounding SNPs, and have conserved long-range haplotypes. Analysis of the sequences flanking the deletion breakpoints revealed an enrichment of microhomology at the breakpoint junctions. More significantly, we found an enrichment of Alu repeat elements, the overwhelming majority of which intersected deletion breakpoints at their poly-A tails. We found no enrichment of LINE elements or segmental duplications, in contrast to other reports. Sequence analysis revealed enrichment of a conserved motif in the sequences surrounding the deletion breakpoints, although whether this motif has any mechanistic role in the formation of some deletions has yet to be determined. Considered together with existing information on more complex inherited variant regions, and reports of de novo variants associated with autism, these data support the presence of different subgroups of CNV in the genome which may have originated through different mechanisms.     

Comparative genomic analysis of human and chimpanzee proteases

Proteolytic enzymes are implicated in multiple physiological and pathological processes. The availability of the sequence of the chimpanzee genome has allowed us to determine that the chimpanzee degradome-the repertoire of protease genes from this organism-is composed of at least 559 protease and protease-like genes and is virtually identical to that of human, containing 561 genes. Despite the high degree of conservation between both genomes, we have identified important differences that vary from deletion of whole genes to small insertion/deletion events or single nucleotide changes that lead to the specific gene inactivation in one species, mostly affecting immune system genes. For example, the genes encoding PRSS33/EOS, a macrophage serine protease conserved in most mammals, and GGTLA1 are absent in chimpanzee, while the gene for metalloprotease MMP23A, located in chromosome 1p36, has been specifically duplicated in the human genome together with its neighbor gene CDC2L1. Other differences arise from single nucleotide changes in protease genes, such as NAPSB and CASP12, resulting in the presence of functional genes in chimpanzee and pseudogenes in human. Finally, we have confirmed that the Trypanosoma lytic factor HPR is inactive in chimpanzee, likely contributing to the susceptibility of chimpanzees to T. brucei infection. This study provides the first analysis of the chimpanzee degradome and might contribute to the understanding of the molecular bases underlying variations in host defense mechanisms between human and chimpanzee.     

Chromosomal Inversions between Human and Chimpanzee Lineages Caused by Retrotransposons

The long interspersed element-1 (LINE-1 or L1) and Alu elements are the most abundant mobile elements comprising 21% and 11% of the human genome, respectively. Since the divergence of human and chimpanzee lineages, these elements have vigorously created chromosomal rearrangements causing genomic difference between humans and chimpanzees by either increasing or decreasing the size of genome. Here, we report an exotic mechanism, retrotransposon recombination-mediated inversion (RRMI), that usually does not alter the amount of genomic material present. Through the comparison of the human and chimpanzee draft genome sequences, we identified 252 inversions whose respective inversion junctions can clearly be characterized. Our results suggest that L1 and Alu elements cause chromosomal inversions by either forming a secondary structure or providing a fragile site for double-strand breaks. The detailed analysis of the inversion breakpoints showed that L1 and Alu elements are responsible for at least 44% of the 252 inversion loci between human and chimpanzee lineages, including 49 RRMI loci. Among them, three RRMI loci inverted exonic regions in known genes, which implicates this mechanism in generating the genomic and phenotypic differences between human and chimpanzee lineages. This study is the first comprehensive analysis of mobile element bases inversion breakpoints between human and chimpanzee lineages, and highlights their role in primate genome evolution.     

Alpha-galactosidase A gene rearrangements causing Fabry disease. Identification of short direct repeats at breakpoints in an Alu-rich gene

Fabry disease, an inborn error of glycosphingolipid catabolism, results from mutations in the X-linked gene encoding the lysosomal enzyme, alpha-galactosidase A (EC 3.2.1.22). Six alpha-galactosidase A gene rearrangements that cause Fabry disease were investigated to assess the role of Alu repetitive elements and short direct and/or inverted repeats in the generation of these germinal mutations. The breakpoints of five partial gene deletions and one partial gene duplication were determined by either cloning and sequencing the mutant gene from an affected hemizygote, or by polymerase chain reaction amplifying and sequencing the genomic region containing the novel junction. Although the alpha-galactosidase A gene contains 12 Alu repetitive elements (representing approximately 30% of the 12-kilobase (kb) gene or approximately 1 Alu/1.0 kb), only one deletion resulted from an Alu-Alu recombination. The remaining five rearrangements involved illegitimate recombinational events between short direct repeats of 2 to 6 base pairs (bp) at the deletion or duplication breakpoints. Of these rearrangements, one had a 3' short direct repeat within an Alu element, while another was unusual having two deletions of 1.7 kb and 14 bp separated by a 151-bp inverted sequence. These findings suggested that slipped mispairing or intrachromosomal exchanges involving short direct repeats were responsible for the generation of most of these gene rearrangements. There were no inverted repeat sequences or alternating purine-pyrimidine regions which may have predisposed the gene to these rearrangements. Intriguingly, the tetranucleotide CCAG and the trinucleotide CAG (or their respective complements, CTGG and CTG) occurred within or adjacent to the direct repeats at the 5' breakpoints in three and four of the five alpha-galactosidase A gene rearrangements, respectively, suggesting a possible functional role in these illegitimate recombinational events. These studies indicate that short direct repeats are important in the formation of gene rearrangements, even in human genes like alpha-galactosidase A that are rich in Alu repetitive elements.