Models for analyzing data from sequencing-by-synthesis operations

• US 10,679,724 B2
• Filed: 04/27/2017
• Issued: 06/09/2020
• Est. Priority Date: 05/11/2012
• Status: Active Grant

• Current Assignee: Life Technologies Corporation (Thermo Fisher Scientific Incorporated)
• Original Assignee: Life Technologies Corporation (Thermo Fisher Scientific Incorporated)
• Inventors: Rearick, Todd
• Primary Examiner(s): Negin, Russell S
• Application Number:
  US15/498,591
• Publication Number:
  US 20170293713A1
• Time in Patent Office: 1,139 Days
• Field of Search: None
• CPC Class Codes: C12Q 1/6874 involving nucleic acid arrays, e.g. sequencing by hybridisation
• CPC Class Codes: C12Q 2537/165 Mathematical modelling, e.g. logarithm, ratio
• CPC Class Codes: C12Q 2563/113 the label being electroactive, e.g. redox labels
• CPC Class Codes: C12Q 2563/116 electrical properties of nucleic acids, e.g. impedance, conductivity or resistance
• CPC Class Codes: G16B 25/00 ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
• CPC Class Codes: G16B 30/00 ICT specially adapted for sequence analysis involving nucleotides or amino acids
• CPC Class Codes: G16B 5/00 ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks

Abstract

A method of modeling a background signal when sequencing a polynucleotide strand using sequencing-by-synthesis includes: flowing a series of nucleotide flows onto a reactor array having multiple reaction confinement regions, one or more copies of the polynucleotide strand being located in a loaded reaction confinement region of the reactor array, the loaded reaction confinement region being located in a vicinity of one or more neighboring reaction confinement regions that may or may not be loaded; receiving output signals from the reactor array; and modeling a background signal for the loaded reaction confinement region using the received output signals and a model adapted to account at least for an exchange of ions between the one or more neighboring reaction confinement regions and a headspace adjacent the loaded reaction confinement region and the one or more neighboring reaction confinement regions.

Claims

1. A method of modeling a background signal when sequencing a polynucleotide strand using sequencing-by-synthesis, comprising:

flowing a series of nucleotide flows onto a reactor array having multiple reaction confinement regions, one or more copies of the polynucleotide strand being located in a loaded reaction confinement region of the reactor array, the loaded reaction confinement region being located in a vicinity of one or more neighboring reaction confinement regions;

receiving output signals from the reactor array; and

estimating a background signal for the loaded reaction confinement region using the received output signals and a background model representing at least an exchange of ions between the one or more neighboring reaction confinement regions and a headspace adjacent to the loaded reaction confinement region and the one or more neighboring reaction confinement regions, wherein the background model is derived using a first characteristic equation, wherein the first characteristic equation comprises a term related to a difference between a proton concentration in each of the neighboring reaction confinement regions and a proton concentration in the headspace.

2. The method of claim 1, wherein the background model represents at least an exchange of ions between at least two neighboring reaction confinement regions and the headspace.

3. The method of claim 1, wherein the background model represents at least an exchange of ions between at least three neighboring reaction confinement regions and the headspace.

4. The method of claim 1, wherein the background model represents at least an exchange of ions between at least four neighboring reaction confinement regions and the headspace.

5. The method of claim 1, wherein the background model represents at least an exchange of ions between six neighboring reaction confinement regions and the headspace, the six neighboring reaction confinement regions being arranged hexagonally around the loaded reaction confinement region.

6. The method of claim 1, wherein the background model simulates a relationship between the loaded reaction confinement region, the neighboring reaction confinement regions, and a bulk fluid, all interacting with each other through the headspace.

7. The method of claim 1, wherein the background model represents at least an exchange of ions between the headspace and a bulk fluid.

8. The method of claim 1, wherein the first characteristic equation represents a conservation of protons in the headspace.

9. The method of claim 8, wherein the first characteristic equation is:

10. The method of claim 8, wherein the first characteristic equation comprises a term related to a difference between a proton concentration in a bulk fluid and a proton concentration in the headspace.

11. The method of claim 8, wherein the first characteristic equation comprises a term related to a difference between a proton concentration in the loaded reaction confinement region and a proton concentration in the headspace.

12. The method of claim 8, wherein the background model is further derived using a second characteristic equation representing a conservation of protons in the loaded reaction confinement region.

13. The method of claim 12, wherein the second characteristic equation is:

14. The method of claim 1, wherein the background model comprises an interaction equation representing a relationship between the loaded reaction confinement region, the neighboring reaction confinement regions, and a bulk fluid, all interacting with each other through the headspace.

15. The method of claim 14, wherein the interaction equation is:

16. The method of claim 1, wherein the background model comprises a simplified equation representing a relationship between the loaded reaction confinement region, the neighboring reaction confinement regions, and a bulk fluid, all interacting with each other through the headspace, the simplified equation being:

17. The method of claim 16, wherein the background model comprises scaling factors a_nn multiplied with proton concentration terms (S_nn-c−

S_nn-e) in the simplified equation corresponding to different neighboring reaction confinement regions, wherein the scaling factors represent unequal contributions of ions exchanged with the headspace by the one or more neighboring reaction confinement regions.

18. A non-transitory machine-readable storage medium comprising instructions which, when executed by a processor, cause the processor to:

receive output signals relating to chemical reactions resulting from the flow of a series of nucleotides onto a reactor array having multiple reaction confinement regions, one or more copies of a polynucleotide strand being located in a loaded reaction confinement region of the reactor array, the loaded reaction confinement region being located in a vicinity of one or more neighboring reaction confinement regions; and

estimate a background signal for the loaded reaction confinement region using the received output signals and a background model representing at least an exchange of ions between the one or more neighboring reaction confinement regions and a headspace adjacent to the loaded reaction confinement region and the one or more neighboring reaction confinement regions, wherein the background model is derived using a first characteristic equation, wherein the first characteristic equation comprises a term related to a difference between a proton concentration in each of the neighboring reaction confinement regions and a proton concentration in the headspace.

19. The non-transitory machine-readable storage medium of claim 18, wherein the background model represents at least an exchange of ions between six neighboring reaction confinement regions and the headspace, the six neighboring reaction confinement regions being arranged hexagonally around the loaded reaction confinement region.

20. The non-transitory machine-readable storage medium of claim 18, wherein the background model simulates a relationship between the loaded reaction confinement region, the neighboring reaction confinement regions, and a bulk fluid, all interacting with each other through the headspace.