本发明涉及生物
技术领域：
：，尤其是涉及癌症特异性tcr及其分析技术，以及所述癌症特异性tcr在癌症治疗中的应用。
背景技术：
：高等生物防御机制高度依赖于包含t细胞和b细胞在内的过继免疫系统。t细胞具有在细胞表面上表达的特异性的受体分子tcr(tcellreceptor)，其能识别并且区分自身或外源抗原。通过由抗原受体反应来传递细胞内信号，促进细胞增殖，进而启动各种免疫应答，如增强炎性细胞因子、趋化因子等的产生。tcr识别与抗原呈递细胞表达的主要组织相容性复合体(mhc)以及抗原肽结合，形成肽-mhc(pmhc)-tcr复合体，将自身和非自身的抗原区分开并且识别抗原肽。tcr是由两个tcr多肽链组成的异源二聚体受体分子，大多数t细胞表达αβtcr，以及少数t细胞表达具有特殊功能的γδtcr。α链和β链tcr分子与多个cd3分子(cd3ζ链、cd3ε链、cd3γ链和cd3δ链)形成复合体，在抗原识别之后传递细胞内信号，并且启动多种免疫应答。内源性抗原，诸如来自癌细胞的癌抗原或在细胞内增殖的病毒抗原，被mhci类分子递呈为抗原肽。抗原呈递细胞通过内吞作用获得并且加工来源于外源性微生物的抗原，然后递呈在mhcii类分子上。这样的抗原被cd8+t细胞和cd4+t细胞各自所表达的tcr识别。tcr基因由基因组中不同区域编码的多个v区(可变区，v)、j区(连接区，j)、d区(多变区，d)和c区(恒定区，c)组成。在t细胞分化过程中，这些基因片段以多种组合进行基因重排。α链和γ链tcr表达由v-j-c组成的基因，而β链和δ链tcr表达由v-d-j-c组成的基因。目前，imgt(internationalimmunogeneticsproject)的数据库具有43种功能性α链tcrv基因片段(trav)、50种tcrj基因片段(traj)，40～42种功能性β链tcrv基因片段(trbv)、2种tcrd基因片段(trbd)、13种tcrj基因片段(trbj)，4～6种功能性γ链tcrv基因片段(trgv)、5种tcrj基因片段(trgj)，3种功能性δ链tcrv基因片段(trdv)、3种tcrd基因片段(trdd)和4种tcrj基因片段(trdj)(nucleicacidresearch，2009，37(suppl1)，d1006-d1012.)。通过这些基因片段的重排产生多样性。此外，在v基因片段和d基因片段之间，或在d基因片段和j基因片段之间，插入或缺失一个或多个碱基，导致随机氨基酸序列的形成，从而产生更多多样性的tcr基因序列。tcr分子与pmhc复合体表面直接结合的区域(tcr足迹)由v区内的三个多样性互补确定区(cdr)cdr1区、cdr2区和cdr3区构成。cdr3区具体包括一部分v区、一部分j区和由随机序列形成的v-d-j区，形成最多样的抗原识别位点。同时，其它区域被称为fr(构架区)，用于形成tcr分子的骨架结构。胸腺中的t细胞在分化和成熟过程中，β链tcr最先进行基因重排，并且缀合ptα分子以形成前tcr复合体分子。然后，α链tcr进行重排形成αβtcr分子，并且当没有形成功能性αβtcr时，则在其他α链tcr基因等位基因中发生重排。已知在胸腺中进行正/负选择之后，选择出具有适当亲和力的tcr以获得抗原特异性。t细胞产生对特定抗原具有高水平特异性的一种tcr。由于在活的生物体中有许多抗原特异性t细胞，可形成多样性tcr库，以有效起到抵御多种病原体的防御机制的作用，tcr库是免疫细胞的特异性或多样性的重要指示物。对tcr库的分析是提升免疫反应效果或治疗自身免疫性疾病的有用方法。如果t细胞应答于抗原而进行增殖，则在多样性库中观察到特定tcr基因的比率增加(克隆性增加)。已经尝试通过对tcr库进行分析，以根据克隆性的增加来检测肿瘤中表达tcr的淋巴细胞的发育(leukemiaresearch，2003，27，305-312)；报道了在暴露于选择性刺激具有特定vβ链的tcr的分子(诸如超抗原)时，特定vβ链的使用频率增加(immunology1999，96，465-72.)。为了研究抗原特异性免疫应答，频繁用于分析由免疫失调诱导的难治的自体免疫性疾病，诸如类风湿性关节炎、系统性红斑狼疮、肖格伦综合征和特发性血小板减少性紫癜，并且已证明了其有用性。当前免疫治疗已经成为肿瘤临床治疗中不可或缺的环节。免疫治疗的药物和方案涉及到机体免疫系统识别和攻击癌细胞的各个阶段。已有的肿瘤免疫药物包括以下多个类型：靶向癌细胞的抗体、过继细胞治疗、溶瘤病毒、树突状细胞相关治疗、dna和蛋白水平的肿瘤疫苗、免疫激活细胞因子以及其他免疫调节化合物。其中针对t细胞检验点抑制蛋白的抗体类药物和肿瘤抗原特异性的t细胞过继疗法近年来取得突破，广受瞩目。过继细胞疗法(adoptivecelltransfer，act)，是将激活的t淋巴细胞输入机体行使免疫功能。一般操作方式为：从病患肿瘤组织中分离肿瘤浸润淋巴细胞群，从中分离出t细胞并辅以t细胞生长因子(如il-2)在体外培养和激活；筛选出具有肿瘤特异性的t细胞进行体外扩大培养，随后输入患者体内，联合化疗或放疗共同治疗。act最具挑战的步骤是t细胞特异地识别肿瘤细胞。嵌合抗原受体-t细胞(chimericantigenreceptort-cell,car-t)疗法和tcr-t疗法是提高其效率的有效途径。car-t需要构建嵌合抗原受体，通常是用编码两个抗体可变区的基因编码一个可变区，再嫁接到tcr的胞内部分，后者可以激活t细胞的免疫活性。而tcr-t疗法则用人类肿瘤抗原刺激表达人类mhcⅰ的小鼠，从而得到特异杀死肿瘤的小鼠t细胞，克隆其tcr并使其在病患t细胞上表达，最终将这些改造好的细胞输回病患体内进行免疫治疗。因此，如何高效分离和分析患者个体化的tcr，并且开发可用于肿瘤治疗的tcr-t细胞有重要的临床价值。技术实现要素：本发明发明人利用单细胞转录组分析技术，通过分析癌症患者的癌症组织中浸润的t细胞，发现、分离并表征了一系列新的克隆性的tcr基因及其序列，这些tcr可能是由肿瘤细胞抗原刺激t细胞产生的，表达这些tcr的t细胞可能具有特异性识别肿瘤细胞抗原，并杀死肿瘤细胞的活性，存在用于过继细胞疗法的前景。进一步的，本发明发明人提供了一种利用单细胞tcr序列进行柔性对接，进而预测tcr，mhc和小肽段的结合能力，并用于预测或发现癌症患者体内新的肿瘤抗原的方法。本发明的一个目的是提供一系列新的tcr。本发明的再一个目的是提供编码这些tcr的核酸序列。本发明的另一个目的是提供携带所述tcr的核酸序列的表达载体。本发明的另一个目的是提供表达所述tcr的t细胞。本发明的另一个目的是提供所述t细胞的制备方法，以及相应t细胞在过继免疫治疗中的应用。本发明的另一个目的是提供一种新的t细胞的单细胞转录组分析方法。本发明的另一个目的是提供一种预测tcr，mhc和小肽段结合能力的计算方法。本发明的再一个目的是提供一种预测新的肿瘤抗原的方法。本发明的技术方案如下：一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表1-46中的每张表所列的tcr的α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表1-46中的每张表所列的对应tcr的α链的vj和β链的vdj。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别如表1-46中的每张表所列的对应tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列，或实质与之相似的序列。优选所述tcr的α链的全长核酸序列和β链的全长核酸序列分别如表1-46中的每张表所列的对应tcr的α链的全长核酸序列和β链的全长核酸序列，或与之基本上同源的核酸序列。具体如下：一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表1中的seqidno.1和seqidno.3，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表1中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表1中的seqidno.2和seqidno.4，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表2中的seqidno.5和seqidno.7，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表2中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表2中的seqidno.6和seqidno.8，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表3中的seqidno.9和seqidno.11，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表3中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表3中的seqidno.10和seqidno.12，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表4中的seqidno.13和seqidno.15，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表4中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表4中的seqidno.14和seqidno.16，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表5中的seqidno.17和seqidno.19，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表5中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表5中的seqidno.18和seqidno.20，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表6中的seqidno.21和seqidno.23，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表6中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表6中的seqidno.22和seqidno.24，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表7中的seqidno.25和seqidno.27，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表7中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表7中的seqidno.26和seqidno.28，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表8中的seqidno.29和seqidno.31，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表8中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表8中的seqidno.30和seqidno.32，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表9中的seqidno.33和seqidno.35，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表9中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表9中的seqidno.34和seqidno.36，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表10中的seqidno.37和seqidno.39，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表10中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表10中的seqidno.38和seqidno.40，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表11中的seqidno.41和seqidno.43，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表11中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表11中的seqidno.42和seqidno.44，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表12中的seqidno.45和seqidno.47，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表12中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表12中的seqidno.46和seqidno.48，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表13中的seqidno.49和seqidno.51，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表13中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表13中的seqidno.50和seqidno.52，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表14中的seqidno.53和seqidno.55，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表14中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表14中的seqidno.54和seqidno.56，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表15中的seqidno.57和seqidno.59，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表15中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表15中的seqidno.58和seqidno.60，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表16中的seqidno.61和seqidno.63，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表16中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表16中的seqidno.62和seqidno.64，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表17中的seqidno.65和seqidno.67，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表17中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表17中的seqidno.66和seqidno.68，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表18中的seqidno.69和seqidno.71，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表18中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表18中的seqidno.70和seqidno.72，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表19中的seqidno.73和seqidno.75，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表19中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表19中的seqidno.74和seqidno.76，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表20中的seqidno.77和seqidno.79，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表20中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表20中的seqidno.78和seqidno.80，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表21中的seqidno.81和seqidno.83，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表21中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表21中的seqidno.82和seqidno.84，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表22中的seqidno.85和seqidno.87，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表22中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表22中的seqidno.86和seqidno.88，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表23中的seqidno.89和seqidno.91，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表23中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表23中的seqidno.90和seqidno.92，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表24中的seqidno.93和seqidno.95，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表24中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表24中的seqidno.94和seqidno.96，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表25中的seqidno.97和seqidno.99，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表25中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表25中的seqidno.98和seqidno.100，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表26中的seqidno.101和seqidno.103，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表26中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表26中的seqidno.102和seqidno.104，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表27中的seqidno.105和seqidno.107，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表27中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表27中的seqidno.106和seqidno.108，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表28中的seqidno.109和seqidno.111，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表28中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表28中的seqidno.110和seqidno.112，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表29中的seqidno.113和seqidno.115，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表29中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表29中的seqidno.114和seqidno.116，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表30中的seqidno.117和seqidno.119，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表30中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表30中的seqidno.118和seqidno.120，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表31中的seqidno.121和seqidno.123，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表31中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表31中的seqidno.122和seqidno.124，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表32中的seqidno.125和seqidno.127，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表32中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表32中的seqidno.126和seqidno.128，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表33中的seqidno.129和seqidno.131，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表33中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表33中的seqidno.130和seqidno.132，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表34中的seqidno.133和seqidno.135，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表34中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表34中的seqidno.134和seqidno.136，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表35中的seqidno.137和seqidno.139，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表35中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表35中的seqidno.138和seqidno.140，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表36中的seqidno.141和seqidno.143，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表36中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表36中的seqidno.142和seqidno.144，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表37中的seqidno.145和seqidno.147，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表37中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表37中的seqidno.146和seqidno.148，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表38中的seqidno.149和seqidno.151，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表38中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表38中的seqidno.150和seqidno.152，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表39中的seqidno.153和seqidno.155，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表39中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表39中的seqidno.154和seqidno.156，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表40中的seqidno.157和seqidno.159，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表40中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表40中的seqidno.158和seqidno.160，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表41中的seqidno.161和seqidno.163，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表41中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表41中的seqidno.162和seqidno.164，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表42中的seqidno.165和seqidno.167，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表42中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表42中的seqidno.166和seqidno.168，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表43中的seqidno.169和seqidno.171，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表43中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表43中的seqidno.170和seqidno.172，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表44中的seqidno.173和seqidno.175，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表44中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表44中的seqidno.174和seqidno.176，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表45中的seqidno.177和seqidno.179，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表45中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表45中的seqidno.178和seqidno.180，或实质与之相似的序列。一种tcr，其α链的cdr3的氨基酸序列和β链的cdr3的氨基酸序列分别为表46中的seqidno.181和seqidno.183，或实质与之相似的序列。优选所述tcr的α链的vj和β链的vdj分别如表46中对应所述。优选所述tcr的α链的v区的氨基酸序列和β链的v区的氨基酸序列分别为表46中的seqidno.182和seqidno.184，或实质与之相似的序列。分离的核酸，其编码上述tcr的α链的cdr3的氨基酸序列，或者β链的cdr3的氨基酸序列，或者α链的v区的氨基酸序列，或者β链的v区的氨基酸序列，或者α链的vj，或者β链的vdj，或者与所述氨基酸序列实质上相似的氨基酸序列。分离的核酸，其编码如表1-46中所述的任一tcr的α链的全长核酸序列，或任一tcr的β链的全长核酸序列，或与之基本上同源的核酸序列。表达载体，包含所述的核酸。根据本发明，所述载体包括但不限于病毒、质粒、粘粒、噬菌体、酵母等。宿主细胞，包含所述的核酸。根据本发明，所述宿主细胞包括但不限于真核细胞、细菌细胞、昆虫细胞或人细胞。例如：vreo细胞、hela细胞、cos细胞、cho细胞、hek293细胞、bhk细胞、mdkii细胞、sf9细胞等。一种具有上述tcr或编码所述tcr的核酸序列的t细胞群、t细胞株或重组表达的t细胞。所述的t细胞群、t细胞株、重组表达的t细胞或编码该tcr的核酸序列在诊断或治疗中是有用的。对于诊断，通过检验是否上述序列仅在癌症的患者中，是否上述序列更多地在癌症患者中观察到，或者上述序列在癌症患者的癌症组织中积累，可发现癌症，或者可预测病理状况或预后。对于癌症的治疗，可以利用具有上述tcr的t细胞群，t细胞株，或重组表达的t细胞。在本发明的一个实施方式中，所述癌症包括肝癌，尤其是肝细胞癌。上述tcr、分离的核酸、表达载体、t细胞群、t细胞株或重组表达的t细胞在制备治疗癌症的药物中的应用。在本发明的一个实施方式中，所述癌症包括肝癌，尤其是肝细胞癌。一种制备含有所述tcr的t细胞的方法。在本发明的一个实施方式中，所述制备方法可以包括如下步骤：(1)基于本发明所述的tcr，确定候选hla和测试肽的氨基酸序列；(2)合成所确定的hla和测试肽，并体外形成复合体；(3)用该hla-肽刺激淋巴细胞。根据本发明，可以基于使用hla结合肽预测算法计算的得分，进行该候选hla-肽的氨基酸序列的确定。例如可以使用bimas、syfpeithi、rankpep或netmhc等，确定所述候选hla-肽。在本发明的另一个实施方式中，所述制备方法可以包括以下步骤：(1)将本发明所述的tcrα或tcrβ基因引入到用于基因表达的逆转录病毒载体中；(2)从表达tcrα和tcrβ基因的逆转录病毒载体，创建引入基因的病毒；(3)使用所述带有tcrα和tcrβ基因的病毒独立且依次地感染从患者收集的淋巴细胞，以进行转染，或者创建包括tcrα和tcrβ基因的基因表达逆转录病毒载体，以一次转化这两种基因；(4)证明tcrα/tcrβ异源二聚体在细胞表面上表达。一种t细胞的单细胞转录组tcr分析方法，所述方法包括如下步骤：(1)获得单个的t细胞；(2)构建每个t细胞的cdna文库并测序，获得每个细胞每个基因的表达量；(3)鉴别单个t细胞的tcr序列和克隆性识别。根据本发明，可采用本领域已知的多种方法获得单个的t细胞，例如对于血液中的单个t细胞，可采用密度梯度离心法；对于组织中的单个t细胞，可采用研磨法。根据本发明，可采用本领域已知的各种单细胞转录组的cdna文库构建方法，构建每个t细胞的cdna文库并测序，获得每个细胞每个基因的表达量，例如：汤富酬2009年创立的方法(tang,f.etal.rna-seqanalysistocapturethetranscriptomelandscapeofasinglecell.nat.protoc.5,516–535(2010))、strt-seq(single-celltaggedreversetranscriptionsequencing)、smart-seq和smart-seq2、cell-seq(cellexpressionbylinearamplificationandsequencing)和pma-seq(phi29-mrnaamplificationandsequencing)等方法。在本发明的一个优选实施方式中，采用smart-seq2构建每个t细胞的cdna文库并测序，获得每个细胞每个基因的表达量。本发明的发明人，对比研究了汤富酬2009年创立的方法和smart-seq2，发现smart-seq2法在保证测序质量的情况下可以检测到更多的基因，包括t淋巴细胞共有的标志物cd3基因；并且smart-seq2法更有利于扩增出完整的cdna，更适用于进行t细胞单细胞转录组扩增。本发明的发明人通过实验研究，进一步改进了smart-seq2中的操作条件，提高了mrna的反转录收率，以及pcr扩增后产物的纯化效率。在本发明的具体实施方式中，采用smart-seq2法反转录时，采用如下的反转录条件，提高了反转录cdna的产率和cdna全长的比例：相比于常用的50℃30分钟的反转录条件，改进后的反转录条件能将cdna产率提高16～23％，cdna全长的平均长度提高20％左右。在本发明的具体实施方式中，采用smart-seq2法对pcr扩增产物纯化时的方法如下，提高了pcr产物的纯度，有利于后续的测序和文库构建质量的提高：用磁珠进行两次纯化，第一次纯化时，所加入的磁珠体积与pcr扩增产物体积相同，第二次纯化时，所加入的磁珠体积为pcr扩增产物体积的2倍。根据本发明，在进行步骤(3)的分析时，对步骤(2)获得的生物信息数据进行比对和质量控制，去除低质量的部分。根据本发明，对于cdna的测序读段(reads)数据质量控制的方法为：保留符合以下条件的测序读段：①未知碱基占给定读段总序列不超过10％，②phred质量值低于5的碱基不超过50％，③不含有接头序列。根据本发明，对于细胞质量控制的方法为，去除数据量和数据质量低的细胞，保留符合以下条件的细胞：①cd3d的tpm大于3；②当分离cd4+t细胞时，cd4的tpm需要大于3，同时cd8的tpm小于30；③分离cd8+t细胞时，cd8的tpm需要大于3，同时cd4的tpm小于30；④线粒体基因上的读段占所有读段的比例不高于10％。其中，tpm值的定义为：其中cij表示为基因i在细胞j中的读段的数量。根据本发明，对用于分析的单细胞的基因表达量的质量控制方法为：当一个基因在所有细胞中被检测到的读段平均数量大于1才用于后续分析。根据本发明，在步骤(3)中使用软件tracer进行单个t细胞的tcr序列识别。根据本发明，在步骤(3)的克隆性识别中，采用如下方法：比较任意两个细胞中的tcrɑ和tcrβ的序列，当两个细胞中至少一种tcrɑ同时至少一种tcrβ的序列完全一致时，且一致的tcrɑ和tcrβ的序列是可以翻译成有效蛋白质的，且tcrɑ的tmp值至少大于10和tcrβ的tmp值至少大于15，这样的两个细胞被认为来自同一克隆。一种预测tcr，mhc和小肽段结合能力的计算方法，其特征在于包括如下步骤：1)获得癌症患者肿瘤免疫细胞的tcr的rna序列，肿瘤患者的mhc类型，以及小肽段的序列，输入rosettadock软件；2)根据已知序列和蛋白质结构数据库，对tcr序列进行蛋白质结构的同源建模；3)确认tcr中cdr的6个loop区(环状区域)，并进行分步模拟，计算出所述6个loop区的结合自由能；4)将mhc，tcr和小肽段结合在一起，分别进行低分辨率以及高分辨率的对接进程计算，达到最大迭代次数终止计算；5)分析结果，rmsd，计算对接自由能和表示结合能力强弱的打分函数值(rosettascore)。根据本发明，步骤1)中，癌症患者肿瘤免疫细胞的tcr的rna序列可以来自于已知各种公共数据库中已经披露的tcr的rna序列，也可以是采用本领域已知的各种测序方法由癌症患者肿瘤免疫细胞中测序获得的。优选，采用本发明提供的t细胞的单细胞转录组tcr分析方法，对采自癌症患者的肿瘤免疫细胞进行单细胞测序和分析，获得大量潜在可用的tcr的rna序列，用于本发明的计算和预测分析。根据本发明，步骤1)中，患者的mhc类型，可以采用本领域已知的外显子测序方法并运行optitype得到mhc类型，例如参照szoleka1,schubertb2,mohrc2,sturmm1,feldhahnm1,kohlbachero1.optitype:precisionhlatypingfromnext-generationsequencingdata.bioinformatics.2014dec1；30(23):3310-6.doi:10.1093/bioinformatics/btu548.epub2014aug20.记载的方法进行实验和分析。根据本发明，步骤1)中，小肽段序列可以采用netmhc和本领域已知的rna测序技术，预测患者体内的小肽段序列，例如，参照andreattam,nielsenm.gappedsequencealignmentusingartificialneuralnetworks:applicationtothemhcclassisystem.bioinformatics(2016)feb15；32(4):511-7；nielsenm,lundegaardc,worningp,lauemollersl,lamberthk,buuss,brunaks,lundo.reliablepredictionoft-cellepitopesusingneuralnetworkswithnovesequencerepresentations.proteinsci.,(2003)12:1007-17记载的方法进行实验和分析。根据本发明，步骤2)中同源建模可采用本领域常用的各种同源建模方法，将得到的tcr的rna序列翻译成氨基酸序列，根据已知序列和蛋白质结构数据库，通过寻找同源蛋白，预测tcr的三维结构。根据本发明，步骤3)中，根据小肽段的氨基酸残基计算出小肽段主链中心，依据rosettadock软件识别出tcr的cdr的loop区，计算各loop区与小肽段主链中心的距离，选择距离最近的6个loop区作为后续分析tcr，mhc和小肽段结合能力的基础。所述小肽段主链中心是指，小肽的主链碳原子的坐标的中点，即，该小肽段第一个氨基酸残基的主链碳原子的坐标和最后一个氨基酸残基的主链碳原子的坐标的平均值。小肽段主链中心的计算及其方法是本领域已知的。根据本发明，步骤3)中，进行分步模拟时，可以释放所有6个loop区；也可以为减少可变结构域对计算结果的干扰，每次模拟只释放6个loop区中的5个、4个、3个、2个或1个，相应固定剩余的1个、2个、3个、4个或5个，如此类推，计算出每个loop区的结合自由能。优选每次模拟只释放6个loop区中的1个，固定剩余的5个，如此类推，计算出每个loop区的结合自由能。根据本发明，步骤4)中首先通过低分辨率搜索探索同源建模得到的构象空间，然后通过蒙特卡罗最小化算法对所有原子进行局部细化。低分辨率对接时，蛋白质表示为主干加上侧链的质心表示，即侧链被表示为一个巨大的原子，以节省cpu时间。在这个阶段，rosettadock试图为高分辨率搜索找到对接对象的大致方向。进行高分辨率对接(局部细化)时，蛋白质中的所有原子都被表示出来，在低分辨率搜索中找到的位置也得到了优化。高分辨率阶段消耗了rosettadock最多的cpu时间。根据本发明，依据步骤5)中获得的评分函数可以预测该患者体内最有可能引起免疫应答的新的抗原。rosettascore越低的构象代表能量状态越低，能量状态越低的tcr-mhc-小肽段结构，越有可能是真实存在于生物体内的结合构象，由此，参与形成这种构象的小肽段越有可能是最易与所述tcr结合的小肽，即越有可能是引起免疫应答的新抗原，而相应的tcr序列也越有可能是结合力最强的tcr，可用于tcr-t疗法。因此，所述计算方法还可用于预测或筛选新的肿瘤抗原，和/或预测或筛选可用于未来开发的tcr序列，提高从高通量测序技术获得的大量的肿瘤免疫细胞tcr序列中得到具有潜在生物活性价值的tcr的效率。在本发明的一个实施方式中，所述患者是肝癌患者，本发明预测tcr，mhc和小肽段结合能力的计算方法，以及预测或筛选新的肿瘤抗原，和/或预测或筛选可用于未来开发的tcr序列，是针对肝癌患者进行的。为了评估tcr对接基准内的结构多样性，本发明的发明人比较了tcr/pmhc结构中主干构象以及更多的环状柔性结构的区域。结合tcrcdr环的叠加显示了较大的结构变异，尤其是tcr的cdr3α和cdr3β环之间，cdr1α和cdr2α环的结构变异程度较小。在pmhcs的重叠中，肽骨架构象表现出很大的多样性，由不同的肽序列和长度以及mhc等位基因和结合tcr驱动。发明人把cdr结构变异程度大的区域(loop区)，根据与小肽段主链中心的距离远近，选择出距离最近的6个区域，在对接的过程中6个loop区分别进行无约束动态调整，从而达到tcr/pmhc复合物间最适合的对接结果。除了分析基准中tcr/pmhc复合物之间的变化外，发明人还计算了tcrs的结合构象变化(根据每个测试用例计算的无约束和结合tcr结构)作为位置的函数，结果发现cdr3α环在结合时表现出最大的平均构象变化，其次是cdr3β和cdr1α，其他cdr和pmhc结合位点的构象变化不太明显。由此说明，采用上述确定用于计算的loop区的方法有效可行，确定出tcr的cdr的6个loop区就足以进行准确而又快捷的计算。由于tcr-pmhc配合物的生产、结晶和结构测定具有挑战性，因此对新型配合物的建模具有重要的意义。本发明人利用了已知复合物中保守的结构特征，如受限的tcr结合位点和一般保守的对角对接模式，提供了前述的快速tcr-pmhc建模方法。在本发明中：当用来描述核酸时，“同源的”表示当两个核酸或其指定的序列在最佳比对和比较时至少80％的核苷酸、并且更优选至少大约98％至99％的核苷酸是一致的，伴有适当的核苷酸插入或缺失。术语“同源物”或“同源的”也指就结构和/或功能而言的同源性。就序列同源性而言，如果多个序列是至少80％以上相同的，例如：至少90％、至少95％、至少97％或至少99％，则它们是同源物。术语“基本上同源的”是指至少90％以上相同，例如至少95％相同、至少97％相同或至少99％相同的序列。当用来描述多肽序列时，术语“实质相似性”表示这种多肽包含在大约10-100个氨基酸残基的比较窗(例如，抗体的重链或轻链可变区，tcr的α链或β链的v区)范围内与参考序列具有至少80％以上序列一致性、或最优选90％一致性、或最优选95％一致性、或最优选99％一致性的序列。在氨基酸序列的背景下，“实质相似性”进一步包括氨基酸的保守性置换。术语“实质一致性”表示当最佳比对(如通过使用默认空位权重的程序gap或bestfit)时，两个肽序列共享至少80％序列一致性，优选至少90％序列一致性，更优选至少95％或更高的序列一致性(例如，至少99％或更高的序列一致性)。优选地，不相同的残基位置因保守性氨基酸置换而不同。本发明的基因或氨基酸序列的同源物的确定可以由技术人员容易地确定。术语“恶性肿瘤”、“肿瘤”和“癌症”可互换使用，是指以细胞不受控制、过度增生性或异常生长或转移为特征的疾病或失调。如无特殊说明，在本发明中，氨基酸序列都是从n端到c端，碱基序列都是从5’端到3’端。附图说明图1肿瘤组织、正常组织和外周血t细胞的流式细胞分析图。第一行：外周血中t细胞；第二行：正常组织中的t细胞；第三行：肿瘤组织中t细胞。图2合格的单细胞cdna示例。96个细胞同时使用lifetechreal-timepcr仪7500检测所得的结果图。对象rt-pcr的ct值低于26(左图黑框所示)；对象溶解曲线的峰值在85℃到90℃之间(右图黑框所示)。图3合格的单细胞cdna示例。基于毛细管电泳的fragmentanalysis的检测结果。1700左右的峰是全长转录组的片段大小，1100左右的峰是作为内参的ercc。图4构建成的cdna文库的fragmentanalysis测定结果图。图5tang2009的方法的扩增偏向统计图。图6smart-seq2的扩增偏向统计图。图7使用rosettadock进行tcr/pmhc对接模拟的流程图图8ttc35-1116的tcr的cdrloop区区分具体实施方式以下结合实施例对本发明做进一步描述。以下实施例是以肝癌患者为例，对t细胞单细胞转录组分析方法的示例说明。需要说明的是，实施例不能作为对本发明保护范围的限制，本领域的技术人员理解，任何在本发明基础上所作的改进和变化都在本发明的保护范围之内。以下实施例所用化学试剂都是常规试剂，均可商购获得。所用分析软件及其来源如下：gsnap(http://research-pub.gene.com/gmap/)；tracer(https://github.com/teichlab/tracer)；统计软件r(https://www.r-project.org/)。实施例1t细胞的单细胞转录组数据获取1、临床样本收集从2014年8月到2015年6月，在北京大学人民医院和北京世纪坛医院采集患者的手术组织和外周血，包括癌组织(直径3～5mm)和癌旁正常组织，并采集外周血(3ml)。患者为肝细胞癌，未经受术前辅助放疗或化疗，共5例。本项研究符合赫尔辛基宣言的医学伦理标准，并通过北京大学医学伦理委员会的审核。血样在手术前采集于edta抗凝管中暂时于冰上保存；癌组织和癌旁正常组织样本在手术中采集，其中癌症组织剔除坏死组织；癌旁组织为远离癌组织至少5cm处的正常组织。癌组织和癌旁组织在离体30分钟以内置于冰上和rnalater(qiagen)溶液中，当日内完成单细胞分离操作。2、单细胞悬液制备外周血：采用密度梯度离心法分离外周血单核细胞。具体操作为将3ml全血缓慢加到3ml分离液(sigma，cat.no.1077)上，400g室温离心30分钟，小心吸取白色层单核细胞，用10mlpbs清洗，4℃离心15分钟，重复上述清洗步骤一次。最后将细胞溶解于0.5mlpbs，并加入1％小牛血清(fbs)。癌组织和癌旁正常组织：采用研磨法获得癌组织和癌旁正常组织的单细胞。首先将手术离体的组织剪成1mm3大小的碎块，浸泡于rpmi-1640培养基中，并加入10％小牛血清。使用铜网快速磨碎组织，通过40μm筛除组织碎片，400g离心10分钟收集单细胞悬液。使用红细胞裂解液进一步去除组织中混入的红细胞。同样用10mlpbs清洗两次，最后将细胞溶解于0.5mlpbs，并加入1％小牛血清。3、目的t细胞单细胞分离分离的目的细胞包括细胞毒t细胞(cd3阳性，cd8阳性)，辅助性t细胞(cd3阳性，cd4阳性，cd25阴性)和调节性t细胞(cd3阳性，cd4阳性，cd25阳性)。这三种细胞分别用不同抗体进行荧光标记，抗体来自ebioscience公司，每106个细胞使用5μl抗体：兔抗cd3抗体(facs，cat#48-0037-41)兔抗cd4抗体(facs，cat#11-0048-41)小鼠抗cd8抗体(facs，cat#17-0086-41)小鼠抗cd25抗体(facs，cat#12-0259-42)7aad(facs，cat#00-6993-50)，7aad用于标记死亡细胞。96孔板的每个孔中预先加入反应溶液：引物序列为：aagcagtggtatcaacgcagagtacttttttttttttttttttttttttttttttvn目的t细胞的分离结果如图1所示。根据细胞表面的分子标记，选择细胞毒t细胞，辅助性t细胞和调节性t细胞，用流式细胞仪分别将单个细胞收集入相应的96孔板的每个孔中。4、mrna反转录和cdna扩增对于分离至96孔板中的单细胞进行反转录获得cdna，操作步骤按照smart-seq2方法(picelli,s.etal.full-lengthrna-seqfromsinglecellsusingsmart-seq2.nat.protoc.9,171–181(2014).)，如下：1)单细胞裂解：将上述溶液中的单细胞旋涡震荡至少10秒。在pcr仪上72℃孵育3分钟。2)加入内参rna(erccrnaspike-inmix，invitrogen，cat.no.4456740)。事先需要稀释350倍，加入1μl。内参rna有助于对基因表达量的定量计算。3)反转录：反应体系为：tso引物的序列为：aagcagtggtatcaacgcagagtacatrgrg+g反应条件为：4)pcr扩增：反应体系为：kapahifihotstartreadymix(2x)12.5μlispcr引物(10μm)0.25μl无核酸酶的超纯水2.25μl。ispcr引物序列为：aagcagtggtatcaacgcagagt反应条件为：经扩增后的pcr产物用agencourtampurexp磁珠(beckman)进行纯化，方法如下：(1)在25μl上一步反应溶液中加入25μl磁珠，吹打混匀；(2)室温放置5分钟；(3)把盛有溶液的试管或板子放于磁力架上5分钟；(4)移除液体；(5)用100μl80％的乙醇清洗磁珠，放置30秒后移除，重复此过程一次；(6)从磁力架上拿下，加入20μleb溶液，吹吸混匀；(7)放置2分钟后，置于磁力架上，放置2分钟后，再吸出液体。研究中发现，上述过程中残留在溶液中的引物会降低文库构建的效率，使文库包含非细胞cdna的成分。为此需要再增加一遍纯化操作，除将磁珠用量变成50μl之外，纯化过程与前述一致。进行质量检测，即通过rt-pcr检测t细胞特别表达的基因cd3，判断扩增的有效性。反应体系为：cd3的引物序列为：tcattgccactctgctcc(正向)和gttcacttgttccgagcc(反向)。反应条件为：判断cdna是否可用的标准有两条：一是对象rt-pcr的ct值低于26；二是对象溶解曲线的峰值在85℃到90℃之间。图2给出了本实施例获得的合格的cdna示例。另一项质量控制手段为fragmentanalysis，基于毛细管电泳检测样本dna的片段大小和浓度。图3给出了本实施例获得的合格的cdna示例。5、测序文库构建cdna文库构建采用truepreptmdnalibraryprepkitv2for试剂盒(vazyme，cat.no.td501/502/503)；分别匹配双端index为truepreptmindexkitv2for(vazyme，cat.no.td202)。按照试剂盒说明书操作，用1μgcdna起始建库。用磁珠做片段大小的选择，获得目的片段大小为400bp～600bp的cdna文库，最后经fragmentanalysis进行质量控制，经测定，文库构建合格，相应的分析结果参见图4。采用illuminahiseq4000测序，测序模式为双端150bp，通常一个细胞的数据量需要1百万条读段。实施例2生物信息分析1、数据比对与质量控制对于从测序仪获得的读段(reads)，首先去除低质量的部分，保留的标准如下：①未知碱基占给定读段总序列不能超过10％，②质量值低于5的碱基不能超过50％，③不能含有接头序列。使用gsnap软件完成比对。计算基因表达量时，使用的参考基因集合来自ucsc的“knowngene.txt”，使用r语言包“findoverlaps”统计读段在基因上的归属，使用tpm值标定每个基因在每个细胞中的表达量，使用的公式为：其中cij表示为基因i在细胞j中的读段的数量。数据量和数据质量低的t细胞需要被过滤掉。保留符合以下标准的细胞：①cd3d的tpm大于3；②当分离cd4+t细胞时，cd4的tpm需要大于3，同时cd8的tpm小于30；③分离cd8+t细胞时，cd8的tpm需要大于3，同时cd4的tpm小于30；④线粒体基因上的读段占所有读段的比例不高于10％。另外，在文库容量(librarysize)和基因表达数量上也设定了一些参考标准。当一个基因在所有细胞中被检测到的读段平均数量大于1才用于后续分析。2、单细胞tcr序列鉴定和克隆性识别使用软件tracer鉴定单个t细胞中的tcr基因的cdr3部分的序列，重排的方式以及tcr基因的表达量。tracer首先定制一个tcr重排的参考序列库，如对α链重排，将把若干个n碱基(对应v基因的先导序列)、一种v基因段的一种等位基因的序列、连接处的若干个n碱基、一种j基因段的一种等位基因的序列以及c区域的序列连在一起，形成一种可能的重排模板序列；所有v基因的等位基因和j基因的等位基因的共约7,000种组合，形成了α链的参考序列库。类似的方法建立一个β链的参考序列库，包括共约2,272种组合的参考序列库。用比对软件bowtie2(langmeadandsalzberg,2012)将测序读段和这个参考序列库作比对，筛选出来自tcr的读段。bowtie2的参数设置使得其对空位或n碱基的比对有较低的罚分。再用软件trinity(grabherretal.,2011)把筛选出来的tcr的读段组装成重叠群；用igblast将重叠群和imgt(lefrancetal.,2009)数据库比对，识别出重叠群是由哪种v基因的等位基因和(d)j基因的等位基因重排产生的，v(d)j重排的连接处的序列是什么；如果重叠群上从v基因序列开始一直到c区域序列的这段序列有完整的阅读框并且没有终止密码子则这个tcr重排是能翻译成完整蛋白的，进而根据基序识别出cdr3序列。同时，我们对同样的数据采用mixcr和vdjpuzzle进行的tcr序列的鉴定。vdjpuzzle与tracer一样，是一种针对单细胞转录组鉴定tcr序列的方法，但在我们的数据中它的敏感性只有80％，而mixcr和tracer分别为98％和99％。敏感性指鉴定出有效tcr序列的细胞比例。mixcr是一个面向群体细胞的tcr序列鉴定软件，在我们的数据中它的假阳性率为20％，而tracer和vdjpuzzle均不超过1％。假阳性的判定方法是在一个t细胞上发现超过两种ɑ链和β链的tcr序列(一个细胞最多分别表达两种ɑ链和β链的tcr基因型)。对三种软件进行交叉验证时，tracer、mixcr和vdjpuzzle三者的交叉验证率分别为93％,96％,96％，没有显著区别。综上，我们使用了tracer检验的结果。tcr序列用于判断t细胞的克隆性。每个t细胞表达1～2种tcrɑ序列和1～2种tcrβ序列。当两个细胞中至少一种tcrɑ同时至少一种tcrβ的序列完全一致时，这两个细胞被认为来自同一克隆。同时，要求一致的tcrɑ和tcrβ的序列是可以翻译成有效蛋白质的。第三，tcrɑ的tmp值至少大于10，tcrβ的tmp值至少大于15。通过tracer软件识别和判定，并与已知数据库进行对比，共发现46个未被报道过的、在肿瘤组织中浸润的t细胞克隆性表达的tcr，每个tcr的基本情况如下表1-46所示：表1表2表3表4表5表6表7表8表9表10表11表12表13表14表15表16表17表18表19表20表21表22表23表24表25表26表27表28表29表30表31表32表33表34表35表36表37表38表39表40表41表42表43表44表45表46实施例3：tcr、mhc和小肽段结合能力的计算，以及新的肿瘤抗原的预测从实施例2中获得的tcr序列中，选择了来自于最大的克隆组群中的7种tcr序列，它们分别来自于前述肝癌患者中的5个，确定出对应患者的mhc类型(参照szoleka1,schubertb2,mohrc2,sturmm1,feldhahnm1,kohlbachero1.optitype:precisionhlatypingfromnext-generationsequencingdata.bioinformatics.2014dec1；30(23):3310-6.doi:10.1093/bioinformatics/btu548.epub2014aug20.的方法)，利用netmhc和每个患者的rna-seq预测得到了126种小肽序列(参照andreattam,nielsenm.gappedsequencealignmentusingartificialneuralnetworks:applicationtothemhcclassisystem.bioinformatics(2016)feb15；32(4):511-7；nielsenm,lundegaardc,worningp,lauemollersl,lamberthk,buuss,brunaks,lundo.reliablepredictionoft-cellepitopesusingneuralnetworkswithnovelsequencerepresentations.proteinsci.,(2003)12:1007-17的方法)，利用rosettadock软件组装了543个由tcr/pmhc复合物组成的对接模型(表47)。表47肝癌患者的tcr、mhc和小肽的数量患者编号新抗原tcrmhctcr/pmhc复合物p0205141456p03224114164p04075713171p0508111444p1116934108合计543对cdrloop区灵活性的观察和分析，把这部分区域分成6个部分，分别独立的进行和对象蛋白的精细对接模拟。mhc，tcr和小肽结合在一起对接模型，分别进行低分辨率以及高分辨率的对接进程。达到最大迭代次数终止计算。上述543个模型是理论数量，经对接计算后能形成有效对接的tcr/pmhc复合物的数量少于543个。分析结果时，我们从rmsd最小的结构中挑选出打分函数值最低的构象，认为是接合自由能最低的构象，同时表示找到了结合能力最强的新抗原。以5个患者之一p1116，针对其编号为ttc35-1116的免疫细胞克隆获得的tcr为例：该tcr是前述表44中编号为p1116_c0080-1的tcr，其α链和β链的全长核酸序列如表44所述，翻译后α链的氨基酸序列为dqqvkqnspslsvqegrisilncdytnsmfdyflwykkypaegptflisissikdknedgrftvflnksakhlslhivpsqpgdsavyfcaakpgssntgklifgqgttlqvkpdiqnpepavyqlkdprsqdstlclftdfdsqinvpktmesgtfitdktvldmkamdsksngaiawsnqtsftcqdifketnatypssdvpcdatlteksfetdmnlnfqnlsvmglrilllkvagfnllmtlrlwss，β链的氨基酸序列为：eagvaqsprykiiekrqsvafwcnpisghatlywyqqilgqgpklliqfqnngvvddsqlpkdrfsaerlkgvdstlkiqpakledsavylcasslggggvqyfgpgtrllvledlrnvtppkvslfepskaeiankqkatlvclargffpdhvelswwvngkevhsgvstdpqaykesnysyclssrlrvsatfwhnp。利用rosettadock软件，得到该tcr的7个loop区的分区信息：fr1-imgt：残基位置1至26号残基，cdr1-imgt：27至38位残基，fr2-imgt：39至55号残基，cdr2-imgt：56至65号残基，fr3-imgt：第66至104号残基，cdr3-imgt：第105位至117号残基，fr4-imgt:118-128号残基。我们选取离小肽主链原子中心最近的6个loop区。基于该loop区信息，利用rosettadock软件构建mhc，tcr和9个小肽对接模型，得到了以rosettascore为打分函数值的结果(表48)。从结果可知，肽段“klahvaaqa”是此患者给定的tcr结合最强的小肽。该小肽可用于后续肿瘤新抗原的进一步确认研究，以及相应抗体和免疫细胞的开发等。而该tcr则可用于后续tcr-t疗法的实证研究。表48患者p1116的ttc35-1116细胞的hla类型，肽段序列和rosettascore其他的6个tcr分析结果如下：tcridhlascore小肽全长p0205_c0010a*02:061522.46tkpevdftvp0205_c0010a*02:061548.7krmcivwgvp0205_c0010a*02:061689.35lvdcevvslp0205_c0010a*02:061720.41gvlvgwnllp0205_c0010a*02:061725.62wsyhfiglip0205_c0010a*02:061767.14nllclqgrvp0205_c0010a*24:021904.04syhfigliwp0205_c0010a*02:061954.86rymwsyhfip0205_c0010a*02:062289.03cqparsdslp0205_c0010a*02:062774.4itkpevdftvp0205_c0010a*02:062826.26ftvevfrsalp0205_c0010a*02:063691.6ilvdcevvslp0205_c0010a*02:064017.51rklglcefpvp0205_c0010c*08:014022.75ftvevfrsalp0205_c0010a*02:064031.15vlvgwnllclp0205_c0010a*02:064150.2ymwsyhfiglp0205_c0010a*02:064255.86aqareptaslp0205_c0010a*24:024423.03rymwsyhfigp0205_c0010a*02:065120.58kipgellrwltcridhlascorepep_fullnamep1116_c0161a*02:034683.35qlnpdvpdap1116_c0161a*02:034695.84avmsllrrvp1116_c0161a*02:034729.57klahvaaqatcridhlascorepep_fullnamep1116_c0014a*02:039375.62avmsllrrvp1116_c0014a*02:039123.01klahvaaqap1116_c0014a*02:038976.5qlnpdvpdatcridhlascorepep_fullnamep0508_c0024-2a*02:018038.82qldyvdvylp0508_c0024-2c*06:028700.85erfgrdlnyp0508_c0024-2a*02:018764.63flhswetptp0508_c0024-2a*30:018923.12carraqneap0508_c0024-2b*57:019058.1rlqdqdpswp0508_c0024-2b*57:019389.24raimempffp0508_c0024-2a*02:019988.57iiflhswetp0508_c0024-2a*02:0110545klqldyvdvp0508_c0024-2b*57:0112149.51leiiflhswp0508_c0024-2a*30:0112602.61rslleqdvkp0508_c0024-2a*02:0112767.34allqkyagi对比实施例：1、cdna建库方法的对比将smart-seq2方法和汤富酬创建的方法(为方便描述，以下简称为“tang2009”，tang,f.etal.rna-seqanalysistocapturethetranscriptomelandscapeofasinglecell.nat.protoc.5,516–535(2010).)在cdna扩增效率上：在完成cdna扩增后，以无核酸的超纯水为阴性对照，利用rt-pcr来检测管家基因β-actin、gapdh或cd3基因的表达来检验扩增效率。判断基因是否表达的标准有两条：一是细胞样品rt-pcr的ct值明显小于阴性对照；二是细胞样品的溶解曲线峰值在85℃到90℃之间(阴性对照约为78℃)。根据实验的结果，用smart-seq2法扩增后，在绝大部分细胞中均可以检测到β-actin、gapdh和cd3。用tang2009法扩增后可在大部分细胞中检测到β-actin和gapdh的表达，仅能在极少数细胞中检测到cd3。在cdna文库质量上：文库的构建不受扩增方法影响，用tang2009法和smart-seq2法扩增均能提供构建文库所需的cdna量，所建文库的片段大小均符合illumina测序仪要求。但如果对比最终构建所得的cdna文库数和用于建库的起始单细胞数的比值，以三个肝癌患者的样本比较来看，用smart-seq2法扩增t细胞的成功率总体上高于tang2009法，如下表所示。构建cdna文库成功率表患者编号建库方法ptcpthptrttctthttrntcnthntr20141202tang200946％24％26％62％52％78％///20150205smart-seq280％80％67％80％60％67％///20150508smart-seq292％74％92％95％92％88％70％86％/注：“/”表示无该类型样本。ptc为外周血中的细胞毒t细胞，pth为外周血中的辅助性t细胞，ptr为外周血中的调节性t细胞。ttc为癌组织中的细胞毒t细胞，tth为癌组织中的辅助性t细胞，ttr为癌组织中的调节性t细胞。ntc为正常肝脏组织中的细胞毒t细胞，nth为正常肝脏组织中的辅助性t细胞，ntr为正常肝脏组织中的调节性t细胞。以下每个表的简写相同。在测序质量上：以三个肝癌患者的样本为例，具体参数见下表。tang2009法(患者20141202)样本平均测序质量表smart-seq2法(患者20150205)样本平均测序质量表细胞类型总读段高质量读段率比对率表达基因数ptc156141997.88％62.47％3042pth161040398.33％79.72％2553ptr148232398.42％60.33％3103ttc144257399.51％51.83％3071tth189006977.07％66.7％2478ttr155868598.49％67.87％3478smart-seq2法(患者20150508)样本平均测序质量表细胞类型总读段比对率高质量读段率表达基因数ntc181983097.10％99.25％2187nth143165693.38％99.44％2267ptc1476204.595.40％99.38％2745pth153059086.35％99.39％2026ptr151319096.09％99.43％2295.5ttc160217595.83％99.38％2823tth1562359.594.52％99.38％2795.5ttr151448696.33％99.43％2907作为单细胞rna-seq的数据，多数细胞的高质量读段率、比对率均较高，说明测序质量好，所得数据适合做表达量和表达序列的分析。从检测到的基因数方面看，smart-seq2法所能检测到的基因数明显高于tang2009法。在扩增偏向上，tang2009法和smart-seq2法扩增cdna时存在不同的倾向性。从测序结果的分析图中不难看出，tang2009法更倾向于扩增出cdna的3’端，这种偏向性不利于整个tcr结构的组装，而smart-seq2法则能更为均一地扩增cdna，有利于组装出tcr的全部序列(附图5和6)。参考文献：1.eltahlaaa,rizzettos,pirozyanmr,betz-stableinbd,venturiv,kedzierskak,lloydar,bullra,lucianif.linkingthetcellreceptortothesinglecelltranscriptomeinantigen-specifichumantcells.immunolcellbiol.2016jul；94(6):604-11.doi:10.1038/icb.2016.16.epub2016feb10.2.bolotinda,poslavskys,mitrophanovi,shugaym,mamedoviz,putintsevaev,chudakovdm.mixcr:softwareforcomprehensiveadaptiveimmunityprofiling.natmethods.2015may；12(5):380-1.doi:10.1038/nmeth.3364.3.stubbingtonmj,t,proserpiov,clares,speakao,dougang,teichmannsa.tcellfateandclonalityinferencefromsingle-celltranscriptomes.natmethods.2016apr；13(4):329-32.doi:10.1038/nmeth.3800.epub2016mar7.4.grabherr,m.g.,haas,b.j.,yassour,m.,levin,j.z.,thompson,d.a.,amit,i.,adiconis,x.,fan,l.,raychowdhury,r.,zeng,q.,etal.(2011).full-lengthtranscriptomeassemblyfromrna-seqdatawithoutareferencegenome.nat.biotechnol.29,644–652.5.langmead,b.,andsalzberg,s.l.(2012).fastgapped-readalignmentwithbowtie2.nat.methods9,357–359.6.lefranc,m.-p.,giudicelli,v.,ginestoux,c.,jabado-michaloud,j.,folch,g.,bellahcene,f.,wu,y.,gemrot,e.,brochet,x.,lane,j.,etal.(2009).imgt,theinternationalimmunogeneticsinformationsystem.nucleicacidsres.37,d1006–d1012.当前第1页12当前第1页12