Performance of gene expression-based single sample predictors for assessment of clinicopathological subgroups and molecular subtypes in cancers: a case comparison study in non-small cell lung cancer.

The development of multigene classifiers for cancer prognosis, treatment prediction, molecular subtypes or clinicopathological groups has been a cornerstone in transcriptomic analyses of human malignancies for nearly two decades. However, many reported classifiers are critically limited by different preprocessing needs like normalization and data centering. In response, a new breed of classifiers, single sample predictors (SSPs), has emerged. SSPs classify samples in an N-of-1 fashion, relying on, e.g. gene rules comparing expression values within a sample. To date, several methods have been reported, but there is a lack of head-to-head performance comparison for typical cancer classification problems, representing an unmet methodological need in cancer bioinformatics. To resolve this need, we performed an evaluation of two SSPs [k-top-scoring pair classifier (kTSP) and absolute intrinsic molecular subtyping (AIMS)] for two case examples of different magnitude of difficulty in non-small cell lung cancer: gene expression-based classification of (i) tumor histology and (ii) molecular subtype. Through the analysis of ~2000 lung cancer samples for each case example (n = 1918 and n = 2106, respectively), we compared the performance of the methods for different sample compositions, training data set sizes, gene expression platforms and gene rule selections. Three main conclusions are drawn from the comparisons: both methods are platform independent, they select largely overlapping gene rules associated with actual underlying tumor biology and, for large training data sets, they behave interchangeably performance-wise. While SSPs like AIMS and kTSP offer new possibilities to move gene expression signatures/predictors closer to a clinical context, they are still importantly limited by the difficultness of the classification problem at hand.

Pre-operative plasma cell-free circulating tumor DNA and serum protein tumor markers as predictors of lung adenocarcinoma recurrence.

Background: Lung cancer patients have a risk of recurrence even after curatively intended surgery. Cell-free circulating tumor DNA (ctDNA) and circulating tumor marker measurements are easily accessible through peripheral blood and could potentially identify patients with worse prognosis. The aim of this study was to examine ctDNA in pre-operative plasma and the role of tumor markers in pre-operative serum for their predictive potential on risk of tumor recurrence. Methods: Mutation analysis by 26-gene targeted sequencing was performed on 157 lung adenocarcinomas (ACs) from patients surgically treated at the Lund University Hospital 2005-2014. Of these, 58 tumors from patients in stages I-IIIA (34 stage I, 14 stage II and 10 stage III) with mutation(s) in EGFR, BRAF or KRAS were included. ctDNA from corresponding plasma (median 1.5 ml, range 1-1.6) was analyzed for one tumor-specific mutation in either of these three oncogenes using ultrasensitive IBSAFE droplet digital PCR (ddPCR). The tumor markers cancer antigen 125 (CA 125) and carbohydrate antigen 19-9 (CA 19-9) were analyzed in corresponding serum with electrochemiluminiscence immunoassay. Results: 6/7 patients with ctDNA and 19/51 without detected ctDNA were diagnosed with recurrence (log-rank test p = .001). 8/10 patients with positive serum tumor markers and 17/47 without tumor markers were diagnosed with recurrence (log-rank test, p = .0002). Fifteen patients had positive ctDNA and/or tumor markers, 12 of these had recurrence (log-rank test, p < .0001). Conclusion: A combination of tumor markers and ctDNA single mutation detection in low-volume pre-operative blood samples is a promising prognostic test. Prediction of recurrent disease in surgically treated early stage lung cancer can likely be further improved by using larger volumes of blood.

A combined gene expression tool for parallel histological prediction and gene fusion detection in non-small cell lung cancer.

