VHL Deficiency Drives Enhancer Activation of Oncogenes in Clear Cell Renal Cell Carcinoma

The roles of nuclear orphan receptor NR2F6 in anti-viral innate immunity

Proper transcription regulation by key transcription factors, such as IRF3, is critical for anti-viral defense. Dynamics of enhancer activity play important roles in many biological processes, and epigenomic analysis is used to determine the involved enhancers and transcription factors. To determine new transcription factors in anti-DNA-virus response, we have performed H3K27ac ChIP-Seq and identified three transcription factors, NR2F6, MEF2D and MAFF, in promoting HSV-1 replication. NR2F6 promotes HSV-1 replication and gene expression in vitro and in vivo, but not dependent on cGAS/STING pathway. NR2F6 binds to the promoter of MAP3K5 and activates AP-1/c-Jun pathway, which is critical for DNA virus replication. On the other hand, NR2F6 is transcriptionally repressed by c-Jun and forms a negative feedback loop. Meanwhile, cGAS/STING innate immunity signaling represses NR2F6 through STAT3. Taken together, we have identified new transcription factors and revealed the underlying mechanisms involved in the network between DNA viruses and host cells.

Integrative multiomics enhancer activity profiling identifies therapeutic vulnerabilities in cholangiocarcinoma of different etiologies

OBJECTIVES

Cholangiocarcinoma (CCA) is a heterogeneous malignancy with high mortality and dismal prognosis, and an urgent clinical need for new therapies. Knowledge of the CCA epigenome is largely limited to aberrant DNA methylation. Dysregulation of enhancer activities has been identified to affect carcinogenesis and leveraged for new therapies but is uninvestigated in CCA. Our aim is to identify potential therapeutic targets in different subtypes of CCA through enhancer profiling.

DESIGN

Integrative multiomics enhancer activity profiling of diverse CCA was performed. A panel of diverse CCA cell lines, patient-derived and cell line-derived xenografts were used to study identified enriched pathways and vulnerabilities. NanoString, multiplex immunohistochemistry staining and single-cell spatial transcriptomics were used to explore the immunogenicity of diverse CCA.

RESULTS

We identified three distinct groups, associated with different etiologies and unique pathways. Drug inhibitors of identified pathways reduced tumour growth in in vitro and in vivo models. The first group (ESTRO), with mostly fluke-positive CCAs, displayed activation in estrogen signalling and were sensitive to MTOR inhibitors. Another group (OXPHO), with mostly BAP1 and IDH-mutant CCAs, displayed activated oxidative phosphorylation pathways, and were sensitive to oxidative phosphorylation inhibitors. Immune-related pathways were activated in the final group (IMMUN), made up of an immunogenic CCA subtype and CCA with aristolochic acid (AA) mutational signatures. Intratumour differences in AA mutation load were correlated to intratumour variation of different immune cell populations.

CONCLUSION

Our study elucidates the mechanisms underlying enhancer dysregulation and deepens understanding of different tumourigenesis processes in distinct CCA subtypes, with potential significant therapeutics and clinical benefits.

Histone 3 lysine 9 acetylation-specific reprogramming regulates esophageal squamous cell carcinoma progression and metastasis

Dysregulation of histone acetylation is widely implicated in tumorigenesis, yet its specific roles in the progression and metastasis of esophageal squamous cell carcinoma (ESCC) remain unclear. Here, we profiled the genome-wide landscapes of H3K9ac for paired adjacent normal (Nor), primary ESCC (EC) and metastatic lymph node (LNC) esophageal tissues from three ESCC patients. Compared to H3K27ac, we identified a distinct epigenetic reprogramming specific to H3K9ac in EC and LNC samples relative to Nor samples. This H3K9ac-related reprogramming contributed to the transcriptomic aberration of targeting genes, which were functionally associated with tumorigenesis and metastasis. Notably, genes with gained H3K9ac signals in both primary and metastatic lymph node samples (common-gained gene) were significantly enriched in oncogenes. Single-cell RNA-seq analysis further revealed that the corresponding top 15 common-gained genes preferred to be enriched in mesenchymal cells with high metastatic potential. Additionally, in vitro experiment demonstrated that the removal of H3K9ac from the common-gained gene MSI1 significantly downregulated its transcription, resulting in deficiencies in ESCC cell proliferation and migration. Together, our findings revealed the distinct characteristics of H3K9ac in esophageal squamous cell carcinogenesis and metastasis, and highlighted the potential therapeutic avenue for intervening ESCC through epigenetic modulation via H3K9ac.

VHL mutation drives human clear cell renal cell carcinoma progression through PI3K/AKT-dependent cholesteryl ester accumulation

BACKGROUND

Cholesteryl ester (CE) accumulation in intracellular lipid droplets (LDs) is an essential signature of clear cell renal cell carcinoma (ccRCC), but its molecular mechanism and pathological significance remain elusive.

METHODS

Enabled by the label-free Raman spectromicroscopy, which integrated stimulated Raman scattering microscopy with confocal Raman spectroscopy on the same platform, we quantitatively analyzed LD distribution and composition at the single cell level in intact ccRCC cell and tissue specimens in situ without any processing or exogenous labeling. Since we found that commonly used ccRCC cell lines actually did not show the CE-rich signature, primary cancer cells were isolated from human tissues to retain the lipid signature of ccRCC with CE level as high as the original tissue, which offers a preferable cell model for the study of cholesterol metabolism in ccRCC. Moreover, we established a patient-derived xenograft (PDX) mouse model that retained the CE-rich phenotype of human ccRCC.

FINDINGS

Surprisingly, our results revealed that CE accumulation was induced by tumor suppressor VHL mutation, the most common mutation of ccRCC. Moreover, VHL mutation was found to promote CE accumulation by upregulating HIFα and subsequent PI3K/AKT/mTOR/SREBPs pathway. Inspiringly, inhibition of cholesterol esterification remarkably suppressed ccRCC aggressiveness in vitro and in vivo with negligible toxicity, through the reduced membrane cholesterol-mediated downregulations of integrin and MAPK signaling pathways.

INTERPRETATION

Collectively, our study improves current understanding of the role of CE accumulation in ccRCC and opens up new opportunities for treatment.

FUNDING

This work was supported by National Natural Science Foundation of China (No. U23B2046 and No. 62027824), National Key R&D Program of China (No. 2023YFC2415500), Fundamental Research Funds for the Central Universities (No. YWF-22-L-547), PKU-Baidu Fund (No. 2020BD033), Peking University First Hospital Scientific and Technological Achievement Transformation Incubation Guidance Fund (No. 2022CX02), and Beijing Municipal Health Commission (No. 2020-2Z-40713).

Integrative epigenome-transcriptome analysis unravels cancer-specific over-expressed genes potentially regulating immune microenvironment in clear cell renal cell carcinoma

BACKGROUND

Clear cell Renal Cell Carcinoma (ccRCC) is the frequently diagnosed histological life-threatening tumor subtype in the urinary system. Integrating multi-omics data is emerging as a tool to provide a comprehensive view of biology and disease for better therapeutic interventions.

METHOD

We have integrated freely available ccRCC data sets of genome-wide DNA methylome, transcriptome, and active histone modification marks, H3K27ac, H3K4me1, and H3K4me3 specific ChIP-seq data to screen genes with higher expression. Further, these genes were filtered based on their effect on survival upon alteration in expression.

RESULTS

The six multi-omics-based identified genes, RUNX1, MSC, ADA, TREML1, TGFA, and VWF, showed higher expression with enrichment of active histone marks and hypomethylated CpG in ccRCC. In continuation, the identified genes were validated by an independent dataset and showed a correlation with nodal and metastatic status. Furthermore, gene ontology and pathway analysis revealed that immune-related pathways are activated in ccRCC patients.

CONCLUSIONS

The network analysis of six overexpressed genes suggests their potential role in an immunosuppressive environment, leading to tumor progression and poor prognosis. Our study shows that the multi-omics approach helps unravel complex biology for patient subtyping and proposes combination strategies with epi-drugs for more precise immunotherapy in ccRCC.

Clinical Significance and Expression Pattern of RIP5 and VGLL4 in Clear Cell Renal Cell Carcinoma Patients Treated with Sunitinib

While clear cell renal cell carcinoma (ccRCC) is curable, advanced metastatic (mRCC) remains a clinical challenge. We analyzed clinical, pathohistological, and molecular data (Receptor Interacting Protein 5-RIP5 and Vestigial Like Family Member 4-VGLL4 expression) of 55 mRCC patients treated with first-line treatment with sunitinib. The trend of linear increase in the protein expression of RIP5 was observed with the progression of tumor grade. Overall, 80% of RIP5-positive cells were in the control kidneys and high-grade mRCC. On the contrary, RIP5 displayed low expression in grade 2 mRCC (5.63%). The trend of linear decrease in the expression of VGLL4 was observed with the progression of tumor grade. The highest protein expression of VGLL4 was observed in grade 2 (87.82%) in comparison to grade 3 and 4 and control. High expression of RIP5 mRNA was associated with longer first-line overall survival and longer progression-free survival in mRCC. In addition, a high VGLL4 mRNA expression showed better overall survival in patients with ccRCC. In conclusion, high mRNA expression of RIP5 and VGLL4 are important markers of better survival rates in mRCC patients.

