Epigenetic genes regulated by the BRAFV600E signaling are associated with alterations in the methylation and expression of tumor suppressor genes and patient survival in melanoma

Mutations in KMT2C, BCOR and KDM5C Predict Response to Immune Checkpoint Blockade Therapy in Non-Small Cell Lung Cancer

Simple Summary

Efficient biomarkers are urgently needed to predict response to immune checkpoint blockade (ICB) therapy for non-small cell lung cancer (NSCLC), particularly NSCLC with low tumor mutational burden (TMB). Here, we show that mutations of three chromatin remodeling-related genes, including KMT2C, BCOR and KDM5C, are associated with the ICB response in NSCLC, including NSCLC with low TMB level. Furthermore, this association is further improved by a combined use of KMT2C/BCOR/KDM5C mutations with TMB or PD-L1 expression. These data suggest that KMT2C/BCOR/KDM5C mutation status has the potential to serve as a predictive biomarker for ICB therapy in NSCLC.

Abstract

Efficient predictive biomarkers are urgently needed to identify non-small cell lung cancer (NSCLC) patients who could benefit from immune checkpoint blockade (ICB) therapy. Since chromatin remodeling is required for DNA repair process, we asked whether mutations in chromatin remodeling genes could increase tumor mutational burden (TMB) and predict response to ICB therapy in NSCLC. Analysis of seven ICB-treated NSCLC cohorts revealed that mutations of three chromatin remodeling-related genes, including KMT2C, BCOR and KDM5C, were significantly associated with ICB response, and combined mutations of these three genes further enhance this association. NSCLC patients with KMT2C/BCOR/KDM5C mutations had comparable clinical outcomes to TMB-high patients in terms of objective response rate, durable clinical benefit and overall survival. Although KMT2C/BCOR/KDM5C mutations were positively correlated with TMB levels in NSCLC, the association of this mutation with better ICB response was independent of tumor TMB and programmed death-ligand 1 (PD-L1) level, and combination of KMT2C/BCOR/KDM5C mutations with TMB or PD-L1 further improve the prediction of ICB response in NSCLC patients. Cancer Genome Atlas (TCGA) pan-cancer analysis suggested that the association of KMT2C/BCOR/KDM5C mutations with ICB response observed here might not result from DNA repair defects. In conclusion, our data indicate that KMT2C/BCOR/KDM5C mutation has the potential to serve as a predictive biomarker, alone or combined with PD-L1 expression or TMB, for ICB therapy in NSCLC.

Identification of gene-specific DNA methylation signature for Colorectal Cancer

BACKGROUND

Colorectal Cancer (CC), a common disease causing approximately million deaths annually, has been the third most frequent type of malignancy. We aimed to identify gene-specific DNA methylation signature to function as prognostic and predictive markers for CC patient survival.

METHODS

Expression profiles of gene-specific DNA methylation and the corresponding clinical information of 201 CC patients were downloaded from The Cancer Genome Atlas (TCGA) dataset and differentially expressed gene-specific DNA methylation was identified after tumor subtype classification. A risk score model was further built by analyzing the expression data of these gene-specific DNA methylations from the training dataset of CC patients.

RESULTS

Totally, 214 gene-specific DNA methylations were found to be expressed significantly between different subtypes of CC, including 150 up-regulated and 64 down-regulated ones. Up-regulated gene-specific DNA methylation accounted for 70.1% and the down-regulated gene-specific DNA methylation accounted for 29.9%. Hereinto, six gene-specific DNA methylations were obtained, including methy_vimentin and methy_ TFPI2, which were found significantly correlated with overall survival status of patients with CC.

CONCLUSIONS

With the six gene-specific DNA methylation signatures, patients in the training set were divided into low-risk and high- risk groups. What's more, gene-specific DNA methylation target genes were highly associated with protein phosphorylation, which indicated that further research on phosphorylation of target gene-coding protein might provide new sight on the treatment of CC.

REC8 is a novel tumor suppressor gene epigenetically robustly targeted by the PI3K pathway in thyroid cancer

