High-resolution profiling of histone h3 lysine 36 trimethylation in metastatic renal cell carcinoma

Deciphering the interplay between SETD2 mediated H3K36me3 and RNA N6-methyladenosine in clear cell renal cell carcinoma (ccRCC)

RNA N6-methyladenosine (m6A) plays diverse roles in RNA metabolism and its deregulation contributes to tumor initiation and progression. Clear cell renal cell carcinoma (ccRCC) is characterized by near ubiquitous loss of VHL followed by mutations in epigenetic regulators PBRM1, SETD2, and BAP1. Mutations in SETD2, a histone H3 lysine 36 trimethylase (H3K36me3), are associated with reduced survival, greater metastatic propensity, and metabolic reprogramming. While m6A and H3K36me3 deregulation are separately implicated in renal tumorigenesis, H3K36me3 may participate directly in m6A targeting, but the m6A-H3K36me3 interplay has not been investigated in the context of ccRCC. Using RCC-relevant SETD2 isogenic knockout and rescue cell line models, we demonstrate a dynamic redistribution of m6A in the SETD2 depleted transcriptome, with a subset of transcripts involved in metabolic reprogramming demonstrating SETD2 dependent m6A and expression level changes. Using a panel of six histone modifications we show that m6A redistributes to regions enriched in gained active enhancers upon SETD2 inactivation. Finally, we demonstrate a reversal of transcriptomic programs involved in SETD2 loss mediated metabolic reprogramming, and reduced cell viability through pharmacologic inhibition or genetic ablation of m6A writer METTL3 specific to SETD2 deficient cells. Thus, targeting m6A may represent a novel therapeutic vulnerability in SETD2 mutant ccRCC.

Extensive genotype-phenotype heterogeneity in renal cell carcinoma - a proof-of-concept study

Background

Renal cell carcinoma (RCC) is characterized by a high degree of genomic but also functional intratumoral heterogeneity (ITH). Mutations in VHL, chromatin remodeling genes such as SETD2 and genes that regulate the PI3K/AKT/mTOR pathway have been identified as recurrent drivers despite genomic ITH. Whether and to what extent these mutations shape functional ITH including the formation of spatial niches is incompletely understood. Herein, we analyze the correlation between mutational drivers and their functional proxies in a spatially defined manner.

Methods

A total of 23 RCCs were analyzed by panel next-generation sequencing followed by immunohistochemistry for five functional proxies for key genetic alterations including the expression of CD31, GLUT1, phospho-mTOR S2448, H3K36me3 and Ki-67. Antibody stainings were scored semiquantitatively in the tumor periphery and the tumor center.

Results

Unexpectedly, the presence of a VHL mutation was not found to correlate with its functional proxies including the expression of CD31/microvessel density or the expression of the glucose transporter GLUT1. Likewise, there was no correlation between the presence of activating mutations in genes of the PI3K/AKT/mTOR pathway and the expression of activated phospho-mTOR S2448. Furthermore, mutations in the methyltransferase gene SETD2 were not found to correlate with the expression level of its downstream target H3K36me3. Lastly, there was no correlation between the expression of the proliferation marker Ki-67 and the number of driver mutations.

Conclusion

This proof-of-concept study adds genotype-phenotype heterogeneity as additional layer of complexity to the known genomic and functional ITH in RCC.

SETD2 loss-of-function uniquely sensitizes cells to epigenetic targeting of NSD1-directed H3K36 methylation

BACKGROUND

SETD2 is the sole epigenetic factor responsible for catalyzing histone 3, lysine 36, tri-methylation (H3K36me3) in mammals. Its role in regulating cellular processes such as RNA splicing, DNA repair, and spurious transcription initiation underlies its broader tumor suppressor function. SETD2 mutation promotes the epithelial-mesenchymal transition and is clinically associated with adverse outcomes highlighting a therapeutic need to develop targeted therapies against this dangerous mutation.

RESULTS

We employ an unbiased genome-wide synthetic lethal screen, which identifies another H3K36me writer, NSD1, as a synthetic lethal modifier in SETD2-mutant cells. Confirmation of this synthetic lethal interaction is performed in isogenic clear cell renal cell carcinoma and immortalized renal epithelial cell lines, in mouse and human backgrounds. Depletion of NSD1 using a CRISPRi targeting approach promotes the loss of SETD2-mutant cells coincident with elevated levels of DNA damage and apoptosis. Surprisingly, only suppression of NSD1, but not related H3K36-methyltransferases, promotes synthetic lethality in these models. Mapping of genomic H3K36me2 targeting by NSD1 and NSD2 individually highlights the independent functions of these epigenetic writers. Furthermore, as a proof-of-principle, we demonstrate the therapeutic feasibility of targeting this synthetic lethal interaction by recapitulating the phenotype using BT5, a first-in-class pharmacologic inhibitor against NSD1.

CONCLUSIONS

These findings unify genome-wide screening approaches with the latest genetic and pharmacologic modeling methodologies to reveal an entirely novel epigenetic approach to individualize therapies against a challenging loss-of-function SETD2 mutation in cancer.

Structure and function of the lysine methyltransferase SETD2 in cancer: From histones to cytoskeleton

SETD2 is known to be the unique histone methyltransferase responsible for the trimethylation of the lysine 36 of histone H3 thus generating H3K36me3. This epigenetic mark is critical for transcriptional activation and elongation, DNA repair, mRNA splicing, and DNA methylation. Recurrent SETD2-inactivating mutations and altered H3K36me3 levels are found in cancer at high frequency and numerous studies indicate that SETD2 acts as a tumor suppressor. Recently, SETD2 was further shown to methylate non-histone proteins particularly the cytoskeletal proteins tubulin and actin with subsequent impacts on cytoskeleton structure, mitosis and cell migration. Herein, we provide a review of the role of SETD2 in different cancers with special emphasis on the structural basis of the functions of this key lysine methyltransferase. Moreover, beyond the role of this enzyme in epigenetics and H3K36me3-dependent processes, we highlight the putative role of "non-epigenetic/H3K36me3" functions of SETD2 in cancer, particularly those involving the cytoskeleton.

Revealing the characteristics of SETD2-mutated clear cell renal cell carcinoma through tumor heterogeneity analysis

Background

Renal cell carcinoma (RCC) is the most prevalent type of malignant kidney tumor in adults, with clear cell renal cell carcinoma (ccRCC) comprising about 75% of all cases. The SETD2 gene, which is involved in the modification of histone proteins, is often found to have alterations in ccRCC. Yet, our understanding of how these SETD2 mutations affect ccRCC characteristics and behavior within the tumor microenvironment is still not fully understood.

Methods

We conducted a detailed analysis of single-cell RNA sequencing (scRNA-seq) data from ccRCC. First, the data was preprocessed using the Python package, "scanpy." High variability genes were pinpointed through Pearson's correlation coefficient. Dimensionality reduction and clustering identification were performed using Principal Component Analysis (PCA) and the Leiden algorithm. Malignant cell identification was conducted with the "InferCNV" R package, while cell trajectories and intercellular communication were depicted using the Python packages "VIA" and "cellphoneDB." We then employed the R package "Deseq2" to determine differentially expressed genes (DEGs) between groups. Using high-dimensional weighted gene correlation network analysis (hdWGCNA), co-expression modules were identified. We intersected these modules with DEGs to establish prognostic models through univariate Cox and the least absolute shrinkage and selection operator (LASSO) method.

Results

We identified 69 and 53 distinctive cell clusters, respectively. These were classified further into 12 unique cell types. This analysis highlighted the presence of an abnormal tumor sub-cluster (MT + group), identified by high mitochondrial-encoded protein gene expression and an indication of unfavorable prognosis. Investigation of cellular interactions spotlighted significant interactions between the MT + group and endothelial cells, macrophaes. In addition, we developed a prognostic model based on six characteristic genes. Notably, risk scores derived from these genes correlated significantly with various clinical features. Finally, a nomogram model was established to facilitate more accurate outcome prediction, incorporating four independent risk factors.

Conclusion

Our findings provide insight into the crucial transcriptomic characteristics of ccRCC associated with SETD2 mutation. We discovered that this mutation-induced subcluster could stimulate M2 polarization in macrophages, suggesting a heightened propensity for metastasis. Moreover, our prognostic model demonstrated effectiveness in forecasting overall survival for ccRCC patients, thus presenting a valuable clinical tool.

SETD2 loss in renal epithelial cells drives epithelial-to-mesenchymal transition in a TGF-β-independent manner

Histone-lysine N-methyltransferase SETD2 (SETD2), the sole histone methyltransferase that catalyzes trimethylation of lysine 36 on histone H3 (H3K36me3), is often mutated in clear cell renal cell carcinoma (ccRCC). SETD2 mutation and/or loss of H3K36me3 is linked to metastasis and poor outcome in ccRCC patients. Epithelial-to-mesenchymal transition (EMT) is a major pathway that drives invasion and metastasis in various cancer types. Here, using novel kidney epithelial cell lines isogenic for SETD2, we discovered that SETD2 inactivation drives EMT and promotes migration, invasion, and stemness in a transforming growth factor-beta-independent manner. This newly identified EMT program is triggered in part through secreted factors, including cytokines and growth factors, and through transcriptional reprogramming. RNA-seq and assay for transposase-accessible chromatin sequencing uncovered key transcription factors upregulated upon SETD2 loss, including SOX2, POU2F2 (OCT2), and PRRX1, that could individually drive EMT and stemness phenotypes in SETD2 wild-type (WT) cells. Public expression data from SETD2 WT/mutant ccRCC support the EMT transcriptional signatures derived from cell line models. In summary, our studies reveal that SETD2 is a key regulator of EMT phenotypes through cell-intrinsic and cell-extrinsic mechanisms that help explain the association between SETD2 loss and ccRCC metastasis.

SETD2 Deficiency Confers Sensitivity to Dual Inhibition of DNA Methylation and PARP in Kidney Cancer

SETD2 deficiency alters the epigenetic landscape by causing depletion of H3K36me3 and plays an important role in diverse forms of cancer, most notably in aggressive and metastatic clear-cell renal cell carcinomas (ccRCC). Development of an effective treatment scheme targeting SETD2-compromised cancer is urgently needed. Considering that SETD2 is involved in DNA methylation and DNA repair, a combination treatment approach using DNA hypomethylating agents (HMA) and PARP inhibitors (PARPi) could have strong antitumor activity in SETD2-deficient kidney cancer. We tested the effects of the DNA HMA 5-aza-2'-dexoxydytidine (DAC), the PARPi talazoparib (BMN-673), and both in combination in human ccRCC models with or without SETD2 deficiency. The combination treatment of DAC and BMN-673 synergistically increased cytotoxicity in vitro in SETD2-deficient ccRCC cell lines but not in SETD2-proficient cell lines. DAC and BMN-673 led to apoptotic induction, increased DNA damage, insufficient DNA damage repair, and increased genomic instability. Furthermore, the combination treatment elevated immune responses, upregulated STING, and enhanced viral mimicry by activating transposable elements. Finally, the combination effectively suppressed the growth of SETD2-deficient ccRCC in in vivo mouse models. Together, these findings indicate that combining HMA and PARPi is a promising potential therapeutic strategy for treating SETD2-compromised ccRCC.

SIGNIFICANCE

SETD2 deficiency creates a vulnerable epigenetic status that is targetable using a DNA hypomethylating agent and PARP inhibitor combination to suppress renal cell carcinoma, identifying a precision medicine-based approach for SETD2-compromised cancers.

The arginine methyltransferase Prmt1 coordinates the germline arginine methylome essential for spermatogonial homeostasis and male fertility

Arginine methylation, catalyzed by the protein arginine methyltransferases (PRMTs), is a common post-translational protein modification (PTM) that is engaged in a plethora of biological events. However, little is known about how the methylarginine-directed signaling functions in germline development. In this study, we discover that Prmt1 is predominantly distributed in the nuclei of spermatogonia but weakly in the spermatocytes throughout mouse spermatogenesis. By exploiting a combination of three Cre-mediated Prmt1 knockout mouse lines, we unravel that Prmt1 is essential for spermatogonial establishment and maintenance, and that Prmt1-catalyzed asymmetric methylarginine coordinates inherent transcriptional homeostasis within spermatogonial cells. In conjunction with high-throughput CUT&Tag profiling and modified mini-bulk Smart-seq2 analyses, we unveil that the Prmt1-deposited H4R3me2a mark is permissively enriched at promoter and exon/intron regions, and sculpts a distinctive transcriptomic landscape as well as the alternative splicing pattern, in the mouse spermatogonia. Collectively, our study provides the genetic and mechanistic evidence that connects the Prmt1-deposited methylarginine signaling to the establishment and maintenance of a high-fidelity transcriptomic identity in orchestrating spermatogonial development in the mammalian germline.

Investigating the Role of FoxP3 in Renal Cell Carcinoma Metastasis with BAP1 or SEDT2 Mutation

Forkhead box protein P3 (FoxP3) primarily functions as the master regulator in regulatory T cells (Tregs) differentiation, but its high level of expression has also been found in tumor cells recently. The aim of our study was to clarify the role of FoxP3 in renal cell carcinoma (RCC) progression and metastasis. We verified the FoxP3 characteristic clinicopathological data from The Cancer Genome Atlas (TCGA) database using bioinformatics tools. Meanwhile, RNA sequencing was performed to determine the FoxP3 biofunction in RCC progression. Our results showed that high expression of FoxP3 was found in BAP1- or SETD2-mutant patients with RCC, and a higher FoxP3 expression was related to worse prognosis. However, there was no statistically significant relationship between the FoxP3 IHC score and RCC malignant progression owning to the limited number of patients in our tissue microarray. Using in vitro FoxP3 loss-of-function assays, we verified that silencing FoxP3 in 786-O and ACHN cells could inhibit the cell migration/invasion capability, which was consistent with the data from RNA sequencing in 786-O cells and from the TCGA datasets. Using an in vivo nude mice orthotopic kidney cancer model, we found that silencing FoxP3 could inhibit tumor growth. In conclusion, our study demonstrated that BAP1 or SEDT2 mutation could lead to higher expression of FoxP3 in RCC patients, and FoxP3 could eventually stimulate RCC cells' invasion and metastasis, which might indicate that FoxP3 could function as a potential oncogene in RCC progression.

Competing endogenous RNA analysis identified lncRNA DSCR9 as a novel prognostic biomarker associated with metastasis and tumor microenvironment in renal cell carcinoma

Distant metastasis is the main cause of death in patients with clear cell renal carcinoma (ccRCC). The dysregulation of the tumor microenvironment is responsible for tumorigenesis and metastasis in ccRCC. The role of long non-coding RNAs in the tumor immune of ccRCC remains unclear. The present study screened differentially expressed protein-coding genes and non-coding genes between ccRCC and normal tissues based on three datasets. The commonly deregulated genes were used to identify distant metastasis-related long non-coding RNAs (lncRNAs) and prognostic lncRNAs. Pearson correlation analysis was used to identify immune-related lncRNAs. A competing endogenous RNA network was constructed and hub lncRNAs were identified. A total of 1650 coding genes, 821 lncRNAs and 62 miRNAs were commonly deregulated in the three datasets. A total of 408 lncRNAs associated with the overall survival of patients with ccRCC were identified. Among them, 82 lncRNAs were distant metastasis-related. Further analysis identified 52 lncRNAs associated with the immune pathway. Functional analyses concordantly demonstrated the role of the 52 lncRNAs in metastasis and tumor immunology. The ceRNA network analysis indicated lncRNA DSCR9 as the key lncRNA regulator. Univariate and multivariate analysis in two independent cohorts validated that DSCR9 could be an independent risk factor for the progression-free survival of patients with ccRCC. Further analyses indicated that DSCR9 might be associated with the immunotherapeutic response. reverse transcription-quantitative PCR demonstrated that the RNA expression level of DSCR9 was upregulated in ccRCC compared with normal kidney samples. The present study demonstrated the potential of LncRNA DSCR9 in assessing the prognosis and developing future immunotherapy for patients with metastatic ccRCC.

Targeted Profiling of Epitranscriptomic Reader, Writer, and Eraser Proteins Regulated by H3K36me3

Trimethylation of lysine 36 on histone H3 (H3K36me3), an epigenetic mark associated with actively transcribed genes, plays an important role in multiple cellular processes, including transcription elongation, DNA methylation, DNA repair, etc. Aberrant expression and mutations of the main methyltransferase for H3K36me3, i.e., SET domain-containing 2 (SETD2), were shown to be associated with various cancers. Here, we performed targeted profiling of 154 epitranscriptomic reader, writer, and eraser (RWE) proteins using a scheduled liquid chromatography-parallel-reaction monitoring (LC-PRM) method coupled with the use of stable isotope-labeled (SIL) peptides as internal standards to investigate how H3K36me3 modulates the chromatin occupancies of epitranscriptomic RWE proteins. Our results showed consistent changes in chromatin occupancies of RWE proteins upon losses of H3K36me3 and H4K16ac and a role of H3K36me3 in recruiting METTL3 to chromatin following induction of DNA double-strand breaks. In addition, protein-protein interaction network and Kaplan-Meier survival analyses revealed the importance of METTL14 and TRMT11 in kidney cancer. Taken together, our work unveiled cross-talks between histone epigenetic marks (i.e., H3K36me3 and H4K16ac) and epitranscriptomic RWE proteins and uncovered the potential roles of these RWE proteins in H3K36me3-mediated biological processes.

SETD2-H3K36ME3: an important bridge between the environment and tumors

Epigenetic regulation plays an important role in the occurrence, development and treatment of tumors. The histone methyltransferase SET-domain-containing 2 (SETD2) plays a key role in mammalian epigenetic regulation by catalyzing histone methylation and interacting with RNA polymerase II to mediate transcription elongation and mismatch repair. As an important bridge between the environment and tumors, SETD2-H3K36me3 plays an important role in the occurrence and development of tumors. Many tumors, including renal cancer, gastric cancer, lung cancer, are closely related to SETD2 gene mutations. As a key component of common tumor suppressor mechanisms, SETD2-H3K36me3is an important target for clinical disease diagnosis and treatment. Here, we reviewed the structure and function of the SETD2 and how SETD2-H3K36me3 functions as a bridge between the environment and tumors to provide an in-depth understanding of its role in the occurrence and development of various tumors, which is of great significance for future disease diagnosis and treatment.

The Chromatin Remodeler ATRX: Role and Mechanism in Biology and Cancer

The alpha-thalassemia mental retardation X-linked (ATRX) syndrome protein is a chromatin remodeling protein that primarily promotes the deposit of H3.3 histone variants in the telomere area. ATRX mutations not only cause ATRX syndrome but also influence development and promote cancer. The primary molecular characteristics of ATRX, including its molecular structures and normal and malignant biological roles, are reviewed in this article. We discuss the role of ATRX in its interactions with the histone variant H3.3, chromatin remodeling, DNA damage response, replication stress, and cancers, particularly gliomas, neuroblastomas, and pancreatic neuroendocrine tumors. ATRX is implicated in several important cellular processes and serves a crucial function in regulating gene expression and genomic integrity throughout embryogenesis. However, the nature of its involvement in the growth and development of cancer remains unknown. As mechanistic and molecular investigations on ATRX disclose its essential functions in cancer, customized therapies targeting ATRX will become accessible.

Genomic-transcriptomic evolution in lung cancer and metastasis

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy1. Here, using paired whole-exome and RNA sequencing data, we investigate intratumour transcriptomic diversity in 354 non-small cell lung cancer tumours from 347 out of the first 421 patients prospectively recruited into the TRACERx study2,3. Analyses of 947 tumour regions, representing both primary and metastatic disease, alongside 96 tumour-adjacent normal tissue samples implicate the transcriptome as a major source of phenotypic variation. Gene expression levels and ITH relate to patterns of positive and negative selection during tumour evolution. We observe frequent copy number-independent allele-specific expression that is linked to epigenomic dysfunction. Allele-specific expression can also result in genomic-transcriptomic parallel evolution, which converges on cancer gene disruption. We extract signatures of RNA single-base substitutions and link their aetiology to the activity of the RNA-editing enzymes ADAR and APOBEC3A, thereby revealing otherwise undetected ongoing APOBEC activity in tumours. Characterizing the transcriptomes of primary-metastatic tumour pairs, we combine multiple machine-learning approaches that leverage genomic and transcriptomic variables to link metastasis-seeding potential to the evolutionary context of mutations and increased proliferation within primary tumour regions. These results highlight the interplay between the genome and transcriptome in influencing ITH, lung cancer evolution and metastasis.

The role of histone methylation in renal cell cancer: an update

Renal cell carcinoma accounts for 2-3% of all cancers. It is difficult to diagnose early. Recently, genome-wide studies have identified that histone methylation was one of the functional classes that is most frequently dysregulated in renal cell cancer. Mutation or mis-regulation of histone methylation, methyltransferases, demethylases are associated with gene expression and tumor progression in renal cell cancer. Herein, we summarize histone methylations, demethylases and their alterations and mechanisms in renal cell cancer.

PBRM1, SETD2 and BAP1 - the trinity of 3p in clear cell renal cell carcinoma

Biallelic inactivation of the tumour suppressor gene Von Hippel-Lindau (VHL) occurs in the vast majority of clear cell renal cell carcinoma (ccRCC) instances, disrupting cellular oxygen-sensing mechanisms to yield a state of persistent pseudo-hypoxia, defined as a continued hypoxic response despite the presence of adequate oxygen levels. However, loss of VHL alone is often insufficient to drive oncogenesis. Results from genomic studies have shown that co-deletions of VHL with one (or more) of three genes encoding proteins involved in chromatin modification and remodelling, polybromo-1 gene (PBRM1), BRCA1-associated protein 1 (BAP1) and SET domain-containing 2 (SETD2), are common and important co-drivers of tumorigenesis. These genes are all located near VHL on chromosome 3p and are often altered following cytogenetic rearrangements that lead to 3p loss and precede the establishment of ccRCC. These three proteins have multiple roles in the regulation of crucial cancer-related pathways, including protection of genomic stability, antagonism of polycomb group (PcG) complexes to maintain a permissive transcriptional landscape in physiological conditions, and regulation of genes that mediate responses to immune checkpoint inhibitor therapy. An improved understanding of these mechanisms will bring new insights regarding cellular drivers of ccRCC growth and therapy response and, ultimately, will support the development of novel translational therapeutics.

RNA mis-splicing drives viral mimicry response after DNMTi therapy in SETD2-mutant kidney cancer

Tumors with mutations in chromatin regulators present attractive targets for DNA hypomethylating agent 5-aza-2'-deoxycytidine (DAC) therapy, which further disrupts cancer cells' epigenomic fidelity and reactivates transposable element (TE) expression to drive viral mimicry responses. SETD2 encodes a histone methyltransferase (H3K36me3) and is prevalently mutated in advanced kidney cancers. Here, we show that SETD2-mutant kidney cancer cells are especially sensitive in vitro and in vivo to DAC treatment. We find that the viral mimicry response are direct consequences of mis-splicing events, such as exon inclusions or extensions, triggered by DAC treatment in an SETD2-loss context. Comprehensive epigenomic analysis reveals H3K9me3 deposition, rather than DNA methylation dynamics, across intronic TEs might contribute to elevated mis-splicing rates. Through epigenomic and transcriptomic analyses, we show that SETD2-deficient kidney cancers are prone to mis-splicing, which can be therapeutically exacerbated with DAC treatment to increase viral mimicry activation and provide synergy with combinatorial immunotherapy approaches.

Histone methyltransferase SETD2: An epigenetic driver in clear cell renal cell carcinoma

SET domain-containing 2 (SETD2) is a lysine methyltransferase that catalyzes histone H3 lysine36 trimethylation (H3K36me3) and has been revealed to play important roles in the regulation of transcriptional elongation, RNA splicing, and DNA damage repair. SETD2 mutations have been documented in several cancers, including clear cell renal cell carcinoma (ccRCC). SETD2 deficiency is associated with cancer occurrence and progression by regulating autophagy flux, general metabolic activity, and replication fork speed. Therefore, SETD2 is considered a potential epigenetic therapeutic target and is the subject of ongoing research on cancer-related diagnosis and treatment. This review presents an overview of the molecular functions of SETD2 in H3K36me3 regulation and its relationship with ccRCC, providing a theoretical basis for subsequent antitumor therapy based on SETD2 or H3K36me3 targets.

SETD2 regulates the methylation of translation elongation factor eEF1A1 in clear cell renal cell carcinoma

BACKGROUND

SET domain-containing protein 2 (SETD2) is commonly mutated in renal cell carcinoma. SETD2 methylates histone H3 as well as a growing list of non-histone proteins.

OBJECTIVE

Initially, we sought to explore SETD2-dependent changes in lysine methylation of proteins in proximal renal tubule cells. Subsequently, we focused on changes in lysine methylation of the translation elongation factor eEF1A1.

METHODS

To accomplish these objectives, we initially performed a systems-wide analysis of protein lysine-methylation and expression in wild type (WT) and SETD2-knock out (KO) kidney cells and later focused our studies on eEF1A1 as well as the expression of lysine methyltransferases that regulate its lysine methylation.

RESULTS

We observed decreased lysine methylation of the translation elongation factor eEF1A1. EEF1AKMT2 and EEF1AKMT3 are known to methylate eEF1A1, and we show here that their expression is dependent on SET-domain function of SETD2. Globally, we observe differential expression of hundreds of proteins in WT versus SETD2-KO cells, including increased expression of many involved in protein translation. Finally, we observe decreased progression free survival and loss of EEF1AKMT2 gene expression in SETD2-mutated tumors predicted to have loss of function of the SET domain.

CONCLUSION

Overall, these data suggest that SETD2-mutated ccRCC, via loss of enzymatic function of the SET domain, displays dysregulation of protein translation as a potentially important component of the transformed phenotype.

The role of histone H3K36me3 writers, readers and erasers in maintaining genome stability

Histone Post-Translational Modifications (PTMs) play fundamental roles in mediating DNA-related processes such as transcription, replication and repair. The histone mark H3K36me3 and its associated methyltransferase SETD2 (Set2 in yeast) are archetypical in this regard, performing critical roles in each of these DNA transactions. Here, we present an overview of H3K36me3 regulation and the roles of its writers, readers and erasers in maintaining genome stability through facilitating DNA double-strand break (DSB) repair, checkpoint signalling and replication stress responses. Further, we consider how loss of SETD2 and H3K36me3, frequently observed in a number of different cancer types, can be specifically targeted in the clinic through exploiting loss of particular genome stability functions.

Mutational landscape and clinical outcome of pediatric acute myeloid leukemia with 11q23/KMT2A rearrangements

BACKGROUND

Alterations of 11q23/KMT2A are the most prevalent cytogenetic abnormalities in acute myeloid leukemia (AML) and the prognostic significance of 11q23/KMT2A-rearranged AML based on various translocation partners varies among different studies. However, few studies evaluated the molecular characteristics of 11q23/KMT2A-rearranged pediatric AML. We aim to analyze the mutational landscape of 11q23/KMT2A-rearranged AML and assess their prognostic value in outcomes.

METHODS

The mutational landscape and clinical prognosis of 105 children with 11q23/KMT2A-rearranged AML in comparison with 277 children with non-11q23/KMT2A-rearranged AML were analyzed using publicly accessible next-generation sequencing data from Therapeutically Applicable Research to Generate Effective Treatments (TARGET) dataset.

RESULTS

Pediatric AML patients with 11q23/KMT2A-rearrangements harbored a low number of mutations (Median, 1 mutation/patient, range, 1-22), 58% of which involved in RAS pathway mutations (KRAS, NRAS, and PTPN11) and 10.5% of which comprised of SETD2 mutations. Compared with non-11q23/KMT2A-rearranged AML, the incidence of KRAS (32.4% vs. 10.1%, P〈0.001) and SETD2 (10.5% vs. 1.4%, P=0.001) gene mutations in 11q23/KMT2A-rearranged AML was significantly higher. Both KRAS and SETD2 mutations occurred more often in t(10;11)(p12;q23). KRAS mutations were correlated with worse 5-year event-free survival [EFS] (Plog-rank = 0.001) and 5-year overall survival [OS] (Plog-rank = 0.009) and the presence of SETD2 mutations increases the 5-year relapse rate (PGray = 0.004). Multivariate analyses confirmed KRAS mutations in 11q23/KMT2A-rearranged AML as an independent predictor for poor EFS (hazard ratio [HR] = 2.10, P=0.05) and OS (HR = 2.39, P=0.054).

CONCLUSION

Our findings show that pediatric patients with 11q23/KMT2A rearrangements have characteristic mutation patterns and varying clinical outcomes depending on different translocation partners, which could be utilized to develop more accurate risk stratification and tailored therapies.

SETD2: from chromatin modifier to multipronged regulator of the genome and beyond

Histone modifying enzymes play critical roles in many key cellular processes and are appealing proteins for targeting by small molecules in disease. However, while the functions of histone modifying enzymes are often linked to epigenetic regulation of the genome, an emerging theme is that these enzymes often also act by non-catalytic and/or non-epigenetic mechanisms. SETD2 (Set2 in yeast) is best known for associating with the transcription machinery and methylating histone H3 on lysine 36 (H3K36) during transcription. This well-characterized molecular function of SETD2 plays a role in fine-tuning transcription, maintaining chromatin integrity, and mRNA processing. Here we give an overview of the various molecular functions and mechanisms of regulation of H3K36 methylation by Set2/SETD2. These fundamental insights are important to understand SETD2’s role in disease, most notably in cancer in which SETD2 is frequently inactivated. SETD2 also methylates non-histone substrates such as α-tubulin which may promote genome stability and contribute to the tumor-suppressor function of SETD2. Thus, to understand its role in disease, it is important to understand and dissect the multiple roles of SETD2 within the cell. In this review we discuss how histone methylation by Set2/SETD2 has led the way in connecting histone modifications in active regions of the genome to chromatin functions and how SETD2 is leading the way to showing that we also have to look beyond histones to truly understand the physiological role of an ‘epigenetic’ writer enzyme in normal cells and in disease.

High expression of SETDB1 mediated by miR-29a-3p associates with poor prognosis and immune invasion in breast invasive carcinoma

Background

Breast invasive carcinoma (BRCA) has a poor prognosis. Numerous studies have shown that SET domain bifurcated histone lysine methyltransferase 1 (SETDB1) is involved in the initiation and progression of many cancers. This study aims to reveal the potential mechanism of SETDB1 in the development and progression of BRCA.

Methods

The ONCOMINE database, TIMER database, UALCAN database and GEPIA database were used to analyze the expression of SETDB1 in human cancers. We evaluated the expression level of SETDB1 in cell lines by quantitative real-time polymerase chain reaction (qPCR), and the survival analysis of SETDB1 was performed on PrognoScan and Kaplan-Meier plotter websites. The upstream regulator was obtained from starBase database.

Results

We confirmed that SETDB1 messenger RNA (mRNA) level showed high expression in breast cell lines, and we also found that SETDB1 showed high expression in many types of cancers. Moreover, SETDB1 overexpression was positively correlated with poor prognosis in BRCA. Furthermore, we first predicted miR-29a-3p was a potential upstream regulator of SETDB1 in BRCA. Our findings indicated that SETDB1 might play a carcinogenic role by increasing the infiltration of immune cell and influencing immune checkpoint expression.

Conclusions

This study suggested that miR-29a-3p can mediate the expression of SETDB1 with poor prognosis and tumor immune infiltration in BRCA.

SETD2 loss perturbs the kidney cancer epigenetic landscape to promote metastasis and engenders actionable dependencies on histone chaperone complexes

SETD2 is a H3K36 trimethyltransferase that is mutated with high prevalence (13%) in clear cell renal cell carcinoma (ccRCC). Genomic profiling of primary ccRCC tumors reveals a positive correlation between SETD2 mutations and metastasis. However, whether and how SETD2 loss promotes metastasis remains unclear. In this study, we used SETD2-mutant metastatic ccRCC patient-derived cell line and xenograft models and showed that H3K36me3 restoration greatly reduced distant metastases of ccRCC in mice in an MMP1-dependent manner. An integrated multi-omics analysis using ATAC-seq, ChIP-seq, and RNA-seq established a tumor suppressor model in which loss of SETD2-mediated H3K36me3 activates enhancers to drive oncogenic transcriptional output through regulation of chromatin accessibility. Furthermore, we uncovered mechanism-based therapeutic strategies for SETD2-deficient cancer through the targeting of specific histone chaperone complexes including ASF1A/B and SPT16. Overall, SETD2 loss creates a permissive epigenetic landscape for cooperating oncogenic drivers to amplify transcriptional output, providing unique therapeutic opportunities.

The impact of gene-body H3K36me3 patterns on gene expression level changes in chronic myelogenous leukemia

Chronic myelogenous leukemia (CML) is a malignant clonal disease of hematopoietic stem cells. Researches have exhibited that the progression of CML is related to histone modifications. Here, we perform the systematic analyses of H3K36me3 patterns and gene expression level changes. We observe that the genes with higher gene-body H3K36me3 levels in normal cells show fewer expression changes during leukemogenesis, while the genes with lower gene-body H3K36me3 levels in normal cells yield obvious expression changes during leukemogenesis (ρ = -0.98, P = 9.30 × 10^-8). These findings are conserved in human lung/breast cancers and mouse CML, regardless of gene expression levels and gene lengths. Regulatory element analysis and Random Forest regression display that Hoxd13, Rara, Scl, Smad3, Smad4 and Tgif1 induce the up-regulation of genes with lower H3K36me3 levels (ρ = 0.97, P = 2.35 × 10^-56). Enrichment analysis shows that the differentially expressed genes with lower H3K36me3 levels are involved in leukemia-related pathways, such as leukocyte migration and regulation of leukocyte activation. Finally, six driver genes (Tp53, Wt1, Dnmt3a, Cacna1b, Phactr1 and Gbp4) with lower H3K36me3 levels are identified. Our analyses indicate that lower gene-body H3K36me3 levels may serve as a biomarker for the progression of CML.

The methyltransferase SETD2 couples transcription and splicing by engaging mRNA processing factors through its SHI domain

Heterogeneous ribonucleoproteins (hnRNPs) are RNA binding molecules that are involved in key processes such as RNA splicing and transcription. One such hnRNP protein, hnRNP L, regulates alternative splicing (AS) by binding to pre-mRNA transcripts. However, it is unclear what factors contribute to hnRNP L-regulated AS events. Using proteomic approaches, we identified several key factors that co-purify with hnRNP L. We demonstrate that one such factor, the histone methyltransferase SETD2, specifically interacts with hnRNP L in vitro and in vivo. This interaction occurs through a previously uncharacterized domain in SETD2, the SETD2-hnRNP Interaction (SHI) domain, the deletion of which, leads to a reduced H3K36me3 deposition. Functionally, SETD2 regulates a subset of hnRNP L-targeted AS events. Our findings demonstrate that SETD2, by interacting with Pol II as well as hnRNP L, can mediate the crosstalk between the transcription and the splicing machinery.

The methylation of Histone 3 at Lysine 36 (H3K36) has been implicated in the regulation of transcription and coupled processes such as mRNA splicing. Here the authors show that the histone methyltransferase SETD2 interacts with hnRNP L to mediate the crosstalk between the transcription and splicing machineries.

Monomorphic Epitheliotropic Intestinal T-Cell Lymphoma in Asia Frequently Shows SETD2 Alterations

Simple Summary

Monomorphic epitheliotropic intestinal T-cell lymphoma (MEITL) is a rare primary T-cell lymphoma of the digestive tract that is characterized by an aggressive clinical course. The aim of this study was to analyze the clinicopathological characteristics and genomic profile of Asian MEITL. In this study, nine cases of Japanese MEITL were analyzed by targeted Next Generation Sequencing and immunohistochemistry and were integrated with previously reported whole-genome copy number microarray-based assay data. All cases showed alterations of the tumor suppressor gene SETD2 and mutations in one or more genes of the JAK/STAT pathway. Therefore, we concluded that the combination of epigenetic deregulation and cell signaling activation may represent a major oncogenic event in the pathogenesis of Asian MEITL, similar to Western MEITL.

Abstract

Monomorphic epitheliotropic intestinal T-cell lymphoma (MEITL) is a rare primary T-cell lymphoma of the digestive tract derived from intraepithelial lymphocytes and characterized by an aggressive clinical course. In this study, nine cases of Japanese MEITL were analyzed by targeted Next Generation Sequencing (NGS) and immunohistochemistry and were integrated with previously reported whole-genome copy number microarray-based assay data. The highlight of our findings is that all cases showed alterations of the tumor suppressor gene SETD2 by mutations and/or loss of the corresponding 3p21 locus. We also demonstrated that all cases showed mutations in one or more genes of JAK/STAT pathway. Therefore, the combination of epigenetic deregulation and cell signaling activation represent major oncogenic events in the pathogenesis of MEITL in Asian MEITL, similar to Western MEITL.

Assessment of hazard immune-related genes and tumor immune infiltrations in renal cell carcinoma

BACKGROUND

The present study aimed to explore and validate a prognostic immune signature, to formulate a prognosis for ccRCC patients combined with immune-infiltration analysis.

METHODS

Public datasets were used as our source of multi-omics data. Differential analysis was performed via the edgeR package. A prognostic immune signature was identified by univariate Cox analysis, and we constructed an integrative tumor-associated immune genes (TAIG) model from the multivariate Cox results. In order to interrogate and identify the related crosstalk, functional analysis was deployed. Significantly, we implemented the CIBERSORT algorithm to estimate the immune cell fractions in the ccRCC samples, and analyzed the differential abundance of tumor-infiltrating immune cells in two TAIG groups, using a Wilcoxon rank-sum test. The prognostic role of differential immune cells was further assessed via a Kaplan-Meier analysis. In addition, we investigated the associations of a single immune signature with specific immune cells.

RESULTS

A total of 628 ccRCC patients were comprised in our integrative analysis, including 537 ccRCC patients in the discovery group and 91 patients in the validation group. Fourteen key immune signatures were subsequently identified. A figure of 0.802 was registered for AUC, and worse prognosis was evinced for those patients with a higher TAIG. Correlation analysis indicated that TAIG correlated closely with both clinical variables and TMB. Moreover, functional analysis implicated the immune-related GO items or crosstalk. Hence, we were able to identify the relationships obtaining between tumor-infiltrating immune cells and TAIG. The differential abundance of immune cells showed a significant prognostic difference and consisted of memory-activated CD4⁺ T cells, T follicular helper cells, T regulatory cells, and so on. Moreover, we also characterized the associations between identified signatures and specific immune cells. Finally, the five-year AUC in the ICGC cohort was 0.72, suggesting the robustness of the TAIG that we constructed.

CONCLUSIONS

Overall, our team characterized the tumor-associated immune signature in ccRCC, and further identified the prognostic tumor-infiltrating immune cells related to TAIG. This in turn provided a solid foundation for investigating individualized immunotherapy, as well as other relevant mechanisms.

Regulation of SETD2 stability is important for the fidelity of H3K36me3 deposition

Background

The histone H3K36me3 mark regulates transcription elongation, pre-mRNA splicing, DNA methylation, and DNA damage repair. However, knowledge of the regulation of the enzyme SETD2, which deposits this functionally important mark, is very limited.

Results

Here, we show that the poorly characterized N-terminal region of SETD2 plays a determining role in regulating the stability of SETD2. This stretch of 1–1403 amino acids contributes to the robust degradation of SETD2 by the proteasome. Besides, the SETD2 protein is aggregate prone and forms insoluble bodies in nuclei especially upon proteasome inhibition. Removal of the N-terminal segment results in the stabilization of SETD2 and leads to a marked increase in global H3K36me3 which, uncharacteristically, happens in a Pol II-independent manner.

Conclusion

The functionally uncharacterized N-terminal segment of SETD2 regulates its half-life to maintain the requisite cellular amount of the protein. The absence of SETD2 proteolysis results in a Pol II-independent H3K36me3 deposition and protein aggregation.

Genome-wide alterations of uracil distribution patterns in human DNA upon chemotherapeutic treatments

Numerous anti-cancer drugs perturb thymidylate biosynthesis and lead to genomic uracil incorporation contributing to their antiproliferative effect. Still, it is not yet characterized if uracil incorporations have any positional preference. Here, we aimed to uncover genome-wide alterations in uracil pattern upon drug treatments in human cancer cell line models derived from HCT116. We developed a straightforward U-DNA sequencing method (U-DNA-Seq) that was combined with in situ super-resolution imaging. Using a novel robust analysis pipeline, we found broad regions with elevated probability of uracil occurrence both in treated and non-treated cells. Correlation with chromatin markers and other genomic features shows that non-treated cells possess uracil in the late replicating constitutive heterochromatic regions, while drug treatment induced a shift of incorporated uracil towards segments that are normally more active/functional. Data were corroborated by colocalization studies via dSTORM microscopy. This approach can be applied to study the dynamic spatio-temporal nature of genomic uracil.

Transcription of intragenic CpG islands influences spatiotemporal host gene pre-mRNA processing

