单细胞细胞通讯全流程分析教程,代做分析和辅导

0. 分析参数文件和细胞通讯的演示数据

0.1 细胞通讯分析总的参数文件,后面部分细胞通讯分析模块会用到

分析参数文件

参数文件名称:total_analysis_params_demo.xlsx ,很多分析模块都是这个总的参数文件,我的这个总的参数文件如果有更新的话,我就上传到百度网盘。包括拟时序分析的所有分析模块都是用的这一个参数文件

从百度网盘下载该参数文件,该参数文件在百度网盘的位置:

32708717df303050ef27689cbd5c912b.jpeg

我们可以把这个参数文件下载到本地的:D:/omics_tools/demo_data这个目录下。下载后,在本地的文件为:

d2db4d9b58c59185e1397359054d9e9a.jpeg

拟时序分析所有模块都用的这个总的参数文件的完整的文件路径就是:D:/omics_tools/demo_data/total_analysis_params_demo.xlsx

细胞通讯分析的演示数据集

bb7dd2e2058150cf4e238d8225a5a15b.jpeg

0.2 分析思路

1b11ce759560498a564be04087791fc9.jpeg

0.3 每个模块下的教学视频

bf3ca0d25e6e7fd164928f88e6d1d201.jpeg

每个分析模块下都放了一个对应的b站教学视频链接,大家都可以看该模块下对应的b站教学视频来学习该模块的使用方法

1.细胞通讯分数计算和差异贡献度通路筛选

1.1对每个分组单独进行所有类型细胞的细胞通讯分析计算

分析模块位置

c05e051a3777da6411f8ab3ba6348096.jpeg

软件界面

33f25f5a4f963761be64f247cc514d03.jpeg

使用我OmicsTools前面单细胞的seurat系列分析注释好细胞名称的rds文件来做这一步的细胞通讯分析,这一步分析耗时比较长,可能会要运行1-2个小时。

运行完成的结果

484e9ba3e42496d15a0931ebec1b2b58.jpeg

这里我用了Tumor和Normal两个分组,每个分组都使用CellChat计算处理细胞通讯网络的两个rds文件,后续的分析需要用到这一步的结果时,直接给出这里细胞通讯分析的这个结果目录的目录路径就可以了。

9245979fa2198a3326216e6920d7e194.jpeg

这是CellChat细胞通讯数据库的信号通路类别占比百分比的汇总图,展示了每种类别的包含的这些信号通路的百分比情况。

1.2对部分感兴趣细胞的细胞通讯分析计算(这一步可以放到后面再做,先把总的细胞通讯做的差不多了,找到感兴趣的细胞再做这一步)

软件界面

c721ed759c020b8e9cb611269da9e4c7.jpeg

使用我OmicsTools前面单细胞的seurat系列分析注释好细胞名称的rds文件同时给出一组感兴趣的细胞名称来做这一步的细胞通讯分析,这一步分析耗时比较长,可能会要运行1-2个小时。

一组感兴趣的细胞名称的示例文件sub_cell_type.csv内容如下:

3a6adc5da5c7713c58c0786a767f97ea.jpeg

运行完成的结果

2e7dfed9a8b65d6581b1066fcf4e4158.jpeg

这里我用了Tumor和Normal两个分组和一组感兴趣的细胞,每个分组都使用CellChat计算处理特定的一组细胞间的细胞通讯网络的两个rds文件,后续的分析需要用到这一步的结果时,直接给出这里细胞通讯分析的这个结果目录的目录路径就可以了。

14a32813883681ab518ed72e74e3c549.jpeg

这是CellChat细胞通讯数据库的信号通路类别占比百分比的汇总图,展示了每种类别的包含的这些信号通路的百分比情况。

1.3对每个分组总的细胞通讯可视化和lr贡献度可视化

软件界面

432d5704c81d99857078af64d530eb66.jpeg

该步的分析结果目录主要是存放在func_save_dir指定的目录路径下

分析结果

结果文件列表

67e361e28db28c30c5a98123d989be5f.jpeg

细胞通讯网络的介绍

CellChat是一个用于分析和可视化细胞间通讯网络的工具。在CellChat的输出中,interaction numbers和interaction weights是两个重要的指标,它们分别代表:

1. Interaction numbers (交互数量):

这个指标表示在细胞间通讯网络中检测到的显著交互的数量。具体来说:

  • 它反映了不同细胞类型之间可能发生的配体-受体对的数量。
  • 每一个检测到的配体-受体对都被计为一个交互。
  • 较高的交互数量可能意味着两种细胞类型之间有更多样化的通讯方式。

2. Interaction weights (交互权重):

这个指标衡量了细胞间交互的强度或重要性。具体含义是:

  • 它通常基于配体和受体的表达水平来计算。
  • 权重值越高,表示该交互在细胞通讯网络中可能越重要。
  • 它考虑了细胞数量、基因表达水平等因素,给出了一个综合的交互强度评分。

这两个指标结合起来可以帮助研究者更全面地理解细胞间通讯的模式:

  • Interaction numbers告诉我们通讯的配体-受体对的数量
  • Interaction weights则反映了这些通讯的相对重要性

在分析CellChat结果时,通常需要同时考虑这两个指标,以获得细胞间通讯网络的完整图景。

这两个指标的数值范围及其含义:

1. Interaction numbers (交互数量):

数值范围:通常是非负整数,从0开始,上限取决于数据集的特性。

  • 0:表示两种细胞类型之间没有检测到显著的交互。
  • 1-5:较低的交互数量,可能表示细胞类型之间的通讯有限。
  • 6-20:中等水平的交互,表示存在多种通讯途径。
  • 20以上:高水平的交互,表明细胞类型之间有丰富的通讯网络。

具体的数值范围可能因研究而异,取决于研究的规模、细胞类型的多样性以及使用的配体-受体数据库。

2. Interaction weights (交互权重):

数值范围:通常是非负实数,CellChat中常见的范围可能在0到1之间,但也可能超过1,具体取决于权重的计算方法。

  • 0-0.2:弱交互,可能不太重要或不太可靠。
  • 0.2-0.5:中等强度的交互,可能有一定的生物学意义。
  • 0.5-0.8:较强的交互,在细胞通讯网络中可能扮演重要角色。
  • 0.8-1.0或更高:非常强的交互,可能是关键的细胞通讯途径。

没有标记通讯交互数量或通讯强度的细胞通讯总图

0b18bf69dee3c4be62850242bc99ae48.jpeg

标记了细胞通讯交互数量的图

f75588e21641c208999dd25fd04f8234.jpeg

标记了细胞通讯交互强度的图

93d2683d5d58710a51300b82b712f449.jpeg

形成闭环的线是细胞的自分泌通讯,没有形成闭环的线是细胞间的旁分泌通讯。

将待研究的基因跟细胞通讯网络联系起来的结果解读,以CABLES1基因为例

CABLES1基因在肾癌的肿瘤vs正常组织中的单细胞所有细胞类型中的差异表达结果

b9edb60c0b50a96b7ad7387f7c4031a2.jpeg

查看使用seurat之前做的差异分析结果中有该基因表达差异的细胞名称跟细胞通讯的通讯互作网络联合起来解读

a540aed526c31780b552615a9f7b26bc.jpeg

1.4 根据通路中受配体的lr贡献度得分提取出贡献度显著差异的通路

软件界面

2823a319a95b4af7c8c3110c8743f673.jpeg

运行结果

8bc9613ab3da3b7b5c035fffe244c039.jpeg

6f61244d8bd6781d5c97c9eb2f2c6515.jpeg

10.1.4.sig_lr_cellchat.Tumor-Normal.sig_pathway_lr_contribution_res.csv该文件是有差异的通路的汇总文件,这些差异的通讯通路一般是细胞通讯中的贡献度排名最高的受体-配体对在正常或肿瘤间改变了名字变成了其它的受体-配体对或者二者的值差距较大。

1.5对贡献度差异通路可视化

软件界面

467e7b680ce63d6c2459135c83ef9b00.jpeg

运行结果

1dd0b2b1ee11b020f7bf105be8097aeb.jpeg

Nomal中的CCL通路

7f02920c5452c837da97b3d1b2c27e91.jpeg

Tumor中的CCL通路

a9adeccea8fbade4b1a67d119760db53.jpeg

在肿瘤组织中,CCL通路下CCL20-CCR6的受体配体的通讯相较于正常组织贡献度占比明显升高了,由正常组织中的第12位上升到了肿瘤组织中的第2位

对要研究的基因互作的通讯通路受体配体对和细胞进行筛选和可视化

2.1.1 获取某个基因的核心ppi互作蛋白列表

软件界面

fd22eb0f2585d5e38bd70fde71ff654f.jpeg

运行结果

76268fe2e1db74e37f49fc5122871614.jpeg

83ffccb55c7e7a5c48afc84b51555e14.jpeg

得到了string数据库中跟我们想研究的CABLES1的一组互作基因列表

2.1.2查看一组基因列表是否在细胞通讯通路中,跟细胞通讯相关

软件界面

fa16d23c7e8635475f8898b54fbe8b21.jpeg

6ffdfce555833e1ede24d10e100d982e.jpeg

在excel中对分析参数进行修改

该分析模块的名称step_module_name是filter_genelist_in_cellchat,即用的是这个总的参数文件的完整的文件:D:/omics_tools/demo_data/total_analysis_params_demo.xlsx 里面的该模块的参数值。

39a166d39e2d0d4c9e8fd80704e43aca.jpeg

修改参数值的时候,只用改这个模块的参数值中我标红的那几个参数值,包括大家电脑中要分析的数据的文件路径和一些要调整和特别指定的重要参数值,没有标红的参数值,大家一律用默认的,不用修改。

运行结果

c4686243d79d179d29602b1bb8d7a238.jpeg

这个里面是用的CABLES1的一些互作蛋白的基因名称列表跟CellChat细胞通讯数据库中的通路中的受体配体列表取交集,得到跟我们想研究的基因的互作蛋白哪些也是属于细胞通讯通路中的受体配体的结果表格文件。

ppi_genes_in_cellchat_pathway_filter_in_cellchat_genes.csv文件中的结果:

ffd5aeb5d02ac85fcd6874bb9423816f.jpeg

只找到了在细胞通讯通路受体配体中只有一个基因是CABLES1的互作蛋白DLK1。

ppi_genes_in_cellchat_pathway_filter_total_target_interactions.csv文件中的结果:

ae52426d905ade58d3810fd2a8b89371.jpeg

这是DLK1配体所属的NOTCH通路以及该通路中跟DLK1互作通讯的受体配体对情况

2.2将有差异的通路中的受配体跟该基因互作蛋白的结果取交集

软件界面

6fecd407cd4639f3b3615755f04f9b34.jpeg

在excel中对分析参数进行修改

该分析模块的名称step_module_name是extract_common_pathway_genes,即用的是这个总的参数文件的完整的文件:D:/omics_tools/demo_data/total_analysis_params_demo.xlsx 里面的该模块的参数值。

7c6282bda5da51e2c5690abb782b2c0a.jpeg

修改参数值的时候,只用改这个模块的参数值中我标红的那几个参数值,包括大家电脑中要分析的数据的文件路径和一些要调整和特别指定的重要参数值,没有标红的参数值,大家一律用默认的,不用修改。

运行结果

c844311da024bc84cd8c37a1d655527f.jpeg

这是将有差异的通路中的受配体跟该基因互作蛋白的结果取交集的结果文件列表

10.2.2common_pathway_genes_sig_lr_contribution_pathway_target_genes.csv文件中的内容:

026366fee28384d40f44d8927e248c99.jpeg

只有DLK1这一个基因既是CABLES1的互作蛋白,又是细胞通讯贡献度有差异通路中的受体配体对。

10.2.2common_pathway_genes_sig_lr_contribution_ppi_pathway_info.csv文件中的内容:

4c892b3a5bdd82086bd0b62f63964f1a.jpeg

这是交集中贡献度有差异的NOTCH通路以及该通路中跟DLK1互作通讯的受体配体对详细信息

10.2.2common_pathway_genes_target_NOTCH_pathway_gene_list.tsv文件中的内容:

fc1628c58731f29af8b763660a617d1b.jpeg

这是交集中贡献度有差异的NOTCH通路以及该通路中跟DLK1互作通讯的受体配体对详细基因列表

2.3将最后筛选到的存在lr贡献度差异的靶基因细胞通讯通路进行细胞通讯网络的可视化

软件界面

a42194170bb62249ed99d78fcd1e6f0c.jpeg

在excel中对分析参数进行修改

该分析模块的名称step_module_name是sig_lr_pathway_cellchat_plot,即用的是这个总的参数文件的完整的文件:D:/omics_tools/demo_data/total_analysis_params_demo.xlsx 里面的该模块的参数值。

c722f076d867af80713933c2674a9da6.jpeg

修改参数值的时候,只用改这个模块的参数值中我标红的那几个参数值,包括大家电脑中要分析的数据的文件路径和一些要调整和特别指定的重要参数值,没有标红的参数值,大家一律用默认的,不用修改。

运行结果

结果文件列表

16b6d9109c65c581bff0499d2a72692b.jpeg

这些图都是肿瘤组和正常组都是贡献度条形图barplot,通讯网络circle图,通讯网络的和弦图chord图各画了一个图,在解读图片结果的时候,需要把正常组织的结果图和肿瘤组织的结果图对比着来看。

贡献度图比较分析

81fd312b6708ea547ca133859f26b088.jpeg

通讯网络circle图,通讯网络的和弦图chord图比较分析

8a841798de22ba03635fef58386f1e4c.jpeg

2.4 感兴趣的基因的互作蛋白与筛选到的细胞通讯通路的受体配体的ppi互作网络

把2.2步将有差异的通路中的受配体跟该基因互作蛋白的结果取交集基因列表放到String数据库网站中就可以得到这样到互作网络图了

34ccc54cdc479ebf7eb9c3653a8e3266.jpeg

2.5 绘制细胞通讯的heatmap热图来确定发送信号和接受信号的重要细胞

软件界面

bc9f12cbb312f6a97694f676bf9311e3.jpeg

在excel中对分析参数进行修改

该分析模块的名称step_module_name是cellchat_heatmap_plot,即用的是这个总的参数文件的完整的文件:D:/omics_tools/demo_data/total_analysis_params_demo.xlsx 里面的该模块的参数值。

adc074bc5b4947adac428ef7c8bc6360.jpeg

修改参数值的时候,只用改这个模块的参数值中我标红的那几个参数值,包括大家电脑中要分析的数据的文件路径和一些要调整和特别指定的重要参数值,没有标红的参数值,大家一律用默认的,不用修改。

运行完成的结果

dea8ec7fed4b68ecf5452b86417542de.jpeg

把Normal和Tumor的热图像这样放在一起比较

0e9743895b7f34ae19a365865f1ba2b6.jpeg

2.6 对筛选好的细胞通讯通路和感兴趣细胞的细胞通讯进行详细的可视化分析

软件界面

623d904e132fedc169f8f9962ef2b898.jpeg

在excel中对分析参数进行修改

该分析模块的名称step_module_name是sig_pathway_cellchat_show,即用的是这个总的参数文件的完整的文件:D:/omics_tools/demo_data/total_analysis_params_demo.xlsx 里面的该模块的参数值。

5d87c7d3cd880a5c825241c8d76b5009.jpeg

修改参数值的时候,只用改这个模块的参数值中我标红的那几个参数值,包括大家电脑中要分析的数据的文件路径和一些要调整和特别指定的重要参数值,没有标红的参数值,大家一律用默认的,不用修改。

运行完成的结果

结果文件列表

f3a278b1cb05db88be37e0227bf2b630.jpeg

这一步会得到特别多的图

该模块结果中的子步骤cellchat1的结果图

911cf55946e7c86de2ce16c0a1d4a0ca.jpeg

该通路在肿瘤和正常组织中的细胞通讯网络图

该模块结果中的子步骤cellchat2的结果图

5bb58bd35c3453d19ae94651e08d27c2.jpeg

该模块结果中的子步骤cellchat3的结果图

8f5b93ad0f09e88e148e91d9a9066585.jpeg

该通路下每个受体-配体对介导的细胞通讯网络

43a65a5b1bc0dd5f920fff6ae7e9d3e4.jpeg