DCAF1-based PROTACs with activity against clinically validated targets overcoming intrinsic- and acquired-degrader resistance

Targeted protein degradation (TPD) mediates protein level through small molecule induced redirection of E3 ligases to ubiquitinate neo-substrates and mark them for proteasomal degradation. TPD has recently emerged as a key modality in drug discovery. So far only a few ligases have been utilized for TPD. Interestingly, the workhorse ligase CRBN has been observed to be downregulated in settings of resistance to immunomodulatory inhibitory drugs (IMiDs). Here we show that the essential E3 ligase receptor DCAF1 can be harnessed for TPD utilizing a selective, non-covalent DCAF1 binder. We confirm that this binder can be functionalized into an efficient DCAF1-BRD9 PROTAC. Chemical and genetic rescue experiments validate specific degradation via the CRL4DCAF1 E3 ligase. Additionally, a dasatinib-based DCAF1 PROTAC successfully degrades cytosolic and membrane-bound tyrosine kinases. A potent and selective DCAF1-BTK-PROTAC (DBt-10) degrades BTK in cells with acquired resistance to CRBN-BTK-PROTACs while the DCAF1-BRD9 PROTAC (DBr-1) provides an alternative strategy to tackle intrinsic resistance to VHL-degrader, highlighting DCAF1-PROTACS as a promising strategy to overcome ligase mediated resistance in clinical settings.

Integrating scRNA and bulk-RNA sequencing develops a cell senescence signature for analyzing tumor heterogeneity in clear cell renal cell carcinoma

Introduction

Cellular senescence (CS) plays a critical role in cancer development, including clear cell renal cell carcinoma (ccRCC). Traditional RNA sequencing cannot detect precise molecular composition changes within tumors. This study aimed to analyze cellular senescence's biochemical characteristics in ccRCC using single RNA sequencing (ScRNA-seq) and traditional RNA sequencing (Bulk RNA-seq).

Methods

Researchers analyzed the biochemical characteristics of cellular senescence in ccRCC using ScRNA-seq and Bulk RNA-seq. They combined these approaches to identify differences between malignant and non-malignant phenotypes in ccRCC across three senescence-related pathways. Genes from these pathways were used to identify molecular subtypes associated with senescence, and a new risk model was constructed. The function of the gene DUSP1 in ccRCC was validated through biological experiments.

Results

The combined analysis of ScRNA-seq and Bulk RNA-seq revealed significant differences between malignant and non-malignant phenotypes in ccRCC across three senescence-related pathways. Researchers identified genes from these pathways to identify molecular subtypes associated with senescence, constructing a new risk model. Different subgroups showed significant differences in prognosis level, clinical stage and grade, immune infiltration, immunotherapy, and drug sensitivity.

Discussion

Senescence signature markers are practical biomarkers and predictors of molecular typing in ccRCC. Differences in prognosis level, clinical stage and grade, immune infiltration, immunotherapy, and drug sensitivity between different subgroups indicate that this approach could provide valuable insights into senescence-related treatment options and prognostic assessment for patients with ccRCC. The function of the gene DUSP1 in ccRCC was validated through biological experiments, confirming its feasibility as a novel biomarker for ccRCC. These findings suggest that targeted therapies based on senescence-related mechanisms could be an effective treatment option for ccRCC.

Unravelling the Role of P300 and TMPRSS2 in Prostate Cancer: A Literature Review

Prostate cancer is one of the most common malignant diseases in men, and it contributes significantly to the increased mortality rate in men worldwide. This study aimed to review the roles of p300 and TMPRSS2 (transmembrane protease, serine 2) in the AR (androgen receptor) pathway as they are closely related to the development and progression of prostate cancer. This paper represents a library-based study conducted by selecting the most suitable, up-to-date scientific published articles from online journals. We focused on articles that use similar techniques, particularly those that use prostate cancer cell lines and immunohistochemical staining to study the molecular impact of p300 and TMPRSS2 in prostate cancer specimens. The TMPRSS2:ERG fusion is considered relevant to prostate cancer, but its association with the development and progression as well as its clinical significance have not been fully elucidated. On the other hand, high p300 levels in prostate cancer biopsies predict larger tumor volumes, extraprostatic extension of disease, and seminal vesicle involvement at prostatectomy, and may be associated with prostate cancer progression after surgery. The inhibition of p300 has been shown to reduce the proliferation of prostate cancer cells with TMPRSS2:ETS (E26 transformation-specific) fusions, and combining p300 inhibitors with other targeted therapies may increase their efficacy. Overall, the interplay between the p300 and TMPRSS2 pathways is an active area of research.

Etiology of super-enhancer reprogramming and activation in cancer

Super-enhancers are large, densely concentrated swaths of enhancers that regulate genes critical for cell identity. Tumorigenesis is accompanied by changes in the super-enhancer landscape. These aberrant super-enhancers commonly form to activate proto-oncogenes, or other genes upon which cancer cells depend, that initiate tumorigenesis, promote tumor proliferation, and increase the fitness of cancer cells to survive in the tumor microenvironment. These include well-recognized master regulators of proliferation in the setting of cancer, such as the transcription factor MYC which is under the control of numerous super-enhancers gained in cancer compared to normal tissues. This Review will cover the expanding cell-intrinsic and cell-extrinsic etiology of these super-enhancer changes in cancer, including somatic mutations, copy number variation, fusion events, extrachromosomal DNA, and 3D chromatin architecture, as well as those activated by inflammation, extra-cellular signaling, and the tumor microenvironment.

mTOR Regulation of N-Myc Downstream Regulated 1 (NDRG1) Phosphorylation in Clear Cell Renal Cell Carcinoma

The mechanistic target of rapamycin (mTOR) kinase is a component of two signaling complexes that are known as mTOR complex 1 (mTORC1) and mTORC2. We sought to identify mTOR-phosphorylated proteins that are differently expressed in clinically resected clear cell renal cell carcinoma (ccRCC) relative to pair-matched normal renal tissue. Using a proteomic array, we found N-Myc Downstream Regulated 1 (NDRG1) showed the greatest increase (3.3-fold) in phosphorylation (on Thr346) in ccRCC. This was associated with an increase in total NDRG1. RICTOR is a required subunit in mTORC2, and its knockdown decreased total and phospho-NDRG1 (Thr346) but not NDRG1 mRNA. The dual mTORC1/2 inhibitor, Torin 2, significantly reduced (by ~100%) phospho-NDRG1 (Thr346). Rapamycin is a selective mTORC1 inhibitor that had no effect on the levels of total NDRG1 or phospho-NDRG1 (Thr346). The reduction in phospho-NDRG1 (Thr346) due to the inhibition of mTORC2 corresponded with a decrease in the percentage of live cells, which was correlated with an increase in apoptosis. Rapamycin had no effect on ccRCC cell viability. Collectively, these data show that mTORC2 mediates the phosphorylation of NDRG1 (Thr346) in ccRCC. We hypothesize that RICTOR and mTORC2-mediated phosphorylation of NDRG1 (Thr346) promotes the viability of ccRCC cells.

Hypoxia induces polycystin-1 expression in the renal epithelium

Mutations in polycystin-1 which is encoded by the PKD1 gene are the main causes for the development of autosomal dominant polycystic kidney disease. However, only little is known about the physiological function of polycystin-1 and even less about the regulation of its expression. Here, we show that expression of PKD1 is induced by hypoxia and compounds that stabilize the hypoxia-inducible transcription factor (HIF) 1α in primary human tubular epithelial cells. Knockdown of HIF subunits confirms HIF-1α-dependent regulation of polycystin-1 expression. Furthermore, HIF ChIP-seq reveals that HIF interacts with a regulatory DNA element within the PKD1 gene in renal tubule-derived cells. HIF-dependent expression of polycystin-1 can also be demonstrated in vivo in kidneys of mice treated with substances that stabilize HIF. Polycystin-1 and HIF-1α have been shown to promote epithelial branching during kidney development. In line with these findings, we show that expression of polycystin-1 within mouse embryonic ureteric bud branches is regulated by HIF. Our finding links expression of one of the main regulators of accurate renal development with the hypoxia signalling pathway and provides additional insight into the pathophysiology of polycystic kidney disease.

PBRM1-deficient PBAF complexes target aberrant genomic loci to activate the NF-κB pathway in clear cell renal cell carcinoma

PBRM1 encodes an accessory subunit of the PBAF SWI/SNF chromatin remodeller, and the inactivation of PBRM1 is a frequent event in kidney cancer. However, the impact of PBRM1 loss on chromatin remodelling is not well examined. Here we show that, in VHL-deficient renal tumours, PBRM1 deficiency results in ectopic PBAF complexes that localize to de novo genomic loci, activating the pro-tumourigenic NF-κB pathway. PBRM1-deficient PBAF complexes retain the association between SMARCA4 and ARID2, but have loosely tethered BRD7. The PBAF complexes redistribute from promoter proximal regions to distal enhancers containing NF-κB motifs, heightening NF-κB activity in PBRM1-deficient models and clinical samples. The ATPase function of SMARCA4 maintains chromatin occupancy of pre-existing and newly acquired RELA specific to PBRM1 loss, activating downstream target gene expression. Proteasome inhibitor bortezomib abrogates RELA occupancy, suppresses NF-κB activation and delays growth of PBRM1-deficient tumours. In conclusion, PBRM1 safeguards the chromatin by repressing aberrant liberation of pro-tumourigenic NF-κB target genes by residual PBRM1-deficient PBAF complexes.

