Next Article in Journal
Melanocortin-4 Receptor PLC Activation Is Modulated by an Interaction with the Monocarboxylate Transporter 8
Previous Article in Journal
Mechanism of Histone Arginine Methylation Dynamic Change in Cellular Stress
Previous Article in Special Issue
The Glioblastoma CircularRNAome
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Boundary of Exploratory Genomics and Translation in Sequential Glioblastoma

1
Szentagothai Research Center, University of Pecs, 20. Ifjusag Street, 7624 Pecs, Hungary
2
Department of Molecular Medicine, Markusovszky University Teaching Hospital, 5. Markusovszky Street, 9700 Szombathely, Hungary
3
Department of Pathology, School of Medicine, University of Pecs, 12. Szigeti Street, 7624 Pecs, Hungary
4
Office of the Dean, School of Medicine, University of Pecs, 20. Ifjusag Street, 7624 Pecs, Hungary
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(14), 7564; https://doi.org/10.3390/ijms25147564
Submission received: 11 June 2024 / Revised: 1 July 2024 / Accepted: 5 July 2024 / Published: 10 July 2024
(This article belongs to the Special Issue Molecular Genomics for Translational Research in Cancer)

Abstract

:
OMICS methods brought significant advancements to the understanding of tumor cell biology, which transformed the treatment and prognosis of several cancers. Clinical practice and outcomes, however, changed significantly less in the case of glioblastoma (GBM). In this study, we aimed to assess the utility of whole exome (WES) sequencing in the clinical setting. Ten pairs of formalin-fixed, paraffin-embedded (FFPE) GBM specimens were obtained at onset (GBM-P) and at recurrence (GBM-R). Histopathological and molecular features of all samples supported the diagnosis of GBM based on WHO CNS5. WES data were filtered, applying a strict and custom-made pipeline, and occurrence of oncogenic and likely oncogenic variants in GBM-P, GBM-R or both were identified by using the VarSeq program version 2.5.0 (Golden Helix, Inc.). Characteristics and recurrence of the variants were analyzed in our own cohort and were also compared to those available in the COSMIC database. The lists of oncogenic and likely oncogenic variants corresponded to those identified in other studies. The average number of these variants were 4 and 5 out of all detected 24 and 34 variants in GBM-P and GBM-R samples, respectively. On average, one shared oncogenic/likely oncogenic variant was found in the pairs. We assessed the identified variants’ therapeutic significance, also taking into consideration the guidelines by the Association for Molecular Pathology (AMP). Our data support that a thorough WES analysis is suitable for identifying oncogenic and likely oncogenic variants in an individual clinical sample or a small cohort of FFPE glioma specimens, which concur with those of comprehensive research studies. Such analyses also allow us to monitor molecular dynamics of sequential GBM. In addition, careful evaluation of data according to the AMP guideline reveal that though therapeutic applicability of the variants is generally limited in the clinic, such information may be valuable in selected cases, and can support innovative preclinical and clinical trials.

1. Introduction

Gliomas are the most common tumors of the central nervous system (CNS) and are traditionally subdivided based on histological grades I–IV, where grade IV glioma represents the most aggressive subtype, glioblastoma (GBM). During the last two decades, comprehensive sequencing of genomic DNA from solid tumors opened a new window for a better understanding of cancer cell biology by identifying characteristic somatic mutations that drive occurrence and progression in various histological entities, including GBM [1,2]. The fourth and fifth editions of the World Health Organization’s (WHO) classification of tumors of the CNS (WHO CNS4 and 5) integrated the most characteristic genomic and genetic changes into the routine histopathological work up [3,4].
While The Cancer Genome Atlas (TCGA) network mostly focused their OMICS analyses on gliomas at onset and proposed molecular subtypes of GBM [1,5,6], the Glioma Longitudinal Analysis Consortium (GLASS) carried out comprehensive genomic and epigenomic analyses to catalogue longitudinal molecular changes [7]. Wang et al. [8] also investigated genomic and transcriptomic data of longitudinal GBM and found that almost two thirds of patients experienced expression-based subtype changes over time. From the branching patterns and the evolutionary rates, the authors inferred that the relapse-associated clone existed years prior to the clinical diagnosis. They also estimated that 15% of tumors may have hypermutation at recurrence in highly expressed genes [8]. Neilsen et al. [9] established that deletions frequently occurred in CDKN2 A and CDKN2 B, activating mutations or amplifications were common in EGFR, and mutations were also frequent in TERT and PI3 K at both time points. However, the amount of the variants significantly varied in samples, suggesting clonal evolution over time [9]. Barthel et al. [10] noted that driver genes found at onset were typically retained at recurrence. Treatments using alkylating agents led to hypermutator phenotypes but at different rates among glioma subtypes. The study established that the strongest selective pressure affects glioma initiation and early development, while therapies affect the molecular profiles in a stochastic manner [10]. Deep whole genome sequencing by Körber et al. [11] revealed that de novo GBMs commonly carry chromosome 7 gain, 9 p loss or 10 loss early during tumor initiation, while TERT promoter mutation may arise later in association with rapid growth [11].
Using a customized enrichment/hybrid capture-based next-generation sequencing (NGS) gene panel, Sahm et al. [12] were able to identify potential treatment targets in various brain tumors using samples in routine neuropathology [12]. In longitudinal whole exome and whole genome analyses of tumors from a proband and a large cohort, Erson-Omay et al. [13] aimed to use genomics-guided personalized treatment to extend the survival of the proband with GBM [13]. Studying common drug targets in primary and recurrent GBM specimens, Schäfer et al. [14] found profound changes in the expressions of 9 out of 10 investigated molecules, with the only exception being FGFR1. This heterogeneity during the course of disease suggests that the molecular treatment design needs to be adjusted and personalized over time [14].
Cho et al. [15] used targeted NGS from primary and recurrent GBM to evaluate the mutational status of six DNA repair-related genes. The authors found both germline and somatic mutations, most frequently in MSH6 and POLE. The presence of MGMT (O6-methylguanine methyltransferase) promoter methylation and tumor mutational burden (TMB) were associated with mismatch repair (MMR) gene alterations, suggesting that evaluation of MMR genes for both biological and therapeutic considerations may be reasonable in GBM [15].
In this study, we aimed to conduct our WES studies using formalin-fixed, paraffin-embedded (FFPE) primary and recurrent GBM pairs from adult patients in the clinical setting, in order to define if this approach has real-life diagnostic value and may reveal information applicable to the care of individual patients. We used deep bioinformatic analyses for an enhanced and reliable identification of oncogenic and likely oncogenic variants in corresponding primary and recurrent GBM sample pairs, while also defining those variants that persist over time. The samples included here were partly subjects of our previous epigenomic studies [16,17,18].

2. Results

2.1. WES Coverage and Mapping Quality in the GBM-P and GBM-R Cohorts

The average coverage was 163 ± 63 SD, and the median coverage was 139 in the 10 primary (GBM-P) samples. The fraction of exome with at least 30× coverage was 77% and the median insert size was 132 bp. The mean mapping quality was 57.5, and the percentage of aligned mapped reads was 99%. The average coverage was 186 ± 89 SD and the median coverage was 160 in the 10 recurrent (GBM-R) samples. The fraction of exome with at least 30× coverage was 77% and the median insert size was 129 bp. The mean mapping quality was 57.5, and the percentage of aligned mapped reads was 98% (Table S1).

2.2. Genomic Variants in the GBM Sample Pairs

First, we determined the sum of variant numbers in affected genes in the 10 GBM sample pairs. We used a 553 gene-containing glioma panel (Table S2) and the COSMIC database filters to focus only on the glioma-associated variants among the other numerous variants detected by NGS. We excluded the benign and likely benign variants. Subsequently, three variant categories were distinguished based on the VarSeq Cancer Classifier algorithm and scoring system: oncogenic (5+ ≥ score), likely oncogenic (3+ ≥ score) and variant of uncertain significance (VUS) (2+ ≤ score). The software analyses were performed as “duos” of GBM-P and GBM-R samples, and variants were simultaneously evaluated and checked as to whether they exclusively appeared in one of the pairs (GBM-P or GBM-R) or in both (GBM-S). This approach resulted in three groups based on the appearance of variants. Variants that exclusively appeared in primary GBM samples, were sorted into the GBM-p, while those that exclusively appeared in recurrent GBM samples were sorted into the GBM-R category. The GBM-S category was created for variants shared between the primary and recurrent samples, thereby representing variants persisting over time. The average number of variants in the primary samples was 24, of which 1 was oncogenic, 3 were likely oncogenic and 20 were VUS. In the recurrent samples, we identified an average of 34 variants, of which 2 were oncogenic, 3 were likely oncogenic and 29 were VUS. The average number of variants occurring in both tumors was five, of which one was oncogenic, one was likely oncogenic and three were VUS (Table 1).

2.3. Oncogenic and Likely Oncogenic Variants in the 10 GBM Sample Pairs

In the next step, we analyzed the distribution of oncogenic and likely oncogenic variants in genes that appeared to be affected in more than one patient (Figure 1A–C). Detailed characteristics of all oncogenic and likely oncogenic variants (affected gene, variant allele frequency or VAF, read depth, alternative read count, cancer classifier, COSMIC mutation ID and functional significance) in the three groups (GBM-P, GBM-R, GBM-S) are summarized in Table S3.
Figure 1A–C show the distribution of oncogenic and likely oncogenic variants by gene (A), sample (B) and sequence ontology (C). Not surprisingly, the highest number of oncogenic and likely oncogenic variants were identified in the TP53 gene, followed by those in the KMT2 D, PTEN, APC, RB1 and CREBBP genes. Of note, RB1 and NF1 had only oncogenic variants, just like BRCA2 and RUNX1, while CREBBP, VHL, DNMT3 A, HNF1 A, NOTCH1 and BTK carried only likely oncogenic variants (Figure 1A). When we examined the distribution of the gene variants in the primary, recurrent and shared categories, curiously PTEN only had variants shared between the primary and recurrent samples. However, there were variants, such as those in APC and NOTCH1, which were only present in the primary samples or variants, such as those in PTPN11, BAP1, STK11 and NF1, which were exclusively present in the recurrent samples (Figure 1B). Regarding sequence ontology, the most common variant category was frameshift mutation followed by missense and splice site mutations. BRCA2, ARID1 A and KMT2 D were only affected by frameshift, while PTPN11, HFN1 A, BTK, EGFR and CREBBP only by missense and RB1 and NOTCH1 only by splice site mutations (Figure 1C).

2.4. Variants That also Frequently Occur in Tumors Other Than GBM

We selected variants that were found at least in 1 of our 10 GBM sample pairs and were identified in + ≥15 tumor samples in the COSMIC database. Of the nine highlighted variants, eight were oncogenic, while one was likely oncogenic. Six different variants of TP53 appeared in our samples, with the most prominent variant being the NM_000546.6:c.742 C>T, NP_000537.3:p.Arg248 Trp variant. This variant had 1044 entries in the COSMIC database. In addition to TP53, APC, ARID1 A and RUNX1 gene variants also frequently appeared in tumors other than glioma in the COSMIC database (Table 2).

2.5. Potential Therapeutic Targets Detected in the 10 GBM Sample Pairs

In GBM, currently no approved genetic variant-based therapy is available. We investigated the genes with recurrent oncogenic and likely oncogenic variants in the 10 GBM sample pairs at three levels. First, we checked if there is an available clinical study in any cancer for the identified recurrent genetic variants. Only the TP53 (NM_000546.6:c.742 C>T, NP_000537.3:p.Arg248 Trp; ClViC ID: EID4879) and the EGFR (NM_005228.5:c.787 A>C; NP_005219.2:p.Thr263 Pro; ClViC ID: EID4187) genes carried variants that have had recent preclinical studies. A subset of 58 cancer cell lines with the TP53 Arg248 Trp variant appeared insensitive to the MDM2 inhibitor AMGMDS3, while the EGFR Thr263 Pro mutant compared to the wild-type Ba/F3 cell line showed some sensitivity to the tyrosine kinase inhibitor (TKI) erlotinib [19,20]. Both aforementioned gene variants have 10+ > COSMIC entries. In our study, the TP53 variant appeared in the second and fifth, while the EGFR variant appeared in the sixth and ninth GBM pairs. The CREBBP gene also had a variant (NM_004380.3:c.4394 G>A, NP_004371.2:p.Gly1465 Glu) that appeared in more than one sample, namely the sixth, eighth and tenth GBM pairs. However, this latter variant had no targeting drug trial in any cancer (Table S4).
Second, we investigated if there were clinical trials targeting genetic variants in GBM, which were not present in our samples. Only studies targeting variants in the POLE (ClViC ID: EID1861) and RB1 (ClViC ID: EID1595) genes were found. In a GBM case study, a patient with the NM_006231.4:c.1270 C>G; NP_006222.2:p.Leu424 Val POLE variant was treated with pembrolizumab, an immune checkpoint inhibitor, which resulted in a radiographic response in the intracranial lesion. In a study on GBM cell lines, samples that carried the intact RB1 gene showed response to palbociclib, a selective inhibitor of CDK4 and CDK6 cyclin-dependent kinases, while samples with homozygous deletions or mutations causing loss of the Rb protein had no palbociclib response [21,22].
Third, for a more comprehensive approach, we considered what the Association for Molecular Pathology (AMP) guidelines [23] suggest for the evaluation of evidence concerning the clinical impact of variants. Successful targeting of certain gene variants falling into Tier IA (FDA-approved therapy available or included in profession guidelines as biomarker) and Tier IB categories (well-powered studies with consensus from field experts shows clinical significance), according to this guideline, may be classified as Tier IIC (FDA-approved therapies in other tumors or investigational therapies available; multiple small studies with consensus) or Tier IID category (preclinical trial or a few case reports without consensus) in gliomas. This consideration would significantly broaden the potential selection of clinical study targets from specific oncogenic or likely oncogenic variants detected in GBM to genes and variants successfully targeted in other tumor types. However, the overview of such extensive potential but not yet established targets is beyond the scope of our study.

2.6. Clonality and Tumor Mutation Rate (TMR) in the 10 GBM Sample Pairs

We assessed tumor heterogeneity by clustering variants based on VAF to estimate univariate density and cluster classification (Figure S1). Tumor heterogeneity is inferred by clustering VAFs (using mclust) in the primary and recurrent sample pairs. The MATH score is a simple quantitative measure of intra-tumor heterogeneity, calculated from the width of the VAF distribution. Higher MATH scores are found to be associated with higher clonal heterogeneity of samples. Based on this approach, clonal heterogeneity is slightly increasing in five and decreasing in another five of our sample pairs over time. Table 3 summarizes the results of the calculated tumor mutation rate per megabase for all samples including the total number of mutations, the number of mutations per megabase and the log10-transformed number of mutations per megabase. With the exception of sample pair 7, increases in TMR indicators can be seen in each sample pair. The TMR values used in this study represent an arbitrary measure to reflect changes in the mutation rates when primary and recurrent GBM pairs are compared and should not be confused with the widely used tumor mutation burden (TMB) (see details in the Materials and Methods section):

3. Discussion

The widespread use of NGS in studying large numbers of tumor samples in consortial research settings has identified the main molecular drivers of GBM occurrence and recurrence [3,4]. Due to great inter-tumor heterogeneity, determination of the individual molecular profiles may also be necessary. Therefore, we assessed what clinical utility in diagnostics and treatment approaches WES could provide. If we determine the top altered genes and their variants with oncogenic or likely oncogenic significance, can we infer druggable targets with direct clinical benefit for GBM in general and for a given patient in particular? Also, can we gain reliable information using FFPE GBM specimens for NGS?
Regarding technical aspects of NGS, studies have previously shown that quality and number of SNV calls in FFPE samples can be similar to those in freshly-frozen tissue and blood samples [24,25], and the Mutect2 program could detect a similar number of SNVs in FFPE and fresh frozen samples [25,26]. In accordance, we had a relatively high sequencing coverage in both GBM-P (163x ± 63 SD) and GBM-R (186x ± 69 SD) with acceptable DNA quality from our FFPE samples (Table S1). In addition to the recommended high mean coverage, other factors such as application of Mutect2 hard filters, the variant allele frequency (VAF) percentage or alternative read count cutoff and panels of altered genes obtained from tumor databases can greatly influence the variant yield and its reliability [12].
To extract reliable genomic information, we built a stringent filtering pipeline to keep the number of artifacts and false positive variants as low as possible. We used the Mutect2 program for SNV calls, created a filtering pipeline with the exclusion of strand bias and fragment-flagged variants, VAF ≥ 15% and a minimum alternative allele count of ≥20. In addition, we applied the COSMIC filter card, identified oncogenic and likely oncogenic variants based on the new VarSeq Cancer Classifier criteria system (2023.11.15, version 1.0) and implemented a custom-made glioma gene panel based on published data [27,28,29,30,31,32,33] (Table S2). This strict pipeline allowed us to extract clinically important and true genetic variants that presumably contribute to the formation, progression and recurrence of GBM.
In our cohort, the TP53 gene carried the most variants that appeared in more than half of the GBM pairs (Table S3). In addition, we detected one of the most frequent TP53 oncogenic variants (indicated in the COSMIC database) in two of our samples, which are currently being targeted in a preclinical study. In this trial, nevertheless, the TP53 Arg248 Trp variant appeared insensitive to the MDM2 inhibitor AMGMDS3 in a subset of 58 cancer lines (ClViC ID: EID4879) [19]. TP53 is one of the most commonly mutated tumor suppressor and oncodriver gene in tumors, including GBM. Its mutations disrupt the p53 pathway and induce oncogenic transformation [1,5,7,34,35]. Therefore, clinical monitoring and potential targeting of TP53 variants may be important elements of personalized medicine in the future.
Following TP53, the second-highest numbers of variants were found in the KMT2 D gene, which is a histone methyltransferase that modifies histone 3 at lysine 4 (H3 K4) by mono-methylation (H3 K4 me1) [36]. KMT2 D is frequently mutated in tumors, including GBM [37], but its normal expression is positively correlated with immune cell infiltration and negatively with tumor mutation burden. Dhar and Lee [38] described that most of the mutations found in KMT2 D are missense and truncating mutations, in concordance with our findings (Figure 1C) [38]. Malfunction of the gene product may lead to decreased protein expression, which can contribute to epigenetic changes and affect the mitochondrial metabolic profile [39]. These changes may facilitate tumor development and progression; therefore, KMT2 D may be a biomarker and a potential therapeutic target.
The third-highest numbers of variants were found in PTEN and APC, followed by SMARCA4, CREBBP, RB1, EGFR, ARID1 A and POLE. Similar to TP53, PTEN is also a known tumor suppressor and oncodriver gene that in a mutated form may contribute to tumor progression and recurrence [40]. In our cohort, the variants in this gene appeared exclusively in shared form (both the primary and recurrent samples) in almost half of the GBM pairs, and were mostly missense. While generally loss of function is the most common form of PTEN alterations, missense mutations can also disrupt its tumor-suppressing role and enhance malignant transformation [40]. The PTEN wild type and mutation status standalone as well as in combination with other molecular signatures can influence tumor survival time [40], suggesting this gene’s inclusion in clinical tumor panels to serve as a prognostic biomarker.
The APC gene contained mostly oncogenic variants, while the SMARCA4 and CREBBP genes most likely contained oncogenic variants. Mutations in CREBBP can promote cancer development and progression [41] and act as a tumor suppressor involved in DNA repair mechanisms (i.e., histone or p53 acetylation) [42]. In addition, the gene is part of the WNT pathway modulating glioma stem cell maintenance, differentiation and proliferation [18]. In contrast to the findings by Ellis et al. [43], where three CREBBP variants were VUS (one of them also in shared form) [43], our investigation revealed only oncogenic variants of solely missense types. Furthermore, variants of CREBBP appeared in all three groups (p,R,S), while one of the variants also appeared in three independent GBM pairs. In contrast to CREBBP, APC variants only appeared in GBM-p, while SMARCA4 variants were predominantly found in GBM-R samples. Altogether, the above results suggest that these genetic variants take part in the initiation and progression of GBM: APC in early tumor development and SMARCA4 in recurrence, while CREBBP in both stages. Similar to the CREBBP gene products, APC is also a member of the WNT pathway, while SMARCA4 can interact with elements of the WNT pathway, including CREBBP [44]. Altogether, variants in the above-mentioned genes may converge in an oncogenic pathway offering simultaneous targets in GBM.
Various mutations in EGFR, a tyrosine kinase receptor gene, are prominent drivers in the development and progression of several cancers, including GBM. In our study, oncogenic and likely oncogenic variants in this gene appeared in all three investigation categories (GBM-p, GBM-R and GBM-S). One of the EGFR variants (Thr263 Pro) recurred in two GBM sample pairs, and is currently under testing in a cell line-based preclinical trial (ClViC ID: EID4187) [21]. In consensus with previous data, all EGFR mutations in our cohort were missense variants, suggesting that this mutation type is the most common somatic EGFR SNV alteration [35,45]. Following the great clinical success in other cancers, experimental targeting of EGFR in gliomas has also been intensively studied in both preclinical and clinical settings, as explained later. It is also important to note that WES is not the ideal approach to detect the most common EGFR alterations, CNVs such as the EGFRvIII (exon 2–exon 7 deletion) variant or gene amplification in FFPE GBM samples. Due to the confounding effect of DNA fragmentation, we did not determine copy number alterations genome-wide in our specimens. We mention here for completeness the EGFR CNVs, along with the gene’s extensive therapeutic targeting in preclinical and clinical studies, though only with limited success thus far [46].
The mutation rate of the cell cycle and apoptosis regulator RB1 [47] was lower in our primary samples than those indicated in the TCGA [1]. Just as in the case of previously detailed genes, RB1 mutations also have tumor-promoting effects and disrupt normal cell biology in many cancers, including GBM [47]. Our findings are in agreement with those by Barthel et al. [10], stating that driver gene mutations in EGFR, NF1 and RB1 usually persist in the recurrent samples [10].
A review paper by the GLASS consortium [7] discerned various patterns of mutational profiles in primary and recurrent glioma samples. In one study, new mutations developed in the recurrent samples as a result of temozolomide treatment, while in another study, some recurrent tumors carried the same TP53 variants as the primary samples, while yet another study showed that variants in the recurrent tumors were driven by branched subclonal evolution and did not appear in the primary sample. When linear evolution was observed, the initial and recurrent tumor profiles matched [7]. Of note, from all the analyzed genes in the 10 GBM pairs, only the seventh patient lacked any oncogenic or likely oncogenic variants in the recurrent sample, or an increase in TMR (Table 1 and Table 3). Moreover, only this patient did not receive temozolomide that is known to alter the molecular evolution pattern, including occurrence of hypermutation in recurrent GBM [48]. Eskilsson and Verhaak [49] observed that some recurrent tumors deviated from their primary mutated gene patterns due to subclonal evolution and tumor heterogeneity, while others showed a linear evolution [49]. While our limited data from a small FFPE sample cohort do not allow us to infer a general clonal evolution pattern, they suggest both decreasing and increasing clonal heterogeneity over time (Figure S1). Variants that exclusively appeared in the primary or recurrent samples (i.e., APC, NOTCH1 in GBM-P and PTPN11, BAP1, STK11 and NF1 in GBM-R) probably are signs of branching clonal evolution patterns influenced, at least in part, by treatment effects [14], while the detection of shared variants aligned more with the linear evolution pattern (i.e., PTEN, EGFR, TP53).
The above data indicate that WES is suitable for the molecular characterization of individual clinical samples and is capable of identifying one or more drivers of tumor development. As these variants align in critical pathways of tumorigenesis, simultaneous targeting of multiple pathway elements might be considered for treatment development. Nonetheless, in this study we asked the question as to whether identification of oncogenic drivers in the clinical setting presently has direct utility in the treatment of a given patient with GBM.
Several approaches have been developed or are under development to block selected oncogenic drivers. The short list includes small molecular inhibitors, monoclonal antibodies (mAbs), specific vaccines and engineered immune cells. Regulatory RNAs (microRNA, siRNA, lncRNA etc.) and mRNA vaccines are relatively new and potent tools in tumor therapy because of the flexibility and adaptability of the techniques [50]. The mutation rate of TP53 is relatively high in GBM [8,11]; therefore, therapeutic targeting of this gene product is constantly under development [51]. Members of the tyrosine kinase receptor (RTK), PI3 K and RAS signaling pathways, such as EGFR, MEK and NF1, can be modulated by small molecular inhibitors and antibodies [50]. EGFR is one of the main driver molecules in GBM and has been targeted by CAR-T cells, vaccines (rindopepimut), small inhibitors, mAbs (cetuximab) and antibody–drug conjugates (Depatux-M, ABBV–221) in phase I, II or III trials. Though some approaches had been quite promising in the early phases, at the end, none of the studies have thus far yielded successful treatment in GBM, in contrast to other tumors [50]. EGFR RTK inhibitors have been particularly successful in non-small cell lung cancer, but appeared ineffective when tested in GBM. One could postulate that the reason might be that the compounds used in lung cancer inhibit the intracellular (tyrosine kinase) domain of the receptor where the mutations are also located, while the extracellular domain is mostly mutated in GBM. However, this mutated (e.g., EGFRvIII) receptor in GBM is also constitutively active through the phosphorylation of the down-stream signaling molecules by the intracellular tyrosine kinase domain. Therefore, the lack of efficacy of tyrosine kinase inhibitors in GBM may be related to reasons other than the different sites of mutations. Similarly, other drugs successfully targeting certain genes or variants in a variety of tumors may not be equally effective in GBM for many reasons related to the blood–brain barrier, the tumor microenvironment, molecular and biochemical properties of different tumors and host factors outside of the tumors. Therefore, until a clinical trial establishes the clinical benefit of targeting any given mutated gene or its oncogenic variant in GBM, inferring success from the promising outcome of trials in other cancers cannot be supported by the available data. Altogether, we may conclude that determination of the mutational profiles and driver genes in clinical GBM specimens is feasible by WES. However, presently it should be restricted to carefully selected patients for whom the standard medication has been exhausted and some benefit from a targeted therapy may be expected based on individual tumor characteristics.
Our study has some strength and weakness. We defined characteristic genomic changes and dynamics in sequential GBM sample pairs in the interface of research and translation. We used deep bioinformatics analyses to extract as many and as reliable variants relevant to gliomagenesis as possible. Our data support the utility of NGS in the diagnostic work up, even of individual GBM samples, but our analyses show as yet limited real-life therapeutic applicability. The use of DNA from FFPE samples and the applied bioinformatic analyses also involve some interpretive difficulties. The DNA modification by formalin greatly contributed to high numbers of false positive artefacts, but we made all efforts to bioinformatically minimize this effect by filtering out strand bias and fragment-flagged variants. As we studied DNA from a small cohort of FFPE samples, we used the COSMIC database and a glioma gene panel to identify the most relevant variants in GBM. This approach obviously prevented the discovery of new pathogenic variants, and made obligatory (if performed from the data of our small cohort) identifications of already known affected molecular pathways; however, our research did not aim for scientific discovery, rather to establish clinical utility of the variant detection. Overall, our analyses provided new information regarding molecular gliomagenesis and diagnostics, and support a careful as well as selective clinical applicability of this approach in clinical settings.

