Skip to main content

Genome-wide analysis of the transcription factor binding preference of human bi-directional promoters and functional annotation of related gene pairs



Bi-directional gene pairs have received considerable attention for their prevalence in vertebrate genomes. However, their biological relevance and exact regulatory mechanism remain less understood. To study the inner properties of this gene organization and the difference between bi- and uni-directional genes, we conducted a genome-wide investigation in terms of their sequence composition, functional association and regulatory motif discovery.


We identified 1210 bi-directional gene pairs based on the GRCh37 assembly data, accounting for 11.6% of all the human genes owning RNAs. CpG islands were detected in 98.42% of bi-directional promoters and 61.07% of unidirectional promoters. Functional enrichment analysis in GO and GeneGO both revealed that bi-directional genes tend to be associated with housekeeping functions in metabolism pathways and nuclear processes, and 46.84% of the pair members are involved in the same biological function. By fold-enrichment analysis, we characterized 73 and 43 putative transcription factor binding sites(TFBS) that preferentially occur in bi-directional promoters from TRANSFAC and JASPAR database respectively. By text mining, some of them were verified by individual experiments and several novel binding motifs were also identified.


Bi-directional promoters feature a significant enrichment of CpG-islands as well as a high GC content. We provided insight into the function constraints of bi-directional genes and found that paired genes are biased toward functional similarities. We hypothesized that the functional association underlies the co-expression of bi-directional genes. Furthermore, we proposed a set of putative regulatory motifs in the bi-directional promoters for further experimental studies to investigate transcriptional regulation of bi-directional genes.


According to the orientation and status of the 5’ end, the adjacently located genes can be arranged in convergent, divergent or tandem configuration[1]. Among these categories, the divergent gene arrangement is found more frequently than expected by chance in the human genome, accounting for about 10% of all human genes[2, 3]. Bi-directional gene pair is defined as two genes arranged in a head-to-head (adjacent 5’ ends) fashion on opposite strands of DNA with less than 1,000 bp between their transcription start sites(TSS)[1]. Accordingly, the entire intervening region between the two TSSs is designated as a putative bi-directional promoter. A gene is termed as uni-directional if no oppositely oriented TSS is found within 10 kb upstream of the given TSS, or if a similarly oriented TSS is found at least 1 kb upstream. Thus the entire 1 kb of 5’ flanking DNA is considered as the uni-directional promoter.

Considerable attention has been focused on bi-directional genes in recent years. Examples including LRRC49/THAP10[4], SURF-1/SURF-2[5], COL4A1/COL4A2[6], PCD10/SERPINI1[7] and HAND2/DEIN[8] have been identified in human through individual experiments. A considerable number of bi-directional gene pairs were found to be conserved among mammalian species[9, 10]. Since evolutionary conservation usually indicates functional implications, we proposed that bi-directional gene organization is under selection to fulfil a specific functional role. Whereas most of the bi-directional gene pairs have been found in the process of studying a single gene, a genome-wide analysis of their function and physiologic consequences is currently insufficient.

The expression data obtained from biotechnologies such as SAGE and microarray indicated a correlated expression profile between bi-directional genes[1113]. Based on the assumption that ‘co-expression implies co-regulation’, the requirement for co-regulation of functionally related genes appears to underlie the observed co-expression. However, it is still under discussion whether the co-expression evolved merely as a consequence of their physical proximity or if function dictated their co-regulation. There are several examples of bi-directional gene pairs that are related by function, e.g. in DNA repair[1, 2], aging[14], de novo purine synthesis[15] and carcinogenesis[5]. Despite this observation, a systematic study on the degree of internal co-function of the bi-directional genes has not been carried out to date.

More recent studies have suggested an intrinsic difference in nucleotide composition of bi-directional promoters compared to uni-directional ones[1, 2, 13, 16]. These characteristic feature lead us to hypothesize that divergent genes will be transcribed with a special set of regulatory signals. Currently our understanding of transcription regulation relies greatly on experimental identification of prospective regulatory regions. Yet many specifics underlying the regulatory design are unknown. Therefore, it seems necessary to re-evaluate the underlying mechanisms and biological relevance of bi-directional promoters systematically.

In the present study, we have undertaken a genome-wide survey of gene organization in the human genome. To reveal functions collectively performed by such bi-directional genes, we mapped them to the Gene Ontology (GO) and GeneGO pathways. We also explored the functional similarity between the genes on the plus strand and those on the minus strand. We devoted our effort into exploring the binding preference of transcription factors on the bi-directional promoters and statistically identified a set of over-represented transcription factor binding sites(TFBS) in bi-directional promoters, the research scheme is shown in figure 1.

