Spanve:一种检测大规模空间转录组学数据中空间变异基因的有效统计方法
文章题目
Spanve: an Effective Statistical Method to Detect Spatially Variable Genes in Large-scale Spatial Transcriptomics Data
中文名:
Spanve:一种检测大规模空间转录组学数据中空间变异基因的有效统计方法
文章地址:
https://doi.org/10.1101/2023.02.08.527623
本文主旨
提出了Spanve(空间邻域可变表达基因),这是一种基于统计的方法,用于从空间转录组学数据中检测空间依赖性表达基因,通过将空间和表达的依赖性建模为两个分布的距离,该方法具有较高的计算效率和准确性。Spanve能够输入空间转录组并改进空间组织区域的识别。
出发点
作者认为 判断空间转录组学数据的空间模式,需要找到 HVG 高可变基因 或者说 找到SV spatially variable gene 。因为SV基因可由细胞类型定位或空间依赖性细胞相互作用引起,因此可用于多个下游分析,例如重建组织结构。SV基因也有助于提供生物学见解,因为不同的空间表达模式反映了不同的基因功能。
目前的找到sv的方法 取决于 使用协方差矩阵或者说使用聚类方法 他们都不适用于大规模数据集。
此外,基因表达建模技术经常依赖于高斯分布或泊松分布,这可能导致较高的假阳性率,因为并非所有基因表达都符合这些假设。
Spanve这个方法考虑到了 上述所说的问题 不考虑将细胞内的表达建模为传统的统计类分布(经典的几种分布)而转换使用细胞与空间邻居之间的差异性。
由于可能存在位置效应,包括细胞类型分布、局部环境和细胞相互作用,简单分布可能无法描述基因表达,而复杂分布则大大增加了参数的数量和拟合时间。没有基因表达分布的假设,Spanve能够准确地识别具有空间表达模式的基因
此外,通过考虑细胞与空间之间的作用,spanve可以感知空间domain
此外 ,还采用了空间插补方法。

评价指标
使用Moran's I指数:

https://www.jianshu.com/p/305c608f2884?ivk_sa=1024320u
使用Geary's C指数:
这个其实和上面的莫兰指数差不多:

二者区别和联系

moran 和 geary本身都代表空间自相关性
Spanve在 Moran 和 Geary上的效果都是最好的:

CH index
CH指数越高,SV基因作为特征的数据越有可能获得基于原始数据的聚类标记,这意味着SV基因更能代表原始数据。

结果

插补的引入
SPANGE同时引入了插补的概念!这个用于检测生成差异表达基因(SV)的工具同时可以用于插补 插补的方法如下

结果说明

插补后数据的AMI和ARI有着显著的提升

Spanve插补后聚类有了很大改善,可以获得更清晰的组织层边界
总结
Spanve作为一种新开发的空间可变基因(SV)识别方法,该方法可提供可缩放的集群感知基因。本文作者估计,随着空间转录组学数据规模的扩大,可扩展性将变得更加重要。与当前方法相比,Spanve在大规模空间转录组学数据中需要最少的时间和计算成本。此外,Spanve鉴定的基因在空间上更加异质,有利于进一步聚类。Spanve作者团队同时还创建了一种空间插补方法,该方法允许使用流行的聚类方法来识别更准确的组织结构。
限制和不足
一个限制是它不能同时处理来自类似组织切片的数据,以检测稳健的SV基因并利用先前的数据。
Spanve:一种检测大规模空间转录组学数据中空间变异基因的有效统计方法的更多相关文章
- ComplexBrowser: a tool for identification and quantification of protein complexes in large-scale proteomics datasets(大规模蛋白组学数据集中鉴定和定量蛋白复合物)
文献名:ComplexBrowser: a tool for identification and quantification of protein complexes in large-scale ...
- Nature Methods | 新软件SAVER-X可对单细胞转录组学数据进行有效降噪
图片来源(Nature Methods) 摘要 单细胞转 ...
- 4种检测是否支持HTML5的方法,你知道几个?
4种检测是否支持HTML5的方法,你知道几个? 1,检查特定的属性是否存在于全局的对象里面,比如说window或navigator. 比如geolocation,它是HTML5新加支持的新特性:它是由 ...
- Spring MVC内置支持的4种内容协商方式【享学Spring MVC】
每篇一句 十个光头九个富,最后一个会砍树 前言 不知你在使用Spring Boot时是否对这样一个现象"诧异"过:同一个接口(同一个URL)在接口报错情况下,若你用rest访问,它 ...
- 关于ubuntu系统boot分区空间不足而又无法卸载旧内核的解决方法
2016年09月03日 14:16:45 萧氏一郎 阅读数:7802 标签: ubuntuboot分区清理更多 个人分类: linux 版权声明:本文为本猿原创文章,转载务必注明出处,多谢. ht ...
- Oracle创建删除用户,角色,表空间,导入导出数据库命令总结(转载)
无意间看到一篇文章,觉得对于ORACLE的新手很实用,特转载,原文出处这里 说明:在创建数据库时输入的密码,是修改系统默认的密码,以system和sysman等系统默认身份登录时要输入的密码就是修改后 ...
- 表空间tablespace,数据文件datafiles,和控制文件control files介绍
https://docs.oracle.com/cd/B19306_01/server.102/b14220/physical.htm#i5919 本文系翻译 表空间tablespace,数据文件da ...
- [20170623]利用传输表空间恢复部分数据.txt
[20170623]利用传输表空间恢复部分数据.txt --//昨天我测试使用传输表空间+dblink,上午补充测试发现表空间设置只读才能执行impdp导入原数据,这个也很好理解.--//这样的操作模 ...
- 5 种使用 Python 代码轻松实现数据可视化的方法
数据可视化是数据科学家工作中的重要组成部分.在项目的早期阶段,你通常会进行探索性数据分析(Exploratory Data Analysis,EDA)以获取对数据的一些理解.创建可视化方法确实有助于使 ...
- EXPDP/IMPDP与EXP/IMP在不同用户和表空间之间迁移数据的实现方法
1. EXPDP/IMPDP方式 SQL> create user zlm identified by zlm; User created. SQL> grant connect,reso ...
随机推荐
- css样式实现平行四边形
强大的css样式实现平行四边形: 啥也不说了,直接上代码 <!DOCTYPE html> <html lang="en"> <head> < ...
- 网络监测工具之Zabbix的搭建与测试方法(三) ---Zabbix Agent
安装客户端 在官方网站下载最新版本zabbix agent v6.2.6,然后默认安装,其中配置服务端的界面如下图: 其他一律默认即可. 启用发现功能 如上图所示,开启发现规则,默认搜索zabbix ...
- js任务队列EventLoop
JS 执行机制 在我们学js 的时候都知道js 是单线程的如果是多线程的话会引发一个问题在同一时间同时操作DOM 一个增加一个删除JS就不知道到底要干嘛了,所以这个语言是单线程的但是随着HTML5到来 ...
- SICP:符号求导、集合表示和Huffman树(Python实现)
绪论 到目前为止,我们已经使用过的所有复合数据,最终都是从数值出发构造起来的(比如我们在上一篇博客<SICP 2.2: 层次性数据和闭包性质(Python实现)>所介绍的链表和树就基于数来 ...
- TensorRT基础笔记
一,概述 TensorRT 是 NVIDIA 官方推出的基于 CUDA 和 cudnn 的高性能深度学习推理加速引擎,能够使深度学习模型在 GPU 上进行低延迟.高吞吐量的部署.采用 C++ 开发,并 ...
- 初探富文本之OT协同实例
初探富文本之OT协同实例 在前边初探富文本之OT协同算法一文中我们探讨了为什么需要协同.为什么仅有原子化的操作并不能实现协同.为什么要有操作变换.如何进行操作变换.什么时候能够应用操作.服务端如何进行 ...
- [Codeforces Round #794 (Div. 2)] D. Linguistics
我是什么东西艹艹艹 <我离正解只差个sort> 首先,观察字符串,可以发现:若存在形似\(AA--BB\)或\(BB--AA\)等有两个相同的字符挨在一起的情况,则我们在它们中间放一块隔板 ...
- 【ASP.NET Core】动态映射MVC路由
ASP.NET Core 中的几大功能模块(Razor Pages.MVC.SignalR/Blazor.Mini-API 等等)都以终结点(End Point)的方式公开.在HTTP管道上调用时,其 ...
- freeswitch号码黑名单
概述 freeswitch是一款简单好用的的VOIP开源软交换平台. 在客户的呼叫过程中,会遇到一些异常的号码,包括高投诉风险号码.敏感号码.特殊号码等. 在客户呼叫流程中,需要针对这些异常号码做呼叫 ...
- 如何用Python实现http客户端和服务器
功能:客户端可以向服务器发送get,post等请求,而服务器端可以接收这些请求,并返回给客户端消息. 客户端: #coding=utf-8import http.clientfrom urllib i ...