一种识别空间转录组空间区域和细胞类型的特征学习方法

1.本发明涉及生物信息学领域，更具体地，涉及一种识别空间转录组空间区域和细胞类型的特征学习方法。


背景技术：

2.空间转录组测序是一种新兴的测序技术，它能够以多细胞、单细胞、甚至亚细胞分辨率提供带有空间坐标位置的转录组信息，为细胞功能、表型和组织微环境中位置的关系提供了重要信息，而解析这些问题的第一步就是正确划分空间区域。
3.划分空间区域的关键在于如何学习单个细胞的特征，以及如何融合空间位置信息。现有的多个方法可以进行空间区域的划分。其中bayesspace(zhao,e.et al.“spatial transcriptomics at subspot resolution with bayesspace.(使用bayesspace的亚斑点分辨率空间转录组学)”nat biotechnol 39,1375-1384,2021)使用带有马尔可夫随机场的贝叶斯模型来实现空间聚类；spagcn(hu,j.et al.“spagcn:integrating gene expression,spatial location and histology to identify spatial domains and spatially variable genes by graph convolutional network.(spagcn：整合基因表达、空间位置和组织学，通过图卷积网络识别空间域和空间可变基因)”nat methods 18,1342-1351,2021)使用图卷积网络整合基因表达，空间位置以及组织图像来识别空间域；stagate(dong,k.&zhang,s，“deciphering spatial domains from spatially resolved transcriptomics with an adaptive graph attention auto-encoder.(使用自适应图注意力自动编码器从空间分辨的转录组学中破译空间域)”nat commun 13,1739,2022)使用自适应的图注意力自编码器学习空间区域。
4.现有方法的主要问题是：(1)对空间位置的利用不够，导致区域划分效果有限。现有方法仅基于每个细胞与相邻细胞的表达相似性聚合邻居信息，这是对表达数据的过度利用，从而造成表达相似细胞的过度平滑。(2)没有数据的测序噪声进行显式的建模，这使得这些方法难以处理稀疏程度高的空间转录组数据。但随着测序技术的发展，数据的分辨率逐渐升高，稀疏程度也随之提高，建模数据的测序噪声也变得尤为重要。(3)对于兼具单细胞分辨率的空间转录组数据，不能同时识别空间区域和细胞类型，使得数据分析流程复杂化。


技术实现要素：

5.本发明提供一种识别空间转录组空间区域和细胞类型的特征学习方法，同时识别空间转录组数据的空间区域和细胞类型，还保留了明确的生物学信号，可以实现良好的低维特征可视化和正确的轨迹推断。
6.为解决上述技术问题，本发明的技术方案如下：
7.一种识别空间转录组空间区域和细胞类型的特征学习方法，包括以下步骤：
8.s1：通过空间转录组测序方法获取定量的基因表达矩阵和空间坐标位置矩阵；
9.s2：对所述基因表达矩阵进行标准化，得到标准化后的基因表达矩阵；
10.s3：利用结合零膨胀负二项分布(zero-inflated negative binomial distribution,简称zinb)的分布假设的自编码器(autoencoder，简称ae)模块进行降噪低维特征的学习，所述自编码器模块的输入为标准化后的基因表达矩阵，所述自编码器模块的输出为细胞类型相关特征(cell type-related embedding，简称ce)，通过对细胞类型相关特征进行聚类得到每个细胞的细胞类型标签；
11.s4：将细胞类型相关特征和空间坐标位置矩阵作为图注意力网络(graph attention network，简称gat)模块的输入，图注意力网络模块输出空间区域相关特征(spatial domain-related embedding，简称se)，通过对空间区域相关特征进行聚类得到每个细胞的空间区域标签。
12.优选地，步骤s1中基因表达矩阵和空间坐标位置矩阵，具体为：
13.基因表达矩阵x中的每行代表一个测序基本单元，包含单个细胞或多个细胞，一共有n个细胞，基因表达矩阵x中的每列代表一个基因，一共有g个基因；
14.空间坐标位置矩阵y的每行代表一个细胞，空间坐标位置矩阵y的列是该细胞的二维或三维位置信息。
15.优选地，步骤s2中对所述基因表达矩阵进行标准化，具体为：
[0016][0017]
式中，xn表示细胞n的基因表达向量，为所述基因表达矩阵的第n行，sn表示细胞n的文库大小，sf表示标准化后所有细胞的文库大小(library size)，表示标准化后的细胞n的基因表达向量。
[0018]
优选地，所述基因表达矩阵中的每个元素都遵循零膨胀负二项分布，其参数为(π
ng
，r
ng
，pg)，其中π
ng
是观察到真实基因表达值为0的概率，(r
ng
，pg)是负二项分布的标准参数，x
ng
是在达到r
ng
次失败时的成功次数，pg是每个伯努利试验的失败概率。
[0019]
优选地，x
ng
的似然函数为：
[0020][0021]
其中，δ0()表示狄拉克函数，γ()表示伽马函数。
[0022]
优选地，步骤s3中自编码器模块包括两层编码器和两层解码器，表示为：
[0023][0024]z′
＝f2(z)
[0025]r′
＝f3(z
′
)
[0026]
π＝f4(z
′
)
[0027]
式中，f1是编码器，编码器包含两层，分别将输入数据从g维降维至m
′
维、从m
′
维降维至m维，表示标准化后的基因表达矩阵，z为降维后的细胞类型相关特征矩阵；f2、f3、f4分别为解码器的三个输出，解码器包括两层，分别将输入数据从m维升维至m
′
维、m
′
维升维至g维，z
′
为解码器的第一个输出，f3、f4分别作用于z
′
以学习参数r
′
和参数π，参数r和参数π分别为n
×
g维的矩阵，分别由元素r
ng
和π
ng
组成，r
′
由元素r
′
ng
组成，r
ng
＝s
′
nr′
ng
，其中s
′n是细胞n的缩放系数，s
′n＝exp(log
10
sn)。
[0028]
优选地，由元素lg组成的参数向量l是通过指数变换从可学习向量获得的，具体为：
[0029][0030]
式中，lg是负二项分布的logit。
[0031]
优选地，所述自编码器模块的损失函数定义为负二项分布的负对数似然函数(negative log-likelihood function，简称nll)：
[0032]
loss1＝nll
zinb
(x；π，r，p)
[0033]
式中，nll
zinb
()表示负二项分布的负对数似然函数。
[0034]
优选地，步骤s4具体为：
[0035]
首先根据各个细胞的位置构建对称的邻接矩阵a，节点的集合写成v＝{v1，v2，...，vn}，那么邻接矩阵a的每个元素a
ij
可以表示为：
[0036][0037]
式中，是节点vj的近邻节点集合，通过knn或距离截断来计算；
[0038]
构建相似性矩阵∑，令细胞之间的空间区域相关特征相关性随着距离的增加而减少，并趋于指数衰减，所以理想情况下∑是基于高斯核的n
×
n维空间相关矩阵：
[0039][0040]
式中，γ表示带宽参数；
[0041]
所述图注意力网络模块包括两个图注意层，令图注意层的输入是n
×
n维的特征矩阵为h＝(h1，h2，...，hn)，输出为n
×n′
维的矩阵h
′
＝(h
′1，h
′2，...，h
′n)，则：
[0042][0043]
式中，w是n
′×
n维的权重矩阵，是vj邻点的集合，α
ij
是使用softmax函数的归一化注意力系数矩阵：
[0044][0045]eij
＝a
t
(whi||whj)
[0046]
其中a是可学习向量，||是连接操作，在图注意层中使用指数线性单元作为激活函数σ；
[0047]
空间区域相关特征u经过学习，可以表示为：
[0048]
u＝gat2(gat1(z))
[0049]
式中，gat1()、gat2()分别表示第一层图注意层和第二层图注意层。
[0050]
优选地，所述图注意力网络模块的目标函数为：
[0051]
loss2＝λ
spatial
*l
spatial
+λ
rec
*l
rec
[0052]
式中，l
spatial
＝mse(uu
t
，∑)，l
rec
＝mse(u，z)，λ
spatial
和λ
rec
为权重参数。
[0053]
与现有技术相比，本发明技术方案的有益效果是：
[0054]
(1)本发明融合了统计分布假设和深度学习算法，使得模型兼具统计建模的表征能力和深度学习的学习效率，显著提高了空间转录组数据的区域划分能力和特征学习效果。
[0055]
(2)本发明的ae模块的建模学习到了精细可靠的细胞类型相关特征，因此对于具有单细胞分辨率的空间转录组数据，可以同时输出细胞类型和空间区域的标签，简化了空间转录组数据的分析步骤。
[0056]
(3)本发明的gat模块充分利用空间位置信息，利用全局的位置信息来防止局部过拟合的出现。同时，可以自适应地更新低维表征，保留了数据的生物学意义，学习到可解读的特征空间和可靠的轨迹推断结果。
附图说明
[0057]
图1为本发明的方法流程示意图。
[0058]
图2为本发明的方法框架示意图。
[0059]
图3为实施例提供的利用lisi指标评估多种数据的区域划分结果示意图。
[0060]
图4为实施例提供的利用ari、nmi、acc评估几种算法在皮层数据中的区域划分的准确性结果示意图。
具体实施方式
[0061]
附图仅用于示例性说明，不能理解为对本专利的限制；
[0062]
为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；
[0063]
对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。
[0064]
下面结合附图和实施例对本发明的技术方案做进一步的说明。
[0065]
实施例1
[0066]
一种识别空间转录组空间区域和细胞类型的特征学习方法，如图1所示，包括以下步骤：
[0067]
s1：通过空间转录组测序方法获取定量的基因表达矩阵和空间坐标位置矩阵；
[0068]
s2：对所述基因表达矩阵进行标准化，得到标准化后的基因表达矩阵；
[0069]
s3：利用结合零膨胀负二项分布的分布假设的自编码器模块进行降噪低维特征的学习，所述自编码器模块的输入为标准化后的基因表达矩阵，所述自编码器模块的输出为细胞类型相关特征，通过对细胞类型相关特征进行聚类得到每个细胞的细胞类型标签；
[0070]
s4：将细胞类型相关特征和空间坐标位置矩阵作为图注意力网络模块的输入，图注意力网络模块输出空间区域相关特征，通过对空间区域相关特征进行聚类得到每个细胞的空间区域标签。
[0071]
实施例2
[0072]
本实施例在实施例1的基础上，继续公开以下内容：
[0073]
步骤s1中基因表达矩阵和空间坐标位置矩阵，具体为：
[0074]
各种空间转录组测序方法都可以得到定量的基因表达矩阵和空间坐标位置矩阵，其中，基因表达矩阵x中的每行代表一个测序基本单元，包含单个细胞或多个细胞，一共有n个细胞，基因表达矩阵x中的每列代表一个基因，一共有g个基因；
[0075]
空间坐标位置矩阵y的每行代表一个细胞，空间坐标位置矩阵y的列是该细胞的二维或三维位置信息。
[0076]
步骤s2中对所述基因表达矩阵进行标准化，具体为：
[0077][0078]
式中，xn表示细胞n的基因表达向量，为所述基因表达矩阵的第n行，sn表示细胞n的文库大小，sf表示标准化后所有细胞的文库大小，表示标准化后的细胞n的基因表达向量，对每个细胞的表达式均进行上式所述的标准化，从而得到标准化后的表达矩阵
[0079]
模型的整个框架如图2所示，结合了零膨胀负二项分布和深度学习方法进行特征学习，模型包含自编码器和图注意力网络两个模块。其中，ae模块结合了zinb的分布假设，使用自编码器框架进行降噪低维特征的学习。此模块的输入为标准化后的基因表达矩阵输出为细胞类型相关特征ce。ce可以通过聚类进而得到每个细胞的细胞类型标签。gat模块将ce矩阵和坐标位置y作为输入，学习空间区域相关特征se矩阵。gat模块首先根据坐标位置构建邻接矩阵和相似矩阵，然后根据邻接矩阵整合邻居细胞的特征从而更新目标细胞的se，同时，为了充分利用空间位置，在训练过程中，约束se的相似性随着细胞物理距离的增加而减小。在学习得到se矩阵后，可以进一步通过聚类得到每个细胞的空间区域标签。
[0080]
所述基因表达矩阵中的每个元素都遵循零膨胀负二项分布，其参数为(π
ng
，r
ng
，pg)，其中π
ng
是观察到真实基因表达值为0的概率，(r
ng
，pg)是负二项分布的标准参数，x
ng
是在达到r
ng
次失败时的成功次数，pg是每个伯努利试验的失败概率。
[0081]
x
ng
的似然函数为：
[0082][0083]
其中，δ0()表示狄拉克函数，f()表示伽马函数。
[0084]
使用ae模块来同时学习z和分布参数，zinb参数是n
×
g维的参数矩阵π、r和g维的参数向量p，分别由元素π
ng
，r
ng
和pg组成。
[0085]
步骤s3中自编码器模块包括两层编码器和两层解码器，表示为：
[0086][0087]z′
＝f2(z)
[0088]r′
＝f3(z
′
)
[0089]
π＝f4(z
′
)
[0090]
式中，f1是编码器，编码器包含两层，分别将输入数据从g维降维至m
′
维、从m
′
维降维至m维，表示标准化后的基因表达矩阵，z为降维后的细胞类型相关特征矩阵(n
×
m维)；
f2、f3、f4分别为解码器的三个输出，解码器包括两层，分别将输入数据从m维升维至m
′
维、m
′
维升维至g维，z
′
为解码器的第一个输出，f3、f4分别作用于z
′
以学习参数r
′
和参数π，参数r和参数π分别为n
×
g维的矩阵，分别由元素r
ng
和π
ng
组成，r
′
由元素r
′
ng
组成，r
ng
＝s
′
nr′
ng
，其中s
′n是细胞n的缩放系数，在此，假设基因表达水平的对数与测序深度的基数10的对数之间存在线性关系，s
′n＝exp(log
10
sn)，所以重建的表达是
[0091]
由元素lg组成的参数向量l是通过指数变换从可学习向量获得的，具体为：
[0092][0093]
式中，lg是负二项分布的logit。
[0094]
所述自编码器模块的损失函数定义为负二项分布的负对数似然函数：
[0095]
loss1＝nll
zinb
(x；π，r，p)
[0096]
式中，nll
zinb
()表示负二顶分布的负对数似然函数。
[0097]
gat模块的输入是ce和空间坐标y，它首先从空间位置构建相邻矩阵和相似性矩阵，然后用两个图注意层进行学习并输出n
×
m维的se特征u。
[0098]
步骤s4具体为：
[0099]
首先根据各个细胞的位置构建对称的邻接矩阵a，节点的集合写成v＝{v1，v2，...，vn}，那么邻接矩阵a的每个元素a
ij
可以表示为：
[0100][0101]
式中，是节点vj的近邻节点集合，通过knn或距离截断来计算；
[0102]
构建相似性矩阵∑，令细胞之间的空间区域相关特征相关性随着距离的增加而减少，并趋于指数衰减，所以理想情况下∑是基于高斯核的n
×
n维空间相关矩阵：
[0103][0104]
式中，γ表示带宽参数，控制空间位置对特征的影响程度，γ越大影响越大；
[0105]
所述图注意力网络模块包括两个图注意层，令图注意层的输入是n
×
n维的特征矩阵为h＝(h1，h2，...，hn)，输出为n
×n′
维的矩阵h
′
＝(h
′1，h
′2，...，h
′n)，则：
[0106][0107]
式中，w是n
′×
n维的权重矩阵，是vj邻点的集合，α
ij
是使用softmax函数的归一化注意力系数矩阵：
[0108][0109]eij
＝a
t
(whi||whj)
[0110]
其中a是可学习向量，||是连接操作，在图注意层中使用指数线性单元作为激活函
数σ；
[0111]
空间区域相关特征u经过学习，可以表示为：
[0112]
u＝gat2(gat1(z))
[0113]
式中，gat1()、gat2()分别表示第一层图注意层和第二层图注意层。
[0114]
gat模块学习的目标是使不同位置的u相关性呈现指数衰减，同时保持z中的信息，即n个位置u之间的相关性uu
t
接近∑，u接近z，因此，目标函数可以写成，所以，所述图注意力网络模块的目标函数为：
[0115]
loss2＝λ
spatial
*l
spatial
+λ
rec
*l
rec
[0116]
式中，l
spatial
＝mse(uu
t
，∑)，l
rec
＝mse(u，z)，λ
spatial
和λ
rec
为权重参数。
[0117]
实施例2
[0118]
本实施例在实施例1和实施例2的基础上，提供以下具体实施例：
[0119]
使用以下设定进行数据的处理和分析。首先，不同层分别使用了不同的非线性激活函数。其中，编码器和解码器的f2均使用了线性整流函数(rectified linear unit，简称relu)，解码器的f3和f4分别使用指数函数(exponential function，简称exp)和逻辑函数(logistic function，通常称作sigmoid)。模型损失函数的λ
spatial
和λ
rec
默认设置为1.0(0.5～2)和0.08(0.05～1)；knn计算邻域点的数量为10；γ为距离矩阵的0.05分位值。在训练过程中，我们使用适应性矩估计(adaptive moment estimation，简称adam)优化器来最小化损失函数，学习率为0.001(0.0001～0.01)，网络中随机丢弃(dropout)概率为0.1(0～0.5)，默认迭代次数为50(40～400)。
[0120]
在实施例中，此模型应用于不同平台的空间转录组数据，包括stereo-seq的嗅球(0lfactory bulb)数据、slide-seq的海马体(hippocampus)数据、具有标准分层的starmap的视觉皮层(visual cortex)数据，以及具有复杂组织结构的stereo-seq全脑(brain)数据，验证此模型的准确性和普适性。图3用local inverse simpson’s index(简称lisi)指标(korsunsky，i.et al.“fast，sensitive and accurate integration of single-cell data with harmony(使用harmony快速、灵敏且准确地集成单细胞数据).”nat methods 16，1289-1296，2019)评估了区域划分效果。lisi用于衡量空间区域的区域聚集性，取值范围为[1，∞)，越低的lisi值代表越高的区域聚集程度。本发明(sece)在所有数据中都具有最低的lisi取值和最小的标准差，表明sece在所有数据中均有最优的表现。图4展示了几种算法在具有金标准分层的cortex数据中的区域划分结果。利用指标adjusted rand index(简称ari)、normalized mutual information(简称nmi)、accuracy(简称acc)评估区域划分的准确性。sece的ari、nmi和acc均达到最高，表明sece得到了与已知区域最相近的区域划分结果。以上证据显示sece比现有方案准确性更好，普适性更好。
[0121]
相同或相似的标号对应相同或相似的部件；
[0122]
附图中描述位置关系的用语仅用于示例性说明，不能理解为对本专利的限制；
[0123]
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。

技术特征：
1.一种识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，包括以下步骤：s1：通过空间转录组测序方法获取定量的基因表达矩阵和空间坐标位置矩阵；s2：对所述基因表达矩阵进行标准化，得到标准化后的基因表达矩阵；s3：利用结合零膨胀负二项分布假设的自编码器模块进行降噪低维特征的学习，所述自编码器模块的输入为标准化后的基因表达矩阵，所述自编码器模块的输出为细胞类型相关特征，通过对细胞类型相关特征进行聚类得到每个细胞的细胞类型标签；s4：将细胞类型相关特征和空间坐标位置矩阵作为图注意力网络模块的输入，图注意力网络模块输出空间区域相关特征，通过对空间区域相关特征进行聚类得到每个细胞的空间区域标签。2.根据权利要求1所述的识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，步骤s1中基因表达矩阵和空间坐标位置矩阵，具体为：基因表达矩阵x中的每行代表一个测序基本单元，包含单个细胞或多个细胞，一共有n个细胞，基因表达矩阵x中的每列代表一个基因，一共有g个基因；空间坐标位置矩阵y的每行代表一个细胞，空间坐标位置矩阵y的列是该细胞的二维或三维位置信息。3.根据权利要求2所述的识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，步骤s2中对所述基因表达矩阵进行标准化，具体为：式中，x
n
表示细胞n的基因表达向量，为所述基因表达矩阵的第n行，s
n
表示细胞n的文库大小，s
f
表示标准化后所有细胞的文库大小，表示标准化后的细胞n的基因表达向量。4.根据权利要求3所述的识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，所述基因表达矩阵中的每个元素都遵循零膨胀负二项分布，其参数为(π
ng
,r
ng
,p
g
)，其中π
ng
是观察到真实基因表达值为0的概率，(r
ng
,p
g
)是负二项分布的标准参数，x
ng
是在达到r
ng
次失败时的成功次数，p
g
是每个伯努利试验的失败概率。5.根据权利要求4所述的识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，x
ng
的似然函数为：其中，δ0()表示狄拉克函数，γ()表示伽马函数。6.根据权利要求5所述的识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，步骤s3中自编码器模块包括两层编码器和两层解码器，表示为：z
′
＝f2(z)r
′
＝f3(z
′
)π＝f4(z
′
)式中，f1是编码器，编码器包含两层，分别将输入数据从g维降维至m
′
维、从m
′
维降维至m
维，表示标准化后的基因表达矩阵，z为降维后的细胞类型相关特征矩阵；f2、f3、f4分别为解码器的三个输出，解码器包括两层，分别将输入数据从m维升维至m
′
维、m
′
维升维至g维，z
′
为解码器的第一个输出，f3、f4分别作用于z
′
以学习参数r
′
和参数π，参数r和参数π分别为n
×
g维的矩阵，分别由元素r
ng
和π
ng
组成，r
′
由元素r
′
ng
组成，r
ng
＝s
′
n
r
′
ng
，其中s
′
n
是细胞n的缩放系数，s
′
n
＝exp(log
10
s
n
)。7.根据权利要求6所述的识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，由元素l
g
组成的参数向量l是通过指数变换从可学习向量获得的，具体为：式中，l
g
是负二项分布的logit。8.根据权利要求7所述的识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，所述自编码器模块的损失函数定义为负二项分布的负对数似然函数：loss1＝nll
zinb
(x；π,r,p)式中，nll
zinb
()表示负二项分布的负对数似然函数。9.根据权利要求8所述的识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，步骤s4具体为：首先根据各个细胞的位置构建对称的邻接矩阵a，节点的集合写成v＝{v1,v2,
…
,v
n
}，那么邻接矩阵a的每个元素a
ij
可以表示为：式中，是节点v
j
的近邻节点集合，通过knn或距离截断来计算；构建相似性矩阵σ，令细胞之间的空间区域相关特征相关性随着距离的增加而减少，并趋于指数衰减，所以理想情况下σ是基于高斯核的n
×
n维空间相关矩阵：式中，γ表示带宽参数；所述图注意力网络模块包括两个图注意层，令图注意层的输入是n
×
n维的特征矩阵为h＝(h1,h2,
…
,h
n
)，输出为n
×
n
′
维的矩阵h
′
＝(h
′1,h
′2,
…
,h
′
n
)，则：式中，w是n
′×
n维的权重矩阵，是v
j
邻点的集合，α
ij
是使用softmax函数的归一化注意力系数矩阵：e
ij
＝a
t
(wh
i
||wh
j
)
其中a是可学习向量，||是连接操作，在图注意层中使用指数线性单元作为激活函数σ；空间区域相关特征u经过学习，可以表示为：u＝gat2(gat1(z))式中，gat1()、gat2()分别表示第一层图注意层和第二层图注意层。10.根据权利要求9所述的识别空间转录组空间区域和细胞类型的特征学习方法，其特征在于，所述图注意力网络模块的目标函数为：loss2＝λ
spatial
*l
spatial
+λ
rec
*l
rec
式中，l
spatial
＝mse(uu
t
,σ)，l
rec
＝mse(u,z)，λ
spatial
和λ
rec
为权重参数。

技术总结
本发明公开一种识别空间转录组空间区域和细胞类型的特征学习方法，融合统计分布假设与深度学习的模型框架，同时识别空间转录组数据的空间区域和细胞类型。本发明适用于多种空间转录组数据，包括各种测序平台产生的具有不同分辨率和通量的数据。尤其是对于稀疏程度极高，且具备高通量和单细胞分辨率的空间转录组数据，效果提升非常大。同时，本发明还保留了明确的生物学信号，可以实现良好的低维特征可视化和正确的轨迹推断。化和正确的轨迹推断。化和正确的轨迹推断。


技术研发人员：谢志 尉园园 何尧
受保护的技术使用者：中山大学中山眼科中心
技术研发日：2023.02.15
技术公布日：2023/9/13