CN108603228B - 通过分析无细胞dna确定肿瘤基因拷贝数的方法 - Google Patents
通过分析无细胞dna确定肿瘤基因拷贝数的方法 Download PDFInfo
- Publication number
- CN108603228B CN108603228B CN201680081723.6A CN201680081723A CN108603228B CN 108603228 B CN108603228 B CN 108603228B CN 201680081723 A CN201680081723 A CN 201680081723A CN 108603228 B CN108603228 B CN 108603228B
- Authority
- CN
- China
- Prior art keywords
- dataset
- genetic loci
- genetic
- baseline
- sequencing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
- C12Q1/6886—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material for cancer
-
- 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/6809—Methods for determination or identification of nucleic acids involving differential detection
-
- 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
- C12Q1/6874—Methods for sequencing involving nucleic acid arrays, e.g. sequencing by hybridisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/10—Ploidy or copy number detection
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/10—Sequence alignment; Homology search
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B99/00—Subject matter not provided for in other groups of this subclass
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Organic Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Genetics & Genomics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Medical Informatics (AREA)
- Evolutionary Biology (AREA)
- Immunology (AREA)
- Molecular Biology (AREA)
- Biochemistry (AREA)
- General Engineering & Computer Science (AREA)
- Microbiology (AREA)
- Pathology (AREA)
- Oncology (AREA)
- Hospice & Palliative Care (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
- Algebra (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
Abstract
本文提供了方法以改进核酸样品中拷贝数变异的自动检测。这些方法提供了用于确定样品内遗传基因座的基线拷贝数的改进方法,减少了由于遗传基因座的特征、样品制备和探针耗尽导致的变异。
Description
交叉引用
本申请要求于2015年12月17日提交的美国临时申请号62/269,051的优先权,其在此通过引用以其整体并入。
背景
癌症是由个体正常细胞内突变的积累引起的,其中至少一些导致细胞分裂调节不当。这种突变通常包括拷贝数变异,其中肿瘤基因组内基因的拷贝数相对于受试者的非癌性细胞增加或减少。
利用检测和表征肿瘤细胞中的拷贝数变异来监测肿瘤的进展,预测患者的结果,和改善治疗选择。然而,常规方法是对细胞样品进行的,细胞样品通常通过痛苦和时间密集的活检获得。这种活检通常也只能检查受试者体内的一小部分肿瘤细胞,因此并不总是代表肿瘤细胞的群体。对于不需要细胞活检、荧光原位杂交(FISH)、比较基因组杂交阵列或定量荧光聚合酶链式反应(PCR)测定的用于肿瘤中拷贝数变异的更简单、更快速的检测存在需求。
使用测序数据确定拷贝数变异的一个特殊挑战是,出于与真实拷贝数无关的原因,遗传基因座的覆盖深度会出现变异。例如,即使对于样品中存在的在同样拷贝数的单独遗传基因座,扩增效率、PCR效率和鸟嘌呤-胞嘧啶含量也可导致不同的覆盖深度。为了改善拷贝数检测,需要消除由于这种效应的偏倚的改进方法。
概述
对在来源于无细胞体液的样品的肿瘤细胞中检测拷贝数变异的改进方法存在相当大的需求。本方面解决了这种需求并且提供了另外的优点。在一方面,本公开内容提供了一种方法,包括:(a)获得受试者的无细胞体液样品的脱氧核糖核酸(DNA)分子的测序读段;(b)从序列读段生成第一数据集,对于多个遗传基因座中的每个遗传基因座,所述第一数据集包含与测序读段覆盖(“读段覆盖”)有关的定量量度;(c)通过进行饱和平衡校正和探针效率校正来校正第一数据集;(d)确定第一数据集的基线读段覆盖,其中基线读段覆盖涉及饱和平衡和探针效率;和(e)确定多个遗传基因座中每个遗传基因座相对于基线读段覆盖的拷贝数状态。在一些实施方案中,对于多个遗传基因座中的每个遗传基因座,第一数据集包含(i)遗传基因座的鸟嘌呤-胞嘧啶含量(“GC含量”)相关的定量量度。在一些实施方案中,该方法包括在(c)之前从第一数据集去除为高变异遗传基因座的遗传基因座,其中去除包括:(i)拟合涉及鸟嘌呤-胞嘧啶含量相关的定量量度和遗传基因座的测序读段覆盖的定量量度的模型;和(ii)从遗传基因座去除至少10%的遗传基因座,其中去除遗传基因座包括去除与该模型最不同的遗传基因座,由此提供基线定线(baselining)遗传基因座的第一数据集。在一些实施方案中,该方法包括去除至少45%的遗传基因座。
在一些实施方案中,进行饱和平衡校正包括通过以下将基线定线数据遗传基因座的第一数据集转换为饱和校正数据集:(i)对于来自基线定线遗传基因座的第一数据集每个遗传基因座,确定来源于遗传基因座的DNA分子链在测序读段内被代表的概率有关的定量量度;(ii)通过将基线定线遗传基因座的第一数据集的读段覆盖与基线定线遗传基因座的第一数据集的GC含量和与来源于基线定线遗传基因座的第一数据集中的每个基因座的DNA链在测序读段内被代表的概率有关的定量量度相关联来确定读段覆盖的第一转换;和(iii)将所述第一转换应用于来自基线定线遗传基因座的第一数据集的每个遗传基因座的读段覆盖以提供饱和校正数据集,其中所述饱和校正数据集包括基线定线遗传基因座的第一数据集的转换的读段覆盖的第一集。
在一些实施方案中,确定第一转换包括(i)确定与基线定线遗传基因座的第一数据集的读段覆盖的中心趋势有关的量度;(ii)基于遗传基因座的GC含量和与来源于遗传基因座的DNA链被代表在测序读段内的概率相关的定量量度,确定拟合与基线定线遗传基因座的第一数据集的读段覆盖的中心趋势有关的量度的函数;和(iii)对于基线定线遗传基因座的第一数据集的每个遗传基因座,确定由所述函数预测的读段覆盖与所述读段覆盖的差异,其中所述差异为转换的读段覆盖。在一些实施方案中,该函数是表面近似。在本文提供的一些实施方案中,表面近似是二维二次多项式。
在一些实施方案中,进行探针效率校正包括通过以下转换饱和校正数据集为探针效率校正数据集:(i)从饱和校正数据集去除相对于转换的读段覆盖的第一集为高变异遗传基因座的遗传基因座,从而提供基线定线遗传基因座的第二数据集;(ii)确定与基线定线遗传基因座的第二数据集的探针效率相关的转换的读段覆盖的第一集的第二转换;和(iii)利用第二转换将基线定线遗传基因座的第二数据集的转换的读段覆盖的第一数据集进行转换,从而提供探针效率校正数据集,其中探针效率校正数据集包括基线定线遗传基因座的第二数据集的转换的读段覆盖的第二集。在一些实施方案中,从第一数据集去除为高变异遗传基因座的遗传基因座包括:(i)拟合涉及GC含量和饱和校正数据集的转换的读段覆盖的第一集的模型;和(ii)从饱和校正数据集去除至少10%的遗传基因座,其中去除遗传基因座包括去除与该模型最不同的遗传基因座,由此提供基线定线遗传基因座的第二数据集。在本文提供的一些实施方案中,去除是至少45%的遗传基因座。
在一些实施方案中,探针效率通过对一种或更多种参考样品进行饱和平衡校正来确定,其中探针效率是通过进行饱和平衡校正获得的转换的读段覆盖。在一些实施方案中,一种或更多种参考样品是来自无癌症的受试者的无细胞体液样品。在本文提供的一些实施方案中,一种或更多种参考样品是来自具有癌症的受试者的无细胞体液样品,其中对应的遗传基因座未经历拷贝数改变。
在一些实施方案中,确定第二转换包括(i)将对来自一种或更多种参考样品的遗传基因座确定的探针效率拟合到来自基线定线遗传基因座的第二数据集的读段覆盖的第一集;(ii)将基线定线遗传基因座的第二数据集的每个遗传基因座的转换的读段覆盖除以基于(i)的拟合的预测的探针效率。在一些实施方案中,该方法还包括:(f)通过将基线定线遗传基因座的第二数据集的转换的读段覆盖与基线定线遗传基因座的第二数据集的GC含量和与来源于基线定线遗传基因座的第二数据集中的每个基因座的DNA链在测序读段内被代表的概率有关的定量量度相关联,确定转换的读段覆盖的第二集的第三转换;和(g)将第三转换应用于转换的读段覆盖的第二集以提供第四数据集,其中第四数据集包括转换的定量读段覆盖的第三集。
在一些实施方案中,无细胞体液样品的DNA使用与来自遗传基因座的集的遗传基因座的至少一部分互补的一种或更多种寡核苷酸探针富集遗传基因座的集。在一些实施方案中,来自遗传基因座的集的每个遗传基因座的GC含量是与来自遗传基因座的集的遗传基因座的至少一部分互补的一种或更多种寡核苷酸探针的鸟嘌呤-胞嘧啶含量的中心趋势相关的量度。在一些实施方案中,遗传基因座的读段覆盖是与对应于一种或更多种寡核苷酸探针的遗传基因座区域的读段覆盖的中心趋势相关的量度。在一些实施方案中,进行饱和平衡校正和进行探针效率校正包括拟合Langmuir模型,其中Langmuir模型包含探针效率(K)和饱和平衡常数(Isat)。在一些实施方案中,K和Isat对于一种或更多种寡核苷酸探针中的每种寡核苷酸探针根据经验确定。在一些实施方案中,进行饱和平衡校正和进行探针校正包括将遗传基因座的读段覆盖拟合于Langmuir模型,假设遗传基因座以相同的拷贝数状态存在,从而提供基线读段覆盖。在一些实施方案中,相同的拷贝数状态是二倍体。在一些实施方案中,基线读段覆盖是依赖于探针效率和饱和平衡的函数。
在一些实施方案中,确定拷贝数状态包括比较遗传基因座的读段覆盖与基线读段覆盖。在一些实施方案中,无细胞体液选自由以下组成的组:血清、血浆、尿和脑脊液。在一些实施方案中,读段覆盖通过将测序读段映射到参考基因组来确定。在一些实施方案中,获得测序读段包括将衔接子与来自受试者的无细胞体液的DNA分子连接。在一些实施方案中,DNA分子是双链DNA分子并且衔接子与双链DNA分子连接,使得每个衔接子将DNA分子的互补链不同地加标签以提供加标签链。在一些实施方案中,确定来源于遗传基因座的DNA链在所述测序读段内被代表的概率相关的定量量度包括将测序读段分选为配对读段和非配对读段,其中(i)每个配对读段对应于所述集中由来源于双链多核苷酸分子的第一加标签链和第二差异加标签互补链生成的序列读段,并且(ii)每个非配对读段代表第一加标签链,并且来源于双链多核苷酸分子的第二差异加标签互补链在所述序列读段集中的所述序列读段之中没有被代表。在一些实施方案中,该方法还包括确定映射到一个或更多个遗传基因座中的每一个的(i)所述配对读段和(ii)所述非配对读段的定量量度,以基于映射到每个基因座的配对读段和非配对读段相关的所述定量量度确定所述样品中映射到所述一个或更多个遗传基因座中的每一个的全部双链DNA分子相关的定量量度。在一些实施方案中,衔接子包括条形码序列。
在一些实施方案中,确定读段覆盖包括基于测序读段映射到参考基因组的位置和条形码序列叠并测序读段。在一些实施方案中,遗传基因座包括一种或更多种致癌基因。在一些实施方案中,一种方法包括通过确定受试者的种系基因组对其为杂合的基线定线遗传基因座内的变体的相对量,确定受试者的肿瘤细胞中基线定线遗传基因座的至少子集已经经历拷贝数改变。在一些实施方案中,变体的相对量不是近似相等的。在一些实施方案中,将变体的相对量对其不是近似相等的基线定线遗传基因座从基线定线遗传基因座去除,从而提供等位基因频率校正的基线定线遗传基因座。在一些实施方案中,等位基因频率校正的基线定线遗传基因座在前述权利要求的任一项的方法中用作基线定线基因座。
在另一方面,本公开内容提供了包括以下的方法:接收受试者的无细胞体液样品的脱氧核糖核酸(DNA)分子的测序读段到存储器中;用计算机处理器执行代码以进行以下步骤:从序列读段生成对于多个遗传基因座中的每个遗传基因座包含与测序读段覆盖(“读段覆盖”)有关的定量量度的第一数据集;通过进行饱和平衡校正和探针效率校正来校正第一数据集;确定第一数据集的基线读段覆盖,其中基线读段覆盖涉及饱和平衡和探针效率;和确定多个遗传基因座中每个遗传基因座相对于基线读段覆盖的拷贝数状态。
在另一方面,本公开内容提供了一种系统,包括:网络;数据库,所述数据库连接到网络,包括配置成存储核酸(例如,DNA)序列数据的计算机存储器;生物信息学计算机,其包括计算机存储器和一个或更多个计算机处理器,所述计算机连接到所述网络;其中计算机还包括机器可执行代码,当由一个或更多个计算机处理器执行时,该机器可执行代码复制存储在数据库上的核酸(例如,DNA)序列数据,将复制的数据写入生物信息学计算机中的存储器,并进行包括以下的步骤:从核酸(例如,DNA)序列数据生成第一数据集,其包含对于多个遗传基因座中的每个遗传基因座,与测序读段覆盖(“读段覆盖”)相关的定量量度;通过进行饱和平衡校正和探针效率校正来校正第一数据集;确定对于第一数据集的基线读段覆盖,其中基线读段覆盖与饱和平衡和探针效率有关;并确定多个遗传基因座中每个遗传基因座相对于基线读段覆盖的拷贝数状态。在一些实施方案中,数据库连接于DNA测序仪。
通过引用并入
本说明书中提及的所有出版物、专利和专利申请通过引用并入本文,其程度如同每一个单独的出版物、专利或专利申请被具体和单独地指明通过引用并入的相同程度。
附图简述
本发明的新颖的特征在所附权利要求书中具体陈述。本发明的特征和优势的更好的理解将通过参考以下详细说明和附图来获得,详细说明陈述了其中利用本发明的原理的说明性实施方案,在附图中:
图1说明示例性致癌基因和用于序列捕获探针的靶。
图2说明基因水平信号相对于跨三个加标(spike-in)的理论拷贝数和跨加标基因的探针水平信号变异。
图3说明将诱饵量与独特分子计数相关联的诱饵优化实验。
图4A和图4B说明p(图4A)和GC含量(图4B)对独特分子计数的非线性影响。
图5说明没有进行饱和或探针效率校正的每个探针的独特分子计数。
图6说明饱和校正后每个探针的独特分子计数。
图7说明饱和校正后和探针效率校正后每个探针的独特分子计数。
图8说明真实拷贝数与探针饱和和探针效率相关的独特分子计数之间的相互作用的提议的Langmuir模型。
图9说明典型临床样品中饱和校正、探针效率校正和第二轮探针效率校正后用于基线定线遗传基因座的探针信号降噪。
图10A和图10B说明将饱和后校正的UMC针对参考样品中确定的探针效率作图,以便进行探针效率校正。图10A是来自肿瘤细胞中无拷贝数改变的受试者。图10B是来自肿瘤细胞中具有拷贝数改变的受试者。
图11说明患者样品中饱和和探针效率校正的拷贝数变异检测的最终报告。样品上方的星号表示基于校正的信号和次等位基因频率校正的基线优化检测的基因扩增。
图12示出了被编程或以其他方式配置为实现本公开内容的方法的计算机系统1201。
图13说明使用本公开内容的方法测量的基因ERBB2的观察到的拷贝数(CN)vs.理论CN。实心圆点表示观察到的~2拷贝数(二倍体样品),空心圆点表示检测到的扩增事件,且粗的水平虚线表示平均基因CN截断值。
图14说明与对照方法(方形)相比使用本公开内容的方法(圆点)测量的基因ERBB2的观察到的拷贝数(CN)vs.理论CN。实心圆点表示观察到的~2拷贝数(二倍体样品),空心圆点表示检测到的扩增事件,且粗的水平虚线表示平均基因CN截断值。
图15说明将探针拷贝数针对用于本公开内容的方法(三角形)vs.对照方法(X)的验证研究中使用的探针作图。
详细说明
定义
如本文所用的术语“遗传变异”通常是指受试者的核酸样品或基因组中的改变、变异或多态性。这样的改变、变异或多态性可以是相对于参考基因组的,该参考基因组可以是所述受试者或其他个体的参考基因组。单核苷酸多态性(SNP)是多态性的一种形式。在一些实例中,一个或更多个多态性包括一个或更多个单核苷酸变异(SNV)、插入、缺失、重复、小的插入、小缺失、小重复、结构变体接合、可变长度串联重复和/或侧翼序列。拷贝数变异(CNV)、颠换和其他重排也是遗传变异的形式。基因组改变可以是碱基变化、插入、缺失、重复、拷贝数变异或颠换。
如本文使用的,术语“多核苷酸”通常指包含一个或更多个核酸亚单位的分子。多核苷酸可包含选自腺苷(A)、胞嘧啶(C)、鸟嘌呤(G)、胸腺嘧啶(T)和尿嘧啶(U)或其变体的一个或更多个亚单位。核苷酸可包含A、C、G、T或U或其变体。核苷酸可包括可被掺入到增长的核酸链中的任何亚单位。此类亚单位可以是A、C、G、T或U,或为一个或更多个互补A、C、G、T或U特异性的或与嘌呤(即,A或G或其变体)或嘧啶(即,C、T或U或其变体)互补的任何其他亚单位。亚单位能够使得个体核酸碱基或碱基组(例如,AA、TA、AT、GC、CG、CT、TC、GT、TG、AC、CA或尿嘧啶-其对应物)得以解析。在一些实例中,多核苷酸为脱氧核糖核酸(DNA)或核糖核酸(RNA)或其衍生物。多核苷酸可以是单链或双链的。
如本文所用的术语“受试者”通常是指动物如哺乳动物物种(例如,人),或禽类(例如,鸟)物种,或其他生物体如植物。更具体地,所述受试者可以是脊椎动物、哺乳动物、小鼠、灵长类动物、猿或人。动物包括,但不限于,农场动物、运动动物和宠物。受试者可以是健康的个体,患有或疑似患有疾病或有患病倾向的个体,或需要治疗或疑似需要治疗的个体。受试者可以是患者。
术语“基因组”通常指生物体的遗传信息的整体。基因组可以以DNA或以RNA来编码。基因组可包含编码蛋白的编码区域和非编码区域。基因组可包含生物体中所有染色体的序列的集合。例如,人类基因组具有总计46条染色体。所有这些的序列共同构成人类基因组。
术语“衔接子(adaptor)”、“衔接子(adaptor)”和“标签”在本说明书全文中作为同义词使用。可通过任何方法,包括连接、杂交或其他方法,使衔接子或标签与待“加标签”的多核苷酸序列偶联。
如本文所用的术语“文库衔接子(library adaptor)”或“文库衔接子(libraryadaptor)”通常是指其身份(例如,序列)可用来区分生物样品(在本文中也称为“样品”)中的多核苷酸的分子(例如,多核苷酸)。
如本文所用的术语“测序衔接子”通常是指适合于允许测序仪器对靶多核苷酸进行测序(例如通过与靶多核苷酸相互作用以实现测序)的分子(例如,多核苷酸)。测序衔接子允许测序仪器对靶多核苷酸进行测序。在一个实例中,测序衔接子包含与附接至测序系统的固体支持物如流通池的捕获多核苷酸杂交或结合的核苷酸序列。在另一个实例中,测序衔接子包含与多核苷酸杂交或结合以生成发夹环的核苷酸序列,该发夹环允许测序系统对靶多核苷酸进行测序。测序衔接子可包含测序仪基序,该测序仪基序可以是与其他分子(例如,多核苷酸)的流通池序列互补并可被测序系统使用以对靶多核苷酸进行测序的核苷酸序列。该测序仪基序还可包括用于在测序如合成测序中使用的引物序列。该测序仪基序可包括使文库衔接子与测序系统偶联和对靶多核苷酸进行测序所需的序列。
除非另外确定,当如本文所用的术语“至少”、“至多”或“约”在系列之前时是指所述系列的每一个成员。
与参考数值相关的术语“约”及其语法上的等同表述可包括从该值加上或减去10%的值的范围。例如,“约10”的量可包括9至11的量。在其他实施方案中,与参考数值相关的术语“约”可包括从该值加上或减去10%、9%、8%、7%、6%、5%、4%、3%、2%或1%的值的范围。
与参考数值相关的术语“至少”及其语法上的等同表述可包括该参考数值和大于该值的值。例如,“至少10”的量可包括值10和大于10的任何数值,如11、100和1,000。
与参考数值相关的术语“至多”及其语法上的等同表述可包括该参考数值和小于该值的值。例如,“至多10”的量可包括值10和10以下的任何数值,如9、8、5、1、0.5和0.1。
术语“定量量度”是指数量的任何量度,包括绝对和相对量度。定量量度可以是,例如,数目(例如计数)、百分比、程度或阈值。
术语“读段覆盖”是指被原始序列读段或由处理序列读段的覆盖,例如从原始序列读段推断的独特分子计数。
术语“基线读段覆盖”是指基于给定探针参数例如GC含量、探针效率、连接效率或下拉效率,包含二倍体基因组环境的样品中探针的预期读段覆盖。
本文所用的“探针”是指包含官能性的多核苷酸。官能性可以是可检测标记(荧光)、结合部分(生物素)、或固体支持物(磁力可吸引的颗粒或芯片)。
“互补性”是指核酸与另一核酸序列通过传统Watson-Crick类型或其他非传统类型形成氢键的能力。互补性百分比指示核酸分子中可与第二核酸序列形成氢键(Watson-Crick碱基配对)的残基的百分比(10个中的5、6、7、8、9、10个分别为50%、60%、70%、80%、90%和100%互补)。“完全互补”是指核酸序列中所有的连续残基将与第二核酸序列中相同数目的连续残基氢键键合。
本文所用的“大体上互补”是指跨越8、9、10、11、12、13、14、15、16、17、18、19、20、21、22、23、24、25、30、35、40、45、50或更多个核苷酸的区域至少60%、65%、70%、75%、80%、85%、90%、95%、97%、98%、99%或100%的互补性程度,或是指在严格条件下杂交的两个核酸。序列同一性,诸如为了评估互补性百分比的目的,可通过任何合适的比对算法测量,所述比对算法包括但不限于Needleman-Wunsch算法(参见例如在万维网网页:ebi.ac.uk/Tools/psa/emboss_needle/nucleotide.html可得的EMBOSS Needle aligner,任选地以默认设置)、BLAST算法(参见例如在blast.ncbi.nlm.nih.gov/Blast.cgi可得的BLAST比对工具,任选地以默认设置)、或Smith-Waterman算法(参见例如在万维网网页:ebi.ac.uk/Tools/psa/emboss water/nucleotide.html可得的EMBOSS Water aligner,任选地以默认设置)。最佳比对可使用选择的算法的任何合适的参数(包括默认参数)评估。
“杂交”是指其中一个或更多个多核苷酸反应形成经由核苷酸残基的碱基之间的氢键被稳定化的复合体的反应。氢键键合可通过Watson Crick碱基配对、Hoogstein结合、或以符合碱基互补性的任何其他序列特异性方式发生。复合体可包括形成双链体结构的两条链、形成多链复合体的三条或更多条链、单个自杂交链、或这些的任何组合。杂交反应可构成更广泛的过程诸如PCR的起始、或内切核酸酶对多核苷酸的酶促裂解中的一个步骤。与第一序列互补的第二序列被称为第一序列的“互补物”。应用于多核苷酸的术语“可杂交的”是指杂交反应中多核苷酸形成经由核苷酸残基的碱基之间的氢键键合被稳定化的复合体的能力。
术语“严格杂交条件”指在其下多核苷酸将优先与其靶子序列杂交,并且在较小程度上与其它序列杂交或根本不与其他序列杂交的条件。在核酸杂交实验的上下文中的“严格杂交”是序列依赖性的,并在不同环境参数下是不同的。核酸杂交的广泛指南见于Tijssen(1993)Laboratory Techniques in Biochemistry and Molecular Biology--Hybridization with Nucleic Acid Probes第I部分第2章“Overview of principles ofhybridization and the strategy of nucleic acid probe assays”,Elsevier,NewYork。
通常,高度严格的杂交和洗涤条件被选择为比在定义的离子强度和pH下特定序列的热解链温度(Tm)低约5℃。Tm为50%的靶序列与完全匹配的探针杂交的温度(在定义的离子强度和pH下)。非常严格的条件被选择为等于特定探针的Tm。
严格杂交条件包括包含水、缓冲剂(在pH 6-9或pH 7-8的磷酸盐、tris、SSPE或SSC缓冲液)、盐(钠或钾)、和变性剂(SDS、甲酰胺或吐温)的缓冲液和37℃-70℃、60℃-65℃的温度。
在DNA印迹或RNA印迹中的滤器上用于具有多于100个互补残基的互补核酸的杂交的严格杂交条件的实例为在42℃,50%福尔马林与1mg肝素,杂交被进行过夜。高度严格的洗涤条件的实例为在72℃,0.15M NaCl,持续约15分钟。严格洗涤条件的实例为在65℃,0.2X SSC洗涤,持续15分钟(参见,Sambrook等人对SSC缓冲液的描述)。通常,高严格洗涤之前进行低严格洗涤以去除背景探针信号。对于多于100个核苷酸的双链体的示例性中度严格洗涤为在45℃,1x SSC,持续15分钟。对于,例如,多于100个核苷酸的双链体的示例性低度严格洗涤为在40℃,4-6x SSC,持续15分钟。通常,比起在特定杂交测定中对于不相关探针观察到的2倍(或更高)的信噪比指示检测到特异性杂交。
在一方面,本公开内容提供了一种方法,包括:(a)获得来源于受试者的无细胞体液样品的脱氧核糖核酸(DNA)分子的测序读段;(b)从测序读段生成第一数据集,所述第一数据集包含对于多个遗传基因座中的每个遗传基因座,与以下相关的定量量度:(i)遗传基因座的鸟嘌呤-胞嘧啶含量和(ii)遗传基因座的测序读段覆盖相关的定量量度;(c)通过以下转换第一数据集为第二数据集:(i)从所述第一数据集去除相对于测序读段覆盖相关的定量量度为高变异遗传基因座的遗传基因座;(ii)对于来自剩余遗传基因座的第一集的每个遗传基因座,确定来源于遗传基因座的样品的DNA链在测序读段内被代表的概率有关的定量量度;(iii)通过将剩余遗传基因座的第一集的测序读段覆盖相关的定量量度与剩余遗传基因座的第一集的GC含量相关的定量量度和与来源于剩余遗传基因座的第一集中的每个基因座的DNA链在测序读段内被代表的概率有关的定量量度相关联来确定测序读段覆盖相关的定量量度的第一转换;和(iv)将所述第一转换应用于来自剩余遗传基因座的第一集的每个遗传基因座的序列读段覆盖以提供第二数据集,其中第二数据集包括剩余遗传基因座的第一集的测序读段覆盖的转换的定量量度的第一集。
在一些实施方案中,该方法还包括通过以下将第二数据集转换为第三数据集:(d)从第二数据集去除相对于测序读段覆盖的转换的定量量度的第一集为高变异遗传基因座的遗传基因座,从而提供剩余遗传基因座的第二集;(e)确定与剩余遗传基因座的第二集的效率相关的测序读段覆盖的转换的定量量度的第一集的第二转换和(f)利用第二转换将剩余遗传基因座的第二集的测序读段覆盖的转换的定量量度的第一集进行转换,从而提供第三数据集,其中第三数据集包括(d,i)的剩余遗传基因座的第二集的测序读段覆盖相关的转换的定量量度的第二集;
获得来自受试者的无细胞体液的DNA分子的测序读段
获得来自受试者的无细胞体液的DNA分子的测序读段可包括获得无细胞体液。示例性的无细胞体液是或可来源于血清、血浆、血液、唾液、尿、滑液、全血、淋巴液、腹水、间质液或细胞外液、细胞之间的间隙中的流体(包括龈沟液)、骨髓、脑脊液、唾液、粘液、痰液、精液、汗液、尿液或任何其他体液。无细胞体液可选自由以下组成的组:血浆、尿或脑脊液。无细胞体液可以是血浆。无细胞体液可以是尿。无细胞体液可以是脑脊液。
核酸分子,包括DNA分子,可从无细胞体液提取。DNA分子可以是基因组DNA。DNA分子可来自受试者的健康组织的细胞。DNA分子可来自已经历体细胞突变的非癌性细胞。DNA分子可来自母体样品中的胎儿。本领域技术人员将理解,在其中DNA分子来自母体样品中的胎儿的实施方案中,受试者可以是指胎儿,尽管样品是母体样品。DNA分子可来自受试者的癌变前细胞。DNA分子可来自受试者的癌性细胞。DNA分子可来自受试者的原代肿瘤内的细胞。DNA分子可来自受试者的继发肿瘤。DNA分子可以是循环DNA。循环DNA可包括循环肿瘤DNA(ctDNA)。DNA分子可以是双链或单链的。可选地,DNA分子可包括双链部分和单链部分的组合。DNA分子不必是无细胞的。在一些情形中,DNA分子可从样品分离。例如,DNA分子可以是从体液例如血清或血浆分离的无细胞DNA。
样品可包含不同量的核酸分子的基因组等同物。例如,约30ng DNA的样品可含有约10,000个单倍体人类基因组等同物,而在cfDNA的情况下,可含有约2000亿个单独的多核苷酸分子。类似地,约100ng DNA的样品可含有约30,000个单倍体人类基因组等同物,而在cfDNA的情况下,可含有约6000亿个单独的分子。
无细胞DNA分子可使用本领域已知的多种技术从体液分离和提取。在一些情况下,可使用商购可得的试剂盒诸如QiagenCirculating Nucleic Acid试剂盒方案分离、提取和制备无细胞核酸。在其他实例中,可使用Qiagen QubitTMdsDNA HS测定试剂盒方案、AgilentTMDNA 1000试剂盒或TruSeqTMSequencing Library Preparation;Low-Throughput(LT)方案来定量核酸。无细胞核酸可以是胎儿来源的(经由从妊娠受试者取得的流体),或可源自受试者自身的组织。无细胞核酸可来源于赘生物(例如肿瘤或腺瘤)。
通常,通过分配步骤(partitioning step)从体液中提取和分离无细胞核酸,在该分配步骤中,如在溶液中发现的无细胞核酸与细胞和体液的其他不可溶性组分分开。分配可包括但不限于技术诸如离心或过滤。在其他情况下,细胞并非首先与无细胞核酸分配,而是经裂解。在一个实例中,完整细胞的基因组DNA通过选择性沉淀来分配。无细胞核酸(包括DNA)可保持可溶性并可与不可溶性基因组DNA分开并被提取。通常,在不同试剂盒特定的添加缓冲液和其他洗涤步骤之后,可使用异丙醇沉淀来沉淀核酸。可使用进一步的清洁步骤诸如基于二氧化硅的柱以去除污染物或盐。可针对特定应用优化一般步骤。例如,可在整个反应中添加非特异性批量(bulk)载体核酸以优化该程序的某些方面诸如收率。
无细胞DNA分子可以是至多500个核苷酸的长度、至多400个核苷酸的长度、至多300个核苷酸的长度、至多250个核苷酸的长度、至多225个核苷酸的长度、至多200个核苷酸的长度、至多190个核苷酸的长度、至多180个核苷酸的长度、至多170个核苷酸的长度、至多160个核苷酸的长度、至多150个核苷酸的长度、至多140个核苷酸的长度、至多130个核苷酸的长度、至多120个核苷酸的长度、至多110个核苷酸的长度、或至多100个核苷酸的长度。
无细胞DNA分子可以是至少500个核苷酸的长度、至少400个核苷酸的长度、至少300个核苷酸的长度、至少250个核苷酸的长度、至少225个核苷酸的长度、至少200个核苷酸的长度、至少190个核苷酸的长度、至少180个核苷酸的长度、至少170个核苷酸的长度、至少160个核苷酸的长度、至少150个核苷酸的长度、至少140个核苷酸的长度、至少130个核苷酸的长度、至少120个核苷酸的长度、至少110个核苷酸的长度、或至少100个核苷酸的长度。尤其地,无细胞核酸可以是在140和180个核苷酸之间的长度。
无细胞DNA可包含来自健康组织和肿瘤的不同量的DNA分子。肿瘤来源的无细胞DNA可以是样品中无细胞DNA总量的至少0.1%、样品中无细胞DNA总量的至少0.2%、样品中无细胞DNA总量的至少0.5%、样品中无细胞DNA总量的至少0.7%、样品中无细胞DNA总量的至少1%、样品中无细胞DNA总量的至少2%、样品中无细胞DNA总量的至少3%、样品中无细胞DNA总量的至少4%、样品中无细胞DNA总量的至少5%、样品中无细胞DNA总量的至少10%、样品中无细胞DNA总量的至少15%、样品中无细胞DNA总量的至少20%、样品中无细胞DNA总量的至少25%、或样品中无细胞DNA总量的至少30%、或更多。
在一些情况中,DNA分子可在提取过程中被修剪并包括在100和400个核苷酸之间的长度的片段。在一些情况中,核酸可在提取后被修剪并可包括在100和400个核苷酸之间的长度的核苷酸。在某些情况下,DNA分子的长度已经在100到400个核苷酸之间,且另外的修剪并非有目的地实施。
受试者可以是动物。受试者可以是哺乳动物,诸如犬、马、猫、小鼠、大鼠或人类。受试者可以是人类。受试者可疑似患有癌症。受试者可此前接受过癌症诊断。受试者的癌症状态可以是未知的。受试者可以是雄性或雌性的。受试者可以是至少20岁、至少30岁、至少40岁、至少50岁、至少60岁或至少70岁。
测序可通过本领域已知的任何方法。例如,测序技术包括经典技术(例如,使用标记的终止子或引物的双脱氧测序反应(Sanger方法)以及平板或毛细管中的凝胶分离)和下一代技术。示例性技术包括使用可逆终止的标记的核苷酸的合成测序、焦磷酸测序、454测序、Illumina/Solexa测序、与标记的寡核苷酸探针文库的等位基因特异性杂交、使用等位基因特异性杂交至标记的克隆文库,随后连接的合成测序、在聚合步骤期间实时监测标记的核苷酸的掺入、聚合酶克隆测序(polony sequencing)、SOLiD测序靶向测序、单分子实时测序、外显子测序、基于电子显微术的测序、组测序(panel sequencing)、晶体管介导的测序、直接测序、随机鸟枪法测序、全基因组测序、杂交测序、毛细管电泳、凝胶电泳、双链体测序、循环测序、单碱基延伸测序、固相测序、高通量测序、大规模平行标志测序、乳液PCR、在较低变性温度共扩增-PCR(COLD-PCR)、多重PCR、可逆染料终止子测序、双末端测序(paired-end sequencing)、近期测序(near-term sequencing)、核酸外切酶测序、连接测序、短读段测序、单分子测序、实时测序、反向终止子测序、纳米孔测序、MS-PET测序及其组合。在一些实施方案中,测序方法是大规模平行测序,即,对至少100、1000、10,000、100,000、100万、1000万、1亿、或10亿个多核苷酸分子的任一种同时(或快速连续)测序。在一些实施方案中,测序可由基因分析仪,诸如例如可从Illumina或Applied Biosystems商业获得的基因分析仪进行。最近通过利用聚合酶或连接酶的连续或单一延伸反应以及通过与探针文库的单一或连续差示杂交已经证明了单独分子的测序。可由DNA测序仪(例如,设计用于进行测序反应的机器)进行测序。在一些实施方案中,DNA测序仪可包含数据库或连接到数据库,例如,包含DNA序列数据的数据库。
可使用的测序技术包括,例如,使用合成测序系统。在第一步骤中,将DNA剪切成约300-800个碱基对的片段,并将片段平端化。寡核苷酸衔接子然后连接到片段的末端。衔接子充当用于片段扩增和测序的引物。可使用例如含有5'-生物素标签的衔接子B将片段附接至DNA捕获珠,例如链霉抗生物素蛋白包被的珠。附着于珠的片段在油-水乳液的液滴内PCR扩增。结果是每个珠上克隆扩增的DNA片段的多个拷贝。在第二步骤中,珠被捕获在孔中(皮升大小)。焦磷酸测序对每个DNA片段平行进行。加入一个或更多个核苷酸生成光信号,其由测序仪器中的CCD相机记录。信号强度与掺入的核苷酸数量成正比。焦磷酸测序利用核苷酸添加时释放的焦磷酸(PPi)。在腺苷5'磷酸硫酸盐存在下,通过ATP硫酸化酶将PPi转化为ATP。萤光素酶使用ATP将萤光素转化为氧化萤光素,并且该反应生成被检测和分析的光。
可使用的DNA测序技术的另一个实例是来自Life Technologies Corporation(Carlsbad,Calif.)的Applied Biosystems的SOLiD技术。在SOLiD测序中,将基因组DNA剪切成片段,并将衔接子附接到片段的5'和3'末端以生成文库片段。可选地,可通过以下引入内部衔接子:连接衔接子到片段的5'和3'末端,环化片段,消化环化片段以生成内部衔接子,并将衔接子附接到所得片段的5'和3'末端以生成配偶-配对文库。接下来,在含有珠、引物、模板和PCR组分的微反应器中制备克隆珠群体。PCR后,模板被变性并富集珠以分离具有延伸的模板的珠。在选定的珠上的模板经历允许与载玻片键合的3'修饰。序列可通过连续杂交和连接具有由特定荧光团识别的中心确定的碱基(或碱基对)的部分随机的寡核苷酸来确定。记录颜色后,去除连接的寡核苷酸,且然后重复该过程。
可使用的DNA测序技术的另一个实例是离子半导体测序,例如使用LifeTechnologies(South San Francisco,Calif.)的Ion Torrent以商标ION TORRENT出售的系统。离子半导体测序描述于,例如Rothberg,等人,An integrated semiconductordevice enabling non-optical genome sequencing,Nature 475:348-352(2011);美国公布号2010/0304982;美国公布号2010/0301398;美国公布号2010/0300895;美国公布号2010/0300559;和美国公布号2009/0026082中,其每一个的内容通过引用以其整体并入。
可使用的测序技术的另一个实例是Illumina测序。Illumina测序是基于使用折叠PCR(fold-back PCR)和锚定引物在固体表面上扩增DNA。基因组DNA被片段化,且衔接子被添加到片段的5'和3'末端。附着于流通池通道表面的DNA片段被延伸并被桥式扩增。片段变成双链,且双链分子变性。多个循环的固相扩增、随后变性可在流通池的每个通道中产生相同模板的单链DNA分子的约1,000个拷贝的几百万个簇。引物、DNA聚合酶和四种荧光团标记的、可逆终止的核苷酸用于进行顺序测序。核苷酸掺入后,使用激光激发荧光团,拍摄图像并记录第一碱基的身份。来自每个掺入的碱基的3'终止子和荧光团被除去,并重复掺入、检测和鉴定步骤。按照这种技术的测序描述在美国专利第7,960,120号;美国专利第7,835,871号;美国专利第7,232,656号;美国专利第7,598,035号;美国专利第6,911,345号;美国专利第6,833,246号;美国专利第6,828,100号;美国专利第6,306,597号;美国专利第6,210,891号;美国公布号2011/0009278;美国公布号2007/0114362;美国公布号2006/0292611;和美国公布号2006/0024681,其每一个通过引用以其整体并入。
可使用的测序技术的另一个实例包括Pacific Biosciences(Menlo Park,Calif.)的单分子、实时(SMRT)技术。在SMRT中,四种DNA碱基中的每一种与四种不同荧光染料中的一种相连。这些染料是磷酸连接的(phospholinked)。单个DNA聚合酶与模板单链DNA的单个分子一起固定在零模波导(ZMW)的底部。将一个核苷酸掺入增长链需要几毫秒。在此期间,荧光标记物被激发并产生荧光信号,且荧光标签被切下。检测染料的相应荧光表明哪个碱基被掺入。重复该过程。
可使用的测序技术的另一个实例是纳米孔测序(Soni&Meller,2007,Progresstoward ultrafast DNA sequence using solid-state nanopores,Clin Chem 53(11):1996-2001)。纳米孔是直径为1纳米级别的小孔。将纳米孔浸入导电流体中并跨其施加电势导致由于离子通过纳米孔的传导而产生轻微的电流。流动的电流量对纳米孔的尺寸敏感。当DNA分子通过纳米孔时,DNA分子上的每个核苷酸以不同程度阻挡纳米孔。因此,当DNA分子通过纳米孔时,通过纳米孔的电流的变化代表DNA序列的读取。
可使用的测序技术的另一个实例涉及使用化学品敏感的场效应晶体管(chemFET)阵列来对DNA测序(例如,如美国公布号2009/0026082中描述的)。在该技术的一个实例中,可将DNA分子置于反应室中,并且模板分子可与结合聚合酶的测序引物杂交。在测序引物的3'末端将一个或更多个三磷酸掺入到新核酸链中可通过chemFET的电流变化来检测。阵列可有多个chemFET传感器。在另一个实例中,可将单个核酸与珠附接,并且可在珠上扩增核酸,并且可将单个珠转移到chemFET阵列上的单个反应室中,每个室具有chemFET传感器,并且核酸可被测序。
可使用的测序技术的另一个实例涉及使用电子显微镜,例如,由Moudrianakis,E.N.和Beer M.,在Base sequence determination in nucleic acids with theelectron microscope,III.Chemistry and microscopy of guanine-labeled DNA,PNAS53:564-71(1965)中所描述的。在该技术的一个实例中,使用可使用电子显微镜区分的金属标记物来标记单个DNA分子。然后将这些分子在平坦的表面上展开并使用电子显微镜成像以测量序列。
在测序之前,衔接子序列可附接到核酸分子,并且可以为感兴趣的特定序列富集核酸。序列富集可在衔接子序列附接之前或之后发生。
核酸分子或富集的核酸分子可附接到适用于本文公开的任何测序平台的任何测序衔接子。例如,序列衔接子可包括流通池序列、样品条形码或两者。在另一个实例中,序列衔接子可以是发夹形衔接子、Y形衔接子、叉状衔接子,和/或包含样品条形码。在某些情况下,衔接子不包含测序引物区域。在某些情况下,衔接子附着的DNA分子被扩增,并且扩增产物如本文所述富集特定序列。在某些情况下,在准备测序文库后,DNA分子富集了特定序列。衔接子可包含条形码序列。不同的条形码可以是至少1、2、3、4、5、6、7、8、9、10、11、12、13、14、15、16、17、18、19、20、21、22、23、24、25或更多个(或如遍及全文描述的任何长度)核酸碱基,例如7个碱基。条形码可以是随机序列、简并序列、半简并序列或定义的序列。在某些情况下,存在足够的条形码多样性,使得大体上(例如,至少70%、至少80%、至少90%或至少99%的)每种核酸分子以不同条形码序列加标签。在某些情况下,存在足够的条形码多样性,使得来自特定遗传基因座的大体上(例如,至少70%、至少80%、至少90%或至少99%的)每种核酸分子以不同条形码序列加标签。
测序衔接子可包含能够与一种或更多种测序引物杂交的序列。测序衔接子还可包含与固体支持物杂交的序列,例如流通池序列。例如,测序衔接子可以是流通池衔接子。该测序衔接子可附接至多核苷酸片段的一端或两端。在另一个实例中,测序衔接子可以是发夹形。例如,发夹形衔接子可包含互补双链部分和环部分,其中该双链部分可附接(例如,连接)至双链多核苷酸。发夹形测序衔接子可附接至多核苷酸片段的两端以生成环形分子,可对该环形分子多次测序。
在一些情况下,所述文库衔接子均不含有样品标识基序(或样品分子条形码)。可经由测序衔接子提供这样的样品标识基序。样品标识基序可包含至少4、5、6、7、8、9、10、20、30或40个核苷酸碱基的测序仪序列,该测序仪序列允许从来自其他样品的多核苷酸分子中鉴定来自给定样品的多核苷酸分子。例如,这可允许在相同池中对来自两个受试者的多核苷酸分子进行测序并随后鉴定所述受试者的序列读段。
测序仪基序包括使文库衔接子与测序系统偶联并对与该文库衔接子偶联的靶多核苷酸进行测序所需的核苷酸序列。该测序仪基序可包括与流通池序列互补的序列和可与用于在测序中使用的引物(或引发序列)选择性地杂交的序列(测序起始序列)。例如,这样的测序起始序列可与在合成测序(例如,Illumina)中使用的引物互补。这样的引物可包含在测序衔接子中。测序起始序列可以是引物杂交位点。
在一些情况下,所述文库衔接子均不含有完整的测序仪基序。所述文库衔接子可不含有或含有部分测序仪基序。在一些情况下,所述文库衔接子包含测序起始序列。所述文库衔接子可包含测序起始序列但不包含流通池序列。该测序起始序列可与用于测序的引物互补。该引物可以是序列特异性引物或通用引物。这样的测序起始序列可位于所述文库衔接子的单链部分上。作为替代方案,这样的测序起始序列可以是引发位点(例如,扭结(kink)或切口)以允许聚合酶在测序过程中与文库衔接子偶联。
衔接子可通过连接附着于DNA分子。在某些情况下,衔接子与双链DNA分子连接,使得每个衔接子将DNA分子的互补链不同地加标签。在一些情况下,衔接子序列可通过PCR附着,其中单链DNA的第一部分与靶序列互补,且第二部分包含衔接子序列。
可通过序列捕获方法进行对感兴趣的特定序列的富集。序列捕获可使用与感兴趣的靶杂交的固定的探针进行。序列捕获可使用附接于官能团例如生物素的探针进行,其允许探针与特定序列杂交以通过下拉(pulldown)从样品富集所述特定序列。在一些情况中,在与官能化的探针杂交之前,来自文库片段的特定序列诸如衔接子序列可通过使互补的、未官能化的多核苷酸序列与片段退火来掩蔽,以减少非特异性或脱靶结合。序列探针可靶向特定基因。序列捕获探针可靶向特定遗传基因座或基因。这样的基因可以是致癌基因。被捕获探针靶向的示例性基因包括图1所示的那些。带有点突变(SNV)的示例性基因包括但不限于,AKT1、ATM、CCNE1、CTNNB1、FGFR1、GNAS、JAK3、MLH1、NPM1、PTPN11、RIT1、TERT、ALK、BRAF、CDH1、EGFR、FGFR2、HNF1A、KIT MPL、NRAS、RAF1、ROS1、TP53、APC、BRCA1、CDK4、ERBB2、FGFR3、HRAS、KRAS、MYC、NTRK1、RB1、SMAD4、TSC1、AR、BRCA2、CDK6、ESR1、GATA3、IDH2、MAP2K2、NFE2L2、PIK3CA、RHEB、SRC、ARID1A、CCND2、CDKN2B、FBXW7、GNAQ、JAK2、MET、NOTCH1、PTEN、RHOA和STK11。带有拷贝数变异的示例性基因包括但不限于,AR、CCNE1、CDK6、ERBB2、FGFR2、KRAS、MYC、PIK3CA、BRAF、CDK4、EGFR、FGFR1、KIT、MET、PDGFRA和RAF1。具有基因融合的示例性基因包括但不限于:ALK、FGFR2、FGFR3、NTRK1、RET和ROS1。具有插入缺失的示例性基因包括但不限于:EGFR(例如,在外显子19和20处)、ERBB2(例如,在外显子19和20处)、和MET(例如,跳读外显子14)。示例性靶可包括CCND1和CCND2。序列捕获探针可平铺在一个基因上(例如,探针可靶向重叠区域)。序列探针可靶向非重叠区域。序列探针可对长度、解链温度和二级结构优化。
鸟嘌呤-胞嘧啶(GC)含量的定量量度
鸟嘌呤-胞嘧啶含量是DNA分子中为鸟嘌呤或胞嘧啶的含氮碱基的百分比。与遗传基因座的GC含量相关的定量量度可以是整个遗传基因座的GC含量。与遗传基因座的GC含量相关的定量量度可以是基因的外显子区域的GC含量。与遗传基因座的GC含量相关的定量量度可以是被映射到遗传基因座的读段覆盖的区域的GC含量。与GC含量相关的定量量度可以是对应于遗传基因座的序列捕获探针的GC含量。与遗传基因座的GC含量相关的定量量度可以是对应于遗传基因座的序列捕获探针的GC含量的中心趋势相关的量度。与中心趋势相关的量度可以是中心趋势的任何量度,如平均值、中位数或模式。与中心趋势相关的量度可以是中位数。给定区域的GC含量可通过将鸟苷和胞嘧啶碱基的数量除以该区域上的碱基总数来测量。
测序读段覆盖的定量量度
与测序读段覆盖相关的定量量度是指示来源于对应于遗传基因座(例如来自参考基因组的特定位置、碱基、区域、基因或染色体)的DNA分子的读段的数量的量度。为了将读段与遗传基因座相关联,读段可被映射于参考或与参考比对。进行映射或比对的软件(例如,Bowtie、BWA、mrsFAST、BLAST、BLAT)可将测序读段与遗传基因座相关联。在映射过程期间,可优化特定的参数。优化映射过程的非限制性示例可包括掩蔽重复区域;采用映射质量(例如,MAPQ)分数截止;使用不同的种子长度来生成比对;并限制基因组位置之间的编辑距离。
与测序读段覆盖相关的定量量度可包括与遗传基因座相关的读段的计数。在某些情况下,将计数转换为新的度量,以减轻不同测序深度、文库复杂度或遗传基因座大小的影响。示例性度量为每百万的每千碱基读段(Read Per Kilobase per Million,RPKM)、每百万的每千碱基片段(Fragments Per Kilobase per Million,FPKM)、M值的修剪平均值(TMM)、方差稳定原始计数和对数转换原始计数。可用于特定应用的其他转换对于本领域技术人员也是已知的。
可使用叠并读段来确定定量量度,其中每个叠并读段对应于初始模板DNA分子。叠并和定量读段家族的方法见于PCT/US2013/058061和PCT/US2014/000048中,其中的每一个通过引用整体并入本文。具体而言,可使用以下叠并方法,所述叠并方法使用条形码和来自测序读段的信息以将读段叠并为家族,使得每个家族共享条形码序列和至少一部分的测序读段序列。于是,对于大多数的家族,每个家族来源于单个初始模板DNA分子。从来自于家族的映射序列得到的计数可称为“独特分子计数”(UMC)。在某些情况下,确定与测序读段覆盖相关的定量量度包括通过与文库大小相关的度量将UMC归一化以提供归一化的UMC(“归一化的UMC”)。示例性方法是将遗传基因座的UMC除以所有UMC的总和;将遗传基因座的UMC除以所有常染色体UMC的总和。当比较多个测序读段数据集时,UMC可例如被两个测序读段数据集的遗传基因座的中值UMC归一化。在某些情况下,与测序读段覆盖相关的定量量度可以是归一化UMC,其如下进一步归一化:(i)从来源于训练样本的测序读段,确定相应的遗传基因座的归一化UMC;(ii)对于每个遗传基因座,样本的归一化UMC通过训练样本在相应基因座处的归一化UMC的中值归一化,由此提供遗传基因座的相对丰度(RA)。
共有序列可基于它们的序列来鉴定,例如通过基于前5个、10个、15个、20个或25个碱基内的相同序列叠并测序读段。在某些情况下,叠并允许另外相同的读段中的1个差异、2个差异、3个差异、4个差异或5个差异。在某些情况下,叠并使用读段的映射位置,例如测序读段的初始碱基的映射位置。在某些情况下,叠并使用条形码,且共享条形码序列的测序读段被叠并成共有序列。在某些情况下,叠并使用条形码和初始模板分子的序列二者。例如,共享条形码和映射到参考基因组中相同位置的所有读段可被叠并。在另一个实例中,可叠并共享条形码和初始模板分子的序列(或与初始模板分子的序列的同一性百分比)的所有读段。
在某些情况下,测序读段覆盖的定量量度是对基因组的特定子区域确定的。区域可以是箱(bin)、感兴趣的基因、外显子、对应于序列探针的区域、对应于引物扩增产物的区域、或对应于引物结合位点的区域。在某些情况下,基因组的子区域是对应于序列捕获探针的区域。如果读段的至少一部分映射到对应于序列捕获探针的区域的至少一部分,则读段可映射到对应于序列捕获探针的区域。如果读段的至少一部分映射到对应于序列捕获探针的区域的大部分,则读段可映射到对应于序列捕获探针的区域。如果读段的至少一部分映射跨对应于序列捕获探针的区域的中心点,则读段可映射到对应于序列捕获探针的区域。在某些情况下,与遗传基因座的测序读段覆盖相关的定量量度是遗传基因座内对应于基因组位置的探针的RA的中位数。例如,如果KRAS被三个探针覆盖,其具有2、3和5的RA,则遗传基因座的RA将是3。
“饱和平衡”校正
通常,本文描述的方法可用于增加核酸样品中变体调用(例如检测拷贝数变体)的特异性和灵敏度。例如,方法可减少数据样本中的噪声或失真的数量,减少检测到的假阳性变体的数量。当噪声和/或失真降低时,特异性和灵敏度增加。噪声可被认为是信号的不希望的随机增加。失真可被认为是信号或信号的一部分的幅值改变。
噪声可通过在拷贝和/或读取多核苷酸中的错误而引入。例如,在测序过程中,单个多核苷酸可首先经历扩增。扩增可能引入错误,致使扩增的多核苷酸的亚组可能在特定的基因座处包含与在该基因座处的原始碱基不同的碱基。此外,在读取过程中,在任何特定基因座处的碱基可能被不正确地读取。因此,序列读段的集合可能在基因座处包含一定百分比的与原始碱基不同的碱基判定。在通常的测序技术中,这种错误率可以是个位数,例如,2%-3%。当对全部被假定为具有相同序列的分子的集合进行测序时,这样的噪声是足够小的,使得人们可以以高可靠性鉴定原始碱基。
然而,如果亲本多核苷酸的集合包括具有在特定基因座处的序列变体的多核苷酸亚组,则噪声可能是一个显著的问题。例如,当无细胞DNA不仅包括种系DNA还包括来自另一个来源的DNA诸如胎儿DNA或来自癌细胞的DNA时,可能是这样的情况。在这种情况下,如果具有序列变体的分子的频率与通过测序过程引入的错误的频率在相同的范围内,则真序列变体可能无法与噪声区别。这可能会干扰,例如,样品中的序列变体的检测。
失真在测序过程中可表现为由亲本群体中相同频率的分子产生的信号强度的差异,例如,序列读段的总数。例如,可通过扩增偏倚、GC偏倚或测序偏倚引入失真。这可能会干扰样品中的拷贝数变异的检测。GC偏倚导致了在序列读取中GC含量丰富或贫乏区域的不均匀呈现。
本文公开的方法包括通过以下确定用于确定基线的遗传基因座的初始集:从数据集去除对其的测序读段覆盖相关的定量量度或测序读段覆盖相关的转化的定量量度与预测模型差异最大的那些遗传基因座(其在本文中可被称为去除高变异遗传基因座),从而提供剩余遗传基因座的第一集。在一些情况下,去除这些遗传基因座包括拟合将测序读段覆盖相关的定量量度与遗传基因座的GC含量相关的定量量度相关联的模型。例如,预测模型可将遗传基因座的RA与基因座的GC含量相关联。在某些情况下,预测模型是回归模型,包括非参数回归模型,例如LOESS和LOWESS回归模型。在某些情况下,基线定线通过去除与预测模型偏离最大的遗传基因座的5%、10%、15%、20%、25%、30%、35%、40%、45%、50%、55%、60%、65%或70%来进行。在某些情况下,基线定线通过去除与预测模型偏离最大的遗传基因座的至少5%、至少10%、至少15%、至少20%、至少25%、至少30%、至少35%、至少40%、至少45%、至少50%、至少55%、至少60%、至少65%或至少70%来进行。在某些情况下,通过测量遗传基因座相对于模型的残差来确定偏离。可以选择确切的截断值以提供排除从剩余遗传基因座的特定数量的变异。
PCT/US2014/072383中公开了确定来源于遗传基因座的来自样品的DNA链在测序读段内被代表的概率相关的定量量度的方法,其全部内容通过引用特此并入。确定定量量度可包括估计样品中存在的来源于基因座的初始模板DNA分子的数量。双链多核苷酸不生成序列读段的概率可以基于代表初始模板DNA分子的两条链的读段和仅代表初始模板DNA分子的单链的读段的相对数量来确定。
样品中未检测到的初始模板DNA分子的数量可以基于代表初始模板DNA分子的两条链的读段和仅代表初始模板DNA分子的单链的读段的相对数量来估计。作为实例,记录了特定遗传基因座即基因座A的计数,其中1000个分子是配对的(例如,两个链都被检测到)并且1000个分子是未配对的(例如,仅检测到单个链)。应注意,这里使用的术语“配对”和“未配对”与有时将其应用于测序读段以指示分子的两端还是单端被测序的这些术语不同。对于单个Watson或Crick链,假设使其通过转化后的过程具有统一的概率p,人们可以如下计算未能使其通过该过程的(未发现的)分子的比例:R,配对与未配对分子之比=1000/1000=1,因此R=1=p2/(2p(1-p))。这意味着p=2/3,并且丢失的分子的量等于(1-p)2=1/9。因此,在该实例中,大约11%的转化分子丢失并且从未被检测到。除了采用二项分布外,估算未发现的分子的数目的其他方法还包括基于观察到的序列读段的冗余的指数、β、γ或经验分布。在后一种情况下,配对和未配对分子的读出计数的分布可来源于这样的冗余,从而推断特定基因座处原始多核苷酸分子的潜在分布。这通常可导致更好地估算未发现的分子的数目。在某些情况下,p是来源于遗传基因座的来自样品的DNA链在测序读段内被代表的概率相关的定量量度。在一些情况下,p被类似地推导,但是使用不同的读段分布模型(例如,二项式、泊松、β、γ和负二项分布)。
对于测序读段覆盖的定量量度的转换可通过以下来确定:将来自去除了高变异遗传基因座的遗传基因座的集的定量量度或转换的测序读段覆盖与GC含量相关的定量量度和来源于遗传基因座的DNA链在测序读段内被代表的概率相关的定量量度相关联。在某些情况下,剩余的遗传基因座被假定为二倍体和/或以相同的拷贝数存在。在某些情况下,转换通过以下来确定:将剩余遗传基因座的测序读段覆盖相关的定量量度的中心趋势相关的量度与GC含量相关的定量量度和来源于遗传基因座的DNA链在测序读段内被代表的概率相关的定量量度拟合。转换可以,例如,(i)将去除高变异遗传基因座后剩余遗传基因座的测序读段覆盖的定量量度的中心趋势与GC含量相关的定量量度和来源于遗传基因座的DNA链在测序读段内被代表的概率相关的定量量度二者拟合。在某些情况下,剩余基因座的测序读段覆盖的定量量度的中心趋势相关的量度是剩余遗传基因座的UMC的中心趋势。在某些情况下,表面近似用于将剩余遗传基因座的UMC的表面或剩余遗传基因座的UMC的中心趋势与以下拟合:(i)GC含量相关的定量量度和(ii)来源于遗传基因座的DNA链在测序读段内被代表的概率相关的定量量度。例如,表面近似可以是初始模板DNA分子相关的量度(例如,UMC)与GC含量和p的定量量度的二维二次多项式表面拟合。在某些情况下,测序覆盖相关的转换的定量量度是从以下计算的基于以上确定的转换预期的值:(i)GC含量相关的定量量度和(ii)来源于遗传基因座的DNA链在测序读段内被代表的概率相关的定量量度。在某些情况下,测序覆盖相关的转换的定量量度是每个遗传基因座的残差(例如,基于表面近似的基因座的测序读段覆盖相关的预期定量量度与样品中遗传基因座的测序读段覆盖相关的观察到的定量量度的差异或商)。任选地,在确定了测序覆盖相关的转换的定量量度之后,可以基于测序读段覆盖的新转换的定量量度再次如上所述去除高变异遗传基因座。
“探针效率”校正
本文公开了使用参考样品确定和去除遗传基因座的偏倚的方法。在某些情况下,参考样品是来自无癌症的受试者的无细胞DNA的测序读段。在某些情况下,参考样品是来自具有在感兴趣的遗传基因座中基本上没有拷贝数变异的癌症细胞的受试者的无细胞DNA的测序读段。在某些情况下,参考样品是来自患有癌症的受试者的无细胞DNA的测序读段,其中怀疑已经历拷贝数变异的区域被排除在分析之外。在某些情况下,参考样品是来自无癌症的受试者的血浆样品。在某些情况下,参考样品是来自患有癌症的受试者的血浆样品。
参考样品的每个遗传基因座可以如以上“饱和平衡校正”中描述的处理,以提供测序读段覆盖的转换的定量量度。在某些情况下,测序覆盖相关的转换的定量量度是从以下计算的基于以上确定的转换预期的值:(i)GC含量相关的定量量度和(ii)来源于来自参考遗传基因座的遗传基因座的DNA链在测序读段内被代表的概率相关的定量量度。在某些情况下,测序覆盖相关的转换的定量量度是每个参考遗传基因座的残差(例如,基于表面近似的基因座的测序读段覆盖相关的预期定量量度与参考样品中遗传基因座的测序读段覆盖相关的观察到的定量量度的差异或商)。参考样品中遗传基因座的测序读段覆盖相关的转换的定量量度可被认为是遗传基因座的“效率”。例如,低效扩增的遗传基因座比非常有效地扩增的遗传基因座(在样品中以相同拷贝数存在)具有更低的UMC。
可以基于来自参考样品的遗传基因座确定的效率来校正与样品的测序读段覆盖相关的转换的定量量度。这种校正可以减少由从样品产生测序读段的过程引入到样品中的变异,这可以与连接效率、下拉效率、PCR效率、流通池聚类损失、解复用损失、叠并损失和比对损失相关。在一个实施方案中,校正包括将样品的测序覆盖的饱和后的转换的定量量度除以或减去与测序覆盖相关的预测的饱和后转换的定量量度。在某些情况下,与遗传基因座的测序覆盖有关的预测的饱和后转换的定量量度通过拟合来自样品的遗传基因座的测序覆盖相关的饱和后转换的定量量度与参考的测序读段覆盖相关的饱和后转换的定量量度之间的关系来确定。在一些情况下,拟合包括对参考的测序读段覆盖相关的饱和后转换的定量量度进行来自样品的遗传基因座的测序覆盖相关的饱和后转换的定量量度的局部回归(例如LOESS或LOWESS)或稳健线性回归。在某些情况下,拟合可以是线性回归、非线性回归或非参数回归。
任选地,来自探针效率校正的转换的定量量度可以是输入到“饱和平衡校正”转换中的输入,以产生变异减小的与测序读段覆盖相关的第三、进一步转换的定量量度。一般而言,为了进一步减少在测序读段覆盖的转换的定量量度内的变异,可以使用本文公开的任何方法来转换测序覆盖的转换的定量量度另外的次数。
基因水平概要
推断的拷贝数的基因水平概要可以基于如本文所公开的确定的测序读段覆盖的转换的定量量度来确定。通过丢弃高变异遗传基因座,可以推断相对于上述操作中选择的基线的拷贝数。例如,如果样品中剩余的遗传基因座被推断为二倍体,那么测序覆盖相关的转换的定量量度与基线不同的遗传基因座可以被推断为在肿瘤细胞中经历了拷贝数改变。在某些情况下,使用观察到的探针信号的基因水平中位数和使用基因中观察到的探针水平标准差计算的估计的标准差以及全基因组正常二倍体探针信号标准差,计算基因水平z-分数。
次等位基因频率基线优化
本文提供了使用测序读段中的变体的次等位基因频率来检测本文所述的拷贝数的基因水平概要中的错误和校正错误的方法。10%和90%之间、20%和80%之间、30%和70%之间、40%和60%之间、或约50%的来自无细胞体液的核酸的测序读段中存在的序列变体可以是受试者的种系序列中存在的杂合变体。在一些情况下,遗传基因座已被确定经历了如上所述的扩增。将变体的数量与推断的拷贝数进行比较以确定变体频率是否与推断的拷贝数不一致。在一个实例中,杂合遗传基因座可以在用于确定基线拷贝数的遗传基因座(例如,排除高变异遗传基因座后剩余的遗传基因座)中检查。在某些情况下,样品中大量的遗传基因座已被扩增,并且该基线可能被错误鉴定。在这种情况下,杂合性可能偏离1:1的比率,且检测并纠正不准确的基线定线。在第二个实例中,例如,基于与测序读段覆盖相关的转换的定量量度,可以推断遗传基因座以三倍体拷贝数存在。如果受试者的种系基因组具有带有遗传基因座的第一等位基因的一个染色体,并且第二染色体具有第二等位基因,则第一或第二等位基因可在癌细胞中已经重复。
Langmuir样饱和模型
不受理论束缚,本文公开了假设为基于探索历史临床数据以及涉及合成的加标模型系统的靶向实验的诱饵-cfDNA相互作用的调控机制的Langmuir样饱和模型。因此,在干扰测定效应(例如连接效率、PCR扩增偏差、测序假象等)的不存在下,诱饵下拉过程可被描述为
本描述中的K是诱饵效率,其依赖于诱饵序列特征以及其与基因组中邻近靶向的诱饵位置的DNA片段的相互作用。Isat是由下拉反应中有限的初始诱饵计数驱动的饱和参数,其是总诱饵池浓度以及复制计数的函数。本文使用的复制计数是指存在的序列捕获探针的相对或绝对量。例如,序列捕获阵列可以提供阵列上不同摩尔量的探针,以考虑不同的探针效率。图8说明了基于诱饵效率K和饱和参数Isat使真实拷贝数和独特分子计数相关联的模型。
诱饵效率K主要由GC含量驱动,而Isat由更复杂的诱饵耗尽机制和RNA二级结构相互作用驱动,可通过研究独特分子计数与总读段计数相互作用来粗略检查。除了非线性下拉反应外,探针信号还可以通过包含以下假设的乘法模型进一步建模:在幼稚模型下,cfDNA片段按基因组位置均匀分布,且随机采样过程是造成覆盖变异的支配因素。然后,可以通过将观察到的UMC与样品中的真实分子计数相关联来对拷贝数信号(例如,UMC)建模,考虑到潜在位置cfDNA特征(profile)、连接效率、下拉效率、PCR效率、流通池聚类损失、解复用损失、叠并损失和比对损失。
除了非线性下拉反应外,探针信号还可以通过简单乘法模型进一步建模,包括以下假设。在幼稚模型下,cfDNA片段按基因组位置均匀分布,且随机采样过程是造成覆盖变异的支配因素。然后,拷贝数信号,即与给定探针相关的读段计数可以被建模为:
观察到的UMC=真实UMC x潜在位置cfDNA特征(诱饵,cfDNA片段)x连接效率(位置、大小、cfDNA片段)x下拉效率(探针,cfDNA片段)x PCR效率(DNA片段)x流通池聚类损失x解复用损失和叠并损失x比对损失(cfDNA片段序列)。
这个模型假设了上述模型的乘性。可以通过一系列步骤,例如本文公开的基线确定方法,从观察到的UMC(例如,给定序列捕获探针的UMC)关于已建立的基线推断潜在的诱饵特异性拷贝数信号。
本文公开的方法提供了用于从样本集和训练集估算探针效率和诱饵饱和的方法。可选地,可以通过进行一组诱饵滴定实验来推断这样的参数,其中对每种探针观察到改变靶序列浓度对UMC的影响。如果K、Isat和UMC是已知的,则确定对应于未经历拷贝数变异的肿瘤细胞的UMC值或范围可以是可能的。例如,假设大部分遗传基因座没有经过拷贝数改变,观察到的UMC将主要来源于二倍体样品。已经经历拷贝数变异的样品将是其UMC落在探针的相应的K和Isat值的期望范围之外的那些遗传基因座。在某些情况下,例如,UMC值或范围将是依赖于每种探针的K和Isat的函数。例如,对应于二倍体拷贝数的UMC在两种探针之间可以是不同的。
计算机控制系统
本公开内容提供了被编程为实现本公开内容的方法的计算机控制系统。图12示出了被编程或以其他方式配置为实现本公开内容的方法的计算机系统1201。计算机系统1201包括中央处理单元(CPU,本文也被称为“处理器”和“计算机处理器”)1205,其可以是单核或多核处理器或用于并行处理的多于一个处理器。计算机系统1201还包括存储器或存储器位置1210(例如,随机存取存储器、只读存储器、闪速存储器)、电子存储单元1215(例如,硬盘)、用于与一个或更多个其他系统进行通信的通信界面1220(例如,网络适配器)和外围设备1225,诸如高速缓冲存储器、其他存储器、数据存储和/或电子显示适配器。存储器1210、储存单元1215、界面1220和外围设备1225与CPU 1205通过通信总线(实线),诸如主板(motherboard)通信。存储单元1215可以是用于存储数据的数据存储单元(或数据储存库)。计算机系统1201可借助于通信界面1220被可操作地耦合至计算机网络(“网络”)1230。网络1230可以是因特网(Internet)、互联网(internet)和/或外联网、或与因特网通信的内联网和/或外联网。在一些情况下,网络1230为电信和/或数据网络。网络1230可包括局部区域网络。网络1230可包括一个或更多个计算机服务器,这可支持分布式计算,诸如云计算。在一些情况下,借助于计算机系统1201,网络1230可实现对等网络(peer-to-peer network),其可使耦合至计算机系统1201的设备能够作为客户端或服务器运行。
CPU 1205可执行一系列的机器可读指令,该机器可读指令可以以程序或软件来体现。指令可被存储于存储器位置,诸如存储器1210中。指令可被导向CPU 1205,其可随后编程或以其他方式配置CPU 1205,以实现本公开内容的方法。由CPU 1205进行的操作的实例可包括读取、解码、执行和写回。
CPU 1205可以是电路诸如集成电路的一部分。系统1201的一个或更多个其他组件可被包含在该电路中。在一些情况下,电路为专用集成电路(ASIC)。
存储单元1215可存储文件,诸如驱动程序、库和保存的程序。存储单元1215可存储用户数据,例如,用户偏好和用户程序。在一些情况下,计算机系统1201可包括一个或更多个另外的数据存储单元,该数据存储单元在计算机系统1201的外部,诸如位于通过内联网或因特网而与计算机系统1201通信的远程服务器上。
计算机系统1201可与一个或更多个远程计算机系统通过网络1230进行通信。例如,计算机系统1201可与用户的远程计算机系统进行通信。远程计算机系统的实例包括个人计算机(例如便携式PC)、板型或平板PC(例如iPad、/>Galaxy Tab)、电话、智能电话(例如/>iPhone、Android支持的设备、/>)或个人数字助理。用户可经由网络1230访问计算机系统1201。
如本文描述的方法可通过机器(例如,计算机处理器)可执行代码的方式实现,该机器可执行代码被存储在计算机系统1201的电子存储位置,诸如,例如存储器1210或电子存储单元1215上。机器可执行代码或机器可读代码可以以软件的形式提供。在使用期间,代码可由处理器1205执行。在一些情况下,代码可从存储单元1215检索并存储在存储器1210上,以用于由处理器1205即时访问。在一些情况下,可排除电子存储单元1215,而将机器可执行指令存储于存储器1210中。
代码可被预编译并配置为用于与具有适用于执行该代码的处理器的机器一起使用,或者可在运行时间期间被编译。代码可以以编程语言的形式提供,该编程语言可被选择以便使得代码能够以预编译的或按编译原样(as-compiled)的方式被执行。
本文所提供的系统和方法的各方面,诸如计算机系统1201,可以以编程来体现。技术的多个方面可被认为是通常呈一种机器可读介质进行或体现的机器(或处理器)可执行代码和/或相关数据的形式的“产品”或“制品(articles of manufacture)”。机器可执行代码可被存储于电子存储单元诸如存储器(例如,只读存储器、随机存取存储器、闪速存储器)或硬盘上。“存储”型介质可包括计算机、处理器等的任何或所有有形存储器,或其相关模块,诸如多种半导体存储器、磁带驱动器、磁盘驱动器等,其可在任何时间为软件编程提供非暂时性存储。该软件的所有或部分有时可通过因特网或多种其他电信网络进行通信。例如,此类通信可使得将软件从一个计算机或处理器加载到另一个计算机或处理器,例如,从管理服务器或主机加载到应用服务器的计算机平台。因此,能够携带软件元件的另一类型的介质包括诸如在本地设备之间的物理界面、通过有线和光纤陆线网络以及在多种空中链路(air-links)上使用的光波、电波和电磁波。携带此类波的物理元件,诸如有线或无线链路、光链路等,也可被认为是携带软件的介质。如本文使用的,除非被限制为非暂时性的、有形的“储存”介质,否则术语诸如计算机或机器“可读介质”指参与将指令提供至处理器用于执行的任何介质。
因此,机器可读介质,诸如计算机可执行代码,可采取多种形式,包括但不限于有形存储介质、载波介质或物理传输介质。非易失性存储介质包括,例如光盘或磁盘,诸如在任何计算机等中的任何存储设备,诸如可用于实现如附图中示出的数据库等。易失性存储介质包括动态存储器,诸如此类计算机平台的主存储器。有形的传输介质包括同轴电缆;铜线和光纤,包括构成计算机系统内的总线的导线。载波传输介质可采取电信号或电磁信号或者声波或光波的形式,诸如在射频(RF)和红外(IR)数据通信期间生成的那些。因此,计算机可读介质的常见形式包括,例如:软盘(floppy disk)、软性磁盘(flexible disk)、硬盘、磁带、任何其他磁介质、CD-ROM、DVD或DVD-ROM、任何其他光学介质、穿孔卡片纸带、具有打孔模式的任何其他物理存储介质、RAM、ROM、PROM和EPROM、FLASH-EPROM、任何其他存储器芯片或盒、传输数据或指令的载波、传输此类载波的缆线或链路,或者计算机可从其读取编程代码和/或数据的任何其他介质。这些计算机可读介质的形式中的许多形式可参与向处理器传送一个或更多个指令的一个或更多个序列以用于执行。
计算机系统1201可包括电子显示器1235或与之通信,该电子显示器1235包括用户界面(UI)1240,用于提供例如报告。UI的实例包括但不限于图形用户界面(GUI)和基于网络的用户界面。
本公开内容的方法和系统可通过一个或更多个算法来实现。算法可在由中央处理器1205执行后通过软件来实现。
实施例
实施例1
检查先前生成的拷贝数变异加标数据揭示了在原始读段计数和UMC二者中的显著探针间信号变异,以及探针/基因水平拷贝数信号对潜在拷贝数变化的响应。参见图2。图3说明了三种基因(CCND1、CCND2和ERBB2)的推断拷贝数相比于理论拷贝数,证明了归一化的覆盖对样品中诱饵量的非线性响应。这些结果表明在下拉期间的诱饵耗竭,这通过在相同基因内在相邻探针中追踪诱饵滴定效应而证实,且独特分子计数有相当大的差异(从而对于具有高初始UMC的探针观察到更快的独特分子计数饱和)。
图4A说明,与每种探针相关的UMC关于探针p具有非线性响应。图4B说明,与每种探针相关的UMC关于探针GC含量具有非线性响应。
图5说明未进行饱和或探针效率校正的探针的UMC。图6显示饱和校正后的相同样品。图7显示探针效率校正后的相同样品。在每个阶段,基因组位置内的变异减少,导致肿瘤细胞潜在拷贝数变异的更清晰图像出现。图7中中位数探针探针效率校正后信号高于1.2的基因被称为已经发生了拷贝数变异。不同水平的探针效率校正后信号可能是由于肿瘤异质性或继发性肿瘤。
图9显示了饱和校正和探针效率校正后基线定线遗传基因座探针信号降噪的典型进程。
图10A示出了x轴上的参考样品的探针效率和来自肿瘤细胞中无拷贝数变异的受试者的样品的饱和校正后信号的图。关系是近似线性的。图10B说明来自肿瘤细胞中带有拷贝数变异的受试者的类似图。响应不像图10A那样线性。由通过确定来自参考样品的探针效率与基线定线遗传基因座(以黑色表示)的饱和校正后的UMC之间的关系来推断的预测的效率校正,将减少推定已经在肿瘤细胞中经历拷贝数扩增的遗传基因座(灰色圆点)中由于不同探针效率的变异。图11示出了基于饱和和探针效率校正后的UMC和MAF优化的基线定线,来自患者样品的拷贝数变异的示例性报告。星号表示指示属于在受试者的肿瘤细胞中已经经历拷贝数变异的遗传基因座的点。
实施例2
无细胞DNA从患有癌症的受试者获得,制备了条形码化的测序文库,通过用探针集序列捕获来富集一组致癌基因,并对条形码化的测序文库进行测序。将测序读段映射到参考基因组,并基于条形码序列和映射位置叠并成家族。对于对应于来自探针组的探针的中点的每个基因组坐标,计数跨越该中点的读段家族的数目以获得每探针UMC。对每种基因确定中位数每探针UMC。为了进行“饱和平衡校正”,将基因按其中位数每探针GC含量分组。去除其中位数每探针UMC与具有相似的中位数每探针GC含量的那些基因显著不同的基因。
对于每种探针,如本文所述确定p和GC含量。来自上一步的剩余的基因用于进行中位数基因水平UMC的二维二次多项式表面拟合以探测p和GC含量。将p和GC含量与预期的UMC相关联的函数用于确定预期的每探针UMC。通过将观察到的每探针UMC除以期望的每探针UMC来确定数据集的残差。每种探针的残留UMC是测序覆盖的转换的定量量度。
将基因再次按其中位数每探针GC含量分组,并且去除了其中位数每探针残留UMC与具有相似中位数每探针GC含量的基因显著不同的基因。“探针效率”校正然后如前述段落中所述的通过获得参考样品的残留UMC进行。然后将来自样品的每种探针的残留UMC除以来自参考的每种相应探针的残留UMC以获得探针效率校正后的UMC。
类似于以上的饱和平衡校正,剩余的基因用于进行探针效率校正后的UMC的二维二次多项式表面拟合以探测p和GC含量。将p和GC含量与预期的探针效率校正后的UMC相关联的函数用于确定预期的每探针探针效率校正后的UMC。通过将观察到的每探针探针效率校正后的UMC除以期望的每探针探针效率校正后的UMC来确定数据集的残差。每种探针的残留探针效率校正后的UMC是探针GC校正后的信号。
将剩余基因按其中位数每探针GC含量分组,并且去除了其中位数每探针GC校正后的信号与具有相似中位数每探针GC含量的那些基因显著不同的基因。
重复该实施例的过程,以探针GC校正后的信号作为起始输入代替初始UMC。
对于每种基因,使用探针GC校正后信号的中位数来总结每种基因。其中位数探针GC校正后信号与其他基因显著不同的基因被认为是肿瘤细胞中经历基因扩增或缺失的候选者。
对于每种基因,确定种系杂合等位基因并定量每种等位基因的相对频率。发现用于基线定线的遗传基因座具有的等位基因比率为约1:1,验证了基线定线遗传基因座的选择。
基于基因水平中位数探针GC校正后的信号和估计的从全基因组正常二倍体探针信号的标准差确定每种基因的Z-分数。Z-分数高于截止值的基因据报道在肿瘤细胞中已经经历了基因扩增。
实施例3
本文所述的方法通过在本公开内容的方法对比于对照方法中测量ERBB2拷贝数来验证。本公开内容的方法产生了观察到的拷贝数(CN)对比于理论拷贝数的线性响应,没有观察到正常(健康)队列中的假阳性CNV结果。参见图13,其显示推断的基因拷贝数相比于理论拷贝数估计值,其中实心圆点表示观察到的拷贝数~2(二倍体样品),开放圆点表示检测到的扩增事件且粗的水平虚线标记平均基因CN截止值。还参见图14,其描绘了图13的数据,其中对照数据由正方形表示。所有的CNV遵循预期的滴定趋势,降至2.15个拷贝。此外,由于变异减小,本公开内容的方法减少了数据中观察到的“噪声”,与对照方法相比,允许容易地区分CNV。参见图15的最右侧;三角形表示本公开内容的方法,而X表示对照方法。
虽然本文已经显示和描述了本发明的优选实施方案,但对于本领域技术人员将明显的是,此类实施方案仅通过示例的方式提供。并非意图将本发明限制于本说明书中提供的具体实例。虽然已参考以上提及的说明书描述了本发明,但本文实施方案的描述和说明并不意图以限制性的意义来解释。在不偏离本发明的情况下,本领域技术人员现在将想到许多变化、改变和替换。此外,应当理解,本发明的所有方面并不限于本文阐述的取决于多种条件和变量的具体描写、配置或相对比例。应当理解,在实践本发明时可采用本文描述的本发明的实施方案的各种替代选择。因此可预期,本发明还应涵盖任何此类的替代选择、修改、变化或等同物。以下权利要求意图界定本发明的范围,并且从而涵盖在这些权利要求范围内的方法和结构及其等同物。
Claims (36)
1.一种计算机系统,其包括存储器和处理器,其中所述存储器用于接收受试者的体液样品的无细胞DNA分子的测序读段;且所述处理器被配置为:
(a)从所述测序读段生成第一数据集,对于多个遗传基因座中的每个遗传基因座,所述第一数据集包含:(i)与测序读段覆盖相关的定量量度(“读段覆盖”);和(ii)所述遗传基因座的鸟嘌呤-胞嘧啶含量相关的定量量度(“GC含量”);
(b)从所述第一数据集去除高变异遗传基因座,其中去除包括:
(i)拟合涉及鸟嘌呤-胞嘧啶含量相关的定量量度和所述遗传基因座的测序读段覆盖的定量量度的模型;和
(ii)从所述第一数据集去除至少10%的遗传基因座,其中去除遗传基因座包括去除与所述模型差异最大的遗传基因座的至少10%,由此提供基线定线遗传基因座的第一数据集;
(c)通过进行饱和平衡校正和探针效率校正来校正所述第一数据集,其中进行饱和平衡校正包括通过以下将基线定线遗传基因座的所述第一数据集转换为饱和校正数据集:
(i)对于每个基线定线遗传基因座,从基线定线遗传基因座的第一数据集确定来源于所述基线定线遗传基因座的来自样品的DNA分子链在测序读段内被代表的概率相关的定量量度;
(ii)通过将基线定线遗传基因座的所述第一数据集中的读段覆盖与基线定线遗传基因座的所述第一数据集的GC含量和与基线定线遗传基因座的所述第一数据集中的来源于每个基因座的DNA链在所述测序读段内被代表的概率有关的定量量度相关联,确定所述读段覆盖的第一转换;和
(iii)将所述第一转换应用于来自基线定线遗传基因座的第一数据集的每个遗传基因座的读段覆盖以提供饱和校正数据集,其中所述饱和校正数据集包括基线定线遗传基因座的第一数据集的转换的读段覆盖的第一集;
并且其中进行探针效率校正包括通过以下转换所述饱和校正数据集为探针效率校正数据集:
(i)从所述饱和校正数据集去除相对于转换的读段覆盖的所述第一集为高变异遗传基因座的遗传基因座,从而提供基线定线遗传基因座的第二数据集;
(ii)确定与基线定线遗传基因座的所述第二数据集的探针效率相关的转换的读段覆盖的第一集的第二转换,其中所述探针效率通过对一种或更多种参考样品进行饱和平衡校正来确定,其中所述探针效率是通过进行所述饱和平衡校正获得的转换的读段覆盖;和
(iii)利用所述第二转换将基线定线遗传基因座的所述第二数据集的转换的读段覆盖的所述第一集进行转换,从而提供探针效率校正数据集,其中所述探针效率校正数据集包括基线定线遗传基因座的所述第二数据集的转换的读段覆盖的第二集;
(d)确定所述第一数据集的基线读段覆盖,其中所述基线读段覆盖与饱和平衡和探针效率有关;并
(e)确定所述多个遗传基因座中每个遗传基因座相对于所述基线读段覆盖的拷贝数状态。
2.根据权利要求1所述的计算机系统,其中所述处理器被配置为去除至少45%的所述遗传基因座。
3.根据权利要求1所述的计算机系统,其中确定所述第一转换包括(i)确定与基线定线遗传基因座的所述第一数据集的读段覆盖的中心趋势有关的量度;(ii)基于所述遗传基因座的GC含量和与来源于所述遗传基因座的DNA链在所述测序读段内被代表的概率相关的定量量度,确定拟合与基线定线遗传基因座的所述第一数据集的读段覆盖的中心趋势有关的量度的函数;和(iii)对于基线定线遗传基因座的所述第一数据集的每个遗传基因座,确定由所述函数预测的读段覆盖与所述读段覆盖的差异,其中所述差异为转换的读段覆盖。
4.根据权利要求3所述的计算机系统,其中所述函数是表面近似。
5.根据权利要求4所述的计算机系统,其中所述表面近似是二维二次多项式。
6.根据权利要求1所述的计算机系统,其中从所述第一数据集去除为高变异遗传基因座的遗传基因座包括:
(i)拟合涉及所述GC含量和所述饱和校正数据集的转换的读段覆盖的第一集的模型;和
(ii)从饱和校正数据集去除至少10%的遗传基因座,其中去除遗传基因座包括去除与所述模型差异最大的遗传基因座,由此提供基线定线遗传基因座的第二数据集。
7.根据权利要求6所述的计算机系统,其中所述处理器被配置为去除至少45%的所述遗传基因座。
8.根据权利要求1所述的计算机系统,其中所述一种或更多种参考样品是来自无癌症的受试者的无细胞体液样品。
9.根据权利要求1所述的计算机系统,其中所述一种或更多种参考样品是来自癌症受试者的无细胞体液样品,其中对应的遗传基因座未经历拷贝数改变。
10.根据权利要求1所述的计算机系统,其中确定所述第二转换包括(i)将对来自所述一种或更多种参考样品的遗传基因座确定的探针效率拟合到来自基线定线遗传基因座的第二数据集的读段覆盖的第一集;(ii)将基线定线遗传基因座的所述第二数据集的每个遗传基因座的转换的读段覆盖除以基于(i)的拟合的预测的探针效率。
11.根据权利要求1所述的计算机系统,其中所述处理器被配置为还:
(iv)通过将基线定线遗传基因座的所述第二数据集的转换的读段覆盖与基线定线遗传基因座的所述第二数据集的GC含量和与来源于基线定线遗传基因座的所述第二数据集中的每个基因座的DNA链在所述测序读段内被代表的概率有关的定量量度相关联,确定转换的读段覆盖的所述第二集的第三转换;
(v)将所述第三转换应用于转换的读段覆盖的第二集以提供第四数据集,其中所述第四数据集包括转换的定量读段覆盖的第三集。
12.根据权利要求1所述的计算机系统,其中所述体液样品的无细胞DNA分子或其扩增子使用与来自遗传基因座的集的遗传基因座的至少一部分互补的一种或更多种寡核苷酸探针富集遗传基因座的集。
13.根据权利要求12所述的计算机系统,其中来自遗传基因座的集的每个遗传基因座的GC含量是与来自遗传基因座的集的遗传基因座的至少一部分互补的所述一种或更多种寡核苷酸探针的鸟嘌呤-胞嘧啶含量的中心趋势相关的量度。
14.根据权利要求12所述的计算机系统,其中所述遗传基因座的读段覆盖是与对应于所述一种或更多种寡核苷酸探针的遗传基因座区域的读段覆盖的中心趋势相关的量度。
15.根据权利要求12所述的计算机系统,其中进行饱和平衡校正和进行探针效率校正包括拟合Langmuir模型,其中所述Langmuir模型包含K和Isat。
16.根据权利要求15所述的计算机系统,其中K和Isat对于所述一种或更多种寡核苷酸探针中的每种寡核苷酸探针根据经验确定。
17.根据权利要求16所述的计算机系统,其中进行饱和平衡校正和进行探针校正包括将所述遗传基因座的读段覆盖拟合于Langmuir模型,假设所述遗传基因座以相同的拷贝数状态存在,从而提供基线读段覆盖。
18.根据权利要求17所述的计算机系统,其中所述相同的拷贝数状态是二倍体。
19.根据权利要求17所述的计算机系统,其中所述基线读段覆盖是依赖于所述探针效率和所述饱和平衡的函数。
20.根据权利要求17所述的计算机系统,其中确定拷贝数状态包括比较所述遗传基因座的读段覆盖与所述基线读段覆盖。
21.根据权利要求8或9所述的计算机系统,其中所述无细胞体液选自由以下组成的组:血液、血清、血浆、尿和脑脊液。
22.根据权利要求1-20任一项所述的计算机系统,其中所述读段覆盖通过将所述测序读段映射于参考基因组来确定。
23.根据权利要求1-20任一项所述的计算机系统,其中获得所述测序读段包括将衔接子与来自受试者的体液的无细胞DNA分子连接。
24.根据权利要求23所述的计算机系统,其中所述无细胞DNA分子是双链DNA分子并且所述衔接子与所述双链DNA分子连接,使得每个衔接子将所述DNA分子的互补链不同地加标签以提供加标签链。
25.根据权利要求24所述的计算机系统,其中确定来源于遗传基因座的DNA链在所述测序读段内被代表的概率相关的定量量度包括将测序读段分选为配对读段和非配对读段,其中(i)每个配对读段对应于所述集中由来源于双链多核苷酸分子的第一加标签链和第二差异加标签互补链生成的测序读段,并且(ii)每个非配对读段代表第一加标签链,并且来源于双链多核苷酸分子的第二差异加标签互补链在所述测序读段之中没有被代表。
26.根据权利要求25所述的计算机系统,其中所述处理器被配置为还确定映射到一个或更多个遗传基因座中的每一个的(i)所述配对读段和(ii)所述非配对读段的定量量度,以基于映射到每个基因座的配对读段和非配对读段相关的所述定量量度确定所述样品中映射到所述一个或更多个遗传基因座中的每一个的全部双链DNA分子相关的定量量度。
27.根据权利要求23所述的计算机系统,其中所述衔接子包含条形码序列。
28.根据权利要求22所述的计算机系统,其中确定所述读段覆盖包括基于所述测序读段映射到所述参考基因组的位置和条形码序列叠并所述测序读段。
29.根据权利要求1-20任一项所述的计算机系统,其中所述遗传基因座包括一种或更多种致癌基因。
30.根据权利要求1-20任一项所述的计算机系统,其中所述处理器被配置为还通过确定受试者的种系基因组对其为杂合的基线定线遗传基因座内的变体的相对量,确定受试者的肿瘤细胞中基线定线遗传基因座的至少子集已经经历拷贝数改变。
31.根据权利要求30所述的计算机系统,其中所述变体的相对量不是近似相等的。
32.根据权利要求31所述的计算机系统,其中将所述变体的相对量对其不是近似相等的基线定线遗传基因座从基线定线遗传基因座去除,从而提供等位基因频率校正的基线定线遗传基因座。
33.根据权利要求32所述的计算机系统,其中所述等位基因频率校正的基线定线遗传基因座在前述权利要求的任一项的计算机系统中用作基线定线基因座。
34.一种系统,所述系统包括:
(a)网络;
(b)数据库,所述数据库被连接到所述网络,包括配置成存储核酸序列数据的计算机存储器;
(c)生物信息学计算机,其包括计算机存储器和一个或更多个计算机处理器,所述计算机被连接到所述网络;
其中所述计算机还包括机器可执行代码,当由所述一个或更多个计算机处理器执行时,所述机器可执行代码复制存储在所述数据库中的所述核酸序列数据,将复制的数据写入所述生物信息学计算机中的存储器,并进行包括以下的步骤:
(i)从所述核酸序列数据生成第一数据集,对于多个遗传基因座中的每个遗传基因座,所述第一数据集包含:(a)与测序读段覆盖相关的定量量度(“读段覆盖”);和(b)所述遗传基因座的鸟嘌呤-胞嘧啶含量相关的定量量度(“GC含量”);
(ii)从所述第一数据集去除高变异遗传基因座,其中去除包括:
(a)拟合涉及鸟嘌呤-胞嘧啶含量相关的定量量度和所述遗传基因座的测序读段覆盖的定量量度的模型;和
(b)从所述第一数据集去除至少10%的遗传基因座,其中去除遗传基因座包括去除与所述模型差异最大的遗传基因座的至少10%,由此提供基线定线遗传基因座的第一数据集;
(iii)通过进行饱和平衡校正和探针效率校正来校正所述第一数据集,其中进行饱和平衡校正包括通过以下将基线定线遗传基因座的所述第一数据集转换为饱和校正数据集:
(a)对于每个基线定线遗传基因座,从基线定线遗传基因座的第一数据集确定来源于所述基线定线遗传基因座的来自样品的DNA分子链在测序读段内被代表的概率相关的定量量度;
(b)通过将基线定线遗传基因座的所述第一数据集中的读段覆盖与基线定线遗传基因座的所述第一数据集的GC含量和与基线定线遗传基因座的所述第一数据集中的来源于每个基因座的DNA链在所述测序读段内被代表的概率有关的定量量度相关联,
确定所述读段覆盖的第一转换;和
(c)将所述第一转换应用于来自基线定线遗传基因座的第一数据集的每个遗传基因座的读段覆盖以提供饱和校正数据集,其中所述饱和校正数据集包括基线定线遗传基因座的第一数据集的转换的读段覆盖的第一集;
并且其中进行探针效率校正包括通过以下转换所述饱和校正数据集为探针效率校正数据集:
(a)从所述饱和校正数据集去除相对于转换的读段覆盖的所述第一集为高变异遗传基因座的遗传基因座,从而提供基线定线遗传基因座的第二数据集;
(b)确定与基线定线遗传基因座的所述第二数据集的探针效率相关的转换的读段覆盖的第一集的第二转换,其中所述探针效率通过对一种或更多种参考样品进行饱和平衡校正来确定,其中
所述探针效率是通过进行所述饱和平衡校正获得的转换的读段覆盖;和
(c)利用所述第二转换将基线定线遗传基因座的所述第二数据集的转换的读段覆盖的所述第一集进行转换,从而提供探针效率校正数据集,其中所述探针效率校正数据集包括基线定线遗传基因座的所述第二数据集的转换的读段覆盖的第二集;
(iv)确定对于所述第一数据集的基线读段覆盖,其中所述基线读段覆盖与饱和平衡和探针效率有关;并
(v)确定所述多个遗传基因座中每个遗传基因座相对于所述基线读段覆盖的拷贝数状态。
35.根据权利要求34所述的系统,其中所述数据库连接于核酸测序仪。
36.一种非暂时性计算机可读介质,其包括机器可执行代码,当由一个或多个计算机处理器执行时实施包括以下的方法:
(a)获得受试者的体液样品的无细胞DNA分子的测序读段;
(b)从所述测序读段生成第一数据集,对于多个遗传基因座中的每个遗传基因座,所述第一数据集包含:(i)与测序读段覆盖相关的定量量度(“读段覆盖”);和(ii)所述遗传基因座的鸟嘌呤-胞嘧啶含量相关的定量量度(“GC含量”);
(c)从所述第一数据集去除高变异遗传基因座,其中去除包括:
(i)拟合涉及鸟嘌呤-胞嘧啶含量相关的定量量度和所述遗传基因座的测序读段覆盖的定量量度的模型;和
(ii)从所述第一数据集去除至少10%的遗传基因座,其中去除遗传基因座包括去除与所述模型差异最大的遗传基因座的至少10%,由此提供基线定线遗传基因座的第一数据集;
(d)通过进行饱和平衡校正和探针效率校正来校正所述第一数据集,其中进行饱和平衡校正包括通过以下将基线定线遗传基因座的所述第一数据集转换为饱和校正数据集:
(i)对于每个基线定线遗传基因座,从基线定线遗传基因座的第一数据集确定来源于所述基线定线遗传基因座的来自样品的DNA分子链在测序读段内被代表的概率相关的定量量度;
(ii)通过将基线定线遗传基因座的所述第一数据集中的读段覆盖与基线定线遗传基因座的所述第一数据集的GC含量和与基线定线遗传基因座的所述第一数据集中的来源于每个基因座的DNA链在所述测序读段内被代表的概率有关的定量量度相关联,确定所述读段覆盖的第一转换;和
(iii)将所述第一转换应用于来自基线定线遗传基因座的第一数据集的每个遗传基因座的读段覆盖以提供饱和校正数据集,其中所述饱和校正数据集包括基线定线遗传基因座的第一数据集的转换的读段覆盖的第一集;
并且其中进行探针效率校正包括通过以下转换所述饱和校正数据集为探针效率校正数据集:
(i)从所述饱和校正数据集去除相对于转换的读段覆盖的所述第一集为高变异遗传基因座的遗传基因座,从而提供基线定线遗传基因座的第二数据集;
(ii)确定与基线定线遗传基因座的所述第二数据集的探针效率相关的转换的读段覆盖的第一集的第二转换,其中所述探针效率通过对一种或更多种参考样品进行饱和平衡校正来确定,其中所述探针效率是通过进行所述饱和平衡校正获得的转换的读段覆盖;和
(iii)利用所述第二转换将基线定线遗传基因座的所述第二数据集的转换的读段覆盖的所述第一集进行转换,从而提供探针效率校正数据集,其中所述探针效率校正数据集包括基线定线遗传基因座的所述第二数据集的转换的读段覆盖的第二集;
(e)确定所述第一数据集的基线读段覆盖,其中所述基线读段覆盖与饱和平衡和探针效率有关;并
(f)确定所述多个遗传基因座中每个遗传基因座相对于所述基线读段覆盖的拷贝数状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311131468.XA CN117174167A (zh) | 2015-12-17 | 2016-12-16 | 通过分析无细胞dna确定肿瘤基因拷贝数的方法 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562269051P | 2015-12-17 | 2015-12-17 | |
US62/269,051 | 2015-12-17 | ||
PCT/US2016/067356 WO2017106768A1 (en) | 2015-12-17 | 2016-12-16 | Methods to determine tumor gene copy number by analysis of cell-free dna |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311131468.XA Division CN117174167A (zh) | 2015-12-17 | 2016-12-16 | 通过分析无细胞dna确定肿瘤基因拷贝数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108603228A CN108603228A (zh) | 2018-09-28 |
CN108603228B true CN108603228B (zh) | 2023-09-01 |
Family
ID=59057742
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311131468.XA Pending CN117174167A (zh) | 2015-12-17 | 2016-12-16 | 通过分析无细胞dna确定肿瘤基因拷贝数的方法 |
CN201680081723.6A Active CN108603228B (zh) | 2015-12-17 | 2016-12-16 | 通过分析无细胞dna确定肿瘤基因拷贝数的方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311131468.XA Pending CN117174167A (zh) | 2015-12-17 | 2016-12-16 | 通过分析无细胞dna确定肿瘤基因拷贝数的方法 |
Country Status (7)
Country | Link |
---|---|
US (3) | US20170240973A1 (zh) |
EP (1) | EP3390668A4 (zh) |
JP (3) | JP2019507585A (zh) |
CN (2) | CN117174167A (zh) |
CA (1) | CA3008651A1 (zh) |
SG (1) | SG11201805119QA (zh) |
WO (1) | WO2017106768A1 (zh) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012142213A2 (en) | 2011-04-15 | 2012-10-18 | The Johns Hopkins University | Safe sequencing system |
US9892230B2 (en) | 2012-03-08 | 2018-02-13 | The Chinese University Of Hong Kong | Size-based analysis of fetal or tumor DNA fraction in plasma |
EP3447495B2 (en) | 2012-10-29 | 2024-03-13 | The Johns Hopkins University | Papanicolaou test for ovarian and endometrial cancers |
US10364467B2 (en) | 2015-01-13 | 2019-07-30 | The Chinese University Of Hong Kong | Using size and number aberrations in plasma DNA for detecting cancer |
US11286531B2 (en) | 2015-08-11 | 2022-03-29 | The Johns Hopkins University | Assaying ovarian cyst fluid |
US20190287645A1 (en) * | 2016-07-06 | 2019-09-19 | Guardant Health, Inc. | Methods for fragmentome profiling of cell-free nucleic acids |
CA3027919C (en) | 2016-09-30 | 2023-02-28 | Guardant Health, Inc. | Methods for multi-resolution analysis of cell-free nucleic acids |
US9850523B1 (en) | 2016-09-30 | 2017-12-26 | Guardant Health, Inc. | Methods for multi-resolution analysis of cell-free nucleic acids |
US20180225413A1 (en) * | 2016-12-22 | 2018-08-09 | Grail, Inc. | Base Coverage Normalization and Use Thereof in Detecting Copy Number Variation |
MX2020001575A (es) | 2017-08-07 | 2020-11-18 | Univ Johns Hopkins | Materiales y métodos para evaluar y tratar el cáncer. |
CN111278993A (zh) * | 2017-09-15 | 2020-06-12 | 加利福尼亚大学董事会 | 从无细胞核酸中检测体细胞单核苷酸变体并应用于微小残留病变监测 |
US11482303B2 (en) | 2018-06-01 | 2022-10-25 | Grail, Llc | Convolutional neural network systems and methods for data classification |
WO2020069350A1 (en) | 2018-09-27 | 2020-04-02 | Grail, Inc. | Methylation markers and targeted methylation probe panel |
EP3874060A1 (en) * | 2018-10-31 | 2021-09-08 | Guardant Health, Inc. | Methods, compositions and systems for calibrating epigenetic partitioning assays |
CN110232951B (zh) * | 2018-12-06 | 2023-08-01 | 苏州金唯智生物科技有限公司 | 判断测序数据饱和的方法、计算机可读介质和应用 |
US11581062B2 (en) | 2018-12-10 | 2023-02-14 | Grail, Llc | Systems and methods for classifying patients with respect to multiple cancer classes |
CN113661249A (zh) | 2019-01-31 | 2021-11-16 | 夸登特健康公司 | 用于分离无细胞dna的组合物和方法 |
WO2020186024A1 (en) * | 2019-03-13 | 2020-09-17 | Grail, Inc. | Systems and methods for enriching for cancer-derived fragments using fragment size |
US11475981B2 (en) | 2020-02-18 | 2022-10-18 | Tempus Labs, Inc. | Methods and systems for dynamic variant thresholding in a liquid biopsy assay |
US11211147B2 (en) | 2020-02-18 | 2021-12-28 | Tempus Labs, Inc. | Estimation of circulating tumor fraction using off-target reads of targeted-panel sequencing |
US11211144B2 (en) | 2020-02-18 | 2021-12-28 | Tempus Labs, Inc. | Methods and systems for refining copy number variation in a liquid biopsy assay |
WO2023282916A1 (en) | 2021-07-09 | 2023-01-12 | Guardant Health, Inc. | Methods of detecting genomic rearrangements using cell free nucleic acids |
JP2023540221A (ja) | 2020-08-25 | 2023-09-22 | ガーダント ヘルス, インコーポレイテッド | バリアントの起源を予測するための方法およびシステム |
EP4305200A1 (en) | 2021-03-09 | 2024-01-17 | Guardant Health, Inc. | Detecting the presence of a tumor based on off-target polynucleotide sequencing data |
WO2023023402A2 (en) | 2021-08-20 | 2023-02-23 | Guardant Health, Inc. | Methods for simultaneous molecular and sample barcoding |
WO2023150633A2 (en) | 2022-02-02 | 2023-08-10 | Guardant Health, Inc. | Multifunctional primers for paired sequencing reads |
US20240052419A1 (en) * | 2022-08-10 | 2024-02-15 | Guardant Health, Inc. | Methods and systems for detecting genetic variants |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103374518A (zh) * | 2012-04-12 | 2013-10-30 | 维里纳塔健康公司 | 拷贝数变异的检测和分类 |
WO2013166517A1 (en) * | 2012-05-04 | 2013-11-07 | Complete Genomics, Inc. | Methods for determining absolute genome-wide copy number variations of complex tumors |
CN103562407A (zh) * | 2011-02-09 | 2014-02-05 | 伯乐生命医学产品有限公司 | 核酸的分析 |
WO2014039556A1 (en) * | 2012-09-04 | 2014-03-13 | Guardant Health, Inc. | Systems and methods to detect rare mutations and copy number variation |
WO2014165596A1 (en) * | 2013-04-03 | 2014-10-09 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
WO2015100427A1 (en) * | 2013-12-28 | 2015-07-02 | Guardant Health, Inc. | Methods and systems for detecting genetic variants |
WO2015184404A1 (en) * | 2014-05-30 | 2015-12-03 | Verinata Health, Inc. | Detecting fetal sub-chromosomal aneuploidies and copy number variations |
Family Cites Families (280)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4725536A (en) | 1985-09-19 | 1988-02-16 | Genetics Institute, Inc. | Reagent polynucleotide complex with multiple target binding regions, and kit and methods |
US6150517A (en) | 1986-11-24 | 2000-11-21 | Gen-Probe | Methods for making oligonucleotide probes for the detection and/or quantitation of non-viral organisms |
US5149625A (en) | 1987-08-11 | 1992-09-22 | President And Fellows Of Harvard College | Multiplex analysis of DNA |
US4942124A (en) | 1987-08-11 | 1990-07-17 | President And Fellows Of Harvard College | Multiplex sequencing |
US5656731A (en) | 1987-10-15 | 1997-08-12 | Chiron Corporation | Nucleic acid-amplified immunoassay probes |
US5124246A (en) | 1987-10-15 | 1992-06-23 | Chiron Corporation | Nucleic acid multimers and amplified nucleic acid hybridization assays using same |
US5424186A (en) | 1989-06-07 | 1995-06-13 | Affymax Technologies N.V. | Very large scale immobilized polymer synthesis |
US5871928A (en) | 1989-06-07 | 1999-02-16 | Fodor; Stephen P. A. | Methods for nucleic acid analysis |
US6040138A (en) | 1995-09-15 | 2000-03-21 | Affymetrix, Inc. | Expression monitoring by hybridization to high density oligonucleotide arrays |
US5143854A (en) | 1989-06-07 | 1992-09-01 | Affymax Technologies N.V. | Large scale photolithographic solid phase synthesis of polypeptides and receptor binding screening thereof |
US5925525A (en) | 1989-06-07 | 1999-07-20 | Affymetrix, Inc. | Method of identifying nucleotide differences |
US6309822B1 (en) | 1989-06-07 | 2001-10-30 | Affymetrix, Inc. | Method for comparing copy number of nucleic acid sequences |
US6551784B2 (en) | 1989-06-07 | 2003-04-22 | Affymetrix Inc | Method of comparing nucleic acid sequences |
US5744101A (en) | 1989-06-07 | 1998-04-28 | Affymax Technologies N.V. | Photolabile nucleoside protecting groups |
US5800992A (en) | 1989-06-07 | 1998-09-01 | Fodor; Stephen P.A. | Method of detecting nucleic acids |
US5200314A (en) | 1990-03-23 | 1993-04-06 | Chiron Corporation | Polynucleotide capture assay employing in vitro amplification |
EP0562047A4 (en) | 1990-12-06 | 1995-11-02 | Affymax Tech Nv | Sequencing by hybridization of a target nucleic acid to a matrix of defined oligonucleotides |
US6582908B2 (en) | 1990-12-06 | 2003-06-24 | Affymetrix, Inc. | Oligonucleotides |
US5981179A (en) | 1991-11-14 | 1999-11-09 | Digene Diagnostics, Inc. | Continuous amplification reaction |
US5424413A (en) | 1992-01-22 | 1995-06-13 | Gen-Probe Incorporated | Branched nucleic acid probes |
US5573905A (en) | 1992-03-30 | 1996-11-12 | The Scripps Research Institute | Encoded combinatorial chemical libraries |
US6020124A (en) | 1992-04-27 | 2000-02-01 | Trustees Of Dartmouth College | Detection of soluble gene sequences in biological fluids |
US5981176A (en) | 1992-06-17 | 1999-11-09 | City Of Hope | Method of detecting and discriminating between nucleic acid sequences |
WO1995000530A1 (en) | 1993-06-25 | 1995-01-05 | Affymax Technologies N.V. | Hybridization and sequencing of nucleic acids |
US5500356A (en) | 1993-08-10 | 1996-03-19 | Life Technologies, Inc. | Method of nucleic acid sequence selection |
US6309823B1 (en) | 1993-10-26 | 2001-10-30 | Affymetrix, Inc. | Arrays of nucleic acid probes for analyzing biotransformation genes and methods of using the same |
US5681697A (en) | 1993-12-08 | 1997-10-28 | Chiron Corporation | Solution phase nucleic acid sandwich assays having reduced background noise and kits therefor |
CH686982A5 (fr) | 1993-12-16 | 1996-08-15 | Maurice Stroun | Méthode pour le diagnostic de cancers. |
US20030017081A1 (en) | 1994-02-10 | 2003-01-23 | Affymetrix, Inc. | Method and apparatus for imaging a sample on a device |
US5714330A (en) | 1994-04-04 | 1998-02-03 | Lynx Therapeutics, Inc. | DNA sequencing by stepwise ligation and cleavage |
US5695934A (en) | 1994-10-13 | 1997-12-09 | Lynx Therapeutics, Inc. | Massively parallel sequencing of sorted polynucleotides |
US6406848B1 (en) | 1997-05-23 | 2002-06-18 | Lynx Therapeutics, Inc. | Planar arrays of microparticle-bound polynucleotides |
US5604097A (en) | 1994-10-13 | 1997-02-18 | Spectragen, Inc. | Methods for sorting polynucleotides using oligonucleotide tags |
US6013445A (en) | 1996-06-06 | 2000-01-11 | Lynx Therapeutics, Inc. | Massively parallel signature sequencing by ligation of encoded adaptors |
US5846719A (en) | 1994-10-13 | 1998-12-08 | Lynx Therapeutics, Inc. | Oligonucleotide tags for sorting and identification |
US6600996B2 (en) | 1994-10-21 | 2003-07-29 | Affymetrix, Inc. | Computer-aided techniques for analyzing biological sequences |
DE69535240T2 (de) | 1994-10-28 | 2007-06-06 | Gen-Probe Inc., San Diego | Zusammensetzungen und Verfahren für die gleichzeitige Detektion und Quantifizierung von einer Mehrheit spezifischer Nuklein Säure Sequenzen |
US5750341A (en) | 1995-04-17 | 1998-05-12 | Lynx Therapeutics, Inc. | DNA sequencing by parallel oligonucleotide extensions |
US5648245A (en) | 1995-05-09 | 1997-07-15 | Carnegie Institution Of Washington | Method for constructing an oligonucleotide concatamer library by rolling circle replication |
US5968740A (en) | 1995-07-24 | 1999-10-19 | Affymetrix, Inc. | Method of Identifying a Base in a Nucleic Acid |
GB9516636D0 (en) | 1995-08-14 | 1995-10-18 | Univ London | In-situ nucleic acid amplification and detection |
US5763175A (en) | 1995-11-17 | 1998-06-09 | Lynx Therapeutics, Inc. | Simultaneous sequencing of tagged polynucleotides |
US5854033A (en) | 1995-11-21 | 1998-12-29 | Yale University | Rolling circle replication reporter systems |
US6156504A (en) | 1996-03-15 | 2000-12-05 | The Penn State Research Foundation | Detection of extracellular tumor-associated nucleic acid in blood plasma or serum using nucleic acid amplification assays |
EP0938320B2 (en) | 1996-03-26 | 2014-06-18 | Michael S. Kopreski | Method enabling use of extracellular rna extracted from plasma or serum to detect, monitor or evaluate cancer |
US6458530B1 (en) | 1996-04-04 | 2002-10-01 | Affymetrix Inc. | Selecting tag nucleic acids |
US6300077B1 (en) | 1996-08-14 | 2001-10-09 | Exact Sciences Corporation | Methods for the detection of nucleic acids |
GB9620209D0 (en) | 1996-09-27 | 1996-11-13 | Cemu Bioteknik Ab | Method of sequencing DNA |
US5935793A (en) | 1996-09-27 | 1999-08-10 | The Chinese University Of Hong Kong | Parallel polynucleotide sequencing method using tagged primers |
US6124092A (en) | 1996-10-04 | 2000-09-26 | The Perkin-Elmer Corporation | Multiplex polynucleotide capture methods and compositions |
US6117631A (en) | 1996-10-29 | 2000-09-12 | Polyprobe, Inc. | Detection of antigens via oligonucleotide antibody conjugates |
US6046005A (en) | 1997-01-15 | 2000-04-04 | Incyte Pharmaceuticals, Inc. | Nucleic acid sequencing with solid phase capturable terminators comprising a cleavable linking group |
WO1999028505A1 (en) | 1997-12-03 | 1999-06-10 | Curagen Corporation | Methods and devices for measuring differential gene expression |
US6054276A (en) | 1998-02-23 | 2000-04-25 | Macevicz; Stephen C. | DNA restriction site mapping |
US6787308B2 (en) | 1998-07-30 | 2004-09-07 | Solexa Ltd. | Arrayed biomolecules and their use in sequencing |
AU5584999A (en) | 1998-08-28 | 2000-03-21 | Invitrogen Corporation | System for the rapid manipulation of nucleic acid sequences |
US6653077B1 (en) | 1998-09-04 | 2003-11-25 | Lynx Therapeutics, Inc. | Method of screening for genetic polymorphism |
US6503718B2 (en) | 1999-01-10 | 2003-01-07 | Exact Sciences Corporation | Methods for detecting mutations using primer extension for detecting disease |
GB9901475D0 (en) | 1999-01-22 | 1999-03-17 | Pyrosequencing Ab | A method of DNA sequencing |
JP2002535999A (ja) | 1999-02-05 | 2002-10-29 | アマシャム バイオサイエンス ユーケイ リミテッド | ゲノム分析方法 |
US6629040B1 (en) | 1999-03-19 | 2003-09-30 | University Of Washington | Isotope distribution encoded tags for protein identification |
WO2000058516A2 (en) | 1999-03-26 | 2000-10-05 | Whitehead Institute For Biomedical Research | Universal arrays |
US6964846B1 (en) | 1999-04-09 | 2005-11-15 | Exact Sciences Corporation | Methods for detecting nucleic acids indicative of cancer |
EP1046717B1 (en) | 1999-04-20 | 2010-10-06 | National Institute of Advanced Industrial Science and Technology | Method and probes for determining a concentration of target nucleic acid molecules and method for analyzing data obtained by the method |
US6355431B1 (en) | 1999-04-20 | 2002-03-12 | Illumina, Inc. | Detection of nucleic acid amplification reactions using bead arrays |
US6242186B1 (en) | 1999-06-01 | 2001-06-05 | Oy Jurilab Ltd. | Method for detecting a risk of cancer and coronary heart disease and kit therefor |
US6818395B1 (en) | 1999-06-28 | 2004-11-16 | California Institute Of Technology | Methods and apparatus for analyzing polynucleotide sequences |
US6326148B1 (en) | 1999-07-12 | 2001-12-04 | The Regents Of The University Of California | Detection of copy number changes in colon cancer |
US6440706B1 (en) | 1999-08-02 | 2002-08-27 | Johns Hopkins University | Digital amplification |
US6849403B1 (en) | 1999-09-08 | 2005-02-01 | Exact Sciences Corporation | Apparatus and method for drug screening |
US6586177B1 (en) | 1999-09-08 | 2003-07-01 | Exact Sciences Corporation | Methods for disease detection |
WO2001023610A2 (en) | 1999-09-29 | 2001-04-05 | Solexa Ltd. | Polynucleotide sequencing |
AU4714201A (en) | 1999-12-07 | 2001-06-18 | Exact Sciences Corporation | Supracolonic aerodigestive neoplasm detection |
US6489114B2 (en) | 1999-12-17 | 2002-12-03 | Bio Merieux | Process for labeling a ribonucleic acid, and labeled RNA fragments which are obtained thereby |
WO2001057269A2 (en) | 2000-02-07 | 2001-08-09 | Illumina, Inc. | Nucleic acid detection methods using universal priming |
US20020072058A1 (en) | 2000-03-24 | 2002-06-13 | Voelker Leroy L. | Method for amplifying quinolone-resistance-determining-regions and identifying polymorphic variants thereof |
US20030207300A1 (en) | 2000-04-28 | 2003-11-06 | Matray Tracy J. | Multiplex analytical platform using molecular tags |
EP1158055A1 (fr) | 2000-05-26 | 2001-11-28 | Xu Qi University of Teaxs Laboratoire de Leucémie Chen | Méthode pour le diagnostic de cancers |
AU2002246612B2 (en) | 2000-10-24 | 2007-11-01 | The Board Of Trustees Of The Leland Stanford Junior University | Direct multiplex characterization of genomic DNA |
US20020142345A1 (en) | 2000-12-22 | 2002-10-03 | Nelsen Anita J. | Methods for encoding and decoding complex mixtures in arrayed assays |
US20030049616A1 (en) | 2001-01-08 | 2003-03-13 | Sydney Brenner | Enzymatic synthesis of oligonucleotide tags |
JP5073888B2 (ja) | 2001-03-28 | 2012-11-14 | 花王株式会社 | 静電荷像現像用トナー |
US6849404B2 (en) | 2001-05-07 | 2005-02-01 | Bioneer Corporation | Polymerase chain reaction of DNA of which base sequence is completely unidentified |
US7406385B2 (en) | 2001-10-25 | 2008-07-29 | Applera Corporation | System and method for consensus-calling with per-base quality values for sample assemblies |
ATE312946T1 (de) | 2002-03-05 | 2005-12-15 | Epigenomics Ag | Verfahren und vorrichtung zur bestimmung der gewebespezifität von freier dna in körperflüssigkeiten |
US20030186251A1 (en) | 2002-04-01 | 2003-10-02 | Brookhaven Science Associates, Llc | Genome sequence tags |
US7727720B2 (en) | 2002-05-08 | 2010-06-01 | Ravgen, Inc. | Methods for detection of genetic disorders |
US7822555B2 (en) | 2002-11-11 | 2010-10-26 | Affymetrix, Inc. | Methods for identifying DNA copy number changes |
US10229244B2 (en) | 2002-11-11 | 2019-03-12 | Affymetrix, Inc. | Methods for identifying DNA copy number changes using hidden markov model based estimations |
EP1578994A2 (en) | 2002-11-11 | 2005-09-28 | Affymetrix, Inc. | Methods for identifying dna copy number changes |
US7704687B2 (en) | 2002-11-15 | 2010-04-27 | The Johns Hopkins University | Digital karyotyping |
EP1606417A2 (en) | 2003-03-07 | 2005-12-21 | Rubicon Genomics Inc. | In vitro dna immortalization and whole genome amplification using libraries generated from randomly fragmented dna |
US20040259118A1 (en) | 2003-06-23 | 2004-12-23 | Macevicz Stephen C. | Methods and compositions for nucleic acid sequence analysis |
JP5183063B2 (ja) | 2003-07-05 | 2013-04-17 | ザ ジョンズ ホプキンス ユニバーシティ | 遺伝的変異の検出および列挙のための方法ならびに組成物 |
ATE435301T1 (de) | 2003-10-16 | 2009-07-15 | Sequenom Inc | Nicht invasiver nachweis fötaler genetischer merkmale |
DE10348407A1 (de) | 2003-10-17 | 2005-05-19 | Widschwendter, Martin, Prof. | Prognostische und diagnostische Marker für Zell-proliferative Erkrankungen von Brustgeweben |
US20070111233A1 (en) | 2003-10-30 | 2007-05-17 | Bianchi Diana W | Prenatal diagnosis using cell-free fetal DNA in amniotic fluid |
US20060024681A1 (en) | 2003-10-31 | 2006-02-02 | Agencourt Bioscience Corporation | Methods for producing a paired tag from a nucleic acid sequence and methods of use thereof |
JP2007524410A (ja) | 2004-01-23 | 2007-08-30 | リングヴィテ エーエス | ポリヌクレオチドライゲーション反応の改良 |
US7217522B2 (en) | 2004-02-12 | 2007-05-15 | Campass Genetics Llc | Genetic analysis by sequence-specific sorting |
US20100216153A1 (en) | 2004-02-27 | 2010-08-26 | Helicos Biosciences Corporation | Methods for detecting fetal nucleic acids and diagnosing fetal abnormalities |
US20060046258A1 (en) | 2004-02-27 | 2006-03-02 | Lapidus Stanley N | Applications of single molecule sequencing |
WO2005111242A2 (en) | 2004-05-10 | 2005-11-24 | Parallele Bioscience, Inc. | Digital profiling of polynucleotide populations |
US7276720B2 (en) | 2004-07-19 | 2007-10-02 | Helicos Biosciences Corporation | Apparatus and methods for analyzing samples |
US20060035258A1 (en) | 2004-08-06 | 2006-02-16 | Affymetrix, Inc. | Methods for identifying DNA copy number changes |
US7937225B2 (en) | 2004-09-03 | 2011-05-03 | New York University | Systems, methods and software arrangements for detection of genome copy number variation |
US20060073506A1 (en) | 2004-09-17 | 2006-04-06 | Affymetrix, Inc. | Methods for identifying biological samples |
US9109256B2 (en) | 2004-10-27 | 2015-08-18 | Esoterix Genetic Laboratories, Llc | Method for monitoring disease progression or recurrence |
US7424371B2 (en) | 2004-12-21 | 2008-09-09 | Helicos Biosciences Corporation | Nucleic acid analysis |
US7393665B2 (en) | 2005-02-10 | 2008-07-01 | Population Genetics Technologies Ltd | Methods and compositions for tagging and identifying polynucleotides |
ITRM20050068A1 (it) | 2005-02-17 | 2006-08-18 | Istituto Naz Per Le Malattie I | Metodo per la rivelazione di acidi nucleici di agenti patogeni batterici o di parassiti nelle urine. |
WO2006099604A2 (en) | 2005-03-16 | 2006-09-21 | Compass Genetics, Llc | Methods and compositions for assay readouts on multiple analytical platforms |
EP1861512A4 (en) | 2005-03-18 | 2009-12-09 | Fluidigm Corp | THERMAL REACTION DEVICE AND USE METHOD THEREFOR |
ES2313143T3 (es) | 2005-04-06 | 2009-03-01 | Maurice Stroun | Metodo para el diagnostico de cancer mediante la deteccion de adn y arn circulantes. |
WO2007145612A1 (en) | 2005-06-06 | 2007-12-21 | 454 Life Sciences Corporation | Paired end sequencing |
US20070020640A1 (en) | 2005-07-21 | 2007-01-25 | Mccloskey Megan L | Molecular encoding of nucleic acid templates for PCR and other forms of sequence analysis |
US7666593B2 (en) | 2005-08-26 | 2010-02-23 | Helicos Biosciences Corporation | Single molecule sequencing of captured nucleic acids |
DK1929039T4 (en) | 2005-09-29 | 2014-02-17 | Keygene Nv | High throughput-screening af mutageniserede populationer |
US7329860B2 (en) | 2005-11-23 | 2008-02-12 | Illumina, Inc. | Confocal imaging methods and apparatus |
US7537897B2 (en) | 2006-01-23 | 2009-05-26 | Population Genetics Technologies, Ltd. | Molecular counting |
US20070172839A1 (en) | 2006-01-24 | 2007-07-26 | Smith Douglas R | Asymmetrical adapters and methods of use thereof |
US8383338B2 (en) | 2006-04-24 | 2013-02-26 | Roche Nimblegen, Inc. | Methods and systems for uniform enrichment of genomic regions |
US7702468B2 (en) | 2006-05-03 | 2010-04-20 | Population Diagnostics, Inc. | Evaluating genetic disorders |
RU2008146868A (ru) | 2006-05-18 | 2010-06-27 | Кэрис МПИ, Инк.445 Норт Фифс Стрит, 3-ий Флор, Феникс, Аризона 85004, США (US) | Система и способ определения персонализированого медицинского вмешательства при болезненном состоянии |
US20080090239A1 (en) | 2006-06-14 | 2008-04-17 | Daniel Shoemaker | Rare cell analysis using sample splitting and dna tags |
FR2904833A1 (fr) | 2006-08-11 | 2008-02-15 | Bioquanta Sarl | Procede de dosage d'acide nuclieque par fluorescence |
US7754429B2 (en) | 2006-10-06 | 2010-07-13 | Illumina Cambridge Limited | Method for pair-wise sequencing a plurity of target polynucleotides |
ES2679996T3 (es) | 2006-11-15 | 2018-09-03 | Biospherex Llc | Secuenciación multi-etiqueta y análisis ecogenómico |
WO2008070144A2 (en) | 2006-12-06 | 2008-06-12 | Duke University | Imprinted genes and disease |
US8262900B2 (en) | 2006-12-14 | 2012-09-11 | Life Technologies Corporation | Methods and apparatus for measuring analytes using large scale FET arrays |
ES2609094T3 (es) | 2007-01-25 | 2017-04-18 | Dana-Farber Cancer Institute, Inc. | Uso de anticuerpos anti-EGFR en el tratamiento de enfermedades mediadas por un EGFR mutante |
WO2008092155A2 (en) | 2007-01-26 | 2008-07-31 | Illumina, Inc. | Image data efficient genetic sequencing method and system |
CL2008000717A1 (es) | 2007-03-13 | 2008-09-22 | Amgen Inc | Metodo para pronosticar si un paciente sera no respondedor al tratamiento con un agente de union especifica a un polipeptido de egfr que comprende determinar la presencia o ausencia de una mutacion de k-ras en un tumor del paciente. |
WO2008148072A2 (en) | 2007-05-24 | 2008-12-04 | The Brigham And Women's Hospital, Inc. | Disease-associated genetic variations and methods for obtaining and using same |
EP2164985A4 (en) | 2007-06-01 | 2014-05-14 | 454 Life Sciences Corp | SYSTEM AND METHOD FOR IDENTIFYING INDIVIDUAL SAMPLES FROM A MULTIPLEX MIXTURE |
WO2008154317A1 (en) | 2007-06-06 | 2008-12-18 | Pacific Biosciences Of California, Inc. | Methods and processes for calling bases in sequence by incorporation methods |
PT2183693E (pt) | 2007-07-23 | 2014-01-14 | Univ Hong Kong Chinese | Diagnóstico de aneuploidia cromossómica fetal utilizando sequenciação genómica |
US20100112590A1 (en) | 2007-07-23 | 2010-05-06 | The Chinese University Of Hong Kong | Diagnosing Fetal Chromosomal Aneuploidy Using Genomic Sequencing With Enrichment |
US20090053719A1 (en) | 2007-08-03 | 2009-02-26 | The Chinese University Of Hong Kong | Analysis of nucleic acids by digital pcr |
AU2008295992B2 (en) | 2007-09-07 | 2014-04-17 | Fluidigm Corporation | Copy number variation determination, methods and systems |
CN102007407A (zh) | 2007-11-21 | 2011-04-06 | 考斯摩斯德公司 | 基因组鉴定系统 |
US20110171640A1 (en) | 2008-02-12 | 2011-07-14 | Novartis Vaccines And Diagnostics, Inc. | Method for isolating cell free apoptotic or fetal nucleic acids |
US8216789B2 (en) | 2008-02-27 | 2012-07-10 | University Of Washington | Diagnostic panel of cancer antibodies and methods for use |
US8206926B2 (en) | 2008-03-26 | 2012-06-26 | Sequenom, Inc. | Restriction endonuclease enhanced polymorphic sequence detection |
BRPI0909212A2 (pt) | 2008-03-28 | 2015-08-18 | Pacific Biosciences California | Composições e método para sequeciamento de ácido nucléico |
US20110160290A1 (en) | 2008-05-21 | 2011-06-30 | Muneesh Tewari | Use of extracellular rna to measure disease |
US20090298709A1 (en) | 2008-05-28 | 2009-12-03 | Affymetrix, Inc. | Assays for determining telomere length and repeated sequence copy number |
DE102008025656B4 (de) | 2008-05-28 | 2016-07-28 | Genxpro Gmbh | Verfahren zur quantitativen Analyse von Nukleinsäuren, Marker dafür und deren Verwendung |
US20100041048A1 (en) | 2008-07-31 | 2010-02-18 | The Johns Hopkins University | Circulating Mutant DNA to Assess Tumor Dynamics |
US20100062494A1 (en) | 2008-08-08 | 2010-03-11 | President And Fellows Of Harvard College | Enzymatic oligonucleotide pre-adenylation |
WO2010021936A1 (en) | 2008-08-16 | 2010-02-25 | The Board Of Trustees Of The Leland Stanford Junior University | Digital pcr calibration for high throughput sequencing |
EP3216874A1 (en) | 2008-09-05 | 2017-09-13 | TOMA Biosciences, Inc. | Methods for stratifying and annotating cancer drug treatment options |
US8383345B2 (en) | 2008-09-12 | 2013-02-26 | University Of Washington | Sequence tag directed subassembly of short sequencing reads into long sequencing reads |
PT2562268T (pt) | 2008-09-20 | 2017-03-29 | Univ Leland Stanford Junior | Diagnóstico não invasivo de aneuploidia fetal por sequenciação |
US8546128B2 (en) | 2008-10-22 | 2013-10-01 | Life Technologies Corporation | Fluidics system for sequential delivery of reagents |
US20100301398A1 (en) | 2009-05-29 | 2010-12-02 | Ion Torrent Systems Incorporated | Methods and apparatus for measuring analytes |
EP2379748A4 (en) | 2008-12-23 | 2012-08-29 | Illumina Inc | MULTIBASE RELEASE FOR LONG READINGS IN SEQUENCING BY SYNTHESIS PROTOCOLS |
US20100323348A1 (en) | 2009-01-31 | 2010-12-23 | The Regents Of The University Of Colorado, A Body Corporate | Methods and Compositions for Using Error-Detecting and/or Error-Correcting Barcodes in Nucleic Acid Amplification Process |
US20120165202A1 (en) | 2009-04-30 | 2012-06-28 | Good Start Genetics, Inc. | Methods and compositions for evaluating genetic markers |
US9085798B2 (en) | 2009-04-30 | 2015-07-21 | Prognosys Biosciences, Inc. | Nucleic acid constructs and methods of use |
US8574835B2 (en) | 2009-05-29 | 2013-11-05 | Life Technologies Corporation | Scaffolded nucleic acid polymer particles and methods of making and using |
US8673627B2 (en) | 2009-05-29 | 2014-03-18 | Life Technologies Corporation | Apparatus and methods for performing electrochemical reactions |
US20130143747A1 (en) | 2011-12-05 | 2013-06-06 | Myriad Genetics, Incorporated | Methods of detecting cancer |
US9524369B2 (en) | 2009-06-15 | 2016-12-20 | Complete Genomics, Inc. | Processing and analysis of complex nucleic acid sequence data |
RU2014144463A (ru) | 2009-06-25 | 2015-06-20 | Фред Хатчинсон Кансэр Рисёч Сентер | Способ измерения адаптивного иммунитета |
US10017812B2 (en) | 2010-05-18 | 2018-07-10 | Natera, Inc. | Methods for non-invasive prenatal ploidy calling |
US20120220478A1 (en) | 2009-07-20 | 2012-08-30 | Bar Harbor Biotechnology, Inc. | Methods for assessing disease risk |
EP2494065B1 (en) | 2009-10-26 | 2015-12-23 | Lifecodexx AG | Means and methods for non-invasive diagnosis of chromosomal aneuploidy |
EP2499262A4 (en) | 2009-11-12 | 2015-01-07 | Esoterix Genetic Lab Llc | ANALYSIS OF THE NUMBER OF COPIES OF GENETIC LOCUS |
US9023769B2 (en) | 2009-11-30 | 2015-05-05 | Complete Genomics, Inc. | cDNA library for nucleic acid sequencing |
US9752187B2 (en) | 2009-12-11 | 2017-09-05 | Nucleix | Categorization of DNA samples |
US8835358B2 (en) | 2009-12-15 | 2014-09-16 | Cellular Research, Inc. | Digital counting of individual molecules by stochastic attachment of diverse labels |
US9315857B2 (en) | 2009-12-15 | 2016-04-19 | Cellular Research, Inc. | Digital counting of individual molecules by stochastic attachment of diverse label-tags |
DK2516680T3 (en) | 2009-12-22 | 2016-05-02 | Sequenom Inc | Method and kits to identify aneuploidy |
US10662474B2 (en) | 2010-01-19 | 2020-05-26 | Verinata Health, Inc. | Identification of polymorphic sequences in mixtures of genomic DNA by whole genome sequencing |
DK3382037T3 (da) | 2010-01-19 | 2021-05-25 | Verinata Health Inc | Fremgangsmåder til bestemmelse af fraktionen af føtale nukleinsyrer i maternelle prøver |
US20120100548A1 (en) * | 2010-10-26 | 2012-04-26 | Verinata Health, Inc. | Method for determining copy number variations |
WO2011090556A1 (en) | 2010-01-19 | 2011-07-28 | Verinata Health, Inc. | Methods for determining fraction of fetal nucleic acid in maternal samples |
US9260745B2 (en) | 2010-01-19 | 2016-02-16 | Verinata Health, Inc. | Detecting and classifying copy number variation |
US10388403B2 (en) | 2010-01-19 | 2019-08-20 | Verinata Health, Inc. | Analyzing copy number variation in the detection of cancer |
US20110177512A1 (en) | 2010-01-19 | 2011-07-21 | Predictive Biosciences, Inc. | Method for assuring amplification of an abnormal nucleic acid in a sample |
EP2536854B1 (en) | 2010-02-18 | 2017-07-19 | The Johns Hopkins University | Personalized tumor biomarkers |
EP2547698B1 (en) | 2010-03-14 | 2015-07-29 | The Translational Genomics Research Institute | Methods of determining susceptibility of tumors to tyrosine kinase inhibitors |
CN101967517B (zh) | 2010-03-19 | 2012-11-07 | 黄乐群 | 一种无需借助pcr的基因检测方法 |
US10047397B2 (en) | 2010-04-16 | 2018-08-14 | Chronix Biomedical | Breast cancer associated circulating nucleic acid biomarkers |
WO2011140510A2 (en) | 2010-05-06 | 2011-11-10 | Bioo Scientific Corporation | Oligonucleotide ligation, barcoding and methods and compositions for improving data quality and throughput using massively parallel sequencing |
PL2576837T3 (pl) | 2010-06-04 | 2018-04-30 | Chronix Biomedical | Krążeniowe biomarkery typu kwasu nukleinowego związane z rakiem prostaty |
US9080210B2 (en) | 2010-06-09 | 2015-07-14 | Keygene N.V. | High throughput screening using combinatorial sequence barcodes |
EP2400035A1 (en) | 2010-06-28 | 2011-12-28 | Technische Universität München | Methods and compositions for diagnosing gastrointestinal stromal tumors |
WO2012006291A2 (en) | 2010-07-06 | 2012-01-12 | Life Technologies Corporation | Systems and methods to detect copy number variation |
SG10201505723UA (en) | 2010-07-23 | 2015-09-29 | Harvard College | Methods for detecting signatures of disease or conditions in bodily fluids |
CN103140288B (zh) | 2010-07-29 | 2015-03-11 | Toto株式会社 | 光催化剂涂装体和光催化剂涂覆液 |
EP2601609B1 (en) | 2010-08-02 | 2017-05-17 | Population Bio, Inc. | Compositions and methods for discovery of causative mutations in genetic disorders |
US20120034603A1 (en) | 2010-08-06 | 2012-02-09 | Tandem Diagnostics, Inc. | Ligation-based detection of genetic variants |
US11031095B2 (en) | 2010-08-06 | 2021-06-08 | Ariosa Diagnostics, Inc. | Assay systems for determination of fetal copy number variation |
EP2426217A1 (en) | 2010-09-03 | 2012-03-07 | Centre National de la Recherche Scientifique (CNRS) | Analytical methods for cell free nucleic acids and applications |
AU2011299088B2 (en) | 2010-09-09 | 2015-12-24 | Traxxsson, Llc | Combination methods of diagnosing cancer in a patient |
ES2690753T3 (es) | 2010-09-21 | 2018-11-22 | Agilent Technologies, Inc. | Aumento de la confianza en las identificaciones de alelos con el recuento molecular |
WO2012042374A2 (en) | 2010-10-01 | 2012-04-05 | Anssi Jussi Nikolai Taipale | Method of determining number or concentration of molecules |
WO2012048341A1 (en) | 2010-10-08 | 2012-04-12 | President And Fellows Of Harvard College | High-throughput single cell barcoding |
US8725422B2 (en) | 2010-10-13 | 2014-05-13 | Complete Genomics, Inc. | Methods for estimating genome-wide copy number variations |
EP3461914A1 (en) | 2010-10-22 | 2019-04-03 | Cold Spring Harbor Laboratory | Varietal counting of nucleic acids for obtaining genomic copy number information |
WO2012066451A1 (en) | 2010-11-15 | 2012-05-24 | Pfizer Inc. | Prognostic and predictive gene signature for colon cancer |
KR102185244B1 (ko) | 2010-11-30 | 2020-12-02 | 더 차이니즈 유니버시티 오브 홍콩 | 암과 연관된 유전적 또는 분자적 이상들의 검출 |
US20120184449A1 (en) | 2010-12-23 | 2012-07-19 | Sequenom, Inc. | Fetal genetic variation detection |
US9163281B2 (en) | 2010-12-23 | 2015-10-20 | Good Start Genetics, Inc. | Methods for maintaining the integrity and identification of a nucleic acid template in a multiplex sequencing reaction |
KR20210131432A (ko) | 2010-12-30 | 2021-11-02 | 파운데이션 메디신 인코포레이티드 | 종양 샘플의 다유전자 분석의 최적화 |
WO2012097053A1 (en) | 2011-01-11 | 2012-07-19 | Via Genomes, Inc. | Methods, systems, databases, kits and arrays for screening for and predicting the risk of and identifying the presence of tumors and cancers |
WO2012103031A2 (en) | 2011-01-25 | 2012-08-02 | Ariosa Diagnostics, Inc. | Detection of genetic abnormalities |
US20120238464A1 (en) | 2011-03-18 | 2012-09-20 | Baylor Research Institute | Biomarkers for Predicting the Recurrence of Colorectal Cancer Metastasis |
WO2012129363A2 (en) | 2011-03-24 | 2012-09-27 | President And Fellows Of Harvard College | Single cell nucleic acid detection and analysis |
EP2691544B1 (en) | 2011-03-30 | 2017-09-13 | Verinata Health, Inc | Method for verifying bioassay samples |
EP2670866A4 (en) | 2011-04-05 | 2015-09-02 | Translational Genomics Res Inst | BIOMARKERS AND METHODS OF USE |
WO2014014498A1 (en) | 2012-07-20 | 2014-01-23 | Verinata Health, Inc. | Detecting and classifying copy number variation in a fetal genome |
WO2012142213A2 (en) | 2011-04-15 | 2012-10-18 | The Johns Hopkins University | Safe sequencing system |
US9411937B2 (en) | 2011-04-15 | 2016-08-09 | Verinata Health, Inc. | Detecting and classifying copy number variation |
CA2834291A1 (en) | 2011-04-25 | 2012-11-01 | Biorad Laboratories, Inc. | Methods and compositions for nucleic acid analysis |
US8697408B2 (en) | 2011-05-06 | 2014-04-15 | New England Biolabs, Inc. | Ligation enhancement |
CN103890245B (zh) | 2011-05-20 | 2020-11-17 | 富鲁达公司 | 核酸编码反应 |
US9752176B2 (en) | 2011-06-15 | 2017-09-05 | Ginkgo Bioworks, Inc. | Methods for preparative in vitro cloning |
KR101454886B1 (ko) | 2011-08-01 | 2014-11-03 | 주식회사 셀레믹스 | 핵산분자의 제조방법 |
US10704164B2 (en) | 2011-08-31 | 2020-07-07 | Life Technologies Corporation | Methods, systems, computer readable media, and kits for sample identification |
WO2013033721A1 (en) | 2011-09-02 | 2013-03-07 | Atreca, Inc. | Dna barcodes for multiplexed sequencing |
US8712697B2 (en) | 2011-09-07 | 2014-04-29 | Ariosa Diagnostics, Inc. | Determination of copy number variations using binomial probability calculations |
US20130079241A1 (en) | 2011-09-15 | 2013-03-28 | Jianhua Luo | Methods for Diagnosing Prostate Cancer and Predicting Prostate Cancer Relapse |
WO2013052907A2 (en) | 2011-10-06 | 2013-04-11 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
US9367663B2 (en) | 2011-10-06 | 2016-06-14 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
US10424394B2 (en) | 2011-10-06 | 2019-09-24 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
US20130102485A1 (en) | 2011-10-19 | 2013-04-25 | Inhan Lee | Method of Determining a Diseased State in a Subject |
EP2768985B1 (en) | 2011-10-21 | 2019-03-20 | Chronix Biomedical | Colorectal cancer associated circulating nucleic acid biomarkers |
NO3051026T3 (zh) | 2011-10-21 | 2018-07-28 | ||
US20130122499A1 (en) | 2011-11-14 | 2013-05-16 | Viomics, Inc. | System and method of detecting local copy number variation in dna samples |
WO2013086352A1 (en) | 2011-12-07 | 2013-06-13 | Chronix Biomedical | Prostate cancer associated circulating nucleic acid biomarkers |
KR101768652B1 (ko) | 2011-12-08 | 2017-08-16 | 파이브3 제노믹스, 엘엘씨 | Mdm2-포함 이중 소염색체들 및 그의 방법들 |
US20130184165A1 (en) | 2012-01-13 | 2013-07-18 | Data2Bio | Genotyping by next-generation sequencing |
ES2665071T3 (es) | 2012-02-17 | 2018-04-24 | Fred Hutchinson Cancer Research Center | Composiciones y métodos para identificar mutaciones de manera precisa |
US11177020B2 (en) | 2012-02-27 | 2021-11-16 | The University Of North Carolina At Chapel Hill | Methods and uses for molecular tags |
ES2663234T3 (es) | 2012-02-27 | 2018-04-11 | Cellular Research, Inc | Composiciones y kits para recuento molecular |
WO2013128281A1 (en) | 2012-02-28 | 2013-09-06 | Population Genetics Technologies Ltd | Method for attaching a counter sequence to a nucleic acid sample |
US9890429B2 (en) | 2012-02-29 | 2018-02-13 | Dana-Farber Cancer Institute, Inc. | Compositions, kits, and methods for the identification, assessment, prevention, and therapy of cancer |
US9892230B2 (en) | 2012-03-08 | 2018-02-13 | The Chinese University Of Hong Kong | Size-based analysis of fetal or tumor DNA fraction in plasma |
EP2825675B1 (en) | 2012-03-13 | 2017-12-27 | Patel, Abhijit Ajit | Measurement of nucleic acid variants using highly-multiplexed error-suppressed deep sequencing |
WO2013142213A1 (en) | 2012-03-20 | 2013-09-26 | Wake Forest University Health Sciences | Methods, systems, and computer readable media for tracking and verifying receipt of contents of a delivery within an organization |
EP4234713A3 (en) | 2012-03-20 | 2024-02-14 | University Of Washington Through Its Center For Commercialization | Methods of lowering the error rate of massively parallel dna sequencing using duplex consensus sequencing |
EP4239081A3 (en) | 2012-03-26 | 2023-11-08 | The Johns Hopkins University | Rapid aneuploidy detection |
US8209130B1 (en) | 2012-04-04 | 2012-06-26 | Good Start Genetics, Inc. | Sequence assembly |
WO2013159035A2 (en) | 2012-04-19 | 2013-10-24 | Medical College Of Wisconsin, Inc. | Highly sensitive surveillance using detection of cell free dna |
CA2873585C (en) | 2012-05-14 | 2021-11-09 | Cb Biotechnologies, Inc. | Method for increasing accuracy in quantitative detection of polynucleotides |
WO2013181170A1 (en) | 2012-05-31 | 2013-12-05 | Board Of Regents, The University Of Texas System | Method for accurate sequencing of dna |
EP2859123A4 (en) | 2012-06-11 | 2015-12-16 | Sequenta Inc | METHOD OF SEQUENCE DETERMINATION USING SEQUENCE TAGS |
US11261494B2 (en) | 2012-06-21 | 2022-03-01 | The Chinese University Of Hong Kong | Method of measuring a fractional concentration of tumor DNA |
WO2014004726A1 (en) | 2012-06-26 | 2014-01-03 | Caifu Chen | Methods, compositions and kits for the diagnosis, prognosis and monitoring of cancer |
US20160040229A1 (en) | 2013-08-16 | 2016-02-11 | Guardant Health, Inc. | Systems and methods to detect rare mutations and copy number variation |
US20140066317A1 (en) * | 2012-09-04 | 2014-03-06 | Guardant Health, Inc. | Systems and methods to detect rare mutations and copy number variation |
SG11201504334WA (en) | 2012-12-10 | 2015-07-30 | Resolution Bioscience Inc | Methods for targeted genomic analysis |
AU2014204058A1 (en) | 2013-01-05 | 2015-07-23 | Foundation Medicine, Inc. | System and method for outcome tracking and analysis |
CA2898326C (en) | 2013-01-18 | 2022-05-17 | Foundation Medicine, Inc. | Methods of treating cholangiocarcinoma |
WO2014152990A1 (en) | 2013-03-14 | 2014-09-25 | University Of Rochester | System and method for detecting population variation from nucleic acid sequencing data |
GB2525568B (en) | 2013-03-15 | 2020-10-14 | Abvitro Llc | Single cell barcoding for antibody discovery |
EP3795696B1 (en) | 2013-03-15 | 2023-04-26 | The Board of Trustees of the Leland Stanford Junior University | Identification and use of circulating nucleic acid tumor markers |
EP2971139A4 (en) | 2013-03-15 | 2016-12-07 | Abbott Molecular Inc | SYSTEMS AND METHOD FOR PROVIDING THE CHANGE OF A GENOMIC COPY NUMBER |
EP2971097B1 (en) | 2013-03-15 | 2018-08-01 | Verinata Health, Inc | Generating cell-free dna libraries directly from blood |
ES2877088T3 (es) | 2013-03-15 | 2021-11-16 | Guardant Health Inc | Procedimiento para detectar cáncer |
SG11201507739TA (en) | 2013-03-19 | 2015-10-29 | Toppan Printing Co Ltd | Method for predicting sensitivity to egfr inhibitor |
EP3470530B1 (en) | 2013-05-23 | 2020-11-25 | The Board of Trustees of the Leland Stanford Junior University | Transposition into native chromatin for analysis ofchromatin |
JP2015096049A (ja) | 2013-11-15 | 2015-05-21 | 凸版印刷株式会社 | Vegf阻害剤長期奏功性予測方法 |
EP4306659A3 (en) | 2014-04-14 | 2024-03-27 | Yissum Research and Development Company of the Hebrew University of Jerusalem Ltd. | A method and kit for determining the tissue or cell origin of dna |
WO2015175705A1 (en) | 2014-05-13 | 2015-11-19 | Board Of Regents, The University Of Texas System | Gene mutations and copy number alterations of egfr, kras and met |
DK3178941T3 (da) | 2014-07-25 | 2022-01-17 | Bgi Genomics Co Ltd | Fremgangsmåde til bestemmelse af fraktionen af cellefrie føtale nukleinsyrer i en prøve af perifert blod fra en gravid kvinde og anvendelse deraf |
WO2016015058A2 (en) | 2014-07-25 | 2016-01-28 | University Of Washington | Methods of determining tissues and/or cell types giving rise to cell-free dna, and methods of identifying a disease or disorder using same |
US20160053301A1 (en) | 2014-08-22 | 2016-02-25 | Clearfork Bioscience, Inc. | Methods for quantitative genetic analysis of cell free dna |
EP3192047A4 (en) | 2014-09-10 | 2018-03-28 | Pathway Genomics Corporation | Health and wellness management methods and systems useful for the practice thereof |
EP3191628B1 (en) | 2014-09-12 | 2022-05-25 | The Board of Trustees of the Leland Stanford Junior University | Identification and use of circulating nucleic acids |
AU2015360298B2 (en) | 2014-12-12 | 2018-06-07 | Verinata Health, Inc. | Using cell-free DNA fragment size to determine copy number variations |
CN113930507A (zh) | 2014-12-31 | 2022-01-14 | 夸登特健康公司 | 疾病的检测和治疗以及用于传送测试结果的系统和方法 |
US10364467B2 (en) | 2015-01-13 | 2019-07-30 | The Chinese University Of Hong Kong | Using size and number aberrations in plasma DNA for detecting cancer |
WO2016127944A1 (en) | 2015-02-10 | 2016-08-18 | The Chinese University Of Hong Kong | Detecting mutations for cancer screening and fetal analysis |
US10844428B2 (en) | 2015-04-28 | 2020-11-24 | Illumina, Inc. | Error suppression in sequenced DNA fragments using redundant reads with unique molecular indices (UMIS) |
PL3387152T3 (pl) | 2015-12-08 | 2022-05-09 | Twinstrand Biosciences, Inc. | Ulepszone adaptory, sposoby i kompozycje do sekwencjonowania dupleksowego |
EP3443066A4 (en) | 2016-04-14 | 2019-12-11 | Guardant Health, Inc. | EARLY DETECTION METHODS FOR CANCER |
-
2016
- 2016-12-16 JP JP2018531350A patent/JP2019507585A/ja not_active Withdrawn
- 2016-12-16 WO PCT/US2016/067356 patent/WO2017106768A1/en active Application Filing
- 2016-12-16 CN CN202311131468.XA patent/CN117174167A/zh active Pending
- 2016-12-16 SG SG11201805119QA patent/SG11201805119QA/en unknown
- 2016-12-16 EP EP16876854.7A patent/EP3390668A4/en active Pending
- 2016-12-16 CA CA3008651A patent/CA3008651A1/en active Pending
- 2016-12-16 CN CN201680081723.6A patent/CN108603228B/zh active Active
-
2017
- 2017-02-27 US US15/442,993 patent/US20170240973A1/en not_active Abandoned
-
2020
- 2020-01-08 US US16/737,819 patent/US11242569B2/en active Active
-
2021
- 2021-04-05 JP JP2021064099A patent/JP2021101732A/ja active Pending
- 2021-12-16 US US17/552,728 patent/US20220356527A1/en active Pending
-
2023
- 2023-06-30 JP JP2023108348A patent/JP2023126874A/ja active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103562407A (zh) * | 2011-02-09 | 2014-02-05 | 伯乐生命医学产品有限公司 | 核酸的分析 |
CN103374518A (zh) * | 2012-04-12 | 2013-10-30 | 维里纳塔健康公司 | 拷贝数变异的检测和分类 |
WO2013166517A1 (en) * | 2012-05-04 | 2013-11-07 | Complete Genomics, Inc. | Methods for determining absolute genome-wide copy number variations of complex tumors |
WO2014039556A1 (en) * | 2012-09-04 | 2014-03-13 | Guardant Health, Inc. | Systems and methods to detect rare mutations and copy number variation |
CN104781421A (zh) * | 2012-09-04 | 2015-07-15 | 夸登特健康公司 | 检测稀有突变和拷贝数变异的系统和方法 |
WO2014165596A1 (en) * | 2013-04-03 | 2014-10-09 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
WO2015100427A1 (en) * | 2013-12-28 | 2015-07-02 | Guardant Health, Inc. | Methods and systems for detecting genetic variants |
WO2015184404A1 (en) * | 2014-05-30 | 2015-12-03 | Verinata Health, Inc. | Detecting fetal sub-chromosomal aneuploidies and copy number variations |
Non-Patent Citations (1)
Title |
---|
"乳腺癌患者血清游离DNA的定量检测";薛建祥等;《肿瘤》;20111231;第31卷(第12期);第1099-1102页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117174167A (zh) | 2023-12-05 |
EP3390668A1 (en) | 2018-10-24 |
WO2017106768A1 (en) | 2017-06-22 |
JP2021101732A (ja) | 2021-07-15 |
US20170240973A1 (en) | 2017-08-24 |
CN108603228A (zh) | 2018-09-28 |
CA3008651A1 (en) | 2017-06-22 |
JP2023126874A (ja) | 2023-09-12 |
US20220356527A1 (en) | 2022-11-10 |
EP3390668A4 (en) | 2020-04-01 |
SG11201805119QA (en) | 2018-07-30 |
US11242569B2 (en) | 2022-02-08 |
JP2019507585A (ja) | 2019-03-22 |
US20200140960A1 (en) | 2020-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108603228B (zh) | 通过分析无细胞dna确定肿瘤基因拷贝数的方法 | |
US20210257048A1 (en) | Methods and systems for calling mutations | |
TWI832482B (zh) | 核酸鹼基修飾的測定 | |
US20180230530A1 (en) | Methods and systems for detecting genetic variants | |
AU2017336153B2 (en) | Methods for multi-resolution analysis of cell-free nucleic acids | |
CA3111019A1 (en) | Methods and systems for high-depth sequencing of methylated nucleic acid | |
US20200327954A1 (en) | Methods and systems for differentiating somatic and germline variants | |
CN110914450A (zh) | 无细胞dna的体细胞来源或种系来源的鉴定 | |
US20240076720A1 (en) | Methods for analyzing nucleic acids | |
US20230360727A1 (en) | Computational modeling of loss of function based on allelic frequency | |
JP2024056939A (ja) | 生体試料のフィンガープリンティングのための方法 | |
JP2024056984A (ja) | エピジェネティック区画アッセイを較正するための方法、組成物およびシステム | |
EP3894600B1 (en) | Method for quantifying gene fusion dna | |
US20200075124A1 (en) | Methods and systems for detecting allelic imbalance in cell-free nucleic acid samples | |
US20220068433A1 (en) | Computational detection of copy number variation at a locus in the absence of direct measurement of the locus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |