Mapping actionable pathways and mutations in brain tumours using targeted RNA next generation sequencing

Many biology-based precision drugs are available that neutralize aberrant molecular pathways in cancer. Molecular heterogeneity and the lack of reliable companion diagnostic biomarkers for many drugs makes targeted treatment of cancer inaccurate for many individuals. Identifying actionable hyperactive biological pathways in individual cancers may improve this situation. To achieve this we applied a novel targeted RNA next generation sequencing (t/RNA-NGS) technique to surgically obtained glioma tissues. The test combines mutation detection with analysis of biological pathway activities that are involved in tumour behavior in many cancer types (e.g. tyrosine kinase signaling, angiogenesis signaling, immune response, metabolism), via quantitative measurement of transcript levels and splice variants of hundreds of genes. We here present proof of concept that the technique, which uses molecular inversion probes, generates a histology-independent molecular diagnosis and identifies classifiers that are strongly associated with conventional histopathology diagnoses and even with patient prognosis. The test not only confirmed known glioma-associated molecular aberrations but also identified aberrant expression levels of actionable genes and mutations that have so far been considered not to be associated with glioma, opening up the possibility of drug repurposing for individual patients. Its cost-effectiveness makes t/RNA-NGS to an attractive instrument to aid oncologists in therapy decision making.


Introduction
Many different cancer types share similar driver mutations. Examples are loss of function mutations in tumour suppressor genes (e.g. p53 and CDKN2A) and DNA repair proteins, leading to genetic instability and loss of cell cycle control) [1][2][3], and activating mutations, amplifications or fusion events in proto-oncogenes [4,5]. Other shared features are related to micro-environmental effectors, such as altered metabolism as a result of hypoxia [6,7], induction of angiogenesis [8] and immune suppression [9].
Additionally, aberrations exist that are tumour-type specific. Examples are expression of hormone receptors in cancers of prostate, ovary and breast [10,11]; mutations affecting metabolism (isocitrate dehydrogenase mutations [p.IDH1-R132H] in glioma and acute myeloid leukemia [12]); and mutations affecting the PI3K and MAPK pathway (PIK3CA, p.KRAS-G12/G13 mutations in adenocarcinomas [13], p.BRAF-V600E in melanoma [14]). Such specificity is however never absolute. As an example, IDH-and BRAF mutations are sporadically found in other cancers too [15][16][17]. Detection of such relatively rare, and therefore a priori unexpected mutations in individual patients could lead to repurposing of precision medicines in basket trials, in which precision drugs are administered to patients based on DNA profiling [18][19][20][21].
A number of actionable biological pathways in cancer involve the products of genes that are not mutated, but epigenetically regulated, for example by altered transcription factor availability, repressor activity or gene methylation, [22,23]. Activity of such pathways cannot be directly inferred from DNA analyses. Whole genome methylation analysis has robust diagnostic power [24] but does not allow analysis of activity of biological pathways, involved in cancer development and progression. An example is angiogenesis, initiated by hypoxiainducable factor (HIF-1α)-induced expression of an abundance of growth factors and followed by extensive crosstalk between tumour cells, tip-and stalk endothelial cells and pericytes [25,26]. DNA analysis also does not provide information on post-transcriptional events. For example, expression of alternative splice variants of vascular endothelial growth factor (VEGF-A) has implications for the regulation of angiogenesis [27] and splice variants of receptor tyrosine kinases can lead to auto-active and oncogenic PI3K signaling (e.g. EGFR VIII and MET Δ7-8 in glioma and MET Δ14 in lung cancer [28][29][30][31]).
A comprehensive overview of gene expression levels and alternative splice variants can be obtained with whole RNA next generation sequencing (w/RNA-NGS), provided sufficient coverage to detect alternative exonexon boundaries. w/RNA-NGS is increasingly performed in a research setting but is cost-wise still not suitable for implementation in routine patient care. There is therefore a huge need for novel and cost-effective methods to obtain clinically actionable and reliable information for individual patients, to be able to implement personalized treatment approaches.
Due to its low incidence (6 per 100,000) and high molecular heterogeneity [32], glioma is a difficult tumour type to organize clinical trials with, although the molecular underpinnings of gliomagenesis and glioma progression are relatively well established [12]. In the absence of alternatives, treatment of its most malignant form, glioblastoma, is still confined to palliative surgery, followed by chemotherapy with temozolomide (TMZ) and radiotherapy [33] which extents median life expectancy with only few months. Surgical cure for this tumour type is not possible due to its diffuse infiltrative nature [34]. Glioma is therefore one of the most challenging tumours for which new treatment strategies are urgently needed.
We here analyzed test and validation cohorts of in total 103 surgically derived brain tumours with quantitative targeted RNA next generation sequencing (t/RNA-NGS) [35][36][37][38]. The technique uses single molecule molecular inversion probes (smMIPs) and sensitively and quantitatively measures expression levels of and mutations in actionable genes. We show that t/RNA-NGS provides a histology-independent molecular diagnosis and identifies classifier transcripts that are closely associated with histopathological diagnosis and prognosis. By measuring hyperactivity of cancer-related pathways the test may also stratify individual patients for treatment with appropriate medicine.

