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

A Bidirectional Non-Coding RNA Promoter Mediates Long-Range Gene Expression Regulation

Recent evidence suggests that human gene promoters display gene expression regulatory mechanisms beyond the typical single gene local transcription modulation. In mammalian genomes, genes with an associated bidirectional promoter are abundant; bidirectional promoter architecture serves as a regulatory hub for a gene pair expression. However, it has been suggested that its contribution to transcriptional regulation might exceed local transcription initiation modulation. Despite their abundance, the functional consequences of bidirectional promoter architecture remain largely unexplored. This work studies the long-range gene expression regulatory role of a long non-coding RNA gene promoter using chromosome conformation capture methods. We found that this particular bidirectional promoter contributes to distal gene expression regulation in a target-specific manner by establishing promoter-promoter interactions. In particular, we validated that the promoter-promoter interactions of this regulatory element with the promoter of distal gene BBX contribute to modulating the transcription rate of this gene; removing the bidirectional promoter from its genomic context leads to a rearrangement of BBX promoter-enhancer interactions and to increased gene expression. Moreover, long-range regulatory functionality is not directly dependent on its associated non-coding gene pair expression levels.

Landscape of Overlapping Gene Expression in the Equine Placenta

Increasing evidence suggests that overlapping genes are much more common in eukaryotic genomes than previously thought. These different-strand overlapping genes are potential sense–antisense (SAS) pairs, which might have regulatory effects on each other. In the present study, we identified the SAS loci in the equine genome using previously generated stranded, paired-end RNA sequencing data from the equine chorioallantois. We identified a total of 1261 overlapping loci. The ratio of the number of overlapping regions to chromosomal length was numerically higher on chromosome 11 followed by chromosomes 13 and 12. These results show that overlapping transcription is distributed throughout the equine genome, but that distributions differ for each chromosome. Next, we evaluated the expression patterns of SAS pairs during the course of gestation. The sense and antisense genes showed an overall positive correlation between the sense and antisense pairs. We further provide a list of SAS pairs with both positive and negative correlation in their expression patterns throughout gestation. This study characterizes the landscape of sense and antisense gene expression in the placenta for the first time and provides a resource that will enable researchers to elucidate the mechanisms of sense/antisense regulation during pregnancy.

Characterization of bidirectional gene pairs in The Cancer Genome Atlas (TCGA) dataset

The “bidirectional gene pair” indicates a particular head-to-head gene organization in which transcription start sites of two genes are located on opposite strands of genomic DNA within a region of one kb. Despite bidirectional gene pairs are well characterized, little is known about their expression profiles and regulation features in tumorigenesis. We used RNA-seq data from The Cancer Genome Atlas (TCGA) dataset for a systematic analysis of the expression profiles of bidirectional gene pairs in 13 cancer datasets. Gene pairs on the opposite strand with transcription end site distance within one kb or on the same strand with the distance of two genes between 1–10 kb and gene pairs comprising two randomly chosen genes were used as control gene pairs (CG1, CG2, and random). We identified and characterized up-/down-regulated genes by comparing the expression level between tumors and adjacent normal tissues in 13 TCGA datasets. There were no consistently significant difference in the percentage of up-/down-regulated genes between bidirectional and control/random genes in most of TCGA datasets. However, the percentage of bidirectional gene pairs comprising two up- or two down-regulated genes was significantly higher than gene pairs from CG1/2 in 12/11 analyzed TCGA datasets and the random gene pairs in all 13 TCGA datasets. Then we identified the methylation correlated bidirectional genes to explore the regulatory mechanism of bidirectional genes. Like the differentially expressed gene pairs, the bidirectional genes in a pair were significantly prone to be both hypo- or hyper-methylation correlated genes in 12/13 TCGA datasets when comparing to the CG2/random gene pairs despite no significant difference between the percentages of hypo-/hyper-methylation correlated genes in bidirectional and CG2/random genes in most of TCGA datasets. Finally, we explored the correlation between bidirectional genes and patient’s survival, identifying prognostic bidirectional genes and prognostic bidirectional gene pairs in each TCGA dataset. Remarkably, we found a group of prognostic bidirectional gene pairs in which the combination of two protein coding genes with different expression level correlated with different survival prognosis in survival analysis for OS. The percentage of these gene pairs in bidirectional gene pair were significantly higher than the gene pairs in controls in COAD datasets and lower in none of 13 TCGA datasets.

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.

