Methods And Compositions For Base Calling Nucleic Acids

• US 20100063743A1
• Filed: 03/17/2009
• Published: 03/11/2010
• Est. Priority Date: 03/19/2008
• Status: Active Grant

• Current Assignee: QIAGEN Sciences, LLC (Qiagen NV)
• Original Assignee: Intelligent Bio-Systems Incorporated. (Qiagen NV)
• Inventors: GORDON, STEVEN, VEATCH, PHILLIP A.
• Granted Patent: US 8,612,161 B2
• Time in Patent Office: Days
• US Class Current: 702/19
• CPC Class Codes: C12Q 1/6806 Preparing nucleic acids for analysis, e.g. for polymerase chain reaction [PCR] assay C12Q1/6804 takes precedence
• CPC Class Codes: C12Q 1/6825 Nucleic acid detection involving sensors
• CPC Class Codes: C12Q 1/6869 Methods for sequencing
• CPC Class Codes: C12Q 2525/117 incorporating modified base
• CPC Class Codes: C12Q 2525/186 incorporating a non-extendable or blocking moiety
• CPC Class Codes: C12Q 2527/137 Concentration of a component of medium
• CPC Class Codes: C12Q 2537/165 Mathematical modelling, e.g. logarithm, ratio
• CPC Class Codes: C12Q 2563/107 fluorescence
• CPC Class Codes: C12Q 2565/601 being a microscope, e.g. atomic force microscopy [AFM]
• CPC Class Codes: G16B 25/00 ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression

Abstract

The invention provides methods and compositions, including, without limitation, algorithms, computer readable media, computer programs, apparatus, and systems for determining the identity of nucleic acids in nucleotide sequences using, for example, data obtained from sequencing by synthesis methods. The methods of the invention include correcting one or more phenomena that are encountered during nucleotide sequencing, such as using sequencing by synthesis methods. These phenomena include, without limitation, sequence lead, sequence lag, spectral crosstalk, and noise resulting from variations in illumination and/or filter responses.

Claims

1. A method for determining an identity of a nucleic acid at an interrogation position in a nucleotide sequence from data acquired from one or more channels, comprising

a) obtaining a data set for one or more probe intensities at one or more nucleic acid positions in said sequence, wherein each probe corresponds to a nucleic acid,

b) determining the ratio contribution to probe intensity at said interrogation position from probe intensities at the interrogation position and at one or both of

i) at least one subsequent nucleic acid positions in said sequence, and

ii) at least one preceding nucleic acid positions in said sequence, and

c) applying said ratio contribution to probe intensity to said data set to arrive at an identity for a nucleic acid at said interrogation position in said nucleotide sequence.

2. The method of claim 1, wherein said determining the ratio contribution to probe intensity comprises measuring the rate at which a lag occurs at one or more nucleotide position in said nucleotide sequence.

3. The method of claim 2, wherein said determining the ratio contribution to probe intensity comprises measuring the rate at which a lag occurs at each nucleotide position in said nucleotide sequence.

4. The method of claim 1, wherein said determining the ratio contribution to probe intensity comprises measuring the rate at which a lead occurs at one or more nucleotide positions in said nucleotide sequence.

5. The method of claim 1, further comprising calling a nucleic acid at said interrogation position in said nucleotide sequence.

6. The method of claim 1, further comprising repeating steps b) and c) to arrive at an identity for a nucleic acid at more than one interrogation position in said nucleotide sequence.

7. The method of claim 1, further comprising

a) applying a sequence lead-lag compensation equation to determine the ratio contribution to probe intensity from probe at

i) said interrogation position,

ii) each position preceding said interrogation position, and

iii) each position subsequent to said interrogation position, and

b) summing up said ratio contribution to probe intensity.

8. The method of claim 1, further comprising applying a sequence lead-lag compensation equation to said ratio contribution to probe intensity at a plurality of positions in said sequence.

9. The method of claim 8, wherein said sequence lead-lag compensation equation is determined by applying equation

10. The method of claim 8, wherein said sequence lead-lag compensation equation is determined by applying equation

11. The method of claim 10, wherein said sequence lead-lag compensation equation is determined by applying equation

12. The method of claim 1, wherein the nucleic acid comprises a base selected from the group consisting of adenine (A), guanine (G), cytosine (C), and thymine (T).

13. The method of claim 1, wherein said probe is fluorescent.

14. The method of claim 1, further comprising field flattening of background data for said data set.

15. The method of claim 14, wherein said field flattening comprises

a) obtaining a first data set for a plurality of pixel intensities of a first raw image of a probe at a first concentration on a solid support, wherein said first raw image is produced using a first spectral filter for detecting a first probe,

b) obtaining a second data set for a plurality of pixel intensities of a second smoothed image of said probe on said solid support, wherein said second smoothed image is produced using a low pass filter,

c) determining a field flattening intensity value for a plurality of pixels of said first raw image, and

d) generating a third field flattened image of said probe on said solid support using said field flattening intensity of said plurality of pixels, wherein the correlation of intensity of a plurality of pixels to their spatial location on said third field flattened image is reduced compared to the intensity of a plurality of pixels at a corresponding location on said first raw image.

16. The method of claim 15, wherein said field flattening intensity value of a pixel is determined by equation

F_x,y=R_x,yM_x0,y0/M_x,y

where

F_x,y is a field flattening intensity value,

R_x,y is the intensity of a pixel of the plurality of pixels on the first raw image,

M_x,y is the intensity of a pixel of the plurality of pixels on the second smoothed image at a corresponding spatial location to the pixel on the first raw image, and

M_x0,y0 is the intensity of a reference pixel on said second smoothed image or is an arbitrary scale factor.

17. The method of claim 15, further comprising repeating steps a) to d), using a second spectral filter for detecting a second probe.

18. The method of claim 15, further comprising repeating steps a) to d), using said probe at a second concentration on said solid support.

19. The method of claim 14, wherein said solid support comprises a microscope slide.

20. The method of claim 14, wherein said solid support comprises a silicon chip.

21. The method of claim 1, further comprising reducing spectral crosstalk in said one or more channels.

22. The method of claim 21, wherein said reducing spectral crosstalk comprises

a) determining spectral crosstalk factors for each of said one or more probes in its corresponding channel from one or more adjacent channels,

b) applying said spectral crosstalk factors to determine a spectral crosstalk matrix, and

c) applying said spectral crosstalk matrix to said data set for said one or more probe intensities.

23. The method of claim 21, wherein said reducing spectral crosstalk comprises

a) determining probe intensity for one or more probes from one or more channels, wherein each channel corresponds to a probe;

b) determining the ratios of said probe intensities in said one or more channels to arrive at a signature ratio for said probe intensity in its corresponding channel,

c) applying said signature ratios in a matrix equation, and

d) inverting said matrix equation to arrive at an inverted matrix.

24. The method of claim 23, further comprising e) applying said inverted matrix to data from said one or more channels.

25. The method of claim 22, wherein said determining spectral crosstalk matrix comprises using equation

34. A computer readable medium containing a computer program for performing the method of claim 1.

26. The method of claim 26, wherein said equation is solved to determine spectral crosstalk matrix K⁻

1 using equation

27. An algorithm for processing data for nucleic acids in a nucleotide sequence, wherein said data is acquired from one or more channels, the algorithm comprising

a) determining the ratio contribution to probe intensity in said one or more channels for one or more interrogation positions, from probe intensities at the interrogation position and at one or both of

iv) at least one subsequent nucleic acid positions in said sequence, and

v) at least one preceding nucleic acid positions in said sequence,

b) processing data from said one or more channels to correct for sequence lead and sequence lag, and

c) reconstructing said data in said one or more channels.

28. The algorithm of claim 27, wherein said processing data comprises applying said ratio contribution to probe intensity to determine, for said probe at said one or more interrogation positions, a sequence lead-lag compensation equation.

29. The algorithm of claim 28, wherein said sequence lead-lag compensation equation is determined by applying equation

30. The algorithm of claim 27, wherein said sequence lead-lag compensation equation is determined by applying equation

31. The algorithm of claim 30, wherein said sequence lead-lag compensation equation is determined by applying equation

32. The algorithm of claim 27, further comprising field flattening of background data.

33. The algorithm of claim 27, further comprising reducing spectral crosstalk between said data comprised in a plurality of channels.

35. A computer program product for processing data for nucleic acids in a nucleotide sequence to determine an identity of a nucleic acid at an interrogation position in said nucleotide sequence, said computer program product comprising

a) computer code that inputs data from one or more channels for one or more probe intensities, wherein each channel corresponds to a probe, and each probe corresponds to a nucleic acid,

b) computer code that applies to the input data a sequence lead-lag compensation equation to correct for sequence lead and sequence lag,

c) computer code that compares probe intensities in said one or more channels that have been corrected for sequence lead and sequence lag,

d) computer code that determines the highest probe intensity of the compared probe intensities, and

e) computer code that identifies a nucleic acid at said interrogation position according to said highest probe intensity.

36. The computer program product of claim 35, further comprising f) computer code that applies field flattening of background data.

37. The computer program product of claim 35, further comprising f) computer code that reduces spectral crosstalk between data comprised in said one or more channels.

38. An apparatus that processes data for nucleic acids in a nucleotide sequence to determine an identity of a nucleic acid at an interrogation position in said nucleotide sequence, said apparatus comprising

a) means for inputting data from one or more channels for one or more probe intensities, wherein each channel corresponds to a probe, and each probe corresponds to a nucleic acid,

b) means for applying to the input data a sequence lead-lag compensation equation to correct for sequence lead and sequence lag,

c) means for comparing probe intensities in said one or more channels that have been corrected for sequence lead and sequence lag,

d) means for determining the highest probe intensity of the compared probe intensities, and

e) means for identifying a nucleic acid at said interrogation position according to said highest probe intensity.

39. The apparatus of claim 38, further comprising f) means for applying field flattening of background data.

40. The apparatus of claim 38, further comprising f) means for reducing spectral crosstalk between data comprised in said one or more channels.

41. A system for processing data to determine an identity of a nucleic acid at an interrogation position in said nucleotide sequence, said system comprising

a) a processor, and

b) a computer readable medium readable by said processor, said computer readable medium storing a computer program that comprisescode that receives as input a plurality of nucleic acid base calls at an interrogation position in a nucleotide sequence,

code that applies to the input data a sequence lead-lag compensation equation to correct for sequence lead and sequence lag, andcode that identifies a nucleic acid at said interrogation position according to the corrected data.

42. The system of claim 41, wherein said computer readable medium further comprises 4) code that applies field flattening of background data.

43. The system of claim 41, said computer readable medium further comprises 4) code that reduces spectral crosstalk between data comprised in said one or more channels.

44. A method for field flattening an image of a probe on a solid support, comprising

a) obtaining a first data set for a plurality of pixel intensities of a first raw image of a probe at a first concentration on a solid support, wherein said first raw image is produced using a first spectral filter for detecting a first probe,

b) obtaining a second data set for a plurality of pixel intensities of a second smoothed image of said probe on said solid support, wherein said second smoothed image is produced using a low-pass filter,

c) determining a field flattening intensity value for a plurality of pixels of said first raw image, and

d) generating a third field flattened image of said probe on said solid support using said field flattening intensity of said plurality of pixels, wherein the correlation of intensity of a plurality of pixels to their spatial location on said third field flattened image is reduced compared to the intensity of a plurality of pixels at a corresponding location on said first raw image.

45. The method of claim 44, wherein said field flattening intensity value of a pixel is determined by equation

F_x,y=R_x,yM_x0,y0/M_x,y

where

F_x,y is a field flattening intensity value,

R_x,y is the intensity of a pixel of the plurality of pixels on the first raw image,

M_x,y is the intensity of a pixel of the plurality of pixels on the second smoothed image at a corresponding spatial location to the pixel on the first raw image, and

M_x0,y0 is the intensity of a reference pixel on said second smoothed image or is any other scale factor of interest.

46. The method of claim 44, further comprising repeating steps a) to d), using a second spectral filter for detecting a second probe.

47. The method of claim 44, further comprising repeating steps a) to d), using said probe at a second concentration on said solid support.

48. The method of claim 44, wherein said probe is fluorescent.

49. The method of claim 48, wherein said probe corresponds to a nucleic acid.

50. The method of claim 49, wherein said nucleic acid comprises a base selected from the group consisting of adenine (A), guanine (G), cytosine (C), and thymine (T).

51. The method of claim 44, wherein said solid support comprises a microscope slide.

52. The method of claim 44, wherein said solid support comprises a silicon chip.