Research Paper Volume 12, Issue 3 pp 2132—2141
EP300 mutation is associated with tumor mutation burden and promotes antitumor immunity in bladder cancer patients
- 1 Department of Urology, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 2 The State Key Laboratory of Biotherapy, West China Hospital, Sichuan University, Chengdu, China
- 3 Department of Cardiothoracic Surgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
received: October 10, 2019 ; accepted: January 2, 2020 ; published: February 3, 2020 ;https://doi.org/10.18632/aging.102728
How to Cite
Copyright © 2020 Zhu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Bladder cancer is a leading cause of morbidity and mortality worldwide. Currently, immunotherapy has become a worthwhile therapy for bladder cancer. Tumor mutation burden (TMB) has been regarded as the most prevalent biomarker to predict immunotherapy. Bladder cancer is reported to have the third highest mutation rate. However, whether these gene mutations are related to TMB and immune response remain unknown. In this study, we downloaded somatic mutation data of bladder cancer from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) datasets, and found 11 frequently mutated genes were covered by both two cohorts including FGFR3, TTN, XIRP2, CREBBP, PIK3CA, TP53, MUC16, EP300 (E1A binding protein P300), ARID1A, ERBB2, and KDM6A. Among them, EP300 mutation was associated with higher TMB and indicated a favorable clinical prognosis. Furthermore, based on Gene set enrichment analysis (GSEA) and CIBERSORT algorithm, we observed that EP300 mutation upregulated signaling pathways involved in immune system and enhanced antitumor immune response. In conclusion, EP300 is frequently mutated in bladder cancer, and its mutation is associated with increased TMB and promotes antitumor immunity, which may serve as a biomarker to predict immune response.
Worldwide, bladder cancer accounts for approximately 3.0% of all malignancies. In 2018, an estimated 549,393 bladder cancer cases were identified and contributed to 199,922 deaths . Although approximately 70% of bladder cancer patients present with the non-muscle invasive type, as many as half of them experience disease recurrence after receiving transurethral resection, and up to quarter of them progress to muscle-invasive bladder cancer (MIBC) after repeated recurrences . Moreover, approximately 30% of patients have muscle invasive disease at the time of initial diagnosis, and radical cystectomy is the current standard of care . Nevertheless, even after surgery, approximately 50% of patients develop metastases . Over the past years, systemic chemotherapy using cisplatin-based regimens has been the predominant treatment for metastatic urothelial carcinoma (mUC). However, approximately 50% of patients with mUC are ineligible to receive cisplatin due to severe adverse effects . Therefore, for patients for whom first-line chemotherapy is not successful, there is an urgent need to find a novel treatment option to improve clinical outcomes.
Currently, with the emergence of immune checkpoint inhibitors targeting programmed cell death protein and its ligand (PD-1/PD-L1), immunotherapy has revolutionized anti-cancer therapy. Bladder cancer has been recognized as an immunoresponsive tumor, and immunotherapy has been reported to be a worthwhile neoadjuvant therapy for treatment of MIBC [6, 7]. However, only a minority of patients benefit from these immunotherapies, and fewer achieve a durable response . Therefore, understanding the molecular determinants of immunotherapeutic responsiveness is one of critical challenges in bladder cancer. Emerging studies have explored the correlation between tumor mutation burden (TMB) and immunotherapy response [9, 10]. Accumulation of somatic mutations is one of the leading cause of tumorigenesis and contributes to the expression of neoantigens . Neoantigens can activate CD8+ CTLs to act against tumors by recognizing target antigens presented on tumor cells . Previous studies have reported that a high TMB predicts a better outcome in non-small cell lung cancer and melanoma [13, 14]. Meanwhile, TMB has also been regarded as a predictive biomarker of tumor behavior and immunological response [9, 15]. Previous study has demonstrated that bladder cancer has the third highest mutation prevalence and is highly antigenic . However, whether these gene mutations are related to TMB and immune response in bladder cancer remains unclear.
In this study, we first identified somatic mutations in bladder cancer patients from US and China using The Cancer Genome Atlas (TCGA) dataset and International Cancer Genome Consortium (ICGC) dataset. Then, we found the common mutant genes in the two cohorts, and further explored the association of these gene mutations with TMB and prognosis. Finally, we investigated whether gene mutations are associated with immune response. Findings emerged from this study may identify a novel biomarker for TMB and suggest a potential immunotherapy for bladder cancer patients.
Landscape of somatic mutations in bladder cancer
We defined 30 frequently mutated genes in American bladder cancer samples from TCGA cohort, and the five most frequently mutated genes were TP53 (47%), TTN (40%), KMT2D (26%), KDM6A (25%), and ARID1A (24%) (Figure 1A). Meanwhile, we also defined 30 frequently mutated genes in Chinese bladder cancer samples from ICGC cohort, and the five most frequently mutated genes were MUC4 (58%), TP53 (24%), FMN2 (23%), PIK3CA (23%), and KDM6A (21%) (Figure 1B). Interestingly, we observed that some genes frequently mutated in both American and Chinese patients. Therefore, we performed comparative analysis of 30 frequently mutated genes between TCGA and ICGC cohorts of bladder cancer. 11 frequently mutated genes reported by TCGA cohort were covered by ICGC cohort, including FGFR3, TTN, XIRP2, CREBBP, PIK3CA, TP53, MUC16, EP300, ARID1A, ERBB2, and KDM6A (Figure 1C). Then, we focused on these common mutated genes in subsequent analysis.
Figure 1. Landscapes of frequently mutated genes in bladder cancer. (A) Oncoplot depicts the frequently mutated genes in bladder cancer from TCGA cohort. The left panel shows mutation frequency, and genes are ordered by their mutation frequencies. The bottom panel presents different mutation types. (B) Waterfall plot displays the frequently mutated genes in bladder cancer from ICGC cohort. The left panel shows the genes ordered by their mutation frequencies. The right panel presents different mutation types. (C) Venn diagram of frequently mutated genes covered by both TCGA and ICGC cohorts.
EP300 mutation associated with TMB and survival prognosis
The TMB score in bladder cancer is ranging from 0.02 to 130 per Mb with a median of 5.06 per Mb. Among common mutated genes, patients with mutation in TTN, XIRP2, PIK3CA, TP53, MUC16, EP300, ARID1A, ERBB2, and KDM6A had significantly higher TMB (Figure 2A). Previous study has reported that higher TMB indicated a favorable overall survival in bladder cancer patients . Therefore, we performed Kaplan-Meier analysis to identify whether these gene mutations associated with increased TMB were also related to the prognosis of patients with bladder cancer. As shown in Figure 2B, only EP300 mutation (HR=0.612; 95% CI, 0.399-0.937; P=0.024) was associated with a positive prognosis. However, EP300 mutation did not remain statistically significance after taking into account age, gender, grade, TNM classification and TMB status in Cox regression model (Table 1), suggesting that EP300 mutation might not be an independent risk factor of prognosis in bladder cancer patients.
Figure 2. Gene mutations are associated with TMB and clinical prognosis. (A) Most gene mutations are associated with a higher TMB. ** p<0.01; *** p<0.001. (B) Kaplan-Meier survival analysis of patients with gene mutations. Only patients with complete clinical information were included (n=402). The p-value is marked in each graph. WT, wild type; MT, mutant type.
Table 1. Univariate and multivariate overall survival analysis of bladder cancer patients by the COX proportional hazards model.
|HR (95% CI)||P-value||HR (95% CI)||P-value|
|Age (year) (<70, ≥70)||1.704 (1.242-2.339)||<0.001||1.765 (1.284-2.425)||<0.001|
|Gender (male, female)||0.922 (0.650-1.307)||0.648|
|Grade (low, high)||21.636 (0.273-1715.747)||0.168|
|TNM classification||2.533 (1.672-3.838)||<0.001||2.479 (1.635-3.758)||<0.001|
|(stage I&II, stage III&IV)|
|TMB (low, high)||0.943 (0.915-0.972)||<0.001||0.947 (0.918-0.977)||<0.001|
|EP300 (wide, mutant)||0.575 (0.338-0.980)||0.042||0.665 (0.385-1.149)||0.144|
Enrichment pathway analysis of EP300 mutation
As TMB is reported to be a biomarker for immunotherapy and EP300 mutation was associated with an increased TMB, we further investigated the relation of EP300 mutation and immune response. GSEA performed with TCGA revealed that Biocarta IL2 Pathway, Biocarta Insulin Pathway, Kegg Natural Killer Cell Mediated Cytotoxicity, Pid IL12 Stat4 Pathway, Reactome Cytokine Signaling in Immune System, and Reactome MHC II Antigen Presentation were significantly enriched in samples with EP300 mutation (Figure 3A–3F). These findings indicated that samples with EP300 mutation upregulated signaling pathways involved in immune system.
Figure 3. Significantly enriched pathways associated with EP300 mutation. Gene set enrichment analysis has been performed with TCGA. Gene enrichment plots shows that a series of gene sets including (A) Biocarta IL2 Pathway, (B) Biocarta Insulin Pathway, (C) Kegg Natural Killer Cell Mediated Cytotoxicity, (D) Pid IL12 Stat4 Pathway, (E) Reactome Cytokine Signaling in Immune System, and (F) Reactome MHC II Antigen Presentation are enriched in EP300-mutant group. NES, normalized enrichment score. The p-value is marked in each plot.
Tumor-infiltrating immune cells associated with EP300 mutation in bladder cancer
We next evaluated the association of EP300 mutation and tumor-infiltrating immune cells in bladder cancer microenvironment using CIBERSORT algorithm. As shown in Figure 4A, the results showed that the composition of 22 immune cell types in each sample varied significantly. Moreover, we observed that activated memory CD4 T cells and resting NK cells were more enriched in EP300 mutant type group, nevertheless, regulatory T cells (Tregs) and gamma delta T cells were enriched in wild-type group (Figure 4B). Furthermore, the result from the correlation matrix revealed that activated memory CD4 T cells had the strongest positive correlation with CD8 T cells, and also positively correlated with resting NK cells. However, activated CD4 memory T cells had negative correlation with Tregs (Figure 4C).
Figure 4. EP300 mutation is correlated with tumor-infiltrating immune cells. (A) Stacked bar chart shows distribution of 22 immune cells in each sample. (B) Violin plot displays the differentially infiltrated immune cells between EP300-mutant groups and EP300-wild group. Blue color represents EP300-wild group, and red color represents EP300-mutant group. (C) Correlation matrix of immune cell proportions. The red color represents positive correlation and the blue color represents negative correlation.
We respectively characterized the somatic mutation landscape of 412 American bladder cancer samples from TCGA dataset and 101 Chinese bladder cancer samples from ICGC dataset. Then, we found that EP300 was frequently mutated in both TCGA and ICGC cohorts, and its mutation was associated with higher TMB and favorable clinical prognosis. Meanwhile, samples with EP300 mutation positive correlated with signaling pathways implicated in immune response. Tumor-infiltrating immune cells results demonstrated that EP300 mutated samples were more infiltrated in activated memory CD4 T cells and resting NK cells, and less infiltrated in Tregs and gamma delta T cells, which supported the previous findings that such immune cells and pathways play predominant roles in the tumor microenvironment and promote the immune response [18, 19].
EP300 is known to encode an adenoviral E1A-binding proteins (also known as p300) involved in multiple cellular processes functioning as a transcriptional co-factor and histone acetyltransferase (KAT) as well as its homologous protein partner cAMP response element binding protein (CBP) . EP300 mutation has been reported in multiple types of human cancer, including breast cancer and squamous cell carcinoma of the head and neck [21, 22]. Specific mutations in EP300 predict a lower risk for relapses and decrease the mortality of breast cancer . Bladder cancer is a genomically heterogeneous tumor with a high mutation rate, and EP300 is also a frequently mutated gene in bladder cancer. In our study, we found that mutation in EP300 indicated a better prognosis and was associated with increased TMB. TMB represents the accumulation of somatic mutations in tumors, and a high TMB contributes to the exposure of more neoantigens, which is likely to elicit a T-cell-dependent immune response . Therefore, we speculated that EP300 mutation might enhance immune response.
Previous study indicated that as one of the three main KAT families, the p300-CBP family, and p300 in particular, could acetylate and stabilize Forkhead box P3 (Foxp3) protein expression, which was the key transcription factor of Tregs . Tregs are known as adaptive immune cells that maintain self-tolerance and impose a critical check on excessive immune response . Increased Tregs infiltrations are associated with poor outcome in several cancer types, including breast, kidney, and lung cancer [26–28]. Therefore, inhibition of p300 may impair Treg homeostasis and thereby promote antitumor immunity. Recent studies have confirmed that deletion of EP300 in Tregs suppressed the proliferation and function of Tregs and limited the ability of cancers to avoid immune destruction, thereby decreasing tumor growth [29, 30]. In our study, we also found that bladder cancer samples with EP300 mutation had lower Tregs loads. Thus, we speculated that EP300 mutation might impair the function of p300 and affect the generation of Tregs. Moreover, we observed that EP300 mutation was associated with increased memory CD4 T cells and NK cells infiltrations. The frequency of CD4+ lymphocytes with a memory phenotype is reported to significantly increase in patients with tumor free lymph node, and the generation of memory response against tumor antigens and may prevent tumor relapse in patients with different cancers, which predict a longer survival [31, 32]. NK cells as a significant part of innate immune system are characterized by secreting cytokines and cytolytic activity against target cells. Recent studies have shown that NK cells can selectively kill cancer stem cells, which indicate that NK cell-based therapy may become an effective treatment to inhibit tumor relapse and metastasis [33, 34]. Therefore, we found that EP300 mutation elicited the variation of infiltrated immune cells contributing to antitumor immunity in bladder cancer.
The main limitation in our study is that the ICGC database lacks corresponding clinical data of Chinese bladder cancer, so we cannot verify whether EP300 mutation is associated with the prognosis of bladder cancer patients in China, and whether it can give rise to the same immune response. Even if EP300 is also frequently mutated in Chinese bladder cancer sample, the effect may be somewhat heterogeneous between different ethnic groups. As a result, the association of EP300 mutation and prognosis, including analysis in infiltrated immune cells and signaling pathways need further validation in Chinese bladder samples.
In summary, our study demonstrated that EP300 was frequently mutated in bladder cancer, and EP300 mutation was associated with higher TMB and indicated a better prognosis. Furthermore, EP300 mutation upregulated the signaling pathways of the immune system and evoked an antitumor immune response. These findings reveal a novel gene whose mutation could be served as a biomarker to predict immune response.
Materials and Methods
Somatic gene mutations for American bladder cancer samples (n=412) and Chinese bladder cancer samples (n=101) were respectively acquired from TCGA portal (http://portal.gdc.cancer.gov/projects) (up to June 10, 2019) and ICGC portal (http://dcc.icgc.org/ releases/current/Projects) (up to April 3, 2019). Clinical data for 412 bladder cancer samples were downloaded from TCGA. For clinical data, only patients with bladder cancer with complete information (n=402) were included, and those with any missing data concerning survival time, status, age, gender, grade, or TNM classification were excluded.
Classification of bladder cancer based on TMB
TMB was defined as the number of somatic, coding, indels mutations and base substitution per megabase of genome examined. All base substitutions and idels in the coding region of targeted genes were counted. Silent mutations failing to contribute to an amino acid change were not counted. To calculate the TMB score of each sample, the total number of mutations counted was divided by the exome size (38 megabase (Mb) was used as the estimate of the exome size) .
MAF files containing somatic variants for American bladder cancer samples were detected using VarScan and visualized with maftools package . TSV files containing somatic variants for Chinese bladder cancer samples were annotated according to hg19 reference genome and visualized with GenVisR package. Gene set enrichment analysis (GSEA) was performed using Broad Institute GSEA software 3.0 according to previously described . The Gene expression data was downloaded from TCGA portal. Patients were available into two groups according to mutation status of E1A binding protein P300 (EP300). The gene set “c2 all.v6.0 symbols.gmt” was downloaded from Molecular Signatures Database (http://software.broadinstitute.org/gsea/msigdb/index.jsp) and was used for the enrichment analysis. Permutations were set to 1000 to obtain a normalized enrichment score (NES). A normal p-value <0.05 was considered significantly enriched. CIBERSORT , a deconvolution algorithm that evaluates the proportions of 22 tumor-infiltrating lymphocyte subsets in a bulk tumor sample, was used to estimate the relative abundance of immune cell infiltration in patients with different EP300 status. The number of permutations was set to 1000, and a threshold p-value of <0.05 was the criterion for successful computation of a sample.
Statistical analyses were performed with R software (version 3.6.0) and GraphPad Prism 5.0 software (San Diego, CA, USA). Kaplan-Meier survival analysis was used to determine survival curves that reflect the association between gene mutations and prognosis, which were evaluated by the log-rank test. Univariate and multivariate Cox regression analyses were used for survival analysis of clinical characteristics of patients, including age, gender, grade, TNM classification, TMB and EP300. Mann-Whitney U test was used to analyze the correlation between mutant genes and TMB. For all comparisons, a two-tailed p-value <0.05 was considered statistically significant.
EP300: E1A binding protein P300; TMB: tumor mutation burden; TCGA: The Cancer Genome Atlas; ICGC: International Cancer Genome Consortium; GSEA: Gene set enrichment analysis; MIBC: muscle-invasive bladder cancer; mUC: metastatic urothelial carcinoma; Mb: megabase; NES: normalized enrichment score; KAT: histone acetyltransferase; CBP: cAMP response element binding protein; Foxp3: Forkhead box P3.
X. Gou designed the research; G. Zhu and X. Gou prepared the figures and drafted the manuscript; G. Zhu, L. Pei, Y. Li analyzed the data; X. Gou contributed analytic tools and finalized the manuscript. All authors have read and approved the final manuscript.
Conflicts of Interest
The authors declare that there are no conflicts of interest
- 1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018; 68:394–424. https://doi.org/10.3322/caac.21492 [PubMed]
- 2. Cookson MS, Herr HW, Zhang ZF, Soloway S, Sogani PC, Fair WR. The treated natural history of high risk superficial bladder cancer: 15-year outcome. J Urol. 1997; 158:62–67. https://doi.org/10.1097/00005392-199707000-00017 [PubMed]
- 3. Milowsky MI, Rumble RB, Booth CM, Gilligan T, Eapen LJ, Hauke RJ, Boumansour P, Lee CT. Guideline on Muscle-Invasive and Metastatic Bladder Cancer (European Association of Urology Guideline): American Society of Clinical Oncology Clinical Practice Guideline Endorsement. J Clin Oncol. 2016; 34:1945–52. https://doi.org/10.1200/JCO.2015.65.9797 [PubMed]
- 4. Ploussard G, Shariat SF, Dragomir A, Kluth LA, Xylinas E, Masson-Lecomte A, Rieken M, Rink M, Matsumoto K, Kikuchi E, Klatte T, Boorjian SA, Lotan Y, et al. Conditional survival after radical cystectomy for bladder cancer: evidence for a patient changing risk profile over time. Eur Urol. 2014; 66:361–70. https://doi.org/10.1016/j.eururo.2013.09.050 [PubMed]
- 5. Dash A, Galsky MD, Vickers AJ, Serio AM, Koppie TM, Dalbagni G, Bochner BH. Impact of renal impairment on eligibility for adjuvant cisplatin-based chemotherapy in patients with urothelial carcinoma of the bladder. Cancer. 2006; 107:506–13. https://doi.org/10.1002/cncr.22031 [PubMed]
- 6. Morales A, Eidinger D, Bruce AW. Intracavitary Bacillus Calmette-Guerin in the Treatment of Superficial Bladder Tumors. J Urol. 2017; 197:S142–45. https://doi.org/10.1016/j.juro.2016.10.101 [PubMed]
- 7. Necchi A, Anichini A, Raggi D, Briganti A, Massa S, Lucianò R, Colecchia M, Giannatempo P, Mortarini R, Bianchi M, Farè E, Monopoli F, Colombo R, et al. Pembrolizumab as Neoadjuvant Therapy Before Radical Cystectomy in Patients With Muscle-Invasive Urothelial Bladder Carcinoma (PURE-01): An Open-Label, Single-Arm, Phase II Study. J Clin Oncol. 2018; 36:JCO1801148. https://doi.org/10.1200/JCO.18.01148 [PubMed]
- 8. Braun DA, Burke KP, Van Allen EM. Genomic Approaches to Understanding Response and Resistance to Immunotherapy. Clin Cancer Res. 2016; 22:5642–50. https://doi.org/10.1158/1078-0432.CCR-16-0066 [PubMed]
- 9. Goodman AM, Kato S, Bazhenova L, Patel SP, Frampton GM, Miller V, Stephens PJ, Daniels GA, Kurzrock R. Tumor Mutational Burden as an Independent Predictor of Response to Immunotherapy in Diverse Cancers. Mol Cancer Ther. 2017; 16:2598–608. https://doi.org/10.1158/1535-7163.MCT-17-0386 [PubMed]
- 10. Rizvi NA, Hellmann MD, Snyder A, Kvistborg P, Makarov V, Havel JJ, Lee W, Yuan J, Wong P, Ho TS, Miller ML, Rekhtman N, Moreira AL, et al. Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science. 2015; 348:124–28. https://doi.org/10.1126/science.aaa1348 [PubMed]
- 11. Gubin MM, Artyomov MN, Mardis ER, Schreiber RD. Tumor neoantigens: building a framework for personalized cancer immunotherapy. J Clin Invest. 2015; 125:3413–21. https://doi.org/10.1172/JCI80008 [PubMed]
- 12. Chen L, Flies DB. Molecular mechanisms of T cell co-stimulation and co-inhibition. Nat Rev Immunol. 2013; 13:227–42. https://doi.org/10.1038/nri3405 [PubMed]
- 13. Chen H, Chong W, Wu Q, Yao Y, Mao M, Wang X. Association of LRP1B Mutation With Tumor Mutation Burden and Outcomes in Melanoma and Non-small Cell Lung Cancer Patients Treated With Immune Check-Point Blockades. Front Immunol. 2019; 10:1113. https://doi.org/10.3389/fimmu.2019.01113 [PubMed]
- 14. Chen Y, Liu Q, Chen Z, Wang Y, Yang W, Hu Y, Han W, Zeng H, Ma H, Dai J, Zhang H. PD-L1 expression and tumor mutational burden status for prediction of response to chemotherapy and targeted therapy in non-small cell lung cancer. J Exp Clin Cancer Res. 2019; 38:193. https://doi.org/10.1186/s13046-019-1192-1 [PubMed]
- 15. Birkbak NJ, Kochupurakkal B, Izarzugaza JM, Eklund AC, Li Y, Liu J, Szallasi Z, Matulonis UA, Richardson AL, Iglehart JD, Wang ZC. Tumor mutation burden forecasts outcome in ovarian cancer with BRCA1 or BRCA2 mutations. PLoS One. 2013; 8:e80023. https://doi.org/10.1371/journal.pone.0080023 [PubMed]
- 16. Alexandrov LB, Nik-Zainal S, Wedge DC, Aparicio SA, Behjati S, Biankin AV, Bignell GR, Bolli N, Borg A, Børresen-Dale AL, Boyault S, Burkhardt B, Butler AP, et al, and Australian Pancreatic Cancer Genome Initiative, and ICGC Breast Cancer Consortium, and ICGC MMML-Seq Consortium, and ICGC PedBrain. Signatures of mutational processes in human cancer. Nature. 2013; 500:415–21. https://doi.org/10.1038/nature12477 [PubMed]
- 17. Aggen DH, Drake CG. Biomarkers for immunotherapy in bladder cancer: a moving target. J Immunother Cancer. 2017; 5:94. https://doi.org/10.1186/s40425-017-0299-1 [PubMed]
- 18. Zuazo M, Arasanz H, Fernández-Hinojal G, García-Granda MJ, Gato M, Bocanegra A, Martínez M, Hernández B, Teijeira L, Morilla I, Lecumberri MJ, Fernández de Lascoiti A, Vera R, et al. Functional systemic CD4 immunity is required for clinical responses to PD-L1/PD-1 blockade therapy. EMBO Mol Med. 2019; 11:e10293. https://doi.org/10.15252/emmm.201910293 [PubMed]
- 19. Hsu J, Hodgins JJ, Marathe M, Nicolai CJ, Bourgeois-Daigneault MC, Trevino TN, Azimi CS, Scheer AK, Randolph HE, Thompson TW, Zhang L, Iannello A, Mathur N, et al. Contribution of NK cells to immunotherapy mediated by PD-1/PD-L1 blockade. J Clin Invest. 2018; 128:4654–68. https://doi.org/10.1172/JCI99317 [PubMed]
- 20. Iyer NG, Ozdag H, Caldas C. p300/CBP and cancer. Oncogene. 2004; 23:4225–31. https://doi.org/10.1038/sj.onc.1207118 [PubMed]
- 21. Bemanian V, Noone JC, Sauer T, Touma J, Vetvik K, Søderberg-Naucler C, Lindstrøm JC, Bukholm IR, Kristensen VN, Geisler J. Somatic EP300-G211S mutations are associated with overall somatic mutational patterns and breast cancer specific survival in triple-negative breast cancer. Breast Cancer Res Treat. 2018; 172:339–51. https://doi.org/10.1007/s10549-018-4927-3 [PubMed]
- 22. Saloura V, Izumchenko E, Zuo Z, Bao R, Korzinkin M, Ozerov I, Zhavoronkov A, Sidransky D, Bedi A, Hoque MO, Koeppen H, Keck MK, Khattri A, et al. Immune profiles in primary squamous cell carcinoma of the head and neck. Oral Oncol. 2019; 96:77–88. https://doi.org/10.1016/j.oraloncology.2019.06.032 [PubMed]
- 23. McGranahan N, Furness AJ, Rosenthal R, Ramskov S, Lyngaa R, Saini SK, Jamal-Hanjani M, Wilson GA, Birkbak NJ, Hiley CT, Watkins TB, Shafi S, Murugaesu N, et al. Clonal neoantigens elicit T cell immunoreactivity and sensitivity to immune checkpoint blockade. Science. 2016; 351:1463–69. https://doi.org/10.1126/science.aaf1490 [PubMed]
- 24. van Loosdregt J, Vercoulen Y, Guichelaar T, Gent YY, Beekman JM, van Beekum O, Brenkman AB, Hijnen DJ, Mutis T, Kalkhoven E, Prakken BJ, Coffer PJ. Regulation of Treg functionality by acetylation-mediated Foxp3 protein stabilization. Blood. 2010; 115:965–74. https://doi.org/10.1182/blood-2009-02-207118 [PubMed]
- 25. Sakaguchi S, Ono M, Setoguchi R, Yagi H, Hori S, Fehervari Z, Shimizu J, Takahashi T, Nomura T. Foxp3+ CD25+ CD4+ natural regulatory T cells in dominant self-tolerance and autoimmune disease. Immunol Rev. 2006; 212:8–27. https://doi.org/10.1111/j.0105-2896.2006.00427.x [PubMed]
- 26. West NR, Kost SE, Martin SD, Milne K, Deleeuw RJ, Nelson BH, Watson PH. Tumour-infiltrating FOXP3(+) lymphocytes are associated with cytotoxic immune responses and good clinical outcome in oestrogen receptor-negative breast cancer. Br J Cancer. 2013; 108:155–62. https://doi.org/10.1038/bjc.2012.524 [PubMed]
- 27. Zhu G, Pei L, Yin H, Lin F, Li X, Zhu X, He W, Gou X. Profiles of tumor-infiltrating immune cells in renal cell carcinoma and their clinical implications. Oncol Lett. 2019; 18:5235–42. https://doi.org/10.3892/ol.2019.10896 [PubMed]
- 28. Marshall EA, Ng KW, Kung SH, Conway EM, Martinez VD, Halvorsen EC, Rowbotham DA, Vucic EA, Plumb AW, Becker-Santos DD, Enfield KS, Kennett JY, Bennewith KL, et al. Emerging roles of T helper 17 and regulatory T cells in lung cancer progression and metastasis. Mol Cancer. 2016; 15:67. https://doi.org/10.1186/s12943-016-0551-1 [PubMed]
- 29. Ghosh S, Taylor A, Chin M, Huang HR, Conery AR, Mertz JA, Salmeron A, Dakle PJ, Mele D, Cote A, Jayaram H, Setser JW, Poy F, et al. Regulatory T Cell Modulation by CBP/EP300 Bromodomain Inhibition. J Biol Chem. 2016; 291:13014–27. https://doi.org/10.1074/jbc.M115.708560 [PubMed]
- 30. Liu Y, Wang L, Predina J, Han R, Beier UH, Wang LC, Kapoor V, Bhatti TR, Akimova T, Singhal S, Brindle PK, Cole PA, Albelda SM, Hancock WW. Inhibition of p300 impairs Foxp3⁺ T regulatory cell function and promotes antitumor immunity. Nat Med. 2013; 19:1173–77. https://doi.org/10.1038/nm.3286 [PubMed]
- 31. Vahidi Y, Faghih Z, Talei AR, Doroudchi M, Ghaderi A. Memory CD4+ T cell subsets in tumor draining lymph nodes of breast cancer patients: A focus on T stem cell memory cells. Cell Oncol (Dordr). 2018; 41:1–11. https://doi.org/10.1007/s13402-017-0352-6 [PubMed]
- 32. Hotta K, Sho M, Fujimoto K, Shimada K, Yamato I, Anai S, Konishi N, Hirao Y, Nonomura K, Nakajima Y. Prognostic significance of CD45RO+ memory T cells in renal cell carcinoma. Br J Cancer. 2011; 105:1191–96. https://doi.org/10.1038/bjc.2011.368 [PubMed]
- 33. Bryceson YT, March ME, Ljunggren HG, Long EO. Synergy among receptors on resting NK cells for the activation of natural cytotoxicity and cytokine secretion. Blood. 2006; 107:159–66. https://doi.org/10.1182/blood-2005-04-1351 [PubMed]
- 34. Dianat-Moghadam H, Rokni M, Marofi F, Panahi Y, Yousefi M. Natural killer cell-based immunotherapy: from transplantation toward targeting cancer stem cells. J Cell Physiol. 2018; 234:259–73. https://doi.org/10.1002/jcp.26878 [PubMed]
- 35. Chalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, Schrock A, Campbell B, Shlien A, Chmielecki J, Huang F, He Y, Sun J, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med. 2017; 9:34. https://doi.org/10.1186/s13073-017-0424-2 [PubMed]
- 36. Mayakonda A, Lin DC, Assenov Y, Plass C, Koeffler HP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 2018; 28:1747–56. https://doi.org/10.1101/gr.239244.118 [PubMed]
- 37. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005; 102:15545–50. https://doi.org/10.1073/pnas.0506580102 [PubMed]
- 38. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015; 12:453–57. https://doi.org/10.1038/nmeth.3337 [PubMed]