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

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,一经查实,立即删除!

相关文章

auto与decltype

auto: 1.定义&#xff1a; 在C中&#xff0c; auto 是一个类型说明符&#xff0c;它让编译器在编译阶段自动推导变量的类型&#xff0c;其类型取决于初始化表达式的类型。auto 在声明变量时使用&#xff0c;编译器会根据变量初始化表达式自动推断类型。 #include<iostrea…

[Code]R2U-Net中的眼部血管分割

DenseUnet.py import torch import torch.nn as nn import torch.nn.functional as F# 定义一个名为Single_level_densenet的类,继承自nn.Module,它构建了一个单层级的DenseNet结构 class Single_level_densenet(nn.Module):def __init__(self, filters, num_conv=4):super…

Java中的“封装“详解

封装&#xff08;Encapsulation&#xff09;是面向对象编程&#xff08;OOP&#xff09;的四大基本特性之一。它通过将数据和操作数据的方法绑定在一起&#xff0c;并隐藏对象的内部实现细节&#xff0c;只提供有限的访问接口来实现。这种机制不仅提高了代码的安全性&#xff0…

深度学习常用方法(一)

1. Dropout 的原理 Dropout 是一种防止神经网络过拟合&#xff08;学习得过于复杂&#xff0c;导致泛化能力差&#xff09;的方法。 原理&#xff1a;在每次训练时&#xff0c;随机“丢弃”一部分神经元&#xff08;即暂时让它们失效&#xff0c;设置为零&#xff09;&#x…

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

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

【Oracle】个人收集整理的Oracle常用SQL及命令

【建表】 create table emp( id number(12), name nvarchar2(20), primary key(id) ); 【充值一】 insert into emp select rownum,dbms_random.string(*,dbms_random.value(6,20)) from dual connect by level<101; 【充值二】 begin for i in 1..100 loop inser…

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

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

Android 编译和使用libheif

项目中需要使用libheif,libde265,libyuv。一下是相应的cmakelist.txt。这里直接使用了静态库。 里面涉及到c包的链接&#xff0c;需要stdc。 ${PROJECT_SOURCE_DIR}/../jniLibs/${ANDROID_ABI}/liblibde265.a这个路径由于操作过程中copy出现问题&#xff0c;多了一层路径&…

【无标题】多模态对齐

第一章 绪论 1.1 研究背景和意义 在现实生活中,人类往往是以事件为单位来认识和理解世界的,一个事件包 含了时间、地点、参与者等多个静态因素。将这些事件信息从非结构化数据中提 取出来,已逐渐成为信息抽取领域的关键任务之一。事件通常并非孤立存在,而 是彼此关联,它…

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中可以直接用类的名字创建数组&#xff0c;如: M…

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

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

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

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

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

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

恶补英语初级第12天,《描述进行中的动作》

对话 Where’s Sally, Jack? She’s in the garden, Jane. What’s she doing? she’s sitting under the tree. Is Tim in the garden, too? Yes, he is. He’s climbing the tree. I beg your pardon? Who’s climbing the tree? Tim is. What about the dog? The dog’…

联表查询,外键

联表查询 使用多表可以减少冗余&#xff0c;便于增删改操作&#xff0c;查询操作复杂 表与表之间通过主键和外键关联 外键 如果一张表的某个列指向另一个表的主键列 &#xff0c;就称之为外键列 -- 添加外键 alter table t_emp add constraint 外键名 foreign key (depart)…

【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;可以访问控制台地址 访问后的…