CN111951883A - 一种特征mRNA表达谱组合及结肠癌早期预测方法 - Google Patents
一种特征mRNA表达谱组合及结肠癌早期预测方法 Download PDFInfo
- Publication number
- CN111951883A CN111951883A CN202010775101.1A CN202010775101A CN111951883A CN 111951883 A CN111951883 A CN 111951883A CN 202010775101 A CN202010775101 A CN 202010775101A CN 111951883 A CN111951883 A CN 111951883A
- Authority
- CN
- China
- Prior art keywords
- mrna
- sample
- prediction
- expression
- characteristic
- 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
- 108020004999 messenger RNA Proteins 0.000 title claims abstract description 188
- 230000014509 gene expression Effects 0.000 title claims abstract description 80
- 206010009944 Colon cancer Diseases 0.000 title claims abstract description 59
- 208000029742 colonic neoplasm Diseases 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000000523 sample Substances 0.000 claims abstract description 89
- 238000012706 support-vector machine Methods 0.000 claims abstract description 22
- 239000002773 nucleotide Substances 0.000 claims abstract description 5
- 125000003729 nucleotide group Chemical group 0.000 claims abstract description 5
- 206010028980 Neoplasm Diseases 0.000 claims description 59
- 238000012795 verification Methods 0.000 claims description 31
- 238000012549 training Methods 0.000 claims description 26
- 238000011156 evaluation Methods 0.000 claims description 14
- 238000012360 testing method Methods 0.000 claims description 14
- 238000002790 cross-validation Methods 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 12
- 238000010200 validation analysis Methods 0.000 claims description 12
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 9
- 201000010099 disease Diseases 0.000 claims description 7
- 238000012216 screening Methods 0.000 claims description 5
- 238000012353 t test Methods 0.000 claims description 5
- 201000009030 Carcinoma Diseases 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000000692 Student's t-test Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000010832 independent-sample T-test Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 102100030716 Ankyrin repeat and SOCS box protein 8 Human genes 0.000 claims description 2
- 102100033471 Cbp/p300-interacting transactivator 2 Human genes 0.000 claims description 2
- 102100040325 E3 ubiquitin-protein ligase RNF185 Human genes 0.000 claims description 2
- 102100028953 Gelsolin Human genes 0.000 claims description 2
- 102100034551 Glycolipid transfer protein Human genes 0.000 claims description 2
- 101710094738 Glycolipid transfer protein Proteins 0.000 claims description 2
- 101000703115 Homo sapiens Ankyrin repeat and SOCS box protein 8 Proteins 0.000 claims description 2
- 101000944098 Homo sapiens Cbp/p300-interacting transactivator 2 Proteins 0.000 claims description 2
- 101001104290 Homo sapiens E3 ubiquitin-protein ligase RNF185 Proteins 0.000 claims description 2
- 101001059150 Homo sapiens Gelsolin Proteins 0.000 claims description 2
- 101001056814 Homo sapiens Integral membrane protein 2C Proteins 0.000 claims description 2
- 101000585618 Homo sapiens Leptin receptor gene-related protein Proteins 0.000 claims description 2
- 101000583839 Homo sapiens Muscleblind-like protein 1 Proteins 0.000 claims description 2
- 101000605432 Homo sapiens Phospholipid phosphatase 1 Proteins 0.000 claims description 2
- 101000995300 Homo sapiens Protein NDRG2 Proteins 0.000 claims description 2
- 101000711928 Homo sapiens RING finger protein 11 Proteins 0.000 claims description 2
- 101000707534 Homo sapiens Serine incorporator 1 Proteins 0.000 claims description 2
- 101000706152 Homo sapiens Syntaxin-12 Proteins 0.000 claims description 2
- 101000610729 Homo sapiens Trafficking kinesin-binding protein 2 Proteins 0.000 claims description 2
- 101000894525 Homo sapiens Transforming growth factor-beta-induced protein ig-h3 Proteins 0.000 claims description 2
- 101000766345 Homo sapiens Tribbles homolog 3 Proteins 0.000 claims description 2
- 101000666896 Homo sapiens V-type immunoglobulin domain-containing suppressor of T-cell activation Proteins 0.000 claims description 2
- 102100025464 Integral membrane protein 2C Human genes 0.000 claims description 2
- 102100029825 Leptin receptor gene-related protein Human genes 0.000 claims description 2
- 102100030965 Muscleblind-like protein 1 Human genes 0.000 claims description 2
- 102100038121 Phospholipid phosphatase 1 Human genes 0.000 claims description 2
- 102100034436 Protein NDRG2 Human genes 0.000 claims description 2
- 102100034186 RING finger protein 11 Human genes 0.000 claims description 2
- 102100031707 Serine incorporator 1 Human genes 0.000 claims description 2
- 102100030403 Signal peptide peptidase-like 2A Human genes 0.000 claims description 2
- 101150068300 Sppl2a gene Proteins 0.000 claims description 2
- 102100031117 Syntaxin-12 Human genes 0.000 claims description 2
- 102100040377 Trafficking kinesin-binding protein 2 Human genes 0.000 claims description 2
- 102100021398 Transforming growth factor-beta-induced protein ig-h3 Human genes 0.000 claims description 2
- 102100026390 Tribbles homolog 3 Human genes 0.000 claims description 2
- 102100038282 V-type immunoglobulin domain-containing suppressor of T-cell activation Human genes 0.000 claims description 2
- 108010036805 rap1 GTP-Binding Proteins Proteins 0.000 claims description 2
- 102000012356 rap1 GTP-Binding Proteins Human genes 0.000 claims description 2
- 230000001225 therapeutic effect Effects 0.000 claims description 2
- 108020004414 DNA Proteins 0.000 description 21
- 238000013399 early diagnosis Methods 0.000 description 5
- 208000001333 Colorectal Neoplasms Diseases 0.000 description 4
- 230000034994 death Effects 0.000 description 4
- 231100000517 death Toxicity 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 2
- 208000035475 disorder Diseases 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 208000035984 Colonic Polyps Diseases 0.000 description 1
- 208000008589 Obesity Diseases 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000001684 chronic effect Effects 0.000 description 1
- 206010009887 colitis Diseases 0.000 description 1
- 210000001072 colon Anatomy 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 210000001035 gastrointestinal tract Anatomy 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000036210 malignancy Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 235000020824 obesity Nutrition 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000004393 prognosis Methods 0.000 description 1
- 238000001243 protein synthesis Methods 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 210000000664 rectum Anatomy 0.000 description 1
- 229920002477 rna polymer Polymers 0.000 description 1
- 210000001599 sigmoid colon Anatomy 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000014616 translation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
- G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/20—Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/10—Sequence alignment; Homology search
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
-
- 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
Abstract
本发明公开了一种特征mRNA表达谱组合及结肠癌早期预测方法,所述特征mRNA表达谱组合的核苷酸探针序列如SEQ ID NO.1‑20所示。本发明的基于mRNA表达谱组合特征评估结肠癌早期风险具有很高的精确度和准确率(ROC曲线下面积AUC=1.000)。只需要获取上述20种mRNA的相对表达量,通过支持向量机模型计算给出结肠癌早期患病概率,可作为结肠癌早期预测的参考依据。
Description
技术领域
本发明属于生物技术和医学领域,具体地说,涉及一种特征mRNA表达谱组合及结肠癌早期预测方法。
背景技术
结肠癌(colon cancer)是常见的发生于结肠部位的消化道恶性肿瘤,多发于直肠与乙状结肠交界处。结肠癌男女患病比率为2-3:1,40-50岁人群发病率最高。慢性结肠炎患者、结肠息肉患者、男性肥胖者等为易感人群。结肠癌早期无明显症状,早期诊断较为困难。全球疾病负担(Global Burden of Disease,GBD)数据显示,2017年全球患有结肠直肠癌的人数超过930万,其中中国患病人数高达235万。2017年全球患有结肠直肠癌的死亡人数约为90万,占总死亡人数的1.60%。中国2017年死亡患者数约为19万,占总死亡人数的1.79%。统计结果显示,从1990年到2017年全球结肠直肠癌患病率和死亡率持续增长。中国结肠直肠癌患病率和死亡率在2010年之前低于全球平均水平,2010年后患病率和死亡率急剧增长高于全球平均水平。
支持向量机(Support Vector Machine,SVM)是一类按监督学习方式对数据进行二元分类的广义线性分类器,其决策边界是对学习样本求解的最大边距超平面。SVM模型是将实例表示为空间中的点,这样映射就使得单独类别的实例被尽可能宽的明显的间隔分开。然后,将新的实例映射到同一空间,并基于它们落在间隔的哪一侧来预测所属类别。当训练数据是线性可分时,SVM通过硬间隔最大化学习进行分类。当训练数据线性不可分时,SVM通过使用核技巧以及软间隔最大化学习进行分类。SVM对于特征含义相似的中等大小的数据集很强大,也适用于小型数据集。通常情况下,对样本量小于1万的数据集SVM都有很好的预测效果。SVM在疾病诊断、肿瘤分类、肿瘤基因识别等有着广泛的应用。
肿瘤早期诊断一直是医学界的难题。现有的早期诊断方法多是观测某一个或一类标志物的表达水平,难以达到理想的诊断效果。由于这些标志物在肿瘤患者和正常人群中的表达分布有部分重叠,难以界定标志物的临界值将肿瘤患者和正常人群较好地分开。因此,利用多个标志物表达特征组合可能是肿瘤早期诊断的一种有效方法。Messenger RNA(mRNA)是由DNA的一条链作为模板转录而来的、携带遗传信息的能指导蛋白质合成的一类单链核糖核酸。肿瘤组织与正常组织相比往往表现出大量mRNA的失调,研究表明这些失调的mRNA跟肿瘤的发生、病理机制和预后状态有密切关系。然而,由于单个mRNA分子在肿瘤和正常人群中表达分布有重叠,难以界定早期预测的临界值。
因此,有必要建立一种有助于结肠癌的早期预测的更稳定的多个差异mRNA表达特征组合的诊断模型。
发明内容
有鉴于此,本发明针对上述的问题,提供了一种特征mRNA表达谱组合及结肠癌早期预测方法,能够准确地进行结肠癌I/II期预测。
为了解决上述技术问题,本发明公开了一种特征mRNA表达谱组合,包括ASB8、CITED2、GLTP、GSN、ITM2C、LEPROT、MBNL1、NDRG2、PLPP1、RAP1A、RNF11、RNF185、SERINC1、SPPL2A、STX12、TGFBI、TRAK2、TRIB3、TXNIP和VSIR,其核苷酸探针序列如SEQ ID NO.1-20所示。
本发明还公开了一种基于上述的特征mRNA表达谱组合的结肠癌早期预测方法,包括以下步骤:
步骤1、获取结肠癌早期患者稳定差异表达的特征mRNA;
步骤2、选取特征mRNA表达数据,对每个样本进行数据标准化;
步骤3、使用支持向量机对标准化后的数据构建早期预测模型;
步骤4、根据患者特征mRNA的表达水平进行早期预测;
该方法为非疾病的诊断和治疗目的。
可选地,所述步骤1中的获取结肠癌早期患者稳定差异表达的特征mRNA具体为:
步骤1.1、从Genomic Data Commons Data Portal数据库中下载结肠癌患者肿瘤组织和癌旁组织转录组数据以及临床数据,获得结肠癌患者肿瘤组织基因表达谱readcounts数值,即为测序读段数值,进行对数转换;
步骤1.2、选取具有一定表达丰度的mRNA,即在所有样本中mRNA的read counts大于等于10;再对所有mRNA的read counts取对数,设样本总数为n,筛选后mRNA总数为m,v为mRNA的read counts,u为取对数之后的表达值,则有;
uij=log2vij,i∈(1,n),j∈(1,m) (1)
其中,i为样本编号,j为mRNA编号,uij为第i个样本、第j个mRNA编号取对数之后的表达值,vij为第i个样本、第j个mRNA编号的read counts数值;
步骤1.3、选取疾病分期为I期和II期的结肠癌患者,将这些患者记为结肠癌早期患者,结肠癌早期患者总数记为n′;
步骤1.4、选取肿瘤和正常样本中稳定表达的mRNA,即在肿瘤和正常样本中变异系数均小于0.1的mRNA,设μ为所有样本中mRNA的表达均值,σ为标准差,变异系数的计算公式为:
其中,j为mRNA编号,cv为变异系数,cvj为第j个样本的变异系数,σj为第j个mRNA编号的标准差,μj为第j个mRNA编号的mRNA的表达均值,设m1为稳定表达的mRNA总数,则有:
步骤1.5、选取肿瘤和正常样本中差异表达的mRNA;使用取对数后的表达值计算肿瘤和正常样本mRNA取对数后的倍数变化f,公式为:
其中,j为mRNA编号,fj为第j个mRNA编号的倍数变化,μ1j为第j个mRNA编号的肿瘤样本的表达均值,μ2j为第j个mRNA编号的正常样本的表达均值;
然后使用独立样本t检验比较肿瘤和正常样本中mRNA的表达差异,独立样本t检验公式为:
对所有t检验得出的p值进行错误发现率(false discovery rate,FDR)校正,定义q为FDR校正后的数值,r为p值在m1个mRNA中排序后的位置,则有:
其中,j为mRNA编号,qj代表第j个mRNA编号的FDR校正后的数值,pj代表第j个mRNA编号的t检验得出的p值,rj代表第j个mRNA编号的p值在m1个mRNA中排序后的位置;
最后选取倍数变化f的绝对值大于1且FDR校正后q值小于等于0.05的mRNA,记为特征mRNA,设特征mRNA总数为m2,则有:
m2=m1{|fj|≥1,qj≤0.05},j∈(1,m1) (7)
。
可选地,所述步骤2中的选取特征mRNA表达数据,对每个样本进行数据标准化,公式为:
其中i为样本编号,j为特征mRNA编号;μi为第i个样本所有特征mRNA表达均值,σi为第i个样本所有特征mRNA标准差,uij为取对数后的特征mRNA表达值,uij′为标准化后的mRNA数值。
可选地,所述步骤3中的使用支持向量机对标准化后的数据构建早期预测模型,具体为:
步骤3.1、先对所有样本进行分组:将全部样本中80%划分为训练集+验证集,余下20%划分为测试集;训练集+验证集用于5折交叉验证,即将训练集+验证集分为相等的5组,按顺序将其中一组作为验证集,其余4组作为训练集。给定参数,训练集用于构建模型,验证集用于检验模型精确度;
步骤3.2、最优参数筛选:SVM中参数gamma控制高斯核的宽度,C是正则化参数,限制每个点的重要性;参数网格设置为:
gamma=[0.001,0.01,0.1,1,10,100] (9)
C=[0.001,0.01,0.1,1,10,100] (10)
在交叉验证中,依次使用每两个参数gamma和C的组合构建模型,然后用验证集检验模型精确度;对每个参数组合,5折交叉验证的每次验证产生1个精确度,共进行5次验证即产生5个精确度;选取5次验证的平均精确度最高的参数组合作为最优参数;
步骤3.3、使用最优参数和训练集+验证集的数据构建模型,最后用测试集对模型进行评估。评估指标包括精确度(accuracy)、准确率(precision)、召回率(recall)、特异性(specificity)、F1分数(F1 score)、马修斯相关系数(Matthews correlationcoefficient,MCC)和受试者工作曲线(receiver operating curve,ROC)下面积(areaunder the curve,AUC);在测试集中,定义实际为肿瘤且预测为肿瘤计数为true positive(TP),实际为正常但预测为肿瘤计数为false positive(FP),实际为肿瘤但预测为正常为false negative(FN),实际为正常且预测为正常为true negative(TN);以上评估指标计算公式为:
以上评估指标中精确度、准确率、召回率、特异性、F1分数和AUC返回介于(0,1)之间的值;精确度越高表示模型总体预测效率越高;准确率越高说明犯I类错误越小;召回率越高说明犯II类错误越小;特异性高说明在预测为正例的样本中很少有负例混入;F1分数是一个综合指标,为准确率和召回率的调和平均;MCC是观察到的和预测的二元分类之间的相关系数,返回介于(-1,1)之间的值,其中1表示完美预测,0表示不比随机预测好,-1表示预测和观察之间的完全不一致;AUC越高表明分类器预测的正实例概率越高;因此,以上指标越接近1表明模型整体的预测效果越好;
步骤3.4、若以上评估指标都大于0.9,说明模型具有较好的预测效果,则使用所有数据,用最优参数组合构建最终预测模型。
可选地,所述步骤4中的根据患者特征mRNA的表达水平进行早期预测,具体为:
步骤4.1、对预测样本的特征mRNA表达数据进行标准化,设u为预测样本特征mRNA表达值,μ为预测样本特征mRNA表达均值,σ为预测样本特征mRNA标准差,公式为:
其中j为特征mRNA编号,uj′为标准化后的mRNA数值;
步骤4.2、将预测样本标准化后的mRNA数值代入最终预测进行预测;预测结果为1表示患有结肠癌,预测结果为0表示正常。
与现有技术相比,本发明可以获得包括以下技术效果:
1)预测速度快:使用本发明构建的预测模型可以对大规模样本进行快速预测,100个样本的预测时间只需要几秒钟。
2)准确度高:本发明构建的预测模型预测精确度和准确率较高,都达到90%以上,ROC曲线下面积AUC可达1.000。
3)平台异质性影响较小:由于不同分析平台测定的mRNA表达值有较大差异,本发明预测使用标准化后的特征mRNA表达值,因此受平台异质性的影响较小。
当然,实施本发明的任一产品并不一定需要同时达到以上所述的所有技术效果。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明数据筛选和模型构建的流程;
图2是本发明支持向量机模型交叉验证参数优化过程;
图3是本发明支持向量机模型测试集评估指标;
图4是本发明支持向量机模型测试集ROC曲线。
具体实施方式
以下将配合实施例来详细说明本发明的实施方式,藉此对本发明如何应用技术手段来解决技术问题并达成技术功效的实现过程能充分理解并据以实施。
一种基于特征mRNA表达谱组合的结肠癌早期预测方法,包括以下步骤:
步骤1、获取结肠癌早期患者稳定差异表达的mRNA(特征mRNA),具体为:
步骤1.1、从Genomic Data Commons Data Portal数据库中下载结肠癌患者肿瘤组织和癌旁组织转录组数据以及临床数据,获得结肠癌患者肿瘤组织基因表达谱readcounts数值,即为测序读段数值,进行对数转换;
步骤1.2、选取具有一定表达丰度的mRNA,即在所有样本中mRNA的read counts大于等于10。再对所有mRNA的read counts取对数,设样本总数为n,筛选后mRNA总数为m,v为mRNA的read counts,u为取对数之后的表达值,则有;
uij=log2vij,i∈(1,n),j∈(1,m) (1)
其中,i为样本编号,j为mRNA编号,uij为第i个样本、第j个mRNA编号取对数之后的表达值,vij为第i个样本、第j个mRNA编号的read counts数值。
步骤1.3、选取疾病分期为I期和II期的结肠癌患者,将这些患者记为结肠癌早期患者,结肠癌早期患者总数记为n′;
步骤1.4、选取肿瘤和正常样本中稳定表达的mRNA,即在肿瘤和正常样本中变异系数均小于0.1的mRNA,设μ为所有样本中mRNA的表达均值,σ为标准差,变异系数的计算公式为:
其中,j为mRNA编号,cv为变异系数,cvj为第j个样本的变异系数,σj为第j个mRNA编号的标准差,μj为第j个mRNA编号的mRNA的表达均值,设m1为稳定表达的mRNA总数,则有:
步骤1.5、选取肿瘤和正常样本中差异表达的mRNA。使用取对数后的表达值计算肿瘤和正常样本mRNA取对数后的倍数变化f,公式为:
其中j为mRNA编号,fj为第j个mRNA编号的倍数变化,μ1j为第j个mRNA编号的肿瘤样本的表达均值,μ2j为第j个mRNA编号的正常样本的表达均值。
然后使用独立样本t检验比较肿瘤和正常样本中mRNA的表达差异,独立样本t检验公式为:
对所有t检验得出的p值进行错误发现率(false discovery rate,FDR)校正,定义q为FDR校正后的数值,r为p值在m1个mRNA中排序后的位置,则有:
其中,j为mRNA编号,qj代表第j个mRNA编号的FDR校正后的数值,pj代表第j个mRNA编号的t检验得出的p值,rj代表第j个mRNA编号的p值在m1个mRNA中排序后的位置。
最后选取倍数变化f的绝对值大于1且FDR校正后q值小于等于0.05的mRNA,记为特征mRNA,设特征mRNA总数为m2,则有:
m2=m1{|fj|≥1,qj≤0.05},j∈(1,m1) (7)
步骤2、选取特征mRNA表达数据,对每个样本进行数据标准化,公式为:
其中i为样本编号,j为特征mRNA编号。μi为第i个样本所有特征mRNA表达均值,σi为第i个样本所有特征mRNA标准差,uij为取对数后的特征mRNA表达值,uij′为标准化后的mRNA数值。
步骤3、使用支持向量机对标准化后的数据构建早期预测模型,具体为:
步骤3.1、先对所有样本进行分组。将全部样本中80%划分为训练集+验证集,余下20%划分为测试集。训练集+验证集用于5折交叉验证,即将训练集+验证集分为相等的5组,按顺序将其中一组作为验证集,其余4组作为训练集。给定参数,训练集用于构建模型,验证集用于检验模型精确度。
步骤3.2、最优参数筛选。SVM中参数gamma控制高斯核的宽度,C是正则化参数,限制每个点的重要性。参数网格设置为:
gamma=[0.001,0.01,0.1,1,10,100] (9)
C=[0.001,0.01,0.1,1,10,100] (10)
在交叉验证中,依次使用每两个参数gamma和C的组合构建模型,然后用验证集检验模型精确度。对每个参数组合,5折交叉验证的每次验证产生1个精确度,共进行5次验证即产生5个精确度。选取5次验证的平均精确度最高的参数组合作为最优参数。
步骤3.3、使用最优参数和训练集+验证集的数据构建模型,最后用测试集对模型进行评估。评估指标包括精确度(accuracy)、准确率(precision)、召回率(recall)、特异性(specificity)、F1分数(F1 score)、马修斯相关系数(Matthews correlationcoefficient,MCC)和受试者工作曲线(receiver operating curve,ROC)下面积(areaunder the curve,AUC)。在测试集中,定义实际为肿瘤且预测为肿瘤计数为true positive(TP),实际为正常但预测为肿瘤计数为false positive(FP),实际为肿瘤但预测为正常为false negative(FN),实际为正常且预测为正常为true negative(TN)。以上评估指标计算公式为:
以上评估指标中精确度、准确率、召回率、特异性、F1分数和AUC返回介于(0,1)之间的值。精确度越高表示模型总体预测效率越高;准确率越高说明犯I类错误越小;召回率越高说明犯II类错误越小;特异性高说明在预测为正例的样本中很少有负例混入;F1分数是一个综合指标,为准确率和召回率的调和平均;MCC是观察到的和预测的二元分类之间的相关系数,返回介于(-1,1)之间的值,其中1表示完美预测,0表示不比随机预测好,-1表示预测和观察之间的完全不一致;AUC越高表明分类器预测的正实例概率越高。因此,以上指标越接近1表明模型整体的预测效果越好。
步骤3.4、若以上评估指标都大于0.9,说明模型具有较好的预测效果。则使用所有数据,用最优参数组合构建最终预测模型。
步骤4、根据患者特征mRNA的表达水平进行早期预测,具体为:
步骤4.1、对预测样本的特征mRNA表达数据进行标准化,设u为预测样本特征mRNA表达值,μ为预测样本特征mRNA表达均值,σ为预测样本特征mRNA标准差,公式为:
其中j为特征mRNA编号,uj′为标准化后的mRNA数值。
步骤4.2、将预测样本标准化后的mRNA数值代入最终预测进行预测。预测结果为1表示患有结肠癌,预测结果为0表示正常。
实施例1
步骤1、获取结肠癌早期患者稳定差异表达的mRNA(特征mRNA),详细流程见图1。
步骤1.1、从Genomic Data Commons Data Portal数据库中下载结肠癌患者肿瘤组织和癌旁组织转录组数据以及临床数据,获得结肠癌患者肿瘤组织基因表达谱readcounts数值,进行对数转换。
步骤1.2、选取具有一定表达丰度的mRNA,即在所有样本中mRNA的read counts大于等于10,详见公式(1)。
步骤1.3、选取疾病分期为I期和II期的结肠癌患者,详见公式(2)-(3),将这些患者记为结肠癌早期患者。
步骤1.4、选取肿瘤和正常样本中稳定表达的mRNA,即在肿瘤和正常样本中变异系数均小于0.1的mRNA。
步骤1.5、选取肿瘤和正常样本中差异表达的mRNA,详见公式(4)-(7)。记为特征mRNA。本例中选取前20个结肠癌特征mRNA(按FDR校正后P值从小到大排序)进行模型构建,见表1。20个结肠癌特征mRNA的核苷酸探针序列见表2。
表1.结肠癌特征mRNA
表2.结肠癌特征mRNA的核苷酸探针序列
步骤2、对每个样本进行数据标准化,详见公式(8)。
步骤3、使用支持向量机对标准化后的数据构建早期诊断模型。
步骤3.1、先对所有样本进行分组。将全部样本中80%划分为训练集+验证集,余下20%划分为测试集。训练集+验证集用于5折交叉验证,即将训练集+验证集分为相等的5组,按顺序将其中一组作为验证集,其余4组作为训练集。给定参数,训练集用于构建模型,验证集用于检验模型精确度。详见图1。
步骤3.2、最优参数筛选。SVM参数网格设置见公式(9)-(10)。在交叉验证中,依次使用每两个参数gamma和C的组合构建模型,然后用验证集检验模型精确度。对每个参数组合,5折交叉验证的每次验证产生1个精确度,共进行5次验证即产生5个精确度。选取5次验证的平均精确度最高的参数组合作为最优参数。图2所示为交叉验证参数优化过程,当参数gamma=0.01,参数C=100时模型交叉验证精确度最高:1.000。因此该模型的最优参数为:gamma=0.01,C=100。
步骤3.3、使用最优参数和训练集+验证集的数据构建模型,最后用测试集对模型进行评估。评估指标包括精确度(accuracy)、准确率(precision)、召回率(recall)、特异性(specificity)、F1分数(F1 score)、马修斯相关系数(Matthews correlationcoefficient,MCC)和受试者工作曲线(receiver operating curve,ROC)下面积(areaunder the curve,AUC)。评估指标详见公式(11)-(17)。
步骤3.4、图3所示为以上评估指标中的精确度、准确率、召回率、特异性、F1分数和MCC,这6个指标均为1.000;图4所示为ROC曲线和AUC,测试集中AUC为1.000。说明以上评估指标说明该模型有很好的预测效果。因此使用所有数据,用最优参数组合构建最终预测模型。
步骤4、根据患者特征mRNA的表达水平进行早期预测:
步骤4.1、对预测样本的特征mRNA表达数据进行标准化,详见公式(18)。本发明随机选取10例样本进行预测,并在构建最终预测模型时将这10例样本剔除。所选取的10例样本编号和标准化后特征mRNA数值见表3。
表3. 10例样本编号和特征mRNA标准化后的数值
步骤4.2、将预测样本标准化后的mRNA数值代入最终预测进行预测。预测结果为1表示患有结肠癌,预测结果为0表示正常。10例样本编号,对应的TCGA编号,实际状态和预测结果见表4。10例样本预测结果与实际状态完全符合,说明本发明可以对结肠癌进行精确的早期预测。
表4. 10例样本编号,对应的TCGA编号,实际和预测的状态
综上所述,本发明的特征mRNA表达谱组合具有很高的预测准确性,能够有效地进行结肠癌的早期预测。此外,本发明没有平台依赖性,能够对多种来源的数据进行预测。
上述说明示出并描述了发明的若干优选实施例,但如前所述,应当理解发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述发明构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离发明的精神和范围,则都应在发明所附权利要求的保护范围内。
SEQUENCE LISTING
<110> 广东省第二人民医院
<120> 一种特征mRNA表达谱组合及结肠癌早期预测方法
<130> 2020
<160> 20
<170> PatentIn version 3.3
<210> 1
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 1
gctctgtaat aacagtaata aatagctctg 30
<210> 2
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 2
ccacccccat atgtacagat gataataggg 30
<210> 3
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 3
atcgatgtca tctacgagac gtatacccag 30
<210> 4
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 4
gtgaagcaag gctttgagcc tccctccttt 30
<210> 5
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 5
cttagcttgt actttggacg cgtttctata 30
<210> 6
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 6
aaccaggatt atcatctata ttggaagtca 30
<210> 7
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 7
atcctttcaa accctcatga ctgacaaaaa 30
<210> 8
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 8
ggagggtgcc ctcccgtctc ccacaacttc 30
<210> 9
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 9
agcatgcttc ttgctaacat ttgggctcag 30
<210> 10
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 10
gaaatgtccc ctcaaactca ttgcagcaga 30
<210> 11
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 11
tgtagtatcc atatgttgct taaatttcct 30
<210> 12
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 12
ttatgactgc atagtttgtg gaaacaaaga 30
<210> 13
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 13
caggtcagaa gaatgatgga atgttttaga 30
<210> 14
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 14
atgttcagac tggttcttct tacatatact 30
<210> 15
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 15
gcacagtggc ccaaagagca gcttcagaga 30
<210> 16
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 16
agctatgagt tgaaatgttc tgtcaaatgt 30
<210> 17
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 17
agaaaatgtt gtgctgtatg ttcttgattt 30
<210> 18
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 18
acacttgggg tccacaatcc caggtccata 30
<210> 19
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 19
tccacccttt tctgagagtt attacagcca 30
<210> 20
<211> 30
<212> DNA
<213> 人工序列(Artificial sequence)
<400> 20
ctgatcttgg cgtgtagtcc tgcacctgtt 30
Claims (6)
1.一种基于特征mRNA表达谱组合,其特征在于,包括ASB8、CITED2、GLTP、GSN、ITM2C、LEPROT、MBNL1、NDRG2、PLPP1、RAP1A、RNF11、RNF185、SERINC1、SPPL2A、STX12、TGFBI、TRAK2、TRIB3、TXNIP和VSIR,其核苷酸探针序列如SEQ ID NO.1-20所示。
2.一种基于权利要求1所述的特征mRNA表达谱组合的结肠癌早期预测方法,其特征在于,包括以下步骤:
步骤1、获取结肠癌早期患者稳定差异表达的特征mRNA;
步骤2、选取特征mRNA表达数据,对每个样本进行数据标准化;
步骤3、使用支持向量机对标准化后的数据构建早期预测模型;
步骤4、根据患者特征mRNA的表达水平进行早期预测;
该方法为非疾病的诊断和治疗目的。
3.根据权利要求2所述的预测方法,其特征在于,所述步骤1中的获取结肠癌早期患者稳定差异表达的特征mRNA具体为:
步骤1.1、从Genomic Data Commons Data Portal数据库中下载结肠癌患者肿瘤组织和癌旁组织转录组数据以及临床数据,获得结肠癌患者肿瘤组织基因表达谱read counts数值,即为测序读段数值,进行对数转换;
步骤1.2、选取具有一定表达丰度的mRNA,即在所有样本中mRNA的read counts大于等于10;再对所有mRNA的read counts取对数,设样本总数为n,筛选后mRNA总数为m,v为mRNA的read counts,u为取对数之后的表达值,则有:
uij=log2vij,i∈(1,n),j∈(1,m) (1)
其中,i为样本编号,j为mRNA编号,uij为第i个样本、第j个mRNA编号取对数之后的表达值,vij为第i个样本、第j个mRNA编号的read counts数值;
步骤1.3、选取疾病分期为I期和II期的结肠癌患者,将这些患者记为结肠癌早期患者,结肠癌早期患者总数记为n′;
步骤1.4、选取肿瘤和正常样本中稳定表达的mRNA,即在肿瘤和正常样本中变异系数均小于0.1的mRNA,设μ为所有样本中mRNA的表达均值,σ为标准差,变异系数的计算公式为:
其中,j为mRNA编号,cv为变异系数,cvj为第j个样本的变异系数,σj为第j个mRNA编号的标准差,μj为第j个mRNA编号的mRNA的表达均值,设m1为稳定表达的mRNA总数,则有:
步骤1.5、选取肿瘤和正常样本中差异表达的mRNA;使用取对数后的表达值计算肿瘤和正常样本mRNA取对数后的倍数变化f,公式为:
其中,j为mRNA编号,fj为第j个mRNA编号的倍数变化,μ1j为第j个mRNA编号的肿瘤样本的表达均值,μ2j为第j个mRNA编号的正常样本的表达均值;
然后使用独立样本t检验比较肿瘤和正常样本中mRNA的表达差异,独立样本t检验公式为:
对所有t检验得出的p值进行错误发现率(false discovery rate,FDR)校正,定义q为FDR校正后的数值,r为p值在m1个mRNA中排序后的位置,则有:
其中,j为mRNA编号,qj代表第j个mRNA编号的FDR校正后的数值,pj代表第j个mRNA编号的t检验得出的p值,rj代表第j个mRNA编号的p值在m1个mRNA中排序后的位置;
最后选取倍数变化f的绝对值大于1且FDR校正后q值小于等于0.05的mRNA,记为特征mRNA,设特征mRNA总数为m2,则有:
5.根据权利要求2所述的预测方法,其特征在于,所述步骤3中的使用支持向量机对标准化后的数据构建早期预测模型,具体为:
步骤3.1、先对所有样本进行分组:将全部样本中80%划分为训练集+验证集,余下20%划分为测试集;训练集+验证集用于5折交叉验证,即将训练集+验证集分为相等的5组,按顺序将其中一组作为验证集,其余4组作为训练集;给定参数,训练集用于构建模型,验证集用于检验模型精确度;
步骤3.2、最优参数筛选:SVM中参数gamma控制高斯核的宽度,C是正则化参数,限制每个点的重要性;参数网格设置为:
gamma=[0.001,0.01,0.1,1,10,100] (9)
C=[0.001,0.01,0.1,1,10,100] (10)
在交叉验证中,依次使用每两个参数gamma和C的组合构建模型,然后用验证集检验模型精确度;对每个参数组合,5折交叉验证的每次验证产生1个精确度,共进行5次验证即产生5个精确度;选取5次验证的平均精确度最高的参数组合作为最优参数;
步骤3.3、使用最优参数和训练集+验证集的数据构建模型,最后用测试集对模型进行评估,评估指标包括精确度(accuracy)、准确率(precision)、召回率(recall)、特异性(specificity)、F1分数(F1 score)、马修斯相关系数(Matthews correlationcoefficient,MCC)和受试者工作曲线(receiver operating curve,ROC)下面积(areaunder the curve,AUC);在测试集中,定义实际为肿瘤且预测为肿瘤计数为true positive(TP),实际为正常但预测为肿瘤计数为false positive(FP),实际为肿瘤但预测为正常为false negative(FN),实际为正常且预测为正常为true negative(TN);以上评估指标计算公式为:
以上评估指标中精确度、准确率、召回率、特异性、F1分数和AUC返回介于(0,1)之间的值;精确度越高表示模型总体预测效率越高;准确率越高说明犯I类错误越小;召回率越高说明犯II类错误越小;特异性高说明在预测为正例的样本中很少有负例混入;F1分数是一个综合指标,为准确率和召回率的调和平均;MCC是观察到的和预测的二元分类之间的相关系数,返回介于(-1,1)之间的值,其中1表示完美预测,0表示不比随机预测好,-1表示预测和观察之间的完全不一致;AUC越高表明分类器预测的正实例概率越高;因此,以上指标越接近1表明模型整体的预测效果越好;
步骤3.4、若以上评估指标都大于0.9,说明模型具有较好的预测效果,则使用所有数据,用最优参数组合构建最终预测模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010775101.1A CN111951883A (zh) | 2020-08-04 | 2020-08-04 | 一种特征mRNA表达谱组合及结肠癌早期预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010775101.1A CN111951883A (zh) | 2020-08-04 | 2020-08-04 | 一种特征mRNA表达谱组合及结肠癌早期预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111951883A true CN111951883A (zh) | 2020-11-17 |
Family
ID=73339497
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010775101.1A Pending CN111951883A (zh) | 2020-08-04 | 2020-08-04 | 一种特征mRNA表达谱组合及结肠癌早期预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111951883A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1512758A1 (en) * | 2003-08-27 | 2005-03-09 | Veridex, LLC | Colorectal cancer prognostics |
US20070237770A1 (en) * | 2001-11-30 | 2007-10-11 | Albert Lai | Novel compositions and methods in cancer |
WO2009047515A1 (en) * | 2007-10-10 | 2009-04-16 | Paul Jenkins | Method of identifying cancer risk markers |
US20140094379A1 (en) * | 2011-01-25 | 2014-04-03 | Almac Diagnostics Limited | Colon Cancer Gene Expression Signatures and Methods of Use |
CN108034708A (zh) * | 2017-12-14 | 2018-05-15 | 中国中医科学院中药研究所 | 通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统 |
CN110714078A (zh) * | 2019-09-29 | 2020-01-21 | 浙江大学 | 一种用于ii期结直肠癌复发预测的标记基因及应用 |
-
2020
- 2020-08-04 CN CN202010775101.1A patent/CN111951883A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070237770A1 (en) * | 2001-11-30 | 2007-10-11 | Albert Lai | Novel compositions and methods in cancer |
EP1512758A1 (en) * | 2003-08-27 | 2005-03-09 | Veridex, LLC | Colorectal cancer prognostics |
WO2009047515A1 (en) * | 2007-10-10 | 2009-04-16 | Paul Jenkins | Method of identifying cancer risk markers |
US20140094379A1 (en) * | 2011-01-25 | 2014-04-03 | Almac Diagnostics Limited | Colon Cancer Gene Expression Signatures and Methods of Use |
CN108034708A (zh) * | 2017-12-14 | 2018-05-15 | 中国中医科学院中药研究所 | 通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统 |
CN110714078A (zh) * | 2019-09-29 | 2020-01-21 | 浙江大学 | 一种用于ii期结直肠癌复发预测的标记基因及应用 |
Non-Patent Citations (1)
Title |
---|
梅端;柴华金;黄江;: "一种改进的SVM算法对miRNA表达谱的分析", 西北师范大学学报(自然科学版), no. 02, pages 51 - 56 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112020565A (zh) | 用于确保基于测序的测定的有效性的质量控制模板 | |
CN111748632A (zh) | 一种特征lincRNA表达谱组合及肝癌早期预测方法 | |
CA2877430C (en) | Systems and methods for generating biomarker signatures with integrated dual ensemble and generalized simulated annealing techniques | |
CN104508670B (zh) | 用于生成生物标志物签名的系统和方法 | |
CN112927757B (zh) | 基于基因表达和dna甲基化数据的胃癌生物标志物识别方法 | |
CN111748633A (zh) | 一种特征miRNA表达谱组合及头颈鳞状细胞癌早期预测方法 | |
WO2023197825A1 (zh) | 多癌种早筛模型构建方法以及检测装置 | |
CN115631789B (zh) | 一种基于泛基因组的群体联合变异检测方法 | |
CN111763738A (zh) | 一种特征mRNA表达谱组合及肝癌早期预测方法 | |
CN111748634A (zh) | 一种特征lincRNA表达谱组合及结肠癌的早期预测方法 | |
CN111944900A (zh) | 一种特征lincRNA表达谱组合及子宫内膜癌早期预测方法 | |
CN111733251A (zh) | 一种特征miRNA表达谱组合及肾透明细胞癌早期预测方法 | |
CN111944902A (zh) | 一种基于lincRNA表达谱组合特征的肾乳头状细胞癌早期预测方法 | |
CN116312800A (zh) | 一种基于血浆中循环rna全转录组测序的肺癌特征识别方法、装置和存储介质 | |
Vishwakarma et al. | Classification algorithm for high‐dimensional protein markers in time‐course data | |
CN111808965A (zh) | 一种特征lincRNA表达谱组合及肾透明细胞癌早期预测方法 | |
CN111733252A (zh) | 一种特征miRNA表达谱组合及胃癌早期预测方法 | |
CN111850124A (zh) | 一种特征lincRNA表达谱组合及肺鳞癌早期预测方法 | |
CN111951883A (zh) | 一种特征mRNA表达谱组合及结肠癌早期预测方法 | |
CN111793692A (zh) | 一种特征miRNA表达谱组合及肺鳞癌早期预测方法 | |
CN115537467A (zh) | 基于深度神经网络的卵巢癌生存预后预测分子模型的建立方法及其应用 | |
CN111944901A (zh) | 一种特征mRNA表达谱组合及肾乳头状细胞癌早期预测方法 | |
CN111718997A (zh) | 一种特征mRNA表达谱组合及胃癌早期预测方法 | |
CN111944898A (zh) | 一种特征mRNA表达谱组合及肾透明细胞癌早期预测方法 | |
CN111793691A (zh) | 一种特征mRNA表达谱组合及肺鳞癌早期预测方法 |
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 |