CN109994155A - 一种基因变异识别方法、装置和存储介质 - Google Patents
一种基因变异识别方法、装置和存储介质 Download PDFInfo
- Publication number
- CN109994155A CN109994155A CN201910251891.0A CN201910251891A CN109994155A CN 109994155 A CN109994155 A CN 109994155A CN 201910251891 A CN201910251891 A CN 201910251891A CN 109994155 A CN109994155 A CN 109994155A
- Authority
- CN
- China
- Prior art keywords
- gene
- genetic mutation
- sequence
- site
- sequencing read
- 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
Links
- 230000035772 mutation Effects 0.000 title claims abstract description 277
- 238000000034 method Methods 0.000 title claims abstract description 55
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 557
- 238000012163 sequencing technique Methods 0.000 claims abstract description 275
- 238000003780 insertion Methods 0.000 claims description 22
- 230000037431 insertion Effects 0.000 claims description 22
- 238000012512 characterization method Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 9
- 230000002708 enhancing effect Effects 0.000 abstract description 2
- 230000010354 integration Effects 0.000 description 27
- 230000008569 process Effects 0.000 description 18
- 238000003062 neural network model Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 15
- 238000012545 processing Methods 0.000 description 14
- 210000001082 somatic cell Anatomy 0.000 description 14
- 238000005516 engineering process Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 101150096316 5 gene Proteins 0.000 description 7
- 230000002159 abnormal effect Effects 0.000 description 7
- 239000000284 extract Substances 0.000 description 7
- 235000013399 edible fruits Nutrition 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 229920002477 rna polymer Polymers 0.000 description 4
- 238000012216 screening Methods 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 3
- 210000004602 germ cell Anatomy 0.000 description 3
- 238000012549 training Methods 0.000 description 3
- 101150028074 2 gene Proteins 0.000 description 2
- 101150090724 3 gene Proteins 0.000 description 2
- 241001269238 Data Species 0.000 description 2
- IQFYYKKMVGJFEH-XLPZGREQSA-N Thymidine Chemical compound O=C1NC(=O)C(C)=CN1[C@@H]1O[C@H](CO)[C@@H](O)C1 IQFYYKKMVGJFEH-XLPZGREQSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012864 cross contamination Methods 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- UYTPUPDQBNUYGX-UHFFFAOYSA-N guanine Chemical compound O=C1NC(N)=NC2=C1N=CN2 UYTPUPDQBNUYGX-UHFFFAOYSA-N 0.000 description 2
- 239000002773 nucleotide Substances 0.000 description 2
- 125000003729 nucleotide group Chemical group 0.000 description 2
- 238000007637 random forest analysis Methods 0.000 description 2
- 229930024421 Adenine Natural products 0.000 description 1
- GFFGJBXGBJISGV-UHFFFAOYSA-N Adenine Chemical compound NC1=NC=NC2=C1N=CN2 GFFGJBXGBJISGV-UHFFFAOYSA-N 0.000 description 1
- DWRXFEITVBNRMK-UHFFFAOYSA-N Beta-D-1-Arabinofuranosylthymine Natural products O=C1NC(=O)C(C)=CN1C1C(O)C(O)C(CO)O1 DWRXFEITVBNRMK-UHFFFAOYSA-N 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 229960000643 adenine Drugs 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- IQFYYKKMVGJFEH-UHFFFAOYSA-N beta-L-thymidine Natural products O=C1NC(=O)C(C)=CN1C1OC(CO)C(O)C1 IQFYYKKMVGJFEH-UHFFFAOYSA-N 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 229940104230 thymidine Drugs 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- 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/1093—General methods of preparing gene libraries, not provided for in other subgroups
-
- 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/20—Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
-
- 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/1082—Preparation or screening gene libraries by chromosomal integration of polynucleotide sequences, HR-, site-specific-recombination, transposons, viral vectors
-
- 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
-
- C—CHEMISTRY; METALLURGY
- C40—COMBINATORIAL TECHNOLOGY
- C40B—COMBINATORIAL CHEMISTRY; LIBRARIES, e.g. CHEMICAL LIBRARIES
- C40B40/00—Libraries per se, e.g. arrays, mixtures
- C40B40/04—Libraries containing only organic compounds
- C40B40/06—Libraries containing nucleotides or polynucleotides, or derivatives thereof
- C40B40/08—Libraries containing RNA or DNA which encodes proteins, e.g. gene libraries
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- 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
- 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
-
- 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
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Biotechnology (AREA)
- Genetics & Genomics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Molecular Biology (AREA)
- Organic Chemistry (AREA)
- Evolutionary Biology (AREA)
- Medical Informatics (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Wood Science & Technology (AREA)
- Zoology (AREA)
- Data Mining & Analysis (AREA)
- Biomedical Technology (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Biochemistry (AREA)
- Microbiology (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Physiology (AREA)
- Bioethics (AREA)
- Immunology (AREA)
- Public Health (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
Abstract
本公开涉及一种基因变异识别方法、装置和存储介质,其中,该方法包括:获取基因变异候选位点对应的至少一个基因测序读段;根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的序列特征和非序列特征,其中,所述序列特征为与位点的位置相关的特征;基于所述序列特征和所述非序列特征,对所述基因变异候选位点的基因变异进行识别。本公开实施例的可以将基因的序列特征和非序列特征相结合,更加全面地分析基因变异位点的特征,增强基因变异识别的准确性。
Description
技术领域
本公开涉及计算机技术领域,尤其涉及一种基因变异识别方法、装置和存储介质。
背景技术
随着生物技术的发展,通过基因测序技术可以测定人类基因的序列,基因序列的分析可以作为进一步基因研究和改造的基础。目前,基因的二代测序技术相比于一代测序技术而言,极大地提高了基因测序的效率,降低了基因测序的成本,并且保持了基因测序的准确行性。第一代测序技术如果完成一个人类基因组的测序可能需要3年的时间,而使用二代测序技术则可以将时间缩短为仅仅1周。
虽然二代测序技术可以生成更大的原始基因测序数据,但是同时还会产生更多的噪声和错误。如何从海量的基因测序数据中识别体细胞基因的变异,筛掉胚系基因变异以及由于噪声和错误带来的干扰,对于二代测序技术的应用具有重要意义。
发明内容
有鉴于此,本公开提出了一种基因变异识别方案。
根据本公开的一方面,提供了一种基因变异识别方法,所述方法包括:
获取基因变异候选位点对应的至少一个基因测序读段;
根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的序列特征和非序列特征,其中,所述序列特征为与位点的位置相关的特征;
基于所述序列特征和所述非序列特征,对所述基因变异候选位点的基因变异进行识别。
在一种可能的实现方式中,所述属性信息包括序列属性信息;根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的序列特征,包括:
根据所述基因变异候选位点的基因位置信息,确定所述基因变异候选位点所在的预设位点区间;
获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息;其中,所述序列属性信息为与位点的位置相关的表征基因属性的信息;
根据所述预设位点区间中每个位点的序列属性信息,生成所述基因变异候选位点的序列特征。
在一种可能的实现方式中,所述获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息,包括:
确定所述至少一个基因测序读段在所述每个位点的基因类型;
统计所述每个位点对应的每种基因类型的基因数量。
在一种可能的实现方式中,所述获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息,包括:
根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,确定每个基因测序读段在所述每个位点的缺失基因的基因类型;
统计所述至少一个基因测序读段在所述每个位点上每种基因类型的缺失基因数量。
在一种可能的实现方式中,所述获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息,包括:
根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,确定每个基因测序读段在所述每个位点的插入基因的基因类型;
统计所述至少一个基因测序读段在所述每个位点上每种基因类型的插入基因数量。
在一种可能的实现方式中,所述序列属性信息包括以下至少一种信息:
参考基因的基因类型;每种基因类型的基因数量;每种基因类型的缺失基因数量;每种基因类型的插入基因数量。
在一种可能的实现方式中,所述属性信息包括非序列属性信息;根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的非序列特征,包括:
获取所述至少一个基因测序读段的非序列属性信息;其中,所述非序列属性信息为与位点的位置不相关的表征基因属性的信息;
根据所述至少一个基因测序读段的非序列属性信息,确定所述基因变异候选位点的非序列特征。
在一种可能的实现方式中,所述非序列信息包括以下至少一种信息:
对比质量;正负链偏好;基因测序读段长度;边缘偏好。
在一种可能的实现方式中,所述根据所述至少一个基因测序读段的非序列属性信息,确定所述基因变异候选位点的非序列特征,包括:
根据每个基因测序读段中每个位点的对比质量,确定每个基因测序读段的对比质量;其中,所述对比质量用于表征基因测序读段中每个基因序列的基因测序的准确性;
根据每个基因测序读段的对比质量,确定所述基因变异候选位点对应的非序列特征。
在一种可能的实现方式中,所述根据所述至少一个基因测序读段的非序列属性信息,确定所述基因变异候选位点的非序列特征,包括:
根据每个基因测序读段所属基因链的正负链信息,确定所述至少一个基因测序读段所属基因链的正负链比例;
根据所述正负链比例,确定所述基因变异候选位点对应的非序列特征。
在一种可能的实现方式中,所述基于所述序列特征和所述非序列特征,对所述基因变异候选位点的基因变异进行识别,包括:
将所述序列特征和所述非序列特征进行特征整合,得到所述基因变异候选位点的整合特征;
基于所述基因变异候选位点的整合特征,对所述基因变异候选位点的基因变异进行识别。
在一种可能的实现方式中,所述基于所述基因变异候选位点的整合特征,对所述基因变异候选位点的基因变异进行识别,包括:
根据所述基因变异候选位点的整合特征,得到所述基因变异候选位点的基因发生变异的变异值;
在所述变异值大于或等于预设阈值的情况下,确定所述基因变异候选位点的基因存在变异。
在一种可能的实现方式中,所述获取基因变异候选位点对应的至少一个基因测序读段,包括:
获取由体细胞基因进行基因测序得到的基因测序读段;
将所述基因测序读段的基因序列与参考基因组的基因序列进行比对,得到比对结果;
根据所述比对结果确定所述体细胞基因的基因存在异常的基因变异候选位点;
获取所述基因变异候选位点对应的至少一个基因测序读段。
根据本公开的另一方面,提供了一种基因变异识别装置,所述装置包括:
获取模块,用于获取基因变异候选位点对应的至少一个基因测序读段;
确定模块,用于根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的序列特征和非序列特征,其中,所述序列特征为与位点的位置相关的特征;
识别模块,用于基于所述序列特征和所述非序列特征,对所述基因变异候选位点的基因变异进行识别。
在一种可能的实现方式中,所述属性信息包括序列属性信息;所述确定模块,包括:
第一确定子模块,用于根据所述基因变异候选位点的基因位置信息,确定所述基因变异候选位点所在的预设位点区间;
第一获取子模块,用于获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息;其中,所述序列属性信息为与位点的位置相关的表征基因属性的信息;
第一生成子模块,用于根据所述预设位点区间中每个位点的序列属性信息,生成所述基因变异候选位点的序列特征。
在一种可能的实现方式中,所述第一获取子模块,具体用于确定所述至少一个基因测序读段在所述每个位点的基因类型;统计所述每个位点对应的每种基因类型的基因数量。
在一种可能的实现方式中,所述第一获取子模块,具体用于根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,确定每个基因测序读段在所述每个位点的缺失基因的基因类型;统计所述至少一个基因测序读段在所述每个位点上每种基因类型的缺失基因数量。
在一种可能的实现方式中,所述第一获取子模块,具体用于根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,确定每个基因测序读段在所述每个位点的插入基因的基因类型;统计所述至少一个基因测序读段在所述每个位点上每种基因类型的插入基因数量。
在一种可能的实现方式中,所述序列属性信息包括以下至少一种信息:
参考基因的基因类型;每种基因类型的基因数量;每种基因类型的缺失基因数量;每种基因类型的插入基因数量。
在一种可能的实现方式中,所述属性信息包括非序列属性信息;所述确定模块,包括:
第二获取子模块,用于获取所述至少一个基因测序读段的非序列属性信息;其中,所述非序列属性信息为与位点的位置不相关的表征基因属性的信息;
第二确定子模块,用于根据所述至少一个基因测序读段的非序列属性信息,确定所述基因变异候选位点的非序列特征。
在一种可能的实现方式中,所述非序列信息包括以下至少一种信息:
对比质量;正负链偏好;基因测序读段长度;边缘偏好。
在一种可能的实现方式中,所述第二确定子模块,具体用于根据每个基因测序读段中每个位点的对比质量,确定每个基因测序读段的对比质量;其中,所述对比质量用于表征基因测序读段中每个基因序列的基因测序的准确性;根据每个基因测序读段的对比质量,确定所述基因变异候选位点对应的非序列特征。
在一种可能的实现方式中,所述第二确定子模块,具体用于根据每个基因测序读段所属基因链的正负链信息,确定所述至少一个基因测序读段所属基因链的正负链比例;根据所述正负链比例,确定所述基因变异候选位点对应的非序列特征。
在一种可能的实现方式中,所述识别模块,包括:
整合子模块,具体用于将所述序列特征和所述非序列特征进行特征整合,得到所述基因变异候选位点的整合特征;
识别子模块,用于基于所述基因变异候选位点的整合特征,对所述基因变异候选位点的基因变异进行识别。
在一种可能的实现方式中,所述识别子模块,具体用于根据所述基因变异候选位点的整合特征,得到所述基因变异候选位点的基因发生变异的变异值;在所述变异值大于或等于预设阈值的情况下,确定所述基因变异候选位点的基因存在变异。
在一种可能的实现方式中,所述获取模块,具体用于,
获取由体细胞基因进行基因测序得到的基因测序读段;
将所述基因测序读段的基因序列与参考基因组的基因序列进行比对,得到比对结果;
根据所述比对结果确定所述体细胞基因的基因存在异常的基因变异候选位点;
获取所述基因变异候选位点对应的至少一个基因测序读段。
根据本公开的另一方面,提供了一种基因变异识别装置,包括:处理器;用于存储处理器可执行指令的存储器;其中,所述处理器被配置为执行上述方法。
根据本公开的另一方面,提供了一种非易失性计算机可读存储介质,其上存储有计算机程序指令,其中,所述计算机程序指令被处理器执行时实现上述方法。
本公开实施例提供获取基因变异候选位点对应的至少一个基因测序读段,可以根据至少一个基因测序读段的属性信息,确定基因变异候选位点的序列特征和非序列特征,从而可以基于确定的序列特征和非序列特征对基因变异候选位点的基因变异进行识别。这里,序列特征可以是与位点的位置相关的特征,非序列特征可以是与位点的位置不相关的特征,从而在基因变异识别过程中,可以将基因的序列特征和非序列特征相结合,更加全面地分析基因变异位点的特征,筛掉胚系基因变异以及由于噪声和错误带来的干扰,更好地对基因变异进行识别,增强基因变异识别的准确性。
根据下面参考附图对示例性实施例的详细说明,本公开的其它特征及方面将变得清楚。
附图说明
包含在说明书中并且构成说明书的一部分的附图与说明书一起示出了本公开的示例性实施例、特征和方面,并且用于解释本公开的原理。
图1示出根据本公开一实施例的基因变异识别方法的流程图。
图2示出根据本公开一实施例的获取基因变异候选位点对应的至少一个基因测序读段的流程图。
图3示出根据本公开一实施例的基因变异候选位点的序列特征过程的流程图。
图4示出根据本公开一实施例的基因变异候选位点的非序列特征过程的流程图。
图5示出根据本公开一实施例的识别基因变异候选位点的基因变异过程的流程图。
图6示出根据本公开一实施例的神经网络模型的框图。
图7示出根据本公开一实施例的基因变异识别装置的框图。
图8示出根据本公开一示例性实施例示出的一种用于基因变异识别的装置的框图。
具体实施方式
以下将参考附图详细说明本公开的各种示例性实施例、特征和方面。附图中相同的附图标记表示功能相同或相似的元件。尽管在附图中示出了实施例的各种方面,但是除非特别指出,不必按比例绘制附图。
在这里专用的词“示例性”意为“用作例子、实施例或说明性”。这里作为“示例性”所说明的任何实施例不必解释为优于或好于其它实施例。
另外,为了更好的说明本公开,在下文的具体实施方式中给出了众多的具体细节。本领域技术人员应当理解,没有某些具体细节,本公开同样可以实施。在一些实例中,对于本领域技术人员熟知的方法、手段、元件和电路未作详细描述,以便于凸显本公开的主旨。
本公开实施例提供的基因变异识别方案,可以获取基因变异候选位点对应的至少一个基因测序读段,从而可以根据至少一个基因测序读段对基因变异候选位点的基因变异进行识别。在基因变异识别过程中,可以根据至少一个基因测序读段的序列属性信息生成序列特征,根据至少一个基因测序读段的非序列属性信息生成非序列特征,然后可以通过序列特征和非序列特征对基因变异候选位点的基因变异进行识别,从而可以整合至少一个基因测序读段的序列属性信息和非序列属性信息,更加全面地利用基因测序读段的序列属性信息。
在相关技术中,通常是利用支持向量机、随机森林等传统随机森林等传统机器学习方法进行基因变异识别,这种方式虽然实现简单,但难以利用基因变异候选位点附近基因序列的序列属性信息,基因变异识别的效果在基因数据量增加到一定程度之后会陷入瓶颈。还有一些相关技术采用深度学习方法,利用神经网络对基因变异进行识别。但是,神经网络难以整合基因序列的非序列信息,无法对基因数据进行更加全面地分析。在本公开实施例中,在基因变异识别过程中可以利用由多模态信息整合的神经网络模型提取基因变异候选位点的序列特征和非序列特征,从而可以综合基因序列的序列属性信息和非序列属性信息,更加全面地对基因数据进行分析,筛掉胚系基因变异以及由于噪声和错误带来的干扰,更好地对基因变异进行识别。下述实施例将会对基因变异识别过程作详细说明。
图1示出根据本公开一实施例的基因变异识别方法的流程图。该基因变异识别方法可以由基因变异识别装置或其它处理设备执行,其中,基因变异识别装置可以为用户设备(User Equipment,UE)、移动设备、用户终端、终端、蜂窝电话、无绳电话、个人数字处理(Personal Digital Assistant,PDA)、手持设备、计算设备、车载设备、可穿戴设备等,或者,基因变异识别装置可以为服务器。在一些可能的实现方式中,该基因变异识别方法可以通过处理器调用存储器中存储的计算机可读指令的方式来实现。
如图1所示,该基因变异识别方法包括:
步骤11,获取基因变异候选位点对应的至少一个基因测序读段。
在本公开实施例中,基因变异识别装置可以获取由基因测序得到的基因测序读段,然后在基因测序得到的基因测序读段中,获取基因变异候选位点对应的至少一个基因测序读段。这里的基因测序读段可以理解为经过基因测序后标注有基因类型的基因序列,每个基因测序读段的长度可以相同也可以不同。在长度不同的情况下,每个基因测序读段的长度可以在预设长度范围内,从而可以保证每个基因测序读段的长度比较接近。基因类型可以理解为碱基类型,基因类型可以包括胞嘧啶(C)、鸟嘌呤(G)、腺嘌呤(A)、胸腺嘧啶(T),从而基因测序读段可以是包括AGCT的基因序列。这里的基因变异候选位点可以是基因序列存在异常的位点。基因序列的位点可以表示基因序列的位置,针对每个位点,可以存在至少一个基因测序读段,即,在同一个位点可以存在由基因测序得到的至少一个基因测序读段。相应地,基因变异候选位点对应至少一个基因测序读段,其中,这至少一个基因测序读段都在这一位点上出现异常。基因变异候选位点可以为至少一个,每个基因变异候选位点可以对应至少一个基因测序读段。为了便于理解,本公开实施例以一个基因变异候选位点进行说明。
步骤12,根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的序列特征和非序列特征,其中,所述序列特征为与位点的位置相关的特征。
在本公开实施例中,在获取基因变异候选位点对应的至少一个基因测序读段之后,可以提取该基因变异候选位点对应的至少一个基因测序读段的属性信息,并根据提取的属性信息生成该基因变异候选位点的序列特征和非序列特征。属性信息可以包括序列属性信息和非序列属性信息。序列属性信息可以是与位点的位置相关的表征基因测序读段的基因属性的信息。非序列属性信息可以是不受到位点的位置限制并且可以表征基因属性的信息。在提取属性信息时,可以随机选择该基因候选位点对应的若干个基因测序读段,提取随机选择的若干个基因测序读段的属性信息;还可以提取该基因候选位点对应的每个基因测序读段的属性信息。
这里,在提取序列属性信息时,可以提取至少一个基因测序读段在该基因变异候选位点的序列属性信息,还可以提取至少一个基因测序读段在该基因变异候选位点附近位点的序列属性信息。这里,在确定基因变异候选位点的序列特征时,可以利用带有卷积层和池化层的神经网络模型,对基因变异候选位点对应的至少一个基因测序读段提取基因变异候选位点的序列特征。该神经网络模型可以包括两个分支结构,其中一个分支可以提取基因测序读段的序列特征,该分支可以包括卷积层和池化层;另一个分支可以提取基因测序读段的非序列特征。该神经网络模型从而可以整合多种模态信息(序列属性信息和非序列属性信息),对基因变异候选位点的基因变异进行识别。在确定基因变异候选位点的非序列特征时,可以利用上述神经网络模型,由该神经网络模型的另一个分支提取至少一个基因测序读段的非序列特征,该分支结构可以包括全连接层,全连接层可以用于提取不受位置限制的非序列特征。
步骤13,基于所述序列特征和所述非序列特征,对所述基因变异候选位点的基因变异进行识别。
在本公开实施方式中,在确定基因变异候选位点的序列特征和非序列特征之后,可将序列特征和非序列特征进行融合,对该基因变异候选位点的基因变异进行识别,例如,可以利用上述神经网络模型判断该基因变异候选位点的基因是否变异,或者,该基因变异候选位点的基因是否是由于噪声等原因而导致的基因序列异常。
本公开实施例中可以根据基因变异候选位点的序列特征和非序列特征对基因变异候选位点的基因变异进行识别,从而可以更加全面地对基因测序数据进行分析。在对基因变异候选位点的基因变异进行识别时,首先需要获取基因变异候选位点对应的至少一个基因测序读段。本公开实例还提供了一种获取基因变异候选位点对应的至少一个基因测序读段的过程。
图2示出根据本公开一实施例的获取基因变异候选位点对应的至少一个基因测序读段的流程图。在一种可能的实现方式中,获取基因变异候选位点对应的至少一个基因测序读段,可以包括以下步骤:
步骤111,获取由体细胞基因进行基因测序得到的基因测序读段。
这里,通过体细胞基因进行基因测序可以得到至少一个基因测序读段,基因测序读段可以是对体细胞基因进行基因类型标注的序列。体细胞基因在进行基因测序之后,不仅可以得到基因测序读段中每个基因的基因类型,还可以得到基因测序读段中每个基因所在位点的基因位置信息。同一个位点可以对应至少一个基因测序读段。
在一种可能的实现方式中,通过体细胞基因进行基因测序可以得到至少一个基因测序读段,可以对基因测序得到的基因测序读段进行预处理,这里的预处理方式可以包括交叉污染筛选、测序质量筛选、比对质量筛选、读段长度异常筛选等。通过预处理,可以筛选掉交叉污染的基因测序读段,以及筛选掉测序质量和比对质量较低、读段长度异常的基因测序读段。
步骤112,将所述基因测序读段的基因序列与参考基因组的基因序列进行比对,得到比对结果。
在本公开实施例中,在获取由体细胞基因进行基因测序得到的基因测序读段之后,可以将获取的基因测序读段的基因序列与相同位点的参考基因组的基因序列的进行比对,得到对比结果。举例来说,可以将每个进行基因测序得到的基因测序读段与相同位点的参考基因组的基因序列进行对比,确定基因测序读段的基因序列与参考基因组的基因序列不同的位点。还可以将具有相同位点的至少一个基因测序读段与相同位点的参考基因组的基因序列进行对比,确定至少一个基因测序读段的基因序列与参考基因组的基因序列不同的位点。
步骤113,根据所述比对结果确定所述体细胞基因的基因存在异常的基因变异候选位点。
在本公开实施例中,可以根据比对结果确定基因测序读段与参考基因组的基因序列不同的位点,如果该位点对应的至少一个基因测序读段中,在该位点发送变异的基因测序读段的比例大于预设比例,则可以确定该位点为基因变异候选位点,否则,可以认为该位点不是基因变异候选位点。基因测序读段在该位点与参考基因组的基因序列不同,可能是因为测序错误导致的不同,通过这种方式,可以减少由于基因测序失误引起的基因序列异常现象。
步骤114,获取所述基因变异候选位点对应的至少一个基因测序读段。
在本公开实施例中,在确定基因变异候选位点之后,可以获取基因变异候选位点对应的至少一个基因测序读段。其中,每个基因变异候选位点对应的至少一个基因测序读段,在该基因变异候选位点的基因序列与相同位点的参考基因组的基因序列可以不同。这里的基因变异候选位点可以为至少一个。
通过上述获取基因变异候选位点对应的至少一个基因测序读段的过程,不仅可以较为准确地确定基因变异候选位点,还可以在基因测序得到的基因测序读段中确定基因变异候选位点对应的至少一个基因测序读段。
本公开实施例中可以根据基因变异候选位点对应的至少一个基因测序读段的序列属性信息,确定该基因变异候选位点的序列特征,从而在对基因变异候选位点的基因变异进行识别时,可以考虑基因变异候选位点所对应的至少一个基因测序读段的序列属性。下面通过一示例对确定基因变异候选位点的序列特征的过程进行详细说明。
图3示出根据本公开一实施例的基因变异候选位点的序列特征过程的流程图。如图3所示,上述步骤12可以包括以下步骤:
步骤121a,根据所述基因变异候选位点的基因位置信息,确定所述基因变异候选位点所在的预设位点区间;
步骤122a,获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息;其中,所述序列属性信息为与位点的位置相关的表征基因属性的信息;
步骤123a,根据所述预设位点区间中每个位点的序列属性信息,生成所述基因变异候选位点的序列特征。
在本公开实施例的示例中,对于每一个基因变异候选位点可以存在至少一个基因测序读段。为了提高基因变异识别的准确度,不仅可以考虑该基因变异候选位点的序列属性信息,还可以考虑该基因变异候选位点附近的位点的序列属性信息。在确定基因变异候选位点的序列特征时,可以根据基因变异候选位点的基因位置信息,确定该基因变异候选位点所在的预设位点区间,例如,可以将基因变异候选位点前后150个碱基对的区间作为基因变异候选位点所在的预设位点区间。然后可以针对该预设位点区间内的每个位点,获取至少一个基因测序读段在该位点的序列属性信息,由该位点的序列属性信息可以生成该位点对应序列特征。序列特征可以用序列特征向量进行表示。由基因变异候选位点所在预设位点区间中至少一个位点对应的至少一个序列特征向量,可以形成基因变异候选位点的序列特征矩阵。举例来说,若基因变异候选位点所在预设位点区间包括3个位点b1、b2、b3,3个位点对应的序列特征向量分别为a1、a2、a3,基因变异候选位点的序列特征矩阵为[a1a2a3],其中,a1、a2、a3的序列特征分别对应b1、b2、b3的序列属性信息。
这里,序列属性信息可以包括但不限于:参考基因组的基因类型;每种基因类型的基因数量;每种基因类型的缺失基因数量;每种基因类型的插入基因数量。参考基因组的基因类型可以是参考基因组在基因变异候选位点的基因类型。每种基因类型的基因数量可以是至少一个基因测序读段在该基因变异候选位点每种基因类型的基因数量,例如,该基因变异候选位点对应5个基因测序读段,每个基因测序读段在该基因变异候选位点的基因类型分别为:A、C、C、G、G,则每种基因类型的基因数量分别为:A为1个;C为2个;G为2个。每种基因类型的缺失基因数量可以是至少一个基因测序读段在该基因变异候选位点每种基因类型的缺失基因数量,例如,每个基因测序读段在该基因变异候选位点缺失的基因类型分别为:A、C、C、G、G,则每种基因类型的缺失基因数量分别为:A为1个;C为2个;G为2个。每种基因类型的插入基因数量可以是至少一个基因测序读段在该基因变异候选位点每种基因类型的插入基因数量,例如,每个基因测序读段在该基因变异候选位点插入的基因类型分别为:A、C、C、G、G,则每种基因类型的插入基因数量分别为:A为1个;C为2个;G为2个。
在一种可能的实现方式中,在获取至少一个基因测序读段在预设位点区间中每个位点的序列属性信息时,可以针对该预设位点区间中的每个位点,确定至少一个基因测序读段在该位点的基因类型,并统计该位点所对应的每种基因类型的基因数量,从而可以确定基因变异候选位点对应的至少一个基因测序读段,在该位点每种基因类型的基因数量。
在一种可能的实现方式中,在获取至少一个基因测序读段在预设位点区间中每个位点的序列属性信息时,可以根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,针对该预设位点区间中的每个位点,确定每个基因测序读段在该位点的缺失基因的基因类型,并统计至少一个基因测序读段在该位点上每种基因类型的缺失基因数量,从而可以确定基因变异候选位点对应的至少一个基因测序读段,在该位点每种基因类型的缺失基因数量。
在一种可能的实现方式中,在获取至少一个基因测序读段在预设位点区间中每个位点的序列属性信息时,可以根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,针对该预设位点区间中的每个位点,确定每个基因测序读段在该位点的缺失基因的基因类型,并统计至少一个基因测序读段在该位点上每种基因类型的插入基因数量,从而可以确定基因变异候选位点对应的至少一个基因测序读段,在该位点每种基因类型的插入基因数量。
举例来说,假设序列属性信息包括参考基因组的基因类型、每种基因类型的基因数量、每种基因类型的缺失基因数量、每种基因类型的插入基因数量,在确定基因变异候选位点的序列特征时,可以针对基因变异候选位点所在的预设位点区间中的每一个位点,提取基因变异候选位点对应的至少一个基因测序读段在该位点的上述四个信息,例如,基因变异候选位点对应的5个基因测序读段,针对预预设位点区间中的某一位点,可以分别确定参考基因组在该位点的基因类型、5个基因测序读段在该位点各基因类型的基因数量、5个基因测序读段在该位点各基因类型的缺失基因数量和5个基因测序读段在该位点各基因类型的插入基因数量。然后综合该位点对应的至少一个序列属性信息,可以得到该位点的序列特征。基因变异候选位点的序列特征可以包括预设位点区间中每个位点的序列特征。
本公开实施例的示例中不仅在对基因变异候选位点的基因变异进行识别时,考虑了基因变异候选位点所对应的至少一个基因测序读段的序列属性,还考虑了至少一个基因测序读段的非序列属性。下面通过一示例对确定基因变异候选位点的非序列特征的过程进行详细说明。
图4示出根据本公开一实施例的基因变异候选位点的非序列特征过程的流程图。如图4所示,上述步骤12可以包括以下步骤:
步骤121b,获取所述至少一个基因测序读段的非序列属性信息;其中,所述非序列属性信息为与位点的位置不相关的表征基因属性的信息;
步骤122b,根据所述至少一个基因测序读段的非序列属性信息,生成所述基因变异候选位点的非序列特征。
在本公开实施例的示例中,为了提高基因变异识别的准确度,不仅可以考虑至少一个基因测序读段的序列属性信息,还可以考虑至少一个基因测序读段的非序列属性信息。这里,非序列信息可以包括以下至少一种信息:对比质量;正负链偏好;基因测序读段长度;边缘偏好。在确定基因变异候选位点的非序列特征时,可以获取至少一个基因属性序列读段的非序列属性信息,然后由获取的非序列属性信息生成基因变异候选位点的非序列特征。
在一种可能的实现方式中,在根据所述至少一个基因测序读段的非序列属性信息,确定所述基因变异候选位点的非序列特征时,可以根据每个基因测序读段中每个位点的对比质量,确定每个基因测序读段的对比质量,然后根据每个基因测序读段的对比质量,确定所述基因变异候选位点对应的非序列特征。这里,对比质量可以用于表征基因测序读段中每个基因序列的基因测序的准确性,如果某个基因序列的对比质量低于预设值,则可以认为该基因序列由基因测序得到的基因类型不准确,从而可以将对比质量作为判断基因变异候选位点的基因是否发生变异的一个参考因素。举例来说,基因变异候选位点对应至少一个基因测序读段,则可以根据每个基因序列的对比质量,确定每个基因测序读段的对比质量,以一个基因测序读段举例,可以将该基因测序读段所包括的基因序列的对比质量的平均值或者中间值,作为该基因测序读段的对比质量,还可以在该基因测序读段随机选择至少一个基因序列,将选择的至少一个基因序列对比质量的平均值或者中间值作为该基因测序读段的对比质量。然后由每个基因测序读段的对比质量得到该基因变异候选位点对应的对比质量,例如,计算该基因变异候选位点对应的至少一个基因测序读段对比质量的平均值或者均值,得到该基因变异候选位点对应的对比质量,从而可以根据该基因变异候选位点对应的对比质量确定基因变异候选位点对应的非序列特征。
在一种可能的实现方式中,在根据至少一个基因测序读段的非序列属性信息,确定基因变异候选位点的非序列特征时,可以根据每个基因测序读段所属基因链的正负链信息,确定至少一个基因测序读段所属基因链的正负链比例,然后根据确定的正负链比例,确定基因变异候选位点对应的非序列特征。这里,正负链偏好可以是基因测序读段所属基因链中正链和负链的比例,基因链可以包括正链和负链,其中,正链可以是与核糖核酸(RNA)的碱基序列相同的脱氧核糖核酸(DNA)单链,负链可以是与核糖核酸(RNA)的碱基序列互补的脱氧核糖核酸(DNA)单链。举例来说,基因变异候选位点对应5个基因测序读段,其中,3个基因测序读段对应基因链的正链,2个基因测序读段对应基因链的负链,则正负链偏好可以是3:2。
在一种可能的实现方式中,在根据至少一个基因测序读段的非序列属性信息,确定基因变异候选位点的非序列特征时,可以根据每个基因测序读段的基因测序读段长度,确定基因变异候选位点的非序列特征。基因测序读段长度可以是每个基因测序读段所具有碱基序列的长度,举例来说,一个基因测序读段包括4个碱基序列,则该基因测序读段的长度为4,可以由每个基因测序读段长度确定基因变异候选位点的非序列特征,还可以由至少一个基因测序读段长度的中间值或者平均值确定基因变异候选位点的非序列特征。
在一种可能的实现方式中,在根据至少一个基因测序读段的非序列属性信息,确定基因变异候选位点的非序列特征时,可以根据每个基因测序读段的边缘偏好,确定基因变异候选位点的非序列特征。这里,边缘偏好可以是某一位点在基因测序读段中位于边缘位置与中间位置的比例。举例来说,可以将基因测序读段平均分为3段,其中,基因测序读段两端的2段可以作为边缘位置,基因测序读段中间的1段可以作为中间位置,基因变异候选位点对应5个基因测序读段,基因变异候选位点如果位于其中3个基因测序读段的边缘位置,位于其中2个基因测序读段的中间位置,该基因变异候选位点的边缘偏好可以为3:2。相应地,可以由基因变异候选位点在每个基因测序读段的边缘偏好,确定基因变异候选位点的非序列特征,还可以由至少一个基因测序读段所对应的边缘偏好的中间值或者平均值,确定基因变异候选位点的非序列特征。
通过上述方式,可以针对至少一个基因测序读段在基因变异候选位点的非序列属性信息生成基因变异候选位点的非序列特征,从而可以在基因变异识别时考虑基因变异候选位点的非序列特征度特征,使基因变异识别更加准确。在确定非序列特征时,可以是由非序列属性信息中任意至少一个信息的组合生成至少一个基因测序读段的非序列特征。
下面通过一示例对基因变异候选位点的基因变异进行识别的过程进行说明。
图5示出根据本公开一实施例的识别基因变异候选位点的基因变异过程的流程图。如图5所示,上述步骤13可以包括以下步骤:
步骤131,将所述序列特征和所述非序列特征进行特征整合,得到所述基因变异候选位点的整合特征;
步骤132,基于所述基因变异候选位点的整合特征,对所述基因变异候选位点的基因变异进行识别。
在本公开实施例中,在确定基因变异候选位点的序列特征和非序列维度特征之后,可以利用神经网络模型对序列特征和非序列特征进行特征整合,将序列特征形成的序列特征矩阵与非序列特征形成的非序列特征矩阵合成为一个特征矩阵,得到由整合特征形成的整合特征矩阵,然后利用神经网络模型根据该整合特征矩阵对变异候选位点的基因变异进行识别。通过这种方式,可以利用神经网络模型整合基因变异候选位点对应的序列属性信息和非序列属性信息,从而可以更加全面地对基因测序数据进行分析,使基因变异识别更加准确。在训练过程中,可以选取存在单核苷酸多态性(Single NucleotidePolymorphism,SNP)的基因测序读段、存在插入/缺失(Insertion/Deletion,InDel)的基因测序读段作为训练样本,从而训练后得到的基因变异识别模型可以有效地对SNP、InDel的基因变异进行识别。
在一种可能的实现方式中,根据所述基因变异候选位点的整合特征,对所述基因变异候选位点的基因变异进行识别,可以包括:根据所述基因变异候选位点的整合特征,得到所述基因变异候选位点的基因发生变异的变异值;在所述变异值大于或等于预设阈值的情况下,确定所述基因变异候选位点的基因存在变异。这里,基因发生变异的变异值可以是表征该基因变异候选位点发生变异的可能性,例如,变异值越大,该基因变异候选位点发生变异的可能性越大。可以利用上述神经网络对二维特征进行处理得到变异值,并根据变异值判断基因变异候选位点的基因是否存在变异。在一种可能的实现方式中,变异值可以在0至1之间。预设阈值可以根据应用场景进行设置,例如,0.3、0.5,如果变异值大于预设阈值,则可以认为该基因变异候选位点的基因发生变异,否则,可以为该基因变异候选位点的基因未发生变异。
本公开实施例中可以利用神经网络模型对基因变异候选位点的基因变异进行识别,该神经网络模型可以提取基因变异候选位点的序列特征和非序列特征。本公开实施例还提供了一种神经网络模型的结构。
图6示出根据本公开一实施例的神经网络模型的框图。如图6所示,神经网络模型可以包括两个分支结构,第一分支和第二分支。第一分支可以用于提取基因变异候选位点对应的至少一个基因测序读段的序列特征,第一分支可以包括卷积层和池化层。第二分支可以用于提取基因变异候选位点对应的至少一个基因测序读段的非序列特征,第二分支可以包括全连接层。神经网络模型提取基因变异候选位点的序列特征和非序列特征之后,可以将序列特征和非序列特征进行整合,例如,将序列特征的序列特征矩阵与非序列特征的非序列特征矩阵进行拼接,得到整合特征的整合特征矩阵,然后再经过全连接层可以得到基因变异候选位点的变异值。
本公开实施例通过提取基因变异候选位点对应的至少一个基因测序读段的序列属性信息和非序列属性信息,利用对序列属性信息和非序列属性信息整合的整合特征对基因变异进行识别,从而综合考虑基因变异候选位点对应的序列属性信息和非序列属性信息,更加全面地分析基因测序信息,更好地对基因候选位点的基因变异进行识别,筛掉胚系基因变异以及由于噪声和错误带来的干扰,提高基因变异识别的准确率。
本领域技术人员可以理解,在具体实施方式的上述方法中,各步骤的撰写顺序并不意味着严格的执行顺序而对实施过程构成任何限定,各步骤的具体执行顺序应当以其功能和可能的内在逻辑确定。
图7示出根据本公开实施例的基因变异识别装置的框图,如图7所示,所述基因变异识别装置包括:
获取模块71,用于获取基因变异候选位点对应的至少一个基因测序读段;
确定模块72,用于根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的序列特征和非序列特征,其中,所述序列特征为与位点的位置相关的特征;
识别模块73,用于基于所述序列特征和所述非序列特征,对所述基因变异候选位点的基因变异进行识别。
在一种可能的实现方式中,所述属性信息包括序列属性信息;所述确定模块72,包括:
第一确定子模块,用于根据所述基因变异候选位点的基因位置信息,确定所述基因变异候选位点所在的预设位点区间;
第一获取子模块,用于获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息;其中,所述序列属性信息为与位点的位置相关的表征基因属性的信息;
第一生成子模块,用于根据所述预设位点区间中每个位点的序列属性信息,生成所述基因变异候选位点的序列特征。
在一种可能的实现方式中,所述第一获取子模块,具体用于确定所述至少一个基因测序读段在所述每个位点的基因类型;统计所述每个位点对应的每种基因类型的基因数量。
在一种可能的实现方式中,所述第一获取子模块,具体用于根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,确定每个基因测序读段在所述每个位点的缺失基因的基因类型;统计所述至少一个基因测序读段在所述每个位点上每种基因类型的缺失基因数量。
在一种可能的实现方式中,所述第一获取子模块,具体用于根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,确定每个基因测序读段在所述每个位点的插入基因的基因类型;统计所述至少一个基因测序读段在所述每个位点上每种基因类型的插入基因数量。
在一种可能的实现方式中,所述序列属性信息包括以下至少一种信息:
参考基因的基因类型;每种基因类型的基因数量;每种基因类型的缺失基因数量;每种基因类型的插入基因数量。
在一种可能的实现方式中,所述属性信息包括非序列属性信息;所述确定模块,包括:
第二获取子模块,用于获取所述至少一个基因测序读段的非序列属性信息;其中,所述非序列属性信息为与位点的位置不相关的表征基因属性的信息;
第二确定子模块,用于根据所述至少一个基因测序读段的非序列属性信息,确定所述基因变异候选位点的非序列特征。
在一种可能的实现方式中,所述非序列信息包括以下至少一种信息:
对比质量;正负链偏好;基因测序读段长度;边缘偏好。
在一种可能的实现方式中,所述第二确定子模块,具体用于根据每个基因测序读段中每个位点的对比质量,确定每个基因测序读段的对比质量;其中,所述对比质量用于表征基因测序读段中每个基因序列的基因测序的准确性;根据每个基因测序读段的对比质量,确定所述基因变异候选位点对应的非序列特征。
在一种可能的实现方式中,所述第二确定子模块,具体用于根据每个基因测序读段所属基因链的正负链信息,确定所述至少一个基因测序读段所属基因链的正负链比例;根据所述正负链比例,确定所述基因变异候选位点对应的非序列特征。
在一种可能的实现方式中,所述识别模块73,包括:
整合子模块,具体用于将所述序列特征和所述非序列特征进行特征整合,得到所述基因变异候选位点的整合特征;
识别子模块,用于基于所述基因变异候选位点的整合特征,对所述基因变异候选位点的基因变异进行识别。
在一种可能的实现方式中,所述识别子模块,具体用于根据所述基因变异候选位点的整合特征,得到所述基因变异候选位点的基因发生变异的变异值;在所述变异值大于或等于预设阈值的情况下,确定所述基因变异候选位点的基因存在变异。
在一种可能的实现方式中,所述获取模块71,具体用于,
获取由体细胞基因进行基因测序得到的基因测序读段;
将所述基因测序读段的基因序列与参考基因组的基因序列进行比对,得到比对结果;
根据所述比对结果确定所述体细胞基因的基因存在异常的基因变异候选位点;
获取所述基因变异候选位点对应的至少一个基因测序读段。
在一些实施例中,本公开实施例提供的装置具有的功能或包含的模块可以用于执行上文方法实施例描述的方法,其具体实现可以参照上文方法实施例的描述,为了简洁,这里不再赘述。
图8是根据一示例性实施例示出的一种用于基因变异识别的装置1900的框图。例如,装置1900可以被提供为一服务器。参照图8,装置1900包括处理组件1922,其进一步包括一个或多个处理器,以及由存储器1932所代表的存储器资源,用于存储可由处理组件1922的执行的指令,例如应用程序。存储器1932中存储的应用程序可以包括一个或一个以上的每一个对应于一组指令的模块。此外,处理组件1922被配置为执行指令,以执行上述方法。
装置1900还可以包括一个电源组件1926被配置为执行装置1900的电源管理,一个有线或无线网络接口1950被配置为将装置1900连接到网络,和一个输入输出(I/O)接口1958。装置1900可以操作基于存储在存储器1932的操作系统,例如Windows ServerTM,MacOS XTM,UnixTM,LinuxTM,FreeBSDTM或类似。
在示例性实施例中,还提供了一种非易失性计算机可读存储介质,例如包括计算机程序指令的存储器1932,上述计算机程序指令可由装置1900的处理组件1922执行以完成上述方法。
本公开可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质,其上载有用于使处理器实现本公开的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身,诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如,通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备,或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令,并转发该计算机可读程序指令,以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码,所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等,以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中,通过利用计算机可读程序指令的状态信息来个性化定制电子电路,例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA),该电子电路可以执行计算机可读程序指令,从而实现本公开的各个方面。
这里参照根据本公开实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本公开的各个方面。应当理解,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机可读程序指令实现。
这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器,从而生产出一种机器,使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时,产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中,这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作,从而,存储有指令的计算机可读介质则包括一个制造品,其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。
也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上,使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤,以产生计算机实现的过程,从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。
附图中的流程图和框图显示了根据本公开的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分,所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中技术的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (10)
1.一种基因变异识别方法,其特征在于,所述方法包括:
获取基因变异候选位点对应的至少一个基因测序读段;
根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的序列特征和非序列特征,其中,所述序列特征为与位点的位置相关的特征;
基于所述序列特征和所述非序列特征,对所述基因变异候选位点的基因变异进行识别。
2.根据权利要求1所述的方法,其特征在于,所述属性信息包括序列属性信息;根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的序列特征,包括:
根据所述基因变异候选位点的基因位置信息,确定所述基因变异候选位点所在的预设位点区间;
获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息;其中,所述序列属性信息为与位点的位置相关的表征基因属性的信息;
根据所述预设位点区间中每个位点的序列属性信息,生成所述基因变异候选位点的序列特征。
3.根据权利要求2所述的方法,其特征在于,所述获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息,包括:
确定所述至少一个基因测序读段在所述每个位点的基因类型;
统计所述每个位点对应的每种基因类型的基因数量。
4.根据权利要求2所述的方法,其特征在于,所述获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息,包括:
根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,确定每个基因测序读段在所述每个位点的缺失基因的基因类型;
统计所述至少一个基因测序读段在所述每个位点上每种基因类型的缺失基因数量。
5.根据权利要求2所述的方法,其特征在于,所述获取所述至少一个基因测序读段在所述预设位点区间中每个位点的序列属性信息,包括:
根据每个基因测序读段的基因序列与参考基因组的基因序列进行比对的比对结果,确定每个基因测序读段在所述每个位点的插入基因的基因类型;
统计所述至少一个基因测序读段在所述每个位点上每种基因类型的插入基因数量。
6.根据权利要求1所述的方法,其特征在于,所述序列属性信息包括以下至少一种信息:
参考基因的基因类型;每种基因类型的基因数量;每种基因类型的缺失基因数量;每种基因类型的插入基因数量。
7.根据权利要求1至6任意一项所述的方法,其特征在于,所述属性信息包括非序列属性信息;根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的非序列特征,包括:
获取所述至少一个基因测序读段的非序列属性信息;其中,所述非序列属性信息为与位点的位置不相关的表征基因属性的信息;
根据所述至少一个基因测序读段的非序列属性信息,确定所述基因变异候选位点的非序列特征。
8.一种基因变异识别装置,其特征在于,所述装置包括:
获取模块,用于获取基因变异候选位点对应的至少一个基因测序读段;
确定模块,用于根据所述至少一个基因测序读段的属性信息,确定所述基因变异候选位点的序列特征和非序列特征,其中,所述序列特征为与位点的位置相关的特征;
识别模块,用于基于所述序列特征和所述非序列特征,对所述基因变异候选位点的基因变异进行识别。
9.一种基因变异识别装置,其特征在于,包括:
处理器;
用于存储处理器可执行指令的存储器;
其中,其中,所述处理器通过调用所述可执行指令实现如权利要求1至7中任意一项所述的方法。
10.一种非易失性计算机可读存储介质,其上存储有计算机程序指令,其特征在于,所述计算机程序指令被处理器执行时实现权利要求1至7中任意一项所述的方法。
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910251891.0A CN109994155B (zh) | 2019-03-29 | 2019-03-29 | 一种基因变异识别方法、装置和存储介质 |
JP2021514554A JP7064654B2 (ja) | 2019-03-29 | 2019-05-31 | 遺伝子変異認識方法、装置および記憶媒体 |
KR1020217020204A KR20210116454A (ko) | 2019-03-29 | 2019-05-31 | 유전자 변이 인식 방법 및 장치 및 기억 매체 |
SG11202011523VA SG11202011523VA (en) | 2019-03-29 | 2019-05-31 | Gene mutation identification method and apparatus, and storage medium |
PCT/CN2019/089499 WO2020199336A1 (zh) | 2019-03-29 | 2019-05-31 | 一种基因变异识别方法、装置和存储介质 |
TW108137265A TWI748263B (zh) | 2019-03-29 | 2019-10-16 | 一種基因變異辨識方法、裝置和儲存介質 |
US17/102,136 US20210082539A1 (en) | 2019-03-29 | 2020-11-23 | Gene mutation identification method and apparatus, and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910251891.0A CN109994155B (zh) | 2019-03-29 | 2019-03-29 | 一种基因变异识别方法、装置和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109994155A true CN109994155A (zh) | 2019-07-09 |
CN109994155B CN109994155B (zh) | 2021-08-20 |
Family
ID=67131990
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910251891.0A Active CN109994155B (zh) | 2019-03-29 | 2019-03-29 | 一种基因变异识别方法、装置和存储介质 |
Country Status (7)
Country | Link |
---|---|
US (1) | US20210082539A1 (zh) |
JP (1) | JP7064654B2 (zh) |
KR (1) | KR20210116454A (zh) |
CN (1) | CN109994155B (zh) |
SG (1) | SG11202011523VA (zh) |
TW (1) | TWI748263B (zh) |
WO (1) | WO2020199336A1 (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111081318A (zh) * | 2019-12-06 | 2020-04-28 | 人和未来生物科技(长沙)有限公司 | 一种融合基因检测方法、系统和介质 |
CN111081314A (zh) * | 2019-12-13 | 2020-04-28 | 北京市商汤科技开发有限公司 | 基因变异的识别方法及装置、电子设备和存储介质 |
CN111081313A (zh) * | 2019-12-13 | 2020-04-28 | 北京市商汤科技开发有限公司 | 基因变异的识别方法及装置、电子设备和存储介质 |
CN111091867A (zh) * | 2019-12-18 | 2020-05-01 | 中国科学院大学 | 基因变异位点筛选方法及系统 |
CN111091873A (zh) * | 2019-12-13 | 2020-05-01 | 北京市商汤科技开发有限公司 | 基因变异的识别方法及装置、电子设备和存储介质 |
CN111304308A (zh) * | 2020-03-02 | 2020-06-19 | 北京泛生子基因科技有限公司 | 一种审核高通量测序基因变异检测结果的方法 |
CN113539357A (zh) * | 2021-06-10 | 2021-10-22 | 阿里巴巴新加坡控股有限公司 | 基因检测方法、模型训练方法、装置、设备及系统 |
CN115620802A (zh) * | 2022-09-02 | 2023-01-17 | 蔓之研(上海)生物科技有限公司 | 一种基因数据的处理方法及系统 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113517022A (zh) * | 2021-06-10 | 2021-10-19 | 阿里巴巴新加坡控股有限公司 | 基因检测方法、特征提取方法、装置、设备及系统 |
CN113299344A (zh) * | 2021-06-23 | 2021-08-24 | 深圳华大医学检验实验室 | 基因测序分析方法、装置、存储介质和计算机设备 |
CN113628683B (zh) * | 2021-08-24 | 2024-04-09 | 慧算医疗科技(上海)有限公司 | 一种高通量测序突变检测方法、设备、装置及可读存储介质 |
CN115458052B (zh) * | 2022-08-16 | 2023-06-30 | 珠海横琴铂华医学检验有限公司 | 基于一代测序的基因突变分析方法、设备和存储介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014129894A1 (en) * | 2013-02-19 | 2014-08-28 | Cergentis B.V. | Sequencing strategies for genomic regions of interest |
CN104462869A (zh) * | 2014-11-28 | 2015-03-25 | 天津诺禾致源生物信息科技有限公司 | 检测体细胞单核苷酸突变的方法和装置 |
KR20160010277A (ko) * | 2014-07-18 | 2016-01-27 | 에스케이텔레콤 주식회사 | 산모의 무세포 dna의 차세대 서열분석을 통한 태아의 단일유전자 유전변이의 예측방법 |
CN105574361A (zh) * | 2015-11-05 | 2016-05-11 | 上海序康医疗科技有限公司 | 一种检测基因组拷贝数变异的方法 |
CN106529211A (zh) * | 2016-11-04 | 2017-03-22 | 成都鑫云解码科技有限公司 | 变异位点的获取方法及装置 |
CN106683081A (zh) * | 2016-12-17 | 2017-05-17 | 复旦大学 | 基于影像组学的脑胶质瘤分子标记物无损预测方法和预测系统 |
CN108021788A (zh) * | 2017-12-06 | 2018-05-11 | 深圳市新合生物医疗科技有限公司 | 基于细胞游离dna的深度测序数据提取生物标记物的方法和装置 |
US20180280391A1 (en) * | 2011-08-23 | 2018-10-04 | Foundation Medicine, Inc. | Novel kif5b-ret fusion molecules and uses thereof |
CN109326316A (zh) * | 2018-09-18 | 2019-02-12 | 哈尔滨工业大学(深圳) | 一种癌症相关SNP、基因、miRNA和蛋白质相互作用的多层网络模型构建方法和应用 |
WO2019242991A1 (en) * | 2018-06-22 | 2019-12-26 | Julius-Maximilians-Universität Würzburg | Method for determining a quantification of old and new rna |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2528205B (en) * | 2013-03-15 | 2020-06-03 | Guardant Health Inc | Systems and methods to detect rare mutations and copy number variation |
CN104293940B (zh) * | 2014-09-30 | 2017-07-28 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
CN105989246B (zh) * | 2015-01-28 | 2018-10-26 | 深圳华大智造科技有限公司 | 一种基于基因组组装的变异检测方法和装置 |
AU2016258914A1 (en) * | 2015-05-01 | 2017-11-09 | Guardant Health, Inc. | Diagnostic methods |
JP6675164B2 (ja) * | 2015-07-28 | 2020-04-01 | 株式会社理研ジェネシス | 変異判定方法、変異判定プログラムおよび記録媒体 |
JP6679065B2 (ja) * | 2015-10-07 | 2020-04-15 | 国立研究開発法人国立がん研究センター | 稀少突然変異の検出方法、検出装置及びコンピュータプログラム |
KR101936933B1 (ko) * | 2016-11-29 | 2019-01-09 | 연세대학교 산학협력단 | 염기서열의 변이 검출방법 및 이를 이용한 염기서열의 변이 검출 디바이스 |
CN106611106B (zh) * | 2016-12-06 | 2019-05-03 | 北京荣之联科技股份有限公司 | 基因变异检测方法及装置 |
KR102035615B1 (ko) * | 2017-08-07 | 2019-10-23 | 연세대학교 산학협력단 | 유전자 패널에 기초한 염기서열의 변이 검출방법 및 이를 이용한 염기서열의 변이 검출 디바이스 |
CN107944228B (zh) * | 2017-12-08 | 2021-06-01 | 广州漫瑞生物信息技术有限公司 | 一种基因测序变异位点的可视化方法 |
-
2019
- 2019-03-29 CN CN201910251891.0A patent/CN109994155B/zh active Active
- 2019-05-31 WO PCT/CN2019/089499 patent/WO2020199336A1/zh active Application Filing
- 2019-05-31 JP JP2021514554A patent/JP7064654B2/ja active Active
- 2019-05-31 SG SG11202011523VA patent/SG11202011523VA/en unknown
- 2019-05-31 KR KR1020217020204A patent/KR20210116454A/ko not_active Application Discontinuation
- 2019-10-16 TW TW108137265A patent/TWI748263B/zh active
-
2020
- 2020-11-23 US US17/102,136 patent/US20210082539A1/en active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180280391A1 (en) * | 2011-08-23 | 2018-10-04 | Foundation Medicine, Inc. | Novel kif5b-ret fusion molecules and uses thereof |
WO2014129894A1 (en) * | 2013-02-19 | 2014-08-28 | Cergentis B.V. | Sequencing strategies for genomic regions of interest |
KR20160010277A (ko) * | 2014-07-18 | 2016-01-27 | 에스케이텔레콤 주식회사 | 산모의 무세포 dna의 차세대 서열분석을 통한 태아의 단일유전자 유전변이의 예측방법 |
CN104462869A (zh) * | 2014-11-28 | 2015-03-25 | 天津诺禾致源生物信息科技有限公司 | 检测体细胞单核苷酸突变的方法和装置 |
CN105574361A (zh) * | 2015-11-05 | 2016-05-11 | 上海序康医疗科技有限公司 | 一种检测基因组拷贝数变异的方法 |
CN106529211A (zh) * | 2016-11-04 | 2017-03-22 | 成都鑫云解码科技有限公司 | 变异位点的获取方法及装置 |
CN106683081A (zh) * | 2016-12-17 | 2017-05-17 | 复旦大学 | 基于影像组学的脑胶质瘤分子标记物无损预测方法和预测系统 |
CN108021788A (zh) * | 2017-12-06 | 2018-05-11 | 深圳市新合生物医疗科技有限公司 | 基于细胞游离dna的深度测序数据提取生物标记物的方法和装置 |
WO2019242991A1 (en) * | 2018-06-22 | 2019-12-26 | Julius-Maximilians-Universität Würzburg | Method for determining a quantification of old and new rna |
CN109326316A (zh) * | 2018-09-18 | 2019-02-12 | 哈尔滨工业大学(深圳) | 一种癌症相关SNP、基因、miRNA和蛋白质相互作用的多层网络模型构建方法和应用 |
Non-Patent Citations (2)
Title |
---|
RYAN POPLIN 等: "Creating a universal SNP and small indel variant caller with deep neural networks", 《NATURE BIOTECHNOLOGY》 * |
赵彦: "基于神经网络的基因缺失变异检测", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111081318A (zh) * | 2019-12-06 | 2020-04-28 | 人和未来生物科技(长沙)有限公司 | 一种融合基因检测方法、系统和介质 |
CN111081314A (zh) * | 2019-12-13 | 2020-04-28 | 北京市商汤科技开发有限公司 | 基因变异的识别方法及装置、电子设备和存储介质 |
CN111081313A (zh) * | 2019-12-13 | 2020-04-28 | 北京市商汤科技开发有限公司 | 基因变异的识别方法及装置、电子设备和存储介质 |
CN111091873A (zh) * | 2019-12-13 | 2020-05-01 | 北京市商汤科技开发有限公司 | 基因变异的识别方法及装置、电子设备和存储介质 |
CN111091873B (zh) * | 2019-12-13 | 2023-07-18 | 北京市商汤科技开发有限公司 | 基因变异的识别方法及装置、电子设备和存储介质 |
CN111091867A (zh) * | 2019-12-18 | 2020-05-01 | 中国科学院大学 | 基因变异位点筛选方法及系统 |
CN111304308A (zh) * | 2020-03-02 | 2020-06-19 | 北京泛生子基因科技有限公司 | 一种审核高通量测序基因变异检测结果的方法 |
CN113539357A (zh) * | 2021-06-10 | 2021-10-22 | 阿里巴巴新加坡控股有限公司 | 基因检测方法、模型训练方法、装置、设备及系统 |
CN113539357B (zh) * | 2021-06-10 | 2024-04-30 | 阿里巴巴达摩院(杭州)科技有限公司 | 基因检测方法、模型训练方法、装置、设备及系统 |
CN115620802A (zh) * | 2022-09-02 | 2023-01-17 | 蔓之研(上海)生物科技有限公司 | 一种基因数据的处理方法及系统 |
CN115620802B (zh) * | 2022-09-02 | 2023-12-05 | 蔓之研(上海)生物科技有限公司 | 一种基因数据的处理方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
KR20210116454A (ko) | 2021-09-27 |
CN109994155B (zh) | 2021-08-20 |
TWI748263B (zh) | 2021-12-01 |
SG11202011523VA (en) | 2020-12-30 |
WO2020199336A1 (zh) | 2020-10-08 |
JP7064654B2 (ja) | 2022-05-10 |
JP2022500773A (ja) | 2022-01-04 |
TW202036582A (zh) | 2020-10-01 |
US20210082539A1 (en) | 2021-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109994155A (zh) | 一种基因变异识别方法、装置和存储介质 | |
CN109979531A (zh) | 一种基因变异识别方法、装置和存储介质 | |
Schraiber et al. | Methods and models for unravelling human evolutionary history | |
US10600217B2 (en) | Methods for the graphical representation of genomic sequence data | |
US20170178017A1 (en) | Systems and methods that formulate embeddings of problems for solving by a quantum processor | |
US9875215B2 (en) | Systems and methods that formulate problems for solving by a quantum processor using hardware graph decomposition | |
CN109979530A (zh) | 一种基因变异识别方法、装置和存储介质 | |
US20170329899A1 (en) | Display of estimated parental contribution to ancestry | |
BR112021004545A2 (pt) | sistema de determinação ancestral global | |
WO2021139316A1 (zh) | 建立表情识别模型方法、装置、计算机设备及存储介质 | |
CN110262957B (zh) | 一种相似程序间测试用例的重用方法及其实现系统 | |
CN111582315A (zh) | 样本数据处理方法、装置及电子设备 | |
JP6665999B2 (ja) | データ処理装置、決定木生成方法、識別装置及びプログラム | |
CN103294932A (zh) | 用于碱基序列分析的参考序列处理系统及方法 | |
CN109783805A (zh) | 一种网络社区用户识别方法及装置 | |
CN111370055B (zh) | 内含子保留预测模型建立方法及其预测方法 | |
WO2023124779A1 (zh) | 基于三代测序数据检测点突变的分析方法和装置 | |
CN110348022A (zh) | 一种相似度分析的方法、装置、存储介质及电子设备 | |
CN107430600A (zh) | 可扩展的web数据提取 | |
CN115167965A (zh) | 交易进度条的处理方法及装置 | |
US11177018B2 (en) | Stable genes in comparative transcriptomics | |
CN107622184A (zh) | 氨基酸可信度和修饰位点定位的评估方法 | |
Häuser | Quantitative Analysis and Characterization of Natural Language Evolution Datasets | |
KUMAR | Analysis of Machine Learning Based Big Data Analytics | |
Mohammed et al. | DNA base-calling using polynomial classifiers |
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: 40007439 Country of ref document: HK |
|
GR01 | Patent grant | ||
GR01 | Patent grant |