Method Article
单细胞测序揭示了生物系统中的基因型异质性, 但目前的技术缺乏对群落组成和功能的深入分析所必需的吞吐量。在这里, 我们描述了一个微流控流程的序列 > 5万单细胞基因组从不同的细胞群。
测序技术经历了从大批量到单细胞分辨率的范式转变, 以响应对生物系统中细胞异质性作用的逐渐理解。然而, 由于处理基因组测序的局限性, 大量人群的单细胞测序受到限制。本文描述了单细胞基因组测序 (SiC 序列) 的方法, 它利用液滴微流体分离、放大和条形码的单细胞基因组。微中的细胞封装允许分离纯化和 tagmentation DNA, 而微流控合并有效地将每个基因组与一个独特的单细胞寡核苷酸条码配对, 允许 > 5万单胞细胞被测序每次运行。测序数据由条形码 demultiplexed, 生成来自单个单元格的读取组。作为单细胞测序的高通量和低偏置方法, SiC 序列将使针对不同细胞群的更广泛的基因组研究得以实现。
基因组作为细胞身份和功能的蓝图, 包含了有机体的编码潜能的整体。在基因组层面上对细胞生物学的理解可以解释异质细胞群中观察到的表型多样性。这种异质性在生物系统中很明显, 对人类健康和疾病有着广泛的影响。例如, 肿瘤细胞中基因复制数的变化与癌症的演变和扩散有关, 如1,2。在细菌感染中, 在一小部分基因组中存在的致病性岛屿可以水平转移, 并导致抗药性细菌的扩散3,4。在单细胞水平上研究基因组的一个主要挑战是可用的 DNA 数量少, 以及需要分析数以千计的细胞来取样基因型的全部多样性。由于这些原因, 实验吞吐量的限制阻碍了单细胞研究的有效性, 使结果偏向于最丰富的细胞。单单元隔离技术, 如流排序5,6, 光学镊子7, 在大容量凝胶8中嵌入, 微流体9能够处理数以百计的细胞进行测序;但是, 这仅代表大多数示例中的一小部分。单细胞基因组测序的方法具有极高的吞吐量, 可以更深入和更完整地分析细胞种群, 从而阐明基因型多样性在这些社区中的作用。
液滴微流体能够在数以百万计的 picoliter 反应器内对细胞和生物试剂进行高通量的操作。到目前为止, microdroplet 技术已经被用来研究不同组织的细胞之间的差异表达模式10,11,12, 深序列长分子13,14 、15, 并对单个单元格16进行染色质免疫沉淀排序 (芯片序列) 分析。事实上, microdroplets 能够进行高通量的、分割的操作, 使它们能够在单细胞基因组学中应用。然而, 这种技术的发展呈现出自己独特的技术挑战。细胞必须裂解, 纯化, 并以最小的偏差放大, 以均匀的样本细胞数量。另外, 与哺乳动物细胞中的 polyadenylated mRNA 记录不同, 基因组中没有可比的分子母题来促进目标核酸的捕获。由于这些原因, 单细胞基因组测序在 microdroplet 平台中难以实现。
在这项工作中, 我们提供了我们以前报告的单细胞微流控方法的详细协议, 能够在单个实验17中对成千上万个细胞的基因组进行排序。使用这种技术, 称为 SiC 序列, 细菌细胞被封装在微米尺度的水凝胶和单独裂解, tagmented, 并合并 microdroplet 含有一个独特的寡核苷酸条码, 这是拼接到细胞的基因组 DNA 通过单重叠扩展聚合酶链反应 (PCR)。水凝胶作为孤立的容器, 其中高分子量的基因组 DNA 被 sterically 包裹, 允许较小的分子, 如洗涤剂和裂解酶, 以获取和纯化 DNA 之前, 条形码18。此协议在几个小时内处理 > 5万个单元格, 导致条形码库准备好进行排序。按照顺序, 读取是根据他们的单单元条码序列 demultiplexed 的, 结果是由数以百万计的读取组成的数据集, 每个都带有一个蜂窝索引。
1. 微流控器件制造
2. 琼脂糖微中细胞的封装
注意: 请参阅图 2A。
3. 打破和洗涤琼脂糖微
4. 通过裂解酶裂解琼脂糖中的细胞
5. 基于洗涤剂的微凝胶体处理
6. 用数字 PCR 法生成条码滴
7. Tagmentation 基因组 DNA 在液滴中的提取
注意: 请参阅图 2B。
8. 单细胞条形码微流控双合并
注意: 请参阅图 2C。
9. 图书馆的准备和排序
10. 单细胞数据分析
注意: 定制 Python 脚本的质量控制和初步分析的 SiC 序列数据可以从 https://www.github.com/AbateLab/SiC-seq 下载。
SiC 序列实验工作流包含3个用软光刻过程制作的微流控芯片 (图 1)。共流 dropmaker (图 3A) 生成25µm 的数字条码水滴, 用于标记基因组 DNA, 并使用唯一的单细胞标识符。条形码寡核苷酸由 15 bp 退化序列, 由 PCR 手柄两侧的放大 (表 2, 棒底漆)。条码被稀释到 femtomolar 浓度, 以实现单分子封装, 所有的水滴接收0或1条码片段 (s)。含有条形码的水滴放大, 产生了许多双绞型条码 amplicons 的拷贝。核酸染色用于验证成功的放大和量化条形码片段的封装速率 (图 4B)。微是由细菌细胞悬浮液和熔融琼脂糖凝胶在同等流速下产生的 (图 2A)。琼脂糖是准备在两次预期的最终浓度, 因为共流 dropmaking 过程有效稀释水溶液的系数为2。当琼脂糖冷却, 它固化成一个25µm 直径微凝胶体占据球形体积的水滴。
一系列的洗涤和裂解步骤净化微中的高分子量基因组 DNA (图 2B)。在破乳剂后, 在大体积内进行水洗, 稀释微量有机溶剂, 抑制下游酶处理。在显微镜下观察洗涤的微, 以验证细胞封装速率 (图 4A)。在微凝胶体悬浮液中加入含有广泛溶解活性的酶的鸡尾酒, 以消化细菌和真核微生物的细胞壁19。第二种蛋白酶 K 和洗涤剂的治疗会降低蛋白质和 solubilizes 细胞的碎片。
纯化 dna 的 Tagmentation 是在水滴中进行的, 以避免在微18之间扩散小 tagmented dna 片段导致潜在的交叉污染。液滴封装设备 (图 3B) 划分每个微凝胶体与一个缓冲和 tagmentation 酶, 它同时碎片双链 DNA, 同时也 "标记" 它与预先加载的寡核苷酸20。微被装入水滴作为接近包装的微粒, 达到封装率接近1微凝胶体为每下落与少量双峰21(图 4C)。
在微流控流的最后一步 (图 2C) 中, 一个设备执行双重合并操作, 组合1个条码降、1个含微凝胶体的降, 以及在受控的两步过程中的放大组合。首先, 含有 PCR 试剂的液滴被配对, 并与黄色 (图 5) 中显示的条形码下落合并。微流控通道中的盐水电极产生高电场梯度, 触发熔滴的合并。以类似的方式, 第一个合并的水滴与微凝胶体液滴配对, 并在以红色显示的区域中第二次合并。在单重叠扩展 (SOE) PCR 中采集液滴和热循环脱离芯片。条形码和 tagmented 基因组 DNA 重叠的互补端允许只有适当的条形码结构的融合和指数放大。
顺序数据首先通过读取质量进行筛选, 然后根据它们的 15 bp 单单元条码序列对读取进行分组分析。要认为条形码组是有效的, 应包含最少的读数数;这种阈值限制了对具有大量排序数据的单元格的分析, 并从数据集中删除了 PCR 变异的条形码 "孤儿"。在这个示例运行中, 最小值被设置为每组 7.5 kbps (每个50读取150个 bp)。条形码计数与组大小的直方图显示, 有效条码组的重要部分刚好高于阈值大小 (图 6A)。
在一个已知微生物群落组成的控制实验中, 用纯度和相对丰度指标来评价 SiC 序列的质量。在这里, 分析了由3革兰氏阴性菌、5革兰氏阳性菌和2株酵母组成的合成10细胞群落。给定条码组的纯度定义为组中最常见的基因组的读数映射数除以组中读取总数的个数。绝大多数条码组的纯度大于 0.95 (图 6B)。单元格类型的相对丰度是通过计算原始读取数和通过计数条形码组来计算的, 其中, 组分配了与其成员读取的一致性相对应的单元格类型 (图 6C)。大量的读取和条码组跟踪大致相同的比例, 表明细胞的人口正在取样, 使某些物种不包含在不成比例的小或大型条码组。从单个物种中绘制所有条码组的总覆盖范围, 表明整个基因组的覆盖率很高, 且很少或没有辍学区域 (图 7)。覆盖范围的均匀性可以通过规范化覆盖率值的频率分布进行验证, 大多数值围绕平均值 (图 7, 嵌入)。
图 1: 通过光刻工艺制备微流控器件.(A) 具有单个特征高度的主模具由自旋涂覆一层 SU-8 光刻胶到硅片上制成。然后用光刻掩模和 UV 光进行图案化, 将暴露的 SU-8 交联。最后, uncrosslinked SU-8 溶解在显影液中。所产生的模具用于铸造与玻璃滑块粘合, 以产生完整的微流控装置。(B) 对于双层设备, 制造过程类似于从自旋涂层和曝光步骤开始。然后重复这些步骤以创建两层设备。请单击此处查看此图的较大版本.
图 2:SiC 序列工作流概述.(A) 微生物悬浮液与 dropmaker 装置中的熔融琼脂糖共同流动, 以封装微中的单个细胞。(B) 微受到一系列清洗, 以净化细菌基因组 DNA。溶解酶消化革兰阳性菌和酵母细胞壁, 洗涤剂 solubilizes 细胞碎片。微被重新封装成水滴, tagmentation 以减少交叉污染。(C) 微流控合并结合了数字 PCR 条码、tagmented 微凝胶体基因组, 以及一个速率 > 1 赫的放大组合。非芯片的国企-PCR 拼接一个独特的单细胞条码到 tagmented 基因组, 并有选择地放大完全条形码结构。请单击此处查看此图的较大版本.
图 3:用于 dropmaking 和微凝胶体重新封装的微流控设备。(A) 此面板显示一个共流 dropmaker (25 µm 的特征高度)。细胞和熔融琼脂糖被引入到设备中, 在相同的流速, 以产生25µm 液滴在25µm x 25 µm 结。对于数字条码 dropmaking, 细胞入口堵塞, 并引入 PCR 组合进入琼脂糖入口。(B) 此面板显示微凝胶体重新封装设备 (25 µm 的特征高度)。微流入漏斗形入口, 以维持其紧密包装的订购, 并在25µm x 30 µm 路口重新封装之前接收 tagmentation 混合体。请单击此处查看此图的较大版本.
图 4: 水滴和微的显微照片.(A) 此面板显示在酶裂解之前洗涤了25µm 微。细菌被荧光染色, 以定量的封装速率。泊松负荷统计数据规定, 单元格的封装速率为 1, 在10滴或更少的情况下, 以最小化多封装事件的频率。(B) 此面板显示了一个荧光显微图像, 25 µm 数字条码水滴处理的核酸染色。含有放大条码片段的水滴产生强烈的荧光信号。(C) 此面板显示微重新封装在50µm 滴中。微的封闭包装使封装速率接近1凝胶每滴, 很少双峰。请单击此处查看此图的较大版本.
图 5: 单细胞基因组条形码的微流控双合并装置.两步合并操作对条形码水滴与 tagmented 基因组在高吞吐量。首先生成一滴 PCR 混合物, 并将其与在黄色中使用盐水电极显示的条形码水滴合并。接下来, 一个包含微凝胶体的水滴在红色显示的区域中第二次被引入并合并。油入口允许精确控制重新注入水滴之间的间距。条形码回注室及其间隔油放置在较短的25µm 层上, 以蓝色着色。所有其他设备功能属于较厚的层, 总高度为45µm。请单击此处查看此图的较大版本.
图 6:10 细胞合成微生物群落的条码组度量.(A) 此面板显示条码组大小的分布。当组大小增加时, 给定大小的组数呈指数递减。每个组的最小阈值为 7.5 kbps 限制对具有足够信息量的组进行分析, 并删除 PCR 突变序列 "孤儿"。(B) 此面板显示条码组纯度的分布。绝大多数 (> 90%) 的群体是非常高纯度 (> 95%)。(C) 此面板显示了在读取和条码组级别计算的10种的相对丰度。2计数方法产生相似的结果, 表明条码组的大小是一致的不同物种。请单击此处查看此图的较大版本.
图 7: 聚合基因组覆盖率枯草芽孢杆菌条码组.将所有条码组的读数从映射到细菌枯草杆菌(N = 9398) 进行汇总并进行分析。圆形覆盖图说明了 SiC 序列读数的覆盖均匀性, 没有可观察到的辍学区域。绕圆周的虚线表示平均覆盖范围 (5.55x)。相对覆盖频率的嵌入直方图表明, 大部分的碱基覆盖在接近全基因组平均的深度, 由虚线表示。请单击此处查看此图的较大版本.
装置 | 第一层高度 (µm) | 第一层旋转速度 (rpm) | 第二层高度 (µm) | 第二层旋转速度 (rpm) |
共流 dropmaker | 25 | 4000 | N/A | N/A |
凝胶再 encapsulator | 25 | 4000 | N/A | N/A |
双重合并 | 25 | 4000 | 20 | 5000 |
表 1: 微流控设备制造参数.本表显示了在 SiC 序列工作流中使用的微流控器件的列表, 其所需的光刻胶自旋涂层的速度 (基于 SU-8 3025 的制造商规格)。
标签 | 序列 (5 ' > 3 ') | ||||
酒吧 | GCAGCTGGCGTAATAGCGAGTACAATCTGCTCTGATGCCGCATAGNNNNNNNNNNNNNNNTAAGCCAGCCCCGACACT | ||||
DNA_BAR | CTGTCTCTTATACACATCTCCGAGCCCACGAGACGTGTCGGGGCTGGCTTA | ||||
P7_BAR | CAAGCAGAAGACGGCATACGAGATCAGCTGGCGTAATAGCG | ||||
P5_DNA | AATGATACGGCGACCACCGAGATCTACACTCGTCGGCAGCGTC | ||||
I7_READ | GCCCACGAGACGTGTCGGGGCTGGCTTA |
表 2: 底漆序列。
补充文件 1:请单击此处下载此文件.
补充文件 2:请单击此处下载此文件.
SiC 序列微流控流产生了数以千计的细菌细胞的单细胞基因组测序数据。将数字条形码拼接到微凝胶体封装的单元格的基因组上, 可以将数据的在硅片反褶积转化为来自同一单元格的条形码读取组。一个已知成分的微生物群落的控制实验是评价条形码组纯度的必要方法。大量低纯度组表明, 细胞封装率过高, 或有明显的水滴交叉污染发生在微流控处理步骤。根据泊松统计数据, 条形码和单元格应封装在1粒子的目标比, 每10滴, 以限制多个封装事件的速率小于5% 的所有非空滴。比这更高的封装速率会以指数的形式增加双峰的速率, 因此在 dropmaking 过程中对封装率的验证至关重要。用户在单个微凝胶体中应特别小心封装多个单元格, 因为共享同一条形码序列的不同单元格的读取不能 bioinformatically 分隔。如果1单元接收2不同的条形码, 条码组的纯度不受影响, 虽然丰度指标是在按条形码顺序计数时倾斜的。
由于合并条件不理想, 液滴交叉污染也可能出现。在成功的操作中, 微流控合并装置 (图 5) 可以 controllably 1 个条码液滴与1微凝胶体和一定量 PCR 试剂。非理想的流速会导致液滴配对不正确的比率: 例如, 1 条码可以与2微配对。协议中列出的所有流量都是估计的, 可能需要根据设备几何和雾滴大小的细微变化来调整。使用高速录音功能 (> 1万帧/秒) 的摄像头的用户应在微流控操作的开始和结束过程中验证正确的水滴合并。用户如果不使用高速相机, 就可以收集少量的合并输出, 并手动测量显微镜下的水滴大小。液滴尺寸应均匀: 过量的并入条形码或微凝胶体下降表明, 应相应减少回注率。
在处理微和 microdroplets 时, 应采取一些一般的预防措施, 以维护其完整性。微, 虽然机械健壮, 必须充分冷却之前, 打破和洗涤步骤, 以确保完整的凝胶。非球形微表明琼脂糖没有足够的时间凝固。清洗微时, 按所需的速度旋转悬浮液以避免产品损失。琼脂糖凝胶具有与水紧密匹配的折射率, 在管22中可能很难看到, 因此用户在吸入前应仔细识别凝胶液的边界。在实验室手套和油管上, 静态力23的累积使水中油滴容易发生聚结。因此, 我们建议在泵浦启动前, 用徒手装入滴注注射器, 并用防静电枪处理所有回灌线。通过缓慢旋转注射器中的乳液, 并手动吸出较大的水滴, 在顶部积聚, 由于其较大的浮力, 可以去除大的合并液滴。
SiC 序列是第一个证明单细胞基因组测序 > 5万细菌细胞的技术。该平台在现有方法的吞吐量方面具有显著优势, 并能够更深入地取样异构微生物群落。迄今为止, 单细胞基因组测序的微流控技术已使用 microchambers9和 microwells24进行细胞隔离和扩增, 但仅有数以上百的细胞的吞吐量。将单个单元格的流排序为 wellplates5,6不需要专门的微流控工具, 但具有类似的低吞吐量。鉴于来自环境的土壤和水样通常在物种级25、26的 > 1000 中具有 alpha 多样性, 因此, 由于它能够取样大量的生物体, 所以 SiC 序列是非常有利的。SiC 序列工作流适用于实验室文化、自然环境或活体宿主的细胞输入。细胞样品只需要在水悬浮和无大颗粒 (> 10 µm) 适合微流控封装。例如, 以前已将该方法应用于海水样本, 使用一系列洗涤和过滤步骤在封装17之前对单元格进行预处理。
SiC 序列协议从每个单元中生成相对稀疏的排序数据量, 可能不适用于所有应用程序。一些生物信息学算法, 如 "从头" 基因组组装或单核苷酸变体 (SNV) 调用需要更高的覆盖率深度才能有效地工作。相反, 可以通过分类 binning 方法27将条码组聚集在硅片中, 以便可以将算法应用于较大的读取集。碳化硅序列工作流的总体条形码效率相对较低, 在输入样本的可用性较低的情况下也可能带来挑战。SiC 序列依赖于泊松分布的条形码封装步骤, 因此大约10% 的细胞接收到一个分子条码, 并在最后的图书馆准备步骤中被放大。虽然这与其他基于 microdroplet 的条形码方案10类似, 但使用珍贵细胞样本的用户可能难以达到适当的库产量, 因此可能需要在最终的时间内增加 PCR 周期的数量。放大步骤。另一种潜在的解决方案为使用微流控技术的用户是在数字 PCR 步骤后排序正条码滴, 从而使整体条形码效率 > 85%28。
SiC 序列技术的一个潜在的未来方向是调整工作流, 用于哺乳动物细胞, 为新的临床单细胞研究铺平道路。例如, 对单癌细胞复制数变异的分析, 可以进一步了解异质性在肿瘤病理学中的作用2。或者, 将 SiC 序列与现有的方法相结合, 以探测和丰富感兴趣的 DNA 顺序29将使目标单细胞序列的亚群或罕见的细胞株进行排序。通过环境样本, 已知的代谢通路中的基因可以被定位, 并与相邻的基因进行分析, 以确定新的基因组岛。从人类宿主环境中, 低效的致病细菌样本可以被分离和序列在单细胞水平, 以更密切地检查他们的基因型的致病力的起源。
与此工作流程有关的专利可能被授权给特派团生物, 其中 Adam r 减弱为股东。
这项工作得到了国家科学基金会通过职业奖 (赠款号 DBI-1253293) 的支持;国立卫生研究院 (NIH) (赠款编号 HG007233-01、R01-EB019453-01、1R21HG007233、DP2-AR068129-01、R01-HG008978);和国防高级研究项目机构生活铸造项目 (合同编号 HR0011-12-C-0065, N66001-12-C-4211, HR0011-12-C-0066)。
Name | Company | Catalog Number | Comments |
3" silicon wafers, P type, virgin test grade | University Wafers | 447 | |
SU-8 3025 photoresist | Microchem | 17030192 | |
Spin coater | Specialty Coating Systems | G3P-8 | |
Photomasks | CadArt Servcies | (custom) | See Supplemental Files for mask designs |
PGMEA developer | Sigma-Aldrich | 484431 | |
Isopropanol | Sigma-Aldrich | 109827 | |
Sylgard 184 silicone elastomer kit | Krayden | 4019862 | |
Degassing chamber | Bel-Art | 42025 | |
0.75 mm biopsy punch | World Precision Instruments | 504529 | |
Glass microscope slides (75 mm x 50 mm) | Corning | 294775X50 | |
Aquapel (hydrophobic glass treatment) | Pittsburgh Glass Works | 47100 | |
PE-2 polyethylene tubing | Scientific Commodities | B31695-PE/2 | |
1 mL syringes | BD | 309628 | |
27 gauge needles | BD | 305109 | |
Syringe pump | New Era Pump Systems | NE-501 | |
Novec HFE-7500 fluorinated oil (HFE) | 3M | 98-0212-2928-5 | |
FC-40 fluorinated oil | Sigma-Aldrich | F9755 | |
PEG-PFPE surfactant | Ran Biotechnologies | 008-FluoroSurfactant | |
Space heater | Lasko | CD09250 | |
Agarose, low gelling temperature | Sigma-Aldrich | a9414 | |
TE (10X) | Rockland | mb-007 | |
PBS 1X, pH 7.4 | E&K Scientific Products | EK-65083 | |
OptiPrep (density gradient medium) | Sigma-Aldrich | d1556 | |
1H,1H,2H-Perfluoro-1-Octanol (PFO) | Sigma-Aldrich | 370533 | |
Span 80 (sorbitane monooleate) | Sigma-Aldrich | s6760 | |
Hexane | Sigma-Aldrich | 139386 | |
Tween 20 (polysorbate 20) | Sigma-Aldrich | p2287 | |
Lysozyme Type IV | MP Biomedicals | 195303 | |
Mutanolysin | Sigma-Aldrich | M9901 | |
Zymolyase (yeast lytic enzyme) | Zymo Research | e1004 | |
Lysostaphin | Sigma-Aldrich | L7386 | |
Sodium chloride | Sigma-Aldrich | S9888 | |
EDTA | Sigma-Aldrich | E6758 | |
Tris-HCl, pH 7.5, 1M | Invitrogen | 15567-027 | |
Dithiothreitol (DTT) | Teknova | d9750 | |
Lithium dodecyl sulfate | Sigma-Aldrich | L9781 | |
Proteinase K | New England Biosciences | P8107S | |
Ethanol, 200 Proof (100%) | Koptec | V1001 | |
SYBR Green I (nucleic acid stain) | Invitrogen | S7563 | |
PEG 6k | Sigma-Aldrich | 81260 | |
Triton X-100 (octylphenol ethoxylate) | Sigma-Aldrich | t8787 | |
Nextera DNA Library Prep Kit | Illumina | FC-121-1030 | |
Phusion Hot Start Flex Master Mix (High-Fidelity Hot Start Master Mix) | New England Biosciences | m05365 | |
Platinum Multiplex PCR Master Mix (Taq Master Mix) | Applied Biosystems | 4464263 | |
Warmstart 2.0 Bst Polymerase (isothermal polymerase) | New England Biosciences | m0538m | |
NT buffer from Nextera XT kit (neutralization buffer) | Illumina | FC-131-1024 | |
Cold cathode fluorescent inverter | (custom) | (custom) | |
DC power supply | Mastech | HY1503D | |
Zerostat 3 anti-static gun | Milty | 5036694022153 | |
3D-printed centrifuge syringe holder | (custom) | (custom) | See Supplemental Files for 3D print file |
Zymo DNA Clean & Concentrator-5 | Zymo Research | D4003 |
请求许可使用此 JoVE 文章的文本或图形
请求许可This article has been published
Video Coming Soon
版权所属 © 2025 MyJoVE 公司版权所有,本公司不涉及任何医疗业务和医疗服务。