• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    熒光壽命數(shù)據(jù)的相量分析及其應(yīng)用*

    2020-08-29 07:37:20林丹櫻牛敬敬劉雄波張瀟張嬌于斌屈軍樂
    物理學(xué)報 2020年16期
    關(guān)鍵詞:端點光子壽命

    林丹櫻 牛敬敬 劉雄波 張瀟 張嬌 于斌 屈軍樂

    (深圳大學(xué)物理與光電工程學(xué)院, 光電子器件與系統(tǒng)教育部/廣東省重點實驗室, 深圳 518060)

    1 引 言

    隨著熒光標記技術(shù)的快速發(fā)展, 熒光顯微成像技術(shù)已成為生物醫(yī)學(xué)研究領(lǐng)域不可或缺的工具之一. 熒光壽命作為熒光的一個重要參量, 反映了熒光分子從激發(fā)態(tài)回到基態(tài)的退激發(fā)速率, 是熒光分子的固有特性, 其變化能非常靈敏地反映熒光分子所處微環(huán)境(如細胞微環(huán)境中的溫度、黏度、pH 值、離子濃度等)的變化情況[1]. 同時由于熒光壽命的測量可不受熒光探針濃度、激發(fā)光強度、光漂白等因素的影響, 相比強度測量更有利于實現(xiàn)定量化, 因此探測樣品中熒光壽命的分布和變化的熒光壽命顯微成像技術(shù)(fluorescence lifetime imaging microscopy, FLIM)常常被應(yīng)用于定量測量細胞內(nèi)的一些生物物理和生物化學(xué)參數(shù)[2,3], 或氣溶膠的黏度測量等[4]. 另一方面, 許多熒光團分子具有相似的熒光光譜, 假如利用它們對樣品中不同的結(jié)構(gòu)進行特異性標記, 直接從光譜上很難分離, 但利用它們熒光壽命不同的特點, 可借助FLIM 技術(shù)對其熒光進行區(qū)分, 從而實現(xiàn)多結(jié)構(gòu)的同時標記和成像. 這些特點使得FLIM 在生命科學(xué)研究中有著越來越廣泛的應(yīng)用.

    一般而言, FLIM 中熒光壽命的測量方法有時域法與頻域法兩類, 其中頻域法使用周期調(diào)制的連續(xù)光激發(fā)樣品, 檢測熒光信號相對激發(fā)光的振幅和相位變化來計算樣品的熒光壽命; 時域法采用高重復(fù)頻率的飛秒脈沖激光激發(fā)樣品, 利用門控技術(shù)、掃描相機技術(shù)、時間相關(guān)單光子計數(shù)(timecorrelated single photon counting, TCSPC)技術(shù)等記錄脈沖過后熒光信號的衰減過程并擬合或計算熒光壽命. 因此, 不管對于哪種類型的FLIM, 熒光壽命數(shù)據(jù)的處理都是非常關(guān)鍵的一環(huán). 例如, 時域法熒光壽命數(shù)據(jù)最常用的處理方法是非線性最小二乘(nonlinear least-squares, NL-LS)擬合法,即選擇一種衰減模型(單指數(shù)、雙指數(shù)或多指數(shù)衰減)對每個像素點的測量數(shù)據(jù)進行曲線擬合, 然后分析出各衰減組分的熒光壽命值. 但是這種方法要求衰減模型的選擇要合適, 同時要求每個像素點要采集足夠多的光子數(shù)(通常 >1000)才能夠進行有效擬合[5]. 這就使得基于這種數(shù)據(jù)處理方法的FLIM 技術(shù)在應(yīng)用上受到了相應(yīng)的限制. 例如,TCSPC-FLIM 通常需要在激發(fā)光較弱的情況下進行成像, 以避免相鄰兩個脈沖之間檢測到的光子數(shù)多于1 而導(dǎo)致的光子堆積問題, 因此要求每個像素點采集足夠多的光子數(shù)就需要延長采集時間, 由此帶來的問題就是獲得一幅熒光壽命圖像的時間相對較長, 成像速度受限, 無法應(yīng)用于一些熒光壽命變化速度較快的場合.

    因此, 針對這些問題, 研究者們開發(fā)出了其他的算法, 例如相量分析(phasor analysis, PA)[6]、極大 似 然 估 計 (maximum likelihood estimate,MLE)[7]、一階矩(the first moment, M1)[8]、貝葉斯 分 析(Bayesian analysis, BA)[9]、壓 縮 感 知(compressed sensing, CS)[10]等, 這些方法可通過降低壽命分析對光子數(shù)的要求從而間接地提高FLIM 技術(shù)的成像速度[11]. 其中, PA 法通過將時域信息變換到頻域, 可直接計算得到熒光壽命值,簡單快速, 無需擬合, 且在低光子數(shù)情況下也能得到比較準確的壽命值, 適合于快速FLIM 成像. 另一方面, PA 法生成的相量圖能將具有相似熒光衰減特性的熒光團分子所對應(yīng)的像素點顯示在相鄰區(qū)域, 形成一定的簇狀分布, 這種特點使得利用該方法可以更方便地對數(shù)據(jù)進行可視化和聚類分析,因此結(jié)合PA 法的FLIM 技術(shù)(phasor-FLIM)越來越受到科研人員的青睞, 在生命科學(xué)和生物醫(yī)學(xué)研究中應(yīng)用越來越廣泛.

    本文首先詳細闡述phasor-FLIM 的基本原理及使用方法, 并在此基礎(chǔ)上介紹該技術(shù)應(yīng)用于細胞代謝狀態(tài)測量、蛋白質(zhì)相互作用研究、細胞微環(huán)境測量以及輔助病理診斷和提高超分辨成像分辨率等方面的最新研究進展, 最后對其發(fā)展前景進行展望.

    2 Phasor-FLIM 的基本原理

    如前所述, FLIM 技術(shù)中熒光壽命的測量通常分為頻域法和時域法[1,6]. 這兩種方法在數(shù)學(xué)意義上互為傅里葉變換, 但它們獲取熒光壽命信息的方式不同, 得到的數(shù)據(jù)內(nèi)容和形式不同, 從而數(shù)據(jù)處理方法一般也不同. 頻域法一般使用正弦調(diào)制的連續(xù)光激發(fā)樣品, 測量得到的是具有相同頻率的熒光信號, 但由于熒光壽命的影響, 熒光信號的振幅和相位相比激發(fā)光均發(fā)生了變化, 因此通過計算熒光信號相對激發(fā)光的振幅調(diào)制度變化和相位延遲可計算得到熒光壽命. 而時域法則需要采用高重復(fù)頻率的飛秒脈沖激光激發(fā)樣品, 利用前面提到的門控技術(shù)、掃描相機或TCSPC 技術(shù)等直接或間接記錄脈沖過后的熒光衰減過程, 得到的是熒光強度(或光子數(shù))隨時間的變化關(guān)系, 因此一般可通過曲線擬合得到熒光壽命.

    PA 法最先被用于處理頻域FLIM 技術(shù)得到的熒光壽命數(shù)據(jù), 其相量由頻域FLIM 測量得到的解調(diào)系數(shù)和相位延遲來構(gòu)建, 是原始數(shù)據(jù)的直接表達[12]. PA 法同樣適用于時域FLIM 數(shù)據(jù)的分析,但需要先將時域的熒光衰減變換到頻域. 由于時域FLIM 中的TCSPC-FLIM 目前應(yīng)用最為廣泛,因此PA 法在該技術(shù)中的應(yīng)用也是報道得最多的.以下分別介紹這兩類技術(shù)中PA 法分析熒光壽命的基本原理, 并結(jié)合熒光相量圖的特點闡述其典型的應(yīng)用思路.

    2.1 頻域法FLIM 及其相量分析

    在頻域FLIM 中, 激發(fā)光常采用正弦調(diào)制的連續(xù)光源, 如激光、氙燈、LED 燈等, 激發(fā)光強可描述為

    其中E(0)和E(t)分別為起始時刻和t時刻的激發(fā)光強度,ME=a/A為激發(fā)光的強度調(diào)制度(a和A分別是激發(fā)光的振幅和平均強度, 如圖1(a)所示),w為調(diào)制的角頻率. 利用該激發(fā)光激發(fā)標記有熒光團的樣品后, 所檢測到的熒光也是正弦調(diào)制的, 且其頻率與激發(fā)光相同, 但強度調(diào)制度會降低,相位也有一定延遲. 如果熒光衰減符合單指數(shù)衰減規(guī)律, 則熒光光強可描述為

    圖1 熒光壽命的測量方法及相量分析(PA)法示意圖 (a)頻域法測量原理示意圖; (b)單指數(shù)衰減的壽命相量示例圖; (c)雙指數(shù)衰減的壽命相量示例圖; (d)時間相關(guān)單光子計數(shù)(TCSPC)測量原理示意圖Fig. 1. Schematic diagram of fluorescence lifetime measurement and phasor analysis (PA): (a) Frequency domain method; (b) lifetime phasor of single-exponential decay; (c) lifetime phasor of bi-exponential decay; (d) time-correlated single photon counting (TCSPC) method.

    相應(yīng)地,I(0)和I(t)分別為起始時刻和t時刻的熒光強度,MF=b/B為熒光的強度調(diào)制度(b和B分別是熒光的振幅和平均強度, 如圖1(a)所示),f為相位延遲或相移. 通過測量相移f和解調(diào)系數(shù)M=MF/ME, 即可計算出相應(yīng)的熒光壽命, 即:

    可以證明[12,13], 當(dāng)熒光的衰減過程符合單指數(shù)衰減規(guī)律(即I(t) =I(0)e–t/t,t為壽命)時, (3)式和(4)式算得的壽命tf和tm是相等的, 理論上只需要計算其中一個即可. 但假如熒光的衰減過程是兩個或多個單指數(shù)衰減過程的組合(即雙指數(shù)衰減或多指數(shù)衰減), 則情況更加復(fù)雜, 通常需要先在不同調(diào)制頻率wi下重復(fù)測量多組fi和Mi, 并通過以下公式先計算出等效相移和解調(diào)系數(shù)再求解壽命, 即:

    其中fi為第i個調(diào)制頻率測得的光強占總光強的比例.

    1984 年, Jameson 等[13]利用相量的概念對頻域法FLIM 得到的數(shù)據(jù)進行幾何表示, 他們利用單個像素點對應(yīng)的解調(diào)系數(shù)M和相移f來構(gòu)建一個相量, 即以M作為該相量的模, 以f作為該相量的輻角, 則可以認為相量與像素點是一一對應(yīng)的, 相量圖上一個相量的端點就代表了一個像素點的全部熒光壽命信息(如圖1(b)所示). 該相量在實軸和虛軸的分量可用Weber 符號表示[14], 即:

    對于單指數(shù)衰減情形, 由(3)式和(4)式可得到cosf=M, 因此可以得到:

    即以坐標(G,S)表示的相量端點被約束在圓心位于(0.5, 0)處、半徑為0.5 的半圓上. 半圓上的每個點表示不同的壽命, 壽命值從左到右遞減, 其中(1, 0)表示接近零的壽命, (0, 0)表示無限長的壽命, 如圖1(b)所示.

    幾何意義上, 多指數(shù)衰減過程對應(yīng)的相量端點應(yīng)位于其各個單指數(shù)衰減組分對應(yīng)的壽命相量端點連接組成的集合內(nèi). 如圖1(c)中, 雙指數(shù)衰減過程對應(yīng)的壽命相量端點(藍色)落在半圓以內(nèi), 位于兩個單指數(shù)衰減組分對應(yīng)的壽命相量端點(綠色)的連線上, 且與兩端點的距離(p1,p2)由兩個組分的占比(a1,a2)決定.

    2.2 時域法FLIM 及其相量分析

    PA 法同樣適用于時域法FLIM 數(shù)據(jù)的分析.這里以目前應(yīng)用最廣泛的TCSPC-FLIM 技術(shù)為例. 如圖1(d)所示, TCSPC 將每一次脈沖信號作為一個信號周期, 每個周期內(nèi)當(dāng)探測到第一個熒光光子時就在其到達時間對應(yīng)的時間通道中進行計數(shù), 經(jīng)過多次累積即可建立一個反映熒光衰減過程的光子數(shù)-時間分布直方圖, 用于求解熒光壽命.

    2008 年, Digman 等[16]將Weber[14]于1981年提出的熒光脈沖響應(yīng)的傅里葉分析方法用于TCSPC 技術(shù), 并使用以下關(guān)系式計算相量端點坐標, 也可以將每個像素的全部熒光衰減信息轉(zhuǎn)換為相量圖上的單個點, 即:

    其中G(w)和S(w)分別是熒光脈沖響應(yīng)傅里葉頻譜的實部和虛部, 這里用作壽命相量的兩個分量;I(t)是脈沖過后t時刻的熒光光強,w是脈沖激光的角頻率,T是周期長度(實際可取動態(tài)范圍),n是諧波的階次, 一般情況下取基頻分量, 即n=1. 因為TCSPC 是在N個離散的時間通道tk中記錄熒光衰減分布直方圖的, 因此上述積分運算可轉(zhuǎn)化為求和計算:

    其中Ck是第k個時間通道記錄的光子數(shù),是所有時間通道記錄的總光子數(shù).

    如果熒光脈沖響應(yīng)呈單指數(shù)衰減規(guī)律, 則代入(11)式可得G(w)和S(w)與熒光壽命t的關(guān)系為:

    可以驗證,G,S仍然滿足(8)式的關(guān)系, 即單指數(shù)衰減情形下以坐標(G,S)表示的相量端點仍被約束在圓心位于(0.5, 0)處、半徑為0.5 的半圓上, 與頻域法得到的相量圖一致.

    而如果熒光脈沖響應(yīng)是呈多指數(shù)衰減的, 則同樣根據(jù)線性疊加關(guān)系可知, 坐標(G,S)表示的相量端點也位于半圓以內(nèi), 其中雙指數(shù)衰減過程對應(yīng)的相量端點也在其兩個單指數(shù)衰減組分對應(yīng)的相量端點連線上, 與頻域法得到的結(jié)論一致.

    2.3 相量圖的特點及其典型應(yīng)用思路

    如前所述, 對于滿足單指數(shù)衰減規(guī)律的情形,熒光壽命相量圖上以坐標(G,S)表示的相量端點被約束在圓心位于(0.5, 0)處、半徑為0.5 的半圓上, 而多指數(shù)衰減過程的相量端點則位于半圓以內(nèi), 具體地說, 是位于各單指數(shù)衰減組分相量端點連接組成的集合內(nèi). 例如, 對于雙指數(shù)衰減情形,其壽命相量端點位于兩個單指數(shù)衰減組分相量端點的連線上, 且位置與兩個組分的占比有關(guān). 相量圖的這些特點, 使得其在熒光壽命數(shù)據(jù)的定量分析、可視化和聚類分析方面有很大優(yōu)勢.

    如前所述, 熒光壽命反映的是熒光分子從激發(fā)態(tài)回到基態(tài)的退激發(fā)速率, 因此當(dāng)處于激發(fā)態(tài)的熒光分子所處的微環(huán)境不同, 或者熒光分子與其他分子發(fā)生相互作用和能量轉(zhuǎn)移時, 熒光壽命會發(fā)生靈敏的變化. 所以, 許多熒光團存在兩個甚至兩個以上的衰減速率, 分別對應(yīng)于熒光團的不同狀態(tài), 例如細胞內(nèi)的還原型煙酰胺腺嘌呤二核苷酸(nicotinamide adenine dinucleotide, NADH)處于自由態(tài)時(主要位于細胞質(zhì)中)的熒光壽命為幾百皮秒(t1), 與蛋白質(zhì)綁定后處于結(jié)合態(tài)時(主要位于線粒體中)的熒光壽命則達到幾納秒(t2)[6].2005 年Brid 等[17]利用干擾NADH/NAD+比例的代謝藥物氰化鉀(Potassium cyanide, KCN)阻斷細胞的呼吸作用, 證明了細胞中自由態(tài)和結(jié)合態(tài)NADH 的占比與NADH / NAD+比例有關(guān), 從而提出利用FLIM 測量自由態(tài)和結(jié)合態(tài)NADH 的比例, 來間接獲得活細胞內(nèi)的NADH/NAD+比例.NADH 與其氧化形式NAD+作為生物體內(nèi)許多氧化還原反應(yīng)的輔酶參與生命活動并相互轉(zhuǎn)化, 它們之間的平衡反映了氧化磷酸化和糖酵解的比率, 因此NADH/NAD+比例及其變化可用于監(jiān)測細胞代謝方式的變化. 對于細胞內(nèi)NADH/NAD+比例的探測, 傳統(tǒng)的方法有酶循環(huán)法、毛細管電泳法、質(zhì)譜分析法等, 但這些方法只適用于對細胞提取物進行檢測. 直接對細胞內(nèi)NADH/NAD+比例進行探測也可以采用熒光強度成像的方法, 這是因為NADH 發(fā)熒光而NAD+不發(fā)熒光. 但由于熒光團濃度的不均一性以及自由態(tài)和結(jié)合態(tài)NADH 熒光量子產(chǎn)率的不同, 基于強度成像的NADH/NAD+測量結(jié)果通常存在較大誤差. Brid 等[17]的工作為定量監(jiān)測活細胞內(nèi)NADH/NAD+的比例提供了一種新的思路, 即通過FLIM 定量測量自由態(tài)和結(jié)合態(tài)NADH 的比例來間接獲得NADH/NAD+比例. 他們的研究還表明, 當(dāng)NADH/NAD+比例發(fā)生變化時, 自由態(tài)和結(jié)合態(tài)NADH 的熒光壽命值(即t1和t2)幾乎保持恒定, 變化的是它們各自的占比(即a1和a2). 多指數(shù)衰減情形的壽命相量與其各個單指數(shù)衰減組分壽命相量之間的線性關(guān)系,使得采用了PA 法分析熒光壽命的phasor-FLIM 特別適合于在已知單指數(shù)衰減組分壽命的前提下, 對多組分尤其是雙組分中不同組分的占比進行定量分析. 所以, phasor-FLIM 的典型應(yīng)用之一, 就是根據(jù)壽命相量端點的位置對NADH 兩種壽命組分的占比進行定量分析, 從而用于細胞代謝狀態(tài)的測量[18?27].

    此外, 在利用熒光共振能量轉(zhuǎn)移(fluorescence resonance energy transfer, FRET)對蛋白質(zhì)相互作用進行研究的FLIM-FRET 實驗, 以及一些利用FLIM 研究細胞微環(huán)境參量變化的場合, 也會涉及到雙指數(shù)衰減甚至多指數(shù)衰減及其組分占比的變化, 因此phasor-FLIM 在這些應(yīng)用中的報道近幾年也多了起來[28?34].

    其次, PA 法生成的相量圖中, 具有相似熒光衰減特征的像素點所對應(yīng)的相量端點會顯示在相近的位置, 因而形成一定的簇狀分布(如圖2(a)和圖2(b)所示), 稱為“相量簇”. 同一個相量簇對應(yīng)于具有相似熒光衰減特性的多個像素點, 但這些像素點在樣品中的分布可以是不連續(xù)的, 和它們的空間分布并沒有直接聯(lián)系. 利用該特點, 可在相量圖上選取感興趣的相量簇, 并通過其對應(yīng)關(guān)系找到相量簇中各相量端點所對應(yīng)的像素點在樣品上的位置, 就可以方便地將樣品中具有相似熒光壽命特征的區(qū)域標記出來(如圖2(c)所示), 還可以進一步對這些熒光壽命特征相似的區(qū)域進行聚類分析, 從而為研究一些生物醫(yī)學(xué)問題帶來極大的便利[18?20].倘若對不同的相量簇對應(yīng)的像素點賦以不同的偽彩色, 則利用該方法還可以方便地實現(xiàn)熒光壽命圖像的多色顯示(如圖2(d)所示), 從而起到輔助病理判斷等作用[35?37]. 最近, 這種利用PA 法對熒光壽命的分布進行聚類分析的方法還被應(yīng)用于輔助提高基于受激輻射耗盡(stimulated emission depletion, STED)的超分辨成像的分辨率[38?40].

    圖2 Phasor-FLIM 的應(yīng)用思路示意圖 (a)包含未處理壽命信息的熒光強度圖; (b)經(jīng)PA 法分析得到的壽命相量圖; (c)對壽命相量直接進行分析; (d)通過相量聚類分析和偽彩色標記得到的熒光壽命圖Fig. 2. Schematic diagram of phasor-FLIM application:(a) Fluorescence intensity image with untreated lifetime information; (b) lifetime phasor plot obtained by PA analysis;(c) direct analysis of lifetime phasors; (d) phasor-mapped FLIM image based on phasor clustering analysis and pseudo-color assignment.

    綜上, 目前phasor-FLIM 的應(yīng)用可以用圖3來表示. 以下將舉例說明phasor-FLIM 的應(yīng)用研究進展.

    圖3 Phasor-FLIM 的應(yīng)用分類示意圖Fig. 3. Application classification diagram of phasor-FLIM.

    3 Phasor-FLIM 的應(yīng)用研究進展

    3.1 基于NADH 的細胞代謝應(yīng)用研究

    細胞的能量主要來自細胞呼吸, 正常細胞或正常分化的細胞在有氧條件下采用糖酵解進行代謝,缺氧條件下通過氧化磷酸化進行代謝. 而高度增殖的細胞(如癌細胞或干細胞)即使在氧氣充足的條件下也多選擇糖酵解作為主要的產(chǎn)能方式, 這種現(xiàn)象被稱為“Warburg 效應(yīng)”[41]. 而這種細胞代謝方式的不同導(dǎo)致兩者所含NADH 的濃度和狀態(tài)存在差異. 因此基于前面提到的NADH 處于兩種不同狀態(tài)具有不同熒光壽命的特性, 以NADH 作為內(nèi)源熒光標志物的雙光子FLIM 成像常被用于研究正常細胞、干細胞、癌細胞以及其他疾病發(fā)生時細胞的代謝差異, 而結(jié)合了PA 法的phasor-FLIM可方便地實現(xiàn)活組織中單細胞代謝表型的觀測, 在細胞分化和增殖、疾病的機理研究和診斷等方面均具有很好的應(yīng)用前景, 目前也取得了一些重要的應(yīng)用研究進展.

    細胞分化和增殖過程會改變糖酵解和氧化磷酸化之間的平衡, 所以代謝變化可用于研究細胞分化和增殖狀態(tài). Stringari 等[18?20]通過對細胞內(nèi)包括NADH 在內(nèi)的各種內(nèi)源性熒光標志物進行FLIM 成像, 并利用PA 法對FLIM 圖像進行分割,可以區(qū)分膠原蛋白、視黃醇、視黃酸, 以及處于自由態(tài)和結(jié)合態(tài)的NADH. 他們對小腸進行雙光子FLIM 成像并對利用PA 法得到的壽命相量簇進行聚類分析, 用來識別高度增殖的小腸干細胞, 通過代謝狀態(tài)對小腸干細胞和分化的后代進行分類分析, 以監(jiān)測與代謝變化相關(guān)的生理(病理)過程.Lee 等[21]對白細胞和白血病細胞進行FLIM 成像,因為白血病細胞快速增殖, 糖酵解占主導(dǎo)地位, 壽命相量簇向短壽命方向移動, 利用PA 法定量分析自由態(tài)和結(jié)合態(tài)NADH 比例的變化, 可從血液中快速篩選和分離白血病細胞. 與傳統(tǒng)的生物分子診斷技術(shù)相比, 這種基于phasor-FLIM 的單細胞篩查方法對細胞友好, 具有臨床篩選血液細胞的潛力. 細胞外基質(zhì)(extracellular matrix, ECM)是所有組織必不可少的動態(tài)組成部分, 并通過提供機械和生化信號來直接影響細胞行為. ECM 的變化可以改變組織的動態(tài)平衡, 從而潛在地促進細胞轉(zhuǎn)化和腫瘤的發(fā)生. Romero-Lopez 等[22]將正常細胞接種在提取自正常人結(jié)腸和轉(zhuǎn)移至肝臟的結(jié)腸腫瘤的ECM 中, 利用phasor-FLIM 技術(shù)定量分析了自由態(tài)和結(jié)合態(tài)NADH 的比例, 結(jié)果發(fā)現(xiàn)接種在腫瘤ECM 中的細胞比接種在正常ECM 中的細胞具有更高的游離NADH 水平, 糖酵解速率較高, 表明ECM 在癌細胞及其相關(guān)脈管系統(tǒng)的生長中起到了重要作用.

    圖4 Phasor-FLIM 用于分析細胞在缺氧和線粒體毒性藥物氰化鉀刺激下NADH/NAD+比例的變化, 研究代謝狀態(tài)的轉(zhuǎn)變[27]Fig. 4. Phasor-FLIM was used to analyze the change of NADH/NAD+ ratio under the stimulation of hypoxia and mitochondrial toxic drug potassium cyanide, for studying the change of metabolic state of cells[27].

    由于許多疾病的發(fā)病機制與細胞代謝也有著密切的關(guān)聯(lián), phasor-FLIM 在疾病的機理研究及診斷方面也有不少應(yīng)用實例. 例如, 亨廷頓病(Huntington’s disease, HD)是一種常染色體神經(jīng)退行性疾病, 能量代謝障礙是HD 的主要發(fā)病機制. Sameni 等[23]使用phasor-FLIM 來定量測量HD 中自由態(tài)和結(jié)合態(tài)NADH 的比例變化, 作為活細胞代謝變化的間接測量, 用以研究HD 發(fā)病機制, 結(jié)果表明HD 的代謝障礙為糖酵解增加, 導(dǎo)致氧化應(yīng)激和細胞死亡. 這種定量分析的方法可用于篩選HD 組織并進行潛在的藥物篩選, 對診斷和治療疾病具有重要意義. 阿茲海默癥(Alzheimer’s disease, AD)是一種老年人多發(fā)的神經(jīng)退行性疾病, 與抗氧化保護降低和線粒體功能障礙有關(guān).Dong 等[24,25]使用雙光子激發(fā)FLIM 技術(shù)結(jié)合PA 法定量分析了老年小鼠海馬區(qū)神經(jīng)元游離NADH 的水平, 結(jié)果發(fā)現(xiàn)隨著年齡的增加, 游離的NADH 濃度降低, 氧化磷酸化占據(jù)主導(dǎo)地位,AD 進一步惡化, 而還原性治療可恢復(fù)老年小鼠及AD 小鼠神經(jīng)元中游離NADH 的水平. Hato 等[26]則采用類似的方法對活體小鼠腎臟進行雙光子FLIM 成像, 使用PA 法分析腎臟中的代謝變化,提供了一種研究腎臟疾病代謝的方法. Datta 等[27]通過雙光子FLIM 對人誘導(dǎo)多能干細胞分化的心肌 細 胞(human induced pluripotent stem cellderived cardiomyocytes, hiPS-CMs)進行成像, 檢測缺氧和線粒體毒性藥物氰化鉀病理刺激下代謝狀態(tài)的變化. 如圖4 所示, 缺氧狀態(tài)下和用藥物刺激時hiPS-CMs 的壽命相量分布向自由態(tài)NADH的方向移動, 代謝方式轉(zhuǎn)變?yōu)樘墙徒? 這種非侵入性成像技術(shù)有助于研究心臟病的發(fā)病機理和治療方法[27].

    3.2 基于FLIM-FRET 的蛋白互作研究

    當(dāng)一個熒光分子(供體)的熒光發(fā)射光譜與另一個熒光分子(受體)的激發(fā)光譜相重疊, 且兩者間的距離合適(一般為6—10 nm)時, 供體的激發(fā)能誘發(fā)受體發(fā)出熒光, 同時其自身的熒光強度發(fā)生衰減, 稱為FRET 效應(yīng). 由于該效應(yīng)發(fā)生的程度(即FRET 效率)與供體和受體的距離緊密相關(guān),該效應(yīng)常被用于研究蛋白質(zhì)間的相互作用(簡稱“蛋白互作”). 具體的做法是將供體熒光分子和受體熒光分子分別標記于兩個蛋白質(zhì)分子上, 通過測量FRET 效率來反映兩個蛋白質(zhì)之間的距離, 從而反映其是否發(fā)生相互作用以及相互作用的程度.從熒光壽命的角度看, 當(dāng)FRET 效應(yīng)發(fā)生時, 由于供體將能量轉(zhuǎn)移到受體上, 供體的熒光壽命將變短. 考慮到FLIM 測量不受光漂白等因素的影響,利用FLIM 進行FRET 效率的測量比測量熒光強度變化的方法要更加準確, 因此FLIM-FRET 已被廣泛應(yīng)用于蛋白互作的研究. 發(fā)生和不發(fā)生FRET(或者說FRET 效率很高和很低)兩種狀態(tài)下供體的壽命可認為是不變的, 相當(dāng)于兩個單指數(shù)衰減組分, 因此根據(jù)PA 法的分析, 兩種狀態(tài)下得到的壽命相量應(yīng)位于半圓上. 而當(dāng)部分供體發(fā)生FRET 時, 總體的熒光衰減滿足雙指數(shù)衰減規(guī)律,壽命相量端點將位于上述兩個單組分壽命相量端點的連線上, 且到兩端的距離與兩種組分的占比有關(guān), 由此可方便地計算出FERT 效率. 2012 年和2013 年, Hinde 等[28,29]利 用 phasor-FLIM 和FRET 監(jiān)測了小G 蛋白與RBD 蛋白在溶血磷脂酸(lysophosphatidic acid, LPA)刺激下相互作用的變化, 結(jié)果表明, 在熒光強度圖像看不出明顯差異的情況下, FLIM 圖像和PA 定量分析均顯示小G 蛋白與RBD 蛋白的標簽蛋白EPCF 和Critine在LPA 刺激下FRET 效率顯著增加(如圖5 所示). Lou 等[30]在 共 表 達H2B-eGFP 和H2BmCherryrry兩種熒光蛋白的Hela 細胞中利用phasor-FLIM 測得的FRET 效率對核小體的緊實度進行了量化, 用于定量反應(yīng)DNA 的損傷情況.

    圖5 Phasor-FLIM 用于定量測量RhoA-kRas 單鏈生物傳感器的熒光共振能量轉(zhuǎn)移(FRET)效率, 研究蛋白互作[28]Fig. 5. Phasor-FLIM was used in quantitative FRET efficiency detection of a RhoA-kRas single chain biosensor, studying interaction between proteins[28].

    3.3 基于敏感探針的細胞微環(huán)境參量測量

    由于熒光壽命與熒光分子所處的微環(huán)境密切相關(guān), 因此采用一些對細胞中某種微環(huán)境參量(如溫度、pH 值、離子濃度、黏度等)特別敏感的熒光探針對這些微環(huán)境參量及其變化進行定量測量是FLIM 的常規(guī)應(yīng)用. Chen 等[31]將PA 法應(yīng)用于溫度的FLIM 測量, 通過熒光分子羅丹明B 的壽命相量圖定量分析比較了Hela 細胞在8 種不同溫度下壽命相量簇的分布變化. Battisti 等[32]采用具有pH 依賴性的熒光團E2GFP, 利用phasor-FLIM 將所有的實驗圖像在同一個相量圖上進行比較, 從而分析細胞在生理狀態(tài)和非生理狀態(tài)下pH 值的變化(如圖6 所示). 深圳大學(xué)Zhou 等[33]使用phasor-FLIM 分析監(jiān)測了pH 敏感的細胞間納米藥物的釋放. 他們合成聚合物納米顆粒PAHCit/DOX, 用于阿霉素(DOX)在癌細胞中的有效釋放, 并對其釋放過程進行了監(jiān)測. 具體來說, 這種PAH-Cit/DOX 納米顆粒在生理pH 下穩(wěn)定, 但在弱酸性條件下能夠有效釋放DOX. 同時, 他們采用phasor-FLIM 對該納米粒子在癌細胞中釋放DOX 的熒光壽命進行了監(jiān)測, 從而提供了一種評估納米載體藥物釋放效率的方法. Ferri 等[34]則開發(fā)了一種基于BODIPY 的黏度敏感熒光分子BoMe, 將BoMe 與phasor-FLIM 相結(jié)合用于評估早衰綜合癥(Hutchinson-Gilford progeria syndrome,HGPS)細胞內(nèi)的黏度.

    3.4 輔助病理診斷

    圖6 Phasor-FLIM 用于分析細胞在正常(靜止)狀態(tài)和氧化應(yīng)激狀態(tài)下pH 值的變化[32]Fig. 6. Phasor-FLIM was used to analyze the changes of pH value of cells in normal state (at rest) and under oxidative stress[32].

    非黑色素瘤皮膚癌(non-melanoma skin cancers, NMSC)是白種人群體中最常見的腫瘤疾病, 光化性角病(actinic keratosis, AK)、博文病(Bowen disease, BD)、基 底 細 胞 癌(basal cell carcinoma, BCC)和浸潤性鱗狀細胞癌(invasive squamous cell carcinoma, SCC)都是其常見的亞型. 該類型皮膚疾病的臨床癥狀很不明顯, 因而采用活組織實現(xiàn)對NMSC 的診斷比較困難, 目前臨床上主要是通過對蘇木精和伊紅(hematoxylin and eosin, H&E)染色切片的形態(tài)學(xué)特征分析來進行診斷. 但這種方法鑒別與診斷上述不同NMSC亞型的準確性通常與病理學(xué)家的專業(yè)知識、技能和經(jīng)驗有關(guān), 可能會導(dǎo)致不同人員鑒別與診斷的結(jié)果存在巨大差異, 主觀性大, 出錯概率較高, 影響后續(xù)治療. 2017 年, 深圳大學(xué)的Luo 等[35]對H&E 染色的皮膚病理切片進行了雙光子FLIM 成像, 并結(jié)合PA 法定量分析了AK, BD, BCC 這三種類型的皮膚腫瘤對應(yīng)的壽命相量簇分布情況, 結(jié)果表明不同皮膚腫瘤疾病的壽命相量簇其坐標值、角斜率、相量面積分布均存在差異, 提示phasor-FLIM 技術(shù)可為鑒別診斷這些皮膚腫瘤疾病提供一種簡便可行的組織病理學(xué)分析方法, 其準確性和客觀性相對于明場H&E 診斷將有顯著提高. 此外,他們還采用PA 法對H&E 染色的BCC 細胞切片的熒光壽命數(shù)據(jù)進行了半定量分析[36], 并通過聚類分析對不同的相量簇賦以不同的偽彩色, 從而使FLIM 圖像上的病理學(xué)特征更加明顯, 可進一步輔助病理診斷和促進細胞病變及發(fā)展機理方面的深入研究(如圖7 所示).

    輸尿管的阻塞會增加腎小管的壓力, 降低腎小球的濾過率, 并通過腎臟激活許多血管活性激素和細胞因子, 導(dǎo)致纖維化. 臨床上常用Picrosirius 紅或Masson Trichrome 染色的組織切片, 依靠病理學(xué)家專業(yè)的知識來診斷間質(zhì)纖維化. Ranjit 等[37]利用phasor-FLIM 研究了小鼠腎臟單側(cè)輸尿管梗阻(unilateral ureteral obstruction, UUO)纖維化的發(fā)展過程, 根據(jù)相量分布的平均位置、形狀、角度和數(shù)目, 區(qū)分了健康組織和病變的纖維化組織,同時可以從FLIM 圖像定量得到纖維化程度, 可用于監(jiān)測疾病發(fā)展或治療效果.

    3.5 提升超分辨成像分辨率

    圖7 Phasor-FLIM 聚類分析和偽彩色標記用于增強H&E 染色基底細胞癌(BCC)切片病理學(xué)特征的可視化, 可輔助病理診斷[36]Fig. 7. Phasor-FLIM clustering analysis and pseudo-color assignment was used to enhance visualization of pathological features of basal cell carcinoma (BCC) sections stained with H&E, assisting pathological diagnosis[36].

    圖8 Phasor-FLIM 聚類分析用于濾除受激輻射耗盡(STED)成像中環(huán)形擦除光區(qū)域的光子, 可輔助提升超分辨成像分辨率[39]Fig. 8. Phasor-FLIM cluster analysis was used to filter out the photons in the annular depletion region in stimulated radiation depletion (STED) imaging, improving the resolution of super-resolution imaging[39].

    STED 技術(shù)是近年來發(fā)展起來的一種超分辨顯微成像技術(shù), 在“突破”分辨率極限和相關(guān)應(yīng)用方面都取得了顯著成果[42]. 該技術(shù)使用高強度環(huán)形擦除光使衍射極限范圍內(nèi)除中心點以外的熒光分子發(fā)生受激輻射而不產(chǎn)生熒光(稱為“擦除”), 從而有效減小系統(tǒng)的點擴散函數(shù)獲得超分辨率. 提高擦除光功率可進一步提高分辨率, 但代價是增加了光漂白、光毒性和反斯托克斯發(fā)射背景, 從而降低了信噪比, 并嚴重限制了STED 技術(shù)在活細胞成像中的應(yīng)用. STED 中擦除光的作用, 是使處于激發(fā)態(tài)的粒子發(fā)生受激輻射而迅速躍遷回基態(tài), 相當(dāng)于縮短了這部分熒光分子的壽命. 因此, 受到擦除光照射的熒光分子與未受到擦除光照射的熒光分子會分布在相量圖中的不同位置, 對壽命相量進行聚類分析, 可將兩種狀態(tài)下的熒光分子分離. 2015年, Lanzanò等[38]利用該原理, 將phasor-FLIM與STED 結(jié)合, 在較低擦除光功率的條件下通過濾除掉環(huán)形區(qū)域的光子實現(xiàn)了與較高擦除功率相當(dāng)?shù)姆直媛? 2018 年, 深圳大學(xué)的Wang 等[39]也將phasor-FLIM 應(yīng)用于STED, 通過只選擇激發(fā)光斑中心未被擦除光照射的目標光子進行圖像重構(gòu), 從而在不增加擦除光能量的前提下成功將 光 子 進 行 分 離, 將 分 辨 率 從150 nm 提 升到80 nm (如圖8 所示). 類似地, Tortarolo 等[40]也利用PA 法實現(xiàn)了光子分離, 提出了pSTEDSPLIT 方法, 在不增加擦除光功率的前提下可有效提高分辨率.

    4 總結(jié)與展望

    FLIM 技術(shù)通過采集熒光探針分子的壽命信息進行成像, 具有特異性強、靈敏度高、可定量測量的優(yōu)勢, 常用于探測細胞微環(huán)境參量分布、能量代謝狀態(tài), 表征材料內(nèi)部結(jié)構(gòu)及性質(zhì)等方面, 因而在生物醫(yī)學(xué)、材料學(xué)等諸多領(lǐng)域有廣泛的應(yīng)用. 為了從采集到的大通量壽命數(shù)據(jù)中快速精準地獲取樣品內(nèi)部的結(jié)構(gòu)與功能信息, 發(fā)展一種簡單快速而且實用的壽命數(shù)據(jù)分析方法在實際生物醫(yī)學(xué)應(yīng)用中具有重要的意義. 近些年興起的PA 法在頻域計算壽命, 無需擬合, 分析簡便, 對科研人員的專業(yè)知識背景要求不高, 因此發(fā)展和應(yīng)用均非常迅速.相比傳統(tǒng)的NL-LS 擬合, 該方法不僅分析速度快,相同壽命精度下要求的采集光子數(shù)少, 間接縮短了數(shù)據(jù)采集時間, 有利于發(fā)展快速FLIM 成像技術(shù),而且能將具有相似衰減特征的像素點以“相量簇”的形式呈現(xiàn), 具有強大的數(shù)據(jù)可視化功能, 也便于對數(shù)據(jù)進行聚類分析等, 因而在基于NADH 的細胞代謝狀態(tài)測量、基于FLIM-FRET 的蛋白互作研究、基于敏感探針的細胞微環(huán)境參量測量和輔助病理診斷、提升超分辨成像分辨率等方面均取得了很好的應(yīng)用研究進展. 基于此, 本文詳細闡述了PA 法的基本原理及其運用方法, 并在此基礎(chǔ)上介紹了運用PA 法的phasor-FLIM 在上述幾個方面的應(yīng)用實例, 著重討論了PA 法在這些FLIM 應(yīng)用實例中的優(yōu)勢所在, 為相關(guān)領(lǐng)域的研究提供有益的參考.

    然而, 盡管近年來phasor-FLIM 在實際應(yīng)用中已經(jīng)取得了許多重要的研究進展, 但與此同時還存在一些局限性. 例如, 雖然PA 處理時對采集光子數(shù)的要求要比傳統(tǒng)擬合算法低很多, 但是在快速FLIM 成像時當(dāng)采集的光子數(shù)較少或信噪比較差時, 其相量簇的分布就會比較雜散, 據(jù)此分析得到的結(jié)果可靠性往往不夠高. 為了解決這一問題,目前已有報道將中值濾波、小波變換等圖像處理方法應(yīng)用于對低信噪比的壽命數(shù)據(jù)進行預(yù)降噪處理,以提高熒光壽命分析的精度. 此外, 為了獲取可靠的熒光壽命計算結(jié)果, PA 法處理前通常需要利用系統(tǒng)的儀器響應(yīng)函數(shù)加以校準, 相比無需校準的衰減曲線尾部擬合的方法要多一個步驟, 對非專業(yè)人員而言有時候不夠便利. 但是即便如此, 相信隨著生命科學(xué)對快速定量成像需求的日益增加以及PA 法的進一步發(fā)展完善, phasor-FLIM 一定會有越來越廣泛的應(yīng)用前景, 在生物醫(yī)學(xué)等領(lǐng)域發(fā)揮越來越重要的作用.

    猜你喜歡
    端點光子壽命
    《光子學(xué)報》征稿簡則
    非特征端點條件下PM函數(shù)的迭代根
    人類壽命極限應(yīng)在120~150歲之間
    中老年保健(2021年8期)2021-12-02 23:55:49
    倉鼠的壽命知多少
    不等式求解過程中端點的確定
    馬烈光養(yǎng)生之悟 自靜其心延壽命
    華人時刊(2018年17期)2018-12-07 01:02:20
    人類正常壽命為175歲
    奧秘(2017年12期)2017-07-04 11:37:14
    參數(shù)型Marcinkiewicz積分算子及其交換子的加權(quán)端點估計
    基丁能雖匹配延拓法LMD端點效應(yīng)處理
    在光子帶隙中原子的自發(fā)衰減
    不卡av一区二区三区| 欧美97在线视频| 制服诱惑二区| 国产成人一区二区在线| 人妻人人澡人人爽人人| 考比视频在线观看| 美国免费a级毛片| 女人高潮潮喷娇喘18禁视频| 一二三四中文在线观看免费高清| 日产精品乱码卡一卡2卡三| 亚洲视频免费观看视频| 久久精品国产综合久久久| 国产精品一国产av| 美女高潮到喷水免费观看| 另类精品久久| 日日爽夜夜爽网站| 视频区图区小说| 亚洲国产日韩一区二区| 成年av动漫网址| 午夜91福利影院| 日日撸夜夜添| 国产精品成人在线| 在线观看一区二区三区激情| 五月天丁香电影| 爱豆传媒免费全集在线观看| 97在线人人人人妻| 午夜91福利影院| 99热网站在线观看| 免费观看av网站的网址| 亚洲国产精品999| 免费黄网站久久成人精品| 天堂8中文在线网| 制服诱惑二区| 久久久久国产网址| 国产精品久久久久久久久免| 久久精品国产综合久久久| 一区在线观看完整版| 欧美日韩综合久久久久久| 亚洲精品一区蜜桃| 久久国产亚洲av麻豆专区| 国产欧美日韩综合在线一区二区| 国产高清不卡午夜福利| 日韩制服丝袜自拍偷拍| 精品国产国语对白av| 免费观看a级毛片全部| www.av在线官网国产| 精品99又大又爽又粗少妇毛片| 久久ye,这里只有精品| 国产精品 国内视频| 日韩一区二区视频免费看| 9热在线视频观看99| 国产男女内射视频| 国产精品 国内视频| 免费看不卡的av| 激情视频va一区二区三区| 免费久久久久久久精品成人欧美视频| 欧美日韩一区二区视频在线观看视频在线| 女人被躁到高潮嗷嗷叫费观| 黄色毛片三级朝国网站| 成人国产麻豆网| 亚洲成色77777| 女人高潮潮喷娇喘18禁视频| 一区二区三区乱码不卡18| 一本—道久久a久久精品蜜桃钙片| 大香蕉久久网| 成年av动漫网址| 日韩欧美精品免费久久| 色吧在线观看| 永久免费av网站大全| 午夜免费观看性视频| 最近2019中文字幕mv第一页| 亚洲精品国产一区二区精华液| 99久久中文字幕三级久久日本| 久久精品国产自在天天线| 久久久久久久久久人人人人人人| 亚洲第一av免费看| 婷婷成人精品国产| 建设人人有责人人尽责人人享有的| 一级a爱视频在线免费观看| 日日爽夜夜爽网站| 黑人猛操日本美女一级片| 免费在线观看黄色视频的| 久久久久久久久久久久大奶| 尾随美女入室| 一个人免费看片子| 在线观看国产h片| 人人妻人人澡人人看| 1024视频免费在线观看| 大香蕉久久网| 亚洲图色成人| 亚洲国产欧美网| 极品人妻少妇av视频| 国产免费福利视频在线观看| 国产精品.久久久| 亚洲图色成人| 夫妻午夜视频| 在线观看三级黄色| 视频区图区小说| 26uuu在线亚洲综合色| 在线精品无人区一区二区三| 久久久a久久爽久久v久久| 免费黄色在线免费观看| 高清欧美精品videossex| 精品第一国产精品| 可以免费在线观看a视频的电影网站 | 最近中文字幕高清免费大全6| 亚洲一区二区三区欧美精品| 岛国毛片在线播放| 久久精品国产a三级三级三级| 国产精品99久久99久久久不卡 | 欧美国产精品va在线观看不卡| 免费观看无遮挡的男女| 精品一区二区三区四区五区乱码 | 我的亚洲天堂| 亚洲精品国产色婷婷电影| 啦啦啦在线观看免费高清www| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲一级一片aⅴ在线观看| 国产国语露脸激情在线看| 亚洲色图 男人天堂 中文字幕| 交换朋友夫妻互换小说| 中文字幕人妻丝袜制服| 午夜福利乱码中文字幕| 18禁动态无遮挡网站| 午夜福利影视在线免费观看| 美女福利国产在线| 亚洲精品国产一区二区精华液| 国产极品天堂在线| 国产高清不卡午夜福利| 亚洲在久久综合| 蜜桃在线观看..| 另类精品久久| 91精品伊人久久大香线蕉| 桃花免费在线播放| 国产一区二区三区av在线| 国产在线免费精品| 国产成人精品一,二区| 老司机亚洲免费影院| 久久人人爽av亚洲精品天堂| 中文字幕色久视频| 天天躁夜夜躁狠狠久久av| 欧美人与善性xxx| 在线观看一区二区三区激情| 国产精品嫩草影院av在线观看| 久久影院123| 成人国产av品久久久| 欧美变态另类bdsm刘玥| 免费黄色在线免费观看| 91精品国产国语对白视频| 国产精品三级大全| 日韩不卡一区二区三区视频在线| 一本大道久久a久久精品| √禁漫天堂资源中文www| 最新的欧美精品一区二区| 在线天堂中文资源库| 考比视频在线观看| 成年av动漫网址| 国产精品蜜桃在线观看| 一级a爱视频在线免费观看| 婷婷色综合www| 久久综合国产亚洲精品| 国产熟女午夜一区二区三区| 日产精品乱码卡一卡2卡三| 你懂的网址亚洲精品在线观看| 日本黄色日本黄色录像| 男女无遮挡免费网站观看| 欧美亚洲日本最大视频资源| 免费在线观看视频国产中文字幕亚洲 | 少妇 在线观看| 国产亚洲一区二区精品| 久久精品熟女亚洲av麻豆精品| 亚洲久久久国产精品| 久久久精品免费免费高清| 飞空精品影院首页| 超碰成人久久| 日本爱情动作片www.在线观看| 欧美日韩精品成人综合77777| 天天躁夜夜躁狠狠久久av| 麻豆av在线久日| 街头女战士在线观看网站| 两性夫妻黄色片| 1024视频免费在线观看| 色94色欧美一区二区| 伦精品一区二区三区| 久久精品国产a三级三级三级| 国产精品国产三级国产专区5o| 美女国产高潮福利片在线看| 精品久久久精品久久久| 午夜福利影视在线免费观看| 午夜老司机福利剧场| 国产综合精华液| 熟女少妇亚洲综合色aaa.| 观看av在线不卡| 免费高清在线观看视频在线观看| 在线看a的网站| 欧美精品一区二区大全| 国产高清不卡午夜福利| 亚洲精品美女久久久久99蜜臀 | 校园人妻丝袜中文字幕| 最近最新中文字幕免费大全7| 亚洲欧美成人综合另类久久久| 国产精品人妻久久久影院| 久久久久国产一级毛片高清牌| 日韩av免费高清视频| www.精华液| 久久久久久久久免费视频了| 天堂8中文在线网| 咕卡用的链子| 午夜久久久在线观看| 99热网站在线观看| 一区在线观看完整版| 中文欧美无线码| 成人手机av| 国产精品 国内视频| 最近中文字幕2019免费版| 人人妻人人爽人人添夜夜欢视频| 日本欧美视频一区| 亚洲欧美色中文字幕在线| 中文字幕精品免费在线观看视频| 亚洲欧美清纯卡通| 亚洲色图综合在线观看| 国产亚洲最大av| 成人国语在线视频| 欧美日韩成人在线一区二区| 99国产综合亚洲精品| 丝袜美腿诱惑在线| 中文字幕另类日韩欧美亚洲嫩草| 人人妻人人澡人人爽人人夜夜| 狠狠婷婷综合久久久久久88av| 欧美日韩综合久久久久久| 蜜桃在线观看..| 午夜免费观看性视频| 国产免费又黄又爽又色| 免费看av在线观看网站| 日韩 亚洲 欧美在线| a级毛片黄视频| 亚洲四区av| 国语对白做爰xxxⅹ性视频网站| 久久 成人 亚洲| 男男h啪啪无遮挡| 高清不卡的av网站| 免费观看性生交大片5| 国产精品久久久久久av不卡| 美女国产视频在线观看| 国产成人免费观看mmmm| 美女主播在线视频| 欧美日韩视频精品一区| 成人亚洲欧美一区二区av| 精品亚洲乱码少妇综合久久| 成人亚洲精品一区在线观看| 亚洲中文av在线| 青草久久国产| 亚洲欧洲精品一区二区精品久久久 | 少妇 在线观看| 涩涩av久久男人的天堂| 黑丝袜美女国产一区| 久久99精品国语久久久| xxx大片免费视频| 午夜日韩欧美国产| 精品国产超薄肉色丝袜足j| 国产精品久久久av美女十八| 精品国产一区二区三区久久久樱花| 欧美日韩精品成人综合77777| 中国国产av一级| 97人妻天天添夜夜摸| 久久精品国产鲁丝片午夜精品| 中文字幕人妻丝袜一区二区 | 你懂的网址亚洲精品在线观看| 色婷婷av一区二区三区视频| 亚洲av电影在线观看一区二区三区| 老司机影院成人| 亚洲少妇的诱惑av| 久久鲁丝午夜福利片| 精品人妻偷拍中文字幕| 成年女人毛片免费观看观看9 | 桃花免费在线播放| 我要看黄色一级片免费的| 国产精品熟女久久久久浪| 国产日韩欧美亚洲二区| 99国产综合亚洲精品| 一级毛片 在线播放| 国产av国产精品国产| 热re99久久精品国产66热6| 亚洲第一青青草原| 日本免费在线观看一区| 日韩一卡2卡3卡4卡2021年| 丁香六月天网| xxxhd国产人妻xxx| 亚洲精品日韩在线中文字幕| 免费av中文字幕在线| 久久97久久精品| 另类精品久久| 宅男免费午夜| 亚洲精品,欧美精品| 亚洲精品av麻豆狂野| 免费看av在线观看网站| 免费看不卡的av| 午夜日本视频在线| 久久久久人妻精品一区果冻| 国产黄色视频一区二区在线观看| 午夜影院在线不卡| 欧美激情高清一区二区三区 | 免费观看av网站的网址| 美女国产视频在线观看| 亚洲国产欧美网| 狠狠精品人妻久久久久久综合| 亚洲第一青青草原| 最新的欧美精品一区二区| 男女边吃奶边做爰视频| 女性被躁到高潮视频| 欧美av亚洲av综合av国产av | 亚洲精品美女久久久久99蜜臀 | www.自偷自拍.com| 黄色一级大片看看| 久久久精品国产亚洲av高清涩受| videossex国产| 亚洲精品久久成人aⅴ小说| 91aial.com中文字幕在线观看| 啦啦啦在线免费观看视频4| 18禁观看日本| 国产免费现黄频在线看| 99香蕉大伊视频| 啦啦啦在线免费观看视频4| 国产女主播在线喷水免费视频网站| 久久精品亚洲av国产电影网| 日本爱情动作片www.在线观看| 欧美亚洲日本最大视频资源| 99re6热这里在线精品视频| 欧美日韩一级在线毛片| 男男h啪啪无遮挡| 国产亚洲av片在线观看秒播厂| 男女边摸边吃奶| 天堂8中文在线网| 国产乱人偷精品视频| 亚洲成色77777| 免费不卡的大黄色大毛片视频在线观看| 色婷婷久久久亚洲欧美| 亚洲国产看品久久| 亚洲一码二码三码区别大吗| 国产成人精品婷婷| 曰老女人黄片| 精品人妻在线不人妻| 亚洲av中文av极速乱| 国产精品三级大全| 97精品久久久久久久久久精品| 久久久精品免费免费高清| 日韩av不卡免费在线播放| 国产精品 欧美亚洲| 老司机亚洲免费影院| 女人高潮潮喷娇喘18禁视频| 十八禁网站网址无遮挡| 最新中文字幕久久久久| 久久久久久久精品精品| 国产精品蜜桃在线观看| www.av在线官网国产| 欧美亚洲 丝袜 人妻 在线| 亚洲欧洲日产国产| 春色校园在线视频观看| 国产一区有黄有色的免费视频| 日本免费在线观看一区| 国产一区亚洲一区在线观看| 欧美国产精品va在线观看不卡| 国产精品麻豆人妻色哟哟久久| 夫妻性生交免费视频一级片| 日韩精品有码人妻一区| 欧美日韩国产mv在线观看视频| 亚洲国产精品国产精品| 国产高清不卡午夜福利| 免费不卡的大黄色大毛片视频在线观看| 久久精品久久久久久久性| 人人妻人人爽人人添夜夜欢视频| 制服诱惑二区| 啦啦啦视频在线资源免费观看| 韩国精品一区二区三区| 亚洲成人av在线免费| 国产97色在线日韩免费| 国产在视频线精品| 亚洲第一av免费看| 国产福利在线免费观看视频| 国产亚洲欧美精品永久| 97精品久久久久久久久久精品| 亚洲,欧美精品.| 热re99久久精品国产66热6| 国产亚洲av片在线观看秒播厂| 免费观看av网站的网址| 777久久人妻少妇嫩草av网站| 久久精品熟女亚洲av麻豆精品| 欧美精品人与动牲交sv欧美| 久久97久久精品| 久久久精品免费免费高清| 哪个播放器可以免费观看大片| 性色avwww在线观看| 少妇精品久久久久久久| 亚洲欧美一区二区三区黑人 | 成人二区视频| 99国产综合亚洲精品| 黑人猛操日本美女一级片| 日韩中文字幕欧美一区二区 | 国产av国产精品国产| 亚洲国产成人一精品久久久| 老司机影院毛片| 国产片内射在线| 久久精品国产亚洲av天美| 亚洲国产精品一区二区三区在线| 久久久精品94久久精品| 午夜福利乱码中文字幕| 亚洲一码二码三码区别大吗| av一本久久久久| 在线观看www视频免费| 久久久国产精品麻豆| 国产精品 欧美亚洲| 国产亚洲精品第一综合不卡| 日韩一本色道免费dvd| 啦啦啦在线免费观看视频4| 天天躁狠狠躁夜夜躁狠狠躁| 91午夜精品亚洲一区二区三区| 国产日韩一区二区三区精品不卡| 丝瓜视频免费看黄片| 大话2 男鬼变身卡| 男女高潮啪啪啪动态图| 久久久久久免费高清国产稀缺| 欧美最新免费一区二区三区| 亚洲国产精品一区三区| 天天躁狠狠躁夜夜躁狠狠躁| 免费人妻精品一区二区三区视频| 一本久久精品| 午夜免费男女啪啪视频观看| 亚洲精品,欧美精品| 美女福利国产在线| 欧美激情 高清一区二区三区| 国产男人的电影天堂91| 建设人人有责人人尽责人人享有的| 国产精品一区二区在线观看99| 国产高清国产精品国产三级| 99久久人妻综合| 丝袜人妻中文字幕| 国产成人a∨麻豆精品| 一区福利在线观看| 日韩精品有码人妻一区| 亚洲欧美中文字幕日韩二区| 看免费成人av毛片| 成人国语在线视频| 久久人人爽人人片av| 亚洲av电影在线进入| 亚洲av.av天堂| 亚洲成人手机| 久久精品熟女亚洲av麻豆精品| 国产国语露脸激情在线看| 国产日韩欧美视频二区| 色网站视频免费| 青春草亚洲视频在线观看| 91成人精品电影| 国产精品无大码| 国产日韩欧美视频二区| 寂寞人妻少妇视频99o| xxxhd国产人妻xxx| 男女无遮挡免费网站观看| 免费在线观看黄色视频的| www日本在线高清视频| 久久久久精品人妻al黑| 日韩中文字幕视频在线看片| 精品国产超薄肉色丝袜足j| 亚洲成人av在线免费| 激情五月婷婷亚洲| 纯流量卡能插随身wifi吗| 国产av精品麻豆| 成人影院久久| 亚洲图色成人| 亚洲激情五月婷婷啪啪| 一区二区三区乱码不卡18| 老汉色av国产亚洲站长工具| 国产亚洲最大av| 久久久欧美国产精品| 高清av免费在线| 亚洲美女搞黄在线观看| 91午夜精品亚洲一区二区三区| 精品国产国语对白av| 国产毛片在线视频| 成人国产av品久久久| 亚洲国产精品一区二区三区在线| 中国国产av一级| 在线观看一区二区三区激情| 亚洲色图 男人天堂 中文字幕| 久久久久网色| 国产成人精品在线电影| 成人18禁高潮啪啪吃奶动态图| 波多野结衣av一区二区av| 精品久久久久久电影网| 久久久久久人妻| 免费观看在线日韩| 国产片内射在线| 久久精品国产a三级三级三级| 狠狠婷婷综合久久久久久88av| 91在线精品国自产拍蜜月| 亚洲精品一区蜜桃| 亚洲av欧美aⅴ国产| 我的亚洲天堂| 国产亚洲av片在线观看秒播厂| 国产精品一国产av| 国产在视频线精品| kizo精华| 大片电影免费在线观看免费| 麻豆精品久久久久久蜜桃| 国产精品一区二区在线观看99| 国产亚洲午夜精品一区二区久久| 丝袜在线中文字幕| 国产精品偷伦视频观看了| 嫩草影院入口| 两个人免费观看高清视频| 国产欧美日韩综合在线一区二区| 最近手机中文字幕大全| 五月开心婷婷网| 国产精品免费大片| 久久久久久人人人人人| 中文字幕人妻丝袜制服| 美女xxoo啪啪120秒动态图| 成年女人在线观看亚洲视频| 黄色配什么色好看| 久久久久视频综合| 女人高潮潮喷娇喘18禁视频| 久久免费观看电影| 又黄又粗又硬又大视频| 欧美激情 高清一区二区三区| 七月丁香在线播放| 999精品在线视频| 中文字幕人妻丝袜制服| 啦啦啦视频在线资源免费观看| h视频一区二区三区| 一区二区av电影网| 一边摸一边做爽爽视频免费| 一级片免费观看大全| 亚洲欧美一区二区三区国产| 男的添女的下面高潮视频| 美国免费a级毛片| 免费高清在线观看日韩| 久久韩国三级中文字幕| 国产毛片在线视频| 国产亚洲最大av| 亚洲精品久久成人aⅴ小说| 日本猛色少妇xxxxx猛交久久| 午夜影院在线不卡| 国产男女超爽视频在线观看| 日本爱情动作片www.在线观看| 熟女av电影| 搡老乐熟女国产| 18禁国产床啪视频网站| 一级毛片电影观看| 精品国产一区二区三区四区第35| 一区二区av电影网| 丝瓜视频免费看黄片| 精品一区二区免费观看| 久久人人97超碰香蕉20202| 日韩精品有码人妻一区| 亚洲美女搞黄在线观看| 亚洲综合色惰| 老女人水多毛片| 少妇的丰满在线观看| 久久久久久久亚洲中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 精品一品国产午夜福利视频| av一本久久久久| 少妇被粗大的猛进出69影院| 亚洲国产精品一区三区| 亚洲,欧美精品.| 日日爽夜夜爽网站| 国产av码专区亚洲av| 一级毛片我不卡| 亚洲经典国产精华液单| 久久97久久精品| 亚洲成色77777| 人人妻人人添人人爽欧美一区卜| 精品视频人人做人人爽| 国产男人的电影天堂91| 久久精品aⅴ一区二区三区四区 | 亚洲综合精品二区| 亚洲av日韩在线播放| 久久久久国产精品人妻一区二区| 2022亚洲国产成人精品| 亚洲成人一二三区av| 少妇被粗大猛烈的视频| 自线自在国产av| 国产成人免费观看mmmm| 免费少妇av软件| 欧美少妇被猛烈插入视频| 少妇的逼水好多| 两个人看的免费小视频| 亚洲在久久综合| 亚洲,一卡二卡三卡| 国产成人精品在线电影| 午夜福利影视在线免费观看| 亚洲av日韩在线播放| 久久精品夜色国产| 一区在线观看完整版| 在线观看人妻少妇| 精品一区在线观看国产| 日韩,欧美,国产一区二区三区| 国产免费一区二区三区四区乱码| 免费女性裸体啪啪无遮挡网站| 女性被躁到高潮视频| 亚洲av国产av综合av卡| 老汉色∧v一级毛片| 亚洲av成人精品一二三区| 日本91视频免费播放| 亚洲精品国产一区二区精华液| 国产精品国产av在线观看| 在线观看美女被高潮喷水网站| 大香蕉久久网| 我要看黄色一级片免费的| 母亲3免费完整高清在线观看 | 丝袜喷水一区| 高清av免费在线| 免费看不卡的av| av国产精品久久久久影院| 欧美最新免费一区二区三区| 亚洲欧美清纯卡通|