Patients
The study described here was performed with brain tumour tissue from newly diagnosed patients who were operated for a glioma between 2013 and 2018 (n = 103). The cohort was separated in a training cohort (n = 75, tumours that were operated between 2013 and 2017) and an independent validation set of 28 tumours (operated in 2017 and 2018). Researchers were blinded to histopathology and clinical outcome. The study protocol was approved by the Ethical Committee for Human Experimentation of the Radboudumc. All patients signed informed consent. Directly after surgery, tissue samples were snap frozen in liquid nitrogen and stored at -80°C until further processing. In retrospect, patient characteristics were extracted from Radboudumc electronic patient files (EPIC) and documented in the electronic data capture system CASTOR. Histopathology and molecular diagnoses were extracted from the Dutch Pathology archive PALGA and summarized in Table 1.

RNA preparation and cDNA synthesis
Cryosections of 10 μm were cut for RNA isolation using TRIzol (ThermoFisher Scientific, Waltham, MA). For every sample a 4 μm serial section was stained with H&E to estimate percentage tumour area by an experienced neuropathologist (BK). RNA was reverse transcribed into cDNA using random hexamer primers and Superscript II (Invitrogen, CA) according to standard protocols. In parallel, tissue was processed to formalin-fixed paraffinembedded (FFPE) tissue blocks for routine diagnosis. Samples from 2014 onwards were also subjected to genetic analysis using targeted DNA next generation sequencing [18].

