CN111370063A - 一种基于Pacbio数据的MSI微卫星不稳定性检测方法和系统 - Google Patents

一种基于Pacbio数据的MSI微卫星不稳定性检测方法和系统 Download PDF

Info

Publication number
CN111370063A
CN111370063A CN202010207692.2A CN202010207692A CN111370063A CN 111370063 A CN111370063 A CN 111370063A CN 202010207692 A CN202010207692 A CN 202010207692A CN 111370063 A CN111370063 A CN 111370063A
Authority
CN
China
Prior art keywords
sequence
msi
sequences
motifs
core
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.)
Granted
Application number
CN202010207692.2A
Other languages
English (en)
Other versions
CN111370063B (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.)
Shanghai Oe Biotech Co ltd
Original Assignee
Shanghai Oe Biotech 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 Shanghai Oe Biotech Co ltd filed Critical Shanghai Oe Biotech Co ltd
Priority to CN202010207692.2A priority Critical patent/CN111370063B/zh
Publication of CN111370063A publication Critical patent/CN111370063A/zh
Application granted granted Critical
Publication of CN111370063B publication Critical patent/CN111370063B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biophysics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Biotechnology (AREA)
  • Evolutionary Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明提出了一种基于Pacbio数据的MSI微卫星不稳定性检测方法,包括以下步骤:步骤一:测序错误的校正,利用ccs软件将Pacbio测序结果的subreads校正为CCS序列,min‑passes>=10;步骤二,检测MSI‑PCR上下游引物的位置,并提取产物序列;步骤三,根据MSI核心区域旁邻20bp的特异性序列,进一步提取产物序列中的核心序列;步骤四,旁邻特异性序列是反向比对的,需要对核心序列进行反向互补转化;步骤五,根据给定的microsatellite motif,检测产物序列或核心序列中motif的重复次数以及motif之间的嵌合序列;步骤六,筛选重复次数>=3的motif,保留它们之间的嵌合序列,作为MSI序列;步骤七,统计全部MSI序列的频次和频率,过滤频率小于5%的MSI序列;考虑到肿瘤的异质性,肿瘤组织样本可以放宽过滤阈值。

Description

