CN108138365B - 一种高通量的单细胞转录组建库方法 - Google Patents

Abstract

本发明公开了一种高通量的单细胞转录组建库方法，包括：采用纳升级的微量分液平台，将具有预定细胞密度的细胞悬液分配到纳升级的微孔芯片的每个微孔中；用细胞裂解液裂解所述细胞；将反转录试剂分配到所述微孔中进行反转录反应；将cDNA扩增试剂分配到所述微孔中进行cDNA扩增反应，同时采用荧光定量PCR仪监测单细胞扩增情况。

Description

一种高通量的单细胞转录组建库方法

技术领域

本发明涉及细胞的转录组文库构建技术领域，尤其涉及一种高通量的单细胞转录组建库方法。

背景技术

随着基因测序技术水平的提高以及重大国际合作项目的相继开展，基因组研究发展越来越快。2011年，《自然方法》杂志(Nature Methods)将单细胞测序列为年度值得期待的技术之一。2013年，《科学》杂志(Science)将单细胞测序列为年度最值得关注的六大领域榜首。这些都表明，单细胞测序已逐渐成为科研热点，有望成为近年最值得关注的测序技术。

在单细胞研究工作中，扩大试验规模是确保采集足够多的生物多样性信息的关键。使用微流控芯片虽然可以将单个细胞分离，并进行样品扩增，但是通量不高。其它针对单细胞的方法均存在很多缺陷，导致检测灵敏度不高、基因表达信息丢失严重、技术噪音高、操作失误率高、重复性差，而且针对数以千计甚至万计的细胞数量时，实验费用又成为了最大的瓶颈，高昂的单细胞扩增的费用导致单细胞扩增仍然处于不成熟的阶段，在技术上和成本上还远没有达到大规模应用的地步。

发明内容

本发明提供一种高通量的单细胞转录组建库方法，具有高通量、低成本、一体化的特点，能够实现单细胞建库的自动化和规模化。

本发明的高通量的单细胞转录组建库方法，包括：采用纳升级的微量分液平台将具有预定细胞密度的细胞悬液分配到纳升级的微孔芯片的每个微孔中，用细胞裂解液裂解上述微孔内的细胞；采用上述微量分液平台将反转录试剂分配到上述微孔中进行反转录反应；采用上述微量分液平台将cDNA扩增试剂分配到上述微孔中进行cDNA扩增反应，同时采用荧光定量PCR仪辅助监测单细胞扩增情况。

附图说明

图1为本发明实施例中部分样品的扩增溶解曲线，其中上图示出了荧光强度与循环数的关系，下图示出了荧光强度与温度的关系，曲线a显示阳性对照，曲线b显示阴性对照，曲线c、d显示为两个单细胞样品；

图2为本发明实施例中部分样品的Agilent 2100检测结果；

图3为本发明实施例中单细胞样品中检测的基因数目，其中文库1、2、3、4分别代表一个单细胞的文库，rpkm表示每1百万个map上的reads中map到外显子的每1K个碱基上的reads个数(Reads Per Kilobase of exon model per Million mapped reads)。

具体实施方式

下面通过具体实施方式结合附图对本发明作进一步详细说明。

本发明的高通量的单细胞转录组建库方法是基于微量分液平台实现的，因为要实现单细胞的分配，样品的微量化非常关键。目前，微量分液平台的分液能力已经实现了纳升级的突破，相比传统的微升级的微量分液平台，是质的飞跃。本发明正是以纳升级的微量分液平台为基础实现单细胞的分配，进而实现一系列的原位细胞裂解、反转录反应和cDNA扩增反应等，而完成单细胞转录组建厍。

相比现有技术而言，本发明的显著优势体现在高通量、低成本和一体化。所谓高通量是指，相比现有技术尤其是基于微流控芯片的技术，本发明能够一次性获得大量的单细胞，已经证实单个芯片一次即可完成大于1000个单细胞的文库构建。所谓低成本是指，由于本发明方法的简单易行及高通量的特点，而显著降低了每个单细胞的文库构建的平均成本。所谓一体化是指，单细胞分配完成以后的所有步骤都集成于一块芯片上进行。

可用于本发明中的纳升级的微量分液平台，例如Wafergen Biosystems公司的SmartChip^TM MultiSample NanoDispenser(MSND)，可以称作SmartChip^TM多样品纳升级分配器。

本发明中配合使用的微孔芯片是纳升级的微孔芯片，这对于单细胞的分离以及其后的各种反应的进行至关重要。可以采用的纳升级的微孔芯片，例如Wafergen Biosystems公司的SmartChip，可以定制成各种纳升级的容量，可用于本发明的典型的微孔芯片的每个微孔具有350nL的容积。

荧光定量PCR仪用于辅助监测单细胞的扩增情况，根据扩增情况来判断单细胞微孔。可用于本发明的荧光定量PCR仪，例如Wafergen Biosystems公司的SmartChip^TM Real-Time PCR Cycler，可以称作SmartChip^TM实时PCR循环仪。

本发明的方法能够在微孔芯片中一步实现单细胞分离、细胞裂解、反转录反应(RT)以及cDNA扩增反应，根据需要还可以任选地使用TN5转座酶试剂进行片段化处理，以及PCR文库扩增等步骤。使用本发明的方法，单个微孔芯片一次即可完成大于1000个单细胞的文库构建。而且每个反应的试剂用量能够限制在纳升级，例如在本发明的一个实施例中，每个反应的试剂用量在300nL左右，这大大降低了实验成本。此外，转录组测序(RNA-seq)表明该方法构建的文库检测灵敏度高，重复性好，每个单细胞文库至少可以检测到10000个以上的基因。本发明的方法的开发成功，可以有效提高单细胞水平的RNA-seq建库速度，减少人工操作，节省高昂的试剂费用，实现单细胞建库的自动化和规模化。

在本发明中，控制细胞悬液的密度，对于高效率地将单细胞分配到微孔芯片的微孔中是比较关键的。发明人发现，在细胞均匀分布的情况下，当细胞密度小于等于4个/μL时，98％以上的微孔中细胞数是不多于1个。

在本发明的一个优选的实施例中，细胞密度为2-4个细胞/微升，如果细胞密度小于2个细胞/微升，可能有太多微孔中分配不到细胞，降低单细胞的有效分配率；当细胞密度大于4个细胞/微升，可能太多微孔中分配到不止一个细胞，也会相应地降低单细胞的有效分配率。

在本发明中，细胞悬液、细胞裂解液、反转录试剂和cDNA扩增试剂均是通过纳升级的微量分液平台分配到微孔中的，经过大量的实验研究发现，几十纳升级别的用量能够取得较好的效果，比如分配到微孔中的细胞悬液、细胞裂解液、反转录试剂和cDNA扩增试剂的体积可以均为35纳升至50纳升。

在本发明的一个优选的实施例中，采用Wafergen Biosystems公司的MSND作为微量分液平台，其最低分液体积为35纳升，因此在该优选的实施例中，分配到微孔中的细胞悬液、细胞裂解液、反转录试剂和cDNA扩增试剂的体积均为35纳升。同时，考虑到后续一些步骤的试剂用量，我们将微孔芯片的每个微孔的容积定制为350纳升。本领域技术人员可以理解，在一定范围内调整微孔的容积，定制不同规格的微孔芯片对实现本发明也是可行的，例如微孔的容积不小于350纳升均可以。我们分析，上述试剂的用量可以在几十纳升级别的用量范围内波动，然而太多或太少的试剂用量可能导致本发明的有效性降低，比如超过100纳升的用量可能不但单细胞分离效果不好，随后的各步反应试剂用量也会增大，从而提高了成本。因此，本发明中，对于细胞悬液、细胞裂解液、反转录试剂和cDNA扩增试剂，不推荐100纳升以上的试剂用量。

由于本发明是对极微量的试剂进行分液操作，其中可能含有气泡，而气泡在极微量的试剂中可能难以除去，因此最好在将细胞悬液、反转录试剂或cDNA扩增试剂分配到微孔中后，通过离心去除气泡。如果cDNA扩增反应完成以后，还进行其它反应，也最好通过离心去除气泡。

本发明采用荧光定量PCR仪辅助监测单细胞扩增情况。如果以总RNA作为阳性对照，以细胞上清液作为阴性对照，典型的单细胞扩增产物的标志是，溶解曲线为单峰，Tm峰值大于等于88℃，同时Ct值介于阳性对照与阴性对照的数值之间。为了进一步提高单细胞选择的准确性，可以去除Ct值偏高或偏低的情况，例如去除最低和最高10％Ct值的样品。

虽然经过细胞裂解、反转录反应和cDNA扩增反应以后，已经实现了单个细胞的转录组文库的构建。但是，在实际应用中，通常还会进行后续反应。在本发明的一个优选的实施例中，在cDNA扩增反应之后：采用微量分液平台将TN5转座酶试剂分配到微孔中以对cDNA扩增反应的产物进行片段化处理，然后将终止反应试剂分配到微孔中以终止反应；之后，采用微量分液平台将PCR扩增反应试剂分配到微孔中以对片段化处理后的产物进行PCR文库扩增。

结合上述细胞悬液、细胞裂解液、反转录试剂和cDNA扩增试剂的用量情况，TN5转座酶试剂的用量可以为100纳升至150纳升；终止反应试剂的用量可以为35纳升至50纳升；PCR扩增反应试剂的用量可以为50纳升至75纳升。

在本发明的一个优选的实施例中，TN5转座酶试剂的用量为100纳升；终止反应试剂的用量为35纳升；PCR扩增反应试剂的用量为50纳升。

考虑到后续还可能进行混库(pooling)，即将不同单细胞的文库混合在一起。为了使混库之后的不同单细胞的转录组能够区别开来，可以在PCR文库扩增过程中，针对不同的微孔(即不同的单细胞)采用不同的引物对。在实际应用中，不同微孔中的两条引物只要其中一条序列不同即可。基于这样的思想，在本发明的一个优选实施方案中，引物对包括5’端引物和3’端引物，不同5’端引物之间的差别或不同3’端引物之间的差别仅在于各自具有一段专一性的标签序列，而其它部分的序列相同。具体地，标签序列可以是一段具有特定长度N(比如N为6-10的自然数，优选8)的随机序列。标签序列可以位于引物的各个位置，优选位于引物的中段位置。

经过上述PCR文库扩增过程，不同的单细胞文库的两端具有了不同的引物序列。为了实现高通量测序等操作需求，可以将不同单细胞的文库混合在一起。因此，本发明的方法还可以包括，在PCR文库扩增之后：吸取根据荧光定量PCR仪辅助监测到的单细胞孔位中的产物并将其混合在一起。进一步，可以对混合在一起的产物，采用核酸纯化试剂盒纯化并选择预定长度的核酸片段。

以下通过实施例详细说明本发明的技术方案，需要说明的是，实施例用于说明本发明方法的可行性，不应当理解为对本发明保护范围的限制。

本实施例包括单细胞制备、单细胞裂解、反转录(RT)反应、cDNA扩增反应、TN5转座酶建库扩增等一系列步骤。试剂加样操作和PCR操作按照Wafergen Biosystems公司相关仪器说明进行操作。TN5转座酶按照Simone Picelli等(Genome Res.2014.24：2033-2040)发表的方法进行制备。

1.单细胞裂解液制备

配置如下裂解液(表1)，混合均匀，MSND平台进行加样，芯片中每孔加样35nL。2600rcf，12℃离心5min，去除气泡。

表1

