Skip to main content

Alu elements contain many binding sites for transcription factors and may play a role in regulation of developmental processes

Abstract

Background

The human genome contains over one million Alu repeat elements whose distribution is not uniform. While metabolism-related genes were shown to be enriched with Alu, in structural genes Alu elements are under-represented. Such observations led researchers to suggest that Alu elements were involved in gene regulation and were selected to be present in some genes and absent from others. This hypothesis is gaining strength due to findings that indicate involvement of Alu elements in a variety of functions; for example, Alu sequences were found to contain several functional transcription factor (TF) binding sites (BSs). We performed a search for new putative BSs on Alu elements, using a database of Position Specific Score Matrices (PSSMs). We searched consensus Alu sequences as well as specific Alu elements that appear on the 5 Kbp regions upstream to the transcription start site (TSS) of about 14000 genes.

Results

We found that the upstream regions of the TSS are enriched with Alu elements, and the Alu consensus sequences contain dozens of putative BSs for TFs. Hence several TFs have Alu-associated BSs upstream of the TSS of many genes. For several TFs most of the putative BSs reside on Alu; a few of these were previously found and their association with Alu was also reported. In four cases the fact that the identified BSs resided on Alu went unnoticed, and we report this association for the first time. We found dozens of new putative BSs. Interestingly, many of the corresponding TFs are associated with early markers of development, even though the upstream regions of development-related genes are Alu-poor, compared with translational and protein biosynthesis related genes, which are Alu-rich. Finally, we found a correlation between the mouse B1 and human Alu densities within the corresponding upstream regions of orthologous genes.

Conclusion

We propose that evolution used transposable elements to insert TF binding motifs into promoter regions. We observed enrichment of biosynthesis genes with Alu-associated BSs of developmental TFs. Since development and cell proliferation (of which biosynthesis is an essential component) were proposed to be opposing processes, these TFs possibly play inhibitory roles, suppressing proliferation during differentiation.

Background

Over 90% of the human and other vertebrate genomes don't have any known functional role [1, 2], and as such were considered until recently as "Junk DNA" or genomic "dark matter". A large fraction of the "non-functional" DNA originated from mobile elements [1] such as Alu, which comprise about 10% of the nucleotides of the human genome [1] with over one million inserted copies. This abundance is somewhat of a surprise, since Alu is a non-autonomous retroelement, i.e. it doesn't encode proteins that assist its mobilization, and therefore it needs to rely on the cell's machinery for its duplication in the genome [3]. Currently Alu is believed to be a parasite of the mobilization machinery of long interspersed elements (LINE) [4].

A typical Alu is a dimer, comprised of a central A-rich region which is flanked by two similar sequence elements of about 130 bp (left and right arms). During evolution the Alu elements have spread in the human genome in several bursts at different times, which allow their classification into several families and subfamilies, on the basis of their insertions, deletions and mutations. The major families are: old AluJ, whose members are AluJb and AluJo (both spread in the genome 80–100 mya); the middle, AluS family, which includes AluSx, AluSg, AluSp AluSc and AluSq (35 –50 mya), and the youngest Alu family, AluY (25 mya) [5, 6].

One of the biggest mysteries associated with Alu concerns the non-random manner in which these elements are distributed throughout the human genome [1, 7]. The distribution of Alu was found to be highly correlated with local GC content and with the density of genes and with intron density [8]. Furthermore, the density of Alu is higher in intragenic (12.5% of the nucleotides) than in intergenic regions (9.6% of the nucleotides) [8], which is surprising if indeed Alu elements are non-functional bits of DNA. In addition, Alu elements are negatively selected in imprinted regions of primate genomes [9]. Transposable elements are underrepresented in the region surrounding the TSS of genes [10]. Another analysis revealed that almost 20% of the genes contain transposable elements in the 3' UTR [11]. Significant differences were observed also when Alu density was compared among full length gene sequences of different biological process classes (this study [12] was limited to chromosomes 21 and 22). Genes that were related to metabolism, transport and signalling were mostly Alu rich, whereas genes that had poor Alu content belonged mostly to structural proteins and to information processing and storage pathways [12]. It was also shown recently that housekeeping genes contain more Alu than tissue specific genes [13].

This non-random of distribution of Alu, in particular the high concentration near and within genes, might be the result of positive feedback; selective pressure drives the organism to use these elements for some new functions, and the functions acquired by Alu generate advantage for phenotypes that have high Alu concentration near genes. Indeed, the possible functional roles of Alu in genomic organization and gene expression has received increasing attention during the last two decades (see [14, 15] for reviews). One of the more fascinating possible consequences of integration of mobile elements near genes is acquisition of a transcriptional regulatory role by the Alu element, as a carrier of different TF binding sites. The Alu element harbours BSs of several transcription factors that were shown to bind to Alu; most of these TFs were nuclear factors, hormones, calcium nuclear factors, as well as other TFs [16–22]. The idea of Alu acting as a carrier of cis regulatory elements was suggested by Britten [23] and by Oei at el. [20] who found YY1 and SP1 BSs on young Alu family members. These experimental results imply that Alu elements have influence on the regulation of expression levels of many genes that have Alu in their promoter regions; however, such a wide ranged regulatory impact has not yet been tested experimentally.

The main aim of our study was to find new putative BSs that reside on Alu. Very recently, several groups have performed similar analysis for all transposable element classes; however, the set of PSSMs studied [24, 25] contained about half of what is presented here. In addition, our analysis is complemented by an extensive literature search of past experimental work on the Alu-associated binding sites that we found, as well as the biological functions of the corresponding TFs and target genes. We started by searching the Alu consensus sequences and found new, previously unnoticed BSs, many of which turned out to be related to developmental processes.

The regulatory role played by Alu for genes that are related to differentiation or development was already established for at least six genes [26]. In these cases Alu was shown to act as enhancer or silencer involved in regulating the expression of six developmental genes. For one of these, the T-cell marker gene CD 8α, the silencing mechanism was shown to act via binding of LYF-1, bHLH and GATA3 to the Alu situated in the last intron of CD 8α [27, 28]. These experimental results demonstrated the regulatory impact of Alu on developmental genes and during developmental stages. To explore further and enhance the possible regulatory roles of Alu we focused on the genomic region most likely to have a regulatory function – the upstream regions of genes. We scanned these regions and searched for BSs on the particular Alu sequences that were identified. Our findings show for the first time that TFs that are known to regulate the developmental processes may bind to Alu elements that are incorporated in the promoter regions of genes that need to be activated or suppressed during development (such binding has already been demonstrated for Gfi-1, PITX2 and Nkx2.5 – see below).

Results and discussion

The Alu consensus sequences are enriched with Binding Motifs

We focused on the eight major Alu subfamilies AluJo, Jb, Sx, Sz, Sc, Sq, Sp and Y [3], whose sequences were downloaded from RepBase [29]. We refer to the ensemble of sequences that bind a TF as it's Binding Motif (BM). In order to find BMs of a TF, we used its Position Specific Score Matrix (PSSM); 410 PSSMs were downloaded from two databases [30, 31]. For each PSSM we set two threshold scores, T5 and T1, selected so that the probability to find a target BM with score higher than T5 in a suitably chosen random sequence (see Methods) is less than 0.05 (0.01 for T1). Furthermore, for each BM we identified Tmax, the maximal score observed on all the above listed Alu subfamilies' consensus sequences. For 66 BMs Tmax was higher than T5 and for 25 TFs – higher than T1 (see Table 1). A good measure for the statistical significance of enrichment of Alu with BMs is the number of BMs whose score exceeds T5; for example, we found on AluSx 35 BMs with scores higher than their respective T5. The probability to find this many BMs in a random model (i.e. p-value) is 0.023 (see Methods).

Table 1 Twenty five binding motifs whose scores pass T1. Those that pass FDR of Q = 0.20 on at least one major Alu subfamily sequence are in bold. The consensus sites of the PSSMs (from the IUPAC convention [33]) is in the second column and the target sequences with the highest scores among all the major Alu consensus sequences – in the third, with nucleotides that agree with the consensus denoted by capital letters. The fourth column contains the locations of the putative target sites on the Alu sequences of the subfamily with the highest score (in case the same sequence appears on several Alu elements, we choose the one with the highest number of copies in the 5 kb upstream regions). Two Alu subsequences serve as putative target sites of several TFs (designated by*, and **). The fifth column contains the p-values (see text and methods) and the number of subfamilies on which the BSs of the third column resides is listed in column 6.

In order to control the problem of multiple testing we performed a false discovery rate (FDR) analysis at 20% FDR[32]. The 14 BMs that passed the FDR test with Q = 0.20 in AluSx are listed in Table 1, where we also compare the IUPAC consensus sequences of these BMs [33] with the target sequence that had the highest score in all the major Alu subfamilies. Most sequences with scores higher than T1 differ at not more than two nucleotides from the IUPAC consensus (see Table 1 and Additional data file 1), and there are several PSSMs, such as PITX2_Q2, LYF-1.01, PAX4.01, SIX3.01, DELTAEF.01 and NKX25.01 (Table 1, and Additional data file 1), whose IUPAC consensus sequences agree with the highest-scoring target sites on at least one of the Alu consensus sequences.

BMs tend to cluster together on specific parts of Alu consensus sequences, that also tend to be more conserved during evolution

We checked the BM's locations along the consensus sequences on both strands of the eight Alu subfamilies. In general, most of the putative BSs were found on the sense strand. It is interesting to note that Alu elements show a 55% to 45% bias towards antisense orientation [34] relative to nearby genes. It would be of considerable interest to check experimentally whether this orientational difference is reflected in functional roles of the TFs. The target sites that passed T5 cover, on average, 50% of the consensus sequence of each subfamily. On the other hand, the total length of all the binding motifs found on one Alu consensus sequence exceeds, on the average, 420 nucleotides (150% of the Alu's length). This redundancy in putative target sites is caused by the fact that several PSSMs have very similar structure, giving top scores to the same DNA sequence.

The short sequences GGTAATCCC on the 5' to 3' strand and its complementary sequence GGGATTACC serve as a putative BS of five TFs, four of which are homeodomain proteins (Table 1 and Additional data file 1, TFs and locations marked by *). This sequence is part of the consensus of all eight Alu subfamilies, and in the five middle AluS subfamilies the BS appears on both the left and right monomers. Moreover, this is precisely the sequence identified by Britten [23] as one of the most conserved subsequences on Alu, which was seen as a surprise since it did not have any role in retrotransportation or for any other function. Later on, it was found that in the germ line the Specific Alu Binding Protein (SABP) can bind to this sequence and prevent methylation [35], but since the importance of this function of is still not understood, and hence the high level of conservation is still a mystery. Another region is a putative harbour for several BMs is the middle poly A- stretch (that bridges the two Alu monomers). We found this sequence to be a putative target site for many binding motifs such as all MEF2 family members, HNF1.03, OC.2 and PAX4 (marked by ** in Table 1 and Additional data file 1). These putative target sites were also found to accumulate less mutations than the expected rate [23]. More recently both the central poly A stretch and the subsequence mentioned above [36–38] were found to be associated with A-to-I editing.

Existing evidence for TF binding to Alu (on DNA in solution)

Having identified BSs of various TFs on Alu, we searched for existing experimental evidence for binding and regulation associated with these motifs. Indeed, several TFs that are known to bind to motifs that passed our thresholds have already been shown to bind to DNA, and the fact that the binding sequences reside on Alu has been noticed (denoted by * on Table 2). We show in Table 2 four additional BMs that we found and report here for the first time; their corresponding TFs were shown to bind to the target sites identified by us, but the fact that they reside on Alu has not been previously noted. These TFs (bold in Table 2) are PITX2, Nkx2.5 Gfi-1 and LUN1. The most striking example was LUN-1; Chu et al. [39] observed (by EMSA) that it binds to 55 DNA sequences, 27 of which are palindromic; we discovered that all of these 27 reside on Alu. The situation is different for Gfi-1; there the PSSM [30] is based on 21 target sites, only five of which are part of Alu (Table 3). As to PITX2, it's PSSM [31] was based on target sites on the promoter of PLOD1, all of which reside on Alu subsequences [40]. Finally, the heart development marker Nkx2.5 is the only one whose PSSM was calculated using only off-Alu target sites – nevertheless, the target site found by us, that resides on Alu, is the one with the highest possible score.

Table 2 TFs for which the binding to the Alu-associated BM has been experimentally verified. Most passed T1 and FDR of 20%. The last two (marked by **) did not pass our promoter regions search, but were reported in the literature as bound to Alu. Four TFs (bold) are new: they were found by us but their being bound to Alu went unnoticed. The first column is the PSSM symbol from MatInspector [30] or TRANSCFAC [31] the second and third columns are the gene symbol and description of the biological process. The last two columns list those target genes of each TF that have on-Alu BSs, and the appropriate references.
Table 3 List of human target genes of (murine) Gfi-1 which contain Alu in their 5' upstream region. In these genes some of the Gfi-1 target sites are on Alu sequences; such sequences are denoted by '+' in the third column. The first column contains the target gene name; the second column is a potential binding sequence in the promoter region of the gene. The fourth column is a result of EMSA binding experiment for the sequence in column 2 [41].

Evidence (in live cells) for regulation by TFs bound to Alu

The regulatory role of single – gene BSs on Alu was demonstrated for several TFs by a transfection assay for hormone ligands [16, 22], for YY1 [20], and for Estrogen Receptors (ER) [18]. We found experimental evidence indicating that Gfi-1, PITX2 and Nkx2.5 regulate a reporter gene's expression [40–42] via BSs that reside on Alu; where the PITX2 function was validated in two different experiments [40, 42]. In addition, Alu was essential for regulation of PLOD1 by PITX2, shown by using two constructs of PLOD1 promoter: one with several Alu elements and one without any Alu. These researchers [40] found a 3-fold change in the induction of a reporter gene with constructs that contained several Alu elements with respect to the expression level of the construct from which Alu elements were excluded. The same procedure of comparing promoters with and without Alu was done again on PLOD1 [42], but this time they checked PITX2 together with Nkx2.5 in order to prove the synergic effect on PLOD1 expression, and again a big difference between the expression levels of the reporter gene was found.

ChIP analysis was done on 21 genes' promoters that contain Gfi-1 putative target sites; five of these promoters contained suspected BSs on Alu (table 3). However, only one of these five genes (JAK3) had its Gfi-1 BS exclusively on Alu [41]. Hence the positive result of the ChIP experiment on this gene shows that Gfi-1 actually binds in vivo to Alu.

Developmental and stress response TFs constitute a majority of the BSs on Alu

Most of the stress response TFs that bind to Alu are nuclear factors and hormone ligands that were validated experimentally, such as ER, RAR, LXRE and SP1 (Tables 2, 4). As was mentioned above, a subset of them are common to all Alu subfamilies while some appear only on the young Alu sequences, or on Alu elements that were mutated [20]. Another interesting group of TFs is developmental (usually referred to either as early markers or regulators of development); more than 30 members of this group passed the T5 threshold on at least one Alu consensuses sequence. The largest set of developmental TF binding sites that appear on Alu belong to hemateopoietic transcription factors [43] such as Gfi-1, LYF-1, GATA3, GATA2 and ONECUT2 (Table 4). Another set of BMs belong to TFs that are related to muscle and heart development [44], such as PITX2, MEF2, NKX2.5 and LUN1 (Tables 2, 4), and some to brain development: PAX4, NRL, OTX2 (Tables 2, 4). Note that even though the PAX6.01 PSSM did pass our thresholds (and hence appears in Table 4), it has been demonstrated [45] that PAX6 actually binds to very specific Alu sequences, that constitute only 0.2% of the genome-wide Alu elements. The PSSM that was used in our analysis (taken from MatInspector) was constructed from BMs that were identified in different experiments and its dominant binding sequences do not coincide with those that [45] identified as the Alu mediated PAX6 BMs.

Table 4 List of biological processes whose regulation involves transcription factors with putative target sites on Alu. Those transcription factors that were experimentally found to bind to Alu (see Table 2) are in bold. These reported binding events involved specific Alu sequences, and were not based on genome-wide screening.

The distribution of transposable elements in the 5' upstream regions of human genes

The distribution of transposable elements in the 5' upstream regions of human genesis not uniform. We computed the density of different transposable elements on 5000 bp long upstream sequences of about 14000 genes. In the region between 2000 bp and 5000 bp from the transcription start site (TSS), 45% of the nucleotides belong to transposable elements. This fraction is the same as the overall density in the entire genome – see Fig 1. In contrast, in the region closer to the TSS, between 0 and 2000 bp – the density of transposable elements is lower. The Alu density exhibits a more interesting behaviour; for short distances – between 0 and 600 bp – it is lower than the genomic average of about 10%. This is not surprising since the immediate proximal region to the TSS is rich with important binding sites for initiating factors and insertion of Alu and other transposable elements in this region is most probably selected against. The density increases steadily till about 2000 bp, where it reaches a level of 18% of the base pairs [Fig. 1], after which it maintains a steady level (and has to decrease to the genomic average as we move further away into the intergenic region). The Alu content in the deep intergenic spacer regions is less than 10%; hence the density of on-Alu BSs in these regions is also lower than in the upstream 5000 bp regions of genes. The presences of less BSs in these regions might imply less functional relevance. Finally, 85% of the genes in our database contained at least one Alu insertion in the 5 Kbp upstream region.

Figure 1
figure 1

The nucleotide density of different repetitive elements in the 5 Kbp upstream to the TSS regions of 13686 genes. The average genomic Alu content is 10% and the nucleotide density associated with all the repetitive elements is about 45% [1].

Non random distribution of Alu elements on promoter regions of different gene groups

We analyzed the distribution of Alu elements in the 5000 bp upstream regions of genes, treating separately groups of genes that belong to different biological processes. Protein biosynthesis genes are enriched with Alu compared to the genome-wide average value of 3.63 copies in the 5000 bp long upstream regions (Table 5). In contrast, genes that are related to various kinds of developmental processes, tissue formation (such as cell communication, central nervous system development genes and muscle development genes) contain a smaller number of Alu elements in their 5000 bp upstream region than the average value (see Table 5). Many of these differences between the Alu content of genes of a particular biological process and the Alu content of the entire genome (in the 5 Kbp upstream regions) were found (by t-test) to be highly significant (Table 5). As mentioned above, similar differences in Alu content were reported [12] when the full length gene sequences on chromosomes 21 and 22 were analyzed.

Table 5 List of biological processes and statistics of the corresponding Alu content in the 5 kbp upstream regions of the associated genes.

Studying BS content in specific (non-consensus) Alu elements in the promoter regions

Hundreds of thousands of putative BSs reside on Alu elements in the upstream regions of genes. For several TFs one observes that nearly all the putative BSs (more than 90%) are located on Alu. The search for BMs presented in the previous sections was done on the consensus sequences of various Alu subfamilies. Since, however, Alu elements have mutated after insertion, there is considerable deviation from the consensus in the specific Alu sequences that are found across the genome. This variation may change the BS content of specific individual Alu elements versus the consensus. Hence we repeated our search of BSs in the actual 5000 bp upstream regions of about 14,000 genes, and if a BS passed a significance filter, we checked whether it resides on a specific Alu. This search was limited to those 66 BMs that passed our filter based on the consensus sequences, i.e. their largest score, Tmax, exceeded T5 for the consensus of at least one of the eight Alu subfamilies studied. For each of these 66 BMs we searched the 14,000 upstream regions for putative BMs whose score exceeds (the consensus-based) Tmax.

We found that for 29 BMs out of the 66 the majority of putative BSs that pass Tmax reside on Alu sequences; the relative on and off Alu BM content is presented, for 19 out of these 66, in Fig 2. For example, for 5 out the 6 nuclear factors (Table 4) whose BMs appear among the 66, between 66% and 99% of the BMs (found in the 5000 bp upstream regions) are on Alu (see Fig. 2). The highest over-representations in Alu were observed for the BS of PITX2; over 99% of the 34000 BSs that pass the threshold were on Alu, and for LXRE's BSs over 99% of it's 9000 BSs that pass the threshold were on Alu, (see Fig. 2). In general, the distributions of putative BSs of TFs that bind to TGTAATCCC (Table 1) exhibit remarkable differences between the number of putative BS on and off Alu; over 90% of the total number of putative BSs of these TFs are located within Alu sequences. More moderate differences were found for developmental TFs such as SIX3.01, PAX4.01, NRL.01, MEF2.04 and NKX25.01 (Fig. 2). In general, for every BM that passed FDR of Q = 0.2 (see Table 1) over 80% of the BSs on the 5 kbp upstream region are located on Alu.

Figure 2
figure 2

Number of BSs of 19 TFs (selected from the list of 66 – see text) on the 5 Kbp region upstream of 13686 genes. A DNA sequence was considered to be a BS if its score was higher than Tmax. The BS were counted separately for on- Alu (blue) and off-Alu (red) appearances. For 16 out of 19 TFs the majority of the putative BSs reside on Alu. (A) developmental TFs (B) other TFs.

On the other hand, for the 37 (out of 66) remaining PSSMs the majority of BSs were not located on Alu, like IRF1.01 (48% out of the total number of 4000 BS were located on Alu) and HNF1.03 (27%). Still, even for these TFs, with mainly off-Alu BMs, the Alu repeats provide between hundreds to thousands of putative BSs in the 5 kbp upstream regions (see Table 6).

Table 6 Amount of binding motifs (for 25 TFs with Tmax > T1) in the 500 bp and 5000 bp regions of 13686 genes. For each BM we show statistics with respect to sequences that passed Tmax and T5 (see text for further details). We also checked whether the BS resides on Alu or not (see methods). and give the number of BSs that reside on Alu and their percentage out of the total number of BS.

As mentioned above, in the 500 bp upstream region, which constitutes the proximal promoter [46] and hence is believed to be most important for regulation, the Alu density is below the genome average (Fig. 1). To study this region in more detail, we repeated the analysis that was done for the 5000 bp region, and compared the distributions of BSs on and off Alu in the first 500 bp upstream to the TSS (see Fig 3 for 19 representative TFs). Only 9 out of the 66 BMs had more on Alu than off Alu BMs. These BMs include almost all homeodomain proteins, LXRE.01, LYF-1 and PAX4.01 (but not SIX3, see Fig. 3). In spite of the relatively low Alu density, on-Alu BMs do constitute a sizeable fraction of the putative BMs that are found in the 500 bp long upstream region (Fig. 3 and Table 6).

Figure 3
figure 3

Number of BSs of 19 TFs (see text) on the 500 bp region upstream of 13686 genes. A DNA sequence was considered to be a BS if its score was higher than Tmax. The BS were counted separately for on (blue) and off Alu (red) 7 appearance. For 7 out of 19 TFs the majority of the putative BSs reside on Alu. (A) developmental TFs (B) other TFs.

Genes of different biological processes contain different amounts of BSs as a result of differences in their Alu content

We used the results of the search for on and off Alu BMs to look for differences between genes that belong to different biological processes. As expected, for those biological processes whose genes were previously found to be enriched with Alu (Table 6), we also found that a high percentage of their regulatory signals were on-Alu. The Alu rich groups in this category include metabolic and bio-synthesis associated genes. Since for several TFs a large majority of BMs reside on Alu, we find that Alu-rich functional families have a large number of these Alu-associated BMs, whereas Alu-poor biological processes (such as developmental genes) will have less of these BMs (see Fig. 4 and Additional data file 2). For example, the average number of putative BSs of PITX2 is almost 4 times greater in protein biosynthetic genes than in CNS developmental genes (Fig. 4). Interestingly, as a result, developmental TFs such as LYF-1, PAX4, MEF2D and PITX2, that have on-Alu BMs, appear less in the (Alu-poor) upstream region of developmental genes than in the Alu-rich non-developmental group of protein biosynthesis related genes. One possible explanation of this may be that the developmental TFs main role is to suppress pathways that are in competition with or exclusive of the developmental process (versus acting directly on development – associated genes).

Figure 4
figure 4

We present for 6 TFs the number of putative BSs in the 5 Kbp region upstream, averaged over the genes that belong to various biological processes. The 6 TFs shown were chosen because their BSs are non-overlapping and their Tmax is greater than T5 (see text) a. for on-Alu BSs b. all BSs(on and off Alu).

One should emphasize that the statements made above are based mainly on in-silico evidence, which must be supplemented by firm experimental observations of binding and functional roles of the bound TFs.

Inter-species comparisons: B1, the murine 'Alu like' elements, contain several BSs that are similar to the primate Alu BSs

We scanned the B1 consensus sequences that were downloaded from RepBase [29]. These sequences were less enriched with BSs that pass the T5 threshold; this might be due to the fact that the length of the B1 sequences is only half that of Alu. Twenty three binding motifs passed T5 on B1 (see Table 7) and we found that 11 of these passed T5 also in Alu, as expected from the similarity between B1 and Alu [47] (bold in Table 7). One of these shared BSs was T(G/T)TAATCCC where G appears in all Alu families and T appears in the consensus of B1s. This BS is the target of OTX2.01, a homeodomain protein in both human and mouse. An additional homeodomain protein, CRX01.01, is also bound to this BS in mouse; it did not pass T5 in Alu. The members of the MEF2 family of BMs appear in mouse (see Table 7) with a weaker score than in Alu, (unpublished data). The PSSM with the most significant score on B1 was PAX6.02, which is not surprising since this PSSM is built from those BSs of PAX6 that were known to reside on B1 [48]. In addition, a group of BMs that are related to nuclear factors and response proteins passed T5 on B1 (TTF1.01 and ERR_01) but not on Alu. This is surprising since there are experiments that validated that the TFs that correspond to these BMs bind to primate Alu [16].

Table 7 BMs with target sites on B1_mM and B1_Mus ('Alu Like' sequence in mouse). The scores of the target sites were above T5 (those for which the corresponding ALU BS had scores greater than T5 are in bold). The IUPAC convention [33] consensus sites are in column 2, and the third column shows the target site with the highest score among all the major consensus sequence (nucleotides that agree with the consensus are denoted with capital letters). The location of the putative target sites is in column 4, given for the subfamily on which the target site of column 3 appears. Column 6 lists the number, n, of B1 subfamilies on which the BS shown in column 3 is located.

The distributions of B1 in mouse promoters and Alu in human promoters are positively correlated

When the mouse genome was first sequenced, regions orthologous to human DNA were identified. The densities of B1 and B2 in mouse and Alu in human were measured in windows of 500,000 bp and the Alu density in humans was found to be correlated with B2 and B1 density in mouse [49]. A high correlation was quite unexpected, since Alu and the mouse SINE elements B1 and B2 spread in the rodent and human genomes after the divergence of rodents and primates. We checked whether this trend is present also in the regulatory regions of genes, of size 5 kbp, and we found a positive correlation between B1 and ALU densities even in this relatively short region. In more than 5400 pairs of orthologous genes, that had available upstream sequences in both human and mouse, we found positive correlation (see methods) of 0.495 (pval<1e-05). In addition we ordered the genes by increasing Alu content in the 5000 upstream bp region, and divided the genes on this basis into 27 sets, with each set containing 200 genes. We then computed the average Alu and B1 associated nucleotide densities for each of the 27 sets, and plotted the two average Alu/B1 contents (Fig 5). Each point in Fig 5 represents a group of 200 genes, with the x (y) coordinates representing the Alu (B1) densities in human (mouse). The monotonic behaviour evident in Fig 5 indicates that the Alu and B1 contents of orthologous genes are largely similar in human and mouse. This similarity is far from being obvious, since spreading of Alu and B1 across the genome is believed to have happened after the divergence of primates and rodents.

Figure 5
figure 5

Scatter plot of nucleotide densities in 5 Kbp upstream of the TSS, that belong to B1 (in mouse genes) and to Alu (in orthologous human genes). For both species the values presented by a point are obtained by averaging over 200 genes that belong to a "gene set". There are 27 gene sets, that were obtained by sorting the 5400 human promoters (that have mouse othologues) according to their Alu content, and dividing them into 27 bins (of 200 genes in each). Each mouse gene set contains the orthologous genes of the corresponding human gene-set.

Conclusion

Summary

This research suggests that evolution used transposable elements to insert modules of transcription factor binding motifs into promoters and, by means of their presence, assemble higher level regulatory networks. In order to explore this question we focused on Alu elements, which are good potential candidates to be part of the building blocks of regulatory networks for two reasons. First, Alu elements are abundant in the upstream region of the TSS of genes, and second, Alu elements contain dozens of putative BSs for TFs. Some of these BSs were found before and their association with Alu was also reported, whereas in some cases although the BSs were found, the fact that they reside on Alu went unnoticed. Finally, we list here also BSs on Alu that were not identified previously. Our findings imply that the biological pathway on which Alu-mediated regulation appears to have the most significant impact is the development process. Many of the TFs that have binding motifs on Alu are associated with development; moreover, some of these BSs were previously demonstrated to be functional in vivo and essential to regulation of some target genes.

We have also shown that a few subsequences on Alu, that were known to be highly conserved, provide many of the binding sites for the transcription factors. We studied the Alu distribution in the upstream region of genes and found that while Alu are "repelled" from the 500 bp long region immediately preceding the transcription starts site, the 4500 bp long sequence further upstream is enriched (compared to the entire genome). The Alu content in the promoters varies between functional families of genes.

Comparison of Alu in human and B1 in mice shows a significant correlation between their distributions (with respect to orthologs), and high similarity between the binding motifs that reside on the two mobile elements.

Several particular cases (TFs and genes) were discussed in detail.

Future directions

Most importantly, more experiments are needed to verify that the BMs that were found by us (see Table 1 and Additional data file 1) indeed bind their corresponding TFs and, furthermore, to verify their regulatory roles. Of particular interest is the role of Alu in developmental processes. We found that promoters of genes involved in biosynthesis and metabolism are enriched (via Alu) in BMs of several developmental TFs (Table 4). It is our belief that since development and cell proliferation (of which biosynthesis is an essential component) are opposing processes [50], these BMs play an inhibitory role, suppressing proliferation during differentiation. This hypothesis is supported by the known suppressory role of GFI1, but the presently known activities on PITX2 and NKX2.5 do not fit this picture. Further experiments that address this issue by measuring the regulatory effects of other Alu-associated developmental TFs are needed. Clearly, many of the genes and their functions date earlier than the rodent/primate divergence, whereas the Alu insertions happened later. This poses an open and debated [20, 51] question about the mechanism by which pre-existing regulatory circuits must have been replaced, modified or supplemented by Alu insertions.

Methods

Downloads

The Alu consensus sequences were downloaded from Repbase [29].

To scan for BMs we downloaded 410 vertebrate PSSMs from the MathInspector [30] database, with two additional PSSMs ($PITX2_Q2 and &LUN1_01) taken from the TRANSFAC database [31].

The coordinates of the upstream regions of 21000 human genes were taken from [52]. The annotations of the U133 Array (Affymetrix, Santa Clara, CA) were used in order to remove multiple upstream regions for different splice variants of the same gene, (a single upstream region was associated with each gene symbol). We also removed from our database those Refseq genes that have duplicates on chromosomes X and Y, and those with upstream regions that were not fully sequenced. All this left 13686 promoters of different Refseq genes for analysis.

The coordinates of Alu elements and other transposable elements in the human genome were extracted from the file chromOut.zip, which was downloaded from [52]. This is the output file of RepeatMasker [53]. The sequences of the Alu that reside in promoters were then extracted from the human genome, by intersecting these Alu coordinates with the upstream sequences of our 13686 genes.

Similarly, for mouse we downloaded from [54] the files upstream5000.zip which contained the coordinates and sequences of the regions upstream to the TSS of mouse genes, and the RepeatMasker output file of the mouse genome, chromOut.zip. We used annotations of Mouse 430_2 Array (Affymetrix, Santa Clara, CA) in order to remove multiple promoters.

In order to annotate genes to biological processes and molecular functions we used the DAVID database (based on the GO annotation) [55].

Calculating scores and probabilities

For each PSSM we calculated the maximal score, on both sense and antisense orientations.

In order to find the probability of getting some maximal score for a BM on an Alu sequence we produced 100,000 random 281 bp long sequences, of di-nucleotide distribution similar to that of the AluSx consensus sequence. The next step was to calculate the maximal score for a given PSSM for each random sequence. The p-value for each given score was the percent of random sequences that contained a target site with higher score. Although, such procedure might seem to be biased towards AluSx, it can be seen in Tables 8 and 9 that the subfamilies on which we focus in this research have similar nucleotide and di- nucleotide contents. Note that since we did not use a twice repeated random sequence of 140 nucleotides, our p-values are an upper bound to the true ones.

Table 8 The nucleotides distribution in the consensus sequences of major Alu subfamilies.
Table 9 Number of di -nucleotide appearances in the consensus sequences of major Alu subfamilies.

Orthologous genes in mouse and human

We compared gene symbols of human and mouse and got about 5400 orthologous gene pairs and their respective promoters in human and mouse. The Pearson's correlation coffeeicient between the Alu and B1 densities on human and mouse genes was calculated over these 5400 pairs of promoters.

References

  1. Lander ES: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.

    Article  PubMed  CAS  Google Scholar 

  2. Kidwell MG, Lisch DR: Perspective: transposable elements, parasitic DNA, and genome evolution. Evolution Int J Org Evolution. 2001, 55: 1 -24.

    Article  CAS  Google Scholar 

  3. Batzer MA, Deininger PL: Alu repeats and human genomic diversity. Nat Rev Genet. 2002, 3: 370-379. 10.1038/nrg798.

    Article  PubMed  CAS  Google Scholar 

  4. Dewannieux M, Esnault C, Heidmann T: LINE-mediated retrotransposition of marked Alu sequences. Nat Genet. 2003, 35: 41-48. 10.1038/ng1223.

    Article  PubMed  CAS  Google Scholar 

  5. Batzer MA, Deininger PL, Hellmann-Blumberg U, Jurka J, Labuda D, Rubin CM, Schmid CW, Zietkiewicz E, Zuckerkandl E: Standardized nomenclature for Alu repeats. Journal of Molecular Evolution (Historical Archive). 1996, 42: 3-6. 10.1007/BF00163204.

    Article  CAS  Google Scholar 

  6. Kapitonov V, Jurka J: The age of Alu subfamilies. Journal of Molecular Evolution. 1996, 42: 59-65. 10.1007/BF00163212.

    Article  PubMed  CAS  Google Scholar 

  7. Medstrand P, van de Lagemaat LN, Dunn CA, Landry JR, Svenback D, Mager DL: Impact of transposable elements on the evolution of mammalian gene regulation. Cytogenet Genome Res. 2005, 110: 342-345. 10.1159/000084966.

    Article  PubMed  CAS  Google Scholar 

  8. Grover D, Mukerji M, Bhatnagar P, Kannan K, Brahmachari SK: Alu repeat analysis in the complete human genome: trends and variations with respect to genomic composition. Bioinformatics. 2004, 20: 813-817. 10.1093/bioinformatics/bth005.

    Article  PubMed  CAS  Google Scholar 

  9. Greally JM: Short interspersed transposable elements (SINEs) are excluded from imprinted regions in the human genome. PNAS. 2002, 99: 327-332. 10.1073/pnas.012539199.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  10. Mariño-Ramirez L, Lewis KC, Landsman D, Jordan IK: Transposable elements donate lineage-specific regulatory sequences to host genomes. Cytogenet Genome Res. 2005, 110: 333-341. 10.1159/000084965.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Jordan IK, Rogozin IB, Glazko GV, Koonin EV: Origin of a substantial fraction of human regulatory sequences from transposable elements. Trends Genet. 2003, 19: 68 -672. 10.1016/S0168-9525(02)00006-9.

    Article  PubMed  CAS  Google Scholar 

  12. Grover D, Majumder PP, Rao CB, Brahmachari SK, Mukerji M: Nonrandom distribution of alu elements in genes of various functional categories: insight from analysis of human chromosomes 21 and 22. Mol Biol Evol. 2003, 20: 1420 -11424. 10.1093/molbev/msg153.

    Article  PubMed  CAS  Google Scholar 

  13. Ganapathi M, Srivastava P, Sutar S, Kumar K, Dasgupta D, Pal Singh G, Brahmachari V, Brahmachari S: Comparative analysis of chromatin landscape in regulatory regions of human housekeeping and tissue specific genes. BMC Bioinformatics. 2005, 6: 126-10.1186/1471-2105-6-126.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Kazazian HHJ: Mobile elements: drivers of genome evolution. Science. 2004, 303: 1626-1632. 10.1126/science.1089670.

    Article  PubMed  CAS  Google Scholar 

  15. Grover D, Kannan K, Brahmachari SK, Mukerji M: ALU-ring elements in the primate genomes. Genetica. 2005, 124: 273-289. 10.1007/s10709-005-3086-8.

    Article  PubMed  CAS  Google Scholar 

  16. Babich V, Aksenov N, Alexeenko V, Oei SL, Buchlow G, Tomilin N: Association of some potential hormone response elements in human genes with the Alu family repeats. Gene. 1999, 239: 341 -3349. 10.1016/S0378-1119(99)00391-1.

    Article  PubMed  CAS  Google Scholar 

  17. Le Goff W, Guerin M, Chapman MJ, Thillet J: A CYP7A promoter binding factor site and Alu repeat in the distal promoter region are implicated in regulation of human CETP gene expression. J Lipid Res. 2003, 44: 902 -9910. 10.1194/jlr.M200423-JLR200.

    Article  PubMed  CAS  Google Scholar 

  18. Norris J, Fan D, Aleman C, Marks JR, Futreal PA, Wiseman RW, Iglehart JD, Deininger PL, McDonnell DP: Identification of a new subclass of Alu DNA repeats which can function as estrogen receptor-dependent transcriptional enhancers. J Biol Chem. 1995, 270: 22777-22782. 10.1074/jbc.270.39.22777.

    Article  PubMed  CAS  Google Scholar 

  19. Blunt T, Finnie NJ, Taccioli GE, Smith GC, Demengeot J, Gottlieb TM, Mizuta R, Varghese AJ, Alt FW, Jeggo PA: Defective DNA-dependent protein kinase activity is linked to V(D)J recombination and DNA repair defects associated with the murine scid mutation. Cell. 1995, 80: 813-823. 10.1016/0092-8674(95)90360-7.

    Article  PubMed  CAS  Google Scholar 

  20. Oei SL, Babich VS, Kazakov VI, Usmanova NM, Kropotov AV, Tomilin NV: Clusters of regulatory signals for RNA polymerase II transcription associated with Alu family repeats and CpG islands in human promoters. Genomics. 2004, 83: 873-882. 10.1016/j.ygeno.2003.11.001.

    Article  PubMed  CAS  Google Scholar 

  21. Piedrafita FJ, Molander RB, Vansant G, Orlova EA, Pfahl M, Reynolds WF: An Alu element in the myeloperoxidase promoter contains a composite SP1-thyroid hormone-retinoic acid response element. J Biol Chem. 1996, 271: 14412 -114420. 10.1074/jbc.271.24.14412.

    Article  PubMed  CAS  Google Scholar 

  22. Vansant G, Reynolds WF: The consensus sequence of a major Alu subfamily contains a functional retinoic acid response element. Proc Natl Acad Sci U S A. 1995, 92: 8229 -88233. 10.1073/pnas.92.18.8229.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  23. Britten RJ: Evolutionary selection against change in many Alu repeat sequences interspersed through primate genomes. Proc Natl Acad Sci U S A. 1994, 91: 5992 -55996. 10.1073/pnas.91.13.5992.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  24. Thornburg BG, Gotea V, Makalowski W: Transposable elements as a significant source of transcription regulating signals. Gene. 2006, 365: 104-110. 10.1016/j.gene.2005.09.036.

    Article  PubMed  CAS  Google Scholar 

  25. Stepanova M, Tiazhelova T, Skoblov M, Baranova A: A comparative analysis of relative occurrence of transcription factor binding sites in vertebrate genomes and gene promoter areas. Bioinformatics. 2005, 21: 1789-1796. 10.1093/bioinformatics/bti307.

    Article  PubMed  CAS  Google Scholar 

  26. Hamdi HK, Nishio H, Tavis J, Zielinski R, Dugaiczyk A: Alu-mediated phylogenetic novelties in gene regulation and development. J Mol Biol. 2000, 299: 931 -9939. 10.1006/jmbi.2000.3795.

    Article  PubMed  CAS  Google Scholar 

  27. Hambor JE, Mennone J, Coon ME, Hanke JH, Kavathas P: Identification and characterization of an Alu-containing, T-cell- specific enhancer located in the last intron of the human CD8 alpha gene. Mol Cell Biol. 1993, 13: 7056-7070.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Hanke JH, Hambor JE, Kavathas P: Repetitive Alu elements form a cruciform structure that regulates the function of the human CD8[alpha] T cell-specific enhancer. Journal of Molecular Biology. 1995, 246: 63-73. 10.1006/jmbi.1994.0066.

    Article  PubMed  CAS  Google Scholar 

  29. Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467. 10.1159/000084979.

    Article  PubMed  CAS  Google Scholar 

  30. Quandt K, Frech K, Karas H, Wingender E, Werner T: MatInd and MatInspector: new fast and versatile tools for detection of consensus matches in nucleotide sequence data. Nucleic Acids Res. 1995, 23: 4878-4884.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  31. Wingender E, Chen X, Hehl R, Karas H, Liebich I, Matys V, Meinhardt T, Pru{beta} M, Reuter I, Schacherer F: TRANSFAC: an integrated system for gene expression regulation. Nucl Acids Res. 2000, 28: 316-319. 10.1093/nar/28.1.316.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  32. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc. 1995, 57: 289-300.

    Google Scholar 

  33. Cavener DR: Comparison of the consensus sequence flanking translational start sites in Drosophila and vertebrates. Nucl Acids Res. 1987, 15: 1353-1361.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Medstrand P, van de Lagemaat LN, Mager DL: Retroelement Distributions in the Human Genome: Variations Associated With Age and Proximity to Genes. Genome Res. 2002, 12: 1483-1495. 10.1101/gr.388902.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  35. Chesnokov IN, Schmid CW: Specific Alu binding protein from human sperm chromatin prevents DNA methylation. J Biol Chem. 1995, 270: 18539-18542. 10.1074/jbc.270.31.18539.

    Article  PubMed  CAS  Google Scholar 

  36. Eisenberg E, Adamsky K, Cohen L, Amariglio N, Hirshberg A, Rechavi G, Levanon EY: Identification of RNA editing sites in the SNP database. Nucl Acids Res. 2005, 33: 4612-4617. 10.1093/nar/gki771.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  37. Kim DDY, Kim TTY, Walsh T, Kobayashi Y, Matise TC, Buyske S, Gabriel A: Widespread RNA editing of embedded Alu elements in the human transcriptome. Genome Res. 2004, 14: 1719-1725. 10.1101/gr.2855504.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  38. Athanasiadis A, Rich A, Maas S: Widespread A-to-I RNA editing of Alu-containing mRNAs in the human transcriptome. PLoS Biology. 2004, 2: e391-10.1371/journal.pbio.0020391.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Chu D, Kakazu N, Gorrin-Rivas MJ, Lu HP, Kawata M, Abe T, Ueda K, Adachi Y: Cloning and characterization of LUN, a novel RING finger protein that Is highly expressed in lung and specifically binds to a palindromic sequence. J Biol Chem. 2001, 276: 14004-14013.

    PubMed  CAS  Google Scholar 

  40. Hjalt TA, Amendt BA, Murray JC: PITX2 regulates procollagen lysyl hydroxylase (PLOD) gene expression: implications for the pathology of Rieger syndrome. J Cell Biol. 2001, 152: 545-552. 10.1083/jcb.152.3.545.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  41. Duan Z, Horwitz M: Targets of the transcriptional repressor oncoprotein Gfi-1. PNAS. 2003, 100: 5932-5937. 10.1073/pnas.1031694100.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  42. Ganga M, Espinoza HM, Cox CJ, Morton L, Hjalt TA, Lee Y, Amendt BA: PITX2 Isoform-specific regulation of atrial natriuretic factor expression: synergism and repression with Nkx2.5. J Biol Chem. 2003, 278: 22437-22445. 10.1074/jbc.M210163200.

    Article  PubMed  CAS  Google Scholar 

  43. Rothenberg EV, Taghon T: Molecular genetics of T cell devlopment. Annual Review of Immunology. 2005, 23: 601-649. 10.1146/annurev.immunol.23.021704.115737.

    Article  PubMed  CAS  Google Scholar 

  44. Brand T: Heart development: molecular insights into cardiac specification and early morphogenesis. Developmental Biology. 2003, 258: 1-19. 10.1016/S0012-1606(03)00112-X.

    Article  PubMed  CAS  Google Scholar 

  45. Zhou YH, Zheng JB, Gu X, Saunders GF, Yung WKA: Novel PAX6 binding sites in the human genome and the role of repetitive elements in the evolution of gene regulation. Genome Res. 2002, 12: 1716-1722. 10.1101/gr.188302.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  46. Zhang MQ: Computational prediction of eukayotic prtoein-coding genes. Nat Rev Genet. 2002, 3: 698-709. 10.1038/nrg890.

    Article  PubMed  CAS  Google Scholar 

  47. Vassetzky NS, Ten OA, Kramerov DA: B1 and related SINEs in mammalian genomes. Gene. 2003, 319: 149-160. 10.1016/S0378-1119(03)00805-9.

    Article  PubMed  CAS  Google Scholar 

  48. Zhou YH, Zheng JB, Gu X, Li WH, Saunders GF: A novel Pax-6 binding site in rodent B1 repetitive elements: coevolution between developmental regulation and repeated elements?. Gene. 2000, 245: 319-328. 10.1016/S0378-1119(00)00019-6.

    Article  PubMed  CAS  Google Scholar 

  49. Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420: 520-562. 10.1038/nature01262.

  50. Zhu L, Skoultchi AI: Coordinating cell proliferation and differentiation. Current Opinion in Genetics & Development. 2001, 11: 91-97. 10.1016/S0959-437X(00)00162-3.

    Article  CAS  Google Scholar 

  51. Brosius J: RNAs from all categories generate retrosequences that may be exapted as novel genes or regulatory elements. Gene. 1999, 238: 115-134. 10.1016/S0378-1119(99)00227-9.

    Article  PubMed  CAS  Google Scholar 

  52. UCSC Genome Browser FTP site. [http://hgdownload.cse.ucsc.edu/goldenPath/hg17/bigZips]

  53. RepeatMasker. [http://www.repeatmasker.org/]

  54. UCSC Genome Browser FTP. [http://hgdownload.cse.ucsc.edu/goldenPath/mm5/bigZips/]

  55. Dennis G, Sherman B, Hosack D, Yang J, Gao W, Lane HC, Lempicki R: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biology. 2003, 4: P3-10.1186/gb-2003-4-5-p3.

    Article  PubMed  Google Scholar 

  56. Li Y, Bolten C, Bhat BG, Woodring-Dietz J, Li S, Prayaga SK, Xia C, Lala DS: Induction of human Liver X Receptor {alpha} gene expression via an autoregulatory loop mechanism. Mol Endocrinol. 2002, 16: 506-514. 10.1210/me.16.3.506.

    PubMed  CAS  Google Scholar 

Download references

Acknowledgements

The authors wish to thank Tal Shay for help identifying the Alu elements and with the intial genomic mapping of Alu elements, Yuval Tabach for useful comments and fruitful discussions, Or Zuk and Assif Yitzhaky for help with writing computer scripts for searching BSs on the genes' promoters. This work was partially supported by the Ridgefield Foundation and by the Wolfson Family Charitable Trust, London, on Tumour Cell Diversity.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Eytan Domany.

Additional information

Authors' contributions

PP participated in the conception of the study, carried out the mapping of the Alu elements and BSs in the upstream of genes, analyzed the data, wrote the initial draft of the manuscript. ED conceived of the study, and participated in its design coordination and supervision, helped to draft and edited the manuscript.

Electronic supplementary material

12864_2006_516_MOESM1_ESM.xls

Additional File 1: Sixty six binding motifs whose scores pass T5 on at least one major Alu subfamily sequence. The consensus sites of the PSSMs (from the IUPAC convention [33]) is in the second column and the the target sequences with the highest scores among all the major Alu consensus sequences – in the third, with nucleotides that agree with the consensus denoted by capital letters. The fourth column contains the locations of the putative target sites on the Alu sequences of the subfamily with the highest score (in case the same sequence appears on several Alu elements, we choose the one with the highest number of copies in the 5 kb upstream regions). Two Alu subsequences serve as putative target sites of several TFs (designated by *, and **). The fifth column contains the p-values (see text and methods) and the number of subfamilies on which the BSs of the third column resides is listed in column 6. (XLS 24 KB)

12864_2006_516_MOESM2_ESM.xls

Additional File 2: We present for 66 PSSMs the number of putative BSs in the 5 Kbp region upstream, averaged over the genes that belong to various biological processes. The PSSMs that are shown are those for which Tmax>T5. The number of genes in each biological process is given in column 2. The Alu density and number of Alu repeats per gene are given in columns 3 and 4 respectively. (XLS 128 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Polak, P., Domany, E. Alu elements contain many binding sites for transcription factors and may play a role in regulation of developmental processes. BMC Genomics 7, 133 (2006). https://doi.org/10.1186/1471-2164-7-133

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-7-133

Keywords