Skip to main content

Mutational signature analysis in non-small cell lung cancer patients with a high tumor mutational burden



Lung cancer is the leading cause of cancer death worldwide. With the growing number of targeted therapies and the introduction of immuno-oncology (IO), personalized medicine has become standard of care in patients with metastatic disease. The development of predictive and prognostic biomarkers is of great importance. Mutational signatures harbor potential clinical value as predictors of therapy response in cancer. Here we set out to investigate particular mutational processes by assessing mutational signatures and associations with clinical features, tumor mutational burden (TMB) and targetable mutations.


In this retrospective study, we studied tumor DNA from patients with non-small cell lung cancer (NSCLC) irrespective of stage. The samples were sequenced using a 2 megabase (Mb) gene panel. On each sample TMB was determined and defined as the total number of single nucleotide mutations per Mb (mut/Mb) including non-synonymous mutations. Mutational signature profiling was performed on tumor samples in which at least 30 somatic single base substitutions (SBS) were detected.


In total 195 samples were sequenced. Median total TMB was 10.3 mut/Mb (range 0–109.3). Mutational signatures were evaluated in 76 tumor samples (39%; median TMB 15.2 mut/Mb). SBS signature 4 (SBS4), associated with tobacco smoking, was prominently present in 25 of 76 samples (33%). SBS2 and/or SBS13, both associated with activity of the AID/APOBEC family of cytidine deaminases, were observed in 11 of 76 samples (14%). SBS4 was significantly more present in early stages (I and II) versus advanced stages (III and IV; P = .005).


In a large proportion of NSCLC patients tissue panel sequencing with a 2 Mb panel can be used to determine the mutational signatures. In general, mutational signature SBS4 was more often found in early versus advanced stages of NSCLC. Further studies are needed to determine the clinical utility of mutational signature analyses.


Lung cancer is the leading cause of cancer-related death in many countries[1]. Its high mortality has urged major efforts to optimize the treatment of lung cancer have been made that have resulted in targeted therapies and immunotherapy. These treatments are based on the presence or absence of specific predictive and prognostic biomarkers. In the last decade several biomarkers predicting immune checkpoint blockade outcomes have been discovered. One of these biomarkers is the expression of PD-L1 on tumor cells, indicating eligibility for immune checkpoint inhibitor therapy [2, 3]. However, PD-L1 expression as a predictive biomarker for response to immune checkpoint blockade (ICB) is fairly unreliable due to dynamic and heterogeneous expression in the tumor microenvironment, divergent assay interpretation and lack of PD-L1 platform standardization [4,5,6,7,8]

The first FDA-approved tumor type-agnostic biomarker for immunotherapy is microsatellite instability (MSI) [9]. MSI is caused by inactivation of the mismatch repair (MMR) machinery, resulting in the accumulation of DNA replication errors in repetitive sequences or microsatellites. Alternatively, TMB is suggested as a biomarker, as the number of mutations observed in a tumor seem to correlate with clinical outcome and effectiveness of immunotherapy [10,11,12]. However, TMB is not yet approved as a predictive biomarker for NSCLC as recent studies demonstrated varying results regarding therapeutic benefits [13, 14].

A relative new strategy in the search for new biomarkers is to study the molecular processes in the cancer cell that cause a specific pattern of mutations; a so-called mutational signature. Somatic mutations in a cancer genome are the cumulative result of mutational processes that started since embryonic development [15]. Different mutational processes, such as exposure to UV-light or tobacco smoking, generate a unique combination of mutation types, that can be detected as mutational signatures. Mutational signatures have shown their applicability in cancer diagnosis and prediction of response to treatment [16,17,18]. Here we investigated mutational signature analyses in a retrospective cohort of NSCLC patients who underwent comprehensive genomic tumor profiling.


Sample and data collection

We included 210 samples from patients diagnosed with NSCLC at the Radboud university medical center (Nijmegen, the Netherlands). Samples were sequenced between March 2019 and March 2020. Seven tumor samples were taken before 2019. We included tumor samples with a final diagnosis of NSCLC of any stage either derived directly from the lungs (lobectomy, pneumonectomy, lung or bronchial biopsy and bronchoalveolar lavage) or derived from a metastatic locus (Additional file 1: Table S1). Only one sample per solitary tumor was included. We excluded samples with low tumor cell percentages (< 20%), low median unique coverage (< 60), duplicate tumors and stage 0 disease (carcinoma in situ). A total of 195 tumor samples from 192 patients were used for further analysis. The study eligibility criteria included patients with histological or cytological confirmation of NSCLC (Fig. 1).

Fig. 1
figure 1

Flowchart of the study. NSCLC non-small cell lung cancer, AIS adenocarcinoma in situ, SNV single nucleotide variants. * In one sample molecular sequencing analysis was already performed in the context of the TSO500 validation trial written by Leonie Kroeze et al. [23]

The study was conducted in accordance with the institutional guidelines and regulations from the Radboud university medical center. Written informed consent was obtained for all patients. We obtained the electronic medical records of all these patients and extracted age, gender, smoking status, stage, Eastern Cooperative Oncology Group (ECOG) performance score at baseline and treatment modality.

Tumor sequencing and analyses

All tumor samples were subjected to sequencing analysis using TSO500 (Illumina), a next-generation sequencing panel containing 523 cancer related genes (total size: 2 Mb), performed either as part of standard care in patients with advanced NSCLC or for specific study purposes (for example the LEMA trial) [19]. Sequencing libraries were prepared using the hybrid capture-based TSO500 library preparation kit following the manufacturer’s protocol. After quantification, normalization and pooling, the libraries were sequenced on a NextSeq 500 (Illumina), with 10 libraries sequenced per run (NextSeq high output). The sequence data were processed and analyzed by the TruSight Oncology 500 Local App version 1.3 or 2.0 (Illumina). Analyses of single nucleotide variants (SNVs), multiple nucleotide variants (MNVs), copy number variants (CNVs), microsatellite instability (MSI) and tumor mutational burden (TMB) were performed. A tumor was considered MSI positive when at least 25% of accessible MSI sites are unstable. TMB was defined as the number of synonymous and nonsynonymous mutations with a variant allele frequency of at least 5% per Mb of sequence. We used 10 mut/Mb as cut-off value based on large-base clinical studies [20, 21]. Patient characteristics were investigated according to TMB status (Additional file 2: Table S2). For tailored reporting purposes, a gene panel of 15 genes relevant for lung cancer was investigated for relevant actionable mutations (Additional file 2: Table S3). Tumor samples with at least 30 somatic single base substitutions were investigated by means of mutational signature analysis using COSMIC mutational signatures v3 (n = 76) as described by Kroeze et al. [22]23.


Differences between groups were calculated using the Fisher's exact test. A p-value of < 0.05 is considered significant.


Clinical and histopathologic features of NSCLC

In total, 195 NSCLC samples eligible for inclusion in this study were molecularly profiled (Fig. 1). Mean age of diagnosis was 67 years (range 34–88) and 54.8% of the patients were male. Most tumors were stage IV (n = 98; 50.3%), followed by stage III (n = 44; 22.5%), stage I (n = 30; 15.4%) and stage II (n = 23; 11.8%). Review of histopathology showed that the majority of tumors represented adenocarcinoma (AC, n = 110; 56.4%), followed by squamous cell carcinoma (SCC, n = 42; 21.5%) and large cell neuroendocrine carcinoma (LCNEC, n = 11; 5.6%). A history of smoking, active or former, was reported for 98.9% of cases (n = 193). For 84% of the tumors PD-L1 status was investigated (n = 164). PD-L1 expression was negative in 80 tumors (48.8%), intermediate (1–50% positive stained tumor cells) in 32 tumors (19.5%) and high (above 50% positive stained tumor cells) in 52 tumors (31.7%) (Additional file 2: Table S2).

High tumor mutational burden is associated with lower tumor stages and specific driver gene mutations

The median unique exon coverage was 474 (range 64–851). MSI could be assessed in 168 tumors (86.2%), but none showed MSI (Additional file 1: Table S1). The median overall total TMB was 10.3 (range 0.0–109.3; Additional file 2: Figure S2A). Median total TMB varied among different disease stages. Stage IV tumors presented with a lower TMB (median TMB 9.5) than tumors in earlier stages of disease (median TMB 10.2; P = 0.013, 13.9; P = 0.003, and 13.4; P = 0.051, for stage I, II and III, respectively; Additional file 2: Table S2; Additional file 2: Figure S2B). Moreover, median total TMB was higher in tumors with PD-L1 expression (> 1% positive tumor cells) than tumors that did not express PD-L1 (median TMB 13.0 vs 9.5, respectively; P = 0.004), of which tumors with high PD-L1 expression (> 50 positive tumor cells) had the highest TMB (median TMB 14.2; P = 0.001; Additional file 2: Figure S2C). Age appeared to be lower in the high TMB group (> 10 mut/Mb) (P = 0.025). No correlations between TMB and sex or histopathological subtypes were observed (Additional file 2: Table S2).

Furthermore, we investigated somatic driver mutations with clinical relevance for NSCLC patients (Additional file 2: Table S3). TP53 mutations were more frequent in the TMB high group versus TMB low group (P = 0.007). We also noted a higher frequency of STK11 mutations in the TMB high group (P = 0.012).

Mutational signature SBS4 is associated with lower tumor stages

We analyzed the mutation types and mutational signatures of tumor samples with at least 30 SBS (n = 76; Table 1). The median TMB of this set of tumors was 15.2 mut/Mb. The mutation types observed in the tumors were mostly C > A and C > T mutations (Fig. 2). By refitting of mutational signatures (COSMIC v3) we could assign the majority of somatic mutations detected in the tumor samples to a known mutational signature (average cosine 0.77; range 0.50–0.96; Additional file 1: Table S1). SBS signature 4 (SBS4), associated with tobacco smoking, contributed to the mutation spectrum in 25 of 76 tumors (33%) with a relative contribution of at least 20%. All patients showing SBS4 had a smoking history, among which one patient only reported passive smoking during decennia. SBS signatures 2 and/or 13 (SBS2/13), both associated with activity of the AID/APOBEC family of cytidine deaminases, cumulatively contributed to the mutation spectrum in 11 of 76 tumors (14%). SBS signature 29 (SBS29), associated with tobacco chewing, contributed to the mutation spectrum in 5 of 76 tumors (7%). The contribution of each of these mutational signatures was almost mutually exclusive, with the exception of one sample in which both SBS4 and SBS2/13 were detected with a relative contribution of more than 20% (Additional file 1: Table S1; Fig. 2). Furthermore, SBS signature 39 (SBS39), a signature with unknown etiology, likely contributed to the mutation spectrum in 11 of 76 tumors (14%). In one tumor (T213, TMB 26 mut/Mb) the mutation spectrum was almost solely explained by signatures associated with exposure to UV light (SBS7a and SBS7b; Fig. 2). Clinical and pathological revision revealed that this patient presented with an ulcerative skin lesion with malignant properties, issued as a lung metastatic lesion. Histopathological analysis performed on endobronchial biopsies pointed in the direction of a keratin positive NSCLC not otherwise specified. However, we found pathological mutations in the PTPN11, NF1, CUX1, IKZF1 gene and TERT promoter, which better fit a diagnosis of a melanoma than of lung cancer.

Table 1 Patient characteristics, according to mutational signature status
Fig. 2
figure 2

Clinicopathological and molecular features of lung tumor samples with ≥ 30 somatic SBS in genes covered by the TSO500 panel. Samples are sorted from left-to-right based on high-to-low total TMB. Every column represents one tumor sample. Numbers indicated in purple and green indicate two independent tumors (not clonally related) from the same individual. Only lung cancer related genes in which 2 or more somatic mutations or amplifications were detected are included. TMB tumor mutational burden, Non-syn non-synonymous, MT mutation, NOS not otherwise specified, LCNEC Large cell neuroendocrine carcinoma, SBS single base substitution

