QPCR A&E
QPCR A&E
QPCR A&E
• Introduction
– Mnemonics once again
– Mathematics of qPCR
• Data Analysis and Evaluation
– Quantification Strategies in qPCR
• Absolute Quantification
• Relative Quantification
– Fidelity in qPCR
• Specificity, Sensitivity, Accuracy, Reproducibility
• Experimental Variations, Replicates,
• Standard Deviation Calculations
• Optimizing qPCR Experiments
– Primer and probe optimization
– Multiplex assay optimization
• Questions & Answers
– I got mail ...
R CT ∆Rn
• Easy, but …
– Standards
• appropriate?
• same/different RT?
– Expensive
– Least accurate method
• quantitative accuracy = f(standards, standard curve)
Rainer B. Lanz, M.S., Ph.D. 7
Relative Quantification: RQ
• An Active Reference
– …is used to determine changes in the amount of a
given sample relative to another - internal - sample.
• a different amplicon in the same PCR reaction as the
amplification of the amplicon for the GOI
– Does not require standards with known concentrations
• Calculation Methods for Relative Quantitations
– Standard Curve method (∆CT)
• Two ‘standard’ curves (relative control & GOI)
• May include a 2nd normalization with an arbitrarily
chosen calibrator sample
– Comparative CT method (∆∆CT)
• no standards, but with amplification of a reference
• contingent upon similar amplification efficiencies of the
amplicons for GOI and reference
• Always relative to a calibrator sample
– Same amplicon:
• EA = EB ⇒ NA/NB = 2-∆CT
For example: if ∆CT between A and B is 5 cycles, then there
is 2-5 = 1/32 as much A than B.
– Different amplicons:
For example: GOI (x) and endogenous control (c):
• EX ≠ EC ⇒ Nx/Nc = K (1+Ec)CTc / (1+Ex)CTx
Rainer B. Lanz, M.S., Ph.D. 9
RQ: Mathematically
– NCT = N0 (1+E)CT and I = k NCT
– The relative Intensities of samples A and B is:
• IA = kA· NCTA = kA· N0A (1+EA)CTA and
• IB = kB· NCTB = kB· N0B (1+EB)CTB
– at threshold: IA = IB thus: kA· NCTA = kB· NCTB
– Solving for constants yields: K = kB/kA = NCTA/ NCTB ,
• inserting NCTA= N0A (1+EA)CTA and NCTB = N0B (1+EB)CTB and
rearranging we get:
– N0A/N0B = K· (1+EB)CTB / (1+EA)CTA (II)
• The fractions of A and B expressed as percentages are:
A = 100·[K·(1+EB)CTB/(1+EA)CTA] /1+K·[(1+EB)CTB/(1+EA)CTA]
B = 100·[1] /1+K·[(1+EB)CTB/(1+EA)CTA]
– Relative Standards:
• For example: the ratio of treatment (t) vs. control (c):
〈GOI〉
〈Ref〉
2nd normalization:
Calibrator = Brain
Liverc-myc
LiverGAPDH
Kidneyc-myc
KidneyGAPDH
(NA/ NB)t
=
(NA/ NB)c
SuperArray Bioscience
Corporation Newsletter 1
24 wells!
Relative standard
method: 48 wells!
EX vs. ER
Efficiencies:
|s| < 0.1
Livak and Schmittgen, 2001,
E max. Methods 25, 402-408
amplification
efficiency: – Comparing important linear regression plots for qPCR :
s = -3.32 Ampl. Efficiency ∆∆CT Validation
CT ∆CT
Log [] Log []
Rainer B. Lanz, M.S., Ph.D. 22
What If EGOI ≠ Eref ?
• Use Efficiency Correction
– Note: Rainer does NOT recommend efficiency corrections
for qPCR data analyses (if you follow the recommendations,
you most likely won’t have this problem)
(Ex)∆CT x(control-sample)
Relative N =
(ER)∆CT R(control-sample)
CV% = 100 SD
〈(1+E)-CT 〉
0.039 / 14.561 x 100 = 0.267%
CT < 36 • Variability
– In general, the intra-assay variation of 10-20% and a
Baseline mean inter-assay variation of 15-50% on molecule
Threshold basis is realistic over the wide dynamic range (of over
a billion fold range).
– Variability is highest at >107 and <102 template copy
ranges
• Cut-off value: cycle 35, i.e. disregard CT values for cycle
numbers 36 and higher.
– For the threshold methods, the precision is
dependent on the proper setting of the threshold,
which itself is dependent on proper base line settings.
– Evaluate only assays with a normal distribution of the
CTs
Rainer B. Lanz, M.S., Ph.D. 35
Integrated Genomics - The Future?
• Real-Time StatMinerTM
– http://www.integromics.com/StatMiner.php
• Maximizing PCR
– Optimizing primer concentrations :-|
• ∆Rn > 2.5 logs (∆Rn = f{quencher})
– Optimizing probe concentrations ;-(
– Optimizing AD experiments ✓
– Optimizing multiplex experiments
• Almost always! ✓
Optimizing Primer Concentrations
• Primer Optimization Matrix
– Maximize ∆Rn :
– Suggested conc.:
• ≤900nM for TaqMan
• ~50nM for SYBR Green
Primer Probe
[nM] [nM]
100/900 50
100/900 125
100/900 250
100/900 500
– Suggested conc.:
• ~250nM
… In your classes you said that we should ignore Ct values larger than
35. Why? I have seen plots that go up to 45 cycles, or higher.
CT 35 is close to the limit of detection for threshold-
based qPCR . By definition, a gene is not detectable
When the average CT > 35.
3) You mentioned in the first class about AmpErase UNG and ROX.
Which one would you (strongly) recommend that we include in our assay
design?
ROX
Rainer B. Lanz, M.S., Ph.D. 46
Hi Dr. Lanz,
These are a few questions I have:
1) If we want to use the 18s rRNA as our endogenous control, could we use
the fabricated arrays like the QuantumRNA™ Classic 18S from Applied
Biosystems? Or is there a particular one you would recommend?
Isn’t this a newer version of Ambion’s Internal Standards?
Do you need RNA standards or 18S for a endogenous
control? I’ve always used the original ABI primers for 18S
X00686: 5'accgcagctaggaataatgga3' 5'gcctcagttccgaaaacca3'
2) In the class today, do you think you could include an example about how to
the data analysis with standard curve, as well as how to normalize data
between plates and at the same time use an endogenous control as
GAPDH?
Done!
3) I run out of my plate-to-plate controls, so I had to make up a new set of
controls. I run a plate with old and new controls (measuring GOI and
endogenous). So, from then on I have used the new controls, but I’m
comparing to data where I was using the old controls. Can you do a new-
control/old-control ratio and then normalize all you data with that number? Is
kind of complicated!
Relative (intra-plate) quantities ⇒ always comparable!
Rainer B. Lanz, M.S., Ph.D. 47
I am taking your qPCR classes right now and have thought of a couple of
questions.
First, when running the delta delta Ct relative quantification you have said
that the Es must be similar. Would it still be possible to analyze data for Es
that are not similar as long as you know both E values? I would think (at
least theoretically) you could calculate the difference between the Es and
know how big of a change to discount in your analysis.
Yes, but in general I don’t recommend this approach.
Rather, I would re-think the assay.
Remember that 0.1 difference in the E values results in
a 5-fold difference in sample quantity at PCR cycle 30.
If the differences in E is due to a ‘not-so-good’ 18S
amplicon, I do have some suggestions (>discussion)
Second, does the primer design change much when designing them for
C. elegans? Since C. elegans have very small introns, I did not know how
this would affect my design
(trying to space across 2 exons).
Put a hybridization probe over the exon junction.