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 ...
随机推荐
- snprintf拼接字符串
例如编辑一个txt文档,不断将字符输入,最终形成一个长句子.可以看成是字符串的不断拼接.snprintf函数具有这个功能. #include<stdio.h> void main(void ...
- [Codeforces Round #841]
[Codeforces Round #841] Codeforces Round #841 (Div. 2) and Divide by Zero 2022 A. Joey Takes Money J ...
- .NET周报【12月第4期 2022-12-31】
祝大家新年快乐! 国内文章 『 再看.NET7』数值类型 https://mp.weixin.qq.com/s/ctiBMPY6Hditk81AzHSRng 在C#中,有int16,用short来定义 ...
- [seaborn] seaborn学习笔记10-绘图实例(2) Drawing example(2)
文章目录 10 绘图实例(2) Drawing example(2) 1. Grouped violinplots with split violins(violinplot) 2. Annotate ...
- 递归实现指数型枚举 (n个可选可不选)
递归实现指数型枚举 从 1∼n 这 n 个整数中随机选取任意多个,输出所有可能的选择方案. 输入格式 输入一个整数 n. 输出格式 每行输出一种方案. 同一行内的数必须升序排列,相邻两个数用恰好 1 ...
- P8855 [POI2002]商务旅行
简要题意 给出一个 \(N\) 个节点的树和一个长度为 \(M\) 的序列 \(S\).你需要从 \(1\) 出发,依次经过 \(S\) 中的所有点,求至少需要经过的边数. \(1 \le N \le ...
- 自定义alert弹框,去掉IP以及端口号提示
最新版例子~~ 如果同时多个弹框,只显示第一个 <!DOCTYPE html> <html lang="en"> <head> <met ...
- antDesign 【NG-ZORRO、Angular】日期选择框时间段nz-range-picker设置默认选择日期及限制日期可选范围
下面的代码包含 1.只可以选择今天以后 2.只可以选择今天开始一周内 3.只能选择今天之前的 import { Component } from '@angular/core'; import dif ...
- 1.31 wlx 魔怔 9 解法交互题小结
参考题解地址 1. 从树上任意一个节点开始,跳到其随机一个后代,跳到叶子的期望次数为 \(H_{siz_u}=\ln(siz_u)\). 证明: 首先考虑一条链的情况.设在第 \(i\) 个点期望次数 ...
- MRS_下载相关问题汇总
解决问题如下: MRS下载编译时,更改生成文件是HEX文件还是BIN文件 关于MounRiver下载时如何选择配置部分擦除 关于MounRiver下载起始地址配置问题 MRS下载编译时,更改生成文件是 ...