WO2013133283A1 - 判定方法、判定装置、判定システム、および、プログラム - Google Patents
判定方法、判定装置、判定システム、および、プログラム Download PDFInfo
- Publication number
- WO2013133283A1 WO2013133283A1 PCT/JP2013/056020 JP2013056020W WO2013133283A1 WO 2013133283 A1 WO2013133283 A1 WO 2013133283A1 JP 2013056020 W JP2013056020 W JP 2013056020W WO 2013133283 A1 WO2013133283 A1 WO 2013133283A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- pixel group
- value
- spot
- median
- determination
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/53—Immunoassay; Biospecific binding assay; Materials therefor
- G01N33/5308—Immunoassay; Biospecific binding assay; Materials therefor for analytes not provided for elsewhere, e.g. nucleic acids, uric acid, worms, mites
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/69—Microscopic objects, e.g. biological cells or cellular parts
- G06V20/695—Preprocessing, e.g. image segmentation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N21/645—Specially adapted constructive features of fluorimeters
- G01N21/6456—Spatial resolved fluorescence measurements; Imaging
- G01N21/6458—Fluorescence microscopy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10056—Microscopic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10064—Fluorescence image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30072—Microarray; Biochip, DNA array; Well plate
Definitions
- the present invention relates to a determination method, a determination device, a determination system, and a program.
- the gene expression level can be estimated by obtaining the detection intensity of the label.
- the selective binding amount of the test substance such as the gene expression level requires reliability as data, but the selective binding substance such as complementary DNA immobilized on the carrier as a spot is not uniform.
- the detection intensity may fluctuate due to non-biological factors, such as when there is unevenness or when dust is attached to the spot.
- the variation is evaluated based on the value of CV (coefficient of variation).
- CV coefficient of variation
- SD standard deviation
- the microarray analysis software GenePix Pro (Molecular Devices) compares each intensity value of the pixels included in the spot with the average value of the intensity of the pixels around the spot, and among the pixels included in the spot, When the predetermined ratio satisfies the condition “(intensity of each pixel included in the spot) ⁇ (average pixel intensity around the spot) ⁇ 0”, the spot is determined to be defective (see Non-Patent Document 1). As a result, a case where the background around the spot becomes high due to some abnormality, for example, dust adhering or dirt on the chip is detected.
- the spot uniformity determination method in the conventional method has a problem that the reliability of the selective binding amount of the substance to be inspected cannot be determined appropriately.
- the uniformity evaluation method described in Patent Document 1 has a problem that the uniformity within a single spot cannot be evaluated although the uniformity between a plurality of spots can be evaluated.
- the standard deviation increases when the intensity of some pixels in the pixel group constituting the spot is extremely high or low due to dust adhesion or the like. Since the CV value exceeds the threshold value, the intensity data can be used for a group of pixels other than some of the pixels, but the data is excluded, which is an overdetection.
- a median (median value) in a pixel group in the spot in the image is used.
- an intensity median (spot median) is obtained in about 70 pixel groups (pixel size: 10 ⁇ m square) included in the circle.
- FIG. 1 is a scatter plot showing the detection intensity when the same test substance is hybridized in two DNA chips.
- the horizontal axis represents the detection intensity of each spot on one DNA chip, and the vertical axis represents the detection intensity of each spot on the other DNA chip. That is, the coordinates (X, Y) of one point are the detection intensity (X) measured on one DNA chip and the detection measured on the other DNA chip for a spot on which the same selective binding substance is fixed. Intensity (Y) is represented.
- FIG. 2 is a diagram showing an intensity image (left figure) of the spot 1 and a fluctuation graph (right figure) of the detected intensity along the vertical line of the spot 1.
- FIG. 3 is a diagram showing a spot 2 intensity image (left figure) and a detected intensity variation graph (right figure) along the vertical line of the spot 2.
- a broken-line circle indicates a spot portion, and the detected intensity is represented by a white gradation value.
- the detected intensity is not a uniform circle but is uneven.
- the fluctuation graph (the right figure) shows the fluctuation of the detection intensity in the line that vertically cuts the spot portion across one horizontal pixel and ten or more horizontal pixels, and it fluctuates greatly where it should ideally be a constant intensity.
- the causes of defects are mainly spot defects (such as defects when the selective binding substance is fixed) and dust adhesion.
- the CV value is equal to or greater than a predetermined reference value
- the data is excluded as a defective spot.
- the CV value becomes equal to or higher than the reference value, and the data is excluded. Is done.
- the representative value of the signal intensity of the spot the median of the pixel group in the spot is usually used, and even if an outlier is included in a part of the spot, it hardly fluctuates. There is a possibility that it can be used sufficiently as data. That is, the conventional evaluation method using the CV value has a problem in that although it can be used sufficiently, the defective spot is overdetected and data is excluded.
- the present invention has been made in view of the above, and in data obtained from a microarray experiment or the like, it is possible to appropriately evaluate the reliability of the selective binding amount of a test substance in consideration of non-biological effects.
- An object is to provide a determination method, a determination apparatus, a determination system, and a program.
- the determination method of the present invention is obtained as the detection intensity of the label by binding the labeled analyte to the selective binding substance immobilized as a spot on the carrier.
- a determination method for determining the reliability of the selective binding amount of the substance to be inspected wherein a position of the spot is determined in image data obtained by imaging the detection intensity on the carrier, and a pixel group corresponding to the spot
- a pixel group extracting step for extracting a pixel group, a median value of the detected intensity in the pixel group extracted in the pixel group extracting step, and a group obtained by excluding pixels having a higher predetermined ratio and / or a lower predetermined ratio from the pixel group
- a median calculating step for calculating a ratio or difference between the detected median and the median value of the detected intensity, and the median calculating step
- the ratio or the difference calculated Te based on a predetermined reference value, characterized in that it comprises a and a reliability determining step of determining the quality of the reliability.
- the median calculation step calculates a value of a ratio obtained by the following Equation 1 and / or Equation 2, and the reliability determination step includes the median If the value of the ratio calculated in the calculation step is equal to or greater than the reference value, it is determined as defective.
- / X (Formula 2) (Where X is the median value of the detected intensity in the pixel group, Xt is the median value of the detected intensity in the group obtained by removing the upper predetermined ratio of pixels from the pixel group, and Xb is the median value of the detected intensity. (This is the median value of the detection intensity in the group obtained by excluding the lower predetermined ratio of pixels from the pixel group.)
- the determination method of the present invention is characterized in that, in the determination method described above, the reference value is a value obtained by the following expression 3.
- S C + Z / X (Formula 3) (Here, S is the reference value, C is a constant, Z is an offset value according to the sensitivity setting of the device that detects the detection intensity of the marker, and X is the detection value in the pixel group) The median value of intensity.)
- the device for detecting the detection intensity of the marker is a photomultiplier
- the offset value is a value obtained by the following expression 4.
- Z X ⁇ (A) * B (Formula 4) (Here, Z is the offset value, X is the gain voltage of the photomultiplier, and A and B are constants.)
- the determination method of the present invention is the determination method described above, wherein the carrier is a microarray, the label is a fluorescent label, the detection intensity is a fluorescence amount, and the reliability determination step includes: The quality of the spot is determined as the quality of the reliability.
- the determination apparatus of the present invention provides a selective detection of the substance to be inspected, which is obtained as the detection intensity of the label by binding the substance to be inspected to the selective binding substance immobilized as a spot on the carrier.
- the control unit determines the position of the spot in the image data obtained by imaging the detection intensity on the carrier and corresponds to the spot
- a pixel group extracting means for extracting a pixel group to be detected, a median value of the detected intensity in the pixel group extracted by the pixel group extracting means, and pixels having a predetermined upper ratio and / or a lower predetermined ratio from the pixel group
- Median calculating means for calculating a ratio or difference between the detected median value and the median value of the detected intensity, and the median calculating method.
- the ratio or the difference calculated by, based on a predetermined reference value, characterized in that and a reliability determining means for determining the acceptability of the reliability.
- the determination system of the present invention includes a detection device that reads the detection intensity of a label obtained by binding a labeled test substance to a selective binding substance immobilized as a spot on a carrier, and the detection intensity described above.
- a determination apparatus comprising at least a control unit that determines the reliability of the selective binding amount of the substance to be inspected, wherein the control unit of the determination device includes the detection Image data acquisition means for acquiring the detected intensity in the carrier read through the apparatus as imaged image data, and determining the position of the spot in the image data acquired by the image data acquisition means.
- Pixel group extracting means for extracting a pixel group corresponding to the spot, and the pixel extracted by the pixel group extracting means
- a median calculating means for calculating a ratio or difference between the median value of the detected intensity in the pixel group and the median value of the detected intensity in a group obtained by removing pixels of the upper predetermined ratio and / or the lower predetermined ratio from the pixel group;
- a reliability determination unit that determines the quality of the reliability based on the ratio or the difference calculated by the median calculation unit and a predetermined reference value.
- the program of the present invention provides the selective binding of the test substance obtained as the detection intensity of the label by binding the labeled test substance to the selective binding substance immobilized as a spot on the carrier.
- a pixel group extraction step for determining and extracting a pixel group corresponding to the spot; a median value of the detection intensity in the pixel group extracted in the pixel group extraction step; And / or the ratio of the detected intensity to the median value in the group excluding the lower predetermined ratio of pixels
- a median calculating step for calculating a difference; and a reliability determining step for determining whether the reliability is good or not based on the ratio or the difference calculated in the median calculating step and a predetermined reference value.
- the present invention also relates to a recording medium, and is characterized by recording the above-described program.
- the reliability of the selective binding amount of the test substance is appropriately considered in consideration of non-biological effects. There is an effect that it can be evaluated.
- FIG. 1 is a scatter plot showing the detection intensity when the same test substance is hybridized on two DNA chips.
- FIG. 2 is a diagram showing an intensity image (left figure) of the spot 1 and a fluctuation graph (right figure) of the detected intensity along the vertical line of the spot 1.
- FIG. 3 is a diagram showing an intensity image (left figure) of the spot 2 and a fluctuation graph (right figure) of the detected intensity along the vertical line of the spot 2.
- FIG. 4 is a block diagram illustrating an example of the overall configuration of the determination system according to the present embodiment.
- FIG. 5 is a flowchart showing an example of basic processing of the determination apparatus 100 in the present embodiment.
- FIG. 5 is a flowchart showing an example of basic processing of the determination apparatus 100 in the present embodiment.
- FIG. 6 is a diagram showing the relationship between the spot median (X), the median (Xt: Top Cut Median) when the upper predetermined ratio is excluded, and the median (Xb: Bottom Cut Median) when the lower predetermined ratio is excluded. It is.
- FIG. 7 is a diagram showing white noise and a correction curve in a region where the signal intensity is weak.
- FIG. 8 is a diagram showing the results of measuring the same carrier provided with the same substance to be inspected when the PMT (photomultiplier tube) gain voltage (%) is set to “40” (40% ⁇ 1 V). is there.
- FIG. 9 is a diagram showing the results of measuring the same carrier provided with the same substance to be inspected when the setting of the PMT gain voltage (%) is “55” (55% ⁇ 1 V).
- FIG. 10 is a graph plotting the offset value and the PMT gain voltage (sensitized control voltage).
- FIG. 11 is a scatter plot showing the detection intensity when the same test substance is hybridized in two DNA chips.
- FIG. 12 is a diagram showing an intensity image (upper figure) of a spot having a non-uniform inner side and a fluctuation graph (lower figure) of detected intensity at a spot crossing line.
- FIG. 13 is a diagram illustrating a result of sorting the values of Equation 1 obtained for each spot in ascending order.
- FIG. 14 is a diagram illustrating a result of sorting the values of Equation 2 obtained for each spot in descending order.
- FIG. 15 is a scatter plot showing the results when the PMT gain voltage setting is 40%.
- FIG. 16 is a scatter plot showing the results when the PMT gain voltage setting is 55%.
- FIG. 17 is a scatter plot showing the results when the PMT gain voltage setting is 70%.
- FIG. 18 is a log-log graph showing the correlation between the PMT gain voltage setting value and the offset value.
- the determination device is connected to an input device such as a measurement device that acquires image data or the like, but the present invention is not limited to this, and the determination device is included in the input device.
- the image data may be stored in the storage unit in advance without being connected, or the image data may be acquired from the outside by communication.
- this embodiment schematically has the following basic features.
- the reliability of the selective binding amount of the test substance obtained as the detection intensity of the label by binding the test substance labeled to the selective binding substance immobilized as a spot on the carrier. Determine.
- the “substance to be inspected” refers to a sample obtained directly or indirectly from a cell or tissue, for example, genomic DNA, RNA, cDNA, aRNA (amplified using cDNA or its complementary sequence as a template). RNA), protein, sugar chain, lipid and the like.
- the “label” refers to a substance that can be detected by a detection means, such as a fluorescent label, a bioluminescent label, or a radioisotope label.
- the “selective binding substance” refers to a substance that selectively binds to a substance, such as complementary DNA to DNA, complementary RNA to DNA, an antibody to an antigen, an enzyme to a chemical substance, and the like. .
- the “carrier” is not only a general substrate such as a DNA chip or microarray, but also a DNA chip substrate made of polymethyl methacrylate having a concavo-convex structure (3D-Gene (trade name), manufactured and sold by Toray Industries, Inc. (company name) JP, 2004-264289, A) may be sufficient, and in that case, a spot is performed on a convex part.
- a DNA chip substrate made of polymethyl methacrylate having a concavo-convex structure (3D-Gene (trade name), manufactured and sold by Toray Industries, Inc. (company name) JP, 2004-264289, A) may be sufficient, and in that case, a spot is performed on a convex part.
- the position of a spot is determined in image data obtained by imaging the detection intensity on the carrier, and a pixel group corresponding to the spot is extracted. That is, a group of pixels constituting one spot location on the image is extracted for each spot.
- Xb) and the ratio or difference is calculated.
- the ratio may be calculated based on the following Expression 1 and / or Expression 2.
- the upper predetermined ratio and the lower predetermined ratio excluded are 30% of the extracted pixel group, but are not limited thereto. Any value in the range of ⁇ 40% may be used.
- this embodiment determines the quality of reliability based on the calculated ratio or difference and a predetermined reference value S. For example, in the present embodiment, if the calculated value of the ratio or difference is equal to or greater than the reference value S, it may be determined as defective.
- the value (%) of the formula 1 was 34%, and the value (%) of the formula 2 was 22.2%.
- the value (%) of Formula 1 was 30%, and the value (%) of Formula 2 was 22.3%. Since sufficient data is obtained in both spots 1 and 2, the reference value S may be about 25% to 30% in order not to exclude these data.
- the reference value S may be a value obtained by Equation 3 below.
- S C + Z / X (Formula 3) (Here, S is a reference value, C is a constant, Z is an offset value according to the sensitivity setting of the device for detecting the detection intensity of the label, and X is the median value of the detection intensity in the pixel group) .)
- the offset value Z may be a value obtained by Expression 4 below.
- Z X ⁇ (A) * B (Formula 4) (Here, Z is an offset value, X is a gain voltage of the photomultiplier, and A and B are constants.)
- FIG. 4 is a block diagram illustrating an example of the overall configuration of the determination system according to the present embodiment.
- the determination system generally includes an input device 112 that functions as detection means for reading the detection intensity of the sign, an output device 114, and a determination device 100.
- the determination apparatus 100 generally includes communication control connected to a control unit 102 such as a CPU that controls the entire determination apparatus 100 in an integrated manner, and a communication apparatus (not shown) such as a router connected to a communication line or the like.
- the interface unit 104, the input / output control interface unit 108 connected to the input device 112 and the output device 114, and the storage unit 106 for storing various databases and tables are configured. It is connected so that it can communicate via a road.
- the determination device 100 may be connected to the network 300 via a communication device such as a router and a wired or wireless communication line such as a dedicated line, and further via the network 300. It may be connected to the external system 200.
- the various databases and tables (image data file 106a and pixel group file 106b) stored in the storage unit 106 may be storage means such as a fixed disk device.
- the storage unit 106 may store various programs, tables, files, databases, web pages, and the like used for various processes.
- the image data file 106a is an image data storage unit that stores image data obtained by imaging the detected intensity on the carrier.
- the image data stored in the image data file 106a may be image information obtained by scanning the detection intensity on the carrier plane by the detection means.
- the image data file 106 a may store image data in advance, or may be image data input from the input device 112 that functions as a detection unit described later, and is received from the external system 200 via the network 300.
- the image data may be the same.
- the image data file 106a stores, as image data, grayscale image data having a detection intensity as a gradation value.
- the pixel group file 106b is pixel group storage means for storing pixel group information (intensity information and the like) for each spot.
- the information stored in the pixel group file 106b is a detection intensity value (tone value) associated with identification information (carrier block number, row number, column number, etc.) that uniquely identifies a spot. Or the median value of a pixel group.
- the pixel group file 106b is individually The position (block number and the row number and column number of the spot) where the selective binding substance (DNA fragment or the like) corresponding to the gene is arranged in association with the gradation value of the pixel group in the spot is stored. .
- the input / output control interface unit 108 controls the input device 112 and the output device 114.
- the output device 114 a monitor (including a home television), a printer, a recording medium output device, and the like can be used.
- the input device 112 may be a detection device (detection means) that reads the detection intensity of the sign in addition to a keyboard, a mouse, and the like.
- the “detection means” is means for reading the detection intensity of the label obtained by binding the labeled test substance to the selective binding substance immobilized as a spot on the carrier.
- the detection means may be an inspection means for identifying the position where the selective binding substance is spotted and acquiring the detection intensity, such as a fluorescence microscope camera.
- the detection means may be a photomultiplier (photomultiplier tube).
- the detection means is not limited to a means for imaging the detection intensity as in a fluorescence microscope camera, but may be any means that reads the detection intensity, and by scanning along the carrier plane, the determination apparatus 100 side You may image by.
- the selective binding substance is DNA
- a small amount of double-stranded DNA binding fluorescent substance may be incorporated to detect the selective binding amount.
- the detection means may detect an absorption wavelength unique to DNA. Even if the selective binding substance is a protein, sugar chain, etc., depending on the properties of the selective binding substance, absorption wavelength, fluorescent substance, radioisotope, hybridization, antigen-antibody reaction, etc. You may make it detect using a technique.
- control unit 102 has a control program such as an OS (Operating System), a program defining various processing procedures, and an internal memory for storing necessary data. Information processing for executing various processes is performed.
- the control unit 102 is functionally conceptually configured to include an image data acquisition unit 102a, a pixel group extraction unit 102b, a median calculation unit 102c, a reliability determination unit 102d, and a reference value determination unit 102e.
- the image data acquisition unit 102a is image data acquisition means for acquiring image data obtained by imaging the detected intensity in the carrier.
- the image data acquisition unit 102a may acquire the detected intensity of the carrier read via the input device 112 functioning as a detection unit as imaged image data.
- the image data acquisition unit 102a may directly acquire image data from the input device 112 such as a fluorescence microscope camera, and by scanning the input device 112 such as a photomultiplier along the carrier plane. The detected intensity for each obtained coordinate may be imaged.
- the image data acquisition unit 102 a may receive image data from the external system 200 via the network 300.
- the image data acquisition unit 102a stores the acquired image data in the image data file 106a.
- the pixel group extraction unit 102b is a pixel group extraction unit that determines a spot position in the image data stored in the image data file 106a and extracts a pixel group corresponding to the spot. For example, the pixel group extraction unit 102b partitions each spot location on the image based on the arrangement (spot center coordinates, pixel radius, etc.) of the position (spot location) where the selective binding substance on the carrier is arranged, You may extract the pixel group in each division.
- the pixel group extraction unit 102b displays the image data and the array pattern data on the output device 114 in a superimposed manner, and causes the user to perform an input for moving the array pattern on the display via the input device 112 such as a mouse. Depending on the situation, positioning may be performed.
- the pixel group extraction unit 102b stores pixel group information (intensity information and the like) for each spot in the pixel group file 106b. For example, the pixel group extraction unit 102b associates the detection intensity value (such as the gradation value) of each pixel with the identification information (such as the carrier block number, row number, and column number) that uniquely identifies the spot. You may store in the group file 106b.
- the median calculation unit 102c detects the median value (X) of the detection intensity in the pixel group stored in the pixel group file 106b and the detection in the group obtained by removing the upper predetermined ratio and / or the lower predetermined ratio from the pixel group. It is a median calculating means for calculating a ratio or difference between the intensity and the median value (Xt / Xb).
- the ratio may be calculated based on the following Expression 1 and / or Expression 2.
- the median calculation unit 102c rearranges the detection intensities (gradation values, etc.) of the pixel group of a certain spot stored in the pixel group file 106b by sorting in ascending or descending order, and the upper predetermined ratio (x% ) And the medium detection intensity in the group excluding the lower predetermined ratio (y%) may be obtained.
- / X (Formula 2) (Where X is the median value of the detected intensity in the extracted pixel group, Xt is the median value of the detected intensity in the group obtained by removing the upper predetermined ratio of pixels from the pixel group, and Xb is the pixel in question. (This is the median value of the detection intensity in the group obtained by removing the lower predetermined ratio of pixels from the group.)
- the reliability determination unit 102d is a reliability determination unit that determines whether the reliability is good or not based on the ratio or difference calculated by the median calculation unit 102c and a predetermined reference value. For example, the reliability determination unit 102d may determine that the ratio is not satisfactory if the ratio or difference calculated by the median calculation unit 102c is greater than or equal to the reference value S. The reliability determination unit 102d may exclude defective spot data from the analysis based on the reliability determination result, and may output the reliability determination result to the output device 114. The control unit 102 may perform processing (exclusion of defective spot data, etc.) based on the reliability determination result by the reliability determination unit 102d and output the data analysis result to the output device 114.
- the output is not limited to being output to the monitor as the output device 114, but may be output to a recording medium or the like via a printer or a recording medium output device.
- the control unit 102 such as the reliability determination unit 102 d may control the communication control interface 104 and transmit the reliability determination result and the analysis result data to the external system 200 via the network 300.
- the reference value determination unit 102e is a reference value determination unit that determines a reference value that is used as a reference for reliability determination by the reliability determination unit 102d.
- the reference value determination unit 102e may determine the reference value S according to the following Expression 3.
- S C + Z / X (Formula 3) (Here, S is a reference value, C is a constant, Z is an offset value according to the sensitivity setting of the device for detecting the detection intensity of the label, and X is the median value of the detection intensity in the pixel group) .
- the reference value determination unit 102e may determine the offset value Z in Expression 3 based on Expression 4 below.
- Z X ⁇ (A) * B (Formula 4) (Here, Z is an offset value, X is a gain voltage of the photomultiplier, and A and B are constants.)
- the communication control interface unit 104 performs communication control between the determination device 100 and the network 300 (or a communication device such as a router). That is, the communication control interface unit 104 has a function of communicating data with other terminals via a communication line.
- the network 300 has a function of connecting the determination apparatus 100 and the external system 200 to each other, such as the Internet.
- the external system 200 is connected to the determination apparatus 100 via the network 300, and has a function of providing an external database, an external program, and the like regarding the detection intensity data of the sign.
- the external system 200 may be configured as a server device such as a WEB server or an ASP server, or a terminal device.
- the hardware configuration of the external system 200 is an information processing device such as a commercially available workstation or personal computer. You may comprise by the attached apparatus.
- Each function of the external system 200 is realized by a CPU, a disk device, a memory device, an input device, an output device, a communication control device, and the like in the hardware configuration of the external system 200 and a program for controlling them.
- the user of the determination apparatus 100 can detect a detection intensity database related to a DNA chip or the like provided by the external system 200, an external database such as a gene arrangement database, a program for executing a determination method, or the like via the network 300. You may make it acquire detection intensity value data, arrangement
- FIG. 5 is a flowchart showing an example of basic processing of the determination apparatus 100 in the present embodiment.
- the image data acquisition unit 102a acquires image data obtained by imaging the detected intensity on the carrier and stores it in the image data file 106a (step SA-1).
- the image data acquisition unit 102a may acquire the detected intensity of the carrier read via the input device 112 functioning as a detection unit as imaged image data.
- the image data acquisition unit 102a may directly acquire image data from the input device 112 such as a fluorescence microscope camera, and by scanning the input device 112 such as a photomultiplier along the carrier plane. The detected intensity for each obtained coordinate may be imaged.
- the image data acquisition unit 102a may read data from the external recording medium storing the detected intensity data or the image data via the input / output control interface unit 108.
- the pixel group extraction unit 102b determines the position of the spot in the image data stored in the image data file 106a, extracts the pixel group corresponding to the spot, and stores the extracted pixel group information in the pixel group file 106b.
- Store step SA-22).
- the pixel group extraction unit 102b partitions each spot location on the image based on the arrangement (spot center coordinates, pixel radius, etc.) of the position (spot location) where the selective binding substance on the carrier is arranged, You may extract the pixel group in each division.
- the pixel group extraction unit 102b stores pixel group information (intensity information and the like) for each spot in the pixel group file 106b.
- the pixel group extraction unit 102b associates the detection intensity value (such as the gradation value) of each pixel with the identification information (such as the carrier block number, row number, and column number) that uniquely identifies the spot. You may store in the group file 106b.
- the median calculation unit 102c acquires the median value (X) of the detection intensity in the pixel group stored in the pixel group file 106b (Step SA-3). For example, the median calculation unit 102c rearranges the detection intensities (gradation values, etc.) of the pixel group of a certain spot by sorting in ascending or descending order, and acquires the detection intensity value at the middle position (center position in the ranking). By doing so, the median value (X) may be obtained.
- the median calculation unit 102c detects the median value (Xt) of the detection intensity in the group obtained by removing the upper predetermined ratio (x%) of pixels from the pixel group stored in the pixel group file 106b and / or the lower predetermined value from the pixel group.
- the median value (Xb) of the detection intensity in the group excluding the pixels of the ratio (y%) is acquired (step SA-4).
- the upper predetermined ratio (x%) and the lower predetermined ratio (y%) may be, for example, 30%.
- the median calculation unit 102c rearranges the detection intensities (gradation values, etc.) of a pixel group of a certain spot by sorting in ascending or descending order, and removes the upper predetermined ratio (x%) Detection intensity (Xt) and the medium detection intensity (Xb) in the group excluding the lower predetermined ratio (y%) may be obtained.
- FIG. 6 shows the relationship between the spot median (X), the median (Xt: Top Cut Median) when the upper predetermined ratio is removed, and the median (Xb: Bottom Cut Median) when the lower predetermined ratio is excluded.
- FIG. One bar graph corresponds to one pixel, and the length is a detection intensity corresponding to the gradation value of the pixel.
- the pixel groups in the spot are sorted in descending order by intensity.
- the spot median (X) is an intermediate value (center position) in the rank rank of the pixel group.
- the median Xt Top Cut Median
- the median Xb Bottom Cut Median
- the median Xb Bottom Cut Median
- the median calculation unit 102c determines that the ratio or difference between the spot median X and the median Xt excluding the upper predetermined ratio and / or the ratio between the spot median X and the median Xb excluding the lower predetermined ratio or The difference is calculated (step SA-5). 6 represents the difference between the spot median X and the median Xt excluding the upper predetermined ratio, and the difference between the spot median X and the median Xb excluding the lower predetermined ratio.
- the median calculation unit 102c may calculate the ratio based on the following Expression 1 and / or Expression 2.
- the reliability determination unit 102d determines whether the reliability is good or not based on the ratio or difference calculated by the median calculation unit 102c and a predetermined reference value (step SA-6). For example, the reliability determination unit 102d may determine that the ratio is not good if the ratio or difference (such as the absolute value of the difference) calculated by the median calculation unit 102c is equal to or greater than a predetermined reference value S (%).
- the reference value S (%) may be a constant value C (%), and when white noise from an electric circuit system cannot be ignored in a region where the signal intensity is weak, the reference C (%) is used as a reference.
- the value may be corrected by the value determining unit 102e (details will be described later).
- the above is an example of the basic processing of the determination apparatus 100.
- FIG. 7 is a diagram showing white noise and a correction curve in a region where the signal intensity is weak.
- the horizontal axis is the spot median X representing signal intensity
- the vertical axis is the ratio value obtained by Equation 1 and Equation 2.
- the values obtained by Equation 2 are plotted for the positive region on the vertical axis, and the values obtained by Equation 1 are plotted for the negative region. Yes.
- the reference value determination unit 102e may determine the reference value S by correcting the basic C (%) based on the following Equation 3.
- the offset value Z is constant in the same setting, but is a value that changes when the gain setting of the photomultiplier at the time of scanning is changed.
- the spraying state also changes, so that the threshold curve needs to be changed accordingly.
- FIG. 10 is a graph plotting the offset value and the gain voltage (sensitized control voltage).
- the horizontal axis represents the logarithmic value of the gain voltage of the photomultiplier (PMT), and the vertical axis represents the logarithmic value of the offset value.
- the offset value changes depending on the gain voltage, and is independent of the types of the signs (cy3 and cy5).
- the logarithmic value of the offset value and the logarithmic value of the gain voltage are expressed in a linear format (the expression in FIG. 10).
- the inclination and intercept were values specific to the scanner body (photomultiplier).
- the reference value determination unit 102e may determine the offset value Z based on the following Equation 4.
- Z X ⁇ (A) * B (Formula 4)
- Z is an offset value
- X is a gain voltage of the photomultiplier
- a and B are constants.
- the slope of the graph is about 5.1639 to 5.12522, and the intercept is about ⁇ 7.1363 to ⁇ 7.2711. That is, the range of A is about ⁇ 7.1363 to ⁇ 7.27111, and the range of B is about 5.1639 to 5.12522.
- Example 1 in which an appropriate threshold value is studied will be described below. First, the same test substance was hybridized in two DNA chips, and it was confirmed that there was no big difference in detection intensity.
- FIG. 11 is a scatter plot showing the detection intensity when the same test substance is hybridized in two DNA chips. The horizontal axis represents the detection intensity of each spot on one DNA chip, and the vertical axis represents the detection intensity of each spot on the other DNA chip.
- the following formulas 1 and 2 are used for a spot having a detected intensity of gain voltage setting value “40” (the photomultiplier tube control voltage is 40% ⁇ 1 V) and a spot median of 2000 or more. The value of was calculated.
- the reason why the spot median of 2000 or more was examined is that, as described above, white noise such as an electric circuit cannot be ignored in a portion where the signal is weak (see, for example, FIGS. 7 and 8).
- X is the median value of the detection intensity in the pixel group
- Xt is the median value of the detection intensity in the group obtained by removing the upper predetermined ratio of pixels from the pixel group
- Xb is from the pixel group. (This is the median value of the detection intensity in the group excluding the lower predetermined ratio of pixels.)
- FIG. 13 is a diagram showing a result of sorting the values of Formula 1 obtained for each spot in ascending order.
- FIG. 14 is a diagram illustrating a result of sorting the values of Equation 2 obtained for each spot in descending order. Note that Block represents a block number, Column represents a row number, and Row represents a column number, and a spot is uniquely identified by these. In the item of S_532_Median, a spot median is indicated.
- the reference value for judging the reliability is preferably about 25% in absolute value.
- Example 2 A second embodiment in which an appropriate offset value is studied will be described below.
- Example 2 two chips of 3D-Gene (registered trademark) Human Ver1.1 (manufactured by Toray) were used as the carrier.
- Stratagene's Human Reference RNA was used as a test substance, and hybridization was performed as described in the protocol.
- 3D-Gene (registered trademark) Scanner manufactured by Toray Industries, Inc. was used as a detection means.
- the 3D-Gene (registered trademark) Scanner (manufactured by Toray) photomultiplier (PMT) setting was changed to 40%, 55%, and 70%.
- the values of Equation 1 and Equation 2 were calculated for each spot.
- the control voltage of the PMT changes in proportion. Therefore, the photomultiplier tube gain also changes.
- FIG. 15 is a scatter plot showing the result when the PMT setting is 40%
- FIG. 16 is a scatter plot showing the result when the PMT setting is 55%
- FIG. It is a scatter plot which shows a result in case of setting 70%.
- the cut value was 30%, and the horizontal axis plotted with the spot median and the vertical axis plotted with two values according to Formula 1 and Formula 2 excluding the absolute value calculation.
- Equation 4 for offset value correction is shown below.
- Z X ⁇ (A) * B (Formula 4)
- Z is an offset value
- X is a gain voltage of the photomultiplier.
- A was ⁇ 7.13 and B was 5.16.
- the reliability of the selective binding amount of the test substance is appropriately evaluated in consideration of the non-biological influence.
- the embodiment of the present invention has been described with respect to an example in which DNA is used as a selective binding substance.
- the selective binding substance includes a protein library such as an antibody, a library of compounds, and the like. You may arrange.
- the material of the “carrier” is not limited to glass, and may be a membrane or plastic.
- a fluorescent chemical substance eg, Cy3, Cy5
- the present invention is not limited to this, and the label includes a dye having no fluorescence property, a radioisotope, GFP / GRP. Proteins such as, His tag, biotinylation and the like can be used.
- processing is performed in response to a request from a client terminal configured with a separate casing from the determination apparatus 100, and the processing result is displayed. You may comprise so that it may return to the said client terminal.
- each illustrated component is functionally schematic and does not necessarily need to be physically configured as illustrated.
- each device of the determination device 100 is interpreted and executed by a CPU (Central Processing Unit) and the CPU. It may be realized by a program or hardware based on wired logic. The program is recorded on a recording medium to be described later, and is mechanically read by the determination apparatus 100 as necessary. That is, the storage unit 106 such as ROM or HD stores a computer program for performing various processes by giving instructions to the CPU in cooperation with an OS (Operating System). This computer program is executed by being loaded into the RAM, and constitutes a control device in cooperation with the CPU.
- OS Operating System
- the computer program may be stored in an application program server connected to the determination apparatus 100 via an arbitrary network 300, and may be downloaded in whole or in part as necessary. is there.
- the program according to the present invention can be stored in a computer-readable recording medium.
- the “recording medium” means a memory card, USB memory, SD card, flexible disk, magneto-optical disk, ROM, EPROM, EEPROM, CD-ROM, MO, DVD, and Blu-ray (registered trademark). It includes any “portable physical medium” such as Disc.
- program is a data processing method described in an arbitrary language or description method, and may be in any form such as source code or binary code.
- program is not necessarily limited to a single configuration, but is distributed in the form of a plurality of modules and libraries, or in cooperation with a separate program typified by an OS (Operating System). Including those that achieve the function.
- OS Operating System
- a well-known configuration and procedure can be used for a specific configuration for reading a recording medium, a reading procedure, an installation procedure after reading, and the like in each device described in the embodiment.
- Various databases and the like (image data file 106a to pixel group file 106b) stored in the storage unit 106 are storage means such as a memory device such as a RAM and a ROM, a fixed disk device such as a hard disk, a flexible disk, and an optical disk. Yes, it stores various programs, tables, databases, web page files, etc. used for various processing and website provision.
- the determination device 100 may be configured as an information processing device such as a known personal computer or workstation, or may be configured by connecting an arbitrary peripheral device to the information processing device.
- the determination apparatus 100 may be realized by installing software (including a program, data, and the like) that causes the information processing apparatus to realize the method of the present invention.
- the specific form of the distribution / integration of the devices is not limited to that shown in the figure, and all or a part of them may be functionally or physically in arbitrary units according to various additions or according to functional loads. Can be distributed and integrated. That is, the above-described embodiments may be arbitrarily combined and may be selectively implemented.
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Immunology (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Hematology (AREA)
- Urology & Nephrology (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Cell Biology (AREA)
- Microbiology (AREA)
- Pathology (AREA)
- Tropical Medicine & Parasitology (AREA)
- Analytical Chemistry (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Biochemistry (AREA)
- Biotechnology (AREA)
- Multimedia (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Apparatus Associated With Microorganisms And Enzymes (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
Description
CV値=(ピクセル群における強度の標準偏差)/(ピクセル群における強度の平均値)
|X-Xt|/X …(式1)
|X-Xb|/X …(式2)
(ここで、Xは、上記ピクセル群における上記検出強度のメジアン値であり、Xtは、当該ピクセル群から上位所定割合のピクセルを除いた群における上記検出強度のメジアン値であり、Xbは、当該ピクセル群から下位所定割合のピクセルを除いた群における上記検出強度のメジアン値である。)
S=C+Z/X …(式3)
(ここで、Sは、上記基準値であり、Cは、定数であり、Zは、上記標識の検出強度を検出する装置の感度設定によるオフセット値であり、Xは、上記ピクセル群における上記検出強度のメジアン値である。)
Z=X^(A)*B …(式4)
(ここで、Zは、上記オフセット値であり、Xは、上記フォトマルチプライヤーのゲイン電圧であり、AおよびBは、定数である。)
以下、本実施形態の概要について説明し、その後、本実施形態の構成および処理等について図面を参照しながら詳細に説明する。本願発明者らは、鋭意検討した結果、一例として以下に示す本実施形態による判定方法を開発するに至った。
|X-Xt|/X …(式1)
|X-Xb|/X …(式2)
(ここで、Xは、抽出したピクセル群における検出強度のメジアン値であり、Xtは、当該ピクセル群から上位所定割合のピクセルを除いた群における検出強度のメジアン値であり、Xbは、当該ピクセル群から下位所定割合のピクセルを除いた群における検出強度のメジアン値である。)
S=C+Z/X …(式3)
(ここで、Sは、基準値であり、Cは、定数であり、Zは、標識の検出強度を検出する装置の感度設定によるオフセット値であり、Xは、ピクセル群における検出強度のメジアン値である。)
Z=X^(A)*B …(式4)
(ここで、Zは、オフセット値であり、Xは、フォトマルチプライヤーのゲイン電圧であり、AおよびBは、定数である。)
まず、本実施形態における判定システムの構成について説明する。図4は、本実施形態における判定システムの全体構成の一例を示すブロック図である。
|X-Xt|/X …(式1)
|X-Xb|/X …(式2)
(ここで、Xは、抽出したピクセル群における検出強度のメジアン値であり、Xtは、当該ピクセル群から上位所定割合のピクセルを除いた群における検出強度のメジアン値であり、Xbは、当該ピクセル群から下位所定割合のピクセルを除いた群における検出強度のメジアン値である。)
S=C+Z/X …(式3)
(ここで、Sは、基準値であり、Cは、定数であり、Zは、標識の検出強度を検出する装置の感度設定によるオフセット値であり、Xは、ピクセル群における検出強度のメジアン値である。)
Z=X^(A)*B …(式4)
(ここで、Zは、オフセット値であり、Xは、フォトマルチプライヤーのゲイン電圧であり、AおよびBは、定数である。)
次に、このように構成された本実施形態における判定装置100の処理の一例について、以下に図5および図6を参照しながら詳細に説明する。図5は、本実施の形態における本判定装置100の基本処理の一例を示すフローチャートである。
|X-Xt|/X …(式1)
|X-Xb|/X …(式2)
(ここで、Xは、抽出したピクセル群における検出強度のメジアン値であり、Xtは、当該ピクセル群から上位所定割合のピクセルを除いた群における検出強度のメジアン値であり、Xbは、当該ピクセル群から下位所定割合のピクセルを除いた群における検出強度のメジアン値である。)
上述の判定装置100の基本処理においては、基準値Sは、一定値Cとしたが、基準値決定部102eが、基準値を補正して決定する補正処理の一例について以下に説明する。ここで、図7は、シグナル強度が弱い領域におけるホワイトノイズと補正曲線を示す図である。図7において、横軸は、シグナル強度を代表するスポットメジアンXであり、縦軸は、式1および式2によって得られた比の値である。なお、この例では、式1および式2において絶対値を求めていないので、縦軸の正の領域は、式2によって求めた値、負の領域は、式1によって求めた値がプロットされている。
S=C+Z/X …(式3)
(ここで、Sは、基準値であり、Cは、定数であり、Zは、標識の検出強度を検出する装置の感度設定によるオフセット値であり、Xは、ピクセル群における検出強度のメジアン値である。)
Z=X^(A)*B …(式4)
ここで、Zは、オフセット値であり、Xは、フォトマルチプライヤーのゲイン電圧であり、AおよびBは、定数である。図10の例では、グラフの傾きは、5.1639~5.12552程度であり、切片は、-7.1363~-7.2711程度である。すなわち、Aの範囲は-7.1363~-7.2711程度であり、Bの範囲は5.1639~5.12552程度である。
適切な閾値についての検討を行った実施例1について以下に説明する。まず、2枚のDNAチップにおいて同一の被検査物質をハイブリダイズさせ、検出強度に大きな差異がないことを確かめた。図11は、2枚のDNAチップにおいて同一の被検査物質をハイブリダイズさせた場合の検出強度を示すスキャッタープロットである。横軸は、一方のDNAチップにおける各スポットの検出強度を表しており、縦軸は、他方のDNAチップにおける各スポットの検出強度を表している。
(X-Xt)/X*100(%) …(式1)
(X-Xb)/X*100(%) …(式2)
(ここで、Xは、ピクセル群における検出強度のメジアン値であり、Xtは、当該ピクセル群から上位所定割合のピクセルを除いた群における検出強度のメジアン値であり、Xbは、当該ピクセル群から下位所定割合のピクセルを除いた群における検出強度のメジアン値である。)
適切なオフセット値についての検討を行った実施例2について以下に説明する。実施例2においては、担体として、3D-Gene(登録商標) Human Ver1.1(東レ製)のチップを2枚使用した。また、被検査物質として、Stratagene社製Human Reference RNAを用いて、プロトコルに記載のとおりにハイブリダイゼーションまでを行った。また、検出手段として、3D-Gene(登録商標) Scanner(東レ製)を用いた。
Z=X^(A)*B …(式4)
ここで、Zは、オフセット値であり、Xは、フォトマルチプライヤーのゲイン電圧である。本実施例2において、Aは、-7.13であり、Bは、5.16であった。
さて、これまで本発明の実施の形態について説明したが、本発明は、上述した実施の形態以外にも、上記特許請求の範囲に記載した技術的思想の範囲内において種々の異なる実施の形態にて実施されてよいものである。
102 制御部
102a 画像データ取得部
102b ピクセル群抽出部
102c メジアン算出部
102d 信頼性判定部
102e 基準値決定部
104 通信制御インターフェース部
106 記憶部
106a 画像データファイル
106b ピクセル群ファイル
108 入出力制御インターフェース部
112 入力装置
114 出力装置
200 外部システム
300 ネットワーク
Claims (8)
- 担体上にスポットとして固定化された選択結合性物質に、標識された被検査物質が結合することにより標識の検出強度として得られる、上記被検査物質の選択的結合量の信頼性を判定する判定方法であって、
上記担体における上記検出強度を画像化した画像データにおいて上記スポットの位置を決定して、当該スポットに対応するピクセル群を抽出するピクセル群抽出ステップと、
上記ピクセル群抽出ステップにて抽出された上記ピクセル群における上記検出強度のメジアン値と、当該ピクセル群から上位所定割合および/または下位所定割合のピクセルを除いた群における上記検出強度のメジアン値と、の比または差を算出するメジアン算出ステップと、
上記メジアン算出ステップにて算出された上記比または上記差と、所定の基準値とに基づいて、上記信頼性の良否を判定する信頼性判定ステップと、
を含むことを特徴とする、判定方法。 - 請求項1に記載の判定方法において、
上記メジアン算出ステップは、
以下の式1および/または式2によって得られる比の値を算出し、
|X-Xt|/X …(式1)
|X-Xb|/X …(式2)
(ここで、Xは、上記ピクセル群における上記検出強度のメジアン値であり、Xtは、当該ピクセル群から上位所定割合のピクセルを除いた群における上記検出強度のメジアン値であり、Xbは、当該ピクセル群から下位所定割合のピクセルを除いた群における上記検出強度のメジアン値である。)
上記信頼性判定ステップは、
上記メジアン算出ステップにて算出された上記比の値が上記基準値以上であれば不良と判定することを特徴とする、判定方法。 - 請求項1または2に記載の判定方法において、
上記基準値は、以下の式3によって得られる値であること
S=C+Z/X …(式3)
(ここで、Sは、上記基準値であり、Cは、定数であり、Zは、上記標識の検出強度を検出する装置の感度設定によるオフセット値であり、Xは、上記ピクセル群における上記検出強度のメジアン値である。)
を特徴とする、判定方法。 - 請求項3に記載の判定方法において、
上記標識の検出強度を検出する装置は、フォトマルチプライヤーであり、
上記オフセット値は、以下の式4によって得られる値であること
Z=X^(A)*B …(式4)
(ここで、Zは、上記オフセット値であり、Xは、上記フォトマルチプライヤーのゲイン電圧であり、AおよびBは、定数である。)
を特徴とする、判定方法。 - 請求項1乃至4のいずれか一つに記載の判定方法において、
上記担体は、マイクロアレイであり、
上記標識は、蛍光標識であり、
上記検出強度は、蛍光量であり、
上記信頼性判定ステップは、上記信頼性の良否としてスポットの良否を判定すること
を特徴とする、判定方法。 - 担体上にスポットとして固定化された選択結合性物質に、標識された被検査物質が結合することにより標識の検出強度として得られる、上記被検査物質の選択的結合量の信頼性を判定する、制御部を少なくとも備えた判定装置において、
上記制御部は、
上記担体における上記検出強度を画像化した画像データにおいて上記スポットの位置を決定して、当該スポットに対応するピクセル群を抽出するピクセル群抽出手段と、
上記ピクセル群抽出手段により抽出された上記ピクセル群における上記検出強度のメジアン値と、当該ピクセル群から上位所定割合および/または下位所定割合のピクセルを除いた群における上記検出強度のメジアン値と、の比または差を算出するメジアン算出手段と、
上記メジアン算出手段により算出された上記比または上記差と、所定の基準値とに基づいて、上記信頼性の良否を判定する信頼性判定手段と、
を備えたことを特徴とする、判定装置。 - 担体上にスポットとして固定化された選択結合性物質に、標識された被検査物質が結合することにより得られる標識の検出強度を読取る検出装置と、上記検出強度として得られる、上記被検査物質の選択的結合量の信頼性を判定する制御部を少なくとも備えた判定装置と、を接続して構成した判定システムにおいて、
上記判定装置の上記制御部は、
上記検出装置を介して読取られた上記担体における上記検出強度を、画像化した画像データとして取得する画像データ取得手段と、
上記画像データ取得手段により取得された上記画像データにおいて上記スポットの位置を決定して、当該スポットに対応するピクセル群を抽出するピクセル群抽出手段と、
上記ピクセル群抽出手段により抽出された上記ピクセル群における上記検出強度のメジアン値と、当該ピクセル群から上位所定割合および/または下位所定割合のピクセルを除いた群における上記検出強度のメジアン値と、の比または差を算出するメジアン算出手段と、
上記メジアン算出手段により算出された上記比または上記差と、所定の基準値とに基づいて、上記信頼性の良否を判定する信頼性判定手段と、
を備えたことを特徴とする、判定システム。 - 担体上にスポットとして固定化された選択結合性物質に、標識された被検査物質が結合することにより標識の検出強度として得られる、上記被検査物質の選択的結合量の信頼性を判定する方法を、制御部を少なくとも備えたコンピュータに実行させるためのプログラムであって、
上記制御部において、
上記担体における上記検出強度を画像化した画像データにおいて上記スポットの位置を決定して、当該スポットに対応するピクセル群を抽出するピクセル群抽出ステップと、
上記ピクセル群抽出ステップにて抽出された上記ピクセル群における上記検出強度のメジアン値と、当該ピクセル群から上位所定割合および/または下位所定割合のピクセルを除いた群における上記検出強度のメジアン値と、の比または差を算出するメジアン算出ステップと、
上記メジアン算出ステップにて算出された上記比または上記差と、所定の基準値とに基づいて、上記信頼性の良否を判定する信頼性判定ステップと、
を含む上記方法を実行させるためのプログラム。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201380012447.4A CN104160265B (zh) | 2012-03-08 | 2013-03-05 | 判定方法、判定装置、判定系统及程序 |
KR1020147023805A KR102048599B1 (ko) | 2012-03-08 | 2013-03-05 | 판정 방법, 판정 장치, 판정 시스템 및 프로그램 |
BR112014021383-6A BR112014021383A2 (ja) | 2012-03-08 | 2013-03-05 | A judgment method, a determining device, a judgment system, and a program |
EP13758079.1A EP2824444B1 (en) | 2012-03-08 | 2013-03-05 | Determination method, determination device, determination system, and program |
US14/383,601 US9519822B2 (en) | 2012-03-08 | 2013-03-05 | Determination method, determination device, determination system, and computer program |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012052311A JP5895613B2 (ja) | 2012-03-08 | 2012-03-08 | 判定方法、判定装置、判定システム、および、プログラム |
JP2012-052311 | 2012-03-08 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013133283A1 true WO2013133283A1 (ja) | 2013-09-12 |
Family
ID=49116758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2013/056020 WO2013133283A1 (ja) | 2012-03-08 | 2013-03-05 | 判定方法、判定装置、判定システム、および、プログラム |
Country Status (7)
Country | Link |
---|---|
US (1) | US9519822B2 (ja) |
EP (1) | EP2824444B1 (ja) |
JP (1) | JP5895613B2 (ja) |
KR (1) | KR102048599B1 (ja) |
CN (1) | CN104160265B (ja) |
BR (1) | BR112014021383A2 (ja) |
WO (1) | WO2013133283A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016073218A (ja) * | 2014-10-03 | 2016-05-12 | ニプロ株式会社 | 判定装置、判定方法及びコンピュータプログラム |
WO2016152159A1 (ja) * | 2015-03-26 | 2016-09-29 | 東洋製罐グループホールディングス株式会社 | Dnaチップ画像のスポット有効性判定装置、dnaチップ画像のスポット有効性判定方法、及びdnaチップ画像のスポット有効性判定プログラム |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3023912B1 (fr) * | 2014-07-18 | 2020-05-08 | Thales | Calcul de performance pour aeronef |
CN110443787B (zh) * | 2019-07-30 | 2023-05-26 | 云谷(固安)科技有限公司 | 矫正装置和矫正方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004264289A (ja) | 2002-11-22 | 2004-09-24 | Toray Ind Inc | 選択結合性物質が固定化された基材 |
JP2004340574A (ja) | 2001-07-31 | 2004-12-02 | Toagosei Co Ltd | アレイにおけるスポットの均一性評価方法 |
JP2006084281A (ja) * | 2004-09-15 | 2006-03-30 | Hitachi Software Eng Co Ltd | スポット画像認識方法及びスポット画像認識システム |
JP2008039584A (ja) | 2006-08-07 | 2008-02-21 | Toray Ind Inc | 帯電防止性カバーを有するマイクロアレイ |
JP2009058356A (ja) * | 2007-08-31 | 2009-03-19 | Mitsubishi Rayon Co Ltd | 画像読み取り方法および装置 |
JP2009236626A (ja) * | 2008-03-26 | 2009-10-15 | Toyo Kohan Co Ltd | マイクロアレイの測定方法 |
JP2011004737A (ja) * | 2009-05-27 | 2011-01-13 | Fujifilm Corp | 核酸マイクロアレイの異常スポットを検出する方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6251601B1 (en) * | 1999-02-02 | 2001-06-26 | Vysis, Inc. | Simultaneous measurement of gene expression and genomic abnormalities using nucleic acid microarrays |
US7130458B2 (en) * | 2000-10-24 | 2006-10-31 | Affymetrix, Inc. | Computer software system, method, and product for scanned image alignment |
US7031844B2 (en) | 2002-03-18 | 2006-04-18 | The Board Of Regents Of The University Of Nebraska | Cluster analysis of genetic microarray images |
US7796804B2 (en) | 2007-07-20 | 2010-09-14 | Kla-Tencor Corp. | Methods for generating a standard reference die for use in a die to standard reference die inspection and methods for inspecting a wafer |
EP2347361B1 (en) * | 2008-10-13 | 2018-11-14 | Roche Diagnostics GmbH | Algorithms for classification of disease subtypes with gene expression profiling |
US20100304997A1 (en) * | 2009-05-27 | 2010-12-02 | Fujifilm Corporation | Method for detecting abnormal spots of nucleic acid microarray |
US20140307931A1 (en) * | 2013-04-15 | 2014-10-16 | Massachusetts Institute Of Technology | Fully automated system and method for image segmentation and quality control of protein microarrays |
-
2012
- 2012-03-08 JP JP2012052311A patent/JP5895613B2/ja active Active
-
2013
- 2013-03-05 US US14/383,601 patent/US9519822B2/en active Active
- 2013-03-05 EP EP13758079.1A patent/EP2824444B1/en active Active
- 2013-03-05 BR BR112014021383-6A patent/BR112014021383A2/ja not_active IP Right Cessation
- 2013-03-05 KR KR1020147023805A patent/KR102048599B1/ko active IP Right Grant
- 2013-03-05 CN CN201380012447.4A patent/CN104160265B/zh active Active
- 2013-03-05 WO PCT/JP2013/056020 patent/WO2013133283A1/ja active Application Filing
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004340574A (ja) | 2001-07-31 | 2004-12-02 | Toagosei Co Ltd | アレイにおけるスポットの均一性評価方法 |
JP2004264289A (ja) | 2002-11-22 | 2004-09-24 | Toray Ind Inc | 選択結合性物質が固定化された基材 |
JP2006084281A (ja) * | 2004-09-15 | 2006-03-30 | Hitachi Software Eng Co Ltd | スポット画像認識方法及びスポット画像認識システム |
JP2008039584A (ja) | 2006-08-07 | 2008-02-21 | Toray Ind Inc | 帯電防止性カバーを有するマイクロアレイ |
JP2009058356A (ja) * | 2007-08-31 | 2009-03-19 | Mitsubishi Rayon Co Ltd | 画像読み取り方法および装置 |
JP2009236626A (ja) * | 2008-03-26 | 2009-10-15 | Toyo Kohan Co Ltd | マイクロアレイの測定方法 |
JP2011004737A (ja) * | 2009-05-27 | 2011-01-13 | Fujifilm Corp | 核酸マイクロアレイの異常スポットを検出する方法 |
Non-Patent Citations (2)
Title |
---|
"GenePix Pro 7.0 Microarray Image Analysis", 2006, INTER MEDICAL CO., LTD. |
See also references of EP2824444A4 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016073218A (ja) * | 2014-10-03 | 2016-05-12 | ニプロ株式会社 | 判定装置、判定方法及びコンピュータプログラム |
WO2016152159A1 (ja) * | 2015-03-26 | 2016-09-29 | 東洋製罐グループホールディングス株式会社 | Dnaチップ画像のスポット有効性判定装置、dnaチップ画像のスポット有効性判定方法、及びdnaチップ画像のスポット有効性判定プログラム |
JP2016183916A (ja) * | 2015-03-26 | 2016-10-20 | 東洋製罐グループホールディングス株式会社 | Dnaチップ画像のスポット有効性判定装置、dnaチップ画像のスポット有効性判定方法、及びdnaチップ画像のスポット有効性判定プログラム |
Also Published As
Publication number | Publication date |
---|---|
EP2824444A1 (en) | 2015-01-14 |
BR112014021383A2 (ja) | 2018-05-08 |
EP2824444A4 (en) | 2016-01-06 |
US20150098611A1 (en) | 2015-04-09 |
US9519822B2 (en) | 2016-12-13 |
KR20140132343A (ko) | 2014-11-17 |
JP2013186007A (ja) | 2013-09-19 |
KR102048599B1 (ko) | 2019-11-25 |
EP2824444B1 (en) | 2018-04-18 |
CN104160265A (zh) | 2014-11-19 |
JP5895613B2 (ja) | 2016-03-30 |
CN104160265B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6970234B2 (ja) | 細胞分泌プロファイルの分析およびスクリーニング | |
CA2825727C (en) | Microarray analysis method and microarray reading device | |
JP2017504307A (ja) | アレイ上のフィーチャーをデジタルカウントするための方法およびシステム | |
JP2008256428A (ja) | マイクロアレイ画像のブロック位置検出方法 | |
JP5895613B2 (ja) | 判定方法、判定装置、判定システム、および、プログラム | |
JP2018529947A (ja) | 生物学的試料の多重分析におけるクロストークの修正 | |
JPWO2002001477A1 (ja) | 遺伝子発現データの処理方法および処理プログラム | |
JP3537752B2 (ja) | バイオチップを用いたハイブリダイゼーション反応の実験結果表示方法及び実験誤差評価方法 | |
JP2019505884A (ja) | デジタル画像内の少なくとも1つのオブジェクトの全体的明るさを判定するための方法 | |
US7068828B2 (en) | Biochip image analysis system and method thereof | |
JP3880361B2 (ja) | 蛍光シグナル処理方法及びハイブリダイゼーション反応結果表示方法 | |
US20050203708A1 (en) | Method and system for microarray gradient detection and characterization | |
US20190012770A1 (en) | Method for Determining the Positions of a Plurality of Objects in a Digital Image | |
JP2013224894A (ja) | マイクロアレイの解析方法 | |
JP4271688B2 (ja) | バイオチップ | |
US20050226535A1 (en) | Method and system for rectilinearizing an image of a microarray having a non-rectilinear feature arrangement | |
US20050026306A1 (en) | Method and system for generating virtual-microarrays | |
JP2004113188A (ja) | 発現レベル情報の解析方法および解析装置 | |
KR20050048729A (ko) | 바이오칩 촬영영상 품질 분석 시스템 및 그 방법 |
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: 13758079 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 20147023805 Country of ref document: KR Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2013758079 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14383601 Country of ref document: US |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01A Ref document number: 112014021383 Country of ref document: BR |
|
ENP | Entry into the national phase |
Ref document number: 112014021383 Country of ref document: BR Kind code of ref document: A2 Effective date: 20140828 |