CN113373234A - Small cell lung cancer molecular typing determination method based on mutation characteristics and application - Google Patents
Small cell lung cancer molecular typing determination method based on mutation characteristics and application Download PDFInfo
- Publication number
- CN113373234A CN113373234A CN202110770114.4A CN202110770114A CN113373234A CN 113373234 A CN113373234 A CN 113373234A CN 202110770114 A CN202110770114 A CN 202110770114A CN 113373234 A CN113373234 A CN 113373234A
- Authority
- CN
- China
- Prior art keywords
- lung cancer
- small cell
- cell lung
- molecular typing
- typing
- 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.)
- Pending
Links
- 206010041067 Small cell lung cancer Diseases 0.000 title claims abstract description 76
- 208000000587 small cell lung carcinoma Diseases 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000035772 mutation Effects 0.000 title claims abstract description 23
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 31
- 238000012163 sequencing technique Methods 0.000 claims abstract description 23
- 230000008859 change Effects 0.000 claims abstract description 11
- 108020004414 DNA Proteins 0.000 claims description 19
- 206010028980 Neoplasm Diseases 0.000 claims description 16
- 201000011510 cancer Diseases 0.000 claims description 12
- 210000001082 somatic cell Anatomy 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000010186 staining Methods 0.000 claims description 11
- 238000012217 deletion Methods 0.000 claims description 10
- 230000037430 deletion Effects 0.000 claims description 10
- 238000005406 washing Methods 0.000 claims description 10
- 230000003321 amplification Effects 0.000 claims description 9
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 9
- WZUVPPKBWHMQCE-UHFFFAOYSA-N Haematoxylin Chemical compound C12=CC(O)=C(O)C=C2CC2(O)C1C1=CC=C(O)C(O)=C1OC2 WZUVPPKBWHMQCE-UHFFFAOYSA-N 0.000 claims description 8
- 201000010099 disease Diseases 0.000 claims description 8
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 8
- 238000011532 immunohistochemical staining Methods 0.000 claims description 8
- 210000001519 tissue Anatomy 0.000 claims description 8
- 238000010219 correlation analysis Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 7
- 108700019961 Neoplasm Genes Proteins 0.000 claims description 6
- 102000048850 Neoplasm Genes Human genes 0.000 claims description 6
- 210000000349 chromosome Anatomy 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 5
- BJTZZQBKPAGZNC-UHFFFAOYSA-N 5,5-diaminocyclohexa-1,3-diene-1-carboxylic acid Chemical compound NC1(CC(C(=O)O)=CC=C1)N BJTZZQBKPAGZNC-UHFFFAOYSA-N 0.000 claims description 4
- 108700028369 Alleles Proteins 0.000 claims description 4
- 108091023037 Aptamer Proteins 0.000 claims description 4
- 108091026890 Coding region Proteins 0.000 claims description 4
- 239000012153 distilled water Substances 0.000 claims description 4
- 238000004043 dyeing Methods 0.000 claims description 4
- 239000012634 fragment Substances 0.000 claims description 4
- 238000003780 insertion Methods 0.000 claims description 4
- 230000037431 insertion Effects 0.000 claims description 4
- 230000003902 lesion Effects 0.000 claims description 4
- 238000007854 ligation-mediated PCR Methods 0.000 claims description 4
- 230000000869 mutational effect Effects 0.000 claims description 4
- 230000007935 neutral effect Effects 0.000 claims description 4
- 230000001575 pathological effect Effects 0.000 claims description 4
- 229920000642 polymer Polymers 0.000 claims description 4
- 230000008439 repair process Effects 0.000 claims description 4
- 230000000717 retained effect Effects 0.000 claims description 4
- 238000002864 sequence alignment Methods 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 238000007482 whole exome sequencing Methods 0.000 claims description 4
- 101100298998 Caenorhabditis elegans pbs-3 gene Proteins 0.000 claims description 3
- 108090000790 Enzymes Proteins 0.000 claims description 3
- 102000004190 Enzymes Human genes 0.000 claims description 3
- 239000003814 drug Substances 0.000 claims description 3
- 238000012239 gene modification Methods 0.000 claims description 3
- 230000033607 mismatch repair Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000000391 smoking effect Effects 0.000 claims description 3
- 210000004027 cell Anatomy 0.000 claims description 2
- 229910019142 PO4 Inorganic materials 0.000 claims 1
- NBIIXXVUZAFLBC-UHFFFAOYSA-K phosphate Chemical compound [O-]P([O-])([O-])=O NBIIXXVUZAFLBC-UHFFFAOYSA-K 0.000 claims 1
- 239000010452 phosphate Substances 0.000 claims 1
- 206010064571 Gene mutation Diseases 0.000 abstract description 5
- 101001117317 Homo sapiens Programmed cell death 1 ligand 1 Proteins 0.000 abstract description 5
- 102100024216 Programmed cell death 1 ligand 1 Human genes 0.000 abstract description 5
- 210000001266 CD8-positive T-lymphocyte Anatomy 0.000 abstract description 4
- 238000011161 development Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 abstract description 2
- 238000011337 individualized treatment Methods 0.000 abstract description 2
- 238000011282 treatment Methods 0.000 description 5
- 230000010558 Gene Alterations Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 239000008363 phosphate buffer Substances 0.000 description 3
- 108010078814 Tumor Suppressor Protein p53 Proteins 0.000 description 2
- 102000015098 Tumor Suppressor Protein p53 Human genes 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002512 chemotherapy Methods 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 101150079831 OR51I2 gene Proteins 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 230000002759 chromosomal effect Effects 0.000 description 1
- 230000018044 dehydration Effects 0.000 description 1
- 238000006297 dehydration reaction Methods 0.000 description 1
- 239000003596 drug target Substances 0.000 description 1
- 238000009169 immunotherapy Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000036210 malignancy Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000036438 mutation frequency Effects 0.000 description 1
- 239000002773 nucleotide Substances 0.000 description 1
- 125000003729 nucleotide group Chemical group 0.000 description 1
- 238000010837 poor prognosis Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
- C12Q1/6886—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material for cancer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/53—Immunoassay; Biospecific binding assay; Materials therefor
- G01N33/574—Immunoassay; Biospecific binding assay; Materials therefor for cancer
- G01N33/57407—Specifically defined cancers
- G01N33/57423—Specifically defined cancers of lung
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/53—Immunoassay; Biospecific binding assay; Materials therefor
- G01N33/574—Immunoassay; Biospecific binding assay; Materials therefor for cancer
- G01N33/57484—Immunoassay; Biospecific binding assay; Materials therefor for cancer involving compounds serving as markers for tumor, cancer, neoplasia, e.g. cellular determinants, receptors, heat shock/stress proteins, A-protein, oligosaccharides, metabolites
- G01N33/57492—Immunoassay; Biospecific binding assay; Materials therefor for cancer involving compounds serving as markers for tumor, cancer, neoplasia, e.g. cellular determinants, receptors, heat shock/stress proteins, A-protein, oligosaccharides, metabolites involving compounds localized on the membrane of tumor or cancer cells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/68—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides 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
-
- 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
- G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
- G16B50/10—Ontologies; Annotations
-
- 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
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/112—Disease subtyping, staging or classification
-
- 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
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/156—Polymorphic or mutational markers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2333/00—Assays involving biological materials from specific organisms or of a specific nature
- G01N2333/435—Assays involving biological materials from specific organisms or of a specific nature from animals; from humans
- G01N2333/705—Assays involving receptors, cell surface antigens or cell surface determinants
- G01N2333/70503—Immunoglobulin superfamily, e.g. VCAMs, PECAM, LFA-3
- G01N2333/70517—CD8
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2333/00—Assays involving biological materials from specific organisms or of a specific nature
- G01N2333/435—Assays involving biological materials from specific organisms or of a specific nature from animals; from humans
- G01N2333/705—Assays involving receptors, cell surface antigens or cell surface determinants
- G01N2333/70503—Immunoglobulin superfamily, e.g. VCAMs, PECAM, LFA-3
- G01N2333/70532—B7 molecules, e.g. CD80, CD86
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Urology & Nephrology (AREA)
- Hematology (AREA)
- General Health & Medical Sciences (AREA)
- Biotechnology (AREA)
- Biomedical Technology (AREA)
- Analytical Chemistry (AREA)
- Cell Biology (AREA)
- Pathology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Biochemistry (AREA)
- Microbiology (AREA)
- Oncology (AREA)
- Hospice & Palliative Care (AREA)
- Food Science & Technology (AREA)
- General Physics & Mathematics (AREA)
- Organic Chemistry (AREA)
- Biophysics (AREA)
- Medicinal Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Genetics & Genomics (AREA)
- Evolutionary Biology (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Wood Science & Technology (AREA)
- Zoology (AREA)
- Databases & Information Systems (AREA)
- Bioethics (AREA)
- General Engineering & Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
The invention belongs to the technical field of genes, and discloses a small cell lung cancer molecule typing determination method based on mutation characteristics and application thereof. The invention comprehensively analyzes the characteristics of gene mutation, copy number change and the like of 178 cases of small cell lung cancer by a Whole Exon Sequencing (WES) method, and analyzes the immune characteristics of TMB, TNB, PDL1, CD8+ T cells and the like. The invention integrates the gene characteristics and the immune characteristics of the small cell lung cancer, and provides three gene change molecular typing of the small cell lung cancer for the first time. The invention is essential for understanding the characteristics of heterogeneity of small cell lung cancer and developing individualized treatment, and can promote the development of individualized clinical experiments of small cell lung cancer.
Description
Technical Field
The invention belongs to the technical field of genes, and particularly relates to a small cell lung cancer molecular typing determination method based on mutation characteristics and application.
Background
At present, Small Cell Lung Cancer (SCLC) has high malignancy and poor prognosis, and chemotherapy or chemotherapy combined radiotherapy is always applied clinically as two important treatment means. Recently, although immunotherapy has made some research progress in small cell lung cancer, the progress is limited. Due to the unique biological characteristics of the small cell lung cancer, the genome characteristics are determined and analyzed, and the method has very important clinical significance. In fact, molecular biology studies have shown that small cell lung cancer can be divided into four subtypes, SCLC-A, SCLC-N, SCLC-P and SCLC-Y, which can provide necessary clues for personalized treatment. However, current molecular typing is derived by analyzing the transcriptome level of foreign data.
At present, relevant molecular typing is carried out aiming at genome mutation characteristics and immune related characteristics of the small cell lung cancer, and no relevant literature report exists at home and abroad. Therefore, the genome alteration characteristics and the immune characteristics of the small cell lung cancer are deeply and comprehensively analyzed, and related subtype classification is carried out, so that the individual treatment of the small cell lung cancer is further promoted.
Through the above analysis, the problems and defects of the prior art are as follows: the prior art has no technology for carrying out subtype typing analysis and determination on the genome characteristics and the immune characteristics of the small cell lung cancer.
The difficulty in solving the above problems and defects is: since TP53 and RB1 are the two most common mutations in small cell lung cancer, it is difficult to re-molecularly type based on gene alteration profiles.
The significance of solving the problems and the defects is as follows: the method is very important for individual treatment of small cell lung cancer and excavation of drug targets.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a small cell lung cancer molecular typing determination method based on mutation characteristics and application thereof.
The invention is realized in such a way that the application of the small cell lung cancer gene modification molecular typing in the preparation of the drugs for treating the small cell lung cancer is provided.
The small cell lung cancer gene alteration molecular typing is as follows: cluster 1 (mainly enriched in Signature 4: smoking-related gene characteristics (leading to the conversion of G.C to T.A, with a tendency to methylate CpG dinucleotides)), Cluster2 (mainly enriched in Signature 6: DNA mismatch repair-related gene characteristics (additional C.G to T.A at NpCpG base sites)), and Cluster 3 (mainly enriched in Signature 5: unidentified related gene characteristics (several characteristics of A.T to G.C)).
Another object of the present invention is to provide a kit, andthe above-mentionedThe kit comprises the small cell lung cancer gene alteration molecular typing in the application.
Another object of the present invention is to provide a method for determining molecular typing of small cell lung cancer based on mutation characteristics, comprising:
and (3) carrying out signature correlation analysis on the genome by analyzing the gene change characteristics of the small cell lung cancer, and carrying out unsupervised clustering to determine the typing of the small cell lung cancer molecules.
Further, the small cell lung cancer molecular typing determination method based on the mutation characteristics comprises the following steps:
step one, constructing a DNA library and sequencing a whole exome; carrying out sequence comparison and variation detection;
step two, performing immunohistochemical staining by adopting an enhanced marking polymer system; analyzing the gene change characteristics of the small cell lung cancer, carrying out signature correlation analysis on the genome, carrying out unsupervised clustering, and determining the typing result of the small cell lung cancer molecules.
Further, in step one, the performing DNA library construction and whole exome sequencing comprises:
cutting genome DNA into 200bp fragments by enzyme method; after end repair and A tailing, adding T-shaped joints at both ends; amplifying the purified DNA by ligation-mediated PCR, sequencing the paired-end multiplex samples using NovaSeq6000 system; sequencing libraries were generated using 96RXN xGen outer Research Panel v 1.0.
Further, the sequencing comprises: the sequencing depth was 200 x per tissue sample and 100 x per control.
Further, in the first step, the sequence alignment and variation detection comprises:
(1) performing FASTP pretreatment on the raw data to obtain an aptamer sequence, and comparing clean readings in Fast Q format with a reference human genome hg19/GRCh37 by using Burrow-Wheeler aligner; sequencing the mapped BAM file by using an SAM tool and Picard and processing a PCR copy;
(2) generating a final BAM file by using the GATK, detecting the mononucleotide mutation of the somatic cell by using the MuTect, and detecting the small insertion and deletion of the somatic cell by using a GATK somatic cell Indel detector; performing variant annotation based on multiple databases using ANNOVAR;
(3) further analyzing the cancer hotspots screened from the disease database by the retained non-synonymous SNVs with the variation allele frequency (namely cut off is more than or equal to 3% or VAF cut off is more than or equal to 1%);
(4) calculating tumor mutation load, and calculating the number of nonsynonymous SNV and the number of included variation of each coding region; identifying important driving genes by combining two methods of MutsigCV and dDNdScv;
(5) the genome of the cancer important target is used for recognizing 2.0 copy number variation.
Further, the plurality of databases includes a variant description, a population frequency database, a variant function prediction database, and a phenotype or disease database.
Further, the identification of 2.0 copy number variation by using cancer important target genome comprises: at the chromosome arm level, meaningful amplifications or deletions were screened for further analysis with FDR based on criteria with cutoff < 10%; at the lesion level, significant amplifications were screened with FDR based on criteria with a cutoff of < 5% and G-Score based on cutoff >0.3, and further analysis was performed with FDR (criteria with a cutoff of < 5% and G-Score based on cutoff of < -0.2).
Further, the performing immunohistochemical staining comprises:
with anti-PD-L1, anti-CD 8+After overnight staining at 4 ℃, washing with phosphate buffer for 3 times, each for 5 min;
staining with corresponding secondary antibody at 37 deg.C for 30min, washing with PBS 3 times for 5min each time;
further dyeing with 3, 3-diaminobenzoic acid, washing with distilled water;
performing hematoxylin staining, dehydrating, removing, and fixing with neutral gum;
and collecting stained tissue images by using a KF-PRO-120 and KF-BIO digital pathological slide scanner.
Further, the anti-PD-L1 is Cst, 13684, 1: 100, respectively; the anti-CD 8+Cst, 85336, 1: 100.
by combining all the technical schemes, the invention has the advantages and positive effects that: the invention comprehensively analyzes the characteristics of gene mutation, copy number change and the like of 178 cases of small cell lung cancer by a Whole Exon Sequencing (WES) method, and analyzes the immune characteristics of TMB, TNB, PDL1, CD8+ T cells and the like. The invention integrates the gene characteristics and the immune characteristics of the small cell lung cancer, and provides three gene change molecular typing of the small cell lung cancer for the first time. The invention is essential for solving the characteristics of small cell lung cancer heterogeneity and developing individualized treatment, and can promote the development of individualized clinical experiments of small cell lung cancer.
Drawings
FIG. 1 is a flow chart of a method for determining molecular typing of small cell lung cancer based on mutation characteristics according to an embodiment of the present invention.
FIG. 2 is a schematic diagram of the genomic molecular typing of small cell lung cancer provided by the embodiments of the present invention.
FIG. 3 is a schematic diagram of the difference of the molecular typing gene mutation of the small cell lung cancer provided by the embodiment of the invention.
FIG. 4 is a schematic diagram of the difference between the copy numbers of the typing genes of the small cell lung cancer molecules according to the embodiment of the present invention.
FIG. 5 is a schematic diagram of the differences of the molecular typing immune-related characteristics of small cell lung cancer provided by the embodiments of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Aiming at the problems in the prior art, the invention provides a small cell lung cancer molecular typing determination method based on mutation characteristics and application thereof, and the invention is described in detail below with reference to the accompanying drawings.
The small cell lung cancer molecular typing determination method based on the mutation characteristics provided by the embodiment of the invention comprises the following steps:
and (3) carrying out signature correlation analysis on the genome by analyzing the gene change characteristics of the small cell lung cancer, and carrying out unsupervised clustering to determine the typing of the small cell lung cancer molecules.
The invention provides an application of small cell lung cancer gene modification molecule typing in preparation of a small cell lung cancer treatment drug.
The small cell lung cancer gene alteration molecular typing is as follows:
cluster 1 (mainly enriched Signature 4: smoking-related gene characteristics (leading to the conversion of G.C to T.A, with a tendency to methylate CpG dinucleotides)), Cluster2 (mainly enriched Signature 6: DNA mismatch repair-related gene characteristics (extra C.G to T.A at NpCpG base sites)) and Cluster 3 (mainly enriched Signature 5: unidentified related gene characteristics (several characteristics of A.T to G.C)).
It is another object of the present invention to provide a kit comprising the small cell lung cancer gene-altering molecular typing for use as described.
As shown in fig. 1, the method for determining the molecular typing of small cell lung cancer based on the mutation characteristics provided by the embodiment of the present invention includes the following steps:
s101, constructing a DNA library and sequencing a whole exome; carrying out sequence comparison and variation detection;
s102, performing immunohistochemical staining by using an enhanced marker polymer system;
s103, analyzing the gene change characteristics of the small cell lung cancer, carrying out signature correlation analysis on the genome, carrying out unsupervised clustering, and determining the typing result of the small cell lung cancer molecules.
The DNA library construction and whole exome sequencing provided by the embodiment of the invention comprises the following steps:
cutting genome DNA into 200bp fragments by enzyme method; after end repair and A tailing, adding T-shaped joints at both ends; amplifying the purified DNA by ligation-mediated PCR, sequencing the paired-end multiplex samples using NovaSeq6000 system; sequencing libraries were generated using 96RXN xGen outer Research Panel v 1.0.
The sequencing provided by the embodiment of the invention comprises the following steps: the sequencing depth was 200 x per tissue sample and 100 x per control.
The sequence alignment and variation detection provided by the embodiment of the invention comprises the following steps:
(1) performing FASTP pretreatment on the raw data to obtain an aptamer sequence, and comparing clean readings in Fast Q format with a reference human genome hg19/GRCh37 by using Burrow-Wheeler aligner; sequencing the mapped BAM file by using an SAM tool and Picard and processing a PCR copy;
(2) generating a final BAM file by using the GATK, detecting the mononucleotide mutation of the somatic cell by using the MuTect, and detecting the small insertion and deletion of the somatic cell by using a GATK somatic cell Indel detector; performing variant annotation based on multiple databases using ANNOVAR;
(3) further analyzing the cancer hotspots screened from the disease database by the retained non-synonymous SNVs with the variation allele frequency (namely cut off is more than or equal to 3% or VAF cut off is more than or equal to 1%);
(4) calculating tumor mutation load, and calculating the number of nonsynonymous SNV and the number of included variation of each coding region; identifying important driving genes by combining two methods of MutsigCV and dDNdScv;
(5) the genome of the cancer important target is used for recognizing 2.0 copy number variation.
The plurality of databases provided by the embodiments of the present invention include a variant description, a population frequency database, a variant function prediction database, and a phenotype or disease database.
The method for identifying 2.0 copy number variation by using the cancer important target genome provided by the embodiment of the invention comprises the following steps: at the chromosome arm level, meaningful amplifications or deletions were screened for further analysis with FDR based on criteria with cutoff < 10%; at the lesion level, significant amplifications were screened with FDR based on criteria with a cutoff of < 5% and G-Score based on cutoff >0.3, and further analysis was performed with FDR (criteria with a cutoff of < 5% and G-Score based on cutoff of < -0.2).
The immunohistochemical staining method provided by the embodiment of the invention comprises the following steps:
with CST, 13684, 1: 100 anti-PD-L1, CST, 85336, 1: 100 anti-CD 8+After overnight staining at 4 ℃, washing with phosphate buffer for 3 times, each for 5 min;
staining with corresponding secondary antibody at 37 deg.C for 30min, washing with PBS 3 times for 5min each time;
further dyeing with 3, 3-diaminobenzoic acid, washing with distilled water;
performing hematoxylin staining, dehydrating, removing, and fixing with neutral gum;
and collecting stained tissue images by using a KF-PRO-120 and KF-BIO digital pathological slide scanner.
The anti-PD-L1 provided by the embodiment of the invention is CST, 13684, 1: 100, respectively; the anti-CD 8+CST, 85336, 1: 100.
the typing result of the small cell lung cancer molecule provided by the embodiment of the invention comprises the following steps:
obtaining 3 molecular typing subclasses which are respectively Cluster 1, Cluster2 and Cluster 3 by carrying out unsupervised clustering.
The technical solution of the present invention is further described with reference to the following specific embodiments.
Example 1:
DNA library construction and Whole Exome Sequencing (WES).
The genomic DNA was enzymatically cleaved into 200 bp-sized fragments. After end repair and a tailing, a T-joint was added at both ends. To construct the library, the purified DNA was amplified by ligation-mediated PCR, followed by generation of the final sequencing library using 96RXN xGen outer Research Panel v1.0(Integrated DNA Technologies, USA) according to the manufacturer's instructions. Paired-end multiplex samples were sequenced using the NovaSeq6000 system. The sequencing depth was 200 x per tissue sample and 100 x per control.
And (3) sequence alignment and variation detection.
The raw data were pre-processed by FASTP to trim the aptamer sequences, and clean reads in Fast Q format were subsequently aligned with the reference human genome (hg19/GRCh37) using a Burrow-Wheeler aligner (BWA, v0.7.15). The SAM tool and Picard (2.12.1) (http:// Picard. sourceform. net) are used to sort the mapped BAM files and process the PCR copies. To calculate sequencing coverage and depth, final BAM files were generated from GATK (genome analysis toolkit 3.8) for local alignment and base quality recalibration. Single Nucleotide Variants (SNV) were detected from somatic cells using MuTect, and small insertions and deletions (Indels) were detected from somatic cells using the GATK somatic cell Indel detector. ANNOVAR software was used for variant annotation based on multiple databases, including variant descriptions (HGV), population frequency databases (1000 genome project, dbSNP, ExAC), variant function prediction databases (Polyphen-2 and SIFT), and phenotypic or disease databases (OMIM, COSMIC, ClinVar). The non-synonymous SNVs retained were further analyzed for screening of cancer hotspots from disease databases using Variant Allele Frequency (VAF) (cut off 3% or more) or VAF (cut off 1% or more). Tumor Mutational Burden (TMB) was calculated, and the number of nonsynonymous SNVs and the number of inclusion variations per coding region were calculated. Important driver genes were identified by combining both MutsigCV and dNdScv methods, as described in previous studies, with false discovery rates (FDR, cut off < 10%). Copy Number Variation (CNV) was identified using the cancer important target genome identification (gist) 2.0 algorithm. At the chromosomal arm level, meaningful amplifications or deletions were screened with FDR (cutoff < 10%) for further analysis. At the lesion level, significant amplifications were screened with FDR (cut off < 5%) and G-Score (cut off >0.3) and significant deletions were screened with FDR (cut off < 5%) and G-Score (cut off < -0.2) for further analysis.
Immunohistochemical staining.
Immunohistochemical staining was performed using an Enhanced Labeling Polymer System (ELPS). After staining with anti-PD-L1 (Cst, 13684, 1: 100), anti-CD 8+ (Cst, 85336, 1: 100) overnight at 4 ℃ washing 3 times with phosphate buffer each for 5 min. The cells were stained with the corresponding secondary antibody at 37 ℃ for 30min and then washed 3 times with PBS for 5min each. Further dyeing with 3, 3-Diaminobenzoic Acid (DAB) and washing with distilled water. Hematoxylin staining was then performed, followed by dehydration, removal, and fixation with neutral gums. And collecting stained tissue images by using a KF-PRO-120 and KF-BIO digital pathological slide scanner.
Results
Genome molecular typing construction of small cell lung cancer
By analyzing the gene change characteristics of the small cell lung cancer, signature correlation analysis is carried out on the genome, unsupervised clustering is carried out, and finally 3 subclasses are obtained, namely Cluster 1, Cluster2 and Cluster 3.
The invention compares the difference of the driving gene mutation of the three molecular typing of Cluster 1, Cluster2 and Cluster 3. The results show that TP53 and RB1 are still the genes with the highest mutation frequency of the small cell lung cancer, and are 93.26 percent and 44.38 percent respectively. The OR51I2 gene mutation is different among three subtypes.
The invention compares the difference between the gene copy number changes of the chromosome level of the three molecular typing of Cluster 1, Cluster2 and Cluster 3. The results indicate that at the gene copy number amplification level, there are some differences in the three molecular typing on chromosomes 1, 2, 3, 5, 6, 15, 16, 20, and 22. In addition, there are some differences in the three molecular typing on chromosomes 5, 7, 9, 10, 11, 14, 15, and 16 in terms of gene copy number deletion. (. 0.05;. 0.01;. 0.001)
The invention compares the differences of the immune related characteristics of the three molecular classifications of Cluster 1, Cluster2 and Cluster 3. The results show that there are some differences among TNB, PDL1(TPS), PDL1(CPS) and CD8+ T cells in all three molecular typing. Importantly, the surprising finding of this study was that the Cluster 1 subtype has higher TNB, and has higher PDL1 expression and infiltration of CD8+ T cells.
The above description is only for the purpose of illustrating the present invention and the appended claims are not to be construed as limiting the scope of the invention, which is intended to cover all modifications, equivalents and improvements that are within the spirit and scope of the invention as defined by the appended claims.
Claims (10)
1. An application of small cell lung cancer gene modification molecule typing in preparing medicine for treating small cell lung cancer is disclosed.
2. The use of claim 1, wherein said small cell lung cancer genetically altered molecular typing comprises:
cluster 1, enriching Signature 4, smoking related gene characteristics;
cluster2, enriching Signature 6, DNA mismatch repair related gene characteristics;
cluster 3, enriching Signature 5, and gathering unknown related gene characteristics.
3. A kit comprising the small cell lung cancer genetically altered molecular typing for use according to any one of claims 1 to 2.
4. A small cell lung cancer molecular typing determination method based on mutation characteristics is characterized by comprising the following steps:
and (3) carrying out signature correlation analysis on the genome by analyzing the gene change characteristics of the small cell lung cancer, and carrying out unsupervised clustering to determine the typing of the small cell lung cancer molecules.
5. The method for determining molecular typing of small cell lung cancer based on mutation characteristics according to claim 4, wherein the method for determining molecular typing of small cell lung cancer based on mutation characteristics comprises the following steps:
step one, constructing a DNA library and sequencing a whole exome; carrying out sequence comparison and variation detection;
step two, performing immunohistochemical staining by adopting an enhanced marking polymer system; analyzing the gene change characteristics of the small cell lung cancer, carrying out signature correlation analysis on the genome, carrying out unsupervised clustering, and determining the typing result of the small cell lung cancer molecules.
6. The method for molecular typing and determination of small cell lung cancer based on mutation characteristics as claimed in claim 5, wherein in step one, the DNA library construction and whole exome sequencing comprises:
cutting genome DNA into 200bp fragments by enzyme method; after end repair and A tailing, adding T-shaped joints at both ends; amplifying the purified DNA by ligation-mediated PCR, sequencing the paired-end multiplex samples using NovaSeq6000 system; generating a sequencing library using 96RXN xGen outer research & reservoir 1.0;
the sequencing comprises the following steps: the sequencing depth was 200 x per tissue sample and 100 x per control.
7. The method for molecular typing of small cell lung cancer according to claim 5, wherein the sequence alignment and variation detection in step one comprises:
(1) performing FASTP pretreatment on the raw data to obtain an aptamer sequence, and comparing clean readings in Fast Q format with a reference human genome hg19/GRCh37 by using Burrow-Wheeler aligner; sequencing the mapped BAM file by using an SAM tool and Picard and processing a PCR copy;
(2) generating a final BAM file by using the GATK, detecting the mononucleotide mutation of the somatic cell by using the MuTect, and detecting the small insertion and deletion of the somatic cell by using a GATK somatic cell Indel detector; performing variant annotation based on multiple databases using ANNOVAR;
(3) further analyzing the cancer hotspots screened from the disease database by the retained non-synonymous SNVs with the variation allele frequency (namely cut off is more than or equal to 3% or VAF cut off is more than or equal to 1%);
(4) calculating tumor mutation load, and calculating the number of nonsynonymous SNV and the number of included variation of each coding region; identifying important driving genes by combining two methods of MutsigCV and dDNdScv;
(5) the genome of the cancer important target is used for recognizing 2.0 copy number variation.
8. The method of claim 7 for molecular typing of small cell lung cancer based on mutational characteristics wherein the plurality of databases comprise a variant description, a population frequency database, a variant function prediction database, and a phenotype or disease database;
the identification of 2.0 identification copy number variation by using cancer important target genome comprises the following steps: at the chromosome arm level, meaningful amplifications or deletions were screened for further analysis with FDR based on criteria with cutoff < 10%; at the lesion level, significant amplifications were screened with FDR based on criteria with a cutoff of < 5% and G-Score based on cutoff >0.3, and further analysis was performed with FDR (criteria with a cutoff of < 5% and G-Score based on cutoff of < -0.2).
9. The method of determining molecular typing of small cell lung cancer based on mutational characteristics according to claim 4, wherein said performing immunohistochemical staining comprises:
with anti-PD-L1, anti-CD 8+After overnight staining at 4 ℃ the cells were washed with phosphate buffer3 times, each time for 5 min;
staining with corresponding secondary antibody at 37 deg.C for 30min, washing with PBS 3 times for 5min each time;
further dyeing with 3, 3-diaminobenzoic acid, washing with distilled water;
performing hematoxylin staining, dehydrating, removing, and fixing with neutral gum;
and collecting stained tissue images by using a KF-PRO-120 and KF-BIO digital pathological slide scanner.
10. The method for molecular typing of small cell lung cancer based on mutational characteristics as claimed in claim 9, wherein the anti-PD-L1 is Cst, 13684, 1: 100, respectively; the anti-CD 8+Cst, 85336, 1: 100.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110770114.4A CN113373234A (en) | 2021-07-07 | 2021-07-07 | Small cell lung cancer molecular typing determination method based on mutation characteristics and application |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110770114.4A CN113373234A (en) | 2021-07-07 | 2021-07-07 | Small cell lung cancer molecular typing determination method based on mutation characteristics and application |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113373234A true CN113373234A (en) | 2021-09-10 |
Family
ID=77581417
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110770114.4A Pending CN113373234A (en) | 2021-07-07 | 2021-07-07 | Small cell lung cancer molecular typing determination method based on mutation characteristics and application |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113373234A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114187964A (en) * | 2021-12-13 | 2022-03-15 | 深圳市海普洛斯生物科技有限公司 | Construction method of lung cancer perioperative molecular residual lesion gene detection panel and detection model |
CN116403648A (en) * | 2023-06-06 | 2023-07-07 | 中国医学科学院肿瘤医院 | Small cell lung cancer immune novel typing method established based on multidimensional analysis |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111223525A (en) * | 2020-01-07 | 2020-06-02 | 广州基迪奥生物科技有限公司 | Tumor exon sequencing data analysis method |
WO2020136133A1 (en) * | 2018-12-23 | 2020-07-02 | F. Hoffmann-La Roche Ag | Tumor classification based on predicted tumor mutational burden |
-
2021
- 2021-07-07 CN CN202110770114.4A patent/CN113373234A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020136133A1 (en) * | 2018-12-23 | 2020-07-02 | F. Hoffmann-La Roche Ag | Tumor classification based on predicted tumor mutational burden |
CN111223525A (en) * | 2020-01-07 | 2020-06-02 | 广州基迪奥生物科技有限公司 | Tumor exon sequencing data analysis method |
Non-Patent Citations (6)
Title |
---|
CHARLES M. RUDIN ET AL.: "Small-cell lung cancer", 《NAT REV DIS PRIMERS》, pages 3 * |
COSMIC: "Mutational Signatures (v2 - March 2015)", pages 1, Retrieved from the Internet <URL:https://cancer.sanger.ac.uk/signatures/signatures_v2/> * |
HAIYONG WANG ET AL.: "Molecular subtyping of small-cell lung cancer based on mutational signatures with different genomic features and therapeutic strategies", 《CANCER SCI》, pages 665 - 679 * |
JULIE GEORGE ET AL.: "Comprehensive genomic profiles of small cell lung cancer", 《NATURE》, pages 47 * |
NAIYER A. RIZVI ET AL.: "Mutational landscape determines sensitivity to PD-1 blockade in non–small cell lung cancer", 《SCIENCE》, vol. 348, no. 6230, pages 125 * |
王安强: "混合型肝细胞癌—胆管细胞癌的克隆进化与起源", 《中国博士学位论文全文数据库 医药卫生科技辑》, pages 8 - 15 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114187964A (en) * | 2021-12-13 | 2022-03-15 | 深圳市海普洛斯生物科技有限公司 | Construction method of lung cancer perioperative molecular residual lesion gene detection panel and detection model |
CN116403648A (en) * | 2023-06-06 | 2023-07-07 | 中国医学科学院肿瘤医院 | Small cell lung cancer immune novel typing method established based on multidimensional analysis |
CN116403648B (en) * | 2023-06-06 | 2023-08-01 | 中国医学科学院肿瘤医院 | Small cell lung cancer immune novel typing method established based on multidimensional analysis |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108753967B (en) | Gene set for liver cancer detection and panel detection design method thereof | |
EP3274475B1 (en) | Alignment and variant sequencing analysis pipeline | |
KR102028375B1 (en) | Systems and methods to detect rare mutations and copy number variation | |
Wadapurkar et al. | Computational analysis of next generation sequencing data and its applications in clinical oncology | |
US20190338349A1 (en) | Methods and systems for high fidelity sequencing | |
CN104884633B (en) | By the high-throughput Genotyping that a small amount of inhereditary material is sequenced | |
KR20200093438A (en) | Method and system for determining somatic mutant clonability | |
CN110343748B (en) | Method for analyzing tumor mutation load based on high-throughput targeted sequencing | |
CN113373234A (en) | Small cell lung cancer molecular typing determination method based on mutation characteristics and application | |
Larson et al. | A clinician’s guide to bioinformatics for next-generation sequencing | |
Abbasi et al. | Significance and limitations of the use of next-generation sequencing technologies for detecting mutational signatures | |
CN116631508B (en) | Detection method for tumor specific mutation state and application thereof | |
EP4150113A1 (en) | Homologous recombination repair deficiency detection | |
Segawa et al. | HLA genotyping by next-generation sequencing of complementary DNA | |
CN115803447A (en) | Detection of structural variation in chromosome proximity experiments | |
Yadav et al. | Next-Generation sequencing transforming clinical practice and precision medicine | |
Romagnoli et al. | Resolving complex structural variants via nanopore sequencing | |
Li et al. | VarBen: generating in silico reference data sets for clinical next-generation sequencing bioinformatics pipeline evaluation | |
CN113674802B (en) | Method and device for performing variation detection based on methylation sequencing data | |
CN112695081B (en) | New susceptibility gene of primary biliary cholangitis and application thereof | |
CN111028888A (en) | Detection method of genome-wide copy number variation and application thereof | |
CN112513292A (en) | Method and device for detecting homologous sequence based on high-throughput sequencing | |
CN109979534B (en) | C site extraction method and device | |
Cheng et al. | Whole genome error-corrected sequencing for sensitive circulating tumor DNA cancer monitoring | |
EP3524695B1 (en) | Method for the enrichment of genomic regions |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210910 |
|
RJ01 | Rejection of invention patent application after publication |