类器官外显子测序数据标准分析
以 Clinical & Translational Immunology 2021; 文章为例,标题是:《Patient-derived organoids of bladder cancer recapitulate antigen expression profiles and serve as a personal evaluation model for CAR-T cells in vitro》,是中国深圳的中科院先进院的团队的研究成果。
跟普通的肿瘤外显子测序数据的分析步骤并没有较大差异,如下所示:
普通的肿瘤外显子测序数据的分析
虽然这个研究仅仅是产生了3个类器官数据,但是需要有9个样品进行外显子测序哦,其中每个病人的正常对照仅仅是数据分析过程使用,并不需要在后续可视化进行呈现。对类器官数据的可视化主要是比较样品培养前的肿瘤组织和培养后的类器官组织的突变结果的一致性, 首先看SNV;
SNV的一致性
可以看到, 绝大部分的SNV在每个样品培养前的肿瘤组织和培养后的类器官组织的突变结果具备非常好的一致性。包括上图显示的各种突变形式:
- onsense mutations (e.g. ARID1A and ERBB2)
- missense mutations
- splice site mutations (e.g. EP300 and FAT1)
- multiple hit mutations
- frame shift indels/dels (e.g. TP53, RB1,CREBBP and CDKN1A)
其实就是maftools的标准图表啦:
maftools的标准图表
还有CNV的对比:
CNV的一致性
样本量多少都是类似的分析
总之就是要说明自己的 organoids recapitulate the genetic alterations in the parental tumors. 既然是外显子测序,那就是从肿瘤外显子的SNV和CNV两个分析结果的角度来说明。
再比如:2021年12月17日,北京大学第一医院泌尿外科联合深圳市第二人民医院在Advanced Science期刊 (IF=16.8) 在线发表了标题为《Patient-derived Upper Tract Urothelial Carcinoma Organoids as a Platform for Drug Screening》的研究成果,也是这样的展现:
首先看:A) Somatic genomic landscape of 15 UTUC organoid lines (_O) and the corresponding parental tumors (_T). 仍然是 maftools的标准图表 :
maftools的标准图表看SNV的一致性
另外,他们在展现CNV的一致性,应该是采用了发表在Gigascience . 2021 May的工具《MesKit: a tool kit for dissecting cancer evolution of multi-region tumor biopsies through somatic alterations》:
展现CNV的一致性
值得玩味的是文章里面并没有提到 这个MesKit工具,我在原文里面看了看:https://onlinelibrary.wiley.com/doi/full/10.1002/advs.202103999 ,其提到的肿瘤外显子数据分析流程里面的软件工具列表如下所示:
- Fastp (v0.12.6) was used to remove adaptors and filter low quality reads.[39]
- Single-nucleotide variant (SNVs) were called by GATK(v4.1.9).[40]
- Sequence reads were mapped against the human reference genome (hg38) using Burrows-Wheeler Alignment with maximal exact matches (BWA-MEM) v0.7.12.[41]
- Read mapping was followed by the marking of duplicates, merging of lanes, and realignment of indels using Picard (v2.23.8).
- Somatic SNVs and indels in the tumors and organoids were called with Mutect2 (default options) using the matched blood samples as the reference.
- Somatic SNVs with a VAF < 0.05 and supported by less than 3 reads were filtered out.
- Mutation effect predictions and annotations were performed by VEP (release 101).[42]
- CNVs were detected using TitanCNA (1.30.0).[43]
- The clone architecture and the VAF of mutated sites in organoids and the corresponding tumors were inferred by sciClone.[44]
相关文章
- Abp + MongoDb 改造默认的审计日志存储位置
- 《CLR Via C#》读书笔记:27.计算限制的异步操作
- 程序员如何开发高效的程序,这项技术值得了解
- 四大技术变革重塑企业数据库
- [Abp 源码分析]七、仓储与 Entity Framework Core
- 用 C# 在 LinqPad 建立 Linked Server 跨服务器数据库操作
- 使用 SharpZipLib 打包数据到 ZIP 文件
- 1行代码就能跑个量子计算!AWS年度巨献 | 狄拉克孙子点赞
- 过去五年20种涨跌势头最猛的技术技能
- 【NCTS峰会回顾】阿里巴巴潘家腾:阿里妈妈在线下测试域的智能化建设
- 不亚于谷歌量子霸权 微软的拓扑量子位即将实现
- 假如谷歌的“量子优越性”是一场革命,我们还应该知道什么?
- Kafka事务流由基础到实践
- 我花10个小时,写出了小白也能看懂的阿里数据中台分析
- 用Rust让numpy、scikit和pandas加速100倍!开源Weld技术揭秘
- 媒体近距离接触谷歌量子计算机:还要解决哪些问题
- 谷歌宣称首次实现量子优越性,IBM“不服”,中国同行咋看?
- 别逗了,你真以为分库分表就能无限扩容吗?
- 没有爬虫,就没有今天的搜索引擎
- Nature发文:量子计算研究应该全球互通,不应筑起壁垒