Methods and systems for digitally counting features on arrays

• US 9,582,877 B2
• Filed: 10/07/2014
• Issued: 02/28/2017
• Est. Priority Date: 10/07/2013
• Status: Active Grant

• Current Assignee: Becton, Dickinson & Co
• Original Assignee: Cellular Research, Inc. (Becton, Dickinson & Co)
• Inventors: Rudoff, Roger, Stern, David, Wu, George T., Fodor, Stephen P. A., Fu, Glenn, Facer, Geoffrey
• Primary Examiner(s): Liew, Alex
• Application Number:
  US14/508,911
• Publication Number:
  US 20160055632A1
• Time in Patent Office: 875 Days
• Field of Search: 382128-134
• US Class Current: 1/1
• CPC Class Codes: G01N 2021/6439 with indicators, stains, dyes, tags, labels, marks
• CPC Class Codes: G01N 21/274 Calibration, base line adjustment, drift correction
• CPC Class Codes: G01N 21/6428 Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" in vivo A61B5/00; immunoassay G01N33/53
• CPC Class Codes: G01N 21/6452 Individual samples arranged in a regular 2D-array, e.g. multiwell plates
• CPC Class Codes: G01N 21/6456 Spatial resolved fluorescence measurements; Imaging
• CPC Class Codes: G01N 2201/129 Using chemometrical methods
• CPC Class Codes: G06F 18/24 Classification techniques
• CPC Class Codes: G06T 2207/10056 Microscopic image
• CPC Class Codes: G06T 2207/10064 Fluorescence image
• CPC Class Codes: G06T 2207/30072 Microarray; Biochip, DNA array; Well plate
• CPC Class Codes: G06T 2207/30242 Counting objects in image
• CPC Class Codes: G06T 7/0012 Biomedical image inspection
• CPC Class Codes: G06T 7/277 involving stochastic approaches, e.g. using Kalman filters

Abstract

Methods, systems and platforms for digital imaging of multiple regions of an array, and detection and counting of the labeled features thereon, are described.

Claims

1. An array reader system comprising

a computer system comprising a digital processor configured to execute computer-executable instructions to at least:

(a) acquiring image data of a region of an array, wherein the array comprises a plurality of stochastically labeled and non-labeled features, and wherein the features comprise oligonucleotide probes;

(b) measuring a signal intensity and a local background intensity for one of the plurality of stochastically labeled and non-labeled features based on the acquired image data;

(c) calculating a local background corrected signal intensity for the one of the plurality of stochastically labeled and non-labeled features using the signal intensity and the local background intensity measured in (b);

(d) transforming the local background corrected signal intensity calculated in (c) for the one of the plurality of stochastically labeled and non-labeled features to determine the number of stochastically labeled and non-labeled features in the of the array; and

(e) calculating the absolute number of a target molecule based on the number of stochastically labeled and non-labeled features determined in (d) within the region of the array.

2. The array reader system of claim 1, further comprising an optical imaging system that comprises one or more image sensors and an illumination system, wherein the illumination system comprises at least one light source;

and

wherein acquiring the image data of the region of the array comprises acquiring the image data of the region of the array from the optical imaging system.

3. The array reader system of claim 1, wherein the array comprises a microarray, microscope slide, or microwell plate.

4. The array reader system of claim 2, wherein the optical imaging system has a magnification of less than 1, equal to 1, or greater than 1.

5. The array reader system of claim 2, wherein the optical imaging system comprises a fluorescence imaging system.

6. The array reader system of claim 2, wherein the optical imaging system comprises a phosphorescence imaging system.

7. The array reader system of claim 2, wherein the optical imaging system comprises an imaging system that operates in a transmitted light, reflected light, or scattered light imaging mode, or combinations thereof.

8. The array reader system of claim 2, wherein the one or more image sensors have a resolution of at least 320×

240 pixels.

9. The array reader system of claim 2, wherein the one or more image sensors are CCD image sensors.

10. The array reader system of claim 2, wherein the one or more image sensors are CMOS Image sensors.

11. The array reader system of claim 2, wherein the one or more image sensors comprise one or more circuit boards.

12. The array reader system of claim 2, wherein the optical imaging system further comprises one or more components selected from the group including, but not limited to, a microscope objective, a camera lens, a finite-conjugate lens, an infinite-conjugate lens, a plano-convex lens, a double convex lens, a plano-concave lens, a double concave lens, an achromatic cemented doublet, or a bandpass filter.

13. The array reader system of claim 5, wherein the fluorescence imaging system is designed for use with fluorescein, Cy3, Cy5, or phycoerythrin fluorophores.

14. The array reader system of claim 2, wherein the at least one light source is an LED or LED assembly.

15. The array reader system of claim 2, wherein the at least one light source is electronically synchronized with the image sensor, the at least one light source being turned on when the image sensor is acquiring an image and turned off when the image sensor is not acquiring an image.

16. The array reader system of claim 2, wherein the illumination system is an off-axis illumination system.

17. The array reader system of claim 16, wherein the off-axis illumination system satisfies the Scheimpflug condition.

18. The array reader system of claim 16, wherein the off-axis illumination system does not satisfy the Scheimpflug condition.

19. The array reader system of claim 16, wherein the off-axis illumination subsystem is a Kohler illumination system.

20. The array reader system of claim 16, wherein the off-axis illumination system is an Abbe illumination system.

21. The array reader system of claim 2, wherein the illumination system is an epi-illumination system.

22. The array reader system of claim 21, wherein the epi-illumination system is a Kohler illumination system.

23. The array reader system of claim 21, wherein the epi-illumination system is an Abbe illumination system.

24. The array reader system of claim 2, wherein the illumination system is a trans-illumination system.

25. The array reader system of claim 24, wherein the trans-illumination system is a Kohler illumination system.

26. The array reader system of claim 24, wherein the trans-illumination system is an Abbe illumination system.

27. The array reader system of claim 2, wherein the optical imaging system further comprises a translation stage.

28. The array reader system of claim 27, wherein the translation stage is a single-axis translation stage.

29. The array reader system of claim 27, wherein the translation stage is a dual-axis translation stage.

30. The array reader system of claim 27, wherein the translation stage is a multi-axis translation stage.

31. The array reader system of claim 1, wherein the digital processor is configured to automatically locates features of the array within the acquired image data.

32. The array reader system of claim 1, wherein calculating the local background corrected signal intensity for the one of the plurality of stochastically labeled and non-labeled features using the signal intensity and the local background intensity measured in (b) comprises:

(i) centering a predefined analysis window on an array feature within an image,

(ii) calculating an intensity value statistic for signal and background pixels according to a predefined pattern of pixels within the array feature, and

(iii) utilizing the intensity value statistic for signal and background pixels to calculate the local background corrected signal intensity for the array feature.

33. The array reader system of claim 32, wherein the digital processor is configured to perform a k-means clustering analysis of the background corrected signal intensity values for the array features to determine a dynamic signal intensity threshold for discrimination between the plurality of stochastically labeled and non-labeled features of the array.

34. The array reader system of claim 32, wherein the digital processor is configured to perform a k-medoids clustering analysis of the background corrected signal intensity values for the array features to determine a dynamic signal intensity threshold for discrimination between the plurality of stochastically labeled and non-labeled features of the array.

35. The array reader system of claim 32, wherein the digital processor is configured to perform a mixture model statistical analysis of the background corrected signal intensity values for the array features to determine a dynamic signal intensity threshold for discrimination between the plurality of stochastically labeled and non-labeled features of the array.

36. The array reader system of claim 32, wherein the digital processor is configured to perform an empirical analysis based on sorting of background corrected signal intensity values for the array features to determine a dynamic signal intensity threshold for discrimination between the plurality of stochastically labeled and non-labeled features of the array.

37. The array reader system of claim 32, wherein the digital processor is configured to perform an empirical analysis based on sorting of pairwise differences in background corrected signal intensity values for the array features to determine a dynamic signal intensity threshold for discrimination between the plurality of stochastically labeled and non-labeled features of the array.

38. The array reader system of claim 32, wherein the digital processor is configured to perform one or more statistical analyses of the background corrected signal intensity values for the array features to determine a dynamic signal intensity threshold for discrimination between the plurality of stochastically labeled and non-labeled features of the array, and wherein the one or more statistical analyses are selected from the list comprising k-means clustering, k-medoids clustering, mixture model statistical analysis, an empirical analysis, or any combination thereof.

39. The array reader system of claim 1, wherein calculating the absolute number of the target molecule based on the number of stochastically labeled and non-labeled features determined in (d) within the region of the array comprises calculate the absolute number of the target molecule based on the number of stochastically labeled features, the number of non-labeled features within the region of the array, and the predictions of the Poisson distribution.

40. The array reader system of claim 1, wherein the digital processor is configured to calculate a confidence interval for the number of target molecules.

41. The array reader system of claim 2, wherein the optical imaging system and computer system are combined within a single, stand-alone instrument.

42. The array reader system of claim 2, wherein the optical imaging system and computer system are configured as separate instrument modules.

43. The array reader system of claim 32, wherein the digital processor is configured to perform an analysis of the local background corrected signal intensities for the array features to determine a dynamic signal intensity threshold, and wherein the analysis comprises fitting a model function to the local background corrected signal intensities by varying model parameters.

44. The array reader system of claim 32, wherein the digital processor is configured to perform an analysis of the local background corrected signal intensities for the array features to determine a dynamic signal intensity threshold, and wherein the analysis comprises maximizing a quality metric relating to a statistical difference between feature intensities above the threshold and feature intensities below the threshold.