WO2013121905A1 - 単一粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム - Google Patents

単一粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム Download PDF

Info

Publication number
WO2013121905A1
WO2013121905A1 PCT/JP2013/052446 JP2013052446W WO2013121905A1 WO 2013121905 A1 WO2013121905 A1 WO 2013121905A1 JP 2013052446 W JP2013052446 W JP 2013052446W WO 2013121905 A1 WO2013121905 A1 WO 2013121905A1
Authority
WO
WIPO (PCT)
Prior art keywords
light
light intensity
time
particle
particles
Prior art date
Application number
PCT/JP2013/052446
Other languages
English (en)
French (fr)
Inventor
田邊哲也
Original Assignee
オリンパス株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by オリンパス株式会社 filed Critical オリンパス株式会社
Priority to CN201380009765.5A priority Critical patent/CN104115001B/zh
Priority to EP13748583.5A priority patent/EP2816344A4/en
Priority to JP2014500168A priority patent/JP5996625B2/ja
Publication of WO2013121905A1 publication Critical patent/WO2013121905A1/ja
Priority to US14/451,021 priority patent/US9494779B2/en

Links

Images

Classifications

    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00—Microscopes
    • G02B21/0004—Microscopes specially adapted for specific applications
    • G02B21/002—Scanning microscopes
    • G02B21/0024—Confocal scanning microscopes (CSOMs) or confocal "macroscopes"; Accessories which are not restricted to use with CSOMs, e.g. sample holders
    • G02B21/008—Details of detection or image processing, including general computer control
    • G02B21/0084—Details of detection or image processing, including general computer control time-scale detection, e.g. strobed, ultra-fast, heterodyne detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Electro-optical investigation, e.g. flow cytometers
    • G01N15/1429—Electro-optical investigation, e.g. flow cytometers using an analyser being characterised by its signal processing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Electro-optical investigation, e.g. flow cytometers
    • G01N15/1434—Electro-optical investigation, e.g. flow cytometers using an analyser being characterised by its optical arrangement
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/6408—Fluorescence; Phosphorescence with measurement of decay time, time resolved fluorescence
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/645—Specially adapted constructive features of fluorimeters
    • G01N21/6456—Spatial resolved fluorescence measurements; Imaging
    • G01N21/6458—Fluorescence microscopy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/6486—Measuring fluorescence of biological material, e.g. DNA, RNA, cells
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00—Microscopes
    • G02B21/0004—Microscopes specially adapted for specific applications
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00—Microscopes
    • G02B21/0004—Microscopes specially adapted for specific applications
    • G02B21/002—Scanning microscopes
    • G02B21/0024—Confocal scanning microscopes (CSOMs) or confocal "macroscopes"; Accessories which are not restricted to use with CSOMs, e.g. sample holders
    • G02B21/0052—Optical details of the image generation
    • G02B21/0076—Optical details of the image generation arrangements using fluorescence or luminescence
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84—Systems specially adapted for particular applications
    • G01N2021/8405—Application to two-phase or mixed materials, e.g. gas dissolved in liquids
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/12—Circuits of general importance; Signal processing
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00—Microscopes
    • G02B21/16—Microscopes adapted for ultraviolet illumination ; Fluorescence microscopes

Definitions

  • the present invention uses an optical system capable of detecting light from a minute region in a solution, such as an optical system of a confocal microscope or a multiphoton microscope, and uses atoms, molecules or aggregates thereof dispersed or dissolved in a solution ( These are hereinafter referred to as “particles”), for example, biomolecules such as proteins, peptides, nucleic acids, lipids, sugar chains, amino acids or aggregates thereof, particulate objects such as viruses and cells, or non-
  • the present invention relates to an optical analysis technique capable of detecting biological particles and obtaining information useful for analyzing or analyzing their states (interaction, binding / dissociation state, etc.).
  • the detected light may be fluorescence, phosphorescence, chemiluminescence, bioluminescence, scattered light, or the like.
  • the particles that emit light may be either particles that emit light themselves, or particles to which an arbitrary luminescent label or luminescent probe is added.
  • Such optical analysis techniques include, for example, fluorescence correlation spectroscopy (FCS; see, for example, Patent Literature 1-3 and Non-Patent Literature 1-3), and fluorescence intensity distribution analysis (Fluorescence-Intensity Distribution Analysis: FIDA, for example, Patent Document 4, Non-Patent Document 4) and Photon Counting Histogram (PCH, for example, Patent Document 5) are known.
  • Patent Documents 6 to 8 propose a method for detecting a fluorescent substance based on the passage of time of a fluorescence signal of a sample solution measured using an optical system of a confocal microscope.
  • the applicant of the present application described in Patent Documents 9 to 11 is an optical analysis technique using an optical system capable of detecting light from a minute region in a solution, such as an optical system of a confocal microscope or a multiphoton microscope.
  • a new optical analysis technique based on a principle different from optical analysis techniques such as FCS and FIDA has been proposed.
  • optical analysis techniques such as FCS and FIDA
  • light detection region light from fluorescent molecules floating in a minute region (hereinafter referred to as “light detection region”) that is a light detection region in a sample solution is continuously emitted.
  • the light intensity data obtained by the measurement is analyzed by a calculation process for calculating statistical fluctuations, and the concentration and / or other characteristics of the fluorescent molecules are detected.
  • the position of the light detection region is moved in the sample solution, that is, the sample solution is scanned by the light detection region.
  • the light detection region includes light emitting particles that are dispersed in the sample solution and move randomly, the light emitted from the light emitting particles is individually detected, so that one light emitting particle in the sample solution is detected. By detecting each of them, it is possible to obtain information on the counting of the luminescent particles and the concentration or number density of the luminescent particles in the sample solution.
  • the sample required for the measurement is a very small amount (for example, about several tens of ⁇ L) as in the optical analysis techniques such as FCS and FIDA.
  • the measurement time is short (measurement of time in the order of seconds is repeated several times in one measurement), and the concentration of particles to be observed is the same as that of optical analysis technology such as FCS and FIDA. It is possible to detect the presence of luminescent particles in a sample solution lower than a well measurable level (about 1 nM) and quantitatively detect its concentration, number density or other characteristics.
  • the “scanning molecule counting method” is used for analysis of rare or expensive samples often used in the field of medical and biological research and development, clinical diagnosis of diseases, screening of bioactive substances, etc.
  • the number of specimens is large, it is possible to perform experiments or tests at a lower cost or more quickly than conventional biochemical methods, and at a lower concentration so that FCS, FIDA, etc. cannot be performed well. It is expected to be a powerful tool capable of detecting the concentration and / or properties of
  • the observation target particle when the observation target particle is a luminescent particle (as will be described later, the observation target particle may be a particle that does not emit light in the detection wavelength band (non-luminescent particle)), time-series light intensity data Of the increase in bell-shaped or pulse-shaped light intensity observed in, conforms to the conditions of the light intensity increase profile (peak intensity, full width at half maximum, etc.) assumed when the luminescent particles pass through the light detection region A change in light intensity is detected as a signal representing the presence of particles. An increase in light intensity that does not meet the conditions of the increase in light intensity of the luminescent particles is determined as noise.
  • the time-series light intensity data is discrete photon counting data (see FIG. 2A).
  • the time-series light intensity data is preferably smoothed with respect to time, and the change in the light intensity of the particles in the smoothed time-series light intensity data. A bell-shaped signal that meets the profile conditions is detected.
  • the signal when the change in the intensity of the particle signal is small, the signal is buried in noise at the stage of the photon counting data, and even after the smoothing process, the change in the light intensity of the particle. It does not exhibit the characteristics of a bell-shaped profile and is not detected as a particle signal.
  • the magnitude of the signal corresponding to the presence of the particle changes depending on the passage route of the particle in the photodetection region, and the intensity change of the signal of the particle passing through the periphery of the photodetection region is small. Accordingly, the signal of particles passing through the periphery of many weak light detection regions is overlooked, or the noise signal is misjudged as a signal of weak particles, which is the sensitivity of the scanning molecule counting method. Or it was one of the obstacles to improving accuracy.
  • a main object of the present invention is to provide a novel particle signal detection means or method that enables further improvement in sensitivity and / or accuracy in the above-described scanning molecule counting method.
  • Another object of the present invention is to provide a novel particle signal detection means or method that enables more accurate discrimination between weak particle signals and noise signals in the scanning molecule counting method. Is to provide.
  • the inventor of the present invention in the time-series light intensity data of the scanning molecule counting method, in the portion corresponding to the particle signal and the noise signal portion, the time change of the light intensity (photon sequence) ) was found to be different. That is, in the case of a noise signal, photons are always detected randomly, whereas in the case of a particle signal, photons are detected in a concentrated manner over time.
  • the detection accuracy of the particle signal in the scanning molecule counting method can be improved by a novel algorithm capable of detecting the difference in the occurrence pattern of the light intensity with time.
  • the occurrence pattern of light intensity with time changes also changes depending on the characteristics of the particles (polarization characteristics and emission wavelength characteristics). Tried.
  • the above-described problem is an optical analyzer that detects single particles that are dispersed in a sample solution and move randomly using an optical system of a confocal microscope or a multiphoton microscope,
  • a light detection region moving unit that moves the position of the light detection region of the optical system of the microscope in the sample solution, a light detection unit that detects light from the light detection region, and a light detection region in the sample solution.
  • a signal processing unit that individually detects, the signal processing unit within the light detection region of the time change of the light intensity value in each of the analysis windows set in time series on the time series light intensity data Assuming the first state where there is no single particle One occurrence probability and a second occurrence probability assuming a second state in which a single particle exists in the light detection region, and based on the first and second occurrence probabilities This is accomplished by an apparatus characterized by detecting a signal representative of the presence of each single particle on time-series light intensity data.
  • the term “single particles dispersed in the sample solution and moving randomly” refers to particles such as atoms, molecules, or aggregates thereof dispersed or dissolved in the sample solution. Any particle may be used as long as it is not fixed to the substrate or the like and freely moves in the solution in Brownian motion.
  • the single particle to be observed may be a particle that emits light (light-emitting particle) or a particle that does not emit light (in the detection wavelength band) (non-light-emitting particle).
  • the “signal indicating the presence of each single particle” is the time series light intensity data corresponding to the light emitted when the luminescent particle passes through the light detection region.
  • the luminescent particles are typically fluorescent particles, but may be particles that emit light by phosphorescence, chemiluminescence, bioluminescence, light scattering, or the like.
  • the observation target particle may be a non-light emitting particle, in which case the light from the light detection region includes significant background light, and the “signal indicating the presence of each single particle” is the background light.
  • Light intensity from the light source (reversed scanning molecule counting method).
  • the “light detection region” of the optical system of the confocal microscope or the multiphoton microscope is a minute region in which light is detected in those microscopes, and when illumination light is given from the objective lens, Corresponds to the area where the illumination light is collected (in a confocal microscope, it is determined in particular by the positional relationship between the objective lens and the pinhole.
  • the light detection unit detects light from the light detection region by photon counting that counts the number of photons that arrive every predetermined measurement unit time (bin time), and in this case, the light intensity in time series The data becomes time-series photon count data.
  • the position of the light detection region in the sample solution is basically the same as the “scanning molecule counting method” described in Patent Documents 9 to 11. While moving the sample, that is, while scanning the sample solution with the light detection region, the light detection and the generation of the time series light intensity data representing the time series light intensity value are sequentially performed. A signal representing the presence of a single particle is detected on the data.
  • the detection of the signal of the single particle is not executed based on whether or not the increase / decrease amount of the light intensity value exceeds a predetermined amount,
  • the pattern of the time change of the light intensity value on the time-series light intensity data is likely to occur when there is no single particle in the light detection region or when there is a single particle in the light detection region. Is evaluated in time series. That is, as already mentioned, in the time-series light intensity data, there is a difference in the temporal change pattern of the light intensity between the part corresponding to the particle signal and the noise signal part.
  • the temporal change pattern of the light intensity value in a certain part above is a pattern that is likely to occur when there is no single particle in the light detection region, it is assumed that the part is only a noise signal part. If the pattern is easy to be generated when there is a single particle in the light detection region, it can be determined that the portion corresponds to the signal of the particle. Then, whether the pattern of the time change of a certain light intensity value is a pattern that is likely to occur when there is a single particle in the light detection region or not, the pattern of the time change of the light intensity value is This can be determined by the probability of occurrence in the case where there is a single particle in the light detection region and in the case where there is no single particle.
  • the apparatus of the present invention in order to detect a portion where particles exist in the light detection region on the time-series light intensity data, that is, a signal indicating the presence of a single particle, Assuming a first state in which no single particle exists in the light detection region of the time change of the light intensity value in each analysis window set in time series on the light intensity data The first occurrence probability and the second occurrence probability when a second state in which a single particle exists in the light detection region is assumed are calculated.
  • the “analysis window” is a region of any time width on the time series light intensity data, and the “analysis window” is set sequentially or in time series on the time series light intensity data. Is done.
  • the “first occurrence probability” is a probability that a time change of the light intensity measured in the “analysis window” occurs when no single particle exists in the light detection region
  • the “second occurrence probability” is a probability that a time change of the light intensity measured in the “analysis window” occurs when a single particle is present in the light detection region. If there is a single particle in the light detection region, the “second occurrence probability” will increase relative to the “first occurrence probability”.
  • the “probability” and the “second occurrence probability” the time domain in which the single particle exists in the light detection region on the time-series light intensity data is determined. It is possible to detect a signal indicating the presence of each.
  • the width of the analysis window may be set to an arbitrary width equal to or longer than the time required for a single particle to pass through the light detection region, for example.
  • the analysis window may be set for each unit time or for every predetermined time interval on the time-series light intensity data (in this case, the analysis windows set in time series overlap each other). Alternatively, it may be set by dividing the time-series light intensity data by a predetermined time width.
  • the “unit time” is a time period that gives one light intensity value in optical measurement, and in the case of photon counting, it may be one or more bin time width times. .
  • the first and second occurrence probabilities are the actual values from the average pattern of the temporal change of the light intensity value assuming the first and second states in the analysis window, respectively. It can be determined based on the deviation of the temporal change pattern of the measured light intensity value. The smaller the deviation of the actual measurement pattern from the average pattern, the higher the probability that the actual measurement pattern will occur. Accordingly, the first and second occurrence probabilities are respectively the detected light intensity value per unit time in the analysis window and the expected value per unit time assuming the first and second states.
  • the unit time occurrence probability of the light intensity value per unit time is calculated as following the Poisson distribution in which the light intensity value per unit time has an expected value per unit time,
  • the first and second occurrence probabilities may be calculated using the corresponding unit time occurrence probabilities.
  • the ratio or odds ratio of the first and second occurrence probabilities is calculated in time series, and the time when a single particle is present in the light detection region is determined based on the ratio of the probability or the odds ratio. Also good.
  • the above scanning molecule counting method it is possible to separately detect at least two different components of light from the light detection region and generate time-series light intensity data of each component.
  • the component to be detected so that the data of the component to be detected reflects any characteristic of the single particle. For example, if components having different polarization directions are selected as a plurality of components to be detected, the polarization characteristics of particles are reflected in the time-series light intensity data of the plurality of components. Further, if light components having different wavelength bands are selected as the plurality of components to be detected, the emission wavelength characteristics of the particles are reflected in the time-series light intensity data of the plurality of components. The above-described characteristics of the observation target particles are also reflected in the temporal change pattern of the light intensity value in the time-series light intensity data of a plurality of components.
  • the light detection unit separately detects at least two different components of light from the light detection region, and the signal processing unit detects the time-series light of each of the components.
  • the intensity data may be generated, and the signal processing unit may be configured to calculate the first and second occurrence probabilities of each of the components.
  • the detection component by appropriately selecting the detection component, the predetermined characteristic value of the single particle to be observed is reflected in the second occurrence probability of each detected light component, that is,
  • the second occurrence probability can be a function of a predetermined characteristic value of the single particle.
  • the second occurrence probability that is a function of the predetermined characteristic value has a predetermined characteristic value in a relatively high region. It can be determined that the particles are present (if there are particles that do not have the predetermined characteristic value, the second occurrence probability that is a function of the predetermined characteristic value is low). That is, by using the occurrence probability of the first and second occurrence probabilities for each component, not only the presence or absence of particles in the light detection region, but also in the light detection region on the time-series light intensity data. The time during which particles having a predetermined characteristic value exist is determined, and the presence of particles having a predetermined characteristic value can be detected.
  • a plurality of types of singles having predetermined characteristic values different from each other as a single particle are used. if it contains particles, for each type of single particles, there was time of the particles in the time series of the light intensity data on is can be determined.
  • the second occurrence probability of each component that is a function of their different predetermined characteristic values is calculated (since no particles exist in the first state, the first The characteristic value of the particle is not reflected in the occurrence probability of the particle.), See the first occurrence probability of each of the plurality of components and the second occurrence probability of each of the components of the plurality of types of single particles Then, it can be estimated that a kind of particle that gives a second occurrence probability having a relatively high value exists in the light detection region.
  • the above-described apparatus of the present invention is a function of different predetermined characteristic values for each single particle type.
  • the characteristic value of the particle to be reflected in the second occurrence probability is an index value indicating the polarization characteristic of a single particle, such as fluorescence anisotropy or the emission wavelength of a single particle. It may be an index value representing a characteristic, for example, a ratio of emission intensities in different emission wavelength bands.
  • the moving speed of the position of the light detection region in the sample solution in the apparatus of the present invention described above can be appropriately changed based on the characteristics of the single particle to be observed or the number density or concentration in the sample solution. It's okay.
  • the moving speed of the light detection region is increased, the amount of light obtained from one luminescent particle is reduced when the single particle is a luminescent particle, and when the single particle is a non-luminescent particle.
  • the amount of reduction in the light intensity value due to the presence of one non-light emitting particle is reduced. Therefore, it is preferable that the moving speed of the light detection region can be appropriately changed so that a change in light intensity value due to a single particle can be measured with high accuracy or sensitivity.
  • the moving speed of the position of the light detection region in the sample solution is preferably set higher than the diffusion moving speed of the single particle to be detected (the average moving speed of the particles due to Brownian motion).
  • the temporal change pattern of the detected light intensity value is assumed to be when there are no particles in the light detection region (first and second states).
  • the second movement is considered in consideration of the movement of the particles by Brownian motion when the particles pass through the light detection region. The calculation of the occurrence probability is complicated.
  • the movement speed of the light detection area is higher than the diffusion movement speed of the single particle to be detected so that the effect of the Brownian motion of the particles while the particles pass through the light detection area can be ignored. It is preferable that the setting is fast. Since the diffusion movement speed varies depending on the single particle, as described above, the apparatus of the present invention is configured such that the movement speed of the light detection region can be appropriately changed according to the characteristics (particularly, the diffusion coefficient). It is preferable that
  • the movement of the position of the light detection region in the sample solution may be performed by an arbitrary method.
  • the position of the light detection region may be changed by changing the optical path of the optical system of the microscope by using a galvanometer mirror adopted in a laser scanning optical microscope, or
  • the position of the sample solution may be moved (for example, by moving a microscope stage) to move the position of the light detection region in the sample solution.
  • the movement trajectory of the position of the light detection region may be arbitrarily set, and may be selected from, for example, a circle, an ellipse, a rectangle, a straight line, and a curve.
  • the number of single particles included in the light detection region may be counted by counting the number of signals (particle counting).
  • information on the number density or concentration of the identified single particles in the sample solution can be obtained by combining the number of detected single particles and the amount of movement of the position of the light detection region.
  • the number density or concentration ratio of a plurality of sample solutions, or the relative number density or concentration ratio with respect to the standard sample solution serving as a reference for the concentration or number density is calculated, or An absolute number density value or concentration value may be determined using a relative number density or concentration ratio relative to a standard sample solution that is a reference for concentration or number density.
  • the total volume of the movement locus of the position of the light detection region is specified by any method, for example, by moving the position of the light detection region at a predetermined speed, the number density or concentration of single particles Can be calculated specifically.
  • photodetection is performed by detecting light while moving the position of the light detection region in the sample solution and individually detecting signals from a single particle.
  • the processing of the optical analysis technique for determining the presence or absence of a single particle with reference to the probability (first and second occurrence probability) of occurrence of a temporal change pattern of the light intensity value can also be realized by a general-purpose computer. is there.
  • a computer program for optical analysis for detecting single particles dispersed and moving randomly in a sample solution using an optical system of a confocal microscope or a multiphoton microscope
  • the procedure for calculating the probability and the first and second occurrence probabilities Computer program, characterized in that to execute the steps of detecting a signal indicative of the presence of each single particle in time series light intensity data on the computer is provided based. Also in this case, like the apparatus of the present invention, typically, in the procedure of detecting light from the light detection region and generating time-series light intensity data, the measurement unit time (bin time) Light from the photodetection region is detected by photon counting that counts the number of photons that arrive every time, and time-series light intensity data is time-series photon count data.
  • the “analysis window” may be set similarly to the case of the apparatus of the present invention.
  • the particle to be observed is a luminescent particle
  • the “temporary increase in the light intensity value” in the time-series light intensity data becomes the “signal indicating the presence of each single particle”.
  • the particles to be observed are non-luminous particles (in the detection wavelength band)
  • the light from the light detection region includes significant background light
  • the light intensity from the background light is temporarily The reduction becomes a “signal that represents the presence of each single particle”.
  • the computer program is stored in a computer-readable storage medium and provided. The computer reads the program stored in the storage medium and execute processing and operation processing of the information, to realize the above procedure.
  • the computer-readable recording medium may be a magnetic disk, a magneto-optical disk, a CD-ROM, a DVD-ROM, a semiconductor memory, or the like.
  • the above-described program may be distributed to a computer via a communication line, and the computer that has received this distribution may execute the program.
  • the first and second occurrence probabilities are assumed when the light intensity value per unit time in the analysis window and the first and second states are assumed, respectively. More specifically, the unit of the light intensity value per unit time is calculated based on the Poisson distribution in which the light intensity value per unit time has the expected value per unit time. The time occurrence probability may be calculated, and the first and second occurrence probabilities may be calculated using the corresponding unit time occurrence probabilities. Then, it is determined that a single particle is present in the light detection region at the time of the analysis window where the second occurrence probability is greater than the first occurrence probability, or the second occurrence probability and the first occurrence probability Based on the odds ratio, the time domain in which a single particle was present in the light detection region may be determined.
  • the above computer program also separately detects at least two different components of light from the light detection region, and generates time-series light intensity data for each of the components, and for each of the components,
  • the second occurrence probability which is a function of the first occurrence probability and a predetermined characteristic value of the single particle to be observed, is calculated, and based on the occurrence probability of the first and second occurrence probabilities for each component
  • a signal representing the presence of a single particle having a predetermined characteristic value on the series of light intensity data may be detected.
  • each of the single particles has a predetermined characteristic value different from each other.
  • a second probability of occurrence of each component that is a function is calculated and based on the first probability of occurrence of each of the components and the second probability of occurrence of each component of each of the plurality of types of single particles.
  • a signal representing the presence of a single particle may be detected for each type of single particle on the series of light intensity data.
  • the moving speed of the position of the light detection region in the sample solution may be appropriately changed based on the characteristics of the single particle, the number density or concentration in the sample solution, and preferably in the sample solution.
  • the moving speed of the position of the light detection region at is set higher than the diffusion moving speed of the single particle to be detected.
  • the position of the light detection region in the sample solution may be moved by an arbitrary method.
  • the position of the light detection region is changed by changing the optical path of the optical system of the microscope or moving the position of the sample solution. May be changed.
  • the movement trajectory of the position of the light detection region may be arbitrarily set, and may be selected from, for example, a circle, an ellipse, a rectangle, a straight line, and a curve.
  • a procedure for counting the number of single particles detected during movement of the position of the light detection region by counting the number of signals from single particles detected individually and A procedure for determining the number density or concentration of single particles in the sample solution based on the number of single particles detected may be included.
  • an optical analysis method for detecting light of individual luminescent particles while moving the position of a light detection region in a sample solution A novel method for determining the presence / absence of a single particle is realized by referring to the probability (first and second occurrence probability) of occurrence of a temporal change pattern of the intensity value.
  • an optical analysis method for detecting single particles dispersed and moving randomly in a sample solution using an optical system of a confocal microscope or a multiphoton microscope, In the process of moving the position of the light detection area of the optical system of the microscope, and measuring the light intensity from the light detection area while moving the position of the light detection area in the sample solution, the light intensity data
  • the generation process and the time variation of the light intensity value in each of the analysis windows set in time series on the time series light intensity data are the first in which no single particle exists in the light detection region.
  • the “analysis window” may be set similarly to the case of the apparatus of the present invention.
  • the particle to be observed is a luminescent particle
  • the “temporary increase in the light intensity value” in the time-series light intensity data becomes the “signal indicating the presence of each single particle”.
  • the particles to be observed are non-luminous particles (in the detection wavelength band)
  • the light from the light detection region includes significant background light
  • the light intensity from the background light is temporarily The reduction becomes a “signal that represents the presence of each single particle”.
  • the first and second occurrence probabilities are assumed when the light intensity value per unit time in the analysis window and the first and second states are assumed, respectively. More specifically, the unit of the light intensity value per unit time is calculated based on the Poisson distribution in which the light intensity value per unit time has the expected value per unit time. The time occurrence probability may be calculated, and the first and second occurrence probabilities may be calculated using the corresponding unit time occurrence probabilities. Then, it is determined that a single particle is present in the light detection region at the time of the analysis window where the second occurrence probability is greater than the first occurrence probability, or the second occurrence probability and the first occurrence probability Based on the odds ratio, the time domain in which a single particle was present in the light detection region may be determined.
  • At least two different components of the light from the light detection region are separately detected, and time-series light intensity data of each of the components is generated.
  • a second occurrence probability that is a function of one occurrence probability and a predetermined characteristic value of a single particle to be observed is calculated, and a time series based on the occurrence probability of the first and second occurrence probabilities for each component
  • a signal representing the presence of a single particle having a predetermined characteristic value on the light intensity data may be detected.
  • a second occurrence probability of each of the components is calculated, and a time series based on the first occurrence probability of each of the components and the second occurrence probability of each of each component of the plurality of types of single particles For each single particle type, a signal indicating the presence of the single particle may be detected on the light intensity data.
  • the moving speed of the position of the light detection region in the sample solution may be appropriately changed based on the characteristics of the single particle, the number density or concentration in the sample solution, and preferably in the sample solution.
  • the moving speed of the position of the light detection region at is set higher than the diffusion moving speed of the single particle to be detected.
  • the position of the light detection region in the sample solution may be moved by an arbitrary method.
  • the position of the light detection region is changed by changing the optical path of the optical system of the microscope or moving the position of the sample solution. May be changed.
  • the movement trajectory of the position of the light detection region may be arbitrarily set, and may be selected from, for example, a circle, an ellipse, a rectangle, a straight line, and a curve.
  • the optical analysis technique of the present invention described above is typically a biological molecule such as a protein, peptide, nucleic acid, lipid, sugar chain, amino acid or aggregate thereof, or a particulate biological object such as a virus or cell. It is used for analysis or analysis of the state of matter in solution, but it may also be used for analysis or analysis of the state of non-biological particles (eg, atoms, molecules, micelles, metal colloids, etc.) in solution. It should be understood that such cases are also within the scope of the present invention.
  • a biological molecule such as a protein, peptide, nucleic acid, lipid, sugar chain, amino acid or aggregate thereof, or a particulate biological object such as a virus or cell. It is used for analysis or analysis of the state of matter in solution, but it may also be used for analysis or analysis of the state of non-biological particles (eg, atoms, molecules, micelles, metal colloids, etc.) in solution. It should be understood that such cases are also within the scope of
  • the scanning molecule counting method in the scanning molecule counting method, a state in which a change pattern of the light intensity value is likely to occur rather than simply referring to the increase / decrease of the light intensity value on the time-series light intensity data.
  • the change in the light intensity value due to the single particle is relatively small, and even if it is difficult to determine whether the signal due to the particle or the noise signal is difficult only with the absolute value of the light intensity value change, It is expected that a signal and a noise signal can be distinguished.
  • the concentration range of single particles in a sample solution that can be detected by the scanning molecule counting method is expected to expand to a lower concentration side.
  • the single particle to be observed and other particles Identification is also possible, and improvement in measurement accuracy or sensitivity is expected even in sample solutions in which impurities and the like are mixed.
  • FIG. 1A is a schematic diagram of the internal structure of an optical analyzer that performs the scanning molecule counting method according to the present invention.
  • FIG. 1B is a schematic diagram of a confocal volume (observation region of a confocal microscope).
  • FIG. 1C is a schematic diagram of a mechanism for changing the direction of the mirror 7 to move the position of the light detection region in the sample solution.
  • FIG. 1D is a schematic diagram of a mechanism for moving the position of the photodetection region in the sample solution by moving the horizontal position of the microplate.
  • FIG. 2A and 2B are a schematic diagram for explaining the principle of detecting the light of the luminescent particles in the scanning molecule counting method to which the present invention is applied, and a schematic diagram of the time change of the measured light intensity, respectively.
  • FIG. 2 (C) and 2 (D) are schematic diagrams for explaining the principle of detecting the presence of a single particle that does not emit light in the detection light wavelength band in the inverted scanning molecule counting method to which the present invention is applied, respectively. It is a schematic diagram of the time change of the light intensity measured.
  • FIG. 3A shows a schematic example of time-series light intensity data.
  • (Left) When luminescent particles with high luminance exist.
  • (Middle) When luminescent particles with low luminance exist.
  • FIG. 3B is a diagram for explaining an analysis window set on the time-series light intensity data in the present invention.
  • FIG. 3C shows a mode of movement of particles when the light-emitting particles cross the light detection region by moving the position of the light detection region in the sample solution at a speed faster than the diffusion movement speed of the light-emitting particles. It is a model figure.
  • FIG. 3D is a schematic diagram of the intensity distribution in the radius r direction of the light detection region of the light emitted and detected from the luminescent particles in the light detection region.
  • 4A is time-series light intensity data schematically represented for the purpose of explanation, and the upper row shows a state where the luminescent particles are present (left: second state) and a state where the luminescent particles are present. Represents the expected value of the light intensity value on the time-series light intensity data when each of the non-existent states (right: first state) is assumed, and the middle row represents the actually measured light intensity value, The lower row shows the occurrence probability of the measured light intensity value when assuming a state where the luminescent particles are present (left) and a state where the luminescent particles are not present (right).
  • FIG. 4B is a model diagram of the mutually polarized components of the light emitted from the luminescent particles (upper), the average temporal change (lower left) of the light intensity value of each component, and the average of the total light intensity value Time change (lower right).
  • FIG. 4C shows the emission wavelength spectra of the light emitting particles D1 and D2 having different emission wavelength characteristics. When light components in different wavelength bands are separately detected, detection is performed in each wavelength band. The light intensity value of the component from each luminescent particle is shown. The area of the shaded portion corresponds to the light intensity value of each component. It is the figure which represented the process sequence of the scanning molecule counting method performed according to this invention in the format of the flowchart.
  • FIG. 4C shows the emission wavelength spectra of the light emitting particles D1 and D2 having different emission wavelength characteristics. When light components in different wavelength bands are separately detected, detection is performed in each wavelength band. The light intensity value of the component from each luminescent particle is shown. The area of the shaded portion corresponds
  • FIG. 6 shows a part (A) of time-series light intensity data (photon count data) obtained according to the scanning molecule counting method according to the present invention and a smoothing obtained by smoothing the time-series light intensity data of (A).
  • A time-series light intensity data
  • B time-series light intensity data
  • A time-series light intensity data
  • C time-series light intensity data
  • FIG. 7 is a view similar to FIG. 6 and shows another part of the time-series light intensity data.
  • FIG. 8 is a diagram showing the relationship between the concentration of luminescent particles (ATTO647N) in the sample solution and the number of detected particles when the luminescent particles are detected and counted according to the scanning molecule counting method of the present invention.
  • FIG. 9A shows a part (upper) of time-series light intensity data (photon count data) obtained according to the inverted scanning molecule counting method according to the present invention and the time-series light intensity data according to the teaching of the present invention.
  • the odds ratio (lower stage) of the occurrence probability (second and first occurrence probabilities) calculated on the assumption that each of the particles is present and not present is shown.
  • FIG. 9B is an enlarged view of a part of FIG.
  • FIG. 10 is an enlarged view of FIG. 9B.
  • 11A is a diagram showing the relationship between the concentration of particles in a sample solution and the number of detected particles when particles are detected and counted according to the inverted scanning molecule counting method according to the present invention.
  • 11 (B) is a diagram showing the relationship between the particle concentration in the sample solution and the number of particles detected by smoothing and fitting of bell-shaped function in the same time-series light intensity data as (A). is there. In the figure, error bars are standard deviations.
  • FIG. 12 shows the fluorescence polarization anisotropy of luminescent particles in time-series light intensity data obtained by separately detecting components of different polarization directions in accordance with the scanning molecule counting method according to the present invention.
  • the left is the discrimination result of the type of luminescent particles detected in the solution containing the plasmid
  • the right is the discrimination result of the type of luminescent particles detected in the solution containing the fluorescent dye (TAMRA). It is.
  • optical analysis apparatus 1 that realizes the optical analysis technique according to the present invention is capable of executing FCS, FIDA, and the like as schematically illustrated in FIG.
  • the apparatus may be a combination of an optical system of a confocal microscope and a photodetector.
  • optical analysis apparatus 1 includes optical systems 2 to 17 and a computer 18 for controlling the operation of each part of the optical system and acquiring and analyzing data.
  • the optical system of the optical analyzer 1 may be the same as the optical system of a normal confocal microscope, in which the laser light (Ex) emitted from the light source 2 and propagated through the single mode fiber 3 is a fiber.
  • the light is emitted as a divergent light at an angle determined by a specific NA at the outgoing end of the light, becomes parallel light by the collimator 4, is reflected by the dichroic mirror 5, the reflection mirrors 6, 7, and the objective lens 8. Is incident on.
  • a microplate 9 in which a sample container or well 10 into which a sample solution of 1 to several tens of ⁇ L is dispensed is typically arranged is emitted from the objective lens 8.
  • the laser light is focused in the sample solution in the sample container or well 10 to form a region with high light intensity (excitation region).
  • the luminescent particle typically a particle to which a luminescent label such as a fluorescent particle or a fluorescent dye is added is dispersed or dissolved in the sample solution.
  • the luminescent particles enter the excitation region, the luminescent particles are excited and light is emitted.
  • the sample solution typically contains dispersed or dissolved particles that do not emit light in the detection wavelength band and any luminescent substance that produces background light.
  • the luminescent material is excited to emit substantially constant light, which becomes background light and particles that do not emit light in the detection wavelength band.
  • the background light is reduced.
  • the light (Em) emitted from the excitation region passes through the objective lens 8 and the dichroic mirror 5, is reflected by the mirror 11, collected by the condenser lens 12, passes through the pinhole 13, and passes through the barrier filter. 14 (here, only a light component in a specific wavelength band is selected), introduced into the multimode fiber 15, reaches the photodetector 16, and is converted into a time-series electrical signal. Thereafter, it is input to the computer 18 and processing for optical analysis is performed in a manner described later.
  • the pinhole 13 is disposed at a position conjugate with the focal position of the objective lens 8, and as a result, as shown in FIG.
  • the focal region of the laser beam illustrated in FIG. 1B is usually a light detection region in the present optical analyzer having an effective volume of about 1 to 10 fL (typically, the light intensity is in the region).
  • the photodetector 16 since light from one luminescent particle, for example, weak light from one fluorescent dye molecule or background light reduction due to the presence of non-luminescent particles is detected, the photodetector 16 is suitable. In this case, an ultrasensitive photodetector that can be used for photon counting is used. When the light detection is based on photon counting, the light intensity is measured in a mode in which the number of photons arriving at the light detector is measured every measurement unit time (BIN TIME) sequentially over a predetermined time.
  • BIN TIME measurement unit time
  • the time-series light intensity data is time-series photon count data.
  • a stage position changing device 17a for moving the horizontal position of the microplate 9 may be provided on a microscope stage (not shown) in order to change the well 10 to be observed.
  • the operation of the stage position changing device 17a may be controlled by the computer 18. With this configuration, it is possible to achieve quick measurement even when there are a plurality of specimens.
  • a mechanism for scanning the sample solution by the light detection region that is, for moving the position of the focal region, that is, the light detection region in the sample solution.
  • a mechanism for moving the position of the light detection region for example, a mirror deflector 17 that changes the direction of the reflection mirror 7 may be employed as schematically illustrated in FIG. A method of moving the absolute position of the light detection area).
  • Such mirror deflector 17 may be similar to the galvanometer mirror apparatus installed in a conventional laser scanning microscope.
  • FIG. 1D the horizontal position of the container 10 (microplate 9) into which the sample solution is injected is moved to detect light in the sample solution.
  • the stage position changing device 17a may be operated to move the relative position of the region (a method of moving the absolute position of the sample solution).
  • the mirror deflector 17 or the stage position changing device 17a cooperates with the light detection by the light detector 16 under the control of the computer 18 in order to achieve a desired movement pattern of the position of the light detection region.
  • the movement trajectory of the position of the light detection region may be arbitrarily selected from a circle, an ellipse, a rectangle, a straight line, a curve, or a combination thereof (so that various movement patterns can be selected by the program in the computer 18). It may be.)
  • the position of the light detection region may be moved in the vertical direction by moving the objective lens 8 or the stage up and down.
  • the above optical system is used as a multiphoton microscope. In that case, since only the focal region of the excitation light (light detection region) is the emission of light, a pinhole 13 may be removed. Further, when the luminescent particles or the substance that gives the background light emits light regardless of the excitation light due to chemiluminescence or bioluminescence phenomenon, the optical systems 2 to 5 for generating the excitation light may be omitted. When the luminescent particles or the substance providing the background light emits light by phosphorescence or scattering, the optical system of the confocal microscope is used as it is.
  • a plurality of excitation light sources 2 may be provided so that the wavelength of the excitation light can be appropriately selected according to the excitation wavelength of the substance that gives the luminescent particles or the background light. It may be.
  • each photodetector 16 may be configured to separately detect different components of light from the light detection region. As will be described in more detail later, it is possible to selectively detect single particles having specific emission characteristics by appropriately selecting the components to be detected.
  • a mechanism for dividing the optical path in an arbitrary manner is provided in the detection light optical path after passing through the pinhole 13. Provided. For example, when the light from the photodetection region is detected by being divided into different polarization components, the polarization beam splitter 14a is inserted into a portion to which the reference light 14a is attached.
  • a polarizer (not shown) is inserted in the excitation light optical path. Further, by inserting a dichroic mirror 14a that reflects light in a specific wavelength band in the detection light path and transmits another wavelength band into the detection light path 14a, light components in different wavelength bands are separated. Can be detected.
  • the computer 18 includes a CPU and a memory, and the CPU executes various arithmetic processes to execute the procedure of the present invention. Each procedure may be configured by hardware. All or a part of the processing described in the present embodiment may be executed by the computer 18 using a computer-readable storage medium storing a program for realizing the processing. That is, the computer 18 may realize the processing procedure of the present invention by reading a program stored in a storage medium and executing information processing / calculation processing.
  • the computer-readable recording medium may be a magnetic disk, a magneto-optical disk, a CD-ROM, a DVD-ROM, a semiconductor memory, or the like.
  • the above-mentioned program is distributed to a computer via a communication line. The computer that has received the distribution may execute the program.
  • the optical analysis technology of the present invention can be briefly described as a scanning molecule counting method or an inverted scanning molecule counting method.
  • the light intensity on the time-series light intensity data is simply Rather than refer to the increase or decrease of the value, the pattern of the temporal change in the measured light intensity value occurs with a higher probability whether the particle is present in the light detection region or not. It is estimated whether it is a pattern.
  • the presence / absence and / or number of signals of a single particle to be observed is detected on the time-series light intensity data.
  • the scanning molecule counting method of the present invention and the principle of detection of a single particle signal to be observed will be described.
  • the observation target particle is a luminescent particle
  • the observation target particle when passing through a region where one luminescent particle exists during the movement of the light detection region CV (time to to t2 in the figure) (t1), Light is emitted from the luminescent particles, and a pulse-like signal having a significant light intensity (Em) appears on the time-series light intensity data as illustrated in FIG. Therefore, the movement of the position of the light detection region CV and the light detection are performed, and pulse-like signals (significant light intensity) appearing in the meantime illustrated in FIG. 2B are detected one by one.
  • the luminescent particles are individually detected, and the number of the luminescent particles is counted to obtain information on the number, concentration, or number density of the luminescent particles present in the measured region.
  • the particle to be observed is a non-luminescent particle
  • background light is emitted from the light detection region (or the light detection region is illuminated with illumination light) in the optical measurement by the scanning molecule counting method.
  • It is possible to detect the presence of non-luminous particles by capturing the decrease in the background light detected when the observation target particle enters the light detection region (inverted scanning molecule counting method).
  • light detection is performed while moving the position of the light detection region CV in the sample solution.
  • the luminescent substance is dispersed in the sample solution, a large number of luminescent substances exist in the light detection region CV, and the light detection region CV is moving (time to to t2 in the figure).
  • the scanning molecule counting method As described above, on the actually measured time-series light intensity data, in addition to the light emitting particle signal, the light intensity value due to thermal noise, stray light, Raman scattering of water, etc. There is an increase (noise signal). In the case of the inverted scanning molecule counting method, there is a reduction in the light intensity value due to the fluctuation of the intensity of the background light (also referred to as a noise signal hereinafter).
  • the intensity value of a single particle signal becomes weak, it becomes difficult to distinguish a single particle signal from a noise signal based on the intensity, time width, and shape of the signal.
  • the observation target particle is a luminescent particle
  • the profile of the light intensity value is actually shown in FIG. It does not become a smooth bell shape as depicted in (B), but becomes discrete as illustrated on the left in FIG. Therefore, when the number of photons emitted from the luminescent particles decreases, the intensity, time width, and shape (in FIG. 3A) of the signal on the light intensity data in that case do not include particles but only noise signals.
  • the intensity value of the light emitted and detected from the luminescent particles in the light detection region decreases as the position of the luminescent particles moves away from the approximate center of the light detection region (see FIG. 3D), and the light since the number of the more luminescent particles away from the approximate center of the detection area is increased, the number of hard particles of signal identification of the noise signal increases.
  • a similar situation occurs in the case of a decrease in the number of photons of background light in the inverted scanning molecule counting method in which the observation target particles are non-light emitting particles.
  • the accuracy and / or sensitivity thereof is improved as the detection accuracy and the detection number of the particle signal per certain moving distance of the light detection region are higher.
  • reducing the signal strength criterion to increase the number of detected particle signals increases the possibility of erroneously determining a noise signal as a particle signal and improves the detection accuracy of the particle signal. Therefore, if the signal strength determination reference value is increased, many low-brightness particle signals will not be detected. Therefore, in the present invention, a novel algorithm is proposed that can more accurately achieve the detection of a single particle signal that is difficult to distinguish from a noise signal depending on the intensity, time width, and shape of the signal.
  • the light intensity value is the number of detected photons per bin time. Therefore, as schematically illustrated in FIG. 3A, the light intensity values are discretely distributed in the time axis direction.
  • the particle signal with high luminance has little influence of the noise signal and has a substantially bell-shaped profile (left figure), but the particle signal with low luminance has the intensity value almost the same as that of the noise signal, and further noise Since signals are superimposed (middle diagram), it is difficult to extract a substantially bell-shaped profile, and it is difficult to distinguish from a time domain (right diagram) where only noise signals without particles are present.
  • the algorithm for detecting a single particle signal light detection is performed for a time change (photon number sequence) of a light intensity value for each predetermined time width on time-series light intensity data. If it is assumed that particles exist in the region, the probability that the light intensity value will change with time (occurrence probability (when particles are present), second occurrence probability) and particles exist in the light detection region When it is assumed that the light intensity value has not been generated, the probability of occurrence of a temporal change in the light intensity value (occurrence probability (when no particle is present), first occurrence probability) is calculated. And it is estimated that the state which gives a high occurrence probability is an actual state.
  • an interval of an arbitrary time width (hereinafter referred to as an analysis window) is set in the time-series light intensity data.
  • the number of photon detection events occurring in each unit time is considered to follow a Poisson distribution having an expected value in the unit time. Therefore, in any unit time ti, the number of detected photons Ci
  • the probability (unit time occurrence probability) becomes Given by.
  • Ei is an expected value of the number of photons per unit time ti.
  • the analysis window includes n + 1 unit times, the probability P (occurrence probability) of occurrence of the detected photon number sequence Ci in the analysis window is Given by.
  • the expected value Ei of the number of occurrences of photon detection events in each unit time ti is determined by the presence or absence of particles in the light detection region in the time region corresponding to the analysis window during optical measurement.
  • Bg is the time average value of the noise signal when the observation target particle is a luminescent particle
  • Bg is the time average value of the background light when the observation target particle is a non-luminescent particle.
  • the unit time occurrence probability Pni in each measurement unit time ti is calculated in time series, and further, in equation (2)
  • the probability Pn (first occurrence probability) of occurrence of the actual detected photon number sequence is calculated when it is assumed that no particle is present.
  • the position of the light detection region CV is moved, so that the particle is in the light detection region CV as schematically illustrated in FIG. It will pass through.
  • the intensity value of the light emitted and detected from the particles in the light detection region or the background light reduction amount due to the presence of the particles in the light detection region is shown in FIG.
  • the position of the particle decreases with increasing distance from the approximate center of the light detection region. Therefore, the expected value Epi for each unit time ti when particles are present in the light detection region is also expressed as a bell-shaped function with time as a variable.
  • the bell-shaped function is approximated by a Gaussian function, the expected value Epi is Given by.
  • the Gaussian function is assumed to have a peak intensity Q at an arbitrary time tc (for example, the center of the analysis window) within the analysis window.
  • the full width at half maximum of the Gaussian function of the equation (4) indicates that the light detection region at the moving speed v is emitted from the particles in the light detection region as illustrated in FIG. Equal to the time d / v passing through the full width at half maximum d of the distribution in the radius r direction of the reduced amount of background light that is reduced by the presence of particles in the intensity value or light detection region, from which w is Given by.
  • the full width at half maximum d in FIG. 3D can be determined from the optical system.
  • the peak intensity Q in equation (4) is assumed that the total number of photons in the analysis window matches the total expected value of equation (4). Given by. Further, when the observation target particle is a non-light emitting particle, the expected value of the absolute value of the reduction amount of the photon number is set according to the equation (4), and the peak intensity Q is Given by.
  • the unit time occurrence probability Ppi at each measurement unit time ti is calculated in time series, and further, in equation (2)
  • the probability Pp (second occurrence probability) that the actual sequence of detected photons is generated when the state in which particles are present is assumed is calculated.
  • FIG. 4A shows an example of the series of processes described above.
  • a photon number sequence having a photon number of [00011100] is detected in each unit time (rectangle) on the time axis (arrow) as shown in the middle of the figure.
  • the expected value of each unit time is 0.1 over the entire area when no particles are present, as shown in the upper right diagram.
  • the expected value of each unit time when it is assumed that particles are present is calculated from the equation (6) using the number of detected photons and the equation (4) using the result, as shown in the upper left figure.
  • the value has a bell-shaped profile.
  • the probability that the number of detected photons is generated per unit time is assumed to be absent. It is calculated for each of the case of assuming that (lower right) and particles were present (lower left). Then, the probability of occurrence of the detected photon number is summed up according to the equation (2), and the occurrence probability Pn of the photon number sequence and the particle are present when it is assumed that no particle exists for the detected photon number sequence. If it is assumed, the occurrence probability Pp of the photon number sequence is calculated. In the illustrated example, since the probability Pp> probability Pn, detected photons sequence is determined as a signal of the particle.
  • whether or not a particle signal exists in the analysis window is determined by the odds ratio OR between the occurrence probability Pp and the occurrence probability Pn.
  • OR Pp (1-Pn) / (1-Pp) Pn (8)
  • the presence of the particle signal in the analysis window may be determined when the magnitude exceeds a predetermined value.
  • the analysis window is preferably set to be longer than the time required for a single particle to pass through the light detection region. If the light detection area of radius r is moving at a speed v, the time width of the analysis window is 2r / v (9) Is set longer.
  • the analysis window may be preferably set in order for each unit time on the time-series light intensity data. According to this setting, the occurrence probability Pp, the occurrence probability Pn, and / or the odds ratio OR is calculated along the time-series light intensity data. However, in this case, since the amount of calculation increases, the analysis window may be set every several unit times. Further, the analysis window may be set by dividing the time series light intensity data by the time width of the analysis window. In this case, the analysis windows are set without overlapping each other.
  • the bin time is set to be equal to or less than the time (Equation 9) required for a single particle to pass through the light detection region. This is because a signal when a single particle passes is captured over a plurality of bin times so that a temporal change pattern of the light intensity value while the single particle passes can be detected. (If the bin time is longer than the time of Equation 9, the temporal change pattern of the light intensity value during the passage of a single particle cannot be captured.)
  • the single particle signal detection method using the generation probabilities Pp and Pn described above separately measures multiple different light components and measures each component separately. It can be extended when generating the sequence light intensity data. In that case, by selecting the detection component so that the emission characteristic of the single particle to be observed is reflected in the time-series light intensity data for each component, the emission characteristic of the single particle to be observed is selected. It is possible to selectively detect the reflected signal.
  • selectively detecting a signal reflecting the emission characteristics of a single particle will be described.
  • Cpi, Csi, Bgp, and Bgs are the numbers of p-wave and s-wave photons per unit time and the intensities of s-wave and p-wave background light, respectively.
  • Epi is a value obtained in Equation (4) by giving Q by Equation (12) and Bg by Bgp + Bgs.
  • the occurrence probabilities Pppi and Ppsi of the number of detected photons per unit time in the fluorescence p-wave and s-wave when it is assumed that particles are present are the same as in equation (1). Given by. Therefore, the occurrence probability Pp of the detected photon number sequence in the analysis window when it is assumed that the luminescent particles having the fluorescence anisotropy R exist is the sum of the occurrence probabilities Pppi and Ppsi. Given by. On the other hand, assuming that no particles are present, the occurrence probability Pn of the detected photon sequence in the analysis window is It becomes.
  • Pnpi and Pnsi are the occurrence probabilities of each component per unit time, respectively, using the background light Bgp and Bgs as expected values, respectively, as in the equation (1) It becomes.
  • each time-series light intensity data of the p-wave and s-wave is generated per unit time according to equations (14) and (17).
  • the occurrence probabilities Pp and Pn in the analysis window according to the equations (15) and (16) are calculated, and these occurrence probabilities are compared, whereby the signal of the particle having the fluorescence anisotropy R is obtained.
  • the fluorescence anisotropy R is a known number. Therefore, as the value of the fluorescence anisotropy R, a value experimentally or theoretically determined by an arbitrary method may be used.
  • the particle signal can be detected while identifying the type. That is, when identifying the type of particle, the occurrence probability when the particle as described above is assumed to be present is calculated for each value of fluorescence anisotropy for the photon sequence of each analysis window. . Then, it is estimated that particles having fluorescence anisotropy that gives the highest occurrence probability among the occurrence probabilities obtained using different values of fluorescence anisotropy exist on the photon number sequence being analyzed. It will be possible.
  • the ratio of the light intensity (number of photons) detected in different detection wavelength bands depending on the emission wavelength spectrum of the luminescent particles as described above is reflected in the number of photons on the time-series light intensity data by the scanning molecule counting method. . Therefore, in calculating the probability of occurrence of the photon number sequence, it is possible to selectively detect signals of luminescent particles having a specific emission wavelength characteristic by considering the ratio of the number of photons as a predetermined characteristic value. It becomes possible.
  • the expected value Epi of each unit time of the light intensity in a certain analysis window when particles are present is the same as that in the equation (4), and the peak intensity Q is Given by.
  • C ⁇ i, C ⁇ i, Bg ⁇ p, and Bg ⁇ are the photon numbers of Ch1 and Ch2 per unit time and the intensity of background light of Ch1 and Ch2, respectively. Therefore, the expected values Ep ⁇ and Ep ⁇ of the photon numbers of Ch1 and Ch2 are The occurrence probabilities Pp ⁇ i and Pp ⁇ i of the number of detected photons per unit time in Ch1 and Ch2 when it is assumed that particles are present are the same as in equation (1) Given by.
  • Pn ⁇ i and Pn ⁇ i are occurrence probabilities of each component per unit time, respectively, using the background light Bg ⁇ and Bg ⁇ as expected values, respectively, as in the equation (1). It becomes.
  • the occurrence probability per unit time according to the equations (22) and (25) in each of the time-series light intensity data of each detection wavelength band After calculating the above, by calculating and comparing the occurrence probabilities Pp and Pn in the analysis window according to the equations (23) and (24), it is possible to detect the presence or absence of a signal of particles having specific emission wavelength characteristics.
  • the light intensity ratio ⁇ : ⁇ between Ch1 and Ch2 is a known number.
  • the intensity ratio a value experimentally or theoretically determined by any method may be used.
  • a signal of particles having a specific emission wavelength characteristic can be selectively detected, so that a plurality of types of emission particles having different emission wavelength characteristics are mixed in the sample solution.
  • the particle signal can be detected while identifying the type. That is, when identifying the type of particle, the probability of occurrence of the photon sequence of each analysis window assuming that such a particle is present is the light intensity ratio ( ⁇ : ⁇ ) of Ch1 and Ch2. ) Every time. Then, it can be estimated that particles having a light intensity ratio that gives the highest generation probability among the generation probabilities obtained by using different light intensity ratios exist on the photon number sequence being analyzed.
  • FIG. 5 shows processing in the present embodiment expressed in the form of a flowchart.
  • the particles to be observed in the optical analysis technique of the present invention are arbitrary as long as they are dispersed particles in the sample solution and move randomly in the solution, such as dissolved molecules.
  • biomolecules such as proteins, peptides, nucleic acids, lipids, sugar chains, amino acids or aggregates thereof, viruses, cells, metal colloids, other non-biological molecules, etc. It's okay.
  • a luminescent label fluorescent molecule, phosphorescent molecule, chemical / bioluminescent molecule
  • the sample solution is typically an aqueous solution, but is not limited thereto, and may be an organic solvent or any other liquid.
  • the inverted scanning molecule counting method it is a non-light emitting particle to be observed, but it may be arbitrary as in the case of the light emitting particle.
  • the luminescent substance that gives background light may be any luminescent molecule, for example, a fluorescent molecule, a phosphorescent molecule, or a chemical / bioluminescent molecule, and the luminescent substance is always present in several or more molecules in the photodetection region. Dissolved or dispersed in the sample solution at a concentration.
  • the sample solution is typically an aqueous solution, but is not limited thereto, and may be an organic solvent or any other liquid.
  • the computer 18 stores a program stored in a storage device (not shown) (procedure for moving the position of the light detection region in the sample solution, and from the light detection region during movement of the position of the light detection region).
  • a program stored in a storage device (not shown) (procedure for moving the position of the light detection region in the sample solution, and from the light detection region during movement of the position of the light detection region).
  • irradiation of excitation light and measurement of light intensity in the light detection region in the sample solution are started.
  • the mirror deflector 17 or the stage position changing device 17a drives the mirror 7 (galvanomirror) or the microplate 9 on the stage of the microscope, and the wells.
  • the position of the light detection region is moved, and at the same time, the light detector 16 sequentially converts the detected light into an electric signal and transmits it to the computer 18.
  • time-series light intensity data is generated from the transmitted signal and stored.
  • the photodetector 16 is an ultra-sensitive photodetector that can detect the arrival of one photon, and therefore when the light detection is based on photon counting, the time-series light intensity data is a time-series data. It may be photon count data.
  • each of the plurality of photodetectors 16 simultaneously detects the light intensity values (number of photons) of those light components, and detects each different component detected. Time series light intensity data is generated.
  • the moving speed of the position of the photodetection region is preferably set to a value faster than the random movement of the single particle, that is, the moving speed due to the Brownian movement.
  • the moving speed of the position of the light detection region is slower than the movement of the particle due to Brownian motion, the particles move randomly in the light detection region. Then, the calculation of the expected value for the occurrence probability calculated assuming the state where particles are present becomes complicated, and the accuracy may be lowered. Therefore, preferably, as shown in FIG. 3C, the particles cross the light detection region in a substantially straight line, and the expected value of the profile of the change in light intensity is drawn in FIG. 3D.
  • the movement speed of the position of the light detection region is the average movement speed due to the Brownian motion of the particle (diffusion movement speed) so that the position of the particle becomes a bell-like shape that decreases as it moves away from the approximate center of the light detection region. ) Is set faster.
  • the moving speed of the position of the photodetection region may be set to 15 mm / s, which is 10 times or more.
  • the profile of the change in light intensity by setting various moving speeds of the position of the light detection region can be expected (typically, the excitation light intensity distribution).
  • a preliminary experiment for finding a condition that is substantially the same as that described above may be repeatedly performed to determine a suitable moving speed of the position of the light detection region.
  • the intensity value of the background light is calculated on the time-series light intensity data (step 110).
  • the intensity value of the background light may be an average value of intensity values (number of photons) in a region where no particle signal exists on the time-series light intensity data.
  • the obtained time-series light intensity data An average value of all intensity values excluding data of a predetermined ratio (for example, 20%) from the higher one among the light intensity values may be adopted as the intensity value of the background light.
  • a predetermined ratio of data is considered to be a particle signal from the higher light intensity value.
  • a predetermined ratio for example, from the lower one of the light intensity values on the obtained time-series light intensity data (for example, An average value of all intensity values excluding 20%) data may be adopted as the intensity value of the background light.
  • a predetermined ratio of data from the lower light intensity value is considered to be a particle signal.
  • analysis window setting is executed on the time-series light intensity data (step 120).
  • the length of one analysis window may be determined from the size of the light detection region and the moving speed (see equation (9)).
  • the analysis window may be preferably set in time series for each bin time. However, for the purpose of reducing the amount of calculation, it may be set every several bin times, or the analysis windows may be set so as not to overlap each other.
  • the occurrence probability Pn of the light intensity value sequence or the photon number sequence in the analysis window when it is assumed that there is no particle in the light detection region is calculated, and the occurrence probability Pn and the occurrence probability Pp are compared. It is then determined whether or not particles are present in the analysis window (step 150). In this determination, the odds ratio between the occurrence probability Pn and the occurrence probability Pp is calculated (see equation (8)), and when the odds ratio is equal to or greater than a predetermined value, it may be determined that particles are present.
  • time-series light intensity data is obtained for a plurality of different components and the occurrence probability Pp is obtained in consideration of a predetermined characteristic value of the particle
  • the presence / absence of a signal of a particle having the predetermined characteristic value is selectively It should be understood that it is detected.
  • the expected intensity value may be assumed to have an intensity peak at the center of the analysis window.
  • the peak of the particle signal does not exist at the center of the analysis window in most cases, but the value of the occurrence probability Pp decreases as the position of the peak of the actual particle signal moves away from the center of the analysis window. However, the value is higher than the value of the occurrence probability Pp when there is no actual particle signal.
  • the calculation of the occurrence probabilities Pp and Pn in the analysis window and the detection of the particle signal in the processing of the above steps 130 to 150 may be executed in all analysis windows set in the light intensity data (step 160). ).
  • the volume value and the number of single particles can be used to calculate the volume in the sample solution.
  • the number density or concentration of a single particle is determined.
  • the total volume of the region through which the light detection region has passed may be theoretically calculated based on the wavelength of the excitation light or detection light, the numerical aperture of the lens, and the adjustment state of the optical system.
  • experimentally for example, with respect to a solution having a known particle concentration (control solution), measurement of light intensity, detection of a single particle, and It may be determined from the number of particles detected by performing the counting and the concentration of particles in the control solution.
  • Vt N / C (29) Given by.
  • a control solution a plurality of solutions having different concentrations of particles are prepared, and measurement is performed on each of them, and the calculated average value of Vt is adopted as the total volume Vt of the region through which the light detection region passes. It may be.
  • the volume of the light detection region and the total volume of the region through which the light detection region has passed are given by any method, for example, using FCS or FIDA, without depending on the above method. It's okay. Further, in the optical analyzer of the present embodiment, the relationship between the concentration C and the number N of various standard particles (formula (29)) for the assumed movement pattern of the light detection region.
  • the information may be stored in advance in a storage device of the computer 18 so that the information on the relationship stored as appropriate when the user of the device performs the optical analysis can be used.
  • the count and concentration of the signal of the particle having the predetermined characteristic value are obtained. It should be understood that it is possible to calculate Therefore, when a plurality of types of particles having different characteristic values are mixed in the sample solution, it is possible to calculate the particle count and concentration for each type while identifying the type of particle.
  • the particle signal was detected based on the generation probabilities Pp and Pn of the photon number sequence on the time-series light intensity data according to the present invention.
  • sample solution solutions containing various concentrations of ATTO647N (Sigma) as luminescent particles in Tris buffer (containing 10 mM Tris-HCl (pH 8.0), 0.05% Tween20) were prepared.
  • Tris buffer containing 10 mM Tris-HCl (pH 8.0), 0.05% Tween20
  • a single molecule fluorescence measurement device MF20 (Olympus Corporation) equipped with an optical system of a confocal fluorescence microscope and a photon counting system is used as an optical analysis device, and the above-mentioned “(2) Sample solution according the manner detailed in the measurement of the light intensity "was obtained for the above sample solution, the time-series light intensity data (photon count data).
  • the acquired time-series photon count data is described in “(3) Individual detection of single particle signal” and steps 110 to 160 in FIG.
  • calculation of background light intensity, setting of analysis window, occurrence probability Pn of photon sequence in analysis window when it is assumed that there is no particle in the light detection area, and case that particle is in the light detection area The calculation of the generation probability Pp of the photon sequence in the analysis window and the calculation of the odds ratio of the generation probability were sequentially executed, and a single particle signal was detected and counted.
  • the time-series light intensity data is only a component in one detection wavelength band
  • the occurrence probability Pn and the occurrence probability Pp are as described in “2. Detection of a signal of a single particle”. It calculated using Formula (1), (2), (6) etc. which were described in the column.
  • FIG. 6 shows a part (A) of time-series light intensity data, data (B) obtained by smoothing (A), and occurrence probabilities Pn and Pp calculated from (A) according to the teaching of the present invention.
  • One example of the odds ratio (C) is shown.
  • 633 nm laser light was used as the excitation light, and light in the wavelength band of 660 to 710 nm was measured using a bandpass filter to generate time-series photon count data.
  • the moving speed of the position of the light detection region in the sample solution was 3000 rpm (7.5 mm / sec), and BIN TIME was 10 ⁇ sec.
  • the intensity of the background light is the total of all data except the higher and lower 20% of the total number of photons of 1 ms around the illustrated time domain.
  • the average number of photons in the data was used.
  • the full width at half maximum (d / v) of the Gaussian function (Expression (4)) representing the distribution of the expected value is 80 ⁇ sec. .
  • the time width of the analysis window was 300 ⁇ sec.
  • the concentration of ATTO647N in the sample solution was 1 pM.
  • the magnitude of the signal intensity in the smoothed data is a particle having a predetermined threshold value.
  • the intensity of the pulsed signal depending on whether or not exceeds a predetermined threshold value, whether the particles of the signal is determined. Since the predetermined threshold is typically set to 1 photon, all of the predetermined threshold values are determined as noise except for the increase in the number of photons shown at (a) in the approximate center of the figure. Since the peak value of the pulse signal shown in (a) is the number of photons ⁇ 1, it is assumed that it is a particle signal, but the determination result changes due to the influence of an error and is determined as a noise signal. In some cases (ie, the result is unstable).
  • FIG. 7 is a diagram showing another part of the time-series light intensity data similar to FIG.
  • the moving speed of the position of the light detection region in the sample solution was set to 6000 rpm (15 mm / second).
  • a significant pulse-like signal is observed in addition to the pulse-like signal shown in (a).
  • the occurrence probabilities Pn, Pp In the odds ratio (C) no significant increase was observed except in the region corresponding to the pulse signal shown in (a). This suggests that the determination based on the occurrence probabilities Pn and Pp calculated according to the teaching of the present invention makes the distinction between the particle signal and the noise signal clear, and improves the stability of particle determination. Yes.
  • FIG. 8 shows the relationship between the luminescent particle concentration and the number of detected particles when the particle signal is determined using the odds ratio of the generation probabilities Pn and Pp according to the present invention for sample solutions containing luminescent particles of various concentrations. ing.
  • 642 nm laser light (3 mW) was used as excitation light, and light in a wavelength band of 560 to 620 nm was measured using a bandpass filter to generate time-series photon count data.
  • the moving speed of the position of the light detection region in the sample solution was 12000 rpm (90 mm / second), the BIN TIME was 10 ⁇ sec, and the light measurement time was 600 seconds.
  • the intensity of the background light is the number of photons of all data except the higher and lower 20% of the 1 msec data around each analysis window. The average value was used.
  • the full width at half maximum (d / v) of the Gaussian function (Expression (4)) representing the distribution of the expected value is 50 ⁇ sec. .
  • the time width of the analysis window was 350 ⁇ sec.
  • the number of detected particles increased with the concentration of the luminescent particles in the concentration range of 1 fM or more. This indicates that the particle signal detection method of the present invention can detect the particle signal on the time-series photon count data, and that the particle concentration in the sample solution can be determined in a concentration region of 1 fM or more. Yes.
  • the particle signal was detected based on the occurrence probabilities Pp and Pn of the photon sequence on the time-series light intensity data according to the present invention.
  • a tris buffer solution containing 10 mM Tris-HCl (pH 8.0), 0.05% Tween 20
  • 1 mM ATTO647N a luminescent substance that generates background light.
  • a single molecule fluorescence measurement device MF20 (Olympus Corporation) equipped with an optical system of a confocal fluorescence microscope and a photon counting system is used as an optical analysis device, and the above-mentioned “(2) Sample solution according the manner detailed in the measurement of the light intensity "was obtained for the above sample solution, the time-series light intensity data (photon count data).
  • the acquired time-series photon count data is described in “(3) Individual detection of single particle signal” and steps 110 to 160 in FIG.
  • calculation of background light intensity, setting of analysis window, occurrence probability Pn of photon sequence in analysis window when it is assumed that there is no particle in the light detection area, and case that particle is in the light detection area The calculation of the generation probability Pp of the photon number sequence in the analysis window and the calculation of the odds ratio of the generation probability were sequentially performed, and the signal of a single particle was counted.
  • the time-series light intensity data is only a component in one detection wavelength band
  • the occurrence probability Pn and the occurrence probability Pp are as described in “2. Detection of a signal of a single particle”. It calculated using the formula (1), (2), (7) etc. which were described in the column.
  • FIGS. 9A, 9B, and 10 show enlarged diagrams (upper stage) of the obtained time-series light intensity data in order and their odds ratios.
  • the excitation light is 633 nm laser light (50 ⁇ W)
  • the bandpass filter is used to measure the light in the wavelength band of 660-710 nm
  • the time series photon count data is obtained. Generated.
  • the moving speed of the position of the light detection region in the sample solution was 6000 rpm (15 mm / sec), and BIN TIME was 50 ⁇ sec.
  • the intensity of the background light is the number of photons of all data except the higher and lower 20% of the 1 msec data around each analysis window.
  • the average value was used.
  • the full width at half maximum (d / v) of the Gaussian function (Expression (4)) representing the distribution of the expected value is 100 ⁇ sec.
  • the time width of the analysis window was 300 ⁇ sec.
  • the concentration of beads was 100 fM.
  • FIG. 11A shows the relationship between the particle concentration and the number of detected particles when the particle signal is determined using the odds ratio of the generation probabilities Pn and Pp according to the present invention for sample solutions containing beads of various concentrations. Show. In this experiment, 633 nm laser light (50 ⁇ W) was used as excitation light, and light in a wavelength band of 660 to 710 nm was measured using a bandpass filter to generate time-series photon count data. The moving speed of the position of the light detection region in the sample solution was 9000 rpm (90 mm / second), the BIN TIME was 10 ⁇ sec, and the light measurement time was 200 seconds.
  • the intensity of the background light is the number of photons of all data except the higher and lower 20% of the 1 msec data around each analysis window. The average value was used.
  • the full width at half maximum (d / v) of the Gaussian function (Expression (4)) representing the distribution of the expected value is 100 ⁇ sec. .
  • the time width of the analysis window was 300 ⁇ sec. Then, it was determined that the particle signal was present in a region where the odds ratio> 10 10 . For comparison, FIG.
  • 11B shows that the same data is smoothed with time-series photon count data as described in Patent Documents 9 to 11, and then fitted with a Gaussian function so that the number of peak photons ⁇ 5. The relationship between the number of detected particles obtained by detecting a certain pulse signal and the particle concentration is shown.
  • the figure shows that the particle concentration in the sample solution can be determined in a concentration range of 50 aM or more.
  • the determinable concentration range is 500 aM or more.
  • a solution containing a fluorescent dye TAMRA (MW 430.45, Sigma-Aldrich Cat. No. C2734) at 100 fM in a phosphate buffer (containing 0.05% Tween 20) (fluorescent dye solution)
  • SYTOX Orange binds to a single plasmid to form a single luminescent particle).
  • a single molecule fluorescence measurement device MF20 (Olympus Corporation) equipped with an optical system of a confocal fluorescence microscope and a photon counting system is used as an optical analysis device, and the above-mentioned “(2) Sample solution According to the procedure described in “Measurement of light intensity”, time-series light intensity data (photon count data) of the s-polarized component and the p-polarized component were obtained simultaneously and separately for the sample solution. At that time, for both the TAMRA solution and the plasmid solution, 543 nm laser light was used as the excitation light, and light in the wavelength band of 560 to 620 nm was detected using a bandpass filter.
  • the polarization direction of the excitation light was set to be the same direction as the p-polarized component of the detection light.
  • the moving speed of the light detection region in the sample solution was set to 6000 rpm (15 mm / second), the BIN TIME was set to 10 ⁇ sec, and the measurement time was set to 2 seconds. Time series light intensity data (photon count data) was obtained for the sample solution.
  • the acquired time-series photon count data is described in “(3) Individual detection of single particle signal” and steps 110 to 160 in FIG.
  • calculation of background light intensity, setting of analysis window, occurrence probability Pn of photon sequence in analysis window when it is assumed that there is no particle in the light detection area, and case that particle is in the light detection area The calculation of the generation probability Pp of the photon number sequence in the analysis window and the calculation of the odds ratio of the generation probability were sequentially performed, and the signal of a single particle was counted.
  • the occurrence probability Pn and the occurrence probability Pp including the fluorescence anisotropy as a known parameter are determined from “3. Calculated using equations (15), (16), etc.
  • the intensity of the background light was defined as an average value of the number of photons of all data except 20% of the higher and lower photons in the data of 1 ms around each analysis window.
  • the full width at half maximum (d / v) of the Gaussian function (Expression (4)) representing the distribution of the expected value is 120 ⁇ sec.
  • the time width of the analysis window was 400 ⁇ sec.
  • occurrence probability Pp a value using 0.4 of the fluorescence anisotropy of the plasmid and a value using 0.32 of the fluorescence anisotropy of TAMRA are calculated (that is, for one analysis window).
  • Two occurrence probabilities Pp are calculated.) For each of them, when the odds ratio with the occurrence probability Pn exceeds 20, it is determined that there is a signal of the type of particle corresponding to the analysis window. (When the odds ratio of the occurrence probability Pp of the plasmid and the odds ratio of the probability Pp of TAMRA exceeds 20, it is assumed that there is a particle signal with a high odds ratio.)
  • FIG. 12 shows the occurrence probability Pp obtained using the fluorescence anisotropy different from the occurrence probability Pn in the time-series light intensity data obtained by the scanning molecule counting method for the plasmid solution and the fluorescent dye solution. using odds ratios, to determine the presence of particles signals shows the results of counting the number of particles for each type of particle.
  • most of the detected particle signals were detected corresponding to the respective types in each solution.
  • 421 particle signals were detected, of which 82 were erroneously detected as TAMRA.
  • the fluorescent dye solution 279 particle signals were detected, of which The number of falsely detected plasmids was 58.
  • the occurrence probability of the light intensity value sequence in the time-series light intensity data is assumed to be present.
  • the detection accuracy of the signal of the particle or the S / S is calculated by determining the presence or absence of the signal of the particle based on the calculated occurrence probability and the case of assuming that the particle does not exist. It has been shown that the N ratio is improved and the detection sensitivity (the range of the particle concentration of the sample solution in which the particle signal can be detected with an acceptable accuracy) can be improved.

Abstract

共焦点顕微鏡又は多光子顕微鏡による光計測を用いた走査分子計数法に於いて、その感度及び／又はＳ／Ｎ比の向上を可能にする技術が提供される。本発明による試料溶液中の単一粒子の存在を検出する光分析技術では、試料溶液内に於いて光検出領域の位置を移動し、光検出領域からの光の強度を測定して光強度データを生成し、光強度データ上の光強度の時間変化の、光検出領域内に単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、光検出領域内に単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出し、それらの発生確率に基づいて光検出領域内に単一粒子の存在した時間を検出して単一粒子の各々の存在を表す信号を検出する。

Description

単一粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
　本発明は、共焦点顕微鏡又は多光子顕微鏡の光学系などの溶液中の微小領域からの光が検出可能な光学系を用いて、溶液中に分散又は溶解した原子、分子又はこれらの凝集体（以下、これらを「粒子」と称する。）、例えば、タンパク質、ペプチド、核酸、脂質、糖鎖、アミノ酸若しくはこれらの凝集体などの生体分子、ウイルス、細胞などの粒子状の対象物、或いは、非生物学的な粒子を検出して、それらの状態（相互作用、結合・解離状態など）の分析又は解析に於いて有用な情報を取得することが可能な光分析技術に係り、より詳細には、上記の如き光学系を用いて単一の粒子の存在による光の変化（単一粒子が発する光又は単一粒子が生ずる影）を個別に検出して種々の光分析を可能にする光分析装置、光分析方法並びに光分析用コンピュータプログラムに係る。なお、本明細書に於いて、検出される光は、蛍光、りん光、化学発光、生物発光、散乱光等であってよい。光を発する粒子（以下、「発光粒子」と称する。）は、それ自身が光を発する粒子、又は、任意の発光標識若しくは発光プローブが付加された粒子のいずれであってもよい。
　近年の光計測技術の発展により、共焦点顕微鏡の光学系とフォトンカウンティング（１光子検出）も可能な超高感度の光検出技術とを用いて、一光子又は蛍光一分子レベルの微弱光の検出・測定が可能となっている。そこで、そのような微弱光の計測技術を用いて、生体分子等の特性、分子間相互作用又は結合・解離反応の検出を行う光分析技術が種々提案されている。そのような光分析技術としては、例えば、蛍光相関分光分析（Fluorescence　Correlation　Spectroscopy：ＦＣＳ。例えば、特許文献１－３、非特許文献１－３参照）、蛍光強度分布分析(Fluorescence-Intensity　Distribution　Analysis：ＦＩＤＡ。例えば、特許文献４、非特許文献４)やフォトンカウンティングヒストグラム（Photon　Counting　Histogram：ＰＣＨ。例えば、特許文献５）などが知られている。また、特許文献６～８には、共焦点顕微鏡の光学系を用いて計測される試料溶液の蛍光信号の時間経過に基づいて蛍光性物質を検出する方法が提案されている。
　更に、本願出願人は、特許文献９～１１に於いて、共焦点顕微鏡又は多光子顕微鏡の光学系などの溶液中の微小領域からの光が検出可能な光学系を用いた光分析技術であって、ＦＣＳ、ＦＩＤＡ等の光分析技術とは異なる原理による新規な光分析技術を提案した。上記のＦＣＳ、ＦＩＤＡ等の光分析技術の場合、試料溶液内に於ける光の検出領域である微小領域（以下、「光検出領域」と称する。）に浮遊する蛍光分子からの光を連続的に計測することにより得られた光強度データが統計的なゆらぎを算出する演算処理により解析され、蛍光分子の濃度及び／又はその他の特性が検出される。これに対し、特許文献９～１１に於いて提案された新規な光分析技術では、試料溶液内に於いて光検出領域の位置を移動させながら、即ち、光検出領域により試料溶液内を走査しながら、光検出領域が試料溶液中に分散してランダムに運動する発光粒子を包含したときに、その発光粒子から発せられる光を個別に検出し、これにより、試料溶液中の発光粒子を一つずつ検出して、発光粒子のカウンティングや試料溶液中の発光粒子の濃度又は数密度に関する情報の取得が可能となる。この光分析技術（以下、「走査分子計数法」と称する。）によれば、測定に必要な試料は、ＦＣＳ、ＦＩＤＡ等の光分析技術と同様に微量（例えば、数十μＬ程度）であってもよく、また、測定時間が短く（一回の測定で秒オーダーの時間の計測が数回繰り返される。）、更に、観測対象となる粒子の濃度が、ＦＣＳ、ＦＩＤＡ等の光分析技術の良好に測定可能なレベル（１ｎＭ程度）よりも低い試料溶液に於いて、発光粒子の存在を検出し、その濃度、数密度又はその他の特性を定量的に検出することが可能となる。
　かくして、「走査分子計数法」は、医学・生物学の研究開発の分野でしばしば使用される希少な或いは高価な試料についての分析を行う場合や、病気の臨床診断や生理活性物質のスクリーニングなど、検体数が多い場合に、従前の生化学的方法に比して、低廉に、或いは、迅速に実験又は検査が実行でき、しかも、ＦＣＳ、ＦＩＤＡ等を良好に実行できない程度のより低い濃度の粒子の濃度及び／又は特性の検出が可能な強力なツールとして期待される。
特開２００５－０９８８７６ 特開２００８－２９２３７１ 特開２００９－２８１８３１ 特許第４０２３５２３号 国際公開２００８－０８０４１７ 特開２００７－２０５６５ 特開２００８－１１６４４０ 特開平４－３３７４４６号公報 国際公開第２０１１／１０８３６９ 国際公開第２０１１／１０８３７０ 国際公開第２０１１／１０８３７１
金城政孝、蛋白質　核酸　酵素　Vol.４４、No.９、１４３１－１４３８頁　１９９９年 エフ・ジェイ・メイヤー・アルムス（F.J.Meyer-Alms）、フルオレセンス・コリレーション・スペクトロスコピー（Fluorescence　Correlation　Spectroscopy）、アール・リグラー編（R.Rigler）、スプリンガー（Springer）、ベルリン、２０００年、２０４－２２４頁 加藤則子外４名、遺伝子医学、Vol.６、No.２、２７１－２７７頁 カスク他３名、米国科学アカデミー紀要　１９９９年、９６巻、13756‐13761頁（P.　Kask,　K.　Palo,　D.　Ullmann,　K.　Gall　PNAS　96,　13756-13761　(1999)） Anal.　Chem.　Vol.70　No.3　６３２頁　１９９８年
　上記の特許文献９～１１に記載の走査分子計数法では、典型的には、光検出領域からの光強度の時系列データ（時系列光強度データ）に於いて、釣鐘状の光強度の変化が観測対象となる単一粒子の存在に対応するとの仮定の下、かかる釣鐘状の光強度の変化の検出が行われる。例えば、観測対象粒子が発光粒子である場合には（後述の如く、観測対象粒子が検出波長帯域の光を発しない粒子（非発光粒子）である場合もあり得る。）、時系列光強度データに観測された釣鐘状又はパルス状の光強度の増大のうち、発光粒子が光検出領域内を通過した際に想定される光強度の増大のプロファイル（ピーク強度、半値全幅等）の条件に適合する光強度の変化が粒子の存在を表す信号として検出される。そして、発光粒子の光強度の増大のプロファイルの条件に適合しない光強度の増大は、ノイズとして判定される。
　しかしながら、時系列光強度データに於いて上記の如き釣鐘状の信号を検出する場合、粒子の信号の大きさが弱くなると、粒子の信号とノイズ信号との識別が困難となる。
　後述の実施形態の欄に於いてより詳細に説明される如く、光強度の測定は、観測対象粒子による微弱な光強度又はその変化を捉えるために、典型的には、フォトンカウンティングにより実行され、時系列光強度データは、離散的なフォトンカウンティングデータとなる（図２（Ａ）参照）。その場合、釣鐘状の信号の検出を容易にするために、好適には、時系列光強度データが時間について平滑化され、平滑化された時系列光強度データに於いて粒子の光強度の変化のプロファイル条件に適合する釣鐘状の信号が検出される。しかしながら、かかる態様に於いては、粒子の信号の強度変化が小さいとき、その信号は、フォトンカウンティングデータの段階でノイズに埋もれてしまい、平滑化処理後に於いても、粒子の光強度の変化の釣鐘状のプロファイルの特徴を呈さず、粒子の信号として検出されないこととなる。特に、粒子の存在に対応する信号の大きさは、光検出領域内での粒子の通過経路によって変化し、光検出領域の周辺部を通過する粒子の信号の強度変化は小さい。従って、微弱な多くの光検出領域の周辺部を通過する粒子の信号が見落とされ、或いは、ノイズ信号が微弱な粒子の信号として誤判定されることとなり、このことが、走査分子計数法の感度又は精度の向上の妨げの一つとなっていた。
　かくして、本発明の主な課題は、上記の走査分子計数法に於いて、その感度及び／又は精度の更なる向上を可能にする新規な粒子信号の検出手段又は手法を提供することである。
　また、本発明のもう一つの課題は、上記の走査分子計数法に於いて、微弱な粒子の信号とノイズ信号とのより精度の良い識別を可能にする新規な粒子信号の検出手段又は手法を提供することである。
　上記の課題に関して、本発明の発明者は、走査分子計数法の時系列光強度データに於いて、粒子の信号に対応する部分と、ノイズ信号の部分とでは、光強度の時間変化（光子数列）の発生パターンに違いがあることを見出した。即ち、ノイズ信号の場合には、常にランダムに光子が検出されるのに対して、粒子の信号の場合には、光子が時間的に集中して検出される。本発明に於いては、かかる光強度の時間変化の発生パターンの違いを検出可能な新規なアルゴリズムにて走査分子計数法に於ける粒子の信号の検出精度の向上が図られる。また、光強度の時間変化の発生パターンは、粒子の特性（偏光特性や発光波長特性）によっても変化するので、光強度の時間変化の発生パターンの違いに基づいて、異なる特性の粒子の識別も試みられる。
　かくして、本発明によれば、上記の課題は、共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出する光分析装置であって、試料溶液内に於ける顕微鏡の光学系の光検出領域の位置を移動する光検出領域移動部と、光検出領域からの光を検出する光検出部と、試料溶液内に於いて光検出領域の位置を移動させながら光検出部にて検出された光検出領域からの光の時系列の光強度データを生成し、時系列の光強度データに於いて単一粒子の各々の存在を表す信号を個別に検出する信号処理部とを含み、信号処理部が、時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、光検出領域内に単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、光検出領域内に単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出し、第一及び第二の発生確率に基づいて時系列の光強度データ上にて単一粒子の各々の存在を表す信号を検出することを特徴とする装置によって達成される。
　上記本発明の構成に於いて、「試料溶液中にて分散しランダムに運動する単一粒子」とは、試料溶液中に分散又は溶解した原子、分子又はそれらの凝集体などの粒子であって、基板などに固定されず、溶液中を自由にブラウン運動している粒子であれば任意の粒子であってよい。観測対象となる単一粒子は、光を発する粒子（発光粒子）であっても、（検出波長帯域に於いて）光を発さない粒子（非発光粒子）であってもよい。観測対象粒子が発光粒子である場合には、「単一粒子の各々の存在を表す信号」は、発光粒子がその光検出領域内の通過時に発する光に対応した時系列光強度データ上に於ける「光強度値の一時的な増大」となる。発光粒子は、典型的には、蛍光性粒子であるが、りん光、化学発光、生物発光、光散乱等により光を発する粒子であってもよい。また、観測対象粒子は非発光粒子であってもよく、その場合には、光検出領域からの光が有意な背景光を含み、「単一粒子の各々の存在を表す信号」が、背景光からの光強度の一時的な低減となる（反転型走査分子計数法）。なお、共焦点顕微鏡又は多光子顕微鏡の光学系の「光検出領域」とは、それらの顕微鏡に於いて光が検出される微小領域であり、対物レンズから照明光が与えられる場合には、その照明光が集光された領域に相当する（共焦点顕微鏡に於いては、特に対物レンズとピンホールとの位置関係により確定される。発光粒子が照明光なしで発光する場合、例えば、化学発光又は生物発光により発光する粒子の場合には、顕微鏡に於いて照明光は要しない。）。また、典型的には、光検出部は、所定の計測単位時間（ビンタイム）毎に到来する光子数を計数するフォトンカウンティングにより光検出領域からの光を検出し、その場合、時系列の光強度データが時系列のフォトンカウントデータとなる。
　上記から理解される如く、本発明の装置に於いては、基本的には、特許文献９～１１に記載の「走査分子計数法」と同様に、試料溶液内に於いて光検出領域の位置を移動しながら、即ち、試料溶液内を光検出領域により走査しながら、逐次的に、光の検出及び時系列の光強度値を表す時系列光強度データの生成が為され、時系列光強度データ上に於いて、単一粒子の存在を表す信号が検出される。かかる単一粒子の存在を表す信号の検出に際して、本発明では、単に光強度値の増減量が所定量を越えたか否かに基づいて単一粒子の信号の検出が実行されるのではなく、時系列光強度データ上の光強度値の時間変化のパターンが、光検出領域内に単一粒子がない場合と光検出領域内に単一粒子が有る場合とのいずれに於いて発生し易いパターンであるかが時系列に評価される。即ち、既に触れた如く、時系列光強度データに於いて、粒子の信号に対応する部分とノイズ信号の部分とでは、光強度の時間変化の発生パターンに違いがあるので、時系列光強度データ上の或る部分に於ける光強度値の時間変化のパターンが、光検出領域内に単一粒子がない場合に発生し易いパターンであるときには、その部分は、ノイズ信号のみの部分であると判定でき、光検出領域内に単一粒子がある場合に発生し易いパターンであるときには、その部分は、粒子の信号に対応する部分であると判定できることとなる。そして、或る光強度値の時間変化のパターンが、光検出領域内に単一粒子がある場合又はない場合のいずれに発生し易いパターンであるかは、その光強度値の時間変化のパターンが、光検出領域内に単一粒子がある場合及びない場合のそれぞれに於いて発生する確率によって判定することができる。
　かくして、本発明の装置では、時系列光強度データ上で光検出領域内に粒子が存在した部分、即ち、単一粒子の存在を表す信号を検出するために、上記の如く、「時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、光検出領域内に単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、光検出領域内に単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率」とが算出される。ここで、「解析窓」とは、時系列光強度データ上に於ける任意の時間幅の領域であり、「解析窓」は、時系列光強度データ上にて逐次的に又は時系列に設定される。また、「第一の発生確率」は、光検出領域内に単一粒子が存在していない場合に、「解析窓」に於いて計測された光強度の時間変化が発生する確率であり、「第二の発生確率」は、光検出領域内に単一粒子が存在している場合に、「解析窓」に於いて計測された光強度の時間変化が発生する確率である。そして、光検出領域内に単一粒子が存在していれば、「第二の発生確率」が「第一の発生確率」に対して相対的に増大することとなるので、「第一の発生確率」と「第二の発生確率」とを参照することによって、時系列の光強度データ上にて光検出領域内に単一粒子の存在した時間領域が決定され、これにより、単一粒子の各々の存在を表す信号の検出が可能となる。
　なお、上記の構成に於いて、解析窓の幅は、例えば、単一粒子が光検出領域を通過するのに要する時間以上の任意の幅に設定されてよい。また、解析窓は、時系列光強度データ上に於いて、単位時間毎に或いは所定の時間間隔毎に設定されてもよく（この場合、時系列に設定される解析窓は互いに重複する。）、或いは、時系列光強度データを所定の時間幅にて分割することにより設定されてもよい。「単位時間」とは、光計測に於いて一つの光強度値を与える幅の時間であり、フォトンカウンティングの場合であれば、１つ又はそれ以上の数のビンタイムの幅の時間であってよい。
　上記の構成に於いて、第一及び第二の発生確率は、それぞれ、解析窓内に於ける第一及び第二の状態を仮定した場合の光強度値の時間変化の平均のパターンからの実際に計測された光強度値の時間変化のパターンのずれに基づいて決定することができる。平均のパターンからの実測値のパターンのずれが小さいほど、その実測値のパターンの発生する確率が高くなる。従って、第一及び第二の発生確率は、それぞれ、解析窓内に於ける単位時間毎の検出された光強度値と第一及び第二の状態を仮定した場合の単位時間毎の期待値とに基づいて算出することが可能である。また、時系列光強度データ上に出現する光強度値は、光検出領域内から放出される光子数又はこれに比例する量であるので、単位時間毎の光強度値は、ポアソン分布に従うと考えられる。そこで、上記の構成に於いて、好適には、単位時間毎の光強度値が単位時間毎の期待値を有するポアソン分布に従うものとして単位時間毎の光強度値の単位時間発生確率が算出され、第一及び第二の発生確率がそれぞれ対応する単位時間発生確率を用いて算出されてよい。
　かくして、上記の装置では、簡単には、第二の発生確率が第一の発生確率より大きい解析窓の時間に光検出領域内に単一粒子が存在したと判定されてよい。また、第一及び第二の発生確率の比又はオッズ比を時系列に算出し、かかる確率の比又はオッズ比の値に基づいて光検出領域内に単一粒子が存在した時間が決定されてもよい。
　ところで、上記の走査分子計数法に於いて、光検出領域からの光の少なくとも二つの互いに異なる成分を別々に検出して成分の各々の時系列の光強度データを生成することが可能である。その場合、検出される成分のデータに於いて、単一粒子の任意の特性が反映されるように、検出される成分を選択することが可能である。例えば、検出される複数の成分として、偏光方向が互いに異なる成分を選択すれば、複数の成分の時系列光強度データに於いて、粒子の偏光特性が反映される。また、検出される複数の成分として、互いに異なる波長帯域の光成分を選択すれば、複数の成分の時系列光強度データに於いて、粒子の発光波長特性が反映される。そして、複数の成分の時系列光強度データに於ける光強度値の時間変化のパターンに於いても、観測対象粒子の前記の如き特性が反映されることとなる。
　そこで、上記の本発明の装置は、更に、光検出部が光検出領域からの光の少なくとも二つの互いに異なる成分を別々に検出し、信号処理部が、前記の成分の各々の時系列の光強度データを生成し、更に、信号処理部が、成分の各々の第一及び第二の発生確率を算出できるよう構成されていてよい。かかる構成に於いては、検出成分を適宜選択することにより、観測対象となる単一粒子の有する所定の特性値を検出された光の成分の各々の第二の発生確率に反映させること、即ち、第二の発生確率を単一粒子の有する所定の特性値の関数とすることが可能となる。そして、成分毎の第一及び第二の発生確率の発生確率を参照して、所定の特性値の関数である第二の発生確率が相対的に高い領域に於いて、所定の特性値を有する粒子が存在したと判定できることとなる（所定の特性値を有していない粒子が存在していた場合、所定の特性値の関数である第二の発生確率は低くなる。）。即ち、成分毎の第一及び第二の発生確率の発生確率を用いることにより、単に光検出領域内の粒子の存在の有無だけでなく、時系列の光強度データ上にて光検出領域内に所定の特性値を有する粒子が存在している時間が決定され、所定の特性値を有する粒子の存在を検出することが可能となる。
　更に、上記の如き、少なくとも二つの互いに異なる成分毎の第一及び第二の発生確率の発生確率を用いる態様によれば、単一粒子として互いに異なる所定の特性値を有する複数の種類の単一粒子が含まれる場合、それぞれの種類の単一粒子について、時系列の光強度データ上にて粒子の存在した時間が決定できることとなる。即ち、単一粒子の種類毎に、それらの互いに異なる所定の特性値の関数である成分の各々の第二の発生確率を算出し（第一の状態では粒子が存在していないので、第一の発生確率には、粒子の特性値は反映されない。）、複数の成分の各々の第一の発生確率と複数の種類の単一粒子の各々の成分の各々の第二の発生確率とを参照したとき、相対的に高い値の第二の発生確率を与える種類の粒子が光検出領域内に存在していると推定することができる。かくして、単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含む場合、上記の本発明の装置は、単一粒子の種類毎に、互いに異なる所定の特性値の関数である成分の各々の第二の発生確率を算出し、成分の各々の第一の発生確率と複数の種類の単一粒子の各々の成分の各々の第二の発生確率とに基づいて時系列の光強度データ上にて単一粒子の種類毎に該単一粒子の存在を表す信号を検出するよう構成されていてよい。要すれば、かかる構成により、試料溶液中に複数の種類の単一粒子が含まれている場合に、種類を識別しながら、粒子の検出が可能となる。
　なお、上記の第二の発生確率に反映させるべき粒子の特性値としては、既に触れた如く、単一粒子の偏光特性を表す指標値、例えば、蛍光異方性や、単一粒子の発光波長特性を表す指標値、例えば、互いに異なる発光波長帯域の発光強度の比などであってよい。
　上記の本発明の装置に於ける試料溶液内での光検出領域の位置の移動速度は、観測対象の単一粒子の特性又は試料溶液中の数密度又は濃度に基づいて適宜変更可能となっていてよい。光検出領域の移動速度が速くなると、単一粒子が発光粒子である場合には、一つの発光粒子から得られる光量は低減することとなり、また、単一粒子が非発光粒子である場合には、一つの非発光粒子の存在による光強度値の低減量が小さくなる。従って、光検出領域の移動速度は、単一粒子による光強度値の変化が精度よく又は感度よく計測できるように、適宜変更可能となっていることが好ましい。また、試料溶液内での光検出領域の位置の移動速度は、好適には、検出対象となる単一粒子の拡散移動速度（ブラウン運動による粒子の平均の移動速度）よりも高く設定される。上記に説明されている如く、本発明の装置では、検出された光強度値の時間変化のパターンを、光検出領域内に粒子が無い場合と有る場合（第一及び第二の状態）を仮定して算出された時間変化パターンの（第一及び第二の）発生確率を用いて評価するところ、粒子の光検出領域内の通過時に粒子がブラウン運動により移動することも考慮すると、第二の発生確率の算出が複雑となる。そこで、本発明では、粒子の光検出領域内の通過中に於ける粒子のブラウン運動による効果が無視できるように、光検出領域の移動速度が検出対象となる単一粒子の拡散移動速度よりも速く設定されていることが好ましい。なお、拡散移動速度は、単一粒子によって変わるので、上記の如く、その特性（特に、拡散係数）に応じて、本発明の装置は、光検出領域の移動速度が適宜変更可能であるよう構成されていることが好ましい。
　試料溶液内での光検出領域の位置の移動は、任意の方式で為されてよい。例えば、レーザー走査型光学顕微鏡に於いて採用されているガルバノミラーを用いるなどして、顕微鏡の光学系の光路を変更して光検出領域の位置が変更されるようになっていてよく、或いは、試料溶液の位置を（例えば、顕微鏡のステージを移動するなどして）動かして、光検出領域の試料溶液内に於ける位置を移動するようになっていてよい。光検出領域の位置の移動軌跡は、任意に設定されてよく、例えば、円形、楕円形、矩形、直線及び曲線のうちから選択可能であってよい。特に、顕微鏡の光学系の光路を変更して光検出領域の位置が変更される場合、光検出領域の移動は、速やかであり、且つ、試料溶液に於いて機械的振動や流体力学的な作用が実質的に発生しないので、検出対象となる単一粒子が力学的な作用の影響を受けることなく安定した状態にて、光の計測が可能である点で有利である。
　上記の本発明の実施の態様の一つに於いては、信号の数を計数して光検出領域に包含された単一粒子の数を計数するようになっていてよい（粒子のカウンティング）。その場合、検出された単一粒子の数と光検出領域の位置の移動量と組み合わせることにより、試料溶液中の同定された単一粒子の数密度又は濃度に関する情報が得られることとなる。具体的には、例えば、複数の試料溶液の数密度若しくは濃度の比、或いは、濃度若しくは数密度の基準となる標準試料溶液に対する相対的な数密度若しくは濃度の比が算出されるか、又は、濃度若しくは数密度の基準となる標準試料溶液に対する相対的な数密度若しくは濃度の比を用いて、絶対的な数密度値又は濃度値が決定されてよい。或いは、任意の手法により、例えば、所定の速度にて光検出領域の位置を移動するなどして、光検出領域の位置の移動軌跡の全体積を特定すれば、単一粒子の数密度又は濃度が具体的に算定できることとなる。
　上記の本発明の装置に於いて試料溶液内に於ける光検出領域の位置を移動させながら光検出を行い、単一粒子からの信号を個別に検出する光分析技術であって、検出された光強度値の時間変化のパターンが発生する確率（第一及び第二の発生確率）を参照して、単一粒子の有無を決定する光分析技術の処理は、汎用のコンピュータによっても実現可能である。
　従って、本発明のもう一つの態様によれば、共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出するための光分析用コンピュータプログラムであって、試料溶液内に於いて顕微鏡の光学系の光検出領域の位置を移動する手順と、試料溶液内に於いて光検出領域の位置を移動させながら光検出領域からの光の強度を測定して光強度データを生成する手順と、時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度の時間変化値の、光検出領域内に単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、光検出領域内に単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出する手順と、第一及び第二の発生確率に基づいて時系列の光強度データ上にて単一粒子の各々の存在を表す信号を検出する手順とをコンピュータに実行させることを特徴とするコンピュータプログラムが提供される。この場合に於いても、本発明の装置と同様に、典型的には、光検出領域からの光を検出して時系列の光強度データを生成する手順に於いて、計測単位時間（ビンタイム）毎に到来する光子数を計数するフォトンカウンティングにより光検出領域からの光が検出され、時系列の光強度データが時系列のフォトンカウントデータである。「解析窓」は、本発明の装置の場合と同様に設定されてよい。観測対象となる粒子が、発光粒子である場合には、時系列光強度データ上に於ける「光強度値の一時的な増大」が「単一粒子の各々の存在を表す信号」となる。また、観測対象となる粒子が、（検出波長帯域に於いて）非発光粒子である場合には、光検出領域からの光が有意な背景光を含み、背景光からの光強度の一時的な低減が「単一粒子の各々の存在を表す信号」となる。なお、コンピュータプログラムは、コンピュータ読み取り可能な記憶媒体に記憶されて提供される。コンピュータは、記憶媒体に記憶されているプログラムを読み出して、情報の加工・演算処理を実行することにより、上記の手順を実現する。ここで、コンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭ、半導体メモリ等であってよい。更に、上述のプログラムは、通信回線によってコンピュータへ配信され、この配信を受けたコンピュータがプログラムを実行するようにしても良い。
　上記の構成に於いても、典型的には、第一及び第二の発生確率は、それぞれ、解析窓内に於ける単位時間毎の光強度値と第一及び第二の状態を仮定した場合の単位時間毎の期待値とに基づいて算出され、より具体的には、単位時間毎の光強度値が単位時間毎の期待値を有するポアソン分布に従うものとして単位時間毎の光強度値の単位時間発生確率が算出され、第一及び第二の発生確率がそれぞれ対応する単位時間発生確率を用いて算出されてよい。そして、第二の発生確率が第一の発生確率より大きい解析窓の時間に光検出領域内に単一粒子が存在したと判定されるか、第二の発生確率と第一の発生確率とのオッズ比に基づいて、光検出領域内に単一粒子が存在した時間領域が決定されてよい。
　また、上記のコンピュータプログラムも、光検出領域からの光の少なくとも二つの互いに異なる成分が別々に検出されて、かかる成分の各々の時系列の光強度データが生成され、更に、成分の各々について、第一の発生確率と、観測対象である単一粒子の所定の特性値の関数である第二の発生確率が算出され、成分毎の第一及び第二の発生確率の発生確率に基づいて時系列の光強度データ上にて所定の特性値を有する単一粒子の存在を表す信号が検出されるよう構成されていてよい。そして、単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含む場合には、上記のコンピュータプログラムに於いて、単一粒子の種類毎に、互いに異なる所定の特性値の関数である成分の各々の第二の発生確率が算出され、成分の各々の第一の発生確率と複数の種類の単一粒子の各々の成分の各々の第二の発生確率とに基づいて時系列の光強度データ上にて単一粒子の種類毎に該単一粒子の存在を表す信号が検出されてよい。
　更に、試料溶液内での光検出領域の位置の移動速度は、単一粒子の特性、試料溶液中の数密度又は濃度に基づいて適宜変更可能となっていてよく、好適には、試料溶液内での光検出領域の位置の移動速度は、検出対象となる単一粒子の拡散移動速度よりも高く設定される。試料溶液内での光検出領域の位置の移動は、任意の方式で為されてよく、好適には、顕微鏡の光学系の光路を変更して或いは試料溶液の位置を動かして光検出領域の位置が変更されるようになっていてよい。光検出領域の位置の移動軌跡は、任意に設定されてよく、例えば、円形、楕円形、矩形、直線及び曲線のうちから選択可能であってよい。
　更に、上記のコンピュータプログラムに於いても、個別に検出された単一粒子からの信号の数を計数して光検出領域の位置の移動中に検出された単一粒子の数を計数する手順及び／又は検出された単一粒子の数に基づいて、試料溶液中の単一粒子の数密度又は濃度を決定する手順が含まれていてよい。
　上記の本発明の装置又はコンピュータプログラムによれば、試料溶液内に於ける光検出領域の位置を移動させながら、個々の発光粒子の光の検出を行う光分析方法であって、検出された光強度値の時間変化のパターンが発生する確率（第一及び第二の発生確率）を参照して、単一粒子の有無を決定する新規な方法が実現される。
　かくして、本発明によれば、更に、共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出する光分析方法であって、試料溶液内に於いて顕微鏡の光学系の光検出領域の位置を移動する過程と、試料溶液内に於いて光検出領域の位置を移動させながら光検出領域からの光の強度を測定して光強度データを生成する過程と、時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、光検出領域内に単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、光検出領域内に単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出する過程と、第一及び第二の発生確率に基づいて時系列の光強度データ上にて単一粒子の各々の存在を表す信号を検出する過程とを含むことを特徴とする方法が提供される。この場合に於いても、本発明の装置と同様に、典型的には、光検出領域からの光を検出して時系列の光強度データを生成する過程に於いて、計測単位時間（ビンタイム）毎に到来する光子数を計数するフォトンカウンティングにより光検出領域からの光が検出され、時系列の光強度データが時系列のフォトンカウントデータである。「解析窓」は、本発明の装置の場合と同様に設定されてよい。観測対象となる粒子が、発光粒子である場合には、時系列光強度データ上に於ける「光強度値の一時的な増大」が「単一粒子の各々の存在を表す信号」となる。また、観測対象となる粒子が、（検出波長帯域に於いて）非発光粒子である場合には、光検出領域からの光が有意な背景光を含み、背景光からの光強度の一時的な低減が「単一粒子の各々の存在を表す信号」となる。
　上記の構成に於いても、典型的には、第一及び第二の発生確率は、それぞれ、解析窓内に於ける単位時間毎の光強度値と第一及び第二の状態を仮定した場合の単位時間毎の期待値とに基づいて算出され、より具体的には、単位時間毎の光強度値が単位時間毎の期待値を有するポアソン分布に従うものとして単位時間毎の光強度値の単位時間発生確率が算出され、第一及び第二の発生確率がそれぞれ対応する単位時間発生確率を用いて算出されてよい。そして、第二の発生確率が第一の発生確率より大きい解析窓の時間に光検出領域内に単一粒子が存在したと判定されるか、第二の発生確率と第一の発生確率とのオッズ比に基づいて、光検出領域内に単一粒子が存在した時間領域が決定されてよい。
　また、上記の方法も、光検出領域からの光の少なくとも二つの互いに異なる成分が別々に検出されて、かかる成分の各々の時系列の光強度データが生成され、更に、成分の各々について、第一の発生確率と、観測対象である単一粒子の所定の特性値の関数である第二の発生確率が算出され、成分毎の第一及び第二の発生確率の発生確率に基づいて時系列の光強度データ上にて所定の特性値を有する単一粒子の存在を表す信号が検出されるよう構成されていてよい。そして、単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含む場合には、上記の方法に於いて、単一粒子の種類毎に、互いに異なる所定の特性値の関数である成分の各々の第二の発生確率が算出され、成分の各々の第一の発生確率と複数の種類の単一粒子の各々の成分の各々の第二の発生確率とに基づいて時系列の光強度データ上にて単一粒子の種類毎に該単一粒子の存在を表す信号が検出されてよい。
　更に、試料溶液内での光検出領域の位置の移動速度は、単一粒子の特性、試料溶液中の数密度又は濃度に基づいて適宜変更可能となっていてよく、好適には、試料溶液内での光検出領域の位置の移動速度は、検出対象となる単一粒子の拡散移動速度よりも高く設定される。試料溶液内での光検出領域の位置の移動は、任意の方式で為されてよく、好適には、顕微鏡の光学系の光路を変更して或いは試料溶液の位置を動かして光検出領域の位置が変更されるようになっていてよい。光検出領域の位置の移動軌跡は、任意に設定されてよく、例えば、円形、楕円形、矩形、直線及び曲線のうちから選択可能であってよい。
　更に、上記の方法に於いても、個別に検出された単一粒子からの信号の数を計数して光検出領域の位置の移動中に検出された単一粒子の数を計数する過程及び／又は検出された単一粒子の数に基づいて、試料溶液中の単一粒子の数密度又は濃度を決定する過程が含まれていてよい。
　上記の本発明の光分析技術は、典型的には、タンパク質、ペプチド、核酸、脂質、糖鎖、アミノ酸若しくはこれらの凝集体などの生体分子、ウイルス、細胞などの粒子状の生物学的な対象物の溶液中の状態の分析又は解析の用途に用いられるが、非生物学的な粒子（例えば、原子、分子、ミセル、金属コロイドなど）の溶液中の状態の分析又は解析に用いられてもよく、そのような場合も本発明の範囲に属することは理解されるべきである。
　かくして、本発明によれば、走査分子計数法に於いて、単に、時系列光強度データ上の光強度値の増減を参照するのではなく、光強度値の時間変化のパターンが発生しやすい状態を決定して光検出領域内に単一粒子が存在しているか否かが推定される。かかる構成によれば、単一粒子による光強度値の変化が比較的小さく、光強度値変化の絶対値のみでは、粒子による信号かノイズ信号かが判断が難しい場合でも、より精度良く、粒子の信号とノイズ信号との識別が可能となることが期待される。また、粒子の信号とノイズ信号との識別の精度が向上されることにより、走査分子計数法に於いて検出可能な試料溶液中の単一粒子の濃度範囲の低濃度側への拡大が期待される。更に、観測対象となる単一粒子の所定の特性値を考慮して算出された光強度値の時間変化の発生確率を用いる場合には、観測対象となる単一粒子とそれ以外の粒子との識別も可能となり、不純物等が混在する試料溶液に於いても計測の精度又は感度の向上が期待される。
　本発明のその他の目的及び利点は、以下の本発明の好ましい実施形態の説明により明らかになるであろう。
図１（Ａ）は、本発明による走査分子計数法を実行する光分析装置の内部構造の模式図である。図１（Ｂ）は、コンフォーカル・ボリューム（共焦点顕微鏡の観察領域）の模式図である。図１（Ｃ）は、ミラー７の向きを変更して試料溶液内に於いて光検出領域の位置を移動する機構の模式図である。図１（Ｄ）は、マイクロプレートの水平方向位置を移動して試料溶液内に於ける光検出領域の位置を移動する機構の模式図である。 図２（Ａ）、（Ｂ）は、それぞれ、本発明が適用される走査分子計数法に於ける発光粒子の光を検出の原理を説明する模式図及び計測される光強度の時間変化の模式図である。図２（Ｃ）、（Ｄ）は、それぞれ、本発明が適用される反転型走査分子計数法に於ける検出光波長帯域で発光しない単一粒子の存在を検出する原理を説明する模式図及び計測される光強度の時間変化の模式図である。 図３（Ａ）は、時系列光強度データの模式的な例を示している。（左）輝度が大きい発光粒子が存在している場合。（中）輝度が小さい発光粒子が存在している場合。（右）発光粒子が存在していない場合。図３（Ｂ）は、本発明に於いて時系列光強度データ上に設定される解析窓を説明する図である。図３（Ｃ）は、試料溶液内の光検出領域の位置を発光粒子の拡散移動速度よりも速い速度にて移動することにより発光粒子が光検出領域を横切る場合の粒子の運動の態様を表すモデル図である。図３（Ｄ）は、光検出領域に於ける発光粒子から放出され検出される光の光検出領域の半径ｒ方向の強度分布の模式図である。 図４（Ａ）は、説明の目的で模式的に表した時系列光強度データであり、上段は、発光粒子が存在している状態（左：第二の状態）及び発光粒子が存在していない状態（右：第一の状態）をそれぞれ仮定した場合の時系列光強度データ上の光強度値の期待値を表しており、中段は、実際に計測された光強度値を表しており、下段は、発光粒子が存在している状態（左）及び発光粒子が存在していない状態（右）をそれぞれ仮定した場合の計測された光強度値の発生確率を示している。図４（Ｂ）は、発光粒子から放出される光の互いに偏光成分のモデル図（上）と、各成分の光強度値の平均的な時間変化（左下）と、合計の光強度値の平均的な時間変化（右下）とを示している。図４（Ｃ）は、互いに異なる発光波長特性を有する発光粒子Ｄ１、Ｄ２の発光波長スペクトルであり、互いに異なる波長帯域の光の成分を別々に検出した場合に於いて、各波長帯域にて検出される各発光粒子からの成分の光強度値を示している。陰影部分の面積が各成分の光強度値に相当する。 本発明に従って実行される走査分子計数法の処理手順をフローチャートの形式で表した図である。 図６は、本発明による走査分子計数法に従って得られた時系列光強度データ（フォトンカウントデータ）の一部（Ａ）と、（Ａ）の時系列光強度データをスムージングして得られた平滑化された時系列光強度データ（Ｂ）と、（Ａ）の時系列光強度データの、本発明の教示に従って粒子が存在する場合と存在しない場合の各々を仮定して算出された発生確率（第二及び第一の発生確率）のオッズ比（Ｃ）を示している。 図７は、図６と同様の図であって、時系列光強度データの別の部分を示している。 図８は、本発明による走査分子計数法に従って発光粒子の検出及び計数を行った場合の試料溶液中の発光粒子（ＡＴＴＯ６４７Ｎ）濃度と検出された粒子数との関係を示した図である。 図９（Ａ）は、本発明による反転型走査分子計数法に従って得られた時系列光強度データ（フォトンカウントデータ）の一部（上段）と、その時系列光強度データの、本発明の教示に従って粒子が存在する場合と存在しない場合の各々を仮定して算出された発生確率（第二及び第一の発生確率）のオッズ比（下段）を示している。図９（Ｂ）は、（Ａ）の一部を拡大した図である。 図１０は、図９（Ｂ）を更に拡大した図である。 図１１（Ａ）は、本発明による反転型走査分子計数法に従って粒子の検出及び計数を行った場合の試料溶液中の粒子濃度と検出された粒子数との関係を示した図であり、図１１（Ｂ）は、（Ａ）と同一の時系列光強度データに於いて、試料溶液中の粒子濃度と平滑化及び釣鐘型関数のフィッティングにより検出された粒子数との関係を示した図である。図中、エラーバーは、標準偏差である。 図１２は、本発明による走査分子計数法に従って光の計測を互いに異なる偏光方向の成分を別々に検出して得られた時系列光強度データに於いて、発光粒子の蛍光偏光異方性を考慮して算出された発生確率を用いて、発光粒子の検出及び発光粒子の種類の識別を行ったときの結果を示している。図中、左は、プラスミドを含む溶液に於いて検出された発光粒子の種類の判別結果であり、右は、蛍光色素（ＴＡＭＲＡ）を含む溶液に於いて検出された発光粒子の種類の判別結果である。
１…光分析装置（共焦点顕微鏡）
２…光源
３…シングルモードオプティカルファイバー
４…コリメータレンズ
５…ダイクロイックミラー
６、７、１１…反射ミラー
８…対物レンズ
９…マイクロプレート
１０…ウェル（試料溶液容器）
１２…コンデンサーレンズ
１３…ピンホール
１４…バリアフィルター
１４ａ…ダイクロイックミラー又はビームスプリッター
１５…マルチモードオプティカルファイバー
１６…光検出器
１７…ミラー偏向器
１７ａ…ステージ位置変更装置
１８…コンピュータ
　以下、本発明の好ましい実施形態について詳細に説明する。
光分析装置の構成
　本発明による光分析技術を実現する光分析装置は、基本的な構成に於いて、図１（Ａ）に模式的に例示されている如き、ＦＣＳ、ＦＩＤＡ等が実行可能な共焦点顕微鏡の光学系と光検出器とを組み合わせてなる装置であってよい。同図を参照して、光分析装置１は、光学系２～１７と、光学系の各部の作動を制御すると共にデータを取得し解析するためのコンピュータ１８とから構成される。光分析装置１の光学系は、通常の共焦点顕微鏡の光学系と同様であってよく、そこに於いて、光源２から放射されシングルモードファイバー３内を伝播したレーザー光（Ｅｘ）が、ファイバーの出射端に於いて固有のＮＡにて決まった角度にて発散する光となって放射され、コリメーター４によって平行光となり、ダイクロイックミラー５、反射ミラー６、７にて反射され、対物レンズ８へ入射される。対物レンズ８の上方には、典型的には、１～数十μＬの試料溶液が分注される試料容器又はウェル１０が配列されたマイクロプレート９が配置されており、対物レンズ８から出射したレーザー光は、試料容器又はウェル１０内の試料溶液中で焦点を結び、光強度の強い領域（励起領域）が形成される。
　観測対象となる単一粒子が発光粒子である場合、試料溶液中には、発光粒子、典型的には、蛍光性粒子又は蛍光色素等の発光標識が付加された粒子が分散又は溶解され、かかる発光粒子が励起領域に進入すると、その間、発光粒子が励起され光が放出される。一方、観測対象となる単一粒子が非発光粒子である場合、試料溶液中には、典型的には、検出波長帯域にて発光しない粒子と、背景光を生ずる任意の発光物質が分散又は溶解され、検出波長帯域にて発光しない粒子が励起領域に進入していないときには、発光物質が励起されて実質的に一定の光が放出されて、背景光となり、検出波長帯域にて発光しない粒子が励起領域に進入すると、背景光が低減することとなる。
　かくして、励起領域から放出された光（Ｅｍ）は、対物レンズ８、ダイクロイックミラー５を通過し、ミラー１１にて反射してコンデンサーレンズ１２にて集光され、ピンホール１３を通過し、バリアフィルター１４を透過して（ここで、特定の波長帯域の光成分のみが選択される。）、マルチモードファイバー１５に導入されて、光検出器１６に到達し、時系列の電気信号に変換された後、コンピュータ１８へ入力され、後に説明される態様にて光分析のための処理が為される。なお、当業者に於いて知られている如く、上記の構成に於いて、ピンホール１３は、対物レンズ８の焦点位置と共役の位置に配置されており、これにより、図１（Ｂ）に模式的に示されている如きレーザー光の焦点領域、即ち、光検出領域内から発せられた光のみがピンホール１３を通過し、光検出領域以外からの光は遮断される。図１（Ｂ）に例示されたレーザー光の焦点領域は、通常、１～１０ｆＬ程度の実効体積を有する本光分析装置に於ける光検出領域であり（典型的には、光強度が領域の中心を頂点とするガウス様分布となる。実効体積は、光強度が中心光強度の１／ｅ２となる面を境界とする略楕円球体の体積である。）、コンフォーカル・ボリュームと称される。また、本発明では、１つの発光粒子からの光、例えば、一つの蛍光色素分子からの微弱光又は非発光粒子の存在による背景光の低減が検出されるので、光検出器１６としては、好適には、フォトンカウンティングに使用可能な超高感度の光検出器が用いられる。光の検出がフォトンカウンティングによる場合、光強度の測定は、所定時間に亘って、逐次的に、計測単位時間（ＢＩＮ　ＴＩＭＥ）毎に、光検出器に到来するフォトンの数を計測する態様にて実行される。従って、この場合、時系列の光強度のデータは、時系列のフォトンカウントデータである。また、顕微鏡のステージ（図示せず）には、観察するべきウェル１０を変更するべく、マイクロプレート９の水平方向位置を移動するためのステージ位置変更装置１７ａが設けられていてよい。ステージ位置変更装置１７ａの作動は、コンピュータ１８により制御されてよい。かかる構成により、検体が複数在る場合にも、迅速な計測が達成可能となる。
　更に、上記の光分析装置の光学系に於いては、光検出領域により試料溶液内を走査する、即ち、試料溶液内に於いて焦点領域、即ち、光検出領域の位置を移動するための機構が設けられる。かかる光検出領域の位置を移動するための機構としては、例えば、図１（Ｃ）に模式的に例示されている如く、反射ミラー７の向きを変更するミラー偏向器１７が採用されてよい（光検出領域の絶対的な位置を移動する方式）。かかるミラー偏向器１７は、通常のレーザー走査型顕微鏡に装備されているガルバノミラー装置と同様であってよい。或いは、別の態様として、図１（Ｄ）に例示されている如く、試料溶液が注入されている容器１０（マイクロプレート９）の水平方向の位置を移動し、試料溶液内に於ける光検出領域の相対的な位置を移動するべくステージ位置変更装置１７ａが作動されてもよい（試料溶液の絶対的な位置を移動する方式）。いずれの方式による場合も、所望の光検出領域の位置の移動パターンを達成するべく、ミラー偏向器１７又はステージ位置変更装置１７ａは、コンピュータ１８の制御の下、光検出器１６による光検出と協調して駆動される。光検出領域の位置の移動軌跡は、円形、楕円形、矩形、直線、曲線又はこれらの組み合わせから任意に選択されてよい（コンピュータ１８に於けるプログラムに於いて、種々の移動パターンが選択できるようになっていてよい。）。なお、図示していないが、対物レンズ８又はステージを上下に移動することにより、光検出領域の位置が上下方向に移動されるようになっていてもよい。
　観測対象物となる発光粒子又は背景光を与える物質が多光子吸収により発光する場合には、上記の光学系は、多光子顕微鏡として使用される。その場合には、励起光の焦点領域（光検出領域）のみで光の放出があるので、ピンホール１３は、除去されてよい。また、発光粒子又は背景光を与える物質が化学発光や生物発光現象により励起光によらず発光する場合には、励起光を生成するための光学系２～５が省略されてよい。発光粒子又は背景光を与える物質がりん光又は散乱により発光する場合には、上記の共焦点顕微鏡の光学系がそのまま用いられる。更に、光分析装置１に於いては、図示の如く、複数の励起光源２が設けられていてよく、発光粒子又は背景光を与える物質の励起波長によって適宜、励起光の波長が選択できるようになっていてよい。
　また更に、光検出器１６も複数個備えられ、各々の光検出器１６が光検出領域からの光のうちの互いに異なる成分を別々に検出するよう構成されていてよい。後により詳細に述べる如く、検出する成分を適宜選択することにより、特定の発光特性を有する単一粒子を選択的に検出することが可能となる。そのような光検出領域からの光のうちの互いに異なる成分を検出する場合には、ピンホール１３の通過後の検出光光路に於いては、光路を任意の態様にて分割するための機構が設けられる。例えば、光検出領域からの光を互いに異なる偏光成分に分けて検出する場合には、検出光光路の符号１４ａが付された部位に偏光ビームスプリッタ１４ａが挿入される。この場合、励起光光路には、ポーラライザ（図示せず）が挿入される。また、検出光光路の部位１４ａに、検出光光路に特定の波長帯域の光を反射し、別の波長帯域を透過するダイクロイックミラー１４ａが挿入されることにより、互いに異なる波長帯域の光成分が別々に検出可能となる。
　コンピュータ１８は、ＣＰＵおよびメモリを備え、ＣＰＵが各種演算処理を実行することにより、本発明の手順を実行する。なお、各手順は、ハードウェアにより構成するようにしてもよい。本実施形態で説明される処理の全て或いは一部は、それらの処理を実現するプログラムを記憶したコンピュータ読み取り可能な記憶媒体を用いて、コンピュータ１８により実行されてよい。即ち、コンピュータ１８は、記憶媒体に記憶されているプログラムを読み出して、情報の加工・演算処理を実行することにより、本発明の処理手順を実現するようになっていてよい。ここで、コンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭ、半導体メモリ等であってよく、或いは、上記のプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータがプログラムを実行するようにしても良い。
本発明の光分析技術の原理
　「発明の概要」の欄に記載されている如く、本発明の光分析技術に於いては、端的に述べれば、走査分子計数法又は反転型走査分子計数法に於いて、光検出領域の位置を移動しながら計測して得られた時系列光強度データから観測対象となる単一粒子の信号を検出する際に、単に、時系列光強度データ上の光強度値の増減を参照するのではなく、計測された光強度値の時間変化のパターンが、光検出領域内に粒子が存在する場合と存在しない場合とのいずれに於いて、より高い確率で発生するパターンであるかが推定される。そして、その推定結果に基づいて、時系列光強度データ上で観測対象となる単一粒子の信号の存在の有無及び／又はその数が検出される。以下、本発明の走査分子計数法及び観測対象となる単一粒子の信号の検出の原理について説明する。
１.走査分子計数法の原理
　走査分子計数法に於いて実行される基本的な処理に於いては、特許文献９～１１に記載されている如く、端的に述べれば、光検出領域の位置を移動するための機構（ミラー偏向器１７）を駆動して光路を変更し、或いは、試料溶液が注入されている容器１０（マイクロプレート９）の水平方向の位置を移動して、図２（Ａ）にて模式的に描かれているように、試料溶液内に於いて光検出領域ＣＶの位置を移動しながら、即ち、光検出領域ＣＶにより試料溶液内を走査しながら、光検出が実行される。そうすると、観測対象粒子が発光粒子である場合、光検出領域ＣＶが移動する間（図中、時間ｔｏ～ｔ２）に於いて１つの発光粒子の存在する領域を通過する際（ｔ１）には、発光粒子から光が放出され、図２（Ｂ）に描かれている如き時系列の光強度データ上に有意な光強度（Ｅｍ）のパルス状の信号が出現することとなる。そこで、上記の光検出領域ＣＶの位置の移動と光検出を実行し、その間に出現する図２（Ｂ）に例示されている如きパルス状の信号（有意な光強度）を一つずつ検出することによって、発光粒子が個別に検出され、その数をカウントすることにより、計測された領域内に存在する発光粒子の数、或いは、濃度若しくは数密度に関する情報が取得できることとなる。
　また、観測対象となる粒子が非発光粒子である場合、上記の走査分子計数法による光計測に於いて、光検出領域から背景光を放出させ（或いは、光検出領域を照明光にて照明し）、観測対象となる粒子が光検出領域に進入した際の検出される背景光の低下を捉えることにより、非発光粒子の存在を検出することが可能となる（反転型走査分子計数法）。この場合、より具体的には、図２（Ｃ）にて模式的に描かれているように、試料溶液内に於いて光検出領域ＣＶの位置を移動しながら、光検出が実行される。ここで、試料溶液中には発光物質が分散させておくと、光検出領域ＣＶ内には、多数の発光物質が存在し、光検出領域ＣＶの移動中（図中、時間ｔｏ～ｔ２）、それらの発光物質からの光が略一様に検出される。しかしながら、光検出領域ＣＶが移動する間に単一の非発光粒子の存在する領域を通過する際（ｔ１）には、発光物質の存在領域の体積が低減するので、発光物質の放出する光の総量が低減し、図２（Ｄ）に描かれている如く、時系列の光強度データ上に、釣鐘型のパルス状の有意な光強度（Ｅｍ）の低下が出現することとなる。かくして、上記の光検出領域ＣＶの位置の移動と光検出を実行し、その間に出現する図２（Ｄ）に例示されている如きパルス状の有意な光強度の低下、即ち、非発光粒子の存在を表す信号を一つずつ検出することによって、非発光粒子が個別に検出され、その数をカウントすることにより、計測された領域内に存在する単一粒子の数、或いは、濃度若しくは数密度に関する情報が取得できることとなる。
　ところで、上記の如き走査分子計数法に於いて、実際に計測された時系列光強度データ上には、発光粒子の信号の他に、熱ノイズ、迷光、水のラマン散乱等による光強度値の増大（ノイズ信号）が存在する。また、反転型走査分子計数法の場合には、背景光の強度のゆらぎによる光強度値の低減（この場合も、以下、ノイズ信号と称する。）が存在する。そこで、時系列光強度データ上に於いて単一粒子の信号の検出を行う際には、まず、パルス状の信号を検出し、しかる後に、検出されたパルス状の信号の強度、時間幅及び形状が検査され、単一粒子が光検出領域を通過した場合に得られるべき信号の強度、時間幅及び形状の条件に適合している信号のみが、単一粒子の信号として判定され、それ以外の信号はノイズ信号として判定される。
　しかしながら、単一粒子の信号の強度値が弱くなると、信号の強度、時間幅及び形状に基づく単一粒子の信号とノイズ信号との識別は、困難となる。特に、観測対象粒子が発光粒子である場合、発光粒子から放出される光子数が微量であり、且つ、光子は確率的に放出されるので、光強度値のプロファイルは、実際には、図２（Ｂ）に描かれている如き滑らかな釣鐘状にはならず、図３（Ａ）左にて例示されている如く、離散的になる。従って、発光粒子から放出される光子数が少なくなると、その場合の光強度データ上の信号の強度、時間幅及び形状（図３（Ａ）中）は、粒子が存在せずノイズ信号のみ存在する光強度データ（図３（Ａ）右）の場合と区別がつき難くなる。また、光検出領域に於ける発光粒子から放出され検出される光の強度値は、発光粒子の位置が光検出領域の略中心から離れるほど低下し（図３（Ｄ）参照）、且つ、光検出領域の略中心から離れるほど発光粒子の数は多くなるので、ノイズ信号との識別が困難な粒子の信号の数が大きくなる。同様の状況は、観測対象粒子が非発光粒子である反転型走査分子計数法に於ける背景光の光子数の低下の場合にも生ずる。
　走査分子計数法及び反転型走査分子計数法に於いて、その精度及び／又は感度は、光検出領域の或る移動距離当たりの粒子信号の検出精度と検出数が高いほど、向上される。しかしながら、上記の如く、粒子信号の検出数を多くすべく、信号の強度の判定基準値を低減すると、ノイズ信号を誤って粒子信号と判定する可能性が増大し、粒子信号の検出精度を向上すべく、信号の強度の判定基準値を増大すると、多くの輝度の低い粒子信号が検出されないこととなる。そこで、本発明に於いては、信号の強度、時間幅及び形状によってはノイズ信号との識別が困難な単一粒子の信号の検出がより精度良く達成可能とする新規なアルゴリズムが提案される。
２．単一粒子の信号の検出
　時系列光強度データがフォトンカウントデータである場合に於いて、光強度値は、ビンタイム当たりの検出光子数となる。従って、図３（Ａ）に模式的に例示されているように、光強度値は、時間軸方向に離散的に分布する。その場合、輝度が大きい粒子信号は、ノイズ信号の影響が少なく、略釣鐘状のプロファイル（左図）を有するが、輝度の小さい粒子信号は、強度値がノイズ信号と略同様となり、更に、ノイズ信号が重畳するので（中図）、略釣鐘状のプロファイルを抽出することが困難となり、粒子のないノイズ信号のみが存在する時間領域（右図）との識別が困難となる。しかしながら、輝度の小さい粒子信号の存在する時間領域とノイズ信号のみが存在する時間領域とでは、光子が検出される事象（光子数が１以上となる事象）の発生頻度とパターンに違いを有する。即ち、図からも理解される如く、ノイズ信号の場合は、常にランダムに光子検出の事象が生ずるのに対し、粒子信号の場合は、光子検出の事象が時間的に集中し、特に、その中央付近の強度値が高くなる傾向を有する。同様の現象は、非発光粒子の信号と背景光のゆらぎとの間にも観察される。
　かくして、本発明では、かかる光子検出の事象の発生頻度とパターンの違いに着目して、粒子信号を選択的に検出することが試みられる。そのために、本発明による新規な単一粒子信号の検出のためのアルゴリズムに於いては、時系列光強度データ上の所定の時間幅毎の光強度値の時間変化（光子数列）について、光検出領域内に粒子が存在していたと仮定した場合に前記の光強度値の時間変化が発生する確率（発生確率（粒子存在時）、第二の発生確率）と、光検出領域内に粒子が存在していなかったと仮定した場合に前記の光強度値の時間変化が発生する確率（発生確率（粒子不在時）、第一の発生確率）とが算出される。そして、高い発生確率を与える状態が実際の状態であると推定される。
　より具体的には、まず、図３（Ｂ）を参照して、時系列光強度データに於いて、任意の時間幅の区間（以下、解析窓と称する。）が設定される。この解析窓は、単位時間（通常は、ビンタイムであってよい。）ｔｉ（ｉ＝１，２，…　以下、同様）毎に、検出された光子数Ｃｉを有する。ところで、各単位時間内に於いて生ずる光子検出の事象の数は、その単位時間に於ける期待値を有するポアソン分布に従うと考えられるので、任意の単位時間ｔｉに於いて、検出光子数Ｃｉとなる確率（単位時間発生確率）は、
Figure JPOXMLDOC01-appb-M000001
 により与えられる。ここで、Ｅｉは、単位時間ｔｉに於ける光子数の期待値である。そして、解析窓にｎ＋１個の単位時間が含まれるとき、解析窓内に於ける検出光子数列Ｃｉが発生する確率Ｐ（発生確率）は、
Figure JPOXMLDOC01-appb-M000002
 により与えられる。
　上記の各単位時間ｔｉに於ける光子検出事象の発生数の期待値Ｅｉは、光計測中の解析窓に対応する時間領域に於ける光検出領域内の粒子の有無によって決定される。かくして、光検出領域内に粒子が存在していない場合、光子検出の事象は、常にランダムに生ずるので、各計測単位時間ｔｉに於ける期待値Ｅｎｉは、
　　　Ｅｎｉ＝Ｂｇ　　　…（３）
と設定されてよい。ここで、Ｂｇは、観測対象粒子が発光粒子である場合は、ノイズ信号の時間平均値であり、観測対象粒子が非発光粒子である場合は、背景光の時間平均値である。従って、式（３）の値を式（１）に代入して、各計測単位時間ｔｉに於ける単位時間発生確率Ｐｎｉが時系列に算定され、更に、式（２）により、光検出領域内に粒子が存在していない状態を仮定した場合に実際の検出光子数列が発生する確率Ｐｎ（第一の発生確率）が算出される。
　一方、光検出領域内に粒子が存在していた場合、光検出領域ＣＶの位置は移動しているので、粒子は、図３（Ｃ）に模式的に描かれている如く、光検出領域ＣＶ内を通過することとなる。かかる過程に於いて、光検出領域に於ける粒子から放出され検出される光の強度値又は光検出領域内に粒子が存在することによって低減する背景光の低減量は、図３（Ｄ）に示されている如く、粒子の位置が光検出領域の略中心から離れるほど低下する。従って、光検出領域内に粒子が存在していた場合の各単位時間ｔｉに於ける期待値Ｅｐｉも、時間を変数とする釣鐘状の関数として表される。ここで、その釣鐘状の関数がガウス関数で近似されるとすると、期待値Ｅｐｉは、
Figure JPOXMLDOC01-appb-M000003
 により与えられる。ここで、ガウス関数は、解析窓内の任意の時間ｔｃ（例えば、解析窓の中心）にてピーク強度Ｑを有すると仮定されている。また、式（４）のガウス関数の半値全幅は、移動速度ｖの光検出領域が、図３（Ｄ）に例示されている如き、光検出領域に於ける粒子から放出され検出される光の強度値又は光検出領域内に粒子が存在することによって低減する背景光の低減量の半径ｒ方向の分布の半値全幅ｄを通過する時間ｄ／ｖに等しく、この条件から、ｗは、
Figure JPOXMLDOC01-appb-M000004
 により与えられる。なお、図３（Ｄ）の半値全幅ｄは、光学系から決定可能である。
　式（４）中のピーク強度Ｑは、観測対象粒子が発光粒子の場合、解析窓内の光子数の総計が式（４）の期待値の総計に一致すると仮定すると、
Figure JPOXMLDOC01-appb-M000005
 により与えられる。また、観測対象粒子が非発光粒子の場合には、光子数の低減量の絶対値の期待値が式（４）に従うと設定して、ピーク強度Ｑは、
Figure JPOXMLDOC01-appb-M000006
 により与えられる。かくして、式（４）の値を式（１）に代入して、各計測単位時間ｔｉに於ける単位時間発生確率Ｐｐｉが時系列に算定され、更に、式（２）により、光検出領域内に粒子が存在していた状態を仮定した場合に実際の検出光子数列が発生する確率Ｐｐ（第二の発生確率）が算出される。そして、粒子が存在していた状態を仮定した場合の検出光子数列の発生確率Ｐｐが粒子が存在していなかった状態を仮定した場合の検出光子数列の発生確率Ｐｎを所定の程度上回るとき、その解析窓に於いて、粒子の信号が存在していると判定される。
　図４（Ａ）は、上記の一連の処理の例を示している。同図を参照して、図の例に於いて、図中中段の如く、時間軸（矢印）上の単位時間（長方形）の各々に順に光子数が［０００１１１００］となる光子数列が検出されているとき、背景光の平均が０．１であるとすると、上段右図の如く、粒子が存在していないと仮定した場合の各単位時間の期待値は全域に亘って０．１となる。一方、粒子が存在していると仮定した場合の各単位時間の期待値は、検出光子数を用いた式（６）と、更にその結果を用いた式（４）とから、上段左図の如き釣鐘状のプロファイルを有する値となる。かくして、単位時間毎の検出光子数と期待値を用いた式（１）により、単位時間毎に、検出光子数が発生する確率（ビンタイム毎の発生確率）が、粒子が存在していないと仮定した場合（下段右）と粒子が存在していたと仮定した場合（下段左）とのそれぞれについて算出される。そして、かかる検出光子数が発生する確率をそれぞれ式（２）に従って総乗して、検出された光子数列について、粒子が存在していないと仮定した場合の光子数列の発生確率Ｐｎと粒子が存在していたと仮定した場合光子数列の発生確率Ｐｐとが算出される。図示の例では、発生確率Ｐｐ＞発生確率Ｐｎとなるので、検出光子数列が粒子の信号として判定される。
　なお、実施の態様に於いては、解析窓内に粒子の信号が存在するか否かの判定は、発生確率Ｐｐと発生確率Ｐｎとのオッズ比ＯＲ
　　ＯＲ＝Ｐｐ（１－Ｐｎ）／（１－Ｐｐ）Ｐｎ　　　…（８）
を算出し、その大きさが所定値を越えたときに、解析窓内に於ける粒子の信号の存在が判定されてよい。
　上記の粒子の信号の検出処理過程に於いて、解析窓は、好適には、単一粒子が光検出領域内を通過するのに要する時間以上に設定される。半径ｒの光検出領域が速度ｖにて移動しているとすれば、解析窓の時間幅は、
　　　２ｒ／ｖ　　　…（９）
よりも長く設定される。また、解析窓は、好適には、時系列光強度データ上にて単位時間毎に順に設定されてよい。かかる設定によれば、発生確率Ｐｐ、発生確率Ｐｎ及び／又はオッズ比ＯＲが時系列光強度データに沿って算出されることとなる。しかしながら、その場合、演算量が多くなるので、数個の単位時間毎に解析窓が設定されるようになっていてもよい。更に、時系列光強度データ上を解析窓の時間幅にて分割して解析窓を設定してもよい。この場合、解析窓は、互いに重複せずに設定されることとなる。
　解析窓が単位時間毎に設定される場合には、一つの粒子信号が存在するとき、粒子の信号の存在の判定が連続する解析窓に於いて継続することとなる。即ち、一つの粒子の信号は、粒子の信号の存在の判定が継続する一つの区間に対応することとなる。従って、粒子の信号の存在の判定が継続する区間の数を計数することにより粒子の信号の計数が可能となる。また、上記の粒子の信号の検出処理過程に於いて、ビンタイムは、単一粒子が光検出領域内を通過するのに要する時間（式９）以下に設定される。これは、単一粒子が通過する際の信号を複数のビンタイムに亘って捉えて、単一粒子が通過する間の光強度値の時間変化のパターンが検出できるようにするためである。（もしビンタイムが式９の時間よりも長いと単一粒子が通過する間の光強度値の時間変化のパターンが捕捉できない。）
３．複数の光成分の計測を用いた単一粒子信号の検出
　上記の発生確率Ｐｐ、Ｐｎを用いた単一粒子信号の検出手法は、互いに異なる複数の光成分を別々に計測して成分毎に時系列光強度データを生成する場合に拡張することが可能である。その場合、観測対象となる単一粒子の発光特性が、成分毎の時系列光強度データに反映されるように、検出成分を選択することにより、かかる観測対象となる単一粒子の発光特性の反映された信号を選択的に検出することが可能となる。以下、単一粒子の発光特性の反映された信号を選択的に検出する例のいくつかについて説明する。
（ｉ）特定の偏光特性を有する単一粒子の信号の検出
　一定の方向に偏光した光を励起光とし、検出光として、発光粒子の蛍光のｐ波とｓ波を別々に検出する場合、特定の偏光特性を有する粒子の信号を選択的に検出することが可能となる（図４（Ｂ）上参照）。まず、或る発光粒子の蛍光異方性Ｒは、蛍光のｐ波とｓ波の光強度Ｃｐ、Ｃｓを用いて下記の如く与えられる。
　　Ｒ＝（Ｃｐ－Ｃｓ）／（Ｃｐ＋Ｃｓ）　　　…（１０）
　従って、全蛍光強度（Ｃｐ＋Ｃｓ）に対する蛍光のｐ波とｓ波の光強度の割合は、それぞれ
　　Ｃｐ／（Ｃｐ＋Ｃｓ）＝１／２＋Ｒ／２　　　…（１１ａ）
　　Ｃｓ／（Ｃｐ＋Ｃｓ）＝１／２－Ｒ／２　　　…（１１ｂ）
により与えられる。ところで、図４（Ｂ）下段に模式的に描かれている如く、蛍光のｓ波とｐ波の光強度の和について、粒子が存在している場合の或る解析窓内に於ける光強度の各単位時間の期待値Ｅｐｉは、式（４）と同様であるから、ピーク強度Ｑは、
Figure JPOXMLDOC01-appb-M000007
 により与えられる。ここで、Ｃｐｉ、Ｃｓｉ、Ｂｇｐ、Ｂｇｓは、それぞれ、単位時間毎のｐ波とｓ波の光子数、ｓ波とｐ波の背景光の強度である。また、蛍光のｐ波とｓ波のそれぞれの各単位時間の期待値Ｅｐｐｉ、Ｅｐｓｉは、
　　Ｅｐｐｉ＝（１／２＋Ｒ／２）Ｅｐｉ＋Ｂｇｐ　　　…（１３ａ）
　　Ｅｐｓｉ＝（１／２－Ｒ／２）Ｅｐｉ＋Ｂｇｓ　　　…（１３ｂ）
により与えられる。ここで、Ｅｐｉは、式（４）に於いて、Ｑを式（１２）により与え、ＢｇをＢｇｐ＋Ｂｇｓにより与えて得られた値である。かくして、粒子が存在していると仮定した場合の蛍光のｐ波とｓ波に於ける単位時間毎の検出光子数の発生確率Ｐｐｐｉ、Ｐｐｓｉは、式（１）と同様に
Figure JPOXMLDOC01-appb-M000008
 により与えられる。従って、蛍光異方性Ｒを有する発光粒子が存在していたと仮定した場合の解析窓内の検出光子数列の発生確率Ｐｐは、発生確率Ｐｐｐｉ、Ｐｐｓｉの総乗となるので、
Figure JPOXMLDOC01-appb-M000009
 により与えられる。一方、粒子が存在していないと仮定した場合、解析窓内の検出光子数列の発生確率Ｐｎは、
Figure JPOXMLDOC01-appb-M000010
 となる。ここで、Ｐｎｐｉ、Ｐｎｓｉは、各成分の単位時間毎の発生確率であり、それぞれ、背景光Ｂｇｐ、Ｂｇｓを期待値として用いて、式（１）と同様に
Figure JPOXMLDOC01-appb-M000011
 となる。
　かくして、発光粒子の蛍光のｐ波とｓ波を別々に検出する場合、ｐ波及びｓ波の時系列光強度データのそれぞれに於いて、式（１４）及び（１７）による単位時間毎の発生確率を算出した後、式（１５）、（１６）による解析窓に於ける発生確率Ｐｐ、Ｐｎを算出して、これらの発生確率を比較することにより、蛍光異方性Ｒを有する粒子の信号の有無が検出できることとなる。なお、この場合、蛍光異方性Ｒは、既知数である。従って、蛍光異方性Ｒの値は、任意の手法にて実験的又は理論的に決定した値が用いられてよい。理論的に蛍光異方性Ｒの値を算出する場合、例えば、非特許文献５に記載されている如く、粒子の分子量Ｍを用いて、下記の如く推定することができる。
　　Ｒ＝Ｒｏ／（１＋τθ）　　　…（１８ａ）
ここで、Ｒｏは、粒子が運動していないときの異方性の値（＝０．３）であり、τは、回転相関時間である。θは、下記の式にて与えられる。
　　θ＝ηＭ（Ｖ＋ｈ）／ＲoＴ　　　…（１８ｂ）
ここで、η、Ｖ＋ｈ、Ｒo、Ｔは、それぞれ、粘度、水和体積、気体定数、温度である。
　また、上記の如く、特定の蛍光異方性を有する単一粒子の信号を検出する態様によれば、試料溶液内に複数の互いに蛍光異方性の異なる複数種類の発光粒子が混在している場合に、種類を識別しつつ粒子の信号の検出が可能となる。即ち、粒子の種類の識別を行う場合には、各解析窓の光子数列について、上記の如き粒子が存在していると仮定した場合の発生確率が、蛍光異方性の値毎に算出される。そして、互いに異なる蛍光異方性の値を用いて得られた発生確率のうち、最も高い発生確率を与える蛍光異方性を有する粒子が、解析されている光子数列上に存在していると推定できることとなる。
（ii）特定の発光波長特性を有する単一粒子の信号の検出
　光検出領域からの光を検出する際に、互いに異なる波長帯域の光成分を別々に計測する場合、各波長帯域の成分の光強度値の大きさは、発光粒子の発光波長特性に依存して変化することとなる。例えば、図４（Ｃ）に示されている如く、互いに異なる発光波長スペクトルを有する発光粒子Ｄ１、Ｄ２について、図中、点線枠にて示された波長帯域の光成分をそれぞれ別々にＣｈ１、Ｃｈ２として検出する場合、Ｃｈ１及びＣｈ２の各々に於いて検出される各発光粒子の光は、各発光粒子の発光波長スペクトルのうち、Ｃｈ１及びＣｈ２の検出波長帯域と重複する部分（陰影のつけられた部分）となり、それぞれの光量は、陰影のつけられた部分の面積となる。従って、図から理解される如く、Ｃｈ１及びＣｈ２に於ける光量の比は、発光粒子の発光波長スペクトルのプロファイルに依存して変化することとなる。
　上記の如き発光粒子の発光波長スペクトルに依存した互いに異なる検出波長帯域にて検出された光強度（光子数）の比は、走査分子計数法による時系列光強度データ上の光子数に反映される。そこで、上記の光子数列の発生確率の算出に於いて、かかる光子数の比を所定の特性値として考慮することにより、特定の発光波長特性を有する発光粒子の信号を選択的に検出することが可能となる。
　具体的には、まず、或る発光粒子について、Ｃｈ１及びＣｈ２の検出波長帯域にて検出される光の成分の強度比α、βは、それぞれ、
　　α＝Ｓα／（Ｓα＋Ｓβ）；β＝Ｓβ／（Ｓα＋Ｓβ）…（１９）
により与えられる。Ｓα、Ｓβは、それぞれ、Ｃｈ１及びＣｈ２の検出波長帯域に於ける発光粒子の発光波長スペクトル強度の積分値（面積）である。
　Ｃｈ１及びＣｈ２の検出光子数の総和について、粒子が存在している場合の或る解析窓内に於ける光強度の各単位時間の期待値Ｅｐｉは、式（４）と同様であり、ピーク強度Ｑは、
Figure JPOXMLDOC01-appb-M000012
 により与えられる。ここで、Ｃαｉ、Ｃβｉ、Ｂｇαｐ、Ｂｇβは、それぞれ、単位時間毎のＣｈ１とＣｈ２の光子数、Ｃｈ１とＣｈ２の背景光の強度である。従って、Ｃｈ１とＣｈ２の光子数のそれぞれの期待値Ｅｐα、Ｅｐβは、
Figure JPOXMLDOC01-appb-M000013
 により与えられ、粒子が存在していると仮定した場合のＣｈ１及びＣｈ２に於ける単位時間毎の検出光子数の発生確率Ｐｐαｉ、Ｐｐβｉは、式（１）と同様に
Figure JPOXMLDOC01-appb-M000014
 により与えられる。従って、Ｃｈ１とＣｈ２の光強度比がα：βとなる発光波長特性を有する発光粒子が存在していたと仮定した場合の解析窓内の検出光子数列の発生確率Ｐｐは、発生確率Ｐｐαｉ、Ｐｐβｉの総乗となるので、
Figure JPOXMLDOC01-appb-M000015
 により与えられる。一方、粒子が存在していないと仮定した場合、解析窓内の検出光子数列の発生確率Ｐｎは、
Figure JPOXMLDOC01-appb-M000016
 となる。ここで、Ｐｎαｉ、Ｐｎβｉは、各成分の単位時間毎の発生確率であり、それぞれ、背景光Ｂｇα、Ｂｇβを期待値として用いて、式（１）と同様に
Figure JPOXMLDOC01-appb-M000017
 となる。
　かくして、発光粒子の互いに異なる発光波長帯域の成分を別々に検出する場合、各検出波長帯域の時系列光強度データのそれぞれに於いて、式（２２）及び（２５）による単位時間毎の発生確率を算出した後、式（２３）、（２４）による解析窓に於ける発生確率Ｐｐ、Ｐｎを算出して比較することにより、特定の発光波長特性を有する粒子の信号の有無が検出できることとなる。なお、この場合、Ｃｈ１とＣｈ２の光強度比α：βは、既知数である。かかる強度比は、任意の手法にて実験的又は理論的に決定した値が用いられてよい。また、上記の手法によれば、特定の発光波長特性を有する粒子の信号を選択的に検出できることとなるので、試料溶液内に複数の互いに発光波長特性の異なる複数種類の発光粒子が混在している場合に、種類を識別しつつ粒子の信号の検出が可能となる。即ち、粒子の種類の識別を行う場合には、各解析窓の光子数列について、上記の如き粒子が存在していると仮定した場合の発生確率が、Ｃｈ１とＣｈ２の光強度比（α：β）毎に算出される。そして、互いに異なる光強度比を用いて得られた発生確率のうち、最も高い発生確率を与える光強度比を有する粒子が、解析されている光子数列上に存在していると推定できることとなる。
走査分子計数法の処理操作過程
　図１（Ａ）に例示の光分析装置１を用いた本発明に従った走査分子計数法の実施形態に於いては、具体的には、（１）単一粒子を含む試料溶液の調製、（２）試料溶液の光強度の測定処理、（３）単一粒子信号の検出処理が実行される。図５は、フローチャートの形式にて表した本実施形態に於ける処理を示している。
（１）試料溶液の調製
　本発明の光分析技術の観測対象物となる粒子は、溶解された分子等の、試料溶液中にて分散し溶液中にてランダムに運動する粒子であれば、任意のものであってよく、例えば、タンパク質、ペプチド、核酸、脂質、糖鎖、アミノ酸若しくはこれらの凝集体などの生体分子、ウイルス、細胞、或いは、金属コロイド、その他の非生物学的分子などであってよい。観測対象物となる粒子が発光粒子であり、その粒子が、元来、光を発する粒子でない場合には、発光標識（蛍光分子、りん光分子、化学・生物発光分子）が観測対象物となる粒子に任意の態様にて付加されたものが用いられる。試料溶液は、典型的には水溶液であるが、これに限定されず、有機溶媒その他の任意の液体であってよい。
　反転型走査分子計数法を行う場合、観測対象物となる非発光粒子であるが、発光粒子の場合と同様に任意のものであってよい。なお、粒径が光検出領域の径の好適には１５％以上、より好適には、３５％以上であることが見出されている。また、背景光を与える発光物質としては、任意の発光分子、例えば、蛍光分子、りん光分子、化学・生物発光分子であってよく、発光物質は、光検出領域内に常に数分子以上存在する濃度にて試料溶液中に溶解又は分散される。試料溶液は、典型的には水溶液であるが、これに限定されず、有機溶媒その他の任意の液体であってもよい。
（２）試料溶液の光強度の測定（図５－ステップ１００）
　本実施形態の走査分子計数法又は反転型走査分子計数法による光分析に於ける光強度の測定は、測定中にミラー偏向器１７又はステージ位置変更装置１７ａを駆動して、試料溶液内での光検出領域の位置の移動（試料溶液内の走査）を行う他は、ＦＣＳ又はＦＩＤＡに於ける光強度の測定過程と同様の態様にて実行されてよい。操作処理に於いて、典型的には、マイクロプレート９のウェル１０に試料溶液を注入して顕微鏡のステージ上に載置した後、使用者がコンピュータ１８に対して、測定の開始の指示を入力すると、コンピュータ１８は、記憶装置（図示せず）に記憶されたプログラム（試料溶液内に於いて光検出領域の位置を移動する手順と、光検出領域の位置の移動中に光検出領域からの光を検出して時系列の光強度データを生成する手順）に従って、試料溶液内の光検出領域に於ける励起光の照射及び光強度の計測が開始される。かかる計測中、コンピュータ１８のプログラムに従った処理動作の制御下、ミラー偏向器１７又はステージ位置変更装置１７ａは、ミラー７（ガルバノミラー）又は顕微鏡のステージ上のマイクロプレート９を駆動して、ウェル１０内に於いて光検出領域の位置の移動を実行し、これと同時に光検出器１６は、逐次的に検出された光を電気信号に変換してコンピュータ１８へ送信し、コンピュータ１８では、任意の態様にて、送信された信号から時系列の光強度データを生成して保存する。なお、典型的には、光検出器１６は、一光子の到来を検出できる超高感度光検出器であるので、光の検出が、フォトンカウンティングによる場合、時系列光強度データは、時系列のフォトンカウントデータであってよい。また、互いに異なる光成分を別々に検出する場合には、複数の光検出器１６の各々が、それらの光成分の光強度値（光子数）をそれぞれ同時に検出し検出された互いに異なる成分毎に時系列光強度データが生成される。
　ところで、光検出領域の位置の移動速度は、好適には、単一粒子のランダムな運動、即ち、ブラウン運動による移動速度よりも速い値に設定される。光検出領域の位置の移動速度が粒子のブラウン運動による移動に比して遅い場合には、粒子が光検出領域内をランダムに移動することとなる。そうすると、粒子が存在している状態を仮定して算定する発生確率のための期待値の演算が複雑となり、また、精度も低下し得る。そこで、好適には、図３（Ｃ）に描かれている如く、粒子が光検出領域を略直線に横切り、光強度の変化のプロファイルの期待値が、図３（Ｄ）に描かれている如き粒子の位置が光検出領域の略中心から離れるほど低下するプロファイルと同様の釣鐘状となるように、光検出領域の位置の移動速度は、粒子のブラウン運動による平均の移動速度（拡散移動速度）よりも速く設定される。
　具体的には、拡散係数Ｄを有する単一粒子がブラウン運動によって半径ｒの光検出領域（コンフォーカルボリューム）を通過するときに要する時間Δτは、平均二乗変位の関係式
　　（２ｒ）２＝６Ｄ・Δτ　　　…（２６）
から、
　　Δτ＝（２ｒ）２／６Ｄ　　　…（２７）
となるので、単一粒子がブラウン運動により移動する速度（拡散移動速度）Ｖdifは、概ね、
　　Ｖdif＝２ｒ／Δτ＝３Ｄ／ｒ　　　…（２８）
となる。そこで、光検出領域の位置の移動速度は、かかるＶdifを参照して、それよりも十分に早い値に設定されてよい。例えば、単一粒子の拡散係数が、Ｄ＝２．０×１０－１０ｍ２／ｓ程度であると予想される場合には、ｒが、０．６２μｍ程度だとすると、Ｖdifは、１．０×１０－３ｍ／ｓとなるので、光検出領域の位置の移動速度は、その１０倍以上の１５ｍｍ／ｓと設定されてよい。なお、単一粒子の拡散係数が未知の場合には、光検出領域の位置の移動速度を種々設定して光強度の変化のプロファイルが、予想されるプロファイル（典型的には、励起光強度分布と略同様）となる条件を見つけるための予備実験を繰り返し実行して、好適な光検出領域の位置の移動速度が決定されてよい。
（３）単一粒子の信号の個別検出（ステップ１１０～１６０）
　時系列光強度データが生成されると、まず、時系列光強度データ上にて背景光の強度値の算出が行われる（ステップ１１０）。背景光の強度値は、時系列光強度データ上に於いて、粒子の信号が存在していない領域の強度値（光子数）の平均値であってよい。かくして、背景光の強度値の算出の一つの手法に於いては、通常の走査分子計数法の場合（観測対象粒子が発光粒子である場合）には、得られた時系列光強度データ上の光強度値のうち、高い方から所定の割合（例えば、２０％）のデータを除いた全強度値の平均値が背景光の強度値として採用されてよい。これは、光強度値の高い方から所定の割合のデータが粒子の信号であると考えられるためである。また、反転型走査分子計数法の場合（観測対象粒子が非発光粒子である場合）には、得られた時系列光強度データ上の光強度値のうち、低い方から所定の割合（例えば、２０％）のデータを除いた全強度値の平均値が背景光の強度値として採用されてよい。これは、光強度値の低い方から所定の割合のデータが粒子の信号であると考えられるためである。時系列光強度データが複数の成分のそれぞれに生成されているときは、成分毎に背景光の強度値の算出が行われる。なお、背景光の強度値は、観測対象の粒子を含まない試料溶液を用いて得られた時系列光強度データ上の光強度値の平均であってもよい。
　次に、本実施形態の処理過程に於いては、時系列光強度データ上に於いて解析窓の設定が実行される（ステップ１２０）。既に触れた如く、一つの解析窓の長さは、光検出領域の大きさと移動速度とから決定されてよい（式（９）参照）。また、解析窓は、好適には、ビンタイム毎に時系列に設定されてよい。しかしながら、演算量を低減する目的で、数個のビンタイム毎に設定されてもよく、或いは、解析窓が互いに重複しないように設定されてもよい。
　かくして、解析窓の設定が為されると、上記に説明された原理に従って、光検出領域内に粒子が無いと仮定した場合の解析窓内の光強度値列又は光子数列の発生確率Ｐｎ（ステップ１３０）、光検出領域内に粒子が有ると仮定した場合の解析窓内の光強度値列又は光子数列の発生確率Ｐｐ（ステップ１４０）がそれぞれ算出され、発生確率Ｐｎと発生確率Ｐｐとが比較され、解析窓内に粒子が存在しているか否かが判定される（ステップ１５０）。かかる判定に於いては、発生確率Ｐｎと発生確率Ｐｐのオッズ比が算出され（式（８）参照）、オッズ比が所定値以上のとき、粒子が存在していると判定されてよい。互いに異なる複数の成分について時系列光強度データが得られ、粒子の所定の特性値を考慮して発生確率Ｐｐが得られている場合、所定の特性値を有する粒子の信号の有無が選択的に検出されることが理解されるべきである。
　なお、光検出領域内に粒子が有ると仮定した場合の発生確率Ｐｐの算出に於いて、強度値の期待値は、解析窓の中心に強度のピークが存在するものと仮定されてよい。実際には、粒子の信号のピークは、解析窓の中心に存在しないことが殆どであるが、実際の粒子の信号のピークの位置が解析窓の中心から離れるほど、発生確率Ｐｐの値が低減するが、実際の粒子の信号が存在していないときの発生確率Ｐｐの値よりは、高い値となる。
　上記のステップ１３０～１５０の処理に於ける解析窓内の発生確率Ｐｐ、Ｐｎの算出及び粒子信号の検出は、光強度データに設定された全ての解析窓に於いて実行されてよい（ステップ１６０）。
（５）粒子濃度その他の分析
　かくして、上記の処理により、各解析窓に於ける粒子信号の有無が判定されると、光強度データ上に於ける粒子信号の計数や濃度の算出等の分析が行われてよい（ステップ１７０）。既に触れた如く、各解析窓に於ける粒子信号の有無が判定されると、一つの粒子信号が存在する解析窓が連続することとなる。従って、連続する粒子信号が存在する解析窓の組を計数することにより、粒子の数が得られることとなる。
　また、粒子の数の決定が為されている場合、任意の手法にて、光検出領域の通過した領域の総体積が算定されれば、その体積値と単一粒子の数とから試料溶液中の単一粒子の数密度又は濃度が決定される。光検出領域の通過した領域の総体積は、励起光又は検出光の波長、レンズの開口数、光学系の調整状態に基づいて理論的に算定されてよい。また、実験的に、例えば、粒子濃度が既知の溶液（対照溶液）について、検査されるべき試料溶液の測定と同様の条件にて、上記に説明した光強度の測定、単一粒子の検出及び計数を行うことにより検出された粒子の数と、対照溶液の粒子の濃度とから決定されるようになっていてよい。具体的には、例えば、粒子濃度Ｃの対照溶液について、対照溶液の粒子の検出数がＮであったとすると、光検出領域の通過した領域の総体積Ｖｔは、
　　　Ｖｔ＝Ｎ／Ｃ　　　…（２９）
により与えられる。また、対照溶液として、粒子の複数の異なる濃度の溶液が準備され、それぞれについて測定が実行されて、算出されたＶｔの平均値が光検出領域の通過した領域の総体積Ｖｔとして採用されるようになっていてよい。そして、Ｖｔが与えられると、粒子の計数結果がｎの試料溶液の発光粒子の濃度ｃは、
　　　ｃ＝ｎ／Ｖｔ　　　…（３０）
により与えられる。なお、光検出領域の体積、光検出領域の通過した領域の総体積は、上記の方法によらず、任意の方法にて、例えば、ＦＣＳ、ＦＩＤＡを利用するなどして与えられるようになっていてよい。また、本実施形態の光分析装置に於いては、想定される光検出領域の移動パターンについて、種々の標準的な粒子についての濃度Ｃと粒子の数Ｎとの関係（式（２９））の情報をコンピュータ１８の記憶装置に予め記憶しておき、装置の使用者が光分析を実施する際に適宜記憶された関係の情報を利用できるようになっていてよい。なお、互いに異なる複数の成分について時系列光強度データが得られ、粒子の所定の特性値を考慮して発生確率Ｐｐが得られている場合、所定の特性値を有する粒子の信号の計数と濃度の算出が可能であることが理解されるべきである。従って、試料溶液中に互いに異なる特性値を有する複数の種類の粒子が混在している場合、粒子の種類を識別しながら、種類毎に粒子の計数と濃度の算出が可能である
　上記に説明した本発明の有効性を検証するために、以下の如き実験を行った。なお、以下の実施例は、本発明の有効性を例示するものであって、本発明の範囲を限定するものではないことは理解されるべきである。
　走査分子計数法に於いて、本発明による時系列光強度データ上の光子数列の発生確率Ｐｐ、Ｐｎに基づいて粒子信号の検出を行った。
　試料溶液として、トリス緩衝液（１０ｍＭ　Tris-HCl(pH8.0)、０．０５％　Tween20を含む）に発光粒子としてＡＴＴＯ６４７Ｎ（シグマ）を種々の濃度にて含む溶液を調製した。光の測定に於いては、光分析装置として、共焦点蛍光顕微鏡の光学系とフォトンカウンティングシステムを備えた１分子蛍光測定装置ＭＦ２０（オリンパス株式会社）を用い、上記の「（２）試料溶液の光強度の測定」にて説明した要領に従って、上記の試料溶液について、時系列光強度データ（フォトンカウントデータ）を取得した。
　光計測後のデータ処理に於いては、まず、取得された時系列フォトンカウントデータに於いて、「（３）単一粒子の信号の個別検出」及び図５のステップ１１０～１６０に記載された要領に従って、背景光強度の算出、解析窓の設定、光検出領域内に粒子が無いと仮定した場合の解析窓内の光子数列の発生確率Ｐｎ及び光検出領域内に粒子が有ると仮定した場合の解析窓内の光子数列の発生確率Ｐｐの算出、発生確率のオッズ比の算出を順に実行し、単一粒子の信号を検出し計数した。なお、本実施例に於いては、時系列光強度データは、一つの検出波長帯域の成分のみであるので、発生確率Ｐｎ、発生確率Ｐｐは、「２．単一粒子の信号の検出」の欄に記載された式（１）、（２）、（６）等を用いて算出した。
　図６は、時系列光強度データの一部（Ａ）と、（Ａ）を平滑化して得られたデータ（Ｂ）と、（Ａ）から本発明の教示に従って算出された発生確率Ｐｎ、Ｐｐのオッズ比（Ｃ）の一つの例を示している。同図のデータの取得に於いては、励起光は、６３３ｎｍのレーザー光を用い、バンドパスフィルターを用いて、６６０－７１０ｎｍの波長帯域の光を測定し、時系列フォトンカウントデータを生成した。試料溶液中に於ける光検出領域の位置の移動速度は、３０００ｒｐｍ（７．５ｍｍ／秒）とし、ＢＩＮ　ＴＩＭＥを１０μ秒とした。発生確率Ｐｎ、Ｐｐの算出に於いては、背景光の強度は、図示の時間領域の周辺１ｍ秒の全データの光子数のうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、８０μ秒とした。また、解析窓の時間幅は、３００μ秒とした。試料溶液中のＡＴＴＯ６４７Ｎの濃度は、１ｐＭとした。
　図６を参照して、平滑化して得られたデータ（Ｂ）に於いて、特許文献９～１１に記載されている如き、平滑化データに於ける信号強度の大きさが所定の閾値による粒子信号の判定処理では、パルス状の信号の強度が所定の閾値を越えたか否かにより、粒子の信号か否かが判定される。所定の閾値は、典型的には、光子数１に設定されるので、図の略中央の（ａ）にて示した光子数の増大の除き、全てノイズとして判定される。（ａ）にて示したパルス状信号のピーク値は、光子数～１であるので、粒子信号であると思料されるところ、判定結果は誤差の影響で変わってしまい、ノイズ信号として判定される場合もあり得る（即ち、結果が不安定となる。）。
　一方、発生確率Ｐｎ、Ｐｐのオッズ比（Ｃ）に於いては、（ａ）にて示した信号に対応する領域に於いて、値の顕著な増大が見られ、従って、安定的に粒子信号であると判定できることとなる。
　図７は、図６と同様の時系列光強度データの別の一部を示した図である。同図のデータの取得に於いては、試料溶液中に於ける光検出領域の位置の移動速度は、６０００ｒｐｍ（１５ｍｍ／秒）とした。図７を参照すると、平滑化して得られたデータ（Ｂ）に於いて、（ａ）にて示したパルス状信号の他に有意なパルス状の信号が観察されるが、発生確率Ｐｎ、Ｐｐのオッズ比（Ｃ）に於いて（ａ）にて示したパルス状信号に対応する領域以外には、有意な増大は見られなかった。このことは、本発明の教示に従って算出された発生確率Ｐｎ、Ｐｐに基づく判定によれば、粒子信号とノイズ信号との識別が明確となり、粒子判定の安定性が向上されることを示唆している。
　図８は、種々の濃度の発光粒子を含む試料溶液について本発明による発生確率Ｐｎ、Ｐｐのオッズ比を用いて粒子信号の判定を行った場合の発光粒子濃度と検出粒子数との関係を示している。なお、本実験に際しては、励起光は、６４２ｎｍのレーザー光（３ｍＷ）を用い、バンドパスフィルターを用いて、５６０－６２０ｎｍの波長帯域の光を測定し、時系列フォトンカウントデータを生成した。試料溶液中に於ける光検出領域の位置の移動速度は、１２０００ｒｐｍ（９０ｍｍ／秒）とし、ＢＩＮ　ＴＩＭＥを１０μ秒とし、光の測定時間を６００秒間とした。発生確率Ｐｎ、Ｐｐの算出に於いては、背景光の強度は、各解析窓の周辺１ｍ秒のデータのうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、５０μ秒とした。また、解析窓の時間幅は、３５０μ秒とした。そして、オッズ比＞１０８を越える領域に粒子信号が存在するものと判定した。図８から理解される如く、１ｆＭ以上の濃度域に於いて、検出粒子数は、発光粒子濃度と共に増大した。このことは、本発明の粒子信号の検出手法により、時系列フォトンカウントデータ上の粒子の信号を検出できること、また、試料溶液中の粒子濃度を１ｆＭ以上の濃度域に於いて決定できることを示している。
　反転型走査分子計数法に於いて、本発明による時系列光強度データ上の光子数列の発生確率Ｐｐ、Ｐｎに基づいて粒子信号の検出を行った。
　試料溶液として、背景光を生成する発光物質として、１ｍＭ　ＡＴＴＯ６４７Ｎを含むトリス緩衝液（１０ｍＭ　Tris-HCl(pH8.0)、０．０５％　Tween20を含む）にポリスチレンビーズ（直径４μｍ）を分散した溶液を調製した。光の測定に於いては、光分析装置として、共焦点蛍光顕微鏡の光学系とフォトンカウンティングシステムを備えた１分子蛍光測定装置ＭＦ２０（オリンパス株式会社）を用い、上記の「（２）試料溶液の光強度の測定」にて説明した要領に従って、上記の試料溶液について、時系列光強度データ（フォトンカウントデータ）を取得した。
　光計測後のデータ処理に於いては、まず、取得された時系列フォトンカウントデータに於いて、「（３）単一粒子の信号の個別検出」及び図５のステップ１１０～１６０に記載された要領に従って、背景光強度の算出、解析窓の設定、光検出領域内に粒子が無いと仮定した場合の解析窓内の光子数列の発生確率Ｐｎ及び光検出領域内に粒子が有ると仮定した場合の解析窓内の光子数列の発生確率Ｐｐの算出、発生確率のオッズ比の算出を順に実行し、単一粒子の信号を計数した。なお、本実施例に於いては、時系列光強度データは、一つの検出波長帯域の成分のみであるので、発生確率Ｐｎ、発生確率Ｐｐは、「２．単一粒子の信号の検出」の欄に記載された式（１）、（２）、（７）等を用いて算出した。
　図９（Ａ）、（Ｂ）、図１０は、得られた時系列光強度データの一部を順に拡大した図（上段）とそのオッズ比を示している。同図のデータの取得に於いては、励起光は、６３３ｎｍのレーザー光（５０μＷ）を用い、バンドパスフィルターを用いて、６６０－７１０ｎｍの波長帯域の光を測定し、時系列フォトンカウントデータを生成した。試料溶液中に於ける光検出領域の位置の移動速度は、６０００ｒｐｍ（１５ｍｍ／秒）とし、ＢＩＮ　ＴＩＭＥを５０μ秒とした。発生確率Ｐｎ、Ｐｐの算出に於いては、背景光の強度は、各解析窓の周辺１ｍ秒のデータのうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、１００μ秒とした。また、解析窓の時間幅は、３００μ秒とした。ビーズの濃度は、１００ｆＭとした。
　これらの図を参照して、オッズ比の増大は、フォトンカウントデータ上の背景光の低減部分と対応していることが理解される。これにより、本発明の粒子信号検出の手法は、反転型走査分子計数法にも適用可能であることが示唆された。
　図１１（Ａ）は、種々の濃度のビーズを含む試料溶液について本発明による発生確率Ｐｎ、Ｐｐのオッズ比を用いて粒子信号の判定を行った場合の粒子濃度と検出粒子数との関係を示している。なお、本実験に際しては、励起光は、６３３ｎｍのレーザー光（５０μＷ）を用い、バンドパスフィルターを用いて、６６０－７１０ｎｍの波長帯域の光を測定し、時系列フォトンカウントデータを生成した。試料溶液中に於ける光検出領域の位置の移動速度は、９０００ｒｐｍ（９０ｍｍ／秒）とし、ＢＩＮ　ＴＩＭＥを１０μ秒とし、光の測定時間を２００秒間とした。発生確率Ｐｎ、Ｐｐの算出に於いては、背景光の強度は、各解析窓の周辺１ｍ秒のデータのうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、１００μ秒とした。また、解析窓の時間幅は、３００μ秒とした。そして、オッズ比＞１０１０を越える領域に粒子信号が存在するものと判定した。なお、比較のため、図１１（Ｂ）は、同一のデータについて、特許文献９～１１に記載の如く時系列フォトンカウントデータを平滑化した後、ガウス関数をフィッティングしてピーク光子数≧５であるパルス状信号を検出して得られた検出粒子数と粒子濃度との関係を示している。
　図１１（Ａ）、（Ｂ）を参照して、図示の濃度域（１ａＭ～１ｆＭ）に於いて、平滑化された時系列フォトンカウントデータに於いてガウス関数をフィッティングして得られた結果（Ｂ）の場合、濃度の増大に対して検出粒子数の増大が観察されるが、ばらつきが比較的大きかった（近似直線の相関係数ｒ２＝０．９３）。これに対し、本発明による粒子信号の検出手法による結果（Ａ）に於いては、検出粒子数は、粒子濃度に概ね比例した（近似直線の相関係数ｒ２＝０．９９）。また、同図から、試料溶液中の粒子濃度を５０ａＭ以上の濃度域に於いて決定できることを示している。（（Ｂ）の場合は、決定可能な濃度範囲は、５００ａＭ以上である。）
　蛍光異方性の異なる二種類の発光粒子について、本発明による走査分子計数法により種類毎に粒子の検出が可能であることを検証した。
　試料溶液として、リン酸緩衝液（０．０５％　Tween20を含む）に蛍光色素TAMRA(Ｍ．Ｗ．４３０．４５　Sigma-Aldrich　Cat.No.C2734)を１００ｆＭにて含む溶液（蛍光色素溶液）と、リン酸緩衝液に１ｐＭのプラスミド（pbr322　２．９ＭＤａ　タカラバイオ株式会社　Cat.No.3035）と１０ｎＭのＤＮＡインターカレータ蛍光色素SYTOX　Orange（インビトロジェン社　Cat.No.S-11368）を含む溶液（プラスミド溶液、SYTOX　Orangeが単一のプラスミドに結合して単一の発光粒子となる。）を調製した。光の測定に於いては、光分析装置として、共焦点蛍光顕微鏡の光学系とフォトンカウンティングシステムを備えた１分子蛍光測定装置ＭＦ２０（オリンパス株式会社）を用い、上記の「（２）試料溶液の光強度の測定」にて説明した要領に従って、上記の試料溶液について、ｓ偏光成分及びｐ偏光成分の時系列光強度データ（フォトンカウントデータ）を同時に且つ別々に取得した。その際、ＴＡＭＲＡ溶液及びプラスミド溶液の両方について、励起光は、５４３ｎｍのレーザー光を用い、バンドパスフィルターを用いて、５６０－６２０ｎｍの波長帯域の光を検出した。また、励起光の偏光方向は、検出光のｐ偏光成分と同じ方向となるよう設定した。試料溶液中に於ける光検出領域の移動速度は、６０００rpm（１５ｍｍ／秒）と設定し、ＢＩＮ　ＴＩＭＥを１０μ秒とし、測定時間は、２秒間とした。上記の試料溶液について、時系列光強度データ（フォトンカウントデータ）を取得した。
　光計測後のデータ処理に於いては、まず、取得された時系列フォトンカウントデータに於いて、「（３）単一粒子の信号の個別検出」及び図５のステップ１１０～１６０に記載された要領に従って、背景光強度の算出、解析窓の設定、光検出領域内に粒子が無いと仮定した場合の解析窓内の光子数列の発生確率Ｐｎ及び光検出領域内に粒子が有ると仮定した場合の解析窓内の光子数列の発生確率Ｐｐの算出、発生確率のオッズ比の算出を順に実行し、単一粒子の信号を計数した。なお、本実施例に於いては、ｐ波とｓ波の二成分の時系列光強度データから、蛍光異方性を既知のパラメータとして含む発生確率Ｐｎ、発生確率Ｐｐを、「３．複数の光成分の計測を用いた単一粒子信号の検出（ｉ）特定の偏光特性を有する単一粒子の信号の検出」の欄に記載された式（１５）、（１６）等を用いて算出した。その際、背景光の強度は、各解析窓の周辺１ｍ秒のデータのうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、１２０μ秒とした。また、解析窓の時間幅は、４００μ秒とした。また、発生確率Ｐｐに関しては、プラスミドの蛍光異方性の０．４を用いた値と、ＴＡＭＲＡの蛍光異方性の０．３２を用いた値とを算出し（即ち、一つの解析窓について、二つの発生確率Ｐｐが算出される。）、それぞれについて、発生確率Ｐｎとのオッズ比が２０を越えたとき、その解析窓に対応する種類の粒子の信号が存在するものと判定した。（プラスミドの発生確率Ｐｐのオッズ比とＴＡＭＲＡの発生確率Ｐｐのオッズ比とが共に２０を越えたときには、オッズ比の高い種類の粒子信号が存在するものとした。）
　図１２は、プラスミド溶液と蛍光色素溶液について上記の走査分子計数法により得られた時系列光強度データに於いて、発生確率Ｐｎと互いに異なる蛍光異方性を用いて得られた発生確率Ｐｐとのオッズ比を用いて、粒子信号の存在を判定し、粒子の種類毎に粒子数を計数した結果を示している。同図を参照して、図から理解される如く、それぞれの溶液に於いて、大半の検出された粒子信号がそれぞれの種類に対応して検出された。（プラスミド溶液中では、４２１個の粒子信号が検出され、そのうち、ＴＡＭＲＡとして誤検出された数は、８２個であった。一方、蛍光色素溶液中では、２７９個の粒子信号が検出され、そのうち、プラスミドとして誤検出された数は、５８個であった。）この結果は、本発明の手法により、種類を識別しつつ試料溶液中の単一粒子の検出が可能であることを示唆している。
　かくして、上記の実施例の結果から理解される如く、本発明の教示に従って、走査分子計数法に於いて時系列光強度データに於ける光強度値列の発生確率を、粒子が存在したと仮定した場合と粒子が存在していなかったと仮定した場合とのそれぞれについて、算出し、算出された発生確率に基づいて、粒子の信号の有無を判定することにより、粒子の信号の検出精度又はＳ／Ｎ比が向上され、また、検出感度（許容可能な精度にて粒子の信号が検出可能な試料溶液の粒子濃度の範囲）が改善できることが示された。

Claims (30)

  1. 　共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出する光分析装置であって、
    　前記試料溶液内に於ける前記光学系の光検出領域の位置を移動する光検出領域移動部と、
    　前記光検出領域からの光を検出する光検出部と、
    　前記試料溶液内に於いて前記光検出領域の位置を移動させながら前記光検出部にて検出された前記光検出領域からの光の時系列の光強度データを生成し、前記時系列の光強度データに於いて前記単一粒子の各々の存在を表す信号を個別に検出する信号処理部とを含み、
    　前記信号処理部が、前記時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、前記光検出領域内に前記単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、前記光検出領域内に前記単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出し、前記第一及び第二の発生確率に基づいて前記時系列の光強度データ上にて前記単一粒子の各々の存在を表す信号を検出することを特徴とする装置。
  2. 　請求項１の装置であって、前記第一及び第二の発生確率が、それぞれ、前記解析窓内に於ける単位時間毎の光強度値と前記第一及び第二の状態を仮定した場合の前記単位時間毎の期待値とに基づいて算出されることを特徴とする装置。
  3. 　請求項２の装置であって、前記単位時間毎の光強度値が前記単位時間毎の前記期待値を有するポアソン分布に従うものとして前記単位時間毎の光強度値の単位時間発生確率が算出され、前記第一及び第二の発生確率がそれぞれ対応する前記単位時間発生確率を用いて算出されることを特徴とする装置。
  4. 　請求項１乃至３の装置であって、前記第二の発生確率が前記第一の発生確率より大きい解析窓の時間に前記光検出領域内に前記単一粒子が存在したと判定されることを特徴とする装置。
  5. 　請求項１の装置であって、前記単一粒子が所定の特性値を有し、前記光検出部が前記光検出領域からの光の少なくとも二つの互いに異なる成分を別々に検出し、前記信号処理部が、前記成分の各々の時系列の光強度データを生成し、更に、前記信号処理部が、前記成分の各々の前記第一及び第二の発生確率を算出し、前記成分の各々の前記第二の発生確率が前記所定の特性値の関数であり、前記成分毎の前記第一及び第二の発生確率の発生確率に基づいて前記時系列の光強度データ上にて前記所定の特性値を有する前記単一粒子の存在を表す信号が検出されることを特徴とする装置。
  6. 　請求項５の装置であって、前記単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含み、前記単一粒子の種類毎に、前記互いに異なる所定の特性値の関数である前記成分の各々の前記第二の発生確率が算出され、前記成分の各々の前記第一の発生確率と前記複数の種類の単一粒子の各々の前記成分の各々の前記第二の発生確率とに基づいて前記時系列の光強度データ上にて前記単一粒子の種類毎に該単一粒子の存在を表す信号が検出されることを特徴とする装置。
  7. 　請求項１乃至６の装置であって、前記単一粒子が発光粒子であり、前記単一粒子の各々の存在を表す信号が光強度の一時的な増大であることを特徴とする装置。
  8. 　請求項５又は６を引用する請求項７の装置であって、前記所定の特性値が前記単一粒子の偏光異方性であることを特徴とする装置。
  9. 　請求項５又は６を引用する請求項７の装置であって、前記所定の特性値が前記単一粒子の互いに異なる発光波長帯域の発光強度の比であることを特徴とする装置。
  10. 　請求項１乃至４の装置であって、前記単一粒子が検出波長帯域にて発光しない粒子であり、前記光検出領域からの光が背景光を含み、前記単一粒子の各々の存在を表す信号が前記背景光からの光強度の一時的な低減であることを特徴とする装置。
  11. 　共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出する光分析方法であって、
    　前記試料溶液内に於いて前記光学系の光検出領域の位置を移動する過程と、
    　前記試料溶液内に於いて前記光検出領域の位置を移動させながら前記光検出領域からの光の強度を測定して光強度データを生成する過程と、
    　前記時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、前記光検出領域内に前記単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、前記光検出領域内に前記単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出する過程と、
    　前記第一及び第二の発生確率に基づいて前記時系列の光強度データ上にて前記単一粒子の各々の存在を表す信号を検出する過程と
    を含むことを特徴とする方法。
  12. 　請求項１１の方法であって、前記第一及び第二の発生確率が、それぞれ、前記解析窓内に於ける単位時間毎の光強度値と前記第一及び第二の状態を仮定した場合の前記単位時間毎の期待値とに基づいて算出されることを特徴とする方法。
  13. 　請求項１２の方法であって、前記単位時間毎の光強度値が前記単位時間毎の前記期待値を有するポアソン分布に従うものとして前記単位時間毎の光強度値の単位時間発生確率が算出され、前記第一及び第二の発生確率がそれぞれ対応する前記単位時間発生確率を用いて算出されることを特徴とする方法。
  14. 　請求項１１乃至１３請求項の方法であって、前記第二の発生確率が前記第一の発生確率より大きい解析窓の時間に前記光検出領域内に前記単一粒子が存在したと判定されることを特徴とする方法。
  15. 　請求項１１の方法であって、前記単一粒子が所定の特性値を有し、前記光検出領域からの光の少なくとも二つの互いに異なる成分が別々に検出されて、前記成分の各々の時系列の光強度データが生成され、更に、前記成分の各々の前記第一及び第二の発生確率が算出され、前記成分の各々の前記第二の発生確率が前記所定の特性値の関数であり、前記成分毎の前記第一及び第二の発生確率の発生確率に基づいて前記時系列の光強度データ上にて前記所定の特性値を有する前記単一粒子の存在を表す信号が検出されることを特徴とする方法。
  16. 　請求項１５の方法であって、前記単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含み、前記単一粒子の種類毎に、前記互いに異なる所定の特性値の関数である前記成分の各々の前記第二の発生確率が算出され、前記成分の各々の前記第一の発生確率と前記複数の種類の単一粒子の各々の前記成分の各々の前記第二の発生確率とに基づいて前記時系列の光強度データ上にて前記単一粒子の種類毎に該単一粒子の存在を表す信号が検出されることを特徴とする方法。
  17. 　請求項１１乃至１６の方法であって、前記単一粒子が発光粒子であり、前記単一粒子の各々の存在を表す信号が光強度の一時的な増大であることを特徴とする方法。
  18. 　請求項１５又は１６を引用する請求項１７の方法であって、前記所定の特性値が前記単一粒子の偏光異方性であることを特徴とする方法。
  19. 　請求項１５又は１６を引用する請求項１７の方法であって、前記所定の特性値が前記単一粒子の互いに異なる発光波長帯域の発光強度の比であることを特徴とする方法。
  20. 　請求項１１乃至１４の方法であって、前記単一粒子が検出波長帯域にて発光しない粒子であり、前記光検出領域からの光が背景光を含み、前記単一粒子の各々の存在を表す信号が前記背景光からの光強度の一時的な低減であることを特徴とする方法。
  21. 　共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出するための光分析用コンピュータプログラムであって、
    　前記試料溶液内に於いて前記光学系の光検出領域の位置を移動する手順と、
    　前記試料溶液内に於いて前記光検出領域の位置を移動させながら前記光検出領域からの光の強度を測定して光強度データを生成する手順と、
    　前記時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度の時間変化値の、前記光検出領域内に前記単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、前記光検出領域内に前記単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出する手順と、
    　前記第一及び第二の発生確率に基づいて前記時系列の光強度データ上にて前記単一粒子の各々の存在を表す信号を検出する手順と
    をコンピュータに実行させることを特徴とするコンピュータプログラム。
  22. 　請求項２１のコンピュータプログラムであって、前記第一及び第二の発生確率が、それぞれ、前記解析窓内に於ける単位時間毎の光強度値と前記第一及び第二の状態を仮定した場合の前記単位時間毎の期待値とに基づいて算出されることを特徴とするコンピュータプログラム。
  23. 　請求項２２のコンピュータプログラムであって、前記単位時間毎の光強度値が前記単位時間毎の前記期待値を有するポアソン分布に従うものとして前記単位時間毎の光強度値の単位時間発生確率が算出され、前記第一及び第二の発生確率がそれぞれ対応する前記単位時間発生確率を用いて算出されることを特徴とするコンピュータプログラム。
  24. 　請求項２１乃至２３請求項のコンピュータプログラムであって、前記第二の発生確率が前記第一の発生確率より大きい解析窓の時間に前記光検出領域内に前記単一粒子が存在したと判定されることを特徴とするコンピュータプログラム。
  25. 　請求項２１のコンピュータプログラムであって、前記単一粒子が所定の特性値を有し、前記光検出領域からの光の少なくとも二つの互いに異なる成分が別々に検出されて、前記成分の各々の時系列の光強度データが生成され、更に、前記成分の各々の前記第一及び第二の発生確率が算出され、前記成分の各々の前記第二の発生確率が前記所定の特性値の関数であり、前記成分毎の前記第一及び第二の発生確率の発生確率に基づいて前記時系列の光強度データ上にて前記所定の特性値を有する前記単一粒子の存在を表す信号が検出されることを特徴とするコンピュータプログラム。
  26. 　請求項２５のコンピュータプログラムであって、前記単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含み、前記単一粒子の種類毎に、前記互いに異なる所定の特性値の関数である前記成分の各々の前記第二の発生確率が算出され、前記成分の各々の前記第一の発生確率と前記複数の種類の単一粒子の各々の前記成分の各々の前記第二の発生確率とに基づいて前記時系列の光強度データ上にて前記単一粒子の種類毎に該単一粒子の存在を表す信号が検出されることを特徴とするコンピュータプログラム。
  27. 　請求項２１乃至２６のコンピュータプログラムであって、前記単一粒子が発光粒子であり、前記単一粒子の各々の存在を表す信号が光強度の一時的な増大であることを特徴とするコンピュータプログラム。
  28. 　請求項２５又は２６を引用する請求項２７のコンピュータプログラムであって、前記所定の特性値が前記単一粒子の偏光異方性であることを特徴とするコンピュータプログラム。
  29. 　請求項２５又は２６を引用する請求項２７のコンピュータプログラムであって、前記所定の特性値が前記単一粒子の互いに異なる発光波長帯域の発光強度の比であることを特徴とするコンピュータプログラム。
  30. 　請求項２１乃至２４のコンピュータプログラムであって、前記単一粒子が検出波長帯域にて発光しない粒子であり、前記光検出領域からの光が背景光を含み、前記単一粒子の各々の存在を表す信号が前記背景光からの光強度の一時的な低減であることを特徴とするコンピュータプログラム。
PCT/JP2013/052446 2012-02-17 2013-02-04 単一粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム WO2013121905A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201380009765.5A CN104115001B (zh) 2012-02-17 2013-02-04 利用单个粒子检测的光分析装置、光分析方法
EP13748583.5A EP2816344A4 (en) 2012-02-17 2013-02-04 OPTICAL ANALYSIS DEVICE USING SINGLE PARTICLE DETECTION TECHNIQUE, OPTICAL ANALYSIS METHOD, AND COMPUTER PROGRAM FOR OPTICAL ANALYSIS
JP2014500168A JP5996625B2 (ja) 2012-02-17 2013-02-04 単一粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
US14/451,021 US9494779B2 (en) 2012-02-17 2014-08-04 Optical analysis device, optical analysis method and computer program for optical analysis using single particle detection

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2012-032421 2012-02-17
JP2012032421 2012-02-17

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US14/451,021 Continuation US9494779B2 (en) 2012-02-17 2014-08-04 Optical analysis device, optical analysis method and computer program for optical analysis using single particle detection

Publications (1)

Publication Number Publication Date
WO2013121905A1 true WO2013121905A1 (ja) 2013-08-22

Family

ID=48984019

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2013/052446 WO2013121905A1 (ja) 2012-02-17 2013-02-04 単一粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム

Country Status (5)

Country Link
US (1) US9494779B2 (ja)
EP (1) EP2816344A4 (ja)
JP (1) JP5996625B2 (ja)
CN (1) CN104115001B (ja)
WO (1) WO2013121905A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015052965A1 (ja) * 2013-10-07 2015-04-16 オリンパス株式会社 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JPWO2013157319A1 (ja) * 2012-04-18 2015-12-21 オリンパス株式会社 光分析を用いた単一粒子検出装置、単一粒子検出方法及び単一粒子検出用コンピュータプログラム
US9494779B2 (en) 2012-02-17 2016-11-15 Olympus Corporation Optical analysis device, optical analysis method and computer program for optical analysis using single particle detection
JP2019523887A (ja) * 2016-06-20 2019-08-29 プレアー ソシエテ・アノニム 流体中の浮遊粒子を検出および／または特性評価するための装置と方法

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150226573A1 (en) * 2014-02-11 2015-08-13 Qualcomm Incorporated Pedometer integrated pedestrian positioning
JP6475069B2 (ja) * 2015-04-23 2019-02-27 アズビル株式会社 粒子検出装置及び粒子の検出方法
CA3018582A1 (en) * 2016-03-25 2017-09-28 Bioceryx Inc. Apparatuses and methods for assessing target sequence numbers
DE102018208647A1 (de) * 2018-05-30 2019-12-05 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Lasermesseinrichtung zur Messung einer Distanz zu einem Objekt sowie Verfahren zum Betreiben derselben
CN112683873B (zh) * 2021-03-22 2021-06-01 泛肽生物科技(浙江)有限公司 液态芯片的检测模型构建方法及装置、分析方法及装置
CN113720744B (zh) * 2021-11-04 2022-01-25 碧兴物联科技(深圳)股份有限公司 一种基于偏振检测技术的大气颗粒物含量实时监测方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04337446A (ja) 1991-05-15 1992-11-25 Hitachi Ltd 微粒子計測方法、定量方法および微粒子計測装置
JPH10512952A (ja) * 1995-01-16 1998-12-08 ソイニ，エルッキ フロー蛍光法及び装置
JP2005098876A (ja) 2003-09-25 2005-04-14 Institute Of Physical & Chemical Research ２成分相互作用分析方法
JP2007020565A (ja) 2005-06-13 2007-02-01 Kansai Bunri Sogo Gakuen 試料中のウイルスを検出する方法およびシステム
JP4023523B2 (ja) 1996-10-12 2007-12-19 オリンパス株式会社 比明度関数の決定によるサンプルの分析方法
JP2008116440A (ja) 2006-10-13 2008-05-22 Shiga Pref Gov 試料中の蛍光性物質を検出する方法およびシステム
WO2008080417A1 (en) 2006-12-28 2008-07-10 Flult Biosystems Gmbh A method of determining characteristic properties of a sample containing particles
JP2008292371A (ja) 2007-05-25 2008-12-04 Institute Of Physical & Chemical Research 蛍光相関分光による複数成分相互作用の分析方法及び相互作用制御化合物のスクリーニング方法
JP2009281831A (ja) 2008-05-21 2009-12-03 Hamamatsu Photonics Kk 蛍光解析装置及び解析方法
JP2010017127A (ja) * 2008-07-10 2010-01-28 Olympus Corp 標的核酸比率推定方法
JP2011508219A (ja) * 2007-12-19 2011-03-10 シンギュレックス・インコーポレイテッド 単一分子検出走査分析器およびその使用方法
WO2011108369A1 (ja) 2010-03-01 2011-09-09 オリンパス株式会社 光分析装置、光分析方法並びに光分析用コンピュータプログラム

Family Cites Families (48)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4251733A (en) 1978-06-29 1981-02-17 Hirleman Jr Edwin D Technique for simultaneous particle size and velocity measurement
US5866336A (en) 1996-07-16 1999-02-02 Oncor, Inc. Nucleic acid amplification oligonucleotides with molecular energy transfer labels and methods based thereon
US6235471B1 (en) 1997-04-04 2001-05-22 Caliper Technologies Corp. Closed-loop biochemical analyzers
US6710871B1 (en) 1997-06-09 2004-03-23 Guava Technologies, Inc. Method and apparatus for detecting microparticles in fluid samples
GB2326229A (en) 1997-06-13 1998-12-16 Robert Jeffrey Geddes Carr Detecting and analysing submicron particles
SE9800360D0 (sv) 1998-02-06 1998-02-06 Goeteborg University Science I Method, apparatus and flow cell for high sensitivity detection of fluorescent molecules
US20030036855A1 (en) 1998-03-16 2003-02-20 Praelux Incorporated, A Corporation Of New Jersey Method and apparatus for screening chemical compounds
KR100560588B1 (ko) 1998-03-16 2006-03-16 에머샘 바이오사이언시즈 코프 공초점 마이크로스코피 영상 시스템
US6388788B1 (en) 1998-03-16 2002-05-14 Praelux, Inc. Method and apparatus for screening chemical compounds
US6376843B1 (en) 1999-06-23 2002-04-23 Evotec Oai Ag Method of characterizing fluorescent molecules or other particles using generating functions
DK1175602T3 (da) 1999-04-29 2003-05-26 Evotec Ag Fremgangsmåde til karakterisering af fluorescente molekyler eller andre partikler, der anvender frembringerfunktioner
US8264680B2 (en) 1999-05-28 2012-09-11 Yokogawa Electric Corporation Biochip reader and electrophoresis system
US6965113B2 (en) 2000-02-10 2005-11-15 Evotec Ag Fluorescence intensity multiple distributions analysis: concurrent determination of diffusion times and molecular brightness
US20010035954A1 (en) 2000-03-10 2001-11-01 Rahn John Richard Method and apparatus for measuring particle size distributions using light scattering
DE10035190C5 (de) 2000-07-20 2009-07-16 MAX-PLANCK-Gesellschaft zur Förderung der Wissenschaften e.V. Verfahren und Vorrichtung zur Fluoreszenzmessung
DE10038528A1 (de) 2000-08-08 2002-02-21 Zeiss Carl Jena Gmbh Verfahren und Anordnung zur Erhöhung der spektralen und räumlichen Detektorauflösung
US6947133B2 (en) 2000-08-08 2005-09-20 Carl Zeiss Jena Gmbh Method for increasing the spectral and spatial resolution of detectors
CN1353305A (zh) * 2000-11-10 2002-06-12 北京理工大学 超高灵敏度荧光显微探测与处理系统
HU226937B1 (en) 2000-11-17 2010-03-29 Mta Szegedi Biolog Koezpont Method and apparatus for determining polarization amount of material by a laser scanning microscope
JP3999662B2 (ja) 2000-12-14 2007-10-31 オリンパス株式会社 蛍光分析装置および蛍光分析方法
US6782297B2 (en) 2001-05-24 2004-08-24 Eric Paul Tabor Methods and apparatus for data smoothing
US6750963B2 (en) 2002-05-21 2004-06-15 Agilent Technologies, Inc. Imaging systems for signals on a surface
AU2003283860B2 (en) 2002-11-07 2009-09-17 Erasmus Universiteit Rotterdam Fret probes and methods for detecting interacting molecules
US7038848B2 (en) 2002-12-27 2006-05-02 Olympus Corporation Confocal microscope
JP4315794B2 (ja) 2002-12-27 2009-08-19 オリンパス株式会社 共焦点顕微鏡
AU2005290314A1 (en) 2004-09-28 2006-04-06 Singulex, Inc. System and method for spectroscopic analysis of single particles
CN101147055A (zh) 2005-01-31 2008-03-19 伊利诺伊大学评议会 用于表征清澈和混浊介质中的颗粒的方法和设备
WO2007010803A1 (ja) 2005-07-15 2007-01-25 Olympus Corporation 光測定装置
WO2007118209A2 (en) 2006-04-07 2007-10-18 Kim Laboratories Apparatus and method for rapid detection of analytes
US8445655B2 (en) 2006-06-16 2013-05-21 Cornell Research Foundation, Inc. Functional nucleic acid ligands to fluorescent proteins
WO2008007580A1 (fr) 2006-07-13 2008-01-17 Olympus Corporation Procédé d'analyse de particules fines
US20080052009A1 (en) 2006-08-23 2008-02-28 Washington, University Of Method for deconvolving single-molecule intensity distributions for quantitative biological measurements
GB0618057D0 (en) 2006-09-14 2006-10-25 Perkinelmer Ltd Improvements in and relating to scanning confocal microscopy
JPWO2008099778A1 (ja) 2007-02-14 2010-05-27 株式会社ニコン スリット走査共焦点顕微鏡
JP2009145242A (ja) 2007-12-14 2009-07-02 Olympus Corp 光測定装置
US7619734B2 (en) 2008-03-03 2009-11-17 Agilent Technologies, Inc. Methods and systems for computing a size distribution of small particles
JP2009288161A (ja) 2008-05-30 2009-12-10 Olympus Corp 光測定装置及び光測定方法
US20100177190A1 (en) 2008-12-16 2010-07-15 Ann-Shyn Chiang Microscopy system with revolvable stage
JP2011002415A (ja) 2009-06-22 2011-01-06 Olympus Corp 蛍光相関分光装置
JP5856962B2 (ja) 2010-09-10 2016-02-10 オリンパス株式会社 単一発光粒子の光強度を用いた光分析方法
EP2615445B1 (en) 2010-10-13 2014-05-21 Olympus Corporation Method of measuring a diffusion characteristic value of a particle
EP2667183A4 (en) 2011-01-20 2017-05-10 Olympus Corporation Photoanalysis method and photoanalysis device using detection of light from single light-emitting particle
JP5885738B2 (ja) * 2011-04-13 2016-03-15 オリンパス株式会社 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JP2013036765A (ja) 2011-08-04 2013-02-21 Olympus Corp 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
WO2013031309A1 (ja) 2011-08-26 2013-03-07 オリンパス株式会社 光分析を用いた単一粒子検出装置、単一粒子検出方法及び単一粒子検出用コンピュータプログラム
WO2013069504A1 (ja) 2011-11-10 2013-05-16 オリンパス株式会社 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JP5996625B2 (ja) 2012-02-17 2016-09-21 オリンパス株式会社 単一粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JP6010114B2 (ja) * 2012-04-18 2016-10-19 オリンパス株式会社 光分析を用いた単一粒子検出装置、単一粒子検出方法及び単一粒子検出用コンピュータプログラム

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04337446A (ja) 1991-05-15 1992-11-25 Hitachi Ltd 微粒子計測方法、定量方法および微粒子計測装置
JPH10512952A (ja) * 1995-01-16 1998-12-08 ソイニ，エルッキ フロー蛍光法及び装置
JP4023523B2 (ja) 1996-10-12 2007-12-19 オリンパス株式会社 比明度関数の決定によるサンプルの分析方法
JP2005098876A (ja) 2003-09-25 2005-04-14 Institute Of Physical & Chemical Research ２成分相互作用分析方法
JP2007020565A (ja) 2005-06-13 2007-02-01 Kansai Bunri Sogo Gakuen 試料中のウイルスを検出する方法およびシステム
JP2008116440A (ja) 2006-10-13 2008-05-22 Shiga Pref Gov 試料中の蛍光性物質を検出する方法およびシステム
WO2008080417A1 (en) 2006-12-28 2008-07-10 Flult Biosystems Gmbh A method of determining characteristic properties of a sample containing particles
JP2008292371A (ja) 2007-05-25 2008-12-04 Institute Of Physical & Chemical Research 蛍光相関分光による複数成分相互作用の分析方法及び相互作用制御化合物のスクリーニング方法
JP2011508219A (ja) * 2007-12-19 2011-03-10 シンギュレックス・インコーポレイテッド 単一分子検出走査分析器およびその使用方法
JP2009281831A (ja) 2008-05-21 2009-12-03 Hamamatsu Photonics Kk 蛍光解析装置及び解析方法
JP2010017127A (ja) * 2008-07-10 2010-01-28 Olympus Corp 標的核酸比率推定方法
WO2011108369A1 (ja) 2010-03-01 2011-09-09 オリンパス株式会社 光分析装置、光分析方法並びに光分析用コンピュータプログラム
WO2011108371A1 (ja) 2010-03-01 2011-09-09 オリンパス株式会社 光分析装置、光分析方法並びに光分析用コンピュータプログラム
WO2011108370A1 (ja) 2010-03-01 2011-09-09 オリンパス株式会社 光分析装置、光分析方法並びに光分析用コンピュータプログラム

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
ANAL. CHEM., vol. 70, no. 3, 1998, pages 632
MASATAKA KINJO, PROTEIN, NUCLEIC ACID, ENZYME, vol. 44, no. 9, 1999, pages 1431 - 1438
NORIKO KATO ET AL., GENE MEDICINE, vol. 6, no. 2, pages 271 - 277
P. KASK; K. PALO; D. ULLMANN; K. GALL, PNAS, vol. 96, 1999, pages 13756 - 13761
R. RIGLER: "Fluorescence Correlation Spectroscopy", 2000, SPRINGER, article F.J.MEYER-ALMS, pages: 204 - 224
See also references of EP2816344A4 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9494779B2 (en) 2012-02-17 2016-11-15 Olympus Corporation Optical analysis device, optical analysis method and computer program for optical analysis using single particle detection
JPWO2013157319A1 (ja) * 2012-04-18 2015-12-21 オリンパス株式会社 光分析を用いた単一粒子検出装置、単一粒子検出方法及び単一粒子検出用コンピュータプログラム
WO2015052965A1 (ja) * 2013-10-07 2015-04-16 オリンパス株式会社 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
CN105593667A (zh) * 2013-10-07 2016-05-18 奥林巴斯株式会社 利用单个发光粒子检测的光分析装置、光分析方法以及光分析用计算机程序
JPWO2015052965A1 (ja) * 2013-10-07 2017-03-09 オリンパス株式会社 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
US9739698B2 (en) 2013-10-07 2017-08-22 Olympus Corporation Optical analysis device, optical analysis method and computer program for optical analysis using single light-emitting particle detection
JP2019523887A (ja) * 2016-06-20 2019-08-29 プレアー ソシエテ・アノニム 流体中の浮遊粒子を検出および／または特性評価するための装置と方法

Also Published As

Publication number Publication date
JPWO2013121905A1 (ja) 2015-05-11
EP2816344A1 (en) 2014-12-24
CN104115001B (zh) 2016-08-24
CN104115001A (zh) 2014-10-22
US9494779B2 (en) 2016-11-15
US20140339444A1 (en) 2014-11-20
JP5996625B2 (ja) 2016-09-21
EP2816344A4 (en) 2015-09-23

Similar Documents

Publication Publication Date Title
JP5996625B2 (ja) 単一粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JP5802653B2 (ja) 光分析装置、光分析方法並びに光分析用コンピュータプログラム
JP5885738B2 (ja) 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JP5904996B2 (ja) 単一発光粒子検出を用いた光分析装置、光分析方法並びに光分析用コンピュータプログラム
WO2013031309A1 (ja) 光分析を用いた単一粒子検出装置、単一粒子検出方法及び単一粒子検出用コンピュータプログラム
JP6010114B2 (ja) 光分析を用いた単一粒子検出装置、単一粒子検出方法及び単一粒子検出用コンピュータプログラム
JP6013338B2 (ja) 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JP5893564B2 (ja) 複数の波長帯域の光計測を用いた光分析方法
JP5941923B2 (ja) 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JP6013337B2 (ja) 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JP5907882B2 (ja) 単一発光粒子の偏光特性を観測する光分析装置、光分析方法及びそのための光分析用コンピュータプログラム
JP6010035B2 (ja) 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
WO2013118519A1 (ja) 光分析を用いた単一粒子検出装置、単一粒子検出方法及び単一粒子検出用コンピュータプログラム
JP5945506B2 (ja) 単一発光粒子の光の波長特性を用いた光分析装置及び光分析方法
JP2013036765A (ja) 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
WO2013021984A1 (ja) 共焦点顕微鏡又は多光子顕微鏡の光学系を用いた光分析装置及び光分析方法
JP2014202550A (ja) 複数の波長帯域に於ける単一発光粒子検出を用いた光分析装置及び光分析方法
JP6157326B2 (ja) 光検出を用いた単一発光粒子検出方法
JP2013019764A (ja) 共焦点顕微鏡又は多光子顕微鏡の光学系を用いた光分析装置及び光分析方法
JP2014149195A (ja) 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
JP2014149196A (ja) 単一発光粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム
WO2016125278A1 (ja) 単一発光粒子検出を用いた光分析装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13748583

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2014500168

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2013748583

Country of ref document: EP