JP2023021090A - 細胞自家蛍光の画像を用いて細胞を分類する方法及び解析装置 - Google Patents
細胞自家蛍光の画像を用いて細胞を分類する方法及び解析装置 Download PDFInfo
- Publication number
- JP2023021090A JP2023021090A JP2022121735A JP2022121735A JP2023021090A JP 2023021090 A JP2023021090 A JP 2023021090A JP 2022121735 A JP2022121735 A JP 2022121735A JP 2022121735 A JP2022121735 A JP 2022121735A JP 2023021090 A JP2023021090 A JP 2023021090A
- Authority
- JP
- Japan
- Prior art keywords
- cells
- autofluorescence
- image
- information
- cell
- 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.)
- Pending
Links
Images
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
- C12M1/00—Apparatus for enzymology or microbiology
- C12M1/34—Measuring or testing with condition measuring or sensing means, e.g. colony counters
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- 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
-
- 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
Abstract
【課題】細胞の自家蛍光を用いて、細胞の分類をすることを課題とする。【解決手段】複数の細胞に関する自家蛍光画像を用いて、前記複数の細胞に関する自家蛍光情報を取得する工程と、前記領域情報を取得する工程と、前記複数の細胞に関する自家蛍光情報と前記領域情報から、前記個々の細胞の自家蛍光情報を抽出する工程と、前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類する工程を含む方法を提供する。【選択図】図8
Description
本発明は、細胞自家蛍光の画像を用いて細胞を分類する方法、及び細胞の解析装置に関する。
細胞内因性成分には、蛍光を発するものがあることが知られる。細胞内の還元状態や代謝状態に関連する補酵素のうちの複数は蛍光を発する。長谷川らは特に生体組織の酸化還元状態により蛍光輝度が上下する補酵素NADHとFADに着目した。そして、臓器の撮影画像において、NADHとFADに対応するそれぞれの蛍光波長の輝度比を算出し、生体組織の酸化還元状態の判断を行うシステムの開発を行った(非特許文献1)。
高解像度顕微鏡とフェムト秒レーザーを組み合わせた蛍光寿命イメージング顕微鏡が開発されている。蛍光寿命イメージング顕微鏡を用いて、個々の細胞内のNADHやFADの絶対値の定量から、常細胞ががん化していく過程における自家蛍光波長成分の変化が測定されている(非特許文献2)。
また、バクテリア等の自家蛍光と形状や動きの情報を画像より細胞単位で取得し、そのデータを教師として機械学習により特徴づけを行うことで、未知サンプルの情報を解析する手段の開発が行われている(特許文献1)。さらに、染色サンプルの蛍光画像について、高解像度の画像から細胞内の蛍光輝度を空間分布測定することで細胞の分類を行う手法が開発されている(特許文献2)。
一方、一般的に同種の細胞であっても自家蛍光の発色要因の補酵素の絶対量は決まっておらず、培養環境、方法や細胞が由来する個体や臓器の違い等により自家蛍光の輝度やスペクトル情報は異なる。そのため個々の細胞の自家蛍光の特定の波長の輝度の絶対量や、それぞれの波長の輝度の演算等では細胞一つ一つの表現型や状態を判定することはできない。したがって、細胞種や細胞の状態を特定したり、推定するためには、細胞の表現型と自家蛍光輝度と相関関係が分かる対照実験や検量データを作成するか、あるいは、特許文献2のような教師データの蓄積等が必要であった。
顕微イメージングでは視野内に観察される細胞の差異は確認できるが、同時性のある大域的な情報が得られない。したがって、例えば、数センチ四方といった広域な組織中や培養シャーレ内に含まれる多数の細胞種を表現型により分類することは難しい。特許文献1の方法では大域的な情報を取得するが、組織全体としての変化や差異を確認するのみであり、個々の細胞情報の取得はできず、表現型による細胞種の分類はできない。特許文献1の方法では、自家蛍光の各種の波長における輝度変化の要因が、例えば細胞数の変化によるものなのか個々の細胞の輝度の変化なのかを確認することができない。
K. Hasegawa, Y. Ogasawara, 可視化情報学会論文集 2010, 30, 11, 73-79
M. C. Skala, K. M. Riching, A. Gendron-Fitzpatrick, J. Eickhoff, K. W. Eliceiri, J. G. White, and N. Ramanujam, PNAS, 2007, 104, 49, 19494-19499
本発明は、複数の細胞の自家蛍光情報を用いて、細胞の分類を行うことを課題とする。
本発明者らは、複数の細胞に関し、それらの個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、複数の細胞を2以上の群に分類することが可能であることを見出した。
すなわち、本発明に係る方法は、複数の細胞に関する自家蛍光画像を用いて、前記複数の細胞に関する自家蛍光情報を取得する工程と、前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得する工程と、前記複数の細胞に関する自家蛍光情報と前記領域情報から、前記個々の細胞の自家蛍光情報を抽出する工程と、前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類する工程、を含む。
本発明によれば、複数の細胞を、個々の細胞の自家蛍光情報を用いて、細胞の分類をすることができる。本発明によれば、染色の必要がなく、非侵襲に細胞を分類することができる。
以下に、より具体的な構成例について記載するが、本発明は下記方法に限定されるものではない。
(第一の実施形態)
本発明は、第一の実施形態として、
複数の細胞に関する自家蛍光画像を用いて、前記複数の細胞に関する自家蛍光情報を取得する工程と、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得する工程と、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、前記個々の細胞の自家蛍光情報を抽出する工程と、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類する工程、を含む方法を提供する。
このように本実施形態によれば、クラスタリングを行うことで複数の細胞の分類を行うので、教師データを必要としない。また、自家蛍光画像を用いるため、染色の必要がなく、非侵襲に細胞を分類することができる。以下、詳細に説明する。
本発明は、第一の実施形態として、
複数の細胞に関する自家蛍光画像を用いて、前記複数の細胞に関する自家蛍光情報を取得する工程と、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得する工程と、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、前記個々の細胞の自家蛍光情報を抽出する工程と、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類する工程、を含む方法を提供する。
このように本実施形態によれば、クラスタリングを行うことで複数の細胞の分類を行うので、教師データを必要としない。また、自家蛍光画像を用いるため、染色の必要がなく、非侵襲に細胞を分類することができる。以下、詳細に説明する。
(細胞)
本実施形態の方法にはあらゆる細胞が用いられるが、好ましくは哺乳細胞が用いられる。細胞は培養細胞、組織切片の細胞、組織、生体試料から抽出された細胞等が対象となり、正常細胞、がん細胞、あるいは、がん細胞と正常細胞、がん細胞同士、または正常細胞同士の混合培養、あるいは、混合組織等が対象となる。また、がん細胞の上皮間葉転換、未分化細胞からの分化誘導、細胞の老化等の同種細胞の形質の有無が予想される細胞が対象となる。
本実施形態の方法にはあらゆる細胞が用いられるが、好ましくは哺乳細胞が用いられる。細胞は培養細胞、組織切片の細胞、組織、生体試料から抽出された細胞等が対象となり、正常細胞、がん細胞、あるいは、がん細胞と正常細胞、がん細胞同士、または正常細胞同士の混合培養、あるいは、混合組織等が対象となる。また、がん細胞の上皮間葉転換、未分化細胞からの分化誘導、細胞の老化等の同種細胞の形質の有無が予想される細胞が対象となる。
本実施形態の方法においては、細胞から、自家蛍光画像、あるいは、自家蛍光画像と領域抽出用画像が取得される。これらの画像は、細胞を、顕微鏡などの拡大レンズを含む光学系を介してカメラ等の画像取得装置によって、撮像して取得される。
細胞は、その画像取得装置によって撮像が可能な状態であれば、どのような状態であってもよい。たとえば、画像取得装置が一般的な顕微鏡である場合、細胞は、シャーレ、プレート、フラスコ等に、懸濁された、付着した、培養された状態であることができ、あるいは、切片中に含まれてもよい。撮像が可能であれば、細胞は、血液や、組織といった、生体中に存在するものでもよい。
本実施形態の方法は、生きたままの細胞に適用できるという特徴を有するが、対象は生細胞に限定されることはなく、例えば、細胞は固定化されたものや、固定化された切片に含まれるものであっても構わない。
(自家蛍光画像)
自家蛍光とは、細胞や組織に内在する物質が、励起波長の照射を受けて、発する蛍光をいう。特に自家蛍光を発する内因性成分としてはNAD(P)H、フラビン類(FAD等)、コラーゲン、フィブロネクチン、トリプトファン、葉酸等、細胞内で産生される自家蛍光分子があげられる。
自家蛍光とは、細胞や組織に内在する物質が、励起波長の照射を受けて、発する蛍光をいう。特に自家蛍光を発する内因性成分としてはNAD(P)H、フラビン類(FAD等)、コラーゲン、フィブロネクチン、トリプトファン、葉酸等、細胞内で産生される自家蛍光分子があげられる。
本実施形態の方法においては、自家蛍光を発している内因性成分が何であるかを特定される必要はなく、自家蛍光は、上述に例を記した蛍光性内因性成分以外の内因性成分が発する蛍光も含む。
自家蛍光画像は、複数の細胞を含む1の試料に対し、1またはそれ以上取得することができる。自家蛍光画像は、特定波長の励起光を照射し、発せられる蛍光から得られる。励起光は1の波長の光であってもよいし、複数の波長の光からなる光でもよい。また、蛍光は、1の波長の光であってもよいし、複数の波長の光からなる光を蛍光としてもよい。
励起光の照射の例として、350nm以上の波長を発する水銀ランプ、キセノンランプ等の白色光源を用いて、測定しようとする蛍光波長領域以上をカットするショートパスフィルターを通して照射することを挙げられる。また、別の例として、中心波長を有する特定の領域波長を出射するLEDを光源とし、測定しようとする蛍光波長領域以上をカットするショートパスフィルターを通して照射することを挙げられるが、これらに限定されるものではない。
自家蛍光画像の撮像の際の波長の例として、365nm励起/430nmロングパスフィルター観察、365nm励起/530nmロングパスフィルター観察、395nm励起/430nmロングパスフィルター観察、365nm励起/530nmロングパスフィルター観察あるいは450nm励起/530nmロングパスフィルター観察を挙げられるが、これらに限定されるものではない。
(自家蛍光情報)
複数の細胞に関する自家蛍光情報は、複数の細胞に関する自家蛍光画像から取得することができる。自家蛍光情報は、自家蛍光の輝度を含み、さらには、輝度の位置情報、面積情報、その他の情報を含むことができる。輝度は、蛍光強度を示す値であれば特に限定されず、どのような指標を用いてもよい。なお、輝度は強度とも換言できる。
複数の細胞に関する自家蛍光情報は、複数の細胞に関する自家蛍光画像から取得することができる。自家蛍光情報は、自家蛍光の輝度を含み、さらには、輝度の位置情報、面積情報、その他の情報を含むことができる。輝度は、蛍光強度を示す値であれば特に限定されず、どのような指標を用いてもよい。なお、輝度は強度とも換言できる。
輝度は、色空間座標上で表し、その座標で得られる値としてもよい。色空間の一例としてRGBの色空間を用いることができる。輝度はRGBの色空間におけるR、G、及びBの少なくともいずれかの指標として算出される値であってもよい。また、輝度は、Lab色空間またはHSV色空間、あるいはHLS色空間における指標として算出される値であってもよい。これら以外の色彩空間を用いて輝度情報を記述してもよい。
本明細書において、色空間におけるRは600nm以上800nm以下、Gは500nm以上600nm未満、Bは350nm以上500nm未満の範囲のいずれかの波長とすることができる。
本明細書において、色空間におけるRは600nm以上800nm以下、Gは500nm以上600nm未満、Bは350nm以上500nm未満の範囲のいずれかの波長とすることができる。
複数の細胞に関する自家蛍光情報と個々の細胞の領域を示す領域情報から、個々の細胞の自家蛍光情報が抽出される。
個々の細胞の自家蛍光情報として、細胞1個につき、1またはそれ以上の自家蛍光情報が決定されることが好ましく、代表的には、細胞1個につき、1の自家蛍光情報が決定される。
個々の細胞の自家蛍光情報は、個々の細胞の自家蛍光の輝度とすることができ、さらには、該当細胞の輝度の代表値をその細胞の自家蛍光情報とすることができる。代表値として、最小値、最大値、平均値、中央値等を用いることができる。
例えば、色空間におけるRを指標として自家蛍光情報を取得するというのは、自家蛍光画像における個々の細胞に関し、600nm以上800nm以下の範囲のいずれかの波長の輝度の値を取得することを意味する。複数の細胞を2以上の群に分類するためには、R、G、Bのうちの少なくとも1つの指標に関する輝度の値をプロットしてクラスタリングを行う必要がある。
個々の細胞の自家蛍光情報として、細胞1個につき、1またはそれ以上の自家蛍光情報が決定されることが好ましく、代表的には、細胞1個につき、1の自家蛍光情報が決定される。
個々の細胞の自家蛍光情報は、個々の細胞の自家蛍光の輝度とすることができ、さらには、該当細胞の輝度の代表値をその細胞の自家蛍光情報とすることができる。代表値として、最小値、最大値、平均値、中央値等を用いることができる。
例えば、色空間におけるRを指標として自家蛍光情報を取得するというのは、自家蛍光画像における個々の細胞に関し、600nm以上800nm以下の範囲のいずれかの波長の輝度の値を取得することを意味する。複数の細胞を2以上の群に分類するためには、R、G、Bのうちの少なくとも1つの指標に関する輝度の値をプロットしてクラスタリングを行う必要がある。
(領域情報)
領域情報は、自家蛍光画像における個々の細胞の領域を示す情報である。複数の細胞に関する自家蛍光情報を、それぞれの細胞ごとに対応させるために、細胞の領域情報が必要である。領域情報は、少なくとも、細胞の位置、大きさ、面積、形状、厚み、細胞の分布のいずれかの情報を含む。
領域情報が細胞の位置の情報である場合、自家蛍光画像を各細胞に対応できれば、位置はどのように特定されてもよく、例えば、顕微鏡の視野上の位置、培養皿や組織等の、試料上の位置等で特定されてもよい。
領域情報は、自家蛍光画像における個々の細胞の領域を示す情報である。複数の細胞に関する自家蛍光情報を、それぞれの細胞ごとに対応させるために、細胞の領域情報が必要である。領域情報は、少なくとも、細胞の位置、大きさ、面積、形状、厚み、細胞の分布のいずれかの情報を含む。
領域情報が細胞の位置の情報である場合、自家蛍光画像を各細胞に対応できれば、位置はどのように特定されてもよく、例えば、顕微鏡の視野上の位置、培養皿や組織等の、試料上の位置等で特定されてもよい。
領域情報は、自家蛍光画像を用いて取得することができる。すなわち、自家蛍光画像から、個々の細胞の位置、大きさ、面積、形状、厚み、細胞の分布等を知ることができ、これを用いて、個々の細胞の領域を知ることができる。
あるいは、領域情報は、後述する領域抽出用画像から取得してもよい。
あるいは、領域情報は、後述する領域抽出用画像から取得してもよい。
領域情報は、例えば、領域抽出用画像あるいは自家蛍光画像から作成された、自家蛍光画像あるいは自家蛍光情報と重ね合わせて用いられるマスク画像の情報であってもよい。
マスク画像は、自家蛍光情報の、個々の細胞に相当する情報を抽出し、それ以外をマスクするように機能する。
マスク画像は、自家蛍光情報の、個々の細胞に相当する情報を抽出し、それ以外をマスクするように機能する。
具体的な一例を次に説明する。領域抽出用画像を二値化し、細胞位置情報を画像マスクとすることができる。これを、あらかじめ、RGB各色情報として表された自家蛍光画像と重ね合わせて、個々の細胞のRGBの値を得ることができる。
上記において、画像マスクは、領域情報の一例であり、RGB各色情報として表された自家蛍光画像は、複数の細胞に関する自家蛍光情報の一例であり、個々の細胞のRGBの値は、個々の細胞の自家蛍光情報の一例である。
(領域抽出用画像)
領域抽出用画像は、自家蛍光画像における個々の細胞の領域情報を得るために自家蛍光画像とは別に取得される画像であり、個々の細胞そのもの、細胞膜、核、その他細胞内小器官を識別できるものでありうる。領域抽出用画像は明視野画像であってもよいし、例えば核を例とする特定の細胞内小器官を染色した上で得られた、明視野画像、あるいは蛍光画像として取得してもよい。ただし、領域抽出用画像が蛍光画像である場合、その励起波長または蛍光波長が、自家蛍光のそれらの波長と重複しないよう選択することが好ましい。領域抽出用画像は1の試料に対し、1またはそれ以上取得することができる。
領域抽出用画像は明視野で拡大結像及び撮影することにより得られる明視野画像であってもよい。あるいは染色された試料もしくは非染色の試料の蛍光撮影画像であってもよい。
領域抽出用画像は、自家蛍光画像における個々の細胞の領域情報を得るために自家蛍光画像とは別に取得される画像であり、個々の細胞そのもの、細胞膜、核、その他細胞内小器官を識別できるものでありうる。領域抽出用画像は明視野画像であってもよいし、例えば核を例とする特定の細胞内小器官を染色した上で得られた、明視野画像、あるいは蛍光画像として取得してもよい。ただし、領域抽出用画像が蛍光画像である場合、その励起波長または蛍光波長が、自家蛍光のそれらの波長と重複しないよう選択することが好ましい。領域抽出用画像は1の試料に対し、1またはそれ以上取得することができる。
領域抽出用画像は明視野で拡大結像及び撮影することにより得られる明視野画像であってもよい。あるいは染色された試料もしくは非染色の試料の蛍光撮影画像であってもよい。
(第3の情報)
自家蛍光情報または領域情報とは別に、第3の情報を取得する工程を含んでもよい。第3の情報は、明視野画像から取得してもよいし、蛍光画像から取得してもよいし、画像以外から取得される情報であってもよい。第3の情報の例としては、細胞中の特定のマーカーの有無の情報や、その他の情報を挙げることができ、例えば、所望のマーカーを特異的に染色した免疫染色画像を挙げることができる。
自家蛍光情報または領域情報とは別に、第3の情報を取得する工程を含んでもよい。第3の情報は、明視野画像から取得してもよいし、蛍光画像から取得してもよいし、画像以外から取得される情報であってもよい。第3の情報の例としては、細胞中の特定のマーカーの有無の情報や、その他の情報を挙げることができ、例えば、所望のマーカーを特異的に染色した免疫染色画像を挙げることができる。
(明視野画像)
領域抽出用画像として、あるいは、第3の情報の取得のために、明視野画像を用いることができる。明視野画像取得方法は次の例が挙げられる。すなわち、可視光波長の光源または、複数の光源波長の混合である白色光源を対物レンズの光軸に平行か、あるいは平行でなくとも光軸に対して垂直ではない角度で照射する。そして、対象試料からの反射光または複屈折による回折光、及びそれらの干渉を検出することができる。
領域抽出用画像として、あるいは、第3の情報の取得のために、明視野画像を用いることができる。明視野画像取得方法は次の例が挙げられる。すなわち、可視光波長の光源または、複数の光源波長の混合である白色光源を対物レンズの光軸に平行か、あるいは平行でなくとも光軸に対して垂直ではない角度で照射する。そして、対象試料からの反射光または複屈折による回折光、及びそれらの干渉を検出することができる。
(蛍光画像)
自家蛍光画像、領域抽出用画像、あるいは、第3の情報は、蛍光画像を用いることができる。蛍光画像取得方法の例を次に説明する。蛍光励起光源波長として紫外から可視光までの特定の波長を中心波長とした光源を、対物レンズの光軸に平行か、あるいは、少なくとも光軸に対して垂直ではない角度で試料に照射し、試料の蛍光を励起させる。励起により生じた蛍光は対物レンズ等の観察側の対物レンズ等の前方もしくは後方に設置したカットフィルターを介して検出される。励起光及び蛍光観察カットフィルターは励起光源波長の一部が観察側の蛍光カットフィルターを透過しないよう両者の波長カットフィルターが選択される。
自家蛍光画像、領域抽出用画像、あるいは、第3の情報は、蛍光画像を用いることができる。蛍光画像取得方法の例を次に説明する。蛍光励起光源波長として紫外から可視光までの特定の波長を中心波長とした光源を、対物レンズの光軸に平行か、あるいは、少なくとも光軸に対して垂直ではない角度で試料に照射し、試料の蛍光を励起させる。励起により生じた蛍光は対物レンズ等の観察側の対物レンズ等の前方もしくは後方に設置したカットフィルターを介して検出される。励起光及び蛍光観察カットフィルターは励起光源波長の一部が観察側の蛍光カットフィルターを透過しないよう両者の波長カットフィルターが選択される。
(分類する工程)
本実施形態において分類する工程は、個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、複数の細胞を2以上の群に分類する。分類された結果が、ユーザに分かりやすいように可視化することができる。例えば、少なくとも個々の細胞の自家蛍光の輝度を1つの軸とする、1次元、2次元あるいは3次元のグラフとして可視化することができる。グラフは、自家蛍光の輝度のみを軸とする1次元のグラフでもよい。あるいは、2以上の自家蛍光の輝度を軸とする2次元、3次元のグラフでもよい。さらには、自家蛍光の輝度を1つの軸とし、自家蛍光の輝度以外の、領域情報、あるいは、第3の情報をその他の軸とする2次元、3次元のグラフでもよい。グラフを用いて、目視で、複数の細胞を2以上の群にすることができる。なお、グラフの軸は、後述する特徴量の一種といえる。
本実施形態において分類する工程は、個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、複数の細胞を2以上の群に分類する。分類された結果が、ユーザに分かりやすいように可視化することができる。例えば、少なくとも個々の細胞の自家蛍光の輝度を1つの軸とする、1次元、2次元あるいは3次元のグラフとして可視化することができる。グラフは、自家蛍光の輝度のみを軸とする1次元のグラフでもよい。あるいは、2以上の自家蛍光の輝度を軸とする2次元、3次元のグラフでもよい。さらには、自家蛍光の輝度を1つの軸とし、自家蛍光の輝度以外の、領域情報、あるいは、第3の情報をその他の軸とする2次元、3次元のグラフでもよい。グラフを用いて、目視で、複数の細胞を2以上の群にすることができる。なお、グラフの軸は、後述する特徴量の一種といえる。
なお、グラフとは、数値情報を可視化するものすべてを指し、グラフは情報を可視化できればよく、1次元であればヒストグラム、2次元・3次元であれば散布図を例示することができるがこれらに限定されるものではない。さらには、カーネル密度推定等の確率密度関数プロットやヒートマップ等をグラフとして用いてもよい。
あるいは、個々の細胞の自家蛍光情報を用いて複数の細胞を2以上の群に分類する工程では、教師なし機械学習により、複数の細胞を2以上の群に分類することができる。教師なし機械学習の例として、主成分分析、k-means法、GMMを挙げられる。
あるいは、個々の細胞の自家蛍光情報を用いて複数の細胞を2以上の群に分類する工程では、教師なし機械学習により、複数の細胞を2以上の群に分類することができる。教師なし機械学習の例として、主成分分析、k-means法、GMMを挙げられる。
機械学習の際には、特徴量が必要である。特徴量は、データの特徴が数値化されたものである。特徴量の、少なくとも1つとして、個々の細胞の自家蛍光情報が挙げられ、より具体的には、個々の細胞の自家蛍光の輝度を挙げられる。さらに、自家蛍光の輝度以外の、領域情報や、第3の情報を特徴量として用いることができる。自家蛍光情報が複数ある場合、それをさらなる特徴量として用いてもよい。
(表現型への対応)
分類した2以上の群は、群ごとに細胞の表現型に対応させることができる。2以上の群が対応する表現型の組合せの例としてがん細胞と非がん細胞、悪性度がそれぞれ異なるがん細胞群、未分化細胞と分化細胞、分化細胞と脱分化細胞、及び同種がん細胞の上皮系と間葉系細胞等を挙げることができる。
したがって、本実施形態は例えば、組織や培養細胞中のがん細胞の悪性度の判別、各種幹細胞の未分化細胞の維持培養、分化誘導等における分化の違いの判定等に用いることができる。
分類した2以上の群は、群ごとに細胞の表現型に対応させることができる。2以上の群が対応する表現型の組合せの例としてがん細胞と非がん細胞、悪性度がそれぞれ異なるがん細胞群、未分化細胞と分化細胞、分化細胞と脱分化細胞、及び同種がん細胞の上皮系と間葉系細胞等を挙げることができる。
したがって、本実施形態は例えば、組織や培養細胞中のがん細胞の悪性度の判別、各種幹細胞の未分化細胞の維持培養、分化誘導等における分化の違いの判定等に用いることができる。
(定量的な情報の取得)
本実施形態において、分類する工程によって分けられた群に関する定量的な情報を取得する工程を含んでいてもよい。ここで、定量的な情報とは、各群に属する細胞の数、各群における細胞の数の比、色空間におけるRGBのうちの少なくとも2つの輝度の比、各群における輝度の代表値の比、などが挙げられる。ここで、代表値として、最小値、最大値、平均値、中央値等を用いることができる。
本実施形態において、分類する工程によって分けられた群に関する定量的な情報を取得する工程を含んでいてもよい。ここで、定量的な情報とは、各群に属する細胞の数、各群における細胞の数の比、色空間におけるRGBのうちの少なくとも2つの輝度の比、各群における輝度の代表値の比、などが挙げられる。ここで、代表値として、最小値、最大値、平均値、中央値等を用いることができる。
(視野・画素数・解像度)
細胞の分類は、培養皿上の細胞の全て、あるいは組織に含まれる細胞の全て等、多数の細胞に適用されることが想定される。一方、従来の顕微鏡の観察では、一視野内で限られた数の細胞しか同時に観察できなかった。一括した情報を得るためには、高画素あるいは広視野な撮像を用いることが好ましい。高画素または高視野の撮像を用いることで、情報を一括して得られるだけでなく、個々の細胞ごとの自家蛍光情報を非常に多数得られ、より、正確な分類が可能となる。
この観点より、複数の細胞に関する自家蛍光画像は、1視野が細胞培養容器全体を代表するのに統計的に十分な複数の細胞を一括撮像できる、60mm2以上の領域を撮像して取得されることが好ましい。
35mmdishのコンフルエントでの細胞数と言われている細胞数1.2×106個が培養されているとき、1視野中に15000個程度の細胞があると、統計的に1%の許容誤差と99%の信頼度で、シャーレ1枚を1視野で代表することができる。このため、撮影時の視野としては例えば哺乳類の細胞サイズ20μmとし、これを15000個敷き詰めたときのサイズである約60mm2の視野が必要である。ただし、この程度まで許容誤差の低減、信頼性が必要ないと判断する場合、より少ない細胞数で充分なため、視野はこれより低くても構わない。
細胞の分類は、培養皿上の細胞の全て、あるいは組織に含まれる細胞の全て等、多数の細胞に適用されることが想定される。一方、従来の顕微鏡の観察では、一視野内で限られた数の細胞しか同時に観察できなかった。一括した情報を得るためには、高画素あるいは広視野な撮像を用いることが好ましい。高画素または高視野の撮像を用いることで、情報を一括して得られるだけでなく、個々の細胞ごとの自家蛍光情報を非常に多数得られ、より、正確な分類が可能となる。
この観点より、複数の細胞に関する自家蛍光画像は、1視野が細胞培養容器全体を代表するのに統計的に十分な複数の細胞を一括撮像できる、60mm2以上の領域を撮像して取得されることが好ましい。
35mmdishのコンフルエントでの細胞数と言われている細胞数1.2×106個が培養されているとき、1視野中に15000個程度の細胞があると、統計的に1%の許容誤差と99%の信頼度で、シャーレ1枚を1視野で代表することができる。このため、撮影時の視野としては例えば哺乳類の細胞サイズ20μmとし、これを15000個敷き詰めたときのサイズである約60mm2の視野が必要である。ただし、この程度まで許容誤差の低減、信頼性が必要ないと判断する場合、より少ない細胞数で充分なため、視野はこれより低くても構わない。
また、以上の領域における複数の細胞を一括撮像して自家蛍光画像を取得するために、15000個の細胞を個別に認識する必要がある。これを満たすために、複数の細胞に関する自家蛍光画像は、少なくとも30000画素以上の画素数と、かつ1画素当り20μm以下の解像度を有することが好ましい。
(第二の実施形態)
本発明は第二の実施形態として、
画像取得装置と画像処理装置を含む解析装置であって、
前記画像取得装置は、
複数の細胞に関する自家蛍光画像を取得し、
前記画像処理装置は、
前記自家蛍光画像を用いて前記複数の細胞に関する自家蛍光情報を取得し、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得し、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、個々の細胞の自家蛍光情報を抽出し、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類することを特徴とする解析装置を提供する。
本発明は第二の実施形態として、
画像取得装置と画像処理装置を含む解析装置であって、
前記画像取得装置は、
複数の細胞に関する自家蛍光画像を取得し、
前記画像処理装置は、
前記自家蛍光画像を用いて前記複数の細胞に関する自家蛍光情報を取得し、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得し、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、個々の細胞の自家蛍光情報を抽出し、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類することを特徴とする解析装置を提供する。
画像処理装置は、個々の細胞の領域を示す領域情報を、複数の細胞に関する自家蛍光画像に基づいて取得することができる。あるいは、画像取得装置は、さらに、複数の細胞に関する領域抽出用画像を取得することができる。その場合、画像処理装置は個々の細胞の領域を示す領域情報を、複数の細胞に関する領域抽出用画像を用いて取得することができる。
(装置の構成)
図1に、本実施形態における解析装置100の装置の構成の一例を示す。
図1に示すように、解析装置100は、画像取得装置1Aと画像処理装置2Aがケーブル3A等のインターフェースを介してデータ送受信可能に接続されて構成されている。なお、画像取得装置1Aと画像処理装置2Aとの接続方式は特に限定されない。例えば、画像取得装置と画像処理装置はLAN(Local Area Network)により接続されることとしてもよいし、無線により接続される構成としてもよい。
図1に、本実施形態における解析装置100の装置の構成の一例を示す。
図1に示すように、解析装置100は、画像取得装置1Aと画像処理装置2Aがケーブル3A等のインターフェースを介してデータ送受信可能に接続されて構成されている。なお、画像取得装置1Aと画像処理装置2Aとの接続方式は特に限定されない。例えば、画像取得装置と画像処理装置はLAN(Local Area Network)により接続されることとしてもよいし、無線により接続される構成としてもよい。
画像取得装置1Aは、載置ステージ上に載置された、細胞培養容器上の細胞あるいはスライドガラス上の組織切片等の画像を取得し、画像処理装置2Aに送信する。
画像取得装置1Aの構成の例を図2に示す。画像取得装置1Aは照射部11、結像部12、撮像部13、制御部14、通信I/F15、操作部16等、を備えて構成されている。照射部11は、蛍光観察用の光源フィルター等と細胞形態観察用の白色光源により構成され、載置ステージに載置された培養細胞あるいは組織切片等に光を照射する。結像部12は、接眼レンズ、対物レンズ、フィルター等により構成され、照射した光により細胞培養容器上の細胞あるいはスライド上の組織切片から発せられる透過光、反射光、または蛍光を結像する。撮像部13は、例えば、CMOS(Complementary MOS)センサー等を備え、結像部により結像面に結像される像を撮像して画像のデジタル画像データ(R、G、Bの画像データ)を生成するカメラである。制御部14は、CPU(Central Processing Unit)、RAM(Random Access Memory)等を備えて構成される。制御部14は内部に備え付けの各種プログラムとの協働により各種処理を実行する。通信I/F15は、生成された画像の画像データを画像処理装置2Aに送信する。本実施形態において、画像取得装置1Aは、明視野観察に適した照射手段及び結像手段を組み合わせた明視野ユニット、蛍光観察に適した照射手段及び結像手段を組み合わせた蛍光ユニットを備えることができる。ユニットを切り替えることにより明視野/蛍光を切り替える可能とすることができる。操作部16は、文字入力キー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードで押下操作されたキーの押下信号とマウスによる操作信号とを、入力信号として制御部14に出力する。
画像取得装置1Aの構成の例を図2に示す。画像取得装置1Aは照射部11、結像部12、撮像部13、制御部14、通信I/F15、操作部16等、を備えて構成されている。照射部11は、蛍光観察用の光源フィルター等と細胞形態観察用の白色光源により構成され、載置ステージに載置された培養細胞あるいは組織切片等に光を照射する。結像部12は、接眼レンズ、対物レンズ、フィルター等により構成され、照射した光により細胞培養容器上の細胞あるいはスライド上の組織切片から発せられる透過光、反射光、または蛍光を結像する。撮像部13は、例えば、CMOS(Complementary MOS)センサー等を備え、結像部により結像面に結像される像を撮像して画像のデジタル画像データ(R、G、Bの画像データ)を生成するカメラである。制御部14は、CPU(Central Processing Unit)、RAM(Random Access Memory)等を備えて構成される。制御部14は内部に備え付けの各種プログラムとの協働により各種処理を実行する。通信I/F15は、生成された画像の画像データを画像処理装置2Aに送信する。本実施形態において、画像取得装置1Aは、明視野観察に適した照射手段及び結像手段を組み合わせた明視野ユニット、蛍光観察に適した照射手段及び結像手段を組み合わせた蛍光ユニットを備えることができる。ユニットを切り替えることにより明視野/蛍光を切り替える可能とすることができる。操作部16は、文字入力キー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードで押下操作されたキーの押下信号とマウスによる操作信号とを、入力信号として制御部14に出力する。
画像処理装置2Aは、画像取得装置1Aから送信された画像を用いて、複数の細胞を2以上の群へ分類する。
図3に、画像処理装置2Aの構成の例を示す。画像処理装置2Aは、制御部21、操作部22、表示部23、通信I/F24、記憶部等25を備えて構成され、各部はバス26を介して接続されている。
図3に、画像処理装置2Aの構成の例を示す。画像処理装置2Aは、制御部21、操作部22、表示部23、通信I/F24、記憶部等25を備えて構成され、各部はバス26を介して接続されている。
制御部21は、CPU等を備えて構成される。記憶部25に記憶されている各種プログラムとの協働により各種処理を実行し、画像処理装置の動作を統括的に制御する。例えば、制御部21は、記憶部25に記憶されているプログラムとの協働により画像処理を実行する。制御部は、画像取得装置1Aから送信された画像を解析する。制御部は、領域情報と、複数の細胞に関する自家蛍光情報の情報から、個々の細胞の自家蛍光情報を抽出する。そして、これを用いて、複数の細胞を2以上の群に分類する。ただし、途中の工程、特に分類には、人が介在してもよい。
操作部22は、文字入力キー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードで押下操作されたキーの押下信号とマウスによる操作信号とを、入力信号として制御部21に出力する。
表示部23は、例えば、CRT(Cathode Ray Tube)やLCD(Liquid Crystal Display)等のモニタを備えて構成されており、制御部21から入力される表示信号の指示にしたがって、各種画面を表示する。本実施形態において、表示部23は、例えば、算出された特徴量、特徴量に基づいたグラフ、分類した結果、あるいは、分類した群を対応させた表現型等の表示のための手段として機能する。
通信I/F24は、画像取得装置1Aをはじめとする外部機器との間でデータ送受信を行なうためのインターフェースである。通信I/F24は、明視野画像と蛍光画像の入力手段として機能する。本実施形態において、通信I/F24は、入力手段として機能する。
記憶部25は、例えばHDD(Hard Disk Drive)や半導体の不揮発性メモリー等で構成されている。記憶部25には、前述のように各種プログラムや各種データ等が記憶されている。例えば、記憶部25は、後述する画像解析処理で使用される倍率テーブルをはじめとする各種データを記憶している。
その他、画像処理装置2Aは、LANアダプターやルーター等を備え、LAN等の通信ネットワークを介して外部機器と接続される構成としてもよい。
本実施形態における画像処理装置2Aは、画像取得装置1Aから送信された自家蛍光画像及び必要に応じ領域抽出用画像を用いて解析を行う。
自家蛍光画像は、対象である複数の細胞に対し、画像取得装置1Aにおいて所定波長の励起光を照射して蛍光物質を発光させ、この蛍光を光源波長以上のカットフィルターを通して拡大結像及び撮影することにより得られる画像である。自家蛍光画像に現れる蛍光は、細胞における、代謝状態等に対応する特定の生体物質の発現と状態を示すものである。
自家蛍光画像は、対象である複数の細胞に対し、画像取得装置1Aにおいて所定波長の励起光を照射して蛍光物質を発光させ、この蛍光を光源波長以上のカットフィルターを通して拡大結像及び撮影することにより得られる画像である。自家蛍光画像に現れる蛍光は、細胞における、代謝状態等に対応する特定の生体物質の発現と状態を示すものである。
画像取得装置1Aと画像処理装置2Aにおける工程の一例が図4に示される。
ステップP1において、画像取得装置は領域抽出用画像を取得する(本工程はなくともよい)。P2において、画像取得装置は、自家蛍光画像を取得する。ステップP1とステップP2はいずれが先であってもよく、前後関係は問わない。
ステップP1において、画像取得装置は領域抽出用画像を取得する(本工程はなくともよい)。P2において、画像取得装置は、自家蛍光画像を取得する。ステップP1とステップP2はいずれが先であってもよく、前後関係は問わない。
通信I/F24により、画像取得装置1Aからの領域抽出用画像及び自家蛍光画像が画像処理装置に入力される。
画像処理装置2AはステップP3において、領域抽出用画像あるいは自家蛍光画像を用いて領域情報を取得する。例えば、領域抽出用画像から領域を抽出し、ラベリング処理を行い、各細胞にラベルが付与された情報が取得される。領域情報が、自家蛍光画像を用いて取得される場合、ステップP1の領域抽出用画像の取得は不要である。
ステップP4において、自家蛍光画像を用いて、自家蛍光情報を取得する。例えば、自家蛍光画像についてRGB成分ごとの分離が行われる。
なお、ステップP3とP4はいずれが先であってもよく、前後関係は問わない。
ステップP5において、細胞ごとの自家蛍光情報の抽出が行われ、続くステップP6において、自家蛍光情報、及びその他の情報から特徴量が決定される。ステップP7において、細胞ごとの特徴量が出力される。この出力結果は画像処理装置2Aの表示部23に表示することが可能である。
画像処理装置2AはステップP3において、領域抽出用画像あるいは自家蛍光画像を用いて領域情報を取得する。例えば、領域抽出用画像から領域を抽出し、ラベリング処理を行い、各細胞にラベルが付与された情報が取得される。領域情報が、自家蛍光画像を用いて取得される場合、ステップP1の領域抽出用画像の取得は不要である。
ステップP4において、自家蛍光画像を用いて、自家蛍光情報を取得する。例えば、自家蛍光画像についてRGB成分ごとの分離が行われる。
なお、ステップP3とP4はいずれが先であってもよく、前後関係は問わない。
ステップP5において、細胞ごとの自家蛍光情報の抽出が行われ、続くステップP6において、自家蛍光情報、及びその他の情報から特徴量が決定される。ステップP7において、細胞ごとの特徴量が出力される。この出力結果は画像処理装置2Aの表示部23に表示することが可能である。
続く工程においては2通りのルートが存在する。解析の目的、用いた画像あるいは得られた特徴量等に応じて、以下の2つのルートのうち、より適切なものを、あるいは両方をコールアップして用いることができる。
ルート1を説明する。ステップP8において、ステップP7で出力した特徴量に対し、適当な特徴量を軸として選択し、各細胞を1つの点として1次元、2次元あるいは3次元のグラフとして可視化する。その後、得られたグラフから、目視で複数の群への細胞を分類する。それらを、さらに、細胞の表現型に対応させ、さらにその結果を表示部23に表示することができる。
ルート2について説明する。ステップP9において、ステップP7で出力した特徴量に関して、適当なものを選択し、教師なし機械学習を用いてクラスタリング処理を行い、細胞を分類する。それらを、さらに、細胞の表現型に対応させ、さらにその結果を表示部23に表示することができる。
ルート1を説明する。ステップP8において、ステップP7で出力した特徴量に対し、適当な特徴量を軸として選択し、各細胞を1つの点として1次元、2次元あるいは3次元のグラフとして可視化する。その後、得られたグラフから、目視で複数の群への細胞を分類する。それらを、さらに、細胞の表現型に対応させ、さらにその結果を表示部23に表示することができる。
ルート2について説明する。ステップP9において、ステップP7で出力した特徴量に関して、適当なものを選択し、教師なし機械学習を用いてクラスタリング処理を行い、細胞を分類する。それらを、さらに、細胞の表現型に対応させ、さらにその結果を表示部23に表示することができる。
以下、実施例により、本発明をさらに詳細に説明するが、本発明は以下に限定されるものではない。
<単一種類の細胞の、がん悪性度による分離>
(細胞の培養)
ヒト肺癌細胞 SK-LU1(ケーエーシー株式会社)を細胞培養用培地104(ケーエーシー株式会社)中で培養した。細胞培養容器としては35mm glass bottom dish(IWAKI)を用いた。
(細胞の培養)
ヒト肺癌細胞 SK-LU1(ケーエーシー株式会社)を細胞培養用培地104(ケーエーシー株式会社)中で培養した。細胞培養容器としては35mm glass bottom dish(IWAKI)を用いた。
(自家蛍光画像観察)
画像取得装置1Aを用いた。細胞培養液をPhosphate buffer saline(以下PBS、ThermoFisher社)に置換した。365nm LED(オプトコード株式会社製)の光源に対して400nmショートパスフィルターを介して照射して励起、430nmロングパスフィルターを通して撮像した自家蛍光画像を撮像した。また、450nm LED(オプトコード株式会社製)の光源に対して470nmショートパスフィルターを介して照射して励起し、530nmロングパスフィルターを通して自家蛍光画像を撮像した。
画像取得装置1Aを用いた。細胞培養液をPhosphate buffer saline(以下PBS、ThermoFisher社)に置換した。365nm LED(オプトコード株式会社製)の光源に対して400nmショートパスフィルターを介して照射して励起、430nmロングパスフィルターを通して撮像した自家蛍光画像を撮像した。また、450nm LED(オプトコード株式会社製)の光源に対して470nmショートパスフィルターを介して照射して励起し、530nmロングパスフィルターを通して自家蛍光画像を撮像した。
(免疫染色)
培養細胞の培養液を除去して直ちに、4%パラホルムアルデヒド水溶液(富士フイルム和光純薬社)を用いて固定化し、続いてTris buffer saline(以下TBS、ThermoFisher社)で3回洗浄をした。Triton-X 100(0.1%)を3分間処理して膜透過処理を行ったのち、TBSで洗浄した。Blocking one(ナカライテスク社)を用いてブロッキング処理したのち、腫瘍悪性化マーカーであるビメンチン(vimentin)に対して、anti-vimentin rabbit IgG(abcam社)で標識を行った。TBSで洗浄後、続いて核をHoechst色素(株式会社 同仁化学研究所)、ビメンチンをAlexa-fluor 555 anti-rabbit IgG(abcam社)を用いて、核とビメンチンを可視化した。
培養細胞の培養液を除去して直ちに、4%パラホルムアルデヒド水溶液(富士フイルム和光純薬社)を用いて固定化し、続いてTris buffer saline(以下TBS、ThermoFisher社)で3回洗浄をした。Triton-X 100(0.1%)を3分間処理して膜透過処理を行ったのち、TBSで洗浄した。Blocking one(ナカライテスク社)を用いてブロッキング処理したのち、腫瘍悪性化マーカーであるビメンチン(vimentin)に対して、anti-vimentin rabbit IgG(abcam社)で標識を行った。TBSで洗浄後、続いて核をHoechst色素(株式会社 同仁化学研究所)、ビメンチンをAlexa-fluor 555 anti-rabbit IgG(abcam社)を用いて、核とビメンチンを可視化した。
(染色画像観察)
核については、365nm LED(オプトコード株式会社製)の光源に対して400nmショートパスフィルターを介して照射して励起し、430nm-490nmバンドパスフィルターを通して撮像した核の蛍光画像を得た。ビメンチンについては、450nm LED(オプトコード株式会社製)の光源に対して470nmショートパスフィルターを介して照射して励起し、530nm-550nmバンドパスフィルターを通して撮像したビメンチン染色画像を撮像した。
核については、365nm LED(オプトコード株式会社製)の光源に対して400nmショートパスフィルターを介して照射して励起し、430nm-490nmバンドパスフィルターを通して撮像した核の蛍光画像を得た。ビメンチンについては、450nm LED(オプトコード株式会社製)の光源に対して470nmショートパスフィルターを介して照射して励起し、530nm-550nmバンドパスフィルターを通して撮像したビメンチン染色画像を撮像した。
(細胞の自家蛍光情報の抽出)
核領域の面積を取得し、ラベリングを行った。各自家蛍光画像について核画像とのマージを行い、個々の細胞ごとに核領域の自家蛍光成分を取得した。各々の自家蛍光画像のRGB成分をそれぞれ分離し、個々の細胞ごとに核領域の自家蛍光のR成分、G成分、B成分を取得した。それら成分ごとの詳細データ(核面積、核面積中の最小値、最大値、平均値、中央値等)を、ラベリングした細胞ごとに出力した。
核領域の面積を取得し、ラベリングを行った。各自家蛍光画像について核画像とのマージを行い、個々の細胞ごとに核領域の自家蛍光成分を取得した。各々の自家蛍光画像のRGB成分をそれぞれ分離し、個々の細胞ごとに核領域の自家蛍光のR成分、G成分、B成分を取得した。それら成分ごとの詳細データ(核面積、核面積中の最小値、最大値、平均値、中央値等)を、ラベリングした細胞ごとに出力した。
(自家蛍光輝度情報のプロット)
自家蛍光画像から取得した各細胞の輝度情報を用いて、以下の縦軸及び横軸で細胞ごとにプロットした。
横軸:350nm励起/430nmロングパスフィルター観察の自家蛍光画像のB成分輝度値の平均値
縦軸:450nm励起/530nmロングパスフィルター観察の自家蛍光画像のG成分の平均値
自家蛍光画像から取得した各細胞の輝度情報を用いて、以下の縦軸及び横軸で細胞ごとにプロットした。
横軸:350nm励起/430nmロングパスフィルター観察の自家蛍光画像のB成分輝度値の平均値
縦軸:450nm励起/530nmロングパスフィルター観察の自家蛍光画像のG成分の平均値
(ビメンチン染色画像の輝度情報抽出、サブタイプ分離、プロット)
自家蛍光画像でなく、染色画像に対しても、自家蛍光画像に対して行うのと同様に、本発明の実施形態の方法を適用した。すなわち、染色画像における個々の細胞の輝度値を記録し、核領域の面積を画像処理で取得した。ビメンチン染色画像について個々の細胞ごとに、核領域の色彩情報を取得した。自家蛍光画像の輝度情報抽出と同様に、細胞の位置データと細胞の詳細データ(核面積、核面積中の最小値、最大値、平均値、中央値)が一致するように記録した。その上で、ビメンチン輝度値により、陽性/陰性の判断を行った。
横軸として、350nm励起/430nmロングパスフィルター観察の自家蛍光画像のB成分輝度値の平均値、縦軸を450nm励起/530nmロングパスフィルター観察の自家蛍光画像のG成分の平均値とし、プロットした。
自家蛍光画像でなく、染色画像に対しても、自家蛍光画像に対して行うのと同様に、本発明の実施形態の方法を適用した。すなわち、染色画像における個々の細胞の輝度値を記録し、核領域の面積を画像処理で取得した。ビメンチン染色画像について個々の細胞ごとに、核領域の色彩情報を取得した。自家蛍光画像の輝度情報抽出と同様に、細胞の位置データと細胞の詳細データ(核面積、核面積中の最小値、最大値、平均値、中央値)が一致するように記録した。その上で、ビメンチン輝度値により、陽性/陰性の判断を行った。
横軸として、350nm励起/430nmロングパスフィルター観察の自家蛍光画像のB成分輝度値の平均値、縦軸を450nm励起/530nmロングパスフィルター観察の自家蛍光画像のG成分の平均値とし、プロットした。
(教師なし学習)
各細胞の輝度を説明変数とし、機械学習(k-means法)によりクラスターに分類することを試みた。
以下、各実験について記述する。なお、各試薬の調製等は上述の通りである。
各細胞の輝度を説明変数とし、機械学習(k-means法)によりクラスターに分類することを試みた。
以下、各実験について記述する。なお、各試薬の調製等は上述の通りである。
(SK-LU1細胞の自家蛍光の輝度値のプロット)
35mm glass bottom dish(IWAKI社製)にSK-LU1細胞(ケーエーシー株式会社)を培養した。画像取得装置を用いてディッシュ上のSK-LU1細胞の自家蛍光画像を撮像した。結果の一例を図5に示す。励起光を照射し、励起光が結像部に入らないよう選定したロングパスフィルターを介してRGBカラーセンサー搭載カメラで撮像した像を図5(a)に、解析のために撮像した画像をRGB別に16ビットグレースケール画像に分解した像を図5(b)に示す。本実施例では、RGB別に分解した像を自家蛍光画像として用いた。
続いて撮像後のディッシュ上のSK-LU1細胞の核とビメンチンを免疫染色により染色し、同様にディッシュ上で画像撮影装置を用いて染色画像を撮像した。ビメンチンは、がん細胞の上皮間葉系転換が生じた際に発現する分子マーカーである。
結果の一例を図6に示す。図6(a)は核染色して得られた画像、図6(b)はビメンチンを免疫染色により染色して得られた画像である。本実施例では、核染色して得られた画像を領域抽出用画像とした。
35mm glass bottom dish(IWAKI社製)にSK-LU1細胞(ケーエーシー株式会社)を培養した。画像取得装置を用いてディッシュ上のSK-LU1細胞の自家蛍光画像を撮像した。結果の一例を図5に示す。励起光を照射し、励起光が結像部に入らないよう選定したロングパスフィルターを介してRGBカラーセンサー搭載カメラで撮像した像を図5(a)に、解析のために撮像した画像をRGB別に16ビットグレースケール画像に分解した像を図5(b)に示す。本実施例では、RGB別に分解した像を自家蛍光画像として用いた。
続いて撮像後のディッシュ上のSK-LU1細胞の核とビメンチンを免疫染色により染色し、同様にディッシュ上で画像撮影装置を用いて染色画像を撮像した。ビメンチンは、がん細胞の上皮間葉系転換が生じた際に発現する分子マーカーである。
結果の一例を図6に示す。図6(a)は核染色して得られた画像、図6(b)はビメンチンを免疫染色により染色して得られた画像である。本実施例では、核染色して得られた画像を領域抽出用画像とした。
自家蛍光画像と領域抽出用画像を用いて、各細胞の核領域において、励起・観察波長別かつRGB別に、細胞の自家蛍光輝度を取得した。すなわち、領域抽出用画像を二値化し、核から推定した細胞位置情報を画像マスクとし、自家蛍光像との重ね合わせにより、個々の細胞の自家蛍光のRGB各色情報を抽出した。
個々の細胞の核領域内における自家蛍光輝度を平均したものを各細胞の平均輝度とした。
本実施例においては、各細胞の自家蛍光情報として、2種類の自家蛍光輝度を採用した。すなわち、365nm励起/430nmロングパスフィルター観察自家蛍光画像のB成分の平均輝度と450nm励起/530nmロングパスフィルター観察自家蛍光画像のG成分の平均輝度を、自家蛍光輝度値とした。
各細胞を以下の縦軸及び横軸でスキャッタープロットした。
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察自家蛍光画像のG成分の平均値
横軸:各細胞における365nm励起/430nmロングパスフィルター観察自家蛍光画像のB成分の平均値
結果を図7に示す。可視化を容易にするため対数プロットとした。
グラフの各点はそれぞれ細胞1つ1つの平均輝度値を示している。
個々の細胞の核領域内における自家蛍光輝度を平均したものを各細胞の平均輝度とした。
本実施例においては、各細胞の自家蛍光情報として、2種類の自家蛍光輝度を採用した。すなわち、365nm励起/430nmロングパスフィルター観察自家蛍光画像のB成分の平均輝度と450nm励起/530nmロングパスフィルター観察自家蛍光画像のG成分の平均輝度を、自家蛍光輝度値とした。
各細胞を以下の縦軸及び横軸でスキャッタープロットした。
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察自家蛍光画像のG成分の平均値
横軸:各細胞における365nm励起/430nmロングパスフィルター観察自家蛍光画像のB成分の平均値
結果を図7に示す。可視化を容易にするため対数プロットとした。
グラフの各点はそれぞれ細胞1つ1つの平均輝度値を示している。
(SK-LU1細胞の自家蛍光の輝度値のプロットの目視によるクラスタリング(図8))
図7で示したプロットは、目視により、図8に示すように少なくとも2つ(点線・実線)のクラスターに分離できる可能性を示した。細胞の発する自家蛍光は代謝に関与する補酵素(NADH,FAD)の蛍光に主に起因すると考えられる。個々の細胞の代謝経路の差によって、各活性化補酵素の割合が異なる。それら活性化補酵素の内に蛍光を発する補酵素があった場合、それぞれの補酵素の蛍光波長が異なるため、異なる代謝経路を持つ細胞では自家蛍光の色合いに差異が生じ、本実施形態のように非侵襲な方法で分類することができると考えられる。
図7で示したプロットは、目視により、図8に示すように少なくとも2つ(点線・実線)のクラスターに分離できる可能性を示した。細胞の発する自家蛍光は代謝に関与する補酵素(NADH,FAD)の蛍光に主に起因すると考えられる。個々の細胞の代謝経路の差によって、各活性化補酵素の割合が異なる。それら活性化補酵素の内に蛍光を発する補酵素があった場合、それぞれの補酵素の蛍光波長が異なるため、異なる代謝経路を持つ細胞では自家蛍光の色合いに差異が生じ、本実施形態のように非侵襲な方法で分類することができると考えられる。
(sk-lu1細胞のビメンチン輝度値によるサブタイプの分離)
上述の通り撮影した図6(b)のビメンチン染色画像とそこから取得した個々の細胞のビメンチン輝度値データから、細胞をビメンチンの陽性/陰性の2つのサブタイプに分類した。それら2つのサブタイプを別の形状として、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察自家蛍光画像のG成分の平均値
横軸:各細胞における365nm励起/430nmロングパスフィルター観察自家蛍光画像のB成分の平均値
結果を図9に示す。図中subtype2として黒四角で示される各点がビメンチン陽性細胞であり、subtype1として黒丸で示される各点がビメンチン陰性細胞である。
上述の通り撮影した図6(b)のビメンチン染色画像とそこから取得した個々の細胞のビメンチン輝度値データから、細胞をビメンチンの陽性/陰性の2つのサブタイプに分類した。それら2つのサブタイプを別の形状として、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察自家蛍光画像のG成分の平均値
横軸:各細胞における365nm励起/430nmロングパスフィルター観察自家蛍光画像のB成分の平均値
結果を図9に示す。図中subtype2として黒四角で示される各点がビメンチン陽性細胞であり、subtype1として黒丸で示される各点がビメンチン陰性細胞である。
(図8、図9のマージ)
図8と図9の結果をマージしたものを図10に示す。図8の2つのクラスターと、図9のビメンチン陽性/陰性の判断による2つのサブタイプの主要なプロット位置が重なっている。図中、実線で囲まれたクラスターに多くのビメンチン陽性細胞のサブタイプが含まれた。一方、図中、点線で囲まれたクラスターにはほとんどビメンチン陽性細胞のサブタイプが含まれていなかった。すなわち、自家蛍光画像による細胞表現型の分類結果はビメンチン陽性/陰性と高く相関することが示された。ビメンチンは、従来は、特異的に染色しなければ、その陽性/陰性の情報を得ることはできない。すなわち、本実施形態の方法により、自家蛍光を用いて、染色することなく、ビメンチン陽性/陰性と同程度に、分類ができたことが示された。さらには、ビメンチンはがん細胞の上皮間葉系転換が生じた際に発現する分子マーカーであり、腫瘍悪性化マーカーの一つの指標になりえる。すなわち、本実施形態の方法により、1視野中に“がん悪性度の異なる細胞”という複数の異なる代謝経路を持つ細胞種が存在していた場合、自家蛍光を用いて、染色することなく、腫瘍悪性化の有無、あるいは、上皮間葉系転換の有無により、細胞が分類できたことが示された。
図8と図9の結果をマージしたものを図10に示す。図8の2つのクラスターと、図9のビメンチン陽性/陰性の判断による2つのサブタイプの主要なプロット位置が重なっている。図中、実線で囲まれたクラスターに多くのビメンチン陽性細胞のサブタイプが含まれた。一方、図中、点線で囲まれたクラスターにはほとんどビメンチン陽性細胞のサブタイプが含まれていなかった。すなわち、自家蛍光画像による細胞表現型の分類結果はビメンチン陽性/陰性と高く相関することが示された。ビメンチンは、従来は、特異的に染色しなければ、その陽性/陰性の情報を得ることはできない。すなわち、本実施形態の方法により、自家蛍光を用いて、染色することなく、ビメンチン陽性/陰性と同程度に、分類ができたことが示された。さらには、ビメンチンはがん細胞の上皮間葉系転換が生じた際に発現する分子マーカーであり、腫瘍悪性化マーカーの一つの指標になりえる。すなわち、本実施形態の方法により、1視野中に“がん悪性度の異なる細胞”という複数の異なる代謝経路を持つ細胞種が存在していた場合、自家蛍光を用いて、染色することなく、腫瘍悪性化の有無、あるいは、上皮間葉系転換の有無により、細胞が分類できたことが示された。
(自家蛍光強度を説明変数とした機械学習によるクラスタリング)
上述のとおり撮影した自家蛍光画像から365nm励起/430nmロングパスフィルター観察のB成分輝度値の最小値、縦軸を450nm励起/530nmロングパスフィルター観察のG成分輝度値の最小値の2つを説明変数として取得した。この説明変数を用いて教師なし機械学習により細胞を2つのクラスターに分離した。2つのクラスターをそれぞれ黒四角、黒丸として、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
横軸:各細胞における365nm励起/430nmロングパスフィルター観察自家蛍光画像のB成分の平均値
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察自家蛍光画像のG成分の平均値
結果を図11に示す。
この結果と図9を比較することで、ビメンチン陽性/陰性から判断した2つのサブタイプと、機械学習により分離した、2つのクラスターの位置はおおむね一致していることが示された。すなわち、機械学習を用いた場合においても、本実施形態の方法により、自家蛍光を用いて、染色することなく、ビメンチン陽性/陰性と同程度に、細胞の分類ができたことが示された。さらに、機械学習を用いた場合において、自家蛍光を用いて、染色することなく、腫瘍悪性化の有無、あるいは、上皮間葉系転換の有無により、細胞が分類できたことが示された。
上述のとおり撮影した自家蛍光画像から365nm励起/430nmロングパスフィルター観察のB成分輝度値の最小値、縦軸を450nm励起/530nmロングパスフィルター観察のG成分輝度値の最小値の2つを説明変数として取得した。この説明変数を用いて教師なし機械学習により細胞を2つのクラスターに分離した。2つのクラスターをそれぞれ黒四角、黒丸として、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
横軸:各細胞における365nm励起/430nmロングパスフィルター観察自家蛍光画像のB成分の平均値
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察自家蛍光画像のG成分の平均値
結果を図11に示す。
この結果と図9を比較することで、ビメンチン陽性/陰性から判断した2つのサブタイプと、機械学習により分離した、2つのクラスターの位置はおおむね一致していることが示された。すなわち、機械学習を用いた場合においても、本実施形態の方法により、自家蛍光を用いて、染色することなく、ビメンチン陽性/陰性と同程度に、細胞の分類ができたことが示された。さらに、機械学習を用いた場合において、自家蛍光を用いて、染色することなく、腫瘍悪性化の有無、あるいは、上皮間葉系転換の有無により、細胞が分類できたことが示された。
<2種類のがん細胞における、細胞種による分離>
はじめに、各手順について説明する。
はじめに、各手順について説明する。
(細胞の培養)
ヒト肺癌細胞 PC-9(ケーエーシー株式会社)と、ヒト肺癌細胞 NCI-H322(ケーエーシー株式会社)を、細胞培養用培地103(ケーエーシー株式会社)中で混合培養した。細胞培養容器としては35mm glass bottom dish(IWAKI)を用いた。
ヒト肺癌細胞 PC-9(ケーエーシー株式会社)と、ヒト肺癌細胞 NCI-H322(ケーエーシー株式会社)を、細胞培養用培地103(ケーエーシー株式会社)中で混合培養した。細胞培養容器としては35mm glass bottom dish(IWAKI)を用いた。
(自家蛍光画像観察)
画像取得装置1Aを用いた。細胞培養液をPhosphate buffer saline(以下PBS、ThermoFisher社)に置換した。365nm LED(オプトコード株式会社製)の光源に対して400nmショートパスフィルターを介して照射して励起、430nmロングパスフィルターを通して撮像した自家蛍光画像を撮像した。また、450nm LED(オプトコード株式会社製)の光源に対して470nmショートパスフィルターを介して照射して励起し、530nmロングパスフィルターを通して自家蛍光画像を撮像した。
画像取得装置1Aを用いた。細胞培養液をPhosphate buffer saline(以下PBS、ThermoFisher社)に置換した。365nm LED(オプトコード株式会社製)の光源に対して400nmショートパスフィルターを介して照射して励起、430nmロングパスフィルターを通して撮像した自家蛍光画像を撮像した。また、450nm LED(オプトコード株式会社製)の光源に対して470nmショートパスフィルターを介して照射して励起し、530nmロングパスフィルターを通して自家蛍光画像を撮像した。
(ミトコンドリア染色、核染色)
培養細胞に対して、ミトコンドリアをMitoPT(商標)-TMREアッセイキット(Immunochemistry Technologies社)、並びに核をHoechst色素(株式会社 同仁化学研究所)で染色し、核とミトコンドリア膜電位差を可視化した。
培養細胞に対して、ミトコンドリアをMitoPT(商標)-TMREアッセイキット(Immunochemistry Technologies社)、並びに核をHoechst色素(株式会社 同仁化学研究所)で染色し、核とミトコンドリア膜電位差を可視化した。
(染色画像観察)
核については、365nm LED(オプトコード株式会社製)の光源に対して400nmショートパスフィルターを介して照射して励起し、430nm-490nmバンドパスフィルターを通して撮像した核の蛍光画像を得た。ミトコンドリア膜電位差については、530nm LED(オプトコード株式会社製)の光源に対して550nmショートパスフィルターを介して照射して励起し、570nm-630nmバンドパスフィルターを通して撮像したミトコンドリア染色画像を得た。
核については、365nm LED(オプトコード株式会社製)の光源に対して400nmショートパスフィルターを介して照射して励起し、430nm-490nmバンドパスフィルターを通して撮像した核の蛍光画像を得た。ミトコンドリア膜電位差については、530nm LED(オプトコード株式会社製)の光源に対して550nmショートパスフィルターを介して照射して励起し、570nm-630nmバンドパスフィルターを通して撮像したミトコンドリア染色画像を得た。
(細胞の自家蛍光情報の抽出)
核領域の面積を取得し、ラベリングを行った。各自家蛍光画像について核画像とのマージを行い、個々の細胞ごとに核領域の自家蛍光成分を取得した。各々の自家蛍光画像のRGB成分をそれぞれ分離し、個々の細胞ごとに核領域の自家蛍光のR成分、G成分、B成分を取得した。また、自家蛍光のRGB色空間を、HLS色空間、HSV色空間、Lab色空間に変換し、色空間ごとに各成分を取得した。それら成分ごとの詳細データ(核面積、核面積中の最小値、最大値、平均値、中央値等)を、ラベリングした細胞ごとに出力した。
核領域の面積を取得し、ラベリングを行った。各自家蛍光画像について核画像とのマージを行い、個々の細胞ごとに核領域の自家蛍光成分を取得した。各々の自家蛍光画像のRGB成分をそれぞれ分離し、個々の細胞ごとに核領域の自家蛍光のR成分、G成分、B成分を取得した。また、自家蛍光のRGB色空間を、HLS色空間、HSV色空間、Lab色空間に変換し、色空間ごとに各成分を取得した。それら成分ごとの詳細データ(核面積、核面積中の最小値、最大値、平均値、中央値等)を、ラベリングした細胞ごとに出力した。
(自家蛍光輝度情報のプロット)
自家蛍光画像から取得した各細胞の輝度情報を用いて、以下の縦軸及び横軸で細胞ごとにプロットした(図12)。
横軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像、HSV色空間におけるH成分の平均値
自家蛍光画像から取得した各細胞の輝度情報を用いて、以下の縦軸及び横軸で細胞ごとにプロットした(図12)。
横軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像、HSV色空間におけるH成分の平均値
(細胞種ごとの自家蛍光輝度情報のプロット)
自家蛍光画像から取得した各細胞の輝度情報について、各細胞がPC-9あるいはNCI-H322のどちらであるかがラベリングされている。このラベルを用いて、以下の縦軸及び横軸で細胞種ごとに異なる色でプロットした(図13)。
横軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像、HSV色空間におけるH成分の平均値
自家蛍光画像から取得した各細胞の輝度情報について、各細胞がPC-9あるいはNCI-H322のどちらであるかがラベリングされている。このラベルを用いて、以下の縦軸及び横軸で細胞種ごとに異なる色でプロットした(図13)。
横軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像、HSV色空間におけるH成分の平均値
(ミトコンドリア染色画像の輝度情報抽出、プロット)
自家蛍光画像でなく、染色画像に対しても、自家蛍光画像に対して行うのと同様に、本発明の実施形態の方法を適用した。すなわち、染色画像における個々の細胞の輝度値を記録し、核領域の面積を画像処理で取得した。ミトコンドリア染色画像について個々の細胞ごとに、核領域の色彩情報を取得した。自家蛍光画像の輝度情報抽出と同様に、細胞の位置データと細胞の詳細データ(核面積、核面積中の最小値、最大値、平均値、中央値)が一致するように記録した。その上で、ミトコンドリア輝度値により、ミトコンドリア活性の大小の判断を行い、ミトコンドリアの活性が高いエリアと低いエリアを異なる色を用いて以下の縦軸及び横軸でプロットした(図14)。
横軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像、HSV色空間におけるH成分の平均値
自家蛍光画像でなく、染色画像に対しても、自家蛍光画像に対して行うのと同様に、本発明の実施形態の方法を適用した。すなわち、染色画像における個々の細胞の輝度値を記録し、核領域の面積を画像処理で取得した。ミトコンドリア染色画像について個々の細胞ごとに、核領域の色彩情報を取得した。自家蛍光画像の輝度情報抽出と同様に、細胞の位置データと細胞の詳細データ(核面積、核面積中の最小値、最大値、平均値、中央値)が一致するように記録した。その上で、ミトコンドリア輝度値により、ミトコンドリア活性の大小の判断を行い、ミトコンドリアの活性が高いエリアと低いエリアを異なる色を用いて以下の縦軸及び横軸でプロットした(図14)。
横軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
縦軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像、HSV色空間におけるH成分の平均値
(教師なし学習)
各細胞の輝度の一部を説明変数とし、機械学習(k-means法)によりクラスターに分類することを試みた(図15)。
以下、各実験について記述する。なお、各試薬の調製等は上述の通りである。
各細胞の輝度の一部を説明変数とし、機械学習(k-means法)によりクラスターに分類することを試みた(図15)。
以下、各実験について記述する。なお、各試薬の調製等は上述の通りである。
(PC-9細胞とNCI-H322細胞の自家蛍光撮像、ミトコンドリア染色像撮像)
35mm glass bottom dish(IWAKI社製)にPC-9細胞(ケーエーシー株式会社)とNCI-H322細胞(ケーエーシー株式会社)を混合培養した。画像取得装置を用いてディッシュ上の細胞の自家蛍光画像を撮像した。具体的には励起光を照射し、励起光が結像部に入らないよう選定したロングパスフィルターを介してRGBカラーセンサー搭載カメラで撮像した。本実施例では画像を、RGB別に分解した像を自家蛍光画像として用いた。
続いて撮像後のディッシュ上の細胞の核とミトコンドリアを染色し、同様にディッシュ上で画像撮影装置を用いて染色画像を撮像した。ミトコンドリア染色試薬は、ミトコンドリア膜電位に応じて染色され、ミトコンドリアにおける代謝状態を反映する。なお、ミトコンドリアが活発な代謝状態であるほど、強く染色される。
本実施例では、核染色して得られた画像を領域抽出用画像とした。
35mm glass bottom dish(IWAKI社製)にPC-9細胞(ケーエーシー株式会社)とNCI-H322細胞(ケーエーシー株式会社)を混合培養した。画像取得装置を用いてディッシュ上の細胞の自家蛍光画像を撮像した。具体的には励起光を照射し、励起光が結像部に入らないよう選定したロングパスフィルターを介してRGBカラーセンサー搭載カメラで撮像した。本実施例では画像を、RGB別に分解した像を自家蛍光画像として用いた。
続いて撮像後のディッシュ上の細胞の核とミトコンドリアを染色し、同様にディッシュ上で画像撮影装置を用いて染色画像を撮像した。ミトコンドリア染色試薬は、ミトコンドリア膜電位に応じて染色され、ミトコンドリアにおける代謝状態を反映する。なお、ミトコンドリアが活発な代謝状態であるほど、強く染色される。
本実施例では、核染色して得られた画像を領域抽出用画像とした。
(自家蛍光の輝度値のプロット)
自家蛍光画像と領域抽出用画像を用いて、各細胞の核領域において、励起・観察波長別かつRGB別に、細胞の自家蛍光輝度を取得した。すなわち、領域抽出用画像を二値化し、核から推定した細胞位置情報を画像マスクとし、自家蛍光像との重ね合わせにより、個々の細胞の自家蛍光のRGB各色情報を抽出した。続いて自家蛍光のRGB各色情報に対して色空間の変換を行うことで、RGB色空間以外に、HSV色空間、HLS色空間、Lab色空間における自家蛍光輝度を取得した。
個々の細胞の核領域内における自家蛍光輝度を平均したものを各細胞の自家蛍光輝度とした。
本実施例においては、各細胞の自家蛍光情報として、12種類の自家蛍光輝度を採用した。特に本実施形態においては、すなわち、365nm励起/430nmロングパスフィルター観察自家蛍光画像において、B成分(RGB色空間)、G成分(RGB色空間)、H成分(HSV色空間)、S成分(HSV色空間)、a成分(Lab色空間)、b成分(Lab色空間)、続いて、450nm励起/530nmロングパスフィルター観察自家蛍光画像において、G成分(RGB色空間)、R成分(RGB色空間)、H成分(HSV色空間)、S成分(HSV色空間)、a成分(Lab色空間)、b成分(Lab色空間)の平均輝度を、自家蛍光輝度値とした。
各細胞について、一例として以下の縦軸及び横軸でスキャッタープロットした。
縦軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
横軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像の、HSV色空間におけるH成分の平均値
結果を図12に示す。
グラフの各点はそれぞれ細胞1つ1つの平均輝度値を示している。
自家蛍光画像と領域抽出用画像を用いて、各細胞の核領域において、励起・観察波長別かつRGB別に、細胞の自家蛍光輝度を取得した。すなわち、領域抽出用画像を二値化し、核から推定した細胞位置情報を画像マスクとし、自家蛍光像との重ね合わせにより、個々の細胞の自家蛍光のRGB各色情報を抽出した。続いて自家蛍光のRGB各色情報に対して色空間の変換を行うことで、RGB色空間以外に、HSV色空間、HLS色空間、Lab色空間における自家蛍光輝度を取得した。
個々の細胞の核領域内における自家蛍光輝度を平均したものを各細胞の自家蛍光輝度とした。
本実施例においては、各細胞の自家蛍光情報として、12種類の自家蛍光輝度を採用した。特に本実施形態においては、すなわち、365nm励起/430nmロングパスフィルター観察自家蛍光画像において、B成分(RGB色空間)、G成分(RGB色空間)、H成分(HSV色空間)、S成分(HSV色空間)、a成分(Lab色空間)、b成分(Lab色空間)、続いて、450nm励起/530nmロングパスフィルター観察自家蛍光画像において、G成分(RGB色空間)、R成分(RGB色空間)、H成分(HSV色空間)、S成分(HSV色空間)、a成分(Lab色空間)、b成分(Lab色空間)の平均輝度を、自家蛍光輝度値とした。
各細胞について、一例として以下の縦軸及び横軸でスキャッタープロットした。
縦軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
横軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像の、HSV色空間におけるH成分の平均値
結果を図12に示す。
グラフの各点はそれぞれ細胞1つ1つの平均輝度値を示している。
(細胞の正解ラベルごとのプロット)
撮影した自家蛍光画像において、自家蛍光情報はPC-9細胞とNCI-H322細胞のラベルとともに取得されている。それら2つの細胞種を別の色として、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
縦軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
横軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像の、HSV色空間におけるH成分の平均値
結果を図13に示す。図中濃いグレーで示される各点がPC-9細胞であり、薄いグレーで示される各点がNCI-H322細胞である。図に示す通り、PC-9細胞とNCI-H322細胞は異なる自家蛍光の傾向を持つことが示されている。
前述のとおり、細胞の発する自家蛍光は主たる要素は代謝に関与する補酵素(NADH,FAD)の蛍光に起因すると考えられる。個々の細胞の代謝経路の差によって、各活性化補酵素の割合が異なる。それら活性化補酵素中に、蛍光を発する補酵素があった場合、それぞれの補酵素の蛍光波長が異なるため、異なる代謝経路を持つ細胞は自家蛍光の色合いに差異が生じる。本実施形態では、1視野中に異なる代謝経路を持つ2種の細胞種を混在させたため、それらの自家蛍光が異なる傾向を持つプロットとして出力されたと考えられる。
撮影した自家蛍光画像において、自家蛍光情報はPC-9細胞とNCI-H322細胞のラベルとともに取得されている。それら2つの細胞種を別の色として、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
縦軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
横軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像の、HSV色空間におけるH成分の平均値
結果を図13に示す。図中濃いグレーで示される各点がPC-9細胞であり、薄いグレーで示される各点がNCI-H322細胞である。図に示す通り、PC-9細胞とNCI-H322細胞は異なる自家蛍光の傾向を持つことが示されている。
前述のとおり、細胞の発する自家蛍光は主たる要素は代謝に関与する補酵素(NADH,FAD)の蛍光に起因すると考えられる。個々の細胞の代謝経路の差によって、各活性化補酵素の割合が異なる。それら活性化補酵素中に、蛍光を発する補酵素があった場合、それぞれの補酵素の蛍光波長が異なるため、異なる代謝経路を持つ細胞は自家蛍光の色合いに差異が生じる。本実施形態では、1視野中に異なる代謝経路を持つ2種の細胞種を混在させたため、それらの自家蛍光が異なる傾向を持つプロットとして出力されたと考えられる。
(ミトコンドリア染色画像の輝度情報抽出、プロット)
上述のとおり撮影したミトコンドリア染色画像とそこから取得した個々の細胞のミトコンドリア染色輝度値データから、細胞をミトコンドリア活性が強いエリアと弱いエリアの2つのエリアに分割した。それら2つのエリアを別の色で、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
縦軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
横軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像の、HSV色空間におけるH成分の平均値
結果を図14に示す。図中濃いグレーで示される各点がミトコンドリアの活性が相対的に高い細胞であり、一方で薄いグレーで示される各点がミトコンドリアの活性が相対的に低い細胞である。
図13と比較すると、PC-9細胞はミトコンドリアの活性が高いエリアに多く含まれている一方、NCI-H322細胞はミトコンドリアの活性が低いエリアに多く含まれている。
すなわち、細胞種によってミトコンドリアの活性状態は異なり、結果として、自家蛍光画像による細胞表現型の分類結果はミトコンドリアの活性と高く相関することが示された。
上述のとおり撮影したミトコンドリア染色画像とそこから取得した個々の細胞のミトコンドリア染色輝度値データから、細胞をミトコンドリア活性が強いエリアと弱いエリアの2つのエリアに分割した。それら2つのエリアを別の色で、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
縦軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、Lab色空間におけるb成分の平均値
横軸:各細胞における450nm励起/530nmロングパスフィルター観察の自家蛍光画像の、HSV色空間におけるH成分の平均値
結果を図14に示す。図中濃いグレーで示される各点がミトコンドリアの活性が相対的に高い細胞であり、一方で薄いグレーで示される各点がミトコンドリアの活性が相対的に低い細胞である。
図13と比較すると、PC-9細胞はミトコンドリアの活性が高いエリアに多く含まれている一方、NCI-H322細胞はミトコンドリアの活性が低いエリアに多く含まれている。
すなわち、細胞種によってミトコンドリアの活性状態は異なり、結果として、自家蛍光画像による細胞表現型の分類結果はミトコンドリアの活性と高く相関することが示された。
(自家蛍光強度を説明変数とした機械学習によるクラスタリング)
上述のとおり撮影した自家蛍光画像から、365nm励起/430nmロングパスフィルター観察自家蛍光画像において、B成分(RGB色空間)、G成分(RGB色空間)、H成分(HSV色空間)、S成分(HSV色空間)、a成分(Lab色空間)、b成分(Lab色空間)、続いて、450nm励起/530nmロングパスフィルター観察自家蛍光画像において、G成分(RGB色空間)、R成分(RGB色空間)、H成分(HSV色空間)、S成分(HSV色空間)、a成分(Lab色空間)、b成分(Lab色空間)の平均輝度、合計12個を説明変数として取得した。これらの説明変数を用いて教師なし機械学習(k-means法)により細胞を2つのクラスターに分離した。2つのクラスターをそれぞれ濃いグレー、薄いグレーとして、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
縦軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、LAB色空間におけるB成分の平均値
横軸:各細胞における450NM励起/530NMロングパスフィルター観察の自家蛍光画像の、HSV色空間におけるH成分の平均値
結果を図15に示す。
図13で示したプロットと図15で示したプロットを比較することで、自家蛍光を用いた教師なし機械学習によって分離した各クラスターと、正解ラベルによるPC-9細胞あるいはNCI-H322細胞ごとのクラスターがおおむね一致している。また、教師なし機械学習による正解率は72%だった。
以上、教師なし機械学習によりクラスター分離を行うことによって、ある程度の精度で、2種類の細胞種の分類ができた。また、ミトコンドリア活性は、従来は、特異的に染色しなければその大小の情報を得ることはできない。すなわち、本実施形態の方法により、機械学習を用いて、自家蛍光を用いて、染色することなく、ミトコンドリア染色と同程度に、細胞種の分類ができたことが示された。
上述のとおり撮影した自家蛍光画像から、365nm励起/430nmロングパスフィルター観察自家蛍光画像において、B成分(RGB色空間)、G成分(RGB色空間)、H成分(HSV色空間)、S成分(HSV色空間)、a成分(Lab色空間)、b成分(Lab色空間)、続いて、450nm励起/530nmロングパスフィルター観察自家蛍光画像において、G成分(RGB色空間)、R成分(RGB色空間)、H成分(HSV色空間)、S成分(HSV色空間)、a成分(Lab色空間)、b成分(Lab色空間)の平均輝度、合計12個を説明変数として取得した。これらの説明変数を用いて教師なし機械学習(k-means法)により細胞を2つのクラスターに分離した。2つのクラスターをそれぞれ濃いグレー、薄いグレーとして、以下を縦軸及び横軸として、それぞれスキャッタープロットした。
縦軸:各細胞における365nm励起/430nmロングパスフィルター観察の自家蛍光画像の、LAB色空間におけるB成分の平均値
横軸:各細胞における450NM励起/530NMロングパスフィルター観察の自家蛍光画像の、HSV色空間におけるH成分の平均値
結果を図15に示す。
図13で示したプロットと図15で示したプロットを比較することで、自家蛍光を用いた教師なし機械学習によって分離した各クラスターと、正解ラベルによるPC-9細胞あるいはNCI-H322細胞ごとのクラスターがおおむね一致している。また、教師なし機械学習による正解率は72%だった。
以上、教師なし機械学習によりクラスター分離を行うことによって、ある程度の精度で、2種類の細胞種の分類ができた。また、ミトコンドリア活性は、従来は、特異的に染色しなければその大小の情報を得ることはできない。すなわち、本実施形態の方法により、機械学習を用いて、自家蛍光を用いて、染色することなく、ミトコンドリア染色と同程度に、細胞種の分類ができたことが示された。
本発明の実施形態は以下の方法及び構成を含む。
(方法1)
複数の細胞に関する自家蛍光画像を用いて、前記複数の細胞に関する自家蛍光情報を取得する工程と、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得する工程と、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、前記個々の細胞の自家蛍光情報を抽出する工程と、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類する工程、を含む方法。
(方法2)
前記領域情報は前記複数の細胞に関する自家蛍光画像を用いて取得することを特徴とする方法1に記載の方法。
(方法3)
さらに、前記複数の細胞に関する領域抽出用画像を取得する工程を含み、
前記領域情報は、前記複数の細胞に関する領域抽出用画像を用いて取得することを特徴とする方法1に記載の方法。
(方法4)
前記自家蛍光情報は、自家蛍光の輝度を含む方法1から3のいずれかに記載の方法。
(方法5)
前記輝度は、RGBの色空間におけるR、G、及びBの少なくともいずれかの指標として算出される値である方法4に記載の方法。
(方法6)
前記輝度は、Lab色空間またはHSV色空間における指標として算出される値である方法4に記載の方法。
(方法7)
前記分類する工程では、少なくとも前記個々の細胞の自家蛍光の輝度を1つの軸とする、1次元、2次元あるいは3次元のグラフとして可視化することを特徴とする方法1から6のいずれか1項に記載の方法。
(方法8)
前記分類する工程では、少なくとも前記個々の細胞の自家蛍光情報を1つの特徴量として、主成分分析、k-means法、GMMからなる群より選ばれるいずれかの教師なし機械学習を行うことを特徴とする方法1から6のいずれか1項に記載の方法。
(方法9)
さらに、前記2以上の群を、群ごとに細胞の表現型に対応させる工程を含む方法1から8のいずれか1項に記載の方法。
(方法10)
前記細胞の表現型は、がん細胞と非がん細胞を含む組合せ、未分化細胞と分化細胞を含む組合せ、分化細胞と脱分化細胞を含む組合せ、及び同種がん細胞の上皮系と間葉系細胞を含む組合せ、から選ばれる少なくとも1の組合せを含む方法9に記載の方法。
(方法11)
前記複数の細胞に関する自家蛍光画像は、60mm2以上の領域を撮像して取得されることを特徴とする方法1から10のいずれかに記載の方法。
(方法12)
前記複数の細胞に関する自家蛍光画像は、少なくとも30000画素以上の画素数と、
かつ1画素あたり20μm以下の解像度を有することを特徴とする方法1から11のいずれかに記載の方法。
(方法13)
さらに、前記分類する工程によって分けられた群に関する定量的な情報を取得する工程を含む、方法1乃至12のいずれかに記載の方法。
(構成1)
画像取得装置と画像処理装置を含む解析装置であって、
前記画像取得装置は、
複数の細胞に関する自家蛍光画像を取得し、
前記画像処理装置は、
前記自家蛍光画像を用いて前記複数の細胞に関する自家蛍光情報を取得し、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得し、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、個々の細胞の自家蛍光情報を抽出し、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類することを特徴とする解析装置。
(構成2)
前記画像処理装置は、前記領域情報を前記複数の細胞に関する自家蛍光画像を用いて取得することを特徴とする構成1に記載の解析装置。
(構成3)
さらに、前記画像取得装置は、
前記複数の細胞に関する領域抽出用画像を取得し、
前記画像処理装置は、
前記領域情報を前記複数の細胞に関する領域抽出用画像を用いて取得することを特徴とする構成1に記載の解析装置。
(構成4)
前記画像処理装置が、前記分類をする際、少なくとも前記個々の細胞の自家蛍光の輝度を1つの軸とする、1次元、2次元あるいは3次元のグラフとして可視化することを特徴とする構成1から3のいずれか1項に記載の解析装置。
(構成5)
前記画像処理装置が、前記分類する際、少なくとも前記個々の細胞の自家蛍光情報を1つの特徴量として、主成分分析、k-means法、GMMからなる群より選ばれるいずれかの教師なし機械学習を行うことを特徴とする構成1から4のいずれかに記載の解析装置。
(構成6)
さらに、前記画像処理装置は、前記2以上の群を、群ごとに細胞の表現型に対応させることを特徴とする構成1から5のいずれかに記載の解析装置。
(構成7)
前記細胞の表現型は、がん細胞と非がん細胞を含む組合せ、未分化細胞と分化細胞を含む組合せ、分化細胞と脱分化細胞を含む組合せ、及び同種がん細胞の上皮系と間葉系細胞を含む組合せ、から選ばれる少なくとも1の組合せを含む構成6に記載の解析装置。
本発明は上記実施の形態に制限されるものではなく、本発明の精神及び範囲から離脱することなく、様々な変更及び変形が可能である。従って、本発明の範囲を公にするために以下の請求項を添付する。
(方法1)
複数の細胞に関する自家蛍光画像を用いて、前記複数の細胞に関する自家蛍光情報を取得する工程と、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得する工程と、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、前記個々の細胞の自家蛍光情報を抽出する工程と、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類する工程、を含む方法。
(方法2)
前記領域情報は前記複数の細胞に関する自家蛍光画像を用いて取得することを特徴とする方法1に記載の方法。
(方法3)
さらに、前記複数の細胞に関する領域抽出用画像を取得する工程を含み、
前記領域情報は、前記複数の細胞に関する領域抽出用画像を用いて取得することを特徴とする方法1に記載の方法。
(方法4)
前記自家蛍光情報は、自家蛍光の輝度を含む方法1から3のいずれかに記載の方法。
(方法5)
前記輝度は、RGBの色空間におけるR、G、及びBの少なくともいずれかの指標として算出される値である方法4に記載の方法。
(方法6)
前記輝度は、Lab色空間またはHSV色空間における指標として算出される値である方法4に記載の方法。
(方法7)
前記分類する工程では、少なくとも前記個々の細胞の自家蛍光の輝度を1つの軸とする、1次元、2次元あるいは3次元のグラフとして可視化することを特徴とする方法1から6のいずれか1項に記載の方法。
(方法8)
前記分類する工程では、少なくとも前記個々の細胞の自家蛍光情報を1つの特徴量として、主成分分析、k-means法、GMMからなる群より選ばれるいずれかの教師なし機械学習を行うことを特徴とする方法1から6のいずれか1項に記載の方法。
(方法9)
さらに、前記2以上の群を、群ごとに細胞の表現型に対応させる工程を含む方法1から8のいずれか1項に記載の方法。
(方法10)
前記細胞の表現型は、がん細胞と非がん細胞を含む組合せ、未分化細胞と分化細胞を含む組合せ、分化細胞と脱分化細胞を含む組合せ、及び同種がん細胞の上皮系と間葉系細胞を含む組合せ、から選ばれる少なくとも1の組合せを含む方法9に記載の方法。
(方法11)
前記複数の細胞に関する自家蛍光画像は、60mm2以上の領域を撮像して取得されることを特徴とする方法1から10のいずれかに記載の方法。
(方法12)
前記複数の細胞に関する自家蛍光画像は、少なくとも30000画素以上の画素数と、
かつ1画素あたり20μm以下の解像度を有することを特徴とする方法1から11のいずれかに記載の方法。
(方法13)
さらに、前記分類する工程によって分けられた群に関する定量的な情報を取得する工程を含む、方法1乃至12のいずれかに記載の方法。
(構成1)
画像取得装置と画像処理装置を含む解析装置であって、
前記画像取得装置は、
複数の細胞に関する自家蛍光画像を取得し、
前記画像処理装置は、
前記自家蛍光画像を用いて前記複数の細胞に関する自家蛍光情報を取得し、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得し、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、個々の細胞の自家蛍光情報を抽出し、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類することを特徴とする解析装置。
(構成2)
前記画像処理装置は、前記領域情報を前記複数の細胞に関する自家蛍光画像を用いて取得することを特徴とする構成1に記載の解析装置。
(構成3)
さらに、前記画像取得装置は、
前記複数の細胞に関する領域抽出用画像を取得し、
前記画像処理装置は、
前記領域情報を前記複数の細胞に関する領域抽出用画像を用いて取得することを特徴とする構成1に記載の解析装置。
(構成4)
前記画像処理装置が、前記分類をする際、少なくとも前記個々の細胞の自家蛍光の輝度を1つの軸とする、1次元、2次元あるいは3次元のグラフとして可視化することを特徴とする構成1から3のいずれか1項に記載の解析装置。
(構成5)
前記画像処理装置が、前記分類する際、少なくとも前記個々の細胞の自家蛍光情報を1つの特徴量として、主成分分析、k-means法、GMMからなる群より選ばれるいずれかの教師なし機械学習を行うことを特徴とする構成1から4のいずれかに記載の解析装置。
(構成6)
さらに、前記画像処理装置は、前記2以上の群を、群ごとに細胞の表現型に対応させることを特徴とする構成1から5のいずれかに記載の解析装置。
(構成7)
前記細胞の表現型は、がん細胞と非がん細胞を含む組合せ、未分化細胞と分化細胞を含む組合せ、分化細胞と脱分化細胞を含む組合せ、及び同種がん細胞の上皮系と間葉系細胞を含む組合せ、から選ばれる少なくとも1の組合せを含む構成6に記載の解析装置。
本発明は上記実施の形態に制限されるものではなく、本発明の精神及び範囲から離脱することなく、様々な変更及び変形が可能である。従って、本発明の範囲を公にするために以下の請求項を添付する。
100 解析装置
1A 画像取得装置
2A 画像処理装置
3A ケーブル
11 照射部
12 結像部
13 撮像部
14 制御部
15 通信I/F
16 操作部
17 バス
2A 画像処理装置
21 制御部
22 操作部
23 表示部
24 通信I/F
25 記憶部
26 バス
1A 画像取得装置
2A 画像処理装置
3A ケーブル
11 照射部
12 結像部
13 撮像部
14 制御部
15 通信I/F
16 操作部
17 バス
2A 画像処理装置
21 制御部
22 操作部
23 表示部
24 通信I/F
25 記憶部
26 バス
Claims (20)
- 複数の細胞に関する自家蛍光画像を用いて、前記複数の細胞に関する自家蛍光情報を取得する工程と、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得する工程と、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、前記個々の細胞の自家蛍光情報を抽出する工程と、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類する工程、を含む方法。 - 前記領域情報は前記複数の細胞に関する自家蛍光画像を用いて取得することを特徴とする請求項1に記載の方法。
- さらに、前記複数の細胞に関する領域抽出用画像を取得する工程を含み、
前記領域情報は、前記複数の細胞に関する領域抽出用画像を用いて取得することを特徴とする請求項1に記載の方法。 - 前記自家蛍光情報は、自家蛍光の輝度を含む請求項1に記載の方法。
- 前記輝度は、RGBの色空間におけるR、G、及びBの少なくともいずれかの指標として算出される値である請求項4に記載の方法。
- 前記輝度は、Lab色空間またはHSV色空間における指標として算出される値である請求項4に記載の方法。
- 前記分類する工程では、少なくとも前記個々の細胞の自家蛍光の輝度を1つの軸とする、1次元、2次元あるいは3次元のグラフとして可視化することを特徴とする請求項1から6のいずれか1項に記載の方法。
- 前記分類する工程では、少なくとも前記個々の細胞の自家蛍光情報を1つの特徴量として、主成分分析、k-means法、GMMからなる群より選ばれるいずれかの教師なし機械学習を行うことを特徴とする請求項1から6のいずれか1項に記載の方法。
- さらに、前記2以上の群を、群ごとに細胞の表現型に対応させる工程を含む請求項1から6のいずれか1項に記載の方法。
- 前記細胞の表現型は、がん細胞と非がん細胞を含む組合せ、未分化細胞と分化細胞を含む組合せ、分化細胞と脱分化細胞を含む組合せ、及び同種がん細胞の上皮系と間葉系細胞を含む組合せ、から選ばれる少なくとも1の組合せを含む請求項9に記載の方法。
- 前記複数の細胞に関する自家蛍光画像は、60mm2以上の領域を撮像して取得されることを特徴とする請求項1から6のいずれか1項に記載の方法。
- 前記複数の細胞に関する自家蛍光画像は、少なくとも30000画素以上の画素数と、
かつ1画素あたり20μm以下の解像度を有することを特徴とする請求項1から6のいずれか1項に記載の方法。 - さらに、前記分類する工程によって分けられた群に関する定量的な情報を取得する工程を含む、請求項1乃至6のいずれか1項に記載の方法。
- 画像取得装置と画像処理装置を含む解析装置であって、
前記画像取得装置は、
複数の細胞に関する自家蛍光画像を取得し、
前記画像処理装置は、
前記自家蛍光画像を用いて前記複数の細胞に関する自家蛍光情報を取得し、
前記自家蛍光画像における個々の細胞の領域を示す領域情報を取得し、
前記複数の細胞に関する自家蛍光情報と前記領域情報から、個々の細胞の自家蛍光情報を抽出し、
前記個々の細胞の自家蛍光情報を用いてクラスタリングを行うことで、前記複数の細胞を2以上の群に分類することを特徴とする解析装置。 - 前記画像処理装置は、前記領域情報を前記複数の細胞に関する自家蛍光画像を用いて取得することを特徴とする請求項14に記載の解析装置。
- さらに、前記画像取得装置は、
前記複数の細胞に関する領域抽出用画像を取得し、
前記画像処理装置は、
前記領域情報を前記複数の細胞に関する領域抽出用画像を用いて取得することを特徴とする請求項14に記載の解析装置。 - 前記画像処理装置が、前記分類をする際、少なくとも前記個々の細胞の自家蛍光の輝度を1つの軸とする、1次元、2次元あるいは3次元のグラフとして可視化することを特徴とする請求項14から16のいずれか1項に記載の解析装置。
- 前記画像処理装置が、前記分類する際、少なくとも前記個々の細胞の自家蛍光情報を1つの特徴量として、主成分分析、k-means法、GMMからなる群より選ばれるいずれかの教師なし機械学習を行うことを特徴とする請求項14から16のいずれか1項に記載の解析装置。
- さらに、前記画像処理装置は、前記2以上の群を、群ごとに細胞の表現型に対応させることを特徴とする請求項14から16のいずれか1項に記載の解析装置。
- 前記細胞の表現型は、がん細胞と非がん細胞を含む組合せ、未分化細胞と分化細胞を含む組合せ、分化細胞と脱分化細胞を含む組合せ、及び同種がん細胞の上皮系と間葉系細胞を含む組合せ、から選ばれる少なくとも1の組合せを含む請求項19に記載の解析装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021126090 | 2021-07-30 | ||
JP2021126090 | 2021-07-30 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2023021090A true JP2023021090A (ja) | 2023-02-09 |
Family
ID=85086953
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2022121735A Pending JP2023021090A (ja) | 2021-07-30 | 2022-07-29 | 細胞自家蛍光の画像を用いて細胞を分類する方法及び解析装置 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP2023021090A (ja) |
WO (1) | WO2023008369A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116453116B (zh) * | 2023-06-20 | 2023-09-08 | 杭州济扶科技有限公司 | 一种高适应性的明场细胞活死分类算法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1821097A1 (en) * | 2002-12-02 | 2007-08-22 | River Diagnostics B.V. | Use of high wavenumber Raman spectroscopy for measuring tissue |
JP5804194B2 (ja) | 2012-03-30 | 2015-11-04 | コニカミノルタ株式会社 | 医用画像処理装置及びプログラム |
CN109073536B (zh) * | 2016-05-06 | 2021-08-06 | 索尼公司 | 信息处理装置、信息处理方法、程序以及信息处理系统 |
US11002678B2 (en) | 2016-12-22 | 2021-05-11 | University Of Tsukuba | Data creation method and data use method |
JP6763589B1 (ja) | 2020-02-15 | 2020-09-30 | 東京港醸造株式会社 | 製麹装置及びそれを用いた製麹方法 |
-
2022
- 2022-07-25 WO PCT/JP2022/028624 patent/WO2023008369A1/ja active Application Filing
- 2022-07-29 JP JP2022121735A patent/JP2023021090A/ja active Pending
Also Published As
Publication number | Publication date |
---|---|
WO2023008369A1 (ja) | 2023-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10025271B2 (en) | Method and system for detecting and/or classifying cancerous cells in a cell sample | |
JP2024045407A (ja) | 身体試料中の実体を検出する方法および装置 | |
US8320655B2 (en) | Process and system for analyzing the expression of biomarkers in cells | |
US10083340B2 (en) | Automated cell segmentation quality control | |
Holzner et al. | High-throughput multiparametric imaging flow cytometry: toward diffraction-limited sub-cellular detection and monitoring of sub-cellular processes | |
JP6031513B2 (ja) | 網赤血球の同定および測定 | |
CN111699510A (zh) | 数字病理学图像的变换 | |
US20070298454A1 (en) | Analysis Of Cell Morphology And Motility | |
WO2017082048A1 (ja) | 分類器構成方法およびこれを用いた細胞の生死判定方法 | |
WO2012060163A1 (ja) | 細胞分析装置 | |
JP2015510592A (ja) | デジタルホログラフィック顕微鏡を備えるフローサイトメータ | |
CN101149327A (zh) | 基于细胞显微图像信息的抗肿瘤药物评价和筛选方法 | |
EP2630492B1 (en) | Internal focus reference beads for imaging cytometry | |
JPWO2011132587A1 (ja) | 細胞観察装置および細胞観察方法 | |
WO2013132998A1 (ja) | 画像処理装置、顕微鏡システム、及び画像処理方法 | |
Krauß et al. | Colocalization of fluorescence and Raman microscopic images for the identification of subcellular compartments: a validation study | |
WO2023008369A1 (ja) | 細胞自家蛍光の画像を用いて細胞を分類する方法及び解析装置 | |
JP2013109119A (ja) | 顕微鏡制御装置およびプログラム | |
Takahama Junior et al. | Usefulness of tissue autofluorescence imaging in actinic cheilitis diagnosis | |
WO2021198279A1 (en) | Methods and devices for virtual scoring of tissue samples | |
Matsui et al. | Label‐free multiphoton excitation imaging as a promising diagnostic tool for breast cancer | |
KR20060022266A (ko) | 산점도 분포를 사용한 슬라이드를 분류하기 위한 시스템 | |
JP2022550916A (ja) | 深紫外線顕微鏡を用いた無標識の血液学的および組織病理学的分析 | |
US20240029458A1 (en) | A method for automated determination of platelet count based on microscopic images of peripheral blood smears | |
Pop et al. | A fast and automated framework for extraction of nuclei from cluttered 3D images in fluorescence microscopy |