The role of the PI3K pathway in human cancer has been well established, but much of its molecular mechanism, particularly the epigenetic aspect, remains to be defined. We hypothesized that aberrant methylation and hence altered expression of certain unknown important genes induced by the genetically activated PI3K pathway signaling is a major epigenetic mechanism in human tumorigenesis. Through a genome-wide search for such genes that were epigenetically controlled by the PI3K pathway in thyroid cancer cells, we found a wide range of genes with broad functions epigenetically targeted by the PI3K pathway. The most prominent among these genes was REC8, classically known as a meiotic-specific gene, which we found to be robustly down-regulated by the PI3K pathway through hypermethylation. REC8 hypermethylation was strongly associated with genetic alterations and activities of the PI3K pathway in thyroid cancer cell lines, thyroid cancer tumors, and some other human cancers; it was also associated with poor clinicopathological outcomes of thyroid cancer, including advanced disease stages and patient mortality. Demethylating the hypermethylated REC8 gene restored its expression in thyroid cancer cells in which the PI3K pathway was genetically over-activated and induced expression of REC8 protein inhibited the proliferation and colony formation of these cells. These findings are consistent with REC8 being a novel major bona fide tumor suppressor gene and a robust epigenetic target of the PI3K pathway. Aberrant inactivation of REC8 through hypermethylation by the PI3K pathway may represent an important mechanism mediating the oncogenic functions of the PI3K pathway.

Gene signatures of estrogen and progesterone receptor pathways predict the prognosis of colorectal cancer

The associations of estrogen receptor (ER) and progesterone receptor (PR) pathways with the prognosis of colorectal cancer (CRC) are still controversial. The aim of this study was to readdress these issues by introducing a gene signature-based approach to semiquantitate pathway activity. In this approach, the ER and PR pathway activities in CRC were computed based on the expression profiles of the signature genes of ER and PR pathways, respectively. The results showed that the ER pathway activity was progressively significantly decreased from normal colorectal mucosa, colorectal adenoma to CRC. ER pathway signaling was a favorable factor for the presence of microsatellite stability (MSS) in CRC in seven cohorts tested, while was an unfavorable factor for cancer recurrence in all four CRC cohorts tested (n = 1122; overall HR: 0.311, 95% CI: 0.199-0.488, P < 0.001). Subset stratification in stage II patients showed that ER pathway remained significantly inversely associated with recurrence. PR pathway was also suppressed in colorectal tumors and inversely associated with recurrence of CRC, but to a much lesser extent than ER pathway. Moreover, the inverse association of PR pathway with cancer recurrence was more likely observed in CRC with high ER pathway activity, suggesting the interactions between the two pathways. PR pathway was not associated with MSS in CRC, but it was more significant than ER pathway associated with advance cancer stages and cancer response to adjuvant chemotherapy. These results suggested the potential application of the gene signatures of ER and PR pathways, especially the former, as novel markers for prognosis and management of CRC.

Clinically relevant genes and regulatory pathways associated with NRASQ61 mutations in melanoma through an integrative genomics approach

Therapies such as BRAF inhibitors have become standard treatment for melanoma patients whose tumors harbor activating BRAFV600 mutations. However, analogous therapies for inhibiting NRAS mutant signaling have not yet been well established. In this study, we performed an integrative analysis of DNA methylation, gene expression, and microRNA expression data to identify potential regulatory pathways associated with the most common driver mutations in NRAS (Q61K/L/R) through comparison of NRASQ61-mutated melanomas with pan-negative melanomas. Surprisingly, we found dominant hypomethylation (98.03%) in NRASQ61-mutated melanomas. We identified 1,150 and 49 differentially expressed genes and microRNAs, respectively. Integrated functional analyses of alterations in all three data types revealed important signaling pathways associated with NRASQ61 mutations, such as the MAPK pathway, as well as other novel cellular processes, such as axon guidance. Further analysis of the relationship between DNA methylation and gene expression changes revealed 9 hypermethylated and down-regulated genes and 112 hypomethylated and up-regulated genes in NRASQ61 melanomas. Finally, we identified 52 downstream regulatory cascades of three hypomethylated and up-regulated genes (PDGFD, ZEB1, and THRB). Collectively, our observation of predominant gene hypomethylation in NRASQ61 melanomas and the identification of NRASQ61-linked pathways will be useful for the development of targeted therapies against melanomas harboring NRASQ61 mutations.

DNA methylation profiles in primary cutaneous melanomas are associated with clinically significant pathologic features

DNA methylation studies have elucidated a methylation signature distinguishing primary melanomas from benign nevi and provided new insights about genes that may be important in melanoma development. However, it is unclear whether methylation differences among primary melanomas are related to tumor pathologic features with known clinical significance. We utilized the Illumina GoldenGate Cancer Panel array to investigate the methylation profiles of 47 primary cutaneous melanomas. Arraywide methylation patterns revealed a positive association of methylation with Breslow thickness and mutated BRAF, a negative association with mitotic rate, and a weak association with ulceration. Hierarchical clustering on CpG sites exhibiting the most variable methylation (n = 235) divided the melanoma samples into three clusters, including a highly methylated cluster that was positively associated with Breslow thickness and an intermediately methylated cluster associated with Breslow thickness and mitotic rate. Our findings provide support for the existence of methylation-defined subsets in melanomas with increased methylation associated with Breslow thickness.

DNA methylation characteristics of primary melanomas with distinct biological behaviour

In melanoma, the presence of promoter related hypermethylation has previously been reported, however, no methylation-based distinction has been drawn among the diverse melanoma subtypes. Here, we investigated DNA methylation changes associated with melanoma progression and links between methylation patterns and other types of somatic alterations, including the most frequent mutations and DNA copy number changes. Our results revealed that the methylome, presenting in early stage samples and associated with the BRAF(V600E) mutation, gradually decreased in the medium and late stages of the disease. An inverse relationship among the other predefined groups and promoter methylation was also revealed except for histologic subtype, whereas the more aggressive, nodular subtype melanomas exhibited hypermethylation as well. The Breslow thickness, which is a continuous variable, allowed for the most precise insight into how promoter methylation decreases from stage to stage. Integrating our methylation results with a high-throughput copy number alteration dataset, local correlations were detected in the MYB and EYA4 genes. With regard to the effects of DNA hypermethylation on melanoma patients' survival, correcting for clinical cofounders, only the KIT gene was associated with a lower overall survival rate. In this study, we demonstrate the strong influence of promoter localized DNA methylation changes on melanoma initiation and show how hypermethylation decreases in melanomas associated with less favourable clinical outcomes. Furthermore, we establish the methylation pattern as part of an integrated apparatus of somatic DNA alterations.

Activities of multiple cancer-related pathways are associated with BRAF mutation and predict the resistance to BRAF/MEK inhibitors in melanoma cells

Drug resistance is a major obstacle in the targeted therapy of melanoma using BRAF/MEK inhibitors. This study was to identify BRAF V600E-associated oncogenic pathways that predict resistance of BRAF-mutated melanoma to BRAF/MEK inhibitors. We took in silico approaches to analyze the activities of 24 cancer-related pathways in melanoma cells and identify those whose activation was associated with BRAF V600E and used the support vector machine (SVM) algorithm to predict the resistance of BRAF-mutated melanoma cells to BRAF/MEK inhibitors. We then experimentally confirmed the in silico findings. In a microarray gene expression dataset of 63 melanoma cell lines, we found that activation of multiple oncogenic pathways preferentially occurred in BRAF-mutated melanoma cells. This finding was reproduced in 5 additional independent melanoma datasets. Further analysis of 46 melanoma cell lines that harbored BRAF mutation showed that 7 pathways, including TNFα, EGFR, IFNα, hypoxia, IFNγ, STAT3, and MYC, were significantly differently expressed in AZD6244-resistant compared with responsive melanoma cells. A SVM classifier built on this 7-pathway activation pattern correctly predicted the response of 10 BRAF-mutated melanoma cell lines to the MEK inhibitor AZD6244 in our experiments. We experimentally showed that TNFα, EGFR, IFNα, and IFNγ pathway activities were also upregulated in melanoma cell A375 compared with its sub-line DRO, while DRO was much more sensitive to AZD6244 than A375. In conclusion, we have identified specific oncogenic pathways preferentially activated in BRAF-mutated melanoma cells and a pathway pattern that predicts resistance of BRAF-mutated melanoma to BRAF/MEK inhibitors, providing novel clinical implications for melanoma therapy.