一种特征mRNA表达谱组合及头颈鳞状细胞癌早期预测方法技术

本发明专利技术公开了一种特征mRNA表达谱组合及头颈鳞状细胞癌早期预测方法，所述mRNA核苷酸探针序列如SEQ ID NO.1‑20所示。本发明专利技术的基于mRNA表达谱组合特征评估头颈鳞状细胞癌早期风险具有很高的精确度和准确率(ROC曲线下面积AUC＝1.000)。只需要获取上述20种mRNA的相对表达量，通过支持向量机模型计算给出头颈鳞状细胞癌早期患病概率，可作为头颈鳞状细胞癌早期预测的参考依据。

【技术实现步骤摘要】
一种特征mRNA表达谱组合及头颈鳞状细胞癌早期预测方法
本专利技术属于生物技术和医学领域，具体地说，涉及一种特征mRNA表达谱组合及头颈鳞状细胞癌早期预测方法。
技术介绍
头颈鳞状细胞癌(headandnecksquamouscellcarcinoma)占头颈癌的90％，是起源于上呼吸道细胞的一类解破学和分布多样的恶性肿瘤，包括唇和口腔、口咽、下咽、喉、鼻旁窦和唾液腺的恶性肿瘤。头颈鳞状细胞癌通常起始于衬于粘膜表面的鳞状细胞，最常见的头颈鳞状细胞癌类型是位于口腔和口咽部的肿瘤。支持向量机(SupportVectorMachine,SVM)是一类按监督学习方式对数据进行二元分类的广义线性分类器，其决策边界是对学习样本求解的最大边距超平面。SVM模型是将实例表示为空间中的点，这样映射就使得单独类别的实例被尽可能宽的明显的间隔分开。然后，将新的实例映射到同一空间，并基于它们落在间隔的哪一侧来预测所属类别。当训练数据是线性可分时，SVM通过硬间隔最大化学习进行分类。当训练数据线性不可分时，SVM通过使用核技巧以及软间隔最大化学习进行分类。SVM本文档来自技高网...

【技术保护点】
1.一种特征mRNA表达谱组合，其特征在于，包括AC011462.1、ARHGEF10L、BMP1、CCM2、CD276、COLGALT1、DCBLD1、GPD1L、GPT2、HOMER3、MPC1、MRGBP、P3H1、PLOD3、PRADC1、SERPINH1、SLC26A6、SMDT1、SNAI2和TPT1，其核苷酸探针序列如SEQ ID NO.1-20所示。/n

【技术特征摘要】
1.一种特征mRNA表达谱组合，其特征在于，包括AC011462.1、ARHGEF10L、BMP1、CCM2、CD276、COLGALT1、DCBLD1、GPD1L、GPT2、HOMER3、MPC1、MRGBP、P3H1、PLOD3、PRADC1、SERPINH1、SLC26A6、SMDT1、SNAI2和TPT1，其核苷酸探针序列如SEQIDNO.1-20所示。


2.一种基于权利要求1所述的特征mRNA表达谱组合的头颈鳞状细胞癌早期预测方法，其特征在于，包括以下步骤：
步骤1、获取头颈鳞状细胞癌早期患者稳定差异表达的特征mRNA；
步骤2、选取特征mRNA表达数据，对每个样本进行数据标准化；
步骤3、使用支持向量机对标准化后的数据构建早期预测模型；
步骤4、根据患者特征mRNA的表达水平进行早期预测；
该方法为非疾病的诊断和治疗目的。


3.根据权利要求2所述的头颈鳞状细胞癌早期预测方法，其特征在于，所述步骤1中的获取头颈鳞状细胞癌早期患者稳定差异表达的特征mRNA，具体为：
步骤1.1、从GenomicDataCommonsDataPortal数据库中下载头颈鳞状细胞癌患者肿瘤组织和癌旁组织转录组数据以及临床数据，获得头颈鳞状细胞癌患者肿瘤组织基因表达谱readcounts数值，即为测序读段数值，进行对数转换；
步骤1.2、选取具有一定表达丰度的mRNA，即在所有样本中mRNA的readcounts大于等于10；再对所有mRNA的readcounts取对数，设样本总数为n，筛选后mRNA总数为m，v为mRNA的readcounts，u为取对数之后的表达值，则有；



其中，i为样本编号，j为mRNA编号，uij为第i个样本、第j个mRNA编号取对数之后的表达值，vij为第i个样本、第j个mRNA编号的readcounts数值；
步骤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检验公式为：



其中n1为肿瘤样本数，n2为正常样本数，μ1为肿瘤样本mRNA表达均值，μ2为正常样本mRNA表达均值，为肿瘤样本mRNA方差，为正常样本mRNA方差；
对所有t检验得出的p值进行错误发现率(falsediscoveryrate,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)。

【专利技术属性】
技术研发人员：刘斐，贺轲，李文兴，安三奇，
申请(专利权)人：广东省第二人民医院广东省卫生应急医院，
类型：发明
国别省市：广东;44