Lung adenocarcinomas (LUADs) lead to the majority of deaths attributable to lung cancer. We performed whole-exome sequencing (WES) and immune profiling analyses of a unique set of clinically annotated early-stage LUADs to better understand the pathogenesis of this disease and identify clinically relevant molecular markers. We performed WES of 108 paired stage I-III LUADs and normal lung tissues using the Illumina HiSeq 2000 platform. Ten immune markers (PD-L1, PD-1, CD3, CD4, CD8, CD45ro, CD57, CD68, FOXP3 and Granzyme B) were profiled by imaging-based immunohistochemistry (IHC) in a subset of LUADs (n = 92). Associations among mutations, immune markers and clinicopathological variables were analyzed using ANOVA and Fisher’s exact test. Cox proportional hazards regression models were used for multivariate analysis of clinical outcome. LUADs in this cohort exhibited an average of 243 coding mutations. We identified 28 genes with significant enrichment for mutation. SETD2-mutated LUADs exhibited relatively poor recurrence- free survival (RFS) and mutations in STK11 and ATM were associated with poor RFS among KRAS-mutant tumors. EGFR, KEAP1 and PIK3CA mutations were predictive of poor response to adjuvant therapy. Immune marker analysis revealed that LUADs in smokers and with relatively high mutation burdens exhibited increased levels of immune markers. Analysis of immunophenotypes revealed that LUADs with STK11 mutations exhibited relatively low levels of infiltrating CD4+/CD8+ T-cells indicative of a muted immune response. Tumoral PD-L1 was significantly elevated in TP53 mutant LUADs whereas PIK3CA mutant LUADs exhibited markedly down-regulated PD-L1 expression. LUADs with TP53 or KEAP1 mutations displayed relatively increased CD57 and Granzyme B levels indicative of augmented natural killer (NK) cell infiltration. Our study highlights molecular and immune phenotypes that warrant further analysis for their roles in clinical outcomes and personalized immune-based therapy of LUAD.

Significance: The knowledge that cancer is an evolutionary process is old, but only recently can sequencing technology provide data for clinically relevant evolutionary analyses of cancer. Approaches developed for evolutionary biology can reveal the relationship among clonal lineages, the ancestral states of gene sequences, and the timing of evolutionary events. We performed whole exome sequencing of cancer tissues from multiple sites of dozens of subjects, demonstrating nonlinear patterns of tumor progression and early origins of metastatic lineages and quantifying the times of occurrence of driver mutations. These findings direct research attention away from the search for genes that induce metastasis toward genes that are mutated early in tumorigenesis, providing therapeutic targets effective against both primary tumors and metastases.

Lung squamous cell carcinoma (LUSC) accounts for 20–30% of non-small cell lung cancers (NSCLCs). There are limited treatment strategies for LUSC in part due to our inadequate understanding of the molecular underpinnings of the disease. We performed whole-exome sequencing (WES) and comprehensive immune profiling of a unique set of clinically annotated early-stage LUSCs to increase our understanding of the pathobiology of this malignancy. Matched pairs of surgically resected stage I-III LUSCs and normal lung tissues (n = 108) were analyzed by WES. Immunohistochemistry and image analysis-based profiling of 10 immune markers were done on a subset of LUSCs (n = 91). Associations among mutations, immune markers and clinicopathological variables were statistically examined using analysis of variance and Fisher’s exact test. Cox proportional hazards regression models were used for statistical analysis of clinical outcome. This early-stage LUSC cohort displayed an average of 209 exonic mutations per tumor. Fourteen genes exhibited significant enrichment for somatic mutation: TP53, MLL2, PIK3CA, NFE2L2, CDH8, KEAP1, PTEN, ADCY8, PTPRT, CALCR, GRM8, FBXW7, RB1 and CDKN2A. Among mutated genes associated with poor recurrence-free survival, MLL2 mutations predicted poor prognosis in both TP53 mutant and wild-type LUSCs. We also found that in treated patients, FBXW7 and KEAP1 mutations were associated with poor response to adjuvant therapy, particularly in TP53-mutant tumors. Analysis of mutations with immune markers revealed that ADCY8 and PIK3CA mutations were associated with markedly decreased tumoral PD-L1 expression, LUSCs with PIK3CA mutations exhibited elevated CD45ro levels and CDKN2A-mutant tumors displayed an up-regulated immune response. Our findings pinpoint mutated genes that may impact clinical outcome as well as personalized strategies for targeted immunotherapies in early-stage LUSC.

Recent studies have revealed the mutational signatures underlying the somatic evolution of cancer, and the prevalences of associated somatic genetic variants. Here we estimate the intensity of positive selection that drives mutations to high frequency in tumors, yielding higher prevalences than expected on the basis of mutation and neutral drift alone. We apply this approach to a sample of 525 head and neck squamous cell carcinoma exomes, producing a rank-ordered list of gene variants by selection intensity. Our results illustrate the complementarity of calculating the intensity of selection on mutations along with tallying the prevalence of individual substitutions in cancer: while many of the most prevalently-altered genes were heavily selected, their relative importance to the cancer phenotype differs from their prevalence and from their P value, with some infrequent variants exhibiting evidence of strong positive selection. Furthermore, we extend our analysis of effect size by quantifying the degree to which mutational processes (such as APOBEC mutagenesis) contributes mutations that are highly selected, driving head and neck squamous cell carcinoma. We calculate the substitutions caused by APOBEC mutagenesis that make the greatest contribution to cancer phenotype among patients. Lastly, we demonstrate via in vitro biochemical experiments that the APOBEC3B protein can deaminate the cytosine bases at two sites whose mutant states are subject to high net realized selection intensities—PIK3CA E545K and E542K. By quantifying the effects of mutations, we deepen the molecular understanding of carcinogenesis in head and neck squamous cell carcinoma.

Purpose: Human papillomavirus (HPV)-negative head and neck squamous cell carcinomas (HNSCC) commonly bear disruptive mutations in TP53, resulting in treatment resistance. In these patients, direct targeting of p53 has not been successful, but synthetic lethal approaches have promise. Although Aurora A kinase (AURKA) is overexpressed and an oncogenic driver, its inhibition has only modest clinical effects in HPV-negative HNSCC. We explored a novel combination of AURKA and WEE1 inhibition to overcome intrinsic resistance to AURKA inhibition. Experimental Design: AURKA protein expression was determined by fluorescence-based automated quantitative analysis of patient specimens and correlated with survival. We evaluated treatment with the AURKA inhibitor alisertib (MLN8237) and the WEE1 inhibitor adavosertib (AZD1775), alone or in combination, using in vitro and in vivo HNSCC models. Results: Elevated nuclear AURKA correlated with worse survival among patients with p16(−) HNSCC. Alisertib caused spindle defects, G₂–M arrest and inhibitory CDK1 phosphorylation, and cytostasis in TP53 mutant HNSCC FaDu and UNC7 cells. Addition of adavosertib to alisertib instead triggered mitotic entry and mitotic catastrophe. Moreover, in FaDu and Detroit 562 xenografts, this combination demonstrated synergistic effects on tumor growth and extended overall survival compared with either vehicle or single-agent treatment. Conclusions: Combinatorial treatment with adavosertib and alisertib leads to synergistic antitumor effects in in vitro and in vivo HNSCC models. These findings suggest a novel rational combination, providing a promising therapeutic avenue for TP53-mutated cancers.

Cancer progression is an evolutionary process. During this process, evolving cancer cell populations encounter restrictive ecological niches within the body, such as the primary tumor, circulatory system, and diverse metastatic sites. Efforts to prevent or delay cancer evolution—and progression—require a deep understanding of the underlying molecular evolutionary processes. Herein we discuss a suite of concepts and tools from evolutionary and ecological theory that can inform cancer biology in new and meaningful ways. We also highlight current challenges to applying these concepts, and propose ways in which incorporating these concepts could identify new therapeutic modes and vulnerabilities in cancer.