8月13日，國際學術期刊Cell Discovery在線發(fā)表了中國科學院上海營養(yǎng)與健康研究所計算生物學重點實驗室（馬普計算生物學研究所）邵振研究組研究論文“MAP: model-based analysis of proteomic data to detect proteins with significant abundance changes”，報道了一種新計算模型MAP，用于統(tǒng)計分析基于同位素標記產(chǎn)生的定量蛋白質組數(shù)據(jù)并鑒定其中差異表達的蛋白質。
基于同位素標記和質譜技術的定量蛋白質組實驗（如iTRAQ、TMT和SILAC等）能同時檢測數(shù)千甚至上萬個蛋白質在不同樣本之間的相對豐度或表達差異。這類數(shù)據(jù)已有的差異表達分析方法大多依賴于對并行或已有的技術重復數(shù)據(jù)進行前期比較來構建實驗的技術誤差模型，并以它為基礎檢驗每個蛋白質在被比較樣本之間表達差異的統(tǒng)計顯著性。該方法占用了有限的實驗通道，也難以保證誤差模型的精確適用性。
針對這一局限，在MAP模型中研究人員發(fā)展了一種新穎的分步回歸（step-by-step regression）分析流程，實現(xiàn)直接對被比較的兩個iTRAQ樣本構建技術誤差模型。在此類研究中，一個常用的經(jīng)驗假設是技術誤差對樣本間每個蛋白質iTRAQ信號log2比率（log2-ratio）的貢獻服從以0為中心的正態(tài)分布N（0, σ2）。其中，方差σ2依賴于該蛋白質的信號強度，并且常被用一個指數(shù)衰減函數(shù)來刻畫其依賴關系，即所要構建的全局誤差函數(shù)。MAP模型首先使用滑動窗口掃描兩個樣本的M-A圖，同時對窗口中0附近的log2比率進行線性建模，以其斜率的平方作為誤差函數(shù)的局域估計。然后，對所得局域估計進行第二輪指數(shù)擬合，獲得被比較樣本的全局誤差函數(shù)，并以它為參照計算每個蛋白質信號差異的顯著性P值（圖一）。
圖一：MAP模型的分步回歸分析流程。
(a)局域線性擬合；(b)全局指數(shù)擬合構建技術誤差模型；(c)計算每個蛋白質信號差異的顯著性P值。
同位素標記定量蛋白質組數(shù)據(jù)長期存在比率壓縮的難題。研究人員使用MAP模型分別比較分析了三個批次產(chǎn)生的小鼠胚胎干細胞分化前后蛋白質組數(shù)據(jù)，發(fā)現(xiàn)蛋白質iTRAQ信號log2比率在不同批次間關聯(lián)很低（圖二a），可能是因為技術誤差對其貢獻所服從的正態(tài)分布N（0, σ2）在批次間各不相同。根據(jù)MAP模型，研究人員提出使用每個批次的全局誤差函數(shù)對其中每個蛋白質iTRAQ信號的log2比率進行重標度（rescaling），使得在不同批次中技術誤差對其貢獻均服從標準正態(tài)分布N(0, 1)，從而發(fā)展了一個新的Z統(tǒng)計量。比較不同批次蛋白質Z統(tǒng)計量之間的關聯(lián)，可以發(fā)現(xiàn)它具有明顯更好的可重復性（圖二b）。
圖二：不同批次蛋白質iTRAQ信號的log2比率(a)和Z統(tǒng)計量(b)的皮爾森關聯(lián)系數(shù)。
為方便蛋白質組領域研究者使用MAP模型進行數(shù)據(jù)分析，研究人員搭建了一個網(wǎng)絡服務平臺（http://bioinfo.sibs.ac.cn/shaolab/MAP）。該平臺額外搭載了一個整合分析模塊，能夠通過整合多個批次生物重復比較結果的次優(yōu)P值或者平均Z統(tǒng)計量來最終鑒定差異表達蛋白質，并新發(fā)展了一種分析方法用于估測基于這些統(tǒng)計量所定義差異表達蛋白質的錯誤發(fā)現(xiàn)率（FDR）。此外，在用于雙樣本比較的MAP模型基礎上，研究人員還通過分別用樣本方差和卡方分布分位數(shù)取代原分步回歸建模流程中所使用的log2比率和標準正態(tài)分布分位數(shù)，進一步發(fā)展了適用于多樣本比較的拓展eMAP模型。
上述研究由中國科學院上海營養(yǎng)與健康所研究助理李木山和博士研究生涂世奇等在邵振研究員的指導下，與中國科學院植物生理生態(tài)研究所、復旦大學上海醫(yī)學院和美國西南醫(yī)學中心等多家單位的研究人員合作完成，得到了國家自然科學基金委、科技部和中國科學院等多項基金的資助。
官方微信
