CN108388773A - A kind of identification method of tumor neogenetic antigen - Google Patents
A kind of identification method of tumor neogenetic antigen Download PDFInfo
- Publication number
- CN108388773A CN108388773A CN201810101691.2A CN201810101691A CN108388773A CN 108388773 A CN108388773 A CN 108388773A CN 201810101691 A CN201810101691 A CN 201810101691A CN 108388773 A CN108388773 A CN 108388773A
- Authority
- CN
- China
- Prior art keywords
- mutation
- sequence
- amino acid
- hla
- comparison result
- 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
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
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Engineering & Computer Science (AREA)
- Biotechnology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Analytical Chemistry (AREA)
- Bioinformatics & Computational Biology (AREA)
- Chemical & Material Sciences (AREA)
- Evolutionary Biology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Genetics & Genomics (AREA)
- Molecular Biology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Peptides Or Proteins (AREA)
Abstract
The present invention provides a kind of methods for the tumour specific antigen that individual specimen can be analyzed from NGS data;The present invention can utilize NGS data, the tumour specific antigen for fast and accurately analyzing individual specimen provides reference frame for the diagnosis and analysis of doctor, and can synchronize the comparison for carrying out normal cell and tumour cell and mankind's reference gene group, analysis time is short, and analysis efficiency is high.
Description
Technical field
The present invention relates to a kind of technical fields for identifying tumor neogenetic antigen using bioinformatics method.
Background technology
In immunotherapy of tumors, the neoantigen (neoantigen) that identification tumor tissue cell generates is to determine that downstream is faced
The committed step of bed treatment.Normal cell is in Carcinogenesis, since inhereditary material changes, leads to its DNA sequence dna and its
His normal cell is variant, therefore will produce tumor associated antigen (expression of tumour cell height) or tumour specific antigen and (only exist
It is expressed in tumour cell).These antigens can theoretically be resisted due to the specific epitope with mark tumour cell
Then original is combined, and then activate T cell in human leukocyte antigen (HLA) identification on delivery cell with TCR, start immune anti-
It answers, is the potential target spot of immunotherapy of tumors.Pass through the analysis to patient's tumor tissues and normal structure NGS data, Ke Yijian
Make somatic mutation of the tissue, such as point mutation, insertion and deletion mutation, structure variation, Gene Fusion etc..From these mutation
Event can predict the variation that tumour cell occurs in downstream transcription and expression process, and then deduce that may be present new anti-
Original provides foundation for clinical treatment.
Invention content
The purpose of the present invention is to provide a kind of tumour specific antigens that individual specimen can be analyzed from NGS data
Method.
The technical solution adopted by the present invention to solve the technical problems is:
A kind of identification method of tumor neogenetic antigen, includes the following steps:
S1:Obtain the normal cell DNA sequencing result and DNA of tumor cell sequencing result of tumor patient;
S2:Respectively by the short read sequence of normal cell sequencing and the short read sequence alignment of tumour cell sequencing to the mankind
Reference gene group obtains the comparison result of the comparison result and tumour cell of normal cell respectively;
S3:All believable mutation are detected from the comparison result of short read and carry out functional annotation, and mutation includes embryonal system
Mutation and somatic mutation,
S4:According to the annotation of mutation as a result, being directed to the corresponding transcript sequence of each catastrophic event editor, a packet is obtained
Polypeptide sequence is combined the length needed to intercept into small peptide section by the polypeptide sequence containing all mutation according to HLA partings;
S5:Identify that the HLA partings of HLA (human leukocyte antigen) partings and tumour cell of normal cell, prediction S4 obtain
HLA partings, the binding ability of the HLA partings of small peptide and tumour cell of the small peptide and normal cell that obtain,;
S6:It is preferred that going out and there is candidate data of the small peptide of binding ability as tumour antigen with tumour cell HLA partings.
Further, after step S2 obtains comparison result, Quality Control is carried out to comparison result, the content of Quality Control includes:Entirely catch
Obtain the coverage in region, average sequencing depth, comparison rate, the proportion of repetitive sequence and the proportion for uniquely comparing short read;If than
Quality Control requirement is met to result, then enters S3;If comparison result does not meet Quality Control requirement, sample is reacquired, repeats step
S1 and S2.
Further, the step S2 comparison results obtained are optimized, optimization includes being directed to insertion and deletion region again
It compares and the amendment of target area base quality, the optimization for obtaining the optimization comparison result and normal cell of tumour cell compares knot
Fruit.
Further, in step S2, the short read sequence alignment of normal cell sequencing result to mankind's reference gene group and swells
The short read sequence alignment of oncocyte sequencing result is synchronized to mankind's reference gene group and is carried out.
Further, the range of step S3 functional annotations includes sudden change region, gene, transcript, base variation, amino acid change
Change, clinvar, thousand human genomes, esp6500, dbsnp, cosmic, polyphen, sift, Membrane protein conformation, cancer is related
Gene.
Further, the method that step S3 detects all believable mutation is:
S3.1, the optimization comparison result for optimizing comparison result and normal cell for obtaining tumour cell, detect tumour cell
With the germline mutation site of normal cell;
S3.2, the believable mutation of acquisition is filtered to obtained germline mutation site and somatic mutation site;
S3.3, the somatic mutation filtered out and germline mutation are annotated.
Further, step S4 includes that the method for the polypeptide sequence of all mutation is:
S4.1, normal nucleotide sequence is found out from ensembl databases according to transcript number, according to the work(of mutation
It can annotate and the nucleotide of corresponding position on nucleotide sequence is made into modification, obtain a nucleotides sequence for including all mutation
Row;That provide such as functional annotation is G100A, then the 100th bases G for splicing nucleotide sequence is revised as A.
S4.2, the transcript nucleotide sequence comprising all mutation is translated into polypeptide sequence;
S4.3, sequence obtain sporting in polypeptide sequence and work as premutation, and centered on being currently mutated, interception n is a forward
Amino acid and backward m amino acid of interception obtain polypeptide sequence, the maximum length that n can present for HLA partings, and m is HLA partings
The maximum length or m that can be presented are from when premutation to the length of first terminator codon;
S4.4, in the small peptide that polypeptide sequence successively intercepted length is N, N is that HLA partings combine the length needed;To include
When N small peptide of premutation is used for step S6.Such as:It is 8, N 8 that HLA partings, which combine the length needed,;Since mutated site
7 amino acid are intercepted forward, form the small peptide that a length is 8 amino acid with the mutated site;Since mutated site to
6 amino acid of preceding interception intercept 1 amino acid backward, and it is 8 amino that this 7 amino acid form the 2nd article of length with mutated site
The small peptide of acid, and so on, 8 small peptides containing the mutation can be obtained altogether.
Further, the method for formation nucleotide sequence is in step S4.1:It is found out from database according to transcript number
CDS sequences, the sequence that upper 3 ' UTR region is spliced after CDS sequences form nucleotide sequence.
Further, the rule of the interception of step S4.3 is:
A, for point mutation:Centered on the position of mutation, n amino acid is intercepted forward, intercepts m amino backward
Acid, the longest peptide fragment that n=m=HLA partings can present, if when leading portion or back segment curtailment, how many cuts how many;If
Point mutation belongs to stop loss, then m is from when premutation to the length of first terminator codon;
B, for non-frameshift mutation:Non- frameshit insertion will intercept n amino forward from first amino acid of insetion sequence
Acid;M amino acid, the longest peptide fragment that n=m=HLA partings can present are intercepted backward from the last one amino acid of insetion sequence;
Centered on insetion sequence, m amino acid after n amino acid and insetion sequence before insetion sequence, insetion sequence and
Collectively constitute polypeptide sequence;
Non- frameshift deletion then centered on deletion segment, respectively forwardly intercepts n amino acid, intercepts m amino backward
Acid, the longest peptide fragment that n=m=HLA partings can present;
C, for frameshift mutation:Centered on first amino acid for starting frameshift mutation, n amino acid, n are intercepted forward
The longest peptide fragment that=HLA partings can present;M amino acid is intercepted backward, and m is from when premutation to first terminator codon
Length, i.e., intercept backward to first terminator codon.
The advantage of the invention is that:
1. entire analytic process, since fastq files, user is without preparing other input files.
2. all analytical procedures have corresponding Quality Control step, the accuracy of result is improved.
3. the identification for mutational site is more comprehensive, not only the case where consideration individual cells, it is also considered that monoblock tissue
Mutation distribution.
4. there is multigroup confirmation for learning analysis result.
5. committed step is all carried out at the same time there are many algorithm, it can both be mutually authenticated, and also reduce the false negative of result.
6. having carried out algorithm optimization and parallel processing for taking longer step, the speed of single sample analysis is accelerated
Degree.
Description of the drawings
Fig. 1 is the flow chart of the present invention.
Specific implementation mode
A kind of identification method of tumor neogenetic antigen, includes the following steps:
S1:Obtain the normal cell DNA sequencing result and DNA of tumor cell sequencing result of tumor patient;
S2:Respectively by the short read sequence ratio of the short read sequence of normal cell sequencing result and tumour cell sequencing result
To arriving mankind's reference gene group, the comparison result of the comparison result and tumour cell of normal cell is obtained respectively;
S3:All believable mutation are detected from the comparison result of short read and carry out functional annotation, and mutation includes embryonal system
Mutation and somatic mutation, the embryo of tumour cell and normal cell is detected using the haplotypecaller functions of GATK respectively
It is mutational site;Detect the somatic mutation position of tumour cell and normal cell respectively using tools such as Mutect (1.1.7)
Point;.
The range of functional annotation include sudden change region, gene, transcript, base variation, amino acid variation, clinvar,
Thousand human genomes, esp6500, dbsnp, cosmic, polyphen, sift, Membrane protein conformation, cancer related gene;
S4:According to mutation annotation as a result, obtain a polypeptide sequence for including all mutation, by polypeptide sequence according to
HLA partings combine the length needed to intercept into small peptide section;
S4.1, it is numbered from ensembl databases according to transcript find out corresponding CDS Region Nucleotides sequence first, so
The sequence for splicing upper 3 ' UTR region afterwards generates the transcript sequence of wild type.According to the functional annotation of mutation by nucleotide sequence
The nucleotide of upper corresponding position makes modification, obtains a nucleotide sequence for including all mutation.As functional annotation provides
It is G100A, then the 100th bases G for splicing nucleotide sequence is revised as A.
S4.2, the nucleotide sequence comprising all mutation is translated into polypeptide sequence;
S4.3, sequence obtain sporting in polypeptide sequence and work as premutation, and centered on being currently mutated, interception n is a forward
Amino acid and backward m amino acid of interception obtain polypeptide sequence, the maximum length that n can present for HLA partings, and m is HLA partings
The maximum length or m that can be presented are from when premutation to the length of first terminator codon;
S4.4, in the small peptide that polypeptide sequence successively intercepted length is N, N is that HLA partings combine the length needed;To include
When the small peptide of premutation is used for step S5.By taking N=11 as an example, 11 will be obtained and include the small peptide when premutation.If polypeptide sequence
Also include other mutation, and have certain small peptides simultaneously comprising current in N small peptide in row except when other than premutation
Mutation and other mutation, then become the small peptide with mutation combination.
S5:Identify that the HLA partings of HLA (human leukocyte antigen) partings and tumour cell of normal cell, prediction S4 obtain
HLA partings, the binding ability of the HLA partings of small peptide and tumour cell of the small peptide and normal cell that obtain.
The HLA of tumour cell is likely to occur change, compares the difference of normal cell HLA partings and tumour cell HLA partings
It is different, compare wild type peptide fragment (normal cell expression) and HLA binding abilities and saltant type peptide fragment (tumor cells expression) and HLA
The difference of binding ability knows the difference of tumour cell and normal cell, is the preferred accumulation foundation of follow-up neoantigen.
The HLA partings of normal structure sample and tumor tissues sample are predicted using tools such as SOAP-HLA, obtain HLA_I types
(A, B, C) and HLA_II types (DR, DQ, DP).
S6:It is preferred that go out has the small peptide of binding ability as the candidate as tumour antigen with tumour cell HLA partings
Data.
Affinity prediction is carried out using multiple HLA partings affinity forecasting softwares respectively, is such as directed to HLA I types
(netMHC4.0 etc.), for HLA II types (netMHCII 2.2 etc.), all parameters are acquiescence.Prediction result can provide often
The affinity of one HLA partings and each nascent polypeptide filters out the relatively strong (NM of wherein affinity<500) polypeptide.
In one embodiment, after step S2 obtains comparison result, Quality Control, the content packet of Quality Control are carried out to comparison result
Contain:The coverage of entire capture region, average sequencing depth, comparison rate, the proportion of repetitive sequence and unique ratio for comparing reads
Weight;If comparison result meets Quality Control requirement, enter S3;If comparison result does not meet Quality Control requirement, sample is reacquired, weight
Multiple step S1 and S2.
In general we require Q30>80%, i.e. sequencing error rate is more than that 80%, GC contains less than millesimal reads
Amount 50% or so, chain skewed popularity are less than 70%, then it is assumed that normal cell sequencing result and tumour cell the sequencing knot currently obtained
The analysis of step after fruit can be used in.
In one embodiment, the step S2 comparison results obtained are optimized, optimization includes being directed to insertion and deletion area
Domain again compare and target area base quality amendment, obtain tumour cell optimization comparison result and normal cell it is excellent
Change comparison result.
In step S2, the short read sequence alignment of normal cell sequencing result is surveyed to mankind's reference gene group and tumour cell
The short read sequence alignment of sequence result is synchronized to mankind's reference gene group and is carried out.
The method that step S3 detects all believable mutation is:
S3.1, the optimization comparison result for optimizing comparison result and normal cell for obtaining tumour cell, detect tumour cell
With the germline mutation site of normal cell;
S3.2, the believable mutation of acquisition is filtered to obtained germline mutation site and somatic mutation site;
S3.3, the somatic mutation filtered out and germline mutation are annotated.
In one embodiment, when mutation type is point mutation, the rule of the interception of step S4.3 is:With mutation
Centered on position, n amino acid is intercepted forward, intercepts m amino acid, the longest peptide that n=m=HLA partings can present backward
Section, if when leading portion or back segment curtailment, how many cuts how many;If point mutation, which belongs to stop loss, (loses termination codon
Son), then m is from when premutation to the length of first terminator codon.
In this way, for each mutation, it can simulate and be possible to the small peptide containing this mutation in tumour cell,
It is predicted finally by with the binding force of patient's HLA partings, filters out the small peptide of binding ability, the candidate as tumour antigen
Data.
Currently, it is that HLA II types combine 16 amino acid needed, therefore, n=that HLA partings, which combine the maximum length needed,
M=16.
Combine the length needed for for 8, N=8 with HLA partings;Intercept 7 amino forward since current point mutation
Acid is mutated the small peptide to form that a length is 8 sites with current point;Intercept 6 amino acid forward since current point mutation,
1 amino acid is intercepted backward, this 7 amino acid are mutated the small peptide to form that the 2nd article of length is 8 sites with current point, with such
It pushes away, 8 small peptides for containing current point mutation can be obtained altogether.
In one embodiment, when mutation type is non-frameshift mutation, the rule of the interception of step S4.3 is:Non- frameshit is inserted
N amino acid will be intercepted forward from first amino acid of insetion sequence by entering;Backward from the last one amino acid of insetion sequence
Intercept m amino acid, the longest peptide fragment that n=m=HLA partings can present;Centered on insetion sequence, insetion sequence is inserted into sequence
M amino acid after n amino acid and insetion sequence before row and collectively constitute polypeptide sequence.
Non- frameshift deletion then centered on deletion segment, respectively forwardly intercepts n amino acid, intercepts m amino backward
Acid, the longest peptide fragment that n=m=HLA partings can present.
Currently, it is that HLA II types combine 16 amino acid needed, therefore, n=that HLA partings, which combine the maximum length needed,
M=16.
In one embodiment, when mutation type is frameshift mutation, the rule of the interception of step S4.3 is:From beginning frameshit
Centered on first amino acid of mutation, n amino acid, the longest peptide fragment that n=HLA partings can present are intercepted forward;Backward
Intercept m amino acid, m is from when premutation is to the length of first terminator codon, i.e., intercepts close to first termination backward
Numeral.
Currently, it is that HLA II types combine 16 amino acid needed, therefore, n=that HLA partings, which combine the maximum length needed,
16。
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
Member, without departing from the inventive concept of the premise, can also make several improvements and modifications, these improvements and modifications also should be regarded as
In the scope of the present invention.
Claims (9)
1. a kind of identification method of tumor neogenetic antigen, which is characterized in that include the following steps:
S1:Obtain the normal cell DNA sequencing result and DNA of tumor cell sequencing result of tumor patient;
S2:The short read sequence alignment of the short read sequence of normal cell sequencing and tumour cell sequencing is referred to the mankind respectively
Genome obtains the comparison result of the comparison result and tumour cell of normal cell respectively;
S3:All believable mutation are detected from the comparison result of short read and carry out functional annotation, and mutation includes germline mutation
And somatic mutation,
S4:According to the annotation of mutation as a result, being directed to the corresponding transcript sequence of each catastrophic event editor, it includes institute to obtain one
Polypeptide sequence is combined the length needed to intercept into small peptide section by the polypeptide sequence for having mutation according to HLA partings;
S5:Identify that the HLA partings of HLA (human leukocyte antigen) partings and tumour cell of normal cell, prediction S4 obtain
The binding ability of the HLA partings of small peptide and normal cell, the HLA partings of small peptide and tumour cell,;
S6:It is preferred that going out and there is candidate data of the small peptide of binding ability as tumour antigen with tumour cell HLA partings.
2. a kind of identification method of tumor neogenetic antigen according to claim 1, it is characterised in that:Step S2 obtains ratio
After result, Quality Control is carried out to comparison result, the content of Quality Control includes:The coverage of entire capture region, average sequencing depth,
Comparison rate, the proportion of repetitive sequence and the proportion for uniquely comparing short read;If comparison result meets Quality Control requirement, enter S3;
If comparison result does not meet Quality Control requirement, sample is reacquired, repeats step S1 and S2.
3. a kind of identification method of tumor neogenetic antigen according to claim 2, it is characterised in that:Step S2 is obtained
Comparison result optimize, optimization include for insertion and deletion region compare again and target area base quality is repaiied
Just, the optimization comparison result of the optimization comparison result and normal cell of tumour cell is obtained.
4. a kind of identification method of tumor neogenetic antigen according to claim 3, it is characterised in that:In step S2, just
Short read sequence of the short read sequence alignment of normal cell sequencing result to mankind's reference gene group and tumour cell sequencing result
It compares mankind's reference gene group and synchronizes progress.
5. a kind of identification method of tumor neogenetic antigen according to claim 1, it is characterised in that:Step S3 functions are noted
The range released include sudden change region, gene, transcript, base variation, amino acid variation, clinvar, thousand human genomes,
Esp6500, dbsnp, cosmic, polyphen, sift, Membrane protein conformation, cancer related gene.
6. a kind of identification method of tumor neogenetic antigen according to claim 5, it is characterised in that:Step S3 detects institute
There is the method for believable mutation to be:
S3.1, obtain tumour cell optimization comparison result and normal cell optimization comparison result, detect tumour cell and just
The germline mutation site of normal cell;
S3.2, the believable mutation of acquisition is filtered to obtained germline mutation site and somatic mutation site;
S3.3, the somatic mutation filtered out and germline mutation are annotated.
7. a kind of identification method of tumor neogenetic antigen according to claim 1, it is characterised in that:Step S4 includes institute
There is the method for the polypeptide sequence of mutation to be:
S4.1, normal nucleotide sequence is found out from ensembl databases according to transcript number, is noted according to the function of mutation
It releases and the nucleotide of corresponding position on nucleotide sequence is made into modification, obtain a nucleotide sequence for including all mutation;Such as
That functional annotation provides is G100A, then the 100th bases G for splicing nucleotide sequence is revised as A.
S4.2, the transcript nucleotide sequence comprising all mutation is translated into polypeptide sequence;
S4.3, sequence obtain sporting in polypeptide sequence and work as premutation, centered on being currently mutated, intercept n amino forward
M amino acid acquisition polypeptide sequence of acid and backward interception, the maximum length that n can present for HLA partings, m is HLA partings institute energy
The maximum length or m of presentation are from when premutation to the length of first terminator codon;
S4.4, the length needed for the combination of HLA partings in the small peptide that polypeptide sequence successively intercepted length is N, N;To include current
N small peptide of mutation is used for step S6.
8. a kind of identification method of tumor neogenetic antigen according to claim 7, it is characterised in that:Shape in step S4.1
It is at the method for nucleotide sequence:CDS sequences are found out from database according to transcript number, and upper 3 ' are spliced after CDS sequences
The sequence of UTR region forms nucleotide sequence.
9. a kind of identification method of tumor neogenetic antigen according to claim 7, it is characterised in that:Step S4.3's cuts
The rule taken is:
A, for point mutation:Centered on the position of mutation, n amino acid is intercepted forward, intercepts m amino acid, n backward
The longest peptide fragment that=m=HLA partings can present, if when leading portion or back segment curtailment, how many cuts how many;As fruit dot is prominent
Change belongs to stop loss, then m is from when premutation to the length of first terminator codon;
B, for non-frameshift mutation:Non- frameshit insertion will intercept n amino acid forward from first amino acid of insetion sequence;From
The last one amino acid of insetion sequence intercepts m amino acid, the longest peptide fragment that n=m=HLA partings can present backward;To insert
Centered on entering sequence, m amino acid after n amino acid and insetion sequence before insetion sequence, insetion sequence and common
Form polypeptide sequence;
Non- frameshift deletion then centered on deletion segment, respectively forwardly intercepts n amino acid, intercepts m amino acid, n=m backward
The longest peptide fragment that=HLA partings can present;
C, for frameshift mutation:Centered on first amino acid for starting frameshift mutation, n amino acid, n=are intercepted forward
The longest peptide fragment that HLA partings can present;M amino acid is intercepted backward, and m is from when premutation to first terminator codon
Length is intercepted backward to first terminator codon.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810101691.2A CN108388773B (en) | 2018-02-01 | 2018-02-01 | A kind of identification method of tumor neogenetic antigen |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810101691.2A CN108388773B (en) | 2018-02-01 | 2018-02-01 | A kind of identification method of tumor neogenetic antigen |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108388773A true CN108388773A (en) | 2018-08-10 |
CN108388773B CN108388773B (en) | 2019-11-19 |
Family
ID=63074863
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810101691.2A Active CN108388773B (en) | 2018-02-01 | 2018-02-01 | A kind of identification method of tumor neogenetic antigen |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108388773B (en) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109584960A (en) * | 2018-12-14 | 2019-04-05 | 上海鲸舟基因科技有限公司 | Predict the method, apparatus and storage medium of tumor neogenetic antigen |
CN109584966A (en) * | 2019-01-08 | 2019-04-05 | 杭州纽安津生物科技有限公司 | A kind of design method of tumour general vaccines and its application in cancer of pancreas |
CN109637585A (en) * | 2018-12-27 | 2019-04-16 | 北京优迅医学检验实验室有限公司 | The antidote and device of depth is sequenced |
CN110179972A (en) * | 2019-07-10 | 2019-08-30 | 杭州纽安津生物科技有限公司 | Self assembly polypeptide neoantigen tumor vaccine and preparation method thereof |
CN110265086A (en) * | 2019-07-04 | 2019-09-20 | 北京肿瘤医院(北京大学肿瘤医院) | Gene detection method and device |
CN110277135A (en) * | 2019-08-10 | 2019-09-24 | 杭州新范式生物医药科技有限公司 | A kind of method and system based on expected effect selection individuation knubble neoantigen |
CN110322925A (en) * | 2019-07-18 | 2019-10-11 | 杭州纽安津生物科技有限公司 | A method of prediction fusion generates neoantigen |
CN110349627A (en) * | 2019-06-27 | 2019-10-18 | 杭州纽安津生物科技有限公司 | The design method and its Automation Design product of polypeptide vaccine sequence |
CN110491450A (en) * | 2019-08-23 | 2019-11-22 | 深圳市新合生物医疗科技有限公司 | Tumor neogenetic antigen prediction platform and its application in neoantigen vaccine development system |
CN110675913A (en) * | 2019-01-16 | 2020-01-10 | 倍而达药业(苏州)有限公司 | Screening method of tumor neoantigen based on HLA typing and structure |
CN110706747A (en) * | 2019-09-17 | 2020-01-17 | 北京橡鑫生物科技有限公司 | Method and device for detecting tumor neoantigen polypeptide |
WO2020187143A1 (en) * | 2019-03-15 | 2020-09-24 | 痕准生物科技有限公司 | Method for identifying neoantigens |
CN113012756A (en) * | 2021-03-08 | 2021-06-22 | 杭州纽安津生物科技有限公司 | Screening method of individual tumor neoantigen peptide and vaccine preparation thereof |
CN113035276A (en) * | 2021-03-11 | 2021-06-25 | 深圳荻硕贝肯精准医学有限公司 | Method and system for analyzing heterozygous deletion of human HLA chromosome region |
CN116825188A (en) * | 2023-06-25 | 2023-09-29 | 北京泛生子基因科技有限公司 | Method, device and computer readable storage medium for identifying tumor neoantigen at multiple groups of chemical layers based on high-throughput sequencing technology |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104059966A (en) * | 2014-05-20 | 2014-09-24 | 吴松 | STAG2 gene mutant sequence and detection method thereof as well as use of STAG2 gene mutation in detecting bladder cancer |
US20170016075A1 (en) * | 2015-07-14 | 2017-01-19 | Personal Genome Diagnostics, Inc. | Neoantigen analysis |
CN107523563A (en) * | 2017-09-08 | 2017-12-29 | 杭州和壹基因科技有限公司 | A kind of Bioinformatics method for Circulating tumor DNA analysis |
-
2018
- 2018-02-01 CN CN201810101691.2A patent/CN108388773B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104059966A (en) * | 2014-05-20 | 2014-09-24 | 吴松 | STAG2 gene mutant sequence and detection method thereof as well as use of STAG2 gene mutation in detecting bladder cancer |
US20170016075A1 (en) * | 2015-07-14 | 2017-01-19 | Personal Genome Diagnostics, Inc. | Neoantigen analysis |
CN107523563A (en) * | 2017-09-08 | 2017-12-29 | 杭州和壹基因科技有限公司 | A kind of Bioinformatics method for Circulating tumor DNA analysis |
Non-Patent Citations (2)
Title |
---|
武雅蓉: "测序数据质控统计--微不足道?意义重大!", 《HTTP://JU.OUTOFMEMORY.CN/ENTRY/330230》 * |
青萍你好: "INDEL的重新比对和碱基质量分数的重新校", 《一周一PAPER,一周一技术》 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109584960A (en) * | 2018-12-14 | 2019-04-05 | 上海鲸舟基因科技有限公司 | Predict the method, apparatus and storage medium of tumor neogenetic antigen |
CN109637585A (en) * | 2018-12-27 | 2019-04-16 | 北京优迅医学检验实验室有限公司 | The antidote and device of depth is sequenced |
CN109584966A (en) * | 2019-01-08 | 2019-04-05 | 杭州纽安津生物科技有限公司 | A kind of design method of tumour general vaccines and its application in cancer of pancreas |
CN109584966B (en) * | 2019-01-08 | 2019-09-20 | 杭州纽安津生物科技有限公司 | A kind of design method and cancer of pancreas general vaccines of tumour general vaccines |
CN110675913B (en) * | 2019-01-16 | 2022-04-12 | 倍而达药业(苏州)有限公司 | Screening method of tumor neoantigen based on HLA typing and structure |
CN110675913A (en) * | 2019-01-16 | 2020-01-10 | 倍而达药业(苏州)有限公司 | Screening method of tumor neoantigen based on HLA typing and structure |
CN113316818B (en) * | 2019-03-15 | 2024-04-02 | 痕准生物医学技术(厦门)有限公司 | Method for identifying neoantigen |
WO2020187143A1 (en) * | 2019-03-15 | 2020-09-24 | 痕准生物科技有限公司 | Method for identifying neoantigens |
CN113316818A (en) * | 2019-03-15 | 2021-08-27 | 痕准生物科技有限公司 | Method for identifying neoantigens |
CN110349627A (en) * | 2019-06-27 | 2019-10-18 | 杭州纽安津生物科技有限公司 | The design method and its Automation Design product of polypeptide vaccine sequence |
CN110349627B (en) * | 2019-06-27 | 2021-08-27 | 杭州纽安津生物科技有限公司 | Design method of polypeptide vaccine sequence and its automatic design product |
CN110265086A (en) * | 2019-07-04 | 2019-09-20 | 北京肿瘤医院(北京大学肿瘤医院) | Gene detection method and device |
CN110179972A (en) * | 2019-07-10 | 2019-08-30 | 杭州纽安津生物科技有限公司 | Self assembly polypeptide neoantigen tumor vaccine and preparation method thereof |
CN110322925A (en) * | 2019-07-18 | 2019-10-11 | 杭州纽安津生物科技有限公司 | A method of prediction fusion generates neoantigen |
CN110322925B (en) * | 2019-07-18 | 2021-09-03 | 杭州纽安津生物科技有限公司 | Method for predicting generation of neoantigen by fusion gene |
CN110277135A (en) * | 2019-08-10 | 2019-09-24 | 杭州新范式生物医药科技有限公司 | A kind of method and system based on expected effect selection individuation knubble neoantigen |
CN110491450A (en) * | 2019-08-23 | 2019-11-22 | 深圳市新合生物医疗科技有限公司 | Tumor neogenetic antigen prediction platform and its application in neoantigen vaccine development system |
CN110706747B (en) * | 2019-09-17 | 2021-09-07 | 北京橡鑫生物科技有限公司 | Method and device for detecting tumor neoantigen polypeptide |
CN110706747A (en) * | 2019-09-17 | 2020-01-17 | 北京橡鑫生物科技有限公司 | Method and device for detecting tumor neoantigen polypeptide |
CN113012756A (en) * | 2021-03-08 | 2021-06-22 | 杭州纽安津生物科技有限公司 | Screening method of individual tumor neoantigen peptide and vaccine preparation thereof |
WO2022188280A1 (en) * | 2021-03-08 | 2022-09-15 | 杭州纽安津生物科技有限公司 | Method for screening individual tumor neoantigen peptide, and vaccine formulation thereof |
CN113035276A (en) * | 2021-03-11 | 2021-06-25 | 深圳荻硕贝肯精准医学有限公司 | Method and system for analyzing heterozygous deletion of human HLA chromosome region |
CN116825188A (en) * | 2023-06-25 | 2023-09-29 | 北京泛生子基因科技有限公司 | Method, device and computer readable storage medium for identifying tumor neoantigen at multiple groups of chemical layers based on high-throughput sequencing technology |
CN116825188B (en) * | 2023-06-25 | 2024-04-09 | 北京泛生子基因科技有限公司 | Method, device and computer readable storage medium for identifying tumor neoantigen at multiple groups of chemical layers based on high-throughput sequencing technology |
Also Published As
Publication number | Publication date |
---|---|
CN108388773B (en) | 2019-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108388773B (en) | A kind of identification method of tumor neogenetic antigen | |
CN108491689B (en) | Tumour neoantigen identification method based on transcript profile | |
CN109033749B (en) | Tumor mutation load detection method, device and storage medium | |
Anagnostou et al. | Evolution of neoantigen landscape during immune checkpoint blockade in non–small cell lung cancer | |
CN109584960B (en) | Method, device and storage medium for predicting tumor neoantigen | |
Izraelson et al. | Comparative analysis of murine T‐cell receptor repertoires | |
CN108701173A (en) | System, composition and method for finding the prediction MSI and new epitope sensitive to checkpoint inhibitor | |
CN106599616B (en) | Ultralow frequency mutational site determination method based on duplex-seq | |
CN107391965A (en) | A kind of lung cancer somatic mutation determination method based on high throughput sequencing technologies | |
CN110277135B (en) | Method and system for selecting individualized tumor neoantigen based on expected curative effect | |
CN111415707B (en) | Prediction method of clinical individuation tumor neoantigen | |
Borden et al. | Cancer neoantigens: challenges and future directions for prediction, prioritization, and validation | |
CN110621785B (en) | Method and device for haplotyping diploid genome based on three-generation capture sequencing | |
CN112164423B (en) | Fusion gene detection method, device and storage medium based on RNAseq data | |
CN111676288B (en) | System for predicting lung adenocarcinoma patient prognosis and application thereof | |
CN110752041A (en) | Method, device and storage medium for predicting neoantigen based on next generation sequencing | |
JPWO2020023420A5 (en) | ||
KR20230165259A (en) | Identification of clonal neoantigens and their use | |
CN114694750B (en) | Single-sample tumor somatic mutation distinguishing and TMB (tumor necrosis factor) detecting method based on NGS (Next Generation broadcasting) platform | |
Cao et al. | Identification and development of a novel 4-gene immune-related signature to predict osteosarcoma prognosis | |
CN114929899A (en) | Method and system for screening new antigen and application thereof | |
CN112210596B (en) | Tumor neoantigen prediction method based on gene fusion event and application thereof | |
CN111192632B (en) | Method and device for extracting gene fusion immunotherapy new antigen by integrating DNA and RNA deep sequencing data | |
WO2024051097A1 (en) | Neoantigen identification method and device for tumor-specific circular rnas, apparatus and medium | |
CN115424740B (en) | Tumor immunotherapy effect prediction system based on NGS and deep learning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |