[发明专利]分析生物基因组基因表达、拷贝数变异的可视化方法

摘要

本发明公开了一种分析生物基因组基因表达、拷贝数变异的可视化方法，每个染色体的所有基因，按照每个基因的起始位置作为横坐标，用matlab软件中的独立样本检验函数ttest2，每个基因会得到一个P值，然后对P值进行处理，然后规定两组样本的颜色，确定Y轴朝向，用matlab中的line函数绘制线，确定百分比值。本发明的可视化方法，在两种癌症类型下进行基因表达量或者拷贝数变异对比的时候，或者在同一个癌症类型下正常类型样本和肿瘤类型样本的基因表达量或者拷贝数变异对比的时候可以直观的反映哪种类型的癌症或者哪种类型样本的值大；对于基因表达量，高级曼哈顿图可以直观显示出基因是过表达呢还是表达不足；对于拷贝数变异，高级曼哈顿图可以直观显示出基因是扩增还是丢失。

主权项

一种分析生物基因组基因表达、拷贝数变异的可视化方法，其特征在于，包括以下步骤：(1)横纵坐标的确定：把每个染色体的所有基因，按照每个基因的起始位置作为横坐标；用matlab软件中的独立样本检验函数ttest2，然后每个基因会得到一个P值，然后对P值进行处理，‑10log10(P)这个值始终是正的，我们将这个值作为纵坐标的幅度值；(2)颜色的确定：事先需要规定两组样本的颜色，然后用matlab中的中位数函数，遍历出这两组样本中每个基因的拷贝数的中位数值，然后比较每个基因的两个中位数，谁的中位数值大，则颜色就取规定的该中位数所属组的颜色；(3)Y轴朝向的确定：在找到每个基因拷贝数中位数值大的所属的样本组后，拿这个大的基因拷贝数值中位数和该组所属的癌症类型中的正常样本的基因拷贝数中位数值做差，做差后的结果，如果值是正值则纵轴方向朝上，如果值是负值则纵轴方向朝下；(4)绘制线：用matlab中的line函数；(5)垂直虚线：因为中心体也有具体的位置，同样可以用line函数，把虚线绘制出来；(6)百分比值的确定：用Bonferonni校正的p值为2*10的‑6次，做水平虚线，然后计算每个染色体中的p臂或者q臂中，P值大于那条水平虚线阈值的基因个数占该染色体整个p臂或者q臂的比例。