Skip to main content

Tumor Mutational Burden: What It Means for Cancer Morphology

Pierre-Antoine Bannier ·

Tumor mutational burden varies by more than 100-fold across cancer types — yet a single threshold of 10 mutations per megabase determines whether a patient qualifies for pembrolizumab, regardless of where their tumor originated. HistoAtlas data across 6,745 TCGA samples and 21 cancer types reveals that this molecular metric leaves a measurable imprint on the physical architecture of tumors: their proliferation rate, nuclear shape, and immune landscape.

What Is Tumor Mutational Burden and Why Does It Matter?

TMB counts the somatic mutations accumulated in a tumor’s genome, typically reported as mutations per megabase (mut/Mb). A higher mutation count increases the probability that some mutations will produce neoantigens — abnormal peptides presented on the cell surface — that the immune system can recognize [1].

This logic underpinned the FDA’s June 2020 accelerated approval of pembrolizumab for all TMB-high solid tumors (≥10 mut/Mb) — only the fourth time the FDA approved a cancer treatment based on a biomarker rather than the tumor’s organ of origin [2]. In the KEYNOTE-158 trial, TMB-high patients showed a 29% overall response rate, with 57% of responses lasting at least 12 months.

But TMB is not a simple on/off switch. Cancer types differ enormously in their baseline mutational rates. Melanoma and lung cancer accumulate mutations from UV exposure and smoking, respectively, while pediatric cancers and, on average, prostate tumors tend to have low TMB. The question for computational pathology is direct: can we see the consequences of high mutational burden in the tissue itself?

Tumor Mutational Burden Range: What Counts as High or Low?

TMB is a continuous variable, but clinical practice generally divides it into ranges. While exact cutoffs vary by assay and cancer type, a commonly used framework is:

  • Low TMB: fewer than 6 mut/Mb — typical of thyroid carcinoma, prostate adenocarcinoma, and pediatric cancers
  • Intermediate TMB: 6–19 mut/Mb — includes many breast, ovarian, and head and neck cancers
  • High TMB: 20 or more mut/Mb — common in melanoma, lung cancers, and MSI-high colorectal cancers

The FDA’s pembrolizumab threshold of ≥10 mut/Mb sits within the intermediate range, reflecting a pragmatic cutoff rather than a sharp biological boundary [5].

Tumor Mutational Burden Testing

TMB is measured through next-generation sequencing of tumor DNA. The two main approaches are whole-exome sequencing (WES), which surveys all protein-coding regions (~30 Mb), and targeted gene panels that sequence a smaller set of cancer-relevant genes (typically 0.8–2.4 Mb) and extrapolate the total burden. The FDA-approved companion diagnostic for pembrolizumab is FoundationOne CDx, a 324-gene panel. Other widely used panels include MSK-IMPACT (468 genes) and TruSight Oncology 500. Panel-based TMB estimates generally correlate well with WES-derived values, though smaller panels can overestimate or underestimate TMB in individual cases [7].

Tumor Mutational Burden Reshapes Morphology Unevenly Across Cancer Types

Not all cancers show equal sensitivity to mutation-associated morphological change. HistoAtlas quantified the number of statistically significant associations between somatic mutations (across 25 frequently mutated genes) and 40 histomic features, for each of the 21 cancer types in the atlas.

Figure 1. Number of statistically significant mutation–morphology associations per cancer type. Colorectal (COAD), stomach (STAD), and uterine (UCEC) cancers — known for high mutational burden and microsatellite instability — show the most mutation-driven morphological variation. Dashed line at 50 highlights the sharp drop-off after the top four.
Data: HistoAtlas / TCGA (n=6,745). BH-adjusted p < 0.05, Mann-Whitney U with Cliff's delta.

The ranking is telling. Colorectal adenocarcinoma (COAD) leads with 330 significant associations, followed by stomach adenocarcinoma (STAD, 221) and uterine endometrial carcinoma (UCEC, 183). These are precisely the cancer types known for microsatellite instability (MSI) and mismatch repair deficiency — molecular states that drive high TMB [4]. Lung adenocarcinoma (LUAD), where smoking-induced mutations dominate, follows with 108.

Below the top four, a sharp drop-off appears. Bladder, head and neck, and rectal cancers show 34–45 associations each. Low-TMB cancers like prostate (PRAD, 21) and thyroid (THCA, 6) show minimal mutation-associated morphological variation — their tissue architecture is shaped by other forces.

This gradient is not simply a sample-size artifact. BRCA has the largest cohort (n=1,037) yet ranks eighth, while COAD (n=441) leads. The signal reflects genuine biology: cancers with higher mutational burden produce more diverse morphological phenotypes.

High Tumor Mutational Burden Is Associated with Both Proliferation and Immune Infiltration

What exactly changes in a TMB-high tumor’s appearance? HistoAtlas compared TTN-mutant (TMB-high proxy) versus wild-type tumors across the full pan-cancer cohort (n=5,304) using Cliff’s delta effect sizes.

Figure 2. Effect of high mutational burden on tumor morphology and immune infiltration (pan-cancer). Bars show Cliff's delta comparing TTN-mutant versus wild-type tumors — a validated proxy for TMB-high status [2]. Mitotic index shows the largest effect (d = 0.316), followed by multiple immune infiltration metrics. All associations BH-significant (padj < 10-10).
Data: HistoAtlas / TCGA (n=5,304 pan-cancer). Mann-Whitney U test, Cliff's delta effect size.

Mitotic index — the density of actively dividing cells — shows the strongest effect 0.316 Cliff's d (p < 10⁻⁷⁷, n=5,304) . TMB-high tumors show higher mitotic rates. This association likely reflects shared upstream biology: defective DNA repair drives both mutation accumulation and, when cell-cycle regulators are affected, increased proliferation [5].

But proliferation is only half the story. Six immune infiltration metrics are significantly elevated in TMB-high tumors. Peritumoral immune richness (d = 0.278), intratumoral eosinophil density (d = 0.275), and neutrophil density (d = 0.230) all show moderate effects at extraordinary statistical significance (all p < 10⁻⁴⁰). Intratumoral lymphocyte density — the canonical TIL metric — is also elevated (d = 0.113, p < 10⁻¹⁰), though with a smaller effect than innate immune cells.

DNA Repair Pathway Activity Mirrors the Proliferation–TMB Axis

In rapidly proliferating tumors, DNA repair pathway activity is correspondingly elevated — a natural consequence of increased DNA replication requiring more repair. HistoAtlas’s pathway-level Spearman correlations confirm this coupling across cancer types.

