CN113593637B - 测序方法及其分析方法和系统、计算机可读存储介质和电子设备 - Google Patents

测序方法及其分析方法和系统、计算机可读存储介质和电子设备 Download PDF

Info

Publication number
CN113593637B
CN113593637B CN202110480833.2A CN202110480833A CN113593637B CN 113593637 B CN113593637 B CN 113593637B CN 202110480833 A CN202110480833 A CN 202110480833A CN 113593637 B CN113593637 B CN 113593637B
Authority
CN
China
Prior art keywords
reads
sequencing
read
quality
base
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110480833.2A
Other languages
English (en)
Other versions
CN113593637A (zh
Inventor
樊济才
金欢
陈方
孙雷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Genemind Biosciences Co Ltd
Original Assignee
Genemind Biosciences Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Genemind Biosciences Co Ltd filed Critical Genemind Biosciences Co Ltd
Publication of CN113593637A publication Critical patent/CN113593637A/zh
Application granted granted Critical
Publication of CN113593637B publication Critical patent/CN113593637B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • G16B30/20—Sequence assembly
    • 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/30—Detection of binding sites or motifs
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C12N15/09—Recombinant DNA-technology
    • C12N15/10—Processes for the isolation, preparation or purification of DNA or RNA
    • C12N15/1034—Isolating an individual clone by screening libraries
    • C12N15/1068—Template (nucleic acid) mediated chemical library synthesis, e.g. chemical and enzymatical DNA-templated organic molecule synthesis, libraries prepared by non ribosomal polypeptide synthesis [NRPS], DNA/RNA-polymerase mediated polypeptide synthesis
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869—Methods for sequencing
    • C—CHEMISTRY; METALLURGY
    • 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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00—Computing arrangements based on specific mathematical models
    • G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
    • 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
    • G16B35/00—ICT specially adapted for in silico combinatorial libraries of nucleic acids, proteins or peptides
    • G16B35/10—Design of libraries
    • 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
    • G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • G16B40/20—Supervised data analysis

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Organic Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Genetics & Genomics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Analytical Chemistry (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Biochemistry (AREA)
  • Evolutionary Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Microbiology (AREA)
  • Immunology (AREA)
  • Data Mining & Analysis (AREA)
  • Biomedical Technology (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Library & Information Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Plant Pathology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Mathematical Physics (AREA)

Abstract

本发明的一个目的在于提出一种有效的测序方法。根据本发明的实施例，该方法包括：(1)对芯片表面的测序模板进行第一测序，以便通过形成第一新生测序链获得第一测序数据，所述测序模板通过测序接头连接在所述芯片表面上；(2)对至少一部分所述第一新生测序链的3’‑末端进行第一封闭处理；和(3)对所述测序模板进行第二测序，以便通过形成第二新生测序链获得第二测序数据。

Description

测序方法及其分析方法和系统、计算机可读存储介质和电子 设备
技术领域
本发明涉及测序领域，具体的，本发明涉及测序方法、测序结果分析方法、测序结果分析系统、计算机可读存储介质和电子设备。
背景技术
据报道，早在上世纪八十年代人们就提出了单分子测序。2003年斯坦福大学生物工程系的教授Stephen Quake博士成功演示了第一个单分子DNA测序实验。2008年Helicos公司的第一台单分子测序仪(HeliScope)上市。2009年Korlach与Turner在《科学》杂志上发表文献介绍了PacBio单分子测序技术原理。随后，2010年PacBio公司推出了PacBio RS测序系统，并于2011年正式商用。2014年Oxford Nanopore公司在AGBT(基因组生物学技术进展年会)上展示了其MinION测序系统。据报道，无论是Helicos、PacBio还是MinION测序平台，其单重测序(Single-Pass)的测序错误率均很高，最高可达30％。许多研究显示，上述测序平台的错误类型主要是InDel，且是随机发生的，可以通过重复读取的方法降低其测序错误率。
有文献报道PacBio可以采用CCS(环形一致序列)克服其SMRT测序技术的高错误率问题。另外，MinION通过2D和1D2的测序方法可以将测序准确率大幅提升，最高可达97％的准确率。
有文献报道，Helicos通过双重测序法(Two-Pass)进行测序可以将其测序中的缺失类型的错误率降低至1％以内，但是其所使用的文库是5内端为特异性接头，3为端为polyA接头，其操作过程比较繁琐复杂；同时其变性洗脱DNA链的步骤使用的是热水，可能会洗脱不完全，从而干扰后续的测序引物杂交以及测序；并且芯片表面的DNA链未经处理也是随后的测序引入错误的来源之一。
由此，现有的测序方法有待进一步改进。
发明内容
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。为此，本发明的一个目的在于提出一种有效的测序方法。
在本发明的第一方面，本发明提出了一种测序方法，根据本发明的实施例，该方法包括：(1)对芯片表面的测序模板进行第一测序，以便通过形成第一新生测序链获得第一测序数据，所述测序模板通过测序接头连接在所述芯片表面上；(2)对至少一部分所述第一新生测序链的3’-末端进行第一封闭处理；和(3)对所述测序模板进行第二测序，以便通过形成第二新生测序链获得第二测序数据。
根据本发明的实施例，通过进行两轮测序，后续可以互相校正，以提高测序结果的准确率，同时，在进行第一轮测序，即第一测序之后，通过对残留在芯片表面的新生测序链的3’-末端进行封闭处理，能够有效地避免在第二轮测序即第二测序过程中产生干扰信号。由此，能够进一步提高测序结果的准确性。
在本发明的第二方面，本发明提出了一种测序结果分析方法，根据本发明的实施例，所述测序结果包括第一测序数据和第二测序数据，所述第一测序数据和所述第二测序数据均由多个读段构成，所述第一测序数据中的至少一部分所述读段在所述第二测序数据中存在对应读段，所述第一测序数据和所述第二测序数据是通过前面所述的方法得到的，所述测序结果分析方法包括：(a)基于所述第一测序数据和所述第二测序数据各自的至少一部分进行相互校正，以便获得最终序列信息。
根据本发明的实施例，通过对两轮测序的结果进行相互校正能够提高测序结果的准确率。另外，如前所述，在进行第一轮测序，即第一测序之后，通过对残留在芯片表面的新生测序链的3’-末端进行封闭处理，能够有效地避免在第二轮测序即第二测序过程中产生干扰信号。由此，能够进一步提高测序结果的准确性。
在本发明的第三方面，本发明还提出了一种测序结果分析系统。根据本发明的实施例，该系统包括：测序设备，所述测序设备适于通过前面所述的方法获得测序结果，所述测序结果包括第一测序数据和第二测序数据，所述第一测序数据和所述第二测序数据均由多个读段构成，所述第一测序数据中的至少一部分所述读段在所述第二测序数据中存在对应读段；分析设备，所述分析设备适于基于所述第一测序数据和所述第二测序数据各自的至少一部分进行相互校正，以便获得最终序列信息。
利用该系统能够有效地实施前面所述的测序结果分析方法，从而通过对两轮测序的结果进行相互校正能够提高测序结果的准确率。另外，如前所述，在进行第一轮测序，即第一测序之后，通过对残留在芯片表面的新生测序链的3’-末端进行封闭处理，能够有效地避免在第二轮测序即第二测序过程中产生干扰信号。由此，能够进一步提高测序结果的准确性。
另外，本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，根据本发明的实施例，该程序被处理器执行时实现前面所述方法的步骤。
本发明还提供了一种电子设备，其包括：前面所述的计算机可读存储介质；以及一个或者多个处理器，用于执行所述计算机可读存储介质中的程序。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是根据本发明一个实施例的测序方法的流程示意图；
图2是根据本发明又一个实施例的测序方法的流程示意图；
图3是根据本发明又一个实施例的测序方法的流程示意图；
图4是根据本发明一个实施例的测序结果分析方法的流程示意图；
图5是根据本发明又一个实施例的测序结果分析方法的流程示意图；
图6是根据本发明又一个实施例的测序结果分析方法的流程示意图；
图7是根据本发明一个实施例的测序结果分析系统的结构示意图；
图8是根据本发明又一个实施例的测序结果分析系统的结构示意图；
图9是根据本发明又一个实施例的测序结果分析系统的结构示意图；
图10是根据本发明一个实施例的获得Reads1和Reads2的测序方法示意图；
图11是根据本发明一个实施例的测序文库构建示意图；
图12是根据本发明一个实施例的获得Consensus Reads(一致序列/共同序列)的分析方法流程示意图。
具体实施方式
本发明中的术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
在本发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系，除非另有明确的限定。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
在本发明的第一方面，本发明提出了一种降低单分子测序平台(例如GenoCareTM单分子测序平台)测序输出序列噪声及错误率的测序与分析方法测序方法，参考图1～3和图10～12，对根据本发明实施例的测序方法进行描述。
根据本发明的实施例，该方法包括：
S10第一测序获得第一测序数据
在该步骤中，对芯片表面的测序模板进行第一测序，以便通过形成第一新生测序链获得第一测序数据，所述测序模板通过测序接头连接在所述芯片表面上。
需要说明的是，在本文中所使用的术语“芯片”是指在测序平台中所使用的测序芯片，只要是通过边合成边测序的原理进行测序，就可以采用本发明的方法进行处理，其中，优选采用单分子测序平台，例如GenoCare单分子测序平台。当然，本领域技术人员能够理解的是，还可以适用于其他的单分子测序平台，在此不再赘述。
参考图2，在步骤S10之前，还可以通过下列步骤获得可以用于单分子测序平台的芯片：
S10a：使测序文库中的文库分子与芯片表面的测序接头进行杂交；
S10b：利用所述文库分子作为初始模板，通过合成互补链形成所述测序模板；和
S10c除去所述初始模板，并对所述芯片表面的核酸分子的3’-末端进行第二封闭处理。
由此，通过第二封闭处理，可以进一步去除残留的活性3’-末端对后续反应的影响。
参考图3，在进行S10c步骤之前，还可以进一步包括S11b：对步骤S10b中延伸不完全的所述互补链的3’-末端进行第三封闭处理。由此，可以进一步提高测序的准确性，降低不期望的测序噪音。
S20对至少一部分所述第一新生测序链的3’-末端进行第一封闭处理在该步骤中，对至少一部分所述第一新生测序链的3’-末端进行第一封闭处理，通过第一封闭处理，可以有效增加有效数据量，降低无效数据对信息分析的干扰。
根据本发明的一个实施例，步骤S20包括去除所述芯片表面的所述第一新生测序链，并对残余在所述芯片表面的所述第一新生测序链的3’-末端进行第一封闭处理。
根据本发明的一个实施例，步骤S20包括对所述第一新生测序链的3’-末端进行第一封闭处理，并去除封闭后的第一新生测序链。
S30第二测序获得第二测序数据
在该步骤中，对所述测序模板进行第二测序，以便通过形成第二新生测序链获得第二测序数据。
根据本发明的实施例，通过进行两轮测序，后续可以互相校正，以提高测序结果的准确率，同时，在进行第一轮测序，即第一测序之后，通过对残留在芯片表面的新生测序链的3’-末端进行封闭处理，能够有效地避免在第二轮测序即第二测序过程中产生干扰信号。由此，能够进一步提高测序结果的准确性。
根据本发明的实施例，所述第一封闭处理、所述第二封闭处理和所述第三封闭处理可以分别独立地通过使3据本末端羟基与延伸反应阻断剂相连而进行的。由此，可以进一步提高封闭效果，从而进一步提高测序的准确性，降低不期望的测序噪音。
根据本发明的实施例，所述延伸反应阻断剂为ddNTP或其衍生物。由此，可以进一步提高封闭效果，从而进一步提高测序的准确性，降低不期望的测序噪音。
根据本发明的实施例，所述第一封闭处理、所述第二封闭处理和所述第三封闭处理分别独立地采用DNA聚合酶和末端转移酶的至少之一进行。由此，可以进一步提高封闭效果，从而进一步提高测序的准确性，降低不期望的测序噪音。
根据本发明的实施例，所述第一封闭处理和所述第三封闭处理分别独立地通过聚合酶连接所述ddNTP或其衍生物，所述第二封闭处理通过所述末端转移酶连接所述ddNTP或其衍生物。由此，可以进一步提高封闭效果，从而进一步提高测序的准确性，降低不期望的测序噪音。
在本发明的第二方面，本发明提出了一种测序结果分析方法，利用该方法能够对前面所述的测序方法所产生的两轮测序数据进行有效分析，从而进一步提高测序的准确性，避免测序误差。
参考图4～图6、以及图12，根据本发明的实施例，所述测序结果包括第一测序数据和第二测序数据，所述第一测序数据和所述第二测序数据均由多个读段构成，所述第一测序数据中的至少一部分所述读段在所述第二测序数据中存在对应读段，所述第一测序数据和所述第二测序数据是通过前面所述的方法得到的，所述测序结果分析方法包括：
基于所述第一测序数据和所述第二测序数据各自的至少一部分进行相互校正，以便获得最终序列信息。
根据本发明的实施例，所述相互校正包括下列步骤：
在所述第一测序数据和所述第二测序数据中选择高质量读段和所述高质量读段的对应读段，所述读段的长度不低于预定长度，所述读段具有不低于预定质量阈值的测序质量；和
将所述高质量读段与所述高质量读段的对应读段进行比对，并基于所述比对结果进行序列信息校正。
根据本发明的实施例，通过对两轮测序的结果进行相互校正能够提高测序结果的准确率。另外，如前所述，在进行第一轮测序，即第一测序之后，通过对残留在芯片表面的新生测序链的3’-末端进行封闭处理，能够有效地避免在第二轮测序即第二测序过程中产生干扰信号。由此，能够进一步提高测序结果的准确性。
参考图5，所述相互校正包括下列步骤：
S100构建第一读段集合
在该步骤中，根据所述读段的长度，基于所述第一测序数据，构建第一读段集合，所述第一读段集合中的每一个读段长度均不低于第一预定长度。
S200构建第二读段集合和第三读段集合
在该步骤中，根据所述对应读段的长度，基于所述第一读段集合，构建第二读段集合和第三读段集合，所述第二读段集合中每一个读段的所述对应读段的长度均不低于第二预定长度，所述第三读段集合中每一个读段的所述对应读段的长度均处于预定长度范围内。
S300构建第四读段集合和第五读段集合
在该步骤中，根据所述第二读段集合中所述读段及其所述对应读段的测序质量，基于所述第二读段集合及其所述对应读段，构建第四读段集合和第五读段集合。
根据本发明的实施例，所述第四读段集合和所述第五读段集合分别是按照下列原则确定的：
将所述第二读段集合中的所述读段与其所述对应读段进行测序质量比较；
选择测序质量高的一方作为所述第四读段集合的元素，选择测序质量低的一方作为所述第五读段集合的元素；
对于测序质量相同的情形，则选择来自所述第二读段集合的所述读段作为所述第四读段集合的元素，则选择所述对应读段作为所述第五读段集合的元素。
S400构建第六读段集合
在该步骤中，利用测序质量，对所述第四读段集合进行过滤处理，以便构建第六读段集合，所述第六读段集合中的所述读段的测序质量均不低于第一预定质量阈值。
S500构建第七读段集合
在该步骤中，利用所述第六读段集合，从所述第五读段集合中选择与所述第六读段集合中的所述读段对应的所述读段，以便构建第七读段集合。
S600将第六读段集合和第七读段集合进行比较，确定第一差异位点
在该步骤中，将所述第六读段集合与所述第七读段集合进行读段比对，并在所述第六读段集合的所述读段上确定第一差异位点
S700对第一差异位点进行校正
在该步骤中，利用预先确定的测序误差预测模型，对所述第一差异位点进行校正，以便确定第一序列信息，所述测序误差预测模型用于确定在测序过程中，差异位点发生插入或者缺失的概率。
参考图6，在得到第一测序信息后，还可以进一步包括：
S400a构建第八读段集合
在该步骤中，利用测序质量，对所述第三读段集合进行过滤处理，以便构建第八读段集合，其中，所述第八读段集合中的所述读段的测序质量均不低于第二预定质量阈值。
S500a构建第九读段集合
在该步骤中，利用所述第八读段集合，从所述第二测序数据中选择与所述第七读段集合中的所述读段对应的所述读段，以便构建第九读段集合。
S600a将第八读段集合和第九读段集合进行比对，确定第二差异位点
在该步骤中，将所述第八读段集合与所述第九读段集合进行读段比对，并在所述第八读段集合的所述读段上确定第二差异位点。
S700a对第二差异位点进行校正确定第二序列信息。
在该步骤中，利用所述测序误差预测模型，对所述第二差异位点进行校正，以便确定第二序列信息。
根据本发明的实施例，所述测序误差预测模型是基于所述第一测序数据和所述第二测序数据与参考基因组的比对结果，对朴素贝叶斯模型进行训练获得的。
根据本发明的实施例，，针对所述第一差异位点和所述第二差异位点：
如果来自所述第六读段集合的读段在所述差异位点存在碱基，来自所述第七读段集合的对应读段在所述差异位点不存在碱基，并且在所述差异位点发生缺失的概率为50％以上，则保留所述第六读段集合的读段在所述差异位点的碱基作为最终测序结果。
如果来自所述第六读段集合的读段在所述差异位点不存在碱基，来自所述第七读段的读段集合的定读段在所述差异位点存在碱基，并且在所述差异位点发生插入的概率为50％以上，则保留所述第六读段集合的读段在所述差异位点的碱基作为最终测序结果。
如果来自所述第六读段集合的读段在所述差异位点存在碱基，来自所述第七读段的读段集合的定读段在所述差异位点也存在碱基，则选择所述第六读段集合的读段在所述差异位点的碱基作为最终测序结果。
根据本发明的实施例，所述第一预定长度和所述第二预定长度分别独立地不低于20bp，优选不低于25bp，所述预定长度范围为10～25bp；所述第一预定质量阈值和所述第二预定质量阈值分别独立地不低于50，优选不低于60。
由此，根据本发明的实施例，本发明提供了一种联合使用接头D7-S1-T/D9-S2和测序引物D7S1T-R2P，利用GenoCare单分子测序平台进行Two-Pass测序获得Reads1和Reads2的测序方法。
本发明进一步的目的在于，提供一种对上述Two-Pass测序方法获得的Reads1和Reads2进行分析获得Consensus Reads的分析方法。该分析方法可以显著降低输出的Consensus Reads(一致序列/共同序列)中的噪声序列以及碱基错误率。
根据本发明的实施例，用于GenoCare单分子Two-Pass测序文库构建的接头以及测序引物，所述接头由寡核苷酸链D7-S1-T和带有5带磷酸基团修饰的D9-S2退火获得，所述测序引物为D7S1T-R2P。其中，所述D7-S1-T的序列为SEQ ID NO：1，所述D9-S2的序列为SEQ IDNO：2，所述D7S1T-R2P的序列为SEQ ID NO：3。
根据本发明的实施例，本发明提供一种使用上述接头和测序引物在GenoCare单分子测序平台进行Two-Pass测序获得Reads1和Reads2的方法，包括：
步骤一：构建Two-Pass测序文库，使用文库制备试剂盒(Universal DNALibrary Prep Kit for Illumina V2(ND606-01))将退火后的D7-S1-T/D9-S2接头与准备好的片段化的人类gDNA进行连接，连接后无需进行PCR扩增，直接使用纯化试剂盒(VAHTSDNA Clean Beads(N411-01))进行纯化获得目的文库；
步骤二：将步骤一获得的文库与测序芯片表面接头进行杂交；
步骤三：对步骤二中杂交于芯片表面的初始模板进行互补链合成；
步骤四(可选地)：对步骤三中延伸不完全的新生链的3'末端进行封闭，减少其对测序过程的干扰；
步骤五：变性去除步骤二中杂交于芯片表面的初始模板；
步骤六：对芯片表面残余接头的3'末端进行封闭，减少其对测序过程的干扰；
步骤七：以步骤三中合成的互补链为模板杂交所述测序引物D7S1T-R2P；
步骤八：以步骤三中合成的互补链为模板，以步骤七中杂交的测序引物D7S1T-R2P为引物进行Read1测序；
步骤九：变性去除步骤八中的新生测序链；
步骤十：对经过步骤九处理后可能残余的步骤八中的新生测序链的3'末端：进行封闭，防止其在Read2测序时继续延伸；
步骤十一：以步骤三中合成的互补链为模板杂交所述测序引物D7S1T-R2P；
步骤十二：以步骤三中合成的互补链为模板，以步骤十一中杂交的测序引物D7S1T-R2P为引物进行Read2测序；
步骤十三：对步骤八和步骤十二获得的测序数据进行拆分获得坐标一一对应的Reads1和Reads2两部分序列。
进一步地，本发明提供一种对上述获得的Reads1和Reads2进行分析获得Consensus Reads的分析方法，包括：
步骤十四：校正模型构建，提取步骤十三获得的Reads1和Reads2序列中同一坐标两次测序读长均≥列中同一坐标的Reads，分别输出为T1(Read1)和T2(Read2)两个文件，将T1和T2中的Reads分别与参考基因组进行比对，通过朴素贝叶斯方法计算在不同的前后碱基组合下，中间碱基发生Deletion或Insertion的概率。在预测过程中，对于不同的前后碱基组合下的中间碱基，根据模型中发生Deletion或Insertion的概率决定是否保留中间碱基。若Deletion的概率大于50％，则保留中间碱基，反之舍弃该中间碱基。
步骤十五：对步骤十三获得的Reads1数据按读长进行过滤，将读长≥25bp的Reads1序列集合名为Fa1。使用25bp长度对短读序列进行过滤可以去除一部分噪声序列，提高测序数据mapping的准确性。
步骤十六：根据步骤十三获得的Read2的序列读长将与其对应的步骤十五获得的Fa1拆分为两个集合，其中Read2合，其中的步骤时对应的Fa1中的Reads的集合命名为Fa2，10bp合，其中的步骤十＜25bp时对应的Fa1中的Reads的集合命名为Fa3。这里将Fa1拆分为两部分进行分析的目的是在提高Consensus Reads的准确率的同时，减少因长度过滤而导致的数据通量的损失。
步骤十七：将步骤十六获得的Fa2中的Reads和与其坐标相对应的步骤十三获得的Reads2中的Reads的Q值进行比较，取Q值较高的Reads(Q值相等时取Fa2中的Reads)的集合命名为Fa4，取Q值较低的Reads(Q值相等时取Reads2中的Reads)的集合命名为Fa5。该步骤的目的是将Reads1和Reads2中的序列划分为测序质量相对较高和较低的两个集合，保证最终输出的Consensus Reads中的序列是两次测序中测序质量较高相对更准确的Reads。
步骤十八：对步骤十七获得的Fa4和Fa5进行进一步的过滤，将Fa4中Q值≥a4的Reads的集合命名为Fa6，将与Fa6中Reads坐标一一对应的Fa5中的Reads的集合命名为Fa7。
步骤十九：将步骤十八获得的Fa6和Fa7中的Reads一一对齐，根据其序列的相似性进行分级，并以Fa7为参考序列对Fa6进行校正，在Fa6序列中标记出与Fa7不同的位置，根据步骤十四构建的校正模型逐个判断不同位置碱基是否为Deletion或者Insertion，从而得到校正后的用于输出的Consensus Reads Part1。
该步骤中的不同的位置表述仅针对某一位置上只在Fa6或者Fa7的一条Reads上测出碱基。这时步骤十四构建的矫正模型将判断该碱基是否应该被保留。对于某一位置上Fa6和Fa7上均测出碱基，但是碱基种类不一致，则以Fa6的碱基为准，本模型不对上述情况进行矫正。
步骤二十：对所述步骤十六获得的Fa3中的Reads进行进一步过滤，将所述Fa3中Q值≥a3的Reads的集合命名为Fa8；
步骤二十一：提取所述步骤十三获得的Reads2中与所述Fa8中坐标一一对应的Reads的集合，命名为Fa9；
步骤二十二：将所述Fa8与Fa9中的Reads一一对齐，根据其序列的相似性进行分级，并以Fa9为参考序列对Fa8进行校正，在Fa8中标记出与Fa8不同的位置，根据步骤十四构建的校正模型逐个判断不同位置碱基是否为Deletion(缺失)或者Insertion(插入)，从而得到校正后的用于输出的Consensus Reads Part2；
步骤二十三：根据不同的应用需求将不同的相似性级别的所述Consensus ReadsPart1与所述Consensus Reads Part2进行合并，获得输出的Consensus Reads。
根据本发明的实施例，对于步骤二所述文库与测序芯片表面接头杂交过程包括(试剂常规)：
1)将用于杂交的文库在90-100文预变性2-5分钟；
2)将从所述步骤1)获得的产物迅速在冰水混合物上冷却2分钟以上，获得变性的杂交文库母液；
3)将从所述步骤2)获得的变性的杂交文库母液使用80％的GenoCare杂交液稀释至合适的浓度，优选地为0.1～2nM，得到稀释的杂交文库；
4)将从所述步骤3)获得的30～50μL稀释的杂交文库通入使用复溶试剂预处理好的测序芯片通道，在40～60文杂交10～30分钟；
5)向芯片通道通入200～1000μL体积的清洗液1，去除步骤4)杂交后剩余的稀释的杂交文库；
6)向芯片通道通入200～1000μL体积的清洗液2，去除步骤5)中的清洗液1，完成文库与测序芯片表面接头杂交。
根据本发明的实施例，所述GenoCare杂交液为3*SSC溶液。
根据本发明的实施例，所述复溶试剂的组分包括：清洗液1，组分包括：150mM的氯化钠，15mM的柠檬酸钠，150mM的4-羟乙基哌嗪乙磺酸，0.1％的十二烷基硫酸钠。
根据本发明的实施例，清洗液3，组分包括：450mM的氯化钠，45mM的柠檬酸钠。
根据本发明的实施例，所述清洗液2的组分包括：150mM的氯化钠，150mM的4-羟乙基哌嗪乙磺酸。
根据本发明的实施例，对于步骤三所述初始模板的互补链合成的过程包括：
1)向芯片通道中通入200～1000μL体积的延伸试剂，在50～70试条件下反应5～10分钟；
2)向芯片道中通入200～1000μL体积的清洗液1，去除步骤1)中反应后的延伸试剂；
3)向芯片道中通入200～1000μL体积的清洗液2，去除步骤2)中的清洗液1，完成初始模板的互补链的合成。
根据本发明的实施例，所述延伸试剂的组分包括：10～100U/mL的DNA聚合酶，优选地，可以是Bst DNA聚合酶、Bsu DNA聚合酶、Klenow DNA聚合酶等，0.2～2mM的dNTP，0.5～2M的甜菜碱，20mM的三羟甲基氨基甲烷，10mM的氯化钠，10mM的氯化钾，10mM的硫酸铵，3mM的氯化镁，0.1％的Triton X-100，pH值为8.3。
根据本发明的实施例，对于步骤四所述对延伸不完全的链的3据本发明进行封闭的过程包括：
1)向芯片通道中通入200～1000μL体积的封闭试剂1，在30～60试条件下反应5～30分钟；
2)向芯片通道中通入200～1000μL体积的清洗液1，去除步骤1)中反应后的封闭试剂1，完成对延伸不完全的链的3完成对延的封闭。
根据本发明的实施例，所述封闭试剂1的组分包括：10～100U/mL的DNA聚合酶，优选地，可以是Klenow DNA聚合酶、Bsu DNA聚合酶、N9 DNA聚合酶等，10～100μM的ddNTP，5mM的氯化锰，20mM的三羟甲基氨基甲烷，10mM的氯化钠，10mM的氯化钾，10mM的硫酸铵，3mM的氯化镁，0.1％的Triton X-100，pH值为8.3。
根据本发明的实施例，对于步骤五所述去除初始模板的过程包括：
1)向芯片通道中通入200～1000μL体积的变性试剂，优选地，所述变性试剂可以是甲酰胺、0.1M的NaOH等，于50～60反应2～5分钟；
2)向芯片通道中通入200～1000μL体积的清洗液1，去除步骤1)中反应后的变性试剂和从芯片上变性分离的初始模板；
重复步骤1)和步骤2)一次，完成初始模板的去除。
根据本发明的实施例，对于步骤六所述封闭芯片表面残余接头3据本发明的过程包括：
1)向芯片通道中通入200～1000μL体积的清洗液2；
2)向芯片通道中通入200～1000μL体积的封闭试剂2，在30～60试条件下反应5～30分钟；
3)向芯片通道中通入200～1000μL体积的清洗液1，去除步骤2)中反应后的封闭试剂2，完成对芯片表面残余接头3完成对芯的封闭；
根据本发明的实施例，所述封闭试剂2的组分包括：100U/mL的末端转移酶(Terminal Transferase(NEB，M0315L))，1×Terminal Transferase Buffer，0.25mM的氯化钴，10～100μμM的ddNTP。
根据本发明的实施例，对于步骤七所述杂交测序引物D7S1T-R2P的过程包括：
1)将所述测序引物D7S1T-R2的母液使用清洗液3稀释至合适的浓度，优选地为0.1～1μM，获得稀释的测序引物杂交液；
2)将从步骤1)获得的200～1000μL稀释的测序引物杂交液通入芯片通道，在50～60引条件下杂交10～30分钟；
3)向芯片通道中通入200～1000μL体积的清洗液1，去除步骤2)杂交后残余的测序引物；
4)向芯片通道中通入200～1000μL体积的清洗液2，去除步骤3)中的清洗液1，完成测序引物的杂交。
根据本发明的实施例，对于步骤八所述Read1的测序过程参照GenoCare单分子双色测序通用试剂盒(备案号：粤深械备20190887号)中的描述进行。
根据本发明的实施例，对于步骤九所述去除新生测序链的过程参照步骤五进行。
根据本发明的实施例，对于步骤十所述残余新生链的3'末端的封闭过程参照步骤四进行。
根据本发明的实施例，对于步骤十一所述杂交测序引物D7S1T-R2P的过程参照步骤七进行。
根据本发明的实施例，对于步骤十二所述Read2的测序过程参照步骤八进行。
根据本发明的实施例，对于步骤十三所述将测序数据进行拆分获得坐标一一对应的Reads1和Reads2两部分序列的过程包括：
使用python语言按测序循环数将BaseCalling输出的“.fa_”文件中的每条Read从中间平均分为两个部分，分别输出为两份序列坐标一致的“.fa_”文件“Reads1.fa_”、“Reads2.fa_”；
使用python语言将从步骤1)获得的“Reads1.fa_”、“Reads2.fa_”文件中所用Reads中的字符“_”移除，输出“Reads1.fa”、“Reads2.fa”文件，完成将测序数据进行拆分获得坐标一一对应的Reads1和Reads2两部分序列。
根据本发明的实施例，对于步骤十四所述校正模型的构建过程包括：
提取步骤十三获得的Reads1和Reads2序列中同一坐标两次测序读长均≥25bpReads，分别输出为T1(Read1)和T2(Read2)两个fast文件；
将从步骤1)获得的T1和T2文件中对应的两条reads进行滑动对齐，标记对齐结果中两条reads相同和不同的Base，得到Common Reads；
将从步骤1)获得的T1和T2文件分别和参考序列做mapping得到Sam1和Sam2文件；
根据从步骤3)获得的Sam1和Sam2中对应的且mapping到同一位置的reads，找到参考序列中的最长公共子串Ref Reads；
将从步骤2)获得的Common Reads中两次测序不同的base与从步骤4)获得的RefReads进行比较，使用朴素贝叶斯方法计算在不同的前后碱基组合下中间碱基发生Deletion或Insertion的概率，完成校正模型的构建。
与现有技术相比，本发明有益效果在于：本发明结合GenoCare单分子测序平台的特点，提供了一套联合使用接头D7-S1-T/D9-S2和测序引物D7S1T-R2P，利用GenoCare单分子测序平台进行Two-Pass测序获得Reads1和Reads2的测序方法。另一方面，本发明提供了一种使用所述Two-Pass测序方法获得的Reads1和Reads2进行分析获得Consensus Reads的分析方法。该分析方法可以显著降低输出的Consensus Reads(一致序列/共同序列)中的噪声序列以及碱基错误率。
在本发明的第三方面，本发明还提出了一种能够实施上述测序结果分析方法的测序结果分析系统。参考图7～9，根据本发明的实施例，该系统包括：测序设备，所述测序设备适于通过前面所述的方法获得测序结果，所述测序结果包括第一测序数据和第二测序数据，所述第一测序数据和所述第二测序数据均由多个读段构成，所述第一测序数据中的至少一部分所述读段在所述第二测序数据中存在对应读段；分析设备，所述分析设备适于基于所述第一测序数据和所述第二测序数据各自的至少一部分进行相互校正，以便获得最终序列信息。
利用该系统能够有效地实施前面所述的测序结果分析方法，从而通过对两轮测序的结果进行相互校正能够提高测序结果的准确率。另外，如前所述，在进行第一轮测序，即第一测序之后，通过对残留在芯片表面的新生测序链的3’-末端进行封闭处理，能够有效地避免在第二轮测序即第二测序过程中产生干扰信号。由此，能够进一步提高测序结果的准确性。
另外，本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现前面所述方法的步骤。
本发明还提供了一种电子设备，其包括：前面所述的计算机可读存储介质；以及一个或者多个处理器，用于执行所述计算机可读存储介质中的程序。
下面参考具体实施例，对本发明进行描述，需要说明的是，这些实施例仅仅是描述性的，而不以任何方式限制本发明。
实施例
本实施例提供了一套降低GenoCare单分子测序平台测序输出序列噪声及错误率的测序与分析方法。其中Genocare单分子测序平台是使用TIRF成像系统检测掺入核苷酸种类的平台。Genocare测序过程有多种方式，第一种方式：四种核苷酸带有同种荧光信号，每轮反应加入一种核苷酸进行信号检测；第二种方式：四种核苷酸带有两种不同的荧光信号，每轮反应加入两种核苷酸进行信号检测；第三种方式：四种核苷酸带有四种不同的荧光信号，每轮反应加入四种核苷酸进行信号检测。具体测序过程可参看文章Single molecμLetargeted sequencing for cancer gene mutation detection，Scientific RepoRts|6:26110|DOI:10.1038/srep26110、专利CN201680047468.3、CN201910907555.7、CN201880077576.4或CN201911331502.1中测序过程的描述。
进一步地，本实施例提供的所述测序与分析方法包括：
一套联合联合使用接头D7-S1-T/D9-S2和测序引物D7S1T-R2P，利用GenoCare单分子测序平台进行Two-Pass测序获得Reads1和Reads2的测序方法。其中，所述接头D7-S1-T/D9-S2由寡核苷酸链D7-S1-T和带有5带磷酸基团修饰的D9-S2组成。所述D7-S1-T的序列为SEQ ID NO：1，所述D9-S2的序列为SEQ ID NO：2，所述测序引物D7S1T-R2P的序列为SEQ IDNO：3。具体地，本发明中所涉及的引物序列和名称总结如表1所示。
表1
一套对上述Two-Pass测序方法获得的Reads1和Reads2进行分析获得ConsensusReads的分析方法。该分析方法可以显著降低输出的Consensus Reads中的噪声序列以及碱基错误率。
进一步地，本实施例提供的所述一套联合联合使用接头D7-S1-T/D9-S2和测序引物D7S1T-R2P，利用GenoCare单分子测序平台进行Two-Pass测序获得Reads1和Reads2的测序方法包括：
步骤一：构建Two-Pass测序文库。使用Universal DNA Library Prep Kitfor Illumina V2(ND606-01)将退火后的D7-S1-T/D9-S2接头与准备好的片段化的人类gDNA进行连接，连接后无需进行PCR扩增，直接使用VAHTS DNA Clean Beads(N411-01)进行纯化获得目的文库。
具体地，本实施例中构建Two-Pass测序文库的步骤包括：
1)人类gDNA片段化：使用Covaris设置参数Peak Power，75；Duty Factor，25；Cycles/Burst，50；Time(s)，250对0.1～1ug人类gDNA进行超声打断获得100～300bp的DNA片段。可选地，本步骤也可以使用酶切的方法实现。
2)DNA片段进行末端修复和加A尾，反应体系如表2所示。
表2
H2O (16.2-X)μL
End Prep Mix 3.8μL
DNA片段(总量50ng) XμL
Total 20μL
反应条件为：20条反应15分钟，接着在65，条件下反应10分钟。
3)末端修复加A产物与接头进行连接，反应体系如表3所示。
表3
末端修复加A产物 20μL
D7-S1-T/D9-S2接头(20μM) 5μL
Ligation Mix 25μL
Total 50μL
反应条件为，混匀后室温放置15min。
4)连接产物纯化
纯化使用VAHTS DNA Clean Beads(N411-01)的试剂和说明书所示步骤进行纯化，对该试剂盒的操作说明书稍作改动，以回收产物10μL，完成测序文库的构建。具体步骤如下：
a)将连接后的PCR体系转移至1.5mL EP管中，加入0.8加(4μL)磁珠，吹打混匀10次，室温放置3分钟；
b)将1.5mL EP管放置在磁力架上，静置2-3分钟，移去上清；
c)用200μ体积80％乙醇洗涤，漂洗磁珠，室温孵育30sec，小心移除上清；
d)开盖干燥磁珠约5-10分钟至残余乙醇完全挥发；
e)加入22μL体积的去离子水从磁力架上去取进行洗脱，充分混匀后室温静置3分钟，置于磁力架上3分钟，待液体澄清后，回收产物20μL，待，再加入1.2×(24μL)磁珠，吹打混匀10次，室温放置3分钟；
f)将1.5mL EP管放置在磁力架上，静置2-3分钟，移去上清；
g)重复步骤c)-d)一次；
h)加入11μL体积的去离子水从磁力架上取下进行洗脱，充分混匀后室温静置3分钟，置于磁力架上3分钟，待液体澄清后，回收产物10μL，完成测序文库构建。
5)定量及检测
使用Qubit 3.0仪器和Qubit dsDNA HS检测试剂盒对所述构建的文库进行浓度检测。
使用Labchip DNA HS检测试剂盒和LabChip仪器对所述构建的文库进行片段分布检测。
步骤二：将步骤一获得的文库与测序芯片表面探针进行杂交。
芯片选择：
1)芯片选择：所用的芯片为环氧基修饰的芯片，通过探针上的氨基和芯片表面的环氧基团反应的方法(专利申请号CN201811191589.2)来固定探针，序列为5’-TTTTTTTTTTTCCTTGATACCTGCGACCATCCAGTTCCACTCAGATGTGTATAAGAGACAGT-3’(SEQ ID NO:4)。
文库与芯片上探针杂交过程如下：
1)取3μL体积20nM浓度的步骤一构建的测序文库，加入3μL的去离子水，混合均匀，于95℃变性5分钟；
2)将从步骤1)获得的变性文库迅速置于冰水混合物冷却2分钟以上；
3)向步骤2)的产物中加入24μL体积的GenoCare杂交液，将文库稀释至2nM的工作浓度。其中杂交液3xSSC缓冲液，3xSSC溶液是将20*SSC缓冲液((西格玛，#S6639-1L))用无核酸酶水(Rnase-free水)稀释而成。
4)将从步骤3)获得的30μL体积稀释的杂交文库通入从芯片的一条通道中，于42℃杂交反应30分钟，然后冷却至室温；
5)向步骤4)获得的杂交通道中通入200μL体积的清洗液1，去除未杂交至芯片表面的文库；
向芯片杂交通道通入200μL体积的清洗液2，替换通道内的清洗液1，完成文库与测序芯片表面接头的杂交。
清洗液1组分包括：150mM的氯化钠，15mM的柠檬酸钠，150mM的4-羟乙基哌嗪乙磺酸，0.1％的十二烷基硫酸钠。
清洗液2的组分包括：150mM的氯化钠，150mM的4-羟乙基哌嗪乙磺酸。
步骤三：初始模板进行互补链合成。
初始模板为步骤二中与探针进行杂家的文库，初始模板互补链合成的具体步骤如下：
1)将步骤二中完成文库杂交的芯片置于GenoCare测序仪；
2)向芯片杂交通道泵入750μL体积的延伸试剂，其中，延伸试剂组分为：120U/mLBst DNA聚合酶(NEB，#M0275M)，0.2mM dNTP(dATP、dTTP、dCTP、dGTP各0.2μM的混合物)，1M甜菜碱，20mM的三羟甲基氨基甲烷，10mM的氯化钠，10mM的氯化钾，10mM的硫酸铵，3mM的氯化镁，0.1％的Triton X-100，pH值为8.3；
3)将芯片升温至60±0.5℃，反应10分钟；
4)向芯片杂交通道泵入220μL体积的清洗液1，去除延伸试剂；
5)向芯片杂交通道泵入440μL体积的清洗液2，去除步骤4)中的清洗液1，完成初始模板互补链的合成。
步骤四(可选地)：对步骤三中延伸不完全的新生链的3骤四(进行封闭,封闭的具体步骤如下：
1)将芯片降温至37±0.5℃，维持90秒；
2)向步骤三所述延伸后的通道中泵入750μL体积的封闭试剂1，反应10分钟。所述封闭试剂1的组分为：100U/mL Klenow DNA聚合酶大片段(3′→5′exo-,NEB，#M0212M)12.5μM ddNTP mix(ddATP、ddTTP、ddCTP、ddGTP各12.5μM的混合物)，5mM的氯化锰，20mM的三羟甲基氨基甲烷，10mM的氯化钠，10mM的氯化钾，10mM的硫酸铵，3mM的氯化镁，0.1％的TritonX-100，pH值为8.3；
3)向步骤2)所述封闭后的通道中通入220μL体积的清洗液1，去除封闭反应后剩余的封闭液，完成对延伸不完全的新生链的3'末端的封闭。
步骤五：变性去除初始模板，去除初始模板的过程如下：
1)将芯片降温至55片±0.5℃
2)向步骤四所述封闭后的通道中通入800μL体积的甲酰胺，变性2分钟；
3)向步骤2)所述变性后通道中通入220μL体积的清洗液1，去除变性后的初始模板；
4)重复步骤2)和步骤3)一次，完成对初始模板的去除。
步骤六：对芯片表面残余接头的3'末端：进行封闭，封闭芯片表面残余接头3闭芯片的过程包括：
1)将芯片降温至37±0.5℃；
2)向步骤五所述封闭后的通道中通入440μL体积的清洗液2，替换通道内剩余的清洗液1；
3)向所述步骤2)处理后的通道中通入750μL体积的封闭试剂2，反应15分钟。其中，封闭试剂2的组分为：100U/mL Terminal Transferase(NEB，M0315L)，1×TerminalTransferase Buffer，0.25mM氯化钴，100μM ddNTP mix(ddATP、ddTTP、ddCTP、ddGTP各100μM的混合物)；
4)向步骤3)所述封闭后的通道中通入220μL体积的清洗液1，完成对芯片表面残余接头3完成对的封闭。
步骤七：杂交测序引物D7S1T-R2P，杂交测序引物D7S1T-R2P的过程如下：
1)将芯片升温至55±0.5℃，保持1分钟；
2)向步骤六所述封闭后的通道中通入800μL体积的稀释的测序引物杂交液，杂交反应30分钟。所述稀释的测序引物杂交液为含有0.1μM引物D7S1T-R2P的清洗液3，清洗液3组分包括：450mM的氯化钠，45mM的柠檬酸钠；
3)将芯片降温至37±0.5℃，保持90秒；
4)向步骤2)所述杂交通道中通入220μL体积的清洗液1，去除通道中未被杂交的测序引物；
5)向所述步骤4)处理后的通道中通入440μL体积的清洗液2，替换通道中剩余的清洗液1，完成测序引物的杂交。
步骤八：进行Read1测序，Read1测序的过程如下：
利用Genocare单分子测序平台进行80个循环的测序，测序过程中采用四种核苷酸带有两种不同的荧光信号，每轮反应加入两种标记不同荧光信号的核苷酸进行信号检测的方式进行测序。
步骤九：去除新生测序链。
去除新生测序链的过程按照步骤五中的步骤进行。
步骤十：封闭残余新生链的3'末端。
封闭残余新生链的3'末端的过程按照步骤四中的步骤进行。
步骤十一：杂交所述测序引物D7S1T-R2P。
杂交测序引物D7S1T-R2P的过程按照步骤七中的步骤进行。
步骤十二：进行Read2测序。
Read2测序的过程按照步骤八中的步骤进行。
步骤十三：将测序数据进行拆分获得坐标一一对应的Reads1和Reads2两部分序列。
具体地，本实施例中所述将测序数据进行拆分获得坐标一一对应的Reads1和Reads2两部分序列的过程包括：
使用python语言将160循环测序BaseCalling输出的“.fa_”文件中的每条Read拆分为前80循环和后80循环两个部分，并将所有Reads中的字符“_”移除，分别输出为两份序列坐标一致的“.fa”文件“Reads1.fa”、“Reads2.fa”，，完成将测序数据进行拆分获得坐标一一对应的Reads1和Reads2两部分序列。
进一步地，本实施例中提供的所述一套对上述Two-Pass测序方法获得的Reads1和Reads2进行分析获得Consensus Reads的分析方法包括：
步骤十四：构建校正模型。
具体地，本实施例中所述构建校正模型的过程包括：
1)使用python语言，提取步骤十三获得的Reads1和Reads2序列中同一坐标两次测序读长均≥列中同一坐标的Reads，分别输出为T1(Read1)和T2(Read2)两个文件。其中同一坐标的对应方法是在生成Reads文件时将同一坐标Reads在不同文件中的Reads ID设置为一致；
2)将T1和T2中位置对应的Reads相互间做Align，在Align结果中标记两条Reads一致和不一致的Base，得到Common Reads。其中位置对应是通过比较两条Reads将的Reads ID是否一致实现；
3)分别将文件T1和T2和Reference做Mapping，得到Sam1和Sam2文件。将Sam1和Sam2中位置对应且mapping到同一位置的Reads，找到Reference中最长公共子串RefReads。公共子串指两条对应的Reads mapping后均覆盖的区域；
4)比较步骤2)中的Common Reads和步骤3)中的Ref Reads。对于Common Reads中不一致的Base，标记其是否真实存在于Reference中。若存在，对于没有测到的Reads则为Deletion。若不存在，对于测到的Reads则为Insertion；
5)统计步骤4)中的Deletion和Insertion情况，同时统计该不一致位置上前后Base的种类。因此得到在不同Base类型前或后引起Insertion或Deletion的概率。
具体地，本实例中运用的朴素贝叶斯模型如下：
其中：P(D|XY)表示对于某碱基在前后分别为X和Y碱基时发生Deletion的概率，X，Y∈[A，C，G，T]。P(D)表示对于某碱基发生Deletion的概率；P(I)表示对于某碱基发生Insertion的概率。
通过统计不同碱基下发生Deletion或Insertion时，前后碱基出现频率即可得到P(XY|D)和P(XY|I)，从而可以计算得到P(D|XY)和P(I|XY)。
步骤十五：过滤读长得到Fa1。
具体地，本实施例中所述读长过滤的过程包括：
使用Python语言逐行读取Reads1文件中所有reads，若Reads长度大于等于25bp，则输出的文本文件Fa1中。
步骤十六：根据Reads2读长，将Fa1中Reads进行分类。
具体地，本实例中所述根据Reads2读长分类Fa1中Reads的过程包括：
将Fa1中所有Reads对应在Reads2中的Reads读出，根据Reads2中Reads的长度，若Read2≥25bp，则将对应的Fa1中Reads保存于Fa2文件中；若10bp≤Read2＜25bp，则将对应的Fa1中Reads保存于Fa3文件中。
步骤十七：根据Q值输出置信Reads。
具体地，本实例中所述根据Q值重新输出置信Reads的过程包括：
1)将步骤十六得到的Fa2中所有Reads取出，并同时取出其对应的Reads2中的Reads。从Reads ID中分割得到该Reads的Quality Score值(简称Q值)。
2)比较两条对应的Reads的Q值，将Q值较大的Reads输出到文件Fa4中，Q值较小的Reads输出到文件Fa5中。若两者Q值相等，则默认将Reads1中的Reads输出至Fa4中，Reads2中的Reads输出至Fa5中。
步骤十八：根据Q值过滤Fa4和Fa5中Reads。
具体地，本实例中所述根据Q值过滤Fa4和Fa5中Reads的过程包括：
取Fa4中Reads，根据其Q值，若大于等于60，则将输出到文件Fa6中，同时将该Reads对应的Fa5中的Reads输出到文件Fa7中。
步骤十九：使用Fa7中Reads矫正Fa6中Reads，得到Consensus Reads Parts1(简称CRP1)。
具体地，本实例中所述使用Fa7中的Reads矫正Fa6中Reads的过程包括：
1)取Fa6中Reads和其对应的在Fa7中的Reads。将两条对应Reads相互配准，得到共同的一致性序列部分。其中两条序列配准使用Smith-Waterman算法，一致性序列指配准后通过在序列中增加、删除或修改部分Base，得到的局部最佳匹配序列。
2)得到一致性序列后，根据步骤十四构建的矫正模型，逐个判断一致性序列中不一致的Base位置。根据该Base位置前后的碱基类型计算该位置出现Deletion或Insertion的概率。若Deletion的概率大于50％，则认为该位置所测Base不应该出现，从而删除该位置Base。反之，保留该位置上的Base。
3)矫正所有不一致Base后，输出矫正后的Reads，即为CRP1。这里的不一致Base特指两条对应Reads中没有同时被测出的Base。若两次均测出该Base，但Base类型不一致，不在本实例矫正的候选范围内，该情况下，最终Base类型以Fa6中Reads的Base类型为准。
步骤二十：根据Q值过滤Fa3中Reads。
具体地，本实例中所述根据Q值过滤Fa3中Reads的过程包括：
取Fa3中所有Reads，分割Fa3中Reads的Reads ID，得到每条Reads的Q值。将Q值≥60的Reads输出到文件Fa8中。
步骤二十一：输出Fa8文件中对应的Reads2中的Reads。
具体地，本实例中所述根据Fa8中Reads输出Reads2中Reads的过程包括：
取Fa8文件中所有Reads，取出其对应的Reads2中的Reads，将其输出到文件Fa9中。
步骤二十二：使用Fa9中Reads矫正Fa8中Reads，得到Consensus Reads Parts2(简称CRP2)。
具体地，本实例中所述的使用Fa9中Reads矫正Fa8中Reads的过程参照所述步骤十九进行。
步骤二十三：根据不同应用对应测序数据准确率的需求，将符合相似度阈值的CRP1和CRP2中的Reads合并输出，得到Consensus Reads。
具体地，本实例中所述不同应用对应测序数据准确率的需求，过滤ConsensusReads Part中Reads并输出的过程包括：
1)根据不同应用对应测序数据准确率的需求，设定对应的相似度阈值。其中对Part1和Part2的相似度阈值可以不同；
2)分别计算CRP1和CRP2中的Reads相似度，相似度是指某Reads在Reads1和Reads2中对应的Reads的相似度。相似度计算步骤是先将两条对应Reads相互配准。再计算配准得到的一致性序列中一致的Base数占总Base数的比值。其中配准方法、一致性序列和不一致Base定义参照步骤十九。
3)根据不同应用对应测序数据准确率的需求，分别将CRP1和CRP2中符合相似度阈值要求的Reads输出到最终的文件中，得到Consensus Reads，参考表4。
表4:不同相似度阈值过滤输出序列与参考基因组mapping分析比较
注：数据损失主要发生在读长过滤步骤，由于Read1和Read2测序是相互独立事件，所以必然存在部分读长不一致的序列
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
序列表
<110> 深圳市真迈生物科技有限公司
<120> 测序方法及其分析方法和系统、计算机可读存储介质和电子设备
<130> PIDC4200099
<160> 4
<170> PatentIn version 3.5
<210> 1
<211> 51
<212> DNA
<213> Artificial Sequence
<220>
<223> 引物
<400> 1
ctcagatcct acaacgacgc tctaccgatg aagatgtgta taagagacag t 51
<210> 2
<211> 51
<212> DNA
<213> Artificial Sequence
<220>
<223> 引物
<400> 2
ctgtctctta tacacatctg agtggaactg gatggtcgca ggtatcaagg a 51
<210> 3
<211> 43
<212> DNA
<213> Artificial Sequence
<220>
<223> 引物
<400> 3
ctacaacgac gctctaccga tgaagatgtg tataagagac agt 43
<210> 4
<211> 62
<212> DNA
<213> Artificial Sequence
<220>
<223> 探针序列
<400> 4
tttttttttt tccttgatac ctgcgaccat ccagttccac tcagatgtgt ataagagaca 60
gt 62

