qPCR 효율성 계산기: 표준 곡선 분석 도구

qPCR 효율 계산기

qPCR 효율 계산기는 정량 PCR(qPCR) 반응이 각 증폭 사이클마다 표적 DNA 서열을 얼마나 잘 복제하는지 추정합니다. 희석 계열에서 형광 신호가 설정된 임계값을 넘는 사이클 번호인 Ct 값을 사용하여 효율, 기울기, 절편 및 R²를 계산합니다.

qPCR 효율이란?

정량 PCR은 반응이 사이클마다 표적 DNA 서열을 복제할 때 형광을 추적합니다. 증폭 효율은 각 사이클마다 표적 DNA의 양이 얼마나 증가하는지를 나타냅니다. 효율이 100%이면 표적 DNA의 양은 매 사이클 두 배가 됩니다.

실제 반응이 정확히 100%에 도달하는 경우는 드뭅니다. 대부분의 실험실에서는 신뢰할 수 있는 정량을 위해 90%에서 110% 사이의 효율을 허용 범위로 봅니다. 90% 미만이나 110% 초과의 값은 부적절한 프라이머 설계, 시료 추출 과정에서 잔류한 PCR 저해제 또는 희석 계열의 피펫팅 오류와 같은 문제를 나타냅니다.

효율은 시료 간 유전자 발현을 비교할 때 특히 중요합니다. 배수 변화를 계산하는 일반적인 ΔΔCt 방법은 모든 반응의 효율이 100%라고 가정합니다. 측정된 효율이 100%와 다르거나 표적 유전자와 참조 유전자의 효율이 서로 다르면 배수 변화 결과가 부정확할 수 있습니다.

qPCR 효율 공식

효율은 표준 곡선의 기울기로 계산합니다. 표준 곡선은 희석 계열에서 시작 DNA 양의 밑이 10인 로그(log10)에 대해 Ct 값을 표시하여 얻은 직선입니다.

E = 10^(-1/기울기) - 1

E는 효율이며 소수로 표시합니다(백분율로 나타내려면 100을 곱합니다). 기울기는 표준 곡선 직선의 기울기입니다. 시작 DNA 양이 감소할수록 Ct가 증가하므로 기울기는 음수입니다.

기울기가 -3.32이면 효율은 100%에 해당하며, 계산은 다음과 같습니다.

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

반응이 매 사이클 DNA 양을 두 배로 늘린다면 10배 희석은 Ct를 약 3.32사이클만큼 이동시켜야 하므로 이 결과가 나옵니다. log2(10)이 약 3.32이기 때문입니다.

실제로는 -3.1에서 -3.6 사이의 기울기를 허용 범위로 봅니다. 이 범위는 대략 90%에서 110%의 효율에 해당합니다. 기울기가 -3.1에 가까울수록 계산된 효율은 높아지고, -3.6에 가까울수록 낮아집니다.

qPCR 효율 계산 방법: 표준 곡선법

희석 계열을 준비합니다. 농도가 가장 높은 시료에서 시작하여 이후 각 시료를 일반적으로 10배와 같이 일정한 배수로 희석합니다. 각 희석액에 qPCR을 수행하고 Ct 값을 기록합니다.
각 희석액을 log10 값으로 변환합니다. 농도가 가장 높은 첫 번째 시료를 위치 0으로 세면, 시료 번호 i의 log10(상대량)은 -i × log10(희석 배수)입니다.
모든 희석액에 대해 Ct를 세로축에, log10 값을 가로축에 표시하고 선형 회귀를 수행합니다.
회귀 직선의 기울기, 절편 및 R²를 읽습니다.
기울기를 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%입니다.

R² = 1.0이고 효율이 93%이면 표준 품질 기준을 충족하므로 유전자 발현 정량에 적합합니다.

코드 예시

이 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 공식에 오차가 생깁니다. 효율이 약 5퍼센트포인트 이상 다를 때 연구자들은 흔히 2001년에 Michael Pfaffl이 설명한 방법과 같은 효율 보정 계산으로 전환합니다. 2009년에 발표된, 널리 인용되는 qPCR 보고 기준인 MIQE 지침은 정량 결과와 함께 계산된 효율, 기울기 및 R²를 보고하도록 권고합니다.

자주 묻는 질문

좋은 qPCR 효율 백분율은 얼마인가?

일반적으로 90%에서 110% 사이를 허용하며, 잘 최적화된 분석법은 약 95–105%에 해당합니다. -3.32의 기울기는 정확히 100%에 해당합니다.

qPCR 효율이 110%보다 높은 이유는 무엇인가?

흔한 원인으로는 희석 계열의 피펫팅 오류, 프라이머-다이머, 농축 시료와 희석 시료에 불균등한 영향을 주는 저해제, 또는 증폭의 지수 증폭 단계 밖에 설정된 형광 임계값이 있습니다. 용융 곡선에서 피크가 두 개 이상 나타나면 프라이머-다이머를 의심할 수 있습니다.

qPCR 효율이 90%보다 낮은 이유는 무엇인가?

가능한 원인으로는 부적절한 프라이머 설계, 시료 추출 과정에서 잔류한 PCR 저해제 또는 최적이 아닌 결합 온도가 있습니다. 프라이머 특이성을 확인하고 온도 구배를 수행하면 원인을 찾는 데 도움이 될 수 있습니다.

R² 값이 낮다는 것은 무엇을 의미하는가?

R²는 Ct 값이 로그량에 대해 직선에 얼마나 잘 맞는지를 나타냅니다. 0.98 미만의 값은 피펫팅의 불일치, 저해된 시료 또는 검출 한계에 가까운 희석 지점을 의미할 수 있습니다.

표준 곡선에는 희석 지점이 몇 개 필요한가?

선형 회귀에 필요한 최소 개수는 세 개입니다. 다섯 개 또는 여섯 개의 지점을 사용하면 기울기를 더 신뢰성 있게 구할 수 있고 이상치를 쉽게 확인할 수 있습니다.

실험마다 효율이 달라질 수 있는가?

그렇습니다. 새로운 시약 로트, 새로운 프라이머 배치, 다른 시료 유형 또는 다른 qPCR 장비로 인해 효율이 달라질 수 있습니다. 많은 실험실에서는 몇 개월마다 또는 새로운 시약 로트를 사용할 때마다 효율을 다시 확인합니다.

참고 문헌

Bustin SA, Benes V, Garson JA 외. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clinical Chemistry. 2009;55(4):611-622.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Research. 2001;29(9):e45.
Svec D, Tichopad A, Novosadova V, Pfaffl MW, Kubista M. PCR 효율 추정치는 얼마나 정확한가: 정밀하고 신뢰할 수 있는 qPCR 효율 평가를 위한 권고. Biomolecular Detection and Quantification. 2015;3:9-16.

qPCR 효율성 계산기: 표준 곡선 분석 도구