CN107250381B - Dna集合的归一化迭代条形码和测序 - Google Patents

Dna集合的归一化迭代条形码和测序 Download PDF

Info

Publication number
CN107250381B
CN107250381B CN201680010540.5A CN201680010540A CN107250381B CN 107250381 B CN107250381 B CN 107250381B CN 201680010540 A CN201680010540 A CN 201680010540A CN 107250381 B CN107250381 B CN 107250381B
Authority
CN
China
Prior art keywords
sample
nucleic acid
pool
samples
labeled
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
CN201680010540.5A
Other languages
English (en)
Other versions
CN107250381A (zh
Inventor
约瑟夫·C·梅勒
杰克·T·莱昂纳德
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.)
Seqwell Inc
Original Assignee
Seqwell Inc
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 Seqwell Inc filed Critical Seqwell Inc
Publication of CN107250381A publication Critical patent/CN107250381A/zh
Application granted granted Critical
Publication of CN107250381B publication Critical patent/CN107250381B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12NMICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C12N15/09Recombinant DNA-technology
    • C12N15/10Processes for the isolation, preparation or purification of DNA or RNA
    • C12N15/1034Isolating an individual clone by screening libraries
    • C12N15/1065Preparation or screening of tagged libraries, e.g. tagged microorganisms by STM-mutagenesis, tagged polynucleotides, gene tags
    • 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
    • C12Q1/6874Methods for sequencing involving nucleic acid arrays, e.g. sequencing by hybridisation
    • 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
    • C40B40/08Libraries containing RNA or DNA which encodes proteins, e.g. gene libraries

Abstract

本发明尤其提供用于由多个原始含核酸的样品制备线性、非选择性扩增的DNA片段的统一库的组合物和方法,其中来自所述多个原始样品中的每一个的所述片段的比例表述为归一化的并且所述库是以高度并行化、基于池方式产生。本发明尤其用于制备如下库:其中特定信息经编码而使由所述库的高通量测序衍生的较短测序读序被分析或组合成较长级序列,所述较长级序列在可能的极大型样品集合中完全可追踪到原始含核酸的样品。本发明的组合物涵盖本文所述的多种构筑体,其可与一种或多种用于本发明方法中的另外试剂和使用说明书一起以不同方式封装。

Description

DNA集合的归一化迭代条形码和测序
技术领域
本发明涉及DNA测序的技术领域。更具体而言,本文所述的组合物和方法用于生成来自多个来源的DNA库且用于下一代测序。
背景技术
在多种不同研究和临床努力中,期望并行获得多个不同样品中所含的核酸序列信息。
通常被称为“Sanger测序”的单独样品或单重DNA测序的普遍存在的实施方案是使用DNA聚合酶以DNA模板依赖性方式将荧光双脱氧核苷酸终止子并入一系列嵌套延伸产物中,之后,在耦接到荧光计的毛细管电泳仪器上分离延伸产物来进行,所述荧光计会输出对应于连续核苷酸并入事件的信号。Sanger测序当前被视为单独DNA读序质量的黄金标准;其尤其在每次从少量DNA模板(<100个)获得序列结果时,并且尤其在关注区域的长度处于约数百个与750个碱基对之间时,具有精确性和可靠性。
大规模复用测序(有时被称为下一代测序(NGS)、大规模平行测序(MPS)或第二代测序)涵盖数种技术平台,其使得一个样品中的数千到数百万个异质DNA分子可以同时得到准备和分析。因此,大型混合物(通常被称为“库”)中的多种核苷酸序列可在单一分析中得到检索。
除了在单一分析中所获得的序列数量方面不同以外,单重(即,Sanger)和复用测序还在其如何准许从集合追踪分离样品方面有所不同。在单重测序中,经特定确认的序列的来源或位置一般为明确的,因为模板DNA的单一分离样品为分析的输入并且可通过对输出序列数据进行直接追踪来联系。举例来说,可以在电泳之前在其中装载有样品的微量滴定板或毛细管上记录样品位置。在复用/NGS测序中,通过使用输入条形码、标记或索引方案,集合中的特定模板DNA的位置或来源可通过将可鉴别的序列标记(即,不同DNA序列)并入由特定孔或位置衍生的库片段中来得到确认。然而,对于测序数千个分离或可追踪样品来说,用于独立样品的分子标记的常规输入-索引模式(有时被简单地称为“条形码”)通常过于昂贵。
Sanger测序仍为各种情况下,且尤其在需要(1)较长的高质量(低错误率)序列,(2)测序模板或DNA克隆的板或位置集合上的位置信息,和/或(3)项目成本、设计、样品数量或与常规NGS方法相关的其它材料约束条件方面的总体实用性的情况下用于进行DNA测序的默认方法。举例来说,Sanger测序在例如测序筛分(在工程改造肽期间的克隆鉴别的情形下进行)、定制的长扩增子组的全长测序或通过化学基因合成生成的全长产物的测序中仍然是选定方法。
NGS已取代Sanger测序,其中通常以数百万或数十亿个碱基读数测量的所获得的原始数据量具有高实用性。举例来说,现在可以在一天之内以足够深度对人类基因组的全部大约60亿碱基进行测序。需要精确计算大量不同核酸的许多其它应用(例如对人类细胞中所表达的mRNA进行测序和定量)也极大地受益于现在可以以便宜和简单的方式获得的大量原始序列数据。
尽管现代测序仪器的容量已有所增长,但使用现代仪器的测序容量来同时测序较大数目的不同样品仍然具有相对挑战性。举例来说,通过常规估算,测序一个二倍体人类基因组(60亿个核苷酸)所需的测序容量的量大致与测序一千个大肠杆菌基因组(五百万个核苷酸)所需的容量相当。然而,与由一千个大肠杆菌样品常规制备用于下一代测序的核酸样品库相关的材料和成本将超过单一人类基因组样品的材料和成本至少100倍。这是构成在同时测序大量样品方面的困难的所谓的“库制备瓶颈”的一个实例。
已知用于由核酸样品制备适用于NGS的库的多种方法。代表性方法包括例如美国专利第8,420,319号;第8,728,737号和第9,040,256号中所述的方法。这些方法的统一特征为随机片段由核酸样品(例如,基因组DNA样品)产生,并且库是通过添加不对称衔接子,产生可通过多种手段(例如PCR)扩增的适合片段而衍生。利用这些方法,核酸样品的量或浓度也必须通常被定量为指定精确程度以便实现库制备方法的效率、质量或产率。
通过常规方法制造的库的复用通常通过将一种或多种可鉴别的序列标记或条码并入连接到每一衍生库中的片段的衔接子中来实现。可随后通过测量多个库中的每一个的量或浓度并将其混合在一起以制得统一库而将许多此类库以期望比例合并。在一些情况下,这一方法可通过使用机器人技术或能够测量浓度并将单独的库适应性地汇集成一个期望混合库的其它自动化方法来得到自动化。
可使用其它流线型归一化策略来提高库制备方法的并行化的效率。举例来说,可在建库之前使输入核酸样品归一化,或可使用具有有限可逆结合能力的样品容器(例如,赛默飞世尔(ThermoFisher)
Figure BDA0001380396140000031
归一化盘)或呈现有限可逆结合能力的核酸结合组合物(参见,例如美国专利第8,536,322号)使单独的库归一化。然而,这些方法增加有消耗额外时间并增加NGS库制备的材料成本的不同步骤。此外,样品仍然被以平行方式而非以共同方式处理,其意指这一不足之处将随样品数量而成直线增长。
因此所属领域已知的现行方法通常缺乏形成适用于NGS而不需要对输入样品或所得库进行单独调整或并行化归一化的复用化、归一化库的能力。
因此本发明的一个目的在于提供可由多个不同样品形成线性、非选择性扩增的DNA分子的统一库的方法,其适用于NGS,同时维持由每一样品产生的DNA分子的归一化百分比。
本发明的另一个目标在于提供可用于以并行化和汇集方式由多个不同样品形成线性、非选择性扩增的DNA分子的统一库的方法,所述方法不需要明确归一化每一样品中所存在的DNA的量,也不需要测量、评定或明确归一化多个不同库来制造统一库。
本发明的另一个目标在于减少将另外与使用当前在所属领域中已知或实践的方法中的任一个由多个不同样品制备DNA分子的归一化、统一库相关的困难和时间。
发明内容
本发明尤其提供由多个原始的含核酸的样品制备线性、非选择性扩增的DNA片段的统一库的组合物和方法,其中来自多个原始样品中的每一个的片段的比例表述为归一化的(即,片段以大致同等方式展现)并且库是以高度并行化、基于池方式产生。本发明尤其用于制备如下库:其中特定信息经编码而使由所述库的高通量测序衍生的较短测序读序被分析或组合成较长级序列,所述较长级序列在可能的极大型样品集合中完全可追踪到原始含核酸的样品。本发明的组合物涵盖本文所述的多种构筑体,其可与一种或多种用于本发明方法中的另外试剂和使用说明书一起以不同方式封装。
通常,在本发明方法中,首先对多个核酸样品中的每一样品进行样品标记步骤,其中每一样品中的多个核分子经修饰而并入可鉴别的序列标记。用包括唯一、可鉴别的序列标记的构筑体处理(例如在本文所述的条件或分子生物学家已知的常规反应条件下接触)待分析的每一样品以使得一个样品可基于所述标记区别于另一个样品。在优选实施例中,样品标记步骤在将所并入标记的数量限制为远低于每一样品中所存在的可利用的修饰位点的数量的数量(例如样品的核酸分子中所存在的每约1000个核苷酸碱基(例如小于约1,000、900、800、700、600或500个核苷酸碱基)小于一个所并入标记)的反应条件下进行,并且优选地,对于全部多个样品都保持相同。
根据本发明方法,可进行样品标记步骤以使得对于全部多个经处理样品来说,每一样品中所并入的可鉴别的序列标记将大致相等,其与每一样品中所存在的核酸的量或浓度无关。这可以在使用本文所述的指定反应条件,并且更一般来说,其中用于并入可鉴别的序列标记的试剂或构筑体相对于一系列可能的样品核酸浓度保持在一定限度内的至少一些实施例中实现。
在样品标记步骤之后,将多个以可鉴别方式标记的样品随后汇集并且在池标记步骤中进一步加以处理,在所述池标记步骤中,对池中的多个经样品标记的核分子进行修饰以并入鉴别池的其它可鉴别的序列标记(有时在本文中被称为“池标记”)。
在优选实施例中,进行池标记步骤所处的反应条件使得池标记被平均添加到来自池中核酸分子上所存在的任何样品标记的一系列X个核苷酸中。在各个实施例中,X可选择为指定值,例如约100、500或1000个核苷酸。当如某些优选实施例中所规定进行池标记步骤时,所述步骤的一个非显而易见的特征在于池标记的併入与池中所存在的任何核酸分子上所存在的任何样品标记的位置无关,并且也与由先前经历样品标记的任何特定样品产生的核酸分子的相对浓度无关。因此可以说,池中所存在的任何样品标记的X个核苷酸中出现池标记的情况得到归一化。
在池标记步骤之后,池中所存在的核酸分子可在其它实施例中得到扩增。扩增可用杂交到对于全部样品标记、全部池标记或这两者交替常见的序列的区域或子区域的引物由例如PCR来辅助。以此方式,获得由多个原始含核酸的样品产生的线性、非选择性扩增的DNA片段的统一、归一化库。
因此,在第一方面中,本发明提供一种由多个样品获得DNA片段的迭代标记库的方法。所述方法可包括以下步骤:(a)使多个样品的第一和第二样品中的靶核酸分子分别与均一和限定量的各自包含可鉴别的样品序列标记的第一和第二构筑体接触,其中第一构筑体中的样品序列标记与第二构筑体中的样品序列标记不同,并且接触会进行一定持续时间并在如下条件下进行:其中第一构筑体中的样品序列标记被整合到第一样品中的核酸分子中,并且第二构筑体中的样品序列标记被整合到第二样品中的核酸分子中;(b)汇集第一和第二样品以产生经样品标记的核酸分子的池;(c)使经样品标记的核酸分子的池与包含池标记的第三构筑体接触,其中接触会进行一定持续时间并在如下条件下进行:其中池标记以按照核苷酸的规定频率被整合到池的核酸分子中,由此生成经样品标记的、经池标记的核酸分子的池;以及(d)用寡核苷酸扩增经样品标记的、经池标记的核酸分子,由此生成经扩增核酸分子,所述寡核苷酸包含(i)与经样品标记的、经池标记的核酸分子中的核酸序列互补的核酸序列,和(ii)适用于下一代测序的核酸序列。当如此处和我们本发明描述中的其它处所述进行所述方法时,所述方法会产生DNA库,其中与样品对应的每一库中的DNA片段的相对丰度得到归一化。在各个实施例中,池标记包括可以按照核苷酸的规定频率(例如,在平均约100-20,000(例如,100-10,000或100-1,000)个核苷酸之间产生池标记的频率)整合到经汇集的、经样品标记的核酸分子中的可鉴别的序列标记,其远离池中的经样品标记的核酸分子上的每一样品序列标记。
多个样品中的至少一个样品可为或可包括基因组DNA。另外,多个样品中的至少一个样品可为或可包括cDNA、合成性DNA或在病毒、细菌、酵母菌、真菌、原虫、昆虫、鱼、禽、哺乳动物或植物中天然存在的DNA。不同类型的DNA或来自多个来源的DNA可存在于任何给定样品中。如所提到,可并行处理多个样品(例如2到约9,600个样品)(例如可并行处理约100个96孔培养盘,其每一孔中都含有样品)。
在本发明方法中,第一和第二构筑体可包括衔接子,所述衔接子随后连接到靶核酸分子;或转座酶或其生物活性变异体,其随后将可鉴别的样品序列标记引入到靶核酸分子中。适合的转座酶可极具活性并且几乎不呈现或不呈现序列偏差。举例来说,用于本发明方法中的构筑体可包括Tn3、Tn5、Tn9、Tn10、γ-δ、或Mu转座酶、其组合、或其生物活性变异体。
在样品序列标记被整合到靶核酸分子中的条件下,样品序列标记的量可以比随后用于接触经样品标记的核酸分子的池的池标记的量低至少两倍。
转而参看第三构筑体,这一构筑体也可包括随后经连接(例如连接到经样品标记的核酸分子)的衔接子或转座酶(包括上文所列的转座酶中的任一种)或其生物活性变异体,其随后将池标记引入到经样品标记的核酸分子中。第一、第二或第三构筑体中的任一种也可独立地包括特异性结合捕获剂的部分(例如,特异性结合捕获剂中的蛋白质或核酸序列的核酸序列)。举例来说,特异性结合捕获剂的部分可为已被生物素化或洋地黄毒苷化(digoxigenylated)的部分。
在本发明的方法中,可以通过将经样品标记的和/或经池标记的核酸分子群暴露于有限量的捕获剂来使所述群归一化。
上文所述的步骤之后可对DNA片段的迭代标记库进行下一代测序。
在另一方面中,本发明提供一个试剂盒,其包含:(a)多个构筑体(例如,2到约9,600个构筑体),其包括可鉴别的样品序列标记,其中多个构筑体的第一构筑体中的可鉴别的样品序列标记不同于第二构筑体中的可鉴别的样品序列标记;和(b)使用说明书。说明书可呈任何形式,包括书写/印刷形式,通过储存于任何合适媒体上的音频或视频录制提供,或通过网站类指导的链接提供。试剂盒可进一步包括池标记试剂(例如如本文所述的构筑体),其用于将可鉴别的序列标记添加到经样品标记的核酸的池中的一种或多种核酸中,和/或用于扩增经样品标记的和/或经池标记的核酸分子的试剂。
附图说明
此处经简要描述的附图意图帮助说明本发明并且并非限制性的。
图1为可用于实施本发明的方案和步骤的代表性实施例的示意性图解。所图示的方法以样品标记开始并且以测序分析结束。
图2图示本发明的一个实施例,其中核酸样品经历样品标记步骤,在所述步骤中,样品的一种或多种核酸标记有可鉴别的序列标记。
图3图示样品标记步骤的归一化方面,其中使用相对固定量的标记试剂将可鉴别的序列标记添加到不同量的DNA中。随后将来自多个样品的经标记的产物合并(汇集)。
图4图示池标记步骤;处理经汇集的、经样品标记的核酸以获得序列标记,其鉴别经样品标记的核酸作为池的成员。
图5图示与不同的和可鉴别的池标记一起经历池标记步骤的经样品标记的核酸分子的多个不同池。
图6图示经并行测序以从库中的多个片段中的每一个获得短读序信息的片段的库。对于这些片段中的每一个来说,一组短读序经产生而包含一种或多种样品标记、一种或多种集合标记和每一分子上的一种或多种其它关注区。
图7为展示在对使用本发明方法生成的库进行测序之后所获得的覆盖区的实例的图式(参见实例1中所述的研究)。“多个样品”为由伯克霍尔德菌(Burkolderia bacteria)的分离物收集的基因组DNA的96孔培养盘中的96个样品。
图8显示如实例2中所示,使用本发明方法由多个千碱基质粒模板获得的组合的连续核苷酸序列的实例。
图9为展示本发明方法的归一化表现的实例的图式,所述效能通过对具有不同原始输入核酸浓度的多个样品进行测序来确定。通过本发明方法实现的归一化作用可与在此处所报告的作用相同(或大致相同)。
具体实施方式
呈现以下描述以使得所属领域的一般技术人员能够进行和使用本发明的各个方面和实例。特定材料、技术和应用的描述仅作为实例而提供。所属领域的一般技术人员将容易明白本文所述的实例的各种修改,并且在不脱离本发明的精神和范围的情况下,本文中所界定的一般原理可以应用于其它实例和应用。因此,本发明并非意图限于所述和所示的实例,而是被赋予与所附权利要求书一致的范围。
在一个方面中,本发明包含由多个原始含核酸样品制备线性、非选择性扩增的DNA片段的统一库的方法,其中来自多个原始样品中的每一个的线性片段的比例表述为归一化的。这些方法尤其用于制备DNA测序库,其中特定信息经编码而提供通过对待分析或待组合成较长级序列的库进行高通量短读序测序而衍生的测序读序,其在可能的大型样品集合中全部可源于原始分离样品。
在另一个方面中,本发明包含一种或多种试剂盒,其提供用于实现所揭示的本发明方法的试剂和说明书。
除非使用术语的上下文另外明确指明,否则本发明的传授内容利用含义以以下方式解释的术语。
在本文中术语“约”用于表明如下值:其包括用于测定值的装置或方法的内在误差变化;或表明所陈述/参考值±10%,无论哪个更高。
我们使用术语“生物活性变异体”来指与所提及部分(例如“母”分子或模板)相似但并不与其一致,并呈现用于本文所述的方法中的一种或多种的足够活性的部分(例如,代替所提及部分)。在一些情况下,所提及部分为天然存在的,而其生物活性变异体并非天然存在的。举例来说,在所提及部分为天然存在的核酸序列的情况下,其生物活性变异体可包括有限数目的非天然存在的核苷酸;可具有与其天然存在的对应部分不同的核酸序列;或可在其它方面不同于其天然存在的对应部分。举例来说,标记化试剂(tagmentationreagent)可包括“转座酶结合位点”(TBS),其与天然存在的TBS有所不同,但仍然保持很好地结合转座酶以使本发明组合物和方法发挥作用的能力。在所提及部分为天然存在的蛋白质的情况下,其生物活性变异体可包括有限数目的非天然存在的氨基酸;可具有与其天然存在的对应部分不同的肽序列;或可在其他方面不同于其天然存在的对应部分(例如借助于翻译后以不同方式修饰(例如其糖基化模式可有所不同)。
术语“解复用”意指通过数种手段中的任一种分析任何经测序的库片段上存在的可鉴别的序列标记信息,并将指配所述片段的来源的所述信息用于经修饰而并入所述可鉴别的序列标记的特定原始样品。
在一些实施例中,所述衔接子带有简并随机、半随机或非随机“可鉴别的序列标记”的子序列。在各个实施例中,可鉴别的序列标记在突触复合体之间可为相同或不同的。在一些实施例中,可鉴别的序列标记可用于“解复用”(即,计算上分组或分离)来源于单独样品的测序读序(即,经样品标记的读序)和/或来自某些池的测序读序(即,经池标记的读序)。
“试剂盒”是指用于递送进行本发明方法用的材料或试剂的任何递送系统。在分析的情形下,此类递送系统包括允许将反应试剂(合适容器中的例如探针、酶等)和/或支持材料(例如,缓冲液、用于进行分析的书面说明书等)从一个位置储存、传输或递送到另一个位置的系统。举例来说,试剂盒包括一个或多个封装体(例如箱),其含有本发明分析用的相关反应试剂和/或支持材料。此类内含物可共同或分别递送到目标接收者。举例来说,第一个容器可含有用于分析的酶,而第二个容器含有一个或多个探针。
“库”或“片段库”为衍生自一种或多种核酸样品的核酸分子的集合,其中已一般通过并入包含一个或多个引物结合位点和可鉴别的序列标记的末端衔接子序列来修饰核酸片段。
在本发明的上下文中且尤其就归一化来说,术语“有限的”、“有限量”等等用于描述构筑体,如包含活性转座酶的构筑体的量,所述构筑体在广泛范围的目标DNA浓度内将不会生成超过期望数量的“标记情况”(即,变位事件)。位点之间的距离暗中控制所产生的可测序片段的尺寸。举例来说,如果位点之间的平均距离为500,那么片段尺寸将平均为500。
当用于本文中时,术语“复用”或“复用化”一般是指多个核酸样品或其所衍生的库分子在一个池、试管或反应中的集合。
术语“归一化”是指一种如下方法:通过所述方法DNA样品之间的DNA测序读序的变化率大体上降低。尽管归一化通常通过利用样品稀释消耗大量劳动测量和调整每一DNA样品在库制备的各个阶段的浓度来实现,但本文中所述的实施例中的归一化是通过标记和汇集技术来实现。类似地,当用于本文中时,术语“归一化的”是指核酸分子或片段的异质混合物的一般期望属性,通过所述属性,由不同原始来源引起的片段丰度的变化率相对于原始来源中的丰度或量已大体上降低。当用于衍生或复用库(例如“归一化库”)的情形下时,术语是指如下属性:通过所述属性,来自多个样品中的每一个的库中的核酸片段的比例量(或百分比表述)小于原始样品的量的变化。尽管归一化通常通过利用样品稀释消耗大量劳动测量和调整每一DNA样品在库制备的各个阶段的浓度来实现,但本文中所述的实施例中的归一化是通过标记和汇集技术来实现。
一般以同义使用的“寡核苷酸”或“聚核苷酸”意指一种通过磷酸二酯键或其类似物连接的天然或修饰核苷单体的线性聚合物。术语“寡核苷酸”通常是指较短的聚合物,例如,包含约3到约100个单体的聚合物,并且术语“聚核苷酸”通常是指较长的聚合物,例如,包含约100个单体到数千个单体,例如10,000个单体或更多的聚合物。寡核苷酸包含长度通常在12到60个核苷酸,并且更通常18到40个核苷酸范围内的探针或引物。寡核苷酸和聚核苷酸可为天然或合成的。寡核苷酸和聚核苷酸包括脱氧核糖核苷、核糖核苷和其非天然类似物,如其异头形式;肽核酸(PNA)等等,其限制条件为其能够借助于单体与单体相互作用的常规模式特异性结合于靶基因组,如沃森-克里克型碱基配对(Watson-Crick type ofbase pairing)、碱基堆积、胡斯坦或反胡斯坦型碱基配对(Hoogsteen or reverseHoogsteen types of base pairing)等等。
核苷单体通常通过磷酸二酯键来连接。只要寡核苷酸由一连串字母,如“ATGC……”表示,那么应理解,除非另外指出,否则核苷酸从左到右呈5'到3'顺序,并且“A”表示脱氧腺苷,“C”表示脱氧胞苷,“G”表示脱氧鸟苷,“T”表示脱氧胸苷,并且“U”表示核糖核苷、尿苷。寡核苷酸通常包含四种天然脱氧核苷酸;然而,其还可包含核糖核苷或非天然的核苷酸类似物。一般熟习此项技术者应清楚,当寡核苷酸具有天然或非天然核苷酸时,其可用于本文所述的方法和工艺中。举例来说,在需要通过酶加以处理的情况下,通常需要仅由天然核苷酸组成的寡核苷酸。同样地,在酶具有针对活性的特定寡核苷酸或聚核苷酸基质要求,例如,单链DNA、RNA/DNA双螺旋体等等的情况下,则适于寡核苷酸或聚核苷酸基质的组合物的选择恰在普通技术人员的知识范围内,尤其在论文,如Sambrook等人,《分子克隆(Molecular Cloning)》,第二版(纽约冷泉港实验室(Cold Spring Harbor Laboratory,New York),1989)和类似参考文献的指导下。寡核苷酸和聚核苷酸可为单链或双链的。
寡核苷酸和聚核苷酸可以任选地包括一种或多种非标准核苷酸、核苷酸类似物和/或经修饰的核苷酸。经修饰的核苷酸的实例包括(但不限于)二氨基嘌呤、S2T、5-氟尿嘧啶、5-溴尿嘧啶、5-氯尿嘧啶、5-碘尿嘧啶、次黄嘌呤、黄嘌呤、4-乙酰基胞嘧啶、5-(羧基羟甲基)尿嘧啶、5-羧甲基氨甲基-2-硫代尿苷、5-羧甲基氨甲基尿嘧啶、二氢尿嘧啶、β-D-半乳糖基Q核苷、肌苷、N6-异戊烯基腺嘌呤、1-甲基鸟嘌呤、1-甲基肌苷、2,2-二甲基鸟嘌呤、2-甲基腺嘌呤、2-甲基鸟嘌呤、3-甲基胞嘧啶、5-甲基胞嘧啶、N6-腺嘌呤、7-甲基鸟嘌呤、5-甲基氨甲基尿嘧啶、5-甲氧基氨甲基-2-硫尿嘧啶、β-D-甘露糖基Q核苷、5'-甲氧基羧甲基尿嘧啶、5-甲氧基尿嘧啶、2-甲硫基-D46-异戊烯基腺嘌呤、尿嘧啶-5-氧基乙酸(v)、丁氧核苷(wybutoxosine)、假尿嘧啶、Q核苷、2-硫胞嘧啶、5-甲基-2-硫尿嘧啶、2-硫尿嘧啶、4-硫尿嘧啶、5-甲基尿嘧啶、尿嘧啶-5-氧基乙酸甲酯、尿嘧啶-5-氧基乙酸(v)、5-甲基-2-硫尿嘧啶、3-(3-氨基-3-N-2-羧丙基)尿嘧啶、(acp3)w、2,6-二氨基嘌呤等等。核酸分子亦可在碱基部分(例如在通常可用于形成具有补充核苷酸的氢键的一个或多个原子处和/或在通常不能够形成具有补充核苷酸的氢键的一个或多个原子处)、糖部分或磷酸骨架处得到修饰。
“寡核苷酸标记”、“序列标记”或“标记”意指连接到聚核苷酸的寡核苷酸或其复合物并用于鉴别和/或追踪反应中的聚核苷酸。通常,寡核苷酸标记连接到聚核苷酸的3'端或5'端以形成线性共轭物,所述有时在本文中被称作“标记聚核苷酸”或等效地被称作“经寡核苷酸标记的聚核苷酸共轭物”或“标记-聚核苷酸共轭物”。寡核苷酸标记可在尺寸和组成方面大幅变化;以下参考文献提供对选择适合于具体实施例的寡核苷酸标记组的指导:Brenner,美国专利第5,635,400号;Brenner等人,《美国国家科学院院刊(Proc.Natl.Acad.Sci.)》,97:1665;Shoemaker等人,《自然·遗传学(Nature Genetics)》,14:450,1996;莫里斯等人,欧洲专利公开案0799897A1;Wallace,美国专利第5,981,179号;和类似文献。
“聚合酶链反应(PCR)”是指通过对DNA的补充链进行同步引物延伸以使特定DNA序列试管内扩增的反应。换句话说,PCR为制得通过引物结合位点侧接的靶核酸的多个拷贝或复制物的反应,此类反应包含重复以下步骤一次或多次:(i)使靶核酸变性,(ii)将引物退火到引物结合位点,和(iii)在三磷酸核苷存在下通过核酸聚合酶延伸引物。通常,反应在热循环仪器中按照针对每一步骤优化的不同温度循环。特定温度、每一步骤的持续时间和步骤之间的变化率取决于本领域普通技术人员熟知的多种因素,例如通过以下参考文献示范的因素:McPherson等人编,《PCR:实用方法(PCR:A Practical Approach)》和《PCR2:实用方法(PCR2:A Practical Approach)》(IRL出版社,牛津(Oxford),分别为1991和1995)。举例来说,在使用Taq DNA聚合酶的常规PCR中,双链靶核酸可在大于90℃的温度下变性,引物在50℃到75℃范围内的温度下退火,并且引物在72℃到78℃范围内的温度下延伸。
PCR涵盖反应的衍生形式,包括(但不限于)RT-PCR、实时PCR、嵌套PCR、定量PCR、复用PCR等等。反应体积在数百纳升(例如200nL)到数百微升(例如200微升)范围内。“逆转录PCR”或“RT-PCR”意指前面是将靶RNA转化成补充单链DNA的逆转录反应的PCR,其随后扩增,例如Tecott等人,美国专利第5,168,038号。“实时PCR”意指随反应进行而监测反应产物,即扩增子的量的PCR。存在实时PCR的多种形式,其主要在用于监测反应产物的检测化学反应方面有所不同例如,Gelfand等人,美国专利第5,210,015号(“Taqman”);Wittwer等人,美国专利第6,174,670号和第6,569,627号(嵌入染料);Tyagi等人,美国专利第5,925,517号(分子信标)。实时PCR的检测化学反应在Mackay等人,《核酸研究(Nucleic Acids Research)》,30:1292-1305(2002)中得到综述。“嵌套PCR”意指两级PCR,其中第一PCR的扩增子变成使用一组新引物的第二PCR的样品,所述引物中的至少一个结合于第一扩增子的内部位置。如本文所用,关于嵌套扩增反应的“初始引物”意指用于生成第一扩增子的引物,并且“第二引物”意指用于生成第二或嵌套扩增子的一个或多个引物。“复用PCR”意指如下PCR:其中多个靶序列(或单一靶序列和一个或多个参考序列)在同一反应混合物中同时进行(参见例如,Bernard等人,《分析生物化学(Anal.Biochem.)》,273:221-228,1999(双色实时PCR))。通常,针对经扩增的每一序列采用不同引物组。“定量PCR”意指被设计成用于测量样品或试样中的一个或多个特定靶序列的丰度的PCR。定量PCR包括此类靶序列的绝对定量和相对定量。定量PCR的技术对于所述领域的普通技术人员为熟知的,如以下参考文献中所示范:Freeman等人(《生物技术(Biotechniques)》,26:112-126,1999;Becker-Andre等人(《核酸研究》,17:9437-9447,1989;Zimmerman等人(《生物技术》,21:268-279,1996;Diviacco等人(《基因(Gene)》,122:3013-3020,1992;Becker-Andre等人(《核酸研究》,17:9437-9446,1989);和类似文献。
术语“经池标记的”或“经集合标记的”是指由通过多种方法中的任一种将可鉴别的序列标记并入到经样品标记的核酸的池中产生的多个核酸分子。然而,在一些实施例中,所述“池标记”带有不可鉴别的序列标记,并实际上充当用于将核酸添加到经样品标记的核酸上的最小序列衔接子,使得池与特定DNA测序平台或扩增策略相容。
“引物”包括天然或合成的寡核苷酸,其能够在形成具有聚核苷酸模板的双螺旋体后充当核酸合成的起始点并从其3'端沿模板延伸以便形成经延伸的双螺旋体。在延伸过程中添加的核苷酸序列通过模板聚核苷酸的序列测定。引物通常通过DNA聚合酶延伸。引物的长度通常在3到36个核苷酸、5到24个核苷酸或14到36个核苷酸范围内。在某些方面中,引物为通用引物或非通用引物。引物对可侧接相关序列或相关序列组。引物和探针在序列中可为简并的。在某些方面中,引物结合邻近于靶序列,无论所述序列是否是待捕获用于分析的序列还是待拷贝的标记。
术语“样品标记的”和“经样品标记的”是指将可鉴别的序列标记并入到包含核酸样品的核酸分子中的过程和结果。
术语“鸟枪片段化(shotgun fragmentation)”描述使用机械DNA片段化方法(例如,声学、雾化或声处理)或利用
Figure BDA0001380396140000161
或其它非特异性核酸内切酶的酶促DNA片段化的各种方法。DNA断裂点通常为随机或半随机分布的,其跨越给定样品中所存在的靶DNA分子分布。除通过标记化生成的片段(单独描述)以外,由大多数鸟枪片段化方法生成的DNA片段在衔接子连接到片段末端之前需要末端修复和另外的酶促修饰,在多步骤方法中,其通常称作“鸟枪调适(shotgun adaptation)”。
我们使用术语“标记(tagment)”和“标记化(tagmentation)”来描述其中纯化转座酶在人工变位反应用作试剂以制备下一代DNA测序(NGS)的库的方法。在标记化中,转座酶首先负载有带有TBS的合成寡核苷酸,随后负载后的转座酶将活跃于分割靶DNA并同时所负载的寡核苷酸“标记”和TBS连接到靶DNA上的新生成的末端。
现较详细地参考本发明所涵盖的方法,在图1中,我们图示可在实践本发明的一个实施例中实施的步骤(又参见图2-6)。
在图1中,我们展示样品标记步骤100,其中多个样品的核酸样品经修饰而通过所属领域中已知的数种方法中的任一种将可鉴别的序列标记并入到由“鸟枪”片段化和调适产生的DNA片段中,吾人还可将所述可鉴别的序列标记称作DNA条码。可实施鸟枪片段化和调适,其优于例如,随机剪切、末端修复和衔接子连接,优于转座酶调适(即,标记化),或优于所属领域中已知的或所属领域中后续研发的任何相似方法。
样品标记步骤的关键特征在于添加到每一核酸样品中的可利用的标记的数量保持不变且与每一样品中的可利用的核酸的质量相比相对较低。这一限制条件可替代地例如通过调节时间、浓度或所属领域中已知的其它参数的反应条件来实现,以影响可由物理剪切、核酸酶消化或转座酶类调适获得的片段化,以便生成超过规定长度的片段。在本发明的情况下,规定或期望长度还可表示为样品中所存在的全部核苷酸的每个核苷酸所添加的标记数量。因此,所添加的标记数量可保持限制于每单元质量或摩尔浓度的1000个核苷酸低于一个或其它规定水平。举例来说,在其它实施例中,所添加的标记数量可保持为每单元质量或摩尔浓度的2000、3000、4000、5000或更多个核苷酸小于一个。
如通过各个示范性实施例将清楚,通过使样品中的样品标记率小于规定的每个核苷酸水平施加的限制条件在一些情况下可指示,并非样品中所存在的全部核酸分子都将必须接受样品标记情况。对本文中所主张的发明特别重要且相关的是,如果在相似条件下标记多个样品,并且允许每一样品中的核酸的量有所变化,而标记情况的数量限制在低于规定水平,那么样品之间的标记情况的数量将仍然大致相等。
仍参考图1,在对多个模板样品进行标记步骤100的情况下,将不同或特有聚核苷酸标记序列并入到多个样品中的每一单独样品中将视为有利的。如步骤100中所示,由此经修饰的模板样品随后被称为已‘经样品标记’或被称为具有‘样品标记’。样品标记指示样品集合中样品的特定位置,如微量滴定板上的微孔位置将另外被视为有利的。此外,将不同序列随机化的子区包括在样品标记(长度上有4-20个碱基)内以在后续DNA序列分析期间计数或检测由同一母分子衍生的DNA序列可为有利的。
仍参考图1,随后显示步骤102,其中已在步骤100中经样品标记的多个样品一起汇集成一种混合物,其随后含有多个经特别标记的模板样品,每一模板样品除了可变长度的另外未知的或部分未知的序列以外还含有多个已知标记序列中的一个,所述未知的序列侧接、邻近或为已知标记序列的上游或下游。
仍参考图1,随后显示步骤104,其中对来自步骤102的经汇集标记的样品进行池标记步骤,以使得包含经汇集的样品的原始核酸分子随机或半随机片段化,通过所属领域的从业者已知的数种方法中的任一种裂解或断裂成较小块,并且随后或同时用已知的单链或双链DNA衔接子序列调适,所述衔接子序列在一端或两端上侧接新衍生的片段。可用于步骤102的伴随的片段化方法可包括例如通过各种酶,如DNA核酸外切酶以及(或不存在)DNA切口酶进行酶促片段化,或通过机械力,如雾化(例如
Figure BDA0001380396140000181
)或声处理(例如
Figure BDA0001380396140000182
)进行剪切,或优选地通过修饰转座酶进行标记化,其准许同时片段化并将已知的寡核苷酸衔接子连接到片段化基质。在通过酶促或机械手段进行的片段化产生不规则核酸末端,即不同3'或5'突出端或不规则暴露的羟基的混合物的情况下,宜“补平(polish)”或3'-腺苷酸化具有酶的DNA,其准许衔接子后续连接到经片段化的模板样品。在此类情况下,经片段化、经补平和经3'-腺苷酸化的DNA可随后通过连接或其它适合方法处理以在多个新衍生片段的一端或两端上含有侧接DNA衔接子。步骤104的结果为由步骤102中的标记样品的池衍生的核酸分子的片段化库,其中在一端或两端上存在多个含有已知衔接子序列的片段。
通过对先前已通过本发明方法样品标记的样品进行池标记而实现的关键特征在于池标记的添加仅视池中所存在的核酸分子的可用质量而定,并且并不取决于经合并而制得池的单独样品分子的相对量,或代表所述同一池中的每一样品的单独核酸分子的所得频率。因此,当进行池标记时,存在使用本发明实施方法中的任一种将池标记添加到处于任何样品标记的规定数量的核苷酸中的位点的相等机率,所述样品标记代表包含池的分子中已经存在的任一原始样品。当使用本发明实施方法进行应用时,包含池的核酸分子的所得特征在于,对于池中任何核酸分子上存在的每一样品标记来说,存在通过池标记步骤添加的池标记中的一个将存在于规定距离的样品标记中的相等机率。
在本发明的各种实施方法中,经池标记的样品中的样品标记和池标记的间距可由调节添加池标记所处的条件来加以控制或限制。举例来说,通过改变时间、浓度或所属领域中已知的其它参数的反应条件来改变样品标记和池标记之间的距离,以影响通过物理剪切、核酸酶消化或转座酶类调适获得的片段化。在将公认为通常有利的是,一些优选的实施条件下,多个样品标记中的每一个和位于同一核酸分子上的相应邻近池标记的距离和相对定向将使得可通过PCR或其它相似手段进行所述同一核酸分子的序列子区域的扩增,所述子区域形成在样品标记和池标记之间。
参考图1的步骤106,通过合适方法扩增和富集步骤104的片段化和池标记。合适方法可包括例如聚合酶链反应(PCR)。所属领域的一般技术人员将认识到使用对并入原始样品标记步骤100中的序列具有特异性,以及对池标记步骤104具有特异性的寡核苷酸引物的扩增可促进对库片段的的扩增和富集,所述库片段除了通过步骤104中所添加的衔接子在一端处侧接的原始模板分子的随机片段以外,在相对端处还含有用于鉴别步骤100的样品的聚核苷酸标记序列。
在其它实施例中,在已由原始多个样品中的每一个产生多个不同池,所述池中的每一个为由所述步骤100、102、104和106产生的经汇集、经片段化、经调适、经扩增和富集的库的情况下,将公认为尤其有利的是,步骤106的扩增和富集中所用的寡核苷酸引物可经修饰而并入一个或多个特别鉴别的聚核苷酸标记序列(即,DNA条码或可鉴别的序列标记),其可进一步与由步骤106产生的合适片段的整个池相关。在寡核苷酸引物由此经修饰的情况下,通过使用所述池-标记序列知晓的所得池可用于特别鉴别针对池中的任何样品片段的来源的池。
在其它实施例中,在已由多个多个样品产生多个不同池的情况下,将公认为尤其有利的是,池标记序列作为任何池标记步骤的一部分用于并入一个或多个特别鉴别的池-标记序列(即,DNA条码或可鉴别的序列标记),其可进一步与所存在的合适核酸分子的整个池相关。
仍参考图1,随后显示步骤108,其中经扩增和富集的库池(每一库池标记有可鉴别的池标记)经汇集而形成一池标记池,即,经汇集的池。
仍参考图1,随后显示步骤112,其中已根据需要另外经尺寸分级以用于进一步处理的作为所得所得步骤10的统一库,经历大规模平行测序(MPS),如
Figure BDA0001380396140000201
HiSeq或MiSeq测序或Ion
Figure BDA0001380396140000202
Proton测序。在一些模式下,步骤100到106中所用的衔接子和寡核苷酸序列的特定设计或选择将影响测序平台的优选选择。在优选模式下,将实现短读序测序以便测定来自步骤110中所衍生的库的多个不同片段的核苷酸序列。此外,同样在优选模式下,以如下方式进行短读序测序:所述方式使得多个片段中的每一片段的序列测定引起测定由步骤100或步骤106产生的每一片段中任何或全部聚核苷酸标记序列的序列,以及衍生所述片段的原始模板的已知或未知区域的序列。
仍参考图1,随后显示步骤114,其中根据来源于样品标记步骤100的聚核苷酸标记序列的组合或池扩增和富集步骤106处的池标记对步骤112中所获得的单独测序读序进行归类或解复用。在一些实施例中,聚核苷酸标记序列的精确组合可用于映射或关联多个单独测序读序,其借助于步骤100处引入的标记针对原始样品集合中的样品的位置具有相同标记,或借助于步骤106处引入的标记针对代表整个集合的片段的池具有相同标记。以此方式,对于展现在多个原始集合中的每一个中的多个原始样品中的每一个来说,可衍生多个单独测序读序,其特别映射或对应于一个原始样品。此外,紧邻的前述多个读序中的单独测序读序中的每一个含有输入到步骤100的原始核酸模板样品的可能未知的区域的其他序列测定。
仍参考图1,随后显示步骤116,其中对由步骤114产生的多个测序读序进行组合、分析或以计算方式处理以衍生对应于原始模板核酸中的区域的组合序列,所述测序读序与模板样品的一个原始集合中的一个特定模板核酸样品相关。组合方法可通过各种现有DNA序列组合工具,如Velvet、MIRA、SSAKE或ALLPATHS以计算方式完成。如果衍生读序的重新组合并非期望的,那么也可使用各种现有DNA序列比对工具,如BWA、BOWTIE、MAQ等对衍生读序与核苷酸序列参考物进行比较和/或映射。
图2显示图1中所示的样品标记步骤100的代表性实施例。具体言之,示意性图解显示多个核酸分子202标记有固定数量的样品标记201。模板DNA 202可由线性或圆形DNA模板,例如cDNA、扩增子或质粒;或较长DNA模板,如基因组DNA、线粒体DNA、细菌人工染色体(BAC)或任何任意的相关DNA分子以不同方式组成。在不同实施例中,样品标记步骤可包含所属领域中已知的一组方法中的任一种,通过所述方法,核酸样品经修饰或处理而并入特定聚核苷酸标记序列;聚核苷酸标记序列本身可包含一个或多个衔接子区域,引物结合位点或可鉴别的序列标记,其又通常称作DNA条码。
仍参考图2,在一些实施例中,在样品标记步骤期间并入的标记204可由同一已知的寡核苷酸序列或不同寡核苷酸序列组成。此外,标记204可由具有相同或不同长度的寡核苷酸序列组成,其中所述长度可在4到10bp、或8到15bp、或10到20bp范围内。
现参考图3,为了强调本发明的归一化方面,显示较详细的过程图,通过所述方法,具有不同原始浓度的样品得到样品标记和汇集。具体言之,显示样品标记反应的三种情况(301、302、303),其中所描绘的每一样品中所存在的核酸(DNA)的量由低(例如1ng/μL)向中等(例如10ng/μL)再向高(100ng/μL)变化。对于三种样品标记反应情况中的每一种,使用相对固定和受限的的样品标记条件以使得每一规定单元质量的核酸的标记数量受到限制。所得经样品标记的样品304、305和306由此具有相似数量的样品标记情况,但仍然具有不同量的核酸分子。多个经标记的DNA样品经汇集而形成样品池31010,其中通过方法(如图2的方法)制造的多个可能以不同方式标记的DNA样品被合并成一种异质混合物。这一池的重要特征在于尽管其包含与每一样品中的核酸分子的质量大致成比例的不同的来自每一样品的核酸分子的质量,但来自池中所存在的每一样品的样品标记情况的数量为固定和归一化的。图3的过程为图1的样品汇集步骤102的实施例。
现参考图4,显示一个方法的实施例,其中通过池标记步骤对由不同样品产生并经汇集而形成经样品标记的池402为池402的多个经样品标记的核酸分子进行汇集标记(标记为403),以衍生多个核酸分子(404),其样品标记、池标记或这两者的不同情况。由于池标记步骤是在经聚集的样品标记的样品上进行,因此池标记的分布为均一的,并且池标记位点产生在样品标记位点的规定距离内的机率与池中的全部不同样品标记相等。
在各个实施例中,如图3和图4中所描绘的池标记可以趋向于形成较大或较小片段的方式,例如通过滴定用于产生所述片段的碎片化酶的量来进行。所述碎片化酶可包括例如DNA核酸外切酶I、或Sau3AI核酸内切酶、或dsDNA
Figure BDA0001380396140000221
或经修饰的DNA转座酶,如
Figure BDA0001380396140000222
酶。通过所述方法,潜在地有利于使不同测序应用能够通过变更本发明的规定方法或任何相关试剂盒来指明,以便形成较大或短的线性dsDNA片段。
现参看图5,显示一种如下方法:其中由图4中所描绘的方法产生的片段的库可针对保留一端上的样品标记序列和通过鸟枪衔接子侧接的另一端上的一部分潜在未知的序列的片段得到扩增和富集。
现参考图5,显示一个代表性实施例,其中可用不同池标记处理经样品标记的核酸分子的多个不同池,以便获得大量的样品的可能性构型和池标记组合,使得可对更多样品进行并行测序。
现参考图6,显示对由图5中所示的步骤衍生的多个库片段进行测序的绘图,其中每一库片段由样品标记602和一个或多个池标记600和606以及未知的插入区域604组成。由样品标记和池标记获得序列信息使得未知的插入区域的序列与原始培养盘或集合中的样品相关联和/或指配于所述样品。
以下为本发明方法的某些实施例,其通过实例说明方法对特定DNA测序任务的应用。
实例
实例1:本实例说明本发明中的方法对来自96孔培养盘中所含的伯克霍尔德菌菌株的基因组DNA(gDNA)样品的集合的应用。进行以下实验以制造针对每一伯克霍尔德菌属分离株的高质量gDNA序列。使用如下文所述的修饰标记化反应设定对每一gDNA样品进行样品标记:
Figure BDA0001380396140000241
SEQ ID NO:1
5'-AGACGTGTGCTCTTCCGATCTCAACCCGAACCGAGTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3'
SEQ ID NO:2:5'-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-3'
全部96个反应都接受同一通用反向引物(SEQ ID NO:2)/用同一通用反向引物(SEQ ID NO:2)进行,并且每一反应还接受不同的样品标记正向引物/用不同的样品标记正向引物进行(即,96种特有的样品标记引物之一,其实例显示在SEQ ID NO:1中)。使用热循环仪在同一反应中进行标记化和利用PCR的样品标记,所述热循环仪编程如下:
Figure BDA0001380396140000242
完成PCR后,将来自96孔PCR培养盘的每一孔的10μL经样品标记的PCR产物(样品标记有可鉴别的测序标记)汇集到单一1.5mL微量离心试管(
Figure BDA0001380396140000251
DNA LoBind)中,并随后根据制造商的说明书使用0.7×体积的
Figure BDA00013803961400002510
XP磁珠(贝克曼库尔特公司(BECKMAN COULTER)
Figure BDA0001380396140000254
)加以纯化。
使用伴随PicoGreen dsDNA分析试剂盒
Figure BDA0001380396140000255
的荧光定量对所汇集的纯化产物进行定量并在
Figure BDA0001380396140000256
Infinite M200微量板读取器上获得读数。随后用10mMTris-HCl,pH 8.0将定量后的经汇集的、经样品标记的产物的浓度调整到0.8ng/μL,并使用以下反应设定条件对其进行第二轮标记化:
Figure BDA0001380396140000257
将反应聚集在8孔PCR条管的0.2mL孔中,并用塑料盖加以密封。使用
Figure BDA0001380396140000258
iCycler DNA热循环仪使标记化反应在55℃下保温10分钟,随后保持在4℃下。
第二标记化之后,使用池-标记引物对产物进行PCR扩增,其使用以下反应条件:
Figure BDA0001380396140000259
SEQ ID NO:3:5'-CAAGCAGAAGACGGCATACGAGATTCGCCTTAGTGACT
GGAGTTCAGA CGTGTGCTCTTCCGATC-3'
SEQ ID NO:4:5'-AATGATACGGCGACCACCGAGATCTACACTAGATCGC
TCGTCGGCAGC GTCAGATGTGTATAAGAGACAG-3'
对于每一PCR反应,我们使用池标记正向引物和索引反向引物(Index ReversePrimer)的不同组合。将PCR反应集中在先前用于标记化步骤的相同8孔PCR条管中,并随后循环如下:
Figure BDA0001380396140000261
扩增之后,随后将样品混合并汇集到一个1.5mL微量离心试管中以形成池。使用制造商所建议的结合和洗涤方案用60μL
Figure BDA0001380396140000262
Figure BDA0001380396140000263
XP磁珠从池混合物中纯化并大小选择100μL。
随后将经纯化的、经大小选择的集合池稀释10,000倍并使用
Figure BDA0001380396140000264
库定量试剂盒通过qPCR加以分析。随后使用
Figure BDA0001380396140000265
测序系统(150循环v3化学反应)以以下读长配置针对总共164次测序循环对定量后的库进行测序:
Figure BDA0001380396140000271
读序的解复用:使用定制处理脚本,分析由
Figure BDA0001380396140000272
测序仪获得的测序数据以获得由库的测序产生的大约1百万个读序中的每一个的集合和样品标记。FASTQ文件由使用这一解复用方法处理的96个样品中的每一个生成。
产生序列的孔的%的评估:图7显示由测序数据衍生的读序的相对百分率(占总体的%),所述指数可指配于使用所述方法处理的原始培养盘上的96个伯克霍尔德菌属样品中的每一个。
实例2
本实例说明本发明方法供用于由质粒克隆的集合衍生长组合读序的用途,所述质粒克隆衍生自非洲爪蟾属开放阅读框(ORF)集合。质粒克隆的十(10)个96孔培养盘处理如下。
使用这些反应设定条件来进行样品鸟枪调适:
Figure BDA0001380396140000273
使用不同的经样品标记的正向引物(例如SEQ ID NO:1)将反应集中于96孔PCR培养盘的每一孔中,但在所有96个反应中相同通用反向引物。随后在如下的热循环仪中使反应保温:
Figure BDA0001380396140000281
完成热循环后,将来自PCR培养盘的每一孔的10μL经样品标记的PCR产物汇集到单一1.5mL微量离心试管(
Figure BDA0001380396140000282
DNA LoBind)中,并使用0.7×体积的
Figure BDA0001380396140000283
XP珠粒加以纯化。
使用伴随PicoGreen dsDNA分析试剂盒
Figure BDA0001380396140000284
的荧光定量对所纯化的经样品标记的PCR产物进行定量并在
Figure BDA0001380396140000285
Infinite M200微量板读取器上获得读数。随后在Tris-HCl,pH 8.0中将定量后的产物的平均浓度调整到0.8ng/μL。
随后使用以下反应设定条件对所纯化的经样品标记的PCR产物进行第二标记化:
Figure BDA0001380396140000286
将反应聚集在8孔PCR条管的0.2mL孔中,并用塑料盖加以密封。使用
Figure BDA0001380396140000291
iCycler DNA热循环仪使标记化反应在55℃下保温10分钟,随后保持在4℃下。
随后使用以下PCR反应设定条件使用经池标记的引物对第二标记化反应的产物进行扩增:
Figure BDA0001380396140000292
为了由第二标记化的产物形成十(10)个经池标记的子库,使用来自克隆非洲爪蟾属ORF的不同的96孔培养盘的96种经样品标记的DNA的池作为模板来设定每一反应,并随后每一反应接受经池标记的正向和反向引物的特有组合以明确地指配每一96孔培养盘的测序读序。将这些十(10)个PCR反应聚集在如用于第二标记化步骤的相同的8孔PCR条管中,并在密封之后,使反应循环如下:
Figure BDA0001380396140000293
扩增之后,将十(10)个反应(现表示不同的经池标记的子库)汇集到单一1.5mL微量离心试管中以形成单一经池标记的库(即,所汇集的样品的池)。使用制造商所建议的结合和洗涤方案用60μL
Figure BDA0001380396140000301
XP(贝克曼库尔特公司
Figure BDA0001380396140000302
)从经池标记的库中大小选择100μL,一般富集大小为300bp及更高的片段。
将经纯化的、经大小选择的、经池标记的库稀释10,000倍并使用
Figure BDA0001380396140000303
库定量试剂盒通过qPCR加以分析。随后使用
Figure BDA0001380396140000304
测序系统以以下读长配置针对总共127次测序循环对定量后的库进行测序:
Figure BDA0001380396140000305
读序的解复用:使用定制处理脚本,分析由
Figure BDA0001380396140000306
测序仪获得的测序数据以获得由库的测序产生的大约1百万个读序中的每一个的池标记和样品标记。FASTQ文件由来自使用这一解复用方法处理的10个培养盘中的每一个的96个样品中的每一个生成。
读序组合:随后使用软件程序
Figure BDA0001380396140000307
对解复用读序进行重新组合。这一分析的输出结果为一系列组合后的连续序列(重叠群),其对应于原始集合中模板DNA中的每一个的全长序列。图8显示源于这一方法的长度为大约3.5kb的新生组合体的实例。
实例3
本实例通过对由96个DNA标准物(全部具有相等DNA质量)生成的DNA测序读序的数量与由同一标准DNA的两倍稀释系列产生的读序的数量进行比较来说明本发明的归一化特性。
如Picelli等人(《基因组研究(Genome Research)》15:2033-2040,2014)所述,将极具活性的Tn5转座酶的储备浓度纯化为35μM。
制备两个DNA标准物的96孔培养盘(2-log DNA梯度0.1-10.0kb,NEW ENGLAND
Figure BDA0001380396140000311
):对照培养盘A(全部96个孔中为10ng DNA标准物)和稀释培养盘B(DNA标准物的两倍连续稀释物,其以列H的所有12个孔的128ng标准输入DNA开始,并以列B中的所有12个孔的2ng DNA标准物结束。使来自培养盘A和B的每一样品独立地标记有可鉴别的序列标记,其使用以下反应设定条件:
将以下组分添加到每一96孔PCR培养盘的每一孔中以单独地标记每一样品:
Figure BDA0001380396140000312
SEQ ID NO:5:5'-CAAGCAGAAGACGGCATACGAGATACTAACTGGTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3'
SEQ ID NO:6:5'-磷酸-CTGTCTCTTATACACATCT-3'
用粘性箔密封培养盘之后,通过标记化将序列标记直接添加到样品中。使用
Figure BDA0001380396140000313
Tetrad DNA热循环仪使反应保温如下:55℃持续15分钟,72℃持续5分钟,并随后保持在4℃下。
随后,每孔添加2μL的60单位/毫升蛋白酶K溶液
Figure BDA0001380396140000314
和3μL的终止溶液(167mM MgCl2;6.7%SDS)。用粘性箔重新密封96孔PCR培养盘;保温于热循环仪中处于50℃下15分钟;并随后在68℃下保温5分钟。
在完成蛋白酶K反应之后,将来自培养盘A(DNA对照物)的全部经样品标记的产物汇集于单一微量离心试管中;并且将来自培养盘B(DNA稀释系列)的全部经样品标记的产物汇集于分开的微量离心试管中。随后根据作业说明书如下分开纯化来自培养盘A和来自培养盘B的所汇集的经样品标记的产物:将相等体积的MAGwiseTM顺磁纯化珠粒(SEQWELLTM)添加到每一样品中并在磁性分离之前结合5分钟;随后,用80%乙醇洗涤珠粒两次,并随后再悬浮于10mM Tris-HCl,pH 8.0中以洗脱所纯化的DNA。最后,在保温5分钟以以磁性方式从悬浮液去除珠粒之后,将约30μL的洗脱后的每一DNA转移到干净的1.5mL微量离心试管(
Figure BDA0001380396140000321
DNA LoBind)。
在纯化两个经样品标记的产物之后,将全部体积的每一洗脱液(30μL)转移到8孔PCR条管中的单独孔,并添加以下组分以设定第二标记化:
Figure BDA0001380396140000322
SEQ ID NO:7:5'-AATGATACGGCGACCACCGAGATCTACACTAGATCGCTCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-3'
SEQ ID NO:8:5-磷酸-CTGTCTCTTATACACATCT-3
用粘性箔密封8孔PCR条管之后,将经池标记的反应如下保温于热循环仪中:
Figure BDA0001380396140000331
随后将12μL终止溶液添加到每一池标记反应中。重新密封之后,使8孔PCR条管返回到热循环仪,并保温如下:
Figure BDA0001380396140000332
将相等体积的MAGwiseTM顺磁纯化珠粒(SEQWELLTM)添加到每一样品中,并在磁性分离之前保温5分钟。随后,根据操作说明书用80%乙醇洗涤珠粒两次,并随后再悬浮于10mM Tris-HCl中以洗脱所纯化的DNA。最后,在保温5分钟以以磁性方式从悬浮液去除珠粒之后,将约30μL的洗脱后的每一DNA转移到干净的1.5mL微量离心试管(
Figure BDA0001380396140000333
DNALoBind)。
如下文所述设定两种衔接子填充反应以延伸标记化之后残留在DNA片段上的内嵌式3'-突出端。填充反应集中于8孔PCR条管中。
Figure BDA0001380396140000334
在密封8孔PCR条管之后,将反应保温于热循环仪中处于72℃下10分钟。随后在衔接子填充之后立即将以下引物直接添加到反应中以使用PCR开始进行库扩增:
Figure BDA0001380396140000341
SEQ ID NO:9:5'-AATGATACGGCGACCACCGA-3'
SEQ ID NO:10:5'-CAAGCAGAAGACGGCATACGA-3'
重新密封8孔PCR条管,随后使其返回到热循环仪以使用以下PCR循环条件进行库扩增:
Figure BDA0001380396140000342
在PCR扩增之后,使两种库(对应于来自对照物和稀释系列的样品)各与0.75×体积的MAGwiseTM顺磁纯化珠粒(SEQWELLTM)混合以对每一库进行大小选择,并在磁性分离之前保温5分钟。随后,根据制造商的操作说明书用80%乙醇洗涤珠粒两次,并随后再悬浮于10mM Tris-HCl中以洗脱所纯化的DNA。最后,在保温5分钟以以磁性方式从悬浮液去除珠粒之后,将约30μL的每一纯化后的库转移到两个干净的1.5mL微量离心试管(
Figure BDA0001380396140000343
DNA LoBind)。通过凝胶电泳分析每一经大小选择的库的等分试样,并且在检测后,片段长度分布在长度上呈现为大约300bp及更高。
随后将所纯化的、经大小选择的库稀释10,000倍并使用
Figure BDA0001380396140000344
Figure BDA0001380396140000345
库定量试剂盒通过qPCR加以分析。随后使用
Figure BDA0001380396140000346
Figure BDA0001380396140000351
500测序系统以高输出流动池配置并使用以下读长配置针对总共318次循环对定量后的库进行测序:
Figure BDA0001380396140000352
读序的解复用:使用定制处理脚本,分析由
Figure BDA0001380396140000353
500测序仪获得的测序数据以获得由测序每一库产生的大约5亿个读序中的每一样品中的样品特异性和培养盘特异性的可鉴别序列标记。FASTQ文件也由使用这一解复用方法处理的2个培养盘中的每一个上的96个样品生成。
评定归一化结果:将来自每一经稀释的DNA样品的解复用测序读序的数量与由固定量的对照DNA(设定为100%)产生的读序的平均数目进行比较。这一分析的输出结果显示在图9中。平均读序数量相对于大于60倍的一系列DNA输入质量(2ng到128ng的DNA输入)属于窄1.4倍范围内。实例3中的结果突显出可使用本发明方法实现显著归一化作用。

Claims (22)

1.一种由多个样品获得DNA片段的迭代标记库的方法,所述方法包含:
(a)使所述多个样品的第一样品和第二样品中的靶核酸分子分别与均一和限定量的各自包含可鉴别的样品序列标记的第一构筑体和第二构筑体接触,其中所述第一构筑体中的所述样品序列标记与所述第二构筑体中的所述样品序列标记不同,并且所述接触会进行一定持续时间并在如下条件下进行:其中所述第一构筑体中的所述样品序列标记被整合到所述第一样品中的所述核酸分子中,并且所述第二构筑体中的所述样品序列标记被整合到所述第二样品中的所述核酸分子中,并且其中所述第一构筑体和第二构筑体的均一和限定量使得每一样品中整合的可鉴别的样品序列标记的数目对于第一样品和第二样品来说是大致相等的,而与每一样品中存在的核酸的量或浓度无关;
(b)汇集所述第一样品和第二样品以产生经样品标记的核酸分子的池;
(c)使所述经样品标记的核酸分子的池与包含池标记的第三构筑体接触,其中所述接触会进行一定持续时间并在如下条件下进行:其中所述池标记以按照核苷酸的规定频率被整合到所述池的核酸分子中,由此生成经样品标记的、经池标记的核酸分子的池;以及
(d)用引物扩增所述经样品标记的、经池标记的核酸分子,由此生成DNA片段的迭代标记库,所述引物包含(i)与所述经样品标记的、经池标记的核酸分子中的核酸序列互补的核酸序列,和(ii)适用于下一代测序的核酸序列。
2.根据权利要求1所述的方法,其中所述池标记包含可鉴别的序列标记。
3.根据权利要求1所述的方法,其中所述多个样品中的至少一个样品包含基因组DNA。
4.根据权利要求1所述的方法,其中所述多个样品中的至少一个样品包含cDNA。
5.根据权利要求1所述的方法,其中所述多个样品中的至少一个样品包含合成性DNA。
6.根据权利要求1所述的方法,其中所述多个样品中的至少一个样品包含在病毒、细菌、酵母菌、真菌、原虫、昆虫、鱼、禽、哺乳动物或植物中天然存在的DNA。
7.根据权利要求1所述的方法,其中所述多个样品包含2-9,600个样品。
8.根据权利要求1所述的方法,其中所述第一构筑体和第二构筑体包含衔接子,所述衔接子随后连接到所述靶核酸分子。
9.根据权利要求1所述的方法,其中所述第一构筑体和/或第二构筑体包含转座酶或其生物活性变异体,所述转座酶或其生物活性变异体随后将所述可鉴别的样品序列标记引入到所述靶核酸分子中。
10.根据权利要求9所述的方法,其中在所述样品序列标记被整合到所述靶核酸分子中的条件下,所述样品序列标记的量比随后用于接触所述经样品标记的核酸分子的池的所述池标记的量低至少两倍。
11.根据权利要求1所述的方法,其中所述第三构筑体包含衔接子,所述衔接子随后连接到所述经样品标记的核酸分子。
12.根据权利要求1所述的方法,其中所述第三构筑体包含转座酶或其生物活性变异体,所述转座酶或其生物活性变异体随后将所述池标记引入到所述经样品标记的核酸分子中。
13.根据权利要求1所述的方法,其中所述第一构筑体、第二构筑体或第三构筑体独立地包含特异性结合捕获剂的部分。
14.根据权利要求13所述的方法,其中所述特异性结合捕获剂的部分为特异性结合所述捕获部分中的蛋白质或核酸序列的核酸序列。
15.根据权利要求13所述的方法,其中所述特异性结合捕获剂的部分为已被生物素化或洋地黄毒苷化的部分。
16.根据权利要求13所述的方法,其进一步包含如下步骤:通过将经样品标记的和/或经池标记的核酸分子群暴露于有限量的所述捕获剂来使所述群归一化。
17.根据权利要求1所述的方法,其中所述池标记以按照核苷酸的规定频率被整合到所述经汇集的、经样品标记的核酸分子中。
18.根据权利要求17所述的方法,其中所述规定频率在平均100-20,000个核苷酸之间产生池标记,其远离所述池中的所述经样品标记的核酸分子上的每一样品序列标记。
19.根据权利要求1所述的方法,其进一步包含对所述DNA片段的迭代标记库进行下一代测序。
20.根据权利要求1所述的方法,其中所述方法产生归一化库。
21.根据权利要求1所述的方法,其中所述第一构筑体和第二构筑体的均一和限定量使得所述多个样品的第一样品和第二样品的核酸分子中对于每500个核苷酸碱基对来说整合了小于一个的样品序列标记。
22.根据权利要求1所述的方法,其中所述第一构筑体和第二构筑体的均一和限定量使得所述多个样品的第一样品和第二样品的核酸分子中对于每1000个核苷酸碱基对来说整合了小于一个的样品序列标记。
CN201680010540.5A 2015-01-16 2016-01-16 Dna集合的归一化迭代条形码和测序 Active CN107250381B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201562104107P 2015-01-16 2015-01-16
US62/104,107 2015-01-16
PCT/US2016/013753 WO2016115550A1 (en) 2015-01-16 2016-01-16 Normalized iterative barcoding and sequencing of dna collections

Publications (2)

Publication Number Publication Date
CN107250381A CN107250381A (zh) 2017-10-13
CN107250381B true CN107250381B (zh) 2022-03-08

Family

ID=56406512

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680010540.5A Active CN107250381B (zh) 2015-01-16 2016-01-16 Dna集合的归一化迭代条形码和测序

Country Status (5)

Country Link
US (1) US20180010120A1 (zh)
EP (1) EP3245304B1 (zh)
CN (1) CN107250381B (zh)
HK (1) HK1247252A1 (zh)
WO (1) WO2016115550A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018081113A1 (en) 2016-10-24 2018-05-03 Sawaya Sterling Concealing information present within nucleic acids
JP7047373B2 (ja) 2017-12-25 2022-04-05 トヨタ自動車株式会社 次世代シーケンサー用プライマー並びにその製造方法、次世代シーケンサー用プライマーを用いたdnaライブラリー並びにその製造方法、及びdnaライブラリーを用いたゲノムdna解析方法
WO2023056366A1 (en) * 2021-09-29 2023-04-06 Seqwell, Inc. Kits and methods for preparation of nucleic acid libraries for sequencing

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102439177A (zh) * 2009-04-02 2012-05-02 弗卢伊蒂格姆公司 用于对目标核酸进行条形码化的多引物扩增方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5763239A (en) * 1996-06-18 1998-06-09 Diversa Corporation Production and use of normalized DNA libraries
US9080211B2 (en) * 2008-10-24 2015-07-14 Epicentre Technologies Corporation Transposon end compositions and methods for modifying nucleic acids
ES2595433T3 (es) * 2010-09-21 2016-12-30 Population Genetics Technologies Ltd. Aumento de la confianza en las identificaciones de alelos con el recuento molecular
EP2675819B1 (en) * 2011-02-18 2020-04-08 Bio-Rad Laboratories, Inc. Compositions and methods for molecular labeling
US9005935B2 (en) * 2011-05-23 2015-04-14 Agilent Technologies, Inc. Methods and compositions for DNA fragmentation and tagging by transposases
US9752176B2 (en) * 2011-06-15 2017-09-05 Ginkgo Bioworks, Inc. Methods for preparative in vitro cloning
EP2898096B1 (en) * 2012-09-21 2024-02-14 The Broad Institute, Inc. Methods for labeling of rnas
EP3597772A1 (en) * 2013-04-17 2020-01-22 Agency For Science, Technology And Research Method for generating extended sequence reads

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102439177A (zh) * 2009-04-02 2012-05-02 弗卢伊蒂格姆公司 用于对目标核酸进行条形码化的多引物扩增方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
利用高效转座质粒构建羽毛分解菌 Streptinomonas maltophilia YHJ-1 的突变体库;黄亦钧等;《农业生物技术学报》;20091231;第17卷(第5期);第931-936页 *

Also Published As

Publication number Publication date
CN107250381A (zh) 2017-10-13
EP3245304A4 (en) 2018-10-24
EP3245304B1 (en) 2021-06-30
EP3245304A1 (en) 2017-11-22
US20180010120A1 (en) 2018-01-11
WO2016115550A1 (en) 2016-07-21
HK1247252A1 (zh) 2018-09-21

Similar Documents

Publication Publication Date Title
US11725241B2 (en) Compositions and methods for identification of a duplicate sequencing read
US20210071171A1 (en) Compositions and methods for targeted nucleic acid sequence enrichment and high efficiency library generation
US10704091B2 (en) Genotyping by next-generation sequencing
AU2013325107C1 (en) Method of producing a normalised nucleic acid library using solid state capture material
JP2017533709A (ja) データの速度および密度を増大させるための多数のプライマーからのシーケンシング
US20170145493A1 (en) Nucleic acid ligation systems and methods
CN107250381B (zh) Dna集合的归一化迭代条形码和测序
JP2022513343A (ja) 次世代シーケンスにおいて低サンプルインプットを扱うための正規化対照
CN114555821B (zh) 检测与dna靶区域独特相关的序列
WO2020232081A2 (en) Capture and analysis of target genomic regions
KR20220130592A (ko) 핵산의 정확한 병렬 정량분석을 위한 고감도 방법
Rao et al. Recent trends in molecular techniques for food pathogen detection
US20240141426A1 (en) Compositions and methods for identification of a duplicate sequencing read
CA3208896A1 (en) Highly sensitive methods for accurate parallel quantification of variant nucleic acids

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