Claims (9)

1.一种测序方法，其特征在于，包括：
(1)对芯片表面的测序模板进行第一测序，以便通过形成第一新生测序链获得第一测序数据，所述测序模板通过测序接头连接在所述芯片表面上；
(2)对至少一部分所述第一新生测序链的3’-末端进行第一封闭处理；和
(3)对所述测序模板进行第二测序，以便通过形成第二新生测序链获得第二测序数据；
其中，步骤(2)包括：去除所述芯片表面的所述第一新生测序链，并对残余在所述芯片表面的所述第一新生测序链的3’-末端进行所述第一封闭处理。
2.根据权利要求1所述的测序方法，其特征在于，在(1)之前，包括：
(1-a)使测序文库中的文库分子与芯片表面的测序接头进行杂交；
(1-b)利用所述文库分子作为初始模板，通过合成互补链形成所述测序模板；
(1-c)除去所述初始模板，并对所述芯片表面的核酸分子的3’-末端进行第二封闭处理；
任选地，在(1-c)之前，进一步包括：
(1-b-1)对步骤(1-b)中延伸不完全的所述互补链的3’-末端进行第三封闭处理；
任选地，所述第一封闭处理、所述第二封闭处理和所述第三封闭处理分别独立地通过使3’-末端羟基与延伸反应阻断剂相连而进行的；
任选地，所述延伸反应阻断剂为ddNTP或其衍生物；
任选地，所述第一封闭处理、所述第二封闭处理和所述第三封闭处理分别独立地采用DNA聚合酶和末端转移酶的至少之一进行；
任选地，所述第一封闭处理和所述第三封闭处理分别独立地通过聚合酶连接所述ddNTP或其衍生物，所述第二封闭处理通过所述末端转移酶连接所述ddNTP或其衍生物。
3.一种测序结果分析方法，其特征在于，包括如下步骤：
(1)通过权利要求1～2任一项所述的方法获得第一测序数据和第二测序数据，所述第一测序数据和所述第二测序数据均由多个读段构成，所述第一测序数据中的至少一部分所述读段在所述第二测序数据中存在对应读段，
(2)所述测序结果分析方法包括：
(a)基于所述第一测序数据和所述第二测序数据各自的至少一部分进行相互校正，以便获得最终序列信息。
4.根据权利要求3所述的方法，其特征在于，所述相互校正包括下列步骤：
在所述第一测序数据和所述第二测序数据中选择高质量读段和所述高质量读段的对应读段，所述读段的长度不低于预定长度，所述读段具有不低于预定质量阈值的测序质量；和
将所述高质量读段与所述高质量读段的对应读段进行比对，并基于所述比对结果进行序列信息校正。
5.根据权利要求3所述的方法，其特征在于，(a)包括：
(a-1)根据所述读段的长度，基于所述第一测序数据，构建第一读段集合，所述第一读段集合中的每一个读段长度均不低于第一预定长度；
(a-2)根据所述对应读段的长度，基于所述第一读段集合，构建第二读段集合和第三读段集合，所述第二读段集合中每一个读段的所述对应读段的长度均不低于第二预定长度，所述第三读段集合中每一个读段的所述对应读段的长度均处于预定长度范围内；
(a-3)根据所述第二读段集合中的所述读段及其所述对应读段的测序质量，基于所述第二读段集合及其所述对应读段，构建第四读段集合和第五读段集合，其中，所述第四读段集合和所述第五读段集合分别是按照下列原则确定的：
将所述第二读段集合中的所述读段与其所述对应读段进行测序质量比较，
选择测序质量高的一方作为所述第四读段集合的元素，选择测序质量低的一方作为所述第五读段集合的元素，
对于测序质量相同的情形，则选择来自所述第二读段集合的所述读段作为所述第四读段集合的元素，则选择所述对应读段作为所述第五读段集合的元素；
(a-4)利用测序质量，对所述第四读段集合进行过滤处理，以便构建第六读段集合，所述第六读段集合中的所述读段的测序质量均不低于第一预定质量阈值；
(a-5)利用所述第六读段集合，从所述第五读段集合中选择与所述第六读段集合中的所述读段对应的所述读段，以便构建第七读段集合；
(a-6)将所述第六读段集合与所述第七读段集合进行读段比对，并在所述第六读段集合的所述读段上确定第一差异位点；和
(a-7)利用预先确定的测序误差预测模型，对所述第一差异位点进行校正，以便确定第一序列信息，所述测序误差预测模型用于确定在测序过程中，差异位点发生插入或者缺失的概率；
任选地，所述测序结果分析方法进一步包括：
(a-4a)利用测序质量，对所述第三读段集合进行过滤处理，以便构建第八读段集合，其中，所述第八读段集合中的所述读段的测序质量均不低于第二预定质量阈值；
(a-5a)利用所述第八读段集合，从所述第二测序数据中选择与所述第七读段集合中的所述读段对应的所述读段，以便构建第九读段集合；
(a-6a)将所述第八读段集合与所述第九读段集合进行读段比对，并在所述第八读段集合的所述读段上确定第二差异位点；和
(a-7a)利用所述测序误差预测模型，对所述第二差异位点进行校正，以便确定第二序列信息；
任选地，所述测序误差预测模型是基于所述第一测序数据和所述第二测序数据与参考基因组的比对结果，对朴素贝叶斯模型进行训练获得的；
任选地，针对所述第一差异位点和所述第二差异位点：
如果来自所述第六读段集合的读段在所述差异位点存在碱基，来自所述第七读段集合的对应读段在所述差异位点不存在碱基，并且在所述差异位点发生缺失的概率为50％以上，则保留所述第六读段集合的读段在所述差异位点的碱基作为最终测序结果；
如果来自所述第六读段集合的读段在所述差异位点不存在碱基，来自所述第七读段的读段集合的定读段在所述差异位点存在碱基，并且在所述差异位点发生插入的概率为50％以上，则保留所述第六读段集合的读段在所述差异位点的碱基作为最终测序结果；和
如果来自所述第六读段集合的读段在所述差异位点存在碱基，来自所述第七读段的读段集合的定读段在所述差异位点也存在碱基，则选择所述第六读段集合的读段在所述差异位点的碱基作为最终测序结果；
任选地，所述第一预定长度和所述第二预定长度分别独立地不低于20bp，优选不低于25bp；
任选地，所述预定长度范围为10～25bp；
任选地，所述第一预定质量阈值和所述第二预定质量阈值分别独立地不低于50，优选不低于60。
6.一种测序结果分析系统，其特征在于，包括
测序设备，所述测序设备适于通过权利要求1～2任一项所述的方法获得测序结果，所述测序结果包括第一测序数据和第二测序数据，所述第一测序数据和所述第二测序数据均由多个读段构成，所述第一测序数据中的至少一部分所述读段在所述第二测序数据中存在对应读段；
分析设备，所述分析设备适于基于所述第一测序数据和所述第二测序数据各自的至少一部分进行相互校正，以便获得最终序列信息；
任选地，所述相互校正包括下列步骤：
在所述第一测序数据和所述第二测序数据中选择高质量读段和所述高质量读段的对应读段，所述读段的长度不低于预定长度，所述读段具有不低于预定质量阈值的测序质量；和
将所述高质量读段与所述高质量读段的对应读段进行比对，并基于所述比对结果进行序列信息校正。
7.根据权利要求6所述的测序结果分析系统，其特征在于，所述分析设备进一步包括：
第一读段集合确定模块，根据所述读段的长度，基于所述第一测序数据，构建第一读段集合，所述第一读段集合中的每一个读段长度均不低于第一预定长度；
第二读段集合和第三读段集合确定模块，根据所述对应读段的长度，基于所述第一读段集合，构建第二读段集合和第三读段集合，所述第二读段集合中每一个读段的所述对应读段的长度均不低于第二预定长度，所述第三读段集合中每一个读段的所述对应读段的长度均处于预定长度范围内；
第四读段集合和第五读段集合确定模块，根据所述第二读段集合中所述读段及其所述对应读段的测序质量，基于所述第二读段集合及其所述对应读段，构建第四读段集合和第五读段集合，其中，所述第四读段集合和所述第五读段集合分别是按照下列原则确定的：
将所述第二读段集合中的所述读段与其所述对应读段进行测序质量比较，
选择测序质量高的一方作为所述第四读段集合的元素，选择测序质量低的一方作为所述第五读段集合的元素，
对于测序质量相同的情形，则选择来自所述第二读段集合的所述读段作为所述第四读段集合的元素，则选择所述对应读段作为所述第五读段集合的元素；
第六读段集合确定模块，利用测序质量，对所述第四读段集合进行过滤处理，以便构建第六读段集合，所述第六读段集合中的所述读段的测序质量均不低于第一预定质量阈值；
第七读段集合确定模块，利用所述第六读段集合，从所述第五读段集合中选择与所述第六读段集合中的所述读段对应的所述读段，以便构建第七读段集合；
第一差异位点确定模块，将所述第六读段集合与所述第七读段集合进行读段比对，并在所述第六读段集合的所述读段上确定第一差异位点；和
第一序列信息确定模块，利用预先确定的测序误差预测模型，对所述第一差异位点进行校正，以便确定第一序列信息，所述测序误差预测模型用于确定在测序过程中，差异位点发生插入或者缺失的概率；
任选地，所述测序结果分析系统进一步包括：
第八读段集合确定模块，利用测序质量，对所述第三读段集合进行过滤处理，以便构建第八读段集合，其中，所述第八读段集合中的所述读段的测序质量均不低于第二预定质量阈值；
第九读段集合确定模块，利用所述第八读段集合，从所述第二测序数据中选择与所述第七读段集合中的所述读段对应的所述读段，以便构建第九读段集合；
第二差异位点确定模块，将所述第八读段集合与所述第九读段集合进行读段比对，并在所述第八读段集合的所述读段上确定第二差异位点；
第二序列信息确定模块，利用所述测序误差预测模型，对所述第二差异位点进行校正，以便确定第二序列信息；
任选地，所述测序误差预测模型是基于所述第一测序数据和所述第二测序数据与参考基因组的比对结果，对朴素贝叶斯模型进行训练获得的；
任选地，针对所述第一差异位点和所述第二差异位点：
如果来自所述第六读段集合的读段在所述差异位点存在碱基，来自所述第七读段集合的对应读段在所述差异位点不存在碱基，并且在所述差异位点发生缺失的概率为50％以上，则保留所述第六读段集合的读段在所述差异位点的碱基作为最终测序结果；
如果来自所述第六读段集合的读段在所述差异位点不存在碱基，来自所述第七读段的读段集合的定读段在所述差异位点存在碱基，并且在所述差异位点发生插入的概率为50％以上，则保留所述第六读段集合的读段在所述差异位点的碱基作为最终测序结果；和
如果来自所述第六读段集合的读段在所述差异位点存在碱基，来自所述第七读段的读段集合的定读段在所述差异位点也存在碱基，则选择所述第六读段集合的读段在所述差异位点的碱基作为最终测序结果；
任选地，所述第一预定长度和所述第二预定长度分别独立地不低于20bp，优选不低于25bp；
任选地，所述预定长度范围为10～25bp；
任选地，所述第一预定质量阈值和所述第二预定质量阈值分别独立地不低于50，优选不低于60。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1～2中任一项所述方法的步骤。
9.一种电子设备，其特征在于，包括：
权利要求8中所述的计算机可读存储介质；以及
一个或者多个处理器，用于执行所述计算机可读存储介质中的程序。
CN202110480833.2A 2020-04-30 2021-04-30 测序方法及其分析方法和系统、计算机可读存储介质和电子设备 Active CN113593637B (zh)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN2020103625876 2020-04-30
CN202010362587 2020-04-30
CN2020108675693 2020-08-25
CN202010867569 2020-08-25

Publications (2)

Publication Number Publication Date
CN113593637A CN113593637A (zh) 2021-11-02
CN113593637B true CN113593637B (zh) 2024-05-03

Family

ID=78243116

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110480833.2A Active CN113593637B (zh) 2020-04-30 2021-04-30 测序方法及其分析方法和系统、计算机可读存储介质和电子设备

Country Status (4)

Country Link
US (1) US20230178183A1 (zh)
EP (1) EP4144745A4 (zh)
CN (1) CN113593637B (zh)
WO (1) WO2021219114A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024093421A1 (zh) * 2022-10-31 2024-05-10 深圳市真迈生物科技有限公司 表面处理方法、测序方法和试剂盒

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105624272A (zh) * 2014-10-29 2016-06-01 深圳华大基因科技有限公司 基因组预定区域核酸测序文库的构建方法及装置
CN108165618A (zh) * 2017-12-08 2018-06-15 东南大学 一种包含核苷酸和3’端可逆封闭核苷酸的dna测序方法
CN109207571A (zh) * 2017-06-30 2019-01-15 深圳华大基因研究院 一种检测核酸内切酶酶切位点的方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0718409A3 (en) * 1994-12-22 1999-02-03 Hitachi, Ltd. DNA preparation method
WO2007123744A2 (en) * 2006-03-31 2007-11-01 Solexa, Inc. Systems and devices for sequence by synthesis analysis
CN102296065B (zh) * 2011-08-04 2013-05-15 盛司潼 用于构建测序文库的系统与方法
CN103543270B (zh) * 2012-07-09 2015-07-08 国家纳米科学中心 一种蛋白质原位表达芯片及其制备方法和应用
EP3434776A1 (en) * 2012-12-12 2019-01-30 The Broad Institute, Inc. Methods, models, systems, and apparatus for identifying target sequences for cas enzymes or crispr-cas systems for target sequences and conveying results thereof
US10435736B2 (en) * 2014-12-18 2019-10-08 Mgi Tech Co., Ltd. Target region enrichment method based on multiplex PCR, and reagent
US10095831B2 (en) * 2016-02-03 2018-10-09 Verinata Health, Inc. Using cell-free DNA fragment size to determine copy number variations
CN107227344B (zh) * 2017-05-16 2019-09-24 北京量化健康科技有限公司 一种高密度、高稳定性的核酸芯片及其制备方法
JP7299169B2 (ja) * 2017-12-01 2023-06-27 イルミナ インコーポレイテッド 体細胞突然変異のクローン性を決定するための方法及びシステム
CN113913945A (zh) 2018-10-12 2022-01-11 深圳市真迈生物科技有限公司 制备芯片的方法、芯片以及固定生物分子的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105624272A (zh) * 2014-10-29 2016-06-01 深圳华大基因科技有限公司 基因组预定区域核酸测序文库的构建方法及装置
CN109207571A (zh) * 2017-06-30 2019-01-15 深圳华大基因研究院 一种检测核酸内切酶酶切位点的方法
CN108165618A (zh) * 2017-12-08 2018-06-15 东南大学 一种包含核苷酸和3’端可逆封闭核苷酸的dna测序方法

Also Published As

Publication number Publication date
CN113593637A (zh) 2021-11-02
EP4144745A1 (en) 2023-03-08
WO2021219114A1 (zh) 2021-11-04
EP4144745A4 (en) 2024-02-21
US20230178183A1 (en) 2023-06-08

Similar Documents

Publication Publication Date Title
CN113593636B (zh) 测序结果分析方法、系统及计算机可读存储介质和电子设备
EP3763825B1 (en) High multiplex pcr with molecular barcoding
CN108300716B (zh) 接头元件、其应用和基于不对称多重pcr进行靶向测序文库构建的方法
US11702690B2 (en) Sequencing library, and preparation and use thereof
EP3885445B1 (en) Methods of attaching adapters to sample nucleic acids
CN107604046B (zh) 用于微量dna超低频突变检测的双分子自校验文库制备及杂交捕获的二代测序方法
JP2022510723A (ja) 遺伝子標的エリアの富化方法及びキット
WO2018184495A1 (zh) 一步法构建扩增子文库的方法
CN111041563A (zh) 一种靶向序列捕获及pcr建库方法
CN112251821A (zh) 一种快速高效的构建二代测序文库的试剂盒
CN112941635A (zh) 一种提高文库转化率的二代测序建库试剂盒及其方法
CN110603327A (zh) Pcr引物对及其应用
CN113593637B (zh) 测序方法及其分析方法和系统、计算机可读存储介质和电子设备
CN111647953A (zh) 用于检测地中海贫血基因突变的高通量文库构建试剂盒及文库构建方法
CN116024308A (zh) 肺癌相关基因高通量扩增子文库的制备方法、多重pcr引物对及应用
CN109686404B (zh) 检测样本混淆的方法及装置
CN112301430B (zh) 一种建库方法及应用
CN114807300A (zh) 单引物多重扩增技术在检测片段化稀有特征核酸分子中的应用及试剂盒
CN110468211B (zh) 膀胱癌肿瘤突变基因特异性引物、试剂盒和文库构建方法
CN110452958B (zh) 一种微量碎片化核酸甲基化检测的接头、引物、试剂盒及其应用
US20180100180A1 (en) Methods of single dna/rna molecule counting
WO2023202030A1 (zh) 一种小分子rna的高通量测序文库构建方法
WO2023066255A1 (zh) 测序方法、测序数据处理方法、设备和计算机设备
CN114277114B (zh) 一种扩增子测序添加唯一性标识符的方法及应用
CN106701949B (zh) 一种减少扩增偏倚的基因突变检测方法和试剂

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
REG Reference to a national code

Ref country code: HK

Ref legal event code: DE

Ref document number: 40064482

Country of ref document: HK

GR01 Patent grant
GR01 Patent grant