Figure 1
figure 1

Diagrammatic Representation of the Research Scheme


Identification of bi-directional and uni-directional genes

We calculated the distances between the transcription start sites of nearest gene neighbors for four kinds of combination of the clusters and the result was showed in table 1. We identified 1210 bi-directional gene pairs based on the curated transcript cluster NMs and NRs, accounting for 11.67% of all the genes owning RNAs, which agrees the view that bidirectional gene pairs are prevalent in the human genome. If only transcripts with conclusive mRNA were reserved, 878 bi-directional gene pairs, in the proportion of 9.31%, were discovered upon the removal of pairs consisting of NMs and NRs. Redundant gene pair entries that share the same intergenic sequence were removed to yield 822 bi-directional gene pairs for the analyses.

Table 1 Distribution of bi-directional gene pairs on each chromosome

CpG islands are preferentially located in bi-directional promoters

There have been contradictory observations on the CpG island frequency in bi-directional promoters. Adachi et al.[1, 16] considered the presence of CpG island to be a common feature of bidirectional promoters. In contrast, Takai, et al. [13] reported that CpG islands are not preferentially associated with bidirectional promoters. The author attributed the discrepancy to the different criteria adopted to define a CpG island. In order to rationalize these controversial observations, we performed genome-wide computational analysis of the bi-directional promoters on the basis of two different definition systems. According to traditional definition by Gardiner-Garden[17], CpG islands were detected in 809 bi-directional promoters, representing 98.42% of a total of 822 pairs. A lower percentage of 61.07% was recorded for unidirectional promoters. Based on more strict criteria[18] (DNA fragment no less than 500 bp with GC-content >= 55% and Obs/Exp value >=0.60), CpG-islands were present in 86.37% of bidirectional promoters compared to 28.48% of uni-directional promoters. In addition, we analyzed pure IG sequence to remove the difference caused by the extended IG region. Invariably the frequency of CpG island in bi-directional promoters is higher than those in unidirectional ones. As shown in Figure 2, the CpG density in bidirectional promoters (histogram in top left) is significantly higher than that in unidirectional promoter (histogram in top right) in all comparisons. Consistent with a significant enrichment of CpG-islands, bidirectional promoters feature a high C+G content.

Figure 2
figure 2

Density Distribution of CpG Islands between Bi-directional Genes and Unidirectional Genes(The figure is reproduced with permission from the rights owner Liu,B.[28])

Functional Enrichment of Bi-directional Genes

Gene ontology associated with bi-directional promoter regulation

Genes regulated by bi-directional promoters were examined for functional classifications and associations. Among the 1,644 genes involved in the 822 human bi-directional gene pairs, 1,121, 1,219, and 1,256 genes were directly annotated by ‘biological process’, ‘molecular function’ and ‘cellular component’ subcategories in GO annotation system, respectively. We found several GO classes significantly over-represented among bi-directional genes. Cellular, metabolic and biosynthetic processes emerged as the most significantly enriched functional class. GO items of cell cycle and its child nodes were also significantly presented. Cellular response to stress or stimulus and their related subclasses of damage response, break repair were also focused. To summarize, the most enriched GO categories correspond to the known physiological roles of the cell, indicating that bi-directional genes are frequently involved in basic cellular metabolic processes. See Additional file 1 for the complete list of enriched GO terms.

Functional similarities among the paired genes

Among 822 annotated bi-directional gene pairs, we found 385 pairs (46.84%) whose plus and minus genes share at least one GO term, with SARS2/MRPS12 having most GO terms in common. Number and percentage of gene pairs that overlapped for various GO terms were showed in Figure 3. (See Additional file 2 for detailed list of shared GO terms). Such shared or related function supports the hypothesis that bi-directional genes are more likely to be functionally associated than uni-directional genes.

Figure 3
figure 3

Number and percentage of gene pairs that overlapped for various GO terms The X-axis represented the number of shared GO terms per gene pair. Number and percentage of gene pairs that overlapped by various GO terms among the total 385 pairs are plotted on the Y-axis.

We also provided separate estimates for each of the Gene Ontologies. We obtained 337 annotated pairs in subcategory "cellular component", 185 pairs in "molecular function" and 146 pairs in "biological process" respectively. It’s observed that, in general, bi-directional gene products are more likely to perform coordinated roles in the same cellular component, compared to the other two subsystems. Figure 4 illustrated the shared GO terms and P values in subcategory "cellular component".

Figure 4
figure 4

The shared GO terms and P value in Biological Process subcategory Red boxes represent GO terms that are occupied exclusively by genes on the plus strand; Green boxes represent GO terms that are occupied exclusively by genes on the minus strand, while yellow ones were common terms shared by plus and minus genes within the bi-directional gene pair. The color darkens with the significance of enrichment.

Then we set out to find out the GO terms that represent coordinated functions of bi-directional pairs. In Biological Process, the GO terms related to metabolic process and its branch such as primary metabolic process, cellular process and biopolymer biosynthetic process topped the list of both gene pair members. Their child nodes were focused on RNA (mRNA, ncRNA) metabolic process, cellular (macromolecule or biopolymer) catabolic process, organelle organization, mitotic cell cycle etc. In molecular function, the GO terms involved in DNA-directed RNA polymerase activity, RNA methyltransferase activity, purine NTP-dependent helicase activity, NAD or NADH binding, NADH dehydrogenase (quinone) activity, etc. are significantly over-represented as compared to others. In Cellular Component, we found that bi-directional genes tend to be tightly associated into the same class of organelle, organelle envelope, nucleus, nucleoplasm, nucleolus, membrane-bounded or non-membrane-bounded organelle, etc. Interestingly, almost all the items shared by the two divergent genes are related to metabolism and energy transfer. We proposed that genes involved in functions including metabolism, are more likely to be organized in the head-to-head configuration.

GeneGO pathway enrichment

Based on the P values from MetaCore™, totally we found 45 pathways that are significantly enriched with bi-directional genes out of the total 451 distinct pathways. According to the different classification criterion, the 45 pathways were assigned to 18 regulatory processes, 8 protein function, 4 disease maps and 15 metabolic maps. Extreme enrichment occurred for, in order of descending significance level, NHEJ mechanisms of DSBs repair, Oxidative phosphorylation, Nucleotide excision repair and GTP-XTP metabolism, Chromosome condensation in prometaphase, Role of Brca1 and Brca2 in DNA repair. Enriched pathways are further clustered into larger functional categories according to the GeneGO annotation. Regulatory processes/Cell cycle and Regulatory processes/DNA-damage ranked among the top enriched functional categories. Table 2 lists some most enriched categories ordered in decreasing level of significance.

Table 2 Statistically enriched GeneGO Pathway categories

Functional enrichment in GeneGO versus GO

So far we have been analyzing the level of gene function enrichment using two function annotation schemes respectively. The GO results show a clear agreement with those derived from the GeneGO pathways. For example, the GO terms that are significantly enriched include genes that are engaged in processes such as DNA metabolic process, which correspond to the Metabolic maps/Metabolic maps (common pathways)/Nucleotide metabolism pathway in GeneGO; Cell cycle, which corresponds to the same pathways in GeneGO; response to DNA damage stimulus, which corresponds to Regulatory processes/DNA-damage in GeneGO. This agreement is also apparent in that ‘‘DNA repair’’ is the most enriched GO term and ‘‘* DNA damage_Nucleotide excision repair,’’ which corresponds to the Regulatory processes/DNA-damage pathway, is one of the top enriched pathways in GeneGO as well.

Bi-directional promoters are characterized by a distinct collection of putative transcription factor binding sites

We characterized the enrichment of known motifs from TRANSFAC and JASPAR in bi-directional promoters relative to background uni-directional promoters. Based on the Jaspar PSSM information, we categorized 43 transcription factors as over-represented and 6 as under-represented. In the TRANSFAC database, 73 TFBSs found increased presence in bi-directional promoters. Complete lists of over-represented motifs and their enrichment folds are provided as Additional file 3. Although there is slight difference between the two databases, a large majority of the TFBSs overlap. The overlapped TFBSs and their over-represented folds were illustrated in Figure 5. We hypothesize that over-represented motifs correspond to transcription factors that are more likely to bind to bidirectional promoters than to unidirectional promoters. In contrast, under-represented motifs correspond to transcription factors that preferentially regulate unidirectional promoters. Shared motifs show no preference.

Figure 5
figure 5

The over-represented TFBSs that overlapped in TRANSFAC and JASPAR

