临床数据 4. 肿瘤克隆进化分析结果解读?

news2024/10/5 13:17:17

62f993c1dbe29c2d8f0fee396cf3573d.gif


临床数据分析方案

桓峰基因公众号推出临床数据分析方案,整理如下:

临床数据 1. 临床基因突变数据如何发高分?

临床数据 2. 基于NGS的胃癌诊疗全过程的临床应用方案

临床数据 3. 肿瘤微小残留病灶(MRD)如何发文章?


克隆进化生信分析教程

我们在桓峰基因公众号也推出克隆进化分析教程,从组织到单细胞的克隆进化分析都有介绍,供大家参考:

Clone 1. 肿瘤克隆进化之前世今生

Clone 2. 肿瘤克隆进化之不同进化模

Topic 1.克隆进化之sciClone

Topic 2.克隆进化之ClonEvol

Topic 3. 克隆进化之 fishplot

Topic 4. 克隆进化之 Pyclone

Topic 5. 克隆进化之 CITUP

Topic 6. 克隆进化之 Canopy

Topic 7. 克隆进化之 Cardelino

Topic 8. 克隆进化之 RobustClone

Topic 9. 克隆进化之 TimeScape

克隆进化分析报告

  • 准备基因组检测结果

1.    CNV 检测的结果

可以通过多种软件获得,来源一般三种方式:

  1. WGS

  2. WES

  3. Panel

Note: 没有也没办法

CNV 格式要求必须有这4列,如下:

1. chromosome

2. segment start position

3. segment stop position

4. copy number value for that segment.

Note: Unrepresented regions are assumed to have a copy number of 2

举例说明,如下:

Chrstartend value

11    315000 1.293159

1315001  50150000 2.269226

1 5015000151060000 2.521881

1 5106000170690000 2.231479

1 7069000170720000 1.930534

1 70720001 121185000 2.215604

2.    SNV 检测的结果

可以通过多种软件获得,来源一般三种方式:

  1. WGS

  2. WES

  3. Panel

Note: 没有就做不了哦!!!

SNV 格式要求必须有这5列,如下:

1. chromosome

2. position

3. reference-supporting read counts

4. variant-supporting read counts

5. variant allele fraction (between 0-100)

其他:Gene, driver gene.

举例说明,如下:

ChrPos RefDepth AltDepth    VAF     Gene is.dirver

chr12488153        1      352 99.44 TNFRSF14     FALSE

chr18073500        0        0  0.00   ERRFI1     FALSE

chr18073519        0        0  0.00   ERRFI1     FALSE

3.    LOH检测的结果

可以通过多种软件获得,来源一般三种方式:

  1. WGS

  2. WES

  3. Panel

  4. 没有也没办法

LOH格式要求必须有这3列,如下:

1. chromosome

2. window start position

3. window stop position

Note:属于这些窗口的单核苷酸变异将不包括在分析中。

举例说明,如下:

18020922 9987323


  • 生信分析流程

克隆进化整个生信分析流程分成三部分:sciClone + clonevol + fishplot

d6d0fef2b1b5924882192d7caf9fca54.jpeg

1.    sciClone 产出结果解读

  1. a)    两两样本比较聚类结果:

  2. 0374e8cb74a8bf343ea14ca4944058a2.png6357d5bd9de78d61ad74579a2a933edd.png

  3. 3ba92c106f92304a98d9bb060ef3f1f5.png

  1. b)    Table1. 聚类后信息表

  2. c302811ec4dff8d5fe9cdbaf6315f919.png

  1. c)    三维样本结果

  2. 081d6c8b75908e402c0da465665ac5b8.png

2.    clonevol 产出结果解读

  1. a)    绘制进化树

  2. a5801281ff2f59489fa7e1b07418f3ba.png

  1. b)   聚类flow图

  2. 6161333a6413305202d083c3c3e58d92.png

  1. c)    样本聚类两两比较

  2. 4bf9068c14d95fe2d6acdc8e744cc953.png

  1. d)   所有克隆的进化树

  2. 8dbfcdb0ff5af1b5e2f92efe43f014c8.png

  1. e)    样本VAF在不同cluster中的分布

  2. d8687c66f29bdfbb34464670f38137cc.png

  1. f)     综合结果分五部分内容

  2. 1. 每个样本在不同cluster中的突变频率;               

  3. 2. 每个样本的贝尔图;

  4. 3. 每个样本的细胞总体分布;                                  

  5. 4. 基于节点的一致树;

  6. 5.基于分支的一致树。

a892480048bf93dafed65352759b6915.png

3.    fishplot 产出结果解读

结果显示肺癌的进化模式:转移性多克隆支型进化模式

d9535b1d4fda6fb6cad2bdb3b79d6205.png

  • 实例分析

这里我们分享一篇2012年五月发表于 Nature(IF:11.059)。作者基于全基因组测序研究血液病AML复发的克隆进化,该文章使用桓峰基因公众号里面生信分享教程即可实现,有需要类似思路的老师可以联系我们!

806d52c7cb139e103b6b1f6cc6d72271.png

摘       要

大多数急性髓系白血病 (AML) 患者在复发后死于进展性疾病,这与细胞遗传学水平的克隆进化有关1,2。确定突变与复发相关的谱系,我们对原发性肿瘤进行了测序以及来自八名 AML 患者的复发基因组,并经过验证使用深度测序的数百个体细胞突变;这允许我们精确地定义克隆性和克隆进化模式复发。除了发现新的、反复突变的基因(例如,WAC,SMC3,DIS3,DDX41和DAXX)在AML中,我们还发现了AML复发期间的两种主要克隆进化模式:

(1)原发性肿瘤中的创始克隆获得突变和进化为复发克隆

(2)原始的亚克隆克隆在初始治疗中幸存下来,获得了额外的突变在复发时扩大。

在所有病例中,化疗都未能根除创始克隆。复发特异性与复发特异性的比较所有八例病例的原发性肿瘤突变均显示增加在转位中,可能是由于细胞毒性引起的 DNA 损伤化学疗法。这些数据表明,AML复发与新突变和克隆进化的增加有关,这部分是由患者的化疗决定的接受以建立和维持缓解。

生信分析流程

  • 相关数据集选择:

    Patients:8例AML患者的原发和复发配对样本

  • 测序:WGS和Capture array design

  • 基因集选择:突变基因结果列表

  • 生信分析方法:

我们从文章的分析流程中提取所有的分析内容,整理出来就 7个分析条目,每个条目都包括分析的内容,这些分析构成了整个文章,本文属于生信分析的文章,下面我们就看看哪些分析可以利用桓峰基因公众号的教程来实现,点击分析条码就会跳转到对应公众号的教程,跟着教程做,您也能发轻松发高分,如下:

1. 使用 BWA 0.5.5  将 PE 比对到 NCBI build36;

2. 使用 SomaticSniper and SAMtools indel caller 检查体细胞突变位点;

3. 使用 BreakDancer 识别 SV;

4. 使用 GATK and Pindel  detecting indels

5. 使用 VarScan 2 识别胚系突变;

6. 原发肿瘤和复发对亚克隆的核密度估计(sciClone)

7. 克隆进化鱼图(fishplot)

研究结果

1. 基于NGS在8个原发和复发配对的急性髓系白血病患者中捕获体细胞突变位点

a,在八例病例中检测到的1-3级突变位点(不是包括易位)。显示的所有突变均使用捕获进行验证其次是深度测序。共享突变为灰色,原发性肿瘤特异性突变为蓝色,复发特异性突变为红色。合计每个病例的 1-3 级突变数显示在浅灰色上方矩形。

