Mean Depth 1000 More Than100×
Total Page:16
File Type:pdf, Size:1020Kb
A B ×20 ×50 ×100 100% 1000 800 90% 600 80% 400 Mean depth 70% 200 0 60% 50% Mean coverage 40% PBL002_T PBL003_T PBL004_T PBL005_T PBL006_T PBL007_T PBL008_T PBL009_T PBL010_T PBL001_P 30% 20% 10% 0% PBL002_T PBL003_T PBL004_T PBL005_T PBL006_T PBL007_T PBL008_T PBL009_T PBL010_T PBL001_P Figure S1. Depths and coverages of targeted deep sequencing. A. Depths of targeted deep sequencing on 10 PBL cases (mean depth = 677). B. Coverages of the targeted sequences analyzed at the indicated depths; 98.8% of the targeted regions were sequenced at a depth more than 100×. A B ×2 ×10 ×20 250 100% 200 90% 150 80% 100 Mean depth 50 70% 0 60% 50% PBL002_T PBL004_T PBL007_T PBL008_T PBL010_T PBL001_P PBL001_N PBL002_N PBL004_N PBL007_N PBL008_N PBL010_N PBL001_M1 PBL001_M2 PBL001_M3 PBL001_M4 PBL001_M5 PBL001_M6 Mean coverage 40% 30% 20% 10% 0% PBL002_T PBL004_T PBL007_T PBL008_T PBL010_T PBL001_P PBL001_N PBL002_N PBL004_N PBL007_N PBL008_N PBL010_N PBL001_M1 PBL001_M2 PBL001_M3 PBL001_M4 PBL001_M5 PBL001_M6 C D 14 80 12 70 10 60 50 8 40 6 30 4 Number of mutations Number of mutations 20 2 10 0 0 PBL002_T PBL004_T PBL007_T PBL008_T PBL010_T PBL001_P PBL001_P PBL001_M1 PBL001_M2 PBL001_M3 PBL001_M4 PBL001_M5 PBL001_M6 Nonsynonymous SNV Stopgain SNV Nonframeshift indels Frameshift indels Splice site Figure S2. Sequencing coverages and numbers of mutations detected by WES. A. Depths of WES on tumor and normal DNA samples from 6 cases (mean depth = 124 for tumor samples and 141 for normal samples). B. Coverages of the exome sequences; 97.5% of the exome sequences were analyzed at a depth more than 20×. C-D. The number of somatic non-silent mutations and indels detected by WES and validated by PCR-based amplicon deep sequencing in 6 primary tumor samples (C) and 7 multiple samples from PBL001 (D). SNV, single-nucleotide variant. P M1 M2 M3 CLEC19A NOTCH1 SIK3 MFSD2B CTNNB1 UBE2E2 GOLM1 PTPRU CTSW KIAA0368 UBE2E2 UBE2E2 PTPRU CCDC88C CTNNB1 KDM6A CTNNB1 MRGPRE BRS3 CCDC88C SYNJ1 OR4D11 KDM6A OR4D11 CCDC88C DSP KDM6A MRGPRE PLXNB1 CTNNB1 CCDC88C SIK3 ZC3H18 OR4D11 ABCA13 PGK2_c.A155T NTNG1_c.A300T HUWE1 OR4D11 LEF1 MAGEB1 PGK2_c.G154A KLKB1 NTNG1_c.C301T LRRIQ3 TMEM132C_c.A2550G UBE2E2 PHKA2 EXPH5 ADGRG2 MUC3A ADCY1 BRWD1 ZBED2 FOXP4 MRGPRE LCN9 TMPRSS5 SIK3 MYBPC1 IQCF1 GRM2 CYP2W1 TMEM132C_c.G2551T MRGPRE ABHD12B KDM6A SERPINB7 FRMD6 HECW1 PTPRU HOXB13 LOC650293 ZNF804B MYO6 XCR1 SIK3 GLRA3 MADD KMT2A CFAP43 NARFL DNAH17 EDN2 CUX1 PCLO HIST1H4D NUDCD3 PTPRU SLC5A10 NLRP10 PCNXL3 DLL1 AOC1 CNBP IFT43 GBF1 CEP295 GSX2 MKI67 CADPS2 ASCC3 KRTAP13-2 GALR1 CELSR1 DMBT1 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 Cellular Frequency Cellular Frequency Cellular Frequency Cellular Frequency M4 M5 M6 UBE2U SIK3 PKP3 KAT6B PTPRU MVK PIK3CD OGN CBS UBE2E2 OR51T1 KDM6A CCDC88C ADGB SIK3 KCNK17 ZBTB38 PPL MRGPRE BTBD16 MGMT OGN CTNNB1 OR51M1 MRGPRE CCDC88C MRGPRE GPR149 CTNNB1 OR51M1 FANCG OGN KDM6A OR51B4 SYNE1 KCTD1 UBE2E2 GLI3 PTPRU OR4D11 MLLT10 NPSR1 ARHGAP29 OR4D11 CUTC CDH3 PABPN1L UBE2E2 AHNAK OR4D11 DGCR14 TNXB OR51T1 ABCD4 SIK3 KDM6A GART PTPRU SETD4 ATR TIGD3 RAB11FIP3_c.C379A ZFHX3 OR51T1 CUL9 OR51M1 STIM1 CDKN2AIP GBP5 CTNNB1 XKRX ZNF516 RAB11FIP3_c.C380A TRNT1 KIAA2026 ADGRG4 RNF40 TRH HNRNPCL2 FA2H NTAN1 SLC17A6 MOXD1 SEMA3E CYSRT1 ZCCHC11 CD109 RPS6KA5 KCNQ5 COL6A5 RBM25 DMBT1 RANBP10 ALG1L TRH DLK1 LAMC3 KIAA1462 SLC34A2 TAF7L RNF38 HECW2 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 Cellular Frequency Cellular Frequency Cellular Frequency Figure S3. Tumor subclonal populations in multiple samples from PBL001 estimated by PyClone. Cellular frequencies of detected somatic mutaions in each sample from PBL001 are shown on each panel. Mutations on chromosomal regions with copy number gains were not used for the analysis because mutant cell fractions may be underestimated. P, primary sample; M1-M6, metastatic samples. B A Probability Probability 0.00 0.05 0.10 0.15 0.00 0.05 0.10 0.15 for primary(left),andmetastaticPBL(right). substitution patterns,the±2flankingbases,and transcriptionstranddirections,plotted immediately 3’and5’tothemutatedbase.B.Mutational signaturesconsistingofthe of substitutions.Eachsubstitutiontypecontains16 barsrepresentingthesequencecontext samples fromPBL001(lowerplot).Theprobability bars arecoloredaccordingtothe6types of 6primaryPBLsamplesinwhichWESwasperformed (upperplot),and6metastatic Figure S4. A C G Metastatic PBL Primary PBL ACA ACA ACC ACC ACG ACG ACT ACT CCA CCA T C>A CCC CCC C>A Primary PBL CCG CCG Mutationalsignaturesofprimaryandmetastatic PBL.A.Mutationalsignatures CCT CCT GCA GCA GCC GCC T C GCG GCG GCT GCT TCA TCA TCC TCC TCG TCG TCT TCT C G Strand directions ACA ACA ACC ACC ACG ACG ACT ACT A CCA CCA C>G C>G + CCC CCC G CCG CCG − CCT CCT GCA GCA GCC GCC GCG GCG GCT GCT TCA TCA TCC TCC TCG TCG TCT TCT ACA ACA C ACC ACC ACG ACG G ACT ACT T CCA CCA C>T CCC CCC C>T CCG CCG CCT CCT C GCA GCA GCC GCC GCG GCG Metastatic PBL T GCT GCT TCA TCA TCC TCC TCG TCG TCT TCT T C ATA ATA ATC ATC ATG ATG ATT ATT C CTA CTA T>A CTC CTC T>A Strand directions CTG CTG CTT CTT GTA GTA GTC GTC GTG GTG + GTT GTT G TTA TTA − TTC TTC T TTG TTG TTT TTT ATA ATA ATC ATC ATG ATG ATT ATT CTA CTA T>C CTC CTC T>C CTG CTG CTT CTT GTA GTA GTC GTC GTG GTG GTT GTT TTA TTA TTC TTC TTG TTG TTT TTT ATA ATA ATC ATC ATG ATG ATT ATT CTA CTA T>G CTC CTC T>G CTG CTG CTT CTT GTA GTA GTC GTC GTG GTG GTT GTT TTA TTA TTC TTC TTG TTG TTT TTT 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 2122 PBL001_P PBL001_M1 PBL001_M2 PBL001_M3 PBL001_M4 PBL001_M5 PBL001_M6 PBL002_T PBL004_T PBL007_T PBL008_T PBL010_T CN Loss Gain Figure S5. SNP array-based copy number analysis of PBL. Copy number alterations in multiple samples from PBL001 and 5 primary tumor samples identified by SNP array analysis are shown. CN, copy number. A PBL001 PBL002 Chr11 Chr11 Hetero Hetero SNPs SNPs 1 1 0.5 0.5 ratio ratio Allelic Allelic 0 0 PBL003 PBL004 Chr11 Chr11 Hetero Hetero SNPs SNPs 1 1 ratio 0.5 0.5 ratio Allelic Allelic 0 0 PBL005 PBL006 Chr11 Chr11 Hetero Hetero SNPs SNPs 1 1 ratio 0.5 ratio 0.5 Allelic Allelic 0 0 PBL007 PBL008 Chr11 Chr11 Hetero Hetero SNPs SNPs 1 1 0.5 0.5 ratio ratio Allelic Allelic 0 0 B PBL010 PBL009 Chr11 Chr11 Hetero Hetero SNPs SNPs 1 1 0.5 0.5 ratio ratio Allelic Allelic 0 0 Figure S6. Targeted deep sequencing based allele-specific copy number analysis of chromosome 11. A. Sequencing-based copy number analysis identified CN-LOH with the minimal common region identical to chromosome 11p15.5 locus in 8 cases. Allelic ratios were calculated by allele frequencies and sequenced depths of SNPs; red/green dots represent major/minor allele frequencies at SNP sites. B. Allelic ratio of PBL009, the case with gain of whole chromosome 11, is also plotted. A PBL002 (tumor purity = 0.57) PBL004 (tumor purity = 0.86) TTCGTTCGTGGAAACGTTTC GGGTTATTTAA GTTACGCGTCGT TTC GTTC GTGGAAACGT TTCGGGTTAT TTAA GTTACGCGTCGT (T) (T) (T) (T) (T) (T) (T) PBL010 (tumor purity = 0.44) TTCGTTCGTGGAAACGTTTCGGGTTATTTAA GTTACGCGTCGT (T) (T) (T) (T) (T) (T) (T) B PBL009 (tumor purity = 0.69) TTCGTTCGTGGAAACGTTTCGGGTTAT TTAA GTTACGCGTCGT (T) (T) (T) (T) (T) (T) (T) Figure S7. Bisulfite Sanger sequencing of ICR1 on chromosome 11p15.5. A. Bisulfite Sanger sequencing of 3 CN-LOH (+) cases (CpG sites are indicated by arrows). From the estimated tumor purity (see Methods), all samples examined were considered homozygously methylated, indicating that their CN-LOH causes paternal uniparental disomy. B. Bisulfite sequencing of PBL009, the case with gain of whole chromosome 11. A B * M 3.0 IGF2 ACTB ** * * 2.5 IGF2 1000 bp 500 bp 2.0 100 bp 1.5 xpression of e e v 1.0 Relati PBL001 PBL005PBL009 PBL001 PBL005PBL009 0.5 Normal pancreas (Negative control) Normal pancreas (Negative control) 0 Normal PBL001 PBL005 PBL009 pancreas Sample Figure S8. Relative expression levels of IGF2 in CN-LOH-negative samples. A. Reverse transcription PCR was performed in normal pancreas (as a normal control), PBL001 (as a CN-LOH-positive control), PBL005, and PBL009. ACTB (encoding β-actin) was used as a reference gene. The lanes of negative control contain water control. M, DNA size markers. B. DNA bands were quantified using Image Lab software (Bio-Rad). Relative intensities of IGF2 compared to ACTB were computed, and mean intensities ± s.d. of 3 technical replicates are plotted with corresponding P-values (Student’s t-test). *P < 0.05, **P < 0.01. A 3 2 ratio 1 Log2 0 CN = 2 -1 -2 1 0.5 ratio Allelic 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19202122 Chromosome B C 3 2 1 ratio 0 Log2 -1 -2 ● 1 ● ● ● ●● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ●● ● ● ● ● ● ● ● ● ● ● ●● ● ● ● ● ● ● ● ● 0.5 ● ●● ● ratio Allelic ● ● ● 0 Chr5 APC Figure S9. Aberrant activation of Wnt signaling pathway in PBL009. A. Sequencing-based copy number analysis identified CN-LOH within chromosome 5q involving APC gene locus. Upper panel shows the genome-wide copy number profile in PBL009. A log2 ratio of zero (horizontal red line) corresponds to a normal, diploid copy number.