作者,Evil Genius
最近收到了很多粉丝的提问,其中对于空间通讯信号的用法需要在加强一点。
其实单独运用COMMOT做通讯信号流的做法并不完美,而实应该和空间轨迹向量场结合起来才更具有意义。
所以这里我们需要做到的是,
1、细胞的空间分布情况
2、细胞的空间分布趋势
3、细胞之间的共定位情况
4、细胞的信号流向
前三点其实分享了很多了,关于空间轨迹向量场也都分享过,这个地方我提醒大家一点,空间轨迹向量场表征的是细胞或者通路的聚集方向,信号流是通讯的方向,两者之间的结合才更具有意义
COMMOT软件特点
a package that infers CCC by simultaneously considering numerous ligand–receptor pairs for either spatial transcriptomics data or spatially annotated scRNA-seq data equipped with spatial distances between cells estimated from paired spatial imaging data; summarizes and compares directions of spatial signaling; identifies downstream effects of CCC on gene expressions using ensemble of trees models; and provides visualization utilities for the various analyses.
软件总览----COMMOT
配体和受体通常在有限的空间范围内与多种复合物相互作用。考虑到这一点,作者提出了具有三个重要特征的collective optimal transport:首先,the use of non-probability mass distributions to control the marginals of the transport plan to maintain comparability between species(需要一点数学背景知识);其次,对CCC实施空间距离约束,以避免连接空间上相距较远的细胞;最后,将多种配体分布结合到多中受体分布以解释多种相互作用。
这里我以一个空间肝癌样本为例,分享一下COMMOT的用法。
第一步看形态学划分,这个看HE图片应该都能看出来。
数据的简单处理
import os
import gc
import ot
import pickle
import anndata
import scanpy as sc
import pandas as pd
import numpy as np
from scipy import sparse
from scipy.stats import spearmanr, pearsonr
from scipy.spatial import distance_matrix
import matplot