IP Library › Granted Patent US 12,423,815
Granted Patent B2
US 12,423,815 · App. 17/858,268 · Granted Sep 23, 2025

Methods and systems for real time extraction of crosstalk in illumination emitted from reaction sites

Inventor: Mohsen Rezaei (Brown Deer, WI)
Assignee: ILLUMINA, INC.
G06T7/0014G01N33/582G06T5/73G16B30/00G06T2207/20224
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 12,423,815
App. No.
17/858,268
Granted
Sep 23, 2025
Kind
B2
Abstract

Biosensor including an array of reaction sites and corresponding light sensors may experience crosstalk in which photons from one reaction site are detected by neighbors of its corresponding light sensor, and such crosstalk may be corrected using sharpening kernels corresponding to the sensors in the array. Such sharpening kernels may be derived from point spread functions, which may be determined in real time analysis based on images captured during sequencing.

Claims (93)

1. A method comprising:

obtaining a plurality of analysis images of light emitted during sequencing of a biological sample;

obtaining noise dependencies by performing acts comprising, for each location in a point spread function which comprises a plurality of locations, calculating a noise correlation for that location, wherein the noise correlation is a correlation between noise in a first plurality of measurements and noise in a second plurality of measurements, wherein:

each measurement from the first plurality of measurements is captured by a sensor whose position relative to a corresponding sensor which captured a measurement from the second plurality of measurements is the same as that location's position in the point spread function relative to a center of the point spread function; and

each measurement from the first plurality of measurements and the second plurality of measurements measures light emitted during sequencing of the biological sample from the plurality of analysis images;

populating the point spread function based on the noise dependencies;

generating a sharpening kernel based on the point spread function;

obtaining a plurality of sharpened images by applying the sharpening kernel to the plurality of analysis images;

one or more times, repeating:

obtaining noise dependencies;

populating the point spread function;

generating the sharpening kernel; and

obtaining the plurality of sharpened images;

wherein, on each repetition, the plurality of analysis images for that repetition of obtaining noise dependencies is the plurality of sharpened images from a most recent preceding application of the sharpening kernel;

identifying a sharpening kernel generated on a repetition as an optimal sharpening kernel; and

applying the optimal sharpening kernel to compensate for crosstalk in images subsequently captured while sequencing the biological sample.

2. The method of claim 1 , wherein repeating populating the point spread function comprises:

generating a dependency matrix having dimensions equal to those of the point spread function by, for each location in the point spread function, populating a corresponding location in the dependency matrix with a most recently obtained dependency for that location;

multiplying the dependency matrix by a scalar constant; and

adding a result of multiplying the dependency matrix by the scalar constant to a most recently populated preceding point spread function.

3. The method of claim 2 , wherein the scalar constant has a value greater than or equal to 0.08, and less than or equal to 0.12.

4. The method of claim 1 , wherein identifying the optimal sharpening kernel comprises:

determining, for each repetition, a signal to noise ratio obtained by applying the sharpening kernel generated on that repetition to the plurality of analysis images; and

identifying the sharpening kernel from which a highest signal to noise ratio is obtained as the optimal sharpening kernel.

5. The method of claim 4 , wherein determining, for each repetition, the signal to noise ratio comprises calculating a sharpness of the plurality of sharpened images obtained on that repetition.

6. The method of claim 1 , wherein obtaining noise dependencies comprises obtaining a noise map based on subtracting a first image from the plurality of analysis images from a second image from the plurality of analysis images.

7. The method of claim 6 , wherein:

obtaining noise dependencies comprises dividing the noise map into a plurality of units, wherein each unit is a matrix of values having dimensions at least as great as those of the point spread function;

for each location in the point spread function, obtaining the noise dependency for that location comprises calculating a correlation between a first set of values and a second set of values, wherein:

the first set of values comprises, for each unit from the plurality of units, a value at a first location in that unit; and

the second set of values comprises, for each unit from the plurality of units, a value at a second location in that unit, wherein the first location for the value from that unit in the first set of values has a position relative to the second location that is the same as that location's position in the point spread function relative to the center of the point spread function.

8. The method of claim 6 , wherein:

the first image from the plurality of analysis images is an image captured from a first sequencing cycle; and

the second image from the plurality of analysis images is an image from a second sequencing cycle.

9. The method of claim 1 , wherein:

each image from the plurality of analysis images comprises an image from a different sequencing cycle;

the plurality of analysis images comprises more than two images;

obtaining noise dependencies comprises:

obtaining a plurality of intermediate correlations, wherein each of the intermediate correlations corresponds to two analysis images from the plurality of analysis images, and wherein each of the plurality of intermediate correlations is obtained based on:

obtaining an intermediate map for that intermediate correlation by subtracting one of the analysis images corresponding to that intermediate correlation from the other analysis image corresponding to that intermediate correlation;

dividing the intermediate map for that intermediate correlation into a plurality of units, wherein each unit is a matrix of values having dimensions at least as great as those of the point spread function;

for each location in the point spread function, calculating a correlation between a first set of values and a second set of values, wherein:

the first set of values comprises, for each unit from the plurality of units, a value at a first location in that unit; and

the second set of values comprises, for each unit from the plurality of units, a value at a second location in that unit, wherein the first location for the value from that unit in the first set of values has a position relative to the second location that is the same as that location's position in the point spread function relative to the center of the point spread function;

determining the set of noise dependencies based on the plurality of intermediate correlations.

10. The method of claim 1 , wherein obtaining noise dependencies, populating the point spread function, generating the sharpening kernel, and obtaining the plurality of sharpened images are repeated between two and eight times.

11. A system comprising:

a sensor array;

a processor to:

obtain a plurality of analysis images of light emitted during sequencing of a biological sample;

obtain noise dependencies by performing acts comprising, for each location in a point spread function which comprises a plurality of locations, calculating a noise correlation for that location, wherein the noise correlation is a correlation between noise in a first plurality of measurements and noise in a second plurality of measurements, wherein:

each measurement from the first plurality of measurements is captured by a sensor from the sensor array whose position relative to a corresponding sensor in the sensor array which captured a measurement from the second plurality of measurements is the same as that location's position in the point spread function relative to a center of the point spread function; and

each measurement from the first plurality of measurements and the second plurality of measurements measures light emitted during sequencing of the biological sample from the plurality of analysis images;

populate the point spread function based on the noise dependencies;

generate a sharpening kernel based on the point spread function;

obtain a plurality of sharpened images by applying the sharpening kernel to the plurality of analysis images;

one or more times, repeat:

obtaining noise dependencies;

populating the point spread function;

generating the sharpening kernel; and

obtaining the plurality of sharpened images;

wherein, on each repetition, the plurality of analysis images for that repetition of obtaining noise dependencies is the plurality of sharpened images from a most recent preceding application of the sharpening kernel;

identify a sharpening kernel generated on a repetition as an optimal sharpening kernel; and

apply the optimal sharpening kernel to compensate for crosstalk in images subsequently captured while sequencing the biological sample.

12. The system of claim 11 , wherein repeating populating the point spread function comprises:

generating a dependency matrix having dimensions equal to those of the point spread function by, for each location in the point spread function, populating a corresponding location in the dependency matrix with a most recently obtained dependency for that location;

multiplying the dependency matrix by a scalar constant; and

adding a result of multiplying the dependency matrix by the scalar constant to a most recently populated preceding point spread function.

13. The system of claim 12 , wherein the scalar constant has a value greater than or equal to 0.08, and less than or equal to 0.12.

14. The system of claim 11 , wherein identifying the optimal sharpening kernel comprises:

determining, for each repetition, a signal to noise ratio obtained by applying the sharpening kernel generated on that repetition to the plurality of analysis images; and

identifying the sharpening kernel from which a highest signal to noise ratio is obtained as the optimal sharpening kernel.

15. The system of claim 14 , wherein determining, for each repetition, the signal to noise ratio comprises calculating a sharpness of the plurality of sharpened images obtained on that repetition.

16. The system of claim 11 , wherein obtaining noise dependencies comprises obtaining a noise map based on subtracting a first image from the plurality of analysis images from a second image from the plurality of analysis images.

17. The system of claim 16 , wherein:

obtaining noise dependencies comprises dividing the noise map into a plurality of units, wherein each unit is a matrix of values having dimensions at least as great as those of the point spread function;

for each location in the point spread function, obtaining the noise dependency for that location comprises calculating a correlation between a first set of values and a second set of values, wherein:

the first set of values comprises, for each unit from the plurality of units, a value at a first location in that unit; and

the second set of values comprises, for each unit from the plurality of units, a value at a second location in that unit, wherein the first location for the value from that unit in the first set of values has a position relative to the second location that is the same as that location's position in the point spread function relative to the center of the point spread function.

18. The system of claim 16 , wherein:

the first image from the plurality of analysis images is an image from a first sequencing cycle; and

the second image from the plurality of analysis images is an image from a second sequencing cycle.

19. The system of claim 11 , wherein:

each image from the plurality of analysis images comprises an image from a different sequencing cycle;

the plurality of analysis images comprises more than two images;

obtaining noise dependencies comprises:

obtaining a plurality of intermediate correlations, wherein each of the intermediate correlations corresponds to two analysis images from the plurality of analysis images, and wherein each of the plurality of intermediate correlations is obtained based on:

obtaining an intermediate map for that intermediate correlation by subtracting one of the analysis images corresponding to that intermediate correlation from the other analysis image corresponding to that intermediate correlation;

dividing the intermediate map for that intermediate correlation into a plurality of units, wherein each unit is a matrix of values having dimensions at least as great as those of the point spread function;

for each location in the point spread function, calculating a correlation between a first set of values and a second set of values, wherein:

the first set of values comprises, for each unit from the plurality of units, a value at a first location in that unit; and

the second set of values comprises, for each unit from the plurality of units, a value at a second location in that unit, wherein the first location for the value from that unit in the first set of values has a position relative to the second location that is the same as that location's position in the point spread function relative to the center of the point spread function;

