Hesabu ufanisi wa PCR kutoka kwa thamani za Ct na sababu za kupunguza. Changanua mipangilio ya kiwango, tambua ufanisi wa uimarishaji, na thibitisha majaribio yako ya qPCR ya kiasi.
Thamani lazima iwe chanya
Thamani lazima iwe chanya
Thamani lazima iwe chanya
Thamani lazima iwe chanya
Thamani lazima iwe chanya
Ingiza data sahihi ili kuunda chati
Ufanisi wa qPCR ni kipimo cha jinsi mchakato wa PCR unavyofanya kazi. Ufanisi wa 100% unamaanisha kuwa kiasi cha bidhaa ya PCR kinadouble kwa kila mzunguko wakati wa awamu ya kuongezeka.
Ufanisi unakokotwa kutoka kwa mteremko wa mwelekeo wa kawaida, ambao unapatikana kwa kupanga thamani za Ct dhidi ya logarithm ya mkusanyiko wa awali wa sampuli (mfululizo wa mabadiliko).
Ufanisi (E) unakokotwa kwa kutumia formula:
E = 10^(-1/slope) - 1
Ufanisi wa Quantitative Polymerase Chain Reaction (qPCR) ni kipimo muhimu ambacho kinahusiana moja kwa moja na usahihi na uaminifu wa majaribio yako ya qPCR. Kihesabu ufanisi wa qPCR husaidia watafiti kubaini jinsi vizuri mchakato wao wa PCR unavyoongeza mfuatano wa DNA lengo na kila mzunguko wa joto. Mchakato mzuri wa qPCR unapaswa kuwa na ufanisi kati ya 90-110%, ukionyesha kwamba kiasi cha bidhaa ya PCR kinadoubishwa kwa karibu na kila mzunguko wakati wa awamu ya kuongezeka.
Ufanisi duni wa kuimarisha unaweza kusababisha upimaji usio sahihi, matokeo yasiyoaminika, na hitimisho potofu ya majaribio. Kwa kuhesabu na kufuatilia ufanisi wako wa qPCR, unaweza kuboresha hali za mchakato, kuthibitisha muundo wa primer, na kuhakikisha ubora wa data yako ya quantitative PCR.
Kihesabu hiki kinatumia njia ya curve ya kiwango, ambayo inachora thamani za mzunguko wa mabadiliko (Ct) dhidi ya logarithm ya mkusanyiko wa sampuli (inayoonyeshwa na mfululizo wa upunguzaji), ili kubaini ufanisi wa jaribio lako la qPCR. Mteremko unaotokana na curve hii ya kiwango unatumika kisha kuhesabu ufanisi wa kuimarisha kwa kutumia formula rahisi ya kihesabu.
Ufanisi wa mchakato wa qPCR unahesabiwa kutoka kwa mteremko wa curve ya kiwango kwa kutumia formula ifuatayo:
Ambapo:
Kwa mchakato mzuri wa PCR wenye ufanisi wa 100% (kuongezeka kwa bidhaa za amplicon kwa ukamilifu na kila mzunguko), mteremko ungekuwa -3.32. Hii ni kwa sababu:
10^{(-1/-3.32)} - 1 = 10^{0.301} - 1 = 2 - 1 = 1.0 \text{ (au 100%)}
Asilimia ya ufanisi inahesabiwa kwa kuzingatia ufanisi wa desimali kwa 100:
\text{Ufanisi (%)} = E \times 100\%
Curve ya kiwango inaundwa kwa kuchora thamani za Ct (axisi ya y) dhidi ya logarithm ya mkusanyiko wa awali wa sampuli au kipengele cha upunguzaji (axisi ya x). Uhusiano kati ya hizi mbili unapaswa kuwa wa moja kwa moja, na ubora wa uhusiano huu wa moja kwa moja unakaguliwa kwa kutumia kipimo cha uamuzi (R²).
Kwa ufanisi wa kuaminika wa qPCR:
Maandalizi ya data: Kihesabu kinachukua thamani zako za Ct kwa kila pointi ya upunguzaji na kipengele cha upunguzaji kama ingizo.
Mabadiliko ya log: Mfululizo wa upunguzaji unabadilishwa kuwa kiwango cha logarithmic (log msingi 10).
Uchambuzi wa mwelekeo: Kihesabu kinafanya uchambuzi wa mwelekeo wa moja kwa moja kwenye data iliyobadilishwa na log ili kubaini mteremko, y-intercept, na thamani ya R².
Hesabu ya ufanisi: Kwa kutumia thamani ya mteremko, ufanisi unahesabiwa kwa kutumia formula E = 10^(-1/slope) - 1.
Ufafanuzi wa matokeo: Kihesabu kinatoa ufanisi kama asilimia, pamoja na mteremko na thamani ya R² ili kukusaidia kutathmini uaminifu wa jaribio lako la qPCR.
Fuata hatua hizi ili kuhesabu ufanisi wako wa qPCR:
Weka idadi ya upunguzaji: Chagua ni pointi ngapi za upunguzaji unao katika curve yako ya kiwango (kati ya pointi 3-7 zinapendekezwa).
Ingiza kipengele cha upunguzaji: Ingiza kipengele cha upunguzaji kilichotumika kati ya sampuli zinazofuatana (kwa mfano, 10 kwa mfululizo wa upunguzaji wa 10-fold, 5 kwa mfululizo wa upunguzaji wa 5-fold).
Ingiza thamani za Ct: Ingiza thamani za Ct kwa kila pointi ya upunguzaji. Kawaida, upunguzaji wa kwanza (Upunguzaji 1) unakuwa na mkusanyiko wa juu wa sampuli, ukisababisha thamani ya Ct kuwa ya chini zaidi.
Tazama matokeo: Kihesabu kitafanya hesabu kiotomatiki na kuonyesha:
Tafsiri matokeo: Kadiria ikiwa ufanisi wako wa qPCR uko ndani ya upeo wa kukubalika (90-110%) na ikiwa thamani ya R² inaonyesha curve ya kiwango ya kuaminika (≥ 0.98).
Nakili matokeo: Tumia kitufe cha "Nakili Matokeo" ili kunakili thamani zote zilizohesabiwa kwa rekodi zako au machapisho.
Hebu tufanye mfano:
Wakati wa kuchora kwenye curve ya kiwango:
Kihesabu kitafanya uchambuzi wa mwelekeo na kubaini:
Kwa kutumia formula ya ufanisi:
Hii inaonyesha ufanisi mzuri wa qPCR wa 93%, ambao uko ndani ya upeo wa kukubalika wa 90-110%.
Kabla ya kutumia jozi mpya za primer kwa majaribio ya quantitative, ni muhimu kuthibitisha utendaji wake. Kuwa na ufanisi wa qPCR husaidia:
Wakati wa kuendeleza majaribio mapya ya qPCR, hesabu za ufanisi ni muhimu kwa:
Katika majaribio ya upimaji wa jamaa, kujua ufanisi wa PCR ni muhimu kwa:
Katika mazingira ya kliniki na uchunguzi, ufanisi wa qPCR ni muhimu kwa:
Kwa maombi ya usalama wa mazingira na chakula, hesabu za ufanisi husaidia:
Ingawa njia ya curve ya kiwango ndiyo njia inayotumika zaidi kwa kuhesabu ufanisi wa qPCR, kuna mbinu mbadala:
Njia hii inahesabu ufanisi kutoka kwa data ya fluorescence ya mzunguko mmoja wa kuimarisha, bila kuhitaji mfululizo wa upunguzaji. Programu kama LinRegPCR inachambua awamu ya kuongezeka ya majibu binafsi ili kubaini ufanisi.
Faida:
Hasara:
Digital PCR (dPCR) inatoa upimaji wa kiasi bila kuhitaji curve ya kiwango au hesabu za ufanisi.
Faida:
Hasara:
Baadhi ya programu za uchambuzi wa qPCR hutoa njia za kulinganisha kiasi ambazo zinakadiria ufanisi bila curve kamili ya kiwango.
Faida:
Hasara:
Maendeleo ya qPCR na hesabu za ufanisi yamekua kwa kiasi kikubwa katika miongo michache iliyopita:
Mchakato wa Polymerase Chain Reaction (PCR) ulivumbuliwa na Kary Mullis mwaka 1983, ukarevolutionize biolojia ya molekuli. Hata hivyo, PCR ya jadi ilikuwa ya ubora au ya nusu-kiasi tu. Mfumo wa kwanza wa qPCR wa wakati halisi ulitengenezwa mwanzoni mwa miaka ya 1990 na Russell Higuchi na wenzake, ambao walionyesha kwamba kufuatilia bidhaa za PCR kadri zinavyokusanywa (kwa kutumia fluorescence ya ethidium bromide) kunaweza kutoa habari ya kiasi.
Kadri teknolojia ya qPCR ilivyoendelea, watafiti walitambua umuhimu wa viwango na uthibitishaji. Wazo la ufanisi wa PCR likawa katikati ya upimaji wa kuaminika:
Uwanja umeendelea kuimarika na:
Leo, kuhesabu na kuripoti ufanisi wa qPCR inachukuliwa kuwa muhimu kwa kuchapisha data za qPCR zinazoaminika, na zana kama hii ya kihesabu husaidia watafiti kuzingatia mbinu bora katika uwanja.
1' Formula ya Excel kwa kuhesabu ufanisi wa qPCR kutoka kwa mteremko
2' Weka kwenye seli B2 ikiwa mteremko uko kwenye seli A2
3=10^(-1/A2)-1
4
5' Formula ya Excel kubadilisha ufanisi kuwa asilimia
6' Weka kwenye seli C2 ikiwa ufanisi wa desimali uko kwenye seli B2
7=B2*100
8
9' Kazi ya kuhesabu ufanisi kutoka kwa thamani za Ct na kipengele cha upunguzaji
10Function qPCR_Efficiency(CtValues As Range, DilutionFactor As Double) As Double
11 Dim i As Integer
12 Dim n As Integer
13 Dim sumX As Double, sumY As Double, sumXY As Double, sumXX As Double
14 Dim logDilution As Double, slope As Double
15
16 n = CtValues.Count
17
18 ' Hesabu mwelekeo wa moja kwa moja
19 For i = 1 To n
20 logDilution = (i - 1) * WorksheetFunction.Log10(DilutionFactor)
21 sumX = sumX + logDilution
22 sumY = sumY + CtValues(i)
23 sumXY = sumXY + (logDilution * CtValues(i))
24 sumXX = sumXX + (logDilution * logDilution)
25 Next i
26
27 ' Hesabu mteremko
28 slope = (n * sumXY - sumX * sumY) / (n * sumXX - sumX * sumX)
29
30 ' Hesabu ufanisi
31 qPCR_Efficiency = (10 ^ (-1 / slope) - 1) * 100
32End Function
33
1# Kazi ya R kuhesabu ufanisi wa qPCR kutoka kwa thamani za Ct na kipengele cha upunguzaji
2calculate_qpcr_efficiency <- function(ct_values, dilution_factor) {
3 # Unda thamani za upunguzaji wa log
4 log_dilutions <- log10(dilution_factor) * seq(0, length(ct_values) - 1)
5
6 # Fanya uchambuzi wa mwelekeo
7 model <- lm(ct_values ~ log_dilutions)
8
9 # Pata mteremko na R-squared
10 slope <- coef(model)[2]
11 r_squared <- summary(model)$r.squared
12
13 # Hesabu ufanisi
14 efficiency <- (10^(-1/slope) - 1) * 100
15
16 # Rudisha matokeo
17 return(list(
18 efficiency = efficiency,
19 slope = slope,
20 r_squared = r_squared,
21 intercept = coef(model)[1]
22 ))
23}
24
25# Mfano wa matumizi
26ct_values <- c(15.0, 18.5, 22.0, 25.5, 29.0)
27dilution_factor <- 10
28results <- calculate_qpcr_efficiency(ct_values, dilution_factor)
29cat(sprintf("Ufanisi: %.2f%%\n", results$efficiency))
30cat(sprintf("Mteremko: %.4f\n", results$slope))
31cat(sprintf("R-squared: %.4f\n", results$r_squared))
32
1import numpy as np
2from scipy import stats
3import matplotlib.pyplot as plt
4
5def calculate_qpcr_efficiency(ct_values, dilution_factor):
6 """
7 Hesabu ufanisi wa qPCR kutoka kwa thamani za Ct na kipengele cha upunguzaji.
8
9 Parameters:
10 ct_values (list): Orodha ya thamani za Ct
11 dilution_factor (float): Kipengele cha upunguzaji kati ya sampuli zinazofuatana
12
13 Returns:
14 dict: Kichwa chenye ufanisi, mteremko, r_squared, na intercept
15 """
16 # Unda thamani za upunguzaji wa log
17 log_dilutions = np.log10(dilution_factor) * np.arange(len(ct_values))
18
19 # Fanya uchambuzi wa mwelekeo
20 slope, intercept, r_value, p_value, std_err = stats.linregress(log_dilutions, ct_values)
21
22 # Hesabu ufanisi
23 efficiency = (10 ** (-1 / slope) - 1) * 100
24 r_squared = r_value ** 2
25
26 return {
27 'efficiency': efficiency,
28 'slope': slope,
29 'r_squared': r_squared,
30 'intercept': intercept
31 }
32
33def plot_standard_curve(ct_values, dilution_factor, results):
34 """
35 Picha curve ya kiwango pamoja na mwelekeo wa mwelekeo.
36 """
37 log_dilutions = np.log10(dilution_factor) * np.arange(len(ct_values))
38
39 plt.figure(figsize=(10, 6))
40 plt.scatter(log_dilutions, ct_values, color='blue', s=50)
41
42 # Unda pointi za mwelekeo wa mwelekeo
43 x_line = np.linspace(min(log_dilutions) - 0.5, max(log_dilutions) + 0.5, 100)
44 y_line = results['slope'] * x_line + results['intercept']
45 plt.plot(x_line, y_line, 'r-', linewidth=2)
46
47 plt.xlabel('Log Upunguzaji')
48 plt.ylabel('Thamani ya Ct')
49 plt.title('Curve ya Kiwango ya qPCR')
50
51 # Ongeza mwelekeo na R² kwenye picha
52 equation = f"y = {results['slope']:.4f}x + {results['intercept']:.4f}"
53 r_squared = f"R² = {results['r_squared']:.4f}"
54 efficiency = f"Ufanisi = {results['efficiency']:.2f}%"
55
56 plt.annotate(equation, xy=(0.05, 0.95), xycoords='axes fraction')
57 plt.annotate(r_squared, xy=(0.05, 0.90), xycoords='axes fraction')
58 plt.annotate(efficiency, xy=(0.05, 0.85), xycoords='axes fraction')
59
60 plt.grid(True, linestyle='--', alpha=0.7)
61 plt.tight_layout()
62 plt.show()
63
64# Mfano wa matumizi
65ct_values = [15.0, 18.5, 22.0, 25.5, 29.0]
66dilution_factor = 10
67results = calculate_qpcr_efficiency(ct_values, dilution_factor)
68
69print(f"Ufanisi: {results['efficiency']:.2f}%")
70print(f"Mteremko: {results['slope']:.4f}")
71print(f"R-squared: {results['r_squared']:.4f}")
72print(f"Intercept: {results['intercept']:.4f}")
73
74# Picha curve ya kiwango
75plot_standard_curve(ct_values, dilution_factor, results)
76
1/**
2 * Hesabu ufanisi wa qPCR kutoka kwa thamani za Ct na kipengele cha upunguzaji
3 * @param {Array<number>} ctValues - Orodha ya thamani za Ct
4 * @param {number} dilutionFactor - Kipengele cha upunguzaji kati ya sampuli zinazofuatana
5 * @returns {Object} Kichwa chenye ufanisi, mteremko, rSquared, na intercept
6 */
7function calculateQPCREfficiency(ctValues, dilutionFactor) {
8 // Unda thamani za upunguzaji wa log
9 const logDilutions = ctValues.map((_, index) => index * Math.log10(dilutionFactor));
10
11 // Hesabu maana za mwelekeo wa moja kwa moja
12 const n = ctValues.length;
13 let sumX = 0, sumY = 0, sumXY = 0, sumXX = 0, sumYY = 0;
14
15 for (let i = 0; i < n; i++) {
16 sumX += logDilutions[i];
17 sumY += ctValues[i];
18 sumXY += logDilutions[i] * ctValues[i];
19 sumXX += logDilutions[i] * logDilutions[i];
20 sumYY += ctValues[i] * ctValues[i];
21 }
22
23 // Hesabu mteremko na intercept
24 const slope = (n * sumXY - sumX * sumY) / (n * sumXX - sumX * sumX);
25 const intercept = (sumY - slope * sumX) / n;
26
27 // Hesabu R-squared
28 const yMean = sumY / n;
29 let totalVariation = 0;
30 let explainedVariation = 0;
31
32 for (let i = 0; i < n; i++) {
33 const yPredicted = slope * logDilutions[i] + intercept;
34 totalVariation += Math.pow(ctValues[i] - yMean, 2);
35 explainedVariation += Math.pow(yPredicted - yMean, 2);
36 }
37
38 const rSquared = explainedVariation / totalVariation;
39
40 // Hesabu ufanisi
41 const efficiency = (Math.pow(10, -1 / slope) - 1) * 100;
42
43 return {
44 efficiency,
45 slope,
46 rSquared,
47 intercept
48 };
49}
50
51// Mfano wa matumizi
52const ctValues = [15.0, 18.5, 22.0, 25.5, 29.0];
53const dilutionFactor = 10;
54const results = calculateQPCREfficiency(ctValues, dilutionFactor);
55
56console.log(`Ufanisi: ${results.efficiency.toFixed(2)}%`);
57console.log(`Mteremko: ${results.slope.toFixed(4)}`);
58console.log(`R-squared: ${results.rSquared.toFixed(4)}`);
59console.log(`Intercept: ${results.intercept.toFixed(4)}`);
60
Ufanisi mzuri wa qPCR kwa kawaida uko kati ya 90% na 110% (0.9-1.1). Ufanisi wa 100% unawakilisha kuongezeka kwa bidhaa za PCR kwa ukamilifu na kila mzunguko. Ufanisi nje ya upeo huu unaweza kuashiria matatizo na muundo wa primer, hali za mchakato, au uwepo wa vizuizi.
Ufanisi zaidi ya 100% unaweza kutokea kutokana na:
Thamani ya R² ya chini (chini ya 0.98) inaonyesha ukosefu wa usawa katika curve yako ya kiwango, ambayo inaweza kusababishwa na:
Kwa hesabu za ufanisi zinazoweza kuaminika, angalau pointi 3 za upunguzaji zinahitajika, lakini pointi 5-6 zinapendekezwa kwa matokeo sahihi zaidi. Pointi hizi zinapaswa kufunika upeo wote wa nguvu za mkusanyiko zinazotarajiwa katika sampuli zako za majaribio.
Katika upimaji wa jamaa kwa kutumia njia ya ΔΔCt, ufanisi sawa kati ya jeni lengwa na jeni za rejea unatarajiwa (kwa kawaida 100%). Wakati ufanisi unapotofautiana kwa kiasi kikubwa:
Hapana, ufanisi unapaswa kuamuliwa kwa kila jozi ya primer na unapaswa kuthibitishwa tena:
Vizuizi vya PCR vinaweza:
Maneno haya mara nyingi hutumika kwa kubadilishana, lakini:
Kulinganisha sampuli zenye ufanisi tofauti kwa kiasi kikubwa hakupendekezwi kwa sababu:
Bustin SA, Benes V, Garson JA, et al. Mwongozo wa MIQE: taarifa ya chini kwa ajili ya kuchapisha majaribio ya quantitative real-time PCR. Clin Chem. 2009;55(4):611-622. doi:10.1373/clinchem.2008.112797
Pfaffl MW. Mfano mpya wa kihesabu kwa upimaji wa jamaa katika real-time RT-PCR. Nucleic Acids Res. 2001;29(9):e45. doi:10.1093/nar/29.9.e45
Svec D, Tichopad A, Novosadova V, Pfaffl MW, Kubista M. Ufanisi wa makadirio: kuunganisha msingi na upotoshaji katika uchambuzi wa data za quantitative PCR. Biomol Detect Quantif. 2015;3:9-16. doi:10.1016/j.bdq.2015.01.005
Taylor SC, Nadeau K, Abbasi M, Lachance C, Nguyen M, Fenrich J. Jaribio la mwisho la qPCR: Kuzalisha data ya ubora wa kuchapishwa, inayoweza kurudiwa mara ya kwanza. Trends Biotechnol. 2019;37(7):761-774. doi:10.1016/j.tibtech.2018.12.002
Ruijter JM, Ramakers C, Hoogaars WM, et al. Ufanisi wa kuimarisha: kuunganisha msingi na upotoshaji katika uchambuzi wa data za quantitative PCR. Nucleic Acids Res. 2009;37(6):e45. doi:10.1093/nar/gkp045
Higuchi R, Fockler C, Dollinger G, Watson R. Uchambuzi wa PCR wa kinetic: kufuatilia wakati halisi wa majibu ya kuimarisha ya DNA. Biotechnology (N Y). 1993;11(9):1026-1030. doi:10.1038/nbt0993-1026
Bio-Rad Laboratories. Mwongozo wa Maombi ya Real-Time PCR. https://www.bio-rad.com/webroot/web/pdf/lsr/literature/Bulletin_5279.pdf
Thermo Fisher Scientific. Mwongozo wa Real-Time PCR. https://www.thermofisher.com/content/dam/LifeTech/global/Forms/PDF/real-time-pcr-handbook.pdf
Kihesabu chetu cha Ufanisi wa qPCR kinatoa zana rahisi lakini yenye nguvu kwa watafiti kuthibitisha na kuboresha majaribio yao ya quantitative PCR. Kwa kuhesabu ufanisi kwa usahihi kutoka kwa curve za kiwango, unaweza kuhakikisha upimaji wa kuaminika, kutatua matatizo ya majaribio, na kuzingatia mbinu bora katika majaribio ya qPCR.
Jaribu kihesabu chetu leo ili kuboresha ubora na uaminifu wa data yako ya qPCR!
Gundua zana zaidi ambazo zinaweza kuwa na manufaa kwa mtiririko wako wa kazi