CN110886021A - 一种单细胞dna文库的构建方法 - Google Patents

一种单细胞dna文库的构建方法 Download PDF

Info

Publication number
CN110886021A
CN110886021A CN201811041751.2A CN201811041751A CN110886021A CN 110886021 A CN110886021 A CN 110886021A CN 201811041751 A CN201811041751 A CN 201811041751A CN 110886021 A CN110886021 A CN 110886021A
Authority
CN
China
Prior art keywords
transposase
genome
library
primer
sequencing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201811041751.2A
Other languages
English (en)
Other versions
CN110886021B (zh
Inventor
王雨舟
吴靓
赵至坤
厉磊
徐礼钦
刘石平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BGI Shenzhen Co Ltd
Original Assignee
BGI Shenzhen Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by BGI Shenzhen Co Ltd filed Critical BGI Shenzhen Co Ltd
Priority to CN201811041751.2A priority Critical patent/CN110886021B/zh
Publication of CN110886021A publication Critical patent/CN110886021A/zh
Application granted granted Critical
Publication of CN110886021B publication Critical patent/CN110886021B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C40COMBINATORIAL TECHNOLOGY
    • C40BCOMBINATORIAL CHEMISTRY; LIBRARIES, e.g. CHEMICAL LIBRARIES
    • C40B50/00Methods of creating libraries, e.g. combinatorial synthesis
    • C40B50/06Biochemical methods, e.g. using enzymes or whole viable microorganisms
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Organic Chemistry (AREA)
  • Biochemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Engineering & Computer Science (AREA)
  • Microbiology (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Molecular Biology (AREA)
  • General Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Biotechnology (AREA)
  • Immunology (AREA)
  • Medicinal Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Genetics & Genomics (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明公开了一种单细胞DNA文库的构建方法。本发明公开的单细胞DNA文库的构建方法包括:利用转座酶复合体打断单细胞基因组DNA,得到片段化基因组,利用片段化基因组构建单细胞全基因组文库;转座酶复合体包括转座酶和转座酶接头,转座酶接头由名称分别为A和C的单链DNA组成;A为所述转座酶的识别序列的反向互补序列;C为含有所述识别序列的单链DNA。本发明方法可以利用皮克级别的单细胞基因组DNA构建单细胞全基因组文库,且基因组覆盖均一性、覆盖度和测序深度高,适合于扩大规模实现高通量的单细胞扩增。

Description

一种单细胞DNA文库的构建方法
技术领域
本发明涉及生物技术领域中,一种单细胞DNA文库的构建方法。
背景技术
多细胞基因组的测序中,对来源于混合细胞的DNA进行高深度的测序。这类方法能够鉴定主要的细胞亚克隆,但受测序错误的限制,对少数细胞群体的鉴定能力弱。多细胞的测序方法会受到细胞数目的限制,且对亚克隆细胞中低频率的DNA拷贝数变异(CNV)和单核苷酸变异(SNV)的检测不灵敏。
单细胞基因组学将基因组的研究延伸到了细胞水平,为对遗传学的理解提供新的观点。单细胞基因组学可用于鉴定和组装不可培养的微生物的基因组,评估遗传镶嵌在正常生理状态和疾病状态下的作用,以及确认肿瘤内的遗传异质性在癌症发展和治疗反应中所起到的作用,单细胞基因组学对珍贵的细胞样品,如胚细胞、循环肿瘤细胞的研究具有重要意义。现有的单细胞基因组技术,主要依赖于测序文库构建前的单细胞全基因组扩增来获得完整的基因组信息,主要有简并寡核苷酸引物PCR、多重置换扩增、DLP直接文库构建、LIANTI技术和TnBC技术。
简并寡核苷酸引物PCR(Degenerate oligonucleotide-primed PCR(DOP-PCR)),利用随机引物对单细胞全基因组进行PCR预扩增,增大起始的DNA量,再进行后续的文库构建。DOP-PCR需要预扩增,会偏好性地扩增基因组上胞嘧啶和鸟苷含量丰富的区域,导致较低的基因组物理覆盖度,不利于单核苷酸变异SNV的检测。
多重置换扩增(Multiple displacement amplification(MDA)),利用随机引物在具有合成能力和链置换活性的聚合酶作用下,对单细胞的基因组DNA进行预扩增,增大起始的DNA量,然后进行文库构建。该方法覆盖度较好,具有较低的错误率,假阳性较低,但均一性较差。该方法需要预扩增,扩增均一性较差,且φ29聚合酶的活性会产生少量嵌合序列副产物,不利于拷贝数变异CNV的检测。
DLP直接文库构建(Direct library preparation(DLP)),利用Tn5酶转座活性直接处理原始单细胞DNA,在定制的微流控仪器上完成文库构建。其中,由于Nextera的Tn5酶的引物设计使用两种不同的接头,只有50%的基因组能够被扩增、测序。此外,该技术需要定制微流控的设备,预先在相应的孔中固定特定的引物,设计较复杂,不利于该技术的应用。
LIANTI技术(Linear Amplification via Transposon Insertion(LIANTI)技术),利用对称性接头打断单细胞基因组DNA,提高转座酶模板利用率,并使用T7 promoter进行转录再逆转录成cDNA,实现线性扩增降低偏好性和PCR错误,可大幅提升单细胞DNA测序的CNV分辨率和碱基准确性。该方法未检测插入片段的大小是否适合于后续的二代测序,片段长度过长会影响测序的覆盖度。此外,转录和反转录的过程时间较长,成本较高。
TnBC技术(Transposon Barcoded(TnBC)技术),利用饱和浓度的Mu转座酶打断原始单细胞DNA,再用与Mu匹配的MuEnd引物进行扩增纯化,纯化后的产物以带有barcode的引物进行文库构建。进行二代测序后,比对到参考基因组,将每个片段起始和末位的几个碱基作为独特片段标记(UFI)。根据UFI去除冗余的序列、连接各DNA片段,以及构成单倍型。该方法经过Mu转座酶打断扩增后,需纯化后再进行带有barcode的引物扩增,易在操作过程中损失部分基因组信息。
发明内容
本发明所要解决的技术问题是如何构建单细胞全基因组文库。
为解决上述技术问题,本发明首先提供了一种单细胞全基因组文库的构建方法,所述方法包括:
利用转座酶复合体打断单细胞基因组DNA,使打断的片段两端的接头相同,得到片段化基因组,利用所述片段化基因组构建单细胞全基因组文库;
所述转座酶复合体包括转座酶和转座酶接头,所述转座酶接头由名称分别为A和C的单链DNA组成;
所述A为所述转座酶的识别序列的反向互补序列;
所述C为含有所述识别序列的单链DNA。
上述方法中,所述C可由所述识别序列与名称为c1的单链DNA组成;所述利用所述片段化基因组构建单细胞全基因组文库包括:利用名称为引物1的引物对所述片段化基因组进行扩增,得到PCR产物,所述PCR产物即为单细胞全基因组文库;所述引物含有所述c1,能识别所述C。
所述c1和所述识别序列可分别位于所述C的5′端和3′端。
所述转座酶复合体具体可由所述转座酶和所述转座酶接头形成的复合体。所述转座酶复合体具体可由所述转座酶和所述转座酶接头在25℃孵育60min得到。
所述识别序列具体为所述转座酶实现转座所依赖的特定的19bp ME(mosaic end)序列。
在利用所述转座酶复合体打断所述样本基因组DNA后,所述片段化基因组的两端分别连接有所述转座酶接头。
上述方法中,所述引物1可由所述C中所述c1及所述识别序列的部分或全部序列组成。
上述方法中,所述c1可为序列表中序列1的第1-15位。
上述方法中,所述转座酶可为Tn5。所述识别序列可为序列表中序列1的第16-34位。
具体的,所述C的序列可为序列表中序列1。
所述引物1的序列具体可为序列表中序列1的第1-22位。
上述方法还可包括对所述PCR产物连接用于测序的接头步骤。
上述方法还可包括:对所述PCR产物连接所述用于测序的接头,得到连接接头的文库;利用识别所述用于测序的接头的引物对所述连接接头的文库进行PCR扩增,得到单细胞全基因组文库。
上述方法中,所述用于测序的接头可为用于高通量测序的接头。
所述高通量测序具体可为利用Illumina Hiseq的高通量测序。
所述用于测序的接头具体可为由序列表中序列2和序列3所示的单链DNA组成的Y型接头。
所述识别所述用于测序的接头的引物可含有标签序列(barcode序列)。所述标签序列的长度可为b1)或b2):
b1)6-16nt;
b2)10nt。
所述标签序列在特定的样本中可为确定的序列,以区分特定样本与其他样本的基因组;所述标签序列在特定的样本中也可为不确定的序列,即其序列为4n种,n为标签序列长度,以在测序过程中对该样本的基因组DNA进行区分。
所述识别所述用于测序的接头的引物具体可为由P5和P7组成的引物对。
本发明还提供了下述X1)-X4)中任一产品:
X1)所述转座酶复合体;
X2)所述转座酶接头;
X3)成套试剂,由X31)和X32)组成;
X31)所述转座酶复合体或所述转座酶接头;
X32)所述引物1;
X4)成套试剂,由X41)、X42)和X43)组成;
X41)所述转座酶复合体或所述转座酶接头;
X42)所述引物1;
X43)所述识别所述用于测序的接头的引物。
X3)和X4)所述成套试剂均可用于构建单细胞全基因组文库或用于单细胞全基因组测序。
本发明还提供了下述任一应用:
Ⅰ、所述单细胞全基因组文库的构建方法在基因组测序中的应用;
Ⅱ、所述产品在基因组测序或构建基因组文库中的应用。
上述应用中,所述基因组测序可为对单细胞的基因组进行的测序。所述基因组文库可为用于单细胞的基因组文库。
本发明以皮克级别的单细胞基因组DNA(<10pg)为模板,在无预扩增的条件下利用只包埋一种接头的转座酶打断基因组,使打断的片段两端的接头相同,提高基因组的覆盖度,实现单细胞全基因组文库构建,减少了PCR扩增的偏向性,提高基因组覆盖均一性和覆盖度;本发明单细胞裂解、扩增等过程反应体系简单可以相互兼容,提高了实验的成功率,降低了测序的成本;实验过程简便、耗时短,减少了样品的损失,使测序信息更完整,适合于扩大规模实现高通量的单细胞扩增,优于LIANTI和TnBC的方法。此外,在测序前扩增又在序列中引入了特异性的Barcode,可将多个样品Pooling到一起再进行低深度的测序,对DNA的拷贝数变异检测灵敏度高、耗时短(整个文库构建过程可在5-6h内完成),精简了单细胞DNA预扩增的过程,在单管单细胞层面实现快速的单细胞全基因组文库构建,适用于发现细胞的亚克隆。本发明还可有效结合流式细胞分选技术和高通量自动化操作平台,获得大量单细胞基因组拷贝数变异信息,应用范围广泛。
附图说明
图1为PCR产物的电泳检测结果。N为阴性对照,100cells表示阳性对照,三个1cell分别为实验组的三个重复。
图2为通过Agilent 2100检测片段分布情况。纵坐标为荧光强度。
图3为两种构建文库方法的覆盖度的比较。单接头表示本发明的单细胞DNA文库的构建方法构建的文库,双接头表示利用双接头方法构建的文库。
图4为两种构建文库方法的测序深度的比较。单接头表示本发明的单细胞DNA文库的构建方法构建的文库,双接头表示利用双接头方法构建的文库。
具体实施方式
下面结合具体实施方式对本发明进行进一步的详细描述,给出的实施例仅为了阐明本发明,而不是为了限制本发明的范围。下述实施例中的实验方法,如无特殊说明,均为常规方法。下述实施例中所用的材料、试剂、仪器等,如无特殊说明,均可从商业途径得到。以下实施例中的定量试验,均设置三次重复实验,结果取平均值。下述实施例中,如无特殊说明,序列表中各核苷酸序列的第1位均为相应DNA的5′末端核苷酸,末位均为相应DNA的3′末端核苷酸。
实施例1、单细胞DNA文库的构建
本实施例提供了一种单细胞DNA文库的构建方法,该方法利用Tn5转座酶包埋含有扩增序列和Tn5识别序列的接头,将单细胞的基因组DNA(<10pg)在无预扩增的情况下,随机打断并引入同一种接头。用引入的接头引物进行无偏向性地扩增,富集DNA。纯化后通过加A和连接反应,在DNA片段两端引入Illumina文库构建的包含P5、P7序列的Y型接头,磁珠纯化后用相应的P5和P7引物扩增,P7引物中带有特异性的Barcode,纯化后即为IlluminaHiseq平台的测序文库。具体操作方法如下:
一、转座酶复合体的制备
1、合成名称分别为A和C的单链DNA,C由名称为c1的扩增接头和转座酶Tn5的识别序列(识别位点)组成,二者直接相连,A的序列为转座酶Tn5的识别序列的反向互补序列。人工合成寡核苷酸时,其5′端核苷酸不含有磷酸,此处需要在其5′端进行磷酸修饰,以使5′末端核苷酸可以与3′末端核苷酸连接。各引物如下:
A:5′-pho-CTGTCTCTTATACACATCT-3′;
C:5′-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3′(序列表中序列1),下划线所示核苷酸序列为c1,加粗核苷酸序列为转座酶Tn5的识别序列。
2、使用Annealing Buffer溶解A和C,浓度均为100μM,其中,Annealing Buffer由1M Tris-HCl溶液pH7.8(T2694-1L,SIGMA)、0.5M EDTA溶液pH8.0(AM9260G,AMBION和5MNaCl氯化钠水溶液(S5150-1L,SIGMA)配制,1L Annealing Buffer中含有10mL 1M Tris-HCl溶液pH7.8,2mL 0.5M EDTA溶液pH8.0和10mL 5M NaCl氯化钠水溶液,余量为水。
配制如下反应体系:
Figure BDA0001792279460000051
将上述反应体系涡旋震荡充分混匀,并短暂离心使溶液回到管底。置于PCR仪内,按照如下反应程序进行反应:
循环数 温度 时间
1 75℃ 15min
1 60℃ 10min
1 50℃ 10min
1 40℃ 10min
1 25℃ 30min
反应结束后混匀,即得到Tn5接头,-20℃保存。
3、Tn5接头包埋
a)在灭菌PCR管中依次添加各反应组分:
成分 体积
步骤2得到的Tn5接头 12.5μL
Tn5转座酶(BGI产品BGE005) 87.5μL
总体积 100μL
b)使用移液器轻轻吹打,充分混匀。
c)将反应置于25℃反应60min,即得到转座酶复合体(命名为Tagment EnzymeAdvanced Mix V5S),置于-20℃存。
二、单细胞DNA文库的构建
1、单细胞裂解液制备
取87μL 30mM Tris-HCl(pH 7.5)和3μL蛋白酶(Qiagen,Cat No./ID:19155QIAGEN Protease(7.5AU))混合即得到的单细胞裂解液,蛋白酶在该溶液中的浓度为0.67mg/mL。
2、单细胞样品制备及裂解
将Hela细胞用PBS洗涤去除外部DNA和细胞碎片,使用流式细胞仪分选单个细胞置入灭菌的装有2μL步骤1的单细胞裂解液的PCR管中,即实验组。以无单细胞的单细胞裂解液作为阴性对照,向单细胞裂解液中添加100个Hela细胞作为阳性对照。将实验组、阴性对照和阳性对照于50℃孵育1h进行反应,然后于75℃孵育20min,80℃孵育5min失活蛋白酶,冷却到10℃,即得到单细胞裂解产物(实验组基因组DNA含量<10pg)。
3、利用Tn5酶对单细胞基因组进行片段化
a)于室温解冻5 X TAG Buffer(华大基因产品,货号为BGE005B01),上下颠倒混匀后备用;
b)按照如下反应体系对单细胞基因组进行片段化:
成分 体积
H<sub>2</sub>0 5.8μL
单细胞裂解产物 2μL
5 X TAG buffer 2μL
步骤1的转座酶复合体 0.2μL
总体积 10μL
c)将上述反应体系于55℃温浴30min,然后向反应体系中添加2.5μL 0.1%SDS水溶液后室温孵育5min,终止反应,置于冰上,即得到单细胞片段化基因组。
按照上述a)-c)的方法,将单细胞裂解产物分别替换为阳性对照或阴性对照进行实验。
4、单细胞片段化基因组的PCR扩增
按照下述方法,对步骤3得到的单细胞片段化基因组进行PCR扩增:
a)引物名称及序列:
所用引物为Partial C primer接头引物序列,人工合成寡核苷酸时,其5′端核苷酸不含有磷酸,此处需要在其5′端进行磷酸修饰,以使5′末端核苷酸可以与3′末端核苷酸连接,引物具体如下:
Partial C primer:5’-pho-GTCTCGTGGGCTCGGAGATGTG-3’
将灭菌PCR管置于冰浴中,按照如下体系配制PCR扩增体系,其中,5 x KAPAFidelity Buffer和KAPA HiFi DNA polymerase均为KAPA Biotech产品:
成分 体积
灭菌蒸馏水 24μL
步骤3的单细胞片段化基因组 12.5μL
5 x KAPA Fidelity Buffer 10μL
10mM each dNTP 1.5μL
Partial C primer(10μM) 1μL
KAPA HiFi DNA polymerase(1U/μL) 1μL
总体积 50μL
将上述PCR扩增体系按照如下扩增程序进行扩增,得到PCR产物:
Figure BDA0001792279460000071
5、PCR产物的磁珠纯化
利用1×Ampure XP bead(Beckman Coulter)纯化步骤4得到的PCR产物。
a)起始PCR产物体积应为50μL。因PCR过程中样品挥发会导致产物体积不足50μL,进行下面操作之前必须使用灭菌蒸馏水将体积补齐至50μL。
b)涡旋震荡混匀AMPure XP beads并吸取50μL体积至含有50μL PCR产物的EP管中,使用移液器轻轻吹打10次充分混匀。室温孵育10分钟。
c)步骤b)结束后进行短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心移除上清。
d)保持EP管始终处于磁力架中,步骤c)结束后,向EP管中加入200μL新鲜配制的80%(体积百分比)乙醇水溶液漂洗磁珠。室温孵育30秒后小心移除上清。
e)重复步骤d),总计漂洗两次。
f)保持EP管始终处于磁力架中,开盖,空气中干燥磁珠3分钟。
g)将EP管从磁力架中取出,向EP管中加入21μL灭菌超纯水,涡旋振荡或使用移液器轻轻吹打充分混匀,将EP管短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心吸取上清至灭菌EP管中,即得到纯化后的PCR产物,于-20℃保存。
6、PCR产物末端加A
将步骤5获得的纯化后的PCR产物在dATP、Klenow(exo-3’→5’)酶(ENZYMATICS,货号P7010-HC-L,50,000U/ml)的作用下,添加一个A碱基,用于后续与Y型接头的连接。配制如下反应体系:
Figure BDA0001792279460000081
37℃反应30min,4℃保持。其中,10X NEBuffer 2.1为NEB产品,货号为B7202S。
7、加A后纯化
将步骤6的产物进行柱纯化,纯化过程参考试剂盒Purelink PCR Micro kit(INVITROGEN,货号A11358),最后用11μl EB洗脱产物。
8、Y型接头制备
a)合成名称分别为Y P5AD和Y P7AD的单链DNA,两序列中有部分序列可互补配对,退火后形成Y型接头。人工合成寡核苷酸时,其5′端核苷酸不含有磷酸,此处需要在其5′端进行磷酸修饰,以使5′末端核苷酸可以与3′末端核苷酸连接。各引物如下:
Y P5AD:5’-pho-TACACTCTTTCCCTACACGACGCTCTTCCGATCT-3′(序列2)
Y P7AD:5’-pho-GATCGGAAGAGCACACGTCTGAACTCCAGTCAC-3′(序列3)
b)使用超纯水溶解Y P5AD和Y P7AD,浓度均为100μM,在10X NEBuffer 2.1缓冲液(NEB,货号B7202S)的作用下制备Y型接头。
配制如下反应体系:
Figure BDA0001792279460000091
将上述反应体系涡旋震荡充分混匀,并短暂离心使溶液回到管底。置于PCR仪内,按照如下反应程序进行反应:
循环数 温度 时间
1 75℃ 15min
1 60℃ 10min
1 50℃ 10min
1 40℃ 10min
1 25℃ 30min
反应结束后混匀,即得到Y型接头,-20℃保存。
9、Y型接头连接
利用T4 DNA连接酶,将步骤7的纯化产物与步骤8获得的Y型接头进行连接反应,得到连接产物。反应体系如下:
Figure BDA0001792279460000092
按上述体系配制后,混匀,23℃,60min连接,4℃保持。
PNK buffer:10x Polynucleotide Kinase Buffer,Enzymatics,货号B904;T4DNA ligase(600U/μl),FERMENTAS,货号EL0011。
10、加接头后磁珠纯化
利用1×Ampure XP bead(Beckman Coulter)纯化步骤9得到的连接产物。
a)起始PCR产物体积应为50μL。因PCR过程中样品挥发会导致产物体积不足50μL,进行下面操作之前必须使用灭菌蒸馏水将体积补齐至50μL。
b)涡旋震荡混匀AMPure XP beads并吸取50μL体积至含有50μL PCR产物的EP管中,使用移液器轻轻吹打10次充分混匀。室温孵育10分钟。
c)步骤b)结束后进行短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心移除上清。
d)保持EP管始终处于磁力架中,步骤c)结束后,向EP管中加入200μL新鲜配制的80%(体积百分比)乙醇水溶液漂洗磁珠。室温孵育30秒后小心移除上清。
e)重复步骤d),总计漂洗两次。
f)保持EP管始终处于磁力架中,开盖,空气中干燥磁珠3分钟。
g)将EP管从磁力架中取出,向EP管中加入21μL灭菌超纯水,涡旋振荡或使用移液器轻轻吹打充分混匀,将EP管短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心吸取上清至灭菌EP管中,即得到纯化后的PCR产物,于-20℃保存。纯化后产物取1μl利用Qubit dsDNA HS Assay Kit(Thermo Fisher Scientific)测定dsDNA浓度。
11、加接头后扩增
用P5(5’-AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT)和P7(5’-CAAGCAGAAGACGGCATACGAGATAACCGCGG GTGACTGGAGTTCAGACGTGT,加粗部分为barcode可根据不同样品更换带有不同barcode的P7引物)组成的引物对对步骤10获得的纯化后产物进行扩增,以获得充足的Illumina测序文库。按如下体系配制:
其中,5 x KAPA Fidelity Buffer和KAPA HiFi DNA polymerase均为KAPABiotech产品:
成分 体积
灭菌蒸馏水 13.5μL
步骤10的纯化产物 20μL
5 x KAPA Fidelity Buffer 10μL
10mM each dNTP 1.5μL
P5(amp)(10μM) 2μL
P7(Tag)(10μM) 2μL
KAPA HiFi DNA polymerase(1U/μL) 1μL
总体积 50μL
将上述PCR扩增体系按照如下扩增程序进行扩增,得到PCR产物:
Figure BDA0001792279460000101
Figure BDA0001792279460000111
12、磁珠纯化
利用1×Ampure XP bead(Beckman Coulter)纯化步骤11得到的PCR产物。
a)起始PCR产物体积应为50μL。因PCR过程中样品挥发会导致产物体积不足50μL,进行下面操作之前必须使用灭菌蒸馏水将体积补齐至50μL。
b)涡旋震荡混匀AMPure XP beads并吸取50μL体积至含有50μL PCR产物的EP管中,使用移液器轻轻吹打10次充分混匀。室温孵育10分钟。
c)步骤b)结束后进行短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心移除上清。
d)保持EP管始终处于磁力架中,步骤c)结束后,向EP管中加入200μL新鲜配制的80%(体积百分比)乙醇水溶液漂洗磁珠。室温孵育30秒后小心移除上清。
e)重复步骤d),总计漂洗两次。
f)保持EP管始终处于磁力架中,开盖,空气中干燥磁珠3分钟。
g)将EP管从磁力架中取出,向EP管中加入21μL灭菌超纯水,涡旋振荡或使用移液器轻轻吹打充分混匀,将EP管短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心吸取上清至灭菌EP管中,即得到纯化后的PCR产物,于-20℃保存。纯化后产物取1μl利用Qubit dsDNA HS Assay Kit(Thermo Fisher Scientific)测定dsDNA浓度。
13.质量与评估
a)质量浓度检测
利用Qubit dsDNA HS Assay Kit(Thermo Fisher Scientific)测定,实验组的浓度大于3.5ng/μl。
b)电泳检测
采用1%琼脂糖凝胶检测纯化后的PCR产物。结果(图1)显示,阴性对照(N)的泳道只有100bp左右的引物二聚体,实验组及阳性对照都有明显的200bp以上的产物条带。
c)Agilent 2100检测
将纯化后的PCR产物稀释到合适的浓度,通过Bioanalyzer High SensitivityDNA Kit(Agilent)测定PCR产物片段大小范围。结果(图2)显示,实验组的PCR产物主要集中在300bp,主要在300bp~500bp之间,符合预期的片段大小。
14、使用Illumina Hiseq平台方法进行测序
对步骤12得到的纯化产物利用Illumina Hiseq平台进行测序,测序结果的覆盖度和测序深度(测序深度为测序得到的总核苷酸数与待测基因组大小的比值)分别如图3和图4所示,单接头方法的覆盖度为3.24%,测序深度为0.05;双接头方法的覆盖度为2.00%,测序深度为0.02。上述结果表明,利用单接头构建文库的方法优于利用双接头构建,覆盖度和测序深度都得到大幅度提高。另外,结果显示,实验组的覆盖度和测序深度均与阳性对照组无显著差异。
其中,利用双接头构建文库的操作步骤如下:
1.接头(Adapter Mix)制备
a)引物名称及序列:
Primer A:5'-[phos]CTGTCTCTTATACACATCT-3'
Primer B:5'-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-3'
Primer C:5'-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3'
b)使用Annealing Buffer溶解Primer A、Primer B、Primer C至100μM。
c)分别配制如下反应体系:
Figure BDA0001792279460000121
注:可根据实际需要等比放大反应体系。
d)分别将反应1和反应2涡旋震荡充分混匀,并短暂离心使溶液回到管底。置于PCR仪内,进行如下反应程序(Hot lid:105℃):
Hold 温度 时间
1 75℃ 15min
1 60℃ 10min
1 50℃ 10min
1 40℃ 10min
1 25℃ 30min
e)反应结束后,将反应1和反应2的产物等体积混合,混匀。命名为Adapter Mix,-20℃保存。
2.接头(Adapter Mix)包埋
a)在灭菌PCR管中依次添加各反应组分:
成分 体积
Adapter Mix 12.5μl
Tn5转座酶(BGI产品BGE005) 87.5μl
总体积 100μl
b)使用移液器轻轻吹打,充分混匀。
c)将反应置于25℃反应60min。反应产物命名为Tagment Enzyme Advanced MixV5S,置于-20℃存。
3.单细胞裂解液制备
取87μL 30mM Tris-HCl(pH 7.5)和3μL蛋白酶(Qiagen,Cat No./ID:19155QIAGEN Protease(7.5AU))混合即得到的单细胞裂解液,蛋白酶在该溶液中的浓度为0.67mg/mL。
4.单细胞样品制备及裂解
Hela细胞用PBS洗涤去除外部DNA和细胞碎片,使用流式细胞仪分选单个细胞进入灭菌的装有2μL步骤1的单细胞裂解液的PCR管中。以无单细胞的单细胞裂解液作为阴性对照,向单细胞裂解液中添加100个Hela细胞作为阳性对照。50℃1h孵育反应。75℃20min,80℃5min失活蛋白酶,冷却到10℃。
5.利用Tn5酶对单细胞基因组进行片段化。
a)于室温解冻5 X TAG Buffer(华大基因产品,货号为BGE005B01),上下颠倒混匀后备用;确认0.1%SDS处于室温,并轻弹管壁确认有无沉淀。如有沉淀,可于37℃加热并剧烈震荡充分混匀沉淀即会溶解。
b)在上述包含裂解产物的PCR管中依次添加各反应组分:
成分 体积
H<sub>2</sub>0 5.8μL
裂解产物 2μL
5 X TAG Buffer 2μL
Fragment Enzyme Advanced Mix V5S 0.2μL
总体积 10μL
c)55℃温浴30min,添加2.5μL 0.1%SDS水溶液室温5min,终止反应,置于冰上。
6.片段化产物PCR扩增
a)参考引物名称及序列:
Illumina文库接头引物序列:引物序列来源于试剂盒(Nextera DNA LibraryPreparation Kit,FC-121-1030),i5,i7部分为barcode。
Nextera P5 primer:5’-AATGATACGGCGACCACCGAGATCTACAC[i5]TCGTCGGCAGCGTC
Nextera P7 primer:5’-CAAGCAGAAGACGGCATACGAGAT[i7]GTCTCGTGGGCTCGG
b)将灭菌PCR管置于冰浴中,依次添加各反应组分:
成分 体积
灭菌蒸馏水 5.25μl
步骤5产物 12.5μl
5 x KAPA Fidelity Buffer 5μl
10mM each dNTP 0.75μl
Nextera P5 primer(10μM) 0.5μl
Nextera P7 primer(10μM) 0.5μl
KAPA HiFi DNA聚合酶(1U/μl) 0.5μl
总体积 25μl
c)使用移液器轻轻吹打10次充分混匀。
d)将PCR管置于PCR仪中,设置如下反应程序(Hot lid:105℃):
Figure BDA0001792279460000141
7.PCR产物的磁珠纯化
利用1×Ampure XP bead(Beckman Coulter)纯化步骤6得到的PCR产物。
a)起始PCR产物体积应为25μL。因PCR过程中样品挥发会导致产物体积不足25μL,进行下面操作之前必须使用灭菌蒸馏水将体积补齐至25μL。
b)涡旋震荡混匀AMPure XP beads并吸取25μL体积至含有25μL PCR产物的EP管中,使用移液器轻轻吹打10次充分混匀。室温孵育10分钟。
c)步骤b)结束后进行短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心移除上清。
d)保持EP管始终处于磁力架中,步骤c)结束后,向EP管中加入200μL新鲜配制的80%(体积百分比)乙醇水溶液漂洗磁珠。室温孵育30秒后小心移除上清。
e)重复步骤d),总计漂洗两次。
f)保持EP管始终处于磁力架中,开盖,空气中干燥磁珠3分钟。
g)将EP管从磁力架中取出,向EP管中加入21μL灭菌超纯水,涡旋振荡或使用移液器轻轻吹打充分混匀,将EP管短暂离心并置于磁力架中分离磁珠和液体。待溶液澄清(约5分钟)小心吸取上清至灭菌EP管中,即得到纯化后的PCR产物,于-20℃保存。
8.质量与评估
a)质量浓度检测
b)利用Qubit dsDNA HS Assay Kit(Thermo Fisher Scientific)测定,实验组的浓度大于3.5ng/μl。
c)电泳检测
d)采用1%琼脂糖凝胶检测纯化后的PCR产物。
e)Agilent 2100检测
f)将纯化后的PCR产物稀释到合适的浓度,通过Bioanalyzer High SensitivityDNA Kit(Agilent)测定PCR产物片段大小范围。
9.使用Illumina平台方法进行测序。
序列表
<110> 深圳华大生命科学研究院
<120> 一种单细胞DNA文库的构建方法
<160> 3
<170> SIPOSequenceListing 1.0
<210> 1
<211> 34
<212> DNA
<213> 人工序列()
<400> 1
gtctcgtggg ctcggagatg tgtataagag acag 34
<210> 2
<211> 34
<212> DNA
<213> 人工序列()
<400> 2
tacactcttt ccctacacga cgctcttccg atct 34
<210> 3
<211> 33
<212> DNA
<213> 人工序列()
<400> 3
gatcggaaga gcacacgtct gaactccagt cac 33