T/RNA-NGS with smMIPs
All enzymes were from NEB (Ipswich, MA) unless stated otherwise. The procedure of smMIP-based targeted RNA Next Generation Sequencing (t/RNA-NGS) to detect expression of metabolic genes has been described before [35,[37][38][39]. SmMIPs were designed by an adjusted version of the MIPgen algorithm [40] and were ordered from Integrated DNA Technologies (Leuven, Belgium). The set of smMIP probes used in [35] was expanded with probes for detection of regions of interest (ROI) in novel transcripts of interest that play a role in a wide variety of cancer types. Furthermore, smMIP probes were included for specific detection of splice   [37]. Probes for detection of Vascular Endothelial Growth Factor (VEGF) isoforms 121, 165 and 189 were designed with extension and ligation probes in exons 5/8, 5/7 and 5/6, respectively. Additional smMIPs were included that were directed against shared exons of transcript variants. For transcripts in which relevant mutations can be expected, smMIPs were selected with extension and ligation probes flanking the respective ROI. At least 5 smMIPs were used for each transcript of interest, evenly distributed along the transcript. cDNA was subjected to smMIP capture as described in [35,36]. In short, a mixture of 936 phophorylated smMIPs was hybridized to 15-50 ng of cDNA. After overnight hybridization and enzymatic gap-filling (KlenTaq polymerase, Epicentre, Madison, WI) and ligation (Ampligase, Epicentre) circular smMIPs are formed, the number of which is linearly related to the number of RNA molecules present in the original sample. After treatment with exonuclease to remove cDNA and linear smMIPs, the libraries of circular smMIPs were subjected to PCR with a unique barcoded primer set per sample. The PCR products of the expected length of 266 bp were purified with Ampure beads (Beckmann Coulter Genomics, High Wycombe, UK) and quantified on a TapeStation 2200 (Agilent Technologies, Santa Clara, CA). PCR libraries were sequenced on the Illumina Nextseq platform (Illumina, San Diego, CA) at the Radboudumc sequencing facility (output 2 × 151 bases). The barcode in the PCR primers allows for pooling of multiple samples into one sequencing library, followed by demultiplexing of reads, generating for all individual samples a FASTQ file.

Data processing
For each patient sample, all reads in the FASTQ file were mapped against reference transcripts (UCSC human genome assembly hg19) using the SeqNext module of JSI SequencePilot version 4.2.2 build 502 (JSI Medical Systems, Ettenheim, Germany). In all smMIP molecules a unique random 8 N sequence is included (unique molecular identifier or UMI), allowing to assign all identical sequencing reads to one originating circularized (unique) smMIP. This excludes PCR amplification bias and makes the assay quantitative. Number of unique reads for each individual smMIP in a sample were normalized to the total of unique reads in that sample and expressed as fragment per million (FPM). Mean FPM values from all  [42]. Mapped reads were filtered and deduplicated using sambamba v0.6.6 and feature quantification was performed using featureCounts v1.5.0-p1 against the RefSeq database (downloaded from the UCSC genome browser on 01/05/2017). Gene and exonlevel fragments per kilobase per million mapped reads (FPKMs) were calculated using a custom python script.

T/RNA-NGS profiles have prognostic value
To investigate the prognostic value of t/RNA-NGS, we profiled a training set of 75 brain tumours, including 69 grade II-IV diffuse gliomas and 6 brain lesions that upon routine histopathology were diagnosed as an ependymoma (n = 1), one dysembryoplastic neuroepithelial tumour (DNET), one pleomorphic xanthoastrocytoma, a variant glioma, a brain metastasis of lung adenocarcinoma and a lymphoproliferative disorder (LPD) ( Table 1). Annotated unique smMIP counts for each tumour sample ranged from 275,000 to 1,111,000 (not shown). Hierarchical unsupervised agglomerative clustering of the gene expression data of the grade II-IV gliomas (excluding the 6 rare cancers) resulted in 3 main clusters A, B and C, comprising of 26, 38 and 5 tumours, respectively (Fig. 1).
In a next step we performed a Wilcoxon Mann-Whitney U test to find genes that were differentially expressed between clusters A and B; A and C; and B and C. A total of 83 genes were differentially expressed between clusters A and B with p < 0.05 and False Discovery Rate (FDR) < 0.05 (Additional file 3: S2a). Among these were transcripts encoding transporters and enzymes involved in metabolism, as described before [39]. Membrane receptor tyrosine kinases (RTKs) NTRK2, ERBB3, ERBB4 were higher expressed in cluster B with average fold changes (FC) of 3.5, 4.8 and 5.2 respectively, while cluster A was characterized by significantly higher expression levels of EGFR vIII (FC = 250) and VEGF-A isoforms VEGF-A 121 , VEGF-A 165 and VEGF-A 189 (all with FC > 11, Additional file 2: Table S1a). Between cluster B and C 69 genes were differentially expressed (Additional file 2: Table S1b). Cluster A and C significantly differed with respect to expression levels of 9 genes (Additional file 2: Table S1c). The functional significance of these differences were not subject of further investigation in this study.
We then coupled the profiles to survival data. As shown in Fig. 2a, cluster B gene expression profiles of gliomas were associated with good prognosis (median survival, defined as time between surgery and death, > 6 years; exact value could not be calculated for available follow-up time) whereas gene expression profiles in cluster A and C were associated with median survival of 467 and 135 days, respectively. Results were highly significant between clusters A and B (p < 0.0001), A and C (p = 0.0078), and C and B (p < 0.0001). A Kaplan-Meier curve with survival data including also the non-glioma patients is presented in Additional file 1: Figure SI. Gene expression levels (in FPM) were transformed to a z-score for each individual transcript. After generating the dendrogram and heatmap using the average clustering method in R, histopathology results and mutation status were added in retrospect (legend and upper annotation bars). Cluster A and C strongly correlate with a wild-type IDH status, while cluster B is strongly associated with the IDH1 R132H mutation (see also Additional file 4: Table S3). IDH2 mutations (R172K, R172M, and R172W) are annotated as IDH2-MT (black in the annotation bar) and all coclustered with the IDH1 R132H mutation. IDH1-other refer to variants other than the IDH1 R132H (see Table 1)

Mutation analysis
We next investigated whether we could identify nucleotide variations that significantly associate with the subgroups. To this end we included all identified sequence variants of potential clinical significance with a coverage of > 10% of the unique reads. This analysis revealed that the well known hotspot mutation in the metabolic enzyme isocitrate dehydrogenase 1 (p.IDH1-R132H) [12] was significantly associated with cluster B (p < 3 × 10 − 5 when compared to cluster A, see Additional file 3: Table  S2). Additionally a duplication in succinate dehydrogenase A (SDHA) of unknown significance was strongly associated with group B (p = 3 × 10 − 5 ). The test also identified known oncogenic mutations in mitochondrial IDH2 in four patients (two p.R172K, one p.R172C and one p.R172M mutation, Table 1). Gene expression profiles from these gliomas co-clustered with those of the p.IDH1-R132H gliomas in cluster B (Fig. 1, black in IDH-status annotation bar). All IDH1 R132 and IDH2 R172 mutations were in retrospect confirmed by routine genetic analyses (not shown). As described before, also other variants of IDH1 (p.V178I, p.Y183C) were identified by the assay [39] (Table 1).
We then analyzed the profiles in relation to survival with treatment and histopathology as additional parameters. Retrospective analysis of histopathology and clinical follow-up data revealed that all patients in clusters A and C received temozolomide (TMZ) and/or radiotherapy upon signs of tumour progression after surgery, suggesting that the differences in survival were not related to treatment. We therefore concentrated on patients in cluster B who were treated with TMZ (n = 9, astrocytomas), procarbazine/lomustine/vincristine (PCV) chemotherapy (n = 6, all 1p/19q codeleted oligodendrogliomas), both adjuvant therapies (n = 3), or did not receive additional treatment (n = 19). Of 2 patients adjuvant therapy status was unknown (n = 2). As shown in Fig. 2b, outcome for oligodendroglioma patients treated with PCV was better than for astrocytoma patients treated with TMZ (p = 0.04).
Whereas group B was dominated by IDH-mutated gliomas, this group also contained 6 IDH wild-type glioblastomas. Comparison of gene expression profiles of these gliomas with those of cluster A revealed that expression levels of all VEGF isoforms were significantly lower in the group B IDH wt gliomas (Additional file 4: Table S3). This suggests that these 6 tumours grouped with the IDH mut gliomas based on the lack of an angiogenic response.

T/RNA-NGS based molecular diagnosis
To investigate whether t/RNA-NGS profiles are associated with histopathology diagnoses, we performed a group-based analysis according to the WHO 2016 grading system [43] and compared grade II/III astrocytomas (n = 12), grade II/III oligodendrogliomas (n = 19) and glioblastomas (n = 38). Wilcoxon Mann-Whitney U tests identified 79 genes that were differentially expressed between diffuse grade II/III oligodendrogliomas and glioblastomas (Additional file 5: Table S4a). In the top 5 of differentially expressed genes were LDHA and BCAT1, genes that are known to be hypermethylated in low grade, IDH mut gliomas [36,44]. Expression levels of 50 genes were significantly different between grade II/III astrocytomas and glioblastomas (Additional file 5: Table  S4b), whereas expression levels of only 1 gene (ALK) differed significantly between astrocytomas and oligodendrogliomas (Additional file 5: Table S4c). As Fig. 2 Kaplan-Meier analyses of unsupervised clusters. a For patients in clusters A (n = 26) and C (n = 5) median survival was 467 days and 135 days, respectively. For cluster B (n = 38) the follow-up period was too short to determine median survival. Survival of patients in cluster B was significantly better than of clusters A and C (both p < 0.0001). Survival in cluster A was significantly better than C (P = 0.008). b Oligodendroglioma patients in cluster B (diagnosed according to WHO 2016 classification, and treated with PCV) performed better than astrocytoma patients in this cluster (all treated with TMZ) (p = 0.04). Two patients (N16-06 and N16-10) were excluded from analysis because no survival data was present expected, both oligodendrogliomas and astrocytomas were distinguished from glioblastomas by the p.IDH1-R132H mutation. EGFR vIII was expressed in 30% of glioblastomas (n = 13; mean FPM = 550, range 6-1450) (Fig.3a) and never in grade II/III gliomas, in good agreement with literature [45]. Expression of wild-type EGFR was found among all gliomas, but was more prominent in glioblastomas than in grade II/III gliomas (mean FPM = 612 vs. 126; p = 0.0061). Levels of EGFR and MET expression in glioblastomas were inversely correlated (Spearman R = -0.59, p < 0.0001; Fig. 3b). Not surprising given the strong association of IDH1 mutations and grade II/III gliomas with cluster B of Fig. 1, the supervised histopathology-based analysis again identified ErbB3, ErbB4 and TrkB (the product of NTRK2) as potentially targetable proteins in grade II/III gliomas (Fig. 3c).

Carbonic anhydrase 12 is a poor prognostic factor
To investigate whether prognostic factors other than IDH1-mutations can be identified by t/RNA-NGS, we performed subgroup analysis on tumour profiles from IDH wt patients with survival < 14 months and > 14 months (the median overall survival time of this group of patients). A Fisher's exact test on the sequence variations in these groups identified two different splice isoforms of CA12: a full-length variant 1 (CA12v1), and a variant 2 (CA12v2) that lacks exon 9. Expression levels of CA12v1 were associated with poor survival (Fig. 4a). Detection of this variant was based on one smMIP with both exon 8-9 and exon 8-10 boundaries in its ROI. This smMIP can detect both isoforms due to the small size of exon 9 of 33 nucleotides. Interestingly, CA12v1 was never detected in IDH mut gliomas. CA12v1 expression values higher than 50 FPM translated in poor prognosis as shown by Kaplan-Meier analysis (272 vs.1002 days from surgery to death, p = 0.0137) (Fig. 4b).

Validation of t/RNA-NGS data
We have previously shown that the quality of t/RNA-NGS data correlates well with w/RNA-NGS [35]. All IDH1 and IDH2 mutations that were found with regular molecular screening, were verified by t/RNA-NGS [36]. Interestingly, the broad set-up of our t/RNA-NGS assay allowed the frequent detection of androgen receptor (AR) in gliomas at relatively high levels. To confirm and validate these and other findings on the proteome level Fig. 3 Actionable gene expression in grade II/III astrocytomas, oligodendrogliomas and glioblastoma. a EGFR vIII was detected at significant levels in 30% of glioblastomas (GB, red) but not in grade II/III astrocytomas (A) and oligodendrogliomas (O) (note the log scale of the Y-axis, expression < 5 FPM was considered EGFR vIII -negative). b Wild-type EGFR expression is inversely correlated with expression of MET in glioblastoma (Spearman R = − 0.59; p < 0.0001). c Mean expression of ErbB3, ErbB4, and NTRK2 in GB (red) and grade II/III O/A (blue). Error bars represent the standard deviation we performed immunohistochemical analyses for EGFR, MET, CA12 and AR on a number of tumours with low or high FPM values for the corresponding genes. Transcript levels for EGFR, MET, CA12, and AR correlated with protein levels (Fig. 5). In one case we found that EGFR-transcript levels from frozen tumour tissue did not match with protein levels in FFPE blocks. This discrepancy could be attributed to tumour heterogeneity since the FFPE block contained both EGFR-positive and negative areas (Fig. 5a).
We also detected prominent expression of FOLH1, the gene encoding prostate specific membrane antigen (PSMA), in approximately 70% of the gliomas. To validate this finding we performed IHC analysis for PSMA on a number of tumours with either low or high FOLH1 FPM values. IHC revealed expression of PSMA on the blood vessels of tumours with high FOLH1 FPM values (Fig. 5f).

Actionable markers in individual brain cancers
We then investigated whether t/RNA-NGS can stratify patient tumours based on relatively high expression of actionable targets (arbitrarily defined as > 2 fold expression compared to the mean of all gliomas). Results of a selection of tyrosine kinases (EGFR, MET, ALK, AXL, KIT, RET, ROS1) for 12 gliomas are presented in Fig. 6a. In only one tumour high expression levels of CD2, CD3, CD4, CD8 CTLA-4 and PD-1 were found, suggesting an inflammatory phenotype (not shown). Many mutations were detected in the assay that are described as somatic mutations in the Cosmic database. Potentially interesting mutations that we found (based on FATHMM score) included p.RET-S1002R, p.EGFR-G719D, p.EGFR-A289D, PTEN mutations ad TP53 mutations (not shown).
In the pleomorphic xanthoastrocyoma (PXA) in our test cohort a BRAF V600E mutation was found which was confirmed via genetic testing (Table 1). This mutation has been described before in PXA [46]. These are mostly slow growing, low-grade gliomas but some patients will undergo malignant progression. Mutations in BRAF have been identified as poor prognostic factors for this tumour type [47,48]. Interestingly, patient 14-03 was operated for a recurrent tumour 4 years later (sample 18-13 in the validation cohort). Samples 14-03 and 18-13 were therefore re-analyzed with t/RNA-NGS in the same sequencing run. The recurrent tumour showed a 70-80 fold increase of expression of transcription factors ASCL1 and HOXA1 and tyrosine kinases FGFR3, HER3, HER4 and MET (Fig. 6b). Another interesting finding was high expression of estrogen receptor 2 (ESR2) in the recurrent, but not the primary tumour. Whereas the p.BRAF-V600E mutation was retained, the recurrent tumour showed an additional p.NF1-E1571* nonsense mutation (not shown). This mutation has been described as pathogenic in the Cosmic database (https://cancer. sanger.ac.uk/cosmic/mutation/overview?id=133079) (FATHMM score = 0.99).
Other observations that were of potential clinical relevance were that the brain metastasis of lung carcinoma expressed relatively high levels of MET, RON and CTLA-4 (Fig. 7C). This may suggest sensitivity to the MET/RON inhibitor BMS777607/ASLAN002, combined with immune checkpoint inhibition [49]. The ependymoma showed high levels of EGFR vIII (not shown). One glioblastoma expressed the MET Δ7-8 variant that we described previously [31]. In the DNET, considered as a benign grade I tumour, we found high expression levels of ErbB3, FGFR2, FOLH1 and AXL.

Validation cohort
To validate these results, we next performed t/RNA-NGS on an independent validation cohort of 26 analyzable glioma samples. p.IDH1-R132H mutations were found in 7 samples, all diagnosed as grade II/III astrocytomas and  (Table 1). Differentially expressed genes between IDH mut and IDH wt gliomas in the training set were confirmed in the validation cohort (Additional file 6: Table S5). EGFR vIII expression was found in 4 of 19 (21%) of IDHwt glioblastomas and not in p.IDH1-R132H cancers, confirming results in the test cohort. Thus, the test is robust and yields batch-independent clinically relevant information. Survival data in the validation set could not be interpreted due to short time of follow-up.

Discussion
We here describe the clinical application of a relatively novel and cost-effective multiplex next generation RNA sequencing assay for cancer pathway profiling. The assay generates expression profiles and sequence information of genes that have been identified by literature surveys as potentially important for diagnosis, prognosis and response to precision medicines in a variety of cancer types. We previously demonstrated that the test allows  [36,38]. To analyze clinically actionable pathways we expanded the assay with smMIPs to quantify expression levels of and detect mutations in the involved genes. Analysis of the cohort of 103 brain cancers, of which 97 well characterized gliomas, allowed robust group-based analysis, and confirmed high inter tumour molecular heterogeneity with respect to tyrosine kinase expression. Indeed, the biology in individual gliomas that drives progression is mostly a black box. Therefore, in the routine clinical setting molecular genetics is increasingly implemented as add-on to histopathology to aid the diagnostic process, but influence on therapeutic decision making is limited. Biomarkers that influence clinical management include IDH mutations, combined with loss of chromosome arms 1p and 19q, leading to a diagnosis of oligodendroglioma that responds well to PCV, whereas methylation of the MGMT promoter in glioblastoma predicts a better response to TMZ [50][51][52]. EGFR status is also determined, amplification of which, in combination with the EGFR vIII mutation, is associated with glioblastoma [53]. We here show that a single t/RNA-NGS test gives the same genetic information that can be derived from DNA analysis, but additionally provides information on gene activity from mutated but also non-mutated alleles. Thus, RNAprofiling is better suited than DNA-profiling to draw conclusions on activity of biological pathways, and thus may be more reliable as a guide to personalized treatment plans [54,55]. It must be noted that for this study we isolated RNA from tumour tissue, snap-frozen directly after surgery. Although t/RNA-NGS can also be performed on RNA isolated from FFPE tissue, we found that whereas the quality of t/RNA-NGS data is good if tissue is formalin-fixed directly after surgery, it decreases in case of large time lags between surgery and formalin fixation (data not shown).
Tumours with overexpression of a targetable tyrosine kinase will be resistant to the corresponding inhibitor if other kinases are present that converge on the same signaling pathways [56]. Insight in the entire repertoire of expressed tyrosine kinases, downstream signaling intermediates and counteracting phosphatases may therefore predict intrinsic resistance or rapid development thereof [57][58][59]. t/RNA-NGS provides such insight for individual glioma patients, and even can identify gene activities that have so far been associated with other cancer types, opening new avenues for targeted therapy. For example, expression of the prostate cancer marker PSMA on microvasculature of glioblastoma as we show here, may provide opportunities for tumour-vascular targeting [60][61][62]. Of interest, it was reported before that high expression levels of AR in glioblastoma translates into sensitivity to the AR antagonist enzalutamide [63]. Notably, Fig. 6 Actionable markers in brain tumours. a overview of expression of potentially actionable targets in a subset of patients, expressed as fold change compared to the mean of the entire group. For clarity only tyrosine kinases are shown for a small selection of patients. b comparative analysis of a primary pleomorphous xantoastrocytoma and its recurrence 4 years later. c shows a selection of actionable targets in sample 13-05, a metastasis of a lung cancer. Note high expression levels of the macrophage receptor RON, MET and CTLA-4. Also VEGF isoforms 121 and 189 are detected our data also confirmed the identification of CA12v1 as a marker of poor prognosis in glioblastoma patients [64]. CA12 is a transmembrane enzyme that is involved in intracellular pH homeostasis and extracellular acidification, providing a biological explanation for the association with poor survival [65]. Overexpression of CA12v1, but not CA12v2, in cell models leads to acidification of the extracellular milieu, whereas overexpression of C12v2 leads to increased oxygen consumption (manuscript in preparation). Because CA12 is expressed at low levels in normal brain, it is a interesting therapeutic target for the group of CA12-positive glioblastomas [66].
In conclusion we here present a novel t/RNA-NGS 'one size fits all' assay that can provide a molecular diagnosis and prognosis for glioma and other cancers, as well as generate a list of (actionable) gene expression levels and mutations for individual cancers. The work described in this manuscript was performed with RNA, isolated from optimally biobanked tissue. The small size of the RNA regions of interest (150 bases) that is required for successful smMIP capture, however also allows testing of RNA isolated from formalin fixed, paraffin embedded tissue blocks, provided that autolysis is prevented by proper and rapid processing. The low costs and high throughput capacity of the assay (up to 400 samples can be processed simultaneously due to barcoding) makes it highly interesting to further analyze biobanked tissues in future studies and identify classifier profiles that can predict prognosis and response to therapies. Thus, t/RNA-NGS may allow repurposing of drugs in an individualized manner.
Additional file 1: Figure S1 Additional file 2: Table S1a. Differential gene expression in cluster A vs. B. A total of 83 genes were differentially expressed between the two clusters. Only genes that were significantly different are shown. Mean gene expression values (FPM) values for the clusters are given. For significance: A Wilcoxon-Mann-Whitney test with multiple testing correction was performed. Values are significant when the p-value is lower than the False Discovery Rate (FDR). The cutoff for the FDR was < 0.05. Table  S1b: Differential gene expression in cluster B vs. C. A total of 69 genes were differentially expressed between the two clusters. Only genes that were significantly different are shown. Mean gene expression values (FPM) values for the clusters are given. For significance: A Wilcoxon-Mann-Whitney test with multiple testing correction was performed. Values are significant when the p-value is lower than the False Discovery Rate (FDR). The cutoff for the FDR was < 0.05. Table S1c: Differential gene expression in cluster A vs. C. A total of 9 genes were differentially expressed between the two clusters. Only genes that were significantly different are shown. Mean gene expression values (FPM) values for the clusters are given. For significance: A Wilcoxon-Mann-Whitney test with multiple testing correction was performed. Values are significant when the p-value is lower than the False Discovery Rate (FDR). The cutoff for the FDR was < 0.05.
Additional file 3: Table S2. Mutation detection in glioma samples. A Fisher's exact test was performed to identify genetic mutations that distinguish the clusters A, B, and C as defined in Fig. 1. Mutations that were present in at least 10% of the unique reads were tested for significance. There were no significant differences between cluster A vs. C. The cutoff for the FDR was set at 0.01. Genes shown here were significantly different between clusters.
Additional file 4: Table S3. Differential gene expression in IDH IDHwt gliomas from cluster B vs cluster A. A total of 4 genes were differentially expressed between the cluster A and the IDH IDHwt gliomas that grouped in cluster B. Only genes that were significantly different are shown. Mean gene expression values (FPM) values for the clusters are given. For significance: A Wilcoxon-Mann-Whitney test with multiple testing correction was performed. Values are significant when the p-value is lower than the False Discovery Rate (FDR). The cutoff for the FDR was < 0.05.
Additional file 5: Table S4a. Differential gene expression in oligodendroglioma (O) vs. glioblastoma (G). A total of 79 genes were differentially expressed between the two histological types. Only genes that were significantly different are shown. Mean gene expression values (FPM) values for the histological type are given. For significance: A Wilcoxon-Mann-Whitney test with multiple testing correction was performed. Values are significant when the p-value is lower than the False Discovery Rate (FDR). The cutoff for the FDR was < 0.05. Table S4b. Differential gene expression in Astrocytomas (A) vs. glioblastoma (G). A total of 50 genes were differentially expressed between the two histological types. Only genes that were significantly different are shown. Mean gene expression values (FPM) values for the histological type are given. For significance: A Wilcoxon-Mann-Whitney test with multiple testing correction was performed. Values are significant when the p-value is lower than the False Discovery Rate (FDR). The cutoff for the FDR was < 0.05. Table S4c. Differential gene expression in Astrocytomas (A) vs. oligodendroglioma (O). One gene was differentially expressed between the two histological types. Only genes that were significantly different are shown. Mean gene expression values (FPM) values for the histological type are given. For significance: A Wilcoxon-Mann-Whitney test with multiple testing correction was performed. Values are significant when the p-value is lower than the False Discovery Rate (FDR). The cutoff for the FDR was < 0.05. Additional file 6: Table S5. Analysis of validation cohort. Shown are fold changes of gene expression values, calculated by dividing mean values of each transcript in the group of IDHwt tumors by that in the group of IDHmut tumors in the test cohort and validation cohort.