Next, we investigated specific mutational signatures in relation to tumor characteristics and driver gene mutations. SBS4 was more present in early stages (I and II) versus advanced stages (III and IV; P = 0.005). Samples harboring SBS4 had a higher median TMB than samples with SBS2/13, SBS29 or other contributions (Additional file 2: Figure S1). PD-L1 status did not differ among SBS4, SBS2/13 and SBS29 (Additional file 2: Table S3). TP53 mutations were more often seen in tumors with a high relative contribution of signature SBS4 compared to tumors in which the contribution of SBS4 was low or absent (P = 0.002) (Additional file 2: Table S4). Other driver mutations were not associated with the presence of SBS4, SBS2/13 or SBS29 in a tumor. No difference was noticed regarding samples with SBS4 versus all other samples regarding sampling locus or histological subtype (AC or SCC).


Here, we investigated tumor mutational burden and mutational signatures in a clinical cohort of NSCLC. Using a 2 Mb gene panel in 39% of tumors sufficient single base substitutions were detected to describe the mutational signatures. SBS4 and SBS2/13 were detected in 33% and 14%. SBS4 appears to be more common in early stages of NSCLC. In addition, mutational signature analyses led to the re-classification of a NSCLC as a metastasis of a tumor that likely originated from the skin.

By mutational signature refitting in samples with a high TMB we observed a major contribution of SBS4, SBS2/13 and SBS29 as one would expect in NSCLC [24]. A major proportion of somatic mutations is known to be formed by exogenous exposures like tobacco smoking in lung cancer. Furthermore, genomic alterations introduced by smoking persist for many years after smoking cessation [25]. Previously, Rizvi and colleagues found a link between signature SBS4 and progression-free survival in patients with NSCLC who received Pembrolizumab. These data have not yet been confirmed in other studies but reveal the potential of qualitative analysis of mutation patterns in relation to treatment outcome. The apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) related signatures SBS2 and SBS13 are frequently observed in our population. The presence of these mutational signatures is markedly associated with response to immune checkpoint blockade (ICB) therapy [26]. It is suggested that the cellular machinery underlying SBS2 and SBS13 is activated by tobacco smoke via direct or indirect pathways [27]. The clinical implications of our study need further assessment in larger series of clinical lung cancer samples in which TMB and mutational signature analyses can be combined with long-term follow-up data on treatment and survival.

In one tumor sample SBS7a and SBS7b could fully explain the mutation spectrum, which suggests that this tumor originates from sun exposed skin [24]. SBS7 mutational signatures are normally seen in skin cancers such as melanoma [28]. Analyzing mutational signatures in NSCLC can as such contribute to assess the primary tumor site of the malignancy.

This is the first study that describes higher TMB levels in early NSCLC stages versus stage IV. In a systematic review of Willes regarding TMB and lung cancer, four publications assessed lung cancer stage, all reporting no significant association with TMB [29]. Within this review two publications referred to patients with SCC or LCNEC/SCLC only [30, 31], whereas the remaining two articles referred to patients with AC [32, 33]. The clonal structure of a tumor varies considerably between primary and metastatic sites, with higher rates of monoclonal structures recorded in metastases due to clonal selection [34]. As such, TMB can be measured from a primary or metastatic tumor sample, causing systematic bias in TMB values.

To our knowledge, this is the first study showing signature SBS4 is more frequently detected in early versus advanced stages of NSCLC. An explanation for this finding could be the higher risk of cardiovascular and pulmonary disease due to tobacco smoking, which often leads to relatively early imaging diagnostics.

Since Alexandrov and Stratton uncovered and catalogued mutational signatures, our understanding of the mutational processes that cause somatic mutations is markedly expanded [15]. However, the understanding of mutational processes in most cancer types is remarkably limited. Refitting of mutational signatures in clinical samples may provide additional insights into the processes underlying cancer development. However, as we have seen in our study this is only feasible when sufficient somatic mutations are detected in a tumor sample. We have applied the criterium of 30 SBSs and a minimal relative contribution of 20% to conclude that a mutational mechanism has been active in a tumor. Using a next gen sequencing panel covering 2 Mb in our population 39% of tumor samples could be evaluated for mutational signature analysis. When whole exome sequencing (WES) or whole genome sequencing (WGS) can be applied, more mutations are detected per tumor resulting in a higher percentage of samples suitable for mutational signature analysis. However, WES or WGS are not yet part of standard of care in cancer diagnostics. Additional studies are needed to assess the clinical utility of mutational signature analyses and to set guidelines to further translate findings of tumor-specific mutational processes into clinical practice.

In conclusion, this is the first report that systematically studies mutational signatures in a set of tumors sequenced in routine clinical practice in NSCLC. With a panel covering about 500 genes mutational signatures can be determined in a significant proportion of NSCLC. Interestingly, mutational signature SBS4 was more common in early versus advanced stages of NSCLC. Furthermore, mutational signature profiling may facilitate the diagnosis of the primary tumor site in a clinical setting. However, studies of greater magnitude and/or with other sequencing strategies using WES/WGS are needed to assess the clinical utility and to translate our findings into clinical practice.

Availability of data and materials

All patients are known in the outpatient clinic of Pulmonology at the Radboud university medical center. Most of the included patients have provided informed consent for Lung cancer Early Molecular Assessment trial (LEMA) ( NCT02894853) or Lung cancer biobank (NCT01084785) or gave otherwise permission for using body specimens for research purposes. Missing informed consents for living patients have been collected.





Activation-induced deaminase/apolipoprotein B mRNA-editing enzyme, catalytic polypeptide


Copy number variants


Catalogue of somatic mutations in cancer


Eastern cooperative oncology group


Electronic case report form


Food and drug administration


Immune checkpoint blockade




Large cell neuroendocrine carcinoma




Multiple nucleotide variants


Microsatellite instability


Mutations per megabase


Non-small cell lung cancer


Programmed cell death- ligand 1


Single base substitute


Squamous cell carcinoma


Small cell lung cancer


Single nucleotide variants


Tumor mutational burden


Whole exome sequencing


Whole genome sequencing


  1. de Koning HJ, van der Aalst CM, de Jong PA, Scholten ET, Nackaerts K, Heuvelmans MA, et al. Reduced lung-cancer mortality with volume CT screening in a randomized trial. NEJM N Engl J Med. 2020;382:503–13.

    Article  Google Scholar 

  2. Yu H, Boyle TA, Zhou C, Rimm DL, Hirsch FR, et al. PD-L1 expression in lung cancer. J Thorac Oncol. 2016;11:964–75.

    Article  Google Scholar 

  3. Reck M, Rodríguez-Abreu D, Robinson AG, Hui R, Csőszi T, Fülöp A, et al. Pembrolizumab versus Chemotherapy for PD-L1–positive non–small-cell lung cancer. N Engl J Med. 2016;375:1823–33.

    Article  CAS  Google Scholar 

  4. Patel SP, Kurzrock R. PD-L1 expression as a predictive biomarker in cancer immunotherapy. Mol Cancer Ther. 2015;14:847–56.

    Article  CAS  Google Scholar 

  5. Shukuya T, Carbone DP. Predictive markers for the efficacy of anti-PD-1/PD-L1 antibodies in lung cancer. J Thorac Oncol. 2016;11:976–88.

    Article  Google Scholar 

  6. Schoenfeld AJ, Rizvi H, Bandlamudi C, Sauter JL, Travis WD, Rekhtman N, et al. Clinical and molecular correlates of PD-L1 expression in patients with lung adenocarcinoma. Ann Oncol. 2000;31(5):599–608.

    Article  Google Scholar 

  7. Arfe A, Fell G, Alexander B, Awad MM, Rodig SJ, Trippa L, Schoenfeld JD, et al. Meta-analysis of PD-L1 expression as a predictor of survival after checkpoint blockade. Immunotherapy. 2016;4:479–88.

    Google Scholar 

  8. Brito AB, Camandaroba MP, de Lima VC. Anti-PD1 versus anti-PD-L1 immunotherapy in first-line therapy for advanced non-small cell lung cander: A systematic review and meta-analysis. Thorac Cancer. 2021;12(7):1058–66.

    Article  CAS  Google Scholar 

  9. Chang L, Chang M, Chang HM, Chang F. Microsatellite instability: a predictive biomarker for cancer immunotherapy. Appl Immunohistochem Mol Morphol. 2018;26(2):e15–21.

    Article  CAS  Google Scholar 

  10. 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–8.

    Article  CAS  Google Scholar 

  11. Carbone DP, Reck M, Paz-Ares L, Creelan B, Horn L, Steins M, et al. First-line Nivolumab in stage IV or recurrent non–small-cell lung cancer. N Engl J Med. 2017;376:2415–3242.

    Article  CAS  Google Scholar 

  12. Hellmann MD, Paz-Ares L, Bernabe Caro R, Zurawski B, et al. Nivolumab plus Ipilimumab in advanced non-small-cell lung cancer. N Engl J Med. 2019;381:2020–31.

    Article  CAS  Google Scholar 

  13. Langer CJ, Gadgeel SM, Borghaei H, et al. KEYNOTE-021: TMB and outcomes for carboplatin and Pemetrexed with or without Pembrolizumab for non-squamous NSCLC. J Thor Oncol. 2019;17(11):1497–508.

    Google Scholar 

  14. Garassino M, Rodriguez-Abreu D, Gadgeel S, Esteban E, et al. Evaluation of TMB in KEYNOTE-189: Pembrolizumab plus chemotherapy vs placebo plus chemotherapy for nonsquamous NSCLC. J Thor Oncol. 2019;14(10):S216-7.

    Article  Google Scholar 

  15. Alexandrov LB, Stratton MR. Mutational signatures: the patterns of somatic mutations hidden in cancer genomes. Curr Opin Genet Dev. 2014;24:52–60.

    Article  CAS  Google Scholar 

  16. Ma J, Setton J, Lee NY, Riaz N, Powell SN. The therapeutic significance of mutational signatures form DNA repair deficiency in cancer. Nat Commun. 2018;9:3292.

    Article  Google Scholar 

  17. Chen H, Chong W, Teng C, Yao Y, Wang X, Li X. The immune response-related mutational signatures and driver genes in non-small-cell lung cancer. Cancer Sci. 2019;110(8):2348–56.

    Article  CAS  Google Scholar 

  18. Anagnostou V, Niknafs N, Marrone K, Bruhm DC, White JR, Naidoo J, et al. Multimodal genomic features predict outcome of immune checkpoint blockade in non-small-cell lung cancer. Nat Cancer. 2020;1:99–111.

    Article  Google Scholar 

  19. Lung Cancer Early Molecular Assessment Trial.

  20. Ready N, Hellmann MD, Awad MM, Otterson GA, Gutierrez M, Gainor JF, et al. First-line Nivolumab plus Ipilimumab in advanced non–small-cell lung cancer (CheckMate 568). J Clin Oncol. 2019;37:992–1000.

    Article  CAS  Google Scholar 

  21. Hellmann MD, Paz-Ares L, Bernabe Caro R, Zurawski B, et al. Nivolumab plus Ipilimumab in advancded non-small-cell lung cancer (CheckMate 227). N Engl J Med. 2019;381:21.

    Article  Google Scholar 

  22. Forbes SA, Beare D, Boutselakis H, Bamford S, Bindal N, Tate J, Cole CG, Ward S, Dawson E, Ponting L, Stefancsik R. COSMIC: somatic cancer genetics at high-resolution. Nucleic Acids Res. 2017;45(D1):D777-83.

    Article  CAS  Google Scholar 

  23. Kroeze LI, de Voer RM, Kamping EJ, von Rhein D, Jansen EAM, Hermsen MJW, et al. J Mol Diagn. 2020;22(6):757–69.

    Article  CAS  Google Scholar 

  24. Alexandrov LB, Kim J, Haradhvala NJ, Huang MN, Ng AW, Wu Y, Boot A, Covington KR, Gordenin DA, Bergstrom EN, Islam SA. The repertoire of mutational signatures in human cancer. Nature. 2020;578(7793):94–101.

    Article  CAS  Google Scholar 

  25. Nicola PA, Rafee S, Burghel G, Wallace A, Schlecht H, et al. Persistence of smoking signature 4 in the non-small cell lung cancer genome. Cancer Res. 2020;80(16 Suppl):3810.

    Google Scholar 

  26. Wang S, Jia M, Liu X-S. APOBEC3B and APOBEC mutational signature as potential predictive markers for immunotherapy response in non-small cell lung cancer. Oncogene. 2018;37:3924–36.

    Article  CAS  Google Scholar 

  27. Alexandrov LB, Ju YS, Haase K, van Loo P, Martincorena I, Nik-Zainal S, et al. Mutational signatures associated with tobacco smoking in human cancer. Science. 2016;354(6312):618–22.

    Article  CAS  Google Scholar 

  28. Yang C, Sanchez-Vega F, Chang JC, Chatila WK, et al. Lung-only melanoma: UV mutational signature supports origin from occult cutaneous primaries and argues against the concept of primary pulmonary melanoma. Modern Pathol. 2020;33:2244–55.

    Article  CAS  Google Scholar 

  29. Willis C, Fiander M, Tran D, Korytowsky B, Thomas J-M, Calderon F, et al. Tumor mutational burden in lung cancer: a systematic literature review. Oncotarget. 2019;10(61):6604–22.

    Article  Google Scholar 

  30. Choi M, Kadara H, Zhangs J, Parra ER, Rodriguez-Canales J, Gaffney SG, et al. Mutation profiles in early-stage lung squamous cell carcinoma with clinical follow-up and correlation with markers of immune function. Ann Oncol. 2017;28:83–9.

    Article  CAS  Google Scholar 

  31. Xiao D, Li F, Pan H, Liang H, Wu K, He J. Integrative analysis of genomic sequencing data reveals higher prevalence of LRP1B mutations in lung adenocarcinoma patients with COPD. Sci Rep. 2017;7(1):2121.

    Article  Google Scholar 

  32. Kim HS, Lee JH, Nam SJ, Ock C-Y, Moon J-W, Yoo CW, et al. Association of PD-L1 expression with tumor-infiltrating immune cells and mutation burden in high-grade neuroendocrine carcinoma of the lung. J Thorac Oncol. 2018;13(5):636–48.

    Article  Google Scholar 

  33. Xiao D, Pan H, Li F, Wu K, Zhang X, He J. Analysis of ultra-deep targeted sequencing reveals mutation burden is associated with gender and clinical outcome in lung adenocarcinoma. Oncotarget. 2016;7(16):22857–64.

    Article  Google Scholar 

  34. Schnidrig D, Turajlic S, Litchfield K. Tumour mutational burden: primary versus metastatic tissue creates systematic bias. Immuno-Oncol Technol. 2019;4:8–14.

    Article  Google Scholar 

Download references


Not applicable.


The authors received no financial support for the research, authorship, or publication of this article.

Author information

Authors and Affiliations



Data were collected, analyzed and interpreted by the authors. All authors reviewed the manuscript and vouch for accuracy and completeness of data and for the adherence of the study to the protocol. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Richarda M. de Voer or Michel M. van den Heuvel.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the principles of the Declaration of Helsinki (version 8, 19 October 2013).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Table S1.

Clinical, histopathological and molecular features of non-small cell lung cancer samples.

Additional file 2.

Additional figures and tables.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

van den Heuvel, G.R.M., Kroeze, L.I., Ligtenberg, M.J.L. et al. Mutational signature analysis in non-small cell lung cancer patients with a high tumor mutational burden. Respir Res 22, 302 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: