Abstract:BACKGROUND:Although endocrine therapy impedes estrogen-ER signaling pathway and thus reduces breast cancer mortality, patients remain at continued risk of relapse after tamoxifen or other endocrine therapies. Understanding the mechanisms of endocrine resistance, particularly the role of transcriptional regulation is very important and necessary.METHODS:We propose a two-step workflow based on linear model to investigate the significant differences between MCF7 and OHT cells stimulated by 17β-estradiol (E2) respect to regulatory transcription factors (TFs) and their interactions. We additionally compared predicted regulatory TFs based on RNA polymerase II (PolII) binding quantity data and gene expression data, which were taken from MCF7/MCF7+E2 and OHT/OHT+E2 cell lines following the same analysis workflow. Enrichment analysis concerning diseases and cell functions and regulatory pattern analysis of different motifs of the same TF also were performed.RESULTS:The results showed PolII data could provide more information and predict more recognizably important regulatory TFs. Large differences in TF regulatory mode were found between two cell lines. Through verified through GO annotation, enrichment analysis and related literature regarding these TFs, we found some regulatory TFs such as AP-1, C/EBP, FoxA1, GATA1, Oct-1 and NF-κB, maintained OHT cells through molecular interactions or signaling pathways that were different from the surviving MCF7 cells. From TF regulatory interaction network, we identified E2F, E2F-1 and AP-2 as hub-TFs in MCF7 cells; whereas, in addition to E2F and E2F-1, we identified C/EBP and Oct-1 as hub-TFs in OHT cells. Notably, we found the regulatory patterns of different motifs of the same TF were very different from one another sometimes.CONCLUSIONS:We inferred some regulatory TFs, such as AP-1 and NF-κB, cooperated with ER through both genomic action and non-genomic action. The TFs that were involved in both protein-protein interactions and signaling pathways could be one of the key resistant mechanisms of endocrine therapy and thus also could be new treatment targets for endocrine resistance. Our flexible workflow could be integrated into an existing analytical framework and guide biologists to further determine underlying mechanisms in human diseases.

A Poisson Mixture Model to Identify Changes in RNA Polymerase II Binding Quantity Using High-Throughput Sequencing Technology

Inference of RNA Polymerase II Transcription Dynamics from Chromatin Immunoprecipitation Time Course Data

Bioinformatic analysis of changes in RNA polymerase II transcription stimulated by estradiol in MCF7 cells.

A signal processing approach for enriched region detection in RNA polymerase II ChIP-seq data

Modeling Expression Plasticity of Genes That Differentiate Drug-Sensitive from Drug-Resistant Cells to Chemotherapeutic Treatment

Differential Expression Analysis for RNAseq Using Poisson Mixed Models.

A Bi-Poisson Model for Clustering Gene Expression Profiles by RNA-seq.

Rna Polymerase Ii Binding Patterns Reveal Genomic Regions Involved in Microrna Gene Regulation

Differential expression analysis for paired RNA-seq data

Using Poisson mixed-effects model to quantify transcript-level gene expression in RNA-Seq.

Identification of potential models for predicting progestin insensitivity in patients with endometrial atypical hyperplasia and endometrioid endometrial cancer based on ATAC-Seq and RNA-Seq integrated analysis

Multi-scale Poisson process approaches for differential expression analysis of high-throughput sequencing data

EBSeq-HMM: a Bayesian Approach for Identifying Gene-Expression Changes in Ordered RNA-seq Experiments

Differential RNA Methylation Analysis for MeRIP-seq Data under General Experimental Design

Transcriptional Regulation Prediction of Antiestrogen Resistance in Breast Cancer Based on RNA Polymerase II Binding Data.

Construction of RNA polymerase II- Chromatin immunoprecipitation- primer extension-mass spectrum and its application in function research of single nucleotide polymorphism

Detecting m6A methylation regions from Methylated RNA Immunoprecipitation Sequencing

COMPUTATIONAL ANALYSIS AND MODELING OF GENOME-SCALE AVIDITY DISTRIBUTION OF TRANSCRIPTION FACTOR BINDING SITES IN CHIP-PET EXPERIMENTS

Spatially Enhanced Differential RNA Methylation Analysis from Affinity-Based Sequencing Data with Hidden Markov Model.

Identifying down and up-regulated chromosome regions using RNA-Seq data

Statistical Models for Detecting Differential Chromatin Interactions Mediated by a Protein.