Gallbladder adenocarcinomas undergo subclonal diversification and selection from precancerous lesions to metastatic tumors

  1. Minsu Kang
  2. Hee Young Na
  3. Soomin Ahn  Is a corresponding author
  4. Ji-Won Kim  Is a corresponding author
  5. Sejoon Lee
  6. Soyeon Ahn
  7. Ju Hyun Lee
  8. Jeonghwan Youk
  9. Haesook T Kim
  10. Kui-Jin Kim
  11. Koung Jin Suh
  12. Jun Suh Lee
  13. Se Hyun Kim
  14. Jin Won Kim
  15. Yu Jung Kim
  16. Keun-Wook Lee
  17. Yoo-Seok Yoon
  18. Jee Hyun Kim
  19. Jin-Haeng Chung
  20. Ho-Seong Han
  21. Jong Seok Lee
  1. Department of Internal Medicine, Seoul National University Bundang Hospital, Republic of Korea
  2. Department of Pathology, Seoul National University Bundang Hospital, Republic of Korea
  3. Department of Pathology and Translational Genomics, Samsung Medical Center, Sungkyunkwan University School of Medicine, Republic of Korea
  4. Genealogy Inc, Republic of Korea
  5. Center for Precision Medicine, Seoul National University Bundang Hospital, Republic of Korea
  6. Medical Research Collaboration Center, Seoul National University Bundang Hospital, Republic of Korea
  7. Department of Data Science, Dana Farber Cancer Institute, Harvard T.H. Chan School of Public Health, United States
  8. Biomedical Research Institute, Seoul National University Bundang Hospital, Republic of Korea
  9. Department of Surgery, Seoul National University Bundang Hospital, Republic of Korea
5 figures, 3 tables and 3 additional files

Figures

Figure 1 with 2 supplements
Clinical history of patients, mutational landscape, and study workflow.

(A) Clinical history of 11 patients is summarized in the swimmer plot. (B) The mutational landscape of 11 primary tumors (GB) is visualized and compared with the two previous studies on gallbladder …

Figure 1—figure supplement 1
Inference of clonal structure by PyClone algorithm.

(A–K) The clonal population structure was statistically inferred for GB-A1 (A), GB-S1 (B), GB-S2 (C), GB-S3 (D), GB-A2 (E), GB-S4 (F), GB-S5 (G), GB-S6 (H), GB-S7 (I), GB-S8 (J), and GB-S9 (K) using …

Figure 1—figure supplement 2
BiIlN and primary gallbladder adenocarcinoma (GBAC) of the GB-S7 patient presumed to be derived from different origins.

(A) H&E staining of BiIlN. (B) H&E staining of GBAC. (C) Variant allele frequency of mutated genes in BilIN and primary GBAC. BilIN, biliary intraepithelial neoplasia; H&E, hematoxylin and eosin.

Spatial and temporal clonal evolution of four patients with gallbladder adenocarcinoma (GBAC) whose precancerous BilIN tissues were analyzed.

(A–D) The most probable phylogenetic trees and MapScape and TimeScape results are visualized for GB-A1 (A), GB-S1 (B), GB-S2 (C), and GB-S3 (D). In MapScape visualization (table), the numbers in the …

Spatial and temporal clonal evolution of additional seven patients with gallbladder adenocarcinoma (GBAC).

(A–G) The most probable phylogenetic trees and MapScape and TimeScape results are visualized for GB-A2 (A), GB-S4 (B), GB-S5 (C), GB-S6 (D), GB-S7 (E), GB-S8 (F), and GB-S9 (G). In MapScape …

Figure 4 with 3 supplements
Mutational signature analysis.

(A–B) The 100% stacked bar plots compare the proportions of known COSMIC Mutational Signatures v2 within our dataset and two public (MSKCC and Shanhai) datasets (A), and each category split …

Figure 4—figure supplement 1
Mutational signatures within each sample from 11 gallbladder adenocarcinoma (GBAC) patients.

The 100% stacked bar plots comparing the proportions of known COSMIC Mutational Signatures v2 within each sample from 11 GBAC patients. BilIN, biliary intraepithelial neoplasia; CBD, common bile …

Figure 4—figure supplement 2
Mutational signatures within each category split according to the type of sample.

The 100% stacked bar plots comparing the proportions of known COSMIC Mutational Signatures v2 within each category split according to the type of sample. The total number of mutations (N) and cosine …

Figure 4—figure supplement 3
Mutational signature analysis validated by two additional tools, Signal and MuSiCa.

