讓ImageJ在Western Blot圖像分析中綻放-2

        《在Western Blot條帶定量分析中圖像灰度值的應(yīng)用實(shí)測-下》一文提到，用ImageJ軟件對圖像定量分析流程中， 分析選區(qū)（ROI）的合理設(shè)置是基礎(chǔ)。但用Image J軟件為Western Blot條帶定量分析設(shè)置ROI的操作方法、操作質(zhì)控和驗(yàn)證的標(biāo)準(zhǔn)，Image J用戶指南并無介紹。
       他山之石，可以攻玉。我們不煩先了解一下，Bio-Rad、Alapha和UVP等專業(yè)凝膠圖像分析軟件操作指南中，有關(guān)電泳凝膠圖像和Western Blot條帶分析選區(qū)ROI設(shè)置經(jīng)驗(yàn)做法。

一、專業(yè)圖像分析軟件關(guān)于ROI創(chuàng)建的指導(dǎo)要求

1. 1 AlphaView Software-講究松緊適度

       Alapha公司的AlphaView Software 在用戶指南【Chapter 5:Image Analysis Tools 的Multiplex Band Analysis】指出，ROI邊框要充分圈選全部待分析對象，實(shí)現(xiàn)“應(yīng)檢盡檢”。

       采用“一帶一景”（Local Background）“因區(qū)施策”法計(jì)算條帶背景時(shí)，將選區(qū)中強(qiáng)度值最低的10個(gè)像素作為采樣背景取其均值作為該條帶背景值（Background region applied using the average of the 10 lowest pixels in a band region）。這就要求創(chuàng)建的ROI要略大于被檢測條帶（slightly larger than the bands），確保選框內(nèi)有足夠數(shù)量（不應(yīng)低于10個(gè)）的條帶背景像素，同時(shí)，避免引入過多背景像素和摻入鄰近條帶像素。

       正如AlphaView指南給出的object示范那樣，ROI上、下、左、右四個(gè)方位選框內(nèi)側(cè)，都應(yīng)該包含不低于1個(gè)像素寬度的背景。

1. 2 Bio-Rad Quantity One-力求完整包容

       Bio-Rad公司Quantity One在指南中對分析選區(qū)ROI創(chuàng)建要求是，繪制的選區(qū)要囊括全部含有數(shù)據(jù)信息的檢測像素。針對背景不均勻的圖像，采用Local background背景扣除法時(shí)，選框內(nèi)應(yīng)包含一圈1個(gè)像素寬但結(jié)構(gòu)連續(xù)、完整的圖像背景地帶。選區(qū)最外圍1圈像素強(qiáng)度值（intensities of the pixels in a 1-pixel border around the volume）之和除以邊框像素總數(shù)（total number of border pixels）所得均值作為該選區(qū)的背景值。

       換句話說，如緊貼ROI選框的像素背景像素缺失或誤入過多強(qiáng)悍的背景干擾像素，軟件依然將外圍一圈像素視作背景，據(jù)此產(chǎn)生背景強(qiáng)度將升高，必會(huì)造成檢測信號(hào)損失而拉低選區(qū)的檢測值。

1. 3 UVP LS Software-安內(nèi)攘外并舉

       UVP公司的LS Software的Define Area Density Background操作，其指南在【Define Background】中指出，在針對復(fù)雜不均勻的Western Blot條帶圖像進(jìn)行Area Density分析時(shí)，軟件默認(rèn)每個(gè)ROI選區(qū)都會(huì)自動(dòng)定義一個(gè)VIP專屬背景（separate background）。這個(gè)用于自動(dòng)計(jì)算圖像背景值的選區(qū)，位于標(biāo)記選區(qū)外，其寬為選區(qū)邊框外緊鄰的三個(gè)像素（three pixels wide），長為三個(gè)像素湊成的狹長環(huán)形區(qū)帶的周長（the total of the perimeter around the region marked）。

       概括起來即，ROI內(nèi)全部像素為檢測目標(biāo)像素，ROI選框外推三個(gè)像素的環(huán)形區(qū)域視為背景帶。 

1. 4 小結(jié)

       以上調(diào)研表明，Bio-Rad、Alapha和UVP等知名凝膠圖像分析軟件在條帶定量分析中，對分析選區(qū)ROI的設(shè)置規(guī)則有所不同。其區(qū)別的關(guān)鍵點(diǎn)在于，分析選區(qū)內(nèi)是否需要引入圖像的背景以及對背景像素范圍的控制。

1）選區(qū)內(nèi)背景分析法

       這種分析模式的代表是AlphaView Software和Bio-Rad公司的Quantity One軟件。

       所建立ROI的最外圍一圈須是背景像素，相當(dāng)于采用有效像素區(qū)外延伸一圈像素，作為該條帶背景的計(jì)算依據(jù)。

       此模式產(chǎn)生的選區(qū)背景的特征是：選區(qū)內(nèi)所有像素均被扣除相同數(shù)量的背景值。其背景去除效果，類似于Photoshop的圖像亮度整體調(diào)節(jié)功能。

       適用對象為：Western Blot實(shí)驗(yàn)圖像上條帶的背景均衡、均勻，且信號(hào)峰前、后均背景本底條件一致或十分接近。

2）選區(qū)外背景分析法

       以UVP品牌LS Software軟件為代表的這種分析模式，分析選區(qū)內(nèi)部全部視為有效信號(hào)像素，從選區(qū)邊框向外延伸3圈像素寬度作為條帶參考背景。

       這種模式產(chǎn)生的選區(qū)背景適用對象與前一種分析模式相同。

      與前一種選區(qū)設(shè)置要求的區(qū)別是，它采用ROI外部的3圈像素的均值計(jì)算背景，而前者是ROI內(nèi)部最外圍1圈像素作為背景。       

       現(xiàn)在的問題是，ImageJ軟件在用于Western Blot條帶定量分析過程中，ROI創(chuàng)建應(yīng)該遵循什么規(guī)則？

       要回答這個(gè)問題，首先要了解ImageJ軟件中ROI的特征。而觀察和分析ROI，最直觀有效的方法之一，就是下面要討論的Plot Profile分析模塊。

二、ImageJ軟件Plot Profile分析模塊功能特點(diǎn)
2.1 Plot Profile模塊工作界面 

       在ImageJ軟件中執(zhí)行【Analyze】→【Plot Profile】命令，即可調(diào)出Plot Profile分析模塊，顯示直線選區(qū)（line selections）或矩形選區(qū) (rectangular selections)的像素強(qiáng)度-水平距離關(guān)系二維圖形。  

       X軸代表垂直方向像素列離開選區(qū)左側(cè)邊界的水平方向距離（distance/pixels），即水平方向跨越的像素?cái)?shù)量（horizontal distance through the selection），用像素?cái)?shù)表示。
       Y軸代表選區(qū)內(nèi)垂直方向像素列的平均像素強(qiáng)度值（vertically averaged pixel intensity）。本示例為Western Blot暗底白帶型效果圖，列強(qiáng)度值為像素灰度值。

       對于矩形選區(qū)(rectangular selections)，當(dāng)選區(qū)為2個(gè)像素或像素列時(shí)顯示為斜線圖；選區(qū)像素列達(dá)到3個(gè)像素顯示為峰線曲線圖。
       用矩形選區(qū)可選定當(dāng)行像素顯示單行或多行像素的峰線曲線走勢圖。在ImageJ軟件為Western Blot實(shí)驗(yàn)條帶創(chuàng)建分析選區(qū)ROI時(shí)，這個(gè)功能可用于預(yù)判斷條帶邊界范圍。

       當(dāng)選區(qū)為1個(gè)像素列時(shí)，在窗口左下角顯示的數(shù)值為該列像素的強(qiáng)度均值；

       當(dāng)選區(qū)僅包含1個(gè)像素時(shí)，顯示的是該像素灰度值。
       點(diǎn)擊菜單中【Analyze】→【Measure】命令，在彈出的【Results】窗口中，顯示的是黃色選框內(nèi)從像素列0到像素列16在內(nèi)全部像素的統(tǒng)計(jì)分析結(jié)果。

2.2 Plot Profile實(shí)時(shí)顯示功能 

       點(diǎn)擊Plot Profile窗口底部的Live按鈕，字體可在紅色-黑色之間隨意切換。

       當(dāng)Live呈現(xiàn)為紅色字體狀態(tài)時(shí)，當(dāng)前活動(dòng)圖像窗口選區(qū)的調(diào)整變化，窗口的選區(qū)像素強(qiáng)度值峰線曲線圖隨之發(fā)生改變。這一功能，對于創(chuàng)建Western Blot條帶分析選區(qū)ROI，具有特別重要的參考價(jià)值。

       此外，在后續(xù)對選區(qū)背景扣除條件設(shè)置方案的評估中，可實(shí)時(shí)預(yù)覽不同背景扣除條件的去背景效果，有利于背景扣除條件的評估優(yōu)化。 

三、運(yùn)用ImageJ軟件為Western Blot實(shí)驗(yàn)條帶建立選區(qū)ROI的基本原理和規(guī)則 

3.1 用ImageJ軟件為Western Blot條帶建立分析選區(qū)的基本原理

       一般認(rèn)為，CCD成像系統(tǒng)生成Western Blot圖像強(qiáng)度信號(hào)，包含三種成份，即CCD成像芯片本身的電子噪音，來自實(shí)驗(yàn)樣品（如轉(zhuǎn)印PVDF膜、化學(xué)試劑和緩沖液的自發(fā)熒光）的光子信號(hào)和條帶發(fā)光產(chǎn)生的檢測光子信號(hào)。對在成像過程，CCD本底噪音電流效應(yīng)、樣品背景光學(xué)能量和連同條帶釋放的光學(xué)能量都被如實(shí)記錄。CCD對感受電子信號(hào)不問出處，來者不拒，是無差別累積存儲(chǔ)。因此，圖像中的像素強(qiáng)度值是由圖像背景、檢測信號(hào)和CCD本底噪音三者相生相伴、相互疊加形成的。圖像的背景與檢測有效信號(hào)這兩股勢力之間，并非此消彼長的博弈和競爭的關(guān)系。
       而ImageJ軟件的Plot Profile分析模塊顯示的圖像信號(hào)分布-位置關(guān)系圖中，任一像素列的強(qiáng)度數(shù)值，同樣是背景-信號(hào)的綜合值。
       眾所周知，WB條帶的蛋白含量有其特殊分布規(guī)律：沿著條帶電泳方向，離開條帶短軸中軸線越遠(yuǎn)，條帶含量越低直至最終消失。因此，在圖像短軸方向上，檢測信號(hào)與背景信號(hào)疊加信號(hào)強(qiáng)度也相應(yīng)響應(yīng)衰減趨勢。表現(xiàn)在Plot的圖像信號(hào)強(qiáng)度分布圖上，是信號(hào)分布曲線隨著條帶信號(hào)分布，呈現(xiàn)出低谷-峰值-衰減-低谷的規(guī)律性變化。
       依據(jù)圖像某個(gè)區(qū)域信號(hào)強(qiáng)度曲線的顯著變化，即可判定此處條帶有效檢測信號(hào)高低和有無，據(jù)此評估條帶有效信號(hào)的分布范圍。

3.2 用Plot Profile分析模塊創(chuàng)建WB條帶ROI的基本規(guī)則

       首先，重視利用Plot Profile模塊分析功能，以圖像信號(hào)曲線走勢，指導(dǎo)研判條帶有效信號(hào)分布范圍；

       其次，Measure分析的是ROI內(nèi)全部有效像素的分析統(tǒng)計(jì)值。因此，ROI須最大程度地囊括條帶的全部有效信號(hào)像素，同時(shí)限制和盡量減少背景像素進(jìn)入選區(qū)。ImageJ圖像分析軟件創(chuàng)建的分析選區(qū)ROI不僅不需要引入背景像素做參考，想法，還應(yīng)避免框選代表圖像背景的像素。這一點(diǎn)，與其他專業(yè)凝膠圖像分析軟件的選區(qū)設(shè)置規(guī)則，存在根本區(qū)別 。

       再次，ROI設(shè)置合理與否判定標(biāo)準(zhǔn)之一是，背景扣除后，ROI統(tǒng)計(jì)值中min值應(yīng)為0。min值大于0，則需適度調(diào)整選區(qū)位置和范圍，重新計(jì)算滾球半徑。

       創(chuàng)建ROI的具體流程，可參考《如何用ImageJ為Western Blot條帶定量建立分析選區(qū)》。

用參考文獻(xiàn)

1. Bio-Rad Quantity One User Guide for Version 4.4

2. AlphaView Software User Guide_94-13934-00 Revision 3a (2011)

3. Life Science Software Installation and User Instructions (81-0254-01 Rev L)

4. ImageJ User Guide (IJ 1.46r)