CN105934522B - 用于多核苷酸测序的组合物和方法 - Google Patents
用于多核苷酸测序的组合物和方法 Download PDFInfo
- Publication number
- CN105934522B CN105934522B CN201480064725.5A CN201480064725A CN105934522B CN 105934522 B CN105934522 B CN 105934522B CN 201480064725 A CN201480064725 A CN 201480064725A CN 105934522 B CN105934522 B CN 105934522B
- Authority
- CN
- China
- Prior art keywords
- polynucleotide
- sequences
- sequence
- target polynucleotide
- translocation
- 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
Links
Images
Classifications
-
- 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/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C07—ORGANIC CHEMISTRY
- C07K—PEPTIDES
- C07K14/00—Peptides having more than 20 amino acids; Gastrins; Somatostatins; Melanotropins; Derivatives thereof
- C07K14/195—Peptides having more than 20 amino acids; Gastrins; Somatostatins; Melanotropins; Derivatives thereof from bacteria
- C07K14/35—Peptides having more than 20 amino acids; Gastrins; Somatostatins; Melanotropins; Derivatives thereof from bacteria from Mycobacteriaceae (F)
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/26—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating electrochemical variables; by using electrolysis or electrophoresis
- G01N27/416—Systems
- G01N27/447—Systems using electrophoresis
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Organic Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Zoology (AREA)
- Engineering & Computer Science (AREA)
- Wood Science & Technology (AREA)
- Genetics & Genomics (AREA)
- Biophysics (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- Biotechnology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- Microbiology (AREA)
- Medicinal Chemistry (AREA)
- Gastroenterology & Hepatology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Apparatus Associated With Microorganisms And Enzymes (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Enzymes And Modification Thereof (AREA)
Abstract
利用用目标多核苷酸通过孔的易位的分步易位步骤来表征目标多核苷酸,包括表征目标多核苷酸的序列的方法和组合物。
Description
相关申请的引用
本专利申请主张2013年11月26日提交的标题为“用于多核苷酸测序的组合物和方法”的美国临时专利申请No.61/909,316的权益,该专利申请的全部内容作为参考并入本文。
序列表
本专利申请包含以ASCII格式电子提交的序列表,并且该序列表以其全部内容作为参考并入本文。创建于2014年11月25日的所述ASCII拷贝的名称为12957-139-228_SL.txt和大小为19,778字节。
发明背景
本发明公开一般地涉及表征目标多核苷酸的方法和组合物,其包括表征目标多核苷酸的序列。
由于多核苷酸(例如,DNA或RNA)中编码的信息对医学和生命科学非常重要,因此需要快速且经济地对多核苷酸测序。目前,商品化测序技术需要样品和文库制备,两者均很费力。此外,对于多种应用,读取数据比预期更慢。因此,通量受限并且成本相对较高。纳米孔测序代表了正在开发的用于快速且廉价地对目标多核苷酸测序的一种新方法。
纳米孔测序利用了可以为离子电流提供通道的纳米孔。电泳驱动多核苷酸通过纳米孔,并且由于多核苷酸穿过纳米孔,因此它降低了通过纳米孔的电流。每个通过的核苷酸或一系列核苷酸获得了特征电流,并且对电流电平的记录对应于多核苷酸序列。由于一些电流电平取决于多个核苷酸(通常3-4个),因此仍需要改善目前工艺水平以改善准确度。由于多核苷酸易位通过纳米孔所获得的电流电平的任何其它信息可以提供优势,如形状和持久时间。
纳米孔测序的常见问题是多核苷酸通过纳米孔的易位过快以至于单个核苷酸的电流电平(current level)过短而难以分辨。纳米孔测序的一种方法涉及在抗电压电位的多核苷酸结合蛋白的引导下控制多核苷酸通过纳米孔的易位,所述结合蛋白如解旋酶、易位酶或聚合酶。尽管有这种控制易位,但是一些测序错误形式仍存在并且导致较差的测序准确度。
因此,仍需要提供进一步控制多核苷酸通过纳米孔的易位并且在核苷酸差异中更好地分辨核苷酸易位的方法和组合物。本发明公开满足了该需要并提供了相关优势。
发明概述
提供了表征目标多核苷酸的方法。所述方法包括:(a)穿过与Hel308解旋酶(解旋酶)和目标多核苷酸接触的孔施加势差;(b)测量目标多核苷酸通过孔的一个或多个分步易位步骤产生的一个或多个信号;和(c)根据分步易位步骤的电信号,表征目标多核苷酸。目标多核苷酸的表征包括鉴别以下中的一个或多个:(1)目标多核苷酸的序列;(2)目标多核苷酸的修饰;(3)目标多核苷酸的长度;(4)目标多核苷酸的身份;(5)目标多核苷酸的来源;或(6)目标多核苷酸的二级结构。还提供了调节目标多核苷酸通过孔的分步易位步骤的方法和表征目标多核苷酸的组合物,其包含孔、Hel308解旋酶和包含在小于1mM ATP溶液或核苷酸类似物溶液中的目标多核苷酸。
附图说明
图1A示出了通过解旋酶的多核苷酸易位的静电尺蠖模型(electrostaticinchworm model)。
图1B示意性示出了根据一些实施方式的第一示例性组合物,其包含与Hel308解旋酶接触的孔。
图1C示意性示出了根据一些实施方式的用于表征目标多核苷酸的示例性方法中的步骤。
图2A示出了根据一些实施方式的Phi29聚合酶和Hel308 Tga解旋酶易位事件的比较。与使用Phi29 DNA聚合酶所观察的易位步骤相比,示出了使用Hel308 Tga解旋酶所观察的分步易位步骤。
图2B示出了根据一些实施方式的Phi29聚合酶和Hel308 Tga解旋酶易位事件的比较。将使用Phi29聚合酶作为分子马达易位通过MspA-M2纳米孔的单链多核苷酸模板所产生的预期电流电平与使用Hel308 Tga解旋酶作为分子马达所观察到的那些相比,示出了用Hel308 Tga解旋酶观察到的分步易位步骤。
图2C示出了根据一些实施方式的Phi29聚合酶和Hel308 Tga解旋酶易位事件的比较。对于简单重复的核苷酸序列(SEQ ID NO:74),与使用Phi29 DNA聚合酶所观察的易位步骤相比,示出了使用Hel308 Tga解旋酶所观察的分步易位步骤。
图3示出了根据一些实施方式对分步易位步骤所提议的“基于紧握”的机制。
图4A和图4B示出了根据一些实施方式ATP浓度对分步易位步骤停留时间的示例性影响。
图5根据一些实施方式,对完整步骤(菱形)和1/2步骤(正方形)绘制了电流迹线电子产生的序列重构准确度(隐马尔可夫模型(HMM))对不同电平的附加噪声的图。
图6A示出了根据一些实施方式HMM解码纳米孔中序列所需的具有非零概率的状态转移,其中多核苷酸是通过马达酶移动的。所述马达为在1个核苷酸的步骤中移动多核苷酸的Phi29 DNAP或类似的酶。
图6B示出了根据一些实施方式HMM解码纳米孔中序列所需的具有非零概率的状态转移,其中多核苷酸是通过马达酶移动的。所述马达为使聚合物能够分步移动的Hel308解旋酶或类似的酶。
图7示出了根据一些实施方式发现电流形状作为高斯偏移的函数的预期准确度。菱形示出了具有完整核苷酸步骤的马达。圆形示出了具有分步易位步骤的马达,而正方形示出了具有与持久时间值合并的分步易位步骤的马达。
图8示出了根据一些实施方式随焦磷酸盐浓度的改变,Hel308解旋酶活性的示例性调节。
图9示出了根据一些实施方式用核苷酸抑制剂原钒酸钠和用核苷酸类似物腺苷5′-(β,γ-亚氨基)三磷酸锂盐水合物对Hel308解旋酶活性的示例性调节。
图10示出了根据一些实施方式使用可以通过电平(level)和电平持续时间(levelduration)从两个独立序列读取获得的额外分步易位步骤所提供的信息的方法的实例。
图11示出了根据一些实施方式使用可以通过电平和电平持续时间从两个并行序列读取获得的额外分步易位步骤所提供的信息的方法的实例。
图12示出了根据一些实施方式使用通过具有或不具有持续时间信息的电流迹线的由额外分步易位步骤所提供的信息的方法的实例。
图13A-图13E示出了根据一些实施方式基于具有Hel308解旋酶3′突出结合位点和胆固醇双层锚点的三元多核苷酸复合物,通过Hel308解旋酶的控制的多核苷酸易位。实心圆形(●)表示5′磷酸盐。实心菱形(◆)表示3′胆固醇。具有凹口的实心半透明圆表示Hel308解旋酶。虚线表示任意长度。大的灰色箭头表示多核苷酸进入或移出孔的多核苷酸移动的方向(通过或抵抗外加场)。大的黑色箭头表示解旋酶沿多核苷酸的易位方向,其为3′至5’。孔(漏斗状圆锥形物体)位于膜(双水平线)中。
图14A-图14D示出了根据一些实施方式基于具有Hel308解旋酶3′突出结合位点和胆固醇双层锚点的三元多核苷酸复合物,通过Hel308解旋酶的控制的多核苷酸易位。实心圆环(●)表示5′磷酸盐。实心菱形(◆)表示3′胆固醇。具有凹口的实心半透明圆表示Hel308解旋酶。虚线表示任意长度。大的灰色箭头表示多核苷酸进入或移出孔的多核苷酸移动的方向(通过或抵抗外加场)。大的黑色箭头表示解旋酶沿多核苷酸的易位方向,其为3′至5’。孔(漏斗状圆锥形物体)位于膜(双水平线)中。符号与图13A-13E中相同。在该图中,存在在多核苷酸“ii”上产生用于Hel308解旋酶结合的3′突出的单个杂交多核苷酸“i”,以及含有任选的胆固醇部分。
图15A-图15C示出了根据一些实施方式的沿梯度力的相同方向的控制的易位。具有凹口的实心半透明圆表示Hel308解旋酶。虚线表示任意长度。大的灰色箭头表示通过外加场多核苷酸移动进入孔的方向。大的黑色箭头表示解旋酶沿多核苷酸的易位方向,其为3′至5’。孔(漏斗状圆锥形物体)位于膜(双水平线)中。
图16示意性示出了已在Hel308作为其成员的SF2家族中鉴别的多个基序(按照出现顺序分别为SEQ ID NO:75-81),例如,DEAD-箱(SEQ ID NO:2)解旋酶(根据Tuteja等人,“Unraveling DNA解旋酶s:Motif,structure,mechanism and function,”EuropeanJournal of Biochemistry 271(10):1849-1863(2004)改编)。
图17A-图17D示出了根据一些实施方式使用某些参数通过Hel308 Mbu解旋酶、Hel308 Tga解旋酶和Phi29聚合酶易位事件产生的示例性信号。
图18示意性示出了根据一些实施方式实施使用表征多核苷酸条型码的分步易位的测定的示例性方法中的步骤。
图19A示意性示出了用于表征来自多核苷酸通过孔的单步易位的信号的示例性隐马尔可夫模型(HMM)的方面。
图19B示意性示出了根据一些实施方式用于表征来自使用Hel308解旋酶的多核苷酸通过孔的分步易位的信号的示例性HMM的方面。
图20A示出了根据一些实施方式使用分步步骤的从头序列的示例性结果。
图20B-图20C示出了根据一些实施方式使用分步步骤的图形匹配的示例性结果。
图21A-21C示意性示出了根据一些实施方式对于多核苷酸通过孔的不同易位,可以作为时间函数产生的信号。
图22A-图22D示出了根据一些实施方式使用多核苷酸通过孔的分布易位所提供的信息的说明性方法中的步骤。
图23示出了根据一些实施方式对于适合用作各自条型码的第一说明性多核苷酸序列(SEQ ID NO:89)和第二说明性多核苷酸序列(SEQ ID NO:90),可以作为时间函数产生的示例性模拟信号。
图24A-图24D示出了根据一些实施方式对于适合用作各自条型码的第一和第二说明性多核苷酸序列可以作为时间函数产生的示例性模拟信号。
图25A和图25B分别示出了根据一些实施方式对于适合用作各自条型码的第一和第二说明性多核苷酸序列可以作为时间函数产生的示例性模拟信号。
图26A-图26D分别示出了根据一些实施方式对于适合用作各自条型码的第一和第二说明性多核苷酸序列作为时间函数产生的示例性测量信号。
发明详述
本发明公开提供了使用目标多核苷酸通过孔易位的一个或多个分步易位步骤来表征目标多核苷酸,包括表征目标多核苷酸的序列的方法和组合物。
在开发纳米孔测序技术时,可以在分子马达,如解旋酶、易位酶或抗势差(例如,抗由势差产生的作用力)的聚合酶的引导下,实现一定水平的多核苷酸通过纳米孔的控制易位。分子马达可以以分步方式移动多核苷酸,通常每步使用一个或多个核苷酸。这种控制的棘轮将多核苷酸通过纳米孔的易位从μsec/核苷酸的天然速率减缓至msec/核苷酸。
分子马达可以使用核苷酸水解的能量来驱动多核苷酸通过纳米孔的易位。解旋酶是其中ATP水解是用于多核苷酸易位的能源的实例。图1中的卡通画示出了通过解旋酶的多核苷酸易位的静电尺蠖模型(electrostatic inchworm model)(参见Frick等人,CurrentPharmaceutical Design,12:1315-1338(2006))。在该模型中,单链多核苷酸保持在带负电荷的槽(cleft)中,其将解旋酶的两个RecA结构域与第三个结构域分开。在不存在ATP的情况下,书夹残基(bookend residue)(例如,HCV解旋酶中的Trp501)和夹钳残基(例如,HCV解旋酶中的Arg393)防止单链多核苷酸滑过槽。一旦ATP结合,RecA结构域旋转,从而移动带正电荷的Arg-夹钳。Arg-夹钳吸引带负电荷的单链多核苷酸,其反过来清除书夹。然后,带负电荷的槽排斥单链多核苷酸,并且单链多核苷酸易位通过解旋酶直至ATP被水解。因此,在该示例性模型中,多核苷酸易位通过解旋酶包括至少两个步骤:其中解旋酶结合至ATP并经历构象变化的第一步骤,和ATP被水解并且多核苷酸易位通过解旋酶的第二步骤。
图1B示意性示出了根据一些实施方式的第一示例性组合物,其包含与Hel308解旋酶接触的孔。在图1B中,具有凹口的实心半透明圆表示如本文所提供的Hel308解旋酶。直线表示多核苷酸,而虚线表示任意长度的多核苷酸。大的灰色箭头表示多核苷酸移动进入或移出孔的方向,而大的黑色箭头表示解旋酶沿多核苷酸易位的方向,其为3′至5′。在所示实施方式中,孔(漏斗状圆锥形物体)位于膜(双水平线)中,尽管可以适合地使用其它孔构型。在图1B所示的实施方式中,多核苷酸运动的方向可以与孔两侧通过势差(说明性地,180V的势差,尽管可以适合地使用其它势差)所产生的外加场相同。为了使多核苷酸移动方向与孔两侧通过势差所产生的外加场相反,因此可以翻转DNA取向,如以下参考图15A-15C更详细地描述的。如本文更详细地提供的,Hel308解旋酶可以导致多核苷酸通过孔的分步易位,其可以有利于核苷酸的表征。例如,基于可以表征的多核苷酸,这种分步易位可以生产一个或多个信号。所述一种或多种信号可以包括如在本文其它处所描述的电信号,或者可以包括如在本文其它处所描述的光信号。示例性电信号可以选自电流、电压、隧穿、电阻、电位、电压、电导率和横向电测量的测量。
说明性地,由于Hel308解旋酶使多核苷酸分步易位通过孔,因此孔内不同核苷酸碱基的通过可以导致通过孔的电流的可测量的变化;这种电流可以被称为“阻塞”电流。如本文更详细地描述的,可以基于信号变化,例如,基于通过孔的电流变化(这些变化基于通过Hel308解旋酶的多核苷酸通过孔的分步易位步骤)确定多核苷酸的一个或多个特征,如多核苷酸序列、多核苷酸修饰、多核苷酸长度、多核苷酸身份、多核苷酸来源或多核苷酸二级结构或其任何适合的组合。在其中孔不对称(例如,包括孔口直径大于孔底)的实施方式中,Hel308解旋酶可以与孔口接触,如图1B中所示。这种构型可以被称为“正向”构型。更一般地说,“正向构型”可以表示不考虑孔是否包括比孔底更宽的孔口,其中分子可以天然地通过孔的方向。作为另外一种选择,可以任意定义“正向方向”。
图1C示意性示出了根据一些实施方式的用于表征目标多核苷酸的示例性方法中的步骤。所述方法可以包括在与Hel308解旋酶和目标多核苷酸接触的孔两侧施加势差的步骤(步骤110)。以与以下参考图13A-13E和14A-14D进一步描述的类似的方式,多核苷酸的易位可以处于施加在易位通过孔的多核苷酸上的由势差所引起的外加力相反的方向,或者多核苷酸的易位可以处于施加在易位通过孔的多核苷酸上的由势差所引起的外加力的方向。任选地,步骤110-130可以重复一次或多次。分步易位步骤(步骤120)可以包括Hel308解旋酶完整易位循环的第一分步易位步骤或可以包括Hel308解旋酶完整易位循环的第二易位步骤。
如本文所使用的,术语“多核苷酸”是指脱氧核糖核酸(DNA)、核糖核酸(RNA)或其类似物。多核苷酸可以是单链的、双链的或含有单链和双链序列两者。多核苷酸分子可以来源于双链DNA(dsDNA)形式(例如,基因组DNA、PCR和扩增产物等),或者可以来源于如DNA(ssDNA)或RNA的单链形式并且其可以转化为dsDNA形式,并且反之亦然。多核苷酸分子的准确序列可以是已知的或未知的。以下是多核苷酸的示例性实例:基因或基因片段(例如,探针、引物、EST或SAGE标签)、基因组DNA、基因组DNA片段、外显子、内含子、信使RNA(mRNA)、转运RNA、核糖体RNA、核糖酶、cDNA、重组多核苷酸、合成多核苷酸、分枝多核苷酸、质粒、载体、任何序列的分离的DNA、任何序列的分离的RNA、任何上述序列的核酸探针、引物或扩增拷贝。
多核苷酸可以包括核苷酸或核苷酸类似物。核苷酸通常含有糖、核碱基和至少一个磷酸基。核苷酸可以是无碱基的(abasic)(即,缺少核碱基)。核苷酸包括脱氧核糖核苷酸、修饰的脱氧核糖核苷酸、核糖核苷酸、修饰的核糖核苷酸、肽核苷酸、修饰的肽核苷酸、修饰磷酸盐糖主链核苷酸及其混合物。核苷酸的实例包括(例如)腺苷一磷酸(AMP)、腺苷二磷酸(ADP)、腺苷三磷酸(ATP)、胸苷一磷酸(TMP)、胸苷二磷酸(TDP)、胸苷三磷酸(TTP)、胞苷酸(CMP)、胞苷二磷酸(CDP)、胞苷三磷酸(CTP)、鸟苷一磷酸(GMP)、鸟苷二磷酸(GDP)、鸟苷三磷酸(GTP)、尿苷一磷酸(UMP)、尿苷二磷酸(UDP)、尿苷三磷酸(UTP)、脱氧腺苷酸(dAMP)、脱氧腺苷二磷酸(dADP)、脱氧腺苷三磷酸(dATP)、脱氧胸腺嘧啶核苷一磷酸(dTMP)、脱氧胸腺嘧啶核苷二磷酸(dTDP)、脱氧胸苷三磷酸(dTTP)、去氧胞二磷(dCDP)、脱氧胞苷三磷酸(dCTP)、脱氧鸟苷一磷酸(dGMP)、脱氧鸟苷二磷酸(dGDP)、脱氧鸟苷三磷酸(dGTP)、脱氧尿苷一磷酸(dUMP)、脱氧尿苷二磷酸(dUDP)和脱氧尿苷三磷酸(dUTP)。还可以在本文所述的方法中使用包含修饰的核碱基的核苷酸类似物。无论是具有天然主链还是类似结构,可以包含在多核苷酸中的示例性修饰的核碱基包括(例如)肌苷、黄嘌呤(xathanine)、次黄嘌呤(hypoxathanine)、异胞嘧啶、异鸟嘌呤、2-氨基嘌呤、5-甲基胞嘧啶、5-羟甲基胞嘧啶、2-氨基腺嘌呤、6-甲基腺嘌呤、6-甲基鸟嘌呤、2-丙基鸟嘌呤、2-丙基腺嘌呤、2-硫脲嘧啶、2-硫胸腺嘧啶、2-硫胞嘧啶、15-卤代脲嘧啶、15-卤代胞嘧啶、5-丙炔基尿嘧啶、5-丙炔基胞嘧啶、6-偶氮尿嘧啶、6-偶氮胞嘧啶、6-偶氮胸腺嘧啶、5-尿嘧啶、4-硫尿嘧啶、8-卤代腺嘌呤或鸟嘌呤、8-氨基腺嘌呤或鸟嘌呤、8-硫腺嘌呤或鸟嘌呤、8-硫烷基腺嘌呤或鸟嘌呤、8-羟基腺嘌呤或鸟嘌呤、5-卤素取代的尿嘧啶或胞嘧啶、7-甲基鸟嘌呤、7-甲基腺嘌呤、8-氮杂鸟嘌呤、8-氮杂腺嘌呤、7-去氮鸟嘌呤、7-去氮腺嘌呤、3-去氮鸟嘌呤、3-去氮腺嘌呤等。如本领域中已知的,某些核苷酸类似物不能引入多核苷酸,例如,核苷酸类似物,如腺苷5′-磷酰硫酸。
如本文所使用的,术语“孔”旨在表示延伸穿过屏障(如膜)的结构,其允许离子和/或水溶性分子从屏障一侧穿过至屏障另一侧。孔可以(但不必)存在于膜中。例如,通常抑制离子或水溶性分子通过的屏障可以包括延伸穿过屏障以允许离子或水溶性分子从屏障一侧穿过至屏障另一侧的多孔结构。孔(例如,跨膜孔)包括(例如)生物孔、固态孔和生物和固态杂交孔。
如本文所使用的,术语“生物孔”旨在表示由生物来源材料制成的孔,其延伸穿过屏障(包括,例如,膜),其允许离子和/或水溶性分子从屏障一侧穿过至屏障另一侧。生物来源是指来源于或分离自生物环境(如生物或细胞)的材料,或生物学可用结构的合成产生形式。生物孔包括(例如)多肽孔和多核苷酸孔。
如本文所使用的,术语“多肽孔”旨在表示延伸穿过屏障(如,例如,膜)并且允许离子和/或水溶性分子从屏障一侧流至屏障另一侧的一个或多个多肽。多肽孔可以是单体、均聚物或杂聚物。多肽孔的结构包括(例如)α螺旋束孔和β-折叠孔以及在本领域中熟知的所有其它结构。示例性多肽孔包括α-溶血素、耻垢分枝杆菌(Mycobacterium smegmatis)孔蛋白A、短杆菌肽A,麦芽糖孔蛋白、OmpF、OmpC、PhoE、Tsx、F-pilus、SP1(Wang等人,Chem.Commun.,49:1741-1743,2013)和线粒体孔蛋白(VDAC)XX、Tom40(美国专利No.6,015,714和Derrington等人,Proc.Natl.Acad.Sci.USA,107:16060(2010))。“耻垢分枝杆菌(Mycobacterium smegmatis)孔蛋白A(MspA)”是分枝杆菌(Mycobacteria)产生的膜孔蛋白,其允许亲水性分子进入细菌。MspA形成紧密互相连接的八聚物和类似于高脚杯的跨膜β-折叠,并且含有中性通道/孔。
如本文所使用的,术语“多核苷酸孔”旨在表示延伸穿过屏障(如,例如,膜)并且允许离子和/或水溶性分子从屏障一侧流至屏障另一侧的一个或多个多核苷酸。多核苷酸孔可以包括(例如)多核苷酸折纸。
如本文所使用的,术语“固态孔”旨在表示由非生物来源材料制成的孔,其延伸穿过屏障(如,例如,膜),并且允许离子和/或水溶性分子从屏障一侧穿过至屏障另一侧。固态旨在表示非生物来源的材料。固态孔可以是无机或有机材料。固态孔包括(例如)氮化硅孔、二氧化硅孔和石墨烯孔。
如本文所使用的,术语“生物和固态杂交孔”旨在表示由生物和非生物来源的材料两者制成的杂交孔,其延伸穿过屏障(如,例如,膜),并且允许水合离子和/或水溶性分子从屏障一侧穿过至屏障另一侧。生物来源的材料如上所定义并且包括(例如)多肽和多核苷酸。生物和固态杂交孔包括(例如)多肽-固态杂交孔和多核苷酸-固态孔。
如本文所使用的,术语“解旋酶”旨在表示具有使用来源于(例如)核苷三磷酸(NTP)水解的能量使双链多核苷酸打开的活性的多核苷酸结合蛋白。打开双链多核苷酸导致多核苷酸沿其活性位点易位。该术语旨在包含具有使单链多核苷酸以及部分双链多核苷酸易位或与之结合的活性的多肽。“Hel308解旋酶”是ATP依赖性DNA解旋酶和超家族2解旋酶。来自果蝇(Drosophila melanogaster)的组成成员Mus308由融合至C末端DNA聚合酶结构域的N末端SF2解旋酶结构域组成。智人中的Hel308作为具有有限前行能力的3′至5′DNA解旋酶SF2起作用。Hel308解旋酶与ski2样解旋酶可互换使用。有用的同源物可以仅包含解旋酶结构域(即缺少聚合酶结构域)。仅有解旋酶的同源物存在于后生动物和古细菌中。后生动物的实例是人Hel308和Mus301。古细菌的实例是Tga和Mbu。
除非本文明确描述,否则如本文所使用的术语“Hel308解旋酶底物”旨在表示核苷酸或核苷酸类似物,其能够被解旋酶水解并提供使双链或部分双链多核苷酸打开或使单链多核苷酸易位的能量。Hel308解旋酶的常见底物包括ATP。然而,该术语含义内的其它Hel308解旋酶底物包括ATP以外的核苷酸,如先前描述的那些,和能够被Hel308解旋酶水解的核苷酸类似物。示例性类似物包括(例如)磷酸盐类似物,如γ巯基类似物、α巯基类似物等、ATPγS、ATPαS、AMP、PNP、ApCpp、AppCp和AppNHp。
如本文所使用的,术语“易位”旨在表示目标多核苷酸沿解旋酶和/或孔(或在它们内部)的运动。
如本文所使用的,当结合解旋酶使用时,术语“完整易位循环”旨在表示目标多核苷酸的一个或多个核苷酸单元沿解旋酶和/或孔运动的完整间隔。完全间隔可以在循环中的任一点开始,并且可以(例如)包括图3中所示的间隔,其包括ATP结合以及所结合的ATP水解的步骤。因此,如本文所使用的完整易位循环可以在核苷酸底物结合开始并在核苷酸底物水解结束。类似地,完整易位循环可以在核苷酸底物水解开始并在核苷酸结合结束。类似地,完整易位循环可以在以上举例说明的两个起始点之间的任一点开始,只要它在起始点前的步骤结束。
如本文所使用的,当结合解旋酶使用时,术语“分步易位步骤”旨在表示表征完整易位循环的一部分的可检测事件。例如,分步易位步骤可以是目标多核苷酸的一个或多个核苷酸单元沿解旋酶和/或孔的部分易位。在具体的实施方式中,当发生构象变化时,可以在ATP结合和水解之间发生分步步骤。构象变化有效地将完整易位循环划分为至少两个部分或分步易位步骤。分步步骤可以或可以不伴有沿解旋酶的核酸运动。
如本文所使用的,术语“信号”旨在表示提供信息的指示物。信号包括(例如)电信号和光信号。
如本文所使用的,术语“电信号”旨在表示提供信息的电学特性的指示物。所述指示物可以是(例如)电流、电压、隧穿、电阻、电位、电压、电导率和横向电测量。“电流”是指电荷流。当穿过孔施加势差时电荷流动。
如本文所使用的,术语“光信号”旨在表示提供信息的光学性能的指示物。光信号包括(例如)荧光信号和拉曼信号。
如本文所使用的,术语“同源性”旨在表示两个多核苷酸之间或两个多肽之间的序列相似性。可以通过比较每条序列中可以出于比较的目的进行比对的位置来确定相似性。序列间的相似度是所述序列共有的匹配或同源位置的数目的函数。可以使用本领域中已知的软件程序进行两序列的比对以确定它们的序列相似性百分比,所述软件程序如(例如)Ausubel等人,Current Protocols in Molecular Biology,John Wiley and Sons,Baltimore,MD(1999)中所述的那些。优选地,使用缺省参数进行比对,它的实例如下所述。可以使用的本领域熟知的一个校直程序是设置为缺省参数的BLAST。具体地,程序为BLASTN和BLASTP,使用以下缺省参数:遗传密码=标准;过滤器=无;链=两条;截止值=60;预期值=10;基质=BLOSUM62;描述=50条序列;分类=高分;数据库=非冗余,GenBank+EMBL+DDBJ+PDB+GenBank CDS translations+SwissProtein+SPupdate+PIR。这些程序的详细内容可见于国家生物技术信息中心(National Center for Biotechnology Information)。
本发明公开提供了表征目标多核苷酸的方法。所述方法包括:(a)横跨与Hel308解旋酶和目标多核苷酸接触的孔施加势差;(b)测量目标多核苷酸通过孔的一个或多个分步易位步骤产生的一个或多个信号;和(c)根据分步易位步骤的电信号,表征目标多核苷酸。
如本文所述,多核苷酸包括脱氧核糖核酸(DNA)、核糖核酸(RNA)或其类似物。多核苷酸通常将含有磷酸二酯键,尽管在一些情况下,多核苷酸还可以具有交替的主链,其包括(例如)磷酰胺(Beaucage等人,Tetrahedron,49(10):1925(1993)及其中的参考文献;Letsinger,J.Org.Chem.,35:3800(1970);Sprinzl等人,Eur.J.Biochem.,81:579(1977);Letsinger等人,Nucl.Acids Res.,14:3487(1986);Sawai等人,Chem.Lett.,805(1984);Letsinger等人,J.Am.Chem.Soc.,110:4470(1988);和Pauwels等人,Chemica Scripta,26:141(1986))、硫代磷酯(Mag等人,Nucleic Acids Res.,19:1437(1991);和美国专利No.5,644,048)、二硫代磷酸酯(Briu等人,J.Am.Chem.Soc.,111:2321(1989))、O-甲基亚磷酰胺键(参见Eckstein,Oligonucleotides and Analogues:A Practical Approach,OxfordUniversity Press)以及肽核酸主链和键(参见Egholm,J.Am.Chem.Soc.,114:1895(1992);Meier等人,Chem.Int.Ed.Engl.,31:1008(1992);Nielsen,Nature,365:566(1993);Carlsson等人,Nature,380:207(1996))。其它多核苷酸包括具有带正电荷的主链(Denpcy等人,Proc.Natl.Acad.Sci.USA,92:6097(1995));非离子主链(美国专利No.5,386,023、5,637,684、5,602,240、5,216,141和4,469,863;Kiedrowshi等人,Angew.Chem.Int.Ed.English,30:423(1991);Letsinger等人,J.Am.Chem.Soc.,110:4470(1988);Letsinger等人,Nucleosides&Nucleotides,13:1597(1994);第2和3章,ASCSymposium Series 580,“Carbohydrate Modificatiohs in Antisense Research”,Y.S.Sanghui和P.Dan Cook主编;Mesmaeker等人,Bioorganic&Medicinal Chem.Lett.,4:395(1994);Jeffs等人,J.Biomolecular NMR,34:17(1994);Tetrahedron Lett.,37:743(1996))和非核糖主链的那些,所述非核糖主链包括美国专利No.5,235,033和5,034,506,以及第6和7章,ASC Symposium Series 580,“Carbohydrate Modifications inAntisense Research”,Y.S.Sanghui主编P.Dan Cook中所述的那些。含有一个或多个碳环糖的多核苷酸分子也包含在多核苷酸的定义内(参见Jenkins等人,Chem.Soc.Rev.,(1995)第169-176页)。在Rawls,C&E News,Jun.2,1997,第35页中描述了一些多核苷酸。
可以根据本发明所述的方法表征目标多核苷酸。示例性多核苷酸包括(例如)基因或基因片段(例如,探针、引物、EST或SAGE标签)、基因组DNA、基因组DNA片段、外显子、内含子、信使RNA(mRNA)、转运RNA、核糖体RNA、核糖酶、cDNA、重组多核苷酸、合成多核苷酸、分枝多核苷酸、质粒、载体、任何序列的分离的DNA、任何序列的分离的RNA、任何上述序列的核酸探针、引物或扩展拷贝。
在本文的具体实施方式中使用的目标多核苷酸可以具有任意多种长度,通常具有足以延伸通过孔并在孔的一侧被解旋酶结合的长度。一般说来,这种长度为至少约10个核苷酸长。然而,对于使用本发明所述的方法的表征来说,比该最小尺寸长的多个长度是适用的。有用的多核苷酸的示例性长度包括(例如)至少约10,20,30,40,50,60,70,80,90,100,150,200,300,400,500,1,000,5,000或10,000、100,000个核苷酸或更长。作为另外一种选择或者另外,所述长度可以不长于1,000,000、100,000、10,000、1,000、100个核苷酸或更少。因此,可以使用本发明所述的方法测序的多核苷酸可以(例如)在短多核苷酸、片段、cDNA、基因和基因组片段的范围内。
在本发明所述的方法中使用的多核苷酸可以是单链的、双链的或含有单链和双链序列两者。多核苷酸分子可以来自于双链多核苷酸(例如,dsDNA)并且可以转化为单链多核苷酸。多核苷酸分子还可以来自于单链多核苷酸(例如,ssDNA、ssRNA),并且所述ssDNA可以转化为双链多核苷酸。在本发明公开的一些方面,双链或部分双链多核苷酸包括封闭性多核苷酸。这类多核苷酸可以包括结合本文中图13A-13E、14A-14D和15A-15C举例说明的那些。在WO 2013/057495中描述了多核苷酸易位通过孔的示例性形式。
在一些方面,本发明公开提供了表征目标多核苷酸的方法。所述方法包括鉴别:(1)目标多核苷酸的序列;(2)目标多核苷酸的修饰;(3)目标多核苷酸的长度;(4)目标多核苷酸的身份;(5)目标多核苷酸的来源;或(6)目标多核苷酸的二级结构。
多核苷酸序列是指多核苷酸的一级结构或者多核苷酸分子中核苷酸的相继次序。通过使用目标多核苷酸通过孔的分步易位步骤所产生的信号来表征目标多核苷酸中的核苷酸,可以确定多核苷酸的序列。
多核苷酸的修饰是指多核苷酸中核苷酸的任何共价或非共价修饰,包括(例如)核苷酸甲基化或羟甲基化。的确,修饰可以包括可以引入多核苷酸单链中的任意数目的核苷酸类似物,包括,例如,8-氧代鸟苷、5-甲酰胞嘧啶和5-羧基胞嘧啶以及在本文其它处所描述的其它类似物。核苷酸的修饰提供了信号中相应的变化。因此,通过使用目标多核苷酸通过孔的分步易位步骤所产生的信号来表征目标多核苷酸中修饰的核苷酸,可以确定多核苷酸的一个或多个修饰。
多核苷酸的长度是指多核苷酸中核苷酸的数目。多核苷酸的长度可以通过(例如)确定多核苷酸的一级序列或通过测量其在孔中的停留时间或者通过对穿过孔的核苷酸数目计数来确定。在一些实施方式中,停留时间对应于电流瞬时变化的持续时间。瞬时变化可以认为是由于多核苷酸的存在所造成的孔电流的任何偏离。在一些实施方式中,所述偏离导致电流幅度减小。这种减小通常可以是原始未封闭的孔电流的至多95%,90%,80%,60%,50%,40%,30%,20%或10%或更小。作为另外一种选择或者另外,所述减小可以是至少10%,20%,30%,40%,50%,60%,70%,80%或90%或以上。在一些情况下,多核苷酸可以导致电流幅度相对于未封闭孔提高。可以通过取决于所使用的实验条件的可重复数学函数来描述持续时间和多核苷酸长度之间的关系。对于给定类型的多核苷酸(例如,DNA或RNA),所述函数可以是线性或非线性的(例如,S形或指数)函数。
多核苷酸的身份是指多核苷酸的类型。身份还可以是指如本领域中已知的多核苷酸的名称。例如,多核苷酸的身份可以是(例如)DNA、RNA、双链多核苷酸、单链多核苷酸和/或部分双链多核苷酸。多核苷酸的身份还可以包括确定多核苷酸的基因产物或结构功能。例如,多核苷酸可以编码多肽或者它可以是结构多核苷酸,如核糖体RNA。可以根据所有或部分多核苷酸的核苷酸序列、与所有或部分多核苷酸互补的第二多核苷酸的序列、所有或部分多核苷酸所编码的RNA序列或者所有或部分多核苷酸所编码的蛋白质序列来确定多核苷酸的身份。在具体的实例中,可以通过构成多核苷酸一部分的“标签”或“条形码”序列鉴别多核苷酸。在这些实例中,可以通过对标签或条型码所预期的信号类型来分配多核苷酸的身份。多核苷酸的来源可以是指多核苷酸的物种来源或是指合成来源。可以使用本领域中熟知的程序,例如,BLASTN,通过在多核苷酸序列数据库中对多核苷酸序列进行比对来确定多核苷酸的身份和来源。
多核苷酸的二级结构是指多核苷酸分子中自互补性区域的分子内碱基配对。示例性二级结构包括(例如)双螺旋、发夹、环、凸起(bulges)、双螺旋、接合(junction)、茎、假结(pseudoknot)、三螺旋、H-DNA、锤头型(hammerhead)和自拼接核糖酶。可以(例如)通过测量其在孔中停留时间的相应变化或测量通过分步易位步骤产生的信号的相应变化来确定多核苷酸的二级结构。
孔是延伸穿过屏障(包括,例如,膜)的结构,其允许离子和/或水溶性分子从屏障一侧穿过至屏障另一侧。孔可以(但不必)存在于膜中。例如,通常抑制离子或水溶性分子通过的屏障可以包括延伸穿过屏障以允许离子或水溶性分子从屏障一侧穿过至屏障另一侧的多孔结构。本发明公开的膜可以是(例如)不透性或半透性屏障,其分成可以具有相同或不同的组成的两个液体室。可以根据本发明公开使用任何膜,只要所述膜可以配置以包含跨膜的孔并维持膜两侧的势差。适合的势差如下所述。
可以在本发明公开的组合物和方法中使用本领域中熟知的多种膜。在本领域中熟知的这些膜包括多种不同的结构和组成。例如,膜可以是单层或多层结构,只要可以引入用于多核苷酸表征的孔。膜中的层是指形成屏障的不透性或半透性材料。以下还描述了单层和多层膜的实例。
形成膜的材料可以是生物或非生物来源的。生物来源的材料是指来源于或分离自生物环境(如生物或细胞)的材料,或生物学可用结构的合成产生形式。由生物来源材料制成的示例性膜包括脂质双分子层。非生物来源的材料也称为固态材料并且可以形成固态膜。
适合的脂质双分子层和制备或获得脂质双分子层的方法在本领域中是熟知的并且在(例如)美国专利公开US 2010/0196203和PCT专利公开WO2006/100484中公开。适合的脂质双分子层包括(例如)细胞膜、细胞器膜、脂质体、平面脂质双分子层和支撑的脂质双分子层。可以(例如)从两个相对的磷脂层形成脂质双分子层,所述磷脂的布置使得它们疏水性尾部基团彼此面对以形成疏水性内部,而脂肪的亲水性头部基团朝外面对双分子层每侧的水相环境。还可以(例如)通过Montal and Mueller的方法(Proc.Natl.Acad.Sci.USA.,1972;69:3561-3566)形成脂质双分子层,其中脂质单分子层从水溶液/空气界面上通过垂直于界面的孔的任一侧。通常通过首先将其溶于有机溶剂,然后使一滴溶剂在孔任一侧上的水溶液表面上蒸发,将脂质加入到电解质水溶液的表面。一旦有机溶剂蒸发,孔任一侧上的溶液/空气界面物理地上下移动通过孔直至形成双层。双分子层形成的其它常规方法包括尖端浸渍(tip-dipping)、涂覆双分子层(painting bilayers)和脂质体双分子层膜片钳(patch-clamping)。获得或产生脂质双分子层的多种其它方法在本领域中是熟知的并且等同地适用于本发明公开的组合物和方法中的用途。
固态膜在本领域中是熟知的并且在(例如)PCT专利公开WO2000/079257中公开。如上所述,固态膜由非生物来源的材料的一个或多个层制成。固态膜可以是单层,如支持基底上的涂层或膜,或者是自立式成分(free-standing element)。固态膜还可以是夹心构造的多层复合材料。对根据本发明公开可以使用的材料没有具体限制,只要可以配制所得固态膜以包含跨膜孔并且在膜两侧建立势差。可以从有机和无机材料制备固态膜,包括(例如)微电子材料、绝缘材料如Si3N4、Al2O3和SiO、有机和无机聚合物如聚酰胺、三嵌段共聚物(例如,两亲PMOXA-PDMS-PMOXA ABA三嵌段共聚物)、塑料材料如或弹性体如二组分加成-固化硅橡胶和玻璃。另外,可以从石墨烯单层(石墨烯是密集挤压成二维蜂房晶格的原子薄层)、石墨烯多层或与一个或多个其它固态材料层混合的一个或多个石墨烯层制备固态膜(PCT专利公开WO 2013/016486)。含有石墨烯的固态膜可以包括石墨烯纳米带或石墨烯纳米间隙中的至少一种石墨烯层,其可以用作电学传感器来表征目标多核苷酸(参见PCT专利公开WO2013/016486)。可以通过本领域中熟知的方法制备固态膜。例如,可以通过化学气相淀积(CVD)或从石墨剥落(PCT专利公开WO 2013/016486)制备石墨烯膜。
本发明公开的组合物和方法可以使用位于屏障中的孔以用于目标多核苷酸的表征。可以从生物或非生物来源的材料制备孔。因此,孔包括(例如)生物孔、固态孔和生物和固态杂交孔。
孔可以具有有利于多核苷酸中核苷酸序列检测的与其相关的功能性。例如,孔可以包含酶(如解旋酶)或者其它功能性,所述酶或其它功能性连接至孔、与孔结合或位于孔附近以控制多核苷酸通过孔的速率。孔可以具有与之结合的检测电路或传感器,包括(例如)膜片钳电路、隧道电极电路或横向电导测定电路(如石墨烯纳米带或石墨烯纳米间隙)。孔还可以包括检测多核苷酸上标记物(包括(例如)荧光部分或拉曼信号产生部分)的光学传感器,其基于片段与孔(例如,将片段通过孔)的相互作用确定核苷酸序列。
在具体的实施方式中,可以在本发明公开的组合物和方法中使用生物孔(包括多肽孔和多核苷酸孔),只要所述孔具有允许多核苷酸通过屏障(例如,膜)的缢缩区。缢缩区是孔腔中的位置,在此分析物(例如,多核苷酸或核苷酸)的阻塞影响了孔所产生的可检测信号。可以在本发明公开的组合物和方法中使用具有多个缢缩区长度的孔,其包括(例如)至少1、2、3、4、5、6、7、8、9或10个核苷酸的长度。作为另外一种选择或者另外,可以使用至多约10、9、8、7、6、5、4、3、2或1个核苷酸长度。然而,缢缩区长度可以影响信号质量。例如,较短的缢缩区可以导致更好的核苷酸易位分辨率或重构准确度。在一个实施方式中,生物孔具有约5个核苷酸或以下的缢缩区,位于缢缩区的5个或小于5个核苷酸调节电信号,其具有比得自大于5个核苷酸的电信号更好的核苷酸易位分辨率。在一些情况下,对于小于2nt的缢缩,信噪比提高不会导致测序准确度改善。当由于酶的随机运动而跳过核苷酸时,如果大于较小缢缩的均聚物不再能被检测并且再读取的缺少降低准确度时,这可以产生。因此,根据本发明公开,可以使用适合的多肽孔和具有5个核苷酸或以下的缢缩区的多核苷酸孔。考虑到本文所提供的教导和指导,本领域技术人员将理解对于具体需要,什么长度的缢缩区是适用的。例如,在需要高质量结果的应用中,本领域技术人员可以使用具有较短缢缩区的孔。
生物孔是由生物来源材料制成,延伸穿过屏障(例如,膜)的孔,其允许离子和/或水溶性分子从屏障一侧穿过至屏障另一侧。如本文所述使用的膜,当提及孔时,生物来源是指来源于或分离自生物环境(如生物或细胞)的结构,或者生物学可用结构的合成产生形式。生物来源材料包括(例如)多肽和多核苷酸。因此,生物孔包括(例如)多肽孔和多核苷酸孔。
重构成屏障(例如,膜),如脂质双分子层的多肽孔可以用于纳米孔测序。存在可以根据本发明公开使用的多种多肽孔,只要多肽可以形成允许目标多核苷酸穿过屏障(例如,膜)的缢缩区。根据所涉及的多肽,多肽孔可以是单体、均聚物或杂聚物。多肽孔可以包括一些重复亚基,如7或8个亚基。因此,所述多肽孔可以是(例如)六聚、七聚或八聚孔。
多肽孔包括(例如)α-螺旋束孔和β-折叠孔以及在本领域中熟知的所有其它孔。α-螺旋束孔包括由α-螺旋形成的孔。适合的α-螺旋束孔包括(例如)内膜蛋白和α外膜蛋白,如WZA和ClyA毒素。β-折叠孔包括通过β-链形成的孔。适合的β-折叠孔包括(例如)β-毒素,如α-溶血素、炭疽毒素和杀白细胞素,和细菌的外膜蛋白/孔蛋白,如耻垢分枝杆菌(Mycobacterium smegmatis)孔蛋白(Msp),包括MspA、外膜孔蛋白F(OmpF)、外膜孔蛋白G(OmpG)、外膜磷脂酶A和奈氏球菌(Neisseria)自转运蛋白脂蛋白(NalP)。其它孔包括(例如)胞溶素(参见,例如,WO 2013 153359)或来自皮疽诺卡氏菌(Norcadia farcinica)的MspA同源物。
α-溶血素多肽是可以在本发明公开的方法和组合物中使用的七聚多肽孔。它由连接到5nm长的β-折叠的3.6nm孔腔(vestibule)组成,其含有允许单链多核苷酸通过但不允许双链多核苷酸通过的1.4nm缢缩。α-溶血素的约5nm长的圆柱形β-折叠孔每次可以容纳多至约10个核苷酸。位于该β-折叠中的核苷酸显著调节孔电流并随后削弱对最窄1.4nm的孔缢缩中的单个核苷酸特异的离子特征,从而降低测序应用中核苷酸易位的整体分辨率。
MspA是可以在本发明公开的组合物和方法中使用的八聚多肽孔。它含有直径~1.2nm,缢缩长度~0.5nm的单个缢缩;与α-溶血素的圆柱形结构相反,内孔形成锥形漏斗形状。Derrington等人表明基因工程MspA区分三核苷酸组(AAA、GGG、TTT、CCC)的能力,其核苷酸分离效率与天然α-溶血素相比具有令人深刻印象的3.5倍提高(Derrington等人,Proc.Natl.Acad.Sci.USA,107:16060(2010))。据报导在涉及固定化单链多核苷酸的实验中,观察到MspA缢缩内或附近少至3个核苷酸有助于孔电流,比已知调节天然α-溶血素中离子电流的~10个核苷酸具有显著改善。作者假设这可以通过定点突变进一步改善至或许单个核苷酸,这是未来MspA突变的目标。
在一些方面,多肽孔为耻垢分枝杆菌(Mycobacterium smegmatis)孔蛋白A(MspA)。在一些方面,所述MspA具有SEQ ID NO:1所示或与SEQ ID NO:1具有至少15%、至少20%、至少25%、至少30%、至少35%、至少40%、至少45%、至少50%、至少55%、至少60%、至少65%、至少70%、至少75%、至少80%、至少85%、至少90%、至少95%或至少99%同源性的氨基酸序列。
MspA是适合的多肽孔。另外,MspA突变体可以在本发明公开的组合物和方法中使用以调节多核苷酸通过孔的易位。在本文的实施方式中使用的MspA孔可以具有SEQ ID NO:1所示的氨基酸序列,其对应于GLDNELSLVDGQDRTLTVQQWDTFLNGVFPLDRNRLTREWFHSGRAKYIVAGPGADEFEGTLELGYQIGFPWSLGVGINFSYTTPNILINNGNITAPPFGLNSVITPNLFPGVSISARLGNGPGIQEVATFSVRVSGAKGGVAVSNAHGTVTGAAGGVLLRPFARLIASTGDSVTTYGEPWNMN,它是具有以下突变的MspA序列:D90N、D91N、D93N、D118R、D134R&E139K。SEQ ID NO:1所示的MspA孔突变体被称为“M2NNN”。可以在本发明公开的组合物和方法中使用其它MspA突变体,其与SEQ ID NO:1具有至少15%、至少20%、至少25%、至少30%、至少35%、至少40%、至少45%、至少50%、至少55%、至少60%、至少65%或至少70%、至少75%、至少80%、至少85%、至少90%、至少95%或至少99%同源性。对另一序列具有特定百分比(例如,50%)的同源性的多肽或多肽区(或者多核苷酸或多核苷酸区)是指当比对时,相比较的两个序列中氨基酸(或核苷酸碱基)的百分比是相同的。如本文所述,可以使用本领域中已知的软件程序来进行确定它们序列同一性百分比的两个序列的比对。可以根据本领域熟知的方法对天然MspA多肽进行突变,包括特定区域或特定氨基酸残基的***、缺失、取代或其它所选修饰,所述突变方法包括编码MspA多肽的核酸的定点突变(Zoller,M.J.,Curr.Opin.Biotechnol.,3:348-354,(1992))。在US2012/0055792A1中还描述了有用的MspA突变体。
可以通过多种本领域熟知的方法分离在本发明公开的组合物和方法中使用的天然或突变MspA多肽,例如,重组表达***、沉淀、凝胶过滤、离子交换、反相和亲和色谱法等。其它熟知的方法描述于Deutscher等人,Guide to Protein Purification:Methods inEnzymology,Vol.182,(Academic Press,(1990))。作为另外一种选择,可以使用熟知的重组方法获得本发明公开的分离的天然或突变MspA多肽。本领域技术人员可以选择本发明公开的天然或突变MspA多肽的生物纯化的方法和条件,并且可以(例如)通过功能测定来监控纯化。
制备天然或突变MspA多肽的一个示例性方法是使用本领域中熟知的方法在适合的宿主细胞,如细菌细胞、酵母细胞或其它适合的细胞中表达编码所述多肽的多核苷酸,并回收所表达的天然或突变体MspA多肽,再使用熟知的纯化方法,如本文所述的那些。可以直接从已用如本文所述的表达载体转化的细胞中分离天然或突变MspA多肽。还可以将重组表达的天然或突变MspA多肽表达为具有适当亲和标签的融合多肽,所述标签如谷胱甘肽S转移酶(GST)或多组氨酸,并亲和纯化。还可以使用本领域技术人员熟知的多肽合成方法,通过化学合成产生天然或突变MspA多肽。
重构成屏障(例如,膜),如脂质双分子层的多核苷酸孔还可以用于纳米孔测序。多核苷酸孔是延伸穿过屏障(例如,膜)并允许离子和/或水溶性分子从屏障一侧穿过至屏障另一侧的一个或多个多核苷酸。可以根据本发明公开的具体实施方式使用任何多核苷酸孔,只要多核苷酸可以形成允许目标多肽穿过屏障(例如,膜)的缢缩区。示例性多核苷酸孔包括(例如)多核苷酸折纸孔(polynucleotide origami pore)。可以使用“折纸法”制备图形在二维或三维空间中延伸的多核苷酸折纸孔,如Rothemund,Nature,440:297-302(2006)中所述。折纸是利用基因组多核苷酸的长链和多个较短的合成“订书针”多核苷酸链产生延伸结构的通用技术。最初的折纸结构基本上是二维结构。此后,折纸技术扩展至三维结构(Douglas等人,Nature 459:414-418(2009);Ke等人,Nano Letters,6:2445-2447(2009);Andersen等人,Nature 459:73-76(2009))。
在本发明公开的组合物和方法中还可以使用固态孔。固态孔是由非生物来源的材料制成的孔,其延伸穿过屏障(例如,膜),允许离子和/或水溶性分子从屏障一侧穿过至屏障另一侧。
可以通过在固态屏障(例如,膜)中产生孔来形成固态孔。因此,与固态膜类似并且如本文所述,可以通过多种材料形成固态孔,包括无机和有机材料两者。
适合的固态孔包括(例如)氧化铝、氧化钽、氧化钛、二氧化硅、二氧化铪、氧化锆、氮化硼、氮化硅、石墨烯或其纳米层合物(例如,石墨烯-Al2O3)或它们的任意组合(PCT专利公开WO 2013016486A1)。可以通过使用定制的反馈控制离子束雕刻工具或使用来自场发射枪(FEG)TEM的聚焦汇聚电子束在膜中分解溅射(decompositionally sputter)纳米孔或本领域中熟知的任何其它方法来制备固态孔(PCT专利公开WO2013016486A1)。例如,可以通过使用来自FEG TEM的聚焦汇聚电子束钻穿石墨烯-Al2O3膜制备石墨烯纳米层合物孔,如石墨烯-Al2O3孔(Venkatesan等人,ACS Nano.,6:441-450(2012))。
可以在本发明公开的组合物和方法中使用生物和固态杂交孔。生物和固态杂交孔是由生物来源和非生物来源的材料制成的杂交孔,其延伸穿过屏障(例如,膜),允许离子和/或水溶性分子从屏障一侧穿过至屏障另一侧。生物来源的材料如上所定义并且包括(例如)多肽和多核苷酸。非生物来源的材料被称为固态材料,如本文所述。
因此,生物和固态杂交孔包括(例如)多肽-固态杂交孔和多核苷酸-固态杂交孔。多肽-固态杂交孔包括一种或多种多肽和固态材料。多核苷酸-固态杂交孔包括一种或多种多核苷酸和固态材料。通过将多肽或多核苷酸孔与固态孔工程设计来制备生物和固态杂交孔(参见PCT专利公开WO2013/016486)。以上描述了适合的多肽孔、多核苷酸孔和固态孔的实例。
纳米孔测序装置可以具有单孔或多孔。多孔可以作为纳米孔阵列使用以表征不止一个具有相同或不同组成的目标核苷酸。本文所使用的多孔的示例性数目包括(例如)至少1、4、16、64、256、512、1028、4096、16384、32768、100000、100万、1000万个孔或以上。在优选的实施方式中,多孔的数目将大于4096。纳米孔阵列在本领域中是已知的并且在(例如)PCT专利公开WO 2013/016486中公开。例如,可以使用电子束平印术和反应性离子蚀刻步骤在SiN/Al2O3膜中制造~15nm直径的固态孔的高密度阵列,从而有利于多核苷酸分子的高通量分析。
本发明所述的方法可以使用穿过屏障(例如,膜)的势差。所述势差可以是势差、化学势差或电化学势差。可以通过向至少一个液池注入或施用电流的电压电源将势差加载到屏障(例如,膜)两侧。可以通过双池的离子组成的差异,将化学势加载到屏障两侧。可以通过双池离子组成的差异并结合电势来建立电化学势差。不同的离子组成可以是(例如)每个池中不同的离子或者每个池中相同离子的不同浓度。
将电势应用于孔两侧以迫使多核苷酸易位通过孔在本领域中是熟知的并且可以根据本发明公开使用(Deamer等人,Trends Biotechnol.,18:147-151(2000);Deamer等人,Ace Chem Res.,35:817-825(2002);和Li等人,Nat Mater.,2(9):611-615(2003))。可以通过应用于孔两侧的电压实施本发明所述的方法。电压范围可以选自40mV至1V以上。通常,本发明所述的方法将在100至200mV的范围内进行。在具体的实例中,所述方法在140mV或180mV进行。在马达运动期间,电压不需要是稳定的。通常应用电压极性,从而通过电泳驱使带负电荷的多核苷酸进入孔中。在一些情况下,可以降低电压或反转极性以有利于马达的适当功能。
在一些情况下,可以使用压力差的应用以迫使多核苷酸易位通过孔。可以在本文举例说明的方法中电势或其它势差的位置使用压力差。
本发明所述的方法生产了对应于一个或多个核苷酸通过孔的易位的一个或多个信号。因此,例如,当目标多核苷酸通过孔时,由于缢缩的碱基依赖性阻塞,穿过屏障的电流改变。可以使用如本文所述的或本领域中已知的任何多种方法测量来自电流变化的信号。如上所述,每个信号对孔中的核苷酸物质是唯一的,从而所得信号可以用于确定多核苷酸的特征。例如,可以确定产生特征信号的一种或多种核苷酸物质的身份。在本发明所述的方法中有用的信号包括(例如)电信号和光信号,以下将进一步描述。在一些方面,电信号可以是电流、电压、隧穿、电阻、电压、电导率的测量;或者横向电测量(PCT专利公开WO 2013/016486)。在一些方面,电信号是通过孔的电流。
在本文所述方法中检测的电信号可以是流动通过孔的电流,它是电荷流(Deamer等人,Trends Biotechnol.,18:147-151(2000);Deamer等人,Ace Chem Res.,35:817-825(2002);和Li等人,NatMater.,2(9):611-615(2003))。如本文所述,可以使用与孔偶联的检测电路测量电信号,例如,膜片钳电路或隧道电极电路。可以检测的电压、隧穿、电阻和电导率信号的实例和用于它们检测的装置在本领域中是已知的,如(例如)Wanunu,Phys LifeRev.,9(2):125-58(2012);和Venkatesan等人,NatNanotechnol.,6(10):615-24(2011)中所述。
在本发明所述的方法中有用的光信号包括(例如)荧光信号和拉曼信号。可以通过将目标核苷酸与光信号产生标记物(例如,荧光部分或拉曼信号产生部分)偶联来产生光信号。例如,在dela Torre等人,Nanotechnology,23(38):385308(2012)中,使用全内反射荧光(TIRF)显微术的光学方案来照射宽的TiO2-涂覆膜区域。在Soni等人,Rev SciInstrum.,81(1):014301(2010)中,使用方法将两个单分子测量方式结合,即全内反射显微术和使用纳米孔的生物分子电学检测。
如本文所述,可以将孔与传感电路偶联,包括(例如)膜片钳电路、隧道电极电路或横向电导测定电路(如石墨烯纳米带或石墨烯纳米间隙)以记录本发明实施方式中的电信号。另外,还可以将孔与检测多核苷酸上的标记物,例如,荧光部分或拉曼信号产生部分的光学传感器偶联。
纳米孔测序方法可以使用减缓目标多核苷酸通过孔的易位的机制。例如,可以连接或掺入多核苷酸结合蛋白,如解旋酶、易位酶或聚合酶以调节易位速率。所述连接可以是(例如)瞬时的或持久的,并且可以通过目标多核苷酸(由于将其从孔中拉出)或者通过在本领域中熟知的多种多肽、化学接头或捕获部分介导。示例性技术描述于Manrao等人,NatBiotechnol.,30(4):349-353(2012)和Cherf等人,Nat Biotechnol.,30(4):344-348(2102)。在具体的实施方式中,解旋酶或其它分子马达可以用于减缓或停止目标多核苷酸通过孔的易位。例如,当使用水解核苷酸以影响易位的马达时,可以从马达漏去核苷酸和/或可以使马达经受抑制剂(例如,不可水解的核苷酸类似物),从而目标多核苷酸仍与马达结合并且不明显易位通过孔。在一些实施方式中,随后可以通过将核苷酸递送至马达和/或除去抑制剂来允许易位发生。本发明所述的方法可以包括将孔与目标多核苷酸和Hel308解旋酶接触的步骤以控制多核苷酸通过孔的易位速率。如以下进一步描述的,Hel308解旋酶可以鉴别为ATP-依赖性DNA解旋酶和超家族2解旋酶。考虑到本文所提供的教导和指导,对于根据本发明的实施方式的用途,本领域技术人员可以适合地选择或改变任何Hel308解旋酶。以下进一步描述了适合的Hel308解旋酶。
在本文所述方法的一些方面,目标多核苷酸的易位处于通过孔的电流方向的相反方向。在其它方面,目标多核苷酸的易位处于通过孔的电流方向的相同方向。
因此,可以在至少两个模式中进行本发明所述的方法,其中目标多核苷酸的易位与(例如)通过孔的电流或其它势的方向相反或相同。该结果可以通过将本发明公开的Hel308解旋酶与目标多核苷酸的5′或3′端结合来实现。当提及双链多核苷酸时,5′或3′取向是指双链多核苷酸内的单链。因此,Hel308解旋酶可以将多核苷酸从孔中拉出或送入孔中,即与电压梯度对多核苷酸产生的作用力的方向相反(参见图13A-13E和14A-14D),或者当多核苷酸在与电压梯度或其它势所产生的作用力相同的方向移动时,使用解旋酶调节易位速度(参见图15A-15C)。
图13A-13E举例说明了根据一些实施方式,针对势(如电压梯度)所产生的作用力,Hel308解旋酶控制的多核苷酸易位,例如,基于具有Hel308解旋酶3′突出结合位点和胆固醇双分子层锚点的三元多核苷酸复合物。实心圆环(●)表示5′磷酸盐。实心菱形()表示3′胆固醇。具有凹口的实心半透明圆表示Hel308解旋酶。虚线表示任意长度。大的灰色箭头表示多核苷酸进入或移出孔的多核苷酸移动的方向(通过或抵抗外加场)。大的黑色箭头表示解旋酶沿多核苷酸的易位方向,其为3′至5’。孔(漏斗状圆锥形物体)位于膜(双水平线)中。
图13A-13E示出了具有Hel308解旋酶3′突出结合位点和胆固醇双分子层锚点的三元多核苷酸复合物用于多核苷酸序列的使用。胆固醇-标记的多核苷酸“i”是任选的并且用于杂交至目标多核苷酸“ii”,目标多核苷酸“ii”将易位通过孔并且有利于整个复合物向脂质双分子层的招募(图13A)。例如,通过电压梯度将5′磷酸盐拉动穿过孔,从而导致目标多核苷酸“ii”的5′端首先进入孔,并导致胆固醇-标记的多核苷酸剥落(图13B)。随着含有磷酸盐的多核苷酸被拉动通过孔至相反一侧,第二杂交多核苷酸“iii”剥落,这是因为孔过窄从而不允许双链多核苷酸易位(13C)。多核苷酸iii的一个目的是产生Hel308解旋酶可以优先结合的Hel308解旋酶结合位点,它通常是约8个核苷酸的3′单链多核苷酸突出。此外,通过迫使Hel308解旋酶分子在易位多核苷酸的3′端结合,将易位通过孔的多核苷酸的长度最大化。复合物的多核苷酸“iii”可以是任何长度,包括任意长度,并且3′端不需要临近多核苷酸“i”的5′端。一旦到达孔口,Hel308解旋酶通过其3′至5′易位酶活性,将多核苷酸逆电压梯度反向拉入顺室(cis chamber)(图13D和13E)。
图14A-14D还示出了根据一些实施方式,针对势(如电压梯度)所产生的作用力,Hel308解旋酶控制的多核苷酸易位,例如,基于具有Hel308解旋酶3′突出结合位点和胆固醇双分子层锚点的三元多核苷酸复合物。实心圆环(●)表示5′磷酸盐。实心菱形()表示3′胆固醇。具有凹口的实心半透明圆表示Hel308解旋酶。虚线表示任意长度。大的灰色箭头表示多核苷酸进入或移出孔的多核苷酸移动的方向(通过或抵抗外加场)。大的黑色箭头表示解旋酶沿多核苷酸的易位方向,其为3′至5’。孔(漏斗状圆锥形物体)位于膜(双水平线)中。然而,该方案举例说明了使用单个杂交多核苷酸“i”在目标多核苷酸“ii”上产生用于Hel308解旋酶结合的3′突出,并且还含有任选的胆固醇部分。有可能Hel308解旋酶可以结合多核苷酸“ii”的单链区上的任何位置。示出了多个Hel308解旋酶分子,并通过“E1”、“E2”和“E3”表示。到达孔口的解旋酶首先将引起返回同侧(cis side)的控制易位过程。如果它脱离,则将发生不受控制的易位直至下一个结合的Hel308解旋酶分子到达孔口并开始控制的易位。
图15A-15C举例说明了当多核苷酸在与势(如电压梯度)所产生的作用力相同的方向移动时,使用Hel308解旋酶调节多核苷酸易位速度。具有凹口的实心半透明圆表示Hel308解旋酶。虚线表示任意长度。大的灰色箭头表示通过外加场多核苷酸移动进入孔的方向。大的黑色箭头表示解旋酶沿多核苷酸的易位方向,其为3′至5’。孔(漏斗状圆锥形物体)位于膜(双水平线)中。在该示例性方案中,目标多核苷酸的3′端首先进入孔。当沿易位多核苷酸的3′至5′端易位时,Hel308解旋酶控制多核苷酸的易位速度。
如上所述,在Hel308解旋酶的背景中,分步易位步骤可以表示目标多核苷酸的一个或多个核苷酸沿解旋酶和/或孔的部分易位。因此,分步易位步骤是指小于完整易位循环的核苷酸步骤的一部分。当发生构象变化时,在ATP结合和水解之间发生分步易位步骤。完整核苷酸步骤可能需要一个或多个分步易位步骤。构象变化有效地将完整易位循环划分为至少两个部分或分步易位步骤。
可以以相同方式使用部分或分步易位步骤以产生表征通过孔的一个或多个核苷酸的唯一信号。因此,由于对于易位通过孔的一个或多个核苷酸的每一个,电流变化对应于每个分步易位步骤,因此本发明所述的方法可以产生至少两个电信号。因此,在一些方面,分步易位步骤包括Hel308解旋酶完整易位循环的第一分步易位步骤。在其它方面,分步易位步骤包括Hel308解旋酶完整易位循环的第二分步易位步骤。可以单独或与其伴侣(例如,分别为第二或第一分步易位步骤)一起使用第一或第二分步易位步骤中的每一个以表征通过孔缢缩区的一个或多个核苷酸。
例如,如实施例I中进一步描述的,Hel308解旋酶可以结合至ATP并经历构象变化,从而提供第一分步易位步骤,并且Hel308解旋酶可以通过ATP水解将目标多核苷酸的一个或多个核苷酸沿解旋酶和/或孔易位,从而提供第二分步易位步骤。可以使用第一和第二分步易位步骤中的任一个或两者来确定(例如)产生信号的核苷酸或一个或多个核苷酸的核苷酸序列。当通过不止一个核苷酸产生信号时,则产生信号的多核苷酸部分被称为单词。因此,该类核苷酸单词可以是至少4、5、6、7、8、9、10个或更多个核苷酸长度并且对应于孔缢缩区的长度。作为另外一种选择或另外,核苷酸单词可以是至多10、9、8、7、6、5或4个或更少的核苷酸长度。
如上所述和如以下实施例III中进一步举例说明的,可以使用得自完整易位循环的两个部分步骤的电信号鉴别多核苷酸中一个或多个核苷酸残基。使用来自两个分步易位步骤的信号对相同的一个或多个核苷酸提供了重复信号并使得能够在单个确定内获得更大的准确度。因此,使用来自两个部分易位步骤的信号可以导致表征准确度提高,与使用得自完整易位循环的单个电信号或其它信号的一个或多个核苷酸的表征相比,其差错率降低了25至50%。类似地,使用来自部分易位步骤的信号可以导致差错率降低至少5%、10%、20%、30%、40%、45%、55%、60%、65%、70%、75%、80%、85%、90%或更高。考虑到本文所提供的教导和指导,本领域技术人员将知道出于给定目的如何调节准确度,如通过(例如)降低如上所述的缢缩区的尺寸来提高核苷酸易位的分辨率。
在其它实施方式中,得自分步易位步骤的其它信息可以用于在许多方面提高纳米孔测序。例如,对于相同核苷酸单词,可以在算法中使用得自分步易位步骤的测量以改善纳米孔碱基调用准确度。对于相同的核苷酸单词,得自分步易位步骤的测量可以用于降低均聚物读取差错率,这是因为在单次确定中对相同核苷酸单词读取了两次。因此,对于相同核苷酸单词,得自分步易位步骤的测量使天然多核苷酸易位反应的可分辨的分辨率加倍,从而导致序列特异性图形的分辨率提高。后者的一个使用是用于检测序列重复或单核苷酸多态性(SNP)的序列特异性图形识别算法。
如上所述,方法可以包括(a)穿过与Hel308解旋酶和目标多核苷酸接触的孔施加势差;(b)测量目标多核苷酸通过孔的一个或多个分步易位步骤产生的一个或多个信号;和(c)根据分步易位步骤的电信号,表征目标多核苷酸。在一些方面,所述方法还包括重复步骤(a)-(c)一次。通过重复步骤(a)-(c),可以表征相邻核苷酸或相邻核苷酸单词。可以根据需要重复步骤(a)-(c)直至表征了一些或所有目标多核苷酸。例如,可以通过步骤(a)-(c)的任何所需迭代次数确定目标多核苷酸的部分或全部的序列。因此,可以确定目标多核苷酸的整体或部分的一个或多个特征。
如本文所述,根据本发明的实施方式,可以使用任何Hel308解旋酶或其变体。下表1和2中提供了示例性hel308解旋酶。
表1.示例性Hel308解旋酶
下表2示出了Hel308解旋酶以及Hel308基序和扩展的Hel308基序的更多实施方式。
表2.示例性Hel308解旋酶、Hel308基序和扩展的Hel308基序。
**:该Hel308解旋酶序列的相关详细信息,参见国际专利公开No.WO2013/057495。
在本发明实施方式中,还可以使用保留了多核苷酸结合和解旋酶酶活性的Hel308解旋酶的变体或突变体。可以根据本领域熟知的方法获得这种变体或突变体,包括编码天然Hel308解旋酶的核酸的定点突变(Zoller,M.J.,Curr.Opin.Biotechnol.,3:348-354,(1992))。
另外,如上所述并且如本领域中已知的,Hel308解旋酶属于SF2家族并且是3′至5′解旋酶(其还可以称为A型解旋酶)。多个解旋酶的核心域可以包括另一种的常见基序,如RecA结合折叠含有参与核苷酸结合和水解的Walker A基序(其还可以称为基序I)和WalkerB基序(其还可以称为基序II)和基序VI。为了解详情,参见Flechsig等人,“In SilicoInvestigation of Conformational Motions in Super家族2解旋酶蛋白s,”PLoS One:6(7):e21809(2011)。另外,SF2家族的解旋酶可以共有9个保守基序,其可以称为Q、I、Ia、Ib、II、III、IV、V和VI。由于基序II的序列(DEAD(SEQ ID NO:2)或DEAH(SEQ ID NO:3)或DEXH),SF2解旋酶家族还可以称为DEAD-箱(SEQ ID NO:2)蛋白或者DEAH-箱(SEQ ID NO:3)解旋酶。包含在SF2家族中的解旋酶包括RecQ-样家族和Snf2-样家族。多种SF2解旋酶是A型,其中有些例外,如XPD家族。SF2家族的X-射线晶体衍射研究表明解旋酶保守基序在蛋白三级结构中密切相关,并且它们可以形成大的功能域。为了解详情,参见Tuteja等人,“Unraveling DNA解旋酶s:Motif,structure,mechanism and function,”EuropeanJournal of Biochemistry271(10):1849-1863(2004)和Hall等人,“解旋酶motifs:theengine that powers DNA unwinding,”Molecular Microbiology34:867-877(1999)。由Tuteja改编而来的图16示意性示出了已在SF2家族中鉴别的多个基序,例如,DEAD-箱(SEQID NO:2)解旋酶,其中Hel308是成员。如Tuteja所述,开放框(open boxes)代表保守基序。每个解旋酶基序的共有序列由单字母编码表示,例如,图16中“C”可以是D、E、H、K或R;图16中“O”可以是S或T;和图16中“X”可以是任何氨基酸。图16中还示出了分配给基序的名称,例如,Q、I、Ia、Ib、II、III、IV、V和VI。如以上进一步提及的,基序I可以称为Walker A基序,并且在Tuteja中称为ATP酶A WalkerI,基序II可以称为Walker B基序,并且在Tuteja中称为ATP酶B Walker II。箭头所指的基序之间的数字是***基序之间的氨基酸残基的典型范围。
另外,如WO 2013/057495中所述,Hel308解旋酶可以包括氨基酸基序Q-X1-X2-G-R-A-G-R(SEQ ID NO:66),其中X1可以是C、M或L;X1可以是C;X2可以是任何残基,包括疏水性或中性残基,如A、F、M、C、V、L、I、S、T、P或R。任选地,上述基序中末端R可以与P偶联。
考虑到本文所提供的教导和指导,本领域技术人员可以通过确定序列身份或者与以上举例说明的Hel308解旋酶中的一个或多个比对来确定参考解旋酶是否是Hel308解旋酶。
另外,考虑到本文所提供的教导和指导,本领域技术人员可以适合地突变Hel308解旋酶以减缓多核苷酸通过孔的部分易位,例如,通过减缓Hel308解旋酶实施的水解步骤,通过以可以减缓水解的方式与另一种蛋白质的同源基序类似地突变Hel308的基序。作为一个实例,Tanaka等人,“ATPase/解旋酶motif mutants of Escherichia coli PriA蛋白essential for recombination-依赖性DNA replication,”Genes to Cells 8:251-261(2003),描述了在其保守ATP酶/DNA解旋酶基序(即Walker A、B和QXXGRXGR基序)中具有氨基酸取代的Pria蛋白(DEXH-型解旋酶)的突变体。根据Tanaka,某些突变体在某些条件下在水解ATP时高度妥协,并且在某些条件下所有的Walker A和Walker B突变体蛋白显示出高度减弱的DNA解旋酶活力。因此,可以预计与Tanaka中公开的那些类似的Hel308解旋酶的Walker A和Walker B基序的突变可以预期会减弱DNA解旋酶活力或减缓ATP水解,可以预期它减缓多核苷酸通过孔的部分易位并因此提高该多核苷酸的表征。作为另一个实例,Hishida等人,“Role of Walker Motif A of RuvB蛋白in Promoting Branch Migrationof Holliday Junctions:Walker motif A mutations affect ATP binding,ATPhydrolyzing,and DNA binding activities of RuvB,”Journal of BiologicalChemistry 274(36):25335-25342(1999),描述了作为ATP-依赖性六聚DNA解旋酶的大肠杆菌(Escherichia coli)RuvB蛋白的突变体。根据Hishida,Walker基序A的某些点突变影响RuvB的ATP水解和ATP结合活性,以及DNA结合、六聚体形成和促进分枝迁移的活性。因此,可以预计与Hishida中公开的那些类似的Hel308解旋酶的Walker A基序的突变可以预期会影响ATP水解和ATP结合,在某些实施方式中,可以预期它减缓多核苷酸通过孔的部分易位并因此可以提高该多核苷酸的表征。
因此,本发明公开提供了表征目标多核苷酸的方法。所述方法可以包括:(a)穿过与Hel308解旋酶和目标多核苷酸接触的孔施加势差;(b)测量通过所述Hel308解旋酶的所述目标多核苷酸通过所述孔的一个或多个分步易位步骤所产生的一个或多个信号;和(c)根据通过所述分步易位步骤产生的所述一个或多个信号,表征所述目标多核苷酸。
本发明公开还提供了表征目标多核苷酸的方法,其中所述势差包括势差。还提供了表征目标多核苷酸的方法,其中所述信号包括电信号或光信号。所述电信号可以选自电流、电压、隧穿、电阻、电位、电压、电导率和横向电测量的测量。所述电信号包括通过孔的电流。
在其它方面,本发明公开提供了表征目标多核苷酸的方法,其中所述分步易位步骤包括Hel308解旋酶完整易位循环的第一分步易位步骤。所述分步易位步骤还可以包括Hel308解旋酶完整易位循环的第二分步易位步骤。目标多核苷酸的易位可以处于施加于通过孔易位的多核苷酸上的作用力的相反的方向,或者处于施加于通过孔易位的多核苷酸上的作用力的方向。
另外提供了表征目标多核苷酸的方法,其中使用得自完整易位循环的两个分步步骤的电信号表征了目标多核苷酸中的一个或多个核苷酸残基,与使用得自完整易位循环的单一电信号表征一个或多个核苷酸相比,准确度大于50%。
还提供了表征目标多核苷酸的方法,其中所述孔为生物孔。所述生物孔可以是多肽孔或多核苷酸孔。在一些方面,多肽孔具有5个核苷酸或以下的缢缩区。在其它方面,多肽孔包括耻垢分枝杆菌(Mycobacterium smegmatis)孔蛋白A(MspA)。MspA可以具有SEQ IDNO:1所示或与SEQ ID NO:1具有至少15%、至少20%、至少25%、至少30%、至少35%、至少40%、至少45%、至少50%、至少55%、至少60%、至少65%或至少70%同源性的氨基酸序列。
还提供了表征目标多核苷酸的方法,其中所述孔为固态孔或者生物和固态杂交孔。所述生物和固态杂交孔包括多肽-固态杂交孔或者多核苷酸-固态杂交孔。
本发明公开另外提供了表征目标多核苷酸的方法,其中所述Hel308解旋酶是表1和2中所示的解旋酶或其变体。还提供了表征目标多核苷酸的方法,其中所述目标多核苷酸选自单链、双链和部分双链多核苷酸。
在一些实施方式中,根据所述分步易位步骤所产生的一个或多个信号表征多核苷酸包括应用改进型维特比算法。
在一些实施方式中,所述方法还包括步骤(c)之后的(D),改变至少一个参数以改变通过所述Hel308解旋酶的所述目标多核苷酸通过所述孔的一个或多个分步易位步骤的时机;和(e)使用改变的至少一个参数,重复步骤(a)-(c)。所述方法还可以包括合并步骤(c)和(e)期间产生的信号并基于所述合并的信号表征所述目标多核苷酸。在一些实施方式中,所述改变的至少一个参数选自温度、盐浓度、辅因子浓度、ATP产物(如无机焦磷酸)浓度、ADP浓度、pH和所使用的特定Hel308解旋酶。
在一些实施方式中,表征包括检测和鉴别一种或多种信号中的电平以及基于检测和识别的电平确定和输出目标多核苷酸的序列。
说明性地,所述检测和鉴别一种或多种信号中的电平包括输出全电平(fulllevel)、分电平(fractional level)、所有电平和电平标识符(level identifier)中的一个或多个。
基于检测和识别的电平,所述确定和输出目标多核苷酸的序列可以包括将全电平、分电平、所有电平和电平标识符中的一个或多个作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的相关置信信息选择并输出至少一个调用序列。
基于检测和鉴别的电平,所述确定和输出目标多核苷酸的序列可以包括将全电平、分电平、所有电平和电平标识符中的一个或多个作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的部分的相关置信信息选择所调用的多个序列的部分并将其彼此连接。
基于所述检测和识别的电平,所述确定和输出目标多核苷酸的序列可以包括将全电平、分电平、所有电平和电平标识符中的一个或多个作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于所调用的序列与模型序列比较的相关置信信息选择和输出至少一个调用序列。
基于所述检测和识别的电平,所述确定和输出目标多核苷酸的序列可以包括将全电平、分电平、所有电平和电平标识符中的一个或多个作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于所调用的多个序列的部分与模型序列比较的相关置信信息选择所调用的多个序列的部分并将其彼此连接。
本发明公开还提供了调节目标多核苷酸通过孔的分步易位步骤的方法。所述方法可以包括:(a)穿过与Hel308解旋酶和目标多核苷酸接触的孔施加势差;(b)将所述Hel308解旋酶与不同于所述底物参比浓度的Hel308解旋酶底物浓度接触,与所述参比浓度相比,所述底物浓度在分步易位步骤持续时间中所产生的变化与所述底物浓度的差异成正比,和(c)测量目标多核苷酸通过孔的分步易位步骤所产生的信号。步骤(b)可以类似地包括使用底物类似物或抑制剂来实现分步易位步骤持续时间的变化。因此,可以在本发明公开的方法中使用本文所述的或本领域中已知的任何底物类似物或核苷酸抑制剂以用于调节Hel308解旋酶底物、用作参比浓度的Hel308底物中的任一个或Hel308解旋酶底物和用作参比浓度的Hel308底物两者的分步易位步骤。
能够调节多核苷酸的分步易位步骤的Hel308解旋酶底物可以是能够被解旋酶水解的核苷酸或核苷酸类似物。核苷酸底物提供了使双链或部分双链多核苷酸解开或者使单链多核苷酸易位通过孔的能量。Hel308解旋酶的常见底物包括(例如)ATP。Hel308解旋酶底物还包括能够被解旋酶水解的核苷酸和核苷酸类似物。
如本文所述,与核苷酸底物结合有关的一个或多个分步易位步骤的停留时间可以与Hel308解旋酶底物浓度成反比。例如,在一些测试条件下,每个核苷酸易位所观察到的两个分步易位步骤中仅有一个的停留时间与Hel308解旋酶底物浓度成反比。因此,一个易位步骤可以对底物浓度敏感,而另一个易位步骤不敏感。
可以通过改变Hel308解旋酶底物的浓度来实现获得不同分步易位步骤长度的分步易位步骤的调节。可以确定调节程度或幅度,从而本领域技术人员可以选择适合于所需目标多核苷酸表征的特定分步易位步骤长度。可以通过将Hel308解旋酶置于不同于底物参比浓度的Hel308解旋酶底物浓度来确定调节程度。与参比浓度相比底物浓度的变化导致了不同的分步易位步骤停留时间,所述停留时间与相比于参比浓度的底物浓度差异成正比。
因此,可以通过使用不同于底物参比浓度的Hel308解旋酶底物浓度来调节目标多核苷酸通过孔的分步易位步骤。解旋酶溶液内的其它组分或反应条件也可以用于改变分步易位步骤的停留时间,并因此改变单个易位循环的分步易位步骤长度。不同的分步易位步骤可以类似地用于获得额外的信号信息以提高目标多核苷酸表征的准确度。
影响(例如)底物结合至Hel308解旋酶的动力学和解旋酶的底物水解的反应组分和反应条件可以用于改变分步易位步骤的停留时间。这些其它因素包括(例如)反应条件的温度、金属浓度,包括二价金属浓度、离子浓度、溶剂粘度。例如,上述因素和条件以及磷酸盐和/或焦磷酸盐浓度可以影响水解步骤。另外,孔两侧的电压可以(例如)影响底物结合和/或构成Hel308解旋酶停留时间的解旋酶暂停(解旋酶pause)。其它因素包括(例如)pH、阳离子类型或二阶阳离子浓度和类型、解旋酶突变等,所有都可以影响停留时间。在这点上,例如,提高焦磷酸盐浓度可以用于减缓Hel308解旋酶的催化速度并因此提高停留时间。此外,例如,原钒酸钠和腺苷5′-(β,γ-亚氨基)三磷酸锂盐水合物也可以用于减缓解旋酶活性。在以下实施例V中举例说明了调节解旋酶活性的焦磷酸盐和核苷酸类似物的使用。
随着顺序步骤之间电流差的提高,使用分步状态进行数据分析的益处也提高。第一次进行近似法时,分步易位步骤将使用相邻完整易位步骤中间的值。如果分步易位步骤比1/2核苷酸(0.3埃)小得多,则在一些情况下或者甚至在大多数情况下,分步值可能难以或甚至不可能观察到。如果分步易位步骤正好是1/2核苷酸长度,那么平均地,所得电流可以最大程度地不同于对应于完整核苷酸步骤的之前和后续电流值。酶的修饰可以允许聚合物亚基进行纳米分数(fractions of a nanometer)的复位。这可以通过提高或降低酶活性水解位点的相对高度至纳米孔的限制缢缩的酶修饰进行。在一些实施方式中,这可以通过添加或移除解旋酶的氨基酸或用具有较大水力学半径的氨基酸的取代来实现。在其它实施方式中,这可以通过改变氨基酸的电荷来实现,所述氨基酸的电荷可以改变对纳米孔边缘的静电斥力或引力。不希望受任何理论束缚,如果“基于紧握”的假设是正确的(如参考图3更详细地描述的),有可能特定突变将影响解旋酶向上推动解旋酶-多核苷酸复合物的程度,其可以转化为核苷酸Z轴易位百分比的变化。
考虑了精细调节分步易位步骤的持续时间:预期解旋酶ATP酶结构域的某些突变将影响ATP水解速率是合理的。反过来,这将预期会影响分步易位步骤之一的停留时间。例如,如果水解速率缓慢,那么预期会提高分步易位步骤之一的停留时间。其它突变可以影响ATP结合至解旋酶的速率(kon)。在这种情况下,随着ATP结合所花费的时间的提高,相应分步易位步骤的停留时间将提高。
Hel308解旋酶的参比浓度可以是(例如)目标多核苷酸表征中常用的底物的量或者它可以是不同的。例如,如果Hel308解旋酶底物的浓度是常用的1.0mM,那么1mM将对应于参比浓度。参比浓度可以经验性地来源于或得自本领域中熟知的报道。在该具体实例中,除1mM之外的底物浓度将是不同于参比浓度的Hel308解旋酶底物浓度。如以下进一步描述的,可以使用Hel308解旋酶底物的多个浓度和参比浓度以调节或确定分步易位步骤的变化量。
可以改变Hel308解旋酶底物浓度和参比底物浓度,只要两浓度不是饱和浓度。说明性地,Hel308解旋酶底物的饱和浓度为约1mM核苷酸底物。因此,如果参比浓度是1mM,那么要改变的Hel308解旋酶底物浓度可以是小于1mM的任何浓度,其包括(例如)0.1μM、1.0μM、10μM、100μM、200μM、300μM、400μM、500μM、600μM、700μM、800μM、900μM。基于Hel308解旋酶底物浓度和/或参比,其它示例浓度可以是(例如)1.0mM、2.0mM、3.0mM、4.0mM和4.9mM或更小。类似地,只要它们不同,Hel308解旋酶底物和参比底物浓度两者的浓度可以是未饱和的浓度。因此,Hel308解旋酶底物和参比浓度可以是以上所列的任何示例性浓度以及它们之间的任何浓度,例如,0.01μM至5mM的任何浓度范围和该范围之间的所有浓度。
相对于表征目标多核苷酸的方法,可以如上所述实施本发明所述的用于调节分步易位步骤的方法。一旦确定了适合于特定需要的Hel308解旋酶底物浓度,则该底物浓度可以在本文所述的用于表征目标多核苷酸的方法中使用。以类似方式,可以以影响(例如)底物结合动力学和水解的反应组分和条件实施类似确定,以确定适合于特定需要的组分浓度或反应条件。然后,可以在本发明所述的用于表征目标多核苷酸的方法中使用该适合的浓度或条件。新的底物浓度、反应组分浓度和/或反应条件将导致不同的停留时间,其可以提供附件信号信息,从而以参考实施例IX的如下所述的方式用于提高测定的准确度。
因此,本发明公开提供了调节目标多核苷酸分步易位步骤的方法,其进一步包括根据一个或多个分步易位步骤的一个或多个信号,表征目标多核苷酸。表征可以包括鉴别以下中的一个或多个:(1)目标多核苷酸的序列;(2)目标多核苷酸的修饰;(3)目标多核苷酸的长度;(4)目标多核苷酸的身份;(5)目标多核苷酸的来源;或(6)目标多核苷酸的二级结构。
本发明公开还提供了调节目标多核苷酸分步易位步骤的方法,其中所述方法使用了势差,所述势差包括势差。还提供了调节目标多核苷酸分步易位步骤的方法,其中分步易位步骤所产生的信号包括电信号或光信号。另外,提供了调节目标多核苷酸分步易位步骤的方法,所述电信号是选自电流、电压、隧穿、电阻、电位、电压、电导率和横向电测量的测量。所述电信号也可以是通过孔的电流。
还提供了调节目标多核苷酸分步易位步骤的方法,其中底物浓度是Hel308解旋酶底物的半饱和浓度。在一些实施方式中,参比浓度是Hel308解旋酶底物的饱和浓度。在其它方面,底物浓度和参比浓度两者是Hel308解旋酶底物的半饱和浓度。还提供了调节目标多核苷酸分步易位步骤的方法,其中Hel308解旋酶底物是三磷腺苷(ATP)。
还提供了调节目标多核苷酸分步易位步骤的方法,其中所述分步易位步骤包括Hel308解旋酶完整易位循环的第一分步易位步骤或Hel308解旋酶完整易位循环的第二分步易位步骤。目标多核苷酸的易位可以处于施加于通过孔易位的多核苷酸上的作用力的相反的方向,或者处于施加于通过孔易位的多核苷酸上的作用力的方向。
本发明公开还提供了调节目标多核苷酸分步易位步骤的方法,其中使用得自完整易位循环的两个分步步骤的电信号表征目标多核苷酸中的一个或多个核苷酸残基,与使用得自完整易位循环的单个电信号的一个或多个核苷酸的表征相比,其准确度大于50%。在本发明公开所述方法的一些方面,与参比浓度相比,在较低底物浓度时,表征目标多核苷酸中的一个或多个核苷酸残基具有更大的准确度。
还提供了调节目标多核苷酸分步易位步骤的方法,其中所述孔是生物孔。所述生物孔可以是多肽孔或多核苷酸孔。在一些方面,多肽孔具有5个核苷酸或以下的缢缩区。在其它方面,多肽孔包括耻垢分枝杆菌(Mycobacterium smegmatis)孔蛋白A(MspA)。MspA可以具有SEQ ID NO:1所示或与SEQ ID NO:1具有至少15%、至少20%、至少25%、至少30%、至少35%、至少40%、至少45%、至少50%、至少55%、至少60%、至少65%、至少70%、至少75%、至少80%、至少85%、至少90%、至少95%或至少99%同源性的氨基酸序列。
还提供了调节目标多肽分步易位步骤的方法,其中所述孔是固态孔或生物和固态杂交孔。所述生物和固态杂交孔可以是多肽-固态杂交孔或者多核苷酸-固态杂交孔。
还提供了调节目标多核苷酸分步易位步骤的方法,其中所述方法中的Hel308解旋酶包括表1和2中所示的解旋酶或其变体。所述目标多核苷酸选自单链、双链和部分双链多核苷酸。
本发明公开还提供了表征目标多核苷酸的组合物。所述组合物包括包含在小于1mM ATP溶液或核苷酸类似物的溶液中的孔、Hel308解旋酶和目标多核苷酸。在所述组合物的一些方面,小于1mM ATP的溶液是0.1μM、1.0μM、10μM、100μM、0.5mM或0.9mM的ATP。
本发明公开的组合物可以包括以上或以下所述的任何组分,所述组分在本发明所述的方法中用于表征多核苷酸或用于调节目标多核苷酸易位的分步易位步骤。例如,组合物可以包含如上所述的孔。根据本文所提供的教导和指导,所述孔可以是(例如)生物孔,如多肽孔或多核苷酸孔。作为另外一种选择,所述孔可以是如上所述的固态孔或杂交孔。
另外,所述组合物将包含用于表征的目标多核苷酸、Hel308解旋酶和Hel308解旋酶底物。如孔一样,所述目标多核苷酸、Hel308解旋酶和Hel308解旋酶底物可以是本文所述的任何示例性多核苷酸、Hel308解旋酶、底物以及变体和类似物以及本领域中熟知的那些。
因此,本发明公开提供了用于表征目标多核苷酸的组合物,其中所述孔是生物孔。所述生物孔可以是多肽孔或多核苷酸孔。所述多肽孔可以具有5个核苷酸或更小的缢缩区并且可以是耻垢分枝杆菌(Mycobacterium smegmatis)孔蛋白A(MspA)。MspA可以具有SEQID NO:1所示或与SEQ ID NO:1具有至少15%、至少20%、至少25%、至少30%、至少35%、至少40%、至少45%、至少50%、至少55%、至少60%、至少65%、至少70%、至少75%、至少80%、至少85%、至少90%、至少95%或至少99%同源性的氨基酸序列。
还提供了用于表征目标多核苷酸的组合物,其中所述孔是固态孔。还提供了用于表征目标多核苷酸的组合物,其中所述孔是生物和固态杂交孔。所述生物和固态杂交孔可以是多肽-固态杂交孔或者多核苷酸-固态杂交孔。
还提供了用于表征目标多核苷酸的组合物,其中所述Hel308解旋酶是表1和2中所示的解旋酶或其变体。还提供了用于表征目标多核苷酸的组合物,其中所述目标多核苷酸选自单链、双链和部分双链多核苷酸。
应理解在本文所提供的本发明公开的定义内还包括了基本上不影响本发明公开的多个实施方式的活动的修改。因此,以下实施例旨在说明但不限制本发明公开。
实施例I
使用Hel308解旋酶的分步易位步骤
实施例I描述了使用示例性Hel308解旋酶所观察到的分步易位步骤。
脂质双分子层由1,2-二植烷酰基-sn-甘油基-3-胆碱磷酸(Avanti极性脂)组成。双分子层在Teflon中横跨水平直径~20微米的孔。以~2.5ng/ml的浓度,将M2-NNN-MspA加入至双分子层的底侧(grounded side)。一旦***单个孔,则隔室内冲入实验缓冲液以避免进一步***。Axopatch-200B膜片钳放大器(Axon Instruments)将180mV的电压施加到双分子层两侧并测量离子电流。使用4-极贝塞耳滤波器以50kHz低通滤波模拟信号,然后以低通滤波频率的5倍数字化。使用以LabWindoWs/CVI写入的常规软件(National Instruments)控制数据采集。双分子层两侧上~60μl的隔室含有0.3MKCl、1mM EDTA、1mM DTT、10mMMgCl2和10mM HEPES/KOH缓冲液,pH 8.0的实验缓冲液。将野生型Hel 308 Tga或野生型Phi29聚合酶中的任一种用作马达。在存在Hel308 Tga的情况下,缓冲液添加1mM ATP。在存在Phi29的情况下,缓冲液添加dCTP、dATP、dTTP和dGTP各100μM。
图2A-2C示出了根据一些实施方式的Phi29聚合酶和Hel308 Tga解旋酶易位事件的比较。图2A示出了与使用phi29 DNA聚合酶(DNAP)所观察的易位步骤相比,使用Hel308Tga解旋酶所观察的分步易位步骤。易位多核苷酸(SEQ ID NO:67:/5Phos/AAACCTTCCXCCCGTACCGTGCCGTACCGTTCCGTTCCGTACCGTA TTTTTTTTTCTCACTATCGCATTCTCATGCAGGTCGTAGCC,其中X=无碱基)杂交至含有胆固醇的多核苷酸(SEQ ID NO:68:CCTGCATGAGAATGCGATAGTGAGA TTTTTTTTTTTTTTTTTTTT-CholTEG)。使用MspA-M2纳米孔。对Hel308 Tga解旋酶多核苷酸易位所观察到的电平数目几乎是对phi29 DNAP所观察到的电平数目的两倍。在示踪之间划线以表明相应电平。phi29示踪(顶部)是共有事件(consensus),而Hel308解旋酶示踪(底部)是所测量的单个易位事件。共有事件可以表示来自对相同序列的多次读取的可靠检测电平的组合。这种组合潜在可能比单个读取更可靠,这是因为它可能不必需包括错误,如单个分子易位所发生的错误,例如,核苷酸“跳跃”或者核苷酸“切换”,如本领域中已知的。
图2B示出了将使用Phi29聚合酶作为分子马达易位通过MspA-M2纳米孔的单链多核苷酸模板所产生的预期电流电平与使用Hel308 Tga解旋酶作为分子马达所观察到的那些相比,用Hel308 Tga解旋酶观察到的分步易位步骤。易位多核苷酸(SEQ ID NO:69:/5Phos/CTCACCTATCCTTCCACTXXCCCCCTTTGGGTTTAAATTTTTTCAGATCTCACTATC TTTTTAAAGTTTTCTCACTATCGCATTCTCATGCAGGTCGTAGCC,其中X=无碱基)杂交至含有胆固醇的多核苷酸(SEQ ID NO:68:CCTGCATGAGAATGCGATAGTGAGATTTTTTTTTTTTTTTTTTTT-CholTEG)。phi29示踪(顶部)是预测的,而Hel308解旋酶示踪(底部)是所测量的单个易位事件。如共有事件一样,预测可能不必需包括错误,如单个分子易位可能发生的错误,例如,核苷酸“跳跃”或者核苷酸“切换”,如本领域中已知的。预测可以表示基于k聚体表的模拟数据,该表是先前采集的。Phi29预测图形基于每个碱基的完整步骤,并且示出了可以使用完整步骤分子马达预期的一类图形。比较地,明显可以看出Hel308 Tga解旋酶具有分步步骤。
图2C示出了与使用phi29 DNAP所观察的易位步骤相比,使用Hel308Tga解旋酶所观察的分步易位步骤。易位多核苷酸序列是:SEQ ID NO:70:/5Phos/CATCATCATCATCATCATXXCCCCCTAAACAAGAATACCACGA CTAGCATTTTTCAGATCTCACTATCGCATTCTCATGCAGGTCGTAGCC。所述易位多核苷酸杂交至含有胆固醇的多核苷酸(SEQ ID NO:68:CCTGCATGAGAATGCGATAGTGAGATTTTTTTTTTTTTTTTTTTT-CholTEG)并通过MspA-M2纳米孔。使用简单重复序列5′-CAT-3′显示重复图形。对Hel308 Tga解旋酶多核苷酸易位所观察到的电平数目是对phi29 DNAP所观察到的电平数目的两倍。phi29和Hel308解旋酶示踪两者均是共有示踪。共有示踪的使用可以有利于不同分子马达之间易位步骤尺寸的比较,并且可以降低或消除潜在可以使解释复杂化的假象,如跳跃和切换。
不希望受任何理论束缚,在分步易位步骤的进一步解释中,提出了“基于紧握”的机制。图3示出了根据一些实施方式对分步易位步骤所提议的“基于紧握”的机制。多核苷酸(黑色实线)与解旋酶(水平线填充的图形)结合。一旦ATP结合(步骤1),解旋酶经历构象变化(步骤2)。由于解旋酶握紧多核苷酸,因此多核苷酸相对于解旋酶的位置不必需改变。解旋酶上的参考点(灰色三角形)相对于解旋酶握紧的多核苷酸(参见握紧的多核苷酸上的参考点,灰色正方形)不移动。解旋酶构象变化将解旋酶-多核苷酸复合物从纳米孔顶部推出,同时也将多核苷酸拉入孔缢缩中(带有箭头的黑线所指出的黑线)。第二多核苷酸参考点(白色圆圈)显示构象变化期间(步骤2)多核苷酸相对于孔缢缩移动,其导致分步步骤所测量的电流改变。最后,ATP水解并且解旋酶沿多核苷酸易位(步骤3)。这导致多核苷酸相对于解旋酶和孔移动完整核苷酸。总之,在第一分步易位步骤中,Hel308解旋酶结合至ATP并经历构象变化,该构象变化拉动解旋酶握紧的多核苷酸并通过一个核苷酸的一部分使多核苷酸移距,这反过来产生了可测量的电流变化。在第二分步易位步骤中,ATP水解并且Hel308解旋酶完成了一个核苷酸通过纳米孔的易位。可以适合地使用其它机制来解释本发明的分步易位步骤的观察结果。
实施例II
ATP浓度和分步易位步骤之间的关系
实施例II描述了ATP浓度对分步易位步骤停留时间的影响。
为了进一步阐明分步易位步骤的生物化学机制,在不同的ATP浓度下检验分步易位步骤的停留时间。首先用由400mM KCl、10mM HEPES,pH8组成的缓冲溶液填充顺式和反式孔。通过在~25μm直径Teflon孔上涂覆十六烷和脂肪的混合物形成由DPhPC组成的脂质双分子层,并进行电导测定以确保脂质双分子层和Teflon孔之间的千兆欧姆密封(Gigaohmseal)。使用连接到一对Ag/AgCl电极(该电极连接至顺式和反式孔)的Axopatch200B膜片钳放大器进行所有电学测量。膜形成后,将MspA纳米孔注入顺式孔中,其中通过电导测定监控纳米孔向脂质双分子层的掺入。一旦将单个纳米孔掺入到双分子层,灌注顺式室以防止多个孔***。然后,以10nM的最终浓度将单链多核苷酸注入顺式室,对膜施加电压,并通过瞬间电流响应检测多核苷酸通过孔的易位。多核苷酸易位检测时,将电压设置为0V,并且将1mM MgCl2、115nM Hel308解旋酶和各个浓度的ATP(10μM、30μM、100μM和1mM)注入到顺式孔内。然后,将电压设置为保持电位(0.01、0.1和1mM ATP为140mV;0.03mM ATP为180mV),记录电流。如图14A-14D所示,在注入顺式孔内之前,易位多核苷酸(SEQ ID NO:71:/5Phos/CATCATCATCATCATCATXXCCCCCTAAACAAGAATACCACGACTAGCATTTTTCAGATCTCACTATCGCATTCTCATGCAGGTCGTAGCC,其中X=无碱基)杂交至含有胆固醇的多核苷酸(SEQ ID NO:68:CCTGCATGAGAATGCGATAGTGAGATTTTTTTTTTTTTTTTTTTT-CholTEG)(如在本文其它处更详细地描述的)。以这种方式,多核苷酸的5′端首先易位通过纳米孔,然后通过Hel308解旋酶的加工,再拉回通过纳米孔。Axopatch放大器以50kHz的采样率并且以10kHz的低通滤波器记录***的电流响应。在该频率范围内,由于多核苷酸的Hel308解旋酶加工通过纳米孔的步骤转换(包括分步易位步骤)是明显可鉴别的。实验后,使用计算机算法鉴别多核苷酸易位事件。通过使用确定邻近值之间统计学显著性的本领域熟知的学生t检验,鉴别这些易位事件内统计学显著的电流电平(有关进一步详细情况,参见Carter等人,在本文其它处引用,或JohnE.Freund,Mathematical Statistics,第5版,Prentice Hall)。对于从该特定序列所观察的电流,与单一步骤分子马达相比,由于存在加工的核苷酸,因此所鉴别的统计学显著的电流电平几乎是两倍,其中拓扑结构(电流电平的峰值和谷值)在每个峰值之间和每个谷值之间具有几乎两倍的水平,如通过直接观察所测量的。
为了减少实验误差,在核苷酸易位的大分辨率区域中进行多核苷酸通过纳米孔易位的持续时间的数据分析。在链测序中,并且具体地,在多核苷酸的纳米孔测序中,由于当与相邻多核苷酸序列相比时,离子流中的显著性差异阻断,因此无碱基区域可以获得相对高的信噪比。出于该原因,与一些不明确的“噪音”影响相比,无碱基区域的附近区域内的统计学显著性水平潜在可以更可能是由于通过纳米孔的核苷酸加工所造成的。出于该原因,选择了围绕并包含无碱基电流峰值的27个电流电平的持续时间来进行基于电流电平持续时间的数据分析。
图4A和图4B示出了根据一些实施方式ATP浓度对分步易位步骤停留时间的示例性影响。在图4A中,将来自包含分步易位步骤的易位通过纳米孔的多核苷酸序列的无碱基电流峰值的电流电平顺序标记为1-27,并对中值持续时间作图。
通过使用学生t检验的算法检测电流电平以确定相邻电流值之间的统计学显著性(参见Carter等人,在本文其它处引用)。其它技术能够进行该检测,包括速度阈值和卡方(chi-squared)最小化,并且在本领域中这些对于核苷酸加工-相关电流变化以及对于图像处理中的步骤检测均是已知的。存在与这些电平中每一个相关的持续时间,并且在与穿越孔的相同序列的多个多核苷酸的相同电平的比较中,计算每个电平的中值持续时间。因此,这些中值持续时间代表了与每个电平相关的典型持续时间。然而,由于停留时间的指数分布,这些停留时间的时间常数更指示了它们的ATP依赖性。出于该原因,在图4B中,用指数衰减曲线(a*e-t/τ)对偶数和奇数电平(其中“偶数”和“奇数”与图4A中的电平指数有关)的持续时间柱状图拟合,并对每一个的时间常数作图。通过将每个多核苷酸易位事件中的每个电平的持续时间引入到相同大小的容器(equivalently-sized bins)中,构建电平持续时间的柱状图。然后,使用商用曲线拟合算法(the Matlab Curve Fitting Toolbox)对这些柱状图拟合,所述算法使用最小二乘法将数据与指数衰减模型拟合。该方法将残差平方和最小化,其中残差定义为数据点和对该点的拟合反应之间的差异。这是将数据与参数模型拟合的标准技术。图4B中的误差线对应于每次拟合的95%置信界限(confidence bounds)。
如图4B所示,偶数电平的停留时间随ATP浓度降低而升高,而奇数电平的停留时间保持恒定。因此,偶数电平的停留时间对应于第一分步易位步骤,其表面上与ATP结合有关并且与ATP浓度成反比,并呈现指数分布,而奇数电平的停留时间对应于第二分步易位步骤,其表面上与ATP水解有关并且是ATP依赖性的。
实施例III
多核苷酸测序中分步易位步骤的使用
实施例III描述了与使用得自完整易位循环的单个电信号相比,通过使用得自完整易位循环的两个分步易位步骤的电信号测序准确度提高。
由于MspA“读取-头(read-head)”对缢缩区内的4个核苷酸(4-聚体)延伸敏感,因此从测量对应于MspA纳米孔中所见的所有4-聚体组合的电流的四聚体图中产生了电流示踪。有关测量对应于4-聚体组合的电流的进一步详细信息,参见Laszlo等人,“Decodinglong nanopore sequencing reads of natural DNA,”Nature Biotechnology 32:829-833(2014)。然而,应理解不同的孔可以对缢缩区内不同数目的核苷酸敏感。在本实施例中,如下所述通过比较使用原始de Bruijn序列的隐马尔可夫模型(HMM)结果,确定了序列准确度。对于典型的实验噪声水平(如图5中所示,即~0.5-2pA,或约0.5至1.5pA),与分步步骤(正方形)相比,使用完整步骤(菱形)的重建准确度降低。
简要地,使用先前所述的方法建立孔(参见Butler等人,Proc.Natl.Acad.Sci.USA,105:20647-20652(2008);Manrao等人,PLoS ONE,6:e25723(2011))。简要地,从1,2-二植烷酰基-sn-甘油基-3-胆碱磷酸(Avanti极性类脂)在Teflon中形成跨过~20微米直径水平孔的脂质双分子层。双分子层两侧上的隔室含有10mMHepes,pH 8.0、400mM KCl、1mM DTT和10mM MgCl2的实验缓冲液。将Axopatch-200B(AxonInstruments)用于在双分子层两侧施加电压(140mV或180mV)并测量离子电流。以~2.5ng/ml的浓度,将MspA加入至底部顺式隔室(cis compartment)。一旦将单个MspA蛋白***Teflon孔,则顺式隔室内冲入实验缓冲液以抑制或避免进一步***。所有实验在23℃进行。使用4-极贝塞耳滤波器以20kHz对模拟离子电流信号低通滤波并使用NationalInstruments 6363数字化仪以100kHz数字化。使用以LabWindows/CVI写入的常规软件(National Instruments)控制数据采集。使用Matlab写入的常规软件(The Mathworks)分析数据。除了ATP滴定实验外,通常以1mM使用ATP,在这种情况下,ATP浓度范围在10uM至1mM。以10nM使用杂交至含有胆固醇的多核苷酸的易位多核苷酸。以115nM的最终浓度使用Hel308 Tga解旋酶。将多核苷酸和ATP加入至顺式室,最后加入Hel308 Tga解旋酶。作为另外一种选择,可以使用在本领域中熟知的ATP再生***。一个示例性***包括2mM ATP、10mM磷酸肌酸二钠盐、3.5U/mL肌酸激酶和0.6U/mL无机焦磷酸酶。
图5根据一些实施方式,对完全步骤(菱形)和步骤(正方形)绘制了电流迹线电子产生的序列重构准确度(隐马尔可夫模型(HMM))对不同水平的附加噪声的图。图5示出了来源于对于de Bruijn序列(256-聚体)的模型电流阻断示踪的HMM/Viterbi算法分析的序列重建准确度。在一些方面,一般HMM算法类似于Timp等人,Biophys J.2012May16;102(10):L37-9.doi:10.1016/j.bpj.2012.04.009中所述的算法。该算法可以从一系列观察测量恢复潜在的M“状态”组。该算法的基本形式依赖于两个通过实验确定的概率组:状态-状态“转移”概率和状态-观察“发出”概率。对于N个测量,在步骤i=1、2、3……N中给出测量。一个概率组是转换矩阵,其描述了对于给定时间i和状态Si(其中S是M状态组中的状态),后续状态Si+1(其中Si+1不必需是Si)的概率。对于纳米孔***,通过对4nt敏感并且检验4个规范核苷酸(A、C、G、T)的纳米孔,导致44=256个状态,其对应于4nt的每个组合。这些状态中的每一个仅可以转换为4个相邻状态之一。
图6A示出了根据一些实施方式HMM解码纳米孔中序列所需的具有非零概率的状态转移,其中多核苷酸是通过马达酶移动的。所述马达为在1个核苷酸的步骤中移动多核苷酸的phi29 DNAP或类似的酶。图6B示出了根据一些实施方式HMM解码纳米孔中序列所需的具有非零概率的状态转移,其中多核苷酸是通过马达酶移动的。所述马达为使聚合物能够分步移动的Hel308解旋酶或类似的酶。
在图6A中对于在单个核苷酸步骤中移动的酶示出了该***转换矩阵的非零转换概率。使用这类酶,每个多核苷酸状态或n聚体必需达到4个相邻n-聚体状态之一。对于采取一个分步易位步骤的酶,将有更多状态。在这点上,在可以观察到另一种完整步骤状态之前,给定完整步骤状态必须进入半步(或分步步骤)状态。因此,通过更多可辩别的途径,存在更多可用的状态,因此辅助了多核苷酸表征的准确度。
通过q*4n+1提供了状态数,其中n是纳米孔的读取大小,而q是完成完整易位循环所需的步骤数。对于q=2和n=4,如对Hel308解旋酶和M2-NNN MspA所观察的,存在2048个状态。在图6B中图示描述了在分步核苷酸步骤中移动的酶的转移概率矩阵。对应于完整状态的每个状态仅可以转移至1024个“半状态”或“分步状态”之一,而每个半(或分步)状态可以转移至4个不同的状态,其对应于纳米孔读取头内的新的状态。对于HMM解码算法,使用另一个概率组:在时间t的电流测定Ct属于状态Si的概率。通过实验确定或根据先前实验观察估计概率组。可以通过比对算法的迭代应用来完成这种估计,如Laszlo等人2014(在本文其它处引用)中所述或者通过HMM的期望最大化来完成。为了评价分步易位步骤的使用,将使用分步核苷酸步骤的酶的测序准确度与使用单个核苷酸步骤的酶的相比较。通过在MATLAB中进行的常规软件实施HMMViterbi译码算法,并对每个条件产生10个蒙特卡洛模拟电子实验。根据这100个蒙特卡洛模拟的平均和标准偏差获得了测序平均值和标准偏差。基于Manrao 2012(在本文其它处引用)的结果,产生电流电平。以图5中X-轴所示的值给出的高斯宽度,加入高斯噪声以改变在序列重构中使用的电子观察电流值。,并且典型的纳米孔测序实验具有约1pArms的平均水平波动。对于所添加的0.5pA宽度的高斯偏移,分步和完整步骤重构两者均获得了100%的测序准确度测量。在具有高于0.5pA的宽度的添加的高斯偏移之上,分步易位步骤的测序准确度大于非分步易位步骤的测序准确度。因此,当将大于0.5pA的高斯噪声加入至平均电流电平时,分步易位步骤的其它信息提供或给予了提高的序列重构准确度。
除噪声波动之外,由酶的随机运动所造成的跳跃水平将或可以预期会降低测序准确度。相邻四聚体中核苷酸图形的重读会或可以部分补偿这种准确度的降低。通过添加分步易位步骤,存在核苷酸图形的其它重读。例如,在相邻分步步骤中包含了给定k聚体的相关信息,因此在那些相邻分步步骤期间重读了k聚体。例如,假设具有序列ATCGTC的多核苷酸正分布易位通过具有4-核苷酸敏感的读取头的纳米孔。不希望受任何理论束缚,对于完整步进马达,仅当C和G之间的区域位于读取头中心时读取4-聚体TCGT(也就是说,先前“步骤”将使“TC”位于读取头中心,并且仅读取ATCG;下一步将使‘GT’位于读取头中心并仅读取CGTC)。因此,如果马达跳过TCGT读取步骤,则不再测量与特定4-聚体有关的信息。然而,不希望受任何理论束缚,通过分步步进马达,如Hel308解旋酶,在完整步骤期间,两个相邻核苷酸之间的区域可以位于读取头中心,而在分步步骤期间,单个核苷酸可以位于读取头中心。因此当上述多核苷酸的‘CG’位于读取头中心时,读取TCGT,如在完整步进情况中一样。先前分步步骤可以仅具有位于读取头中心的C,并且可以读取ATCGT相关信息;后一分步步骤仅具有位于读取头中心的G,并读取TCGTC相关信息。由于在分步步骤情况中,‘TCGT’相关信息可以读取3次,而在完整步进情况下,仅读取1次,因此4-聚体的这种额外的“重读”可以使得能够获得TCGT相关信息,所述马达跳过了与之相关的步骤-这在完整步进情况中可能是错误的。对于移除水平的所有部分,测序准确度存在2-8%的改善。这通过实施电流电平随机去除的额外的电子蒙特卡罗模拟进行显示。总之,对于纳米孔测序实验中所见的错误形式,测序准确度存在稳健提高。图10、11和12示出了方案,通过该方案,额外的分步易位步骤信息可以用于改善序列准确度。在以下实施例VI中还举例说明了这些方案。所述方案是示例性使用,并且不意欲限制。
实施例IV
分步易位步骤在图形匹配中的使用
实施例IV描述了使用已知算法的分步易位步骤的示例性使用以识别电平。使用动态规划算法,如Needleman-Wunsch比对,额外的水平为在多个水平内准确找到图形提供了帮助。有关Needleman Wunsche比对算法的进一步详细信息,参见Durbin等人,BiologicalSequence Ahalysis,第11版(Cambridge University Press,Cambridge,UK 2006)。除电平电流平均值之外或作为它的替代,使用电平持续时间、电平电流标准偏差或电平分布可以进一步提高图形匹配准确度。在本实施例中,使用Needleman Wunsch比对算法来鉴别对应于包含在对应于1000个碱基的序列的电平内的15个碱基的序列的电平。比较了下列使用:(1)对应于完整核苷酸移动的电平、(2)对应于2个半(或分步)步骤移动的电平、(3)对应于2个半步骤(或分步步骤)移动的电平和持续时间。通过从不同宽度的高斯分布中产生的随机值所提供的值迁移的电平,通过10个蒙特卡罗模拟电子产生观察结果。图7示出了结果,其示出了根据一些实施方式发现作为高斯偏移的函数的电流图形的预期准确度。根据10个蒙特卡罗模拟的平均和标准偏差产生比对准确度的平均和标准偏差。图7中,菱形示出了具有完整核苷酸步骤的马达,而圆圈示出了具有分步易位步骤的马达,正方形示出了具有分步易位步骤结合持续时间值的马达。简要地,将对应于15个核苷酸的电平图形包含在对应于随机1000个核苷酸的序列的电平图形内。电平对应于具有完整核苷酸步骤(菱形)的马达,如phi29 DNAP(仅完整易位步骤)或对应于具有分步易位步骤的马达(圆形),如Hel308解旋酶。除了电流值之外,使用持续时间以进一步改善匹配(正方形)。根据图7中的结果,可以理解对于提高的噪音,对于使用分步易位步骤移动的算法,匹配质量显著更好。如果还使用持续时间值,则将进一步改善匹配质量。为了将电平与Needleman Wunsch算法相匹配,使用输入水平相似性测量或得分比较水平。在测试中,使用学生t检验比较电流电平。为了比较两个持续时间的相似性(或对其打分),确定了持续时间的自然对数差异,并将其加入到学生t检验给出的得分中。可以在Needleman Wunsch算法的术语表中定义术语“得分”。这些打分功能代表了可以用于比较信号电平(例如,电流值)和持续时间的方法的非限制性实例。
实施例V
调节解旋酶分步步骤
实施例V举例说明了改变反应组分对改变Hel308解旋酶停留时间的使用。
图8示出了根据一些实施方式随焦磷酸盐浓度的改变,Hel308解旋酶活性的示例性调节。图9示出了根据一些实施方式用核苷酸抑制剂原钒酸钠和用核苷酸类似物腺苷5′-(β,γ-亚氨基)三磷酸锂盐水合物对Hel308解旋酶活性的示例性调节。
通过提高焦磷酸盐浓度调节Hel308解旋酶活性。简要地,反应条件为实施例III中所述的那些。通过包含0至50mM范围内的多个浓度的焦磷酸盐,例如,0mM(对照)、5mM、10mM、20mM、30mM、40mM和50mM。图8示出了结果,并且示出了与在不存在焦磷酸盐(对照)的情况下的解旋酶活性相比的解旋酶活性的百分比。5和10mM焦磷酸盐浓度导致解旋酶活性降低至对照的大于75%。大于10mM的焦磷酸盐浓度导致解旋酶活性进一步降低,并因此降低了解旋酶的停留时间。使用荧光测定监控解旋酶打开双螺旋DNA的能力。49-nt FRET多核苷酸(50nM最终浓度)包含5′荧光素基团(/FAM/)。含有40-nt猝灭剂的多核苷酸(50nM最终浓度)包含荧光猝灭剂Black Hole Quencher(/BHQ1/)。使用本领域熟知的方法,通过将它们加热超过它们的解链温度至75℃并缓慢冷却至室温将这两条多核苷酸杂交在一起。双螺旋包括3′至5′解旋酶可以结合的9-碱基3′突出。与含有40-nt猝灭剂的多核苷酸100%互补的互补FRET 40-nt多核苷酸以至少10倍摩尔过量存在。由于起初猝灭剂和荧光团非常接近,因此荧光淬灭。基于打开双螺旋DNA的解旋酶,含有40-nt猝灭剂的多核苷酸与再结合至49-ntFRET多核苷酸相比,更可能结合至互补FRET 40-nt多核苷酸。因此,在存在适合的激发光源的情况下,新的单链49-nt FRET多核苷酸发荧光。测定缓冲液包含10mM HEPES,pH 8.0、400mM KCl、1mM MgCl2、1mM DTT、1mM ATP。读取荧光前,使反应在室温下进行20分钟。
在分别存在核苷酸抑制剂或类似物原钒酸钠和腺苷5′-(β,γ-亚氨基)三磷酸锂盐水合物的情况下,Hel308解旋酶活性并且因此停留时间也显示降低。简要地,反应条件是实施例III中所述的那些。在5mM抑制剂或类似物浓度下,包含原钒酸钠(图9中的“原钒酸钠”)或腺苷5′-(β,γ-亚氨基)三磷酸锂盐水合物(图9中的“AMP-PNP”)。图9示出了结果,并且示出了与在不存在核苷酸抑制剂或类似物(对照)的情况下的解旋酶活性相比的解旋酶活性的百分比。5mM抑制剂或类似物浓度导致解旋酶活性降低至大于对照的85%,并因此可以预期会提高解旋酶停留时间或解旋酶沿DNA移动所花的时间。例如,停留时间增加可以延长分步步骤的时间,因此允许更长的信号采集时间。
实施例VI
处理分步步骤信息以改善测序准确度的方法
实施例VI举例说明了处理得自分步易位步骤的额外信息以改善测序准确度的三种方法。
图10示出了使用电流电平和持续时间信息处理得自分步易位步骤的额外信息的示例性方法。所述方法可以应用于两个独立序列读取。使用该方案,对电流示踪进行步骤检测算法,其中发现了电流电平和那些电平的持续时间。至少部分基于所述电平的持续时间,两状态HMM将电平鉴别为完整步骤(长)或者半(或分步)步骤(短)或者观察结果内的电位跳跃。然后,通过HMM、Viterbi或图形匹配算法或其适当组合使用这些鉴别的长和短步骤以及跳跃信息来分别对两个状态类型重建多核苷酸序列(长和短分别对应于完整步骤和半(或分步)步骤电平)。然后,比较调用序列并用于改善多核苷酸测序准确度,例如,通过调节HMM、Viterbi或图形匹配算法。比对可以用于较差地鉴别两个独立序列读取的匹配位置。
图11示出了使用电流电平和持续时间信息处理得自分步易位步骤的额外信息的示例性方法。所述方法可以应用于两个并行的序列读取。在该方法中,首先对电流示踪进行步骤检测算法以发现电平。然后,将平均(或中值)电平电流值以及每个电平的持续时间成对输入到检验持续时间和电流值并对半(或分步)状态以及完整状态估计或调用最优序列的二维HMM、Viterbi或图形匹配算法或其适合的组合中。在该技术中,HMM发射概率不是二维的:Pi(发射t)=Pi(curt,durt)=Pi(curt)*Pi(durt),其中i是对应于多核苷酸的长步骤或短(分步)易位步骤的“状态”,curt和durt分别是电平数t的电平电流和持续时间。二维HMM可以作为输入采取长电平(完整状态)的共有图谱和概率分布以及短电平(半(或分步)状态)的共有图谱和概率分布。二维HMM可以作为输出提供核苷酸序列的调用。
图12示出了直接使用电流示踪处理得自分步易位步骤的额外信息的示例性方法。可以通过或不通过使用持续时间信息应用所述方法。通过参考持续时间信息的使用,在该方法中,使用持续时间-依赖性HMM直接分析电流示踪。在这种HMM形式中,同时将电平持续时间确定为最可能的序列和完整或者半(或分步)步骤状态。如果在两个时间迭代之间状态保持不变,则将提高给定状态的持续时间。然后,该持续时间用于改善状态是处于完整状态还是分步状态的评价。
实施例VII
处理分步步骤信息以改善测序准确度的其它方法
实施例VII描述了处理分步步骤信息以改善测序准确度的其它示例性方法
基于来自使用单步分子马达易位通过纳米孔的多核苷酸的信号,先前已将隐马尔可夫模型(HMM)和Viterbi算法用于碱基调用(base-calling)。为了解详情,参见Timp等人,“DNA Base-Calling from a Nanopore Using a Viterbi Algorithm,”BiophysicalJournal 102:L37-L39(May 2012)。图19A示意性示出了来自用于表征多核苷酸通过孔的单步易位的信号的示例性隐马尔可夫模型(HMM)的方面,例如,其中给定信号电平对应于一个核苷酸通过孔的易位,例如,通过聚合酶或解旋酶通过孔的易位。如本文在其它处所提及的,信号电平可以不必需对应于孔缢缩内单个核苷酸的存在,但是可以对应于包含多个核苷酸,例如,2、3、4、5、6、7、8、9、10或10个以上的核苷酸的“单词”的存在。这个“单词”也可以称为“k-聚体”。在图19A所示的实施方式中,所述“单词”或“k-聚体”是4个核苷酸长或者是“四聚体”或“4-聚体”,其对应于基于孔缢缩中4个核苷酸存在的信号电平。
在图19A中,可以看出对于易位通过孔的多核苷酸的给定位置i,孔缢缩中给定的四聚体可以包括四种核苷酸的任何可能组合,例如,AAAA、AAAC、AAAG、AAAT……TTTT。基于对应于该四聚体的信号电平来唯一地鉴别该四聚体可以不必是可能的。例如,预期序列中两个不同的四聚体,例如,彼此临近的两个不同的四聚体潜在地可以具有与另一个相同的信号电平。Timp公开了DNA三联体(3-聚体)的示例性电流值,基于此可以看出某些三联体可以具有与另一个相同的信号电平,因此抑制了碱基的核苷酸碱基调用,其中所述三联体仅基于对应于该三联体的电流电平。应理解某些四聚体(并且,更一般地,某些k-聚体)可以具有彼此难以区分的信号电平,因此抑制了碱基的核苷酸碱基调用,其中四聚体或k-聚体仅基于对应于该4-聚体或k-聚体的电流电平。因此,使用HMM的术语,基于信号电平的观察,彼此难以区分的这类四聚体或k-聚体中的碱基可以作为“隐状态”建模。
可以使用基于孔缢缩中多核苷酸的其它单步位置的观察结果的其它信息,从而提高准确鉴别该四聚体或k-聚体中的碱基,并因此准确鉴别“隐状态”的可能性。例如,在图19A中,还可以看出对于易位通过孔的多核苷酸的下一位置i+1,孔缢缩中给定四聚体可以仅具有4个核苷酸的某些可能组合,这是因为位置i的最后三个核苷酸对应于位置i+1的前3个核苷酸。照此,对于i和i+1状态的信号测量可以用于提高正确鉴别在i和i+1位置(或者,同样地,i-1和i位置)中的一个或两个处存在的四聚体的可能性。例如,基于对应于多核苷酸的位置i的序列AAAA,对于位置i+1,仅四个序列AAAA、AAAC、AAAG和AAAT是可用的。对于位置i的每个可能序列,可以容易地鉴别位置i+1的可用的4个序列。类似地,基于多核苷酸位置i+1处的序列,可以容易地鉴别多核苷酸位置i+2处可用的4个序列。用于单步马达的Viterbi算法——其中信号电平和位置i、i+1、i+2、……、i+n之间存在一对一的对应,其中n是多核苷酸中核苷酸的数目,可以表示来自有序电平组L={l1,l2,...ln}的信号。对应于多核苷酸位置i的每个电平li可以表示为该信号电平的平均值(平均值i)、该信号电平的标准偏差(stdi)或该信号电平的持续时间(duri)中的一个或多个。当对应于多核苷酸现有位置(i位置)的四聚体是q时,可能的四聚体组可以表示为定义对应于多核苷酸前一位置(i-1位置)的四聚体的可能值的prev(q)={q1,q2,...q4}。例如,prev(AACC)={AAAC,CAAC,GAAC,TAAC}。基于对应于jth观察电平的观察信号电平Oji,位置i处存在的给定四聚体q的可能性得分可以表示为:
其中s(li|q)对应于表示观察水平li,给定的四聚体q的可能性的奖励得分,InsPen是***罚分(罚分对应于所观察到的但不对应于多核苷酸中四聚体的信号电平),DelPen是缺失罚分(罚分对应于多核苷酸中不具有对应信号电平的四聚体)。
图19B示意性示出了根据一些实施方式用于表征来自使用Hel308解旋酶的多核苷酸通过孔的分步易位的信号的示例性HMM的方面。在图19B中,还可以看出对于易位通过孔的多核苷酸的给定位置i,孔缢缩中给定的四聚体可以包括四种核苷酸的任何可能组合,例如,AAAA、AAAC、AAAG、AAAT……TTTT。可以使用基于孔缢缩中多核苷酸的分步步骤位置以及其他单步位置的观察结果的其它信息,从而以相对于仅使用单独的单步位置的改善的准确度,提高准确鉴别该四聚体或k-聚体中的碱基,并因此准确鉴别“隐状态”的可能性。
例如,在图19B中,还可以看出对于分步步骤马达,易位通过孔的多核苷酸的下一位置是“i分步”并且其中孔缢缩中给定四聚体可以仅具有4个核苷酸的某些可能组合,这是因为位置i的最后三个核苷酸对应于位置“i分步”的前3个核苷酸。照此,i和i分步状态的信号测量可以用于提高正确鉴别所提供的四聚体的可能性。例如,基于对应于多核苷酸位置i的序列AAAA,对于位置i分步,仅4个序列AAAA、AAAC、AAAG和AAAT是可用的。对于位置i处每个可能的序列,可以容易地鉴别位置i分步处可用的4个序列。
另外,在图19B中,还可以看出对于易位通过孔的多核苷酸的下一位置i+1完整(该位置紧挨着i分步),孔缢缩中给定四聚体可以仅具有1个可能序列,这是因为位置i+1完整的4个核苷酸对应于位置i分步的相同核苷酸。照此,对应于i、i分步和i+1完整位置的信号测量可以用于提高正确鉴别i、i分步和i+1完整位置(或者,同样地,i-1和i位置)中一些或全部处存在的四聚体的可能性。例如,基于对应于多核苷酸位置i的序列AAAA,对于位置i分步和对于i+1完整,仅4个序列AAAA、AAAC、AAAG和AAAT是可用的。用于分步步骤马达的改进型Viterbi算法-其中信号电平和分步步骤以及完整步骤位置i、i分步、i+1完整、i+1分步、i+2完整、i+2分步、……、i+n分步、i+n完整之间存在对应关系,其中n是多核苷酸中核苷酸的数目,可以作为电平组表示信号电平l。类似地,如以上参考图19A所讨论的,对应于i完整或i分步位置的每个信号电平li可以表示为该信号电平的平均值(平均值i)、该信号电平的标准偏差(stdi)或该信号电平的持续时间(duri)中的一个或多个。考虑到四聚体q,在电流分步易位步骤中,对应于先前完整位移步骤的可能四聚体组可以定义为prev(q)={q1,q2,...q4}。例如,prev(AACC)={AAAC,CAAC,GAAC,TAAC}。
基于对应于位置i的观察信号电平li,对于对应于完整易位状态的位置,位置i处存在的给定四聚体q的可能性得分scoref以及对于对应于半(或分步)易位状态的位置,位置i处存在的给定四聚体q的可能性得分scoreh可以如下表示为:
其中sf(li|q)对应于表示观察电平li,完整易位状态的给定四聚体q的可能性的奖励得分,sh(li|p,q)对应于表示观察电平li,分步易位状态的给定四聚体q和前一四聚体p的可能性的奖励得分,InsPen是***罚分(罚分对应于所观察到的但不对应于多核苷酸中四聚体的信号电平),DelPen是缺失罚分(罚分对应于多核苷酸中不具有对应信号电平的四聚体)。
另外,动态规划可以用于分步步骤分子马达(如Hel308解旋酶)的图形匹配。在Laszlo等人,“Decoding long nanopore sequencing reads of natural DNA,”NatureBiotechnology 32:829-833(2014)中对于单步分子马达描述了动态图形匹配。例如,对于单步分子马达,信号电平l可以表示为电平组L={l1,l2,…ln},其中对应于多核苷酸的完整易位步骤位置的每个信号电平li可以表示为该信号电平的平均值(平均值i)、该信号电平的标准偏差(stdi)或该信号电平的持续时间(duri)中的一个或多个。基于观察的信号电平li,测量的给定四聚体qj的可能性得分可以表示为:
其中i代表电平序列中的位置;j代表DNA序列中的位置,四聚体qj的最后一个碱基将是位置j处的碱基;得分(i,j)代表电平l1…li和四聚体q1…qj之间匹配的良好程度;s(li|qj)对应于表示观察电平li,给定四聚体qj的可能性的奖励得分;InsPen是***罚分(罚分对应于所观察到的但不对应于多核苷酸中四聚体的信号电平);DelPen是缺失罚分(罚分对应于多核苷酸中不具有对应信号电平的四聚体)。
对于分步步骤分子马达(如Hel308),信号电平l可以表示为电平组L={l1l2,…ln},其中对应于多核苷酸的i完整或i分步位置的每个信号电平li可以表示为该信号电平的平均值(平均值i)、该信号电平的标准偏差(stdi)或该信号电平的持续时间(duri)中的一个或多个。基于观察的信号电平li,对应于完整易位状态的测量的给定四聚体qj的可能性得分scoref以及对应于半(或分步)易位状态的测量的给定四聚体qj的可能性得分scoreh可以如下表示为:
其中i代表电平序列中的位置;j代表DNA序列中的位置,四聚体qj的最后一个碱基将是位置j处的碱基;scoref(i,j)和scoreh(i,j)代表分别假定为完整或分步状态的电平l1…li和四聚体q1…qj之间的匹配良好程度;sf(li|qj)和sh(li|qj)分别对应于表示观察电平li,处于完整或分步状态的给定四聚体的可能性的奖励得分;InsPen是***罚分(罚分对应于所观察到的但不对应于多核苷酸中四聚体的信号电平);DelPen是缺失罚分(罚分对应于多核苷酸中不具有对应信号电平的四聚体)。
现将参考图20A描述使用分步步骤的一些示例性从头测序结果。基于人DNA产生了75个500-聚体多核苷酸文库,并且如在本文其它地方所述的参考实施例II和III类似地采集纳米孔数据。使用以上方程(2)和(3),使用所述的改进型Viterbi算法,分析基于数据的核苷酸碱基调用。然后,将碱基调用序列与150个500-聚体组比对,其中75个是正确的500-聚体,并且其中75个是“诱饵”或“虚假”500-聚体序列。在图20A中,该图示出了作为比对准确度函数的读取长度(使用LASTAL aligner,如Kielbasa等人,“Adaptive seeds tamegenomic sequence comparison,”Genome Research21:487-493(2011)中所述),空心菱形对应于其中碱基调用序列与正确(“目标”)序列比对的结果,而实心菱形对应于其中碱基调用序列与“诱饵”或“虚假”序列比对的结果。根据图20A可以理解对于大于约200个基底对的读取长度,可以获得大于约60%的准确度。还可以使用已知的技术,如读取DNA的两条链来潜在地提高准确度。
现将参考图20B-20C描述使用分步步骤的一些示例性图形匹配结果。如上所述,参考图20A使用75个500-聚体多核苷酸的相同文库和相同实验规程。使用以上方程(5)和(6),使用所述的用于图形匹配的动态规化,分析基于数据的核苷酸碱基调用。然后,将碱基调用序列与150个500-聚体组比对,其中75个是正确的500-聚体,并且其中75个是“诱饵”或“虚假”500-聚体序列。在图20B中,该图示出了作为比对得分的函数的比对尺寸,空心菱形对应于其中碱基调用序列与正确(“目标”)序列比对的结果,而实心菱形对应于其中碱基调用序列与“诱饵”或“虚假”序列比对的结果。根据图20B可以理解对于大于约200个碱基对的比对尺寸可以获得大于约40的比对得分。在图20C中,该图也示出了作为比对得分的函数的比对尺寸,空心菱形对应于其中碱基调用序列与正确(“目标”)序列比对的结果,而实心菱形对应于其中碱基调用序列与“诱饵”或“虚假”序列比对的结果。根据图20C可以理解对于大于约50个碱基对的比对尺寸可以获得大于约20的比对得分。可以观察到分步步骤模型可以比单步模型准确鉴别到更多事件。
另外,可以看出对于使用1332个电平的易位事件,针对80kb数据组的图像匹配(方程5和6)对单链花费了145秒,而使用1332个电平的相同事件,针对该数据组的从头测序(方程2和3)对单链花费了69秒。可以看出图形匹配复杂性随核苷酸数据组线性增加,而从头序列的复杂性独立于数据组。观察图形匹配以准确鉴别较短的事件,这是从头测序不能鉴别的。另外,可以看出图形匹配的分步步骤模型比单步模型产生了更多真阳性结果,这表明分步步骤模型可以是更好的模型以解释解旋酶数据。
实施例VIII
使用其它Hel308解旋酶的分步易位步骤
实施例VIII描述了使用示例性Hel308解旋酶作为分子马达所观察到的分步易位步骤。
如上所述,使用DphPC脂质双分子层中的单个2NNN MspA纳米孔并且使用下表3中所列的参数,参考实施例I类似地实施实施例VIII的实验,其中“Hel308 Mbu(A)”是指用于使用Hel308 Mbu的第一实验的一组参数,并且“Hel308 Mbu(B)”是指用于使用Hel308 Mbu的第二实验的一组参数。脂质双分子层由1,2-二植烷酰基-sn-甘油基-3-胆碱磷酸(Avanti极性脂)形成。双分子层在Teflon中横跨水平直径~20微米的孔。以~2.5ng/ml的浓度,将M2-NNN-MspA加入至双分子层的底侧(grounded side)。一旦***单个孔,则隔室内冲入实验缓冲液以避免进一步***。Axopatch-200B膜片钳放大器(Axon Instruments)将180mV的电压施加到双分子层两侧并测量离子电流。使用4-极贝塞耳滤波器以50kHz低通滤波模拟信号,然后以低通滤波频率的5倍数字化。使用以LabWindows/CVI写入的常规软件(National Instruments)控制数据采集。双分子层两侧~60μl的隔室含有适当浓度的KCl、1mM EDTA、1mM DTT、1mM ATP、5mM MgCl2和10mM HEPES/KOH缓冲液,pH 8.0的实验缓冲液。以对分子马达指明的浓度使用野生型Mbu Hel 308解旋酶。
在Hel308 Mbu实验和Hel308 Tga实验中,以3′至5′的方向读取DNA,而在Phi29聚合酶实验中,以5′至3′的方向读取DNA。
表3
图17A-17D示出了根据一些实施方式使用某些参数的Hel308 Mbu解旋酶、Hel308Tga解旋酶和Phi29聚合酶易位事件的比较。图17A示出了使用表3中所示的“Hel308 Mbu(A)”参数,通过Hel308 Mbu解旋酶观察到的易位步骤。易位多核苷酸(SEQ ID NO:72:/5Phos/AAACCTTCCXCCCGTACCGTGCCGTACCGTTCCGTTCCGTACCGTATTTTTTTTTCTCACTATCGCATTCTCATGCAGGTCGTAGCC,其中X=无碱基)杂交至含有胆固醇的多核苷酸(SEQ ID NO:73:AAAAAAAATACGGTACGGAACGGAACGGTACGGCACGGTACGGGTTTTTTTTTTTTTTTT/3CholTEG)。图17B示出了使用表3中所示的“Hel308 Mbu(B)”参数和使用图17A中相同的多核苷酸序列通过Hel308Mbu解旋酶所观察到的易位步骤。图17C示出了使用表3中所示的“Hel308Tga”参数和使用图17A中相同的多核苷酸序列通过Hel308 Tga解旋酶所观察到的易位步骤。图17D示出了使用表3中所示的“phi29”参数和使用图17A中相同的多核苷酸序列通过Phi29聚合酶所观察到的易位步骤;对纵轴反映图17D中的phi29图,以有利于图17A、17B、17C和17D之间的比较。
从图17A-17D可以看出,对于使用每种解旋酶的测序,纳米孔检测了通常表示为“a”(对应于信号谷值)和“b”(对应于信号峰值)的特征。还可以看出对于使用Hel308 Tga解旋酶的测序(图17C),相比于phi29解旋酶(图17D),观察到了约两倍的电平数。还可以看出对于在“Hel308Mbu(B)”条件下使用Hel308 Mbu解旋酶的测序(图17B),相比于在“Hel308Mbu(A)”条件下使用Hel308 Mbu解旋酶的测序(图17A),观察到了更多的电平数。还可以看出对于在“Hel308 Mbu(B)”条件下使用Hel308 Mbu解旋酶的测序(图17B),相比于Hel308 Tga解旋酶(图17C),观察到了较少的电平,而相比于phi29解旋酶(图17D),观察到了较多的电平。图17A-17D可以理解为表示(1)Hel308解旋酶(例如,Tga和Mbu两者)的多个变体示出了分步步骤,而对于聚合酶Phi29未观察到分步步骤;和(2)通过改变环境变量或参数,例如,KCl浓度,可以阐明分步步骤。另外,其它数据表明当使用Mbu时,电平持续时间随ATP浓度的降低而增加,例如,Mbu中的持续时间以及因此所造成的分布步骤的物理机制也可以是ATP-依赖性的。
实施例IX
应激物的使用,任选地与多种方式结合
基于本文所提供的发明公开应清楚多种环境变量或参数可以影响纳米孔***如何读取或如何基于特定多核苷酸序列产生信号。可以提供这种影响的示例性变量或参数可以包括温度、盐浓度(例如,Mg、Cl)、辅因子(例如,ATP)浓度、ATP产物如焦磷酸盐的浓度、pH、所使用的具体的分子马达(例如,所使用的具体的Hel308解旋酶)、压力等。
例如,如上所述参考实施例II和图4A和4B,ATP的浓度可以影响对应于某些易位步骤的电平的停留时间。例如,观察到第一分步易位步骤的停留时间随ATP浓度的降低的增加,并且表面上与ATP结合有关并与ATP浓度成反比。作为另一实例,如上所述参考实施例V和图8,焦磷酸盐的浓度可以影响Hel308解旋酶的活性。例如,观察到Hel308解旋酶的活性随焦磷酸盐浓度的升高而降低,因此提高了解旋酶的停留时间。作为另一实例,如上所述参考实施例V和图9,核苷酸抑制剂或类似物的浓度可以影响Hel309解旋酶的活性。例如,观察到基于原钒酸钠或腺苷5′-(β,γ-亚氨基)-三磷酸锂盐水合物(AMP-PNP)的存在,Hel308解旋酶活性降低,因此提高了解旋酶的停留时间。作为另一实例,如上所述参考实施例VIII和图17A和17B,盐浓度可以影响所观察到的电平的数目。例如,在使用Hel308 Mbu解旋酶测序期间,观察到盐(例如,KCl)浓度升高提高了所观察到的电平的数目。本领域技术人员将能够容易地设想调整任何适合的参数以调整基于多核苷酸序列产生信号的方式。
另外,应理解这些参数的不同组合可以影响测序准确度以及测序通量。例如,解旋酶停留时间的提高可以提高准确度,例如,可以提高所观察到的电平的数目,但是潜在地可以降低测序通量。对于基于分步步骤观察的测序,与另一组步骤相比,一些步骤潜在地可以受特定变量影响更大。变量-独立步骤可以用于设置准确度基线,同时可以调整其它步骤以满足特定测序需要(例如,提高准确度而降低通量,或者提高通量而降低准确度)。在一些实施方式中,多方式装置可以通过基于测序仪的需要来调整准确度和通量而进行利用,例如,通过在测序期间调整一个或多个参数。作为一个非限制性、说明性实例并且如上所述,已观察到使用Hel308 Tga,ATP浓度的降低可以提高分步状态的持续时间。分步状态持续时间的增加可以提高测序准确度,例如,通过改善分步状态读取的信噪比(SNR)或允许应用低频滤波器,但是它可以降低通量。多方式装置可以通过以高浓度ATP开始测序,从而相对快地确定序列的大致“结构(scaffold)”,然后可以降低ATP浓度以通过更高质量,尽管更慢的读取来“填充结构中的空隙”。
另外,注意可以顺序或彼此平行地使用任何适合个数的不同参数以提高通过Hel308解旋酶的目标多核苷酸通过孔的易位所产生的一个或多个信号的分辨率。图21A-21C示意性示出了根据一些实施方式对于多核苷酸通过孔的不同易位,可以作为时间函数产生的信号。图21A-21C分别示出了虚线曲线,其对应于其中在仅有势差所施加的外加力的情况下,而不是通过分子马达,多核苷酸易位通过孔的情况下所产生的理想化信号,并且其具有无穷大的信号分辨率。在这些条件下,当它们穿过孔时,信号是核苷酸位置和序列的连续改变函数。
图21A也示出了仅使用在垂直虚线所指明的时间发生的完整易位步骤所产生的示例性信号(粗实线)。所述信号可以是电信号或光信号,如本文其它处所述。另外,所述信号可以包括这种电信号或光信号的任何适合的特征,如平均信号电平、信号持续时间或标准偏差(例如,宽带噪声或限带噪声)。在图21A中可以看出通过单个步骤信号从相对低电平改变为相对高电平,然后再通过单个步骤改变为相对低电平,其对应于多核苷酸易位通过孔的完整易位步骤间发生的转换。在图21A中还可以看出信号在不同时间在点(a)、(b)和(c)处与理想化信号相交,并因此可以认为在这些点对理想化信号“采样”。然而,由于有效采样率相对较低,因此所述信号对理想化信号的采样相对较差。例如,点(a)和(b)处的值彼此是相同的,其对应于不同易位步骤的简并信号电平。由于信号没有对点(a)和(b)之间存在的理想化曲线部分充分采样,因此对应于点(a)和(b)的物理易位步骤可以是彼此难以区分的,从而导致多核苷酸序列相关信息的损失。另外,由于信号没有对点(a)和(b)之间存在的理想化曲线部分充分采样,因此对应于点(b)和(c)之间理想化曲线的下降斜率的物理易位步骤潜在地仅可以部分表征这些步骤期间易位通过孔的多核苷酸部分。
除了如上所述由虚线曲线表示的理想化信号之外,图21B还示出了使用在垂直虚线所指明的时间发生的时间分隔的完整易位步骤的组合或完整和分步易位步骤的组合所产生的示例性信号(粗实线)。时间分隔的完整易位步骤可以对应于分别易位多核苷酸,但在彼此相对改变(例如,彼此相对改变了完整易位循环持续时间的约50%)的时间的两个分子马达所产生的信号。完整和分步易位步骤的组合可以对应于通过部分和完整易位步骤分步易位多核苷酸的单个分子马达(例如,Hel308解旋酶)所产生的信号,如其中所述分步易位步骤在完整易位循环持续时间的约50%发生。所述信号可以参考图21A如上所述。在图21B中可以看出通过步骤序列信号从相对低电平改变为相对高电平,然后再通过另一步骤序列改变为相对低电平,其对应于多核苷酸易位通过孔时的时间分隔的完整易位步骤之间或通过完整和分步易位步骤的组合所发生的转换。在图21B中还可以看出信号在与图21A相比显著更多的点(和时间)与理想化信号相交,并因此可以认为在这些点对理想化信号“采样”。由于比图21A中的有效采样率相对更高,因此所述信号对理想化信号的采样比图21A中相对更好。例如,点(a)和(b)处的值彼此是相同的,其对应于不同易位步骤的简并信号电平。由于图21A中的信号也对点(a)和(b)之间存在的理想化曲线部分充分采样,因此对应于点(a)和(b)的物理易位步骤可以是彼此难以区分的,从而导致产生了多核苷酸序列的相关附加信息。另外,由于在图21B中信号对点(a)和(b)之间存在的理想化曲线部分更充分地采样,因此对应于点(b)和(c)之间理想化曲线的下降斜率的物理易位步骤可以比使用图21A中的信号可以潜在地实现的更好地表征这些步骤期间易位通过孔的多核苷酸部分。
除了如上所述由虚线曲线表示的理想化信号之外,图21C还示出了使用在垂直虚线所指明的时间发生的时间分隔的完整易位步骤的组合或完整和分步易位步骤的组合所产生的另一种示例性信号(粗实线)。时间分隔的完整易位步骤可以对应于分别易位多核苷酸,但在彼此相对改变(例如,彼此相对改变了完整易位循环持续时间的约25%、50%和75%)的时间的多个分子马达所产生的信号。完整和分步易位步骤的组合可以对应于通过部分和完整易位步骤分步易位多核苷酸的单个分子马达(例如,Hel308解旋酶)所产生的信号,如其中所述分步易位步骤在完整易位循环持续时间的约25%、50%和75%发生。所述信号可以参考图21A如上所述。在图21C中可以看出通过比图21B中更多个数的步骤序列信号从相对低电平改变为相对高电平,然后再通过比图21B中更多个数的步骤另一序列改变为相对低电平,其对应于多核苷酸易位通过孔时的时间分隔的完整易位步骤之间或通过完整和分步易位步骤的组合所发生的转换。在图21C中还可以看出信号在与图21B相比显著更多的点(和时间)与理想化信号相交,并因此可以认为在这些点对理想化信号“采样”。由于比图21B中的有效采样率相对更高,因此所述信号对理想化信号的采样比图21B中相对更好,并因此可以比使用图21A或21B中的信号潜在可以实现的更好地表征这些步骤期间易位通过孔的多核苷酸。
应理解可以使用任何适合的参数选择以提高理想化样品曲线的任何选择部分的采样。例如,如上所述,可以使用来自不同分子马达的时间-改变(相-改变)的完整易位步骤的组合。在这点上,尽管图21B描述了分子马达彼此之间时间改变了完整易位步骤时间的50%,并且尽管图21B描述了分子马达彼此之间时间改变了完整易位步骤时间的25%、50%和75%,但是这些值是纯说明性的,并且分子马达可以彼此之间时间-改变任何适合的时间量,例如,可以彼此之间改变完整易位步骤时间的5%至95%中的任一个量,例如,可以彼此之间改变完整易位步骤时间的10%至90%中的任一个量,例如,可以彼此之间改变完整易位步骤时间的25%至75%中的任一个量,例如,可以彼此之间改变完整易位步骤时间的40%至60%中的任一个量。作为另一实例,完整和分步易位步骤的组合可以对应于通过部分和完整易位步骤分步易位多核苷酸的单个分子马达(例如,Hel308解旋酶)所产生的信号。尽管图21B描述了彼此之间在完整易位步骤时间的50%发生的分步易位步骤,并且尽管图21C描述了彼此之间在完整易位步骤时间的25%、50%和75%发生的分步易位步骤,但是这些值是纯说明性的,并且分步易位步骤可以在相对于完整易位步骤的任何适合的时间发生,例如,在完整易位步骤时间的5%至95%,例如,在完整易位步骤时间的10%至90%,例如,在完整易位步骤时间的25%至75%,例如,在完整易位步骤时间的40%至60%。
另外,应理解可以通过改变任何适合的参数来适合地调整完整或分步步骤发生的相对时间和因此信号对理想化信号采样的时间。例如,如上所述。可以影响信号产生的示例性变量或参数可以包括温度、盐浓度(例如,Mg、Cl)、辅因子(例如,ATP)浓度、ATP产物如焦磷酸盐的浓度、pH、所使用的具体的分子马达等。在一些实施方式中,可以基于第一参数组产生第一信号,从而在第一离散时间组对理想化信号采样,并且可以基于第二参数组(其在至少一个方面不同于第一参数组)产生第二信号,从而在第二离散时间组对理想化信号采样。可以合并第一和第二信号以提供通过比单独的第一或第二信号更大的分辨率对理想化信号采样的信号曲线。应理解可以以类似的方式合并任何适合个数的信号以提供通过比那些信号中任何单独一个更大的分辨率对理想化信号采样的信号曲线。
实施例X
用于序列鉴别的其它方法
参考实施例X描述了用于序列鉴别的一些其它方法。
在一些实施方式中,可以单独或彼此结合地使用某些类型的信息以获得序列特异性信息:(A)单独的完整步骤响应信息、(B)单独的分步步骤响应信息、(C)不包括鉴别器的完整步骤和分步步骤响应信息,和(D)与鉴别器一起的完整步骤和分步步骤响应信息。
“响应信息”表示得自***对给定多核苷酸序列(k聚体)的数据,它对k聚体或k聚体的亚组(包括所关心的k聚体在内)是唯一的。响应信息的实例包括平均电流电平、中值电流电平、电流噪声的带宽电平、电流噪声的带限电平、电平持续时间等。
“鉴别器”表示多核苷酸(k聚体)与纳米孔环境相互作用所获得的数据,其中根据“理想响应”,特定电平相对于其它电平存在。例如,在存在相对高或相对低电平ATP浓度的情况下使用Hel308 Tga解旋酶作为分子马达的***可以分别对于每隔一个电平显示出相对短或者相对长的持续时间,其中每隔一个电平为根据来自相邻电平的理想响应的约50%。在本实施例中,电平持续时间可以用作鉴别因素,这是因为它可以用于鉴别沿理想响应的序列位置(相对于相邻层)。
“理想响应”表示***对通过纳米孔的特定多核苷酸易位的响应,其具有足够高的分辨率从而可以分辨多核苷酸足够小的运动。例如,理想响应是易位通过纳米孔的DNA的无限高的分辨率的连续电流示踪。
再次提及以上在本实施例中进一步提到的项目(A)-(D),可以独立使用项目(A)-(D)中的每一个或与项目(A)-(D)中一个或多个其它项目结合使用以鉴别多核苷酸序列。例如,由于(例如)计算资源限制、时间限制、优化方法的先验知识等,可以根据项目(A)-(D)中的任何其它项目独立计算项目(A)-(D)中的一个或多个。基于要计算的项目(A)-(D)中不止一个项目,可以仅使用项目(A)-(D)中的一个项目的结果。对于要使用该类计算中哪一个的确定可以基于结果的可信程度。例如,结果的可信性可以基于以下中的一种或多种:(a)响应信息本身(例如,高电平的ATP可以缩短Hel308 Tga中的分布步骤长度,其可以降低项目(B)相对于项目(A)的可信性);(b)测序算法本身(例如,Viterbi算法可以对它所提出的最佳序列产生可能性得分,其可以用于确定所提出的序列的可信性水平);(c)测序算法所产生的序列(例如,可以根据算法所提出的序列与序列对照表和/或要测序的多核苷酸的任何先验知识之间的比较来分配可信性);或者(d)项目(a)-(c)的任何适合的组合。
注意在一些情况下,通过使用根据项目(A)-(D)中不止一个项目所提出的序列来确定真实序列可以是有益的。例如,可以基于一些或所有这些所提出的序列确定共有序列。可以基于所有或一些所提出的序列来确定共有序列。可以将共有序列全局应用于整个多核苷酸序列或者局部应用于序列的一部分。可以基于来自一些或所有项目(A)-(D)的可信性值来确定共有序列。可信性值可以是以上在本实施例中进一步描述的那些。可信性值可以局部应用于序列的一部分,或者全局应用于整个序列。可以通过多轮上述方法来确定最终的共有序列,其中可以将每轮所获得的共有序列用作所提议的序列,并且每轮之间用于确定每轮的可信性的方法可以是不同的。
作为一个实例,Viterbi算法可以用于通过仅对易位通过纳米孔的DNA的完整步骤或仅对分步步骤测序来确定两个不同的所提议的序列(项目(A)和(B),如本实施例中上述的)。将该算法对每个DNA片段的可能性得分用于确定所提议的序列的每个区域的置信度,并且每个区域的置信度的集合体可以导致产生第一轮提议的共有序列。然后,可以就已知序列的对照表将该共有序列与两个初始提议的序列相比较。对照表和这三个提议的序列之间的相似性可以导致对三个提议序列中每个序列的每个区域产生置信值。这三个提议的序列之间的第二轮基于置信度的比较可以导致产生最终提议的共有序列。
在一些实施方式中,图22A-22D示出了根据一些实施方式使用多核苷酸通过孔的分步易位所提供的信息的说明性方法中的步骤。图22A示出了根据一些实施方式使用多核苷酸通过孔的分步易位所提供的信息的方法的高电平概述。图22A中所示的方法包括获得信号(步骤2210),如通过Hel308解旋酶的目标多核苷酸通过孔的一个或多个分步易位步骤所产生的一个或多个信号,如本文其它处更详细地描述的。图22A中所示的方法还包括电平检测和识别(步骤2220),例如,检测和鉴别信号中不同的信号电平,例如,检测和鉴别对应于多核苷酸通过孔的分步易位步骤的电平,以及检测和鉴别对应于多核苷酸通过孔的完整易位步骤的电平。图22A中所示的方法还包括序列确定(步骤2230),例如,基于检测和鉴别信号中不同的信号电平来表征多核苷酸的序列。图22A所示的方法还包括输出序列(步骤2240),例如,基于序列调用结果输出实际核苷酸的可能核苷酸序列。
图22B-22D示出了图22A中所示的一个或多个步骤的任选的子步骤。例如,图22B示出了图22A所示的步骤2210和2220的一种可能实施的额外细节。图22B中所示的方法还包括获得信号(步骤2210),如通过Hel308解旋酶的目标多核苷酸通过孔的一个或多个分步易位步骤所产生的一个或多个信号,如本文其它处更详细地描述的。图22B所示的方法还可以任选地包括获得输入参数(2211)。这些输入参数可以包括(但不限于)定义应检测和确定以对应于信号的那些特征信号的特征的参数。例如,输入参数可以定义信号值中的阈值幅度变化,可以将高于该信号幅度变化检测为对应于电平。或者,例如,输入参数可以限定应检测仅对应于完整易位步骤的信号电平,或仅对应于分步易位步骤的信号电平,或者对应于完整和分步易位步骤两者的信号电平。输入参数还可以含有与错误形式(例如,核苷酸跳跃或核苷酸切换)有关的信息,其可能包括某些错误形式的倾向和/或程度,当确定电平时,可以对其加以考虑。输入参数还可以包括与特定环境有关的信息,在此之间纳米孔、分子马达和多核苷酸相互作用(例如,温度、盐度、pH、辅因子浓度等),这可以用于确定给定信号的电平。图22B中所示的方法还包括电平检测,例如,检测对应于多核苷酸通过孔的分步易位步骤的信号中不同的信号电平(步骤2221)。例如,基于步骤2210获得的信号和步骤2211获得的输入参数,该电平检测可以检测充分统计学显著不同于其它信号区域的信号区域以对应于电平。电平检测的示例性方法(其还可以称为边缘检测或步骤检测)在本领域中是已知的,并且包括学生t检验和卡方最大化。对于可以适合地调整以用于在步骤2221的检测电平中使用的步骤检测算法的一些实例,参见Carter等人,“A Comparison of Step-DetectionMethods:How Well Can You Do?,”Biophysical Journal 94:306-308(January 2008)。
图22B中所示的方法还包括基于步骤2221的电平检测的输出电平信息(步骤2222)。电平信息可以包括对于给定电平检测的平均、中值、形式、分步、持续时间、最大和/或最小电流,或者这些值的任意组合,或者属于给定电平的电流值的亚组的这些值(例如,在首先除去与错误形式有关的电流信息后,可以使用平均电流)。电平信息还可以包括电流的标准偏差或电流的频带限制亚组(例如,在使用低通、高通、带通或带阻滤波器或这些滤波器的任何组合后,获得电流)。电平信息还可以包括与电平持续时间有关的信息以及与电平有关的错误形式信息。图22B中所示方法还包括电平识别(步骤2223),例如,确定在步骤2222输出的电平信息的在步骤2221检测的哪个电平对应于目标多核苷酸的完整或分步易位步骤。例如,步骤2223可以包括分析在步骤2222输出电平信息的在步骤2221检测的不同电平的持续时间,并且基于这些持续时间,鉴别对应于完整易位步骤的某些电平,并且鉴别对应于分步易位步骤的其它某些电平。作为一个实例,可以认为具有比第一阈值更短的持续时间的信号电平对应于噪音并因此弃去,而可以认为具有比第一阈值更长但比第二阈值更短的持续时间的信号电平对应于分步易位步骤并因此鉴别为该信号,而可以认为具有比第二阈值更长但比第三阈值更短的持续时间的信号电平对应于完整易位步骤并因此鉴别为该信号,而可以认为具有比第三阈值更长的持续时间的信号电平对应于错误或不存在多核苷酸并且因此弃去。
图22B中所示的方法还包括输出以下输出中的一种或多种:全电平、分电平、全部电平和电平识别符(level identifiers)。例如,如上所提及的,在步骤2211获得的输入参数可以限定应检测仅对应于完整易位步骤的信号电平,或仅对应于分步易位步骤的信号电平,或者对应于完整和分步易位步骤两者的信号电平(例如,“全部电平”)。注意,在一些实施方式中,通过输入参数选择“全部电平”可以对应于绕过电平识别步骤,从而电平检测步骤2221直接输出全部电平。作为另外一种选择,基于电平检测结果2223和输入参数2211,可以输出所期望的信号的鉴别电平,例如,用于进一步加工,如以下参考图22C和22D所述的。电平识别符可以包括有利于电平进一步分析的任何适合的信息,例如,表明步骤2223期间使用的完整或分步步骤的持续时间以表示所识别的电平对应的转换类型的指数。
再参考图22A,可以使用图22B中所示方法或者使用另一种适合的方法产生的全电平、分电平、全部电平和电平鉴别符中的一个或多个可以用作进行序列确定的输入(图22A中的步骤2230)。例如,图22C示出了基于这些全电平、分电平、全部电平和电平识别符中的一个或多个,例如,作为输入的全电平、分电平、全部电平和电平识别符中的一个或多个进行序列确定的第一示例性方法。图22C中所示的方法包括基于全电平、分电平、全部电平和电平识别符中的一个或多个的输入的序列调用步骤(步骤2231)。序列调用可以包括基于输入信号电平的基于目标多核苷酸的哪个核苷酸碱基可以调用的任何适合的方法。序列调用的示例性方法包括(但不限于)Viterbi算法,如实施例VII中参考图19A所述,改进型Viterbi算法,如实施例VII中参考图19B所述,或类似于实施例XI中所述的图形匹配。可以适合地使用序列调用的其它方法。序列调用的输出(步骤2231)可以包括多种调用序列,例如,序列A、序列B、...序列N,以及对于每个该调用序列的置信度信息。不同的调用序列可以基于彼此不同的对步骤2231的输入。例如,第一调用序列(例如,序列A)可以基于对步骤2231的输入,其中基于在步骤2210获得的给定信号,仅鉴别完整易位电平,第二调用序列(例如,序列B)可以基于对步骤2231的输入,其中仅鉴别分步易位电平,和第三调用序列(例如,序列N)可以基于对步骤2231的输入,其中鉴别所有易位电平(例如,完整和分步易位电平两者)。作为另外一种选择或另外,其它调用序列可以基于其它电平,所述其它电平是基于在步骤2211获得的替代输入参数鉴别的,如限定应检测和确定哪种特征信号的特征以对应于信号的参数的不同值,如信号值的不同阈值幅度变化,可以将高于此的信号幅度变化检测为对应于电平。每个不同的调用序列可以具有相关置信度信息,例如,代表调用序列对应于目标核苷酸真实序列的可能性的值。
在图22C所示的实施方式中,序列选择步骤(步骤2232)可以选择一个或多个调用序列并作为输出提供选择的序列(步骤2240)。作为一个实例,序列选择步骤(步骤2232)可以包括对多个调用序列比较置信度信息,并且可以在步骤2240选择并输出具有高置信度的调用序列,例如,对应于真实序列的最高可能性。作为另一个实例,给定调用序列的置信度信息可以包括分别代表调用序列的相应部分对应于目标多核苷酸的该部分的真实序列的可能性的多个置信度值。对于调用序列的不同部分(例如,10个碱基对长,或50个碱基对长,或100个碱基对长,或10-100个碱基对长,或10-50个碱基对长,或50-100个碱基对长的部分),序列选择步骤(步骤2232)可以包括不同调用序列在该部分的置信度值的比较,以及选择对该部分具有最高值的调用序列的部分。所选部分可以与对这些部分具有最高值的其它调用序列的所选部分连结或进行比对。
图22D示出了可以用于序列确定的替代方法(2230)。图22D中所示的方法可以包括作为输入获得多个调用序列,例如,序列A、序列B、...序列N,以及对于每个该调用序列的置信度信息,其可以类似于参考图22C的上述那些。在这点上,尽管未具体说明,但是图22D中所示的方法可以包括序列调用步骤2231:类似于参考图22C的上述那些的接收输入,类似于参考图22C的上述那些的提供输出和与步骤2231类似的操作。作为另外一种选择,图22D中所示的方法可以从任何其它适合来源获得多个调用序列。
图22D所示的方法还可以包括获得模型序列(步骤2234)。例如,这些序列可以包括对于一个或多个不同物种,如一个或多个不同的病原体,先前已知的序列。说明性地,可以将模型序列保存在查找表、数据库或保存在永久计算机可读介质上的其它适合的数据结构中。图22D所示的方法还可以包括序列选择步骤(步骤2233)。在图22D所示的实施方式中,序列选择步骤可以基于在步骤2234获得的一个或多个模型序列,选择作为输入接受的一个或多个调用序列,并作为输出提供所提议的序列和新的置信度信息。作为一个实例,序列选择步骤(图22D中的步骤2233)可以包括将多个调用序列中的一个或多个与在步骤2234获得的一个或多个模型序列相比较,并且可以选择和输出所提议的序列,其可以对应于具有最高的新的置信度信息的调用序列,例如,对应于模型序列的最高可能性。输入置信度信息可以对序列(或者该序列内的区域)的可能性至关重要,以匹配模型序列(模型序列内的区域)以确定最可能的序列,并将其作为提议的序列输出。例如,对于与模型序列Z最佳比对的输入序列A和与模型序列Y最佳比对的输入序列B,基于A和Z之间的比对优于B和Y之间的比对,所提议的序列可以是模型序列Z。然而,在B和Y具有最佳比对的情况下,则A比B具有更高置信度值的情况可以使Z成为所提议的序列。另外,在另一种情况中,可以比较序列区域,从而使输出的提议序列包含来自A、B、Z和Y的序列信息。作为另外一种选择,给定调用序列可以包括新的置信度信息,例如,分别代表调用序列的相应部分对于该部分对应于一个或多个模型序列的部分的可能性的多个新的置信度值。对于调用序列的不同部分(例如,10个碱基对长,或50个碱基对长,或100个碱基对长,或10-100个碱基对长,或10-50个碱基对长,或50-100个碱基对长的部分),序列选择步骤(步骤2233)可以包括将不同调用序列在该部分的新的置信度值与模型序列进行比较,和选择对该部分具有最高的新的置信度值的调用序列的部分。所选部分可以与对这些部分具有最高的新的置信度值的其它调用序列的所选部分连结或进行比对。
图22D中所示的方法还可以包括基于步骤2233的新的置信度信息输出,确定对于所提议的序列新的置信度信息以及步骤2233的输出是否满足要求(步骤2235)。作为一个实例,步骤2235可以将新的置信度信息(其可以是新的置信度值)与阈值置信度值相比较,在该阈值置信度值或之上,所提议的序列可以确定为与模型充分匹配,而低于阈值置信度值,所提议的序列可以确定为与模型未充分匹配。新的置信度信息可以包括输入置信度信息的结果,所提议的序列和输入序列之间的关系,所提议的序列和模型序列之间的关系和/或输入序列和模型序列之间的关系。例如,在所提议的序列只不过是输入序列之一的情况下,新的置信度信息可以是输入序列的输入置信度值的加权平均数以及与最佳比对模型序列的比对得分。在其它情况下,如当所提议的序列是输入序列区域的组合时,新的置信度信息可以包括输入可信度值的加权平均值的加权平均值以及所提议序列的区域中(与模型序列)的比对得分。基于在步骤2235确定新的置信度信息满足要求(“是”),步骤2235作为输出提供所提议的序列(步骤2240)。基于在步骤2235确定新的置信度信息不满足要求,步骤2235返回至2233,此时继续序列选择,例如,通过调用序列与模型序列的进一步比较。序列选择算法或模型序列组可以依赖于参数,其可以包括所提议的序列、新的置信度信息、实施序列选择算法的次数和已使用的模型序列中的一个或多个。例如,初始通过序列选择算法可以使用相对少的模型序列(例如,为了通量)。然而,如果输入序列和模型序列之间的比对相对较差,则新的置信度信息可能不满足要求,并因此一旦返回步骤2233,则可以与新的一个或多个模型序列的精炼组进行比较。
实施例XI
图形识别,任选地用于SNP鉴别
在一些实施方式中,本文所公开的方法和组合物可以与多路核酸检测、基因分型和扩增方法组合使用。用于多路核酸检测、基因分型和扩增的方法在本领域中是熟知的并且技术人员可以容易地选择并使用上述方法。例如,在一个实施方式中,本文所公开的方法和组合物可以与以下美国专利中所述的多路核酸检测、基因分型和扩增方法组合使用:美国专利No.6,890,741、6,913,884、7,955,794、7582,420和8,288,103和美国专利公开2013-0244882,以上专利作为参考并入本文。
在一些实施方式中,可以与本文所公开的方法和组合物组合的用于多路核酸检测、基因分型和扩增的方法包括在固体载体上实施或与固体载体组合实施的方法,所述固体载体如阵列(随机和有序阵列)或珠。例如,在一些方面,可以将要测定的目标多核苷酸,如基因组DNA固定至固体载体。可以对这些固定的目标多核苷酸进行在本领域中熟知的多路核酸检测和基因分型方法。可以使用本文所公开的方法表征所得目标多核苷酸。
在一些实施方式中,用于表征目标多核苷酸的方法还可以包括产生要测定的目标多核苷酸所必需的步骤。因此,在一些实施方式中,所述方法可以包括以下步骤:(a)提供从3′至5′分别包含第一、第二和第三目标结构域的多个目标核酸序列,其中第一目标结构域包含检测位置、第二目标结构域为至少一个核苷酸;(b)将目标核酸序列与每个目标序列的探针组接触以形成一组第一杂交复合物,每个探针组包括:第一探针,其从5′至3′包含通用引物序列和与目标序列的第一目标结构域基本互补的序列,以及适合于和检测位置碱基配对的探查位置(interrogation position)(例如,在3′端4个末端碱基内),和第二探针,其从5′至3′包含与目标序列的第三目标结构域基本互补的序列和通用引物序列,其中任选地至少一个探针含有位点鉴别序列(例如,标签或条型码);(c)在一定情况下将杂交复合物与延伸酶和dNTP接触,其中如果探查位置处的碱基与检测位置的碱基完美互补,则第一探针发生延伸通过第二目标结构域以形成第二杂交复合物;和(d)将延伸的第一探针与第二探针连接以形成扩增模板。在该方法的一些方面,探针组的第一或第二探针可以包括等位基因鉴别序列(例如,标签或条型码)。
在一些实施方式中,用于表征目标多核苷酸的方法还可以包括以下步骤:(a)提供从3′至5′分别包含第一、第二和第三目标结构域的多个目标核酸序列,其中第一目标结构域包含检测位置、第二目标结构域为至少一个核苷酸;(b)将目标核酸序列与探针接触,所述探针分别从5′至3′包含通用引发序列和与目标序列的第一目标结构域基本互补的序列,以及适合于和检测位置碱基配对的探查位置(例如,在3′端4个末端碱基内),其中任选地所述探针含有位点鉴别序列(例如,标签或条型码);(c)在一定情况下将杂交复合物与延伸酶和dNTP接触,其中如果探查位置处的碱基与检测位置的碱基完美互补,则探针发生延伸通过第二和第三目标结构域以形成可以用作扩增模板的延伸探针。
产生用于在本文所述的方法中测定的目标多核苷酸的方法还可以包括将扩增模板扩增以产生扩增子。在一些方面,基于已引入扩增模板的鉴别序列,包含第一或第二探针的通用引发序列的引物还包含等位基因鉴别序列或位点鉴别序列(例如,标签或条型码)。可以使用本文所公开的方法表征这些扩增子,所述扩增子可以包含位点鉴别序列和等位基因鉴别序列两者。目标序列的表征可以表明样品的基因分型基于位点和等位基因鉴别序列的存在。
在一些实施方式中,用于产生扩增子的引物包括不允许延伸酶在扩增期间使用以穿过残基的一个或多个修饰的残基。例如,在一些方面,一个引物包括无碱基位点(无嘌呤/无嘧啶位点)、C3间隔臂亚磷酰胺(Int C3间隔臂)、三甘醇间隔臂(Int间隔臂9)或18原子的六乙二醇间隔臂(Int间隔臂18)以防止延伸酶继续引物延伸。应理解本领域的技术人员可以选择可以实施这种相同功能的其它修饰的残基。一个或多个修饰的残基可以位于等位基因鉴别序列内或等位基因鉴别序列的任一侧,只要对于使用本文所公开的方法的目标多核苷酸的表征产生了足够长度的5′突出。例如,5′突出具有足够长度以允许扩增子的固定化。
在一些实施方式中,通过上述方法产生的扩增子还与切口核酸内切酶接触,从而在第二探针序列中或附近产生3′突出。这些切口酶可以是序列特异性的,从而仅切割双链产物中的一条链。多种切口核酸内切酶在本领域中是熟知的并且认识到基于探针和引发序列,技术人员可以容易地选择适合的核酸内切酶。为了在切口核酸内切酶切割后产生3′突出,可以使用本领域中已知的一些方法,包括(例如)扩增子的部分变性,从而从扩增子释放较小部分的切口链,而扩增子的剩余部分仍保持杂交在一起。为了辅助除去扩增子的较小部分,可以加入所述较小部分的反向互补序列以杂交至不希望的链。
在一些实施方式中,可以通过在以上方法中所述的第二探针序列中包含一个或多个尿嘧啶残基来产生3′突出,并将扩增子与在尿嘧啶位置处特异地产生单个核苷酸缺口的尿嘧啶特异的酶接触。这类尿嘧啶特异的酶的非限制性实例是Uracil-Specific ExcisionReagent(USERTM)Enzyme(New England Biolabs)。因此,使用熟知的方法,可以容易地将通过所述酶产生的较小交替片段与扩增子变性开。
在某些方面,所产生的3′突出具有足够长度以有利于本文所述的解旋酶的结合。因此,在一些方面,3′突出包括至少4个核苷酸长度。在其它方面,3′突出包括4-20个核苷酸长度,或在某些方面,8-16个核苷酸长度,或在其它方面,包括10-16个核苷酸长度。
短语“位点鉴别序列”是指已分配或已知连接至目标多核苷酸上特定位点的核酸残基序列(例如,标签或条型码)。目标多核苷酸的位点可以是(例如)基因组上接近于要测定的等位基因的基因、基因的一部分(例如,外显子或内含子)或非编码区(例如,启动子或增强子)。位点鉴别序列可以是对所关心的目标序列的位点特异的天然存在的序列和/或对所关心的目标序列非天然的合成序列。可以通过标签或条型码的预期信号图案分配位点鉴别序列。
短语“等位基因鉴别序列”是指已分配给目标多核苷酸检测位置中特定核酸残基的核酸残基序列(例如,标签或条型码)。等位基因鉴别序列可以表明检测位置中核酸残基(例如,A、T、C或G)的存在。还可以通过标签或条型码的预期信号图案来分配等位基因鉴别序列。
在另一个实施方式中,用于表征目标多核苷酸的方法还可以包括图18中所述的步骤。该方法可以包括以下步骤:(a)提供具有不同的所关心的目标核酸序列的样品,其中所述不同的目标核酸序列任选地固定在固体载体上;(b)将样品与每个不同的所关心的目标核酸序列的探针组接触以形成杂交复合物,每个探针组包括:第一探针,其从5′至3′包含第一通用引发序列和与第一目标结构域基本互补并且具有适合于和检测位置碱基配对的探查位置(interrogation position)的序列;和第二探针,其从5′至3′包含与第三目标结构域基本互补的序列和第二通用引物序列,其中至少一个探针含有对所关心的目标序列非天然的位点鉴别序列(例如,标签或条型码);(c)将杂交复合物与延伸酶和dNTP接触,其中对于每个杂交复合物,如果探查位置处的碱基与检测位置的碱基完美互补,则第一探针沿第二目标结构域延伸;(d)将延伸的第一探针与第二探针连接以形成扩增模板;(e)用第一和第二通用引物扩增所述扩增模板以产生扩增子,其中至少一个引物包括等位基因鉴别序列(例如,标签或条形码),其中等位基因鉴别序列包括无碱基位点;(f)将扩增子与切口核酸内切酶接触,从而在第二引物序列中产生3′突出;和(g)使用本文所述的目标多核苷酸表征方法,检测不同扩增子的位点鉴别序列和等位基因鉴别序列的存在,借此表明样品中所关心的不同靶标序列的存在。
如本文所使用的,短语“多路”或语法等价形式是指不止一个所关心的目标序列的检测、分析或扩增。在一个实施方式中,多路是指至少100或200个不同的目标序列,同时至少500个不同的靶标序列是优选的。更优选地,至少1000个,特别优选地,5000或10,000个,并且最优选地,大于50,000或100,000个。检测可以在如本文所述的多个平台上进行。
在一些方面,本文的发明公开提供了检测样品中核酸目标序列的方法。如本领域技术人员将理解的,样品溶液可以包含多种东西,其包括(但不限于)体液(包括(但不限于)实际上任何生物的血液、尿液、血清、淋巴、唾液、***和***分泌物、汗液和***,并且哺乳动物样品是优选的,并且人样品是特别优选的);环境样品(包括(但不限于)空气、农业水和土壤样品);生物战剂样品;研究样品;纯化样品,如纯化的基因组DNA、RNA、蛋白质等;粗样品(细菌、病毒、基因组DNA等)。如本领域技术人员将理解的,实际上可以对样品进行任何实验操作。
如果需要,使用已知技术制备目标多核苷酸。例如,可以使用已知的裂解缓冲液、超声、电穿孔等处理样品使细胞裂解,并根据需要,如以下所列进行纯化和扩增,如本领域技术人员将理解的。另外,可以以多种方式进行本文所列的反应,如本领域技术人员将理解的。可以以任何顺序同时或顺序添加反应组分,其中下文列出了优选实施方式。另外,反应可以包括在测定中可以包含的多种其它试剂。这些包括试剂,如盐、缓冲剂、中性蛋白质,例如白蛋白、去污剂等,其可以用于辅助最优化杂交和检测,和/或降低非特异性或背景相互作用。另外,根据样品制备方法和目标纯度,可以使用改善测定效率的试剂,如蛋白酶抑制剂、核酸酶抑制剂、抗微生物剂等。
另外,在大部分实施方式中,将双链目标多核苷酸变性以使它们成为单链,从而允许引物和本文所述的其它探针杂交。一个实施方式使用了热步骤,通常通过将反应温度升高至约95℃,尽管还可以使用pH变化及其它技术。
如本文所列,目标多核苷酸可以是反应产物,如来自反应的检测序列、连接的探针、来自PCR反应的延伸探针或PCR扩增产物(“扩增子”)等。
在一些实施方式中,目标多核苷酸包括需要序列信息的位置,其在本文中一般地称为“检测位置”。在某一实施方式中,检测位置是单个核苷酸,尽管在一些实施方式中,它可以包括彼此邻接或被一个或多个核苷酸分开的多个核苷酸。如本文所使用的“多个”表示至少两个。如本文所使用的,杂交中与检测位置碱基碱基配对的碱基被称为“读取位置”或“审查位置”;因此,本发明的第一或第二步骤探针中的多个包含审查位置。
本文所公开的方法可以采取多种构型,如附图中所示和如本文中更详细描述的。通常,这些部分包括复杂性减小部分(complexity reduction component)、特异性部分和扩增部分。可以以如下所公开的多种方式配置部分。也就是说,在一个实施方式中,首先执行复杂性减小步骤。随后是扩增或特异性步骤。作为另外一种选择,首先执行特异性步骤。随后可以是复杂性减小或扩增步骤。作为另外一种选择,首先执行扩增。随后是复杂性和特异性步骤。
然而,以上表明可以以任何顺序执行三个部分中的每一个。本领域技术人员将理解当首先执行扩增时,将有可能涉及一定程度的复杂性减小或特异性。另外,当首先执行特异性部分时,将有一定程度的复杂性减小。另外,在一些实施方式中,当首先执行扩增时,将有一定程度的特异性和复杂性减小。然而,如下所述,所述方法通常包括三个部分。
探针和引物
如本领域技术人员所理解的,存在一些可以在本文所公开的方法中使用的探针或引物。这些探针/引物可以采取多种构型,并且可以具有以下更详细描述的多种结构成分。第一步骤探针可以是等位基因特异性探针或位点特异性探针。“等位基因特异性”探针或引物是指与目标序列杂交并且区分等位基因或与目标序列杂交并且以等位基因特异的方式修饰的探针或引物。“位点特异性”探针或引物是指以位点特异性方式杂交至目标序列,但不必需区分等位基因的探针或引物。如下所述,还可以修饰(即延伸)位点特异性引物,从而它包含了特定等位基因相关的信息,但是位点特异性引物不区分等位基因。
在多种实施方式中,探针或引物包含一个或多个通用引发位点和/或鉴别序列。例如,在一种构型中,四个等位基因碱基中的每一个与不同序列(即等位基因鉴别序列(例如,标签或条型码))结合,每个序列具有类似的扩增效率。在另一个构型中,探针之一包含位点鉴别序列(例如,标签或条型码)。
引物和探针核酸的大小可以不同,如本领域技术人员将理解的,探针的每个部分和探针的总长度一般在5至500个核苷酸长度改变。根据用途和扩增技术,每个部分可以在10至300、15至250或10至35个核苷酸长度之间。因此,例如,探针的通用引发位点可以在15-20个核苷酸长度之间,在某些实施方式中,使用了18个核苷酸长度。探针的位点和/或等位基因鉴别序列可以在10-300个核苷酸长度之间,在某些实施方式中,使用了20-100个核苷酸长度。探针的目标特异性部分可以在15-50个核苷酸长度。另外,引物可以包含额外的扩增引发位点。
在一个实施方式中,等位基因或位点特异性探针包含与目标序列的第一结构域基本互补的目标结构域。一般地,可以将探针设计成与目标序列(样品的目标序列或其它探针序列,如本文所述)互补,从而发生目标和本文所述的探针的杂交。这种互补性不必需是完全的;可以存在多个防碍本发明的目标序列和单链核酸之间杂交的碱基对错配。然而,如果突变数目过大从而即使在最不严格的杂交条件下也不能发生杂交时,所述序列与目标序列不互补。因此,本文中“基本互补”表示探针与目标序列足够互补,从而在所选反应条件下杂交。
另外,可以构建在本文所述的方法中使用的探针以含有必需的引发位点或用于后续扩增方案的位点。在某些实施方式中,所述引发位点是通用引发位点。本文中“通用引发位点”或“通用引发序列”是指将结合用于扩增的引物的探针序列。
如本领域技术人员将理解的,一般地,可以实施高度多路反应,对于所有反应,所有通用引发位点相同。作为另外一种选择,可以同时或顺序使用通用引发位点“组”和相应探针。通用引发位点用于扩增修饰的探针以形成多个扩增子,然后以多种方式检测扩增子,如本文所列。
因此,本文所述的方法提供了第一目标探针组。本文中“探针组”是指在特定多路测定中使用的多个目标探针。在本文中,多个表示至少2个,其中根据测试的测定、样品和目的,优选大于10个。在一个实施方式中,探针组包括大于100个探针,优选大于500个探针,并且更优选大于1000个探针。在特别优选的实施方式中,每个探针组含有至少5000个探针,最优选大于10,000个探针。
复杂性减小部分
复杂性减小可以是本文所述的多路方案的部分。通常,复杂性减小是富集特定目标或位点的方法。也就是说,认为复杂性减小是导致非目标核酸从样品中移除或未正确杂交或根本未杂交至目标核酸的探针/引物的移除的方法。另外,复杂性减小包括酶促步骤期间未修饰的探针的移除。也就是说,复杂性减小包括在酶促步骤(即扩增或特异性步骤或两者)之前移除非目标核酸,即富集目标核酸,或移除未杂交的探针或引物。
存在多个包含复杂性减小步骤的方法。这些包括(但不限于)以目标特异性方式修饰的目标核酸或探针/引物的选择性固定化,非目标核酸的选择性移除和非目标核酸的选择性破坏。这种破坏包括(但不限于)非目标核酸的变性、降解或切割。另外,复杂性减小可以包括如目标选择性扩增的部分,尽管这还包括扩增和部分。
在某些实施方式中,复杂性减小伴随着以目标特异性方式选择性固定修饰的引物。也就是说,位点特异性或等位基因特异性引物与目标杂交。目标可以固定化或处于溶液中。杂交后,引物在引物延伸反应中延伸。在一些方面,引物或NTP包含纯化标签,其使得能够从反应混合物中移除或纯化延伸产物。一旦延伸,通常可以将修饰的引物固定在固体载体上。修饰的引物固定化后,可以清洗载体以除去未修饰(即延伸)的非目标核酸和引物两者。因此,固定化的引物包括目标位点相关信息,其包括特定等位基因信息。这导致目标核酸的富集或非目标核酸的移除。
在另一个实施方式中,复杂性减小部分包括目标多核苷酸的选择性固定化。也就是说,优先地将目标多核苷酸而不是非目标核酸固定在固体载体上。
在一个实施方式中,将目标多核苷酸、探针或引物,包括修饰的引物连接至固体载体。本文中“固体载体”或其它语法等价形式是指适合于或者可以修饰以适合于与目标序列连接的任何材料。如本领域技术人员将理解的,可能的底物数目非常大。可能的底物包括(但不限于)玻璃和修饰或官能化的玻璃、塑料(包括丙烯酸类、聚苯乙烯以及苯乙烯共聚物及其它材料、聚丙烯、聚乙烯、聚丁烯、聚氨脂、TeflonTM等)、多糖、尼龙或硝化纤维素、陶瓷、树脂、二氧化硅或二氧化硅基材料,包括硅和修饰的硅、碳、金属、无机玻璃、塑料、光纤束和多种其它聚合物。磁珠和高通量微滴板是特别优选的。
固体载体的组成和几何形状根据其用途而改变。在某些实施方式中,包括微球或珠的载体可以用于固体载体。本文中“微球”或“珠”或语法等价形式表示小的分离颗粒。根据生物活性剂的种类和合成方法,珠的成分将改变。适合的珠成分包括在肽、核酸和有机部分合成中使用的那些,其包括(但不限于)塑料、陶瓷、玻璃、聚苯乙烯、甲基苯乙烯、丙烯酸聚合物、顺磁性物质、二氧化钍溶胶、石墨化碳、二氧化钛、乳胶或交联右旋糖酐,如琼脂糖、纤维素、尼龙、交联胶束和聚四氟乙烯以及本文所列的用于固体载体的任何其它材料,所有材料均可以使用。来自Bangs Laboratories,Fishers IN的“Microsphere DetectionGuide”是有帮助的指南。优选地,在该实施方式中,当执行复杂性减小时,所述微球是磁性微球或珠。
一旦连接到固体载体,对目标序列、探针或引物进行如本文所述的分析。
在本发明中可以使用多种杂交或清洗条件,其包括高、中和低严格性条件;参见,例如Maniatis等人,Molecular Cloning:A Laboratory Manual,第2版,1989和ShortProtocols in Molecular Biology,Ausubel等人主编,以上文献作为参考并入本文。严格性条件是序列依赖性的并且在不同情况下将是不同的。较长的序列在较高的温度下特异性杂交。核酸杂交的广泛指南可见于Tijssen,Techniques in Biochemistry and MolecularBiology-Hybridization with Nucleic Acid Probes,“Overview of principles ofhybridization and the strategy of nucleic acid assays”(1993)。通常,严格性条件选择为比限定离子强度和pH下特定序列的热熔点(Tm)低约5-10℃。Tm是(限定离子强度、pH和核酸浓度下)50%与目标互补的探针杂交至目标序列达到平衡时的温度(由于在Tm,目标序列过量存在,因此平衡时50%的探针被占据)。严格性条件将是其中pH 7.0至8.3,盐浓度小于约1.0M钠离子,通常约0.01至1.0M钠离子浓度(或其它盐),并且对于短探针(例如,10至50个核苷酸)温度至少约30℃,对于长探针(例如,大于50个核苷酸),至少约60℃的那些条件。还可以通过添加螺旋去稳定剂,如甲酰胺来实现严格性条件。
本文中“延伸酶”是指通过加入NTP使序列延伸的酶。如本领域中熟知的,存在多种适合的延伸酶,其中聚合酶(RNA和DNA两者,基于目标序列组成和预环化探针(precircleprobe))是优选。优选的聚合酶是缺乏链置换活性,从而它们将能够仅在探针末端添加必要的碱基,而不会进一步延伸探针以包含与目标结构域互补的核苷酸并因此防止环化的那些。适合的聚合酶包括(但不限于)DNA和RNA聚合酶两者,其包括DNA聚合酶I的克列诺片段、SEQUENASE 1.0和SEQUENASE 2.0(U.S.Biochemical)、T5DNA聚合酶、Phi29 DNA聚合酶和多种RNA聚合酶,如来自栖热菌属(Thermus sp.)的RNA聚合酶,或来自噬菌体的Qβ复制酶,以及SP6、T3、T4和T7RNA聚合酶等均可使用。
聚合酶还可以包括本质上缺少5′至3′核酸外切酶活性以确保探针将不会延伸通过探针的5′端的那些。缺少5′至3′核酸外切酶活性的示例性酶包括DNA聚合酶的克列诺片段和DNAPTaq聚合酶的Stoffel片段。例如,Taq DNA聚合酶的Stoffel片段由于基因操作而缺少5′至3′核酸外切酶活性,这导致缺少N末端289个氨基酸的截短蛋白的产生。(参见,例如,Lawyer等人,J.Biol.Chem.,264:6427-6437(1989);和Lawyer等人,PCR Meth.Appl.,2:275-287(1993))。对于来源于海栖热袍菌(T.maritima)的聚合酶、Tsps17、TZ05、Tth和Taf,已产生了类似的突变体聚合酶。
其它聚合酶是缺少3′至5′核酸外切酶活性的那些,这种活性通常被称为校对活性并且它移除了引物模板双螺旋3′端的错配碱基。尽管3′至5′核酸外切酶活性的存在提供了链合成中提高的保真度,但是热稳定的DNA聚合酶(如Tma,包括缺少5′至3′核酸外切酶活性的Tma突变形式)中存在的3′至5′核酸外切酶活性也降解单链DNA,如PCR中使用的引物、单链模板和单链PCR产物。引物延伸过程中使用的寡核苷酸引物的3′端的完整性至关重要,这只因为正是从该末端开始新生链的延伸。3′端的降解导致寡核苷酸缩短,这反过来导致引发反应中特异性的丧失(即引物越短,则越有可能发生伪引发或非特异性引发)。
其它聚合酶是热稳定性聚合酶。耐热酶可以包括在最优条件下,在40℃保持1小时后仍保留其大部分活性的任何酶。缺少5′至3′核酸外切酶和3′至5′核酸外切酶的热稳定性聚合酶的实例包括Taq DNA聚合酶的Stoffel片段。该聚合酶由于基因操作而缺少5′至3′核酸外切酶活性,并且由于Taq聚合酶天然缺少3′至5′核酸外切酶活性,因此不存在3′至5′活性。Tth DNA聚合酶衍生自嗜热栖热菌(Thermus thermophilus),并且可得自EpicentreTechnologies,Molecular Biology Resource Inc.或Perkin-Elmer Corp.。其它有用的缺少3′核酸外切酶活性的DNA聚合酶包括Vent[R](exo-),其得自New England Biolabs,Inc.(从带有来自古细菌Thermococcus litoralis的DNA聚合酶基因的大肠杆菌(E.coli)菌株纯化得到)和来源于黄栖热菌(Thermus flavus)并且得自Amersham Corporation的HotTub DNA聚合酶。其它优选的热稳定并且缺少5′至3′核酸外切酶活性和3′至5′核酸外切酶活性的酶包括AmpliTaq Gold。可以使用至少基本相同的其它DNA聚合酶,如其它N-末端截短的水生栖热菌(Thermus aquaticus)(Taq)DNA聚合酶I。称为KlenTaq I和KlenTaq LA的聚合酶非常适合该目的。当然,根据本发明还可以使用具有这些特性的任何其它聚合酶。
在探针3′端添加一个或多个核苷酸的条件将取决于所使用的特定的酶,并且通常将按照所使用的酶的生产商所推荐的条件进行。
特异性部分
一般在复杂性减小步骤之后,在本文所述的方法中包含特异性步骤。“特异性部分”是指区分目标核酸的步骤,优选地,在等位基因电平。也就是说,特异性部分是等位基因特异性步骤(例如,基因分型或SNP分析)。尽管可以通过将等位基因特异性探针与模板(即以上复杂性减小步骤的产物)简单杂交来达到一定特异性电平,但是在优选的实施方式中,特异性步骤包括酶促步骤。也就是说,酶促步骤的保真度改善了等位基因区分的特异性。优选的酶包括如本文中更详细地描述的DNA聚合酶、RNA聚合酶和连接酶。
如上所述的聚合酶还可以适合于特异性步骤。
多种连接酶是已知的并且适合于在本文所述的方法中使用。在Lehman,Science,186:790-797(1974);Engler等人,DNA Ligases,3-30页,在Boyer主编,The Enzymes,Vol.15B(Academic Press,New York,1982)中;等中描述了示例性连接酶。优选的连接酶包括T4DNA连接酶、T7DNA连接酶、大肠杆菌(E.coil)DNA连接酶、Taq连接酶、Pfu连接酶和Tth连接酶。它们的使用规程是熟知的,例如,Sambrook等人(以上引用);Barany,PCR Methodsan Applications,1:5-16(1991);Marsh等人,Strategies,5:73-76(1992);等。通常,连接酶需要存在用于连接至邻接链的3′羟基的5′磷酸基。优选的连接酶包括热稳定的或(嗜热的)连接酶,如pfu连接酶、Tth连接酶、Taq连接酶和Ampligase TM DNA连接酶(EpicentreTechnologies,Madison,Wis.)。Ampligase具有低钝端连接活性。
在某些实施方式中,连接酶是具有最小错配连接的连接酶。可以通过用特异性更强的NAD+-依赖性连接酶,如大肠杆菌(E.coli)连接酶和(热稳定)Taq连接酶替换特异性较低的T4DNA连接酶来提高连接酶的特异性。连接反应中NAD类似物的使用进一步提高了连接反应的特异性。参见,授权于Wallace等人的美国专利No.5,508,179。
在一个实施方式中,使用固定化目标进行特异性部分。也就是说,将复杂性减小步骤的产物固定在固体载体上,如本文所列。如本文所讨论的,特异性反应的目标被称为“特异性目标”。也就是说,复杂性减小步骤的产物是特异性目标。
在一个实施方式中,所述载体是与初始复杂性减小步骤中相同的载体。在该实施方式中,在特异性测定之前从固体载体移除目标核酸。可以通过使杂交复合物变性的任何方法除去目标核酸,从而导致目标核酸的释放。如本领域技术人员所理解的,在该实施方式中,目标核酸不共价结合至固体载体。也就是说,目标探针稳定连接到载体。也就是说,尽管探针的连接不必需是共价的,但是它足够稳定以经受得住杂交复合物的变性和未连接的目标核酸的移除。
在替代实施方式中,特异性目标处于溶液中。也就是说,在复杂性减小步骤之后,固定的目标核酸和目标探针之间的杂交复合物变性并且从杂交复合物洗脱修饰的目标探针。在某一实施方式中,在溶液中分析特异性目标。在替代实施方式中,将液相特异性目标固定在连续固体载体上。
这些特异性测定(即基因分型技术)分成五个一般种类:(1)依赖于常规杂交方法的技术,其使用严格性条件的改变(温度、缓冲条件等)来辨别检测位置的核苷酸;(2)添加碱基(“the base”),从而与检测位置的核苷酸形成碱基对的延伸技术;(3)连接技术,其依赖于连接酶的特异性(或者,在一些情况下,依赖于化学技术的特异性),从而如果在检测位置存在完全互补(perfect complementarity),则优先发生连接反应;(4)切割技术,其也依赖于酶促或化学特异性,从而如果存在完全互补,则优先发生切割;和(5)组合这些方法的技术。一般地,参见专利No.6,890,741、6,913,884、7,955,794、7582,420和8,288,103以及美国专利公开2013-0244882,以上专利作为参考并入本文。
在某些实施方式中,进行延伸基因分型。在该实施方式中,许多技术可以用于将核苷酸添加至与检测位置临近的目标序列杂交的探针的读取位置。通过依赖于酶促特异性,优先添加完全互补碱基。本文所述的一些方法依赖于检测位置处核苷酸的酶促掺入。可以使用多种本领域熟知的方法进行,如单碱基延伸或多碱基延伸。在某些实施方式中,基因分型伴随着不使用链终止核苷酸的引物延伸。照此,认为这种基因分型是多碱基延伸。所述方法包括提供设计用于检测给定SNP的一个等位基因的询问寡核苷酸(interrogatoroligonucleotide)。通过探查到的不同的SNP等位基因的数目确定寡核苷酸的数目。例如,如果探查到1000个SNP,则每个具有两个等位基因,因此2000个寡核苷酸将是必需的。询问子(interrogator)与含有SNP的DNA延伸互补,其中每个询问子的末端碱基对应于SNP位置,或者其中SNP-特异性位置位于询问子的最后1、2、3或4个核苷酸内。在一些实施方式中,询问子不是引物的末端位置,而是位于距引物3′端1、2、3、4、5或6个核苷酸的位置处。例如,当SNP具有A和C等位基因时,提供了以T和G结尾的询问子,并且在一些实施方式中,可以将其固定在单独的元件(珠)上以检测所述两个等位基因。尽管匹配和错配均将杂交至给定等位基因,但是仅匹配可以用作DNA聚合酶延伸反应的引物。因此,在探针与目标DNA杂交后,进行聚合酶反应。这导致在存在dNTP的情况下,与DNA聚合酶的杂交物的延伸。
在某些实施方式中,由于未延伸的探针或引物可以与延伸引物在捕获探针结合中竞争,因此期望从测定混合物,并且具体地,从固体载体中除去未延伸或未反应的探针或引物。相对于延伸引物,未延伸的引物的浓度可以相对较高,这是因为通常需要大量过量的引物以产生有效的引物退火。因此,可以使用一些不同的技术来辅助未延伸的探针或引物的除去。这些通常包括基于通过结合至固体载体,保护反应的引物并降解未延伸的引物,和分离未反应和反应的引物除去未反应的引物的方法。
扩增部分
在该实施方式中,本文提供了可以在表征多核苷酸的方法中使用的包括多核苷酸扩增的方法以及核酸扩增反应的产物(即扩增子)。适合的扩增方法包括目标扩增和信号扩增。目标扩增包括要检测的目标序列的扩增(即复制),从而导致目标分子数目的显著增加。目标扩增策略包括(但不限于)聚合酶链反应(PCR)、链置换扩增(SDA)、核酸序列基扩增(NASBA)和滚环扩增(RCA)。这些扩增策略对本领域技术人员来说是熟知的,并且可以容易地选择用于所述方法。
作为另外一种选择,替代技术使用目标作为模板来复制信号探针,而不是扩增目标,这使得少量目标分子能够导致产生大量信号检验,然后可以进行检测。信号扩增策略包括连接酶链反应(LCR)、循环探针技术(CPT)、侵入性裂解技术,如Invader技术、Q-β复制酶(QR)技术和“扩增探针”,如“分支DNA”的使用,其导致产生了结合单个目标序列的多重标记探针。
所有这些方法都可以包含引物核酸(包括核酸类似物),其杂交至目标序列以形成杂交复合物,并且以修饰引物形成修饰的引物的方式添加酶。例如,PCR通常需要两种引物,dNTP和DNA聚合酶;LCR需要两种引物,其邻近杂交目标序列和连接酶;CPT需要一种可切割引物和切割酶;侵入性裂解需要两种引物和裂解酶;等。因此,一般地,将目标核酸加入至包含必要扩增部分的反应混合物,并形成修饰的引物。
一般地,将修饰的引物用作二次反应的目标序列,然后其产生多个扩增链,并且可以如本文所列进行检测。根据需要,以多种方式除去未反应的引物,如本领域技术人员将理解的和如本文所列的。因此,反应从向目标序列中加入引物核酸开始,其形成了杂交复合物。一旦引物和目标序列之间形成杂交复合物,则使用酶(有时称为“扩增酶”)修饰引物。如对于本文所列的所有方法,可以在测定期间的任意点(在加入引物之前、期间或之后)加入酶。酶的身份将取决于所使用的扩增技术。类似地,修饰将取决于扩增技术。
在某些实施方式中,目标扩增技术是聚合酶链反应(PCR)。PCR是广泛使用和描述的,并且包括引物延伸与热循环结合使用以扩增目标序列;参见美国专利No.4,683,195和4,683,202和PCR Essential Data,J.W.Wiley&sons,Ed.C.R.Newton,1995,以上全部文献作为参考并入本文。另外,存在多个PCR的变化,其也在本发明中得到应用,包括“定量竞争性PCR”或“QC-PCR”、“任意引物PCR”或“AP-PCR”、“免疫-PCR”、“Alu-PCR”、“PCR单链构象多态性”或“PCR-SSCP”、“反转录酶PCR”或“RT-PCR”、“生物素俘获PCR”、“小载体PCR(vectorette PCR)”、“锅饼式(panhandle)PCR”和“PCR选择cDNA减少(subtraction)”、“等位基因-特异性PCR”等。应理解本领域技术人员可以容易地选择可以在本文所述的方法中使用的PCR的适当改变。
在某些实施方式中,扩增反应是如本文所述的多路扩增反应。在一个实施方式中,扩增反应使用多个PCR引物来扩增多个靶标序列。在该实施方式中,使用多个扩增引物对同时扩增多个目标序列。
在替代实施方式中,多路PCR反应使用如本文所述的通用引物。也就是说,通用PCR引物杂交至目标序列上的通用引发位点并借此扩增多个靶标序列。该实施方式潜在地是优选的,这是因为它仅需要有限的PCR引物。也就是说,少至一个引物对可以扩增多个靶标序列。
如前所述,使用人DNA作为模板产生了Golden Gate扩增子(Cold Spring HarbSymp Quant Biol.2003;68:69-78.Highly parallel SNP genotyping.Fan JB等人)。基于等位基因,所得扩增子具有表示为P1和P2的两个引物之一。此外,在所有扩增子上存在通用反向引物(“反向P3”)。
P1:TCTCGTCGCTCATCAACT(SEQ ID NO:82)
P2:GAGTCGAGGTCATATCGT(SEQ ID NO:83)
反向P3:GTCTGCCTATAGTGAGTC(SEQ ID NO:84)
将使用16个循环的第二轮PCR用于添加称为“P1_barcode_A”和“P1_barcode_B”的等位基因条型码引物。使用含有多个脱氧尿嘧啶残基的延伸的通用反向引物(“UniversaldU Reverse”)。
P1_barcode_A:
/5phos/TTTTTTTTTTTTTTTCCTTCCXXTTTTCTTCTTCTTCAAGAAGAAGATCTCGTCGCTCATCAACT(SEQ ID NO:85)
P2_barcode_B:
/5phos/TTTTTTTTTTTTTTTCCTTCCXXTTTTTTAATTAATTTTGTTGTTGTGAGTCGAGGTCATATCGT(SEQ ID NO:86)
Universal dU Reverse:
ATACGGCG/dU/CCACCGACC/dU/CAGCGTC/dU/GCCTATAG/dU/G AGTC(SEQ ID NO:87)
其中/5phos/表示5′磷酸盐,/du/是脱氧尿嘧啶残基,并且X是无碱基部分。
PCR之后,将样品与USER酶(New England Biolabs,Ipswich,MA)在37℃培育2.5小时以产生单链缺口,在此存在dU残基。将样品加热至65℃,10min以除去碎片化的DNA并产生3′突出。使用PCR Cleanup试剂盒(Qiagen)纯化样品。
通过加热至65℃并缓慢冷却,将样品与含有胆固醇的低聚“P3_Chol”以1∶1的摩尔比退火。
P3_Chol:
ACCGACACTGCGTCTGCCTATAGTGAGTC/iSp9//3CholTEG/(SEQ ID NO:88)
其中/iSp9/表示9-原子三甘醇间隔臂,并且/3CholTEG/表示3′胆固醇TEG(三甘醇)部分。
脂质双分子层由1,2-二植烷酰基-sn-甘油基-3-胆碱磷酸(Avanti极性脂)形成。双分子层在Teflon中横跨电平直径~20微米的孔。以~2.5ng/ml的浓度,将M2-NNN-MspA加入至双分子层的底侧(grounded side)。一旦***单个孔,则隔室内冲入实验缓冲液以避免进一步***。Axopatch-200B膜片钳放大器(Axon Instruments)将180mV的电压施加到双分子层两侧并测量离子电流。使用4-极贝塞耳滤波器以50kHz低通滤波模拟信号,然后以低通滤波频率的5倍数字化。使用以LabWindows/CVI写入的常规软件(National Instruments)控制数据采集。
双分子层两侧~60μl的隔室含有0.4M KCl、1mM EDTA、1mM DTT、1mM ATP、10mMMgCl2和10mM HEPES/KOH缓冲液,pH 8.0的实验缓冲液。以150nM,将Hel308野生型Tga用作马达。
图23示出了根据一些实施方式对于适合用作各自条型码的第一说明性多核苷酸序列(SEQ ID NO:89)和第二说明性多核苷酸序列(SEQ ID NO:90),可以作为时间函数产生的示例性模拟信号。在图24中可以看出(1)对应于第一说明性多核苷酸序列通过Hel308解旋酶通过孔的分步易位的模拟信号具有随时间的特征图形,其包括相对高信号电平的两个“峰”,然后下降,然后是相对低信号电平的另外两个“峰”,而(2)对应于第二说明性多核苷酸序列通过Hel308解旋酶通过孔的分步易位的模拟信号包括相对低信号电平的两个“峰”,然后升高,然后是相对高信号电平的另外两个“峰”。因此,可以预期可以(例如)使用图形匹配来容易地彼此区分包括不同特征的真实信号,如模拟信号(1)和(2),并因此可以有利于区分彼此的测定结果。例如,图24A-24D示出了根据一些实施方式对于适合用作各自条型码的第一和第二说明性多核苷酸序列可以作为时间函数产生的示例性模拟信号。可以看出在图24A和24B中通常表示为“条形码”的模拟信号的部分包括相对高信号电平的两个“峰”,然后下降,然后是相对低信号电平的另外两个“峰”,并因此可以理解为对应于第一说明性多核苷酸序列。还可以看出在图24C和24D中通常表示为“条形码”的模拟信号的部分包括相对低信号电平的两个“峰”,然后升高,然后是相对高信号电平的另外两个“峰”,并因此可以理解为对应于第二说明性多核苷酸序列。
在另一个实例中,以与本实施例中如上所述的类似方式,将2NNN MspA孔***DPhPC脂质双分子层。缓冲液包括400mM KCl、10mM HEPES pH 8、5mM MgCl2和1mM EDTA。所述试剂包括1mM DTT和1mM ATP。酶包括约150mM的Hel308 Tga。DNA为约10nM,并且将测序的单链(表示RS1801131 SNP1和SNP2)杂交至含有胆固醇的多核苷酸。使用后处理解释这些链测序期间所获得的信号,其包括电平查找(level-finding)和使用算法与预期序列比对,如在本文其它处所述的。
图25A和25B分别示出了根据一些实施方式对于适合用作各自条型码的第一和第二说明性多核苷酸序列可以作为时间函数产生的示例性模拟信号。图25A和25B中的虚线框内的模拟信号部分对于分别表示为rs1801131SNP1和rs1801131 SNP2的序列分别包括了不同的图像,并因此可以用作各自的条型码。用作条型码的序列与图23中所示的相同。
图26A-26D分别示出了根据一些实施方式对于适合用作各自条型码的第一和第二说明性多核苷酸序列作为时间函数产生的示例性测量信号。图26A和26B中的虚线框内的测量信号部分分别包括不同的图形,可以看出其对应于表示为rs1801131 SNP1的序列条型码,而图26C和26D中的虚线框内的测量信号部分分别包括不同的图形,可以看出其对应于表示为rs1801131 SNP2的序列条型码,并且该条形码与表示为rs1801131 SNP1的序列的条型码是容易区分的。
作为参考并入
在整个发明申请中,使用或未使用括号引用了多个出版物。在本发明申请中,这些出版物的公开内容出于所有目的以它们的全部内容作为参考并入本文,其包括(但不限于)为了更充分地描述本发明公开所述领域目前的工艺水平。
其它替代实施方式
应注意可以使用执行指令(例如,软件命令)以实施本文所公开的操作的多种类型的数据处理器环境(例如,在一个或多个数据处理器上)来实施本文所提供的***和方法。非限制性实例包括在单个通用计算机或工作站上,或者在网络***上,或者在客户-服务器配置中,或者在应用业务供应商配制中实施。例如,可以在多种不同类型的处理装置上,通过包含设备处理子***可执行的程序命令的程序代码来实施本文所述的方法和***。软件程序命令可以包括源编码、结果代码、机器代码或者可操作以导致处理***执行本文所述的方法和操作的任何其它存储数据。然而,还可以使用其它装置,如配置用于执行本文所述的方法和***的固件或甚至适当设计的硬件。
还应注意所述***和方法可以包括通过用于一个或多个数据处理装置通讯的网络(例如,局域网、广域网、因特网、它们的组合等)、光学纤维介质、载波、无线网络等传输的数据信号。数据信号可以具有任何或所有本文所公开的提供给装置或装置所提供的数据。
可以在一种或多种不同类型的计算机实施数据存储中存储和实施***和方法数据(例如,关联、数据输入、数据输出、中间数据结果、最终数据结果等),如不同类型的存储设备和编程构造体(例如,RAM、ROM、闪存、平面文件、数据库、编程数据结构、编程变量、IF-THEN(或类似类型)声明构造体等)。注意数据结构描述了在供计算机程序使用的数据库、程序、存储器或其它计算机可读介质中组织和存储数据中使用的格式。
还可以在多种不同类型的计算机可读存储介质上提供***和方法,其包括含有在通过处理器的执行中用于实施方法操作和执行本文所述的***的命令(例如,软件)的计算机存储机制(例如,永久性介质,如CD-ROM、磁盘、RAM、闪存、计算机硬盘等)。
此外,本文所提供的计算机组件、软件模块、功能、数据存储和数据结构可以直接或间接彼此连接以使得能够进行它们的操作所需的数据流动。还注意模块或处理器包括(但不限于)执行软件操作的代码单元,并且可以(例如)作为代码的子例程单元或作为代码的软件功能单元或作为对象(如在面向对象的范例中那样)或作为小程序或以计算机脚本语言或作为另一类型的计算机代码实施。基于实际情形,可以将软件部件和/或功能定位在单个计算机上或分布在多个计算机上。
尽管已参考所公开的实施方式描述了本发明公开,但是本领域技术人员将容易地理解以上详述的具体实例和研究仅是本发明公开的说明。应理解在不背离本发明公开的精神的情况下,可以做出多种改变。因此,本发明公开仅限于以下权利要求。
Claims (71)
1.表征目标多核苷酸的方法,所述方法包括:
(a)施加跨越与Hel308解旋酶和目标多核苷酸接触的孔的电势差,其中所述Hel308解旋酶与包含ATP的底物接触,并且其中所述Hel308解旋酶具有选自SEQ ID NO:4、SEQ IDNO:5、SEQ ID NO:22和SEQ ID NO:23的Hel308基序,并且所述ATP具有从10μM至1mM的浓度;
(b)测量在所述Hel308解旋酶的完整易位循环期间通过所述孔的所述目标多核苷酸的每个核苷酸产生的至少两个电信号;和
(c)根据在步骤(b)测量的所述至少两个信号表征所述目标多核苷酸。
2.根据权利要求1所述的方法,其中,表征所述目标多核苷酸包含鉴别以下中的一种或多种:所述目标多核苷酸的序列、所述目标多核苷酸的修饰、所述目标多核苷酸的长度、所述目标多核苷酸的身份、所述目标多核苷酸的来源和所述目标多核苷酸的二级结构。
3.根据权利要求2所述的方法,还包括重复步骤(a)-(c)一次或多次。
4.根据权利要求1所述的方法,其中,所述目标多核苷酸的易位处于通过势差施加于易位通过所述孔的多核苷酸上的外加作用力的相反方向。
5.根据权利要求1所述的方法,其中,所述目标多核苷酸的易位处于通过势差施加于易位通过所述孔的多核苷酸上的外加作用力的方向。
6.根据权利要求1所述的方法,其中,所述电信号为选自电流、电压、隧穿、电阻、电位、电导率和横向电测量的测量值。
7.根据权利要求6所述的方法,其中,所述电信号包括穿过所述孔的电流。
8.根据权利要求1所述的方法,其中,使用得自完整易位循环的两个分步易位步骤的电信号表征所述目标多核苷酸中的一个或多个核苷酸残基,与使用得自完整易位循环的单个电信号的一个或多个核苷酸的表征相比,其准确度大于50%。
9.根据权利要求1所述的方法,其中,所述孔为生物孔。
10.根据权利要求9所述的方法,其中,所述生物孔为多肽孔。
11.根据权利要求9所述的方法,其中,所述生物孔为多核苷酸孔。
12.根据权利要求10所述的方法,其中,所述多肽孔具有5个核苷酸或以下的缢缩区。
13.根据权利要求10所述的方法,其中,所述多肽孔包含耻垢分枝杆菌孔蛋白A(MspA)。
14.根据权利要求13所述的方法,其中,所述MspA由SEQ ID NO:1的氨基酸序列组成。
15.根据权利要求1所述的方法,其中,所述孔为固态孔。
16.根据权利要求1所述的方法,其中,所述孔为生物和固态杂交孔。
17.根据权利要求16所述的方法,其中,所述生物和固态杂交孔为多肽-固态杂交孔。
18.根据权利要求16所述的方法,其中,所述生物和固态杂交孔为多核苷酸-固态杂交孔。
19.根据权利要求1所述的方法,其中,所述目标多核苷酸选自由单链多核苷酸和双链多核苷酸组成的组。
20.根据权利要求1所述的方法,其中,所述目标多核苷酸是部分双链多核苷酸。
21.调节目标多核苷酸通过孔的分步易位步骤的方法,所述方法包括:
(a)施加跨越与Hel308解旋酶和所述目标多核苷酸接触的孔的电势差,其中所述Hel308解旋酶具有选自SEQ ID NO:4、SEQ ID NO:5、SEQ ID NO:22和SEQ ID NO:23的Hel308基序;
(b)将所述Hel308解旋酶与不同于底物参比浓度的Hel308解旋酶底物浓度接触,与所述参比浓度相比,所述底物浓度在分步易位步骤持续时间中所产生的变化与所述底物浓度的差异成正比,其中所述底物包含ATP并且所述ATP具有从10μM至1mM的浓度,和
(c)测量在所述Hel308解旋酶的完整易位循环期间通过所述孔的所述目标多核苷酸的每个核苷酸产生的至少两个电信号。
22.根据权利要求21所述的方法,还包括根据所述至少两个电信号表征所述目标多核苷酸。
23.根据权利要求22所述的方法,其中,表征所述目标多核苷酸包括鉴别以下中的一种或多种:所述目标多核苷酸的序列、所述目标多核苷酸的修饰、所述目标多核苷酸的长度、所述目标多核苷酸的身份、所述目标多核苷酸的来源和所述目标多核苷酸的二级结构。
24.根据权利要求21所述的方法,其中,所述底物浓度为所述Hel308解旋酶底物的半饱和浓度。
25.根据权利要求21所述的方法,其中,所述参比浓度为所述Hel308解旋酶底物的饱和浓度。
26.根据权利要求21所述的方法,其中,所述底物浓度和所述参比浓度两者均不是所述底物的饱和浓度。
27.根据权利要求21所述的方法,其中,所述底物浓度和所述参比浓度是所述Hel308解旋酶底物的半饱和浓度。
28.根据权利要求21所述的方法,其中,所述Hel308解旋酶底物为三磷腺苷(ATP)。
29.根据权利要求21所述的方法,其中,所述目标多核苷酸的易位处于通过势差施加于易位通过所述孔的多核苷酸上的外加作用力的相反方向。
30.根据权利要求21所述的方法,其中,所述目标多核苷酸的易位处于通过势差施加于易位通过所述孔的多核苷酸上的外加作用力的方向。
31.根据权利要求21所述的方法,其中,所述电信号为选自电流、电压、隧穿、电阻、电位、电导率和横向电测量的测量。
32.根据权利要求31所述的方法,其中,所述电信号包括穿过所述孔的电流。
33.根据权利要求21所述的方法,其中,使用得自完整易位循环的两个分步步骤的电信号表征所述目标多核苷酸中的一个或多个核苷酸残基,与使用得自完整易位循环的单个电信号的一个或多个核苷酸的表征相比,其准确度大于50%。
34.根据权利要求21所述的方法,其中,与所述参比浓度相比,在较低底物浓度时,表征所述目标多核苷酸中的一个或多个核苷酸残基具有更大的准确度。
35.根据权利要求21所述的方法,其中,所述孔为生物孔。
36.根据权利要求35所述的方法,其中,所述生物孔为多肽孔。
37.根据权利要求35所述的方法,其中,所述生物孔为多核苷酸孔。
38.根据权利要求36所述的方法,其中,所述多肽孔具有5个核苷酸或以下的缢缩区。
39.根据权利要求36所述的方法,其中,所述多肽孔包含耻垢分枝杆菌孔蛋白A(MspA)。
40.根据权利要求39所述的方法,其中,所述MspA由SEQ ID NO:1所示的氨基酸序列组成。
41.根据权利要求21所述的方法,其中,所述孔为固态孔。
42.根据权利要求21所述的方法,其中,所述孔为生物和固态杂交孔。
43.根据权利要求42所述的方法,其中,所述生物和固态杂交孔为多肽-固态杂交孔。
44.根据权利要求42所述的方法,其中,所述生物和固态杂交孔为多核苷酸-固态杂交孔。
45.根据权利要求21所述的方法,其中,所述目标多核苷酸选自由单链多核苷酸和双链多核苷酸组成的组。
46.根据权利要求21所述的方法,其中,所述目标多核苷酸是部分双链多核苷酸。
47.根据权利要求1所述的方法,其中,所述表征包括应用改进型维特比算法。
48.根据权利要求1所述的方法,还包括:
(d)在步骤(c)后,改变至少一种参数从而改变所述目标多核苷酸通过所述孔的由所述Hel308解旋酶造成的一个或多个分步易位步骤的时机;和
(e)使用所述改变的至少一种参数,重复步骤(a)-(c)。
49.根据权利要求48所述的方法,还包括合并步骤(c)和(e)期间产生的信号并基于合并的电信号表征所述目标多核苷酸。
50.根据权利要求49所述的方法,其中,所述改变的至少一种参数选自温度、盐浓度、辅因子浓度、ATP产物浓度、pH和所使用的具体的Hel308解旋酶。
51.根据权利要求1所述的方法,其中,所述表征包括检测和鉴别所述至少两个电信号中的电平以及基于所述检测和鉴别的电平确定和输出所述目标多核苷酸的序列。
52.根据权利要求51所述的方法,其中,所述检测和鉴别所述至少两个电信号中的电平包括输出全电平。
53.根据权利要求51所述的方法,其中,所述检测和鉴别所述至少两个电信号中的电平包括输出分电平。
54.根据权利要求51所述的方法,其中,所述检测和鉴别所述至少两个电信号中的电平包括输出所有电平。
55.根据权利要求51所述的方法,其中,所述检测和鉴别所述至少两个电信号中的电平包括输出电平标识符。
56.根据权利要求52所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将全电平作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的相关置信信息选择并输出至少一个所调用的序列。
57.根据权利要求53所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将分电平作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的相关置信信息选择并输出至少一个所调用的序列。
58.根据权利要求54所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将所有电平作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的相关置信信息选择并输出至少一个所调用的序列。
59.根据权利要求55所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将电平标识符作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的相关置信信息选择并输出至少一个所调用的序列。
60.根据权利要求52所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将全电平作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的部分的相关置信信息选择多个所调用序列的部分并将其彼此连接。
61.根据权利要求53所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将分电平作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的部分的相关置信信息选择多个所调用序列的部分并将其彼此连接。
62.根据权利要求54所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将所有电平作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的部分的相关置信信息选择多个所调用序列的部分并将其彼此连接。
63.根据权利要求55所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将电平标识符作为输入,基于所述输入调用多个序列,和基于所调用的多个序列的部分的相关置信信息选择多个所调用序列的部分并将其彼此连接。
64.根据权利要求52所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将全电平作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于所调用的序列与模型序列比较的相关置信信息选择和输出至少一个调用序列。
65.根据权利要求53所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将分电平作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于所调用的序列与模型序列比较的相关置信信息选择和输出至少一个调用序列。
66.根据权利要求54所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将所有电平作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于所调用的序列与模型序列比较的相关置信信息选择和输出至少一个调用序列。
67.根据权利要求55所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将电平标识符作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于所调用的序列与模型序列比较的相关置信信息选择和输出至少一个调用序列。
68.根据权利要求52所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将全电平作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于多个调用序列的部分与模型序列比较的相关置信信息选择多个调用序列的部分并将其彼此连接。
69.根据权利要求53所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将分电平作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于多个调用序列的部分与模型序列比较的相关置信信息选择多个调用序列的部分并将其彼此连接。
70.根据权利要求54所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将所有电平作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于多个调用序列的部分与模型序列比较的相关置信信息选择多个调用序列的部分并将其彼此连接。
71.根据权利要求55所述的方法,其中,基于所述检测和鉴别的电平,所述确定和输出所述目标多核苷酸的序列包括将电平标识符作为输入,基于所述输入调用多个序列,将所调用的序列与模型序列比较,和基于多个调用序列的部分与模型序列比较的相关置信信息选择多个调用序列的部分并将其彼此连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110859489.8A CN113528632A (zh) | 2013-11-26 | 2014-11-26 | 用于多核苷酸测序的组合物和方法 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361909316P | 2013-11-26 | 2013-11-26 | |
US61/909,316 | 2013-11-26 | ||
PCT/US2014/067582 WO2015081178A1 (en) | 2013-11-26 | 2014-11-26 | Compositions and methods for polynucleotide sequencing |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110859489.8A Division CN113528632A (zh) | 2013-11-26 | 2014-11-26 | 用于多核苷酸测序的组合物和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105934522A CN105934522A (zh) | 2016-09-07 |
CN105934522B true CN105934522B (zh) | 2021-07-20 |
Family
ID=52424089
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201480064725.5A Active CN105934522B (zh) | 2013-11-26 | 2014-11-26 | 用于多核苷酸测序的组合物和方法 |
CN202110859489.8A Pending CN113528632A (zh) | 2013-11-26 | 2014-11-26 | 用于多核苷酸测序的组合物和方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110859489.8A Pending CN113528632A (zh) | 2013-11-26 | 2014-11-26 | 用于多核苷酸测序的组合物和方法 |
Country Status (12)
Country | Link |
---|---|
US (5) | US9689033B2 (zh) |
EP (3) | EP4282983A3 (zh) |
JP (3) | JP6800749B2 (zh) |
CN (2) | CN105934522B (zh) |
AU (3) | AU2014354726B2 (zh) |
CA (2) | CA2926871C (zh) |
DK (2) | DK3556869T3 (zh) |
ES (2) | ES2958715T3 (zh) |
FI (1) | FI3556869T3 (zh) |
LT (1) | LT3074534T (zh) |
SI (1) | SI3074534T1 (zh) |
WO (2) | WO2015081211A2 (zh) |
Families Citing this family (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6457811B2 (ja) * | 2011-09-23 | 2019-01-23 | オックスフォード ナノポール テクノロジーズ リミテッド | ポリマー単位を含むポリマーの解析 |
WO2013121224A1 (en) | 2012-02-16 | 2013-08-22 | Oxford Nanopore Technologies Limited | Analysis of measurements of a polymer |
CN112646019B (zh) | 2012-04-10 | 2022-08-16 | 牛津纳米孔科技公开有限公司 | 突变胞溶素孔 |
GB201222928D0 (en) | 2012-12-19 | 2013-01-30 | Oxford Nanopore Tech Ltd | Analysis of a polynucleotide |
GB201313477D0 (en) | 2013-07-29 | 2013-09-11 | Univ Leuven Kath | Nanopore biosensors for detection of proteins and nucleic acids |
CN105934522B (zh) * | 2013-11-26 | 2021-07-20 | 伊鲁米那股份有限公司 | 用于多核苷酸测序的组合物和方法 |
EP3097210B2 (en) | 2014-01-22 | 2022-04-20 | Oxford Nanopore Technologies Limited | Method for attaching one or more polynucleotide binding proteins to a target polynucleotide |
SG10202111841TA (en) | 2014-02-19 | 2021-12-30 | Univ Washington | Nanopore-based analysis of protein characteristics |
CN106574300A (zh) | 2014-05-02 | 2017-04-19 | 牛津纳米孔技术公司 | 改善目标多核苷酸相对于跨膜孔移动的方法 |
JP6709213B2 (ja) * | 2014-10-16 | 2020-06-10 | オックスフォード ナノポール テクノロジーズ リミテッド | ポリマーの解析 |
WO2017027518A1 (en) * | 2015-08-10 | 2017-02-16 | Stratos Genomics, Inc. | Single molecule nucleic acid sequencing with molecular sensor complexes |
CN108521782A (zh) * | 2015-10-21 | 2018-09-11 | 豪夫迈·罗氏有限公司 | 含氟聚合物作为疏水层以支持用于纳米孔的脂质双层形成的用途 |
WO2017149316A1 (en) | 2016-03-02 | 2017-09-08 | Oxford Nanopore Technologies Limited | Mutant pore |
CN109071618A (zh) | 2016-04-06 | 2018-12-21 | 牛津纳米孔技术公司 | 突变体孔 |
WO2017180141A1 (en) | 2016-04-14 | 2017-10-19 | Lockheed Martin Corporation | Selective interfacial mitigation of graphene defects |
EP3475447B1 (en) * | 2016-06-23 | 2021-04-21 | F. Hoffmann-La Roche AG | Formation and calibration of nanopore sequencing cells |
EP3497233B1 (en) | 2016-08-08 | 2021-11-10 | F. Hoffmann-La Roche AG | Basecalling for stochastic sequencing processes |
CA3038347A1 (en) | 2016-10-19 | 2018-04-26 | Illumina, Inc. | Methods for chemical ligation of nucleic acids |
GB201707122D0 (en) | 2017-05-04 | 2017-06-21 | Oxford Nanopore Tech Ltd | Pore |
SG11201913174PA (en) | 2017-06-30 | 2020-01-30 | Vib Vzw | Novel protein pores |
CA3232118A1 (en) * | 2017-09-15 | 2019-03-21 | Illumina, Inc. | Sequence-detection system |
EP3682023A4 (en) * | 2017-09-15 | 2021-06-02 | Illumina, Inc. | TUNING AND CALIBRATION CHARACTERISTICS OF A SEQUENCE DETECTION SYSTEM |
US10310959B2 (en) | 2017-11-07 | 2019-06-04 | Bank Of America Corporation | Pre-deployment validation system using intelligent databases |
WO2019243421A1 (en) * | 2018-06-21 | 2019-12-26 | F. Hoffmann-La Roche Ag | Tunneling junctions for sequencing |
US20190390268A1 (en) * | 2018-06-26 | 2019-12-26 | Electronic Biosciences, Inc. | Controlled nanopore translocation utilizing extremophilic replication proteins |
JP2020031557A (ja) * | 2018-08-28 | 2020-03-05 | 株式会社日立ハイテクノロジーズ | 生体分子分析装置 |
WO2020084705A1 (ja) * | 2018-10-24 | 2020-04-30 | 株式会社日立ハイテク | 生体ポリマ分析デバイス及びそれを用いた分析装置、並びに分析方法 |
WO2021056598A1 (zh) | 2019-09-29 | 2021-04-01 | 北京齐碳科技有限公司 | 一种Mmup单体变体及其应用 |
CN114807317A (zh) * | 2021-01-22 | 2022-07-29 | 上海羿鸣生物科技有限公司 | 一种优化的dna线性扩增方法及试剂盒 |
GB202107192D0 (en) * | 2021-05-19 | 2021-06-30 | Oxford Nanopore Tech Ltd | Method |
CN113322180B (zh) * | 2021-05-31 | 2022-07-15 | 中国科学院物理研究所 | 基于纳米孔测序的力谱分析方法和分析装置 |
US20230298693A1 (en) * | 2022-02-28 | 2023-09-21 | Rajant Health Incorporated | Alignment-free variant calling |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013057495A3 (en) * | 2011-10-21 | 2013-08-22 | Oxford Nanopore Technologies Limited | Method of characterizing a target polynucleotide using a pore and a hel308 helicase |
Family Cites Families (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4469863A (en) | 1980-11-12 | 1984-09-04 | Ts O Paul O P | Nonionic nucleic acid alkyl and aryl phosphonates and processes for manufacture and use thereof |
US5235033A (en) | 1985-03-15 | 1993-08-10 | Anti-Gene Development Group | Alpha-morpholino ribonucleoside derivatives and polymers thereof |
US5034506A (en) | 1985-03-15 | 1991-07-23 | Anti-Gene Development Group | Uncharged morpholino-based polymers having achiral intersubunit linkages |
US4683202A (en) | 1985-03-28 | 1987-07-28 | Cetus Corporation | Process for amplifying nucleic acid sequences |
US4683195A (en) | 1986-01-30 | 1987-07-28 | Cetus Corporation | Process for amplifying, detecting, and/or-cloning nucleic acid sequences |
US5216141A (en) | 1988-06-06 | 1993-06-01 | Benner Steven A | Oligonucleotide analogs containing sulfur linkages |
US5602240A (en) | 1990-07-27 | 1997-02-11 | Ciba Geigy Ag. | Backbone modified oligonucleotide analogs |
US5386023A (en) | 1990-07-27 | 1995-01-31 | Isis Pharmaceuticals | Backbone modified oligonucleotide analogs and preparation thereof through reductive coupling |
US5644048A (en) | 1992-01-10 | 1997-07-01 | Isis Pharmaceuticals, Inc. | Process for preparing phosphorothioate oligonucleotides |
US5637684A (en) | 1994-02-23 | 1997-06-10 | Isis Pharmaceuticals, Inc. | Phosphoramidate and phosphorothioamidate oligomeric compounds |
US5508179A (en) | 1994-03-18 | 1996-04-16 | Bio-Rad Laboratories, Inc. | Use of deoxyribose nicotinamide adenine dinucleotide to enhance the specificity of NAD+ -dependent ligation reactions |
US5795782A (en) * | 1995-03-17 | 1998-08-18 | President & Fellows Of Harvard College | Characterization of individual polymer molecules based on monomer-interface interactions |
EP2383776B1 (en) | 1999-06-22 | 2015-02-25 | President and Fellows of Harvard College | Solid state nanopore device for evaluating biopolymers |
WO2001057269A2 (en) | 2000-02-07 | 2001-08-09 | Illumina, Inc. | Nucleic acid detection methods using universal priming |
US7955794B2 (en) | 2000-09-21 | 2011-06-07 | Illumina, Inc. | Multiplex nucleic acid reactions |
US6913884B2 (en) | 2001-08-16 | 2005-07-05 | Illumina, Inc. | Compositions and methods for repetitive use of genomic DNA |
US7582420B2 (en) | 2001-07-12 | 2009-09-01 | Illumina, Inc. | Multiplex nucleic acid reactions |
GB0505971D0 (en) | 2005-03-23 | 2005-04-27 | Isis Innovation | Delivery of molecules to a lipid bilayer |
EP2126588A1 (en) | 2007-02-20 | 2009-12-02 | Oxford Nanopore Technologies Limited | Formation of lipid bilayers |
EP2682460B1 (en) * | 2008-07-07 | 2017-04-26 | Oxford Nanopore Technologies Limited | Enzyme-pore constructs |
EP2344891B1 (en) | 2008-09-22 | 2016-03-16 | University Of Washington | Msp nanopores and related methods |
GB0905140D0 (en) * | 2009-03-25 | 2009-05-06 | Isis Innovation | Method |
EP2614156B1 (en) * | 2010-09-07 | 2018-08-01 | The Regents of The University of California | Control of dna movement in a nanopore at one nucleotide precision by a processive enzyme |
US10175195B2 (en) | 2011-07-27 | 2019-01-08 | The Board Of Trustees Of The University Of Illinois | Nanopore sensors for biomolecular characterization |
EP2798083B1 (en) | 2011-12-29 | 2017-08-09 | Oxford Nanopore Technologies Limited | Method for characterising a polynucelotide by using a xpd helicase |
CA2861808C (en) | 2011-12-29 | 2021-02-23 | Oxford Nanopore Technologies Limited | Enzyme method |
JP5809572B2 (ja) | 2012-01-30 | 2015-11-11 | ルネサスエレクトロニクス株式会社 | 半導体装置 |
CN112646019B (zh) | 2012-04-10 | 2022-08-16 | 牛津纳米孔科技公开有限公司 | 突变胞溶素孔 |
WO2013185137A1 (en) * | 2012-06-08 | 2013-12-12 | Pacific Biosciences Of California, Inc. | Modified base detection with nanopore sequencing |
AU2013291763B2 (en) * | 2012-07-19 | 2019-06-20 | Oxford Nanopore Technologies Limited | Modified helicases |
CN105934522B (zh) * | 2013-11-26 | 2021-07-20 | 伊鲁米那股份有限公司 | 用于多核苷酸测序的组合物和方法 |
-
2014
- 2014-11-26 CN CN201480064725.5A patent/CN105934522B/zh active Active
- 2014-11-26 DK DK19161148.2T patent/DK3556869T3/da active
- 2014-11-26 WO PCT/US2014/067639 patent/WO2015081211A2/en active Application Filing
- 2014-11-26 SI SI201431264T patent/SI3074534T1/sl unknown
- 2014-11-26 EP EP23188555.9A patent/EP4282983A3/en active Pending
- 2014-11-26 WO PCT/US2014/067582 patent/WO2015081178A1/en active Application Filing
- 2014-11-26 LT LTEP14830903.2T patent/LT3074534T/lt unknown
- 2014-11-26 EP EP19161148.2A patent/EP3556869B1/en active Active
- 2014-11-26 CA CA2926871A patent/CA2926871C/en active Active
- 2014-11-26 JP JP2016521718A patent/JP6800749B2/ja active Active
- 2014-11-26 CN CN202110859489.8A patent/CN113528632A/zh active Pending
- 2014-11-26 ES ES19161148T patent/ES2958715T3/es active Active
- 2014-11-26 ES ES14830903T patent/ES2735015T3/es active Active
- 2014-11-26 DK DK14830903.2T patent/DK3074534T3/da active
- 2014-11-26 CA CA3157586A patent/CA3157586A1/en active Pending
- 2014-11-26 AU AU2014354726A patent/AU2014354726B2/en active Active
- 2014-11-26 FI FIEP19161148.2T patent/FI3556869T3/fi active
- 2014-11-26 EP EP14830903.2A patent/EP3074534B1/en active Active
- 2014-11-26 US US14/554,741 patent/US9689033B2/en active Active
-
2017
- 2017-05-26 US US15/606,354 patent/US10364462B2/en active Active
-
2019
- 2019-06-12 US US16/439,022 patent/US11041196B2/en active Active
-
2020
- 2020-04-14 AU AU2020202489A patent/AU2020202489B2/en active Active
- 2020-11-25 JP JP2020195453A patent/JP7349972B2/ja active Active
-
2021
- 2021-05-25 US US17/329,482 patent/US11879155B2/en active Active
-
2022
- 2022-09-16 AU AU2022231786A patent/AU2022231786A1/en active Pending
-
2023
- 2023-05-30 JP JP2023088480A patent/JP2023116553A/ja active Pending
- 2023-12-08 US US18/533,863 patent/US20240167085A1/en active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013057495A3 (en) * | 2011-10-21 | 2013-08-22 | Oxford Nanopore Technologies Limited | Method of characterizing a target polynucleotide using a pore and a hel308 helicase |
Non-Patent Citations (3)
Title |
---|
DNA Base-Calling from a Nanopore Using a Viterbi Algorithm;Winston Timp,et al.;《Biophysical Journa》;20120516;第102卷;摘要,L37页左栏第1段至L39页右栏第1段 * |
Nanopore sensors for nucleic acid analysis;Bala Murali Venkatesan,et al.;《NATURE NANOTECHNOLOGY》;20110918;第6卷;第615-624页 * |
Stochastic sensors inspired by biology;Hagan Bayley,et al.;《NATURE》;20010913;第413卷;第226-230页 * |
Also Published As
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105934522B (zh) | 用于多核苷酸测序的组合物和方法 | |
EP1225234B1 (en) | A method for sequencing nucleic acids | |
CN109312401B (zh) | 连接的核酸的纳米孔测序方法 | |
US9267168B2 (en) | Methods and compositions for isolating template nucleic acids | |
EP4063521B1 (en) | Method of nanopore sequencing | |
US10697006B2 (en) | Hairpin-mediated amplification method | |
US20170037465A1 (en) | Methods for Nucleic Acid Base Determination | |
CN113039285A (zh) | 用于纳米孔测序的液体样品工作流程 | |
EP4362028A1 (en) | Mutant aerolysin and uses thereof | |
Chu | Progress Toward Single-Molecule Nanopore DNA Sequencing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |