Core transcriptional regulatory circuitries in cancer

BIT: Bayesian Identification of Transcriptional Regulators

BIT is a novel Bayesian hierarchical model capable of predicting transcriptional regulators (TRs) from the input of user-provided epigenomic regions. TRs are critical molecules in transcriptional regulation. Many diseases and cancers are linked to the dysfunction of TRs. Knowing TRs in certain biological process can help find new biomarkers or therapeutic targets. Thus, BIT formulates a novel Bayesian hierarchical model with the Pólya-gamma data augmentation strategy. Based on collected ChIP-seq datasets, BIT can identify TRs responsible for the genome-wide binding pattern within the user-provided epigenomic regions. BIT has been validated by using a simulation study and three applications.

Ligandability assessment of the C-terminal Rel-homology domain of NFAT1

Transcription factors are generally considered challenging, if not "undruggable", targets but they promise new therapeutic options due to their fundamental involvement in many diseases. In this study, we aim to assess the ligandability of the C-terminal Rel-homology domain of nuclear factor of activated T cells 1 (NFAT1), a TF implicated in T-cell regulation. Using a combination of experimental and computational approaches, we demonstrate that small molecule fragments can indeed bind to this protein domain. The newly identified binder is the first small molecule binder to NFAT1 validated with biophysical methods and an elucidated binding mode by X-ray crystallography. The reported eutomer/distomer pair provides a strong basis for potential exploration of higher potency binders on the path toward degrader or glue modalities.

Targeted Degradation of CDK9 Potently Disrupts the MYC Transcriptional Network

Cyclin-dependent kinase 9 (CDK9) coordinates signaling events that regulate RNA polymerase II (Pol II) pause-release states. It is an important co-factor for transcription factors, such as MYC, that drive aberrant cell proliferation when their expression is deregulated. CDK9 modulation offers an approach for attenuating dysregulation in such transcriptional programs. As a result, numerous drug development campaigns to inhibit CDK9 kinase activity have been pursued. More recently, targeted degradation has emerged as an attractive approach. However, comprehensive evaluation of degradation versus inhibition is still critically needed to assess the biological contexts in which degradation might offer superior therapeutic benefits. We validated that CDK9 inhibition triggers a compensatory mechanism that dampens its effect on MYC expression and found that this feedback mechanism was absent when the kinase is degraded. Importantly, CDK9 degradation is more effective than its inhibition for disrupting MYC transcriptional regulatory circuitry likely through the abrogation of both enzymatic and scaffolding functions of CDK9.

Highlights

- KI-CDK9d-32 is a highly potent and selective CDK9 degrader. - KI-CDK9d-32 leads to rapid downregulation of MYC protein and mRNA transcripts levels. - KI-CDK9d-32 represses canonical MYC pathways and leads to a destabilization of nucleolar homeostasis. - Multidrug resistance ABCB1 gene emerged as the strongest resistance marker for the CDK9 PROTAC degrader.

Transcription factors in chimeric antigen receptor T-cell development

Chimeric antigen receptor (CAR) T-cell therapy is a new and innovative approach to treating cancers that has shown promising results in the treatment of lymphoma. However, it has been found to be less effective in the treatment of solid tumors. To overcome the limitation, researchers have explored the use of combined CAR-T therapy with other complementary regimens that target specific genes or biomarkers, which would enhance the synergistic therapeutic effects. Transcription factors (TFs) have been identified as potential markers that can regulate gene expression in CAR-T cells to enhance their cytotoxicity and safety. TFs are known to bind DNA specifically and recruit cofactor proteins to regulate the expression of target genes. By targeting TFs, it is possible to improve the anti-tumor response of CAR-T cells by altering their phenotype and transcriptional map, thereby increasing their effector function, such as reducing the exhaustion, enhancing the survival, and cytotoxicity of CAR-T cells. This review summarizes the application of transcription factors in CART therapy to enhance the synergistic therapeutic effect of CAR-T cells in the treatment of solid tumors and improve their anti-tumor responses.

CDK5RAP3 is a novel super-enhancer-driven gene activated by master TFs and regulates ER-Phagy in neuroblastoma

Super enhancers (SEs) are genomic regions comprising multiple closely spaced enhancers, typically occupied by a high density of cell-type-specific master transcription factors (TFs) and frequently enriched in key oncogenes in various tumors, including neuroblastoma (NB), one of the most prevalent malignant solid tumors in children originating from the neural crest. Cyclin-dependent kinase 5 regulatory subunit-associated protein 3 (CDK5RAP3) is a newly identified super-enhancer-driven gene regulated by master TFs in NB; however, its function in NB remains unclear. Through an integrated study of publicly available datasets and microarrays, we observed a significantly elevated CDK5RAP3 expression level in NB, associated with poor patient prognosis. Further research demonstrated that CDK5RAP3 promotes the growth of NB cells, both in vitro and in vivo. Mechanistically, defective CDK5RAP3 interfered with the UFMylation system, thereby triggering endoplasmic reticulum (ER) phagy. Additionally, we provide evidence that CDK5RAP3 maintains the stability of MEIS2, a master TF in NB, and in turn, contributes to the high expression of CDK5RAP3. Overall, our findings shed light on the molecular mechanisms by which CDK5RAP3 promotes tumor progression and suggest that its inhibition may represent a novel therapeutic strategy for NB.

Reactivation of the G1 enhancer landscape underlies core circuitry addiction to SWI/SNF

Several cancer core regulatory circuitries (CRCs) depend on the sustained generation of DNA accessibility by SWI/SNF chromatin remodelers. However, the window when SWI/SNF is acutely essential in these settings has not been identified. Here we used neuroblastoma (NB) cells to model and dissect the relationship between cell-cycle progression and SWI/SNF ATPase activity. We find that SWI/SNF inactivation impairs coordinated occupancy of non-pioneer CRC members at enhancers within 1 hour, rapidly breaking their autoregulation. By precisely timing inhibitor treatment following synchronization, we show that SWI/SNF is dispensable for survival in S and G2/M, but becomes acutely essential only during G1 phase. We furthermore developed a new approach to analyze the oscillating patterns of genome-wide DNA accessibility across the cell cycle, which revealed that SWI/SNF-dependent CRC binding sites are enriched at enhancers with peak accessibility during G1 phase, where they activate genes involved in cell-cycle progression. SWI/SNF inhibition strongly impairs G1-S transition and potentiates the ability of retinoids used clinically to induce cell-cycle exit. Similar cell-cycle effects in diverse SWI/SNF-addicted settings highlight G1-S transition as a common cause of SWI/SNF dependency. Our results illustrate that deeper knowledge of the temporal patterns of enhancer-related dependencies may aid the rational targeting of addicted cancers.

Defining super-enhancers by highly ranked histone H4 multi-acetylation levels identifies transcription factors associated with glioblastoma stem-like properties

BACKGROUND

Super-enhancers (SEs), which activate genes involved in cell-type specificity, have mainly been defined as genomic regions with top-ranked enrichment(s) of histone H3 with acetylated K27 (H3K27ac) and/or transcription coactivator(s) including a bromodomain and extra-terminal domain (BET) family protein, BRD4. However, BRD4 preferentially binds to multi-acetylated histone H4, typically with acetylated K5 and K8 (H4K5acK8ac), leading us to hypothesize that SEs should be defined by high H4K5acK8ac enrichment at least as well as by that of H3K27ac.

RESULTS

Here, we conducted genome-wide profiling of H4K5acK8ac and H3K27ac, BRD4 binding, and the transcriptome by using a BET inhibitor, JQ1, in three human glial cell lines. When SEs were defined as having the top ranks for H4K5acK8ac or H3K27ac signal, 43% of H4K5acK8ac-ranked SEs were distinct from H3K27ac-ranked SEs in a glioblastoma stem-like cell (GSC) line. CRISPR-Cas9-mediated deletion of the H4K5acK8ac-preferred SEs associated with MYCN and NFIC decreased the stem-like properties in GSCs.

CONCLUSIONS

Collectively, our data highlights H4K5acK8ac's utility for identifying genes regulating cell-type specificity.

Super-enhancers complexes zoom in transcription in cancer

Super-enhancers (SEs) consist of multiple typical enhancers enriched at high density with transcription factors, histone-modifying enzymes and cofactors. Oncogenic SEs promote tumorigenesis and malignancy by altering protein-coding gene expression and noncoding regulatory element function. Therefore, they play central roles in the treatment of cancer. Here, we review the structural characteristics, organization, identification, and functions of SEs and the underlying molecular mechanism by which SEs drive oncogenic transcription in tumor cells. We then summarize abnormal SE complexes, SE-driven coding genes, and noncoding RNAs involved in tumor development. In summary, we believe that SEs show great potential as biomarkers and therapeutic targets.

SOX9 is a key component of RUNX2-regulated transcriptional circuitry in osteosarcoma

BACKGROUND

The absence of prominent, actionable genetic alternations in osteosarcomas (OS) implies that transcriptional and epigenetic mechanisms significantly contribute to the progression of this life-threatening form of cancer. Therefore, the identification of potential transcriptional events that promote the survival of OS cells could be key in devising targeted therapeutic approaches for OS. We have previously shown that RUNX2 is a transcription factor (TF) essential for OS cell survival. Unfortunately, the transcriptional network or circuitry regulated by RUNX2 in OS cells is still largely unknown.

METHODS

The TFs that are in the RUNX2 transcriptional circuitry were identified by analyzing RNAseq and ChIPseq datasets of RUNX2. To evaluate the effect of SOX9 knockdown on the survival of osteosarcoma cells in vitro, we employed cleaved caspase-3 immunoblotting and propidium iodide staining techniques. The impact of SOX9 and JMJD1C depletion on OS tumor growth was examined in vivo using xenografts and immunohistochemistry. Downstream targets of SOX9 were identified and dissected using RNAseq, pathway analysis, and gene set enrichment analysis. Furthermore, the interactome of SOX9 was identified using BioID and validated by PLA.

RESULT

Our findings demonstrate that SOX9 is a critical TF that is induced by RUNX2. Both in vitro and in vivo experiments revealed that SOX9 plays a pivotal role in the survival of OS. RNAseq analysis revealed that SOX9 activates the transcription of MYC, a downstream target of RUNX2. Mechanistically, our results suggest a transcriptional network involving SOX9, RUNX2, and MYC, with SOX9 binding to RUNX2. Moreover, we discovered that JMJD1C, a chromatin factor, is a novel binding partner of SOX9, and depletion of JMJD1C impairs OS tumor growth.

CONCLUSION

The findings of this study represent a significant advancement in our understanding of the transcriptional network present in OS cells, providing valuable insights that may contribute to the development of targeted therapies for OS.

Super-enhancer-driven lncRNA LIMD1-AS1 activated by CDK7 promotes glioma progression

Long non-coding RNAs (lncRNAs) are tissue-specific expression patterns and dysregulated in cancer. How they are regulated still needs to be determined. We aimed to investigate the functions of glioma-specific lncRNA LIMD1-AS1 activated by super-enhancer (SE) and identify the potential mechanisms. In this paper, we identified a SE-driven lncRNA, LIMD1-AS1, which is expressed at significantly higher levels in glioma than in normal brain tissue. High LIMD1-AS1 levels were significantly associated with a shorter survival time of glioma patients. LIMD1-AS1 overexpression significantly enhanced glioma cells proliferation, colony formation, migration, and invasion, whereas LIMD1-AS1 knockdown inhibited their proliferation, colony formation, migration, and invasion, and the xenograft tumor growth of glioma cells in vivo. Mechanically, inhibition of CDK7 significantly attenuates MED1 recruitment to the super-enhancer of LIMD1-AS1 and then decreases the expression of LIMD1-AS1. Most importantly, LIMD1-AS1 could directly bind to HSPA5, leading to the activation of interferon signaling. Our findings support the idea that CDK7 mediated-epigenetically activation of LIMD1-AS1 plays a crucial role in glioma progression and provides a promising therapeutic approach for patients with glioma.

Genomic and Epigenetic Changes Drive Aberrant Skeletal Muscle Differentiation in Rhabdomyosarcoma

Rhabdomyosarcoma (RMS), the most common soft-tissue sarcoma in children and adolescents, represents an aberrant form of skeletal muscle differentiation. Both skeletal muscle development, as well as regeneration of adult skeletal muscle are governed by members of the myogenic family of regulatory transcription factors (MRFs), which are deployed in a highly controlled, multi-step, bidirectional process. Many aspects of this complex process are deregulated in RMS and contribute to tumorigenesis. Interconnected loops of super-enhancers, called core regulatory circuitries (CRCs), define aberrant muscle differentiation in RMS cells. The transcriptional regulation of MRF expression/activity takes a central role in the CRCs active in skeletal muscle and RMS. In PAX3::FOXO1 fusion-positive (PF+) RMS, CRCs maintain expression of the disease-driving fusion oncogene. Recent single-cell studies have revealed hierarchically organized subsets of cells within the RMS cell pool, which recapitulate developmental myogenesis and appear to drive malignancy. There is a large interest in exploiting the causes of aberrant muscle development in RMS to allow for terminal differentiation as a therapeutic strategy, for example, by interrupting MEK/ERK signaling or by interfering with the epigenetic machinery controlling CRCs. In this review, we provide an overview of the genetic and epigenetic framework of abnormal muscle differentiation in RMS, as it provides insights into fundamental mechanisms of RMS malignancy, its remarkable phenotypic diversity and, ultimately, opportunities for therapeutic intervention.

The Functional Role and Regulatory Mechanism of FTO m6A RNA Demethylase in Human Uterine Leiomyosarcoma

Uterine leiomyosarcoma (uLMS) is the most frequent subtype of uterine sarcoma that presents a poor prognosis and high rates of recurrence and metastasis. The origin and molecular mechanism underlying and driving its clinical and biological behavior remain largely unknown. Recently, we and others have revealed the role of microRNAs, DNA methylation, and histone modifications in contributing to the pathogenesis of uLMS. However, the connection between reversible m6A RNA methylation and uLMS pathogenesis remains unclear. In this study, we assessed the role and mechanism of FTO m6A RNA demethylase in the pathogenesis of uLMS. Immunohistochemistry analysis revealed that the levels of RNA demethylases FTO and ALKBH5 were aberrantly upregulated in uLMS tissues compared to adjacent myometrium with a significant change by histochemical scoring assessment (p < 0.01). Furthermore, the inhibition of FTO demethylase with its small, potent inhibitor (Dac51) significantly decreased the uLMS proliferation dose-dependently via cell cycle arrest. Notably, RNA-seq analysis revealed that the inhibition of FTO with Dac51 exhibited a significant decrease in cell-cycle-related genes, including several CDK members, and a significant increase in the expression of CDKN1A, which correlated with a Dac51-exerted inhibitory effect on cell proliferation. Moreover, Dac51 treatment allowed the rewiring of several critical pathways, including TNFα signaling, KRAS signaling, inflammation response, G2M checkpoint, and C-Myc signaling, among others, leading to the suppression of the uLMS phenotype. Moreover, transcription factor (TF) analyses suggested that epitranscriptional alterations by Dac51 may alter the cell cycle-related gene expression via TF-driven pathways and epigenetic networks in uLMS cells. This intersection of RNA methylation and other epigenetic controls and pathways provides a framework to better understand uterine diseases, particularly uLMS pathogenesis with a dysregulation of RNA methylation machinery. Therefore, targeting the vulnerable epitranscriptome may provide an additional regulatory layer for a promising and novel strategy for treating patients with this aggressive uterine cancer.

Profiling and integrated analysis of transcriptional addiction gene expression and prognostic value in hepatocellular carcinoma

Transcriptional dysregulation caused by genomic and epigenetic alterations in cancer is called "transcriptional addiction". Transcriptional addiction is an important pathogenic factor of tumor malignancy. Hepatocellular carcinoma (HCC) genomes are highly heterogeneous, with many dysregulated genes. Our study analyzed the possibility that transcriptional addiction-related genes play a significant role in HCC. All data sources for conducting this study were public cancer databases and tissue microarrays. We identified 38 transcriptional addiction genes, and most were differentially expressed genes. Among patients of different groups, there were significant differences in overall survival rates. Both nomogram and risk score were independent predictors of HCC outcomes. Transcriptional addiction gene expression characteristics determine the sensitivity of patients to immunotherapy, cisplatin, and sorafenib. Besides, HDAC2 was identified as an oncogene, and its expression was correlated with patient survival time. Our study conclusively demonstrated that transcriptional addiction is crucial in HCC. We provided biomarkers for predicting the prognosis of HCC patients, which can more precisely guide the patient's treatment.

SWI/SNF Blockade Disrupts PU.1-Directed Enhancer Programs in Normal Hematopoietic Cells and Acute Myeloid Leukemia

In acute myeloid leukemia (AML), SWI/SNF chromatin remodeling complexes sustain leukemic identity by driving high levels of MYC. Previous studies have implicated the hematopoietic transcription factor PU.1 (SPI1) as an important target of SWI/SNF inhibition, but PU.1 is widely regarded to have pioneer-like activity. As a result, many questions have remained regarding the interplay between PU.1 and SWI/SNF in AML as well as normal hematopoiesis. Here we found that PU.1 binds to most of its targets in a SWI/SNF-independent manner and recruits SWI/SNF to promote accessibility for other AML core regulatory factors, including RUNX1, LMO2, and MEIS1. SWI/SNF inhibition in AML cells reduced DNA accessibility and binding of these factors at PU.1 sites and redistributed PU.1 to promoters. Analysis of nontumor hematopoietic cells revealed that similar effects also impair PU.1-dependent B-cell and monocyte populations. Nevertheless, SWI/SNF inhibition induced profound therapeutic response in an immunocompetent AML mouse model as well as in primary human AML samples. In vivo, SWI/SNF inhibition promoted leukemic differentiation and reduced the leukemic stem cell burden in bone marrow but also induced leukopenia. These results reveal a variable therapeutic window for SWI/SNF blockade in AML and highlight important off-tumor effects of such therapies in immunocompetent settings.

SIGNIFICANCE

Disruption of PU.1-directed enhancer programs upon SWI/SNF inhibition causes differentiation of AML cells and induces leukopenia of PU.1-dependent B cells and monocytes, revealing the on- and off-tumor effects of SWI/SNF blockade.

Expression of paired box 9 defines an aggressive subset of lung adenocarcinoma preferentially occurring in smokers

AIMS

