3D基因组工具（HiC可视化）trackc--bioinfomatics tools 35

让学习成为一种生活方式2024-12-04 13:08
01 3D genome data analysis guides

茶树三维基因组-文献精读19

复制代码 
https://trackc.readthedocs.io/en/latest/install.html    #官网
https://github.com/seqyuan/trackc           #官网

https://trackc.readthedocs.io/en/latest/analysis_guide/index.html   #HiC可视化案例
GridSpec

在基因组区域中可视化多组学数据通常需要结合多个轨迹（tracks）。为了创建适合发表的高质量图像，精细调整每个轨迹的比例非常重要。快速调整轨迹比例以及轻松添加或移除轨迹的能力是一个重要问题。

trackc 提供了两种方法来绘制布局：

  1. trackc.tenon

  2. trackc.make_spec

这些函数为用户提供了一种快速灵活的方法，用于调整每个轨迹的比例，同时方便地在多组学可视化布局中添加或移除轨迹。

1. trackc.tenon

  • trackc.tenon 返回一个对象。例如：

    复制代码 
    ten = trackc.tenon(width=5, height=1)

    在此示例中，ten 是一个表示虚拟图形的对象，其宽度为 5，高度为 1。用户可以通过以下方法向该图形添加轨迹：

    复制代码 
    ten.add(pos='bottom', height=4, hspace=0.1)

    • height 参数 表示新添加轨迹的相对高度。例如，如果虚拟图形的高度为 1，而新轨迹的高度为 4，则子图的实际高度将为 1 * 4。

    • 使用 ten.add 方法 类似于搭建积木，用户可以通过 pos 参数 选择将新轨迹添加到顶部（top）或底部（bottom）。

    • hspace=0.1 参数 控制新轨迹与相邻轨迹之间的间距。

2. trackc.make_spec
  • trackc.make_spec 定义整个图像的大小，允许用户从上到下或从左到右设置一组子图，可以通过 height_ratios 或 width_ratios 控制各部分的比例。

可用轨迹

本节包括多个快速教程，展示了使用 trackc 进行组学数据可视化的功能。

02 安装

trackc 运行需要 Python 版本 >= 3.8。

通过 PyPI 安装

复制代码 
使用以下命令安装 trackc：

pip install trackc
更新 trackc 的命令：

pip install --upgrade trackc
指定版本和安装源的命令：

pip install -i https://pypi.org/simple trackc==版本号

开发版本

复制代码 
从 GitHub 安装 trackc 的命令：

pip install git+https://github.com/seqyuan/trackc@main
03 heatmap
复制代码 
conf-hicmap_1.yml等很多
复制代码 
trackc:
  - ax: t1
    height: 2
    track_type: hicmap
    track_para:
      mat:
        method: extractContactRegions
        clr: /path/GSM4417639_SK-N-DZ_C6BC81F2_b38d5.mcool::/resolutions/5120000
        row_regions:
          - chr6
          - chr8
      mapC:
        cmap: plasma
复制代码 
运行
复制代码 
trackc cli conf-hicmap_1.yml -s 4,1 -o chr6_chr8.pdf

conf-hicmap_4.yml

复制代码 
trackc:
  - ax: t1
    height: 1
    track_type: hicmap
    track_para:
      mat:
        method: extractCisContact
        clr: ./github/seqyuan/trackc_data/examples/HPDE6C7.chr18.mcool::/resolutions/25000
      mapC:
        map_type: tri
        height: 50
        label: tri
  - ax: t2
    height: 1
    hspace: 0.3
    track_type: hicmap
    track_para:
      mat2:
        method: extractCisContact
        clr: ./github/seqyuan/trackc_data/examples/HPDE6C7.chr18.mcool::/resolutions/25000
      mapC:
        map_type: tri
        height: 50
        label: tri-mat2
  - ax: t3
    height: 2
    hspace: 0.3
    track_type: hicmap
    track_para:
      mat:
        method: extractCisContact
        clr: ./github/seqyuan/trackc_data/examples/HPDE6C7.chr18.mcool::/resolutions/25000
      mapC:
        map_type: tri
        height: 50
        label: tri-symmetric
        symmetric: True

  - ax: t4
    height: 1
    hspace: 0.3
    track_type: hicmap
    track_para:
      mat:
        method: extractCisContact
        clr:./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000
        extend: 50
      mapC:
        cmap: PuBu
        map_type: rec
        height: 50
        label: rec
  - ax: t5
    height: 1
    hspace: 0.3
    track_type: hicmap
    track_para:
      mat2:
        method: extractCisContact
        clr: ./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000
        extend: 50
      mapC:
        cmap: cividis
        map_type: rec
        height: 50
        label: rec-mat2
  - ax: t6
    height: 2
    hspace: 0.3
    track_type: hicmap
    track_para:
      mat2:
        method: extractCisContact
        clr: ./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000
        extend: 50
      mapC:
        cmap: magma
        map_type: rec
        height: 50
        symmetric: True
        logdata: True
        label: rec-symmetric

  - ax: t1
    track_type: scale_track
    track_para:
      scale_adjust: Mb
      tick_pos: top
      ratio2ax: 1

运行

复制代码 
trackc cli conf-hicmap_4.yml -r 18:10000000-14500000
04 rearranged_interactions

ectopic_interactions.yaml

复制代码 
trackc:
  - ax: t1
    height: 2
    track_type: hicmap
    track_para:
      mat:
        method: extractContactRegions
        clr: ./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/50000
        balance: False
        row_regions: 18:45000000-78077248
      mapC:
        map_type: tri
        label: tumor res=50k
        cmap: PuBu
        maxrange: 70
        ax_on: False
  - ax: t1
    track_type: scale_track
    track_para:
      region: 18:45000000-78077248
      scale_adjust: Mb
      tick_pos: top
      ratio2ax: 0.3
      space: 0.01
      tick_fontsize: 6
      label_fontsize: 7
  - ax: t2
    height: 0.2
    track_type: zoomin
    track_para:
      raw_regions: 18:45000000-78077248
      zoomin_regions:
        - 18:47340000-50370000
        - 18:61140000-63630000
        - 18:74030000-77560000
      line_on: False
      fill: True
      alpha: 0.5

  - ax: t3
    height: 0.4
    track_type: multi_scale_track
    track_para:
      regions:
        - 18:47340000-50370000
        - 18:61140000-63630000
        - 18:74030000-77560000
      scale_adjust: Mb
      intervals: 2
      tick_fontsize: 7

  - ax: t4
    height: 2
    track_type: hicmap
    track_para:
      mat2:
        method: extractContactRegions
        clr: ./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000
        balance: False
        row_regions:
          - 18:47340000-50370000
          - 18:61140000-63630000
          - 18:74030000-77560000
      mapC:
        map_type: tri
        label: ectopic interactions
        cmap: PuBu
        ax_on: False
        maxrange: 100

  - ax: t5
    height: 0.2
    track_type: zoomin
    track_para:
      raw_regions:
        - 18:47340000-50370000
        - 18:61140000-63630000
        - 18:74030000-77560000
      zoomin_regions:
        - 18:47400000-48280000
        - 18:75280000-74030000
      line_on: False
      fill: True
      alpha: 0.5

  - ax: t6
    height: 0.4
    track_type: multi_scale_track
    track_para:
      regions:
        - 18:47400000-48280000
        - 18:75280000-74030000
      scale_adjust: Mb
      intervals: 2
      tick_fontsize: 7
  - ax: t7
    height: 1.25
    track_type: hicmap
    track_para:
      mat:
        method: extractContactRegions
        clr: ./github/seqyuan/trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000
        balance: False
        row_regions:
          - 18:47400000-48280000
          - 18:75280000-74030000
      mapC:
        map_type: tri
        label: neo tad
        cmap: PuBu
        height: 40
        ax_on: False
        maxrange: 200
        minrange: 10

  - ax: t8
    height: 0.5
    track_type: bw_track
    track_para:
      bw: ./github/seqyuan/trackc_data/examples/GSM3178671_BxPC3_H3K27ac.bigwig
      regions:
        - 18:47400000-48280000
        - 18:75280000-74030000
      label: H3K27ac
      binsize: 2000
      color:
       - '#5B7695'

  - ax: t9
    height: 2.8
    track_type: gene_track
    track_para:
      bed12: ./github/seqyuan/trackc_data/examples/hg19_chr18.gene.bed12
      regions:
        - 18:47400000-48280000
        - 18:75280000-74030000
      line: 12
      gene_fontszie: 8

运行

复制代码 
trackc cli ectopic_interactions.yaml -s 5,1
05 Virtual 4C

4C.yaml

复制代码 
trackc:
  - ax: t1
    height: 4
    track_type: hicmap
    track_para:
      mat:
        method: extractContactRegions
        clr: ./github/seqyuan/trackc_data/tutorials/4C/GSM4604287_1360.sub.cool
        balance: False
      mapC:
        map_type: tri
        label: AML 1360
        cmap: PuBu
        maxrange: 6
        minrange: -6

  - ax: t2
    height: 0.7
    track_type: gene_track
    track_para:
      bed12: ./github/seqyuan/trackc_data/tutorials/4C/GRCh38.gene.chr8chr14.bed12
      line: 3
      gene_fontszie: 10
      show_label:
        - MYC
  - ax: t3
    height: 0.7
    track_type: bw_track
    track_para:
      bw: ./github/seqyuan/trackc_data/tutorials/4C/GSM4604189_H3K27ac.bw
      maxrange: 10
      label: H3K27ac
      binsize: 10000
      invert_y: False
  - ax: t4
    height: 0.6
    track_type: virtual4C
    track_para:
      clr: ./github/seqyuan/trackc_data/tutorials/4C/GSM4604287_1360.sub.cool
      maxrange: 10
      label: Virtual 4C
      target_color: tab:red
      track_type: line
      target: chr8:127735434-127735435

  - ax: t5
    height: 0.6
    track_type: multi_scale_track
    track_para:
      scale_adjust: Mb
      intervals: 1
      tick_fontsize: 8

运行

复制代码 
trackc cli 4C.yaml -s 6,0.75 -r 'chr8:127000000-129200000 chr14:96500000-99300000'
06 常用命令行
复制代码 
import trackc as tc

ten = tc.tenon(figsize=(8,1))
ten.add(pos='bottom', height=3.6)
ten.add(pos='bottom', height=1)

chr6_len = clr.chromsizes['chr6']
chr8_len = clr.chromsizes['chr8']

mat = tc.tl.extractContactRegions(clr='GSM4417639_SK-N-DZ_C6BC81F2_b38d5.mcool::/resolutions/5120000', row_regions=['chr6', 'chr8'])
tc.pl.mapC(ax=ten.axs(0), mat=mat.cmat, map_type='triangle', maxrange=5000)

tc.pl.multi_scale_track(ax=ten.axs(1), regions=['chr6:0-{0}'.format(chr6_len), 'chr8:0-{0}'.format(chr8_len)],
                        scale_adjust='Mb', intervals=2, tick_rotation=0)
# tc.savefig('hicmap.pdf')
相关推荐
PH_modest9 分钟前
【Qt跬步积累】—— 初识Qt
开发语言·qt
咔咔一顿操作29 分钟前
第七章 Cesium 3D 粒子烟花效果案例解析：从原理到完整代码
人工智能·3d·信息可视化·cesium
怀旧,36 分钟前
【C++】18. 红⿊树实现
开发语言·c++
多恩Stone1 小时前
【3DV 进阶-2】Hunyuan3D2.1 训练代码详细理解下-数据读取流程
人工智能·python·算法·3d·aigc
xiaopengbc1 小时前
在 Python 中实现观察者模式的具体步骤是什么？
开发语言·python·观察者模式
Python大数据分析@1 小时前
python用selenium怎么规避检测？
开发语言·python·selenium·网络爬虫
ThreeAu.1 小时前
Miniconda3搭建Selenium的python虚拟环境全攻略
开发语言·python·selenium·minicoda·python环境配置
偷心伊普西隆1 小时前
Python EXCEL 理论探究：格式转换时处理缺失值方法
python·excel
zhangfeng11331 小时前
R 语法高亮为什么没有，是需要安装专用的编辑软件，R语言自带的R-gui 功能还是比较简单
开发语言·r语言
热门推荐
01UV安装并设置国内源02A股预测还能更准？开源大模型Kronos带你跑通预测+回测全流程032025 年高教社杯全国大学生数学建模竞赛C 题 NIPT 的时点选择与胎儿的异常判定 完整成品思路模型代码分享，全网首发高质量！！！042025年数学建模国赛C题超详细解题思路05UV 工具安装与国内镜像源配置指南06不再让Windows更新！&Edge游戏助手卸载及关闭自动更新07KGG转MP3工具|非KGM文件|解密音频08突破百度网盘的下载限速，两种方法教会你【超详细】09Linux下V2Ray安装配置指南10教你如何认证 Gemini 教育优惠的二次验证，薅个 1年的 Gemini Pro 会员