WO2020235283A1 - 分析装置および分析方法 - Google Patents

分析装置および分析方法 Download PDF

Info

Publication number
WO2020235283A1
WO2020235283A1 PCT/JP2020/017202 JP2020017202W WO2020235283A1 WO 2020235283 A1 WO2020235283 A1 WO 2020235283A1 JP 2020017202 W JP2020017202 W JP 2020017202W WO 2020235283 A1 WO2020235283 A1 WO 2020235283A1
Authority
WO
WIPO (PCT)
Prior art keywords
points
emission
light emitting
fluorescence
time
Prior art date
Application number
PCT/JP2020/017202
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 SG11202112784UA priority Critical patent/SG11202112784UA/en
Priority to CN202080036799.3A priority patent/CN113840902A/zh
Priority to US17/612,046 priority patent/US20220236183A1/en
Priority to GB2116497.5A priority patent/GB2597202B/en
Priority to JP2021520664A priority patent/JP7282880B2/ja
Priority to DE112020002033.9T priority patent/DE112020002033T5/de
Publication of WO2020235283A1 publication Critical patent/WO2020235283A1/ja

Links

Images

Classifications

    • 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/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/59—Transmissivity
    • 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
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
    • C12M1/00—Apparatus for enzymology or microbiology
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
    • C12M1/00—Apparatus for enzymology or microbiology
    • C12M1/34—Measuring or testing with condition measuring or sensing means, e.g. colony counters
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869—Methods for sequencing
    • 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
    • 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
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/26—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating electrochemical variables; by using electrolysis or electrophoresis
    • G01N27/416—Systems
    • G01N27/447—Systems using electrophoresis
    • 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
    • G01N2021/6417—Spectrofluorimetric devices
    • G01N2021/6421—Measuring at two or more wavelengths
    • 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/6428—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • G01N2021/6439—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" with indicators, stains, dyes, tags, labels, marks
    • G01N2021/6441—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" with indicators, stains, dyes, tags, labels, marks with two or more labels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/26—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating electrochemical variables; by using electrolysis or electrophoresis
    • G01N27/416—Systems
    • G01N27/447—Systems using electrophoresis
    • G01N27/44704—Details; Accessories
    • G01N27/44717—Arrangements for investigating the separated zones, e.g. localising zones
    • G01N27/44721—Arrangements for investigating the separated zones, e.g. localising zones by optical means

