AbstractExtreme responders to anticancer therapy are rare among advanced breast cancer patients. Researchers, however, have yet to investigate treatment responses therein on the whole exome level. We performed whole exome analysis to characterize the genomic landscape of extreme responders among metastatic breast cancer patients. Clinical samples were obtained from breast cancer patients who showed exceptional responses to anti-HER2 therapy or hormonal therapy and from those who did not. Matched breast tumor tissue (somatic DNA) and blood samples (germline DNA) were collected from a total of 30 responders and 15 non-responders. Whole exome sequencing using Illumina HiSeq2500 was performed for all 45 patients (90 samples). Somatic single nucleotide variants (SNVs), indels, and copy number variants (CNVs) were identified for the genomes of each patient. Group-specific somatic variants and mutational burden were statistically analyzed. Sequencing of cancer exomes for all patients revealed 1839 somatic SNVs (1661 missense, 120 nonsense, 43 splice-site, 15 start/stop-lost) and 368 insertions/deletions (273 frameshift, 95 in-frame), with a median of 0.7 mutations per megabase (range, 0.08 to 4.2 mutations per megabase). Responders harbored a significantly lower nonsynonymous mutational burden (median, 26 vs. 59, P = 0.02) and fewer CNVs (median 13.6 vs. 97.7, P = 0.05) than non-responders. Multivariate analyses of factors influencing progression-free survival showed that a high mutational burden and visceral metastases were significantly related with disease progression. Extreme responders to treatment for metastatic breast cancer are characterized by fewer nonsynonymous mutations and CNVs.