Abstract

Giant Cell Tumour of Bone (GCTB) is a rare bone tumour. Locally aggressive and recurrent, it might evolve into pulmonary metastases. Our present work is aimed at investigating the involvement of the epidermal growth factor receptor (ErbB) family and its downstream effectors in the development and recurrence of GCTB. For this purpose, we used a cohort of 32 GCTB patients and we evaluated the clinicohistological features and the expression of RANKL, EGFR, and HER2. The mutation status of KRAS, PI3KCA, and PTEN gene as potential oncogene involved in GCTB was also evaluated. We found a significant correlation between advanced histological stages, overexpression of EGFR/HER2, and tumour recurrence. Moreover, two mutations were found in the PIK3CA gene: a missense mutation, 1634A>C, detected for the first time in GCTB patients, without influencing the stability of the protein, and a frameshift mutation, c.1658_1659delGTinsC, causing the loss of the protein kinase domain. Altogether, these results suggest that overexpression of HER2/EGFR, Campanacci, and histological stages could be used as a novel prognostic marker for GCTB recurrence.

1. Introduction

Giant Cell Tumour of Bone (GCTB) is an intermediate locally aggressive primary bone tumour. Due to its local recurrence and rare pulmonary metastases, GCTB is considered to have a low malignant potential [1, 2]. Aggressiveness characteristics include a high mitotic rate, necrosis, and recurrence after resection [3]. Local recurrence after surgery is estimated to occur in 15–50% of patients, with a higher prevalence after curettage than after en bloc resection [4]. Recurrence of GCTB is not fatal in most cases but can lead to disability and to a poor quality of life as a result of repeated and radical operations, loss of bone stock, and secondary arthritis of the joints [5]. In United States, GCTB accounts for approximately 3 to 5% of all primary bone tumours and 15 to 20% of all benign bone tumours [6, 7]. A slightly higher incidence was suggested in a population-based series from the Swedish Cancer Registry; of the 4625 bone tumours diagnosed over a 53-year period, 505 (11%) were GCTB [8]. Asian populations have a significantly higher incidence than Western populations. In China, GCTB represents approximately 20% of all primary bone tumours [9, 10].

GCTB usually occurs in patients between 20 and 40 years old with a slightly higher rate in women than men (1.5 to 1 ratio) [11]. GCTB is formed by multinuclear osteoclast-like cells. These cells are originated essentially from macrophages and monocytes stimulated by various cytokines [12]. GCTB is also formed by stromal cells which have been identified as the neoplastic cell population and evolved from mesenchymal stem cells [13, 14]. Identifying risk factors influencing recurrence in GCTB is still nonwell developed. Previous studies reported that recurrence is influenced by clinical features, such as gender, age, location, tumour size, and surgery method [15]. Moreover, certain factors were proposed as predictors of recurrence and aggressiveness markers in GCTB. Indeed, the overexpression of interleukin-17A (IL-17A) and β-catenin is closely associated with GCTB progression and recurrence [16]. In addition, it has been reported that Ezrin, VEGF, and CD44 play a role in invasion, metastasis, and recurrence in GCTB [17]. A recent study has indicated that postsurgery fluids are an adjuvant in the mechanism of tumour recurrence, increasing stem cell growth and AKT/mTOR activity [18].

A strong correlation exists between solid cancers and the overexpression of the ErbB receptors [19, 20]. However, this correlation was not studied in GCTB. The ErbB family members, especially EGFR and HER2, activate the proliferative signalling through the MAPkinase and the PI3kinase pathways [21]. ErbB downstream signalling such as PI3K/AKT and mitogen-activated protein kinase (MAPK) pathways was involved in bone loss. In fact, during osteolysis, the activation of nuclear factor-κB (NF-κB), MAPK, and PI3K/AKT signalling pathways is triggered by RANKL, inflammatory factors, and oxidative stress [22]. It has been demonstrated that PI3K and its downstream target AKT (serine/threonine kinase) play a significant role in the differentiation of progenitor cells into mature osteoblasts [23]. Moreover, the PI3K/AKT and MAPK pathways ensure specific roles in osteoblast differentiation and in RANKL-induced osteoclastogenesis [2426].

PI3K/AKT and RAS/RAF/MEK/ERK pathways regulate certain cellular functions, including cell proliferation, differentiation, survival, and apoptosis [27]. PTEN (phosphatase and TENsin homolog), a tumour suppressor gene, exerts a negative control of the phosphoinositide 3-kinase signalling pathway and cell proliferation/survival by dephosphorylating the phosphatidylinositol 3,4,5-triphosphate [2830]. Deregulation of these pathways leads to the occurrence and progression of certain cancers mainly by somatic mutations [31, 32]. The most common PIK3CA mutations had been reported in exons 9 and 20, representing different hotspot mutations [33]. There are 3 hotspot mutations in PIK3CA gene: E542K, E545K at exon 9 (helix domain), and H1047R at exon 20 (kinase domain) [34]. The 3 hotspots represent almost 80% of PIK3CA mutations and lead to constitutive PI3K activity by different mechanisms [35]. KRAS is an oncogene playing a significant role in downstream signalling of the EGFR pathway. KRAS mutations are observed in many human malignancies, such as non-small-cell lung cancer, colorectal cancer, and pancreatic ductal adenocarcinoma [3639]. KRAS mutations commonly arise in codon 12 and 13 of exon 2 [40]. Considering all this background, a possible overexpression of RANK/RANKL and members of the ErbB receptor family in GCTB, in addition to the involvement of PIK3CA, PTEN, and KRAS mutations, is suggested to influence GCTB pathogenesis.

For this purpose, we collected a cohort of GCTB patients within the region of Sfax, South Tunisia, and we performed several investigations, including clinical and histological characterisation, and mutational analysis. The main goal was to propose a novel prognostic marker for GCTB recurrence.

2. Materials and Methods

2.1. Clinical Cases and Samples

32 patients with GCTB (11 men, 21 women) were recruited in the Orthopaedics and Traumatology Department, Habib Bourguiba University Hospital, Sfax. From 2015 to 2017, complete clinical and pathological data were generated for each patient. Biopsies were collected after surgery and prepared for further analysis. Informed consent was obtained for all individual participants included in the study.

2.2. Tumour Morphology

Formalin-Fixed Paraffin-Embedded (FFPE) tissue sections (4–5 mm) were deparaffinised, hydrated, and stained with hematoxylin and eosin (H&E). Baseline and on-study tissue samples from each patient were evaluated for changes in overall tissue composition and architecture by light microscopy. Then, the collected tumours were graded according to the “Campanacci” and “Jaffe-Lichtenstein” grading.

2.3. Campanacci Grading

GCTBs were classified by Campanacci et al. based on radiographic appearance: stage I: latent, stage II: active, and stage III: aggressive [41].

2.4. Jaffe-Lichtenstein Grading

The histological grading system of Jaffe et al. can be summarized as follows: For grade 1, tumours are the least aggressive and show a moderate loss, with vascular stroma composed of spindle and ovoid cells, with few or no mitotic figures and numerous giant cells. For grade 2, tumours manifest, cytologically, a very compact cellular stroma showing definite evidence of atypia, tend strongly toward recurrence, and in some cases eventually undergo malignant transformation. For grade 3, tumours show a compact stroma in which the stromal cells are universally pleomorphic and the giant cells are smaller and less numerous [42].

2.5. Immunohistochemistry

For each sample, paraffin-embedded sections were cut 3 μm thin, dewaxed in xylene, rehydrated, and submitted to antigen retrieval in EDTA buffer (pH 9.0) or citrate (pH 6.0). Each slide was blocked for endogenous peroxidase activity through incubation in H2O2 3% for 5 min, followed by 60 min incubation at ambient temperature with anti-Her2 (mouse monoclonal antibody, clone CB11, dilution 1 : 400, Leica Biosystems) and anti-EGFR (mouse monoclonal antibody, clone EGFR.113, dilution 1 : 20, Leica Biosystems). For RANKL detection, an overnight incubation at 4°C with a 5 μg/ml working solution of an anti-RANKL (mouse monoclonal antibody IgG2B, clone: 70525, MAB626, antibody, R&D Systems) was adopted. Post Primary (rabbit anti-mouse IgG) was then used to detect mouse antibodies. Sections were incubated with the Novolink™ Polymer that recognizes rabbit immunoglobulins (Novolink™ Polymer Detection System, RE7290-K) and further incubated with the substrate/chromogen, 3,3-diaminobenzidine (DAB), prepared from DAB chromogenanaNovolink™ DAB substrate Buffer (polymer). Reactions with the peroxidase produced a visible brown precipitate at the antigen site. Sections were counterstained with hematoxylin and coverslipped. Immunohistochemical staining was evaluated according to intensity; intensities were scored as 0, no staining; +, weak staining; ++, moderate staining; or +++, strong staining. Positive and negative controls were used to validate antibodies.

2.6. Mutational Analysis
2.6.1. DNA Samples

Fresh-frozen surgically resected samples were collected from GCTB patients. Genomic DNA was extracted from tumour tissues using phenol-chloroform standard procedures.

(1) Sample Preparation and Phenol/Chloroform Extraction of DNA. Frozen GCTB tissues were cut up to 3 μm into small chips (2–4) using a cryostat and placed in a sterile tube. Tissues are disaggregated and treated with PK-buffer (2 M Tris-HCl, 0.5 M EDTA, 1.5 M NaCl) without Sodium Dodecyl Sulphate (SDS). The tissue was homogenized while it is submerged in PK-buffer using a syringe until complete homogenization. 20 μl SDS was added then vortexed and resuspended in 50 μl proteinase K (10 mg/ml stock; HiMedia). An equal volume of phenol : chloroform : isoamyl alcohol was added and mixed by inverting the tube for 20 min and centrifuging at 2500 tr/min (at 4°C) for 10 min. Chloroform is then used to facilitate the removal of phenol. The sample was incubated overnight at -20°C. DNA is subsequently purified by precipitation in ethanol. Finally, DNA is resolubilized in Tris-EDTA buffer.

2.6.2. Sequencing of the PIK3CA, PTEN, and KRAS Genes

The exons 9 and 20 of PIK3CA, exons 5 and 7 of PTEN, and exon 2 of KRAS genes were amplified by PCR, using primers and PCR condition described previously. After PCR amplification, each PCR product was purified and subsequently analysed by direct sequencing in an ABI PRISM 3100-avant automated DNA sequencer, by means of the Big Dye Terminator Cycle Sequencing reaction kit v1.1. The resulted sequences were compared to human genomic DNA sequence using blast, and the mutations were identified using ENSEMBL and Cosmic database.

2.7. Computational Analysis
2.7.1. The Sequence Alignment

