本發(fā)明涉及一種微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法及應(yīng)用，屬于檢驗醫(yī)學(xué)。

背景技術(shù)：

1、基質(zhì)輔助激光解析電離飛行時間質(zhì)譜(matrix-assisted?laser?desorptionionization-time?of?flight?mas?spectrometry，maldi-tof?ms)是近年來發(fā)展的全新用于微生物鑒定和分型的生物質(zhì)譜技術(shù)。其通過離子的質(zhì)荷比(m/z)與離子的飛行時間成正比的原理來分析離子，測定樣品分子量并分析混合生物大分子的軟電離技術(shù)。其通過采集2000～20000m/z區(qū)間的微生物蛋白分子指紋圖譜并與數(shù)據(jù)庫中的蛋白圖譜進行識別比對，從而判斷細菌的種屬。與基于免疫和核酸的檢測相比，maldi-tof-ms具有靈敏度高、準確度高及分辨率高等特點，而且使用的試劑成本相對較低，不需要專門培訓(xùn)的人員操作，是一種高通量篩選鑒定方法，一次運行能識別大量菌株并在數(shù)分鐘內(nèi)得出結(jié)論。maldi-tof-ms在微生物檢測鑒定中的應(yīng)用越來越深入，質(zhì)譜檢測相關(guān)的技術(shù)指南和專家共識也已經(jīng)出版發(fā)行，質(zhì)譜檢測技術(shù)日趨成熟，在對微生物進行菌種鑒定的過程中獲得了海量的核糖體蛋白質(zhì)組數(shù)據(jù)，但是如何對獲得的海量數(shù)據(jù)進行有效的數(shù)據(jù)挖掘是一個技術(shù)難題。

2、在臨床抗感染治療中，細菌耐藥是一個世界性的難題，如何通過菌種鑒定的途徑同時檢測細菌是否耐藥，在公布日為2023年02月03日，公布號為cn115678958a的中國發(fā)明專利中通過4個質(zhì)譜峰進行區(qū)分cskp(碳青霉烯類敏感的肺炎克雷伯菌)和crkp(碳青霉烯類耐藥的肺炎克雷伯菌)；也有研究通過40個特征峰構(gòu)建kpc-crkp自建質(zhì)譜庫的方法解決該問題。自2020年以來，不斷有文獻報道細菌如：肺炎克雷伯菌中已有耐藥表型和毒力表型的融合，那么能否在利用maldi-tof-ms進行菌種鑒定和耐藥檢測的同時，有效挖掘核糖體蛋白質(zhì)組數(shù)據(jù)以實現(xiàn)包括細菌毒力狀況的檢測，逐漸成為目前亟待解決的問題。

技術(shù)實現(xiàn)思路

1、本發(fā)明的第一個目的是提供一種微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法，以解決現(xiàn)有技術(shù)中無法有效挖掘maldi-tof-ms質(zhì)譜獲得的微生物核糖體蛋白質(zhì)組數(shù)據(jù)問題。

2、本發(fā)明的第二個目的是提供一種微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法在鑒定耐碳青霉烯類肺炎克雷伯菌毒力水平中的應(yīng)用，為現(xiàn)有技術(shù)提供一種快速、簡便的鑒定耐碳青霉烯類肺炎克雷伯菌毒力水平的方法。

3、為了實現(xiàn)上述目的，本發(fā)明中一種微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法所采用的技術(shù)方案是：

4、一種微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法，包括以下步驟：

5、s1、樣本前處理：將待分析的微生物樣本進行培養(yǎng)并根據(jù)預(yù)設(shè)表型分組，然后進行基質(zhì)輔助激光解析電離飛行時間質(zhì)譜檢測獲得核糖體蛋白質(zhì)譜譜圖數(shù)據(jù)集；

6、s2、質(zhì)譜數(shù)據(jù)集預(yù)處理：對s1獲得的核糖體蛋白質(zhì)譜譜圖數(shù)據(jù)集進行預(yù)處理，獲得微生物樣本的共同質(zhì)譜譜圖數(shù)據(jù)集；

7、s3、質(zhì)譜數(shù)據(jù)集基礎(chǔ)分析：對s1和s2獲得的質(zhì)譜譜圖進行質(zhì)量控制分析、種群代表性質(zhì)譜峰生成及主成分分析，獲取微生物樣本的種群核糖體蛋白表達譜；

8、s4、質(zhì)譜數(shù)據(jù)集差異性分析：對s2中所述微生物樣本的共同質(zhì)譜譜圖數(shù)據(jù)集進行矩陣分析及數(shù)據(jù)處理篩選出差異性質(zhì)譜峰，結(jié)合s3中所述微生物樣本的種群核糖體蛋白表達譜進行層次聚類分析輸出可視化熱圖，進一步做簇類分析；

9、s5、質(zhì)譜數(shù)據(jù)集分子標記篩選：對s4獲得差異性質(zhì)譜峰篩選分析得到候選分子標記，對不同候選分子標記進行組合分析及臨床診斷效能評價確定微生物預(yù)設(shè)表型分組的連鎖分子標記。

10、上述方案的有益效果在于：本發(fā)明的微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法為開拓性發(fā)明。本發(fā)明以耐碳青霉烯類肺炎克雷伯菌為具體檢測對象，以其毒力為具體的分析表型，開發(fā)出一種微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法。本發(fā)明將耐碳青霉烯類肺炎克雷伯菌依據(jù)毒力水平分為無毒力組、中間毒力組和高毒力組，將微生物質(zhì)譜檢測中的蛋白質(zhì)數(shù)據(jù)集經(jīng)質(zhì)量控制分析及數(shù)據(jù)處理后共獲取389個共有質(zhì)譜峰，進一步通過矩陣分析和聚類分析發(fā)現(xiàn)130個差異性蛋白峰，繼續(xù)進行數(shù)據(jù)分析篩選出3個聯(lián)合標記質(zhì)譜峰用于區(qū)分不同毒力水平的細菌，形成了一種大樣本量、高重現(xiàn)性、完整可靠的核糖體蛋白質(zhì)組數(shù)據(jù)分析方法。

11、作為進一步地改進，s1中所述預(yù)設(shè)表型包括微生物的毒力或耐藥性。

12、上述方案的有益效果在于：致病微生物的毒力和耐藥性為是臨床醫(yī)學(xué)研究和病原微生物研究中，對于致病微生物重點關(guān)注的領(lǐng)域，而且研究的比較深入，對于相關(guān)毒力和耐藥性有明確且清晰的劃分。

13、優(yōu)選地，當(dāng)預(yù)設(shè)表型為微生物毒力或耐藥性時，可以根據(jù)現(xiàn)有技術(shù)中公開的相應(yīng)毒力或耐藥性檢測的方法，如：熒光定量pcr、藥物敏感性鑒定，實驗檢測后進行分組。

14、具體地，在進行實際分析時，所述預(yù)設(shè)表型還可以為細菌的st分型、表面抗原分型(血清分型)等。

15、作為進一步地改進，s1中所述數(shù)據(jù)集的子集要求鑒定百分比≥75％、質(zhì)譜峰數(shù)量為100-200。

16、作為進一步地改進，s2中所述微生物樣本的共同質(zhì)譜譜圖數(shù)據(jù)集來自于不同的生物學(xué)重復(fù)和技術(shù)重復(fù)。

17、作為進一步地改進，s3中所述質(zhì)量控制分析包含原始質(zhì)譜譜圖數(shù)據(jù)集和共同質(zhì)譜譜圖數(shù)據(jù)集。

18、作為進一步地改進，s4中所述差異性質(zhì)譜峰的篩選如下：①排除所有菌株中出現(xiàn)頻率100％的質(zhì)譜峰；②排除所有菌株內(nèi)出現(xiàn)頻率＜20％的質(zhì)譜峰。

19、作為進一步地改進，s4中所述簇類分析從預(yù)設(shè)表型分組和質(zhì)譜峰分布上雙向綜合分析。

20、作為進一步地改進，s5中所述候選分子標記的篩選如下：①待分析質(zhì)譜峰在預(yù)設(shè)表型分組的目標組內(nèi)pop＜20％；②待分析質(zhì)譜峰在預(yù)設(shè)表型分組其他組內(nèi)表達相對較高。

21、作為進一步地改進，s5中所述組合分析使用排列組合的方法，公式如下：

22、

23、式中，m表示候選蛋白聯(lián)合標記組合種類數(shù)量，n表示可供篩選的候選分子標記數(shù)量。

24、作為進一步地改進，s5中所述組合分析后使用分子標記的臨床診斷效能評價指標進行評估。

25、具體地，包括靈敏性、特異度、準確度、陽性預(yù)測值和陰性預(yù)測值的評估。

26、作為進一步地改進，s1中待分析的微生物樣本培養(yǎng)至至少f2代。

27、上述方案的有益效果在于：在待分析的微生物樣本培養(yǎng)至f2代后，取同一樣本的不同后代(如：同一種菌株的不同菌落)進行分組檢測，獲得質(zhì)譜分析數(shù)據(jù)，能夠減小實驗分析過程中的誤差。

28、為了實現(xiàn)上述目的，本發(fā)明中一種微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法在鑒定耐碳青霉烯類肺炎克雷伯菌毒力水平中的應(yīng)用所采用的技術(shù)方案是：

29、一種微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法在鑒定耐碳青霉烯類肺炎克雷伯菌毒力水平中的應(yīng)用。

30、上述方案的有益效果在于：本發(fā)明微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法在鑒定耐碳青霉烯類肺炎克雷伯菌毒力水平中的應(yīng)用為開拓性發(fā)明。通過使用本發(fā)明的深入分析方法篩選出了耐碳青霉烯類肺炎克雷伯菌高毒力菌株的聯(lián)合標記質(zhì)譜峰、中間毒力菌株的聯(lián)合標記質(zhì)譜峰和無毒力菌株的聯(lián)合標記質(zhì)譜峰。經(jīng)過分析驗證，由上述三個聯(lián)合標記質(zhì)譜峰組成的蛋白分子標記可以準確、快速地鑒定出待測耐碳青霉烯類肺炎克雷伯菌的毒力水平。這樣既能充分利用在菌種鑒定中獲得的大量核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)，又能直接將分析的結(jié)果用于臨床試驗檢測數(shù)據(jù)的延伸化分析，這將對臨床抗感染的精準治療帶來極大提升，并提早了精準治療的時間窗和降低了患者的治療成本，甚至挽救患者生命。而且相較于蛋白和核酸的檢測，利用三個聯(lián)合標記質(zhì)譜峰組成的蛋白分子標記進行耐碳青霉烯類肺炎克雷伯菌毒力的鑒定，無需額外實驗操作及檢驗成本，無檢測時間上的延遲，適合大規(guī)模、高通量檢測鑒定。

31、作為進一步地改進，對待測耐碳青霉烯類肺炎克雷伯菌的核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集進行深度分析，根據(jù)不同毒力的聯(lián)合標記質(zhì)譜峰對待測耐碳青霉烯類肺炎克雷伯菌的毒力水平進行判定。

32、作為進一步地改進，高毒力耐碳青霉烯類肺炎克雷伯菌的聯(lián)合標記質(zhì)譜峰(m/z)為2368.64-2525.54-3130.11-3611.80、中間毒力耐碳青霉烯類肺炎克雷伯菌的聯(lián)合標記質(zhì)譜峰(m/z)為2761.26-11178.83-11212.18-11873.99、無毒力耐碳青霉烯類肺炎克雷伯菌的聯(lián)合標記質(zhì)譜峰(m/z)為2502.75-2580.52-4401.36-11178.83。

33、作為進一步地改進，所述高、中間和無毒力耐碳青霉烯類肺炎克雷伯菌聯(lián)合標記質(zhì)譜峰中每個質(zhì)譜峰檢測數(shù)據(jù)的容忍度為800ppm。

34、與現(xiàn)有技術(shù)相比，本發(fā)明的有益效果在于：

35、1、本發(fā)明微生物核糖體蛋白質(zhì)組學(xué)數(shù)據(jù)集的深度分析方法中篩選的聯(lián)合質(zhì)譜峰蛋白分子標記通過在微生物質(zhì)譜檢測的基礎(chǔ)上分析該菌種鑒定的質(zhì)譜峰數(shù)據(jù)來判定其毒力狀況，是臨床實驗檢測數(shù)據(jù)的延伸化分析。

36、2、本發(fā)明中的三組聯(lián)合質(zhì)譜峰蛋白標記可直接應(yīng)用于臨床診療，無需額外實驗操作及檢驗成本，無檢測時間上的延遲，直接在臨床常規(guī)微生物質(zhì)譜菌種鑒定的同時比對其檢測的質(zhì)譜峰數(shù)據(jù)：2368.64-2525.54-3130.11-3611.80、2761.26-11178.83-11212.18-11873.99、2502.75-2580.52-4401.36-11178.8，即可判定耐碳青霉烯類肺炎克雷伯菌的毒力水平。

37、3、本發(fā)明的目的在于提供一種大樣本量、高重現(xiàn)性、完整可靠且操作簡便的微生物核糖體蛋白質(zhì)組數(shù)據(jù)集分析方法。與現(xiàn)有的轉(zhuǎn)錄組學(xué)、代謝組學(xué)及蛋白質(zhì)組學(xué)技術(shù)相比較，本發(fā)明中用于研究的每株菌株都有多個生物學(xué)試驗，將多個質(zhì)譜譜圖合并后產(chǎn)生的代表性質(zhì)譜譜圖作為該菌株的質(zhì)譜數(shù)據(jù)集，從而保證了后續(xù)數(shù)據(jù)分析結(jié)果的高重現(xiàn)性和高穩(wěn)定性；本發(fā)明可進行大樣本量的數(shù)據(jù)分析處理，在發(fā)現(xiàn)差異表達分子的基礎(chǔ)上直接篩選出區(qū)分相關(guān)表型的分子標記。