当前位置：首页 > news >正文

Wig/BigWig 格式实战：3步从BAM文件生成UCSC基因组浏览器轨迹

news 2026/7/6 11:59:52

从BAM到BigWig：基因组浏览器轨迹生成全流程实战指南

在ChIP-seq或RNA-seq等高通量测序数据分析中，可视化是理解数据分布特征的关键环节。本文将手把手带您完成从原始BAM文件到UCSC基因组浏览器兼容的BigWig格式的完整转换流程，提供可直接复用的Shell脚本和关键参数优化建议。

1. 环境准备与工具安装

开始前需要确保Linux系统已安装以下工具：

samtools：用于处理BAM文件和计算测序深度
UCSC工具集：包含wigToBigWig等格式转换工具
基因组染色体尺寸文件：记录各染色体长度信息

安装UCSC工具集的命令如下：

# 下载wigToBigWig工具 wget http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/wigToBigWig chmod +x wigToBigWig # 获取染色体尺寸文件（以hg38为例） fetchChromSizes hg38 > hg38.chrom.sizes

提示：对于不同基因组版本，需替换hg38为相应版本号（如mm10、rn6等）。染色体尺寸文件也可从UCSC官网直接下载。

2. 从BAM到Wig：测序深度提取与格式转换

2.1 使用samtools计算测序深度

核心命令通过管道组合实现深度计算与格式转换：

samtools depth -aa your_file.bam | \ awk 'BEGIN{OFS="\t"; print "variableStep chrom="$1" span=1"} {print $2,$3}' > output.wig

参数说明：

-aa：输出所有位点（包括零覆盖区域）
awk部分实现格式转换：
- BEGIN块输出Wig头信息
- 主块输出位置和深度值

2.2 处理大文件的优化策略

对于全基因组数据，建议采用分区处理：

# 按染色体并行处理 for chr in $(cut -f1 hg38.chrom.sizes); do samtools view -b your_file.bam $chr | \ samtools depth -aa | \ awk -v chr=$chr 'BEGIN{OFS="\t"; print "variableStep chrom="chr" span=1"} {print $2,$3}' > ${chr}.wig & done wait

3. Wig到BigWig：二进制压缩与索引生成

3.1 格式转换基础命令

wigToBigWig input.wig hg38.chrom.sizes output.bw

关键参数对比：

参数	作用	推荐值
-clip	修剪超出染色体范围的数据	建议启用
-fixedSummaries	生成固定窗口摘要	大文件时建议使用
-uncached	禁用本地缓存	服务器环境建议启用

3.2 高级应用：窗口化处理

对于全基因组数据，可先生成窗口化Wig文件：

# 生成100bp窗口的Wig samtools depth -aa your_file.bam | \ awk 'BEGIN{OFS="\t"; window=100; last_chr=""; pos=0} { if($1!=last_chr){pos=0; last_chr=$1; print "variableStep chrom="$1" span="window} if($2>=pos+window){ pos=int($2/window)*window; print pos,$3 } }' > windowed.wig

4. 可视化优化与实战技巧

4.1 UCSC基因组浏览器轨迹配置

创建自定义轨迹描述文件（track.txt）：

track type=bigWig name="My ChIP-seq Data" description="ChIP-seq peaks from experiment X" bigDataUrl=http://yourdomain.com/path/to/output.bw color=50,150,200 altColor=200,100,50 maxHeightPixels=100:50:20 visibility=full autoScale=on

4.2 常见问题排查

问题1：坐标超出染色体范围

# 检查染色体最大坐标 cut -f1,2 hg38.chrom.sizes | while read chr len; do awk -v len=$len '$1=="'$chr'" && $2>len{print "Error:"$0}' input.wig done

问题2：内存不足处理

# 使用split分割大文件 split -l 1000000 bigfile.wig chunk_ for f in chunk_*; do wigToBigWig $f hg38.chrom.sizes ${f}.bw done

5. 自动化脚本集成

完整处理脚本示例：

#!/bin/bash BAM=$1 GENOME=$2 # hg38/mm10 etc OUTDIR=$3 # 获取工具路径 WIGTOBW=$(which wigToBigWig) SAMTOOLS=$(which samtools) # 创建临时目录 mkdir -p ${OUTDIR}/tmp # 处理流程 echo "Processing $BAM..." $SAMTOOLS depth -aa $BAM | \ awk 'BEGIN{OFS="\t"; print "variableStep chrom="$1" span=1"} {print $2,$3}' > ${OUTDIR}/tmp/raw.wig # 格式转换 $WIGTOBW ${OUTDIR}/tmp/raw.wig \ <(fetchChromSizes $GENOME) \ ${OUTDIR}/result.bw 2>&1 | tee ${OUTDIR}/conversion.log # 清理临时文件 rm -rf ${OUTDIR}/tmp echo "Done. Output: ${OUTDIR}/result.bw"

使用方式：