裂解液组分	芯片孔用量(nL)
		10％Triton X-100	0.5
40U/μl RNase抑制剂	1.25
		10μM Oligo-dT30VN引物	12.5
10mm dNTP	12.5
		ddH<sub>2</sub>O	8.25
总量	35

Oligo-dT30VN引物序列：

5′-AAGCAGTGGTATCAACGCAGAGTACT30VN-3′(SEQ ID NO：1)。

2.单细胞制备及裂解

2.1.使用Hela S3细胞系进行单细胞文库制备，使用胰酶对贴壁型Hela细胞进行消化，将消化掉的细胞悬浮液取1mL。将细胞在1000rpm条件下离心5min，用1×PBS重悬浮，重复一次。Percoll原液与10×PBS按体积比9∶1混合成细胞的等渗溶液。将重悬的细胞加入配置的20％的等渗溶液中，细胞密度调节至4个/μL，完成细胞样品制备。

2.2.取10pg总RNA作为阳性对照，细胞悬浮液离心获得的上清液作为阴性对照。

2.3.采用MSND平台进行加样，芯片中每孔35nL体积。用膜封好，在2600rcf，12℃条件下离心5min，去气泡。

2.4.转移至可放置芯片的热循环仪中，72℃反应5min进行细胞裂解。裂解后将芯片在2600rcf，12℃条件下离心5min。

3.反转录反应

3.1.配置反转录试剂(表2)，混合均匀，MSND平台进行加样，芯片中每孔加样35nL。2600rcf，12℃离心5min，去除气泡。

表2

反转录试剂组分	芯片孔用量(nL)
		200U/μL SSII	4.2
5×SuperScript II First-Strand Buffer	11.2
		5M Betaine(甜菜碱)	11.2
100mM MgCl<sub>2</sub>	5.04
		100mM DTT	1.4
100μM TSO	0.56
		40U/μL RNAse抑制剂	1.4
总量	35

TSO序列：5′-AAGCAGTGGTATCAACGCAGAGTACATrGrG+G-3′(SEQ ID NO：2，rG表示riboguanosines，+G表示锁核酸修饰)。

3.2.转移芯片至热循环仪中，42℃90min；50℃2min，42℃2min，重复2次；72℃5min。反应后将芯片在2600rcf，12℃条件下离心5min。

4.cDNA扩增反应

4.1.配置cDNA扩增试剂(表3)，混合均匀，MSND平台进行加样，芯片中每孔加样35nL。2600rcf，12℃离心5min，去除气泡。加样完毕后，转移至荧光定量PCR仪中，进行PCR反应及溶解曲线测定。

表3

IS PCR引物：5′-AAGCAGTGGTATCAACGCAGAGT-3′(SEQ ID NO：3)。

4.2.根据Ct值及溶解曲线判定单细胞扩增成功

选取溶解曲线为单峰，Tm峰值＞＝88℃，同时Ct值介于阳性和阴性数值之间的样品。为了保证是单细胞扩增产物，去除最低和最高10％Ct值的样品。图1示出了部分样品的PCR反应及溶解曲线测定结果。

5.TN5转座酶建库

根据Simone Picelli等(Genome Res.2014.24：2033-2040)发表的方法进行制备，TN5酶浓度调至0.75mg/ml，配置5×打断缓冲液(50mM TAPS-NaOH，pH8.5(@RT)，25mMMgCl₂，50％DMF)。

5.1接头(Adapter Mix)制备

5.1.1参考引物名称及序列：

引物A：5′-CTGTCTCTTATACACATCT-3′(SEQ ID NO：4)；

引物B：5′-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-3′(SEQ ID NO：5)；

引物C：5′-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3′(SEQ ID NO：6)。

5.1.2使用退火缓冲液溶解引物A、引物B、引物C至100μM。

5.1.3分别配制如下反应体系，见表4：

表4

5.1.4分别将反应1和反应2涡旋震荡充分混匀，并短暂离心使溶液回到管底。置于PCR仪内，进行如下反应程序：75℃15min；60℃10min；50℃10min；40℃10min；25℃30min。

5.1.5反应结束后，将反应1和反应2等体积混合，混匀。命名为接头(AdapterMix)，-20℃保存。

5.2接头(Adapter Mix)包埋

5.2.1在灭菌PCR管中依次添加各反应组分，见表5：

表5

接头(Adapter Mix)	12.5μL
		TN5	87.5μL

5.2.2使用移液器轻轻吹打，充分混匀。

5.2.3将反应置于25℃反应60min。反应产物命名为片段化酶高级混合物V5S(Tagment Enzyme Advanced Mix V5S)，置于-20℃存。

5.3 DNA片段化

5.3.1配置TN5片段化试剂(表6)，混合均匀，MSND平台进行加样，芯片中每孔加样100nL。2600rcf，12℃离心5min，去除气泡。

表6

片段化试剂组分	芯片孔用量(nL)
		5×缓冲液	45
Tagment Enzyme Advanced Mix V5S	55
		总量	100

5.3.2 55℃热循环仪中温浴7min。

5.4终止反应

配置打断终止反应试剂(表7)，其中384孔板中每孔添加一种N5XX引物，混合均匀，MSND平台选择(72×72模式)加样，芯片中每孔加样35nL，2600rcf，25℃离心5min，去除气泡。室温放置5min。

表7

打断终止反应试剂组分	芯片孔用量(nL)
		2.25％SDS	1.45
N5XX(0.5μM)	6.5
		引物1(10μM)	6.5
引物2(10μM)	6.5
		5×KAPA Fidelity Buffer	14.05
总量	35

引物N5XX序列：

AATGATACGGCGACCACCGAGATCTACACXXXXXXXXTCGTCGGCAGCGTC(SEQ ID NO：7，X表示随机碱基)；

引物1序列：5′-AATGATACGGCGACCACCGA-3′(SEQ ID NO：8)；

引物2序列：5′-CAAGCAGAAGACGGCATACGA-3′(SEQ ID NO：9)。

5.5PCR文库扩增

配置扩增反应试剂(表8)，384孔板中每孔添加一种N7XX引物混合均匀，MSND平台选择(72×72模式)加样，芯片中每孔加样50nL，2600rcf，12℃离心5min，去除气泡。

表8

引物N7XX序列：

CAAGCAGAAGACGGCATACGAGATXXXXXXXXGTCTCGTGGGCTCGG(SEQ ID NO：10，X表示随机碱基)。

将芯片放入热循环仪中进行扩增。

设置如下反应程序：

热盖105℃；72℃3min；95℃3min；(98℃20sec；60℃15sec；72℃25sec)×15-20循环；72℃5min；4℃保温。

6.根据4.2，将判定单细胞成功扩增出来的孔位，使用直径小于200μm的玻璃针或者特制的吸量仪吸取混库(pooling)在一起。

7.将混库后的产物使用Agencourt AMPure XP beads，进行选择性纯化。

7.1涡旋震荡混匀AMPure XP beads并吸取40μL体积至50μL PCR混库后的产物中，使用移液器轻轻吹打10次充分混匀。室温孵育5分钟。

7.2将反应管短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心转移上清至干净EP管中，丢弃磁珠。

7.3涡旋震荡混匀AMPure XP beads并吸取7.5μL体积至上清中，使用移液器轻轻吹打10次充分混匀。室温孵育5分钟。

7.4将反应管短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心移除上清。

7.5保持EP管始终处于磁力架中，加入200μL新鲜配制的80％乙醇漂洗磁珠。室温孵育30秒后小心移除上清。

7.6重复上步，总计漂洗两次。

7.7保持EP管始终处于磁力架中，开盖空气干燥磁珠10分钟。

7.8将EP管从磁力架中取出，加入15μL灭菌超纯水洗脱。涡旋振荡或使用移液器轻轻吹打充分混匀。将反应管短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心吸取上清至灭菌EP管中，于-20℃保存。如需获得长度分布更集中的文库，扩增产物可使用胶回收试剂盒进行片段长度分选和纯化。

8 Agilent 2100检测文库构建结果。结果如图2。

9根据纯化结果选择不同测序类型进行illumina上机测序。测序结果如图3。

以上内容是结合具体的实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (12)

1.一种高通量的单细胞转录组建库方法，其特征在于，所述方法包括：

采用纳升级的微量分液平台将具有预定细胞密度的细胞悬浮液分配到纳升级的微孔芯片的每个微孔中，用细胞裂解液裂解所述微孔内的细胞；

采用所述微量分液平台将反转录试剂分配到所述微孔中进行反转录反应；

采用所述微量分液平台将cDNA扩增试剂分配到所述微孔中进行cDNA扩增反应，同时采用荧光定量PCR仪辅助监测单细胞扩增情况；

所述预定细胞密度是指2-4个细胞/微升；

分配到所述微孔中的细胞悬浮液的体积为35纳升至50纳升；

所述微孔芯片的每个微孔的容积为不小于350纳升。

2.根据权利要求1所述的方法，其特征在于，分配到所述微孔中的细胞裂解液、反转录试剂和cDNA扩增试剂的体积均为35纳升至50纳升。

3.根据权利要求1所述的方法，其特征在于，分配到所述微孔中的细胞悬浮液、细胞裂解液、反转录试剂和cDNA扩增试剂的体积均为35纳升。

4.根据权利要求1所述的方法，其特征在于，将所述细胞悬浮液、反转录试剂或cDNA扩增试剂分配到所述微孔中后，通过离心去除气泡。

5.根据权利要求1所述的方法，其特征在于，所述荧光定量PCR仪辅助监测单细胞扩增情况中，以溶剂曲线为单峰，Tm峰值大于等于88℃，同时Ct值介于阳性对照与阴性对照的数值之间，作为单细胞扩增产物的标志；其中所述阳性对照是总RNA，所述阴性对照是细胞上清液。

6.根据权利要求5所述的方法，其特征在于，所述方法还包括：去除最低和最高10％Ct值的样品。

7.根据权利要求1所述的方法，其特征在于，所述方法还包括，在所述cDNA扩增反应之后：

采用所述微量分液平台将TN5转座酶试剂分配到所述微孔中以对所述cDNA扩增反应的产物进行片段化处理，然后将终止反应试剂分配到所述微孔中以终止反应；

采用所述微量分液平台将PCR扩增反应试剂分配到所述微孔中以对所述片段化处理后的产物进行PCR文库扩增。

8.根据权利要求7所述的方法，其特征在于，所述TN5转座酶试剂的体积为100纳升至150纳升；所述终止反应试剂的体积为35纳升至50纳升；所述PCR扩增反应试剂的体积为50纳升至75纳升。

9.根据权利要求7所述的方法，其特征在于，所述TN5转座酶试剂的体积为100纳升；所述终止反应试剂的体积为35纳升；所述PCR扩增反应试剂的体积为50纳升。

10.根据权利要求7所述的方法，其特征在于，所述PCR文库扩增过程中，不同的微孔中采用不同的引物对。

11.根据权利要求10所述的方法，其特征在于，所述引物对包括5’端引物和3’端引物，不同5’端引物之间的差别或不同3’端引物之间的差别仅在于各自具有一段专一性的标签序列，而其它部分的序列相同。

12.根据权利要求7所述的方法，其特征在于，所述方法还包括，在所述PCR文库扩增之后：

吸取根据所述荧光定量PCR仪辅助监测到的单细胞空位中的产物并将其混合在一起；

对所述混合在一起的产物，采用核酸纯化试剂盒纯化并选择预定长度的核酸片段。

Images