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

推荐订阅源

V
V2EX
小众软件
小众软件
GbyAI
GbyAI
B
Blog RSS Feed
月光博客
月光博客
A
About on SuperTechFans
Microsoft Security Blog
Microsoft Security Blog
Blog — PlanetScale
Blog — PlanetScale
U
Unit 42
云风的 BLOG
云风的 BLOG
P
Proofpoint News Feed
Recent Announcements
Recent Announcements
Cyber Security Advisories - MS-ISAC
Cyber Security Advisories - MS-ISAC
Vercel News
Vercel News
F
Fortinet All Blogs
Apple Machine Learning Research
Apple Machine Learning Research
freeCodeCamp Programming Tutorials: Python, JavaScript, Git & More
雷峰网
雷峰网
Google DeepMind News
Google DeepMind News
MongoDB | Blog
MongoDB | Blog
Engineering at Meta
Engineering at Meta
aimingoo的专栏
aimingoo的专栏
博客园_首页
G
Google Developers Blog

博客园 - 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 

参考: