Tumor mutational burden quantification from targeted gene panels: major advancements and challenges
===================================================================================================

* Laura Fancello
* Sara Gandini
* Pier Giuseppe Pelicci
* Luca Mazzarella

## Abstract

Tumor mutational burden (TMB), the total number of somatic coding mutations in a tumor, is emerging as a promising biomarker for immunotherapy response in cancer patients. TMB can be quantitated by a number of NGS-based sequencing technologies. Whole Exome Sequencing (WES) allows comprehensive measurement of TMB and is considered the gold standard. However, to date WES remains confined to research settings, due to high cost of the large genomic space sequenced. In the clinical setting, instead, targeted enrichment panels (gene panels) of various genomic sizes are emerging as the routine technology for TMB assessment. This stimulated the development of various methods for panel-based TMB quantification, and prompted the multiplication of studies assessing whether TMB can be confidently estimated from the smaller genomic space sampled by gene panels. In this review, we inventory the collection of available gene panels tested for this purpose, illustrating their technical specifications and describing their accuracy and clinical value in TMB assessment. Moreover, we highlight how various experimental, platform-related or methodological variables, as well as bioinformatic pipelines, influence panel-based TMB quantification. The lack of harmonization in panel-based TMB quantification, of adequate methods to convert TMB estimates across different panels and of robust predictive cutoffs, currently represents one of the main limitations to adopt TMB as a biomarker in clinical practice. This overview on the heterogeneous landscape of panel-based TMB quantification aims at providing a context to discuss common standards and illustrates the strong need of further validation and consolidation studies for the clinical interpretation of panel-based TMB values.

*   Tumor mutational burden
*   TMB
*   Gene panels
*   Targeted enrichment sequencing
*   Immunotherapy

## Tumor mutational burden: an emerging biomarker for cancer immunotherapy

Immunotherapy with immune checkpoint inhibitors targeting cytotoxic T lymphocyte associated 4 (CTLA-4) or programmed cell death 1 (PD-1) or its ligand (PD-L1) can provide important clinical benefit to patients affected by multiple cancers, most notably lung cancer [1, 2], melanoma [3], renal cancer [4] and urothelial carcinoma [5]. However, only a fraction of patients currently treated by immune checkpoint inhibitors derive benefit from it, while a minority of them suffers from severe side effects. Given the significant cost and non-negligible toxicity of these therapies, the identification of strategies to adequately select those patients most likely to show a favorable response is recognized as an urgent medical need. A few potential biomarkers have been identified up to now, such as PD-L1 gene expression, microsatellite instability (MSI), mismatch repair deficiency (dMMR), POLE or JAK1/2 mutations, immune cell infiltration, IFNγ expression, tumor mutational burden (TMB) or neoantigen burden [6, 7].

TMB is a measure of the total amount of somatic coding mutations in a tumor and it is currently investigated as a potential biomarker in non-small cell lung carcinoma (NSCLC) [8–10]. Accumulating evidence, however, suggests its potential usefulness also in melanoma [8, 11–14], urothelial cancer [5, 15, 16], mismatch-repair deficient colorectal tumors [17] and other cancer types [18]. Its pattern and distribution is highly variable across different cancer types, with over 1000-fold difference between cancer types with the lowest mutational burden and those with the highest mutational burden, such as those associated with DNA environmental damage, i.e. by exposure to tobacco smoke or UVs [19, 20]. Increased TMB was also observed in tumors with defects in DNA mismatch repair and DNA replication or in tumors characterized by microsatellite instability, as in colorectal cancer [21, 22]. Highly mutated tumors are more likely to produce abundance of tumor-specific mutant epitopes, which may function as neoantigens recognized as non-self by the immune system. Therefore, increased activation of immune cells by treatment with immune checkpoint inhibitors may lead to improved immune-mediated tumor-cell clearance and clinical response in these tumors (Fig. 1). A significant association between neoantigen production and immune-mediated clinical response was indeed observed in several studies [9, 11, 14, 23]. Measurement of this neoantigen production, though, is expensive and time-consuming. Tumor neoantigens can be generated by mutations or by gene fusions, especially out-of-frame fusions. Although some pipelines have recently been developed for the identification of neoantigens derived from gene fusions [24], most research up to now has estimated overall neoantigen load based only on somatic nonsynonymous coding mutations, called by Whole Exome Sequencing (WES). Briefly, somatic nonsynonymous coding mutations are identified by WES and, if RNA sequencing is also available, only mutations occurring in expressed genes are retained. Peptides containing selected mutations are then identified in silico and the efficiency of their presentation to the immune system may be evaluated by mass spectrometry or by algorithms that consider their predicted affinity to the MHC class I complex and patient-specific HLA class I alleles [14, 25]. In comparison with overall neoantigen load, TMB is easier to measure and correlates with it. Although not all mutations can give rise to tumor immunogenic peptides, their number influences the amount of neoantigens potentially produced. High TMB correlates with long-term clinical benefit from immune checkpoint inhibitors in patients with melanoma [14], NSCLC [9, 26–28] and urothelial carcinoma [5, 15, 16, 29]. In addition to that, patients with mismatch repair (MMR) deficient tumors are more responsive to immunotherapy, probably due to their high tumor mutational burden [17]. Therefore, although not always capable to explain the clinical benefit alone, TMB is a good approximation for neoantigen load assessment [14], is technically less challenging and less expensive and may represent a better suited predictive biomarker for immunotherapy response. 

![Fig. 1](http://jitc.bmj.com/https://jitc.bmj.com/content/jitc/7/1/183/F1.medium.gif)

[Fig. 1](http://jitc.bmj.com/content/7/1/183/F1)

Fig. 1 
Tumor mutational burden as immunotherapy biomarker. Interaction between tumor mutational burden, neoantigen production and immune checkpoints. Hyper-mutated tumors (bottom) are more likely than hypo-mutated tumors (top) to generate tumor-specific peptides (neoantigens) recognized by the immune system. However, immune surveillance can be restrained by simultaneous high expression of PD-L1, which delivers a suppressive signal to T cells. PD-L1/PD-1 interaction and other immune checkpoints can be inhibited by immune checkpoint inhibitors, restoring immune response

TMB may also represent a relevant prognostic biomarker. In BRCA-1/2 mutated ovarian cancers, TMB correlates with improved overall survival [30, 31]. In breast cancer patients, tumors with high TMB and favorable immune-infiltrate (“hot tumors”) are associated with prolonged survival [32]. Consistently, basal cell carcinoma, which is characterized by very high TMB, presents with slow growth rates and rare metastases. Although not definitively demonstrated, we can speculate that this less aggressive phenotype may be due to hypermutation, which would trigger, via increased neoantigen production, a more effective immune response of the host [33].

## Quantification of tumor mutational burden from gene panels: “yes we can”

Initial studies showing a correlation between TMB and enhanced response to immunotherapy were based on WES datasets for TMB quantification [9, 14, 17]. WES allows a direct measurement of TMB, yet it remains unsuitable as routine technology in clinical practice, because expensive, labor-intensive and time-consuming. Therefore, several studies explored the possibility to provide equally accurate and clinically predictive TMB estimates from targeted enrichment sequencing, using various gene panels (Table 1, Additional file 1: Table S1). 

View this table:
[Table1](http://jitc.bmj.com/content/7/1/183/T1)

Overview of the main published studies on TMB quantification from gene panels

The main challenge for accurate panel-based TMB quantification is the ability to extrapolate the global mutational burden from the narrow sequencing space targeted by a gene panel. In silico analyses were performed to test the concordance between panel-based and WES-based TMB, which is considered the reference for TMB quantification. Publicly available WES datasets were downsampled to the subset of genes targeted in the panel under consideration and TMB values from such simulated gene panels were compared with TMB values from the original WES (Additional file 7: Figure S1), finding high correlation between the two (Additional file 2: Table S2, Additional file 8: Figure S2). Most of these in silico analyses were performed using publicly available WES datasets from TCGA, with the exception of the Oncomine Tumor Mutation Load Assay or NovoPM and CANCERPLEX gene panels, for which WES datasets from COSMIC or from other sources were used. Regardless, similar correlation values were reported for the different gene panels tested (Additional file 2: Table S2, Additional file 8: Figure S2). For some of these gene panels (FoundationOne, Trusight170, Oncomine Tumor Mutation Load Assay, Oncomine Comprehensive Assay V3 and MSK-IMPACT gene panels), an empirical approach was also used to test the concordance between panel-based and WES-based TMB quantification, based on matched sequencing by gene panel and WES of the same tumor sample and comparison of matched TMB values (Additional file 3: Table S3, Additional file 9: Figure S3). Accuracy of panel-based TMB quantification is influenced by statistical sampling effects and small panels provide less precise TMB estimates [22, 34–36]. It was demonstrated that TMB values from the FoundationOne gene panel, which targets 1.1 Mb of genomic space, are similar to those from WES, whereas accuracy drops importantly when sequencing less than 0.5 Mb [22]. Another study simulated sequencing of theoretical gene panels of different sizes and identified 1.5 to 3 Mb as the best suited targeted genomic size to confidently estimate TMB [35]. Moreover, the deviation between WES- and panel-based TMB appears more relevant for samples with low to moderate underlying TMB levels, compared to samples with high TMB [22, 35, 36]. Another retrospective study on a commercial panel of 248 genes likewise cautions against small gene panels which would lead to TMB overestimation [37].

Besides the accuracy of panel-based TMB quantification, it is critical to assess its capability to discriminate between immunotherapy responders and non-responders, as previously observed for WES-based TMB. Several exploratory analyses demonstrated that panel-based TMB, as simulated in silico by downsampling a WES dataset to only include genes targeted by the FoundationOne gene panel, associates with immunotherapy response [8, 26] or with signatures of immune checkpoint activation [38]. Comparable results were observed in similar in silico analyses for other gene panels, such as the Trusight170 [39, 40] or MSK-IMPACT [26] (Additional file 4: Table S4). Notably, direct measurement of TMB from the Oncomine Tumor Load Assay shows that this panel-based TMB value allows to classify colorectal cancer cases based on their MSI status [39, 41]. Since in this cancer type MSI positively correlates with immunotherapy response, this is a further, yet indirect evidence, of the capability to predict immunotherapy response, using a panel-based TMB estimate. Most importantly, a few clinical studies demonstrated that TMB directly estimated using gene panels is higher in those patients who benefit more from immune checkpoint blockade treatment, thus providing “real-life” evidence for its potential clinical predictive value (Fig. 2, Additional file 5: Table S5). A direct association with immunotherapy response was shown for the MSK-IMPACT [42, 43] and the Guardant360 gene panels [44] but most of the reported studies utilized the FoundationOne gene panel (Fig. 2, Additional file 5: Table S5). In particular, in the CheckMate 227 trial, NSCLC patients with high TMB (> 10 mutations per Mb, measured by FoundationOne) presented increased progression-free survival after combined anti-CTLA-4 plus anti-PD-1/PD-L1 therapy [45]. Interestingly, TMB was predictive of anti-PD-L1 monotherapy response in NSCLC (POPLAR trial, [27]) and metastatic urothelial carcinoma patients [5, 15, 16], independently from the PD-L1 expression status. Analysis of archival tumor samples from melanoma patients treated by anti-PD-1/PD-L1 monotherapy also showed superior response rates, progression-free survival and overall survival in high TMB cases [12]. Moreover, a retrospective study on 151 patients across diverse tumor types showed that cancer patients with higher TMB, benefit more from anti-PD-1/PD-L1, anti-CTLA-4 or high dose IL2 monotherapy [18]. The same was not observed for combined anti-PD-1/PD-L1 plus anti-CTLA-4 therapy but the available number of samples may be too small to draw conclusions [18]. 

![Fig. 2](http://jitc.bmj.com/https://jitc.bmj.com/content/jitc/7/1/183/F2.medium.gif)

[Fig. 2](http://jitc.bmj.com/content/7/1/183/F2)

Fig. 2 
TMB association with progression-free survival. Forest plot of hazard ratios (HR) comparing progression-free survival (PFS) between patients with high or low TMB, as indicated in the “Comparison” column. If not specified otherwise, TMB is reported as number of mutations per Mb. All patients were treated with immune checkpoint inhibitors (*ICI*). Bars represent the 95% confidence intervals. Size of the box is proportional to precision. Reference to the study and the analyzed cancer type are also reported together with the log-rank *p*-value. *Q1-Q4*: quartiles; *VUS*: variants of unknown significance. ***: TMB quantified from blood; ****: Cox proportional hazards model adjusted for age, gender, disease stage and prior therapy by ipilimumab

## Need for standardization of TMB quantification and reporting

Despite the increasing number of studies showing the potential clinical relevance of panel-based TMB as a predictive biomarker for immunotherapy response, its use in the clinical setting is currently limited by the absence of standard methods of quantification and the lack of a robust and universal cutoff to identify immunotherapy responders.

Panel-based TMB quantification is influenced by various experimental factors affecting library construction and sequencing, by the pipeline used to call mutations and by the capability to extrapolate TMB values from the restricted genomic space sampled by gene panels to the whole genome (Fig. 3a). Experimental factors (e.g. tumor purity or sequencing depth) and the variant calling pipeline (e.g. the variant calling algorithm and the method to remove germline variants) can significantly affect the number of called somatic mutations and have a similar impact on both panel-based and WES-based TMB quantification. Indeed, the adoption of a well-documented standard pipeline was already claimed for WES analyses as an urgent need to allow data interoperability between different platforms [46]. The same applies to panel sequencing for TMB quantification. In this context, an important factor investigated for its influence on the number of called somatic variants is the method chosen to identify and remove germline variants. Indeed, since only somatic mutations can potentially produce tumor neoantigens recognized as non-self by the immune system, it is important to remove germline variants in TMB quantification. It was observed that the use of an in silico method for somatic variant calling instead of matched tumor-normal samples, leads to increased false positive somatic variants, which has an important influence on the accuracy of TMB quantification, especially for small gene panels [34]. To avoid this, it was proposed to perform TMB quantification using only high-confidence regions [47] (e.g. regions of the genome, devoid of potential systematic biases or structural variants, where mutations can be confidently called), as defined by Zook et al. [48]. It was also observed that increased somatic false positives are generated by the in silico germline filtering method for patients with non-caucasian ancestry compared to caucasian patients, as the former are less represented in public databases used for germline variant filtering [34]. The use of ExAC, the largest and more representative public germline WES database, in addition to dbSNP and 1000 Genomes, is recommended to reduce this difference [34]. 

![Fig. 3](http://jitc.bmj.com/https://jitc.bmj.com/content/jitc/7/1/183/F3.medium.gif)

[Fig. 3](http://jitc.bmj.com/content/7/1/183/F3)

Fig. 3 
Differences in the workflow for panel-based TMB quantification. **a**. Overview of the factors influencing panel-based TMB quantification. Several variables in library construction, sequencing and in the pipeline to call mutations influence panel-based TMB quantification. Furthermore, panel-based TMB quantification is influenced by differences in the bioinformatic method to extrapolate global TMB from mutations identified in the narrow genomic region targeted by the gene panel. **b**. Differences across various studies in panel-based TMB quantification: gene panel technical specifications, preanalytical factors and the bioinformatics workflow used to extrapolate from the genomic space targeted by gene panels global TMB are described. *FM1*: Foundation Medicine’s FoundationOne panel (v1: 185 genes, v2: 236 genes, v3: 315 genes, v4: 405 genes); *NA:* not available; *±*: algorithm developed by Sun et al. for in silico removal of germline variants [74]

In addition to these factors, which similarly influence WES- and panel-based analyses, panel-based TMB quantification also requires to extrapolate the number of somatic coding mutations observed in the targeted genomic space to the number that would be observed across the whole genome. Extrapolation methods may differ for various choices in variant filtering, such as removal of known cancer mutations or synonymous mutations (Fig. 3b). Standard gene panels are commonly enriched in known cancer genes, which are more likely to be mutated in a tumor and expectedly enriched in mutations. Therefore, it was proposed to remove known cancer variants of targeted genes when performing TMB quantification, to avoid overestimation of TMB when extrapolating it across the whole genome [12, 22]. Buchhalter et al., showed that removal of cancer mutational hotspots slightly decreases the number of high TMB tumors identified but does not change the general picture [35]. However, the importance of this filtering, routinely performed only for Foundation Medicine panels, depends on panel size and composition: some gene panels may be larger and less enriched in cancer genes by including, for example, pharmacogenomic variants. As for synonymous mutations, it is claimed that, although not biologically meaningful, their inclusion may reduce sampling noise and improve the approximation of TMB across the whole genome. Indeed, several works compared TMB quantification with or without synonymous variants and observed that, when including synonymous variants, panel-based TMB shows increased correlation with WES-based TMB values [35, 49] and stronger association with clinical response [9]. Starting from the above observations, we can infer that some recommendations to build a standardized and robust analysis pipeline for TMB quantification are starting to emerge at least for the following points: i. germline variants can be most efficiently identified and removed by matched normal sequencing and, if this is not possible, the largest available germline variant databases should be used for in silico filtering, especially for non-caucasian ancestries; ii. TMB extrapolation to the whole genome is accurately performed by counting all somatic mutations, including synonymous mutations, to enlarge the sampling space and better approximate global TMB across the whole genome.

As already mentioned, panel-based TMB quantification is affected by the genomic size targeted by the panel and by its gene composition. Notably, gene panels tested up to now widely differ for number of targeted genes (from 73 to 710) and size (from 0.39 to 2.8 Mb of targeted genomic space) (Additional file 1: Table S1). These considerations raise the question of how to convert TMB estimates between different gene panels to allow cross-platform comparability. Indeed, although the majority of studies correlating TMB to immunotherapy response are currently based on FoundationOne, several other types of gene panels exist and the offer is steadily increasing (Additional file 1: Table S1). Moreover, we still need to standardize the conversion of the reference WES-based TMB values to panel-based TMB, since the lower sequencing coverage and higher sequencing depth of gene panels, as compared to WES, may lead to decreased accuracy of TMB values and increased sensitivity in variant calling. For cross-panels or panel-to-WES TMB conversion, an in silico approach was proposed, where TMB distributions derived from two different technologies were interpolated and aligned and TMB cutoffs were mapped across distributions [38]. However, a consensus on a standard method to convert TMB values is still missing.

Not surprisingly, in this heterogeneous landscape a robust TMB cutoff to discriminate between immunotherapy responders and non-responders is still to be defined. Moreover, the adopted cutoffs sometimes differ across different studies on the same gene panel (Table 1). Up to now, the TMB cutoff of 10 mutations per Mb, measured by the FoundationOne gene panel and found to best discriminate between responders and non-responders to immunotherapy in NSCLC patients, is the only one which has been validated in a separate further study [28, 50, 51]; this cutoff was also observed, but not yet validated, in melanoma [38] and in metastatic urothelial carcinoma [15] (Table 1). Interestingly, these cancer types present a TMB distribution similar to that of NSCLC [52]. Indeed, due to the diversity of TMB distribution across different cancer types, the adoption of cancer-specific TMB cutoffs was proposed [35, 43]. TMB cutoff was initially most commonly established using distribution-based stratification, which can be heavily influenced by outliers, while it is now often identified based on statistically sound methods, such as ROC curves. In the future, the use of ROC curves as a common method of TMB cutoff determination will greatly help to converge to a robust TMB cutoff and will facilitate comparison across different platforms.

The heterogeneity in experimental and analytical protocols, in the extrapolation of panel-based TMB values and in gene panel technical specifications is currently limiting the potential use of TMB in a clinical setting (Fig. 3). For this reason, a common standard for TMB quantification and a consensus on a clinically useful TMB cutoff are urgently needed. Some efforts in this sense are ongoing by the Quality Assurance Initiative Pathology (QuIP) in Germany ([https://quip.eu](https://quip.eu)) and by the European Society of Pathology. Moreover, in the US, governmental organizations, health-sector industries, several NGS gene panel manufacturers and academic institutions set up a TMB Harmonization Working Group ([https://www.focr.org/tmb](https://www.focr.org/tmb)) and planned a 3-phase project for TMB harmonization. In the first phase, they performed in silico analyses of publicly available TCGA data to identify sources of variability in TMB quantification between WES and gene panels. Recently concluded, this work established that panel-based TMB is comparable between different gene panels for TMB values ranging 0 to 40 mutations per Mb, that it strongly correlates with WES-based TMB and it is possibly influenced by the type of cancer under investigation. It also found that the observed variance across gene panels stems from their different gene composition and technical specifications, as well as from the bioinformatic pipeline adopted [53]. The second step of the project empirically validates TMB estimates from different gene panels by mapping them to WES-based TMB values, used as a gold standard, whereas the last step will define best practices for TMB use as immunotherapy biomarker in clinical routine. Following preliminary results on the influence of the bioinformatic workflow and of gene panel size and composition on TMB quantification, the working group recommends the use of gene panels larger than 1 Mb and the standardization of the bioinformatic algorithms, in addition to standardization of sample processing. Moreover, it suggests the inclusion of actionable genes, genes associated with mutagenesis and negative predictors of response in these gene panels and the alignment of panel-based TMB values to WES-based ones to allow interoperability across different assays [54].

## TMB quantification beyond tissue biopsies and current gene panels

Most studies on TMB as a predictive biomarker for immunotherapy response were performed on bioptical or surgical specimens from solid tumors. Since obtaining tissue biopsies may be challenging and invasive for patients, it would be critical for the clinical routine to assess TMB using cell-free DNA (cfDNA) from blood, which includes circulating tumor DNA (ctDNA), as a surrogate specimen to biopsy. High throughput molecular profiling of ctDNA remains technically challenging but increasing efforts are being made in this direction. A few studies previously investigated the feasibility of WES on ctDNA and highlighted some inherent limitations, such as the low amount of available ctDNA, which reduces sensitivity, or ctDNA being more associated with metastases rather than with primary tumors [55–59]. In one of the largest studies attempting optimization of WES-based TMB quantification from liquid biopsy, WES was performed in parallel on DNA from tissue biopsies and on cfDNA from liquid biopsies of 32 metastatic patients and comparable sequencing depth and coverage were obtained [60]. Performance of variant detection was dependent on the fraction of tumor DNA within the analyzed cfDNA, as previously described. In those samples positive for the presence of ctDNA, variant detection sensitivity of cfDNA-WES compared to tDNA-WES was 94%, regardless of the tumor type (2 cholangiosarcoma and 19 lung, 5 head and neck, 2 prostate, 2 colorectal, 1 breast and 1 bladder cancer were analyzed). Most importantly, in ctDNA positive samples, TMB values from WES on liquid biopsies were robust and consistent with those from WES on tissue biopsies, which demonstrates for the first time the feasibility of TMB quantification from liquid biopsies, using WES.

More recently, it was demonstrated that targeted enrichment sequencing by gene panels is another valid approach for TMB quantification from liquid biopsies. In particular, Gandara et al. developed, tested and analytically validated a novel gene panel for TMB quantification from blood [61]. The panel is based on hybridization-capture enrichment and targets 394 genes, corresponding to 1.1 Mb of genomic space (Fig. 3). Its clinical utility was evaluated via a retrospective study on 259 NSCLC samples from patients treated with immunotherapy or chemotherapy in the OAK and POPLAR clinical trials. Blood-derived TMB (bTMB) calculated using this novel gene panel correlated well with tissue-derived TMB (tTMB) measured by FoundationOne. Moreover, measured TMB was found to be significantly associated with response to anti-PD-L1 immunotherapy in the POPLAR trial and this was further confirmed on patient samples from the OAK trial. A prospective validation is also currently ongoing in the BFAST trial ([NCT03178552](http://jitc.bmj.com/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT03178552&atom=%2Fjitc%2F7%2F1%2F183.atom)) on advanced and metastatic NSCLC patients. Interestingly, it was observed that the capability of TMB, as measured by this panel, to predict anti-PD-1/PD-L1 immunotherapy response is independent from PD-L1 expression levels [61]. One of the main pitfalls of the panel, though, is its limit of detection, defined as a minimum of 1% tumor content in at least 20 ng of cell-free DNA input, and its dependency on the overall tumor burden, which influences the likelihood of detecting ctDNA. The exclusive use of single nucleotide variants (SNVs) for TMB quantification represents another limitation, although future versions of the algorithm are planned to be released, which will also use indels. The commercial Guardant360 and GuardantOMNI gene panels were also designed for blood-based TMB quantification [62]. Their limit of detection was defined as a minimum of 0.3% tumor content in at least 5 ng of cell-free DNA input. They were validated in silico by subsetting TCGA WES datasets to only include genes targeted by the panels. Panel accuracy in TMB quantification was then evaluated by correlation of TMB values obtained from the simulated gene panels with those from WES. Their predictive value was similarly evaluated in silico on 30 lung cancer samples with matched information on immunotherapy response. The performance showed by TMB from the simulated gene panel in responder identification was comparable to that of WES-based TMB (Additional file 4: Table S4). The Guardant360 panel was further tested in a small retrospective study on 69 patients with various tumor types [44]. No comparison with tissue-based TMB has been reported yet, but a significant correlation between high blood-derived TMB measured by Guardant360 and immunotherapy response was observed [44]. Finally, a further gene panel for bTMB quantification was recently developed in China. Consistency between panel-based bTMB values and WES-based tTMB values, tested in silico and empirically by matched blood and tissue samples, was comparable to that of the panels described above. Similar results were also found for its predictive value, based on in silico analyses. Interestingly, the authors also raised the issue of the different frequency of oncogenic driver mutations, such as EGFR or KRAS, between Asian and white population. For this reason, they compare TCGA WES-based TMB with panel-based TMB with or without inclusion of EGFR and/or KRAS mutations. Although similar results are yielded, the raised issue is an important point to be further investigated in panel-based TMB quantification [49]. TMB quantification from liquid biopsies suffers from ctDNA detection limits, which also depend on tumor size and number of cancer cells, but these results encourage to further explore and more extensively validate this approach.

Besides new technologies to estimate TMB from liquid biopsies, another significant step towards routine use of TMB in clinical practice is TMB quantification from an even smaller set of genes than in targeted enrichment gene panels. Although panel size is known to affect accuracy of TMB quantification, the use of a highly customized set of genes may represent a valid and even less expensive approach. In this view, Lyu et al., proposed a computational framework to assess the best and smallest subset of genes necessary to estimate TMB as a biomarker for lung adenocarcinoma [63]. They were able to identify a model of only 24 genes which predicted in silico immunotherapy response with 93% specificity and 85% sensitivity and they suggested that other small custom sequencing gene panels may be designed in a cancer-specific way to assess TMB with further reduced costs.

## Future perspectives and recommendations

TMB is one of the most rapidly developing biomarkers for immunotherapy response, with about 37 ongoing clinical trials currently registered in [ClinicalTrials.gov](http://clinicaltrials.gov) that use TMB as stratification biomarker [64]. Several gene panels were recently optimized to estimate TMB at reduced sequencing costs, and emerging evidence supports the feasibility of TMB quantification from liquid biopsies. However, harmonization in TMB quantification and reporting remains the main challenge for the near future: standard procedures are required to allow interoperability between different gene panels, compare results across studies and define a universal cutoff to confidently identify patients most likely to benefit from immunotherapy.

Even an accurate TMB value is an imperfect predictor of immunotherapy response and further studies are needed to enhance its value as clinically useful immunotherapy biomarker. TMB is used as an approximation of neoantigen burden, upon the assumption that the higher the mutational burden, the higher the probability for immunogenic peptides to be generated, which leads to stronger immune response upon inhibition of immune checkpoints. Interestingly, neoantigen clonality, in addition to the overall amount of neoantigens, influences immunotherapy response in NSCLC patients [65]. In particular, tumors enriched in clonal neoantigens (e.g. present in all tumor cells) are more sensitive to immune checkpoint inhibitors than tumors enriched in subclonal neoantigens (e.g. present only in a subset of tumor cells), in advanced NSCLC and melanoma patients [65]. Indeed, clonality of produced neoantigens seems to be associated with a more effective immune surveillance. On the other hand, enrichment in subclonal neoantigens may activate T cells against only a subset of tumor cells, leading to less effective tumor control. Based on these observations, it would be interesting to investigate if information on mutation clonality (e.g. variant allele frequency) improves the predictive power of TMB. Evaluation of mutation clonality from gene panels is not trivial though: the reduced genomic space targeted by gene panels may not be representative of the overall clonal architecture and the mutations sampled herein may not be those generating neoantigens. Interestingly, McGranahan et al. observed a relationship between subclonal mutations and mutational signatures associated with alkylating agents and, in NSCLC, between clonal mutations and mutational signatures associated with smoking [65]. Mutational signatures associated with smoking were also found to be significantly associated with high tumor mutational burden and with response to immunotherapy [9]. Therefore, although the extraction of mutational signatures from gene panels may be hampered by the small number of sampled mutations, these observations suggest that they may prove helpful to infer neoantigen clonality and enhance TMB predictive value.

Integration of TMB with other potential immunotherapy biomarkers represents another promising way to refine prediction of immunotherapy responders. For example, TMB, defects in DNA mismatch-repair pathway and the MSI status all are measures of genomic instability that can provide indirect assessment of tumor antigenicity, whereas PD-L1 expression, immune cell infiltration and inflammatory signatures represent biomarkers of the T cell-inflamed tumor microenvironment. Therefore, their integration can refine prediction of immunotherapy outcome by combining information on tumor complexity and on the immune response. Indeed, emerging evidence suggests that, at least in NSCLC, TMB and PD-L1 expression are independent predictors and TMB may complement or even outperform PD-L1 expression [10, 26, 50, 66]. Moreover, it was observed that most tumors with high MSI also present elevated TMB, whereas the opposite does not hold true. The combination of TMB with MSI and PD-L1 expression in gastrointestinal tumors significantly improved the identification of immunotherapy responders [67]. In another study, it was observed that TMB is an independent predictor and only weakly correlates with T cell-inflamed gene expression profiles (GEP) or PD-L1 expression. Thus, TMB and T cell-inflamed GEP were jointly used to identify immunotherapy responders: patients with both high TMB and high T cell-inflamed GEP were those with the highest objective response rates on tumors from four KEYNOTE clinical trials across 22 cancer types. Similarly, in melanoma patients, a response score based on the combination of TMB, infiltration of CD8+ T cells and gene expression profiles for PD-L1, CD8 and a set of 394 immune genes demonstrated higher sensitivity and similar specificity than each biomarker alone [68]. To date, the FoundationOne and Guardant360 gene panels allow to measure both TMB and MSI but no other potential immunotherapy biomarker. Moreover, they do not provide the user any combinatorial model to integrate them. Although further validation in prospective clinical studies is required for all these potential biomarkers, several observations suggest that simultaneous profiling of both TMB and other immunotherapy biomarkers currently under investigation may represent the next step forward in the design of new gene panels for clinical use. The Friends and QuIP initiatives for TMB harmonization recommended to include as much relevant genetic and molecular information as possible in these panels, to avoid the need to re-biopsy the patient for further information. In line with this recommendation, we propose to also include in gene panels for TMB quantification other potential immunotherapy biomarkers but also negative predictors of immunotherapy response [69, 70] and variants predisposing to adverse reaction to immunotherapy [71, 72]. These and other recommendations which emerge from the studies reviewed here, including the one from the TMB Harmonization Working Group, are summed up in Additional file 6: Table S6.

## Funding

LF work was supported by the Italian Ministry of Health (Ricerca Corrente, “Alleanza Contro il Cancro” - ACC network).

## Acknowledgments

We thank all researchers, clinicians and organizations working in this field for their contributions and we apologize to those whose work we did not report or cite.

## Additional files

[[40425\_2019\_647\_MOESM1\_ESM.xlsx]](pending:yes)

[[40425\_2019\_647\_MOESM2\_ESM.xlsx]](pending:yes)

[[40425\_2019\_647\_MOESM3\_ESM.xlsx]](pending:yes)

[[40425\_2019\_647\_MOESM4\_ESM.xlsx]](pending:yes)

[[40425\_2019\_647\_MOESM5\_ESM.xlsx]](pending:yes)

[[40425\_2019\_647\_MOESM6\_ESM.xlsx]](pending:yes)

[[40425\_2019\_647\_MOESM7\_ESM.pdf]](pending:yes)

[[40425\_2019\_647\_MOESM8\_ESM.pdf]](pending:yes)

[[40425\_2019\_647\_MOESM9\_ESM.pdf]](pending:yes)

## Authors’ contributions

LF designed, wrote and revised the manuscript. PGP, SG and LM revised the manuscript. All authors read and approved the final manuscript.

## Availability of data and materials

Not applicable.

## Notes

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

LF, SG, PGP, and LM declare that they have no competing interest.

## Footnotes

*   Electronic supplementary material The online version of this article (10.1186/s40425-019-0647-4) contains supplementary material, which is available to authorized users.

ACC
:   Adrenocortical carcinoma
AUC
:   Area under the curve
BLCA
:   Bladder urothelial carcinoma
BRCA
:   Breast invasive carcinoma
CESC
:   Cervical squamous cell carcinoma and endocervical adenocarcinoma
cfDNA
:   Circulating free DNA
CHOL
:   Cholangiosarcoma
COADREAD
:   Colon adenocarcinoma
CRC
:   Colorectal cancer
ctDNA
:   Circulating tumor DNA
DLBC
:   Lymphoid neoplasm diffuse large B-cell lymphoma
ESCA
:   Esophageal carcinoma
FDA
:   Food and Drug Administration
GBM
:   Glioblastoma
GEP
:   Gene expression profile
HLA
:   Human Leukocyte Antigen
HNSC
:   Head and neck squamous cell carcinoma
KICH
:   Kidney chromophobe
KIRC
:   Kidney renal clear cell carcinoma
KIRP
:   Kidney renal papillary cell carcinoma
LAML
:   Acute myeloid leukemia
LGG
:   Brain lower grade glioma
LIHC
:   Liver hepatocellular carcinoma
LUAD
:   Lung adenocarcinoma
LUSC
:   Lung squamous carcinoma
Mb
:   Megabase
mCRPC
:   Metastatic castration-resistant prostate cancer
MESO
:   Mesothelioma
MHC
:   Major histocompatibility complex
MMR
:   Mismatch repair
MSI
:   Microsatellite instability
NSCLC
:   Non-small cell lung cancer
ORR
:   Objective response rates
OS
:   Overall survival
OV
:   Ovarian serous cystadenocarcinoma
PAAD
:   Pancreatic adenocarcinoma
PCPG
:   Pheochromocytoma and paraganglioma
PFS
:   Progression free survival
PRAD
:   Prostate adenocarcinoma
ROC
:   Receiver operating characteristic
SARC
:   Sarcoma
SCLC
:   Small cell lung cancer
SKCM
:   Skin cutaneous melanoma
SNV
:   Single nucleotide variant
STAD
:   Stomach adenocarcinoma
TCGA
:   The cancer genome atlas
TCR
:   T cell receptor
TGCT
:   Testicular germ cell tumors
THYM
:   Thymoma
TMB
:   Tumor mutational burden
UCEC
:   Uterine corpus endometrial carcinoma
UCS
:   Uterine carcinosarcoma
UVM
:   Uveal melanoma
WES
:   Whole exome sequencing

*   Received December 20, 2018.
*   Accepted June 19, 2019.


*   © The Author(s).

**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/](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/](http://creativecommons.org/publicdomain/zero/1.0/)) applies to the data made available in this article, unless otherwise stated.

## References

1.  1. Antonia SJ , Villegas A , Daniel D , Vicente D , Murakami S , Ri H , et al . Durvalumab after chemoradiotherapy in stage III non-small cell lung cancer. N Engl J Med. 2017;377:1919–1929. [doi:10.1056/NEJMoa1709937](http://dx.doi.org/10.1056/NEJMoa1709937) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1056/NEJMoa1709937&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=28885881&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

2.  2. Borghaei H , Paz-Ares L , Horn L . Nivolumab versus docetaxel in advanced nonsquamous non-small-cell lung cancer. N Engl J Med. 2015;373:1627–1639. pmcid:5705936 [doi:10.1056/NEJMoa1507643](http://dx.doi.org/10.1056/NEJMoa1507643) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1056/NEJMoa1507643&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=26412456&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

3.  3. Hodi F , O’Day S , McDermott D , Weber R , Sosman J , Haanen J , et al . Improved survival with ipilimumab in patients with metastatic melanoma. N Engl J Med. 2010;363:711–723. pmcid:20525992 [doi:10.1056/NEJMoa1003466](http://dx.doi.org/10.1056/NEJMoa1003466) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1056/NEJMoa1003466&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=20525992&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

4.  4. Motzer R , Tannir N , McDermott D , Frontera O , Melichar B , Choueiri T , et al . Nivolumab plus Ipilimumab versus Sunitinib in advanced renal-cell carcinoma. N Engl J Med. 2018;378:14 1277–1290. pmcid:5972549 [doi:10.1056/NEJMoa1712126](http://dx.doi.org/10.1056/NEJMoa1712126) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1056/NEJMoa1712126&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=29562145&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

5.  5. Rosenberg Jonathan E , Hoffman-Censits Jean , Powles Tom , van der Heijden Michiel S , Balar Arjun V , Necchi Andrea , Dawson Nancy , O'Donnell Peter H , Balmanoukian Ani , Loriot Yohann , Srinivas Sandy , Retz Margitta M , Grivas Petros , Joseph Richard W , Galsky Matthew D , Fleming Mark T , Petrylak Daniel P , Perez-Gracia Jose Luis , Burris Howard A , Castellano Daniel , Canil Christina , Bellmunt Joaquim , Bajorin Dean , Nickles Dorothee , Bourgon Richard , Frampton Garrett M , Cui Na , Mariathasan Sanjeev , Abidoye Oyewale , Fine Gregg D , Dreicer Robert . Atezolizumab in patients with locally advanced and metastatic urothelial carcinoma who have progressed following treatment with platinum-based chemotherapy: a single-arm, multicentre, phase 2 trial. The Lancet. 2016;387:10031 1909–1920. [doi:10.1016/S0140-6736(16)00561-4](http://dx.doi.org/10.1016/S0140-6736(16)00561-4) 
    
    

6.  6. Li Xiaolei , Song Wenhui , Shao Changshun , Shi Yufang , Han Weidong . Emerging predictors of the response to the blockade of immune checkpoints in cancer therapy. Cellular & Molecular Immunology. 2018;16:1 28–39. [doi:10.1038/s41423-018-0086-z](http://dx.doi.org/10.1038/s41423-018-0086-z) 
    
    

7.  7. Galluzzi L , Chan TA , Kroemer G , Wolchok JD , López-Soto A . The hallmarks of successful anticancer immunotherapy. Sci Transl Med. 2018;10:459 1–15. [doi:10.1126/scitranslmed.aat7807](http://dx.doi.org/10.1126/scitranslmed.aat7807) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

8.  8. Campesato LF , Barroso-sousa R , Jimenez L , Camargo AA . Comprehensive cancer-gene panels can be used to estimate mutational load and predict clinical benefit to PD-1 blockade in clinical practice. Oncotarget.. 2015;6:33 34221. pmcid:4741447 [doi:10.18632/oncotarget.5950](http://dx.doi.org/10.18632/oncotarget.5950) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.18632/oncotarget.5950&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

9.  9. Rizvi NA , Hellmann MD , Snyder A , Kvistborg P , Makarov V , Havel JJ , et al . Mutational landscape determines sensitivity to PD-1 blockade in non – small cell lung cancer. Science.. 2015;348:6230 124–129. pmcid:4993154 [doi:10.1126/science.aaa1348](http://dx.doi.org/10.1126/science.aaa1348) 
    
    [Abstract/FREE Full Text](http://jitc.bmj.com/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNDgvNjIzMC8xMjQiO3M6NDoiYXRvbSI7czoxODoiL2ppdGMvNy8xLzE4My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

10. 10. Carbone David P. , Reck Martin , Paz-Ares Luis , Creelan Benjamin , Horn Leora , Steins Martin , Felip Enriqueta , van den Heuvel Michel M. , Ciuleanu Tudor-Eliade , Badin Firas , Ready Neal , Hiltermann T. Jeroen N. , Nair Suresh , Juergens Rosalyn , Peters Solange , Minenza Elisa , Wrangle John M. , Rodriguez-Abreu Delvys , Borghaei Hossein , Blumenschein George R. , Villaruz Liza C. , Havel Libor , Krejci Jana , Corral Jaime Jesus , Chang Han , Geese William J. , Bhagavatheeswaran Prabhu , Chen Allen C. , Socinski Mark A. . First-Line Nivolumab in Stage IV or Recurrent Non–Small-Cell Lung Cancer. New England Journal of Medicine. 2017;376:25 2415–2426. [doi:10.1056/NEJMoa1613493](http://dx.doi.org/10.1056/NEJMoa1613493) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1056/NEJMoa1613493&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=28636851&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

11. 11. 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–211. pmcid:5054517 [doi:10.1126/science.aad0095](http://dx.doi.org/10.1126/science.aad0095) 
    
    [Abstract/FREE Full Text](http://jitc.bmj.com/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNTAvNjI1Ny8yMDciO3M6NDoiYXRvbSI7czoxODoiL2ppdGMvNy8xLzE4My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

12. 12. Johnson D. B. , Frampton G. M. , Rioth M. J. , Yusko E. , Xu Y. , Guo X. , Ennis R. C. , Fabrizio D. , Chalmers Z. R. , Greenbowe J. , Ali S. M. , Balasubramanian S. , Sun J. X. , He Y. , Frederick D. T. , Puzanov I. , Balko J. M. , Cates J. M. , Ross J. S. , Sanders C. , Robins H. , Shyr Y. , Miller V. A. , Stephens P. J. , Sullivan R. J. , Sosman J. A. , Lovly C. M. . Targeted Next Generation Sequencing Identifies Markers of Response to PD-1 Blockade. Cancer Immunology Research. 2016;4:11 959–967. pmcid:5134329 [doi:10.1158/2326-6066.CIR-16-0143](http://dx.doi.org/10.1158/2326-6066.CIR-16-0143) 
    
    [Abstract/FREE Full Text](http://jitc.bmj.com/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FuaW1tIjtzOjU6InJlc2lkIjtzOjg6IjQvMTEvOTU5IjtzOjQ6ImF0b20iO3M6MTg6Ii9qaXRjLzcvMS8xODMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

13. 13. Eroglu Z , Zaretsky JM , Hu-Lieskovan S , Kim DW , Algazi A , Johnson DB , et al . High response rate to PD-1 blockade in desmoplastic melanomas. Nature.. 2018;553:7688 347–350. pmcid:5773412 [doi:10.1038/nature25187](http://dx.doi.org/10.1038/nature25187) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1038/nature25187&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=29320474&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

14. 14.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. 2014;2189–2199.
    
    

15. 15. Powles Thomas , Durán Ignacio , van der Heijden Michiel S , Loriot Yohann , Vogelzang Nicholas J , De Giorgi Ugo , Oudard Stéphane , Retz Margitta M , Castellano Daniel , Bamias Aristotelis , Fléchon Aude , Gravis Gwenaëlle , Hussain Syed , Takano Toshimi , Leng Ning , Kadel Edward E , Banchereau Romain , Hegde Priti S , Mariathasan Sanjeev , Cui Na , Shen Xiaodong , Derleth Christina L , Green Marjorie C , Ravaud Alain . Atezolizumab versus chemotherapy in patients with platinum-treated locally advanced or metastatic urothelial carcinoma (IMvigor211): a multicentre, open-label, phase 3 randomised controlled trial. The Lancet. 2018;391:10122 748–757. [doi:10.1016/S0140-6736(17)33297-X](http://dx.doi.org/10.1016/S0140-6736(17)33297-X) 
    
    

16. 16. Balar AV , Galsky MD , Rosenberg JE , Powles T , Petrylak DP , Bellmunt J , et al . Atezolizumab as first-line treatment in cisplatin-ineligible patients with locally advanced and metastatic urothelial carcinoma: a single-arm, multicentre, phase 2 trial. Lancet.. 2017;389:10064 67–76. [doi:10.1016/S0140-6736(16)32455-2](http://dx.doi.org/10.1016/S0140-6736(16)32455-2) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1016/S0140-6736(16)32455-2&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=27939400&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

17. 17.Le DT, Uram JN, H W, R BB, Kemberling H, Eyring AD, et al. PD-1 blockade in tumors with mismatch-repair deficiency. NEJM. 2015:2509–20.
    
    

18. 18. Goodman Aaron M. , Kato Shumei , Bazhenova Lyudmila , Patel Sandip P. , Frampton Garrett M. , Miller Vincent , Stephens Philip J. , Daniels Gregory A. , Kurzrock Razelle . Tumor Mutational Burden as an Independent Predictor of Response to Immunotherapy in Diverse Cancers. Molecular Cancer Therapeutics. 2017;16:11 2598–2608. pmcid:5670009 [doi:10.1158/1535-7163.MCT-17-0386](http://dx.doi.org/10.1158/1535-7163.MCT-17-0386) 
    
    [Abstract/FREE Full Text](http://jitc.bmj.com/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6Im1vbGNhbnRoZXIiO3M6NToicmVzaWQiO3M6MTA6IjE2LzExLzI1OTgiO3M6NDoiYXRvbSI7czoxODoiL2ppdGMvNy8xLzE4My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

19. 19. Lawrence MS , Stojanov P , Polak P , Kryukov GV , Cibulskis K , Sivachenko A , et al . Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature.. 2013;499:7457 214–218. pmcid:3919509 [doi:10.1038/nature12213](http://dx.doi.org/10.1038/nature12213) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1038/nature12213&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=23770567&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 
    
    [Web of Science](http://jitc.bmj.com/lookup/external-ref?access_num=000321557600063&link_type=ISI) 

20. 20. Alexandrov LB , Nik-Zainal S , Wedge DC , Aparicio SAJR , Behjati S , Biankin AV , et al . Signatures of mutational processes in human cancer. Nature.. 2013;500:7463 415–421. pmcid:3776390 [doi:10.1038/nature12477](http://dx.doi.org/10.1038/nature12477) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1038/nature12477&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=23945592&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 
    
    [Web of Science](http://jitc.bmj.com/lookup/external-ref?access_num=000323316100026&link_type=ISI) 

21. 21.Quiroga D, Lyerly HK, Morse MA. Deficient mismatch repair and the role of immunotherapy in metastatic colorectal cancer. Curr Treat Options in Oncol. 2016. 10.1007/s11864-016-0414-4.
    
    

22. 22. 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 Med.. 2017;9:1 1–14. [doi:10.1186/s13073-017-0424-2](http://dx.doi.org/10.1186/s13073-017-0424-2) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1186/s13073-017-0403-7&link_type=DOI) 

23. 23. Rooney Michael S. , Shukla Sachet A. , Wu Catherine J. , Getz Gad , Hacohen Nir . Molecular and Genetic Properties of Tumors Associated with Local Immune Cytolytic Activity. Cell. 2015;160:1-2 48–61. pmcid:4856474 [doi:10.1016/j.cell.2014.12.033](http://dx.doi.org/10.1016/j.cell.2014.12.033) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1016/j.cell.2014.12.033&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=25594174&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 
    
    [Web of Science](http://jitc.bmj.com/lookup/external-ref?access_num=000347923200008&link_type=ISI) 

24. 24.Zhang J, Mardis ER, Maher CA. Genome analysis INTEGRATE-neo : a pipeline for personalized gene fusion neoantigen discovery. Bioinformatics. 2017; doi:10.1093/bioinformatics/btw674.
    
    

25. 25. Yadav M , Jhunjhunwala S , Phung QT , Lupardus P , Tanguay J , Bumbaca S , et al . Predicting immunogenic tumour mutations by combining mass spectrometry and exome sequencing. Nature.. 2014;515:7528 572–576. [doi:10.1038/nature14001](http://dx.doi.org/10.1038/nature14001) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1038/nature14001&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=25428506&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 
    
    [Web of Science](http://jitc.bmj.com/lookup/external-ref?access_num=000346247600055&link_type=ISI) 

26. 26. Hellmann M , Nathanson T , Rizvi H . Genomic features of response to combination immunotherapy in patients with advanced non- small-cell lung cancer. Cancer Cell. 2018;33:843–852. pmcid:5953836 [doi:10.1016/j.ccell.2018.03.018](http://dx.doi.org/10.1016/j.ccell.2018.03.018) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1016/j.ccell.2018.03.018&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

27. 27. Kowanetz M . Tumor mutation load assessed by FoundationOne (FM1) is associated with improved efficacy of atezolizumab (atezo) in patients with advanced NSCLC. Ann Oncol. 2016;27:6 15–42. 
    
    

28. 28. Ready Neal , Hellmann Matthew D. , Awad Mark M. , Otterson Gregory A. , Gutierrez Martin , Gainor Justin F. , Borghaei Hossein , Jolivet Jacques , Horn Leora , Mates Mihaela , Brahmer Julie , Rabinowitz Ian , Reddy Pavan S. , Chesney Jason , Orcutt James , Spigel David R. , Reck Martin , O’Byrne Kenneth John , Paz-Ares Luis , Hu Wenhua , Zerba Kim , Li Xuemei , Lestini Brian , Geese William J. , Szustakowski Joseph D. , Green George , Chang Han , Ramalingam Suresh S. . First-Line Nivolumab Plus Ipilimumab in Advanced Non–Small-Cell Lung Cancer (CheckMate 568): Outcomes by Programmed Death Ligand 1 and Tumor Mutational Burden as Biomarkers. Journal of Clinical Oncology. 2019;37:12 992–1000. [doi:10.1200/JCO.18.01042](http://dx.doi.org/10.1200/JCO.18.01042) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

29. 29.Galsky M, Saci A, Szabo P, Azrilevich A, Horak C, Lambert A, et al. Impact of tumor mutation burden on Nivolumab efficacy in second-line urothelial carcinoma patients: exploratory analysis of the phase II CheckMate 275. Ann Oncol. 2017. 10.1093/annonc/mdx371.
    
    

30. 30. Strickland KC , Howitt BE , Shukla SA , Rodig S , Ritterhouse LL , Liu JF , et al . Association and prognostic significance of BRCA1/2-mutation status with neoantigen load, number of tumor-infiltrating lymphocytes and expression of PD-1/PD-L1 in high grade serous ovarian cancer. Oncotarget.. 2016;7:12 13587–13598. pmcid:4924663 [doi:10.18632/oncotarget.7277](http://dx.doi.org/10.18632/oncotarget.7277) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=26871470&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

31. 31. Birkbak Nicolai Juul , Kochupurakkal Bose , Izarzugaza Jose M. G. , Eklund Aron C. , Li Yang , Liu Joyce , Szallasi Zoltan , Matulonis Ursula A. , Richardson Andrea L. , Iglehart J. Dirk , Wang Zhigang C. . Tumor Mutation Burden Forecasts Outcome in Ovarian Cancer with BRCA1 or BRCA2 Mutations. PLoS ONE. 2013;8:11 e80023. pmcid:3827141 [doi:10.1371/journal.pone.0080023](http://dx.doi.org/10.1371/journal.pone.0080023) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1371/journal.pone.0080023&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=24265793&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

32. 32. Thomas A , Routh ED , Pullikuth A , Jin G , Su J , Chou JW , et al . Tumor mutational burden is a determinant of immune-mediated survival in breast cancer. Oncoimmunology.. 2018;7:10 1–12. [doi:10.1080/2162402X.2018.1490854](http://dx.doi.org/10.1080/2162402X.2018.1490854) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1080/2162402X.2018.1501137&link_type=DOI) 

33. 33. Jayaraman Shyam S. , Rayhan David J. , Hazany Salar , Kolodney Michael S. . Mutational Landscape of Basal Cell Carcinomas by Whole-Exome Sequencing. Journal of Investigative Dermatology. 2014;134:1 213–220. [doi:10.1038/jid.2013.276](http://dx.doi.org/10.1038/jid.2013.276) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1038/jid.2013.276&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=23774526&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

34. 34.Garofalo A, Sholl L, Reardon B, Taylor-Weiner A, Amin-Mansour A, Miao D, et al. The impact of tumor profiling approaches and genomic data strategies for cancer precision medicine. Genome Med. 2016. 10.1186/s13073-016-0333-9.
    
    

35. 35. Buchhalter Ivo , Rempel Eugen , Endris Volker , Allgäuer Michael , Neumann Olaf , Volckmar Anna-Lena , Kirchner Martina , Leichsenring Jonas , Lier Amelie , von Winterfeld Moritz , Penzel Roland , Christopoulos Petros , Thomas Michael , Fröhling Stefan , Schirmacher Peter , Budczies Jan , Stenzinger Albrecht . Size matters: Dissecting key parameters for panel-based tumor mutational burden analysis. International Journal of Cancer. 2018;144:4 848–858. [doi:10.1002/ijc.31878](http://dx.doi.org/10.1002/ijc.31878) 
    
    

36. 36.Qiu P, Poehlein CH, Marton MJ, Laterza OF, Levitan D. Measuring tumor mutational burden (TMB) in plasma from mCRPC patients using two commercial NGS assays. Sci Rep. 2019. 10.1038/s41598-018-37128-y.
    
    

37. 37. Nguyen Andrew , Garner Chad , Reddy Sandeep K. , Sanborn John Z , Benz Stephen Charles , Seery Tara Elisabeth , Rabizadeh Shahrooz . Three-fold overestimation of tumor mutation burden using 248 gene panel versus whole exome. Journal of Clinical Oncology. 2018;36:15_suppl 12117–12117. [doi:10.1200/JCO.2018.36.15\_suppl.12117](http://dx.doi.org/10.1200/JCO.2018.36.15_suppl.12117) 
    
    

38. 38.Panda A, Betigeri A, Subramanian K, Ross JS, Pavlick DC, Ali S, et al. Identifying a clinically applicable mutational burden threshold as a potential biomarker of response to immune checkpoint therapy in solid tumors. JCO Precis Oncol. 2017. 10.1200/PO.17.00146.
    
    

39. 39.Endris V, Buchhalter I, Allgäuer M, Rempel E, Lier A, Volckmar A-L, et al. Measurement of tumor mutational burden (TMB) in routine molecular diagnostics: *in-silico* and real-life analysis of three larger gene panels. Int J Cancer. 2019. 10.1002/ijc.32002.
    
    

40. 40.Zhang S, So AS, Kaplan S, KK M. Comprehensive evaluation of Illumina’s TruSight® tumor 170 panel to estimate tumor mutational burden. Cancer Res. 2017. 10.1158/1538-7445.AM2017-5358.
    
    

41. 41.Chaudhary R, Quagliata L, Martin JP, Alborelli I, Cyanam D, Mittal V, et al. A scalable solution for tumor mutational burden from formalin- fixed , paraffin-embedded samples using the Oncomine tumor mutation load assay. TLCR. 2018:1–15.
    
    

42. 42. Rizvi H , Sanchez-Vega F , La K , Chatila W , Jonsson P , Halpenny D , et al . Molecular determinants of response to anti-programmed cell death (PD)-1 and anti-programmed death-ligand 1 (PD-L1) blockade in patients with non-small-cell lung cancer profiled with targeted next-generation sequencing. J Clin Oncol. 2018;36:7 633–641. pmcid:6075848 [doi:10.1200/JCO.2017.75.3384](http://dx.doi.org/10.1200/JCO.2017.75.3384) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1200/JCO.2017.75.3384&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

43. 43. Samstein Robert M. , Lee Chung-Han , Shoushtari Alexander N. , Hellmann Matthew D. , Shen Ronglai , Janjigian Yelena Y. , Barron David A. , Zehir Ahmet , Jordan Emmet J. , Omuro Antonio , Kaley Thomas J. , Kendall Sviatoslav M. , Motzer Robert J. , Hakimi A. Ari , Voss Martin H. , Russo Paul , Rosenberg Jonathan , Iyer Gopa , Bochner Bernard H. , Bajorin Dean F. , Al-Ahmadie Hikmat A. , Chaft Jamie E. , Rudin Charles M. , Riely Gregory J. , Baxi Shrujal , Ho Alan L. , Wong Richard J. , Pfister David G. , Wolchok Jedd D. , Barker Christopher A. , Gutin Philip H. , Brennan Cameron W. , Tabar Viviane , Mellinghoff Ingo K. , DeAngelis Lisa M. , Ariyan Charlotte E. , Lee Nancy , Tap William D. , Gounder Mrinal M. , D’Angelo Sandra P. , Saltz Leonard , Stadler Zsofia K. , Scher Howard I. , Baselga Jose , Razavi Pedram , Klebanoff Christopher A. , Yaeger Rona , Segal Neil H. , Ku Geoffrey Y. , DeMatteo Ronald P. , Ladanyi Marc , Rizvi Naiyer A. , Berger Michael F. , Riaz Nadeem , Solit David B. , Chan Timothy A. , Morris Luc G. T. . Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nature Genetics. 2019;51:2 202–206. [doi:10.1038/s41588-018-0312-8](http://dx.doi.org/10.1038/s41588-018-0312-8) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1038/s41588-018-0312-8&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=30643254&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

44. 44.Khagi Y, Goodman AM, Daniels GA, Patel SP, Sacco AG, Randall JM, et al. Hypermutated Circulating Tumor DNA : Correlation with Response to Checkpoint Inhibitor – Based Immunotherapy. 2017;5729–5737.
    
    

45. 45. Hellmann Matthew D. , Callahan Margaret K. , Awad Mark M. , Calvo Emiliano , Ascierto Paolo A. , Atmaca Akin , Rizvi Naiyer A. , Hirsch Fred R. , Selvaggi Giovanni , Szustakowski Joseph D. , Sasson Ariella , Golhar Ryan , Vitazka Patrik , Chang Han , Geese William J. , Antonia Scott J. . Tumor Mutational Burden and Efficacy of Nivolumab Monotherapy and in Combination with Ipilimumab in Small-Cell Lung Cancer. Cancer Cell. 2018;33:5 853-861.e46750707
    
    

46. 46. Qiu Ping , Pang Ling , Arreaza Gladys , Maguire Maureen , Chang Ken , Marton Matthew , Levitan Diane . Data Interoperability of Whole Exome Sequencing (WES) Based Mutational Burden Estimates from Different Laboratories. International Journal of Molecular Sciences. 2016;17:5 651. pmcid:4881477 [doi:10.3390/ijms17050651](http://dx.doi.org/10.3390/ijms17050651) 
    
    

47. 47.Mola N, Schu M, Stiegelmeyer S, Jones W, Weigman V. Tumor mutational burden: guidelines for derivation and robustness of measurement. Cancer Res. 2018. 10.1158/1538-7445.AM2018-2250.
    
    

48. 48. Zook JM , Chapman B , Wang J , Mittelman D , Hofmann O , Hide W , et al . Integrating human sequence data sets provides a resource of benchmark SNP and indel genotype calls. Nat Biotechnol. 2014;32:3 246–251. [doi:10.1038/nbt.2835](http://dx.doi.org/10.1038/nbt.2835) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1038/nbt.2835&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=24531798&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

49. 49.Wang Z, Dua J, Cai S, Han M, Dong H, Zhao J, et al. Assessment of blood tumor mutational burden as a potential biomarker for immunotherapy in patients with non–small cell lung cancer with use of a next-generation sequencing cancer gene panel. JAMA Oncol. 2019:1–7.
    
    

50. 50. Hellmann Matthew D. , Ciuleanu Tudor-Eliade , Pluzanski Adam , Lee Jong Seok , Otterson Gregory A. , Audigier-Valette Clarisse , Minenza Elisa , Linardou Helena , Burgers Sjaak , Salman Pamela , Borghaei Hossein , Ramalingam Suresh S. , Brahmer Julie , Reck Martin , O’Byrne Kenneth J. , Geese William J. , Green George , Chang Han , Szustakowski Joseph , Bhagavatheeswaran Prabhu , Healey Diane , Fu Yali , Nathan Faith , Paz-Ares Luis . Nivolumab plus Ipilimumab in Lung Cancer with a High Tumor Mutational Burden. New England Journal of Medicine. 2018;378:22 2093–2104. [doi:10.1056/NEJMoa1801946](http://dx.doi.org/10.1056/NEJMoa1801946) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1056/NEJMoa1801946&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

51. 51.Ramalingam SS. Tumor mutational burden (TMB) as a biomarker for clinical benefit from dual immune checkpoint blockade with nivolumab (nivo) + ipilimumab (ipi) in first-line (1L) non-small cell lung cancer (NSCLC): identification of TMB cutoff from CheckMate 568. Cancer Res. 2018. 10.1158/1538-7445.AM2018-CT078.
    
    

52. 52. Chen Yu-Pei , Zhang Yu , Lv Jia-Wei , Li Ying-Qin , Wang Ya-Qin , He Qing-Mei , Yang Xiao-Jing , Sun Ying , Mao Yan-Ping , Yun Jing-Ping , Liu Na , Ma Jun . Genomic Analysis of Tumor Microenvironment Immune Types across 14 Solid Cancer Types: Immunotherapeutic Implications. Theranostics. 2017;7:14 3585–3594. pmcid:5596445 [doi:10.7150/thno.21471](http://dx.doi.org/10.7150/thno.21471) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.7150/thno.21471&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=28912897&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

53. 53.Fabrizio D, Chen S-J, Xie M, Chen W, Quinn KJ, Zhao C, et al. In silico assessment of variation in TMB quantification across diagnostic platforms: phase 1 of the friends of Cancer research harmonization project. J Immunother Cancer. 2018. 10.1186/s40425-018-0434-7.
    
    

54. 54.Stenzinger A, Allen J, Maas J, Stewart M, Merino D, Dietel M. Tumor mutational burden (TMB) standardization initiative: establishing a consistent methodology for TMB measurement in clinical samples. Ann Oncol. 2018. 10.1093/annonc/mdy269.139.
    
    

55. 55. Murtaza Muhammed , Dawson Sarah-Jane , Tsui Dana W. Y. , Gale Davina , Forshew Tim , Piskorz Anna M. , Parkinson Christine , Chin Suet-Feung , Kingsbury Zoya , Wong Alvin S. C. , Marass Francesco , Humphray Sean , Hadfield James , Bentley David , Chin Tan Min , Brenton James D. , Caldas Carlos , Rosenfeld Nitzan . Non-invasive analysis of acquired resistance to cancer therapy by sequencing of plasma DNA. Nature. 2013;497:7447 108–112. [doi:10.1038/nature12065](http://dx.doi.org/10.1038/nature12065) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1038/nature12065&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=23563269&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 
    
    [Web of Science](http://jitc.bmj.com/lookup/external-ref?access_num=000318221500043&link_type=ISI) 

56. 56.Murtaza M, Dawson SJ, Pogrebniak K, Rueda OM, Provenzano E, Grant J, et al. Multifocal clonal evolution characterized using circulating tumour DNA in a case of metastatic breast cancer. Nat Commun. 2015. 10.1038/ncomms9760.
    
    

57. 57. Butler TM , Johnson-Camacho K , Peto M , Wang NJ , Macey TA , Korkola JE , et al . Exome sequencing of cell-free DNA from metastatic cancer patients identifies clinically actionable mutations distinct from primary disease. PLoS One. 2015;10:8 1–14. [doi:10.1371/journal.pone.0136407](http://dx.doi.org/10.1371/journal.pone.0136407) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1371/journal.pone.0131080&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=25919488&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

58. 58. Klevebring Daniel , Neiman Mårten , Sundling Simon , Eriksson Louise , Darai Ramqvist Eva , Celebioglu Fuat , Czene Kamila , Hall Per , Egevad Lars , Grönberg Henrik , Lindberg Johan . Evaluation of Exome Sequencing to Estimate Tumor Burden in Plasma. PLoS ONE. 2014;9:8 e104417. pmcid:4136786 [doi:10.1371/journal.pone.0104417](http://dx.doi.org/10.1371/journal.pone.0104417) 
    
    

59. 59. Chan K. C. A. , Jiang P. , Zheng Y. W. L. , Liao G. J. W. , Sun H. , Wong J. , Siu S. S. N. , Chan W. C. , Chan S. L. , Chan A. T. C. , Lai P. B. S. , Chiu R. W. K. , Lo Y. M. D. . Cancer Genome Scanning in Plasma: Detection of Tumor-Associated Copy Number Aberrations, Single-Nucleotide Variants, and Tumoral Heterogeneity by Massively Parallel Sequencing. Clinical Chemistry. 2012;59:1 211–224. [doi:10.1373/clinchem.2012.196014](http://dx.doi.org/10.1373/clinchem.2012.196014) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=23065472&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

60. 60. Koeppel Florence , Blanchard Steven , Jovelet Cécile , Genin Bérengère , Marcaillou Charles , Martin Emmanuel , Rouleau Etienne , Solary Eric , Soria Jean-Charles , André Fabrice , Lacroix Ludovic . Whole exome sequencing for determination of tumor mutation load in liquid biopsy from advanced cancer patients. PLOS ONE. 2017;12:11 e0188174. pmcid:5697854 [doi:10.1371/journal.pone.0188174](http://dx.doi.org/10.1371/journal.pone.0188174) 
    
    

61. 61. Gandara DR , Paul SM , Kowanetz M , Schleifman E , Zou W , Li Y , et al . Blood-based tumor mutational burden as a predictor of clinical benefit in non-small-cell lung cancer patients treated with atezolizumab. Nat Med. 2018;24:September 1441 [doi:10.1038/s41591-018-0134-3](http://dx.doi.org/10.1038/s41591-018-0134-3) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

62. 62.Quinn K, Helman E, Nance T, Artieri C, Yen J, Zhao J, et al. Development and analytical validation of a plasma-based tumor mutational burden (TMB) score from next-generation sequencing panels. Ann Oncol. 2018; doi: 0.1093/annonc/mdy269.
    
    

63. 63.Lyu G, Yeh Y, Yeh Y, Wang Y. Mutation load estimation model as a predictor of the response to cancer immunotherapy. npj Genomic Med. 2018. 10.1038/s41525-018-0051-x.
    
    

64. 64.Chan TA, Yarchoan M, Jaffee E, Swanton C, Quezada SA, Stenzinger A, et al. Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic. Ann Oncol. 2018. 10.1093/annonc/mdy495/5160130%0A.
    
    

65. 65. McGranahan N , Furness AJS , Rosenthal R , Ramskov S , Lyngaa R , Saini SK , et al . Clonal neoantigens elicit T cell immunoreactivity and sensitivity to immune checkpoint bloackade. Science (80- ). 2016;351:6280 1463–1470. pmcid:4984254 [doi:10.1126/science.aaf1490](http://dx.doi.org/10.1126/science.aaf1490) 
    
    [Abstract/FREE Full Text](http://jitc.bmj.com/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNTEvNjI4MC8xNDYzIjtzOjQ6ImF0b20iO3M6MTg6Ii9qaXRjLzcvMS8xODMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

66. 66. Kazmi SM . A retrospective analysis to evaluate prevalence and correlation between PD-L1 score and tumor mutational burden (TMB) levels in patients with solid tumor malignancies. Ann Oncol. 2017;28:Suppl 11 xi6–xi29. 
    
    

67. 67.Salem ME, Puccini A, Grothey A, Raghavan D, Goldberg RM, Xiu J, et al. Landscape of Tumor Mutation Load , Mismatch Repair Deficiency, and PD-L1 Expression in a Large Patient Cohort of Gastrointestinal Cancers. 2018;805–813.
    
    

68. 68. Morrison C , Pabla S , Conroy JM , Nesline MK , Glenn ST , Dressman D , et al . Predicting response to checkpoint inhibitors in melanoma beyond PD-L1 and mutational burden. J Immunother Cancer.. 2018;6:1 1–12. [doi:10.1186/s40425-018-0344-8](http://dx.doi.org/10.1186/s40425-018-0344-8) 
    
    [FREE Full Text](http://jitc.bmj.com/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaml0YyI7czo1OiJyZXNpZCI7czoxMToiNi9TdXBwbF8yLzEiO3M6NDoiYXRvbSI7czoxODoiL2ppdGMvNy8xLzE4My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

69. 69. Xiao Wenjing , Du Nan , Huang Taoyuan , Guo Jinan , Mo Xingkui , Yuan Tao , Chen Yong , Ye Ting , Xu Chunwei , Wang Wenxian , Wang Guoqiang , Cai Shangli , Chen Jing . TP53 Mutation as Potential Negative Predictor for Response of Anti-CTLA-4 Therapy in Metastatic Melanoma. EBioMedicine. 2018;32:119–124. pmcid:6020711 [doi:10.1016/j.ebiom.2018.05.019](http://dx.doi.org/10.1016/j.ebiom.2018.05.019) 
    
    

70. 70.Ock C-Y, Hwang J-E, Keam B, Kim S-B, Shim J-J, Jang H-J, et al. Genomic landscape associated with potential response to anti-CTLA-4 treatment in cancers. Nat Commun. 2017. 10.1038/s41467-017-01018-0.
    
    

71. 71.Kirchhoff T, Ferguson R, Simpson D, Kazlow E, Martinez C, Vogelsang M, et al. Germline determinants of immune related adverse events (irAEs) in melanoma immunotherapy response. Ann Oncol. 2017. 10.1093/annonc/mdx376.
    
    

72. 72. Hasan Ali O , Berner F , Bomze D , Fässler M , Diem S , Cozzio A , et al . Human leukocyte antigen variation is associated with adverse events of checkpoint inhibitors. Eur J Cancer. 2019;107:8–14. [doi:10.1016/j.ejca.2018.11.009](http://dx.doi.org/10.1016/j.ejca.2018.11.009) 
    
    

73. 73. Zehir Ahmet , Benayed Ryma , Shah Ronak H , Syed Aijazuddin , Middha Sumit , Kim Hyunjae R , Srinivasan Preethi , Gao Jianjiong , Chakravarty Debyani , Devlin Sean M , Hellmann Matthew D , Barron David A , Schram Alison M , Hameed Meera , Dogan Snjezana , Ross Dara S , Hechtman Jaclyn F , DeLair Deborah F , Yao JinJuan , Mandelker Diana L , Cheng Donavan T , Chandramohan Raghu , Mohanty Abhinita S , Ptashkin Ryan N , Jayakumaran Gowtham , Prasad Meera , Syed Mustafa H , Rema Anoop Balakrishnan , Liu Zhen Y , Nafa Khedoudja , Borsu Laetitia , Sadowska Justyna , Casanova Jacklyn , Bacares Ruben , Kiecka Iwona J , Razumova Anna , Son Julie B , Stewart Lisa , Baldi Tessara , Mullaney Kerry A , Al-Ahmadie Hikmat , Vakiani Efsevia , Abeshouse Adam A , Penson Alexander V , Jonsson Philip , Camacho Niedzica , Chang Matthew T , Won Helen H , Gross Benjamin E , Kundra Ritika , Heins Zachary J , Chen Hsiao-Wei , Phillips Sarah , Zhang Hongxin , Wang Jiaojiao , Ochoa Angelica , Wills Jonathan , Eubank Michael , Thomas Stacy B , Gardos Stuart M , Reales Dalicia N , Galle Jesse , Durany Robert , Cambria Roy , Abida Wassim , Cercek Andrea , Feldman Darren R , Gounder Mrinal M , Hakimi A Ari , Harding James J , Iyer Gopa , Janjigian Yelena Y , Jordan Emmet J , Kelly Ciara M , Lowery Maeve A , Morris Luc G T , Omuro Antonio M , Raj Nitya , Razavi Pedram , Shoushtari Alexander N , Shukla Neerav , Soumerai Tara E , Varghese Anna M , Yaeger Rona , Coleman Jonathan , Bochner Bernard , Riely Gregory J , Saltz Leonard B , Scher Howard I , Sabbatini Paul J , Robson Mark E , Klimstra David S , Taylor Barry S , Baselga Jose , Schultz Nikolaus , Hyman David M , Arcila Maria E , Solit David B , Ladanyi Marc , Berger Michael F . Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10,000 patients. Nature Medicine. 2017;23:6 703–713. pmcid:5461196 [doi:10.1038/nm.4333](http://dx.doi.org/10.1038/nm.4333) 
    
    [CrossRef](http://jitc.bmj.com/lookup/external-ref?access_num=10.1038/nm.4333&link_type=DOI) 
    
    [PubMed](http://jitc.bmj.com/lookup/external-ref?access_num=28481359&link_type=MED&atom=%2Fjitc%2F7%2F1%2F183.atom) 

74. 74. Sun James X. , He Yuting , Sanford Eric , Montesion Meagan , Frampton Garrett M. , Vignot Stéphane , Soria Jean-Charles , Ross Jeffrey S. , Miller Vincent A. , Stephens Phil J. , Lipson Doron , Yelensky Roman . A computational approach to distinguish somatic vs. germline origin of genomic alterations from deep sequencing of cancer specimens without a matched normal. PLOS Computational Biology. 2018;14:2 e1005965. pmcid:5832436 [doi:10.1371/journal.pcbi.1005965](http://dx.doi.org/10.1371/journal.pcbi.1005965)