A hidden Markov model to identify combinatorial epigenetic regulation patterns for estrogen receptor α target genes

The crucial role of epigenetic regulation in breast cancer anti-estrogen resistance: Current findings and future perspectives

Breast cancer (BC) cell de-sensitization to Tamoxifen (TAM) or other selective estrogen receptor (ER) modulators (SERM) is a complex process associated with BC heterogeneity and the transformation of ER signalling. The most influential resistance-related mechanisms include modifications in ER expression and gene regulation patterns. During TAM/SERM treatment, epigenetic mechanisms can effectively silence ER expression and facilitate the development of endocrine resistance. ER status is efficiently regulated by specific epigenetic tools including hypermethylation of CpG islands within ER promoters, increased histone deacetylase activity in the ER promoter, and/or translational repression by miRNAs. Over-methylation of the ER α gene (ESR1) promoter by DNA methyltransferases was associated with poor prognosis and indicated the development of resistance. Moreover, BC progression and spreading were marked by transformed chromatin remodelling, post-translational histone modifications, and expression of specific miRNAs and/or long non-coding RNAs. Therefore, targeted inhibition of histone acetyltransferases (e.g. MYST3), deacetylases (e.g. HDAC1), and/or demethylases (e.g. lysine-specific demethylase LSD1) was shown to recover and increase BC sensitivity to anti-estrogens. Indicated as a powerful molecular instrument, the administration of epigenetic drugs can regain ER expression along with the activation of tumour suppressor genes, which can in turn prevent selection of resistant cells and cancer stem cell survival. This review examines recent advances in the epigenetic regulation of endocrine drug resistance and evaluates novel anti-resistance strategies. Underlying molecular mechanisms of epigenetic regulation will be discussed, emphasising the utilization of epigenetic enzymes and their inhibitors to re-program irresponsive BCs.

Three-dimensional analysis reveals altered chromatin interaction by enhancer inhibitors harbors TCF7L2-regulated cancer gene signature

Distal regulatory elements influence the activity of gene promoters through chromatin looping. Chromosome conformation capture (3C) methods permit identification of chromatin contacts across different regions of the genome. However, due to limitations in the resolution of these methods, the detection of functional chromatin interactions remains a challenge. In the current study, we employ an integrated approach to define and characterize the functional chromatin contacts of human pancreatic cancer cells. We applied tethered chromatin capture to define classes of chromatin domains on a genome-wide scale. We identified three types of structural domains (topologically associated, boundary, and gap) and investigated the functional relationships of these domains with respect to chromatin state and gene expression. We uncovered six distinct sub-domains associated with epigenetic states. Interestingly, specific epigenetically active domains are sensitive to treatment with histone acetyltransferase (HAT) inhibitors and decrease in H3K27 acetylation levels. To examine whether the subdomains that change upon drug treatment are functionally linked to transcription factor regulation, we compared TCF7L2 chromatin binding and gene regulation to HAT inhibition. We identified a subset of coding RNA genes that together can stratify pancreatic cancer patients into distinct survival groups. Overall, this study describes a process to evaluate the functional features of chromosome architecture and reveals the impact of epigenetic inhibitors on chromosome architecture and identifies genes that may provide insight into disease outcome.

Transcription factor-associated combinatorial epigenetic pattern reveals higher transcriptional activity of TCF7L2-regulated intragenic enhancers

Background

Recent studies have suggested that combinations of multiple epigenetic modifications are essential for controlling gene expression. Despite numerous computational approaches have been developed to decipher the combinatorial epigenetic patterns or “epigenetic code”, none of them has explicitly addressed the relationship between a specific transcription factor (TF) and the patterns.

Methods

Here, we developed a novel computational method, T-cep, for annotating chromatin states associated with a specific TF. T-cep is composed of three key consecutive modules: (i) Data preprocessing, (ii) HMM training, and (iii) Potential TF-states calling.

Results

We evaluated T-cep on a TCF7L2-omics data. Unexpectedly, our method has uncovered a novel set of TCF7L2-regulated intragenic enhancers missed by other software tools, where the associated genes exert the highest gene expression. We further used siRNA knockdown, Co-transfection, RT-qPCR and Luciferase Reporter Assay not only to validate the accuracy and efficiency of prediction by T-cep, but also to confirm the functionality of TCF7L2-regulated enhancers in both MCF7 and PANC1 cells respectively.

Conclusions

Our study for the first time at a genome-wide scale reveals the enhanced transcriptional activity of cell-type-specific TCF7L2 intragenic enhancers in regulating gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3764-9) contains supplementary material, which is available to authorized users.

A computational model for predicting fusion peptide of retroviruses

As a pivotal domain within envelope protein, fusion peptide (FP) plays a crucial role in pathogenicity and therapeutic intervention. Taken into account the limited FP annotations in NCBI database and absence of FP prediction software, it is urgent and desirable to develop a bioinformatics tool to predict new putative FPs (np-FPs) in retroviruses. In this work, a sequence-based FP model was proposed by combining Hidden Markov Method with similarity comparison. The classification accuracies are 91.97% and 92.31% corresponding to 10-fold and leave-one-out cross-validation. After scanning sequences without FP annotations, this model discovered 53,946 np-FPs. The statistical results on FPs or np-FPs reveal that FP is a conserved and hydrophobic domain. The FP software programmed for windows environment is available at https://sourceforge.net/projects/fptool/files/?source=navbar.

Nuclear-encoded mitochondrial MTO1 and MRPL41 are regulated in an opposite epigenetic mode based on estrogen receptor status in breast cancer

BACKGROUND

MTO1 and MRPL41 are nuclear-encoded mitochondrial genes encoding a mitochondrial tRNA-modifying enzyme and a mitochondrial ribosomal protein, respectively. Although both genes have been known to have potential roles in cancer, little is known about their molecular regulatory mechanism, particularly from an epigenetic approach. In this study, we aimed to address their epigenetic regulation through the estrogen receptor (ER) in breast cancer.

METHODS

Digital differential display (DDD) was conducted to identify mammary gland-specific gene candidates including MTO1 and MRPL41. Promoter CpG methylation and expression in breast cancer cell lines and tissues were examined by methylation-specific PCR and real time RT-PCR. Effect of estradiol (E2), tamoxifen, and trichostatin A (TSA) on gene expression was examined in ER + and ER- breast cancer cell lines. Chromatin immunoprecipitation and luciferase reporter assay were performed to identify binding and influencing of the ER to the promoters.

RESULTS

Examination of both cancer tissues and cell lines revealed that the two genes showed an opposite expression pattern according to ER status; higher expression of MTO1 and MRPL41 in ER- and ER+ cancer types, respectively, and their expression levels were inversely correlated with promoter methylation. Tamoxifen, E2, and TSA upregulated MTO1 expression only in ER+ cells with no significant changes in ER- cells. However, these chemicals upregulated MRPL41 expression only in ER- cells without significant changes in ER+ cells, except for tamoxifen that induced downregulation. Chromatin immunoprecipitation and luciferase reporter assay identified binding and influencing of the ER to the promoters and the binding profiles were differentially regulated in ER+ and ER- cells.

CONCLUSIONS

These results indicate that different epigenetic status including promoter methylation and different responses through the ER are involved in the differential expression of MTO1 and MRPL41 in breast cancer.

LOcating non-unique matched tags (LONUT) to improve the detection of the enriched regions for ChIP-seq data

One big limitation of computational tools for analyzing ChIP-seq data is that most of them ignore non-unique tags (NUTs) that match the human genome even though NUTs comprise up to 60% of all raw tags in ChIP-seq data. Effectively utilizing these NUTs would increase the sequencing depth and allow a more accurate detection of enriched binding sites, which in turn could lead to more precise and significant biological interpretations. In this study, we have developed a computational tool, LOcating Non-Unique matched Tags (LONUT), to improve the detection of enriched regions from ChIP-seq data. Our LONUT algorithm applies a linear and polynomial regression model to establish an empirical score (ES) formula by considering two influential factors, the distance of NUTs to peaks identified using uniquely matched tags (UMTs) and the enrichment score for those peaks resulting in each NUT being assigned to a unique location on the reference genome. The newly located tags from the set of NUTs are combined with the original UMTs to produce a final set of combined matched tags (CMTs). LONUT was tested on many different datasets representing three different characteristics of biological data types. The detected sites were validated using de novo motif discovery and ChIP-PCR. We demonstrate the specificity and accuracy of LONUT and show that our program not only improves the detection of binding sites for ChIP-seq, but also identifies additional binding sites.

Comparative annotation of functional regions in the human genome using epigenomic data

Epigenetic regulation is dynamic and cell-type dependent. The recently available epigenomic data in multiple cell types provide an unprecedented opportunity for a comparative study of epigenetic landscape. We developed a machine-learning method called ChroModule to annotate the epigenetic states in eight ENCyclopedia Of DNA Elements cell types. The trained model successfully captured the characteristic histone-modification patterns associated with regulatory elements, such as promoters and enhancers, and showed superior performance on identifying enhancers compared with the state-of-art methods. In addition, given the fixed number of epigenetic states in the model, ChroModule allows straightforward illustration of epigenetic variability in multiple cell types. Using this feature, we found that invariable and variable epigenetic states across cell types correspond to housekeeping functions and stimulus response, respectively. Especially, we observed that enhancers, but not the other regulatory elements, dictate cell specificity, as similar cell types share common enhancers, and cell-type-specific enhancers are often bound by transcription factors playing critical roles in that cell type. More interestingly, we found some genomic regions are dormant in cell type but primed to become active in other cell types. These observations highlight the usefulness of ChroModule in comparative analysis and interpretation of multiple epigenomes.