A distinct subset of lung adenocarcinomas (LADs), arising from a series of peripheral lung cells defined as the terminal respiratory unit (TRU), is characterised by thyroid transcription factor 1 (TTF-1) expression. The clinical relevance of transcription factors (TFs) other than TTF-1 remains unknown in LAD and was explored in the present study.

METHODS AND RESULTS

Seventy-one LAD samples were subjected to high-throughput transcriptome screening of LAD using cap analysis gene expression (CAGE) sequencing data; CAGE provides genome-wide expression levels of the transcription start sites (TSSs). In total, 1083 invasive LAD samples were subjected to immunohistochemical examination for paired box 9 (PAX9) and TTF-1 expression levels. PAX9 is an endoderm development-associated TF that most strongly and inversely correlates with the expression of TTF-1 TSS subsets. Immunohistochemically, PAX9 expression was restricted to the nuclei of ciliated epithelial and basal cells in the bronchi and bronchioles and the nuclei of epithelial cells of the bronchial glands; moreover, PAX9 expression was observed in 304 LADs (28%). PAX9-positive LADs were significantly associated with heavy smoking, non-lepidic subtype, EGFR wild-type tumours and PD-L1 expression (all P < 0.0001). All these characteristics were opposite to those of TRU-type LADs with TTF-1 expression. PAX9 expression was an independent prognostic factor for decreased overall survival (P = 0.022).

CONCLUSIONS

Our results revealed that PAX9 expression defines an aggressive subset of LADs preferentially occurring in smokers that may arise from bronchial or bronchiolar cells.

Leukemia core transcriptional circuitry is a sparsely interconnected hierarchy stabilized by incoherent feed-forward loops

Lineage-defining transcription factors form densely interconnected circuits in chromatin occupancy assays, but the functional significance of these networks remains underexplored. We reconstructed the functional topology of a leukemia cell transcription network from the direct gene-regulatory programs of eight core transcriptional regulators established in pre-steady state assays coupling targeted protein degradation with nascent transcriptomics. The core regulators displayed narrow, largely non-overlapping direct transcriptional programs, forming a sparsely interconnected functional hierarchy stabilized by incoherent feed-forward loops. BET bromodomain and CDK7 inhibitors disrupted the core regulators' direct programs, acting as mixed agonists/antagonists. The network is predictive of dynamic gene expression behaviors in time-resolved assays and clinically relevant pathway activity in patient populations.

Transcriptional regulatory mechanisms and signaling networks in cancer

Cancer is a general term that refers to a wide range of illnesses that are characterized by the development of aberrant cells that have the capacity to divide uncontrollably, invade, and harm healthy tissue. It is caused by both genetic and epigenetic changes that suppress abnormal proliferation and prevent cells from surviving outside of their normal niches. Complex protein networks are responsible for the development of a suitable environment via multiple cells signaling pathways. The study of these pathways is essential for analysing network context and developing novel cancer therapies. Transcription factors (TFs) are actively involved in gene expression and maintain the combinatorial on-and-off states of the gene. In addition, the TFs regulate cell identity and state; these TFs cooperate to establish cell-type-specific gene expression. In this chapter, we describe the number of transcription factors and their role in the progression of cancer. The knowledge of transcriptional factors and their network is crucial for emphasizing the specific transcriptional addiction and for designing new anticancer therapies.

TRANSPARENT: a Python tool for designing transcription factor regulatory networks

Transcription factors are proteins able to selectively bind DNA short traits, namely transcription factors binding sites, in order to regulate gene expression in terms of both repression and activation. Despite plenty of studies focusing on transcription factors and on the role they play in specific biological tasks or diseases, is available in the literature, to our knowledge there is no tool able to automatically provide a list of transcription factors involved in this task and the associated interaction network through a solid computational analysis. TRANScriPtion fActor REgulatory NeTwork (TRANSPARENT) is a user-friendly Python tool designed to help researchers in studying given biological tasks or given diseases in human, by identifying transcription factors controlling and regulating the expression of genes associated with that task or disease. The tool takes in input a list of gene IDs and provides (1) a set of transcription factors that are significantly associated with the input genes, (2) the correspondent P values (i.e., the probability that this observed association was driven by chance) and (3) a transcription factor network that can be directly visualized through STRING database. The effectiveness and reliability of the tool were assessed by applying it to two different test cases: schizophrenia and autism disorders. The obtained results clearly show that identified TFs, for both datasets, are significantly associated with those disorders, in terms of both gene enrichment and coherence with the literature. The proposed tool TRANSPARENT can be a useful instrument to investigate transcription factor networks and unveil the role that TFs play in given biological tasks and diseases.

DeepCIA: a novel deep-learning model for cancer type identification using class activation map via transcription factor expression

Deep learning methods are powerful analytical tools for large-scale data analysis. Here, we introduce DeepCIA as a novel diagnostic deep-learning model for cancer type identification using a class activation map via transcription factor expression. Although many deep learning researches attempts have recently been made in relation to cancer diagnosis, there are difficulties in using cancer data due to a large-scale problem. Therefore, From The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) public databases, we selected transcription factor expression profiles of eight cancer types. TCGA included 3496 samples and divided the train and validation sets in an 8:2 ratio. ICGC included 552 samples and was used as a test set for external validation. To compare the performance of 1D-CNN models, we also used SVM and KNN from machine learning. In external validation, 1D-CNN showed a high average accuracy of 98% and was superior to support vector machine (SVM) and k-nearest neighbor (KNN) with a difference in the accuracy of 10-12%. Also, 1D-CNN performed very well in several performance metrics (98.2% Recall, 98.1% Precision, 98.2% F score, 99.8% Specificity, 99.8% AUC, and 99.0% Balanced Accuracy). In each data set evaluation, 1-network, 5-network, and 2-network with high accuracy were selected and visualized through the Class Activation Map. We identified the Cys2Hys2 zinc finger group with the highest distribution across all cancer types. Collectively, DeepCIA can be used as a decision support system for cancer and a classifier for diagnosing unknown primary cancer, while emphasizing its usefulness in cancer diagnosis.

Targeting Class I Histone Deacetylases in Human Uterine Leiomyosarcoma

Uterine leiomyosarcoma (uLMS) is the most frequent subtype of uterine sarcoma that presents a poor prognosis, high rates of recurrence, and metastasis. Currently, the molecular mechanism of the origin and development of uLMS is unknown. Class I histone deacetylases (including HDAC1, 2, 3, and 8) are one of the major classes of the HDAC family and catalyze the removal of acetyl groups from lysine residues in histones and cellular proteins. Class I HDACs exhibit distinct cellular and subcellular expression patterns and are involved in many biological processes and diseases through diverse signaling pathways. However, the link between class I HDACs and uLMS is still being determined. In this study, we assessed the expression panel of Class I HDACs in uLMS and characterized the role and mechanism of class I HDACs in the pathogenesis of uLMS. Immunohistochemistry analysis revealed that HDAC1, 2, and 3 are aberrantly upregulated in uLMS tissues compared to adjacent myometrium. Immunoblot analysis demonstrated that the expression levels of HDAC 1, 2, and 3 exhibited a graded increase from normal and benign to malignant uterine tumor cells. Furthermore, inhibition of HDACs with Class I HDACs inhibitor (Tucidinostat) decreased the uLMS proliferation in a dose-dependent manner. Notably, gene set enrichment analysis of differentially expressed genes (DEGs) revealed that inhibition of HDACs with Tucidinostat altered several critical pathways. Moreover, multiple epigenetic analyses suggested that Tucidinostat may alter the transcriptome via reprogramming the oncogenic epigenome and inducing the changes in microRNA-target interaction in uLMS cells. In the parallel study, we also determined the effect of DL-sulforaphane on the uLMS. Our study demonstrated the relevance of class I HDACs proteins in the pathogenesis of malignant uLMS. Further understanding the role and mechanism of HDACs in uLMS may provide a promising and novel strategy for treating patients with this aggressive uterine cancer.

Dissecting super-enhancer driven transcriptional dependencies reveals novel therapeutic strategies and targets for group 3 subtype medulloblastoma

BACKGROUND

Medulloblastoma is the most common malignant pediatric brain tumor and group 3 subtype medulloblastoma (G3-MB) exhibits the worst prognosis. Super enhancers (SEs) are large clusters of enhancers that play important roles in cancer through transcriptional control of cell identity genes, oncogenes and tumor-dependent genes. Dissecting SE-driven transcriptional dependencies of cancer leads to identification of novel oncogenic mechanisms, therapeutic strategies and targets.

METHODS

Integrative SE analyses of primary tissues and patient-derived tumor cell lines of G3-MB were performed to extract the conserved SE-associated gene signatures and their oncogenic potentials were evaluated by gene expression, tumor-dependency and patient prognosis analyses. SE-associated subtype-specific upregulated tumor-dependent genes, which were revealed as members of SE-driven core transcriptional regulatory network of G3-MB, were then subjected to functional validation and mechanistic investigation. SE-associated therapeutic potential was further explored by genetic or pharmaceutical targeting of SE complex components or SE-associated subtype-specific upregulated tumor-dependent genes individually or in combination, and the underlying therapeutic mechanisms were also examined.

RESULTS

The identified conserved SE-associated transcripts of G3-MB tissues and cell lines were enriched of subtype-specifically upregulated tumor-dependent genes and MB patients harboring enrichment of those transcripts exhibited worse prognosis. Fourteen such conserved SE-associated G3-MB-specific upregulated tumor-dependent genes were identified to be members of SE-driven core transcriptional regulatory network of G3-MB, including three well-recognized TFs (MYC, OTX2 and CRX) and eleven newly identified downstream effector genes (ARL4D, AUTS2, BMF, IGF2BP3, KIF21B, KLHL29, LRP8, MARS1, PSMB5, SDK2 and SSBP3). An OTX2-SE-ARL4D regulatory axis was further revealed to represent a subtype-specific tumor dependency and therapeutic target of G3-MB via contributing to maintaining cell cycle progression and inhibiting neural differentiation of tumor cells. Moreover, BET inhibition with CDK7 inhibition or proteasome inhibition, two combinatory strategies of targeting SE complex components (BRD4, CDK7) or SE-associated effector gene (PSMB5), were shown to exhibit synergistic therapeutic effects against G3-MB via stronger suppression of SE-associated transcription or higher induction of ER stress, respectively.

CONCLUSIONS

Our study verifies the oncogenic role and therapeutic potential of SE-driven transcriptional dependencies of G3-MB, resulting in better understanding of its tumor biology and identification of novel SE-associated therapeutic strategies and targets.

Recurrent noncoding somatic and germline WT1 variants converge to disrupt MYB binding in acute promyelocytic leukemia

Genetic alternations can occur at noncoding regions, but how they contribute to cancer pathogenesis is poorly understood. Here, we established a mutational landscape of cis-regulatory regions (CREs) in acute promyelocytic leukemia (APL) based on whole-genome sequencing analysis of paired tumor and germline samples from 24 patients and epigenetic profiling of 16 patients. Mutations occurring in CREs occur preferentially in active enhancers bound by the complex of master transcription factors in APL. Among significantly enriched mutated CREs, we found a recurrently mutated region located within the third intron of WT1, an essential regulator of normal and malignant hematopoiesis. Focusing on noncoding mutations within this WT1 intron, an analysis on 169 APL patients revealed that somatic mutations were clustered into a focal hotspot region, including one site identified as a germline polymorphism contributing to APL risk. Significantly decreased WT1 expression was observed in APL patients bearing somatic and/or germline noncoding WT1 variants. Furthermore, biallelic WT1 inactivation was recurrently found in APL patients with noncoding WT1 variants, which resulted in the complete loss of WT1. The high incidence of biallelic inactivation suggested the tumor suppressor activity of WT1 in APL. Mechanistically, noncoding WT1 variants disrupted MYB binding on chromatin and suppressed the enhancer activity and WT1 expression through destroying the chromatin looping formation. Our study highlights the important role of noncoding variants in the leukemogenesis of APL.

Differential Gene Expression and Weighted Correlation Network Dynamics in High-Throughput Datasets of Prostate Cancer

Precision oncology is an absolute need today due to the emergence of treatment resistance and heterogeneity among cancerous profiles. Target-propelled cancer therapy is one of the treasures of precision oncology which has come together with substantial medical accomplishment. Prostate cancer is one of the most common cancers in males, with tremendous biological heterogeneity in molecular and clinical behavior. The spectrum of molecular abnormalities and varying clinical patterns in prostate cancer suggest substantial heterogeneity among different profiles. To identify novel therapeutic targets and precise biomarkers implicated with prostate cancer, we performed a state-of-the-art bioinformatics study, beginning with analyzing high-throughput genomic datasets from The Cancer Genome Atlas (TCGA). Weighted gene co-expression network analysis (WGCNA) suggests a set of five dysregulated hub genes (MAF, STAT6, SOX2, FOXO1, and WNT3A) that played crucial roles in biological pathways associated with prostate cancer progression. We found overexpressed STAT6 and SOX2 and proposed them as candidate biomarkers and potential targets in prostate cancer. Furthermore, the alteration frequencies in STAT6 and SOX2 and their impact on the patients' survival were explored through the cBioPortal platform. The Kaplan-Meier survival analysis suggested that the alterations in the candidate genes were linked to the decreased overall survival of the patients. Altogether, the results signify that STAT6 and SOX2 and their genomic alterations can be explored in therapeutic interventions of prostate cancer for precision oncology, utilizing early diagnosis and target-propelled therapy.

PRRX1 is a master transcription factor of stromal fibroblasts for myofibroblastic lineage progression

Although stromal fibroblasts play a critical role in cancer progression, their identities remain unclear as they exhibit high heterogeneity and plasticity. Here, a master transcription factor (mTF) constructing core-regulatory circuitry, PRRX1, which determines the fibroblast lineage with a myofibroblastic phenotype, is identified for the fibroblast subgroup. PRRX1 orchestrates the functional drift of fibroblasts into myofibroblastic phenotype via TGF-β signaling by remodeling a super-enhancer landscape. Such reprogrammed fibroblasts have myofibroblastic functions resulting in markedly enhanced tumorigenicity and aggressiveness of cancer. PRRX1 expression in cancer-associated fibroblast (CAF) has an unfavorable prognosis in multiple cancer types. Fibroblast-specific PRRX1 depletion induces long-term and sustained complete remission of chemotherapy-resistant cancer in genetically engineered mice models. This study reveals CAF subpopulations based on super-enhancer profiles including PRRX1. Therefore, mTFs, including PRRX1, provide another opportunity for establishing a hierarchical classification system of fibroblasts and cancer treatment by targeting fibroblasts.

Cancer associated fibroblasts are an important and highly heterogeneous component of the tumor microenvironment. Here the authors identify PRRX1 as a master transcription factor determining a fibroblast lineage with myofibroblastic phenotype, associated with unfavourable prognosis in several cancer types.

Intrinsic transcriptional heterogeneity in neuroblastoma guides mechanistic and therapeutic insights

Cell state is controlled by master transcription factors (mTFs) that determine the cellular gene expression program. Cancer cells acquire dysregulated gene expression programs by mutational and non-mutational processes. Intratumoral heterogeneity can result from cells displaying distinct mTF-regulated cell states, which co-exist within the tumor. One archetypal tumor associated with transcriptionally regulated heterogeneity is high-risk neuroblastoma (NB). Patients with NB have poor overall survival despite intensive therapies, and relapsed patients are commonly refractory to treatment. The cellular populations that comprise NB are marked by different cohorts of mTFs and differential sensitivity to conventional therapies. Recent studies have highlighted mechanisms by which NB cells dynamically shift the cell state with treatment, revealing new opportunities to control the cellular response to treatment by manipulating cell-state-defining transcriptional programs. Here, we review recent advances in understanding transcriptionally defined cancer heterogeneity. We offer challenges to the field to encourage translation of basic science into clinical benefit.

Identifying common transcriptome signatures of cancer by interpreting deep learning models

Background

Cancer is a set of diseases characterized by unchecked cell proliferation and invasion of surrounding tissues. The many genes that have been genetically associated with cancer or shown to directly contribute to oncogenesis vary widely between tumor types, but common gene signatures that relate to core cancer pathways have also been identified. It is not clear, however, whether there exist additional sets of genes or transcriptomic features that are less well known in cancer biology but that are also commonly deregulated across several cancer types.

Results

Here, we agnostically identify transcriptomic features that are commonly shared between cancer types using 13,461 RNA-seq samples from 19 normal tissue types and 18 solid tumor types to train three feed-forward neural networks, based either on protein-coding gene expression, lncRNA expression, or splice junction use, to distinguish between normal and tumor samples. All three models recognize transcriptome signatures that are consistent across tumors. Analysis of attribution values extracted from our models reveals that genes that are commonly altered in cancer by expression or splicing variations are under strong evolutionary and selective constraints. Importantly, we find that genes composing our cancer transcriptome signatures are not frequently affected by mutations or genomic alterations and that their functions differ widely from the genes genetically associated with cancer.

Conclusions

Our results highlighted that deregulation of RNA-processing genes and aberrant splicing are pervasive features on which core cancer pathways might converge across a large array of solid tumor types.

Supplementary Information

The online version contains supplementary material available at (10.1186/s13059-022-02681-3).

POU6F1 cooperates with RORA to suppress the proliferation of lung adenocarcinoma by downregulation HIF1A signaling pathway

Lung adenocarcinoma (LUAD) represents the most frequently diagnosed histological subtype of non-small cell lung cancer with the highest mortality worldwide. Transcriptional dysregulation is a hallmark of nearly all kinds of cancers. In the study, we identified that the POU domain, class 6, transcription factor 1 (POU6F1), a member of the POU family of transcription factors, was closely associated with tumor stage and death in LUAD. We revealed that POU6F1 was downregulated in LUAD tissues and downregulated POU6F1 was predictive of an unfavorable prognosis in LUAD patients. In vitro assays, including CCK8, soft agar, transwell, clone formation, wound-healing assay, and nude mouse xenograft model all revealed that POU6F1 inhibited the growth and invasion of LUAD cells. Mechanistically, POU6F1 bound and stabilized retinoid-related orphan receptor alpha (RORA) to exert the transcriptional inhibition of hypoxia-inducible factor 1-alpha (HIF1A) and alter the expression of HIF1A signaling pathway-associated genes, including ENO1, PDK1, and PRKCB, thereby leading to the suppression of LUAD cells. Collectively, these results demonstrated the suppressive role of POU6F1/RORA in the progression of LUAD and may potentially be used as a target for the treatment of LUAD.

A feedforward circuit between KLF5 and lncRNA KPRT4 contributes to basal-like breast cancer

Basal-like breast cancer (BLBC) is the most aggressive subtype of breast cancer with a poor prognosis. Long noncoding RNAs (lncRNAs) play critical roles in human cancers. Krüppel-like Factor 5 (KLF5) is a key oncogenic transcription factor in BLBC. However, the underlying mechanism of mutual regulation between KLF5 and lncRNA remains largely unknown. Here, we demonstrate that lncRNA KPRT4 promotes BLBC cell proliferation in vitro and in vivo. Mechanistically, KLF5 directly binds to the promoter of KPRT4 to promote KPRT4 transcription. Reciprocally, KPRT4 recruits the YB-1 transcription factor to the KLF5 promoter by interacting with YB-1 at its 5' domain and forming an RNA-DNA-DNA triplex structure at its 3' domain, resulting in enhanced transcription of KLF5 and ultimately establishing a feedforward circuit to promote cell proliferation. Moreover, the antisense oligonucleotide (ASO)-based therapy targeting KPRT4 substantially attenuated tumor growth in vivo. Clinically, the expression levels of YB-1, KLF5 and KPRT4 are positively correlated in clinical breast specimens. Together, our data suggest that KPRT4 is a major molecule for BLBC progression and that the feedforward circuit between KLF5 and KPRT4 may represent a potential therapeutic target in BLBC.

A distinct core regulatory module enforces oncogene expression in KMT2A-rearranged leukemia

In this study, Harada et al. identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. This study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Acute myeloid leukemia with KMT2A (MLL) rearrangements is characterized by specific patterns of gene expression and enhancer architecture, implying unique core transcriptional regulatory circuitry. Here, we identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. Rapid transcription factor degradation followed by measurements of genome-wide transcription rates and superresolution microscopy revealed that MEF2D and IRF8 form a distinct core regulatory module with a narrow direct transcriptional program that includes activation of the key oncogenes MYC, HOXA9, and BCL2. Our study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Distinct properties of pure- and mixed-type high-grade fetal lung adenocarcinomas by genetic profiling and transcription factor expression

The clinicopathological differences among high-grade fetal lung adenocarcinomas completely comprising tumor cells that resemble fetal lung epithelium (pure type) and those with fetal lung-like components admixed with conventional adenocarcinoma cells (mixed type) remain undetermined. Here, we examined the clinicopathological, immunohistochemical, and molecular features of 11 lung adenocarcinomas with fetal lung-like morphology among 3895 consecutive cases of primary lung cancer based on the expression pattern of transcription factors. According to the current WHO classification, two cases (0.05%) were categorized as low-grade fetal adenocarcinoma, two cases (0.05%) were pure-type high-grade fetal adenocarcinoma, five cases (0.1%) were mixed-type high-grade fetal adenocarcinoma, and the remaining two cases (0.05%) were lung adenocarcinoma with high-grade fetal features (fetal lung-like morphology occupied less than 50%). CTNNB1 mutations were exclusively identified in low-grade fetal adenocarcinomas. In contrast, mixed-type high-grade fetal adenocarcinoma or lung adenocarcinoma with high-grade fetal features frequently harbored mitogenic drivers including EGFR mutations. Furthermore, almost all tumor cells expressed CDX2 and HNF4α in both cases of pure-type high-grade fetal lung adenocarcinoma, but lacked TTF-1 positivity. In contrast, TTF-1 was frequently expressed in mixed-type high-grade fetal lung adenocarcinoma and in lung adenocarcinoma with high-grade fetal features. Our data suggest similar prevalence of low-grade fetal lung adenocarcinoma and pure-type high-grade fetal lung adenocarcinoma, and indicate that pure- and mixed-type high-grade fetal lung adenocarcinomas are distinct, with the former akin to low-grade fetal adenocarcinoma with respect to purely embryonic morphology and absence of common lung adenocarcinoma mitogenic drivers, and the latter being genetically and transcriptionally related to conventional lung adenocarcinoma.

Targeting Oncogenic Transcriptional Networks in Neuroblastoma: From N-Myc to Epigenetic Drugs

Neuroblastoma (NB) is one of the most frequently occurring neurogenic extracranial solid cancers in childhood and infancy. Over the years, many pieces of evidence suggested that NB development is controlled by gene expression dysregulation. These unleashed programs that outline NB cancer cells make them highly dependent on specific tuning of gene expression, which can act co-operatively to define the differentiation state, cell identity, and specialized functions. The peculiar regulation is mainly caused by genetic and epigenetic alterations, resulting in the dependency on a small set of key master transcriptional regulators as the convergence point of multiple signalling pathways. In this review, we provide a comprehensive blueprint of transcriptional regulation bearing NB initiation and progression, unveiling the complexity of novel oncogenic and tumour suppressive regulatory networks of this pathology. Furthermore, we underline the significance of multi-target therapies against these hallmarks, showing how novel approaches, together with chemotherapy, surgery, or radiotherapy, can have substantial antineoplastic effects, disrupting a wide variety of tumorigenic pathways through combinations of different treatments.

Loss of synergistic transcriptional feedback loops drives diverse B-cell cancers

Background

The most common B-cell cancers, chronic lymphocytic leukemia/lymphoma (CLL), follicular and diffuse large B-cell (FL, DLBCL) lymphomas, have distinct clinical courses, yet overlapping “cell-of-origin”. Dynamic changes to the epigenome are essential regulators of B-cell differentiation. Therefore, we reasoned that these distinct cancers may be driven by shared mechanisms of disruption in transcriptional circuitry.

Methods

We compared purified malignant B-cells from 52 patients with normal B-cell subsets (germinal center centrocytes and centroblasts, naïve and memory B-cells) from 36 donor tonsils using >325 high-resolution molecular profiling assays for histone modifications, open chromatin (ChIP-, FAIRE-seq), transcriptome (RNA-seq), transcription factor (TF) binding, and genome copy number (microarrays).

Findings

From the resulting data, we identified gains in active chromatin in enhancers/super-enhancers that likely promote unchecked B-cell receptor signaling, including one we validated near the immunoglobulin superfamily receptors FCMR and PIGR. More striking and pervasive was the profound loss of key B-cell identity TFs, tumor suppressors and their super-enhancers, including EBF1, OCT2(POU2F2), and RUNX3. Using a novel approach to identify transcriptional feedback, we showed that these core transcriptional circuitries are self-regulating. Their selective gain and loss form a complex, iterative, and interactive process that likely curbs B-cell maturation and spurs proliferation.

Interpretation

Our study is the first to map the transcriptional circuitry of the most common blood cancers. We demonstrate that a critical subset of B-cell TFs and their cognate enhancers form self-regulatory transcriptional feedback loops whose disruption is a shared mechanism underlying these diverse subtypes of B-cell lymphoma.

Funding

National Institute of Health, Siteman Cancer Center, Barnes-Jewish Hospital Foundation, Doris Duke Foundation.

Reporter gene systems for the identification and characterization of cancer stem cells

Cancer stem cells (CSCs) are tumor cells that share functional characteristics with normal and embryonic stem cells. CSCs have increased tumor-initiating capacity and metastatic potential and lower sensitivity to chemo- and radiotherapy, with important roles in tumor progression and the response to therapy. Thus, a current goal of cancer research is to eliminate CSCs, necessitating an adequate phenotypic and functional characterization of CSCs. Strategies have been developed to identify, enrich, and track CSCs, many of which distinguish CSCs by evaluating the expression of surface markers, the initiation of specific signaling pathways, and the activation of master transcription factors that control stemness in normal cells. We review and discuss the use of reporter gene systems for identifying CSCs. Reporters that are under the control of aldehyde dehydrogenase 1A1, CD133, Notch, Nanog homeobox, Sex-determining region Y-box 2, and POU class 5 homeobox can be used to identify CSCs in many tumor types, track cells in real time, and screen for drugs. Thus, reporter gene systems, in combination with in vitro and in vivo functional assays, can assess changes in the CSCs pool. We present relevant examples of these systems in the evaluation of experimental CSCs-targeting therapeutics, demonstrating their value in CSCs research.

Super-enhancer-mediated core regulatory circuitry in human cancer

Super-enhancers (SEs) are congregated enhancer clusters with high level of loading of transcription factors (TFs), cofactors and epigenetic modifications. Through direct co-occupancy at their own SEs as well as each other's, a small set of so called "master" TFs form interconnected core regulatory circuitry (CRCs) to orchestrate transcriptional programs in both normal and malignant cells. These master TFs can be predicted mathematically using epigenomic methods. In this Review, we summarize the identification of SEs and CRCs in cancer cells, the mechanisms by which master TFs and SEs cooperatively regulate cancer-type-specific expression programs, and the cancer-type- and subtype-specificity of CRC and the significance in cancer biology.

CDK7 Inhibitor THZ1 Induces the Cell Apoptosis of B-Cell Acute Lymphocytic Leukemia by Perturbing Cellular Metabolism

B-cell acute lymphocytic leukemia (B-ALL) is a malignant blood cancer that develops in children and adults and leads to high mortality. THZ1, a covalent cyclin-dependent kinase 7 (CDK7) inhibitor, shows anti-tumor effects in various cancers by inhibiting cell proliferation and inducing apoptosis. However, whether THZ1 has an inhibitory effect on B-ALL cells and the underlying mechanism remains obscure. In this study, we showed that THZ1 arrested the cell cycle of B-ALL cells in vitro in a low concentration, while inducing the apoptosis of B-ALL cells in vitro in a high concentration by activating the apoptotic pathways. In addition, RNA-SEQ results revealed that THZ1 disrupted the cellular metabolic pathways of B-ALL cells. Moreover, THZ1 suppressed the cellular metabolism and blocked the production of cellular metabolic intermediates in B-ALL cells. Mechanistically, THZ1 inhibited the cellular metabolism of B-ALL by downregulating the expression of c-MYC-mediated metabolic enzymes. However, THZ1 treatment enhanced cell apoptosis in over-expressed c-MYC B-ALL cells, which was involved in the upregulation of p53 expression. Collectively, our data demonstrated that CDK7 inhibitor THZ1 induced the apoptosis of B-ALL cells by perturbing c-MYC-mediated cellular metabolism, thereby providing a novel treatment option for B-ALL.

Topography of transcriptionally active chromatin in glioblastoma

Molecular profiling of the most aggressive brain tumor glioblastoma (GBM) on the basis of gene expression, DNA methylation, and genomic variations advances both cancer research and clinical diagnosis. The enhancer architectures and regulatory circuitries governing tumor-intrinsic transcriptional diversity and subtype identity are still elusive. Here, by mapping H3K27ac deposition, we analyze the active regulatory landscapes across 95 GBM biopsies, 12 normal brain tissues, and 38 cell line counterparts. Analyses of differentially regulated enhancers and super-enhancers uncovered previously unrecognized layers of intertumor heterogeneity. Integrative analysis of variant enhancer loci and transcriptome identified topographies of transcriptional enhancers and core regulatory circuitries in four molecular subtypes of primary tumors: AC1-mesenchymal, AC1-classical, AC2-proneural, and AC3-proneural. Moreover, this study reveals core oncogenic dependency on super-enhancer-driven transcriptional factors, long noncoding RNAs, and druggable targets in GBM. Through profiling of transcriptional enhancers, we provide clinically relevant insights into molecular classification, pathogenesis, and therapeutic intervention of GBM.

Tribbles homolog 2 (Trib2), a pseudo serine/threonine kinase in tumorigenesis and stem cell fate decisions

The family of Tribbles proteins play many critical nonenzymatic roles and regulate a wide range of key signaling pathways. Tribbles homolog 2 (Trib2) is a pseudo serine/threonine kinase that functions as a scaffold or adaptor in various physiological and pathological processes. Trib2 can interact with E3 ubiquitin ligases and control protein stability of downstream effectors. This protein is induced by mitogens and enhances the propagation of several cancer cells, including myeloid leukemia, liver, lung, skin, bone, brain, and pancreatic. Thus, Trib2 can be a predictive and valuable biomarker for the diagnosis and treatment of cancer. Recent studies have illustrated that Trib2 plays a major role in cell fate determination of stem cells. Stem cells have the capacity to self-renew and differentiate into specific cell types. Stem cells are important sources for cell-based regenerative medicine and drug screening. Trib2 has been found to increase the self-renewal ability of embryonic stem cells, the reprogramming efficiency of somatic cells, and chondrogenesis. In this review, we will focus on the recent advances of Trib2 function in tumorigenesis and stem cell fate decisions.