Claims (10)

1.一种单细胞全基因组文库的构建方法,包括:
利用转座酶复合体打断单细胞基因组DNA,使打断的片段两端的接头相同,得到片段化基因组,利用所述片段化基因组构建单细胞全基因组文库;
所述转座酶复合体包括转座酶和转座酶接头,所述转座酶接头由名称分别为A和C的单链DNA组成;
所述A为所述转座酶的识别序列的反向互补序列;
所述C为含有所述识别序列的单链DNA。
2.根据权利要求1所述的方法,其特征在于:所述C由所述识别序列与名称为c1的单链DNA组成;所述利用所述片段化基因组构建单细胞全基因组文库包括:利用名称为引物1的引物对所述片段化基因组进行扩增,得到PCR产物,所述PCR产物即为单细胞全基因组文库;所述引物含有所述c1。
3.根据权利要求2所述的方法,其特征在于:所述引物1由所述C中所述c1及所述识别序列的部分或全部序列组成。
4.根据权利要求2或3所述的方法,其特征在于:所述c1为序列表中序列1的第1-15位。
5.根据权利要求1-4中任一所述的方法,其特征在于:所述转座酶为Tn5;和/或,所述识别序列为序列表中序列1的第16-34位。
6.根据权利要求2-5中任一所述的方法,其特征在于:所述方法还包括对所述PCR产物连接用于测序的接头步骤。
7.根据权利要求6所述的方法,其特征在于:所述方法还包括:对所述PCR产物连接所述用于测序的接头,得到连接接头的文库;利用识别所述用于测序的接头的引物对所述连接接头的文库进行PCR扩增,得到单细胞全基因组文库。
8.根据权利要求6或7所述的方法,其特征在于:所述用于测序的接头为用于高通量测序的接头。
9.下述X1)-X4)中任一产品:
X1)权利要求1-5中任一所述转座酶复合体;
X2)权利要求1-5中任一所述转座酶接头;
X3)成套试剂,由X31)和X32)组成;
X31)权利要求1-5中任一所述转座酶复合体或所述转座酶接头;
X32)权利要求1-5中任一所述引物1;
X4)成套试剂,由X41)、X42)和X43)组成;
X41)权利要求1-5中任一所述转座酶复合体或所述转座酶接头;
X42)权利要求1-5中任一所述引物1;
X43)权利要求7或8所述识别所述用于测序的接头的引物。
10.下述任一应用:
Ⅰ、权利要求1-8中任一所述的方法在基因组测序中的应用;
Ⅱ、权利要求9所述的产品在基因组测序或构建基因组文库中的应用。
CN201811041751.2A 2018-09-07 2018-09-07 一种单细胞dna文库的构建方法 Active CN110886021B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811041751.2A CN110886021B (zh) 2018-09-07 2018-09-07 一种单细胞dna文库的构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811041751.2A CN110886021B (zh) 2018-09-07 2018-09-07 一种单细胞dna文库的构建方法

Publications (2)

Publication Number Publication Date
CN110886021A true CN110886021A (zh) 2020-03-17
CN110886021B CN110886021B (zh) 2023-08-15

Family

ID=69744444

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811041751.2A Active CN110886021B (zh) 2018-09-07 2018-09-07 一种单细胞dna文库的构建方法

Country Status (1)

Country Link
CN (1) CN110886021B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111893576A (zh) * 2020-07-08 2020-11-06 广州再生医学与健康广东省实验室 一种微量细胞基因组测序文库构建方法
CN113337590A (zh) * 2021-06-03 2021-09-03 深圳华大基因股份有限公司 一种二代测序方法和文库构建方法
CN114150046A (zh) * 2020-09-08 2022-03-08 北京大学 一种低深度高通量单细胞基因变异检测方法
WO2022161294A1 (zh) * 2021-02-01 2022-08-04 南方医科大学 一种中通量单细胞拷贝数文库构建的方法及其应用
CN115386622A (zh) * 2022-10-26 2022-11-25 北京寻因生物科技有限公司 一种转录组文库的建库方法及其应用
CN115948621A (zh) * 2023-01-18 2023-04-11 珠海舒桐医疗科技有限公司 一种基于月经血dna的hpv筛查方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104389026A (zh) * 2014-10-30 2015-03-04 北京诺禾致源生物信息科技有限公司 单细胞转录组测序文库的构建方法及其应用
WO2016058134A1 (zh) * 2014-10-14 2016-04-21 深圳华大基因科技有限公司 一种接头元件和使用其构建测序文库的方法
CN105525357A (zh) * 2014-09-30 2016-04-27 深圳华大基因股份有限公司 一种测序文库的构建方法及试剂盒和应用
CN105524920A (zh) * 2016-01-21 2016-04-27 深圳华大基因研究院 用于Ion Proton测序平台的TN5建库的引物组、试剂盒和建库方法
WO2017084023A1 (zh) * 2015-11-17 2017-05-26 深圳华大基因研究院 一种高通量的单细胞转录组建库方法
CN107586835A (zh) * 2017-10-19 2018-01-16 东南大学 一种基于单链接头的下一代测序文库的构建方法及其应用

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105525357A (zh) * 2014-09-30 2016-04-27 深圳华大基因股份有限公司 一种测序文库的构建方法及试剂盒和应用
WO2016058134A1 (zh) * 2014-10-14 2016-04-21 深圳华大基因科技有限公司 一种接头元件和使用其构建测序文库的方法
CN104389026A (zh) * 2014-10-30 2015-03-04 北京诺禾致源生物信息科技有限公司 单细胞转录组测序文库的构建方法及其应用
WO2017084023A1 (zh) * 2015-11-17 2017-05-26 深圳华大基因研究院 一种高通量的单细胞转录组建库方法
CN105524920A (zh) * 2016-01-21 2016-04-27 深圳华大基因研究院 用于Ion Proton测序平台的TN5建库的引物组、试剂盒和建库方法
CN107586835A (zh) * 2017-10-19 2018-01-16 东南大学 一种基于单链接头的下一代测序文库的构建方法及其应用

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111893576A (zh) * 2020-07-08 2020-11-06 广州再生医学与健康广东省实验室 一种微量细胞基因组测序文库构建方法
CN111893576B (zh) * 2020-07-08 2021-08-17 生物岛实验室 一种微量细胞基因组测序文库构建方法
CN114150046A (zh) * 2020-09-08 2022-03-08 北京大学 一种低深度高通量单细胞基因变异检测方法
WO2022161294A1 (zh) * 2021-02-01 2022-08-04 南方医科大学 一种中通量单细胞拷贝数文库构建的方法及其应用
CN113337590A (zh) * 2021-06-03 2021-09-03 深圳华大基因股份有限公司 一种二代测序方法和文库构建方法
CN115386622A (zh) * 2022-10-26 2022-11-25 北京寻因生物科技有限公司 一种转录组文库的建库方法及其应用
CN115386622B (zh) * 2022-10-26 2023-10-27 北京寻因生物科技有限公司 一种转录组文库的建库方法及其应用
CN115948621A (zh) * 2023-01-18 2023-04-11 珠海舒桐医疗科技有限公司 一种基于月经血dna的hpv筛查方法

Also Published As

Publication number Publication date
CN110886021B (zh) 2023-08-15

Similar Documents

Publication Publication Date Title
CN110886021B (zh) 一种单细胞dna文库的构建方法
CN110268059B (zh) 单细胞全基因组文库及制备其的组合索引方法
US20180080021A1 (en) Simultaneous sequencing of rna and dna from the same sample
JP2020522243A (ja) 核酸のマルチプレックス末端タギング増幅
US20230250476A1 (en) Deep Sequencing Profiling of Tumors
JP7058502B2 (ja) ゲノム適用および治療適用のための、核酸分子のクローン複製および増幅のためのシステムおよび方法
JP2003527867A (ja) ポリヌクレオチド配列改変のマイクロアレイベースの分析
JP2021511794A (ja) 核酸増幅のための方法
EP3728636B1 (en) Particles associated with oligonucleotides
WO2004081224A2 (en) Reca-assisted allele specific oligonucleotide extension method
JP4669614B2 (ja) 多型dnaフラグメントおよびその使用
CN111690748A (zh) 使用高通量测序检测微卫星不稳定的探针组、试剂盒及微卫星不稳定的检测方法
WO2022237035A1 (zh) 一种超高通量单细胞测序方法
CN116635535A (zh) 单细胞dna和rna的同时扩增
CN112680796A (zh) 一种靶标基因富集建库方法
CN114875118B (zh) 确定细胞谱系的方法、试剂盒和装置
US20230095295A1 (en) Phi29 mutants and use thereof
CN117015603A (zh) 使用基于转座子的技术与用于误差校正的独特分子标识符制备定向标签化测序文库的方法
EP3433379B1 (en) Primers with self-complementary sequences for multiple displacement amplification
Chetverin et al. Molecular colony technique: a new tool for biomedical research and clinical practice
US20240316556A1 (en) High-throughput analysis of biomolecules
US20240124930A1 (en) Diagnostic and/or Sequencing Method and Kit
CN117305466B (zh) 一种能够识别单碱基甲基化状态的检测方法
EP3828283A1 (en) An improved sequencing method and kit
WO2023215524A2 (en) Primary template-directed amplification and methods thereof

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant