Functional Enhancers As Master Regulators of Tissue-Specific Gene Regulation and Cancer Development

Genetically regulated eRNA expression predicts chromatin contact frequency and reveals genetic mechanisms at GWAS loci

The biological functions of extragenic enhancer RNAs and their impact on disease risk remain relatively underexplored. In this work, we develop in silico models of genetically regulated expression of enhancer RNAs across 49 cell and tissue types, characterizing their degree of genetic control. Leveraging the estimated genetically regulated expression for enhancer RNAs and canonical genes in a large-scale DNA biobank (N > 70,000) and high-resolution Hi-C contact data, we train a deep learning-based model of pairwise three-dimensional chromatin contact frequency for enhancer-enhancer and enhancer-gene pairs in cerebellum and whole blood. Notably, the use of genetically regulated expression of enhancer RNAs provides substantial tissue-specific predictive power, supporting a role for these transcripts in modulating spatial chromatin organization. We identify schizophrenia-associated enhancer RNAs independent of GWAS loci using enhancer RNA-based TWAS and determine the causal effects of these enhancer RNAs using Mendelian randomization. Using enhancer RNA-based TWAS, we generate a comprehensive resource of tissue-specific enhancer associations with complex traits in the UK Biobank. Finally, we show that a substantially greater proportion (63%) of GWAS associations colocalize with causal regulatory variation when enhancer RNAs are included.

Super-enhancer omics in stem cell

The hallmarks of stem cells, such as proliferation, self-renewal, development, differentiation, and regeneration, are critical to maintain stem cell identity which is sustained by genetic and epigenetic factors. Super-enhancers (SEs), which consist of clusters of active enhancers, play a central role in maintaining stemness hallmarks by specifically transcriptional model. The SE-navigated transcriptional complex, including SEs, non-coding RNAs, master transcriptional factors, Mediators and other co-activators, forms phase-separated condensates, which offers a toggle for directing diverse stem cell fate. With the burgeoning technologies of multiple-omics applied to examine different aspects of SE, we firstly raise the concept of "super-enhancer omics", inextricably linking to Pan-omics. In the review, we discuss the spatiotemporal organization and concepts of SEs, and describe links between SE-navigated transcriptional complex and stem cell features, such as stem cell identity, self-renewal, pluripotency, differentiation and development. We also elucidate the mechanism of stemness and oncogenic SEs modulating cancer stem cells via genomic and epigenetic alterations hijack in cancer stem cell. Additionally, we discuss the potential of targeting components of the SE complex using small molecule compounds, genome editing, and antisense oligonucleotides to treat SE-associated organ dysfunction and diseases, including cancer. This review also provides insights into the future of stem cell research through the paradigm of SEs.

Placental co-transcriptional activator Vestigial-like 1 (VGLL1) drives tumorigenesis via increasing transcription of proliferation and invasion genes

Introduction

Vestigial-like 1 (VGLL1) is a co-transcriptional activator that binds to TEA domain-containing transcription factors (TEADs). Its expression is upregulated in a variety of aggressive cancer types, including pancreatic and basal-like breast cancer, and increased transcription of VGLL1 is strongly correlated with poor prognosis and decreased overall patient survival. In normal tissues, VGLL1 is most highly expressed within placental trophoblast cells, which share the common attributes of rapid cellular proliferation and invasion with tumor cells. The impact of VGLL1 in cancer has not been fully elucidated and no VGLL1-targeted therapy currently exists.

Methods

The aim of this study was to evaluate the cellular function and downstream genomic targets of VGLL1 in placental, pancreatic, and breast cancer cells. Functional assays were employed to assess the role of VGLL1 in cellular invasion and proliferation, and ChIP-seq and RNAseq assays were performed to identify VGLL1 target genes and potential impact using pathway analysis.

Results

ChIP-seq analysis identified eight transcription factors with a VGLL1-binding motif that were common between all three cell types, including TEAD1-4, AP-1, and GATA6, and revealed ~3,000 shared genes with which VGLL1 interacts. Furthermore, increased VGLL1 expression led to an enhancement of cell invasion and proliferation, which was supported by RNAseq analysis showing transcriptional changes in several genes known to be involved in these processes.

Discussion

This work expands our mechanistic understanding of VGLL1 function in tumor cells and provides a strong rationale for developing VGLL1-targeted therapies for treating cancer patients.

Identification of drug responsive enhancers by predicting chromatin accessibility change from perturbed gene expression profiles

Individual may response to drug treatment differently due to their genetic variants located in enhancers. These variants can alter transcription factor's (TF) binding strength, affect enhancer's chromatin activity or interaction, and eventually change expression level of downstream gene. Here, we propose a computational framework, PERD, to Predict the Enhancers Responsive to Drug. A machine learning model was trained to predict the genome-wide chromatin accessibility from transcriptome data using the paired expression and chromatin accessibility data collected from ENCODE and ROADMAP. Then the model was applied to the perturbed gene expression data from Connectivity Map (CMAP) and Cancer Drug-induced gene expression Signature DataBase (CDS-DB) and identify drug responsive enhancers with significantly altered chromatin accessibility. Furthermore, the drug responsive enhancers were related to the pharmacogenomics genome-wide association studies (PGx GWAS). Stepping on the traditional drug-associated gene signatures, PERD holds the promise to enhance the causality of drug perturbation by providing candidate regulatory element of those drug associated genes.

HPV16-miRNAs exert oncogenic effects through enhancers in human cervical cancer

BACKGROUND

Cervical cancer is a human papillomavirus (HPV)-related disease. HPV type 16 (HPV16), which is the predominant cause of cervical cancer, can encode miRNAs (HPV16-miRNAs). However, the role of HPV16-miRNAs in the pathogenesis of cervical cancer remains unclear.

METHODS

Human cervical cancer cell lines SiHa (HPV16-positive) and C33A (HPV-negative), and cervical cancer tissues were collected to investigate the expression levels of two HPV16-miRNAs (HPV16-miR-H1 and HPV16-miR-H6). The overexpression and knockdown of HPV16-miR-H1 and HPV16-miR-H6 were performed using the lentiviral vector system and miRNA inhibitors, respectively. RNA-sequencing (RNA-seq) analysis and H3K27ac chromatin immunoprecipitation and sequencing (CHIP-seq) experiments were utilized to explore the roles of HPV16-miR-H1 and HPV16-miR-H6 facilitated by enhancers. CCK8, EdU, transwell, and wound healing assays were performed to verify the effects of HPV16-miR-H1 and HPV16-miR-H6 on cell proliferation and migration.

RESULTS

HPV16-miR-H1 and HPV16-miR-H6 were highly expressed in both SiHa cells and tissue samples from HPV16-positive cervical cancer patients. RNA-seq analysis showed that HPV16-miR-H1 and HPV16-miR-H6 induced the upregulation of numerous tumor progression-associated genes. H3K27ac CHIP-seq experiments further revealed that HPV16-miR-H1 and HPV16-miR-H6 modulated the expression of critical genes by regulating their enhancer activity. The functional study demonstrated that HPV16-miR-H1 and HPV16-miR-H6 increased the migratory capacity of SiHa cells.

CONCLUSIONS

Our data shed light on the role of HPV16-encoded miRNAs in cervical cancer, particularly emphasizing their involvement in the miRNA-enhancer-target gene system. This novel regulatory mechanism of HPV16-miRNAs provides new insights and approaches for the development of therapeutic strategies by targeting HPV16-positive cervical cancer.

Kidney-specific methylation patterns correlate with kidney function and are lost upon kidney disease progression

BACKGROUND

Chronological and biological age correlate with DNA methylation levels at specific sites in the genome. Linear combinations of multiple methylation sites, termed epigenetic clocks, can inform us the chronological age and predict multiple health-related outcomes. However, why some sites correlating with lifespan, healthspan, or specific medical conditions remain poorly understood. Kidney fibrosis is the common pathway for chronic kidney disease, which affects 10% of European and US populations.

RESULTS

Here we identify epigenetic clocks and methylation sites that correlate with kidney function. Moreover, we identify methylation sites that have a unique methylation signature in the kidney. Methylation levels in majority of these sites correlate with kidney state and function. When kidney function deteriorates, all of these sites regress toward the common methylation pattern observed in other tissues. Interestingly, while the majority of sites are less methylated in the kidney and become more methylated with loss of function, a fraction of the sites are highly methylated in the kidney and become less methylated when kidney function declines. These methylation sites are enriched for specific transcription-factor binding sites. In a large subset of sites, changes in methylation patterns are accompanied by changes in gene expression in kidneys of chronic kidney disease patients.

CONCLUSIONS

These results support the information theory of aging, and the hypothesis that the unique tissue identity, as captured by methylation patterns, is lost as tissue function declines. However, this information loss is not random, but guided toward a baseline that is dependent on the genomic loci.

SIGNIFICANCE STATEMENT

DNA methylation at specific sites accurately reflects chronological and biological age. We identify sites that have a unique methylation pattern in the kidney. Methylation levels in the majority of these sites correlate with kidney state and function. Moreover, when kidney function deteriorates, all of these sites regress toward the common methylation pattern observed in other tissues. Thus, the unique methylation signature of the kidney is degraded, and epigenetic information is lost, when kidney disease progresses. These methylation sites are enriched for specific and methylation-sensitive transcription-factor binding sites, and associated genes show disease-dependent changes in expression. These results support the information theory of aging, and the hypothesis that the unique tissue identity, as captured by methylation patterns, is lost as tissue function declines.

Quantitative analysis of the UDP-glucuronosyltransferase transcriptome in human tissues

UDP-glucuronosyltransferases (UGTs) are phase II drug metabolizing enzymes that play important roles in the detoxification of endogenous and exogenous substrates. The 22 human UGTs belong to four families (UGT1, UGT2, UGT3, and UGT8) and differ in their expression, substrate specificity, UDP-sugar preference, and physiological functions. Differential expression/activity of the UGTs contributes to interperson variability in drug responses and toxicity, hormone homeostasis, and disease/cancer risks. However, in normal tissues, the tissue-specific expression profiles and transcriptional regulation of the UGTs are still not fully understood. In this study, we comprehensively analyzed the transcriptome of 22 UGTs in 54 human tissues/regions using RNAseq data from GTEx. We then validated the findings in the liver and small intestine samples using real-time PCR. Our results showed large interindividual variability across tissues in the expression of each UGT and the overall composition of UGT pools, consisting of different UGTs and their splice isoforms. Our results also revealed coexpression of the UGTs, Cytochrome P450s, and many transcription factors in the liver, suggesting potential coregulation or functional coordination. Our results provide the groundwork for future studies to detail further the regulation of the expression and activity of the UGTs.

Dynamic modulation of genomic enhancer elements in the suprachiasmatic nucleus, the site of the mammalian circadian clock

The mammalian suprachiasmatic nucleus (SCN), located in the ventral hypothalamus, synchronizes and maintains daily cellular and physiological rhythms across the body, in accordance with environmental and visceral cues. Consequently, the systematic regulation of spatiotemporal gene transcription in the SCN is vital for daily timekeeping. So far, the regulatory elements assisting circadian gene transcription have only been studied in peripheral tissues, lacking the critical neuronal dimension intrinsic to the role of the SCN as central brain pacemaker. By using histone-ChIP-seq, we identified SCN-enriched gene regulatory elements that associated with temporal gene expression. Based on tissue-specific H3K27ac and H3K4me3 marks, we successfully produced the first-ever SCN gene-regulatory map. We found that a large majority of SCN enhancers not only show robust 24-h rhythmic modulation in H3K27ac occupancy, peaking at distinct times of day, but also possess canonical E-box (CACGTG) motifs potentially influencing downstream cycling gene expression. To establish enhancer-gene relationships in the SCN, we conducted directional RNA-seq at six distinct times across the day and night, and studied the association between dynamically changing histone acetylation and gene transcript levels. About 35% of the cycling H3K27ac sites were found adjacent to rhythmic gene transcripts, often preceding the rise in mRNA levels. We also noted that enhancers encompass noncoding, actively transcribing enhancer RNAs (eRNAs) in the SCN, which in turn oscillate, along with cyclic histone acetylation, and correlate with rhythmic gene transcription. Taken together, these findings shed light on genome-wide pretranscriptional regulation operative in the central clock that confers its precise and robust oscillation necessary to orchestrate daily timekeeping in mammals.

Cell-specific expression of the FAP gene is regulated by enhancer elements

Fibroblast activation protein (FAP) is an integral membrane serine protease that acts as both dipeptidyl peptidase and collagenase. In recent years, FAP has attracted considerable attention due to its specific upregulation in multiple types of tumor cell populations, including cancer cells in various cancer types, making FAP a potential target for therapy. However, relatively few papers pay attention to the mechanisms driving the cell-specific expression of the FAP gene. We found no correlation between the activities of the two FAP promoter variants (short and long) and the endogenous FAP mRNA expression level in several cell lines with different FAP expression levels. This suggested that other mechanisms may be responsible for specific transcriptional regulation of the FAP gene. We analyzed the distribution of known epigenetic and structural chromatin marks in FAP-positive and FAP-negative cell lines and identified two potential enhancer-like elements (E1 and E2) in the FAP gene locus. We confirmed the specific enrichment of H3K27ac in the putative enhancer regions in FAP-expressing cells. Both the elements exhibited enhancer activity independently of each other in the functional test by increasing the activity of the FAP promoter variants to a greater extent in FAP-expressing cell lines than in FAP-negative cell lines. The transcription factors AP-1, CEBPB, and STAT3 may be involved in FAP activation in the tumors. We hypothesized the existence of a positive feedback loop between FAP and STAT3, which may have implications for developing new approaches in cancer therapy.

E2F1 combined with LINC01004 super-enhancer to promote hepatocellular carcinoma cell proliferation and metastasis

INTRODUCTION

Super-enhancer-associated lncRNAs play important roles in the occurrence and development of malignant tumors, including hepatocellular carcinoma (HCC).

OBJECTIVES

The current work aimed to identify and characterize super-enhancer-associated lncRNAs in the pathogenesis of HCC.

METHODS

H3K27ac ChIP-seq data from HepG2 cell line and two HCC tissues were used to identify super-enhancer-associated lncRNAs in HCC. JQ-1 treatment and CRISPR-dCas9 system were performed to confirm super-enhancer activity. Quantitative real-time PCR (qPCR), ChIP-qPCR, and dual-luciferase reporter system assay demonstrated the regulation of E2F1 on super-enhancer. Functional loss experiment was used to identify the function of LINC01004.

RESULTS

In this study, we identified and characterized LINC01004, a novel super-enhancer-associated lncRNA, as a crucial oncogene in HCC. LINC01004 was upregulated in liver cancer tissues and was associated with poor patient prognosis. Moreover, LINC01004 promoted cell proliferation and metastasis of HCC. The binding of E2F1 to the super-enhancer could promote the transcription of LINC01004, while the inhibition of super-enhancer activity decreased LINC01004 expression.

CONCLUSION

This finding might provide mechanistic insights into the molecular mechanisms underlying hepatocarcinogenesis and the biological function of super-enhancer. LINC01004 can serve as a potential therapeutic target for HCC patient.

RANKL-responsive epigenetic mechanism reprograms macrophages into bone-resorbing osteoclasts

Monocyte/macrophage lineage cells are highly plastic and can differentiate into various cells under different environmental stimuli. Bone-resorbing osteoclasts are derived from the monocyte/macrophage lineage in response to receptor activator of NF-κB ligand (RANKL). However, the epigenetic signature contributing to the fate commitment of monocyte/macrophage lineage differentiation into human osteoclasts is largely unknown. In this study, we identified RANKL-responsive human osteoclast-specific superenhancers (SEs) and SE-associated enhancer RNAs (SE-eRNAs) by integrating data obtained from ChIP-seq, ATAC-seq, nuclear RNA-seq and PRO-seq analyses. RANKL induced the formation of 200 SEs, which are large clusters of enhancers, while suppressing 148 SEs in macrophages. RANKL-responsive SEs were strongly correlated with genes in the osteoclastogenic program and were selectively increased in human osteoclasts but marginally presented in osteoblasts, CD4+ T cells, and CD34+ cells. In addition to the major transcription factors identified in osteoclasts, we found that BATF binding motifs were highly enriched in RANKL-responsive SEs. The depletion of BATF1/3 inhibited RANKL-induced osteoclast differentiation. Furthermore, we found increased chromatin accessibility in SE regions, where RNA polymerase II was significantly recruited to induce the extragenic transcription of SE-eRNAs, in human osteoclasts. Knocking down SE-eRNAs in the vicinity of the NFATc1 gene diminished the expression of NFATc1, a major regulator of osteoclasts, and osteoclast differentiation. Inhibiting BET proteins suppressed the formation of some RANKL-responsive SEs and NFATc1-associated SEs, and the expression of SE-eRNA:NFATc1. Moreover, SE-eRNA:NFATc1 was highly expressed in the synovial macrophages of rheumatoid arthritis patients exhibiting high-osteoclastogenic potential. Our genome-wide analysis revealed RANKL-inducible SEs and SE-eRNAs as osteoclast-specific signatures, which may contribute to the development of osteoclast-specific therapeutic interventions.

Construction of a transposase accessible chromatin landscape reveals chromatin state of repeat elements and potential causal variant for complex traits in pigs

Background

A comprehensive landscape of chromatin states for multiple mammalian tissues is essential for elucidating the molecular mechanism underlying regulatory variants on complex traits. However, the genome-wide chromatin accessibility has been only reported in limited tissue types in pigs.

Results

Here we report a genome-wide landscape of chromatin accessibility of 20 tissues in two female pigs at ages of 6 months using ATAC-seq, and identified 557,273 merged peaks, which greatly expanded the pig regulatory element repository. We revealed tissue-specific regulatory elements which were associated with tissue-relevant biological functions. We identified both positive and negative significant correlations between the regulatory elements and gene transcripts, which showed distinct distributions in terms of their strength and distances from corresponding genes. We investigated the presence of transposable elements (TEs) in open chromatin regions across all tissues, these included identifications of porcine endogenous retroviruses (PERVs) exhibiting high accessibility in liver and homology of porcine specific virus sequences to universally accessible transposable elements. Furthermore, we prioritized a potential causal variant for polyunsaturated fatty acid in the muscle.

Conclusions

Our data provides a novel multi-tissues accessible chromatin landscape that serve as an important resource for interpreting regulatory sequences in tissue-specific and conserved biological functions, as well as regulatory variants of loci associated with complex traits in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00767-3.

SNP-to-gene linking strategies reveal contributions of enhancer-related and candidate master-regulator genes to autoimmune disease

We assess contributions to autoimmune disease of genes whose regulation is driven by enhancer regions (enhancer-related) and genes that regulate other genes in trans (candidate master-regulator). We link these genes to SNPs using several SNP-to-gene (S2G) strategies and apply heritability analyses to draw three conclusions about 11 autoimmune/blood-related diseases/traits. First, several characterizations of enhancer-related genes using functional genomics data are informative for autoimmune disease heritability after conditioning on a broad set of regulatory annotations. Second, candidate master-regulator genes defined using trans-eQTL in blood are also conditionally informative for autoimmune disease heritability. Third, integrating enhancer-related and master-regulator gene sets with protein-protein interaction (PPI) network information magnified their disease signal. The resulting PPI-enhancer gene score produced >2-fold stronger heritability signal and >2-fold stronger enrichment for drug targets, compared with the recently proposed enhancer domain score. In each case, functionally informed S2G strategies produced 4.1- to 13-fold stronger disease signals than conventional window-based strategies.

Mapping cis-regulatory elements in the midgestation mouse placenta

The placenta is a temporary organ that provides the developing fetus with nutrients, oxygen, and protection in utero. Defects in its development, which may be caused by misregulated gene expression, can lead to devastating outcomes for the mother and fetus. In mouse, placental defects during midgestation commonly lead to embryonic lethality. However, the regulatory mechanisms controlling expression of genes during this period have not been thoroughly investigated. Therefore, we generated and analyzed ChIP-seq data for multiple histone modifications known to mark cis-regulatory regions. We annotated active and poised promoters and enhancers, as well as regions generally associated with repressed gene expression. We found that poised promoters were associated with neuronal development genes, while active promoters were largely associated with housekeeping genes. Active and poised enhancers were associated with placental development genes, though only active enhancers were associated with genes that have placenta-specific expression. Motif analysis within active enhancers identified a large network of transcription factors, including those that have not been previously studied in the placenta and are candidates for future studies. The data generated and genomic regions annotated provide researchers with a foundation for future studies, aimed at understanding how specific genes in the midgestation mouse placenta are regulated.

Nuclear Dishevelled targets gene regulatory regions and promotes tumor growth

Dishevelled (DVL) critically regulates Wnt signaling and contributes to a wide spectrum of diseases and is important in normal and pathophysiological settings. However, how it mediates diverse cellular functions remains poorly understood. Recent discoveries have revealed that constitutive Wnt pathway activation contributes to breast cancer malignancy, but the mechanisms by which this occurs are unknown and very few studies have examined the nuclear role of DVL. Here, we have performed DVL3 ChIP-seq analyses and identify novel target genes bound by DVL3. We show that DVL3 depletion alters KMT2D binding to novel targets and changes their epigenetic marks and mRNA levels. We further demonstrate that DVL3 inhibition leads to decreased tumor growth in two different breast cancer models in vivo. Our data uncover new DVL3 functions through its regulation of multiple genes involved in developmental biology, antigen presentation, metabolism, chromatin remodeling, and tumorigenesis. Overall, our study provides unique insight into the function of nuclear DVL, which helps to define its role in mediating aberrant Wnt signaling.

MOSTWAS: Multi-Omic Strategies for Transcriptome-Wide Association Studies

Traditional predictive models for transcriptome-wide association studies (TWAS) consider only single nucleotide polymorphisms (SNPs) local to genes of interest and perform parameter shrinkage with a regularization process. These approaches ignore the effect of distal-SNPs or other molecular effects underlying the SNP-gene association. Here, we outline multi-omics strategies for transcriptome imputation from germline genetics to allow more powerful testing of gene-trait associations by prioritizing distal-SNPs to the gene of interest. In one extension, we identify mediating biomarkers (CpG sites, microRNAs, and transcription factors) highly associated with gene expression and train predictive models for these mediators using their local SNPs. Imputed values for mediators are then incorporated into the final predictive model of gene expression, along with local SNPs. In the second extension, we assess distal-eQTLs (SNPs associated with genes not in a local window around it) for their mediation effect through mediating biomarkers local to these distal-eSNPs. Distal-eSNPs with large indirect mediation effects are then included in the transcriptomic prediction model with the local SNPs around the gene of interest. Using simulations and real data from ROS/MAP brain tissue and TCGA breast tumors, we show considerable gains of percent variance explained (1-2% additive increase) of gene expression and TWAS power to detect gene-trait associations. This integrative approach to transcriptome-wide imputation and association studies aids in identifying the complex interactions underlying genetic regulation within a tissue and important risk genes for various traits and disorders.

ARID1A Mutations Promote P300-Dependent Endometrial Invasion through Super-Enhancer Hyperacetylation

Endometriosis affects 1 in 10 women and is characterized by the presence of abnormal endometrium at ectopic sites. ARID1A mutations are observed in deeply invasive forms of the disease, often correlating with malignancy. To identify epigenetic dependencies driving invasion, we use an unbiased approach to map chromatin state transitions accompanying ARID1A loss in the endometrium. We show that super-enhancers marked by high H3K27 acetylation are strongly associated with ARID1A binding. ARID1A loss leads to H3K27 hyperacetylation and increased chromatin accessibility and enhancer RNA transcription at super-enhancers, but not typical enhancers, indicating that ARID1A normally prevents super-enhancer hyperactivation. ARID1A co-localizes with P300 at super-enhancers, and genetic or pharmacological inhibition of P300 in ARID1A mutant endometrial epithelia suppresses invasion and induces anoikis through the rescue of super-enhancer hyperacetylation. Among hyperactivated super-enhancers, SERPINE1 (PAI-1) is identified as an essential target gene driving ARID1A mutant endometrial invasion. Broadly, our findings provide rationale for therapeutic strategies targeting super-enhancers in ARID1A mutant endometrium.

Aberrant activation of super enhancer and choline metabolism drive antiandrogen therapy resistance in prostate cancer

Next generation antiandrogens such as enzalutamide (Enz) are effective initially for the treatment of castration-resistant prostate cancer (CRPC). However, the disease often relapses and the underlying mechanisms remain elusive. By performing H3-lysine-27 acetylation (H3K27ac) ChIP-seq in Enz-resistant CRPC cells, we identified a group of super enhancers (SEs) that are abnormally activated in Enz-resistant CRPC cells and associated with enhanced transcription of a subset of tumor promoting genes such as CHPT1, which catalyzes phosphatidylcholine (PtdCho) synthesis and regulates choline metabolism. Increased CHPT1 conferred CRPC resistance to Enz in vitro and in mice. While androgen receptor (AR) primarily binds to a putative CHPT1 enhancer and mediates androgen-dependent expression of CHPT1 gene in Enz-sensitive prostate cancer cells, AR binds to a different enhancer within the CHPT1 SE locus and facilities androgen-independent expression of CHPT1 in Enz-resistant cells. We further identified a long-non coding RNA transcribed at CHPT1 enhancer (also known as enhancer RNA) that binds to the H3K27ac reader BRD4 and participates in regulating CHPT1 SE activity and CHPT1 gene expression. Our findings demonstrate that aberrantly activated SE upregulates CHPT1 expression and confers Enz resistance in CRPC, suggesting that SE-mediated expression of downstream effectors such as CHPT1 can be viable targets to overcome Enz resistance in PCa.

Polymorphic mobile element insertions contribute to gene expression and alternative splicing in human tissues

BACKGROUND

Mobile elements are a major source of structural variants in the human genome, and some mobile elements can regulate gene expression and transcript splicing. However, the impact of polymorphic mobile element insertions (pMEIs) on gene expression and splicing in diverse human tissues has not been thoroughly studied. The multi-tissue gene expression and whole genome sequencing data generated by the Genotype-Tissue Expression (GTEx) project provide a great opportunity to systematically evaluate the role of pMEIs in regulating gene expression in human tissues.

RESULTS

Using the GTEx whole genome sequencing data, we identify 20,545 high-quality pMEIs from 639 individuals. Coupling pMEI genotypes with gene expression profiles, we identify pMEI-associated expression quantitative trait loci (eQTLs) and splicing quantitative trait loci (sQTLs) in 48 tissues. Using joint analyses of pMEIs and other genomic variants, pMEIs are predicted to be the potential causal variant for 3522 eQTLs and 3717 sQTLs. The pMEI-associated eQTLs and sQTLs show a high level of tissue specificity, and these pMEIs are enriched in the proximity of affected genes and in regulatory elements. Using reporter assays, we confirm that several pMEIs associated with eQTLs and sQTLs can alter gene expression levels and isoform proportions, respectively.

CONCLUSION

Overall, our study shows that pMEIs are associated with thousands of gene expression and splicing variations, indicating that pMEIs could have a significant role in regulating tissue-specific gene expression and transcript splicing. Detailed mechanisms for the role of pMEIs in gene regulation in different tissues will be an important direction for future studies.

Emerging Roles of Estrogen-Regulated Enhancer and Long Non-Coding RNAs

Genome-wide RNA sequencing has shown that only a small fraction of the human genome is transcribed into protein-coding mRNAs. While once thought to be “junk” DNA, recent findings indicate that the rest of the genome encodes many types of non-coding RNA molecules with a myriad of functions still being determined. Among the non-coding RNAs, long non-coding RNAs (lncRNA) and enhancer RNAs (eRNA) are found to be most copious. While their exact biological functions and mechanisms of action are currently unknown, technologies such as next-generation RNA sequencing (RNA-seq) and global nuclear run-on sequencing (GRO-seq) have begun deciphering their expression patterns and biological significance. In addition to their identification, it has been shown that the expression of long non-coding RNAs and enhancer RNAs can vary due to spatial, temporal, developmental, or hormonal variations. In this review, we explore newly reported information on estrogen-regulated eRNAs and lncRNAs and their associated biological functions to help outline their markedly prominent roles in estrogen-dependent signaling.

Ontology-driven integrative analysis of omics data through Onassis

Public repositories of large-scale omics datasets represent a valuable resource for researchers. In fact, data re-analysis can either answer novel questions or provide critical data able to complement in-house experiments. However, despite the development of standards for the compilation of metadata, the identification and organization of samples still constitutes a major bottleneck hampering data reuse. We introduce Onassis, an R package within the Bioconductor environment providing key functionalities of Natural Language Processing (NLP) tools. Leveraging biomedical ontologies, Onassis greatly simplifies the association of samples from large-scale repositories to their representation in terms of ontology-based annotations. Moreover, through the use of semantic similarity measures, Onassis hierarchically organizes the datasets of interest, thus supporting the semantically aware analysis of the corresponding omics data. In conclusion, Onassis leverages NLP techniques, biomedical ontologies, and the R statistical framework, to identify, relate, and analyze datasets from public repositories. The tool was tested on various large-scale datasets, including compendia of gene expression, histone marks, and DNA methylation, illustrating how it can facilitate the integrative analysis of various omics data.

Enhancer RNAs in cancer: regulation, mechanisms and therapeutic potential

Enhancers are distal genomic elements critical for gene regulation and cell identify control during development and diseases. Many human cancers were found to associate with enhancer malfunction, due to genetic and epigenetic alterations, which in some cases directly drive tumour growth. Conventionally, enhancers are known to provide DNA binding motifs to recruit transcription factors (TFs) and to control target genes. However, recent progress found that most, if not all, active enhancers pervasively transcribe noncoding RNAs that are referred to as enhancer RNAs (eRNAs). Increasing evidence points to functional roles of at least a subset of eRNAs in gene regulation in both normal and cancer cells, adding new insights into the action mechanisms of enhancers. eRNA expression was observed to be widespread but also specific to tumour types and individual patients, serving as opportunities to exploit them as potential diagnosis markers or therapeutic targets. In this review, we discuss the brief history of eRNA research, their functional mechanisms and importance in cancer gene regulation, as well as their therapeutic and diagnostic values in cancer. We propose that further studies of eRNAs in cancer will offer a promising 'eRNA targeted therapy' for human cancer intervention.

Investigating the role of super-enhancer RNAs underlying embryonic stem cell differentiation

BACKGROUND

Super-enhancer RNAs (seRNAs) are a kind of noncoding RNA transcribed from super-enhancer regions. The regulation mechanism and functional role of seRNAs are still unclear. Although super-enhancers play a critical role in the core transcriptional regulatory circuity of embryonic stem cell (ESC) differentiation, whether seRNAs have similar properties should be further investigated.

RESULTS

We analyzed cap analysis gene expression sequencing (CAGE-seq) datasets collected during the differentiation of embryonic stem cells (ESCs) to cardiomyocytes to identify the seRNAs. A non-negative matrix factorization algorithm was applied to decompose the seRNA profiles and reveal two hidden stages during the ESC differentiation. We further identified 95 and 78 seRNAs associated with early- and late-stage ESC differentiation, respectively. We found that the binding sites of master regulators of ESC differentiation, including NANOG, FOXA2, and MYC, were significantly observed in the loci of the stage-specific seRNAs. Based on the investigation of genes coexpressed with seRNA, these stage-specific seRNAs might be involved in cardiac-related functions such as myofibril assembly and heart development and act in trans to regulate the co-expressed genes.

CONCLUSIONS

In this study, we used a computational approach to demonstrate the possible role of seRNAs during ESC differentiation.

Super-Enhancer-Associated Hub Genes In Chronic Myeloid Leukemia Identified Using Weighted Gene Co-Expression Network Analysis

Purpose

Super-enhancer (SE)-associated oncogenes extensively potentiate the uncontrolled proliferation capacity of cancer cells. In this study, we aimed to identify the SE-associated hub genes associated with the clinical characteristics of chronic myeloid leukemia (CML).

Methods

Eigengenes from CML clinical modules were determined using weighted gene co-expression network analysis (WGCNA). Overlapping genes between eigengenes and SE-associated genes were used to construct protein–protein interaction (PPI) networks and annotate for pathway enrichment analysis. Expression patterns of the top-ranked SE-associated hub genes were further determined in CML patients and healthy controls via real-time PCR. After treatment of K562 cells with the BRD4 inhibitor, JQ1, for 24 hrs, mRNA and protein levels of SE-associated hub genes were evaluated using real-time PCR and Western blotting, respectively. H3K27ac, H3K4me1 and BRD4 ChIP-seq signal peaks were used to predict and identify SEs visualized by the Integrative Genomics Viewer.

Results

The yellow module was significantly related to the status and pathological phase of CML. SE-associated hub candidate genes were mainly enriched in the cell cycle pathway. Based on the PPI networks of hub genes and the top rank of degree, five SE-associated genes were identified: specifically, BUB1, CENPO, KIF2C, ORC1, and RRM2. Elevated expression of these five genes was not only related to CML status and phase but also positively regulated by SE and suppressed by the BRD4 inhibitor, JQ1, in K562 cells. Strong signal peaks of H3K27ac, H3K4me1 and BRD4 ChIP-seq of the five genes were additionally observed close to the predicted SE regions.

Conclusion

This is the first study to characterize SE-associated genes linked to clinical characteristics of CML via weighted gene co-expression network analysis. Our results support a novel mechanism involving aberrant expression of hub SE-associated genes in CML patients and K562 cells, and these genes will be potential new therapeutic targets for human leukemia.

Retroelement-Linked H3K4me1 Histone Tags Uncover Regulatory Evolution Trends of Gene Enhancers and Feature Quickly Evolving Molecular Processes in Human Physiology

Background: Retroelements (REs) are mobile genetic elements comprising ~40% of human DNA. They can reshape expression patterns of nearby genes by providing various regulatory sequences. The proportion of regulatory sequences held by REs can serve a measure of regulatory evolution rate of the respective genes and molecular pathways. Methods: We calculated RE-linked enrichment scores for individual genes and molecular pathways based on ENCODE project epigenome data for enhancer-specific histone modification H3K4me1 in five human cell lines. We identified consensus groups of molecular processes that are enriched and deficient in RE-linked H3K4me1 regulation. Results: We calculated H3K4me1 RE-linked enrichment scores for 24,070 human genes and 3095 molecular pathways. We ranked genes and pathways and identified those statistically significantly enriched and deficient in H3K4me1 RE-linked regulation. Conclusion: Non-coding RNA genes were statistically significantly enriched by RE-linked H3K4me1 regulatory modules, thus suggesting their high regulatory evolution rate. The processes of gene silencing by small RNAs, DNA metabolism/chromatin structure, sensory perception/neurotransmission and lipids metabolism showed signs of the fastest regulatory evolution, while the slowest processes were connected with immunity, protein ubiquitination/degradation, cell adhesion, migration and interaction, metals metabolism/ion transport, cell death, intracellular signaling pathways.

The therapeutic and diagnostic potential of regulatory noncoding RNAs in medulloblastoma

Medulloblastoma, a central nervous system tumor that predominantly affects children, always requires aggressive therapy. Nevertheless, it frequently recurs as resistant disease and is associated with high morbidity and mortality. While recent efforts to subclassify medulloblastoma based on molecular features have advanced our basic understanding of medulloblastoma pathogenesis, optimal targets to increase therapeutic efficacy and reduce side effects remain largely undefined. Noncoding RNAs (ncRNAs) with known regulatory roles, particularly long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), are now known to participate in medulloblastoma biology, although their functional significance remains obscure in many cases. Here we review the literature on regulatory ncRNAs in medulloblastoma. In providing a comprehensive overview of ncRNA studies, we highlight how different lncRNAs and miRNAs have oncogenic or tumor suppressive roles in medulloblastoma. These ncRNAs possess subgroup specificity that can be exploited to personalize therapy by acting as theranostic targets. Several of the already identified ncRNAs appear specific to medulloblastoma stem cells, the most difficult-to-treat component of the tumor that drives metastasis and acquired resistance, thereby providing opportunities for therapy in relapsing, disseminating, and therapy-resistant disease. Delivering ncRNAs to tumors remains challenging, but this limitation is gradually being overcome through the use of advanced technologies such as nanotechnology and rational biomaterial design.

A distal super enhancer mediates estrogen-dependent mouse uterine-specific gene transcription of Igf1 (insulin-like growth factor 1)

Insulin-like growth factor 1 (IGF1) is primarily synthesized in and secreted from the liver; however, estrogen (E2), through E2 receptor α (ERα), increases uterine Igf1 mRNA levels. Previous ChIP-seq analyses of the murine uterus have revealed a potential enhancer region distal from the Igf1 transcription start site (TSS) with multiple E2-dependent ERα-binding regions. Here, we show E2-dependent super enhancer-associated characteristics and suggest contact between the distal enhancer and the Igf1 TSS. We hypothesized that this distal super-enhancer region controls E2-responsive induction of uterine Igf1 transcripts. We deleted 430 bp, encompassing one of the ERα-binding sites, thereby disrupting interactions of the enhancer with gene-regulatory factors. As a result, E2-mediated induction of mouse uterine Igf1 mRNA is completely eliminated, whereas hepatic Igf1 expression remains unaffected. This highlights the central role of a distal enhancer in the assembly of the factors necessary for E2-dependent interaction with the Igf1 TSS and induction of uterus-specific Igf1 transcription. Of note, loss of the enhancer did not affect fertility or uterine growth responses. Deletion of uterine Igf1 in a PgrCre;Igf1f/f model decreased female fertility but did not impact the E2-induced uterine growth response. Moreover, E2-dependent activation of uterine IGF1 signaling was not impaired by disrupting the distal enhancer or by deleting the coding transcript. This indicated a role for systemic IGF1, suggested that other growth mediators drive uterine response to E2, and suggested that uterine-derived IGF1 is essential for reproductive success. Our findings elucidate the role of a super enhancer in Igf1 regulation and uterine growth.

Dysregulated Transcriptional Control in Prostate Cancer

Recent advances in whole-genome and transcriptome sequencing of prostate cancer at different stages indicate that a large number of mutations found in tumors are present in non-protein coding regions of the genome and lead to dysregulated gene expression. Single nucleotide variations and small mutations affecting the recruitment of transcription factor complexes to DNA regulatory elements are observed in an increasing number of cases. Genomic rearrangements may position coding regions under the novel control of regulatory elements, as exemplified by the TMPRSS2-ERG fusion and the amplified enhancer identified upstream of the androgen receptor (AR) gene. Super-enhancers are increasingly found to play important roles in aberrant oncogenic transcription. Several players involved in these processes are currently being evaluated as drug targets and may represent new vulnerabilities that can be exploited for prostate cancer treatment. They include factors involved in enhancer and super-enhancer function such as bromodomain proteins and cyclin-dependent kinases. In addition, non-coding RNAs with an important gene regulatory role are being explored. The rapid progress made in understanding the influence of the non-coding part of the genome and of transcription dysregulation in prostate cancer could pave the way for the identification of novel treatment paradigms for the benefit of patients.

The anti-cancer drugs curaxins target spatial genome organization

Recently we characterized a class of anti-cancer agents (curaxins) that disturbs DNA/histone interactions within nucleosomes. Here, using a combination of genomic and in vitro approaches, we demonstrate that curaxins strongly affect spatial genome organization and compromise enhancer-promoter communication, which is necessary for the expression of several oncogenes, including MYC. We further show that curaxins selectively inhibit enhancer-regulated transcription of chromatinized templates in cell-free conditions. Genomic studies also suggest that curaxins induce partial depletion of CTCF from its binding sites, which contributes to the observed changes in genome topology. Thus, curaxins can be classified as epigenetic drugs that target the 3D genome organization.

Curaxins are a recently discovered class of anti-cancer agents that disturbs DNA/histone interactions within. Here the authors provide evidence that curaxins affect the spatial genome organization and compromise enhancer-promoter communication necessary for expression of several oncogenes, including MYC.

MicroRNA-196a is regulated by ER and is a prognostic biomarker in ER+ breast cancer

Background

MicroRNAs are potent post-transcriptional regulators involved in all hallmarks of cancer. Mir-196a is transcribed from two loci and has been implicated in a wide range of developmental and pathogenic processes, with targets including Hox, Fox, Cdk inhibitors and annexins. Genetic variants and altered expression of MIR196A are associated with risk and progression of multiple cancers including breast cancer, however little is known about the regulation of the genes encoding this miRNA, nor the impact of variants therein.

Methods

Genomic data and chromatin interaction analysis were used to discover functional promoter and enhancer elements for MIR196A. Expression data were used to associate MIR196A with mechanisms of resistance, breast cancer subtypes and prognosis.

Results

Here we demonstrate that MIR196A displays complex and dynamic expression patterns, in part controlled by long-range transcriptional regulation between promoter and enhancer elements bound by ERα. Expression of this miRNA is significantly increased in drug-resistant models of hormone-receptor positive disease. The expression of MIR196A also proves to be a robust prognostic factor for patients with advanced and post-menopausal ER+ disease.

Conclusion

This work sheds light on the normal and abnormal regulation of MIR196A and provides a novel stratification method for therapeutically resistant breast cancer.

BC200 RNA: An Emerging Therapeutic Target and Diagnostic Marker for Human Cancer

One of the most interesting findings from genome-wide expression analysis is that a considerable amount of noncoding RNA (ncRNA) is present in the cell. Recent studies have identified diverse biological functions of ncRNAs, which are expressed in a much wider array of forms than proteins. Certain ncRNAs associated with diseases, in particular, have attracted research attention as novel therapeutic targets and diagnostic markers. BC200 RNA, a 200-nucleotide ncRNA originally identified as a neuron-specific transcript, is abnormally over-expressed in several types of cancer tissue. A number of recent studies have suggested mechanisms by which abnormal expression of BC200 RNA contributes to the development of cancer. In this article, we first provide a brief review of a recent progress in identifying functions of BC200 RNA in cancer cells, and then offer examples of other ncRNAs as new therapeutic targets and diagnostic markers for human cancer. Finally, we discuss future directions of studies on BC200 RNA for new cancer treatments.

Bromodomain protein inhibition: a novel therapeutic strategy in rheumatic diseases

The reading of acetylation marks on histones by bromodomain (BRD) proteins is a key event in transcriptional activation. Small molecule inhibitors targeting bromodomain and extra-terminal (BET) proteins compete for binding to acetylated histones. They have strong anti-inflammatory properties and exhibit encouraging effects in different cell types in vitro and in animal models resembling rheumatic diseases in vivo. Furthermore, recent studies that focus on BRD proteins beyond BET family members are discussed.

Enhancer and super-enhancer: Positive regulators in gene transcription

Enhancer is a positive regulator for spatiotemporal development in eukaryotes. As a cluster, super-enhancer is closely related to cell identity- and fate-determined processes. Both of them function tightly depending on their targeted transcription factors, cofactors, and genes through distal genomic interactions. They have been recognized as critical components and played positive roles in transcriptional regulatory network or factory. Recent advances of next-generation sequencing have dramatically expanded our ability and knowledge to interrogate the molecular mechanism of enhancer and super-enhancer for transcription. Here, we review the history, importance, advances and challenges on enhancer and super-enhancer field. This will benefit our understanding of their function mechanism for transcription underlying precise gene expression.

Targeting Transcription Factors for Cancer Treatment

Transcription factors are involved in a large number of human diseases such as cancers for which they account for about 20% of all oncogenes identified so far. For long time, with the exception of ligand-inducible nuclear receptors, transcription factors were considered as “undruggable” targets. Advances knowledge of these transcription factors, in terms of structure, function (expression, degradation, interaction with co-factors and other proteins) and the dynamics of their mode of binding to DNA has changed this postulate and paved the way for new therapies targeted against transcription factors. Here, we discuss various ways to target transcription factors in cancer models: by modulating their expression or degradation, by blocking protein/protein interactions, by targeting the transcription factor itself to prevent its DNA binding either through a binding pocket or at the DNA-interacting site, some of these inhibitors being currently used or evaluated for cancer treatment. Such different targeting of transcription factors by small molecules is facilitated by modern chemistry developing a wide variety of original molecules designed to specifically abort transcription factor and by an increased knowledge of their pathological implication through the use of new technologies in order to make it possible to improve therapeutic control of transcription factor oncogenic functions.

Targeting Super-Enhancers for Disease Treatment and Diagnosis

The transcriptional regulation of genes determines the fate of animal cell differentiation and subsequent organ development. With the recent progress in genome-wide technologies, the genomic landscapes of enhancers have been broadly explored in mammalian genomes, which led to the discovery of novel specific subsets of enhancers, termed super-enhancers. Super-enhancers are large clusters of enhancers covering the long region of regulatory DNA and are densely occupied by transcription factors, active histone marks, and co-activators. Accumulating evidence points to the critical role that super-enhancers play in cell type-specific development and differentiation, as well as in the development of various diseases. Here, I provide a comprehensive description of the optimal approach for identifying functional units of super-enhancers and their unique chromatin features in normal development and in diseases, including cancers. I also review the recent updated knowledge on novel approaches of targeting super-enhancers for the treatment of specific diseases, such as small-molecule inhibitors and potential gene therapy. This review will provide perspectives on using super-enhancers as biomarkers to develop novel disease diagnostic tools and establish new directions in clinical therapeutic strategies.

Dissecting super-enhancer hierarchy based on chromatin interactions

Recent studies have highlighted super-enhancers (SEs) as important regulatory elements for gene expression, but their intrinsic properties remain incompletely characterized. Through an integrative analysis of Hi-C and ChIP-seq data, here we find that a significant fraction of SEs are hierarchically organized, containing both hub and non-hub enhancers. Hub enhancers share similar histone marks with non-hub enhancers, but are distinctly associated with cohesin and CTCF binding sites and disease-associated genetic variants. Genetic ablation of hub enhancers results in profound defects in gene activation and local chromatin landscape. As such, hub enhancers are the major constituents responsible for SE functional and structural organization.

Alteration of Epigenetic Regulation by Long Noncoding RNAs in Cancer

Long noncoding RNAs (lncRNAs) are important regulators of the epigenetic status of the human genome. Besides their participation to normal physiology, lncRNA expression and function have been already associated to many diseases, including cancer. By interacting with epigenetic regulators and by controlling chromatin topology, their misregulation may result in an aberrant regulation of gene expression that may contribute to tumorigenesis. Here, we review the functional role and mechanisms of action of lncRNAs implicated in the aberrant epigenetic regulation that has characterized cancer development and progression.