We further investigated the experimental evidence supporting the roles of these transcriptional factors in regulating certain bi-directional genes. Table 3 lists the experimentally validated TFBS that occurred in bi-directional promoters. Some of the reported physiological functions are consistent with our functional enrichment analysis. For example, previous work[19] has demonstrated that GABPA regulates genes involved in a variety of cellular processes including adipocyte differentiation, mitochondrial respiration, and neuromuscular signaling, corresponding to enriched GO terms of cell cycle, cellular and metabolic processes and their child nodes. E2F1 are observed to regulate cell growth during the G0/G1-S phase transition, and over-expression of E2F1 induces apoptosis and DNA synthesis in quiescent fibroblasts[20]. These are in agreement with the significantly enriched GeneGO pathways such as Regulatory processes/Cell cycle and Regulatory processes/DNA-damage.

Table 3 The experimentally validated TFBS that occurred in bi-directional promoters

Interestingly, the over-represented recognition sequence for MYC, ELK1, NF-Y, SP1, ATF, GABPA, SREBP-1, NF-E2, STAT5A, NF-1 as well as SOX-9 rank among the most conserved motifs found in human promoters[21].

Given the enrichment of these motifs in bi-directional promoters and their strong evolutionary conservation across mammalian promoters, we assume that the predicted TFBSs located within bi-directional promoters are more likely to be functional in co-regulation than other TFBSs. Interestingly, it would appear that TFs within the same family tend to have similar binding preference. A TFBS is either over-represented or under-represented in parallel with other family members. These observations suggest a common mode of expression across the family members of transcription factors.


In this study, 11.6% of the human genes were shown to be arranged in a head-to-head fashion, and this proportion is slightly larger than most of the previous report[2], except that Piontkivska et al.[3] reported a number of 1,369 bi-directional promoters. The inconsistency was partly due to the update of TSS coordinates during the accumulation of EST and mRNA evidence. In addition, we used the much more highly curated RefGene track instead of the spliced human ESTs collection, because the large and complicated ESTs data containing thousands of transcripts captured by oligo-capping techniques will lead to an overestimation of the frequency of transcripts, and then introduce false positive result. What’s more, our work focus on the pure mRNA gene pairs and a large part of non-coding RNA, transcribed RNA and miscRNA are excluded from further analysis. Herein we provided a solid evidence for the previous observation[1] that bi-directional promoters had a significant enrichment of CpG-islands as well as a high GC content. Since CpG island is usually the targets of regulation by methylation, it may induce changes in chromatin structure that can confer either positive or negative effects on transcription. Misregulation of bi-directional promoters elicited by mutation or hypermethylation will simultaneously silence genes on both sides. Loss of their vital biological function well explains the role bi-directional genes in the development of human diseases such as aging[14], brain disease[7] and oncogenesis[4].

Our study provided insight into the function constraints of bi-directional genes. Functional enrichment analysis in GO and GeneGO both revealed that bi-directional genes are often associated with housekeeping functions. GO terms, including metabolic process, such as DNA, RNA, biopolymer or macromolecule metabolism, as well as nuclear processes, such as DNA repair and replication or cell cycle regulation are significantly enriched. The GeneGO pathways that are involved in growth or proliferation, such as those engaged in Energy metabolism, Nucleotide metabolism, Vitamin and cofactor metabolism, tend to be more enriched with bi-directional genes. Pathways in genetic information processing (transcription, translation and DNA repair) and cell cycle tend to be enriched as well. To summarize, bi-directional genes are significantly enriched in housekeeping functions such as metabolism pathways and nuclear processes.

Further analyses revealed that the significant functional categories are more likely to be shared by bi-directional genes. This indicated that the bi-directional genes are strongly biased toward functional similarities and coordinated regulation. We postulate that for bi-directional genes involved in basic biological processes, coordinated regulation ensures their synchronized action and thus minimizes transcriptional error. In contrast, genes with less coordinated regulation may be involved in pathways that are more flexible in responding to environmental changes.

We compared the TFBSs between bi- and uni-directional promoters according to their rate of occurrence. We discovered several transcription factors that preferentially regulate bi-directional promoters. Some of the TFBSs matched well with experimentally determined ones and several novel binding motifs were also identified. These bi-directional gene associated motifs may be envisaged as the best candidates for functional regulatory elements. In addition, the motif search result could help identify novel genes, which is linked to a known gene via a bi-directional promoter. And these genes probably perform important conserved functions.

