CN112251422B - 含独特分子标签序列的转座酶复合体及其应用 - Google Patents

含独特分子标签序列的转座酶复合体及其应用 Download PDF

Info

Publication number
CN112251422B
CN112251422B CN202011129375.XA CN202011129375A CN112251422B CN 112251422 B CN112251422 B CN 112251422B CN 202011129375 A CN202011129375 A CN 202011129375A CN 112251422 B CN112251422 B CN 112251422B
Authority
CN
China
Prior art keywords
sequence
umi
sequencing
transposase
sequences
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.)
Active
Application number
CN202011129375.XA
Other languages
English (en)
Other versions
CN112251422A (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.)
Huazhong Agricultural University
Original Assignee
Huazhong Agricultural University
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 Huazhong Agricultural University filed Critical Huazhong Agricultural University
Priority to CN202011129375.XA priority Critical patent/CN112251422B/zh
Publication of CN112251422A publication Critical patent/CN112251422A/zh
Application granted granted Critical
Publication of CN112251422B publication Critical patent/CN112251422B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N9/00—Enzymes; Proenzymes; Compositions thereof; Processes for preparing, activating, inhibiting, separating or purifying enzymes
    • C12N9/10—Transferases (2.)
    • C12N9/12—Transferases (2.) transferring phosphorus containing groups, e.g. kinases (2.7)
    • C12N9/1241—Nucleotidyltransferases (2.7.7)
    • C—CHEMISTRY; METALLURGY
    • C07—ORGANIC CHEMISTRY
    • C07K—PEPTIDES
    • C07K19/00—Hybrid peptides, i.e. peptides covalently bound to nucleic acids, or non-covalently bound protein-protein complexes
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C12N15/09—Recombinant DNA-technology
    • C12N15/10—Processes for the isolation, preparation or purification of DNA or RNA
    • C12N15/1034—Isolating an individual clone by screening libraries
    • C12N15/1093—General methods of preparing gene libraries, not provided for in other subgroups
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6806—Preparing nucleic acids for analysis, e.g. for polymerase chain reaction [PCR] assay
    • C—CHEMISTRY; METALLURGY
    • C40—COMBINATORIAL TECHNOLOGY
    • C40B—COMBINATORIAL CHEMISTRY; LIBRARIES, e.g. CHEMICAL LIBRARIES
    • C40B50/00—Methods of creating libraries, e.g. combinatorial synthesis
    • C40B50/06—Biochemical methods, e.g. using enzymes or whole viable microorganisms
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
    • Y02A50/30—Against vector-borne diseases, e.g. mosquito-borne, fly-borne, tick-borne or waterborne diseases whose impact is exacerbated by climate change

Landscapes

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

Abstract

本发明公开了一种含独特分子标签序列的转座酶复合体及其应用，所述转座酶复合体为由寡聚核苷酸接头和转座酶组装形成的复合体；利用本发明所述的含有UMI标签的转座酶复合体构建的二代测序文库，能够有效区分转座酶插入重复事件和PCR导致的重复读段，提高了测序数据的利用率、序列变异及其频率鉴定的可靠性、基因表达和染色质可及性定量的准确性，以及鉴定DNA结合蛋白的足迹的敏感性。

Description

含独特分子标签序列的转座酶复合体及其应用
技术领域
本发明涉及测序技术领域，具体涉及一种含独特分子标签序列(UMI)的转座酶复合体及其应用。
背景技术
新一代测序(Next Generation Sequencing,NGS)技术深刻影响了生物学研究和临床诊断。构建测序文库是进行NGS的基础。从2010年开始，转座酶Tn5被广泛用于构建基于Illumina公司的测序平台的DNA测序文库。其原理是：Tn5是一种细菌转座子，最早在大肠杆菌中发现。经改造后，Tn5转座酶可以与长度为19bp的ME(mosaic end)双链序列结合，并将其近乎随机地插入双链DNA大片段中。因此，只要将接头与ME序列一起合成，然后与Tn5转座酶孵育，形成转座复合体，再在适宜的体系和条件下与DNA样本孵育，就可一步实现DNA样本片段化和接头连接，极大地减少了所需样本量，提高了建库的效率。
2013年，Jason D Buenrostro等提出ATAC-seq(assay for transposase-accessible chromatin using sequencing)方法，将转座酶Tn5与真核生物的细胞核一起孵育，通过检测Tn5插入基因组各区段的频率，可以评估染色质开放程度(即“染色质可及性”)。该方法具有快速、灵敏的特点，已被广泛用于表观基因组研究。在ATAC-seq结果中，Tn5插入事件的密度反映了一个区域的染色质的可及性；而在一个高度可及的区域中，一些插入事件突然下降的小区域可能表征了结合有转录因子，称为“足迹(Footprint)”。
2020年，Lin Di等和Bo Lu等相继发现转座酶Tn5可插入RNA/DNA杂合双链，提出其可用于快速构建RNA测序文库。
与大多数NGS建库测序的方法一样，上述基于Tn5转座酶的应用，均需要使用PCR扩增，以确保获得足够的DNA片段在Illumina平台上测序。由于不同片段的扩增效率不同，在最终的测序文库中可能某些片段的数量会过多，为后续的一些分析带来偏差。为了消除这种偏差，我们通常将与基因组比对后具有相同的基因组坐标的读段(read)认为是PCR引起的重复拷贝，并在各种分析之前将其去除。然而，这种重复数据删除模式在某些情况下是有缺陷的。例如，随着测序深度的增加，不同的Tn5插入事件产生相同片段的概率会增加。对于小基因组或具有许多串联重复的物种，这种现象会加剧，因为它们在PCR扩增之前就产生相同片段的概率要高得多。在这种情况下，基于基因组坐标删除PCR重复将会误删有信息的片段，也会引入偏差。
UMI(unique molecular identifiers)是一种分子标签，由随机合成的序列组成，常用来对mRNA转录本进行准确定量。在单细胞RNA-seq中，mRNA文库通过逆转录过程，在5’端加入带有UMI的引物标签，再进行PCR扩增和测序。测序结果中不同的mRNA转录本将带有不同的UMI标签，而通过PCR产生的拷贝，将带有相同的UMI标签。因此，在生物信息分析过程中，UMI作为内部验证对照，可以有效去除由PCR引起的重复拷贝，从而做到准确定量。目前，尚未有将UMI技术与Tn5转座酶建库技术结合的报道。
最后，构建好测序文库后，上机测序时，需要选择合适的、能跟文库片段匹配的测序引物。目前Illumina测序平台的通用测序引物有两种类型，分别叫TruSeq和Nextera。其中基于转座酶Tn5构建的文库，使用的测序引物是Nextera，其主要序列是转座酶结合序列ME，其原因是ME序列在Tn5转座时插入DNA样本，其3′与DNA样本相连。因此，测序时将从ME序列3′开始读取信号。为了不影响转座反应，只能将UMI与ME的5′相连，此时用Nextera测序引物将无法获得UMI序列。另一方面，只有样本建库时包含与通用测序引物匹配的序列，才能方便与不同研究者的测序文库混合，一起进行商业化测序，降低成本、提高效率；否则，需要使用自定义测序引物，承包流动槽(lane)测序，存在灵活性和成本的问题。因此，精巧地设计Tn5转座酶复合体的接头序列，是实现UMI技术与Tn5转座酶建库技术结合并商业化推广应用的关键。
专利CN103938277B和CN110628889A公开了两种含随机碱基(即UMI)的测序文库的构建方法，但其读取随机碱基的方式是通过索引引物(index primer)进行的，占据了区分不同样品的索引序列，给商业化测序带来限制。
发明内容
本发明第一目的是提供了一种含独特分子标签序列(UMI)的转座酶复合体，该转座酶复合体具有以下特点：
1)UMI不占用样品索引序列，便于使用商业化测序服务；
2)读取UMI的方式不同(通过read1primer或read2 primer读取)，从而允许更长的UMI。
为实现上述目的，本发明所设计一种含独特分子标签序列的转座酶复合体，所述转座酶复合体为由寡聚核苷酸接头和转座酶组装形成复合体；其中，
寡聚核苷酸接头的第一链内含有一个UMI标签序列，
或，含有多个依次连接的UMI标签序列组成UMI复合标签序列且UMI标签序列为一种或多种；
或，含有多个由间隔序列分隔的UMI标签序列组成UMI复合标签序列且UMI标签序列为一种或多种；其中，间隔序列由长度为0～20bp的可替换的非随机的序列组成，GC含量为25％-75％。
进一步地，所述寡聚核苷酸接头的第一链从5′端到3′端依次包括通用接头序列、间隔序列、UMI标签序列/UMI复合标签序列、间隔序列、条码序列和ME序列；间隔序列为0-20bp的序列；
所述寡聚核苷酸接头的第二链为包含ME序列的互补序列(第二链的长度不超过第一链的长度)。
再进一步地，所述通用接头序列为10-50bp且序列(DNA序列)中，GC含量为25％-75％(其作用为与之后建库所用PCR引物结合)；
UMI标签序列为2-100bp的序列(用于区分转座酶插入重复事件和PCR重复)，条码序列为0-100bp的序列(用于区分不同样品或保持碱基组成平衡)，ME序列用于与转座酶结合的序列(不同的酶选用不同的ME序列)，间隔序列为0-20bp的序列(用于间隔通用接头、多个UMI及条码序列，以及用于区分不同样品或保持碱基组成平衡)。
再进一步地，所述通用接头序列为含有“TTCCGATCT”的序列。
再进一步地，所述寡聚核苷酸接头的第一链的表达式为：
[]TTCCGATCT{}<>AGATGTGTATAAGAGACA，
其中，[]TTCCGATCT为通用接头序列，其长度为10～50个核苷酸；
{}为UMI标签序列或UMI复合标签序列，
若{}为UMI标签序列，其长度为2～50个核苷酸；
或，若{}为UMI复合标签序列，其长度为4～50个核苷酸；
<>为条码序列，其长度为0～50个核苷酸；
AGATGTGTATAAGAGACA为ME序列；
再进一步地，所述寡聚核苷酸接头的第一链选自下列序列中任意一条：
其中，CGACGCTCTTCCGATCT为通用接头序列，AGATGTGTATAAGAGACAG为ME序列；
含有“N”的序列均为{}<>的组合。
作为最优选方案，所述寡聚核苷酸接头的第一链选自下列序列中任意一条：
ME-A6 cgacgctcttccgatctnnnnnnagatgtgtataagagacag
ME-A13 cgacgctcttccgatctnnnnnnnnnnnnnagatgtgtataagagacag
ME-A20 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A25 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnnctgctagatgtgtataagagacag
再进一步地，所述转座酶(商品化或自主表达的酶)选自Tn5转座酶、Tn5变体、Mu、Mu E392Q、Vibhar、RAG和Tn552。
本发明第二目的是提供了上述转座酶复合体在RNA测序文库构建、DNA测序文库构建和染色质可及性测序(ATAC-seq)文库的构建中的应用。
本发明提供了一种用于构建测序文库的试剂盒，它包括复合体组合物、转座反应缓冲液、PCR引物、PCR酶和PCR扩增缓冲液；其中，所述复合体组合物包含一种上述的转座酶复合体/多种上述的转座酶复合体。
利用上述试剂盒构建含UMI的染色质可及性测序(UMI-ATAC-seq)文库的方法包括以下步骤：
1)材料准备和提取细胞核；
2)利用试剂盒转座、纯化和PCR，得到PCR产物；
3)PCR产物测序：
将PCR产物磁珠纯化后，采用Illumina系列平台进行PE150测序；其中，PCR产物的序列结构如下：
P5—index 5—Universal connector A—UMI—barcode—ME—DNA insert—ME—Universal connector B—index 7—P7；
公式中，P5序列和P7序列是Illumina标准文库要求的序列，与测序芯片(flowcell)上固定的序列互补或一致；
index 5和index 7序列为6-8bp可替换的DNA序列(用于多个文库混合测序时区分不同的文库，可以对多个样品建库后混池测序)；
测序时根据universal connector A和universal connector B来设计对应的read1和read2的测序引物，index 5和index 7的测序引物同Illumina测序仪所用；
3)对测序结果进行生物信息分析，得到数据质控结果。
本发明的有益效果：
1.利用本发明所述的含有UMI标签的转座酶复合体构建的二代测序文库，能够有效区分转座酶插入重复事件和PCR导致的重复读段，提高了测序数据的利用率，提高测序读段的利用率、序列变异及其频率鉴定的可靠性、基因表达和染色质可及性定量的准确性，以及鉴定DNA结合蛋白的足迹(Footprint)的敏感性。
2.利用本发明构建的含有UMI标签的转座酶复合体，可以构建兼容主流测序仪的文库，便于商业化推广应用。
3.本发明设计的UMI标签长度为2-100bp。通过不同长度UMI组合，有利于文库碱基平衡性，可以有效提高测序质量和兼容性。
4.本发明设计的兼容标准TruSeq文库的接头，可将所构建的文库与标准IlluminaTruSeq或Nextera文库一起混合上机测序，不需要使用自定义测序引物，可以降低成本、提高效率。
附图说明
图1为含独特分子标签序列的转座酶复合体及其构建文库和应用的流程图；
图2为含有UMI和条码的寡核苷酸接头的第一链的结构示意图；
图3为含有UMI接头的Illumina测序文库的结构示意图；
图4为实施例用含UMI的转座酶复合体进行的ATAC-seq实验的流程图；
图中，图4a为组装含UMI的转座酶复合体的示意图；
图4b为ATAC-seq实验的流程图，流程步骤依次为获取细胞核、与转座酶复合体孵育、纯化DNA、PCR构建文库和上机测序；
图5为样本C019中具有相同比对位点但UMI不同的被挽救读段的分布图；通常的测序文库中，比对位点相同的读段通常会被当成PCR重复而被去除，而如果使用含UMI的转座酶复合体构建文库，根据其所含UMI不同，部分读段可以被“挽救”；
图中，框图说明了UMI可以在高染色质可及性的基因组区段挽救更多的读段。将基因组划分为1kb的区块，并计算每个区块中被挽救的读段的比例，然后将所有的区块按25％、50％和75％的量级进行分组，获得此图。本分析中只有比对质量大于三十的读段被使用；
图6为不同方法去除PCR重复的定量结果相关性图；
图中，通过散点图对比了无重复数据删除(ND)、基于坐标的重复数据删除(CD)和基于UMI的重复数据删除(UD)数据集的量化结果；每个点是250bp区块(对数基数2)中Tn5插入的数量；
图7为UMI的使用提高了转录因子足迹识别的灵敏度图；
图中，图7a为显示了数据经CD和UD方式处理后，由pyDNase确定的足迹(P值<10-30)的维恩图，其中，共同的足迹被定义为至少有一个碱基重叠的足迹；
图7b数据经CD和UD方式处理后，pyDNase计算获得的足迹分数的关系图；其中，足迹分数的计算方式为log(-log10(P值))，点的颜色与a中的维恩图相同，小图显示的是仅CD方式处理后鉴定的足迹；
图7c为数据经CD(蓝色)、UD(红色)或ND(灰色)方式处理后，一些鉴定的足迹周围的Tn5插入频率；阴影区域对应的基因组坐标显示在每个图的上方；
图7d为足迹深度(FPD)计算的示意图，其中，"Avg"代表平均值函数；
图7e为数据经CD和UD方式处理后，计算的FPD值的差异的分布图，所涉及的足迹是指数据经CD方式处理后鉴定的足迹。
具体实施方式
下面结合具体实施例对本发明作进一步的详细描述，以便本领域技术人员理解。
实施例1
含独特分子标签序列的转座酶复合体构建(组装)：
1.设计寡聚核苷酸接头
ME-A6 cgacgctcttccgatctnnnnnnagatgtgtataagagacag
MErev [phos]ctgtctcttatacacatct
注：ME-A6为寡聚核苷酸接头的第一链，其中CGACGCTCTTCCGATCT为通用接头序列，NNNNNN为UMI序列，AGATGTGTATAAGAGACAG为ME序列，与形成完整接头所需的第二链MErev反向互补。
2.转座酶复合体组装
a.将MErev和ME-A6按照1:1混合，一次混合20μl；
b.将混合的引物经过高温变性(95℃3min)后，逐渐降温退火至12℃，获得5′粘性末端的接头MEDS-A6，即为组装所用接头。
c.按照如下比例，将上述合成的接头与Tn5酶(Tn5酶可通过商业化途径购买获得)混合：30μL Tn5酶，4.2μL MEDS-A6(100μM)；25℃，孵育60min，即得到转座酶复合体Tn5-UMI-A6，-20℃保存。
按上述方法将寡聚核苷酸接头的第一链替换为下述表中的寡聚核苷酸，可组装成各种转座酶复合体：
表1
实施例2
试剂盒包括实施例2获得的复合体组合物，缓冲液，PCR引物，PCR酶和PCR扩增缓冲液；
转座反应缓冲液体系为50mM Tris-HCl、25mM MgCl2和50％v/v DMF或其他适合Tn5转座反应的缓冲液。
PCR引物为：
其中下划线部分为可替换的条码序列。
PCR酶为高保真酶，如High-Fidelity 2X PCR Master Mix。
上述试剂盒中的复合体组合物制备方法如下：
1.设计寡聚核苷酸接头集合
合成融合UMI的Tn5组装所需接头序列，见表2。
表2.Tn5接头序列引物
注：ME-A25中的小写字符碱基为调节碱基平衡的条码序列。ME-B为文献(Jason DBuenrostro等，Nature Methods，2013，10:1213)所述序列，与MErev形成另一种寡聚核苷酸接头，便于PCR构建测序文库。
2.复合体组合物的组装(图4a)
a.将上述引物按照如下组合1:1混合，MEDS-A6(MErev+ME-A6)、MEDS-A13(MErev+ME-A13)、MEDS-A20(MErev+ME-A20)、MEDS-A25(MErev+ME-A25)、MEDS-B(MErev+ME-B)，每组引物对应一个离心管，分别混合20μl。
b.将混合的引物经过高温变性(95℃3min)后，逐渐降温退火至12℃，获得组装所用5′粘性末端的接头；
c.将上述四种接头MEDS-A6、MEDS-A13、MEDS-A20、MEDS-A25等摩尔混合为MEDS-A；
d.按照如下比例，将合成的接头与Tn5酶混合：30μL Tn5，2.1μL MEDS-A(100μM)，2.1μL MEDS-B(100μM)。25℃，孵育60min，得到复合体混合物，-20℃保存。
实施例3
以水稻叶片为材料，利用实施例2的试剂盒构建含有UMI接头的ATAC-seq文库(以下称为UMI-ATAC-seq)：
步骤S1，材料准备和文库构建(图4b)
a.水稻叶片细胞核制备：取新鲜水稻叶片，加入500μL细胞裂解液，用刀片快速剁碎，过滤后加入0.1％DAPI染色。用流式细胞仪分选100,000个细胞核至1.5mL离心管中，离心收集细胞核，冰上备用。
b.转座反应和纯化：在冰上按照如下比例配置反应体系：2μL复合体组合物，8μL5x tagmentation buffer，30μL ddH2O。吸打混匀后，37℃，孵育30min。反应完成后用TaKaRa MiniBEST DNAFragment Purification Kit(No.9761)纯化。
c.PCR扩增：采用High-Fidelity 2X PCR Master Mix(M0541)进行，反应体系如下：纯化DNA 20μL，2X PCR Master Mix25μL，PCR Primer 1(10μM)2.5μL，Barcoded PCR Primer 2(10μM)2.5μL。PCR反应程序为：72℃3min，98℃30sec；98℃10sec，63℃30sec，72℃1min，11cycles；72℃5min。
PCR所用引物序列见表3。
表3.PCR扩增引物
注：[i7]为标签序列，长度为8bp，本例中使用的标签序列共有7个，为TGCCTCTT、TCCTCTAC、ATCACGAC、ACAGTGGT、CAGATCCA、ACAAACGG、ACCCAGCA。用于针对多样品混合测序时区分不同样品的测序读段。
d.PCR产物测序：将PCR产物磁珠纯化后，采用Illumina HiSeq X-ten平台进行PE150测序；其中，PCR产物的序列结构如下：
P5—index 5—Universal connector A—UMI—barcode—ME—DNA insert—ME—Universal connector B—index 7—P7(图3)。
上述结构中，P5序列、P7序列是Illumina标准文库的要求，与测序芯片(flowcell)上固定的序列互补或一致。
index 5和index 7序列为6-8bp DNA固定序列，用于多个文库混合测序时区分不同的文库，可以对多个样品建库后混池测序。
测序时根据universal connector A和universal connector B来设计对应的read1和read2的测序引物，index 5和index 7的测序引物同Illumina测序仪所用。
本实施例中，read1测序引物为
“ACACTCTTTCCCTACACGACGCTCTTCCGATCT”，即Illumina测序试剂盒中包含的TruSeq标准测序引物。
步骤S2，对测序结果进行生物信息分析，得到数据质控结果。用fastqc对测序数据进行质量统计；去除测序接头后，利用bwa将测序数据比对至水稻参考基因组，统计比对率；samtools去除重复及低质量比对的序列并统计重复率。
步骤S3，评估UMI的应用对数据利用率的影响
表4 UMI-ATAC-seq测序数据质控结果
注：读段数量：左右两端均比对上的且在合理距离范围内的读段数目。UMI重复率：按照比对至同一位点，片段长度一致且UMI一致的片段仅保留一条的方式去除重复读段的比例。坐标重复率：按照比对至同一位点且片段长度一致的片段仅保留一条的方式去除重复读段的比例。UMI拯救率：考虑UMI后“拯救”的读段占去除的所有重复(不考虑UMI)读段的比例。FRiP：计算落在读段富集区域中的读段比例。TSS score：计算TSS周围(上下游3000bp)读段最高富集倍数。
通过实施例可以获得7个水稻样品ATAC-seq信息，在混库测序后，成功拆分得到7个品种的测序读段(如表3)。将UMI标签进行分类，发现通过重复插入事件造成的重复序列占比情况如表4。
结果表明，UMI-ATAC-seq可以挽救标准ATAC-seq中约20％被误认为PCR重复的读段(或约占总比对到基因组的读段的6％)。进一步将基因组分成250个bp的bins，并计算每个bin中的被挽救的读段数。发现染色质上Tn5越容易接近的区域，通过UMI方法可以挽救的读段数就越多。在样本C019中，81.8％被挽救的读段数来自水稻基因组中最易接近的前2.0％的bins。并且，被挽救的读段数在易受拷贝数变异或基因组组装错误影响的区域中富集，例如rDNA区域、转座子、逆转录转座子和端粒区(图5)。
为了评估加入UMI对定量的影响，分别通过基于相同基因组坐标进行的PCR重复数据消除(CD)、基于UMI的重复数据消除(UD)和不进行重复数据消除(ND)来量化每个250bpbin的Tn5插入，并相互比较。结果表明，UD和ND的相关性略高于UD和CD(R2分别为0.995和0.961)；并得到了与其他六个样品相似的结果(图6a、b)。随着染色质区域可及性的增加，CD的量化更偏离UD(图6c)。这些结果表明，在ATAC-seq的定量分析中，仅基于比对坐标的PCR重复的去除可能比不去除重复结果更糟糕。
使用pyDNase4对通过CD和UD鉴定到的转录因子足迹进行进一步的验证。大多数足迹都是在高度易接近的区域被识别出来的，因此在去除PCR重复的方法中受到的影响更大。UD识别的足迹数量明显大于CD(图7a)。虽然只有约6％的读段被挽救，但它们却促成了50％或更多的足迹被发现。比较重叠足迹的得分，发现UD计算的足迹得分几乎总是高于CD(图7b)。
对于仅由CD或UD识别的足迹，使用相同的足迹区间计算UD和CD的足迹得分。结果表明，在大多数情况下，UD足迹得分也高于CD。使用IGV5(整合基因组学观察者)检查已识别足迹周围Tn5插入的频率，发现足迹样分布在UD中更为明显(图7c)。对于仅限于UD的足迹，许多独立的Tn5插入作为PCR重复被移除，导致CD无法识别。此外，足迹深度(FPD)值越大，可以解释为足迹区的Tn5插入量较其侧翼区减少(图7d)。通过计算FPD值，结果显示UD的FPD值普遍高于CD，说明UD模式下足迹模式更为明显(图7e)。因此，UMI-ATAC-seq可以在染色质高度可及的区域挽救了大量独立的Tn5插入事件，从而提高转录因子足迹绘制的灵敏度和准确性。
其它未详细说明的部分均为现有技术。尽管上述实施例对本发明做出了详尽的描述，但它仅仅是本发明一部分实施例，而不是全部实施例，人们还可以根据本实施例在不经创造性前提下获得其他实施例，这些实施例都属于本发明保护范围。
序列表
<110> 华中农业大学
<120> 含独特分子标签序列的转座酶复合体及其应用
<160> 25
<170> SIPOSequenceListing 1.0
<210> 1
<211> 38
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 1
cgacgctctt ccgatctnna gatgtgtata agagacag 38
<210> 2
<211> 39
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 2
cgacgctctt ccgatctnnn agatgtgtat aagagacag 39
<210> 3
<211> 40
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 3
cgacgctctt ccgatctnnn nagatgtgta taagagacag 40
<210> 4
<211> 41
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 4
cgacgctctt ccgatctnnn nnagatgtgt ataagagaca g 41
<210> 5
<211> 42
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 5
cgacgctctt ccgatctnnn nnnagatgtg tataagagac ag 42
<210> 6
<211> 43
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 6
cgacgctctt ccgatctnnn nnnnagatgt gtataagaga cag 43
<210> 7
<211> 44
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 7
cgacgctctt ccgatctnnn nnnnnagatg tgtataagag acag 44
<210> 8
<211> 45
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 8
cgacgctctt ccgatctnnn nnnnnnagat gtgtataaga gacag 45
<210> 9
<211> 46
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 9
cgacgctctt ccgatctnnn nnnnnnnaga tgtgtataag agacag 46
<210> 10
<211> 47
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 10
cgacgctctt ccgatctnnn nnnnnnnnag atgtgtataa gagacag 47
<210> 11
<211> 48
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 11
cgacgctctt ccgatctnnn nnnnnnnnna gatgtgtata agagacag 48
<210> 12
<211> 49
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 12
cgacgctctt ccgatctnnn nnnnnnnnnn agatgtgtat aagagacag 49
<210> 13
<211> 50
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 13
cgacgctctt ccgatctnnn nnnnnnnnnn nagatgtgta taagagacag 50
<210> 14
<211> 51
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 14
cgacgctctt ccgatctnnn nnnnnnnnnn nnagatgtgt ataagagaca g 51
<210> 15
<211> 52
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 15
cgacgctctt ccgatctnnn nnnnnnnnnn nnnagatgtg tataagagac ag 52
<210> 16
<211> 53
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 16
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnagatgt gtataagaga cag 53
<210> 17
<211> 54
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 17
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnnagatg tgtataagag acag 54
<210> 18
<211> 55
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 18
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnnnagat gtgtataaga gacag 55
<210> 19
<211> 56
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 19
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnnnnaga tgtgtataag agacag 56
<210> 20
<211> 57
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 20
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnnnnnag atgtgtataa gagacag 57
<210> 21
<211> 58
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 21
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnnnnnna gatgtgtata agagacag 58
<210> 22
<211> 59
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 22
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnnnnnnn agatgtgtat aagagacag 59
<210> 23
<211> 60
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 23
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnnnnnnn nagatgtgta taagagacag 60
<210> 24
<211> 61
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 24
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnnnnnnn nnagatgtgt ataagagaca 60
g 61
<210> 25
<211> 61
<212> DNA
<213> 人工合成序列(Synthetic sequence)
<400> 25
cgacgctctt ccgatctnnn nnnnnnnnnn nnnnnnnctg ctagatgtgt ataagagaca 60
g 61

Claims (5)

1.一种含独特分子标签序列的转座酶复合体，其特征在于：所述转座酶复合体为由寡聚核苷酸接头和转座酶组装形成的复合体；转座酶为Tn5转座酶；且所述寡聚核苷酸接头的第一链从5′端到3′端依次包括通用接头序列、间隔序列、UMI标签序列或者UMI复合标签序列、间隔序列、条码序列和ME序列，间隔序列由长度为0～20bp的可替换的非随机的序列组成，GC含量为25％-75％；ME序列为用于与转座酶结合的序列；
所述寡聚核苷酸接头的第一链的表达式为：
[]cgacgctcttccgatct{}<>AGATGTGTATAAGAGACA，
其中，cgacgctcttccgatct为通用接头序列的末端，通用接头序列[]cgacgctcttccgatct的长度为10～50个核苷酸；
{}为UMI标签序列或UMI复合标签序列，UMI复合标签由多个依次连接的UMI标签序列组成且UMI标签序列为一种或多种；
或者由多个由间隔序列分隔的UMI标签序列组成且UMI标签序列为一种或多种；
UMI标签序列长度为2～50个核苷酸；
UMI复合标签序列长度为4～50个核苷酸；
<>为条码序列，其长度为0～50个核苷酸，用于调节碱基平衡；
AGATGTGTATAAGAGACA为ME序列；
所述寡聚核苷酸接头的第二链为包含ME序列的互补序列。
2.根据权利要求1所述含独特分子标签序列的转座酶复合体，其特征在于：所述寡聚核苷酸接头的第一链选自下列序列中任意一条：
ME-A2 cgacgctcttccgatctnnagatgtgtataagagacag
ME-A3 cgacgctcttccgatctnnnagatgtgtataagagacag
ME-A4 cgacgctcttccgatctnnnnagatgtgtataagagacag
ME-A5 cgacgctcttccgatctnnnnnagatgtgtataagagacag
ME-A6 cgacgctcttccgatctnnnnnnagatgtgtataagagacag
ME-A7 cgacgctcttccgatctnnnnnnnagatgtgtataagagacag
ME-A8 cgacgctcttccgatctnnnnnnnnagatgtgtataagagacag
ME-A9 cgacgctcttccgatctnnnnnnnnnagatgtgtataagagacag
ME-A10 cgacgctcttccgatctnnnnnnnnnnagatgtgtataagagacag
ME-A11 cgacgctcttccgatctnnnnnnnnnnnagatgtgtataagagacag
ME-A12 cgacgctcttccgatctnnnnnnnnnnnnagatgtgtataagagacag
ME-A13 cgacgctcttccgatctnnnnnnnnnnnnnagatgtgtataagagacag
ME-A14 cgacgctcttccgatctnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A15 cgacgctcttccgatctnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A16 cgacgctcttccgatctnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A17 cgacgctcttccgatctnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A18 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A19 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A20 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A21 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A22 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A23 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A24 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A25-2 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnnnnnnnagatgtgtataagagacag
ME-A25 cgacgctcttccgatctnnnnnnnnnnnnnnnnnnnnctgctagatgtgtataagagacag
其中，CGACGCTCTTCCGATCT为通用接头序列，AGATGTGTATAAGAGACAG为ME序列；
“n”是所述权利要求1中{}<>的组合。
3.一种如权利要求1所述的转座酶复合体在RNA测序文库构建、DNA测序文库构建和染色质可及性测序文库的构建中的应用。
4.一种用于构建测序文库的试剂盒，其特征在于，它包括复合体组合物、转座反应缓冲液、PCR引物、PCR酶和PCR扩增缓冲液；其中，所述复合体组合物包括一种或多种权利要求1所述的转座酶复合体。
5.利用权利要求4所述试剂盒构建含UMI的染色质可及性测序文库的方法，其特征在于：包括以下步骤：
1)材料准备和提取细胞核；
2)利用试剂盒转座、纯化和PCR，得到PCR产物；
3)PCR产物测序：
将PCR产物磁珠纯化后，采用Illumina系列平台进行测序；其中，PCR产物的序列结构如下：P5—index 5—通用接头A—UMI—条码序列—ME—DNA insert—ME—通用接头B—index 7—P7；其中，P5序列和P7序列是Illumina标准文库要求的序列，与测序芯片上固定的序列互补或一致；
index 5和index 7序列为6-8bp可替换的DNA序列；
测序时根据通用接头A和通用接头B来设计对应的read1和read2的测序引物，index 5和index 7的测序引物为Illumina测序仪所用的测序引物；
4)对测序结果进行生物信息分析，得到数据质控结果。
CN202011129375.XA 2020-10-21 2020-10-21 含独特分子标签序列的转座酶复合体及其应用 Active CN112251422B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011129375.XA CN112251422B (zh) 2020-10-21 2020-10-21 含独特分子标签序列的转座酶复合体及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011129375.XA CN112251422B (zh) 2020-10-21 2020-10-21 含独特分子标签序列的转座酶复合体及其应用

Publications (2)

Publication Number Publication Date
CN112251422A CN112251422A (zh) 2021-01-22
CN112251422B true CN112251422B (zh) 2024-04-19

Family

ID=74263695

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011129375.XA Active CN112251422B (zh) 2020-10-21 2020-10-21 含独特分子标签序列的转座酶复合体及其应用

Country Status (1)

Country Link
CN (1) CN112251422B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113604545B (zh) * 2021-08-09 2022-04-29 浙江大学 一种超高通量单细胞染色质转座酶可及性测序方法
CN113517026B (zh) * 2021-06-16 2022-08-19 苏州拉索生物芯片科技有限公司 应用于生物制品的标签序列的生成方法、系统、智能终端及计算机可读存储介质
CN114645074A (zh) * 2022-03-07 2022-06-21 华南农业大学 一种单花蕾微量rna建库方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103938277A (zh) * 2014-04-18 2014-07-23 中国科学院北京基因组研究所 以痕量dna为基础的二代测序文库构建方法
CN107075509A (zh) * 2014-05-23 2017-08-18 数字基因公司 通过数字化转座子的单倍体组测定
CN109154013A (zh) * 2016-03-31 2019-01-04 安捷伦科技有限公司 转座酶和y衔接子用于片段化和标签化dna的用途
CN109526228A (zh) * 2017-05-26 2019-03-26 10X基因组学有限公司 转座酶可接近性染色质的单细胞分析
WO2019076768A1 (en) * 2017-10-16 2019-04-25 Tervisetehnoloogiate Arenduskeskus As METHOD AND KIT FOR PREPARING DNA BANK
CN110628889A (zh) * 2019-08-12 2019-12-31 远辰生物科技(苏州)有限公司 Illumina二代测序平台文库构建引入分子标签的方法、接头序列及应用
CN110678547A (zh) * 2017-05-23 2020-01-10 生物辐射实验室股份有限公司 分子条码化
WO2020014303A1 (en) * 2018-07-11 2020-01-16 X Gen Us Co. Transposome enabled dna/rna-sequencing (ted rna-seq)
CN110997937A (zh) * 2017-09-15 2020-04-10 伊鲁米那股份有限公司 具有可变长度非随机独特分子标识符的通用短衔接子
CN111201329A (zh) * 2018-05-17 2020-05-26 伊鲁米纳公司 具有减少的扩增偏倚的高通量单细胞测序
CN112824534A (zh) * 2019-11-20 2021-05-21 武汉华大医学检验所有限公司 对核酸的目标区域进行扩增的方法、建库和测序方法及试剂盒

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018195224A1 (en) * 2017-04-18 2018-10-25 Fred Hutchinson Cancer Research Center Barcoded transposases to increase efficiency of high-accuracy genetic sequencing

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103938277A (zh) * 2014-04-18 2014-07-23 中国科学院北京基因组研究所 以痕量dna为基础的二代测序文库构建方法
CN107075509A (zh) * 2014-05-23 2017-08-18 数字基因公司 通过数字化转座子的单倍体组测定
CN109154013A (zh) * 2016-03-31 2019-01-04 安捷伦科技有限公司 转座酶和y衔接子用于片段化和标签化dna的用途
CN110678547A (zh) * 2017-05-23 2020-01-10 生物辐射实验室股份有限公司 分子条码化
CN109526228A (zh) * 2017-05-26 2019-03-26 10X基因组学有限公司 转座酶可接近性染色质的单细胞分析
CN110997937A (zh) * 2017-09-15 2020-04-10 伊鲁米那股份有限公司 具有可变长度非随机独特分子标识符的通用短衔接子
WO2019076768A1 (en) * 2017-10-16 2019-04-25 Tervisetehnoloogiate Arenduskeskus As METHOD AND KIT FOR PREPARING DNA BANK
CN111201329A (zh) * 2018-05-17 2020-05-26 伊鲁米纳公司 具有减少的扩增偏倚的高通量单细胞测序
WO2020014303A1 (en) * 2018-07-11 2020-01-16 X Gen Us Co. Transposome enabled dna/rna-sequencing (ted rna-seq)
CN110628889A (zh) * 2019-08-12 2019-12-31 远辰生物科技(苏州)有限公司 Illumina二代测序平台文库构建引入分子标签的方法、接头序列及应用
CN112824534A (zh) * 2019-11-20 2021-05-21 武汉华大医学检验所有限公司 对核酸的目标区域进行扩增的方法、建库和测序方法及试剂盒

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Georgia Giannoukos等.UDiTaS™, a genome editing detection method for indels and genome rearrangements.《BMC Genomics》.2018, *
UDiTaS™, a genome editing detection method for indels and genome rearrangements;Georgia Giannoukos等;《BMC Genomics》;20180321;图1和表2 *
我国科学家开发出转录组快速建库新方法 有望助力新型冠状病毒测序;新华网;《微循环学》;20200331;第55页下半部分 *

Also Published As

Publication number Publication date
CN112251422A (zh) 2021-01-22

Similar Documents

Publication Publication Date Title
CN112251422B (zh) 含独特分子标签序列的转座酶复合体及其应用
CN110997937B (zh) 具有可变长度非随机独特分子标识符的通用短衔接子
CN105506125B (zh) 一种dna的测序方法及一种二代测序文库
CN105696088B (zh) 一种双酶切简化基因组二代测序文库构建方法及配套试剂盒
EP1546345B1 (en) Genome partitioning
CN110546272B (zh) 将衔接子附接至样品核酸的方法
CN106048009B (zh) 一种用于超低频基因突变检测的标签接头及其应用
US20220389416A1 (en) COMPOSITIONS AND METHODS FOR CONSTRUCTING STRAND SPECIFIC cDNA LIBRARIES
CN110628891B (zh) 一种对胚胎进行基因异常筛查的方法
JP2009536817A (ja) 制限断片に基づく分子マーカーのハイスループットな検出
CN113463202B (zh) 一种新的rna高通量测序的方法、引物组和试剂盒及其应用
CN111321202A (zh) 基因融合变异文库构建方法、检测方法、装置、设备及存储介质
EP2531610A1 (en) Complexitiy reduction method
CN112359093B (zh) 血液中游离miRNA文库制备和表达定量的方法及试剂盒
CN108192893B (zh) 基于转录组测序开发艾纳香ssr引物的方法
JP2020529833A (ja) モジュラー核酸アダプター
WO2024093961A1 (zh) 一种简化基因组测序的方法与相关应用
CN115747208B (zh) 一种dna/rna混合物的处理方法
CN111304299A (zh) 一种用于检测常染色体拷贝数变异的引物组合、试剂盒和方法
CN109355289A (zh) 用于从预文库中富集靶标序列的试剂盒、制备方法和应用
US20190218606A1 (en) Methods of reducing errors in deep sequencing
CN111763668B (zh) 测序引物组及基于pcr的全基因组测序方法
US20200190567A1 (en) Method For Detecting Activity Change Of Transposon In Plant Before And After Stress Treatment
CN114854825A (zh) 一种适用于dnbseq技术的简化基因组测序的建库接头及方法
CN107447050B (zh) 一种利用转录组测序全面快速检测大蒜rna病毒的方法

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