Inflammatory Networks in Renal Cell Carcinoma

Cancer-associated inflammation has been established as a hallmark feature of almost all solid cancers. Tumor-extrinsic and intrinsic signaling pathways regulate the process of cancer-associated inflammation. Tumor-extrinsic inflammation is triggered by many factors, including infection, obesity, autoimmune disorders, and exposure to toxic and radioactive substances. Intrinsic inflammation can be induced by genomic mutation, genome instability and epigenetic remodeling in cancer cells that promote immunosuppressive traits, inducing the recruitment and activation of inflammatory immune cells. In RCC, many cancer cell-intrinsic alterations are assembled, upregulating inflammatory pathways, which enhance chemokine release and neoantigen expression. Furthermore, immune cells activate the endothelium and induce metabolic shifts, thereby amplifying both the paracrine and autocrine inflammatory loops to promote RCC tumor growth and progression. Together with tumor-extrinsic inflammatory factors, tumor-intrinsic signaling pathways trigger a Janus-faced tumor microenvironment, thereby simultaneously promoting or inhibiting tumor growth. For therapeutic success, it is important to understand the pathomechanisms of cancer-associated inflammation, which promote cancer progression. In this review, we describe the molecular mechanisms of cancer-associated inflammation that influence cancer and immune cell functions, thereby increasing tumor malignancy and anti-cancer resistance. We also discuss the potential of anti-inflammatory treatments, which may provide clinical benefits in RCCs and possible avenues for therapy and future research.

Epigenetic and transcriptomic characterization reveals progression markers and essential pathways in clear cell renal cell carcinoma

Identifying tumor-cell-specific markers and elucidating their epigenetic regulation and spatial heterogeneity provides mechanistic insights into cancer etiology. Here, we perform snRNA-seq and snATAC-seq in 34 and 28 human clear cell renal cell carcinoma (ccRCC) specimens, respectively, with matched bulk proteogenomics data. By identifying 20 tumor-specific markers through a multi-omics tiered approach, we reveal an association between higher ceruloplasmin (CP) expression and reduced survival. CP knockdown, combined with spatial transcriptomics, suggests a role for CP in regulating hyalinized stroma and tumor-stroma interactions in ccRCC. Intratumoral heterogeneity analysis portrays tumor cell-intrinsic inflammation and epithelial-mesenchymal transition (EMT) as two distinguishing features of tumor subpopulations. Finally, BAP1 mutations are associated with widespread reduction of chromatin accessibility, while PBRM1 mutations generally increase accessibility, with the former affecting five times more accessible peaks than the latter. These integrated analyses reveal the cellular architecture of ccRCC, providing insights into key markers and pathways in ccRCC tumorigenesis.

Aberrant JAK-STAT signaling-mediated chromatin remodeling impairs the sensitivity of NK/T-cell lymphoma to chidamide

BACKGROUND

Natural killer/T-cell lymphoma (NKTL) is a rare type of aggressive and heterogeneous non-Hodgkin's lymphoma (NHL) with a poor prognosis and limited therapeutic options. Therefore, there is an urgent need to exploit potential novel therapeutic targets for the treatment of NKTL. Histone deacetylase (HDAC) inhibitor chidamide was recently approved for treating relapsed/refractory peripheral T-cell lymphoma (PTCL) patients. However, its therapeutic efficacy in NKTL remains unclear.

METHODS

We performed a phase II clinical trial to evaluate the efficacy of chidamide in 28 relapsed/refractory NKTL patients. Integrative transcriptomic, chromatin profiling analysis and functional studies were performed to identify potential predictive biomarkers and unravel the mechanisms of resistance to chidamide. Immunohistochemistry (IHC) was used to validate the predictive biomarkers in tumors from the clinical trial.

RESULTS

We demonstrated that chidamide is effective in treating relapsed/refractory NKTL patients, achieving an overall response and complete response rate of 39 and 18%, respectively. In vitro studies showed that hyperactivity of JAK-STAT signaling in NKTL cell lines was associated with the resistance to chidamide. Mechanistically, our results revealed that aberrant JAK-STAT signaling remodels the chromatin and confers resistance to chidamide. Subsequently, inhibition of JAK-STAT activity could overcome resistance to chidamide by reprogramming the chromatin from a resistant to sensitive state, leading to synergistic anti-tumor effect in vitro and in vivo. More importantly, our clinical data demonstrated that combinatorial therapy with chidamide and JAK inhibitor ruxolitinib is effective against chidamide-resistant NKTL. In addition, we identified TNFRSF8 (CD30), a downstream target of the JAK-STAT pathway, as a potential biomarker that could predict NKTL sensitivity to chidamide.

CONCLUSIONS

Our study suggests that chidamide, in combination with JAK-STAT inhibitors, can be a novel targeted therapy in the standard of care for NKTL.

TRIAL REGISTRATION

ClinicalTrials.gov, NCT02878278. Registered 25 August 2016, https://clinicaltrials.gov/ct2/show/NCT02878278.

Regulatory enhancer profiling of mesenchymal-type gastric cancer reveals subtype-specific epigenomic landscapes and targetable vulnerabilities

OBJECTIVE

Gastric cancer (GC) comprises multiple molecular subtypes. Recent studies have highlighted mesenchymal-subtype GC (Mes-GC) as a clinically aggressive subtype with few treatment options. Combining multiple studies, we derived and applied a consensus Mes-GC classifier to define the Mes-GC enhancer landscape revealing disease vulnerabilities.

DESIGN

Transcriptomic profiles of ~1000 primary GCs and cell lines were analysed to derive a consensus Mes-GC classifier. Clinical and genomic associations were performed across >1200 patients with GC. Genome-wide epigenomic profiles (H3K27ac, H3K4me1 and assay for transposase-accessible chromatin with sequencing (ATAC-seq)) of 49 primary GCs and GC cell lines were generated to identify Mes-GC-specific enhancer landscapes. Upstream regulators and downstream targets of Mes-GC enhancers were interrogated using chromatin immunoprecipitation followed by sequencing (ChIP-seq), RNA sequencing, CRISPR/Cas9 editing, functional assays and pharmacological inhibition.

RESULTS

We identified and validated a 993-gene cancer-cell intrinsic Mes-GC classifier applicable to retrospective cohorts or prospective single samples. Multicohort analysis of Mes-GCs confirmed associations with poor patient survival, therapy resistance and few targetable genomic alterations. Analysis of enhancer profiles revealed a distinctive Mes-GC epigenomic landscape, with TEAD1 as a master regulator of Mes-GC enhancers and Mes-GCs exhibiting preferential sensitivity to TEAD1 pharmacological inhibition. Analysis of Mes-GC super-enhancers also highlighted NUAK1 kinase as a downstream target, with synergistic effects observed between NUAK1 inhibition and cisplatin treatment.

CONCLUSION

Our results establish a consensus Mes-GC classifier applicable to multiple transcriptomic scenarios. Mes-GCs exhibit a distinct epigenomic landscape, and TEAD1 inhibition and combinatorial NUAK1 inhibition/cisplatin may represent potential targetable options.

Epigenomic charting and functional annotation of risk loci in renal cell carcinoma

While the mutational and transcriptional landscapes of renal cell carcinoma (RCC) are well-known, the epigenome is poorly understood. We characterize the epigenome of clear cell (ccRCC), papillary (pRCC), and chromophobe RCC (chRCC) by using ChIP-seq, ATAC-Seq, RNA-seq, and SNP arrays. We integrate 153 individual data sets from 42 patients and nominate 50 histology-specific master transcription factors (MTF) to define RCC histologic subtypes, including EPAS1 and ETS-1 in ccRCC, HNF1B in pRCC, and FOXI1 in chRCC. We confirm histology-specific MTFs via immunohistochemistry including a ccRCC-specific TF, BHLHE41. FOXI1 overexpression with knock-down of EPAS1 in the 786-O ccRCC cell line induces transcriptional upregulation of chRCC-specific genes, TFCP2L1, ATP6V0D2, KIT, and INSRR, implicating FOXI1 as a MTF for chRCC. Integrating RCC GWAS risk SNPs with H3K27ac ChIP-seq and ATAC-seq data reveals that risk-variants are significantly enriched in allelically-imbalanced peaks. This epigenomic atlas in primary human samples provides a resource for future investigation.

VHL Ser65 mutations enhance HIF2α signaling and promote epithelial-mesenchymal transition of renal cancer cells

Background

Von Hippel-Lindau (VHL) disease is an autosomal dominant genetic neoplastic disorder caused by germline mutation or deletion of the VHL gene, characterized by the tendency to develop multisystem benign or malignant tumors. The mechanism of VHL mutants in pathogenicity is poorly understand.

Results

Here we identified heterozygous missense mutations c.193T > C and c.194C > G in VHL in several patients from two Chinese families. These mutations are predicted to cause Serine (c.193T > C) to Proline and Tryptophan (c.194C > G) substitution at residue 65 of VHL protein (p.Ser65Pro and Ser65Trp). Ser65 residue, located within the β-domain and nearby the interaction sites with hypoxia-inducing factor α (HIFα), is highly conserved among different species. We observed gain of functions in VHL mutations, thereby stabilizing HIF2α protein and reprograming HIF2α genome-wide target gene transcriptional programs. Further analysis of independent cohorts of patients with renal carcinoma revealed specific HIF2α gene expression signatures in the context of VHL Ser65Pro or Ser65Trp mutation, showing high correlations with hypoxia and epithelial-mesenchymal transition signaling activities and strong associations with poor prognosis.

Conclusions

Together, our findings highlight the crucial role of pVHL-HIF dysregulation in VHL disease and strengthen the clinical relevance and significance of the missense mutations of Ser65 residue in pVHL in the familial VHL disease.

A local tumor microenvironment acquired super-enhancer induces an oncogenic driver in colorectal carcinoma

Tumors exhibit enhancer reprogramming compared to normal tissue. The etiology is largely attributed to cell-intrinsic genomic alterations. Here, using freshly resected primary CRC tumors and patient-matched adjacent normal colon, we find divergent epigenetic landscapes between CRC tumors and cell lines. Intriguingly, this phenomenon extends to highly recurrent aberrant super-enhancers gained in CRC over normal. We find one such super-enhancer activated in epithelial cancer cells due to surrounding inflammation in the tumor microenvironment. We restore this super-enhancer and its expressed gene, PDZK1IP1, following treatment with cytokines or xenotransplantation into nude mice, thus demonstrating cell-extrinsic etiology. We demonstrate mechanistically that PDZK1IP1 enhances the reductive capacity CRC cancer cells via the pentose phosphate pathway. We show this activation enables efficient growth under oxidative conditions, challenging the previous notion that PDZK1IP1 acts as a tumor suppressor in CRC. Collectively, these observations highlight the significance of epigenomic profiling on primary specimens.

Development and validation of a novel necroptosis-related score to improve the outcomes of clear cell renal cell carcinoma

Necroptosis has been indicated as a key regulator of tumor progression. However, the prognostic regulatory role of necroptosis in clear cell renal cell carcinoma (ccRCC) needs to be further investigated. In this study, necroptosis-related subtypes were identified by mining the public cohort (n = 530) obtained from The Cancer Genome Atlas. By applying Principal Component Analysis (PCA), the necroptosis-related scores (N-Score) were developed to assess the prognosis procession of ccRCC. The results were further validated by an external clinical cohort (n = 116) obtained from the First Affiliated Hospital of Wenzhou Medical University. It has been found that N-Score could precisely distinguish the prognostic outcomes of patients as an independent risk factor (Hazard ratio = 4.990, 95% confidence interval (CI) = 2.007–12.403, p < 0.001). In addition, changes in N-Score were associated with differences in tumor mutational burden as well as immune infiltration characterization. Moreover, higher N-Scores were also correlated significantly molecular drug sensitivity and stronger immune checkpoint activity. Notably, the prognosis of ccRCC could be effectively guided by combining the N-Scores and external clinical indicators. In conclusion, N-Scores could be served as a robust and effective biomarker to improve the prognosis outcomes and targeted therapy of ccRCC.

An oncogenic JMJD6-DGAT1 axis tunes the epigenetic regulation of lipid droplet formation in clear cell renal cell carcinoma

Characterized by intracellular lipid droplet accumulation, clear cell renal cell carcinoma (ccRCC) is resistant to cytotoxic chemotherapy and is a lethal disease. Through an unbiased siRNA screen of 2-oxoglutarate (2-OG)-dependent enzymes, which play a critical role in tumorigenesis, we identified Jumonji domain-containing 6 (JMJD6) as an essential gene for ccRCC tumor development. The downregulation of JMJD6 abolished ccRCC colony formation in vitro and inhibited orthotopic tumor growth in vivo. Integrated ChIP-seq and RNA-seq analyses uncovered diacylglycerol O-acyltransferase 1 (DGAT1) as a critical JMJD6 effector. Mechanistically, JMJD6 interacted with RBM39 and co-occupied DGAT1 gene promoter with H3K4me3 to induce DGAT1 expression. JMJD6 silencing reduced DGAT1, leading to decreased lipid droplet formation and tumorigenesis. The pharmacological inhibition (or depletion) of DGAT1 inhibited lipid droplet formation in vitro and ccRCC tumorigenesis in vivo. Thus, the JMJD6-DGAT1 axis represents a potential new therapeutic target for ccRCC.

Identification of epigenetic dysregulation gene markers and immune landscape in kidney renal clear cell carcinoma by comprehensive genomic analysis

Kidney cancer is one the most lethal cancers of the urinary system, but current treatments are limited and its prognosis is poor. This study focused on kidney renal clear cell carcinoma (KIRC) and analyzed the relationship between epigenetic alterations and KIRC prognosis, and explored the prognostic significance of these findings in KIRC patients. Based on multi-omics data, differentially expressed histone-modified genes were identified using the R package limma package. Gene enhancers were detected from data in the FANTOM5 database. Gene promoters were screened using the R package ChIPseeker, and the Bumphunter in the R package CHAMP was applied to screen differentially methylated regions (DMR). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and Gene Ontology (GO) functional enrichment analysis of genes was performed using the R package clusterProfiler. We identified 51 dysregulated epigenetic protein coding genes (epi-PCGs) from 872 epi-PCGs, and categorized three molecular subtypes (C1, C2, and C3) of KIRC samples with significantly different prognosis. Notably, among the three molecular subtypes, we found a markedly differential immune features in immune checkpoints, cytokines, immune signatures, and immune cell distribution. C2 subtype had significantly lower enrichment score of IFNγ, cytotoxic score (CYT), and angiogenesis. In addition, an 8-gene signature containing 8 epi-PCGs (ETV4, SH2B3, FATE1, GRK5, MALL, HRH2, SEMA3G, and SLC10A6) was developed for predicting KIRC prognosis. Prognosis of patients with a high 8-gene signature score was significantly worse than those with a low 8-gene signature score, which was also validated by the independent validation data. The 8-gene signature had a better performance compared with previous signatures of KIRC. Overall, this study highlighted the important role of epigenetic regulation in KIRC development, and explored prognostic epi-PCGs, which may provide a guidance for exploiting further pathological mechanisms of KIRC and for developing novel drug targets.

Loss of PBRM1 Alters Promoter Histone Modifications and Activates ALDH1A1 to Drive Renal Cell Carcinoma

Subunits of SWI/SNF chromatin remodeling complexes are frequently mutated in human malignancies. The PBAF complex is composed of multiple subunits, including the tumor-suppressor protein PBRM1 (BAF180), as well as ARID2 (BAF200), that are unique to this SWI/SNF complex. PBRM1 is mutated in various cancers, with a high mutation frequency in clear cell renal cell carcinoma (ccRCC). Here, we integrate RNA-seq, histone modification ChIP-seq, and ATAC-seq data to show that loss of PBRM1 results in de novo gains in H3K4me3 peaks throughout the epigenome, including activation of a retinoic acid biosynthesis and signaling gene signature. We show that one such target gene, ALDH1A1, which regulates a key step in retinoic acid biosynthesis, is consistently upregulated with PBRM1 loss in ccRCC cell lines and primary tumors, as well as non-malignant cells. We further find that ALDH1A1 increases the tumorigenic potential of ccRCC cells. Using biochemical methods, we show that ARID2 remains bound to other PBAF subunits after loss of PBRM1 and is essential for increased ALDH1A1 after loss of PBRM1, whereas other core SWI/SNF components are dispensable, including the ATPase subunit BRG1. In total, this study uses global epigenomic approaches to uncover novel mechanisms of PBRM1 tumor suppression in ccRCC.

IMPLICATIONS

This study implicates the SWI/SNF subunit and tumor-suppressor PBRM1 in the regulation of promoter histone modifications and retinoic acid biosynthesis and signaling pathways in ccRCC and functionally validates one such target gene, the aldehyde dehydrogenase ALDH1A1.

Histone acetyltransferases CBP/p300 in tumorigenesis and CBP/p300 inhibitors as promising novel anticancer agents

The histone acetyltransferases CBP and p300, often referred to as CBP/p300 due to their sequence homology and functional overlap and co-operation, are emerging as critical drivers of oncogenesis in the past several years. CBP/p300 induces histone H3 lysine 27 acetylation (H3K27ac) at target gene promoters, enhancers and super-enhancers, thereby activating gene transcription. While earlier studies indicate that CBP/p300 deletion/loss can promote tumorigenesis, CBP/p300 have more recently been shown to be over-expressed in cancer cells and drug-resistant cancer cells, activate oncogene transcription and induce cancer cell proliferation, survival, tumorigenesis, metastasis, immune evasion and drug-resistance. Small molecule CBP/p300 histone acetyltransferase inhibitors, bromodomain inhibitors, CBP/p300 and BET bromodomain dual inhibitors and p300 protein degraders have recently been discovered. The CBP/p300 inhibitors and degraders reduce H3K27ac, down-regulate oncogene transcription, induce cancer cell growth inhibition and cell death, activate immune response, overcome drug resistance and suppress tumor progression in vivo. In addition, CBP/p300 inhibitors enhance the anticancer efficacy of chemotherapy, radiotherapy and epigenetic anticancer agents, including BET bromodomain inhibitors; and the combination therapies exert substantial anticancer effects in mouse models of human cancers including drug-resistant cancers. Currently, two CBP/p300 inhibitors are under clinical evaluation in patients with advanced and drug-resistant solid tumors or hematological malignancies. In summary, CBP/p300 have recently been identified as critical tumorigenic drivers, and CBP/p300 inhibitors and protein degraders are emerging as promising novel anticancer agents for clinical translation.

The renal lineage factor PAX8 controls oncogenic signalling in kidney cancer

Large-scale human genetic data1-3 have shown that cancer mutations display strong tissue-selectivity, but how this selectivity arises remains unclear. Here, using experimental models, functional genomics and analyses of patient samples, we demonstrate that the lineage transcription factor paired box 8 (PAX8) is required for oncogenic signalling by two common genetic alterations that cause clear cell renal cell carcinoma (ccRCC) in humans: the germline variant rs7948643 at 11q13.3 and somatic inactivation of the von Hippel-Lindau tumour suppressor (VHL)4-6. VHL loss, which is observed in about 90% of ccRCCs, can lead to hypoxia-inducible factor 2α (HIF2A) stabilization6,7. We show that HIF2A is preferentially recruited to PAX8-bound transcriptional enhancers, including a pro-tumorigenic cyclin D1 (CCND1) enhancer that is controlled by PAX8 and HIF2A. The ccRCC-protective allele C at rs7948643 inhibits PAX8 binding at this enhancer and downstream activation of CCND1 expression. Co-option of a PAX8-dependent physiological programme that supports the proliferation of normal renal epithelial cells is also required for MYC expression from the ccRCC metastasis-associated amplicons at 8q21.3-q24.3 (ref. 8). These results demonstrate that transcriptional lineage factors are essential for oncogenic signalling and that they mediate tissue-specific cancer risk associated with somatic and inherited genetic variants.

Hypoxia classifier for transcriptome datasets

Molecular gene signatures are useful tools to characterize the physiological state of cell populations, but most have developed under a narrow range of conditions and cell types and are often restricted to a set of gene identities. Focusing on the transcriptional response to hypoxia, we aimed to generate widely applicable classifiers sourced from the results of a meta-analysis of 69 differential expression datasets which included 425 individual RNA-seq experiments from 33 different human cell types exposed to different degrees of hypoxia (0.1–5%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {O}_{2}$$\end{document}O2) for 2–48 h. The resulting decision trees include both gene identities and quantitative boundaries, allowing for easy classification of individual samples without control or normoxic reference. Each tree is composed of 3–5 genes mostly drawn from a small set of just 8 genes (EGLN1, MIR210HG, NDRG1, ANKRD37, TCAF2, PFKFB3, BHLHE40, and MAFF). In spite of their simplicity, these classifiers achieve over 95% accuracy in cross validation and over 80% accuracy when applied to additional challenging datasets. Our results indicate that the classifiers are able to identify hypoxic tumor samples from bulk RNAseq and hypoxic regions within tumor from spatially resolved transcriptomics datasets. Moreover, application of the classifiers to histological sections from normal tissues suggest the presence of a hypoxic gene expression pattern in the kidney cortex not observed in other normoxic organs. Finally, tree classifiers described herein outperform traditional hypoxic gene signatures when compared against a wide range of datasets. This work describes a set of hypoxic gene signatures, structured as simple decision tress, that identify hypoxic samples and regions with high accuracy and can be applied to a broad variety of gene expression datasets and formats.

MEOX2 Regulates the Growth and Survival of Glioblastoma Stem Cells by Modulating Genes of the Glycolytic Pathway and Response to Hypoxia

Simple Summary

Glioblastoma is the most common incurable primary brain tumor in adults, typically leading to death within 15 months of diagnosis. Although there is an ongoing debate in the scientific community about the precise cellular origin of this tumor, glioblastoma stem cells (GSCs), which are able to self-renew, yield a full tumor mass, and determine chemo- and radio-resistance, are recognized to have a pivotal role. Our research aims to understand the role of the mesenchyme homeobox 2 (MEOX2) transcription factor in GSCs where it is strongly and specifically expressed. We have found that MEOX2 is indeed important for the survival of these cells. In fact, when we reduce its expression in two different GSC lines, they undergo a massive death accompanied by the inhibition of key genes of the glycolytic metabolism, the main source of energy for these cells. Our results reveal a novel function for MEOX2 in glioblastoma and suggest a mechanism through which GSCs may survive even in unfavorable conditions.

Abstract

The most widely accepted hypothesis for the development of glioblastoma suggests that glioblastoma stem-like cells (GSCs) are crucially involved in tumor initiation and recurrence as well as in the occurrence of chemo- and radio-resistance. Mesenchyme homeobox 2 (MEOX2) is a transcription factor overexpressed in glioblastoma, whose expression is negatively correlated with patient survival. Starting from our observation that MEOX2 expression is strongly enhanced in six GSC lines, we performed shRNA-mediated knock-down experiments in two different GSC lines and found that MEOX2 depletion resulted in the inhibition of cell growth and sphere-forming ability and an increase in apoptotic cell death. By a deep transcriptome analysis, we identified a core group of genes modulated in response to MEOX2 knock-down. Among these genes, the repressed ones are largely enriched in genes involved in the hypoxic response and glycolytic pathway, two strictly related pathways that contribute to the resistance of high-grade gliomas to therapies. An in silico study of the regulatory regions of genes differentially expressed by MEOX2 knock-down revealed that they mainly consisted of GC-rich regions enriched for Sp1 and Klf4 binding motifs, two main regulators of metabolism in glioblastoma. Our results show, for the first time, the involvement of MEOX2 in the regulation of genes of GSC metabolism, which is essential for the survival and growth of these cells.

Nutritional Epigenetics in Cancer

Alterations in the epigenome are well known to affect cancer development and progression. Epigenetics is highly influenced by the environment, including diet, which is a source of metabolic substrates that influence the synthesis of cofactors or substrates for chromatin and RNA modifying enzymes. In addition, plants are a common source of bioactives that can directly modify the activity of these enzymes. Here, we review and discuss the impact of diet on epigenetic mechanisms, including chromatin and RNA regulation, and its potential implications for cancer prevention and treatment.

A Novel Nomogram for Prediction and Evaluation of Lymphatic Metastasis in Patients With Renal Cell Carcinoma

Background

Lymphatic metastasis is an important mechanism of renal cell carcinoma (RCC) dissemination and is an indicator of poor prognosis. Therefore, we aimed to identify predictors of lymphatic metastases (LMs) in RCC patients and to develop a new nomogram to assess the risk of LMs.

Methods

This study included patients with RCC from 2010 to 2018 in the Surveillance, Epidemiology, and Final Results (SEER) database into the training cohort and included the RCC patients diagnosed during the same period in the Second Affiliated Hospital of Dalian Medical University into the validation cohort. Univariate and multivariate logistic regression analysis were performed to identify risk factors for LM, constructing a nomogram. The receiver operating characteristic (ROC) curves were generated to assess the nomogram’s performance, and the concordance index (C-index), area under curve value (AUC), and calibration plots were used to evaluate the discrimination and calibration of the nomogram. The nomogram’s clinical performance was evaluated by decision curve analysis (DCA), probability density function (PDF) and clinical utility curve (CUC). Furthermore, Kaplan-Meier curves were performed in the training and the validation cohort to evaluate the survival risk of the patients with lymphatic metastasis or not. Additionally, on the basis of the constructed nomogram, we obtained a convenient and intuitive network calculator.

Results

A total of 41837 patients were included for analysis, including 41,018 in the training group and 819 in the validation group. Eleven risk factors were considered as predictor variables in the nomogram. The nomogram displayed excellent discrimination power, with AUC both reached 0.916 in the training group (95% confidence interval (CI) 0.913 to 0.918) and the validation group (95% CI 0.895 to 0.934). The calibration curves presented that the nomogram-based prediction had good consistency with practical application. Moreover, Kaplan-Meier curves analysis showed that RCC patients with LMs had worse survival outcomes compared with patients without LMs.

Conclusions

The nomogram and web calculator (https://liwenle0910.shinyapps.io/DynNomapp/) may be a useful tool to quantify the risk of LMs in patients with RCC, which may provide guidance for clinicians, such as identifying high-risk patients, performing surgery, and establishing personalized treatment as soon as possible.

RNA-binding protein MEX3A controls G1/S transition via regulating the RB/E2F pathway in clear cell renal cell carcinoma

MEX3A is an RNA-binding protein that mediates mRNA decay through binding to 3′ untranslated regions. However, its role and mechanism in clear cell renal cell carcinoma remain unknown. In this study, we found that MEX3A expression was transcriptionally activated by ETS1 and upregulated in clear cell renal cell carcinoma. Silencing MEX3A markedly reduced clear cell renal cell carcinoma cell proliferation in vitro and in vivo. Inhibiting MEX3A induced G1/S cell-cycle arrest. Gene set enrichment analysis revealed that E2F targets are the central downstream pathways of MEX3A. To identify MEX3A targets, systematic screening using enhanced cross-linking and immunoprecipitation sequencing, and RNA-immunoprecipitation sequencing assays were performed. A network of 4,000 genes was identified as potential targets of MEX3A. Gene ontology analysis of upregulated genes bound by MEX3A indicated that negative regulation of the cell proliferation pathway was highly enriched. Further assays indicated that MEX3A bound to the CDKN2B 3′ untranslated region, promoting its mRNA degradation. This leads to decreased levels of CDKN2B and an uncontrolled cell cycle in clear cell renal cell carcinoma, which was confirmed by rescue experiments. Our findings revealed that MEX3A acts as a post-transcriptional regulator of abnormal cell-cycle progression in clear cell renal cell carcinoma.

The renal cancer risk allele at 14q24.2 activates a novel hypoxia-inducible transcription factor-binding enhancer of DPF3 expression

Evolution of clear cell renal cell carcinoma is guided by dysregulation of hypoxia-inducible transcription factor (HIF) pathways following loss of the von Hippel-Lindau tumor suppressor protein. Renal cell carcinoma (RCC)-associated polymorphisms influence HIF–DNA interactions at enhancers of important oncogenes thereby modulating the risk of developing renal cancer. A strong signal of genome-wide association with RCC was determined for the single nucleotide polymorphism (SNP) rs4903064, located on chr14q.24.2 within an intron of DPF3, encoding for Double PHD Fingers 3, a member of chromatin remodeling complexes; however, it is unclear how the risk allele operates in renal cells. In this study, we used tissue specimens and primary renal cells from a large cohort of RCC patients to examine the function of this polymorphism. In clear cell renal cell carcinoma tissue, isolated tumor cells as well as in primary renal tubular cells, in which HIF was stabilized, we determined genotype-specific increases of DPF3 mRNA levels and identified that the risk SNP resides in an active enhancer region, creating a novel HIF-binding motif. We then confirmed allele-specific HIF binding to this locus using chromatin immunoprecipitation of HIF subunits. Consequentially, HIF-mediated DPF3 regulation was dependent on the presence of the risk allele. Finally, we show that DPF3 deletion in proximal tubular cells retarded cell growth, indicating potential roles for DPF3 in cell proliferation. Our analyses suggest that the HIF pathway differentially operates on a SNP-induced hypoxia-response element at 14q24.2, thereby affecting DPF3 expression, which increases the risk of developing renal cancer.

Master lineage transcription factors anchor trans mega transcriptional complexes at highly accessible enhancer sites to promote long-range chromatin clustering and transcription of distal target genes

The term 'super enhancers' (SE) has been widely used to describe stretches of closely localized enhancers that are occupied collectively by large numbers of transcription factors (TFs) and co-factors, and control the transcription of highly-expressed genes. Through integrated analysis of >600 DNase-seq, ChIP-seq, GRO-seq, STARR-seq, RNA-seq, Hi-C and ChIA-PET data in five human cancer cell lines, we identified a new class of autonomous SEs (aSEs) that are excluded from classic SE calls by the widely used Rank Ordering of Super-Enhancers (ROSE) method. TF footprint analysis revealed that compared to classic SEs and regular enhancers, aSEs are tightly bound by a dense array of master lineage TFs, which serve as anchors to recruit additional TFs and co-factors in trans. In addition, aSEs are preferentially enriched for Cohesins, which likely involve in stabilizing long-distance interactions between aSEs and their distal target genes. Finally, we showed that aSEs can be reliably predicted using a single DNase-seq data or combined with Mediator and/or P300 ChIP-seq. Overall, our study demonstrates that aSEs represent a unique class of functionally important enhancer elements that distally regulate the transcription of highly expressed genes.

A lysosome-targeted dextran-doxorubicin nanodrug overcomes doxorubicin-induced chemoresistance of myeloid leukemia

The hypoxic microenvironment is presumed to be a sanctuary for myeloid leukemia cells that causes relapse following chemotherapy, but the underlying mechanism remains elusive. Using a zebrafish xenograft model, we observed that the hypoxic hematopoietic tissue preserved most of the chemoresistant leukemic cells following the doxorubicin (Dox) treatment. And hypoxia upregulated TFEB, a master regulator of lysosomal biogenesis, and increased lysosomes in leukemic cells. Specimens from relapsed myeloid leukemia patients also harbored excessive lysosomes, which trapped Dox and prevented drug nuclear influx leading to leukemia chemoresistance. Pharmaceutical inhibition of lysosomes enhanced Dox-induced cytotoxicity against leukemic cells under hypoxia circumstance. To overcome lysosome associated chemoresistance, we developed a pH-sensitive dextran-doxorubicin nanomedicine (Dex-Dox) that efficiently released Dox from lysosomes and increased drug nuclear influx. More importantly, Dex-Dox treatment significantly improved the chemotherapy outcome in the zebrafish xenografts transplanted with cultured leukemic cells or relapsed patient specimens. Overall, we developed a novel lysosome targeting nanomedicine that is promising to overcome the myeloid leukemia chemoresistance.

Genome-wide profiling in colorectal cancer identifies PHF19 and TBC1D16 as oncogenic super enhancers

Colorectal cancer is one of the most common cancers in the world. Although genomic mutations and single nucleotide polymorphisms have been extensively studied, the epigenomic status in colorectal cancer patient tissues remains elusive. Here, together with genomic and transcriptomic analysis, we use ChIP-Seq to profile active enhancers at the genome wide level in colorectal cancer paired patient tissues (tumor and adjacent tissues from the same patients). In total, we sequence 73 pairs of colorectal cancer tissues and generate 147 H3K27ac ChIP-Seq, 144 RNA-Seq, 147 whole genome sequencing and 86 H3K4me3 ChIP-Seq samples. Our analysis identifies 5590 gain and 1100 lost variant enhancer loci in colorectal cancer, and 334 gain and 121 lost variant super enhancer loci. Multiple key transcription factors in colorectal cancer are predicted with motif analysis and core regulatory circuitry analysis. Further experiments verify the function of the super enhancers governing PHF19 and TBC1D16 in regulating colorectal cancer tumorigenesis, and KLF3 is identified as an oncogenic transcription factor in colorectal cancer. Taken together, our work provides an important epigenomic resource and functional factors for epigenetic studies in colorectal cancer.

Sirtuin 7 super-enhancer drives epigenomic reprogramming in hepatocarcinogenesis

Hepatocellular carcinoma (HCC) is a major cancer burden worldwide with increasing incidence in many developed countries. Super-enhancers (SEs) drive gene expressions required for cell type-specificity and tumor cell identity. However, their roles in HCC remain unclear because of data scarcity from primary tumors. Herein, chromatin profiling of non-alcoholic fatty liver disease (NAFLD)-associated HCCs and matched liver tissues uncovered an average of ∼500 somatically-acquired SEs per patient. The identified SE-target genes were functionally enriched for aberrant metabolism and cancer phenotypes, especially chromatin regulators including deacetylases and Polycomb repressive complexes. Notably, all examined tumors exhibited SE activation of Sirtuin 7 (SIRT7), genome-wide promoter H3K18 deacetylation and concurrent H3K27me3, as well as tumor-suppressor gene silencing. Depletion of SIRT7 SE in hepatoma cells induced global H3K18 acetylation and reactivated key metabolic and immune regulators, leading to marked suppression of tumorigenicity in vitro and in vivo. In concordance, SIRT7 physically interacted with the methyltransferase EZH2, and they were co-expressed in primary HCCs. In summary, our integrative analysis establishes a compendium of SEs in NAFLD-associated HCCs and uncovers SIRT7-driven chromatin regulatory network as potential druggable vulnerability of this increasingly prevalent cancer.

Identification of Hub Genes Associated With Clear Cell Renal Cell Carcinoma by Integrated Bioinformatics Analysis

Background

Clear cell renal cell carcinoma (ccRCC) is a common genitourinary cancer type with a high mortality rate. Due to a diverse range of biochemical alterations and a high level of tumor heterogeneity, it is crucial to select highly validated prognostic biomarkers to be able to identify subtypes of ccRCC early and apply precision medicine approaches.

Methods

Transcriptome data of ccRCC and clinical traits of patients were obtained from the GSE126964 dataset of Gene Expression Omnibus and The Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma (TCGA-KIRC) database. Weighted gene co-expression network analysis (WGCNA) and differentially expressed gene (DEG) screening were applied to detect common differentially co-expressed genes. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes analysis, survival analysis, prognostic model establishment, and gene set enrichment analysis were also performed. Immunohistochemical analysis results of the expression levels of prognostic genes were obtained from The Human Protein Atlas. Single-gene RNA sequencing data were obtained from the GSE131685 and GSE171306 datasets.

Results

In the present study, a total of 2,492 DEGs identified between ccRCC and healthy controls were filtered, revealing 1,300 upregulated genes and 1,192 downregulated genes. Using WGCNA, the turquoise module was identified to be closely associated with ccRCC. Hub genes were identified using the maximal clique centrality algorithm. After having intersected the hub genes and the DEGs in GSE126964 and TCGA-KIRC dataset, and after performing univariate, least absolute shrinkage and selection operator, and multivariate Cox regression analyses, ALDOB, EFHD1, and ESRRG were identified as significant prognostic factors in patients diagnosed with ccRCC. Single-gene RNA sequencing analysis revealed the expression profile of ALDOB, EFHD1, and ESRRG in different cell types of ccRCC.

Conclusions

The present results demonstrated that ALDOB, EFHD1, and ESRRG may act as potential targets for medical therapy and could serve as diagnostic biomarkers for ccRCC.

Targeting enhancer reprogramming to mitigate MEK inhibitor resistance in preclinical models of advanced ovarian cancer

Ovarian cancer is characterized by aberrant activation of the mitogen-activated protein kinase (MAPK), highlighting the importance of targeting the MAPK pathway as an attractive therapeutic strategy. However, the clinical efficacy of MEK inhibitors is limited by intrinsic or acquired drug resistance. Here, we established patient-derived ovarian cancer models resistant to MEK inhibitors and demonstrated that resistance to the clinically approved MEK inhibitor trametinib was associated with enhancer reprogramming. We also showed that enhancer decommissioning induced the downregulation of negative regulators of the MAPK pathway, leading to constitutive ERK activation and acquired resistance to trametinib. Epigenetic compound screening uncovered that HDAC inhibitors could alter the enhancer reprogramming and upregulate the expression of MAPK negative regulators, resulting in sustained MAPK inhibition and reversal of trametinib resistance. Consequently, a combination of HDAC inhibitor and trametinib demonstrated a synergistic antitumor effect in vitro and in vivo, including patient-derived xenograft mouse models. These findings demonstrated that enhancer reprogramming of the MAPK regulatory pathway might serve as a potential mechanism underlying MAPK inhibitor resistance and concurrent targeting of epigenetic pathways and MAPK signaling might provide an effective treatment strategy for advanced ovarian cancer.

Integrative epigenomic and high-throughput functional enhancer profiling reveals determinants of enhancer heterogeneity in gastric cancer

BACKGROUND

Enhancers are distal cis-regulatory elements required for cell-specific gene expression and cell fate determination. In cancer, enhancer variation has been proposed as a major cause of inter-patient heterogeneity-however, most predicted enhancer regions remain to be functionally tested.

METHODS

We analyzed 132 epigenomic histone modification profiles of 18 primary gastric cancer (GC) samples, 18 normal gastric tissues, and 28 GC cell lines using Nano-ChIP-seq technology. We applied Capture-based Self-Transcribing Active Regulatory Region sequencing (CapSTARR-seq) to assess functional enhancer activity. An Activity-by-contact (ABC) model was employed to explore the effects of histone acetylation and CapSTARR-seq levels on enhancer-promoter interactions.

RESULTS

We report a comprehensive catalog of 75,730 recurrent predicted enhancers, the majority of which are GC-associated in vivo (> 50,000) and associated with lower somatic mutation rates inferred by whole-genome sequencing. Applying CapSTARR-seq to the enhancer catalog, we observed significant correlations between CapSTARR-seq functional activity and H3K27ac/H3K4me1 levels. Super-enhancer regions exhibited increased CapSTARR-seq signals compared to regular enhancers, even when decoupled from native chromatin contexture. We show that combining histone modification and CapSTARR-seq functional enhancer data improves the prediction of enhancer-promoter interactions and pinpointing of germline single nucleotide polymorphisms (SNPs), somatic copy number alterations (SCNAs), and trans-acting TFs involved in GC expression. We identified cancer-relevant genes (ING1, ARL4C) whose expression between patients is influenced by enhancer differences in genomic copy number and germline SNPs, and HNF4α as a master trans-acting factor associated with GC enhancer heterogeneity.

CONCLUSIONS

Our results indicate that combining histone modification and functional assay data may provide a more accurate metric to assess enhancer activity than either platform individually, providing insights into the relative contribution of genetic (cis) and regulatory (trans) mechanisms to GC enhancer functional heterogeneity.

Functional succinate dehydrogenase deficiency is a common adverse feature of clear cell renal cancer

This study demonstrates that underexpression of succinate dehydrogenase (SDH) subunits resulting in accumulation of oncogenic succinate is a common, adverse, epigenetic modulating feature in clear cell renal cell carcinoma (ccRCC), during pathogenesis and progression. The study sheds light on the mechanisms of down-regulation of SDH subunits in ccRCC and deciphers the consequent oncogenic effects. It shows that functional SDH deficiency is a common feature of ccRCC (∼80% of all kidney cancers), and not just limited to the 0.05 to 0.5% of kidney cancers with germline SDH mutations.

Reduced succinate dehydrogenase (SDH) activity resulting in adverse succinate accumulation was previously considered relevant only in 0.05 to 0.5% of kidney cancers associated with germline SDH mutations. Here, we sought to examine a broader role for SDH loss in kidney cancer pathogenesis/progression. We report that underexpression of SDH subunits resulting in accumulation of oncogenic succinate is a common feature in clear cell renal cell carcinoma (ccRCC) (∼80% of all kidney cancers), with a marked adverse impact on survival in ccRCC patients (n = 516). We show that SDH down-regulation is a critical brake in the TCA cycle during ccRCC pathogenesis and progression. In exploring mechanisms of SDH down-regulation in ccRCC, we report that Von Hippel-Lindau loss-induced hypoxia-inducible factor–dependent up-regulation of miR-210 causes direct inhibition of the SDHD transcript. Moreover, shallow deletion of SDHB occurs in ∼20% of ccRCC. We then demonstrate that SDH loss-induced succinate accumulation contributes to adverse loss of 5-hydroxymethylcytosine, gain of 5-methylcytosine, and enhanced invasiveness in ccRCC via inhibition of ten-eleven translocation (TET)-2 activity. Intriguingly, binding affinity between the catalytic domain of recombinant TET-2 and succinate was found to be very low, suggesting that the mechanism of succinate-induced attenuation of TET-2 activity is likely via product inhibition rather than competitive inhibition. Finally, exogenous ascorbic acid, a TET-activating demethylating agent, led to reversal of the above oncogenic effects of succinate in ccRCC cells. Collectively, our study demonstrates that functional SDH deficiency is a common adverse feature of ccRCC and not just limited to the kidney cancers associated with germline SDH mutations.

Loss of Von Hippel-Lindau (VHL) Tumor Suppressor Gene Function: VHL-HIF Pathway and Advances in Treatments for Metastatic Renal Cell Carcinoma (RCC)

Renal cell carcinoma (RCC) is a malignancy of the kidney originating from the tubular epithelium. Inactivation of the von Hippel–Lindau tumor-suppressor gene (VHL) is found in most clear cell renal cell carcinomas (ccRCCs). The VHL–HIF–VEGF/VEGFR pathway, which involves the von Hippel–Lindau tumor suppressor protein (VHL), hypoxia-inducible factor (HIF), vascular endothelial growth factor (VEGF), and its receptor (VEGFR), is a well-studied therapeutic target for metastatic ccRCC. Therefore, over the past decade, anti-angiogenic agents targeting VEGFR have served as the standard treatment for metastatic RCC. Recently, based on the immunomodulatory effect of anti-VEGFR therapy, anti-angiogenic agents and immune checkpoint inhibitor combination strategies have also emerged as therapeutic strategies. These advances were made possible by the improved understanding of the VHL–HIF pathway. In this review, we summarize the historical evolution of ccRCC treatments, with a focus on the involvement of the VHL–HIF pathway.

Mutational Analysis of PBRM1 and Significance of PBRM1 Mutation in Anti-PD-1 Immunotherapy of Clear Cell Renal Cell Carcinoma

Renal cell carcinoma is a common solid tumor. PBRM1 is one of the most mutation-prone genes in clear cell renal cell carcinoma (ccRCC) with the occurrence of mutation in 40% of ccRCC patients. Mutations in PBRM1 have been correlated with the efficacy of immunotherapy. However, the mutation types of PBRM1 are not well characterized. The effects of PBRM1 expression levels in the tumor microenvironment are not well studied. In addition, the mechanism and effect of anti-PD-1 immunotherapy in ccRCC tumor microenvironments are not well clarified. In this study, using bioinformatics methods we analyzed the alternation frequency and expression levels of PBRM1 in various tumors. Next, we experimentally validated their expression levels in ccRCC tissues from human and mouse models. We attempted to clarify the mechanisms of anti-PD-1 immunotherapy in ccRCC with various PBRM1 expression levels. Our results showed that deficiency of PBRM1 protein is correlated with CD4 T cell reduction in human and mouse ccRCC tissues. We also showed that anti-PD-1 Immunotherapy can increase the infiltration of T cells in both PBRM1 high and PBRM1 low tumors but to different degrees. Our study indicates that the reduction of CD4 cells in tumor tissues with low expression of PBRM1 may explain the compromised efficacy of anti-PD-1 immunotherapy in patients with PBRM1 mutated ccRCC. Our study sheds light on the potential of PBRM1 as a therapeutic target in ccRCC.

Genome-Wide Enhancer Analysis Reveals the Role of AP-1 Transcription Factor in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is one of the most common cancers in the world, but its epigenomic features have not been determined. Here, we studied the chromatin landscape of active enhancers of HNSCC head tumor tissues by performing H3K27ac and H3K4me1 ChIP-Seq with a Tgfbr1/Pten double conditional knockout HNSCC mouse model. We identified 1,248 gain variant enhancer loci (VELs) and 2,188 lost VELs, as well as 153 gain variant super enhancer loci (VSELs) and 234 lost VSELs. Potentially involved transcription factors were predicted with motif analysis, and we identified AP-1 as one of the critical oncogenic transcription factors in HNSCC and many other types of cancer. Combining transcriptomic and epigenomic data, our analysis also showed that AP-1 and histone modifications coordinately regulate target gene expression in HNSCC. In conclusion, our study provides important epigenomic information for enhancer studies in HNSCC and reveals new mechanism for AP-1 regulating HNSCC.

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.

A selective HDAC8 inhibitor potentiates antitumor immunity and efficacy of immune checkpoint blockade in hepatocellular carcinoma

Insufficient T cell infiltration into noninflamed tumors, such as hepatocellular carcinoma (HCC), restricts the effectiveness of immune-checkpoint blockade (ICB) for a subset of patients. Epigenetic therapy provides further opportunities to rewire cancer-associated transcriptional programs, but whether and how selective epigenetic inhibition counteracts the immune-excluded phenotype remain incompletely defined. Here, we showed that pharmacological inhibition of histone deacetylase 8 (HDAC8), a histone H3 lysine 27 (H3K27)-specific isozyme overexpressed in a variety of human cancers, thwarts HCC tumorigenicity in a T cell-dependent manner. The tumor-suppressive effect of selective HDAC8 inhibition was abrogated by CD8⁺ T cell depletion or regulatory T cell adoptive transfer. Chromatin profiling of human HDAC8-expressing HCCs revealed genome-wide H3K27 deacetylation in 1251 silenced enhancer-target gene pairs that are enriched in metabolic and immune regulators. Mechanistically, down-regulation of HDAC8 increased global and enhancer acetylation of H3K27 to reactivate production of T cell-trafficking chemokines by HCC cells, thus relieving T cell exclusion in both immunodeficient and humanized mouse models. In an HCC preclinical model, selective HDAC8 inhibition increased tumor-infiltrating CD8⁺ T cells and potentiated eradication of established hepatomas by anti-PD-L1 therapy without evidence of toxicity. Mice treated with HDAC8 and PD-L1 coblockade were protected against subsequent tumor rechallenge as a result of the induction of memory T cells and remained tumor-free for greater than 15 months. Collectively, our study demonstrates that selective HDAC8 inhibition elicits effective and durable responses to ICB by co-opting adaptive immunity through enhancer reprogramming.

ACSNI: An unsupervised machine-learning tool for prediction of tissue-specific pathway components using gene expression profiles

Determining the tissue- and disease-specific circuit of biological pathways remains a fundamental goal of molecular biology. Many components of these biological pathways still remain unknown, hindering the full and accurate characterization of biological processes of interest. Here we describe ACSNI, an algorithm that combines prior knowledge of biological processes with a deep neural network to effectively decompose gene expression profiles (GEPs) into multi-variable pathway activities and identify unknown pathway components. Experiments on public GEP data show that ACSNI predicts cogent components of mTOR, ATF2, and HOTAIRM1 signaling that recapitulate regulatory information from genetic perturbation and transcription factor binding datasets. Our framework provides a fast and easy-to-use method to identify components of signaling pathways as a tool for molecular mechanism discovery and to prioritize genes for designing future targeted experiments (https://github.com/caanene1/ACSNI).

Dysregulation of Zinc Finger Protein 395 Contributes to the Pathogenesis of Chondrosarcoma

Introduction

The transcription factor zinc finger protein 395 (ZNF395) is involved in several cellular responses and tumorigenesis. However, the potential role and clinical significance of ZNF395 in chondrosarcoma are not well investigated. This study determines the expression profile, prognostic value and biological function of ZNF395 in human chondrosarcoma.

Methods

The mRNA and protein expressions of ZNF395 in fresh chondrosarcomas and the matched adjacent non-tumor tissues were assessed using real-time PCR and immunoblotting, respectively. The protein expression of ZNF395 in chondrosarcoma specimens was evaluated by immunohistochemistry, and the relationships among its protein level, clinicopathological parameters and prognosis were further detected. Cell viability, colony formation, migration, invasion and apoptosis assay were evaluated in chondrosarcoma cells with depletion of ZNF395.

Results

The mRNA and protein expressions of ZNF395 in chondrosarcoma tissues were significantly higher than those in the matched adjacent non-tumor tissues and benign cartilage tumors. Clinical analysis displayed that ZNF395 was expressed at higher levels in chondrosarcoma patients with higher histological grade and advanced MSTS stage. Furthermore, we demonstrated that high expression of ZNF395 correlated with a worse overall survival of chondrosarcoma patients. Multivariate Cox regression analysis indicated that ZNF395 was an independent prognostic marker in chondrosarcoma patients. Functional studies revealed that depletion of ZNF395 markedly inhibited cell viability, colony formation, migration and invasion, and promoted apoptosis in chondrosarcoma.

Conclusion

These findings suggest that dysregulation of ZNF395 contributes to chondrosarcoma development, and ZNF395 may act as a potent oncogene and serve as a independently prognostic factor, highlight the potential of ZNF395 as a novel biomarker and therapeutic target for chondrosarcoma.

Transcription Factors BARX1 and DLX4 Contribute to Progression of Clear Cell Renal Cell Carcinoma via Promoting Proliferation and Epithelial-Mesenchymal Transition

Dysregulation of transcription factors contributes to the carcinogenesis and progression of cancers. However, their roles in clear cell renal cell carcinoma remain largely unknown. This study aimed to evaluate the clinical significance of TFs and investigate their potential molecular mechanisms in ccRCC. Data were accessed from the cancer genome atlas kidney clear cell carcinoma cohort. Bioinformatics algorithm was used in copy number alterations mutations, and differentially expressed TFs’ analysis. Univariate and multivariate Cox regression analyses were performed to identify clinically significant TFs and construct a six-TF prognostic panel. TFs’ expression was validated in human tissues. Gene set enrichment analysis (GSEA) was utilized to find enriched cancer hallmark pathways. Functional experiments were conducted to verify the cancer-promoting effect of BARX homeobox 1 (BARX1) and distal-less homeobox 4 (DLX4) in ccRCC, and Western blot was performed to explore their downstream pathways. As for results, many CNAs and mutations were identified in transcription factor genes. TFs were differentially expressed in ccRCC. An applicable predictive panel of six-TF genes was constructed to predict the overall survival for ccRCC patients, and its diagnostic efficiency was evaluated by the area under the curve (AUC). BARX1 and DLX4 were associated with poor prognosis, and they could promote the proliferation and migration of ccRCC. In conclusion, the six-TF panel can be used as a prognostic biomarker for ccRCC patients. BARX1 and DLX4 play oncogenic roles in ccRCC via promoting proliferation and epithelial–mesenchymal transition. They have the potential to be novel therapeutic targets for ccRCC.

HIF-1 directly induces TET3 expression to enhance 5-hmC density and induce erythroid gene expression in hypoxia

In mammalian cells, cytosines found within cytosine guanine dinucleotides can be methylated to 5-methylcytosine (5-mC) by DNA methyltransferases and further oxidized by the Ten-eleven translocation dioxygenase (TET) enzymes to 5-hydroxymethylcytosine (5-hmC). We have previously shown that hematopoietic stem and progenitor cells (HSPCs) with TET2 mutations have aberrant 5-hmC distribution and less erythroid differentiation potential. However, these experiments were performed under standard tissue culture conditions with 21% oxygen (O2), whereas HSPCs in human bone marrow reside in ∼1% O2. Therefore, to model human erythropoiesis more accurately, we compared 5-hmC distribution and gene expression in hypoxic vs normoxic conditions. Despite TET enzymes having limited O2 as a substrate in hypoxia, 5-hmC peaks were more numerous and pronounced than in normoxia. Among the TET genes, TET3 was upregulated specifically in hypoxia. We identified 2 HIF-1 binding sites in TET3 by chromatin immunoprecipitation of HIF-1α followed by sequencing, and TET3 upregulation was abrogated with deletion of both sites, indicating that TET3 is a direct HIF-1 target. Finally, we showed that loss of one or both of these HIF-1 binding sites in K562 cells disrupted erythroid differentiation in hypoxia and lowered cell viability. This work provides a molecular link between O2 availability, epigenetic modification of chromatin, and erythroid differentiation.