CN107704727A - Neoantigen Activity Prediction and sort method based on tumour neoantigen characteristic value - Google Patents

Neoantigen Activity Prediction and sort method based on tumour neoantigen characteristic value Download PDF

Info

Publication number
CN107704727A
CN107704727A CN201711071334.8A CN201711071334A CN107704727A CN 107704727 A CN107704727 A CN 107704727A CN 201711071334 A CN201711071334 A CN 201711071334A CN 107704727 A CN107704727 A CN 107704727A
Authority
CN
China
Prior art keywords
neoantigen
tumour
prediction
mhc
activity
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
CN201711071334.8A
Other languages
Chinese (zh)
Other versions
CN107704727B (en
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.)
Wang Yong
Original Assignee
Hangzhou Wind Intelligent Technology 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 Hangzhou Wind Intelligent Technology Co Ltd filed Critical Hangzhou Wind Intelligent Technology Co Ltd
Priority to CN201711071334.8A priority Critical patent/CN107704727B/en
Publication of CN107704727A publication Critical patent/CN107704727A/en
Application granted granted Critical
Publication of CN107704727B publication Critical patent/CN107704727B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/50Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
    • G01N33/53Immunoassay; Biospecific binding assay; Materials therefor
    • 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
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Urology & Nephrology (AREA)
  • Hematology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biotechnology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Biophysics (AREA)
  • Evolutionary Biology (AREA)
  • Theoretical Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Cell Biology (AREA)
  • Genetics & Genomics (AREA)
  • Microbiology (AREA)
  • Medical Informatics (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
  • Peptides Or Proteins (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

The invention discloses a kind of neoantigen immunocompetence marking based on tumour neoantigen characteristic value and sort method, comprise the following steps:The WGS/WES, RNA seq sequencing datas of tumour normal sample input, the prediction of tumour somatic mutation and annotation, associated eigenvalue calculate;MHC I combinations neoantigen prediction based on tumour somatic mutation, associated eigenvalue calculate;The extraction of neoantigen associated eigenvalue;The setting of neoantigen activity scoring functions;Neoantigen sequence based on neoantigen activity scoring functions.Present invention analysis first calculates the mutation of tumour body cell and completes mutation annotation, calculates partial feature value, then predicts MHC I combination neoantigens, calculate partial feature value;All associated eigenvalues of tumour neoantigen are extracted, and then set neoantigen activity scoring functions, neoantigen is ranked up finally by neoantigen activity scoring functions.This method is highly efficient accurate compared with traditional screening technique, has important application value for tumour immunotherapy.

Description

Neoantigen Activity Prediction and sort method based on tumour neoantigen characteristic value
Technical field
The present invention relates to immunotherapy of tumors field, in particular it relates to a kind of new anti-based on tumour neoantigen characteristic value Former activity marking and sort method.
Background technology
In recent years, immunotherapy of tumors yields unusually brilliant results, and clinical test constantly obtains spy and broken, and cure rate and effective remission rate are held Continuous lifting.The efficient precisely screening of tumour neoantigen is work of crucial importance and basic in tumour immunotherapy, particularly It is particularly important to tumour immunotherapies such as TCR-T/TIL, individuation vaccines.
At present, the scheme passed through at present for the screening technique of tumour neoantigen is two steps:Step 1: based on swollen The WGS/WES data of knurl-normal structure, the instruments such as Mutect/Varscan are called to calculate the gene mutation of tumour cell;Step 2nd, NetMHCpan scheduling algorithms prediction MHC-I combination neoantigens are called.
There is presently no activity of the effective method based on antigen to be ranked up, to lift antigen selection efficiency.Above-mentioned side Due to not carrying out active sequence to the MHC-I combinations neoantigen that prediction obtains in case, therefore can be huge to bringing for experimental verification Big workload, cause tumour neoantigen screening efficiency it is low.
The content of the invention
The present invention is directed to above-mentioned the deficiencies in the prior art, there is provided a kind of based on the new of tumour neoantigen characteristic value Antigen active is given a mark and sort method, can substantially reduce the workload of experimental verification, and further realize the height of tumour neoantigen Effect precisely screening.
The technical proposal of the invention is realized in this way:
Based on tumour neoantigen characteristic value neoantigen immunocompetence prediction and sort method, it is characterised in that including with Lower step:
(1), the input of WGS/WES, RNA-SEQ sequencing data of tumour-normal sample:Input tumour-normal sample Full genome sequencing data WGS or full sequencing of extron group data WES, transcript profile sequencing data RNA-SEQ;
(2), the prediction of tumour somatic mutation and annotation, the calculating of associated eigenvalue:Sequencing based on step (1) input Data, call Varscan or Mutect tool analysis to calculate tumour somatic mutation, call VEP (Variant Effect Prediction) instrument complete mutation annotation, call PyClone, Kallisto, Varscan or Mutect instrument calculate as Lower eigenvalue:Mutator clone ratio, mutator expression value TPM, allelic mutation frequency VAF;
(3) prediction of MHC-I combinations neoantigen, the calculating of associated eigenvalue based on tumour somatic mutation:Based on step (2) tumour somatic mutation and annotation data in, NetMHCpan, Netchop, OptiType instrument prediction MHC-I knots are called Neoantigen is closed, and is calculated such as lower eigenvalue:It is mutated peptide fragment and MHC affinity sequence percentage, unmutated peptide fragment and MHC is affine Power sequence percentage, peptide fragment shearing present efficiency;
(4) extraction of all associated eigenvalues of neoantigen:For the MHC-I combination neoantigens of prediction in step (3), extraction Go out all associated eigenvalues of tumour neoantigen;
(5) setting of neoantigen activity scoring functions:For the neoantigen characteristic value of extraction in step (4), set new anti- Former active scoring functions;
(6) the neoantigen sequence based on neoantigen activity scoring functions:By neoantigen activity scoring functions to neoantigen It is ranked up.
Preferably, in step (4), neoantigen associated eigenvalue includes Rm、A、Rn, E, NC, CL, wherein:
Rm- mutation peptide fragment and MHC affinity sequence percentage, are calculated by NetMHCpan;
A-allelic mutation frequency VAF, is calculated by Varscan/Mutect/Strelka2;
Rn- unmutated peptide fragment and MHC affinity sequence percentage, are calculated by NetMHCpan;
E-mutator expression value TPM, is calculated by Kallisto;
NC-peptide fragment shearing presents efficiency, is calculated by netchop;
CL-mutator clone's ratio, is calculated by pyclone.
Preferably, in step (5), the neoantigen Activity Prediction scoring functions of proposition are:
Neo_Score=abundancedissimilarityclonality;
Clonality=NCCL;
Abundance=L (Rm) Atanh (E/k);
Dissimilarity=(1-L (Rn)/2));
Wherein:L (x)=1/ (1+e5(x-2)), tanh (x) is hyperbolic tangent function;
K is transcript gene expression abundance threshold value, default value 1.
Preferably, in step (6), algorithmic procedure is ranked up such as by neoantigen Activity Prediction function pair neoantigen Under:
A), for the MHC-I combination neoantigens of all predictions, neoantigen Activity Prediction function Neo_score is called to calculate Go out the predicted value of neoantigen activity;
B), the predicted value based on neoantigen activity, is ranked up using quick sorting algorithm to neoantigen;
C), neoantigen ranking results are exported.
The design philosophy of the invention and beneficial effect for employing above-mentioned technical proposal be:
Technical scheme, it is proposed that a kind of tumour neoantigen Activity Prediction scoring functions, based on neoantigen activity Anticipation function carries out active sequence to the MHC-I combinations neoantigen predicted, so as to realize efficiently, accurately tumour neoantigen sieves Choosing.
This function is based on the generation of tumour neoantigen, cutting is transported, neoantigen is combined the design of this complete procedure with MHC, beats Function is divided to be divided into 3 parts, wherein, it is to the efficiency NC of small peptide and newly anti-by being mutated that Clonal (clonality) weighs neoantigen The ratio CL that original is distributed in all tumour cells, it is an important factor for influenceing tumor vaccine curative effect;Abundance (abundance) weighs Amount neoantigen expression quantity and neoantigen are combined with MHC-I and form the efficiency of pMHC compounds, neoantigen mutator expression quantity E Higher, allelic mutation frequency A is higher, binding affinity R between compoundmHigher (IC50 values are smaller), immunogenicity is then got over By force;Dissimilar degree (dissimilarity) weighs difference R of the mutation peptide fragment with corresponding normal peptide fragment affinityn, because maincenter is resistance to Existed by mechanism, the two difference is bigger, and the specificity of neoantigen is stronger, and the side effect for treatment is also smaller.Two in function Individual mapping function is used to normalize (0 to 1) calculated value, and threshold value 2 is that peptide fragment-MHC binding affinities screen threshold value, tanh in L (x) (x/k) when ensureing that neoantigen gene expression abundance exceedes given threshold k, function value changes tend towards stability.
This Activity Prediction function considers the combined influence factor during neoantigen produces, and the antigen after sequence is more intentional Justice and application value.
Brief description of the drawings
Fig. 1 is the neoantigen activity methods and sort method based on tumour neoantigen characteristic value described in the embodiment of the present invention Schematic diagram.
Embodiment
The preferred embodiments of the present invention are illustrated below in conjunction with accompanying drawing, it will be appreciated that described herein preferred real Apply example to be merely to illustrate and explain the present invention, be not intended to limit the present invention.
Embodiment one:
As shown in figure 1, a kind of neoantigen activity methods and sort method based on tumour neoantigen characteristic value, including it is following Step:
Step 101:The input of WGS/WES, RNA-seq sequencing data of tumour-normal sample (uses melanoma patient Sample one mel_21, Science 2015:Carreno B M,Magrini V,Beckerhapak M,et al.Cancer immunotherapy.A dendritic cell vaccine increases the breadth and diversity of melanoma neoantigen-specific T cells.[J].Science,2015,348(6236):803-8.)
Step 102:The prediction of tumour somatic mutation and annotation, the calculating of associated eigenvalue:Based on tumour-normal sample WGS/WES, RNA-seq sequencing data, call the tool analysis such as Varscan/Mutect to calculate tumour somatic mutation, adjust Mutation is completed with VEP (Variant Effect Prediction) instrument to annotate, calling PyClone, Kallisto, Varscan/Mutect instruments are calculated such as lower eigenvalue:Mutator is cloned than CL, mutator expression value TPM, equipotential base Because of frequency of mutation VAF.By taking characteristic value corresponding to active peptide fragment in document as an example, the mel_21 of patient's sample one 3 peptide fragments pair The E (TPM) and A (VAF) value that should be calculated are as follows:
Table 1
Step 103:The prediction of MHC-I combinations neoantigen, the calculating of associated eigenvalue based on tumour somatic mutation:It is based on Tumour somatic mutation and annotation data in step (1), call the prediction of NetMHCpan, Netchop, OptiType instrument MHC-I combination neoantigens, and calculate such as lower eigenvalue:It is mutated peptide fragment and MHC affinity sequence percentage Rm, unmutated peptide fragment With MHC affinity sequence percentage Rn, peptide fragment shearing present efficiency NC.Using in document characteristic value corresponding to active peptide fragment as Example, the R that the mel_21 of patient's sample one 3 peptide fragments correspondingly calculatem,RnIt is as follows with NC values:
Table 2
Step 104:The extraction of neoantigen associated eigenvalue:For the MHC-I combination neoantigens predicted in step 103, carry All associated eigenvalues of tumour neoantigen are taken out, the mel_21 of patient's sample one neoantigen and characteristic value are shown in Table 1, table 2;
Step 105:The setting of neoantigen activity scoring functions:For the neoantigen characteristic value extracted in step 104, setting Neoantigen activity scoring functions;
Step 106:Neoantigen sequence based on neoantigen activity scoring functions:By neoantigen activity scoring functions to new Antigen is ranked up, and table 3 gives the marking (Neo_ of the neoantigen through the checking of PMHC activity experiments in the mel_21 of patient's sample one Score the sequence (Rank)) and in complete or collected works.
Table 3
In the neoantigen that 3 identify, CLNEYHLFL can to stimulate DC cells just to show before using tumor vaccine The immunocompetence of CD8+T cells is activated, remaining 2 are then after using vaccine enhancing immune system ability, are possessed in various degree Immunocompetence.We have seen that KMIGNHLWV and CLNEYHLFL is in sequence Top 3 (total candidate's neoantigen number 94) position. And AMFWSVPTS is either tested in still human tumor microenvironment in vitro, because expression quantity is relatively low, its immunogenicity is determined It is on the weak side with immune response, be in 33 in our ranking results, experimental result and our prediction sequence quite it is identical.
Embodiment two:
As shown in figure 1, a kind of neoantigen activity methods and sort method based on tumour neoantigen characteristic value, including it is following Step:
Step 101:The input of WGS/WES, RNA-seq sequencing data of tumour-normal sample (uses melanoma patient Sample two mel_38, Science 2015:Carreno B M,Magrini V,Beckerhapak M,et al.Cancer immunotherapy.A dendritic cell vaccine increases the breadth and diversity of melanoma neoantigen-specific T cells.[J].Science,2015,348(6236):803-8.)
Step 102:The prediction of tumour somatic mutation and annotation, the calculating of associated eigenvalue:Based on tumour-normal sample WGS/WES, RNA-seq sequencing data, call the tool analysis such as Varscan/Mutect to calculate tumour somatic mutation, adjust Mutation is completed with VEP (Variant Effect Prediction) instrument to annotate, calling PyClone, Kallisto, Varscan/Mutect instruments are calculated such as lower eigenvalue:Mutator is cloned than CL, mutator expression value TPM, equipotential base Because of frequency of mutation VAF;By taking characteristic value corresponding to active peptide fragment in document as an example, the mel_38 of patient's sample two 3 peptide fragments pair The E (TPM) and A (VAF) value that should be calculated are as follows:
Table 4
Step 103:The prediction of MHC-I combinations neoantigen, the calculating of associated eigenvalue based on tumour somatic mutation:It is based on Tumour somatic mutation and annotation data in step (1), call the prediction of NetMHCpan, Netchop, OptiType instrument MHC-I combination neoantigens, and calculate such as lower eigenvalue:It is mutated peptide fragment and MHC affinity sequence percentage Rm, unmutated peptide fragment With MHC affinity sequence percentage Rn, peptide fragment shearing present efficiency NC;Using in document characteristic value corresponding to active peptide fragment as Example, the R that the mel_38 of patient's sample two 3 peptide fragments correspondingly calculatem,RnIt is as follows with NC values:
Table 5
Step 104:The extraction of neoantigen associated eigenvalue:For the MHC-I combination neoantigens of prediction in step (2), carry All associated eigenvalues of tumour neoantigen are taken out, the mel_38 of patient's sample two neoantigen and characteristic value are shown in Table 4, table 5;
Step 105:The setting of neoantigen activity scoring functions:For the neoantigen characteristic value of extraction in step (3), setting Neoantigen activity scoring functions;
Step 106:Neoantigen sequence based on neoantigen activity scoring functions:By neoantigen activity scoring functions to new Antigen is ranked up, and table 6 gives the marking (Neo_ of the neoantigen through the checking of PMHC activity experiments in the mel_38 of patient's sample two Score the sequence (Rank)) and in complete or collected works
Table 6
In the neoantigen that 3 identify, FLYNLLTRVY is to stimulate just to show before DC cells using tumor vaccine The immunocompetence of CD8+T cells can be activated, remaining 2 are then after using vaccine enhancing immune system ability, possess different journeys The immunocompetence of degree.We have seen that QLSCISTYV and FLYNLLTRVY is in Top 20 (total candidate's neoantigen number 117).And KLMNIQQKL either in vitro experiment or human tumor microenvironment in, because expression quantity is relatively low, determine its immunogenicity and Immune response is on the weak side, is in 66 in our ranking results, experimental result and our prediction sequence quite it is identical.
In summary, it is proposed that neoantigen immunocompetence scoring functions, can be exempted from effectively measurement neoantigen Epidemic disease activity, help is provided with immunization therapy for clinical trial and tumor research.

Claims (4)

1. the prediction of neoantigen immunocompetence and sort method based on tumour neoantigen characteristic value, it is characterised in that including following Step:
(1), the input of WGS/WES, RNA-SEQ sequencing data of tumour-normal sample:Input the full base of tumour-normal sample Because of sequencing data WGS or full sequencing of extron group data WES, transcript profile sequencing data RNA-SEQ;
(2), the prediction of tumour somatic mutation and annotation, the calculating of associated eigenvalue:Sequencing number based on step (1) input According to calling Varscan or Mutect tool analysis calculates tumour somatic mutation, calls VEP (Variant Effect Prediction) instrument complete mutation annotation, call PyClone, Kallisto, Varscan or Mutect instrument calculate as Lower eigenvalue:Mutator clone ratio, mutator expression value TPM, allelic mutation frequency VAF;
(3) prediction of MHC-I combinations neoantigen, the calculating of associated eigenvalue based on tumour somatic mutation:Based in step (2) Tumour somatic mutation and annotation data, call NetMHCpan, Netchop, OptiType instrument prediction MHC-I to combine new Antigen, and calculate such as lower eigenvalue:Peptide fragment is mutated to arrange with MHC affinity sequence percentage, unmutated peptide fragment and MHC affinity Sequence percentage, peptide fragment shearing present efficiency;
(4) extraction of all associated eigenvalues of neoantigen:For the MHC-I combination neoantigens of prediction in step (3), extract swollen All associated eigenvalues of knurl neoantigen;
(5) setting of neoantigen activity scoring functions:For the neoantigen characteristic value of extraction in step (4), setting neoantigen is lived Property scoring functions;
(6) the neoantigen sequence based on neoantigen activity scoring functions:Neoantigen is carried out by neoantigen activity scoring functions Sequence.
2. the method according to claim 11, it is characterized in that:
In step (4), neoantigen associated eigenvalue includes Rm、A、Rn, E, NC, CL, wherein:
Rm- mutation peptide fragment and MHC affinity sequence percentage, are calculated by NetMHCpan;
A-allelic mutation frequency VAF, is calculated by Varscan/Mutect/Strelka2;
Rn- unmutated peptide fragment and MHC affinity sequence percentage, are calculated by NetMHCpan;
E-mutator expression value TPM, is calculated by Kallisto;
NC-peptide fragment shearing presents efficiency, is calculated by netchop;
CL-mutator clone's ratio, is calculated by pyclone.
3. the method according to claim 11, it is characterized in that:
In step (5), the neoantigen Activity Prediction scoring functions of proposition are:
Neo_Score=abundancedissimilarityclonality;
Clonality=NCCL;
Abundance=L (Rm) Atanh (E/k);
Dissimilarity=(1-L (Rn)/2));
Wherein:L (x)=1/ (1+e5(x-2)), tanh (x) is hyperbolic tangent function;
K is transcript gene expression abundance threshold value, default value 1.
4. the method according to claim 11, it is characterized in that:In step (6), newly resisted by neoantigen Activity Prediction function pair It is as follows that original is ranked up algorithmic procedure:
A), for the MHC-I combination neoantigens of all predictions, neoantigen Activity Prediction function Neo_score is called to calculate newly The predicted value of antigen active;
B), the predicted value based on neoantigen activity, is ranked up using quick sorting algorithm to neoantigen;
C), neoantigen ranking results are exported.
CN201711071334.8A 2017-11-03 2017-11-03 Neoantigen activity prediction and sequencing method based on tumor neoantigen characteristic value Active CN107704727B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711071334.8A CN107704727B (en) 2017-11-03 2017-11-03 Neoantigen activity prediction and sequencing method based on tumor neoantigen characteristic value

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711071334.8A CN107704727B (en) 2017-11-03 2017-11-03 Neoantigen activity prediction and sequencing method based on tumor neoantigen characteristic value

Publications (2)

Publication Number Publication Date
CN107704727A true CN107704727A (en) 2018-02-16
CN107704727B CN107704727B (en) 2020-01-31

Family

ID=61177577

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711071334.8A Active CN107704727B (en) 2017-11-03 2017-11-03 Neoantigen activity prediction and sequencing method based on tumor neoantigen characteristic value

Country Status (1)

Country Link
CN (1) CN107704727B (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108796055A (en) * 2018-06-12 2018-11-13 深圳裕策生物科技有限公司 Tumor neogenetic antigen detection method, device and storage medium based on the sequencing of two generations
CN109021062A (en) * 2018-08-06 2018-12-18 倍而达药业(苏州)有限公司 A kind of screening technique of tumour neoantigen
WO2018183980A3 (en) * 2017-03-31 2018-12-20 Pei Jia Yang Ranking system for immunogenic cancer-specific epitopes
CN109880894A (en) * 2019-03-05 2019-06-14 杭州西合森医学检验实验室有限公司 The construction method of tumour immunity microenvironment prediction model based on RNAseq
WO2019206725A1 (en) * 2018-04-25 2019-10-31 Koninklijke Philips N.V. Tumor functional mutation and epitope loads as improved predictive biomarkers for immunotherapy response
CN110706742A (en) * 2019-09-30 2020-01-17 中生康元生物科技(北京)有限公司 Pan-cancer tumor neoantigen high-throughput prediction method and application thereof
CN111415707A (en) * 2020-03-10 2020-07-14 四川大学 Prediction method of clinical individualized tumor neoantigen
WO2020187143A1 (en) * 2019-03-15 2020-09-24 痕准生物科技有限公司 Method for identifying neoantigens
CN113762416A (en) * 2021-10-15 2021-12-07 南京澄实生物科技有限公司 Antigen immunogenicity prediction method and system based on multi-mode depth coding

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11118792A (en) * 1997-10-17 1999-04-30 Ajinomoto Co Inc Method for sorting macrophage, assay for immunological disease and screening method of remedy for immunological disease
CN103257236A (en) * 2012-02-20 2013-08-21 南京市疾病预防控制中心 Application of specific CTL cell epitope peptide of HLA-A24 restrictive Mycobacterium tuberculosis

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11118792A (en) * 1997-10-17 1999-04-30 Ajinomoto Co Inc Method for sorting macrophage, assay for immunological disease and screening method of remedy for immunological disease
CN103257236A (en) * 2012-02-20 2013-08-21 南京市疾病预防控制中心 Application of specific CTL cell epitope peptide of HLA-A24 restrictive Mycobacterium tuberculosis

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵恢准: "构建小鼠肉瘤随机表达文库筛选肿瘤新抗原", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018183980A3 (en) * 2017-03-31 2018-12-20 Pei Jia Yang Ranking system for immunogenic cancer-specific epitopes
US11485784B2 (en) 2017-03-31 2022-11-01 Act Genomics (Ip) Co., Ltd. Ranking system for immunogenic cancer-specific epitopes
CN112292464A (en) * 2018-04-25 2021-01-29 皇家飞利浦有限公司 Tumor function mutations and epitope burden as improved predictive biomarkers of immunotherapy response
WO2019206725A1 (en) * 2018-04-25 2019-10-31 Koninklijke Philips N.V. Tumor functional mutation and epitope loads as improved predictive biomarkers for immunotherapy response
CN108796055B (en) * 2018-06-12 2022-04-08 深圳裕策生物科技有限公司 Method, device and storage medium for detecting tumor neoantigen based on second-generation sequencing
CN108796055A (en) * 2018-06-12 2018-11-13 深圳裕策生物科技有限公司 Tumor neogenetic antigen detection method, device and storage medium based on the sequencing of two generations
CN109021062A (en) * 2018-08-06 2018-12-18 倍而达药业(苏州)有限公司 A kind of screening technique of tumour neoantigen
CN109880894A (en) * 2019-03-05 2019-06-14 杭州西合森医学检验实验室有限公司 The construction method of tumour immunity microenvironment prediction model based on RNAseq
CN113316818A (en) * 2019-03-15 2021-08-27 痕准生物科技有限公司 Method for identifying neoantigens
WO2020187143A1 (en) * 2019-03-15 2020-09-24 痕准生物科技有限公司 Method for identifying neoantigens
CN113316818B (en) * 2019-03-15 2024-04-02 痕准生物医学技术(厦门)有限公司 Method for identifying neoantigen
CN110706742A (en) * 2019-09-30 2020-01-17 中生康元生物科技(北京)有限公司 Pan-cancer tumor neoantigen high-throughput prediction method and application thereof
CN111415707A (en) * 2020-03-10 2020-07-14 四川大学 Prediction method of clinical individualized tumor neoantigen
CN113762416A (en) * 2021-10-15 2021-12-07 南京澄实生物科技有限公司 Antigen immunogenicity prediction method and system based on multi-mode depth coding

Also Published As

Publication number Publication date
CN107704727B (en) 2020-01-31

Similar Documents

Publication Publication Date Title
CN107704727A (en) Neoantigen Activity Prediction and sort method based on tumour neoantigen characteristic value
CN108388773B (en) A kind of identification method of tumor neogenetic antigen
CN108796055B (en) Method, device and storage medium for detecting tumor neoantigen based on second-generation sequencing
Gameiro et al. Treatment-naïve HPV+ head and neck cancers display a T-cell-inflamed phenotype distinct from their HPV-counterparts that has implications for immunotherapy
KR102011440B1 (en) Systems and Methods for High-Accuracy Variation Determination
US20220122690A1 (en) Attention-based neural network to predict peptide binding, presentation, and immunogenicity
CN111415707B (en) Prediction method of clinical individuation tumor neoantigen
US7756644B2 (en) System and method for automated selection of T-cell epitopes
CN110706742B (en) Pan-cancer tumor neoantigen high-throughput prediction method and application thereof
CN110752041A (en) Method, device and storage medium for predicting neoantigen based on next generation sequencing
KR102278586B1 (en) Method and System for Screening Neoantigens, and Use thereof
JP6710004B2 (en) Monitoring or diagnosis for immunotherapy and design of therapeutic agents
CN110799196A (en) System for ranking immunogenic cancer-specific epitopes
CN114446389B (en) Tumor neoantigen feature analysis and immunogenicity prediction tool and application thereof
CN110464840A (en) A kind of preparation method of tumor vaccine and the tumor vaccine prepared using this method
CN115747327A (en) Novel antigen prediction methods involving frameshift mutations
CN114882951A (en) Method and device for detecting MHC II tumor neoantigen based on next generation sequencing data
Finnigan Jr et al. Mutation-derived tumor antigens: novel targets in cancer immunotherapy
CN115104156A (en) Methods and systems for optimizing vaccine design
CN112210596B (en) Tumor neoantigen prediction method based on gene fusion event and application thereof
Xu et al. A highly effective system for predicting MHC-II epitopes with immunogenicity
US20210033608A1 (en) Methods and Systems for Identification of Human Leukocyte Antigen Peptide Presentation and Applications Thereof
Hao et al. Improvement of neoantigen identification through convolution neural network
CN116580771A (en) Method and device for predicting tumor neoantigen
Zapata et al. dN/dS dynamics quantify tumour immunogenicity and predict response to immunotherapy

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
TR01 Transfer of patent right

Effective date of registration: 20210926

Address after: 310000 room 2-201-6, xinwan science and Technology Innovation Park, xinwan street, Qiantang new area, Hangzhou, Zhejiang Province

Patentee after: Hangzhou zhuotong Internet Co.,Ltd.

Address before: 310019 room 301-3, floor 2, building 3, No. 9, Jiuhuan Road, Jianggan District, Hangzhou City, Zhejiang Province

Patentee before: HANGZHOU FENGQI INTELLIGENT TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20221103

Address after: 231261 Hongxing Villager Group, Fengxiang Village, Guanting Town, Feixi County, Hefei City, Anhui Province

Patentee after: Wang Yong

Address before: 310000 room 2-201-6, xinwan science and Technology Innovation Park, xinwan street, Qiantang new area, Hangzhou, Zhejiang Province

Patentee before: Hangzhou zhuotong Internet Co.,Ltd.

TR01 Transfer of patent right