Successful gene expression studies using validated qPCR assays Jan - - PowerPoint PPT Presentation

successful gene expression studies using validated qpcr
SMART_READER_LITE
LIVE PREVIEW

Successful gene expression studies using validated qPCR assays Jan - - PowerPoint PPT Presentation

Successful gene expression studies using validated qPCR assays Jan Hellemans, CEO Biogazelle webinar October 28 th , 2015 Agenda Requirements for high quality qPCR assays Approaches for qPCR assay validation How good a result


slide-1
SLIDE 1

Successful gene expression studies using validated qPCR assays

Jan Hellemans, CEO Biogazelle webinar October 28th, 2015

slide-2
SLIDE 2

Agenda

  • Requirements for high quality qPCR assays
  • Approaches for qPCR assay validation
  • How good a result can be achieved with highly optimized design

tools

  • How to easily discover additional genes of interest
  • Best practice in qPCR
slide-3
SLIDE 3

Online poll

How would you rank your current knowledge about gene expression assay design? a) basic b) advanced c) expert

slide-4
SLIDE 4

qPCR is reference technology for nucleic acid quantification

  • sensitivity and specificity
  • wide dynamic range
  • speed
  • relatively low cost
  • conceptual and practical simplicity

qPCR is easy to perform ≠ easy to do it right

  • many steps involved
  • all need to be right

Introduction

slide-5
SLIDE 5

R

relative quantification quality control statistical analysis

C

Prepare – cycle – report

P

experiment design samples assays

prepare cycle report

slide-6
SLIDE 6

R

relative quantification quality control statistical analysis

C

Prepare – cycle – report

P

experiment design samples assays

prepare cycle report

slide-7
SLIDE 7

Assay design & validation

design

  • amplicon length
  • primer positions (exonic or intron-spanning)

Considerations & best practices

slide-8
SLIDE 8

Assay design & validation

design

  • amplicon length
  • primer positions (exonic or intron-spanning)

Considerations & best practices

gene exonic intron-spanning

slide-9
SLIDE 9

Assay design & validation

design

  • amplicon length
  • primer positions (exonic or intron-spanning)
  • transcript coverage

Considerations & best practices

gene transcript 1 transcript 2 transcript 3

2 3 1 2 3 2

coverage

slide-10
SLIDE 10

Assay design & validation

design

  • amplicon length
  • primer positions (exonic or intron-spanning)
  • transcript coverage

in silico verification

  • specificity prediction (retropseudogenes and other

homologues)

  • secondary structure analysis

wet lab validation (experimental)

  • specificity assessment (gel, melt, amplicon sequencing)
  • Cq of NTC (for SYBR assays)
  • amplification efficiency determination (slope, E, SE(E), r2)

Considerations & best practices

slide-11
SLIDE 11

Assay design & validation

design

  • amplicon length
  • primer positions (exonic or intron-spanning)
  • transcript coverage

in silico verification

  • specificity prediction (retropseudogenes and other

homologues)

  • secondary structure analysis

wet lab validation (experimental)

  • specificity assessment (gel, melt, amplicon sequencing)
  • Cq of NTC (for SYBR assays)
  • amplification efficiency determination (slope, E, SE(E), r2)

Considerations & best practices

slide-12
SLIDE 12

Properties of the perfect assay

  • specific for the gene of interest => no off-target amplification
  • detection of all transcript variants
  • detection not affected by polymorphisms => no allelic bias or

drop out

  • amplification efficiency ~100%
  • no gDNA co-amplification
  • no primer dimer formation
slide-13
SLIDE 13

The perfect assay

slide-14
SLIDE 14

Online poll

How do you currently obtain your ‘perfect’ qPCR assay? a) using your own home brewed assays b) buying pre-designed assays (commercial) c) currently not designing any gene expression assays

slide-15
SLIDE 15

The perfect assay

  • For some genes, there is no perfect assay
  • no unique sequence (homology with other genes –

pseudogenes)

... or the best possible

slide-16
SLIDE 16

Gene homology in olfactory receptor genes prevents perfect designs

1532 / 2043 (75%) of genes without perfect design have homologous genes that differ less than 12.5% (2 variations per 16 bases)

0% 5% 10% 15% 20% 25% 30% 35% 40% 45% 50%

1 101 201 301 401 501 601 701 801 901 1001 1101 1201 1301 1401 1501 1601 1701 1801 1901 2001

distances (clustalW) between all genes without perfect design

slide-17
SLIDE 17

The perfect assay

  • For some genes, there is no perfect assay
  • no unique sequence (homology with other genes –

pseudogenes)

  • no common sequence among all transcripts
  • regions are excluded because of repeats, secondary

structures, SNPs, homology, ...

  • Make the best possible compromise and report potential issues
  • Design  in silico quality control  wet lab validation

... or the best possible

slide-18
SLIDE 18

Assay design using primerXL

  • database of genomic information (transcripts, SNPs, ...)
  • tools for target region selection (maximize transcript coverage)
  • primer3 design engine
  • analysis of secondary structures and SNPs in primer & probe

annealing regions

  • specificity prediction (BiSearch, bowtie)
  • relaxation cascade (from perfect to best possible)
slide-19
SLIDE 19

Impact of primer mismatches on qPCR assay performance

Lefever, Clin Chem 2013

slide-20
SLIDE 20

BiSearch specificity prediction

  • BiSearch loose
  • 1222222222222222
  • BiSearch strict
  • 1233333333333
slide-21
SLIDE 21

BiSearch specificity prediction

  • BiSearch loose
  • 1222222222222222
  • nly the gene of interest

(FFAR2)

  • BiSearch strict
  • 1233333333333

reads seq gene_list

  • fficial_symbol

location 2843 CATGGCAGTCACCATCTTCTGCTACTGGCGTTTTGTGTGGATCATGCTCTCCCAGCCC

CTTGTGGGGGCCCAGAGGCGGCGCCGAGCCGTGGGGCTGGCTGTGGTGACGC TGCTCAATTTCCTGGTGTGCTTCGGACCTTACAGATCGGAA

ENSG00000126262 FFAR2 19:35940617- 35942667 1897 GTAAGGTCCGAAGCACACCAGGAAATTGAGCAGCGTCACCACAGCCAGCCCC

ACGGCTCGGCGCCGCCTCTGGGCCCCCACAAGGGGCTGGGAGAGCATGATCC ACACAAAACGCCAGTAGCAGAAGATGGTGACTGCCATGAGATCGGAA

ENSG00000126262 FFAR2 19:35940617- 35942667 1535 GTAAGGTCCGAAGCACACCGAGAGCTGGGAGCAGGAGCTACACAGTCTGCTGG

CCTCACTGCACACCCTGCTGGGGGCCCTGTACGAGGGAGCAGAGACTGCTCCT GTGCAGAATGAAGGCCCTGGGGTGGAGATGCTGCTGTCCTCAGAA

ENSG00000141456 AC091153.1 17:4574680- 4607632 1097 CATGGCAGTCACCATCTTCTGAGGACAGCAGCATCTCCACCCCAGGGCCTTCATT

CTGCACAGGAGCAGTCTCTGCTCCCTCGTACAGGGCCCCCAGCAGGGTGTGCA GTGAGGCCAGCAGACTGTGTAGCTCCTGCTCCCAGCTCTCGG

ENSG00000141456 AC091153.1 17:4574680- 4607632 1091 CATGGCAGTCACCATCTTCTGAGGACAGCAGCATCTCCACCCCAGGGCCTTCATT

CTGCACAGGAGCAGTCTCTGCTCCCTCGTACAGGGCCCCCAGCAGGGTGTGCA GTGAGGCCAGCAGACTGTGTAGCTCCTGCTCCCAGCTCTCGGT

ENSG00000141456 AC091153.1 17:4574680- 4607632

slide-22
SLIDE 22

Wet lab validation

PCR composition

  • total volume: 5 μl
  • instrument: Bio-Rad CFX384 (with CFX Automation System)
  • mastermix: Bio-Rad SsoAdvanced SYBR
  • primer conc: 250 nM each

PCR program

  • default cycling protocol for SsoAdvanced SYBR (Ta=60°C)

Samples

  • cDNA: 25 ng (total RNA equivalents – Agilent Universal human reference

RNA = MAQC A)

  • gDNA: 2.5 ng (Roche)
  • NTC: water + carrier (5 ng/μl yeast transfer RNA)
  • synthetic template (pooled 60-mers in concentration range: 20 M – 20

copies)

setup

slide-23
SLIDE 23

Wet lab validation

  • lab validation of 103 053 assays

(human, mouse and rat coding genes)

  • 1 456 142 reactions
  • 3 822 PCR plates (384-well)
  • equivalent to 15 288 PCR plates (96-well)

some numbers

305 m

slide-24
SLIDE 24

Amplification efficiency

  • initial publication: Vermeulen et al., Nucleic Acids Research, 2009
  • Biogazelle approach (easy & cost effective)
  • 60-mer
  • no modifications, standard desalted
  • 7 points dilution series: 20 000 000 > 20 molecules
  • equivalent to full length double stranded template
  • limitation: behavior of first cycles amplifying from cDNA are not

evaluated

synthetic templates

30 nt 3’ 30 nt 5’

ds template ss oligo r2<0.99 1 1 median E 2.00 2.01 average E 2.00 2.01 count E <> [1.90-2.10] 1 3 paired t-test p-value 0.14

slide-25
SLIDE 25

Amplification efficiency

distribution (n = 50 133)

89%

slide-26
SLIDE 26

Amplification efficiency

distribution (n = 50 133)

89% redesign redesign

slide-27
SLIDE 27

Specificity

amplicon sizing ( + melt analysis for SYBR assays)

  • limited sensitivity for detecting low level non-specific

coamplification

  • failure to observe non-specific amplification of sequences

with similar size and/or Tm e.g. expressed pseudogenes or homologous genes next level of specificity assessment

  • in silico specificity predictions by BiSearch
  • massively parallel sequencing of pooled PCR products
  • average coverage > 1000-fold 

lab specificity > 99.9%

  • 50 – 200 times more sensitive than size analysis and Sanger

sequencing

NGS for increased sensitivity

slide-28
SLIDE 28

Specificity

most assays are 100% on-target

slide-29
SLIDE 29

Specificity

0% 25% 50% 75% 100% % on-target

2/3 of non-specific assays may go unnoticed without NGS

0% 20% 40% 60%

0 < x < 0.1 0.1 < x < 0.2 0.2 < x < 0.3 0.3 < x < 0.4 0.4 < x < 0.5 0.5 < x < 0.6 0.6 < x < 0.7 0.7 < x < 0.8 0.8 < x < 0.9 0.9 < x < 1

slide-30
SLIDE 30

Specificity

perfect 60 293 86% acceptable (<10% non-specific) 5 866 8% predicted non-specificity (no specific design found) 1 204 2% failing specificity QC criteria 2 467 4%

the power of in silico verification

slide-31
SLIDE 31

Online poll

How do you validate your assay’s specificity? a) Melt curves b) Size analysis (gel or capillary) c) Restriction digestion with gel analysis d) Sequencing of PCR products

slide-32
SLIDE 32

Online poll

Do you know the MIQE initiative? a) Yes b) No

slide-33
SLIDE 33

MIQE compliant PrimePCR assay

validation data sheet for human, mouse & rat

slide-34
SLIDE 34
  • building on the confidence validated on > 100,000 assays
  • skip wet lab validation
  • 9 organisms
  • assays for SYBR or with probe

27,155 19,762 19,310 25,006 15,307 20,184 19,049 6,572 21,360

PrimePCR assay for 9 extra organisms

slide-35
SLIDE 35

Screening vs targeted

Transcriptome

  • microarray or RNA-seq
  • hypothesis generating
  • few samples (high cost per sample)

Gene panels

  • qPCR (high sensitivity, low cost)
  • selection of genes based on pathway, cell type or disease
  • convenient balance between blind screening and targeted

analysis Individual gene(s)

  • qPCR
  • high flexibility, customization
  • low cost per sample
slide-36
SLIDE 36

0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% ribo-depletion RNAseq - detection poly-A RNAseq - detection qPCR - detection ribo-depletion RNAseq - quantification poly-A RNAseq - quantification qPCR - quantification

Saturation analysis

MAQC A - % of PrimePCR

slide-37
SLIDE 37

Questions?

slide-38
SLIDE 38

Thank you for your attention!