Consensus of gene expression phenotypes and prognostic risk predictors in primary lung adenocarcinoma

Proteogenomic analysis of lung adenocarcinoma reveals tumor heterogeneity, survival determinants, and therapeutically relevant pathways

We present a deep proteogenomic profiling study of 87 lung adenocarcinoma (LUAD) tumors from the United States, integrating whole-genome sequencing, transcriptome sequencing, proteomics and phosphoproteomics by mass spectrometry, and reverse-phase protein arrays. We identify three subtypes from somatic genome signature analysis, including a transition-high subtype enriched with never smokers, a transversion-high subtype enriched with current smokers, and a structurally altered subtype enriched with former smokers, TP53 alterations, and genome-wide structural alterations. We show that within-tumor correlations of RNA and protein expression associate with tumor purity and immune cell profiles. We detect and independently validate expression signatures of RNA and protein that predict patient survival. Additionally, among co-measured genes, we found that protein expression is more often associated with patient survival than RNA. Finally, integrative analysis characterizes three expression subtypes with divergent mutations, proteomic regulatory networks, and therapeutic vulnerabilities. This proteogenomic characterization provides a foundation for molecularly informed medicine in LUAD.

SRIQ clustering: A fusion of Random Forest, QT clustering, and KNN concepts

Gene expression profiling together with unsupervised analysis methods, typically clustering methods, has been used extensively in cancer research to unravel, e.g., new molecular subtypes that hold promise of disease refinement that may ultimately benefit patients. However, many of the commonly used methods require a prespecified number of clusters to extract and frequently require some type of feature pre-selection, e.g. variance filtering. This introduces subjectivity to the process of cluster discovery and the definition of putative novel tumor subtypes. Here, we introduce SRIQ, a novel unsupervised clustering method that could circumvent some of the issues in commonly used unsupervised analysis methods. SRIQ incorporates concepts from random forest machine learning as well as quality threshold- and k-nearest neighbor clustering. It is implemented as a Java and Python pipeline including data pre-processing, differential expression analysis, and pathway analysis. Using 434 lung adenocarcinomas profiled by RNA sequencing, we demonstrate the technical reproducibility of SRIQ and benchmark its performance compared to the commonly used consensus clustering method. Based on differential gene expression analysis and auxiliary molecular data we show that SRIQ can define new tumor subsets that appear biologically relevant and consistent compared and that these new subgroups seem to refine existing transcriptional subtypes that were defined using consensus clustering. Together, this provides support that SRIQ may be a useful new tool for unsupervised analysis of gene expression data from human malignancies.

Gene Expression Profiling as a Potential Tool for Precision Oncology in Non-Small Cell Lung Cancer

Recent technological advances and the application of high-throughput mutation and transcriptome analyses have improved our understanding of cancer diseases, including non-small cell lung cancer. For instance, genomic profiling has allowed the identification of mutational events which can be treated with specific agents. However, detection of DNA alterations does not fully recapitulate the complexity of the disease and it does not allow selection of patients that benefit from chemo- or immunotherapy. In this context, transcriptional profiling has emerged as a promising tool for patient stratification and treatment guidance. For instance, transcriptional profiling has proven to be especially useful in the context of acquired resistance to targeted therapies and patients lacking targetable genomic alterations. Moreover, the comprehensive characterization of the expression level of the different pathways and genes involved in tumor progression is likely to better predict clinical benefit from different treatments than single biomarkers such as PD-L1 or tumor mutational burden in the case of immunotherapy. However, intrinsic technical and analytical limitations have hindered the use of these expression signatures in the clinical setting. In this review, we will focus on the data reported on molecular classification of non-small cell lung cancer and discuss the potential of transcriptional profiling as a predictor of survival and as a patient stratification tool to further personalize treatments.

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.

Pan-cancer application of a lung-adenocarcinoma-derived gene-expression-based prognostic predictor

