
摘要
非靶向宏基因组(又称鸟枪宏基因组,whole-metagenome sequencing,WMS)测序可全面解析微生物群落的分类与功能组成,并支持不依赖培养地重构宏基因组组装基因组(metagenome-assembled genomes,MAGs)。然而,WMS 分析流程复杂、工具割裂、缺乏标准化,严重制约结果的可重复性。作者介绍了 MOSHPIT(MOdular SHotgun metagenome Pipelines with Integrated provenance Tracking)——一套构建于 QIIME 2/rachis 框架(QIIME 2/rachis framework,Q2F)之上的模块化鸟枪宏基因组分析套件。MOSHPIT 整合了在 CAMI II 与 LEMMI 基准中表现优异的组装、分箱、注释与统计工具,并内置完整的数据溯源(provenance)跟踪、检查点可视化、并行化与流程续跑能力。该套件支持读(read)、重叠群(contig)与 MAG 三种尺度的分类及功能分析,使不同经验水平的研究者均可在统一生态系统中完成端到端、可重复的 WMS 分析,从而降低宏基因组学前沿研究的门槛并加速发现。
keywords:鸟枪宏基因组;QIIME 2;宏基因组组装基因组;数据溯源;可重复分析;CAMI II;LEMMI;分类与功能注释
文献信息
- Ziemski, M., Gehret, L., Simard, A., Castro Dau, S., Risch, V., Grabocka, D., Matzzoros, C., Wood, C., Momo Cabrera, P., Hernandez-Velázquez, R., Scharer, M. R., Sandberg, F., Herman, C., Evans, K., Robeson, M. S. II, Bolyen, E., Caporaso, J. G., & Bokulich, N. A. (2026). MOSHPIT: Accessible, reproducible metagenome data science on the QIIME 2 framework. Cell Genomics, 6, 101366.
- 期刊:Cell Genomics(Cell Press 旗下,IF≈10–11,2025–2026 年 JCR 基因组学 Q1)
- 在线发表:2026 年 11 月 11 日;接收:2026 年 9 月 1 日
研究总结
#mermaid-svg-fJZHFZY5X9mohVNu{font-family:\”trebuchet ms\”,verdana,arial,sans-serif;font-size:16px;fill:#333;}@keyframes edge-animation-frame{from{stroke-dashoffset:0;}}@keyframes dash{to{stroke-dashoffset:0;}}#mermaid-svg-fJZHFZY5X9mohVNu .edge-animation-slow{stroke-dasharray:9,5!important;stroke-dashoffset:900;animation:dash 50s linear infinite;stroke-linecap:round;}#mermaid-svg-fJZHFZY5X9mohVNu .edge-animation-fast{stroke-dasharray:9,5!important;stroke-dashoffset:900;animation:dash 20s linear infinite;stroke-linecap:round;}#mermaid-svg-fJZHFZY5X9mohVNu .error-icon{fill:#552222;}#mermaid-svg-fJZHFZY5X9mohVNu .error-text{fill:#552222;stroke:#552222;}#mermaid-svg-fJZHFZY5X9mohVNu .edge-thickness-normal{stroke-width:1px;}#mermaid-svg-fJZHFZY5X9mohVNu .edge-thickness-thick{stroke-width:3.5px;}#mermaid-svg-fJZHFZY5X9mohVNu .edge-pattern-solid{stroke-dasharray:0;}#mermaid-svg-fJZHFZY5X9mohVNu .edge-thickness-invisible{stroke-width:0;fill:none;}#mermaid-svg-fJZHFZY5X9mohVNu .edge-pattern-dashed{stroke-dasharray:3;}#mermaid-svg-fJZHFZY5X9mohVNu .edge-pattern-dotted{stroke-dasharray:2;}#mermaid-svg-fJZHFZY5X9mohVNu .marker{fill:#333333;stroke:#333333;}#mermaid-svg-fJZHFZY5X9mohVNu .marker.cross{stroke:#333333;}#mermaid-svg-fJZHFZY5X9mohVNu svg{font-family:\”trebuchet ms\”,verdana,arial,sans-serif;font-size:16px;}#mermaid-svg-fJZHFZY5X9mohVNu p{margin:0;}#mermaid-svg-fJZHFZY5X9mohVNu .edge{stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .section–1 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section–1 path,#mermaid-svg-fJZHFZY5X9mohVNu .section–1 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section–1 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section–1 path{fill:hsl(240, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section–1 text{fill:#ffffff;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon–1{font-size:40px;color:#ffffff;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge–1{stroke:hsl(240, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth–1{stroke-width:17;}#mermaid-svg-fJZHFZY5X9mohVNu .section–1 line{stroke:hsl(60, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-0 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-0 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-0 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-0 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-0 path{fill:hsl(60, 100%, 73.5294117647%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-0 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-0{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-0{stroke:hsl(60, 100%, 73.5294117647%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-0{stroke-width:14;}#mermaid-svg-fJZHFZY5X9mohVNu .section-0 line{stroke:hsl(240, 100%, 83.5294117647%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-1 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-1 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-1 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-1 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-1 path{fill:hsl(80, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-1 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-1{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-1{stroke:hsl(80, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-1{stroke-width:11;}#mermaid-svg-fJZHFZY5X9mohVNu .section-1 line{stroke:hsl(260, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-2 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-2 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-2 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-2 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-2 path{fill:hsl(270, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-2 text{fill:#ffffff;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-2{font-size:40px;color:#ffffff;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-2{stroke:hsl(270, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-2{stroke-width:8;}#mermaid-svg-fJZHFZY5X9mohVNu .section-2 line{stroke:hsl(90, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-3 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-3 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-3 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-3 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-3 path{fill:hsl(300, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-3 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-3{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-3{stroke:hsl(300, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-3{stroke-width:5;}#mermaid-svg-fJZHFZY5X9mohVNu .section-3 line{stroke:hsl(120, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-4 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-4 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-4 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-4 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-4 path{fill:hsl(330, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-4 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-4{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-4{stroke:hsl(330, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-4{stroke-width:2;}#mermaid-svg-fJZHFZY5X9mohVNu .section-4 line{stroke:hsl(150, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-5 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-5 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-5 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-5 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-5 path{fill:hsl(0, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-5 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-5{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-5{stroke:hsl(0, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-5{stroke-width:-1;}#mermaid-svg-fJZHFZY5X9mohVNu .section-5 line{stroke:hsl(180, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-6 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-6 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-6 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-6 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-6 path{fill:hsl(30, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-6 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-6{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-6{stroke:hsl(30, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-6{stroke-width:-4;}#mermaid-svg-fJZHFZY5X9mohVNu .section-6 line{stroke:hsl(210, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-7 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-7 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-7 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-7 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-7 path{fill:hsl(90, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-7 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-7{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-7{stroke:hsl(90, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-7{stroke-width:-7;}#mermaid-svg-fJZHFZY5X9mohVNu .section-7 line{stroke:hsl(270, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-8 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-8 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-8 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-8 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-8 path{fill:hsl(150, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-8 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-8{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-8{stroke:hsl(150, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-8{stroke-width:-10;}#mermaid-svg-fJZHFZY5X9mohVNu .section-8 line{stroke:hsl(330, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-9 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-9 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-9 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-9 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-9 path{fill:hsl(180, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-9 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-9{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-9{stroke:hsl(180, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-9{stroke-width:-13;}#mermaid-svg-fJZHFZY5X9mohVNu .section-9 line{stroke:hsl(0, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-10 rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-10 path,#mermaid-svg-fJZHFZY5X9mohVNu .section-10 circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-10 polygon,#mermaid-svg-fJZHFZY5X9mohVNu .section-10 path{fill:hsl(210, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-10 text{fill:black;}#mermaid-svg-fJZHFZY5X9mohVNu .node-icon-10{font-size:40px;color:black;}#mermaid-svg-fJZHFZY5X9mohVNu .section-edge-10{stroke:hsl(210, 100%, 76.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .edge-depth-10{stroke-width:-16;}#mermaid-svg-fJZHFZY5X9mohVNu .section-10 line{stroke:hsl(30, 100%, 86.2745098039%);stroke-width:3;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled,#mermaid-svg-fJZHFZY5X9mohVNu .disabled circle,#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:lightgray;}#mermaid-svg-fJZHFZY5X9mohVNu .disabled text{fill:#efefef;}#mermaid-svg-fJZHFZY5X9mohVNu .section-root rect,#mermaid-svg-fJZHFZY5X9mohVNu .section-root path,#mermaid-svg-fJZHFZY5X9mohVNu .section-root circle,#mermaid-svg-fJZHFZY5X9mohVNu .section-root polygon{fill:hsl(240, 100%, 46.2745098039%);}#mermaid-svg-fJZHFZY5X9mohVNu .section-root text{fill:#ffffff;}#mermaid-svg-fJZHFZY5X9mohVNu .section-root span{color:#ffffff;}#mermaid-svg-fJZHFZY5X9mohVNu .section-2 span{color:#ffffff;}#mermaid-svg-fJZHFZY5X9mohVNu .icon-container{height:100%;display:flex;justify-content:center;align-items:center;}#mermaid-svg-fJZHFZY5X9mohVNu .edge{fill:none;}#mermaid-svg-fJZHFZY5X9mohVNu .mindmap-node-label{dy:1em;alignment-baseline:middle;text-anchor:middle;dominant-baseline:middle;text-align:center;}#mermaid-svg-fJZHFZY5X9mohVNu :root{–mermaid-font-family:\”trebuchet ms\”,verdana,arial,sans-serif;}
MOSHPIT
背景
16S只看潜能
WMS功能更强但流程碎
QIIME2缺WMS原生支持
设计
Q2F插件体系
组装分箱注释一体化
溯源与可重放
多用户界面
核心插件
q2-assembly
q2-mag
q2-annotate
辅助插件
q2-fastp
q2-skani
q2-gunc
q2-amrfinderplus
q2-viromics
验证
CAMI II人体微生态
TARA海洋
多研究人肠
可可发酵
意义
降低WMS门槛
端到端可重复
读重叠群MAG三尺度
背景介绍
标记基因扩增子测序(如 16S rRNA)因可扩展性强而被广泛使用,但受引物选择影响,且仅能间接推断功能潜力。相比之下,WMS 可同时获得功能基因与分类组成、达到亚种级分辨率,并支持不依赖培养地重构 MAG。问题在于:WMS 工具链极其碎片化——组装、分箱、质量控制、分类、功能注释、统计与可视化往往由不同团队维护,参数、数据库版本与中间格式互不兼容,导致“同一数据不同实验室结果难复现”。
自 2018 年发布以来,QIIME 2 以可访问性、互操作性与可重复性获得广泛用户群,但 Q2F 长期缺乏对非靶向 WMS 的原生支持。MOSHPIT 正是为弥补这一缺口而设计:它在 Q2F 上封装经过基准验证的最佳工具,用统一Artifact(.qza)与溯源机制把读级、重叠群级和 MAG 级分析串成一条可审计、可重放、可扩展的管道。
重要结果
MOSHPIT 的工作流与插件架构
MOSHPIT 包含三个核心插件与多个可选附属插件。核心插件中,q2-assembly 负责重叠群组装、质量评估与读回贴;q2-mag 负责分箱(binning)、MAG 质量控制与去冗余;q2-annotate 负责分类(Kraken 2、Kaiju)与功能注释(eggNOG-mapper、HUMANN)。附属插件覆盖质控(q2-fastp)、平均核苷酸一致性分析(q2-skani)、嵌合体检测(q2-gunc)、差异丰度(q2-deseq2)、病毒检测(q2-viromics、q2-virsorter2)以及抗性基因注释(q2-amrfinderplus、q2-rgi)。新数据类型、格式与转换器被加入 Q2F,使这些插件可无缝衔接。

图 1. MOSHPIT 主要输入输出类型与插件架构。A 面板展示原始读可由 Kaiju 做分类注释,Kraken 2 可作用于原始读、重叠群或去冗余 MAG;eggNOG-mapper 可用于重叠群或去冗余 MAG 的功能注释,HUMANN 用于读级功能谱。B 面板列出与 MOSHPIT 工作流兼容的插件,星号表示社区插件。
基于 CAMI II 的组装与 MAG 恢复基准
作者用 CAMI II 人体微生态玩具数据集(49 个样本,每样本 16.6M 读)评估标准 MOSHPIT 流程。组装连续性因身体部位而异:气道(airways)中位最长重叠群约 0.45 Mbp,胃肠道(gastrointestinal tract)约 0.9 Mbp;但所有部位的组装准确度都很高,中位重叠群 ANI 超过 95%,胃肠样本接近 99%。去冗余 MAG 与金标准比对显示,分类学恢复力强:属级(level 7)R² 超过 0.85,科级以上 R² 高于 0.9;分类学准确率(taxon accuracy rate,TAR)与检出率(taxon detection rate,TDR)在更细分辨率下自然下降,但整体群落重构保持高精度。

图 2. CAMI II 基准结果。A 面板为各身体部位(A 气道、G 胃肠、O 口腔、S 皮肤、U 泌尿生殖)的组装质量指标。B 面板为去冗余 MAG 与金标准数据集的群落恢复评估,误差棒表示均值标准误。
TARA Oceans 海洋原核微生物分析
为展示大规模 WMS 处理能力,作者用 q2-fondue 从 NCBI SRA 拉取 TARA Oceans 原核组分数据,再用 q2-annotate 中 Kraken 2 做读级分类,并接入 Q2F 的多样性分析。地图面板展示全球样品的香农多样性,四个放大区显示不同深度样品的分类指派;Bray-Curtis 主坐标散点图突出深海样品在组成上的相似性。该案例证明 MOSHPIT 可与既有 QIIME 2 插件互操作,从数据获取到多样性分析都在同一框架内完成。

图 3. TARA Oceans 分析概览。地图描绘全球采集样品的香农多样性,并放大四个位点展示分类指派;Bray-Curtis 主坐标散点图显示深海样品间的组成相似性。
多研究人肠宏基因组中的工业化梯度与耐药组
作者整合 1,578 个人肠样本、覆盖采猎、渔猎、农耕与不同工业化程度人群。重叠群用 MEGAHIT 组装,eggNOG-mapper 注释 KEGG KO 术语,读回贴后算每样本 KO 丰度。PERMANOVA(adonis2)显示人群类型显著解释样本间方差(R²=0.15,p=0.001)。ANCOM-BC2 发现工业化大城市与采猎人群间差异最显著的 KO 之一为 K08169(log-fold change=3.57,p<0.001),该基因编码主要易化超家族(major facilitator superfamily)外排转运蛋白,赋予多重耐药表型。结果与“工业化人群肠道耐药组更富集”的既往结论一致。

图 4. 人肠工业化梯度功能分析。主坐标散点图按城市化水平人群的功能谱聚类;箱线图显示多重耐药转运蛋白 K08169 丰度随城市化水平上升而增加。
可可发酵中的读级与 MAG 级一致性分析
对 14 个可可发酵时间序列样本(两个品种各 7 个),MOSHPIT 同时跑读级(Kraken 2、Kaiju、mOTUs 3)与 MAG 级流程。所有分类器都再现了发酵过程中分类多样性的下降,并检出 Acetobacter 与 Pantoea 丰度降低。MAG 功能注释给出每样本 CAZyme(碳水化合物活性酶)计数:发酵进程中糖基转移酶(glycosyltransferase,GT)类上升,糖苷水解酶(glycoside hydrolase,GH)类下降,表明发酵后期微生物群落在碳水化合物代谢上发生功能倾斜。该案例验证了 MOSHPIT 在食品发酵生境中的端到端可用性。

图 5. 可可发酵读级与 MAG 级结果。A 与 B 面板显示:读级(A)与 MAG 级(B)分析均表明发酵过程多样性下降,并伴随功能基因谱转变;GT 类 CAZyme 增加、GH 类减少。
方法学参考
#mermaid-svg-glIAH1FxezlWrsUi{font-family:\”trebuchet ms\”,verdana,arial,sans-serif;font-size:16px;fill:#333;}@keyframes edge-animation-frame{from{stroke-dashoffset:0;}}@keyframes dash{to{stroke-dashoffset:0;}}#mermaid-svg-glIAH1FxezlWrsUi .edge-animation-slow{stroke-dasharray:9,5!important;stroke-dashoffset:900;animation:dash 50s linear infinite;stroke-linecap:round;}#mermaid-svg-glIAH1FxezlWrsUi .edge-animation-fast{stroke-dasharray:9,5!important;stroke-dashoffset:900;animation:dash 20s linear infinite;stroke-linecap:round;}#mermaid-svg-glIAH1FxezlWrsUi .error-icon{fill:#552222;}#mermaid-svg-glIAH1FxezlWrsUi .error-text{fill:#552222;stroke:#552222;}#mermaid-svg-glIAH1FxezlWrsUi .edge-thickness-normal{stroke-width:1px;}#mermaid-svg-glIAH1FxezlWrsUi .edge-thickness-thick{stroke-width:3.5px;}#mermaid-svg-glIAH1FxezlWrsUi .edge-pattern-solid{stroke-dasharray:0;}#mermaid-svg-glIAH1FxezlWrsUi .edge-thickness-invisible{stroke-width:0;fill:none;}#mermaid-svg-glIAH1FxezlWrsUi .edge-pattern-dashed{stroke-dasharray:3;}#mermaid-svg-glIAH1FxezlWrsUi .edge-pattern-dotted{stroke-dasharray:2;}#mermaid-svg-glIAH1FxezlWrsUi .marker{fill:#333333;stroke:#333333;}#mermaid-svg-glIAH1FxezlWrsUi .marker.cross{stroke:#333333;}#mermaid-svg-glIAH1FxezlWrsUi svg{font-family:\”trebuchet ms\”,verdana,arial,sans-serif;font-size:16px;}#mermaid-svg-glIAH1FxezlWrsUi p{margin:0;}#mermaid-svg-glIAH1FxezlWrsUi .label{font-family:\”trebuchet ms\”,verdana,arial,sans-serif;color:#333;}#mermaid-svg-glIAH1FxezlWrsUi .cluster-label text{fill:#333;}#mermaid-svg-glIAH1FxezlWrsUi .cluster-label span{color:#333;}#mermaid-svg-glIAH1FxezlWrsUi .cluster-label span p{background-color:transparent;}#mermaid-svg-glIAH1FxezlWrsUi .label text,#mermaid-svg-glIAH1FxezlWrsUi span{fill:#333;color:#333;}#mermaid-svg-glIAH1FxezlWrsUi .node rect,#mermaid-svg-glIAH1FxezlWrsUi .node circle,#mermaid-svg-glIAH1FxezlWrsUi .node ellipse,#mermaid-svg-glIAH1FxezlWrsUi .node polygon,#mermaid-svg-glIAH1FxezlWrsUi .node path{fill:#ECECFF;stroke:#9370DB;stroke-width:1px;}#mermaid-svg-glIAH1FxezlWrsUi .rough-node .label text,#mermaid-svg-glIAH1FxezlWrsUi .node .label text,#mermaid-svg-glIAH1FxezlWrsUi .image-shape .label,#mermaid-svg-glIAH1FxezlWrsUi .icon-shape .label{text-anchor:middle;}#mermaid-svg-glIAH1FxezlWrsUi .node .katex path{fill:#000;stroke:#000;stroke-width:1px;}#mermaid-svg-glIAH1FxezlWrsUi .rough-node .label,#mermaid-svg-glIAH1FxezlWrsUi .node .label,#mermaid-svg-glIAH1FxezlWrsUi .image-shape .label,#mermaid-svg-glIAH1FxezlWrsUi .icon-shape .label{text-align:center;}#mermaid-svg-glIAH1FxezlWrsUi .node.clickable{cursor:pointer;}#mermaid-svg-glIAH1FxezlWrsUi .root .anchor path{fill:#333333!important;stroke-width:0;stroke:#333333;}#mermaid-svg-glIAH1FxezlWrsUi .arrowheadPath{fill:#333333;}#mermaid-svg-glIAH1FxezlWrsUi .edgePath .path{stroke:#333333;stroke-width:2.0px;}#mermaid-svg-glIAH1FxezlWrsUi .flowchart-link{stroke:#333333;fill:none;}#mermaid-svg-glIAH1FxezlWrsUi .edgeLabel{background-color:rgba(232,232,232, 0.8);text-align:center;}#mermaid-svg-glIAH1FxezlWrsUi .edgeLabel p{background-color:rgba(232,232,232, 0.8);}#mermaid-svg-glIAH1FxezlWrsUi .edgeLabel rect{opacity:0.5;background-color:rgba(232,232,232, 0.8);fill:rgba(232,232,232, 0.8);}#mermaid-svg-glIAH1FxezlWrsUi .labelBkg{background-color:rgba(232, 232, 232, 0.5);}#mermaid-svg-glIAH1FxezlWrsUi .cluster rect{fill:#ffffde;stroke:#aaaa33;stroke-width:1px;}#mermaid-svg-glIAH1FxezlWrsUi .cluster text{fill:#333;}#mermaid-svg-glIAH1FxezlWrsUi .cluster span{color:#333;}#mermaid-svg-glIAH1FxezlWrsUi div.mermaidTooltip{position:absolute;text-align:center;max-width:200px;padding:2px;font-family:\”trebuchet ms\”,verdana,arial,sans-serif;font-size:12px;background:hsl(80, 100%, 96.2745098039%);border:1px solid #aaaa33;border-radius:2px;pointer-events:none;z-index:100;}#mermaid-svg-glIAH1FxezlWrsUi .flowchartTitleText{text-anchor:middle;font-size:18px;fill:#333;}#mermaid-svg-glIAH1FxezlWrsUi rect.text{fill:none;stroke-width:0;}#mermaid-svg-glIAH1FxezlWrsUi .icon-shape,#mermaid-svg-glIAH1FxezlWrsUi .image-shape{background-color:rgba(232,232,232, 0.8);text-align:center;}#mermaid-svg-glIAH1FxezlWrsUi .icon-shape p,#mermaid-svg-glIAH1FxezlWrsUi .image-shape p{background-color:rgba(232,232,232, 0.8);padding:2px;}#mermaid-svg-glIAH1FxezlWrsUi .icon-shape .label rect,#mermaid-svg-glIAH1FxezlWrsUi .image-shape .label rect{opacity:0.5;background-color:rgba(232,232,232, 0.8);fill:rgba(232,232,232, 0.8);}#mermaid-svg-glIAH1FxezlWrsUi .label-icon{display:inline-block;height:1em;overflow:visible;vertical-align:-0.125em;}#mermaid-svg-glIAH1FxezlWrsUi .node .label-icon path{fill:currentColor;stroke:revert;stroke-width:revert;}#mermaid-svg-glIAH1FxezlWrsUi :root{–mermaid-font-family:\”trebuchet ms\”,verdana,arial,sans-serif;}
从SRA用q2-fondue取读
q2-fastp质控与去人源读
q2-assembly用MEGAHIT或metaSPAdes组装
Bowtie2建索引并读回贴
q2-mag用MetaBAT2分箱
BUSCO与GUNC做质量与嵌合体控制
q2-sourmash去冗余得非冗余MAG
q2-annotate用Kraken2或Kaiju分类
eggNOG-mapper或HUMANN功能注释
读回贴估MAG丰度并整合功能表
Q2F多样性统计与ANCOM-BC2等下游分析
可复用要点:用 Q2F Artifact 统一中间产物,避免脚本拼凑;用 CAMI II/LEMMI 上榜工具保证方法质量;用 provenance replay 自动生成可重跑文档;用 Parsl 并行在笔记本到 HPC 间扩展;用 sourmash MinHash 做 MAG 去冗余;用读回贴 RPKM/TPM 把 MAG 丰度与功能注释表点乘得到样本级功能丰度。
总结
MOSHPIT 的本质不是“再写一个组装器”,而是把经过独立基准验证的 WMS 工具,规范化为 QIIME 2 生态中的第一类公民。它补上了 QIIME 2 在非靶向宏基因组上的短板:研究者可以用同一套 Artifact、同一套溯源日志、同一套用户界面,完成从原始读、重叠群到 MAG 的分类与功能分析,并直接对接 QIIME 2 已有的多样性与统计模块。
作者也明确其边界:核心版只封装经 CAMI II 验证的成熟算法,长读长原生组装、深度学习分箱器等尚未原生纳入;Q2F action 抽象会略微限制极端参数微调;案例研究用了公开/子集数据而非全生物群穷举基准;.qza 缓存与格式转换会带来磁盘开销。这些是为“严格可重复”让位于“极致灵活”的主动取舍。
未来路线包括新增分箱细化工具、新分类器、机器学习 MAG 质量评估、生物合成基因簇检测,以及真核宏基因组支持。对微生物生态、临床宏基因组、食品发酵与海洋组学用户而言,MOSHPIT 提供的是一条“低门槛入场、高可审计产出”的 WMS 分析路径。
参考文献
- Ziemski, M., Gehret, L., Simard, A., Castro Dau, S., Risch, V., Grabocka, D., Matzzoros, C., Wood, C., Momo Cabrera, P., Hernandez-Velázquez, R., Scharer, M. R., Sandberg, F., Herman, C., Evans, K., Robeson, M. S. II, Bolyen, E., Caporaso, J. G., & Bokulich, N. A. (2026). MOSHPIT: Accessible, reproducible metagenome data science on the QIIME 2 framework. Cell Genomics, 6, 101366. https://doi.org/10.1016/j.xgen.2026.101366
- Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857 (2019).
- Meyer, F. et al. Critical Assessment of Metagenome Interpretation: the second round of challenges. Nat. Methods 19, 429–440 (2022).
- Seppey, M. et al. LEMMI: A continuous benchmarking platform for metagenomics classifiers. Genome Res. 30, 1208–1216 (2020).
- MOSHPIT 文档:https://moshpit.qiime2.org
- q2-assembly:https://github.com/bokulich-lab/q2-assembly
- q2-mag:https://github.com/bokulich-lab/q2-mag
- q2-annotate:https://github.com/bokulich-lab/q2-annotate
- nf-moshpit:https://github.com/bokulich-lab/nf-moshpit
- Cell Genomics 期刊主页:https://www.cell.com/cell-genomics/home
网硕互联帮助中心







评论前必须登录!
注册