CN108130372A - 一种用于急性髓系白血病药物指示的方法及装置 - Google Patents
一种用于急性髓系白血病药物指示的方法及装置 Download PDFInfo
- Publication number
- CN108130372A CN108130372A CN201810048989.1A CN201810048989A CN108130372A CN 108130372 A CN108130372 A CN 108130372A CN 201810048989 A CN201810048989 A CN 201810048989A CN 108130372 A CN108130372 A CN 108130372A
- Authority
- CN
- China
- Prior art keywords
- gene
- sample
- prognosis
- aml
- expression
- 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
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
- 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
-
- 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/106—Pharmacogenomics, i.e. genetic variability in individual responses to drugs and drug metabolism
-
- 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/118—Prognosis of disease development
-
- 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
-
- 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/158—Expression markers
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Analytical Chemistry (AREA)
- Genetics & Genomics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Physics & Mathematics (AREA)
- Organic Chemistry (AREA)
- Pathology (AREA)
- Zoology (AREA)
- Biophysics (AREA)
- Immunology (AREA)
- Biotechnology (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Wood Science & Technology (AREA)
- Microbiology (AREA)
- Biochemistry (AREA)
- Oncology (AREA)
- General Engineering & Computer Science (AREA)
- Hospice & Palliative Care (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明实施例提供了一种用于急性髓系白血病药物指示的方法及装置,方法包括:从预设的目标样本中筛选出目标评估基因,所述目标评估基因为STAB1基因;获取正常核型急性髓系白血病CN‑AML患者的目标评估基因;根据CN‑AML患者目标评估基因的基因表达量对所述CN‑AML患者进行分层,将所述CN‑AML患者分为预后良好组及预后不良组;对处于预后不良组的CN‑AML患者指示相应的药物信息;如此,因目标评估基因只有一个,因此能简单地对CN‑AML患者进行预后评估和疾病分层;另外所述目标评估基因为STAB1基因,因STAB1基因为膜蛋白基因,因此能利用现有的检测设备快速检测出,提高了检测效率。
Description
技术领域
本发明属于分子生物技术领域,尤其涉及一种用于急性髓系白血病药物指示的方法及装置。
背景技术
急性髓性白血病(AML,Acute Myeloid Leukemia)是一种以造血干细胞分化受阻导致不同阶段未成熟细胞异常增殖而正常造血组织减少为特征的恶性克隆性疾病,严重影响人类的健康。
根据细胞遗传学分层可将AML预后分为低危组、中危组和高危组,其中,被划分在中危组的核型正常的(CN-AML,Cytogenetically normal Acute Myeloid Leukemia)患者约占全部AML患者50%,这类患者存在明显异质性,预后存在明显差异。
目前可在CN-AML患者中检测到分子遗传学的变化,包括与CN-AML预后相关的基因突变和基因表达的改变,然后利用分子遗传学并结合联用多项临床指征,比如:基因突变、患者年龄、细胞遗传学水平及基因表达等进行联合打分可以实现CN-AML患者的预后评估。但是这种评估方法一般需要筛选出多个评估基因、导致基因筛选过程、联合打分的过程都非常繁琐,且不能有效地进行预后评估的疾病分层,更不能对相应层级的患者指示对应的药物信息,导致临床不能根据相应的药物指示信息快速制定治疗方案,延误治疗时机。
发明内容
针对现有技术存在的问题,本发明实施例提供了一种用于正常核型急性髓系白血病药物指示的方法及装置,用于解决现有技术中因评估基因太多,导致筛选基因的过程繁琐,且不能有效、简易地利用评估基因对CN-AML患者进行预后评估和疾病分层,不能对相应层级的患者指示对应的药物信息,导致临床不能根据相应的药物指示快速制定治疗方案,延误治疗时机的技术问题。
本发明提供一种用于急性髓系白血病药物指示的方法,所述方法包括:
从数据库中获取目标样本的基因表达信息,所述目标样本为正常核型急性髓系白血病CN-AML患者的样本;
根据预设的生存时间为分类标准,将所述目标样本分为第一类样本及第二类样本,所述第一类样本为生存时间小于2年的样本,所述第二类样本为生存时间大于2年的样本;
根据预设的第一筛选条件对所述第一类样本及所述第二类样本的基因进行筛选,获取多个差异表达基因;
根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;
从所述数据库中获取影响所述CN-AML预后的临床因子,根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的经验预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因;
根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定出所述目标评估基因,所述目标评估基因为STAB1基因;
获取正常核型急性髓系白血病CN-AML患者的目标评估基因;
根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组;
对处于预后不良组的CN-AML患者指示相应的药物信息。
上述方案中,所述从数据库中获取目标样本的基因表达信息,包括:
从癌症基因数据库TCGA中获取预设数量样本的基因表达信息;
根据目标样本的样本标识,从所述预设数量样本中提取目标样本的基因表达信息。
上述方案中,所述根据预设的第一筛选条件对所述第一类样本及所述第二类样本的基因进行筛选,获取多个差异表达基因,包括:
根据所述第一筛选条件对所述第一类样本中的所有基因及所述第二类样本中的所有基因进行筛选,获取多个所述差异表达基因;其中,所述第一筛选条件为假阳性基因的假阳性率FDR<0.05且基因的差异倍数fold-change>1.5。
上述方案中,所述根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因,包括:
获取每个差异表达基因分别在所述目标样本中的生存资料,所述生存资料包括:每个差异表达基因的基因表达量、每个差异表达基因对应样本的生存时间及生存状态;
基于所述每个差异表达基因的生存资料,利用KM统计工具中的曲线函数对每个差异表达基因进行生存分析,生成各差异表达基因的第一KM生存曲线;
从所述各第一KM生存曲线获取每个差异表达基因的第一显著值;
根据预设的第二筛选条件对每个差异表达基因的第一显著值进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;其中,所述第二筛选条件为P≤0.05;所述P为显著值。
上述方案中,所述根据影响所述CN-AML预后的临床因子及影响AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因,包括:
根据影响所述CN-AML预后的临床因子及影响所述AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取每个所述预后相关差异表达基因的第二显著值;
根据预设的第二筛选条件对所述各所述预后相关差异表达基因的第二显著值进行筛选,获取各独立预后基因;其中,所述第二筛选条件为P≤0.05;所述P为显著值;
所述影响所述CN-AML预后的临床因子包括:年龄Age、FMS样的酪氨酸激酶3(FLT3)突变、DNA甲基转移酶3A(DNMT3A)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变、RUNT相关转录因子1(RUNX1)突变及线粒体基因B(MTCYB)突变、核磷蛋白(NPM1)突变及威廉姆斯肿瘤抑制因子1(WT1)突变;
所述影响AML生存的经验预后因子包括:NPM1突变、IDH1突变、IDH2突变及WT1突变。
上述方案中,所述根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定一个目标评估基因,包括:
获取每个独立预后基因分别在所述验证样本集中的各基因表达量,所述验证样本集中包括多个验证样本;
获取每个独立预后基因分别在每个验证样本中的生存资料,所述生存资料包括:每个独立预后基因的基因表达量、每个独立预后基因对应样本的生存时间及生存状态;
基于所述每个独立预后基因的生存资料,利用KM统计工具中的曲线函数对每个独立预后基因进行生存分析,生成各独立预后基因分别在各所述验证样本中的第二KM生存曲线;
从所述各第二KM生存曲线获取每个独立预后基因在各验证样本中的第三显著值;
获取符合所述第二筛选条件的每个独立预后基因各验证样本中的第三显著值的数量,符合第二筛选条件的数量最多的第三显著值对应的独立预后基因为所述目标评估基因。
上述方案中,所述根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者的预后生存时间进行分层,包括:
获取所述目标评估基因在目标样本中的各基因表达量;
确定各基因表达量的中位数;
当所述CN-AML患者目标评估基因的基因表达量大于所述中值时,确定所述CN-AML患者的目标评估基因为高表达基因,将所述当前CN-AML患者分为预后不良组;
当所述CN-AML患者目标评估基因的基因表达量小于所述中值时,确定所述CN-AML患者的目标评估基因为低表达基因,将所述当前CN-AML患者分为预后良好组。
上述方案中,所述对处于预后不良组的CN-AML患者指示相应的药物信息,包括:
获取对目标药物敏感或耐药的基因以及癌症药物数据库中多个药物的拮抗剂的半抑制浓度IC50值,所述目标药物为临床上的经验药物;
利用皮尔森函数确定所述STAB1基因的表达量与所述对药物敏感或者耐药的基因表达量之间的相关性,若所述STAB1基因的表达量与所述对药物敏感或者耐药的基因表达量之间的相关性为正相关,则将所述目标药物作为处于预后不良组的CN-AML患者的指示药物;
利用所述皮尔森函数确定癌症药物数据库中的药物拮抗剂的半抑制浓度IC50值与所述STAB1基因表达量的相关性,若所述IC50值与所述STAB1基因表达量的呈正相关性,则将所述药物作为处于预后不良组的CN-AML患者的指示药物。
上述方案中,所述对处于预后不良组的CN-AML患者指示相应的药物信息,包括:
对处于预后不良组的CN-AML患者指示阿糖胞苷药物、小分子抑制剂NVP-BHG712、小分子抑制剂GSK-J4、小分子抑制剂BRD-K30748066和陶扎塞替Tozasertib药物信息。
本发明还提供一种用于急性髓系白血病药物指示的装置,所述装置包括:
第一获取单元,用于从数据库中获取目标样本的基因表达信息,所述目标样本为正常核型急性髓系白血病CN-AML患者的样本;
分类单元,用于根据预设的生存时间为分类标准,将所述目标样本分为第一类样本及第二类样本,所述第一类样本为生存时间小于2年的样本,所述第二类样本为生存时间大于2年的样本;
第一筛选单元,用于根据预设的第一筛选条件对所述第一类样本及所述第二类样本的基因进行筛选,获取多个差异表达基因;
第二筛选单元,用于根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;
分析单元,用于从所述数据库中获取影响所述CN-AML预后的临床因子,根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的经验预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因;
验证单元,用于根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定出所述目标评估基因,所述目标评估基因为STAB1基因;
第二获取单元,用于获取正常核型急性髓系白血病CN-AML患者的目标评估基因,所述目标评估基因为STAB1基因;
分层单元,用于根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组;
指示单元,用于对处于预后不良组的CN-AML患者指示相应的药物信息。
上述方案中,所述分层单元具体用于:
获取所述目标评估基因在目标样本中的各基因表达量;
确定各基因表达量的中位数;
当所述CN-AML患者目标评估基因的基因表达量大于所述中值时,确定所述CN-AML患者的目标评估基因为高表达基因,将所述当前CN-AML患者分为预后不良组;
当所述CN-AML患者目标评估基因的基因表达量小于所述中值时,确定所述CN-AML患者的目标评估基因为低表达基因,将所述当前CN-AML患者分为预后良好组。
本发明实施例提供了一种用于急性髓系白血病药物指示的方法及装置,所述方法包括:从数据库中获取目标样本的基因表达信息,所述目标样本为正常核型急性髓系白血病CN-AML患者的样本;根据预设的生存时间为分类标准,将所述目标样本分为第一类样本及第二类样本,所述第一类样本为生存时间小于2年的样本,所述第二类样本为生存时间大于2年的样本;根据预设的第一筛选条件对所述第一类样本及所述第二类样本的基因进行筛选,获取多个差异表达基因;根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;从所述数据库中获取影响所述CN-AML预后的临床因子,根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的经验预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因;根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定出所述目标评估基因,所述目标评估基因为STAB1基因;获取正常核型急性髓系白血病CN-AML患者的目标评估基因;根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组;对处于预后不良组的CN-AML患者指示相应的药物信息;如此,因目标评估基因只有一个,因此筛选过程比较简单;然后根据CN-AML患者的目标评估基因把CN-AML患者划分为预后良好组和预后不良组,确定患者的疾病分层;然后根据对预后不良组的CN-AML患者指示相应的药物信息,以辅助临床快速制定预后不良组的治疗方案,避免延误治疗时机;这里,因目标评估基因只有一个,因此能简单地对CN-AML患者进行预后评估和疾病分层;另外所述目标评估基因为STAB1基因,因STAB1基因为膜蛋白基因,因此能利用现有的检测设备快速检测出,提高了检测效率。
附图说明
图1为本发明实施例一提供的用于急性髓系白血病药物指示的方法流程示意图;
图2为本发明实施例二提供的用于急性髓系白血病药物指示的装置结构示意图;
图3为本发明实施例三提供的STAB1基因在TCGACN-AML预测样本中的KM生存曲线示意图;
图4为本发明实施例三提供的STAB1基因在GSE12417A预测样本中的KM生存曲线示意图;
图5为本发明实施例三提供的STAB1基因在GSE71014预测样本中的KM生存曲线示意图;
图6为本发明实施例三提供的STAB1基因在GSE6891预测样本中的KM生存曲线示意图。
具体实施方式
为了解决现有技术中因评估基因太多,导致筛选基因的过程繁琐,且不能有效、简易地利用评估基因对CN-AML患者进行预后的疾病分层及药物指示,导致临床不能根据相应的药物指示尽快制定治疗方案,延误治疗时机的技术问题,本发明实施例提供了一种用于急性髓系白血病药物指示的方法及装置,所述方法包括:获取CN-AML患者的目标评估基因,所述目标评估基因为STAB1基因;根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组;对处于预后不良组的CN-AML患者指示相应的药物信息。
下面通过附图及具体实施例对本发明的技术方案做进一步的详细说明。
实施例一
本实施例提供一种用于急性髓系白血病药物指示的方法,如图1所示,所述方法包括:
S110,从数据库中获取目标样本的基因表达信息,所述目标样本为正常核型急性髓系白血病CN-AML患者的样本;
本步骤中,从癌症基因数据库(TCGA,The Cancer Genome Atlas)下载获取预设数量样本的基因表达信息,预设样本包括:正常核型急性髓系白血病CN-AML患者的和异常核型急性髓系白血病患者的基因表达信息,所述基因表达信息包括:基因表达量。
这里,因在数据中,CN-AML患者的和异常核型急性髓系白血病患者的样本标识是不同的,因此可以根据目标样本的样本标识,从所述预设数量样本中提取目标样本的基因表达信息。所述目标样本为CN-AML患者的基因表达信息。
S111,从数据库中获取目标样本的基因表达信息,所述目标样本为正常核型急性髓系白血病CN-AML患者的样本;
获取到目标样本后,根据预设的生存时间为分类标准,将所述目标样本分为第一类样本及第二类样本;所述预设的生存时间为临床上CN-AML完全缓解的经验指标,具体为2年。本实施中所述第一类样本为生存时间小于2年的样本,所述第二类样本为生存时间大于2年的样本。
S112,根据预设的第一筛选条件对所述第一类样本及所述第二类样本的基因进行筛选,获取多个差异表达基因;
获取到第一类样本及第二类样本后,根据所述第一筛选条件,利用统计建模工具R包的Deseq函数对所述第一类样本中的所有基因及所述第二类样本中的所有基因进行筛选,获取多个所述差异表达基因;其中,所述第一筛选条件为假阳性基因的假阳性率(FDR)<0.05且基因的差异倍数(fold-change)>1.5。
S113,根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;
获取到差异表达基因后,根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因。
具体地,获取每个差异表达基因分别在所述目标样本中的生存资料,所述生存资料包括:每个差异表达基因的基因表达量、每个差异表达基因对应样本的生存时间及生存状态;所述生存状态为生或死的状态,生的状态可以以1对应,死的状态可以以0对应。
基于所述每个差异表达基因的生存资料,利用R包中的曲线函数survival对每个差异表达基因进行生存分析,生成各差异表达基因的第一KM生存曲线;
从所述各第一KM生存曲线获取每个差异表达基因的第一显著值;
根据预设的第二筛选条件对每个差异表达基因的第一显著值进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;其中,所述第二筛选条件为P≤0.05;所述P为显著值。
比如,以基因A为例,基因A在每个目标样本中都存在一个基因表达量,确定目标样本中基因A的表达量取的中位数(中值),利用KM中的Log Rank test函数将表达量大于中值的样本及表达量小于中值的样本进行区分,生成表格;表达量大于中值的样本为高表达量样本,表达量小于中值的样本为低表达量样本。
然后根据基因A的生存资料,利用R包survival函数生成基因A的第一KM生存曲线,并从所述基因A的第一KM生存曲线上读取基因A的显著值P值,当基因A的P值≤0.05时,确定基因A为与所述生存时间相关的预后相关差异表达基因。
S114,从所述数据库中获取影响所述CN-AML预后的临床因子,根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的经验预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因;
确定出所述生存时间相关的多个预后相关差异表达基因后,根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各预后基因。
具体地,首先需要对影响所述CN-AML预后的临床因子进行筛选。
这里还需要从所述数据库中获取影响所述CN-AML预后的临床因子。具体如下:从数据库中获取目标样本的临床信息,利用R包的survival函数,结合目标样本的临床信息,筛选出统计学意义显著的临床因子。当所述临床因子的p≤0.1时,视为该临床因子具有统计学意义,即可以作为影响CN-AML预后的临床因子。本实施例中影响CN-AML预后的临床因子包括:年龄Age、FMS样的酪氨酸激酶3(FLT3)突变、DNA甲基转移酶3A(DNMT3A)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变、RUNT相关转录因子1(RUNX1)突变及线粒体基因B(MTCYB)突变、核磷蛋白(NPM1)突变及威廉姆斯肿瘤抑制因子1(WT1)突变。
然后再结合文献及临床上认定的影响AML生存的经验预后因子,所述影响AML生存的预后因子包括:核磷蛋白(NPM1)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变及威廉姆斯肿瘤抑制因子1(WT1)突变。
最后根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取每个所述预后相关差异表达基因的第二显著值。
当第二显著值P≤0.05时,则说明该基因为独立于Age、FLT3突变、DNMT3A突变、IDH1突变、RUNX1突变、MT_CYB突变、NPM1突变、IDH2突变及WT1突变的独立预后基因。
S115,根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定出所述目标评估基因,所述目标评估基因为STAB1基因;
获取到独立预后基因后,根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定一个目标评估基因。所述验证样本集可以从NCBI网站(https://www.ncbi.nlm.nih.gov/)的GEO数据平台获取。
具体地,获取每个独立预后基因分别在所述验证样本集中的各基因表达量,所述验证样本集中包括多个验证样本;获取每个独立预后基因分别在每个验证样本中的生存资料,所述生存资料包括:每个独立预后基因的基因表达量、每个独立预后基因对应样本的生存时间及生存状态;基于所述每个独立预后基因在各验证样本中的生存资料,利用KM统计工具中的曲线函数对每个预后基因进行生存分析,生成各预后基因分别在各所述验证样本中的第二KM生存曲线;从所述各第二KM生存曲线获取每个预后基因在各验证样本中的第三显著值;获取符合所述第二筛选条件的每个独立预后基因各验证样本中的第三显著值的数量,符合第二筛选条件的数量最多的第三显著值对应的独立预后基因为所述目标评估基因;所述目标评估基因为STAB1基因。
比如,独立预后基因包括基因B、C和D;验证样本包括:a、b和c;以基因B为例来说,获取每个基因B分别在各所述验证样本中的各基因表达量,同时获取在各验证样本中的生存资料。
然后根据基因B在各验证样本中的生存资料,利用R包survival函数分别生成基因B的第三KM生存曲线,并从所述基因B在不同验证样本中的第三KM生存曲线上读取基因B在不同验证样本中的P值,再统计符合第二筛选条件的P值数量,所述第二筛选条件为显著值P≤0.05。
然后以同样的方法统计基因C和D在各验证样本中P值,确定符合第二筛选条件的P值数量。
假设基因B在各验证样本中的P值符合第二筛选条件的数量是3,基因C在各验证样本中的P值符合第二筛选条件的数量是2,基因D在各验证样本中的P值符合第二筛选条件的数量是1,那么就确定基因B为目标评估基因。
这里,因目标评估基因为STAB1基因,可以利用流式检测仪快速获取CN-AML患者的目标评估基因,因STAB1基因为膜蛋白基因,因此能利用现有的检测设备快速检测出。
S116,获取正常核型急性髓系白血病CN-AML患者的目标评估基因,所述目标评估基因为STAB1基因;
信息,利用R包的survival函数,结合目标样本的临床信息,筛选出统计学意义显著的临床因子。当所述临床因子的p≤0.1时,视为该临床因子具有统计学意义,即可以作为影响CN-AML预后的临床因子。本实施例中影响CN-AML预后的临床因子包括:年龄Age、FMS样的酪氨酸激酶3(FLT3)突变、DNA甲基转移酶3A(DNMT3A)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变、RUNT相关转录因子1(RUNX1)突变及线粒体基因B(MTCYB)突变、核磷蛋白(NPM1)突变及威廉姆斯肿瘤抑制因子1(WT1)突变。
然后再结合文献及临床上认定的影响AML生存的经验预后因子,所述影响AML生存的预后因子包括:核磷蛋白(NPM1)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变及威廉姆斯肿瘤抑制因子1(WT1)突变。
最后根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取每个所述预后相关差异表达基因的第二显著值。
当第二显著值P≤0.05时,则说明该基因为独立于Age、FLT3突变、DNMT3A突变、IDH1突变、RUNX1突变、MT_CYB突变、NPM1突变、IDH2突变及WT1突变的独立预后基因。
获取到独立预后基因后,根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定一个目标评估基因。所述验证样本集可以从NCBI网站(https://www.ncbi.nlm.nih.gov/)的GEO数据平台获取。
具体地,获取每个独立预后基因分别在所述验证样本集中的各基因表达量,所述验证样本集中包括多个验证样本;获取每个独立预后基因分别在每个验证样本中的生存资料,所述生存资料包括:每个独立预后基因的基因表达量、每个独立预后基因对应样本的生存时间及生存状态;基于所述每个独立预后基因在各验证样本中的生存资料,利用KM统计工具中的曲线函数对每个预后基因进行生存分析,生成各预后基因分别在各所述验证样本中的第二KM生存曲线;从所述各第二KM生存曲线获取每个预后基因在各验证样本中的第三显著值;获取符合所述第二筛选条件的每个独立预后基因各验证样本中的第三显著值的数量,符合第二筛选条件的数量最多的第三显著值对应的独立预后基因为所述目标评估基因;所述目标评估基因为STAB1基因。
比如,独立预后基因包括基因B、C和D;验证样本包括:a、b和c;以基因B为例来说,获取每个基因B分别在各所述验证样本中的各基因表达量,同时获取在各验证样本中的生存资料。
然后根据基因B在各验证样本中的生存资料,利用R包survival函数分别生成基因B的第三KM生存曲线,并从所述基因B在不同验证样本中的第三KM生存曲线上读取基因B在不同验证样本中的P值,再统计符合第二筛选条件的P值数量,所述第二筛选条件为显著值P≤0.05。
然后以同样的方法统计基因C和D在各验证样本中P值,确定符合第二筛选条件的P值数量。
假设基因B在各验证样本中的P值符合第二筛选条件的数量是3,基因C在各验证样本中的P值符合第二筛选条件的数量是2,基因D在各验证样本中的P值符合第二筛选条件的数量是1,那么就确定基因B为目标评估基因。
这里,因目标评估基因为STAB1基因,可以利用流式检测仪快速获取CN-AML患者的目标评估基因,因STAB1基因为膜蛋白基因,因此能利用现有的检测设备快速检测出。
S117,根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组。
当目标评估基因确定出之后,获取CN-AML患者的目标评估基因,根据CN-AML患者的目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组。
这里,还需要先评估目标评估基因预测生存时间的准确度。具体地,确定目标评估基因在各预测样本中的各基因表达量,确定目标评估基因在各预测样本中的基因表达量的中值,将基因表达量大于该中值的样本确定为评估基因高表达样本,将基因表达量小于该中值的基因确定评估基因低表达样本。
以预设的生存时间为标准,再分别统计出各预测样本中评估基因高表达样本的第一数量,及各预测样本中评估基因低表达样本的第二数量;根据所述第一数量及所述第二数量确定所述目标评估基因在各预测样本中的准确率。其中,所述预测样本包括目标样本及验证样本。
比如,目标评估基因为基因B,在确定基因B在目标样本中的准确率时,先将样本根据基因B分为评估基因高表达样本及评估基因低表达样本,统计出评估基因高表达样本的生存时间小于2年的数量为m,统计出评估基因低表达样本的生存时间大于2年的数量为n,那么在计算目标评估基因B在目标样本的准确率为:(m+n)/S;所述S为目标样本的数量。
当确定出准确率后,当准确率为60%以上时,认为该目标评估基因的准确度是可行的。
那么可认为根据所述目标评估基因的基因表达量评估CN-AML患者的预后生存时间准确度较高。
然后就可以根据CN-AML患者的目标评估基因的基因表达量对所述CN-AML患者进行分层。当所述CN-AML患者目标评估基因的基因表达量大于所述中值时,确定所述CN-AML患者的目标评估基因为高表达基因,将所述当前CN-AML患者分为预后不良组;
当所述CN-AML患者目标评估基因的基因表达量小于所述中值时,确定所述CN-AML患者的目标评估基因为低表达基因,将所述当前CN-AML患者分为预后良好组。
S118,对处于预后不良组的CN-AML患者指示相应的药物信息。
根据目标评估基因表达量将CN-AML患者分为预后良好组和预后不良组别后,可以对处于预后不良组的CN-AML患者指示相应的药物信息,以辅助临床指定精准的分层治疗方案。
具体如下:先需要获取对目标药物敏感或耐药的基因,所述目标药物为临床上的经验药物,比如阿糖胞苷;获得对阿糖胞苷敏感或者耐药的基因后,采用R包的皮尔森函数计算STAB1和对阿糖胞苷敏感或者耐药基因的表达量相关性,发现STAB1高表达与对阿糖胞苷耐药基因表达呈正相关性,即可说明STAB1高表达样本(预后不良组)能揭示这部分患者对阿糖胞苷耐药,即说明在STAB1高表达患者群体中需要加大阿糖胞苷剂量或者选择其他治疗方案。
同时分别计算GDSC(Drug Sensitivity in Cancer)和CTRP(CancerTherapeutics Response Portal)数据库中药物拮抗剂的半抑制浓度IC50值与STAB1表达量的皮尔森相关性,获得p值<0.05(即相关性具有统计学意义)的药物,包括IC50值与STAB1表达量呈现正相关性和负相关性的药物。IC50值与STAB1呈正相关性表明,STAB1高表达样本中,对该类药物具有敏感性,即说明STAB1高表达患者(预后不良组)对该类药物敏感,可供临床考虑为候选药物/小分子抑制剂。这里STAB1高表达患者对小分子抑制剂NVP-BHG712、小分子抑制剂GSK-J4、小分子抑制BRD-K30748066和陶扎塞替(Tozasertib)敏感。
综上,CN-AML患者的预后生存时间小于2年时,为STAB1高表达的患者,该类型患者对阿糖胞苷具有抗药性,在STAB1高表达患者中需要提高阿糖胞苷药物剂量或选择其他替代药物。
同时,STAB1高表达患者中,对小分子抑制剂NVP-BHG712、小分子抑制剂GSK-J4、小分子抑制剂BRD-K30748066和陶扎塞替(Tozasertib)较为敏感,预示着可能这几类药物能在STAB1高表达患者中有一定的作用。这样就对STAB1高表达患者指示出相应的药物信息,以辅助临床指定精准的分层治疗方案,避免延误治疗时机。
实施例二
相应于实施例一,本实施例提供一种用于急性髓系白血病药物指示的装置,如图2所示,所述装置包括:第一获取单元21、分类单元22、第一筛选单元23、第二筛选单元24、分析单元25、验证单元26、第二获取单元27、分层单元28及指示单元29;其中,
第一获取单元21在获取正常核型急性髓系白血病CN-AML患者的目标评估基因之前,需要先确定出目标评估基因是哪一个基因,具体的确定方式如下:
第一获取单元21用于从癌症基因数据库TCGA中下载获取预设数量样本的基因表达信息,预设样本包括:正常核型急性髓系白血病CN-AML患者的和异常核型急性髓系白血病患者的基因表达信息,所述基因表达信息包括:基因表达量。
因在数据中,CN-AML患者的和异常核型急性髓系白血病患者的样本标识是不同的,因此可以根据目标样本的样本标识,从所述预设数量样本中提取目标样本的基因表达信息。所述目标样本即为CN-AML患者的基因表达信息。
当获取到目标样本后,分类单元22用于根据预设的生存时间为分类标准,将所述目标样本分为第一类样本及第二类样本;所述预设的生存时间为临床上CN-AML完全缓解的经验指标,具体为2年。本实施中所述第一类样本为生存时间小于2年的样本,所述第二类样本为生存时间大于2年的样本。
获取到第一类样本及第二类样本后,第一筛选单元23用于根据预设的第一筛选条件,利用R包的Deseq函数对所述第一类样本中的所有基因及所述第二类样本中的所有基因进行筛选,获取多个所述差异表达基因;其中,所述第一筛选条件为差异表达基因的假阳性率(FDR)<0.05且差异表达基因的差异倍数(fold-change)>1.5。
获取到差异表达基因后,第二筛选单元24用于根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因。
具体地,获取每个差异表达基因分别在所述目标样本中的生存资料,所述生存资料包括:每个差异表达基因的基因表达量、每个差异表达基因对应样本的生存时间及生存状态;所述生存状态为生或死的状态,生的状态可以以1对应,死的状态可以以0对应。
基于所述每个差异表达基因的生存资料,利用R包中的曲线函数survival对每个差异表达基因进行生存分析,生成各差异表达基因的第一KM生存曲线;
从所述各第一KM生存曲线获取每个差异表达基因的第一显著值;
根据预设的第二筛选条件对每个差异表达基因的第一显著值进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;其中,所述第二筛选条件为P≤0.05;所述P为显著值。
比如,以基因A为例,基因A在每个目标样本中都存在一个基因表达量,确定目标样本中基因A的表达量取的中值(中位数),利用KM中的Log Rank test函数将表达量大于中值的样本及表达量小于中值的样本进行区分,生成表格;表达量大于中值的样本为高表达量样本,表达量小于中值的样本为低表达量样本。
然后根据基因A的生存资料,利用R包survival函数生成基因A的第一KM生存曲线,并从所述基因A的第一KM生存曲线上读取基因A的显著值P值,当基因A的P值≤0.05时,确定基因A为与所述生存时间相关的差异表达基因。
确定出所述生存时间相关的多个预后相关差异表达基因后,分析单元25用于根据影响所述CN-AML预后的临床因子及影响AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各预后基因。
具体地,分析单元25首先需要对影响所述CN-AML预后的临床因子进行筛选。
这里还需要从数据库中获取目标样本的临床信息,利用R包的survival函数,结合目标样本的临床信息,筛选出统计学意义显著的临床因子。当所述临床因子的p≤0.1时,视为该临床因子具有统计学意义,即可以作为影响CN-AML预后的临床因子。本实施例中影响CN-AML预后的临床因子包括:年龄Age、FMS样的酪氨酸激酶3(FLT3)突变、DNA甲基转移酶3A(DNMT3A)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变、RUNT相关转录因子1(RUNX1)突变及线粒体基因B(MTCYB)突变、核磷蛋白(NPM1)突变及威廉姆斯肿瘤抑制因子1(WT1)突变。
然后再结合文献及临床上认定的影响AML生存的预后因子,所述影响AML生存的预后因子包括:核磷蛋白(NPM1)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变及威廉姆斯肿瘤抑制因子1(WT1)突变。
最后分析单元25根据影响所述CN-AML预后的临床因子及影响AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取每个所述预后相关差异表达基因的第二显著值。
当第二显著值P≤0.05时,则说明该基因为独立于Age、FLT3突变、DNMT3A突变、IDH1突变、RUNX1突变、MTCYB突变、NPM1突变、IDH2突变及WT1mutation突变的独立预后基因。
获取到独立预后基因后,验证单元26用于根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定一个目标评估基因。所述验证样本集可以从NCBI网站(https://www.ncbi.nlm.nih.gov/)的GEO数据平台获取。
具体地,验证单元26获取每个独立预后基因分别在所述验证样本集中的各基因表达量,所述验证样本集中包括多个验证样本;获取每个独立预后基因分别在每个验证样本中的生存资料,所述生存资料包括:每个独立预后基因的基因表达量、每个独立预后基因对应样本的生存时间及生存状态;基于所述每个独立预后基因在各验证样本中的生存资料,利用KM统计工具中的曲线函数对每个预后基因进行生存分析,生成各预后基因分别在各所述验证样本中的第二KM生存曲线;从所述各第二KM生存曲线获取每个预后基因在各验证样本中的第三显著值;获取符合所述第二筛选条件的每个独立预后基因各验证样本中的第三显著值的数量,符合第二筛选条件的数量最多的第三显著值对应的独立预后基因为所述目标评估基因;所述目标评估基因为STAB1基因。
比如,独立预后基因包括基因B、C和D;验证样本包括:a、b和c;以基因B为例来说,验证单元26获取每个基因B分别在各所述验证样本中的各基因表达量,同时获取在各验证样本中的生存资料。
然后根据基因B在各验证样本中的生存资料,利用R包survival函数分别生成基因B的第三KM生存曲线,并从所述基因B在不同验证样本中的第三KM生存曲线上读取基因B在不同验证样本中的P值,再统计符合第二筛选条件的P值数量,所述第二筛选条件为显著值P≤0.05。
然后以同样的方法统计基因C和D在各验证样本中P值,确定符合第二筛选条件的P值数量。
假设基因B在各验证样本中的P值符合第二筛选条件的数量是3,基因C在各验证样本中的P值符合第二筛选条件的数量是2,基因D在各验证样本中的P值符合第二筛选条件的数量是1,那么就确定基因B为目标评估基因。
当确定目标评估基因后,第二获取单元27还用于获取正常核型急性髓系白血病CN-AML患者的目标评估基因,所述目标评估基因为STAB1基因;分层单元27用于根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组。
这里,评估单元30需要先评估目标评估基因预测生存时间的准确度。具体地,确定目标评估基因在各预测样本中的各基因表达量,确定目标评估基因在各预测样本中的基因表达量的中值,将基因表达量大于该中值的样本确定为评估基因高表达样本,将基因表达量小于该中值的基因确定评估基因低表达样本。
评估单元30以预设的生存时间为标准,再分别统计出各预测样本中评估基因高表达样本的第一数量,及各预测样本中评估基因低表达样本的第二数量;根据所述第一数量及所述第二数量确定所述目标评估基因在各预测样本中的准确率。其中,所述预测样本包括目标样本及验证样本。
比如,目标评估基因为基因B,在确定基因B在目标样本中的准确率时,先将样本根据基因B分为评估基因高表达样本及评估基因低表达样本,统计出评估基因高表达样本的生存时间小于2年的数量为m,统计出评估基因低表达样本的生存时间大于2年的数量为n,那么在计算目标评估基因B在目标样本的准确率为:(m+n)/S;所述S为目标样本的数量。
当评估单元30确定出准确率后,当准确率为60%以上时,认为该目标评估基因的准确度是可行的。
那么可认为根据所述目标评估基因的基因表达量评估CN-AML患者的预后生存时间准确度较高。
然后分层单元27就可以根据CN-AML患者的目标评估基因的基因表达量对所述CN-AML患者进行分层。当所述CN-AML患者目标评估基因的基因表达量大于所述中值时,确定所述CN-AML患者的目标评估基因为高表达基因,将所述当前CN-AML患者分为预后不良组;
当所述CN-AML患者目标评估基因的基因表达量小于所述中值时,确定所述CN-AML患者的目标评估基因为低表达基因,将所述当前CN-AML患者分为预后良好组。
分层单元27确定出CN-AML患者的分组后,指示单元28可以对处于预后不良组的CN-AML患者指示相应的药物信息,以辅助临床指定精准的分层治疗方案,避免延误治疗时机。
具体如下:先需要获取对目标药物敏感或耐药的基因,所述目标药物为临床上的经验药物,比如阿糖胞苷;获得对阿糖胞苷敏感或者耐药的基因后,采用R包的皮尔森函数计算STAB1和对阿糖胞苷敏感或者耐药基因的表达量相关性,发现STAB1高表达与对阿糖胞苷耐药基因表达呈正相关性,即可说明STAB1高表达样本(预后不良组)能揭示这部分患者对阿糖胞苷耐药,即说明在STAB1高表达患者群体中需要加大阿糖胞苷剂量或者选择其他治疗方案。
同时分别计算GDSC(Drug Sensitivity in Cancer)和CTRP(CancerTherapeutics Response Portal)数据库中药物拮抗剂的半抑制浓度IC50值与STAB1表达量的皮尔森相关性,获得p值<0.05(即相关性具有统计学意义)的药物,包括IC50值与STAB1表达量呈现正相关性和负相关性的药物。IC50值与STAB1呈正相关性表明,STAB1高表达样本中,对该类药物具有敏感性,即说明STAB1高表达患者(预后不良组)对该类药物敏感,可供临床考虑为候选药物/小分子抑制剂。这里STAB1高表达患者对小分子抑制剂NVP-BHG712、小分子抑制剂GSK-J4、小分子抑制BRD-K30748066和陶扎塞替(Tozasertib)敏感。
综上,CN-AML患者的预后生存时间小于2年时,为STAB1高表达的患者,该类型患者对阿糖胞苷具有抗药性,在STAB1高表达患者中需要提高阿糖胞苷药物剂量或选择其他替代药物。
同时,STAB1高表达患者中,对小分子抑制剂NVP-BHG712、小分子抑制剂GSK-J4、小分子抑制BRD-K30748066和陶扎塞替(Tozasertib)较为敏感,预示着可能这几类药物能在STAB1高表达患者中有一定的作用。这样就对STAB1高表达患者指示出相应的药物信息,以辅助临床指定精准的分层治疗方案,避免延误治疗时机。
实施例三
实际应用中,可以根据上述方法及装置确定CN-AML的目标评估基因,并利用该基因对CN-AML进行预后分层,指示相应的药物信息,具体如下:
首先从TCGA数据库中下载200例样本的基因表达信息及临床信息,然后根据CN-AML样本的样本标识,从所述预设数量样本中提取CN-AML样本的基因表达信息。所述CN-AML样本的数量为79例。
根据预设的生存时间为分类标准,将所述CN-AML样本分为第一类样本及第二类样本;所述预设的生存时间为临床上CN-AML完全缓解的经验指标,具体为2年。本实施中所述第一类样本为生存时间小于2年的CN-AML样本,所述第二类样本为生存时间大于2年的CN-AML样本。
获取到第一类样本及第二类样本后,用R包的Deseq函数对所述第一类样本中的所有基因及所述第二类样本中的所有基因进行筛选,获取多个所述差异表达基因;其中,所述第一筛选条件为差异表达基因的假阳性率(FDR)<0.05且基因的差异倍数(fold-change)>1.5。这里,所述差异表达基因的数量为353。
获取到差异表达基因后,利用R包survival对每个差异表达基因进行生存分析,生成各差异表达基因的第一KM生存曲线;基于第一KM生存曲线获取每个差异表达基因的第一显著值,筛选出显著值P≤0.05的与所述生存时间相关预后相关差异表达基因,所述预后相关差异表达基因的数量为15个。
然后根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因,所述独立预后基因为6个。
本实施例中影响CN-AML预后的临床因子包括:年龄Age、FMS样的酪氨酸激酶3(FLT3)突变、DNA甲基转移酶3A(DNMT3A)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变、RUNT相关转录因子1(RUNX1)突变及线粒体基因B(MTCYB)突变、核磷蛋白(NPM1)突变及威廉姆斯肿瘤抑制因子1(WT1)突变。
所述影响AML生存的预后因子包括:核磷蛋白(NPM1)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变及威廉姆斯肿瘤抑制因子1(WT1)突变。
获取到独立预后基因后,根据预设的基因芯片验证样本集分别对每个独立预后基因进行验证,根据验证结果确定一个目标评估基因。所述验证样本集可以从GEO平台获取。所述验证样本集包括四组,各验证样本的样本数量分别为:79、163、104和187。
获取每个独立预后基因在各验证样本中的显著值P,确定符合第二筛选条件的P值数量;所述第二筛选条件为显著值P≤0.05。符合第二筛选条件的P值数量最多的独立预后基因即为目标评估基因。本实施例中的目标评估基因为STAB1。
然后利用预测样本集来确定STAB1的准确率,具体地,确定STAB1在各预测样本中的各基因表达量,确定STAB1在各预测样本中的基因表达量的中值,将基因表达量大于该中值的基因确定为高表达评估基因,将基因表达量大于该中值的基因确定为低表达基因。
以预设的生存时间为标准,再分别统计出高表达评估基因在各预测样本中的第一数量,及低表达评估基因在各预测样本中的第二数量;根据所述第一数量及所述第二数量确定所述STAB1在各预测样本中的准确率。其中,所述预测样本包括目标样本及验证样本。其中,STAB1在各预测样本中的准确率如表1所示:
表1
在表1中,TCGACN-AML(79)为目标样本,其余四组为预测样本,以目标样本为例,在计算STAB1在目标样本中的准确率时,具体为:(35+22)/79=0.72。
然后再利用同样的方法计算STAB1在其他验证样本中的准确率,由表1可以看出,STAB1在预测样本中的准确率为60%以上,证明STAB1的准确度是可行的。
进一步地,还可根据STAB1在各预测样本中的KM生存曲线来预测STAB1的可行性,STAB1基因在TCGA CN-AML样本KM生存曲线如图3所示,STAB1基因在GSE12417A(79)样本KM生存曲线如图4所示,STAB1基因在GSE71014(104)样本KM生存曲线如图5所示,STAB1基因在GSE6891(187)样本KM生存曲线如图6所示。其中,在图3、图4、图5及图6中,上方的曲线代表的是低表达STAB1对应的KM生存曲线,下方的曲线代表的是高表达STAB1对应的KM生存曲线,对应的n分别代表的是低表达样本的数量及高表达样本的数量。
需要说明的是,图3、图4、图5及图6是根据STAB1的中值来确定低表达样本的数量及高表达样本的数量。而表1中同样以STAB1的中值来确定STAB1高表达组和STAB1低表达组,但是统计的是STAB1低表达样本中生存时间大于2年的样本数量以及STAB1高表达样本中生存时间小于2年的样本数量,因此数量会出现不一致的情况。
由图3、图4、图5及图6可以看出,在四个预测样本中,STAB1基因的P值均小于0.05,也进一步说明STAB1基因的准确度是可行的。
然后就可以根据CN-AML患者的目标评估基因的基因表达量对所述CN-AML患者进行分层。当所述CN-AML患者目标评估基因的基因表达量大于所述中值时,确定所述CN-AML患者的目标评估基因为高表达基因,将所述当前CN-AML患者分为预后不良组;预后不良组CN-AML患者的预后生存时间小于2年;
当所述CN-AML患者目标评估基因的基因表达量小于所述中值时,确定所述CN-AML患者的目标评估基因为低表达基因,将所述当前CN-AML患者分为预后良好组;预后良好组CN-AML患者的预后生存时间大于2年。
根据目标评估基因表达量将CN-AML患者分为预后良好组和预后不良组别后,可以对处于预后不良组的CN-AML患者指示相应的药物信息,以辅助临床指定精准的分层治疗方案,避免延误治疗时机。
具体如下:先需要获取对目标药物敏感或耐药的基因,所述目标药物为临床上的经验药物,比如阿糖胞苷;获得对阿糖胞苷敏感或者耐药的基因后,采用R包的皮尔森函数计算STAB1和对阿糖胞苷敏感或者耐药基因的表达量相关性,发现STAB1高表达与对阿糖胞苷耐药基因表达呈正相关性,即可说明STAB1高表达样本(预后不良组)能揭示这部分患者对阿糖胞苷耐药,即说明在STAB1高表达患者群体中需要加大阿糖胞苷剂量或者选择其他治疗方案。
同时分别计算GDSC(Drug Sensitivity in Cancer)和CTRP(CancerTherapeutics Response Portal)数据库中药物拮抗剂的半抑制浓度IC50值与STAB1表达量的皮尔森相关性,获得p值<0.05(即相关性具有统计学意义)的药物,包括IC50值与STAB1表达量呈现正相关性和负相关性的药物。IC50值与STAB1呈正相关性表明,STAB1高表达样本中,对该类药物具有敏感性,即说明STAB1高表达患者(预后不良组)对该类药物敏感,可供临床考虑为候选药物/小分子抑制剂。这里STAB1高表达患者对小分子抑制剂NVP-BHG712、小分子抑制剂GSK-J4、小分子抑制BRD-K30748066和陶扎塞替(Tozasertib)敏感。
综上,CN-AML患者的预后生存时间小于2年时,为STAB1高表达的患者,该类型患者对阿糖胞苷具有抗药性,在STAB1高表达患者中需要提高阿糖胞苷药物剂量或选择其他替代药物。
同时,STAB1高表达患者中,对小分子抑制剂NVP-BHG712、GSK-J4,BRD-K30748066和陶扎塞替(Tozasertib)较为敏感,预示着可能这几类药物能在STAB1高表达患者中有一定的作用。这样就对STAB1高表达患者指示出相应的药物信息,以辅助临床指定精准的分层治疗方案,避免延误治疗时机。
本发明实施例提供的用于急性髓系白血病药物指示的方法及装置能带来的有益效果至少是:
本发明实施例提供了一种用于急性髓系白血病药物指示的方法及装置,所述方法包括:从数据库中获取目标样本的基因表达信息,所述目标样本为正常核型急性髓系白血病CN-AML患者的样本;根据预设的生存时间为分类标准,将所述目标样本分为第一类样本及第二类样本,所述第一类样本为生存时间小于2年的样本,所述第二类样本为生存时间大于2年的样本;根据预设的第一筛选条件对所述第一类样本及所述第二类样本的基因进行筛选,获取多个差异表达基因;根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;从所述数据库中获取影响所述CN-AML预后的临床因子,根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的经验预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因;根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定出所述目标评估基因,所述目标评估基因为STAB1基因;获取正常核型急性髓系白血病CN-AML患者的目标评估基因;根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组;对处于预后不良组的CN-AML患者指示相应的药物信息;如此,根据筛选出目标样本中的差异表达基因,再对差异表达基因进行筛选,获取与生存时间相关的多个预后相关差异表达基因,再结合临床信息筛选出多个预后基因,然后利用验证样本集中的多组样本数据分别对每个预后基因进行验证,确定出一个目标评估基因,根据CN-AML患者的目标评估基因把CN-AML患者划分为预后良好组和预后不良组,确定患者的疾病分层;然后根据对预后不良组的CN-AML患者指示相应的药物信息,以辅助临床快速制定精准分层的治疗方案,避免延误治疗时机;这里,因确定的目标评估基因只有一个,因此筛选过程简单,并且能简单地对CN-AML患者进行预后评估和疾病分层;另外所述目标评估基因为STAB1基因,因STAB1基因为膜蛋白基因,因此能利用现有的RT-PCR或流式细胞仪易于快速检测出,提高了检测效率,进一步提高了预后分层过程的简便性。
以上所述,仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种用于急性髓系白血病药物指示的方法,其特征在于,所述方法包括:
从数据库中获取目标样本的基因表达信息,所述目标样本为正常核型急性髓系白血病CN-AML患者的样本;
根据预设的生存时间为分类标准,将所述目标样本分为第一类样本及第二类样本,所述第一类样本为生存时间小于2年的样本,所述第二类样本为生存时间大于2年的样本;
根据预设的第一筛选条件对所述第一类样本及所述第二类样本的基因进行筛选,获取多个差异表达基因;
根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;
从所述数据库中获取影响所述CN-AML预后的临床因子,根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的经验预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因;
根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定出所述目标评估基因,所述目标评估基因为STAB1基因;
获取正常核型急性髓系白血病CN-AML患者的目标评估基因;
根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组;
对处于预后不良组的CN-AML患者指示相应的药物信息。
2.如权利要求1所述的方法,其特征在于,所述从数据库中获取目标样本的基因表达信息,包括:
从癌症基因数据库TCGA中获取预设数量样本的基因表达信息;
根据目标样本的样本标识,从所述预设数量样本中提取目标样本的基因表达信息。
3.如权利要求1所述的方法,其特征在于,所述根据预设的第一筛选条件对所述第一类样本及所述第二类样本的基因进行筛选,获取多个差异表达基因,包括:
根据所述第一筛选条件对所述第一类样本中的所有基因及所述第二类样本中的所有基因进行筛选,获取多个所述差异表达基因;其中,所述第一筛选条件为假阳性基因的假阳性率FDR<0.05且基因的差异倍数fold-change>1.5。
4.如权利要求1所述的方法,其特征在于,所述根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因,包括:
获取每个差异表达基因分别在所述目标样本中的生存资料,所述生存资料包括:每个差异表达基因的基因表达量、每个差异表达基因对应样本的生存时间及生存状态;
基于所述每个差异表达基因的生存资料,利用KM统计工具中的曲线函数对每个差异表达基因进行生存分析,生成各差异表达基因的第一KM生存曲线;
从所述各第一KM生存曲线获取每个差异表达基因的第一显著值;
根据预设的第二筛选条件对每个差异表达基因的第一显著值进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;其中,所述第二筛选条件为P≤0.05;所述P为显著值。
5.如权利要求1所述的方法,其特征在于,所述根据影响所述CN-AML预后的临床因子及影响AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因,包括:
根据影响所述CN-AML预后的临床因子及影响所述AML生存的预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取每个所述预后相关差异表达基因的第二显著值;
根据预设的第二筛选条件对所述各所述预后相关差异表达基因的第二显著值进行筛选,获取各独立预后基因;其中,所述第二筛选条件为P≤0.05;所述P为显著值;
所述影响所述CN-AML预后的临床因子包括:年龄Age、FMS样的酪氨酸激酶3(FLT3)突变、DNA甲基转移酶3A(DNMT3A)突变、异磷酸脱氢酶1(IDH1)突变、异磷酸脱氢酶2(IDH2)突变、RUNT相关转录因子1(RUNX1)突变及线粒体基因B(MTCYB)突变、核磷蛋白(NPM1)突变及威廉姆斯肿瘤抑制因子1(WT1)突变;
所述影响AML生存的经验预后因子包括:NPM1突变、IDH1突变、IDH2突变及WT1突变。
6.如权利要求1所述的方法,其特征在于,所述根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定一个目标评估基因,包括:
获取每个独立预后基因分别在所述验证样本集中的各基因表达量,所述验证样本集中包括多个验证样本;
获取每个独立预后基因分别在每个验证样本中的生存资料,所述生存资料包括:每个独立预后基因的基因表达量、每个独立预后基因对应样本的生存时间及生存状态;
基于所述每个独立预后基因的生存资料,利用KM统计工具中的曲线函数对每个独立预后基因进行生存分析,生成各独立预后基因分别在各所述验证样本中的第二KM生存曲线;
从所述各第二KM生存曲线获取每个独立预后基因在各验证样本中的第三显著值;
获取符合所述第二筛选条件的每个独立预后基因各验证样本中的第三显著值的数量,符合第二筛选条件的数量最多的第三显著值对应的独立预后基因为所述目标评估基因。
7.如权利要求1所述的方法,其特征在于,所述根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者的预后生存时间进行分层,包括:
获取所述目标评估基因在目标样本中的各基因表达量;
确定各基因表达量的中位数;
当所述CN-AML患者目标评估基因的基因表达量大于所述中值时,确定所述CN-AML患者的目标评估基因为高表达基因,将所述当前CN-AML患者分为预后不良组;
当所述CN-AML患者目标评估基因的基因表达量小于所述中值时,确定所述CN-AML患者的目标评估基因为低表达基因,将所述当前CN-AML患者分为预后良好组。
8.如权利要求1所述的方法,其特征在于,所述对处于预后不良组的CN-AML患者指示相应的药物信息,包括:
获取对目标药物敏感或耐药的基因以及癌症药物数据库中多个药物的拮抗剂的半抑制浓度IC50值,所述目标药物为临床上的经验药物;
利用皮尔森函数确定所述STAB1基因的表达量与所述对药物敏感或者耐药的基因表达量之间的相关性,若所述STAB1基因的表达量与所述对药物敏感或者耐药的基因表达量之间的相关性为正相关,则将所述目标药物作为处于预后不良组的CN-AML患者的指示药物;
利用所述皮尔森函数确定癌症药物数据库中的药物拮抗剂的半抑制浓度IC50值与所述STAB1基因表达量的相关性,若所述IC50值与所述STAB1基因表达量的呈正相关性,则将所述药物作为处于预后不良组的CN-AML患者的指示药物。
9.如权利要求1所述的方法,其特征在于,所述对处于预后不良组的CN-AML患者指示相应的药物信息,包括:
对处于预后不良组的CN-AML患者指示阿糖胞苷药物、小分子抑制剂NVP-BHG712、小分子抑制剂GSK-J4、小分子抑制剂BRD-K30748066和陶扎塞替Tozasertib药物信息。
10.一种用于急性髓系白血病药物指示的装置,其特征在于,所述装置包括:
第一获取单元,用于从数据库中获取目标样本的基因表达信息,所述目标样本为正常核型急性髓系白血病CN-AML患者的样本;
分类单元,用于根据预设的生存时间为分类标准,将所述目标样本分为第一类样本及第二类样本,所述第一类样本为生存时间小于2年的样本,所述第二类样本为生存时间大于2年的样本;
第一筛选单元,用于根据预设的第一筛选条件对所述第一类样本及所述第二类样本的基因进行筛选,获取多个差异表达基因;
第二筛选单元,用于根据预设的第二筛选条件对所述多个差异表达基因进行筛选,获取与所述生存时间相关的多个预后相关差异表达基因;
分析单元,用于从所述数据库中获取影响所述CN-AML预后的临床因子,根据影响所述CN-AML预后的临床因子及影响急性髓系白血病AML生存的经验预后因子对所述多个预后相关差异表达基因进行多因素比例风险Cox回归分析,获取各独立预后基因;
验证单元,用于根据预设的验证样本集分别对每个独立预后基因进行验证,根据验证结果确定出所述目标评估基因,所述目标评估基因为STAB1基因;
第二获取单元,用于获取正常核型急性髓系白血病CN-AML患者的目标评估基因,所述目标评估基因为STAB1基因;
分层单元,用于根据CN-AML患者目标评估基因的基因表达量对所述CN-AML患者进行分层,将所述CN-AML患者分为预后良好组及预后不良组;
指示单元,用于对处于预后不良组的CN-AML患者指示相应的药物信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810048989.1A CN108130372A (zh) | 2018-01-17 | 2018-01-17 | 一种用于急性髓系白血病药物指示的方法及装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810045606.5A CN107974506A (zh) | 2018-01-17 | 2018-01-17 | 一种用于急性髓系白血病药物指示的方法及装置 |
CN201810048989.1A CN108130372A (zh) | 2018-01-17 | 2018-01-17 | 一种用于急性髓系白血病药物指示的方法及装置 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810045606.5A Division CN107974506A (zh) | 2018-01-17 | 2018-01-17 | 一种用于急性髓系白血病药物指示的方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108130372A true CN108130372A (zh) | 2018-06-08 |
Family
ID=62005994
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810048989.1A Pending CN108130372A (zh) | 2018-01-17 | 2018-01-17 | 一种用于急性髓系白血病药物指示的方法及装置 |
CN201810048986.8A Pending CN108277278A (zh) | 2018-01-17 | 2018-01-17 | 一种用于正常核型急性髓系白血病预后分层的方法及装置 |
CN201810045606.5A Pending CN107974506A (zh) | 2018-01-17 | 2018-01-17 | 一种用于急性髓系白血病药物指示的方法及装置 |
Family Applications After (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810048986.8A Pending CN108277278A (zh) | 2018-01-17 | 2018-01-17 | 一种用于正常核型急性髓系白血病预后分层的方法及装置 |
CN201810045606.5A Pending CN107974506A (zh) | 2018-01-17 | 2018-01-17 | 一种用于急性髓系白血病药物指示的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (3) | CN108130372A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110229902A (zh) * | 2019-06-24 | 2019-09-13 | 至本医疗科技(上海)有限公司 | 用于胃癌预后预测的评估基因群的确定方法 |
CN110942808A (zh) * | 2019-12-10 | 2020-03-31 | 山东大学 | 一种基于基因大数据的预后预测方法及预测*** |
US11568213B2 (en) | 2018-10-29 | 2023-01-31 | Hitachi, Ltd. | Analyzing apparatus, analysis method and analysis program |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108676885A (zh) * | 2018-06-26 | 2018-10-19 | 华中科技大学鄂州工业技术研究院 | 早期肾癌诊断标记物 |
CN110656172A (zh) * | 2019-01-14 | 2020-01-07 | 南方医科大学珠江医院 | 一种预测小细胞肺癌对ep化疗方案敏感性的分子标志物及试剂盒 |
CN112708675A (zh) * | 2020-12-25 | 2021-04-27 | 中山大学肿瘤防治中心 | 骨髓nk细胞联合mcl1抑制剂在抗白血病中的应用 |
CN113345592B (zh) * | 2021-06-18 | 2022-08-23 | 山东第一医科大学附属省立医院(山东省立医院) | 一种急性髓细胞样白血病预后风险模型的构建及诊断设备 |
CN114093524A (zh) * | 2021-11-02 | 2022-02-25 | 深圳市儿童医院 | 一种儿童抗菌药物使用评估***、计算机可读存储介质及终端 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006048263A3 (en) * | 2004-11-04 | 2007-04-12 | Roche Diagnostics Gmbh | Gene expression profiling in acute promyelocytic leukemia |
CN105046094A (zh) * | 2015-08-26 | 2015-11-11 | 深圳谱元科技有限公司 | 肠道菌群的检测***及其方法和动态式数据库 |
-
2018
- 2018-01-17 CN CN201810048989.1A patent/CN108130372A/zh active Pending
- 2018-01-17 CN CN201810048986.8A patent/CN108277278A/zh active Pending
- 2018-01-17 CN CN201810045606.5A patent/CN107974506A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006048263A3 (en) * | 2004-11-04 | 2007-04-12 | Roche Diagnostics Gmbh | Gene expression profiling in acute promyelocytic leukemia |
CN105046094A (zh) * | 2015-08-26 | 2015-11-11 | 深圳谱元科技有限公司 | 肠道菌群的检测***及其方法和动态式数据库 |
Non-Patent Citations (2)
Title |
---|
MING-KAI CHUANG等: ""An mRNA expression signature for prognostication in de novo acute myeloid leukemia patients with normal karyotype"", 《ONCOTARGET》 * |
郭安源等: ""基于新一代高通量技术的个性化医疗研究进展"", 《重庆医学》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11568213B2 (en) | 2018-10-29 | 2023-01-31 | Hitachi, Ltd. | Analyzing apparatus, analysis method and analysis program |
CN110229902A (zh) * | 2019-06-24 | 2019-09-13 | 至本医疗科技(上海)有限公司 | 用于胃癌预后预测的评估基因群的确定方法 |
CN110942808A (zh) * | 2019-12-10 | 2020-03-31 | 山东大学 | 一种基于基因大数据的预后预测方法及预测*** |
Also Published As
Publication number | Publication date |
---|---|
CN107974506A (zh) | 2018-05-01 |
CN108277278A (zh) | 2018-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108130372A (zh) | 一种用于急性髓系白血病药物指示的方法及装置 | |
Read et al. | Evaluating cell-of-origin subtype methods for predicting diffuse large B-cell lymphoma survival: a meta-analysis of gene expression profiling and immunohistochemistry algorithms | |
CN103299188B (zh) | 用于癌症的分子诊断试验 | |
CN109880910A (zh) | 一种肿瘤突变负荷的检测位点组合、检测方法、检测试剂盒及*** | |
Kornblau et al. | Dynamic single-cell network profiles in acute myelogenous leukemia are associated with patient response to standard induction therapy | |
CN107423578A (zh) | 检测体细胞突变的装置 | |
CN111676288B (zh) | 用于预测肺腺癌患者预后的***及其应用 | |
JP2015527057A5 (zh) | ||
CN113234829B (zh) | 结肠癌预后评估基因集及其构建方法 | |
Zheng et al. | Identification and validation of immunotherapy for four novel clusters of colorectal cancer based on the tumor microenvironment | |
CN111440869A (zh) | 一种用于预测原发性乳腺癌发生风险的dna甲基化标记物及其筛选方法和应用 | |
CN110060733A (zh) | 基于单样本的二代测序肿瘤体细胞变异检测装置 | |
Bouhaddou et al. | Caveolin-1 and Sox-2 are predictive biomarkers of cetuximab response in head and neck cancer | |
CN108642568A (zh) | 一种家犬全基因组低密度品种鉴定专用snp芯片设计方法 | |
Sarıman et al. | Investigation of gene expressions of myeloma cells in the bone marrow of multiple myeloma patients by transcriptome analysis | |
Men et al. | A prognostic 11 genes expression model for ovarian cancer | |
Deng et al. | Immune-related genes to construct a novel prognostic model of breast cancer: a chemosensitivity-based study | |
Xie et al. | Accuracy of matrix-assisted LASER desorption ionization–time of flight mass spectrometry for identification of Candida | |
KR102217272B1 (ko) | 유전체 변이 정보를 이용한 질병 진단 바이오마커 추출 방법 | |
CN114093421B (zh) | 一种判别淋巴瘤分子亚型的方法、装置和存储介质 | |
CN113257353B (zh) | 基于reads深度进行目的基因外显子水平缺失检测的方法及装置 | |
CN105969656B (zh) | 一种单细胞外显子测序肿瘤体细胞突变检测及分析平台 | |
CN114540488B (zh) | 一种用于高通量靶向测序检测肿瘤突变负荷的基因组合、检测装置、检测试剂盒及应用 | |
CN108913773A (zh) | 一种临床评价卵巢癌铂类药物化疗敏感性的多分子标志物及其装置与评价方法 | |
CN105177130B (zh) | 用来评估艾滋病人发生免疫重建炎性综合症的标志物 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180608 |