拷贝数变异分析的python实现及R语言对比

拷贝数变异学习手册(Python版本)

拷贝数变异(Copy number alterations,以下简称CNA) 是一种重要的基因组变异,在癌症的发生和发展过程中起着至关重要的作用。确定肿瘤细胞中CNA的特征对早期肿瘤检测、划分肿瘤异质性、了解肿瘤进展模式和揭示耐药机制具有重要意义。随着单细胞RNA测序(scRNA-seq)的蓬勃发展,研究人员开发出了多种计算方法来从scRNA-seq研究中推断CNA。目前从scRNA-seq数据中推断CNA的计算方法大致可分为两类。一类是基于表达的,即直接从基因表达推断CNA。这类工具包括inferCNV、copykat和SCEVAN。这些工具的工作原理是,CN扩增(AMPs)或缺失(DELs)将分别导致受影响基因组区域内基因的上调或下调。InferCNV 利用基因窗口的校正移动平均数;copykat 采用综合贝叶斯分割方法;SCEVAN采用多通道分割算法。第二类将基因表达与等位基因信息相结合,相关工具有Numbat和CaSpER等。

往期我们 也基于R版本验证了其中部分软件的功能,感兴趣的小伙伴可以 点击**R语言scRNA-seq拷贝数变异分析** 学习哦~

学习手册部分图片集锦:

优惠订阅资料 或参与课程 请联系**Biomamba_zhushou**

更多生信教程可见: 为什么我们敢开放终身订阅

-Biomamba

所见即所得

以html格式提供学习手册,可在浏览器中打开并翻阅学习:‍

教程目录:

零、欢迎关注

一、前言

1.1 研究背景

1.2 主要结果

1.3 小结

二、准备工作

2.1 环境准备

2.2 库准备

三、infercnvpy

3.1 数据下载

3.2 参数详解

3.3 运行infercnvpy

3.4 计算每个细胞的cnv分数

3.5 可视化

四、通过CNV谱进行聚类并识别肿瘤细胞

4.1 cnv聚类分群

4.2 cnv展示umap

4.3 scanpy展示cnv结果

4.4 肿瘤细胞分类

4.5 肿瘤和正常细胞分开绘制热图

五、内嵌copykat

5.1 数据读入

5.2 参数详解

5.3 运行copycat

5.4 计算每个细胞的cnv分数

5.5 可视化

5.6 通过cnv谱进行聚类并识别肿瘤细胞

5.7 肿瘤和正常细胞分开绘制热图

六、对比R版本的infercnv

七、如何获取基因在染色体上的位置

八、软件版本

九、单细胞教程全收录

十、欢迎致谢

文末下载链接中包含以下内容:

测试文件,总大小10.62GB

往期推送

Biomamba 生信基地现已开放终身订阅,现已包含几十篇订阅资料,未来37年持续更新:为什么我们敢开放终身订阅

2025年我们零差评完成订阅任务(2026订阅须知),2026已更新订阅内容:

单细胞Figure1的常见图表(Python版)

告别R依赖!一文搞定Python单细胞基础分析

风靡一时的scTenifoldKnk单细胞虚拟基因敲除,准不准?

scRNA拟时序分析三大Python包: palantir、PAGA、RNAvelocity

正刊都在用的AI驱动蛋白质从头设计流程:RFdiffusion+ProteinMPNN+AlphaFold3

CellphoneDB搞定细胞通讯分析

scRNA-seq肿瘤细胞药物敏感性预测分析

fastXXX、rapidXXX系列加速你的单细胞数据分析

转录因子、基因共表达网络分析手册:pySCENIC在脑、黑色素瘤中的应用

让你的scRNA数据变得更干净~

保姆级视频教程| 空间转录组分析手册(基于Seurat)

解锁前需阅读:2026订阅须知

相关推荐
caoerzhong18 分钟前
中小企业上 WMS 该先上哪几块:JeeWMS 开源 Java 仓库管理系统的分批上线清单
java·python·开源
正在走向自律33 分钟前
AI数据分析与可视化:从基础到应用实践
服务器·人工智能·python·机器学习·数据分析·pandas
databook1 小时前
从手动检查到自动监控:一个数据质量工作流的实现
后端·python·数据分析
北冥有鱼被烹3 小时前
砷化镓与磷化铟:光模块的地基——从化学键到AI缺货潮的底层逻辑
人工智能·python
for_ever_love__4 小时前
线性回归与梯度下降——从零手写一个模型
python·机器学习·线性回归·梯度下降
code_whiter5 小时前
7.自动化测试常用函数
python·功能测试
小园子的小菜5 小时前
Python协程深度解析:从原理演进到实战避坑
后端·python
砚底藏山河5 小时前
量化实战:截面因子有效性检验(IC 分析与分层回测)
java·python·金融·maven
GitFun6 小时前
给策略加了条-8%止损线,11年只触发1次
python·股票·量化交易
马六六i6 小时前
市面上正规的IP驱动产业新场景新工具有哪些
大数据·网络·python·tcp/ip