Gene-expression profiling can be used to classify human tumors into molecular subtypes or risk groups, representing potential future clinical tools for treatment prediction and prognostication. However, it is less well-known how prognostic gene signatures derived in one malignancy perform in a pan-cancer context. In this study, a gene-rule-based single sample predictor (SSP) called classifier for lung adenocarcinoma molecular subtypes (CLAMS) associated with proliferation was tested in almost 15 000 samples from 32 cancer types to classify samples into better or worse prognosis. Of the 14 malignancies that presented both CLAMS classes in sufficient numbers, survival outcomes were significantly different for breast, brain, kidney and liver cancer. Patients with samples classified as better prognosis by CLAMS were generally of lower tumor grade and disease stage, and had improved prognosis according to other type-specific classifications (e.g. PAM50 for breast cancer). In all, 99.1% of non-lung cancer cases classified as better outcome by CLAMS were comprised within the range of proliferation scores of lung adenocarcinoma cases with a predicted better prognosis by CLAMS. This finding demonstrates the potential of tuning SSPs to identify specific levels of for instance tumor proliferation or other transcriptional programs through predictor training. Together, pan-cancer studies such as this may take us one step closer to understanding how gene-expression-based SSPs act, which gene-expression programs might be important in different malignancies, and how to derive tools useful for prognostication that are efficient across organs.

A gene expression-based single sample predictor of lung adenocarcinoma molecular subtype and prognosis

Disease recurrence in surgically treated lung adenocarcinoma (AC) remains high. New approaches for risk stratification beyond tumor stage are needed. Gene expression-based AC subtypes such as the Cancer Genome Atlas Network (TCGA) terminal-respiratory unit (TRU), proximal-inflammatory (PI) and proximal-proliferative (PP) subtypes have been associated with prognosis, but show methodological limitations for robust clinical use. We aimed to derive a platform independent single sample predictor (SSP) for molecular subtype assignment and risk stratification that could function in a clinical setting. Two-class (TRU/nonTRU=SSP2) and three-class (TRU/PP/PI=SSP3) SSPs using the AIMS algorithm were trained in 1655 ACs (n = 9659 genes) from public repositories vs TCGA centroid subtypes. Validation and survival analysis were performed in 977 patients using overall survival (OS) and distant metastasis-free survival (DMFS) as endpoints. In the validation cohort, SSP2 and SSP3 showed accuracies of 0.85 and 0.81, respectively. SSPs captured relevant biology previously associated with the TCGA subtypes and were associated with prognosis. In survival analysis, OS and DMFS for cases discordantly classified between TCGA and SSP2 favored the SSP2 classification. In resected Stage I patients, SSP2 identified TRU-cases with better OS (hazard ratio [HR] = 0.30; 95% confidence interval [CI] = 0.18-0.49) and DMFS (TRU HR = 0.52; 95% CI = 0.33-0.83) independent of age, Stage IA/IB and gender. SSP2 was transformed into a NanoString nCounter assay and tested in 44 Stage I patients using RNA from formalin-fixed tissue, providing prognostic stratification (relapse-free interval, HR = 3.2; 95% CI = 1.2-8.8). In conclusion, gene expression-based SSPs can provide molecular subtype and independent prognostic information in early-stage lung ACs. SSPs may overcome critical limitations in the applicability of gene signatures in lung cancer.

Simplified molecular classification of lung adenocarcinomas based on EGFR, KRAS, and TP53 mutations

Background

Gene expression profiling has consistently identified three molecular subtypes of lung adenocarcinoma that have prognostic implications. To facilitate stratification of patients with this disease into similar molecular subtypes, we developed and validated a simple, mutually exclusive classification.

Methods

Mutational status of EGFR, KRAS, and TP53 was used to define seven mutually exclusive molecular subtypes. A development cohort of 283 cytology specimens of lung adenocarcinoma was used to evaluate the associations between the proposed classification and clinicopathologic variables including demographic characteristics, smoking history, fluorescence in situ hybridization and molecular results. For validation and prognostic assessment, 63 of the 283 cytology specimens with available survival data were combined with a separate cohort of 428 surgical pathology specimens of lung adenocarcinoma.

Results

The proposed classification yielded significant associations between these molecular subtypes and clinical and prognostic features. We found better overall survival in patients who underwent surgery and had tumors enriched for EGFR mutations. Worse overall survival was associated with older age, stage IV disease, and tumors with co-mutations in KRAS and TP53. Interestingly, neither chemotherapy nor radiation therapy showed benefit to overall survival.

Conclusions

The mutational status of EGFR, KRAS, and TP53 can be used to easily classify lung adenocarcinoma patients into seven subtypes that show a relationship with prognosis, especially in patients who underwent surgery, and these subtypes are similar to classifications based on more complex genomic methods reported previously.

A Meta-Analysis of the Performance of a Blood-Based Exposure Response Gene Signature Across Clinical Studies on the Tobacco Heating System 2.2 (THS 2.2)

As part of emerging tobacco harm reduction strategies, modified risk tobacco products (MRTP) are being developed to offer alternatives that have the potential to reduce the individual risk and population harm compared with smoking cigarettes for adult smokers who want to continue using tobacco and nicotine products. MRTPs are defined as any tobacco products that are distributed for use to reduce harm or the risk of tobacco-related disease associated with commercially marketed tobacco products. One such candidate MRTP is the Tobacco Heating System (THS) 2.2, which does not burn tobacco but instead heats it, thus producing significantly reduced levels of harmful and potentially harmful constituents compared with cigarettes. The clinical assessment of candidate MRTPs requires the development of exposure-response markers to distinguish current smokers from either nonsmokers or former smokers with high specificity and sensitivity. Toward this end, a whole blood-derived gene signature was previously developed and reported. Four randomized, controlled, open-label, three-arm parallel group reduced exposure clinical studies have been conducted with subjects randomized to three arms: switching from cigarettes to THS 2.2, continuous use of cigarettes, or smoking abstinence. These clinical studies had an investigational period of 5 days in confinement, which was followed by an 85-day ambulatory period in two studies. Here we tested the previously developed blood-derived signature on the samples derived from those clinical studies. We showed that in all four studies, the signature scores were reduced consistently in subjects who either stopped smoking or switched to THS 2.2 compared with subjects who continued smoking cigarettes.

Prognostic value of cancer antigen -125 for lung adenocarcinoma patients with brain metastasis: A random survival forest prognostic model

Using random survival forest, this study was intended to evaluate the prognostic value of serum markers for lung adenocarcinoma patients with brain metastasis (BM), and tried to integrate them into a prognostic model. During 2010 to 2015, the patients were retrieved from two medical centers. Besides the Cox proportional hazards regression, the random survival forest (RSF) were also used to develop prognostic model from the group A (n = 142). In RSF of the group A, the factors, whose minimal depth were greater than the depth threshold or had a negative variable importance (VIMP), were firstly excluded. Subsequently, C-index and Akaike information criterion (AIC) were used to guide us finding models with higher prognostic ability and lower overfitting possibility. These RSF models, together with the Cox, modified-RPA and lung-GPA index were validated and compared, especially in the group B (CAMS, n = 53). Our data indicated that the KSE125 model (KPS, smoking, EGFR-20 (exon 18, 19 and 21) and Ca125) was the best in survival prediction, and performed well in internal and external validation. In conclusions, for lung adenocarcinoma patients with brain metastasis, a validated prognostic nomogram (KPS, smoking, EGFR-20 and Ca125) can more accurately predict 1-year and 2-year survival of the patients.

Gene Expression Profiling of Large Cell Lung Cancer Links Transcriptional Phenotypes to the New Histological WHO 2015 Classification

INTRODUCTION

Large cell lung cancer (LCLC) and large cell neuroendocrine carcinoma (LCNEC) constitute a small proportion of NSCLC. The WHO 2015 classification guidelines changed the definition of the debated histological subtype LCLC to be based on immunomarkers for adenocarcinoma and squamous cancer. We sought to determine whether these new guidelines also translate into the transcriptional landscape of lung cancer, and LCLC specifically.

METHODS

Gene expression profiling was performed by using Illumina V4 HT12 microarrays (Illumina, San Diego, CA) on samples from 159 cases (comprising all histological subtypes, including 10 classified as LCLC WHO 2015 and 14 classified as LCNEC according to the WHO 2015 guidelines), with complimentary mutational and immunohistochemical data. Derived transcriptional phenotypes were validated in 199 independent tumors, including six WHO 2015 LCLCs and five LCNECs.

RESULTS

Unsupervised analysis of gene expression data identified a phenotype comprising 90% of WHO 2015 LCLC tumors, with characteristics of poorly differentiated proliferative cancer, a 90% tumor protein p53 gene (TP53) mutation rate, and lack of well-known NSCLC oncogene driver alterations. Validation in independent data confirmed aggregation of WHO 2015 LCLCs in the specific phenotype. For LCNEC tumors, the unsupervised gene expression analysis suggested two different transcriptional patterns corresponding to a proposed genetic division of LCNEC tumors into SCLC-like and NSCLC-like cancer on the basis of TP53 and retinoblastoma 1 gene (RB1) alteration patterns.

CONCLUSIONS

Refined classification of LCLC has implications for diagnosis, prognostics, and therapy decisions. Our molecular analyses support the WHO 2015 classification of LCLC and LCNEC tumors, which herein follow different tumorigenic paths and can accordingly be stratified into different transcriptional subgroups, thus linking diagnostic immunohistochemical staining-driven classification with the transcriptional landscape of lung cancer.