Primary analysis in next generation sequencing
Image data analysis, and particularly identifying cluster or polony locations for performing base-calling in a digital image of a flow cell during DNA sequencing is described. A method may include generating a first plurality of flow cell images of a cellular sample immobilized on a support by conducting one or more cycles of sequencing reactions. The cellular sample may include a plurality of concatemer molecules therewithin. For the first plurality of flow cell image, pixel intensities, and a respective color purity of each of the pixel intensities may be determined. A base calling template may include base calling locations based on the pixel intensities and the respective color purity of the pixel intensities. The base calling template may be for registering a second plurality of flow cell images of the support in one or more subsequent cycles of the one or more cycles.
1 . A method comprising:
generating, by a sequencing system, a first plurality of flow cell images of a cellular sample immobilized on a support by conducting one or more cycles of sequencing reactions, the cellular sample comprising a plurality of concatemer molecules therewithin, wherein a first concatemer molecule of the plurality of concatemer molecules corresponds to a first target RNA molecule of the cellular sample, and a second concatemer molecule of the plurality of concatemer molecules corresponds to a second target RNA molecule of the cellular sample;
determining, by a processor, for pixels of the first plurality of flow cell images, pixel intensities and a respective color purity of each of the pixel intensities; and
determining, by the processor and based on the pixel intensities and the respective color purity of the pixel intensities, a base calling template comprising base calling locations of the cellular sample at different axial locations along an axial axis, wherein the base calling template is configured for registering a second plurality of flow cell images of the cellular sample in one or more subsequent cycles of the one or more cycles.
2 . The method of claim 1 , wherein the first or second plurality of flow cell images of the cellular sample are generated at different axial locations along an axial axis.
3 . The method of claim 2 , wherein the axial axis is orthogonal to an image plane, and wherein a field of view of the first or second plurality of flow cell images is within the image plane.
4 . The method of claim 1 , further comprising:
registering, by the processor, the second plurality of flow cell images of the one or more subsequent cycles of the one or more cycles to the base calling template.
5 . The method of claim 1 , further comprising:
performing, by the processor, base calling of the second plurality of flow cell images based on the base calling locations of the base calling template.
6 . The method of claim 1 , wherein the first or second plurality of flow cell images are acquired at a plurality of predetermined axial locations and wherein at least some of the different axial locations of the base calling locations are different from any of the plurality of predetermined axial locations.
7 . The method of claim 1 , wherein each of the plurality of concatemer molecules immobilized on the support corresponds to a base calling location.
8 . The method of claim 1 , wherein the first or second plurality of flow cell images comprises an identical image resolution.
9 . The method of claim 1 , wherein the first plurality of flow cell images or the second plurality of flow cell images comprises optical signals emitted from nucleotide reagents bound to an unbalanced diversity of nucleotide bases of A, G, C and T/U among the plurality of concatemer molecules immobilized on the support.
10 . The method of claim 9 , wherein the unbalanced diversity of nucleotide bases of A, G, C and T/U among the plurality of concatemer molecules comprises:
a percentage of (1) a number of at least one type of nucleotide bases to (2) a total number of all four types of bases is less than 20%, 15%, 10%, or 5% in the one or more cycles.
11 . The method of claim 1 , wherein each of the pixel intensities comprises one or more subpixel intensities of a corresponding pixel.
12 . The method of claim 1 , wherein the respective color purity of each of the pixel intensities comprises:
one or more color purities corresponding to one or more subpixel intensities of a corresponding pixel.
13 . The method of claim 1 , wherein the respective color purity of each of the pixel intensities comprise: the respective color purity for one or more color channels.
14 . The method of claim 13 , wherein determining the respective color purity of each of the pixel intensities comprises: determining a ratio of (1) a signal corresponding to a particular type of nucleotide base to (2) a total amount of signal for other types of nucleotide bases.
15 . The method of claim 13 , wherein determining, for the pixels of the first plurality of flow cell images, the pixel intensities comprises:
determining, based on a comparison of a set of channel intensities at a corresponding pixel or a subpixel location, each channel intensity in the set of channel intensities corresponding to a respective different fluorescent wavelength.
16 . The method of claim 13 , wherein determining the base calling template based on the pixel intensities and the respective color purity of the pixel intensities comprises:
for each of the pixels or subpixels, determining whether or not the respective color purity is greater than a color purity of other pixels or subpixels within a threshold distance;
in response to determining that the respective color purity is greater than the color purity of the other pixels or subpixels within the threshold distance, adding a location of a corresponding pixel or subpixel to the base calling template; and
in response to determining that the respective color purity is no greater than the color purity of the other pixels or subpixels within the threshold distance, making no changes to the base calling template.
17 . The method of claim 16 , wherein the threshold distance is in three dimensions.
18 . The method of claim 1 , wherein the support comprises a flow cell.
19 . A system comprising:
one or more hardware processors;
one or more data storage devices storing instructions executable by the one or more hardware processors that, when executed, cause the one or more hardware processors to perform operations, the operations comprising:
generating, by a sequencing system, a first plurality of flow cell images of a cellular sample immobilized on a support by conducting one or more cycles of sequencing reactions, the cellular sample comprising a plurality of concatemer molecules therewithin, wherein a first concatemer molecule of the plurality of concatemer molecules corresponds to a first target RNA molecule of the cellular sample, and a second concatemer molecule of the plurality of concatemer molecules corresponds to a second target RNA molecule of the cellular sample;
determining, by a processor, for pixels of the first plurality of flow cell image, pixel intensities and a respective color purity of each of the pixel intensities; and
determining, by the processor and based on the pixel intensities and the respective color purity of the pixel intensities, a base calling template comprising base calling locations, wherein the base calling template is configured for registering a second plurality of flow cell images of the cellular sample in one or more subsequent cycles of the one or more cycles.
20 . One or more non-transitory computer storage media encoded with instructions executable by one or more hardware processors that, when executed, cause the one or more hardware processors to perform operations, the operations comprising:
generating, by a sequencing system, a first plurality of flow cell images of a cellular sample immobilized on a support by conducting one or more cycles of sequencing reactions, the cellular sample comprising a plurality of concatemer molecules therewithin, wherein a first concatemer molecule of the plurality of concatemer molecules corresponds to a first target RNA molecule of the cellular sample, and a second concatemer molecule of the plurality of concatemer molecules corresponds to a second target RNA molecule of the cellular sample;
determining, by a processor, for pixels of the first plurality of flow cell image, pixel intensities and a respective color purity of each of the pixel intensities; and
determining, by the processor and based on the pixel intensities and the respective color purity of the pixel intensities, a base calling template comprising base calling locations, wherein the base calling template is configured for registering a second plurality of flow cell images of the cellular sample in one or more subsequent cycles of the one or more cycles.