pyGenomeTracks 使用教程

1. 项目介绍

pyGenomeTracks 是一个用于绘制高质量、高度可定制的基因组浏览器轨迹的 Python 模块。它旨在生成美观且可重复的基因组浏览器快照，支持多种基因组数据类型，包括 bigwig、bed/gtf、bedgraph、epilogos、narrow peaks、links、Hi-C 矩阵等。pyGenomeTracks 不仅适用于科研人员，还可以集成到自动化工作流中，提供强大的基因组数据可视化功能。

2. 项目快速启动

安装

pyGenomeTracks 支持 Python 3.8 及以上版本。推荐使用 conda 进行安装：

conda create -n pygenometracks -c bioconda -c conda-forge pygenometracks

也可以使用 pip 进行安装：

pip install pygenometracks

使用

首先，使用 make_tracks_file 创建一个配置文件：

make_tracks_file --trackFiles <file1.bed> <file2.bw> -o tracks.ini

然后，使用 pyGenomeTracks 绘制基因组区域：

pyGenomeTracks --tracks tracks.ini --region chr2:10,000,000-11,000,000 --outFileName nice_image.pdf

3. 应用案例和最佳实践

案例1：绘制 Hi-C 矩阵

假设你有一个 Hi-C 矩阵文件 hic_matrix.hic，你可以使用以下命令绘制特定区域的 Hi-C 矩阵：

pyGenomeTracks --tracks tracks.ini --region chr1:10,000,000-20,000,000 --outFileName hic_image.pdf

案例2：绘制多个基因组轨迹

你可以将多个基因组轨迹文件（如 bigwig、bed 文件）合并到一个配置文件中，并一次性绘制多个轨迹：

make_tracks_file --trackFiles file1.bw file2.bed file3.bedgraph -o multi_tracks.ini
pyGenomeTracks --tracks multi_tracks.ini --region chr3:5,000,000-15,000,000 --outFileName multi_tracks_image.pdf

4. 典型生态项目

CoolBox

CoolBox 是一个用于 Jupyter Notebooks 的交互式基因组数据浏览器，可以与 pyGenomeTracks 结合使用，提供更丰富的交互式基因组数据可视化体验。

Galaxy 集成

Galaxy 是一个开源的生物信息学平台，pyGenomeTracks 可以集成到 Galaxy 中，提供图形用户界面来创建基因组浏览器轨迹图，并可以将其纳入自动化工作流。

通过这些生态项目的支持，pyGenomeTracks 不仅可以在本地环境中使用，还可以扩展到更广泛的生物信息学工作流中，提供强大的基因组数据可视化能力。