b,来自验证突变的突变等位基因频率分布原发性肿瘤的 1-3 级和病例 UPN 933124的复发(左)。突变体获得五个原发肿瘤特异性突变的等位基因频率来自 454 深度读取计数实验。鉴定出4个突变簇在原发性肿瘤中,复发时发现一个。五个低水平突变在原发性肿瘤和复发中(包括 4 例已知拷贝数量变量区域)被排除在聚类分析之外。来自 AML 中反复突变的基因的非同义突变是显示。突变的突变等位基因频率的变化来自五个突变簇状物显示在原发肿瘤和复发之间(右)。

c,在来自另外七名 AML 患者的原发性肿瘤和复发样本。这原发肿瘤簇与复发样本之间的关系由连接它们的线表示。

89b9726aae3e0df273617e3b80cebf35.png

2. 在八个原发性肿瘤-复发对中鉴定出的编码突变

faba8e6d4687cd373c6b707149f63f09.png

3. 从初级克隆进化的图形表示UPN 933124中肿瘤复发以及肿瘤演变模式

在八个原发性肿瘤和复发配对中观察到。

a,创始克隆UPN 933124的原发肿瘤包含DNMT3A、NPM1、PTPRT、SMC3 和 FLT3 在 AML 和 AML 中反复出现可能与发病机制有关;创始克隆中的一个子克隆

进化为通过获得额外的突变,包括 ETV6 和 MYO18B 的复发性突变,以及WNK1-WAC融合基因。

b, 示例AML肿瘤演变的两种主要模式。模型 1 显示占主导地位原发性肿瘤中的克隆通过获得复发特异性突变进化为复发克隆;这种模式在三个原发性肿瘤和复发对(UPN 804168、UPN 573988 和 UPN 400220)。模型 2 显示携带绝大多数原发性肿瘤突变的小克隆复发时存活并扩大。在五个主要肿瘤和复发对(UPN 933124、UPN 452198、UPN 758168、UPN426980 和 UPN 869586)。

169a07ca5e7a2071149edf18d19b1454.png

4. 原发性肿瘤和复发样本。

a,原发肿瘤的比例和复发特异性每个转换和转化类别中的突变。

b, 转换原发肿瘤的频率和复发特异性突变来自八个AML 肿瘤和复发对。456 个复发特异性突变和 3,590 个8例原发性肿瘤突变用于评估统计学使用比例检验的显著性

3341ab6b610192e38abbd3e4b6608490.png

Reference:

1. Miller CA, White BS, Dees ND, et al. SciClone: inferring clonal architecture and tracking the spatial and temporal patterns of tumor evolution. PLoS Comput Biol. 2014;10(8):e1003665. Published 2014 Aug 7.

2. Dang, H. X., White, B. S., Foltz, S. M., Miller, C. A., Luo, J., Fields, R. C., & Maher, C. A. (2017). ClonEvol: clonal ordering and visualization in cancer sequencing. Annals of Oncology, 28(12), 3076-3082.

3. Miller CA, McMichael J, Dang HX, Maher CA, Ding L, Ley TJ, Mardis ER, Wilson RK. Visualizing tumor evolution with the fishplot package for R. BMC Genomics. 2016 Nov 7;17(1):880. 

4. Ding, L., Ley, T., Larson, D. et al. Clonal evolution in relapsed acute myeloid leukaemia revealed by whole-genome sequencing. Nature 481, 506–510 (2012).

这篇文章利用桓峰基因公众号上的教程完全可以实现这些图表的绘制,过程需要自己细致的探究,同时也需要有一定的生信分析和结果解读的能力,有需要此类分析的老师可以联系桓峰基因!

桓峰基因,铸造成功的您!

未来桓峰基因公众号将不间断的推出单细胞系列生信分析教程,

敬请期待!!

桓峰基因官网正式上线,请大家多多关注,还有很多不足之处,大家多多指正!

http://www.kyohogene.com/

桓峰基因和投必得合作,文章润色优惠85折,需要文章润色的老师可以直接到网站输入领取桓峰基因专属优惠券码:KYOHOGENE,然后上传,付款时选择桓峰基因优惠券即可享受85折优惠哦!https://www.topeditsci.com/

有想进生信交流群的老师可以扫最后一个二维码加微信

备注“单位+姓名+目的”!

b429e1d158bd6b2a180d8d595dea1d6a.png

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

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

相关文章

【雕爷学编程】Arduino动手做(93)--- 0.96寸OLED液晶屏模块11

37款传感器与执行器的提法,在网络上广泛流传,其实Arduino能够兼容的传感器模块肯定是不止这37种的。鉴于本人手头积累了一些传感器和执行器模块,依照实践出真知(一定要动手做)的理念,以学习和交流为目的&am…

UI设计工具都有哪些好用的推荐?

对于UI设计的初学者来说,掌握一个实用且易于使用的界面UI软件是非常重要的。今天,我整理了四个易于使用的界面UI软件。让我们看看。 即时设计 即时设计是一款免费的在线 UI 设计工具,无系统限制,浏览器打开即可使用,…

gin框架内容(二)

上一篇过于gin的内容 https://mp.csdn.net/mp_blog/creation/editor/131953861 CSDNhttps://mp.csdn.net/mp_blog/creation/editor/131953861 一、路由组 为了管理具有相同前缀的URL, 将拥有URL共同前缀的路由划分为一组 为了代码的阅读性,使用{}包裹相同组的路由…

SDN系统方法 | 7. 叶棘网络

随着互联网和数据中心流量的爆炸式增长,SDN已经逐步取代静态路由交换设备成为构建网络的主流方式,本系列是免费电子书《Software-Defined Networks: A Systems Approach》的中文版,完整介绍了SDN的概念、原理、架构和实现方式。原文: Softwar…

【WIN系统】创建任务计划程序:因为空格问题而无法创建的解决办法

在计算机管理中创建任务计划程序如下图: 【问题】 1. 点击“创建任务”或“创建基本任务” 2. 在操作一栏中,我们可以发现,之前已经存在的任务中,“程序或脚本”不带引号,而且允许存在空格。但自己输入时,…

Cesium 实战 - Blender调整模型组件原点,实现直升机尾翼旋转

Cesium 实战 - Blender调整模型组件原点,实现直升机尾翼旋转 1.模型原点问题2.导入模型(zhisheng.glb)3.导出模型4. 通过 czml 调试代码 某个项目需求,在操作直升机模型的时候,希望直升机机翼和尾翼旋转起来。 机翼旋…

多项目进度把控:选择合适的项目管理系统助力无忧

我相信很多企业都在使用项目管理系统,但你有没有想过为什么要使用这些工具?想象一下,如果我们使用传统的而不是这些项目管理系统,excel或work一些简单的在线文档记录能满足我们的需求吗?当需求发生变化时,这…

前端学习——Vue (Day5)

自定义指令 <template><div><h1>自定义指令</h1><input v-focus ref"inp" type"text" /></div> </template><script> export default {// mounted(){// this.$ref.inp.focus()// }// 2. 局部注册指令di…

桥梁安全监测系统中数据采集上传用 什么?

背景 2023年7月6日凌晨时分&#xff0c;G5012恩广高速达万段230公里加80米处6号大桥部分桥面发生垮塌&#xff0c;导致造成2车受损后自燃&#xff0c;3人受轻伤。目前&#xff0c;四川省公安厅交通警察总队高速公路五支队十四大队民警已对现场进行双向管制。 作为世界第一桥梁…

亚马逊云科技纽约峰会官宣生成式AI产品策略升级

作为云计算领域的领导者和创新者&#xff0c;生成式AI&#xff08;Generative AI&#xff09;一直是亚马逊云科技持续关注和投入的主要方向。在今年4月&#xff0c;亚马逊云科技发布了以Amazon Bedrock为代表的生成式AI产品全家桶&#xff0c;正式入局该赛道并宣布了亚马逊云科…

进程信息查看脚本

一、shell 脚本 该脚本可以根据进程名和进程pid查询进程信息。 输出的进程信息包括&#xff1a;进程 PID、进程命令、进程所属用户、CPU占用率、内存占用率等 查询可选择模式&#xff1a;仅一次还是不限次。 #! /bin/bashprintProcessInfo(){P$1echo "------------------…

15 文本编辑器vim

15.1 建立文件命令 如果file.txt就是修改这个文件&#xff0c;如果不存在就是新建一个文件。 vim file.txt 使用vim建完文件后&#xff0c;会自动进入文件中。 15.2 切换模式 底部要是显示插入&#xff0c;是编辑模式&#xff1b; 按esc&#xff0c;底部要是空白的&#xff0…

NetCore 使用 Swashbuckle 搭建 SwaggerHub

什么是SwaggerHub? Hub 谓之 中心, 所以 SwaggerHub即swagger中心. 什么时候需要它? 通常, 公司都拥有多个服务, 例如商品服务, 订单服务, 用户服务, 等等, 每个服务都有自己的environment, endpoint, swagger schema. 然而这些信息都分散在各处, 如果能集中在一个地方展示…

自主AI代理:未来的生产力引擎

摘要 文章介绍了自主AI代理的概念&#xff0c;AI代理由AI驱动&#xff0c;能够自我创建、优先处理和完成任务。自主AI代理可以执行任何数量的任务&#xff0c;包括内容创建、个人助手、个人财务管理、研究和数据分析等。文章强调了知识、记忆和学习在构建成功的自主AI代理中的重…

任务与项目的巧妙运用:项目管理软件的有效实践方法

如果您不熟悉项目管理或以前依赖任务管理系统来管理您的项目&#xff0c;那么任务和项目之间的区别可能会令人困惑。任务和项目是项目管理软件中的主要构建块&#xff0c;使您能够跟踪和组织您的工作。 任务是需要在项目中完成的单个工作单元。项目是需要一起完成以实现单个结果…

Linux数据处理三剑客

目录&#xff1a; Linux 三剑客之 greplinux三剑客之awklinux三剑客之sedlinux三剑客与管道使用【实战】三剑客实战之nginx日志分析实战【实战】三剑客实战之性能、网络统计实战linux进阶命令linux环境配置Linux与Bash编程实战 1.Linux 三剑客之 grep 内容检索:&#xff08;…

高电压放大器ATA-2021B在无损检测领域中的应用

超声无损检测&#xff08;UltrasonicNondestructiveTesting&#xff0c;简称UT&#xff09;是一种常用的材料及构件内部缺陷检测技术。它利用超声波在材料中传播的特性&#xff0c;通过接收回波信号来检测材料内部的缺陷情况。超声无损检测具有精度高、快速、非破坏性等优点&am…

k8s Webhook 使用java springboot实现webhook 学习总结

k8s Webhook 使用java springboot实现webhook 学习总结 大纲 基础概念准入控制器&#xff08;Admission Controllers&#xff09;ValidatingWebhookConfiguration 与 MutatingWebhookConfiguration准入检查&#xff08;AdmissionReview&#xff09;使用Springboot实现k8s-Web…

【录用案例】2区SCI,仅14天见刊!

近日新增两本2区SCI录用&#xff0c;期刊表现质量优&#xff0c;均是隶属于世界前列的出版社&#xff0c;好刊版面紧俏&#xff0c;可放心投稿&#xff01;录用案例如下&#xff0c;可重点参考&#xff1a; 计算机科学类SCI&EI 【影响因子】IF&#xff08;2022&#xff0…

Microsoft todo 数据导出

文章目录 官方说明&#xff1a; https://support.microsoft.com/zh-cn/office/导出您的-microsoft-待办事项帐户-d286b243-affb-4db4-addc-162e16588943 由于 微软待办 会自动与 Outlook 中的任务同步&#xff0c;因此您可以从 Outlook 中导出所有列表和任务。 若要导出列表和…