Inviting instability: Transposable elements, double-strand breaks, and the maintenance of genome integrity.

The ubiquity of mobile elements in mammalian genomes poses considerable challenges for the maintenance of genome integrity. The predisposition of mobile elements towards participation in genomic rearrangements is largely a consequence of their interspersed homologous nature. As tracts of nonallelic sequence homology, they have the potential to interact in a disruptive manner during both meiotic recombination and DNA repair processes, resulting in genomic alterations ranging from deletions and duplications to large-scale chromosomal rearrangements. Although the deleterious effects of transposable element (TE) insertion events have been extensively documented, it is arguably through post-insertion genomic instability that they pose the greatest hazard to their host genomes. Despite the periodic generation of important evolutionary innovations, genomic alterations involving TE sequences are far more frequently neutral or deleterious in nature. The potentially negative consequences of this instability are perhaps best illustrated by the >25 human genetic diseases that are attributable to TE-mediated rearrangements. Some of these rearrangements, such as those involving the MLL locus in leukemia and the LDL receptor in familial hypercholesterolemia, represent recurrent mutations that have independently arisen multiple times in human populations. While TE-instability has been a potent force in shaping eukaryotic genomes and a significant source of genetic disease, much concerning the mechanisms governing the frequency and variety of these events remains to be clarified. Here we survey the current state of knowledge regarding the mechanisms underlying mobile element-based genetic instability in mammals. Compared to simpler eukaryotic systems, mammalian cells appear to have several modifications to their DNA-repair ensemble that allow them to better cope with the large amount of interspersed homology that has been generated by TEs. In addition to the disruptive potential of nonallelic sequence homology, we also consider recent evidence suggesting that the endonuclease products of TEs may also play a key role in instigating mammalian genomic instability.

Human retroelements may introduce intragenic polyadenylation signals.

In the human genome, the insertion of LINE-1 and Alu elements can affect genes by sequence disruption, and by the introduction of elements that modulate the gene's expression. One of the modulating sequences retroelements may contribute is the canonical polyadenylation signal (pA), AATAAA. L1 elements include these within their own sequence and AATAAA sequences are commonly created in the A-rich tails of both SINEs and LINEs. Computational analysis of 34 genes randomly retrieved from the human genome draft sequence reveals an orientation bias, reflected as a lower number of L1s and Alus containing the pA in the same orientation as the gene. Experimental studies of Alu-based pA sequences when placed in pol II or pol III transcripts suggest that the signal is very weak, or often not used at all. Because the pA signal is highly affected by the surrounding sequence, it is likely that the Alu constructs evaluated did not provide the required recognition signals to the polyadenylation machinery. Although the effect of pA signals contributed by Alus is individually weak, the observed reduction of "sense" oriented pA-containing L1 and Alu elements within genes reflects that even a modest influence causes a change in evolutionary pressure, sufficient to create the biased distribution.

Worldwide genetic variation at the 3'-UTR region of the LDLR gene: possible influence of natural selection.

The low density lipoprotein receptor gene (LDLR) contains many Alu insertions, and is especially Alu-rich at its 3'-untranslated region (3'-UTR). Previous studies suggested that the LDLR 3'-UTR could regulate gene expression by the stabilization of its mRNA. Given the faster Alu evolutionary rate, and wondering about its consequences in a possibly regulatory locus, we have studied approximately 800 bp of 222 chromosomes from individuals of African, Asian, Caucasian and Amerind ancestry, to better understand the evolution of the worldwide genetic diversity at this locus. Twenty-one polymorphic sites, distributed in 15 haplotypes, were found. High genetic diversity was observed, concentrated in one Alu insertion (Alu U), which also shows a fast evolutionary rate. Genetic diversity is similar in all populations except Amerinds, suggesting a bottleneck during the peopling of the American continent. Three haplotype clusters (A, B, C) are distinguished, cluster A being the most recently formed (approximately 500,000 years ago). No clear geographic structure emerges from the haplotype network, the global F(st) (0.079) being lower than the average for the human genome. When ancestral population growth is taken into account, neutrality statistics are higher than expected, possibly suggesting the action of balancing selection worldwide.

Alu insertion polymorphisms for the study of human genomic diversity.

Genomic database mining has been a very useful aid in the identification and retrieval of recently integrated Alu elements from the human genome. We analyzed Alu elements retrieved from the GenBank database and identified two new Alu subfamilies, Alu Yb9 and Alu Yc2, and further characterized Yc1 subfamily members. Some members of each of the three subfamilies have inserted in the human genome so recently that about a one-third of the analyzed elements are polymorphic for the presence/absence of the Alu repeat in diverse human populations. These newly identified Alu insertion polymorphisms will serve as identical-by-descent genetic markers for the study of human evolution and forensics. Three previously classified Alu Y elements linked with disease belong to the Yc1 subfamily, supporting the retroposition potential of this subfamily and demonstrating that the Alu Y subfamily currently has a very low amplification rate in the human genome.

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.

Phylogenetic analysis of the Friedreich ataxia GAA trinucleotide repeat.

Friedreich ataxia is an autosomal recessive neurodegenerative disorder associated with a GAA repeat expansion in the first intron of the gene (FRDA) encoding a novel, highly conserved, 210 amino acid protein known as frataxin. Normal variation in repeat size was determined by analysis of more than 600 DNA samples from seven human populations. This analysis showed that the most frequent allele had nine GAA repeats, and no alleles with fewer than five GAA repeats were found. The European and Syrian populations had the highest percentage of alleles with 10 or more GAA repeats, while the Papua New Guinea population did not have any alleles carrying more than 10 GAA repeats. The distributions of repeat sizes in the European, Syrian, and African American populations were significantly different from those in the Asian and Papua New Guinea populations (p < 0.001). The GAA repeat size was also determined in five nonhuman primates. Samples from 10 chimpanzees, 3 orangutans, 1 gorilla, 1 rhesus macaque, 1 mangabey, and 1 tamarin were analyzed. Among those primates belonging to the Pongidae family, the chimpanzees were found to carry three or four GAA repeats, the orangutans had four or five GAA repeats, and the gorilla carried three GAA repeats. In primates belonging to the Cercopithecidae family, three GAA repeats were found in the mangabey and two in the rhesus macaque. However, an AluY subfamily member inserted in the poly(A) tract preceding the GAA repeat region in the rhesus macaque, making the amplified sequence approximately 300 bp longer. The GAA repeat was also found in the tamarin, suggesting that it arose at least 40 million years ago and remained relatively small throughout the majority of primate evolution, with a punctuated expansion in the human genome.

Cis-acting influences on Alu RNA levels.

The human short interspersed repeated element (SINE), Alu, amplifies through a poorly understood RNA-mediated mechanism, termed retroposition. There are over one million copies of Alu per haploid human genome. The copies show some internal variations in sequence and are very heterogeneous in chromosomal environment. However, very few Alu elements actively amplify. The amplification rate has decreased greatly in the last 40 million years. Factors influencing Alu transcription would directly affect an element's retroposition capability. Therefore, we evaluated several features that might influence expression from individual Alu elements. The influence of various internal sequence variations and 3' unique flanks on full-length Alu RNA steady-state levels was determined. Alu subfamily diagnostic mutations do not significantly alter the amount of Alu RNA observed. However, sequences containing random mutations throughout the right half of selected genomic Alu elements altered Alu RNA steady-state levels in cultured cells. In addition, sequence variations at the 3' unique end of the transcript also significantly altered the Alu RNA levels. In general, sequence mutations and 3' end sequences contribute to Alu RNA levels, suggesting that the master Alu element(s) have a multitude of individual differences that collectively gives them a selective advantage over other Alu elements.

Potential gene conversion and source genes for recently integrated Alu elements.

Alu elements comprise >10% of the human genome. We have used a computational biology approach to analyze the human genomic DNA sequence databases to determine the impact of gene conversion on the sequence diversity of recently integrated Alu elements and to identify Alu elements that were potentially retroposition competent. We analyzed 269 Alu Ya5 elements and identified 23 members of a new Alu subfamily termed Ya5a2 with an estimated copy number of 35 members, including the de novo Alu insertion in the NF1 gene. Our analysis of Alu elements containing one to four (Ya1-Ya4) of the Ya5 subfamily-specific mutations suggests that gene conversion contributed as much as 10%-20% of the variation between recently integrated Alu elements. In addition, analysis of the middle A-rich region of the different Alu Ya5 members indicates a tendency toward expansion of this region and subsequent generation of simple sequence repeats. Mining the databases for putative retroposition-competent elements that share 100% nucleotide identity to the previously reported de novo Alu insertions linked to human diseases resulted in the retrieval of 13 exact matches to the NF1 Alu repeat, three to the Alu element in BRCA2, and one to the Alu element in FGFR2 (Apert syndrome). Transient transfections of the potential source gene for the Apert's Alu with its endogenous flanking genomic sequences demonstrated the transcriptional and presumptive transpositional competency of the element.

Upstream flanking sequences and transcription of SINEs.

SINEs, short interspersed repeated DNA elements, undergo amplification through retroposition and subsequent integration into a new location in the genome. Each new SINE insertion will be located in a new chromosomal environment, with different flanking sequences. Modulation of transcription by different flanking sequences may play an important role in determining which SINE elements are preferentially active in a genome. We evaluated the ability of upstream flanking sequences to regulate the transcription of three different SINEs (Alu, B2 and ID) by constructing chimeric constructs with known 5' flanking sequences of RNA polymerase III-transcribed genes. Upstream sequences from the 7SL RNA gene, U6 RNA gene, vault RNA gene, and BC1 gene increase transcription of Alu, B2 and BC1 in transient transfections of NIH3T3, HeLa, Neuro2a and C6 glioma cell lines. The 7SL sequence proved most efficient in increasing SINE transcription. The 7SL upstream fused to the BC1 RNA gene (an ID element) was used to create a transgenic mouse line. In contrast to the tissue-specific endogenous BC1 transcription, BC1 transgene transcripts were detected in all tissues tested. However, expression was much higher in those tissues that express the endogenous gene, demonstrating both transcriptional and post-transcriptional regulation. The BC1 RNA was detected in a similar ribonucleoprotein complex in the different tissues.

Alu repeats and human disease.

Alu elements have amplified in primate genomes through a RNA-dependent mechanism, termed retroposition, and have reached a copy number in excess of 500,000 copies per human genome. These elements have been proposed to have a number of functions in the human genome, and have certainly had a major impact on genomic architecture. Alu elements continue to amplify at a rate of about one insertion every 200 new births. We have found 16 examples of diseases caused by the insertion of Alu elements, suggesting that they may contribute to about 0.1% of human genetic disorders by this mechanism. The large number of Alu elements within primate genomes also provides abundant opportunities for unequal homologous recombination events. These events often occur intrachromosomally, resulting in deletion or duplication of exons in a gene, but they also can occur interchromosomally, causing more complex chromosomal abnormalities. We have found 33 cases of germ-line genetic diseases and 16 cases of cancer caused by unequal homologous recombination between Alu repeats. We estimate that this mode of mutagenesis accounts for another 0.3% of human genetic diseases. Between these different mechanisms, Alu elements have not only contributed a great deal to the evolution of the genome but also continue to contribute to a significant portion of human genetic diseases.

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.

Simultaneous analysis of multiple gene expression patterns as a function of development, injury or senescence.

Concurrent changes in expression of eight genes were examined following cryogenic rat brain injury. Cortical RNA levels were catalogued at time 0, and at 1 h and 1 week following injury. The genes include thymidine kinase (TK), c-fos, renin, myelin basic protein (MBP), proteolipid protein (PLP), glial fibrillary acidic protein (GFAP), insulin-like growth factor-1 (IGF-1), and somatostatin. All demonstrate increased expression following injury. Renin and c-fos exhibit detectable changes as early as 1 h post-injury.

Identification and mutation analysis of a cochlear-expressed, zinc finger protein gene at the DFNB7/11 and dn hearing-loss loci on human chromosome 9q and mouse chromosome 19.

The DFNB7/11 locus for autosomal recessive non-syndromic hearing loss (ARNSHL) has been mapped to an approx. 1.5 Mb interval on human chromosome 9q13-q21. We have determined the cDNA sequence and genomic structure of a novel cochlear-expressed gene, ZNF216, that maps to the DFNB7/11 interval. The mouse orthologue of this gene maps to the murine dn (deafness) locus on mouse chromosome 19. The ZNF216 gene is highly conserved between human and mouse, and contains two regions that show homology to the putative zinc linger domains of other proteins. To determine it mutations in ZNF216 might be the cause of hearing loss at the DFNB7/11 locus, we screened the coding region of this gene in DFNB7/11 families by direct sequencing. No potential disease-causing mutations were found. In addition, Northern blot analysis showed no difference in ZNF216 transcript size or abundance between dn and control mice. These data Suggest that the ZNF216 gene is unlikely to be responsible for hearing loss at the DFNB7/11 and dn loci.

The mouse deafness locus (dn) is associated with an inversion on chromosome 19.

Recombination data for the mouse deafness locus (dn) on chromosome 19 are consistent with the presence of an inversion for which one of the breakpoints is between D19Mit14 and D19Mit96, a distance of less than 226 kb. Fluorescence in situ hybridization studies using a bacterial artificial chromosome on interphase (G1) nuclei provide additional support for the presence of an inversion. The dn gene is probably the orthologue of the human DFNB7/DFNB11 gene on chromosome 9.

Novel variant of the P2X2 ATP receptor from the guinea pig organ of Corti.

ATP functions as a neurotransmitter and a neuromodulator in various tissues by acting on metabotropic (P2Y) and ionotropic (P2X) receptors. Evidence suggests that ATP activates P2X receptors on several cell types in the organ of Corti of guinea pig including outer hair cells (OHCs), Deiters' cells, Hensen's cells, pillar cells and inner hair cells (IHCs). Determining the sequence and structure of P2X receptors in guinea pig organ of Corti is important for understanding the function of ATP in the cochlea. We screened a guinea pig organ of Corti cDNA library for P2X2 ATP receptors using rat P2X2 cDNA as a probe. We sequenced three P2X2 variants which were found to be abundant in this library. One is a novel P2X2 isoform (P2X2-3) created by a retained intron coding for an additional 27 amino acids (81 bp) in the putative extracellular domain. We have also sequenced a variant (P2X2-2) that lacks both the 81-bp sequence and a 192-bp sequence in the 3' intracellular domain. A third variant (P2X2-1) contains the intracellular 192-bp sequence but not the extracellular 81-bp sequence found in P2X2-3. The multiple transcripts arise from alternative intron and exon splicing events. In situ hybridization with a probe common to the three variants localized P2X2 to many of the cells of the organ of Corti.

Temporal changes in gene expression following cryogenic rat brain injury.

Expression of 18 genes was examined at 8 different time points between 1 h and 28 days following cryogenic rat brain injury. The genes include thymidine kinase (TK), p53 tumor suppressor, c-fos, renin, myelin basic protein (MBP), proteolipid protein (PLP), transferrin, transferrin receptor, platelet-derived growth factor A (PDGF A), platelet-derived growth factor B (PDGF B), platelet-derived growth factor receptor alpha (PDGF alpha receptor), platelet-derived growth factor receptor beta (PDGF beta receptor), glial fibrillary acidic protein (GFAP), transforming growth factor-beta 1 (TGF-beta 1), basic fibroblast growth factor (bFGF), fibroblast growth factor receptor-1 (FGF-R1), insulin-like growth factor-1 (IGF-1), and somatostatin. Time courses of gene expression were determined for RNAs derived from hippocampus and cortex. Genes were divided into categories based upon those in which statistically significant changes in expression were first observed at or before 24 h (early genes) and those in which changes were first observed at or after 72 h (late genes). In the present model, many genes demonstrate elevated RNA levels in the cortex prior to hippocampus, following injury. RNAs transcribed from late genes tend to be elevated concurrently in cortex and hippocampus.

Assembly of a high-resolution map of the Acadian Usher syndrome region and localization of the nuclear EF-hand acidic gene.

Usher syndrome type 1C (USH1C) occurs in a small population of Acadian descendants from southwestern Louisiana. Linkage and linkage disequilibrium analyses localize USH1C to chromosome 11p between markers D11S1397 and D11S1888, an interval of less than 680 kb. Here, we refine the USH1C linkage to a region less than 400 kb, between genetic markers D11S1397 and D11S1890. Using 17 genetic markers from this interval, we have isolated a contiguous set of 60 bacterial artificial chromosomes (BACs) that span the USH1C critical region. Exon trapping of BAC clones from this region resulted in the recovery of an exon of the nuclear EF-hand acidic (NEFA) gene. However, DNA sequence analysis of the NEFA cDNA from lymphocytes of affected individuals provided no evidence of mutation, making structural mutations in the NEFA protein unlikely as the cellular cause of Acadian Usher syndrome.

Polymorphic Alu insertions and the Asian origin of Native American populations.

A rapid PCR-based assay was used to study the distribution of 5 polymorphic Alu insertions in 895 unrelated individuals from 30 populations, 24 from North, Central, and South America. Although a significant level of interpopulation variability was detected, the variability was less than that observed in a worldwide population survey. This is consistent with the bottleneck effect and genetic drift forces that may have acted on the migrating founder groups. The results corroborate the Asian origin of native American populations but do not support the multiple-wave migration hypothesis supposedly responsible for the tri-partite Eskaleut, Nadene, and Amerind linguistic groups. Instead, these populations exhibit three major identifiable clusters reflecting geographic distribution. Close similarity between the Chinese and native Americans suggests recent gene flow from Asia.

Recent B2 element insertions in the mouse genome.

B2 elements are a family of short interspersed repeats that have amplified within rodent genomes. Recent mobility of only two individual B2 elements has been reported to date. We identified an additional recent B2 insertion occurring within intron 4 of the murine beta-glucuronidase gene (Gus-s) of the BalbC strain of mouse by analyzing orthologous loci of a nonrandomly selected B2 element. The basis of selection for the B2 element was its high level of sequence identity to the B2 consensus. The selected B2 element was amplified by the polymerase chain reaction (PCR) using primers to the unique flanking sequences from genomic DNA of several species and laboratory strains of mice. Our results demonstrated the presence of the selected B2 element only in the genome of Mus musculus BalbC strain. Cloning and sequencing of a representative sample of the products obtained confirmed the absence of the B2 element within this intron in addition to other variations in the sequence. The detection of the B2 element only in the BalbC strain suggests that the element recently inserted within this mouse population when the initial laboratory colony was formed. Sequence comparison of the two previously identified recent B2 inserts also shows a low divergence in relation to the B2 type II consensus. The data presented confirms that recently inserted B2 elements closely match their consensus sequence, potentially allowing for their identification.

Identification and characterization of two polymorphic Ya5 Alu repeats.

Two new polymorphic Alu elements (HS2.25 and HS4.14) belonging to the young (Ya5/8) subfamily of human-specific Alu repeats have been identified. DNA sequence analysis of both Alu repeats revealed that each Alu repeat had a long 3'-oligo-dA-rich tail (41 and 52 nucleotides in length) and a low level of random mutations. HS2.25 and HS4.14 were flanked by short precise direct repeats of 8 and 14 nucleotides in length, respectively. HS2.25 was located on human chromosome 13, and HS4.14 on chromosome 1. Both Alu elements were absent from the orthologous positions within the genomes of non-human primates, and were highly polymorphic in a survey of twelve geographically diverse human groups.