Alternative splicing (AS) and alternative polyadenylation (APA) generate diverse transcripts in mammalian genomes during development and differentiation. Epigenetic marks such as trimethylation of histone H3 lysine 36 (H3K36me3) and DNA methylation play a role in generating transcriptome diversity. Intragenic CpG islands (iCGIs) and their corresponding host genes exhibit dynamic epigenetic and gene expression patterns during development and between different tissues. We hypothesise that iCGI-associated H3K36me3, DNA methylation and transcription can influence host gene AS and/or APA. We investigate H3K36me3 and find that this histone mark is not a major regulator of AS or APA in our model system. Genomewide, we identify over 4000 host genes that harbour an iCGI in the mammalian genome, including both previously annotated and novel iCGI/host gene pairs. The transcriptional activity of these iCGIs is tissue- and developmental stage-specific and, for the first time, we demonstrate that the premature termination of host gene transcripts upstream of iCGIs is closely correlated with the level of iCGI transcription in a DNA-methylation independent manner. These studies suggest that iCGI transcription, rather than H3K36me3 or DNA methylation, interfere with host gene transcription and pre-mRNA processing genomewide and contributes to the spatiotemporal diversification of both the transcriptome and proteome.

Transcription and splicing: A two-way street

RNA synthesis by RNA polymerase II and RNA processing are closely coupled during the transcription cycle of protein-coding genes. This coupling affords opportunities for quality control and regulation of gene expression and the effects can go in both directions. For example, polymerase speed can affect splice site selection and splicing can increase transcription and affect the chromatin landscape. Here we review the many ways that transcription and splicing influence one another, including how splicing "talks back" to transcription. We will also place the connections between transcription and splicing in the context of other RNA processing events that define the exons that will make up the final mRNA. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Alternative Lengthening of Telomeres and Differential Expression of Endocrine Transcription Factors Distinguish Metastatic and Non-metastatic Insulinomas

Insulin-producing pancreatic neuroendocrine tumors (PanNETs)/insulinomas are generally considered to be indolent tumors with an excellent prognosis after complete resection. However, some insulinomas have a poor prognosis due to relapses and metastatic disease. Recently, studies in non-functional PanNETs indicated that behavior can be stratified according to alpha- and beta-cell differentiation, as defined by expression of the transcription factors ARX and PDX1, respectively. It is unknown whether similar mechanisms play a role in insulinomas. Therefore, we determined ARX and PDX1 expression in a cohort of 35 sporadic primary insulinomas and two liver metastases of inoperable primary insulinomas. In addition, WHO grade and loss of ATRX or DAXX were determined by immunohistochemistry, and alternative lengthening of telomeres (ALT) and CDKN2A status by fluorescence in situ hybridization. These findings were correlated with tumor characteristics and clinical follow-up data. In total, five out of 37 insulinoma patients developed metastatic disease. Metastatic insulinomas were all larger than 3 cm, whereas the indolent insulinomas were smaller (p value < 0.05). All three primary insulinomas that metastasized showed ARX expression, 2/3 showed ALT, and 1/3 had a homozygous deletion of CDKN2A as opposed to absence of ARX expression, ALT, or CDKN2A deletions in the 32 non-metastatic cases. The two liver metastases also showed ARX expression and ALT (2/2). The presence of ARX expression, which is usually absent in beta-cells, and genetic alterations not seen in indolent insulinomas strongly suggest a distinct tumorigenic mechanism in malignant insulinomas, with similarities to non-functional PanNETs. These observations may inform future follow-up strategies after insulinoma surgery.

Histone methyltransferase SETD2: a potential tumor suppressor in solid cancers

Epigenetic regulation plays an important role in the occurrence, development and treatment of malignant tumors; and a great deal of attention has been paid to the histone methylation level in recent years. As a 230-kD epigenetic regulator, the histone H3 lysine 36 histone (H3K36) methyltransferase SETD2 is a key enzyme of the nuclear receptor SET domain-containing (NSD) family, which is associated with a specific hyperphosphorylated domain, a large subunit of RNA polymerase II (RNAPII), named RNAPII subunit B1 (RPB1), and SETD2 which methylates the ly-36 position of dimethylated histone H3 (H3K36me2) to generate trimethylated H3K36 (H3K36me3). SETD2 is involved in various cellular processes, including transcriptional regulation, DNA damage repair, non-histone protein-related functions and some other processes. Great efforts of high-throughput sequencing have revealed that SETD2 is mutated or its function is lost in a range of solid cancers, including renal cancer, gastrointestinal cancer, lung cancer, pancreatic cancer, osteosarcoma, and so on. Mutation, or functional loss, of the SETD2 gene produces dysfunction in corresponding tumor tissue proteins, leading to tumorigenesis, progression, chemotherapy resistance, and unfavorable prognosis, suggesting that SETD2 possibly acts as a tumor suppressor. However, its underlying mechanism remains largely unexplored. In the present study, we summarized the latest advances of effects of SETD2 expression at the mRNA and protein levels in solid cancers, and its potential molecular and cellular functions as well as clinical applications were also reviewed.

H3K36me3, message from chromatin to DNA damage repair

Histone marks control many cellular processes including DNA damage repair. This review will focus primarily on the active histone mark H3K36me3 in the regulation of DNA damage repair and the maintenance of genomic stability after DNA damage. There are diverse clues showing H3K36me3 participates in DNA damage response by directly recruiting DNA repair machinery to set the chromatin at a “ready” status, leading to a quick response upon damage. Reduced H3K36me3 is associated with low DNA repair efficiency. This review will also place a main emphasis on the H3K36me3-mediated DNA damage repair in the tumorigenesis of the newly found oncohistone mutant tumors. Gaining an understanding of different aspects of H3K36me3 in DNA damage repair, especially in cancers, would share the knowledge of chromatin and DNA repair to serve to the drug discovery and patient care.

Inferring gene expression networks with hubs using a degree weighted Lasso approach

MOTIVATION

Genome-scale gene networks contain regulatory genes called hubs that have many interaction partners. These genes usually play an essential role in gene regulation and cellular processes. Despite recent advancements in high-throughput technology, inferring gene networks with hub genes from high-dimensional data still remains a challenging problem. Novel statistical network inference methods are needed for efficient and accurate reconstruction of hub networks from high-dimensional data.

RESULTS

To address this challenge we propose DW-Lasso, a degree weighted Lasso (least absolute shrinkage and selection operator) method which infers gene networks with hubs efficiently under the low sample size setting. Our network reconstruction approach is formulated as a two stage procedure: first, the degree of networks is estimated iteratively, and second, the gene regulatory network is reconstructed using degree information. A useful property of the proposed method is that it naturally favors the accumulation of neighbors around hub genes and thereby helps in accurate modeling of the high-throughput data under the assumption that the underlying network exhibits hub structure. In a simulation study, we demonstrate good predictive performance of the proposed method in comparison to traditional Lasso type methods in inferring hub and scale-free graphs. We show the effectiveness of our method in an application to microarray data of Escherichia coli and RNA sequencing data of Kidney Clear Cell Carcinoma from The Cancer Genome Atlas datasets.

AVAILABILITY AND IMPLEMENTATION

Under the GNU General Public Licence at https://cran.r-project.org/package=DWLasso.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Multi-omics analysis reveals epithelial-mesenchymal transition-related gene FOXM1 as a novel prognostic biomarker in clear cell renal carcinoma

Identification of novel clinical biomarker in clear cell renal carcinoma (ccRCC) is warranted. Integrating transcriptome (n=1669), DNA methylation (n=577) and copy number data (n=832), we developed a method to identify driver biomarkers by analyzing the omics-level dynamics of Epithelial-Mesenchymal Transition (EMT)-related genes in ccRCC. We first identified 504 expression dynamic changed genes involved in ccRCC-associated key pathways such as EMT, cell cycle, EGFR and PI3K/AKT signaling. Further analysis identified 229 (90 gene promoters) aberrant expression quantitative trait methylation (eQTM) and 256 genes with expression quantitative trait copy number (eQTCN) alterations. Among them, FOXM1 was affected by both eQTM and eQTCN. FOXM1 copy number amplification (115/500, 23% of patients), occurred in an amplified peak in chromosome 12q13.3, was enriched in late-stage ccRCC samples and was associated with worse survival. FOXM1-overexpressed pT3 patients with distant metastasis showed ~25% shorter overall survival in both training (log-rank P=0.006) and validation (log-rank P=0.018) cohorts. The eQTM-gene hybrid signature (cg00044170 and FOXM1), superior to either gene expression or DNA methylation alone, showed great potential in diagnosing localized ccRCC in training (area under curve = 0.958) and validation datasets. FOXM1 could be a novel prognostic biomarker and shed light for early diagnosis at molecular level in ccRCC.

Loss of SETD2 Induces a Metabolic Switch in Renal Cell Carcinoma Cell Lines toward Enhanced Oxidative Phosphorylation

SETD2, a histone H3 lysine trimethyltransferase, is frequently inactivated and associated with recurrence of clear cell renal cell carcinoma (ccRCC). However, the impact of SETD2 loss on metabolic alterations in ccRCC is still unclear. In this study, SETD2 null isogenic 38E/38F clones derived from 786-O cells were generated by zinc finger nucleases, and subsequent metabolic, genomic, and cellular phenotypic changes were analyzed by targeted metabolomics, RNA sequencing, and biological methods, respectively. Our results showed that compared with parental 786-O cells, 38E/38F cells had elevated levels of MTT/Alamar blue levels, ATP, glycolytic/mitochondrial respiratory capacity, citrate synthase (CS) activity, and TCA metabolites such as aspartate, malate, succinate, fumarate, and α-ketoglutarate. The 38E/38F cells also utilized alternative sources beyond pyruvate to generate acetyl-CoA for the TCA cycle. Moreover, 38E/38F cells showed disturbed gene networks mainly related to mitochondrial metabolism and the oxidation of fatty acids and glucose, which was associated with increased PGC1α, mitochondrial mass, and cellular size/complexity. Our results indicate that SETD2 deficiency induces a metabolic switch toward enhanced oxidative phosphorylation in ccRCC, which can be related to PGC1α-mediated metabolic networks. Therefore, this current study lays the foundation for the further development of a global metabolic analysis of cancer cells in individual patients, which ultimately will have significant potential for the discovery of novel therapeutics and precision medicine in SETD2-inactivated ccRCC.

SETD2 mutations in primary central nervous system tumors

Mutations in SETD2 are found in many tumors, including central nervous system (CNS) tumors. Previous work has shown these mutations occur specifically in high grade gliomas of the cerebral hemispheres in pediatric and young adult patients. We investigated SETD2 mutations in a cohort of approximately 640 CNS tumors via next generation sequencing; 23 mutations were detected across 19 primary CNS tumors. Mutations were found in a wide variety of tumors and locations at a broad range of allele frequencies. SETD2 mutations were seen in both low and high grade gliomas as well as non-glial tumors, and occurred in patients greater than 55 years of age, in addition to pediatric and young adult patients. High grade gliomas at first occurrence demonstrated either frameshift/truncating mutations or point mutations at high allele frequencies, whereas recurrent high grade gliomas frequently harbored subclones with point mutations in SETD2 at lower allele frequencies in the setting of higher mutational burdens. Comparison with the TCGA dataset demonstrated consistent findings. Finally, immunohistochemistry showed decreased staining for H3K36me3 in our cohort of SETD2 mutant tumors compared to wildtype controls. Our data further describe the spectrum of tumors in which SETD2 mutations are found and provide a context for interpretation of these mutations in the clinical setting.

Loss of Chromatin-Remodeling Proteins and/or CDKN2A Associates With Metastasis of Pancreatic Neuroendocrine Tumors and Reduced Patient Survival Times

Despite prognostic grading and staging systems, it is a challenge to predict outcomes for patients with pancreatic neuroendocrine tumors (PanNETs). Sequencing studies of PanNETs have identified alterations in death domain-associated protein (DAXX) and alpha-thalassemia/mental retardation X-linked chromatin remodeler (ATRX). In tumors, mutations in DAXX or ATRX and corresponding loss of protein expression correlate with shorter times of disease-free survival and disease-specific survival of patients. However, DAXX or ATRX proteins were lost in only 50% of distant metastases analyzed. We performed whole-exome sequencing analyses of 20 distant metastases from 20 patients with a single nonsyndrome, nonfunctional PanNET. We found distant metastases contained alterations in multiple endocrine neoplasia type 1 (MEN1) (n = 8), ATRX (n = 5), DAXX (n = 5), TSC2 (n = 3), and DEP domain containing 5 (DEPDC5) (n = 3). We found copy number loss of cyclin dependent kinase inhibitor 2A (CDKN2A) in 15 metastases (75%) and alterations in genes that regulate chromatin remodeling, including set domain containing 2 (SETD2) (n = 4), AT-rich interaction domain 1A (ARID1A) (n = 2), chromodomain helicase DNA binding protein 8 (CHD8) (n = 2), and DNA methyl transferase 1 (DNMT1) (n = 2). In a separate analysis of 347 primary PanNETs, we found loss or deletion of DAXX and ATRX, disruption of SETD2 function (based on loss of H3 lysine 36 trimethylation), loss of ARID1A expression or deletions in CDKN2A in 81% of primary PanNETs with distant metastases. Among patients with loss or deletion of at least 1 of these proteins or genes, 39% survived disease-free for 5 years and 44% had disease-specific survival times of 10 years. Among patients without any of these alterations, 98% survived disease-free for 5 years and 95% had disease-specific survival times of 10 years. Therefore, primary PanNETs with loss of DAXX, ATRX, H3 lysine 36 trimethylation, ARID1A, and/or CDKN2A associate with shorter survival times of patients. Our findings indicate that alterations in chromatin-remodeling genes and CDKN2A contribute to metastasis of PanNETs.

Global isoform-specific transcript alterations and deregulated networks in clear cell renal cell carcinoma

Extensive genome-wide analyses of deregulated gene expression have now been performed for many types of cancer. However, most studies have focused on deregulation at the gene-level, which may overlook the alterations of specific transcripts for a given gene. Clear cell renal cell carcinoma (ccRCC) is one of the best-characterized and most pervasive renal cancers, and ccRCCs are well-documented to have aberrant RNA processing. In the present study, we examine the extent of aberrant isoform-specific RNA expression by reporting a comprehensive transcript-level analysis, using the new kallisto-sleuth-RATs pipeline, investigating coding and non-coding differential transcript expression in ccRCC. We analyzed 50 ccRCC tumors and their matched normal samples from The Cancer Genome Altas datasets. We identified 7,339 differentially expressed transcripts and 94 genes exhibiting differential transcript isoform usage in ccRCC. Additionally, transcript-level coexpression network analyses identified vasculature development and the tricarboxylic acid cycle as the most significantly deregulated networks correlating with ccRCC progression. These analyses uncovered several uncharacterized transcripts, including lncRNAs FGD5-AS1 and AL035661.1, as potential regulators of the tricarboxylic acid cycle associated with ccRCC progression. As ccRCC still presents treatment challenges, our results provide a new resource of potential therapeutics targets and highlight the importance of exploring alternative methodologies in transcriptome-wide studies.

SETD2 Haploinsufficiency for Microtubule Methylation Is an Early Driver of Genomic Instability in Renal Cell Carcinoma

Loss of the short arm of chromosome 3 (3p) occurs early in >95% of clear cell renal cell carcinoma (ccRCC). Nearly ubiquitous 3p loss in ccRCC suggests haploinsufficiency for 3p tumor suppressors as early drivers of tumorigenesis. We previously reported methyltransferase SETD2, which trimethylates H3 histones on lysine 36 (H3K36me3) and is located in the 3p deletion, to also trimethylate microtubules on lysine 40 (αTubK40me3) during mitosis, with αTubK40me3 required for genomic stability. We now show that monoallelic, Setd2-deficient cells retaining H3K36me3, but not αTubK40me3, exhibit a dramatic increase in mitotic defects and micronuclei count, with increased viability compared with biallelic loss. In SETD2-inactivated human kidney cells, rescue with a pathogenic SETD2 mutant deficient for microtubule (αTubK40me3), but not histone (H3K36me3) methylation, replicated this phenotype. Genomic instability (micronuclei) was also a hallmark of patient-derived cells from ccRCC. These data show that the SETD2 tumor suppressor displays a haploinsufficiency phenotype disproportionately impacting microtubule methylation and serves as an early driver of genomic instability.Significance: Loss of a single allele of a chromatin modifier plays a role in promoting oncogenesis, underscoring the growing relevance of tumor suppressor haploinsufficiency in tumorigenesis. Cancer Res; 78(12); 3135-46. ©2018 AACR.

HistoneH3 demethylase JMJD2A promotes growth of liver cancer cells through up-regulating miR372

Changes in histone lysine methylation status have been observed during cancer formation. JMJD2A protein is a demethylase that is overexpressed in several tumors. Herein, our results demonstrate that JMJD2A accelerates malignant progression of liver cancer cells in vitro and in vivo. Mechanistically, JMJD2A promoted the expression and mature of pre-miR372 epigenetically. Notably, miR372 blocks the editing of 13th exon-introns-14th exon and forms a novel transcript(JMJD2AΔ) of JMJD2A. In particular, JMJD2A inhibited P21(WAF1/Cip1) expression by decreasing H3K9me3 dependent on JMJD2AΔ. Thereby, JMJD2A could enhance Pim1 transcription by suppressing P21(WAF1/Cip1). Furthermore, through increasing the expression of Pim1, JMJD2A could facilitate the interaction among pRB, CDK2 and CyclinE which prompts the transcription and translation of oncogenic C-myc. Strikingly, JMJD2A may trigger the demethylation of Pim1. On the other hand, Pim1 knockdown and P21(WAF1/Cip1) overexpression fully abrogated the oncogenic function of JMJD2A. Our observations suggest that JMJD2A promotes liver cancer cell cycle progress through JMJD2A-miR372-JMJD2AΔ-P21WAF1/Cip1-Pim1-pRB-CDK2-CyclinE-C-myc axis. This study elucidates a novel mechanism for JMJD2A in liver cancer cells and suggests that JMJD2A can be used as a novel therapeutic targets of liver cancer.

SETD2: an epigenetic modifier with tumor suppressor functionality

In the past decade important progress has been made in our understanding of the epigenetic regulatory machinery. It has become clear that genetic aberrations in multiple epigenetic modifier proteins are associated with various types of cancer. Moreover, targeting the epigenome has emerged as a novel tool to treat cancer patients. Recently, the first drugs have been reported that specifically target SETD2-negative tumors. In this review we discuss the studies on the associated protein, Set domain containing 2 (SETD2), a histone modifier for which mutations have only recently been associated with cancer development. Our review starts with the structural characteristics of SETD2 and extends to its corresponding function by combining studies on SETD2 function in yeast, Drosophila, Caenorhabditis elegans, mice, and humans. SETD2 is now generally known as the single human gene responsible for trimethylation of lysine 36 of Histone H3 (H3K36). H3K36me3 readers that recruit protein complexes to carry out specific processes, including transcription elongation, RNA processing, and DNA repair, determine the impact of this histone modification. Finally, we describe the prevalence of SETD2-inactivating mutations in cancer, with the highest frequency in clear cell Renal Cell Cancer, and explore how SETD2-inactivation might contribute to tumor development.

Whole-Exome Sequencing in Two Extreme Phenotypes of Response to VEGF-Targeted Therapies in Patients With Metastatic Clear Cell Renal Cell Carcinoma

Advances in next-generation sequencing have provided a unique opportunity to understand the biology of disease and mechanisms of sensitivity or resistance to specific agents. Renal cell carcinoma (RCC) is a heterogeneous disease and highly variable clinical responses have been observed with vascular endothelial growth factor (VEGF)-targeted therapy (VEGF-TT). We hypothesized that whole-exome sequencing analysis might identify genotypes associated with extreme response or resistance to VEGF-TT in metastatic (mRCC). Patients with mRCC who had received first-line sunitinib or pazopanib and were in 2 extreme phenotypes of response were identified. Extreme responders (ERs) were defined as those with partial response or complete response for 3 or more years (n=13) and primary refractory patients (PRPs) were defined as those with progressive disease within the first 3 months of therapy (n=14). International Metastatic RCC Database Consortium prognostic scores were not significantly different between the groups (P=.67). Considering the genes known to be mutated in RCC at significant frequency, PBRM1 mutations were identified in 7 ERs (54%) versus 1 PRP (7%) (P=.01). In addition, mutations in TP53 (n=4) were found only in PRPs (P=.09). Our data suggest that mutations in some genes in RCC may impact response to VEGF-TT.

Multicenter Validation of Enhancer of Zeste Homolog 2 Expression as an Independent Prognostic Marker in Localized Clear Cell Renal Cell Carcinoma

Purpose

Enhancer of zeste homolog 2 (EZH2), a chromatin remodeler, is implicated in the pathogenesis of clear cell renal cell carcinoma (ccRCC). However, the effect of EZH2 on outcomes in localized ccRCC is unclear, and molecular biomarkers are not currently integrated into prognostic models or adjuvant therapy trials.

Methods

We performed Cox regression to evaluate the association of tumor-based EZH2 gene and protein expression with survival in three independent cohorts: a cohort from The Cancer Genome Atlas (n = 532), a cohort from University of Texas Southwestern Medical Center (n = 122), and a cohort from Mayo Clinic (n = 1,338). Analyses were adjusted for the prognostic stage, size, grade, and necrosis (SSIGN) score as well as within low-, intermediate-, and high-risk SSIGN groups.

Results

Patients in The Cancer Genome Atlas cohort with EZH2-high gene expression were 1.5 times more likely to experience overall death than patients with EZH2-low expression (95% CI, 1.1 to 2.3; P = .028). Patients in the University of Texas Southwestern Medical Center cohort with EZH2-high protein expression were two times more likely to experience overall death than patients with EZH2-low expression (95% CI, 1.1 to 4.4; P = .034). Similarly, patients in the Mayo Clinic cohort with EZH2-high protein expression were 1.4 times more likely to experience overall death (95% CI, 1.2 to 1.7; P < .001). Patients in the Mayo Clinic cohort with EZH2-high protein expression were nearly two times more likely to experience RCC-specific death (95% CI, 1.5 to 2.6; P < .001); EZH2 protein expression was particularly prognostic among patients with low-risk SSIGN tumors (HR, 6.1; 95% CI, 3.4 to 11.1; P < .001).

Conclusion

EZH2 expression accurately predicts risk of RCC death beyond existing clinicopathologic models, particularly in low- and intermediate-risk SSIGN tumors. Further studies are required to incorporate molecular biomarkers into surveillance guidelines and adjuvant clinical trials.

Integration of Recurrent Somatic Mutations with Clinical Outcomes: A Pooled Analysis of 1049 Patients with Clear Cell Renal Cell Carcinoma

BACKGROUND

Analyses of associations between clinicopathologic outcomes and recurrent somatic mutations in clear cell renal cell carcinoma (ccRCC) have been limited to individual cohorts.

OBJECTIVE

To define clinicopathologic associations between specific mutations and ccRCC disease characteristics.

DESIGN, SETTING, AND PARTICIPANTS

DNA sequencing data were pooled from three collaborative genomic cohorts (n=754) and our institutional database (n=295). All patients had clinical data and identification of somatic mutations from their primary tumors.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Analysis of gene mutations for associations with maximal tumor size (linear regression) and pathologic stage (logistic regression). Cancer-specific survival (CSS) and recurrence-free survival (RFS) were calculated using competing risks methods. Analyses were adjusted for cohort site, and results were adjusted for multiple testing (q value). Relevant genes were used in multivariable models that included confounding variables and the validated Mayo Clinic Stage, Size, Grade, and Necrosis (SSIGN) score.

RESULTS AND LIMITATIONS

Association with tumor size was found for mutations in BAP1 (q=0.013). No mutations were found to be associated with stage after adjusted analysis. Mutations in BAP1 (q=0.004) and TP53 (q=0.001) were associated with decreased CSS in a multivariable model; only TP53 (q=0.005) remained significant when SSIGN score was included. SETD2 mutations (q=0.047) were associated with decreased RFS in multivariable models, including models with SSIGN score.

CONCLUSIONS

In >1000 patients with ccRCC, pooled analysis and multivariable modeling demonstrated that three mutated genes have statistically significant associations with poor clinical outcomes. This included the more commonly mutated BAP1 and SETD2 and the less frequently mutated TP53. After adjustment for clinical confounders, mutations of TP53 and SETD2 were associated with decreased CSS and RFS, respectively.

PATIENT SUMMARY

Using rigorous statistical methods, this study affirmed that certain mutations in clear cell renal cell carcinoma may portend inferior survival and an increased risk of recurrence.