TFM-Explorer: mining cis-regulatory regions in genomes

Identification of upstream transcription factor binding sites in orthologous genes using mixed Student’s t-test statistics

Background

Transcription factor (TF) regulates the transcription of DNA to messenger RNA by binding to upstream sequence motifs. Identifying the locations of known motifs in whole genomes is computationally intensive.

Methodology/Principal findings

This study presents a computational tool, named “Grit”, for screening TF-binding sites (TFBS) by coordinating transcription factors to their promoter sequences in orthologous genes. This tool employs a newly developed mixed Student’s t-test statistical method that detects high-scoring binding sites utilizing conservation information among species. The program performs sequence scanning at a rate of 3.2 Mbp/s on a quad-core Amazon server and has been benchmarked by the well-established ChIP-Seq datasets, putting Grit amongst the top-ranked TFBS predictors. It significantly outperforms the well-known transcription factor motif scanning tools, Pscan (4.8%) and FIMO (17.8%), in analyzing well-documented ChIP-Atlas human genome Chip-Seq datasets.

Significance

Grit is a good alternative to current available motif scanning tools.

Locating transcription factor-binding (TF-binding) site in the genome and identification their function is fundamental in understanding various biological processes. Improve the performance of the prediction tools is important because accurate TF-binding site prediction can save cost and time for wet-lab experiments. Also, genome wide TF-binding site prediction can provide new insights for transcriptome regulation in system biology perspective. This study developed a new TF-binding site prediction tool based on mixed Student’s t-test statistical method. The tool is amongst the top-ranked TF-binding site predictors, as such, it can help the researchers in TF-binding site identification and transcriptional regulation mechanism interpretation of genes.

Network mapping of primary CD34+ cells by Ampliseq based whole transcriptome targeted resequencing identifies unexplored differentiation regulatory relationships

With the exception of a few master transcription factors, regulators of neutrophil maturation are poorly annotated in the intermediate phenotypes between the granulocyte-macrophage progenitor (GMP) and the mature neutrophil phenotype. Additional challenges in identifying gene expression regulators in differentiation pathways relate to challenges wherein starting cell populations are heterogeneous in lineage potential and development, are spread across various states of quiescence, as well as sample quality and input limitations. These factors contribute to data variability make it difficult to draw simple regulatory inferences. In response we have applied a multi-omics approach using primary blood progenitor cells primed for homogeneous proliferation and granulocyte differentiation states which combines whole transcriptome resequencing (Ampliseq RNA) supported by droplet digital PCR (ddPCR) validation and mass spectrometry-based proteomics in a hypothesis-generation study of neutrophil differentiation pathways. Primary CD34+ cells isolated from human cord blood were first precultured in non-lineage driving medium to achieve an active, proliferating phenotype from which a neutrophil primed progenitor was isolated and cultured in neutrophil lineage supportive medium. Samples were then taken at 24-hour intervals over 9 days and analysed by Ampliseq RNA and mass spectrometry. The Ampliseq dataset depth, breadth and quality allowed for several unexplored transcriptional regulators and ncRNAs to be identified using a combinatorial approach of hierarchical clustering, enriched transcription factor binding motifs, and network mapping. Network mapping in particular increased comprehension of neutrophil differentiation regulatory relationships by implicating ARNT, NHLH1, PLAG1, and 6 non-coding RNAs associated with PU.1 regulation as cell-engineering targets with the potential to increase total neutrophil culture output. Overall, this study develops and demonstrates an effective new hypothesis generation methodology for transcriptome profiling during differentiation, thereby enabling identification of novel gene targets for editing interventions.

Identification and Functional Annotation of Genes Related to Horses' Performance: From GWAS to Post-GWAS

Simple Summary

