本公開涉及生物信息領(lǐng)域，特別是涉及一種細(xì)胞通訊分析方法、裝置、電子設(shè)備和存儲介質(zhì)。

背景技術(shù)：

1、多細(xì)胞生物的細(xì)胞與細(xì)胞之間往往會通過細(xì)胞因子和膜蛋白等進(jìn)行通訊，從而調(diào)節(jié)生命活動。因此，受體與配體介導(dǎo)的細(xì)胞間通訊對細(xì)胞分化、組織和器官內(nèi)的細(xì)胞類型相互作用、與免疫反應(yīng)的研究上起到重要的作用。相關(guān)技術(shù)中，為了高效且準(zhǔn)確的對細(xì)胞中的細(xì)胞互作關(guān)系進(jìn)行分析，涌現(xiàn)了大量的細(xì)胞通訊分析應(yīng)用。細(xì)胞通訊分析應(yīng)用可以基于配受體數(shù)據(jù)集進(jìn)行細(xì)胞通訊分析，但是目前，細(xì)胞通訊分析結(jié)果可靠性的較低。

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

1、本公開實施例提供了一種細(xì)胞通訊分析方法、裝置、電子設(shè)備和存儲介質(zhì)，該方法可以提高細(xì)胞通訊分析結(jié)果的可靠性。

2、根據(jù)本公開的一方面，提供了一種細(xì)胞通訊分析方法，包括：

3、獲取目標(biāo)配受體數(shù)據(jù)集，所述目標(biāo)配受體數(shù)據(jù)集為對多個預(yù)設(shè)配受體數(shù)據(jù)集與自建配受體數(shù)據(jù)集進(jìn)行整合得到的配受體數(shù)據(jù)集，所述自建配受體數(shù)據(jù)集為基于多個基因庫中的潛在配受體基因構(gòu)建的配受體數(shù)據(jù)集；

4、基于所述目標(biāo)配受體數(shù)據(jù)集與待分析細(xì)胞群的空間轉(zhuǎn)錄數(shù)據(jù)及細(xì)胞注釋數(shù)據(jù)生成多個細(xì)胞通訊分析應(yīng)用對應(yīng)的細(xì)胞通訊分析數(shù)據(jù)，并將每一細(xì)胞通訊分析數(shù)據(jù)輸入至相應(yīng)的細(xì)胞通訊分析應(yīng)用中進(jìn)行細(xì)胞通訊分析，得到多個子分析結(jié)果；

5、對所述多個子分析結(jié)果進(jìn)行并集處理，得到融合分析結(jié)果；

6、基于多個所述子分析結(jié)果與所述融合分析結(jié)果，計算所述融合分析結(jié)果的第一可信度分?jǐn)?shù)與所述目標(biāo)配受體數(shù)據(jù)集中每一配受體對的第二可信度分?jǐn)?shù)；

7、當(dāng)所述第一可信度分?jǐn)?shù)大于預(yù)設(shè)閾值時，基于所述第二可信度分?jǐn)?shù)篩選出目標(biāo)配受體對，并根據(jù)所述目標(biāo)配受體對從所述融合分析結(jié)果中確定目標(biāo)細(xì)胞分析結(jié)果。

8、根據(jù)本公開的一方面，提供了一種細(xì)胞通訊分析裝置，包括：

9、第一獲取單元，用于獲取目標(biāo)配受體數(shù)據(jù)集，所述目標(biāo)配受體數(shù)據(jù)集為對多個預(yù)設(shè)配受體數(shù)據(jù)集與自建配受體數(shù)據(jù)集進(jìn)行整合得到的配受體數(shù)據(jù)集，所述自建配受體數(shù)據(jù)集為基于多個基因庫中的潛在配受體基因構(gòu)建的配受體數(shù)據(jù)集；

10、分析單元，用于基于所述目標(biāo)配受體數(shù)據(jù)集與待分析細(xì)胞群的空間轉(zhuǎn)錄數(shù)據(jù)及細(xì)胞注釋數(shù)據(jù)生成多個細(xì)胞通訊分析應(yīng)用對應(yīng)的細(xì)胞通訊分析數(shù)據(jù)，并將每一細(xì)胞通訊分析數(shù)據(jù)輸入至相應(yīng)的細(xì)胞通訊分析應(yīng)用中進(jìn)行細(xì)胞通訊分析，得到多個子分析結(jié)果；

11、合并處理單元，用于對所述多個子分析結(jié)果進(jìn)行并集處理，得到融合分析結(jié)果；

12、計算單元，用于基于多個所述子分析結(jié)果與所述融合分析結(jié)果，計算所述融合分析結(jié)果的第一可信度分?jǐn)?shù)與所述目標(biāo)配受體數(shù)據(jù)集中每一配受體對的第二可信度分?jǐn)?shù)；

13、確定單元，用于當(dāng)所述第一可信度分?jǐn)?shù)大于預(yù)設(shè)閾值時，基于所述第二可信度分?jǐn)?shù)篩選出目標(biāo)配受體對，并根據(jù)所述目標(biāo)配受體對從所述融合分析結(jié)果中確定目標(biāo)細(xì)胞分析結(jié)果。

14、可選地，每一所述配受體對與對應(yīng)的配受體可信度分?jǐn)?shù)在所述目標(biāo)配受體數(shù)據(jù)集中對應(yīng)存儲，所述配受體可信度分?jǐn)?shù)指示了所述配受體對在通用場景中的可信度；

15、所述計算單元具體用于：

16、基于多個所述子分析結(jié)果與所述融合分析結(jié)果，計算所述融合分析結(jié)果的第一可信度分?jǐn)?shù)與每一配受體對在所述融合分析結(jié)果中的第三可信度分?jǐn)?shù)；

17、基于每一所述配受體對的所述第三可信度分?jǐn)?shù)與所述配受體可信度分?jǐn)?shù)，計算每一所述配受體對的第二可信度分?jǐn)?shù)。

18、可選地，所述融合分析結(jié)果中包含所述待分析細(xì)胞群中每一細(xì)胞基于所述目標(biāo)配受體數(shù)據(jù)集中的每一所述配受體對進(jìn)行細(xì)胞通訊的互作結(jié)果；

19、所述計算單元具體用于：

20、針對所述融合分析結(jié)果中的每一所述互作結(jié)果，確定多個所述子分析結(jié)果中包含所述互作結(jié)果的所述子分析結(jié)果的第一數(shù)目；

21、基于所述第一數(shù)目與所述子分析結(jié)果的總數(shù)目，計算所述互作結(jié)果對應(yīng)的互作結(jié)果可信度分?jǐn)?shù)；

22、基于所述融合分析結(jié)果中的多個所述互作結(jié)果對應(yīng)的所述互作結(jié)果可信度分?jǐn)?shù)，計算所述融合分析結(jié)果的第一可信度分?jǐn)?shù)與每一配受體對在所述融合分析結(jié)果中的第三可信度分?jǐn)?shù)。

23、可選地，所述計算單元具體用于：

24、計算所述融合分析結(jié)果中多個所述互作結(jié)果對應(yīng)的所述互作結(jié)果可信度分?jǐn)?shù)的平均值，作為所述融合分析結(jié)果的第一可信度分?jǐn)?shù)；

25、針對每一配受體對，獲取多個所述細(xì)胞與所述配受體對的多個所述互作結(jié)果對應(yīng)的多個所述互作結(jié)果可信度分?jǐn)?shù)，并計算多個所述互作結(jié)果可信度分?jǐn)?shù)的平均值，作為所述配受體對在所述融合分析結(jié)果中的第三可信度分?jǐn)?shù)。

26、可選地，所述目標(biāo)配受體數(shù)據(jù)集通過以下方式構(gòu)建：

27、獲取多個所述預(yù)設(shè)配受體數(shù)據(jù)集與所述自建配受體數(shù)據(jù)集；

28、將所述預(yù)設(shè)配受體數(shù)據(jù)集中的配受體對與所述自建配受體數(shù)據(jù)集中的配受體對整合；

29、針對整合后的配受體數(shù)據(jù)集中的每一所述配受體對生成配受體注釋，將所述配受體對與所述配受體注釋對應(yīng)存儲，構(gòu)建所述目標(biāo)配受體數(shù)據(jù)集。

30、可選地，所述針對整合后的配受體數(shù)據(jù)集中的每一所述配受體對生成配受體注釋，將所述配受體對與所述配受體注釋對應(yīng)存儲，構(gòu)建所述目標(biāo)配受體數(shù)據(jù)集，包括：

31、獲取第一參考配受體數(shù)據(jù)集；

32、針對整合后的配受體數(shù)據(jù)集中的每一所述配受體對，基于所述第一參考配受體數(shù)據(jù)集獲取所述配受體對的配受體可信度分?jǐn)?shù)；

33、將所述配受體可信度分?jǐn)?shù)作為所述配受體對的配受體注釋，將所述配受體對與所述配受體注釋對應(yīng)存儲，構(gòu)建所述目標(biāo)配受體數(shù)據(jù)集。

34、可選地，所述自建配受體數(shù)據(jù)集通過以下方式構(gòu)建：

35、獲取多個所述基因庫，并從所述基因庫中獲取多個所述潛在配受體基因；

36、基于多個所述潛在配受體基因構(gòu)建潛在配受體對；

37、獲取第二參考配受體數(shù)據(jù)集，并基于所述第二參考配受體數(shù)據(jù)集從所述潛在配受體對中篩選得到目標(biāo)潛在配受體對；