Definitions

  • the present disclosure relates to an analysis method and an analyzer that detect fluorescence emitted from a plurality of types of phosphors at a plurality of emission points while identifying each of them.
  • the measurement mode is a mode in which the lamp light is irradiated to the absorption point of each capillary to detect the absorption of the lamp light generated when the analysis target passes through the absorption point, or the laser light is applied to the emission point of each capillary.
  • modes such as irradiating and detecting fluorescence or scattered light generated when the analysis target passes through the light emitting point.
  • Patent Document 1 all the capillaries around the A light emitting points on the capillaries of A (A is an integer of 2 or more) are arranged on the same plane, and the laser beam is introduced from the side of the arrangement plane.
  • the emission points of all the capillaries are irradiated at once, and the fluorescence generated at each emission point is wavelength-dispersed from the direction perpendicular to the array plane for batch detection.
  • the fluorescence emitted from the A light emitting points is collectively collimated by one condenser lens, transmitted through one transmission type diffraction grating, and the primary diffracted light of each fluorescence is one.
  • the imaging lens collectively forms an image on one two-dimensional sensor.
  • B color detection is possible by setting B (B is an integer of 1 or more) detection regions of any wavelength band for the wavelength dispersion image of each capillary.
  • B 1 is called single color detection
  • B ⁇ 2 is called multicolor detection.
  • DNA sequencing of different DNA samples by the Sanger method can be performed in each capillary.
  • DNA fragments contained in DNA samples are labeled with four types of phosphors according to the terminal base species A, C, G, and T, and the emission fluorescence of each is identified by multicolor detection. ..
  • Patent Document 2 all the capillaries around the A light emitting points on the capillaries of A (A is an integer of 2 or more) are arranged on the same plane, and the laser beam is introduced from the side of the arrangement plane to all.
  • the emission points of the capillary are collectively irradiated, and the fluorescence generated at each emission point is divided according to the wavelength component from the direction perpendicular to the arrangement plane and collectively detected.
  • the emission fluorescence from A emission points is individually collimated with A condenser lenses to obtain A luminous flux, and B (B is an integer of 1 or more) dichroic mirrors are arranged1.
  • Each luminous flux is incident in parallel on a set of dichroic mirror arrays, each is divided into B luminous fluxes in different wavelength bands, and the generated total A ⁇ B luminous flux is incident in parallel on one two-dimensional sensor. , A ⁇ B divided images are generated on the image.
  • the A ⁇ B divided images overlap each other on the image. It will be possible to set A ⁇ B detection areas. This makes it possible to detect the B color of each capillary. Therefore, in the multi-capillary electrophoresis apparatus of Patent Document 2, for example, as in the case of Patent Document 1, DNA sequencing of different DNA samples by the Sanger method can be performed in each capillary.
  • the emission fluorescence of the C type (C is an integer of 1 or more) phosphor emission and the wavelength band of the B type (B is an integer of 1 or more) are detected.
  • B color detection is performed at each time in the area. However, B ⁇ C.
  • the B color detection result is subjected to color conversion to obtain the concentration of the C type phosphor.
  • the largest element should be 1, and the other elements should be shown as a ratio to the maximum value.
  • the entire row of the matrix Y can be determined by performing the above steps in order for all the C type phosphors D (c).
  • the matrix Y is determined only by the characteristics of the fluorophore D (c) and the detection regions W (b) of different wavelength bands and does not change during the electrophoretic analysis.
  • the matrix Y is kept constant for different electrophoretic analyzes. Therefore, for each emission point, the concentration Z (c) of the phosphor D (c) at each time can be obtained from the signal intensity X (b) of the detection region W (b) at each time by the following equation.
  • Equation 1 is a simultaneous equation showing the relationship between the concentration of the unknown C type phosphor and the known B color fluorescence intensity, and (Equation 6) corresponds to finding the solution. Therefore, in general, the condition B ⁇ C is required as described above. If B ⁇ C, the solution cannot be uniquely obtained (that is, there can be multiple solutions), so the color conversion cannot be performed as in (Equation 6).
  • the emission fluorescence of the phosphor D (c) (c is 1, 2, 3, or 4) is in the wavelength band due to spectral crosstalk.
  • W (b) (b is 1, 2, 3, or 4) represents the intensity ratio detected.
  • the elements Y (b) (c) of the matrix Y are obtained by electrophoretic analysis of a sample in which the phosphor D (c) (c is 1, 2, 3, or 4) fluoresces independently. Can be decided.
  • the four-color fluorescence intensities X (1), X (2), X (3), and X (4) when the phosphor D (1) fluoresces independently are the elements Y (1), respectively.
  • Y (b) (c) is a fixed value determined only by the characteristics of the phosphor D (c) and the wavelength band W (b), and does not change during electrophoresis.
  • the concentrations of the phosphors D (1), D (2), D (3), and D (4) at each time are the four-color fluorescence intensities X (1), X (2) at each time. ), X (3), and X (4), it is obtained by the following equation that embodies (Equation 6).
  • the above color conversion process is performed independently for each of the A capillaries.
  • This is based on the premise that the crosstalk between capillaries (referred to as spatial crosstalk in this disclosure) is sufficiently small. In other words, it means that the ratio of the signal intensity obtained by detecting the fluorescence emitted from any one capillary to the signal intensity derived from the fluorescence emitted from the other capillary is sufficiently small. Therefore, the step of determining the matrix Y by sequentially performing the above-mentioned steps of fluorescing one type of phosphor D (c 0 ) independently for the C type phosphor D (c) is a capillary of A capillaries. Are performed substantially at the same time, and the matrix Y of each of the A capillaries is obtained in parallel. This is an inevitable process for deriving the matrix Y in a short time and without any hassle.
  • Spatial crosstalk is basically reduced by devising the optical system, but there are also attempts to reduce it by calculation processing.
  • Patent Document 3 not a plurality of capillaries, but emission fluorescence from a plurality of light emitting points randomly arranged on a plane is collectively imaged on one two-dimensional sensor by one condenser lens. ..
  • the emission fluorescence from each emission point is obtained as the signal intensity of the detection region provided at each imaging position on the two-dimensional sensor.
  • the ratio of spatial crosstalk that the signal intensity of emission fluorescence from any one emission point gives to the signal intensity of emission fluorescence from another emission point is the distance between those two emission points, or the corresponding two detections.
  • the above function is obtained in advance, and in the fluorescence images of a plurality of emission points randomly arranged on a plane, the distance between two arbitrary detection regions and the mutual spatial crosstalk are obtained from the above function, and the original fluorescence image is used. By subtracting the obtained spatial crosstalk, we are trying to reduce the spatial crosstalk.
  • the present inventors tried to reduce spatial crosstalk by the method of Patent Document 3 for the optical system of Patent Document 1 or Patent Document 2, but they did not work well.
  • the signal intensity X ( ⁇ ) ( ⁇ ) of the detection region ⁇ for the emission point ⁇ which is any one emission point of A emission points, is any one emission point ⁇ '( ⁇ ⁇ ) other than the above.
  • the spatial crosstalk given to the signal strength X ( ⁇ ') ( ⁇ ) of the detection region ⁇ for ⁇ ') tends to decrease with the distance between the two detection regions, one of the above distances It became clear that it could not be represented by the function of.
  • the two detection regions are located near the central axis of the optical system and far from the central axis of the optical system.
  • the ratio of spatial crosstalk was different.
  • the ratio of spatial crosstalk from one of the above two detection regions to the other was different from the ratio of spatial crosstalk from the other to one.
  • the method of Patent Document 3 has stopped working.
  • the signal intensity X ( ⁇ ) ( ⁇ ) of the detection region ⁇ , which is any one detection region of the B detection regions, is the signal intensity X ( ⁇ ) ( ⁇ ) of the emission point ⁇ , which is any one emission point of A emission points.
  • any one detection region ⁇ '( ⁇ ⁇ 'or ⁇ ⁇ ⁇ 'of B detection regions.
  • ⁇ ⁇ ⁇ ' spatial crosstalk and spectral crosstalk coexist, so the crosstalk between the above two detection regions is essentially impossible to express as a function of the distance between them. there were. This subject will be described in detail in [Modes for Carrying Out the Invention].
  • the spatial crosstalk between a plurality of light emitting points which occurs in an arbitrary optical system that detects the light emitted from the plurality of light emitting points, is reduced by calculation processing, so that the light emitted from the plurality of light emitting points is emitted.
  • spatial crosstalk between multiple light emitting points which occurs in an arbitrary optical system that detects fluorescence of multiple types of phosphors emitted from multiple light emitting points, and multiple types of each light emitting point.
  • the absorption at multiple absorption points can be identified and each can be identified.
  • spatial crosstalk between multiple absorption points and between multiple types of absorbers for each absorption point occur in an arbitrary optical system that detects the absorption of multiple types of absorbers at multiple absorption points.
  • the concentration of the phosphor D (a, c) at the emission point P (a) at an arbitrary time is Z (a, c), and the signal in the detection region W (a', b) for the emission point P (a').
  • X is an A ⁇ B row and 1 column matrix whose elements are X (a', b)
  • Z is an A ⁇ C row and 1 column matrix whose elements are Z (a, c)
  • Y is Y
  • (Equation 9) is the same formula as (Equation 1), but when comparing (Equation 1) to (Equation 4) and (Equation 9) to (Equation 12), it can be seen that they are completely different.
  • (Equation 9) to (Equation 12) are not only the relational expressions of b and c, but also the relational expression of a, and different emission points P (a) are related to each other. That is, (Equation 9) to (Equation 12) collectively consider spatial crosstalk and spectral crosstalk between different emission points P (a) in addition to spectral crosstalk for the same emission point P (a). It is made possible, and is essentially different from (Equation 1) to (Equation 4).
  • A' ⁇ a that is, the emission fluorescence of phosphor D (a, c) at emission point P (a) is arbitrary detection of P (a') by spatial crosstalk and spectral crosstalk for different emission points.
  • it is generally difficult to control the concentration of the phosphor D (a 0 , c 0 ) it is convenient to standardize the row Y (a, b) (a 0 , c 0 ). ..
  • the largest element should be 1, and the other elements should be shown as a ratio to the maximum value.
  • the entire row of the matrix Y can be determined by sequentially performing the above steps for all combinations of the C-type phosphors D (a, c) at the A light emitting points P (a).
  • the matrix Y is determined only by the characteristics of the emission point P (a), the phosphor D (a, c), and the detection region W (a, b), and does not change during the electrophoretic analysis.
  • the matrix Y is constant for different electrophoretic analyzes. Is kept in. Therefore, for each emission point, the concentration Z (a, c) of the phosphor D (a, c) at each time is determined from the signal intensity X (a, b) of the detection region W (a, b) at each time. It is calculated by the following equation.
  • Equation 14 is the same equation as (Equation 6), but comparing (Equation 2) to (Equation 4) and (Equation 6) with (Equation 10) to (Equation 12) and (Equation 14), It can be seen that the two are completely different.
  • the spectral crosstalk and the spatial crosstalk are Both can be eliminated and the concentration Z (a, c) of A ⁇ C phosphors can be obtained.
  • Y - is a matrix of (A ⁇ C) rows (A ⁇ B) columns.
  • eliminating spectral crosstalk is referred to as color conversion
  • eliminating spatial crosstalk is referred to as spatial correction.
  • the (A ⁇ C) ⁇ (A ⁇ B) elements of Y ⁇ include both the part that executes color conversion and the part that executes spatial correction.
  • the time series data of Z (a, c) obtained by executing (Equation 14) for each time is referred to as color conversion + spatial correction data in the present disclosure.
  • Each element of the matrix X may be a value obtained by subtracting the background light in advance, or may be a value subjected to appropriate noise reduction processing. Similarly, for each element of the matrix Y, the value may be obtained by subtracting the background light in advance, or an appropriate change may be added.
  • Equation 9 and (Equation 14) to (Equation 17) are formally the same as the general formulas used in Patent Document 3, for example, [Equation 7] to [Equation 9] in Patent Document 3.
  • the contents are significantly different (hereinafter, the mathematical formulas used in Patent Document 3 are described by []).
  • the element ⁇ ij of the matrix A of [Equation 7] is a function of the distance dij between the detection regions of the emission points ⁇ (meas) i and ⁇ (meas) j, and is specifically shown in [Equation 10]. As you can see, it is represented by the sum of the exponential functions of dij, and ⁇ ij decays as dij increases.
  • Patent Document 3 different fluorescence images are obtained for different samples, and the positions of the detection regions of a plurality of emission points change randomly in each fluorescence image. Therefore, in Patent Document 3, the above function is obtained in advance, and for each fluorescence image, the distance between the detection regions of any two emission points among the plurality of emission points is the combination of all the two emission points. Is obtained, and ⁇ ij is derived by substituting it into the above function. That is, ⁇ ij, that is, the matrix A changes for each sample or each fluorescence image.
  • the elements Y (a') (a) of the matrix Y in (Equation 16) are the light emitting points P ( We have found that it cannot be expressed as a function of the distance between each detection area of a) and P (a'). It has also been found that it is impossible to calculate the elements Y (a') (a) from the configuration of the optical system. Therefore, in the present disclosure, the elements Y (a') and (a) are actually obtained under the condition that the positions of the detection regions of a plurality of light emitting points do not change even for different samples.
  • the concentration of the illuminant D (a, c) at the emission point P (a) is Z (a, c)
  • the emission intensity of W (a', b) for P (a') is X ( Let a', b).
  • (Equation 9) to (Equation 17) are established, and similarly, both spectral crosstalk and spatial crosstalk are eliminated, and the concentration Z (a, c) of A ⁇ C illuminants is determined.
  • the light emission includes fluorescence, phosphorescence, scattered light, and the like.
  • the concentration of the absorber D (a, c) at the absorption point P (a) at an arbitrary time is Z (a, c)
  • the absorbance of W (a', b) at the absorption point P (a') is X.
  • the concentration of the signal generator D (a, c) at the signal generation point P (a) at an arbitrary time is Z (a, c), and the concentration of W (a', b) for the signal generation point P (a').
  • the signal strength be X (a', b).
  • (Equation 9) to (Equation 17) are established, and similarly, both spectral crosstalk and spatial crosstalk are eliminated, and the concentration Z (a, c) of A ⁇ C signal generators is satisfied. Can be sought.
  • spatial crosstalk between a plurality of light emitting points which occurs in an arbitrary optical system for detecting light emission from a plurality of light emitting points, can be eliminated or reduced by calculation processing from a plurality of light emitting points. It is possible to identify the light emission of and detect each of them independently.
  • spatial crosstalk between multiple light emitting points which occurs in an arbitrary optical system that detects the fluorescence of multiple types of phosphors emitted from multiple light emitting points, and multiple types of each light emitting point.
  • Raw data from electrophoretic analysis of unknown samples Color conversion data obtained by performing color conversion on the raw data shown in FIG.
  • Raw data obtained by electrophoretic analysis of samples containing unknown phosphors Color conversion + spatial correction data obtained by performing color conversion + spatial correction on the raw data in Fig. 17.
  • Color conversion data obtained by performing color conversion on the raw data shown in FIG. Color conversion + spatial correction data obtained by performing color conversion + spatial correction on the raw data in Fig. 19.
  • Example 1 In order to investigate the characteristics of spatial crosstalk in detail, we constructed the simple optical system shown in Fig. 1.
  • the optical system in Fig. 1 is a pinhole plate 1-1, a light emitting point side opening plate 1-2, a condenser lens 1-3, a sensor side opening plate 1-4, a colored glass filter 1-5, and a two-dimensional sensor 1-. 6.
  • Halogen lamp A halogen lamp (light source) that irradiates light 1-7 is provided.
  • the optical system shown in FIG. 1 was configured as follows.
  • a light emitting point 1-8 having a diameter of 0.05 mm was formed by irradiating a pinhole plate 1-1 having a pinhole having a diameter of 0.05 mm with halogen lamp light 1-7 from below.
  • 3 was placed at a position 1.54 mm away from the light emitting point 1-8 in the upward direction, and the sensor-side aperture plate 1-4 having an aperture of ⁇ 0.7 mm was placed directly above the condenser lens 1-3.
  • the two-dimensional sensor 1-6 was placed at a position 15 mm above the condenser lens 1-3, and the colored glass filter 1-5 was placed directly under the two-dimensional sensor 1-6.
  • the pinhole plate 1-1, the light emitting point side opening plate 1-2, the sensor side opening plate 1-4, the colored glass filter 1-5, and the two-dimensional sensor 1-6 are arranged in parallel with each other.
  • the light 1-9 emitted from the light emitting points 1-8 is transmitted through the aperture of ⁇ 0.2 mm, condensed by the condenser lens 1-3, transmitted through the aperture of ⁇ 0.7 mm, and the colored glass filter 1-
  • a light emission image 1-10 with a diameter of 0.5 mm was formed on the two-dimensional sensor 1-6 through 5.
  • the light emitting points 1-8 were imaged at an image magnification of 10 times with the two-dimensional sensor 1-6 in focus.
  • Fig. 2 is a light emitting image including images 1-10 obtained by the two-dimensional sensor 1-6 in the simple optical system of Fig. 1.
  • the sensor size of the two-dimensional sensor 1-6 is 13 ⁇ 13 mm, and the signal range of each pixel is 0 to 65536.
  • Fig. 2 (a) and Fig. 2 (b) are the same emission image, but the signal display scale (gradation) of Fig. 2 (a) is set to 0 to 50000, and the signal display scale of Fig. 2 (b) is set to 0 to 5. It is set to 500.
  • the maximum signal value of the luminescent image is about 50,000, and while the signals of all pixels are displayed almost not saturated in FIG. 2 (a), the luminescent image is displayed saturated in FIG. 2 (b).
  • FIGS. 3 (a) to 3 (c) shows the signal intensity distribution in the horizontal direction passing through the center of the luminescent image in the luminescent image of FIG. Fig. 3 (a), Fig. 3 (b), and Fig. 3 (c) use the same data, but the vertical axis is changed.
  • the horizontal axis is common and indicates the distance from the center of the luminescent image.
  • the plus and minus on the horizontal axis indicate the right side and the left side of the luminescent image in the luminescent image, respectively.
  • the vertical axis of FIGS. 3 (a) and 3 (b) shows the absolute signal strength, the vertical axis scale of FIG. 3 (a) is 0 to 60000, and the vertical axis scale of FIG.
  • FIGS. 3 (a) to 3 (c) show the case where the output intensity of the halogen lamp light 1-7 is gradually reduced under the conditions for acquiring the emission image of FIG. Similar to ⁇ , it shows the signal intensity distribution.
  • the maximum signal intensities of the emission images of ⁇ , ⁇ , and ⁇ are about 50,000, about 25,000, and about 10,000, respectively.
  • the signal intensity tends to decrease as the distance from the center of the luminescent image increases, but it is compared with the size of the luminescent image (width of about 0.5 mm) seen in Fig. 3 (a). It can be seen that the hem is much larger. It can also be seen that as the maximum signal strength decreases, so does the strength of the hem. Therefore, looking at Fig. 3 (c), it can be seen that ⁇ , ⁇ , and ⁇ overlap. This is a new discovery and brings some important findings. For example, at a position ⁇ 1 mm away from the center of the luminescent image, that is, at a position outside the luminescent image seen in Fig.
  • Fig. 4 (a) and Fig. 5 (a) are the results derived from Fig. 3.
  • Fig. 4 (a) shows positions -1, -2, -3, -4, -5, -6 mm away from the center of the imaging point
  • Fig. 5 (a) shows +1 and + from the center of the imaging point.
  • the graph is a graph in which the absolute signal intensity at positions 2, +3, +4, +5, +6 mm is on the vertical axis, and the absolute signal intensity at the center of the emission image is on the horizontal axis in each case.
  • the above absolute signal strength is the average value of the absolute signal strength with a width of ⁇ 0.1 mm centered at each position in Fig. 3.
  • the absolute signal intensity of spatial crosstalk at an arbitrary position away from the center of the imaging point is linear with respect to the absolute signal intensity of the center of the imaging point, and the spatial crosstalk ratio is constant. It shows that it becomes. That is, the spatial crosstalk at an arbitrary position is obtained by subtracting the value obtained by multiplying the absolute signal intensity at the center of the imaging point by the spatial crosstalk ratio from the absolute signal intensity at an arbitrary position away from the center of the imaging point. It was newly found that it is possible to eliminate or reduce. This corresponds to subtracting the value of the corresponding approximate straight line at the same position from the absolute signal strength of each plot in FIGS. 4 (a) and 5 (a).
  • Fig. 4 (b) and Fig. 5 (b) are the results of performing the above deduction operation on the results of Fig. 4 (a) and Fig. 5 (a), respectively.
  • the types of plots showing each position and the horizontal axis are common to Fig. 4 (a) and Fig. 4 (b), and Fig. 5 (a) and Fig. 5 (b).
  • the vertical axes of FIGS. 4 (b) and 5 (b) show the absolute signal strength after the above subtraction calculation, and the scale is enlarged as compared with FIGS. 4 (a) and 5 (a). There is.
  • a straight line with a spatial crosstalk of ⁇ 0.01% is overlaid with a dotted line.
  • Fig. 6 and Fig. 7 are the same experimental results as in Fig. 4 and Fig. 5.
  • the pinhole plate 1-1 was once removed and reattached, and then the same results as in Fig. 2 and Fig. 3 were obtained, and the results in Fig. 6 and Fig. 7 were derived by the same method.
  • the mounting position of the pinhole plate 1-1 is almost the same as the original position, but it is not exactly the same position.
  • the results in FIGS. 6 and 7 are similar to the results in FIGS. 4 and 5, indicating that the method of eliminating or reducing spatial crosstalk at an arbitrary position is highly reproducible by subtraction calculation processing.
  • Patent Document 3 it was found that it is impossible to express the spatial crosstalk ratio that occurs between the center of the imaging point and an arbitrary position as a function of the distance between the two. This means that the method shown in Patent Document 3 does not work effectively.
  • the conditions of the optical system and a plurality of imaging points, that is, a plurality of detection regions are fixed, and the mutual spatial crosstalk ratio is shown in each of FIGS. 4 to 7 regardless of the mutual distance. As described above, it differs from Patent Document 3 in that it is derived experimentally. Whenever the conditions are changed, it is necessary to reacquire the spatial crosstalk ratio.
  • Example 2 A model experiment was conducted on the simplest case of an optical system that detects the fluorescence of multiple types of phosphors emitted from multiple emission points in the detection regions of multiple wavelength bands.
  • A 2 emission points
  • B 2 wavelength bands, respectively.
  • C the same effect can be obtained by the same method when each of A, B, and C has other numerical values. Needless to say.
  • Figure 8 shows a model experimental system.
  • the model experimental system in FIG. 8 is a light source that irradiates a laser beam LB in the arrangement direction of two capillarys Cap (1) and Cap (2) and two capillarys Cap (1) and Cap (2) (not shown). ) And. Emission points P (1) and P (2) are provided at positions where the laser beam LB is irradiated on the two capillary caps (1) and Cap (2), respectively.
  • Two types of phosphors D (1,1) and D (1,2) labeled substances are electrophoresed on the capillary Cap (1), and two types of phosphors D (2,1) are electrophoresed on the capillary Cap (2).
  • D (2,2) labeled substances are electrophoresed.
  • D (1,1) and D (2,1) are the same type of phosphor
  • D (1,2) and D (2,2) are the same type of phosphor.
  • the same type of phosphor means that the fluorescence spectra of the emitted fluorescence are equal.
  • the phosphors D (1,1) and D (1,2) emit fluorescence by irradiation with the laser beam LB at the emission point P (1), and the detection region W (1,1) provided in the sensor S. Detected by and W (1,2).
  • the detection region W (1,1) mainly detects the emission fluorescence of the phosphor D (1,1)
  • the detection region W (1,2) mainly detects the emission fluorescence of the phosphor D (1,2).
  • the detection wavelength band of is designed, mutual spectral crosstalk is significantly present.
  • the phosphors D (2,1) and D (2,2) emit fluorescence by irradiation with the laser beam LB at the emission point P (2), and the detection region W (2) provided in the sensor S is provided. Detected by, 1) and W (2,2).
  • the detection region W (2,1) mainly detects the emission fluorescence of the phosphor D (2,1)
  • the detection region W (2,2) mainly detects the emission fluorescence of the phosphor D (2,2).
  • the detection wavelength band of is designed, mutual spectral crosstalk is significantly present.
  • the fluorescence emission of the phosphors D (1,1) and D (1,2) at the emission point P (1) is also detected in the detection regions W (2,1) and W (2,2), and similarly.
  • Fluorescence emission of phosphors D (2,1) and D (2,2) at the emission point P (2) is also detected in the detection regions W (1,1) and W (1,2), and spatial crosstalk Is significantly present.
  • the signal intensities (fluorescence intensities) in the detection regions W (1,1), W (1,2), W (2,1), and W (2,2) are X (1,1), X (1,2). ), X (2,1), and X (2,2), as shown in the lower part of Fig. 8, electrophoresis for each of the emission point P (1) and emission point P (2).
  • time series data of signal intensities X (1,1) and X (1,2), X (2,1) and X (2,2) are obtained.
  • the detection areas W (1,1), W (1,2), W (2,1), and W (2,2) are provided on one sensor S. It is drawn, but it is not always necessary.
  • the detection areas W (1,1) and W (1,2) may be provided on one sensor, and the detection areas W (2,1) and W (2,2) may be provided on the other sensor. However, the detection areas W (1,1), W (1,2), W (2,1), and W (2,2) may be provided on four different sensors, respectively.
  • some imaging means and spectroscopic means are required between the light emitting points P (1) and P (2) and the sensor S. What kind of imaging means and spectroscopic means are used in this embodiment? In FIG. 8, the imaging means and the spectroscopic means are omitted.
  • FIG. 9 is a schematic diagram showing the relationship between spectral crosstalk and spatial crosstalk in FIG.
  • the fluorescence emission of the phosphor D (1,1) at the emission point P (1) is mainly detected in (i) the detection region W (1,1), and (ii) the detection region W (ii). It is detected as a secondary in 1 and 2).
  • the emission fluorescence was detected in (iii) detection region W (2,1) with much lower intensity than in (i) and (ii), and (iv) detection region W (2,2). However, it is detected with a lower intensity than (iii).
  • (i) and (ii) show spectral crosstalk
  • (iii) and (iv) show spatial crosstalk and spectral crosstalk.
  • FIG. 9 shows the relationship between the strength and weakness of the above detection intensity and the relationship between crosstalk.
  • FIG. 9 (b) shows the fluorescence emission of the phosphor D (1,2) at the emission point P (1)
  • FIG. 9 (c) shows the fluorescence D at the emission point P (2).
  • the fluorescence emission of (2,1) and FIG. 9 (d) are drawn for the fluorescence emission of the phosphor D (2,2) at the emission point P (2).
  • FIGS. 9 (a) to 9 (d) occur all at once at the same time, and their fluorescence intensities are different.
  • Fig. 10 shows a model experimental system similar to Fig. 8.
  • the lamp light LL is applied to the two capillary Caps (1) and Caps (2) from the light source perpendicular to the arrangement direction. Therefore, the emission fluorescence associated with the irradiation of the laser beam LB is not detected, but the absorbance or absorbance associated with the transmission of the lamp light LL is detected.
  • Absorption points P (1) and P (2) are provided on the two capillaries Cap (1) and Cap (2) at positions where the lamp light LL is irradiated.
  • D (1,1) and D (1,2) are electrophoresed on the capillary Cap (1)
  • two types of absorbers D (2,1) and D (2) are electrophoresed on the capillary Cap (2).
  • , 2) are electrophoresed.
  • D (1,1) and D (2,1) are the same type of absorber
  • D (1,2) and D (2,2) are the same type of absorber.
  • the same type of absorber means that the light absorption spectra are equal.
  • the absorbers D (1,1) and D (1,2) absorb the light at the absorption point P (1) by the irradiation of the lamp light LL, and the unabsorbed light is the detection region provided in the sensor S. Detected at W (1,1) and W (1,2).
  • the detection region W (1,1) mainly detects the absorption of the absorber D (1,1)
  • the detection region W (1,2) mainly detects the absorption of the absorber D (1,2).
  • the wavelength band is designed, there is significant mutual spectral crosstalk.
  • the absorbers D (2,1) and D (2,2) absorb the light at the absorption point P (2) by the irradiation of the lamp light LL, and the unabsorbed light is provided in the sensor S. It is detected in the detection areas W (2,1) and W (2,2).
  • the detection region W (2,1) mainly detects the absorption of the absorber D (2,1)
  • the detection region W (2,2) mainly detects the absorption of the absorber D (2,2).
  • the wavelength band is designed, there is significant mutual spectral crosstalk.
  • the absorption of the absorbers D (1,1) and D (1,2) at the absorption point P (1) was also detected in the detection regions W (2,1) and W (2,2).
  • the absorption of the absorbers D (2,1) and D (2,2) at the absorption point P (2) was also detected in the detection regions W (1,1) and W (1,2), and the spatial crosstalk was significant. Exists in. This is because a part of the transmitted light of the lamp light LL transmitted through the absorption point P (1) is also detected in the detection regions W (2,1) and W (2,2) and is transmitted through the absorption point P (2).
  • Figure 11 shows the raw data obtained by electrophoresis analysis by injecting known samples into the capillaries Cap (1) and Cap (2) using the model experimental system of FIG. In both cases, the horizontal axis is the electrophoresis time (arbitrary unit) and the vertical axis is the fluorescence intensity (arbitrary unit).
  • Figures 11 (a) and 11 (b) are the same raw data obtained at the emission point P (1), and (b) is an enlargement of the vertical scale of (a).
  • FIGS. 11 (c) and 11 (d) are the same raw data obtained at emission point P (2), and FIG. 11 (d) is an enlargement of the vertical scale of FIG. 11 (c). It is a thing.
  • FIGS. 11 (a) to 11 (d) indicate times 0 to 500
  • the vertical axis scales in FIGS. 11 (a) and 11 (c) are fluorescence intensities 0 to 250
  • the axis scale is fluorescence intensity -0.2 to 0.8
  • X (1,1) and X (2,1) are represented by solid lines
  • X (1,2) and X (2,2) are represented by dotted lines.
  • Z (1,1), Z (1,2), Z (2,1), and Z (2,2) are the phosphors D (1,1), D (1,2) at each time. ), D (2,1), and D (2,2).
  • the 2-by-2 matrix Y and the inverse matrix Y -1 are the same at the emission point P (1) and the emission point P (2). It was. However, in general, the matrix Y and the inverse matrix Y -1 may differ between the light emitting point P (1) and the light emitting point P (2).
  • Each element of the matrix Y for the emission point P (1) is the fluorescence intensity ratio of X (1,1) to X (1,2) by the phosphor D (1,1) at time 100, and the fluorescence at time 300. It was determined by the intensity ratio of fluorescence intensity X (1,1) to X (1,2) by body D (1,2).
  • Each element of the matrix Y for the emission point P (2) is the fluorescence intensity ratio of the fluorescence intensities X (2,1) and X (2,2) by the phosphor D (2,1) at time 200, and the fluorescence at time 400. It was determined by the intensity ratio of the fluorescence intensities X (2,1) and X (2,2) by the body D (2,2).
  • FIG. 12 is color conversion data obtained by performing color conversion of (Equation 18) and (Equation 19) at each time on the raw data of FIG.
  • the notation method is the same as in FIG.
  • the spectral crosstalk of the four large peaks in FIGS. 11 (a) and 11 (c) is resolved by the four large peaks in FIGS. 12 (a) and 12 (c).
  • the spectral crosstalk of the four small peaks in FIGS. 11 (b) and 11 (d) is also shown by the four small peaks in FIGS. 12 (c) and 12 (d) (also indicated by arrows).
  • these peaks themselves still remain, indicating that spatial crosstalk has not been resolved. This is the problem of the conventional method.
  • the matrix Y and the inverse matrix Y -1 are extended from 2 rows and 2 columns of (Equation 18) and (Equation 19) to 4 rows and 4 columns. Further, in (Equation 18) and (Equation 19), the color conversion was performed independently at the emission point P (1) and the emission point P (2), whereas in (Equation 20), the emission point P (1) was performed. ) And the light emitting point P (2) are mixed and color conversion + spatial correction is performed at once.
  • Each element of the matrix Y of (Equation 20) has fluorescence intensities X (1,1), X (1,2), X (2,1) by the phosphor D (1,1) at time 100 in FIG.
  • the upper left 2 rows and 2 columns of the matrix Y of (Equation 20) and the lower right 2 rows and 2 columns correspond to the matrices Y of (Equation 18) and (Equation 19), and the values of each element are almost the same as each other. equal.
  • the upper left 2 rows and 2 columns of the matrix Y -1 of (Equation 20) and the lower right 2 rows and 2 columns correspond to the matrices Y -1 of (Equation 18) and (Equation 19).
  • the values of each element are approximately equal to each other. In other words, these elements are responsible for color conversion that eliminates spectral crosstalk for the same emission point.
  • (Equation 21) Right-hand side 1st term and (Equation 22) Right-hand side 2nd term are responsible for the color conversion of the conventional method, and correspond to (Equation 18) and (Equation 19), respectively.
  • (Equation 21) the second term on the right side and (Equation 22) the first term on the right side are responsible for spatial correction and color conversion, which are not handled by the conventional method.
  • the color conversion, the spatial correction, and the color conversion may be separated and executed individually, or only one of them may be performed.
  • the processing can be changed depending on the light emitting point, for example, the light emitting point P (1) is subjected to spatial correction and color conversion, and the light emitting point P (2) is subjected to only color conversion.
  • FIG. 13 shows the color conversion + spatial correction data obtained by executing the color conversion + spatial correction of (Equation 20) at each time on the raw data of FIG.
  • the notation method is the same as in FIG.
  • the spectral crosstalk of the four large peaks in FIGS. 11 (a) and 11 (c) is resolved by the four large peaks in FIGS. 13 (a) and 13 (c). You can see that it is done.
  • the spatial and spectral crosstalks of the four small peaks in FIGS. 11 (b) and 11 (d) were also eliminated in FIGS. 13 (c) and 13 (d), resulting in peaks. Can be seen to have disappeared (indicated by the arrow).
  • spectral crosstalk for the same emission point and a plurality of emission points generated in an arbitrary optical system for detecting fluorescence of a plurality of types of phosphors emitted from a plurality of emission points. It was shown that spatial crosstalk and spectral crosstalk between them can be eliminated or reduced by computational processing.
  • FIG. 14 shows the raw data obtained as a result of performing electrophoretic analysis by injecting an unknown sample into the capillaries Cap (1) and Cap (2) using the model experimental system of FIG. Since the same model experimental system as in [Example 2] is used, (Equation 18) to (Equation 22) can be used as they are.
  • the notation method is the same as in FIG. However, the vertical scale in Figures 11 (b) and 11 (d) is reduced to fluorescence intensity -0.5 to 2.0. At time 100, 200, 300, and 400, a large peak was observed at emission point P (1), while a small peak (indicated by an arrow) was observed at emission point P (2).
  • FIG. 15 is color conversion data obtained by performing color conversion of (Equation 18) and (Equation 19) at each time on the raw data of FIG.
  • the notation method is the same as in FIG.
  • the emission fluorescence of the phosphor D (1,1) was detected independently at times 100 and 300, and the emission fluorescence of the phosphor D (1) was detected at times 200 and 400.
  • 2) Emission fluorescence was detected alone.
  • the identity of the four small peaks (indicated by arrows) detected at times 100, 200, 300, and 400 at the emission point P (2) is unknown. It was.
  • FIG. 16 shows the color conversion + spatial correction data obtained by executing the color conversion + spatial correction of (Equation 20) at each time on the raw data of FIG.
  • the notation method is the same as in FIG.
  • the emission fluorescence of the phosphor D (1,1) was independently detected at times 100 and 300, as in FIG. 15 (a).
  • Emission fluorescence of phosphor D (1,2) was detected alone at 200 and 400.
  • the fluorescence emission detection regions W (2,1) and W of the phosphors D (1,1) and D (1,2) at the emission point P (1) are shown in FIG. 16 (d).
  • the weak emission fluorescence of phosphor D (2,2) was independently detected at time 100 and 300 at the emission point P (2).
  • the weak emission fluorescence of phosphor D (2,1) was detected alone at times 200 and 400.
  • the peak intensity of these weak emission fluorescence was less than 1% of the peak intensity of the emission fluorescence observed at the emission point P (1).
  • the spatial crosstalk ratio generated in the model experimental system in Fig. 8 is also less than 1%. Therefore, in Fig. 15, the true weak emission fluorescence and the false weak emission fluorescence due to spatial crosstalk coexist, making it impossible to distinguish between the two.
  • the above results show that spatial crosstalk can push up the lower limit of detection in detecting luminescence from each emission point.
  • the spatial crosstalk ratio is 1%, even if the detection limit is 0.1% when there is one light emitting point, is the signal of 1% or less a true signal? Or, because it is not possible to distinguish whether it is a fake signal due to spatial crosstalk, the effective lower limit of detection rises to 1%.
  • the detection sensitivity and dynamic range are both orders of magnitude lower in the emission detection of multiple emission points than in the emission detection of one emission point. The present disclosure can solve such a problem and avoid a decrease in detection sensitivity and dynamic range in light emission detection of a plurality of light emitting points.
  • Fluorescent material D (1,1) at emission point P (1) at time 100, phosphor D (1,2) at emission point P (1) at time 200, and phosphor at emission point P (1) at time 300 Samples are prepared so that the phosphor D (1,3) fluoresces independently at the emission point P (1) at D (1,1) and time 400. Then, the sample is prepared so that the fluorescence other than these is not emitted. As shown in Fig. 17 (a), only four large peaks corresponding to the above fluorescence emission were observed. It was found that the spectral crosstalk ratio of fluorophore D (1,3) was similar to that of fluorophore D (1,2), but slightly different. On the other hand, the four small peaks shown by the arrows shown in Fig. 17 (d) are the results of the spatial crosstalk and spectral crosstalk of the emission fluorescence that are the sources of the four large peaks shown in Fig. 17 (a), respectively. It can be judged that.
  • FIG. 18 shows the color conversion + spatial correction data obtained by executing the color conversion + spatial correction of (Equation 20) at each time on the raw data of FIG.
  • the notation method is the same as in FIG.
  • the emission fluorescence of the phosphor D (1,1) was detected independently at time 100 and 300, and the emission fluorescence of the phosphor D (1,2) was detected at time 200.
  • Emission fluorescence was detected alone.
  • the emission fluorescence of phosphor D (1,3) was detected independently at time 400, but the spectral crosstalk ratios were between phosphor D (1,1) and phosphor D (1,2). ), So the spectral crosstalk was not resolved.
  • the phosphor D (1,1) fluoresces independently at the emission point P (1), and the concentration and fluorescence intensity of the phosphor D (1,1) are stepwise. Samples have been prepared to rise to. Then, the sample is prepared so that the fluorescence other than these is not emitted. As shown in FIG. 19 (a), four large peaks corresponding to the above fluorescence emission were observed, but the sensor S used in the model experimental system of FIG. 8 is saturated at the fluorescence intensity of 200. , Times 200, 300, and 400 were detected saturated.
  • FIG. 20 is color conversion data obtained by performing color conversion of (Equation 18) and (Equation 19) at each time on the raw data of FIG.
  • the notation method is the same as in FIG.
  • the emission fluorescence of the phosphor D (1,1) at time 100 at the emission point P (1) was detected independently after the spectral crosstalk was eliminated.
  • the spectral crosstalk was not resolved because the emission fluorescence of the phosphor D (1,1) at times 200, 300, and 400 was saturated. This is because, as is clear from Fig. 19 (a), when the fluorescence intensity is saturated, the spectral crosstalk ratio derived from the intensity ratio of X (1,1) and X (1,2) changes (number).
  • FIG. 21 is color conversion + spatial correction data obtained by executing (Equation 20) color conversion + spatial correction on the raw data of FIG. 19 at each time.
  • the notation method is the same as in FIG. In Fig. 21 (a), equivalent results are obtained for the same reason as in Fig. 20 (a).
  • Fig. 21 (d) both spatial crosstalk and spectral crosstalk were eliminated for the emission fluorescence of phosphor D (2,1) at time 100.
  • spatial crosstalk and spectral crosstalk were not resolved for the emission fluorescence of phosphor D (2,1) at times 200, 300, and 400. This is because the fluorescence intensity of X (1,1) or X (1,2) shown in FIG.
  • Example 6 Here, a method for determining the matrix Y of (Equation 20) will be described by a method different from [Example 2].
  • the phosphors D (1,1), D (2,1), D (1,2), and D (2,2) are used in this order.
  • the matrix Y was determined by fluorescing each of them independently, that is, by fluorescing them alternately at the emission point P (1) and the emission point P (2).
  • the electrophoresis speed of either Capillary Cap (1) or Cap (2) deviates from the assumption, fluorescence may be emitted simultaneously at the emission point P (1) and the emission point P (2). Can no longer determine the matrix Y.
  • Figure 22 shows a more realistic way to avoid the above problems.
  • the sample was injected only into the capillary Cap (1), electrophoretic analysis was performed with the capillaries Cap (1) and Cap (2), and the fluorophore D (1,1) was obtained. Fluorescence is emitted from each of D (1,2) independently, and fluorescence is detected in the detection regions W (1,1), W (1,2), W (2,1), and W (2,2).
  • the sample was injected only into the capillary Cap (2), and the electrophoretic analysis was performed with the capillaries Cap (1) and Cap (2) to obtain the phosphor D (2,1).
  • Fluorescence is emitted from each of D (2,2) independently, and fluorescence is detected in the detection regions W (1,1), W (1,2), W (2,1), and W (2,2).
  • the matrix Y can be determined as in [Example 2].
  • the matrix Y can be determined more easily and reliably than in [Example 2].
  • 2 and 2) are drawn so as to emit fluorescence multiple times, it is sufficient to emit fluorescence once for each.
  • the next electrophoresis analysis is performed after the first electrophoresis analysis is completed, but it is not necessary to extend the interval of the electrophoresis analysis to that extent.
  • the same effect as in FIG. 22 can be obtained.
  • the electrophoresis may be performed for a short time, then the sample may be injected into the capillary Cap (2), and then the electrophoresis may be restarted. This can reduce the time required to determine the matrix Y.
  • the sample injected into the capillary Cap (1) and the sample injected into the capillary Cap (2) can have the same composition or the same composition. This simplifies the sample to be prepared and reduces the cost for it.
  • FIG. 23 shows a method for further simplifying the determination of the matrix Y.
  • the phosphor D (1,1) in the capillary Cap (1) was similar to that in FIG. ) And D (1,2), and the phosphors D (2,1) and D (2,2) in the capillary Cap (2) are prepared so that they fluoresce at different times.
  • the composition of the sample injected into the capillary Cap (1) and the sample injected into the capillary Cap (2) may be changed so that the electrophoresis rate of the substance labeled with each phosphor is different.
  • a sample injected into Capillary Cap (1) has a 50-base-long DNA fragment labeled with fluorescent substance D (1,1) and a 60-base-long DNA fragment labeled with fluorescent substance D (1,2).
  • Samples containing DNA fragments and injected into Capillary Cap (2) were labeled with a 70-base-long DNA fragment labeled Fluorescent D (2,1) and Fluorescent D (2,2). It may be prepared so as to contain a DNA fragment having a length of 80 bases.
  • FIG. 23 shows that even if samples of the same composition are injected into the capillaries Cap (1) and Cap (2) at the same timing, different electrophoresis conditions are set for the capillaries Cap (1) and Cap (2). Can also be realized. For example, during electrophoresis, the applied voltage of Capi (2) is temporarily lowered, the temperature of Capi (2) during electrophoresis is lowered, and so on.
  • FIG. 24 is a block diagram of a multi-capillary electrophoresis device, which is an example of an analyzer.
  • the multi-capillary electrophoresis device is widely used as an analyzer for analyzing DNA sequences and DNA fragments.
  • the multi-capillary electrophoresis equipment includes capillary 24-1, cathode 24-4, anode 24-5, cathode side buffer solution 24-6, anode side buffer solution 24-7, and pump block 24-9. , Syringe 24-11, laser light source 24-12.
  • four capillaries 24-1 were used, and DNA sequences of different samples were performed in each capillary 24-1. Samples of DNA sequences consist of DNA fragments labeled with four fluorophore.
  • each capillary 24-1 near the light emitting point 24-14 is removed in advance, each capillary 24-1 near the light emitting point 24-14 is arranged on the same plane, and the laser beams 24-13 are collected. After illuminating, it was introduced along the array plane from the side of the above array plane. (5) Then, DNA fragments labeled with four types of phosphors are electrophoresed inside each capillary 24-1, and are excited by irradiation with a laser beam 24-13 as they pass through emission points 24-14. , It emitted fluorescence. In other words, four types of phosphors emitted fluorescence from the four emission points, and the fluorescence intensity of each changed from moment to moment with electrophoresis.
  • Figure 25 shows Raman emitting light from four emission points P (1) to P (4) due to laser beam irradiation when four capillary Caps (1) to Cap (4) are filled with a standard solution. This is the acquisition result of a two-dimensional sensor image including a wavelength dispersion image of scattered light.
  • the horizontal axis direction is the arrangement direction of the four capillaries, and the vertical axis direction is the wavelength direction.
  • the four vertically elongated streaky images indicated by the arrows are the wavelength dispersion images of the Raman scattered light of the capillaries Cap (1) to Cap (4), respectively.
  • the wavelength of the wavelength dispersion image on the two-dimensional sensor image was calibrated for each light emitting point, and the relationship between the pixel position and the wavelength was obtained.
  • B 20 different wavelength band detection regions are set for each of the emission points P (1) to P (4) on the two-dimensional sensor image. did.
  • Each of the 20 detection regions was set to detect light emission of 500 to 700 nm in an evenly divided wavelength band at 10 nm intervals.
  • the detection regions W (1,1), W (1,2), ..., W (1,20) set for the emission point P (1) are 500 to 510 nm, 510 to 520 nm, ..., respectively.
  • 690-700 nm wavelength band emission was set to be detected, and the signal intensities were set to X (1,1), X (1,2), ..., X (1,20). The same applies to the light emitting points P (2) to P (4).
  • the detection area set in 26 is also valid for multiple different emission detections or for multiple different analysis sessions.
  • the emission fluorescence of dR110, dR6G, dTAMRA, and dROX at the emission point P (1) is in the detection regions W (1,5), W (1,7), W (1,10), and W (1), respectively. , 12), it is mainly detected.
  • each emission fluorescence is also detected in the other detection regions of emission points P (1) by spectral crosstalk, and in the detection regions of emission points P (2) to P (3) by spatial crosstalk and spectral crosstalk. Is also detected with weak intensity.
  • dR110, dR6G, dTAMRA, and dROX at the emission point P (1) are set to D (1,1), D (1,2), D (1,3), and D (1,4), respectively.
  • Their concentrations were Z (1,1), Z (1,2), Z (1,3), and Z (1,4).
  • X is a matrix of 80 rows and 1 column
  • Y is a matrix of 80 rows and 16 columns
  • Z is a matrix of 16 rows and 1 column.
  • the method for determining the matrix Y is as described above. However, it is not always necessary to set all 80 x 16 elements. For example, an element whose absolute value is sufficiently small compared to other elements may be replaced with zero to simplify the related calculation. If the range covered by spatial crosstalk is limited, (Equation 9) to (Equation 12) may be defined within the limited range. In this case, a wider range of analysis can be performed by sequentially sliding the limited range, if necessary. For example, if the range of spatial crosstalk can be limited to adjacent capillaries, spatial crosstalk and spectral crosstalk for two or more distant capillaries may be ignored and related calculations may be omitted.
  • FIG. 27 shows the configuration of this optical system and a quadrant image of four light emitting points acquired by this optical system.
  • the optical system shown in Fig. 27 includes a condenser lens array 27-1, a long-pass filter 27-3, a dichroic mirror 27-4, 27-5, 27-6, and 27-7, and a two-dimensional sensor 27-8. As shown in FIGS.
  • each luminous flux 27-9 is incidentally incident on a set of dichroic mirror arrays consisting of four dichroic mirrors 27-4, 27-5, 27-6, and 27-7.
  • the divided luminous fluxes 27-10, 27-11, 27-12, and 27-13 are 520 to 550 nm, 550 to 580 nm, and 580 to 610 nm, respectively.
  • the arrangement directions of the four capillarys Cap (1) to Cap (4) that is, the arrangement directions of the light emitting points P (1) to P (4) and the division direction by the dichroic mirror array, that is, the wavelength directions are mutual. Since it is vertical, as shown in Fig. 27 (c), the divided images W (1,1) to W (4,4) are aligned in the two-dimensional sensor image 27-14 without overlapping each other. Was detected in. For example, the fluorescence emitted from the emission point P (1) is detected as a divided image of W (1,1), W (1,2), W (1,3), and W (1,4).
  • the divided images W (1,1) to W (4,4) in the two-dimensional sensor image 27-14 were set in the detection area, respectively.
  • the signal intensities detected in each detection area were set to X (1,1) to X (4,4), respectively.
  • the detection area set in Fig. 27 (c) is also valid for multiple different emission detections or for multiple different analysis sessions.
  • the DNA fragments of A and G were labeled. Since the maximum fluorescence emission wavelengths of dR110, dR6G, dTAMRA, and dROX are 541 nm, 568 nm, 595 nm, and 618 nm, respectively, the emission fluorescence is 520 to 550 nm, 550 to 580 nm, and 580 to, respectively.
  • the emission fluorescence of dR110, dR6G, dTAMRA, and dROX at the emission point P (1) is the detection regions W (1,1), W (1,2), W (1,3), and W (1), respectively. , 4) are mainly detected.
  • each emission fluorescence is also detected in the other detection regions of emission points P (1) by spectral crosstalk, and in the detection regions of emission points P (2) to P (3) by spatial crosstalk and spectral crosstalk. Is also detected.
  • dR110, dR6G, dTAMRA, and dROX at the emission point P (1) are the phosphors D (1,1), D (1,2), D (1,3), and D (1,4), respectively.
  • concentrations were Z (1,1), Z (1,2), Z (1,3), and Z (1,4). The same applies to the light emitting points P (2) to P (4).
  • X is a 16-by-1 matrix
  • Y is a 16-by-16 matrix
  • Z is a 16-by-1 matrix.
  • Figure 28 shows four types of multi-channel multi-channel reaction tanks arranged on a plane, in which a complementary strand extension reaction using a different DNA fragment in each channel as a template is performed in units of one base and incorporated into the complementary strand in the extension reaction.
  • a complementary strand extension reaction using a different DNA fragment in each channel as a template is performed in units of one base and incorporated into the complementary strand in the extension reaction.
  • the multi-channel extension reactor includes a laser light source 28-1, a dichroic mirror 28-3, a lens 28-4, a dichroic mirror 28-8, 28-9 and 28-10, a lens 28-11, and a first two-dimensional sensor 28. -12, lens 28-13, second 2D sensor 28-14, lens 28-15, third 2D sensor 28-16, lens 28-17 and 4th 2D sensor 28-18.
  • the maximum fluorescence emission wavelengths of the phosphors labeled with the base types T, C, A, and G were 535 nm, 565 nm, 595 nm, and 625 nm, respectively.
  • the laser beam 28-2 oscillated from the laser light source 28-1 was transmitted through the dichroic mirror 28-3, then focused by the lens 28-4 and irradiated to the multi-channel sample 28-5.
  • the fluorescence emission 28-6 of each phosphor excited by the laser beam irradiation was collectively collimated with the lens 28-4, and then the obtained luminous flux 28-7 was transferred to the dichroic mirror 28-3. It was reflected by.
  • the dichroic mirror 28-3 has a spectral characteristic that transmits laser beam light and reflects emission fluorescence.
  • the luminous flux 28-7 was divided into four luminous fluxes having four different wavelength components by three types of dichroic mirrors 28-8, 28-9, and 28-10.
  • the first divided light beam is imaged on the first two-dimensional sensor 28-12 by the lens 28-11, and the second divided light beam is formed on the second two-dimensional sensor 28-14 by the lens 28-13.
  • An image is formed, the third divided light beam is formed on the third two-dimensional sensor 28-16 by the lens 28-15, and the fourth divided light beam is formed on the fourth two-dimensional sensor 28-18 by the lens 28-17. The image was formed on the top.
  • Sample 28-19 is a schematic diagram of sample 28-5 observed from the front.
  • A 5 channels each formed emission points P (1) to P (5).
  • the emission point P (3) and the emission point P (3) of interest Is directly affected by spatial crosstalk, and only the emission points P (1), P (2), P (4), and P (5) existing around the emission point P (3) are analyzed. Therefore, only the emission points P (1) to P (5) are simply drawn on the sample 28-19.
  • the same analysis may be performed for the light emitting point and the light emitting points around it.
  • the first two-dimensional sensor image 28-20 is a divided image of the sample 28-19 acquired by the first two-dimensional sensor 28-12, and each of the light emitting points P (1) to P (5) is imaged.
  • the detection areas W (1,1) to W (5,1) were set at the points, and the respective signal intensities were set to X (1,1) to X (5,1).
  • fluorescent components in the wavelength band of 520 to 550 nm were detected.
  • the second two-dimensional sensor image 28-21 is a divided image of the sample 28-19 acquired by the second two-dimensional sensor 28-14, and each of the light emitting points P (1) to P (5) is imaged.
  • the detection areas W (1,2) to W (5,2) were set at the points.
  • the third two-dimensional sensor image 28-22 is a divided image of the sample 28-19 acquired by the third two-dimensional sensor 28-16, and each of the light emitting points P (1) to P (5) is imaged.
  • the detection areas W (1,3) to W (5,3) were set at the points.
  • fluorescent components in the wavelength band of 580 to 610 nm were detected.
  • the fourth two-dimensional sensor image 28-23 is a divided image of the sample 28-19 acquired by the fourth two-dimensional sensor 28-18, and each of the light emitting points P (1) to P (5) is imaged.
  • the detection areas W (1,4) to W (5,4) were set at the points.
  • the emission fluorescence of the four types of phosphors at the emission point P (3) is in the detection regions W (3,1), W (3,2), W (3,3), and W (3, 3, respectively, in order of wavelength. It is mainly detected in 4).
  • each emission fluorescence is also detected in the other detection regions of emission point P (3) by spectral crosstalk, and emission points P (1), P (2), P (4) by spatial crosstalk and spectral crosstalk.
  • P (5) is also detected in the detection area.
  • the four types of phosphors at the emission point P (3) are divided into D (3,1), D (3,2), D (3,3), and D (3,4), respectively, in order of wavelength. Their concentrations were Z (3,1), Z (3,2), Z (3,3), and Z (3,4). The same applies to the emission points P (1), P (2), P (4), and P (5).
  • X is a matrix of 20 rows and 1 column
  • Y is a matrix of 20 rows and 20 columns
  • Z is a matrix of 20 rows and 1 column.
  • FIG. 29 is a flowchart showing one analysis session of the conventional method.
  • the traditional analysis session is performed by an analytical system that has an analyzer, computer, display, and database that analyzes sample 29-1.
  • the analyzer has a sensor (not shown) that receives light from sample 29-1.
  • the sample 29-1 of type A (A is an integer of 2 or more) labeled with the phosphor of type C (C is an integer of 1 or more) is input to the analyzer.
  • each sample was analyzed in parallel by step 29-2, and the emission fluorescence was detected for each analysis and at each time in the wavelength band of type B (B is an integer of 1 or more).
  • a ⁇ B fluorescence intensity X (a, b) time-series raw data are acquired, and these time-series raw data are sent to the computer.
  • step 29-3 analyzed for (a 0), the matrix Y of C rows B string stored in the database - spectral crosstalk by color conversion using (a 0) (a 0) Is solved, and the concentration Z (a 0 , c) of the C type phosphor is obtained from the fluorescence intensity X (a 0 , b) at each time.
  • the time series color conversion data of Z (a 0 , c) is output in step 29-4. Steps 29-3 and 29-4 are performed for all a 0 . As shown in Fig.
  • the conventional method is characterized in that A analyzes and analyzes are performed independently.
  • the fluorescence intensity X (a, b) is expressed as X (b)
  • the concentration Z (a 0 , c) is expressed as Z (c).
  • FIG. 30 is a flowchart showing one analysis session of this method. Similar to FIG. 29, the analysis session of this example shown in FIG. 30 is performed by an analysis system having an analyzer, a computer, a display device, and a database for analyzing sample 29-1.
  • the analyzer has a sensor (not shown) that receives light from sample 29-1.
  • the sample 30-1 of type A (A is an integer of 2 or more) labeled with the phosphor of type C (C is an integer of 1 or more) is input to the analyzer.
  • each sample was analyzed in parallel by step 30-2, and emission fluorescence was detected for each analysis and at each time in the wavelength band of type B (B is an integer of 1 or more).
  • a ⁇ B fluorescence intensity X (a, b) time-series raw data are acquired, and these time-series raw data are sent to the computer.
  • Fig. 29 as shown by the dotted arrow in step 30-2, there is a significant spatial crosstalk of luminescence fluorescence between different analyzes.
  • the step 30-3, all analyzed (a), the matrix of which is stored in the database (A ⁇ C) line (A ⁇ B) column Y - color conversion + space correction using Spectral crosstalk and spatial crosstalk are eliminated by, and at each time, the fluorescence intensity X (a, b) is combined with the concentration Z (a, c) of the C type phosphor for the A type sample or analysis.
  • the time-series color conversion + spatial correction data of Z (a, c) is output for each analysis (a) in step 30-4.
  • this method is characterized in that the color conversion + spatial correction in step 30-3 is performed collectively for all analyzes.
  • common matrix Y color conversion + space correction to all analyzes - is also characteristic points used.
  • Figure 31 is performed prior to FIG. 30, the matrix Y stored in the database - is a flowchart illustrating a method for obtaining the.
  • the above is realized by the sample 31-1 for determining the matrix Y, but the above may be realized by the setting of the analyzer.
  • the analyzer in step 31-2, the above A ⁇ C single fluorescent emission was detected in the wavelength band of B type (B is an integer of 1 or more), and A ⁇ B fluorescence intensity X.
  • the time-series raw data of (a, b) is acquired, and these time-series raw data are sent to the computer.
  • a matrix of (A ⁇ B) rows (A ⁇ C) columns is derived from the (A ⁇ B) ⁇ (A ⁇ C) data obtained in step 31-3.
  • FIG. 32 is a flowchart in which the analysis sessions of steps 30-2 to 30-4 for the A type sample 30-1 shown in FIG. 30 are repeated a plurality of times while utilizing the matrix Y ⁇ obtained in FIG. 31. is there.
  • the A-type samples 30-1 to be analyzed in each analysis session are different from each other.
  • the same matrix Y stored in the database - that are utilizing a is a feature of the present method.
  • Figure 33 shows an example of a computer configuration.
  • the calculator is connected to the analyzer.
  • the computer not only analyzes the data but also controls the analyzer.
  • the database and display device are drawn outside the computer, but in Fig. 33, they are drawn inside the computer.
  • Data analysis conditions and analyzer control conditions are set through the keyboard, which is the input unit.
  • the time-series raw data of the fluorescence intensity X (a, b) output from the analyzer is sequentially stored in the memory. Further, it stored in the database internal to the HDD (A ⁇ C) line (A ⁇ B) column of matrix Y - is stored in the memory.
  • the CPU the fluorescence intensity X (a, b) stored in the memory and Y - the product calculated in the phosphor concentration Z (a, c) deriving a time-series color conversion + space correction data, sequentially memory
  • the analysis results can be collated with the information on the network through the network interface NIF.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Engineering & Computer Science (AREA)
  • Immunology (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Organic Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Molecular Biology (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Biotechnology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Genetics & Genomics (AREA)
  • General Engineering & Computer Science (AREA)
  • Microbiology (AREA)
  • Biomedical Technology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Sustainable Development (AREA)
  • Medicinal Chemistry (AREA)
  • Biophysics (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

複数の発光点からの，複数種類の蛍光体の発光を識別するため，各発光点からの蛍光をそれぞれ複数の波長帯で検出する分析方法および分析装置において，複数の発光点および複数の波長帯の信号強度の間に空間クロストークおよびスペクトルクロストークが存在し，上記識別の性能が低下する。予め定められた演算式に，複数の発光点それぞれの，複数の波長帯での検出信号のすべてを入力することによって，空間クロストークおよびスペクトルクロストークを解消し，各発光点における各発光体の濃度を導出する。

Description

分析装置および分析方法
　本開示は，複数の発光点において，複数種類の蛍光体から発光される蛍光を，それぞれ識別しながら検出する分析方法および分析装置に関する。
　近年，複数のキャピラリに，電解質溶液，あるいは高分子ゲルやポリマを含む電解質溶液等の電気泳動分離媒体を充填し，並列に電気泳動分析を行うマルチキャピラリ電気泳動装置が広く用いられている。分析対象は，低分子から，タンパク質，核酸等の高分子まで，幅広い。また，計測モードには，ランプ光を各キャピラリの吸光点に照射し，分析対象が吸光点を通過する際に生じるランプ光の吸収を検出するモード，あるいは，レーザ光を各キャピラリの発光点に照射し，分析対象が発光点を通過する際に生じる蛍光あるいは散乱光を検出するモード等，多数ある。
　例えば，特許文献1では，A本（Aは2以上の整数）のキャピラリ上のA個の発光点の周辺の全キャピラリを同一平面上に配列し，配列平面の側方よりレーザビームを導入して全キャピラリの発光点を一括照射し，配列平面に垂直方向より各発光点で発生する蛍光を波長分散させて一括検出している。検出装置では，A個の発光点から発光される蛍光を1個の集光レンズで一括してコリメートし，1個の透過型の回折格子を透過させ，各蛍光の1次回折光を1個の結像レンズで1個の2次元センサ上に一括して結像させている。ここで，A個の発光点の配列方向と，回折格子による波長分散方向が互いに垂直になるようにすることにより，2次元センサ上で各キャピラリからの発光蛍光の波長分散像が互いに重なり合わないようにしている。各キャピラリの波長分散像について，B個（Bは1以上の整数）の任意の波長帯の検出領域を設定することでB色検出が可能になる。B＝1の場合を単色検出，B≧2の場合を多色検出と呼ぶ。特許文献1のマルチキャピラリ電気泳動装置では，例えば，各キャピラリで異なるDNAサンプルのサンガー法によるDNAシーケンスを行うことができる。サンガー法では，DNAサンプルに含まれるDNA断片に，末端塩基種A，C，G，およびTに応じて，4種類の蛍光体を標識し，それぞれの発光蛍光を多色検出によって識別している。
　特許文献2では，A本（Aは2以上の整数）のキャピラリ上のA個の発光点の周辺の全キャピラリを同一平面上に配列し，配列平面の側方よりレーザビームを導入して全キャピラリの発光点を一括照射し，配列平面に垂直方向より各発光点で発生する蛍光を波長成分に応じて分割させて一括検出している。検出装置では，A個の発光点からの発光蛍光をそれぞれA個の集光レンズで個別にコリメートしてA本の光束とし，B個（Bは1以上の整数）のダイクロイックミラーを配列した1組のダイクロイックミラーアレイに各光束を並列に入射してそれぞれをB本の異なる波長帯の光束に分割し，生成された合計A×B本の光束を1個の2次元センサに並列に入射し，A×B個の分割像を画像上に生成する。ここで，A個の発光点の配列方向と，ダイクロイックミラーアレイによる各光束のB個の分割方向が互いに垂直になるようにすることにより，A×B個の分割像が画像上で互いに重なり合わないようになり，A×B個の検出領域を設定できるようになる。これにより，各キャピラリのB色検出が可能になる。したがって，特許文献2のマルチキャピラリ電気泳動装置では，例えば，特許文献1の場合と同様に，各キャピラリで異なるDNAサンプルのサンガー法によるDNAシーケンスを行うことができる。
　しかしながら，一般に，多色検出を行っただけでは，複数種類の蛍光体の発光蛍光を識別することができない。何故ならば，各蛍光体の蛍光スペクトルは互いに重なり合っているため，任意のひとつの波長帯に，複数種類の蛍光体の蛍光が入り混じる（本開示では，スペクトルクロストークと呼ぶ）からである。また，異なる濃度の，複数種類の蛍光体が同時に蛍光を発光する場合があるからである。そこで，次の工程（本開示では，色変換と呼ぶ）によって，スペクトルクロストークを解消し，上記識別を可能にしている。
　A個（Aは2以上の整数）の発光点のそれぞれについて，C種類（Cは1以上の整数）の蛍光体の発光蛍光を，B種類（Bは1以上の整数）の波長帯の検出領域でB色検出を各時刻で行う。ただし，B≧Cとする。各発光点について，かつ各時刻について，B色検出結果に対して色変換を施し，C種類の蛍光体の濃度を取得する。発光点P(a)（a＝1，2，…，A）それぞれについて，異なる波長帯の検出領域W(b)（b＝1，2，…，B）で，蛍光体D(c)（c＝1，2，…，C）の発光蛍光を検出する。任意の時刻において，発光点P(a)における蛍光体D(c)の濃度をZ(c)，発光点P(a)についての検出領域W(b)の信号強度をX(b)とする。ここで，Xが信号強度X(b)を要素とするB行1列の行列，Zが濃度Z(c)を要素とするC行1列の行列，YがY(b)(c)を要素とするB行C列の行列として，次式が成り立つ。（数1）～（数4）は，bおよびcの関係式であるが，aの関係式ではなく，各発光点P(a)で独立に成立する。B＝1の単色検出の場合は，B≧Cにより，C=1となり，X，Y，Zはいずれも行列ではなくなる。
Figure JPOXMLDOC01-appb-M000007
Figure JPOXMLDOC01-appb-M000008
Figure JPOXMLDOC01-appb-M000009
Figure JPOXMLDOC01-appb-M000010
　ここで，B行C列の行列Yの要素Y(b)(c)は，スペクトルクロストークによって，蛍光体D(c)の発光蛍光が，異なる波長帯の検出領域W(b)で検出される信号強度比率を表す。いずれか1種類の蛍光体D(c0)を単独に蛍光発光させることにより，行列Yの1列Y(b)(c0)（b＝1，2，…，B）を決定することができる。ここで，蛍光体D(c0)の濃度を制御することは一般に困難であるため，1列Y(b)(c0)を規格化すると便利である。例えば，B個の要素の内，最大の要素を1として，他の要素を最大値に対する比率で示すのが良い。あるいは，B個の要素の合計が1になるように，各要素の比率を決めるのが良い。つまり，
Figure JPOXMLDOC01-appb-M000011
 とする。そして，上記工程をC種類の蛍光体D(c)すべてについて順番に行うことによって行列Yの全列を決定することができる。行列Yは，蛍光体D(c)および異なる波長帯の検出領域W(b)の特性のみで決まり，電気泳動分析の最中に変化しない。また，光学系，蛍光体D(c)，検出領域W(b)，等の条件を固定している限り，異なる電気泳動分析についても行列Yは一定に保たれる。したがって，各発光点について，各時刻における蛍光体D(c)の濃度Z(c)は，各時刻における検出領域W(b)の信号強度X(b)から，次式によって求められる。
Figure JPOXMLDOC01-appb-M000012
 　ここで，Y－はC行B列の，Yの一般逆行列であり，Y－＝(YT×Y)－1×YT)によって求められる。行列YがB＝Cの正方行列の場合は，Y－は逆行列Y－1と等しい。
　（数1）は，未知であるC種類の蛍光体の濃度と，既知であるB色蛍光強度の関係を示す連立方程式であり，（数6）はその解を求めることに相当する。したがって一般には，上述の通り，B≧Cの条件が必要である。仮に，B＜Cであると，解を一意に求めることができないため（つまり，複数の解が存在し得るため），（数6）のように色変換を実行できない。
　例として，サンガー法のC=4，4色検出のB=4の場合について詳細に説明する。サンガー反応により，鋳型DNAに対する様々な長さのDNA断片のコピーを作製する際に，末端の塩基種A，C，G，およびTに応じて4種類の蛍光体D(1)，D(2)，D(3)，およびD(4)で標識し，これらを電気泳動によって長さ分離しながら順番に，レーザビームを照射して蛍光発光させる。蛍光体D(1)，D(2)，D(3)，およびD(4)のそれぞれの極大発光波長に合わせた4種類の波長帯の検出領域W(1)，W(2)，W(3)，およびW(4)で発光蛍光を4色検出し，それらの信号強度X(1)，X(2)，X(3)，およびX(4)の時系列データ（本開示では，生データと呼ぶ。ただし，B＝4，C=4の場合に限らない）を取得する。各時刻における蛍光体D(1)，D(2)，D(3)，およびD(4)の濃度をZ(1)，Z(2)，Z(3)，およびZ(4)とすると，（数1）は次式となる。
Figure JPOXMLDOC01-appb-M000013
　ここで，4行4列の行列Yの要素Y(b)(c)は，スペクトルクロストークによって，蛍光体D(c)（cは1，2，3，または4）の発光蛍光が波長帯W(b)（bは1，2，3，または4）で検出される強度比率を表す。各キャピラリで，蛍光体D(c)（cは1，2，3，または4）がそれぞれ単独に蛍光発光するサンプルを電気泳動分析することにより，行列Yの要素Y(b)(c)を決定することができる。例えば，蛍光体D(1)が単独で蛍光発光している時の4色蛍光強度X(1)，X(2)，X(3)，およびX(4)はそれぞれ，要素Y(1)(1)，Y(2)(1)，Y(3)(1)，およびY(4)(1)を与える。また，蛍光体D(2)が単独で蛍光発光している時の4色蛍光強度X(1)，X(2)，X(3)，およびX(4)はそれぞれ，要素Y(1)(2)，Y(2)(2)，Y(3)(2)，およびY(4)(2)を与える。蛍光体D(3)，D(4)についても同様である。Y(b)(c)は，蛍光体D(c)および波長帯W(b)の特性のみで決まる固定値であり，電気泳動中に変化しない。したがって，各キャピラリについて，各時刻における蛍光体D(1)，D(2)，D(3)，およびD(4)の濃度は，各時刻における4色蛍光強度X(1)，X(2)，X(3)，およびX(4)から，（数6）を具体化した次式によって求められる。
Figure JPOXMLDOC01-appb-M000014
　このように，逆行列Y－1を，4色蛍光強度に乗じることによってスペクトルクロストークを解消し，4種類の蛍光体の濃度，すなわち末端が4種類の塩基のDNA断片の濃度の時系列データ（本開示では，色変換データと呼ぶ。ただし，B＝4，C=4の場合に限らない）を取得する。
　以上の色変換の工程は，既に述べている通り，A本のキャピラリそれぞれについて独立に行う。これは，キャピラリ間のクロストーク（本開示では，空間クロストークと呼ぶ）が十分に小さいという前提の上で成り立つ。つまり，任意のひとつのキャピラリから発光される蛍光を検出して得られる信号強度に，他のキャピラリから発光される蛍光に由来する信号強度が混在する比率が十分に小さいことを意味している。したがって，上述した，1種類の蛍光体D(c0)を単独に蛍光発光させる工程をC種類の蛍光体D(c)について順番に行うことによって行列Yを決定する工程は，A本のキャピラリについて実質的に同時に行い，A本のキャピラリそれぞれの行列Yを並列に求めている。これは，行列Yの導出を短時間に，かつ手間をかけずに実行するために必然的な工程である。
　空間クロストークは，光学系の工夫によって低減することが基本であるが，計算処理によって低減する試みもある。特許文献3では，複数のキャピラリではなく，平面上にランダムに配列した複数の発光点からの発光蛍光を1個の集光レンズによって1個の2次元センサ上に一括して結像させている。各発光点からの発光蛍光は，2次元センサ上のそれぞれの結像位置に設けられた検出領域の信号強度として得られる。任意のひとつの発光点からの発光蛍光の信号強度が，他の発光点からの発光蛍光の信号強度に与える空間クロストークの比率は，それら2つの発光点間の距離，あるいは対応する2つの検出領域の距離の関数として表すことができ，上記距離とともに減少することを見出している。上記関数を予め求めておき，平面上にランダムに配列した複数の発光点の蛍光画像において，任意の2つの検出領域の距離と上記関数から相互の空間クロストークを求め，元の蛍光画像から，求めた空間クロストークを差引くことで，空間クロストークの低減を図っている。
特許第3897277号公報 特許第6456983号公報 特表2018-529947号公報
　特許文献1，および特許文献2のいずれの場合も，A本のキャピラリ上のA個の発光点から発光した蛍光は，2次元センサ上で互いに分離して結像されるため，空間クロストークは本質的に低く抑えられている。空間クロストークの要因には，（1）レンズの収差，（2）マルチキャピラリ，レンズ，フィルタ，ダイクロイックミラー，および2次元センサ等の要素から構成される光学系の内部で生じる，蛍光の要素間の多重反射，（3）2次元センサの画素間のブルーミング，等が考えられる。空間クロストークを低減するには，上記の要因（1）～（3）の影響が最小になるように光学系を構築する必要がある。例えば，要因（2）の影響を抑えるため，レンズに施す減反射コートとして，より反射率の小さいものを選択する等が考えられる。しかしながら，空間クロストークを完全にゼロにすることはできない。空間クロストークは，任意のひとつのキャピラリから発光される蛍光の検出において，実効的に，検出下限を上昇させるため，検出感度を低下させると同時に，検出ダイナミックレンジを狭める可能性がある。したがって，空間クロストークを少しでも低減することは，上記の課題を緩和，または解決する上で極めて重要である。
　本発明者らは、特許文献1，または特許文献2の光学系について，特許文献3の方法によって，空間クロストークを低減させることを試みたが，良好に機能しなかった。まず，B＝1の単色検出の場合について説明する。A個の発光点の任意の1個の発光点である発光点αについての検出領域βの信号強度X(α)(β)が，上記以外の任意の1個の発光点α'（α≠α'）についての検出領域βの信号強度X(α')(β)に与える空間クロストークが，上記2個の検出領域の間の距離とともに減少する傾向が認められたものの，ひとつの上記距離の関数で表すことが出来ないことが明らかになった。例えば，上記2個の検出領域の間の距離が一定であっても，上記2個の検出領域が，光学系の中心軸の近くに位置する場合と，光学系の中心軸から遠くに位置する場合で，空間クロストークの比率が異なった。また，上記2個の検出領域の一方から他方への空間クロストークの比率と，他方から一方への空間クロストークの比率が異なった。B≧2の多色検出の場合は，特許文献3の方法が一層機能しなくなった。A個の発光点の任意の1個の発光点である発光点αについての，B個の検出領域の任意の1個の検出領域である検出領域βの信号強度X(α)(β)が，上記以外の任意の1個の発光点α'（α≠α'）についての，B個の検出領域の任意の1個の検出領域β'（β＝β'またはβ≠β'であり，β＝β'は同じ波長帯の検出領域，β≠β'は異なる波長帯の検出領域を意味する）の信号強度X(α')(β')に与える空間クロストークが，上記2個の検出領域の間の距離の関数で表すことが出来なかった。特に，β≠β'の場合，空間クロストークとスペクトルクロストークが混在するため，上記2個の検出領域の間のクロストークは，本質的に，両者の距離の関数で表すことが不可能であった。本課題については，［発明を実施するための形態］で詳細に説明する。
　そこで，本開示では，複数の発光点からの発光をそれぞれ検出する任意の光学系において発生する，複数の発光点間の空間クロストークを計算処理によって低減することによって，複数の発光点からの発光を識別し，それぞれを独立に検出する方法を提案する。また，複数の発光点から発光される，複数種類の蛍光体の蛍光を，それぞれ検出する任意の光学系において発生する，複数の発光点間の空間クロストーク，および各発光点についての複数種類の蛍光体間のスペクトルクロストークを計算処理によって低減することによって，複数の発光点からの，複数種類の蛍光体の蛍光発光を識別し，それぞれを独立に検出する方法を提案する。あるいは，複数の吸光点における吸光をそれぞれ検出する任意の光学系において発生する，複数の吸光点間の空間クロストークを計算処理によって低減することによって，複数の吸光点における吸光を識別し，それぞれを独立に検出する方法を提案する。また，複数の吸光点における，複数種類の吸光体の吸光を，それぞれ検出する任意の光学系において発生する，複数の吸光点間の空間クロストーク，および各吸光点についての複数種類の吸光体間のスペクトルクロストークを計算処理によって低減することによって，複数の吸光点における，複数種類の吸光体の吸光を識別し，それぞれを独立に検出する方法を提案する。
　特許文献1または特許文献2等の光学系において，A個（Aは2以上の整数）の発光点それぞれについて，B個（Bは1以上の整数）の波長帯の検出領域で，C種類（Cは1以上の整数）の蛍光体の発光蛍光を検出する場合について説明する。ここで，異なる波長帯の検出領域はそれぞれ，蛍光の異なる波長成分を検出するものとする。また，異なる種類の蛍光体はそれぞれ，異なる蛍光スペクトルを有する蛍光を発光するものとする。発光点P(a)（a＝1，2，…，A）それぞれについて，異なる波長帯の検出領域W(a,b)（b＝1，2，…，B）において，発光点P(a)に存在する蛍光体D(a,c)（c＝1，2，…，C）の発光蛍光を検出する。任意の時刻における，発光点P(a)における蛍光体D(a,c)の濃度をZ(a,c)，発光点P(a')についての検出領域W(a',b)の信号強度をX(a',b)とする。ここで，XがX(a',b)を要素とするA×B行1列の行列，ZがZ(a,c)を要素とするA×C行1列の行列，YがY(a',b)(a,c)を要素とする(A×B)行(A×C)列の行列として，次式が成り立つことを本発明者らは初めて見出した。
Figure JPOXMLDOC01-appb-M000015
Figure JPOXMLDOC01-appb-M000016
Figure JPOXMLDOC01-appb-M000017
Figure JPOXMLDOC01-appb-M000018
　（数9）は（数1）と同じ式であるが，（数1）～（数4）と（数9）～（数12）を比較すると，両者が全く異なることが分かる。（数9）～（数12）は，bおよびcの関係式であるだけでなく，aの関係式であり，異なる発光点P(a)が相互に関連している。すなわち，（数9）～（数12）は，同一発光点P(a)についてのスペクトルクロストークに加えて，異なる発光点P(a)間の空間クロストークおよびスペクトルクロストークを一括して考慮できるようにしたものであり，（数1）～（数4）と本質的に異なる。
　ここで，(A×B)行(A×C)列の行列Yの要素Y(a',b)(a,c)は，（1）a'＝a，すなわち同じ発光点についてのスペクトルクロストークによって，P(a')における蛍光体D(a',c)の発光蛍光が，任意の検出領域W(a',b)で検出される信号強度比率を表すことに加えて，（2）a'≠a，すなわち異なる発光点についての空間クロストークおよびスペクトルクロストークによって，発光点P(a)における蛍光体D(a,c)の発光蛍光が，P(a')の任意の検出領域W(a',b)で検出される信号強度比率を表す。いずれか1個の発光点P(a0)において，いずれか1種類の蛍光体D(a0,c0)を単独に蛍光発光させることにより，行列Yの1列Y(a,b)(a0,c0)（a＝1，2，…，A，および，b＝1，2，…，B）を決定することができる。ここで，蛍光体D(a0,c0)の濃度を制御することは一般に困難であるため，1列Y(a,b)(a0,c0)を規格化することが便利である。例えば，A×B個の要素の内，最大の要素を1として，他の要素を最大値に対する比率で示すのが良い。あるいは，A×B個の要素の合計が1になるように，各要素の比率を決めるのが良い。つまり，
Figure JPOXMLDOC01-appb-M000019
 とする。そして，上記工程をA個の発光点P(a)におけるC種類の蛍光体D(a,c)のすべての組み合わせについて順番に行うことによって行列Yの全列を決定することができる。行列Yは，発光点P(a)，蛍光体D(a,c)，および検出領域W(a,b)の特性のみで決まり，電気泳動分析の最中に変化しない。また，光学系，発光点P(a)，蛍光体D(a,c)，検出領域W(a,b)，等の条件を固定している限り，異なる電気泳動分析についても行列Yは一定に保たれる。したがって，各発光点について，各時刻における蛍光体D(a,c)の濃度Z(a,c)は，各時刻における検出領域W(a,b)の信号強度X(a,b)から，次式によって求められる。
Figure JPOXMLDOC01-appb-M000020
　（数14）は（数6）と同じ式であるが，（数2）～（数4）および（数6）と，（数10）～（数12）および（数14）を比較すると，両者が全く異なることが分かる。（数14）によって，各時刻について，予め求めた行列Yの一般逆行列Y－と，A×B個の信号強度X(a,b)を掛け合わせることによって，スペクトルクロストークと空間クロストークの両方を解消して，A×C個の蛍光体の濃度Z(a,c)を求めることができる。Y－は(A×C)行(A×B)列の行列である。本開示では，スペクトルクロストークを解消させることを色変換と呼ぶことに対応させて，空間クロストークを解消させることを空間補正と呼ぶ。つまり，Y－の(A×C)×(A×B)個の要素の中には，色変換を実行する部分と，空間補正を実行する部分の両方が含まれている。そして，各時刻について（数14）を実行して得られるZ(a,c)の時系列データを，本開示では，色変換+空間補正データと呼ぶ。
　行列Xの各要素は，予め背景光を差引いた値としても良いし，適当なノイズ低減処理を施した値としても良い。また，行列Yの各要素についても，同様に，予め背景光を差引いた値としても良いし，適当な変化を加えても構わない。
　B＝1，およびC=1の単色検出の場合，上記は次のように簡略化される。発光点P(a)（a＝1，2，…，A）それぞれについて，1種類の蛍光体の発光蛍光を1個の検出領域で検出する。任意の時刻における，発光点P(a)における蛍光体の濃度をZ(a)，発光点P(a')についての信号強度をX(a')とする。ここで，XがX(a')を要素とするA行1列の行列，ZがZ(a)を要素とするA行1列の行列，YがY(a')(a)を要素とするA行A列の行列とすると，（数10）～（数12）は次式のように簡略化される。（数9）および（数14）はそのまま成立する。
Figure JPOXMLDOC01-appb-M000021
Figure JPOXMLDOC01-appb-M000022
Figure JPOXMLDOC01-appb-M000023
　（数9）および（数14）～（数17）は，特許文献3で用いられている一般式，例えば，特許文献3の［数7］～［数9］と形式的には同じであるが，内容は大きく異なる（以降，特許文献3で用いられている数式を［　］で記す）。［数7］の行列Aの要素αijは，発光点φ(meas)iとφ(meas)jそれぞれの検出領域の間の距離dijの関数であり，具体的には［数10］に示されている通り，dijの指数関数の和で表され，dijの増大に伴ってαijが減衰する。特許文献3では，異なるサンプルに対して異なる蛍光画像が得られ，各蛍光画像において複数の発光点の検出領域の位置がランダムに変化する。そこで特許文献3では，予め上記関数を求めておき，各蛍光画像について，複数の発光点の内の任意の2個の発光点の検出領域の間の距離をすべての2個の発光点の組み合わせについて求め，上記関数に代入してαijを導出している。つまり，サンプル毎，あるいは蛍光画像毎にαij，すなわち行列Aが変化する。これに対して本開示では，上述した通り，あるいは［発明を実施するための形態］で後述する通り，（数16）の行列Yの要素Y(a')(a)は，発光点P(a)とP(a')それぞれの検出領域の間の距離の関数で表せないことを見出している。また，光学系の構成等から，要素Y(a')(a)を計算で求めることも不可能であることを見出している。そこで本開示では，異なるサンプルに対しても，複数の発光点の検出領域の位置が変化しない条件として，その条件下で要素Y(a')(a)を実測で求めている。具体的には，既に述べている通り，いずれか1個の発光点P(a0)においてのみ，単独に蛍光発光させることにより，行列Yの1列Y(a')(a0)（a'＝1，2，…，A）を決定する。そして，上記工程をA個の発光点P(a)のすべてについて順番に行うことによって行列Yの全列を決定する。行列Yは，発光点P(a)の特性のみで決まり，分析の最中に変化しない。また，光学系，発光点P(a)とその検出領域，等の条件を固定している限り，異なるサンプルの分析についても行列Yは一定に保たれる。以上の工程を，特許文献3の光学系，あるいは装置構成で実施することは不可能である。何故ならば，複数の発光点の検出領域の位置を固定することができないし，仮に固定できたとしても，各発光点を単独で，かつ順番に発光させることができないからである。したがって，特許文献3から，本開示の方法を想到することはできない。
　以上は次のように読み替えることができる。すなわち，A個（Aは2以上の整数）の発光点それぞれについて，B個（Bは1以上の整数）の任意の波長帯の検出領域で，C種類（Cは1以上の整数）の発光体の発光を検出する。発光点P(a)（a＝1，2，…，A）それぞれについて，異なる波長帯の検出領域W(a,b)（b＝1，2，…，B）において，発光点P(a)に存在する発光体D(a,c)（c＝1，2，…，C）の発光を検出する。任意の時刻における，発光点P(a)における発光体D(a,c)の濃度をZ(a,c)，P(a')についてのW(a',b)の発光強度をX(a',b)とする。この場合も，（数9）～（数17）が成立し，同様に，スペクトルクロストークと空間クロストークの両方を解消して，A×C個の発光体の濃度Z(a,c)を求めることができる。ここで，発光とは，蛍光，燐光，散乱光などが当てはまる。
　また，以上は次のように読み替えても良い。すなわち，A個（Aは2以上の整数）の吸光点それぞれについて，B個（Bは1以上の整数）の任意の波長帯の検出領域で，C種類（Cは1以上の整数）の吸光体の吸光を検出する。吸光点P(a)（a＝1，2，…，A）それぞれについて，異なる波長帯の検出領域W(a,b)（b＝1，2，…，B）において，吸光点P(a)に存在する吸光体D(a,c)（c＝1，2，…，C）の吸光を検出する。任意の時刻における，吸光点P(a)における吸光体D(a,c)の濃度をZ(a,c)，吸光点P(a')についてのW(a',b)の吸光度をX(a',b)とする。この場合も，（数9）～（数17）が成立し，同様に，スペクトルクロストークと空間クロストークの両方を解消して，A×C個の吸光体の濃度Z(a,c)を求めることができる。
　あるいは，以上を光計測以外の多点検出に読み替えることもできる。すなわち，A個（Aは2以上の整数）の信号発生点それぞれについて，B個（Bは1以上の整数）の任意の周波数帯の検出領域で，C種類（Cは1以上の整数）の信号発生体の信号を検出する。信号発生点P(a)（a＝1，2，…，A）それぞれについて，異なる周波数帯の検出領域W(a,b)（b＝1，2，…，B）において，信号発生点P(a)に存在する信号発生体D(a,c)（c＝1，2，…，C）の信号を検出する。任意の時刻における，信号発生点P(a)における信号発生体D(a,c)の濃度をZ(a,c)，信号発生点P(a')についてのW(a',b)の信号強度をX(a',b)とする。この場合も，（数9）～（数17）が成立し，同様に，スペクトルクロストークと空間クロストークの両方を解消して，A×C個の信号発生体の濃度Z(a,c)を求めることができる。
　以上では，主として数式を用いた説明を行ったが，これは本開示の内容の理解を容易にするためのものである。本開示の技術を実施する際には，本開示の内容を踏まえた方法を用いれば良く，必ずしも数式の通りに行う必要はなく，数式を変形しても構わないし，数式を用いなくても構わない。また，本開示において，「発光体の濃度」との記載は「発光体の発光強度」と読み替えることができ，「蛍光体の濃度」との記載は「蛍光体の蛍光強度」と読み替えることができ，「吸光体の濃度」との記載は「吸光体の吸光度」と読み替えることができる。
　本開示に関連する更なる特徴は、本明細書の記述、添付図面から明らかになるものである。また、本開示の態様は、要素及び多様な要素の組み合わせ及び以降の詳細な記述と添付される請求の範囲の様態により達成され実現される。
　本明細書の記述は典型的な例示に過ぎず、本開示の請求の範囲又は適用例を如何なる意味に於いても限定するものではない。
　本開示に関連する更なる特徴は、本明細書の記述、添付図面から明らかになるものである。また、本開示の態様は、要素及び多様な要素の組み合わせ及び以降の詳細な記述と添付される請求の範囲の様態により達成され実現される。
　本明細書の記述は典型的な例示に過ぎず、本開示の請求の範囲又は適用例を如何なる意味に於いても限定するものではない。
　本開示によれば，複数の発光点からの発光をそれぞれ検出する任意の光学系において発生する，複数の発光点間の空間クロストークを計算処理によって解消または低減することによって，複数の発光点からの発光を識別し，それぞれを独立に検出することが可能となる。また，複数の発光点から発光される，複数種類の蛍光体の蛍光を，それぞれ検出する任意の光学系において発生する，複数の発光点間の空間クロストーク，および各発光点についての複数種類の蛍光体間のスペクトルクロストークを計算処理によって解消または低減することによって，複数の発光点からの，複数種類の蛍光体の蛍光発光を識別し，それぞれを独立に検出することが可能となる。
　さらに，空間クロストークおよびスペクトルクロストークを解消または低減することによって，空間クロストークおよびスペクトルクロストークに伴って発生する検出感度の低下，あるいは検出ダイナミックレンジの低下を回避することが可能となる。
　上記した以外の，課題，構成及び効果は，以下の実施形態の説明により明らかにされる。
単純な光学系の模式図 単純な光学系で2次元センサで取得した発光画像 発光画像における発光像の信号強度分布 発光像の中心の絶対信号強度と，発光像の中心から左側に離れた位置の絶対信号強度の関係 発光像の中心の絶対信号強度と，発光像の中心から右側に離れた位置の絶対信号強度の関係 発光像の中心の絶対信号強度と，発光像の中心から左側に離れた位置の絶対信号強度の関係 発光像の中心の絶対信号強度と，発光像の中心から右側に離れた位置の絶対信号強度の関係 モデル実験系の模式図 モデル実験系における，スペクトルクロストークと空間クロストークの模式図 モデル実験系の模式図 既知サンプルの電気泳動分析で得られた生データ 図11の生データに対して色変換を実施した色変換データ 図11の生データに対して色変換+空間補正を実施した色変換+空間補正データ 未知サンプルの電気泳動分析で得られた生データ 図14の生データに対して色変換を実施した色変換データ 図14の生データに対して色変換+空間補正を実施した色変換+空間補正データ 未知の蛍光体を含むサンプルの電気泳動分析で得られた生データ 図17の生データに対して色変換+空間補正を実施した色変換+空間補正データ 既知の蛍光体の濃度を段階的に上昇させたサンプルの電気泳動分析で得られた生データ 図19の生データに対して色変換を実施した色変換データ 図19の生データに対して色変換+空間補正を実施した色変換+空間補正データ 複数のキャピラリにタイミングをずらしてサンプル注入する方法の模式図 複数のキャピラリでタイミングをずらして蛍光発光を得る方法の模式図 マルチキャピラリ電気泳動装置の模式図 特許文献1の光学系による，4本のキャピラリのラマン散乱光の波長分散像 図25における，各キャピラリの20個の波長帯の検出領域の設定を示す図 特許文献2の光学系による，4本のキャピラリからの発光蛍光を個別に4個の波長帯に分割して結像する模式図 平面上に配列する5個の発光点からの発光蛍光を一括して4個の波長帯に分割して結像する光学系と模式図 従来法のフローチャート 本法のフローチャート 本法の色変換+空間補正を行う行列を決定するためのフローチャート 分析のセッションを複数回繰り返す場合のフローチャート 計算機の構成図
［実施例1］
　空間クロストークの特性を詳細に調べるため，図1に示す単純な光学系を構築した。図1の光学系は，ピンホール板1-1，発光点側開口板1-2，集光レンズ1-3，センサ側開口板1-4，色ガラスフィルタ1-5，2次元センサ1-6，ハロゲンランプ光1-7を照射するハロゲンランプ（光源）を備える。具体的には、以下のように図1の光学系を構成した。φ0.05 mmのピンホールを有するピンホール板1-1に，下方向からハロゲンランプ光1-7を照射することによって，φ0.05 mmの発光点1-8を形成した。また，φ0.2 mmの開口を有する発光点側開口板1-2を発光点1-8の上方向に0.2 mm離れた位置に配置し，焦点距離がf＝1.4 mmの集光レンズ1-3を発光点1-8の上方向に1.54 mm離れた位置に配置し，φ0.7 mmの開口を有するセンサ側開口板1-4を集光レンズ1-3の直上に配置した。さらに，2次元センサ1-6を集光レンズ1-3の上方向に15 mm離れた位置に配置し，色ガラスフィルタ1-5を2次元センサ1-6の直下に配置した。以上で，ピンホール板1-1，発光点側開口板1-2，センサ側開口板1-4，色ガラスフィルタ1-5，2次元センサ1-6は互いに平行に配置した。発光点1-8から発光した光1-9は，φ0.2 mmの開口を透過し，集光レンズ1-3で集光され，φ0.7 mmの開口を透過し，色ガラスフィルタ1-5を透過し，2次元センサ1-6上に，φ0.5 mmの発光像1-10を形成した。ここで，発光点1-8は，2次元センサ1-6上にピントが合った状態で，像倍率10倍で結像された。
　図2は，図1の単純な光学系において2次元センサ1-6で得られた像1-10を含む発光画像である。2次元センサ1-6のセンササイズは13×13 mmであり，各画素の信号範囲は0～65536である。図2(a)と図2(b)は同じ発光画像であるが，図2(a)の信号表示スケール（階調）を0～50000とし，図2(b)の信号表示スケールを0～500としている。発光像の最大信号値は約50000であり，図2(a)では全画素の信号がほぼ飽和せずに表示されているのに対して，図2(b)では発光像が飽和して表示されている。図2(a)によれば，期待通りにφ0.5 mmの発光像が得られており，φ0.5 mmの発光像の外側は信号ゼロに見える。しかしながら，図2(b)によれば，φ0.5 mmの発光像の外側に低い強度の裾が広がっていることが確認できる。
　図3(a)～(c)のα（点線）は，図2の発光画像において，発光像の中心を通る，水平方向の信号強度分布を示している。図3(a)，図3(b)，および図3(c)は同じデータを用いているが，縦軸を変更している。横軸は共通であり，発光像の中心からの距離を示している。横軸のプラスおよびマイナスは，発光画像において，発光像の右側と左側をそれぞれ示している。図3(a)と(b)の縦軸は絶対信号強度を示し，図3(a)は縦軸スケールを0～60000とし，図3(b)は縦軸スケールを0～100としている。一方，図3(c)は，縦軸をログスケールとして，かつ最大信号強度を100％とする相対信号強度として，縦軸スケールを0.001％～100％としている。図3(a)～(c)のβ（実線），γ（破線）は，図2の発光画像の取得条件下で，ハロゲンランプ光1-7の出力強度を段階的に下げた場合について，αと同様に信号強度分布を示したものである。図3(a)から分かるように，α，β，γの発光像の最大信号強度はそれぞれ約50000，約25000，約10000である。一方，図3(b)を見ると，発光像の中心から離れるに従って信号強度が減少する傾向があるが，図3(a)で見られる発光像の大きさ（幅0.5 mm程度）と比較して，ずっと大きな裾が広がっていることが分かる。また，最大信号強度が減少すると，裾の強度も減少していることも分かる。そこで，図3(c)を見ると，α，β，γがひとつに重なることが分かる。これは新たな発見であり，いくつかの重要な知見をもたらす。例えば，発光像の中心から±1 mm離れた位置，つまり，図2(a)で見られる発光像の外側の位置において，最大信号強度の約0.1％の信号強度が観察されている。仮に，±1 mm離れた位置に隣接する発光像の検出領域が存在した場合，これは約0.1％の空間クロストークが存在することを意味する。図1に示すように，極めて単純化された光学系を用いて，発光点をジャストフォーカスしているにも関わらず，このように無視できない大きさの空間クロストークが存在することは驚きである。さらに驚くべきは，α，β，γがひとつに重なったことから，光学系および発光点の位置が固定されている限り，発光点の発光強度，あるいは結像点の信号強度によらずに，結像点の相対信号強度分布が一定であると判明したことである。
　図4(a)および図5(a)は，図3から導いた結果である。図4(a)は結像点の中心から-1，-2，-3，-4，-5，-6 mm離れた位置，図5(a)は結像点の中心から+1，+2，+3，+4，+5，+6 mm離れた位置における絶対信号強度を縦軸，いずれの場合も発光像の中心の絶対信号強度を横軸とするグラフである。上記の絶対信号強度は，図3における，各位置を中心とする±0.1 mm幅の絶対信号強度の平均値とした。その結果，各位置における3個のプロット（図3のα，β，γに対応）はいずれも，原点を通る近似直線に乗ることが判明した。各近似直線の傾きは，対応する位置における空間クロストーク比率を示す。つまり，結像点の中心から離れるに従い，近似直線の傾きが小さくなり，空間クロストーク比率が低下することが分かる。比較のため，空間クロストークが0.1％および0.01％の直線を点線で重ね表示してある。例えば，結像点の中心から-1 mm離れた位置の空間クロストーク比率は0.1％強であるのに対して，-3 mm離れた位置の空間クロストーク比率は0.01％程度であることが読み取れる。以上の結果は，結像点の中心から離れた任意の位置における空間クロストークの絶対信号強度は，結像点の中心の絶対信号強度に対して直線性があり，空間クロストーク比率が一定になることを示している。すなわち，結像点の中心から離れた任意の位置における絶対信号強度から，結像点の中心の絶対信号強度に空間クロストーク比率を乗じた値を差引くことによって，任意の位置における空間クロストークを解消または減少させることが可能であることが新たに見出された。これは，図4(a)および図5(a)における各プロットの絶対信号強度から，同じ位置における対応する近似直線の値を差引くことに相当する。
　図4(b)および図5(b)はそれぞれ，図4(a)および図5(a)の結果に対して，上記の差引演算を行った結果である。各位置を示すプロットの種類，および横軸は，図4(a)と図4(b)，図5(a)と図5(b)で共通である。図4(b)および図5(b)の縦軸は，上記の差引演算を行った後の絶対信号強度を示し，図4(a)および図5(a)と比較してスケールを拡大してある。比較のため，空間クロストークが±0.01％の直線を点線で重ね表示してある。その結果，各位置における差引演算後の絶対信号強度が，像中心の絶対信号強度によらずに，ほぼゼロになり，各位置における空間クロストークが±0.01％よりも小さくなることが判明した。結像点の中心から-1 mm離れた位置で比較すると，本法によって，空間クロストークを少なくとも1桁以上低減させることができた。各位置における差引演算後の絶対信号強度が完全にゼロにならないのは，図4(a)および図5(a)における各プロットと対応する近似直線のずれ，誤差が存在するためである。各位置における，原点を通る近似直線の決定係数（R2）が高いほど（1に近いほど），つまり，結像点の中心の絶対信号強度に対する，各位置の絶対信号強度の直線性が高いほど，上記の誤差は小さくなり，差引演算後の絶対信号強度はゼロに近くなる。当然ながら，差引演算後の絶対信号強度は負になることもある。しかしながら，上記の直線性が存在する限り，各位置の絶対信号強度の大きさは差引演算によって少なくとも小さくなる。また，近似直線の傾きが小さいほど，つまり，差引演算前の空間クロストークが小さいほど，各位置の絶対信号強度と近似直線の差分の絶対値も小さくなるため，上記の誤差も小さくなる。
　図6および図7は，図4および図5と同様の実験結果である。図1の光学系において，ピンホール板1-1を一旦取り外し，再度取り付けた後，図2と図3と同様の結果を取得し，同様の手法により図6および図7の結果を導出した。ピンホール板1-1の取り付け位置は，ほぼ元の位置と同じであるが，厳密に全く同じ位置ではない。図6および図7の結果は，図4および図5の結果と同等であり，差引演算処理によって，任意の位置における空間クロストークを解消または減少させる手法の再現性が高いことが示された。
　一方で，図4～図7の比較から，次に示すように，もうひとつ重要な知見が得られた。図4(a)の-1 mm，図5(a)の+1 mm，図6(a)の-1 mm，および図7(a)の+1 mmの近似直線を比較する。これら4個の近似直線の傾きは，いずれも結像点の中心から1 mm離れた位置における空間クロストーク比率を示すが，相互に異なっていることが判明した。したがって，例えば，図4(a)の-1 mmの近似直線を用いて，図5(a)の+1 mm，図6(a)の-1 mm，または図7(a)の+1 mmの各プロットの差引演算処理を行うと，上記の誤差が大きくなり，空間クロストークの低減が不十分になることが明らかになった。この現象は，他の位置についても同様に生じており，場合によっては，差引演算処理が逆効果，つまり，空間クロストークを増大させてしまうこともあることが分かった。このような現象が発生する要因は，（1）図1のような単純な光学系においてすら，結像が及ぼす空間クロストークの点対称性が必ずしも高くないこと（例えば図2で，結像点の左方向と右方向で空間クロストーク比率が異なる），（2）発光点の位置がわずかにずれる等，光学系の微妙な変化に伴って，結像が及ぼす空間クロストークが有意に変化すること，等が考えられる。したがって，結像点の中心と任意の位置の間に発生する空間クロストーク比率を，両者の間の距離の関数で表すことは不可能であることが判明した。これは，特許文献3に示される手法が有効に機能しないことを意味している。本開示では，光学系，および複数の結像点，すなわち複数の検出領域の条件を固定し，相互の距離とは無関係に，相互の空間クロストーク比率を，図4～図7のそれぞれで示すように，実験で導出する点で特許文献3と異なる。条件に変更が加われば，その都度，空間クロストーク比率を取得し直す必要がある。
　以上の実験で得られた知見，すなわち，結像点の中心から離れた任意の位置における空間クロストークの絶対信号強度は，結像点の中心の絶対信号強度に対して直線性があり，空間クロストーク比率が一定になる事実を拡張し，一般化すると（数9）～（数17）が導出される。
［実施例2］
　複数の発光点から発光される，複数種類の蛍光体の蛍光を，それぞれ複数の波長帯の検出領域で検出する光学系の最もシンプルな場合についてモデル実験を行った。本実施例，および同じモデル実験系を用いる他の実施例においては，例として，（数9）～（数12）において，A＝2個の発光点それぞれについて，B＝2個の波長帯の検出領域で，C=2種類の蛍光体の発光蛍光を検出する場合を選んだが，A，B，Cのそれぞれが他の数値である場合についても同様の方法によって同様の効果が得られることは言うまでもない。
　図8はモデル実験系を示す。図8のモデル実験系は，2本のキャピラリCap(1)およびCap(2)と、2本のキャピラリCap(1)およびCap(2)の配列方向にレーザビームLBを照射する光源（不図示）とを備える。2本のキャピラリCap(1)およびCap(2)上に，レーザビームLBが照射される位置に，発光点P(1)およびP(2)がそれぞれ設けられている。キャピラリCap(1)で2種類の蛍光体D(1,1)およびD(1,2)で標識された物質が電気泳動され，キャピラリCap(2)で2種類の蛍光体D(2,1)およびD(2,2)で標識された物質が電気泳動されている。ここで，D(1,1)とD(2,1)は同じ種類の蛍光体であり，D(1,2)とD(2,2)は同じ種類の蛍光体である。同じ種類の蛍光体とは，発光される蛍光の蛍光スペクトルが等しいことを意味する。蛍光体D(1,1)およびD(1,2)は，発光点P(1)において，レーザビームLBの照射によって蛍光を発光し，センサSに設けられた検出領域W(1,1)およびW(1,2)で検出される。検出領域W(1,1)は蛍光体D(1,1)の発光蛍光を，検出領域W(1,2)は蛍光体D(1,2)の発光蛍光を主として検出するように，それぞれの検出波長帯が設計されているが，相互のスペクトルクロストークが有意に存在している。同様に，蛍光体D(2,1)およびD(2,2)は，発光点P(2)において，レーザビームLBの照射によって蛍光を発光し，センサSに設けられた検出領域W(2,1)およびW(2,2)で検出される。検出領域W(2,1)は蛍光体D(2,1)の発光蛍光を，検出領域W(2,2)は蛍光体D(2,2)の発光蛍光を主として検出するように，それぞれの検出波長帯が設計されているが，相互のスペクトルクロストークが有意に存在している。また，発光点P(1)における蛍光体D(1,1)およびD(1,2)の蛍光発光は検出領域W(2,1)およびW(2,2)においても検出され，同様に，発光点P(2)における蛍光体D(2,1)およびD(2,2)の蛍光発光は検出領域W(1,1)およびW(1,2)においても検出され，空間クロストークが有意に存在している。検出領域W(1,1)，W(1,2)，W(2,1)，およびW(2,2)における信号強度（蛍光強度）をX(1,1)，X(1,2)，X(2,1)，およびX(2,2)とすると，図8下側に示されているように，発光点P(1)，発光点P(2)のそれぞれについて，電気泳動に伴い，信号強度X(1,1)およびX(1,2)，X(2,1)およびX(2,2)の時系列データが得られる。
　尚，図8では，検出領域W(1,1)，W(1,2)，W(2,1)，およびW(2,2)が1個のセンサS上に設けられているように描かれているが，必ずしもその必要はない。検出領域W(1,1)とW(1,2)がひとつのセンサ上に，検出領域W(2,1)とW(2,2)がもうひとつのセンサ上に設けられていても良いし，検出領域W(1,1)，W(1,2)，W(2,1)，およびW(2,2)が異なる4個のセンサ上にそれぞれ設けられていても良い。また，発光点P(1)およびP(2)と，センサSの間に，何らかの結像手段および分光手段が必要であるが，本実施例ではどのような結像手段および分光手段であるかを問わないため，図8では結像手段および分光手段を省略している。
　図9は，図8におけるスペクトルクロストークおよび空間クロストークの関係を示す模式図である。図9(a)では，発光点P(1)における蛍光体D(1,1)の蛍光発光が，(i)検出領域W(1,1)で主として検出され，(ii)検出領域W(1,2)で副として検出されている。また，上記発光蛍光は，(iii)検出領域W(2,1)でも(i)および(ii)と比較してずっと小さい強度で検出され，さらに，(iv)検出領域W(2,2)でも(iii)よりさらに小さい強度で検出される。(i)と(ii)はスペクトルクロストークを示しており，(iii)と(iv)は空間クロストークおよびスペクトルクロストークを示している。つまり，従来法では(i)と(ii)だけを考慮していたが，本開示では(i)～(iv)のすべてを考慮する。図9では，以上の検出強度の強弱の関係，およびクロストークの関係を，矢印の太さ，実線または点線，(i)～(iv)の番号で共通的に示している。図9(a)と同様に，図9(b)は発光点P(1)における蛍光体D(1,2)の蛍光発光，図9(c)は発光点P(2)における蛍光体D(2,1)の蛍光発光，図9(d)は発光点P(2)における蛍光体D(2,2)の蛍光発光の場合についてそれぞれ描いている。実際の分析においては，図9(a)～(d)が同時刻に一斉に発生し，しかもそれぞれの蛍光強度が異なっている。
　図10は，図8と同様なモデル実験系を示す。ただし，図10のモデル実験系では，2本のキャピラリCap(1)およびCap(2)に対し，その配列方向と垂直に光源からランプ光LLが照射される。したがって，レーザビームLBの照射に伴う発光蛍光が検出されるのではなく，ランプ光LLの透過に伴う吸光，あるいは吸光度が検出される。2本のキャピラリCap(1)およびCap(2)上に，ランプ光LLが照射される位置に，吸光点P(1)およびP(2)がそれぞれ設けられている。キャピラリCap(1)で2種類の吸光体D(1,1)およびD(1,2)が電気泳動され，キャピラリCap(2)で2種類の吸光体D(2,1)およびD(2,2)が電気泳動されている。ここで，D(1,1)とD(2,1)は同じ種類の吸光体であり，D(1,2)とD(2,2)は同じ種類の吸光体である。同じ種類の吸光体とは，光の吸光スペクトルが等しいことを意味する。吸光体D(1,1)およびD(1,2)は，吸光点P(1)において，ランプ光LLの照射によって光を吸収し，吸収されなかった光がセンサSに設けられた検出領域W(1,1)およびW(1,2)で検出される。検出領域W(1,1)は吸光体D(1,1)の吸光を，検出領域W(1,2)は吸光体D(1,2)の吸光を主として検出するように，それぞれの検出波長帯が設計されているが，相互のスペクトルクロストークが有意に存在している。同様に，吸光体D(2,1)およびD(2,2)は，吸光点P(2)において，ランプ光LLの照射によって光を吸収し，吸収されなかった光がセンサSに設けられた検出領域W(2,1)およびW(2,2)で検出される。検出領域W(2,1)は吸光体D(2,1)の吸光を，検出領域W(2,2)は吸光体D(2,2)の吸光を主として検出するように，それぞれの検出波長帯が設計されているが，相互のスペクトルクロストークが有意に存在している。また，吸光点P(1)における吸光体D(1,1)およびD(1,2)の吸光は検出領域W(2,1)およびW(2,2)においても検出され，同様に，吸光点P(2)における吸光体D(2,1)およびD(2,2)の吸光は検出領域W(1,1)およびW(1,2)においても検出され，空間クロストークが有意に存在している。これは，吸光点P(1)を透過するランプ光LLの透過光の一部が検出領域W(2,1)およびW(2,2)においても検出され，吸光点P(2)を透過するランプ光LLの透過光の一部が検出領域W(1,1)およびW(1,2)においても検出されるためである。検出領域W(1,1)，W(1,2)，W(2,1)，およびW(2,2)における吸光度をX(1,1)，X(1,2)，X(2,1)，およびX(2,2)とすると，図8下側に示されているように，吸光点P(1)，吸光点P(2)のそれぞれについて，電気泳動に伴い，吸光度X(1,1)およびX(1,2)，X(2,1)およびX(2,2)の時系列データが得られる。以降では，図8を用いた場合についての説明を行うが，図10を用いた場合にも同様の効果が得られることは言うまでもない。
　図11は，図8のモデル実験系を用いて，既知のサンプルをキャピラリCap(1)およびCap(2)に注入して電気泳動分析で得られた生データである。いずれも，横軸を電気泳動時間（任意単位），縦軸を蛍光強度（任意単位）としている。図11(a)と図11(b)は発光点P(1)で得られた同一の生データであり，(b)は(a)の縦軸スケールを拡大したものである。同様に，図11(c)と図11(d)は発光点P(2)で得られた同一の生データであり，図11(d)は図11(c)の縦軸スケールを拡大したものである。図11(a)～(d)の横軸スケールは時刻0～500を示し，図11(a)(c)の縦軸スケールは蛍光強度0～250，図11(b)(d)の縦軸スケールは蛍光強度-0.2～0.8としている。X(1,1)およびX(2,1)は実線，X(1,2)およびX(2,2)は点線で表されている。図11に示されている通り，時刻100において発光点P(1)で蛍光体D(1,1)（図9(a)の状態），時刻200において発光点P(2)で蛍光体D(2,1)（図9(c)の状態），時刻300において発光点P(1)で蛍光体D(1,2)（図9(b)の状態），時刻400において発光点P(2)で蛍光体D(2,2)（図9(d)の状態）がそれぞれ単独で蛍光発光するようにサンプルが調製されている。そして，これら以外の蛍光は発光されないようにサンプルが調製されている。図11(a)と図11(c)を見る限り，上記の蛍光発光に対応する4個の大きなピークだけが観察されているが（図9の(i)(ii)に対応），図11(b)と図11(d)を見ると，矢印で示す4個の小さなピーク，すなわち，時刻100において発光点P(2)で蛍光体D(2,1)，時刻200において発光点P(1)で蛍光体D(1,1)，時刻300において発光点P(2)で蛍光体D(2,2)，時刻400において発光点P(1)で蛍光体D(1,2)が観察された（図9の(iii)(iv)に対応）。したがって，これら4個の小さなピークは，空間クロストークおよびスペクトルクロストークによるものであると判断できる。例えば，時刻100の発光点P(2)における蛍光体D(2,1)の小さなピークは，時刻100の発光点P(1)における蛍光体D(1,1)の蛍光発光の一部が空間クロストークおよびスペクトルクロストークによってX(2,1)とX(2,2)として検出された結果である（図9(a)の(iii)(iv)に対応）。
　空間クロストークを考慮しない場合，つまり，図11(c)および図11(d)の4個の小さなピークを無視する場合（図9で(i)(ii)のみを考慮する場合），図11(a)および図11(c)の4個の大きなピークのスペクトルクロストークの解消，すなわち従来法による色変換を実行すれば良い。このとき，（数6）は次式で表される。
Figure JPOXMLDOC01-appb-M000024
Figure JPOXMLDOC01-appb-M000025
　ここで，Z(1,1)，Z(1,2)，Z(2,1)，およびZ(2,2)は，各時刻における蛍光体D(1,1)，D(1,2)，D(2,1)，およびD(2,2)の濃度を示す。（数18）および（数19）に示す通り，本実施例では，2行2列の行列Yおよび逆行列Y-1は，発光点P(1)と発光点P(2)で同じであった。ただし，一般には，発光点P(1)と発光点P(2)で行列Yおよび逆行列Y-1が異なっていることがある。発光点P(1)についての行列Yの各要素は，時刻100の蛍光体D(1,1)による蛍光強度X(1,1)とX(1,2)の強度比，時刻300の蛍光体D(1,2)による蛍光強度X(1,1)とX(1,2)の強度比により決定した。発光点P(2)についての行列Yの各要素は，時刻200の蛍光体D(2,1)による蛍光強度X(2,1)とX(2,2)の強度比，時刻400の蛍光体D(2,2)による蛍光強度X(2,1)とX(2,2)の強度比により決定した。
　図12は，図11の生データに対して，各時刻で（数18）および（数19）の色変換を実行して得られた色変換データである。表記方法は図11と同等である。期待通りに，図11(a)および図11(c)の4個の大きなピークのスペクトルクロストークが，図12(a)および図12(c)の4個の大きなピークで解消されていることが分かる。一方で，図11(b)および図11(d)の4個の小さなピークのスペクトルクロストークも，図12(c)および図12(d)の4個の小さなピーク（同様に矢印で示す）で解消されているが，これらのピークそのものは依然として残っており，空間クロストークは解消されていないことが分かる。これが従来法の課題である。
　そこで次に，図11について，スペクトルクロストークと空間クロストークの両方を考慮し（図9で(i)～(iv)をすべて考慮し），色変換と空間補正の両方を実行する。このとき，（数14）は次式で表される。
Figure JPOXMLDOC01-appb-M000026
　行列Yおよび逆行列Y-1は，（数18）と（数19）の2行2列から，4行4列に拡張されている。また，（数18）と（数19）では発光点P(1)と発光点P(2)でそれぞれ独立に色変換を行っていたのに対して，（数20）では発光点P(1)と発光点P(2)を混在させて一括で色変換＋空間補正を行っている。（数20）の行列Yの各要素は，図11における，時刻100の蛍光体D(1,1)による蛍光強度X(1,1)，X(1,2)，X(2,1)，およびX(2,2)の強度比，時刻200の蛍光体D(2,1)による蛍光強度X(1,1)，X(1,2)，X(2,1)，およびX(2,2)の強度比，時刻300の蛍光体D(1,2)による蛍光強度X(1,1)，X(1,2)，X(2,1)，およびX(2,2)の強度比，および時刻400の蛍光体D(2,2)による蛍光強度X(1,1)，X(1,2)，X(2,1)，およびX(2,2)の強度比により決定した。（数20）の行列Yの左上の2行2列，および右下の2行2列は，（数18）および（数19）の行列Yと対応しており，各要素の値が互いにほぼ等しい。同様に，（数20）の行列Y-1の左上の2行2列，および右下の2行2列は，（数18）および（数19）の行列Y-1と対応しており，各要素の値が互いにほぼ等しい。つまり，これらの要素は，同一発光点についてのスペクトルクロストークを解消する色変換を担っている。ここで，これらの要素の値が（数20）と（数18）および（数19）で完全に一致しない理由は，次に述べる空間クロストークの考慮の有無の違いである。これに対して，（数20）の行列Yおよび行列Y-1の右上の2行2列，および左下の2行2列は，（数18）および（数19）には存在しない，異なる発光点間の空間クロストークおよびスペクトルクロストークを解消する空間補正および色変換を担っている。（数20）を次式のように変形すると，上記がより分かり易くなる。
Figure JPOXMLDOC01-appb-M000027
Figure JPOXMLDOC01-appb-M000028
　（数21）右辺第1項，および（数22）右辺第2項は従来法の色変換を担っており，それぞれ（数18）および（数19）と対応している。これに対して，（数21）右辺第2項，および（数22）右辺第1項は，従来法では取り扱っていない，空間補正および色変換を担っている。このように，色変換と，空間補正および色変換を切り離して個別に実行しても良いし，どちらか一方のみを行っても良い。あるいは，例えば，発光点P(1)には空間補正および色変換を施し，発光点P(2)には色変換のみを施す，等，発光点によって処理を変化させることも可能である。
　図13は，図11の生データに対して，各時刻で（数20）の色変換+空間補正を実行して得られた色変換+空間補正データである。表記方法は図11と同等である。まず，図12と同様に，図11(a)および図11(c)の4個の大きなピークのスペクトルクロストークが，図13(a)および図13(c)の4個の大きなピークで解消されていることが分かる。さらに，期待通りに，図11(b)および図11(d)の4個の小さなピークの空間クロストークおよびスペクトルクロストークも，図13(c)および図13(d)で解消されて，ピークが消失していること（矢印で示す）が分かる。
　以上により，本開示によって，複数の発光点から発光される，複数種類の蛍光体の蛍光を，それぞれ検出する任意の光学系において発生する，同一発光点についてのスペクトルクロストーク，および複数の発光点間の空間クロストーク，およびスペクトルクロストークを計算処理によって解消または低減することが可能であることが示された。
［実施例3］
　次に，［実施例2］の実験結果を踏まえて，今度は未知のサンプルの分析を行った。図14は，図8のモデル実験系を用いて，未知のサンプルをキャピラリCap(1)およびCap(2)に注入して電気泳動分析を行った結果得られた生データである。［実施例2］と同一のモデル実験系を用いているため，（数18）～（数22）をそのまま用いることができる。表記方法は図11と同等である。ただし，図11(b)(d)の縦軸スケールを蛍光強度-0.5～2.0に縮小している。時刻100，200，300，および400それぞれにおいて，発光点P(1)で大きなピークが観察される一方で，発光点P(2)で小さなピーク（矢印で示す）が観察された。
　図15は，図14の生データに対して，各時刻で（数18）および（数19）の色変換を実行して得られた色変換データである。表記方法は図14と同等である。図15(a)に示すように，発光点P(1)において，時刻100および300で蛍光体D(1,1)の発光蛍光が単独で検出され，時刻200および400で蛍光体D(1,2)の発光蛍光が単独で検出された。一方で，図15(d)に示すように，発光点P(2)において，時刻100，200，300，および400で検出された4個の小さなピーク（矢印で示す）の正体は不明であった。すなわち，それらの由来が，（1）発光点P(2)における蛍光体D(2,1)とD(2,2)の混合物の蛍光発光なのか，（2）発光点P(2)における蛍光体D(2,1)とD(2,2)以外の不純物の蛍光発光なのか，あるいは（3）発光点P(1)における蛍光体D(1,1)およびD(1,2)の蛍光発光の空間クロストークなのか，が不明であった。
　図16は，図14の生データに対して，各時刻で（数20）の色変換+空間補正を実行して得られた色変換+空間補正データである。表記方法は図14と同等である。図16(a)に示すように，図15(a)と同様に，発光点P(1)において，時刻100および300で蛍光体D(1,1)の発光蛍光が単独で検出され，時刻200および400で蛍光体D(1,2)の発光蛍光が単独で検出された。それに加えて，図16(d)に示すように，発光点P(1)における蛍光体D(1,1)およびD(1,2)の蛍光発光の検出領域W(2,1)およびW(2,2)に対する空間クロストークおよびスペクトルクロストークが解消された結果，発光点P(2)において，時刻100および300で蛍光体D(2,2)の微弱な発光蛍光が単独で検出され，時刻200および400で蛍光体D(2,1)の微弱な発光蛍光が単独で検出されていることが判明した。これらの微弱な発光蛍光のピーク強度は，発光点P(1)で観察された発光蛍光のピーク強度の1％弱であった。一方，図12の結果から分かるように，図8のモデル実験系で発生する空間クロストーク比率も1％弱である。そのため，図15では，真の微弱な発光蛍光と，空間クロストークによる偽の微弱な発光蛍光とが混在し，両者の区別がつかなくなったのであった。
　上記の結果は，空間クロストークが各発光点からの発光の検出における検出下限を押し上げ得ることを示している。上記の例のように，空間クロストーク比率が1％だとすると，仮に発光点が1個の場合の検出下限が0.1％であったとしても，1％以下の信号は，それが真の信号なのか，あるいは空間クロストークによる偽の信号なのかを区別できないため，実効的な検出下限は1％に上昇してしまう。つまり，1個の発光点の発光検出と比較して，複数の発光点の発光検出では，検出感度とダイナミックレンジがいずれも1桁低くなってしまう。本開示は，そのような課題を解決し，複数の発光点の発光検出における検出感度とダイナミックレンジの低下を回避することができる。
［実施例4］
　続いて，［実施例2］の実験結果を踏まえて，発光点P(1)で，蛍光体D(1,1)と蛍光体D(1,2)と一緒に，未知の蛍光体D(1,3)で標識された物質を電気泳動分析した。図17は，図8のモデル実験系を用いて，キャピラリCap(1)にのみサンプルを注入し，キャピラリCap(2)にはサンプルを注入せずに，キャピラリCap(1)とCap(2)の両方について電気泳動分析して得られた生データである。表記方法は図14と同等である。時刻100において発光点P(1)で蛍光体D(1,1)，時刻200において発光点P(1)で蛍光体D(1,2)，時刻300において発光点P(1)で蛍光体D(1,1)，時刻400において発光点P(1)で蛍光体D(1,3)がそれぞれ単独で蛍光発光するようにサンプルが調製されている。そして，これら以外の蛍光は発光されないようにサンプルが調製されている。図17(a)に示されている通り，上記の蛍光発光に対応する4個の大きなピークだけが観察された。蛍光体D(1,3)のスペクトルクロストーク比率は，蛍光体D(1,2)のスペクトルクロストーク比率と似ているが，若干異なっていることが分かった。一方で，図17(d)に示される，矢印で示す4個の小さなピークはそれぞれ，図17(a)の4個の大きなピークの元となる発光蛍光の空間クロストークおよびスペクトルクロストークの結果であると判断できる。
　図18は，図17の生データに対して，各時刻で（数20）の色変換+空間補正を実行して得られた色変換+空間補正データである。表記方法は図14と同等である。図18(a)に示すように，発光点P(1)において，時刻100および300で蛍光体D(1,1)の発光蛍光が単独で検出され，時刻200で蛍光体D(1,2)の発光蛍光が単独で検出された。しかしながら，矢印で示す通り，時刻400で蛍光体D(1,3)の発光蛍光が単独で検出されたが，スペクトルクロストーク比率が蛍光体D(1,1)と蛍光体D(1,2)のいずれとも異なるため，スペクトルクロストークの解消はなされなかった。一方で，図18(d)に示すように，時刻100，200，および300における，発光点P(1)での蛍光体D(1,1)およびD(1,2)の蛍光発光による空間クロストークおよびスペクトルクロストークが解消されているにも関わらず，矢印で示す通り，時刻400における，発光点P(1)での蛍光体D(1,3)の蛍光発光による空間クロストークおよびスペクトルクロストークは解消されずに，ピークが残っていることが判明した。これは，［実施例2］で，発光点P(1)での蛍光体D(1,1)，D(1,2)，および発光点P(2)での蛍光体D(2,1)，D(2,2)の発光蛍光について（数20）の導出を行ったためであり，発光点P(1)でのD(1,3)の発光蛍光は，空間クロストークおよびスペクトルクロストークの特性が（数20）と異なるためである。したがって，発光点P(1)での蛍光体D(1,3)の発光蛍光について，（数20）を取得し直せば，上記の問題は解決される。以上の現象は本開示の一側面を表している。
［実施例5］
　今度は，［実施例2］の実験結果を踏まえて，発光点P(1)で，蛍光体D(1,1)で標識された物質の濃度を段階的に引き上げて蛍光発光させる実験を行った。図19は，図8のモデル実験系を用いて，キャピラリCap(1)にのみサンプルを注入し，キャピラリCap(2)にはサンプルを注入せずに，キャピラリCap(1)とCap(2)の両方について電気泳動分析して得られた生データである。表記方法は図14と同等である。ただし，図19(a)(c)の縦軸スケールを蛍光強度0～400に縮小している。時刻100，200，300，400のそれぞれにおいて発光点P(1)で蛍光体D(1,1)がそれぞれ単独で蛍光発光し，蛍光体D(1,1)の濃度および蛍光強度が段階的に上昇するようにサンプルが調製されている。そして，これら以外の蛍光は発光されないようにサンプルが調製されている。図19(a)に示されている通り，上記の蛍光発光に対応する4個の大きなピークが観察されたが，図8のモデル実験系で用いているセンサSは蛍光強度200で飽和するため，時刻200，300，400の3個のピークが飽和して検出された。これに対して，図19(d)に示される，発光点P(1)における蛍光体D(1,1)の発光蛍光の空間クロストークおよびスペクトルクロストークによって生じている4個の小さなピークは，蛍光強度が低いため，いずれも飽和せずに観察された。
　図20は，図19の生データに対して，各時刻で（数18）および（数19）の色変換を実行して得られた色変換データである。表記方法は図19と同等である。図20(a)に示すように，発光点P(1)における，時刻100の蛍光体D(1,1)の発光蛍光は，スペクトルクロストークが解消され，単独で検出された。しかしながら，時刻200，300，400の蛍光体D(1,1)の発光蛍光は飽和しているため，スペクトルクロストークが解消されなかった。これは，図19(a)で明らかなように，蛍光強度が飽和すると，X(1,1)とX(1,2)の強度比から導出されるスペクトルクロストーク比率が変化し，（数18）および（数19）で定義されているスペクトルクロストーク比率から乖離したためである。これに対して，図20(d)に示されるように，蛍光体D(2,1)の蛍光強度が飽和していないため，矢印で示す4個の小さなピークそれぞれについてスペクトルクロストークが解消された。図20(d)に示される4個のピークの強度比が，時刻100，200，300，400における発光点P(1)での蛍光体D(1,1)の発光蛍光の強度比を表している。
　図21は，図19の生データに対して，各時刻で（数20）の色変換+空間補正を実行して得られた色変換+空間補正データである。表記方法は図19と同等である。図21(a)では，図20(a)と同様な理由によって同等な結果が得られている。これに対して，図21(d)に示すように，時刻100における蛍光体D(2,1)の発光蛍光については，空間クロストークとスペクトルクロストークの両方が解消された。しかしながら，時刻200，300，400における蛍光体D(2,1)の発光蛍光については，空間クロストークとスペクトルクロストークが解消されなかった。これは，時刻200，300，400において，図19(a)に示されているX(1,1)またはX(1,2)の蛍光強度が飽和した結果，それらのX(1,1)およびX(1,2)の蛍光強度に対する，図19(d)に示されているX(2,1)およびX(2,2)の蛍光強度の比率が，（数20）で定義されている空間クロストーク比率から乖離したためである。また，その乖離の度合いは，飽和の度合いに応じて増大するため，図21(d)において，時刻200よりも時刻300，時刻300よりも時刻400の蛍光体D(2,1)の蛍光強度が増大した。
［実施例6］
　ここでは，［実施例2］と異なる手法により，（数20）の行列Yを決定する手法を説明する。［実施例2］では，図11(a)のように，蛍光体D(1,1)，D(2,1)，D(1,2)，およびD(2,2)の順番で，それぞれを単独で蛍光発光させることによって，つまり，発光点P(1)と発光点P(2)で交互に蛍光発光させることによって，行列Yを決定した。しかしながら，キャピラリCap(1)とCap(2)のどちらかの電気泳動速度が想定からずれると，発光点P(1)と発光点P(2)で同時に蛍光発光する可能性があり，その場合は行列Yを決定できなくなる。
　図22は，上記の問題を回避するための，より現実的な方法を示している。最初に，図22(a)のように，キャピラリCap(1)にのみサンプルを注入し，キャピラリCap(1)とCap(2)で電気泳動分析を行い，蛍光体D(1,1)とD(1,2)をそれぞれ単独で蛍光発光させて，検出領域W(1,1)，W(1,2)，W(2,1)，W(2,2)で蛍光検出を行う。次に，図22(b)のように，キャピラリCap(2)にのみサンプルを注入し，キャピラリCap(1)とCap(2)で電気泳動分析を行い，蛍光体D(2,1)とD(2,2)をそれぞれ単独で蛍光発光させて，検出領域W(1,1)，W(1,2)，W(2,1)，W(2,2)で蛍光検出を行う。このような方法によって，［実施例2］と同様に，行列Yを決定することができる。また，この方法によれば，［実施例2］よりも簡便かつ確実に行列Yを決定することができる。図22では，分かり易くするために，キャピラリCap(1)における蛍光体D(1,1)とD(1,2)，およびキャピラリCap(2)における蛍光体D(2,1)とD(2,2)をそれぞれ複数回，蛍光発光させるように描いているが，それぞれ1回ずつ蛍光発光させれば十分である。また，上記では，最初の電気泳動分析が終わってから，次の電気泳動分析を行っているが，そこまで電気泳動分析の間隔を広げなくても良い。キャピラリCap(1)のサンプル注入とキャピラリCap(2)のサンプル注入のタイミングを適当にずらすことによって，図22と同等の効果を得ることができる。例えば，キャピラリCap(1)にサンプル注入した後に短い時間だけ電気泳動し，次いでキャピラリCap(2)にサンプル注入してから電気泳動を再開すれば良い。これによって，行列Yを決定するのに要する時間を削減することができる。また，この方法によれば，キャピラリCap(1)に注入するサンプルと，キャピラリCap(2)に注入するサンプルを同じ組成，あるいは同一にすることができる。これによって，準備するサンプルを単純化でき，そのためのコストを削減することができる。
　図23は，行列Yの決定をさらに簡便にする方法を示している。ここでは，キャピラリCap(1)のサンプル注入とキャピラリCap(2)のサンプル注入を同じタイミングで行ったにも関わらず，図22と同様に，キャピラリCap(1)における蛍光体D(1,1)とD(1,2)，およびキャピラリCap(2)における蛍光体D(2,1)とD(2,2)がそれぞれ異なる時刻に蛍光発光するようにサンプルが調製されている。つまり，キャピラリCap(1)に注入するサンプルと，キャピラリCap(2)に注入するサンプルの組成を変え，各蛍光体が標識された物質の電気泳動速度に違いが出るようにすれば良い。例えば，キャピラリCap(1)に注入するサンプルには，蛍光体D(1,1)が標識された50塩基長のDNA断片，および蛍光体D(1,2)が標識された60塩基長のDNA断片が含まれ，キャピラリCap(2)に注入するサンプルには，蛍光体D(2,1)が標識された70塩基長のDNA断片，および蛍光体D(2,2)が標識された80塩基長のDNA断片が含まれるように調製すれば良い。
　あるいは，図23は，キャピラリCap(1)とCap(2)に同一組成のサンプルを同じタイミングで注入したとしても，キャピラリCap(1)とCap(2)で異なる電気泳動条件を設定することによっても実現できる。例えば，電気泳動の最中で，Capi(2)の印加電圧を一時的に低下させる，電気泳動時のCapi(2)の温度を低下させる，等である。
　以上では，行列Yを決定するための専用サンプルを予め準備することを想定していたが，必ずしもその必要はない。分析対象となる実サンプルの電気泳動分析の生データにおいて，各キャピラリにおける各蛍光体の蛍光発光が単独で発生する状態が存在し，それを特定することができれば，その時刻の生データを用いて行列Yを決定することができる。
［実施例7］
　図24は，分析装置の一例であるマルチキャピラリ電気泳動装置の構成図である。マルチキャピラリ電気泳動装置は，DNAシーケンスやDNAフラグメント解析を行う分析装置として広く用いられている。図24に示すように，マルチキャピラリ電気泳動装置は，キャピラリ24-1，陰極24-4，陽極24-5，陰極側緩衝液24-6，陽極側緩衝液24-7，ポンプブロック24-9，シリンジ24-11，レーザ光源24-12を備える。本実施例では，4本のキャピラリ24-1を用い，各キャピラリ24-1で異なるサンプルのDNAシーケンスを実施した。DNAシーケンスのサンプルは，4種類の蛍光体で標識されたDNA断片から構成される。以下の（1）～（6）の工程によって，1回の分析セッションを実行した。（1）まず，4本のキャピラリ24-1の試料注入端24-2を陰極側緩衝液24-6に浸し，試料溶出端24-3をポリマブロック24-9を介して陽極側緩衝液24-7に浸した。（2）次に，ポンプブロック24-9のバルブ24-10を閉じ，ポンプブロック24-9に接続されたシリンジ24-11のピストンを押し下げることにより内部のポリマ溶液に加圧し，ポリマ溶液を各キャピラリ24-1の内部に，サンプル溶出端24-3からサンプル注入端24-2に向かって充填した。（3）続いて，バルブ24-10を開け，各キャピラリ24-1にサンプル注入端24-2から異なるサンプルを電界注入した後，陰極24-4と陽極24-5の間に電源24-8により高電圧を印加することにより，キャピラリ電気泳動を開始した。4種類の蛍光体で標識されたDNA断片は，サンプル注入端24-2からサンプル溶出端24-3に向かって電気泳動された。(4）並行して，各キャピラリ24-1の，サンプル注入端24-2から一定距離電気泳動された位置を発光点24-14とし，レーザ光源24-12から発振された，波長505 nmのレーザビーム24-13を各発光点24-14に一括照射した。ここで，発光点24-14近傍の各キャピラリ24-1の被覆を予め除去し，発光点24-14近傍の各キャピラリ24-1を同一平面上に配列し，レーザビーム24-13を，集光してから，上記の配列平面の側方より，配列平面に沿って導入した。（5）そして，4種類の蛍光体で標識されたDNA断片が，各キャピラリ24-1の内部を電気泳動し，発光点24-14を通過する際にレーザビーム24-13の照射によって励起され，蛍光を発光した。つまり，4個の発光点から4種類の蛍光体が蛍光発光し，電気泳動に伴い，それぞれの蛍光強度が時々刻々と変化した。（6）最後に，各発光点から発光される蛍光をセンサ（不図示）により多色検出し，計算機（不図示）により得られた時系列データを解析することによって，各キャピラリに注入されたサンプルのDNAシーケンスを行った。以上の（1）～（6）の工程からなる分析セッションは複数回，繰り返すことができる。例えば，1回目の分析セッションではサンプル(1)～(4)を分析し，2回目の分析セッションではサンプル(5)～(8)を分析し，…，とすることによって，多数の異なるサンプルを分析することができる。
　本実施例では，上記（6）の多色検出を，特許文献1の光学系を用いて行った。すなわち，A＝4個の発光点P(1)～P(4)から発光される蛍光を1個の集光レンズでコリメートし，1個の透過型の回折格子を透過させ，各蛍光の1次回折光を1個の結像レンズで1個の2次元センサ上に結像させた。図25は，4本のキャピラリCap(1)～Cap(4)に標準溶液を充填した場合の，レーザビーム照射による，4個の発光点P(1)～P(4)から発光されるラマン散乱光の波長分散像を含む2次元センサ画像の取得結果である。横軸方向は4本のキャピラリの配列方向，縦軸方向は波長方向である。矢印で示す，縦方向に伸びた4本の筋状の像がそれぞれキャピラリCap(1)～Cap(4)のラマン散乱光の波長分散像である。この結果に基づいて，各発光点について，2次元センサ画像上における波長分散像の波長校正を行い，画素位置と波長の関係を求めた。そして，その結果に基づいて，図26に示すように，2次元センサ画像上における，発光点P(1)～P(4)のそれぞれについて，B＝20個の異なる波長帯の検出領域を設定した。20個の検出領域はそれぞれ，500～700 nmの発光を，10 nm間隔で等分割した波長帯で検出するように設定した。例えば，発光点P(1)について設定された検出領域W(1,1)，W(1,2)，…，W(1,20)はそれぞれ，500～510 nm，510～520 nm，…，690～700 nmの波長帯の発光を検出するように設定し，それぞれの信号強度をX(1,1)，X(1,2)，…，X(1,20)とした。発光点P(2)～P(4)についても同様である。ここで，各発光点P(1)～P(4)の位置および光学系が固定されている限り，2次元センサ画像における，各キャピラリについての画素位置と波長の関係は維持されるため，図26で設定された検出領域は，複数の異なる発光検出，あるいは複数の異なる分析セッションについても有効である。
　本実施例では，C=4種類の蛍光体として，dR110，dR6G，dTAMRA，およびdROXを用い，それぞれをサンガー反応によって調製される末端塩基種がT，C，A，およびGのDNA断片に標識した。dR110，dR6G，dTAMRA，およびdROXの極大蛍光発光波長はそれぞれ，541 nm，568 nm，595 nm，および618 nmであるため，それぞれの発光蛍光は，540～550 nm，560～570 nm，590～600 nm，および610～620 nmの波長帯の検出領域で最も強い強度で検出される。例えば，発光点P(1)におけるdR110，dR6G，dTAMRA，およびdROXの発光蛍光はそれぞれ，検出領域W(1,5)，W(1,7)，W(1,10)，およびW(1,12)で，主として検出される。ただし，各発光蛍光は，スペクトルクロストークによって発光点P(1)のその他の検出領域においても検出され，空間クロストークおよびスペクトルクロストークによって発光点P(2)～P(3)の検出領域においても弱い強度で検出される。以降，発光点P(1)におけるdR110，dR6G，dTAMRA，およびdROXをそれぞれ，D(1,1)，D(1,2)，D(1,3)，およびD(1,4)，また，それらの濃度をZ(1,1)，Z(1,2)，Z(1,3)，およびZ(1,4)とした。発光点P(2)～P(4)についても同様である。
　このとき，（数9）における（数10）～（数12）は次式となる。
Figure JPOXMLDOC01-appb-M000029
Figure JPOXMLDOC01-appb-M000030
Figure JPOXMLDOC01-appb-M000031
　ここで，Xは80行1列，Yは80行16列，Zは16行1列の行列である。（数23）～（数25）を（数14）に代入して得られる空間補正と色変換を実施し，空間クロストークとスペクトルクロストークの両方を解消し，各キャピラリにおける4種類の蛍光体の濃度の時系列データを取得することができた。その結果，各キャピラリに注入された異なるサンプルのDNAシーケンスを実施することができた。尚，本実施例に対応する従来法は，各キャピラリ毎に（数8）の色変換を実施することである。従来法では空間クロストークを解消できない。
　（数23）～（数25），およびこれらから導出される（数14）は，各発光点の位置，光学系，および用いる蛍光体が固定されている限り，複数の異なる発光検出，あるいは複数の異なる分析セッションについても有効である。
　行列Yの決定方法は既に述べた通りである。ただし，80×16個の要素をすべて設定することは必ずしも必要ない。例えば，他の要素と比較して十分に絶対値が小さい要素をゼロに置き換え，関連する計算を簡略化しても構わない。また，空間クロストークが及ぶ範囲が限定的である場合は，限定された範囲で（数9）～（数12）を定義しても構わない。この場合，必要に応じて，限定された範囲を順次スライドさせ，より広範囲の分析を行うこともできる。例えば，空間クロストークが及ぶ範囲が両隣のキャピラリに限定できる場合は，2本以上離れたキャピラリについての空間クロストークおよびスペクトルクロストークは無視して，関連する計算を省いても構わない。この場合さらに，行列Yを決定するために各発光点で単独に発光させる工程を，互いに空間クロストークの影響を及ぼし合わない複数の発光点について同時に発光させる工程に簡略化することも可能である。
［実施例8］
　本実施例では，［実施例7］と同様に，図24のマルチキャピラリ電気泳動装置を用いたDNAシーケンスを行った。ただし，4個の発光点から発光される，4種類の蛍光体の蛍光発光を，各時刻において，特許文献2の光学系によって多色検出した。図27は，本光学系の構成と，本光学系により取得した4個の発光点の4分割像を示す。図27の光学系は，集光レンズアレイ27-1，ロングパスフィルタ27-3，ダイクロイックミラー27-4，27-5，27-6，および27-7，並びに2次元センサ27-8を備える。図27(a)および図27(b)に示すように，まず，A＝4個の発光点P(1)～P(4)から発光される蛍光をそれぞれ，集光レンズアレイ27-1を構成する4個の集光レンズ27-2で個別にコリメートして4本の光束27-9とし，1個のロングパスフィルタ27-3を一括して透過させることで，レーザビーム光をカットした。次に，4個のダイクロイックミラー27-4，27-5，27-6，および27-7で構成される1組のダイクロイックミラーアレイに各光束27-9を一括して入射することで，各光束27-9を，B=4個の分割光束27-10，27-11，27-12，および27-13にそれぞれ分割し，ダイクロイックミラーアレイから出射させた。ここで，各ダイクロイックミラーの分光特性を調整することによって，分割光束27-10，27-11，27-12，および27-13がそれぞれ，520～550 nm，550～580 nm，580～610 nm，および610～640 nmの波長帯の光を有するようにした。最後に，生成されたA×B＝4×4＝16個の分割光束を2次元センサ27-8に入射させ，図27(c)に示すように，16個の分割像W(1,1)～W(4,4)を得た。ここで，4本のキャピラリCap(1)～Cap(4)の配列方向，すなわち，発光点P(1)～P(4)の配列方向と，ダイクロイックミラーアレイによる分割方向，すなわち波長方向は互いに垂直であるため，図27(c)に示すように，2次元センサ画像27-14の中に，分割像W(1,1)～W(4,4)が互いに重なることなく，整列した状態で検出された。例えば，発光点P(1)から発光される蛍光は，W(1,1)，W(1,2)，W(1,3)，およびW(1,4)の分割像として検出され，それぞれ520～550 nm，550～580 nm，580～610 nm，および610～640 nmの波長帯の光成分を有した。発光点P(2)～P(4)についても同様である。そこで，2次元センサ画像27-14の中の分割像W(1,1)～W(4,4)をそれぞれ検出領域に設定した。また，各検出領域で検出した信号強度をそれぞれX(1,1)～X(4,4)とした。ここで，各発光点P(1)～P(4)の位置および光学系が固定されている限り，2次元センサ画像27-14における，各キャピラリについての画素位置と波長の関係は維持されるため，図27(c)で設定された検出領域は，複数の異なる発光検出，あるいは複数の異なる分析セッションについても有効である。
　本実施例では，［実施例7］と同様に，C=4種類の蛍光体として，dR110，dR6G，dTAMRA，およびdROXを用い，それぞれをサンガー反応によって調製される末端塩基種がT，C，A，およびGのDNA断片に標識した。dR110，dR6G，dTAMRA，およびdROXの極大蛍光発光波長はそれぞれ，541 nm，568 nm，595 nm，および618 nmであるため，それぞれの発光蛍光は，520～550 nm，550～580 nm，580～610 nm，および610～640 nmの波長帯の光成分を有する分割像の検出領域で最も強い強度で検出される。例えば，発光点P(1)におけるdR110，dR6G，dTAMRA，およびdROXの発光蛍光はそれぞれ，検出領域W(1,1)，W(1,2)，W(1,3)，およびW(1,4)で，主として検出される。ただし，各発光蛍光は，スペクトルクロストークによって発光点P(1)のその他の検出領域においても検出され，空間クロストークおよびスペクトルクロストークによって発光点P(2)～P(3)の検出領域においても検出される。以降，発光点P(1)におけるdR110，dR6G，dTAMRA，およびdROXをそれぞれ，蛍光体D(1,1)，D(1,2)，D(1,3)，およびD(1,4)，また，それらの濃度をZ(1,1)，Z(1,2)，Z(1,3)，およびZ(1,4)とした。発光点P(2)～P(4)についても同様である。
　このとき，（数9）における（数10）～（数12）は次式となる。
Figure JPOXMLDOC01-appb-M000032
Figure JPOXMLDOC01-appb-M000033
Figure JPOXMLDOC01-appb-M000034
　ここで，Xは16行1列，Yは16行16列，Zは16行1列の行列である。（数26）～（数28）を（数14）に代入することによって，空間クロストークとスペクトルクロストークの両方を解消し，各キャピラリにおける4種類の蛍光体の濃度の時系列データを取得することができた。その結果，各キャピラリに注入された異なるサンプルのDNAシーケンスを実施することができた。尚，本実施例に対応する従来法は，各キャピラリ毎に（数8）の色変換を実施することである。従来法では空間クロストークを解消できない。
　（数26）～（数28），およびこれらから導出される（数14）は，各発光点の位置，光学系，および用いる蛍光体が固定されている限り，複数の異なる発光検出，あるいは複数の異なる分析セッションについても有効である。
［実施例9］
　図28は，複数の反応槽を平面上に配列させたマルチチャンネルの，各チャンネルで異なるDNA断片を鋳型とする相補鎖伸長反応を1塩基単位で行い，伸長反応で相補鎖に取り込まれる4種類の塩基に標識されたC=4種類の蛍光体の発光蛍光を多色検出することによって，各DNA断片のDNAシーケンスを行うマルチチャンネル伸長反応装置の構成図，および複数の2次元センサで得られる2次元センサ画像の模式図を示す。マルチチャンネル伸長反応装置は，レーザ光源28-1，ダイクロイックミラー28-3，レンズ28-4，ダイクロイックミラー28-8，28-9および28-10，レンズ28-11，第1の2次元センサ28-12，レンズ28-13，第2の2次元センサ28-14，レンズ28-15，第3の2次元センサ28-16，レンズ28-17並びに第4の2次元センサ28-18を備える。塩基種T，C，A，およびGに標識された蛍光体のの極大蛍光発光波長はそれぞれ，535 nm，565 nm，595 nm，および625 nmとした。まず，レーザ光源28-1から発振されたレーザビーム28-2を，ダイクロイックミラー28-3を透過させた後，レンズ28-4で集光させてマルチチャンネルのサンプル28-5に照射した。次に，各チャンネルで，レーザビーム照射で励起された各蛍光体の蛍光発光28-6を，レンズ28-4で一括してコリメートした後，得られた光束28-7をダイクロイックミラー28-3で反射させた。ダイクロイックミラー28-3は，レーザビーム光を透過させ，発光蛍光を反射させる分光特性を有する。続いて，3種類のダイクロイックミラー28-8，28-9，および28-10によって，光束28-7を4個の異なる波長成分を有する4個の光束に分割した。そして，第1の分割光束をレンズ28-11によって第1の2次元センサ28-12上に結像させ，第2の分割光束をレンズ28-13によって第2の2次元センサ28-14上に結像させ，第3の分割光束をレンズ28-15によって第3の2次元センサ28-16上に結像させ，第4の分割光束をレンズ28-17によって第4の2次元センサ28-18上に結像させた。
　サンプル28-19は，サンプル28-5を正面から観察した模式図である。A=5個のチャンネルがそれぞれ発光点P(1)～P(5)を形成した。実際には，より多数のチャンネルがサンプル上に存在し，各チャンネルからの発光蛍光を上記の光学系によって一括検出したが，ここでは，注目する発光点P(3)と，発光点P(3)に直接的に空間クロストークの影響を与える，発光点P(3)の周囲に存在する発光点P(1)，P(2)，P(4)，P(5)のみを解析対象とするため，サンプル28-19には簡略的に発光点P(1)～P(5)のみを描いた。発光点P(3)以外の発光点に注目する場合は，その発光点と，その周辺の発光点を対象として同様に解析すれば良い。いずれの発光点に注目する場合も，蛍光検出の後で，同じ蛍光検出データに対して個別に解析することができる。第1の2次元センサ画像28-20は，第1の2次元センサ28-12によって取得したサンプル28-19の分割像であり，発光点P(1)～P(5)のそれぞれの結像点に検出領域W(1,1)～W(5,1)を設定し，それぞれの信号強度をX(1,1)～X(5,1)とした。第1の2次元センサ画像28-20では，波長帯520～550 nmの蛍光成分を検出した。第2の2次元センサ画像28-21は，第2の2次元センサ28-14によって取得したサンプル28-19の分割像であり，発光点P(1)～P(5)のそれぞれの結像点に検出領域W(1,2)～W(5,2)を設定した。第2の2次元センサ画像28-21では，波長帯550～580 nmの蛍光成分を検出した。第3の2次元センサ画像28-22は，第3の2次元センサ28-16によって取得したサンプル28-19の分割像であり，発光点P(1)～P(5)のそれぞれの結像点に検出領域W(1,3)～W(5,3)を設定した。第3の2次元センサ画像28-22では，波長帯580～610 nmの蛍光成分を検出した。第4の2次元センサ画像28-23は，第4の2次元センサ28-18によって取得したサンプル28-19の分割像であり，発光点P(1)～P(5)のそれぞれの結像点に検出領域W(1,4)～W(5,4)を設定した。第1の2次元センサ画像28-23では，波長帯610～640 nmの蛍光成分を検出した。すなわち，各発光点からの発光蛍光を4個の2次元センサを用いてB=4色検出した。例えば，発光点P(3)における4種類の蛍光体の発光蛍光は波長順にそれぞれ，検出領域W(3,1)，W(3,2)，W(3,3)，およびW(3,4)で，主として検出される。ただし，各発光蛍光は，スペクトルクロストークによって発光点P(3)のその他の検出領域においても検出され，空間クロストークおよびスペクトルクロストークによって発光点P(1)，P(2)，P(4)，P(5)の検出領域においても検出される。以降，発光点P(3)における4種類の蛍光体を波長順にそれぞれ，D(3,1)，D(3,2)，D(3,3)，およびD(3,4)，また，それらの濃度をZ(3,1)，Z(3,2)，Z(3,3)，およびZ(3,4)とした。発光点P(1)，P(2)，P(4)，P(5)についても同様である。
　このとき，（数9）における（数10）～（数12）は次式となる。
Figure JPOXMLDOC01-appb-M000035
Figure JPOXMLDOC01-appb-M000036
Figure JPOXMLDOC01-appb-M000037
　ここで，Xは20行1列，Yは20行20列，Zは20行1列の行列である。（数29）～（数31）を（数14）に代入することによって，空間クロストークとスペクトルクロストークの両方を解消し，各チャンネルにおける4種類の蛍光体の濃度の時系列データを取得することができた。ただし，本実施例では，発光点P(3)のみに注目しているため，行列Zの内，Z(3,1)，Z(3,2)，Z(3,3)，およびZ(3,4)の要素のみを抜き出した。その結果，発光点P(3)のチャンネルにおいてDNAシーケンスを実施することができた。他の発光点に注目して，対応するチャンネルでDNAシーケンスを行う場合も，同様な方法を用いることができる。
［実施例10］
　本実施例では，従来の方法と，本開示の方法をフローチャートで纏める。図29は，従来法の1回の分析セッションを示すフローチャートである。図29に示すように，従来法の分析セッションは，サンプル29-1を分析する分析装置，計算機，表示装置及びデータベースを有する分析システムにより実行される。分析装置は、サンプル29-1からの光が入射するセンサ（不図示）を有する。まず，C種類（Cは1以上の整数）の蛍光体で標識されたA種類（Aは2以上の整数）のサンプル29-1を分析装置に入力する。次に，分析装置において，工程29-2により，各サンプルを並列に分析し，それぞれの分析について，かつ各時刻において，B種類（Bは1以上の整数）の波長帯で発光蛍光を検出し，A×B個の蛍光強度X(a,b)の時系列生データを取得し，これらの時系列生データが計算機に送られる。ここで，a＝1，2，…，およびA，b＝1，2，…，およびBである。続いて，計算機において，工程29-3により，分析(a0)について，データベースに格納されているC行B列の行列Y-(a0)を用いた色変換(a0)によってスペクトルクロストークを解消し，各時刻において，蛍光強度X(a0,b)からC種類の蛍光体の濃度Z(a0,c)を求める。ここで，a0＝1，2，…，またはA，c＝1，2，…，およびCである。最後に，表示装置において，工程29-4により，Z(a0,c)の時系列色変換データを出力する。工程29-3および工程29-4は，すべてのa0について実施される。図29に示されている通り，従来法は，A個の分析および解析が，それぞれ独立に行われていることが特徴である。尚，（数1）～（数4）では，蛍光強度X(a,b)はX(b)，濃度Z(a0,c)はZ(c)と表記している。
　図30は，本法の1回の分析セッションを示すフローチャートである。図29と同様に，図30に示す本実施例の分析セッションは，サンプル29-1を分析する分析装置，計算機，表示装置及びデータベースを有する分析システムにより実行される。分析装置は、サンプル29-1からの光が入射するセンサ（不図示）を有する。まず，C種類（Cは1以上の整数）の蛍光体で標識されたA種類（Aは2以上の整数）のサンプル30-1を分析装置に入力する。次に，分析装置において，工程30-2により，各サンプルを並列に分析し，それぞれの分析について，かつ各時刻において，B種類（Bは1以上の整数）の波長帯で発光蛍光を検出し，A×B個の蛍光強度X(a,b)の時系列生データを取得し，これらの時系列生データが計算機に送られる。ここで，a＝1，2，…，およびA，b＝1，2，…，およびBである。ただし，図29と異なり，工程30-2に点線矢印で示すように，異なる分析の間で発光蛍光の空間クロストークが有意に存在している。続いて，計算機において，工程30-3により，すべての分析(a)について，データベースに格納されている(A×C)行(A×B)列の行列Y-を用いた色変換+空間補正によってスペクトルクロストークおよび空間クロストークを解消し，各時刻において，蛍光強度X(a,b)から，A種類のサンプルまたは分析についての，C種類の蛍光体の濃度Z(a,c)を一括で求める。ここで，a＝1，2，…，およびA，c＝1，2，…，およびCである。最後に，表示装置において，工程30-4により，分析(a)それぞれについて，Z(a,c)の時系列色変換+空間補正データを出力する。図30に示されている通り，本法は，工程30-3における色変換+空間補正を，すべての分析について一括で行うことが特徴である。さらに，色変換+空間補正をすべての分析に共通の行列Y-を用いる点も特徴である。
　図31は，図30より前に行う，データベースに格納する行列Y-の求め方を示すフローチャートである。まず，C種類（Cは1以上の整数）の蛍光体で標識されたA種類（Aは2以上の整数）の行列Y決定用サンプル31-1を分析装置に入力する。このサンプルは，分析(a)における蛍光体(c)の蛍光発光が単独で発生するように調製されている。つまり，a＝1，2，…，およびA，c＝1，2，…，およびCであるが，これらのA×C個のすべての組み合わせについて蛍光発光させるが，2つ以上の組み合わせの蛍光発光が同時に発生しないようなされている。ここでは，行列Y決定用サンプル31-1によって上記を実現するが，分析装置の設定によって，上記を実現しても構わない。次に，分析装置において，工程31-2により，上記のA×C個の単独の蛍光発光をB種類（Bは1以上の整数）の波長帯で検出し，A×B個の蛍光強度X(a,b)の時系列生データを取得し，これらの時系列生データが計算機に送られる。続いて，計算機において，工程31-3により，得られた(A×B)個×(A×C)個のデータから，(A×B)行(A×C)列の行列を導出し，さらに，Yの一般逆行列である，(A×C)行(A×B)列の行列Y-を求める。そして，求めた行列Y-をデータベースに格納し，図30等の以降の分析において行列Y-を活用する。
　図32は，図31で求めた行列Y-を活用しながら，図30に示すA種類のサンプル30-1についての工程30-2～30-4の分析セッションを，複数回繰り返す場合のフローチャートである。ここで，各分析セッションにおける分析対象のA種類のサンプル30-1は相互に異なっている。図32に示されている通り，複数の異なる分析セッションに対して，データベースに格納された同じ行列Y-を活用している点が本法の特徴である。
　図33は，計算機の構成例を示す。計算機は分析装置と接続されている。計算機は，データ解析だけでなく，分析装置の制御も行う。図29～図31ではデータベースおよび表示装置が計算機の外部に描かれているが，図33ではそれらを計算機の内部に描いている。入力部であるキーボードを通じて，データ解析の条件設定や分析装置制御の条件設定を行う。分析装置から出力される蛍光強度X(a,b)の時系列生データが順次メモリに格納される。また，HDDの内部にあるデータベースに格納されている(A×C)行(A×B)列の行列Y-がメモリに格納される。CPUは，メモリに格納された蛍光強度X(a,b)とY-の積を計算し，蛍光体濃度Z(a,c)の時系列色変換+空間補正データを導出し，順次メモリに格納すると同時に，表示部であるモニタに表示する。また，解析結果をネットワークインターフェースNIFを通じてネットワーク上の情報と照合することができる。
［変形例］
　本開示は、上述した実施形態に限定されるものでなく、様々な変形例を含んでいる。例えば、上述した実施形態は、本開示を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備える必要はない。また、ある実施形態の一部を他の実施形態の構成に置き換えることができる。また、ある実施形態の構成に他の実施形態の構成を加えることもできる。また、各実施形態の構成の一部について、他の実施形態の構成の一部を追加、削除又は置換することもできる。
1-1　ピンホール板
1-2　発光点側開口板
1-3　集光レンズ
1-4　センサ側開口板
1-5　色ガラスフィルタ
1-6　2次元センサ
1-7　ハロゲンランプ光
1-8　発光点
1-9　光
1-10　発光像
Cap(1)　キャピラリ1
Cap(2)　キャピラリ2
P(1)　Cap(1)上の発光点
P(2)　Cap(2)上の発光点
D(1,1)　P(1)上の蛍光体1または吸光体1
D(1,2)　P(1)上の蛍光体2または吸光体2
D(2,1)　P(2)上の蛍光体1または吸光体1
D(2,2)　P(2)上の蛍光体2または吸光体2
D(1,3)　P(1)上の蛍光体3または吸光体3
W(1,1)　D(1,1)を主として検出する検出領域
W(1,2)　D(1,2)を主として検出する検出領域
W(2,1)　D(2,1)を主として検出する検出領域
W(2,2)　D(2,2)を主として検出する検出領域
X(1,1)　W(1,1)の信号強度
X(1,2)　W(1,2)の信号強度
X(2,1)　W(2,1)の信号強度
X(2,2)　W(2,2)の信号強度
Z(1,1)　D(1,1)の濃度
Z(1,2)　D(1,2)の濃度
Z(2,1)　D(2,1)の濃度
Z(2,2)　D(2,2)の濃度
S　センサ
LB　レーザビーム
LL　ランプ光
24-1　キャピラリ
24-2　試料注入端
24-3　試料溶出端
24-4　陰極
24-5　陽極
24-6　陰極側緩衝液
24-7　陽極側緩衝液
24-8　電源
24-9　ポリマブロック
24-10　バルブ
24-11　シリンジ
24-12　レーザ光源
24-13　レーザビーム
24-14　発光点
Cap(a)　キャピラリa（a＝1，2，3，および4）
P(a)　キャピラリa上の発光点（a＝1，2，3，および4）
W(a,b)　P(a)からの発光の波長帯bの検出領域（a＝1，2，…，および4）（b＝1，2，…，および20）
27-1　集光レンズアレイ
27-2　集光レンズ
27-3　ロングパスフィルタ
27-4，27-5，27-6，および27-7　ダイクロイックミラー
27-8　2次元センサ
27-9　光束
27-10，27-11，27-12，および27-13　分割光束
27-14　2次元センサ画像
P(a)　発光点a（a＝1，2，3，4，および5）
W(a,b)　P(a)からの発光の波長帯bの検出領域（a＝1，2，3，4，および5）（b＝1，2，3，および4）
28-1　レーザ光源
28-2　レーザビーム
28-3，28-8，28-9，および28-10　ダイクロイックミラー
28-4，28-11，28-13，28-15，および28-17　レンズ
28-5，および28-19　サンプル
28-6　蛍光発光
28-7　光束
28-12，28-14，28-16，および28-18　2次元センサ
28-20，28-21，28-22，および28-23　2次元センサ画像
29-1　サンプル
29-2　分析装置における工程
29-3　計算機における工程
29-4　表示装置における工程
30-1　サンプル
30-2　分析装置における工程
30-3　計算機における工程
30-4　表示装置における工程
31-1　サンプル
31-2　分析装置における工程
31-3　計算機における工程

Claims (40)

  1. 　位置が固定されており，複数種類の発光体が発光する複数の発光点と，
    　前記複数の発光点からの発光をそれぞれ複数の波長帯で検出するセンサと，
    　前記センサの検出信号を処理する計算機と，を備え，
    　前記計算機は，
    　前記複数の発光点および前記複数の波長帯の信号強度を用いた演算によって，前記複数の発光点および前記複数の波長帯の信号強度の間に存在する空間クロストークおよびスペクトルクロストークを低減し，前記複数の発光点のそれぞれにおける前記複数種類の発光体のそれぞれの濃度を導出することを特徴とする分析装置。
  2. 　請求項1に記載の分析装置において，
    　前記計算機は，
    　前記複数の発光点のそれぞれにおける前記複数種類の発光体のそれぞれの濃度が時間変化するとき，各時刻において，その都度得られた信号強度を用いた前記演算を行うことによって，前記複数の発光点のそれぞれにおける前記複数種類の発光体のそれぞれの濃度の時系列を導出することを特徴とする分析装置。
  3. 　請求項1または2に記載の分析装置において，
    　前記計算機は，
    　異なるタイミングで異なる分析を行うとき，その都度得られた信号強度を用いた前記演算を行うことによって，前記複数の発光点のそれぞれにおける前記複数種類の発光体のそれぞれの濃度または濃度の時系列を導出することを特徴とする分析装置。
  4. 　それぞれC種類（Cは1以上の整数）の蛍光体D(a,c)（c＝1，2，…，C）が蛍光を発光するA個（Aは2以上の整数）の発光点P(a)（a＝1，2，…，A）と，
    　前記A個の発光点P(a)からの前記蛍光をそれぞれB種類（Bは1以上の整数）の波長帯の検出領域W(a,b)（b＝1，2，…，B）で検出するセンサと，
    　前記センサの検出信号を処理する計算機と，を備え，
    　前記計算機は，
    　任意の時刻における前記検出領域W(a,b)の信号強度をX(a,b)，前記蛍光体D(a,c)の濃度をZ(a,c)とするとき，任意の時刻において，予め定められた演算式によって，すべてのaおよびbの組み合わせの信号強度X(a,b)から，すべてのaおよびcの組み合わせの濃度Z(a,c)を導出し，空間クロストークおよびスペクトルクロストークを低減することを特徴とする分析装置。
  5. 　請求項4に記載の分析装置において，
    　前記信号強度X(a,b)を要素とする(A×B)行1列の行列Xを，
    Figure JPOXMLDOC01-appb-M000001
     とし，
    　Z(a,c)を要素とする(A×C)行1列の行列Zを，
    Figure JPOXMLDOC01-appb-M000002
     とし，
    　X=Y×Zの関係を満たす，Y(a,b)(a,c)を要素とする(A×B)行(A×C)列の行列Yを
    Figure JPOXMLDOC01-appb-M000003
     とするとき，前記計算機は，(A×C)行(A×B)列の行列Yの一般逆行列Y-を予め求めておき，Z=Y-×Xを上記の演算式とすることを特徴とする分析装置。
  6. 　請求項5に記載の分析装置において，
    　前記計算機は，
    　行列Zのひとつの要素Z(a,c)が単独で存在し，他の要素がゼロと見なせる状態において，行列Xのすべての要素X(a,b)を取得し，X(a,b)に比例した値を行列Yのひとつの列の要素Y(a,b)(a,c)とする工程を，行列Zのすべてaおよびcの組み合わせについて行うことで，行列Yを取得することを特徴とする分析装置。
  7. 　請求項5に記載の分析装置において，
    　前記計算機は，
    　異なるタイミングで異なる分析を行うとき，同一の行列YおよびY-を用いることを特徴とする分析装置。
  8. 　請求項4～7に記載の分析装置において，
    　前記A個の発光点がそれぞれ設けられるA本のキャピラリと，
    　前記A本のキャピラリにレーザビームを照射する光源と，をさらに備え，
    　前記A個の発光点が前記レーザビームで連続的に照射され，前記C種類の蛍光体で標識されたサンプルが前記A本のキャピラリのそれぞれの内部を電気泳動されることで前記A個の発光点P(a)のそれぞれを通過し，その際に，前記C種類の蛍光体のそれぞれが前記レーザビームで励起されて前記蛍光を発光することを特徴とする分析装置。
  9. 　位置が固定されており，複数種類の蛍光体が蛍光を発光する複数の発光点と，
    　前記複数の発光点からの前記蛍光をそれぞれ複数の波長帯で検出するセンサと，
    　前記センサの検出信号を処理する計算機と，を備え，
    　前記計算機は，
    　1個の発光点で1種類の蛍光体が単独で発光する状態を，前記複数の発光点と前記複数種類の蛍光体のすべての組み合わせについて，異なる時刻で引き起こす処理を実行することを特徴とする分析装置。
  10. 　請求項9に記載の分析装置において，
    　前記計算機は，
    　前記複数の発光点および前記複数の波長帯の信号強度，および前記処理で得られた情報を用いた演算によって，前記複数の発光点および前記複数の波長帯の信号強度の間に存在する空間クロストークおよびスペクトルクロストークを低減し，前記複数の発光点のそれぞれにおける前記複数種類の蛍光体のそれぞれの濃度を導出することを特徴とする分析装置。
  11. 　それぞれC種類（Cは1以上の整数）の蛍光体D(a,c)（c＝1，2，…，C）が蛍光を発光するA個（Aは2以上の整数）の発光点P(a)（a＝1，2，…，A）と，
    　前記A個の発光点P(a)からの前記蛍光をそれぞれB種類（Bは1以上の整数）の波長帯の検出領域W(a,b)（b＝1，2，…，B）で検出するセンサと，
    　前記センサの検出信号を処理する計算機と，を備え，
    　前記計算機は，
    　1個の発光点P(a)で1種類の蛍光体D(a,c)が単独で発光する状態を，すべてのaおよびcの組み合わせについて，異なる時刻で引き起こす処理を実行することを特徴とする分析装置。
  12. 　請求項11に記載の分析装置において，
    　前記計算機は，
    　任意の時刻における前記検出領域W(a,b)の信号強度をX(a,b)，前記蛍光体D(a,c)の濃度をZ(a,c)とするとき，すべてのaおよびbの組み合わせの信号強度X(a,b)，および前記処理で得られた情報を用いた演算によって，すべての前記A個の発光点P(a)について，任意の時刻において，空間クロストークおよびスペクトルクロストークを低減し，前記濃度Z(a,c)を導出することを特徴とする分析装置。
  13. 　位置が固定されており，光を透過する複数の吸光点と，
    　前記複数の吸光点のそれぞれに存在する複数種類の吸光体の光吸収によって発生する吸光度を，それぞれ複数の波長帯で検出するセンサと，
    　前記センサの検出信号を処理する計算機と，を備え，
    　前記計算機は，
    　前記複数の吸光点および前記複数の波長帯の吸光度を用いた演算によって，前記複数の吸光点および前記複数の波長帯の吸光度の間に存在する空間クロストークおよびスペクトルクロストークを低減し，前記複数の吸光点のそれぞれにおける前記複数種類の吸光体のそれぞれの濃度を導出することを特徴とする分析装置。
  14. 　請求項13に記載の分析装置において，
    　前記計算機は，
    　前記複数の吸光点のそれぞれにおける前記複数種類の吸光体のそれぞれの濃度が時間変化するとき，各時刻において，その都度得られた吸光度を用いた前記演算を行うことによって，前記複数の吸光点のそれぞれにおける前記複数種類の吸光体のそれぞれの濃度の時系列を導出することを特徴とする分析装置。
  15. 　位置が固定されており，発光体が発光する複数の発光点と，
    　前記複数の発光点からの発光をそれぞれ検出するセンサと，
    　前記センサの検出信号を処理する計算機と，を備え，
    　前記計算機は，
    　前記複数の発光点の信号強度を用いた演算によって，前記複数の発光点の信号強度の間に存在する空間クロストークを低減し，前記複数の発光点のそれぞれにおける前記発光体の濃度を導出することを特徴とする分析装置。
  16. 　請求項15に記載の分析装置において，
    　前記計算機は，
    　前記複数の発光点のそれぞれにおける前記発光体の濃度が時間変化するとき，各時刻において，その都度得られた信号強度を用いた前記演算を行うことによって，前記複数の発光点のそれぞれにおける前記発光体の濃度の時系列を導出することを特徴とする分析装置。
  17. 　請求項15または16に記載の分析装置において，
    　前記計算機は，
    　異なるタイミングで異なる分析を行うとき，その都度得られた信号強度を用いた前記演算を行うことによって，前記複数の発光点のそれぞれにおける前記発光体の濃度または前記濃度の時系列を導出することを特徴とする分析装置。
  18. 　位置が固定されており，光を透過する複数の吸光点と，
    　前記複数の吸光点のそれぞれに存在する吸光体の光吸収によって発生する吸光度をそれぞれ検出するセンサと，
    　前記センサの検出信号を処理する計算機と，を備え，
    　前記計算機は，
    　前記複数の吸光点の吸光度を用いた演算によって，前記複数の吸光点の吸光度の間に存在する空間クロストークを低減し，前記複数の吸光点のそれぞれにおける前記吸光体の濃度を導出することを特徴とする分析装置。
  19. 　請求項18に記載の分析装置において，
    　前記計算機は，
    　前記複数の吸光点のそれぞれにおける前記吸光体の濃度が時間変化するとき，各時刻において，その都度得られた吸光度を用いた前記演算を行うことによって，前記複数の吸光点のそれぞれにおける前記吸光体の濃度の時系列を導出することを特徴とする分析装置。
  20. 　請求項18または19に記載の分析装置において，
    　前記計算機は，
    　異なるタイミングで異なる分析を行うとき，その都度得られた吸光度を用いた前記演算を行うことによって，前記複数の吸光点のそれぞれにおける前記吸光体の濃度または前記濃度の時系列を導出することを特徴とする分析装置。
  21. 　位置が固定されており，複数種類の発光体が発光する複数の発光点を準備することと，
    　センサにより，前記複数の発光点からの発光をそれぞれ複数の波長帯で検出することと，
    　計算機により，前記センサの検出信号を処理することと，を含み，
    　前記処理することは，
    　前記複数の発光点および前記複数の波長帯の信号強度を用いた演算によって，前記複数の発光点および前記複数の波長帯の信号強度の間に存在する空間クロストークおよびスペクトルクロストークを低減することと，
    　前記複数の発光点のそれぞれにおける前記複数種類の発光体のそれぞれの濃度を導出することと，を含む分析方法。
  22. 　請求項21に記載の分析方法において，
    　前記処理することは，
    　前記複数の発光点のそれぞれにおける前記複数種類の発光体のそれぞれの濃度が時間変化するとき，各時刻において，その都度得られた信号強度を用いた前記演算を行うことによって，前記複数の発光点のそれぞれにおける前記複数種類の発光体のそれぞれの濃度の時系列を導出することを含む分析方法。
  23. 　請求項21または22に記載の分析方法において，
    　前記処理することは，
    　異なるタイミングで異なる分析を行うとき，その都度得られた信号強度を用いた前記演算を行うことによって，前記複数の発光点のそれぞれにおける前記複数種類の発光体のそれぞれの濃度または濃度の時系列を導出することを含む分析方法。
  24. 　それぞれC種類（Cは1以上の整数）の蛍光体D(a,c)（c＝1，2，…，C）が蛍光を発光するA個（Aは2以上の整数）の発光点P(a)（a＝1，2，…，A）を準備することと，
    　センサにより，前記A個の発光点P(a)からの前記蛍光をそれぞれB種類（Bは1以上の整数）の波長帯の検出領域W(a,b)（b＝1，2，…，B）で検出することと，
    　計算機により，前記センサの検出信号を処理することと，を含み，
    　前記処理することは，
    　任意の時刻における前記検出領域W(a,b)の信号強度をX(a,b)，前記蛍光体D(a,c)の濃度をZ(a,c)とするとき，任意の時刻において，予め定められた演算式によって，すべてのaおよびbの組み合わせの信号強度X(a,b)から，すべてのaおよびcの組み合わせの濃度Z(a,c)を導出し，空間クロストークおよびスペクトルクロストークを低減することを含む分析方法。
  25. 　請求項24に記載の分析方法において，
    　前記信号強度X(a,b)を要素とする(A×B)行1列の行列Xを，
    Figure JPOXMLDOC01-appb-M000004
     とし，
    　Z(a,c)を要素とする(A×C)行1列の行列Zを，
    Figure JPOXMLDOC01-appb-M000005
     とし，
    　X=Y×Zの関係を満たす，Y(a,b)(a,c)を要素とする(A×B)行(A×C)列の行列Yを
    Figure JPOXMLDOC01-appb-M000006
     とするとき，前記処理することは，(A×C)行(A×B)列の行列Yの一般逆行列Y-を予め求めておき，Z=Y-×Xを上記の演算式とすることを含む分析方法。
  26. 　請求項25に記載の分析方法において，
    　前記処理することは，
    　行列Zのひとつの要素Z(a,c)が単独で存在し，他の要素がゼロと見なせる状態において，行列Xのすべての要素X(a,b)を取得し，X(a,b)に比例した値を行列Yのひとつの列の要素Y(a,b)(a,c)とする工程を，行列Zのすべてaおよびcの組み合わせについて行うことで，行列Yを取得することを含む分析方法。
  27. 　請求項25に記載の分析方法において，
    　前記処理することは，
    　異なるタイミングで異なる分析を行うとき，同一の行列YおよびY-を用いることを含む分析方法。
  28. 　請求項24～27に記載の分析方法において，
    　前記A個の発光点は，それぞれA本のキャピラリに設けられ，
    　前記分析方法は，
    　光源により，前記A本のキャピラリにレーザビームを照射することをさらに含み，
    　前記A個の発光点が前記レーザビームで連続的に照射され，前記C種類の蛍光体で標識されたサンプルが前記A本のキャピラリのそれぞれの内部を電気泳動されることで前記A個の発光点P(a)のそれぞれを通過し，その際に，前記C種類の蛍光体のそれぞれが前記レーザビームで励起されて前記蛍光を発光することを特徴とする分析方法。
  29. 　位置が固定されており，複数種類の蛍光体が蛍光を発光する複数の発光点を準備することと，
    　センサにより，前記複数の発光点からの前記蛍光をそれぞれ複数の波長帯で検出することと，
    　計算機により，前記センサの検出信号を処理することと，を含み，
    　前記処理することは，
    　1個の発光点で1種類の蛍光体が単独で発光する状態を，前記複数の発光点と前記複数種類の蛍光体のすべての組み合わせについて，異なる時刻で引き起こす処理を実行することを含む分析方法。
  30. 　請求項29に記載の分析方法において，
    　前記処理することは，
    　前記複数の発光点および前記複数の波長帯の信号強度，および前記処理を実行することで得られた情報を用いた演算によって，前記複数の発光点および前記複数の波長帯の信号強度の間に存在する空間クロストークおよびスペクトルクロストークを低減することと，
    　前記複数の発光点のそれぞれにおける前記複数種類の蛍光体のそれぞれの濃度を導出することと，を含む分析方法。
  31. 　それぞれC種類（Cは1以上の整数）の蛍光体D(a,c)（c＝1，2，…，C）が蛍光を発光するA個（Aは2以上の整数）の発光点P(a)（a＝1，2，…，A）を準備することと，
    　センサにより，前記A個の発光点P(a)からの前記蛍光をそれぞれB種類（Bは1以上の整数）の波長帯の検出領域W(a,b)（b＝1，2，…，B）で検出することと，
    　計算機により，前記センサの検出信号を処理することと，を含み，
    　前記処理することは，
    　1個の発光点P(a)で1種類の蛍光体D(a,c)が単独で発光する状態を，すべてのaおよびcの組み合わせについて，異なる時刻で引き起こす処理を実行することを含む分析方法。
  32. 　請求項31に記載の分析方法において，
    　前記処理することは，
    　任意の時刻における前記検出領域W(a,b)の信号強度をX(a,b)，前記蛍光体D(a,c)の濃度をZ(a,c)とするとき，すべてのaおよびbの組み合わせの信号強度X(a,b)，および前記処理を実行することで得られた情報を用いた演算によって，すべての前記A個の発光点P(a)について，任意の時刻において，空間クロストークおよびスペクトルクロストークを低減することと，
    　前記濃度Z(a,c)を導出することと，を含む分析方法。
  33. 　位置が固定されており，光を透過する複数の吸光点を準備することと，
    　センサにより，前記複数の吸光点のそれぞれに存在する複数種類の吸光体の光吸収によって発生する吸光度を，それぞれ複数の波長帯で検出することと，
    　計算機により，前記センサの検出信号を処理することと，を含み，
    　前記処理することは，
    　前記複数の吸光点および前記複数の波長帯の吸光度を用いた演算によって，前記複数の吸光点および前記複数の波長帯の吸光度の間に存在する空間クロストークおよびスペクトルクロストークを低減することと，
    　前記複数の吸光点のそれぞれにおける前記複数種類の吸光体のそれぞれの濃度を導出することと，を含む分析方法。
  34. 　請求項33に記載の分析方法において，
    　前記処理することは，
    　前記複数の吸光点のそれぞれにおける前記複数種類の吸光体のそれぞれの濃度が時間変化するとき，各時刻において，その都度得られた吸光度を用いた前記演算を行うことによって，前記複数の吸光点のそれぞれにおける前記複数種類の吸光体のそれぞれの濃度の時系列を導出することを含む分析方法。
  35. 　位置が固定されており，発光体が発光する複数の発光点を準備することと，
    　センサにより，前記複数の発光点からの発光をそれぞれ検出することと，
    　計算機により，前記センサの検出信号を処理することと，を含み，
    　前記処理することは，
    　前記複数の発光点の信号強度を用いた演算によって，前記複数の発光点の信号強度の間に存在する空間クロストークを低減することと，
    　前記複数の発光点のそれぞれにおける前記発光体の濃度を導出することと，を含む分析方法。
  36. 　請求項35に記載の分析方法において，
    　前記処理することは，
    　前記複数の発光点のそれぞれにおける前記発光体の濃度が時間変化するとき，各時刻において，その都度得られた信号強度を用いた前記演算を行うことによって，前記複数の発光点のそれぞれにおける前記発光体の濃度の時系列を導出することを含む分析方法。
  37. 　請求項35または36に記載の分析方法において，
    　前記処理することは，
    　異なるタイミングで異なる分析を行うとき，その都度得られた信号強度を用いた前記演算を行うことによって，前記複数の発光点のそれぞれにおける前記発光体の濃度または前記濃度の時系列を導出することを含む分析方法。
  38. 　位置が固定されており，光を透過する複数の吸光点を準備することと，
    　センサにより，前記複数の吸光点のそれぞれに存在する吸光体の光吸収によって発生する吸光度をそれぞれ検出することと，
    　計算機により，前記センサの検出信号を処理することと，を含み，
    　前記処理することは，
    　前記複数の吸光点の吸光度を用いた演算によって，前記複数の吸光点の吸光度の間に存在する空間クロストークを低減することと，
    　前記複数の吸光点のそれぞれにおける前記吸光体の濃度を導出することと，を含む分析方法。
  39. 　請求項38に記載の分析方法において，
    　前記処理することは，
    　前記複数の吸光点のそれぞれにおける前記吸光体の濃度が時間変化するとき，各時刻において，その都度得られた吸光度を用いた前記演算を行うことによって，前記複数の吸光点のそれぞれにおける前記吸光体の濃度の時系列を導出することを含む分析方法。
  40. 　請求項38または39に記載の分析方法において，
    　前記処理することは，
    　異なるタイミングで異なる分析を行うとき，その都度得られた吸光度を用いた前記演算を行うことによって，前記複数の吸光点のそれぞれにおける前記吸光体の濃度または前記濃度の時系列を導出することを含む分析方法。
PCT/JP2020/017202 2019-05-22 2020-04-21 分析装置および分析方法 WO2020235283A1 (ja)

Priority Applications (6)

Application Number Priority Date Filing Date Title
SG11202112784UA SG11202112784UA (en) 2019-05-22 2020-04-21 Analysis system and analysis method
CN202080036799.3A CN113840902A (zh) 2019-05-22 2020-04-21 分析装置以及分析方法
US17/612,046 US20220236183A1 (en) 2019-05-22 2020-04-21 Analysis device and analysis method
GB2116497.5A GB2597202B (en) 2019-05-22 2020-04-21 Analysis system and analysis method
JP2021520664A JP7282880B2 (ja) 2019-05-22 2020-04-21 分析装置および分析方法
DE112020002033.9T DE112020002033T5 (de) 2019-05-22 2020-04-21 Analysesystem und Analyseverfahren

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2019096011 2019-05-22
JP2019-096011 2019-05-22

Publications (1)

Publication Number Publication Date
WO2020235283A1 true WO2020235283A1 (ja) 2020-11-26

Family

ID=73458607

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/017202 WO2020235283A1 (ja) 2019-05-22 2020-04-21 分析装置および分析方法

Country Status (7)

Country Link
US (1) US20220236183A1 (ja)
JP (2) JP7282880B2 (ja)
CN (1) CN113840902A (ja)
DE (1) DE112020002033T5 (ja)
GB (5) GB2619179B (ja)
SG (1) SG11202112784UA (ja)
WO (1) WO2020235283A1 (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004532969A (ja) * 2000-10-12 2004-10-28 アムニス コーポレイション レポーターラベルビードを読み取るための方法および装置
JP2005524051A (ja) * 2001-08-28 2005-08-11 ベイラー　カレッジ　オブ　メディスン カラーブラインド蛍光のためのパルスマルチライン励起法
JP2006515428A (ja) * 2003-02-27 2006-05-25 コンビセップ インコーポレーテッド 電気泳動計器のためのロボットフレンドリな外部ローディングシステムおよび方法
US20140303016A1 (en) * 2007-06-06 2014-10-09 Pacific Biosciences Of California, Inc. Methods and Processes for Calling Bases in Sequence by Incorporation Methods
WO2018151843A2 (en) * 2017-02-17 2018-08-23 Life Technologies Corporation Automated quality control and spectral error correction for sample analysis instruments

Family Cites Families (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5096807A (en) * 1985-03-06 1992-03-17 Murex Corporation Imaging immunoassay detection system with background compensation and its use
JPS62260652A (ja) 1985-12-27 1987-11-12 コニカ株式会社 写真処理液容器
WO2000055882A1 (en) * 1999-03-18 2000-09-21 Cambridge Research & Instrumentation Inc. High-efficiency multiple probe imaging system
JP4175735B2 (ja) * 1999-05-12 2008-11-05 独立行政法人理化学研究所 マルチキャピラリー電気泳動装置
CA2436098C (en) * 2001-01-26 2014-04-01 Biocal Technology, Inc. Optical detection in a multi-channel bio-separation system
JP3897277B2 (ja) 2001-04-02 2007-03-22 株式会社日立製作所 キャピラリアレイ及びキャピラリアレイ光検出装置
EP1389956B1 (en) * 2001-04-25 2012-10-31 Amnis Corporation Method and apparatus for correcting crosstalk and spatial resolution for multichannel imaging
US7222059B2 (en) * 2001-11-15 2007-05-22 Siemens Medical Solutions Diagnostics Electrophoretic trace simulator
WO2006014494A2 (en) * 2004-07-07 2006-02-09 Corcoran Timothy C Multiple-label fluorescence imaging using excitation-emission matrices
WO2006015251A2 (en) * 2004-07-29 2006-02-09 The Research Foundation Of State University Of New York System and method for cross-talk cancellation in a multilane fluorescence detector
US20060094109A1 (en) * 2004-11-02 2006-05-04 Immunivest Corporation Device and method for analytical cell imaging
US7715001B2 (en) * 2006-02-13 2010-05-11 Pacific Biosciences Of California, Inc. Methods and systems for simultaneous real-time monitoring of optical signals from multiple sources
US7692783B2 (en) * 2006-02-13 2010-04-06 Pacific Biosciences Of California Methods and systems for simultaneous real-time monitoring of optical signals from multiple sources
US20080033677A1 (en) * 2006-06-27 2008-02-07 Applera Corporation Methods And System For Compensating For Spatial Cross-Talk
EP2044432B1 (en) * 2006-06-28 2023-03-29 Life Technologies Corporation Minimizing effects of dye crosstalk
JP4964568B2 (ja) * 2006-11-24 2012-07-04 浜松ホトニクス株式会社 蛍光検出装置、蛍光検出方法および蛍光検出プログラム
EP2155855B1 (en) * 2007-06-06 2016-10-12 Pacific Biosciences of California, Inc. Methods and processes for calling bases in sequence by incorporation methods
US8475640B2 (en) * 2009-04-20 2013-07-02 Advanced Analytical Technologies, Inc. Method for multiplexed capillary electrophoresis signal cross-talk correction
US8817115B1 (en) * 2010-05-05 2014-08-26 Amnis Corporation Spatial alignment of image data from a multichannel detector using a reference image
WO2011143507A2 (en) * 2010-05-12 2011-11-17 Li-Cor, Inc. High dynamic range scanning with reduced channel cross-talk
US20120015825A1 (en) * 2010-07-06 2012-01-19 Pacific Biosciences Of California, Inc. Analytical systems and methods with software mask
IT1403792B1 (it) * 2010-12-30 2013-10-31 St Microelectronics Srl Analizzatore per analisi biochimiche e metodo per la determinazione di concentrazioni di sostanze fluorescenti in una soluzione
EP2972207B1 (en) * 2013-03-14 2022-09-21 Cytonome/ST, LLC Assemblies and methods for reducing optical crosstalk in particle processing systems
CN105143469B (zh) * 2013-05-24 2018-05-22 株式会社日立高新技术 核酸分析装置及使用其的核酸分析方法
US9170198B2 (en) * 2013-12-19 2015-10-27 Luminex Corporation Crosstalk reduction
JP6351765B2 (ja) * 2015-01-27 2018-07-04 株式会社日立ハイテクノロジーズ 多色蛍光分析装置
JP6975704B2 (ja) * 2015-02-03 2021-12-01 株式会社日立ハイテク 多色検出装置
EP3147650A1 (en) 2015-09-22 2017-03-29 MyCartis N.V. Cross-talk correction in multiplexing analysis of biological sample
WO2018150559A1 (ja) * 2017-02-20 2018-08-23 株式会社日立ハイテクノロジーズ 分析システム及び分析方法
US11237111B2 (en) * 2020-01-30 2022-02-01 Trustees Of Boston University High-speed delay scanning and deep learning techniques for spectroscopic SRS imaging
DE112020006681T5 (de) * 2020-04-16 2022-12-01 Hitachi High-Tech Corporation Elektrophoresesystem

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004532969A (ja) * 2000-10-12 2004-10-28 アムニス コーポレイション レポーターラベルビードを読み取るための方法および装置
JP2005524051A (ja) * 2001-08-28 2005-08-11 ベイラー　カレッジ　オブ　メディスン カラーブラインド蛍光のためのパルスマルチライン励起法
JP2006515428A (ja) * 2003-02-27 2006-05-25 コンビセップ インコーポレーテッド 電気泳動計器のためのロボットフレンドリな外部ローディングシステムおよび方法
US20140303016A1 (en) * 2007-06-06 2014-10-09 Pacific Biosciences Of California, Inc. Methods and Processes for Calling Bases in Sequence by Incorporation Methods
WO2018151843A2 (en) * 2017-02-17 2018-08-23 Life Technologies Corporation Automated quality control and spectral error correction for sample analysis instruments

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YANG, T. H. ET AL.: "Comparative microRNA detection from precursor-microRNA-transfected hepatocellular carcinoma cells by capillary electrophoresis with dual-color laser -induced fluorescence", ELECTROPHORESIS, vol. 33, September 2012 (2012-09-01), pages 2769 - 2776, XP055762842, DOI: 10.1002/elps.201200246 *

Also Published As

Publication number Publication date
GB202312065D0 (en) 2023-09-20
JPWO2020235283A1 (ja) 2020-11-26
GB2597202A (en) 2022-01-19
SG11202112784UA (en) 2021-12-30
GB2619179B (en) 2024-03-27
DE112020002033T5 (de) 2022-01-13
GB2597202B (en) 2023-12-13
GB202311980D0 (en) 2023-09-20
JP2023101563A (ja) 2023-07-21
GB2619177B (en) 2024-03-27
CN113840902A (zh) 2021-12-24
GB202116497D0 (en) 2021-12-29
GB202312067D0 (en) 2023-09-20
GB2619178B (en) 2024-05-01
JP7282880B2 (ja) 2023-05-29
US20220236183A1 (en) 2022-07-28
GB2619176A (en) 2023-11-29
GB202312063D0 (en) 2023-09-20
GB2619179A (en) 2023-11-29
GB2619176B (en) 2024-03-27
GB2619177A (en) 2023-11-29
GB2619178A (en) 2023-11-29

Similar Documents

Publication Publication Date Title
Karger et al. Multiwavelength fluorescence detection for DNA sequencing using capillary electrophoresis
EP1837647B1 (en) Multiplexed capillary electrophoresis method
US10295466B2 (en) Fluorescence intensity correcting method, fluorescence intensity calculating method, and fluorescence intensity calculating apparatus
Galievsky et al. “Getting the best sensitivity from on-capillary fluorescence detection in capillary electrophoresis”–A tutorial
DE112015005826B4 (de) Mehrfarbige-fluoreszenz-analysevorrichtung
EP1677097A1 (en) Method and system for measuring the concentrations of fluorescent dyes
US20240085329A1 (en) Analysis system and analysis method
US20170269032A1 (en) Pixelated 2-dimensional fluorescence detection system
JP2011521247A (ja) 累積時間分解発光二次元ゲル電気泳動法
Li et al. Simple two-color base-calling schemes for DNA sequencing based on standard four-label Sanger chemistry
WO2020235283A1 (ja) 分析装置および分析方法
JPH05118991A (ja) 塩基配列決定方法およびその装置
JP7340095B2 (ja) 電気泳動システム
CN115963094A (zh) 用于荧光检测的多通道光学系统及补偿方法、存储介质
JPWO2020235283A5 (ja)
JP2020180983A (ja) キャピラリ電気泳動装置
Huo et al. Pulsed multi-wavelength excitation using fiber-in-capillary light emitting diode induced fluorescence detection in capillary electrophoresis
JP7016957B2 (ja) 生体ポリマー分析方法、および生体ポリマー分析装置
US20100264028A1 (en) Method for multiplexed capillary electrophoresis signal cross-talk correction
CN110088597A (zh) 用于检测分析物的方法和装置
US20220113253A1 (en) Method of intrinsic spectral analysis and applications thereof
WO2021229700A1 (ja) 電気泳動装置及び分析方法
Li et al. On‐capillary raman spectroscopy of nucleoside/borate complexes separated by isotachophoresis
JP4078324B2 (ja) 電気泳動装置およびキャピラリーアレイ
WO2013181362A1 (en) Multi-color detection system for multiplexed capillary electrophoresis

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: 20809448

Country of ref document: EP

Kind code of ref document: A1

DPE1 Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101)
ENP Entry into the national phase

Ref document number: 2021520664

Country of ref document: JP

Kind code of ref document: A

122 Ep: pct application non-entry in european phase

Ref document number: 20809448

Country of ref document: EP

Kind code of ref document: A1