1. 单细胞测序文库制备工程导论

单细胞测序(single-cell sequencing)把“组织平均”拆解为“逐个细胞”的分子读数,而文库制备工程正是把细胞悬液转化为可上机测序文库的全部物理、化学与信息学前置工作。它决定了后续数据的上限:细胞捕获率、双细胞率、基因检出数、转录本覆盖均一性、批次可重复性,几乎都在文库构建阶段被锁定。本课程的目标不是罗列试剂盒说明书,而是建立一套“可设计、可计算、可验证、可排障”的工程化思维:给定样本类型与生物学问题,能反推所需的捕获平台、细胞投入量、条码策略、扩增循环数与质控阈值,并能解释每一步参数变化对最终数据的影响。

技术发展脉络:从孔板到液滴

单细胞测序的起点是低通量的手工分选。2009年Tang等用显微操作挑取单个卵裂球,配合末端加尾与体外转录扩增完成转录组测序,通量仅个位数细胞。2012年前后,SMART-seq与CEL-seq引入模板转换(template switching)与体外转录线性扩增,把单细胞转录组推进到百细胞量级,但仍是孔板操作,成本高、周期长。真正的转折来自微流控与条码(barcoding)思想:2015年Drop-seq与inDrop分别用液滴包裹单细胞与凝胶珠,把“一个细胞一个孔”变成“一个细胞一个液滴”,一次运行即可捕获数千至数万细胞。2016年后,10x Genomics的GemCode/Chromium体系把液滴微流控、凝胶珠条码与油相分隔标准化为商业平台,使单细胞转录组成为常规实验。同期,Smart-seq2、Smart-seq3保留了全长读长与高灵敏度优势,用于低丰度转录本、可变剪接与等位基因分析。近年,单细胞多组学(RNA+ATAC、RNA+蛋白、RNA+V(D)J)与空间转录组进一步扩展了文库制备的对象:不仅要做mRNA,还要在同一细胞或同一空间位置同时捕获表观、蛋白与免疫受体信息。

文库制备在单细胞研究中的核心地位

单细胞实验的误差来源可分为三类:样本解离引入的应激与细胞丢失、捕获过程引入的双细胞与空液滴、文库构建引入的扩增偏倚与条码错误。其中后两类直接由文库制备工程控制。以液滴法为例,细胞悬液与凝胶珠、油相在微流控芯片中汇合,理想情况下每个液滴包含一个细胞与一个凝胶珠。凝胶珠上的捕获探针通常由三部分组成:一段用于后续PCR扩增的通用序列、一段细胞条码(cell barcode,标识来源细胞)、一段UMI(unique molecular identifier,标识原始mRNA分子)。细胞裂解后,mRNA的poly(A)尾与探针上的oligo(dT)杂交,反转录生成带条码的cDNA。此后破乳、纯化、PCR扩增、片段化与接头连接,最终形成P5/P7双端测序文库。这里的关键工程参数包括:细胞悬液浓度与活力、凝胶珠与细胞的比例、油相流速、反转录效率、PCR循环数、片段化时间。任何一个参数偏离,都会在数据中留下可识别的痕迹,例如细胞浓度过高导致双细胞率上升,PCR循环过多导致UMI重复与饱和,片段化不足导致插入片段偏长、测序有效数据率下降。

工程视角的核心判断:文库制备不是“按说明书操作”,而是在给定样本与预算下,对捕获率、双细胞率、基因检出灵敏度、成本四者做权衡。例如,目标细胞数1万、双细胞率控制在5%以内时,液滴法的细胞投入浓度通常需根据平台标定的泊松分布反推;若样本本身细胞数极少(如稀有循环肿瘤细胞),则应优先选择低损失的分选与孔板方案,而非追求高通量。

整体工程流程概览

一个完整的单细胞文库制备工程可拆为六个阶段,各阶段之间存在明确的输入输出与质控点。第一阶段是样本与细胞悬液制备:实体组织需酶解或机械解离,血液或培养细胞需洗涤与计数,关键质控是活力>85%、结团率低、无大量碎片。第二阶段是捕获与条码标记:按平台选择液滴、纳米孔或孔板,完成细胞与条码珠的配对。第三阶段是反转录与cDNA扩增:将带条码的mRNA转为cDNA并扩增到可检测量。第四阶段是文库构建:片段化、末端修复、加A、接头连接、索引PCR,引入样本索引以支持多样本混样。第五阶段是文库质控:用Qubit测浓度、Bioanalyzer或TapeStation测插入片段分布、qPCR测可扩增文库比例。第六阶段是上机测序与数据拆分:测序后按细胞条码与UMI拆分,生成表达矩阵。每个阶段都有对应的失败模式与排查路径,例如cDNA产量低可能源于反转录效率差或细胞活力不足,插入片段偏大可能源于片段化时间不足,索引PCR后出现引物二聚体则需调整纯化比例。

单细胞文库制备工程:六阶段流程与质控锚点 ① 样本与悬液 解离 / 洗涤 / 计数 QC:活力>85% ② 捕获与条码 液滴 / 孔板 / 纳米孔 QC:双细胞率 ③ 反转录扩增 RT + PCR / IVT QC:cDNA产量 ④ 文库构建 片段化 / 接头 / 索引 QC:插入片段 ⑤ 文库质控 Qubit / qPCR / 片段 QC:可扩增比例 ⑥ 上机测序与数据拆分 按细胞条码 + UMI 拆分 → 表达矩阵 关键工程参数(示例范围) • 细胞悬液浓度:100–1000 cells/μL(按平台标定) • 目标双细胞率:<5%(泊松分布反推投入量) • 反转录效率:10–30%(决定基因检出数) • PCR循环数:10–14(避免UMI饱和) • 插入片段主峰:400–600 bp(片段化时间控制) • 文库上机浓度:1–4 nM(qPCR定量后稀释) 常见失败模式与排查方向 • cDNA产量低 → 细胞活力 / RT效率 / 条码珠 • 双细胞率高 → 细胞浓度过高 / 结团 • 插入片段偏大 → 片段化时间不足 • 引物二聚体 → 纯化比例 / 循环数过高 • 基因检出少 → UMI饱和 / 测序深度不足 • 批次效应 → 索引分配 / 混样策略

把上述流程放在工程视角下,最重要的不是记住每个试剂盒的步骤,而是理解“输入—操作—中间变化—输出—判断”的闭环。以液滴法为例,输入是细胞悬液与凝胶珠,操作是微流控共包裹与反转录,中间变化是液滴内mRNA被条码化并转为cDNA,输出是带条码的cDNA文库,判断依据是cDNA产量、插入片段分布与测序后的细胞条码有效性。若cDNA产量低于预期,应先回看细胞活力与计数是否准确,再检查反转录体系与条码珠批次;若插入片段主峰偏大,应缩短片段化时间或调整酶量;若测序后有效细胞条码比例低,应检查油相流速与破乳效率。这些判断不需要等到测序结束才做,在文库质控阶段就能通过片段分布与qPCR结果提前预警。

学习建议:本课程后续章节会逐一展开每个阶段的原理、参数与排障方法。建议在学习时始终带着三个问题:这一步的输入是什么?关键参数如何影响输出?出现异常时,最早能在哪个质控点发现?把这三个问题贯穿30章,就能把零散的实验步骤串成可复用的工程能力。

最后需要强调,单细胞文库制备工程不是孤立的技术操作,而是与生物学问题、测序策略、数据分析深度紧密耦合的系统工程。选择全长还是3'端、选择液滴还是孔板、选择单组学还是多组学,都会反向决定文库制备的路线。课程后续将围绕这些决策点,给出可计算、可验证、可复现的工程方法。