Accurate histological classification and identification of fusion genes represent two cornerstones of clinical diagnostics in non-small cell lung cancer (NSCLC). Here, we present a NanoString gene expression platform and a novel platform-independent, single sample predictor (SSP) of NSCLC histology for combined, simultaneous, histological classification and fusion gene detection in minimal formalin fixed paraffin embedded (FFPE) tissue. The SSP was developed in 68 NSCLC tumors of adenocarcinoma (AC), squamous cell carcinoma (SqCC) and large-cell neuroendocrine carcinoma (LCNEC) histology, based on NanoString expression of 11 (CHGA, SYP, CD56, SFTPG, NAPSA, TTF-1, TP73L, KRT6A, KRT5, KRT40, KRT16) relevant genes for IHC-based NSCLC histology classification. The SSP was combined with a gene fusion detection module (analyzing ALK, RET, ROS1, MET, NRG1, and NTRK1) into a multicomponent NanoString assay. The histological SSP was validated in six cohorts varying in size (n = 11-199), tissue origin (early or advanced disease), histological composition (including undifferentiated cancer), and gene expression platform. Fusion gene detection revealed five EML4-ALK fusions, four KIF5B-RET fusions, two CD74-NRG1 fusion and three MET exon 14 skipping events among 131 tested cases. The histological SSP was successfully trained and tested in the development cohort (mean AUC = 0.96 in iterated test sets). The SSP proved successful in predicting histology of NSCLC tumors of well-defined subgroups and difficult undifferentiated morphology irrespective of gene expression data platform. Discrepancies between gene expression prediction and histologic diagnosis included cases with mixed histologies, true large cell carcinomas, or poorly differentiated adenocarcinomas with mucin expression. In summary, we present a proof-of-concept multicomponent assay for parallel histological classification and multiplexed fusion gene detection in archival tissue, including a novel platform-independent histological SSP classifier. The assay and SSP could serve as a promising complement in the routine evaluation of diagnostic lung cancer biopsies.

NF1-mutated melanoma tumors harbor distinct clinical and biological characteristics.

In general, melanoma can be considered as a UV-driven disease with an aggressive metastatic course and high mutational load, with only few tumors (acral, mucosal, and uveal melanomas) not induced by sunlight and possessing a lower mutational load. The most commonly activated pathway in melanoma is the mitogen-activated protein kinase (MAPK) pathway. However, the prognostic significance of mutational stratification is unclear and needs further investigation. Here, in silico we combined mutation data from 162 melanomas subjected to targeted deep sequencing with mutation data from three published studies. Tumors from 870 patients were grouped according to BRAF, RAS, NF1 mutation or triple-wild-type status and correlated with tumor and patient characteristics. We found that the NF1-mutated subtype had a higher mutational burden and strongest UV mutation signature. Searching for co-occurring mutated genes revealed the RASopathy genes PTPN11 and RASA2, as well as another RAS domain-containing gene RASSF2 enriched in the NF1 subtype after adjustment for mutational burden. We found that a larger proportion of the NF1-mutant tumors were from males and with older age at diagnosis. Importantly, we found an increased risk of death from melanoma (disease-specific survival, DSS; HR, 1.9; 95% CI, 1.21-3.10; P = 0.046) and poor overall survival (OS; HR, 2.0; 95% CI, 1.28-2.98; P = 0.01) in the NF1 subtype, which remained significant after adjustment for age, gender, and lesion type (DSS P = 0.03, OS P = 0.06, respectively). Melanoma genomic subtypes display different biological and clinical characteristics. The poor outcome observed in the NF1 subtype highlights the need for improved characterization of this group.

Multiregion Whole-Exome Sequencing Uncovers the Genetic Evolution and Mutational Heterogeneity of Early-Stage Metastatic Melanoma.

Cancer genome sequencing has shed light on the underlying genetic aberrations that drive tumorigenesis. However, current sequencing-based strategies, which focus on a single tumor biopsy, fail to take into account intratumoral heterogeneity. To address this challenge and elucidate the evolutionary history of melanoma, we performed whole-exome and transcriptome sequencing of 41 multiple melanoma biopsies from eight individual tumors. This approach revealed heterogeneous somatic mutations in the range of 3%-38% in individual tumors. Known mutations in melanoma drivers BRAF and NRAS were always ubiquitous events. Using RNA sequencing, we found that the majority of mutations were not expressed or were expressed at very low levels, and preferential expression of a particular mutated allele did not occur frequently. In addition, we found that the proportion of ultraviolet B (UVB) radiation-induced C>T transitions differed significantly (P < 0.001) between early and late mutation acquisition, suggesting that different mutational processes operate during the evolution of metastatic melanoma. Finally, clinical history reports revealed that patients harboring a high degree of mutational heterogeneity were associated with more aggressive disease progression. In conclusion, our multiregion tumor-sequencing approach highlights the genetic evolution and non-UVB mutational signatures associated with melanoma development and progression, and may provide a more comprehensive perspective of patient outcome. Cancer Res; 76(16); 4765-74. ©2016 AACR.

