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

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.mzph.cn/bicheng/62518.shtml

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

相关文章

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

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

RAG (Retrieval Augmented Generation) 检索增强和生成

1 RAG技术简介 1.1 RAG技术概述 RAG(Retrieval Augmented Generation) 是一种结合了检索(Retrieval)和生成(Generation)的技术,旨在通过利用外部知识库来增强大型语言模型(LLMs&am…

DIY-Tomcat part 3 实现对动态资源的请求

实现ServletRequest package connector;import javax.servlet.RequestDispatcher; import javax.servlet.ServletInputStream; import javax.servlet.ServletRequest; import java.io.BufferedReader; import java.io.IOException; import java.io.InputStream; import java.i…

黑马程序员Java笔记整理(day05)

1.面向对象编程 2.用法 3.对象是什么 4.对象在计算机中是啥 5.无参与有参构造器 小结: 6.this的作用 7.小结 8.封装 9.小结 10.实体类 11.小结 12.static 13.小结 14.static修饰方法 15.static应用前景 16.几个注意事项 17.java中可以直接用类的名字创建数组,如: M…

Perforce SAST专家详解:自动驾驶汽车的安全与技术挑战,Klocwork、Helix QAC等静态代码分析成必备合规性工具

自动驾驶汽车安全吗?现代汽车的软件包含1亿多行代码,支持许多不同的功能,如巡航控制、速度辅助和泊车摄像头。而且,这些嵌入式系统中的代码只会越来越复杂。 随着未来汽车的互联程度越来越高,这一趋势还将继续。汽车越…

《Python基础》之数据加密模块hashlib的用法

目录 一、简介 二、用法 步骤一、导入hashlib库 步骤二、创建哈希对象 步骤三、往哈希对象中传值 1、可以在创建对象的时候传值 2、使用updata传值 步骤四、获取经过哈希对象加密后的值 三、注意事项 1、编码问题 2、安全性 3、多次传值 四、总结 一、简介 hashli…

QT实战--qt各种按钮实现

本篇介绍qt一些按钮的实现,包括正常按钮;带有下拉箭头的按钮的各种实现;按钮和箭头两部分分别响应;图片和按钮大小一致;图片和按钮大小不一致的处理;文字和图片位置的按钮 效果图如下: 详细实现…

【OJ】前K个高频单词和单词识别和两个数组的交集

个人主页 &#xff1a; zxctscl 如有转载请先通知 题目 1. 692. 前K个高频单词1.1 分析1.2 代码 2. KY264 单词识别2.1 分析2.2 代码 3. 349. 两个数组的交集3.1 分析3.2 代码 1. 692. 前K个高频单词 1.1 分析 先试用map来统计每个单词出现的次数&#xff1a; map<string,i…

分布式协同 - 分布式锁一二事儿

文章目录 导图Pre概述概述1. 分布式互斥和临界资源的协调2. 分布式锁的基本原理3. 分布式锁的实现方式a. 基于数据库实现的分布式锁b. 基于Redis实现的分布式锁c. 基于Zookeeper实现的分布式锁 4. 高并发场景下的分布式锁优化a. 分段锁&#xff08;Sharded Locks&#xff09;b.…

Vue-01

Vue框架 Vue官网&#xff1a; Vue.js 框架 数据模型和view的通信就是依靠viewmodel的关键。 目前主流版本仍然是vue2版本。 Vue快速入门 1.新建一个HTML文件&#xff0c;引入Vue.js文件。Vue.js文件是官方引入的一个文件&#xff0c;我们如果要使用Vue就必须引入这个文件。…

【微服务】Nacos

一、安装 1、官网地址&#xff1a;https://nacos.io/download/nacos-server/ 2、启动&#xff1a;找到bin目录下的startup.cmd双击启动&#xff0c;或者打开一个命令窗口输入&#xff1a; startup.cmd -m standalone双击启动后如下&#xff1a;可以访问控制台地址 访问后的…

2025 - AIDD - python的autodock vina 批量分子对接改进版本2.0-全自动对接,完全全自动对接

2025 - AIDD - python的autodock vina 批量分子对接改进版本2.0-全自动对接&#xff0c;完全全自动对接 import warnings from pathlib import Path import subprocess from itertools import product import numpy as np import pandas as pd from MDAnalysis import Univers…

【uniapp】轮播图

前言 Uniapp的swiper组件是一个滑块视图容器组件&#xff0c;可以在其中放置多个轮播图或滑动卡片。它是基于微信小程序的swiper组件进行封装&#xff0c;可以在不同的平台上使用&#xff0c;如微信小程序、H5、App等。 效果图 前端代码 swiper组件 <template><vi…

由于导包而引发的错误

今天在调试时发现删除功能无论如何都无法实现&#xff0c;于是调试找到了mapper层的错误但不知道为什么报错。以下是报错信息。 Caused by: org.apache.ibatis.binding.BindingException: Parameter userIds not found. Available parameters are [arg0, collection, list]at o…

STM32 ADC --- 知识点总结

STM32 ADC — 知识点总结 文章目录 STM32 ADC --- 知识点总结cubeMX中配置注解单次转换模式、连续转换模式、扫描模式单通道采样的情况单次转换模式&#xff1a;连续转换模式&#xff1a; 多通道采样的情况禁止扫描模式&#xff08;单次转换模式或连续转换模式&#xff09;单次…

Web day02 Js Vue Ajax

目录 1.javascript: 1.js的引入方式&#xff1a; 2.js变量 & 数据类型 & 输出语句&#xff1a; 模板字符串&#xff1a; 3.函数 & 自定义对象&#xff1a; 4. json 字符串 & DOM操作&#xff1a; 5. js事件监听&#xff1a; 6.js的模块化导入或者导出&a…

Kylin Server V10 下 RocketMQ 主备自动切换模式部署

一、NameServer简介 NameServer 是一个注册中心,提供服务注册和服务发现的功能。NameServer 可以集群部署,集群中每个节点都是对等的关系,节点之间互不通信。 服务注册 Broker 启动的时候会向所有的 NameServer 节点进行注册,注意这里是向集群中所有的 NameServer 节点注册…

芯片测试-RF中的S参数,return loss, VSWR,反射系数,插入损耗,隔离度等

RF中的S参数&#xff0c;return loss, VSWR&#xff0c;反射系数&#xff0c;插入损耗&#xff0c;隔离度 &#x1f4a2;S参数&#x1f4a2;&#x1f4a2;S11与return loss&#xff0c;VSWR&#xff0c;反射系数&#x1f4a2;&#x1f4a2;S21&#xff0c;插入损耗和增益&#…

一个开源轻量级的服务器资源监控平台,支持告警推送

大家好&#xff0c;今天给大家分享一款开源的轻量级服务器资源监控工具Beszel&#xff0c;提供历史数据记录、Docker容器统计信息监控以及多种警报功能&#xff0c;用于监控服务器资源。 项目介绍 Beszel由hub&#xff08;中心服务器端应用&#xff0c;基于PocketBase构建&…

linux centos nginx编译安装

编译安装nginx&#xff08;Centos&#xff09; 编译需要的基础环境yum -y install pcre-devel openssl openssl-devel gd-devel gcc gcc-c1.下载nginx源码包 Nginx源码包下载地址&#xff1a;nginx源码包下载 2. 上传nginx源码包到服务器 我上传的地址是/home/chenhao/nginx…