(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-09-07
(45)【発行日】2022-09-15
(54)【発明の名称】ブリルアン光散乱に基づくラベルフリーサイトメトリーのシステム及び方法
(51)【国際特許分類】
G01N 21/17 20060101AFI20220908BHJP
C12Q 1/04 20060101ALI20220908BHJP
C12M 1/34 20060101ALI20220908BHJP
G01N 21/64 20060101ALI20220908BHJP
G01N 33/48 20060101ALI20220908BHJP
G01N 21/49 20060101ALI20220908BHJP
【FI】
G01N21/17 A
C12Q1/04
C12M1/34 B
G01N21/64 Z
G01N33/48 M
G01N21/49 Z
(21)【出願番号】P 2018553041
(86)(22)【出願日】2016-12-22
(86)【国際出願番号】 US2016068405
(87)【国際公開番号】W WO2017112896
(87)【国際公開日】2017-06-29
【審査請求日】2019-12-16
(32)【優先日】2015-12-22
(33)【優先権主張国・地域又は機関】US
(32)【優先日】2016-04-15
(33)【優先権主張国・地域又は機関】US
(32)【優先日】2016-05-20
(33)【優先権主張国・地域又は機関】US
(32)【優先日】2016-11-21
(33)【優先権主張国・地域又は機関】US
(73)【特許権者】
【識別番号】507028217
【氏名又は名称】キヤノン ユー.エス. ライフ サイエンシズ, インコーポレイテッド
【氏名又は名称原語表記】CANON U.S. LIFE SCIENCES, INC.
(73)【特許権者】
【識別番号】518223823
【氏名又は名称】ユニバーシティ オブ メリーランド,カレッジ パーク
(74)【代理人】
【識別番号】100094112
【氏名又は名称】岡部 讓
(74)【代理人】
【識別番号】100128668
【氏名又は名称】齋藤 正巳
(74)【代理人】
【識別番号】100096943
【氏名又は名称】臼井 伸一
(72)【発明者】
【氏名】スカセリ,グリアノ
(72)【発明者】
【氏名】ジャンクー,ジタオ
(72)【発明者】
【氏名】フィオリ,アントニオ
(72)【発明者】
【氏名】キム,ハニオプ
【審査官】赤木 貴則
(56)【参考文献】
【文献】特開平07-072066(JP,A)
【文献】米国特許出願公開第2009/0323056(US,A1)
【文献】米国特許出願公開第2012/0302862(US,A1)
【文献】米国特許出願公開第2009/0004670(US,A1)
【文献】LINDSAY, S. M. et al.,Correction of Brillouin linewidths measured by multipass Fabry-Perot spectroscopy,Applied Optics,1977年,Vo. 16, No. 5,pp. 1404-1407
(58)【調査した分野】(Int.Cl.,DB名)
G01N 21/00-G01N 21/74
G01N 33/48-G01N 33/98
C12M 1/00-C12M 3/10
C12Q 1/00-C12Q 3/00
JSTPlus/JMEDPlus/JST7580(JDreamIII)
(57)【特許請求の範囲】
【請求項1】
サンプル内の複数の点においてブリルアン散乱スペクトルに関連付けられる1つ以上の指標を同時に取得する方法であって、
第1の方向に沿って照射光ビームによって前記サンプル
に照射光を照
射することと、
前記照
射光ビームに応答して前記サンプルから放射されるブリルアン散乱光を収集することと、
前記ブリルアン散乱光を、スペクトル分散を誘発する光学配置に送り、前記ブリルアン散乱光の空間スペクトルパターンを生成する検出ユニットに送ることであって、前記光学配置及び前記検出ユニットは第2の方向に沿って位置決めされることと、
前記検出ユニット上で前記ブリルアン散乱光の前記空間スペクトルパターンを検出することであって、前記照
射光ビームに沿った前記サンプルの複数の点が同時に測定されることと、
前記検出ユニットにおいて、各空間点における前記空間スペクトルパターンを較正することと、
前記検出された空間スペクトルパターンに基づいて、各測定サンプル点において前記1つ以上のブリルアン指標を計算することと、
前記サンプルの1以上の特性を同定するための前記サンプルの複数点から計算された1つ以上のブリルアン指標を使用することとを含む、方法。
【請求項2】
前記1つ以上のブリルアン指標は、ブリルアン周波数シフト、ブリルアンスペクトル線幅、ブリルアン利得又は損失スペクトル、及びその組み合わせからなる群から選択される、請求項1に記載の方法。
【請求項3】
前記方法は、各測定サンプル点において計算された前記1つ以上のブリルアン指標に基づいて、前記サンプルの画像を生成することを更に含み、結像中に前記サンプルが前記照
射光ビームに対して動いているか、又は結像中に前記照
射光ビームが静止した前記サンプルに対して動いている、請求項1又は2に記載の方法。
【請求項4】
前記第1の方向と前記第2の方向との間の角度は0度より大きい任意の角度である、請求項1~3のいずれか一項に記載の方法。
【請求項5】
スペクトル分散を誘発する前記光学配置は、仮想画像化フェーズドアレイ(VIPA)、ファブリ-ペローエタロン又はエシェル格子を含む、請求項1~4のいずれか一項に記載の方法。
【請求項6】
スペクトル分散を誘発する前記光学配置は、前記サンプルから前記検出ユニットへの光路内の前記空間スペクトルパターンのサイズ、形状及び/又は角度の広がりを変更する光学要素を更に備える、請求項5に記載の方法。
【請求項7】
前記サンプルは、生物有機体、組織又は生体細胞を含む、生体サンプルである、請求項1~6のいずれか一項に記載の方法。
【請求項8】
前記生体サンプルはマイクロ流体デバイスの1つ以上のチャネル内にある、請求項7に記載の方法。
【請求項9】
前記生体細胞は懸濁されるか、2D基体に付着しているか、又は3D細胞外マトリックス内で培養される、請求項7又は8に記載の方法。
【請求項10】
前記生体細胞は生細胞である、請求項7~9のいずれか一項に記載の方法。
【請求項11】
各測定サンプル点における前記1つ以上のブリルアン指標は、前記サンプルの物理的特性を示す、請求項1~10のいずれか一項に記載の方法。
【請求項12】
前記サンプルの物理的特性は、粘弾性率、密度、屈折率、電歪、及びその組み合わせからなる群から選択される、請求項11に記載の方法。
【請求項13】
前記サンプルの前記物理的特性は、前記サンプルから生成される前記ブリルアン散乱光の前記空間スペクトルパターンに基づいて取得される、請求項11又は12に記載の方法。
【請求項14】
前記サンプルを含む容器内にペンシルビームを生成するために第1のレンズによって前記照
射光ビームを再整形することを更に含む、請求項1~13のいずれか一項に記載の方法。
【請求項15】
第1のレンズに入射する前記照
射光ビームを再整形する空間光変調器又は変形可能ミラー、請求項1~14のいずれか一項に記載の方法。
【請求項16】
サンプル容器内に生成される前記ブリルアン散乱光を第2のレンズによって収集することと、第3のレンズの焦点面において中間像を生成することとを更に含む、請求項1~15のいずれか一項に記載の方法。
【請求項17】
前記中間像の倍率が、前記測定サンプル点の角分散を最小化するように最適化される、請求項16に記載の方法。
【請求項18】
中間像面内で空間フィルター又はアパーチャを用いて、前記サンプルから到来する焦点外の光を除去することを更に含む、請求項16又は17に記載の方法。
【請求項19】
前記方法は、第2の結像系を用いて、前記検出ユニット上に前記中間像を投影することを更に含み
、VIPAを含む、スペクトル分散を誘発する前記光学配置は、前記第2の結像系の無限遠空間内にある、請求項16~18のいずれか一項に記載の方法。
【請求項20】
前記方法は、球面レンズによって前記ブリルアン散乱光をコリメートすることと、前記コリメートされた光を円柱レンズによってz方向において前記VIPA上に合焦させることとを更に含み、異なる位置からのコリメートされたビームが前記VIPAに関するxy平面において異なる角度を有する、請求項5又は6に記載の方法。
【請求項21】
幾つかの円柱レンズ及び球面レンズによって変更される前記コリメートされた光から空間的に前記VIPAが離れており、前記空間スペクトルパターンの焦点を前記検出ユニット上に厳密に合わせることを更に含む、請求項20に記載の方法。
【請求項22】
前記方法は、第1のレンズアレイによって前記ブリルアン散乱光をコリメートすることを更に含み、前記第1のレンズアレイの各レンズレットが前記中間像の一部からの光のみを受光し、像全体が複数の部分に分割され、各部分がレンズレットによって独立してコリメートされる、請求項16に記載の方法。
【請求項23】
前記方法は、一対の円柱レンズによって、前記VIPA出力光を、第2のレンズアレイのレンズレットの口径に合うように圧縮することを更に含み、前記空間スペクトルパターンは、前記第2のレンズアレイの前方焦点面において生成され、前記検出ユニット上に結像される、請求項21に記載の方法。
【請求項24】
前記照
射光は、その中心波長に固定されるか、又はその中心波長付近において調整可能である、UV、可視光又はIR方式の単一波長を有する照
射源によって与えられる、請求項1~23のいずれか一項に記載の方法。
【請求項25】
前記検出ユニットの各空間点での前記空間スペクトルパターンを較正することが、前記検出ユニット上で、前記光学配置の異なる回折次数によって生成された異なるレーザーライン又は弾性散乱ライン間の距離を測定することを更に含む、請求項1~24のいずれか一項に記載の方法。
【請求項26】
前記検出ユニットの各空間点での前記空間スペクトルパターンを較正することが、既知のブリルアン特性の基準物質を用いて、スペクトル分散を誘発する前記光学配置のスペクトル分散特性を計算することを更に含む、請求項1~25のいずれか一項に記載の方法。
【請求項27】
前記照
射光ビームはレーザー源によって生成され、前記検出ユニットは、CCD、CMOSカメラ、又は検出器のアレイである、請求項1~26のいずれか一項に記載の方法。
【請求項28】
前記方法は、狭帯域フィルターを用いて前記空間スペクトルパターン内の前記照射光ビームを吸収することを更に含み、前記狭帯域フィルターは、吸収ガスセル及びファブリ-ペローエタロンデバイスからなる群から選択される、請求項1~27のいずれか一項に記載の方法。
【請求項29】
前記照
射光ビームは単一の波長を有し、及び前記狭帯域フィルターは単一の波長を吸収し、前記照射光ビームの前記波長と前記狭帯域フィルターによって吸収される前記波長は互いにロックされる、請求項28に記載の方法。
【請求項30】
サンプル内の複数の点においてブリルアン散乱スペクトルに関連付けられる1つ以上のブリルアン指標を同時に取得するシステムであって、
第1の方向に沿って照射光ビームによって前記サンプル
に照射光を照
射する照
射源と、
前記照
射光ビームに応答して、前記サンプルから放射されるブリルアン散乱光を収集する1つ以上のレンズと、
前記1つ以上のレンズから前記ブリルアン散乱光を受光する光学配置であって、該光学配置はスペクトル分散を誘発する、光学配置と、
前記ブリルアン散乱光の空間スペクトルパターンを検出する検出ユニットであって、前記1つ以上のレンズ、前記光学配置及び前記検出ユニットは、第2の方向に沿って位置決めされ、前記照
射光ビームに沿った前記サンプルの複数の点が同時に測定される、検出ユニットと、
プロセッサであって、
前記検出ユニットにおいて、各空間点における前記空間スペクトルパターンを較正するための命令と、
前記検出された空間スペクトルパターンに基づいて、各測定サンプル点において前記1つ以上のブリルアン指標を計算するための命令と、前記サンプルの複数点から計算された1つ以上のブリルアン指標に基づき前記サンプルの1以上の特性を同定するための命令と
を実行する、プロセッサと、
を備える、システム。
【請求項31】
前記1つ以上のブリルアン指標は、ブリルアン周波数シフト、ブリルアンスペクトル線幅、ブリルアン利得又は損失スペクトル、及びその組み合わせからなる群から選択される、請求項30に記載のシステム。
【請求項32】
各測定サンプル点において計算された前記1つ以上のブリルアン指標に基づいて、前記サンプルの画像が生成され、結像中に前記サンプルが前記照
射光ビームに対して動いているか、又は結像中に前記照
射光ビームが静止した前記サンプルに対して動いている、請求項30又は31に記載のシステム。
【請求項33】
前記第1の方向と前記第2の方向との間の角度は0度より大きい任意の角度である、請求項30~32のいずれか一項に記載のシステム。
【請求項34】
スペクトル分散を誘発する前記光学配置は、仮想画像化フェーズドアレイ(VIPA)、ファブリ-ペローエタロン又はエシェル格子を含む、請求項30~33のいずれか一項に記載のシステム。
【請求項35】
スペクトル分散を誘発する前記光学配置は、前記サンプルから前記検出ユニットへの光路内の前記空間スペクトルパターンのサイズ、形状及び/又は角度の広がりを変更する光学要素を更に含む、請求項34に記載のシステム。
【請求項36】
前記サンプルは、生物有機体、組織又は生体細胞を含む、生体サンプルである、請求項30~35のいずれか一項に記載のシステム。
【請求項37】
前記生体サンプルはマイクロ流体デバイスの1つ以上のチャネル内にある、請求項36に記載のシステム。
【請求項38】
前記生体細胞は懸濁されるか、2D基体に付着しているか、又は3D細胞外マトリックス内で培養される、請求項36又は37に記載のシステム。
【請求項39】
前記生体細胞は生細胞である、請求項36~38のいずれか一項に記載のシステム。
【請求項40】
各測定サンプル点における前記1つ以上のブリルアン指標は、前記サンプルの物理的特性を示す、請求項30~39のいずれか一項に記載のシステム。
【請求項41】
前記サンプルの物理的特性は、粘弾性率、密度、屈折率、電歪及びその組み合わせからなる群から選択される、請求項40に記載のシステム。
【請求項42】
前記サンプルの前記物理的特性は、前記サンプルから生成される前記ブリルアン散乱光の前記空間スペクトルパターンに基づいて取得される、請求項40又は41に記載のシステム。
【請求項43】
前記照
射光ビームは、前記サンプルを含む容器内にペンシルビームを生成するために第1のレンズによって再整形される、請求項30~42のいずれか一項に記載のシステム。
【請求項44】
前記光学配置は、第1のレンズに入射する前記照
射光ビームを再整形する空間光変調器又は変形可能ミラーを含む、請求項30~43のいずれか一項に記載のシステム。
【請求項45】
サンプル容器内に生成される前記ブリルアン散乱光は、第2のレンズによって収集され、第3のレンズの焦点面において中間像を生成する、請求項30~44のいずれか一項に記載のシステム。
【請求項46】
前記中間像の倍率が、前記測定サンプル点の角分散を最小化するように最適化される、請求項45に記載のシステム。
【請求項47】
中間像面内に、前記サンプルから到来する焦点外の光を除去する空間フィルター又はアパーチャを更に備える、請求項45又は46に記載のシステム。
【請求項48】
前記システムは、前記検出ユニット上に前記中間像を投影する第2の結像系を更に備え
、VIPAを含む、スペクトル分散を誘発する前記光学配置は、前記第2の結像系の無限遠空間内にある、請求項45~47のいずれか一項に記載のシステム。
【請求項49】
前記ブリルアン散乱光は球面レンズによってコリメートされ、前記コリメートされた光を円柱レンズによってz方向において前記VIPA上に合焦させ、異なる位置からのコリメートされたビームが前記VIPAに関するxy平面において異なる角度を有する、請求項34に記載のシステム。
【請求項50】
前記空間スペクトルパターンの焦点を前記検出ユニット上に厳密に合わせるために、幾つかの円柱レンズ及び球面レンズにより変更される前記コリメートされた光から空間的に前記VIPAが離れている、請求項49に記載のシステム。
【請求項51】
前記ブリルアン散乱光は第1のレンズアレイによってコリメートされ、前記第1のレンズアレイの各レンズレットが前記中間像の一部からの光のみを受光し、像全体が複数の部分に分割され、各部分がレンズレットによって独立してコリメートされる、請求項45に記載のシステム。
【請求項52】
前記VIPA出力光は、一対の円柱レンズによって、第2のレンズアレイのレンズレットの口径に合うように圧縮され、前記空間スペクトルパターンは、前記第2のレンズアレイの前方焦点面において生成され、前記検出ユニット上に結像される、請求項50に記載のシステム。
【請求項53】
前記照
射光は、その中心波長に固定されるか、又はその中心波長付近において調整可能である、UV、可視光又はIR方式の単一波長を有する照
射源によって与えられる、請求項30~52のいずれか一項に記載のシステム。
【請求項54】
前記検出ユニットの各空間点での前記空間スペクトルパターンを較正することが、前記検出ユニット上で、前記光学配置の異なる回折次数によって生成された異なるレーザーライン又は弾性散乱ライン間の距離を測定することを更に含む、請求項30~53のいずれか一項に記載のシステム。
【請求項55】
前記検出ユニットの各空間点での前記空間スペクトルパターンを較正することが、既知のブリルアン特性の基準物質を用いて、スペクトル分散を誘発する前記光学配置のスペクトル分散特性を計算することを更に含む、請求項30~54のいずれか一項に記載のシステム。
【請求項56】
前記照
射光ビームはレーザー源によって生成され、前記検出ユニットは、CCD、CMOSカメラ、又は検出器のアレイである、請求項30~55のいずれか一項に記載のシステム。
【請求項57】
前記システムは、前記空間スペクトルパターン内の前記照
射光ビームを吸収する狭帯域フィルターを更に備え、前記狭帯域フィルターは、吸収ガスセル及びファブリ-ペローエタロンデバイスからなる群から選択される、請求項30~56のいずれか一項に記載のシステム。
【請求項58】
前記照
射光ビームは単一の波長を有し、及び前記狭帯域フィルターは単一の波長を吸収し、前記照射光ビームの前記波長と前記狭帯域フィルターによって吸収される前記波長は互いにロックされる、請求項57に記載のシステム。
【発明の詳細な説明】
【技術分野】
【0001】
[関連出願の相互参照]
本出願は、2015年12月22日に出願された米国仮特許出願第62/270,982号、2016年5月20日に出願された米国仮特許出願第62/339,512号、2016年4月15日に出願された米国仮特許出願第62/323,176号及び2016年11月21日に出願された米国仮特許出願第62/425,070号の恩典を主張し、それらの特許出願は、引用することによりその全体が本明細書の一部をなす。
【0002】
本発明は、ブリルアン光散乱技法に基づいてラベルフリー細胞分析を実行する方法及びシステムに関する。具体的には、本発明による方法及びシステムは、照明方向に沿って複数のサンプル及び/又は一サンプル内の複数の点を同時に分析できるようにするブリルアン分光法に基づく。本発明によるシステム及び方法をマイクロ流体システムとともに利用して、分析を妨害しないように、接触することなく細胞の物理的特性を分析することができる。さらに、細胞のブリルアンスペクトルを分析することによって、細胞内(subcellular:亜細胞)物理的特性を取得することができる。
【背景技術】
【0003】
ブリルアン散乱は、物質の音響フォノンによって誘発される非弾性光散乱現象である。弾性的に散乱した光から小さな(通常、GHzオーダーの)ブリルアン周波数シフトを分離するために、従来のブリルアン分光法では一般に、マルチパス走査型ファブリ-ペロー干渉計のような高分解能分光計が使用される(非特許文献1)。音響フォノンの動態は物質の粘弾性特性に直接結び付けられるので、散乱光のブリルアン周波数シフトを測定することによって、機械的情報を取得することができる(非特許文献2)。しかしながら、この方法は、スペクトルを逐点走査することに起因して、かなりの時間を要する。さらに、スループット効率はエタロンのフィネスに制限される。このボトルネックは、大きな角分散を生成することができ(非特許文献3)、高いスループットを有するCCDカメラによって全てのスペクトル成分を同時に測定できるようにする仮想画像化フェーズドアレイ(VIPA:virtually imaged phased array)を使用することによって克服された。このタイプの分光計を使用して、生物組織内のレーザー走査共焦点ブリルアン顕微鏡法が低い照明電力及び短い積分時間において実行され、その方法では、サンプル内の各点が順次に照明される。ブリルアンに基づく弾性マップを作成するために、ブリルアンスペクトルが分析された(非特許文献4及び非特許文献5)。
【0004】
ブリルアン分光法は、サンプル内の音響的に誘発された光散乱の周波数スペクトルを測定することによって、機械的特性を非侵襲的に測定することが可能なものである。サブミクロン分解能におけるブリルアン顕微鏡法は、細胞の物理化学的特性の測定を可能にした(非特許文献6)。実際に、細胞の機械的特性は、数多くの細胞機能、例えば、増殖、遊走、遺伝子発現だけでなく、組織レベルの挙動、例えば、組織形態形成、転移及び血管形成も決定的に制御する。結果として、「機械的フェノタイピング(mechanical phenotyping)」、すなわち、その機械的特性に基づいて細胞を分類する能力が、細胞状態及び生理学的/病理学的条件を特徴付ける強力な手法として現れた。例えば、細胞剛性(cell stiffness)の減少が、転移能の向上と相関することがわかっており、それゆえ、腫瘍の検出及び病期分類のための新規のラベルフリーマーカーとして提案されてきた(非特許文献7及び非特許文献8)。
【0005】
ブリルアン顕微鏡法の現在の構成は一般にエピ構成(epi-configuration)に基づいており、その構成では、照明用の同じ対物レンズによって、後方散乱光が収集される。この構成の1つの短所は、強い後方反射が背景雑音として分光計の中に結合され、それにより、高い消光比を有する分光計(一般には2ステージVIPA)が必要とされることである(非特許文献9)。さらに、各時点において対物レンズの焦点面にある1つの点しか測定することができないので、結像のために逐点走査が必要とされる。
【0006】
細胞機構に関して、高い分解能、高密度マッピング及び高いスループットで、非接触及び非侵襲的に細胞内機構の探査を達成することが非常に重要である。実際には、細胞は非常に不均質であり、それゆえ、細胞内情報は、細胞の生物活性を完全に理解するのに極めて重要である。例えば、細胞核及び細胞骨格の異なる挙動は、細胞遊走、増殖及び転移進行中の癌性細胞遊走等の幾つかのプロセスに非常に関連性がある。そのサイズ及び高い剛性(rigidity)に起因して、細胞核は、細胞遊走のための主な物理的障壁を与える。一方、細胞骨格は、細胞質を満たし、それゆえ、細胞の形状全体を保持する動的で、適応的な3D構造である。3D組織の中を遊走中に、細胞核の運動は細胞骨格の動態と協調されなければならず、このため、位置、剛性及び形状に関して複雑な変化を受け、それにより、遊走効率に影響を及ぼすことになる。正常細胞が癌性細胞に変化するとき、細胞骨格及び細胞核両方の変化が、細胞が付着し、移動する能力に影響を及ぼすことになる。それゆえ、細胞骨格及び細胞核の剛性を特徴付ける能力は、結果として、癌性細胞の転移能に関する有力な機械的バイオマーカーを生成することができる。
【0007】
したがって、ブリルアン顕微鏡法をマイクロ流体技術と組み合わせて、対比機構として細胞の固有の物理的特性を使用することによって、蛍光標識を用いることなく細胞を分類するサイトメトリープラットフォームを作り出す必要がある。具体的には、細胞内部の音響的特性に関連するブリルアン光散乱のスペクトル分析から、細胞内機械的情報を取得することが望ましい。また、細胞内から生成されたブリルアン散乱スペクトルを分析することによって取得された細胞内機械的特性に基づいて、生細胞を分類できることも望ましい。
【0008】
さらに、サンプル内の数百の点を、散乱スペクトルのライン走査並列検出を使用して同時に測定できるようにする新たなブリルアン顕微鏡法構成が必要とされている。言い換えると、1次元においてのみサンプルを走査することによって高解像度2次元ブリルアン画像を瞬時に取得するブリルアン分光法構成を提供することが極めて望ましい。
【先行技術文献】
【非特許文献】
【0009】
【文献】Lindsay S M、Burgess S及びShepherd I W「Correction of Brillouin linewidths measured by multipass Fabry-Perot spectroscopy」Appl. Opt. 16(5), 1404-1407 (1977)
【文献】Dil J G「Brillouin scattering in condensed matter」Rep. Prog. Phys. 45, 286-334 (1982)
【文献】Xiao S、Weiner A M及びLin C「A Dispersion Law for Virtually Imaged Phased- Array Spectral Dispersers Based on Paraxial Wave Theory」IEEE J. Quantum Electronics 40(4), 420-426 (2004)
【文献】Scarcelli G及びYun S H「Confocal Brillouin microscopy for three-dimensional mechanical imaging」Nature Photonics 2, 39-43 (2007)
【文献】Girard M J A、Dupps W J、Baskaran M、Scarcelli G、Yun S H、Quigley H A、Sigal I A及びStrouthidis N G「Translating Ocular Biomechanics into Clinical Practice: Current State and Future Prospects」Curr. Eye Res. 40(1), 1-18 (2015)
【文献】Scarcelli G、Polacheck WJ、Nia HT、Patel K、Grodzinsky AJ、Kamm RD、Yun SH「Noncontact three-dimensional mapping of intracellular hydromechanical properties by Brillouin microscopy」Nat Methods. 12, 1132-1134 (2015)
【文献】Swaminathan, Vら「Mechanical Stiffness Grades Metastatic Potential in Patient Tumor Cells and in Cancer Cell Lines」Cancer Research 71, 5075-5080 (2011)
【文献】Cross, S.E.、Jin, Y.S.、Rao, J. 及びGimzewski, J.K.「Nanomechanical analysis of cells from cancer patients」Nature Nanotechnology 2, 780-783 (2007)
【文献】Scarcelli, G.及びYun, S.H.「Multistage VIPA etalons for high-extinction parallel Brillouin spectroscopy」Optics Express 19, 10913-10922 (2011)
【発明の概要】
【0010】
本発明の一態様において、サンプル内の複数の点においてブリルアン散乱スペクトルに関連付けられる1つ以上のブリルアン指標を同時に取得する方法が提供される。具体的には、サンプルは、第1の方向に沿って光ビームによって照明される。照明光ビームに応答してサンプルから放射されるブリルアン散乱光が収集され、スペクトル分散を誘発する光学配置に送られ、その後、検出ユニットに送られる。光学配置及び検出ユニットは、第1の方向とは異なる第2の方向に沿って位置決めされる。照明光ビームに沿ったサンプルの複数の点が同時に測定されるように、ブリルアン散乱光の空間スペクトルパターンが検出ユニット上で検出される。次に、検出ユニットにおいて各空間点における空間スペクトルパターンが較正される。最終ステップにおいて、検出された空間スペクトルパターンに基づいて、各測定サンプル点において1つ以上のブリルアン指標が計算される。
【0011】
本発明の更に別の態様において、サンプル内の複数の点においてブリルアン散乱スペクトルに関連付けられる1つ以上のブリルアン指標を同時に取得するシステムが提供される。具体的には、サンプルは、光源によって与えられる光ビームによって第1の方向に沿って照明される。1つ以上のレンズが、照明光ビームに応答してサンプルから放射されるブリルアン散乱光を収集し、ブリルアン散乱光を、スペクトル分散を誘発する光学配置に導光し、その後、検出ユニットに導光する。光学配置、1つ以上のレンズ及び検出ユニットは、第1の方向とは異なる第2の方向に沿って位置決めされる。照明光ビームに沿ったサンプルの複数の点が同時に測定されるように、ブリルアン散乱光の空間スペクトルパターンが検出ユニット上で検出される。検出ユニットと通信するプロセッサが、VIPAにおける角分散を考慮に入れるように、各空間点における空間スペクトルパターンを較正する。さらに、プロセッサは、検出された空間スペクトルパターンに基づいて、各測定サンプル点において1つ以上のブリルアン指標を計算する。一実施の形態において、ブリルアン指標は、ブリルアン周波数シフト、ブリルアンスペクトル線幅、ブリルアン利得又は損失スペクトル、及びその組み合わせからなる群から選択される。更に別の実施の形態において、サンプルの物理的特性は、粘弾性率、密度、屈折率、電歪、及びその組み合わせからなる群から選択される。
【0012】
本発明の更に別の態様において、生体細胞を分類する方法及びシステムが提供される。そのシステムは、媒質内に生体細胞を含む生体サンプルを有する容器を含む。次に、生体サンプルが光源によって照明され、生体細胞及び媒質内からブリルアン散乱光が生成される。プロセッサが、各生体細胞内の複数の点においてブリルアン散乱スペクトルを計算する。測定されたブリルアン散乱スペクトルに基づいて、プロセッサによって、生体細胞内の異なる空間点における細胞内物理的特性に関連する1つ以上のブリルアン指標が計算される。さらに、プロセッサは、生体細胞内の異なる空間点における物理的特性に基づいて、生体細胞を分類する。一実施の形態において、細胞内物理的特性に関連する指標はブリルアン周波数シフトである。プロセッサは、生体細胞内の各測定点を含むブリルアン周波数シフトに関するヒストグラムを生成するように更に構成される。ヒストグラムを当てはめるために、プロセッサによって、ガウス分布の線形重ね合わせが適用される。次に、プロセッサはヒストグラム内の各ピークを特定し、ピークは、細胞内の異なる領域からの機械的シグネチャを表す。生体細胞内の異なる空間点における物理的特性は、特定された機械的シグネチャと相関がある。
【0013】
本明細書に組み込まれ、明細書の一部を形成する添付の図面は、本開示の主題に関する種々の実施形態を示す。図面において、同じ参照符号は、同一の、又は機能的に類似の要素を示す。
【図面の簡単な説明】
【0014】
【
図1】逐点走査モードのためのブリルアン分光法構成を示す図である。
【
図2】本発明の一実施形態による、多重化走査モードのためのブリルアン分光法構成を示す図である。
【
図3】本発明の別の実施形態による、多重化走査モードのためのブリルアン分光法構成を示す図である。
【
図4】
図2の構成を使用することによって取得されるスペクトルパターンの原画像、及び、スペクトルパターンの原画像に関連付けられる測定値に関する空間分解能の特性評価を例示する図であり、正方形が測定データを表し、青線が直線当てはめデータを表す。
【
図5】マイクロ流体チャネルを通って流れる2つの異なる液体に関連付けられるブリルアン周波数を例示する図である。
【
図6】本発明による多重化ブリルアン分光法と、従来の逐点ブリルアン分光法との間の信号対雑音比(SNR)の比較を示す図である。
【
図7】本発明の一実施形態による、多重化走査モードのためのブリルアン分光法構成を示す図である。
【
図8a】カメラによって取り込まれた元の分散パターンを例示する図である。
【
図8b】単一の点におけるスペクトルパターンを例示する図である。
【
図8c】SNR対照明エネルギーの対数プロットを例示する図である。
【
図8d】時間トレースデータの推定ブリルアン周波数の分布を例示する図である。
【
図9】
図7に従って取得された測定値に関する空間分解能の特性評価を例示する図であり、正方形は測定データを表し、赤線は直線当てはめデータを表す。
【
図10a】キュベット内のポリメチルメタクリレート(PMMA)レンズの写真を例示する図であり、その写真はキュベットを屈折率整合液で満たす前に撮影された。
【
図10b】カメラによって取得されたブリルアンスペクトルパターンのスナップショットを例示する図である。
【
図10c】屈折率整合液のブリルアンスペクトルを例示する図であり、点及び実曲線はそれぞれ、測定データ及び当てはめデータを表す。
【
図10d】PMMAレンズのブリルアンスペクトルを例示する図である。
【
図11a】
図7による構成によって走査されたPMMAレンズの断面を例示する図である。
【
図11b】
図10a中の点線によって示される場所におけるPMMAレンズの断面ブリルアン画像を例示する図である。
【
図11c】PMMAレンズの3D走査の図解を例示する図である。
【
図11d】PMMAレンズの3Dブリルアン結像を例示する図であり、各スライスは、z軸に沿った場所におけるPMMAレンズの断面画像を表す。
【
図12】本発明の一実施形態による、逐点走査モードのためのブリルアン分光法構成を示す図である。
【
図13】
図12によるブリルアン分光法構成に関連付けられる制御インターフェースのスナップショットを示す図である。
【
図14a】異なる時点におけるブリルアン信号及び明視野像の両方のスナップショットを例示する図である。
【
図14b】異なる時点におけるブリルアン信号及び明視野像の両方のスナップショットを例示する図である。
【
図14c】異なる時点におけるブリルアン信号及び明視野像の両方のスナップショットを例示する図である。
【
図14d】
図14aの未処理のブリルアン信号から計算されたブリルアン周波数の時間トレースである。
【
図14e】
図14bの未処理のブリルアン信号から計算されたブリルアン周波数の時間トレースである。
【
図14f】
図14cの未処理のブリルアン信号から計算されたブリルアン周波数の時間トレースである。
【
図15a】ブリルアン周波数に及ぼすシースフローの影響を示す図である。
【
図15b】ブリルアン周波数に及ぼすシースフローの影響を示す図である。
【
図16】本発明による、流動している細胞をレーザービームの合焦スポットに位置合わせする原理を示す図である。
【
図17a】
図12によるブリルアン分光法構成を使用することによって実施される細胞スペクトル分析を示す図である。
【
図17b】
図12によるブリルアン分光法構成を使用することによって実施される細胞スペクトル分析を示す図である。
【
図17c】
図12によるブリルアン分光法構成を使用することによって実施される細胞スペクトル分析を示す図である。
【
図18a】
図12によるブリルアン分光法構成を使用することによって取得された測定値から生じる細胞スペクトル分析を示す図である。
【
図18b】2Dブリルアン細胞画像から生じるスペクトル分析を示す図である。
【
図19d】同じ細胞のブリルアン画像を示す図である。
【
図19e】
図19cの画像に基づいて一緒にプロットされた細胞核及び細胞質両方のヒストグラムを示す図である。
【
図20a】細胞核に及ぼすサイトカラシンD(Cytochalasyn D)の影響を示す細胞スペクトル分析を示す図である。
【
図20b】細胞核に及ぼすサイトカラシンDの影響を示す細胞スペクトル分析を示す図である。
【
図20c】細胞核に及ぼすサイトカラシンDの影響を示す細胞スペクトル分析を示す図である。
【
図20d】細胞核に及ぼすサイトカラシンDの影響を示す細胞スペクトル分析を示す図である。
【
図20e】細胞核に及ぼすサイトカラシンDの影響を示す細胞スペクトル分析を示す図である。
【
図21】本発明の別の実施形態による、逐点走査モードのためのブリルアン分光法構成を示す図である。
【
図22】
図21のブリルアン分光法構成によって取得された細胞画像を示す図であり、画像内の各点は細胞内の或る特定の領域において取り込まれた1つの測定値を表し、異なる濃淡はデータ点の相対密度を符号化する。
【
図23a】
図21のブリルアン分光法構成によって取得された、細胞核に及ぼすサイトカラシンDの影響を示す図である。
【
図23b】
図21のブリルアン分光法構成によって取得された、細胞核に及ぼすサイトカラシンDの影響を示す図である。
【
図24】
図7のブリルアン分光法構成を使用することによって2D画像を取得する概念図である。
【
図25】ルビジウムガスを含むガスチャンバーの測定された吸収スペクトルを例示する図である。
【
図26】本発明の別の実施形態(例4において使用される)による、逐点走査モードのためのブリルアン分光法構成を示す図である。
【
図27a】異なる時点においてマイクロ流体チャネルを通って流れている細胞のスナップショットであり、スケールバーは30μmである。
【
図27b】異なる時点においてマイクロ流体チャネルを通って流れている細胞のスナップショットであり、スケールバーは30μmである。
【
図27c】異なる時点においてマイクロ流体チャネルを通って流れている細胞のスナップショットであり、スケールバーは30μmである。
【
図27d】対応するブリルアン信号の時間トレースを例示する図であり、点は測定データを表し、実線は見やすくするために引いた線である。
【
図28a】リン酸緩衝生理食塩水(PBS)及び細胞の両方のシグネチャを含む元のデータを例示する図である。
【
図28b】PBS自体のシグネチャを例示する図である。
【
図29a】細胞集団(29個の細胞)の2D画像を例示する図である。
【
図29b】2D画像からのデータ点の分布を示す図である。
【
図29c】細胞内の複数の点においてブリルアン周波数を取得することに基づく流動実験からのデータ点の分布を示す図である。
【
図30】細胞の蛍光画像によって検証される細胞のブリルアンシグネチャを例示する図であり、ヒストグラムはそれぞれ細胞核及び細胞質からのブリルアンシグネチャを表し、挿入写真は、(i)明視野画像、(ii)蛍光画像、(iii)(i)及び(ii)の融合画像、(iv)ブリルアン画像である。
【
図31a】クロマチン脱凝縮による細胞核の軟化の実験的観測を例示する図である。
【
図31b】クロマチン脱凝縮による細胞核の軟化の実験的観測を例示する図である。
【
図31c】クロマチン脱凝縮による細胞核の軟化の実験的観測を例示する図である。
【
図32】本発明の一実施形態によるシステムを示すブロック図である。
【発明を実施するための形態】
【0015】
本発明は幾つかの実施形態を有し、当該技術分野において既知の細部に関しては種々の特許、特許出願及び他の参考文献に頼る。それゆえ、本明細書において或る特許、特許出願又は他の参考文献が引用されるか、又は繰り返されるとき、その特許、特許出願又は他の参考文献は、あらゆる目的のために、そして説明される提案のために、引用することによりその全体が本明細書の一部をなすことは理解されたい。
【0016】
本発明は、ブリルアン分光法に基づくラベルフリー細胞分析の方法及びシステムに関する。細胞サンプルから放射されるブリルアン散乱光によって、サンプルを特徴付ける物理的情報を取得できるようになる。本発明の一態様では、これを果たすために、サンプルが第1の方向に沿って光ビームによって照明され、サンプルによって放射されたブリルアン散乱光が第2の方向に沿って収集され、検出される。第1の方向に沿ったサンプルの複数の位置を同時に測定することができ、それにより、測定スループットを実効的に改善することができる。さらに、ブリルアンスペクトルから取得された細胞内情報に基づいて、サンプル内の生体細胞を分類することができる。
【0017】
等方性物質のブリルアン周波数シフトf
bは以下のように表すことができる。
【数1】
ただし、nは物質の屈折率であり、λはレーザーの波長であり、
【数2】
は音響速度であり、Eは弾性率であり、ρは密度であり、θは散乱角(ここでは、θ=90度)である。
【0018】
一サンプル内の複数の点におけるブリルアン散乱スペクトルに関連付けられる1つ以上の指標を測定することによって、測定サンプルの機械的及び/又は物理的特性を直接読み出すことができる。一実施形態において、指標としてブリルアン周波数シフトが使用される。サンプルの機械的及び/又は物理的特性は、粘弾性率、密度、屈折率及び電歪を含むことができる。測定サンプルの機械的及び/又は物理的特性を特定するために、ブリルアンスペクトル線幅、ブリルアン利得又は損失スペクトル及びその組み合わせ等の他の指標を使用することができる。
【0019】
図1は、マイクロ流体チャネル111とともにエピ検出を用いる逐点走査モードのためのブリルアン分光法構成を例示する。具体的には、到来するレーザービーム110が、最初に、ビームスプリッター113によって導光され、その後、対物レンズ112によってマイクロ流体チャネル111上に合焦する。励起されたブリルアン散乱光が、同じ対物レンズ112によって収集され、ビームスプリッター113及びミラー114によってコリメーター115に導光される。コリメーター115に後続するコリメートされたビームは、最初に、円柱レンズ116によって合焦し、その後、VIPA117に送り込まれる。到来する散乱光の異なるスペクトル成分がVIPA117によって空間分離され、球面レンズ118によって合焦し、スペクトルパターンが生成される。スペクトルパターンは、その後、球面レンズ119によって検出ユニット120上に結像される。一例として、検出ユニット120はカメラとすることができる。幾つかの実施形態において、球面レンズ119は不要な場合があり、その場合、カメラ120は、球面レンズ118の前方焦点面に直接配置することができる。マイクロ流体チャネル111の中にサンプルが流されるとき、そのブリルアンシグネチャを識別し、その固有の物理的特性に相関させることができる。したがって、各時点において、対物レンズ112の焦点面において1つの点しか測定することができないので、結像のために逐点走査が必要とされる。
【0020】
図2は、マイクロ流体チップとともに多重化ブリルアン分光法を用いる本発明の一実施形態に関連する。具体的には、到来する照明光ビーム210が対物レンズ211によって再整形され、サンプルを有するマイクロ流体チップ212内にペンシルビームが生成される。検出ユニット227に通じる検出経路(x方向に沿っている)は、照明経路(y方向に沿っている)に直交する。一例として、検出ユニット227はカメラとすることができる。チップ212内に生成される散乱光は、最初に、第2の対物レンズ220によって収集され、チューブレンズ221の焦点面において中間像が生成される。一実施形態において、中間像面内に空間フィルター又はアパーチャ228が位置決めされ、サンプルから到来する焦点外の光が除去される。中間像がその後方焦点面上にあるように球面レンズ222が配置され、それゆれ、光が球面レンズ222によってコリメートされる。コリメートされた光は、その後、円柱レンズ223によって、z方向においてVIPA224上に合焦する。到来する散乱光の異なるスペクトル成分がVIPA224によって空間分離され、球面レンズ225によって合焦し、スペクトルパターンが生成される。スペクトルパターンは、その後、球面レンズ226によってカメラ227上に結像される。検出経路において、ビームラインの3つのグループが、サンプルの異なる位置を示す。対物レンズ220の視野内で、サンプルの複数の位置を同時に測定できることが、
図2から明らかである。幾つかの実施形態において、カメラ227が球面レンズ225の前方焦点面に直接配置されるように、球面レンズ226を除くことができる。
【0021】
照明ビームに沿った全ての測定位置から収集された散乱光が同じ球面レンズ222によってコリメートされるので、異なる位置からのコリメートされたビームは、VIPA224に関するxy平面において異なる角度を有する。この結果として、軸上の点に対して、軸外の点に関するスペクトルに更なる空間的なシフトが生じることになり、それゆえ、空間次元(すなわち、y方向)に沿ったスペクトルは、カメラ227において一直線ではなく、湾曲する。したがって、照明ラインを横切る異なる点からの散乱光の入射角は異なる。これは、サンプル内の異なる点からの散乱がVIPAエタロン内でわずかに異なる経路長及び分散を有することを意味する。角分散を、それゆえ、スペクトル線の曲がりを除去するために、レンズ221の中間結像面において生じる結像倍率が最小化された。VIPA224においてスペクトル分散が生じた後に、カメラ227上のスペクトル線の湾曲に影響を及ぼすことなく、結像倍率を復元することができる。
【0022】
さらに、角分散を補償するために、スペクトル較正が必要とされる場合がある。カメラ227上で検出されたスペクトル内に照明ビームラインが存在するとき、VIPAエタロンの種々の異なる回折次数におけるシフトされない照明ビームラインの場所を用いて、空間ラインの各ピクセル(サンプルの1つの点に対応する)においてスペクトル較正を実行することができる。照明ビームラインが利用できないとき、各空間点における較正のために必要とされる2つの未知のパラメーター、すなわち、自由スペクトル領域(FSR)及びスペクトル分散係数がある。それゆえ、未知のパラメーターを求めるのに、2つ以上の既知のサンプルの結果を組み合わせれば十分である。
【0023】
代替的には、
図3において例示されるような多重化ブリルアン分光法に基づく本発明の別の実施形態を使用することによって、空間較正を回避することができる。
【0024】
図3は、本発明の別の実施形態による、レンズアレイを使用する多重化ブリルアン分光法構成を示す。照明経路は
図2に示される照明経路と同じである。対物レンズ311を用いて再整形された光310のビームによって、マイクロ流体チップ312が照明される。検出経路(x方向に沿っている)は、照明経路(y方向に沿っている)の法線方向にある。チップ312内に生成される散乱光は、最初に、第2の対物レンズ320によって収集される。チューブレンズ321の焦点面において中間像が生成される。一実施形態において、中間像面内に空間フィルター又はアパーチャ330が位置決めされ、サンプルから到来する焦点外の光が除去される。
図2の球面レンズ222が、コリメーションのためのレンズアレイ322に置き換えられる。像全体が数多くの部分に分割され、それぞれがレンズレットによって独立してコリメートされるように、レンズアレイ322の各レンズレットが、中間像の小さな部分からの光のみを受光する。各小部分がレンズレットに対して軸上にあると近似され、全てのレンズレットが十分に位置合わせされるので、コリメートされたビーム全体が平行であり、対応するスペクトルは一直線になるであろう。コリメートされた光は、その後、VIPA224の中に送り込まれる。第2のレンズアレイ327のレンズレットの口径を適合させるために、一対の円柱レンズ325及び326を用いて、VIPA224の出力光を圧縮する。レンズアレイ327の前方焦点面においてスペクトルパターンが生成され、球面レンズ328によって検出ユニット329上に結像される。一例として、検出ユニット329はカメラとすることができる。幾つかの実施形態において、レンズアレイ327の前方焦点面にカメラ329を直接配置できるように、球面レンズ328を除くことができる。
【0025】
例1
本発明の検証の一例として、
図2の構成を用いて多重化ブリルアン分光法の空間分解能を実証するために実験が実施された。この実験において、光源として、単一モード532nm連続波レーザーが使用され、サンプルとして、メタノールを含むプラスチックキュベットが使用された。照明経路(y軸に沿っている)内の対物レンズ211は、0.0175の開口数(NA)を有し、検出経路(x軸に沿っている)内の対物レンズ220は0.1のNAを有する。VIPA224は、17GHzの自由スペクトル領域(FSR)と、20mmの入射窓とを有する。カメラ227は、電子増倍電荷結合素子(EMCCD:electron multiplying coupled charge device)である。サンプル(プラスチックキュベット)と対物レンズ220との間にナイフエッジが配置された。ナイフエッジは、並進ステージを用いて、y方向に沿って動かすことができる。ナイフエッジをy方向に向かって動かすとき、キュベットからの散乱信号が部分的に遮断される。ナイフエッジの変位と、カメラの記録された信号とを比較することによって、分光法の空間分解能を取得することができる。
図4a~
図4eは、ナイフエッジが異なる位置にあったときにカメラによって取り込まれた原画像を示す。これらの画像において、スペクトル領域内の複数のスペクトル次数が示された。各次数は3つの縞のグループを含み、両側に2つのブリルアン周波数(ストークス及び反ストークスシフト)と、その間にレーザー周波数とを有する。
図4fは、ナイフエッジの動きと測定された変位との間の線形関係を示す。カメラの1ピクセルはサンプルの2.45μmのサイズに対応することが示された。また、
図2に示されるような現在の構成は、800μm程度の大きいサンプルを測定することができ、結果として、300点より多くの点を同時に測定できることを示す。
【0026】
図2に示されるような本発明の実施形態を使用するマイクロ流体チャネル内の流動液体の測定が
図5において例示された。溶融石英ガラスから形成されるマイクロ流体チャネルは、100μmの幅と、250μmの深さとを有する。チャネル層は、500μmの幅を有する2つのガラス窓によって挟持される。チャネルは、光ビームが通過できるようにするためのその両側が研磨された1インチ×2インチのチップの中央に位置した。測定中に、液体がz方向に沿ってチャネル内で流動できるように、チップは直立するように配置された。レーザービーム210が、チップの側方からチャネルを照明する。水(5.2GHzのブリルアン周波数シフトを有する)及びメタノール(3.9GHzのブリルアン周波数シフトを有する)が、ハンドヘルドシリンジによって管を通して交互にチャネルの中に注入された。カメラの積分時間は、測定シーケンスの時間でもある0.2秒であった。
図5は、多重化ブリルアン分光法によって、マイクロ流体チャネル内の流動液体を明確に識別できることを示す。この実験は、フローサイトメトリーの適用例における本発明の能力を実証する。
【0027】
図2において例示される実施形態による、本発明の利点の別の証拠として、本発明(
図2)と従来のエピ検出構成(
図1)との間の信号対雑音比(SNR)が比較された。各構成のSNRは、光の異なる入力電力密度及びカメラの異なる露光時間において測定された。
図6は、本発明がショット雑音限界であり(理論的な当てはめ)、同じ入力エネルギー密度において従来のエピ検出よりはるかに高いSNRを有することを示す。言い換えると、本発明によれば、同じSNRレベルにおいて、従来のエピ検出より速く測定できるようになる。
【0028】
上記の方法以外に、測定スループットを更に改善する他の方法がある。一実施形態において、マイクロ流体チャネルの真下に励起源としてトランスデューサーを配置し、サンプルによって生成される音響信号を高めることができる。さらに、刺激されたブリルアン散乱は、一般に、自然発生のブリルアン散乱よりはるかに強い。単波長光源を使用する代わりに、可変の波長を有する2つのレーザーを用いて、マイクロ流体チャネル内に刺激されたブリルアン散乱を生成し、信号を更に高めることができる。さらに、1つの超短パルスレーザーを用いて、インパルス性の刺激されたブリルアン散乱を生成することもできる。
【0029】
例2
例2は、
図7に示されるような本発明の実施形態による、多重化(ライン走査)ブリルアン分光法の特性評価に関連する。
図7は、ミラー724と、対物レンズ711、713及び714と、ビームコリメーション用の球面レンズ715と、円柱レンズ716、718及び719と、結像用の球面レンズ722とを含む、ライン走査ブリルアン分光法構成を例示する。一実施形態において、中間像面内に空間フィルター又はアパーチャ725が位置決めされ、サンプルから到来する焦点外の光が除去される。
【0030】
図7による分光計のスペクトル分解能が、この例において特性評価された。1つの非限定的な実施形態において、光源710は、単一モードの532nm cw(連続波)レーザー(Torus、LaserQuantum社)であった。レーザーヘッドからの光は対物レンズ711(NA=0.0175)によって合焦し、サンプル712を照明するためのラインビームが生成された。90度において、ビームラインに沿った散乱光が最初に、一対の対物レンズ713及び714(いずれも4X/0.1NAであった)によって結像され、その後、球面レンズ715(f=400mm)によってコリメートされた。コリメートされた光は、その後、円柱レンズ716(f=200mm)によって、VIPA717(FSR=17GHz、フィネス=35)の入口において合焦した。VIPA717に後続する分散光は、一対の円柱レンズ718(f=1000mm)及び719(f=400mm)によってスリット720の平面上に結像された。この平面は、球面レンズ722(f=60mm)によってカメラ723(iXon、Andor社)上に再び結像された。スリット720の隙間を調整することによって、望ましくない周波数を遮断することができ、ブリルアン周波数成分のみを通過させることができた。
【0031】
サンプルとして、純メタノールを含むプラスチックキュベット712が使用される。カメラ723によって取得された元の分散パターンが
図8aに示される。上側及び下側の2つのグレーラインは弾性散乱周波数(レーザーの周波数と同じ)であり、その間にある2つの明るいラインはブリルアン周波数(ストークス及び反ストークス成分)である。横軸及び縦軸はそれぞれ、空間領域及びスペクトル領域に対応する。
【0032】
図8bは、
図8aにおいて点線によって示される場所におけるブリルアンスペクトルを示す。VIPAの自由スペクトル領域は17GHzであったので、スペクトル分解能がピクセルあたり0.177GHzであったと計算することができた。レーザー周波数の連続するピークと半値全幅との空間比を計算することによって、約35のフィネスも確認された。推定ブリルアン周波数シフトの精度は、空間領域においてピクセルごとの信号対雑音比(SNR)を計算することによって特性評価された。
図8cは、SNR対照明エネルギーの対数プロットを示す。
図8cのSNRは、異なる照明エネルギーで計算された。測定データは、概ね平方根依存性を示し、ショット雑音限界挙動を示した。
図8dは、単一の点において100回測定した場合の推定ブリルアン周波数の典型的な分布を示しており、その分布はガウス形状を有する。この分布の標準偏差を用いて、推定ブリルアン周波数の精度を推定した。この場合、それは8.5MHzであり、0.22%の相対的不確実性に相当する。
【0033】
空間分解能を特性評価するために、キュベット712に後続してすぐに、ただし検出対物レンズ713に先行して、ナイフエッジが配置される。一実施形態において、ナイフエッジを搬送する並進ステージが25.4μmのステップサイズでx方向に動かされ、ナイフエッジの対応する像が分光計のカメラ723によって監視された。その結果が
図9に示される。測定データが直線に当てはめられ、それにより、空間分解能はピクセルあたり3.29μmと計算された。カメラ723は全部で512×512ピクセルを有し、それはサンプル面内の概ね1.68mmに対応する。
図8aに示されるように、スペクトル線は空間次元にわたって完全に一直線ではなく、VIPA717における軸外の点の入射角の偏差に起因してわずかな湾曲を有する。この湾曲の結果として、カメラにわたってスペクトル分解能が変動する場合があるが、その変動はピクセルあたり0.2MHz未満まで、すなわち、1mmの空間的視野を超える範囲内で10MHz以下まで最小化された。幾つかの実施形態において、この誤差は、必要な場合に各点においてスペクトル分解能を較正することによって回避することができる。
【0034】
例3
2次元及び3次元結像
例3は、
図7による多重化ブリルアン分光法に基づく2次元及び3次元結像の特性評価に関連する。
【0035】
非球面PMMAレンズのブリルアンシフトが、
図10aに示される構成において測定された。
図10aは、屈折率整合液でキュベットを満たす前に撮影されたPMMAレンズの写真を示す。下はPMMAレンズを拡大撮影した写真であり、スケールバーは500μmの長さを有する。PMMAレンズはプラスチックキュベットの中に配置され、サンプルの表面における散乱を低減するために、キュベットは屈折率整合液であらかじめ満たされた。キュベットは、垂直に配置された電動並進ステージによって搬送され、その並進ステージによって、サンプルをy方向に走査できるようにした。一実施形態において、キュベット(
図7のサンプル712)は、10mm×10mmのサイズを有する。非球面ポリメチルメタクリレート(PMMA)レンズをキュベット712の中央領域に浮かせておくために、PMMAレンズの縁が光学接着剤を用いて注射器の針の先端に取り付けられ、その後、針の端部がキュベットの壁に固定された。一例として、限定はしないが、電動並進ステージ(T-LSM025A、Zaber社)は、25mm移動範囲とRS-232制御とを有した。並進ステージの最大速度は7mm/s程度に速くすることができる。
【0036】
一実施形態において、走査を自動的に実行できるように、Lab Viewプログラムを用いて、ステージの移動とカメラによる取得とを同期させる。レーザー電力は70mWであり、カメラの露光時間は0.1秒であった。ステージの速度は、30秒以内に300フレームが取り込まれるように50μm/secと設定され、それは全部で1.5mmの変位に対応した。サンプル(PMMAレンズ)が液体の中に浸漬されたとき、屈折率整合のために、裸眼ではほとんど見ることができなかった。しかしながら、サンプル及び整合液の剛性が異なったので、それらはブリルアンスペクトルによって容易に識別することができた。
図10bは、カメラ723によって取得された典型的な信号のスナップショットであり、スリット720によって弾性周波数が遮断されたために、ブリルアン周波数成分のみが示された。PMMAレンズが整合液によって包囲されたことが
図10bにおいて明確に示される。
図10c及び
図10dはそれぞれ、単一の点における屈折率整合液及びPMMAのブリルアンスペクトルを示す。点及び実線はそれぞれ、測定データ及び当てはめられたデータに対応する。スペクトル分解能の特性評価データを使用するとき、PMMA及び整合液のブリルアン周波数シフトはそれぞれ、11.32GHz及び9.11GHzと特定された。
【0037】
図11は、PMMAレンズの2Dブリルアン画像及び3Dブリルアン画像を示す。
図11aは、
図10a内の点線に沿って走査された断面を示す。矢印は、照明ビームを示す。
図11bは、測定されたブリルアン周波数シフトに基づく、整合液内に浸漬されたPMMAレンズの2D結像を示す。一実施形態において、画像の全体サイズは約1.1mm×1.5mmであった。PMMAレンズと整合液との界面は、画像から明確に確認することができ、PMMAレンズの内側領域は極めて均一である。異なる物質の界面において、各物質に対応する2つのブリルアンシグネチャ間にクロストークがある。2つの物質が類似のブリルアンシフトを有する場合には、このクロストーク効果が界面において不明瞭な領域を導入することになる。この実験において、PMMAレンズ及び屈折率整合液は異なるブリルアンシフトを有し、それにより、不明瞭な領域を、2ピクセル、すなわち、6.58μm以内にあると定量化できるようになる。
図11cに示される走査方法を使用することによって、この分光法が迅速な3D結像を実施できることが実証された。
図11dは、別の並進ステージを用いてPMMAレンズがz軸に沿って動かされるときに取得された断面画像の5つのスライスを示す。
【0038】
要するに、
図1~
図11を参照すると、本発明の一実施形態において、照明ビームはレーザービームとすることができ、検出ユニット120、227、329及び723は、CCD、CMOSカメラ、又は検出器のアレイとすることができる。更に別の実施形態において、照明光は、その中心値に固定されるか、又はその中心値付近において調整可能である、UV、可視光又はIR方式における単一波長を有する照明源によって与えられる。各測定サンプル点において計算された1つ以上のブリルアン指標に基づいてサンプルの画像を生成するために、サンプルは、結像又は取得中に照明光ビームに対して移動する可動プラットフォーム上に配置することができる。更に別の実施形態において、照明光ビームは、結像中に静止したサンプルに対して移動している場合がある。一例として、限定はしないが、サンプルは、生物有機体、組織又は生体細胞を含む、生体サンプルとすることができる。一実施形態において、生体細胞は生細胞である。生体細胞は、ブリルアン指標を測定している間に、懸濁される場合があるか、2D基体に付着している場合があるか、又は3D細胞外マトリックス内で培養される場合がある。
【0039】
一実施形態において、カメラ120、227、329及び723上で検出されたスペクトルは、ブリルアン散乱光の空間スペクトルパターンである。検出された空間スペクトルパターンに基づいて、各測定サンプル点において1つ以上のブリルアン指標を計算することができる。一例として、限定はしないが、ブリルアン指標はブリルアン周波数シフト、ブリルアンスペクトル線幅、ブリルアン利得又は損失スペクトル、及びその組み合わせとすることができる。サンプル点において測定されたブリルアン指標はそれぞれ、この時点におけるサンプルの物理的特性を示す。一例として、限定はしないが、サンプルの物理的特性は、粘弾性率、密度、屈折率、電歪及びその組み合わせとすることができる。一実施形態において、ブリルアン散乱光に関連付けられる指標は、ブリルアン周波数シフトである。
【0040】
各測定サンプル点における1つ以上のブリルアン指標に基づいて、サンプルの画像を生成することができる。
図2、
図3及び
図7は90度の角度をなすように配置される照明経路及び検出経路を示すが、いずれの構成も、互いに対して0度より大きい任意の角度に配置される照明経路及び検出経路を用いて実現することができる。
【0041】
図2、
図3及び
図7を参照すると、スペクトル分散を誘発する光学配置は、球面レンズ、円柱レンズ、VIPAを備えることができる。具体的には、その光学配置は、
図2の要素222、223、224、225及び226、
図3の要素322、323、324、325、326、327及び328、並びに
図7の要素715、716、717、719、720及び722を備える。一実施形態において、スペクトル分散を誘発する光学配置は、仮想画像化フェーズドアレイ(VIPA)、ファブリ-ペローエタロン又はエシェル格子を備える。更に別の実施形態において、光学配置は、サンプルから検出ユニットまでの光路内の空間スペクトルパターンのサイズ、形状及び/又は角度の広がりを変更する光学要素を更に備える。具体的には、照明光ビームは、サンプルを含む容器内にペンシルビームを生成するために、第1のレンズ211、311及び711によって再整形される。空間光変調器又は変形可能ミラー(
図7、ミラー724)を用いて、第1のレンズ(
図7、レンズ711)に入射する照明光ビームを再整形することができる。更に別の実施形態において、サンプル容器内で生成されるブリルアン散乱光は、第2のレンズ220、320、713を備え、第3のレンズ221、321、715の焦点面において中間像を生成する第1の結像系によって収集される。中間像の倍率は、VIPA224、324及び717における測定サンプル点の角分散を最小化するように最適化される。サンプルから到来する焦点外の光を除去するために、中間像面において空間フィルター又はアパーチャ228、330及び725を使用することができる。
【0042】
図2の要素222、225、226、
図3の要素322、326、327及び328、
図7の要素715、718及び722の組み合わせである第2の結像系が中間像を検出ユニット227、329、723上に投影し、VIPAを含む光学配置は、第2の結像系の無限遠空間(infinity space)内にある。本発明の更なる態様において、ブリルアン散乱光は、球面レンズ(
図2、レンズ222;
図7、レンズ715)によってコリメートされ、円柱レンズ(
図2、レンズ223;
図7、レンズ716)によってz方向においてVIPA上に合焦し、異なる位置からのコリメートされたビームは、VIPAに関するxy平面において異なる角度を有する。ブリルアン空間スペクトルパターンの焦点を検出ユニット227、329、723上に厳密に合わせるために、VIPAの出力光は幾つかの円柱レンズ及び球面レンズ225、226、718、719、722によって変更される。さらに、本発明の別の態様において、ブリルアン散乱光は第1のレンズアレイ322によってコリメートされ、第1のレンズアレイの各レンズレットが中間像の一部からの光のみを受光し、像全体は複数の部分に分割され、各部分はレンズレットによって独立してコリメートされる。VIPAの出力光は、第2のレンズアレイ327のレンズレットの口径を適合させるために、一対の円柱レンズ325及び326によって補償される。ブリルアン空間スペクトルパターンは、第2のレンズアレイ327の前方焦点面において生成され、検出ユニット329上に結像される。
【0043】
更に別の実施形態において、空間スペクトルパターン内のレーザーラインを吸収する狭帯域フィルターであり、狭帯域フィルターは吸収ガスセル及びファブリ-ペローエタロンデバイスからなる群から選択される。照明源の波長及び狭帯域フィルターによって吸収される波長は互いにロック(lock)される。
【0044】
一実施形態において、ブリルアン空間スペクトルパターンは、光学配置の異なる回折次数によって生成される異なるレーザーライン又は弾性散乱ライン間の距離を検出ユニット227、329、723上で測定することによって較正される。更に別の実施形態において、レーザーラインが狭帯域フィルターによって吸収され、較正のために利用できないとき、既知のブリルアン特性の基準物質を用いて、光学配置のスペクトル分散特性を計算する。
【0045】
図7において例示されるような構成の利点及び制約を理解するために、角度のある幾何学的配置のスペクトル効率が、
図1において例示されるような共焦点エピ検出構成を使用する従来の構成と比較された。照明ビームが垂直偏光され、それゆえ、散乱面に対して垂直である場合には、微分断面積は散乱角から独立しており(非特許文献5)、それゆえ、直交する照明経路-検出経路を備える構成は、エピ検出と同じ微分断面積を有する。しかしながら、角度のある幾何学的配置は結果として、単一の点あたりの幾何学的効率が減少することになる(すなわち、測定の実質的な並列化に起因する利点を排除する)。収集された散乱電力はP=I
ill・V・Ω・Rと書くことができる。ただし、I
illは照明光の強度であり、Vは散乱の相互作用体積であり、Ωは収集立体角であり、Rは散乱係数であり、それはm
-1の単位を有し、ここでは定数と見なすことができる。直交構成において、相互作用体積は、半径r=0.61λ/NA
col及び長さl=0.61λ/NA
illを有する円柱形によって近似することができる。ただし、NA
col及びNA
illはそれぞれ、照明対物レンズ及び収集対物レンズのNAである。収集立体角は、収集開口数Ω=πNA
col
2によって決まる。それゆえ、収集散乱電力は、P
90=I
ill・R・π
2・0.61
3・λ
3/NA
illである。エピ構成の場合、代わりに、相互作用体積は約0.61
2・λ
3/NA
epi
4である。ただし、NA
epiは対物レンズのNAであり、収集立体角はΩ=πNA
epi
2である。したがって、収集散乱電力は、P
epi=I
ill・R・π
2・0.61
2・λ
3/NA
epi
2である。同じ照明強度の場合、2つの構成間の収集電力の比は、結局、η=P
90/P
epi=0.61・NA
epi
2/NA
illになる。直交構成では、視野にわたって均一な長い照明ビームラインを生成するために、一般に低いNA
illが好ましい。例えば、1.737mmの使用可能なレイリー範囲に対応するNA
ill=0.0175が使用された。本発明による特定の実験構成をNA
epi=0.1を有するエピ構成(x方向において同じ分解能を与える)と比較すると、その比は34.8%である。この計算は、対物レンズの異なる組み合わせが使用された他の研究(Scarcelli, G.、Kim, P.及びYun, S.H.「In vivo measurement of age-related stiffening in the crystalline lens by Brillouin optical microscopy」Biophysical Journal 101, 1539-1545 (2011))による予測と一致する。重要なことに、この計算は、
図2の実験結果、そして他のエピ検出結果(Guilluy, Cら「Isolated nuclei adapt to force and reveal a mechanotransduction pathway in the nucleus」Nature Cell Biology 16, 376 (2014))と一致する。直交及びエピ検出における散乱効率は、照明体積と検出体積との間の重なりによって決定され、それゆえ、並列化の利点を最大化するために、角度のある幾何学的配置を、単一点測定あたりエピ検出と同じ効率になるように設計することができる。ライン並列検出を含むとき、多重ブリルアン分光法は、エピ検出における1時間超の長さと比べて、数十秒内に数ミクロンの分解能でmmサイズのサンプルの走査を成し遂げることができる。ビームラインに沿って測定可能なサイズは、照明NAと、カメラのピクセル数との両方によって決定される。空間分解能は、検出NA(x方向、ビームラインに沿っている)及び照明NA(y方向及びz方向)によって決定される。
【0046】
背景除去に関して、角度のある構成は、エピ検出より本質的に低い背景雑音を有する。一般の共焦点構成では、照明光の後方反射が分光計の中に容易に結合され、背景雑音の一因になる。ライン走査構成では、照明経路及び検出経路が直交して配置されるので、後方反射を完全に回避することができる。さらに、ラインビームによって照明される領域のみが励起されるので、ライン走査構成は、光学的分割(optical sectioning)に関して、共焦点構成と類似の機能を有する。一方、1ステージVIPAのみが使用されるので、分光計の消光比は制限される。これは、界面、又はイントラリピッド媒質若しくは組織等の光学的に不透明なサンプルを測定するのを難しくする。機器の全体的な消光比を改善するために、ライン走査構成を、アポディゼーション、狭帯域バンドパスフィルタリング、及びガスチャンバー狭帯域吸収フィルタリング等の、背景雑音を抑圧することができる既存の方法と組み合わせることができる。
【0047】
本発明の別の態様は、細胞内物理的(機械的)特性に基づいて、生細胞を分類するシステム及び方法に関する。一例として、限定はしないが、分類されることになる細胞は、ヒト細胞、動物細胞又は植物細胞とすることができる。細胞内物理的情報は、細胞内部の音響特性に関連するブリルアン光散乱のスペクトル分析から取得される。分類は、細胞集団又は個々の細胞のいずれかに基づくことができる。一実施形態において、細胞は、懸濁された状態で、及び/又は2D基体に付着している状態で、及び/又は3D合成/天然細胞外マトリックス内で培養された状態で分析することができる。更に別の実施形態において、細胞は、静止状態において、又は細胞流動を介して分析することができる。
【0048】
一実施形態において、細胞はin vitroで培養及び調製することができ、マイクロ流体デバイスを用いてプローブビームに送達することができる。細胞の異なる領域における細胞内情報を、標準又はカスタムマイクロ流体デバイスを利用する幾つかの構成において取得することができる。細胞は、フローサイトメトリー又は他の細胞分析技法に関連付けられる概ね全ての構成において分析することができる。一例として、限定はしないが、細胞は、マイクロ流体チャネルの中を流れている間に、及び/又は細胞が懸濁液内に静止したままであるときにマップされた状態で、及び/又は既知の方法に従って3Dゲル又は他の細胞外マトリックスの小滴内に封入された状態で、及び/又は細胞凝集体を形成している間に測定することができる。代替的には、マイクロ流体デバイスを用いて、又は用いることなく、細胞は、基体上に付着している間に、又は3Dゲル内にある間に、又は凝集体を形成している間に、明視野若しくは共焦点結像手法で細胞内情報のブリルアン取得を誘導することによって、他の構成(皿、又はフラスコ、又はプレートアレイ内等)において分析することができる。取得された物理的フェノタイピングに基づいて、細胞周期の異なる段階における異なるタイプの細胞又は同じタイプの細胞を区別することができる。一実施形態において、薬物への細胞応答を分析することができる。分析は、細胞の集団、又は個々の細胞のいずれかにおいて実行することができる。
【0049】
図12は、ブリルアン信号を励起/収集するエピ検出経路、監視結像手法、マイクロ流体デバイス(又は他の細胞播種(cell plating)デバイス)、ブリルアン分光計、及びこれら全てのサブ構成要素をリンクする、ソフトウェアを備えるコンピューターの新規の組み合わせから生じる本発明の例示的な実施形態を示す。第1の到来する照明ビーム1210が、ビームスプリッター1211及び対物レンズ1212を使用することによって、マイクロ流体デバイス1230内部の1つにスポット上に合焦する。照明レーザービーム1210はレーザービームとすることができる。一実施形態において、第1の到来するレーザービーム1210は1μm~10μmのスポットの中に合焦する。後方ブリルアン散乱光が同じ対物レンズ1212によって収集され、ビームスプリッター1211、反射ミラー1213及びカップリングレンズ1214を介して、ブリルアン分光計1215に導光される。ブリルアン分光計1215は、ブリルアン周波数シフトを正確に測定することができるデバイスである。ブリルアン分光計1215の一例が非特許文献6(Scarcelli Polachech 2015 Nature Methods)において開示されている。マイクロ流体デバイス1230の上方に、第2の光源1220(一般にレーザービーム1210とは異なる波長を有する)、ビームスプリッター1221、チューブレンズ1222、対物レンズ1224及び2D画像記録デバイス1225からなる明視野顕微鏡が存在する。この顕微鏡構成の目的は、マイクロ流体デバイス1230内部の細胞の流動状況を監視することである。一実施形態において、その顕微鏡構成はスペクトルフィルター1223を更に含み、スペクトルフィルターは、レーザービーム1210からの光を遮断するが、光源1220からの光が通過できるようにする。対物レンズ1212からの合焦スポットが対物レンズ1224の像面上にあるのを確実にするように、対物レンズ1212及び1224が調整される。
【0050】
一実施形態において、ソフトウェアインターフェースを用いて、2D画像記録デバイス1225によって取得された明視野像と、ブリルアン分光計1215のブリルアン信号とを同期させることができる。これにより、ブリルアン信号を測定サンプルの正確な場所に割り当てることができる。
【0051】
図13は、
図12のブリルアン分光法構成とともに使用されるインターフェースのスナップショットを示す。インターフェースの主な機能は、その中に細胞を有するマイクロ流体デバイスの未処理のブリルアン信号1310(空間スペクトルパターン)及び明視野像1320を同時に記録することである。2つの白色の点1310はブリルアン信号(空間スペクトルパターン)であり、点間の距離はブリルアンピークを表す。レーザーラインからブリルアンピークまでの距離は、所望の周波数シフトである。この表現において、ここで適切な較正が行われた場合、2つのブリルアンピーク間の距離を計算することによって、ブリルアンシフトを測定することができる。この場合に、2つのピーク間の距離が短いほど、大きな周波数シフト、それゆえ、高い剛性に対応する。
【0052】
図14において、ブリルアン信号検出を実際にいかに実施できるかについての非限定的な例が示される。具体的には、
図14a~
図14cはそれぞれ、0.85sec、1.35sec及び1.65secにおいて同時に取得された未処理のブリルアン信号(空間スペクトルパターン)及び明視野像の両方のスナップショットである。そのスナップショットは、上から下にチャネル内部を流れる細胞を示す。丸1402は、細胞が横切っている照明(プローブ)ビームスポットの場所を示す。
図14d~
図14fは、
図14a~
図14cの未処理のブリルアン信号から計算されたブリルアン周波数シフトの時間トレースを示す。
図14a~
図14cの同期した画像から、ブリルアン信号の元の場所を、曖昧にすることなく特定することができる。ブリルアン周波数シフトを特定の細胞内の場所に割り当てること以外に、画像処理技法を組み込むことと、細胞容器を並進させること、又はブリルアンプローブビームの位置を制御することのいずれかとによって、この能力を用いてブリルアンプローブビームを特定の場所に導光することもできる。
【0053】
1つの非限定的な実施形態において、細胞を合焦したプローブビームスポットに位置合わせするために、
図12に示されるようなマイクロ流体デバイス1230を2D並進ステージ上に配置することができる。一例として、限定はしないが、マイクロ流体デバイス1230は、長方形の断面を有するか、又はフローサイトメータの大部分と同様に円形の断面を有する一直線のチャネルとすることができる。チャネルのサイズ(一例として、100μm×100μm)が、細胞のサイズ(一例として、直径が10μm~20μm)よりはるかに大きい場合には、細胞は、マイクロ流体デバイス1230を通ってランダムな場所に流れる場合があるので、測定精度、そして潜在的には測定スループットが減少する。一実施形態において、全ての細胞を照明ビーム1210によって探査できるように、細胞をチャネル内部の同じトレースに沿って位置合わせするために、シースフローと呼ばれる流体力学的合焦技法を実施することができる。これは、各細胞の中心がプローブビーム1210の合焦スポットに絶えず位置合されるのを確実にする実効的な方法である。
【0054】
図15a及び
図15bは、シースフロー技法を使用することによって取得された実験データ(
図15b)と、シースフロー技法を使用することなく取得された実験データ(
図15a)との比較を示す。流量及び分光計1215の取得時間を調整することによって、細胞が光ビーム1210を横切って進むときに、1つの細胞の幾つかの場所を探査することができる。
図15a及び
図15bのプロットに関して、単一の細胞の平均剛性を記述するために、細胞内で取得された全てのブリルアン周波数シフトの平均が計算された。それゆえ、
図15a及び
図15bにおいて、ヒストグラムの各データ点は1つの細胞を表す。それらの結果は、非シースフローによって引き起こされる人工的な広がりが取り除かれるので、測定データの広がり(分布の線幅)が、シースフロー技法を使用することによって約40%狭くされる。1つの非限定的な実施形態において、NIH 3T3細胞株が実験において使用された。
【0055】
更に別の実施形態において、画像誘導位置合わせ(image-guided alignment)と呼ばれる第2の手法を用いて、流動する細胞をレーザービーム1210の合焦スポットに位置合わせすることができる。その概念の概要が
図16に示される。この手法も
図12に示される構成内で実施することができるが、以下の変更が加えられる。
図12において、レーザービーム1210の合焦スポットは、顕微鏡1212の視野の中心に位置する。しかしながら、
図16において、合焦スポット1620は、マイクロ流体チャネル1610の底部に動かされ、x方向に沿った中間点に、すなわち、x
0に設定された。細胞1630がy方向に(上から下に)流れていると仮定して、細胞が視野内に現れると、その場所(x
1)及び合焦スポット1620に対する水平シフト(x
0-x
1)を特定することができる。その後、マイクロ流体チャネル1610を搬送する並進ステージを自動的に、又は手動で動かして、細胞1630を合焦スポット1620に完全に位置合わせすることができる。代替的には、更に別の実施形態において、ブリルアンプローブビームは、レーザー走査補助部品(例えば、検流計スキャナー、ポリゴン)又は他の既知のビーム位置決め技法を用いて調整することができる。この手法の成功は、2つの条件に基づく。第一は、流動実験の全てが層流であることであり、それは、細胞がチャネルに入ると、横方向にドリフトするのではなく、そのトレースを維持するのを確実にする。第二は、細胞及びプローブの場所を調整するフィードフォワードアルゴリズムが、チャネルの位置調整を行うだけの十分な時間があるように流速より速いことである。
図12のブリルアン分光法構成を使用することによって取得された典型的な実験結果が、
図17aに示される。細胞が培養され、実験直前に緩衝液の中に懸濁された。一実施形態において、緩衝液として、リン酸緩衝生理食塩水(PBS)が使用された。流量を調整することによって、細胞が光ビームを横切って進むときに、各細胞内の5~10の位置が探査される。この場合に、測定データ点の全てを用いて、
図17a内のヒストグラムグラフをプロットした。
図17aのヒストグラムは3つの部分からなり、それらの部分の輪郭が3つの実曲線1702、1704及び1706によって描かれる。第1のピーク1704はPBSのシグネチャであり、それは
図17bに示されるように、容易に認識し、位置を特定し、除去することができる。光ビーム1210は、2つの隣接する細胞間の時間間隔内でのみ緩衝液を探査することになるので、測定データ内のPBSのシグネチャの出現が予想される。1つの非限定的な実施形態において、上記のように、細胞流動実験のために、NIH 3T3細胞株が使用される場合がある。
【0056】
細胞内情報を抽出するために、曲線当てはめ法を用いて、
図17aから細胞のシグネチャを抽出することができる。3つのガウス分布1702、1704及び1706の線形重ね合わせを用いて、元のヒストグラムを当てはめた。PBS溶液のシグネチャ1704は、
図17bのデータをあらかじめ当てはめることによって十分に特定することができるので、この情報は、当てはめの既知のパラメーターとして使用することができる。
図17aの元のヒストグラムに最も適合するように破曲線を当てはめた後に、
図17cに示されるように、PBS溶液1704のシグネチャを差し引くことによって、細胞の情報が取得された。
図17cのヒストグラムは2つのピークからなり、それらのピークは、2つのガウス分布1708及び1710の重ね合わせである。実際には、ガウス分布1708及び1710は、細胞内の異なる領域(すなわち、細胞質及び細胞核)からの機械的シグネチャを表す。
【0057】
図17a内のフローデータが、細胞集団の機構を正しく表すことを検証するために、細胞集団のフローデータと、細胞の組の2D画像とが直接比較された。最初に、細胞が
図12のマイクロ流体デバイス1230の中に送達され、その後、2~3分間、チャネル内部で細胞が懸濁されるのを確実にするために停止された。比較として、同じ細胞内の複数の場所において、同一の方法で調製された同じ細胞株を用いてブリルアン周波数を取得することに基づく流動実験が、
図12のブリルアン分光法構成を使用することによって行われた。
図18aは、チャネルを通って流れている細胞がレーザーを横切るラインに沿った一群の点に対して測定されたブリルアン周波数シフトの分布を示すヒストグラムである。一実施形態において、577個の細胞(3546個の細胞内の点)が測定された。
図18aのヒストグラムは、2つのガウス分布1802及び1804の重ね合わせによって当てはめられる。
図18bは、2つのガウス分布1806及び1808の重ね合わせによって当てはめられた流動細胞の2Dブリルアン画像から生じるヒストグラムである。一実施形態において、全部で4378個の細胞内の点を含む、29個の細胞が測定された。
図18a及び
図18bのヒストグラムは、2ピークプロファイル及び各ピークの周波数シフトを含む、同じ特徴を共有することがわかる。上記で規定されたような細胞及び測定値の数は、限定するものではなく、任意の数の細胞及び測定値を用いて、
図18a及び
図18bのヒストグラムを取得できることに留意されたい。
【0058】
さらに、
図17cの2つのピークが、細胞質/細胞骨格及び細胞核の機械的シグネチャに対応することが検証された。細胞の細胞核のみを染色するために蛍光染料が使用され、その後、同じ細胞の明視野/蛍光画像(
図19a及び
図19b)並びに2Dブリルアン画像(
図19d)の両方が取得された。
図19a及び
図19bの画像を融合することによって取得された
図19cの画像を
図19dのブリルアン2D画像と比較することによって、細胞質から細胞核を分離することができる。細胞核及び細胞質の両方のヒストグラムが
図19eに合わせてプロットされており、ピーク1904は細胞核に由来し(
図19bの蛍光領域に対応する)、ピーク1902は細胞質に由来する(
図19bの非蛍光領域に対応する)。
図19eと
図17cとの間の比較に基づいて、細胞がプローブビームを横切るときに、
図12のブリルアン分光法構成を用いて細胞内の複数の点においてブリルアン周波数を取得する流動実験によって、細胞集団内で細胞核及び細胞質の機械的シグネチャを区別できると結論付けられる。蛍光/ブリルアンを組み合わせた構成を用いて、蛍光で標識される特定の細胞内成分のブリルアン取得を正確に割り当てることができるか、又は画像誘導することができる。
【0059】
本発明の別の態様において、細胞内物理的(機械的)フェノタイピングに基づいて細胞を分類する能力が例示される。例えば、細胞に及ぼすサイトカラシンD(CytoD)の影響が分析された。細胞の全体剛性(変形能)は、細胞骨格及び細胞核の両方によって決まることが広く受け入れられている(Fletcher DA、Mullins RD「Cell mechanics and the cytoskeleton」Nature 463, 485-492 (2010))。細胞骨格は3つの主な重合体、すなわち、アクチン、微小管及び中間径フィラメントからなる。F-アクチンはアクチンの重合した状態であり、細胞形状及び機械的特性を制御する際に主要な役割を有する。それゆえ、F-アクチンの重合及び解重合が、細胞形状及び剛性の変化を推進する有向力を生成する。CytoDは、F-アクチンの重合を抑制する薬物である。F-アクチンはCytoDで処置された付着細胞において著しく低減されることになり、結果として、細胞全体がその剛性を低下させることが報告された(Wakatsuki T、Schwab B、Thompson NC、Elson EL.「Effects of cytochalasin D and latrunculin B on mechanical properties of cells.」J Cell Sci. 114, 1025-1036 (2001))。細胞内レベルにおいて、細胞核及び細胞骨格はいずれも、その両方内にアクチンが存在するときに、剛性を低下させることになる。さらに、細胞核及び細胞骨格は機械的に接続されるので、細胞核は細胞骨格からのプレストレスを受ける(Wang N、Tytell JD、Ingber DE.「Mechanotransduction at a distance: mechanically coupling the extracellular matrix with the nucleus」Nat Rev Mol Cell Biol. 10, 75- 82 (2009))。したがって、F-アクチンの崩壊は、最終的に細胞核を軟化させるようなプレストレスも低減するはずである(Chalut KJ、Hoepfler M、Lautenschlaeger F、Boyde L、Chan CJ、Ekpenyong A、Martinez-Arias A、Guck J.「Chromatin decondensation and nuclear softening accompany Nanog downregulation in embryonic stem cells」Biophys J. 103, 2060-2070 (2012))。
【0060】
一実施形態において、CytoDは細胞上に塗布され、その後、トリプシンを用いて剥離される。その後、処置された細胞はPBSに再懸濁され、
図12の構成において測定される。投与量依存性を観測するために、CytoDの3つの投与量が使用された。一例として、限定はしないが、CytoD投与量は、0.5ug/ml、2ug/ml及び5ug/mlであった。処置された細胞(279、232及び577)は投与量ごとに探査され、対照のために、処置されない細胞(268)が探査された。細胞全体にわたる平均値を用いて、細胞の機械的フェノタイピングを表した。
図20aの結果は、細胞の剛性へのCytoDの投与量依存効果を示す。したがって、細胞全体にわたる平均ブリルアン周波数シフト値を用いて、高スループットの機械的フェノタイピングを取得することができる。1つの非限定的な実施形態において、NIH 3T3細胞株を用いて、実験を実行した。
【0061】
平均された細胞値の他に、
図17a~
図17cを参照しながら説明された方法は、細胞内情報を取得することができる。CytoDの異なる投与量に対応する
図20b~
図20eに示されるように、ピーク2004、2008、2012及び2016のブリルアン周波数シフトは小さくなり、それゆえ、
図20bから
図20eに左に移動した。これは、高い投与量のCytoDによる処置ほど細胞核が軟化していることを示している。一方、細胞質を表すピーク2002、2006、2010、2014は、この実験において概ね変化しないことに留意されたい。これは、フラスコ基体から細胞を剥離するために使用されるトリプシン処理ステップが全てのF-アクチンストレス線維を崩壊させるので、結果として、懸濁された細胞に関する細胞骨格に測定可能な影響をほとんど及ぼさないためである。付着細胞において、細胞骨格剛性が著しく低減されることが検証された。
【0062】
図21は、本発明の別の実施形態による、ブリルアン分光法構成を示す。この実施形態において、細胞がブリルアンビームによって探査されると同時に、一緒に記録される細胞の共焦点蛍光画像を取得するために、第2の光学手法が追加される。第2の光学手法は、ダイクロイックミラー2127と、チューブレンズ2122と、ピンホールマスク2126と、光検出器2125とを備えることができる。ダイクロイックミラー2127は、レーザービーム2110及びブリルアン信号を通過させるが、蛍光を反射する。
図12の明視野顕微鏡と比べて、共焦点構成は、焦点外れの面からの蛍光を除去することによって、より良好な横断及び横方向光学分解能を有する。実際には、蛍光標識された細胞を用いて、細胞のサブ領域(例えば、細胞核)を識別することができ、レーザービーム2110の合焦スポットを用いて、蛍光放射を直接励起することができる。その後、共焦点構成を用いて、カメラ2125によって蛍光放射光を検出する。蛍光放射光の強度に基づいて、ブリルアン信号を目標点に割り当てることができるように、細胞内部の合焦スポット(それゆえ、ブリルアン信号)の場所を正確に特定することができる。この手法を用いて、測定中にブリルアンビームを細胞の特定の領域に導光することもできる。更に別の実施形態において、反転構成に加えて、共焦点顕微鏡の直立構成を使用することもできる。この場合、共焦点構成は、マイクロ流体デバイス2130の上に存在することになる。また、蛍光情報を用いて、細胞の縁部における測定値を較正するために、3次元結像ボクセル内に存在する非細胞媒質(例えば、PBS)の量を推定することができる。
【0063】
図22の画像は
図21の実施形態によるブリルアン分光法構成を使用することによって取得されたが、共焦点顕微鏡法ではなく、明視野蛍光が使用された。蛍光染料を用いて、マイクロ流体デバイス2130の中に後に送達される細胞の細胞核を染色した。細胞が懸濁液内に静止したままになると、マイクロ流体デバイスを動かすことによって、単一細胞の2Dマッピングが実行された。ブリルアンビームスポット及び蛍光が空間内で互いに重ね合わせられるので、ブリルアン信号を測定点ごとに蛍光強度に結び付けることができる。その結果を表現する1つの取り得る方法が
図22に示される。各点が細胞内の或る特定の領域において取り込まれた1つの測定値を表し、異なる色は、データ点の相対密度を符号化する。散乱スポットから、横軸上にプロットされた蛍光強度に従って、各ブリルアン信号の場所を明確に見定めることができる。例えば、破線の丸はそれぞれ、細胞核(蛍光領域)及び細胞質(非蛍光領域)からの測定点を示す。
【0064】
図21の実施形態を使用して細胞を分類する一例として、細胞に及ぼすCytoDの影響が分析された。処置する場合、及び処置しない場合の細胞のブリルアン及び蛍光両方の2D結像が行われ、
図22の方法を用いて、データを分析した。その結果が
図23に示されており、(a)及び(b)はそれぞれ、対照群及びCytoD処置群の結果である。細胞がCytoDで処置された後に、細胞核ははるかに軟化することが観測された。
【0065】
図12及び
図21の実施形態において、ブリルアン信号が逐点走査モードにおいて取得され、そのモードは全体的なスループットの制限である。その場合に、個々の細胞を分析するのに十分な、細胞あたりの細胞内点を取得するのは難しい場合がある。したがって、
図12及び
図21の実施形態は、細胞集団分析に関して有用であるが、単一細胞の機械的フェノタイピングに関して強力というほどでない場合があり、それゆえ、細胞選別の適用例には適していない。この問題への1つの解決策は、単一の細胞を分析するのに十分な細胞内点を有するように流速を下げることである。別の解決策は、2つ以上の分光計を利用して、一度に2つ以上の場所を測定することである。代替的には、
図7によるブリルアン分光法構成が、細胞の複数の点を同時に測定することができる。
図7の光学配置は、
図12及び
図21の配置に比べて決定的に重要な変更を有する。この場合の主な変更は、CCDカメラによって記録された湾曲したブリルアンスペクトルパターンを一直線にするために、一対の対物レンズ713及び714が高分解能において利用されることである。これは、取得されたブリルアン信号が照明ラインビームに沿った対応する場所に正確に割り当てられるのを確実にするので決定的に重要である。具体的には、レーザー源710からのビームは最初にミラー724によって導光され、その後、x方向において照明ラインビームになるように対物レンズ711によって合焦する。ラインビームは、マイクロ流体デバイス712の縁部から、そのチャネルの中に差し込む。細胞は、y方向に沿って流動することによって、ラインビームを通り抜けることができる。ラインビームによって励起されたブリルアン散乱信号は、ブリルアン分光計の並列構成によって検出され、その構成は、一対の対物レンズ713及び714と、コリメートレンズ715と、円柱レンズ716と、分散デバイス(すなわち、VIPA)717と、円柱レンズ718及び719と、マスク720と、球面レンズ722と、高感度CCDカメラ723とからなる。
【0066】
図12及び
図21の分光計とは異なり、
図7の多重化ブリルアン分光計は、照明ラインビームに沿って数多くの点からのブリルアン信号を同時に取得することができる。例えば、適切な対物レンズを選択し、分光計内の光学系を調整することによって、x方向において1μm未満の分解能を達成することができる。細胞の通常のサイズは約10μm~15μmであるので、細胞が流動しているときに、1ショットの場合に約10個の点を、単一の細胞の場合に最大で100個を超える点を取り込むことができる。このデータセットは、単一の細胞の細胞内機構を取得するのに十分な情報を確実に与えることになる。それゆえ、
図7の実施形態は、その細胞内機械的フェノタイピングに基づいて、単一の細胞を分類する、より実効的な方法を提供する。幾つかの実施形態において、照明経路及び検出経路が直角に配置されることに加えて、
図7の構成は、45度及び60度の配置角度を含む、0度より大きい任意の角度において機能することもできる。
【0067】
図24は、
図7の実施形態の概念図を示す。照明レーザーライン2410は左縁部(x方向)からマイクロ流体チャネル2430を照明し、細胞2420は、チャネルの内部をy方向に沿って流れる。細胞2420がレーザーライン2410に当たると、細胞2420の照明された領域からのブリルアン信号が同時に取得されることになる。細胞が通り抜けるとき、そのラインビームによって細胞体全体を走査することができる。最終的に、
図24に示されるように、細胞2420の2D画像2440を再構成することができる。
【0068】
図7の実施形態は、細胞又はチャネルから反射されるか、又は散乱するレーザーからのスプリアス光を除去するために、フィルタリング要素を追加することによって更に改善することができる。このために、ファブリ-ペローフィルターを使用することができる。ファブリ-ペローフィルターは、実際には、選択された周波数を有する光が通過できるようにし、他の光を遮断できるようにするエタロンである。フィルターを適切に設計することによって、フィルターの透過ピークを、ブリルアン信号と重なり合うように調整することができる。スプリアス光は、ブリルアン信号とは異なる周波数を有し、抑圧されることになる。しかしながら、
図7の方法は、位置合わせの問題及び潜在的なアーティファクトを導入する場合があり、正確に較正することによって解決されるべきである。
【0069】
代替的には、吸収ガスチャンバーを使用することができる(Piironen P、Eloranta EW.「Demonstration of a high-spectral-resolution lidar based on an iodine absorption filter」Opt. Lett. 19, 234-236 (1994))。ガスチャンバーは、ビームの方向変更(redirection)によるスペクトル分散を導入せず、それゆえ、並列分光計に極めて適している。ガスチャンバーは一般にガラスから製造され、明確に規定された吸収スペクトルを有する特定の原子元素からの蒸気を含む。光は入口を介してガスチャンバーを通り抜け、いかなる光学歪みも導入することなく、窓から出る。吸収ガスは、特定の光学波長において非常に狭い(1GHz未満の)分子吸収スペクトルを有する。吸収ガスを適切に選択することによって、その吸収ラインが光源710の波長と完全に重なり合うガスチャンバーを設計することができる。例えば、ガスチャンバーは、レンズ712と713との間に配置することができ、ブリルアン信号が通り抜ける間に、スプリアス光(レーザー光源710と同じ波長を有する)が吸収されることになる。吸収スペクトルの線幅は非常に狭いので、レーザーの周波数の任意の小さなドリフトがあると、ガスチャンバーの実効性がなくなる。ガスチャンバーを
図7の実施形態とともに十分に機能させるために、レーザー源の波長をガスチャンバーの吸収ラインにロックする必要がある。一例として、
図25は、ルビジウムガスを含む、ガスチャンバーの測定吸収スペクトルを示す。このガスチャンバーは780.24nmの波長において機能し、その結果は、チャンバーが94.5℃まで加熱されたときに、ピーク吸収は約34dBであり、挿入損失が2dBであることを示す。得られたスペクトル減衰によって、強い背景雑音を有する生体サンプルを測定できるようになる。
【0070】
例4
本発明の例示的な構成が
図26に示される。垂直偏光レーザービーム2628が最初に偏光ビームスプリッター(PBS)2634によって反射され、対物レンズ2632によって底部側からマイクロ流体チャネル2624の中央に合焦する。一実施形態において、ブリルアン散乱を励起するための光源として、11mW単一モード直線偏光532nm cwレーザー(Torus、Laser Quantum社)が使用された。4分の1波長板(QWP)2626に起因して、合焦したビームは円偏光される。後方散乱光が同じ対物レンズ2632によって収集され、QWP2626に再び通される。円偏光した光は最終的に、PBS2634によって完全に透過する水平偏光した光に変換される。散乱光は、その後、ミラー2630によって導光され、ファイバーカップリングレンズ(FC)2636によってブリルアン分光計2638に結合される。ブリルアン分光計2638は、交差軸構成において、標準的な2ステージ仮想画像化フェーズドアレイ(VIPA)によって構成される。マイクロ流体チャネル2624の上側に、チャネル2624内部の細胞の流動状況を監視するために、明視野顕微鏡が構築された。破線矢印は、チャネル2624内部の細胞の流動方向を示す。照明のために赤色LED2616が使用され、ブリルアンビームを遮断するためにロングパスフィルターが使用された。一例として、限定はしないが、LED2616の波長は660nmとすることができ、ロングパスフィルターはFEL0600、Thorlabs社とすることができる。合焦したブリルアンビームの場所は、ロングパスフィルター(LP)2612を取り外すことによってあらかじめ決定された。ブリルアン信号に及ぼすLED光2616の影響を排除するために、ブリルアン分光計2638の前方に狭帯域フィルターが配置された。一実施形態において、分光計のデータを測定サンプルの対応する場所に割り当てることができるように、CCDカメラ2610及びブリルアン分光計2638をLab VIEWプログラムによって同期させた。一実施形態において、実験中に、分光計のサンプリング時間が50msと設定された。
図27a~
図27dは、細胞内の複数の点においてブリルアン周波数を取得することに基づく、細胞流動実験の測定手順を示す。実験が行われているとき、ブリルアン分光計2638は絶えず信号を記録する。細胞がビームの場所を通り抜けるとき(矢印によって示される白色の点)、そのブリルアンシフトが測定される。
【0071】
マイクロ流体デバイス
図26の構成において使用されるチップは、入口においてシースフロー構成を備える単一の一直線のチャネル2624を有する。一実施形態において、チャネルの断面は150μm(幅)×50μm(深さ)のサイズを有し、長さは約50mmである。チップは顕微鏡の2Dステージ上に組み付けられ、細胞懸濁液をチャネルの中に送り込むためのシリンジポンプに接続される。シースフロー技法を用いて、細胞をチャネルの中央に位置合わせすることができる。一実施形態において、流速は約20μm/s~40μm/sである。全体的なスループットは一時間あたり細胞約160個から300個まで変化し、平均では一時間あたり細胞230個である。
【0072】
サンプル調製
細胞が培養され、Mg2+及びCa2+不含のリン酸緩衝生理食塩水(PBS)の溶液内に再懸濁される。一実施形態において、最終的な細胞濃度は、100万細胞/ml~1000万細胞/mlである。クロマチン脱凝縮実験の場合、古い培地が吸引され、完全な成長培地内の100ng/mLトリコスタチンA(TSA)溶液と置き換えられる。その後、細胞は2時間、インキュベートされる。インキュベーション後に、細胞が、0.25%トリプシンEDTAを用いて採取され、250×gにおいて5分間、遠心分離され、約250万細胞/mLにおいて、Mg2+及びCa2+不含のPBSに再懸濁された。細胞核を蛍光染色する場合、Vybrant DyeCycle Ruby染料(ThermoFisher社)が使用され、推奨プロトコルに従った。要するに、50万細胞/mLの濃度を有する細胞懸濁液に1μL染料が添加され、十分に混合される。最終的な染料濃度は5μMである。その後、その混合物は37℃で15分~30分間、インキュベートされ、光から保護される。その際、細胞は、洗浄することなく、蛍光実験の準備ができている。一実施形態において、実験のために、NIH 3T3細胞株が使用される。
【0073】
細胞の2次元画像
本発明の一実施形態において、細胞懸濁液が、シリンジポンプを用いて、マイクロ流体チャネルの中に徐々に送り込まれる。チャネル内に細胞が満たされると、ポンプが停止される。数分後に、細胞はチャネルの底部に沈殿するが、依然として、丸みを帯びた形状を保持する。その後、レーザービームスポットの高さが細胞の中心に調整され、顕微鏡のステージを走査することによって、水平面において2D走査が実行される。走査のステップサイズは、両方の次元において0.5μmと設定される。一実施形態において、実験のために、NIH 3T3細胞株が使用される。
【0074】
元のフローデータからの細胞のシグネチャの抽出
緩衝液内に懸濁された細胞は一般に丸い形状を有する。一実施形態において、細胞の直径は10μmから20μmに及ぶ。合焦したブリルアンビームは約0.5μmのスポットサイズを有するので、細胞がビームスポットを横切って流れるときに複数の位置を測定することができる。一実施形態において、流動実験中に、ブリルアン分光計は、20Hzの速度において絶えずデータを取得する。
図27a~
図27dは、測定手順を例示する。具体的には、
図27a~
図27cは、異なる時点(t1、t2及びt3)において細胞がマイクロ流体チャネルを通って流れる場合を示す。スケールバーは30μmである。矢印によって示される白色の点は、レーザービームスポットの場所である。
図27dにおいて、ブリルアン信号の時間トレースが提示される。点は測定データであり、実線は見やすくするために引いた線である。
【0075】
データの後処理において、明視野画像の助けを借りて、時間窓を用いて元のデータストリームから細胞事象を選出することができる。
図28a~
図28cは、577個の細胞について測定されたブリルアン周波数シフトの分布を表す。各データ点は、ブリルアン周波数シフトを測定するブリルアン分光計の1つの測定事象に対応する。細胞の機械的情報を保持するために、時間窓は一般に、各細胞事象の継続時間より広い。それゆえ、細胞及びPBS緩衝液両方の機械的シグネチャが記録され、
図28aに示される。ヒストグラムの第1のピーク2802はPBS緩衝液のシグネチャであり、それは、そのブリルアン周波数シフトに基づいて容易に識別することができる。その際、このシグネチャは、PBS緩衝液のみを流すことによって特徴付けられ、その結果が
図28bに示される。そのフローデータは正規分布に十分に当てはめることができ、それは、ブリルアンシフトが7.51GHzであり、0.047GHzの幅を有することを示す。
図28aからPBS緩衝液のシグネチャ2804を除去することにより、直接、細胞の機械的シグネチャを抽出することができ、その結果が
図28cに示される。細胞のシグネチャは2つのピークを明らかにする。測定データを当てはめるために、2つの正規分布2806及び2808の線形重ね合わせが使用された。それらの結果が
図28cにおいて実曲線2806及び2808によって示され、それぞれ7.57GHz及び7.80GHzの中心周波数を有する。細胞核は細胞の最も剛性が高い細胞小器官であることが当該技術分野において既知である(Friedl, P.、Wolf, K.及びLammerding, J.「Nuclear mechanics during cell migration」Current Opinion in Cell Biology 23, 55-64 (2011);非特許文献5)。それゆえ、実曲線はそれぞれ、細胞質及び細胞核のシグネチャである。これは、細胞核を染色し、蛍光画像をそのブリルアンデータと比較することによって確認することができる。さらに、3T3細胞の場合の細胞核の体積比は、30%より大きいことが報告されており、それは、細胞質と細胞核との直径比が約3:7であることを意味する。それゆえ、細胞内の複数の点においてブリルアン周波数を取得することに基づく流動実験中に、細胞核のシグネチャは、細胞質より、ブリルアン分光計によって測定される機会が2倍より大きく、それも
図28cの2つのピーク2806及び2808の振幅によって確認される。
【0076】
フローデータが細胞集団機構を表す
懸濁された細胞は丸みを帯びた形状を有し、それは簡単にするために球として近似することができる。シースフローの助けを借りて、細胞がマイクロ流体チャネルを通って流れるときに、細胞の中心を幅方向においてレーザービームスポットの場所に位置合わせすることができる。チャネルの深さが細胞サイズより大きいので、細胞の中心は、深さ方向において細胞ごとに異なることになる。このような場合には、細胞がレーザービームスポットを横切って進むときに、各単一細胞の直径を含む断面内の1つのラインがサンプリングされることになる。細胞流動中のブリルアン周波数シフトの測定値が、細胞集団の機構を表すのを確実にするために、それらの測定値が、細胞の中心を通る断面の画像と比較される。具体的には、細胞の2D画像が、
図29aに示されるように取得された。
図29bは、2つの曲線2902及び2904の重ね合わせによって当てはめられる細胞集団の全体的な機構を表す、
図29aの画像データの全ての分布を示す。細胞がレーザーを横切るラインに沿った一群の点に対応するブリルアンフローデータの分布が
図29cにプロットされ、曲線2908及び2906の重ね合わせによって当てはめられる。
図29bの2Dブリルアン画像データを
図27cのブリルアンフローデータと比較することによって、両方の分布が2ピークシグネチャを有することが観測される。より重要なことに、両方の分布に関するこれら2つのピークのブリルアンシフトは同じであり、すなわち、7.57GHz及び7.80GHzである。したがって、フローデータが、細胞集団の機械的特性を正確に明らかにすると結論付けられる。
【0077】
ブリルアンシグネチャが細胞核を細胞質から区別する
懸濁された細胞を球であると見なすと、固定されたビームスポットによって1つの直径が検知されるように、細胞がシースフローを介して位置合わせされる。しかしながら、実際の実験において、垂直方向における細胞の場所はわずかに上下変動するので、ビームスポットは実際には、その直径を通る断面をマッピングしている。
図28cの2ピークシグネチャを更に検証するために、蛍光染料を用いて、細胞の細胞核のみを染色した。
図30に示されるように、明視野(I)、蛍光(II)及びブリルアン2D(IV)の画像が取得された。明視野画像及び蛍光画像を含む融合画像(III)をブリルアン画像と比較することによって、ヒストグラム3002及び3004によってそれぞれ示されるように、細胞核のブリルアンシグネチャを細胞骨格のブリルアンシグネチャから容易に分離することができる。集団内の細胞間の小さな差を考えるとき、細胞内の複数の点においてブリルアン周波数を取得することに基づいて流動実験から取得されたデータの分布は、細胞結像からの分布に類似しているはずであり、それが、
図30と
図28cとの間の高い類似度によって確認される。したがって、
図28cの流動データの2つのピークがそれぞれ細胞核及び細胞骨格のシグネチャに対応することが検証される。
【0078】
クロマチン脱凝縮が細胞核を軟化させる
クロマチンはDNA及びタンパク質の複合物であり、細胞核内に染色体を形成する。その機能は、細胞の細胞核に収まるように、かつDNA構造及び配列を保護するように、DNAを小さな体積の中に効率的に詰め込むことである。研究の積み重ねが、クロマチンの凝縮レベルが細胞核の剛性に密接に関連することを示唆する。例えば、クロマチン脱凝縮の結果として、細胞核の剛性が低下する。逆に、クロマチン凝縮の増加が、細胞核の硬化につながる。一実施形態において、機械的サイトメトリーを用いて、TSA、ヒストンデアセチラーゼ阻害剤によって引き起こされるクロマチン脱凝縮の影響を特徴付けた。
【0079】
図31a~
図31cは、クロマチン脱凝縮による細胞核の軟化の実験的な観測結果を与える。1つの非限定的な実施形態において、356個及び353個の細胞が、TSA処置群及び対照(非TSA処置)群からそれぞれ収集される。
図31a及び
図31bに示されるように、細胞がTSAで処置された後に、細胞核のブリルアンシフトが7.80GHz(曲線3108)から7.71GHz(曲線3104)に動かされ、それは細胞核の明確な軟化を示す。一方、細胞質のブリルアンシフトは概ね同じままであり(7.56GHz(曲線3102及び3106))、それは細胞質がTSAによって影響を及ぼされなかったことを意味する。
図31cに示されるように、TSA処置細胞及び対照細胞の細胞質に関する平均ブリルアンシフトにそれぞれ関連付けられる棒グラフ3110及び3112は、同じ高さを有し、ブリルアン周波数に変化がないことを示す。対照的に、TSA処置細胞及び対照細胞の細胞核に関する平均ブリルアンシフトにそれぞれ関連付けられる棒グラフ3114及び3116は、異なる高さを有し、ブリルアン周波数に変化があることを示す。
【0080】
したがって、本発明による機械的サイトメトリーは、オンチップブリルアン技法を用いて、細胞の細胞核シグネチャを識別することができる。それは、細胞集団の細胞核の機械的特性を測定するラベルフリーの非接触、非侵襲的方法を提供する。一実施形態において、総合スループットは一時間あたり約230細胞であり、それは、原子間力顕微鏡法(AFM)及びマイクロピペット吸引等の既存の技法よりはるかに効率的である。さらに、本発明による機械的サイトメトリーを用いて、細胞核の剛性に及ぼすクロマチン脱凝縮の影響を評価することができる。TSA処置後に明確に軟化することがわかっている。細胞核剛性の変化は、遊走、分化及び悪性転換等の細胞の数多くの重要な活動に関与するので、本発明による機械的サイトメトリーは、細胞集団の細胞核機構の効率的で、ラベルフリーの特性評価が必要とされる場合に適用される可能性がある。
【0081】
要するに、
図12~
図31を参照すると、本発明による生体細胞を分類する方法は、各生体細胞内の複数の点においてブリルアン散乱スペクトルを測定することに基づく。測定されたブリルアン散乱スペクトルに基づいて、生体細胞内の異なる空間点において細胞内物理的特性に関連する1つ以上の指標が特定される。生体細胞は、生体細胞内の異なる空間点における細胞内物理的特性に基づいて分類される。
【0082】
一実施形態において、生体細胞内の各測定点を含む、ブリルアン周波数シフトに関するヒストグラムが生成される。そのヒストグラムを当てはめるために、ガウス分布の線形重ね合わせが適用される。次に、ヒストグラム内の各ピークが特定され、各ピークは、細胞内の異なる領域からの機械的シグネチャを表す。生体細胞内の異なる空間点における物理的特性は、特定された機械的シグネチャと相関がある。
【0083】
生体サンプルは、生物有機体、組織又は生体細胞を含むことができる。一実施形態において、生体細胞は生細胞である。更に別の実施形態において、指標としてブリルアン周波数シフトが選択される。さらに、一実施形態において、生体細胞内の異なる空間点において細胞内物理的特性を特定するステップは、画像を形成することと、特定された物理的シグネチャの空間による違いに基づいて、パラメーターを分割することとを含む。サンプルの物理的特性は、粘弾性率、密度、屈折率、電歪及びその組み合わせからなる群から選択される。
【0084】
一実施形態において、生体サンプルを有する容器は、マイクロ流体デバイスのマイクロ流体チャネルである。生体細胞は、スペクトルが測定される間に、マイクロ流体チャネルを通って流れる。
【0085】
一実施形態において、ブリルアン散乱の元の場所を識別するために、そして照明光ビームをマイクロ流体チャネル内の特定の場所に導光するために、ブリルアン光スペクトルパターンと同時に明視野2D画像が取得される。更に別の実施形態において、多次元ヒストグラムを作成して生体細胞を分類するために、蛍光散乱、ラマン散乱、前方散乱及び側方散乱とともに1つ以上のブリルアン指標が使用される。一実施形態において、細胞核を細胞質から分離するために、2D明視野細胞画像及び2D蛍光細胞画像を含む融合画像が、ブリルアン周波数シフトに基づく細胞画像と比較される。
【0086】
ブリルアン周波数シフトは、点走査モードにおいて、又は多重化走査モードにおいて測定することができる。細胞内から取得された抽出された物理的特性は、単一の生体細胞又は生体細胞の集団の分析結果を指す場合がある。
【0087】
一実施形態において、ブリルアンスペクトルを表すヒストグラムは、細胞質及び細胞核にそれぞれ対応する2つのピークを有する。更に別の実施形態において、ブリルアンスペクトルによって規定された細胞の細胞内機械的特性に基づいて、改変された生体細胞が無傷の生体細胞から区別され、細胞は、細胞骨格又は細胞核等の細胞内成分を標的にする薬物によって改変される。細胞内物理的特性を用いて、細胞を選別する。
【0088】
図32は、本発明の一実施形態によるシステムを示すブロック図である。ブリルアン分光計3202が、ブリルアン散乱光の空間スペクトルパターン3404と、ブリルアン光スペクトルとをユーザーインターフェース3208に与える。プロセッサ3206が、ブリルアン分光計3202、ユーザーインターフェース3208、及びメモリ3210と通信する。プロセッサ3206は、メモリ3210内の命令を実行するように構成される。一実施形態において、メモリ3210は、検出された空間スペクトルパターンに基づいて、各測定サンプル点における1つ以上のブリルアン指標を計算するための命令を含む。更に別の実施形態において、メモリ3210は、測定されたブリルアン散乱スペクトルに基づいて、生体細胞内の異なる空間点において細胞内物理的特性に関連する1つ以上の指標を抽出する命令と、生体細胞内の異なる空間点における物理的特性に基づいて、生体細胞を分類するための命令とを含む。さらに、メモリ3210は、ブリルアン周波数シフトに関するヒストグラムをプロットし、ヒストグラムをユーザーインターフェース3208に通信するための命令、ヒストグラムを当てはめるためにガウス分布の線形重ね合わせを適用するための命令、ヒストグラム内の各ピークを特定するための命令であって、ピークは細胞内の異なる領域からの機械的シグネチャを表す、命令、及び媒質に関連付けられるデータをヒストグラムから除去するための命令を含むことができる。
【0089】
本発明を説明する文脈において(特に添付の特許請求の範囲の文脈において)「1つの(a、an)」及び「その(the)」という用語並びに類似の指示物を使用することは、本明細書において他に指示されない限り、又は文脈によって明確に否定されない限り、単数形及び複数形の両方を含むと解釈されるべきである。「備える」、「有する」「含む」及び「含有する」という用語は、他に言及されない限り、オープンエンドの用語(すなわち、「含むが、それに限定されない」を意味する)と解釈されるべきである。本明細書において値の範囲を列挙することは、本明細書において他に指示されない限り、その範囲内に入る別々の各値を個々に参照する簡潔な方法としての役割を果たすことを意図するにすぎず、本明細書において個々に列挙されたかのように、別々の各値が本明細書に組み込まれる。例えば、範囲10~15が開示される場合には、11、12、13及び14も開示される。本明細書において説明される全ての方法は、本明細書において他に指示されない限り、又は文脈によって明確に否定されない限り、任意の適切な順序において実行することができる。本明細書において与えられるありとあらゆる例、又は例示する用語(例えば、「~等の(such as)」)の使用は、本発明の理解をより容易にすることを意図するにすぎず、別に特許請求されない限り、本発明の範囲に制限を課すものではない。本明細書における文言は、任意の特許請求されない要素を、本発明を実施する上で不可欠であると示すものと解釈されるべきでない。
【0090】
本発明の方法及び構成は、様々な実施形態の形で組み込むことができ、その幾つかのみが本明細書において開示されることは理解されよう。上記の説明を読むと、それらの実施形態の変形が当業者に明らかになる場合がある。本発明者らは、当業者が必要に応じてそのような変形を利用することを期待しており、本発明者らは、本発明が本明細書において具体的に説明されたのとは別のやり方で実施されることを意図している。
【0091】
したがって、本発明は、適用法令によって許されるような添付の特許請求の範囲において列挙される主題の全ての変更及び均等物を含む。さらに、本明細書において他に指示されない限り、又はそうでなくても文脈によって明確に否定されない限り、本発明の全ての取り得る変形における上記の要素の任意の組み合わせが本発明によって包含される。