Detection of low-frequent mitochondrial DNA variants using SMRT sequencing

Transcription

1 Detection of low-frequent mitochondrial DNA variants using SMRT sequencing Marjolein J.A. Weerts SMRT Leiden 2018 June 13

2 Content Mitochondrial DNA & liquid biopsy in oncology Pitfalls when studying human mitochondrial DNA Approach to detect low-frequent mitochondrial DNA variants using SMRT sequencing Application in a pilot study 2

3 Human mitochondrial DNA Inheritance Structure Nuclear DNA Maternal & Paternal Double stranded lineair Mitochondrial DNA Maternal Double stranded circular Size 3.1 Gb 16.6 kb Cellular content Copies/cell Two* Hundreds to Thousands Gene density ~ 1 / 120 kb ~ 1 / 0.45 kb Mutation rate ~10-fold higher * With the exception of aneuploid tumor cells and germ cells 1

4 Liquid biopsy in oncology Less invasive than a regular biopsy! Blood-circulating cell-free DNA (cfdna) Released by normal as well as tumor cells Genetic information on the patient s cancer (primary tumor ánd metastases) 3 JCO 32(6): , 2014

5 Studying human mtdna in cancer HYPOTHESIS The high mtdna copy number per cell as well as the high mtdna mutation rate make it worthwhile to explore the potential of tumorspecific cf-mtdna variants as cancer marker in the blood of cancer patients 4

6 Studying human mtdna in cancer HYPOTHESIS The high mtdna copy number per cell as well as the high mtdna mutation rate make it worthwhile to explore the potential of tumorspecific cf-mtdna variants as cancer marker in the blood of cancer patients AIMS Detection of low-frequent mtdna variants Highly mtdna-specific detection of mtdna variants 4

7 Pitfalls when studying human mtdna High mutation rate Heteroplasmy Nuclear insertions of mitochondrial origin: NUMTs 5

8 Pitfalls when studying human mtdna High mutation rate Several orders of magnitude higher than that of ndna Nearly 10,000 variable positions reported in databases (e.g. mitomap, mtdb) 6

9 Pitfalls when studying human mtdna Heteroplasmy Genetically different mtdna molecules within a single cell Heteroplasmy patterns differ between tissues within an individual 7

10 Pitfalls when studying human mtdna Nuclear insertions of mitochondrial origin: NUMTs Each 175 bp mtdna segment: ~9.5 NUMT copies Cancer cells contain somatic NUMT insertion events NUMT similarity to mtdna interferes with variant detection 9

11 Studying human mtdna in cancer HYPOTHESIS The high mtdna copy number per cell as well as the high mtdna mutation rate make it worthwhile to explore the potential of tumorspecific cf-mtdna variants as cancer marker in the blood of cancer patients AIMS Detection of low-frequent mtdna variants Highly mtdna-specific detection of mtdna variants 10

12 Highly mtdna-specific detection of mtdna variants Minimize the interference of NUMTs: pure mtdna sample (devoid of nuclear DNA) sequence large fragments SMRT sequencing (~80 of human NUMTs are < 500 bp mtdna segments) 11 MJA Weerts et al., Scientific Reports 2018; 8:2261

13 SMRT sequencing approach for mtdna Amplification of mtdna (15-18 PCR cycles) 9x bp amplicons covering the entire mtdna Purify & pool amplicons (equimolar) per sample Barcode each sample (5 PCR cycles), purify & pool samples Generate SMRTbell library 5 Sequencing on RSII or Sequel 12 MJA Weerts et al., Scientific Reports 2018; 8:2261

14 SMRT sequencing approach for mtdna Amplification of mtdna (15-18 PCR cycles) 9x bp amplicons covering the entire mtdna Purify & pool amplicons (equimolar) per sample Barcode each sample (5 PCR cycles), purify & pool samples Generate SMRTbell library 5 Sequencing on RSII or Sequel 12 MJA Weerts et al., Scientific Reports 2018; 8:2261

15 SMRT sequencing approach for mtdna Amplification of mtdna (15-18 PCR cycles) 9x bp amplicons covering the entire mtdna Purify & pool amplicons (equimolar) per sample Barcode each sample (5 PCR cycles), purify & pool samples Generate SMRTbell library 5 Sequencing on RSII or Sequel 12 MJA Weerts et al., Scientific Reports 2018; 8:2261

16 SMRT sequencing approach for mtdna Amplification of mtdna (15-18 PCR cycles) 9x bp amplicons covering the entire mtdna Purify & pool amplicons (equimolar) per sample Barcode each sample (5 PCR cycles), purify & pool samples Generate SMRTbell library 5 Sequencing on RSII or Sequel 12 MJA Weerts et al., Scientific Reports 2018; 8:2261

17 SMRT sequencing approach for mtdna Amplification of mtdna (15-18 PCR cycles) 9x bp amplicons covering the entire mtdna Purify & pool amplicons (equimolar) per sample Barcode each sample (5 PCR cycles), purify & pool samples Generate SMRTbell library 5 Sequencing on RSII or Sequel 12 MJA Weerts et al., Scientific Reports 2018; 8:2261

18 SMRT sequencing approach for mtdna Generate circular consensus reads (CCS2 algorithm) Attribute reads using sample-specific barcode (TSSV) Select highly accurate CCS reads (QV>=99, >= 5 passes) Trim reads (Cutadapt) and align to extended version or rcrs (BWA-MEM) Call positions alternative to the reference sequence in pileup (Rsamtools, min_nucleotide_depth=5) 13 MJA Weerts et al., Scientific Reports 2018; 8:2261

19 SMRT sequencing approach for mtdna Generate circular consensus reads (CCS2 algorithm) Attribute reads using sample-specific barcode (TSSV) Select highly accurate CCS reads (QV>=99, >= 5 passes) Trim reads (Cutadapt) and align to extended version or rcrs (BWA-MEM) Call positions alternative to the reference sequence in pileup (Rsamtools, min_nucleotide_depth=5) 13 MJA Weerts et al., Scientific Reports 2018; 8:2261

20 SMRT sequencing approach for mtdna Generate circular consensus reads (CCS2 algorithm) Attribute reads using sample-specific barcode (TSSV) Select highly accurate CCS reads (QV>=99, >= 5 passes) Trim reads (Cutadapt) and align to extended version or rcrs (BWA-MEM) Call positions alternative to the reference sequence in pileup (Rsamtools, min_nucleotide_depth=5) 13 MJA Weerts et al., Scientific Reports 2018; 8:2261

21 SMRT sequencing approach for mtdna 6 7 Generate circular consensus reads (CCS2 algorithm) Attribute reads using sample-specific barcode (TSSV) Select highly accurate CCS reads (QV>=99, >= 5 passes) Trim reads (Cutadapt) and align to extended version or rcrs (BWA-MEM) Call positions alternative to the reference sequence in pileup (Rsamtools, min_nucleotide_depth=5) compensate for mapping bias due to circularity of the mitochondrial genome 13 MJA Weerts et al., Scientific Reports 2018; 8:2261

22 SMRT sequencing approach for mtdna Generate circular consensus reads (CCS2 algorithm) Attribute reads using sample-specific barcode (TSSV) Select highly accurate CCS reads (QV>=99, >= 5 passes) Trim reads (Cutadapt) and align to extended version or rcrs (BWA-MEM) Call positions alternative to the reference sequence in pileup (Rsamtools, min_nucleotide_depth=5) 13 MJA Weerts et al., Scientific Reports 2018; 8:2261

23 Detection of low-frequent mtdna variants Determine detection limits empirically: mixtures of the cell lines containing different mtdna variants 14 MJA Weerts et al., Scientific Reports 2018; 8:2261

24 Detection of low-frequent mtdna variants Determine detection limits empirically: mixtures of the cell lines containing different mtdna variants Allele frequency Digital PCR 0/2 0/2 1/2 1/2 2/2 2/2 UltraSEEK 0/7 0/7 0/7 2/7 7/7 7/7 SMRT seq 0/23 0/23 0/23 21/23 23/23 23/23 positions 14 MJA Weerts et al., Scientific Reports 2018; 8:2261

25 Studying human mtdna in cancer HYPOTHESIS The high mtdna copy number per cell as well as the high mtdna mutation rate make it worthwhile to explore the potential of tumorspecific cf-mtdna variants as cancer marker in the blood of cancer patients 10

26 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Pilot study: SMRT sequence the entire mtdna of eight cancer patients Tumor, matched normal (tissue of origin!) and cfdna 15 MJA Weerts et al., Neoplasia. 2018; 20(7):

27 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Patient 1 Primary breast cancer (2 cm) with lymph nodes involved Neoadjuvant endocrine therapy Surgery, radiation Prolonged endocrine therapy Disease recurrence after 4 months (bone, lung and liver metastases) Second-line endocrine therapy Disease progression Chemotherapy 16 MJA Weerts et al., Neoplasia. 2018; 20(7):

28 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Patient 1 Primary breast cancer (2 cm) with lymph nodes involved Neoadjuvant endocrine therapy Surgery, radiation Prolonged endocrine therapy Disease recurrence after 4 months (bone, lung and liver metastases) Second-line endocrine therapy Disease progression Chemotherapy cfdna Primary tumor Normal mammary cfdna 16 MJA Weerts et al., Neoplasia. 2018; 20(7):

29 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Patient G>A 5351 A>G 7368 T>C 664 G>A 6255 G>A G>A T>C Primary tumor a Primary tumor b Mammary tissue a Mammary tissue b cfdna (pre surgery) cfdna (post chemo) MJA Weerts et al., Neoplasia. 2018; 20(7):

30 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Patient 1 Primary breast cancer (2 cm) with lymph nodes involved cfdna cfdna Neoadjuvant endocrine therapy Surgery, radiation Prolonged endocrine therapy Disease recurrence after 4 months (bone, lung and liver metastases) Second-line endocrine therapy Disease progression Chemotherapy cfdna Primary tumor Normal mammary cfdna cfdna 18 MJA Weerts et al., Neoplasia. 2018; 20(7):

31 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Patient G>A 5351 A>G 7368 T>C 664 G>A 6255 G>A G>A T>C Primary tumor a Primary tumor b Mammary tissue a Mammary tissue b cfdna (pre surgery) cfdna (pre 2nd line) na na na na na cfdna (pre chemo) na na na na na cfdna (post chemo) na na na 0 0 na na cfdna (post chemo) MJA Weerts et al., Neoplasia. 2018; 20(7):

32 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Patient 2 Colorectal cancer with synchronous hepatic metastases Surgery 20 MJA Weerts et al., Neoplasia. 2018; 20(7):

33 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Patient 2 Colorectal cancer with synchronous hepatic metastases cfdna Surgery Primary tumor Normal colon Metastasis 1 Normal liver Metastasis 2 20 MJA Weerts et al., Neoplasia. 2018; 20(7):

34 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Patient G>A T>C C>T G>A 66 G>T 152 T>C 189 A>G C>T C>T 1924 T>C 2305 T>C 4048 G>A T>C 60 T>C 72 T>C Primary tumor Metastasis Metastasis Colon tissue Liver tissue cfdna MJA Weerts et al., Neoplasia. 2018; 20(7):

35 Tumor-specific cf-mtdna variants as cancer marker in the blood of cancer patients? Patient G>A T>C C>T G>A 66 G>T 152 T>C 189 A>G C>T C>T 1924 T>C 2305 T>C 4048 G>A T>C 60 T>C 72 T>C Primary tumor Metastasis Metastasis Colon tissue Liver tissue cfdna MJA Weerts et al., Neoplasia. 2018; 20(7):

36 Studying human mtdna in cancer HYPOTHESIS The high mtdna copy number per cell as well as the high mtdna mutation rate make it worthwhile to explore the potential of tumorspecific cf-mtdna variants as cancer marker in the blood of cancer patients 22

37 Studying human mtdna in cancer HYPOTHESIS The high mtdna copy number per cell as well as the high mtdna mutation rate make it worthwhile to explore the potential of tumorspecific cf-mtdna variants as cancer marker in the blood of cancer patients CONCLUSIONS SMRT sequencing suitable for low-frequent mtdna single-nucleotide variant detection. Tumor-specific mtdna variants rarely detected as cfdna. 22

38 Acknowledgements Erasmus MC Medical Oncology John Foekens John Martens Stefan Sleijfer LUMC Leiden Genome Technology Center Rolf Vossen Yahya Anvar contact: Philips Research Precision and Decentralized Diagnostics Dianne van Strijp Pieter-Jan van der Zaag Eveline den Biezen Timmermans Anja van de Stolpe These results have been published Weerts et al., Scientific Reports 2018; 8:2261 Weerts et al., Neoplasia. 2018; 20(7):