Mutational and gene fusion analyses of primary large cell and large cell neuroendocrine lung cancer.

Large cell carcinoma with or without neuroendocrine features (LCNEC and LC, respectively) constitutes 3-9% of non-small cell lung cancer but is poorly characterized at the molecular level. Herein we analyzed 41 LC and 32 LCNEC (including 15 previously reported cases) tumors using massive parallel sequencing for mutations in 26 cancer-related genes and gene fusions in ALK, RET, and ROS1. LC patients were additionally subdivided into three immunohistochemistry groups based on positive expression of TTF-1/Napsin A (adenocarcinoma-like, n = 24; 59%), CK5/P40 (squamous-like, n = 5; 12%), or no marker expression (marker-negative, n = 12; 29%). Most common alterations were TP53 (83%), KRAS (22%), MET (12%) mutations in LCs, and TP53 (88%), STK11 (16%), and PTEN (13%) mutations in LCNECs. In general, LCs showed more oncogene mutations compared to LCNECs. Immunomarker stratification of LC revealed oncogene mutations in 63% of adenocarcinoma-like cases, but only in 17% of marker-negative cases. Moreover, marker-negative LCs were associated with inferior overall survival compared with adenocarcinoma-like tumors (p = 0.007). No ALK, RET or ROS1 fusions were detected in LCs or LCNECs. Together, our molecular analyses support that LC and LCNEC tumors follow different tumorigenic paths and that LC may be stratified into molecular subgroups with potential implications for diagnosis, prognostics, and therapy decisions.

Molecular stratification of metastatic melanoma using gene expression profiling: Prediction of survival outcome and benefit from molecular targeted therapy.

Melanoma is currently divided on a genetic level according to mutational status. However, this classification does not optimally predict prognosis. In prior studies, we have defined gene expression phenotypes (high-immune, pigmentation, proliferative and normal-like), which are predictive of survival outcome as well as informative of biology. Herein, we employed a population-based metastatic melanoma cohort and external cohorts to determine the prognostic and predictive significance of the gene expression phenotypes. We performed expression profiling on 214 cutaneous melanoma tumors and found an increased risk of developing distant metastases in the pigmentation (HR, 1.9; 95% CI, 1.05-3.28; P=0.03) and proliferative (HR, 2.8; 95% CI, 1.43-5.57; P=0.003) groups as compared to the high-immune response group. Further genetic characterization of melanomas using targeted deep-sequencing revealed similar mutational patterns across these phenotypes. We also used publicly available expression profiling data from melanoma patients treated with targeted or vaccine therapy in order to determine if our signatures predicted therapeutic response. In patients receiving targeted therapy, melanomas resistant to targeted therapy were enriched in the MITF-low proliferative subtype as compared to pre-treatment biopsies (P=0.02). In summary, the melanoma gene expression phenotypes are highly predictive of survival outcome and can further help to discriminate patients responding to targeted therapy.

A role for antizyme inhibitor in cell proliferation.