该模块结果中的子步骤cellchat4的结果图

da7552f1b4f8df92cbee41a5e0c278a5.jpeg

aae431a9a2cda7347483806135c14f4f.jpeg

该模块结果中的子步骤cellchat5的结果图

4fd75dccf61c96c7b0af9c025ac17fa9.jpeg

d4406327f518960365486491401f06b5.jpeg

该通路受体配体在肿瘤和正常组织中的表达水平小提琴图

该模块结果中的子步骤cellchat6的结果图

d066508bdc2e7191159fcea65825936c.jpeg

e5026a2f0a295466378f7ab859fbf886.jpeg

可以提供的科研服务清单

6011e883e2eba5bccdf4e4e0f46b139f.jpeg

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/481148.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Flume 与 Kafka 整合实战

目录 一、Kafka 作为 Source【数据进入到kafka中,抽取出来】 (一)环境准备与配置文件创建 (二)创建主题 (三)测试步骤 二、Kafka 作为 Sink数据从别的地方抽取到kafka里面】 (…

从零开始理解JVM:对象的生命周期之对象销毁(垃圾回收)

一、JVM参数 在学垃圾回收器之前,我们先要知道,jvm参数是怎么回事。因为配置各种回收器,必须对应各种参数设置。 标准参数(-) 所有的JVM实现都必须实现这些参数的功能,而且向后兼容 -help-version 非标准参…

win10中使用ffmpeg的filter滤镜

1 给视频加文字水印 1.1 添加播放时间 ffmpeg -i input.mp4 -vf "drawtextfontfileC\\:/Windows/fonts/consola.ttf:fontsize30:fontcolorwhite:timecode00\:00\:00\:00:rate25:textTCR\::boxcolor0x000000AA:box1:x20:y20" -y output.mp4 在视频的x20:y20位置添加t…

模拟器快速上手,助力HarmonyOS应用/服务高效开发

文章目录 1 创建模拟器1)打开设备管理界面2)设置本地模拟器实例存储路径3)创建一个模拟器(1)选择模拟器设备(2)创建模拟器(3)启动模拟器(4)关闭模…

HarmonyOS(61) 组件间状态共享的分类以及状态选择器的选取优先级

状态共享 状态共享的分类状态共享选择器State与prop\Link\ObservedObjectLink组合的区别合理选择装饰器的顺序参考资料 状态共享的分类 HarmonyOS的组件之间是可以共享状态数据了,不同的组件之间,状态共享的场景也不一样,根据共享范围从小到…

高德地图 Readme GT 定制版 10.25.0.3249 | 极致简洁

这款定制版高德地图去除了广告,运行速度更快。虽然没有车道级导航、打车功能和红绿灯倒计时等功能,但支持正常登录和收藏功能。检测更新始终为最新版本。 大小:82.5M 下载地址: 百度网盘:https://pan.baidu.com/s/1Y…

去中心化物理基础设施网络(DePIN):重塑未来的基石

一、引言:DePIN的定义与背景 什么是DePIN? 去中心化物理基础设施网络(DePIN,Decentralized Physical Infrastructure Networks)是利用区块链和去中心化技术管理、优化和激励物理资源分配的一种新兴模式。与传统集中式…

模型 布鲁姆法则

系列文章 分享 模型,了解更多👉 模型_思维模型目录。分层提升思维力。 1 布鲁姆法则的应用 1.1 布鲁姆法则在产品开发流程中的应用 背景: 在产品开发领域,创新和效率是关键。布鲁姆法则可以帮助产品经理和设计师系统地提升产品开…

恒创科技:服务器操作系统和客户端操作系统之间的区别

客户端操作系统和服务器操作系统是两种不同的操作系统,旨在满足计算机网络环境中的特定目的。虽然每种类型的操作系统在基本功能方面都有一些相似之处,但它们针对不同的用例进行了优化,并具有针对其特定角色量身定制的特定功能。 什么是服务器…

Flink的双流join理解

如何保证Flink双流Join准确性和及时性、除了窗口join还存在哪些实现方式、究竟如何回答才能完全打动面试官呢。。你将在文中找到答案。 1 引子 1.1 数据库SQL中的JOIN 我们先来看看数据库SQL中的JOIN操作。如下所示的订单查询SQL,通过将订单表的id和订单详情表ord…

