Comparative analyses of bidirectional promoters in vertebrates
© Yang et al; licensee BioMed Central Ltd. 2008
Published: 28 May 2008
Orthologous genes with deep phylogenetic histories are likely to retain similar regulatory features. In this report we utilize orthology assignments for pairs of genes co-regulated by bidirectional promoters to map the ancestral history of the promoter regions.
Our mapping of bidirectional promoters from humans to fish shows that many such promoters emerged after the divergence of chickens and fish. Furthermore, annotations of promoters in deep phylogenies enable detection of missing data or assembly problems present in higher vertebrates. The functional importance of bidirectional promoters is indicated by selective pressure to maintain the arrangement of genes regulated by the promoter over long evolutionary time spans. Characteristics unique to bidirectional promoters are further elucidated using a technique for unsupervised classification, known as ESPERR.
Results of these analyses will aid in our understanding of the evolution of bidirectional promoters, including whether the regulation of two genes evolved as a consequence of their proximity or if function dictated their co-regulation.
Bidirectional promoters are defined as the regulatory regions that are shared between two genes, when those two genes are transcribed away from one another . The genes are said to be in a head-to-head arrangement, with their Transcription Start Sites (TSSs) positioned nearby one another. By definition, the intergenic distance between these genes (i.e. the promoter length) can be no greater than 1000 bp . This distance is measured from the TSS of the gene on the left of the promoter to the TSS of the gene on the right of the promoter. Head-to-head genes are spaced at this distance more frequently than expected in the human genome , suggesting a regulatory theme in gene expression. We recently showed that the human genome contains more bidirectional promoters than previously recognized [3, 4]. Here we map the orthologous regions of bidirectional promoters in seven additional species.
Using the "chains and nets" data from the UCSC Human Genome Browser and the Liftover tool , we are able to use the identity of genes on each side of a bidirectional promoter to find the corresponding functional location in other species. The use of these orthologous genes, which were present in the last common ancestor to the species being compared, is important because bidirectional promoters themselves often do not show a strong signal for conserved sequences. This fact makes the assignment of ancestral relationships difficult in these regulatory regions. Because genes flank both sides of bidirectional promoters, they provide markers of the ancestral history. The presence of the same head-to-head genes over long evolutionary time spans facilitates the assignment of orthology at their intervening promoter regions. Our approach complements and extends the work of Li et al. 2006  who examined bidirectional promoters in multiple species, because we are able to explore the ancestral history of bidirectional promoters in eight species simultaneously using whole genome orthology information. Our method uses orthologous genes as anchors that can be traced across vertebrate genomes. We use this information to track the appearance of bidirectional promoters in vertebrate evolution, predict when gene annotations are missing in higher vertebrates and determine which gene functions regulated by these promoters are the oldest.
Results and discussion
Orthologous gene pairs identify ancestral patterns of gene regulation
Other sets of bidirectional promoters showed a lineage-specific history. For instance, a large group of mammal-specific genes was not present in chickens or fish. A smaller group was only present in primates. In contrast, genes that were present in all species except chimp were likely to be missing from chimp due to assembly problems. Nearly twenty pairs of genes were found only in the human genome.
Intergenic distance at bidirectional promoters
The second graph shows the cumulative percentage of bidirectional promoters mapped in human and a second species, where the human dataset is limited to a 1000 bp distance. The most complete annotations were found in the human-mouse comparison. This result is illustrated by the similar curves for the cumulative percentage of orthologous bidirectional promoters in mouse that fall within 1000 bp. Up to 80% of all human bidirectional promoters were identified in mouse at this similar distance. In comparison, 75% of the human promoters were present in chimp within 1000 bp. The high levels of orthology found in mouse and chimp suggest that the 1000 bp distance will capture similar gene sets in other species. Thus we predict that the gene annotations of chimp, rhesus and dog will improve to represent a minimum of 80% of the bidirectional promoters the human genome.
Evolutionary comparison of head-to-head and tail-to-tail gene pairs
For chicken datasets the head-to-head gene sets were found more frequently than the tail-to-tail sets at 1000 bp, indicating that tail-to-tail arrangements of genes had been allowed to change in both distance and arrangement more often than head-to-head genes. These results indicate that selective pressure acts more strongly over evolutionary time to keep head-to-head genes together at the 1000 bp distance compared to tail-to-tail genes.
The data from the fish genomes indicated that very long distances were necessary to capture a majority of the human gene pairs. Given the compact nature of the fish genomes, it is unlikely that many of these long distance associations are biologically relevant. However, the preservation of tightly associated genes indicated the presence of important regulatory or functional roles that cannot be disturbed.
Gene ontology associated with bidirectional promoter regulation
Functions associated with orthologous genes regulated by bidirectional promoters were examined for those conserved in all seven species, or in the four mammals. Sixty pairs of genes were conserved across all seven species. These genes were examined for functional classifications. Four groups emerged: intracellular membrane bound organelle, macromolecule metabolism, chaperone, and mitochondrion. The p-values on these groups ranged from 10E-3 to 10E-1, and remained statistically significant following Benjamini correction for false discovery rate (i.e. ~2.7E-1).
Genes that were conserved across the four mammals had a much larger range of functional activities. Of 342 pairs of genes, catalytic activity emerged as the most significant enrichment in any functional class (6.1E-4 after Benjamini correction). Thus bidirectional promoters are regulating many enzymes in mammalian genomes. In total, 58 functional classes were significantly enriched in this dataset compared to a random collection of genes. These data indicate that the regulatory domain of bidirectional promoters has expanded to encompass a much larger set of gene functions in mammals.
Training ESPERR to discriminate bidirectional promoters
Our previous work indicated that sequence-based characteristics were different in bidirectional promoters and non-bidirectional promoters . However the size of the datasets was quite disparate (1,005 bidirectional, 17,613 non-bidirectional). Therefore for training ESPERR  we sampled equal size subsets of 800 elements from each class (keeping the remaining elements in each class as test sets for verification). For each training interval we then extracted genomic alignments of six species (human, chimpanzee, macaque, mouse, rat, and dog) from the 17 species alignments available in the UCSC Genome Browser. Regions of the training data overlapping coding exons (from UCSC Known Genes) were masked out. We first performed an unsupervised encoding selection (the first stage of the ESPERR procedure) to create an encoding in 10 symbols. Leave-one-out cross validation on the training data using this encoding yielded a success rate of 76%. On the bidirectional test set, the model trained using this mapping correctly classified 404 elements (89%) and incorrectly classified 50 (17 elements were not included due to insufficient alignment). On the non-bidirectional test set it successfully classified 11,150 elements (70%), and incorrectly classified 4,845 (687 elements not included). Next, we performed the full ESPERR procedure, using the first stage reduction to produce an encoding of size 75, which was then refined with the heuristic search yielding an encoding of size 10. The resulting encoding gave a modest improvement in cross-validation, with a success rate of 82%. However, on the bidirectional test set, the model using this encoding classified 405 elements (89%) and incorrectly classified 49 (17 elements were again not included due to insufficient alignment). On the non-bidirectional test set it successfully classified 10,900 elements (68%), and incorrectly classified 5,095 (687 elements not included). Thus, using the ESPERR heuristic search gives no improvement for classifying this dataset. It is noteworthy that these classification rates, though modest, indicate that there are sequence and evolutionary patterns that can be captured to characterize bidirectional promoters. Particularly interesting is the substantially greater generalization rate for the bidirectional test set, suggesting that there are more characteristic signals for these elements that can be captured. This is consistent with the result of the ESPERR heuristic search – optimizing the encoding using the training data gives a slight improvement in recognizing the bidirectional test elements, but at the cost of poorer performance on the non-bidirectional test set.
Our study of bidirectional promoters across orthologous regions of eight species provides a foundation for optimized annotations of these regulatory regions in higher vertebrates, including chimp, rhesus and dog. Furthermore the functional analyses of genes regulated by these promoters show that a small subset of specialized functions in chickens and fish was expanded in mammals to include a wide breadth of activities. A common regulatory mechanism is likely to exist that coordinately regulates genes in these functional pathways. We continue to investigate the features associated with bidirectional promoters using a classification procedure containing supervised and unsupervised techniques. The results are promising in that they indicate that bidirectional promoters have features consistent with learnable patterns.
Assigning orthologous regions
A multi-stage approach to mapping orthology at bidirectional promoters was developed. Because orthology assignments are strongest in coding regions, we mapped single human genes from head-to-head gene pairs to a second species. To identify orthologous DNA we used "chains and nets" data from the UCSC Genome Browser mysql tables. Chains in the Genome Browser represent sequences of gapless aligned blocks. Nets provide a hierarchical ordering of those chains. Level 1 chains contain the longest, best scoring sequence chains that span any selected region. Gaps in the level 1 chains are recorded in the level 2 chains (of the Browser mysql tables). This ordering process is repeated until all aligned sequences are assigned to a homologous human region. Odd number levels represent aligned regions and even number levels correspond to gapped regions separating the best scoring chained alignments.
We used orthologous regions present in only level 1 and excluded any other levels, which contain both paralogous (duplicated during evolution) and orthologous sequences. Level 1 alignments contain extremely long stretches of genes in conserved synteny (i.e. same gene identity and location) between species. The regions of conserved sequence forming these alignment blocks are separated by gaps that provide spacers between them. Frequently the aligned regions correspond to exons and gapped regions correspond to introns and intergenic regions. Given a human gene, our approach examined whether it fell within an orthologous region defined by level 1 alignment data without knowledge of the exact position within an alignment or on which side of a gap. In a subsequent step, we intersected the positions of gaps and exons of each gene to (1) identify the orthologous gene in the second species and (2) to determine how well the exons align between species.
Mapping orthologous genes
After determining the orthology assignments using the UCSC alignment data, we used the Ensembl annotations  to search the identity of genes within each corresponding alignment. The appearance of several genes in the same region was handled by choosing the candidate with the closest transcription start site to its neighboring gene. This technique identified the most likely gene pair for regulation by a bidirectional promoter. Orthology assignments were checked for each human gene individually, and subsequently checked to see if the pairs from human also formed pairs in the other species. Once orthologous genes were identified for both human genes forming a pair, the orthology assignments were checked in the reverse direction from the other species to human.
Count of Ensembl gene annotations used to find orthologous bidirectional promoters
Orthology mapping for 1642 human genes regulated by bidirectional promoters
Number ofOrthologous Genes
Number of Orthologous Genes
Dealing with special cases
Some orthologous regions aligned perfectly over very long distances in the second species without any gaps. This situation occurred most frequently in close evolutionary comparisons, such as human to chimp or rhesus. The perfect alignments complicated our mapping of individual gene orthology, because there were no breaks to separate genes within the chained alignment region. To circumvent this problem, we used the Liftover tool available at the UCSC genome Browser. The Liftover software converted the genomic coordinates from human to the second species using the genome alignment information at the nucleotide level. This approach converted perfectly aligned genes between human (hg17) and chimp, rhesus, dog, mouse, chicken or zebrafish. Although Liftover data is not available for Fugu, it appears that only 11 such genes align perfectly to Fugu, which is a small enough number to handle manually.
Although a singular approach using the Liftover tool could replace the procedure of using chains and nets to more precisely identify each orthologous coordinate, we chose to use it as a second phase in the mapping process. In this way more information was retained from chains and nets regarding how well each exon aligned and how well all genes aligned when gapped regions were present. The overall procedure appears in Figure 1.
LE and MQY are funded by the Intramural Research Program of the National Human Genome Research Institute, US National Institutes of Health.
This article has been published as part of BMC Bioinformatics Volume 9 Supplement 6, 2008: Symposium of Computations in Bioinformatics and Bioscience (SCBB07). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/9?issue=S6.
- 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-4View ArticlePubMedGoogle Scholar
- 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.1982804PubMed CentralView ArticlePubMedGoogle Scholar
- Yang MQ, Elnitski L: A computational study of bidirectional promoters in the human genome. Springer Lecture Series: Notes in Bioinformatics 2007.Google Scholar
- Yang MQ, Koehly LM, Elnitski LL: Comprehensive annotation of human bidirectional promoters identifies co-regulation among breast and ovarian cancer genes. PLOS Comput Biol 2007, 20: e72. 10.1371/journal.pcbi.0030072View ArticleGoogle Scholar
- Kuhn RM, Karolchik D, Zweig AS, Trumbower H, Thomas DJ, Thakkapallayil A, Sugnet CW, Stanke M, Smith KE, Siepel A, et al.: The UCSC genome browser database: update 2007. Nucleic Acids Res 2007, (35 Database):D668–673. 10.1093/nar/gkl928Google Scholar
- Li YY, Yu H, Guo ZM, Guo TQ, Tu K, Li YX: Systematic analysis of head-to-head gene organization: evolutionary conservation and potential biological relevance. PLoS Comput Biol 2006,2(7):e74. 10.1371/journal.pcbi.0020074PubMed CentralView ArticlePubMedGoogle Scholar
- Yang MQ, Elnitski L: Orthology of Bidirectional Promoters Enables Use of a Multiple Class Predictor for Discriminating Functional Elements in the Human Genome. Proceedings of BIOCOMP 2007 2007, in press.Google Scholar
- Taylor J, Tyekucheva S, King DC, Hardison RC, Miller W, Chiaromonte F: ESPERR: learning strong and weak signals in genomic sequence alignments to identify functional elements. Genome Res 2006,16(12):1596–1604. 10.1101/gr.4537706PubMed CentralView ArticlePubMedGoogle Scholar
- Hubbard TJ, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, et al.: Ensembl 2007. Nucleic Acids Res 2007, (35 Database):D610–617. 10.1093/nar/gkl996Google Scholar
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.