WO2023182169A1 - 細胞評価方法および細胞評価装置 - Google Patents

細胞評価方法および細胞評価装置 Download PDF

Info

Publication number
WO2023182169A1
WO2023182169A1 PCT/JP2023/010411 JP2023010411W WO2023182169A1 WO 2023182169 A1 WO2023182169 A1 WO 2023182169A1 JP 2023010411 W JP2023010411 W JP 2023010411W WO 2023182169 A1 WO2023182169 A1 WO 2023182169A1
Authority
WO
WIPO (PCT)
Prior art keywords
wavelength
refractive index
cell
labeling
image
Prior art date
Application number
PCT/JP2023/010411
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 浜松ホトニクス株式会社
Publication of WO2023182169A1 publication Critical patent/WO2023182169A1/ja

Links

Images

Classifications

    • 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/02—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
    • C12Q1/04—Determining presence or kind of microorganism; Use of selective media for testing antibiotics or bacteriocides; Compositions containing a chemical indicator therefor
    • C12Q1/06—Quantitative determination
    • 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/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
    • 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/41—Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
    • 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/483—Physical analysis of biological material

Definitions

  • the cell can be evaluated (for example, quality is evaluated) based on the phase image or refractive index distribution image of the cell obtained by the above-mentioned conventional technique.
  • the type of cell is unknown, it may not be clear what information about the cell the phase image and refractive index distribution image of the cell obtained by the above-mentioned conventional techniques represent, and these images It is difficult to evaluate cells based on
  • FIG. 28 shows images of Example 1B, showing (a) a phase image at a center wavelength of 433 nm, and (b) a phase image at a center wavelength of 694 nm.
  • FIG. 29 shows each image of Example 1B, showing (a) an intensity image at a center wavelength of 530 nm, and (b) a refractive index variation image.
  • FIG. 30 shows images of Example 2A, showing (a) a phase image at a center wavelength of 530 nm, and (b) a phase image at a center wavelength of 694 nm.
  • FIG. 31 shows images of Example 2A, showing (a) an intensity image at a center wavelength of 530 nm, and (b) a refractive index variation image.
  • FIG. 29 shows each image of Example 1B, showing (a) an intensity image at a center wavelength of 530 nm, and (b) a refractive index variation image.
  • FIG. 30 shows images of Example 2A, showing (a) a
  • a lens 12, a polarizer 13, a first prism 14, and a condenser lens 15 are provided in this order on the optical path of an irradiation optical system that condenses and irradiates the evaluation object S with the light output from the light source 11.
  • An objective lens 21, a second prism 22, a lens 24, a lens 31, a quarter wavelength plate 33, and a lens 34 are provided in this order.
  • the condenser lens 15 is optically connected to the first prism 14.
  • the condenser lens 15 condenses and irradiates the evaluation object S with two mutually orthogonal linearly polarized lights output from the first prism 14 . Since the two linearly polarized lights that are orthogonal to each other output from the first prism 14 travel in different directions from each other after being output from the first prism 14, the condensing position on the evaluation target S by the condenser lens 15 is small. different.
  • the calculation unit 39 and the analysis unit 40 include a processing unit including a CPU that performs various calculation processes, a storage unit including a hard disk drive, RAM, ROM, etc. for storing data and programs, and a liquid crystal display for displaying processing results, etc. and an input section including a keyboard, a mouse, etc. that accepts input of various conditions for acquiring the interference image and displaying the image.
  • the calculation unit 39 and the analysis unit 40 may be configured by a smart device such as a tablet terminal that includes a touch panel or the like as an input unit. Further, the calculation units and storage units of the calculation unit 39 and the analysis unit 40 may be configured by an FPGA (field-programmable gate array) or a microcomputer.
  • the beam splitter 57 is optically connected to the dichroic mirror 56.
  • the beam splitter 57 inputs the light that has arrived from the dichroic mirror 56 and splits it into two, outputs one branched light to the lens 61 and outputs the other branched light to the lens 62 .
  • the lens 71 is optically connected to the light output end 63b of the optical fiber 63, and collimates the light output as divergent light from the light output end 63b.
  • Mirror 72 is optically connected to lens 71 and reflects light arriving from lens 71 to lens 73. The orientation of the reflective surface of mirror 72 is variable.
  • a fluorescent pigment protein may be used as the pigment protein.
  • fluorescent pigment proteins include red fluorescent protein (eforRed) derived from coral (Echinopora forskaliana).
  • the distribution at other specific parts of the cell is also evaluated.
  • the first wavelength or the second wavelength may be equal to the third or fourth wavelength.
  • the number of wavelengths of light used in the refractive index distribution acquisition step S2 may be four, three, or two.
  • FIG. 7 is a diagram schematically illustrating the refractive index distribution of a cell whose nucleus as a specific region is labeled with a labeling substance.
  • FIG. 7(a) shows the refractive index distribution of the cell at the first wavelength ⁇ 1 .
  • FIG. 7(b) shows the refractive index distribution of the cell at the second wavelength ⁇ 2 .
  • FIG. 34(a) is a phase image at a center wavelength of 433 nm.
  • FIG. 34(b) is a phase image at a center wavelength of 620 nm.
  • FIG. 35(a) is an intensity image at a center wavelength of 530 nm.
  • FIG. 35(b) is a refractive index variation image.
  • the refractive index variation image (FIG. 35(b)) was created in the same manner as in Example 1A. Comparing the intensity image (FIG. 35(a)) and the refractive index variation image (FIG. 35(b)) shows that in Example 3, the region labeled with the labeling substance was similarly detected as a change in refractive index. I know that there is.
  • Example 38 and 39 are images of Example 4.
  • Janus Green B was used as the labeling substance
  • light with center wavelengths of 530 nm and 694 nm was used when measuring the refractive index distribution
  • A549 cell line was used as the cell to be evaluated.
  • Comparative Example 4 differs from Example 4 in that no labeling substance was used.
  • the distribution of specific regions labeled with the labeling substance can also be seen in the intensity images (such as FIG. 27(a)).
  • the intensity image is affected not only by light absorption by the labeling substance but also by light scattering and the like.
  • Intensity images may show reduced light intensity due not only to absorption effects but also to refraction and scattering of light at the boundaries of structures with different refractive indices (this is also observed in unlabeled cells). ).

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Organic Chemistry (AREA)
  • Pathology (AREA)
  • General Physics & Mathematics (AREA)
  • Zoology (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Biomedical Technology (AREA)
  • Wood Science & Technology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Toxicology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • Microbiology (AREA)
  • Hematology (AREA)
  • Urology & Nephrology (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • General Engineering & Computer Science (AREA)
  • Genetics & Genomics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

細胞評価方法は、標識ステップＳ１、屈折率分布取得ステップＳ２および解析ステップＳ３を備える。標識ステップＳ１では、第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で、評価対象である細胞の特定部位を標識する。屈折率分布取得ステップＳ２では、標識ステップＳ１で特定部位が標識された細胞の屈折率分布を、第１波長および第２波長それぞれにおいて取得する。解析ステップＳ３では、第１波長および第２波長それぞれにおける屈折率分布を比較することで、細胞における特定部位の分布を評価する。これにより、未知の細胞であっても該細胞を容易に評価することができる方法が実現される。

Description

細胞評価方法および細胞評価装置
　本開示は、細胞における特定部位の分布を評価する方法および装置に関するものである。
　例えば細胞等の透明な観察対象物は通常の顕微鏡により観察することが困難であることから、一般に観察対象物の位相画像または屈折率分布画像を取得することができる観察装置が用いられる。このような観察装置には幾つかの実現方法が知られている。
　従来技術の一つは、２光束の分岐および合波を行うマイケルソン干渉計またはマッハツェンダ干渉計を用い、２光束のうちの一方の光路上に観察対象物を配置する。そして、光路上のミラー等を移動させることで２光束の間の光路長差（位相差）を各値に順次に設定して、観察対象物の複数の干渉画像を取得し、これら複数の干渉画像に基づいて観察対象物の位相画像を作成することができる。
　非特許文献１に記載された技術は、微分干渉顕微鏡に１／４波長板および偏光カメラを組み合わせた構成を用いる。そして、偏光カメラにより取得された複数の偏光成分それぞれについての干渉画像に基づいて、観察対象物の位相画像を作成することができる。
　観察対象物の厚さをｄとし、光の波長をλとし、屈折率をｎとすると、位相φは、φ＝２π・ｎ・ｄ／λ なる関係式で表される。したがって、厚さｄおよび波長λが既知であれば、位相画像から屈折率分布画像を作成することができる。
　非特許文献２に記載された光回折トモグラフィ（Optical Diffraction Tomography: ＯＤＴ）技術は、２光束の分岐および合波を行うマッハツェンダ干渉計を用い、２光束のうちの一方の光路上に観察対象物を配置する。そして、観察対象物への光入射方向を様々に設定して、観察対象物の複数の干渉画像を取得し、これら複数の干渉画像に基づいて観察対象物の屈折率分布画像を作成することができる。
　特許文献１には、位相画像に基づいて幹細胞を評価する発明が開示されている。この文献に開示された発明は、幹細胞の位相画像において細胞核の部分の位相と細胞質の部分の位相とを比較して、その比較結果に基づいて幹細胞の品質の良否を評価する。
特許第５７４５９１９号公報
O. Yasuhiko, K. Takeuchi, H. Yamada, and Y. Ueda, "Single-shot quantitative phase imaging as an extension of differential interference contrast microscopy", Genes to Cells Vol.26, pp.596-610, 2021 O. Yasuhiko, K. Takeuchi, H. Yamada, and Y. Ueda, "Multiple-scattering suppressive refractive index tomography for the label-free quantitative assessment of multicellular spheroids", Biomed. Opt. Express Vol.13, pp.962-979, 2022 D. Fu et al., "Quantitative dispersion microscopy", Biomed. Opt. Express Vol.1, pp.347-353, 2010 T. Sai et al., "Designing refractive index fluids using the Kramers-Kronig relations", Faraday Discussions Vol.223, pp.136-144, 2020
　細胞の種類が既知である場合には、上記の従来技術により取得される細胞の位相画像または屈折率分布画像に基づいて、細胞を評価（例えば品質の良否を評価）することができる。しかし、細胞の種類が未知である場合には、上記の従来技術により取得される細胞の位相画像および屈折率分布画像は、細胞の如何なる情報を表しているのか明確でない場合があり、これらの画像に基づいて細胞を評価することは困難である。
　実施形態は、未知の細胞であっても該細胞を容易に評価することができる方法および装置を提供することを目的とする。
　実施形態は、細胞評価方法である。細胞評価方法は、第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で細胞の特定部位を標識する標識ステップと、標識ステップで特定部位が標識された細胞の屈折率分布を第１波長および第２波長それぞれにおいて取得する屈折率分布取得ステップと、第１波長および第２波長それぞれにおける屈折率分布を比較することで細胞における特定部位の分布を評価する解析ステップと、を備える。
　実施形態は、細胞評価装置である。細胞評価装置は、第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で特定部位が標識された細胞の屈折率分布を第１波長および第２波長それぞれにおいて取得する屈折率分布取得部と、第１波長および第２波長それぞれにおける屈折率分布を比較することで細胞における特定部位の分布を評価する解析部と、を備える。
　実施形態は、プログラムである。プログラムは、第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で特定部位が標識された細胞の屈折率分布を第１波長および第２波長それぞれにおいて取得する屈折率分布取得ステップと、第１波長および第２波長それぞれにおける屈折率分布を比較することで細胞における特定部位の分布を評価する解析ステップと、をコンピュータに実行させるためのプログラムである。
　実施形態の細胞評価方法および細胞評価装置によれば、未知の細胞であっても該細胞を容易に評価することができる。
図１は、細胞評価方法のフローチャートである。 図２は、細胞評価装置１の構成を示す図である。 図３は、細胞評価装置２Ａの構成を示す図である。 図４は、細胞評価装置２Ｂの構成を示す図である。 図５は、光吸収性標識物質の吸光スペクトルおよび屈折率スペクトルを模式的に示す図であり、（ａ）吸光スペクトル、及び（ｂ）屈折率スペクトルを示している。 図６は、（ａ），（ｂ）光吸収性標識物質の屈折率スペクトルを模式的に示す図である。 図７は、特定部位としての細胞核が光吸収性標識物質で標識された細胞の屈折率分布を模式的に説明する図であり、（ａ）第１波長λ１での細胞の屈折率分布、及び（ｂ）第２波長λ２での細胞の屈折率分布を示している。 図８は、解析ステップＳ３における第１解析例を説明する図であり、（ａ）屈折率分布画像を模式的に示す図、（ｂ）屈折率変動画像を模式的に示す図、及び（ｃ）屈折率分布画像（図８（ａ））と屈折率変動画像（図８（ｂ））とを重ねた図を示している。 図９は、解析ステップＳ３における第２解析例を説明する図であり、（ａ）細胞核を標識物質で標識して取得された屈折率変動画像を模式的に示す図、（ｂ）細胞全体を標識物質で標識して取得された屈折率変動画像を模式的に示す図、及び（ｃ）屈折率変動画像（図９（ａ））と屈折率変動画像（図９（ｂ））とを重ねた図を示している。 図１０は、解析ステップＳ３における第３解析例を説明する図であり、一定時間間隔で取得された複数の屈折率変動画像を模式的に示している。 図１１は、解析ステップＳ３における第３解析例を説明する図であり、一定時間間隔で取得された複数の屈折率変動画像（図１０）を重ねた図を示している。 図１２は、解析ステップＳ３における第３解析例を説明する図であり、一定時間間隔で取得された複数の屈折率変動画像（図１０）と屈折率分布画像とを重ねた図を示している。 図１３は、解析ステップＳ３における第４解析例を説明する図であり、（ａ）細胞全体の蛍光画像または強度画像を模式的に示す図、（ｂ）細胞核を標識物質で標識して取得された屈折率変動画像を模式的に示す図、及び（ｃ）細胞全体の蛍光画像または強度画像（図１３（ａ））と屈折率変動画像（図１３（ｂ））とを重ねた図を示している。 図１４は、解析ステップＳ３における第５解析例を説明する図であり、（ａ）細胞核の蛍光画像または強度画像を模式的に示す図、（ｂ）標識対象が未知である標識物質を用いて標識して取得された屈折率変動画像を模式的に示す図、及び（ｃ）細胞核の蛍光画像または強度画像（図１４（ａ））と屈折率変動画像（図１４（ｂ））とを重ねた図を示している。 図１５は、解析ステップＳ３における第６解析例を説明する図であり、（ａ）特定の小胞を標識物質αで標識して取得された屈折率変動画像を模式的に示す図、（ｂ）他の特定の小胞を標識物質βで標識して取得された屈折率変動画像を模式的に示す図、（ｃ）屈折率分布画像を模式的に示す図、及び（ｄ）屈折率変動画像（図１５（ａ））と屈折率変動画像（図１５（ｂ））と屈折率分布画像（図１５（ｃ））とを重ねた図を示している。 図１６は、細胞塊を模式的に示す図である。 図１７は、解析ステップＳ３における第７解析例を説明する図であり、（ａ）細胞塊を構成する各細胞の特定部位を標識物質αで標識して取得された屈折率変動画像を模式的に示す図、（ｂ）屈折率分布画像を模式的に示す図、及び（ｃ）屈折率変動画像（図１７（ａ））と屈折率分布画像（図１７（ｂ））とを重ねた図を示している。 図１８は、解析ステップＳ３における第８解析例を説明する図であり、（ａ）細胞塊を構成する各細胞の特定部位を標識物質αで標識して取得された屈折率変動画像を模式的に示す図、（ｂ）細胞塊を構成する各細胞の他の特定部位を他の標識物質βで標識して取得された屈折率変動画像を模式的に示す図、（ｃ）屈折率分布画像を模式的に示す図、及び（ｄ）屈折率変動画像（図１８（ａ））と屈折率変動画像（図１８（ｂ））と屈折率分布画像（図１８（ｃ））とを重ねた図を示している。 図１９は、各実施例において用いた光吸収性標識物質、中心波長および細胞を纏めた表である。 図２０は、０.００２％ブリリアントクレシルブルー溶液の測定された吸光スペクトルおよび推定された屈折率スペクトルを示すグラフである。 図２１は、０.００３％ニューメチレンブルー溶液の測定された吸光スペクトルおよび推定された屈折率スペクトルを示すグラフである。 図２２は、０.０６６％中性赤溶液の測定された吸光スペクトルおよび推定された屈折率スペクトルを示すグラフである。 図２３は、０.００４％ヤヌスグリーンＢ溶液の測定された吸光スペクトルおよび推定された屈折率スペクトルを示すグラフである。 図２４は、実施例１Ａの各画像であり、（ａ）中心波長４３３ｎｍでの位相画像、及び（ｂ）中心波長６９４ｎｍでの位相画像を示している。 図２５は、実施例１Ａの各画像であり、（ａ）中心波長４３３ｎｍでの屈折率分布画像、及び（ｂ）中心波長６９４ｎｍでの屈折率分布画像を示している。 図２６は、実施例１Ａの各画像であり、（ａ）図２５（ａ）中の破線矩形で示した範囲を拡大して示す図、及び（ｂ）図２５（ｂ）中の破線矩形で示した範囲を拡大して示す図を示している。 図２７は、実施例１Ａの各画像であり、（ａ）中心波長５３０ｎｍでの強度画像、及び（ｂ）屈折率変動画像を示している。 図２８は、実施例１Ｂの各画像であり、（ａ）中心波長４３３ｎｍでの位相画像、及び（ｂ）中心波長６９４ｎｍでの位相画像を示している。 図２９は、実施例１Ｂの各画像であり、（ａ）中心波長５３０ｎｍでの強度画像、及び（ｂ）屈折率変動画像を示している。 図３０は、実施例２Ａの各画像であり、（ａ）中心波長５３０ｎｍでの位相画像、及び（ｂ）中心波長６９４ｎｍでの位相画像を示している。 図３１は、実施例２Ａの各画像であり、（ａ）中心波長５３０ｎｍでの強度画像、及び（ｂ）屈折率変動画像を示している。 図３２は、実施例２Ｂの各画像であり、（ａ）中心波長５３０ｎｍでの位相画像、及び（ｂ）中心波長６９４ｎｍでの位相画像を示している。 図３３は、実施例２Ｂの各画像であり、（ａ）中心波長５３０ｎｍでの強度画像、及び（ｂ）屈折率変動画像を示している。 図３４は、実施例３の各画像であり、（ａ）中心波長４３３ｎｍでの位相画像、及び（ｂ）中心波長６２０ｎｍでの位相画像を示している。 図３５は、実施例３の各画像であり、（ａ）中心波長５３０ｎｍでの強度画像、及び（ｂ）屈折率変動画像を示している。 図３６は、比較例３の各画像であり、（ａ）中心波長４３３ｎｍでの位相画像、及び（ｂ）中心波長６２０ｎｍでの位相画像を示している。 図３７は、比較例３の各画像であり、（ａ）中心波長５３０ｎｍでの強度画像、及び（ｂ）屈折率変動画像を示している。 図３８は、実施例４の各画像であり、（ａ）中心波長５３０ｎｍでの位相画像、及び（ｂ）中心波長６９４ｎｍでの位相画像を示している。 図３９は、実施例４の各画像であり、（ａ）中心波長５３０ｎｍでの強度画像、及び（ｂ）屈折率変動画像を示している。 図４０は、比較例４の各画像であり、（ａ）中心波長５３０ｎｍでの位相画像、及び（ｂ）中心波長６９４ｎｍでの位相画像を示している。 図４１は、比較例４の各画像であり、（ａ）中心波長５３０ｎｍでの強度画像、及び（ｂ）屈折率変動画像を示している。
　以下、添付図面を参照して、細胞評価方法および細胞評価装置の実施の形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。本発明は、これらの例示に限定されるものではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
　図１は、細胞評価方法のフローチャートである。細胞評価方法は、標識ステップＳ１、屈折率分布取得ステップＳ２および解析ステップＳ３を備える。標識ステップＳ１では、第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で、評価対象である細胞の特定部位を標識する。
　屈折率分布取得ステップＳ２では、標識ステップＳ１で特定部位が標識された細胞の屈折率分布を、第１波長および第２波長それぞれにおいて取得する。解析ステップＳ３では、第１波長および第２波長それぞれにおける屈折率分布を比較することで、細胞における特定部位の分布を評価する。
　この細胞評価方法における屈折率分布取得ステップＳ２および解析ステップＳ３の各処理を実施することができる細胞評価装置は、例えば非特許文献１または非特許文献２に記載された技術を採用することで構成され得る。
　図２は、細胞評価装置１の構成を示す図である。細胞評価装置１は、屈折率分布取得部１０および解析部４０を備える。屈折率分布取得部１０は、屈折率分布取得ステップＳ２の処理を行うものであり、非特許文献１に記載された技術に基づくものである。解析部４０は、解析ステップＳ３の処理を行うものである。
　屈折率分布取得部１０は、光源１１、レンズ１２、偏光子１３、第１プリズム１４、コンデンサレンズ１５、対物レンズ２１、第２プリズム２２、レンズ２４、ミラー２７、光学モジュール３０および演算部３９を備える。光学モジュール３０は、レンズ３１、１／４波長板３３、レンズ３４、偏光カメラ３５およびミラー３６～３８を含む。
　屈折率分布取得部１０のうち光学モジュール３０および演算部３９を除く構成は、微分干渉顕微鏡の構成と同様の構成とすることができる。すなわち、微分干渉顕微鏡のカメラポートに光学モジュール３０を装着するとともに、演算部３９を設けることにより、屈折率分布取得部１０を構成することができる。
　光源１１から出力された光を評価対象物Ｓに集光照射する照射光学系の光路上に、順に、レンズ１２、偏光子１３、第１プリズム１４およびコンデンサレンズ１５が設けられている。また、照射光学系による評価対象物Ｓへの光の照射に応じて評価対象物Ｓで生じた光を入力して偏光カメラ３５の撮像面上に結像する結像光学系の光路上に、順に、対物レンズ２１、第２プリズム２２、レンズ２４、レンズ３１、１／４波長板３３およびレンズ３４が設けられている。
　光源１１は、空間的にインコヒーレントな光を出力する。光源１１から出力される光は、時間的にコヒーレントであってもよいし、時間的にインコヒーレントであってもよい。光源１１から出力される光は、直線偏光であってもよいし、無偏光であってもよい。光源１１は、例えばハロゲンランプ、発光ダイオードまたはレーザダイオードを含み、さらに、その後段に拡散板を含むのも好適である。
　光源１１は、互いに異なる中心波長を有する光を出力することができる。例えば、出力波長が互いに異なる複数のレーザダイオードのうちから使用するレーザダイオードを選択する構成とすることで、出力光の波長を可変としてもよい。また、広帯域の光を出力する発光ダイオードと、透過帯域が互いに異なる複数のフィルタとを用いて、使用するフィルタを切り替えることで、出力光の波長を可変としてもよい。
　レンズ１２は、光源１１と光学的に接続されている。レンズ１２は、光源１１から出力された光をコリメートして偏光子１３へ出力する。
　偏光子１３は、レンズ１２と光学的に接続されている。偏光子１３は、レンズ１２によりコリメートされて出力された光を入力し、その偏光子１３の光学軸の方位に応じた偏光面（振動面）を有する直線偏光の光を第１プリズム１４へ出力する。
　第１プリズム１４は、偏光子１３と光学的に接続されている。第１プリズム１４は、偏光子１３から出力された直線偏光の光を入力し、互いに直交する二つの直線偏光の光を互いに異なる方向へ出力する。第１プリズム１４から出力される互いに直交する二つの直線偏光の光それぞれの偏光面は、偏光子１３から第１プリズム１４に入力される直線偏光の光の偏光面に対して４５°だけ傾いている。
　コンデンサレンズ１５は、第１プリズム１４と光学的に接続されている。コンデンサレンズ１５は、第１プリズム１４から出力される互いに直交する二つの直線偏光の光を評価対象物Ｓに集光照射する。第１プリズム１４から出力された互いに直交する二つの直線偏光の光は、第１プリズム１４から出力された後に互いに異なる方向へ進むので、コンデンサレンズ１５による評価対象物Ｓ上の集光位置は僅かに異なる。
　対物レンズ２１は、照射光学系による評価対象物Ｓへの光の照射に応じて評価対象物Ｓから出力された互いに直交する二つの直線偏光の光を入力し、これらの光をコリメートして第２プリズム２２へ出力する。
　第２プリズム２２は、対物レンズ２１と光学的に接続されている。第２プリズム２２は、対物レンズ２１によりコリメートされて出力された互いに直交する二つの直線偏光の光を入力し、第１プリズム１４でずらした光路におけるこれらの光を合波してレンズ２４へ出力する。
　第１プリズム１４および第２プリズム２２は、通常の微分干渉顕微鏡において用いられるウォラストンプリズム（Wollaston prism）であってもよいし、ノマルスキープリズム（Nomarski prism）であってもよい。
　第２プリズム２２と偏光カメラ３５との間の光路上にあるレンズ２４、３１，３４は、第２プリズム２２から出力された光を入力し、偏光カメラ３５の撮像面上に結像する。レンズ２４は、第２プリズム２２から出力された光を入力し、この光を収斂させる。レンズ３１は、レンズ２４から出力されミラー２７で反射された光を入力し、この光をコリメートする。レンズ３４は、レンズ３１から出力されミラー３６，３７で反射された光を入力し、この光を収斂させる。レンズ３１およびレンズ３４はリレー光学系を構成している。
　レンズ３１とレンズ３４との間の光路上に設けられている１／４波長板３３は、第２プリズム２２から出力されレンズ２４，３１を経てコリメートされた光を入力し、互いに回転方向が異なる二つの円偏光の光を出力する。
　偏光カメラ３５は、レンズ３４と光学的に接続されている。偏光カメラ３５は、結像光学系により像が形成される位置に配置された撮像面を有する。偏光カメラ３５は、１／４波長板３３により互いに異なる回転方向の円偏光とされた二つの光を入力し、３以上の偏光成分それぞれについて撮像面上の干渉画像を取得する。
　偏光カメラ３５として、例えば、ソニー株式会社により商品化されているイメージセンサ（Polarsens（登録商標））や、Teledyne DALSA社により商品化されているイメージセンサ（Area Scan Polarization Sensor）やイメージセンサ（Line Scan Polarization Sensor）が用いられ得る。
　演算部３９は、偏光カメラ３５と電気的に接続されている。演算部３９は、偏光カメラ３５により取得された３以上の偏光成分それぞれについての干渉画像に基づいて評価対象物Ｓの複素振幅画像（振幅画像および位相画像）を作成する。解析部４０は、演算部３９により取得された位相画像に基づいて、評価対象物Ｓを評価する。
　演算部３９および解析部４０は、プログラムに従って屈折率分布取得ステップＳ２および解析ステップＳ３の各計算処理を実行するものであり、例えばコンピュータである。両者は一体のものであってもよい。
　演算部３９および解析部４０は、様々な演算処理を行うＣＰＵを含む処理部と、データやプログラムを記憶するハードディスクドライブ、ＲＡＭおよびＲＯＭ等を含む記憶部と、処理結果等を表示する液晶ディスプレイを含む表示部と、干渉画像の取得および画像の表示の際の諸条件の入力を受け付けるキーボードおよびマウス等を含む入力部とを備える。演算部３９および解析部４０は、タッチパネル等を入力部として備えるタブレット端末等のスマートデバイスによって構成されてもよい。また、演算部３９および解析部４０の演算部や記憶部は、ＦＰＧＡ（field-programmable gate array）やマイコンによって構成されていてもよい。
　図３は、細胞評価装置２Ａの構成を示す図である。細胞評価装置２Ａは、屈折率分布取得部５０Ａおよび解析部９０を備える。屈折率分布取得部５０Ａは、屈折率分布取得ステップＳ２の処理を行うものであり、非特許文献２に記載された技術に基づくものである。解析部９０は、解析ステップＳ３の処理を行うものである。
　屈折率分布取得部５０Ａは、レーザ光源５１，５２、ミラー５５、ダイクロイックミラー５６、ビームスプリッタ５７、レンズ６１，６２、光ファイバ６３，６４、レンズ７１、ミラー７２、レンズ７３、コンデンサレンズ７４、対物レンズ７５、ビームスプリッタ８１、レンズ８２、カメラ８３Ａおよび演算部８４を備える。
　レーザ光源５１，５２は、互いに異なる中心波長のレーザ光を出力する。ミラー５５は、レーザ光源５１と光学的に接続されている。ミラー５５は、レーザ光源５１から出力されたレーザ光をダイクロイックミラー５６へ反射させる。ダイクロイックミラー５６は、レーザ光源５２およびミラー５５と光学的に接続されている。ダイクロイックミラー５６は、レーザ光源５２から出力されたレーザ光を入力するとともに、ミラー５５から到達したレーザ光を入力して、これらの入力したレーザ光を合波してビームスプリッタ５７へ出力する。
　ビームスプリッタ５７は、ダイクロイックミラー５６と光学的に接続されている。ビームスプリッタ５７は、ダイクロイックミラー５６から到達した光を入力して２分岐し、一方の分岐光をレンズ６１へ出力し、他方の分岐光をレンズ６２へ出力する。
　レンズ６１は、ビームスプリッタ５７と光学的に接続されている。レンズ６１は、ビームスプリッタ５７から到達した一方の分岐光を入力して、その光を光ファイバ６３の光入射端６３ａに集光して、その光を光入射端６３ａに入射させる。光ファイバ６３は、光入射端６３ａに入射した光を導光し、光出射端６３ｂから発散光として出射する。
　レンズ６２は、ビームスプリッタ５７と光学的に接続されている。レンズ６２は、ビームスプリッタ５７から到達した他方の分岐光を入力して、その光を光ファイバ６４の光入射端６４ａに集光して、その光を光入射端６４ａに入射させる。光ファイバ６４は、光入射端６４ａに入射した光を導光し、光出射端６４ｂから発散光として出射する。
　レンズ７１は、光ファイバ６３の光出射端６３ｂと光学的に接続されており、光出射端６３ｂから発散光として出力された光をコリメートする。ミラー７２は、レンズ７１と光学的に接続されており、レンズ７１から到達した光をレンズ７３へ反射させる。ミラー７２の反射面の方位は可変である。
　レンズ７３は、ミラー７２と光学的に接続されている。コンデンサレンズ７４は、レンズ７３と光学的に接続されている。レンズ７３およびコンデンサレンズ７４は、好適には４ｆ光学系を構成している。レンズ７３およびコンデンサレンズ７４は、ミラー７２の反射面の方位に応じた光照射方向から評価対象物Ｓに対して光を照射する。
　対物レンズ７５は、コンデンサレンズ７４と光学的に接続されている。対物レンズ７５とコンデンサレンズ７４との間に評価対象物Ｓが配置される。対物レンズ７５は、コンデンサレンズ７４から出力されて評価対象物Ｓを経た光（物体光）を入力し、その光をビームスプリッタ８１へ出力する。
　ビームスプリッタ８１は、対物レンズ７５と光学的に接続され、また、光ファイバ６４の光出射端６４ｂとも光学的に接続されている。ビームスプリッタ８１は、対物レンズ７５から出力されて到達した光（物体光）を入力するとともに、光出射端６４ｂから出力されて到達した光（参照光）を入力して、両光をレンズ８２へ出力する。レンズ８２は、ビームスプリッタ８１と光学的に接続されており、ビームスプリッタ８１から到達した物体光および参照光それぞれをコリメートしてカメラ８３Ａへ出力する。
　カメラ８３Ａは、カラーカメラであり、レンズ８２と光学的に接続されている。カメラ８３Ａは、レンズ８２から到達した物体光と参照光との干渉による干渉縞像（干渉強度画像）を撮像する。カメラ８３Ａの撮像面への物体光の入射方向に対して参照光の入射方向は傾斜している。ビームスプリッタ８１により物体光と参照光とが合波される位置は、結像レンズより後段であってもよいが、収差の影響を考慮すると、図に示されるように対物レンズ７５とレンズ８２との間であるのが望ましい。
　演算部８４は、カメラ８３Ａと電気的に接続されている。演算部８４は、カメラ８３Ａにより取得された干渉画像（カラー画像）を解析して、レーザ光源５１，５２それぞれの出力波長での干渉画像を作成する。さらに、演算部８４は、ミラー７２の反射面の方位の変化により評価対象物Ｓへの光入射方向を様々に設定して得られた干渉画像を処理することにより、各波長における評価対象物Ｓの３次元屈折率分布を算出する。
　解析部９０は、演算部８４により取得された屈折率分布画像に基づいて、評価対象物Ｓを評価する。演算部８４および解析部９０も、屈折率分布取得ステップＳ２および解析ステップＳ３の各計算処理を実行するものであり、例えばコンピュータである。両者は一体のものであってもよい。
　図４は、細胞評価装置２Ｂの構成を示す図である。細胞評価装置２Ｂは、屈折率分布取得部５０Ｂおよび解析部９０を備える。屈折率分布取得部５０Ａ（図３）と比較すると、屈折率分布取得部５０Ｂ（図４）は、シャッタ５３，５４を更に備える点で相違し、カラーカメラ８３Ａに替えてモノクロカメラ８３Ｂを備える点で相違する。
　シャッタ５３は、レーザ光源５１とミラー５５との間の光路上に設けられており、レーザ光源５１からミラー５５へのレーザ光の伝搬を許可または禁止する。シャッタ５４は、レーザ光源５２とダイクロイックミラー５６との間の光路上に設けられており、レーザ光源５２からダイクロイックミラー５６へのレーザ光の伝搬を許可または禁止する。シャッタ５３およびシャッタ５４は、同時にレーザ光の伝搬を許可することはない。カメラ８３Ｂは、レーザ光源５１の出力波長での干渉画像と、レーザ光源５２の出力波長での干渉画像と、を個別に撮像する。
　本実施形態の細胞評価方法の屈折率分布取得ステップＳ２の処理は、これらの細胞評価装置の屈折率分布取得部１０，５０Ａ，５０Ｂを用いて実施することができるが、これらの構成に限られるものではなく、評価対象物Ｓの位相画像または屈折率分布画像を取得することができる他の構成を用いて実施してもよい。
　次に、細胞評価方法の標識ステップＳ１、屈折率分布取得ステップＳ２および解析ステップＳ３それぞれについて、より詳細に説明する。
　一般に、光の波長が長いほど物質の屈折率は低くなることが知られている。例えば、非特許文献３には、タンパク質およびＤＮＡそれぞれについて、波長３１０ｎｍでの屈折率と比べて、波長４００ｎｍでの屈折率が低いことが示されている。また、非標識部位では、光の波長が長いほど、非標識部位全体の屈折率が低くなってしまう。そのため、標識（非標識部位の屈折率の変化よりも大きな屈折率の変化を示す標識）を行わないと、細胞内における特定部位を観察することは容易ではない。
　標識ステップＳ１において、評価対象である細胞の特定部位を標識する標識物質として、後の屈折率分布取得ステップＳ２で用いる光の波長である第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質を用いる。特定部位は、特定分子または特定構造であってよく、例えば、小胞構造、核小体、ミトコンドリア、特定のタンパク質などであり、細胞全体であってもよい。
　標識ステップＳ１において特定部位を標識する標識物質は、第１波長および第２波長のうちの長波長における屈折率が短波長における屈折率より高い標識物質であってもよく、逆に、長波長における屈折率が短波長における屈折率より低い標識物質であってもよい。
　また、標識ステップＳ１において特定部位を標識する標識物質は、特定波長λ０において吸収ピークを有する光吸収性標識物質であるのが好適である。この場合、屈折率分布取得ステップＳ２において、特定波長λ０に対して短波長側または長波長側にある第１波長λ１および第２波長λ２それぞれにおいて細胞の屈折率分布を取得するのが好適である。また、屈折率分布取得ステップＳ２において、第１波長λ１および第２波長λ２のうちいずれか一方を特定波長λ０に対して短波長側とし、他方を特定波長λ０に対して長波長側として、第１波長λ１および第２波長λ２のそれぞれにおいて細胞の屈折率分布を取得するのも好適である。
　図５は、光吸収性標識物質の吸光スペクトルおよび屈折率スペクトルを模式的に示す図である。図５（ａ）は吸光スペクトルを示す。図５（ｂ）は屈折率スペクトルを示す。この図に示されるように、光吸収性標識物質は、特定波長λ０において吸収ピークを有し、この特定波長λ０に対して短波長側において屈折率が極小となり、特定波長λ０に対して長波長側において屈折率が極大となる。一般に、吸光スペクトルと屈折率スペクトルとの間の関係は、クラマース・クローニッヒの関係式で表される（非特許文献４参照）。
　そこで、特定波長λ０において吸収ピークを有する光吸収性標識物質を用いる場合には、図５（ｂ）に示されるように、特定波長λ０に対して短波長側において屈折率が極小となる波長付近に第１波長λ１を設定し、特定波長λ０に対して長波長側において屈折率が極大となる波長付近に第２波長λ２を設定することができる。
　図６（ａ）の光吸収性標識物質の屈折率スペクトルに示されるように、屈折率が極小となる波長付近から短波長側において屈折率の波長依存性が大きい波長範囲があるので、この波長範囲に第１波長λ１および第２波長λ２を設定してもよい。また、図６（ｂ）の光吸収性標識物質の屈折率スペクトルに示されるように、屈折率が極大となる波長付近から長波長側においても屈折率の波長依存性が大きい波長範囲があるので、この波長範囲に第１波長λ１および第２波長λ２を設定してもよい。
　このように第１波長λ１および第２波長λ２を設定して、屈折率分布取得ステップＳ２では、光吸収性標識物質で特定部位が標識された細胞の屈折率分布を第１波長λ１および第２波長λ２それぞれにおいて取得する。
　生きた細胞の特定部位を標識することができる光吸収性標識物質の例として、ブリリアントクレシルブルー、ニューメチレンブルー、中性赤、ヤヌスグリーンＢ、ナイルブルー、ビスマルクブラウン、ナイルレッド、オイルレッド、等が挙げられる。死細胞の特定部位を標識することができる光吸収性標識物質の例として、トリパンブルー等が挙げられる。
　また、細胞の特定部位を標識することができる光吸収性標識物質の例として、細胞の固定処理が必要ではあるが、ヘマトキシリン、エオシン、ピクリン酸、オレンジＧ、アゾカルミンＧ、フクシン、アニリン青、アズール青、アズールII、ズダンIII、オイズダン黒、レゾルシン、ケルンエヒトロート、等も挙げられる。
　細胞の特定部位を標識することができる光吸収性標識物質の例として、天然の色素タンパク質を用いることもできる。このような天然の色素タンパク質の例として、イソギンチャク（Cnidopus japonicus）由来の青色タンパク質（cjBlue）、クラゲ（Rhizostoma pulmo）由来の青色タンパク質（rpulFKz1）、サンゴ（Galaxea fascicularis）由来の紫色タンパク質（gfasCP）、イソギンチャク（Anemonia sulcata）由来のピンク色タンパク質（asFP595）、サンゴ（Montipora efflorescens）由来の青色タンパク質（Rtms5）、イソギンチャク（Actinia equina）由来の青色タンパク質（aeCP597）等が挙げられる。天然の色素タンパク質を改変した人工色素タンパク質を用いてもよい。
　また、色素タンパク質として、蛍光色素タンパク質を用いてもよい。このような蛍光色素タンパク質の例として、サンゴ（Echinopora forskaliana）由来の赤色蛍光タンパク質（eforRed）等が挙げられる。
　細胞の特定部位を標識する標識物質は、色素タンパク質Ｖ（光吸収性ポリペプチド）を指令するＤＮＡ配列が任意タンパク質Ｗを指令するＤＮＡ配列のＮ端末またはＣ端末に接続された細胞内発現用コンストラクトであってもよい。ＣＭＶ（Cytomegalovirus）プロモータ等の制御下で発現するように細胞内発現用コンストラクトを設計してもよい。標識ステップＳ１では、リポフェクション法やエレクトロポレーション法等によって、この細胞内発現用コンストラクトを細胞に導入して、細胞の特定部位を標識する。この導入によって人工融合タンパク質が発現する。
　複数種類の標識物質を同時に用いてもよい。すなわち、標識ステップＳ１において、第３波長および第４波長それぞれにおける屈折率が互いに異なる他の標識物質で細胞の他の特定部位を標識してもよい。この場合、屈折率分布取得ステップＳ２において、標識ステップＳ１で標識された細胞の屈折率分布を、第１波長および第２波長それぞれにおいて取得するとともに、第３波長および第４波長それぞれにおいても取得する。
　解析ステップＳ３において、第３波長および第４波長それぞれにおける屈折率分布を比較することで、細胞における他の特定部位の分布をも評価する。第１波長または第２波長は、第３波長または第４波長と等しくてもよい。２種類の標識物質を同時に用いる場合、屈折率分布取得ステップＳ２において用いる光の波長数は、４であってもよいし、３であってもよいし、２であってもよい。
　第１波長および第２波長は例えば次のようにして決めることができる。光吸収性標識物質の吸光スペクトルが既知である場合には、クラマース・クローニッヒの関係式に基づいて吸光スペクトルから屈折率スペクトルを推定し、この屈折率プロファイルにおける極小値付近・極大値付近に第１波長・第２波長を設定すればよい。
　光吸収性標識物質の吸光スペクトルが未知である場合には、吸光スペクトルを計測した後に、クラマース・クローニッヒの関係式に基づいて吸光スペクトルから屈折率スペクトルを推定し、この屈折率プロファイルにおける極小値付近・極大値付近に第１波長・第２波長を設定すればよい。
　第１波長・第２波長を設定した後に、その設定した第１波長・第２波長について屈折率分布取得ステップＳ２の処理を行ってもよい。また、屈折率分布取得ステップＳ２において幾つかの波長について屈折率分布を取得した後に、それらの波長のうちから屈折率分布の差異が認められる第１波長・第２波長を設定してもよい。
　第１波長・第２波長を自由に設定することができる場合には、屈折率プロファイルにおける極小値・極大値に第１波長・第２波長を設定すればよい。幾つかの固定の波長のうちから選択して第１波長・第２波長を選択する場合には、それらの固定の波長のうちから最も屈折率差が大きい２つの波長を第１波長・第２波長としてもよい。また、第１波長・第２波長のうち一方の波長については常に屈折率分布を取得し、他方の波長については光吸収性標識物質の色に基づいて設定した波長で屈折率分布を取得してもよい。
　解析ステップＳ３では、第１波長λ１および第２波長λ２それぞれにおける屈折率分布を比較することで、細胞における特定部位の分布を評価する。図７は、特定部位としての細胞核が標識物質で標識された細胞の屈折率分布を模式的に説明する図である。図７（ａ）は第１波長λ１での細胞の屈折率分布を示す。図７（ｂ）は第２波長λ２での細胞の屈折率分布を示す。
　第１波長λ１での細胞の屈折率分布（図７（ａ））と比べて、第２波長λ２での細胞の屈折率分布（図７（ｂ））では、標識物質で標識されている細胞核の領域（図７（ｂ）中のハッチング領域）の屈折率を高くすることができる。第１波長λ１での細胞の屈折率分布（図７（ａ））と第２波長λ２での細胞の屈折率分布（図７（ｂ））とを比較することで、細胞における細胞核の位置、形状および大きさ等を評価することができる。このように、解析ステップＳ３では、第１波長λ１および第２波長λ２それぞれにおける屈折率分布を比較することで、細胞における特定部位の分布を評価することができる。
　解析ステップＳ３では、第１波長λ１および第２波長λ２それぞれでの細胞の屈折率分布の差に基づいて、細胞における特定部位の分布を評価してもよい。また、これら二つの屈折率分布の比に基づいて、細胞における特定部位の分布を評価してもよい。細胞の位相画像が取得された後に屈折率分布画像が作成される場合には、波長の違いに基づく補正をした後に二つの屈折率分布の差または比を求めてもよいし、波長の違いに基づく補正をすることなく二つの屈折率分布の比を求めてもよい。
　第１波長λ１および第２波長λ２それぞれにおける屈折率分布の比較の結果を細胞の強度画像または蛍光画像に重ねることで、細胞の特定部位に関する解析を行うのも好適である。以下では、第１波長λ１および第２波長λ２それぞれでの細胞の屈折率分布画像の差または比を表す画像を「屈折率変動画像」という。
　解析ステップＳ３では、屈折率分布画像を用いて（さらに他の画像をも用いて）細胞の特定部位に関して様々な解析を行うことができる。
　図８は、解析ステップＳ３における第１解析例を説明する図である。図８（ａ）は、屈折率分布画像を模式的に示す図である。図８（ｂ）は、屈折率変動画像を模式的に示す図である。図８（ｃ）は、屈折率分布画像（図８（ａ））と屈折率変動画像（図８（ｂ））とを重ねた図である。これらの図は、細胞内の小胞構造や細胞構造のうちでも特定の小胞構造を標識物質で標識した例を示している。
　屈折率分布画像（図８（ａ））では、細胞内において標識された特定の小胞構造と標識されていない小胞構造や他の細胞構造とを区別することができない。一方、屈折率変動画像（図８（ｂ））では、細胞内において標識された特定の小胞構造と標識されていない小胞構造や他の細胞構造とを区別することができる。
　したがって、屈折率分布画像（図８（ａ））と屈折率変動画像（図８（ｂ））とを比較することにより、或いは、屈折率分布画像（図８（ａ））と屈折率変動画像（図８（ｂ））とを重ねて表示することにより、細胞における小胞構造や他の細胞構造の分布を求めることができ、さらに、そのうちでも標識された特定の小胞構造の分布をも求めることができる。
　第１解析例では、標識された特定の小胞構造の平均屈折率等を解析することができ、また、この値に基づいて、特定の小胞構造を評価することができる。例えば、小胞内に含まれる物質を推定することができ、空胞や分泌小胞を識別することができる。
　標識物質の局在性が未知であって、屈折率分布画像から大凡の細胞内構造が特定できる場合、標識物質がどこに局在するかを評価することができ、例えば、特定の小胞構造の細胞内の局在性（細胞核周辺への局在傾向等）が分かる。標識物質の局在性が既知であって、屈折率分布画像から細胞内構造を特定することが難しい場合、屈折率分布画像の区分化（セグメンテーション）を行うことができ、例えば、どのような屈折率の特徴を有するものが特定の小胞構造であると判別しうるか評価することができる。
　なお、屈折率分布画像の取得の際に用いる光の波長は、屈折率変動画像の取得の際に用いる光の何れかの波長と等しくてもよいし、異なってもよい。例えば、標識部位と非標識部位との間で屈折率を比較したい場合等では、屈折率分布画像の取得の際に用いる光の波長として、標識物質の屈折率スペクトルが大きく変化する波長領域ではないものを使う等の工夫をしてもよい。
　図９は、解析ステップＳ３における第２解析例を説明する図である。図９（ａ）は、細胞核を標識物質で標識して取得された屈折率変動画像を模式的に示す図である。図９（ｂ）は、細胞全体を標識物質で標識して取得された屈折率変動画像を模式的に示す図である。図９（ｃ）は、屈折率変動画像（図９（ａ））と屈折率変動画像（図９（ｂ））とを重ねた図である。
　第２解析例では、屈折率変動画像（図９（ａ））に基づいて細胞核の面積を求めとともに、屈折率変動画像（図９（ｂ））に基づいて細胞全体の面積を求めて、両者の面積の比を求めることで、細胞核／細胞質の比を解析することができる。また、細胞全体で細胞核以外の領域を細胞質と特定することもでき、細胞全体における細胞核の位置を特定することもできる。細胞核の標識物質のみで細胞を標識し、細胞核の数をカウントすることで、細胞塊に含まれる細胞数を推定することもできる。
　図１０～図１２は、解析ステップＳ３における第３解析例を説明する図である。図１０は、一定時間間隔で取得された複数の屈折率変動画像を模式的に示す図である。図１１は、一定時間間隔で取得された複数の屈折率変動画像（図１０）を重ねた図である。図１２は、一定時間間隔で取得された複数の屈折率変動画像（図１０）と屈折率分布画像とを重ねた図である。
　第３解析例では、標識物質で標識された特定部位の位置が細胞内において経時変化している様子を解析することができる。蛍光画像を取得する場合には、高パワーの励起光を照射する必要があることから光毒性や光退色の問題が生じるが、本実施形態では、高パワーの光を照射する必要がないので光毒性や光退色の問題が少なく、長時間に亘る細胞のタイムラプス観察が可能である。
　図１３は、解析ステップＳ３における第４解析例を説明する図である。図１３（ａ）は、細胞全体の蛍光画像または強度画像を模式的に示す図である。図１３（ｂ）は、細胞核を標識物質で標識して取得された屈折率変動画像を模式的に示す図である。図１３（ｃ）は、細胞全体の蛍光画像または強度画像（図１３（ａ））と屈折率変動画像（図１３（ｂ））とを重ねた図である。
　第４解析例では、細胞全体の蛍光画像または強度画像（図１３（ａ））に基づいて細胞全体の面積を求めるとともに、屈折率変動画像（図１３（ｂ））に基づいて細胞核の面積を求めて、両者の面積の比を求めることで、細胞核／細胞質の比を解析することができる。また、細胞全体で細胞核以外の領域を細胞質と特定することもでき、細胞全体における細胞核の位置を特定することもできる。
　図１４は、解析ステップＳ３における第５解析例を説明する図である。図１４（ａ）は、細胞核の蛍光画像または強度画像を模式的に示す図である。図１４（ｂ）は、標識対象が未知である標識物質を用いて標識して取得された屈折率変動画像を模式的に示す図である。図１４（ｃ）は、細胞核の蛍光画像または強度画像（図１４（ａ））と屈折率変動画像（図１４（ｂ））とを重ねた図である。
　第５解析例では、細胞核の蛍光画像または強度画像（図１４（ａ））と屈折率変動画像（図１４（ｂ））とが互いに一致すれば、標識対象が未知であるとしていた標識物質が細胞核を標識し得るものであると評価することができる。また、細胞核の蛍光画像または強度画像（図１４（ａ））と屈折率変動画像（図１４（ｂ））とを比較することで、屈折率分布画像に基づいて屈折率変動画像を作成する際の補正を行うこともできる。
　なお、蛍光画像または強度画像の取得の際に用いる光の波長は、屈折率変動画像の取得の際に用いる光の何れかの波長と等しくてもよいし、異なってもよい。
　２種類の標識物質を用いる場合、解析ステップＳ３では、２つの屈折率分布画像を用いて（さらに他の画像をも用いて）細胞の２種類の特定部位に関して更に様々な解析を行うことができる。
　図１５は、解析ステップＳ３における第６解析例を説明する図である。図１５（ａ）は、特定の小胞を標識物質αで標識して取得された屈折率変動画像を模式的に示す図である。図１５（ｂ）は、他の特定の小胞を標識物質βで標識して取得された屈折率変動画像を模式的に示す図である。図１５（ｃ）は、屈折率分布画像を模式的に示す図である。図１５（ｄ）は、屈折率変動画像（図１５（ａ））と屈折率変動画像（図１５（ｂ））と屈折率分布画像（図１５（ｃ））とを重ねた図である。
　第６解析例では、標識物質αで標識される特定の小胞が存在している領域と、標識物質βで標識される他の特定の小胞が存在している領域とが互いに異なることを知ることができ、例えば、前者の小胞が細胞核に近い領域に存在し、後者の小胞が細胞核から遠い領域に存在する、といったことを知ることができる。
　単一の細胞だけでなく、複数の細胞が３次元状に集合した細胞塊についても、特定部位の分布を評価することができる。
　図１６は、細胞塊を模式的に示す図である。細胞塊の任意の断面（図中で破線で示した断面）における位相画像または屈折率分布画像を求めることができ、これから該断面における屈折率変動画像を取得することができる。細胞塊の断面だけでなく断層や３次元の屈折率分布画像を求めることもできる。
　図１７は、解析ステップＳ３における第７解析例を説明する図である。図１７（ａ）は、細胞塊を構成する各細胞の特定部位を標識物質αで標識して取得された屈折率変動画像を模式的に示す図である。図１７（ｂ）は、屈折率分布画像を模式的に示す図である。図１７（ｃ）は、屈折率変動画像（図１７（ａ））と屈折率分布画像（図１７（ｂ））とを重ねた図である。
　第７解析例では、細胞塊において、標識物質αで標識される特定部位を有する細胞が存在する領域を知ることができ、例えば、周縁領域に存在することを知ることができる。
　２種類の標識物質を用いて、細胞塊について２種類の特定部位の分布を評価することもできる。
　図１８は、解析ステップＳ３における第８解析例を説明する図である。図１８（ａ）は、細胞塊を構成する各細胞の特定部位を標識物質αで標識して取得された屈折率変動画像を模式的に示す図である。図１８（ｂ）は、細胞塊を構成する各細胞の他の特定部位を他の標識物質βで標識して取得された屈折率変動画像を模式的に示す図である。図１８（ｃ）は、屈折率分布画像を模式的に示す図である。図１８（ｄ）は、屈折率変動画像（図１８（ａ））と屈折率変動画像（図１８（ｂ））と屈折率分布画像（図１８（ｃ））とを重ねた図である。
　第８解析例では、細胞塊において、標識物質α，βそれぞれで標識される特定部位を有する細胞が存在する領域を知ることができ、例えば、標識物質αで標識される特定部位を有する細胞が周縁領域に存在すること、および、標識物質βで標識される特定部位を有する細胞が中央領域に存在すること、を知ることができる。
　次に実施例について説明する。株式会社ニコン製の微分干渉顕微鏡（ECLIPSE Ti-2）のカメラポートに光学モジュール３０を装着して細胞評価装置１（図２）を構成し、これを用いて、標識物質で特定部位を標識した細胞の位相画像を取得し、さらに屈折率分布画像および屈折率変動画像を求めた。標識物質の種類に応じて、４３３±１２ｎｍ、 ５３０±２１.５ｎｍ、６２０±２６ｎｍおよび６９４±２２ｎｍのうちの何れか２波長の光を用いて位相画像を取得した。発光ダイオードとSemrock社製バンドパスフィルタとを組み合わせることで、各波長の光を選択して出力した。
　標識物質として、ブリリアントクレシルブルー、ニューメチレンブルー、ヤヌスグリーンＢ（以上の３種類の標識物質は武藤化学株式会社から購入）および中性赤（シグマアルドリッチ社から購入）を用いた。評価対象の細胞として、ヒト肺がん由来Ａ５４９細胞株およびヒト肝がん由来ＨｅｐＧ２細胞株を用いた。
　図１９は、各実施例において用いた標識物質、中心波長および細胞を纏めた表である。この表には、標識物質を用いなかった比較例３，４についても中心波長および細胞が示されている。
　実施例１Ａ，１Ｂで用いたブリリアントクレシルブルーは、生体染色色素（塩基性色素）の一つであり、エンドサイトーシス（ピノサイトーシス；飲作用）によって細胞内に取込まれることが知られており、また、核小体の染色色素としても知られている。実施例２Ａ，２Ｂで用いたニューメチレンブルーは、ブリリアントクレシルブルーと同様の性質を有する試薬として知られている。
　実施例３で用いた中性赤は、生体染色色素の一つであり、細胞取込みアッセイに使用される。この中性赤の取込みは、ＡＴＰ依存的なｐＨ勾配の維持に関係する。実施例４で用いたヤヌスグリーンＢは、生体染色色素の一つであり、ミトコンドリアを青色に染色することが知られている他、その他の様々なタンパク質を染色することも知られている。また、ヤヌスグリーンＢは、培養細胞で細胞核や細胞質を染色することが知られている。
　標識物質の吸光スペクトルは、株式会社日立製作所製の分光光度計U-4100を用いて測定した。標識物質の屈折率スペクトルは、クラマース・クローニッヒの関係式に基づいて吸光スペクトルから推定した。
　図２０は、０.００２％ブリリアントクレシルブルー溶液の測定された吸光スペクトルおよび推定された屈折率スペクトルを示すグラフである。図２１は、０.００３％ニューメチレンブルー溶液の測定された吸光スペクトルおよび推定された屈折率スペクトルを示すグラフである。図２２は、０.０６６％中性赤溶液の測定された吸光スペクトルおよび推定された屈折率スペクトルを示すグラフである。図２３は、０.００４％ヤヌスグリーンＢ溶液の測定された吸光スペクトルおよび推定された屈折率スペクトルを示すグラフである。これらの屈折率プロファイルに基づいて、屈折率分布を測定する際に用いる光の第１波長・第２波長を図１９のとおりに設定した。
　図２４～図２７は、実施例１Ａの各画像である。実施例１Ａでは、標識物質としてブリリアントクレシルブルーを用い、屈折率分布を測定する際に中心波長４３３ｎｍ，６９４ｎｍの光を用い、評価対象の細胞としてＡ５４９細胞株を用いた。
　図２４（ａ）は、中心波長４３３ｎｍでの位相画像である。図２４（ｂ）は、中心波長６９４ｎｍでの位相画像である。図中において矢印(a)で指し示す位置は、標識物質で標識されている部位であって、中心波長４３３ｎｍで位相が１.５１であって屈折率が１.３５４であり、中心波長６９４ｎｍで位相が２.０５であって屈折率が１.３７８であり、長波長で屈折率が高くなっている。なお、細胞の厚さを５μｍとし、背景培地の屈折率を１.３３３として、屈折率の値を推定した。
　図中において矢印(b)で指し示す位置は、標識物質で標識されていない部位であって、中心波長４３３ｎｍでの位相が１.５２であって屈折率が１.３５４であり、中心波長６９４ｎｍでの位相が０.７２であって屈折率が１.３４９であり、長波長で屈折率が低くなっている。
　図２５（ａ）は、中心波長４３３ｎｍでの屈折率分布画像である。図２５（ｂ）は、中心波長６９４ｎｍでの屈折率分布画像である。これらの屈折率分布画像は、位相画像の各画素値を波長に基づいて補正することで作成したものである。図２５（ａ）の屈折率分布画像は、図２４（ａ）の位相画像に係数４.３３を乗じたものである。図２５（ｂ）の屈折率分布画像は、図２４（ｂ）の位相画像に係数６.９４を乗じたものである。各位相画像に乗じる係数の値は、その位相画像を取得した際に用いた光の波長の一定倍の値であればよい。
　図２６（ａ）は、図２５（ａ）中の破線矩形で示した範囲を拡大して示す図である。図２６（ｂ）は、図２５（ｂ）中の破線矩形で示した範囲を拡大して示す図である。図２７（ａ）は、中心波長５３０ｎｍでの強度画像である。図２７（ｂ）は、図２６（ａ）および図２６（ｂ）から作成した屈折率変動画像である。
　屈折率変動画像（図２７（ｂ））は、中心波長６９４ｎｍでの屈折率分布画像（図２６（ｂ））から中心波長４３３ｎｍでの屈折率分布画像（図２６（ａ））を減算して、背景（位相０）より位相が大きい領域（位相が正の値の領域）を抽出したものである。なお、屈折率変動画像は、２波長での屈折率分布画像の比に基づいて作成してもよいし、２波長での屈折率分布画像の差または比と閾値との比較に基づいて作成してもよい。強度画像（図２７（ａ））と屈折率変動画像（図２７（ｂ））とを比較すると、標識物質で標識された部位が屈折率の変化として検出されていることが分かる。
　図２８および図２９は、実施例１Ｂの各画像である。実施例１Ｂでは、標識物質としてブリリアントクレシルブルーを用い、屈折率分布を測定する際に中心波長４３３ｎｍ，６９４ｎｍの光を用い、評価対象の細胞としてＨｅｐＧ２細胞株を用いた。
　図２８（ａ）は、中心波長４３３ｎｍでの位相画像である。図２８（ｂ）は、中心波長６９４ｎｍでの位相画像である。図中において矢印で指し示す位置は、標識物質で標識されている部位であって、中心波長４３３ｎｍで位相が０.１４であって屈折率が１.３３５であり、中心波長６９４ｎｍで位相が３.５９であって屈折率が１.４１２であり、長波長で屈折率が高くなっている。
　図２９（ａ）は、中心波長５３０ｎｍでの強度画像である。図２９（ｂ）は、屈折率変動画像である。屈折率変動画像（図２９（ｂ））は、実施例１Ａと同様にして作成したものである。強度画像（図２９（ａ））と屈折率変動画像（図２９（ｂ））とを比較すると、実施例１Ｂでも、実施例１Ａと同様に、標識物質で標識された部位が屈折率の変化として検出されていることが分かる。
　図３０および図３１は、実施例２Ａの各画像である。実施例２Ａでは、標識物質としてニューメチレンブルーを用い、屈折率分布を測定する際に中心波長５３０ｎｍ，６９４ｎｍの光を用い、評価対象の細胞としてＡ５４９細胞株を用いた。
　図３０（ａ）は、中心波長５３０ｎｍでの位相画像である。図３０（ｂ）は、中心波長６９４ｎｍでの位相画像である。図３１（ａ）は、中心波長５３０ｎｍでの強度画像である。図３１（ｂ）は、屈折率変動画像である。屈折率変動画像（図３１（ｂ））は、実施例１Ａと同様にして作成したものである。強度画像（図３１（ａ））と屈折率変動画像（図３１（ｂ））とを比較すると、実施例２Ａでも、同様に、標識物質で標識された部位が屈折率の変化として検出されていることが分かる。
　図３２および図３３は、実施例２Ｂの各画像である。実施例２Ｂでは、標識物質としてニューメチレンブルーを用い、屈折率分布を測定する際に中心波長５３０ｎｍ，６９４ｎｍの光を用い、評価対象の細胞としてＨｅｐＧ２細胞株を用いた。
　図３２（ａ）は、中心波長５３０ｎｍでの位相画像である。図３２（ｂ）は、中心波長６９４ｎｍでの位相画像である。図３３（ａ）は、中心波長５３０ｎｍでの強度画像である。図３３（ｂ）は、屈折率変動画像である。屈折率変動画像（図３３（ｂ））は、実施例１Ａと同様にして作成したものである。強度画像（図３３（ａ））と屈折率変動画像（図３３（ｂ））とを比較すると、実施例２Ｂでも、同様に、標識物質で標識された部位が屈折率の変化として検出されていることが分かる。
　図３４および図３５は、実施例３の各画像である。実施例３では、標識物質として中性赤を用い、屈折率分布を測定する際に中心波長４３３ｎｍ，６２０ｎｍの光を用い、評価対象の細胞としてＨｅｐＧ２細胞株を用いた。
　図３４（ａ）は、中心波長４３３ｎｍでの位相画像である。図３４（ｂ）は、中心波長６２０ｎｍでの位相画像である。図３５（ａ）は、中心波長５３０ｎｍでの強度画像である。図３５（ｂ）は、屈折率変動画像である。屈折率変動画像（図３５（ｂ））は、実施例１Ａと同様にして作成したものである。強度画像（図３５（ａ））と屈折率変動画像（図３５（ｂ））とを比較すると、実施例３でも、同様に、標識物質で標識された部位が屈折率の変化として検出されていることが分かる。
　図３６および図３７は、比較例３の各画像である。比較例３では、実施例３に対し、標識物質を用いなかった点で相違する。
　図３６（ａ）は、中心波長４３３ｎｍでの位相画像である。図３６（ｂ）は、中心波長６２０ｎｍでの位相画像である。図３７（ａ）は、中心波長５３０ｎｍでの強度画像である。図３７（ｂ）は、屈折率変動画像である。屈折率変動画像（図３７（ｂ））は、実施例１Ａと同様にして作成したものである。標識物質を用いなかった比較例３では、殆ど細胞内の特定部位を検出することができなかった。
　図３８および図３９は、実施例４の各画像である。実施例４では、標識物質としてヤヌスグリーンＢを用い、屈折率分布を測定する際に中心波長５３０ｎｍ，６９４ｎｍの光を用い、評価対象の細胞としてＡ５４９細胞株を用いた。
　図３８（ａ）は、中心波長５３０ｎｍでの位相画像である。図３８（ｂ）は、中心波長６９４ｎｍでの位相画像である。図３９（ａ）は、中心波長５３０ｎｍでの強度画像である。図３９（ｂ）は、屈折率変動画像である。屈折率変動画像（図３９（ｂ））は、実施例１Ａと同様にして作成したものである。強度画像（図３９（ａ））と屈折率変動画像（図３９（ｂ））とを比較すると、実施例４でも、同様に、標識物質で標識された部位が屈折率の変化として検出されていることが分かる。
　図４０および図４１は、比較例４の各画像である。比較例４では、実施例４に対し、標識物質を用いなかった点で相違する。
　図４０（ａ）は、中心波長５３０ｎｍでの位相画像である。図４０（ｂ）は、中心波長６９４ｎｍでの位相画像である。図４１（ａ）は、中心波長５３０ｎｍでの強度画像である。図４１（ｂ）は、屈折率変動画像である。屈折率変動画像（図４１（ｂ））は、実施例１Ａと同様にして作成したものである。標識物質を用いなかった比較例４では、殆ど細胞内の特定部位を検出することができなかった。
　以上のとおり、本実施形態では、第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で細胞の特定部位を標識し、その細胞の屈折率分布を第１波長および第２波長それぞれにおいて取得し、第１波長および前記第２波長それぞれにおける屈折率分布を比較することで、細胞における特定部位の分布を評価する。本実施形態によれば、未知の細胞であっても該細胞を容易に評価することができる。
　強度画像（図２７（ａ）等）でも、標識物質で標識された特定部位の分布が分かる。しかし、強度画像には、標識物質による光吸収だけでなく、光散乱等の影響が現れる。強度画像には、吸収の影響だけでなく、屈折率が異なる構造の境界での光の屈折や散乱により、光強度が下がって観察される場合がある（これは、非標識の細胞でも観察される場合がある）。
　光散乱等の影響はノイズの原因になりえる。１波長での強度画像でもこれが問題になるが、複数波長での強度画像を比較する際にも問題になる可能性がある。一方、屈折率分布画像では、屈折率を定量することができるので、この問題は回避しやすい。したがって、屈折率変動画像は、光散乱等に起因するノイズが小さいものとなりえる。
　光吸収性標識物質を用いる場合、光吸収がある波長領域（図５（ａ）の特定波長λ０付近の領域）では、光が透過しにくいことから、屈折率分布画像の取得に影響する可能性がある。ＯＤＴにより複数方向からの光を使用して細胞塊を観察する際には、この影響が大きくなると予想され、深部観察が難しくなる。
　しかし、光吸収性標識物質の屈折率スペクトルにおいて、特定波長λ０に対して短波長側において屈折率が極小となる波長付近、または、特定波長λ０に対して長波長側において屈折率が極大となる波長付近に、第１波長・第２波長を設定することにより（図５（ｂ），図６）、光吸収が小さい波長の光により位相画像または屈折率分布画像を取得することができる。
　また、光吸収性標識物質を用いる場合、互いに比較的近い第１波長・第２波長を設定することができる（図５（ｂ），図６）。互いに比較的近い第１波長・第２波長での屈折率分布画像を比較することで、屈折率変動画像の抽出が容易になり、ＳＮ比の高いデータが得られる。
　複数波長の光で屈折率分布画像を取得する前提で、波長の強度画像を利用して対象部位を特定する場合（以下「ケース１」という。）と、複数波長の屈折率分布画像を利用して対象部位を特定する場合（以下「ケース２」という。）と、を比較すると、次のようなことが言える。
　ケース２では、低屈折率のピーク（第１波長）と高屈折率のピーク（第２波長）とが互いに比較的近くにある場合、第１波長・第２波長での屈折率分布画像を比較することで、屈折率変動領域の抽出が容易になり、ＳＮ比の高いデータが得やすくなるメリットがある。これは、ケース２では、非標識部位の屈折率分布の変化が小さい状況で、標識部位の屈折率分布の大きな変化をとらえることができることによるものであり、ケース１に対して優位な点である。
　一方、標識物質の吸光スペクトルや屈折率スペクトルの関係上、低屈折率のピーク（第１波長）と高屈折率のピーク（第２波長）とが互いに比較的近くにない場合でも、ケース２では、屈折率スペクトルにおいて低屈折率のピークが明確であることが多いと想定されるため、第１波長・第２波長での屈折率分布画像を比較することで、屈折率変動領域の抽出が容易になり、ＳＮ比の高いデータが得やすくなる。
　このようなケース１に対するケース２の優位性は、ケース１では光吸収スペクトルの吸光度の高いピークに着目する必要があるのに対して、ケース２では屈折率スペクトルの高屈折率のピークに加えて低屈折率のピークにも着目することに因るものである。
　蛍光画像を取得する場合には、高パワーの励起光を照射する必要があることから光毒性や光退色の問題が生じる。これに対して、本実施形態では、高パワーの光を照射する必要がないので光毒性や光退色の問題が少なく、長時間に亘る細胞のタイムラプス観察が可能である。また、光毒性の問題が少ないことから、細胞本来の機能を保った状態での観察が可能であり、また、光退色の問題が少ないことから、定量的な観察が可能である。
　また、光吸収性標識物質を用いる場合、光吸収がある波長領域（図５（ａ）における特定波長λ０付近の領域）の光を照射すると、吸収された光のエネルギーが熱などに変換されて細胞毒性を引き起こす場合がある。これに対して、本実施形態では、光吸収がある波長領域を避けることができるため（図５（ｂ））、細胞毒性を低減し、細胞本来の機能を保った状態での観察が可能である。
　また、１台の装置で、屈折率分布画像や屈折率変動画像を取得することができるとともに、強度画像（図２７（ａ）等）をも取得することができる。細胞の屈折率変動画像と強度画像とを比較することで、生物学的に意義のある観察が可能となる。
　さらに、本実施形態では、屈折率の絶対値を求めるため、強度画像のみを取得する場合に比べて、定量性に優れた観察が可能となる。また、本実施形態によって、標識のために好適な光吸収性標識物質をスクリーニングしてもよい。
　細胞評価方法および細胞評価装置は、上述した実施形態及び構成例に限定されるものではなく、種々の変形が可能である。
　上記実施形態による細胞評価方法は、第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で細胞の特定部位を標識する標識ステップと、標識ステップで特定部位が標識された細胞の屈折率分布を第１波長および第２波長それぞれにおいて取得する屈折率分布取得ステップと、第１波長および第２波長それぞれにおける屈折率分布を比較することで細胞における特定部位の分布を評価する解析ステップと、を備える。
　上記の細胞評価方法では、標識ステップにおいて、第１波長および第２波長のうちの長波長における屈折率が短波長における屈折率より高い標識物質で細胞の特定部位を標識する構成としてもよい。また、標識ステップにおいて、第１波長および第２波長のうちの長波長における屈折率が短波長における屈折率より低い標識物質で細胞の特定部位を標識する構成としてもよい。
　上記の細胞評価方法では、標識ステップにおいて、特定波長において吸収ピークを有する光吸収性標識物質で細胞の特定部位を標識し、屈折率分布取得ステップにおいて、標識ステップで特定部位が標識された細胞の屈折率分布を、特定波長に対して短波長側または長波長側にある第１波長および第２波長それぞれにおいて取得する構成としてもよい。
　上記の細胞評価方法では、標識ステップにおいて、特定波長において吸収ピークを有する光吸収性標識物質で細胞の特定部位を標識し、屈折率分布取得ステップにおいて、標識ステップで特定部位が標識された細胞の屈折率分布を、第１波長および第２波長のうちいずれか一方を特定波長に対して短波長側とし、他方を特定波長に対して長波長側として、第１波長および第２波長のそれぞれにおいて取得する構成としてもよい。
　上記の細胞評価方法では、標識ステップにおいて、色素タンパク質を指令するＤＮＡ配列が任意タンパク質を指令するＤＮＡ配列のＮ端末またはＣ端末に接続された細胞内発現用コンストラクトを細胞に導入して、この細胞内発現用コンストラクトを標識物質として用いて特定部位を標識する構成としてもよい。
　上記の細胞評価方法では、解析ステップにおいて、第１波長および第２波長それぞれにおける屈折率分布の比または差に基づいて細胞における特定部位の分布を評価する構成としてもよい。また、解析ステップにおいて、第１波長および第２波長それぞれにおける屈折率分布の比較の結果を細胞の強度画像または蛍光画像に重ねることで、細胞の特定部位に関する解析を行う構成としてもよい。
　上記の細胞評価方法では、標識ステップにおいて、第３波長および第４波長それぞれにおける屈折率が互いに異なる他の標識物質で細胞の他の特定部位を標識し、屈折率分布取得ステップにおいて、標識ステップで標識された細胞の屈折率分布を第３波長および第４波長それぞれにおいて取得し、解析ステップにおいて、第３波長および第４波長それぞれにおける屈折率分布を比較することで細胞における他の特定部位の分布を評価する構成としてもよい。また、この場合、第１波長または第２波長は、第３波長または第４波長と等しい構成としてもよい。
　上記実施形態による細胞評価装置は、第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で特定部位が標識された細胞の屈折率分布を第１波長および第２波長それぞれにおいて取得する屈折率分布取得部と、第１波長および第２波長それぞれにおける屈折率分布を比較することで細胞における特定部位の分布を評価する解析部と、を備える。
　上記実施形態によるプログラムは、第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で特定部位が標識された細胞の屈折率分布を第１波長および第２波長それぞれにおいて取得する屈折率分布取得ステップと、第１波長および第２波長それぞれにおける屈折率分布を比較することで細胞における特定部位の分布を評価する解析ステップと、をコンピュータに実行させるためのプログラムである。
　実施形態は、未知の細胞であっても該細胞を容易に評価することができる細胞評価方法および細胞評価装置として利用可能である。
　１，２Ａ，２Ｂ…細胞評価装置、１０…屈折率分布取得部、１１…光源、１２…レンズ、１３…偏光子、１４…第１プリズム、１５…コンデンサレンズ、２１…対物レンズ、２２…第２プリズム、２４…レンズ、２７…ミラー、３０…光学モジュール、３１…レンズ、３３…１／４波長板、３４…レンズ、３５…偏光カメラ、３６～３８…ミラー、３９…演算部、４０…解析部、５０Ａ，５０Ｂ…屈折率分布取得部、５１，５２…レーザ光源、５３，５４…シャッタ、５５…ミラー、５６…ダイクロイックミラー、５７…ビームスプリッタ、６１，６２…レンズ、６３，６４…光ファイバ、７１…レンズ、７２…ミラー、７３…レンズ、７４…コンデンサレンズ、７５…対物レンズ、８１…ビームスプリッタ、８２…レンズ、８３Ａ，８３Ｂ…カメラ、８４…演算部、９０…解析部。

Claims (12)

  1. 　第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で細胞の特定部位を標識する標識ステップと、
    　前記標識ステップで前記特定部位が標識された前記細胞の屈折率分布を前記第１波長および前記第２波長それぞれにおいて取得する屈折率分布取得ステップと、
    　前記第１波長および前記第２波長それぞれにおける屈折率分布を比較することで前記細胞における前記特定部位の分布を評価する解析ステップと、
    を備える、細胞評価方法。
  2. 　前記標識ステップにおいて、前記第１波長および前記第２波長のうちの長波長における屈折率が短波長における屈折率より高い標識物質で前記細胞の前記特定部位を標識する、請求項１に記載の細胞評価方法。
  3. 　前記標識ステップにおいて、前記第１波長および前記第２波長のうちの長波長における屈折率が短波長における屈折率より低い標識物質で前記細胞の前記特定部位を標識する、請求項１に記載の細胞評価方法。
  4. 　前記標識ステップにおいて、特定波長において吸収ピークを有する光吸収性標識物質で前記細胞の特定部位を標識し、
    　前記屈折率分布取得ステップにおいて、前記標識ステップで前記特定部位が標識された前記細胞の屈折率分布を、前記特定波長に対して短波長側または長波長側にある前記第１波長および前記第２波長それぞれにおいて取得する、請求項１～３の何れか１項に記載の細胞評価方法。
  5. 　前記標識ステップにおいて、特定波長において吸収ピークを有する光吸収性標識物質で前記細胞の特定部位を標識し、
    　前記屈折率分布取得ステップにおいて、前記標識ステップで前記特定部位が標識された前記細胞の屈折率分布を、前記第１波長および前記第２波長のうちいずれか一方を前記特定波長に対して短波長側とし、他方を前記特定波長に対して長波長側として、前記第１波長および前記第２波長のそれぞれにおいて取得する、請求項１～３の何れか１項に記載の細胞評価方法。
  6. 　前記標識ステップにおいて、色素タンパク質を指令するＤＮＡ配列が任意タンパク質を指令するＤＮＡ配列のＮ端末またはＣ端末に接続された細胞内発現用コンストラクトを前記細胞に導入して、この細胞内発現用コンストラクトを前記標識物質として用いて前記特定部位を標識する、請求項１～３の何れか１項に記載の細胞評価方法。
  7. 　前記解析ステップにおいて、前記第１波長および前記第２波長それぞれにおける屈折率分布の比または差に基づいて前記細胞における前記特定部位の分布を評価する、請求項１～６の何れか１項に記載の細胞評価方法。
  8. 　前記解析ステップにおいて、前記第１波長および前記第２波長それぞれにおける屈折率分布の比較の結果を前記細胞の強度画像または蛍光画像に重ねることで、前記細胞の前記特定部位に関する解析を行う、請求項１～７の何れか１項に記載の細胞評価方法。
  9. 　前記標識ステップにおいて、第３波長および第４波長それぞれにおける屈折率が互いに異なる他の標識物質で前記細胞の他の特定部位を標識し、
    　前記屈折率分布取得ステップにおいて、前記標識ステップで標識された前記細胞の屈折率分布を前記第３波長および前記第４波長それぞれにおいて取得し、
    　前記解析ステップにおいて、前記第３波長および前記第４波長それぞれにおける屈折率分布を比較することで前記細胞における前記他の特定部位の分布を評価する、請求項１～８の何れか１項に記載の細胞評価方法。
  10. 　前記第１波長または前記第２波長は、前記第３波長または前記第４波長と等しい、請求項９に記載の細胞評価方法。
  11. 　第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で特定部位が標識された細胞の屈折率分布を前記第１波長および前記第２波長それぞれにおいて取得する屈折率分布取得部と、
    　前記第１波長および前記第２波長それぞれにおける屈折率分布を比較することで前記細胞における前記特定部位の分布を評価する解析部と、
    を備える、細胞評価装置。
  12. 　第１波長および第２波長それぞれにおける屈折率が互いに異なる標識物質で特定部位が標識された細胞の屈折率分布を前記第１波長および前記第２波長それぞれにおいて取得する屈折率分布取得ステップと、
    　前記第１波長および前記第２波長それぞれにおける屈折率分布を比較することで前記細胞における前記特定部位の分布を評価する解析ステップと、
    をコンピュータに実行させるためのプログラム。
PCT/JP2023/010411 2022-03-22 2023-03-16 細胞評価方法および細胞評価装置 WO2023182169A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2022-045496 2022-03-22
JP2022045496A JP7390423B2 (ja) 2022-03-22 2022-03-22 細胞評価方法および細胞評価装置

Publications (1)

Publication Number Publication Date
WO2023182169A1 true WO2023182169A1 (ja) 2023-09-28

Family

ID=88101525

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2023/010411 WO2023182169A1 (ja) 2022-03-22 2023-03-16 細胞評価方法および細胞評価装置

Country Status (2)

Country Link
JP (2) JP7390423B2 (ja)
WO (1) WO2023182169A1 (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011112548A (ja) * 2009-11-27 2011-06-09 Sony Corp 生体サンプル解析方法、生体サンプル解析装置及び生体サンプル解析プログラム
JP2012055267A (ja) * 2010-09-10 2012-03-22 Kanagawa Acad Of Sci & Technol 細胞観察装置および細胞観察方法
JP5745919B2 (ja) 2011-04-28 2015-07-08 浜松ホトニクス株式会社 細胞解析方法、細胞解析装置、および細胞解析プログラム
WO2020013325A1 (ja) * 2018-07-13 2020-01-16 国立大学法人東京大学 画像生成装置及び画像生成方法
WO2022054305A1 (ja) * 2020-09-11 2022-03-17 浜松ホトニクス株式会社 観察装置および観察方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011112548A (ja) * 2009-11-27 2011-06-09 Sony Corp 生体サンプル解析方法、生体サンプル解析装置及び生体サンプル解析プログラム
JP2012055267A (ja) * 2010-09-10 2012-03-22 Kanagawa Acad Of Sci & Technol 細胞観察装置および細胞観察方法
JP5745919B2 (ja) 2011-04-28 2015-07-08 浜松ホトニクス株式会社 細胞解析方法、細胞解析装置、および細胞解析プログラム
WO2020013325A1 (ja) * 2018-07-13 2020-01-16 国立大学法人東京大学 画像生成装置及び画像生成方法
WO2022054305A1 (ja) * 2020-09-11 2022-03-17 浜松ホトニクス株式会社 観察装置および観察方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
BACZEWSKA MARIA, EDER KAI, KETELHUT STEFFI, KEMPER BJÖRN, KUJAWIŃSKA MAŁGORZATA: "Refractive Index Changes of Cells and Cellular Compartments Upon Paraformaldehyde Fixation Acquired by Tomographic Phase Microscopy", CYTOMETRY A, WILEY-LISS, HOBOKEN, USA, vol. 99, no. 4, 1 April 2021 (2021-04-01), Hoboken, USA, pages 388 - 398, XP093095101, ISSN: 1552-4922, DOI: 10.1002/cyto.a.24229 *
D. FU ET AL.: "Quantitative dispersion microscopy", BIOMED. OPT. EXPRESS, vol. 1, 2010, pages 347 - 353
O. YASUHIKOK. TAKEUCHIH. YAMADAY. UEDA: "Multiple-scattering suppressive refractive index tomography for the label-free quantitative assessment of multicellular spheroids", BIOMED. OPT. EXPRESS, vol. 13, 2022, pages 962 - 979
O. YASUHIKOK. TAKEUCHIH. YAMADAY. UEDA: "Single-shot quantitative phase imaging as an extension of differential interference contrast microscopy", GENES TO CELLS, vol. 26, 2021, pages 596 - 610
T. SAI ET AL.: "Designing refractive index fluids using the Kramers-Kronig relations", FARADAY DISCUSSIONS, vol. 223, 2020, pages 136 - 144

Also Published As

Publication number Publication date
JP2023139787A (ja) 2023-10-04
JP7390423B2 (ja) 2023-12-01
JP2024019188A (ja) 2024-02-08

Similar Documents

Publication Publication Date Title
Aknoun et al. Living cell dry mass measurement using quantitative phase imaging with quadriwave lateral shearing interferometry: an accuracy and sensitivity discussion
Huang et al. Linear birefringence of the retinal nerve fiber layer measured in vitro with a multispectral imaging micropolarimeter
Haseda et al. Significant correlation between refractive index and activity of mitochondria: single mitochondrion study
Strupler et al. Second harmonic imaging and scoring of collagen in fibrotic tissues
JP4670031B2 (ja) 試料内で励起および／または後方散乱を経た光ビームの光学的検出のための装置
Ghosh et al. Dynamic imaging of homo-FRET in live cells by fluorescence anisotropy microscopy
Williams et al. Full spectrum fluorescence lifetime imaging with 0.5 nm spectral and 50 ps temporal resolution
CN111328370A (zh) 一种表面缺陷检测系统及方法
Ströhl et al. Total internal reflection fluorescence anisotropy imaging microscopy: setup, calibration, and data processing for protein polymerization measurements in living cells
Bohannon et al. Refractive index imaging of cells with variable-angle near-total internal reflection (TIR) microscopy
Guo et al. Limited-angle tomographic phase microscopy utilizing confocal scanning fluorescence microscopy
US20130230864A1 (en) Method for and use of digital holographic microscopy and imaging on labelled cell samples
Hsiao et al. Spinning disk interferometric scattering confocal microscopy captures millisecond timescale dynamics of living cells
US11971357B2 (en) In situ measurement of absolute concentrations by normalized Raman imaging
Miyazaki et al. Label-free dynamic imaging of mitochondria and lysosomes within living cells via simultaneous dual-pump photothermal microscopy
WO2015004917A1 (ja) 観察方法、観察装置、細胞シート製造方法および細胞シート製造装置
Rothe et al. Angular and spectrally resolved investigation of single particles by darkfield scattering microscopy
de Kernier et al. Large field-of-view phase and fluorescence mesoscope with microscopic resolution
Blanke et al. Multiscale label-free imaging of myelin in human brain tissue with polarization-sensitive optical coherence tomography and birefringence microscopy
Yang et al. A quantitative method for measurement of HL-60 cell apoptosis based on diffraction imaging flow cytometry technique
WO2023182169A1 (ja) 細胞評価方法および細胞評価装置
Mehta et al. Sample-less calibration of the differential interference contrast microscope
Huang et al. Diattenuation and polarization preservation of retinal nerve fiber layer reflectance
Cogswell et al. Colour confocal reflection microscopy using red, green and blue lasers
Arnold et al. Development of inspection system for the detection and analysis of solid particles and oil droplets in process water of the petrochemical industry using hyperspectral imaging and fluorescence imaging

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

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2023774765

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2023774765

Country of ref document: EP

Effective date: 20241016