حاسبة كفاءة التفاعل التسلسلي للبوليميراز الكمي: أداة تحليل المنحنى القياسي

حاسبة كفاءة تفاعل PCR الكمي (qPCR)

تقدّر حاسبة كفاءة تفاعل PCR الكمي مدى جودة نسخ تفاعل PCR الكمي (qPCR) لتسلسل DNA المستهدف مع كل دورة تضخيم. وتستخدم قيم Ct، أي رقم الدورة التي تتجاوز فيها الإشارة الفلورية عتبة محددة، من سلسلة تخفيف لحساب الكفاءة والميل والجزء المقطوع وR².

ما كفاءة تفاعل PCR الكمي؟

يتتبع تفاعل PCR الكمي التألق أثناء نسخ التفاعل لتسلسل DNA مستهدف، دورة بعد أخرى. وتصف كفاءة التضخيم مقدار نمو كمية DNA المستهدف مع كل دورة. عند كفاءة قدرها 100%، تتضاعف كمية DNA المستهدف في كل دورة.

نادرًا ما تصل التفاعلات الفعلية إلى 100% تمامًا. وتعد معظم المختبرات الكفاءة بين 90% و110% مقبولة لإجراء قياس كمي موثوق. وتشير قيمة أقل من 90% أو أعلى من 110% إلى مشكلة، مثل سوء تصميم البادئات، أو مثبطات PCR المحمولة من استخلاص العينة، أو أخطاء السحب بالماصة في سلسلة التخفيف.

تهم الكفاءة بدرجة أكبر عند مقارنة التعبير الجيني بين العينات. وتفترض طريقة ΔΔCt الشائعة لحساب التغير النسبي أن كل تفاعل يعمل بكفاءة قدرها 100%. وعندما تختلف الكفاءة المقاسة عن 100%، أو عندما تختلف كفاءة الجين المستهدف عن كفاءة الجين المرجعي، قد تكون نتائج التغير النسبي غير دقيقة.

معادلة كفاءة تفاعل PCR الكمي

تُحسب الكفاءة من ميل المنحنى القياسي، وهو خط يتكون من رسم قيم Ct مقابل اللوغاريتم ذي الأساس 10 (log10) لكمية DNA الابتدائية في سلسلة التخفيف.

E = 10^(-1/الميل) - 1

تمثل E الكفاءة، وتُعبَّر عنها في صورة عدد عشري (تُضرب في 100 للحصول على النسبة المئوية). والميل هو ميل خط المنحنى القياسي. ويكون سالبًا، لأن Ct يرتفع مع انخفاض كمية DNA الابتدائية.

يقابل ميل قدره -3.32 كفاءة قدرها 100%:

10^(-1/-3.32) - 1 = 10^0.301–1 ≈ 1.00

وتصح هذه النتيجة لأن التخفيف بمقدار 10 أضعاف ينبغي أن يغيّر Ct بنحو 3.32 دورة إذا كان التفاعل يضاعف DNA في كل دورة، إذ إن log2(10) يساوي تقريبًا 3.32.

عمليًا، يُعد الميل الواقع بين -3.1 و-3.6 مقبولًا. ويقابل هذا النطاق كفاءة تبلغ تقريبًا من 90% إلى 110%. ويعطي الميل الأقرب إلى -3.1 كفاءة محسوبة أعلى، بينما يعطي الميل الأقرب إلى -3.6 كفاءة أقل.

كيفية حساب كفاءة تفاعل PCR الكمي: طريقة المنحنى القياسي

1. حضّر سلسلة تخفيف. ابدأ بالعينة الأعلى تركيزًا وخففها بعامل ثابت، غالبًا 10 أضعاف، لكل عينة تالية. أجرِ تفاعل PCR الكمي على كل تخفيف وسجّل قيمة Ct الخاصة به.
2. حوّل كل تخفيف إلى قيمة log10. وباعتبار العينة الأولى، وهي الأعلى تركيزًا، في الموضع 0، تكون قيمة log10(الكمية النسبية) للعينة رقم i مساوية لـ -i × log10(عامل التخفيف).
3. ارسم Ct (على المحور الرأسي) مقابل قيمة log10 (على المحور الأفقي) لكل تخفيف، وأجرِ انحدارًا خطيًا.
4. استخرج الميل والجزء المقطوع وR² لخط الانحدار.
5. عوّض عن الميل في E = 10^(-1/الميل) - 1 للحصول على الكفاءة.

معايير جودة المنحنى القياسي

• قيمة R² لا تقل عن 0.98. يقيس ذلك مدى اقتراب النقاط من خط مستقيم. وغالبًا ما تشير القيمة الأقل إلى أخطاء السحب بالماصة أو عدم اتساق التضخيم.
• ميل بين -3.1 و-3.6. ويقابل هذا النطاق كفاءة قدرها 90–110%.
• ثلاث نقاط تخفيف على الأقل. وتعطي خمس أو ست نقاط خطًا أكثر موثوقية، كما تسهّل اكتشاف القيمة الشاذة.
• ينبغي أن تعطي العينة الأعلى تركيزًا قيمة Ct تقارب 15-20. وينبغي أن تبقى العينة الأكثر تخفيفًا تحت Ct البالغ 35، لأن القيم الأعلى من ذلك تقترب من حد الكشف.

مثال محلول

تعطي سلسلة تخفيف من خمس نقاط، بمقدار 10 أضعاف، قيم Ct التالية، من الأعلى إلى الأقل تركيزًا: 15.0، 18.5، 22.0، 25.5، 29.0.

قيم log10(الكمية النسبية) هي 0، -1، -2، -3، -4، لأن كل خطوة تخفيف تضيف -log10(10)، أي -1. ويرتفع Ct بمقدار 3.5 دورة في كل خطوة، لذلك تقع النقاط على خط مستقيم، ويعطي الانحدار الخطي ما يلي:

• الميل: -3.5
• الجزء المقطوع: 15.0
• R²: 1.0

الكفاءة: E = 10^(-1/-3.5) - 1 = 10^0.286–1 ≈ 0.931، أي نحو 93%.

كفاءة قدرها 93% مع R² = 1.0 تستوفي معايير الجودة القياسية، وتناسب قياس التعبير الجيني كميًا.

مثال على الشيفرة

تتبع دالة JavaScript هذه الطريقة نفسها التي تتبعها الحاسبة. ويُعكس مؤشر كل قيمة Ct قبل حساب اللوغاريتم، لأن العينة الأولى، وهي الأعلى تركيزًا، لها أعلى كمية نسبية، وكل عينة لاحقة تكون أكثر تخفيفًا.

1function qpcrEfficiency(ctValues, dilutionFactor) {
2  // Sample 0 is most concentrated; each later sample is diluted further,
3  // so its log10(relative quantity) is negative.
4  const x = ctValues.map((_, i) => -i * Math.log10(dilutionFactor));
5  const n = ctValues.length;
6  let sumX = 0, sumY = 0, sumXY = 0, sumXX = 0;
7  for (let i = 0; i < n; i++) {
8    sumX += x[i];
9    sumY += ctValues[i];
10    sumXY += x[i] * ctValues[i];
11    sumXX += x[i] * x[i];
12  }
13  const slope = (n * sumXY - sumX * sumY) / (n * sumXX - sumX * sumX);
14  const intercept = (sumY - slope * sumX) / n;
15  const efficiency = (Math.pow(10, -1 / slope) - 1) * 100;
16  return { efficiency, slope, intercept };
17}
18
19qpcrEfficiency([15.0, 18.5, 22.0, 25.5, 29.0], 10);
20// => { efficiency: 93.07, slope: -3.5, intercept: 15 }
21

لماذا تهم الكفاءة في التعبير الجيني؟

تفترض حسابات التغير النسبي، بما فيها طريقة ΔΔCt، أن كل تفاعل يضاعف هدفه في كل دورة. وعندما تختلف الكفاءتان المقاسوَتان لجين مستهدف وجين مرجعي، تُدخل معادلة ΔΔCt البسيطة خطأً في النتيجة. وغالبًا ما ينتقل الباحثون إلى حساب مصحح للكفاءة، مثل الطريقة التي وصفها مايكل بفافل في 2001، عندما تختلف الكفاءتان بأكثر من نحو خمس نقاط مئوية. وتوصي إرشادات MIQE، وهي مجموعة واسعة الاستشهاد من معايير الإبلاغ عن تفاعل PCR الكمي المنشورة في 2009، بأن تورد الدراسات الكفاءة المحسوبة والميل وR² إلى جانب نتائج القياس الكمي.

الأسئلة الشائعة

ما النسبة المئوية الجيدة لكفاءة تفاعل PCR الكمي؟

تُقبل عمومًا قيمة بين 90% و110%، بينما تحقق الاختبارات المحسّنة جيدًا نحو 95–105%. ويقابل ميل قدره -3.32 كفاءة تساوي تمامًا 100%.

لماذا تتجاوز كفاءة تفاعل PCR الكمي لدي 110%؟

تشمل الأسباب الشائعة أخطاء السحب بالماصة في سلسلة التخفيف، أو ثنائيات البادئات، أو المثبطات التي تؤثر في العينات المركزة والمخففة بدرجات غير متساوية، أو ضبط عتبة التألق خارج الطور الأسي للتضخيم. وقد يشير منحنى الانصهار ذو القمة المتعددة إلى وجود ثنائيات البادئات.

لماذا تقل كفاءة تفاعل PCR الكمي لدي عن 90%؟

تشمل الأسباب المحتملة سوء تصميم البادئات، أو مثبطات PCR المحمولة من استخلاص العينة، أو درجة حرارة التهجين غير المثلى. وقد يساعد التحقق من نوعية البادئات وإجراء تدرج حراري في تحديد السبب.

ماذا تعني قيمة R² المنخفضة؟

تقيس R² مدى توافق قيم Ct مع خط مستقيم مقابل الكمية اللوغاريتمية. وتشير قيمة أقل من 0.98 إلى عدم اتساق السحب بالماصة، أو عينة مثبطة، أو نقطة تخفيف قريبة من حد الكشف.

كم عدد نقاط التخفيف التي يحتاج إليها المنحنى القياسي؟

ثلاث نقاط على الأقل، وهو الحد الأدنى اللازم لإجراء انحدار خطي. وتعطي خمس أو ست نقاط ميلًا أكثر موثوقية، كما تسهّل اكتشاف القيم الشاذة.

هل يمكن أن تتغير الكفاءة بين التجارب؟

نعم. فقد يغيّر تشغيلة جديدة من الكواشف، أو دفعة جديدة من البادئات، أو نوع مختلف من العينات، أو جهاز qPCR مختلف الكفاءة. ويعيد كثير من المختبرات فحص الكفاءة كل بضعة أشهر أو مع كل تشغيلة جديدة من الكواشف.

المراجع

1. بوستن SA، بينيش V، غارسون JA، وآخرون. إرشادات MIQE: الحد الأدنى من المعلومات لنشر تجارب PCR الكمي الآني. الكيمياء السريرية. 2009؛55(4):611-622.
2. بفافل MW. نموذج رياضي جديد للقياس الكمي النسبي في RT-PCR الآني. أبحاث الأحماض النووية. 2001؛29(9):e45.
3. سفِتس D، تيخوباد A، نوفوسادوفا V، بفافل MW، كوبيستا M. ما مدى جودة تقدير كفاءة PCR: توصيات لإجراء تقييمات دقيقة وموثوقة لكفاءة qPCR. الكشف الجزيئي والقياس الكمي. 2015؛3:9-16.