杜慶華，李慶山，許艷麗

(廣州醫(yī)科大學(xué)附屬?gòu)V州市第一人民醫(yī)院血液內(nèi)科，廣東 廣州 510180)

四種坐標(biāo)函數(shù)對(duì)流式細(xì)胞術(shù)數(shù)據(jù)可視化的影響

杜慶華，李慶山，許艷麗

(廣州醫(yī)科大學(xué)附屬?gòu)V州市第一人民醫(yī)院血液內(nèi)科，廣東 廣州 510180)

目的 研究流式細(xì)胞術(shù)數(shù)據(jù)可視化中四種常見坐標(biāo)函數(shù)的特點(diǎn)，探討這四種坐標(biāo)函數(shù)在數(shù)據(jù)可視化中的應(yīng)用。方法對(duì)比分析線性、對(duì)數(shù)、HyperLog以及Logical四個(gè)坐標(biāo)函數(shù)曲線的特點(diǎn)，使用不同坐標(biāo)函數(shù)對(duì)同一個(gè)數(shù)據(jù)使用散點(diǎn)圖進(jìn)行顯示，對(duì)比圖形顯示的差異。結(jié)果線性函數(shù)對(duì)數(shù)據(jù)呈等比例顯示，但動(dòng)態(tài)范圍不足。對(duì)數(shù)函數(shù)能拉伸1附近的數(shù)值，壓縮顯示數(shù)值大的數(shù)值，但不能顯示小于1的值。HyperLog與Logical函數(shù)數(shù)值小的時(shí)候近似線性顯示，數(shù)值大時(shí)近似對(duì)數(shù)顯示。結(jié)論不同的坐標(biāo)函數(shù)有不同的特點(diǎn)，在進(jìn)行流式細(xì)胞術(shù)數(shù)據(jù)分析時(shí)，要根據(jù)需要選擇相應(yīng)的坐標(biāo)函數(shù)進(jìn)行顯示。

流式細(xì)胞術(shù)；坐標(biāo)函數(shù)；數(shù)據(jù)可視化

數(shù)據(jù)可視化是流式細(xì)胞術(shù)數(shù)據(jù)分析中的一個(gè)重要環(huán)節(jié)，數(shù)據(jù)通過圖形顯示，然后才能通過設(shè)門操作進(jìn)行進(jìn)一步分析。流式數(shù)據(jù)的顯示樣式有直方圖、散點(diǎn)圖、等高圖或者密度圖等。每個(gè)細(xì)胞的熒光強(qiáng)度數(shù)據(jù)要在圖中顯示均需要通過坐標(biāo)函數(shù)計(jì)算其位置，因此坐標(biāo)函數(shù)的選擇會(huì)直接影響圖形的形狀。流式細(xì)胞術(shù)中常用的坐標(biāo)有線性坐標(biāo)、對(duì)數(shù)坐標(biāo)以及其他對(duì)數(shù)衍生的坐標(biāo)。由于細(xì)胞表面及內(nèi)部分子表達(dá)的差異很大，這些數(shù)據(jù)很難通過線性坐標(biāo)顯示[1-2]，因此一直以來使用流式細(xì)胞術(shù)進(jìn)行免疫表型分析基本都使用對(duì)數(shù)坐標(biāo)以顯示數(shù)據(jù)。但因?yàn)閷?duì)數(shù)無法恰當(dāng)?shù)仫@示小于1的數(shù)據(jù)，故Bagwell及Parks建議使用對(duì)數(shù)衍生函數(shù)進(jìn)行坐標(biāo)轉(zhuǎn)換[3-4]。而這些坐標(biāo)函數(shù)間特性的比較暫時(shí)未見詳細(xì)報(bào)道。本文將深入探討不同坐標(biāo)函數(shù)的特點(diǎn)，為坐標(biāo)函數(shù)的選擇提供指引。

1 材料與方法

1.1 儀器及軟件 流式細(xì)胞儀為FACS Canto (美國(guó)BD公司)，數(shù)據(jù)為常規(guī)淋巴細(xì)胞亞群臨床檢測(cè)的流式數(shù)據(jù)，軟件使用本實(shí)驗(yàn)室自行研發(fā)的流式數(shù)據(jù)分析軟件CFCS(軟件著作權(quán)登記號(hào)：2010SR064983)。

1.2 方法

1.2.1 函數(shù)曲線的繪制 對(duì)數(shù)函數(shù)及其反函數(shù)為：

其中r為數(shù)據(jù)的分辨率，d為動(dòng)態(tài)范圍的數(shù)量級(jí)。HypherLog為隱函數(shù)[3]，故僅能寫出其反函數(shù)：

因此在進(jìn)行HypherLog函數(shù)運(yùn)算時(shí)，必須通過其反函數(shù)進(jìn)行求根運(yùn)算。函數(shù)中b為相關(guān)系數(shù)，是影響線性范圍與對(duì)數(shù)范圍的參數(shù)。當(dāng)b=0時(shí)曲線最接近對(duì)數(shù)[3]。

Logical與HypherLog同樣是隱函數(shù)[4-5]，其反函數(shù)為：

其中M為圖形顯示寬寬的數(shù)量級(jí)，而W=2p log (p)/(p+1)，為線性范圍的寬度，其值可以根據(jù)以下公式求出：

其中n為負(fù)值范圍參考點(diǎn)的值。

1.2.2 坐標(biāo)函數(shù)曲線的繪制 使用MATLAB7.0繪制線性、對(duì)數(shù)函數(shù)、HyperLog函數(shù)以及Logical函數(shù)的函數(shù)曲線。

1.2.3 不同坐標(biāo)函數(shù)對(duì)圖形的影響 同一個(gè)數(shù)據(jù)在CFCS上分別使用線性坐標(biāo)、對(duì)數(shù)坐標(biāo)、HyperLog以及Logical顯示，并對(duì)比四者圖形的區(qū)別。

2 結(jié)果

2.1 不同的坐標(biāo)函數(shù)曲線的特點(diǎn) 從函數(shù)曲線發(fā)現(xiàn)線性坐標(biāo)能顯示整個(gè)實(shí)數(shù)域的數(shù)據(jù)，且具有對(duì)稱性。因其等比例的特性，當(dāng)數(shù)據(jù)動(dòng)態(tài)范圍較大時(shí)會(huì)使數(shù)值小的數(shù)據(jù)被壓縮，造成顯示效果不佳。對(duì)數(shù)坐標(biāo)能保證數(shù)值小數(shù)據(jù)的顯示空間，但其缺點(diǎn)是僅僅能顯示≥1的數(shù)據(jù)，對(duì)于＜1的數(shù)值全按0來處理。HyperLog與Logical能顯示實(shí)數(shù)域的數(shù)據(jù)，且二者在數(shù)值較小的時(shí)候函數(shù)曲線接近線性，在數(shù)值較大的時(shí)候接近對(duì)數(shù)曲線(見圖1)。

圖1 線性、對(duì)數(shù)、HyperLog以及Logical的函數(shù)曲線

2.2 不同的坐標(biāo)函數(shù)對(duì)流式數(shù)據(jù)顯示的影響 通過使用不同坐標(biāo)顯示相同的免疫熒光數(shù)據(jù)，發(fā)現(xiàn)線性坐標(biāo)顯示時(shí)細(xì)胞群體往往聚集在數(shù)軸的一端，難以對(duì)群體進(jìn)行區(qū)分。對(duì)數(shù)坐標(biāo)顯示免疫熒光實(shí)驗(yàn)的數(shù)據(jù)時(shí)，細(xì)胞群體分布比較清晰，但是進(jìn)行多色熒光分析進(jìn)行補(bǔ)償后，有8%的細(xì)胞堆疊在基線上，造成視覺誤差，容易對(duì)分析結(jié)果造成影響。Logical與HypherLog兩者顯示效果無顯著差異，均能很好地控制堆疊在基線上細(xì)胞的數(shù)量，細(xì)胞群體分布清晰，補(bǔ)償后群體離散程度較對(duì)數(shù)坐標(biāo)明顯(見圖2)。

圖2 使用4種坐標(biāo)函數(shù)對(duì)同一個(gè)數(shù)據(jù)生成散點(diǎn)圖的比較

3 討論

線性坐標(biāo)就是把數(shù)據(jù)直接按比例顯示在圖形上，其優(yōu)點(diǎn)是簡(jiǎn)單直接，能顯示實(shí)數(shù)域的所有數(shù)據(jù)，一般用于DNA含量的的測(cè)定，如細(xì)胞周期分析。但細(xì)胞表面分子表達(dá)情況差異相當(dāng)大，需要極寬的動(dòng)態(tài)范圍才能顯示這類型的數(shù)據(jù)，而且這類型的數(shù)據(jù)低熒光強(qiáng)度群體離散程度小并近似正態(tài)分布，高熒光強(qiáng)度群體離散程度大接近對(duì)數(shù)正態(tài)分布。線性坐標(biāo)顯示動(dòng)態(tài)范圍寬的數(shù)據(jù)時(shí)，會(huì)壓縮數(shù)值小數(shù)據(jù)的顯示范圍，造成細(xì)胞群體往往聚集在數(shù)軸的一端。而對(duì)數(shù)坐標(biāo)能拉伸數(shù)值小數(shù)據(jù)的顯示范圍，壓縮數(shù)值大的數(shù)據(jù)，擁有較寬的動(dòng)態(tài)范圍，因此對(duì)數(shù)坐標(biāo)比線性坐標(biāo)更適合顯示免疫熒光的數(shù)據(jù)。但是對(duì)數(shù)坐標(biāo)的缺陷是無法對(duì)＜1的值進(jìn)行變換，因此一般對(duì)于＜1的值轉(zhuǎn)換為0。

在實(shí)際應(yīng)用中為了消除通道間的熒光滲漏，往往需要對(duì)數(shù)據(jù)進(jìn)行補(bǔ)償運(yùn)算，該運(yùn)算是多個(gè)通道間的減法運(yùn)算。運(yùn)算過程中會(huì)把某一通道的檢測(cè)誤差引入到其他通道中，因此會(huì)造成數(shù)據(jù)離散增大，如果細(xì)胞群體被影響通道的熒光強(qiáng)度較低時(shí)，數(shù)據(jù)的離散將造成＜1甚者負(fù)值數(shù)據(jù)的產(chǎn)生，這些數(shù)據(jù)在對(duì)數(shù)坐標(biāo)中堆疊在基線下，使該處細(xì)胞密度異常增高，造成假群體的產(chǎn)生[4，6]。從圖2中我們可看到8%的細(xì)胞堆在x坐標(biāo)上。如果設(shè)門分析時(shí)忽略了這些細(xì)胞，將對(duì)結(jié)果造成一定的影響。為了克服對(duì)數(shù)坐標(biāo)的這些缺陷，Parks提出理想的坐標(biāo)函數(shù)應(yīng)具備以下特點(diǎn)：(1)函數(shù)可對(duì)不同的數(shù)據(jù)進(jìn)行顯示的優(yōu)化；(2)為了提供足夠?qū)挼膭?dòng)態(tài)范圍，更好地顯示高熒光強(qiáng)度的對(duì)數(shù)正態(tài)分布數(shù)據(jù)，函數(shù)應(yīng)在隨數(shù)值增大逼近對(duì)數(shù)曲線；(3)函數(shù)應(yīng)在接近0的區(qū)域?qū)ΨQ并近似線性分布，這樣更適合顯示低熒光強(qiáng)度正態(tài)分布的數(shù)據(jù)；(4)線性區(qū)域應(yīng)盡可能平滑地向?qū)?shù)區(qū)域過渡，以避免扭曲地顯示數(shù)據(jù)；(5)線性化強(qiáng)度應(yīng)隨線性范圍增大而增大[4]。

Bagwell及Parks分別提出了HypherLog與Logical函數(shù)。這兩個(gè)函數(shù)的特點(diǎn)是能顯示負(fù)值的數(shù)據(jù)，且小數(shù)值區(qū)域接近線性分布，高數(shù)值區(qū)接近對(duì)數(shù)分布，且線性區(qū)與對(duì)數(shù)區(qū)能平滑過渡。在實(shí)際顯示中通過HypherLog及Logical兩個(gè)坐標(biāo)與對(duì)數(shù)坐標(biāo)比較，發(fā)現(xiàn)對(duì)于熒光強(qiáng)度大的細(xì)胞群體差異不大。HypherLog及Logical兩個(gè)坐標(biāo)下，熒光強(qiáng)度低的群體比對(duì)數(shù)坐標(biāo)更趨于聚集，更重要的是這兩個(gè)坐標(biāo)下能顯示對(duì)數(shù)坐標(biāo)不能顯示的數(shù)值＜1的值。此外，HypherLog及Logical能識(shí)別一些因過補(bǔ)償而造成細(xì)胞群體落在坐標(biāo)軸上的數(shù)據(jù)，并將之顯示出來。不同數(shù)據(jù)＜0的值往往是不一樣的，故HypherLog與Logical需要適當(dāng)調(diào)整參數(shù)使得數(shù)據(jù)能恰當(dāng)?shù)仫@示。HypherLog通過相關(guān)系數(shù)b來控制線性區(qū)的斜率，b越大線性區(qū)域越寬，b為0時(shí)曲線最接近對(duì)數(shù)曲線。但Bagwell發(fā)現(xiàn)b太小會(huì)使陰性區(qū)域的原來為一群的群體在視覺上分裂為兩個(gè)，因此使用時(shí)要注意b的取值。Logical使用負(fù)值大小作為參數(shù)對(duì)曲線的先行范圍進(jìn)行控制，通過公式5計(jì)算W的值，-W與W之間即為線性區(qū)域。

HypherLog的二階導(dǎo)數(shù)不為0而Logical為0，因此在0附近Logical比HypherLog更接近線性，且更快靠近對(duì)數(shù)曲線。從圖2實(shí)際顯示中看，這差異并不會(huì)對(duì)數(shù)據(jù)顯示分析造成很大的影響。

由于不同的參數(shù)會(huì)對(duì)顯示效果造成一定影響，因此HypherLog與Logical兩個(gè)函數(shù)的穩(wěn)定性不如對(duì)數(shù)。當(dāng)進(jìn)行表達(dá)模式分析(如白血病免疫學(xué)分型)的時(shí)候，細(xì)胞群體的位置及走向是極其重要的信息。當(dāng)數(shù)據(jù)中含有個(gè)別極端數(shù)值時(shí)，要將其恰當(dāng)顯示則需要調(diào)整HypherLog及Logical函數(shù)的參數(shù)，細(xì)胞群體的形狀及位置可能因此而發(fā)生扭曲改變(特別是熒光強(qiáng)度低的群體)，容易導(dǎo)致分析錯(cuò)誤。

對(duì)數(shù)函數(shù)能直接運(yùn)算得到結(jié)果，而HypherLog與Logical無法顯式寫出其函數(shù)形式而無法直接進(jìn)行運(yùn)算，因此必須對(duì)方程進(jìn)行求根運(yùn)算。在分析軟件中，一般使用牛頓迭代法進(jìn)行求根運(yùn)算，運(yùn)算結(jié)果的精度隨迭代次數(shù)增加而增加。HypherLog與Logical二者運(yùn)算效率遠(yuǎn)低于對(duì)數(shù)運(yùn)算。Logical運(yùn)算的復(fù)雜程度大于HypherLog，故其運(yùn)算效率最低。在進(jìn)行大數(shù)據(jù)量的分析時(shí)我們可以預(yù)先建立一個(gè)函數(shù)表，進(jìn)行運(yùn)算時(shí)可通過二分查找法進(jìn)行檢索，從而極大地提高運(yùn)算的效率。

變異范圍小和正態(tài)分布的數(shù)據(jù)(如細(xì)胞周期分析)宜使用線性坐標(biāo)顯示。免疫熒光這類動(dòng)態(tài)范圍大的數(shù)據(jù)不宜使用線性顯示。如果這類數(shù)據(jù)無需考慮群體形狀走勢(shì)等信息，僅僅進(jìn)行群體劃分，則直接使用HypherLog或Logical進(jìn)行顯示設(shè)門即可。如果需要進(jìn)行表達(dá)模式分析，建議先使用對(duì)數(shù)顯示以了解群體分布情況，再使用HypherLog或Logical進(jìn)行設(shè)門，以避免負(fù)值數(shù)據(jù)的丟失?？傊煌淖鴺?biāo)函數(shù)有其不同的特性，對(duì)數(shù)據(jù)的顯示及分析有著重要的影響。因此我們必須對(duì)不同坐標(biāo)函數(shù)的特性有所了解，分析時(shí)使用不同的坐標(biāo)函數(shù)顯示數(shù)據(jù)，才能對(duì)數(shù)據(jù)有更全面的認(rèn)識(shí)，以避免坐標(biāo)函數(shù)選擇不當(dāng)而造成的視覺誤差，進(jìn)而影響分析結(jié)果。

[1]Muirhead KA,Schmitt TC,Muirhead AR.Muirhead,determination of linear fluorescence intensities from flow cytometric data accumulated with logarithmic amplifiers[J].Cytometry,1983,3(4):251-256.

[2]Watson JV,Chambers SH,Smith PJ.A pragmatic approach to the analysis of DNA histograms with a definable G1 peak[J].Cytometry,1987,8(1):1-8.

[3]Bagwell CB.Hyperlog-a flexible log-like transform for negative,zero,and positive valued data[J].CytometryA,2005,64(1):34-42.

[4]Parks DR,Roederer M,Moore WA.A new"Logicle"display method avoids deceptive effects of logarithmic scaling for low signals and compensated data[J].CytometryA,2006,69(6):541-551.

[5]Moore WA,Parks DR.Update for the logicle data scale including operational code implementations[J].Cytometry A,2012,81(4): 273-277.

[6]Finak G,Perez JM,Weng A,et al.Optimizing transformations for automated,high throughput analysis of flow cytometry data[J]. BMC Bioinformatics,2010,11:546.

Effects of four scaling functions on flow cytometry data visualization.

DU Qing-hua,LI Qing-shan,XU Yan-li.
Department of Hematology,Guangzhou First People's Hospital,Guangzhou Medical University,Guangzhou 510180, Guangdong,CHINA

ObjectiveTo explore the features of four scaling functions,and investigate the selection of 4 scaling functions in data visualization for flow cytometry.MethodsWe compared the features of 4 function curves, and dot plots from a same data file in 4 different scaling functions.ResultsLinear function displays data in a manner of proportion,but it is insufficient in dynamic range.Logarithmic function can stretch displaying the data value near zero,and can compress displaying the large data value,but it can not display data value less than zero.HyperLog and Logical function can display small data value similar to linear distribution,and display large data value similar to near logarithm.ConclusionThere are different features in different scaling functions.We should choose correct scaling function to display in data analyzing of flow cytometry.

Flow cytometry;Scaling function;Data visualization

R446

A

1003—6350（2015）15—2259—03

10.3969/j.issn.1003-6350.2015.15.0814

2015-03-02)

2013年廣東省自然科學(xué)基金(編號(hào)：S2013010016726)

李慶山。E-mail:qingshanli@hotmail.com