We are also aware of some limitations in our analysis. The motifs for the identification of TFBSs are still incomplete, and the evolutionary importance of the over-representation of TFBS remains to be elucidated. Although some of their function are indicated by functional categories (GO terms) of experimental verified motifs, conclusive evidence of the role played by regulatory factors in the co-regulation of the two genes will be tested in experiments. Eventually, the combination of computational and experimental approaches will permit us to construct mechanistic models of regulatory transcription networks of bi-directional genes. It would be interesting, as a future endeavor, to examine these regulatory elements in other species in a similar fashion and compare the results to those obtained herein. Comparative analyses of these regulators across multiple species will validate our predictions by their appearance in another species. A related work is still in progress.


In this work, we conducted a systematic investigation of bi-directional gene organization focusing on sequence features, functional association and regulatory motif discovery. We confirmed known properties of bi-directional gene organization and also provided new observations. We found that bi-directional gene pairs show a higher probability to be functionally associated, formulating hypotheses that the requirement for co-regulation of functionally related genes is a possible cause for the observed co-expression of bi-directional genes. We also proposed that a special set of motifs in the bi-directional promoters play a role in transcriptional regulation of bi-directional genes. Our data also provide the putative regulatory motifs for experimental studies to investigate how the expression of divergent gene pairs is regulated.


Identification of bi-directional and uni-directional genes in human genome

