惯性聚合 高效追踪和阅读你感兴趣的博客、新闻、科技资讯
阅读原文 在惯性聚合中打开

推荐订阅源

A
About on SuperTechFans
有赞技术团队
有赞技术团队
人人都是产品经理
人人都是产品经理
月光博客
月光博客
美团技术团队
博客园 - 聂微东
阮一峰的网络日志
阮一峰的网络日志
WordPress大学
WordPress大学
让小产品的独立变现更简单 - ezindie.com
让小产品的独立变现更简单 - ezindie.com
博客园_首页
爱范儿
爱范儿
G
Google Developers Blog
aimingoo的专栏
aimingoo的专栏
T
The Blog of Author Tim Ferriss
MongoDB | Blog
MongoDB | Blog
小众软件
小众软件
OSCHINA 社区最新新闻
OSCHINA 社区最新新闻
IT之家
IT之家
I
InfoQ
B
Blog
H
Hackread – Cybersecurity News, Data Breaches, AI and More
大猫的无限游戏
大猫的无限游戏
T
Tailwind CSS Blog
F
Fortinet All Blogs

博客园 - Life·Intelligence

SSHFS + VS Code 挂载集群代码目录(macOS)| 集群vibe coding OpenClaw 多 Channel 实战总结(Windows 环境) OpenClaw Windows 安装与 Debug 最终版教程(适用于 MiniPC i3-N305 / 无 GPU) miniconda转miniforge | conda | license LightDock | 蛋白质-多肽对接 | peptide-protein docking 蛋白结构预测 | alphafold | colabfold | docking 全栈生信 | PyMol使用教程 细胞通讯推断 | CCI | CellChat | CellphoneDB | iTALK | NicheNet 极简 | GRN | SCENIC | pySCENIC | 安装使用最新版scenicplus Linux下载zenodo数据 共定位 | colocalization 分析 | 表观因子 findOverlappingPeaks | peak取交集操作 根据基因名批量查询下载PDB蛋白结构数据库 R小技巧汇总 Signac处理bulk ATAC-seq数据 Differential motif enrichment | CentriMo | meme 亚马逊云 | AWS S3 | 基本操作 ATAC-seq | TOBIAS | footprint分析 TCGA+GTEx基因表达数据合并 | 多癌种表达分析 药物筛选 | drug screen
ChromHMM教程
Life·Intelligence · 2025-03-28 · via 博客园 - Life·Intelligence

2025年03月27日

先看一篇高分文章:Chromatin state dynamics confers specific therapeutic strategies in enhancer subtypes of colorectal cancer

With the use of the ChromHMM algorithm,14 we trained a 10-state, 12-state and 15-state model (for
samples having all the six marks; online supplemental table S1)
consisting of combinatorial patterns of specific histone modifications known as chromatin states, which are linked with the
transcriptional output of the associated genomic loci and regulate the recruitment of various components of the transcriptional
machinery. We chose a 10-state model for in-depth investigation
because it represented biologically interpretable combinatorial
patterns that were not repetitive or ambiguous (figure 1A and
online supplemental figure S1A). We annotated these chromatin states based on the content of each state and its average
relative distribution around known genomic features (online
supplemental figure S1B).

因为一个基因组locus可以有很多歌histone state,某些histone mark之间具有协同或拮抗作用,所以这里相当于做个一个(maybe线性)变换,输入是所有histone mark,输出是统一的chromatin state。

如果有超过5个以上的histone marks,那这个就属于必做的分析。

代码

ChromHMM.sh BinarizeBam -gzip -b 200 -f 2 -paired chrom.sizes.hg38 3.bam design_sheet.txt binarization
# 容易爆内存,我猜一个线程会占用很多内存
default_jvm_mem_opts="-Xms20m -Xmx30g"
ChromHMM.sh LearnModel -gzip -d 0.001 -color 0,0,255 -p 4 -i chrhmm binarization learnmodel 10 hg38

参考项目:/home/zz950/projects/BAF_SOX9/diffbind/chromhmm 

参考: