[发明专利]基于转录组的肿瘤新抗原鉴定方法

摘要

本发明公开了基于转录组的肿瘤抗原鉴定方法，包括获取患者肿瘤组织的RNA样本，对RNA样本建库和扩增，获得肿瘤组织的RNA样本测序结果；将RNA样本测序结果的短读段比对到人类参考基因组，获得RNA比对结果；根据RNA比对结果计算基因表达量，根据RNA比对结果进行突变检测和预测融合基因事件；根据比对结果预测转录组HLA分型；计算基因表达量、突变检测和预测融合基因事件；将转录组样本的基因表达量，转录组突变位点在全外测序样本中的深度，新生短肽与患者HLA分型的结合力作为分析结果交给下游分析人员四个步骤。本发明提供了一种能够从肿瘤患者转录组NGS数据中鉴定出个体样本的肿瘤特异性抗原的方法。

主权项

1.基于转录组的肿瘤抗原鉴定方法，其特征在于，包括以下步骤：S1：获取患者肿瘤组织的RNA样本，对RNA样本建库和扩增，获得肿瘤组织的RNA样本测序结果，根据测序结果鉴定患者的HLA分型；S2：将RNA样本测序结果的短读段比对到人类参考基因组，获得RNA比对结果；S3、根据RNA比对结果计算基因表达量，根据RNA比对结果进行突变检测，根据RNA比对结果预测融合基因事件；根据比对结果预测转录组HLA分型；计算基因表达量、突变检测和预测融合基因事件按指定顺序进行，或者同步进行；S3A：计算基因表达量包含：去除RNA比对结果中的内含子短读段，形成仅由外显子短读段组成的RNA比对结果，计算基因表达量，计算基因表达量包含以下步骤：S3A1、去除建库过程中由PCR扩增产生的重复的短读段序列；S3A2、将比对到内含子的短读段去除,获得由外显子短读段组成的RNA比对结果；S3A3、计算每个外显子短读段数，用以下公式换算成基因的表达量：其中total exon reads表示比对到某个外显子区域的总短读段数，mapped reads(millions)表示每1000000条短读段中比对到该区域的短读段数，exon length(KB)表示这段外显子的长度；S3B：突变检测包含：S3B1：从RNA比对结果中检测RNA样本中的突变，去掉正常组织的全外显子样本中具有的突变，剩余的突变作为RNA水平的体细胞突变；对RNA的突变进行功能注释；S3B2：计算突变位点的深度：针对每个突变位点，计算当前突变位点在全外显子测序样本中的深度，全外显子样本包括肿瘤组织的全外显子样本和正常组织的全外显子样本；S3B3：HLA分型亲和力预测：根据突变的注释结果，针对每一个突变事件编辑相应转录本序列，获得一条包含所有突变的多肽序列，去掉在野生型蛋白序列中可以匹配到的肽段形成转录组特有的新生多肽序列，将转录组特有的多肽序列按照HLA分型结合需要的长度截取成短肽，将短肽与步骤S1的RNA样本测序结果中鉴定的患者HLA分型做亲和力预测；其中，计算突变位点的深度和HLA分型亲和力预测按指定顺序进行或同步进行；S3C、预测融合基因事件：根据步骤S2的RNA比对结果预测转录组的融合基因事件，从中筛选出可信的融合基因后，翻译成蛋白序列，将蛋白序列按照HLA分型结合需要的长度截取成短肽；将短肽与步骤S1的RNA样本测序结果中鉴定的患者HLA分型做亲和力预测；步骤S3B3和S3C中截取短肽的方法为：SI、顺序获取多肽序列上的突变为当前突变，以当前突变为中心、向前截取n个氨基酸和向后截取m个氨基酸获得多肽序列，n为HLA分型所能呈递的最大长度，m为HLA分型所能呈递的最大长度、或者m为从当前突变到第一个终止密码子的长度；SII，在多肽序列上依次截取长度为N的短肽，N为HLA分型结合需要的长度；获得N条包含当前突变的短肽；S4：将转录组样本的基因表达量，转录组突变位点在全外测序样本中的深度，新生短肽与患者HLA分型的亲和力作为分析结果交给下游分析人员。