Human genome assembly GRCh37, released as NCBI Build36 and Ensemble release 55, was downloaded from Genome Reference Consortium ( Gene annotation (NCBI Build36) was retrieved from the NCBI Entrez Gene ftp site ( The transcript annotation including transcription orientation, strand, starting site (hg19) was downloaded from hg19 RefGene table from UCSC Genome Browser ( A total of 45,408 genes (excluding mitochondrial genome) and 31,357 transcripts were collected and filtered for redundancy. This resulted in 44,293 non-redundant items of RefSeqs transcripts. Genes without clear mRNA information (NR, XR and XM) were filtered to ensure the exact transcription of all the genes. The 28520 mRNAs were collapsed into 21757 unique and non-overlapping clusters, which were further ranked according to their chromosome position and TSS coordinates to determine the adjacent gene pairs. Discrimination of bi-directional gene pairs and uni-directional genes was performed by a perl script according to the definition by Trinklein. et al[2]. Redundant gene pair entries that share the same intergenic sequence were removed.

Extraction of bi-directional promoter region

Based on the mapping information of gene and its transcripts, possible multiple TSSs were assessed. The intergenic regions between bi-directional genes’ TSS were taken as bi-directional promoters. For uni-directional genes the region of 1000 bp upstream of the TSS were extracted as promoter. Promoter regions were extracted from the chromosome fasta files of the latest GRCh37 version genome assembly datasets. (

Analysis of Promoter Sequences

The intergenic sequences of bidirectional genes were extended in both sides symmetrically into 1000 bp to meet the definition of a CpG island length. CpG island finder script[22] was run with two types of parameter criteria, %GC>=50, Obs/Exp>=0.60, length 500 and %GC >=55, Obs/Exp>=0.60, length 500 respectively. CpG frequency within the bi-directional and uni-bidirectional promoters was calculated.

Evaluation of Functional Enrichment

We utilized Gene Ontology (GO) categories ( and a commercial software MetaCore-GeneGO Pathway Maps ( to group functionally related genes and to contrast the functional distribution of bi-directional genes to the average distribution in the whole genome. The analysis of over-represented GO terms for bi-directional genes was performed by the GOEAST[23]. Statistical enrichment of a category was quantified using the Hypergeometric test method. Yekutieli multi-test adjustment method was applied to correct for multiple testing.

Genes were then mapped to GeneGO database by MetaCore™ tools to infer pathways preferentially targeted by bi-directional genes. In MetaCore™, the statistical significance of the enriched pathways is indicated by a P value yielded from the Fisher’s exact test. The False discovery rate (FDR) is also applied to correct for multiple testing.

Discovery of over-represented transcription factor binding sites

Putative TFBS in promoter regions were searched for matches to the position-weight matrix(PWM) in the JASPAR[24, 25] and TRANSFAC[26] database. Predetermined PWMs for 73 and 87 vertebrate TFBSs were extracted from TRANSFAC(public version 7.0) and JASPAR PSSM, respectively. Alignment of PWMs on genomic sequence was performed with COTRASIF[27] ( TFBSs within bi-directional promoters were categorized as over-represented, shared or under-represented at 2-fold threshold. Over-represented TFBS was defined as whose normalized number of binding sites in bidirectional promoters is 2-fold larger than those in unidirectional ones while under-represented means the normalized number of binding sites in bidirectional promoters is 2-fold smaller than the number of sites in a single unidirectional promoter. Shared motif is the intermediate state. A total of 18840 uni-directional promoters was used to give a contrast of bi-directional genes.


  1. Adachi N, Lieber MR: Bidirectional gene organization: a common architectural feature of the human genome. Cell. 2002, 109 (7): 807-809. 10.1016/S0092-8674(02)00758-4.

    Article  CAS  PubMed  Google Scholar 

  2. Trinklein ND, Aldred SF, Hartman SJ, Schroeder DI, Otillar RP, Myers RM: An abundance of bidirectional promoters in the human genome. Genome Res. 2004, 14 (1): 62-66. 10.1101/gr.1982804.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Piontkivska H, Yang MQ, Larkin DM, Lewin HA, Reecy J, Elnitski L: Cross-species mapping of bidirectional promoters enables prediction of unannotated 5' UTRs and identification of species-specific transcripts. BMC Genomics. 2009, 10: 189-10.1186/1471-2164-10-189.

    Article  PubMed Central  PubMed  Google Scholar 

  4. De Souza Santos E, De Bessa SA, Netto MM, Nagai MA: Silencing of LRRC49 and THAP10 genes by bidirectional promoter hypermethylation is a frequent event in breast cancer. Int J Oncol. 2008, 33 (1): 25-31.

    CAS  PubMed  Google Scholar 

  5. Lennard AC, Fried M: The bidirectional promoter of the divergently transcribed mouse Surf-1 and Surf-2 genes. Mol Cell Biol. 1991, 11 (3): 1281-1294.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Pollner R, Schmidt C, Fischer G, Kuhn K, Poschl E: Cooperative and competitive interactions of regulatory elements are involved in the control of divergent transcription of human Col4A1 and Col4A2 genes. FEBS Lett. 1997, 405 (1): 31-36. 10.1016/S0014-5793(97)00152-X.

    Article  CAS  PubMed  Google Scholar 

  7. Chen PY, Chang WS, Chou RH, Lai YK, Lin SC, Chi CY, Wu CW: Two non-homologous brain diseases-related genes, SERPINI1 and PDCD10, are tightly linked by an asymmetric bidirectional promoter in an evolutionarily conserved manner. BMC Mol Biol. 2007, 8: 2-10.1186/1471-2199-8-2.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Voth H, Oberthuer A, Simon T, Kahlert Y, Berthold F, Fischer M: Co-regulated expression of HAND2 and DEIN by a bidirectional promoter with asymmetrical activity in neuroblastoma. BMC Mol Biol. 2009, 10: 28-10.1186/1471-2199-10-28.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Lennard A, Gaston K, Fried M: The Surf-1 and Surf-2 genes and their essential bidirectional promoter elements are conserved between mouse and human. DNA Cell Biol. 1994, 13 (11): 1117-1126. 10.1089/dna.1994.13.1117.

    Article  CAS  PubMed  Google Scholar 

  10. Schuettengruber B, Doetzlhofer A, Kroboth K, Wintersberger E, Seiser C: Alternate activation of two divergently transcribed mouse genes from a bidirectional promoter is linked to changes in histone modification. J Biol Chem. 2003, 278 (3): 1784-1793. 10.1074/jbc.M204843200.

    Article  CAS  PubMed  Google Scholar 

  11. Williams EJ, Bowles DJ: Coexpression of neighboring genes in the genome of Arabidopsis thaliana. Genome Res. 2004, 14 (6): 1060-1067. 10.1101/gr.2131104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Kruglyak S, Tang H: Regulation of adjacent yeast genes. Trends Genet. 2000, 16 (3): 109-111. 10.1016/S0168-9525(99)01941-1.

    Article  CAS  PubMed  Google Scholar 

  13. Takai D, Jones PA: Origins of bidirectional promoters: computational analyses of intergenic distance in the human genome. Mol Biol Evol. 2004, 21 (3): 463-467. 10.1093/molbev/msh040.

    Article  CAS  PubMed  Google Scholar 

  14. Bellizzi D, Dato S, Cavalcante P, Covello G, Di Cianni F, Passarino G, Rose G, De Benedictis G: Characterization of a bidirectional promoter shared between two human genes related to aging: SIRT3 and PSMD13. Genomics. 2007, 89 (1): 143-150. 10.1016/j.ygeno.2006.09.004.

    Article  CAS  PubMed  Google Scholar 

  15. Gavalas A, Zalkin H: Analysis of the chicken GPAT/AIRC bidirectional promoter for de novo purine nucleotide synthesis. J Biol Chem. 1995, 270 (5): 2403-2410. 10.1074/jbc.270.5.2403.

    Article  CAS  PubMed  Google Scholar 

  16. Carninci P, Sandelin A, Lenhard B, Katayama S, Shimokawa K, Ponjavic J, Semple CA, Taylor MS, Engstrom PG, Frith MC, Forrest AR, Alkema WB, Tan SL, Plessy C, Kodzius R, Ravasi T, Kasukawa T, Fukuda S, Kanamori-Katayama M, Kitazume Y, Kawaji H, Kai C, Nakamura M, Konno H, Nakano K, Mottagui-Tabar S, Arner P, Chesi A, Gustincich S, Persichetti F, Suzuki H, Grimmond SM, Wells CA, Orlando V, Wahlestedt C, Liu ET, Harbers M, Kawai J, Bajic VB, Hume DA, Hayashizaki Y: Genome-wide analysis of mammalian promoter architecture and evolution. Nat Genet. 2006, 38 (6): 626-635. 10.1038/ng1789.

    Article  CAS  PubMed  Google Scholar 

  17. Gardiner-Garden M, Frommer M: CpG islands in vertebrate genomes. J Mol Biol. 1987, 196 (2): 261-282. 10.1016/0022-2836(87)90689-9.

    Article  CAS  PubMed  Google Scholar 

  18. Takai D, Jones PA: Comprehensive analysis of CpG islands in human chromosomes 21 and 22. Proc Natl Acad Sci U S A. 2002, 99 (6): 3740-3745. 10.1073/pnas.052410099.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Patton J, Block S, Coombs C, Martin ME: Identification of functional elements in the murine Gabp alpha/ATP synthase coupling factor 6 bi-directional promoter. Gene. 2006, 369: 35-44. 10.1016/j.gene.2005.10.009.

    Article  CAS  PubMed  Google Scholar 

  20. Kowalik TF, DeGregori J, Leone G, Jakoi L, Nevins JR: E2F1-specific induction of apoptosis and p53 accumulation, which is blocked by Mdm2. Cell Growth Differ. 1998, 9 (2): 113-118.

    CAS  PubMed  Google Scholar 

  21. Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M: Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals. Nature. 2005, 434 (7031): 338-345. 10.1038/nature03441.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Sujuan Y, Asaithambi A, Liu Y: CpGIF: an algorithm for the identification of CpG islands. Bioinformation. 2008, 2 (8): 335-338.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Burbelo PD, Martin GR, Yamada Y: Alpha 1(IV) and alpha 2(IV) collagen genes are regulated by a bidirectional promoter and a shared enhancer. Proc Natl Acad Sci U S A. 1988, 85 (24): 9679-9682. 10.1073/pnas.85.24.9679.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B: JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004, 32 (Database issue): D91-94. 10.1093/nar/gkh012.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Vlieghe D, Sandelin A, De Bleser PJ, Vleminckx K, Wasserman WW, van Roy F, Lenhard B: A new generation of JASPAR, the open-access repository for transcription factor binding site profiles. Nucleic Acids Res. 2006, 34 (Database issue): D95-97. 10.1093/nar/gkj115.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, Reuter I, Chekmenev D, Krull M, Hornischer K, Voss N, Stegmaier P, Lewicki-Potapov B, Saxel H, Kel AE, Wingender E: TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006, 34 (Database issue): D108-110. 10.1093/nar/gkj143.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Tokovenko B, Golda R, Protas O, Obolenskaya M, El'skaya A: COTRASIF: conservation-aided transcription-factor-binding site finder. Nucleic Acids Res. 2009, 37 (7): e49-10.1093/nar/gkp084.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Liu B, Chen J, Shen B: Computationally Comparative Analysis of CpG islands in Human Bidirectional and Unidirectional Promoters. IEEE Proc. 2010,

    Google Scholar 

  29. Huang CC, Chang WS: Cooperation between NRF-2 and YY-1 transcription factors is essential for triggering the expression of the PREPL-C2ORF34 bidirectional gene pair. BMC Mol Biol. 2009, 10: 67-10.1186/1471-2199-10-67.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Zanotto E, Hakkinen A, Teku G, Shen B, Ribeiro AS, Jacobs HT: NF-Y influences directionality of transcription from the bidirectional Mrps12/Sarsm promoter in both mouse and human cells. Biochim Biophys Acta. 2009, 1789 (5): 432-442.

    Article  CAS  PubMed  Google Scholar 

  31. Lin JM, Collins PJ, Trinklein ND, Fu Y, Xi H, Myers RM, Weng Z: Transcription factor binding and modified histones in human bidirectional promoters. Genome Res. 2007, 17 (6): 818-827. 10.1101/gr.5623407.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Ikeda S, Ayabe H, Mori K, Seki Y, Seki S: Identification of the functional elements in the bidirectional promoter of the mouse O-sialoglycoprotein endopeptidase and APEX nuclease genes. Biochem Biophys Res Commun. 2002, 296 (4): 785-791. 10.1016/S0006-291X(02)00939-7.

    Article  CAS  PubMed  Google Scholar 

  33. Zhang Y, Koushik S, Dai R, Mivechi NF: Structural organization and promoter analysis of murine heat shock transcription factor-1 gene. J Biol Chem. 1998, 273 (49): 32514-32521. 10.1074/jbc.273.49.32514.

    Article  CAS  PubMed  Google Scholar 

  34. Byrd PJ, Cooper PR, Stankovic T, Kullar HS, Watts GD, Robinson PJ, Taylor MR: A gene transcribed from the bidirectional ATM promoter coding for a serine rich protein: amino acid sequence, structure and expression studies. Hum Mol Genet. 1996, 5 (11): 1785-1791. 10.1093/hmg/5.11.1785.

    Article  CAS  PubMed  Google Scholar 

  35. Xu CF, Brown MA, Nicolai H, Chambers JA, Griffiths BL, Solomon E: Isolation and characterisation of the NBR2 gene which lies head to head with the human BRCA1 gene. Hum Mol Genet. 1997, 6 (7): 1057-1062. 10.1093/hmg/6.7.1057.

    Article  CAS  PubMed  Google Scholar 

  36. Gavalas A, Dixon JE, Brayton KA, Zalkin H: Coexpression of two closely linked avian genes for purine nucleotide synthesis from a bidirectional promoter. Mol Cell Biol. 1993, 13 (8): 4784-4792.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Seki Y, Ikeda S, Kiyohara H, Ayabe H, Seki T, Matsui H: Sequencing analysis of a putative human O-sialoglycoprotein endopeptidase gene (OSGEP) and analysis of a bidirectional promoter between the OSGEP and APEX genes. Gene. 2002, 285 (1-2): 101-108. 10.1016/S0378-1119(02)00429-8.

    Article  CAS  PubMed  Google Scholar 

  38. Ikeda S, Mochizuki A, Sarker AH, Seki S: Identification of functional elements in the bidirectional promoter of the mouse Nthl1 and Tsc2 genes. Biochem Biophys Res Commun. 2000, 273 (3): 1063-1068. 10.1006/bbrc.2000.3071.

    Article  CAS  PubMed  Google Scholar 

  39. Otte DM, Schwaab U, Luers GH: The Pxmp2 and PoleI genes are linked by a bidirectional promoter in an evolutionary conserved fashion. Gene. 2003, 313: 119-126. 10.1016/S0378-1119(03)00667-X.

    Article  CAS  PubMed  Google Scholar 

Download references


This research work has been supported in part by the National 973 Program of China (No. 2007CB947002) and the National Nature Science Foundation of China (20872107)

This article has been published as part of BMC Systems Biology Volume 5 Supplement 1, 2011: Selected articles from the 4th International Conference on Computational Systems Biology (ISB 2010). The full contents of the supplement are available online at

Author information

Authors and Affiliations


Corresponding author

Correspondence to Bairong Shen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

BL carried out the promoter compositional analysis, participated in the TFBS categorization and drafted the manuscript. JC participated in the functional enrichment analysis, performed the statistical analysis and draft the manuscript. BS conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.

Bingchuan Liu, Jiajia Chen contributed equally to this work.

Electronic supplementary material

Additional file 1 TableS1:Identification and statistics of GO terms enriched with bi-directional genes.(XLS 20 KB)


Additional file 2 TableS2:Functional similarities for plus and minus genes within annotated bi-directional gene pairs.(XLS 782 KB)


Additional file 3 TableS3:The enrichment of TFBS in bi-directional promoters relative to uni-directional promoters.(XLS 31 KB)

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), 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

Liu, B., Chen, J. & Shen, B. Genome-wide analysis of the transcription factor binding preference of human bi-directional promoters and functional annotation of related gene pairs. BMC Syst Biol 5 (Suppl 1), S2 (2011).

Download citation

  • Published:

  • DOI: