CN113668068A - 基因组甲基化文库及其制备方法和应用 - Google Patents

基因组甲基化文库及其制备方法和应用 Download PDF

Info

Publication number
CN113668068A
CN113668068A CN202110819010.8A CN202110819010A CN113668068A CN 113668068 A CN113668068 A CN 113668068A CN 202110819010 A CN202110819010 A CN 202110819010A CN 113668068 A CN113668068 A CN 113668068A
Authority
CN
China
Prior art keywords
methylation
linker
sequence
library
transposable
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
CN202110819010.8A
Other languages
English (en)
Other versions
CN113668068B (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.)
Guangzhou Dina Biotechnology Co ltd
Original Assignee
Guangzhou Dina Biotechnology 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 Guangzhou Dina Biotechnology Co ltd filed Critical Guangzhou Dina Biotechnology Co ltd
Priority to CN202110819010.8A priority Critical patent/CN113668068B/zh
Publication of CN113668068A publication Critical patent/CN113668068A/zh
Application granted granted Critical
Publication of CN113668068B publication Critical patent/CN113668068B/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
    • C40B40/00Libraries per se, e.g. arrays, mixtures
    • C40B40/04Libraries containing only organic compounds
    • C40B40/06Libraries containing nucleotides or polynucleotides, or derivatives thereof
    • 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/6806Preparing nucleic acids for analysis, e.g. for polymerase chain reaction [PCR] assay
    • 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
    • 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

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Organic Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Biochemistry (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Microbiology (AREA)
  • General Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Biotechnology (AREA)
  • Immunology (AREA)
  • Biophysics (AREA)
  • Medicinal Chemistry (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进行片段化,得到加上了所述转座接头的片段化产物;使用以mdCTP替代了dCTP的dNTPs对所述片段化产物进行缺口补齐;将经所述缺口补齐后的片段化产物进行转化处理,使未被甲基化修饰的胞嘧啶转化为尿嘧啶;将经所述转化处理后的片段化产物进行PCR扩增,得到所述基因组甲基化文库。本发明可以解决甲基化文库比对分析的问题,实现更准确的甲基化位点分析,可以进行基因突变位点分析,也可大大提高文库碱基平衡性,大幅提高测序质量。

Description

基因组甲基化文库及其制备方法和应用
技术领域
本发明涉及基因检测技术领域,特别是涉及一种基因组甲基化文库及其制备方法和应用。
背景技术
DNA甲基化是DNA化学修饰的一种形式,能够在不改变DNA序列的前提下,改变遗传表现。所谓DNA甲基化是指在DNA甲基化转移酶的作用下,在基因组CpG二核苷酸的胞嘧啶5号碳位共价键结合一个甲基基团。大量研究表明,DNA甲基化能引起染色质结构、DNA构象、DNA稳定性及DNA与蛋白质相互作用方式的改变,从而控制基因表达。甲基化异常是癌症中一种常见的表观遗传学改变,大量研究表明甲基化异常与肿瘤的发生、发展、细胞癌变有着密切的联系。
目前的甲基化研究方法中,主要包含甲基化特异性PCR、高分辨率溶解曲线法、重亚硫酸盐测序法、基因组直接测序法等。最常用的方法是基于重亚硫酸盐的方法,该方法中,在重亚硫酸盐的作用下,DNA中未被甲基化修饰的胞嘧啶(C)转化为尿嘧啶(U),通过高通量测序检测出未被转化的胞嘧啶(C),即为被甲基化修饰的胞嘧啶。随着高通量测序技术的发展,DNA甲基化检测也由单基因水平逐渐向整个基因组水平发展。
全基因组亚硫酸盐测序(whole genome bisulfite sequencing,WGBS)技术是一种基于亚硫酸盐转化的DNA甲基化检测技术。这种传统的WGBS技术是对基因组DNA片段化后,末端修复加A,然后加上已知甲基化修饰序列的接头,通过重亚硫酸盐处理,将未被甲基化修饰的C脱氨基成为U,再进行PCR扩增后转变为T,甲基化修饰的C则不变,最后将PCR得到的产物进行测序。WGBS可绘制单碱基分辨率的全基因组DNA甲基化图谱,但要达到要求的测序深度,测序成本比较高。
考虑到WGBS成本较高,还有一种跟WGBS十分相似,但是只测部分基因组的方法,即RRBS,也可以称为“代表性精简的亚硫酸盐测序”。RRBS主要覆盖启动子区域,可以检测到人类基因组中85%的CpG岛。RRBS法基于酶切原理,利用限制性内切酶MspI对基因组进行酶切富集CpG区域,再经重亚硫酸盐处理后进行高通量测序。该技术显著提高了CpG区域的测序深度,降低了测序成本,但受内切酶位点限制,不能完全覆盖CpG区域。
而且,传统的基于重亚硫酸盐转化的甲基化测序,转化后基因组中大多数的胞嘧啶(C)都会转变为胸腺嘧啶(T),造成碱基不平衡,在生物信息学分析过程中大大降低含有非甲基化胞嘧啶的reads比对率,导致不能得到完整的全基因组甲基化图谱,限制了实际应用。
发明内容
基于此,有必要提供一种可提高测序时含有非甲基化胞嘧啶的reads比对率的基因组甲基化文库的制备方法。
一种基因组甲基化文库的制备方法,包括以下步骤:
提供转座酶复合体,所述转座酶复合体由转座接头和转座酶组成,所述转座接头含有转座子序列;
使用所述转座酶复合体对基因组DNA进行片段化,得到加上了所述转座接头的片段化产物;
使用以mdCTP替代了dCTP的dNTPs对所述片段化产物进行缺口补齐;
将经所述缺口补齐后的片段化产物进行转化处理,使未被甲基化修饰的胞嘧啶转化为尿嘧啶;
将经所述转化处理后的片段化产物进行PCR扩增,得到所述基因组甲基化文库。
本发明通过转座酶复合体对基因组DNA进行片段化并同时加上接头,然后利用mdCTP(甲基化修饰的胞嘧啶)替代了dCTP(未被甲基化修饰的胞嘧啶)的dNTPs对片段化的DNA进行缺口补齐,经转化处理后的文库序列含有原始链序列和转化链序列,从而可以解决甲基化文库比对分析的问题,实现更准确的甲基化位点分析,可以进行基因突变位点分析,也可大大提高文库碱基平衡性,大幅提高测序质量。
在其中一个实施例中,所述转座接头为单链接头、双链接头和茎环接头中的一种或多种。
在其中一个实施例中,所述转座接头含有甲基化修饰的胞嘧啶。
在其中一个实施例中,所述转座接头包括序列如SEQ ID NO:1所示的第一接头和序列如SEQ ID NO:2所示的第二接头,以及序列如SEQ ID NO:3所示的第三接头或序列如SEQ ID NO:4所示的第四接头。
在其中一个实施例中,在所述缺口补齐的步骤之前还包括以下步骤:使用具有3’端到5’端核酸外切功能的酶对所述片段化产物进行酶切处理以扩大所述转座酶复合体形成的缺口。
在其中一个实施例中,所述转座酶选自Tn1、Tn2、Tn3、Tn4、Tn5、Tn6、Tn7、Tn9、Tn10、Tn551、Tn971、Tn916、Tn1545、Tn1681、Tgf2、Tol2、Himar1和HARBI1中的一种或多种。
在其中一个实施例中,所述转化处理包括以下步骤:使用TET2酶和氧化增强剂将甲基化修饰的胞嘧啶氧化成羧基胞嘧啶,然后对未被甲基化修饰的胞嘧啶进行脱氨基处理。
在其中一个实施例中,所述转化处理包括以下步骤:使用亚硫酸氢盐将未被甲基化修饰的胞嘧啶磺化和水解脱氨,生成尿嘧啶磺酸盐中间产物,然后在碱性条件下对所述尿嘧啶磺酸盐中间产物进行脱磺酸基处理。
本发明还提供了一种基因组甲基化文库,其根据如上所述的制备方法制备得到。
本发明还提供了一种测序方法,包括以下步骤:
根据如上所述的制备方法制备得到基因组甲基化文库;
对所述基因组甲基化文库进行双端测序,获得双端测序数据;
将所述双端测序数据与参考基因组数据进行比对,测序reads中通过所述缺口补齐形成的部分用于序列定位,测序reads中其他部分用于判断甲基化位点,根据比对到参考基因组序列上的信息,获得所述基因组DNA的甲基化位点。
附图说明
图1为本发明一实施例的转座酶复合体的工作原理示意图;
图2为本发明实施例1的实验流程示意图;
图3为本发明实施例2的实验流程示意图;
图4为本发明实施例3的实验流程示意图;
图5为本发明一实施例的序列比对原理示意图。
具体实施方式
为了便于理解本发明,下面将对本发明进行更全面的描述,并给出了本发明的较佳实施例。但是,本发明可以以许多不同的形式来实现,并不限于本文所描述的实施例。相反地,提供这些实施例的目的是使对本发明的公开内容的理解更加透彻全面。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。本文所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
本发明一实施例的基因组甲基化文库的制备方法,包括以下步骤S1~S5:
S1、提供转座酶复合体(transposome),转座酶复合体由转座接头和转座酶组成,转座接头含有转座子序列;
S2、使用转座酶复合体对基因组DNA(genomic DNA)进行片段化(tagmentation),得到加上了转座接头且具有缺口的片段化产物;
S3、使用以mdCTP替代了dCTP的dNTPs对片段化产物进行缺口补齐(gap filling);
S4、将经缺口补齐后的片段化产物进行转化处理(convert),使未被甲基化修饰的胞嘧啶转化为尿嘧啶;
S5、将经转化处理后的片段化产物进行PCR扩增(amplification),得到基因组甲基化文库。
本发明通过转座酶复合体对基因组DNA进行片段化并同时加上接头,然后利用mdCTP(甲基化修饰的胞嘧啶)替代了dCTP(未被甲基化修饰的胞嘧啶)的dNTPs对片段化的DNA进行缺口补齐,经转化处理后的文库序列含有原始链序列和转化链序列,从而可以解决甲基化文库比对分析的问题,实现更准确的甲基化位点分析,可以进行基因突变位点分析,也可大大提高文库碱基平衡性,大幅提高测序质量。
在一个具体示例中,转座接头为单链接头、双链接头和茎环接头中的一种或多种。
在一个具体示例中,转座接头含有甲基化修饰的胞嘧啶,从而可以避免转座接头在后续的转化处理中序列发生巨大变化,提高比对分析的准确性。可选地,转座接头还含有单分子标签序列,从而便于区分不同样本。
在一个具体示例中,转座接头包括序列如SEQ ID NO:1所示的第一接头和序列如SEQ ID NO:2所示的第二接头,以及序列如SEQ ID NO:3所示的第三接头或序列如SEQ IDNO:4所示的第四接头。采用第一接头、第二接头和第四接头时,制备的基因组甲基化文库是将两条互补链连接在一起,可以自身比对,不需要借助参考序列,序列信息更全面。
在一个具体示例中,在缺口补齐的步骤之前还包括以下步骤:使用具有3’端到5’端核酸外切功能的酶对片段化产物进行酶切处理以扩大转座酶复合体形成的缺口。例如Tn5转座酶复合体处理后会形成9bp的缺口,可以通过酶切处理形成95~105bp左右的缺口。如此,通过缺口补齐形成的部分更长,则未被转化的原始链序列更长,能够更加准确地与参考基因组进行比对分析从而定位其位置。
在一个具体示例中,上述酶切处理的步骤包括:在片段化产物中加入T4DNA聚合酶(polymerase),于10℃~14℃酶切40~80min,然后在70℃~74℃孵育5~15min。
在一个具体示例中,转座酶选自Tn1、Tn2、Tn3、Tn4、Tn5、Tn6、Tn7、Tn9、Tn10、Tn551、Tn971、Tn916、Tn1545、Tn1681、Tgf2、Tol2、Himar1和HARBI1中的一种或多种,但不限于此。
在一个具体示例中,转化处理采用酶法转化,例如可包括以下步骤:使用TET2酶和氧化增强剂将甲基化修饰的胞嘧啶氧化成羧基胞嘧啶,保护经过修饰的胞嘧啶,使其免受下游的脱氨基作用,然后使用脱氨酶例如APOBEC等对未被甲基化修饰的胞嘧啶进行脱氨基处理,但不影响羧基胞嘧啶。
在一个具体示例中,转化处理采用重亚硫酸盐转化,例如可包括以下步骤:使用亚硫酸氢盐将未被甲基化修饰的胞嘧啶磺化和水解脱氨,生成尿嘧啶磺酸盐中间产物,然后在碱性条件下对所述尿嘧啶磺酸盐中间产物进行脱磺酸基处理,但并不影响甲基化修饰的胞嘧啶。
在一个具体示例中,上述PCR扩增时所用的扩增引物含有测序接头序列和/或标签序列。可选地,制备方法还包括以下步骤:对上述PCR扩增所得的产物进行纯化(purification)。可选地,纯化方法为膜柱子法或磁珠法等。
本发明一实施例的基因组甲基化文库,其根据如上所述的制备方法制备得到。
本发明一实施例的测序方法,包括以下步骤S1~S3:
S1、根据如上所述的制备方法制备得到基因组甲基化文库;
S2、对基因组甲基化文库进行双端测序,获得双端测序数据;
S3、将双端测序数据与参考基因组数据进行比对,测序reads中通过缺口补齐形成的部分用于序列定位,测序reads中其他部分用于判断甲基化位点,根据比对到参考基因组序列上的信息,获得基因组DNA的甲基化位点。通过缺口补齐形成的部分与原核酸序列相同,用于序列定位,测序reads中其他部分中未被甲基化修饰的胞嘧啶已转变为胸腺嘧啶,用于判断甲基化位点。
在一个具体示例中,在进行比对之前还包括以下步骤:去除双端测序数据中的接头、测序引物(sequencing primer)和低质量数据,获得干净的双端测序数据以用于比对,从而筛选得到高质量的比对序列。
在一个具体示例中,在进行双端测序之前还包括以下步骤:对基因组甲基化文库进行目的基因片段捕获和纯化,从而获得目的基因的甲基化文库,如此以用于靶向基因甲基化测序。
下面主要结合具体实施方式和附图对本发明作进一步详细的说明。
实施例1采用Tn5转座酶介导的基因组甲基化文库构建
(1)样本准备:取HEK293细胞,使用细胞基因组DNA提取试剂盒进行DNA提取。所得的核酸分别用50μL TE buffer溶解,并使用Qubit 3.0荧光计进行浓度测定。最后取50ngDNA按照下述的操作步骤进行文库构建。
(2)Tn5处理标记:
a、合成含有19bp转座酶识别序列的核酸片段(详细序列如下述SEQ ID NO:1、SEQID NO:2和SEQ ID NO:3所示),其中序列B和序列C中部分胞嘧啶为甲基化修饰的胞嘧啶(C)。序列B和序列C中,“N”代表任意碱基A、T、C或者G,其中胞嘧啶为甲基化修饰的胞嘧啶。
接头序列A:CTGTCTCTTATACACATCT(SEQ ID NO:1);
接头序列B:TCGTCGGCAGCGTCNNNNNNAGATGTGTATAAGAGACAG(SEQ ID NO:2);
接头序列C:TCTCGTGGGCTCGGNNNNNNAGATGTGTATAAGAGACAG(SEQ ID NO:3)。
b、将接头序列A、B和C分别稀释到100μM,序列A+序列B按1:1混合,序列A+序列C按1:1混合,充分混合后离心,于PCR仪中按如下程序(表1)退火得到转座接头(储存于-20℃),用于转座酶复合体的制备。
表1
Figure BDA0003171195200000071
c、按照如下体系(表2)将转座接头与转座酶包埋成转座酶复合体,轻轻吹打20次混合,30℃孵育1小时后完成复合体包埋,该转座酶复合体储存于-20℃。
表2
组分 体积(μL)
转座酶 85μL
转座接头 30μL
偶联缓冲液 85μL
合计 200μL
d、按照如下体系(表3)将50ng的高质量基因组DNA和转座酶复合体进行混合,轻轻吹打20次混合,55℃孵育10分钟后降温至4℃,完成基因组的打断,其工作原理如图1所示。
表3
组分 体积(μL)
5μL
5×打断缓冲液 2μL
gDNA(50ng/μL) 1μL
转座酶复合体 2μL
合计 10μL
e、终止反应:反应完成后立即向产物中加入5μL 5×终止液,使用移液器轻轻吹打充分混匀,置于室温放置5min。
(3)修补缺口:经Tn5处理后会形成9bp的缺口,向上述产物中加入15μL的PCRMaster Mix(mdCTP替代dCTP的dNTP),然后在72℃孵育10min。
(4)转化:使用QIAGEN公司的
Figure BDA0003171195200000081
Fast DNA Bisulfite Kit处理上述产物,将30μL产物加入110μL转换试剂中(85μL Bisulfite Solution,25μL DNA ProtectBuffer),瞬离后放入PCR仪,设定反应程序如下:95℃,5min;60℃,20min;95℃,5min;60℃,20min;20℃,∞;反应结束后将产物转移至1.5mL离心管中,加入310μL Buffer BL和250μL无水乙醇,充分混匀后瞬时离线,将全部液体转移至膜柱子
Figure BDA0003171195200000082
DNA spin columns上离心1min,加入500μL清洗液Buffer BW后离心1min,加入500μL脱磺液Buffer BD后室温(20~30℃)静置15~20min,离心1min,加入500μL清洗液Buffer BW离心1min,重复清洗一次;加入250μL无水乙醇离心1min,开盖后于60℃放置5min以去除残余液体,加入25μL洗脱液Buffer EB后室温(20~30℃)静置1min,离心1min收集洗脱液。
(5)扩增:按照下表4进行反应体系的配制:
表4
组分 体积(μL)
2×HiFidelity PCR Mix 25μL
处理后的DNA样品 23μL
Primer N5 1μL
Primer N7 1μL
合计 50μL
涡旋混匀后瞬时离心,在PCR仪中按如下表5程序进行反应:
表5
Figure BDA0003171195200000091
(6)纯化:使用Beckman公司的AMPure XP磁珠进行文库片段分选,选用的分选比例为第一轮0.7×,第二轮为0.2×,所得文库大小为400bp左右。文库质检:文库通过Qubit3.0核酸定量仪定量,并进行琼脂糖凝胶电泳查看文库分布。本实施例的实验流程如图2所示。
实施例2Tn5转座酶介导和缺口酶切处理的基因组甲基化文库构建
(1)样本准备:取HEK293细胞,使用细胞基因组DNA提取试剂盒进行DNA提取。所得的核酸分别用50μL TE buffer溶解,并使用Qubit 3.0荧光计进行浓度测定。最后取50ngDNA按照下述的操作步骤进行文库构建。
(2)Tn5处理标记:
a、合成含有19bp转座酶识别序列的核酸片段(详细序列如下述SEQ ID NO:1、SEQID NO:2和SEQ ID NO:3所示),其中序列B和序列C中部分胞嘧啶为甲基化修饰的胞嘧啶(C)。序列B和序列C中,“N”代表任意碱基A、T、C或者G,其中胞嘧啶为甲基化修饰的胞嘧啶。
接头序列A:CTGTCTCTTATACACATCT(SEQ ID NO:1);
接头序列B:TCGTCGGCAGCGTCNNNNNNAGATGTGTATAAGAGACAG(SEQ ID NO:2);
接头序列C:TCTCGTGGGCTCGGNNNNNNAGATGTGTATAAGAGACAG(SEQ ID NO:3)。
b、将接头序列A、B和C分别稀释到100μM,序列A+序列B按1:1混合,序列A+序列C按1:1混合,充分混合后离心,于PCR仪中按表1退火得到转座接头(储存于-20℃),用于转座酶复合体的制备。
c、按照表2将转座接头与转座酶包埋成转座酶复合体,轻轻吹打20次混合,30℃孵育1小时后完成复合体包埋,该转座酶复合体储存于-20℃。
d、按照表3将50ng的高质量基因组DNA和转座酶复合体进行混合,轻轻吹打20次混合,55℃孵育10分钟后降温至4℃,完成基因组的打断。
e、终止反应:反应完成后立即向产物中加入5μL终止液,使用移液器轻轻吹打充分混匀,置于室温放置5min。
(3)核酸外切酶处理:向上述产物中加入1μL的T4 DNA Polymerase(5U/μL)和4μL的5×Reaction Buffer,充分混匀,瞬离离心,12℃酶切60min,然后在72℃孵育10min。
(4)修补缺口:经酶切处理后会形成100bp左右的缺口,向上述产物中加入终浓度为10mM的dNTP(其中mdCTP替代dCTP),充分混匀,瞬离,37℃反应15min,70℃孵育10分钟终止反应。
(5)使用QIAGEN公司的
Figure BDA0003171195200000101
Fast DNA Bisulfite Kit处理上述产物,将上述产物补水至30μL加入110μL转换试剂中(85μL Bisulfite Solution,25μLDNA ProtectBuffer),瞬离后放入PCR仪,设定反应程序如下:95℃,5min;60℃,20min;95℃,5min;60℃,20min;20℃,∞;反应结束后将产物转移至1.5mL离心管中,加入310μL Buffer BL和250μL无水乙醇,充分混匀后瞬时离线,将全部液体转移至膜柱子
Figure BDA0003171195200000111
DNA spin columns上离心1min,加入500μL清洗液Buffer BW后离心1min,加入500μL脱磺液Buffer BD后室温(20~30℃)静置15~20min,离心1min,加入500μL清洗液Buffer BW离心1min,重复清洗一次;加入250μL无水乙醇离心1min,开盖后于60℃放置5min以去除残余液体,加入25μL洗脱液Buffer EB后室温(20~30℃)静置1min,离心1min收集洗脱液。
(6)扩增:按照表4进行反应体系的配制,涡旋混匀后瞬时离心,在PCR仪中按如表5程序进行反应。
(7)纯化:使用Beckman公司的AMPure XP磁珠进行文库片段分选,选用的分选比例为第一轮0.7×,第二轮为0.2×,所得文库大小为400bp左右。文库质检:文库通过Qubit3.0核酸定量仪定量,并进行琼脂糖凝胶电泳查看文库分布。本实施例的实验流程如图3所示。
实施例3不同接头组合及Tn5转座酶介导的基因组甲基化文库构建
(1)样本准备:取HEK293细胞,使用细胞基因组DNA提取试剂盒进行DNA提取。所得的核酸分别用50μL TE buffer溶解,并使用Qubit 3.0荧光计进行浓度测定。最后取50ngDNA按照下述的操作步骤进行文库构建。
(2)Tn5处理标记:
a、合成含有19bp转座酶识别序列的核酸片段(详细序列如下述SEQ ID NO:1、SEQID NO:2和SEQ ID NO:4所示),其中序列B中部分胞嘧啶为甲基化修饰胞嘧啶(C)。序列B和序列D中,“N”代表任意碱基A、T、C或者G,其中胞嘧啶为甲基化修饰的胞嘧啶。
接头序列A:CTGTCTCTTATACACATCT(SEQ ID NO:1);
接头序列B:TCGTCGGCAGCGTCNNNNNNAGATGTGTATAAGAGACAG(SEQ ID NO:2);
接头序列D:
CTGTCTCTTATACACATCTNNNNNNNNAGATGTGTATAAGAGACAG(SEQ ID NO.4)。
b、将接头序列A、B和D分别稀释到100μM,序列A+序列B按1:1混合,序列D可自身互补形成茎环,充分混匀后离心,于PCR仪中按表1退火得到转座接头(储存于-20℃),用于转座酶复合体的制备。
c、按照表2将转座接头与转座酶包埋成转座酶复合体,轻轻吹打20次混合,30℃孵育1小时后完成复合体包埋,该复合体储存于-20℃。
d、按照表3将50ng的高质量基因组DNA和转座酶复合体进行混合,轻轻吹打20次混合,55℃孵育10分钟后降温至4℃,完成基因组的打断。
e、终止反应:反应完成后立即向产物中加入5μL 5×终止液,使用移液器轻轻吹打充分混匀,置于室温放置5min。
(3)修补缺口:向上述产物中加入1μL T4 DNA聚合酶(5U/μL)、4μL的5×ReactionBuffer、终浓度为10mM的dNTP(其中mdCTP替代dCTP)37℃反应15min,70℃孵育10min后,加入1μL T4 DNA链接酶(5U)和5μL 5×Buffer,然后在22℃孵育10min,然后在72℃孵育5min。
(4)转化:使用QIAGEN公司的
Figure BDA0003171195200000121
Fast DNA Bisulfite Kit处理上述产物,将上述产物补水至30μL加入110μL转换试剂中(85μL Bisulfite Solution,25μL DNAProtect Buffer),瞬离后放入PCR仪,设定反应程序如下:95℃,5min;60℃,20min;95℃,5min;60℃,20min;20℃,∞;反应结束后将产物转移至1.5mL离心管中,加入310μL BufferBL和250μL无水乙醇,充分混匀后瞬时离线,将全部液体转移至膜柱子
Figure BDA0003171195200000122
DNAspin columns上离心1min,加入500μL清洗液Buffer BW后离心1min,加入500μL脱磺液Buffer BD后室温(20~30℃)静置15~20min,离心1min,加入500μL清洗液Buffer BW离心1min,重复清洗一次;加入250μL无水乙醇离心1min,开盖后于60℃放置5min以去除残余液体,加入25μL洗脱液Buffer EB后室温(20~30℃)静置1min,离心1min收集洗脱液。
(5)扩增:按照表4进行反应体系的配制,涡旋混匀后瞬时离心,在PCR仪中按如表5程序进行反应。
(6)纯化:使用Beckman公司的AMPure XP磁珠进行文库片段分选,选用的分选比例为第一轮0.7×,第二轮为0.2×,所得文库大小为400bp左右。文库质检:文库通过Qubit3.0核酸定量仪定量,并进行琼脂糖凝胶电泳查看文库分布。本实施例的实验流程如图4所示。
实施例4甲基化文库在基因组甲基化测序中的应用
(1)文库构建:按照实施例1~3中的建库流程分别进行基因组甲基化文库建库。
(2)高通量测序:将质检好的文库分别在Illumina测序仪上进行高通量测序,获得文库序列信息。
(3)数据分析:实施例1的文库所得测序reads的3’端有9bp的原始链序列信息,可与参考基因组进行比对分析,根据3’端定位reads在基因组的位置,从而确定reads甲基化位点信息,序列比对原理如图5所示;
实施例2的文库所得测序reads的3’端有100bp左右的原始链序列信息,可与参考基因组进行比对分析,根据3’端定位reads在基因组的位置,从而确定reads甲基化位点信息;
实施例3的文库所得测序reads的中间段有接头序列信息,接头两端序列为原互补原始链序列,可自身进行比对分析,从而确定甲基化位点信息。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
序列表
<110> 广州滴纳生物科技有限公司
<120> 基因组甲基化文库及其制备方法和应用
<160> 4
<170> SIPOSequenceListing 1.0
<210> 1
<211> 19
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 1
ctgtctctta tacacatct 19
<210> 2
<211> 39
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 2
tcgtcggcag cgtcnnnnnn agatgtgtat aagagacag 39
<210> 3
<211> 39
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 3
tctcgtgggc tcggnnnnnn agatgtgtat aagagacag 39
<210> 4
<211> 46
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 4
ctgtctctta tacacatctn nnnnnnnaga tgtgtataag agacag 46

Claims (10)

1.一种基因组甲基化文库的制备方法,其特征在于,包括以下步骤:
提供转座酶复合体,所述转座酶复合体由转座接头和转座酶组成,所述转座接头含有转座子序列;
使用所述转座酶复合体对基因组DNA进行片段化,得到加上了所述转座接头的片段化产物;
使用以mdCTP替代了dCTP的dNTPs对所述片段化产物进行缺口补齐;
将经所述缺口补齐后的片段化产物进行转化处理,使未被甲基化修饰的胞嘧啶转化为尿嘧啶;
将经所述转化处理后的片段化产物进行PCR扩增,得到所述基因组甲基化文库。
2.根据权利要求1所述的制备方法,其特征在于,所述转座接头为单链接头、双链接头和茎环接头中的一种或多种。
3.根据权利要求1所述的制备方法,其特征在于,所述转座接头含有甲基化修饰的胞嘧啶。
4.根据权利要求1所述的制备方法,其特征在于,所述转座接头包括序列如SEQ ID NO:1所示的第一接头和序列如SEQ ID NO:2所示的第二接头,以及序列如SEQ ID NO:3所示的第三接头或序列如SEQ ID NO:4所示的第四接头。
5.根据权利要求1~4任一项所述的制备方法,其特征在于,在所述缺口补齐的步骤之前还包括以下步骤:使用具有3’端到5’端核酸外切功能的酶对所述片段化产物进行酶切处理以扩大所述转座酶复合体形成的缺口。
6.根据权利要求1~4所述的制备方法,其特征在于,所述转座酶选自Tn1、Tn2、Tn3、Tn4、Tn5、Tn6、Tn7、Tn9、Tn10、Tn551、Tn971、Tn916、Tn1545、Tn1681、Tgf2、Tol2、Himar1和HARBI1中的一种或多种。
7.根据权利要求1~4任一项所述的制备方法,其特征在于,所述转化处理包括以下步骤:使用TET2酶和氧化增强剂将甲基化修饰的胞嘧啶氧化成羧基胞嘧啶,然后对未被甲基化修饰的胞嘧啶进行脱氨基处理。
8.根据权利要求1~4任一项所述的制备方法,其特征在于,所述转化处理包括以下步骤:使用亚硫酸氢盐将未被甲基化修饰的胞嘧啶磺化和水解脱氨,生成尿嘧啶磺酸盐中间产物,然后在碱性条件下对所述尿嘧啶磺酸盐中间产物进行脱磺酸基处理。
9.一种基因组甲基化文库,其特征在于,根据权利要求1~8任一项所述的制备方法制备得到。
10.一种测序方法,其特征在于,包括以下步骤:
根据权利要求1~8任一项所述的制备方法制备得到基因组甲基化文库;
对所述基因组甲基化文库进行双端测序,获得双端测序数据;
将所述双端测序数据与参考基因组数据进行比对,测序reads中通过所述缺口补齐形成的部分用于序列定位,测序reads中其他部分用于判断甲基化位点,根据比对到参考基因组序列上的信息,获得所述基因组DNA的甲基化位点。
CN202110819010.8A 2021-07-20 2021-07-20 基因组甲基化文库及其制备方法和应用 Active CN113668068B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110819010.8A CN113668068B (zh) 2021-07-20 2021-07-20 基因组甲基化文库及其制备方法和应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110819010.8A CN113668068B (zh) 2021-07-20 2021-07-20 基因组甲基化文库及其制备方法和应用

Publications (2)

Publication Number Publication Date
CN113668068A true CN113668068A (zh) 2021-11-19
CN113668068B CN113668068B (zh) 2024-07-19

Family

ID=78539635

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110819010.8A Active CN113668068B (zh) 2021-07-20 2021-07-20 基因组甲基化文库及其制备方法和应用

Country Status (1)

Country Link
CN (1) CN113668068B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113969307A (zh) * 2021-11-24 2022-01-25 竹石生物科技(苏州)有限公司 Dna甲基化测序文库及制备方法和dna甲基化检测方法
CN114411266A (zh) * 2022-02-22 2022-04-29 深圳华大基因股份有限公司 一种基于启动子构建文库的方法及其文库
CN115386966A (zh) * 2022-10-26 2022-11-25 北京寻因生物科技有限公司 Dna表观修饰的建库方法、测序方法及其建库试剂盒
WO2023130670A1 (zh) * 2022-01-07 2023-07-13 无锡臻和生物科技有限公司 基于游离dna的基因组癌变信息检测系统和检测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1829524A (zh) * 2003-05-08 2006-09-06 波特·W·安德森 基于gp160和人cd4蛋白共有的保守氨基酸序列的抗hiv-1化合物
CN101855242A (zh) * 2007-09-14 2010-10-06 阿迪马布公司 合理设计的合成抗体文库及其用途
WO2013090588A1 (en) * 2011-12-13 2013-06-20 Oslo Universitetssykehus Hf Methods and kits for detection of methylation status
CN108300766A (zh) * 2018-01-16 2018-07-20 四川大学 利用转座酶对染色质开放区和线粒体甲基化研究的方法
CN108753771A (zh) * 2018-06-20 2018-11-06 深圳市海普洛斯生物科技有限公司 一种Hemi-M甲基化修饰引物及其应用

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1829524A (zh) * 2003-05-08 2006-09-06 波特·W·安德森 基于gp160和人cd4蛋白共有的保守氨基酸序列的抗hiv-1化合物
CN101855242A (zh) * 2007-09-14 2010-10-06 阿迪马布公司 合理设计的合成抗体文库及其用途
WO2013090588A1 (en) * 2011-12-13 2013-06-20 Oslo Universitetssykehus Hf Methods and kits for detection of methylation status
CN108300766A (zh) * 2018-01-16 2018-07-20 四川大学 利用转座酶对染色质开放区和线粒体甲基化研究的方法
CN108753771A (zh) * 2018-06-20 2018-11-06 深圳市海普洛斯生物科技有限公司 一种Hemi-M甲基化修饰引物及其应用

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113969307A (zh) * 2021-11-24 2022-01-25 竹石生物科技(苏州)有限公司 Dna甲基化测序文库及制备方法和dna甲基化检测方法
WO2023130670A1 (zh) * 2022-01-07 2023-07-13 无锡臻和生物科技有限公司 基于游离dna的基因组癌变信息检测系统和检测方法
CN114411266A (zh) * 2022-02-22 2022-04-29 深圳华大基因股份有限公司 一种基于启动子构建文库的方法及其文库
CN115386966A (zh) * 2022-10-26 2022-11-25 北京寻因生物科技有限公司 Dna表观修饰的建库方法、测序方法及其建库试剂盒
CN115386966B (zh) * 2022-10-26 2023-03-21 北京寻因生物科技有限公司 Dna表观修饰的建库方法、测序方法及其建库试剂盒

Also Published As

Publication number Publication date
CN113668068B (zh) 2024-07-19

Similar Documents

Publication Publication Date Title
CN113668068B (zh) 基因组甲基化文库及其制备方法和应用
CN108103055B (zh) 一种单细胞rna逆转录与文库构建的方法
JP6438126B2 (ja) 核酸一本鎖環状ライブラリを構築するための方法および試薬キット
US9518295B2 (en) High-throughput sequencing method for methylated DNA and use thereof
CN103088433B (zh) 全基因组甲基化高通量测序文库的构建方法及其应用
EP2725125B1 (en) High throughput methylation detection method
CN108138364B (zh) 一种核酸单链环状文库的构建方法和试剂
CN111808854B (zh) 带有分子条码的平衡接头及快速构建转录组文库的方法
CN109593757B (zh) 一种探针及其适用于高通量测序的对目标区域进行富集的方法
CN106554955B (zh) 构建pkhd1基因突变的测序文库的方法和试剂盒及其用途
CN107937985A (zh) 一种微量碎片化dna甲基化检测文库的构建方法和检测方法
CN108070658B (zh) 检测msi的非诊断方法
CN107858409B (zh) 一种微量降解基因组dna甲基化建库测序方法及其试剂盒
CN107475779A (zh) 适用于单细胞rrbs测序的文库建立方法及其应用
JP2021521786A (ja) 核酸を結合、修飾、および切断する物質の基質選択性および部位のためのin vitroでの高感度アッセイ
CN108103164B (zh) 一种利用多重荧光竞争性pcr检测拷贝数变异的方法
CN106011230A (zh) 用于检测碎片化dna目标区域的引物组合物及其应用
CN111748637A (zh) 一种用于亲缘关系分析鉴定的snp分子标记组合、多重复合扩增引物组、试剂盒及方法
CN110699425B (zh) 基因目标区域的富集方法及体系
WO2022199242A1 (zh) 一组条码接头以及中通量多重单细胞代表性dna甲基化建库和测序方法
CN115715323A (zh) 一种高兼容性的PCR-free建库和测序方法
CN109295500B (zh) 一种单细胞甲基化测序技术及其应用
CN112080555A (zh) Dna甲基化检测试剂盒及检测方法
CN113215663B (zh) 一种基于高通量测序的胃癌靶向治疗基因组文库的构建方法及引物
CN114808148A (zh) 一种dna文库构建试剂盒、文库构建方法和应用

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