In a flow cytometer, a focused light source such as a laser beam or LED is arranged to shine on particles flowing within a fluid stream. The particles scatter the light and also cause fluorescent light of a different wavelength than the incident beam to be emitted.
The focused spot however is usually not a bright point that decays away monotonically in all directions. Rather, the beam from the light source such as a laser can have “lobes.” The far field of the beam is typically the Fourier transform of the laser aperture. Thus, for example a square aperture will give a far field which is a 2 dimensional Sync function as shown in
While the power in a laser minor side lobe 4 is typically 1% or less of the main lobe 3, a difficulty arises when very large particles 20 and very small particles 21 as shown in
This invention relates to the field of cytometry. It is sometimes necessary to detect particles of widely varying size, typically from 50 nm to 50 μm a dynamic range of 1000×. According to many embodiments, particle parameter measurement which may be recalculated instantaneously or at a later time, dependent on the determination as to if the particle was a large or small particle.
Many embodiments relate to an analog or digital deconvolution filter that can operate on a light intensity data pattern received from light emissions in a flow channel, so as to mitigate the effect of the laser lobes or other stray light in the flow channel caused by optical reflections or imperfections.
Many embodiments relate to a particle detection trigger mechanism that may be retriggered or restarted depending on measured values of data reaching fulfilling certain conditions requiring restart.
Many embodiments relate to a detection system in which laser light is directed at a flow path of particles of a flow cytometer. The laser light results in a light pattern having a plurality of lobes. A first signal is detected exceeding a first threshold. A second signal exceeding a second threshold is detected, wherein the second threshold is greater than the first threshold. Based on detecting the second trigger after a particular time point after detecting the first trigger, is determined that the first and second signals were created by a single relatively large particle.
In most cases, flow cytometers operate by flowing particles inside a small channel 17 formed in a middle of typically a quartz cuvette 12 as shown in
For simplistic analysis let us first consider a “Particle A” having an impulse as the spatial distribution of 0,1,0. It can be seen that if we convolve this particle pattern with a “Laser A” having an impulse spatial distribution of 0,1000,0 the resulting signal is 0,1000,0. Similarly if we convolve a “Laser B” with spatial distribution 0,1,0 with a “Particle B” of impulse spatial distribution 0,1000,0, the resultant signal is 0,1000,0. In this simple example, both cases cause an identical and indistinguishable signature of emitted light. In a more realistic case a particle that scatters more may be larger and have a wider spatial signature such as 0,0,1,10,100,1000,100,10,1,0,0.
A problem occurs when the laser spot spatial intensity variation in the cuvette is not an impulse, but instead is non-monotonically decaying undulating waveform as shown in
Minor laser lobes or their equivalent can also be created by near normal incidence of a laser beam 5 to an optical element 2 as shown in
Typically as shown in
Considering a flow cytometry system with a single laser and nearby side-lobes: First consider a “Particle A” with a spatial pattern that is an impulse of 0,1,0 convolved with a “Laser C” of spatial pattern 0,1,0,1000,0,1,0 (the ones on either side of the 1000 being the side lobes), the convolution yields 0,1,0,1000,0,1,0. Now consider a “Particle B” with a spatial pattern of 0,1000,0, convolved with the same laser pattern for “Laser C” the convolution yields 0,1000,0,1000000,0,1000,0. With a triggering threshold of 1000, one can detect the particle with the spatial pattern of 0,1,0 correctly, however, one will falsely detect the response to the lobe only for the particle with a spatial pattern of 0,1000,0, further, these particles may be “aborted” completely with a low trigger level because they may be classified as “doublets”—i.e. the presence of two particles in the laser beam simultaneously. In certain applications, it can be necessary to be able to distinguish between small particle entering the main lobe and a “larger” particle entering a side lobe. As used herein, the term “larger” not necessarily meaning physically larger but instead could be a particle of the same size but with more fluorescence or scattering.
Certainly a physically larger “Particle C” can have a spatial pattern like 0,0,1,10,100,1000,100,10,1,0,0 which indicates a relatively large physical size. Convolution of this particle with a “Laser C” with pattern 0,1,0,1000,0,1,0, yields 1,10,1100,11000,100101,1000000,100101,11000,1100,10,1 in which case the output signal is monotonically increasing and decreasing with a single peak even though the laser lobes are non-monotonic.
Depending on the spatial distribution of the laser and the spatial distribution of the particle, and the relative sizes of the light lobes compared to the particle size, it is possible for a single particle to create several blips that may be mistakenly interpreted as several different particles.
With knowledge of the laser lobe structure within the cuvette it is possible to mitigate the effect of the lobes. The essence of this method is “deconvolution” of the spatial distribution of the scattering or fluorescent emission of the particle with the spatial distribution of the laser light pattern.
In an analog cytometer system 500 as shown at
The logic 508 will classify any trigger event on the lower trigger, only when followed by a trigger on the higher trigger at a time separated by the distance of the laser lobes, to be a false first trigger of a large particle entering a side lobe, and requiring an immediate restart of the integrators based on the higher trigger level, since the higher trigger indicates a large particle in the main lobe of the laser, and not a side lobe. Comparatively, a smaller particle would not trigger on the high trigger on the side lobe, because the energy provided by the side lobe is not high enough to trigger the lower trigger, but only trigger at the main lobe.
Typically reset of integration logic 510 for determination of peak height and area is not required because the side lobe energy is so low that the main lobe can be integrated right on top of the side lobe integration, which started based on the lower trigger just previously, with less than 1% error. However, calculations such as width, which is determined by integrating a constant current onto a capacitor for the duration of the pulse, can require restart of the integration.
Alternatively, for a digital cytometer system where the entire data stream is digitized it is possible to process captured data through an inverse deconvolution filter by way of a processor shown at
At operation 604, relatively small and large particles are flowed through a flow path. Relatively small particles can be, for example, as small as 50 nm while relatively large particles can be as large as 50 μm particles, which is a dynamic range of 1000×. A detection system monitors for energy signals (fluorescence and/or scatter) resultant from particles interfacing with the light pattern of lobes. Generally, the center lobe provides greater energy than the shown side lobes. Hence, a relatively small particle will trigger a given threshold only when excited by the center lobe, while a relatively large particle will trigger the same threshold when passing through a side and center lobes, thus, causing potential measurement errors.
At operation 606, a first signal is detected that has an energy level exceeding a first threshold. The first threshold can be set low enough to detect a relatively small sized particle passing through the main lobe.
At operation 608, a second signal is detected (successively after the first signal) that has an energy level exceeding a second threshold that is higher than the first threshold. The second threshold can be set at a level to detect a relatively large particle passing through the main lobe. Hence, a small particle would not trigger the second threshold. If the second threshold was not triggered by the second signal, the particle would be processed as a small particle.
With both thresholds being triggered, it is still unknown whether the triggers occurred from a single large particle passing through multiple lobes, or from a large particle passing through a side lobe directly following a small particle passing through a main lobe. Accordingly, at operation 610, the detection system analyzes whether the second signal was received at a predetermined time after the first signal according to the time of particle flight between adjacent lobes. A large particle would trigger at a side lobe and main lobe between an expected time differential, since the flow rate and distance between lobes are both constant. While different particles would more than likely be separated at a distance that is not exactly coincident with the distance between lobes.
When it is determined that the first and second signals are resultant from a relatively large particle according to the time information, a requisite integration is performed to measure characteristics of the large particle at operation 612. When it is determined that the first and second signals are resultant from a relatively large particle directly following a relatively small particle, requisite integrations are performed to measure characteristics of the both particles at operation 614.
At operation 620, relatively small and large particles are flowed through a flow path. Relatively small particles can be, for example, as small as 50 nm while relatively large particles can be as large as 50 μm particles, which is a dynamic range of 1000×. A detection system monitors for energy signals (fluorescence and/or scatter) resultant from particles interfacing with the light pattern of lobes.
At operation 622, a time domain signal having a spatial emission pattern is detected from one or more particles interacting with the plurality of lobes. However, size of the one or more particles is not discernible in the emission pattern due to convolution between the spatial light pattern (plurality of lobes) of the laser and the spatial emission pattern of the particle.
At operation 624, the spatial emission pattern is deconvoluted (e.g., using a deconvolution filter as described above) from the time domain signal, based on dimensions of the plurality of lobes, to discern the size or spacing of the one or more particles.
It should be understood that the present invention as described above can be implemented in the form of control logic using computer software in a modular or integrated manner. Based on the disclosure and teachings provided herein, a person of ordinary skill in the art can know and appreciate other ways and/or methods to implement the present invention using hardware and a combination of hardware and software.
Any of the software components, user interfaces, or methods described in this application, may be implemented as software code to be executed by a processor using any suitable computer language such as, for example, Java, C++ or Perl using, for example, conventional or object-oriented techniques. The software code may be stored as a series of instructions, or commands on a computer readable medium, such as a random access memory (RAM), a read only memory (ROM), a magnetic medium such as a hard-drive or a floppy disk, or an optical medium such as a CD-ROM. Any such computer readable medium may reside on or within a single computational apparatus, and may be present on or within different computational apparatuses within a system or network.
This application is a divisional of U.S. Non-Provisional application Ser. No. 14/332,347, filed Jul. 15, 2014, and also claims the benefit of U.S. Provisional Application No. 61/957,887, filed Jul. 15, 2013, which is incorporated by reference.
Number | Date | Country | |
---|---|---|---|
61957887 | Jul 2013 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 14332347 | Jul 2014 | US |
Child | 15791199 | US |