nanozoo/rsem

RSEM v1.11 Docker镜像文档

镜像概述

本镜像包含RSEM（RNA-Seq by Expectation-Maximization）v1.11版本，RSEM是一款用于RNA测序数据定量分析的工具，主要功能是估计基因和转录本的表达水平。该镜像由Martin维护，基于https://github.com/deweylab/RSEM%E6%9E%84%E5%BB%BA%EF%BC%8C%E9%81%B5%E5%BE%AAGPL3%EF%BC%88GNU General Public License v3）许可证。

核心功能与特性

• 转录组参考索引构建：支持从基因组序列和注释文件构建转录组参考索引，为后续定量分析提供基础
• RNA-Seq数据定量：通过Expectation-Maximization算法估计基因和转录本的表达水平（如FPKM、TPM等）
• 多输入格式支持：兼容FASTQ等常用测序数据格式，支持单端和双端测序数据
• 标准化输出：生成标准化的表达矩阵，便于下游差异表达分析等生物信息学流程
• 版本稳定性：基于v1.11稳定版本，确保分析结果的可靠性和可重复性

使用场景与适用范围

• 转录组学研究：用于RNA-Seq实验中基因和转录本表达水平的定量分析
• 差异表达分析：作为差异表达基因分析前的表达定量步骤，为DESeq2、edgeR等工具提供输入数据
• 多组学整合：在多组学研究中，提供RNA层面的表达数据，支持与基因组、蛋白质组数据的整合分析
• 生物信息学流程集成：可作为自动化分析流程（如RNA-Seq分析管道）中的定量模块

使用方法与配置说明

基本使用

通过以下命令可查看RSEM帮助信息：

bash
docker run --rm [镜像名称] rsem-help

数据挂载

分析时需将数据目录挂载至容器内，建议使用-v参数挂载本地数据目录：

bash
docker run --rm -v /本地数据目录:/data [镜像名称] [RSEM命令] --input /data/输入文件 --output /data/输出目录

典型操作示例

1. 构建转录组参考索引

bash
docker run --rm -v /本地参考数据:/ref -v /输出索引目录:/index [镜像名称] \
  rsem-prepare-reference \
  --bowtie2 \  # 使用Bowtie2作为比对工具
  --reference-name hg38 \  # 参考名称
  /ref/genome.fa \  # 基因组序列文件
  /ref/transcripts.gtf \  # 转录本注释文件
  /index/rsem_index  # 输出索引前缀

2. RNA-Seq数据定量分析（双端数据）

bash
docker run --rm -v /本地测序数据:/fastq -v /索引目录:/index -v /结果目录:/result [镜像名称] \
  rsem-calculate-expression \
  --bowtie2 \
  --paired-end \  # 双端测序数据
  /fastq/sample_R1.fastq.gz \  # R1 reads
  /fastq/sample_R2.fastq.gz \  # R2 reads
  /index/rsem_index \  # 参考索引前缀
  /result/sample_expression  # 输出结果前缀

主要配置参数说明

输出文件说明

定量分析后将在输出目录生成以下主要文件：

• sample_expression.genes.results：基因水平表达结果
• sample_expression.isoforms.results：转录本水平表达结果
• sample_expression.bam：比对结果BAM文件（若指定--output-genome-bam）

许可证信息

本镜像遵循GPL3（GNU General Public License v3）许可证，详细信息可参考GPL3官方文档。