pbsv 项目使用教程

1. 项目介绍

pbsv 是由 Pacific Biosciences 开发的一套工具，用于从 PacBio 单分子实时测序（SMRT）读取中调用和分析二倍体基因组中的结构变异（SV）。pbsv 工具集支持插入、删除、倒位、重复和易位等结构变异的调用，并且支持单样本调用和多样本联合调用。

pbsv 的主要功能包括：

• 调用结构变异（SV）
• 分析结构变异
• 支持单样本和多样本联合调用

2. 项目快速启动

2.1 安装 pbsv

pbsv 可以通过 bioconda 进行安装。以下是安装步骤：

# 安装 miniconda3
module load miniconda3

# 创建并激活 conda 环境
conda create -n mypy_py27 python=2.7
source activate mypy_py27

# 安装 pbsv
conda install -c bioconda pbsv

2.2 使用 pbsv 进行结构变异调用

以下是一个简单的 pbsv 工作流程示例：

1. 对齐 PacBio 读取到参考基因组

   使用 pbmm2 对齐工具将 PacBio 读取对齐到参考基因组。

   pbmm2 align ref.fa movie1.subreads.bam ref_movie1.bam --sort --median-filter --sample sample1
   

2. 发现结构变异的特征

   使用 pbsv discover 工具从对齐的 BAM 文件中发现结构变异的特征。

   pbsv discover ref_movie1.bam ref_sample1.svsig.gz
   

3. 调用结构变异并分配基因型

   使用 pbsv call 工具从结构变异特征文件中调用结构变异并分配基因型。

   pbsv call ref.fa ref_sample1.svsig.gz ref_var.vcf
   

3. 应用案例和最佳实践

3.1 应用案例

pbsv 广泛应用于基因组学研究中，特别是在结构变异的检测和分析方面。例如，研究人员可以使用 pbsv 来检测和分析人类基因组中的结构变异，以研究基因组多样性和疾病相关性。

3.2 最佳实践

• 数据质量控制：在使用 pbsv 之前，确保输入数据的质量，包括测序读取的准确性和覆盖度。
• 参数优化：根据具体的研究需求，调整 pbsv 的参数以获得最佳的结构变异调用结果。
• 结果验证：使用其他工具或实验方法验证 pbsv 调用的结构变异结果，以确保结果的准确性。

4. 典型生态项目

pbsv 作为 PacBio 结构变异调用和分析工具的一部分，与其他 PacBio 工具和生态系统项目紧密集成。以下是一些典型的生态项目：

• SMRT Link：PacBio 的图形用户界面（GUI）工具，集成了 pbsv 和其他 PacBio 分析工具，方便用户进行结构变异调用和分析。
• pbmm2：PacBio 的读取对齐工具，与 pbsv 配合使用，用于将 PacBio 读取对齐到参考基因组。
• bioconda：生物信息学软件包管理工具，用于安装和管理 pbsv 及其依赖项。

通过这些生态项目，pbsv 可以更好地与其他工具和资源集成，提供更全面的基因组分析解决方案。