determining the set of noise dependencies based on the plurality of intermediate correlations.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Dec 6, 2023
From: REZAEI, MOHSEN
To: ILLUMINA, INC.
Reel/Frame 065781/0616 →
Continuity (2)
Provisional Application 63221236 · Jul 13, 2021
Related Publication 20230027409A1 · Jan 26, 2023
References Cited (62)
US 5641658A · Adams et al. · 1997 [cited by applicant]
US 6090592A · Adams et al. · 2000 [cited by applicant]
US 7057026B2 · Barnes et al. · 2006 [cited by applicant]
US 7115400B1 · Adessi et al. · 2006 [cited by applicant]
US 7211414B2 · Hardin et al. · 2007 [cited by applicant]
US 7315019B2 · Turner et al. · 2008 [cited by applicant]
US 7329492B2 · Hardin et al. · 2008 [cited by applicant]
US 7405281B2 · Xu et al. · 2008 [cited by applicant]
US 7414116B2 · Milton et al. · 2008 [cited by applicant]
US 7427673B2 · Balasubramanian et al. · 2008 [cited by applicant]
US 7541444B2 · Milton et al. · 2009 [cited by applicant]
US 7566537B2 · Balasubramanian et al. · 2009 [cited by applicant]
US 7592435B2 · Milton et al. · 2009 [cited by applicant]
US 7595883B1 · El Gamal et al. · 2009 [cited by applicant]
US 8396320B2 · Allebach et al. · 2013 [cited by applicant]
US 8638096B2 · Zhang et al. · 2014 [cited by applicant]
US 8830363B2 · Jang · 2014 [cited by applicant]
US 9453258B2 · Kain et al. · 2016 [cited by applicant]
US 9696408B2 · Eikenberry et al. · 2017 [cited by applicant]
US 9747672B2 · Oniki et al. · 2017 [cited by applicant]
US 10304189B2 · Garcia et al. · 2019 [cited by applicant]
US 11935311B2 · Egertson · 2024 [cited by examiner]
US 20020055100A1 · Kawashima et al. · 2002 [cited by applicant]
US 20040002090A1 · Mayer et al. · 2004 [cited by applicant]
US 20040096853A1 · Mayer · 2004 [cited by applicant]
US 20050100900A1 · Kawashima et al. · 2005 [cited by applicant]
US 20050244870A1 · Chee et al. · 2005 [cited by applicant]
US 20060188901A1 · Barnes et al. · 2006 [cited by applicant]
US 20060240439A1 · Smith et al. · 2006 [cited by applicant]
US 20060281109A1 · Barr Ost et al. · 2006 [cited by applicant]
US 20070128624A1 · Gormley et al. · 2007 [cited by applicant]
US 20070166705A1 · Milton et al. · 2007 [cited by applicant]
US 20080009420A1 · Schroth et al. · 2008 [cited by applicant]
US 20080108082A1 · Rank et al. · 2008 [cited by applicant]
US 20080234136A1 · Drmanac et al. · 2008 [cited by applicant]
US 20080242560A1 · Gunderson et al. · 2008 [cited by applicant]
US 20110059865A1 · Smith et al. · 2011 [cited by applicant]
US 20180060635A1 · Li et al. · 2018 [cited by applicant]
US 20190005351A1 · Zhou et al. · 2019 [cited by applicant]
US 20190333197A1 · Kask · 2019 [cited by examiner]
US 20200080142A1 · Langlois et al. · 2020 [cited by applicant]
CN 103927740A · 2014 [cited by applicant]
CN 104574423A · 2015 [cited by applicant]
EP 2360638A1 · 2011 [cited by examiner]
WO WO1991006678A1 · 1991 [cited by applicant]
WO WO2004018497A2 · 2004 [cited by applicant]
WO WO2005065814A1 · 2005 [cited by applicant]
WO WO2006064199A1 · 2006 [cited by applicant]
WO WO2007010251A2 · 2007 [cited by applicant]
WO WO2007123744A2 · 2007 [cited by applicant]
WO WO2007135368A2 · 2007 [cited by applicant]
WO WO2012058096A1 · 2012 [cited by applicant]
WO 2020061237A1 · 2020 [cited by applicant]
Clenmens, “Flow Imaging”, 2002. (Year: 2002). [cited by examiner]
Bentley, David R., et al. “Accurate whole human genome sequencing using reversible terminator chemistry.” [cited by applicant]
U.S. Appl. No. 61/538,294, entitled “Methods and Compositions for Nucleic Acid Sequencing,” filed Sep. 23, 2011. [cited by applicant]
U.S. Appl. No. 61/619,878, entitled “Methods and Compositions for Nucleic Acid Sequencing,” filed Apr. 3, 2012. [cited by applicant]
U.S. Appl. No. 62/216,125, entitled “Methods and Systems to Correct Crosstalk in Illumination Emitted from Reaction Sites,” filed Jun. 29, 2021. [cited by applicant]
Brauers, Johannes, Claude Seiler, and Til Aach. “Direct PSF estimation using a random noise target.” [cited by applicant]
Salama, Khaled. [cited by applicant]
Zhan, Dazhi, et al. “PSF estimation method of simple-lens camera using normal sinh-arcsinh model based on noise image pairs.” [cited by applicant]
International Search Report and Written Opinion dated Dec. 5, 2022, for International Application No. PCT/US2022/036193, 15 pages. [cited by applicant]