一种基于Pacbio数据的MSI微卫星不稳定性检测方法和系统
技术领域
本发明涉及数据处理技术领域,涉及人类微卫星不稳定性状态MSI的检测,利用Pacbio测序数据分析人类MSI的状态。
背景技术
PacBio平台(PacBio核酸测序平台)采用的是单分子实时测序(SingleMoleculeReal-Time)。PacBio平台的测序基于两个关键技术,第一,Pacific Biosciences公司发明的一种直径只有几十纳米的纳米孔,该纳米孔内只能容许单分子核酸在DNA聚合酶作用下进行反应,检测纳米孔中合成过程的A、T、C、G这四种荧光标记的脱氧核苷酸,即可实现核酸测序;第二,利用共聚焦显微镜实时、快速地对集成在SMRT cell上的无数的纳米孔同时进行记录。PacBio核酸测序平台,其平均读长可达10-40k,且测序不受AT和CG含量的影响,能够对高GC含量或低GC含量的区域进行测序,相比二代测序有较大优势。此外,Pacbio技术可以对一段DNA进行循环测序,从而拿到准确度高的测序结果,解决了三代测序碱基准确性低的问题。
微卫星(microsatellite)是指DNA基因组中小于10个核苷酸的简单重复序列,遍布于整个人类基因组,一般为2-6个碱基重复。这些高度重复序列在机体自我复制时容易出现碱基错配或跳读的现象,因此需要机体的错配修复系统来进行校正。错配修复系统是指一系列可以矫正碱基错配的蛋白,可维持基因组的稳定及降低自发性突变;当错配修复系统发生失调,则发生在微卫星区域高发的体细胞突变不能被校正,进而形成微卫星重复单位的减少或增加而造成微卫星长度改变,形成新的微卫星等位基因,造成微卫星不稳定(microsatellite instability,MSI)。
MSI的发生显著提高了基因的随机突变率,是基因组不稳定的重要指标,并与林奇综合征(Lynch Syndrome)等遗传病相关。MSI的状态也是多种癌症类型的独立预后因素,可指导临床用药及治疗方案的选择,在结直肠癌和子宫内膜癌患者中表现尤为明显。
MSI-PCR通过对肿瘤组织和正常组织中选定的微卫星位点进行PCR扩增及凝胶电泳,通过比较两组电泳结果的差异来确定MSI的状态。然而这种方式只能看到MSI的长度差异,准确的microsatellite重复类型和重复次数无法统计。而且基因组中有数以万计的微卫星位点,不同的位点对于检测MSI的敏感性和准确性也各不相同,其中包含重复次数极高MSI,这类MSI的PCR产物长度远超过一代和二代测序技术的测序长度,它们MSI状态的鉴定成为难题。
发明内容
本发明的目的在于提供一种基于Pacbio数据的MSI微卫星不稳定性检测方法,以解决长片段MSI状态无法精确鉴定的问题。
为了实现上述目的,本发明提出了一种基于Pacbio数据的MSI微卫星不稳定性检测方法,包括以下步骤:
步骤一,测序错误的校正,利用ccs软件将Pacbio测序数据校正为CCS序列;为保证校正后的测序错误率低于0.1%,设置参数min-passes>=10。
步骤二,检测MSI-PCR上下游引物的位置,并提取MSI-PCR产物序列;为保证步骤一获取的CCS序列的可靠性,需要检测CCS序列中是否包含引物序列,通过bwa aln软件默认参数将引物序列比对到CCS序列上,检测引物位置,根据引物位置提取产物序列,过滤不含有引物序列的CCS序列。
步骤三,根据MSI核心区域上下游20bp特异性序列,进一步提取核心序列;所述核心序列是MSI-PCR产物中microsatellite motif重复序列和motif之间的嵌合序列;通过核心序列两侧相邻的20bp特异性序列,利用bwa aln软件默认参数将特异性序列比对到步骤二获得的MSI-PCR产物序列上,检测特异性序列位置,根据特异性序列位置提取核心序列,过滤不含有特异性序列的MSI-PCR产物序列。
步骤四,如果步骤三旁邻特20bp异性序列是反向比对的,需要对所述核心序列进行反向互补转化为正向序列。
步骤五,根据给定的microsatellite motif,检测产物序列或核心序列中motif的重复次数以及motif之间的嵌合序列;检测microsatellite motif重复序列的重复次数,通过匹配核心序列中存在的motif,统计motif的重复次数以及motif之间的嵌合序列。
步骤六,筛选重复次数>=3的motif,保留它们之间的嵌合序列,作为MSI序列;
步骤七,统计全部MSI序列的频次和频率,过滤频率小于5%的MSI序列;考虑到肿瘤组织的异质性,肿瘤组织样本可以放宽过滤阈值,可大于5%。
基于以上方法,本发明还提出了一种基于Pacbio数据的MSI微卫星不稳定性检测系统,其包括:
测序校正模块,其利用ccs软件将Pacbio测序数据校正为CCS序列;
产物序列提取模块,其用于检测MSI-PCR上下游引物的位置,并提取MSI-PCR产物序列;
核心序列提取模块,其用于根据MSI核心区域上下游20bp特异性序列,进一步提取核心序列;
反向序列转化模块,其用于对所述核心序列进行反向互补转化;
统计模块,其用于根据给定的microsatellite motif,检测产物序列或核心序列中motif的重复次数以及motif之间的嵌合序列;
筛选模块,其用于筛选重复次数>=3的motif,保留它们之间的嵌合序列,作为MSI序列;过滤模块,其用于统计全部MSI序列的频次和频率,过滤频率小于阈值的MSI序列。
本发明提出的基于Pacbio数据的MSI微卫星不稳定性检测方法,针对MSI长度<=20kb的序列,都可以通过Pacbio测序结果,分析出准确的MSI序列状态,解决了相关技术中MSI检测敏感性低,检测长度受限的问题。
附图说明
图1是根据本发明实施例的用于Pacbio的数据处理的示意图。
图2是本发明基于Pacbio数据的MSI微卫星不稳定性检测系统的示意图。
具体实施方式
结合以下具体实施例和附图,对发明作进一步的详细说明。实施本发明的过程、条件、实验方法等,除以下专门提及的内容之外,均为本领域的普遍知识和公知常识,本发明没有特别限制内容。
本发明的一个方面,提供了一种PacBio测序数据的处理方法,该方法包括:
步骤一,测序错误的校正,由于Pacbio测序数据存在15%左右的测序错误,测序错误会影响数据的比对和MSI的识别,需要利用ccs软件将Pacbio测序数据校正为CCS序列,为保证校正后的测序错误率低于0.1%,设置参数min-passes>=10;
步骤二,检测MSI-PCR上下游引物的位置,提取MSI-PCR产物序列;为保证步骤一获取的CCS序列的可靠性,需要检测CCS序列中是否包含引物序列,通过bwa aln软件默认参数将引物序列比对到CCS序列上,检测引物位置,根据引物位置提取产物序列,过滤不含有引物序列的CCS序列;
步骤三,根据MSI核心区域上下游20bp特异性序列,进一步提取核心序列;核心序列是MSI-PCR产物中microsatellite motif重复序列和motif之间的嵌合序列。通过核心序列两侧相邻的20bp特异性序列,利用bwa aln软件默认参数将特异性序列比对到步骤二获得的MSI-PCR产物序列上,检测特异性序列位置,根据特异性序列位置提取核心序列,过滤不含有特异性序列的MSI-PCR产物序列;
步骤四,反向序列的转化,如果步骤三核心序列两侧相邻的20bp特异性序列的比对方向为负,则需要将提取的核心序列进行反向互补,转化为正向序列;
本发明的另一个方面,提供了一种检测核心序列中MSI类型和数目的数据处理方法,该方法包括:
步骤一,检测microsatellite motif重复序列的重复次数,通过匹配核心序列中存在的motif,统计motif的重复次数以及motif之间的嵌合序列;
步骤二,MSI序列的筛选,筛选重复次数>=3的motif,保留它们之间的嵌合序列,作为MSI序列;
步骤三,MSI序列的过滤,由于PCR错位,测序错误等原因,MSI序列中依然存在假阳性的结果,需要统计全部MSI序列的频率,过滤频率小于阈值的MSI序列。所述阈值优选为5%;考虑到肿瘤组织的异质性,肿瘤组织样本可以放宽过滤阈值。
本发明提供的基于Pacbio数据的MSI微卫星不稳定性检测系统,其包括:
测序校正模块,其利用ccs软件将Pacbio测序数据校正为CCS序列;
产物序列提取模块,其用于检测MSI-PCR上下游引物的位置,并提取MSI-PCR产物序列;
核心序列提取模块,其用于根据MSI核心区域上下游20bp特异性序列,进一步提取核心序列;
反向序列转化模块,其用于对所述核心序列进行反向互补转化;
统计模块,其用于根据给定的microsatellite motif,检测产物序列或核心序列中motif的重复次数以及motif之间的嵌合序列;
筛选模块,其用于筛选重复次数>=3的motif,保留它们之间的嵌合序列,作为MSI序列;过滤模块,其用于统计全部MSI序列的频次和频率,过滤频率小于阈值的MSI序列。
实施例1,人类HTT基因CAG,CCG motif检测
本实施例测序数据包括1个MSI-PCR的50M Pacbio测序数据。
(1)利用ccs软件合并Pacbio测序数据,共获得5347条CCS序列
(2)通过bwa aln比对引物序列,提取MSI-PCR产物序列,共获得5104条MSI-PCR产物序列
(3)通过bwa aln比对核心序列两侧相邻的20bp特异性序列,提取核心序列,共获得5101条核心序列
(4)特异性序列比对方向为负的核心序列,转化为反向互补序列
(5)匹配CAG,CCG类型的重复序列位置,统计重复次数。
(6)筛选重复次数>=3的motif,保留它们之间的嵌合序列,作为MSI序列。
(7)统计MSI序列的频次和频率,结果如下:
Figure BDA0002421717620000041
Figure BDA0002421717620000051
过滤频率<=5%的结果,保留结果如下:
Figure BDA0002421717620000052
基于上述结果,本发明可以有效检测多种类型的嵌合MSI motif,且能准确预测杂合基因型。
实施例2,人类TBP基因CAG,CAA motif检测
本实施案例测序数据包括1个MSI-PCR的80M Pacbio测序数据。
(1)利用ccs软件合并Pacbio测序数据,共获得6290条CCS序列
(2)通过bwa aln比对引物序列,提取MSI-PCR产物序列,共获得6005条MSI-PCR产物序列
(3)通过bwa aln比对核心序列两侧相邻的20bp特异性序列,提取核心序列,共获得5992条核心序列
(4)特异性序列比对方向为负的核心序列,转化为反向互补序列
(5)匹配CAG,CCG类型的重复序列,统计重复次数
(6)筛选重复次数>=3的motif,保留它们之间的嵌合序列,作为MSI序列。
(7)统计MSI序列的频次和频率,结果如下:
Figure BDA0002421717620000053
Figure BDA0002421717620000061
过滤频率<=5%的结果,保留结果如下:
Figure BDA0002421717620000062
基于上述结果,本发明可以有效检测多种类型的嵌合MSI motif,且能准确预测杂合基因型
本发明的保护内容不局限于以上实施例。在不背离发明构思的精神和范围下,本领域技术人员能够想到的变化和优点都被包括在本发明中,并且以所附的权利要求书为保护范围。

Claims (8)

1.一种基于Pacbio数据的MSI微卫星不稳定性检测方法,其特征在于,包括以下步骤:
步骤一,测序错误的校正,利用ccs软件将Pacbio测序数据校正为CCS序列;
步骤二,检测MSI-PCR上下游引物的位置,并提取MSI-PCR产物序列;
步骤三,根据MSI核心区域上下游20bp特异性序列,进一步提取核心序列;
步骤四,旁邻特异性序列是反向比对的,需要对所述核心序列进行反向互补转化;
步骤五,根据给定的microsatellite motif,检测产物序列或核心序列中motif的重复次数以及motif之间的嵌合序列;
步骤六,筛选重复次数>=3的motif,保留它们之间的嵌合序列,作为MSI序列;
步骤七,统计全部MSI序列的频次和频率,过滤频率小于阈值的MSI序列。
2.如权利要求1所述的基于Pacbio数据的MSI微卫星不稳定性检测方法,其特征在于,步骤一中,为保证校正后的测序错误率低于0.1%,设置参数min-passes>=10。
3.如权利要求1所述的基于Pacbio数据的MSI微卫星不稳定性检测方法,其特征在于,步骤二中,为保证步骤一获取的CCS序列的可靠性,需要检测CCS序列中是否包含引物序列,通过bwa aln软件默认参数将引物序列比对到CCS序列上,检测引物位置,根据引物位置提取产物序列,过滤不含有引物序列的CCS序列。
4.如权利要求1所述的基于Pacbio数据的MSI微卫星不稳定性检测方法,其特征在于,步骤三中,所述核心序列是MSI-PCR产物中microsatellite motif重复序列和motif之间的嵌合序列;通过核心序列两侧相邻的20bp特异性序列,利用bwa aln软件默认参数将特异性序列比对到步骤二获得的MSI-PCR产物序列上,检测特异性序列位置,根据特异性序列位置提取核心序列,过滤不含有特异性序列的MSI-PCR产物序列。
5.如权利要求1所述的基于Pacbio数据的MSI微卫星不稳定性检测方法,其特征在于,步骤四中,若步骤三中核心序列两侧相邻的20bp特异性序列的比对方向为负,则需要将提取的核心序列进行反向互补,转化为正向序列。
6.如权利要求1所述的基于Pacbio数据的MSI微卫星不稳定性检测方法,其特征在于,步骤五中,检测microsatellite motif重复序列的重复次数,通过匹配核心序列中存在的motif,统计motif的重复次数以及motif之间的嵌合序列。
7.如权利要求1所述的基于Pacbio数据的MSI微卫星不稳定性检测方法,其特征在于,步骤七中,所述阈值为5%;或,若为肿瘤组织样本,则阈值的取值范围为大于5%。
8.一种基于Pacbio数据的MSI微卫星不稳定性检测系统,其特征在于,采用如权利要求1-7之任一项所述的基于Pacbio数据的MSI微卫星不稳定性检测方法,所述系统包括:
测序校正模块,其利用ccs软件将Pacbio测序数据校正为CCS序列;
产物序列提取模块,其用于检测MSI-PCR上下游引物的位置,并提取MSI-PCR产物序列;
核心序列提取模块,其用于根据MSI核心区域上下游20bp特异性序列,进一步提取核心序列;
反向序列转化模块,其用于对所述核心序列进行反向互补转化;
统计模块,其用于根据给定的microsatellite motif,检测产物序列或核心序列中motif的重复次数以及motif之间的嵌合序列;
筛选模块,其用于筛选重复次数>=3的motif,保留它们之间的嵌合序列,作为MSI序列;
过滤模块,其用于统计全部MSI序列的频次和频率,过滤频率小于阈值的MSI序列。
CN202010207692.2A 2020-03-23 2020-03-23 一种基于Pacbio数据的MSI微卫星不稳定性检测方法和系统 Active CN111370063B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010207692.2A CN111370063B (zh) 2020-03-23 2020-03-23 一种基于Pacbio数据的MSI微卫星不稳定性检测方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010207692.2A CN111370063B (zh) 2020-03-23 2020-03-23 一种基于Pacbio数据的MSI微卫星不稳定性检测方法和系统