38、將所述目標(biāo)潛在配受體對構(gòu)建為所述自建配受體數(shù)據(jù)集。

39、可選地，所述細(xì)胞通訊分析裝置還包括：

40、第二獲取單元，用于基于所述目標(biāo)配受體數(shù)據(jù)集與多個所述預(yù)設(shè)配受體數(shù)據(jù)集，獲取多個待評估配受體數(shù)據(jù)集；

41、篩選單元，用于從多個所述待評估配受體數(shù)據(jù)集中篩選出評估后配受體數(shù)據(jù)集；

42、所述分析單元具體用于：

43、基于所述評估后配受體數(shù)據(jù)集與待分析細(xì)胞群的空間轉(zhuǎn)錄數(shù)據(jù)及細(xì)胞注釋數(shù)據(jù)生成多個細(xì)胞通訊分析應(yīng)用對應(yīng)的細(xì)胞通訊分析數(shù)據(jù)，并將每一細(xì)胞通訊分析數(shù)據(jù)輸入至相應(yīng)的細(xì)胞通訊分析應(yīng)用中進(jìn)行細(xì)胞通訊分析，得到多個子分析結(jié)果；

44、所述計算單元具體用于：

45、基于多個所述子分析結(jié)果與所述融合分析結(jié)果，計算所述融合分析結(jié)果的第一可信度分?jǐn)?shù)與所述評估后配受體數(shù)據(jù)集中每一配受體對的第二可信度分?jǐn)?shù)。

46、可選地，所述篩選單元具體用于：

47、基于每一所述待評估配受體數(shù)據(jù)集中包含的多個配受體對，確定所述待評估配受體數(shù)據(jù)集的第四可信度分?jǐn)?shù)；

48、根據(jù)多個所述待評估配受體數(shù)據(jù)集的所述第四可信度分?jǐn)?shù)，篩選出評估后配受體數(shù)據(jù)集。

49、可選地，所述融合分析結(jié)果中包含所述待分析細(xì)胞群中的細(xì)胞與所述目標(biāo)配受體數(shù)據(jù)集中配受體對之間的第一互作矩陣、細(xì)胞類型組與所述目標(biāo)配受體數(shù)據(jù)集中所述配受體對之間的第二互作矩陣與不同所述細(xì)胞類型之間的第三互作矩陣；

50、所述計算單元具體用于：

51、基于多個所述子分析結(jié)果與所述融合分析結(jié)果中的所述第一互作矩陣，計算所述第一互作矩陣的第一子可信度分?jǐn)?shù)與所述目標(biāo)配受體數(shù)據(jù)集中每一配受體對的第二可信度分?jǐn)?shù)；基于多個所述子分析結(jié)果與所述融合分析結(jié)果中的所述第二互作矩陣，計算所述第二互作矩陣的第二子可信度分?jǐn)?shù)；基于多個所述子分析結(jié)果與所述融合分析結(jié)果中的所述第三互作矩陣，計算所述第三互作矩陣的第三子可信度分?jǐn)?shù)；

52、根據(jù)所述第一子可信度分?jǐn)?shù)、所述第二子可信度分?jǐn)?shù)與所述第三子可信度分?jǐn)?shù)，計算所述融合分析結(jié)果的第一可信度分?jǐn)?shù)。

53、可選地，所述確定單元具體用于：

54、基于所述第一互作矩陣中獲取所述目標(biāo)配受體對的配受體互作強(qiáng)度與所述待分析細(xì)胞群中每一所述細(xì)胞的細(xì)胞互作強(qiáng)度；

55、基于所述第二互作矩陣，生成所述細(xì)胞類型組與所述目標(biāo)配受體對的所述第一互作關(guān)系圖；

56、基于所述第三互作矩陣，生成不同所述細(xì)胞類型之間的所述第二互作關(guān)系圖；

57、將所述細(xì)胞的細(xì)胞互作強(qiáng)度、所述目標(biāo)配受體對的所述配受體互作強(qiáng)度、所述第一互作關(guān)系圖與所述第二互作關(guān)系圖作為目標(biāo)細(xì)胞分析結(jié)果。

58、根據(jù)本公開的一方面，提供了一種計算機(jī)設(shè)備，包括存儲器和處理器，所述存儲器存儲有計算機(jī)程序，所述處理器執(zhí)行所述計算機(jī)程序時實現(xiàn)如上所述的細(xì)胞通訊分析方法。

59、根據(jù)本公開的一方面，提供了一種計算機(jī)可讀存儲介質(zhì)，所述存儲介質(zhì)存儲有計算機(jī)程序，所述計算機(jī)程序被處理器執(zhí)行時實現(xiàn)如上所述的細(xì)胞通訊分析方法。

60、根據(jù)本公開的一方面，提供了一種計算機(jī)程序產(chǎn)品，該計算機(jī)程序產(chǎn)品包括計算機(jī)程序，所述計算機(jī)程序被計算機(jī)設(shè)備的處理器讀取并執(zhí)行，使得該計算機(jī)設(shè)備執(zhí)行如上所述的細(xì)胞通訊分析方法。

61、本公開實施例中的細(xì)胞通訊分析方法，獲取目標(biāo)配受體數(shù)據(jù)集，目標(biāo)配受體數(shù)據(jù)集為對多個預(yù)設(shè)配受體數(shù)據(jù)集與自建配受體數(shù)據(jù)集進(jìn)行整合得到的配受體數(shù)據(jù)集，自建配受體數(shù)據(jù)集為基于多個基因庫中的潛在配受體基因構(gòu)建的配受體數(shù)據(jù)集；基于目標(biāo)配受體數(shù)據(jù)集與待分析細(xì)胞群的空間轉(zhuǎn)錄數(shù)據(jù)及細(xì)胞注釋數(shù)據(jù)生成多個細(xì)胞通訊分析應(yīng)用對應(yīng)的細(xì)胞通訊分析數(shù)據(jù)，并將每一細(xì)胞通訊分析數(shù)據(jù)輸入至相應(yīng)的細(xì)胞通訊分析應(yīng)用中進(jìn)行細(xì)胞通訊分析，得到多個子分析結(jié)果；對多個子分析結(jié)果進(jìn)行并集處理，得到融合分析結(jié)果；基于多個子分析結(jié)果與融合分析結(jié)果，計算融合分析結(jié)果的第一可信度分?jǐn)?shù)與目標(biāo)配受體數(shù)據(jù)集中每一配受體對的第二可信度分?jǐn)?shù)；當(dāng)?shù)谝豢尚哦确謹(jǐn)?shù)大于預(yù)設(shè)閾值時，基于第二可信度分?jǐn)?shù)篩選出目標(biāo)配受體對，并根據(jù)目標(biāo)配受體對從融合分析結(jié)果中確定目標(biāo)細(xì)胞分析結(jié)果。

62、以此，本公開實施例的細(xì)胞通訊分析方法中，目標(biāo)配受體數(shù)據(jù)集中既包括預(yù)設(shè)配受體數(shù)據(jù)集中的配受體對，也包括從基因庫中獲得的潛在配受體對。因此，目標(biāo)配受體數(shù)據(jù)集中的配受體對較為全面?；谶@個目標(biāo)配受體數(shù)據(jù)集進(jìn)行細(xì)胞通訊分析可以提高分析結(jié)果的可靠性。多個細(xì)胞通訊分析應(yīng)用根據(jù)對應(yīng)的細(xì)胞通訊分析數(shù)據(jù)進(jìn)行細(xì)胞通訊分析，得到多個子分析結(jié)果。將多個子分析結(jié)果取并集進(jìn)行整合后得到融合分析結(jié)果。利用多個子分析結(jié)果與融合分析結(jié)果來計算融合分析結(jié)果的第一可信度分?jǐn)?shù)，因此，第一可信度分?jǐn)?shù)越高，就表示融合分析結(jié)果得到了越多細(xì)胞通訊分析應(yīng)用的支持，融合分析結(jié)果的可靠性越高。因此，通過設(shè)定預(yù)設(shè)閾值來判定融合分析結(jié)果的可靠性。細(xì)胞通訊分析結(jié)果中對于目標(biāo)配受體集中配受體對的分析是十分重要的。因此，利用多個子分析結(jié)果與融合分析結(jié)果來計算每一配受體對的第二可信度分?jǐn)?shù)。第二可信度分?jǐn)?shù)越高也就表示了配受體對得到越多細(xì)胞通訊分析應(yīng)用的支持，配受體對的可靠性越高。因此，利用第二可信度分?jǐn)?shù)篩選出可靠的配受體對作為目標(biāo)配受體對，并根據(jù)目標(biāo)配受體從融合分析結(jié)果中確定目標(biāo)細(xì)胞分析結(jié)果。基于此，目標(biāo)細(xì)胞分析結(jié)果是基于對融合分析結(jié)果的可靠性與配受體對的可靠性的雙重篩選得到的。因此，利用本公開實施例的細(xì)胞通訊分析方法得到的目標(biāo)細(xì)胞分析結(jié)果的可靠性較高。

63、本公開的其他特征和優(yōu)點將在隨后的說明書中闡述，并且，部分地從說明書中變得顯而易見，或者通過實施本公開而了解。本公開的目的和其他優(yōu)點可通過在說明書、權(quán)利要求書以及附圖中所特別指出的結(jié)構(gòu)來實現(xiàn)和獲得。