The polyamines are important for a variety of cellular functions, including cell growth. Their intracellular concentrations are controlled by a complex network of regulatory mechanisms, in which antizyme (Az) has a key role. Az reduces the cellular polyamine content by down-regulating both the enzyme catalysing polyamine biosynthesis, ornithine decarboxylase (ODC), and the uptake of polyamines. The activity of Az is repressed by the binding of a protein, named Az inhibitor (AzI), which is an enzymatically inactive homologue of ODC. Two forms of AzI have been described: AzI1, which is ubiquitous, and AzI2 which is expressed in brain and testis. In the present study, we have investigated the role of AzI1 in polyamine homeostasis and cell proliferation in breast cancer cells. The results obtained showed that the cellular content of AzI increased transiently after induction of cell proliferation by diluting cells in fresh medium. Inhibition of polyamine biosynthesis induced an even larger increase in the cellular AzI content, which remained significantly elevated during the 7-day experimental period. However, this increase was not a consequence of changes in cell cycle progression, as demonstrated by flow cytometry. Instead, the increase appeared to correlate with the cellular depletion of polyamines. Moreover, induced overexpression of AzI resulted in an increased cell proliferation with a concomitant increase in ODC activity and putrescine content. During mitosis, AzI1 was localised in a pattern that resembled that of the two centrosomes, confirming earlier observations. Taken together, the results indicate that AzI fulfils an essential regulatory function in polyamine homeostasis and cell proliferation.

Loss of CITED1, an MITF regulator, drives a phenotype switch in vitro and can predict clinical outcome in primary melanoma tumours.

CITED1 is a non-DNA binding transcriptional co-regulator whose expression can distinguish the 'proliferative' from 'invasive' signature in the phenotype-switching model of melanoma. We have found that, in addition to other 'proliferative' signature genes, CITED1 expression is repressed by TGFß while the 'invasive' signature genes are upregulated. In agreement, CITED1 positively correlates with MITF expression and can discriminate the MITF-high/pigmentation tumour molecular subtype in a large cohort (120) of melanoma cell lines. Interestingly, CITED1 overexpression significantly suppressed MITF promoter activation, mRNA and protein expression levels while MITF was transiently upregulated following siRNA mediated CITED1 silencing. Conversely, MITF siRNA silencing resulted in CITED1 downregulation indicating a reciprocal relationship. Whole genome expression analysis identified a phenotype shift induced by CITED1 silencing and driven mainly by expression of MITF and a cohort of MITF target genes that were significantly altered. Concomitantly, we found changes in the cell-cycle profile that manifest as transient G1 accumulation, increased expression of CDKN1A and a reduction in cell viability. Additionally, we could predict survival outcome by classifying primary melanoma tumours using our in vitro derived 'CITED1-silenced' gene expression signature. We hypothesize that CITED1 acts a regulator of MITF, functioning to maintain MITF levels in a range compatible with tumourigenesis.

Genome-Wide DNA Methylation Analysis in Melanoma Reveals the Importance of CpG Methylation in MITF Regulation.

The microphthalmia-associated transcription factor (MITF) is a key regulator of melanocyte development and a lineage-specific oncogene in melanoma; a highly lethal cancer known for its unpredictable clinical course. MITF is regulated by multiple intracellular signaling pathways, although the exact mechanisms that determine MITF expression and activity remain incompletely understood. In this study, we obtained genome-wide DNA methylation profiles from 50 stage IV melanomas, normal melanocytes, keratinocytes, and dermal fibroblasts and utilized The Cancer Genome Atlas data for experimental validation. By integrating DNA methylation and gene expression data, we found that hypermethylation of MITF and its co-regulated differentiation pathway genes corresponded to decreased gene expression levels. In cell lines with a hypermethylated MITF-pathway, overexpression of MITF did not alter the expression level or methylation status of the MITF pathway genes. In contrast, however, demethylation treatment of these cell lines induced MITF-pathway activity, confirming that gene regulation was controlled via methylation. The discovery that the activity of the master regulator of pigmentation, MITF, and its downstream targets may be regulated by hypermethylation has significant implications for understanding the development and evolvement of melanoma.

Molecular and genetic diversity in the metastatic process of melanoma.

Diversity between metastatic melanoma tumours in individual patients is known; however, the molecular and genetic differences remain unclear. To examine the molecular and genetic differences between metastatic tumours, we performed gene-expression profiling of 63 melanoma tumours obtained from 28 patients (two or three tumours/patient), followed by analysis of their mutational landscape, using targeted deep sequencing of 1697 cancer genes and DNA copy number analysis. Gene-expression signatures revealed discordant phenotypes between tumour lesions within a patient in 50% of the cases. In 18 of 22 patients (where matched normal tissue was available), we found that the multiple lesions within a patient were genetically divergent, with one or more melanoma tumours harbouring 'private' somatic mutations. In one case, the distant subcutaneous metastasis of one patient occurring 3 months after an earlier regional lymph node metastasis had acquired 37 new coding sequence mutations, including mutations in PTEN and CDH1. However, BRAF and NRAS mutations, when present in the first metastasis, were always preserved in subsequent metastases. The patterns of nucleotide substitutions found in this study indicate an influence of UV radiation but possibly also DNA alkylating agents. Our results clearly demonstrate that metastatic melanoma is a molecularly highly heterogeneous disease that continues to progress throughout its clinical course. The private aberrations observed on a background of shared aberrations within a patient provide evidence of continued evolution of individual tumours following divergence from a common parental clone, and might have implications for personalized medicine strategies in melanoma treatment.

The antiproliferative effect of dietary fiber phenolic compounds ferulic acid and p-coumaric acid on the cell cycle of Caco-2 cells.

Epidemiological and animal studies have shown that dietary fiber is protective against the development of colon cancer. Dietary fiber is a rich source of the hydroxycinnamic acids ferulic acid (FA) and p-coumaric acid (p-CA), which both may contribute to the protective effect. We have investigated the effects of FA and p-CA treatment on global gene expression in Caco-2 colon cancer cells. The Caco-2 cells were treated with 150 µM FA or p-CA for 24 h, and gene expression was analyzed with cDNA microarray technique. A total of 517 genes were significantly affected by FA and 901 by p-CA. As we previously have found that FA or p-CA treatment delayed cell cycle progression, we focused on genes involved in proliferation and cell cycle regulation. The expressions of a number of genes involved in centrosome assembly, such as RABGAP1 and CEP2, were upregulated by FA treatment as well as the gene for the S phase checkpoint protein SMC1L1. p-CA treatment upregulated CDKN1A expression and downregulated CCNA2, CCNB1, MYC, and ODC1. Some proteins corresponding to the affected genes were also studied. Taken together, the changes found can partly explain the effects of FA or p-CA treatment on cell cycle progression, specifically in the S phase by FA and G(2)/M phase by p-CA treatment.

Reduction of the putative CD44+CD24- breast cancer stem cell population by targeting the polyamine metabolic pathway with PG11047.

Cancer stem cells (CSCs) are considered to be of particular concern in cancer as they possess inherent properties of self-renewal and differentiation, along with expressing certain genes related to a mesenchymal phenotype. These features favour the promotion of tumour recurrence and metastasis in cancer patients. Thus, the optimal chemotherapeutic treatment should target the CSC population, either by killing these cells and/or by inducing their transition to a more differentiated epithelial-like phenotype. Experiments were carried out on the trastuzumab-resistant human epidermal growth factor receptor 2-overexpressing breast cancer cell line JIMT-1 to unravel the chemotherapeutic effects of the polyamine analogue [1N,12N]bis(ethyl)-cis-6,7-dehydrospermine (PG11047) and of the polyamine biosynthetic inhibitor 2-difluoromethylornithine (DFMO) on the CD44+CD24- CSC population. Furthermore, effects on the properties of self-renewal and epithelial/mesenchymal markers were also investigated. Treatment with PG11047 reduced the CD44+CD24- subpopulation of JIMT-1 cells by approximately 50%, inhibited and/or reduced self-renewal capability of the CSC population, decreased cell motility and induced expression of mesenchymal to epithelial transition-associated proteins that are involved in promoting an epithelial phenotype. By contrast, DFMO slightly increased the CD44+CD24- subpopulation, increased cell motility and the level of mesenchymal-related proteins. DFMO treatment reduced the self-renewal capability of the CSC population. Both PG11047 and DFMO reduced the expression of the human epidermal growth factor receptor 2 protein, which is correlated to malignancy and resistance to trastuzumab in JIMT-1 cells. Our findings indicate that treatment with PG11047 targeted the CSC population by interfering with several stem cell-related properties, such as self-renewal, differentiation, motility and the mesenchymal phenotype.