LUADLIHCBRCABLCAMESOPAADUCECCESCTHCACOADLUSCMitotic indexMitotic index × LUAD: ρ = 0.54 *0.54Mitotic index × LIHC: ρ = 0.43 *0.43Mitotic index × BRCA: ρ = 0.41 *0.41Mitotic index × BLCA: ρ = 0.38 *0.38Mitotic index × MESO: ρ = 0.34 *0.34Mitotic index × PAAD: ρ = 0.30 *0.30Mitotic index × UCEC: ρ = 0.27 *0.27Mitotic index × CESC: ρ = 0.24 *0.24Mitotic index × THCA: ρ = 0.19 *0.19Mitotic index × COAD: ρ = 0.18 *0.18Mitotic index × LUSC: ρ = 0.15 *0.15PleomorphismPleomorphism × LUAD: ρ = 0.27 *0.27Pleomorphism × LIHC: ρ = 0.23 *0.23Pleomorphism × BRCA: ρ = 0.28 *0.28Pleomorphism × BLCA: ρ = 0.36 *0.36Pleomorphism × MESO: ρ = 0.250.25Pleomorphism × PAAD: ρ = 0.140.14Pleomorphism × UCEC: ρ = 0.090.09Pleomorphism × CESC: ρ = 0.000.00Pleomorphism × THCA: ρ = -0.02-0.02Pleomorphism × COAD: ρ = 0.040.04Pleomorphism × LUSC: ρ = 0.030.03TIL densityTIL density × LUAD: ρ = -0.06-0.06TIL density × LIHC: ρ = 0.120.12TIL density × BRCA: ρ = 0.09 *0.09TIL density × BLCA: ρ = 0.25 *0.25TIL density × MESO: ρ = -0.04-0.04TIL density × PAAD: ρ = 0.120.12TIL density × UCEC: ρ = 0.13 *0.13TIL density × CESC: ρ = -0.02-0.02TIL density × THCA: ρ = 0.15 *0.15TIL density × COAD: ρ = -0.03-0.03TIL density × LUSC: ρ = -0.01-0.01Peritumoral immunePeritumoral immune × LUAD: ρ = 0.19 *0.19Peritumoral immune × LIHC: ρ = 0.19 *0.19Peritumoral immune × BRCA: ρ = 0.14 *0.14Peritumoral immune × BLCA: ρ = 0.100.10Peritumoral immune × MESO: ρ = 0.110.11Peritumoral immune × PAAD: ρ = 0.270.27Peritumoral immune × UCEC: ρ = 0.040.04Peritumoral immune × CESC: ρ = 0.110.11Peritumoral immune × THCA: ρ = 0.19 *0.19Peritumoral immune × COAD: ρ = -0.03-0.03Peritumoral immune × LUSC: ρ = -0.19 *-0.19Stromal lymphocytesStromal lymphocytes × LUAD: ρ = 0.070.07Stromal lymphocytes × LIHC: ρ = 0.040.04Stromal lymphocytes × BRCA: ρ = 0.30 *0.30Stromal lymphocytes × BLCA: ρ = 0.100.10Stromal lymphocytes × MESO: ρ = -0.03-0.03Stromal lymphocytes × PAAD: ρ = 0.090.09Stromal lymphocytes × UCEC: ρ = 0.100.10Stromal lymphocytes × CESC: ρ = 0.040.04Stromal lymphocytes × THCA: ρ = 0.13 *0.13Stromal lymphocytes × COAD: ρ = -0.01-0.01Stromal lymphocytes × LUSC: ρ = 0.050.050.00.10.20.30.40.5Spearman ρ
Figure 3. DNA repair pathway activity correlates with morphological features across cancer types. Heatmap shows Spearman correlations between the DNA_REPAIR pathway score and five histomic features. Bold values indicate BH-adjusted significance (p < 0.05). Mitotic index shows the strongest and most consistent positive correlation across all 11 cancer types (rho up to 0.54 in LUAD). Pleomorphism and immune features show more cancer-type-specific patterns.
Data: HistoAtlas / TCGA. Spearman rank correlation, BH-adjusted p < 0.05.

The DNA_REPAIR pathway score correlates most strongly with mitotic index — and the magnitude varies by cancer type in a biologically coherent pattern. Lung adenocarcinoma (LUAD) shows the strongest correlation (rho = 0.54, p < 10⁻³, n=437), consistent with LUAD’s known proliferative biology. Liver (LIHC, rho = 0.43, n=348) and breast cancer (BRCA, rho = 0.40, n=958) follow. In breast cancer, this is consistent with the high proliferative and mutational burden characteristic of the triple-negative subtype. The correlation is significant across all 11 cancer types tested, from LUAD down to LUSC (rho = 0.14).

Beyond mitotic index, tumor pleomorphism shows a more cancer-type-specific pattern: bladder cancer (BLCA, rho = 0.37) and breast cancer (BRCA, rho = 0.29) show the strongest DNA repair–pleomorphism correlations, while cervical and thyroid cancers show near-zero effects. Immune infiltration features show scattered significance across cancer types — notably BRCA for stromal lymphocytes (rho = 0.30) and BLCA for intratumoral lymphocytes (rho = 0.25). These variable patterns confirm that mutational burden leaves morphological signatures, but the specific features affected depend on the cancer type. This may help explain why deep learning models have shown the ability to predict TMB from H&E-stained tissue sections [6], though cross-institutional reproducibility remains an active research challenge.

What Does “Low Tumor Mutational Burden” Mean for Morphology?

A common clinical question is whether low TMB is “good.” The HistoAtlas data provides a nuanced answer. Low-TMB tumors generate fewer neoantigens, resulting in less immune infiltration and generally a reduced likelihood of immunotherapy response — though individual responses vary. But they also show lower mitotic rates and less nuclear irregularity — meaning the tumor grows in a more ordered, less chaotic fashion.

Thyroid carcinoma (THCA) exemplifies this: only 6 significant mutation-morphology associations across all 25 genes and 40 features. These tumors are morphologically homogeneous and genomically quiet — and while they are less likely to respond to checkpoint inhibitors, they also carry a favorable prognosis for most patients.

The clinical takeaway: TMB is not inherently “good” or “bad.” It defines a biological axis — from genomically quiet, morphologically uniform tumors to mutation-rich, proliferative, immune-infiltrated tumors — and the optimal therapy depends on where a patient’s tumor falls along that axis.

Explore Tumor Mutational Burden Data on HistoAtlas

The data behind this article is fully interactive on HistoAtlas. Explore mutation–morphology associations across all 21 cancer types and 40 histomic features — filter by gene, cancer type, and effect size to identify the strongest signals.

For gene-specific mutation data, visit the individual mutation pages: TP53 mutations, KRAS mutations, BRAF mutations, and EGFR mutations each show how specific driver mutations associate with tumor architecture. The HistoAtlas main explorer lets you visualize all 6,745 slides in a UMAP embedding, colored by mutation status, morphology cluster, or any of the 40 histomic features.

For related analysis methods, see our guides to gene set enrichment analysis in cancer and copy number variation across cancer types.


Frequently Asked Questions

What is tumor mutational burden?

Tumor mutational burden (TMB) is the total number of somatic (non-inherited) mutations per megabase of tumor DNA. It is measured by whole-exome sequencing or targeted gene panels such as FoundationOne CDx. TMB-high tumors (≥10 mutations/Mb) generate more neoantigens and are more likely to respond to immune checkpoint inhibitors like pembrolizumab [2].

Is low tumor mutational burden good?

Low TMB is not inherently good or bad. Low-TMB tumors typically show less immune infiltration and are less likely to respond to immunotherapy. However, they often have lower proliferation rates and more favorable prognoses in certain cancer types like thyroid carcinoma. The clinical significance depends on the cancer type and therapeutic context.


References

  1. Wu HX et al. Beyond Tumor Mutation Burden: Tumor Neoantigen Burden as a Biomarker for Immunotherapy and Other Types of Therapy. Frontiers in Oncology, 2021.
  2. Marcus L et al. FDA Approval Summary: Pembrolizumab for the Treatment of Tumor Mutational Burden–High Solid Tumors. Clinical Cancer Research, 2021.
  3. Oh S et al. Spontaneous mutations in the single TTN gene represent high tumor mutation burden. npj Genomic Medicine, 2020.
  4. Hause RJ et al. Classification and characterization of microsatellite instability across 18 cancer types. Nature Medicine, 2016.
  5. Chalmers ZR et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Medicine, 2017.
  6. Jain MS, Massoud TF. Predicting tumour mutational burden from histopathological images using multiscale deep learning. Nature Machine Intelligence, 2020.
  7. Merino DM et al. Establishing guidelines to harmonize tumor mutational burden (TMB): in silico assessment of variation in TMB quantification across diagnostic platforms. Journal for ImmunoTherapy of Cancer, 2020.