查看“Mutect2”的源代码

<div style="padding: 0 4%; line-height: 1.8; color: #1e293b; font-family: 'Helvetica Neue', Helvetica, 'PingFang SC', Arial, sans-serif; background-color: #ffffff; max-width: 1200px; margin: auto;">

    <div style="margin-bottom: 30px; border-bottom: 1.2px solid #e2e8f0; padding-bottom: 25px;">
        <p style="font-size: 1.1em; margin: 10px 0; color: #334155; text-align: justify;">
            <strong>Mutect2</strong> 是由 <strong>Broad Institute</strong> 开发的一款用于探测<strong>体细胞突变</strong>（Somatic Mutations）的生物信息学工具，属于 <strong>[[GATK]]</strong>（Genome Analysis Toolkit）软件包的核心组件。
            <br>它在第一代 MuTect（仅支持 SNV）的基础上进行了重大升级，能够同时检测<strong>单核苷酸变异</strong>（[[SNV]]）和<strong>插入缺失</strong>（[[Indel]]）。
            <br>Mutect2 采用<strong>贝叶斯分类器</strong>和<strong>单倍体组装</strong>（Haplotype Assembly）技术，通过对比“肿瘤-正常组织”配对样本（Tumor-Normal Pair），有效剔除生殖系变异和测序噪音，是目前癌症基因组分析（如 [[WES]]、[[WGS]]）中检测体细胞突变的“金标准”工具之一。
        </p>
    </div>

    <div class="medical-infobox mw-collapsible mw-collapsed" style="width: 100%; max-width: 320px; margin: 0 auto 35px auto; border: 1.2px solid #bae6fd; border-radius: 12px; background-color: #ffffff; box-shadow: 0 8px 20px rgba(0,0,0,0.05); overflow: hidden;">
        
        <div style="padding: 15px; color: #1e40af; background: linear-gradient(135deg, #e0f2fe 0%, #bae6fd 100%); text-align: center; cursor: pointer;">
            <div style="font-size: 1.2em; font-weight: bold; letter-spacing: 1.2px;">Mutect2</div>
            <div style="font-size: 0.7em; opacity: 0.85; margin-top: 4px; white-space: nowrap;">Somatic Variant Caller (点击展开)</div>
        </div>
        
        <div class="mw-collapsible-content">
            <div style="padding: 25px; text-align: center; background-color: #f8fafc;">
                 
                <div style="font-size: 0.8em; color: #64748b; margin-top: 12px; font-weight: 600;">癌症突变的“捕手”</div>
            </div>

            <table style="width: 100%; border-spacing: 0; border-collapse: collapse; font-size: 0.85em;">
                <tr>
                    <th colspan="2" style="padding: 8px 12px; background-color: #e0f2fe; color: #1e40af; text-align: left; font-size: 0.9em; border-top: 1px solid #bae6fd;">软件档案</th>
                </tr>
                <tr>
                    <th style="text-align: left; padding: 6px 12px; background-color: #f8fafc; color: #475569; border-bottom: 1px solid #e2e8f0; width: 40%;">开发者</th>
                    <td style="padding: 6px 12px; border-bottom: 1px solid #e2e8f0; color: #0f172a;">Broad Institute</td>
                </tr>
                <tr>
                    <th style="text-align: left; padding: 6px 12px; background-color: #f8fafc; color: #475569; border-bottom: 1px solid #e2e8f0;">所属套件</th>
                    <td style="padding: 6px 12px; border-bottom: 1px solid #e2e8f0; color: #1e40af;">[[GATK]] 4.x</td>
                </tr>
                <tr>
                    <th style="text-align: left; padding: 6px 12px; background-color: #f8fafc; color: #475569; border-bottom: 1px solid #e2e8f0;">检测对象</th>
                    <td style="padding: 6px 12px; border-bottom: 1px solid #e2e8f0; color: #0f172a;">体细胞 SNVs & Indels</td>
                </tr>
                <tr>
                    <th style="text-align: left; padding: 6px 12px; background-color: #f8fafc; color: #475569; border-bottom: 1px solid #e2e8f0;">输入格式</th>
                    <td style="padding: 6px 12px; border-bottom: 1px solid #e2e8f0; color: #0f172a;">[[BAM]] / CRAM</td>
                </tr>
                
                <tr>
                    <th colspan="2" style="padding: 8px 12px; background-color: #e0f2fe; color: #1e40af; text-align: left; font-size: 0.9em; border-top: 1px solid #bae6fd;">核心算法</th>
                </tr>
                <tr>
                    <th style="text-align: left; padding: 6px 12px; background-color: #f8fafc; color: #475569; border-bottom: 1px solid #e2e8f0;">组装策略</th>
                    <td style="padding: 6px 12px; border-bottom: 1px solid #e2e8f0; color: #1e40af;">局部单倍体组装</td>
                </tr>
                <tr>
                    <th style="text-align: left; padding: 6px 12px; background-color: #f8fafc; color: #475569; border-bottom: 1px solid #e2e8f0;">噪音过滤</th>
                    <td style="padding: 6px 12px; border-bottom: 1px solid #e2e8f0; color: #0f172a;">[[PON]] (Panel of Normals)</td>
                </tr>
                <tr>
                    <th style="text-align: left; padding: 6px 12px; background-color: #f8fafc; color: #475569; border-bottom: 1px solid #e2e8f0;">生殖系扣除</th>
                    <td style="padding: 6px 12px; border-bottom: 1px solid #e2e8f0; color: #0f172a;">配对 Normal + [[gnomAD]]</td>
                </tr>
                 <tr>
                    <th style="text-align: left; padding: 6px 12px; background-color: #f8fafc; color: #475569;">输出格式</th>
                    <td style="padding: 6px 12px; color: #e11d48;">[[VCF]] (需 FilterMutectCalls)</td>
                </tr>
            </table>
        </div>
    </div>

    <h2 style="background: #f1f5f9; color: #0f172a; padding: 10px 18px; border-radius: 0 6px 6px 0; font-size: 1.25em; margin-top: 40px; border-left: 6px solid #0f172a; font-weight: bold;">从 MuTect 到 Mutect2 的进化</h2>
    <p style="margin: 15px 0; text-align: justify;">
        Mutect2 并非简单的升级，而是算法逻辑的重构。它继承了 <strong>[[HaplotypeCaller]]</strong> 的“活性区域组装”（Active Region Assembly）引擎，但针对体细胞突变的特点进行了调整。
    </p>
    <div style="overflow-x: auto; margin: 20px auto;">
        <table style="width: 100%; border-collapse: collapse; border: 1.2px solid #cbd5e1; font-size: 0.9em; text-align: left;">
            <tr style="background-color: #f1f5f9; border-bottom: 2px solid #0f172a;">
                <th style="padding: 12px; border: 1px solid #cbd5e1; color: #0f172a; width: 25%;">特性</th>
                <th style="padding: 12px; border: 1px solid #cbd5e1; color: #1e40af; width: 35%;">MuTect (v1)</th>
                <th style="padding: 12px; border: 1px solid #cbd5e1; color: #475569; width: 40%;">Mutect2 (GATK4)</th>
            </tr>
            <tr>
                <td style="padding: 10px; border: 1px solid #cbd5e1; font-weight: 600;">变异类型</td>
                <td style="padding: 10px; border: 1px solid #cbd5e1;">仅 <strong>SNV</strong></td>
                <td style="padding: 10px; border: 1px solid #cbd5e1;"><strong>SNV + Indel</strong></td>
            </tr>
            <tr>
                <td style="padding: 10px; border: 1px solid #cbd5e1; font-weight: 600;">倍体假设</td>
                <td style="padding: 10px; border: 1px solid #cbd5e1;">二倍体 (Diploid)</td>
                <td style="padding: 10px; border: 1px solid #cbd5e1;"><strong>可变倍体</strong> (适应肿瘤非整倍性)</td>
            </tr>
            <tr>
                <td style="padding: 10px; border: 1px solid #cbd5e1; font-weight: 600;">分析模式</td>
                <td style="padding: 10px; border: 1px solid #cbd5e1;">必须成对 (Tumor-Normal)</td>
                <td style="padding: 10px; border: 1px solid #cbd5e1;">支持成对，也支持 <strong>Tumor-only</strong> 模式</td>
            </tr>
             <tr>
                <td style="padding: 10px; border: 1px solid #cbd5e1; font-weight: 600;">线粒体</td>
                <td style="padding: 10px; border: 1px solid #cbd5e1;">不支持</td>
                <td style="padding: 10px; border: 1px solid #cbd5e1;">支持 (Mitochondria mode)</td>
            </tr>
        </table>
    </div>

    <h2 style="background: #f1f5f9; color: #0f172a; padding: 10px 18px; border-radius: 0 6px 6px 0; font-size: 1.25em; margin-top: 40px; border-left: 6px solid #0f172a; font-weight: bold;">标准分析流程 (Best Practices)</h2>
    <p style="margin: 15px 0; text-align: justify;">
        GATK 团队推荐的体细胞突变分析流程包含极其严格的过滤步骤，以去除假阳性。
    </p>
    <div style="background-color: #f0f9ff; border-left: 5px solid #1e40af; padding: 15px 20px; margin: 20px 0; border-radius: 4px;">
        <ul style="margin: 0; padding-left: 20px; color: #334155;">
            <li style="margin-bottom: 12px;"><strong>Step 1: Create PON (构建正常样本库)：</strong> 使用一组（推荐 > 40个）健康人的样本生成 <strong>[[PON]]</strong> (Panel of Normals)。这是 Mutect2 的核心优势，用于捕捉测序仪特有的系统性错误和罕见生殖系变异。</li>
            <li style="margin-bottom: 12px;"><strong>Step 2: Mutect2 Calling：</strong> 输入肿瘤 BAM、配对正常 BAM（可选）、PON 和生殖系资源库（如 [[gnomAD]]）。Mutect2 会输出一个原始的 VCF 文件。</li>
            <li style="margin-bottom: 12px;"><strong>Step 3: Estimate Contamination：</strong> 使用 <code>GetPileupSummaries</code> 和 <code>CalculateContamination</code> 计算样本间的交叉污染率。</li>
            <li style="margin-bottom: 0;"><strong>Step 4: FilterMutectCalls：</strong> 综合考虑污染率、方向偏差（Orientation Bias, 如 FFPE 样本的 OxoG 损伤）和 PON 数据，给每个突变打上 <code>PASS</code> 或过滤标签。</li>
        </ul>
    </div>

    <div style="font-size: 0.92em; line-height: 1.6; color: #1e293b; margin-top: 50px; border-top: 2px solid #0f172a; padding: 15px 25px; background-color: #f8fafc; border-radius: 0 0 10px 10px;">
        <span style="color: #0f172a; font-weight: bold; font-size: 1.05em; display: inline-block; margin-bottom: 15px;">关键相关概念 [Key Concepts]</span>
        
        <p style="margin: 12px 0; border-bottom: 1px solid #e2e8f0; padding-bottom: 10px;">
            <strong>1. PON (Panel of Normals)：</strong> Mutect2 的“黑名单”。如果一个突变在多个正常样本中都出现了，说明它很可能是测序错误（Technical Artifact）而非真正的体细胞突变，Mutect2 会据此将其过滤。
        </p>

        <p style="margin: 12px 0; border-bottom: 1px solid #e2e8f0; padding-bottom: 10px;">
            <strong>2. AF (Allele Fraction)：</strong> 突变等位基因频率。与生殖系突变通常为 0.5 (杂合) 或 1.0 (纯合) 不同，体细胞突变的 AF 波动巨大（可能低至 1%），这取决于肿瘤纯度（Purity）和异质性（Heterogeneity）。Mutect2 专为检测低 AF 变异而优化。
        </p>

        <p style="margin: 12px 0;">
            <strong>3. Germline Resource：</strong> 如 [[gnomAD]] 或 [[1000 Genomes]]。Mutect2 利用这些公共数据库来估算一个变异是生殖系的可能性。如果患者没有配对的正常样本（Tumor-only），这一步至关重要。
        </p>
    </div>

    <div style="font-size: 0.92em; line-height: 1.6; color: #1e293b; margin-top: 20px; border-top: 2px solid #0f172a; padding: 15px 25px; background-color: #ffffff;">
        <span style="color: #0f172a; font-weight: bold; font-size: 1.05em; display: inline-block; margin-bottom: 15px;">学术参考文献 [Academic Review]</span>
        
        <p style="margin: 12px 0; border-bottom: 1px solid #e2e8f0; padding-bottom: 10px;">
            [1] <strong>Cibulskis K, Lawrence MS, Carter SL, et al. (2013).</strong> <em>Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples.</em> <strong>[[Nature Biotechnology]]</strong>. <br>
            <span style="color: #475569;">[点评]：初代 MuTect 的奠基论文。详细阐述了基于贝叶斯统计的体细胞突变检测原理，成为该领域的引用经典。</span>
        </p>

        <p style="margin: 12px 0; border-bottom: 1px solid #e2e8f0; padding-bottom: 10px;">
            [2] <strong>Benjamin D, Sato T, Cibulskis K, et al. (2019).</strong> <em>Calling Somatic SNVs and Indels with Mutect2.</em> <strong>[[bioRxiv]]</strong>. <br>
            <span style="color: #475569;">[点评]：GATK 团队发布的关于 Mutect2 具体算法的技术预印本。解释了从 MuTect 到 Mutect2 的算法演进，特别是对 Indel 的支持和 PON 的改进。</span>
        </p>
    </div>

    <div style="margin: 40px 0; border: 1px solid #e2e8f0; border-radius: 8px; overflow: hidden; font-family: 'Helvetica Neue', Arial, sans-serif; font-size: 0.9em;">
        <div style="background-color: #eff6ff; color: #1e40af; padding: 8px 15px; font-weight: bold; text-align: center; border-bottom: 1px solid #dbeafe;">
            生物信息学 · 知识图谱
        </div>
        <table style="width: 100%; border-collapse: collapse; background-color: #ffffff;">
            <tr style="border-bottom: 1px solid #f1f5f9;">
                <td style="width: 85px; background-color: #f8fafc; color: #334155; font-weight: 600; padding: 10px 12px; text-align: right; vertical-align: middle;">上级分类</td>
                <td style="padding: 10px 15px; color: #334155;">[[变异检测]] (Variant Calling) • 癌症基因组学</td>
            </tr>
            <tr style="border-bottom: 1px solid #f1f5f9;">
                <td style="width: 85px; background-color: #f8fafc; color: #334155; font-weight: 600; padding: 10px 12px; text-align: right; vertical-align: middle;">输入数据</td>
                <td style="padding: 10px 15px; color: #334155;">[[BAM]] • [[PON]] • 参考基因组</td>
            </tr>
            <tr>
                <td style="width: 85px; background-color: #f8fafc; color: #334155; font-weight: 600; padding: 10px 12px; text-align: right; vertical-align: middle;">竞品软件</td>
                <td style="padding: 10px 15px; color: #334155;">[[Strelka2]] • [[VarScan2]] • [[MuSE]]</td>
            </tr>
        </table>
    </div>

</div>