IP Library Granted Patent US 10,273,540
Granted Patent B2
US 10,273,540 · App. 13/967,665 · Granted Apr 30, 2019

Methods and apparatuses for estimating parameters in a predictive model for use in sequencing-by-synthesis

Inventors: Melville Davey (Westbrook, CT); Michael Meyer (Ithaca, NY); Marcin Sikora (Foster City, CA); Simon Cawley (Oakland, CA); Kirk Pastorian (Temecula, CA)
Assignee: Life Technologies Corporation
C12Q1/6874G06F19/10
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 10,273,540
App. No.
13/967,665
Granted
Apr 30, 2019
Kind
B2
Abstract

A method of estimating a parameter related to sequencing of a sample nucleic acid template includes: receiving signal data relating to nucleotide incorporation events resulting from a series of flows of nucleotides onto an array of wells including (i) a first well containing the sample nucleic acid template and (ii) a plurality of other sample-containing wells; determining sequence information for the sample nucleic acid template using signal data from the first well; and constructing a phase-state model for a set of nucleotide flows that contributed at least in part to the sequence information, wherein the model includes a signal correction parameter that is determined using signal data from the plurality of other sample-containing wells.

Claims (25)

1. A method of estimating a parameter related to sequencing of a sample nucleic acid template, comprising:

(a) measuring signal data relating to nucleotide incorporation events resulting from a series of flows of nucleotides onto a sensor array comprising a plurality of regions of wells, at least one of the regions of wells comprising: (i) a first set of wells including a first well containing the sample nucleic acid template and (ii) a second set of wells including a plurality of other sample-containing wells, wherein the first set of wells and the second set of wells are physically distinguishable from each other by shape or dimension;

(b) determining sequence information for the sample nucleic acid template using signal data from the first well containing the sample nucleic acid template;

(c) constructing a phase-state model for a set of nucleotide flows that contributed at least in part to the sequence information, wherein the model includes a signal correction parameter that is determined using signal data from the plurality of other sample-containing wells, and wherein the model is stored in a machine-readable memory;

(d) calculating, using the phase-state model, predicted signals for the plurality of other sample-containing wells resulting from the set of nucleotide flows;

(e) comparing the predicted signals to the signal data from the plurality of other sample-containing wells;

(f) fitting the signal correction parameter of the phase-state model based on the comparison of the predicted signals to the signal data from the plurality of other sample-containing wells; and

(g) storing the fitted signal correction parameter in the memory.

2. The method of claim 1 , wherein the signal correction parameter is obtained using signal data obtained from at least a portion of the plurality of other sample-containing wells and without using signal data obtained from the first well.

3. The method of claim 1 , wherein the signal correction parameter is obtained using signal data obtained from at least a portion of the plurality of other sample-containing wells and signal data obtained from the first well.

4. The method of claim 1 , further comprising:

performing steps (b) through (g) for each of the obtained signal data from each or some of the plurality of other sample-containing wells to obtain multiple fitted signal correction parameters, wherein each of the multiple fitted signal correction parameters is determined for a given well without using signal data from that given well.

5. The method of claim 4 , wherein the comparing step comprises calculating a fitting metric that measures the fit between the predicted signals and the signal data.

6. The method of claim 1 , wherein the phase-state model includes two or more signal correction parameters, including a carry forward rate and an incomplete extension rate.

7. The method of claim 1 , wherein the comparing step comprises calculating a fitting metric that measures a fit between the predicted signals and the signal data.

8. The method of claim 7 , wherein the fitting step comprises determining a value of the signal correction parameter that optimizes the fitting metric.

9. The method of claim 7 , wherein the fitting step comprises determining a value of the signal correction parameter using Nelder-Mead optimization.

10. The method of claim 8 , wherein the fitting metric is calculated using only nucleotide flows that result in nucleotide non-incorporation or single nucleotide incorporations.

11. The method of claim 1 , further comprising performing a base calling analysis of the signal data using the fitted signal correction parameter.

12. The method of claim 1 , wherein the set of nucleotide flows is a first set of nucleotide flows and the sequence information is a first sequence information, and further comprising:

applying the phase-state model using the fitted signal correction parameter;

calculating, using the phase-state model and the fitted signal correction parameter obtained using signal data from the plurality of other sample-containing wells, predicted signals for the first well resulting from a second set of nucleotide flows that includes nucleotide flows that are not in the first set of nucleotide flows;

making base calls by comparing the signal data from the first well to the predicted signals for the first well; and

obtaining a second sequence information about the sample nucleic acid template, wherein the second sequence information includes sequence information not contained in the first sequence information.

13. The method of claim 12 , further comprising repeating steps (d) through (g) using the second sequence information to obtain a further fitted signal correction parameter.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jan 21, 2014
From: DAVEY, MELVILLE; MEYER, MICHAEL; SIKORA, MARCIN; CAWLEY, SIMON; PASTORIAN, KIRK
To: LIFE TECHNOLOGIES CORPORATION
Reel/Frame 032010/0200 →
Continuity (5)
Continuation In Part 13283320 · Oct 27, 2011
Provisional Application 61775322 · Mar 8, 2013
Provisional Application 61684221 · Aug 17, 2012
Provisional Application 61407377 · Oct 27, 2010
Related Publication 20140051584A1 · Feb 20, 2014