Dataset Information

Integrated analysis of motif activity and gene expression changes of transcription factors.

ABSTRACT: The ability to predict transcription factors based on sequence information in regulatory elements is a key step in systems-level investigation of transcriptional regulation. Here, we have developed a novel tool, IMAGE, for precise prediction of causal transcription factors based on transcriptome profiling and genome-wide maps of enhancer activity. High precision is obtained by combining a near-complete database of position weight matrices (PWMs), generated by compiling public databases and systematic prediction of PWMs for uncharacterized transcription factors, with a state-of-the-art method for PWM scoring and a novel machine learning strategy, based on both enhancers and promoters, to predict the contribution of motifs to transcriptional activity. We applied IMAGE to published data obtained during 3T3-L1 adipocyte differentiation and showed that IMAGE predicts causal transcriptional regulators of this process with higher confidence than existing methods. Furthermore, we generated genome-wide maps of enhancer activity and transcripts during human mesenchymal stem cell commitment and adipocyte differentiation and used IMAGE to identify positive and negative transcriptional regulators of this process. Collectively, our results demonstrate that IMAGE is a powerful and precise method for prediction of regulators of gene expression.

SUBMITTER: Madsen JGS

PROVIDER: S-EPMC5793788 | biostudies-literature | 2018 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Integrated analysis of motif activity and gene expression changes of transcription factors.

Madsen Jesper Grud Skat JGS Rauch Alexander A Van Hauwaert Elvira Laila EL Schmidt Søren Fisker SF Winnefeld Marc M Mandrup Susanne S

Genome research 20171212 2

The ability to predict transcription factors based on sequence information in regulatory elements is a key step in systems-level investigation of transcriptional regulation. Here, we have developed a novel tool, IMAGE, for precise prediction of causal transcription factors based on transcriptome profiling and genome-wide maps of enhancer activity. High precision is obtained by combining a near-complete database of position weight matrices (PWMs), generated by compiling public databases and syste ...[more]

PMID: 29233921

Similar Datasets

Project description:BackgroundThe identification of transcription factors (TFs) associated with a biological process is fundamental to understanding its regulatory mechanisms. From microarray data, however, the activity changes of TFs often cannot be directly observed due to their relatively low expression levels, post-transcriptional modifications, and other complications. Several approaches have been proposed to infer TF activity changes from microarray data. In some models, a linear relationship between gene expression and TF-gene binding strength is assumed. In some other models, the target genes of a TF are first determined by a significance cutoff to binding affinity scores, and then expression differentiation is checked between the target and other genes.ResultsWe propose a novel method, referred to as BASE (binding association with sorted expression), to infer TF activity changes from microarray expression profiles with the help of binding affinity data. It searches the maximum association between bind affinity profile of a TF and expression change profile along the direction of sorted differentiation. The method does not make hard target gene selection, rather, the significances of TF activity changes are evaluated by permutation tests of binding association at the end. To show the effectiveness of this method, we apply it to three typical examples using different kinds of binding affinity data, namely, ChIP-chip data, motif discovery data, and positional weighted matrix scanning data, respectively. The implications obtained from all three examples are consistent with established biological results. Moreover, the inferences suggest new and biological meaningful hypotheses for further investigation.ConclusionThe proposed method makes transcription inference from profiles of expression and binding affinity. The same machinery can be used to deal with various kinds of binding affinity data. The method does not require a linear assumption, and has the desirable property of scale-invariance with respect to TF-specific binding affinity. This method is easy to implement and can be routinely applied for transcriptional inferences in microarray studies.

Project description:Background. Papillary thyroid carcinoma (PTC) is the commonest thyroid malignancy originating from the follicle cells in the thyroid. Despite a good overall prognosis, certain high-risk cases as in those with lymph node metastasis (LNM) have progressive disease and poorer prognosis. MicroRNAs are a class of non-protein-coding, 19-24 nucleotides single-stranded RNAs which regulate gene expression and these molecules have been shown to play a role in LNM. The integrated analysis of miRNAs and gene expression profiles together with transcription factors (TFs) has been shown to improve the identification of functional miRNA-target gene-TF relationships, providing a more complete view of molecular events underlying metastasis process. Objectives. We reanalyzed The Cancer Genome Atlas (TCGA) datasets on PTC to identify differentially expressed miRNAs/genes in PTC patients with LNM-positive (LNM-P) versus lymph node negative (LNN) PTC patients and to investigate the miRNA-gene-TF regulatory circuit that regulate LNM in PTC. Results. PTC patients with LNM (PTC LNM-P) have a significantly shorter disease-free survival rate compared to PTC patients without LNM (PTC LNN) (Log-rank Mantel Cox test, p = 0.0049). We identified 181 significantly differentially expressed miRNAs in PTC LNM-P versus PTC LNN; 110 were upregulated and 71 were downregulated. The five topmost deregulated miRNAs were hsa-miR-146b, hsa-miR-375, hsa-miR-31, hsa-miR-7-2 and hsa-miR-204. In addition, 395 miRNAs were differentially expressed between PTC LNM-P and normal thyroid while 400 miRNAs were differentially expressed between PTC LNN and normal thyroid. We found four significant enrichment pathways potentially involved in metastasis to the lymph nodes, namely oxidative phosphorylation (OxPhos), cell adhesion molecules (CAMs), leukocyte transendothelial migration and cytokine-cytokine receptor interaction. OxPhos was the most significantly perturbed pathway (p = 4.70E-06) involving downregulation of 90 OxPhos-related genes. Significant interaction of hsa-miR-301b with HLF, HIF and REL/NFkB transcription factors were identified exclusively in PTC LNM-P versus PTC LNN. Conclusion. We found evidence of five miRNAs differentially expressed in PTC LNM-P. Alteration in OxPhos pathway could be the central event in metastasis to the lymph node in PTC. We postulate that hsa-miR-301b might be involved in regulating LNM in PTC via interactions with HLF, HIF and REL/NFkB. To the best of our knowledge, the roles of these TFs have been studied in PTC but the precise role of this miRNA with these TFs in LNM in PTC has not been investigated.

Project description:BackgroundDNA methylation is a fundamental component of epigenetic modification, which is intimately involved in the regulation of gene expression. One important DNA methylation pathway reduces the abilities of transcription factors to bind to gene promoter regions. Although many experiments have been designed to measure genome-wide DNA methylation levels at high resolution, the meaning of these different DNA methylation levels on transcription factor binding abilities remains poorly understood. We have, therefore, developed a method to quantitatively explore the extent to which DNA methylation levels can significantly reduce or even abolish the binding of certain transcription factors, resulting in reduced or non-expression of flanking genes. This method allows transcription factors that are functionally active in gene expression to be investigated.ResultsThe method is based on a general model that depicts the relationship between DNA methylation and transcription factor binding ability based on intrinsic component properties, and the model parameters can be optimized through relative analysis of recognized transcription factor binding status and gene expression profiling. With fixed models, transcription factors functionally active in the regulation of gene expression and affected by epigenetic DNA methylation can be identified and subsequently confirmed. The method identified eleven apparently functionally active transcriptional factors in SH-SY5Y neuroblastoma cells.ConclusionsCompared with gene regulatory elements, epigenetic modifications are able to change to dynamically respond to signals from physical, biological and social environments. Our proposed method is therefore designed to provide a dynamic assessment to investigate functionally active transcription factors. With the information deduced from our method, we can predict transcription factor binding status in promoter regions to further investigate how a particular gene is regulated by a specific group of transcription factors organized in a particular pattern. This will be helpful in the diagnosis and development of treatment for numerous diseases, including cancer. Although the method only investigates DNA methylation, it has the potential to be applied to more epigenetic factors, such as histone modification.

Project description:BackgroundEndometriosis is recognized as an estrogen-dependent inflammation disorder, estimated to affect 8%-15% of women of childbearing age. Currently, the etiology and pathogenesis of endometriosis are not completely clear. Underlying mechanism for endometriosis is still under debate and needs further exploration. The involvement of transcription factors and immune mediations may be involved in the pathophysiological process of endometriosis, but the specific mechanism remains to be explored. This study aims to investigate the underlying molecular mechanisms in endometriosis.MethodsThe gene expression profile of endometriosis was obtained from the gene expression omnibus (GEO) database. Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were applied to the endometriosis GSE7305 datasets. Cibersort and MCP-counter were used to explore the immune response gene sets, immune response pathway, and immune environment. Differentially expressed genes (DEGs) were identified and screened. Common biological pathways were being investigated using the kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis. Transcription factors were from The Human Transcription Factors. The least absolute shrinkage and selection operator (Lasso) model identified four differential expressions of transcription factors (AEBP1, HOXB6, KLF2, and RORB). Their diagnostic value was calculated by receiver operating characteristic (ROC) curve analysis and validated in the validation cohort (GSE11691, GSE23339). By constructing the interaction network of crucial transcription factors, weighted gene coexpression network analysis (WGCNA) was used to search for key module genes. Metascape was used for enrichment analysis of essential module genes and obtained HOXB6, KLF2. The HOXB6 and KLF2 were further verified as the only two intersection genes according to Support Vector Machine Recursive Feature Elimination (SVM-RFE) and random forest models. We constructed ceRNA (lncRNA-miRNA-mRNA) networks with four potential transcription factors. Finally, we performed molecular docking for goserelin and dienogest with four transcription factors (AEBP1, HOXB6, KLF2, and RORB) to screen potential drug targets.ResultsImmune and metabolic pathways were enriched in GSVA and GSEA. In single sample gene set enrichment analysis (ssGSEA), most immune infiltrating cells, immune response gene sets, and immune response pathways are differentially expressed between endometriosis and non-endometriosis. Twenty-seven transcription factors were screened from differentially expressed genes. Most of the twenty-seven transcription factors were correlated with immune infiltrating cells, immune response gene sets and immune response pathways. Furthermore, Adipocyte enhancer binding protein 1 (AEBP1), Homeobox B6 (HOXB6), Kruppel Like Factor 2 (KLF2) and RAR Related Orphan Receptor B (RORB) were selected out from twenty-seven transcription factors. ROC analysis showed that the four genes had a high diagnostic value for endometriosis. In addition, KLF2 and HOXB6 were found to play particularly important roles in multiple modules (String, WGCNA, SVM-RFE, random forest) on the gene interaction network. Using the ceRNA network, we found that NEAT1 may regulate the expressions of AEBP1, HOXB6 and RORB, while X Inactive Specific Transcript (XIST) may control the expressions of HOXB6, RORB and KLF2. Finally, we found that goserelin and dienogest may be potential drugs to regulate AEBP1, HOXB6, KLF2 and RORB through molecular docking.ConclusionsAEBP1, HOXB6, KLF2, and RORB may be potential biomarkers for endometriosis. Two of them, KLF2 and HOXB6, are critical molecules in the gene interaction network of endometriosis. Discovered by molecular docking, AEBP1, HOXB6, KLF2, and RORB are targets for goserelin and dienogest.

Dataset Information

Integrated analysis of motif activity and gene expression changes of transcription factors.

Publications

Integrated analysis of motif activity and gene expression changes of transcription factors.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets