Deficiency of PTEN and CDKN2A Tumor-Suppressor Genes in Conventional and Chondroid Chordomas: Molecular Characteristics and Clinical Relevance.

INTRODUCTION: Chordoma is a malignant tumor predominantly involving the skull base and vertebral column. This study aimed to investigate the molecular characteristics of PTEN and CDKN2A in conventional and chondroid chordomas and their correlation with clinical prognosis.
MATERIALS AND METHODS: A total of 42 patients were enrolled, including 26 patients with conventional chordoma and 16 patients with chondroid chordoma. Clinicopathological profiles and tissue specimens were collected. Gene sequencing and fluorescence in situ hybridization were performed to identify genetic alterations in the PTEN and CDKN2A genes. Immunohistochemical staining was used for semiquantitative evaluation of PTEN and CDKN2A expression.
RESULTS: Gene sequencing revealed an intron SNP (c.80-96A>G) and a missense mutation (c.10G>A; p.Gly4Arg) in the PTEN gene and a missense mutation (c.442G>A; p.Ala148Thr) in the CDKN2A gene. Loss of the PTEN locus was identified in 25 (59.5%) cases, and loss of the CDKN2A locus was found in 28 (66.7%) cases. There was no significant correlation between progression-free survival (PFS)/overall survival (OS) and loss of PTEN or CDKN2A. The patients with lower PTEN expression showed significantly shorter PFS and OS than those with higher expression, while there was no significant difference in PFS or OS between patients with lower CDKN2A expression and those with higher CDKN2A expression.
CONCLUSION: Our findings delineated the genetic landscape and expression of PTEN and CDKN2A in chordomas. PTEN expression may serve as a prognostic and predictive biomarker for chordomas.

Introduction

Chordoma is a rare, slow-growing, and locally invasive malignant skeletal tumor that occurs most commonly at the skull base and vertebral column.1 This entity is presumed to derive from remnants of the embryonic notochord, and the annual incidence is estimated to range from 0.51 to 8 per million population.2 Surgical resection is the mainstay of treatment modalities for chordoma, as chordoma is resistant to conventional chemotherapy and radiotherapy. However, despite aggressive surgical resection, chordoma still has a high local recurrence rate, and the clinical prognosis remains poor. The reported 5-year survival rate of patients with chordoma varies from 47% to 80%,2 and the overall median survival is approximately only 6 years.3 To date, the molecular and genetic mechanisms of chordoma have not been fully elucidated; in particular, there is a paucity of research on the genetic aberrations in distinct histological subtypes (conventional, chondroid, and dedifferentiated). Identification of molecular hallmarks in these subtypes may provide insight into potential novel therapies.

In 2001, Scheil et al raised the first human chordoma cell line (U‐CH1), and they found that U-CH1 had del(9) and del(10) chromosomal abnormalities.4 Subsequently, Kuźniacka et al also noted loss of 9p and chromosome 10 in chordomas.5 Chromosome 9p harbors the tumor-suppressor gene CDKN2A (cyclin-dependent kinase inhibitor 2A; located on chromosome 9p21), which encodes the cell cycle regulatory protein p16; PTEN (phosphatase and tensin homolog), another tumor-suppressor gene, is located on chromosome 10q23.3. In the last decade, accumulating evidence indicates that deficiency of PTEN and CDKN2A is common in chordomas.6–9 However, the genetic alterations of PTEN and CDKN2A in distinct histological subtypes of chordoma remain unclear, and the correlation between loss of these tumor-suppressor genes and clinical prognosis is still controversial. In this study, we aimed to investigate the molecular characteristics of PTEN and CDKN2A in conventional and chondroid chordomas and their clinical relevance.

Materials and Methods

Patients and Tumor Specimens

This study was approved by the Institutional Review Board (IRB) of Peking University Third Hospital. A total of 42 patients with chordoma were enrolled. All patients underwent surgical resection at our department, and the histopathological diagnosis was conventional chordoma in 26 patients and chondroid chordoma in 16 patients. Individual clinical and pathological profiles were collected, and histopathological sections were reviewed by two pathologists independently. Tissue specimens were obtained from all patients.

Immunohistochemical (IHC) Staining

IHC staining for PTEN and CDKN2A was performed in formalin-fixed and paraffin-embedded tissue specimens. After baking at 60°C for 30 min, the sections (4 μm) were deparaffinized with dimethylbenzene and hydrated in a graded ethanol series. Antigen retrieval was performed in Tris-EDTA buffer, and then the slides were treated with 0.3% H2O2 to block endogenous peroxide. After incubation with the primary antibody (anti-PTEN antibody, ab31392, 1:100 dilution; anti-CDKN2A antibody, ab54210, 1:500 dilution; Abcam, Cambridge, MA, USA) at 4°C overnight, visible staining was developed using the DAB Horseradish Peroxidase Color Development Kit (Beyotime Institute of Biotechnology, Shanghai, China) according to the manufacturer’s protocol. Finally, the slides were counterstained with hematoxylin, dehydrated, cleared and mounted. Breast and ovarian carcinoma tissues with known positivity for the antibodies were used as positive controls. For negative controls, the primary antibodies were replaced with phosphate-buffered saline.

IHC results were evaluated according to the scoring method proposed by Sinicrope et al,10 based on both the percentage of immunoreactive tumor cells and the intensity of staining. The percentage of tumor cells positive for PTEN or CDKN2A was scored as follows: 0 (no specific staining or staining in <10% of tumor cells); 1 (staining in 11%~25% of tumor cells); 2 (staining in 26%~50% of tumor cells); 3 (staining in 51%~75% of tumor cells); and 4 (staining in >75% of tumor cells). The intensity of immunoreactivity was graded as follows: 0 (no staining); 1 (weak staining); 2 (moderate staining); and 3 (strong staining). The final immunoreactivity score was the product of these two indices, with values ranging from 0 to 12. The median immunoreactivity score was set as the cutoff value for the classification of low and high expression.

Gene Sequencing

Genomic DNA was extracted from frozen tissue samples using a QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) or from paraffin-embedded samples using a GeneRead DNA FFPE Kit (Qiagen) according to the manufacturer’s instructions. Primers were designed using OLIGO Primer Analysis Software (version 7, Molecular Biology Insights, Cascade, CO, USA) and synthesized by Sangon Biological Technology (Shanghai, China). All annotated exons and adjacent introns of PTEN (OMIM: 601728) and CDKN2A (OMIM: 600160) were amplified, and Sanger sequencing was performed. The single-nucleotide polymorphism (SNP) information of detected genetic variations was retrieved from the NCBI dbSNP database (). Known disease-associated mutations in the PTEN and CDKN2A genes were retrieved from the Human Gene Mutation Database (HGMD; ).

Fluorescence in situ hybridization (FISH)

Dual-color interphase FISH was conducted on 4-μm paraffin-embedded tissue sections. The sections were pretreated with the FFPE FISH PreTreatment Kit (Abnova, Taipei, Taiwan). FISH was performed using the PTEN/CEN10p or CDKN2A/CEN9q probe (Abnova) according to the instructions of the manufacturer. Then, the sections were counterstained with 0.1 μM 4.6-diamidino-2-phenylindole (DAPI) for fluorescence microscopy observation, and a total of 200 interphase nuclei were analyzed for each specimen.

Follow-Up

Follow-up data for all patients were obtained during individual office visits or telephone interviews, with a mean follow-up time of 71.5±44.0 months (range 8–156 months). Radiological images, clinical symptoms, and survival status were evaluated. Progression-free survival (PFS) was defined as the interval between the first surgery to the date of tumor progression or tumor-related death. Overall survival (OS) was defined as the interval between the first surgery and the date of death from any cause.

Statistical Analysis

All statistical analyses were performed using SPSS software 25.0 (IBM, Armonk, NY, USA) and GraphPad Prism 8.0 (GraphPad, La Jolla, CA, USA). Continuous variables were expressed as the mean ± standard deviation (SD) when normally distributed or medians (interquartile ranges, IQR) when nonnormally distributed. Categorical values were presented as frequencies (percentage, %). Median PFS and OS were estimated using the Kaplan–Meier method and compared using the Log rank test. Normal distribution was tested using a D’Agostino and Pearson omnibus normality test. Chi-square tests and independent-sample Student’s t tests were performed to determine the correlation between molecular characteristics and clinicopathological features. All analyses were two-sided, and the statistically significant level was set at 0.05.

Results

Clinicopathological Characteristics

The demographic and clinical characteristics of patients with chordoma are summarized in Table 1. Histopathologically, conventional chordomas were composed of physaliferous cells arranged in nests or clusters embedded in a myxoid matrix, and the tumor cells were separated by fibrous septa into lobules (Figure 1A); chondroid chordomas exhibited islands of hyaline-type chondroid or cartilaginous matrix in addition to conventional chordoma appearance (Figure 1B). There was no significant difference in age (P=0.618), sex (P=0.835), maximal diameter of the tumor (P=0.273), surrounding tissue invasion (P=0.636), or recurrence rate (P=1.000) between patients with conventional chordoma and those with chondroid chordoma. Conventional chordoma was more common in the cervical spine, and chondroid chordoma was more common in the sacrococcygeal region (P=0.031). Kaplan–Meier survival curves revealed no significant difference in PFS or OS between patients with conventional chordoma and those with chondroid chordoma (median PFS: 120 months vs 109 months, P=0.605; median OS: 120 months vs 118 months, P=0.652).

Table 1

Clinical Characteristics of Patients with Chordomas

Characteristics	Conventional Chordoma (n=26)	Chondroid Chordoma (n=16)	Statistical Value (Chi-Square or t value)	P value
Gender
Male (n)	18 (69.2%)	13 (81.2%)	0.249a	0.618
Female (n)	8 (30.8%)	3 (18.8%)
Age (years; mean ± SD)	50.85 ± 13.94	49.94 ± 13.15	0.210b	0.835
Location
Cervical (n)	17 (65.4%)	5 (31.3%)	4.627c	0.031*
Sacrococcygeal (n)	9 (34.6%)	11 (68.7%)
Maximal diameter (mm; mean ± SD)	64.73 ± 23.97	75.88 ± 41.14	1.113b	0.273
Surrounding tissue invasion
Presence (n)	16 (61.5%)	11 (68.7%)	0.224c	0.636
Absence (n)	10 (38.5%)	5 (31.3%)
Recurrence
Presence (n)	13 (50.0%)	8 (50.0%)	0.000c	1.000
Absence (n)	13 (50.0%)	8 (50.0%)

Notes: aYates’ corrected Chi-square test; bIndependent-sample t test; cChi-square test. *P<0.05.

Figure 1

Histopathology of chordomas. (A) Conventional chordomas were composed of physaliferous cells that were separated by fibrous septa into lobules. (B) Chondroid chordomas exhibited islands of hyaline-type chondroid or cartilaginous matrix in addition to conventional chordoma appearance.

Genetic Alterations of PTEN and CDKN2A

Gene sequencing revealed an intron SNP (c.80–96A>G) and a missense mutation (c.10G>A; p.Gly4Arg) in the PTEN gene (Figure 2A and B). The intron SNP was identified in 15 (35.7%) patients; the corresponding allele frequency reported in the dbSNP database is ~38.9% (TOPMed Project; not adjusted for ethnicity). The missense mutation c.10G>A was identified in 18 (42.9%) patients, which was significantly higher than the reported allele frequency (~13.5%; TOPMed Project; not adjusted for ethnicity) in the dbSNP database. Additionally, a missense mutation (c.442G>A; p.Ala148Thr) was identified in the CDKN2A gene in three patients (7.1%) (Figure 2C); the corresponding allele frequency in the dbSNP database was only ~2.1% (TOPMed Project; not adjusted for ethnicity). The mutation c.442G>A in CDKN2A has been identified to be disease-associated in the HGMD database. The genetic landscape was presented in Figure 3.

Figure 2

Gene sequence chromatograms. Gene sequencing reveals a missense mutation (c.10G>A; p.Gly4Arg) in the PTEN gene (A), an intron SNP (c.80–96A>G) in the PTEN gene (B), and a missense mutation (c.442G>A; p.Ala148Thr) in the CDKN2A gene (C).

Figure 3

Genetic landscape and immunoreactivity. Genetic landscape map shows the distribution of genetic alterations of PTEN and CDKN2A in chordomas as well as the immunoreactivity for PTEN, CDKN2A and Ki-67.

Copy Number Loss of PTEN and CDKN2A and Clinical Relevance

FISH analysis identified heterozygous loss of the PTEN locus in 20 (47.6%) cases and homozygous deletion of PTEN in 5 (11.9%) cases (Figure 4A-C). Moreover, 22 (52.4%) cases displayed heterozygous loss of the CDKN2A locus, and 6 (14.3%) showed a homozygous deletion of CDKN2A (Figure 4D-F). Copy loss of PTEN and CDKN2A was highly correlated with the low expression level of corresponding proteins (P<0.001 for PTEN and P=0.010 for CDKN2A; Table 2). There was no significant difference in the frequency of PTEN or CDKN2A loss between conventional chordomas and chondroid chordomas (P=0.218 for PTEN and P=0.615 for CDKN2A; Table 3). Kaplan–Meier curves with the Log rank test showed that there was no significant difference in PFS or OS between patients with intact PTEN and patients with PTEN loss (P=0.570 for PFS, and P=0.330 for OS) or between patients with intact CDKN2A and patients with CDKN2A loss (P=0.318 for PFS, and P=0.114 for OS) (Figure 4G-J).

Figure 4

Fluorescence in situ hybridization and Kaplan–Meier survival analysis. Fluorescence in situ hybridization showed intact PTEN loci (A), heterozygous loss of the PTEN locus (B), homozygous loss of PTEN (C), intact CDKN2A loci (D), heterozygous loss of the CDKN2A locus (E), and homozygous deletion of CDKN2A (F); (A–C) PTEN and CEN10p were indicated in red and green, respectively; (D–F) CDKN2A and CEN9q were indicated in red and green, respectively. Kaplan–Meier curves showed that there was no significant difference in progression-free survival time (G; P=0.570) or overall survival time (H; P=0.330) between patients with intact PTEN and patients with PTEN loss, and there was no significant difference in progression-free survival time (I; P=0.318) or overall survival time (J; P=0.114) between patients with intact CDKN2A and patients with CDKN2A loss. The dashed lines indicate 95% confidence interval.

Table 2

Correlation Between Genetic Variations of PTEN and CDKN2A and Immunoreactive Scores in Chordomas

Genetic Variation	PTEN Immunoreactive Score (Mean ± SD)	t value	P value	CDKN2A Immunoreactive Score (mean ± SD)	t value	P value
PTEN intact	10.53 ± 2.13	6.328	<0.001*	–	–	–
PTEN deletion	4.92 ± 3.20			–
With c.80–96A>G in PTEN	7.40 ± 2.87	0.254	0.800	–	–	–
Without c.80–96A>G in PTEN	7.07 ± 4.47			–
With c.10G>A in PTEN	6.78 ± 2.60	0.583	0.563	–	–	–
Without c.10G>A in PTEN	7.50 ± 4.74			–
CDKN2A intact	–	–	–	9.29 ± 3.38	2.688	0.010*
CDKN2A deletion	–			5.96 ± 3.95
With c.442G>A in CDKN2A	–	–	–	5.33 ± 1.16	N.A.a	N.A.a
Without c.442G>A in CDKN2A	–			7.21 ± 4.17

Notes: aThe statistical comparison was not performed due to a low power, as the missense mutation (c.442G>A) was only identified in 3 patients. *P<0.05.

Table 3

Copy Number Loss of PTEN and CDKN2A and Immunoreactive Scores in Chordomas

Copy Number Variation	Conventional Chordoma (n=26)	Chondroid Chordoma (n=16)	Statistical Value (Chi-Square or t value)	P value
PTEN
Intact (n)	8 (30.8%)	9 (56.3%)	3.051a	0.218
Heterozygous loss (n)	15 (57.7%)	5 (31.2%)
Homozygous loss (n)	3 (11.5%)	2 (12.5%)
CDKN2A
Intact (n)	10 (38.5%)	4 (25.0%)	0.973a	0.615
Heterozygous loss (n)	13 (50.0%)	9 (56.3%)
Homozygous loss (n)	3 (11.5%%)	3 (18.7%)
PTEN immunoreactive score (mean ± SD)	6.85 ± 3.82	7.75 ± 4.19	0.718b	0.477
CDKN2A immunoreactive score (mean ± SD)	7.46 ± 3.73	6.44 ± 4.58	0.792b	0.433

Notes: aChi-square test; bIndependent-sample t test.

Expression of PTEN and CDKN2A and Clinical Relevance

According to the semiquantitative immunoreactive scores, the patients were divided into two groups based on the average PTEN expression cutoff value of 7.19. The group with lower PTEN expression (immunoreactive score ≤ 6) contained 20 patients (Figure 5A), and the group with higher expression (immunoreactive score ≥ 8) contained 22 patients (Figure 5B). The detailed clinicopathological data are summarized in Table 4. We compared the PFS and OS between these two groups using Kaplan–Meier curves with the Log rank test. The group with lower PTEN expression showed significantly shorter PFS and OS than those of the group with higher expression (median PFS: 75 months vs 134 months, P=0.023; median OS: 109 months vs 156 months, P=0.018) (Figure 5C and D). A low expression level of PTEN was an independent risk factor for PFS (hazard ratio [HR], 2.722; 95% confidence interval [CI], 1.156–6.407; P=0.023) and OS (HR, 2.966; 95% CI, 1.140–7.715; P=0.018).

Figure 5

Immunohistochemical staining and Kaplan–Meier survival analysis. According to the expression level of PTEN, the patients were categorized into the low expression group (A) and the high expression group (B). Kaplan–Meier curves showed that the patients with low expression of PTEN had significantly shorter progression-free survival time (C; P=0.023) and overall survival time (D; P=0.018) than the patients with high expression of PTEN. Based on the expression level of CDKN2A, the patients were classified into the low CDKN2A expression group (E) and the high CDKN2A expression group (F). Kaplan–Meier curves revealed that there was no significant difference in progression-free survival time (G; P=0.418) or overall survival time (H; P=0.250) between these two groups. The dashed lines indicate 95% confidence interval.

Table 4

Clinicopathological Characteristics of Patients with Different PTEN Expressions

Characteristics	Low Expression of PTEN (n=20)	High Expression of PTEN (n=22)	Statistical Value (Chi-Square Value or t value)	P value
Gender
Male (n)	18 (90.0%)	13 (59.1%)	5.178a	0.023*
Female (n)	2 (10.0%)	9 (40.9%)
Age (years; mean ± SD)	50.80 ± 13.05	50.23 ± 14.18	0.136b	0.893
Histopathological subtype
Conventional (n)	14 (70.0%)	12 (54.5%)	1.061a	0.303
Chondroid (n)	6 (30.0%)	10 (45.5%)
Location
Cervical (n)	11 (55.0%)	11 (50.0%)	0.105a	0.746
Sacrococcygeal (n)	9 (45.0%)	11 (50.0%)
Maximal diameter (mm; mean ± SD)	73.35 ± 37.42	65.00 ± 25.49	0.852b	0.399
Ki-67 labeling index
< 5%	6	15	10.786a	0.005*
5–10%	7	7
>10%	7	0
Surrounding tissue invasion
Presence (n)	14 (70.0%)	13 (59.1%)	0.543a	0.461
Absence (n)	6 (30.0%)	9 (40.9%)
Recurrence
Presence (n)	15 (75.0%)	6 (27.3%)	9.545a	0.002*
Absence (n)	5 (25.0%)	16 (72.7%)

Notes: aChi-square test; bIndependent-sample t test. *P<0.05.

Then, based on the average CDKN2A expression cutoff value of 7.07, the patients were classified into the low CDKN2A expression group (immunoreactive score ≤ 6; n=19) and the high CDKN2A expression group (immunoreactive score ≥ 8; n=23) (Figure 5E and F). The detailed clinicopathological data are summarized in Table 5. Survival analysis revealed that there was no significant difference in PFS or OS between these two groups (median PFS: 109 months vs 120 months, P=0.418; median OS: 109 months vs 156 months, P=0.250) (Figure 5G and H). Low expression levels of CDKN2A could not serve as an independent risk factor for PFS (HR, 1.423; 95% CI, 0.598–3.390; P=0.418) or OS (HR, 1.758; 95% CI, 0.672–4.601; P=0.250). The detailed statistical data are summarized in Table 6.

Table 5

Clinicopathological Characteristics of Patients with Different CDKN2A Expressions

Characteristics	Low Expression of CDKN2A (n=19)	High Expression of CDKN2A (n=23)	Statistical Value (Chi-Square Value or t value)	P value
Gender
Male (n)	15 (78.9%)	16 (69.6%)	5.178a	0.023*
Female (n)	4 (21.1%)	7 (30.4%)
Age (years; mean ± SD)	51.95 ± 13.52	49.30 ± 13.64	0.627b	0.534
Histopathological subtype
Conventional (n)	10 (52.6%)	16 (69.6%)	1.265a	0.261
Chondroid (n)	9 (47.4%)	7 (30.4%)
Location
Cervical (n)	8 (42.1%)	14 (60.9%)	1.469a	0.226
Sacrococcygeal (n)	11 (57.9%)	9 (39.1%)
Maximal diameter (mm; mean ± SD)	73.84 ± 39.42	64.96 ± 23.54	0.905b	0.371
Ki-67 labeling index
< 5%	5	16	9.323a	0.010*
5–10%	8	6
>10%	6	1
Surrounding tissue invasion
Presence (n)	11 (57.9%)	16 (69.6%)	0.617a	0.432
Absence (n)	8 (42.1%)	7 (30.4%)
Recurrence
Presence (n)	14 (73.7%)	7 (30.4%)	7.785a	0.005*
Absence (n)	5 (26.3%)	16 (69.6%)

Notes: aChi-square test; bIndependent-sample t test. *P<0.05.

Table 6

Kaplan–Meier Analysis of Progression-Free Survival and Overall Survival in Patients with Chordoma

Biomarker	Progression-Free Survival			Overall Survival
	Hazard Ratio	95% Confidence Interval	P value	Hazard Ratio	95% Confidence Interval	P value
Copy number loss
PTEN gene (loss vs intact)	1.290	0.537 ~ 3.097	0.570	1.568	0.604 ~ 4.070	0.330
CDKN2A gene (loss vs intact)	0.632	0.218 ~ 1.835	0.318	0.457	0.125 ~ 1.672	0.114
Protein expression
PTEN expression (low vs high)	2.722	1.156 ~ 6.407	0.023*	2.966	1.140 ~ 7.715	0.018*
CDKN2A expression (low vs high)	1.423	0.598 ~ 3.390	0.418	1.758	0.672 ~ 4.601	0.250

Note: *P<0.05.

Discussion

We performed a systematic review on genetic alterations at the PTEN and CDKN2A loci in chordomas. PubMed, Embase, Web of Science, and Google Scholar databases were used for the literature search, and a total of 19 relevant articles were retrieved. The detailed findings are summarized in Table 7. Karyotypic abnormalities involving chromosomes 9 and 10 in chordoma were first reported by Mertens et al in 1994 in which del(9) and –10 were identified in one case with sacral chordoma.11 In 2001, Scheil et al analyzed chromosomal imbalances in a series of 16 chordomas using comparative genomic hybridization in which loss of 9p and 10 chromosome segments was identified in 31% and 19% of chordoma samples, respectively. Additionally, they also established the first chordoma cell line U-CH1 from a recurrent sacrococcygeal chordoma, and U-CH1 exhibited del(9) and del(10) karyotypic abnormalities.4 In 2014, Kuźniacka et al used the conventional G-banding method to analyze chromosomal abnormalities in seven chordoma samples, and they identified loss of 9p and 10 chromosome segments in two and three cases, respectively.5 However, due to the technical limitations of conventional chromosome banding analysis, the locus of chromosomal deletion/loss was not accurately mapped in these studies. In 2008, Hallor et al first used array-based comparative genomic hybridization (aCGH) to analyze the copy number changes in chordomas; these scholars further mapped the chromosomal deletion to the 9p21 segment, and they proposed that this region covered the CDKN2A and CDKN2B loci.12 In the following years, chromosome segments 9p and 10 have become the hotspot regions for studying genetic alterations in chordomas. Many scholars analyzed these chromosomal abnormalities in chordoma tissues and cell lines using various technical platforms, such as aCGH,6,7,13 FISH,13–15 PCR-based microsatellite loss of heterozygosity (LOH) analysis,16 GeneChip array,9 or genome-wide high-resolution SNP-array.17,18 Accumulating evidence has revealed a frequency of PTEN deletion in chordomas ranging from 19% to 80%, and the frequency of CDKN2A deletion in chordomas has been reported to range from 12% to 80%.4–6,12-14,16–18 In our study, loss of PTEN and CDKN2A was identified in 59.5% and 66.7% of all chordoma cases, respectively.

Table 7

Systematic Review of Genetic Alterations at the PTEN and CDKN2A Loci in Chordomas

No.	Author (Year)	Materials/Objects	Method	Findings
1	Mertens et al (1994)11	Human sacral chordoma samples (n=3)	Chromosome banding analysis	Karyotype included del(9) and –10 in one case.
2	Scheil et al (2001)4	Human chordoma samples (n=16); U-CH1 cell line	Comparative genomic hybridization; chromosome banding analysis; FISH	Loss of 9p was identified in 5 of 16 (31%) cases, and loss of chromosome 10 segment was noted in 3 of 16 (19%) cases; U-CH1 had del(9) and del(10) chromosomal abnormalities.
3	Kuźniacka et al (2004)5	Human chordoma samples (n=7)	Chromosome banding analysis	Loss of chromosome 10 segment was noted in 3 of 7 (43%) cases, and loss of 9p sequences was identified in 2 of 7 (29%) cases.
4	Hallor et al (2008)	Human chordoma samples (n=26)	Array comparative genomic hybridization; FISH	Heterozygous deletion at 9p21 covering the CDKN2A locus was observed in 15 of 26 (58%) tumor samples, homozygous deletion was found in 3 of 26 (12%) samples.
5	Han et al (2009)34	Human chordoma samples (n=11), and U-CH1 cell line	Immunohistochemistry; western-blotting	Weak PTEN staining was demonstrated in 4 of 11 (36%) cases, and 6 of 11 (55%) cases were negative for PTEN staining, whereas 1 (9%) case showed positive PTEN staining; PTEN expression was not observed in U-CH1 cells.
6	Presneau et al (2009)35	Human chordoma samples (n=50)	Immunohistochemistry	PTEN displayed no immunoreactivity in 7 of 43 (16%) chordoma samples, and CDKN2A showed no immunoreactivity in 46 of 48 (96%) chordoma samples.
7	Sommer et al (2010)36	Human chordoma samples (n=27)	Immunohistochemistry	CDKN2A expression was absent in 20 of 27 (74%) chordomas. There was no significant difference in median disease-free survival between CDKN2A-positive chordomas (90 months) and negative chordomas (24 months).
8	Horbinski et al (2010)13	Human clival chordoma samples (n=28)	FISH; PCR-based microsatellite LOH analysis	9p LOH was identified in 3 of 25 (12%) samples, and 10q LOH was found in 12 of 21 (57%) samples; FISH identified homozygous loss of the 9p21 locus in 5 of 23 (22%) cases.
9	Brüderlein et al (2010)7	U-CH1 and U-CH2 chordoma cell lines	Array comparative genomic hybridization	Both U-CH1 and U-CH2 showed biallelic loss of the CDKN2A locus on chromosome 9p21; loss of the PTEN locus was noted in U-CH1, but not in U-CH2.
10	Shalaby et al (2010)37	Human chordoma samples (n=147)	Immunohistochemistry	PTEN expression was absent in 19 of 147 (13%) analyzed chordomas.
11	Le et al (2011)6	Human chordoma samples (n=20)	Array comparative genomic hybridization	Loss of CDKN2A was observed in 16/20 (80%) cases, six of which showed homozygous deletion; loss of PTEN was found in 16/20 (80%) cases; no mutations were found in any of the 21 chordoma samples.
12	Dewaele et al (2011)14	Human chordoma samples (n=42)	FISH	Loss of PTEN was observed in 7 of 18 (39%) analyzed samples.
13	Rinner et al (2012)9	MUG-Chor1 chordoma cell line	Affymetrix GeneChip Human Mapping SNP 6.0 array	The cell line displayed chromosomal loss at 10p15.3-q23.32 (includes the PTEN locus), and a homozygous loss was seen at 9p24.3-p13.1 (includes the CDKN2A locus).
14	Diaz et al (2012)17	Human chordoma samples (n=22)	Genome-wide high-resolution SNP-array	Loss of 9p involving CDKN2A was identified in 4 of 18 (22%) chordoma samples.
15	Choy et al (2014)8	Human chordoma samples (n=45)	Matrix-assisted laser desorption/ionization time-of-flight mass spectrometer	A nonsense mutation (p.R130*) in PTEN was identified in 1 of 45 (2.2%) chordoma samples; a nonsense mutation (p.R58*) was identified in CDKN2A in 2 of 45 (4.4%) chordoma samples.
16	Chen et al (2014)38	Human sacral chordoma samples (n=40)	Immunohistochemistry	Positive expression of PTEN was observed in 10/40 (25%) chordoma samples; patients with negative PTEN expression had a poorer prognosis than those with positive PTEN expression (P=0.011).
17	Lee et al (2015)16	Human chordoma samples (n=23)	PCR-based microsatellite LOH analysis	LOH was observed at 10q23 (includes the PTEN locus) in 6 of 23 (26%) specimens; chordoma tissues with LOH at 10q23 showed decreased or absent PTEN expression, and those without LOH at 10q23 demonstrated positive PTEN expression.
18	Wang et al (2016)18	Human chordoma samples (n=24)	Genome-wide high-resolution SNP-array	Copy number loss of CDKN2A was identified in 12 of 22 (55%) chordoma samples, and homozygous deletion of CDKN2A was detected in 5 of 22 (23%) samples.
19	Cottone et al (2020)15	Human chordoma samples (n=384)	Immunohistochemistry; FISH	Complete loss of chromosome 9 was observed in 3 cases; Positive CDKN2A protein expression was noted in 141 of 303 (47%) samples; copy number loss of CDKN2A gene was found in 138 of 274 (50%) samples, including heterozygous loss in 44 samples and homozygous loss in 48 samples.

Abbreviations: FISH, fluorescent in situ hybridization; LOH, loss of heterozygosity.

Point mutations in PTEN and CDKN2A have rarely been studied in the existing literature. In the retrieved articles, only one mutation (p.R130*) in PTEN and one mutation (p.R58*) in CDKN2A were identified by mass spectrometry in chordomas.8 In the current study, three genetic variants were detected, including a missense mutation (c.10G>A; p.Gly4Arg) and an intron mutation (c.80–96A>G) in the PTEN gene and a missense mutation (c.442G>A; p.Ala148Thr) in the CDKN2A gene. There has been no functional validation for these variants. Slattery et al reported that the intron SNP (c.80–96A>G; rs1903858) in PTEN may be associated with breast cancer risk (odds ratio [OR], 0.88; 95% CI, 0.80–0.97; P=0.032),19 while other studies found that this SNP was not associated with primary ovarian insufficiency or prostate cancer.20,21 The missense mutation c.10G>A (rs12573787) in PTEN was found to be associated with an increased risk of chronic myeloid leukemia (OR, 1.71; 95% CI, 1.11–2.63; P=0.016),22 but not with endometrial cancer risk.23 According to the literature, the missense mutation c.442G>A (rs3731249) in CDKN2A was considered to have pathogenicity, and it has been documented as a disease-associated mutation in the HGMD database. A large number of studies have suggested that c.442G>A may contribute to acute lymphoblastic leukemia,24–27 breast cancer,28,29 melanoma,30,31 ovarian cancer,32 and colorectal cancer.33 We found that the frequency of this missense mutation in chordomas (7.1%) was markedly higher than that in the general population (2.2%), which was consistent with the data documented in the above studies. Due to the limited sample size, an accurate OR value could not be calculated in this study.

In addition to chromosomal aberrations in the PTEN and CDKN2A loci, some studies evaluated the expression levels of the PTEN and CDKN2A proteins using immunohistochemistry or Western blotting.15,34-38 Based on the reported results in these previous studies, PTEN staining was negative in 13%~55% of chordomas,34,35,37,38 and immunoreactivity for CDKN2A was absent in 53%~96% of chordomas.15,35,36 However, only two studies analyzed the correlation between the protein expression levels and clinical prognosis.36,38 Sommer et al found no significant difference in median disease-free survival between CDKN2A-positive chordomas (90 months) and CDKN2A-negative chordomas (24 months).36 Chen et al observed that patients with negative PTEN expression had a poorer prognosis than those with positive PTEN expression.38 Notably, a recent study proposed that the absence of CDKN2A protein expression could not be explained by copy number loss of the CDKN2A gene.15 In particular, there has been a paucity of studies concerning genetic alterations of PTEN and CDKN2A in distinct histological subtypes of chordoma. Therefore, we designed this study to investigate the molecular characteristics and expression levels of PTEN and CDKN2A in conventional and chondroid chordomas, and the clinical relevance was analyzed. We found that the copy loss of PTEN and CDKN2A was closely correlated with the low expression level of the corresponding proteins. Additionally, our results indicated that loss of PTEN expression was associated with a poorer prognosis, but loss of CDKN2A expression was not correlated with clinical prognosis, which was highly consistent with the above two studies.36,38

The tumor-suppressor PTEN was first identified in 1997,39 and it has been found to be a negative regulator of the phosphatidylinositol-3-kinase (PI3K)/AKT/mTOR signaling cascade.40 PTEN participates in the regulation of cell division, proliferation, apoptosis, migration, adhesion to surrounding tissues, and neovascularization.41 Increasing evidence indicates that PTEN expression is associated with the prognosis and drug response of several tumors, such as prostate cancer, breast cancer, and endometrial carcinoma.42 Our results suggest that PTEN expression may be a prognostic and predictive biomarker for chordomas; nevertheless, we also noted that there was no significant correlation between copy number loss of PTEN and clinical prognosis. The expression of PTEN can be controlled at the genetic, transcriptional, and posttranscriptional levels,43 which may explain the inconsistent findings regarding the prognostic value of PTEN loss and that of protein expression, while the definitive mechanism still needs further research. CDKN2A, the gene encoding the cell-cycle inhibitor p16CDKN2A (also known as p16, p16,INK4 p16,INK4A or CDK4I), was first identified in 1994.44 The CDKN2A/p16 protein plays an executional role in the cell cycle and senescence by regulating cyclin-dependent kinase (CDK) 4/6 and cyclin D complexes. Previous studies have shown that genetic and epigenetic aberrations of CDKN2A can promote tumorigenesis, recurrence, and metastasis in various cancers, such as lymphoma, melanoma, lung cancer, and gastric cancer.45 In the present study, we failed to find a causal link between CDKN2A deficiency and chordoma prognosis.

There are some limitations to this study. First, although several genetic variants were identified, no functional validation was performed. Second, as chordomas are solitary tumors, no peritumoral parenchyma specimens were available, and peripheral blood samples were not preserved during the hospitalization; we did not classify the genetic variations as somatic or germline. Third, the sample size is still limited. As there were no significant differences in survival time between conventional chordomas and chondroid chordomas, we combined all the patients for further survival analyses. In future studies, we will construct chordoma cell models with site-directed mutagenesis and analyze the potential pathogenic functions of genetic alterations.

Conclusion

This study delineated the genetic characteristics and expression of PTEN and CDKN2A in chordomas. Copy number loss of PTEN and CDKN2A is common and evenly distributed in conventional and chondroid chordomas. PTEN expression may serve as a prognostic and predictive biomarker for chordomas, while no correlation between CDKN2A expression and chordoma prognosis was noted.