Genome-wide promoter extraction and analysis in human, mouse, and rat

Minimum entropy framework identifies a novel class of genomic functional elements and reveals regulatory mechanisms at human disease loci

We introduce CoRE-BED, a framework trained using 19 epigenomic features in 33 major cell and tissue types to predict cell-type-specific regulatory function. CoRE-BED identifies nine functional classes de-novo, capturing both known and new regulatory categories. Notably, we describe a previously undercharacterized class that we term Development Associated Elements (DAEs), which are highly enriched in cell types with elevated regenerative potential and distinguished by the dual presence of either H3K4me2 and H3K9ac (an epigenetic signature associated with kinetochore assembly) or H3K79me3 and H4K20me1 (a signature associated with transcriptional pause release). Unlike bivalent promoters, which represent a transitory state between active and silenced promoters, DAEs transition directly to or from a non-functional state during stem cell differentiation and are proximal to highly expressed genes. CoRE-BED's interpretability facilitates causal inference and functional prioritization. Across 70 complex traits, distal insulators account for the largest mean proportion of SNP heritability (~49%) captured by the GWAS. Collectively, our results demonstrate the value of exploring non-conventional ways of regulatory classification that enrich for trait heritability, to complement existing approaches for cis-regulatory prediction.

Epigenetic aberrations of gene expression in a rat model of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is mostly triggered by environmental and life-style factors and may involve epigenetic aberrations. However, a comprehensive documentation of the link between the dysregulated epigenome, transcriptome, and liver carcinogenesis is lacking. In the present study, Fischer-344 rats were fed a choline-deficient (CDAA, cancer group) or choline-sufficient (CSAA, healthy group) L-amino acid-defined diet. At the end of 52 weeks, transcriptomic alterations in livers of rats with HCC tumours and healthy livers were investigated by RNA sequencing. DNA methylation and gene expression were assessed by pyrosequencing and quantitative reverse-transcription PCR (qRT-PCR), respectively. We discovered 1,848 genes that were significantly differentially expressed in livers of rats with HCC tumours (CDAA) as compared with healthy livers (CSAA). Upregulated genes in the CDAA group were associated with cancer-related functions, whereas macronutrient metabolic processes were enriched by downregulated genes. Changes of highest magnitude were detected in numerous upregulated genes that govern key oncogenic signalling pathways, including Notch, Wnt, Hedgehog, and extracellular matrix degradation. We further detected perturbations in DNA methylating and demethylating enzymes, which was reflected in decreased global DNA methylation and increased global DNA hydroxymethylation. Four selected upregulated candidates, Mmp12, Jag1, Wnt4, and Smo, demonstrated promoter hypomethylation with the most profound decrease in Mmp12. MMP12 was also strongly overexpressed and hypomethylated in human HCC HepG2 cells as compared with primary hepatocytes, which coincided with binding of Ten-eleven translocation 1 (TET1). Our findings provide comprehensive evidence for gene expression changes and dysregulated epigenome in HCC pathogenesis, potentially revealing novel targets for HCC prevention/treatment.

MicroRNA 195-5p Targets Foxo3 Promoter Region to Regulate Its Expression in Granulosa Cells

Forkhead box O3 (Foxo3) is a member of the FOXO subfamily within the forkhead box (FOX) family, which has been shown to be essential for ovarian follicular development and maturation. Previous studies have shown the abundant expression of miR-195-5p in the nuclei of porcine granulosa cells (GCs), suggesting its potential role during ovarian follicle growth. In this study, a conditional immortalized porcine granulosa cell (CIPGC) line was used to determine whether the expression of Foxo3 could be regulated by the nuclear-enriched miR-195-5p. Through silico target prediction, we identified a potential binding site of miR-195-5p within the Foxo3 promoter. The over-expression of miR-195-5p increased Foxo3 expression at both mRNA and protein levels, while the knockdown of miR-195-5p decreased the expression of Foxo3. Furthermore, driven by the Foxo3 promoter, luciferase reporter activity was increased in response to miR-195-5p, while the mutation of the miR-195-5p binding site in the promoter region abolished this effect. In addition, the siRNA knockdown of Argonaute (AGO) 2, but not AGO1, significantly decreased Foxo3 transcript level. However, miR-195-5p failed to upregulate Foxo3 expression when AGO2 was knocked down. Moreover, chromatin immunoprecipitation (CHIP) assay showed that anti-AGO2 antibody pulled down both AGO2 and the Foxo3 promoter sequence, suggesting that AGO2 may be required for miR-195-5p to regulate Foxo3 expression in the nucleus. Additionally, Foxo3 expression was significantly increased by valproic acid (VPA), the inhibitor of deacetylase, as well as by methyltransferase inhibitor BIX-01294, indicating the involvement of histone modification. These effects were further enhanced in the presence of miR-195-5p and were decreased when miR-195-5p was knocked down. Overall, our results suggest that nuclear-enriched miR-195-5p regulates Foxo3 expression, which may be associated with AGO2 recruitment, as well as histone demethylation and acetylation in ovarian granulosa cells.

LEF1 supports metastatic brain colonization by regulating glutathione metabolism and increasing ROS resistance in breast cancer

More than half of all brain metastases show infiltrating rather than displacing growth at the macro-metastasis/organ parenchyma interface (MMPI), a finding associated with shorter survival. The lymphoid enhancer-binding factor-1 (LEF1) is an epithelial-mesenchymal transition (EMT) transcription factor that is commonly overexpressed in brain-colonizing cancer cells. Here, we overexpressed LEF1 in an in vivo breast cancer brain colonization model. It shortened survival, albeit without engaging EMT at the MMPI. By differential proteome analysis, we identified a novel function of LEF1 as a regulator of the glutathione (GSH) system, the principal cellular redox buffer. LEF1 overexpression also conferred resistance against therapeutic GSH depletion during brain colonization and improved management of intracellular ROS. We conclude that besides EMT, LEF1 facilitates metastasis by improving the antioxidative capacity of epithelial breast cancer cells, in particular during colonization of the brain parenchyma.

Benzyl butyl phthalate (BBP) triggers the malignancy of acute myeloid leukemia cells via upregulation of PDK4

Acute Myeloid Leukemia (AML) is a cancer of hematopoietic stem cells with a rapid progression. Recent studies indicated that endocrine disruptor chemicals (EDCs) are potential risk factors for AML progression. Our present data showed that an industrial endocrine disrupting chemical, Benzyl butyl phthalate (BBP), can promote the proliferation of AML cells and decrease their sensitivity to daunorubicin (DNR) and cytarabine (Ara-C) treatments. Further, BBP can increase the glucose consumption, lactate generation, and ATP levels of AML cells. Among the measured glycolysis-related genes, BBP can increase the expression of pyruvate dehydrogenase lipoamide kinase isozyme 4 (PDK4), a mitochondrial protein that regulates the tricarboxylic acid cycle (TCA) cycle. The inhibitor of PDK4 or its specific siRNA can attenuate BBP-induced cell proliferation and ATP generation, which suggested the essential roles of PDK4 in BBP-induced glycolysis and proliferation. Further, BBP can increase the mRNA stability of PDK4, while had no effect on its transcription and protein stability. miR-15b-5p can bind with the 3'UTR of PDK4 to decrease its mRNA stability, while BBP can decrease the expression of miR-15b-5p in AML cells. Collectively, our data showed that BBP can trigger the malignancy of AML cells via regulation of miR-15b-5p/PDK4 signals.

Epigenetic silencing of miR-520c leads to induced S100A4 expression and its mediated colorectal cancer progression

The S100 calcium-binding protein A4 (S100A4) induces epithelial mesenchymal transition, migration, invasion, angiogenesis and metastasis. Its induced expression in several cancer types correlates with poor prognosis. Apart from the functional and transcriptional regulatory aspects of S100A4, its post-transcriptional regulation is not yet clearly elucidated. In this study, we show that microRNAs (miR) miR-505-5p and miR-520c-3p target the 3′-UTR of S100A4 and inhibits its expression and its mediated migration and invasion. 5-Aza treatment significantly increased miR-520c-3p expression and reduced the S100A4 protein amounts. The upstream promoter region of miR-520c is hypermethylated irrespective of the metastasis status of colorectal cancer (CRC) patient tissues and in all analyzed CRC cell lines. Moreover, in a cohort of CRC patient specimen (n = 59), miR-520c-3p was significantly downregulated. miR-520c-3p stably expressing HCT116 cells showed a reduced metastasis formation in livers after implanting in mice spleen. Taken together, our findings demonstrate that S100A4 is post-transcriptionally regulated by tumor suppressor miRs, miR-505c-5p and miR-520c-3p, and particularly miR-520c-3p expression is epigenetically silenced in CRC.

1,25(OH)2 vitamin D(3) contributes to osteoclast-like trans-differentiation of malignant plasma cells

1,25-dihydroxyvitamin D (1,25(OH)₂D) exerts pleiotropic effects including bone turnover and immune system regulation. It inhibits both T and B cell proliferation while decreasing the secretion of inflammatory cytokines and immunoglobulins. 1,25(OH)₂D also modulates monocyte-macrophage and osteoclast (OC) maturation. Since we have previously described that malignant plasma cells may trans-differentiate towards the myeloid lineage participating to skeletal devastation in multiple myeloma (MM), we here evaluated in vitro the role of 1,25(OH)₂D in this lineage switch. We investigated the gene and protein expression of vitamin D receptor (VDR) in MM cell lines. Thus, after cell treatment with 1,25(OH)₂D, we analyzed their morphology and the expression of myeloid and OC markers. Finally, we assessed their bone resorption property on calcium phosphate slices. All MM cells expressed VDR in nuclear and perinuclear sites. Treatment with 1,25(OH)₂D altered their morphology from round to fusiform, while inducing paxillin focalization. 1,25(OH)₂D administration also up-regulated myeloid and OC genes, including C/EBPα, RANK, M-CSFR and V-ATPase, whose promoters contain potential 1,25(OH)₂D responsive elements. Finally, 1,25(OH)₂D increased MM cell capability to generate pits of erosion on calcium phosphate discs. This data suggest that myeloma cells may undergo a functional trans-differentiation into OCs and, under appropriate experimental conditions, 1,25(OH)₂D triggers this lineage switch.

Adaptation of the targeted capture Methyl-Seq platform for the mouse genome identifies novel tissue-specific DNA methylation patterns of genes involved in neurodevelopment

Methyl-Seq was recently developed as a targeted approach to assess DNA methylation (DNAm) at a genome-wide level in human. We adapted it for mouse and sought to examine DNAm differences across liver and 2 brain regions: cortex and hippocampus. A custom hybridization array was designed to isolate 99 Mb of CpG islands, shores, shelves, and regulatory elements in the mouse genome. This was followed by bisulfite conversion and sequencing on the Illumina HiSeq2000. The majority of differentially methylated cytosines (DMCs) were present at greater than expected frequency in introns, intergenic regions, near CpG islands, and transcriptional enhancers. Liver-specific enhancers were observed to be methylated in cortex, while cortex specific enhancers were methylated in the liver. Interestingly, commonly shared enhancers were differentially methylated between the liver and cortex. Gene ontology and pathway analysis showed that genes that were hypomethylated in the cortex and hippocampus were enriched for neuronal components and neuronal function. In contrast, genes that were hypomethylated in the liver were enriched for cellular components important for liver function. Bisulfite-pyrosequencing validation of 75 DMCs from 19 different loci showed a correlation of r = 0.87 with Methyl-Seq data. We also identified genes involved in neurodevelopment that were not previously reported to be differentially methylated across brain regions. This platform constitutes a valuable tool for future genome-wide studies involving mouse models of disease.

Directing an artificial zinc finger protein to new targets by fusion to a non-DNA-binding domain

Transcription factors are often regarded as having two separable components: a DNA-binding domain (DBD) and a functional domain (FD), with the DBD thought to determine target gene recognition. While this holds true for DNA binding in vitro, it appears that in vivo FDs can also influence genomic targeting. We fused the FD from the well-characterized transcription factor Krüppel-like Factor 3 (KLF3) to an artificial zinc finger (AZF) protein originally designed to target the Vascular Endothelial Growth Factor-A (VEGF-A) gene promoter. We compared genome-wide occupancy of the KLF3FD-AZF fusion to that observed with AZF. AZF bound to the VEGF-A promoter as predicted, but was also found to occupy approximately 25 000 other sites, a large number of which contained the expected AZF recognition sequence, GCTGGGGGC. Interestingly, addition of the KLF3 FD re-distributes the fusion protein to new sites, with total DNA occupancy detected at around 50 000 sites. A portion of these sites correspond to known KLF3-bound regions, while others contained sequences similar but not identical to the expected AZF recognition sequence. These results show that FDs can influence and may be useful in directing AZF DNA-binding proteins to specific targets and provide insights into how natural transcription factors operate.

Skp2 deletion unmasks a p27 safeguard that blocks tumorigenesis in the absence of pRb and p53 tumor suppressors

pRb and p53 are two major tumor suppressors. Here, we found that p53 activates expression of Pirh2 and KPC1, two of the three ubiquitin ligases for p27. Loss of p53 in the absence of Skp2, the third ubiquitin ligase for p27, shrinks the cellular pool of p27 ubiquitin ligases to accumulate p27 protein. In the absence of pRb and p53, p27 was unable to inhibit DNA synthesis in spite of its abundance, but could inhibit division of cells that maintain DNA replication with rereplication. This mechanism blocked pRb/p53 doubly deficient pituitary and prostate tumorigenesis lastingly coexistent with bromodeoxyuridine-labeling neoplastic lesions, revealing an unconventional cancer cell vulnerability when pRb and p53 are inactivated.

Deciphering the calcitriol-induced transcriptomic response in keratinocytes: presentation of novel target genes

Numerous studies to date have been aimed at unraveling the large suite of calcitriol (1α,25-dihydroxyvitamin D(3)) response genes in diverse tissues including skin, where this hormone is involved in regulating keratinocyte proliferation, differentiation, permeability barrier formation, innate immunity promotion, antimicrobial peptide production, and wound healing. However, the various approaches differ considerably in probed cell types, scale, throughput, and statistical reliability and do, of note, not reveal much overlap. To further expand our knowledge on presently elusive targets and characterize the extent of fragmentation of existing datasets, we have performed whole-transcriptome microarray examinations of calcitriol-treated human primary keratinocytes. Out of 28,869 genes investigated, we uncovered 86 differentially expressed (67 upregulated and 19 downregulated) candidates that were functionally clustered into five annotation categories: response to wounding, protease inhibition, secondary metabolite biosynthesis, cellular migration, and amine biosynthetic processes. A complementary RTq-PCR study of 78 nominees selected thereof demonstrated significant differential expression of 55 genes (48 upregulated and seven downregulated) within biological replicates. Our hit list contains nine previously authenticated targets (16.36%, proof of concept) and 46 novel genes (83.6%) that have not yet been explicitly described as being differentially regulated within human primary keratinocytes. Direct vitamin D receptor response element predictions within the regulatory promoter regions of 50 of the RTq-PCR-validated targets agreed with known biological functionality and corroborated our stringent data validation pipeline. Altogether, our results indicate the value of continuing these kinds of gene expression studies, which contribute to an enhanced comprehension of calcitriol-mediated processes that may be dysregulated in human skin pathophysiology.

Application of an integrative computational framework in trancriptomic data of atherosclerotic mice suggests numerous molecular players

Atherosclerosis is a multifactorial disease involving a lot of genes and proteins recruited throughout its manifestation. The present study aims to exploit bioinformatic tools in order to analyze microarray data of atherosclerotic aortic lesions of ApoE knockout mice, a model widely used in atherosclerosis research. In particular, a dynamic analysis was performed among young and aged animals, resulting in a list of 852 significantly altered genes. Pathway analysis indicated alterations in critical cellular processes related to cell communication and signal transduction, immune response, lipid transport, and metabolism. Cluster analysis partitioned the significantly differentiated genes in three major clusters of similar expression profile. Promoter analysis applied to functional related groups of the same cluster revealed shared putative cis-elements potentially contributing to a common regulatory mechanism. Finally, by reverse engineering the functional relevance of differentially expressed genes with specific cellular pathways, putative genes acting as hubs, were identified, linking functionally disparate cellular processes in the context of traditional molecular description.

GPMiner: an integrated system for mining combinatorial cis-regulatory elements in mammalian gene group

Background

Sequence features in promoter regions are involved in regulating gene transcription initiation. Although numerous computational methods have been developed for predicting transcriptional start sites (TSSs) or transcription factor (TF) binding sites (TFBSs), they lack annotations for do not consider some important regulatory features such as CpG islands, tandem repeats, the TATA box, CCAAT box, GC box, over-represented oligonucleotides, DNA stability, and GC content. Additionally, the combinatorial interaction of TFs regulates the gene group that is associated with same expression pattern. To investigate gene transcriptional regulation, an integrated system that annotates regulatory features in a promoter sequence and detects co-regulation of TFs in a group of genes is needed.

Results

This work identifies TSSs and regulatory features in a promoter sequence, and recognizes co-occurrence of cis-regulatory elements in co-expressed genes using a novel system. Three well-known TSS prediction tools are incorporated with orthologous conserved features, such as CpG islands, nucleotide composition, over-represented hexamer nucleotides, and DNA stability, to construct the novel Gene Promoter Miner (GPMiner) using a support vector machine (SVM). According to five-fold cross-validation results, the predictive sensitivity and specificity are both roughly 80%. The proposed system allows users to input a group of gene names/symbols, enabling the co-occurrence of TFBSs to be determined. Additionally, an input sequence can also be analyzed for homogeneity of experimental mammalian promoter sequences, and conserved regulatory features between homologous promoters can be observed through cross-species analysis. After identifying promoter regions, regulatory features are visualized graphically to facilitate gene promoter observations.

Conclusions

The GPMiner, which has a user-friendly input/output interface, has numerous benefits in analyzing human and mouse promoters. The proposed system is freely available at http://GPMiner.mbc.nctu.edu.tw/.

Histone modification profiles are predictive for tissue/cell-type specific expression of both protein-coding and microRNA genes

Background

Gene expression is regulated at both the DNA sequence level and through modification of chromatin. However, the effect of chromatin on tissue/cell-type specific gene regulation (TCSR) is largely unknown. In this paper, we present a method to elucidate the relationship between histone modification/variation (HMV) and TCSR.

Results

A classifier for differentiating CD4+ T cell-specific genes from housekeeping genes using HMV data was built. We found HMV in both promoter and gene body regions to be predictive of genes which are targets of TCSR. For example, the histone modification types H3K4me3 and H3K27ac were identified as the most predictive for CpG-related promoters, whereas H3K4me3 and H3K79me3 were the most predictive for nonCpG-related promoters. However, genes targeted by TCSR can be predicted using other type of HMVs as well. Such redundancy implies that multiple type of underlying regulatory elements, such as enhancers or intragenic alternative promoters, which can regulate gene expression in a tissue/cell-type specific fashion, may be marked by the HMVs. Finally, we show that the predictive power of HMV for TCSR is not limited to protein-coding genes in CD4+ T cells, as we successfully predicted TCSR targeted genes in muscle cells, as well as microRNA genes with expression specific to CD4+ T cells, by the same classifier which was trained on HMV data of protein-coding genes in CD4+ T cells.

Conclusion

We have begun to understand the HMV patterns that guide gene expression in both tissue/cell-type specific and ubiquitous manner.

Atlas of gene expression in the mouse kidney: new features of glomerular parietal cells

To gain molecular insight into kidney function, we performed a high-resolution quantitative analysis of gene expression in glomeruli and nine different nephron segments dissected from mouse kidney using Serial Analysis of Gene Expression (SAGE). We also developed dedicated bioinformatics tools and databases to annotate mRNA tags as transcripts. Over 800,000 mRNA SAGE tags were sequenced corresponding to >20,000 different mRNA tags present at least twice in at least one library. Hierarchical clustering analysis of tags demonstrated similarities between the three anatomical subsegments of the proximal tubule, between the cortical and medullary segments of the thick ascending limb of Henle's loop, and between the three segments constituting the aldosterone-sensitive distal nephron segments, whereas the glomerulus and distal convoluted tubule clusterized independently. We also identified highly specific mRNA markers of each subgroup of nephron segments and of most nephron segments. Tag annotation also identified numbers of putative antisense mRNAs. This database constitutes a reference resource in which the quantitative expression of a given gene can be compared with that of other genes in the same nephron segment, or between different segments of the nephron. To illustrate possible applications of this database, we performed a deeper analysis of the glomerulus transcriptome that unexpectedly revealed expression of several ion and water carriers; within the glomerulus, they were found to be preferentially expressed in the parietal sheet. It also revealed the major role of the zinc finger transcription factor Wt1 in the specificity of gene expression in the glomerulus. Finally, functional annotation of glomerulus-specific transcripts suggested a high proliferation activity of glomerular cells. Immunolabeling for PCNA confirmed a high percentage of proliferating cells in the glomerulus parietal sheet.

Identification of cis-regulatory elements in the mammalian genome: the cREMaG database

Background

A growing number of gene expression-profiling datasets provides a reliable source of information about gene co-expression. In silico analyses of the properties shared among the promoters of co-expressed genes facilitates the identification of transcription factors (TFs) involved in the co-regulation of those genes. Our previous experience with microarray data led to the development of a database suitable for the examination of regulatory motifs in the promoters of co-expressed genes.

Methodology

We introduce the cREMaG (cis-Regulatory Elements in the Mammalian Genome) system designed for in silico studies of the promoter properties of co-regulated mammalian genes. The cREMaG system offers an analysis of data obtained from human, mouse, rat, bovine and canine gene expression-profiling studies. More than eight analysis parameters can be utilized in user-defined combinations. The selection of alternative transcription start sites and information about CpG islands are also available.

Conclusions

Using the cREMaG system, we successfully identified TFs mediating transcriptional responses in reference gene sets. The cREMaG system facilitates in silico studies of mammalian transcriptional gene regulation. The resource is freely available at http://www.cremag.org.

Cross-species comparison of genomewide gene expression profiles reveals induction of hypoxia-inducible factor-responsive genes in iron-deprived intestinal epithelial cells

Molecular mechanisms mediating the induction of metal ion homeostasis-related genes in the mammalian intestine during iron deficiency remain unknown. To elucidate relevant regulatory pathways, genomewide gene expression profiles were determined in fully differentiated human intestinal epithelial (Caco-2) cells. Cells were deprived of iron (or not) for 6 or 18 h, and Gene Chip analyses were subsequently performed (Affymetrix). More than 2,000 genes were differentially expressed; genes related to monosaccharide metabolism, regulation of gene expression, hypoxia, and cell death were upregulated, while those related to mitotic cell cycle were downregulated. A large proportion of induced genes are hypoxia responsive, and promoter enrichment analyses revealed a statistical overrepresentation of hypoxia response elements (HREs). Immunoblot experiments demonstrated a >60-fold increase in HIF2α protein abundance in iron-deprived cells; HIF1α levels were unchanged. Furthermore, comparison of the Caco-2 cell data set with a Gene Chip data set from iron-deficient rat intestine revealed 29 common upregulated genes; the majority are hypoxia responsive, and their promoters are enriched for HREs. We conclude that the compensatory response of the intestinal epithelium to iron deprivation relates to hypoxia and that stabilization of HIF2α may be the primary event mediating metabolic and morphological changes observed during iron deficiency.

A transcriptomic computational analysis of mastic oil-treated Lewis lung carcinomas reveals molecular mechanisms targeting tumor cell growth and survival

Background

Mastic oil from Pistacia lentiscus variation chia, a blend of bioactive terpenes with recognized medicinal properties, has been recently shown to exert anti-tumor growth activity through inhibition of cancer cell proliferation, survival, angiogenesis and inflammatory response. However, no studies have addressed its mechanisms of action at genome-wide gene expression level.

Methods

To investigate molecular mechanisms triggered by mastic oil, Lewis Lung Carcinoma cells were treated with mastic oil or DMSO and RNA was collected at five distinct time points (3-48 h). Microarray expression profiling was performed using Illumina mouse-6 v1 beadchips, followed by computational analysis. For a number of selected genes, RT-PCR validation was performed in LLC cells as well as in three human cancer cell lines of different origin (A549, HCT116, K562). PTEN specific inhibition by a bisperovanadium compound was applied to validate its contribution to mastic oil-mediated anti-tumor growth effects.

Results

In this work we demonstrated that exposure of Lewis lung carcinomas to mastic oil caused a time-dependent alteration in the expression of 925 genes. GO analysis associated expression profiles with several biological processes and functions. Among them, modifications on cell cycle/proliferation, survival and NF-κB cascade in conjunction with concomitant regulation of genes encoding for PTEN, E2F7, HMOX1 (up-regulation) and NOD1 (down-regulation) indicated some important mechanistic links underlying the anti-proliferative, pro-apoptotic and anti-inflammatory effects of mastic oil. The expression profiles of Hmox1, Pten and E2f7 genes were similarly altered by mastic oil in the majority of test cancer cell lines. Inhibition of PTEN partially reversed mastic oil effects on tumor cell growth, indicating a multi-target mechanism of action. Finally, k-means clustering, organized the significant gene list in eight clusters demonstrating a similar expression profile. Promoter analysis in a representative cluster revealed shared putative cis-elements suggesting a common regulatory transcription mechanism.

Conclusions

Present results provide novel evidence on the molecular basis of tumor growth inhibition mediated by mastic oil and set a rational basis for application of genomics and bioinformatic methodologies in the screening of natural compounds with potential cancer chemopreventive activities.

Computational design and application of endogenous promoters for transcriptionally targeted gene therapy for rheumatoid arthritis

The promoter regions of genes that are differentially regulated in the synovial membrane during the course of rheumatoid arthritis (RA) represent attractive candidates for application in transcriptionally targeted gene therapy. In this study, we applied an unbiased computational approach to define proximal-promoters from a gene expression profiling study of murine experimental arthritis. Synovium expression profiles from progressing stages of collagen-induced arthritis (CIA) were classified into six distinct groups using k-means clustering. Using an algorithm based on local over-representation and comparative genomics, we identified putatively functional transcription factor-binding sites (TFBS) in TATA-dependent proximal-promoters. Applying a filter based on spacing between TATA box and transcription start site (TSS) combined with the presence of over-represented nuclear factor kappaB (NFkappaB), AP-1, or CCAAT/enhancer-binding protein beta (C/EBPbeta) sites, 382 candidate murine and human promoters were reduced to 66, corresponding to 45 genes. In vitro, 9 out of 10 computationally defined promoter regions conferred cytokine-inducible expression in murine cells and human synovial fibroblasts. Under these conditions, the serum amyloid A3 (Saa3) promoter showed the strongest transcriptional induction and strength. We applied this promoter for driving therapeutically efficacious levels of the interleukin-1 receptor antagonist (Il1rn) in a disease-regulated fashion. These results demonstrate the value of bioinformatics for guiding the selection of endogenous promoters for transcriptionally targeted gene therapy.

FactorY, a bioinformatic resource for genome-wide promoter analysis

The interpretation of the complex molecular descriptions generated by high-throughput gene expression technologies is still challenging. The development of new tools to identify common regulatory mechanisms involved in the control of the expression of a set of co-expressed genes, might enhance our capacity to extract functional information from genomic data sets. Here we present FactorY, a website that allows identification of enriched transcription factor binding sites (TFBSs) in the proximal promoter of a cluster of genes, as well as functional interpretation, and intuitive visualization of the results.

Systematic analysis of alternative promoters correlated with alternative splicing in human genes

Interactions between various events are essential for complex and delicate transcriptional regulation. To delineate the features and potential roles of alternative promoters (APs) correlated with alternative splicing (AS), we have systematically analyzed 9908 putative alternative promoters (PAPs) from 3797 human genes. Our results showed that approximately 65% of AS events are associated with PAPs. Intriguingly, PAPs per human AS gene only averaged 2.6 for our dataset, which was significantly lower than previously reported. This seems to imply that the human genome contains a small pool of appropriable PAPs for AS genes. Exploration of the characteristics of PAPs such as CpG islands, TATA boxes, GC-content, transcription factor binding sites (TFBSs) and repetitive elements suggested that, respectively, 87% and 90% of PAPs of human AS genes are CpG- and TATA box-poor. The GC-content is significantly higher in the downstream of transcription start sites (TSSs) than upstream (58% vs. 53%), and there is a strong negative correlation between the GC-content and the number of PAPs. These suggested that GC-content around the TSSs plays an important role in the regulation of AS. Moreover, different APs contain distinct densities of repetitive elements and TFBSs, indicating that such sequences have an intrinsic role in the divergent regulation of PAPs and AS. Substantial difference was also found between human AS genes in terms of PAP numbers. A close connection between PAPs and AS may play a critical role in the choice of APs and regulation of AS genes. Furthermore, the distribution of AS genes on different human chromosomes also influences the numbers of PAPs and isoforms of AS genes. Our results may provide important clues for further studies on regulatory network of transcription-related events.

High-resolution human core-promoter prediction with CoreBoost_HM

Correctly locating the gene transcription start site and the core-promoter is important for understanding transcriptional regulation mechanism. Here we have integrated specific genome-wide histone modification and DNA sequence features together to predict RNA polymerase II core-promoters in the human genome. Our new predictor CoreBoost_HM outperforms existing promoter prediction algorithms by providing significantly higher sensitivity and specificity at high resolution. We demonstrated that even though the histone modification data used in this study are from a specific cell type (CD4+ T-cell), our method can be used to identify both active and repressed promoters. We have applied it to search the upstream regions of microRNA genes, and show that CoreBoost_HM can accurately identify the known promoters of the intergenic microRNAs. We also identified a few intronic microRNAs that may have their own promoters. This result suggests that our new method can help to identify and characterize the core-promoters of both coding and noncoding genes.

Characterization of mouse synaptic vesicle-2-associated protein (Msvop) specifically expressed in the mouse central nervous system

In this study, we identified and characterized a mouse brain-enriched unigene (msvop) after performing digital differential display. Msvop was the mouse ortholog of Xenopus synaptic vesicle-2-associated protein (svop), the molecular characteristics of which were unknown. The 3125-bp full-length cDNA encoded a 548-aa protein of approximately 60 kDa. A strong promoter element was found in the -200 to -100 bp region in both NG108-15 and HEK293 cells. RT-PCR and in situ hybridization analysis confirmed that msvop was strictly expressed in the mouse central nervous system. In adult brains, msvop was highly expressed in the hippocampus and cerebellum. When the gene was transfected into HEK293 cells in a GFP fusion vector, the protein was specifically localized in the cytosol. These results indicate that msvop is a central nervous system-specific gene and could be utilized for elucidating gene regulation in neuronal cells.

Regulation of the PDK4 isozyme by the Rb-E2F1 complex

Loss of the transcription factor E2F1 elicits a complex metabolic phenotype in mice underscored by reduced adiposity and protection from high fat diet-induced diabetes. Here, we demonstrate that E2F1 directly regulates the gene encoding PDK4 (pyruvate dehydrogenase kinase 4), a key nutrient sensor and modulator of glucose homeostasis that is chronically elevated in obesity and diabetes and acutely induced under the metabolic stress of starvation or fasting. We show that loss of E2F1 in vivo blunts PDK4 expression and improves myocardial glucose oxidation. The absence of E2F1 also corresponds to lower blood glucose levels, improved plasma lipid profile, and increased sensitivity to insulin stimulation. Consistently, enforced E2F1 expression up-regulates PDK4 levels and suppresses glucose oxidation in C(2)C(12) myoblasts. Furthermore, inactivation of Rb, the repressor of E2F-dependent transcription, markedly induces PDK4 and triggers the enrichment of E2F1 occupancy onto the PDK4 promoter as detected by chromatin immunoprecipitation analysis. Two overlapping E2F binding sites were identified on this promoter. Transactivation assays later verified E2F1 responsiveness of this promoter element in C(2)C(12) myoblasts and IMR90 fibroblasts, an effect that was completely abrogated following mutation of the E2F sites. Taken together, our data illustrate how the E2F1 mitogen directly regulates PDK4 levels and influences cellular bioenergetics, namely mitochondrial glucose oxidation. These results are relevant to the pathophysiology of chronic diseases like obesity and diabetes, where PDK4 is dysregulated and could have implications pertinent to the etiology of tumor metabolism, especially in cancers with Rb pathway defects.

Distance conservation of transcription regulatory motifs in human promoters

To understanding the interaction network among transcription-regulation elements in human is an immediate challenge for modern molecular biology. Here a central problem is how to extract evolutionary information and search the evolutionary conservation from the comparison of promoters of closely related species. Through the comparative studies of k-mer distribution in human and mouse transcription factor binding site (TFBS) sequences we have discovered that the average distance between a pair of transcription regulatory 7-mer motifs is conservative in human-mouse promoters. The distance conservation is a new kind of evolutionary conservation, not based on the strict location of bases in genome sequence. By utilizing the conservation of k-mer distance it will be helpful to propose a non-alignment-based approach for fast genome-wide discovery of transcription regulatory motifs. We demonstrated the distance conservation by genome-wide searching of conservative regulatory 7-mer motifs with successful rate 90%. Then, after defining human-mouse pair-distance divergence parameter we studied the tissue-specific motif pairs and found that the parameter for motif pairs is 11-16 times smaller than for their controls for 28 tissues and these pairs can be clearly differentiated on two-dimensional parameter plane. Finally, the mechanism of distance conservation was discussed briefly which is supposed to be related to the module structure of TFBSs.

In silico regulatory analysis for exploring human disease progression

Background

An important goal in bioinformatics is to unravel the network of transcription factors (TFs) and their targets. This is important in the human genome, where many TFs are involved in disease progression. Here, classification methods are applied to identify new targets for 152 transcriptional regulators using publicly-available targets as training examples. Three types of sequence information are used: composition, conservation, and overrepresentation.

Results

Starting with 8817 TF-target interactions we predict an additional 9333 targets for 152 TFs. Randomized classifiers make few predictions (~2/18660) indicating that our predictions for many TFs are significantly enriched for true targets. An enrichment score is calculated and used to filter new predictions.

Two case-studies for the TFs OCT4 and WT1 illustrate the usefulness of our predictions:

• Many predicted OCT4 targets fall into the Wnt-pathway. This is consistent with known biology as OCT4 is developmentally related and Wnt pathway plays a role in early development.

• Beginning with 15 known targets, 354 predictions are made for WT1. WT1 has a role in formation of Wilms' tumor. Chromosomal regions previously implicated in Wilms' tumor by cytological evidence are statistically enriched in predicted WT1 targets. These findings may shed light on Wilms' tumor progression, suggesting that the tumor progresses either by loss of WT1 or by loss of regions harbouring its targets.

• Targets of WT1 are statistically enriched for cancer related functions including metastasis and apoptosis. Among new targets are BAX and PDE4B, which may help mediate the established anti-apoptotic effects of WT1.

• Of the thirteen TFs found which co-regulate genes with WT1 (p ≤ 0.02), 8 have been previously implicated in cancer. The regulatory-network for WT1 targets in genomic regions relevant to Wilms' tumor is provided.

Conclusion

We have assembled a set of features for the targets of human TFs and used them to develop classifiers for the determination of new regulatory targets. Many predicted targets are consistent with the known biology of their regulators, and new targets for the Wilms' tumor regulator, WT1, are proposed. We speculate that Wilms' tumor development is mediated by chromosomal rearrangements in the location of WT1 targets.

Reviewers

This article was reviewed by Trey Ideker, Vladimir A. Kuznetsov(nominated by Frank Eisenhaber), and Tzachi Pilpel.

Lindernia brevidens: a novel desiccation-tolerant vascular plant, endemic to ancient tropical rainforests

A particular adaptation to survival under limited water availability has been realized in the desiccation-tolerant resurrection plants, which tend to grow in a habitat with seasonal rainfall and long dry periods. One of the best-studied examples is Craterostigma plantagineum. Here we report an unexpected finding: Lindernia brevidens, a close relative of C. plantagineum, exhibits desiccation tolerance, even though it is endemic to the montane rainforests of Tanzania and Kenya, where it never experiences seasonal dry periods. L. brevidens has been found exclusively in two fragments of the ancient Eastern Arc Mountains, which were protected from the devastating Pleistocene droughts by the stable Indian Ocean temperature. Analysis of the microhabitat reveals that L. brevidens is found in the same habitat as hygrophilous plant species, which further indicates that the plant never dries out completely. The objective of this investigation was to address whether C. plantagineum and L. brevidens have desiccation-related pathways in common, or whether L. brevidens has acquired novel pathways. A third, closely related, desiccation-sensitive species, Lindernia subracemosa, has been included for comparison. Mechanisms that confer cellular protection during extreme water loss are well conserved between C. plantagineum and L. brevidens, including the interconversion of 2-octulose to sucrose within the two desiccation-tolerant species. Furthermore, transcriptional control regions of desiccation-related genes belonging to the late embryogenesis abundant (LEA) protein family are also highly conserved. We propose that L. brevidens is a neoendemic species that has retained desiccation tolerance through genome stability, despite tolerance being superfluous to environmental conditions.

FOXOs, cancer and regulation of apoptosis

Forkhead box O (FOXO) transcription factors are involved in multiple signaling pathways and play critical roles in a number of physiological and pathological processes including cancer. The importance of FOXO factors ascribes them under multiple levels of regulation including phosphorylation, acetylation/deacetylation, ubiquitination and protein-protein interactions. As FOXO factors play a pivotal role in cell fate decision, mounting evidence suggests that FOXO factors function as tumor suppressors in a variety of cancers. FOXOs are actively involved in promoting apoptosis in a mitochondria-independent and -dependent manner by inducing the expression of death receptor ligands, including Fas ligand and tumor necrosis factor-related apoptosis-inducing ligand, and Bcl-2 family members, such as Bim, bNIP3 and Bcl-X(L), respectively. An understanding of FOXO proteins and their biology will provide new opportunities for developing more effective therapeutic approaches to treat cancer.

CREB binding and activity in brain: regional specificity and induction by electroconvulsive seizure

BACKGROUND

The transcription factor cyclic adenosine monophosphate response element binding protein (CREB) orchestrates diverse neurobiological processes including cell differentiation, survival, and plasticity. Alterations in CREB-mediated transcription have been implicated in numerous central nervous system (CNS) disorders including depression, anxiety, addiction, and cognitive decline. However, it remains unclear how CREB contributes to normal and aberrant CNS function, as the identity of CREB-regulated genes in brain and the regional and temporal dynamics of CREB function remain largely undetermined.

METHODS

We combined microarray and chromatin immunoprecipitation technology to analyze CREB-DNA interactions in brain. We compared the occupancy and activity of CREB at gene promoters in rat frontal cortex, hippocampus, and striatum before and after a rodent model of electroconvulsive therapy.

RESULTS

Our analysis identified >860 CREB binding sites in rat brain. We identified multiple genomic loci enriched with CREB binding sites and find that CREB-occupied transcripts interact extensively to promote cell proliferation, plasticity, and resiliency. We discovered regional differences in CREB occupancy and activity that explain, in part, the diverse biological and behavioral outputs of CREB activity in frontal cortex, hippocampus, and striatum. Electroconvulsive seizure rapidly increased CREB occupancy and/or phosphorylation at select promoters, demonstrating that both events contribute to the temporal regulation of the CREB transcriptome.

CONCLUSIONS

Our data provide a mechanistic basis for CREB's ability to integrate regional and temporal cues to orchestrate state-specific patterns of transcription in the brain, indicate that CREB is an important mediator of the biological responses to electroconvulsive seizure, and provide global mechanistic insights into CREB's role in psychiatric and cognitive function.

CD36-dependent regulation of muscle FoxO1 and PDK4 in the PPAR delta/beta-mediated adaptation to metabolic stress

The transcription factor FoxO1 contributes to the metabolic adaptation to fasting by suppressing muscle oxidation of glucose, sparing it for glucose-dependent tissues. Previously, we reported that FoxO1 activation in C(2)C(12) muscle cells recruits the fatty acid translocase CD36 to the plasma membrane and increases fatty acid uptake and oxidation. This, together with FoxO1 induction of lipoprotein lipase, would promote the reliance on fatty acid utilization characteristic of the fasted muscle. Here, we show that CD36-mediated fatty acid uptake, in turn, up-regulates protein levels and activity of FoxO1 as well as its target PDK4, the negative regulator of glucose oxidation. Increased fatty acid flux or enforced CD36 expression in C(2)C(12) cells is sufficient to induce FoxO1 and PDK4, whereas CD36 knockdown has opposite effects. In vivo, CD36 loss blunts fasting induction of FoxO1 and PDK4 and the associated suppression of glucose oxidation. Importantly, CD36-dependent regulation of FoxO1 is mediated by the nuclear receptor PPARdelta/beta. Loss of PPARdelta/beta phenocopies CD36 deficiency in blunting fasting induction of muscle FoxO1 and PDK4 in vivo. Expression of PPARdelta/beta in C(2)C(12) cells, like that of CD36, robustly induces FoxO1 and suppresses glucose oxidation, whereas co-expression of a dominant negative PPARdelta/beta compromises FoxO1 induction. Finally, several PPRE sites were identified in the FoxO1 promoter, which was responsive to PPARdelta/beta. Agonists of PPARdelta/beta were sufficient to confer responsiveness and transactivate the heterologous FoxO1 promoter but not in the presence of dominant negative PPARdelta/beta. Taken together, our findings suggest that CD36-dependent FA activation of PPARdelta/beta results in the transcriptional regulation of FoxO1 as well as PDK4, recently shown to be a direct PPARdelta/beta target. FoxO1 in turn can regulate CD36, lipoprotein lipase, and PDK4, reinforcing the action of PPARdelta/beta to increase muscle reliance on FA. The findings could have implications in the chronic abnormalities of fatty acid metabolism associated with obesity and diabetes.

Life-history traits drive the evolutionary rates of mammalian coding and noncoding genomic elements

A comprehensive phylogenetic framework is indispensable for investigating the evolution of genomic features in mammals as a whole, and particularly in humans. Using the ENCODE sequence data, we estimated mammalian neutral evolutionary rates and selective pressures acting on conserved coding and noncoding elements. We show that neutral evolutionary rates can be explained by the generation time (GT) hypothesis. Accordingly, primates (especially humans), having longer GTs than other mammals, display slower rates of neutral evolution. The evolution of constrained elements, particularly of nonsynonymous sites, is in agreement with the expectations of the nearly neutral theory of molecular evolution. We show that rates of nonsynonymous substitutions (dN) depend on the population size of a species. The results are robust to the exclusion of hypermutable CpG prone sites. The average rate of evolution in conserved noncoding sequences (CNCs) is 1.7 times higher than in nonsynonymous sites. Despite this, CNCs evolve at similar or even lower rates than nonsynonymous sites in the majority of basal branches of the eutherian tree. This observation could be the result of an overall gradual or, alternatively, lineage-specific relaxation of CNCs. The latter hypothesis was supported by the finding that 3 of the 20 longest CNCs displayed significant relaxation of individual branches. This observation may explain why the evolution of CNCs fits the expectations of the nearly neutral theory less well than the evolution of nonsynonymous sites.

Genome-wide mapping and analysis of active promoters in mouse embryonic stem cells and adult organs

By integrating genome-wide maps of RNA polymerase II (Polr2a) binding with gene expression data and H3ac and H3K4me3 profiles, we characterized promoters with enriched activity in mouse embryonic stem cells (mES) as well as adult brain, heart, kidney, and liver. We identified approximately 24,000 promoters across these samples, including 16,976 annotated mRNA 5' ends and 5153 additional sites validating cap-analysis of gene expression (CAGE) 5' end data. We showed that promoters with CpG islands are typically non-tissue specific, with the majority associated with Polr2a and the active chromatin modifications in nearly all the tissues examined. By contrast, the promoters without CpG islands are generally associated with Polr2a and the active chromatin marks in a tissue-dependent way. We defined 4396 tissue-specific promoters by adapting a quantitative index of tissue-specificity based on Polr2a occupancy. While there is a general correspondence between Polr2a occupancy and active chromatin modifications at the tissue-specific promoters, a subset of them appear to be persistently marked by active chromatin modifications in the absence of detectable Polr2a binding, highlighting the complexity of the functional relationship between chromatin modification and gene expression. Our results provide a resource for exploring promoter Polr2a binding and epigenetic states across pluripotent and differentiated cell types in mammals.

Promoter analysis of intestinal genes induced during iron-deprivation reveals enrichment of conserved SP1-like binding sites

Background

Iron-deficiency leads to the induction of genes related to intestinal iron absorption and homeostasis. By analyzing a large GeneChip^® dataset from the rat intestine, we identified a large cluster of 228 genes that was induced by iron-deprivation. Only 2 of these genes contained 3' iron-response elements, suggesting that other regulation including transcriptional may be involved. We therefore utilized computational methods to test the hypothesis that some of the genes within this large up-regulated cluster are co-ordinately regulated by common transcriptional mechanisms. We thus identified promoters from the up-regulated gene cluster from rat, mouse and human, and performed enrichment analyses with the Clover program and the TRANSFAC database.

Results

Surprisingly, we found a strong statistical enrichment for SP1 binding sites in our experimental promoters as compared to background sequences. As the TRANSFAC database cannot distinguish among SP/KLF family members, many of which bind similar GC-rich DNA sequences, we surmise that SP1 or an SP1-like factor could be involved in this response. In fact, we detected induction of SP6/KLF14 in the GeneChip^® studies, and confirmed it by real-time PCR. Additional computational analyses suggested that an SP1-like factor may function synergistically with a FOX TF to regulate a subset of these genes. Furthermore, analysis of promoter sequences identified many genes with multiple, conserved SP1 and FOX binding sites, the relative location of which within orthologous promoters was highly conserved.

Conclusion

SP1 or a closely related factor may play a primary role in the genetic response to iron-deficiency in the mammalian intestine.

MetaProm: a neural network based meta-predictor for alternative human promoter prediction

BACKGROUND

De novo eukaryotic promoter prediction is important for discovering novel genes and understanding gene regulation. In spite of the great advances made in the past decade, recent studies revealed that the overall performances of the current promoter prediction programs (PPPs) are still poor, and predictions made by individual PPPs do not overlap each other. Furthermore, most PPPs are trained and tested on the most-upstream promoters; their performances on alternative promoters have not been assessed.

RESULTS

In this paper, we evaluate the performances of current major promoter prediction programs (i.e., PSPA, FirstEF, McPromoter, DragonGSF, DragonPF, and FProm) using 42,536 distinct human gene promoters on a genome-wide scale, and with emphasis on alternative promoters. We describe an artificial neural network (ANN) based meta-predictor program that integrates predictions from the current PPPs and the predicted promoters' relation to CpG islands. Our specific analysis of recently discovered alternative promoters reveals that although only 41% of the 3' most promoters overlap a CpG island, 74% of 5' most promoters overlap a CpG island.

CONCLUSION

Our assessment of six PPPs on 1.06 x 109 bps of human genome sequence reveals the specific strengths and weaknesses of individual PPPs. Our meta-predictor outperforms any individual PPP in sensitivity and specificity. Furthermore, we discovered that the 5' alternative promoters are more likely to be associated with a CpG island.

Boosting with stumps for predicting transcription start sites

CoreBoost applies a boosting technique to select important features for predicting core promoters with diverse patterns.

Promoter prediction is a difficult but important problem in gene finding, and it is critical for elucidating the regulation of gene expression. We introduce a new promoter prediction program, CoreBoost, which applies a boosting technique with stumps to select important small-scale as well as large-scale features. CoreBoost improves greatly on locating transcription start sites. We also demonstrate that by further utilizing some tissue-specific information, better accuracy can be achieved.

Multiple non-collinear TF-map alignments of promoter regions

Background

The analysis of the promoter sequence of genes with similar expression patterns is a basic tool to annotate common regulatory elements. Multiple sequence alignments are on the basis of most comparative approaches. The characterization of regulatory regions from co-expressed genes at the sequence level, however, does not yield satisfactory results in many occasions as promoter regions of genes sharing similar expression programs often do not show nucleotide sequence conservation.

Results

In a recent approach to circumvent this limitation, we proposed to align the maps of predicted transcription factors (referred as TF-maps) instead of the nucleotide sequence of two related promoters, taking into account the label of the corresponding factor and the position in the primary sequence. We have now extended the basic algorithm to permit multiple promoter comparisons using the progressive alignment paradigm. In addition, non-collinear conservation blocks might now be identified in the resulting alignments. We have optimized the parameters of the algorithm in a small, but well-characterized collection of human-mouse-chicken-zebrafish orthologous gene promoters.

Conclusion

Results in this dataset indicate that TF-map alignments are able to detect high-level regulatory conservation at the promoter and the 3'UTR gene regions, which cannot be detected by the typical sequence alignments. Three particular examples are introduced here to illustrate the power of the multiple TF-map alignments to characterize conserved regulatory elements in absence of sequence similarity. We consider this kind of approach can be extremely useful in the future to annotate potential transcription factor binding sites on sets of co-regulated genes from high-throughput expression experiments.

Context specific transcription factor prediction

One of the goals of systems biology is the identification of regulatory mechanisms that govern an organism's response to external stimuli. Transcription factors have been hypothesized as a major contributor to an organism's response to various outside stimuli, and a great deal of work has been done to predict the set of transcription factors which regulate a given gene. Most of the current methods seek to identify possible binding sites from genomic sequence. Initial attempts at predicting transcription factors from genomic sequences suffered from the problem of false positives. Making the problem more difficult, it has also been shown that while predicted binding sites might be false positives, they can be shown to bind to their corresponding sequences in vitro. One method for rectifying this is through the use of phylogenetic analysis in which only regions which show high evolutionary conservation are analyzed. However such an approach may be too stringent because of the level of degeneracy shown in transcription factor binding site position weight matrices. Due to the degeneracy, there may be only a few bases that need to be conserved across species. Therefore, while a sequence may not show a high level of evolutionary conservation, these sequences may still show high affinity for the same transcription factor. In predicting transcription factor binding we explore the notion that "Co-expression implies co-regulation" [Allocco et al. BMC Bioinformatics 5:18, 2004]. With multiple genes requiring similar transcription factors binding sites, there exists a basis for eliminating false positives. This method allows for the selection of transcription factors binding sites that are active under a given experimental paradigm, thereby allowing us to indirectly incorporate the effects of chromosome and recognition site presentation upon transcription factor binding prediction. Rather than having to rationalize that a few transcription factors binding sites are over-represented in a cluster of genes, one can show that a few transcription factors are active in the cluster of genes that have been grouped together. Although the method focuses on predicting experiment-specific transcription factor binding sites, it is possible that if such a methodology were used in an iterative process where different experiments were analyzed, one could obtain a comprehensive set of transcription factors binding sites which regulate the various dynamic responses shown by biological systems under a variety of conditions hence building a more comprehensive model of transcriptional regulation.

TRED: a transcriptional regulatory element database, new entries and other development

Transcriptional factors (TFs) and many of their target genes are involved in gene regulation at the level of transcription. To decipher gene regulatory networks (GRNs) we require a comprehensive and accurate knowledge of transcriptional regulatory elements. TRED (http://rulai.cshl.edu/TRED) was designed as a resource for gene regulation and function studies. It collects mammalian cis- and trans-regulatory elements together with experimental evidence. All the regulatory elements were mapped on to the assembled genomes. In this new release, we included a total of 36 TF families involved in cancer. Accordingly, the number of target promoters and genes for TF families has increased dramatically. There are 11,660 target genes (7479 in human, 2691 in mouse and 1490 in rat) and 14,908 target promoters (10,225 in human, 2985 in mouse and 1698 in rat). Additionally, we constructed GRNs for each TF family by connecting the TF-target gene pairs. Such interaction data between TFs and their target genes will assist detailed functional studies and help to obtain a panoramic view of the GRNs for cancer research.

Statistical significance of cis-regulatory modules

BACKGROUND

It is becoming increasingly important for researchers to be able to scan through large genomic regions for transcription factor binding sites or clusters of binding sites forming cis-regulatory modules. Correspondingly, there has been a push to develop algorithms for the rapid detection and assessment of cis-regulatory modules. While various algorithms for this purpose have been introduced, most are not well suited for rapid, genome scale scanning.

RESULTS

We introduce methods designed for the detection and statistical evaluation of cis-regulatory modules, modeled as either clusters of individual binding sites or as combinations of sites with constrained organization. In order to determine the statistical significance of module sites, we first need a method to determine the statistical significance of single transcription factor binding site matches. We introduce a straightforward method of estimating the statistical significance of single site matches using a database of known promoters to produce data structures that can be used to estimate p-values for binding site matches. We next introduce a technique to calculate the statistical significance of the arrangement of binding sites within a module using a max-gap model. If the module scanned for has defined organizational parameters, the probability of the module is corrected to account for organizational constraints. The statistical significance of single site matches and the architecture of sites within the module can be combined to provide an overall estimation of statistical significance of cis-regulatory module sites.

CONCLUSION

The methods introduced in this paper allow for the detection and statistical evaluation of single transcription factor binding sites and cis-regulatory modules. The features described are implemented in the Search Tool for Occurrences of Regulatory Motifs (STORM) and MODSTORM software.

Tissue-specific regulatory elements in mammalian promoters

Transcription factor-binding sites and the cis-regulatory modules they compose are central determinants of gene expression. We previously showed that binding site motifs and modules in proximal promoters can be used to predict a significant portion of mammalian tissue-specific transcription. Here, we report on a systematic analysis of promoters controlling tissue-specific expression in heart, kidney, liver, pancreas, skeletal muscle, testis and CD4 T cells, for both human and mouse. We integrated multiple sources of expression data to compile sets of transcripts with strong evidence for tissue-specific regulation. The analysis of the promoters corresponding to these sets produced a catalog of predicted tissue-specific motifs and modules, and cis-regulatory elements. Predicted regulatory interactions are supported by statistical evidence, and provide a foundation for targeted experiments that will improve our understanding of tissue-specific regulatory networks. In a broader context, methods used to construct the catalog provide a model for the analysis of genomic regions that regulate differentially expressed genes.

A scaffold for the Chinese hamster genome

Chinese hamster ovary (CHO) cells are a prevalent tool in biological research and are among the most widely used host cell lines for production of recombinant therapeutic proteins. While research in other organisms has been revolutionized through the development of DNA sequence-based tools, the lack of comparable genomic resources for the Chinese hamster has impeded similar work in CHO cell lines. A comparative genomics approach, based upon the completely sequenced mouse genome, can facilitate genomic work in this important organism. Using chromosome synteny to define regions of conserved linkage between Chinese hamster and mouse chromosomes, a working scaffold for the Chinese hamster genome has been developed. Mapping CHO and Chinese hamster sequences to the mouse genome creates direct access to relevant information in public databases. Additionally, mapping gene expression data onto a chromosome scaffold affords the ability to interpret information in a genomic context, potentially revealing important structural and regulatory features in the Chinese hamster genome. Further development of this genomic scaffold will provide opportunities to use biomolecular tools for research in CHO cell lines today and will be an asset to future efforts to sequence the Chinese hamster genome.

Unraveling transcriptional control and cis-regulatory codes using the software suite GeneACT

Deciphering gene regulatory networks requires the systematic identification of functional cis-acting regulatory elements. We present a suite of web-based bioinformatics tools, called GeneACT http://promoter.colorado.edu, that can rapidly detect evolutionarily conserved transcription factor binding sites or microRNA target sites that are either unique or over-represented in differentially expressed genes from DNA microarray data. GeneACT provides graphic visualization and extraction of common regulatory sequence elements in the promoters and 3'-untranslated regions that are conserved across multiple mammalian species.

Adaptively inferring human transcriptional subnetworks

Although the human genome has been sequenced, progress in understanding gene regulation in humans has been particularly slow. Many computational approaches developed for lower eukaryotes to identify cis-regulatory elements and their associated target genes often do not generalize to mammals, largely due to the degenerate and interactive nature of such elements. Motivated by the switch-like behavior of transcriptional responses, we present a systematic approach that allows adaptive determination of active transcriptional subnetworks (cis-motif combinations, the direct target genes and physiological processes regulated by the corresponding transcription factors) from microarray data in mammals, with accuracy similar to that achieved in lower eukaryotes. Our analysis uncovered several new subnetworks active in human liver and in cell-cycle regulation, with similar functional characteristics as the known ones. We present biochemical evidence for our predictions, and show that the recently discovered G2/M-specific E2F pathway is wider than previously thought; in particular, E2F directly activates certain mitotic genes involved in hepatocellular carcinomas. Additionally, we demonstrate that this method can predict subnetworks in a condition-specific manner, as well as regulatory crosstalk across multiple tissues. Our approach allows systematic understanding of how phenotypic complexity is regulated at the transcription level in mammals and offers marked advantage in systems where little or no prior knowledge of transcriptional regulation is available.

A clustering property of highly-degenerate transcription factor binding sites in the mammalian genome

Transcription factor binding sites (TFBSs) are short DNA sequences interacting with transcription factors (TFs), which regulate gene expression. Due to the relatively short length of such binding sites, it is largely unclear how the specificity of protein–DNA interaction is achieved. Here, we have performed a genome-wide analysis of TFBS-like sequences for the transcriptional repressor, RE1 Silencing Transcription Factor (REST), as well as for several other representative mammalian TFs (c-myc, p53, HNF-1 and CREB). We find a nonrandom distribution of inexact sites for these TFs, referred to as highly-degenerate TFBSs, that are enriched around the cognate binding sites. Comparisons among human, mouse and rat orthologous promoters reveal that these highly-degenerate sites are conserved significantly more than expected by random chance, suggesting their positive selection during evolution. We propose that this arrangement provides a favorable genomic landscape for functional target site selection.