CN105950760A - 同时检测单细胞mRNA表达水平和端粒长度的方法 - Google Patents

Abstract

一种同时检测单细胞mRNA表达水平和端粒长度的方法。包括：(1)单细胞的捕获和裂解；(2)单细胞mRNA和gDNA的分离；(3)单细胞mRNA的扩增和端粒特异性序列的扩增；(4)qPCR检测单细胞基因表达；(5)单细胞相对端粒长度的检测。该方法可以有效地对单细胞进行裂解并释放完整的mRNA和gDNA，然后用生物素标记的引物捕获mRNA，磁珠吸附生物素后将mRNA和gDNA进行物理分离，对单细胞的基因表达和长度分别进行检测，最后可获得同一个细胞的基因表达水平和端粒长度并进行关联分析。本发明所述检测方法准确性和灵敏度高，方法简便，可在单细胞水平对哺乳动物各种细胞进行平行的基因表达和端粒长度检测。

Description

同时检测单细胞mRNA表达水平和端粒长度的方法

技术领域

本发明属于生物分子技术领域，具体涉及同时检测单细胞mRNA表达水平和端粒长度的方法。

背景技术

单细胞技术是近几年发展起来的新的方法，单细胞技术是相对于我们之前的以整个细胞群体为研究对象而言的。以往我们都是以一个细胞群体为单位从DNA和RNA水平对生物学现象进行研究，但是由于群体细胞异质性的存在，每个细胞群中的细胞不都是一模一样的。以肿瘤细胞为例，一个肿瘤组织混杂有基质细胞、血细胞、免疫细胞、肿瘤干细胞等，这其中肿瘤干细胞往往在肿瘤的发生发展中起着关键作用，但是肿瘤干细胞的数目非常稀少，这样以传统的群体对象对肿瘤组织进行研究，就会使得我们不能把握关键的信息。另外对于稀少的珍贵样品的研究，以往的群体研究手段都需要大量的细胞，而单细胞技术可以很好的克服这一瓶颈。

某一特定时段，由细胞以基因组DNA为模板转录出来的所有的转录产物称为转录组，其中mRNA可以进一步翻译成蛋白质，从而发挥重要功能。正如以上所述，以往的转录学技术也都是以群体为研究对象，这样将大量细胞的基因表达平均化进行分析，使得我们无法揭示不同细胞之间基因表达的异质性，而对于胚胎细胞这些稀有珍贵的细胞也无能为力。

端粒是存在于染色体末端的重复序列，由TTAGGG重复片段构成，其长度在不同组织和细胞类型中不同，平均长度分布在3-21kb，端粒结构稳定和长度的维持是细胞进行持续增殖的必要条件之一。端粒长度异质性与细胞的功能密切相关，目前群体细胞端粒长度的测量方式主要有TRF(Kimura et al.,2010)、qFISH(Lansdorp et al.,1996)、Flow-FISH(Hultdin et al.,1998)、T/S qPCR(Cawthon,2002)，但是这些方法都是以细胞群体为研究对象，而且需要大量的起始细胞数目。另外对于稀有样品也没有办法进行分析，并且不能对单个细胞的端粒长度进行准确测量。

目前主要存在的技术问题：

目前测量端粒的方法需要大量的gDNA或者大量的细胞数目，对单个细胞和微量gDNA不适用，虽然目前单独的检测单细胞mRNA表达水平的方法已经建立，我们实验室也建立了分析单个细胞端粒长度的方法，但是还没有办法做到同时对同一个细胞的mRNA表达水平和端粒长度进行检测。而仅仅依靠直接高速离心的方法来分离单细胞的细胞核和细胞，由于细胞核中也存在mRNA，这样就会造成较大的误差。

发明内容

本发明目的是解决在单细胞水平同时进行端粒长度和mRNA基因表达谱关联分析的问题，提供一种同时检测单细胞mRNA表达水平和端粒长度的方法。

本发明将在实验室之前建立的单细胞端粒长度检测技术的基础上，通过有效分离单细胞mRNA和gDNA，结合高效的单细胞Smart-seq2转录组扩增技术来用一个单细胞对端粒长度和mRNA基因表达谱关联进行分析，为更加深入研究单细胞异质性提供参考。

本发明技术方案

为了达到同时检测同一个细胞的端粒长度和转录组表达谱，本发明设计了如下方案(见图1)：

第1步、同时获得完整的gDNA和RNA，采用Qiagen RLT Plus裂解单细胞以达到同时释放gDNA和RNA的目的，使用生物素标记的Oligo-dT对含有Ploy(A)的mRNA进行杂交捕获；

第2步、使用能够特异性吸附生物素的磁珠将mRNA和gDNA进行分离，上清中即为游离的gDNA，管壁磁珠中为mRNA；

第3步、将管壁磁珠中分离获得的单细胞mRNA进行反转录成第一链cDNA，并对其进行扩增，扩增产物纯化后，对看家基因GAPDH和Actin表达情况进行qPCR调查，排除看家基因未被扩增的单细胞样品；

第4步、将上清中的gDNA全部转移至新的PCR管中，纯化后用EB溶解gDNA，然后直接向管中加入端粒和内参基因Alu引物进行预扩增反应，对扩增产物进行纯化后分别用于端粒和Alu的qPCR定量分析。

本发明的优点和积极效果：

本发明能够同时对单个细胞的mRNA表达水平和端粒长度进行检测，操作简便，结果可靠，通量高耗费低。

附图说明

图1是单细胞同时检测端粒长度和RNA-seq的基本流程。

图2是单细胞同时检测端粒长度和mRNA表达技术的质量控制。(A)进行18个循环的单细胞Tel和Alu预扩增后，对Tel和Alu进行qPCR定量分析，同时以PBS和Water作为阴性对照。(B)使用Smat-seq2扩增同一个单细胞的mRNA，对扩增获得的dscDNA进行纯化和浓度测定用于qPCR定量分析GAPDH和OCT4基因的表达，同时以PBS和Water作为阴性对照。(C)单细胞端粒qPCR定量分析的端粒扩增曲线。(D)单细胞端粒qPCR定量分析的端粒溶解曲线。(E)单细胞端粒qPCR定量分析的端粒扩增曲线。(F)用于评价同一次实验Tel和Alu扩增效率的标准品Ct值与DNA上样量的标准曲线。S1-S5为标准曲线，总gDNA上样量依次为100ng，20ng，4ng，0.8ng，0.16ng。

图3是scTel&RNA-seq技术分析六种人细胞系的端粒长度异质性和与群体端粒长度相关性分析。(A)-(F)分别为人胚胎干细胞WA26、RUES2、骨肉瘤细胞系U2OS、结肠癌细胞系HCT116、宫颈癌细胞系HeLa、人胚胎成纤维细胞系HEF和Mean±SD。(G)Southernblot杂交TRF定量人不同细胞系的端粒长度。(H)T/S ratio qPCR的方法定量人不同细胞系的相对端粒长度。(I)scTel&RNA-seq单细胞端粒长度的平均值。(J)TRF杂交定量的群体细胞端粒实际平均长度与T/S ratio qPCR定量相对长度的相关性。(K)单细胞端粒长度与TRF杂交定量的群体细胞端粒的相关性。(L)单细胞端粒长度与T/S ratio qPCR定量相对端粒长度的相关性。

图4是单细胞Smart-seq2扩增后qPCR验证单细胞GAPDH基因表达水平，qPCR循环数为40，对扩增产物进行2％凝胶电泳分析，条带均为单一条带。

图5是scTel&RNA-seq检测人胚胎干细胞WA26的OCT4和TRF1的基因表达水平与对应的端粒长度。(A)单细胞GAPDH的Ct值。(B)单细胞OCT4的Ct值。(C)单细胞TRF1的Ct值。(D)单细胞OCT4相对基因表达水平，以GAPDH为内参。(E)单细胞TRF1相对基因表达水平，以GAPDH为内参。总共扩增94个WA26单细胞，以GAPDH判断有9个单细胞mRNA未被扩增，不用于后续分析，扩增成功率为91.5％。(F)qPCR检测的单细胞端粒Ct值。(G)qPCR检测的单细胞Alu基因Ct值。(H)单细胞端粒相对长度的定量分析。S1-S5为标准曲线，总gDNA上样量依次为100ng，20ng，4ng，0.8ng，0.16ng。以Tel Ct值判断有5个单细胞未被扩增，不用于后续分析，扩增成功率为94.7％。总共扩增94个WA26单细胞，NTC1为PBS，NTC2为Water。

图6是单细胞基因表达水平与单细胞端粒长度的相关性和异质性分析。(A)and(B)将85个有效扩增的细胞分成高中低三个群，每个群体间的平均值具有显著性差异。不管是OCT4还是TRF1，相对表达水平高的细胞数目都是很少的，而大部分细胞都集中在中和低水平。(C)将单细胞端粒长度分成4个不同的区间，4个区间均显示出了显著性差异。(D)OCT4与TRF1的共表达情况分析。(E)单细胞端粒长度与OCT4的共表达情况分析。(F)单细胞端粒长度与TRF1的共表达情况分析。(G-J)不同端粒长度区间[0.0-0.4]，[0.4-0.8]，[0.8-2.1]，[2.1-5.0]中单细胞端粒长度与对应OCT4基因表达相关性的分析。(K-N)不同端粒长度区间中单细胞端粒长度与对应TRF1基因表达相关性的分析。

具体实施方式

实施例1、同时检测单细胞mRNA表达水平和端粒长度的方法

本发明方法的具体操作步骤如下：

步骤1、磁珠分离单细胞基因组gDNA和mRNA

注意：本步骤中的用于孵育的PCR仪热盖温度为105℃。使用RnaseZap和DNA-OFF擦拭PCR级别超净工作台。RnaseZap擦拭所使用的移液器。保证所用的试剂盒物品都是Rnase free和无DNA污染。实验时穿戴白大褂，发罩，口罩，鞋罩，手套。

1)首先向PCR样品管中分装好4μL的裂解液，裂解液的配方为：2.5μL的BufferRLT plus裂解液，1μL浓度为10μM的Biotin-oligo-dT，0.03μL的40UμL^-1RNA酶抑制剂RRI，0.975μL的水；

2)将所需要的细胞重悬成单细胞，并将其稀释至10倍物镜下能够找到2-3个细胞的浓度，使用口吸管、流式分选仪或者其他可取的设备分选单个细胞到步骤1)中准备好的PCR管中，加入的体积最大为1μL；

3)迅速小心涡旋，室温700g离心10秒，立即置于冰上；

4)在带热盖的PCR仪中72℃孵育3分钟，结束后立即置于冰上5分钟；

5)室温700g离心10秒，立即置于冰上；此时Biotin-oligo-dT引物已经杂交到含有poly(A)的mRNA上；

6)按照1：1的比例加入5μL的链霉亲和素磁珠MyOne^TM Streptavidin C1室温孵育20分钟；

7)结束后将PCR管置于96孔贝克曼磁体上静置8分钟，沉淀吸附mRNA，将上清转到新的PCR管中，对含mRNA的PCR管使用10μL以下缓冲液清洗2次后，都转到装gDNA的PCR管中，最终即分别获得了用一个单细胞的mRNA和gDNA样品；清洗缓冲液配方是：2.5ml的1M Tris-HCl，25μL的3M KCl，150μL的1M MgCl₂，500μL的1M DTT，250μL的100％Tween-20，总体积为50ml，临用前按照1：40加入RRI；

步骤2、单细胞mRNA的扩增

1)向步骤1-7)得到的mRNA磁珠沉淀中加入10μL的单细胞反转录混合液：0.5μL的200UμL^-1SuperScript II reverse transcriptase，0.25μL的40UμL^-1RRI，2μL的5×Superscript II first-strand buffer，0.5μL的100mM DTT，2.0μL的5M Betaine，0.06μL的1M MgCl₂，0.1μL的100μM TSO引物，4.59μL水，总体积10μL；

2)室温700g离心10秒，置于PCR仪上进行如下循环，42℃90min，50℃2min、42℃2min，10个循环，70℃15min，4℃保存；所得产物即为第一链cDNA；

3)反转结束后，立即进行以下的扩增步骤，直接向以上步骤2-2)的反应液中加入15μL的扩增体系；扩增体系为：12.5μL 2×KAPA HiFi HotStart Ready Mix，0.25μL的10μM ISPCR引物，2.25μL水；

4)使用以下程序进行PCR：98℃3min，98℃20s、67℃15s、72℃6min，18-22个循环，72℃5min，4℃保存；此步骤的循环数根据细胞类型需要调整，结束后即获得双链dscDNA；

5)反应结束后进行双链cDNA的纯化，涡旋震荡混匀AMPure XP beads，室温放置30min，按照1：1的比例向以上步骤2-4)PCR产物中加入25μL磁珠，室温孵育5分钟；磁珠用200μL新鲜配制的80％乙醇漂洗两次，最后加入25μL灭菌超纯水洗脱，于-20℃保存备用；

6)对纯化产物使用QuBit测定浓度，将所有样品调整成浓度为0.5ng/μL，体积为50μL，用于qPCR检测，20μL qPCR体系中加入5μL的dscDNA，最后以GAPDH作为内参，按照2^{- ΔΔCt}计算基因相对表达水平；

步骤3、单细胞端粒长度的检测

1)将以上步骤1-7)分离获得的gDNA按1：0.6比例加入AMPure XP beads进行纯化，最后用10μL水洗脱；

2)向洗脱液中直接加入15μL以下的单细胞端粒长度检测预扩增体系：2.5μL 10XiTaq buffer，0.75μL的50mM MgCl₂，0.5μL的10μM dNTP混合液，0.5μL的10μM Tel-F/R引物混合液，0.5μL的10μM Alu-F/R引物混合液，10.125μL水，0.125μL的5U iTaq DNApolymerase，总反应体系为25μL；

3)按照以下的程序进行PCR预扩增：95℃3分钟，95℃30秒、60℃30秒，18个循环，72℃30秒，10℃保存；

4)预扩增结束后按照1：1.8比例加入步骤3-3)AMPure XP beads纯化预扩增产物，最后加入64μL的水洗脱预扩增产物；

5)按照以下的20μL体系进行qPCR检测Tel和Alu，10μL 2×SYBR green Supermix，上游下游引物各0.5μM，4μL水，步骤3-4)扩增产物5μL；

6)按照以下qPCR程序进行单细胞端粒长度的定量分析，95℃10min，95℃15s、60℃1min，40个循环，溶解曲线为55℃到95℃80个循环；

7)单细胞端粒长度的计算方法:首先根据标准曲线，判断Tel和Alu的扩增效率，当两者效率均接近100％时，则按照Ct值的2^-ΔΔCt方法计算相对的端粒长度。

实施例2：效果杨峥

同时检测单细胞端粒长度和RNA-seq(scT&R-seq)的方法

为了达到同时检测同一个细胞的端粒长度和转录组表达谱，我们设计了图1中所示的方案：

首先同时获得完整的gDNA和RNA，我们采用Qiagen RLT Plus裂解单细胞以达到同时释放gDNA和RNA的目的，使用生物素(Biotin)标记的Oligo-dT对含有Ploy(A)的mRNA进行杂交捕获。

随后，使用Biotin特异性吸附的磁珠将mRNA和gDNA进行分离，上清中即为游离的gDNA，管壁磁珠中为mRNA。

最后，将分离获得的mRNA按照Smart-seq2的方法进行转录组的扩增，对扩增产物进行看家基因表达(GAPDH，Actin)情况进行qPCR调查，后续实验中剔除内参基因未被扩增的样品，对合格的样品进行选择性的qPCR分析或者直接进行二代测序分析。

将上清中的gDNA全部转移至新的PCR管中，并用磁珠进行两次纯化去除影响后续PCR反应的残留试剂，最后用EB溶解gDNA，直接向管中加入同时扩增端粒和Alu序列的PCR预混液，对扩增产物进行纯化后分别用于端粒和Alu的qPCR定量分析。

具体过程如图1所示，使用口吸管或者FACS等任何可以有效分离单细胞的装置将单个细胞分装进单独的管子中，管中提取分装好裂解液，裂解液中含有生物素标记的OligodT，PCR中杂交5分钟后，Oligo dT即结合到含有Poly(A)尾巴的mRNA上，随后向其中加入含有链霉亲和素标记的磁珠，特异性的将生物素标记的mRNA捕获，再将管子放在磁体上进行孵育，mRNA则被吸附在管壁，而gDNA不被吸附游离在上清中，最后将上清转移至一个新的管子中进行单细胞端粒长度方法扩增，剩下的mRNA则立即进行Smart-seq2的单细胞扩增。对最终扩增获得双链cDNA进行质量评价，即可用于qPCR分析或者直接进行高通量的测序分析。对最终获得T/R扩增产物进行纯化后用于单细胞相对端粒长度的检测。

一、scT&R-seq技术的质量控制

我们首先用一种端粒重组机制阳性的肿瘤细胞系U2OS(具有非常长的端粒长度)和一种hESC(WA26)进行实验。每个细胞系挑取6个单细胞，同时以PBS和Water作为对照，按照以上的方法进行操作，Smart-seq2的dscDNA扩增循环数为20，T/R ratio预扩增循环数为18(图2A)。对照(PBS和Water)的Tel和Alu的Ct值均在20以上，对于Ct值在20以上的细胞样品则判断为端粒未被扩增，而能够被扩增的样品，其Ct值要在20以下。对同一个样品同时扩增后的dsDNA进行内参基因GAPDH和多能性基因OCT4的qPCR定量分析，结果表明GAPDH的Ct均接近15，而OCT4Ct在肿瘤细胞U2OS中很高，超过35，而在人胚胎干细胞WA26中Ct值在12左右，这与细胞群体的qPCR结果一致((图2B)。我们使用群体细胞HEF的gDNA作为标准品，来评价每次的扩增效率，结果表明单细胞的端粒扩增Ct值均落在标准曲线内，而对照组(PBS，Water)的扩增曲线则偏离标准品((图2C)，且端粒Tel的溶解曲线与标准品的曲线完全一致((图2D)。内参基因Alu的扩增情况也与端粒leis((图2E)，标准品的Tel和Alu扩增的Ct值与gDNA上样量相关性均为0.99，效率均在100％左右((图2F)。

二、scT&R-seq分析人不同细胞系的端粒长度和异质性

在初步建立同时检测端粒长度和mRNA表达水平方法后，我们首先选取人的不同类型细胞系来验证方法的可靠性(图.3)。结果表明不同的细胞系表现了不同的端粒长度异质性(图.3A-F),端粒重组机制的U2OS细胞间端粒长度差异最大(SD＝0.5)，结肠癌细胞系HCT116和胚胎成纤维细胞系HEF差异最小。而且U2OS的单细胞端粒长度普遍比其他细胞系单细胞端粒长(图.3G)，并且通过端粒长度检测金标准TRF和qPCR的方法获得的群体结构也表明U2OS具有最长的单细胞长度(图.3H and 3I)。相对于其他的细胞，hESC具有较长的端粒长度，其次是HEF和HeLa，结肠癌细胞系HCT116的端粒最短，TRF显示其平均端粒长度只有3.5kb(图.3H)。相关性分析显示，TRF检测结果与qPCR检测端粒长度(T/S ratio)的结果具有非常高的相关性，达到了0.918(P＝0.0026)(图.3J)，而TRF与T/R ratio的也有极显著的相关性为0.933(P＝0.0017)(图.3K)，T/R ratio与T/S ratio相关性为0.787(P＝0.0183)(图.3L)，以上结果表明scT&R-seq检测单细胞端粒长度是可靠的和准确的。

三、scTel&RNA-seq检测人胚胎干细胞WA26的OCT4和TRF1的基因表达和端粒长度异质性。

我们的结果表明端粒长度的异质性普遍存在于人的不同细胞系中，我们进一步使用scTel&RNA-seq来研究端粒长度异质性与hESC WA26多能性基因OCT4和TRF1表达的相关性(图5)。总共挑选了94个WA26单细胞，使用scTel&RNA-seq方法同时扩增单细胞的mRNA和端粒序列。以GAPDH的扩增情况来判断mRNA是否被扩增(Ct值大于26)，结果有85个单细胞获得好的扩增效果(图5A)，并且对qPCR产物进行凝胶电泳判断条带大小和单一性，结果表明与判断的情况一致，对照无目的条带，判定为未被扩增的样品无目的条带或大小不一致(图4)，扩增成功率为91.5％(图5A)。同时进行OCT4和TRF1基因的表达分析(图5B and C)，最后以GAPDH作为内参，计算每个细胞的相对基因表达水平(图5D and E)，结果表明OCT4和TRF1在不同的细胞中显示出了极大的异质性。以Tel或Alu的Ct值大于20判断gDNA未被扩增(图5F and G)，结果显示有9个细胞的gDNA未被扩增，同样的单细胞长度也显示了高的异质性(图5H)。同时mRNA和gDNA均未被扩增的细胞有1个，在后续分析中，这些未被扩增的细胞将被过滤掉。

四.单细胞基因表达水平与单细胞端粒长度异质性的相关性分析

对OCT4和TRF1基因表达水平进行分析后发现，这85个有效扩增的细胞可以明显分成高中低三个群(细胞数分别为6，23，56)，且每个群体间的具有显著性差异。不管是OCT4还是TRF1，相对表达水平高的的细胞数目都是很少的(7％左右)，而大部分细胞都集中在中低水平(图6A and B)。同样的单细胞的端粒长度也可以划分成4个明显差异的细胞群体，大部分的细胞的端粒长度集中在0.4-2.1范围(53个)，短端粒的细胞数有18个，特别长端粒的有14个(图6C)。OCT4和TRF1的基因共表达分析结果表明，WA26细胞系中OCT4与TRF1具有低的相关性，相关性系数R为0.2781(P＝0.01)，以Low表达值进行截断将细胞划分为4个区域，结果表明同时高表达OCT4和TRF1的细胞数为10个(占11.76％)，OCT4和TRF1单独高表达的分别占到20％(17个)和15.29(13个)(图6D)。同时我们发现单细胞端粒长度与OCT4基因表达水平没有显著相关性，相关系数R为0.08771(P＝0.4481)(图6E)。但是单细胞端粒长度与TRF1基因表达水平显示出了极显著的相关性，相关系数R为0.745(P<0.0001)(图6E)。进一步对不同端粒长度的细胞群体中OCT4和TRF1基因表达的相关性进行分析(图6G-N)，结果表明不同端粒长度的细胞OCT4与端粒长度的正相关性不高，特别是端粒相对较短的细胞群体中，呈现除了负相关的趋势(图6I and J)。TRF1基因表达与单细胞端粒长度有一定相关性，特别是长端粒的细胞，与端粒长度范围在2.1-5.0的相关性为0.4624(P＝0.0960)(图6K)，与范围在0.8-2.1的相关性0.5365(P＝0.0015)(图6L)，而与短端粒范围(0.0-0.8)的相关性不高(图6M and N)。

参考文献

Angermueller,C.,Clark,S.J.,Lee,H.J.,MacaμLay,I.C.,Teng,M.J.,Hu,T.X.,Krueger,F.,Smallwood,S.A.,Ponting,C.P.,Voet,T.,et al.(2016).Parallel single-cell sequencing linkstranscriptional and epigenetic heterogeneity.Nature methods 13,229-232.

Cawthon,R.M.(2002).Telomere measurement by quantitative PCR.Nucleic acids research 30,e47.Dey,S.S.,Kester,L.,Spanjaard,B.,Bienko,M.,and van Oudenaarden,A.(2015).Integrated genomeand transcriptome sequencing of the same cell.Nature biotechnology 33,285-289.

Han,A.,Glanville,J.,Hansmann,L.,and Davis,M.M.(2014).Linking T-cell receptor sequenceto functional phenotype at the single-cell level.Nature biotechnology 32,684-692.

Hou,Y.,Guo,H.,Cao,C.,Li,X.,Hu,B.,Zhu,P.,Wu,X.,Wen,L.,Tang,F.,Huang,Y.,et al.(2016).Single-cell triple omics sequencing reveals genetic,epigenetic,and transcriptomicheterogeneity in hepatocellμLar carcinomas.Cell research 26,304-319.

HμLtdin,M.,Gronlund,E.,Norrback,K.,Eriksson-Lindstrom,E.,Just,T.,and Roos,G.(1998).Telomere analysis by fluorescence in situ hybridization and flow cytometry.Nucleic acids research26,3651-3656.

Kimura,M.,Stone,R.C.,Hunt,S.C.,Skurnick,J.,Lu,X.,Cao,X.,Harley,C.B.,and Aviv,A.(2010).Measurement of telomere length by the Southern blot analysis of ter分钟al restrictionfragment lengths.Nature protocols 5,1596-1607.

Lansdorp,P.M.,Verwoerd,N.P.,van de Rijke,F.M.,Dragowska,V.,Little,M.T.,Dirks,R.W.,Raap,A.K.,and Tanke,H.J.(1996).Heterogeneity in telomere length of human chromosomes.HumanmolecμLar genetics 5,685-691.

MacaμLay,I.C.,Haerty,W.,Kumar,P.,Li,Y.I.,Hu,T.X.,Teng,M.J.,Goolam,M.,Saurat,N.,Coupland,P.,Shirley,L.M.,et al.(2015).G&T-seq:parallel sequencing of single-cell genomesand transcriptomes.Nature methods 12,519-522.

Patel,A.P.,Tirosh,I.,Trombetta,J.J.,Shalek,A.K.,Gillespie,S.M.,Wakimoto,H.,Cahill,D.P.,Nahed,B.V.,Curry,W.T.,Martuza,R.L.,et al.(2014).Single-cell RNA-seq highlightsintratumoral heterogeneity in primary glioblastoma.Science 344,1396-1401.

PetropoμLos,S.,Edsgard,D.,Reinius,B.,Deng,Q.,PanμLa,S.P.,Codeluppi,S.,Plaza Reyes,A.,Linnarsson,S.,Sandberg,R.,and Lanner,F.(2016).Single-Cell RNA-Seq Reveals Lineage andX Chromosome Dynamics in Human Preimplantation Embryos.Cell.

Shin,J.,Berg,D.A.,Zhu,Y.,Shin,J.Y.,Song,J.,Bonaguidi,M.A.,Enikolopov,G.,Nauen,D.W.,Christian,K.M.,分钟g,G.L.,etal.(2015).Single-Cell RNA-Seq with Waterfall Reveals MolecμLarCascades underlying AdμLt Neurogenesis.Cell stem cell 17,360-372.

Tang,F.,Barbacioru,C.,Bao,S.,Lee,C.,Nordman,E.,Wang,X.,Lao,K.,and Surani,M.A.(2010a).Tracing the derivation of embryonic stem cells from the inner cell mass by single-cell RNA-Seqanalysis.Cell stem cell 6,468-478.

Tang,F.,Barbacioru,C.,Nordman,E.,Li,B.,Xu,N.,Bashkirov,V.I.,Lao,K.,and Surani,M.A.(2010b).RNA-Seq analysis to capture the transcriptome landscape of a single cell.Nature protocols5,516-535.

Wang,F.,Pan,X.,Kalmbach,K.,Seth-Smith,M.L.,Ye,X.,Antumes,D.M.,Yin,Y.,Liu,L.,Keefe,D.L.,and Weissman,S.M.(2013).Robust measurement of telomere length in single cells.Proceedingsof the National Academy of Sciences of the United States of America 110,E1906-1912.

Zeisel,A.,Munoz-Manchado,A.B.,Codeluppi,S.,Lonnerberg,P.,La Manno,G.,Jureus,A.,Marques,S.,Munguba,H.,He,L.,Betsholtz,C.,et al.(2015).Brain structure.Cell types in the mousecortex and hippocampus revealed by single-cell RNA-seq.Science 347,1138-1142.

Claims (2)

1.一种同时检测单细胞mRNA表达水平和端粒长度的方法，该方法包括：

第1步、同时获得完整的gDNA和RNA，采用Qiagen RLT Plus裂解单细胞以达到同时释放gDNA和RNA的目的，使用生物素标记的Oligo-dT对含有Ploy(A)的mRNA进行杂交捕获；

第2步、使用能够特异性吸附生物素的磁珠将mRNA和gDNA进行分离，上清中即为游离的gDNA，管壁磁珠中为mRNA；

第3步、将管壁磁珠中分离获得的单细胞mRNA进行反转录成第一链cDNA，并对其进行扩增，扩增产物纯化后，对看家基因GAPDH和Actin表达情况进行qPCR调查，排除看家基因未被扩增的单细胞样品；

第4步、将上清中的gDNA全部转移至新的PCR管中，纯化后用EB溶解gDNA，然后直接向管中加入端粒和内参基因Alu引物进行预扩增反应，对扩增产物进行纯化后分别用于端粒和Alu的qPCR定量分析。

2.根据权利要求1所述的同时检测单细胞mRNA表达水平和端粒长度的方法，其特征在于该方法的具体操作步骤如下：

步骤1、磁珠分离单细胞基因组gDNA和mRNA

1)首先向PCR样品管中分装好4μL的裂解液，裂解液的配方为：2.5μL的BufferRLT plus裂解液，1μL浓度为10μM的Biotin-oligo-dT，0.03μL的40UμL^-1RNA酶抑制剂RRI，0.975μL的水；

2)将所需要的细胞重悬成单细胞，并将其稀释至10倍物镜下能够找到2-3个细胞的浓度，使用口吸管、流式分选仪或者其他可取的设备分选单个细胞到步骤1)中准备好的PCR管中，加入的体积最大为1μL；

3)迅速小心涡旋，室温700g离心10秒，立即置于冰上；

4)在带热盖的PCR仪中72℃孵育3分钟，结束后立即置于冰上5分钟；

5)室温700g离心10秒，立即置于冰上；此时Biotin-oligo-dT引物已经杂交到含有poly(A)的mRNA上；

6)按照1：1的比例加入5μL的链霉亲和素磁珠MyOne^TMStreptavidin C1室温孵育20分钟；

7)结束后将PCR管置于96孔贝克曼磁体上静置8分钟，沉淀吸附mRNA，将上清转到新的PCR管中，对含mRNA的PCR管使用10μL以下缓冲液清洗2次后，都转到装gDNA的PCR管中，最终即分别获得了用一个单细胞的mRNA和gDNA样品；清洗缓冲液配方是：2.5ml的1M Tris-HCl，25μL的3M KCl，150μL的1M MgCl2，500μL的1M DTT，250μL的100％Tween-20，总体积为50ml，临用前按照1：40加入RRI；

步骤2、单细胞mRNA的扩增

1)向步骤1-7)得到的mRNA磁珠沉淀中加入10μL的单细胞反转录混合液：0.5μL的200UμL^-1SuperScript II reverse transcriptase，0.25μL的40UμL^-1RRI，2μL的5×Superscript II first-strand buffer，0.5μL的100mM DTT，2.0μL的5M Betaine，0.06μL的1M MgCl₂，0.1μL的100μM TSO引物，4.59μL水，总体积10μL；

2)室温700g离心10秒，置于PCR仪上进行如下循环，42℃90min，50℃2min、42℃2min，10个循环，70℃15min，4℃保存；所得产物即为第一链cDNA；

3)反转结束后，立即进行以下的扩增步骤，直接向以上步骤2-2)的反应液中加入15μL的扩增体系；扩增体系为：12.5μL 2×KAPA HiFi HotStart Ready Mix，0.25μL的10μM ISPCR引物，2.25μL水；

4)使用以下程序进行PCR：98℃3min，98℃20s、67℃15s、72℃6min，18-22个循环，72℃5min，4℃保存；此步骤的循环数根据细胞类型需要调整，结束后即获得双链dscDNA；

5)反应结束后进行双链cDNA的纯化，涡旋震荡混匀AMPure XP beads，室温放置30min，按照1：1的比例向以上步骤2-4)PCR产物中加入25μL磁珠，室温孵育5分钟；磁珠用200μL新鲜配制的80％乙醇漂洗两次，最后加入25μL灭菌超纯水洗脱，于-20℃保存备用；

6)对纯化产物使用QuBit测定浓度，将所有样品调整成浓度为0.5ng/μL，体积为50μL，用于qPCR检测，20μL qPCR体系中加入5μL的dscDNA，最后以GAPDH作为内参，按照2^{- ΔΔCt}计算基因相对表达水平；

步骤3、单细胞端粒长度的检测

1)将以上步骤1-7)分离获得的gDNA按1：0.6比例加入AMPure XP beads进行纯化，最后用10μL水洗脱；

2)向洗脱液中直接加入15μL以下的单细胞端粒长度检测预扩增体系：2.5μL 10XiTaq buffer，0.75μL的50mM MgCl₂，0.5μL的10μM dNTP混合液，0.5μL的10μM Tel-F/R引物混合液，0.5μL的10μM Alu-F/R引物混合液，10.125μL水，0.125μL的5U iTaq DNApolymerase，总反应体系为25μL；

3)按照以下的程序进行PCR预扩增：95℃3分钟，95℃30秒、60℃30秒，18个循环，72℃30秒，10℃保存；

4)预扩增结束后按照1：1.8比例加入步骤3-3)AMPure XP beads纯化预扩增产物，最后加入64μL的水洗脱预扩增产物；

5)按照以下的20μL体系进行qPCR检测Tel和Alu，10μL 2×SYBR green Supermix，上游下游引物各0.5μM，4μL水，步骤3-4)扩增产物5μL；

6)按照以下qPCR程序进行单细胞端粒长度的定量分析，95℃10min，95℃15s、60℃1min，40个循环，溶解曲线为55℃到95℃80个循环；

7)单细胞端粒长度的计算方法:首先根据标准曲线，判断Tel和Alu的扩增效率，当两者效率均接近100％时，则按照Ct值的2^-ΔΔCt方法计算相对的端粒长度。