The 100% stacked bar plots comparing the proportions of known COSMIC Mutational Signatures v2 within each category split according to the timing of development during clonal evolution. The total …

ERBB2 copy number variation during neoplastic transformation of BilIN in GB-A1.

(A–G) ERBB2 gene amplification (A and B), HER2 SISH (C and D), and HER2 IHC (F and G) were compared between BilIN (A, C) and (F) and primary GBAC (B), (D), and (G) samples and the mean ERBB2/CEP17 …

Tables

Table 1
Baseline characteristics of 11 patients with GBAC.
Patient ID*SexAge at diagnosisECOG PS at diagnosisStage at diagnosisDifferentiation
GB-A1F701IVPD
GB-S1F661IVMD
GB-S2F661IVMD
GB-S3F751IIIMD
GB-A2M611IVMD
GB-S4M721IVMD
GB-S5F701IVPD
GB-S6M700IVMD
GB-S7F741IVPD
GB-S8M671IIIWD
GB-S9M590IVMD
  1. *

    In patient ID, ‘A’ indicates autopsy cases whereas ‘S’ indicates surgery cases.

  2. GBAC, gallbladder adenocarcinoma; ECOG PS, Eastern Cooperative Oncology Group performance status; M, male; F, female; PD, poorly differentiated; MD, moderately differentiated; WD, well-differentiated.

Table 2
List of subclones expanding during metastasis.
Patient IDSubcloneNo. of mutationsPutative driver mutationsClonal prevalence in primary tumorClonal prevalence in metastasisMetastatic organ
GB-A1E121JARID2 p.K603Q (LOH)0.2%20.5%CBD
F54SMAD4 p.L414fs (LOH)0.4%49.8%Omentum 1
27.0%Omentum 2
GB-S1E6ROBO1 p.P1360Q0.0%59.2%Distant LN
GB-S2H61.2%28.1%Regional LN
GB-A2F14PRKCD p.I153L0.3%73.1%Liver
68.3%Mesentery
G3DICER1 p.T519A2.8%39.9%Lung
61.5%Chest wall
GB-S4F32FBXW2 p.W450C0.0%34.1%Colon wall
H360.0%34.9%Distant LN
GB-S5B28KIAA0100 p.F5S1.2%54.4%Liver
CSMD2 p.E411K
GB-S6B7OSCP1 p.R351X1.1%60.2%Lung
GB-S7C53.3%24.7%Regional LN
  1. Putative driver mutations are indicated, and a full list of mutated genes is specified in Supplementary file 2.

  2. GB, gallbladder; LN, lymph node; CBD, common bile duct.

Author response table 1
Proportion of mutations satisfying the conditions: (1) C:G > T:A variants, and (2) VAF < 5% in sample-specific clusters.
Patient IDSampleType of sampleSample specific clusters*Number of mutations in sample-specific clustersNumber of filtered mutations in sample-specific clustersPercentage (%)
GB-A1BilINFFPE23438.8
CBDFFPE31211310.7
GB-S1BilINFFPE29333.3
GBFFPE418738.9
Distant LNFFPE7600.0
GB-S2BilINFFPE68337.5
GBFFPE07228.6
GB-S3BilINFFPE37114.3
GBFFPE2300.0
GB-A2GBFresh-frozen24000.0
Lung (old)FFPE101715.9
LungFresh-frozen51500.0
Chest wallFresh-frozen32300.0
LiverFresh-frozen43100.0
MesenteryFresh-frozen68225.0
GB-S4Colon wallFFPE33600.0
Distant LNFFPE19222.2
GB-S5GBFFPE31300.0
GB-S6GBFFPE1500.0
LungFFPE2700.0
GB-S7GBFFPE01100.0
Regional LNFFPE3500.0
GB-S8GBFFPE235411.4
Regional LNFFPE08133.7

Additional files

Supplementary file 1

Tables.

(A) Baseline characteristics of samples are summarized. (B) Mutational signatures of our dataset are analyzed by three different tools, Mutalisk, Signal, and MuSiCa.

https://cdn.elifesciences.org/articles/78636/elife-78636-supp1-v3.docx
Supplementary file 2

Detected somatic alternations in 11 patients with gallbladder adenocarcinoma (GBAC).

Full list of mutations called in our cohort.

https://cdn.elifesciences.org/articles/78636/elife-78636-supp2-v3.xlsx
MDAR checklist
https://cdn.elifesciences.org/articles/78636/elife-78636-mdarchecklist1-v3.docx

Download links