4. Materials and Methods

4.1. Subjects of the Study

Surgically removed FFPE GBM specimens were obtained between 1999 and 2017 at the Department of Pathology, University of Pecs (UP). Blocks left over from clinicopathological work up were used for the present molecular analyses after receiving approval (number: 7517 PTE 2018) from the Regional Clinical Research Ethics Committee of UP. After quality assessment, 10 pairs of immunohistochemistry screened, isocitrate dehydrogenase (IDH)−1 R132 H mutation negative initial (GBM-P) and recurrent (GBM-R) tumor blocks were included from patients with late onset disease [age range: 39–66 years] (Table 4).

4.2. Sample Characteristics

The diagnosis of primary (de novo) GBM was established based on standard clinical and histopathological criteria at the time when the study was initiated [3] and later upgraded according to the 2021 WHO CNS5 criteria [4]. Primary tumor specimens were taken before and recurrent specimens after chemoradiation.
In addition to histological features of GBM (mitotic index, microvascular proliferation, necrosis, atypia), molecular characteristics inferred from WES also supported the WHO CNS5-based diagnosis of GBM [4]. All GBM-P and -R specimens were IDH1/2 wild type by WES. Of the 10 GBM-P samples, 8 had TERT promoter mutation and 6 showed EGFR amplification. Of the 10 GBM-R samples, 9 had TERT promoter mutation and 7 had EGFR amplification. We identified homozygous CKDN2 A deletions in two and two of the ten (20%) GBM-P and GBM-R samples, respectively. In one case, where the primary tumor had a homozygous deletion, the corresponding recurrent sample contained a heterozygous deletion, while in another sample pair, the reverse of these changes was noted. Heterozygous deletions of the CDKN2 A gene were noted in one GBM-P and three additional GBM-R samples (Table S5). All samples satisfied the latest WHO CNS5 criteria for GBM diagnosis considering the histopathological and molecular characteristics [4] (Figure 2).

4.3. Sample Preparation and Quality Check

DNA was isolated from FFPE tissue sections received from the Department of Pathology of the University of Pecs (UP). The samples were deparaffinized with xylene and, after washing with alcohol, the QIAamp DNA FFPE tissue kit (Qiagen®, Hilden, Germany) was used for the isolation of DNA. In the final step, DNA was eluted from the columns in 50 µL of buffer and stored at −20 °C until further use. The concentration of the eluate was determined using the Qubit™ 1 X dsDNA HS assay kit on a Qubit 4 fluorimeter (Invitrogen, Carlsbad, CA, USA). Fragment analysis of the eluted DNA was performed using the Agilent Genomic DNA ScreenTape assay kit on an Agilent 4200 TapeStation system (Agilent Technologies, Santa Clara, CA, USA). In the primary samples, the average rate was 50.34 ± 13.91% SD between 200 and 2000 bp fragments, while the average rate was 36.76 ± 8.75% between 2000 and 60,000 bp fragments. In the recurrent samples, the average rate was 58.13 ± 11.75% SD between 200 and 2000 bp fragments, while the average rate was 29.61 ± 9.12% between 2000 and 60,000 bp fragments (Figure 2, Table S6).

4.4. TERT Promoter Sequencing

TERT promoter PCR products were sequenced using NGS to determine the mutational status of two well-known hotspots (C228 T and C250 T). The hotspot-containing region was amplified using GoTaq® Long PCR Master Mix (Promega, Madison, WI, USA) and primers based on the paper by Hafezi et al. [52]. Amplicon libraries were prepared using the NEBNext Ultra II DNA library prep kit (NEB, Ipswitch, MA, USA) and sequencing was performed on a NovaSeq 6000 instrument (Illumina, San Diego, CA, USA) (Figure 2).

4.5. Library Preparation for WES

Library preparation and exome capture were performed using Sure Select Human All Exon V7 (Agilent Technologies, Santa Clara, CA, USA). Briefly, 50 ng of DNA was enzymatically fragmented, end prepped, adaptor ligated and amplified. Hybridization was performed using Human All Exon V7 probes (Agilent Technologies, Santa Clara, CA, USA). Libraries were qualitatively and quantitatively assessed using TapeStation 4200 (Agilent Technologies, Santa Clara, CA, USA) and Qubit 3.0 (Invitrogen, Carlsbad, CA, USA) and were sequenced on a NovaSeq 6000 machine (Illumina Inc. San Diego, CA, USA) for 2 × 150 paired-end reads (Figure 2).

4.6. Bioinformatics

Raw FASTQ file generation, including basecall and demultiplexing, was carried out using bcl2 fastq (v2.20.0.422) [53] on a local HPC cluster. First, raw data quality was checked using FastQC (v0.11.9) [54]. Based on the results, datasets were filtered and adapters and low-quality bases/reads were removed with fastp (v0.21.0) [55]. Clean, high-quality sequences were aligned to the human reference genome (GRCh37), applying the bwa mem algorithm (v0.7.17) [56]. Default parameters were used for the chromosomes. The SAM/BAM file modifications (e.g., sorting, adding read group information, indexing) were carried out using various Picard Tools (v2.23.3) subcommands [57]. Before variant calling, base quality scores were recalibrated by the GATK BSQR module (v4.1.4.1) [58]. Somatic SNVs and short insertions and deletions (INDELs) were identified with the GATK Mutect2 algorithm [59]. Finally, raw variants were pre-filtered using the GATK FilterMutectCalls module [59] (Figure 2).

4.7. Filtering Pipeline for Variant Selection

VarSeq software version 2.5.0 (Golden Helix, Inc., Bozeman, MT, USA) was used for annotation and variant filtering (SNVs and INDELs). The GBM pairs were uploaded as duos.
First level of filtering: Hard filter fragment and strand bias-flagged variants were filtered out using Mutect2. Read depth was set to equal or higher than 50, and VAF was set to equal or higher than 15%. Variants with alternative (Alt) read counts lower than 20 were filtered out. Variants were removed with equal or higher than 2% allele frequencies (in 1 kG Phase3–Variants 5 a with genotype counts), equal or higher than 2% alternative allele frequencies (in gnomAD Exome Variant frequencies 2.1.1) and equal or higher than 5% minor allele frequencies (in NHLBI ESP6500 SI-V2-SSA137 Exome Variant Frequencies 0.0.30). Variants with dbSNP (dbSNP Common 155, NCBI) entries were also filtered out.
Second level of filtering: Variants flagged as initiator codon, intragenic, intron and synonymous were filtered out (Sequence Ontology, RefSeq Genes 105.20220307, NCBI). Functional interpretation of the variants was provided by dbNSFP in silico prediction tools (dbNSFP Functional Prediction Voting) (i.e., SIFT, PolyPhen2 HVAR, MutationTaster, MutationAssessor, FATHMM and FATHMM MKL), and variants that were predicted to be damaging, at least in one, were filtered out. Clinically nonsignificant variants (Benign, Likely Benign, VUS/Weak Benign) were filtered out based on the VarSeq Clinical Classification criteria system. Variants with loss of function (LoF), missense mutation, Other or Unknown predicted effect (RefSeq Genes 105.20220307, NCBI) and with at least one record in the COSMIC database (Cosmic Mutations 96, GHI) were retained. As a last step, we used the Mendel error filter card to extract variants with Shared, Unshared and de Novo Allele tags. The detectability of the variants depended on many factors including DNA fragmentation caused by the FFPE conservation technique, and the strict filter pipeline intended for avoiding artifacts. As a result of the latter, variants of decreased quality were filtered out in the primary or recurrent samples. However, in order to better capture disease-relevant variants, Varseq analyses were performed with “duos” of GBM-P and GBM-R variants, meaning that variants that passed the filters in one sample were specifically searched for in the corresponding sample pair. Using this approach, several variants were identified to be shared between the sample pairs (designated as GBM-S).
Third level of filtering: Next, variants were filtered using a custom-made glioma gene panel based on prior scientific publications [27,28,29,30,31,32,33]. Fusion, Oncogen and Tumor Suppressor variants were retained based on the COSMIC Cancer Gene Consensus 96, GHI database.
Fourth level of filtering: We manually curated and interpreted the remaining variants with VarSeq, VSClinical AMP (Figure 2, Table 5).

4.8. Identification of Oncogenic and Likely Oncogenic Variants in GBM-p, GBM-R and GBM-S

GBM-P and GBM-R NGS data obtained at filter level 2 were loaded into VarSeq as related duo samples. By this approach, we were able to search for and evaluate variants that appeared exclusively in the GBM-P and GBM-R categories, or appeared in both designated as shared category (GBM-S) (in other words, variants that persisted over time from primary to recurrent stages).

4.9. Tumor Mutation Rate (TMR)

TMR values were estimated by defining all variants (with at least 20× depth and after level 1 filtering in Table 5) per megabase. The capture size was set up to 42 MB. In addition to the total numbers of mutations used in the calculation, the numbers of mutations per megabase and the log10-transformed numbers of mutations per megabase were defined. The determination of TMR in this study differs from that generally used for the determination of tumor mutation burden (TMB). Due to the low number of samples and data points, the TMB values (typically calculated from oncogenic, likely oncogenic and VUS/weak oncogenic variants) would be too low in this study and would not be comparable with those calculated from data of large research cohorts. The purpose of TMR generated from the variants defined above is to reveal the changes in mutation rates from the primary to the recurrent stages of disease.

4.10. Tumor Heterogeneity

Variant allele frequency (VAF) values of all variants after level 1 filtering in Table 5 were used to calculate tumor heterogeneity by applying the default parametric finite mixture model (mclust) for clustering [60]. To measure intra-tumor heterogeneity, Mutant-Allele Tumor Heterogeneity (MATH) scores were calculated as well. The MATH score is a non-biased, quantitative method to assess intra-tumor genetic heterogeneity based on NGS-generated data. It is calculated as the ratio of the width to the center of distribution of mutant allele fractions among tumor-specific mutated loci [60,61]. MATH values of our GBM samples were calculated from the median absolute deviation (MAD) and the median of its mutant allele fractions at tumor-specific mutated loci using the built-in math.score function (https://rdrr.io/bioc/maftools/src/R/mathScore.R, accessed on 4 July 2024) from the maftools package [62], based on the following formula: MATH = 100 * MAD/median. Higher MATH scores are associated with higher heterogeneity of samples.

5. Conclusions

Our analyses reflect that the WES method is suitable, but with caution, for studying DNA from FFPE glioma specimens, as it is capable of providing valuable information regarding somatic oncogenic variants in smaller cohorts and individual samples in the clinical setting. In addition to supporting the latest WHO classifications of gliomas, this NGS-based method also allows us to monitor molecular tumor changes and identify the main molecular drivers at onset and at recurrence. The gained information can support patient selection for preclinical and clinical studies, or allow complementing standard treatment protocols in selected patients. Along with the published data, our results suggest that combined and personalized therapeutic interventions may be needed to successfully treat GBM. Therefore, clinical application of comprehensive (WES- or panel-based) molecular diagnostic approaches will likely gain wider roles in testing individual clinical GBM samples in the near future.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25147564/s1.

Author Contributions

Conceptualization: B.K. and M.T.; methodology: M.T., P.U., B.G., B.I.K., Z.K. and B.K.; software: B.G., M.T., A.G.; validation: P.U., B.G., B.I.K., A.G., B.K. and A.M.; formal analysis: B.G., M.T., B.K.; resources: B.K., B.G., A.G.; writing—original draft preparation: M.T. and B.K.; writing—review and editing: all coauthors. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Dr. Janos Szolcsanyi Research Fund by the University of Pecs, School of Medicine (Nr. KA-2019-42) granted to Bernadette Kalman. Project no. TKP2020–4.1.1 has been implemented with the support provided from the National Research, Development and Innovation Fund of Hungary, financed under the Thematic Excellence Program 2020 National Excellence Sub-program supporting the research of Bernadette Kalman. Bioinformatics analyses and infrastructure were supported by ELIXIR Converge (Nr. 871075). BG was supported by funding from the Horizon 2020 research and innovation program under the H2020-MSCA-COFUND Marie Skłodowska-Curie grant agreement (Nr. 754432) and by the Polish Ministry of Science and Higher Education, from financial resources for science in 2018–2023.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Regional Clinical Research Ethics Committee of University of Pecs (Number: 7517 PTE 2018). Informed consent was waived since all patients passed away by the time the study was initiated. The tumor specimens were FFPE blocks left over from clinical histopathological evaluation.

Data Availability Statement

Raw data are available at the European Nucleotide Archive (ENA) under PRJEB63610 accession.

Acknowledgments

The authors are most grateful to the patients whose postmortem tumors could be used for the research.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. The Cancer Genome Atlas Research Network. Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature 2008, 455, 1061–1068. [Google Scholar] [CrossRef] [PubMed]
  2. International Cancer Genome Consortium/The Cancer Genome Atlas Research Network. Pan-cancer analysis of whole genomes. Nature 2020, 578, 82–93. [Google Scholar] [CrossRef] [PubMed]
  3. Louis, D.N.; Perry, A.; Reifenberger, G.; Von Deimling, A.; Figarella-Branger, D.; Cavenee, W.K.; Ohgaki, H.; Wiestler, O.D.; Kleihues, P.; Ellison, D.W. The 2016 World Health Organization classification of tumors of the central nervous system: A summary. Acta Neuropathol. 2016, 131, 803–820. [Google Scholar] [CrossRef] [PubMed]
  4. Louis, D.N.; Perry, A.; Wesseling, P.; Brat, D.J.; Cree, I.A.; Figarella-Branger, D.; Hawkins, C.; Ng, H.K.; Pfister, S.M.; Reifenberger, G.; et al. The 2021 WHO Classification of Tumors of the Central Nervous System: A summary. Neuro-Oncology 2021, 23, 1231–1251. [Google Scholar] [CrossRef] [PubMed]
  5. Verhaak, R.G.W.; Hoadley, K.A.; Purdom, E.; Wang, V.; Wilkerson, M.D.; Miller, C.R.; Ding, L.; Golub, T.; Jill, P.; Alexe, G.; et al. Cancer Genome Atlas Research Network. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell 2010, 17, 98–110. [Google Scholar] [CrossRef] [PubMed]
  6. Noushmehr, H.; Weisenberger, D.J.; Diefes, K.; Phillips, H.S.; Pujara, K.; Berman, B.P.; Pan, F.; Pelloski, C.E.; Sulman, E.P.; Bhat, K.P.; et al. Cancer Genome Atlas Research Network. Identification of a CpG island methylator phenotype that defines a distinct subgroup of glioma. Cancer Cell 2010, 17, 510–522. [Google Scholar] [CrossRef] [PubMed]
  7. The GLASS Consortium. Glioma through the looking GLASS: Molecular evolution of diffuse gliomas and the Glioma Longitudinal Analysis Consortium. Neuro-Oncology 2018, 20, 873–884. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, J.; Cazzato, E.; Ladewig, E.; Frattini, V.; Rosenbloom, D.I.S.; Zairis, S.; Abate, F.; Liu, Z.; Elliott, O.; Shin, Y.-J.; et al. Clonal evolution of glioblastoma under therapy. Nat. Genet. 2016, 48, 768–776. [Google Scholar] [CrossRef] [PubMed]
  9. Neilsen, B.K.; Sleightholm, R.; McComb, R.; Ramkissoon, S.H.; Ross, J.S.; Corona, R.J.; Miller, V.A.; Cooke, M.; Aizenberg, M.R. Comprehensive genetic alteration profiling in primary and recurrent glioblastoma. J. Neurooncol. 2019, 142, 111–118. [Google Scholar] [CrossRef] [PubMed]
  10. Barthel, F.P.; Johnson, K.C.; Varn, F.S.; Moskalik, A.D.; Tanner, G.; Kocakavuk, E.; Anderson, K.J.; Abiola, O.; Aldape, K.; Alfaro, K.D.; et al. Longitudinal trajectories of diffuse glioma in adults. Nature 2019, 576, 112–120. [Google Scholar] [CrossRef] [PubMed]
  11. Körber, V.; Yang, J.; Barah, P.; Wu, Y.; Stichel, D.; Gu, Z.; Fletcher, M.N.C.; Jones, D.; Hentschel, B.; Lamszus, K.; et al. Evolutionary trajectories of IDHWT glioblastomas reveal a common path of early tumorigenesis instigated years ahead of initial diagnosis. Cancer Cell 2019, 35, 692–704. [Google Scholar] [CrossRef]
  12. Sahm, F.; Schrimpf, D.; Jones, D.T.; Meyer, J.; Kratz, A.; Reuss, D.; Capper, D.; Koelsche, C.; Korshunov, A.; Wiestler, B.; et al. Next-generation sequencing in routine brain tumor diagnostics enables an integrated diagnosis and identifies actionable targets. Acta Neuropathol. 2016, 131, 903–910. [Google Scholar] [CrossRef] [PubMed]
  13. Erson-Omay, E.Z.; Henegariu, O.; Omay, S.B.; Harmancı, A.S.; Youngblood, M.W.; Mishra-Gorur, K. Longitudinal analysis of treatment-induced genomic alterations in gliomas. Genome Med. 2017, 9, 12. [Google Scholar] [CrossRef]
  14. Schäfer, N.; Gielen, G.H.; Rauschenbach, L.; Kebir, S.; Till, A.; Reinartz, R.; Simon, M.; Niehusmann, P.; Kleinschnitz, C.; Herrlinger, U.; et al. Longitudinal heterogeneity in glioblastoma: Moving targets in recurrent versus primary tumors. J. Transl. Med. 2019, 17, 96. [Google Scholar] [CrossRef]
  15. Cho, Y.A.; Kim, D.; Lee, B.; Shim, J.H.; Suh, Y.-L. Incidence, clinicopathologic, and genetic characteristics of mismatch repair gene-mutated glioblastomas. J. Neurooncol. 2021, 153, 43–53. [Google Scholar] [CrossRef] [PubMed]
  16. Kraboth, Z.; Galik, B.; Tompa, M.; Kajtar, B.; Urban, P.; Gyenesei, A.; Miseta, A.; Kalman, B. DNA CpG methylation in sequential glioblastoma specimens. J. Cancer Res. Clin. Oncol. 2020, 146, 2885–2896. [Google Scholar] [CrossRef] [PubMed]
  17. Kraboth, Z.; Kajtár, B.; Gálik, B.; Gyenesei, A.; Miseta, A.; Kalman, B. Involvement of the Catecholamine Pathway in Glioblastoma Development. Cells 2021, 10, 549. [Google Scholar] [CrossRef]
  18. Tompa, M.; Kajtar, B.; Galik, B.; Gyenesei, A.; Kalman, B. DNA methylation and protein expression of Wnt pathway markers in progressive glioblastoma. Pathol. Res. Pract. 2021, 222, 153429. [Google Scholar] [CrossRef] [PubMed]
  19. Saiki, A.Y.; Caenepeel, S.; Cosgrove, E.; Su, C.; Boedigheimer, M.; Oliner, J.D. Identifying the determinants of response to MDM2 inhibition. Oncotarget 2015, 6, 7701. [Google Scholar] [CrossRef] [PubMed]
  20. Michaud, K.; Solomon, D.A.; Oermann, E.; Kim, J.-S.; Zhong, W.-Z.; Prados, M.D.; Ozawa, T.; James, C.D.; Waldman, T. Pharmacologic inhibition of cyclin-dependent kinases 4 and 6 arrests the growth of glioblastoma multiforme intracranial xenografts. Cancer Res. 2010, 70, 3228–3238. [Google Scholar] [CrossRef] [PubMed]
  21. Lee, J.C.; Vivanco, I.; Beroukhim, R.; Huang, J.H.Y.; Feng, W.L.; DeBiasi, R.M.; Yoshimoto, K.; King, J.C.; Nghiemphu, P.; Yuza, Y.; et al. Epidermal growth factor receptor activation in glioblastoma through novel missense mutations in the extracellular domain. PloS Med. 2006, 3, e485. [Google Scholar] [CrossRef] [PubMed]
  22. Johanns, T.M.; Miller, C.A.; Dorward, I.G.; Tsien, C.; Chang, E.; Perry, A.; Uppaluri, R.; Ferguson, C.; Schmidt, R.E.; Dahiya, S.; et al. Immunogenomics of hypermutated glioblastoma: A patient with germline POLE deficiency treated with checkpoint blockade immunotherapy. Cancer Discov. 2016, 6, 1230–1236. [Google Scholar] [CrossRef] [PubMed]
  23. Li, M.M.; Datto, M.; Duncavage, E.J.; Kulkarni, S.; Lindeman, N.I.; Roy, S.; Tsimberidou, A.M.; Vnencak-Jones, C.L.; Wolff, D.J.; Younes, A.; et al. Standards and guidelines for the interpretation and reporting of sequence variants in cancer: A joint consensus recommendation of the Association for Molecular Pathology, American Society of Clinical Oncology, and College of American Pathologists. J. Mol. Diagn. 2017, 19, 4–23. [Google Scholar] [CrossRef]
  24. Oh, E.; Choi, Y.-L.; Kwon, M.J.; Kim, R.N.; Kim, Y.J.; Song, J.-Y.; Jung, K.S.; Shin, Y.K. Comparison of accuracy of whole-exome sequencing with formalin-fixed paraffin-embedded and fresh frozen tissue samples. PLoS ONE 2015, 10, e0144162. [Google Scholar] [CrossRef] [PubMed]
  25. Robbe, P.; Popitsch, N.; Knight, S.J.; Antoniou, P.; Becq, J.; He, M.; Kanapin, A.; Samsonova, A.; Vavoulis, D.V.; Ross, M.T.; et al. Clinical whole-genome sequencing from routine formalin-fixed, paraffin-embedded specimens: Pilot study for the 100,000 Genomes Project. Genet. Med. 2018, 20, 1196–1205. [Google Scholar] [CrossRef] [PubMed]
  26. Kerick, M.; Isau, M.; Timmermann, B.; Sültmann, H.; Herwig, R.; Krobitsch, S.; Schaefer, G.; Verdorfer, I.; Bartsch, G.; Klocker, H.; et al. Targeted high throughput sequencing in clinical cancer settings: Formaldehyde fixed-paraffin embedded (FFPE) tumor tissues, input amount and tumor heterogeneity. BMC Med. Genom. 2011, 4, 68. [Google Scholar] [CrossRef] [PubMed]
  27. Jin, H.; Yu, Z.; Tian, T.; Shen, G.; Chen, W.; Fan, M.; He, Q.; Dai, D.; Zhang, X.; Liu, D. Integrative Genomic and Transcriptomic Analysis of Primary Malignant Gliomas Revealed Different Patterns Between Grades and Somatic Mutations Related to Glioblastoma Prognosis. Front. Mol. Biosci. 2022, 9, 873042. [Google Scholar] [CrossRef] [PubMed]
  28. Saadeh, F.S.; Morsi, R.Z.; El-Kurdi, A.; Nemer, G.; Mahfouz, R.; Charafeddine, M.; Khoury, J.; Najjar, M.W.; Khoueiry, P.; Assi, H.I. Correlation of genetic alterations by whole-exome sequencing with clinical outcomes of glioblastoma patients from the Lebanese population. PLoS ONE 2020, 15, e0242793. [Google Scholar] [CrossRef] [PubMed]
  29. Stasik, S.; Juratli, T.A.; Petzold, A.; Richter, S.; Zolal, A.; Schackert, G.; Dahl, A.; Krex, D.; Thiede, C. Exome sequencing identifies frequent genomic loss of TET1 in IDH-wild-type glioblastoma. Neoplasia 2020, 22, 800–808. [Google Scholar] [CrossRef] [PubMed]
  30. Noorani, I.; de la Rosa, J.; Choi, Y.H.; Strong, A.; Ponstingl, H.; Vijayabaskar, M.S.; Lee, J.; Lee, E.; Richard-Londt, A.; Friedrich, M.; et al. PiggyBac mutagenesis and exome sequencing identify genetic driver landscapes and potential therapeutic targets of EGFR-mutant gliomas. Genome Biol. 2020, 21, 26. [Google Scholar] [CrossRef]
  31. Garrett, A.-M.; Lastakchi, S.; McConville, C. The personalisation of glioblastoma treatment using whole exome sequencing: A pilot study. Genes 2020, 11, 173. [Google Scholar] [CrossRef] [PubMed]
  32. Shin, H.; Sa, J.K.; Bae, J.S.; Koo, H.; Jin, S.; Cho, H.J.; Choi, S.W.; Kyoung, J.M.; Kim, J.Y.; Seo, Y.J.; et al. Clinical targeted next-generation sequencing panels for detection of somatic variants in gliomas. Cancer Res. Treat. Off. J. Korean Cancer Assoc. 2020, 52, 41. [Google Scholar] [CrossRef]
  33. Shi, Z.; Li, K.K.; Kwan, J.S.H.; Yang, R.R.; Aibaidula, A.; Tang, Q.; Bao, Y.; Mao, Y.; Chen, H.; Ng, H. Whole-exome sequencing revealed mutational profiles of giant cell glioblastomas. Brain Pathol. 2019, 29, 782–792. [Google Scholar] [CrossRef] [PubMed]
  34. Kraboth, Z.; Kalman, B. Longitudinal characteristics of glioblastoma in genome-wide studies. Pathol. Oncol. Res. 2019, 26, 2035–2047. [Google Scholar] [CrossRef] [PubMed]
  35. McNulty, S.N.; Cottrell, C.E.; Vigh-Conrad, K.A.; Carter, J.H.; Heusel, J.W.; Ansstas, G.; Dahiya, S. Beyond sequence variation: Assessment of copy number variation in adult glioblastoma through targeted tumor somatic profiling. Hum. Pathol. 2019, 86, 170–181. [Google Scholar] [CrossRef] [PubMed]
  36. Dhar, S.S.; Zhao, D.; Lin, T.; Gu, B.; Pal, K.; Wu, S.J.; Alam, H.; Lv, J.; Yun, K.; Gopalakrishnan, V.; et al. MLL4 is required to maintain broad H3 K4 me3 peaks and super-enhancers at tumor suppressor genes. Mol. Cell 2018, 70, 825–841. [Google Scholar] [CrossRef] [PubMed]
  37. Luo, G.; Chen, G.; Chen, P.; Zhou, J. Pan-cancer analysis of histone methyltransferase KMT2 D with potential implications for prognosis and immunotherapy in human cancer. Comb. Chem. High Throughput Screen. 2023, 26, 83–92. [Google Scholar] [CrossRef] [PubMed]
  38. Dhar, S.S.; Lee, M.G. Cancer-epigenetic function of the histone methyltransferase KMT2 D and therapeutic opportunities for the treatment of KMT2 D-deficient tumors. Oncotarget 2021, 12, 1296. [Google Scholar] [CrossRef] [PubMed]
  39. Pacelli, C.; Adipietro, I.; Malerba, N.; Squeo, G.M.; Piccoli, C.; Amoresano, A.; Pinto, G.; Pucci, P.; Lee, J.-E.; Ge, K.; et al. Loss of function of the gene encoding the histone methyltransferase KMT2 D leads to deregulation of mitochondrial respiration. Cells 2020, 9, 1685. [Google Scholar] [CrossRef] [PubMed]
  40. Zhang, P.; Meng, X.; Liu, L.; Li, S.; Li, Y.; Ali, S.; Li, S.; Xiong, J.; Liu, X.; Li, S.; et al. Identification of the prognostic signatures of glioma with different PTEN status. Front. Oncol. 2021, 11, 633357. [Google Scholar] [CrossRef]
  41. Kung, A.L.; Rebel, V.I.; Bronson, R.T.; Ch’Ng, L.-E.; Sieff, C.A.; Livingston, D.M.; Yao, T.-P. Gene dose-dependent control of hematopoiesis and hematologic tumor suppression by CBP. Genes Dev. 2000, 14, 272–277. [Google Scholar] [CrossRef]
  42. Reed, S.M.; Quelle, D.E. p53 acetylation: Regulation and consequences. Cancers 2014, 7, 30–69. [Google Scholar] [CrossRef] [PubMed]
  43. Ellis, H.P.; McInerney, C.E.; Schrimpf, D.; Sahm, F.; Stupnikov, A.; Wadsley, M.; Wragg, C.; White, P.; Prise, K.M.; McArt, D.G.; et al. Clinically actionable insights into initial and matched recurrent glioblastomas to inform novel treatment approaches. J. Oncol. 2019, 2019, 4878547. [Google Scholar] [CrossRef] [PubMed]
  44. DiRenzo, J.; Shang, Y.; Phelan, M.; Sif, S.; Myers, M.; Kingston, R.; Brown, M. BRG-1 is recruited to estrogen-responsive promoters and cooperates with factors involved in histone acetylation. Mol. Cell. Biol. 2000, 20, 7541–7549. [Google Scholar] [CrossRef] [PubMed]
  45. Brastianos, P.K.; Nayyar, N.; Rosebrock, D.; Leshchiner, I.; Gill, C.M.; Livitz, D.; Bertalan, M.S.; D’Andrea, M.; Hoang, K.; Aquilanti, E.; et al. Resolving the phylogenetic origin of glioblastoma via multifocal genomic analysis of pre-treatment and treatment-resistant autopsy specimens. NPJ Precis. Oncol. 2017, 1, 33. [Google Scholar] [CrossRef] [PubMed]
  46. An, Z.; Aksoy, O.; Zheng, T.; Fan, Q.W.; Weiss, W.A. Epidermal growth factor receptor and EGFRvIII in glioblastoma: Signaling pathways and targeted therapies. Oncogene 2018, 37, 1561–1575. [Google Scholar] [CrossRef] [PubMed]
  47. Chkheidze, R.; Raisanen, J.; Gagan, J.; Richardson, T.E.; Pinho, M.C.; Raj, K. Alterations in the RB pathway with inactivation of RB1 characterize glioblastomas with a primitive neuronal component. J. Neuropathol. Exp. Neurol. 2021, 80, 1092–1098. [Google Scholar] [CrossRef]
  48. Yu, Y.; Villanueva-Meyer, J.; Grimmer, M.R.; Hilz, S.; A Solomon, D.; Choi, S.; Wahl, M.; Mazor, T.; Hong, C.; Shai, A.; et al. Temozolomide-induced hypermutation is associated with distant recurrence and reduced survival after high-grade transformation of low-grade IDH-mutant gliomas. Neuro-Oncology 2021, 23, 1872–1884. [Google Scholar] [CrossRef] [PubMed]
  49. Eskilsson, E.; Verhaak, R.G.W. Longitudinal genomic characterization of brain tumors for identification of therapeutic vulnerabilities. Neuro-Oncology 2016, 18, 1037–1039. [Google Scholar] [CrossRef] [PubMed]
  50. Liu, E.K.; Sulman, E.P.; Wen, P.Y.; Kurz, S.C. Novel therapies for glioblastoma. Curr. Neurol. Neurosci. Rep. 2020, 20, 19. [Google Scholar] [CrossRef] [PubMed]
  51. Kumar, A.; Chandna, S. Evidence for a radiation-responsive ‘p53 gateway’contributing significantly to the radioresistance of lepidopteran insect cells. Sci. Rep. 2018, 8, 2. [Google Scholar] [CrossRef] [PubMed]
  52. Hafezi, F.; Jaxel, L.; Lemaire, M.; Turner, J.D.; Perez-Bercoff, D. TERT promoter mutations increase sense and antisense transcription from the TERT promoter. Biomedicines 2021, 9, 1773. [Google Scholar] [CrossRef] [PubMed]
  53. bcl2 fastq (RRID:SCR_015058). Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed on 27 June 2023).
  54. Andrews, S. FASTQC: A Quality Control Tool for High Throughput Sequence Data. 2010. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed on 27 June 2023).
  55. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef] [PubMed]
  56. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
  57. Picard. Available online: http://broadinstitute.github.io/picard (accessed on 27 June 2023).
  58. Konstantin, O.; Ana, C.; Fernando, G.A. Qualimap 2: Advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics 2016, 32, 292–294. [Google Scholar] [CrossRef]
  59. Van der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ Data to High-Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline. Curr. Protoc. Bioinform. 2013, 43, 11.10.1–11.10.33. [Google Scholar] [CrossRef] [PubMed]
  60. Mroz, E.A.; Rocco, J.W. MATH, a novel measure of intratumor genetic heterogeneity, is high in poor-outcome classes of head and neck squamous cell carcinoma. Oral Oncol. 2013, 49, 211–215. [Google Scholar] [CrossRef] [PubMed]
  61. Rocco, J.W. Mutant allele tumor heterogeneity (MATH) and head and neck squamous cell carcinoma. Head. Neck Pathol. 2015, 9, 1–5. [Google Scholar] [CrossRef] [PubMed]
  62. Mayakonda, A.; Lin, D.-C.; Assenov, Y.; Plass, C.; Koeffler, H.P. Maftools: Efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 2018, 28, 1747–1756. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Oncogenic and likely oncogenic variants in our GBM samples. (A) depicts oncogenic and likely oncogenic variants in genes as determined by the VarSeq cancer classifier. (B) depicts the oncogenic and likely oncogenic variant distributions in the primary, recurrent and shared categories. (C) depicts the sequence ontology of detected variants within the affected genes.
Figure 1. Oncogenic and likely oncogenic variants in our GBM samples. (A) depicts oncogenic and likely oncogenic variants in genes as determined by the VarSeq cancer classifier. (B) depicts the oncogenic and likely oncogenic variant distributions in the primary, recurrent and shared categories. (C) depicts the sequence ontology of detected variants within the affected genes.
Ijms 25 07564 g001aIjms 25 07564 g001b
Figure 2. Workflow of laboratory preparations and bioinformatic analyses.
Figure 2. Workflow of laboratory preparations and bioinformatic analyses.
Ijms 25 07564 g002
Table 1. Variant numbers with different oncogenicity category across the 10 GBM pairs. GBM-P: Primary sample variants, GBM-R: Recurrent sample variants, GBM-S: Shared variants that appeared in both GBM-P and -R samples. O = Oncogenic, LO = Likely oncogenic, VUS = Variant of uncertain significance.
Table 1. Variant numbers with different oncogenicity category across the 10 GBM pairs. GBM-P: Primary sample variants, GBM-R: Recurrent sample variants, GBM-S: Shared variants that appeared in both GBM-P and -R samples. O = Oncogenic, LO = Likely oncogenic, VUS = Variant of uncertain significance.
GBM PairsGBM-P GBM-RGBM-S
SUMOLOVUSSUMOLOVUSSUMOLOVUS
12933232702259216
2120392610254202
3130585535473012
41210112434172002
56115556152545104
61311113534287124
73625291700171001
83432293428246303
9100191410134112
101611144132366114
Mean2413203423295113
Table 2. Variants in GBM also frequently detected in other tumors. The table shows variants that are frequent in the COSMIC database. Missense variant = M, Frameshift = F, Splice variant = SV, O = Oncogenic, LO = Likely oncogenic.
Table 2. Variants in GBM also frequently detected in other tumors. The table shows variants that are frequent in the COSMIC database. Missense variant = M, Frameshift = F, Splice variant = SV, O = Oncogenic, LO = Likely oncogenic.
Variant InfoGeneCancer ClassifierCOSMIC Mutations 96
Chr:PosRef/AltGene NameHGVS
cDot
HGVS
pDot
Seq. Ont.ScoreClass.Mutation IDCount
17:7578205C/TTP53NM_000546.6:c.644 G>ANP_000537.3:p.
Ser215 Asn
M9OCOSM4409329
17:7578262C/GTP53NM_000546.6:c.587 G>CNP_000537.3:p.
Arg196 Pro
M9OCOSM4381436
17:7577539G/ATP53NM_000546.6:c.742 C>TNP_000537.3:p.
Arg248 Trp
M9OCOSM106561044
17:7578291T/ATP53NM_000546.6:c.560–2 A>Tp.?SV4LOCOSM4502615
17:7578550G/TTP53NM_000546.6:c.380 C>ANP_000537.3:p.
Ser127 Tyr
M8OCOSM4397037
21:36171607G/ARUNX1NM_001754.5:c.958 C>TNP_001745.2:p.
Arg320 Ter
F7OCOSM4169921
1:27087503C/TARID1 ANM_006015.6:c.2077 C>TNP_006006.3:p.
Arg693 Ter
F6OCOSM18423637
5:112164616C/TAPCNM_000038.6:c.1690 C>TNP_000029.2:p.
Arg564 Ter
F9OCOSM1884882
17:7577548C/TTP53NM_000546.6:c.733 G>ANP_000537.3:p.
Gly245 Ser
M9OCOSM6932670
Table 3. Tumor mutation rates in GBM-P and GBM-R samples. The table shows the tumor mutation rates (TMRs) in the 10 GBM pairs. The TMR values were calculated by defining all variants (with at least 20× depth and after level 1 filtering, see details in Section 4.3) per megabase. GBM-P: Primary sample variants, GBM-R: Recurrent sample variants.
Table 3. Tumor mutation rates in GBM-P and GBM-R samples. The table shows the tumor mutation rates (TMRs) in the 10 GBM pairs. The TMR values were calculated by defining all variants (with at least 20× depth and after level 1 filtering, see details in Section 4.3) per megabase. GBM-P: Primary sample variants, GBM-R: Recurrent sample variants.
SampleTotal MutationsTotal Mutations per MbTotal Mutations per Mb (log10)
1. GBM-P15,227362.552.56
1. GBM-R25,619609.992.79
2. GBM-P11,898283.292.45
2. GBM-R16,759399.022.6
3. GBM-P6281149.552.17
3. GBM-R39,456939.432.97
4. GBM-P11,858282.332.45
4. GBM-R22,414533.672.73
5. GBM-P35,395842.742.93
5. GBM-R45,0601072.863.03
6. GBM-P9111216.932.34
6. GBM-R20,473487.452.69
7. GBM-P20,877497.072.70
7. GBM-R10,648253.522.40
8. GBM-P17,404414.382.62
8. GBM-R22,285530.62.72
9. GBM-P6119145.692.16
9. GBM-R7341174.792.24
10. GBM-P11,005262.022.42
10. GBM-R29,776708.952.85
Table 4. Clinical characteristics of patients. The table shows patients’ gender, age of onset of GBM, applied treatment and time to relapse in weeks.
Table 4. Clinical characteristics of patients. The table shows patients’ gender, age of onset of GBM, applied treatment and time to relapse in weeks.
GBM-P SamplesGBM-R SamplesGenderAge at Onset TreatmentTime to Relapse (Weeks)
UPL22–003804UPL22–003816man61Surgery + radio + TMZ49
UPL22–003805UPL22–003817man39Surgery + radio + TMZ40
UPL22–003806UPL22–003818man62Surgery + radio + TMZ58
UPL22–003807UPL22–003819woman61Surgery + radio + TMZ31
UPL22–003810UPL22–003822man66Surgery + radio + TMZ56
UPL22–003811UPL22–003823woman53Surgery + radio + TMZ55
UPL22–003812UPL22–003824woman63Surgery+ irradiation30
UPL22–003813UPL22–003825woman45Surgery + radio + TMZ143
UPL22–003814UPL22–003826man43Surgery + radio + TMZ135
UPL22–003815UPL22–003827woman56Surgery + radio + TMZ199
Table 5. VarSeq filtering pipeline. This table summarizes filter cards, databases and settings applied in four sequential steps when creating our analysis pipeline.
Table 5. VarSeq filtering pipeline. This table summarizes filter cards, databases and settings applied in four sequential steps when creating our analysis pipeline.
Filter Cards and DatabasesSettings
First level of filteringGATK Mutect2 hard filtersFragment and SB variants were filtered out
Read Depth (DP)≥50
Variant Allele Frequency (VAF)≥15%
Alternative Read Count≥20
Allele Freq (1 kG Phase3)≥2% or missing
Alternative Allele Freq ≥2% or missing
(gnomAD Exome Variant frequencies 2.1.1)
All Minor Allele Frequency (NHLBI 0.0.30)≥5% or missing
dbSNP Common 155 (NCBI)False
Second level of filteringSequence Ontology (RefSeq Genes 105.20220307, NCBI)Initiation codon, intragenic and synonymous variants were filtered out
dbNSFP Functional Prediction VotingFunctional interpretation of variants
Cancer ClassifierBenign, Likely Benign, VUS/Weak benign and missing variants were filtered out
Effect (RefSeq Genes 105.20220307, NCBI)LoF, missense, other and missing variants were retained
COSMIC (Cosmic Mutations 96, GHI) True
VarSeq built-in flag system Technically hard to detect variant extraction
Third level of filteringGeneID (Aux Fields RefSeq Genes 105.20220307, NCBI)553 gene glioma-specific panel
Fourth level of filteringVSClinical, AMP Manual variant interpretation
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tompa, M.; Galik, B.; Urban, P.; Kajtar, B.I.; Kraboth, Z.; Gyenesei, A.; Miseta, A.; Kalman, B. On the Boundary of Exploratory Genomics and Translation in Sequential Glioblastoma. Int. J. Mol. Sci. 2024, 25, 7564. https://doi.org/10.3390/ijms25147564

AMA Style

Tompa M, Galik B, Urban P, Kajtar BI, Kraboth Z, Gyenesei A, Miseta A, Kalman B. On the Boundary of Exploratory Genomics and Translation in Sequential Glioblastoma. International Journal of Molecular Sciences. 2024; 25(14):7564. https://doi.org/10.3390/ijms25147564

Chicago/Turabian Style

Tompa, Marton, Bence Galik, Peter Urban, Bela Istvan Kajtar, Zoltan Kraboth, Attila Gyenesei, Attila Miseta, and Bernadette Kalman. 2024. "On the Boundary of Exploratory Genomics and Translation in Sequential Glioblastoma" International Journal of Molecular Sciences 25, no. 14: 7564. https://doi.org/10.3390/ijms25147564

APA Style

Tompa, M., Galik, B., Urban, P., Kajtar, B. I., Kraboth, Z., Gyenesei, A., Miseta, A., & Kalman, B. (2024). On the Boundary of Exploratory Genomics and Translation in Sequential Glioblastoma. International Journal of Molecular Sciences, 25(14), 7564. https://doi.org/10.3390/ijms25147564

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop