Improved prognostication of glioblastoma beyond molecular subtyping by transcriptional profiling of the tumor microenvironment.

Glioblastoma (GBM), the most aggressive form of brain cancer, is characterized by a high level of molecular heterogeneity, and infiltration by various immune and stromal cell populations. Important advances have been made in deciphering the microenvironment of GBMs, but its association with existing molecular subtypes and its potential prognostic role remain elusive. We have investigated the abundance of infiltrating immune and stromal cells in silico, from gene expression profiles. Two cohorts, including in-house normal brain and glioma samples (n = 70) and a large sample set from TCGA (n = 393), were combined into a single exploratory dataset. A third independent cohort (n = 124) was used for validation. Tumors were clustered based on their microenvironment infiltration profiles, and associations with known GBM molecular subtypes and patient outcome were tested a posteriori in a multivariable setting. We identified a subset of GBM samples with significantly higher abundances of most immune and stromal cell populations. This subset showed increased expression of both immune suppressor and immune effector genes compared to other GBMs and was enriched for the mesenchymal molecular subtype. Survival analyses suggested that tumor microenvironment infiltration pattern was an independent prognostic factor for GBM patients. Among all, patients with the mesenchymal subtype with low immune and stromal infiltration had the poorest survival. By combining molecular subtyping with gene expression measures of tumor infiltration, the present work contributes with improving prognostic models in GBM.

Protein expression of BIRC5, TK1, and TOP2A in malignant peripheral nerve sheath tumours--A prognostic test after surgical resection.

No consensus treatment regime exists beyond surgery for malignant peripheral nerve sheath tumours (MPNST), and the purpose of the present study was to find new approaches to stratify patients with good and poor prognosis and to better guide therapeutic intervention for this aggressive soft tissue cancer. From a total of 67 MPNSTs from Scandinavian patients with and without neurofibromatosis type 1, 30 MPNSTs were investigated by genome-wide RNA expression profiling and 63 MPNSTs by immunohistochemical (IHC) analysis, and selected genes were submitted to analyses of disease-specific survival. The potential drug target genes survivin (BIRC5), thymidine kinase 1 (TK1), and topoisomerase 2-alpha (TOP2A), all encoded on chromosome arm 17q, were up-regulated in MPNST as compared to benign neurofibromas. Each of them was found to be independent prognostic markers on the gene expression level, as well as on the protein level. A prognostic profile was identified by combining the nuclear expression scores of the three proteins. For patients with completely resected tumours only 15% in the high risk group were alive after two years, as compared to 78% in the low risk group. In conclusion, we found a novel protein expression profile which identifies MPNST patients with inferior prognosis even after assumed curative surgery. The tested proteins are drug targets; therefore the expression profile may provide predictive information guiding the design of future clinical trials. Importantly, as the effect is seen on the protein level using IHC, the biomarker panel can be readily implemented in routine clinical testing.

Anticipating the clinical use of prognostic gene expression-based tests for colon cancer stage II and III: is Godot finally arriving?

PURPOSE: According to current recommendations for adjuvant treatment, patients with colon cancer stage II are not routinely offered chemotherapy, unless considered to have a high risk of relapse based on specific clinicopathological parameters. Following these criteria, it is challenging to identify the subgroup of patients that will benefit the most from adjuvant treatment. Contrarily, patients with colon cancer stage III are routinely offered chemotherapy, but due to expected adverse effects and frailty, elderly patients are often excluded from standard protocols. Colon cancer is a disease of the elderly and accordingly, there is a large subgroup of patients for which guidelines for adjuvant treatment remain less clear. In these two clinical settings, improved risk stratification has great potential impact on patient care, anticipating that high-risk patients will benefit from chemotherapy. However, microsatellite instability is the only molecular prognostic marker recommended for clinical use. EXPERIMENTAL DESIGN: In this perspective, we provide an updated view on the status and clinical potential of the many proposed prognostic gene expression-based tests for colon cancer stage II and III. RESULTS: The main limitation for clinical implementation is lack of prospective validation. For patients with stage II, highly promising tests have been identified and clinical trials are ongoing. For elderly patients with stage III, the value of such tests has received less focus, but promising early results have been shown. CONCLUSION: Although awaiting results from prospective trials, improved risk assessment for patients with stage II and III is likely to be achieved in the foreseeable future.

NEDD4 is overexpressed in colorectal cancer and promotes colonic cell growth independently of the PI3K/PTEN/AKT pathway.

Ubiquitination controls multiple cellular processes relevant to cancer pathogenesis. Using Gene Set Enrichment Analysis of an mRNA transcriptome dataset, we have identified genes encoding components of the ubiquitin system that are differentially expressed in colorectal cancers as compared to normal colonic mucosa. Among the significantly overexpressed genes was NEDD4 (neural precursor cell-expressed developmentally down-regulated 4), the prototype member of the HECT (homologous to E6AP C-terminus) E3 ubiquitin ligase family. Previous studies have shown that NEDD4 may act as an oncoprotein by inducing ubiquitination and degradation of the tumor suppressor protein PTEN (phosphatase and tensin homolog). To investigate its functional importance in colorectal cancer, HCT-15 and LoVo colon cancer cells were depleted of NEDD4 by small interfering RNA. The depletion resulted in reduced growth and altered cell morphology in both cell lines. However, NEDD4 depletion did not affect the PTEN protein level or PI3K/AKT signaling pathway activation. Moreover, ectopic expression of NEDD4 did not influence the PTEN subcellular localization or protein level. Collectively, these data demonstrate that NEDD4 is overexpressed in colorectal cancers, and suggest that NEDD4 promotes growth of colon cancer cells independently of PTEN and PI3K/AKT signaling.

Survival meta-analyses for >1800 malignant peripheral nerve sheath tumor patients with and without neurofibromatosis type 1.

There are conflicting reports as to whether malignant peripheral nerve sheath tumor (MPNST) patients with neurofibromatosis type 1 (NF1) have worse prognosis than non-NF1 MPNST patients. Large clinical studies to address this problem are lacking due to the rareness of MPNST. We have performed meta-analyses testing the effect of NF1 status on MPNST survival based on publications from the last 50 years, including only nonoverlapping patients reported from each institution. In addition, we analyzed survival characteristics for 179 MPNST patients from 3 European sarcoma centers. The meta-analyses including data from a total of 48 studies and >1800 patients revealed a significantly higher odds ratio for overall survival (OR(OS)) and disease-specific survival (OR(DSS)) in the non-NF1 group (OR(OS) = 1.75, 95% confidence interval [CI] = 1.28-2.39, and OR(DSS) = 1.68, 95% CI = 1.18-2.40). However, in studies published in the last decade, survival in the 2 patient groups has been converging, as especially the NF1 group has shown improved prognosis. For our own MPNST patients, NF1 status had no effect on overall or disease-specific survival. The compiled literature from 1963 to the present indicates a significantly worse outcome of MPNST in patients with NF1 syndrome compared with non-NF1 patients. However, survival for the NF1 patients has improved in the last decade, and the survival difference is diminishing. These observations support the hypothesis that MPNSTs arising in NF1 and non-NF1 patients are not different per se. Consequently, we suggest that the choice of treatment for MPNST should be independent of NF1 status.

MiR-9, -31, and -182 deregulation promote proliferation and tumor cell survival in colon cancer.

Several microRNAs (miRNAs) are known to be deregulated in colon cancer, but the mechanisms behind their potential involvement on proliferation and tumor cell survival are unclear. The present study aimed to identify miRNAs with functional implications for development of colon cancer. The cell proliferation and apoptosis were examined following perturbations of miRNA levels by employing a comprehensive miRNA library screen. miRNAs nominated for relevance to colon cancer were validated on expression and functional levels. By integrating the effect of miRNA up-regulation with the endogenous miRNA expression levels within the HT29, HCT116, and SW480 colon cancer cell lines, we identified miRNAs controlling cell proliferation (n = 53) and apoptosis (n = 93). From these functionally nominated miRNAs, we narrowed the list to 10 oncogene- and 20 tumor suppressor-like miRNAs that were also differentially expressed between colon cancer (n = 80) and normal colonic mucosa (n = 20). The differential expressions of miR-9, miR-31, and miR-182 were successfully validated in a series of colon carcinomas (n = 30) and polyps (n = 10) versus normal colonic mucosa (n = 10), whereas the functional effect was confirmed in an in-depth validation using different cell viability and apoptotic markers. Several transcription factors and genes regulating cell proliferation were identified as putative target genes by integrative miRNA/mRNA expression analysis obtained from the same colon cancer patient samples. This study suggests that deregulated expression of miR-9, miR-31, and miR-182 during carcinogenesis plays a significant role in the development of colon cancer by promoting proliferation and tumor cell survival.

ColoGuidePro: a prognostic 7-gene expression signature for stage III colorectal cancer patients.

PURPOSE: Improved prognostic stratification of patients with stage II and III colorectal cancer is warranted for postoperative clinical decision making. This study was conducted to develop a clinically feasible and robust prognostic classifier for these patients independent of adjuvant treatment. EXPERIMENTAL DESIGN: Global gene expression profiles from altogether 387 stage II and III colorectal cancer tissue samples from three independent patient series were included in the study. ColoGuidePro, a seven-gene prognostic classifier, was developed from a selected Norwegian learning series (n = 95; no adjuvant treatment) using lasso-penalized multivariate survival modeling with cross-validation. RESULTS: The expression signature significantly stratified patients in a consecutive Norwegian test series, in which patients were treated according to current standards [HR, 2.9 (1.1-7.5); P = 0.03; n = 77] and an external validation series [HR, 3.7 (2.0-6.8); P < 0.001; n = 215] according to survival. ColoGuidePro was also an independent predictor of prognosis in multivariate models including tumor stage in both series (HR, ≥ 3.1; P ≤ 0.03). In the validation series, which consisted of patients from other populations (United States and Australia), 5-year relapse-free survival was significantly predicted for stage III patients only (P < 0.001; n = 107). Here, prognostic stratification was independent of adjuvant treatment (P = 0.001). CONCLUSIONS: We present ColoGuidePro, a prognostic classifier developed for patients with stage II and III colorectal cancer. The test is suitable for transfer to clinical use and has best prognostic prediction potential for stage III patients.

ColoGuideEx: a robust gene classifier specific for stage II colorectal cancer prognosis.

BACKGROUND AND AIMS: Several clinical factors have an impact on prognosis in stage II colorectal cancer (CRC), but as yet they are inadequate for risk assessment. The present study aimed to develop a gene expression classifier for improved risk stratification of patients with stage II CRC. METHODS: 315 CRC samples were included in the study. Gene expression measurements from 207 CRC samples (stage I-IV) from two independent Norwegian clinical series were obtained using Affymetrix exon-level microarrays. Differentially expressed genes between stage I and stage IV samples from the test series were identified and used as input for L1 (lasso) penalised Cox proportional hazards analyses of patients with stage II CRC from the same series. A second validation was performed in 108 stage II CRC samples from other populations (USA and Australia). RESULTS: An optimal 13-gene expression classifier (PIGR, CXCL13, MMP3, TUBA1B, SESN1, AZGP1, KLK6, EPHA7, SEMA3A, DSC3, CXCL10, ENPP3, BNIP3) for prediction of relapse among patients with stage II CRC was developed using a consecutive Norwegian test series from patients treated according to current standard protocols (n=44, p<0.001, HR=18.2), and its predictive value was successfully validated for patients with stage II CRC in a second Norwegian CRC series collected two decades previously (n=52, p=0.02, HR=3.6). Further validation of the classifier was obtained in a recent external dataset of patients with stage II CRC from other populations (n=108, p=0.001, HR=6.5). Multivariate Cox regression analyses, including all three sample series and various clinicopathological variables, confirmed the independent prognostic value of the classifier (p≤0.004). The classifier was shown to be specific to stage II CRC and does not provide prognostic stratification of patients with stage III CRC. CONCLUSION: This study presents the development and validation of a 13-gene expression classifier, ColoGuideEx, for prognosis prediction specific to patients with stage II CRC. The robustness was shown across patient series, populations and different microarray versions.

The exon-level biomarker SLC39A14 has organ-confined cancer-specificity in colorectal cancer.

An alternative transcript variant of SLC39A14, caused by pre-mRNA splicing, was recently suggested as a biomarker for colorectal cancer (CRC). In our study, we have validated the cancer-specific splicing pattern of the mutually exclusive exons 4A and 4B in altogether 244 colorectal tissue samples. Exon-specific quantitative RT-PCR analyses across 136 Stage I-IV CRC samples and 44 normal colonic mucosa samples showed complete cancer-specificity, as well as 94% sensitivity of SLC39A14-exon4B relative to SLC39A14-exon4A expression. However, across 20 samples from a range of healthy tissues, 18 expressed the CRC variant. This was true also for ten benign lymph nodes, demonstrating that the cancer-specificity is mainly confined to the colon and rectum. Hence, clinical use of SLC39A14-exon4B as a detection marker for CRC other than in samples taken from the bowel wall is diminished. Prognostic value by detection of metastasis to lymph nodes is also abated, elucidating an important pitfall to biomarker discovery. However, analyses of ten nondysplastic biopsies from patients with active inflammatory bowel disease showed negative results in seven samples and only weakly positive results in three samples, suggesting value of SLC39A14-exon4B as a marker to distinguish CRC from other pathologic conditions of the colon. In conclusion, the SLC39A14-exon4B transcript variant is a CRC biomarker with high sensitivity and organ-confined specificity. Further use of the transcript and its encoded protein isoform should be explored in an organ-confined context.

Phospholipase C isozymes are deregulated in colorectal cancer--insights gained from gene set enrichment analysis of the transcriptome.

Colorectal cancer (CRC) is one of the most common cancer types in developed countries. To identify molecular networks and biological processes that are deregulated in CRC compared to normal colonic mucosa, we applied Gene Set Enrichment Analysis to two independent transcriptome datasets, including a total of 137 CRC and ten normal colonic mucosa samples. Eighty-two gene sets as described by the Kyoto Encyclopedia of Genes and Genomes database had significantly altered gene expression in both datasets. These included networks associated with cell division, DNA maintenance, and metabolism. Among signaling pathways with known changes in key genes, the "Phosphatidylinositol signaling network", comprising part of the PI3K pathway, was found deregulated. The downregulated genes in this pathway included several members of the Phospholipase C protein family, and the reduced expression of two of these, PLCD1 and PLCE1, were successfully validated in CRC biopsies (n = 70) and cell lines (n = 19) by quantitative analyses. The repression of both genes was found associated with KRAS mutations (P = 0.005 and 0.006, respectively), and we observed that microsatellite stable carcinomas with reduced PLCD1 expression more frequently had TP53 mutations (P = 0.002). Promoter methylation analyses of PLCD1 and PLCE1 performed in cell lines and tumor biopsies revealed that methylation of PLCD1 can contribute to reduced expression in 40% of the microsatellite instable carcinomas. In conclusion, we have identified significantly deregulated pathways in CRC, and validated repression of PLCD1 and PLCE1 expression. This illustrates that the GSEA approach may guide discovery of novel biomarkers in cancer.

Transcriptome instability in colorectal cancer identified by exon microarray analyses: Associations with splicing factor expression levels and patient survival.

BACKGROUND: Colorectal cancer (CRC) is a heterogeneous disease that, on the molecular level, can be characterized by inherent genomic instabilities; chromosome instability and microsatellite instability. In the present study we analyze genome-wide disruption of pre-mRNA splicing, and propose transcriptome instability as a characteristic that is analogous to genomic instability on the transcriptome level. METHODS: Exon microarray profiles from two independent series including a total of 160 CRCs were investigated for their relative amounts of exon usage differences. Each exon in each sample was assigned an alternative splicing score calculated by the FIRMA algorithm. Amounts of deviating exon usage per sample were derived from exons with extreme splicing scores. RESULTS: There was great heterogeneity within both series in terms of sample-wise amounts of deviating exon usage. This was strongly associated with the expression levels of approximately half of 280 splicing factors (54% and 48% of splicing factors were significantly correlated to deviating exon usage amounts in the two series). Samples with high or low amounts of deviating exon usage, associated with overall transcriptome instability, were almost completely separated into their respective groups by hierarchical clustering analysis of splicing factor expression levels in both sample series. Samples showing a preferential tendency towards deviating exon skipping or inclusion were associated with skewed transcriptome instability. There were significant associations between transcriptome instability and reduced patient survival in both sample series. In the test series, patients with skewed transcriptome instability showed the strongest prognostic association (P = 0.001), while a combination of the two characteristics showed the strongest association with poor survival in the validation series (P = 0.03). CONCLUSIONS: We have described transcriptome instability as a characteristic of CRC. This transcriptome instability has associations with splicing factor expression levels and poor patient survival.

DNA sequence profiles of the colorectal cancer critical gene set KRAS-BRAF-PIK3CA-PTEN-TP53 related to age at disease onset.

The incidence of colorectal cancer (CRC) increases with age and early onset indicates an increased likelihood for genetic predisposition for this disease. The somatic genetics of tumor development in relation to patient age remains mostly unknown. We have examined the mutation status of five known cancer critical genes in relation to age at diagnosis, and compared the genomic complexity of tumors from young patients without known CRC syndromes with those from elderly patients. Among 181 CRC patients, stratified by microsatellite instability status, DNA sequence changes were identified in KRAS (32%), BRAF (16%), PIK3CA (4%), PTEN (14%) and TP53 (51%). In patients younger than 50 years (n = 45), PIK3CA mutations were not observed and TP53 mutations were more frequent than in the older age groups. The total gene mutation index was lowest in tumors from the youngest patients. In contrast, the genome complexity, assessed as copy number aberrations, was highest in tumors from the youngest patients. A comparable number of tumors from young (<50 years) and old patients (>70 years) was quadruple negative for the four predictive gene markers (KRAS-BRAF-PIK3CA-PTEN); however, 16% of young versus only 1% of the old patients had tumor mutations in PTEN/PIK3CA exclusively. This implies that mutation testing for prediction of EGFR treatment response may be restricted to KRAS and BRAF in elderly (>70 years) patients. Distinct genetic differences found in tumors from young and elderly patients, whom are comparable for known clinical and pathological variables, indicate that young patients have a different genetic risk profile for CRC development than older patients.

Distinct high resolution genome profiles of early onset and late onset colorectal cancer integrated with gene expression data identify candidate susceptibility loci.

BACKGROUND: Estimates suggest that up to 30% of colorectal cancers (CRC) may develop due to an increased genetic risk. The mean age at diagnosis for CRC is about 70 years. Time of disease onset 20 years younger than the mean age is assumed to be indicative of genetic susceptibility. We have compared high resolution tumor genome copy number variation (CNV) (Roche NimbleGen, 385 000 oligo CGH array) in microsatellite stable (MSS) tumors from two age groups, including 23 young at onset patients without known hereditary syndromes and with a median age of 44 years (range: 28-53) and 17 elderly patients with median age 79 years (range: 69-87). Our aim was to identify differences in the tumor genomes between these groups and pinpoint potential susceptibility loci. Integration analysis of CNV and genome wide mRNA expression data, available for the same tumors, was performed to identify a restricted candidate gene list. RESULTS: The total fraction of the genome with aberrant copy number, the overall genomic profile and the TP53 mutation spectrum were similar between the two age groups. However, both the number of chromosomal aberrations and the number of breakpoints differed significantly between the groups. Gains of 2q35, 10q21.3-22.1, 10q22.3 and 19q13.2-13.31 and losses from 1p31.3, 1q21.1, 2q21.2, 4p16.1-q28.3, 10p11.1 and 19p12, positions that in total contain more than 500 genes, were found significantly more often in the early onset group as compared to the late onset group. Integration analysis revealed a covariation of DNA copy number at these sites and mRNA expression for 107 of the genes. Seven of these genes, CLC, EIF4E, LTBP4, PLA2G12A, PPAT, RG9MTD2, and ZNF574, had significantly different mRNA expression comparing median expression levels across the transcriptome between the two groups. CONCLUSIONS: Ten genomic loci, containing more than 500 protein coding genes, are identified as more often altered in tumors from early onset versus late onset CRC. Integration of genome and transcriptome data identifies seven novel candidate genes with the potential to identify an increased risk for CRC.

Genomic changes in chromosomes 10, 16, and X in malignant peripheral nerve sheath tumors identify a high-risk patient group.

PURPOSE: The purpose of this study was to identify genetic aberrations contributing to clinical aggressiveness of malignant peripheral nerve sheath tumors (MPNSTs). PATIENTS AND METHODS: Samples from 48 MPNSTs and 10 neurofibromas were collected from 51 patients with (n = 31) or without (n = 20) neurofibromatosis type 1 (NF1). Genome-wide DNA copy number changes were assessed by chromosomal and array-based comparative genomic hybridization (CGH) and examined for prognostic significance. For a subset of 20 samples, RNA microarray data were integrated with the genome data to identify potential target genes. RESULTS: Forty-four (92%) MPNSTs displayed DNA copy number changes (median, 18 changes per tumor; range, 2 to 35 changes). Known frequent chromosomal gains at chromosome arms 8q (69%), 17q (67%), and 7p (52%) and losses from 9p (50%), 11q (48%), and 17p (44%) were confirmed. Additionally, gains at 16p or losses from 10q or Xq identified a high-risk group with only 11% 10-year disease-specific survival (P = .00005). Multivariate analyses including NF1 status, tumor location, size, grade, sex, complete remission, and initial metastatic status showed that the genomic high-risk group was the most significant predictor of poor survival. Several genes whose expression was affected by the DNA copy number aberrations were identified. CONCLUSION: The presence of specific genetic aberrations was strongly associated with poor survival independent of known clinical risk factors. Conversely, within the total patient cohort with 34% 10-year disease-specific survival, a low-risk group was identified: without changes at chromosomes 10q, 16p, or Xq in their MPNSTs, the patients had 74% 10-year survival.