Publications (2)

Publication Number Publication Date
CN111370063A true CN111370063A (zh) 2020-07-03
CN111370063B CN111370063B (zh) 2021-09-10

Family

ID=71209026

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010207692.2A Active CN111370063B (zh) 2020-03-23 2020-03-23 一种基于Pacbio数据的MSI微卫星不稳定性检测方法和系统

Country Status (1)

Country Link
CN (1) CN111370063B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112037859A (zh) * 2020-09-02 2020-12-04 迈杰转化医学研究(苏州)有限公司 一种微卫星不稳定性的分析方法和分析装置

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101501193A (zh) * 2006-08-11 2009-08-05 普罗森那技术公司 用于治疗与dna重复不稳定性相关的遗传病的方法和手段
CN102230004A (zh) * 2011-06-08 2011-11-02 北京阅微基因技术有限公司 肿瘤细胞微卫星不稳定状态的复合扩增体系及检测试剂盒
CN103555843A (zh) * 2013-11-05 2014-02-05 上海赛安生物医药科技有限公司 结直肠癌微卫星不稳定性扩增体系及其检测试剂盒
CN104878086A (zh) * 2009-02-11 2015-09-02 卡里斯Mpi公司 肿瘤的分子谱分析
CN105916983A (zh) * 2013-10-25 2016-08-31 塞勒克提斯公司 用于高效且特异性靶向包含高度重复基序的dna序列的稀有切割核酸内切酶的设计
CN107058551A (zh) * 2017-05-04 2017-08-18 北京诺禾致源科技股份有限公司 检测微卫星位点不稳定性的方法及装置
CN107217103A (zh) * 2017-07-14 2017-09-29 常州桐树生物科技有限公司 检测微卫星不稳定性的多重荧光pcr扩增试剂及试剂盒
CN107475385A (zh) * 2017-08-21 2017-12-15 上海派森诺生物科技股份有限公司 一种基于smrt高通量测序技术的菌群多样性组成谱数据分析方法
CN108588239A (zh) * 2018-06-05 2018-09-28 浙江大学 朱鹮微卫星多态位点、鉴定方法及特异引物序列
CN108588240A (zh) * 2018-06-05 2018-09-28 浙江大学 扬子鳄微卫星多态位点、鉴定方法及特异引物序列
CN109182525A (zh) * 2018-09-29 2019-01-11 广州燃石医学检验所有限公司 一种微卫星生物标志物组合、检测试剂盒及其用途
CN109207594A (zh) * 2018-09-29 2019-01-15 广州燃石医学检验所有限公司 一种基于二代测序的通过血浆检测微卫星稳定状态和基因组变化的方法
CN109584961A (zh) * 2018-12-03 2019-04-05 元码基因科技(北京)股份有限公司 基于二代测序技术检测血液微卫星不稳定的方法
CN109949862A (zh) * 2019-03-13 2019-06-28 拓普基因科技(广州)有限责任公司 一种血液ctDNA的微卫星不稳定性检测方法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101501193A (zh) * 2006-08-11 2009-08-05 普罗森那技术公司 用于治疗与dna重复不稳定性相关的遗传病的方法和手段
CN104878086A (zh) * 2009-02-11 2015-09-02 卡里斯Mpi公司 肿瘤的分子谱分析
CN102230004A (zh) * 2011-06-08 2011-11-02 北京阅微基因技术有限公司 肿瘤细胞微卫星不稳定状态的复合扩增体系及检测试剂盒
CN105916983A (zh) * 2013-10-25 2016-08-31 塞勒克提斯公司 用于高效且特异性靶向包含高度重复基序的dna序列的稀有切割核酸内切酶的设计
CN103555843A (zh) * 2013-11-05 2014-02-05 上海赛安生物医药科技有限公司 结直肠癌微卫星不稳定性扩增体系及其检测试剂盒
CN107058551A (zh) * 2017-05-04 2017-08-18 北京诺禾致源科技股份有限公司 检测微卫星位点不稳定性的方法及装置
CN107217103A (zh) * 2017-07-14 2017-09-29 常州桐树生物科技有限公司 检测微卫星不稳定性的多重荧光pcr扩增试剂及试剂盒
CN107475385A (zh) * 2017-08-21 2017-12-15 上海派森诺生物科技股份有限公司 一种基于smrt高通量测序技术的菌群多样性组成谱数据分析方法
CN108588239A (zh) * 2018-06-05 2018-09-28 浙江大学 朱鹮微卫星多态位点、鉴定方法及特异引物序列
CN108588240A (zh) * 2018-06-05 2018-09-28 浙江大学 扬子鳄微卫星多态位点、鉴定方法及特异引物序列
CN109182525A (zh) * 2018-09-29 2019-01-11 广州燃石医学检验所有限公司 一种微卫星生物标志物组合、检测试剂盒及其用途
CN109207594A (zh) * 2018-09-29 2019-01-15 广州燃石医学检验所有限公司 一种基于二代测序的通过血浆检测微卫星稳定状态和基因组变化的方法
CN109584961A (zh) * 2018-12-03 2019-04-05 元码基因科技(北京)股份有限公司 基于二代测序技术检测血液微卫星不稳定的方法
CN109949862A (zh) * 2019-03-13 2019-06-28 拓普基因科技(广州)有限责任公司 一种血液ctDNA的微卫星不稳定性检测方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112037859A (zh) * 2020-09-02 2020-12-04 迈杰转化医学研究(苏州)有限公司 一种微卫星不稳定性的分析方法和分析装置
CN112037859B (zh) * 2020-09-02 2023-12-19 迈杰转化医学研究(苏州)有限公司 一种微卫星不稳定性的分析方法和分析装置

Also Published As

Publication number Publication date
CN111370063B (zh) 2021-09-10

Similar Documents

Publication Publication Date Title
JP5389638B2 (ja) 制限断片に基づく分子マーカーのハイスループットな検出
JP7369492B2 (ja) 核酸の塩基修飾の決定
US9394567B2 (en) Detection and quantification of sample contamination in immune repertoire analysis
EP2663655B1 (en) Paired end random sequence based genotyping
AU2013246050B2 (en) Detection and quantitation of sample contamination in immune repertoire analysis
CN107077537A (zh) 用短读测序数据检测重复扩增
AU2021269294A1 (en) Validation methods and systems for sequence variant calls
CN111108218B (zh) 使用压缩的分子标记的核酸序列数据检测融合的方法
AU2013274366A1 (en) Method of sequence determination using sequence tags
WO2011074960A1 (en) Restriction enzyme based whole genome sequencing
CN105385755A (zh) 一种利用多重pcr技术进行snp-单体型分析的方法
CA3075932A1 (en) Methods and systems for differentiating somatic and germline variants
CN111370063B (zh) 一种基于Pacbio数据的MSI微卫星不稳定性检测方法和系统
CN108359723B (zh) 一种降低深度测序错误的方法
EP2510114B1 (en) Rna analytics method
CN112513292B (zh) 基于高通量测序检测同源序列的方法和装置
EP3474168B1 (en) Method for measuring mutation rate
US20200075124A1 (en) Methods and systems for detecting allelic imbalance in cell-free nucleic acid samples
EP3795685A1 (en) Methods for dna library generation to facilitate the detection and reporting of low frequency variants
EP3596229A1 (en) Method and system for nucleic acid sequencing
CN109817277B (zh) 基于PacBio全长转录组测序数据的质控方法
CN115896303B (zh) 用于大熊猫性别鉴定的snp引物组及应用
WO2017136606A1 (en) Apparatus, systems, and methods for dna amplification with post-sequencing data filtering and cell isolation
CN116705153A (zh) 确定snp检测区域的方法和对测序样本进行校正的方法
CN115631801A (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
CB03 Change of inventor or designer information

Inventor after: Gu Yincong

Inventor after: Xiao Yunping

Inventor after: Wang Shuwei

Inventor after: Shi Xianjun

Inventor after: Lin Bo

Inventor after: Liu Yuchuan

Inventor before: Wang Shuwei

Inventor before: Xiao Yunping

Inventor before: Shi Xianjun

Inventor before: Lin Bo

Inventor before: Gu Yincong

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A method and system for MSI microsatellite instability detection based on Pacbio data

Effective date of registration: 20220831

Granted publication date: 20210910

Pledgee: Shanghai Rural Commercial Bank Co.,Ltd. Xuhui sub branch

Pledgor: SHANGHAI OE BIOTECH CO.,LTD.

Registration number: Y2022310000222

PE01 Entry into force of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Granted publication date: 20210910

Pledgee: Shanghai Rural Commercial Bank Co.,Ltd. Xuhui sub branch

Pledgor: SHANGHAI OE BIOTECH CO.,LTD.

Registration number: Y2022310000222

PC01 Cancellation of the registration of the contract for pledge of patent right