Abstract
Background/Aim: Better diagnostic and prognostic markers are required for a more accurate diagnosis and an earlier detection of glioma progression and for suggesting better treatment strategies. This retrospective study aimed to identify actionable gene variants to define potential markers of clinical significance. Materials and Methods: 56 glioblastomas (GBM) and 44 grade 2-3 astrocytomas were profiled with next generation sequencing (NGS) as part of routine diagnostic workup and bioinformatics analysis was used for the identification of variants. CD34 immunohistochemistry (IHC) was used to measure microvessel density (MVD) and Log-rank test to compare survival and progression in the presence or absence of these variants. Results: Bioinformatic analysis highlighted frequently occurring variants in genes involved in angiogenesis regulation (KDR, KIT, TP53 and PIK3CA), with the most common ones being KDR (rs1870377) and KIT (rs3822214). The KDR variant was associated with increased MVD and shorter survival in GBM. We did not observe any correlation between the KIT variant and MVD; however, there was an association with tumour grade. Conclusion: This study highlights the role of single-nucleotide variants (SNVs) that may be considered non-pathogenic and suggests the prognostic significance for survival of KIT rs3822214 and KDR rs1870377 and potential importance in planning new treatment strategies for gliomas.
- Glioma
- microvascular proliferation
- KIT
- KDR
- artificial intelligence
- polymorphism
- next-generation sequencing
Diffuse gliomas are the most frequent primary malignant tumours of the central nervous system (CNS) in adults, with an annual incidence of approximately 5/100,000 per year (1). The majority consists of astrocytic tumours, which include more differentiated and slower-growing diffuse astrocytomas, isocitrate dehydrogenase (IDH) mutant (mut) (WHO grade 2) and high-grade astrocytomas, IDH mut (WHO grade 3 and 4), IDH wildtype (WT) GBM (WHO grade 4), and other genetically defined but rare subtypes (2, 3). IDH WT GBM is the most aggressive and most common type of glioma in adults, corresponding to 50% of all primary CNS malignancies (2). Lower grade (grades 2 and 3) IDH mut astrocytomas follow a natural progression to grade 4 IDH mut astrocytomas, or secondary GBM (2). GBM is characterised by strong resistance to chemotherapy (4, 5) and radiotherapy (6), whilst the often-precarious location and infiltrating growth of the tumour results in largely ineffective surgical removal. All these factors contribute to the dismal prognosis of grade 4 gliomas with a 5-year survival of 5% (7). The current diagnosis of gliomas relies on magnetic resonance imaging, histological assessment of tumour tissue (8, 9) and genetic profiling of tumour cells (2, 10). The histological hallmarks of GBM in comparison to lower grade astrocytomas are: i) high proliferation activity and ii) microvascular proliferation with or without necrosis (2, 11). GBM is known as one of the most highly angiogenic tumours and the extensive vascular network is associated with an elevated presence of pro-angiogenic factors released by the tumour cells (12-14). The genetic profile of gliomas is assessed using a combination of techniques, such as IHC, sequencing and copy number variation (CNV) studies. An epigenetic marker, methylguanine-DNA-methyltransferase (MGMT) promoter methylation is also routinely assessed as a crucial predictive parameter for patients' management (15). NGS is increasingly used to assess the presence of mutations in relevant genes (16, 17) as it offers a rapid and comprehensive molecular analysis. Using NGS, the mutation status of only a subset of panel genes is routinely considered for diagnostic work-up. Hence, the diagnostic and prognostic significance of the mutation status in many of the other oncogenes in the sequencing panels is not fully characterised in the context of glioma.
In this study, we investigated the correlation between these non-routinely assessed oncogenes from a panel of 50 genes together with clinical, morphological and molecular features of patients with diffuse gliomas, so as to isolate gene variants, which might hold a diagnostic or prognostic significance as well as potential relevance for treatment. We analysed tumour tissues from 100 patients diagnosed with diffuse astrocytic glioma ranging from WHO grade 2 to grade 4 and identified frequent and potentially actionable SNVs in genes, KDR and KIT. The KDR variant rs1870377 in particular, has also been reported in other malignancies and has been variably associated with angiogenesis modulation and shorter survival (18-22). In this study, we explored the role of these KDR and KIT variants as well as TP53 mutations in vascular proliferation and aggressiveness in diffuse gliomas.
Materials and Methods
Patients and tumour samples. For this study, a total of 100 cases were assessed. These included biopsies and resected tumour samples of patients who were operated on at Imperial College Healthcare NHS Trust, London, from 2016 to 2019. The histological and integrated molecular diagnoses were established according to the 2016 WHO Classification of CNS tumours and cIMPACT guidelines using standard protocols of the Department of Cell Pathology of Imperial College (2, 3). Additionally, the MGMT promoter methylation status was assessed in all GBM cases by methylation – sensitive high-resolution melting analysis at the Department of Neuropathology of the Institute of Neurology, University College London. We identified 4 IDH mut GBM and these were excluded from MVD and survival analyses as mut IDH significantly influences prognosis in GBM (23). Patient samples from 2016-2017 did not have NGS analysis as the routine genetic profile was obtained using Sanger sequencing coupled with PCR and qPCR to identify variants of interest. In these cases, SNVs of interest were assessed as follows (see genotyping). Follow-up data were obtained from clinical databases. The survival time was defined as the time from the first surgery to the date of death or last contact. Time to progression was defined as the time between initial diagnosis and progression on neuroimaging. The data was provided from samples used for diagnostic purposes and were de-identified to the researchers.
Sequencing. Samples from 2018 were sequenced using the IonTorrent S5 sequencer (ThermoFisher Scientific, Waltham, MA, USA) and cases from 2019 were profiled using MiSeq instrument (Illumina, San Diego, CA, USA). IonTorrent sequencing was performed at the Molecular Pathology Laboratory of Imperial College using the Ion AmpliSeq™ Cancer Hotspot Panel v2 (ThermoFisher Scientific, Catalogue no. 4475346), according to manufacturer's instruction. Illumina sequencing was performed at the clinical genomics lab of the Royal Marsden Hospital, London, using a QIAseq (Qiagen, Hilden, Germany) targeted DNA from a panel of 33 genes developed in-house according to manufacturer's instructions; see Table II for sequencing panel details.
NGS data analysis. NGS generated by IonTorrent sequencing is subject to in-house quality control parameters using IonReporter whilst NGS data generated by Illumina sequencing underwent quality control at the clinical genomics lab of the Royal Marsden Hospital. Unfiltered NGS reports were collected and a looser cut-off frequency of 0.1% was applied. COSMIC (24), PolyPhen (25), CKB BOOST (26) and OncoKB (27) were used to probe known and unknown variants of genes in the panel as well as the significance of the variant itself. The softwares Bioconductor (28), FactoMineR (29) and RStudio were used for all subsequent bioinformatic analysis.
Genotyping. Cases that had not received NGS or sequencing information was unavailable were genotyped for KDR and KIT with restriction enzyme digestion and qPCR respectively. The following primers were used to amplify the segment of DNA containing the KDR variant (rs18170377): KDR-F: 5’-GGAAGTCCTCCACACTTCTCC-3’ and KDR-R: 5’-GGTAGGCTGCGTTGGAAGTTA-3’. Touchdown PCR from 65°C to 55°C was used to amplify the region of interest. PCR products were visualised using QIAXcel (Qiagen). PCR products were digested with AluI (New England Biolabs, Ipswich, MA, USA) at 37°C for 1h. AluI specifically digests variant KDR (rs1870377). Digested PCR products were visualised using QIAXcel Advanced system (Qiagen). TaqMan™ SNP genotyping qPCR assay (ThermoFisher Scientific) was performed according to manufacturer's instructions using a QuantStudio 5 Real-Time PCR System (ThermoFisher Scientific) to detect KIT variant (rs3822214).
Artificial neural network (ANN). ANN were constructed using Python 3.7, TensorFlow 2.0 (30) and the functional Keras API (31) using the complete targeted NGS data of 48 GBMs. To predict MGMT promoter methylation we approached it as a binary classification problem where 0 and 1 represent the unmethylated and methylated states, respectively. Our ANN consists of a 34-node input layer sequentially followed by 3 hidden layers containing 17, 9 and 5 nodes, with ReLU activations and a final single node output layer with a sigmoidal activation. Drop-out layers with a value of 0.2 were included in-between each hidden layer to prevent overfitting. NGS data were augmented by rearranging the order of gene mutation and variant frequencies for each patient. Processed NGS data were split into an 80:20, training and testing dataset. The model was compiled with a binary cross-entropy loss function and an Adaptive Moment Estimation optimizer with a learning rate of 0.025. The number of nodes for each hidden layer and the number of hidden layers were determined through a series of optimisations and were used to design a feed-forward neural network. We systematically removed the mutation information of genes and measured whether the accuracy of the ANN changes in response to this and determined whether the model begins to overfit upon losing these data, by comparing loss values.
Immunohistochemistry. Formalin-fixed paraffin-embedded sections were stained using a ready-to-use CD34 antibody (Leica Biosystems, Newcastle upon Tyne, UK) using a Leica BOND autostainer and was visualised using a BOND Polymer Refine Detection kit (Leica Biosystems). For each biopsy, the three tumour areas showing the highest MVD were selected and the number of blood vessels was counted at 20× magnification, as previously described (32, 33). Biopsies in which the tumour area was small (<1 cm2) were excluded from the analysis.
Statistical analysis. Statistical analysis was carried out using R, version 4.0.0 and graphics were generated using the Bioconductor (28) and FactoMineR (29) packages in RStudio and Graph Pad Prism 7. We used one-way ANOVA and unpaired t-tests to determine the significance of gene variants on MVD and Log-rank test for survival analysis. Fisher's exact test was used to determine whether gene variants influences MGMT promoter methylation. X2 test was used to assess the distribution of variants across tumour grades.
Results
A total of 100 cases were included in this study. There were 52 IDH WT GBMs, the aforementioned 4 IDH mut GBMs, and 44 lower grade (26 grade 2 and 18 grade 3) IDH mut astrocytomas. Following exclusion due to suboptimal DNA quality and availability of NGS data, 48 GBM and 26 samples of grade 2 and 3 astrocytomas were qualified for generating a mutation landscape. Disparities in the number of cases for further analyses were due to the availability of information for certain criteria (e.g. TP53 mutation status). Table I provides an overview of the patients included in the study.
Glioma mutation landscape. The mutation landscape was produced using archived data from NGS of 2018 and 2019 glioma cases. There was a correlation between age and mutation burden, which is expected given the accumulation of mutations over time (34) (Figures 1A and B). KDR (50%), TP53 (47.9%), KIT (43.75%), EGFR (43.75%) and PIK3CA (31.25%) were the most frequently mutated genes in both GBM and grade 2-3 glioma groups (Figures 1A and B). All patients with the KDR variants carry the Q472H variant (rs1870377) and those with the KIT variants carry the M541L variant (rs3822214). Gender and MGMT promoter methylation appear to have no discernible impact on the mutations and vice versa.
Relationship between KDR Q472H and PIK3CA/KIT mutations. The associations described by the Manhattan clustering of gene mutations in Figure 1A can be clearly visualised through the performing Principle Component Analysis (PCA) (Figure 2). We opted to use Manhattan clustering as it is less influenced by outlier data (35). TP53, KDR, PIK3CA and KIT contribute most to the variance in mutation status between patients; all other genes are clustered towards the centre and have very little contribution, if any, to the variance between patients. The co-occurrence of PIK3CA and KIT is shown by their effect in the same direction in both principal components, while the effect of KDR in the opposite direction of both components, relative to PIK3CA and KIT, suggests that mutations in KDR tend not to co-occur with PIK3CA and KIT. This does not completely rule out the possibility of co-mutating, as seen in Figure 1A, albeit with low mutation frequencies (<5%) of PIK3CA.
KDR Q472H increases MVD in GBMs, and is associated with shorter survival. KDR is also referred to as vascular endothelial growth factor receptor 2 and has a defined role in promoting angiogenesis in both normal and cancer contexts, including gliomas (36, 37). We investigated whether the Q472H variant of KDR would affect angiogenesis by measuring MVD using IHC for the endothelial marker CD34 (Figure 3A and B). As expected, MVD was increased with tumour grade. The mean MVD for GBM, grade 3 and grade 2 astrocytoma were 191.1, 120.3 and 76.1, respectively (Figure 3C). In our grade 2-3 population, we did not find a significant difference in MVD in patients with KDR and KIT variants (Figures 3D and F) and, accordingly, we did not find any significant impact on the time to progression in patients harbouring these variants (Figures 3G and I). We observed a 67% increase in MVD only in patients with grade 3 astrocytomas and pathogenic TP53 variants (Figure 3E) though this does not appear to affect overall survival (Figure 3H). Conversely, in our GBM population, we found that KDR Q472H increased MVD by 35% (Figure 3J) and significantly reduced the probability of overall survival (Figure 3L); with Fisher's exact test, we found that this is not influenced by MGMT promoter methylation. KIT M541L had no impact on either parameter (Figure 3K and M). Enhanced angiogenesis in tumours is associated with irregular patterns of blood vessel formation (38). Two patterns of microvascular proliferation can be recognised, a classical one showing endothelial multilayering and a bizarre one with formation of highly atypical glomeruloid vascular structures (39, 40). We did not observe any significant changes in the frequency of bizarre vs. classical angiogenesis in relation to tumour grade, pathogenic TP53/PIK3CA or KDR Q472H and KIT M541L status.
Increased frequency of KIT M541L in grade 3 and 4 glioma. KIT M541L appears to have no significant association with survival in GBM or grade 2-3 astrocytoma patient populations (Figures 3F, I, K, and M), we observed a clear and significant shift in the frequency of the KIT M541L variant in our lower grade glioma group. Only 8% of our grade 2 patients were positive for KIT M541L vs. 44% of grade 3 and 26% GBM patients, as shown in Figure 4.
Mutation status of angiogenesis-related genes predict MGMT promoter methylation status. We next wanted to investigate whether the SNV detected with the NGS panel could be used for predictive measures. To this end, we first sought to use the sequencing information to predict the methylation status of the MGMT promoter. We constructed an ANN using the targeted NGS data of all 48 GBM cases which received NGS. Our ANN could predict the MGMT methylation status of a patient with 86% accuracy using sequencing information from NGS, after the model iterated over the data 14 times (Figure 5A). The model did not perform overfitting in any case as the loss values between training and test sets converged by the end of the training (Figure 5B). Next, we sought to investigate what information was necessary for the model to be able to accurately predict MGMT methylation status. We found that the removal of PIK3CA, KIT and KDR has the most significant reduction in accuracy, reducing it to 16% from 86% (Figure 5C) and resulted in overfitting (Figure 5J). Removal of solely PIK3CA did not affect accuracy suggesting insignificance, while removal of KIT or KDR reduced the accuracy to 72% and 57%, respectively (Figures 5D-F). Removing the mutation status of PIK3CA/KIT and PIK3CA/KDR reduced accuracy to 57% and 72%, respectively (Figures 5G and I). Reductions in model accuracy were not observed when removing the mutation status of the 31 other genes presented in Figure 1A. These findings suggest there is some level of correlation between PIK3CA, KIT, KDR and MGMT promoter methylation, previously unseen without deep learning techniques.
Discussion
This study set out to elucidate the importance and possible significance of non-routinely assessed oncogenes in the context of glioma. We found that the most frequent variants occurred in genes that are known to be involved in angiogenesis regulation: KDR, KIT, PIK3CA and TP53 (36, 37, 41-44). This clustering, as well as the heatmap, also suggested a redundant relationship between mutations in angiogenesis genes, KDR and KIT (45) in the GBM cohort, wherein KDR variants were predominantly present in the absence of KIT mutations or in low frequency of a KIT mutation.
Further inspection of the NGS data showed that all variants of KDR detected were of the Q472H variant while the KIT variants were KIT M541L. KDR Q472H is already known to be a germline variant (18, 46), and we observed it in 50% of our GBM group and 47% of grade 2-3 astrocytomas, greater than the 21% of individuals reported positive in the general population (47). Our cohort was too small for population analysis. In addition, we did not have a matching non-tumour cohort available to determine whether KIT M541L is a germline variant, given the high frequency it presented (43.75%).
VEGF signalling is central to tumour angiogenesis and its regulation is complex (48, 49). The KDR Q472H variant is regarded as non-pathogenic but has been previously shown to increase MVD in melanoma (18) and lung (19) cancer and to be associated with shorter survival in gastric carcinoma (20). The variant has also been reported in colorectal cancer (21) and head and neck squamous cell carcinomas (HNSCC) (22). We found a significant increase in MVD in the GBM group but did not observe this effect in gliomas grade 2 and 3. Hence, the Q472H variant was enriched in GBMs with increased MVD though it did not seem to affect vascular morphology. The mechanisms of KDR variants in this context are not clear, however, it has recently been demonstrated in vitro (50, 51) that GBM cells can secrete factors to increase human brain endothelial cell migration. If KDR Q472H represents increased activity or sensitivity of the receptor, then this may explain the enhanced MVD in GBM patients who tested positive for KDR Q472H. To put this in to context, this elevated MVD reflects abnormal angiogenesis, which represents a growth advantage permitting the increased access to nutrients (52) and tumour cell invasion (53, 54). Moreover, an association between this variant of KDR and the risk of glioma has been described previously, hence this abnormal angiogenesis may be a factor involved in initial tumorigenesis as commented by the original study (55). It is reasonable to hypothesise that patients with more angiogenic tumours may have a worse survival rate (56), hence patients with the Q472H variant may have poorer prognosis. We showed that GBM patients with KDR Q472H have an overall shorter survival, and this may be related to increases in MVD. MGMT methylation does not influence survival in this case as we did not observe any clear correlations between MGMT promoter methylation and Q472H status.
However, increases in MVD have been shown to act as a double-edged sword in HNSCC (22), and similar effects may be observed in gliomas. Changes in MVD may influence the efficacy of radiotherapy as a result of changes in the oxygenation of the tumour (52, 57). In fact, Tinhofer et al., (22) has demonstrated this by showing that KDR Q472H-positive HNSCC patients had a significantly better survival when receiving chemo-radiation vs. KDR WT patients. Chemo-radiation is standard practice of the treatment of high-grade gliomas, hence identifying KDR Q472H patients may identify those who may benefit greater from these regimens, or from anti-angiogenic therapies (14, 58).
Interestingly, we also found that pathogenic TP53 variants enhance MVD only in grade 3 astrocytomas. WT TP53 is a known suppressor of angiogenesis (44, 59) and, therefore, pathogenic TP53 may enhance angiogenesis as we showed in our astrocytoma grade 3 cases, though in our cohort, this enhancement did not affect survival. TP53 mutations are common in IDH mut astrocytomas (60). It is, thus, possible that pathogenic TP53 primes IDH mut astrocytomas for high level MVD to contribute to progression to GBM.
Transcriptomic analysis (46) has identified an immuno-suppressive phenotype in KDR Q472H tumours; subsequent blockade of KDR augmented the anti-tumour immune response in a mouse melanoma model. This may be relevant to the increasing number of studies surrounding the use of immunotherapy in gliomas (61) and may, in the future, identify those patients who would benefit greater from immunotherapy strategies. Although we did not assess the inflammatory infiltration in this study, it could be important to examine whether KDR Q472H is associated with a reduced inflammatory response in gliomas in the future.
Gain- and loss-of-function mutations in KIT have been implicated in numerous types of cancer (62) in tumour cell proliferation, including cancer stem cell proliferation (63), as well as endothelial proliferation in gliomas, supporting tumour-associated angiogenesis (41, 64). In our study, patients with KIT mutations carry M541L, a variant that has been shown to increase the affinity of the receptor to its ligand, stem cell factor (SCF) (65, 66). This effect on SCF binding should in theory increase the sensitivity to imatinib, yet there are conflicting reports of this observation in clinical settings (67, 68). In our study, we observed that the presence of KIT M541L does not significantly affect MVD in the context of gliomas. Interestingly we observed the majority of KIT M541L positive cases in grade 3 and 4 gliomas. This variant has been found in patients with gliomas (69), however, its significance is not known. It is possible that the KIT M514L mutation affects aggressiveness with a different mechanism that does not involve angiogenesis regulation. It is unknown whether KIT M541L could serve as a potential diagnostic marker of high-grade glioma or faster progression to higher grade for earlier intervention strategies, and this should be verified in a larger cohort of low-grade glioma patients.
We used artificial intelligence (AI) to search for further insights into the relationship between the mutation status of oncogenes in the panel and clinically relevant features. We decided to build an ANN to predict MGMT promoter methylation in GBM as it is used to plan treatment strategies. Temozolomide is the gold-standard therapy of gliomas (70) but its activity is circumvented through MGMT promoter methylation (71). We proposed the use of targeted NGS data fed into an AI model to predict MGMT methylation, and we achieved this with an 86% accuracy. The model is unable to progress with further than 86% accuracy, but issues regarding accuracy could be resolved by providing a larger set of data to train the system on. Interestingly, in our model the variant status of KIT, KDR and PIK3CA – a signal transduction component in the activation of KIT (42) and KDR (43), are necessary to predict MGMT promoter methylation. Even though the link between angiogenesis-related genes and MGMT promoter methylation is not clear, a similar association between MGMT expression and decreased angiogenesis has been described experimentally in GBM, in vitro (72).
In conclusion, using NGS we identified the presence of KDR Q472H in a large proportion of our glioma patient cohort and linked it to an increase in MVD and poorer survival. As a germline mutation, this variant may predispose individuals to more aggressive tumours. Based on this study we suggest that the KDR Q472H status should also be included during reporting after these data have been verified in a larger study cohort. Reporting KDR Q472H status may identify patients who may incur greater benefits from more aggressive radiotherapy treatment regimens or from the inclusion of anti-angiogenic therapies and immunotherapy strategies. Our study also highlights the potential diagnostic significance of a KIT variant, to assess grade as well as monitor the grade change of astrocytomas. Our findings support the importance of incorporating germline variants into the genetic profile of neoplasms and extending it to the diagnosis of gliomas.
Acknowledgements
The Authors would like to thank the members of the molecular genetics laboratory at Hammersmith Hospital, Imperial College Healthcare NHS Trust for their technical assistance with this work.
Footnotes
Authors' Contributions
JK and CL designed and supervised the project. PDP, SM, LG, LT performed sequencing and contributed to sequencing data analysis. NZ and SSD conducted genotyping for the validation cohort. SSD and CL performed histological examination. NZ, SSD, JK and CL analyzed data. NZ constructed the AI model and wrote the manuscript. All Authors have read and approved the manuscript.
This article is freely accessible online.
Conflicts of Interest
The Authors report no conflicts of interest.
- Received August 24, 2020.
- Revision received September 17, 2020.
- Accepted September 21, 2020.
- Copyright© 2020, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved