Zi-Ming Zhao

Early and multiple origins of metastatic lineages within primary tumors

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. Many aspects of the evolutionary process of tumorigenesis that are fundamental to cancer biology and targeted treatment have been challenging to reveal, such as the divergence times and genetic clonality of metastatic lineages. To address these challenges, we performed tumor phylogenetics using molecular evolutionary models, reconstructed ancestral states of somatic mutations, and inferred cancer chronograms to yield three conclusions. First, in contrast to a linear model of cancer progression, metastases can originate from divergent lineages within primary tumors. Evolved genetic changes in cancer lineages likely affect only the proclivity toward metastasis. Single genetic changes are unlikely to be necessary or sufficient for metastasis. Second, metastatic lineages can arise early in tumor development, sometimes long before diagnosis. The early genetic divergence of some metastatic lineages directs attention toward research on driver genes that are mutated early in cancer evolution. Last, the temporal order of occurrence of driver mutations can be inferred from phylogenetic analysis of cancer chronograms, guiding development of targeted therapeutics effective against primary tumors and metastases.

Whole-exome sequencing and immune profiling of early-stage lung adenocarcinoma with fully annotated clinical follow-up

Background 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. Methods 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 Fishers exact test. Cox proportional hazards regression models were used for multivariate analysis of clinical outcome. Results 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. Conclusion(s) Our study highlights molecular and immune phenotypes that warrant further analysis for their roles in clinical outcomes and personalized immune-based therapy of LUAD.BACKGROUND 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. METHODS 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 Fishers exact test. Cox proportional hazards regression models were used for multivariate analysis of clinical outcome. RESULTS 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. CONCLUSION(S) Our study highlights molecular and immune phenotypes that warrant further analysis for their roles in clinical outcomes and personalized immune-based therapy of LUAD.

Mutation profiles in early-stage lung squamous cell carcinoma with clinical follow-up and correlation with markers of immune function.

Background 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. Methods 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. Results 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. Conclusion(s) Our findings pinpoint mutated genes that may impact clinical outcome as well as personalized strategies for targeted immunotherapies in early-stage LUSC.

Detection of Regional Variation in Selection Intensity within Protein-Coding Genes Using DNA Sequence Polymorphism and Divergence

Numerous approaches have been developed to infer natural selection based on the comparison of polymorphism within species and divergence between species. These methods are especially powerful for the detection of uniform selection operating across a gene. However, empirical analyses have demonstrated that regions of protein-coding genes exhibiting clusters of amino acid substitutions are subject to different levels of selection relative to other regions of the same gene. To quantify this heterogeneity of selection within coding sequences, we developed Model Averaged Site Selection via Poisson Random Field (MASS-PRF). MASS-PRF identifies an ensemble of intragenic clustering models for polymorphic and divergent sites. This ensemble of models is used within the Poisson Random Field framework to estimate selection intensity on a site-by-site basis. Using simulations, we demonstrate that MASS-PRF has high power to detect clusters of amino acid variants in small genic regions, can reliably estimate the probability of a variant occurring at each nucleotide site in sequence data and is robust to historical demographic trends and recombination. We applied MASS-PRF to human gene polymorphism derived from the 1,000 Genomes Project and divergence data from the common chimpanzee. On the basis of this analysis, we discovered striking regional variation in selection intensity, indicative of positive or negative selection, in well-defined domains of genes that have previously been associated with neurological processing, immunity, and reproduction. We suggest that amino acid-altering substitutions within these regions likely are or have been selectively advantageous in the human lineage, playing important roles in protein function.

Heterogeneity and mutation in KRAS and associated oncogenes: evaluating the potential for the evolution of resistance to targeting of KRAS G12C

Activating mutations in RAS genes are associated with approximately 20% of all human cancers. New targeted therapies show preclinical promise in inhibiting the KRAS G12C variant. However, concerns exist regarding the effectiveness of such therapies in vivo given the possibilities of existing intratumor heterogeneity or de novo mutation leading to treatment resistance. We performed deep sequencing of 27 KRAS G12-positive lung tumors to determine the prevalence of other oncogenic mutations within KRAS or within commonly mutated downstream genes that could confer resistance at the time of treatment. We also passaged patient-derived xenografts to assess the potential for novel KRAS mutation to arise during subsequent tumor evolution. Furthermore, we estimate the de novo mutation rate in KRAS position 12 and in genes downstream of KRAS. Finally, we present an approach for estimation of the selection intensity for these point mutations that explains their high prevalence in tumors. We find no evidence of heterogeneity that may compromise KRAS G12C targeted therapy within sequenced lung tumors or passaged xenografts. We find that mutations that confer resistance are even less likely to occur downstream of KRAS than to occur within KRAS. Our approach predicts that BRAF V600E would provide the highest fitness advantage for de novo-resistant subclones. Overall, our findings suggest that resistance to targeted therapy of KRAS G12C-positive tumors is unlikely to be present at the time of treatment and, among the de novo mutations likely to confer resistance, mutations in BRAF, a gene with targeted inhibitors presently available, result in subclones with the highest fitness advantage.

The likelihood of heterogeneity or additional mutation in KRAS or associated oncogenes to compromise targeting of oncogenic KRAS G12C

Activating mutations in RAS genes are associated with approximately 20% of all human cancers. New targeted therapies show preclinical promise in inhibiting the KRAS G12C variant, however, concerns exist regarding the effectiveness of such therapies in vivo given the possibilities of existing intratumor heterogeneity or de novo mutation leading to treatment resistance. We performed deep sequencing of 27 KRAS G12 positive lung tumors and found no evidence of other oncogenic mutations within KRAS or within commonly mutated downstream genes that could confer resistance at the time of treatment. Furthermore, we estimate the de novo mutation rate in KRAS position 12 and in genes downstream of KRAS. We find that mutations that confer resistance are about as likely to occur downstream of KRAS as within KRAS. Moreover, we present an approach for estimation of the selection intensity for these point mutations that explains their high prevalence in tumors. Our approach predicts that BRAF V600E would provide the highest fitness advantage for de novo resistant subclones. Overall, our findings suggest that resistance to targeted therapy of KRAS G12C positive tumors is unlikely to be present at the time of treatment and, among the de novo mutations likely to confer resistance, mutations in BRAF, a gene with targeted inhibitors presently available, result in subclones with the highest fitness advantage. One Sentence Summary Mutations conferring resistance to KRAS G12C targeted therapy are unlikely to be present at the time of resection, and the likely mechanisms of evolved resistance are predicted be ones that are responsive to therapies that are in development or that are already available.

Abstract 142: Mutation and immune profiles in early-stage lung squamous cell carcinoma

PURPOSE: 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 sought to perform whole-exome sequencing (WES), comprehensive immune profiling and clinicopathological analysis of early-stage LUSCs to increase our understanding of the pathobiology of this malignancy. METHODS: 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 was done on a subset of LUSCs (n = 91). Associations among mutations, immune markers and clinicopathological variables were statistically examined using ANOVA and Fisher tests. Cox proportional hazards regression models were used for statistical analysis of clinical outcome. RESULTS: This early-stage LUSC cohort displayed an average of 209 exonic mutations per tumor. Fourteen genes exhibited significant enrichment for 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. CONCLUSION: Our findings pinpoint mutated genes that may impact clinical outcome as well as personalized strategies for targeted immunotherapies in early-stage LUSC. Citation Format: Murim Choi, Humam Kadara, Jiexin Zhang, Edwin Parra Cuentas, Jaime Rodriguez Canales, Stephen G. Gaffney, Zi-Ming Zhao, Carmen Behrens, Junya Fujimoto, Chi-Wan Chow, Neda Kalhor, Cesar Moran, David Rimm, Stephen Swisher, Don L. Gibbons, John V. Heymach, Edward Kaftan, Jeffrey Townsend, Thomas J. Lynch, Joseph Schlessinger, J. Jack Lee, Richard Lifton, Roy S. Herbst, Ignacio I. Wistuba. Mutation and immune profiles in early-stage lung squamous cell carcinoma. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 142.

Abstract 89: Whole-exome sequencing and immune profiling of early-stage lung adenocarcinoma

PURPOSE: Lung adenocarcinomas (LUADs) lead to the preponderance of deaths attributable to lung cancer. We performed whole-exome sequencing (WES), comprehensive immune profiling and clinicopathological analysis of LUADs to better understand the molecular pathogenesis of this disease and identify clinically relevant molecular markers. METHODS: We performed WES of 108 paired surgically resected stage I-III LUADs and normal lung tissues using the Illumina HiSeq 2000 platform. Additionally, ten immune related markers (PD-L1, PD-1, CD3, CD4, CD8, CD45ro, CD57, CD68, FOXP3 and Granzyme B) were profiled by imaging-based immunohistochemistry in a subset of LUADs (n = 92). Associations among mutations, immune markers and clinicopathological variables were analyzed using ANOVA and Fishers Exact tests. Cox proportional hazards regression models were employed for multivariate analysis of clinical outcome. RESULTS: LUADs in this cohort exhibited an average of 243 coding mutations per tumor. We identified 28 genes with significant enrichment for mutation. SETD2-mutant LUADS exhibited relatively poor recurrence-free survival (RFS) and mutations in STK11 and ATM were associated with poor RFS in KRAS-mutant tumors. EGFR, KEAP1 and PIK3CA mutations were predictive of poor response to adjuvant therapy. Immune marker analysis demonstrated that PD-L1 expression was increased in smoker compared to non-smoker LUADs and, along with other immune markers, was positively correlated with somatic mutation burden. Moreover, immune marker levels including PD-L1 were elevated in TP53-mutant LUADs. In contrast, STK11 and U2AF1 mutant tumors exhibited a suppressed immune response and LUADs with PIK3CA mutations exhibited markedly decreased tumoral PD-L1 expression. CONCLUSION: Our study highlights mutations that may impact clinical outcome and personalized strategies for immune-based therapy of early-stage LUAD patients. Citation Format: Humam Kadara, Murim Choi, Jiexin Zhang, Edwin Parra Cuentas, Jaime Rodriguez Canales, Stephen Gaffney, Zi-Ming Zhao, Carmen Behrens, Junya Fujimoto, Chi-Wan Chow, Neda Kalhor, Cesar Moran, David Rimm, Stephen G. Swisher, Don L. Gibbons, John V. Heymach, Edward Kaftan, Jeffrey Townsend, Thomas J. Lynch, Joseph Schlessinger, J. Jack Lee, Richard Lifton, Ignacio I. Wistuba, Roy S. Herbst. Whole-exome sequencing and immune profiling of early-stage lung adenocarcinoma. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 89.