JPWO2019176427A1 - 判定方法 - Google Patents

判定方法 Download PDF

Info

Publication number
JPWO2019176427A1
JPWO2019176427A1 JP2020505699A JP2020505699A JPWO2019176427A1 JP WO2019176427 A1 JPWO2019176427 A1 JP WO2019176427A1 JP 2020505699 A JP2020505699 A JP 2020505699A JP 2020505699 A JP2020505699 A JP 2020505699A JP WO2019176427 A1 JPWO2019176427 A1 JP WO2019176427A1
Authority
JP
Japan
Prior art keywords
sphere
phase difference
determined
correlation
index value
Prior art date
Legal status (The legal status 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 status listed.)
Granted
Application number
JP2020505699A
Other languages
English (en)
Other versions
JP6995975B2 (ja
Inventor
崇市郎 中村
崇市郎 中村
龍介 大▲崎▼
龍介 大▲崎▼
祥 小野澤
祥 小野澤
彰洋 朝倉
彰洋 朝倉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Corp
Original Assignee
Fujifilm Corp
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 Fujifilm Corp filed Critical Fujifilm Corp
Publication of JPWO2019176427A1 publication Critical patent/JPWO2019176427A1/ja
Application granted granted Critical
Publication of JP6995975B2 publication Critical patent/JP6995975B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
    • C12M41/00—Means for regulation, monitoring, measurement or control, e.g. flow regulation
    • C12M41/30—Means for regulation, monitoring, measurement or control, e.g. flow regulation of concentration
    • C12M41/36—Means for regulation, monitoring, measurement or control, e.g. flow regulation of concentration of biomass, e.g. colony counters or by turbidity measurements
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/06—Investigating concentration of particle suspensions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1429—Signal processing
    • G01N15/1433—Signal processing using image recognition
    • 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
    • G01N21/453—Holographic interferometry
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • 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
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/06—Investigating concentration of particle suspensions
    • G01N15/075—Investigating concentration of particle suspensions by optical means
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N2015/1006—Investigating individual particles for cytology
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434—Optical arrangements
    • G01N2015/1454—Optical arrangements using phase shift or interference, e.g. for improving contrast
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N2015/1493—Particle size
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N2015/1497—Particle shape
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/06—Illumination; Optics
    • G01N2201/061—Sources
    • G01N2201/06113—Coherent sources; lasers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/06—Illumination; Optics
    • G01N2201/063—Illuminating optical parts
    • G01N2201/0633—Directed, collimated illumination
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00—Microscopes
    • G02B21/16—Microscopes adapted for ultraviolet illumination ; Fluorescence microscopes
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/0005—Adaptation of holography to specific applications
    • G03H2001/0033—Adaptation of holography to specific applications in hologrammetry for measuring or analysing
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/0005—Adaptation of holography to specific applications
    • G03H2001/005—Adaptation of holography to specific applications in microscopy, e.g. digital holographic microscope [DHM]
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • G03H2001/0875—Solving phase ambiguity, e.g. phase unwrapping
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • G03H2001/0883—Reconstruction aspect, e.g. numerical focusing
    • 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/30024—Cell structures in vitro; Tissue sections in vitro

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Organic Chemistry (AREA)
  • Wood Science & Technology (AREA)
  • Zoology (AREA)
  • Pathology (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biotechnology (AREA)
  • Quality & Reliability (AREA)
  • Medical Informatics (AREA)
  • Genetics & Genomics (AREA)
  • General Engineering & Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Dispersion Chemistry (AREA)
  • Microbiology (AREA)
  • Geometry (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Sustainable Development (AREA)
  • Biomedical Technology (AREA)
  • Signal Processing (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Apparatus Associated With Microorganisms And Enzymes (AREA)
  • Holo Graphy (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

複数の細胞の凝集体であるスフェアの状態を非破壊且つ簡便に判定することができる判定方法を提供する。複数の細胞の凝集体であるスフェアを撮像したホログラムからスフェアの位相差画像を生成し、位相差画像と前記スフェアの形状に応じた形状指標値とに基づいて、スフェアの状態を判定する。

Description

開示の技術は、複数の細胞の凝集体であるスフェアの状態を判定する判定方法に関する。
細胞の状態を評価または判定する技術として、例えば、以下の技術が知られている。特許文献１には、細胞群に対して照明光を照射し、照明光の照射によって細胞群を透過した透過光を検出して細胞群の画像を検出し、細胞群における散乱光を検出するための検出光を細胞群に対して照射し、検出した散乱光の強度の空間的な広がりまたは時間的な揺らぎに基づいて細胞群を評価する細胞評価方法が記載されている。
特許文献２には、神経細胞分化過程の細胞を撮像した第１の撮像画像を入力する画像入力部と、細胞の厚さ方向の位置において焦点の合った状態の全焦点画像を、少なくとも第１の撮像画像に基づく第１の原画像として生成する全焦点画像生成部と、第１の原画像において一定以下の輝度分布を有する領域と、第１の原画像において濃度差が一定以下の領域とで共通する領域を、分化過程に出現するロゼッタとして抽出するロゼッタ抽出部と、抽出されたロゼッタの状態を判定するロゼッタ対応判定部と、を備える細胞評価装置が記載されている。
特許文献３には、細胞厚さに基づいて、細胞の分化度を判定する判定ステップを備えることを特徴とする細胞判定方法が記載されている。
特開２０１６−００５４３７号公報 特開２０１７−０００１６３号公報 特開２０１５−１４６７４７号公報
細胞の大量生産が可能な培養手法として、細胞の凝集体であるスフェアを培地中に浮遊させた状態で培養する三次元培養法が知られている。三次元培養による細胞の製造工程においては、工程管理の容易性の観点から、細胞の品質をスフェアの状態のまま非破壊且つ簡便に評価する技術が求められる。しかしながら、現時点においては、三次元空間中にランダムに存在する様々なサイズのスフェアを評価する手法が確立されておらず、特にスフェアの内部における細胞の密度及び生存状況について直接観察することは困難である。このため、従来の二次元培養の手法を適用した評価がなされているが、培養する細胞の数の増加に伴って、評価工数が増えるため、多くの人手及び多くの時間が必要となる。従来の二次元培養の手法を適用した評価では、スフェアを単一細胞に分解したり、蛍光色素を添加したりするといった細胞の破壊を伴う処理が必要とされる。
開示の技術は、複数の細胞の凝集体であるスフェアの状態を非破壊且つ簡便に判定することを目的とする。
開示の技術に係る判定方法は、複数の細胞の凝集体であるスフェアを撮像したホログラムからスフェアの位相差画像を生成し、位相差画像とスフェアの形状に応じた形状指標値とに基づいて、スフェアの状態を判定することを含む。
開示の技術に係る判定方法によれば、スフェアの状態を非破壊且つ簡便に判定することが可能となる。
開示の技術に係る判定方法において、位相差画像と形状指標値とに基づいて、スフェアに含まれる複数の細胞の生存率、密度、均質性、未分化状態逸脱、及びスフェアの外形形状の少なくとも１つに関する判定を行ってもよい。
開示の技術に係る判定方法において、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を導出し、位相差量総和と形状指標値との相関性を用いてスフェアの状態を判定してもよい。
開示の技術に係る判定方法において、位相差量総和と形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象のスフェアについての位相差量総和と形状指標値との相関性とを比較し、判定対象のスフェアについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の程度に応じて判定対象のスフェアの状態を判定してもよい。例えば、判定対象のスフェアについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の幅が閾値を超える場合、判定対象のスフェアに含まれる細胞の生存率、密度、均質性及びスフェアの外形形状の少なくとも１つについて異常ありと判定してもよい。
開示の技術に係る判定方法において、位相差量総和と形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象の培養ロットに属する複数のスフェアの各々についての位相差量総和と形状指標値との相関性とを比較し、判定対象の培養ロットについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて判定対象の培養ロットの良否判定を行ってもよい。
開示の技術に係る判定方法において、形状指標値は、スフェアの体積、断面積、粒径及び周の長さのいずれかであってもよい。
開示の技術に係る判定方法において、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を、スフェアの体積で除算して位相差量密度を導出し、位相差量密度に基づいて、スフェアの状態を判定してもよい。
開示の技術に係る判定方法において、スフェアの状態の判定に用いる位相差画像は、ホログラムから生成可能な複数の位相差画像のうち、位相差画像を構成する複数の画素間の位相差量のばらつきが最大となる位相差画像であることが好ましい。
開示の技術によれば、複数の細胞の凝集体であるスフェアの状態を非破壊且つ簡便に判定することが可能となる。
開示の技術の実施形態に係る判定方法の実施に用いる撮像システムの構成の一例を示す図である。 開示の技術の実施形態に係る判定方法の実施に用いるホログラムの一例を示す図である。 スフェアのフーリエ変換画像の一例を示す図である。 スフェアのアンラッピング前の位相差画像の一例を示す図である。 スフェアのアンラッピング後の位相差画像の一例を示す図である。 開示の技術の実施形態に係る位相差画像の概念を示す図である。 開示の技術の実施形態に係る位相差画像の焦点合わせに関する説明図である。 開示の技術の実施形態に係るオートフォーカス処理を行うコンピュータのハードウェア構成の一例である。 開示の技術の実施形態に係るオートフォーカス処理の流れの一例を示すフローチャートである。 開示の技術の実施形態に係る、スフェアの位相差画像における焦点位置と位相差量のばらつきとの関係の一例を示すグラフである。 開示の技術の実施形態に係る、スフェアの体積と位相差量総和との相関特性の一例を示すグラフである。 開示の技術の実施形態に係る、良ロット及び不良ロットの各々について導出したスフェアの粒径のヒストグラムである。 開示の技術の実施形態に係る、良ロット及び不良ロットの各々について取得した、スフェアの体積と位相差量総和との相関特性の一例を示すグラフである。 開示の技術の実施形態に係る、良ロット及び不良ロットの各々について取得した、位相差量密度のヒストグラムである。 開示の技術の実施形態に係る、細胞の生存率と位相差量密度との相関特性の一例を示すグラフである。 Ｈ２Ｏ２を添加した後、０分経過した時点におけるコロニーの蛍光顕微鏡画像である。 Ｈ２Ｏ２を添加した後、１５分経過した時点におけるコロニーの蛍光顕微鏡画像である。 Ｈ２Ｏ２を添加した後、３７分経過した時点におけるコロニーの蛍光顕微鏡画像である。 Ｈ２Ｏ２を添加した後、０分経過した時点におけるコロニーの位相差画像である。 Ｈ２Ｏ２を添加した後、１５分経過した時点におけるコロニーの位相差画像である。 Ｈ２Ｏ２を添加した後、３７分経過した時点におけるコロニーの位相差画像である。 開示の技術の実施形態に係る、ネクローシスに誘導したコロニー、アポトーシスに誘導したコロニー及び健全なコロニーのそれぞれの位相差量総和の相対値を示すグラフである。 開示の技術の実施形態に係る、未分化状態を維持したスフェア及び未分化状態を逸脱したスフェアのそれぞれの位相差量総和の平均値を示すグラフである。
以下、本発明の実施形態について図面を参照しつつ説明する。尚、各図面において、実質的に同一又は等価な構成要素又は部分には同一の参照符号を付している。
開示の技術の実施形態に係る判定方法は、複数の細胞の凝集体であるスフェアを撮像したホログラムからスフェアの位相差画像を生成し、位相差画像とスフェアの形状に応じた形状指標値とに基づいて、スフェアの状態を判定する、というものである。
図１は、開示の技術の実施形態に係る判定方法の実施に用いる撮像システム１の構成の一例を示す図である。撮像システム１は、公知のデジタルホログラフィ技術を用いてスフェアのホログラムを取得するためのホログラム光学系１０と、スフェアの蛍光顕微鏡観察を行うための蛍光顕微鏡光学系３０と、を含んで構成されている。
デジタルホログラフィ技術は、物体を透過または反射した物体光と、物体光に対してコヒーレントである参照光との干渉によって生じる像をイメージセンサーを用いて撮像し、撮像によって得られた画像について、光伝搬に基づく数値計算を実施することにより、物体からの光波の波面を復元する技術である。デジタルホログラフィ技術によれば、物体の位相分布を定量化し、また、焦点位置を機械的に移動させることなく、物体の三次元情報を取得することができる。
ホログラム光学系１０は、レーザ光源１１、ビームスプリッタ１２及び１８、コリメートレンズ１３及び２１、対物レンズ１５、結像レンズ１７並びにＣＭＯＳ（Complementary Metal Oxide Semiconductor）カメラ１９を含んで構成されている。サンプルステージにセットされる試料１４としてのスフェアは、コリメートレンズ１３と対物レンズ１５との間に配置される。
レーザ光源１１には、例えば波長６３２．８ｎｍのＨｅＮｅレーザを用いることが可能である。レーザ光源１１から出射されたレーザ光は、ビームスプリッタ１２により、２つのレーザ光に分割される。２つのレーザ光の一方は物体光とされ、他方は参照光とされる。物体光は、コリメートレンズ１３によって平行光とされた後、サンプルステージにセットされた試料１４であるスフェアに照射される。スフェアを透過した物体光による像は、対物レンズ１５によって拡大される。対物レンズ１５を透過した物体光は、結像レンズ１７によって再び平行光とされた後、ビームスプリッタ１８を介してＣＭＯＳカメラ１９の撮像面に結像される。一方、参照光は、光ファイバ２０によってコリメートレンズ２１の手前まで導かれる。光ファイバ２０から出射した参照光は、コリメートレンズ２１によって平行光とされ、ビームスプリッタ１８を介してＣＭＯＳカメラ１９の撮像面に入射する。物体光と参照光との干渉によって生じるホログラムが、ＣＭＯＳカメラ１９によって記録される。なお、ＣＭＯＳカメラ１９の撮像面に入射する物体光と参照光の光軸方向が互いに異なったオフアキシャル光学系が構成されていてもよい。
一方、蛍光顕微鏡光学系３０は、励起光源３１、励起用フィルタ３２、ＯＮ／ＯＦＦ切り替えミラー３３、ダイクロイックミラー３４、対物レンズ１５、分光フィルタ３５及びＣＭＯＳカメラ３６を含んで構成されている。
励起光源３１には、例えば水銀ランプを用いることが可能である。励起用フィルタ３２は、励起光源３１から出射された励起光のうち、例えば４５０ｎｍ〜４９０ｎｍの波長域の光を透過させる。励起用フィルタ３２を透過した励起光は、ＯＮ／ＯＦＦ切り替えミラー３３、ダイクロイックミラー３４及び対物レンズ１５を介してサンプルステージにセットされた試料１４であるスフェアに照射される。スフェアに添加された蛍光色素に励起光が照射されることにより、蛍光色素から発せられる蛍光は、対物レンズ１５、ダイクロイックミラー３４、ＯＮ／ＯＦＦ切り替えミラー３３及び分光フィルタ３５を介してＣＭＯＳカメラ３６の撮像面に入射する。分光フィルタ３５は、蛍光色素から発せられた蛍光のうち、例えば５１５ｎｍ以上の波長域の光を透過させる。
本実施形態に係る撮像システム１によれば、ホログラム光学系１０によるホログラム観察と、蛍光顕微鏡光学系３０による蛍光顕微観察とを並行して行うことが可能である。なお、上記した撮像システム１の構成は、一例にすぎず、上記の構成に限定されない。開示の技術に係る判定方法の実施には、デジタルホログラム技術を用いてホログラムを取得することができる、あらゆる撮像システムを利用することが可能である。
以下に、撮像システム１を用いて取得したスフェアのホログラムから、スフェアの位相差画像を取得する方法の一例について説明する。
はじめに、ＣＭＯＳカメラ１９によって取得した図２Ａに例示するホログラムを、二次元フーリエ変換することにより、物体光のみの複素振幅成分を抜き出す。図２Ｂは、この処理によって得られるスフェアのフーリエ変換画像の一例である。
次に、例えば角スペクトル法を適用して任意の空間位置のスフェアの位相を示す画像を復元する。図２Ｃは、この処理によって得られるスフェアのアンラッピング前の位相差画像の一例である。この時点におけるスフェアの位相は、０〜２πの値に畳みこまれている。そこで、例えば、Unweighted Least Squares（重みなし最小２乗法）またはFlynn's Algorithm（フリンのアルゴリズム）などの位相接続（アンラッピング）手法を適用して２π以上の部分も接合していくことにより、図２Ｄに例示するような最終的なスフェアの位相差画像を得ることができる。なお、アンラッピングの手法は数多く提案されており、位相不整合を生じない適切な手法を適宜選択すればよい。
図３は、位相差画像ＩＰの概念を示す図である。図３の下段には、位相差画像ＩＰの各画素ｋにおける位相差量を３次元表示する。図３の上段には、位相差画像ＩＰの各画素ｋにおける位相差量を平面上にグレースケールで示す。
ここで、位相差画像ＩＰにおける位相差量θは、位相差画像ＩＰの同一焦点面内に存在するバックグラウンド（スフェアの存在しない領域）の位相をθＢとし、スフェアの存在する領域の位相をθＳとしたとき、下記の（１）式によって表わされる。また、本明細書中における「位相」という用語は、光を電磁波とみなした場合の電場振幅の位相であり、より一般的な意味で使用される。
また、位相差画像ＩＰの各画素ｋにおける位相差量θｋは、下記（２）式によって表わすことができる。但し、ｎｋは位相差画像ＩＰの各画素ｋに対応する部位におけるスフェアの屈折率であり、ｄｋは位相差画像ＩＰの各画素ｋに対応する部位におけるスフェアの厚さであり、λはホログラム光学系１０における物体光の波長である。
スフェアの位相差画像は、スフェアを透過した物体光の光路長分布を示した画像である。スフェア内における光路長は、スフェアの屈折率とスフェアの厚さの積に相当することから、スフェアの位相差画像は、（２）式にも示されているように、スフェアの屈折率及び厚さ（形状）の情報を含んでいる。
スフェアに対して焦点が合っていない位相差画像からは、回折による広がりの影響によりスフェアの実態に合致した正確な情報が得られない。従って、ＣＭＯＳカメラ１９によって取得したホログラムから位相差画像を取得する際に、スフェアに焦点を合わせることが好ましい。ここで、「スフェアに焦点を合わせる」とは、球形状のスフェアの中央付近においてスライスした位相差画像を得ることを意味する。スフェアに焦点が合った位相差画像を用いてスフェアの状態を判定することにより、より正確な判定結果を得ることができる。なお、ユーザが、取得したスフェアの位相差画像に基づいて、スフェアの状態を判定してもよい。または、判定条件を学習したコンピュータが、取得したスフェアの位相差画像に基づいて、スフェアの状態を判定してもよい。
位相差画像の焦点合わせは、人手によらず、自動化することが好ましい。焦点合わせを自動化することにより、作業者による任意性を排除し、更に、処理時間の短縮を図ることができる。本発明者らは、以下に説明する自動化可能な焦点合わせの手法を見出した。
図４の左のグラフは、スフェアの位相差画像における平面方向の位置と位相差量との関係の一例を示すグラフであり、実線がスフェアに焦点が合っている状態に対応し、点線がスフェアに焦点が合っていない状態に対応する。スフェアに焦点が合っている場合、位相差画像における特定の位置に急峻なピークが現れる。一方、スフェアに焦点が合っていない場合、焦点が合っている場合と比較してピークが低く且つなだらかになる。
図４の右のグラフは、スフェアの位相差画像における位相差量のヒストグラムの一例であり、実線がスフェアに焦点が合っている状態に対応し、点線がスフェアに焦点が合っていない状態に対応する。スフェアに焦点が合っている場合、カーブの幅ｗ（位相差量のばらつき）は、相対的に大きくなり、スフェアに焦点が合っていない場合、カーブの幅ｗ（位相差量のばらつき）は、相対的に小さくなる。
従って、互いに異なる焦点位置（スライス位置）毎にスフェアの位相差画像を取得し、取得した位相差画像の各々について、位相差量のヒストグラムにおけるカーブの幅ｗ（位相差量のばらつき）を求め、求めた幅ｗのうち、最大の幅ｗを有する位相差画像を、スフェアに焦点が合った位相差画像として抽出することにより焦点合わせを実現できる。
上記の焦点合わせは、コンピュータを用いて自動化することが可能である。図５は、上記の焦点合わせを自動で行うオートフォーカス処理を行うコンピュータ５００のハードウェア構成の一例である。
コンピュータ５００は、ＣＰＵ（Central Processing Unit）５０１、一時記憶領域としての主記憶装置５０２、不揮発性の補助記憶装置５０３、ＣＭＯＳカメラ１９との間で通信を行うための通信Ｉ／Ｆ（InterFace）５０４、及び液晶ディスプレイ等の表示部５０５を含んで構成されている。ＣＰＵ５０１、主記憶装置５０２、補助記憶装置５０３、通信Ｉ／Ｆ５０４、及び表示部５０５は、それぞれ、バス５０７に接続されている。補助記憶装置５０３には、上記のオートフォーカス処理の手順を記述したオートフォーカスプログラム５０６が格納されている。コンピュータ５００は、ＣＰＵ５０１がオートフォーカスプログラム５０６を実行することにより、オートフォーカス処理を行う。
図６は、コンピュータ５００において実施されるオートフォーカス処理の流れの一例を示すフローチャートである。
ステップＳ１において、ＣＰＵ５０１は、ＣＭＯＳカメラ１９からスフェアのホログラムを取得する。
ステップＳ２において、ＣＰＵ５０１は、取得したホログラムから、焦点位置（スライス位置）が互いに異なる複数の位相差画像を生成する。
ステップＳ３において、ＣＰＵ５０１は、焦点位置（スライス位置）毎の位相差画像の各々について、位相差量のばらつきを導出する。ＣＰＵ５０１は、例えば、位相差画像における位相差量の最大値と最小値との差分を、当該位相差画像における位相差量のばらつきとして導出してもよい。
ステップＳ４において、ＣＰＵ５０１は、焦点位置（スライス位置）が互いに異なる複数の位相差画像のうち、ステップＳ３にて導出した位相差量のばらつきが最大となる位相差画像を、焦点が合った位相差画像として抽出する。
図７は、スフェアの位相差画像における焦点位置（スライス位置）と位相差量のばらつきとの関係の一例を示したグラフである。図７には、焦点位置が、−４００μｍ、−２００μｍ、０μｍ、＋２００μｍ、及び＋４００μｍに対応するスフェアの位相差画像、並びにグラフが例示されている。なお、図７では、位相差量のばらつきが最大となる焦点位置を０μｍとしている。上記のオートフォーカス処理によれば、位相差量のばらつきが最大となる焦点位置０μｍに対応する位相差画像が、焦点が合った位相差画像として抽出される。位相差量のばらつきが最大となる焦点位置０μｍに対応する位相差画像において、スフェアの輪郭が最も鮮明となる。
開示の技術の実施形態に係る判定方法は、スフェアの位相差画像と、スフェアの形状指標値とに基づいてスフェアの状態を判定することを含む。すなわち、スフェアの位相差画像を、スフェアの形状指標値と関連付けて分析することにより、スフェアの状態が判定される。判定されるスフェアの状態には、例えば、スフェアに含まれる複数の細胞の生存率、密度、均質性、未分化状態逸脱、及びスフェアの外形形状などが含まれる。
スフェアの形状指標値として、例えば、スフェアの体積、断面積、粒径及び周の長さなどを用いることができる。スフェアの粒径及び周の長さは、スフェアに焦点を合わせた位相差画像（すなわち、スフェアの中央付近においてスライスした位相差画像）におけるスフェアの像から直接的に求めることが可能である。スフェアの断面積は、例えば、スフェアの粒径を直径とする円の断面積として導出することが可能である。スフェアの体積は、例えばスフェアの粒径を直径とする球の体積として導出することが可能である。
開示の技術の実施形態に係る判定方法は、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和θＡを導出し、位相差量総和θＡと形状指標値との相関性を用いてスフェアの状態を判定することを含み得る。
位相差量総和θＡは、下記の（３）式によって表わされる。但し、ｓは位相差画像の各画素ｋの面積であり、ｖｋは位相差画像の各画素ｋに対応する部位におけるスフェアの体積である。（３）式に示されるように、位相差量総和θＡは、スフェアの位相差画像の画素毎の位相差量θｋを、全ての画素ｋについて積算した値に相当する。なお、（３）式におけるｄｋは画素ｋに投影されたスフェア部分の厚みを示し、ｎｋはバックグラウンドである培養液とスフェア内部との屈折率差を表す。（３）式では、ｖｋ＝ｄｋ・sを用いた。ここで、（３）式により、位相差量総和θＡの単位は面積のスケールとなって、例えば[μｍ２]であるが、イメージセンサー間で比較を行わない場合には単に１ピクセルあたりの画素毎の位相差量θｋの和として位相差量総和θＡの単位を[pixel]とし、つまりs ＝ 1 [pixel]として扱ってもよい。
開示の技術の実施形態に係る判定方法は、位相差量総和θＡとスフェアの形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象のスフェアについての位相差量総和θＡと形状指標値との相関性とを比較し、判定対象のスフェアについての相関性の、基準相関トレンドラインからの乖離の程度に応じて当該判定対象のスフェアの状態を判定することを含み得る。この場合において、判定対象のスフェアについての相関性の、基準相関トレンドラインからの乖離の幅が閾値を超える場合、当該判定対象のスフェアに含まれる細胞の生存率、密度、均質性及び前記スフェアの外形形状の少なくとも１つについて異常ありと判定してもよい。
例えば、一方の軸に位相差量総和θＡをとり、他方の軸に形状指標値をとったグラフ上にプロットされた、判定対象のスフェアについてのプロットの、基準相関トレンドラインからの乖離の程度に応じて、当該判定対象のスフェアの状態を判定してもよい。この場合において、上記プロットの、基準相関トレンドラインからの乖離の幅が閾値を超える場合、当該判定対象のスフェアに含まれる細胞の生存率、密度、均質性及びスフェアの外形形状の少なくとも１つについて異常ありと判定してもよい。
開示の技術の実施形態に係る判定方法は、基準相関トレンドラインと、判定対象の培養ロットに属する複数のスフェアの各々についての位相差量総和θＡと形状指標値との相関性とを比較し、判定対象の培養ロットについての相関性の、基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて当該判定対象の培養ロットの良否判定を行うことを含み得る。
例えば、一方の軸に位相差量総和θＡをとり、他方の軸に形状指標値をとったグラフ上にプロットされた、判定対象の培養ロットに属する複数のスフェアの各々についてのプロットの、基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて当該培養ロットの良否判定を行ってもよい。
例えば、健全なスフェア（標準サンプル）について予め取得した、位相差量総和θＡとスフェアの形状指標値との相関ラインを、基準相関トレンドラインとして用いることが可能である。また、判定対象のスフェアについての各プロットから、例えば最小二乗法を用いて導出される回帰線を、基準相関トレンドラインとして用いてもよい。
開示の技術の実施形態に係る判定方法は、位相差量総和θＡを当該スフェアの体積で除算して位相差量密度ＤＰを導出し、導出した位相差量密度ＤＰに基づいて、当該スフェアの状態を判定することを含み得る。
位相差量密度ＤＰは、下記の（４）式によって表わされる。但し、Ｖはスフェアの体積である。（４）式に示されるように、位相差量密度ＤＰは、位相差量総和θＡを当該スフェアの体積Ｖで除算した値に相当する。健全な細胞は、その恒常性から内部の屈折率は、媒質の屈折率とは異なる一定の値を維持すると考えられる。一方、死細胞は、恒常性を喪失し、内部の屈折率が媒質の屈折率と略同じになると考えられる。従って、位相差量密度ＤＰを、細胞の状態を示す指標として用いることが可能であると考えられる。例えば、判定対象のスフェアについて取得した位相差量密度ＤＰが、閾値以上である場合に当該スフェアの状態が良好であると判定することができ、位相差量密度ＤＰが、閾値未満である場合に当該スフェアの状態が異常であると判定することができる。なお、２π／λは、定数として扱うことができるので、位相差量密度ＤＰの導出に際し、２π／λの乗算を省略してもよい。ここで、スフェアの体積平均屈折率差ＮａｖｅをＮａｖｅ=Σｎｋ ・ （ｖｋ／Ｖ）とすると、（４）式は、ＤＰ = ( ２π／ λ) ×Ｎａｖｅとなることから、位相差密度は体積平均したスフェアの屈折率差を波長の長さで規格化した値である。本明細書において、Ｖはスフェアの位相像の断面画像から球相当径を算出して求めた。より正確に楕円球とすることも可能である。
開示の技術に係る判定方法の実施例を以下に示す。
［実施例１］
三次元培養法により培養したｉＰＳ細胞（induced pluripotent stem cells）のスフェアを、図１に示す撮像システム１のサンプルステージにセットし、複数のスフェアのホログラムをＣＭＯＳカメラ１９によって撮影した。取得した各スフェアのホログラムに対してコンピュータによる数値計算を実施することにより、スフェアの中央付近においてスライスした位相差画像を取得した。
得られた各スフェアの位相差画像から、各スフェアの形状指標値として体積を導出した。更に、各スフェアの位相差画像について、（３）式によって示される位相差量総和θＡを導出した。横軸にスフェアの体積をとり、縦軸に位相差量総和θＡをとったグラフ上に、各スフェアについてプロットを形成することにより、スフェアの体積と位相差量総和θＡの相関特性を取得した。図８は、スフェアの体積と位相差量総和θＡとの相関特性を示すグラフである。図８に示すように、位相差量総和θＡとスフェアの体積とは、比例関係であることが確認された。図８には、スフェアの体積と位相差量総和θＡとの間の相関性の基準を示す基準相関トレンドラインＬＳ及びプロットが示されている。図８に示す各プロットから、導出した回帰線を、基準相関トレンドラインＬＳとして適用した。
図８には、基準相関トレンドラインＬＳ上に存在するプロットａ１及びａ２に対応するスフェアの位相差画像、及び基準相関トレンドラインＬＳから乖離した位置に存在するプロットａ３、ａ４、ａ５に対応するスフェアの位相差画像が示されている。基準相関トレンドラインＬＳ上に存在するプロットａ１及びａ２に対応するスフェアについては、スフェアの全体に亘って輝度が均一な位相差画像が得られた。このことは、スフェアを構成する複数の細胞が均質であること、及びスフェア内における細胞の密度が均一であること等を示している。一方、基準相関トレンドラインＬＳから乖離した位置に存在するプロットａ３及びａ４に対応するスフェアについては、中心部の輝度が他の部位と比較して低下している位相差画像が得られた。このことは、スフェアを構成する複数の細胞が不均質であること、及びスフェア内における細胞の密度が不均一であること等を示している。また、基準相関トレンドラインＬＳから乖離した位置に存在するプロットａ５に対応するスフェアについては、スフェアの輪郭線の凹凸が顕著な位相差画像が得られた。このことは、スフェアを構成する細胞に異常が生じていること等を示している。
以上の結果から、位相差量総和θＡと、スフェアの形状指標値の一例である体積との相関性を用いてスフェアの状態を判定することが可能であるといえる。また、位相差量総和θＡとスフェアの体積との相関性を示す基準相関トレンドラインＬＳと、判定対象のスフェアについての位相差量総和θＡと形状指標値との相関性とを比較し、判定対象のスフェアについての相関性の、基準相関トレンドラインＬＳからの乖離の程度に応じて当該スフェアの状態を判定することが可能であるといえる。具体的には、一方の軸にスフェアの体積をとり、他方の軸に位相差量総和θＡをとったグラフ上にプロットした、判定対象のスフェアについてのプロットの、基準相関トレンドラインＬＳからの乖離の程度に応じて、各スフェアの状態を判定することが可能であるといえる。従って、例えば、位相差量総和θＡの基準相関トレンドラインＬＳからのマイナス幅が閾値以上であるスフェアについては、当該スフェアに含まれる複数の細胞の密度、均質性、及びスフェアの外形形状のうちの少なくとも１つに異常があると判定することができる。
なお、本実施例では、スフェアの形状指標値としてスフェアの体積を用いたが、これに代えて、スフェアの断面積、粒径または周の長さを用いることも可能である。これらいずれの形状指標値を用いる場合でも、位相差量総和θＡとの相関性を用いてスフェアの状態を判定することが可能である。
［実施例２］
本実施例では、複数のスフェアを含む培養ロットを判定対象とした。判定対象の培養ロットは、下記の表１に示す良ロット及び不良ロットである。良ロットに属するスフェア及び不良ロットに属するスフェアは、それぞれ、三次元培養法により培養したｉＰＳ細胞のスフェアである。培養期間の初日から５日目までの細胞の増殖率は、良ロットに属するスフェアおいては２０倍であり、不良ロットに属するスフェアにおいては３．７倍であった。培養５日目における細胞の生存率は、良ロットに属するスフェアおいては９０．２％であり、不良ロットに属するスフェアにおいては６４．１％であった。
良ロット及び不良ロットの各々について、試料１４であるスフェアを図１に示す撮像システム１のサンプルステージにセットし、複数のスフェアのホログラムをＣＭＯＳカメラ１９によって撮影した。取得した各スフェアのホログラムに対してコンピュータによる数値計算を実施することにより、スフェアの中央付近においてスライスした位相差画像を取得した。
得られた各スフェアの位相差画像から、各スフェアの粒径を導出し、スフェアの粒径のヒストグラムを取得した。図９Ａは、良ロット及び不良ロットの各々について導出したスフェアの粒径のヒストグラムである。図９Ａに示すように、良ロットと不良ロットとの間でヒストグラムに有意な差はみられなかった。このことは、スフェアの状態は、当該スフェアの粒径に反映されにくく、スフェアの粒径のみに基づいてスフェアの状態を判定することは困難であることを示している。
次に、良ロット及び不良ロットの各々について得られた各スフェアの位相差画像から、各スフェアの形状指標値として体積を導出し、更に、各スフェアの位相差画像について、（３）式によって示される位相差量総和θＡを導出した。横軸に、スフェアの体積をとり、縦軸に位相差量総和θＡをとったグラフ上に、各培養ロットの各スフェアについてプロットを形成することにより、スフェアの体積と位相差量総和θＡの相関特性を取得した。図９Ｂは、良ロット及び不良ロットの各々について取得した、スフェアの体積と位相差量総和θＡの相関特性を示すグラフである。図９Ｂには、基準相関トレンドラインＬＳ及びプロットが示されている。図９Ｂに示すように、良ロットにおいては、位相差量総和θＡが基準相関トレンドラインＬＳを下回るスフェアが少数であるのに対し、不良ロットにおいては、位相差量総和θＡが基準相関トレンドラインＬＳを下回るスフェアが、良ロットと比べて、多数存在した。
以上の結果から、位相差量総和θＡと、スフェアの形状指標値の一例である体積との相関性を用いて培養ロットの良否判定を行うことが可能であるといえる。また、位相差量総和θＡとスフェアの体積との相関性を示す基準相関トレンドラインＬＳと、判定対象の培養ロットに属するスフェアについての相関性とを比較し、判定対象の培養ロットについての相関性の、基準相関トレンドラインＬＳからの乖離の程度に応じて当該培養ロットの良否判定が可能であるといえる。具体的には、一方の軸にスフェアの体積をとり、他方の軸に位相差量総和θＡをとったグラフ上にプロットした判定対象の培養ロットに属する各スフェアについてのプロットの、基準相関トレンドラインＬＳからの乖離の状態（基準相関トレンドラインＬＳを下回っているか否か）及び乖離の程度の少なくとも一方に基づいて当該培養ロットの状態を判定することが可能であるといえる。従って、例えば、位相差量総和θＡの基準相関トレンドラインＬＳからのマイナス幅が第１の閾値以上であるスフェアの含有率が、第２の閾値以上である場合、当該培養ロットが不良ロットであると判定することができ、上記含有率が、第２の閾値未満である場合、当該培養ロットが良ロットであると判定することができる。
次に、良ロット及び不良ロットの各々について得られた各スフェアの位相差画像から、（４）式によって示される位相差量密度ＤＰを導出した。図９Ｃは、良ロット及び不良ロットの各々について取得した、位相差量密度ＤＰのヒストグラムである。図９Ｃに示すように、良ロットと不良ロットとの間で、ヒストグラムに有意な差が見られた。具体的には、良ロットにおいては、位相差量密度ＤＰが相対的に高い位置に分布のピークが存在し、不良ロットにおいては、位相差量密度ＤＰが相対的に低い位置に分布のピークが存在した。なお、位相差量密度ＤＰのヒストグラムにおける良ロットと不良ロットとの間の有意差の信頼性を示す指標であるｐ値（有意確率）は、０．０５未満であった。
以上の結果から、スフェアの位相差量密度ＤＰに基づいて、培養ロットの良否判定を行うことが可能であるといえる。従って、例えば、判定対象の培養ロットに属する複数のスフェアについて、位相差量密度ＤＰの平均値を求め、上記平均値が閾値以上である場合に、当該培養ロットを良ロットと判定することができ、上記平均値が閾値未満である場合に、当該培養ロットを不良ロットと判定することができる。
なお、本実施例では、位相差量密度ＤＰを培養ロットの良否判定に用いる場合を例示したが、位相差量密度ＤＰをスフェア単体の良否判定に用いることも可能である。例えば、判定対象のスフェアについて取得した位相差量密度ＤＰが閾値以上である場合に、当該スフェアの状態が良好であると判定することができ、位相差量密度ＤＰが閾値未満である場合に、当該スフェアの状態が異常であると判定することができる。
［実施例３］
本実施例では、互いに異なる複数の培養ロットに属する複数のスフェアを判定対象とした。判定対象のスフェアは、いずれも三次元培養法により培養したｉＰＳ細胞のスフェアである。
複数の培養ロットの各々について、試料１４であるスフェアを図１に示す撮像システム１のサンプルステージにセットし、複数のスフェアのホログラムをＣＭＯＳカメラ１９によって撮影した。取得した各スフェアのホログラムに対してコンピュータによる数値計算を実施することにより、スフェアの中央付近においてスライスした位相差画像を取得した。得られた各スフェアの位相差画像について、（４）式によって示される位相差量密度ＤＰを培養ロット毎に導出した。
次に、複数の培養ロットの各々について、細胞の生存率を取得した。生存率の測定手順は以下のとおりである。遠沈管に収容したスフェアを４００Ｇ（Ｇ：重力加速度）、３分間の遠心分離処理を実施した。遠沈管の底に沈殿したスフェアを回収し、細胞解離酵素であるTrypLE（登録商標） Selectを用いて単一細胞に分解した。分解した細胞に対してTrypan Blueによる死細胞染色処理を施した後、細胞を、市販されている一般的なセルカウンタ（Countess（登録商標）） を用いて計数することにより、染色された死細胞の個数をカウントした。死細胞のカウント数から細胞の生存率を培養ロット毎に導出した。
横軸に細胞の生存率をとり、縦軸に位相差量密度ＤＰをとったグラフ上に、各培養ロットについてプロットを形成することにより、細胞の生存率と位相差量密度ＤＰ（ロット平均値）との相関特性を取得した。図１０は、細胞の生存率と位相差量密度ＤＰ（ロット平均値）との相関特性を示すグラフである。図１０に示すように、細胞の生存率と、位相差量密度ＤＰとは比例関係であることが確認された。
以上の結果より、位相差量密度ＤＰから細胞の生存率を推定できるといえる。位相差量密度ＤＰから細胞の生存率を推定することにより、遠心分離処理及び染色処理等の細胞の破壊を伴う処理を行うことなく（すなわち非破壊で）、細胞の生存率を把握することができる。例えば、細胞の生存率と位相差量密度ＤＰとの相関性を示す相関ラインを予め取得しておき、判定対象の培養ロットに属するスフェアについて取得した位相差量密度ＤＰの平均値と、上記の相関ラインとから当該培養ロットにおける細胞の生存率を推定することが可能である。
なお、本実施例では、位相差量密度ＤＰから培養ロット内における細胞の生存率を推定する場合を例示したが、位相差量密度ＤＰから単一のスフェア内における細胞の生存率を推定することも可能である。例えば、細胞の生存率と位相差量密度ＤＰとの相関性を示す相関ラインを予め取得しておき、判定対象のスフェアについて取得した位相差量密度ＤＰと、上記の相関ラインとから当該スフェア内における細胞の生存率を推定することが可能である。
［実施例４］
細胞の生存率が、位相差画像に与える影響について確認した。二次元培養法（接着培養法）により培養したｉＰＳ細胞をネクローシスに誘導するべく、細胞のコロニーに３ｗｔ％のＨ２Ｏ２（過酸化水素）を添加した後、蛍光色素であるSYTOX（登録商標） Greenによる染色処理を施した。続いて、上記の処理を施したコロニーを図１に示す撮像システム１のサンプルステージにセットし、コロニーのホログラムをＣＭＯＳカメラ１９によって撮影した。取得したホログラムに対してコンピュータによる数値計算を実施することにより、コロニーの位相差画像を取得した。ホログラム撮影と並行して、蛍光顕微鏡画像をＣＭＯＳカメラ３６によって撮影した。
図１１Ａ、図１１Ｂ、及び図１１Ｃは、それぞれ、Ｈ２Ｏ２を添加した後、０分経過した時点、１５分経過した時点、３７経過した時点におけるコロニーの蛍光顕微鏡画像である。図１２Ａ、図１２Ｂ、及び図１２Ｃは、それぞれ、Ｈ２Ｏ２を添加した後、０分経過した時点、１５分経過した時点、３７分経過した時点におけるコロニーの位相差画像である。
図１１Ａ〜図１１Ｃに示すように、Ｈ２Ｏ２を添加した時点からの時間経過に伴って、発光箇所が増加した。このことは、Ｈ２Ｏ２を添加した時点からの時間経過に伴って、死細胞が増加したこと（換言すれば、細胞の生存率が低下したこと）を意味する。また、図１２Ａ〜図１２Ｃに示すように、Ｈ２Ｏ２を添加した時点からの時間経過に伴って、位相差画像における低輝度領域の面積が増加した。このことは、細胞の生存率の低下に伴って、位相差画像における位相差量総和θＡが低下することを意味する。このように、細胞の生存率が、位相差量総和θＡに影響を及ぼすことが確認できた。
［実施例５］
二次元培養法（接着培養法）により培養したｉＰＳ細胞のコロニーにＨ２Ｏ２を添加することにより、ネクローシスに誘導した。また、別のｉＰＳ細胞のコロニーについて、PromoKine社のApoptosis Inducer Setを用いてアポトーシスに誘導した。なお、アポトーシスは、多細胞生物の細胞において増殖制御機構として管理・調節された能動的な細胞死である。一方、ネクローシスは、栄養不足、毒物、外傷などの外的環境要因により起こる受動的な細胞死であり、細胞死に至る過程がアポトーシスと異なる。
上記の処理を施したコロニー及び健全なコロニーを図１に示す撮像システム１のサンプルステージにセットし、各コロニーのホログラムをＣＭＯＳカメラ１９によって撮影した。取得したホログラムに対してコンピュータによる数値計算を実施することにより、各コロニーの位相差画像を取得した。各コロニーの位相差画像について、（３）式によって示される位相差量総和θＡを導出した。
図１３は、ネクローシスに誘導したコロニー（死細胞）、アポトーシスに誘導したコロニー（死細胞）及び健全なコロニー（生細胞）のそれぞれの位相差量総和θＡの相対値を示すグラフである。なお、健全なコロニー（生細胞）の位相差量総和θＡを１とした。図１３に示すように、健全なコロニー（生細胞）と、ネクローシスに誘導したコロニー（死細胞）及びアポトーシスに誘導したコロニー（死細胞）との間で、位相差量総和θＡに有意な差が見られた。有意差の信頼性を示す指標であるｐ値（有意確率）は、０．０５未満であった。
以上の結果から細胞死に至る過程が、ネクローシスであるかアポトーシスであるかにかかわらず、位相差量総和θＡを、細胞の生死の判別に利用できるといえる。例えば、判定対象のスフェアについて取得した位相差量総和θＡが、閾値以上である場合、当該スフェアは、生細胞をより多く含むと判定することができ、判定対象のスフェアについて取得した位相差量総和θＡが、閾値未満である場合、当該スフェアは、死細胞をより多く含むと判定することができる。なお、位相差量総和θＡに代えて位相差量密度ＤＰを、細胞の生死の判別に利用することも可能である。
［実施例６］
未分化状態を維持したｉＰＳ細胞のスフェア及び未分化状態を逸脱したｉＰＳ細胞のスフェアを、図１に示す撮像システム１のサンプルステージにセットし、各スフェアのホログラムをＣＭＯＳカメラ１９によって撮影した。取得した各スフェアのホログラムに対してコンピュータによる数値計算を実施することにより、スフェアの中央付近においてスライスした位相差画像を取得した。各スフェアの位相差画像について、（３）式によって示される位相差量総和θＡを導出した。なお、未分化状態を逸脱したスフェアについては、未分化状態の維持に必要な成分を添加しない基礎培地を用いることにより未分化逸脱状態に誘導した。
図１４は、未分化状態を維持したスフェア及び未分化状態を逸脱したスフェアのそれぞれの位相差量総和θＡの平均値を示すグラフである。図１４に示すように、未分化状態を維持したスフェアと、未分化状態を逸脱したスフェアとの間で、位相差量総和θＡに有意な差が見られた。
以上の結果から、位相差量総和θＡを、幹細胞における未分化状態を維持しているか、未分化状態を逸脱しているかの判定に利用できるといえる。例えば、判定対象のスフェアについて取得した位相差量総和θＡが、閾値以上である場合に当該スフェアが未分化状態を維持していると判定することができ、判定対象のスフェアについて取得した位相差量総和θＡが閾値未満である場合に当該スフェアが未分化状態を逸脱したと判定することができる。なお、位相差量総和θＡに代えて位相差量密度ＤＰを、未分化状態を維持しているか、未分化状態を逸脱しているかの判定に利用することも可能である。
以上の説明から明らかなように、開示の技術の実施形態に係る判定方法によれば、スフェアのホログラムから生成された位相差画像と当該スフェアの形状に応じた形状指標値とに基づいて、スフェアの状態が判定されるので、スフェアの状態を非破壊且つ簡便に判定することが可能である。
すなわち、開示の技術の実施形態に係る判定方法によれば、従来の光学顕微鏡では把握し得ないスフェアの内部における細胞の密度及び生存率などについても判定することができる。また、本判定方法において実施されるホログラムの撮影、画像解析及び統計処理などについては、コンピュータによる自動化が容易であり、従来の手法と比較して、人手及び処理時間を削減することができる。また、作業者による手技のばらつきが影響することもなく、安定した判定結果を得ることができる。また、スフェアを単一細胞に分解したり、蛍光色素を添加したりするといった細胞の破壊を伴う処理を行うことなく、スフェアの状態を非破壊で判定することが可能である。また、ホログラム撮影は、三次元空間中にランダムに存在するスフェアについても容易に行うことが可能であり、三次元培養法による培養中のスフェアの状態を、その場で判定することも可能である。
また、スフェアの状態の判定において、球形状のスフェアの中央付近においてスライスした位相差画像を用いることにより、より正確な判定結果を得ることができる。球形状のスフェアの中央付近においてスライスした１枚の位相差画像を用いることにより、スライス位置が互いに異なる複数の位相差画像を用いる場合と比較して、判定に要する処理を簡略化することができる。
開示の技術に係る判定方法は、複数の細胞の凝集体であるスフェアを撮像したホログラムからスフェアの位相差画像を生成し、位相差画像とスフェアの形状に応じた形状指標値とに基づいて、スフェアの状態を判定することを含む。
開示の技術に係る判定方法によれば、スフェアの状態を非破壊且つ簡便に判定することが可能となる。
開示の技術に係る判定方法において、位相差画像と形状指標値とに基づいて、スフェアに含まれる複数の細胞の生存率、密度、均質性、未分化状態逸脱、及びスフェアの外形形状の少なくとも１つに関する判定を行ってもよい。
開示の技術に係る判定方法において、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を導出し、位相差量総和と形状指標値との相関性を用いてスフェアの状態を判定してもよい。
スフェアの状態の判定に位相差量総和と形状指標値との相関性を用いることにより、スフェアの状態について適確な判定を行うことが可能となる。
開示の技術に係る判定方法において、位相差量総和と形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象のスフェアについての位相差量総和と形状指標値との相関性とを比較し、判定対象のスフェアについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の程度に応じて判定対象のスフェアの状態を判定してもよい。例えば、判定対象のスフェアについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の幅が閾値を超える場合、判定対象のスフェアに含まれる細胞の生存率、密度、均質性及びスフェアの外形形状の少なくとも１つについて異常ありと判定してもよい。
この態様によれば、スフェアの状態判定の適確性を更に高めることが可能となる。
開示の技術に係る判定方法において、位相差量総和と形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象の培養ロットに属する複数のスフェアの各々についての位相差量総和と形状指標値との相関性とを比較し、判定対象の培養ロットについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて判定対象の培養ロットの良否判定を行ってもよい。
この態様によれば、複数のスフェアを含む培養ロットについての良否判定を適確に行うことが可能となる。
開示の技術に係る判定方法において、形状指標値は、スフェアの体積、断面積、粒径及び周の長さのいずれかであってもよい。
開示の技術に係る判定方法において、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を、スフェアの体積で除算して位相差量密度を導出し、位相差量密度に基づいて、スフェアの状態を判定してもよい。
スフェアの状態を判定に位相差量密度を用いることにより、スフェアの状態について適確な判定を行うことが可能となる。
開示の技術に係る判定方法において、スフェアの状態の判定に用いる位相差画像は、ホログラムから生成可能な複数の位相差画像のうち、位相差画像を構成する複数の画素間の位相差量のばらつきが最大となる位相差画像であることが好ましい。
スフェアの状態の判定に用いる位相差画像として、上記の画像を用いることで、当該位相差画像から、スフェアの実態に合致したより正確な情報を得ることができ、スフェアの状態判定の適確性を更に高めることが可能となる。
１ 撮像システム
１０ ホログラム光学系
１１ レーザ光源
１２ ビームスプリッタ
１３ コリメートレンズ
１４ 試料
１５ 対物レンズ
１７ 結像レンズ
１８ ビームスプリッタ
１９、３６ ＣＭＯＳカメラ
２０ 光ファイバ
２１ コリメートレンズ
３０ 蛍光顕微鏡光学系
３１ 励起光源
３２ 励起用フィルタ
３３ ＯＮ／ＯＦＦ切り替えミラー
３４ ダイクロイックミラー
３５ 分光フィルタ
５００ コンピュータ
５０２ 主記憶装置
５０３ 補助記憶装置
５０４ 通信インターフェース
５０５ 表示部
５０６ オートフォーカスプログラム
５０７ バス
ＤＰ 位相差量密度
ＩＰ 位相差画像
ＬＳ 基準相関トレンドライン
θ 位相差量
θＢ バックグランドの位相
θＳ スフェアの存在する領域の位相
θＡ 位相差量総和
θｋ １画素の位相差量
Ｖ スフェアの体積
ｖｋ 位相差画像の各画素ｋに対応する部位におけるスフェアの体積
ａ１、ａ２、ａ３、ａ４、ａ５ プロット
ｋ 画素
ｗ カーブの幅

Claims (9)

  1. 複数の細胞の凝集体であるスフェアを撮像したホログラムから前記スフェアの位相差画像を生成し、
    前記位相差画像と前記スフェアの形状に応じた形状指標値とに基づいて、前記スフェアの状態を判定する
    判定方法。
  2. 前記位相差画像と前記形状指標値とに基づいて、前記スフェアに含まれる複数の細胞の生存率、密度、均質性、未分化状態逸脱、及び前記スフェアの外形形状の少なくとも１つに関する判定を行う
    請求項１に記載の判定方法。
  3. 前記位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を導出し、
    前記位相差量総和と前記形状指標値との相関性を用いて前記スフェアの状態を判定する
    請求項１または請求項２に記載の判定方法。
  4. 前記位相差量総和と前記形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象のスフェアについての位相差量総和と形状指標値との相関性とを比較し、
    前記判定対象のスフェアについての位相差量総和と形状指標値との相関性の、前記基準相関トレンドラインからの乖離の程度に応じて前記判定対象のスフェアの状態を判定する
    請求項３に記載の判定方法。
  5. 前記判定対象のスフェアについての位相差量総和と形状指標値との相関性の、前記基準相関トレンドラインからの乖離の幅が閾値を超える場合、前記判定対象のスフェアに含まれる細胞の生存率、密度、均質性及び前記スフェアの外形形状の少なくとも１つについて異常ありと判定する
    請求項４に記載の判定方法。
  6. 前記位相差量総和と前記形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象の培養ロットに属する複数のスフェアの各々についての位相差量総和と形状指標値との相関性とを比較し、
    前記判定対象の培養ロットについての位相差量総和と形状指標値との相関性の、前記基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて前記判定対象の培養ロットの良否判定を行う
    請求項３に記載の判定方法。
  7. 前記形状指標値は、前記スフェアの体積、断面積、粒径及び周の長さのいずれかである
    請求項１から請求項６のいずれか１項に記載の判定方法。
  8. 前記位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を、前記スフェアの体積で除算して位相差量密度を導出し、
    前記位相差量密度に基づいて、前記スフェアの状態を判定する
    請求項１または請求項２に記載の判定方法。
  9. 前記スフェアの状態の判定に用いる位相差画像は、前記ホログラムから生成可能な複数の位相差画像のうち、位相差画像を構成する複数の画素間の位相差量のばらつきが最大となる位相差画像である
    請求項１から請求項８のいずれか１項に記載の判定方法。
JP2020505699A 2018-03-12 2019-02-13 判定方法 Active JP6995975B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2018044594 2018-03-12
JP2018044594 2018-03-12
PCT/JP2019/005112 WO2019176427A1 (ja) 2018-03-12 2019-02-13 判定方法

Publications (2)

Publication Number Publication Date
JPWO2019176427A1 true JPWO2019176427A1 (ja) 2020-12-10
JP6995975B2 JP6995975B2 (ja) 2022-02-04

Family

ID=67908147

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020505699A Active JP6995975B2 (ja) 2018-03-12 2019-02-13 判定方法

Country Status (4)

Country Link
US (1) US11893728B2 (ja)
EP (1) EP3766984A4 (ja)
JP (1) JP6995975B2 (ja)
WO (1) WO2019176427A1 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7046211B2 (ja) * 2018-09-28 2022-04-01 富士フイルム株式会社 判定方法
JP7123155B2 (ja) * 2018-09-28 2022-08-22 富士フイルム株式会社 判定方法
JPWO2021085033A1 (ja) * 2019-10-28 2021-05-06
WO2021085034A1 (ja) * 2019-10-28 2021-05-06 富士フイルム株式会社 多能性幹細胞の選別方法、分化誘導結果の予測方法及び細胞製品の製造方法
CN112765624B (zh) * 2021-01-19 2022-05-27 浙江科技学院 基于相位优化和稀疏约束的可认证纯相位全息图生成方法
US20220254490A1 (en) * 2021-02-05 2022-08-11 Anvitha Addanki Data processing system for estimating disease progression rates
CN117098987A (zh) 2021-04-09 2023-11-21 富士胶片株式会社 全息装置及细胞的评价方法
JPWO2022224722A1 (ja) 2021-04-21 2022-10-27
WO2023233827A1 (ja) * 2022-06-03 2023-12-07 富士フイルム株式会社 位相画像取得方法及び定量データ取得方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080309944A1 (en) * 2006-04-21 2008-12-18 Consiglio Nazionale Delle Ricerche Quantitative phase-contrast digital holography method for the numerical reconstruction of images, and relevant apparatus
JP2009521216A (ja) * 2005-12-22 2009-06-04 フェイズ ホログラフィック イメージング ペーハーイー アーベー 細胞サンプル分析のための方法と装置
JP2015192644A (ja) * 2014-03-26 2015-11-05 株式会社Ｓｃｒｅｅｎホールディングス スフェロイドの評価方法およびスフェロイド評価装置
US20170261930A1 (en) * 2011-07-19 2017-09-14 Ovizio Imaging Systems NV/SA Method and system for detecting and/or classifying cancerous cells in a cell sample

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011160068A1 (en) * 2010-06-17 2011-12-22 Purdue Research Foundation Digital holographic method of measuring cellular activity and measuring apparatus with improved stability
US8531747B2 (en) * 2009-06-24 2013-09-10 Canon Kabushiki Kaisha Hologram, hologram data generation method, and exposure apparatus
CN105009002B (zh) * 2013-03-01 2019-03-05 西铁城时计株式会社 光束分割元件
JP2015146747A (ja) 2014-02-05 2015-08-20 浜松ホトニクス株式会社 細胞判定方法
JP6284832B2 (ja) 2014-06-20 2018-02-28 富士フイルム株式会社 細胞評価装置および方法並びにプログラム
WO2016084420A1 (ja) 2014-11-27 2016-06-02 株式会社島津製作所 ディジタルホログラフィ装置及びディジタルホログラム生成方法
WO2017159367A1 (ja) * 2016-03-18 2017-09-21 株式会社村田製作所 金属製多孔膜、それを用いた分級方法、および分級装置
US10846849B2 (en) * 2016-06-16 2020-11-24 Hitachi High-Tech Corporation Method for analyzing state of cells in spheroid
JP6376196B2 (ja) 2016-09-29 2018-08-22 株式会社ニコン 細胞評価装置、細胞評価方法およびプログラム
JP7053853B2 (ja) * 2018-08-20 2022-04-12 富士フイルム株式会社 細胞培養方法及び細胞培養装置
JP7159334B2 (ja) * 2018-09-28 2022-10-24 富士フイルム株式会社 細胞構造体、細胞構造体の製造方法、細胞培養方法及びマイクロ流路

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009521216A (ja) * 2005-12-22 2009-06-04 フェイズ ホログラフィック イメージング ペーハーイー アーベー 細胞サンプル分析のための方法と装置
US20080309944A1 (en) * 2006-04-21 2008-12-18 Consiglio Nazionale Delle Ricerche Quantitative phase-contrast digital holography method for the numerical reconstruction of images, and relevant apparatus
US20170261930A1 (en) * 2011-07-19 2017-09-14 Ovizio Imaging Systems NV/SA Method and system for detecting and/or classifying cancerous cells in a cell sample
JP2015192644A (ja) * 2014-03-26 2015-11-05 株式会社Ｓｃｒｅｅｎホールディングス スフェロイドの評価方法およびスフェロイド評価装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
APPLIED OPTICS (2011) VOL.50, NO.20, PP.3589-3597, JPN6019015489, ISSN: 0004495718 *
OPTICS LETTERS (2005) VOL.30, NO.5, PP.468-470, JPN6019015487, ISSN: 0004495717 *

Also Published As

Publication number Publication date
WO2019176427A1 (ja) 2019-09-19
JP6995975B2 (ja) 2022-02-04
EP3766984A1 (en) 2021-01-20
US11893728B2 (en) 2024-02-06
EP3766984A4 (en) 2021-05-12
US20200342599A1 (en) 2020-10-29

Similar Documents

Publication Publication Date Title
JP6995975B2 (ja) 判定方法
Bian et al. Autofocusing technologies for whole slide imaging and automated microscopy
Sikkel et al. High speed sCMOS‐based oblique plane microscopy applied to the study of calcium dynamics in cardiac myocytes
US20210215923A1 (en) Microscope system
Seo et al. Multi-color LUCAS: Lensfree on-chip cytometry using tunable monochromatic illumination and digital noise reduction
CN108254295B (zh) 一种定位与表征球形微粒的方法及其装置
JP2018502283A (ja) 生物学的粒子の位置の決定を含む分析方法
WO2021188839A1 (en) Single-shot autofocusing of microscopy images using deep learning
WO2018158901A1 (ja) 細胞解析方法及び細胞解析装置
JP2020514704A (ja) レンズレス撮像によるサンプル内の粒子の計数方法
EP3532822B1 (en) Trans-illumination imaging with use of interference fringes to enhance contrast and find focus
US11480920B2 (en) Image processing apparatus, evaluation system, image processing program, and image processing method
JP6867731B2 (ja) 微粒子観察装置及び微粒子観察方法
WO2018186120A1 (ja) 細胞塊の評価方法及び細胞塊の状態解析装置
Go et al. Learning-based automatic sensing and size classification of microparticles using smartphone holographic microscopy
US11869181B2 (en) Determination method
CN108593528B (zh) 基于激光干涉的非球形粗糙粒子形状和尺寸测量方法
US20240102912A1 (en) Plaque counting assay method
US20150285785A1 (en) Monitoring and/or characterising biological or chemical material
WO2021085033A1 (ja) 多能性幹細胞の選別方法、分化誘導結果の予測方法及び細胞製品の製造方法
WO2021085034A1 (ja) 多能性幹細胞の選別方法、分化誘導結果の予測方法及び細胞製品の製造方法
JP7046211B2 (ja) 判定方法
US20230341748A1 (en) Method and apparatus for displaying cultured cells
WO2023059764A1 (en) Method and apparatus for searching and analyzing cell images
Vandewiele et al. Low Coherence Digital Holography Microscopy Enhanced by Optical Trapping

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200511

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210427

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210628

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20211116

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20211215

R150 Certificate of patent or registration of utility model

Ref document number: 6995975

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150