Skip to main content

Advertisement

Tumor genomic alterations in severe-combined immunodeficiency bare-lymphocyte syndrome genes are associated with high mutational burden and disproportional neo-antigen rates

Article metrics

Abstract

The progression of cancer requires mutational adaptation to permit unrestrained proliferation. A fraction of cancer mutations are oncogenic drivers, while others are putative ‘passengers’ that do not contribute to oncogenesis. However, altered peptides arising from passenger mutations may bind MHCs and activate non-self immunologic signals (i.e. neoantigens), thus requiring immunoediting for cancer persistence. Disruption of antigen processing machinery in tumor cells may diminish this requirement. Here, we show that rare mutations in antigen processing machinery are associated with high mutational burden and increased predicted neoantigen load, providing insights into the mechanisms of high mutation burden in some patients.

Background

Cancer-expressed neoantigens are a major determinant of immunological activity and response to immunotherapy [1,2,3,4,5,6,7,8,9]. Tumor types with high tumor mutational burden (TMB), generally corresponding to high predicted neoantigen loads (e.g. POLE-mutant [10,11,12] and mismatch repair deficient/microsatellite-instable tumors; MSI) [13] demonstrate less aggressive disease courses and higher responses to immunotherapies targeting CTLA-4 and PD-1/L1 checkpoints [14,15,16]. Recently, it has been demonstrated that tumors with acquired resistance to immunotherapies can undergo immunoediting (genomic loss of neoantigens [1]) or lose key antigen presentation molecules such as beta-2-microglobulin (B2M) [17], a prerequisite for peptide-major histocompatibility complex class I (MHC-I) formation at the cell surface. Thus, immunoediting and genomic alterations in antigen presentation may work together or independently to circumvent adaptive immunity. However, other mechanisms whereby tumors may avoid T cell recognition through disruption of antigen presentation have not been thoroughly explored. Further, the process whereby tumors acquire high mutation burden while avoiding immunoediting is not well characterized in many cases.

Taking a lesson from severe combined immunodeficiency (SCID) hereditary disorders in humans, we asked whether somatic alterations in genes involved in bare lymphocyte syndrome (BLS) may be a possible mechanism for tumors to circumvent adaptive immunity and antigen presentation. BLS is a rare hereditary immunodeficiency syndrome characterized by loss-of-function germline alterations in MHC machinery (i.e. those molecules required for expression of MHC as well as peptide/antigen loading and transport) [18]. Type I BLS is caused by germline alterations in TAP1, TAP2, or TAPBP, genes encoding proteins involved in transport of peptides into the endoplasmic reticulum for loading onto MHC-I. Type II BLS results from alterations in CIITA, RFX5, RFKAP, or RFXANK, which are required for transcription of MHC-II, leading to lymphocytes (and other cell types) which lack MHC-II antigen presentation. Both Type I and Type II BLS lead to poor or nonexistent antigen presentation (via MHC-I or MHC-II) and immunodeficiency (i.e. ‘immune ignorance’) [18], resulting in susceptibility to infections and early mortality.

It is conceivable that tumor cells could somatically alter these genes in order to avoid antigen presentation synonymous to Type I or II BLS. Interestingly, somatic variants in Type-II BLS genes identified in a subset of gastric cancers were found to be associated with lack of recurrence [19]. To explore BLS genes as a target for somatic disruption of antigen presentation, we investigated four cancer types with high TMB in TCGA (colorectal, melanoma, gastric and uterine). We predicted MHC-I neoantigens using a computational approach with publicly available next generation sequencing databases. High mutational burden and increased neoantigen load are associated with rare mutations in BLS genes. In addition, immune surveillance is believed to be a selection force shaping cancer evolution, potentially acting to remove neoantigen-carrying clones by purifying (negative) selection. Consistently, reduced purifying selection on putative neoantigen-carrying genes are detected in patients with Type I BLS gene alterations, suggesting loss of immunoediting due to disrupted antigen presentation. Our findings show an insight into the mechanisms of high mutation burden, and potentially a loss of immunoediting in patients with somatic BLS alterations.

Results and discussion