Novel Biomarker MicroRNAs for Subtyping of Acute Coronary Syndrome: A Bioinformatics Approach

Acute coronary syndrome (ACS) is a life-threatening disease that affects more than half a million people in United States. We currently lack molecular biomarkers to distinguish the unstable angina (UA) and acute myocardial infarction (AMI), which are the two subtypes of ACS. MicroRNAs play significant roles in biological processes and serve as good candidates for biomarkers. In this work, we collected microRNA datasets from the Gene Expression Omnibus database and identified specific microRNAs in different subtypes and universal microRNAs in all subtypes based on our novel network-based bioinformatics approach. These microRNAs were studied for ACS association by pathway enrichment analysis of their target genes. AMI and UA were associated with 27 and 26 microRNAs, respectively, nine of them were detected for both AMI and UA, and five from each subtype had been reported previously. The remaining 22 and 21 microRNAs are novel microRNA biomarkers for AMI and UA, respectively. The findings are then supported by pathway enrichment analysis of the targets of these microRNAs. These novel microRNAs deserve further validation and will be helpful for personalized ACS diagnosis.

Positive Selection on a Regulatory Insertion-Deletion Polymorphism in FADS2 Influences Apparent Endogenous Synthesis of Arachidonic Acid

Long chain polyunsaturated fatty acids (LCPUFA) are bioactive components of membrane phospholipids and serve as substrates for signaling molecules. LCPUFA can be obtained directly from animal foods or synthesized endogenously from 18 carbon precursors via the FADS2 coded enzyme. Vegans rely almost exclusively on endogenous synthesis to generate LCPUFA and we hypothesized that an adaptive genetic polymorphism would confer advantage. The rs66698963 polymorphism, a 22-bp insertion–deletion within FADS2, is associated with basal FADS1 expression, and coordinated induction of FADS1 and FADS2 in vitro. Here, we determined rs66698963 genotype frequencies from 234 individuals of a primarily vegetarian Indian population and 311 individuals from the US. A much higher I/I genotype frequency was found in Indians (68%) than in the US (18%). Analysis using 1000 Genomes Project data confirmed our observation, revealing a global I/I genotype of 70% in South Asians, 53% in Africans, 29% in East Asians, and 17% in Europeans. Tests based on population divergence, site frequency spectrum, and long-range haplotype consistently point to positive selection encompassing rs66698963 in South Asian, African, and some East Asian populations. Basal plasma phospholipid arachidonic acid (ARA) status was 8% greater in I/I compared with D/D individuals. The biochemical pathway product–precursor difference, ARA minus linoleic acid, was 31% and 13% greater for I/I and I/D compared with D/D, respectively. This study is consistent with previous in vitro data suggesting that the insertion allele enhances n-6 LCPUFA synthesis and may confer an adaptive advantage in South Asians because of the traditional plant-based diet practice.

Characteristics of Antisense Transcript Promoters and the Regulation of Their Activity

Recently, an increasing number of studies on natural antisense transcripts have been reported, especially regarding their classification, temporal and spatial expression patterns, regulatory functions and mechanisms. It is well established that natural antisense transcripts are produced from the strand opposite to the strand encoding a protein. Despite the pivotal roles of natural antisense transcripts in regulating the expression of target genes, the transcriptional mechanisms initiated by antisense promoters (ASPs) remain unknown. To date, nearly all of the studies conducted on this topic have focused on the ASP of a single gene of interest, whereas no study has systematically analyzed the locations of ASPs in the genome, ASP activity, or factors influencing this activity. This review focuses on elaborating on and summarizing the characteristics of ASPs to extend our knowledge about the mechanisms of antisense transcript initiation.

The AVR2-SIX5 gene pair is required to activate I-2-mediated immunity in tomato

Plant-invading microbes betray their presence to a plant by exposure of antigenic molecules such as small, secreted proteins called 'effectors'. In Fusarium oxysporum f. sp. lycopersici (Fol) we identified a pair of effector gene candidates, AVR2-SIX5, whose expression is controlled by a shared promoter. The pathogenicity of AVR2 and SIX5 Fol knockouts was assessed on susceptible and resistant tomato (Solanum lycopersicum) plants carrying I-2. The I-2 NB-LRR protein confers resistance to Fol races carrying AVR2. Like Avr2, Six5 was found to be required for full virulence on susceptible plants. Unexpectedly, each knockout could breach I-2-mediated disease resistance. So whereas Avr2 is sufficient to induce I-2-mediated cell death, Avr2 and Six5 are both required for resistance. Avr2 and Six5 interact in yeast two-hybrid assays as well as in planta. Six5 and Avr2 accumulate in xylem sap of plants infected with the reciprocal knockouts, showing that lack of I-2 activation is not due to a lack of Avr2 accumulation in the SIX5 mutant. The effector repertoire of a pathogen determines its host specificity and its ability to manipulate plant immunity. Our findings challenge an oversimplified interpretation of the gene-for-gene model by showing requirement of two fungal genes for immunity conferred by one resistance gene.

The landscape of antisense gene expression in human cancers

High-throughput RNA sequencing has revealed more pervasive transcription of the human genome than previously anticipated. However, the extent of natural antisense transcripts’ (NATs) expression, their regulation of cognate sense genes, and the role of NATs in cancer remain poorly understood. Here, we use strand-specific paired-end RNA sequencing (ssRNA-seq) data from 376 cancer samples covering nine tissue types to comprehensively characterize the landscape of antisense expression. We found consistent antisense expression in at least 38% of annotated transcripts, which in general is positively correlated with sense gene expression. Investigation of sense/antisense pair expressions across tissue types revealed lineage-specific, ubiquitous and cancer-specific antisense loci transcription. Comparisons between tumor and normal samples identified both concordant (same direction) and discordant (opposite direction) sense/antisense expression patterns. Finally, we provide OncoNAT, a catalog of cancer-related genes with significant antisense transcription, which will enable future investigations of sense/antisense regulation in cancer. Using OncoNAT we identified several functional NATs, including NKX2-1-AS1 that regulates the NKX2-1 oncogene and cell proliferation in lung cancer cells. Overall, this study provides a comprehensive account of NATs and supports a role for NATs' regulation of tumor suppressors and oncogenes in cancer biology.

Regulation of Interferon-Stimulated Gene BST2 by a lncRNA Transcribed from a Shared Bidirectional Promoter

Recent genome-wide studies have revealed the presence of thousands of long non-protein-coding RNAs (lncRNAs), some of which may play critical roles in the cell. We have previously shown that a large number of lncRNAs show differential expression in response to interferon (IFN)α stimulation in primary human cells. Here, we show that a subset of IFN-induced lncRNAs are positioned in proximity of protein-coding IFN-stimulated genes (ISGs). The majority of gene pairs originated from bidirectional promoters and showed positively correlated expression. We focused our analysis on a pair consisting of the known protein-coding ISG, BST2, and an un-studied putative lncRNA originating from the promoter region of BST2 in a divergent orientation. We showed that this transcript was a multi-exonic, polyadenylated long RNA that lacked protein-coding capacity. BST2 and the lncRNA were both induced in response to IFNα in diverse cell types. The induction of both genes was mediated through the JAK-STAT pathway, suggesting that IFN-stimulated response elements within the shared promoter activated the transcription of both genes. RNAi-mediated knock-down of the lncRNA resulted in down-regulation of BST2, and we could show that this down-regulation occurred at the level of transcription. Forced overexpression of this lncRNA, which we named BST2 IFN-Stimulated Positive Regulator (BISPR), resulted in up-regulation of BST2, indicating that the regulation of expression of BST2 by BISPR is mediated through interactions involving BISPR RNA itself, rather than the impact of its transcription from an adjacent locus. Importantly, upon IFN stimulation, transcriptional activation of BISPR preceded the induction of BST2, suggesting that expression of BISPR facilitated the initiation of transcription in its paired protein-coding gene. The lncRNA-mediated transcriptional regulation described in this study may help govern the expression of additional protein-coding RNAs involved in IFN response and other cellular processes.

Deep sequencing reveals a novel class of bidirectional promoters associated with neuronal genes

Background

Comprehensive annotation of transcripts expressed in a given tissue is a critical step towards the understanding of regulatory and functional pathways that shape the transcriptome.

Results

Here, we reconstructed a cumulative transcriptome of the human prefrontal cortex (PFC) based on approximately 300 million strand-specific RNA sequence (RNA-seq) reads collected at different stages of postnatal development. We find that more than 50% of reconstructed transcripts represent novel transcriptome elements, including 8,343 novel exons and exon extensions of annotated coding genes, 11,217 novel antisense transcripts and 29,541 novel intergenic transcripts or their fragments showing canonical features of long non-coding RNAs (lncRNAs). Our analysis further led to a surprising discovery of a novel class of bidirectional promoters (NBiPs) driving divergent transcription of mRNA and novel lncRNA pairs and displaying a distinct set of sequence and epigenetic features. In contrast to known bidirectional and unidirectional promoters, NBiPs are strongly associated with genes involved in neuronal functions and regulated by neuron-associated transcription factors.

Conclusions

Taken together, our results demonstrate that large portions of the human transcriptome remain uncharacterized. The distinct sequence and epigenetic features of NBiPs, as well as their specific association with neuronal genes, further suggest existence of regulatory pathways specific to the human brain.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-457) contains supplementary material, which is available to authorized users.

Identification of microRNA as sepsis biomarker based on miRNAs regulatory network analysis

Sepsis is regarded as arising from an unusual systemic response to infection but the physiopathology of sepsis remains elusive. At present, sepsis is still a fatal condition with delayed diagnosis and a poor outcome. Many biomarkers have been reported in clinical application for patients with sepsis, and claimed to improve the diagnosis and treatment. Because of the difficulty in the interpreting of clinical features of sepsis, some biomarkers do not show high sensitivity and specificity. MicroRNAs (miRNAs) are small noncoding RNAs which pair the sites in mRNAs to regulate gene expression in eukaryotes. They play a key role in inflammatory response, and have been validated to be potential sepsis biomarker recently. In the present work, we apply a miRNA regulatory network based method to identify novel microRNA biomarkers associated with the early diagnosis of sepsis. By analyzing the miRNA expression profiles and the miRNA regulatory network, we obtained novel miRNAs associated with sepsis. Pathways analysis, disease ontology analysis, and protein-protein interaction network (PIN) analysis, as well as ROC curve, were exploited to testify the reliability of the predicted miRNAs. We finally identified 8 novel miRNAs which have the potential to be sepsis biomarkers.

Bioinformatic tools in Arabidopsis research

Bioinformatic tools are an increasingly important resource for Arabidopsis researchers. With them, it is possible to rapidly query the large data sets covering genomes, transcriptomes, proteomes, epigenomes, and other "omes" that have been generated in the past decade. Often these tools can be used to generate quality hypotheses at the click of a mouse. In this chapter, we cover the use of bioinformatic tools for examining gene expression and coexpression patterns, performing promoter analyses, looking for functional classification enrichment for sets of genes, and investigating protein-protein interactions. We also introduce bioinformatic tools that allow integration of data from several sources for improved hypothesis generation.

Bidirectional promoters in seed development and related hormone/stress responses

BACKGROUND

Bidirectional promoters are common in genomes but under-studied experimentally, particularly in plants. We describe a targeted identification and selection of a subset of putative bidirectional promoters to identify genes involved in seed development and to investigate possible coordinated responses of gene pairs to conditions important in seed maturation such as desiccation and ABA-regulation.

