Calculateur d'Efficacité qPCR : Outil d'Analyse de Courbe Standard

Calculateur d’efficacité de la qPCR

Un calculateur d’efficacité de la qPCR estime dans quelle mesure une réaction de PCR quantitative (qPCR) copie la séquence d’ADN cible à chaque cycle d’amplification. Il utilise les valeurs de Ct (le numéro du cycle auquel le signal fluorescent franchit un seuil défini) d’une série de dilutions pour calculer l’efficacité, la pente, l’ordonnée à l’origine et le R².

Qu’est-ce que l’efficacité de la qPCR ?

La PCR quantitative suit la fluorescence au fur et à mesure qu’une réaction copie une séquence d’ADN cible, cycle après cycle. L’efficacité de l’amplification décrit l’augmentation de la quantité d’ADN cible à chaque cycle. Avec une efficacité de 100 %, la quantité d’ADN cible double à chaque cycle.

Les réactions réelles atteignent rarement exactement 100 %. La plupart des laboratoires considèrent qu’une efficacité comprise entre 90 % et 110 % est acceptable pour une quantification fiable. Une valeur inférieure à 90 % ou supérieure à 110 % indique un problème, par exemple une mauvaise conception des amorces, des inhibiteurs de PCR entraînés lors de l’extraction de l’échantillon ou des erreurs de pipetage dans la série de dilutions.

L’efficacité est particulièrement importante pour comparer l’expression génique entre les échantillons. La méthode courante ΔΔCt de calcul de la variation relative suppose que chaque réaction a une efficacité de 100 %. Lorsque l’efficacité mesurée diffère de 100 %, ou lorsque l’efficacité d’un gène cible et celle d’un gène de référence diffèrent, les résultats de variation relative peuvent être inexacts.

Formule de l’efficacité de la qPCR

L’efficacité se calcule à partir de la pente de la courbe étalon, une droite obtenue en représentant les valeurs de Ct en fonction du logarithme en base 10 (log10) de la quantité d’ADN initiale dans une série de dilutions.

E = 10^(-1/pente) - 1

E est l’efficacité, exprimée sous forme décimale (multiplier par 100 pour obtenir un pourcentage). La pente est celle de la droite de la courbe étalon. Elle est négative, car le Ct augmente lorsque la quantité d’ADN initiale diminue.

Une pente de -3,32 correspond à une efficacité de 100 % :

10^(-1/-3,32) - 1 = 10^0,301–1 ≈ 1,00

Cela s’explique par le fait qu’une dilution par un facteur 10 devrait décaler le Ct d’environ 3,32 cycles si la réaction doublait la quantité d’ADN à chaque cycle, puisque log2(10) est approximativement égal à 3,32.

En pratique, une pente comprise entre -3,1 et -3,6 est considérée comme acceptable. Cette plage correspond à une efficacité d’environ 90 % à 110 %. Une pente plus proche de -3,1 donne une efficacité calculée plus élevée ; une pente plus proche de -3,6 donne une efficacité plus faible.

Comment calculer l’efficacité de la qPCR : méthode de la courbe étalon

Préparer une série de dilutions. Commencer par l’échantillon le plus concentré et le diluer selon un facteur fixe, généralement 10, pour chaque échantillon suivant. Effectuer la qPCR sur chaque dilution et noter sa valeur de Ct.
Convertir chaque dilution en valeur log10. En comptant le premier échantillon, le plus concentré, comme position 0, l’échantillon numéro i a une valeur log10(quantité relative) égale à -i × log10(facteur de dilution).
Représenter le Ct (axe vertical) en fonction de la valeur log10 (axe horizontal) pour chaque dilution, puis effectuer une régression linéaire.
Relever la pente, l’ordonnée à l’origine et le R² de la droite de régression.
Introduire la pente dans E = 10^(-1/pente) - 1 pour obtenir l’efficacité.

Critères de qualité de la courbe étalon

R² d’au moins 0,98. Cette valeur mesure à quel point les points s’ajustent à une droite. Une valeur plus faible indique souvent des erreurs de pipetage ou une amplification irrégulière.
Pente comprise entre -3,1 et -3,6. Cette plage correspond à une efficacité de 90–110%.
Au moins trois points de dilution. Cinq ou six points donnent une droite plus fiable et facilitent la détection d’une valeur aberrante.
L’échantillon le plus concentré devrait donner un Ct d’environ 15-20. L’échantillon le plus dilué devrait rester inférieur à un Ct de 35, car les valeurs supérieures s’approchent de la limite de détection.

Exemple détaillé

Une série de dilutions en cinq points, au 10, donne les valeurs de Ct suivantes, de l’échantillon le plus concentré au moins concentré : 15,0, 18,5, 22,0, 25,5, 29,0.

Les valeurs log10(quantité relative) sont 0, -1, -2, -3, -4, car chaque étape de dilution ajoute -log10(10), soit -1. Le Ct augmente de 3,5 cycles à chaque étape ; les points se trouvent donc sur une droite et la régression linéaire donne :

Pente : -3,5
Ordonnée à l’origine : 15,0
R² : 1,0

Efficacité : E = 10^(-1/-3,5) - 1 = 10^0,286–1 ≈ 0,931, soit environ 93 %.

Une efficacité de 93 % avec un R² = 1,0 satisfait les critères de qualité habituels et convient à la quantification de l’expression génique.

Exemple de code

Cette fonction JavaScript suit la même méthode que le calculateur. L’indice de chaque valeur de Ct est pris avec son signe opposé avant le calcul du logarithme, car le premier échantillon, le plus concentré, possède la quantité relative la plus élevée et chaque échantillon suivant est plus dilué.

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

Pourquoi l’efficacité est importante pour l’expression génique

Les calculs de variation relative, y compris la méthode ΔΔCt, supposent que chaque réaction double sa cible à chaque cycle. Lorsque les efficacités mesurées d’un gène cible et d’un gène de référence diffèrent, la formule ΔΔCt simple introduit une erreur. Les chercheurs utilisent souvent un calcul corrigé de l’efficacité, comme la méthode décrite par Michael Pfaffl en 2001, lorsque les efficacités diffèrent de plus d’environ cinq points de pourcentage. Les recommandations MIQE, un ensemble de normes de présentation largement cité pour la qPCR, publiées en 2009, préconisent de communiquer l’efficacité calculée, la pente et le R² avec les résultats de quantification.

Foire aux questions

Quel est le bon pourcentage d’efficacité de la qPCR ?

Une valeur comprise entre 90 % et 110 % est généralement acceptée ; les essais bien optimisés se situent autour de 95–105%. Une pente de -3,32 correspond exactement à 100 %.

Pourquoi l’efficacité de ma qPCR est-elle supérieure à 110 % ?

Les causes courantes comprennent des erreurs de pipetage dans la série de dilutions, la formation de dimères d’amorces, des inhibiteurs qui affectent de manière inégale les échantillons concentrés et dilués, ou un seuil de fluorescence défini en dehors de la phase exponentielle de l’amplification. Une courbe de fusion présentant plus d’un pic peut indiquer la présence de dimères d’amorces.

Pourquoi l’efficacité de ma qPCR est-elle inférieure à 90 % ?

Les causes probables comprennent une mauvaise conception des amorces, des inhibiteurs de PCR entraînés lors de l’extraction de l’échantillon ou une température d’hybridation sous-optimale. Vérifier la spécificité des amorces et effectuer un gradient de température peut aider à déterminer la cause.

Que signifie une faible valeur de R² ?

Le R² mesure dans quelle mesure les valeurs de Ct suivent une droite en fonction du logarithme de la quantité. Une valeur inférieure à 0,98 suggère une irrégularité du pipetage, un échantillon inhibé ou un point de dilution proche de la limite de détection.

De combien de points de dilution une courbe étalon a-t-elle besoin ?

Au moins trois, ce qui est le minimum nécessaire pour une régression linéaire. Cinq ou six points donnent une pente plus fiable et facilitent la détection des valeurs aberrantes.

L’efficacité peut-elle varier d’une expérience à l’autre ?

Oui. Un nouveau lot de réactifs, un nouveau lot d’amorces, un type d’échantillon différent ou un autre instrument de qPCR peuvent modifier l’efficacité. De nombreux laboratoires vérifient à nouveau l’efficacité tous les quelques mois ou avec chaque nouveau lot de réactifs.

Références

Bustin SA, Benes V, Garson JA, et al. Les recommandations MIQE : informations minimales pour la publication d’expériences de PCR quantitative en temps réel. Clinical Chemistry. 2009;55(4):611-622.
Pfaffl MW. Un nouveau modèle mathématique pour la quantification relative en RT-PCR en temps réel. Nucleic Acids Research. 2001;29(9):e45.
Svec D, Tichopad A, Novosadova V, Pfaffl MW, Kubista M. Quelle est la fiabilité d’une estimation de l’efficacité de la PCR : recommandations pour une évaluation précise et fiable de l’efficacité de la qPCR. Biomolecular Detection and Quantification. 2015;3:9-16.

Calculateur d'Efficacité qPCR : Outil d'Analyse de Courbe Standard