Cross-species mapping of bidirectional promoters enables prediction of unannotated 5' UTRs and identification of species-specific transcripts

Genetic variability of the activity of bidirectional promoters: a pilot study in bovine muscle

Bidirectional promoters are regulatory regions co-regulating the expression of two neighbouring genes organized in a head-to-head orientation. In recent years, these regulatory regions have been studied in many organisms; however, no investigation to date has been done to analyse the genetic variation of the activity of this type of promoter regions. In our study, we conducted an investigation to first identify bidirectional promoters sharing genes expressed in bovine Longissimus thoracis and then to find genetic variants affecting the activity of some of these bidirectional promoters. Combining bovine gene information and expression data obtained using RNA-Seq, we identified 120 putative bidirectional promoters active in bovine muscle. We experimentally validated in vitro 16 of these bidirectional promoters. Finally, using gene expression and whole-genome genotyping data, we explored the variability of the activity in muscle of the identified bidirectional promoters and discovered genetic variants affecting their activity. We found that the expression level of 77 genes is correlated with the activity of 12 bidirectional promoters. We also identified 57 single nucleotide polymorphisms associated with the activity of 5 bidirectional promoters. To our knowledge, our study is the first analysis in any species of the genetic variability of the activity of bidirectional promoters.

Chorion genes: a landscape of their evolution, structure, and regulation

Differential regulation at the level of transcription provides a means for controlling gene expression in eukaryotes, especially during development. Insect model systems have been extensively used to decipher the molecular basis of such regulatory cascades, and one of the oldest such model systems is the regulation of chorion gene expression during ovarian follicle maturation. Recent experimental and technological advances have shed new light onto the system, allowing us to revisit it. Thus, in this review we try to summarize almost 40 years' worth of studies on chorion gene regulation while-by comparing Bombyx mori and Drosophila melanogaster models-attempting to present a comprehensive, unified model of the various regulatory aspects of choriogenesis that takes into account the evolutionary conservation and divergence of the underlying mechanisms.

Bidirectional promoters in the transcription of mammalian genomes

In the genomes of humans and other mammals a large number of closely spaced pairs of genes that are transcribed in opposite directions were revealed. Their transcription is directed by so-called bidirectional promoters. This review is devoted to the characteristics of bidirectional promoters and features of their structure. The composition of "core" promoter elements in conventional unidirectional and bidirectional promoters is compared. Data on binding sites of transcription factors that are primarily specific for bidirectional promoters are discussed. The examples of promoters that share protein-coding genes transcribed by RNA polymerase II and the non-coding RNA genes transcribed by RNA polymerase III are described. Data obtained from global transcriptome analysis about the existence of short noncoding antisense RNA associated with the promoters in the context of the hypothesis of bidirectional transcription initiation as an inherent property of eukaryotic promoters are discussed.

Volatile evolution of long noncoding RNA repertoires: mechanisms and biological implications

Thousands of genes encoding long noncoding RNAs (lncRNAs) have been identified in all vertebrate genomes thus far examined. The list of lncRNAs partaking in arguably important biochemical, cellular, and developmental activities is steadily growing. However, it is increasingly clear that lncRNA repertoires are subject to weak functional constraint and rapid turnover during vertebrate evolution. We discuss here some of the factors that may explain this apparent paradox, including relaxed constraint on sequence to maintain lncRNA structure/function, extensive redundancy in the regulatory circuits in which lncRNAs act, as well as adaptive and non-adaptive forces such as genetic drift. We explore the molecular mechanisms promoting the birth and rapid evolution of lncRNA genes, with an emphasis on the influence of bidirectional transcription and transposable elements, two pervasive features of vertebrate genomes. Together these properties reveal a remarkably dynamic and malleable noncoding transcriptome which may represent an important source of robustness and evolvability.

Construction and validation of a novel dual reporter vector for studying mammalian bidirectional promoters