RESULTS

We combined a search for 100-600 bp intergenic regions in the Arabidopsis genome with a cis-element based selection for those containing multiple copies of the G-box motif, CACGTG. One of the putative bidirectional promoters identified also contained a CE3 coupling element 5 bp downstream of one G-box and is identical to that characterized previously in the HVA1 promoter of barley. CE3 elements are significantly under-represented and under-studied in Arabidopsis. We further characterized the pair of genes associated with this promoter and uncovered roles for two small, previously uncharacterized, plant-specific proteins in Arabidopsis seed development and stress responses.

CONCLUSIONS

Using bioinformatics we identified putative bidirectional promoters involved in seed development and analysed expression patterns for a pair of plant-specific genes in various tissues and in response to hormones/stress. We also present preliminary functional analysis of these genes that is suggestive of roles in seed development.

Functional networking of human divergently paired genes (DPGs)

Divergently paired genes (DPGs), also known as bidirectional (head-to-head positioned) genes, are conserved across species and lineages, and thus deemed to be exceptional in genomic organization and functional regulation. Despite previous investigations on the features of their conservation and gene organization, the functional relationship among DPGs in a given species and lineage has not been thoroughly clarified. Here we report a network-based comprehensive analysis on human DPGs and our results indicate that the two members of the DPGs tend to participate in different biological processes while enforcing related functions as modules. Comparing to randomly paired genes as a control, the DPG pairs have a tendency to be clustered in similar “cellular components” and involved in similar “molecular functions”. The functional network bridged by DPGs consists of three major modules. The largest module includes many house-keeping genes involved in core cellular activities. This module also shows low variation in expression in both CNS (central nervous system) and non-CNS tissues. Based on analyses of disease transcriptome data, we further suggest that this particular module may play crucial roles in HIV infection and its disease mechanism.

GSTCD and INTS12 regulation and expression in the human lung

Genome-Wide Association Study (GWAS) meta-analyses have identified a strong association signal for lung function, which maps to a region on 4q24 containing two oppositely transcribed genes: glutathione S-transferase, C-terminal domain containing (GSTCD) and integrator complex subunit 12 (INTS12). Both genes were found to be expressed in a range of human airway cell types. The promoter regions and transcription start sites were determined in mRNA from human lung and a novel splice variant was identified for each gene. We obtained the following evidence for GSTCD and INTS12 co-regulation and expression: (i) correlated mRNA expression was observed both via Q-PCR and in a lung expression quantitative trait loci (eQTL) study, (ii) induction of both GSTCD and INTS12 mRNA expression in human airway smooth muscle cells was seen in response to TGFβ1, (iii) a lung eQTL study revealed that both GSTCD and INTS12 mRNA levels positively correlate with percent predicted FEV1, and (iv) FEV1 GWAS associated SNPs in 4q24 were found to act as an eQTL for INTS12 in a number of tissues. In fixed sections of human lung tissue, GSTCD protein expression was ubiquitous, whereas INTS12 expression was predominantly in epithelial cells and pneumocytes. During human fetal lung development, GSTCD protein expression was observed to be highest at the earlier pseudoglandular stage (10-12 weeks) compared with the later canalicular stage (17-19 weeks), whereas INTS12 expression levels did not alter throughout these stages. Knowledge of the transcriptional and translational regulation and expression of GSTCD and INTS12 provides important insights into the potential role of these genes in determining lung function. Future work is warranted to fully define the functions of INTS12 and GSTCD.

LuFLA1PRO and LuBGAL1PRO promote gene expression in the phloem fibres of flax (Linum usitatissimum)

Cell type-specific promoters were identified that drive gene expression in an industrially important product. To identify flax (Linum usitatissimum) gene promoters, we analyzed the genomic regions upstream of a fasciclin-like arabinogalactan protein (LuFLA1) and a beta-galactosidase (LuBGAL1). Both of these genes encode transcripts that have been found to be highly enriched in tissues bearing phloem fibres. Using a beta-glucuronidase (GUS) reporter construct, we found that a 908-bp genomic sequence upstream of LuFLA1 (LuFLA1PRO) directed GUS expression with high specificity to phloem fibres undergoing secondary cell wall development. The DNA sequence upstream of LuBGAL1 (LuBGAL1PRO) likewise produced GUS staining in phloem fibres with developing secondary walls, as well as in tissues of developing flowers and seed bolls. These data provide further evidence of a specific role for LuFLA1 in phloem fibre development, and demonstrate the utility of LuFLA1PRO and LuBGAL1PRO as tools for biotechnology and further investigations of phloem fibre development.

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.

The preservation of bidirectional promoter architecture in eukaryotes: what is the driving force?

Background

The bidirectional gene architecture has been studied in many organisms, and the conservation of bidirectional arrangement has received considerable attention. However, the explanation for the evolutionary conservation about this genomic structure is still insufficient. In this study the large scale identification and pathway enrichment analysis for bidirectional genes were performed in several eukaryotes and the comparative analysis of this arrangement between human and mouse were dissected for the purpose of discovering the driving force of the preservation of this genomic structure.

Results

We identified the bidirectional gene pairs in eight different species and found this structure to be prevalent in eukaryotes. The pathway enrichment analysis indicated the bidirectional genes at the genome level are conserved in certain pathways, such as the DNA repair and some other fundamental cellular pathways. The comparative analysis about the gene expression, function, between human and mouse bidirectional genes were also performed and we observed that the selective force of this architecture doesn't derive from the co-regulation between paired genes, but the functional bias of bidirectional genes at whole genome level is observed strengthened during evolution.

Conclusions

Our result validated the coexpression of bidirectional genes; however failed to support their functional relevance. The conservation of bidirectional promoters seems not the result of functional connection between paired genes, but the functional bias at whole genome level, which imply that the genome-wide functional constraint is important for the conservation of bidirectional structure.

A study on promoter characteristics of head-to-head genes in Saccharomyces cerevisiae

Background

Head-to-head (h2h) genes are prone to have association in expression and in functionality and have been shown conserved in evolution. Currently there are many studies on such h2h gene pairs. We found that the previous studies extremely focused on human genome. Furthermore, they only focused on analyses that require only gene or protein sequences but not conducted a systematic investigation on other promoter features such as the binding evidence of specific transcription factors (TFs). This is mainly because of the incomplete resources of higher organisms, though they are relatively of interest, than model organisms such as Saccharomyces cerevisiae. The authors of this study recently integrated nine promoter features of 6603 genes of S. cerevisiae from six databases and five papers. These resources are suitable to conduct a comprehensive analysis of h2h genes in S. cerevisiae.

Results

This study analyzed various promoter features, including transcription boundaries (TSS, 5'UTR and 3'UTR), TATA box, TF binding evidence, TF regulation evidence, DNA bendability and nucleosome occupancy. The expression profiles and gene ontology (GO) annotations were used to measure if two genes are associated. Based on these promoter features, we found that i) the frequency of h2h genes was close to the expectation, namely they were not relatively frequent in genome; ii) the distance between the TSSs of most h2h genes fell into the range of 0-600 bps and was more centralized in 0-200 bps of the highly associated ones; iii) the number of TFs that regulate both h2h genes influenced the co-expression and co-function of the genes, while the number of TFs that bind both h2h genes influenced only the co-expression of the genes; iv) the association of two h2h genes was influenced by the existence of specific TFs such as STP2; v) the association of h2h genes whose bidirectional promoters have no TATA box was slightly higher than those who have TATA boxes; vi) the association of two h2h genes was not influenced by the DNA bendability and nucleosome occupancy.

Conclusions

This study analyzed h2h genes with various promoter features that have not been used in analyzing h2h genes. The results can be applied to other genomes to confirm if the observations of this study are limited to S. cerevisiae or universal in most organisms.