The prevalence of early-stage non-small cell lung cancer (NSCLC) is expected to increase with recent implementation of annual screening programs. Reliable prognostic biomarkers are needed to identify ...patients at a high risk for recurrence to guide adjuvant therapy.
To develop a robust, individualized immune signature that can estimate prognosis in patients with early-stage nonsquamous NSCLC.
This retrospective study analyzed the gene expression profiles of frozen tumor tissue samples from 19 public NSCLC cohorts, including 18 microarray data sets and 1 RNA-Seq data set for The Cancer Genome Atlas (TCGA) lung adenocarcinoma cohort. Only patients with nonsquamous NSCLC with clinical annotation were included. Samples were from 2414 patients with nonsquamous NSCLC, divided into a meta-training cohort (729 patients), meta-testing cohort (716 patients), and 3 independent validation cohorts (439, 323, and 207 patients). All patients underwent surgery with a negative surgical margin, received no adjuvant or neoadjuvant therapy, and had publicly available gene expression data and survival information. Data were collected from July 22 through September 8, 2016.
Overall survival.
Of 2414 patients (1205 men 50%, 1111 women 46%, and 98 of unknown sex 4%; median age range, 64 15-90 years), a prognostic immune signature of 25 gene pairs consisting of 40 unique genes was constructed using the meta-training data set. In the meta-testing and validation cohorts, the immune signature significantly stratified patients into high- vs low-risk groups in terms of overall survival across and within subpopulations with stage I, IA, IB, or II disease and remained as an independent prognostic factor in multivariate analyses (hazard ratio range, 1.72 95% CI, 1.26-2.33; P < .001 to 2.36 95% CI, 1.47-3.79; P < .001) after adjusting for clinical and pathologic factors. Several biological processes, including chemotaxis, were enriched among genes in the immune signature. The percentage of neutrophil infiltration (5.6% vs 1.8%) and necrosis (4.6% vs 1.5%) was significantly higher in the high-risk immune group compared with the low-risk groups in TCGA data set (P < .003). The immune signature achieved a higher accuracy (mean concordance index C-index, 0.64) than 2 commercialized multigene signatures (mean C-index, 0.53 and 0.61) for estimation of survival in comparable validation cohorts. When integrated with clinical characteristics such as age and stage, the composite clinical and immune signature showed improved prognostic accuracy in all validation data sets relative to molecular signatures alone (mean C-index, 0.70 vs 0.63) and another commercialized clinical-molecular signature (mean C-index, 0.68 vs 0.65).
The proposed clinical-immune signature is a promising biomarker for estimating overall survival in nonsquamous NSCLC, including early-stage disease. Prospective studies are needed to test the clinical utility of the biomarker in individualized management of nonsquamous NSCLC.
Single-cell RNA-sequencing has emerged as a powerful technique for characterizing cellular heterogeneity, but it is currently impractical on large sample cohorts and cannot be applied to fixed ...specimens collected as part of routine clinical care. We previously developed an approach for digital cytometry, called CIBERSORT, that enables estimation of cell type abundances from bulk tissue transcriptomes. We now introduce CIBERSORTx, a machine learning method that extends this framework to infer cell-type-specific gene expression profiles without physical cell isolation. By minimizing platform-specific variation, CIBERSORTx also allows the use of single-cell RNA-sequencing data for large-scale tissue dissection. We evaluated the utility of CIBERSORTx in multiple tumor types, including melanoma, where single-cell reference profiles were used to dissect bulk clinical specimens, revealing cell-type-specific phenotypic states linked to distinct driver mutations and response to immune checkpoint blockade. We anticipate that digital cytometry will augment single-cell profiling efforts, enabling cost-effective, high-throughput tissue characterization without the need for antibodies, disaggregation or viable cells.
We introduce CIBERSORT, a method for characterizing cell composition of complex tissues from their gene expression profiles. When applied to enumeration of hematopoietic subsets in RNA mixtures from ...fresh, frozen and fixed tissues, including solid tumors, CIBERSORT outperformed other methods with respect to noise, unknown mixture content and closely related cell types. CIBERSORT should enable large-scale analysis of RNA mixtures for cellular biomarkers and therapeutic targets (http://cibersort.stanford.edu/).
Growing evidence demonstrates that circulating tumor DNA (ctDNA) minimal residual disease (MRD) following treatment for solid tumors predicts relapse. These results suggest that ctDNA MRD could ...identify candidates for adjuvant therapy and measure response to such treatment. Importantly, factors such as assay type, amount of ctDNA release, and technical and biological background can affect ctDNA MRD results. Furthermore, the clinical utility of ctDNA MRD for treatment personalization remains to be fully established. Here, we review the evidence supporting the value of ctDNA MRD in solid cancers and highlight key considerations in the application of this potentially transformative biomarker.
ctDNA analysis enables detection of MRD and predicts relapse after definitive treatment for solid cancers, thereby promising to revolutionize personalization of adjuvant and consolidation therapies.
Molecular profiles of tumors and tumor-associated cells hold great promise as biomarkers of clinical outcomes. However, existing data sets are fragmented and difficult to analyze systematically. Here ...we present a pan-cancer resource and meta-analysis of expression signatures from ∼18,000 human tumors with overall survival outcomes across 39 malignancies. By using this resource, we identified a forkhead box MI (FOXM1) regulatory network as a major predictor of adverse outcomes, and we found that expression of favorably prognostic genes, including KLRB1 (encoding CD161), largely reflect tumor-associated leukocytes. By applying CIBERSORT, a computational approach for inferring leukocyte representation in bulk tumor transcriptomes, we identified complex associations between 22 distinct leukocyte subsets and cancer survival. For example, tumor-associated neutrophil and plasma cell signatures emerged as significant but opposite predictors of survival for diverse solid tumors, including breast and lung adenocarcinomas. This resource and associated analytical tools (http://precog.stanford.edu) may help delineate prognostic genes and leukocyte subsets within and across cancers, shed light on the impact of tumor heterogeneity on cancer outcomes, and facilitate the discovery of biomarkers and therapeutic targets.
Metabolic reprogramming of the tumor microenvironment is recognized as a cancer hallmark. To identify new molecular processes associated with tumor metabolism, we analyzed the transcriptome of bulk ...and flow-sorted human primary non-small cell lung cancer (NSCLC) together with
FDG-PET scans, which provide a clinical measure of glucose uptake. Tumors with higher glucose uptake were functionally enriched for molecular processes associated with invasion in adenocarcinoma and cell growth in squamous cell carcinoma (SCC). Next, we identified genes correlated to glucose uptake that were predominately overexpressed in a single cell-type comprising the tumor microenvironment. For SCC, most of these genes were expressed by malignant cells, whereas in adenocarcinoma, they were predominately expressed by stromal cells, particularly cancer-associated fibroblasts (CAF). Among these adenocarcinoma genes correlated to glucose uptake, we focused on glutamine-fructose-6-phosphate transaminase 2 (
), which codes for the glutamine-fructose-6-phosphate aminotransferase 2 (GFAT2), a rate-limiting enzyme of the hexosamine biosynthesis pathway (HBP), which is responsible for glycosylation.
was predictive of glucose uptake independent of GLUT1, the primary glucose transporter, and was prognostically significant at both gene and protein level. We confirmed that normal fibroblasts transformed to CAF-like cells, following TGFβ treatment, upregulated HBP genes, including
, with less change in genes driving glycolysis, pentose phosphate pathway, and TCA cycle. Our work provides new evidence of histology-specific tumor stromal properties associated with glucose uptake in NSCLC and identifies
as a critical regulator of tumor metabolic reprogramming in adenocarcinoma.
These findings implicate the hexosamine biosynthesis pathway as a potential new therapeutic target in lung adenocarcinoma.
.
Because cancer is caused by an accumulation of genetic mutations, mutant DNA released by tumors can be used as a highly specific biomarker for cancer. Although this principle was described decades ...ago, the advent and falling costs of next-generation sequencing have made the use of tumor DNA as a biomarker increasingly practical. This review surveys the use of cellular and cell-free DNA for the detection of cancer, with a focus on recent technological developments and applications to solid tumors. It covers (
) key principles and technology enabling the highly sensitive detection of tumor DNA; (
) assessment of tumor DNA in plasma, including for genotyping, minimal residual disease detection, and early detection of localized cancer; (
) detection of tumor DNA in body cavity fluids, such as urine or cerebrospinal fluid; and (
) challenges posed to the use of tumor DNA as a biomarker by the phenomenon of benign clonal expansions.
Accurate prediction of antigen presentation by human leukocyte antigen (HLA) class II molecules would be valuable for vaccine development and cancer immunotherapies. Current computational methods ...trained on in vitro binding data are limited by insufficient training data and algorithmic constraints. Here we describe MARIA (major histocompatibility complex analysis with recurrent integrated architecture; https://maria.stanford.edu/ ), a multimodal recurrent neural network for predicting the likelihood of antigen presentation from a gene of interest in the context of specific HLA class II alleles. In addition to in vitro binding measurements, MARIA is trained on peptide HLA ligand sequences identified by mass spectrometry, expression levels of antigen genes and protease cleavage signatures. Because it leverages these diverse training data and our improved machine learning framework, MARIA (area under the curve = 0.89-0.92) outperformed existing methods in validation datasets. Across independent cancer neoantigen studies, peptides with high MARIA scores are more likely to elicit strong CD4
T cell responses. MARIA allows identification of immunogenic epitopes in diverse cancers and autoimmune disease.
Circulating tumor-derived DNA (ctDNA) is an emerging biomarker for many cancers, but the limited sensitivity of current detection methods reduces its utility for diagnosing minimal residual disease. ...Here we describe phased variant enrichment and detection sequencing (PhasED-seq), a method that uses multiple somatic mutations in individual DNA fragments to improve the sensitivity of ctDNA detection. Leveraging whole-genome sequences from 2,538 tumors, we identify phased variants and their associations with mutational signatures. We show that even without molecular barcodes, the limits of detection of PhasED-seq outperform prior methods, including duplex barcoding, allowing ctDNA detection in the ppm range in participant samples. We profiled 678 specimens from 213 participants with B cell lymphomas, including serial cell-free DNA samples before and during therapy for diffuse large B cell lymphoma. In participants with undetectable ctDNA after two cycles of therapy using a next-generation sequencing-based approach termed cancer personalized profiling by deep sequencing, an additional 25% have ctDNA detectable by PhasED-seq and have worse outcomes. Finally, we demonstrate the application of PhasED-seq to solid tumors.
Outcomes for patients with diffuse large B-cell lymphoma remain heterogeneous, with existing methods failing to consistently predict treatment failure. We examined the additional prognostic value of ...circulating tumor DNA (ctDNA) before and during therapy for predicting patient outcomes.
We studied the dynamics of ctDNA from 217 patients treated at six centers, using a training and validation framework. We densely characterized early ctDNA dynamics during therapy using cancer personalized profiling by deep sequencing to define response-associated thresholds within a discovery set. These thresholds were assessed in two independent validation sets. Finally, we assessed the prognostic value of ctDNA in the context of established risk factors, including the International Prognostic Index and interim positron emission tomography/computed tomography scans.
Before therapy, ctDNA was detectable in 98% of patients; pretreatment levels were prognostic in both front-line and salvage settings. In the discovery set, ctDNA levels changed rapidly, with a 2-log decrease after one cycle (early molecular response EMR) and a 2.5-log decrease after two cycles (major molecular response MMR) stratifying outcomes. In the first validation set, patients receiving front-line therapy achieving EMR or MMR had superior outcomes at 24 months (EMR: EFS, 83% v 50%; P = .0015; MMR: EFS, 82% v 46%; P < .001). EMR also predicted superior 24-month outcomes in patients receiving salvage therapy in the first validation set (EFS, 100% v 13%; P = .011). The prognostic value of EMR and MMR was further confirmed in the second validation set. In multivariable analyses including International Prognostic Index and interim positron emission tomography/computed tomography scans across both cohorts, molecular response was independently prognostic of outcomes, including event-free and overall survival.
Pretreatment ctDNA levels and molecular responses are independently prognostic of outcomes in aggressive lymphomas. These risk factors could potentially guide future personalized risk-directed approaches.