A BLAST homology search was performed using the program BLAST2SEQ available at the NCBI (National Center for Biotechnology Information) website (https://www.ncbi.nlm.nih.gov/blast/bl2seq) to compare individual nucleotide sequences with the human gene wild-type sequence.

The sequence alignment was performed using the ClustalW program (http://bioinfo.hku.hk/services/analyseq/cgi-bin/clustalw_in.pl).

2.7.2. Protein Structure Prediction

Structures are visualized, and figures were generated using the Biovia Studio visualizer (Accelrys Software Inc., Discovery Studio Modeling Environment, Release 2017 R2, San Diego: Accelrys Software Inc., 2007). The structure of the complex between P110 and nSH2-iSH2 domains of regulatory subunit taken as wild-type protein has the pdb code 4L23 (https://www.rcsb.org/). The mutant E545A was generated using the Biovia Studio visualizer by replacing E545 in the wild-type structure by alanine. The same software was used to optimize A545 configuration and to check the absence of any clashes with the surrounding residues.

2.8. Statistical Analysis

The SPSS statistical package (version 20.0; SPSS, USA) was employed for statistical analysis. Pearson’s chi-squared or Fisher tests were used to evaluate the association between Her2, EGFR, and RANKL biomarker overexpression and GCTB recurrence. A value of <0.05 was considered significant.

3. Results

3.1. Case Selection

A total of 32 patients, 11 men and 21 women, were recruited within our Orthopaedics and Traumatology Department. Clinical and pathological examination confirmed the GCTB diagnosis. Patient age ranged between 12 and 74 years, with a majority between 20 and 40 years (56.2%). Tumours were mainly located in the radius 10 (31%), tibia 9 (28%), femur 7 (22%), fibula 3 (9.5%), and other rare locations 3 (9.5%). We applied 2 grading systems in order to classify the tumours. According to Campanacci’s grading system, 6 tumours were stage 1, 19 stage 2, and 7 cases stage 3, respectively (Figure 1(a) and Table 1). According to the Jaffe-Lichtenstein grading system, 3 tumours were classified grade 1, 15 grade 2, and 14 grade 3 (Figure 1(b) and Table 1). Time of recurrence was quite disparate and ranged from 6 to 48 months and rated in 47% of the studied cases. Moreover, curettage has a significantly higher recurrence rate than wide resection. Indeed, 21 of the 32 patients studied underwent curettage as a treatment. Also, we have shown that the type of surgery is significantly associated with recurrence; in fact, among the 21 patients treated with curettage, 13 have recurrence (), and 2 among 7 patients treated with wide resection have recurrence (). Interestingly, the majority of patients with high radiological and histological grades and 10 over 15 recurrent patients with high histological grade underwent curettage as treatment. Therefore, there was no significant correlation between age, gender, and recurrence. However, according to both systems of grading and surgery methods, higher grades were significantly associated with recurrence (Table 1).

3.2. Expression Level of HER2, EGFR, and RANKL and GCTB Tumour Recurrence Correlation

Using immunohistochemistry, we analysed the expression level of HER2, EGFR, and RANKL in GCTB tumours. We found that HER2 was highly expressed, mainly in the cytoplasm of multinucleated giant cells corresponding to the nonneoplastic compartment (Figure 2(a)). Indeed, HER2 was expressed moderately (++) in 6 patients (18.75%) and strongly (+++) in 26 patients (81.25%). GCTB recurrence was found in 4 over 6 patients with moderate expression and 14 over 26 in patients with strong expression. However, there was no significant correlation between HER2 expression and recurrence (Table 2). We analysed then the expression level of EGFR, a strong expression level as confirmed by the high signal we found. EGFR is localized in both types of cells, within the neoplastic and the nonneoplastic compartments (Figure 2(b)). Distribution over tumour samples was as follows: weak expression (+) in 3 patients (9.375%), moderate expression (++) in 16 patients (50%), and strong expression (+++) in 13 patients (40.625%). Interestingly, we found that 70% of patients with strong staining are getting tumour recurrence, and thus, a significant positive correlation between EGFR overexpression and recurrence was confirmed. Moreover, 13 patients displayed an associated overexpression of EGFR and HER2 with a statistically significant tumour recurrence (Table 2). A contrario of HER2, RANKL was essentially observed in the neoplastic compartment of GCTB. We found that RANKL is weakly expressed in 9 patients (28.125%), moderately in 14 (53.125%), and strongly in 9 patients (18.75%) (Figure 2(c)). In addition, 77% of patients with strong expression were getting tumour recurrence (Table 2).

3.3. Mutational and Bioinformatics Analysis of GCTB Tumours

Based on key oncogenes described in Introduction, we decided to perform a genetic study focusing on PIK3CA, PTEN, and KRAS genes. We sequenced these genes in all our GCTB patients. Interestingly, we did not found any mutations in exons 5 and 7 of PTEN gene, in exon 2 of KRAS gene, and in exon 20 of PIK3CA gene. However, we found that the PIK3CA gene harbours two mutations in exon 9. First, a missense mutation, c.1634A>C encoding p.E545A, substitutes the highly conserved “Glutamic acid, E” by “alanine, A” in the position 545, p.E545A within the helical domain (Figures 3(a) and 3(b)). As the missense mutation could affect the PIK3CA protein structure, we decided to perform a bioinformatics analysis on both wild-type and mutant protein. For this purpose, we used the program Biovia Studio visualizer (Accelrys Software Inc., Discovery Studio Modeling Environment, Release 2017 R2, San Diego: Accelrys Software Inc., 2007). In the wild-type protein, E545 belonging to the PIK3CA (p110) is establishing ionic interaction with K379 of the nSH2 domain of the regulatory subunit (p85) (Figure 3(c), i and ii). E545 is also ion pairing with the close K548 of P110. The NH group of E545 is stabilized by hydrogen bonding with the carbonyl group of L380. A replacement of E545 by alanine disrupted the ion pairing with K379 without disrupting the complex. This would lead to lower affinity between the two proteins, but the interaction would remain due to the large patch of amino acids involved in the complex stabilization. Like the wild-type protein, E545A mutant is also stabilized by hydrogen bonding with L380. Furthermore, the nondisturbance of the interaction between the regulatory and the catalytic subunits of PI3K by this mutation is very likely due to the small side chain of alanine residue preventing clashes with surrounding residues. A contrario, Ala545 in the mutated protein is even stabilized through hydrophobic contacts with Leu380 and Ile381 belonging to the regulatory subunit (Figure 3(c), i and iii). Altogether, these observations explain a lack of effect of the mutation E545A on the PIK3CA regulation by the regulatory subunit.

The second mutation we found is a frameshift mutation, c.1658_1659delGTinsC encoding p.S553Tfs7. It is a deletion of the dinucleotide GT plus an insertion of a C resulting in a serine-to-threonine substitution at codon 553, followed by frameshift in a stop codon, 7 amino acids later. These changes induced the creation of a new premature stop codon (TAA) at position 559 and cause the loss of the PI3K protein catalytic domain (Figures 4(a)4(c)).

4. Discussion

In this study, we aimed to determine the risk factors for GCTB recurrence in a cohort of patients and to investigate the HER2, EGFR, and RANKL expression within the GCTB, in relationship with tumour recurrence. We aimed also to evaluate the involvement of PI3K/AKT and MAPK pathways in GCTB development.

Within our cohort, we found a feminine predominance that confirms previous reports [43]. This predominance is not correlated with GCTB tumour recurrence. However, a recent study suggested that recurrence should be more considered in men with GCTB than women [44]. In our study, no significant correlation was found between age and recurrence despite the study of Hu et al. that revealed a positive association between age and recurrence (2016) [45]. Similarly, our analysis revealed that location is not associated with recurrence while others reported that localization in the radius and proximal tibia was more frequently associated with tumour recurrence [46]. Grading systems Campanacci and Jaffe revealed a strong correlation between tumour grade and recurrence in our cohort and are in line with previous studies [47, 48].

We studied then the expression of HER2 and found that over 80% of the analysed tumours overexpressed the protein. HER2 expression was exclusively cytoplasmic, and no concomitant membranous immunostaining has been detected in giant cells. Multinucleated giant cells are “reactive” components and originate from osteoclast precursors (bone marrow-derived monocytes) based on a resemblance in differentiation markers. Despite their nonneoplastic origin, multinuclear giant cells caused lacunar bone resorption, the major characteristics of this tumour [49]. In addition in our case, these cells are overexpressing the most oncogenic receptor in solid cancer (EGFR/HER2) [50]. Despite the fact that literature may qualify the absence of HER2 membrane staining, as an artefact [51], in numerous studies, overexpression of HER2 within the cytoplasm of cancer cells was a fact [52, 53]. Another study suggested that cytoplasmic immunostaining of HER2 was related to the presence of a truncated form “p95her2 [54]. Furthermore, cytoplasmic HER2 expression was described in osteosarcoma and Ewing and Synovial sarcomas [55, 56], and in 2005, Scotlandi et al. reported an association of HER2 expression with a worse prognosis in osteosarcoma (2005) [57].

We studied also the expression of EGFR and found that GCTB tumours exacerbate its expression in both types of compartments, in multinuclear giant cells and mononuclear neoplastic cells. Infante et al. suggested that EGFR could play as a RANK downstream regulator and/or a class of RANK coreceptors in osteoclasts (2008) [58]. It is well known that RANKL stimulation triggers the recruitment of TRAFs to RANK and the subsequent activation of the NF-κB and MAPK cascades (ERK, JNK, and p38) with the induction of Src- and PI3K-dependent AKT activation [59]. EGFR activation could be an extra component of RANK downstream signalling, and activated EGFR could reinforce the assembly of multiprotein signalling complexes including RANK and Gab2 leading to a further activation of JNK, AKT/PKB, and NF-κB [58]. We found also a significant correlation between overexpression of HER2/EGFR with tumour recurrence. These results confirm what has been reported in terms of correlation of HER2/EGFR to highly aggressive tumours [50, 59, 60].

We showed also that RANKL is overexpressed in our cohort, confirming what has been published previously in different tumours, including breast cancer, colorectal cancer, osteosarcoma, and GCTB [61, 62]. We showed that RANKL is predominantly expressed in mononuclear stromal cells corresponding to the neoplastic compartment. The stromal cells regulate the formation of osteoclast-like giant cells in the neoplasm [63]. In fact, GCTB is composed of “neoplastic” cells with fibroblast-like spindled morphology in harmony with “reactive” macrophage-like round or osteoclast-like multinucleated giant cells [64, 65]. Previous studies confirmed that RANK/RANKL interaction and macrophage colony-stimulating factor (M-CSF) are key players in osteoclastogenesis by recruiting mononuclear osteoclast precursor cells that differentiate into multinucleated osteoclast-like giant cells leading to bone resorption. The reciprocal and orchestrated actions between mononuclear stromal cells and multinucleated giant cells help in the understanding of the molecular pathogenesis and tumour biology of GCTB [66]. Furthermore, it was previously shown that blocking RANKL binding to RANK arrests tumour progression, improves survival, and also inhibits pulmonary metastasis [6769].

PI3K/AKT and MAPK constitute important downstream signalling pathways of EGFR and HER2 [70]. We showed in our study that no mutations in PTEN and KRAS genes were detected in all GCTB patients. However, in PIK3CA gene, we detected two mutations in exon 9 among the helical domain: a frameshift mutation 1658_1659delGTinsC; p.S553Tfs7 and a missense mutation c.1634A>C encoding p.E545A, and in all patients, no mutations were detected in exon 20 (encoding for the kinase domain) of PIK3CA gene.

The frameshift mutation shortens the protein and causes the loss of the PI3K protein catalytic domain. As a consequence, this mutation will interfere with the conversion of phosphatidylinositol 4,5-bisphosphate (PIP2) to phosphatidylinositol-3,4,5-triphosphate (PIP3). This PIP3 acts as a membrane anchor for the downstream proteins like PDK1 and AKT, and it seems that this mutation interferes with the PI3K/AKT antiapoptotic pathway through the inhibition of the AKT [71]. This frameshift mutation was described in a patient with endometrial cancer; it was associated to long-term nonprogressive disease as reported by Olivier Tredan (2013). IHC analysis confirmed the absence of phospho-AKT in the corresponding patient. No mutations have been proven in KRAS, PTEN, or AKT1 genes (Table 3) [72]. We can conclude that this frameshift mutation does not have any pathogenic potential.

The missense mutation does not change the protein activity, and we could speculate that E545A mutation in PIK3CA gene is not oncogenic in GCTB despite great evidences previously reported in other cancer models. Indeed, previous studies were pretending that E545A mutation is described as an oncogenic mutation in different types of cancer, such as breast and ovarian (Table 3) [73, 74]. Furthermore, Gymnopoulos et al. showed that mutations in PIK3CA, including E545A, gained function-induced oncogenic transformation (Table 3) (2007) [75].

However, some studies revealed that the coexistence of both mutations of PIK3CA, E545A, and S553T fs7 has a pathogenic potential in the Cowden syndrome, which is known to be a cancer predisposition syndrome. In fact, the heterozygosis for these two mutations in exon 9 of the PIK3CA gene led to an increase in kinase activity [76]. This pathogenic potential induced by the coexistence of both mutations was described also in non-small-cell lung cancer (NSCLC) with a locally advanced stage and with stage IV (Table 3). Therefore, the poor prognosis and the resistance to treatment in these patients might be related to other genetic alterations and factors [77]. In our patients therefore, no pathogenicity was recorded despite the coexistence of the two mutations in PIK3CA gene.

5. Conclusions

The present study is an attempt to identify possible criteria for predicting recurrence of Giant Cell Tumour of Bone. In addition to generally known parameters, like localization of the lesion, clear destruction of cortex with the presence of extra compartmental lesion, and histological criteria for malignancy, we have found the following independent predictors, namely, the overexpression of HER2 and EGFR in the recurrent group. The correlation between HER2/EGFR overexpression with tumour recurrence suggested that HER2/EGFR signalling plays a crucial role in aggressiveness and recurrence in GCTB. However, despite the occurrence of mutations in the protooncogenic PIK3CA gene, PI3K and MAPK pathways may not play a crucial role in the pathogenesis of GCTB.

Data Availability

All data generated or analyzed during this study are included in this article. All data and materials are presented in Materials and Methods and Results as shown in figures and tables.

Conflicts of Interest

All authors declare that there is no conflict of interest.

Acknowledgments

This work was funded by the Tunisian Ministry of Higher Education and Scientific Research. The authors thank the patients and their families for their cooperation in the present study.