We hypothesized that somatic alteration of BLS genes in tumors could abrogate the necessity of immuno-editing by reducing antigen presentation. To test this hypothesis, we determined the rate of mutations in Type I and Type II BLS genes across human tumor datasets in cBioPortal [20]. We identified sporadic mutations in BLS genes (MHC-I and MHC-II) across cancer types, some of which were recurrent (Additional file 1: Figure S1). We specifically focused on tumor types with single-nucleotide variants and indels. Interestingly, we noted a high preponderance of truncating and early-stop codon alterations in BLS genes.

Next, four TCGA tumor databases with a sizeable population of BLS-altered tumors (colorectal, melanoma, gastric and uterine) were selected. Intriguingly, within each database, tumors harboring BLS mutations were strongly associated with a higher TMB (p = 1.42e-17 to 8.82e-75 for all comparisons; Fig. 1a) compared to non-BLS-altered, suggesting that putative disruption of MHC-I or MHC-II antigen presentation may eliminate or reduce the requirement for immune-editing of the tumor by permiting the generation of neoantigens without an immune-mediated consequence.

Fig. 1
figure1

BLS-altered primary human tumors demonstrate high TMB in the context of POLE-mutator phenotypes. a) Left panels demonstrate association of BLS-altered tumors for each TCGA dataset (colorectal, melanoma, gastric/stomach, and uterine). BLS-alterations considered were nonsynonymous variants and deep deletions. The types of mutations are shown in the oncoprint on the right panels for each tumor type: light blue (homologous deletion), green (missense mutation), black (multi hit), blue (frame shift deletion), purple (frame shift insertion), brown (in frame deletion), and red (nonsense mutation). b) Association of TMB with the presence of a BLS class I/II alteration, a POLE mutation, or the combination of both. All comparisons were highly significant by negative binomial model test

Recently, proofreading domain mutations in polymerase-epsilon (Polε/POLE) have been described as associated with high TMB, immunogenic phenotypes, and response to immunotherapies [11, 12, 21]. Consistent with our hypothesis, BLS-altered tumors acquire greatly-increased TMB in the context of a POLE-deficient ‘mutator phenotype’ [21] (Fig. 1b). In all 4 tumor types tested, both POLE-mutant and BLS-mutant tumors had elevated TMBs. However, tumors with mutations in both POLE and BLS genes had exceptionally high mutation rates (p ≤ 0.01 for all tumor types compared to POLE-mutant alone), consistent with the notion that within a mutator-phenotype background, disruption of BLS genes may permit the acquisition of higher TMB without necessity of immunoediting.

Not all mutations bind efficiently to patient-specific MHCs, and thus computational methods based on HLA-haplotype and known binding affinities are required to improve prediction accuracy for putative neoantigens. By in silico prediction of MHC class I specific neoantigens (which is more precise than MHC-II epitope prediction) using transcriptomic and whole-exome next generation sequencing data from TCGA, we asked whether BLS type I altered tumors within this cohort demonstrated a greater degree of predicted neoantigens than non-BLS type I altered counterparts (Additional file 2: Data Set). For both TMB and at each cut off of predicted binding affinity for neoantigens (50-500 nM) (Additional file 3: Table S1), we observed a strong association with BLS type I altered tumor status in all four cancer types (Fig. 2a).

Fig. 2
figure2

BLS altered tumors demonstrate disproportionally high neoantigen rates consistent with lack of immunoediting. a) TMB and neoantigen load across predicted binding affinity cutoffs in 4 cancer types. Type I BLS altered tumors demonstrated higher TMB and higher neoantigen rates across all cutoffs applied: 50 nM (NEO50), 100 (NEO100), 250 nM (NEO250), 500 nM (NEO500). b) Proportion of neoantigens carrying genes with ω (dN/dS ratio) < 1 between Type I BLS altered tumors and the Normal samples. dN/dS based on missense mutations and nonsense mutations were separately calculated and shown. Error bars show the 95% confidence interval

