- 综合排序
- 最热优先
- 最新优先
- 来自专栏作图丫
deconstructSigs 突变 Signature 分析
面对大量的SNV突变数据你是否还觉得无从下手，不知道怎么分析合适？今天给大家介绍一个R包-deconstructSigs。这款R包是基于大样本量预测的signature解析突变特征。 sample.mut.ref) Sample chr pos ref alt 1 1 chr1 905907 A T 2 1 chr1 1192480 C A 3 sample.id # 突变文件中的样品列名 chr # 突变文件中的染色体列名 pos # 突变文件中的突变位置列名 ref # 突变文件中的参考基因组碱基列名 alt # 突变文件中的突变碱基列名 TRUE, tri.counts.method = 'default') whichSignatures参数释义 tumor.ref # 上一步生成突变文件 ，数据为数据框或文本，横行是样本，纵行是突变碱基上下文序列 sample.id # 样品名称，tumor.ref文件的行名 signatures.ref # 预测的已知signatures参考文件
1.9K30编辑于 2022-03-29 - 来自专栏生信技能树
突变位点生存分析
一个简单突变位点做生存分析居然拖了一两个月才有人提交笔记！ 加油哈，广大粉丝们 1 主要流程 1.本次选用BRCA的maf数据和临床数据，主要使用其中的varscan数据 2.使用R包maftools读取maf文件，并可视化top10突变基因 3.选取两基因对BRCA 临床样本进行分组 所选取两基因都未发生突变的样本为一组 剩余样本为一组 4.使用logrank进行生存分析 2 代码及结果图 1.读取maf文件并对数据进行可视化 options(download.file.method ，进行生存分析 #整理生存分析的输入数据---- #1.由随访时间和死亡时间计算生存时间(月) is.empty.chr = function(x){ ifelse(stringr::str_length (x)==0,T,F) } is.empty.chr(meta[2,3]) meta[,3][is.empty.chr(meta[,3])]=0 meta[,4][is.empty.chr(meta[,
2.1K31发布于 2020-06-19 - 来自专栏用户7627119的专栏
肿瘤新抗原突变负荷分析
用R包limma对定义的新抗原负荷亚群进行差异突变分析，FDR < 0.05。根据hg19人参考基因组，分析了30个突变特征，并比较了所识别的亚组之间的突变特征。 新抗原负荷亚组的不同功能通路 对三个亚组的差异表达分析识别出5个显著差异表达基因（CXCL9、CXCL13、IGLL5、AGR3和TFF3）。 接下来评估了30个突变特征，以更好地理解复杂的突变过程。得到5个差异显著的特征，分别是特征1、特征3、特征6、特征13和特征30（图2B）。 NAL-M在特征3、13和30中富集，特征3与乳腺癌、胰腺癌和卵巢癌的生殖系和体细胞BRCA1和BRCA2突变密切相关。特征13表示胞苷脱氨酶的AID/APOBEC家族活性。 图3.新抗原负荷亚组的生存分析及风险因素分析 表1.根据临床参数进行样本比较 ?
3.7K41发布于 2020-08-05 - 来自专栏生信菜鸟团
使用 sigminer 进行突变模式分析
突变模式分析（Mutual Signature Analysis）已经逐步成为变异检测后一个通用分析，本文简单介绍如何使用sigminer进行突变模式分析，以解决2大分析任务： 从头发现签名 已知一些参考 如果你会使用maftools读入突变数据，那么就会使用sigminer读入突变数据，支持 data.frame 和MAF文件。 这里我们使用maftools内置数据集，maftools其实本身也可以做signature分析（但不是它的核心）。 上面展示了最核心的分析和可视化功能，sigminer还支持很多功能，我就不再出现了。 当然，读者完全可以基于上面的分析的结果值进行各种个性化分析。
2.5K21发布于 2020-06-10 - 来自专栏作图丫
MutationalPatterns--进行肿瘤突变分析！
背景介绍 突变过程在基因组 DNA 中留下特征足迹。 install.packages("BiocManager") BiocManager::install("MutationalPatterns") library(MutationalPatterns) 分析流程 突变谱显示了每个突变类型在碱基替换目录中的相对贡献。 ncol = 3, widths = c(3, 3, 1.75)) p4 <- plot_spectrum(type_occurrences, by = tissue, CT = TRUE, legend plot_spectrum_region(type_occurrences_region, mode = "relative_sample") 小编总结 MutationalPatterns的新版本在多种突变分析中均有强大的功能
1.7K21编辑于 2022-03-29 - 来自专栏生物信息学、python、R、linux
SCmut||分析单细胞数据突变
单细胞RNA测序(scRNA-seq)和DNA测序(scDNA-seq)都可以应用于细胞水平基因组分析。对于突变分析，scDNA-seq似乎更常见。 但目前DNA测序中检测突变的方法多种多样，尚不清楚这些方法是否可以用于scRNA-seq数据。 在scRNA-seq乳腺癌数据集中，SCmut可以识别许多高度可信的细胞水平突变，这些突变在许多细胞中都反复出现，并且在不同样品中保持一致。 简而言之，该方法首先从肿瘤和匹配的种系组织的大细胞DNA测序（bcDNA-seq）中收集体细胞突变。 在（i）中，发现的细胞水平突变在肿瘤细胞和非肿瘤细胞之间被很好地分开，在（ii）中，突变被同时在两个独立的数据集中发现。
1.4K10发布于 2020-09-03 单细胞空间突变信息分析导论
今天周六，昨天相亲的91年大姐要求结婚必须再买一个婚房，自然也就没有然后了，10点了还躺着， 我们来简单汇总一下单细胞空间分析突变的内容吧。 10X空间转录组数据中研究并可视化体细胞突变单细胞、空间、外显子解析TP53突变重构肺腺癌细胞图谱单细胞、空间、外显子分析方法更新单细胞、空间、外显子多组学分析探讨多组学（单细胞、空间转录+蛋白、外显子 、甲基化）揭示神经母细胞瘤异质性图谱2024空间转录组分析方向和文献分享单细胞、空间、外显子多组学分析揭示了早期肺腺癌（LUAD）发展过程中上皮细胞的状态和转录组特征系统整理10X单细胞空间数据中可检测到的有害突变位点 ，3’数据和空转数据只能检测到基因末尾98bp的范围，5‘数据能检测到基因开始的98bp范围。 X 基因位置突变的矩阵（以单细胞3'数据为例）,row是基因 + 染色体 + 位置 + ref +alt，column是barcode。
45820编辑于 2024-04-13- 来自专栏作图丫
几项泛癌driver突变分析汇总
通过将突变异质性纳入分析，MutSigCV能够消除大多数明显的人工发现（造成的偏倚），并能够识别真正与癌症相关的基因。 分析了特定癌症类型患者的异质性。 最近发现，在表达DNA胞苷脱氨酶APOBEC3B酶的乳腺癌中，APOBEC3B占突变负荷的一半。该研究讨论了APOBEC3B是否对多种肿瘤类型的突变负责。 分析了19种不同癌症类型的基因表达数据和突变模式、分布和负荷，其中有超过4800个外显体和1000000个体细胞突变。 结合该工作提供的全面的表达的数据分析(图9)、CG碱基对突变频率(图10)、局部胞嘧啶突变特征(图11)、整体突变负荷(图12)和kataegis(图12 c)，所有可用的数据得出结论APOBEC3B是一个多种人类癌症的突变的主要来源 上面这些文章都是从泛癌角度多肿瘤突变进行分析描述，建议大家可以去浏览一下感兴趣的文章哦，里面的绘图和分析都很值得我们借鉴！ 引用： [1].
1.6K11编辑于 2022-03-29 - 来自专栏R语言交流中心
R语言中的突变signature分析工具
新的一年开始了，今天给大家介绍一款用于发现、预测和探索突变特征的综合分析工具包musicatk。此包主要基于COSMIC突变数据中的最新数据进行肿瘤突变模式的探索。 参考基因组的选择 g <- select_genome("hg38") 3. 项目操作 ##创建项目。 突变基序包括单碱基替换(SBS)、双碱基替换(DBS)、插入(INS)和删除(DEL)。 和加速版本非负矩阵分解（nmf） result <- discover_signatures(musica = musica, table_name = "SBS96", num_signatures = 3, 绘制单个样本的signature samples <- sample_names(musica) plot_sample_counts(musica, sample_names = samples[c(3,4,5
2K31编辑于 2022-03-14 - 来自专栏作图丫
突变signature分析你不能错过的R包！
今天给大家带来的是signature分析的R包“YAPSA”，让大家在分析signature的时候多一个选择，增加绘图展示的多样性，最重要的是让你的老板知道你有多优秀。 大家知道前文分享的“deconstructSigs”（点击看往期详情）也是突变signature的绘图软件，好奇的小伙伴就有疑问了，这俩要是一样就没必要再学了吧？ 是利用非负矩阵分解（Nonnegative Matrix Factor）的方法预测；另外，“YAPSA”方法展示可通过阈值和颜色进行调整，较为灵活；再者，“YAPSA”运行更快和所需的资源更少（个人感觉，在突变不是很多的情况下 set(echo=TRUE) opts_chunk$set(fig.show='asis') library(BSgenome.Hsapiens.UCSC.hg19)#注意基因组版本 2、准备数据 3、 <- 3 data1 <- create_mutation_catalogue_from_df( data, this_seqnames.field = "CHROM
1.2K10编辑于 2022-03-29 - 来自专栏生信菜鸟团
新工具：基因组分析和突变检测
数据：文件基于千人基因组 1kGP的数据，通过分析 WGS 种系变异来展示 DRAGEN 的功能和结果。 DRAGEN* *的WGS种系分析：简单来说就是分析速度快、全面的变异检测、可扩展性、简化的工作流程 示例代码：在 illumina 的 ICA 平台上，使用DRAGEN v.4.2流程进行比对和变异检测 ：约检测到 492 万个小突变，其中包括 3,956,307 个 SNV， 960,908 个 插入缺失， 13,886 个SV （≥50 bp），1,156 个 CNV（图2a）。 DRAGEN 在其他样本或数据集的突变检测结果如图3a所示，DRAGEN 在全基因组数据的SNV和Indel突变检测中，召回率在 99.76% 和 99.87% 之间，精确率在 99.90% 和 99.93% 与其他工具或流程相比较，DRAGEN 的性能也更好（图3cde）。
66410编辑于 2024-12-05 - 来自专栏生信修炼手册
使用MutationalPatterns进行肿瘤突变频谱分析
MutationalPatterns是一个bioconductor上的R包，可以用于肿瘤突变频谱的分析。 肿瘤突变频谱针对点突变进行定义，A,T,C,G四种碱基两两突变，共有4X3=12种排列，考虑到正负链碱基配对原则，正链上的A->C突变，对应负链上为T->G, 所以进一步转换成了一个组合的问题，所以某个位点的突变可以划分为以下 ，即上下游各取一个碱基再加上突变位点的碱基，组成了3个碱基的motif, 可以有4X4X6=96种模式，每种模式的频率分布就是突变频谱。 3. 通过这个R包，可以轻松实现突变频谱的常见分析内容。
3.6K20发布于 2019-12-19 - 来自专栏R语言交流中心
R语言实现突变信号（Mutational Signatures）分析
每一个突变过程都会留下一个不同的基因组标记，也就称为突变信号。 ? 实例： sample_names <- c ( "colon1","colon2", "colon3","intestine1","intestine2", "intestine3","liver1", "liver2", "liver3" ) # We assemble a list of files we want toload. 3. RDS格式数据读取源自base包readRDS函数，其实就是R语言自己的数据格式，此格式在可以输出到文件。 ? 3. plot_compare_profiles 多组数据之间的对比情况可视化 实例： plot_compare_profiles(mut_mat[,1],
9.9K55发布于 2019-07-31 - 来自专栏分析工具
maftools癌症体细胞变异(突变)分析工具学习
Maftools 是一个专门用于分析和可视化突变数据的 R 包。 Maftools 的主要功能包括：突变数据解析：可以方便地加载 MAF 文件，并将其解析为 R 中的对象，便于进一步分析。 突变数据可视化：提供了多种强大的可视化工具，如瀑布图（Waterfall plot）、突变频率图、共突变和互斥突变分析、突变频谱分析等，帮助研究人员更直观地理解突变数据的特征。 临床数据整合：可以将临床数据与突变数据整合，分析不同临床亚群的突变特征和生存分析等。统计分析：包括突变负荷分析、通路富集分析、TCGA数据整合分析等功能。 首先是根据突变情况联合各种临床参数；其次现已经进入了免疫治疗时代，可以将其与免疫治疗的数据联合分析；或者如果重点要研究某一个分子的突变情况，研究者可以先来看看这里的突变数据~ 总之，用途还是非常多的！
1.2K10编辑于 2024-09-02 课程丰富----单细胞联合突变的信息分析
作者，Evil Genius多组学是分析的潮流，也有助于我们深刻的认识生物学问题，甚至于找到解决问题的方法，单一组学都有局限性，一叶障目并不可取。当然了我希望大家能抓住空窗期。 汇总一下突变的用药指南卵巢癌的运用卵巢癌(HGSOC)的主要特征是拷贝数改变和基因组重排形式的深刻结构变异，这是在几乎普遍存在的TP53突变的遗传背景下产生的。 骨髓瘤的运用，研究肿瘤异质性肺腺癌的运用全外显子组测序(WES)、单细胞RNA测序和单细胞TCR测序基因组 + 空间组学H3-mutant gliomas生活很好，有你更好
27220编辑于 2025-04-03- 来自专栏用户7627119的专栏
maftools包分析突变数据，绘制瀑布图
☞ 如何从TCGA数据库下载体细胞突变数据（somatic mutation） 今天我们来讲讲，怎么用R的maftools包来分析MAF格式的突变数据，并用瀑布图来展示结果。 maftools这个包的主要分为两部分功能，分析和可视化。下图列出了，这个包中相应的函数的名字。 我们先用maftools包自带的数据，给大家讲解这个包的使用方法。 = laml,addStat = 'median') dev.off() 我们会得到下面的一张突变的汇总图，包括各种突变分类统计图，突变类型统计图，碱基改变统计图，每个样本包含突变数统计图，样本中各种突变分类的箱型图 ，突变最多的10个基因所包含的突变类型，以及样本占比情况。 接下来我们就可以来绘制瀑布图了，我们可以通过top来控制展示多少个突变最多的基因，这里展示20个突变最多的基因。
1.9K40编辑于 2022-09-21 - 来自专栏作图丫
几项泛癌基于突变的网络模块分析汇总
最后，为了使聚类的分配具有鲁棒性，使用一致性聚类分析1000次重复抽样（图1d）。 也分析了癌症网络中反复发生的突变和富集的信号通路，以及在已知的癌症基因以及新的候选癌症基因和通路中发现了影响信号传导的突变。 新的细胞周期driver候选包括ATR，一种磷酸化p53和其他蛋白的激酶，以及PIK3-AKT信号通路中的PIK3CG和PIK3CB，似乎补充了PIK3CA的致瘤作用。 小编总结 对于样本网络的研究，是生信分析中非常重要的一部分，而突变产生的影响可以由网络进行传播破坏正常的分子功能，了解这种因果机制在研究中非常重要。 Sci Rep. 2013;3:2651.
57521编辑于 2022-03-29 - 来自专栏生信技能树
带临床信息的肿瘤突变maf文件分析维度更多
上次我们分享了 指定病人的指定基因的突变全景瀑布图，主要是讲解了maftools这个包的个性化操作，这个教程里面我们仅仅是使用了 TCGA-Clinical Data Resource (CDR) Outcome 而带临床信息的肿瘤突变maf文件分析维度更多。 每个癌症都去找各自的肿瘤突变maf文件很麻烦，所以我们才会选择 PanCanAtlas Publications Scalable Open Science Approach for Mutation 首先病人按照临床指标分组后看突变差异 这里我们仅仅是按照stage举例，首先把stage给格式化： table(phe$ajcc_pathologic_tumor_stage) phe$new_stage 下期预告：有生存统计信息就可以做生存分析，根据突变位点， 可以在各个癌症内部批量进行生存分析检验，挖掘有意义的基因或者基因列表。
1.7K10发布于 2021-10-12 - 来自专栏作图丫
Cancer3D*v2 ：癌症突变的三维模式的互作分析
Res IF: 11.147 Published: 2015；v2.0 2019 Cancer3D数据库是最早将错义突变等基因组数据与包含这些突变的蛋白质的结构域和三维结构信息相结合的数据库之一，可以分析在患者临床数据中发现的和与患者临床数据相关的蛋白质结构中的癌症错义突变 你可以在这个数据库中为特定的子群体找到新的候选驱动区域，特别是在基因水平上进行类似的分析时找不到的。还能够对不同癌症类型和分期、不同性别和年龄、病人的种族中的突变分布进行可视化。 使用BLAST来匹配数据库中所有基因的三维结构序列，并将突变的位置映射到蛋白质结构上。 3. 以BRAF（致癌基因）为例分析蛋白突变，下图是主注释页面，缩放查看整体或局部详细信息。 ①突变直方图。V600突变最频繁的点。 ②蓝色方框表示Pfam域，紫色方框表示PDB的位置。 (3) Tutorial 数据库提供的使用说明。 小编总结 Cancer3D数据库允许用户分析癌症突变的三维结构及其临床相关信息背景下的突变分布模式。
64220编辑于 2022-03-29 课后补充--关于单细胞突变联合inferCNV的分析
作者，Evil Genius我们的课程已经上了4节了，包括单细胞突变、轨迹、CNV和通讯，这一篇做一个联合分析的补充。 1、根据肿瘤和正常组织的特征基因表达来计算肿瘤评分2、通过使用scRNA-seq数据估计基因组拷贝数变化的算法infercnv计算细胞的拷贝数变异(CNV)分数3、单细胞数据中检测肿瘤特异性突变其中CNV 和突变都已经讲完了，大家要活学活用。
26400编辑于 2024-06-30
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号