JP2024519990A - 誘導ラマン散乱トモグラフィシステム及び方法 - Google Patents
誘導ラマン散乱トモグラフィシステム及び方法 Download PDFInfo
- Publication number
- JP2024519990A JP2024519990A JP2023572624A JP2023572624A JP2024519990A JP 2024519990 A JP2024519990 A JP 2024519990A JP 2023572624 A JP2023572624 A JP 2023572624A JP 2023572624 A JP2023572624 A JP 2023572624A JP 2024519990 A JP2024519990 A JP 2024519990A
- Authority
- JP
- Japan
- Prior art keywords
- light beam
- sample
- input light
- srst
- imaging
- 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
- 238000001069 Raman spectroscopy Methods 0.000 title claims abstract description 26
- 238000003325 tomography Methods 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 title claims description 37
- 238000003384 imaging method Methods 0.000 claims description 94
- 239000000523 sample Substances 0.000 claims description 67
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 5
- 239000012472 biological sample Substances 0.000 claims description 4
- 238000007654 immersion Methods 0.000 claims description 2
- 239000011324 bead Substances 0.000 description 48
- 210000001519 tissue Anatomy 0.000 description 44
- 230000003287 optical effect Effects 0.000 description 17
- 229920003229 poly(methyl methacrylate) Polymers 0.000 description 17
- 239000004926 polymethyl methacrylate Substances 0.000 description 17
- 238000010009 beating Methods 0.000 description 15
- 238000000386 microscopy Methods 0.000 description 15
- 208000005017 glioblastoma Diseases 0.000 description 13
- 239000004793 Polystyrene Substances 0.000 description 11
- 210000004027 cell Anatomy 0.000 description 9
- 210000003763 chloroplast Anatomy 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 239000011543 agarose gel Substances 0.000 description 8
- 239000013078 crystal Substances 0.000 description 8
- 238000009826 distribution Methods 0.000 description 8
- 150000002632 lipids Chemical class 0.000 description 8
- 239000002904 solvent Substances 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 8
- 238000000701 chemical imaging Methods 0.000 description 7
- 230000035515 penetration Effects 0.000 description 7
- 241000196324 Embryophyta Species 0.000 description 6
- 230000035559 beat frequency Effects 0.000 description 6
- 238000005286 illumination Methods 0.000 description 6
- 241000234282 Allium Species 0.000 description 5
- 235000002732 Allium cepa var. cepa Nutrition 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 230000005284 excitation Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 229920002223 polystyrene Polymers 0.000 description 5
- 108090000623 proteins and genes Proteins 0.000 description 5
- 102000004169 proteins and genes Human genes 0.000 description 5
- 238000013135 deep learning Methods 0.000 description 4
- 239000000499 gel Substances 0.000 description 4
- 229920000642 polymer Polymers 0.000 description 4
- 240000001307 Myosotis scorpioides Species 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000008832 photodamage Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- HMQFJYLWNWIYKQ-UHFFFAOYSA-N 1,4-diphenylbutadiyne Chemical compound C1=CC=CC=C1C#CC#CC1=CC=CC=C1 HMQFJYLWNWIYKQ-UHFFFAOYSA-N 0.000 description 2
- 101100459910 Arabidopsis thaliana NCS1 gene Proteins 0.000 description 2
- 102000006386 Myelin Proteins Human genes 0.000 description 2
- 108010083674 Myelin Proteins Proteins 0.000 description 2
- 241001323321 Pluto Species 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 229930002875 chlorophyll Natural products 0.000 description 2
- 235000019804 chlorophyll Nutrition 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000002601 intratumoral effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 210000005012 myelin Anatomy 0.000 description 2
- 238000012014 optical coherence tomography Methods 0.000 description 2
- 238000012634 optical imaging Methods 0.000 description 2
- 235000015277 pork Nutrition 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000003174 Brain Neoplasms Diseases 0.000 description 1
- 241000282341 Mustela putorius furo Species 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 238000010870 STED microscopy Methods 0.000 description 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- 230000004323 axial length Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000032823 cell division Effects 0.000 description 1
- 230000012292 cell migration Effects 0.000 description 1
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 description 1
- 239000001752 chlorophylls and chlorophyllins Substances 0.000 description 1
- 238000001218 confocal laser scanning microscopy Methods 0.000 description 1
- 238000004624 confocal microscopy Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000003210 demyelinating effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000011331 genomic analysis Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000007102 metabolic function Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000000399 optical microscopy Methods 0.000 description 1
- 210000003463 organelle Anatomy 0.000 description 1
- 230000001991 pathophysiological effect Effects 0.000 description 1
- 230000029553 photosynthesis Effects 0.000 description 1
- 238000010672 photosynthesis Methods 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 210000001732 sebaceous gland Anatomy 0.000 description 1
- 238000012174 single-cell RNA sequencing Methods 0.000 description 1
- 150000003384 small molecules Chemical class 0.000 description 1
- 238000010869 super-resolution microscopy Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- 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/65—Raman scattering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/4795—Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/0004—Microscopes specially adapted for specific applications
- G02B21/002—Scanning microscopes
- G02B21/0024—Confocal scanning microscopes (CSOMs) or confocal "macroscopes"; Accessories which are not restricted to use with CSOMs, e.g. sample holders
- G02B21/008—Details of detection or image processing, including general computer control
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/33—Immersion oils, or microscope systems or objectives for use with immersion fluids
-
- 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/65—Raman scattering
- G01N2021/653—Coherent methods [CARS]
- G01N2021/655—Stimulated Raman
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/06—Illumination; Optics
- G01N2201/061—Sources
- G01N2201/06113—Coherent sources; lasers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/06—Illumination; Optics
- G01N2201/067—Electro-optic, magneto-optic, acousto-optic elements
- G01N2201/0675—SLM
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/06—Illumination; Optics
- G01N2201/069—Supply of sources
- G01N2201/0696—Pulsed
- G01N2201/0697—Pulsed lasers
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Microscoopes, Condenser (AREA)
Abstract
【要約】誘導ラマン散乱トモグラフィシステムは、第1の入力光ビームを生成する手段と、第2の入力光ビームを生成する手段と、対物レンズと、集光器と、検出器とを含む。第1の入力光ビームは位相変調され、第2の入力光ビームは振幅変調される。対物レンズは、第1および第2の入力光ビームを試料上に向けるように構成される。集光器は、試料からの出力光ビームを集光するように構成される。検出器は、第1の入力光ビームに対応する出力光ビームの少なくとも一部を検出するように構成される。システムは、前記出力光ビームの前記検出された部分から前記試料の深度分解画像を形成するための手段をさらに含む。
Description
本開示は、広義には、誘導ラマン散乱トモグラフィシステム及び方法に関するが、これに限定されない。
三次元(3D)光学顕微鏡法(例えば、共焦点顕微鏡法、多光子顕微鏡法、および高調波発生顕微鏡法)は、生体組織および細胞についての時空間情報(例えば、構造およびアーキテクチャ、代謝機能、ニューロンネットワーク、細胞分裂および遊走など)に関する重要なデータを送達するための強力な画像化ツールである。従来の3D体積データは、一連の2次元(2D)断面画像を取得するために対物レンズまたは試料ステージのいずれかを走査すること(例えば、ラスタ走査)を通して、z方向を機械的に横断する点毎または線走査を介して取得されることができる。従来の3D顕微鏡撮像技術の主な制限は、点走査に典型的に使用される密に集束されたガウスビームが、混濁媒体(例えば、組織)における屈折率不連続性に起因する強い光散乱に遭遇し、制限された光透過(約100~200μm)を被り、したがって、体積測定深部組織撮像が不可能であることである。散乱効果を軽減するために、その伝播経路に沿って散乱が生じた後に固有の自己再構成能力を有する非回折ベッセルビームが、高度な3D撮像技術(例えば、ライトシート顕微鏡法、2光子顕微鏡法、光コヒーレンス顕微鏡法)において回折及び散乱が支配的である混濁組織の深い領域への浸透を助けるために利用されてきた。しかしながら、ベッセルビーム顕微鏡法は、一般に、深さ情報が失われた試料の投影画像を提供する。
誘導ラマン散乱(SRS)顕微鏡法は、生体分子特異性を有する新興の無標識化学撮像技術であり、生物学的および生物医学的システムにおいて広範な用途が見出されている。従来のガウスビームSRS 3D撮像の撮像深度を改善するために、ベッセルビーム励起と結合された光学投影トモグラフィ(OPT)が、複数の投影角度における一連の2D画像から3D SRS画像を取り出すために使用されてきたが、OPTベースのSRS 3D撮像は、依然として、試料ステージまたは入射光ビームを機械的に回転させる必要があり、これは、迅速なインビボ生物学的/生物医学的用途には好適ではない。
従って、上記の問題の少なくとも一部に対処することができる断層撮影システム及び方法を提供する必要がある。
第1の態様によれば、第1の入力光ビームを生成する手段であって、第1の入力光ビームが位相変調される、生成する手段と、第2の入力光ビームを生成する手段であって、第2の入力光ビームが振幅変調される、生成する手段と、第1および第2の入力光ビームを試料上に向けるように構成された対物レンズと、試料からの出力光ビームを集光するように構成された集光器と、第1の入力光ビームに対応する出力光ビームの少なくとも一部を検出するように構成された検出器と、出力光ビームの検出された部分から試料の深さ分解画像を形成する手段とを備える誘導ラマン散乱トモグラフィシステムが提供される。
第2の態様によれば、位相変調された第1の入力光ビームを生成することと、振幅変調された第2の入力光ビームを生成することと、第1および第2の入力光ビームを試料上に向けることと、試料から出力光ビームを収集することと、第1の入力光ビームに対応する出力光ビームの少なくとも一部を検出することと、出力光ビームの検出された部分から試料の深さ分解画像を形成することとを含む誘導ラマン散乱トモグラフィ方法が提供される。
第2の態様の方法を含む3次元ボリューム撮像方法も開示される。
本発明の実施形態は、単なる例として、図面と併せて、以下の記載から、当業者にはよりよく理解され、容易に明らかになるであろう。
本開示は、非回折ゼロ次ベッセルビームに関連する光ビーティング技術(OBT)によって可能になるzスキャンフリーSRSトモグラフィ(SRST)に関し、試料または対物レンズの機械的走査を必要とせずに3D SRSのより深い化学撮像を実現する。SRSTでは、試料は、空間周波数領域において試料の深さ分解された化学情報を符号化する空間光変調器(SLM)を使用することによって生成される調整可能な光ビートパターンを有するベッセルビームによって照射される。3D SRS断層像は、逆高速フーリエ変換(IFFT)を実施することによって迅速に取り戻すことができる。本開示はまた、独特のベッセルビームOBTベースのSRST撮像方法の理論的な導出および分析を示し、また、ベッセルビームOBTベースのSRST撮像システムの開発を示して、様々な試料(例えば、ポリマービーズファントム、ラマン活性結晶、植物細胞、および生物組織)に対する無標識の体積測定によるより深い化学撮像のためのSRSTの能力を実験的に実証する。
図1は、例示的な実施形態によるSRST撮像システム100の概略図を示す。2つのレーザ出力(1041 nmの1つはストークスビーム102として使用され、他の調整可能な出力(680~1300 nm)はポンプビーム104として機能する)を有するフェムト秒(fs)レーザ源(Insight DS dual、Spectra-Physics-図示せず)が、SRST撮像における組織励起のために使用される。ポンプビーム104およびストークスビーム102は、80 MHzの繰り返し率で動作する約100 fsのパルス幅を有する。ポンプビーム104は、空間光変調器(SLM)106(PLUTO BB、Holoeye)によって位相変調される。ストークス光102は、アキシコン108(AX251-B、Thorlabs)を用いてベッセルビームに変換され電気光学変調器(EOM)(APE-ベルリン-図示せず)。ポンプビーム104とストークスビーム102は、ダイクロイックミラー110を用いて結合される。ポンプビーム104およびストークスビーム102を受け取るレンズL1(f=100 mm)の焦点面上のリングパターンは、4-fシステムを通して多光子走査顕微鏡(MPM-4R、Thorlabs)のガルボミラー上に結像され、次いで、さらに、レンズL1の背面開口上に投影される。水浸顕微鏡対物レンズ112(Apo LWD 25X、1.10w、Nikon)から、顕微鏡の別の内部4-fシステムを通して得た。順方向では、コンデンサ114(CC Achromat/Aplanat、N.A.=1.4、Nikon)を使用して、透過ポンプビームおよび透過ストークスビームを収集する。透過ポンプビームは、バンドパスフィルタセット116(Semrock)を使用して透過ストークスビームからスペクトル的に分離され(すなわち、透過ストークスビームは遮断される)、大面積フォトダイオード118(FDS1010、Thorlabs)によって検出される。ロックイン増幅器120(APE-Berlin)を使用して、フィルタリングされたポンプビームを復調し、試料からSRS信号(すなわち、誘導ラマン損失(SRL))を取得する。深度分解SRST画像は、例えば、逆高速フーリエ変換を実施することによって、SRS信号から迅速に取り戻すことができる。
SRS撮像では、試料は、共線ポンプビーム(ωp)、および高開口数(N.A.)対物レンズを通して厳密に焦点を合わせたもとでのストークスビーム(ωs)( ωs<ωp)によって照射される。生成されるSRS信号は、ポンプビーム強度とストークスビーム強度との積IpIsに比例する。ポンプ励起ビームとストークス励起ビームの両方が、試料を励起するベッセルビームである場合、それらの強度分布は、以下のように表すことができる。
ここで、ρは、ρ(pump)またはρ(Stokes)であり得る。Iρ(z)は、軸方向に沿った強度分布であり、J0(krρr)は、第1種のゼロ次ベッセル関数である。krρは、横方向波数ベクトルであり、rは、横方向座標である。ベッセルビームは非常に細いニードルビームであるので、その横方向の寸法を無視することによって、3D試料f(x,y,z)の2D投影SRS画像F(x,y)は、次のように表される。
光学ビーティング技術(OBT)は、ベッセルポンプビームが空間周波数Δkzを有するビーティングパターンを含むことを必要とする。すなわち、Ip(z)がIp(z)[cos(Δkzz)+1]dzとなり、次いで、2D投影F(x,y,Δkz)は、以下のように空間周波数領域における試料の深さ情報を反映する。
Δkz空間では、式(2)の余弦関数の周波数は、zに比例する。すなわち、より深い位置zにある物体は、Δkz空間でより速く振動する。その逆フーリエ変換は、以下のようにあらわされる。
ここで、Δkzは、Δkz,minからΔkz,maxまでに範囲であり、
であり、“*”は畳み込みを示す。方程式(3)の右辺の3つの項は、それぞれ、取り戻された画像f(x,y,z’)Ip(z’)Is(z’)、鏡像f(x,y,-z’)Ip(-z’)Is(-z’)、直流(DC)成分Cδ(z’)を表す。再構成された3D画像f(x,y,z’)は、照明関数Ip(z’)Is(z’)へ正規化し、しかし、鏡像とDC成分を省くことによって得ることができる。
軸方向分解能は、Δkz,max-Δkz,minのビート周波数範囲に反比例する。
図1を参照すると、SLM 106は、ベッセルビーム生成のためにポンプビーム104に位相パターンを課すために使用される。光学ビーティングパターンを有するベッセルポンプビーム104は、SLM 106の後に形成され得、顕微鏡対物レンズの後部開口上に現れる2つの同心リングをもたらす。最終的に、ベッセルポンプビームの光ビーティングパターンは、4-fシステム(L1および顕微鏡対物レンズ112)を通して試料上に投影される。対物レンズ112の後部開口上で、ベッセルストークスビーム102は単一のリングを形成する。これら2つのビームの積は、依然として、試料内にビーティングパターンを有する。
SRSTのためのOBTにおけるビーティングパターン生成のメカニズムを解明するために、図2aは、SLM 106上の位相パターンのいくつかの例を表示する。例示的な実施態様では、位相パターンは、異なる収束角を有する2つのアキシコン位相から構成される。顕微鏡対物レンズ112の後部開口上の対応する強度分布も示されている(図2bのリングを参照)。外側リングおよび内側リングの半径は、それぞれ、r1とr2である。各リングは、顕微鏡対物レンズを通過した後に、試料上にベッセルビームを生成する。2つの重畳されたベッセルビームは、異なる軸方向の波数ベクトル、すなわち、
および
を有し、ここで、λpは、ポンプ波長であり、nは、サンプルの屈折率であり、
および
であり、Rは、対物レンズの後部開口の半径である。それらの干渉は、試料内のベッセルビームに沿って光ビーティングパターンをもたらす。ビート周波数Δkzは、
である。
式(5)によれば、ビート周波数Δkzはr1またはr2を変化させることにより可変である。ビーティングパターンを有する対応するビーティングベッセルビームが図2cに示されている。
周波数領域光コヒーレンストモグラフィ(FD-OCT)と同様に、式(3)は、取り戻された画像、鏡像、およびDC成分を含む。FD-OCTでは、これらが重なり合って、取得された画像が判別できなくなることがある。しかし、SRSTは、そのような重複の問題がないという利点を有する。レンズL1および対物レンズは、4-fシステム(図1参照)を形成し、したがって、サンプル上のビームは、SLM 106の後のビームの像である。式(2)によれば、試料上のz=0平面は、任意のΔkzに対して位相φ=Δkzzが0の平面である。SLM 106上の位相パターンの中心画素の位相は、任意のΔkzに対してゼロであるので、料上のz=0平面は、ちょうどSLM106の像平面である。ベッセルビームは、SLM106の後に形成され、したがって、z>0のときのみ照明関数Ip(z)>0である。したがって、例示的な実施形態によるSRSTで取り戻された画像は、鏡像およびDC成分と区別することができる。さらに、インターフェログラムを生成するために試料の外側に参照アームを必要とし、その撮像を反射モードのみに制限する(ほとんどの生体撮像に適用可能であるが)OCTとは異なり、例示的な実施形態によるOBTベースのSRSTは、インターフェログラムがSRSTにおける3D組織撮像のために試料の内側で直接生成されるので、外部参照アームを必要とせずに透過モードと反射モードの両方で生体撮像に適用され得る汎用性がある。
(i)Δkz空間における撮像範囲およびステップサイズ:例示的な実施形態によるSRST撮像では、141個のフレームの位相パターンがSLM106上に1つずつ表示され、これはステップサイズ0.0049μm-1で0から0.68μm-1までの範囲でのベッセルポンプビーム104の141個のビート周波数Δkzを生成する。141個の生SRS画像が、画像スタックを形成するために、各位相パターンに対して連続的に取得される。次いで、3D深度分解SRST画像が、逆高速フーリエ変換(IFFT)によってSRS原画像のスタックから迅速に取り出される。
(ii)試料内のベッセルビームの長さ: 例示的なSRSTシステムでは、ポンプビームとストークスベッセルビームの両方の長さは、約200μmである。
(iii) レーザ出力: レーザ出力:試料上の平均レーザ出力は、典型的には、例示的なSRSTシステムに対してPSRST=100mWであり、従来のSRSに対してPcon-SRS=10mWである一方、対応するピーク電力密度は、SRSTに対して
であり、従来のSRSに対して
である。ここで、f=80MHzはレーザ繰り返し率であり、τ=100fsはパルス幅であり、A=π(d/2)2=0.7μm2は焦点スポット面積であり(半径dは
である)、C=1/40はガウスビームとベッセルビームとの間の軸方向に沿った長さの比であるスケーリング係数である。ベッセルビームのレーザパワーは軸方向に沿って広がるため、ピークパワー密度はSRST撮像では比較的低い。
撮像時間:一例では、典型的な画素滞留時間は、SRSTに対してΔτSRST=45μsであり、従来のSRSに対してΔτcon-SRS=11μsである。そのようなピクセル滞留時間は、試料におけるSRST撮像の良好な品質のための適切な信号対雑音比(SNR)レベルを維持するために選定される。画素数512×512に対して、2Dラスタ走査速度は、SRST撮像に対して12sであり、従来のSRS撮像に対して3sである。使用されるSLM106(PLUTO BB、Holoeye)のフレームレートは約60Hzであり、顕微鏡のzモータ走査は50Hzである。
例示的な撮像用途で使用される化学的および生物学的試料は、ポリメチルメタクリレート(PMMA)およびポリスチレン(PS)ビーズ、1,4-ジフェニルブタ-1,3-ジイン(DPBD)結晶、春玉葱、ブタ組織、およびマウスの耳を含む。
撮像用の第1の試料は、ポリメチルメタクリレート(PMMA)およびポリスチレン(PS)を含む。この実施例では、10μmポリメチルメタクリレート(PMMA)およびポリスチレン(PS)ビーズの混合物がSRST撮像用の2wt%硬化アガロースゲルファントムに包埋する。設計されたビーズゲルファントムの寸法は2×2×3cm3であり、ビーズファントム内のポンプ(800nm)およびストークス(1040nm)レーザビームの散乱平均自由経路は、それぞれ65μmおよび55μmである。PMMAビーズは2950cm-1においてシグナルを与える(O-CH3中のC-Hの対称振動およびCH2の非対称振動)一方で、PSビーズは2950cm-1および3050cm-1の両方のラマンシフトにおいてSRSシグナルを生成する。2950cm-1および3050cm-1におけるPSのSRS強度比は約0.8:1である。アガロースゲルのSRSシグナルは、使用されるその低濃度(約2重量%)のために、これらの2つのラマンシフトにおいてPSおよびPMMAよりもはるかに弱い。PMMAビーズのSRS画像は、2950cm-1~3050cm-1におけるSRS画像を差し引くことによって得られ、これらはSRS強度比へ正規化される。
例えば、上述のようなポリマービーズアガロースゲルにおけるポンプ光ビームおよびストークス光ビームの散乱平均自由行路は、ミー散乱モデルに基づいて推定することができる。ファントム内の混合ビーズの濃度Nsphereは約105sphere/mm3である。800nm(ポンプビーム)におけるPMMA及びPSの屈折率は、それぞれ、1.48及び1.58である([https://refractiveindex.info/]). ファントム内のアガロースゲルの屈折率は約1.34である。ミー散乱モデル化を単純化するために、全てのビーズの平均屈折率は約1.53である([https://refractiveindex.info/])。ファントム内のアガロースゲルの屈折率は約1.34である。Mie散乱モデル化を単純化するために、全てのビーズの平均屈折率は約1.53であると仮定する。
上記のパラメータをMie散乱計算モデル[https://omlc.org/calc/mie_calc.html]に与えると、得られたビーズのMie散乱断面積σは157.77μm2であり、散乱係数μsはσNsphere=157.77cm-1である。したがって、ビーズゲルファントムにおける800nmでのポンプビームの散乱平均自由経路
は、約65μmである。同様に、ビーズゲルファントムにおける1040nm(ストークス波長)において散乱係数μsは181.61.cm-1であり、したがってビーズゲルファントムにおける1040nmでのストークスビームの散乱平均自由経路
は約55μmである。
図3は、アガロースゲルファントムのマトリックスに埋め込まれた混合された10μmのポリスチレン(PS)及びポリメチルメタクリレート(PMMA)ビーズのSRSTに基づく例示的なSRSTシステムの性能評価を示す。混合された10μmのポリスチレン(PS)及びポリメチルメタクリレート(PMMA)ビーズを上記のようにアガロースゲルファントムに包埋したものの化学選択的SRST画像および従来の点スキャンガウシアンベームSRS画像を、それぞれ図3aおよび3bに示す。ファントム内のビーズのz深さは、疑似カラースケールを使用することによって示される。例示的な実施形態によるSRSTシステムは、異なるファントム深さにわたって異なるポリマービーズに関する深さ分解化学情報を正確に識別し、これは、従来のポイントスキャンSRS撮像とほぼ同一である。ビーズの深さ情報がSRSTにおいてどのように決定されるかを示すために、一例として、図3aのより深いビーズ(ビーズ1)およびより浅いビーズ(ビーズ2)が観察のために選択される。図2cは、これら2つのビーズを横切って引かれた線に沿った空間周波数領域における測定データを示す。これは、より深いビーズはΔkz軸に沿ってより速く振動することを明確に示し、これは、式(2)の予測と一致する。図3dは、x-z平面上のビーズ1およびビーズ2の取り戻されたSRST画像を示し、zスキャンフリー3d撮像のためのSRSTの光学切片化能力を確認する。
より深い3D化学撮像のための例示的な実施形態によるSRST技術の利点を実証するために、SRST及び従来のSRSにおけるビーズのクラスタが、それぞれ図3e及び図3fに表示される。このクラスタの最上層は約z=21μmである。z=32μmにおけるSRST画像とSRS画像を比較すると、はるかに強いSRS信号が提供される。より深い深さz=46μmでは、ビーズは、SRSTにおいて依然として明確に観察することができ、より深い領域におけるビーズのSRST画像輝度レベルは、z=21μmでのより浅いビーズのそれよりわずかに薄暗いのみである。対照的に、z=46μmでのビーズの従来のガウスビームSRS信号レベルは、より浅いビーズのものよりもはるかに弱い(図3fを参照)。結果は、SRST撮像で使用されるベッセルビームが、散乱障害物を迂回した後に自己再構築され、それにより、例示的な実施形態によるSRST技術が、試料内の他のより浅いビーズの影の下にあるより深いビーズにより良好に到達することが可能になり、より深い組織撮像がもたらされることを確認する。
ベッセルビームOBTベースのSRSTによる撮像深度の改善を定量的に分析するために、図3gは、従来のガウスビームSRS顕微鏡法(線304)と比較した、SRSTにおける異なる深度でのビーズの正規化されたSRS強度(線302)を示す。プロットは、指数関数的減衰関数
によって適合される。ここで、z0は、浸透深さである。SRSTについてはz0=130μm(95%信頼区間:[77μm, 183μm])である一方、従来のSRSについてはz0=58μm(95%信頼区間:[42μm, 73μm])である。例示的な実施形態によるSRST技術は、従来のSRS顕微鏡法と比較して、透過深度において2倍を超える改善を提供することができ、より深い組織の3D撮像のためのSRSTの可能性をさらに確認することができることが分かる。この増強効果は、図3aを図3bと比較することによっても明確に観察することができ、SRSTにおいてより深い深さに位置するビーズ(例えばビーズ1)は、従来のSRS画像におけるものよりもはるかに明るく見え、試料におけるガウスビームよりもベッセルビームのより深い光透過を証明する。
例示的な実施形態によるベッセルビームOBTベースのSRSTシステムは、生体試料における深部組織撮像にも使用することができる。例えば、図4aおよび4bは、異なる深度でのマウスの耳の皮膚のSRSTおよび従来のSRS画像(タンパク質および脂質のCH3ストレッチの2935 cm-1でのラマンシフト)の比較を示す。皮膚表面の近く(z=20μm)ではSRSTおよびSRS撮像の両方で脂質の豊富な皮脂腺がはっきりと見られたが、組織のより深くに進むと、SRSTは、従来のSRS撮像と比較して、はるかに強いSRS信号を示す(例えば、図4aおよび4bにおける矢印によってマークされるそれらの場所)。図4cは、図4aおよび図4bの選択された領域(楕円)における強度を比較する。異なる深さでのSRST画像からの測定強度は、表面に近い(z=20μm)強度に正規化される。従来のSRS画像からの強度は、同じ方法で正規化される。図4cは、SRSTの強度が、組織深度の増加に伴って、従来のSRSのものよりはるかにゆっくりと減衰することを明確に示し、SRSTが、より深い組織撮像のためのはるかに優れた能力を有することを確認する。図4dおよび図4eは、SRSTおよび従来のSRSの3D画像を比較し、SRSTおよび従来のSRSの両方が、皮膚表面により近い組織領域において非常に類似した画像を与えるが、SRSTは、より深い組織深度においてはるかに鮮明な画像を提供することを反映している。
例示的な実施形態によるSRSTシステムは、様々な撮像ターゲット(例えば、生体組織、ラマン活性結晶、及び植物細胞)上の3D体積測定深部分子撮像にも使用される(図5)。より良好な視覚化のために、図5a~図5cに示されるように、深度情報は、疑似カラースケールを使用することによって提示される。
図5aは、1,4-ジフェニルブタ-1,3-ジイン(DPBD)結晶
のSRST 3D画像を示す。DPBDは、生体試料中の小分子の追跡及び細胞小器官の標的化に広く適用される有用なラマンタグである。図5dは、例示的な実施形態によるシステムに基づくSRST画像を、2つの異なる深度における従来のポイントスキャンSRS画像と比較する。より浅い深さ(例えば、z=40μm)での結晶子は、SRSTと従来のSRS撮像との間で同様のSRS信号レベルを生成する一方で、より深い結晶子(矢印でマークされる)(例えばz=133μm)は、SRSTにおいて明確に見出すことができるが、従来のSRSにおいてはほとんど不可視である。SRSTは、はるかに良好な撮像深度を示す。
図5bは、春玉葱の葉緑体の2935cm-1におけるSRST 3D画像(2935cm-1におけるSRST信号は以下に記載するようにクロロフィルにおけるC-H伸縮および2光子吸収の両方を含む)を示す。葉緑体の大部分は植物組織表面のより近くに分布し、それによって、暗い環境下でさえも有効な光合成を可能にする。図5eは、例示的な実施形態によるシステムに基づくSRST画像を、図5bのボックス内で強調表示された選択された領域の従来のSRS画像と比較する。いくつかの葉緑体(矢印でマークされる)は、SRST(例えば、z=45μm) において明確に観察され得る。これは、おそらく、SRS撮像におけるガウスビーム照明下では、それらのより深い葉緑体が、いくつかのより浅い葉緑体の影の下にある(例えば、z=27μm)という理由による。一方、本発明の実施形態によるSRSTシステムで使用されるベッセルビームは、植物細胞内の浅い葉緑体をバイパスすることによって、深い葉緑体に到達することができる。
本実施例では、スペクトルフォーカシングSRSシステムを用いて、春玉葱葉のクロロプラストのSRSスペクトルを測定する。各測定は、166cm-1をカバーすることができる。したがって、ポンプビームの中心波長は、4回シフトされて、2786~3045cm-1のスペクトル範囲全体をカバーする。2900~2950cm-1におけるラマンピークは4つの測定の全てにおいて明確に観察され(図6を参照)、CH3振動に対応する。CH3基は、異なるタイプのクロロフィルにおいて一般的に観察される。加えて、強い非化学的に特異的な2光子吸収バックグラウンド(ラマンピーク強度の約70%)も観察される。3000cm-1を超えて、水からのSRS信号により信号が増加する。
図5cは、豚肉組織中の脂肪の(脂質のCH2の2845cm-1における)SRST 3D撮像の別の例を示す。細胞および組織生理機能(例えば、エネルギー貯蔵)において重要な役割を果たす脂質液滴の3D形態および分布は、SRST撮像において、より深い組織領域においてさえも明確に観察することができる。図5fは、例示的な実施形態によるシステムに基づくSRST画像を、図5cにおいてボックスによってラベル付けされた選択されたエリアの従来のSRS画像と比較する。例えば、15μmにおいて矢印でマークされた3つの小さな脂肪滴がある。ガウスビーム照明を用いた従来のSRS撮像では、光の影が、より深い脂肪滴(例えば、z=15μm) 内にアーチファクト(矢印によってもマークされる)を生成するが、組織におけるベッセルビーム伝播の自己再構成特性の利点のために、SRST撮像においてはそのようなアーチファクトは観察されない。上記の結果は、例示的な実施形態によるベッセルビームOBTベースのSRSTが、生物学的および生物医学的システムにおける体積測定のより深い組織撮像のための強力な無標識3D分子撮像ツールであることをさらに実証する。
例示的な実施形態によるSRSTと従来のSRS撮像との間の軸方向分解能及び横方向分解能も比較される。SRSTでは、深度分解能は、式(4)によって決定され、Δkzは、式(5)によって決定される。本発明の実施形態に係るSRST撮像において、r1は固定され、
に対応し、r2は0.5r1からr1まで変化し、0.27から0.54までの
に対応する。ここで、n=1.33が、使用される組織の屈折率(水の屈折率に近い)である。したがって、λp=800nmに対してΔkzの範囲は0~0.6μmである。F(x,y,Δkz)はΔkzの偶関数あるから、撮像取り戻しのためのΔkzの有効ビート周波数はー0.68~0.68μm―1まで変化し、これはSRST撮像における5.6μmの深さ分解能に対応する。この予測されたSRST深さ分解能は、測定された深さ分解能5.49μmと一致する(以下の計算参照)。N.A.=0.54を使用した従来の点走査SRSとの比較において、軸方向に沿ったエアリーディスクの半値全幅(FWHM)(2λ/N.A.2)は、λp=800nmおよびλs=1041nmに対して、それぞれ、5.5μmおよび7.1μmである。深さ分解能は、従来のSRS撮像では、
である。一方、横方向の非回折ベッセルビームサイズは、試料中のガウシアンビーム(λp=800nm、N.A.=0.54)の横方向分解能(約0.90μm)の約1.33分の1に絞ることができる。したがって、例示的な実施形態によるSRST撮像は、従来のSRSと同様の軸方向分解能を共有するが、細胞内分解能を有する無標識3D容積測定深部組織撮像を達成するためのより高い横方向分解能を有する。
図7aは、x-z平面における、水中に浸漬されたDPBD結晶のSRST画像を示す。図7bは、図7aに示される白線に沿った結晶と水との界面にわたる強度プロファイルを示す。データは、誤差関数
を使用してフィッティングされる。ここで、A1、A2、およびz0は、フィッティングパラメータである。フィッティング結果は、曲線702としてプロットされる。軸方向分解能は、フィッティング関数h(x)内のガウス関数のFWHMである。すなわち、強度プロファイルに示され、上で述べたように、分解能は
である。
例示的な実施形態によるepi検出SRST技術は、腫瘍内不均一性の迅速な無標識分子評価、及びサブミクロン分解能での全膠芽腫(GBM)組織標本における分子サブタイプ化にも適用することができる。図8 A~8 Cは、例示的な実施形態によるSRSTシステムを使用して、ほんの一瞬の間に2850cm-1および2940cm―1のラマンシフトで取得された3つのGBMサブタイプ(それぞれ、古典的、間葉系および前神経)の代表的なSRS画像を示す。各ラマンシフトは、2850 cm-1で、GBM中の脂質分布が脂質分子のCH2結合の対称伸縮により可視化される、明確な生体分子分布を提供し(図8の第1列)、一方、脂質およびタンパク質の両方におけるCH3結合の伸縮による2940 cm-1のラマンシフトへの共鳴は、組織の均一な輝度を生じる(第2列)。タンパク質に対する特異性を高めるために、2850 cm-1でのSRS画像を2940 cm-1から減算して、タンパク質分布を強調する(第3列)。さらに、2850 cm-1および減算画像を重ね合わせて、ミエリンを含む細胞形態を示す(最後のカラム)。一般に、ミエリンの濃度および統合性は、前神経-古典的-間葉の順に減少していることが観察され得、これは、迅速なGBM分子サブタイプ化のためのSRST画像化技術の染色のない組織学的可能性を示唆する。
SRS撮像の結果は、染色なしのSRST組織学的画像および2Dサブタイプマップの両方が20~30分以内に得られることを示す。代替実施形態では、例えば、より高速のスキャナを使用することによって、かかる時間をさらに短縮することができることが理解されるであろう。このような性能は、従来の単一細胞RNA配列決定よりも優れている。SRS組織学的結果は、脱髄状態を新しい診断特徴として評価するが、不均一性マッピングは、腫瘍内不均一性への新規の洞察を明らかにする。GBM組織の大部分はゲノム解析の診断結果と一致するが、標本中の残りの画像タイルの重要な部分は他の分子サブタイプに属することが見出され、無標識SRS撮像によって明らかにされたGBM不均一性の実質的な程度を示唆する。
さらなる実施態様では、各GBM組織上のハイパースペクトルSRS撮像から取り出されたSRSスペクトルは、GBM分子サブタイプ化のためのディープラーニングアルゴリズムを用いて調査される。組織画像全体の各画像タイルは、サブタイプ間の比較のために512×512ピクセルを平均することによって単一のスペクトルを生成する。一般に、MES-CL-PNの順にCH2およびCH3伸縮ピーク(2850および2940 cm-1)の減少がある。これに対して、GBMサブタイプのタンパク質スペクトルは互いに類似しており、比較的低い診断情報が埋め込まれていることを示唆している。SRS撮像からのSRS分光情報がGBMサブタイピングのための有意な診断性能をもたらすことができるかどうかを検証するために、SRSスペクトルが二次サポートベクターマシンモデルに提供される。訓練セットおよび検証セットについて80/20の比を有するハンドアウト検証は、混同行列(図9A)に示されるように、80.6%の診断精度を与える。ディープラーニング診断モデルのロバスト性は、約0.89の平均面積積分を有する受信者動作特性(ROC)曲線によってさらに確認される(図9C~9E)。したがって、この実装形態で開発されたディープラーニングアルゴリズムは、無染色脳腫瘍診断のためのSD-SRST 3D撮像、および精密神経手術のための腫瘍マージン画定にさらに拡張され得る。
図10は、例示的な実施形態による誘導ラマン散乱トモグラフィ方法のフローチャートを示す。ステップ1002において、第1の入力光ビームが生成される。第1の入力光ビームは位相変調される。ステップ1004では、第2の入力光ビームが生成される。第2の入力光ビームは振幅変調される。ステップ1006では、第1および第2の入力光ビームが試料上に方向付けられる。ステップ1008では、試料からの出力光ビームが収集される。ステップ1010では、第1の入力光ビームに対応する出力光ビームの少なくとも一部が検出される。ステップ1012では、出力光ビームの検出された部分から試料の深度分解画像が形成される。
上述のように、従来の3D体積撮像における課題、例えば、(i)生細胞及び組織における病態生理学的環境を乱す可能性がある蛍光標識(共焦点蛍光顕微鏡法)の必要性、及び(ii)混濁組織における集束ガウスビーム照明の散乱によるより短い撮像深度は、生物系における迅速な動的及び機能的3D撮像のためのそれらの広範な適用を妨げてきた。本発明の実施形態によるベッセルビーム-光ビーティング技法(OBT)に基づくz-スキャンフリー誘導ラマン散乱トモグラフィ(SRST)は、上記の問題を解決して、無標識の体積測定によるより深い組織の撮像を達成することができる。例示的な実施形態では、組織深さにわたってレーザ集束の機械的走査を必要とせずに、SLMを通してポンプビームに投影される位相パターンを電子的に変化させることによって、試料内の2つの重なり合うベッセルビームの干渉により、光学切片化のためのポンプビームの光学ビーティングパターンを迅速に生成することができる。試料内のビーティングベッセルポンプビームとベッセルストークスビームとの重ね合わせは、3D SRS撮像のための深度符号化SRSを生成する。したがって、細胞内分解能を有する体積測定試料についての深度分解3D化学分布は、IFFTを実施することによって迅速に取り出すことができる。
上述の例は、本ベッセルビームOBTベースのSRST技術の有用性を実証しており、アガロースゲルに埋め込まれた多数の試料(例えば、ポリマービーズ(PMMAおよびポリスチレン)、ラマン活性結晶、植物細胞、および生物組織(マウスの耳の皮膚、ブタの組織))における体積測定のより深い化学撮像を可能にする。3D SRS撮像深度を改善するために、ポンプビーム及びストークスビームとしての非回折ゼロ次ベッセルビームが、OBTベースのSRST撮像に用いられ、これは、光路に沿って遭遇する障害物を超えて自己再構成特性を示し、混濁媒体における散乱効果に対して顕著な復元力を有する。ベッセルビーム伝搬の固有の自己再構成能力は、SRSTによるより深い3D化学撮像について実証されており、これは、上述の例において明確に示されている。これらの例は、障害物の後のベッセルビーム伝播の自己回復を確認し、それによって障害物の影の下の分子を照射し、SRS撮像における影のアーチファクトを排除し、散乱弾性特性を有する高散乱媒体中を伝播するときのベッセルビームによるエネルギー損失が比較的低いことを確認する。結果として、例示的な実施形態によるベッセルビームOBTベースのSRST撮像は、より深い生物医学的撮像における進歩をもたらすことができ、混濁媒体における回折及び散乱は、従来のガウスビームSRS顕微鏡法を使用することによるタイトな焦点の形成を妨げる。
SRSTは、本質的に、多重化検出技術(すなわち、生画像データが、層ごとの走査ではなく、投影画像として取得される)であり、ノイズが試料とは無関係であるという条件で、フェレットの利点により、取り戻された画像のSNRを改善することができることが理解されるであろう。実際に、SRSは、レーザ源の固有の特性である励起レーザビームのショットノイズによって制限されるので、例示的な実施形態によるSRSのノイズは、試料とは無関係である。図3に示される例では、SNRは、生データから取り出されたSRSTまで5.3倍改善される。取り戻されたSRST画像の最終的なSNRは、SRST撮像における試料中の溶媒バックグラウンドによって影響されない。
例示として、図11は、SRSTにおける生データの信号対雑音比(SNR)、取り戻されたSRST画像(図3a)、およびビーズの従来のSRS(図3b)を比較する。これは、Fellgettの利点により、SNRが、生データ(SNR=3)から取り戻されたSRST(SNR=16)まで約5.3倍改善されることを示す。取り戻されたSRSTのSNRは、従来のSRS撮像のSNR(SNR=21)に匹敵するが、SRSTでは、4倍低いピーク電力密度があり、16倍弱い局所信号をもたらし、10倍高い平均電力があり、
高いノイズをもたらす。
溶媒バックグラウンド信号がSRST撮像のSNRを劣化させないことを示すために、SRSシミュレーション結果を図12に示す。具体的には、シミュレーションで使用される2つの異なる試料が、それぞれ図12a~図12bに示されている。試料Iは溶媒なしでz=100μmにおけるビーズであり、試料IIは溶媒中の同じビーズである(溶媒40~160μmの範囲で、溶媒バックグラウンドシグナルレベルはビーズの60%である)。SRSTにおける生データ(投影画像)は、それぞれ図12c~図12dにプロットされているように、Δkz空間(フーリエ領域)における両方の試料のSRS信号を記録する。図12 e~図12fは、それぞれ図12c~図12dにプロットされた曲線の逆高速フーリエ変換(IFFT)によって得られる、z軸に沿った取得されたSRST画像の強度プロファイルを示す。図12c~図12fにおいて、曲線の一方のセットは、純粋なビーズSRS信号を表し、曲線の他方のセットは、純粋なビーズSRS信号+ノイズを表す。
示された生データ中のノイズは、試料とは無関係であるべきである。いくつかの実装形態では、図12cおよび図12dにおいて同じポアソンノイズ(7×10―4nJ/pixel)が加えられている。信号レベルは図12d(約7×10―4nJ/pixel)における方が図12c(約5×10―4nJ/pixel)における方よりはるかに高い。なぜならば、溶媒は投影画像においてはるかに強いバックグラウンドシグナルを生成するからである。最後に、取り戻されたSRST画像(図12e~図12f)のSNRは、ほぼ同じであり、強いSRSバックグラウンド信号を有する溶媒がSRST撮像の最終SNRに影響を及ぼさないことを証明している。
SRSTにおけるベッセルビーム照明下では、レーザパワーは、軸方向に沿って伸びる((約100μmの例示的な実施形態)に沿って広がることが理解されるであろう。したがって、SRSTは、従来のガウスSRS撮像と比較して、試料上で同じ局所出力密度を達成するために、はるかに高い総レーザ出力を有する必要がある。光損傷は、主に、試料上に与えられる局所的な電力密度に依存するので、SRST撮像においてより高い総電力を使用することは、試料を損傷しない。本SRSTの例では、SRST(100mW)における試料上の総レーザパワーは、試料に光損傷を引き起こすことなく、従来のSRS(10mW)のそれより高い。一方、従来のSRSでは、総パワーが約30mWを超える場合、明らかな光損傷が観察される。、
SRSTは、修正された撮像モード、いわゆる周波数領域SRST(FD-SRST)において撮像速度を著しく改善することができることが期待される。これは、以下のように説明することができる。式(5)によれば、複数の波長成分Δλpを有する広帯域レーザ源λpは、複数のΔkzを生成するために使用することができる。したがって、FD-SRSTでは、SRS信号は、スペクトルが各ピクセル内に記録される分光計によって検出されることができ、深度情報は、IFFTを用いて取得されたスペクトルを通して容易に明らかにされることができる。3D画像は、2Dラスタスキャンを1回だけ実行することによって迅速に取得することができる。したがって、3D撮像速度は、2D撮像と同じくらい速くなり得る。このような実施形態では、ポンプビームとストークスベッセルビームの両方が広帯域であるべきである。同時に、両方のビームは、高速3D SRS撮像のためにスペクトル集束技術が広帯域励起の下でさえ特定のラマンピークを選択するために適用され得るように、チャープされるべきである。
例示的な実施形態によるベッセルビーム-OBT法のzスキャンフリー光学切片化特性は、SRST 3D撮像のみに固有ではなく、普遍的であることが理解されるであろう。記載されるOBT法は、高速3D組織撮像のための実質的に任意の他の非線形光学撮像モダリティに容易に適合され得る。例えば、現在のOBTベースのSRSTシステムは、フォトダイオードを光電子増倍管(PMT)に置き換えて、より深い組織領域から3D CARS信号を収集することによって、コヒーレント反ストークスラマン散乱(CARS)トモグラフィの準備ができている。SRST撮像システムにおいてポンプビームまたはストークスビームのみを使用する場合、OBT-SRST技術は、第2/第3高調波発生(SHG/THG)トモグラフィ、蛍光トモグラフィ、および多光子トモグラフィに簡略化することができる。さらに、SRST撮像におけるベッセルビーム-OBT法は、高速超解像3D深部組織撮像を実現するために、誘導放出抑制(STED)顕微鏡法、飽和誘導ラマン散乱顕微鏡法、および高次コヒーレントラマン散乱顕微鏡法(高N.A. 対物レンズが使用される場合)などの超解像顕微鏡法技法とも適合する。
要約すると、例示的な実施形態は、細胞内分解能を有するラベルフリーの体積測定化学撮像においてより深い浸透を達成するために、自己再構成ベッセルビームと結合された光ビーティング技術(OBT)を使用することによって可能にされるzスキャンフリーの誘導ラマン散乱トモグラフィ(SRST)を提供する。軸方向に沿った機械的走査を必要とせずに、体積組織に関する深さ分解SRS信号は、位相変調器を用いてベッセルポンプビームの光ビート周波数を電子的に調整することによって空間周波数領域において符号化され、したがって、深さ分解SRSTは、3D SRS撮像のためにIFFTを実施することによって取り出すことができる。例示的な実施形態はまた、ベッセルビームOBTベースのSRST撮像が、従来のポイントスキャンガウスビームSRS顕微鏡法と比較して、高散乱媒体における撮像深度の少なくとも2倍の改善を提供することを示す。無標識の容積測定のより深い分子撮像のためのSRST技術の有用性は、様々な試料(例えば、ラマン活性結晶、生体組織、および植物細胞)に対して実証され、これは、透過深度の点で従来のSRS顕微鏡法よりも優れている。SRSTにおけるベッセルビーム-OBT法のzスキャンフリー光学切片化能力の一般性は、生物学的および生物医学的システムにおける深部組織体積3D撮像のための実質的に任意の他の非線形光学撮像モダリティに容易に拡張することができる。したがって、zスキャンフリー光学切片化のためのSRSTにおける強力なベッセルビーム-OBT法は、高度な3D顕微鏡撮像用途の分野全体に対して著しい影響を有し得る。
当業者であれば、広範に記載された本発明の範囲から逸脱することなく、特定の実施形態に示された本発明に対して多数の変形および/または修正を行うことができることを理解するであろう。したがって、本実施形態は、あらゆる点で例示的であり、限定的ではないと考えられるべきである。
Claims (25)
- 誘導ラマン散乱トモグラフィシステムであって、
第1の入力光ビームを生成する手段であって、前記第1の入力光ビームが位相変調される、生成する手段と、
第2の入力光ビームを生成する手段であって、前記第2の入力光ビームが振幅変調される、生成する手段と、
前記第1および第2の入力光ビームを試料上に向けるように構成された対物レンズと、
前記試料からの出力光ビームを集光するように構成された集光器と、
前記第1の入力光ビームに対応する前記出力光ビームの少なくとも一部を検出するように構成された検出器と、
前記出力光ビームの検出された部分から前記試料の深さ分解画像を形成する手段と、を備える、誘導ラマン散乱トモグラフィシステム。 - 前記第1の入力光ビームを生成する手段は、レーザ光源の第1の出力を含み、前記第2の入力光ビームを生成する手段は、前記レーザ光源の第2の出力を含む、請求項1に記載のシステム。
- 前記レーザ光源は、広帯域フェムト秒レーザ光源を含む、請求項2に記載のシステム。
- 前記第1の入力光ビームはポンプビームを含み、前記第2の入力光ビームはストークスビームを含む、請求項1に記載のシステム。
- 前記第1の入力光ビームを生成する手段は、所定のパターンに基づく前記ポンプビームを位相変調するための空間的光変調器をさらに含む、請求項4に記載のシステム。
- 前記第2の入力光ビームを生成する手段は、所定の周波数において前記ストークスビームを振幅変調するための電気光学的変調器をさらに含む、請求項4に記載のシステム。
- 前記ポンプビームとストークスビームは、共線ベッセルビームを含む、請求項4に記載のシステム。
- 前記対物レンズは、水浸顕微鏡対物レンズを含む、請求項1に記載のシステム。
- コンデンサの後に配置され、前記第1の入力光ビームに対応する前記出力光ビームの前記部分をスペクトル的に分離するように構成された帯域通過フィルタセットをさらに含む、請求項1に記載のシステム。、
- 前記検出器は、フォトダイオードを含む、請求項1に記載のシステム。
- 前記試料の前記深さ分解画像を形成するための前記手段は、前記試料からの前記出力光ビームの前記検出された部分を復調するためのロックイン増幅器を含む、請求項1に記載のシステム。
- 前記ロックイン増幅器は、逆高速フーリエ変換に基づいて前記出力光ビームの前記検出された部分を復調するように構成される、請求項11に記載のシステム。
- 前記出力光ビームは、前記試料からの反射ビームを含む、請求項1に記載のシステム。
- 前記出力光ビームは、前記試料からの透過ビームを含む、請求項1に記載のシステム。
- 誘導ラマン散乱トモグラフィ方法であって、
第1の入力光ビームを生成することであって、前記第1の入力光ビームは位相変調される、生成することと、
第2の入力光ビームを生成することであって、前記第2の入力光ビームは振幅変調される、生成することと、
前記第1および第2の入力光ビームを試料上に向けることと、
前記試料から出力光ビームを収集することと、
前記第1の入力光ビームに対応する前記出力光ビームの少なくとも一部分を検出することと、
前記出力光の前記検出された一部分から前記試料の深さ分解画像を形成することと、を含む、方法。 - 前記第1の入力光ビームはポンプビームを含み、前記第2の入力光ビームはストークスビームを含む、請求項15に記載の方法。
- 前記ポンプビームは、空間光変調器を用いて所定のパターンに基づいて位相変調される、請求項16に記載の方法。
- 前記ストークスビームは、電気光学変調器を用いて所定の周波数で振幅変調される、請求項16に記載の方法。
- 前記ポンプビーム及びストークスビームを生成することは、共線ベッセルビームを形成することを含む、請求項16に記載の方法。
- 前記試料の前記深さ分解画像を形成することは、ロックイン増幅器を使用して、前記試料からの前記出力光ビームの前記検出された部分を復調することを含む、請求項15に記載の方法。
- 前記出力光ビームの前記検出された部分を復調することは、逆高速フーリエ変換を適用することを含む、請求項20に記載の方法。
- 請求項15に記載の方法を含む、3次元体積撮像方法。
- 前記試料は、化学試料又は生物試料を含む、請求項22に記載の3次元体積撮像方法。
- 前記試料は、ラベルフリーである、請求項22に記載の3次元体積撮像方法。
- 前記第1および第2の入力光ビームを前記試料上に向けることは、前記第1または第2の入力光ビームの焦点を前記試料の深さにわたって走査することなく、前記第1および第2の入力光ビームを向けることを含む、請求項22に記載の3次元体積撮像方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SG10202105493U | 2021-05-24 | ||
SG10202105493U | 2021-05-24 | ||
PCT/SG2022/050340 WO2022250610A1 (en) | 2021-05-24 | 2022-05-23 | Stimulated raman scattering tomography system and method |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2024519990A true JP2024519990A (ja) | 2024-05-21 |
Family
ID=84230326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2023572624A Pending JP2024519990A (ja) | 2021-05-24 | 2022-05-23 | 誘導ラマン散乱トモグラフィシステム及び方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20240241057A1 (ja) |
JP (1) | JP2024519990A (ja) |
CN (1) | CN117413171A (ja) |
WO (1) | WO2022250610A1 (ja) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8441633B2 (en) * | 2009-10-29 | 2013-05-14 | California Institute Of Technology | Multiple-photon excitation light sheet illumination microscope |
WO2011099269A1 (ja) * | 2010-02-10 | 2011-08-18 | 国立大学法人大阪大学 | 顕微鏡及び観察方法 |
CN110208241B (zh) * | 2019-05-11 | 2021-10-26 | 复旦大学 | 基于受激拉曼散射的大气单颗粒物快速三维化学成像方法 |
EP4089400B1 (en) * | 2020-01-09 | 2024-03-06 | National Univeristy Corporation Tokyo University Of Agriculture And Technology | Light detection device, and light detection method |
CN212489863U (zh) * | 2020-02-27 | 2021-02-09 | 浙江大学 | 一种快速高效自适应光学补偿的受激拉曼散射成像系统 |
-
2022
- 2022-05-23 CN CN202280038733.7A patent/CN117413171A/zh active Pending
- 2022-05-23 JP JP2023572624A patent/JP2024519990A/ja active Pending
- 2022-05-23 US US18/563,631 patent/US20240241057A1/en active Pending
- 2022-05-23 WO PCT/SG2022/050340 patent/WO2022250610A1/en active Application Filing
Also Published As
Publication number | Publication date |
---|---|
CN117413171A (zh) | 2024-01-16 |
US20240241057A1 (en) | 2024-07-18 |
WO2022250610A1 (en) | 2022-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ogien et al. | Dual-mode line-field confocal optical coherence tomography for ultrahigh-resolution vertical and horizontal section imaging of human skin in vivo | |
Hu et al. | Harmonic optical tomography of nonlinear structures | |
US9360660B2 (en) | Methods and apparatus for laser scanning structured illumination microscopy and tomography | |
CN101802675B (zh) | 荧光焦点调制显微系统和方法 | |
Soliman et al. | Combining microscopy with mesoscopy using optical and optoacoustic label-free modes | |
US8848199B2 (en) | Tomographic phase microscopy | |
Vogler et al. | Multimodal imaging spectroscopy of tissue | |
Gong et al. | Stimulated Raman scattering tomography enables label‐free volumetric deep tissue imaging | |
Yoshitake et al. | Direct comparison between confocal and multiphoton microscopy for rapid histopathological evaluation of unfixed human breast tissue | |
CN105852800B (zh) | 三维光学相干断层扫描装置及其应用 | |
CA2882784C (en) | Three-dimensional optical coherence tomography apparatus and its application | |
Stockton et al. | Single pixel quantitative phase imaging with spatial frequency projections | |
US9791682B2 (en) | Quantitative nonlinear optical microscopy using a shaped beam | |
Lin et al. | Volumetric chemical imaging in vivo by a remote-focusing stimulated Raman scattering microscope | |
CN113390827A (zh) | 一种基于相干探测的多光子超深度组织成像设备及其检测成像方法 | |
Wang et al. | Methods and applications of full-field optical coherence tomography: a review | |
Quan et al. | Multimodal microscopy: fast acquisition of quantitative phase and fluorescence imaging in 3D space | |
Lee et al. | In vivo sub-femtoliter resolution photoacoustic microscopy with higher frame rates | |
Soltani et al. | Label-free detection of brain tumors in a 9L gliosarcoma rat model using stimulated Raman scattering-spectroscopic optical coherence tomography | |
Heuke et al. | Coherent Stokes Raman scattering microscopy (CSRS) | |
Lin et al. | Super‐Resolution Two‐Photon Fluorescence Tomography Through the Phase‐Shifted Optical Beatings of Bessel Beams for High‐Resolution Deeper Tissue 3D Imaging | |
Wang et al. | Stimulated Raman scattering tomography for rapid three-dimensional chemical imaging of cells and tissue | |
Hu et al. | Noncontact photoacoustic lipid imaging by remote sensing on first overtone of the CH bond | |
US20240241057A1 (en) | Stimulated Raman Scattering Tomography System And Method | |
Nteroli et al. | Two octaves spanning photoacoustic microscopy |