ctc-seq: accurate mutation detection from single …

1
SUMMARY • Circulating tumor cells (CTCs) harbor genetic signatures that may inform the clinical management of patients with targeted therapies. • Advances in single-cell isolation, genome amplification and low input mutation detection (1- 20 cells) are needed to enable genetic profiling of liquid biopsy samples. • Here we describe initial evaluation of a comprehensive solution for targeted next-generation sequencing (NGS) of CTCs, or CTC-Seq, using cancer cell lines spiked into blood as a model. • Accurate mutation detection can be achieved from as little as a single cell, with more robust results when using 5 or more cells. RESULTS CTC-SEQ: ACCURATE MUTATION DETECTION FROM SINGLE CANCER CELLS USING THE CELLSEARCH ® SYSTEM AND QUANTIDEX ® TARGETED NGS PANELS Stephanie Bridger 1 , Charles Saginario 2 , Christian Kazanecki 2 , Huiping Zhu 1 , Robert Zeigler 1 , Diane Ilsley 1 , and Gary J Latham 1 1 Asuragen Inc., Austin, TX; 2 Janssen Diagnostics, LLC, Raritan, NJ INTRODUCTION Liquid biopsies can forecast outcomes in cancer patients and report shifts in tumor progression in real- time. Enumeration of CTCs in peripheral blood has been clinically validated to predict overall survival of patients with metastatic breast, colorectal and prostate cancers (CELLSEARCH ® CTC Test). Yet CTCs also harbor genetic signatures that may inform the clinical management of patients with targeted therapies. NGS is a compelling option to profile “druggable” mutations in CTCs, but advances are needed to accommodate low CTC counts. Here we describe initial evaluation of a comprehensive solution for targeted NGS of CTCs, or CTC-Seq. This approach combines sample qualification, targeted enrichment of clinically-actionable cancer genes, and a customized bioinformatics pipeline tuned for high positive predictive value (PPV). METHODS Four cancer cell lines were spiked into 7.5 mL of normal human blood, enriched using the CELLSEARCH ® System (Janssen Diagnostics, JDx) and further purified by individual cell picking. Control cells were directly purified from culture by individual cell picking. The cell lines represent 6 “driver” mutations in 4 oncogenes (BRAF, EGFR, KRAS and PIK3CA). Phi29-based whole genome amplification (WGA) (REPLI-g SC, Qiagen) was used to amplify genetic material from low input cell samples (1, 3, 5, 10 or 20 picked cells) and high input control gDNA (10 ng). Amplifiable DNA inputs into NGS were qualified with the QuantideX ® qPCR DNA QC Assay (Asuragen, Asu). PCR-based target enrichment, library preparation, and variant analysis were performed using the QuantideX ® NGS Pan Cancer Kit* and QuantideX ® NGS Reporter software (Asu) and sequenced on the MiSeq ® System (Illumina). CONCLUSIONS • Our results demonstrate proof-of-concept for CTC-Seq mutation detection from 1-20 cancer cells captured by the CELLSEARCH ® System and further purified by single cell picking. • QuantideX qPCR DNA QC Assay was predictive of WGA sample performance on NGS. • Adequate sequence coverage was reached for >90% of the QuantideX NGS Pan Cancer target genes enabling accurate mutation detection from as little as a single cell, with more robust results when using 5 or more cells. • EDTA sample collection tubes support more robust NGS analysis and mutation detection than CellSave tubes. • Further optimization of the CTC-Seq workflow and its application to clinical CTC specimens are in progress. Figure 1. CTC-Seq proof-of-concept study design and NGS panel description. (A) Establish feasibility for detecting known variants from limited cell inputs (1 – 20 cells) into WGA as compared to purified genomic DNA templates. Determine cell input requirements for accuracy of mutation detection. Evaluate sample collection tubes to assess compatibility with CELLSEARCH ® and QuantideX workflows. Perform stability study to determine maximum limits for sample storage prior to analysis. (B) The QuantideX NGS Pan Cancer Panel consists of 46 amplicons and covers >1,600 known COSMIC variants in 21 genes of clinical significance in human cancers including: colorectal, breast, melanoma and NSCLC. Expected cell line mutations are represented on 6 of the QuantideX NGS Pan Cancer amplicons (circled above). Figure 2. Integrated CTC-Seq workflow. Cancer cell lines were enriched by the CELLSEARCH ® System and further purified by cell picking (JDx). WGA from as little as a single cell was conducted by two independent sites (Asu and JDx). QuantideX NGS Pan Cancer Kit and NGS Reporter software was used for sample QC, NGS library prep and data reporting (Asu). Sequencing was performed on the MiSeq System (Illumina). Figure 5. Mutation detection and time course stability in CellSave. (A) Accuracy (PPV) and sensitivity of mutation detection is summarized by site where WGA was conducted (JDx and Asu). Cellular inputs ranging from 1 to 20 cells (CellSave < 24 hrs) were tested in duplicate (triplicate for the single cell inputs). Expected mutations were identified from as few as a single cell in all 4 cell lines, with the most accurate results (86 - 100% PPV and analytical sensitivity) achieved with 5 or more cells. (B) Sample stability was evaluated using 3 and 5 cell inputs into CellSave blood for 24-96 hrs prior to enrichment and cell picking. Amplifiable copy number of the resulting WGA products were predictive of successful variant detection. PPV decreased with length of time in CellSave, while ~80% sensitivity was maintained across all points measured. Figure 6. Sample collection comparison (CellSave vs. EDTA tubes) and implications on NGS performance. NCl-H1975 and MCF-7 cell lines at 3 and 5 cell inputs were spiked into CellSave blood (in duplicate) and EDTA blood (in triplicate) for up to 24 hrs. EDTA samples consistently outperformed CellSave samples in terms of usable samples (96% EDTA, 75% CellSave) (A), mutation detection (Sensitivity/PPV for NCl-H1975: 100%/92% EDTA, 91%/91% CellSave) (B) and coverage uniformity (up to 20% increase for EDTA) (C). Time course stability in EDTA suggests number of usable samples and mutation detection is comparable out to 96 hrs (data not shown). Successful attempts were made to recover samples in the “risk” DNA QC category (10 – 100 copies/uL); however samples with < 10 copies/uL consistently failed to yield usable sequence data. Figure 4. Coverage uniformity is a function of sample quality and input. The QuantideX NGS Pan Cancer Panel is optimized for low-quality FFPE DNA with > 90% coverage uniformity across all 46 amplicons. (A) Intact control gDNA and WGA-treated gDNA samples display similar uniformity results suggesting minimal WGA-induced bias; however, (B) combinatorial effects associated with low input WGA (1 - 20 cells) and DNA fixation in CellSave produced characteristic MDA-induced amplification bias resulting in allelic dropout in AKT1, FGFR3, HRAS, and RET regions (CellSave Cells). Of note, FGFR3 and HRAS loci are each GC-rich (60-73%) which may contribute to reduced amplification uniformity in the WGA step. (C) Coverage uniformity in CellSave improved as a function of cell number and amplifiable DNA inputs, where 80% of 20 cell inputs showed comparable uniformity to control cell samples. Figure 3. QuantideX qPCR DNA QC Assay predicts WGA sample performance on NGS. (A) WGA products from control cells and cells spiked into CellSave blood for < 24 hr prior to enrichment were QC’d. Samples with < 100 amplifiable copies/uL were flagged at risk for poor performance on NGS. 98% of the WGA samples from control cells (excluding a single 1 cell input) yielded sufficient amplifiable DNA. A 17% decrease in number of “usable” CellSave samples compared to controls (JDx site) was observed, especially below 5 cell input. Asu CellSave samples showed lower pass rate (46%) compared to JDx (78%), consistent with a batch handling effect. (B) Median amplicon coverage > 500X indicates successful performance on NGS; however high-sensitivity variant calling is possible down to ~200X. CELLSEARCH ® SYSTEM & CELL PICKING Blood Collection Rare Cell Isolation & Staining Cell Picking & Pooling CellSave Preservative Tubes or EDTA Tubes CELLSEARCH ® CXC Kit CELLTRACKS ® AUTOPREP ® System Individual Cell Picking Whole Genome Amplification (WGA) Multiple Displacement Amplification (MDA) Analytics & Reporting QuantideX Reporter Push-Button Analytics & Reporting Suite QUANTIDEX NGS PAN CANCER KIT COMPONENTS Sample QC & Quantification Single-Tube Gene-Specific PCR Library Prep Sample Indexing & Sequencing Primers Library Purification & Quantification QuantideX DNA Assay QuantideX NGS Pan Cancer Panel QuantideX Codes QuantideX Library Pure Prep & Quant Quant Readout Inhibition Readout Gene-Specific PCR FFPE & Multi- Variant Control Barcode Tagging (Dual index Illumina or PGM Chemistries) Sequencing Primers Included for Post-Library Purification Bead Library Purification SEQUENCING* *not included QUANTIDEX NGS PAN CANCER PANEL Gene Codon Range Gene Codon Range ABL1 249-258 HRAS 9-20 303-319 59-76 AKT1 16-27 113-121 AKT2 16-26 IDH1 122-134 ALK1 1174-1196 IDH2 122-134 1274-1278 163-174 BRAF 465-474 JAK2 607-620 591-612 KIT 557-579 EGFR 486-493 815-826 709-722 KRAS 9-20 737-761 55-65 767-798 104-118 849-861 137-148 ERBB2 755-769 MET 1245-1256 774-788 NRAS 9-20 839-847 55-67 877-883 110-119 FGFR1 123-136 144-150 250-262 PDGFRA 560-572 FGFR3 247-260 840-852 363-374 PIK3CA 540-551 638-653 1038-1049 FLT3 829-840 RET 916-926 Accuracy of Mutation Detection Cell Input Requirements Sample Collection Sample Stability Cancer Cell Lines Control gDNA & WGA-gDNA Cell # Input Spiked Blood Control Cells Incubation Time (hr) NCI-H1975 MCF-7 SK-Mel-28 MDA-MB-231 1, 3, 5, 10, 20 CellSave EDTA From Culture 0, 24, 48, 96 CTC-SEQ PROOF-OF-CONCEPT STUDY 1A 3A 5A 6A 4A 1B 3B 5B 6B Control Cells CellSave Cells Sample Description # Cell Asu JDx Asu JDx H1975 1 0 21474 1 13212 0 998 26586 1 3467 6300 0 0 3 12718 7604 0 18405 5 8410 4259 0 11481 10 6408 5139 4476 4271 20 2786 5086 0 11446 MCF-7 1 0 1285 1 43719 2227 0 32 1 5702 5449 8 0 3 9091 4371 5056 33 5 4798 4904 0 1208 10 6652 7167 1 2507 20 10862 4829 11619 4011 SK-Mel-28 1 0 1 1 8197 3238 0 541 1 10608 2480 1231 550 3 20662 3147 22853 618 5 3232 4769 0 4328 10 8118 5285 25295 4317 20 12206 6341 10799 2964 MDA-MB-231 1 5 737 1 7132 20728 237 2660 1 14956 15324 0 0 3 13729 9032 2394 9071 5 13920 8596 5788 17324 10 10342 10729 5074 6813 20 10509 8857 17015 12762 WGA NTC control 0 0 0 # Samples At Risk 0 1 15 6 6C 4B 4C *Research Use Only – Not For Use In Diagnostic Procedures Preliminary research data. The performance characteristics of this assay have not yet been established. Presented at AMP 2015

Upload: others

Post on 28-Jul-2022

4 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: CTC-SEQ: ACCURATE MUTATION DETECTION FROM SINGLE …

SUMMARY• Circulating tumor cells (CTCs) harbor genetic signatures that may inform the clinical management

of patients with targeted therapies. • Advances in single-cell isolation, genome amplification and low input mutation detection (1- 20

cells) are needed to enable genetic profiling of liquid biopsy samples.• Here we describe initial evaluation of a comprehensive solution for targeted next-generation

sequencing (NGS) of CTCs, or CTC-Seq, using cancer cell lines spiked into blood as a model.• Accurate mutation detection can be achieved from as little as a single cell, with more robust

results when using 5 or more cells.

RESULTS

CTC-SEQ: ACCURATE MUTATION DETECTION FROM SINGLE CANCER CELLS USING THE CELLSEARCH® SYSTEM AND QUANTIDEX® TARGETED NGS PANELSStephanie Bridger1, Charles Saginario2, Christian Kazanecki2, Huiping Zhu1, Robert Zeigler1, Diane Ilsley1, and Gary J Latham1

1Asuragen Inc., Austin, TX; 2Janssen Diagnostics, LLC, Raritan, NJ

INTRODUCTIONLiquid biopsies can forecast outcomes in cancer patients and report shifts in tumor progression in real-time. Enumeration of CTCs in peripheral blood has been clinically validated to predict overall survival of patients with metastatic breast, colorectal and prostate cancers (CELLSEARCH® CTC Test). Yet CTCs also harbor genetic signatures that may inform the clinical management of patients with targeted therapies. NGS is a compelling option to profile “druggable” mutations in CTCs, but advances are needed to accommodate low CTC counts. Here we describe initial evaluation of a comprehensive solution for targeted NGS of CTCs, or CTC-Seq. This approach combines sample qualification, targeted enrichment of clinically-actionable cancer genes, and a customized bioinformatics pipeline tuned for high positive predictive value (PPV).

METHODSFour cancer cell lines were spiked into 7.5 mL of normal human blood, enriched using the CELLSEARCH® System (Janssen Diagnostics, JDx) and further purified by individual cell picking. Control cells were directly purified from culture by individual cell picking. The cell lines represent 6 “driver” mutations in 4 oncogenes (BRAF, EGFR, KRAS and PIK3CA). Phi29-based whole genome amplification (WGA) (REPLI-g SC, Qiagen) was used to amplify genetic material from low input cell samples (1, 3, 5, 10 or 20 picked cells) and high input control gDNA (10 ng). Amplifiable DNA inputs into NGS were qualified with the QuantideX® qPCR DNA QC Assay (Asuragen, Asu). PCR-based target enrichment, library preparation, and variant analysis were performed using the QuantideX® NGS Pan Cancer Kit* and QuantideX® NGS Reporter software (Asu) and sequenced on the MiSeq® System (Illumina).

CONCLUSIONS• Our results demonstrate proof-of-concept for CTC-Seq mutation detection from 1-20 cancer

cells captured by the CELLSEARCH® System and further purified by single cell picking. • QuantideX qPCR DNA QC Assay was predictive of WGA sample performance on NGS.• Adequate sequence coverage was reached for >90% of the QuantideX NGS Pan Cancer target

genes enabling accurate mutation detection from as little as a single cell, with more robust results when using 5 or more cells.

• EDTA sample collection tubes support more robust NGS analysis and mutation detection than CellSave tubes.

• Further optimization of the CTC-Seq workflow and its application to clinical CTC specimens are in progress.

Figure 1. CTC-Seq proof-of-concept study design and NGS panel description. (A) Establish feasibility for detecting known variants from limited cell inputs (1 – 20 cells) into WGA as compared to purified genomic DNA templates. Determine cell input requirements for accuracy of mutation detection. Evaluate sample collection tubes to assess compatibility with CELLSEARCH® and QuantideX workflows. Perform stability study to determine maximum limits for sample storage prior to analysis. (B) The QuantideX NGS Pan Cancer Panel consists of 46 amplicons and covers >1,600 known COSMIC variants in 21 genes of clinical significance in human cancers including: colorectal, breast, melanoma and NSCLC. Expected cell line mutations are represented on 6 of the QuantideX NGS Pan Cancer amplicons (circled above).

Figure 2. Integrated CTC-Seq workflow. Cancer cell lines were enriched by the CELLSEARCH® System and further purified by cell picking (JDx). WGA from as little as a single cell was conducted by two independent sites (Asu and JDx). QuantideX NGS Pan Cancer Kit and NGS Reporter software was used for sample QC, NGS library prep and data reporting (Asu). Sequencing was performed on the MiSeq System (Illumina).

Figure 5. Mutation detection and time course stability in CellSave. (A) Accuracy (PPV) and sensitivity of mutation detection is summarized by site where WGA was conducted (JDx and Asu). Cellular inputs ranging from 1 to 20 cells (CellSave < 24 hrs) were tested in duplicate (triplicate for the single cell inputs). Expected mutations were identified from as few as a single cell in all 4 cell lines, with the most accurate results (86 - 100% PPV and analytical sensitivity) achieved with 5 or more cells. (B) Sample stability was evaluated using 3 and 5 cell inputs into CellSave blood for 24-96 hrs prior to enrichment and cell picking. Amplifiable copy number of the resulting WGA products were predictive of successful variant detection. PPV decreased with length of time in CellSave, while ~80% sensitivity was maintained across all points measured.

Figure 6. Sample collection comparison (CellSave vs. EDTA tubes) and implications on NGS performance. NCl-H1975 and MCF-7 cell lines at 3 and 5 cell inputs were spiked into CellSave blood (in duplicate) and EDTA blood (in triplicate) for up to 24 hrs. EDTA samples consistently outperformed CellSave samples in terms of usable samples (96% EDTA, 75% CellSave) (A), mutation detection (Sensitivity/PPV for NCl-H1975: 100%/92% EDTA, 91%/91% CellSave) (B) and coverage uniformity (up to 20% increase for EDTA) (C). Time course stability in EDTA suggests number of usable samples and mutation detection is comparable out to 96 hrs (data not shown). Successful attempts were made to recover samples in the “risk” DNA QC category (10 – 100 copies/uL); however samples with < 10 copies/uL consistently failed to yield usable sequence data.

Figure 4. Coverage uniformity is a function of sample quality and input. The QuantideX NGS Pan Cancer Panel is optimized for low-quality FFPE DNA with > 90% coverage uniformity across all 46 amplicons. (A) Intact control gDNA and WGA-treated gDNA samples display similar uniformity results suggesting minimal WGA-induced bias; however, (B) combinatorial effects associated with low input WGA (1 - 20 cells) and DNA fixation in CellSave produced characteristic MDA-induced amplification bias resulting in allelic dropout in AKT1, FGFR3, HRAS, and RET regions (CellSave Cells). Of note, FGFR3 and HRAS loci are each GC-rich (60-73%) which may contribute to reduced amplification uniformity in the WGA step. (C) Coverage uniformity in CellSave improved as a function of cell number and amplifiable DNA inputs, where 80% of 20 cell inputs showed comparable uniformity to control cell samples.

Figure 3. QuantideX qPCR DNA QC Assay predicts WGA sample performance on NGS. (A) WGA products from control cells and cells spiked into CellSave blood for < 24 hr prior to enrichment were QC’d. Samples with < 100 amplifiable copies/uL were flagged at risk for poor performance on NGS. 98% of the WGA samples from control cells (excluding a single 1 cell input) yielded sufficient amplifiable DNA. A 17% decrease in number of “usable” CellSave samples compared to controls (JDx site) was observed, especially below 5 cell input. Asu CellSave samples showed lower pass rate (46%) compared to JDx (78%), consistent with a batch handling effect. (B) Median amplicon coverage > 500X indicates successful performance on NGS; however high-sensitivity variant calling is possible down to ~200X.

CELLSEARCH® SYSTEM & CELL PICKING

BloodCollection

Rare Cell Isolation

& StainingCell Picking& Pooling

CellSavePreservative

Tubes orEDTA Tubes

CELLSEARCH®

CXC Kit

CELLTRACKS®

AUTOPREP®

System

IndividualCell Picking

Whole GenomeAmplification (WGA)

Multiple DisplacementAmplification

(MDA)

Analytics &ReportingQuantideXReporter

Push-ButtonAnalytics &

Reporting Suite

QUANTIDEX NGS PAN CANCER KIT COMPONENTS

Sample QC &Quantification

Single-Tube Gene-Specific

PCR Library Prep

Sample Indexing & Sequencing

Primers

Library Purification & Quantification

QuantideXDNA Assay

QuantideXNGS Pan

Cancer Panel

QuantideXCodes

QuantideXLibrary Pure

Prep & Quant

Quant Readout

Inhibition Readout

Gene-Specific PCR

FFPE & Multi-Variant Control

Barcode Tagging (Dual index Illumina

or PGM Chemistries)

Sequencing Primers Included for Post-Library

Purification

Bead Library Purification

SE

QU

EN

CIN

G*

*not included

QUANTIDEX NGS PAN CANCER PANEL

Gene Codon Range Gene Codon Range

ABL1249-258

HRAS

9-20

303-319 59-76

AKT1 16-27 113-121

AKT2 16-26 IDH1 122-134

ALK11174-1196

IDH2122-134

1274-1278 163-174

BRAF465-474 JAK2 607-620

591-612KIT

557-579

EGFR

486-493 815-826

709-722

KRAS

9-20

737-761 55-65

767-798 104-118

849-861 137-148

ERBB2

755-769 MET 1245-1256

774-788

NRAS

9-20

839-847 55-67

877-883 110-119

FGFR1123-136 144-150

250-262PDGFRA

560-572

FGFR3

247-260 840-852

363-374PIK3CA

540-551

638-653 1038-1049

FLT3 829-840 RET 916-926

Accuracyof MutationDetection

Cell InputRequirements

SampleCollection

SampleStability

Cancer Cell LinesControl gDNA &

WGA-gDNA

Cell # Input

Spiked BloodControl Cells

IncubationTime (hr)

NCI-H1975MCF-7

SK-Mel-28MDA-MB-231

1, 3, 5, 10, 20

CellSaveEDTA

From Culture

0, 24, 48, 96CTC

-SE

Q P

RO

OF-

OF-

CO

NC

EP

T S

TUD

Y

1A

3A

5A

6A

4A

1B

3B

5B

6B

Control Cells CellSave Cells

Sample Description # Cell Asu JDx Asu JDx

H1975

1 0 214741 13212 0 998 265861 3467 6300 0 03 12718 7604 0 184055 8410 4259 0 1148110 6408 5139 4476 427120 2786 5086 0 11446

MCF-7

1 0 12851 43719 2227 0 321 5702 5449 8 03 9091 4371 5056 335 4798 4904 0 120810 6652 7167 1 250720 10862 4829 11619 4011

SK-Mel-28

1 0 11 8197 3238 0 5411 10608 2480 1231 5503 20662 3147 22853 6185 3232 4769 0 432810 8118 5285 25295 431720 12206 6341 10799 2964

MDA-MB-231

1 5 7371 7132 20728 237 26601 14956 15324 0 03 13729 9032 2394 90715 13920 8596 5788 1732410 10342 10729 5074 681320 10509 8857 17015 12762

WGA NTC control 0 0 0# Samples At Risk 0 1 15 6

6C

4B

4C

* Research Use Only – Not For Use In Diagnostic ProceduresPreliminary research data. The performance characteristics of this assay have not yet been established.Presented at AMP 2015