DE102020202360B4 - Verfahren und Vorrichtung zur Durchführung eines qPCR-Verfahrens - Google Patents

Verfahren und Vorrichtung zur Durchführung eines qPCR-Verfahrens Download PDF

Info

Publication number
DE102020202360B4
DE102020202360B4 DE102020202360.3A DE102020202360A DE102020202360B4 DE 102020202360 B4 DE102020202360 B4 DE 102020202360B4 DE 102020202360 A DE102020202360 A DE 102020202360A DE 102020202360 B4 DE102020202360 B4 DE 102020202360B4
Authority
DE
Germany
Prior art keywords
qpcr
curve
detected
dna strand
strand section
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
DE102020202360.3A
Other languages
English (en)
Other versions
DE102020202360A1 (de
Inventor
Volker Fischer
Christoph Faigle
Torsten Sachse
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Robert Bosch GmbH
Original Assignee
Robert Bosch GmbH
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Robert Bosch GmbH filed Critical Robert Bosch GmbH
Priority to DE102020202360.3A priority Critical patent/DE102020202360B4/de
Priority to CN202180016607.7A priority patent/CN115103915A/zh
Priority to PCT/EP2021/053656 priority patent/WO2021170441A1/de
Priority to EP21706218.1A priority patent/EP4110947A1/de
Priority to US17/904,301 priority patent/US20230094386A1/en
Publication of DE102020202360A1 publication Critical patent/DE102020202360A1/de
Application granted granted Critical
Publication of DE102020202360B4 publication Critical patent/DE102020202360B4/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6844Nucleic acid amplification reactions
    • C12Q1/6851Quantitative amplification
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/20ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the management or administration of healthcare resources or facilities, e.g. managing hospital staff or surgery rooms
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • G16B25/20Polymerase chain reaction [PCR]; Primer or probe design; Probe optimisation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • G16B40/10Signal processing, e.g. from mass spectrometry [MS] or from PCR
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • G16B40/20Supervised data analysis

Abstract

Computer-implementiertes Verfahren zum Betreiben eines quantitativen Polymerase-Kettenreaktions (qPCR)-Verfahren, mit folgenden Schritten:- Zyklisches Ausführen (S1-S3) von qPCR-Zyklen;- Messen (S11) eines Intensitätswertes einer Fluoreszenz zu jedem qPCR-Zyklus, um eine qPCR-Kurve aus Intensitätswerten zu erhalten;- Auswerten (S12) des Verlaufs der qPCR-Kurve mithilfe eines datenbasierten Klassifikationsmodells, das trainiert ist, um ein Klassifikationsergebnis abhängig von dem Verlauf der qPCR-Kurve bereitzustellen (S13);- Betreiben (S14) des qPCR-Verfahrens abhängig von dem Klassifikationsergebnis der Auswertung des Verlaufs der qPCR-Kurve.

Description

  • Technisches Gebiet
  • Die Erfindung betrifft die Anwendung eines Polymerase-Kettenreaktionsverfahrens (PCR-Verfahrens), insbesondere zur Detektion des Vorliegens eines Erregers. Weiterhin betrifft die vorliegende Erfindung die Auswertung von qPCR-Messungen.
  • Technischer Hintergrund
  • Zur Detektion von DNA-Strangabschnitten in einer zu untersuchenden Substanz, wie beispielsweise einem Serum oder dergleichen, werden in automatisierten Systemen PCR-Verfahren durchgeführt. Die PCR-Systeme ermöglichen es, bestimmte zu detektierende DNA-Strangabschnitte, die z.B. einem Erreger zuzuordnen sind, zu amplifizieren und zu detektieren. Ein PCR-Verfahren umfasst generell eine zyklische Anwendung der Schritte der Denaturierung, des Annealing und der Elongation. Insbesondere wird beim PCR-Prozess ein DNA-Doppelstrang in Einzelstränge aufgespalten und diese jeweils durch Anlagerung von Nukleotiden wieder vervollständigt, um die DNA-Strangabschnitte in jedem Zyklus zu vervielfältigen.
  • Das qPCR-Verfahren ermöglicht die Erregerlast, nachgewiesen mit diesem Prozess zu quantifizieren. Dazu sind die Nukleotide zumindest teilweise mit Fluoreszenzmolekülen versehen, die bei Anbinden an den Einzelstrang des zu detektierenden DNA-Strangabschnitts eine Fluoreszenzeigenschaft aktivieren.
  • Nach dem Aufbau der Doppelstränge kann nach jedem Zyklus ein Fluoreszenzwert ermittelt werden, der von der Anzahl der erzeugten DNA-Strangabschnitte abhängt.
  • Im Laufe der Amplifikation kann dann aus den ermittelten Fluoreszenzwerten (Intensitätswerten) eine qPCR-Kurve ermittelt werden, die bei Vorliegen des zu detektierenden DNA-Strangabschnitts in der zu untersuchenden Substanz einen sigmoidförmigen Verlauf hat. Die gemessenen qPCR-Kurven können in der Realität mit Artefakten behaftet sein, so dass in der Regel mehrere parallele Messungen durchgeführt werden, um durch eine Mittelwertbildung der Messwerte eine genauere Auswertung der qPCR-Kurven zu ermöglichen.
  • Offenbarung der Erfindung
  • Erfindungsgemäß sind ein Verfahren zum Durchführen eines qPCR-Verfahrens gemäß Anspruch 1 sowie eine Vorrichtung und ein qPCR-System gemäß den nebengeordneten Ansprüchen vorgesehen.
  • Weitere Ausgestaltungen sind in den abhängigen Ansprüchen angegeben.
  • Gemäß einem ersten Aspekt ist ein Verfahren zum Betreiben eines quantitativen Polymerase-Kettenreaktions (qPCR)-Verfahren vorgesehen, mit folgenden Schritten:
    • - Zyklisches Ausführen von qPCR-Zyklen
    • - Messen eines Intensitätswertes einer Fluoreszenz nach oder während jedem qPCR-Zyklus, um eine qPCR-Kurve aus Intensitätswerten zu erhalten;
    • - Auswerten des Verlaufs der qPCR-Kurve mithilfe eines datenbasierten Klassifikationsmodells, das trainiert ist, um ein Klassifikationsergebnis abhängig von dem Verlauf der qPCR-Kurve bereitzustellen;
    • - Betreiben des qPCR-Verfahrens abhängig von dem Klassifikationsergebnis der Auswertung des Verlaufs der qPCR-Kurve.
  • Das qPCR-Verfahren weist eine zyklische Wiederholung der Schritte der Denaturierung, des Annealing und der Elongation auf. Bei der Denaturierung wird bei einer hohen Temperatur die gesamte doppelsträngige DNA in der zu untersuchenden Substanz in zwei Einzelstränge aufgespalten. In dem Schritt des Annealing wird an die Einzelstränge einer der Substanz zugegebenen Primer angebunden, die den Startpunkt der Amplifizierung der zu detektierenden DNA-Strangabschnitte vorgeben. Bei dem Schritt der Elongation wird ein zweiter komplementärer DNA-Strangabschnitt aus freien Nukleotiden auf den mit dem Primer versehenen Einzelsträngen aufgebaut. Nach jedem dieser Zyklen hat sich also idealerweise die DNA-Menge der zu detektierenden DNA-Strangabschnitte verdoppelt.
  • Durch die Anwendung des qPCR-Verfahrens werden in die zu detektierenden DNA-Strangabschnitte Fluoreszenzmoleküle als Marker eingebaut, so dass über eine Messung der Intensität der Fluoreszenz nach jedem Elongationsschritt ein zeitlicher Verlauf der Intensitätswerte ermittelt werden kann. Die so erhaltene qPCR-Kurve weist dabei drei distinkte Phasen auf, nämlich eine Baseline, in der die Intensität der Fluoreszenz des von eingebauten Markern emittierten Fluoreszenzlichts noch nicht von der Hintergrund-Fluoreszenz zu unterscheiden ist, einer exponentiellen Phase, in der die Fluoreszenzintensität über die Baseline ansteigt, also sichtbar wird, wobei durch die Verdoppelung der DNA-Stränge in jedem Zyklus das Fluoreszenzsignal proportional zur Menge der zu detektierenden DNA-Strangabschnitte exponentiell ansteigt, und einer Plateau-Phase, in der die Reagenzien, d. h. der Primer und die freien Nukleotide, nicht mehr in der benötigten Konzentration vorhanden sind und keine weitere Verdoppelung stattfindet.
  • Für den Nachweis eines vorgegebenen zu detektierenden DNA-Strangabschnitts, der beispielsweise einem Erreger entsprechen kann, ist hierbei der sogenannte ct (cycle threshold)-Wert maßgeblich. Der ct-Wert bestimmt den Beginn der exponentiellen Phase und wird durch Überschreiten eines spezifischen Grenzwerts, der für den jeweils zu detektierenden DNA-Strangabschnitt festgelegt ist und der für alle Proben für den zu detektierenden DNA-Strangabschnitt identisch ist, ermittelt oder rechnerisch durch die zweite Ableitung der qPCR-Kurve in der exponentiellen Phase ermittelt und dem Intensitätswert des steilsten Anstiegs der qPCR-Kurve entspricht. Ist der Zielwert bekannt, kann durch Rückrechnung die Anfangskonzentration des zu detektierenden DNA-Strangabschnitts in der zu untersuchenden Substanz bestimmt werden.
  • In der Realität sind die qPCR-Kurven sehr ungenau und erheblichen Schwankungen unterworfen. Es kann ein Baseline-Drift auftreten, mit der das Ansteigen der Hintergrundfluoreszenz über die Messzyklen bezeichnet wird. Das heißt, dass, selbst wenn keine Amplifikation stattfindet, das Fluoreszenzsignal ansteigt. Weitere Einflussfaktoren, die sich negativ auf die Genauigkeit der qPCR-Kurve auswirken, können beispielsweise aus thermischem Rauschen, Schwankungen bzw. Zumesstoleranzen in der Reagenzkonzentration, Lufteinschlüsse und Artefakten im Fluoreszenzvolumen resultieren.
  • In herkömmlichen qPCR-Systemen findet einerseits eine softwarebasierte Korrektur der qPCR-Kurven statt, andererseits kann vorgesehen sein, eine Probe unter gleichen Bedingungen mehrfach zu messen und die resultierenden qPCR-Kurven durch Mittelwertbildung zu glätten. Dies erfordert jedoch einen erhöhten Aufwand.
  • Eine Idee des obigen Verfahrens besteht darin, mithilfe eines datenbasierten Klassifikationsmodells, insbesondere in Form eines künstlichen neuronalen Netzes, insbesondere eines tiefen neuronalen Netzes oder eines rekurrenten neuronalen Netzes, insbesondere eines LSTM, oder in Form einer Support Vector Machine, zu entscheiden, ob eine gemessene qPCR-Kurve angibt, dass der zu detektierende DNA-Strangabschnitt in der zu untersuchenden Substanz vorkommt oder nicht. Da gemessene qPCR-Kurven in der Regel stark rauschbehaftet sind, kann durch die Benutzung des trainierten datenbasierten Klassifikationsmodells der qPCR-Kurvenverlauf auch in Grenzfällen zuverlässiger dahingehend bewertet werden, ob der zu detektierende DNA-Strangabschnitt in der zu untersuchenden Substanz vorkommt oder nicht.
  • Das Klassifikationsmodell kann trainiert sein, um abhängig von der qPCR-Kurve ein Klassifikationsergebnis bereitzustellen, das ein Vorhandensein eines zu detektierenden DNA-Strangabschnitts angibt oder nicht. Insbesondere ist es möglich, das datenbasierte Klassifikationsmodell auf gelabelte qPCR-Kurven als Trainingsdaten zu trainieren, so dass eine gemessene qPCR-Kurve in entsprechender Weise klassifiziert werden kann. Dies erfordert jedoch ein umfangreiches Training des Klassifikationsmodells mit einer hohen Anzahl von Datensätzen aus manuell klassifizierten bzw. gelabelten qPCR-Kurven, da keine Vorannahmen berücksichtigt werden.
  • Weiterhin können Restfehlerverläufe zwischen der gemessenen qPCR-Kurve und einer parametrierten Vorhandenseinsfunktion und einer parametrierten Nichtvorhandenseinsfunktion bestimmt werden, wobei das Klassifikationsmodell trainiert ist, um abhängig von mindestens einer der Restfehlerverläufe ein Klassifikationsergebnis bereitzustellen, das ein Vorhandensein eines zu detektierenden DNA-Strangabschnitts angibt oder nicht. Somit ist es alternativ möglich, die gemessene qPCR-Kurve an parametrisierte Idealverläufe einer qPCR-Kurve bei Vorhandensein eines zu detektierenden DNA-Strangabschnitts (Vorhandenseins-Kurvenverlauf) in der zu untersuchenden Substanz und bei Nichtvorhandensein eines zu detektierenden DNA-Strangabschnitts (Nicht-Vorhandenseins-Kurvenverlauf) in der zu untersuchenden Substanz zu fitten. Dies erfolgt durch Parametrisierung des Vorhandenseins-Kurvenverlauf und des Nicht-Vorhandenseins-Kurvenverlauf an die gemessene qPCR-Kurve. Die Restfehlerverläufe zwischen dem parametrisierten Vorhandenseins-Kurvenverlauf bzw. dem parametrisierten Nicht-Vorhandenseins-Kurvenverlauf und der tatsächlichen gemessenen qPCR-Kurve können dann mithilfe eines trainierten datenbasierten Klassifikationsmodells ausgewertet werden. Das Klassifikationsmodell ist dabei trainiert, zu entscheiden, ob der jeweilige Kurvenfit auf einem parametrisierten Idealverlauf basiert, der dem gemessenen Verlauf der qPCR-Kurve am nächsten kommt. Das heißt, mithilfe des Klassifikationsmodells wird festgestellt, ob die gemessene qPCR-Kurve eher einer VorhandenseinqPCR-Kurve oder Nichtvorhandensein-qPCR-Kurve entspricht.
  • Dabei wird ein Vorhandensein eines zu detektierenden DNA-Strangabschnitts festgestellt, wenn das Klassifikationsergebnis basierend auf dem Restfehlerverlauf aus der parametrierten Vorhandenseinsfunktion das Vorhandensein des zu detektierenden DNA-Strangabschnitts anzeigt. Entsprechend kann ein Nichtvorhandensein des zu detektierenden DNA-Strangabschnitts festgestellt werden, wenn das Klassifikationsergebnis basierend auf dem Restfehlerverlauf aus der parametrierten Nichtvorhandenseinsfunktion das Nichtvorhandensein des zu detektierenden DNA-Strangabschnitts anzeigt.
  • Die Auswertung kann erfolgen, indem auf ein Vorhandensein des zu detektierenden DNA-Strangabschnitt geschlossen wird, wenn das Klassifikationsergebnis einen Kurvenfit mit dem Vorhandenseins-Kurvenverlauf bestätigt und einen Kurvenfit mit dem Nicht-Vorhandenseins-Kurvenverlauf nicht bestätigt. Analog wird auf ein Nicht-Vorhandensein des zu detektierenden DNA-Strangabschnitt geschlossen, wenn das Klassifikationsergebnis einen Kurvenfit mit dem Nicht-Vorhandenseins-Kurvenverlauf bestätigt und einen Kurvenfit mit dem Vorhandenseins-Kurvenverlauf nicht bestätigt.
  • Während die erste oben beschriebene Variante mögliche fehlerhafte oder ungenaue Grundannahmen über die zugrundeliegenden typischen Kurvenverläufe vermeidet und auch unbekannte Zusammenhänge und Dynamiken abbilden kann, kann die zweite oben beschriebene Variante aufgrund des vorhandenen Domänenwissens mit einem Klassifikationsmodell verwendet werden, das für sein Training einen weniger umfangreichen Trainingsdatensatz erfordert. Insgesamt führt die Nutzung eines datenbasierten Klassifikationsmodells zu einer geringeren Anzahl von Fehlklassifikationen im Vergleich zu herkömmlichen Verfahren.
  • Es kann vorgesehen sein, dass das qPCR-Verfahren betrieben wird, indem bei Feststellen eines Vorhandenseins des zu detektierenden DNA-Strangabschnitts,
    • - signalisiert wird, dass ein ct-Wert ermittelbar ist, und/oder
    • - der ct-Wert aus der parametrierten Vorhandenseinsfunktion ermittelt wird.
  • Gemäß einem weiteren Aspekt ist eine Vorrichtung zum Betreiben eines quantitativen Polymerase-Kettenreaktions (qPCR)-Verfahren, wobei die Vorrichtung ausgebildet ist, um folgende Schritte auszuführen:
    • - Zyklisches Ausführen von qPCR-Zyklen
    • - Messen eines Intensitätswertes einer Fluoreszenz zu jedem qPCR-Zyklus, um eine qPCR-Kurve aus Intensitätswerten zu erhalten;
    • - Auswerten des Verlaufs der qPCR-Kurve mithilfe eines datenbasierten Klassifikationsmodells, das trainiert ist, um ein Klassifikationsergebnis abhängig von dem Verlauf der qPCR-Kurve bereitzustellen;
    • - Betreiben des qPCR-Verfahrens abhängig von dem Klassifikationsergebnis der Auswertung des Verlaufs der qPCR-Kurve.
  • Kurzbeschreibung der Zeichnungen
  • Ausführungsformen werden nachfolgend anhand der beigefügten Zeichnungen näher erläutert. Es zeigen:
    • 1 eine schematische Darstellung eines Zyklus eines PCR-Verfahrens;
    • 2 eine schematische Darstellung einer typischen qPCR-Kurve mit einem Intensitätswerteverlauf;
    • 3 ein gemessener Verlauf einer qPCR-Kurve;
    • 4a und 4b ideale Verläufe der qPCR-Kurve für den Fall einer nicht nachweisbaren Substanz bzw. einer nachweisbaren Substanz; und
    • 5 ein Flussdiagramm zur Veranschaulichung eines Verfahrens zum Betreiben einer qPCR-Messung; und
    • 6 ein Flussdiagramm zur Veranschaulichung eines weiteren Verfahrens zum Betreiben einer qPCR-Messung.
  • Beschreibung von Ausführungsformen
  • 1 zeigt eine schematische Darstellung eines an sich bekannten PCR-Verfahrens mit den Schritten der Denaturierung des Annealing und der Elongation.
  • In dem Schritt des Annealing S1 wird bei einer hohen Temperatur von beispielsweise über 90°C die doppelsträngige DNA in einer Substanz in zwei Einzelstränge aufgebrochen. In einem nachfolgenden Annealing-Schritt S2 wird an die Einzelstränge an einer bestimmten DNA-Position, die einen Beginn eines zu detektierenden DNA-Strangabschnitts markiert, ein sogenannter Primer angebunden. Dieser Primer stellt den Startpunkt einer Amplifizierung des DNA-Strangabschnitts dar. In einem Elongationsschritt S3 wird beginnend an der von dem Primer markierten Position aus der Substanz zugesetzten freien Nukleotiden der komplementäre DNA- Strangabschnitts an den Einzelsträngen aufgebaut, so dass am Ende des Elongationsschrittes die zuvor aufgespaltenen Einzelstränge zu vollständigen Doppelsträngen ergänzt worden sind.
  • Durch Versehen der freien Nukleotide bzw. des Primers mit Fluoreszenzmolekülen, die nur in dem an den DNA-Strangabschnitt gebundenen Zustand Fluoreszenzeigenschaften aufweisen, ist es möglich, durch Ermitteln einer Intensität einer Fluoreszenz im Anschluss an den Elongationsschritt S3 einen Intensitätswertes der Fluoreszenz durch eine geeignete Messung zu erhalten. Der gemessenen Intensität des Fluoreszenzlichts wird ein Intensitätswert zugeordnet.
  • Das Verfahren der Schritte S1 bis S3 wird zyklisch ausgeführt und die Intensitätswerte aufgezeichnet, um einen Intensitätswerteverlauf als eine qPCR-Kurve zu erhalten.
  • Der Intensitätswerteverlauf hat im Idealfall einen Verlauf, der in 2 dargestellt ist. 2 zeigt einen Verlauf einer normalisierten Intensität I über dem Zyklusindex z. Dieser Verlauf unterteilt sich in drei Abschnitte, nämlich einen Baseline-Abschnitt B, in dem die Fluoreszenz der eingebauten Fluoreszenzmoleküle noch nicht von einer Hintergrundfluoreszenz zu unterscheiden ist, einen exponentiellen Abschnitt E, in dem die Intensitätswerte sichtbar sind und exponentiell ansteigen, und in einem Plateau-Abschnitt P, in dem sich der Anstieg der Intensitätswerte abflacht, da die verwendeten Reagenzien (Lösung mit Nukleotiden) aufgebraucht sind und keine weitere Anbindung an aufgebrochenen Einzelsträngen stattfindet.
  • In 3 ist beispielhaft ein bei einer realen Messung erhaltener Verlauf der Intensitätswerte als qPCR-Kurve dargestellt. Man erkennt starke Schwankungen, die sich aus einer Hintergrundfluoreszenz, thermisches Rauschen, Schwankungen in den Reagenzkonzentrationen sowie Bläschen und Artefakte im Fluoreszenzvolumen ergeben können. Man erkennt, dass die Bestimmung des Baseline-Abschnitts, des exponentiellen Abschnitts und des Plateau-Abschnitts der qPCR-Kurve nicht ohne weiteres möglich ist.
  • 4a und 4b zeigen ideale Verläufe einer qPCR-Kurve ohne Vorliegen eines nachzuweisenden DNA-Strangabschnitts bzw. mit Vorliegen eines nachzuweisenden DNA-Strangabschnitts.
  • In 5 ist ein Flussdiagramm zur Veranschaulichung eines Verfahrens zum Betreiben eines qPCR-Messprozesses dargestellt. Das Verfahren kann auf einer Datenverarbeitungseinrichtung ausgeführt werden, die einen qPCR-Prozess auf einem PCR-System steuert und die von einem PCR-System in jedem Zyklus einen Intensitätswert, der die Intensität einer Fluoreszenz abhängig von einem während der qPCR amplifizierten DNA-Strangabschnitt angibt, bereitstellt. In der Datenverarbeitungseinrichtung kann das nachfolgend beschriebene Verfahren in Software und/oder Hardware implementiert sein.
  • In Schritt S11 wird die qPCR-Messung durchgeführt, um Intensitätswerte in aufeinanderfolgenden Zyklen einer qPCR-Messung zu empfangen. Die Anzahl der Zyklen für die qPCR-Messung liegt bei etwa 30 bis 60 Zyklen, vorzugsweise bei 40 Zyklen. Man erhält eine qPCR-Kurve, die die Intensitätswerte (oder davon abgeleitete Werte) über einen Zyklusindex abbildet.
  • In Schritt S12 werden die Intensitätswerte der qPCR-Kurve einem trainierten Klassifikationsmodell zugeführt. Das Klassifikationsmodell ist als datenbasiertes Modell ausgebildet, wie z. B. eine SVM (Support Vector Machine) oder ein tiefes neuronales Netz. Alternativ kann das datenbasierte Klassifikationsmodell auch mit einem neuronalen Netz aus zeitlichen Konvolutionsschichten gebildet sein.
  • Das Klassifikationsmodell kann mit Datensätzen aus real gemessenen qPCR-Verläufen trainiert sein, denen jeweils ein Label zugeordnet ist, das angibt, ob der Datensatz (qPCR-Kurve) einer Messung einer Substanz, die den zu detektierenden DNA-Strangabschnitt enthält, entspricht oder nicht.
  • In Schritt S13 wird entsprechend dem Ergebnis der Klassifikation durch das Klassifikationsmodell bestimmt, ob die qPCR-Kurve ein Vorhandensein oder ein Nichtvorhandensein des zu detektierenden DNA-Strangabschnitts entspricht, d. h. angibt, ob der zu detektierenden DNA-Strangabschnitt in der Substanz enthalten ist oder nicht.
  • In Schritt S14 wird das PCR-Verfahren entsprechend dem Klassifikationsergebnis ausgeführt. Insbesondere kann das qPCR-Verfahren betrieben werden, indem bei Feststellen eines Vorhandenseins des zu detektierenden DNA-Strangabschnitts, signalisiert wird, dass ein ct-Wert ermittelbar ist, und der ct-Wert aus der parametrierten Vorhandenseinsfunktion ermittelt.
  • In 6 ist ein Flussdiagramm zur Veranschaulichung eines alternativen Verfahrens zum Betreiben der PCR-Messung und zur Auswertung einer gemessenen qPCR-Kurve dargestellt.
  • In Schritt S21 wird mithilfe eines qPCR-Verfahrens eine qPCR-Messung durchgeführt und eine qPCR-Kurve durch aufeinanderfolgendes Messen von Intensitätswerten bestimmt.
  • In Schritt S22 wird zunächst eine vorgegebene parametrische Nichtvorhandenseinsfunktion parametrisiert, indem die gemessene qPCR-Kurve an die Nichtvorhandenseinsfunktion gefittet wird. Beispielsweise kann die Nichtvorhandenseinsfunktion eine lineare Funktion, wie sie in 4a dargestellt ist, sein, die im Wesentlichen der Verlaufscharakteristik einer Baseline entspricht. Das Fitten an die Nichtvorhandenseinsfunktion kann beispielsweise mithilfe einer Least-Squares-Minimierung erfolgen.
  • In Schritt S23 wird die gemessene qPCR-Kurve an eine Vorhandenseinsfunktion gefittet. Die Vorhandenseinsfunktion entspricht einer parametrisierten Funktion, die im Wesentlichen einer Verlaufscharakteristik, wie sie in 4b dargestellt ist, entspricht. Die Vorhandenseinsfunktion kann eine Sigmoid-Funktion aufweisen. Durch Parametrisieren der Vorhandenseinsfunktion kann die gemessene qPCR-Kurve an die Vorhandenseinsfunktion gefittet werden.
  • In Schritt S24 werden Restfehlerverläufe der gemessenen qPCR-Kurve zu der parametrisierten Vorhandenseinsfunktion und zu der parametrisierten Nichtvorhandenseinsfunktion ermittelt.
  • In Schritt S25 werden die Restfehlerverläufe einem datenbasierten Klassifikationsmodell zugeführt. Das Klassifikationsmodell ist basierend auf Trainingsdaten trainiert, die Restfehlerverläufe den entsprechenden Vorhandenseins- oder Nichtvorhandenseinsfunktion zuordnet. Das heißt, die Trainingsdaten geben an, ob ein Restfehlerverlauf aus der parametrisierten Vorhandenseinsfunktion das Vorhandensein des zu detektierenden Strangabschnitt bestätigt oder nicht. Weiterhin geben die Trainingsdaten an, ob der Restfehlerverlauf aus der parametrisierten Nichtvorhandenseinsfunktion das Nicht-Vorhandensein des zu detektierenden Strangabschnitt bestätigt oder nicht.
  • Entsprechend wird in Schritt S25 mithilfe des Restfehlerverlaufs festgestellt, dass eine gemessene qPCR-Kurve ein Vorhandensein des zu detektierenden DNA-Strangabschnitts in der Substanz angibt, wenn das Klassifikationsmodell den Restfehlerverlauf aus der Vorhandenseinsfunktion bestätigt. Analog wird mithilfe des Restfehlerverlaufs festgestellt, dass eine gemessene qPCR-Kurve ein Nichtvorhandensein des zu detektierenden DNA-Strangabschnitts in der Substanz angibt, wenn das Klassifikationsmodell den Restfehlerverlauf aus der NichtVorhandenseinsfunktion bestätigt. Das heißt, ergibt sich allgemein aus dem Restfehlerverlauf basierend auf der parametrierten Vorhandenseinsfunktion oder der parametrierten Nichtvorhandenseinsfunktion, dass das Klassifikationsergebnis den Restfehlerverlauf der entsprechenden parametrierten Vorhandenseinsfunktion oder der parametrierten Nichtvorhandenseinsfunktion, die dem Restfehlerverlauf zugrunde liegt, so wird festgestellt, dass der zu detektierende DNA-Strangabschnitt vorhanden bzw. nicht vorhanden ist.
  • Ergibt sich aus dem Klassifikationsergebnis bezüglich des Restfehlerverlaufs ein dem Restfehlerverlauf zugrunde liegenden parametrisierten Vorhandenseins- oder Nichtvorhandenseinsfunktion entgegengesetztes Ergebnis, so kann entschieden werden, die qPCR-Messung zu verwerfen, da eine eindeutige Entscheidung über ein Vorhandensein oder Nichtvorhandensein des zu detektierenden DNA-Strangabschnitts nicht getroffen werden kann.

Claims (10)

  1. Computer-implementiertes Verfahren zum Betreiben eines quantitativen Polymerase-Kettenreaktions (qPCR)-Verfahren, mit folgenden Schritten: - Zyklisches Ausführen (S1-S3) von qPCR-Zyklen; - Messen (S11) eines Intensitätswertes einer Fluoreszenz zu jedem qPCR-Zyklus, um eine qPCR-Kurve aus Intensitätswerten zu erhalten; - Auswerten (S12) des Verlaufs der qPCR-Kurve mithilfe eines datenbasierten Klassifikationsmodells, das trainiert ist, um ein Klassifikationsergebnis abhängig von dem Verlauf der qPCR-Kurve bereitzustellen (S13); - Betreiben (S14) des qPCR-Verfahrens abhängig von dem Klassifikationsergebnis der Auswertung des Verlaufs der qPCR-Kurve.
  2. Verfahren nach Anspruch 1, wobei das datenbasierte Klassifikationsmodell ein neuronales Netz, insbesondere ein tiefes neuronales Netz oder ein rekurrentes neuronales Netz, insbesondere ein LSTM, oder ein Support Vector Machine Modell umfasst.
  3. Verfahren nach einem der Ansprüche 1 bis 2, wobei das Klassifikationsmodell trainiert ist, um abhängig von der qPCR-Kurve das Klassifikationsergebnis bereitzustellen, das ein Vorhandensein eines zu detektierenden DNA-Strangabschnitts angibt oder nicht.
  4. Verfahren nach einem der Ansprüche 1 bis 2, wobei Restfehlerverläufe zwischen der gemessenen qPCR-Kurve und einer parametrierten Vorhandenseinsfunktion und einer parametrierten Nichtvorhandenseinsfunktion bestimmt werden (S24), wobei das Klassifikationsmodell trainiert ist, um abhängig von mindestens einer der Restfehlerverläufe ein Klassifikationsergebnis bereitzustellen (S25), das ein Vorhandensein eines zu detektierenden DNA-Strangabschnitts angibt oder nicht.
  5. Verfahren nach Anspruch 4, wobei ein Vorhandensein eines zu detektierenden DNA-Strangabschnitts festgestellt wird, wenn das Klassifikationsergebnis basierend auf dem Restfehlerverlauf aus der parametrierten Vorhandenseinsfunktion das Vorhandensein des zu detektierenden DNA-Strangabschnitts anzeigt.
  6. Verfahren nach Anspruch 4 oder 5, wobei ein Nichtvorhandensein des zu detektierenden DNA-Strangabschnitts festgestellt wird (S25), wenn das Klassifikationsergebnis basierend auf dem Restfehlerverlauf aus der parametrierten Nichtvorhandenseinsfunktion das Nichtvorhandensein des zu detektierenden DNA-Strangabschnitts anzeigt.
  7. Verfahren nach einem der Ansprüche 1 bis 6, wobei das qPCR-Verfahren betrieben wird (S14), indem bei Feststellen eines Vorhandenseins des zu detektierenden DNA-Strangabschnitts, - signalisiert wird, dass ein ct-Wert ermittelbar ist, - der ct-Wert aus der parametrierten Vorhandenseinsfunktion ermittelt wird.
  8. Vorrichtung zum Betreiben eines quantitativen Polymerase-Kettenreaktions (qPCR)-Verfahren, wobei die Vorrichtung ausgebildet ist, um folgende Schritte auszuführen: - Zyklisches Ausführen von qPCR-Zyklen - Messen eines Intensitätswertes einer Fluoreszenz zu jedem qPCR-Zyklus, um eine qPCR-Kurve aus Intensitätswerten zu erhalten; - Auswerten des Verlaufs der qPCR-Kurve mithilfe eines datenbasierten Klassifikationsmodells, das trainiert ist, um ein Klassifikationsergebnis abhängig von dem Verlauf der qPCR-Kurve bereitzustellen; - Betreiben des qPCR-Verfahrens abhängig von dem Klassifikationsergebnis der Auswertung des Verlaufs der qPCR-Kurve.
  9. Computerprogramm, welches dazu eingerichtet ist, alle Schritte eines Verfahrens nach einem der Ansprüche 1 bis 7 auszuführen.
  10. Elektronisches Speichermedium, auf welchem ein Computerprogramm nach Anspruch 9 gespeichert ist.
DE102020202360.3A 2020-02-25 2020-02-25 Verfahren und Vorrichtung zur Durchführung eines qPCR-Verfahrens Active DE102020202360B4 (de)

Priority Applications (5)

Application Number Priority Date Filing Date Title
DE102020202360.3A DE102020202360B4 (de) 2020-02-25 2020-02-25 Verfahren und Vorrichtung zur Durchführung eines qPCR-Verfahrens
CN202180016607.7A CN115103915A (zh) 2020-02-25 2021-02-15 用于执行qPCR方法的方法和设备
PCT/EP2021/053656 WO2021170441A1 (de) 2020-02-25 2021-02-15 Verfahren und vorrichtung zur durchführung eines qpcr-verfahrens
EP21706218.1A EP4110947A1 (de) 2020-02-25 2021-02-15 Verfahren und vorrichtung zur durchführung eines qpcr-verfahrens
US17/904,301 US20230094386A1 (en) 2020-02-25 2021-02-15 Method and Device for Carrying out a qPCR Process

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102020202360.3A DE102020202360B4 (de) 2020-02-25 2020-02-25 Verfahren und Vorrichtung zur Durchführung eines qPCR-Verfahrens

Publications (2)

Publication Number Publication Date
DE102020202360A1 DE102020202360A1 (de) 2021-08-26
DE102020202360B4 true DE102020202360B4 (de) 2024-03-21

Family

ID=74666703

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102020202360.3A Active DE102020202360B4 (de) 2020-02-25 2020-02-25 Verfahren und Vorrichtung zur Durchführung eines qPCR-Verfahrens

Country Status (5)

Country Link
US (1) US20230094386A1 (de)
EP (1) EP4110947A1 (de)
CN (1) CN115103915A (de)
DE (1) DE102020202360B4 (de)
WO (1) WO2021170441A1 (de)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060009916A1 (en) 2004-07-06 2006-01-12 Xitong Li Quantitative PCR data analysis system (QDAS)
WO2014053467A1 (en) 2012-10-02 2014-04-10 Roche Diagnostics Gmbh Universal method to determine real-time pcr cycle threshold values
WO2017025589A1 (de) 2015-08-13 2017-02-16 Cladiac Gmbh Verfahren und testsystem zum nachweis und/oder quantifizieren einer ziel-nukleinsäure in einer probe
DE102018213026A1 (de) 2018-08-03 2020-02-06 Robert Bosch Gmbh Verfahren zur Durchführung einer Echtzeit-PCR

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3411819A4 (de) * 2016-02-05 2019-10-23 Seegene, Inc. Verfahren zur reduzierung des rauschpegels eines datensatzes für einen zielanalyten

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060009916A1 (en) 2004-07-06 2006-01-12 Xitong Li Quantitative PCR data analysis system (QDAS)
WO2014053467A1 (en) 2012-10-02 2014-04-10 Roche Diagnostics Gmbh Universal method to determine real-time pcr cycle threshold values
WO2017025589A1 (de) 2015-08-13 2017-02-16 Cladiac Gmbh Verfahren und testsystem zum nachweis und/oder quantifizieren einer ziel-nukleinsäure in einer probe
DE102018213026A1 (de) 2018-08-03 2020-02-06 Robert Bosch Gmbh Verfahren zur Durchführung einer Echtzeit-PCR

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GÜNAY, M. [u.a.]: Machine Learning for Optimum CT-Prediction for qPCR. 15th IEEE International Conference on Machine Learning and Applications (ICMLA) (2016): 588-592
TICHOPAD, A. [u.a.]: Quality control for quantitative PCR based on amplification compatibility test. Methods (2010) 50 (4) 308-12

Also Published As

Publication number Publication date
US20230094386A1 (en) 2023-03-30
DE102020202360A1 (de) 2021-08-26
WO2021170441A1 (de) 2021-09-02
CN115103915A (zh) 2022-09-23
EP4110947A1 (de) 2023-01-04

Similar Documents

Publication Publication Date Title
DE102004046388A1 (de) Quantifizierung von Nukleinsäuren unter Verwendung von Wachstumskurven
DE102016202428B4 (de) Messsystem für colorimetrische Assays
DE102018122961A1 (de) Chromatographische datensystemverarbeitungseinrichtung
EP3159681B1 (de) Verfahren und vorrichtung zur automatisierbaren ermittlung der bestimmungsgrenze und des relativen fehlers bei der quantifizierung der konzentration einer zu untersuchenden substanz in einer messprobe
EP3329284A1 (de) Verfahren zur bestimmung von fibrinogen
DE102020202360B4 (de) Verfahren und Vorrichtung zur Durchführung eines qPCR-Verfahrens
EP3839958B1 (de) Verfahren und vorrichtung zur qualitativen auswertung von real-time pcr daten
DE60013303T2 (de) Verfahren zum entfernen systematischer fehler und abweichungen und zum abschätzen zufälliger fehler in chemischen und biologischen testverfahren
DE4331018C2 (de) Verfahren zur Bewertung von Blutproben
EP2288916B1 (de) Immunochromatographisches verfahren und testsystem zur bestimmung von wenigstens einem analyten in einer zu untersuchenden testlösung
DE102020202358B4 (de) Verfahren und Vorrichtung zur Durchführung eines qPCR-Verfahrens
DE102016217948A1 (de) Verfahren zur Prognose der Qualität von Klebverbindungen
DE102020202361A1 (de) Verfahren und Vorrichtung zur Auswertung einer qPCR-Kurve
DE102020202363A1 (de) Verfahren und Vorrichtung zur Durchführung eines qPCR-Verfahrens
DE102014106918A1 (de) Verfahren und Vorrichtung zur automatisierten Durchführung von affinitätsbasierten Assays
EP3926636A2 (de) Verfahren zum erkennen einer amplifikationsphase in einer amplifikation
DE102020215815A1 (de) Verfahren und Vorrichtung zum Trainieren eines Klassifikators für molekularbiologische Untersuchungen
DE102022207161A1 (de) Verfahren zum Quantifizieren einer Lösung und mikrofluidische Analysevorrichtung
EP0861687A1 (de) Verfahren zur Auswertung von Reaktionsverläufen
EP2589950A1 (de) Probenanalysevorrichtung zum Bestimmen von Proben in einer Probenmatrix und Verfahren zum Bestimmen von Proben in einer oder mehreren Probenmatrizen
AT507056B1 (de) Immunochromatographisches verfahren und testsystem zur bestimmung von wenigstens einem analyten in einer zu untersuchenden testlösung
EP3136102A1 (de) Verfahren zur bestimmung eines antikörperbasierten analyten in körperflüssigkeiten eines säugers
DE102019220020A1 (de) Verfahren und Vorrichtung zum Ermitteln einer in einem Fluid enthaltenen Kopienanzahl einer DNA-Sequenz
DE102013113170A1 (de) Verfahren zur Bestimmung eines Messwerts und Analysator zur Ausführung des Verfahrens
DE102017206878A1 (de) Verfahren, Steuergerät und Vorrichtung zur Detektion eines gasförmigen Stoffes in einem Gasgemisch

Legal Events

Date Code Title Description
R163 Identified publications notified
R012 Request for examination validly filed
R016 Response to examination communication
R018 Grant decision by examination section/examining division