Frequent genetic aberrations in the cell cycle related genes in mucosal melanoma indicate the potential for targeted therapy
Journal of Translational Medicine volume 17, Article number: 245 (2019)
Melanoma is one of the most aggressive cancers with extremely poor prognosis, and the median survival time for stage IV patients is approximately 6 to 8 months. Unlike cutaneous melanoma, mucosal melanoma is a rare melanoma subtype among Caucasian patients but its incidence remains as high as 22.6% among Chinese patients. Screening specific genetic variations is the guideline to select targeted drugs for the treatment of advanced melanoma, whereas the genetic variation spectrum and potential therapeutic targets for mucosal melanoma are largely unclear. It is urgent to identify promising genetic variants for mucosal melanoma so as to develop effective targeted therapies for this disease.
Tumor samples from 213 Chinese mucosal melanoma patients were involved in this study. P16INK4a/Cyclin D1/CDK4 copy number was examined using the QuantiGene Plex DNA assay and the correlation between abnormal copy number and clinicopathological parameters was analyzed. Patient-derived xenograft models (PDX) were performed to detect the effects of CDK4/6 inhibitors on the proliferation of mucosal melanoma cells with altered copy number of CDK4 pathway (CDK4, Cyclin D1 and P16INK4a). The molecular mechanisms of CDK4/6 inhibitors on the proliferation of mucosal melanoma were analyzed by RNAseq.
Among the 213 samples, the amplification rate of CDK4 and CCND1 was 47.0% and 27.7%, respectively, and the deletion rate of P16INK4a was 57.7%. Patients with more than one genetic abnormality were up to 81.7%. CDK4 pathway gene copy number variation was not associated with the prognosis of patients with mucosal melanoma (P > 0.05). Drug sensitivity tests showed that AT7519, a broad-spectrum CDK inhibitor, and PD0332991, a specific CDK4/6 inhibitor, exhibited higher inhibitory effect on CDK4 signaling pathway abnormal mucosal melanoma cells-derived PDX tumors growth than CDK4 signaling pathway normal ones. RNA-seq analysis showed that CDK4 inhibitors may affect tumor proliferation through multiple signaling pathways.
Abnormal copy number of cell cycle related genes is frequently found in mucosal melanoma. CDK4/6 inhibitors significantly suppress the PDX tumor growth with abnormal CDK4 pathway. CDK4 signaling variations predict the effectiveness of CDK4 inhibitors in mucosal melanoma.
Melanoma is one of the most aggressive cancers with extremely poor prognosis. Approximately 100,000 new cases of cutaneous melanoma (CM) are diagnosed and 7000 death occurs in the US in 2019 . In China, the incidence of malignant melanoma is gradually growing, with an average of more than 20,000 new cases . Because of the increased number of melanoma patients, more attention has been paid to the prevention and treatment of melanoma in China. Mucosal melanoma (MM), an aggressive subtype of melanoma, is extremely rare among Caucasian melanoma patients, with the incidence rate less than 2% . However, the incidence of MM is very high, accounting for 22.6% among Chinese melanoma patients . Compared with other subtypes, MM is characterized by occult location, late stage of initial diagnosis, high recurrence and metastasis rate. Due to the difficulty of diagnosis and treatment, the prognosis of MM is significantly poor, with a 14% 5-year survival. For Stage IV patients, the average survival time was only 6–8 months, and the 5-year survival rate was less than 5% [4, 5]. It is important to develop promising diagnostic biomarkers and effective treatment strategies for MM.
Genomic analysis shows that the major genetic variation of MM is the large-scale structural amplification or deletion of chromatin, distinct from the point mutation hotspot in CM [6,7,8,9]. Therefore, MM is different from CM in its pathogenesis, diagnosis, treatment and response to drugs. For instance, BRAF active mutation is common (> 50%) in CM, whereas the mutation rate is only about 10% in MM [6, 10]. Moreover, the effect of such targeted therapy on the treatment for MM was unclear due to the limited numbers of MM patients included in clinical trials. Copy number aberrations study may be helpful to investigate the pathogenesis of MM and to screen effective drug targets.
Abnormal cell cycle progression, caused by the mutation or amplification of CDK or cell cycle families, is one of the hallmarks of tumor cells . CDKN2A/CDK4/6/CCND1 complex is the essential regulator of cell cycle. CDK4/6/CCND1 complex phosphorylates the retinoblastoma (Rb) and inhibits its activity, leading to the transition of G1 to S phase. By contrast, CDKN2A negatively regulates the progression of cell cycle. The combination of Palbociclib (a specific inhibitor of CDK4/6) and MLN0128 (an inhibitor of mTOR) significantly inhibits the proliferation of patient derived xenograft (PDX) model of ER-positive breast cancer . High throughput sequencing in nasopharyngeal carcinoma reveals an increase in CCND1 copy number combined with CDKN2A gene deletion. Palbociclib significantly inhibits CDK4 signaling pathway activation in patient derived xenograft (PDX) model . Palbociclib has antitumor activity for NRAS-mutant melanomas in a preclinical mouse model, indicating the CDK4 pathway as a potential therapeutic target . A large cohort analysis in 2017 found that CDK4 gene amplification was higher in acral melanoma (AM), suggesting that CDK4 may be a therapeutic target for AM [15,16,17].
It is not clear whether there are variations in cell cycle related signaling pathways in MM and whether they can be used as effective therapeutic targets. To address this question, we collected 213 MM samples and investigated the status of CDK4 signaling related genes (CDK4, CCND1 and CDKN2A) and its relationship with clinical pathologic parameters. We also explored the inhibitory effects of CDK4/6 inhibitors on PDX tumor growth containing normal or aberrant CDK4 pathway. The results of this study may provide a new idea for the clinical treatment of MM.
Materials and methods
Patients and tissue samples
This study involved samples from primary lesions of 213 MM patients, who were hospitalized during January 2007 and October 2018 at the Peking Cancer Hospital & Institute. We obtained all the clinical and pathological data by medical record review, including age, gender, ulceration, depth of invasion, anatomic site, treatment, tumor-node-metastases stage, mutation status of therapeutic targets (such as KIT, BRAF and NRAS), follow-up time and survival (follow-up persisted until June 2018, or until the missing of follow-up or death of patients). An informed consent was obtained from all participants that were enrolled in clinical trials. This study was approved by the Medical Ethics Committee of the Beijing Cancer Hospital & Institute and was conducted according to the Declaration of Helsinki Principles.
QuantiGenePlex DNA assay
Tissue homogenates were prepared according to the protocols in the user manual of QuantiGene Sample Processing Kit for Formalin-Fixed, Paraffin-Embedded Tissues (FFPE; Panomics of Affymetrics, Santa Clara, CA). Briefly, 5 to 8 pieces of deparaffinized sections (4 to 10 μm) were incubated with 150 μl homogenizing solution supplemented with 1.5 μl of proteinase K (50 μg/μl) at 65 °C for 6 h. The tissue homogenates were separated from debris by brief centrifugation and transferred to a new tube.
The branched DNA (bDNA) assay was performed according to the procedures described in the user manual of QuantiGenePlex DNA Assay (Panomics). Briefly, the homogenate DNA was sheared by the Covaris S2 (Covaris, Woburn, MA) with the following settings: duty cycle 5%, intensity 3, cycles/burst 200, 80 s. For each assay well, 40 μl homogenate was denatured with 2.5 M NaOH (final concentration 0.18 M) in the presence of DNA probe. Neutralized tissue homogenate was transferred to each well of the Hybridization Plate containing Working Bead Mix. All the samples were in duplicates. Hybridization Plate was sealed and incubated at 54 °C ± 1 °C in a shaking incubator (600 rpm) for 18–22 h. The unbound samples were washed away using the Bio-plex pro II wash station (Bio-Rad, Hercules, CA). Then the beads were sequentially hybridized with the DNA Pre-Amplifier, DNA Amplifier, Label Probe and SAPE (Streptavidin-conjugated R-phycoerythrin). Fluorescence intensities were measured by the Bio-plex 100 system (Bio-Rad).
The mean fluorescence intensities of the duplicates were calculated for all genes. The background values were subtracted from each probe set signal. Values of tested genes were normalized to the geometric means of Rpph1, Rpp30 and Rplp0. For each test sample, normalized signal was divided by the reference DNA sample (G1521, Promega, Madison, WI) for each test gene, and the values were multiplied by the known copy number (usually 2 copies) of each gene in the reference genome. The Cdk4, Ccnd1 and P16INK4a copy numbers for the samples were calculated by dividing the sample values tothe control values: No gain referred to copy number ≤ 2, gain referred to copy number > 2, loss referred to copy number < 1.
DNA preparation and TaqMan copy number assays
Genomic DNA was extracted from FFPE sections using a QIAamp DNA FFPE Tissue Kit (Qiagen, Hilden, Germany). To validate the results of QuantiGenePlex DNA Assay, the copy numbers of Cdk4, Ccnd1 and P16INK4a were further quantified by TaqMan Copy Number Assays (Applied Biosystems of ThermoFisher, Waltham, MA). A TaqMan probe targeted on the Rnasep gene was used as a control. Quantitative real-time PCR was performed using the ABI 7500 FAST real-time PCR system (Applied Biosystems). Copy numbers were then determined by CopyCaller v2.0 software (Applied Biosystems) using the comparative Ct (ΔΔCt) method.
Immunohistochemistry analyses were performed using antibodies against Ki67 (dilution 1:400) (Abcam, Cambridge, UK) as described (11, 17). The Ki-67+ cells under 5 random fields were counted and were presented as mean ± SD of three sections.
The HMVII cell line was obtained from Sigma-Aldirch and was cultured at 37 °C in Ham’s F-10 medium supplemented with 1% penicillin and streptomycin (Invitrogen) and 10% fetal bovine serum (HyClone of GE Healthcare, Logan, UT). The GAK cell line was obtained from JCRB Cell Bank and was cultured at 37 °C in Ham’s F-12 medium supplemented with 1% penicillin and streptomycin (Invitrogen) and 10% fetal bovine serum (HyClone of GE Healthcare, Logan, UT).
Cell proliferation assays
CDK4/6 inhibitors including PD0332991 (#S1116), and pan-CDK inhibitor AT7519 (#S1524) were purchased from Selleck Chemicals (Houston, TX). All inhibitors were dissolved at 10 mM in dimethylsulfoxide (DMSO) as stock solutions. After treatment with various concentrations of inhibitors or DMSO for 24 h, proliferation of the cells was evaluated using the Cell Titer-Glo Luminescent Cell Viability Assay (Promega) according to the instructions.
Patient-derived xenograft (PDX) model and treatment
Fragments of patient-derived MM tissues bearing typical CDK4 pathway aberrations were cut into fragments and then subcutaneously inoculated into a 6 week-old NOD/SCID (non-obese diabetic and severe combined immunodeficiency) female mouse (4–6 week-old, 18–22 g-weight) to establish the PDX model. The established PDX model was called passage 0 (P0). When the tumor size reached approximately 500 mm3, the mice were sacrificed, and tumor tissues were separated and re-inoculated into new mice to obtain the subsequent passages called P1, P2, P3 and forth. 10 PDX models containing typical CDK4 pathway aberrations (Table 1) were finally established.
Mice (P2) were randomized (treatment arm versus control arm; n = 4) and treated with control (sodium lactate buffer, pH 4.5) or PD0332991 and with control (saline solution, pH 4.0) or AT7519. For PD0332991 treatment, mice received PD0332991 (50 mg/kg in pH 4.5 sodium lactate buffer) via oral gavage daily. For AT7519 treatment, mice received AT7519 (12 mg/kg in saline solution) via intraperitoneal injection daily. Tumor sizes were measured every 3 days and tumor volume was calculated using the formula: volume = length*width2/2. Percentage of tumor volume on day of treatment (baseline volume) was used as the end-point of study. The treatment lasted for 14 days, after which the mice were sacrificed and the tumors were fixed in 10% formalin for histological and immunohistological analysis. All animal care and experimental procedures were performed in consistent with the Animal Care Ethics approved by the Medical Ethics Committee of the Beijing Cancer Hospital & Institute.
RNA-sequencing was performed at the Shanghai Biotechnology Company. cDNA library was built according to the standard manufacturer’s protocol. Paired-end 2 × 100 bp read sequencing was performed using the Illumina HiSeq X-ten (Illumina, USA). The FASTX-Toolkit (v0.0.13) was used to trim low-quality bases. High-quality reads were aligned to the human GRCh38 reference genome with two mismatches using spliced mapping alignment in Hisat2 (version: 2.0.4). After genome mapping, Stringtie (version: 1.3.0) was run with a reference annotation to generate FPKM values for known gene models. Differentially expressed genes were identified using edgeR. The P value significance threshold in multiple tests was set by the false discovery rate (FDR). The fold-changes were also estimated according to the FPKM (Fragments Per Kilobase of exon model per Million mapped read) in each sample. The differentially expressed genes were selected using the following filter criteria: FDR ≤ 0.05 and fold-change ≥ 2.
Statistical analyses were performed using SPSS 22.0 software. Continuous data such as age was described using mean ± SD for normally distributed data. The correlations between aberration status and clinical parameters were evaluated by Chi square test or Fisher’s exact test. Kaplan–Meier estimates of time-to-event overall survival (OS) and follow-up time were calculated. Log-rank tests were used to estimate the statistical significance between the time-dependent outcomes of OS. All statistical analyses were two-sided, and P < 0.05 was considered as statistically significant.
Aberrations of Cdk4, Ccnd1 and P16INK4a in MM
Firstly, we detected the copy number variation of CDK4, CCND1 and P16INK4a genes in paraffin sections of 213 cases of MMs through QGP method (Table 2). Among the 213 samples, 100 cases (47.0%), 59 cases (27.7%) and 123 cases (57.7%) showed Cdk4 gain, Ccnd1 gain, and P16INK4a loss, respectively. Moreover, 36.6% of MMs contained more than two concurrent aberrations, and 8.0% of MMs contained three aberrations. The overall frequency of MM containing any CNV (≥ 1 CNV) was 81.7%, whereas 39 cases harbored no CNV aberrations in these three genes. 38 cases had Ccnd1 loss, and one case harbored P16INK4a gain (Table 2). Collectively, Cdk4 gain, Ccnd1 gain and/or P16INK4a loss were observed in most MMs.
Correlation of CDK4 pathway aberrations to other driver genes
We also detected the mutation status of c-Kit (exon 9, 11, 13, 17 and 18), BRAF (exon 11 and 15) and NRAS (exon 1 and 2) by Sanger sequencing. The result showed that the mutation rate of c-Kit, BRAF and NRAS was 7.4%, 5.2% and 15.7% in MM patients carrying Cdk4 gain, respectively. For MM patients carrying Ccnd1 gain, the mutation rate of c-Kit, BRAF and NRAS was 1.8%, 8.9%, 10.4%, respectively. For MM patients carrying P16INK4a, the mutation rate of c-Kit, BRAF and NRAS was 9.6%, 4.3%, 10.5%, respectively (Table 2).
Correlation of Anatomic site to TNM stage of mucosal melanoma
We next explored the anatomic sites of MM in our 213 MM patient’s cohort. The prevalent anatomic sites were head and neck (44.2%), anorectum (21.1%), genitourinary (25.8%) and oesophagus (8.9%) (Table 4). The proportion of patients with TNM stage I, II, III, and IV diseases were 2.8%, 53.5%, 28.6% and 15.1%, respectively. The percentages of patients with stage IV of MM from head and neck were significantly lesser than other anatomic sites (P < 0.001).
Correlation of CDK4 pathway aberrations to clinicopathological features
Last, we investigated the correlations between the aberration of Cdk4, Ccnd1 and P16INK4a and the clinicopathological features of MM. Statistical analysis was divided into 8 groups: Cdk4 aberration, Ccnd1 aberration, P16INK4a aberration, ≥ 1 CNV (copy number variation), Cdk4 gain + Ccnd1gain, Cdk4 gain + P16INK4a loss, Ccnd1gain + P16INK4a loss and Cdk4 gain + Ccnd1gain + P16INK4a loss.
In our cohort, both age, gender, anatomic site and follow-up time were not significantly different between patients with or without any CNVs for Cdk4, Ccnd1, P16INK4a or other indicated stochastic combinations. However, the other clinical features were at least significantly associated with a molecular variant of the CDK4 signaling pathway: compared with P16INK4a loss, patients with normal were more likely to ulcerate;The depth of invasion of P16INK4a normal group tended to be deeper than that of group with P16INK4a loss, while there was no statistical difference (P = 0.053). The percentages of patients with stage I–IV of MM were significantly different between patients with CDK4 gain and CDK4 normal (P = 0.004; Table 3). The TNM stage of ≥ 1 CNV group tended to be greater than that of group without any CDK4 pathway aberrations, while there was no statistical difference (P = 0.059);Similar result was observed in the group of CDK4 gain + P16INK4a loss (P = 0.088). The treatment groups for MM patients was significantly different between patients with CDK4 gain and CDK4 normal (P = 0.003; Additional file 1: Table S1).
The overall survival of MM patients with Cdk4 gain, Ccnd1 gain, P16INK4a loss or other indicated stochastic combinations were comparable. Univariate Cox analysis suggested thatCdk4 gain, Ccnd1 gain, P16INK4a loss or other stochastic combinations might not be of prognostic significance for melanoma patients (Table 3).
Sensitivity of MM cells to CDK4/6 inhibitors
To evaluate the effect of CDK4/6 inhibition on the proliferation of MM cells, GAK and HMV II cell lines were treated with different concentrations of pan-CDK inhibitor AT7519 and CDK4/6 specific inhibitor PD0332991. The genetic variations of key genes for these cells are listed in Table 4. Both pan-CDK inhibitor AT7519 and CDK4/6 specific inhibitor PD0332991 significantly inhibited the viability of GAK and HMV II cells. For AT7519, the inhibitory rate in GAK and HMV II cells was obvious when higher than 4 µ mol/l. For PD0332991, GAK and HMV II were strikingly sensitive at a concentration higher than 10 µ mol/l (Fig. 1).
Sensitivity of PDX models to CDK4/6 inhibitors
To analyze the sensitivity of MM containing typical CDK4 pathway aberrations to CDK4/6 inhibitors, we selected 10 PDX models with 5 different CDK4 pathway aberrations. The basic information of the PDX models is shown in Table 1. Two models (PDX-001 and PDX-002) are CDK4 pathway normal. Four models (PDX-003, PDX-004, PDX-005 and PDX-006) are Cdk4 gain + Ccnd1 gain + P16INK4a loss. Two models (PDX-007 and PDX-008) are Cdk4 gain + Ccnd1 gain. One model (PDX-009) is Ccnd1 gain + P16INK4a loss. One model (PDX-010) is Ccnd1 gain.
Then we treated the PDX models with AT7519 and PD0332991. As compared to the vehicle-treated group, AT7519 and PD0332991 showed no inhibitory effect on tumor growth in PDX-001 and PDX-002 model. Interestingly, AT7519 and PD0332991 significantly retarded the tumor growth of PDX-003 to PDX-010 (Fig. 2). To confirm our findings, we performed immunohistochemical staining of Ki-67 in these tumors. Consistently, the number of Ki-67+ cells was significantly decreased after AT7519 and PD0332991 treatments in PDX models with CDK4 pathway aberration (Fig. 3). Taken together, these data indicate that CDK4 aberration dictates the sensitivity of MM-derived PDX tumor growth to CDK inhibitors (Table 5).
Gene expression profile of PDX005 tumor tissue to CDK4/6 inhibitors
To identify the molecular alterations after CDK4/6 inhibitors treatment, mRNA sequencing was performed on PDX tumor cells isolated from three paired tissue samples. On average, 7.6 million reads (between 7.0 and 8.6 million per sample) were obtained using Illumina HiSeq 2500 platform of which 97% had high quality scores (≥ Q20). Compared with vehicle group, AT7519 treated PDX-005 tumors had 1345 DEGs (FC ≥ 2, adjusted P-values < 0.05), among which 725 genes were down-regulated and 620 genes were up-regulated. In PD0332991-treated tumors, a total of 919 DEGs (FC ≥ 2, adjusted P-values < 0.05) were observed, among which 648 genes were down-regulated and 271 genes were up-regulate. The volcanic map of down- and up-regulated genes was listed in Fig. 4.
Functional annotation analysis of all DEGs utilizing: Profiler software revealed “Ribosome” (58 genes, P < 7.8E−24), “Oxidative phosphorylation” (32 genes, P < 4.8E−08), “Focal adhesion” (29 genes, P < 0.001) and “Antigen processing and presentation” (15 genes, P < 0.002) were the most enriched KEGG pathways for AT7519 vs vehicle. For PD0332991 vs vehicle group, “Cell cycle” (25 genes, P < 7.5E−9), “DNA replication” (9 genes, P-values 3.1E−05), “Focal adhesion” (25 genes, P < 5.6E−05) and “p53 signaling pathway” (11 genes, P < 0.0004) were the most enriched KEGG pathways. The top 10 pathways were showed in Fig. 5.
Melanoma is the fifth most common cancer in the Caucasian population . CM is the major subtype melanoma in Caucasian population. Epidemiological statistics showed more than 5000 new cases of MM in China every year, and the incidence of MM is increasing . Unlike CM, MM exhibits distinct biological and clinical features. Targeted therapies, such as BRAF inhibitors and c-KIT inhibitors, have greatly revolutionized the treatment of CM. However, the effect of such targeted therapy on the treatment for MM remained obscure due to following reasons: 1, the number of MM patients included in clinical trials is very limited. 2, the mutation rates of the BRAF and c-KIT genes were only about 13.9 and 9.6%, respectively, in MM patients [18, 19]. Thus, potential targets should be identified for MM patients for developing new and effective targeted therapies.
Copy number variation of CDK signaling pathway genes plays an important role in many tumors. Reducing the activity of CDK signaling pathway can significantly inhibit the growth of tumor. CDK4/6 and CCND1 genes were amplified while CDKN2A, CDKN2B and CDKN1B were deleted in primary and metastatic breast cancers . In patients with type 1 neurofibromatosis associated with breast cancer, the expansion of CDK4 copy number may increase the expression of the NF1 gene and then up-regulate the expression of Her2 gene in breast cancer cells . CDK4 is amplified and its protein expression is increased in esophageal squamous cell carcinoma (ESCC), and both of them are associated with the poor survival of ESCC. CDK4 silencing or PD0332991 treatment significantly inhibits the proliferation of esophageal cancer cells . Similar results are observed that CDK4/6 inhibitors suppress the proliferation of thyroid cancer cells . In addition, research for malignant glioma, neuroblastoma and malignant peripheral nerve sheath tumors demonstrated that CDK4 copy number was amplified and was associated with tumor prognosis [23,24,25]. Moreover, our previous study of 514 cases of acral melanoma showed that 87% of patients had at least one copy number variant of the CDK4 signaling pathway genes. Cdk4 gain, Ccnd1 gain and the combination of Cdk4 gain and Ccnd1 gain were associated with the poor prognosis of acral melanoma, respectively . In this study, copy number amplification of CDK4 and CCND1 genes and copy number deletion of CDKN2A gene were observed in most MM cases, suggesting that CDK4 signaling pathway aberrations maybe the driver for MM. Further analysis showed that CNV was significantly associated with the TNM stage in MM. However, the aberration of CDK4 signaling pathway genes was not significantly associated with the OS of MM patients. Therefore, it is necessary to expand the sample for further verification. Richard et al. found that CDKN2A copy number loss was a frequent event in patients with CM, and the prognosis is worse in patients with increased Ccnd1 copy number and Cdk4 copy number. However, after stratification of the status of NRAS and Braf mutations in patients, copy number variation in the CDK4 signaling pathway genes was not significantly associated with the survival of CM . This is similar to our results because the mutation rates of NRAS and Braf is extremely low in MM patients, and the patients analyzed tend to be independent of NRAS and BRAF mutation.
Previous studies have shown that CDK inhibitors can inhibit the growth of a variety of tumors. CDK4/6 inhibitor Palbociclib combined with mTOR inhibitor MLN0128 suppresses the proliferation of ER negative breast cancer cells and Glioma [12, 27]. Similarly, triple-negative breast cancers with MT4-MMP, EGFR and RB-positive are sensitive to Erlotinib in combination with Palbociclib in PDX model . Our previous study found that both CDK inhibitors and specific CDK4/6 inhibitors effectively inhibited the tumor growth in PDX models of acral melanoma containing CDK4 pathway aberration . Zhou et al. found that CDK4 gene was amplified in 65 MM patients and Palbociclib, an inhibitor of CDK4/6, effectively blunted the tumors harboring Cdk4 copy number gain in a PDX model . Here, we found that both AT7519, a pan CDK inhibitor, and PD0332991, a selective CDK4/6 inhibitor, significantly inhibited the proliferation of MM cell lines in vitro. In vivo, both inhibitors obviously reduced the growth of tumors which harbored abnormal CDK4 signaling pathway in PDX models. However, tumors with normal CDK4 signaling pathway exhibited minimal sensitivity to both inhibitors. These results suggest that CDK4 aberration is an indicator for CDK4/6 inhibitors applied in MM. We also found in clinic that Palbociclib, an inhibitor of CDK4/6, can prolong the survival of MM patients with copy number variations of CDK4 pathway . Therefore, CDK4 pathway genes copy number variation is a therapeutic target for MM.
To explore the potential molecular mechanisms, we subjected a PDX model (PDX-005) containing Cdk4 gain + Ccnd1 gain + P16INK4a loss to RNA sequencing. Analysis of GO and KEGG revealed significant changes in several signaling pathways, including cell cycle and immunity. It is worth noting that in the enrichment analysis of the biological process, cellular component, molecular function and KEGG Pathway, the immune-related signaling pathways were significantly altered during the four above-mentioned analyses. These results suggest that CDK4/6 inhibitors may play an important role in the development of MM by regulating the cellular immune signaling pathway. Alterations in immune-related signaling pathways suggest that CDK4/6 inhibitors may affect the cellular immune system, thus affecting the efficacy of immunotherapies to tumor. Immunotherapy, including PD-1 antibody, has revolutionized the treatment of many tumors. In breast cancer, Shom et al. found that CDK4/6 inhibitors not only induce tumor cell cycle arrest, but also promote anti-tumor immunity in two main way: 1, CDK4/6 inhibitors alter tumor cell expression of endogenous retroviral components, hence increasing intracellular levels of double-stranded RNA. This in turn stimulates production of type III interferons and therefore enhances tumor antigen presentation. 2, CDK4/6 inhibitors markedly check the proliferation of regulatory T cells . In general, the above-mentioned results suggest that CDK4/6 inhibitors are common in inducing antitumor immunity.
In this study, we found that the copy number of the components of CDK4 signaling pathway is altered in MM. Inhibition of CDK signaling pathway and CDK4 signaling pathway alone obviously suppresses the proliferation of MM cells and the tumor growth in PDX models harboring CDK4 pathway abnormity. This study provides theoretical basis for targeting CDK4 pathway in MM.
Availability of data and materials
All the data and materials supporting the conclusion were included in the main paper.
copy number variation
false discovery rate
Formalin-Fixed, Paraffin-Embedded Tissues
Fragments Per Kilobase of exon model per Million mapped read
Kyoto Encyclopedia of Genes and Genomes
non-obese diabetic and severe combined immunodeficiency
patient-derived xenograft model
QuantiGenePlex DNA assay
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69:7–34.
Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66:115–32.
McLaughlin CC, Wu XC, Jemal A, Martin HJ, Roche LM, Chen VW. Incidence of noncutaneous melanomas in the U.S. Cancer. 2005;103:1000–7.
Chi Z, Li S, Sheng X, Si L, Cui C, Han M, Guo J. Clinical presentation, histology, and prognoses of malignant melanoma in ethnic Chinese: a study of 522 consecutive cases. BMC Cancer. 2011;11:85.
Lian B, Si L, Cui C, Chi Z, Sheng X, Mao L, Li S, Kong Y, Tang B, Guo J. Phase II randomized trial comparing high-dose IFN-alpha2b with temozolomide plus cisplatin as systemic adjuvant therapy for resected mucosal melanoma. Clin Cancer Res. 2013;19:4488–98.
Furney SJ, Turajlic S, Stamp G, Nohadani M, Carlisle A, Thomas JM, Hayes A, Strauss D, Gore M, van den Oord J, et al. Genome sequencing of mucosal melanomas reveals that they are driven by distinct mechanisms from cutaneous melanoma. J Pathol. 2013;230:261–9.
Krauthammer M, Kong Y, Ha BH, Evans P, Bacchiocchi A, McCusker JP, Cheng E, Davis MJ, Goh G, Choi M, et al. Exome sequencing identifies recurrent somatic RAC1 mutations in melanoma. Nat Genet. 2012;44:1006–14.
Curtin JA, Fridlyand J, Kageshita T, Patel HN, Busam KJ, Kutzner H, Cho KH, Aiba S, Brocker EB, LeBoit PE, et al. Distinct sets of genetic alterations in melanoma. N Engl J Med. 2005;353:2135–47.
Curtin JA, Busam K, Pinkel D, Bastian BC. Somatic activation of KIT in distinct subtypes of melanoma. J Clin Oncol. 2006;24:4340–6.
Flaherty KT, Puzanov I, Kim KB, Ribas A, McArthur GA, Sosman JA, O’Dwyer PJ, Lee RJ, Grippo JF, Nolop K, Chapman PB. Inhibition of mutated, activated BRAF in metastatic melanoma. N Engl J Med. 2010;363:809–19.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646–74.
Yamamoto T, Kanaya N, Somlo G, Chen S. Synergistic anti-cancer activity of CDK4/6 inhibitor palbociclib and dual mTOR kinase inhibitor MLN0128 in pRb-expressing ER-negative breast cancer. Breast Cancer Res Treat. 2019;174(3):615–25.
Hsu CL, Lui KW, Chi LM, Kuo YC, Chao YK, Yeh CN, Lee LY, Huang Y, Lin TL, Huang MY, et al. Integrated genomic analyses in PDX model reveal a cyclin-dependent kinase inhibitor Palbociclib as a novel candidate drug for nasopharyngeal carcinoma. J Exp Clin Cancer Res. 2018;37:233.
Kwong LN, Costello JC, Liu H, Jiang S, Helms TL, Langsdorf AE, Jakubosky D, Genovese G, Muller FL, Jeong JH, et al. Oncogenic NRAS signaling differentially regulates survival and proliferation in melanoma. Nat Med. 2012;18:1503–10.
Hayes TK, Luo F, Cohen O, Goodale AB, Lee Y, Pantel S, Bagul M, Piccioni F, Root DE, Garraway LA, et al. A functional landscape of resistance to MEK1/2 and CDK4/6 inhibition in NRAS mutant melanoma. Cancer Res. 2019;79(9):2352–66.
Kong Y, Sheng X, Wu X, Yan J, Ma M, Yu J, Si L, Chi Z, Cui C, Dai J, et al. Frequent genetic aberrations in the CDK4 pathway in acral melanoma indicate the potential for CDK4/6 inhibitors in targeted therapy. Clin Cancer Res. 2017;23:6946–57.
Sheppard KE, McArthur GA. The cell-cycle regulator CDK4: an emerging therapeutic target in melanoma. Clin Cancer Res. 2013;19:5320–8.
Kong Y, Si L, Zhu Y, Xu X, Corless CL, Flaherty KT, Li L, Li H, Sheng X, Cui C, et al. Large-scale analysis of KIT aberrations in Chinese patients with melanoma. Clin Cancer Res. 2011;17:1684–91.
Si L, Kong Y, Xu X, Flaherty KT, Sheng X, Cui C, Chi Z, Li S, Mao L, Guo J. Prevalence of BRAF V600E mutation in Chinese melanoma patients: large scale analysis of BRAF and NRAS mutations in a 432-case cohort. Eur J Cancer. 2012;48:94–100.
Basudan A, Priedigkeit N, Hartmaier RJ, Sokol ES, Bahreini A, Watters RJ, Boisen MM, Bhargava R, Weiss KR, Karsten MM, et al. Frequent ESR1 and CDK pathway copy-number alterations in metastatic breast cancer. Mol Cancer Res. 2019;17:457–68.
Wang X, Kallionpaa RA, Gonzales PR, Chitale DA, Tousignant RN, Crowley JP, Chen Z, Yoder SJ, Blakeley JO, Acosta MT, et al. Germline and somatic NF1 alterations are linked to increased HER2 expression in breast cancer. Cancer Prev Res (Phila). 2018;11:655–64.
Ismail A, Bandla S, Reveiller M, Toia L, Zhou Z, Gooding WE, Kalatskaya I, Stein L, D’Souza M, Litle VR, et al. Early G(1) cyclin-dependent kinases as prognostic markers and potential therapeutic targets in esophageal adenocarcinoma. Clin Cancer Res. 2011;17:4513–22.
Valenciaga A, Saji M, Yu L, Zhang X, Bumrah C, Yilmaz AS, Knippler CM, Miles W, Giordano TJ, Cote GJ, Ringel MD. Transcriptional targeting of oncogene addiction in medullary thyroid cancer. JCI Insight. 2018;3:16.
Aoki K, Nakamura H, Suzuki H, Matsuo K, Kataoka K, Shimamura T, Motomura K, Ohka F, Shiina S, Yamamoto T, et al. Prognostic relevance of genetic alterations in diffuse lower-grade gliomas. Neuro Oncol. 2018;20:66–77.
Yu J, Deshmukh H, Payton JE, Dunham C, Scheithauer BW, Tihan T, Prayson RA, Guha A, Bridge JA, Ferner RE, et al. Array-based comparative genomic hybridization identifies CDK4 and FOXM1 alterations as independent predictors of survival in malignant peripheral nerve sheath tumor. Clin Cancer Res. 2011;17:1924–34.
Young RJ, Waldeck K, Martin C, Foo JH, Cameron DP, Kirby L, Do H, Mitchell C, Cullinane C, Liu W, et al. Loss of CDKN2A expression is a frequent event in primary invasive melanoma and correlates with sensitivity to the CDK4/6 inhibitor PD0332991 in melanoma cell lines. Pigment Cell Melanoma Res. 2014;27:590–600.
Olmez I, Brenneman B, Xiao A, Serbulea V, Benamar M, Zhang Y, Manigat L, Abbas T, Lee J, Nakano I, et al. Combined CDK4/6 and mTOR inhibition is synergistic against glioblastoma via multiple mechanisms. Clin Cancer Res. 2017;23:6958–68.
Foidart P, Yip C, Radermacher J, Blacher S, Lienard M, Montero-Ruiz L, Maquoi E, Montaudon E, Chateau-Joubert S, Collignon J, et al. Expression of MT4-MMP, EGFR, and RB in triple-negative breast cancer strongly sensitizes tumors to erlotinib and palbociclib combination therapy. Clin Cancer Res. 2018;25:1838–50.
Zhou R, Shi C, Tao W, Li J, Wu J, Han Y, Yang G, Gu Z, Xu S, Wang Y, et al. Analysis of mucosal melanoma whole-genome landscapes reveals clinically relevant genomic aberrations. Clin Cancer Res. 2019;25:3548–60.
Tang B, Sheng X, Kong Y, Chi Z, Si L, Cui C, Yan X, Mao L, Lian B, Li S, et al. Palbociclib for treatment of metastatic melanoma with copy number variations of CDK4 pathway: case report. Chin Clin Oncol. 2018;7:62.
Goel S, DeCristo MJ, Watt AC, BrinJones H, Sceneay J, Li BB, Khan N, Ubellacker JM, Xie S, Metzger-Filho O, et al. CDK4/6 inhibition triggers anti-tumour immunity. Nature. 2017;548:471–5.
We thank the Dr. Zhiguo Zhao for English language editing.
This work was supported by grants from Natural Science Foundation of China (81672696, 81772912), Fostering Young Scholars of Peking University Health Science Center,Clinical Medicine Plus X—Young Scholars Project (Peking University), the Fundamental Research Funds for the Central Universities, Beijing Baiqianwan Talents Project and Beijing Municipal Administration of Hospitals’ Ascent Plan (No. DFL20181101)
Ethics approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the Medical Ethics Committee of the Beijing Cancer Hospital & Institute and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. For this type of study formal consent is not required.
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.
About this article
Cite this article
Xu, L., Cheng, Z., Cui, C. et al. Frequent genetic aberrations in the cell cycle related genes in mucosal melanoma indicate the potential for targeted therapy. J Transl Med 17, 245 (2019). https://doi.org/10.1186/s12967-019-1987-z