cnvkit被设计来处理同一个批次的多个肿瘤配对样本测序情况，首先对所有的normal数据进行bin处理拿到背景值，然后就这个背景值来处理所有的tumor测序数据计算拷贝数变异情况。该软件使用比较复杂，建议读一读官网教程。所有的命令都被包装到一个python脚本里面，使用该脚本调用一系列字命令，如下：官网教程Copy number calling pipelinebatchtarge...阅读全文... 阅 读 全 部 >
这一节里，我们来看一下独立资源的几个例子。一个很重要的资源是蛋白三维结构数据库 Protein DataBank （PDB）。我们这门课里对于蛋白三维结构的计算和预测讲的很少，但这其实也是非常重要的资源和研究方向。Protein DataBank现在主要是由Rutgers和UC Sandiego共同来管理。它目前共有9万多，将近10万的三维结构，包括8万多的X-r...阅读全文>>... 阅 读 全 部 >
首先说说GATK可以做什么。它主要用于从sequencing 数据中进行variant calling，包括SNP、INDEL。比如现在风行的exome sequencing找variant，一般通过BWA+GATK的pipeline进行数据分析。要run GATK，首先得了解它的网站(http://www.broadinstitute.org/gatk/)。一.准备工作...阅读全文>&g... 阅 读 全 部 >
在 6 月份，基因组学领域的权威期刊Genome Biology发表了一篇方法学评估类文章，在这篇文章中作者系统地分类和评估了目前主要的69种基于全基因组测序（whole genome sequencing，WGS）数据分析结构变异的算法/方法：Comprehensive evaluation of structural variation detection algorithms for ...... 阅 读 全 部 >
在特定的细胞或组织中，发挥作用的转录因子在调控网络中的地位也是不同的。来自胚胎干细胞和其他细胞模型的研究表明，一小组细胞类型特异性或谱系特异性转录因子（TFs）形成了相互连接的自动调节环，以控制特定细胞类型的转录程序。将某种细胞或者组织中的核心转录因子及其调控回路称之为core transcription regulatory circuitry, 简称CRC。核心TF及其相互连接的...阅读全文... 阅 读 全 部 >
本篇文章给大家分享的是有关怎么使用ROSE鉴定超级增强子，小编觉得挺实用的，因此分享给大家学习，希望大家阅读完这篇文章后可以有所收获，话不多说，跟着小编一起来看看吧。ROSE是最经典的超级增强子预测软件，由Richard A. Young大牛团队开发，源代码的网址如下http://younglab.wi.mit.edu/super_enhancer_code.html...阅读全文>>... 阅 读 全 部 >
超级增强子概念 超级增强子 (Super enhancer, SE) 是一类具有超强转录激活特性的顺式调控元件，2013年由美国学者Richard A. Young首次提出。与普通增强子 (Typical enhancer, TE) 相比，超级增强子区域跨度范围通常可达8~20Kb，远高于普通增强子...阅读全文&g... 阅 读 全 部 >
表达谱分析中，经常会使用到相关性分析，探索一组基因间的共表达特征。例如，这些基因间的表达是否存在较强的协同性，一个基因表达值的改变是否与另一个基因表达值改变显著相关，它们之间是共激活还是抑制关系等。对于相关性分析结果的可视化，通常有多种方法，例如相关性散点图、弦图、共表达网络图。本篇教程则主要带大家了解如何基于基因表达值的相关性，绘制弦图。如下弦图来自文献“CUL4B pro...阅读全文>... 阅 读 全 部 >
对于Hi-C互作分析而言，不同分辨率水平可以分析不同的结构特征，其中，分析A/B compartment结构需要的分辨率在100kb，分析TAD结构需要的分辨率在40kb，分析loop结构的分辨率在10kb。随着研究的深入，越来越多的研究者希望通过Hi-C技术去分析DNA loop结构。在染色体的折叠过程中，DNA loop环的形成让线性距离上相隔很远的位点有机会靠近并发生互作，这些loop的..... 阅 读 全 部 >