B Cell Transcriptional Coactivator POU2AF1 (BOB-1) Is an Early Transcription Factor Modulating the Protein Synthesis and Ribosomal Biogenesis in Multiple Myeloma: With Therapeutic Implication

Multiple Myeloma (MM) is a malignancy driven by numerous genetic and epigenetic alterations. Recurrent IgH translocations, somatic mutations and copy number abnormalities all contribute to myelomagenesis, however true drivers of the disease have not been well defined. To identify new targetable dependencies in MM, we generated high-quality active enhancer landscape using large cohort of primary patient myeloma cells (n=70), MM cell lines, and normal plasma cells. We integrated this data with an in-house curated atlas of 600+ active enhancer profile across a wide range of tumor types and normal tissues. Combining these data with gene expression and genetic dependency (CRISPR KO) enabled a multidimensional integration of how transcriptional regulation intersects with tumor specific dependencies. We identified that many of the specific and potent dependencies in MM are transcription factors, especially those establishing plasma cell identity. Among these, the POU2AF1 gene, which encodes the OCA-B/BOB-1, a B cell transcriptional coactivator protein, represented the most striking dependency in MM. Although BOB-1 is expressed throughout B-cell development, we found it to be highly expressed in CD138+ plasma cells from patients with precursor conditions (MGUS and SMM) as well as symptomatic MM compared to normal plasma cells. To functionally validate the role of BOB-1 in MM, we performed loss-of-function studies using shRNA, siRNAs and antisense GapMers specific for BOB-1 and observed a significant impact on MM cell viability and cell cycle arrest. Transcriptomic analysis upon BOB-1 depletion by RNA-sequencing revealed a small set of genes commonly modulated in all 3 MM cell lines tested including the plasma cell differentiation related transcription factor XBP1 and heme oxygenase (HMOX1). Importantly, we observed ribosome biogenesis, RNA polymerase 1A transcription and mRNA translation and elongation processes to be significantly enriched among genes modulated by BOB-1 depletion in MM cells. Bob1 KD resulted in a rapid and robust decrease in the level of transcription of rDNA by RNA polymerase I as determined by qRT-PCR quantification of pre-rRNA (47S). In addition, ChiP assay revealed decreased binding of RNA polymerase 1A to the 18S ribosomal DNA promoter region in BOB-1 depleted cells compared to control. These data indicate that BOB1 downregulation results in the suppression of RNA-polymerase I activity in MM cells. RNA Pol I-dependent transcription governs abundance of rRNA and directly regulates cellular translational and proliferative capacity. Since high protein load is a feature of MM, we evaluated the role of BOB-1 in the translational efficiency of MM cells. We observed that in MM cells compared to control cells, BOB-1 KD decreased, while its overexpression significantly enhanced de novo protein synthesis. As MM is characterized by excess production of monoclonal immunoglobulins, we evaluated impact of BOB-1 perturbation on intracellular light chains (kappa or lambda) production. We observed changes in the intracellular abundance of the light chains with BOB-1 modulation in all MM cell lines tested. As a result of decreased protein production, BOB-1 depletion was associated with induction of resistance to proteasome inhibition suggesting that high expression of BOB-1 may be one the factors driving the exquisite sensitivity of MM cells to proteasome inhibitor. Interestingly, mass spectrometry analysis revealed BOB1 in a protein complex with mTOR, Raptor and mLST8 proteins which are members of mTORC1 complex which is also involved in ribosomal function and may suggest the mechanism of action of Bob-1 at molecular level. In conclusion, here we report BOB1 as a specific proximal dependency in MM cells with potential role in modulating the protein load/capacity balance and ribosomal biogenesis essential for MM cell protein production function and therefore their sensitivity to proteasome inhibition. Disclosures Shirasaki: FIMECS: Consultancy. Mitsiades: BMS: Research Funding; Nurix: Research Funding; H3 Biomedicine: Research Funding; Novartis: Research Funding; Abbvie: Research Funding; Arch Oncology: Research Funding; Janssen/Johnson & Johnson: Research Funding; Fate Therapeutics: Consultancy, Honoraria; Karyopharm: Research Funding; Sanofi: Research Funding; TEVA: Research Funding; EMD Serono: Research Funding; Adicet Bio: Membership on an entity's Board of Directors or advisory committees; FIMECS: Consultancy, Honoraria; Ionis Pharmaceuticals: Consultancy, Honoraria. Hajek: BMS: Consultancy, Honoraria, Research Funding; AbbVie: Consultancy, Honoraria; Celgene: Consultancy, Honoraria, Research Funding; Novartis: Consultancy, Research Funding; Janssen: Consultancy, Honoraria, Research Funding; Takeda: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Amgen: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Pharma MAR: Consultancy, Honoraria. Munshi: Abbvie: Consultancy; Oncopep: Consultancy, Current equity holder in publicly-traded company, Other: scientific founder, Patents & Royalties; Pfizer: Consultancy; Legend: Consultancy; Bristol-Myers Squibb: Consultancy; Janssen: Consultancy; Karyopharm: Consultancy; Celgene: Consultancy; Adaptive Biotechnology: Consultancy; Takeda: Consultancy; Amgen: Consultancy; Novartis: Consultancy.

TET2 mutations are associated with hypermethylation at key regulatory enhancers in normal and malignant hematopoiesis

Mutations in the epigenetic modifier TET2 are frequent in myeloid malignancies and clonal hematopoiesis of indeterminate potential (CHIP) and clonal cytopenia of undetermined significance (CCUS). Here, we investigate associations between TET2 mutations and DNA methylation in whole blood in 305 elderly twins, 15 patients with CCUS and 18 healthy controls. We find that TET2 mutations are associated with DNA hypermethylation at enhancer sites in whole blood in CHIP and in both granulocytes and mononuclear cells in CCUS. These hypermethylated sites are associated with leukocyte function and immune response and ETS-related and C/EBP-related transcription factor motifs. While the majority of TET2-associated hypermethylation sites are shared between CHIP and in AML, we find a set of AML-specific hypermethylated loci at active enhancer elements in hematopoietic stem cells. In summary, we show that TET2 mutations is associated with hypermethylated enhancers involved in myeloid differentiation in both CHIP, CCUS and AML patients.

High Enhancer Activity is an Epigenetic Feature of HPV Negative Atypical Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous disease with significant morbidity and mortality and frequent recurrence. Pre-NGS efforts to transcriptionally classify HNSCC into groups of varying prognosis have identified four accepted molecular subtypes of disease: Atypical (AT), Basal (BA), Classical (CL), and Mesenchymal (MS). Here, we investigated the active enhancer landscapes of these subtypes using representative HNSCC cell lines and identified samples belonging to the AT subtype as having increased enhancer activity compared to the other 3 HNSCC subtypes. Cell lines belonging to atypical subtype were more resistant to bromodomain inhibitors (BETi). PRO-Seq experiments that both TCGA tumors and AT cell lines showed higher eRNA transcripts for enhancers controlling BETi resistance pathways, such as lipid metabolism and MAPK signaling. Additionally, HiChIP experiments suggested higher enhancer-promoter (E-P) contacts in the AT subtype, including on genes identified in the eRNA analysis. Consistently, known BETi resistance pathways were upregulated upon exposure to these inhibitors. Together, our results identify that the AT subtype of HNSCC is associated with high enhancer activity, resistance to BET inhibition, and signaling pathways that could serve as future targets for sensitizing HNSCC to BET inhibition.

Liver enhancer signature and regulatory network of non-alcoholic steatohepatitis resistance in humans with obesity

Non-alcoholic fatty liver disease (NAFLD), which often co-occurs with obesity and type 2 diabetes, is the most common cause of chronic liver disease worldwide. A subset of patients with NAFLD progress to the severe form known as non-alcoholic steatohepatitis (NASH), which increases the risk of developing hepatic fibrosis, cirrhosis, and hepatocellular carcinoma. The molecular underpinnings of the progression from NAFLD to NASH in patients is poorly understood. Active enhancer landscapes are known to determine cell states and can point to key transcription factors (TF), genes, and pathways in disease pathogenesis or resistance. Also, while super-enhancers have helped reveal key disease drivers in several cancer types, they remain undefined in NASH. To define the enhancer signature of NASH-prone (NP) and NASH-resistant (NR) phenotypes in humans, we performed chromatin run-on sequencing (ChRO-seq) analysis on liver biopsies of individuals with severe obesity who were stratified into either the NP or NR group. We defined the active enhancer signature, super-enhancer linked genes, and the candidate TF networks in human NP and NR livers. Notably, we showed that NR-activated genes that are linked to NR-specific enhancers/super-enhancers are involved in serine and glycine biosynthesis (SHMT, BHMT) as well as glycine utilization (GLYAT, GATM). Overall, this study has defined for the first time the active enhancer/super-enhancer landscape and the underlying TF programs of NASH resistance in humans with obesity.

Oncogenic enhancers drive esophageal squamous cell carcinogenesis and metastasis

The role of cis-elements and their aberrations remains unclear in esophageal squamous cell carcinoma (ESCC, further abbreviated EC). Here we survey 28 H3K27ac-marked active enhancer profiles and 50 transcriptomes in primary EC, metastatic lymph node cancer (LNC), and adjacent normal (Nor) esophageal tissues. Thousands of gained or lost enhancers and hundreds of altered putative super-enhancers are identified in EC and LNC samples respectively relative to Nor, with a large number of common gained or lost enhancers. Moreover, these differential enhancers contribute to the transcriptomic aberrations in ECs and LNCs. We also reveal putative driver onco-transcription factors, depletion of which diminishes cell proliferation and migration. The administration of chemical inhibitors to suppress the predicted targets of gained super-enhances reveals HSP90AA1 and PDE4B as potential therapeutic targets for ESCC. Thus, our epigenomic profiling reveals a compendium of reprogrammed cis-regulatory elements during ESCC carcinogenesis and metastasis for uncovering promising targets for cancer treatment.

LncRNA AC007255.1, an immune-related prognostic enhancer RNA in esophageal cancer

Background Growing evidence has suggested that enhancer RNAs (eRNAs), a set of long non-coding RNAs (lncRNAs) that were derived from active enhancer regions, play critical roles in regulating gene expression in human cancers. Nevertheless potential functions of eRNAs in esophageal cancer ESCA have not yet been expounded. Here, this study aimed to explore key prognostic eRNAs in ESCA. Methods LncRNAs that were transcribed from active enhancer regions were analyzed utilizing the PreSTIGE algorithm, followed by prediction of their target genes. Based on the ESCA RNA-seq data from the TANRIC database, overall survival (OS)-related eRNAs were determined. The correlation between AC007255.1 expression and various clinical traits of ESCA was calculated. Functional enrichment analysis was presented based on its co-expressed genes. Based on the TIMER database, we analyzed correlations between AC007255.1 expression and immune infiltration levels. qRT-PCR was utilized to validate the expression of AC007255.1 and PRR15 in ESCA and normal tissues. Results Totally, 2,695 lncRNAs were transcribed from active enhancer regions. Among them, 33 were significantly related to OS. AC007255.1 was a key eRNA. PRR15 was a target gene of AC007255.1 (correlation coefficient r = 0.936). Patients with high AC007255.1 expression indicated poor OS time. There were significant correlations between AC007255.1 expression and clinical characteristics like pathological TNM, grade and stage. AC007255.1 was closely related to tight junction and neutrophil activation involved in immune response. Moreover, AC007255.1 expression was related to the infiltration levels of B cell, dendritic cell and neutrophil. qRT-PCR results confirmed that AC007255.1 and PRR15 were both up-regulated in ESCA tissues, and there was a positive correlation between the two. Conclusion Our findings identified a novel immune-related eRNA AC007255.1 in ESCA, which could be a promising prognostic factor for ESCA.

The MLL3/4 H3K4 methyltransferase complex in establishing an active enhancer landscape

Enhancers are cis-regulatory elements that play essential roles in tissue-specific gene expression during development. Enhancer function in the expression of developmental genes requires precise regulation, while deregulation of enhancer function could be the main cause of tissue-specific cancer development. MLL3/KMT2C and MLL4/KMT2D are two paralogous histone modifiers that belong to the SET1/MLL (also named COMPASS) family of lysine methyltransferases and play critical roles in enhancer-regulated gene activation. Importantly, large-scale DNA sequencing studies have revealed that they are amongst the most frequently mutated genes associated with human cancers. MLL3 and MLL4 form identical multi-protein complexes for modifying mono-methylation of histone H3 lysine 4 (H3K4) at enhancers, which together with the p300/CBP-mediated H3K27 acetylation can generate an active enhancer landscape for long-range target gene activation. Recent studies have provided a better understanding of the possible mechanisms underlying the roles of MLL3/MLL4 complexes in enhancer regulation. Moreover, accumulating studies offer new insights into our knowledge of the potential role of MLL3/MLL4 in cancer development. In this review, we summarize recent evidence on the molecular mechanisms of MLL3/MLL4 in the regulation of active enhancer landscape and long-range gene expression, and discuss their clinical implications in human cancers.

Functional in vivo characterization of sox10 enhancers in neural crest and melanoma development

The role of a neural crest developmental transcriptional program, which critically involves Sox10 upregulation, is a key conserved aspect of melanoma initiation in both humans and zebrafish, yet transcriptional regulation of sox10 expression is incompletely understood. Here we used ATAC-Seq analysis of multiple zebrafish melanoma tumors to identify recurrently open chromatin domains as putative melanoma-specific sox10 enhancers. Screening in vivo with EGFP reporter constructs revealed 9 of 11 putative sox10 enhancers with embryonic activity in zebrafish. Focusing on the most active enhancer region in melanoma, we identified a region 23 kilobases upstream of sox10, termed peak5, that drives EGFP reporter expression in a subset of neural crest cells, Kolmer-Agduhr neurons, and early melanoma patches and tumors with high specificity. A ~200 base pair region, conserved in Cyprinidae, within peak5 is required for transgenic reporter activity in neural crest and melanoma. This region contains dimeric SoxE/Sox10 dimeric binding sites essential for peak5 neural crest and melanoma activity. We show that deletion of the endogenous peak5 conserved genomic locus decreases embryonic sox10 expression and disrupts adult stripe patterning in our melanoma model background. Our work demonstrates the power of linking developmental and cancer models to better understand neural crest identity in melanoma.

Systematic comparison of experimental assays and analytical pipelines for identification of active enhancers genome-wide

Mounting evidence supports the idea that transcriptional patterns serve as more specific identifiers of active enhancers than histone marks; however, the optimal strategy to identify active enhancers both experimentally and computationally has not been determined. In this study, we compared 13 genome-wide RNA sequencing assays in K562 cells and showed that the nuclear run-on followed by cap-selection assay (namely, GRO/PRO-cap) has significant advantages in eRNA detection and active enhancer identification. We also introduced a new analytical tool, Peak Identifier for Nascent-Transcript Sequencing (PINTS), to identify active promoters and enhancers genome-wide and pinpoint the precise location of the 5′ transcription start sites (TSSs) within these regulatory elements. Finally, we compiled a comprehensive enhancer candidate compendium based on the detected eRNA TSSs available in 120 cell and tissue types. To facilitate the exploration and prioritization of these enhancer candidates, we also built a user-friendly web server (https://pints.yulab.org) for the compendium with various additional genomic and epigenomic annotations. With the knowledge of the best available assays and pipelines, this large-scale annotation of candidate enhancers will pave the road for selection and characterization of their functions in a time-, labor-, and cost-effective manner in future.