KMT2C methyltransferase domain regulated INK4A expression suppresses prostate cancer metastasis
Molecular Cancer volume 21, Article number: 89 (2022)
Frequent truncation mutations of the histone lysine N-methyltransferase KMT2C have been detected by whole exome sequencing studies in various cancers, including malignancies of the prostate. However, the biological consequences of these alterations in prostate cancer have not yet been elucidated.
To investigate the functional effects of these mutations, we deleted the C-terminal catalytic core motif of Kmt2c specifically in mouse prostate epithelium. We analysed the effect of Kmt2c SET domain deletion in a Pten-deficient PCa mouse model in vivo and of truncation mutations of KMT2C in a large number of prostate cancer patients.
We show here for the first time that impaired KMT2C methyltransferase activity drives proliferation and PIN formation and, when combined with loss of the tumour suppressor PTEN, triggers loss of senescence, metastatic dissemination and dramatically reduces life expectancy. In Kmt2c-mutated tumours we show enrichment of proliferative MYC gene signatures and loss of expression of the cell cycle repressor p16INK4A. In addition, we observe a striking reduction in disease-free survival of patients with KMT2C-mutated prostate cancer.
We identified truncating events of KMT2C as drivers of proliferation and PIN formation. Loss of PTEN and KMT2C in prostate cancer results in loss of senescence, metastatic dissemination and reduced life expectancy. Our data demonstrate the prognostic significance of KMT2C mutation status in prostate cancer patients. Inhibition of the MYC signalling axis may be a viable treatment option for patients with KMT2C truncations and therefore poor prognosis.
Prostate cancer (PCa) ranks as the second most frequently diagnosed malignancy in men worldwide and is expected to surpass even lung cancer incidence levels within the next decade [1, 2]. Diagnosis and therapy are challenged by enormous inter-tumour heterogeneity regarding clinical, morphological, and molecular features . While patients with localized or regional disease have an excellent prognosis, metastatic PCa remains largely incurable . Therefore, therapeutic strategies must be tailored to the individual risk of the patient to avoid overtreatment of low-risk tumours while ensuring rapid and decisive intervention in high-risk cases. To better stratify PCa and to advance the development of new therapies, a deeper understanding of the genetic and epigenetic events responsible for the progression and metastatic spread of PCa is urgently needed.
Multiple key pathways of prostate tumorigenesis have already been identified. Inactivating mutations of the phosphatase and tensin homolog (PTEN) tumour suppressor gene rank among the most common alterations observed in PCa . Loss of PTEN results in the aberrant activation of the phosphoinositide 3-kinase (PI3K) – AKT signalling pathway, which uncouples proliferation, survival and metabolism from external growth-stimulatory signals . However, loss of PTEN has also been shown to induce cellular senescence. This form of cell cycle arrest, which can be triggered upon oncogenic stress, is usually mediated via the p16INK4A-RB and the p14ARF-p53-p21CIP1 pathways and has previously been shown to act as a barrier to metastatic transformation in PCa . Besides inactivation of tumour suppressor genes, the amplification and overexpression of oncogenes is similarly known to play a crucial role in prostate tumorigenesis. The frequently altered androgen receptor (AR) signalling axis is the most well-studied pathway in the context of PCa. However, other key effectors, such as the activation of the proto-oncogene MYC, have also been found to be fundamental to PCa progression .
Besides genetic mutations, several epigenetic alterations, including DNA and histone modifications, have been identified in clinical PCa samples. Only recently has the pivotal importance of epigenetic reprogramming as a driver of carcinogenesis been widely recognized [9, 10]. Large cancer genome sequencing projects have revealed a substantial number of alterations in epigenetic modulators . The histone lysine N-methyltransferase KMT2C, an enzymatically active scaffold protein within the COMPASS (Complex Proteins Associated with Set1) multi-subunit complex, is the most frequently mutated gene within this group, predominantly presenting with frameshift and nonsense mutations [12, 13]. Previous studies have shown evidence of tumour suppressive roles for KMT2C and its close paralogue KMT2D, and have proposed their involvement in cellular growth, stemness and epithelial differentiation [13, 14]. However, despite the tremendous prevalence of mutations in these genes, the molecular mechanisms contributing to carcinogenic processes are still poorly understood . In PCa, alterations of epigenetic modifiers and chromatin-remodelling genes occur in about 20% of tumours with a strong overrepresentation of KMT2C mutations (7%), suggestive of a crucial role in carcinogenesis . A substantial number of mutations detected in KMT2C lead to loss of function of the encoded protein, whereby the catalytic domain is lost.
In this study, we investigated the effects of loss of the catalytic domain of KMT2C on the development and progression of PCa. We established a transgenic mouse model with prostate-specific deletion of either the catalytic core motif of KMT2C alone or in combination with loss of the tumour suppressor PTEN. In these models we observed that mutant KMT2C drives proliferation in vivo and triggers PCa metastasis when co-deleted with PTEN. Our data show that loss of the KMT2C catalytic core motif, mimicking the scenario in patients, results in an enrichment of the proliferative MYC gene signature and impairs p16INK4A-mediated cell cycle arrest in both our model system and human prostate cancers. Importantly, we show that mutated KMT2C significantly correlates with reduced disease-free survival (DFS) for PCa patients. Taken together, we identify the SET domain deletion of KMT2C as a novel driver of prostate carcinogenesis in murine models and suggest that the presence of mutated forms is a biomarker for poor outcome in PCa patients. Furthermore, our data are indicative of a possible therapeutic application through blockade of the MYC pathway.
KMT2C SET Domain Deletion Initiates Formation of Prostatic Intraepithelial Neoplasia In Vivo
Whole exome sequencing studies of various human cancers have identified frequent somatic mutations in the gene encoding the histone-methyltransferase KMT2C . Like other KMT2 proteins, KMT2C acts as a scaffold for the multi-subunit COMPASS complex where it regulates enhancer elements mainly through mono-methylation of lysine 3 on histone 4 (H3K4me1) via the enzymatically active SET domain located at the C-terminal end of the protein [13, 16] (Fig. 1a). To gain insight into the mutational spectrum of human PCa we analysed a cohort of 1013 patients with either localized or metastatic disease (MSKCC/DFCI cohort)  and found KMT2C to be the 7th most frequently mutated gene (Fig. 1b, left panel). In contrast to previously published data of the mutational pattern of KMT2C in different human cancers , KMT2C mutations in the analysed PCa dataset were distributed along the gene with no apparent mutational hotspot (Supplementary Fig. 1a) in keeping with the types of mutations commonly observed in tumour suppressor genes. While the functional significance of individual missense mutations is difficult to discern, truncation mutations (nonsense, frameshift, and splice mutations), which account for the majority of the changes we detected in KMT2C in the MSKCC/DFCI cohort, are primarily predicted to negatively impact the C-terminal SET domain, and are twice as common in metastatic samples in this dataset (Fig. 1b, right panel). Thus, we hypothesized that mutations leading to loss of the methyltransferase activity of KMT2C play a functional role in the initiation and/or progression of PCa. To investigate the impact of impaired KMT2C methyltransferase activity in vivo we established a mouse model with prostate-specific deletion of the Kmt2c SET domain-encoding exons (Kmt2cSETΔ/Δ). Specifically, we crossed mice carrying loxP sites flanking exons 57 and 58, which encode the catalytic core motif of KMT2C , to mice carrying the Cre recombinase transgene under the control of the androgen-regulated prostate-specific probasin promoter (PbCre4)  (Fig. 1c). Deletion of the Kmt2c SET domain-encoding exons was confirmed after the onset of puberty at both genetic and transcriptional levels by Polymerase Chain Reaction (PCR) and quantitative reverse transcription—PCR (RT-qPCR), respectively (Supplementary Fig. 1b-c). RT-qPCR analysis of total Kmt2c mRNA expression, using primers detecting a region outside of the deleted locus, revealed comparable levels of the mutated mRNA transcripts in Kmt2cSETΔ/Δ mouse prostates to wild type Kmt2c expression in control animals (Supplementary Fig. 1d). Interestingly, we found no compensatory up-regulation of the closely related paralogue Kmt2d (Supplementary Fig. 1e). Upon macroscopic evaluation, we discovered that the total prostate weight was significantly increased in mutant animals at 19 and 90 weeks postpartum (p.p.) compared to wild type controls (Fig. 1d, Supplementary Fig. 1f). Analysis of haematoxylin and eosin (H&E) stained sections of mouse prostates revealed focal areas of prostatic intraepithelial neoplasia (PIN) in Kmt2cSET∆/∆ mice sacrificed as early as 19 weeks p.p. increasing to full penetrance by 90 weeks p.p. (Fig. 1e-f). Immunohistochemistry (IHC) showed a significant increase in the percentage of Ki-67 positive proliferating cells in Kmt2cSET∆/∆ prostates (Fig. 1g-h, Supplementary Fig. 1g) and clusters of cells highly positive for AR expression, the primary driver of prostate cancer initiation and progression (Fig. 1i-j). These data suggest that loss of KMT2C catalytic activity drives proliferation and initiates transcriptional programs involved in prostate tumorigenesis evidenced by increased AR expression. However, even at 90 weeks of age mice showed only PIN, but no signs of PCa, suggesting that secondary events are required for full cancer progression. Thus, similar to many other prominent genetic alterations in PCa (e.g., affecting TP53 or ERG), KMT2C inactivation is insufficient to initiate malignant transformation.
Mutant KMT2C Drives Metastatic Transformation in a Pten-null Mouse Model of Prostate Cancer
To identify genes with a potential additive effect to KMT2C mutations in PCa progression, we analysed the MSKCC/DFCI PCa patient cohort and found a significant co-occurrence of alterations in PTEN amongst other genes (Supplementary Fig. 2a). PTEN is the most prominent tumour suppressor gene in PCa . Deleterious alterations of this gene are found in ~ 14% of primary and over 30% of metastatic PCa . Previous studies have shown that prostate-specific deletion of Pten is sufficient to induce tumorigenesis in mice and that disease progression closely mimics early human PCa [22,23,24]. Therefore, we back-crossed a mouse model in which conditional deletion of Pten is possible  to our Kmt2cSET∆/∆ animals to generate a double transgenic line with prostate-specific loss of both Kmt2c-SET and Pten (Pten∆/∆Kmt2cSET∆/∆) (Fig. 2a).
Prostate-specific deletion of these genes was verified at the genetic level by PCR (Supplementary Fig. 2b). Efficient abrogation of Pten mRNA expression levels were found for both Pten∆/∆ and Pten∆/∆Kmt2cSET∆/∆ prostates (Supplementary Fig. 2c). Deletion of the Kmt2c SET domain resulted in successful depletion of its full-length transcript. Expression levels of the mutant Kmt2c gene in Pten∆/∆Kmt2cSET∆/∆ double transgenic mice were comparable to full-length Kmt2c in Pten∆/∆ mice (Supplementary Fig. 2e). We observed no compensatory upregulation of Kmt2d upon Kmt2c SET domain deletion (Supplementary Fig. 2f). In line with the H3K4 mono-methyltransferase activity of KMT2C we found a global reduction of the enhancer mark H3K4me1 but not H3K27ac in Pten∆/∆Kmt2cSET∆/∆ prostate tissue (Supplementary Fig. 2g-h). Prostates of mice sacrificed at 19 weeks p.p. were significantly enlarged in Pten∆/∆Kmt2cSET∆/∆ animals compared to Pten∆/∆ controls (Fig. 2b-c). Interestingly, heterozygous deletion of Kmt2c-SET in combination with loss of Pten (Pten∆/∆Kmt2cSET∆/+) was sufficient to induce increased prostate weight indicative of tumour development (Fig. 2c).
We found a striking difference in the gross morphology between Pten∆/∆Kmt2cSET∆/∆ and Pten∆/∆ tumours. Deletion of Pten alone in the prostate epithelium led to formation of cystic tumours comprised of ~ 40% fluid, while Pten∆/∆Kmt2cSET∆/∆ animals developed solid tumours with cyst formation observed only in a minority of mice (Fig. 2d). Histopathological analysis of the primary tumour tissue revealed locally invasive areas in Pten∆/∆Kmt2cSET∆/∆ animals compared to age matched Pten∆/∆ mice. In PtenΔ/ΔKmt2cSETΔ/Δ prostates the basal cell layer was focally disrupted, and isolated tumour cells infiltrated into the surrounding stroma (Fig. 2e, Supplementary Fig. 2i). Double transgenic tumour cells infiltrating into blood vessels or nerve sheaths were observed at localised areas (Supplementary Fig. 2j). Furthermore, we found significantly more Ki-67 positive cells in Pten∆/∆Kmt2cSET∆/∆ tumours compared to PtenΔ/Δ controls (Fig. 2e-f). Strikingly, Pten∆/∆Kmt2cSET∆/∆ and even Pten∆/∆Kmt2cSET∆/+ mice showed a significant reduction in median life expectancy to 30 and 59 weeks, respectively, compared to > 80 weeks for Pten∆/∆ animals (p < 0.0001; Fig. 2g). Pten∆/∆ and Pten∆/∆Kmt2cSET∆/+ mice developed large tumours and had to be sacrificed due to tumour burden. In contrast, the general health of Pten∆/∆Kmt2cSET∆/∆ mice deteriorated rapidly after sudden disease onset despite significantly smaller tumour sizes. Necropsies of Pten∆/∆Kmt2cSET∆/∆ mice showed invasion of tumour cells into the urethra and urinary bladder, leading to severe obstructive uropathy with parenchymal reduction in the kidneys due to hydronephrosis (Fig. 3a, Supplementary Fig. 3a). This was accompanied by a sharp increase in serum creatinine levels, suggestive of renal failure (Supplementary Fig. 3b). These findings indicate a switch from organ confined PCa towards aggressive disease spreading beyond the prostate upon ablation of KMT2C activity in combination with PTEN loss.
To further investigate the invasive nature of Pten∆/∆Kmt2cSET∆/∆ prostate tumours we analysed draining lymph nodes and organs to identify metastatic tumour cells in mice showing clinical signs over 40 weeks of age. The prostatic origin of potential metastatic cells was determined by IHC staining of AR and keratin 8 (KRT8) (Fig. 3b). Lymph node metastases were found in more than 50% of analysed Pten∆/∆Kmt2cSET∆/∆ and 20% of Pten∆/∆Kmt2cSET∆/+ mice (Fig. 3b-c). Alongside lymphatic metastasis, prostatic tumour cells were detected in the blood vessels demonstrating metastatic dissemination via the circulatory system (Fig. 3b). Remarkably, 3 out of 10 Pten∆/∆Kmt2cSET∆/∆ mice had further progressed to develop lung metastases (Fig. 3b, d). In contrast, metastatic dissemination could not be found in lymph nodes or in distant organs in Pten∆/∆ mice up to 90 weeks p.p., which is in agreement with previously published data [25, 26] (Fig. 3c-d). Kaplan Meier analysis, whereby an event is defined as the presence of metastatic disease detected on necropsy, showed a median time of metastasis-free survival of 52.6 weeks in Pten∆/∆Kmt2cSET∆/∆ animals (Fig. 3e). This suggests that a significant number of animals likely die due to renal failure from bladder obstruction, resulting from invasive growth of PCa into the local surrounding tissues, before distant metastases of a detectable size can form. Taken together, these data reveal that loss of KMT2C histone methyltransferase activity in Pten-null prostate tumours not only drives proliferation, but also confers invasive properties allowing tumour cells to invade surrounding tissues and spread to distant organs.
Transcriptional Profiling of KMT2C Mutated Prostate Epithelial Cells Reveals Enrichment of the Proliferative MYC Gene Signature
We performed RNA sequencing (RNA-Seq) of prostate tissue from wild type, Kmt2cSETΔ/Δ, PtenΔ/Δ and PtenΔ/ΔKmt2cSETΔ/Δ mice at 19 weeks p.p. following the onset of tumour formation to gain insights into alterations in gene expression that are associated with the observed phenotypes. Cells expressing the epithelial cell adhesion molecule (EpCAM) were isolated from prostate tissue using a magnetic bead-based cell sorting approach. Successful enrichment of EpCAMpos cells was confirmed using fluorescence activated cell sorting (FACS) (Fig. 4a) before RNA isolation and RNA-Seq was conducted. Principal component analysis revealed a clear separation of tumour samples (PtenΔ/Δ and PtenΔ/ΔKmt2cSETΔ/Δ) according to their genotypes. In contrast, Kmt2cSETΔ/Δ samples clustered closely with wild type samples, indicating only minor changes in overall gene expression between these groups (Supplementary Fig. 4a). These results substantiate our previous data, highlighting the importance of KMT2C in tumour progression rather than initiation. Whereas differential gene expression analysis of the transcriptomes of PtenΔ/ΔKmt2cSETΔ/Δ with Pten-null prostate epithelial cells showed significant up-regulation of 252 and down-regulation of 943 genes (Fig. 4b) the number of significantly deregulated genes between the wild type and Kmt2cSETΔ/Δ cells was considerably smaller with only 5 up- and 80 downregulated transcripts in Kmt2cSETΔ/Δ mice (Supplementary Fig. 4b). To gain further insight into the pathways affected by these gene expression changes, we performed fast pre-ranked gene set enrichment analysis (fGSEA) using hallmark gene sets derived from MSigDB [27, 28]. In accordance with our previous results, we found a significant enrichment for genes driving proliferation upon loss of the KMT2C methyltransferase domain (Fig. 4c, Supplementary Fig. 4c). Of note, we found the proto-oncogene Myc, which is a known master regulator of proliferation, and its downstream targets Ccnd1 and Ccnd2 to be upregulated in PtenΔ/ΔKmt2cSETΔ/Δ samples (Fig. 4d, Supplementary Fig. 4d). In line with this observation, MYC target genes were overrepresented in Kmt2c SET domain deleted samples (Fig. 4e, Supplementary Fig. 4e). Overexpression of MYC in the progression of PCa has previously been observed [8, 29]. Additionally, MYC expression has been shown to correlate with increased disease severity and to positively regulate the well described PCa driver AR [30, 31]. In line with this dependency, we observed an upregulation of androgen response genes in PtenΔ/ΔKmt2cSETΔ/Δ samples (Supplementary Fig. 4f). Besides its role in cellular proliferation, MYC has been recognized as regulating epithelial-to-mesenchymal (EMT) transition and in the general acquisition of invasive properties of cancer cells . Accordingly, an enrichment of genes involved in EMT was observed upon Kmt2c SET domain deletion (Fig. 4f, Supplementary Fig. 4 g). To further confirm the invasive nature of PtenΔ/ΔKmt2cSETΔ/Δ prostate epithelial cells we analysed two additional sets of genes described in general to be upregulated in metastatic tumours (RAMASWAMY Metastasis_UP)  or lymph node metastasis in PCa (PCa LN Metastasis_UP)  and found both signatures to be enriched in the double transgenic samples (PtenΔ/ΔKmt2cSETΔ/Δ) compared to controls (Fig. 4g-h). Taken together, changes in the transcriptome observed upon disruption of KMT2C methyltransferase activity verify the activation of proliferative signalling pathways and indicate alterations in processes involved in disease aggressiveness.
Mutant KMT2C Impairs Oncogene-Induced Expression of the Cell Cycle Repressor p16INK4A
While sustained growth-stimulatory signalling is a fundamental trait of cancer cells, the replicative stress induced by excessive proliferation may elicit a counteractive growth arrest. This tumour suppressive process, known as oncogene-induced senescence (OIS), can present a barrier to malignant transformation of precursor lesions and must be overcome for tumours to progress to lethal, metastatic disease [35, 36]. Elevated signalling by strong oncogenic drivers, as well as loss of potent tumour suppressors, such as PTEN, have been described to trigger OIS, highlighting the potentially crucial role of this phenomenon in our model system [35, 37]. The unrestrained proliferation and metastatic dissemination of PtenΔ/ΔKmt2cSETΔ/Δ prostate tumour cells described so far indicate an escape from PTEN-loss induced cellular senescence (PICS) in KMT2C mutant tumours. Another key attribute of senescent cells is the extensive change in expression of inflammatory cytokines and other signalling molecules known as the senescence-associated secretory phenotype (SASP) . Interestingly, analysis of the RNA-Seq data (Fig. 4) revealed a significant deregulation of immune response pathways between the transcriptomes of PtenΔ/Δ and PtenΔ/ΔKmt2cSETΔ/Δ prostate epithelial cells (Fig. 4c). This might reflect changes to the SASP as there is considerable overlap of genes involved in those processes. To further investigate potential loss of the SASP gene signature in PtenΔ/ΔKmt2cSETΔ/Δ prostate cells we performed fGSEA using a set of genes previously described to be induced upon PICS (“core SASP of PICS”)  and found them to be strongly downregulated in the PtenΔ/ΔKmt2cSETΔ/Δ double transgenic group compared to PtenΔ/Δ mouse prostate cells (Fig. 5a). Next, we performed IHC staining of mouse prostate tissue taken from mice at 19 weeks p.p. for the senescence markers β-galactosidase (GLB1) and p16INK4A. PtenΔ/Δ prostate cells showed an accumulation of GLB1 and pronounced induction of p16INK4A which corroborates previously published data describing senescence in the PtenΔ/Δ PCa mouse model . In contrast, expression of both markers was lower in PtenΔ/ΔKmt2cSETΔ/Δ prostate samples (Fig. 5b-d). As senescence is known to be mediated by the two isoforms of Cdkn2a, p16INK4A and p19ARF, and the Cdkn1a gene product p21CIP1  we analysed expression of these genes via RT-qPCR. In line with the induction of senescence upon loss of Pten, expression of Cdkn2a transcripts and, to a lesser extent, Cdkn1a was induced in PtenΔ/Δ samples. Combined loss of Pten with the Kmt2c SET domain led to the loss of p16INK4A expression and a similar but less pronounced effect for p19ARF (Fig. 5e). Western blot analysis for p16INK4A and the cyclin-dependent kinase CDK4, which is known to be inhibited by p16INK4A, confirmed the loss of p16INK4A detected in PtenΔ/ΔKmt2cSETΔ/Δ protein lysates and an upregulation of CDK4. In addition, we analysed p53 expression, the key mediator of p19ARF-p53-p21CIP1-mediated senescence but did not detect a deregulation of this axis at the protein level nor at the transcript level as indicated by no significant change in p21 levels (Fig. 5e-g). We therefore conclude that mutant KMT2C impairs the induction of OIS in Pten-null mice by circumventing p16INK4A-mediated cell cycle arrest.
KMT2C Truncation Mutations are Associated with Increased Proliferation and Reduced Disease-Free Survival for Prostate Cancer Patients
We have identified KMT2C SET domain deletion as a driving event of metastatic PCa in our murine model system and have observed an enrichment of the MYC gene signature accompanied by concurrent loss of p16INK4A as mechanistic downstream effects. Therefore, we hypothesized that KMT2C truncation mutations may be associated with aggressive disease and poor prognosis in PCa patients. Therefore, we stratified sample data derived from the Cancer Genome Atlas Prostate Adenocarcinoma (TCGA-PRAD)  patient cohort into KMT2C wildtype and KMT2C truncated tumours for further pathway analysis. Analogous to our murine data, fGSEA of hallmark gene sets revealed strong upregulation of proliferative pathways, with MYC target genes ranking as the most highly enriched gene sets (Fig. 6a-b). Furthermore, mRNA expression of MYC itself was upregulated in patient samples carrying KMT2C truncation mutations (Fig. 6c). In line with the capacity of MYC to induce transcription and enhance the stability of the AR, we detected enrichment of androgen response genes in the KMT2C truncated patient group (Supplementary Fig. 5a). As our murine data show that activation of proliferative signalling pathways occurs concurrently with evasion of p16INK4A-mediated growth arrest, we performed fGSEA on the “core SASP of PICS” gene set and found a significant reduction of this transcriptional signature in KMT2C truncated tumours (Supplementary Fig. 5b), potentially indicating loss of senescence features. We further observed a reduction in expression of the gene encoding p16INK4A, CDKN2A, in these samples (Fig. 6d). As p16INK4A inhibits cell cycle progression at the G1/S transition, we performed fGSEA of two sets of genes upregulated during this process (REACTOME G1/S, FISCHER G1/S) and found these to be highly enriched in the KMT2C truncated patient group, further highlighting the likely depletion of cell cycle repressors in this cohort (Fig. 6e). Collectively, these data indicate hyperactivation of proliferation in PCa carrying truncated forms of KMT2C, which is remarkably similar to the effect observed in our model system. Thus, we hypothesize that KMT2C mutations, in agreement with our murine data, also influence prognosis in PCa. Therefore, we analysed data derived from the International Cancer Genome Consortium (ICGC)  and found that KMT2C truncation mutations (nonsense and frameshift mutations) correlate with reduced disease-free survival (DFS) in PCa (Fig. 6f). Taken together our data show that KMT2C truncating events drive cancer progression by upregulating MYC target genes and circumventing p16INK4A-mediated growth arrest. These effects may ultimately facilitate the formation of lethal, metastatic disease in PCa patients. Therefore, KMT2C truncation mutations represent a biomarker for aggressive disease and indicate that inhibitors of the MYC signalling axis might be beneficial to these patients.
Although KMT2C is the most frequently mutated epigenetic regulator in PCa , and loss of histone methylation correlates with a poorer clinical outcome , functional studies investigating the impact of these mutations on prostate carcinogenesis are scarce. Here, we describe for the first time that loss of the KMT2C methyltransferase domain accelerates tumour growth and promotes a switch from indolent to lethal, metastatic disease in vivo when combined with other PCa associated mutational events, specifically loss of PTEN. Our study also reveals that this exacerbation of tumour progression is accompanied by an enrichment of the proliferative MYC gene signature and an escape of p16INK4A-mediated oncogene-induced senescence.
Previous studies have shown correlations between expression levels of KMT2C and cancer progression, while others have instead highlighted the importance of genetic alterations affecting protein function [17, 42,43,44]. To date, the contribution of KMT2C to carcinogenesis has most extensively been studied in breast cancer (BCa), which is biologically similar to PCa due to its dependence on steroid hormone signalling. Reduced expression of KMT2C or KMT2D results in decreased oncogenic estrogen receptor (ER) signalling, possibly due to impaired interaction with the ER pioneer factor Forkhead box protein A1 (FOXA1) [44, 45]. In contrast, low KMT2C expression correlates with bad prognosis in BCa  and both truncation mutations as well as mutations in the plant homeodomain (PHD) of KMT2C have been shown to be tumour-promoting events [17, 44]. Thus, the contribution of KMT2C to carcinogenesis appears to be influenced by cellular context, alteration type, as well as hormone-dependence even within a single malignancy. Our analysis of the mutational spectrum in PCa patients showed a high prevalence of nonsense and frameshift mutations resulting in truncated forms of KMT2C, especially in metastatic samples. We therefore focused our work on the genetic loss of the C-terminal methyltransferase domain rather than the analysis of expressional changes.
To study the complex biological consequences of KMT2C alterations in vivo, we established a genetically engineered mouse model (GEMM) with a prostate epithelium specific deletion of the catalytic core motif of Kmt2c. Kmt2c-mutated animals developed PIN lesions but never progressed to develop carcinomas. This observation resembles other GEMMs modelling key alterations observed in human PCa, such as deletion of TP53 or overexpression of ERG or AR in the prostate epithelium [25, 46, 47]. Regardless, loss of the Kmt2c SET domain resulted in the activation of the oncogenic AR signalling pathway and a significant increase in cellular proliferation, albeit at low levels. These data are in contrast to the ability of wildtype KMT2C/D in facilitating the closely related ER signalling axis in BCa [44, 45]. A potential difference in interaction of ER and AR with their co-regulators, including FOXA1, might result in diverging effects of KMT2C on these transcriptional programs. However, the varied effects of KMT2C observed across different BCa model systems hints at a multitude of factors capable of influencing the biological outcome of KMT2C alterations in carcinogenesis.
Previous studies have identified loss of PTEN or overexpression of MYC to be sufficient to induce carcinogenesis in vivo. However, neither alteration leads to metastatic disease suggesting a requirement for other contributory genetic events [23, 48]. Formation of distant metastasis has previously been achieved in the Pten-null model in combination with select additional genetic alterations [26, 39, 49,50,51]. For example, it was recently shown that the combined overexpression of MYC and loss of Pten in prostate luminal epithelial cells can induce lethal metastatic PCa . Our analysis of the mutational spectrum present in PCa patients revealed a significant co-occurrence of PTEN and KMT2C alterations, pointing towards a synergistic effect. In accordance, we demonstrate in our study that Kmt2c SET domain deletion represents an as-yet unknown aberration capable of initiating the important switch from localized to metastatic PCa in a Pten-null GEMM.
A frequent observation in relation to metastatic transformation is the circumvention of the senescent phenotype initiated upon Pten-deficiency , or upregulation of the proto-oncogene MYC [50,51,52]. There is ample evidence for the fundamental role of senescence as an abortive stage of cancer development. However, its impact on tumorigenesis appears to be diverse and highly dependent on the context in which it arises. A multitude of inducers of senescence have been described, yet the molecular effectors activated downstream mainly converge on two major tumour suppressive pathways: the p16INK4A-RB and the p14ARF-p53-p21CIP1 signalling axes . Interestingly, genetic deletion of p19Arf, the murine homolog to p14ARF, in Pten-null mice does not impact p53 expression in vivo and is not sufficient to abrogate senescence or promote tumorigenesis . Furthermore, while senescent features are lost in Pten / p53 double-knockout prostate epithelial cells in a GEMM, the formation of distant metastases was not observed . These findings indicate that circumvention of OIS may be insufficient to drive metastatic transformation and needs to occur in parallel with the induction of an additional driver of cellular dissemination. In our model system, combined deletion of Pten and the Kmt2c catalytic core motif led to loss of p16INK4A in double transgenic mice but failed to convincingly show any change in the p19ARF-p53-p21CIP1 axis. Still, we observed a loss of senescence features and, more importantly, a drastic reduction in life expectancy due to obstructive tumour cell infiltration into the urethra and the development of lymph node and lung metastases. Our data additionally suggest a possible involvement of the MYC signalling axis in the switch to metastatic PCa in our model system based on increased expression of MYC target genes and an upregulation of Myc mRNA levels. However, as other drivers of proliferation might activate a similar set of downstream targets and MYC protein levels and activity are tightly controlled not only on transcriptional level but also through RNA stability, protein turnover and posttranslational modifications  the transcriptional changes observed in our model system may not exclusively point towards MYC but drivers of proliferation in general.
Based on the functional role of KMT2C as an H3K4 mono-methyltransferase, loss of p16INK4A might directly be regulated by impaired placement of activating H3K4me1 marks at enhancer elements associated with Cdkn2a. On the other hand, an upregulation of Myc is unlikely to be a direct effect of H3K4me1 loss and possibly mediated by the downregulation of a transcriptional inhibitor of Myc instead. While our data convincingly show a global deregulation of H3K4me1 upon impaired KMT2C methyltransferase activity deciphering site specific changes in the enhancer landscape in our model system and how they correlate with transcriptional changes remains to be investigated.
The relevance of our findings to human patients is supported by the fact that we observe both enrichment of the MYC gene signature and downregulation of CDKN2A, which encodes p16INK4A, in PCa patients with truncation mutations of KMT2C. Remarkably, we found a significant correlation between DFS and truncated KMT2C in patients, demonstrating the prognostic significance of KMT2C mutation status in PCa. Based on our findings, we speculate that inhibition of MYC transcriptional activity may be a viable treatment option for patients with KMT2C truncation mutations. Furthermore, elucidating the complex transcriptional networks that are altered in response to KMT2C mutations uncovered in this study may facilitate the identification of novel pro-metastatic pathways and promote the development of new clinical therapeutics that can counteract PCa metastases.
Taken together, this work reveals new insights into the previously poorly understood transition from local to lethal metastatic PCa. Furthermore, data presented in this study provide a rationale for the inclusion of KMT2C mutation status in standard diagnostics of patients with suspected PCa. Based on our findings, inhibition of MYC-associated transcriptional activity could represent a strategy for treating PCa patients with deleterious KMT2C mutations and thus a poor prognosis. As many other human cancers also show a high KMT2C mutation prevalence our work may have translational relevance to an array of malignancies similarly affected.
Analysis of mutation data from publicly available datasets
The DFS Rate of PCa patients dependent on KMT2C mutational status was derived from the ICGC Data Portal. For “KMT2C wildtype” all patients with KMT2C mutations of any kind, including those of unknown significance, were removed. Mutations of the consequence type “Frameshift” and “Stop Gained” were grouped as truncation mutations. Consequence type “Missense” was selected for missense mutations. All available PCa projects (PRAD-US, PRAD-CA, PRAD-UK, EOPC-DE, PRAD-CN, PRAD-FR) were included in this analysis.
Pb-Cre4 (Tg(Pbsn-cre)4Prb, RRID: IMSR_JAX:026,662), Ptenflox/flox (Ptentm2Mak, RRID: IMSR_RBRC02300) and Kmt2cflox/flox (Kmt2ctm1.2Jwle, RRID: MGI:4,948,101) mice have previously been described [18, 19, 24]. All animals were maintained on a C57Bl/6-Sv/129 mixed background. Male mice were allocated to experimental groups based on their genotype and sacrificed at 19 or 90 weeks of age as indicated in the individual figure legends. All animal studies were reviewed and approved by the Federal Ministry Republic of Austria for Education, Science and Research and conducted according to regulatory standards (BMBWF GZ 66.009/0135-WF/V/3b/2016).
Sections were reviewed by two independent board-certified pathologists with specific expertise in laboratory animals (L.K. and S.H.). All analyses were performed blinded to genotype.
Immunohistochemistry and histological analysis
Immunohistochemistry (IHC) and haematoxylin/eosin stains (H&E) were performed with formalin-fixed paraffin-embedded (FFPE) tissue using standard protocols. The following antibodies were used for IHC: Ki67 (Cell Signaling, #12,202, 1:400), AR (Abcam, ab133273, 1:100), p16 (Abcam, ab211542, 1:500), GLB1 (Novus Biologicals, NBP2-45,731, 1:120), KRT8 (Abcam, ab53280, 1:100) and p63 (Abcam, ab735, 1:80).
All images were taken with a Zeiss AxioImager M2. Quantification was performed using QuPath 0.1.2 . If automated quantification was not possible, semi-quantitative analysis was performed by a board-certified pathologist with specific expertise in mouse models of prostate cancer (L.K.). Analyses were performed blinded to genotype.
Western blot analysis
Protein extraction from frozen prostate samples for Western blot analysis has previously been described . Western blots were prepared using 15–20 µg of protein lysate, blocked with 5% BSA in 1 × TBS / 0.1% Tween-20 for 1 h and incubated at 4 °C overnight with primary antibodies against AR (1:1000, ab133273, Abcam), p16INK4A (1:2000, ab211542, Abcam), CDK4 (1:200, sc-260, Santa Cruz), p53 (1:1000, CST#2524, Cell Signaling), β-Actin (1:1000, CST#4967, Cell Signaling) and β-Tubulin (1:2000, CST#2146, Cell Signaling). Western blots were quantified using ImageJ2.
RNA Isolation and RT-qPCR
RNA from mouse tissue was extracted using TRI Reagent (Merck) and purified using the ReliaPrep RNA Tissue Miniprep kit (Promega). DNase digestion was performed on a column. For RT-qPCR, 1 µg of total RNA was reverse transcribed to cDNA using the iScript cDNA Synthesis Kit (Bio-Rad). RT-qPCR was performed in triplicates with the Luna Universal qPCR Master Mix (NEB). All procedures were performed according to the manufacturers’ instructions. Real-time monitoring of PCR amplification was performed using the CFX96 Real-Time System C1000 Thermal Cycler (Bio-Rad). mRNA levels were calculated using the Pfaffl analysis method  and normalized to the geometric mean of β-actin and cyclophilin A. RT-qPCR primer sequences are listed in the Supplementary Materials and Methods.
RNA-Seq of mouse prostate epithelial cells
Mice were sacrificed at 19 weeks p.p. and prostates were dissected. Prostate tissue was digested as previously described . Cell suspensions were strained through a 40 µm cell strainer and washed twice using MACS buffer (1 × PBS + 2 mM EDTA + 2% FCS). The cell solution was subject to centrifugation for 5 min at 150 g, and cells were resuspended in 800 µL MACS buffer. Cells were passed 5 times through a 27G needle and counted. 1–3 * 107 cells were collected by centrifugation at 300 g for 5 min and resuspended in 200 µL MACS buffer. 100 ng anti-CD326-biotin (13–5791-82, eBioscience) were added to the cell solution and incubated at room temperature for 10 min. 1 mL MACS buffer was added to the solution and cells were collected by centrifugation at 300 g for 5 min. The cell pellet was resuspended in 200 µL MACS buffer and transferred to a 12 × 75 mm FACS tube. 60 µL streptavidin positive selection beads (T9424, Sigma) were added to the tubes and incubated for 10 min at room temperature. 2.5 mL MACS buffer were added to the cell-bead suspension and the tube was transferred to a magnetic stand. Beads were allowed to adhere to the magnet for 5 min at room temperature before MACS buffer was discarded. Beads were washed twice more in an identical manner using 2.5 mL MACS buffer for each wash and resuspending the cell-bead solution with a P1000 pipette between washes. The bound fraction containing CD326-positive cells was then collected in 5 mL MACS buffer and subject to centrifugation for 5 min at 150 g at 4 °C. Cells were resuspended in 200 µL MACS buffer, counted, and transferred to a fresh (RNase-free) Eppendorf tube. 5 µL of sorted prostate epithelial cells were removed for FACS analysis. The remaining cells were pelleted at 300 g for 5 min at 4 °C, supernatant was removed, and the cell pellet was snap frozen in liquid nitrogen and stored at -80 °C until further use.
RNA isolation was performed using the ReliaPrep RNA Tissue Miniprep kit (Promega). RNA quality was assessed using the 4200 TapeStation System (Agilent). For library preparation, the Ultra II Directional RNA Library Prep Kit (E7760, NEB) was used in combination with a poly(A) mRNA magnetic isolation module (E7490) and multiplex oligos for Illumina (E7600). Library preparation was performed according to the manufacturer’s instructions (optimized for an insert size of 200 nt, input: 250 ng, adaptor dilution: fivefold, 11 amplification cycles). The library quality was analysed using the 4200 TapeStation System.
RNA-Seq data analysis
Raw FASTQ files from murine RNA-Seq were subject to a quality check with FastQC (v0.11.5)  and MultiQC (v1.4) . Adapters and low-quality read ends were trimmed using Trimmomatic (v0.36)  and reads shorter than 35 nt were discarded. The pre-processed reads were mapped to the reference mouse genome (Ensembl GRCm38) utilizing gene annotation (Ensembl v91) with STAR [64, 65]. The quality of mapping was evaluated with RSeQC (v2.6.4)  and Picard (v2.10.6)  and rRNA content was checked with FastQ Screen (v0.13.0) . Gene quantification was performed on uniquely mapped reads only, with featureCounts (v1.5.2) . Differential expression analysis was carried out in R (version 3.5.1)  with DESeq2 package (v1.20.0)  and limma package (v3.38.2) . Genes with a FDR-adjusted P value < 0.05 and fold change ≥ 2 or fold change ≤ 0.5 were considered significantly differentially expressed. Heatmaps of differentially expressed genes (DEG) were generated using unsupervised hierarchical clustering with the pheatmap package (v1.0.12) . fGSEA of gene sets was done with the fgsea package (v1.14.0) . Gene sets were derived from MSigDB (version 6.1.1) (hallmark gene sets, Ramaswamy_metastasis_up, Reactome_G1_S_transition, Fischer G1_S_cell_cycle) or previously published works (core SAPS of PICS , PCa LN metastasis ). Gene sets based on previously published works are detailed in the Supplementary Materials and Methods. Human gene symbols present in custom gene sets were converted to orthologous mouse genes using the biomaRt package (v2.44.4) .
TCGA data analysis
Harmonized TCGA PRAD  RNA-seq data were acquired as HTSeq‐Counts via R package TCGAbiolinks (v2.16.4) . Differentially expressed genes between patients with truncated and wild-type forms of KMT2C were identified with DESeq2 package (v1.20.0) . fGSEA of hallmark gene sets (MSigDB version 6.1.1) and custom gene sets was performed with the fgsea package (v1.14.0) . Detailed steps and parameters of the analysis are described in the Supplementary Materials and Methods.
Data in figures are plotted as individual replicates with their mean and standard deviation for analyses with n ≤ 20 per group or as boxplots depicting the 25th to 75th percentile (box) and min to max (whiskers) for n > 20 per group. DFS rate for data derived from the ICGC database was calculated using the “cohort comparison” tool provided by the platform. All other statistical analyses were performed using GraphPad Prism 8. The significance level of differences between groups was determined by two-tailed unpaired Student’s t-tests for 2 groups or ordinary one-way ANOVA with Tukey’s multiple comparisons tests for 3 or more groups. For Kaplan–Meier analysis the log rank (Mantel-Cox) test was performed.
Availability of data and materials
The RNA-Seq dataset supporting the conclusions of this article is available in the GEO repository, GSE186413 and is publicly available as of date of publication. The following publicly available datasets were used: The Cancer Genome Atlas—Prostate Adenocarcinoma (https://portal.gdc.cancer.gov/projects/TCGA-PRAD). The long tail of oncogenic drivers in prostate cancer , ICGC Prostate Cancer Project (PRAD-US, PRAD-CA, PRAD-UK, EOPC-DE, PRAD-CN, PRAD-FR) .
Complex of proteins associated with Set1
Differentially expressed genes
Epithelial cell adhesion molecule
Fluorescence activated cell sorting
Fast pre-ranked gene set enrichment analysis
Genetically engineered mouse model
Haematoxylin and eosin
Mono-methylation of lysine 4 on histone 3
International Cancer Genome Consortium
Prostatic intraepithelial neoplasia
PTEN-loss induced cellular senescence
Probasin Cre promoter
Polymerase Chain Reaction
Quantitative reverse transcription-PCR
Phosphate and Tensin Homologue
Senescence-associated secretory phenotype
Cancer Genome Atlas Prostate Adenocarcinoma
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F: Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.
Cancer Tomorrow [https://gco.iarc.fr/]
Haffner MC, Zwart W, Roudier MP, True LD, Nelson WG, Epstein JI, De Marzo AM, Nelson PS, Yegnasubramanian S. Genomic and phenotypic heterogeneity in prostate cancer. Nat Rev Urol. 2021;18:79–92.
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer Statistics, 2021. CA Cancer J Clin. 2021;71:7–33.
Jamaspishvili T, Berman DM, Ross AE, Scher HI, De Marzo AM, Squire JA, Lotan TL. Clinical implications of PTEN loss in prostate cancer. Nat Rev Urol. 2018;15:222–34.
Hoxhaj G, Manning BD. The PI3K-AKT network at the interface of oncogenic signalling and cancer metabolism. Nat Rev Cancer. 2020;20:74–88.
Parisotto M, Grelet E, El Bizri R, Dai Y, Terzic J, Eckert D, Gargowitsch L, Bornert JM, Metzger D. PTEN deletion in luminal cells of mature prostate induces replication stress and senescence in vivo. J Exp Med. 2018;215:1749–63.
Koh CM, Bieberich CJ, Dang CV, Nelson WG, Yegnasubramanian S, De Marzo AM. MYC and prostate cancer. Genes Cancer. 2010;1:617–28.
Feinberg AP, Koldobskiy MA, Gondor A. Epigenetic modulators, modifiers and mediators in cancer aetiology and progression. Nat Rev Genet. 2016;17:284–99.
Perdigoto CN. Epigenetic cancer evolution, one cell at a time. Nat Rev Genet. 2019;20:434–5.
Chatterjee A, Rodger EJ, Eccles MR. Epigenetic drivers of tumourigenesis and cancer metastasis. Semin Cancer Biol. 2018;51:149–59.
Kandoth C, McLellan MD, Vandin F, Ye K, Niu B, Lu C, Xie M, Zhang Q, McMichael JF, Wyczalkowski MA, et al. Mutational landscape and significance across 12 major cancer types. Nature. 2013;502:333–9.
Rao RC, Dou Y. Hijacked in cancer: the KMT2 (MLL) family of methyltransferases. Nat Rev Cancer. 2015;15:334–46.
Fagan RJ, Dingwall AK. COMPASS Ascending: Emerging clues regarding the roles of MLL3/KMT2C and MLL2/KMT2D proteins in cancer. Cancer Lett. 2019;458:56–65.
Armenia J, Wankowicz SAM, Liu D, Gao J, Kundra R, Reznik E, Chatila WK, Chakravarty D, Han GC, Coleman I, et al. The long tail of oncogenic drivers in prostate cancer. Nat Genet. 2018;50:645–51.
Lavery WJ, Barski A, Wiley S, Schorry EK, Lindsley AW. KMT2C/D COMPASS complex-associated diseases [KCDCOM-ADs]: an emerging class of congenital regulopathies. Clin Epigenetics. 2020;12:10.
Wang L, Zhao Z, Ozark PA, Fantini D, Marshall SA, Rendleman EJ, Cozzolino KA, Louis N, He X, Morgan MA, et al. Resetting the epigenetic balance of Polycomb and COMPASS function at enhancers for cancer therapy. Nat Med. 2018;24:758–69.
Lee S, Lee DK, Dou Y, Lee J, Lee B, Kwak E, Kong YY, Lee SK, Roeder RG, Lee JW. Coactivator as a target gene specificity determinant for histone H3 lysine 4 methyltransferases. Proc Natl Acad Sci U S A. 2006;103:15392–7.
Wu X, Wu J, Huang J, Powell WC, Zhang J, Matusik RJ, Sangiorgi FO, Maxson RE, Sucov HM, Roy-Burman P. Generation of a prostate epithelial cell-specific Cre transgenic mouse model for tissue-specific gene ablation. Mech Dev. 2001;101:61–9.
Cancer Genome Atlas Research Network. The Molecular taxonomy of primary prostate cancer. Cell. 2015;163:1011–25.
Nguyen B, Mota JM, Nandakumar S, Stopsack KH, Weg E, Rathkopf D, Morris MJ, Scher HI, Kantoff PW, Gopalan A, et al. Pan-cancer analysis of CDK12 alterations identifies a subset of prostate cancers with distinct genomic and clinical characteristics. Eur Urol. 2020;78:671–9.
Wang S, Gao J, Lei Q, Rozengurt N, Pritchard C, Jiao J, Thomas GV, Li G, Roy-Burman P, Nelson PS, et al. Prostate-specific deletion of the murine Pten tumor suppressor gene leads to metastatic prostate cancer. Cancer Cell. 2003;4:209–21.
Trotman LC, Niki M, Dotan ZA, Koutcher JA, Di Cristofano A, Xiao A, Khoo AS, Roy-Burman P, Greenberg NM, Van Dyke T, et al. Pten dose dictates cancer progression in the prostate. PLoS Biol. 2003;1:E59.
Suzuki A, Yamaguchi MT, Ohteki T, Sasaki T, Kaisho T, Kimura Y, Yoshida R, Wakeham A, Higuchi T, Fukumoto M, et al. T cell-specific loss of Pten leads to defects in central and peripheral tolerance. Immunity. 2001;14:523–34.
Chen Z, Trotman LC, Shaffer D, Lin HK, Dotan ZA, Niki M, Koutcher JA, Scher HI, Ludwig T, Gerald W, et al. Crucial role of p53-dependent cellular senescence in suppression of Pten-deficient tumorigenesis. Nature. 2005;436:725–30.
Aytes A, Mitrofanova A, Kinkade CW, Lefebvre C, Lei M, Phelan V, LeKaye HC, Koutcher JA, Cardiff RD, Califano A, et al. ETV4 promotes metastasis in response to activation of PI3-kinase and Ras signaling in a mouse model of advanced prostate cancer. Proc Natl Acad Sci U S A. 2013;110:E3506-3515.
Liberzon A, Subramanian A, Pinchback R, Thorvaldsdottir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27:1739–40.
Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov JP, Tamayo P. The molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1:417–25.
Gurel B, Iwata T, Koh CM, Jenkins RB, Lan F, Van Dang C, Hicks JL, Morgan J, Cornish TC, Sutcliffe S, et al. Nuclear MYC protein overexpression is an early alteration in human prostate carcinogenesis. Mod Pathol. 2008;21:1156–67.
Hawksworth D, Ravindranath L, Chen Y, Furusato B, Sesterhenn IA, McLeod DG, Srivastava S, Petrovics G. Overexpression of C-MYC oncogene in prostate cancer predicts biochemical recurrence. Prostate Cancer Prostatic Dis. 2010;13:311–5.
Bai S, Cao S, Jin L, Kobelski M, Schouest B, Wang X, Ungerleider N, Baddoo M, Zhang W, Corey E, et al. A positive role of c-Myc in regulating androgen receptor and its splice variants in prostate cancer. Oncogene. 2019;38:4977–89.
Meskyte EM, Keskas S, Ciribilli Y. MYC as a multifaceted regulator of tumor microenvironment leading to metastasis. Int J Mol Sci. 2020;21(20):7710.
Ramaswamy S, Ross KN, Lander ES, Golub TR. A molecular signature of metastasis in primary solid tumors. Nat Genet. 2003;33:49–54.
Arias CR, Yeh HY, Soo VW. Biomarker identification for prostate cancer and lymph node metastasis from microarray data and protein interaction network using gene prioritization method. ScientificWorldJournal. 2012;2012:842727.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646–74.
Collado M, Serrano M. Senescence in tumours: evidence from mice and humans. Nat Rev Cancer. 2010;10:51–7.
Alimonti A, Nardella C, Chen Z, Clohessy JG, Carracedo A, Trotman LC, Cheng K, Varmeh S, Kozma SC, Thomas G, et al. A novel type of cellular senescence that can be enhanced in mouse models and human tumor xenografts to suppress prostate tumorigenesis. J Clin Invest. 2010;120:681–93.
Sharpless NE, Sherr CJ. Forging a signature of in vivo senescence. Nat Rev Cancer. 2015;15:397–408.
Guccini I, Revandkar A, D’Ambrosio M, Colucci M, Pasquini E, Mosole S, Troiani M, Brina D, Sheibani-Tezerji R, Elia AR, et al. Senescence reprogramming by TIMP1 deficiency promotes prostate cancer metastasis. Cancer Cell. 2021;39:68-82 e69.
Zhang J, Bajari R, Andric D, Gerthoffert F, Lepsa A, Nahal-Bose H, Stein LD, Ferretti V. The international cancer genome consortium data portal. Nat Biotechnol. 2019;37:367–9.
Bianco-Miotto T, Chiam K, Buchanan G, Jindal S, Day TK, Thomas M, Pickering MA, O’Loughlin MA, Ryan NK, Raymond WA, et al. Global levels of specific histone modifications and an epigenetic gene signature predict prostate cancer progression and development. Cancer Epidemiol Biomarkers Prev. 2010;19:2611–22.
Dawkins JB, Wang J, Maniati E, Heward JA, Koniali L, Kocher HM, Martin SA, Chelala C, Balkwill FR, Fitzgibbon J, Grose RP. Reduced expression of histone methyltransferases KMT2C and KMT2D correlates with improved outcome in pancreatic ductal adenocarcinoma. Cancer Res. 2016;76:4861–71.
Sato K, Akimoto K. Expression levels of KMT2C and SLC20A1 identified by information-theoretical analysis are powerful prognostic biomarkers in estrogen receptor-positive breast cancer. Clin Breast Cancer. 2017;17:e135–42.
Gala K, Li Q, Sinha A, Razavi P, Dorso M, Sanchez-Vega F, Chung YR, Hendrickson R, Hsieh JJ, Berger M, et al. KMT2C mediates the estrogen dependence of breast cancer through regulation of ERalpha enhancer function. Oncogene. 2018;37:4692–710.
Toska E, Osmanbeyoglu HU, Castel P, Chan C, Hendrickson RC, Elkabets M, Dickler MN, Scaltriti M, Leslie CS, Armstrong SA, Baselga J. PI3K pathway regulates ER-dependent transcription in breast cancer through the epigenetic regulator KMT2D. Science. 2017;355:1324–30.
Carver BS, Tran J, Gopalan A, Chen Z, Shaikh S, Carracedo A, Alimonti A, Nardella C, Varmeh S, Scardino PT, et al. Aberrant ERG expression cooperates with loss of PTEN to promote cancer progression in the prostate. Nat Genet. 2009;41:619–24.
Stanbrough M, Leav I, Kwan PW, Bubley GJ, Balk SP. Prostatic intraepithelial neoplasia in mice expressing an androgen receptor transgene in prostate epithelium. Proc Natl Acad Sci U S A. 2001;98:10823–8.
Ellwood-Yen K, Graeber TG, Wongvipat J, Iruela-Arispe ML, Zhang J, Matusik R, Thomas GV, Sawyers CL. Myc-driven murine prostate cancer shares molecular features with human prostate tumors. Cancer Cell. 2003;4:223–38.
Ding Z, Wu CJ, Chu GC, Xiao Y, Ho D, Zhang J, Perry SR, Labrot ES, Wu X, Lis R, et al. SMAD4-dependent barrier constrains prostate cancer growth and metastatic progression. Nature. 2011;470:269–73.
Wang J, Kobayashi T, Floc’h N, Kinkade CW, Aytes A, Dankort D, Lefebvre C, Mitrofanova A, Cardiff RD, McMahon M, et al. B-Raf activation cooperates with PTEN loss to drive c-Myc expression in advanced prostate cancer. Cancer Res. 2012;72:4765–76.
Arriaga JM, Panja S, Alshalalfa M, Zhao J, Zou M, Giacobbe A, Madubata CJ, Yeji Kim J, Rodriguez A, Coleman I, et al. A MYC and RAS co-activation signature in localized prostate cancer drives bone metastasis and castration resistance. Nat Cancer. 2020;1:1082–96.
Hubbard GK, Mutton LN, Khalili M, McMullin RP, Hicks JL, Bianchi-Frias D, Horn LA, Kulac I, Moubarek MS, Nelson PS, et al. Combined MYC activation and pten loss are sufficient to create genomic instability and lethal metastatic prostate cancer. Cancer Res. 2016;76:283–92.
Faget DV, Ren Q, Stewart SA. Unmasking senescence: context-dependent effects of SASP in cancer. Nat Rev Cancer. 2019;19:439–53.
Chen Z, Carracedo A, Lin HK, Koutcher JA, Behrendt N, Egia A, Alimonti A, Carver BS, Gerald W, Teruya-Feldstein J, et al. Differential p53-independent outcomes of p19(Arf) loss in oncogenesis. Sci Signal. 2009;2:ra44.
Madden SK, de Araujo AD, Gerhardt M, Fairlie DP, Mason JM. Taking the Myc out of cancer: toward therapeutic strategies to directly inhibit c-Myc. Mol Cancer. 2021;20:3.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6:p11.
Bankhead P, Loughrey MB, Fernandez JA, Dombrowski Y, McArt DG, Dunne PD, McQuaid S, Gray RT, Murray LJ, Coleman HG, et al. QuPath: Open source software for digital pathology image analysis. Sci Rep. 2017;7:16878.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Drost J, Karthaus WR, Gao D, Driehuis E, Sawyers CL, Chen Y, Clevers H. Organoid culture systems for prostate epithelial and cancer tissue. Nat Protoc. 2016;11:347–58.
Babraham Bioinformatics [https://www.bioinformatics.babraham.ac.uk/projects/fastqc/]
Ewels P, Magnusson M, Lundin S, Kaller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32:3047–8.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Howe KL, Achuthan P, Allen J, Allen J, Alvarez-Jarreta J, Amode MR, Armean IM, Azov AG, Bennett R, Bhai J, et al. Ensembl 2021. Nucleic Acids Res. 2021;49:D884–91.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.
Wang L, Wang S, Li W. RSeQC: quality control of RNA-seq experiments. Bioinformatics. 2012;28:2184–5.
Broad Institute. [http://broadinstitute.github.io/picard/].
Wingett SW, Andrews S. Fast Screen: A tool for multi-genome mapping and quality control. F1000Res. 2018;7:1338.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.
R: A language and environment for statistical computing. R Foundation for Statistical Computing [https://www.R-project.org/]
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.
pheatmap: Pretty Heatmaps. R package version 1.0.12. [https://CRAN.R-project.org/package=pheatmap]
Gennady K. Fast gene set enrichment analysis. Biorxiv. 2021. https://doi.org/10.1101/060012.
Durinck S, Spellman PT, Birney E, Huber W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protoc. 2009;4:1184–91.
Colaprico A, Silva TC, Olsen C, Garofano L, Cava C, Garolini D, Sabedot TS, Malta TM, Pagnotta SM, Castiglioni I, et al. TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 2016;44:e71.
We thank Jae W. Lee (Oregon Health and Sciences University, Portland, OR) for generously providing the Kmt2ctm1.1Jwle mouse strain. Core Facility Bioinformatics and Core Facility Genomics of CEITEC Masaryk University supported by the NCLG research infrastructure (LM2018132 funded by MEYS CR) is gratefully acknowledged for the obtaining of the sequencing data presented in this paper.
LK was supported by the BM Fonds (15142), the Margaretha Hehberger Stiftung (15142), and the Christian-Doppler Lab for Applied Metabolomics. LK and TL were supported by the Austrian Science Fund, FWF grant P 26011 and P 29251 and by the COMET Competence Center CBmed—Center for Biomarker Research in Medicine (FA791A0906.FFG). The COMET Competence Center CBmed is funded by the Austrian Federal Ministry for Transport, Innovation and Technology (BMVIT); the Austrian Federal Ministry for Digital and Economic Affairs (BMDW); Land Steiermark (Department 12, Business and Innovation); the Styrian Business Promotion Agency (SFG); and the Vienna Business Agency (WAW). The COMET program is executed by the FFG. LK was recipient of funds from a European Union Horizon 2020 Marie Sklodowska-Curie Innovative Training Network (ITNETN) grant, award n. 675712.
Ethics approval and consent to participate
Institutional Review Board Statement: The use of clinical material was approved by the Research Ethics Committee of the Medical University Vienna, Austria (1877/2016) and conducted in adherence to the Declaration of Helsinki protocols. Patient consent was waived due to the completely anonymized, retrospective nature of the study.
All animal studies were reviewed and approved by the Federal Ministry Republic of Austria for Education, Science and Research and conducted according to regulatory standards (BMBWF GZ 66.009/0135-WF/V/3b/2016).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Figure 1. related to Figure 1. Supplementary Figure 2. related to Figure 2. Supplementary Figure 3. related to Figure 3. Supplementary Figure 4. related to Figure 4. Supplementary Figure 5. related to Figure 6.
Genotyping Primers. Supplementary Table 2. Genotyping Protocols. Supplementary Table 3. RT-qPCR Primers. Supplementary Table 4. Geneset “PCa LN Metastasis UP”. Supplementary Table 5. Geneset “core SASP of PICS”.
About this article
Cite this article
Limberger, T., Schlederer, M., Trachtová, K. et al. KMT2C methyltransferase domain regulated INK4A expression suppresses prostate cancer metastasis. Mol Cancer 21, 89 (2022). https://doi.org/10.1186/s12943-022-01542-8
- Prostate cancer