[发明专利]一种isomiR分子标志物的筛选方法及装置

权利要求书

1.一种isomiR分子标志物的筛选方法，所述方法至少包括以下步骤：

S1：利用数据库，获取待检测疾病样本和对应健康样本的isomiR的表达数据；

S2：对所述待检测疾病样本和对应健康样本的isomiR的表达数据进行差异表达分析，得到以对应健康样本为对照，待检测疾病样本的isomiR的变化率；

S3：根据变化率选取isomiR；

S4：对S3中选取的isomiR进行靶基因预测，获得靶基因；

S5：对所述靶基因进行功能分析，筛选与待检测疾病相关的靶基因，所述靶基因及其对应的isomiR即为所述待检测疾病的分子标志物。

2.如权利要求1所述的isomiR分子标志物的筛选方法，其特征在于，还包括以下特征中的一项或多项：

a.步骤S1中，所述数据库选自GEO数据库，和/或，获取isomiR的表达数据发方法包括以下步骤：获取待检测疾病样本和对应健康样本的small RNA测序原始数据，并进行比对、定量及注释，得到isomiR的表达数据；

b.步骤S2中，所述差异表达分析方法选自倍数法，优选的，所述倍数法选取1.5倍或者2倍的差异倍数；

c.步骤S3中，根据变化率选取isomiR是指，在步骤S2得到的结果进行变化率排名，选取排名在前的isomiR；

d.步骤S4中，利用靶基因预测网站或预测软件作为预测靶基因，优选的，所述预测软件选自miRanda；

e.步骤S5中，所述功能分析包括通路分析和/或调控网络的构建；

f.所述待检测疾病选自肿瘤疾病。

3.如权利要求2所述的isomiR分子标志物的筛选方法，其特征在于，还包括以下特征中的一项或多项：

g.特征a中，步骤S1中，获取待检测疾病样本和对应健康样本的small RNA测序原始数据的方法包括如下步骤：

S1.1，获取待检测疾病样本和对应健康样本的isomiR的测序数据；

S1.2，获得原始测序数据的SRA数据链接；

S1.3，利用aspera批量下载所需的测序原始数据；

S1.4，将上述步骤得到的数据转变成fastq格式；

h.特征b中，所述差异表达分析方法选用Benjamini-Hochberg方法、FDR方法或者Bonforroni方法校正P-value后得到差异表达的isomiR；

i.特征e中，基于clusterProfiler进行通路分析；

j.特征f中，所述肿瘤疾病选自骨髓瘤。

4.如权利要求3所述的isomiR分子标志物的筛选方法，其特征在于，特征i中，利用KEGG数据库进行通路分析。

5.如权利要求1所述的isomiR分子标志物的筛选方法，其特征在于，步骤S2中，在差异表达分析之前，还包括对所述待检测疾病样本和对应健康样本的isomiR的表达数据进行过滤。

6.一种isomiR分子标志物的筛选装置，其特征在于，所述装置至少包括：

获取模块，用于利用数据库，获取待检测疾病样本和对应健康样本的isomiR的表达数据；

差异表达分析模块，用于对所述待检测疾病样本和对应健康样本的isomiR的表达数据进行差异表达分析，得到以对应健康样本为对照，待检测疾病样本的isomiR的变化率；

选取模块，用于根据变化率选取isomiR；

靶基因预测模块，用于对选取模块中选取的isomiR进行靶基因预测，获得靶基因；

功能分析模块，用于对所述靶基因进行功能分析，筛选出与待检测疾病相关的靶基因，所述靶基因及其对应的isomiR即为所述待检测疾病的分子标志物。