It is assumed that the athletic performance of horses is influenced by a large number of genes; however, to date, not many genomic studies have been performed to identify candidate genes. In this study we performed a systematic review of genome-wide association studies followed by functional analyses aiming to identify the most candidate genes for horse performance. We were successful in identifying 669 candidate genes, from which we built biological process networks. Regulatory elements (transcription factors, TFs) of these genes were identified and used to build a gene–TF network. Genes and TFs presented in this study are suggested to play a role in the studied traits through biological processes related with exercise performance, for example, positive regulation of glucose metabolism, regulation of vascular endothelial growth factor production, skeletal system development, cellular response to fatty acids and cellular response to lipids. In general, this study may provide insights into the genetic architecture underlying horse performance in different breeds around the world.

Abstract

Integration of genomic data with gene network analysis can be a relevant strategy for unraveling genetic mechanisms. It can be used to explore shared biological processes between genes, as well as highlighting transcription factors (TFs) related to phenotypes of interest. Unlike other species, gene–TF network analyses have not yet been well applied to horse traits. We aimed to (1) identify candidate genes associated with horse performance via systematic review, and (2) build biological processes and gene–TF networks from the identified genes aiming to highlight the most candidate genes for horse performance. Our systematic review considered peer-reviewed articles using 20 combinations of keywords. Nine articles were selected and placed into groups for functional analysis via gene networks. A total of 669 candidate genes were identified. From that, gene networks of biological processes from each group were constructed, highlighting processes associated with horse performance (e.g., regulation of systemic arterial blood pressure by vasopressin and regulation of actin polymerization and depolymerization). Transcription factors associated with candidate genes were also identified. Based on their biological processes and evidence from the literature, we identified the main TFs related to horse performance traits, which allowed us to construct a gene–TF network highlighting TFs and the most candidate genes for horse performance.

RNA-seq reveals downregulated osteochondral genes potentially related to tibia bacterial chondronecrosis with osteomyelitis in broilers

Background

Bacterial chondronecrosis with osteomyelitis (BCO) develops in the growth plate (GP) of the proximal femur and tibia and is initiated by damage to the less mineralized chondrocytes followed by colonization of opportunistic bacteria. This condition affects approximately 1% of all birds housed, being considered one of the major causes of lameness in fast growing broilers. Although several studies have been previously performed aiming to understand its pathogenesis, the molecular mechanisms involved with BCO remains to be elucidated. Therefore, this study aimed to generate a profile of global differential gene expression involved with BCO in the tibia of commercial broilers, through RNA sequencing analysis to identity genes and molecular pathways involved with BCO in chickens.

Results

Our data showed 192 differentially expressed (DE) genes: 63 upregulated and 129 downregulated in the GP of the tibia proximal epiphysis of BCO-affected broilers. Using all DE genes, six Biological Processes (BP) were associated with bone development (connective tissue development, cartilage development, skeletal system development, organ morphogenesis, system development and skeletal system morphogenesis). The analyses of the upregulated genes did not indicate any significant BP (FDR < 0.05). However, with the downregulated genes, the same BP were identified when using all DE genes in the analysis, with a total of 26 coding genes explaining BCO in the tibia: ACAN, ALDH1A2, CDH7, CHAD, CHADL, COL11A1, COMP, CSGALNACT1, CYR61, FRZB, GAL3ST1, HAPLN1, IHH, KIF26B, LECT1, LPPR1, PDE6B, RBP4A, SERINC5, SFRP1, SOX8, SOX9, TENM2, THBS1, UCHL1 and WFIKKN2. In addition, seven transcription factors were also associated to BCO: NFATC2, MAFB, HIF1A-ARNT, EWSR1-FLI1, NFIC, TCF3 and NF-KAPPAB.

Conclusions

Our data show that osteochondral downregulated genes are potential molecular causes of BCO in broilers, and the bacterial process seems to be, in fact, a secondary condition. Sixteen genes responsible for bone and cartilage formation were downregulated in BCO-affected broilers being strong candidate genes to trigger this disorder.

A system view and analysis of essential hypertension

OBJECTIVES

The goal of this study was to investigate genes associated with essential hypertension from a system perspective, making use of bioinformatic tools to gain insights that are not evident when focusing at a detail-based resolution.

METHODS

Using various databases (pathways, Genome Wide Association Studies, knockouts etc.), we compiled a set of about 200 genes that play a major role in hypertension and identified the interactions between them. This enabled us to create a protein-protein interaction network graph, from which we identified key elements, based on graph centrality analysis. Enriched gene regulatory elements (transcription factors and microRNAs) were extracted by motif finding techniques and knowledge-based tools.

RESULTS

We found that the network is composed of modules associated with functions such as water retention, endothelial vasoconstriction, sympathetic activity and others. We identified the transcription factor SP1 and the two microRNAs miR27 (a and b) and miR548c-3p that seem to play a major role in regulating the network as they exert their control over several modules and are not restricted to specific functions. We also noticed that genes involved in metabolic diseases (e.g. insulin) are central to the network.

CONCLUSION

We view the blood-pressure regulation mechanism as a system-of-systems, composed of several contributing subsystems and pathways rather than a single module. The system is regulated by distributed elements. Understanding this mode of action can lead to a more precise treatment and drug target discovery. Our analysis suggests that insulin plays a primary role in hypertension, highlighting the tight link between essential hypertension and diseases associated with the metabolic syndrome.

Evolution of the human cold/menthol receptor, TRPM8

Genes showing versatile functions or subjected to fast expansion and contraction during the adaptation of species to specific ecological conditions, like sensory receptors for odors, pheromones and tastes, are characterized by a great plasticity through evolution. One of the most fascinating sensory receptors in the family of TRP channels, the cold and menthol receptor TRPM8, has received significant attention in the literature. Recent studies have reported the existence of TRPM8 channel isoforms encoded by alternative mRNAs transcribed from alternative promoters and processed by alternative splicing. Since the first draft of the human genome was accomplished in 2000, alternative transcription, alternative splicing and alternative translation have appeared as major sources of gene product diversity and are thought to participate in the generation of complexity in higher organisms. In this study, we investigate whether alternative transcription has been a driving force in the evolution of the human forms of the cold receptor TRPM8. We identified 33 TRPM8 alternative mRNAs (24 new sequences) and their associated protein isoforms in human tissues. Using comparative genomics, we described the evolution of the human TRPM8 sequences in eight ancestors since the origin of Amniota, and estimated in which ancestors the new TRPM8 variants originated. In order to validate the estimated origins of this receptor, we performed experimental validations of predicted exons in mouse tissues. Our results suggest a first diversification event of the cold receptor in the Boreoeutheria ancestor, and a subsequent divergence at the origin of Simiiformes.

Screening of Tumor Suppressor Genes in Metastatic Colorectal Cancer

Most tumor suppressor genes are commonly inactivated in the development of colorectal cancer (CRC). The activation of tumor suppressor genes may be beneficial to suppress the development and metastasis of CRC. This study analyzed genes expression and methylation levels in different stages of CRC. Genes with downregulated mRNA expression and upregulated methylation level in advanced CRC were screened as the potential tumor suppressor genes. After comparing the methylation level of screened genes, we found that MBD1 gene had downregulated mRNA expression and upregulated methylation levels in advanced CRC and continuously upregulated methylation level in the progression of CRC. Enrichment analysis revealed that genes expression in accordance with the elevated expression of MBD1 mainly located on chromosomes 17p13 and 17p12 and 8 tumor suppressor genes located on chromosome 17p13. Further enrichment analysis of transcription factor binding site identified that SP1 binding site had higher enrichment and could bind with MBD1. In conclusion, MBD1 may be a tumor suppressor gene in advanced CRC and affect the development and metastasis of CRC by regulating 8 tumor suppressor genes through binding with SP1.

HSD1 and AQP7 short-term gene regulation by cortisone in 3T3-L1 adipocytes

Adipose Tissue (AT) is a complex organ with a crucial regulatory role in energy metabolism and in the development of obesity and the Metabolic Syndrome (MS). Modified responses and the metabolism of hormones have been observed in visceral adiposity during obesity, specifically as related with cortisone. The objective of this study was to assess, in the 3T3-L1 adipocyte cell line, the short-term effect of cortisone on the expression of 11β-Hydroxysteroid dehydrogenase 1 (Hsd1), which is responsible for activation of cortisone into cortisol, and for Aquaporin 7 (Aqp7), involved in glycerol transport through the cell membrane. Total RNA (tRNA) and complementary DNA (cDNA) were obtained from cell samples treated with cortisone (0.1, 1, and 10 μM) during different times (0, 5, 10, 15, and 20 min, and 48 h) to quantify the expression of the aforementioned genes by real time PCR employing MnSOD and Ppia as housekeeping genes. There was a time-dependent response of Aqp7, a dose-dependent response of Hsd1, and an increase observed in the expression of both genes during min 1 of treatment (5- and 6-fold, respectively), followed by a decrease during the following 5-10 min (P < 0.05). With the 1-μM cortisone treatment, both genes showed cubic tendencies in their expression; the Hsd1 tendency is described by the equation y = 0.18×(3)-1.65×(2)+3.59x+1.31, while the Aqp7 tendency is described by y = 0.33×(3)-2.67×(2)+4.93x+1.84. There are immediate and quantitatively important actions of cortisone on the expression of Aqp7 and Hsd1 in 3T3-L1 adipocytes.

Bisphenol A affects placental layers morphology and angiogenesis during early pregnancy phase in mice

Bisphenol A (BPA) is a widespread endocrine disrupter mainly used in food contact plastics. Much evidence supports the adverse effects of BPA, particularly on susceptible groups such as pregnant women. The present study considered placental development - relevant for pregnancy outcomes and fetal nutrition/programming - as a potential target of BPA. Pregnant CD-1 mice were administered per os with vehicle, 0.5 (BPA05) or 50 mg kg(-1) (BPA50) body weight day(-1) of BPA, from gestational day (GD) 1 to GD11. At GD12, BPA50 induced significant degeneration and necrosis of giant cells, increased vacuolization in the junctional zone in the absence of glycogen accumulation and reduction of the spongiotrophoblast layer. In addition, BPA05 induced glycogen depletion as well as significant nuclear accumulation of β-catenin in trophoblasts of labyrinthine and spongiotrophoblast layers, supporting the activation of the Wnt/β-catenin pathway. Transcriptomic analysis indicated that BPA05 promoted and BPA50 inhibited blood vessel development and branching; morphologically, maternal vessels were narrower in BPA05 placentas, whereas embryonic and maternal vessels were irregularly dilated in the labyrinth of BPA50 placentas. Quantitative polymerase chain reaction evidenced an estrogen receptor β induction by BPA50, which did not correspond to downstream genes activation; indeed, the transcription factor binding sites analysis supported the AhR/Arnt complex as regulator of BPA50-modulated genes. Conversely, Creb appeared as the main transcription factor regulating BPA05-modulated genes. Embryonic structures (head, forelimb) showed divergent perturbations upon BPA05 or BPA50 exposure, potentially related to unbalanced embryonic nutrition and/or to modulation of genes involved in embryo development. Our findings support placenta as an important target of BPA, even at environmentally relevant dose levels.

The Human Blood Metabolome-Transcriptome Interface

Biological systems consist of multiple organizational levels all densely interacting with each other to ensure function and flexibility of the system. Simultaneous analysis of cross-sectional multi-omics data from large population studies is a powerful tool to comprehensively characterize the underlying molecular mechanisms on a physiological scale. In this study, we systematically analyzed the relationship between fasting serum metabolomics and whole blood transcriptomics data from 712 individuals of the German KORA F4 cohort. Correlation-based analysis identified 1,109 significant associations between 522 transcripts and 114 metabolites summarized in an integrated network, the 'human blood metabolome-transcriptome interface' (BMTI). Bidirectional causality analysis using Mendelian randomization did not yield any statistically significant causal associations between transcripts and metabolites. A knowledge-based interpretation and integration with a genome-scale human metabolic reconstruction revealed systematic signatures of signaling, transport and metabolic processes, i.e. metabolic reactions mainly belonging to lipid, energy and amino acid metabolism. Moreover, the construction of a network based on functional categories illustrated the cross-talk between the biological layers at a pathway level. Using a transcription factor binding site enrichment analysis, this pathway cross-talk was further confirmed at a regulatory level. Finally, we demonstrated how the constructed networks can be used to gain novel insights into molecular mechanisms associated to intermediate clinical traits. Overall, our results demonstrate the utility of a multi-omics integrative approach to understand the molecular mechanisms underlying both normal physiology and disease.

The role of melanogenesis in regulation of melanoma behavior: melanogenesis leads to stimulation of HIF-1α expression and HIF-dependent attendant pathways

To study the effect of melanogenesis on HIF-1α expression and attendant pathways, we used stable human and hamster melanoma cell lines in which the amelanotic vs. melanotic phenotypes are dependent upon the concentration of melanogenesis precursors in the culture media. The induction of melanin pigmentation led to significant up-regulation of HIF-1α, but not HIF-2α, protein in melanized cells for both lines. Similar upregulation of nuclear HIF-1α was observed in excisions of advanced melanotic vs. amelanotic melanomas. In cultured cells, melanogenesis also significantly stimulated expression of classical HIF-1-dependent target genes involved in angiogenesis and cellular metabolism, including glucose metabolism and stimulation of activity of key enzymes in the glycolytic pathway. Several other stress related genes containing putative HRE consensus sites were also upregulated by melanogenesis, concurrently with modulation of expression of HIF-1-independent genes encoding for steroidogenic enzymes, cytokines and growth factors. Immunohistochemical studies using a large panel of pigmented lesions revealed that higher levels of HIF-1α and GLUT-1 were detected in advanced melanomas in comparison to melanocytic nevi or thin melanomas localized to the skin. However, the effects on overall or disease free survival in melanoma patients were modest or absent for GLUT-1 or for HIF-1α, respectively. In conclusion, induction of the melanogenic pathway leads to robust upregulation of HIF-1-dependent and independent pathways in cultured melanoma cells, suggesting a key role for melanogenesis in regulation of cellular metabolism.

Putative role of HIF transcriptional activity in melanocytes and melanoma biology

Hypoxia-inducible factor-1α (HIF-1α) is a highly oxygen sensitive bHLH protein that is part of the heterodimeric HIF-1 transcription factor. Under hypoxic stress, HIF-1 activity is induced to control expression of multiple downstream target genes, including vascular endothelial growth factor (VEGF). The normal epidermis exists in a constant mild hypoxic microenvironment and constitutively expresses HIF-1α and HIF-2α. Expression of HIF-1α and/or HIF-2α has been suggested to correlate with the increased malignant potential of melanocytes, therefore, failures of melanoma therapies may be partially linked to high HIF activity. Notably, melanomas that have the V600E BRAF mutation exhibit increased HIF-1α expression. We have utilized a bioinformatics approach to identify putative hypoxia response elements (HREs) in a set of genes known to participate in the process of melanogenesis (includingTRPM1, SLC45A2, HRAS, C-KIT, PMEL and CRH). While some of the mechanistic links between these genes and the HIF pathway have been previously explored, others await further investigation. Although agents targeting HIF activity have been proposed as novel treatment modalities for melanoma, there are currently no clinical trials in progress to test their efficacy in melanoma.

Screening and regulatory network analysis of survival-related genes of patients with colorectal cancer

The purpose of this study was to screen key survival-related genes from patients with colorectal cancer and explore signal transduction network of the involved genes. In a previous study, survival-related genes of patients with colorectal cancer were selected by colorectal cancer-related expression data GSE17538 using the Significance Analysis of Microarrays (SAM3.01) software, and 235 genes related to the survival of patients with colorectal cancer were obtained. Therefore, the following screening and analysis were conducted on these 235 genes in this study. First, the enrichment analysis of transcription factor binding sites was conducted on the 235 genes. Genes with more than seven transcription factor binding sites were screened. Then, these genes and upregulated genes in colorectal cancer were intersected. Finally, survival analysis and regulatory network analysis were conducted on the screened genes. This allowed clarification of the relationship between these genes and the survival of patients with colorectal cancer and the signaling network involving these genes in the cell signal transduction network of colorectal cancer. Through the above analysis, six upregulated genes in colorectal cancer related to the survival of colorectal cancer patients and highly regulated by transcription factors were selected, namely STX2, PODXL, KLK6, GRB10, EHBP1 and CREB5. These genes are involved in signal regulatory networks related to colorectal cancer metastasis-related signaling pathways. Therefore, the survival of patients with colorectal cancer is closely correlated with colorectal cancer metastasis. The six survival-related genes affect the survival of patients by regulating colorectal cancer metastasis-associated signaling pathways.

Landscape of transcriptional deregulations in the preeclamptic placenta

Preeclampsia is a pregnancy disease affecting 5 to 8% of pregnant women and a leading cause of both maternal and fetal mortality and morbidity. Because of a default in the process of implantation, the placenta of preeclamptic women undergoes insufficient vascularization. This results in placental ischemia, inflammation and subsequent release of placental debris and vasoactive factors in the maternal circulation causing a systemic endothelial activation. Several microarray studies have analyzed the transcriptome of the preeclamptic placentas to identify genes which could be involved in placental dysfunction. In this study, we compared the data from publicly available microarray analyses to obtain a consensus list of modified genes. This allowed to identify consistently modified genes in the preeclamptic placenta. Of these, 67 were up-regulated and 31 down-regulated. Assuming that changes in the transcription level of co-expressed genes may result from the coordinated action of a limited number of transcription factors, we looked for over-represented putative transcription factor binding sites in the promoters of these genes. Indeed, we found that the promoters of up-regulated genes are enriched in putative binding sites for NFkB, CREB, ANRT, REEB1, SP1, and AP-2. In the promoters of down-regulated genes, the most prevalent putative binding sites are those of MZF-1, NFYA, E2F1 and MEF2A. These transcriptions factors are known to regulate specific biological pathways such as cell responses to inflammation, hypoxia, DNA damage and proliferation. We discuss here the molecular mechanisms of action of these transcription factors and how they can be related to the placental dysfunction in the context of preeclampsia.

D-Light on promoters: a client-server system for the analysis and visualization of cis-regulatory elements

BACKGROUND

The binding of transcription factors to DNA plays an essential role in the regulation of gene expression. Numerous experiments elucidated binding sequences which subsequently have been used to derive statistical models for predicting potential transcription factor binding sites (TFBS). The rapidly increasing number of genome sequence data requires sophisticated computational approaches to manage and query experimental and predicted TFBS data in the context of other epigenetic factors and across different organisms.

RESULTS

We have developed D-Light, a novel client-server software package to store and query large amounts of TFBS data for any number of genomes. Users can add small-scale data to the server database and query them in a large scale, genome-wide promoter context. The client is implemented in Java and provides simple graphical user interfaces and data visualization. Here we also performed a statistical analysis showing what a user can expect for certain parameter settings and we illustrate the usage of D-Light with the help of a microarray data set.

CONCLUSIONS

D-Light is an easy to use software tool to integrate, store and query annotation data for promoters. A public D-Light server, the client and server software for local installation and the source code under GNU GPL license are available at http://biwww.che.sbg.ac.at/dlight.

Proteomic profiling in Drosophila reveals potential Dube3a regulation of the actin cytoskeleton and neuronal homeostasis