Finally, we tested whether BLS type I altered tumors disproportionally develop putative neoantigens with reduced purifying selection, potentially due to disrupted immunoediting. Using a specifically designed method adapted from molecular evolution theory for cancer somatic evolution, we determined the selection pressure of neoantigen-carrying genes in different cohorts by calculating dN/dS ratios, a classic indicator for positive or purifying selection detection [22] (Additional file 3: Table S2). A consistent pattern of positive selection in majority of the neoantigen carrying genes were observed in all 4 cancer types as reported previously for the somatic mutations in all TCGA projects [22]. Interestingly, among the significantly, positively selected genes, a Type I BLS gene -TAP2, which encodes an ABC-transport protein required for transport of class-I peptides into the endoplasmic reticulum [23], was identified as a particular driver in colorectal cancer (q = 8.6e-4). In general, purifying selection showed stronger pressure on neoantigen carrying genes in non-BLS-I altered samples comparing with the BLS-I altered samples (p values of proportion test for percentages of neoantigen carrying genes with ω < 1 in all between groups are all under 0.05; Fig. 2b and Additional file 1: Figure S2). Thus, BLS-I altered samples appear to have reduced purifying selection pressure for developing neoantigens, possibly due to reduced immunoediting restrictions.

The effects of BLS alterations on clonal mutation selection and therapeutic response to immunotherapies are a potential area of importance to study. There have been several studies where data on somatic mutation load/neoantigen load have been collected on/post immunotherapy in patients [24,25,26]. We identified the cases with BLS I/II mutations in post therapy tumors, and compared the mutation load or neoantigen load, if available, with the cases without BLS mutations.

In the Riaz et al. cohort [25], which examined tumor mutational evolution during anti-PD-1 therapy, only 3 patients carried BLS gene mutations in their on-therapy tumors, two of which carried BLS type-II mutations in CIITA, and responded to anti-PD-1 therapy. Since the neoantigens predicted are MHC-I related, only one patient (who experienced stable disease [SD]) with a BLS type-I mutation in TAP2 gene was appropriate to evaluate the change, which was not substantially different compared with all SD patients without BLS mutations.

For the Roh et al. cohort [21], exome sequencing data from longitudinal tissue samples from metastatic melanoma patients treated with sequential immune checkpoint blockade (CTLA-4 blockade followed by PD-1 blockade at time of progression) were available. As shown and described in Additional file 1: Figure S3A, this dataset contained two BLS type-I altered patients and five BLS type-II altered patients. Based on the limited information of the BLS altered cases in this dataset, there is a trend that BLS-altered patients retained high mutation load during or post immunotherapy, particularly in non-responders. However, these are highly exploratory data and conclusions cannot be drawn.

Finally, in the Le et al. cohort [15], mutation profiles and clinical information of 30 patients with mCRC treated with anti-PD-1 therapy, including a substantial proportion that were mismatch repair-deficient/microsatellite-instable, were investigated. There were 7 cases of BLS-altered tumors, which demonstrated higher mutation load, consistent with our prior data (Additional file 1: Figure S3B). However, all patients with BLS mutations were also microsatellite-instable in this dataset, which have an exceptionally high rate of response to immunotherapy [13, 26] (Additional file 1: Figure S3C). As such, 6 out of 7 BLS-altered cases were responders, complicating purifying selection analysis, as tumor response is typically accompanied by a substantial contraction of the tumor mutation burden due to elimination of tumor cells. No differences in outcomes were observed in this small population (data not shown). Since MSI is a major factor influencing genome stability, dysfunction of antigen presentation may be muted by the diversity and frequency of mutations generated in MSI tumors. This could lead to field-effects on anti-tumor immunity and immunoediting, or a substantial response to MHC-I restricted antigens (in BLS type-II-altered tumors, or vice versa). In brief, MSI status seems to predominate over mutation load in this case and may overcome the influences of BLS gene alterations, but it is difficult to conclude this without a larger sample size.

Conclusions

In conclusion, we provide in silico support for the hypothesis that rare functionally-disruptive alterations in BLS genes may permit acquisition of high tumor mutation burdens, particularly in the context of mutation-prone genetic backgrounds such as POLE mutants. Furthermore, after adjusting for mutation burden, there was a significant gain in predicted neoantigens arising from BLS-altered tumors versus other tumors due to weakened constraints of purifying selection. These data support further investigation into the functionality of BLS alterations in cancer both molecularly and in clinical outcomes. Importantly, however, such investigations must include a detailed functional analysis, particularly in cases where the effect of the somatic alteration is unclear, such as SNVs.

Whether or not defects in BLS-associated genes are sufficient to drive immune evasion and their influence on immunotherapy outcomes remains to be more thoroughly explored. It is currently unclear in such cases whether activation of inflammation pathways (e.g. through Type I or II interferon activation/stimulation, such as STING activation or others) could overcome somatic defects in antigen presentation. The answers to such questions will likely need to be explored in murine models and may be subject to the type of alteration and clonality within the tumor cell population. Furthermore, clonal evolution studies of metastatic tumors harboring BLS-like alterations could offer additional insights into their role in immunoediting and driving immune evasion in certain contexts. Future work should extend these findings toward models testing causative association (e.g. mutator-phenotype murine cell lines harboring defects in BLS genes and testing acquisition of tumor mutations over time), as well as their impact on immunotherapy outcomes in larger analyses, which has already been demonstrated with B2M loss in some longitudinal studies [17, 26].

Methods

Genomic databases

Genomic variants and gene expression data, clinical factors (e.g. MSI status in colorectal tumors) and total mutational burden were accessed via the GDC Data Portal (https://portal.gdc.cancer.gov/). Copy number variation information was accessed from the Broad Institute TCGA GDAC (http://firebrowse.org/, GISTIC2 version, [27,28,29,30]). Data from immunotherapy-treated patients were accessed through their associated supplementary data files [24, 25].

Neoantigen prediction

Patients specific HLA class I genotype was predicted using the whole exome sequencing data from matched normal tissue with OptiType [31]. GDC VCF files from the Mutect2 pipeline and the reads count based expression profiles for each case were achieved. Only non-synonymous SNPs or Indels were kept for somatic mutation load calculation and for neoantigen prediction. Raw reads count was transformed to transcript per million (TPM) as described in Bo et al. [32]. Combined with patient specific HLA genotypes, non-synonymous mutation list, and gene expression profile, MHC class I specific binding neoantigens were predicted in MuPeXI [33]. Binding affinity of MHC class I with candidate peptides was evaluated with netMHCpan 4.0 [34]. Mutated peptides with IC50 score < 500 nm while the according normal peptides with IC50 score > 500 nm were considered as candidate neoantigen peptides [35].

Screen for selection signal

We used dNdScv, a set of maximum-likelihood dN/dS methods designed to quantify selection in cancer and somatic evolution [22], to study the selection pressure of the neoantigens in our study. dN/dS is the normalized ratio of non-synonymous to synonymous substitutions and widely used for selection detection in molecular evolution theory [36]. dNdScv utilized a modified dN/dS measurement, as described in Martincorena et al. [22], for the analysis of somatic mutations in cancer genomes with a few critical refinements, including a more comprehensive substitution model by considering the context-dependent mutational processes, addition of other non-synonymous mutation types beyond missense, careful elimination of germline variant contamination, and consideration of the variable mutation rate along the genome. This approach has been successfully applied for genome-wide scan of the somatic mutation selection patterns of different cancer types [22]. In Fig. 2b, samples were separated into “BLS I only” group if any non-synonymous mutations or deep deletions were occurred on TAP1, TAP2 or TAPBP gene and “Others” if no non-synomymous mutations were detected in any BLS genes. For colorectal, gastric and uterine cancer, default parameters were used in dNdScv, while for melanoma, ultra-hypermutator mode was switched on due to large proportion of hypermutator tumors.

Statistics

All statistical analyses were performed in R (www.r-project.org). For mutation load and neoantigen load comparison, negative binomial model was used to test the differences between different sample groups. A two-sample proportion test was used to test the difference of ω (dN/dS ratio) < 1 between BLS I only group and Others group and genes with ω = 0 were excluded. For all comparisons, p value was shown in plots. Bar graphs show mean ± standard error of the mean (SEM), unless otherwise stated in the figure legend.

References

  1. 1.

    Anagnostou V, Smith KN, Forde PM, Niknafs N, Bhattacharya R, White J, et al. Evolution of Neoantigen landscape during immune checkpoint blockade in non-small cell lung Cancer. Cancer Discov. 2017;7(3):264–76.

  2. 2.

    Desrichard A, Snyder A, Chan TA. Cancer Neoantigens and applications for immunotherapy. Clin Cancer Res. 2016;22(4):807–12.

  3. 3.

    Gubin MM, Artyomov MN, Mardis ER, Schreiber RD. Tumor neoantigens: building a framework for personalized cancer immunotherapy. J Clin Invest. 2015;125(9):3413–21.

  4. 4.

    Lu YC, Robbins PF. Cancer immunotherapy targeting neoantigens. Semin Immunol. 2016;28(1):22–7.

  5. 5.

    Lu YC, Robbins PF. Targeting neoantigens for cancer immunotherapy. Int Immunol. 2016;28(7):365–70.

  6. 6.

    McGranahan N, Furness AJ, Rosenthal R, Ramskov S, Lyngaa R, Saini SK, et al. Clonal neoantigens elicit T cell immunoreactivity and sensitivity to immune checkpoint blockade. Science. 2016;351(6280):1463–9.

  7. 7.

    Rooney MS, Shukla SA, Wu CJ, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. 2015;160(1–2):48–61.

  8. 8.

    Snyder A, Makarov V, Merghoub T, Yuan J, Zaretsky JM, Desrichard A, et al. Genetic basis for clinical response to CTLA-4 blockade in melanoma. N Engl J Med. 2014;371(23):2189–99.

  9. 9.

    Van Allen EM, Miao D, Schilling B, Shukla SA, Blank C, Zimmer L, et al. Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science. 2015;350(6257):207–11.

  10. 10.

    Gong J, Wang C, Lee PP, Chu P, Fakih M. Response to PD-1 blockade in microsatellite stable metastatic colorectal Cancer harboring a POLE mutation. J Natl Compr Cancer Netw. 2017;15(2):142–7.

  11. 11.

    Howitt BE, Shukla SA, Sholl LM, Ritterhouse LL, Watkins JC, Rodig S, et al. Association of Polymerase e-mutated and microsatellite-instable endometrial cancers with Neoantigen load, number of tumor-infiltrating lymphocytes, and expression of PD-1 and PD-L1. JAMA Oncol. 2015;1(9):1319–23.

  12. 12.

    Mehnert JM, Panda A, Zhong H, Hirshfield K, Damare S, Lane K, et al. Immune activation and response to pembrolizumab in POLE-mutant endometrial cancer. J Clin Invest. 2016;126(6):2334–40.

  13. 13.

    Le DT, Uram JN, Wang H, Bartlett BR, Kemberling H, Eyring AD, et al. PD-1 blockade in tumors with mismatch-repair deficiency. N Engl J Med. 2015;372(26):2509–20.

  14. 14.

    Johnson DB, Frampton GM, Rioth MJ, Yusko E, Xu Y, Guo X, et al. Targeted next generation sequencing identifies markers of response to PD-1 blockade. Cancer Immunol Res. 2016;4(11):959–67.

  15. 15.

    Rizvi NA, Hellmann MD, Snyder A, Kvistborg P, Makarov V, Havel JJ, et al. Cancer immunologyMutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science. 2015;348(6230):124–8.

  16. 16.

    Chalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome medicine. 2017;9(1):34.

  17. 17.

    Zaretsky JM, Garcia-Diaz A, Shin DS, Escuin-Ordinas H, Hugo W, Hu-Lieskovan S, et al. Mutations associated with acquired resistance to PD-1 blockade in melanoma. N Engl J Med. 2016;375(9):819–29.

  18. 18.

    Reith W, Mach B. The bare lymphocyte syndrome and the regulation of MHC expression. Annu Rev Immunol. 2001;19:331–73.

  19. 19.

    Yavorski JM, Blanck G. MHC class II associated stomach cancer mutations correlate with lack of subsequent tumor development. Molecular and clinical oncology. 2017;7(6):1119–21.

  20. 20.

    Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.

  21. 21.

    Heitzer E, Tomlinson I. Replicative DNA polymerase mutations in cancer. Curr Opin Genet Dev. 2014;24:107–13.

  22. 22.

    Martincorena I, Raine KM, Gerstung M, Dawson KJ, Haase K, Van Loo P, et al. Universal patterns of selection in Cancer and somatic tissues. Cell. 2017;171(5):1029–41 e21.

  23. 23.

    Townsend A, Trowsdale J. The transporters associated with antigen presentation. Semin Cell Biol. 1993;4(1):53–61.

  24. 24.

    Roh W, Chen PL, Reuben A, Spencer CN, Prieto PA, Miller JP, et al. Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance. Sci Transl Med. 2017;9(379).

  25. 25.

    Riaz N, Havel JJ, Makarov V, Desrichard A, Urba WJ, Sims JS, et al. Tumor and microenvironment evolution during immunotherapy with Nivolumab. Cell. 2017;171(4):934–49 e16.

  26. 26.

    Le DT, Durham JN, Smith KN, Wang H, Bartlett BR, Aulakh LK, et al. Mismatch repair deficiency predicts response of solid tumors to PD-1 blockade. Science. 2017;357(6349):409–13.

  27. 27.

    Beroukhim R, Getz G, Nghiemphu L, Barretina J, Hsueh T, Linhart D, et al. Assessing the significance of chromosomal aberrations in cancer: methodology and application to glioma. Proc Natl Acad Sci U S A. 2007;104(50):20007–12.

  28. 28.

    Beroukhim R, Mermel CH, Porter D, Wei G, Raychaudhuri S, Donovan J, et al. The landscape of somatic copy-number alteration across human cancers. Nature. 2010;463(7283):899–905.

  29. 29.

    McCarroll SA, Kuruvilla FG, Korn JM, Cawley S, Nemesh J, Wysoker A, et al. Integrated detection and population-genetic analysis of SNPs and copy number variation. Nat Genet. 2008;40(10):1166–74.

  30. 30.

    Mermel CH, Schumacher SE, Hill B, Meyerson ML, Beroukhim R, Getz G. GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol. 2011;12(4):R41.

  31. 31.

    Szolek A, Schubert B, Mohr C, Sturm M, Feldhahn M, Kohlbacher O. OptiType: precision HLA typing from next-generation sequencing data. Bioinformatics. 2014;30(23):3310–6.

  32. 32.

    Li B, Ruotti V, Stewart RM, Thomson JA, Dewey CN. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics. 2010;26(4):493–500.

  33. 33.

    Bjerregaard AM, Nielsen M, Hadrup SR, Szallasi Z, Eklund AC. MuPeXI: prediction of neo-epitopes from tumor sequencing data. Cancer Immunol Immunother. 2017;66(9):1123–30.

  34. 34.

    Jurtz V, Paul S, Andreatta M, Marcatili P, Peters B, Nielsen M. NetMHCpan-4.0: improved peptide-MHC class I interaction predictions integrating eluted ligand and peptide binding affinity data. J Immunol. 2017;199(9):3360–8.

  35. 35.

    Rajasagi M, Shukla SA, Fritsch EF, Keskin DB, DeLuca D, Carmona E, et al. Systematic identification of personal tumor-specific neoantigens in chronic lymphocytic leukemia. Blood. 2014;124(3):453–62.

  36. 36.

    Yang Z, Nielsen R. Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol. 2000;17(1):32–43.

Download references

Acknowledgements

Not applicable.

Funding

Funding for this work was provided by Susan G. Komen Career Catalyst Grant CCR14299052 (JMB), NIH/NCI R00CA181491 (JMB), and the Vanderbilt-Ingram Cancer Center Support Grant P30 CA68485.

Availability of data and materials

The datasets analyzed in the current study are available in Genomic Data Commons Data Portal and associated raw data repositories.

Author information

YW and YX performed data analysis, biostatistics, and contributed to study design. DBJ, YX, and JMB generated study concept and oversight. SL and LAD provided data and analyses of MSI tumors. All authors read and approved the final manuscript.

Correspondence to Yaomin Xu or Justin M. Balko.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests relevant to this work. DBJ is on advisory boards for Array, BMS, Genoptix, Incyte, Merck, and has received research funding from BMS and Incyte. JMB receives consulting fees from Novartis and grant funding from Incyte and Genentech. Luis A. Diaz (LAD) is a member of the board of directors of Personal Genome Diagnostics (PGDx) and Jounce Therapeutics. LAD holds equity in PapGene, PGDx and Neophore. He is a paid consultant for PGDx and Phoremost. He is an uncompensated consultant for Merck. LAD is an inventor of multiple patents related to technology for circulating tumor DNA analyses and mismatch repair deficiency for diagnosis and therapy (WO2016077553A1) from Johns Hopkins University. These licenses and relationships are associated with equity or royalty payments to LAD. The terms of all these arrangements are being managed by Johns Hopkins and Memorial Sloan Kettering in accordance with their conflict of interest policies. In addition, in the past 5 years, LAD has participated as a paid consultant for Merck and for one-time engagements with Caris, Lyndra, Genocea Biosciences, Illumina and Cell Design Labs.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Rates and types of BLS alterations across human tumor datasets. A) cBioportal was utilized to access rates of deep deletion (GISTIC) and nonsynonymous genomic alterations in BLS genes across all datasets. B) Lollipop plot of individual mutations across human tumor datasets arranged by gene. Figure S2. ω (dN/dS ratio) of neoantigens carrying genes between Type I BLS altered tumors and the Normal samples. dN/dS based on missense mutations and nonsense mutations were separately calculated and shown. Figure S3. Changes in mutation burden and neoantigen load with immunotherapy. A) Sequencing data from longitudinal tissue samples from metastatic melanoma patients treated with sequential immune checkpoint blockade (CTLA-4 blockade followed by PD-1 blockade at time of progression) were extracted [24]. As shown in, there are two BLS I altered patients (Pt16: post-CTLA4_pre-PD1 and on-PD1; Pt22: on-CTLA4) and five BLS type-II altered patients (Pt21: post-CTLA4_pre-PD1; Pt26: post-CTLA4_pre-PD1 and post-PD1; Pt47: pre-CTLA4; Pt49: post-PD1; Pt53: pre-CTLA4 and on-PD1). Pt16 and Pt26 demonstrated relatively high and stable mutation load during anti-PD1 therapy. Pt22 had no BLS mutations before therapy but gained a new BLS type-I mutation during anti-CTLA4 treatment. Since this is a responder, total mutation burden was reduced dramatically, as expected. Only a single sample was available for both Pt21 and Pt49, however, the mutation load was relatively high in both cases. Interestingly, for Pt47, subclones harboring a BLS type-II mutation were eliminated during immunotherapy, accompanying the dramatic reduction in mutation load. Mutation load was slightly increased during the series of treatments for Pt53. B) Comparison of nonsynonymous mutation load between BLS mutated patients and control group in Le et al. 2017 Cohort. Only 30 mCRC patients with tumors post anti-PD1 therapy were investigated. BLS-altered tumors show significant or near significant higher mutation load than normal tumors, based on Wilcoxon test (p values are shown in plot). and C) Strata by MSI/MSS, the higher mutation load is associated with MSI. There are no significant differences between BLS-altered tumors and normal ones in MSI. Since all BLS-altered tumors are all MSI, we are unable to investigate how BLS mutations influence the mutation load without the interference of MSI. (PDF 674 kb)

Additional file 2:

Data Set 1. Neoantigens predicted based on patient-specific HLA genotypes in colorectal, melanoma, gastric, and uterine cancers. (ZIP 63724 kb)

Additional file 3:

Table S1. Mutation profiles of BLS genes in colorectal, melanoma, gastric, and uterine cancers. Table S2. ω (dN/dS ratio) and neutrality tests for selection signal of neoantigen carrying genes in BLS I only and Others group. Syn: synonymous mutation; mis: missense mutation; non: nonsense mutation. w: ω (dN/dS ratio); p: P-values for substitutions are obtained by Likelihood-Ratio Tests as described in (Martincorena et al., 2017); q: q-values are obtained by Benjamini-Hodgberg’s multiple testing correction. (XLS 12882 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, Y., Johnson, D.B., Lu, S. et al. Tumor genomic alterations in severe-combined immunodeficiency bare-lymphocyte syndrome genes are associated with high mutational burden and disproportional neo-antigen rates. j. immunotherapy cancer 7, 123 (2019) doi:10.1186/s40425-019-0584-2

Download citation

Keywords

  • Bare lymphocyte syndrome
  • Severe combined immunodeficiency syndrome
  • Antigen presentation
  • Neoantigens
  • Cancer
  • Immunotherapy