2024学习之前端微信小程序开发教程,从入门到精通-含基础+实战+源码code

目录 一、简单介绍 二、课程需知 三、内容编排 1、小程序基础  起步式 目录结构 小程序框架 场景值  逻辑层 视图层 组件 视图容器 基础内容 表单组件 导航 媒体组件 Api 路由 界面 交互 网络 数据缓存 自定义组件 2、项目实战 …

HarmonyOS4+NEXT星河版入门与项目实战(23)------组件转场动画

文章目录 1、控件图解2、案例实现1、代码实现2、代码解释3、实现效果4、总结1、控件图解 这里我们用一张完整的图来汇整 组件转场动画的用法格式、属性和事件,如下所示: 2、案例实现 这里我们对上一节小鱼游戏进行改造,让小鱼在游戏开始的时候增加一个转场动画,让小鱼自…

数据预处理方法—特征选择、特征缩放、特征构造

特征选择 1.1 原理 特征选择是选择对模型训练最重要的特征,减少数据维度,去除冗余或不相关特征,提高模型性能的性能和训练速度,减少过拟合。 1.2 核心公式 可以使用基于树模型的特征重要性度量,如在随机森林中计算特…

【C++/Qt 】使用QCustomplot类打造一款数学函数图像生成工具(支持latex公式渲染+Python连接AI大模型)

✨✨ Rqtz 个人主页 : 点击✨✨ 🌈Qt系列专栏:点击 软件介绍 基于Qt的开源项目QCustomplot类的一款在线的数学函数图像生成工具,涉及到了数学的latex公式渲染,如何将latex语法转换为Python的函数,和如何在Qt中使用QCustomplot类进…

分页查询功能

EmployeeController /** * 员工分页查询 * * param employeePageQueryDTO * return */ GetMapping("/page") ApiOperation("员工分页查询") public Result<PageResult> page(EmployeePageQueryDTO employeePageQueryDTO) { log.info("…

mp4视频流推送的学习

一、依赖引入&#xff1a; ①使用 CDN 的播放器代码 <!-- 引入 xgplayer 核心 --> <script src"https://unpkg.byted-static.com/xgplayer/3.0.10/dist/index.min.js" charset"utf-8"></script><!-- 引入 xgplayer mp4 插件 -->…

C++趣味编程:基于树莓派Pico的模拟沙漏-倾斜开关与LED的互动实现

沙漏,作为一种古老的计时工具,利用重力让沙子通过狭小通道,形成了计时效果。在现代,我们可以通过电子元件模拟沙漏的工作原理。本项目利用树莓派Pico、倾斜开关和LED,实现了一个电子沙漏。以下是项目的详细技术解析与C++代码实现。 一、项目概述 1. 项目目标 通过倾斜开关…

PyG教程:MessagePassing基类

PyG教程&#xff1a;MessagePassing基类 一、引言二、如何自定义消息传递网络1.构造函数2.propagate函数3.message函数4.aggregate函数5.update函数 三、代码实战1.图数据定义2.实现GNN的消息传递过程3.完整代码4.完整代码的精简版本 四、总结1.MessagePassing各个函数的执行顺…

Linux—进程学习—04(进程地址空间学习)

目录 Linux—进程学习—41.程序地址空间1.1虚拟地址空间的现象1.2虚拟地址空间的理解(感性) 2.进程地址空间2.0 mm_struct结构体2.1 mm_struct结构体的源代码2.2分页&虚拟地址空间解释前面的实验现象 2.3进程地址空间存在的原因2.3.1第一个原因2.3.2第二个原因2.3.3第三个原…

信息安全实验--密码学实验工具:CrypTool

1. CrypTool介绍&#x1f4ad; CrypTool 1的开源教育工具&#xff0c;用于密码学研究。通过CrypTool 1&#xff0c;可以实现加密和解密操作&#xff0c;数字签名。CrypTool1和2有很多区别的。 2. CrpyTool下载&#x1f527; 在做信息安全实验--密码学相关实验时&#xff0c;发…