Regulation of gene expression plays important role in cellular functions. With the development of sequencing techniques, more and more genomes are available and genome-wide analyses of genomic structures that may affect gene expression regulation are now possible. Analyses of several genomes have found a class of regulatory regions that contain elements that initiate transcription of two different genes positioned with a head-to-head arrangement in two opposite directions. These regulatory regions are known as bidirectional promoters. Although bidirectional promoters have been known for years, recent genome-scale studies have shown that the regulation of the expression of up to 10% of the genes are controlled by bidirectional promoters. These findings are based mostly on computational work and only a limited number of putative bidirectional promoters have been experimentally validated. Developing methods to study bidirectional promoters will allow researchers to understand how these regions are regulated and the roles that divergent transcription plays in the expression of genes. Here, we have developed a novel dual-fluorescence reporter gene vector to study the transcriptional output of mammalian bidirectional promoters. We demonstrate that this vector is capable of expressing reporter genes under the control of bidirectional promoters, using the known human OSGEP/APEX bidirectional promoter.

Divergent transcription: a driving force for new gene origination?

The mammalian genome is extensively transcribed, a large fraction of which is divergent transcription from promoters and enhancers that is tightly coupled with active gene transcription. Here, we propose that divergent transcription may shape the evolution of the genome by new gene origination.

Bidirectional promoters as important drivers for the emergence of species-specific transcripts

The diversification of gene functions has been largely attributed to the process of gene duplication. Novel examples of genes originating from previously untranscribed regions have been recently described without regard to a unifying functional mechanism for their emergence. Here we propose a model mechanism that could generate a large number of lineage-specific novel transcripts in vertebrates through the activation of bidirectional transcription from unidirectional promoters. We examined this model in silico using human transcriptomic and genomic data and identified evidence consistent with the emergence of more than 1,000 primate-specific transcripts. These are transcripts with low coding potential and virtually no functional annotation. They initiate at less than 1 kb upstream of an oppositely transcribed conserved protein coding gene, in agreement with the generally accepted definition of bidirectional promoters. We found that the genomic regions upstream of ancestral promoters, where the novel transcripts in our dataset reside, are characterized by preferential accumulation of transposable elements. This enhances the sequence diversity of regions located upstream of ancestral promoters, further highlighting their evolutionary importance for the emergence of transcriptional novelties. By applying a newly developed test for positive selection to transposable element-derived fragments in our set of novel transcripts, we found evidence of adaptive evolution in the human lineage in nearly 3% of the novel transcripts in our dataset. These findings indicate that at least some novel transcripts could become functionally relevant, and thus highlight the evolutionary importance of promoters, through their capacity for bidirectional transcription, for the emergence of novel genes.

Identification of regulatory regions of bidirectional genes in cervical cancer

BACKGROUND

Bidirectional promoters are shared promoter sequences between divergent gene pair (genes proximal to each other on opposite strands), and can regulate the genes in both directions. In the human genome, > 10% of protein-coding genes are arranged head-to-head on opposite strands, with transcription start sites that are separated by < 1,000 base pairs. Many transcription factor binding sites occur in the bidirectional promoters that influence the expression of 2 opposite genes. Recently, RNA polymerase II (RPol II) ChIP-seq data are used to identify the promoters of coding genes and non-coding RNAs. However, a bidirectional promoter with RPol II ChIP-Seq data has not been found.

RESULTS

In some bidirectional promoter regions, the RPol II forms a bi-peak shape, which indicates that 2 promoters are located in the bidirectional region. We have developed a computational approach to identify the regulatory regions of all divergent gene pairs using genome-wide RPol II binding patterns derived from ChIP-seq data, based upon the assumption that the distribution of RPol II binding patterns around the bidirectional promoters are accumulated by RPol II binding of 2 promoters. In HeLa S3 cells, 249 promoter pairs and 1094 single promoters were identified, of which 76 promoters cover only positive genes, 86 promoters cover only negative genes, and 932 promoters cover 2 genes. Gene expression levels and STAT1 binding sites for different promoter categories were therefore examined.

CONCLUSIONS

The regulatory region of bidirectional promoter identification based upon RPol II binding patterns provides important temporal and spatial measurements regarding the initiation of transcription. From gene expression and transcription factor binding site analysis, the promoters in bidirectional regions may regulate the closest gene, and STAT1 is involved in primary promoter.

Human PSENEN and U2AF1L4 genes are concertedly regulated by a genuine bidirectional promoter

Head-to-head genes with a short distance between their transcription start sites may constitute up to 10% of all genes in the genomes of various species. It was hypothesized that this intergenic space may represent bidirectional promoters which are able to initiate transcription of both genes, but the true bidirectionality was proved only for a few of them. We present experimental evidence that, according to several criteria, a 269 bp region located between the PSENEN and U2AF1L4 human genes is a genuine bidirectional promoter regulating a concerted divergent transcription of these genes. Concerted transcription of PSENEN and U2AF1L4 can be necessary for regulation of T-cell activity.

RhesusBase: a knowledgebase for the monkey research community

Although the rhesus macaque is a unique model for the translational study of human diseases, currently its use in biomedical research is still in its infant stage due to error-prone gene structures and limited annotations. Here, we present RhesusBase for the monkey research community (http://www.rhesusbase.org). We performed strand-specific RNA-Seq studies in 10 macaque tissues and generated 1.2 billion 90-bp paired-end reads, covering >97.4% of the putative exon in macaque transcripts annotated by Ensembl. We found that at least 28.7% of the macaque transcripts were previously mis-annotated, mainly due to incorrect exon–intron boundaries, incomplete untranslated regions (UTRs) and missed exons. Compared with the previous gene models, the revised transcripts show clearer sequence motifs near splicing junctions and the end of UTRs, as well as cleaner patterns of exon–intron distribution for expression tags and cross-species conservation scores. Strikingly, 1292 exon–intron boundary revisions between coding exons corrected the previously mis-annotated open reading frames. The revised gene models were experimentally verified in randomly selected cases. We further integrated functional genomics annotations from >60 categories of public and in-house resources and developed an online accessible database. User-friendly interfaces were developed to update, retrieve, visualize and download the RhesusBase meta-data, providing a ‘one-stop’ resource for the monkey research community.

Modeling bidirectional transcription using silkmoth chorion gene promoters

Bidirectional transcription is an interesting feature of eukaryotic genomes; yet not all aspects of its mechanism are understood. Silkmoth choriogenesis is a model system for studying transcriptional regulation at the initiation level. As chorion genes comprise a large group of divergently transcribed gene pairs, we are presented with the possibility of investigating the intricacies of bidirectional transcription. Their well characterized 5' regulatory regions and expression profiles lay the foundation for investigating protein:protein and protein:DNA interactions, and RNA polymerase function during oocyte development. In this article we summarize current knowledge on chorion gene regulation and propose an approach to modeling bidirectional transcription using chorion promoters.

CNS effects of CB2 cannabinoid receptors: beyond neuro-immuno-cannabinoid activity

There are two well characterized cannabinoid receptors (CBRs), CB1-Rs and CB2-Rs, with other candidates, such as GPR55, PPARs and vanilloid TRPV1 (VR1) receptors, which are either activated by cannabinoids and/or endocannabinoids (eCBs). The neuronal and functional expression of CB2-Rs in the brain has been much less well characterized in comparison with the expression of the ubiquitous CB1-Rs. CB2-Rs were previously thought to be predominantly expressed in immune cells in the periphery and were traditionally referred to as peripheral CB2-Rs. We and others have now demonstrated the expression of CB2-Rs in neuronal, glial and endothelial cells in the brain, and this warrants a re-evaluation of the CNS effects of CB2-Rs. In the present review we summarize our current understanding of CNR2 genomic structure, its polymorphic nature, subtype specificity, from mice to human subjects, and its variants that confer vulnerabilities to neuropsychiatric disorders beyond neuro-immuno-cannabinoid activity.

NF-κB activity in perinatal brain during infectious and hypoxic-ischemic insults revealed by a reporter mouse

Infants suffering from infection or hypoxia-ischemia around the time of birth can develop brain damage resulting in life-long impairment such as cerebral palsy, epilepsy and cognitive disability. Inflammation appears to be an important contributor irrespective of whether the primary event is infection or hypoxia-ischemia. Activation of the transcription factor NF-κB is a hallmark of inflammation. To study perinatal brain inflammation, we developed a transgenic reporter mouse for imaging NF-κB activity in live animals and tissue samples. The reporter genes firefly luciferase and a destabilized version of enhanced GFP (dEGFP) were regulated by common NF-κB sites using a bidirectional promoter. Luciferase activity was imaged in vivo, while dEGFP was detected at cellular level in tissue sections. In newborn mice subjected to experimental models of infections or hypoxia-ischemia; luciferase signal increased in brains of live animals. In brain sections dEGFP expression, revealing NF-κB activation was observed in the endothelial cells of the blood-brain barrier in all disease models. In meningitis and hypoxia-ischemia expression of dEGFP was also induced in perivascular astrocytes. In conclusion, by using this transgenic reporter mouse in experimental models of perinatal complications, we could assess NF-κB activity in vivo and subsequently determine the cellular origin in the tissues.

A bidirectional promoter reporter vector for the analysis of the p53/WDR79 dual regulatory element

Analysis of numerous genomes has identified a class of regulatory regions that contain a head-to-head arrangement (5' to 5') on opposite strands of DNA. Often these regulatory regions have fewer than 1000 base pairs separating their corresponding transcription start sites and have been termed as being "bidirectional". This bidirectional arrangement and the divergent gene pairs under the control of these regulatory regions appear to be a common feature within genomes. Establishing methods to study these bidirectional transcriptional promoters, and understanding how they are regulated will allow researchers to gain more insight into the roles that divergent transcription plays in the expression and maintenance of protein coding genes. Recently, the p53 tumor suppressor gene was shown to have a bidirectional gene partner, WDR79. The transcription start sites (TSSs) of human and murine p53 and WDR79 genes are separated by approximately 800 and 930bp, respectively, in a head-to-head fashion, and fit the criteria of what is designated to be a putative bidirectional regulatory region. However, further testing is needed to demonstrate that the region between these genes contains a functional bidirectional promoter. Here, we have developed a bidirectional reporter vector, termed pLucRLuc, to study the transcriptional output of each promoter. This bidirectional reporter vector will allow researchers to determine the output of transcripts mediated by the bidirectional promoters. By focusing our studies on the transcriptional regulation of p53 and its bidirectional gene partner, WDR79, we hope to elucidate key factors that can control and regulate the expression of the p53 and WDR79 genes. Here, we demonstrate that pLucRLuc is a vector capable of expressing reporter genes under the control of bidirectional promoters in multiple human and murine cell lines and that the regulatory region upstream of the p53 and WDR79 TSSs is a bidirectional promoter controlled by common regulatory factors.

Efficient expression from one CMV enhancer controlling two core promoters

The human CMV promoter/enhancer is one of the strongest promoters for recombinant protein expression in mammalian cells, making the promoter very popular for production of recombinant antibodies. We used an antibody vector design where the antibody heavy and light chain genes were transcribed from a promoter complex consisting of two promoters arranged divergently with the 5' ends of the promoters in close proximity. However, when two identical CMV promoters constituted this promoter complex, the antibody expression observed was lower than expected based on the strength of the individual promoters. To optimize expression we prepared truncated promoter complexes where only one CMV enhancer controlled the initiation of transcription from two divergent minimal CMV core promoters. Antibody expression from the truncated promoter complexes was analyzed both when transiently transfected and upon stable site-specific integration into a CHO DG44 derived cell line. The data showed that it was possible for one enhancer to drive the expression of two core promoters. However, efficient expression from both divergent core promoters was seen only when the unique region upstream of the CMV enhancer was removed. Notably, a 12-fold increase in expression was found from the best of the truncated promoter complexes after stable site-specific integration when compared to the full-length double CMV promoter complex.

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

Background

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.

Results

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.

Conclusions

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.

eGOB: eukaryotic Gene Order Browser

A large number of genomes have been sequenced, allowing a range of comparative studies. Here, we present the eukaryotic Gene Order Browser with information on the order of protein and non-coding RNA (ncRNA) genes of 74 different eukaryotic species. The browser is able to display a gene of interest together with its genomic context in all species where that gene is present. Thereby, questions related to the evolution of gene organization and non-random gene order may be examined. The browser also provides access to data collected on pairs of adjacent genes that are evolutionarily conserved.

AVAILABILITY

eGOB as well as underlying data are freely available at http://egob.biomedicine.gu.se

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CONTACT

tore.samuelsson@medkem.gu.se.

Analysis of gene order conservation in eukaryotes identifies transcriptionally and functionally linked genes

The order of genes in eukaryotes is not entirely random. Studies of gene order conservation are important to understand genome evolution and to reveal mechanisms why certain neighboring genes are more difficult to separate during evolution. Here, genome-wide gene order information was compiled for 64 species, representing a wide variety of eukaryotic phyla. This information is presented in a browser where gene order may be displayed and compared between species. Factors related to non-random gene order in eukaryotes were examined by considering pairs of neighboring genes. The evolutionary conservation of gene pairs was studied with respect to relative transcriptional direction, intergenic distance and functional relationship as inferred by gene ontology. The results show that among gene pairs that are conserved the divergently and co-directionally transcribed genes are much more common than those that are convergently transcribed. Furthermore, highly conserved pairs, in particular those of fungi, are characterized by a short intergenic distance. Finally, gene pairs of metazoa and fungi that are evolutionary conserved and that are divergently transcribed are much more likely to be related by function as compared to poorly conserved gene pairs. One example is the ribosomal protein gene pair L13/S16, which is unusual as it occurs both in fungi and alveolates. A specific functional relationship between these two proteins is also suggested by the fact that they are part of the same operon in both eubacteria and archaea. In conclusion, factors associated with non-random gene order in eukaryotes include relative gene orientation, intergenic distance and functional relationships. It seems likely that certain pairs of genes are conserved because the genes involved have a transcriptional and/or functional relationship. The results also indicate that studies of gene order conservation aid in identifying genes that are related in terms of transcriptional control.

Mouse Dfa is a repressor of TATA-box promoters and interacts with the Abt1 activator of basal transcription

Our study of the mouse Ate1 arginyltransferase, a component of the N-end rule pathway, has shown that Ate1 pre-mRNA is produced from a bidirectional promoter that also expresses, in the opposite direction, a previously uncharacterized gene (Hu, R. G., Brower, C. S., Wang, H., Davydov, I. V., Sheng, J., Zhou, J., Kwon, Y. T., and Varshavsky, A. (2006) J. Biol. Chem. 281, 32559-32573). In this work, we began analyzing this gene, termed Dfa (divergent from Ate1). Mouse Dfa was found to be transcribed from both the bidirectional P(Ate1/Dfa) promoter and other nearby promoters. The resulting transcripts are alternatively spliced, yielding a complex set of Dfa mRNAs that are present largely, although not exclusively, in the testis. A specific Dfa mRNA encodes, via its 3'-terminal exon, a 217-residue protein termed Dfa(A). Other Dfa mRNAs also contain this exon. Dfa(A) is sequelogous (similar in sequence) to a region of the human/mouse HTEX4 protein, whose physiological function is unknown. We produced an affinity-purified antibody to recombinant mouse Dfa(A) that detected a 35-kDa protein in the mouse testis and in several cell lines. Experiments in which RNA interference was used to down-regulate Dfa indicated that the 35-kDa protein was indeed Dfa(A). Furthermore, Dfa(A) was present in the interchromatin granule clusters and was also found to bind to the Ggnbp1 gametogenetin-binding protein-1 and to the Abt1 activator of basal transcription that interacts with the TATA-binding protein. Given these results, RNA interference was used to probe the influence of Dfa levels in luciferase reporter assays. We found that Dfa(A) acts as a repressor of TATA-box transcriptional promoters.