The molecular defects associated with Angelman syndrome (AS) and 15q duplication autism are directly correlated to expression levels of the E3 ubiquitin ligase protein UBE3A. Here we used Drosophila melanogaster to screen for the targets of this ubiquitin ligase under conditions of both decreased (as in AS) or increased (as in dup(15)) levels of the fly Dube3a or human UBE3A proteins. Using liquid phase isoelectric focusing of proteins from whole fly head extracts we identified a total of 50 proteins that show changes in protein, and in some cases transcriptional levels, when Dube3a fluctuates. We analyzed head extracts from cytoplasmic, nuclear and membrane fractions for Dube3a regulated proteins. Our results indicate that Dube3a is involved in the regulation of cellular functions related to ATP synthesis/metabolism, actin cytoskeletal integrity, both catabolism and carbohydrate metabolism as well as nervous system development and function. Sixty-two percent of the proteins were >50% identical to homologous human proteins and 8 have previously be shown to be ubiquitinated in the fly nervous system. Eight proteins may be regulated by Dube3a at the transcript level through the transcriptional co-activation function of Dube3a. We investigated one autism-associated protein, ATPα, and found that it can be ubiquitinated in a Dube3a dependent manner. We also found that Dube3a mutants have significantly less filamentous actin than wild type larvae consistent with the identification of actin targets regulated by Dube3a. The identification of UBE3A targets is the first step in unraveling the molecular etiology of AS and duplication 15q autism.

CBS: an open platform that integrates predictive methods and epigenetics information to characterize conserved regulatory features in multiple Drosophila genomes

Background

Information about the composition of regulatory regions is of great value for designing experiments to functionally characterize gene expression. The multiplicity of available applications to predict transcription factor binding sites in a particular locus contrasts with the substantial computational expertise that is demanded to manipulate them, which may constitute a potential barrier for the experimental community.

Results

CBS (Conserved regulatory Binding Sites, http://compfly.bio.ub.es/CBS) is a public platform of evolutionarily conserved binding sites and enhancers predicted in multiple Drosophila genomes that is furnished with published chromatin signatures associated to transcriptionally active regions and other experimental sources of information. The rapid access to this novel body of knowledge through a user-friendly web interface enables non-expert users to identify the binding sequences available for any particular gene, transcription factor, or genome region.

Conclusions

The CBS platform is a powerful resource that provides tools for data mining individual sequences and groups of co-expressed genes with epigenomics information to conduct regulatory screenings in Drosophila.

Caffeic acid phenethyl ester suppresses the proliferation of human prostate cancer cells through inhibition of p70S6K and Akt signaling networks

Caffeic acid phenethyl ester (CAPE) is a bioactive component derived from honeybee hive propolis. CAPE has been shown to have antimitogenic, anticarcinogenic, and other beneficial medicinal properties. Many of its effects have been shown to be mediated through its inhibition of NF-κB signaling pathways. We took a systematic approach to uncover the effects of CAPE from hours to days on the signaling networks in human prostate cancer cells. We observed that CAPE dosage dependently suppressed the proliferation of LNCaP, DU-145, and PC-3 human prostate cancer cells. Administration of CAPE by gavage significantly inhibited the tumor growth of LNCaP xenografts in nude mice. Using LNCaP cells as a model system, we examined the effect of CAPE on gene expression, protein signaling, and transcriptional regulatory networks using micro-Western arrays and PCR arrays. We built a model of the impact of CAPE on cell signaling which suggested that it acted through inhibition of Akt-related protein signaling networks. Overexpression of Akt1 or c-Myc, a downstream target of Akt signaling, significantly blocked the antiproliferative effects of CAPE. In summary, our results suggest that CAPE administration may be useful as an adjuvant therapy for prostate and potentially other types of cancers that are driven by the p70S6K and Akt signaling networks.

Genome-wide profiling of methylation identifies novel targets with aberrant hypermethylation and reduced expression in low-risk myelodysplastic syndromes

Gene expression profiling signatures may be used to classify the subtypes of Myelodysplastic syndrome (MDS) patients. However, there are few reports on the global methylation status in MDS. The integration of genome-wide epigenetic regulatory marks with gene expression levels would provide additional information regarding the biological differences between MDS and healthy controls. Gene expression and methylation status were measured using high-density microarrays. A total of 552 differentially methylated CpG loci were identified as being present in low-risk MDS; hypermethylated genes were more frequent than hypomethylated genes. In addition, mRNA expression profiling identified 1005 genes that significantly differed between low-risk MDS and the control group. Integrative analysis of the epigenetic and expression profiles revealed that 66.7% of the hypermethylated genes were underexpressed in low-risk MDS cases. Gene network analysis revealed molecular mechanisms associated with the low-risk MDS group, including altered apoptosis pathways. The two key apoptotic genes BCL2 and ETS1 were identified as silenced genes. In addition, the immune response and micro RNA biogenesis were affected by the hypermethylation and underexpression of IL27RA and DICER1. Our integrative analysis revealed that aberrant epigenetic regulation is a hallmark of low-risk MDS patients and could have a central role in these diseases.

ReLA, a local alignment search tool for the identification of distal and proximal gene regulatory regions and their conserved transcription factor binding sites

Motivation: The prediction and annotation of the genomic regions involved in gene expression has been largely explored. Most of the energy has been devoted to the development of approaches that detect transcription start sites, leaving the identification of regulatory regions and their functional transcription factor binding sites (TFBSs) largely unexplored and with important quantitative and qualitative methodological gaps.

Results: We have developed ReLA (for REgulatory region Local Alignment tool), a unique tool optimized with the Smith–Waterman algorithm that allows local searches of conserved TFBS clusters and the detection of regulatory regions proximal to genes and enhancer regions. ReLA's performance shows specificities of 81 and 50% when tested on experimentally validated proximal regulatory regions and enhancers, respectively.

Availability: The source code of ReLA's is freely available and can be remotely used through our web server under http://www.bsc.es/cg/rela.

Contact:david.torrents@bsc.es

Supplementary information:Supplementary data are available at Bioinformatics online.

Transcriptomic coordination in the human metabolic network reveals links between n-3 fat intake, adipose tissue gene expression and metabolic health

Understanding the molecular link between diet and health is a key goal in nutritional systems biology. As an alternative to pathway analysis, we have developed a joint multivariate and network-based approach to analysis of a dataset of habitual dietary records, adipose tissue transcriptomics and comprehensive plasma marker profiles from human volunteers with the Metabolic Syndrome. With this approach we identified prominent co-expressed sub-networks in the global metabolic network, which showed correlated expression with habitual n-3 PUFA intake and urinary levels of the oxidative stress marker 8-iso-PGF(2α). These sub-networks illustrated inherent cross-talk between distinct metabolic pathways, such as between triglyceride metabolism and production of lipid signalling molecules. In a parallel promoter analysis, we identified several adipogenic transcription factors as potential transcriptional regulators associated with habitual n-3 PUFA intake. Our results illustrate advantages of network-based analysis, and generate novel hypotheses on the transcriptomic link between habitual n-3 PUFA intake, adipose tissue function and oxidative stress.

Inflammatory gene regulatory networks in amnion cells following cytokine stimulation: translational systems approach to modeling human parturition

A majority of the studies examining the molecular regulation of human labor have been conducted using single gene approaches. While the technology to produce multi-dimensional datasets is readily available, the means for facile analysis of such data are limited. The objective of this study was to develop a systems approach to infer regulatory mechanisms governing global gene expression in cytokine-challenged cells in vitro, and to apply these methods to predict gene regulatory networks (GRNs) in intrauterine tissues during term parturition. To this end, microarray analysis was applied to human amnion mesenchymal cells (AMCs) stimulated with interleukin-1β, and differentially expressed transcripts were subjected to hierarchical clustering, temporal expression profiling, and motif enrichment analysis, from which a GRN was constructed. These methods were then applied to fetal membrane specimens collected in the absence or presence of spontaneous term labor. Analysis of cytokine-responsive genes in AMCs revealed a sterile immune response signature, with promoters enriched in response elements for several inflammation-associated transcription factors. In comparison to the fetal membrane dataset, there were 34 genes commonly upregulated, many of which were part of an acute inflammation gene expression signature. Binding motifs for nuclear factor-κB were prominent in the gene interaction and regulatory networks for both datasets; however, we found little evidence to support the utilization of pathogen-associated molecular pattern (PAMP) signaling. The tissue specimens were also enriched for transcripts governed by hypoxia-inducible factor. The approach presented here provides an uncomplicated means to infer global relationships among gene clusters involved in cellular responses to labor-associated signals.