CN115266962B - 蛋白标志物在制备弥漫型胃癌分子分型的产品中的应用及弥漫型胃癌分子分型的分类器 - Google Patents

蛋白标志物在制备弥漫型胃癌分子分型的产品中的应用及弥漫型胃癌分子分型的分类器 Download PDF

Info

Publication number
CN115266962B
CN115266962B CN202210562342.7A CN202210562342A CN115266962B CN 115266962 B CN115266962 B CN 115266962B CN 202210562342 A CN202210562342 A CN 202210562342A CN 115266962 B CN115266962 B CN 115266962B
Authority
CN
China
Prior art keywords
proteins
gastric cancer
cancer
protein
diffuse
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210562342.7A
Other languages
English (en)
Other versions
CN115266962A (zh
Inventor
秦钧
丁琛
刘明伟
葛赛
夏夏
周泉
冯晋文
甄蓓
汪宜
沈琳
贺福初
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Guhai Tianmu Biomedical Technology Co ltd
Original Assignee
Beijing Guhai Tianmu Biomedical Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Guhai Tianmu Biomedical Technology Co ltd filed Critical Beijing Guhai Tianmu Biomedical Technology Co ltd
Priority to CN202210562342.7A priority Critical patent/CN115266962B/zh
Publication of CN115266962A publication Critical patent/CN115266962A/zh
Application granted granted Critical
Publication of CN115266962B publication Critical patent/CN115266962B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • G01N30/04Preparation or injection of sample to be analysed
    • G01N30/06Preparation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • G01N30/62Detectors specially adapted therefor
    • G01N30/72Mass spectrometers

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明提供蛋白标志物在用于制备弥漫型胃癌分子分型的产品中的应用及弥漫型胃癌分子分型的分类器,提供的分类器中包含蛋白标志物的表达量数值,可将Lauren分型中的弥漫型胃癌进一步分为三种亚型,分别为细胞周期型弥漫型胃癌PX1、上皮‑间充质转化型弥漫型胃癌PX2和免疫富集型弥漫型胃癌PX3。本发明可以为弥漫型胃癌的分子分型、预后以及设计靶向药物奠定基础。

Description

蛋白标志物在制备弥漫型胃癌分子分型的产品中的应用及弥 漫型胃癌分子分型的分类器
本申请为申请日2018年02月13日、申请号201810149737.8、发明名称“弥漫型胃癌的分子分型、用于分型的蛋白标志物及其筛选方法和应用”的分案申请,并享有其优先权。
技术领域
本发明属于医药生物领域,涉及新的肿瘤分子分型以及分型中生物标志物的筛选技术,特别涉及利用蛋白质组定量技术对胃癌患者的胃癌组织和配对的癌旁组织的蛋白质组进行定量检测,通过数据分析筛选获得Lauren分型中弥漫型胃癌新的亚型及筛选分型相关的蛋白标志物的方法。
背景技术
胃癌在全球肿瘤死亡原因中排名第三,我国是胃癌高发地区之一,每年新发病例接近全世界的一半左右。目前应用较广的胃癌分型有基于形态学的Lauren分型。Lauren分型将胃癌分为仍保留腺腔结构的肠型,无腺腔结构、肿瘤细胞弥漫散在的弥漫型,以及介于两者之间的混合型。弥漫型胃癌约占胃癌总数的30%,预后最差。此外,应用较广的还有TNM分期(tumor-node-metastasis)。TNM分期按照肿瘤的大小、侵袭能力、淋巴结节个数以及是否发生近端或者远端转移,对胃癌病人进行分期。Lauren分型的预后能力有限,而TNM分期不能发现肿瘤的分子特征,对治疗的指导意义有限。TCGA研究机构曾经根据胃癌基因组数据和转录组数据对胃癌进行分子分型,可分为4个亚型,分别为:EBV阳性、微卫星不稳定(MSI)、基因组稳定和染色体不稳定型。其中75%的基因组稳定型是弥漫型胃癌。ACRG研究机构将胃癌分为4个亚型,分别为:MSI型、微卫星稳定/上皮细胞向间充质细胞转化型(MSS/EMT)、微卫星稳定/TP53阳性型和微卫星稳定/TP53阴性型。MSS/EMT型主要为弥漫型胃癌。TCGA分型和ACRG分型都是基于基因组和转录组数据进行的分型。而作为生命活动的最终执行者,蛋白质组数据进行分子分型并提出潜在有效的药物靶点具有重大的临床价值。目前尚未有采用蛋白质组学技术对胃癌进行分子分型的技术,针对胃癌分型中的弥漫型胃癌尚无进一步研究其亚型的手段。
发明内容
本发明对弥漫型胃癌进行了进一步分型,得到三组临床特征和分子特征均有明显差异的亚型,并提供了蛋白生物标志物在弥漫型胃癌分子分型中的应用,还构建了用于弥漫型胃癌分子分型的分类器。
本发明的一个方面提供了蛋白标志物在制备弥漫型胃癌进行分子分型的分类器中的应用,其中,所述弥漫型胃癌的分子分型包括弥漫型胃癌PX1亚型、弥漫型胃癌PX2亚型及弥漫型胃癌PX3亚型,
所述弥漫型胃癌PX1亚型的蛋白标志物选自如下任意多种蛋白:RPL23、RPL21、RPS15、RPL38、RPL27、RPS27A、CAMP、RPL36、RPL27A、AZU1、RPL7A、RPL30、MPO、RPL9、PRTN3、RPS20、SBDS、ELANE、MMP9、RPL31、FKBP3、RPL24、RPL10、S100A8、SNRPD2、RPL37A、RPS18、NDUFA6、LTF、RPL22、NDUFS7、BCCIP、S100A9、ATP5C1、HINT2、RBBP9、SRP14、ATP5O、C4BPA、C1R、RPL8、ISOC2;
所述弥漫型胃癌PX2亚型的蛋白标志物选自如下任意多种蛋白:SLC4A2、ATP4B、ATP4A、DBT、CLIC6、VIL1、CKM、GIF、AUH、NNT、IDH2、GLUL、ECI2、BCKDHA、PDIA2、FABP1、C14orf159、ECHDC3、CKMT2、ACSS1、BCKDHB、NDUFB7、CYB5R1、MPC2、MFSD4、NDUFB11、REG3A、APOO、SLC25A4、MLYCD、TXNRD2、ACAT1、SMIM20、TMEM126A、ACADSB、SERPINB5、PCCB、C14orf2、GOT1、CA9、ALDH4A1、CHCHD4、APOOL、FABP3、ADCK3、ATP5F1、ME3、DCXR、Marc2、TPD52L1、NDUFB5、AQP4、CASP1、COA3、FAM162A、LDHD、COX1、NDUFA4、HIP1R、COX4I1、CRAT、PDHX、COQ10B、ATP5L、ADA、ETFDH、BCAT2、MTFP1、NDUFB6、ACADS、TGFBI、GBAS、ND4、FAHD2A、NDUFS7、TMEM214、CTH、COQ9、NDUFA11、ND1、NDUFA3、ALDH6A1、COX17、DLAT、NDUFA9、NUDT19、PPIF、SDHD、PCCA、ALDH5A1、ND5、FHL2、COX6C、ATP5D、PDHB、ALDH1L2、ATPAF2、MCCC2、BDH1、NDUFS2、CLCA1、HINT2、ACO2、FBP2、NDUFS4、OGDH、TPPP3、MTX2、ATP5A1、FGD4、CARKD、PDLIM1、RAP1GAP、COX3、SUOX、REG4、ATP5B、MCCC1、COX7A2L、COX2、NDUFC2、OPLAH、NDUFA7、L2HGDH、AIFM1、PDLIM3、VCAN、OXCT1、PC、CYTB、PYCARD、NADK2、NDUFAB1、C21orf33、SLC25A20、ANPEP、ECHDC2;
所述弥漫型胃癌PX3亚型的蛋白标志物选自如下多种蛋白:ATP4A、ATP4B、CLIC6、SLC4A2、GIF、GLUL、CKM、MPO、PRTN3、NNT、AZU1、IDH2、MLYCD、DBT、MMP9、C14orf2、CAMP、C14orf159、PDIA2、NDUFB11、ATP5F1、CKMT2、AUH、MPC2、ALDH4A1、ATP5L、NDUFB5、FABP3、MFSD4、ND4、COX4I1、CYB5R1、ND1、SLC25A4、BCKDHA、SDHD、MMP8、CA9、NDUFB6、S100A12、APOO、BCKDHB、COA3、TXNRD2、TPD52L1、BPI、ME3、HIP1R、ELANE、MTFP1、COX1、ITGAM、PGLYRP1、NDUFA11、TMEM126A、AQP4、VBP1、RPL27A、EPCAM、REG4、SF3A2、RPL22、GSDMB、RBP2、RPL30、SRP14、SNRPE、EIF2AK2、RPL9、RPS27A、SBDS、ACAT2、ME1、TPPP3、DCTPP1、RPL27、UGDH、FKBP3、CNN1、LSM4、OTC、RBM39、REG3A、RBM14、RPL23、HSD17B2、SERPINB5、FABP1、VIL1;
其中PX1为细胞周期型弥漫型胃癌,PX2为上皮-间充质转化型弥漫型胃癌,PX3为免疫富集型弥漫型胃癌。
在一些实施方式中,所述弥漫型胃癌PX1亚型、所述弥漫型胃癌PX2亚型及所述弥漫型胃癌PX3亚型存在以下临床特征差异或分子特征差异:
所述临床特征方面差异包括:(1)三个亚型的肿瘤原发部位差异显著(Fisher’s精确检验P=0.000),一半以上的PX3样本原发部位在胃窦区;(2)三个亚型的癌旁组织性质差异显著(Fisher’s精确检验P=0.001),3/4以上的PX2样本的癌旁只是有浅表性胃炎,不是慢性萎缩性胃炎或肠化生;(3)三个亚型的总生存时间差异显著(Kaplan–Meier分析P=0.038);
所述分子特征方面的差异为:
(A)PX1和PX2中癌和癌旁差异表达的蛋白富集在细胞周期和DNA复制相关通路,而PX3中癌和癌旁差异表达的蛋白未富集在细胞周期或DNA复制相关通路中;
(B)PX1-PX3中癌和癌旁差异表达的蛋白均富集在E2F靶基因和MYC靶基因的基因集中,但PX1在这些通路中的富集程度更高;
(C)PX1中癌和癌旁差异表达的蛋白富集在RNA处理和DNA修复相关通路,而PX2和PX3中癌和癌旁差异表达的蛋白未富集在这些通路中;
(D)PX2和PX3中癌和癌旁差异表达的蛋白富集在上皮-间充质转化过程、补体、细胞外基质和干扰素gamma响应的相关通路中,而PX1中癌和癌旁差异表达的蛋白未富集在这些通路中;
(E)PX3中癌和癌旁差异表达的蛋白富集在NFKB介导的TNF-alpha信号通路、免疫响应和细胞粘附的相关通路中,而PX1和PX2中癌和癌旁差异表达的蛋白未富集在这些通路中。
本发明再一方面还提供一种弥漫型胃癌分子分型的分类器的构建方法,其包括:
1)采集数量为1/2S对来自弥漫型胃癌患者的癌与癌旁组织样本的蛋白质组表达谱数据集,包括全部样本S归集的蛋白质的数量、种类以及各种蛋白的表达量FOT值;S为包括胃癌组织样本和胃癌旁组织样本全部样本的数量,S/2为患者人数,即样本人数;
2)对全部样本S对应的各种蛋白的表达量FOT值进行筛查,要求定量的蛋白含有至少两条专属肽段(unique peptide),并且离子分值(ion score)大于等于20,获得初筛数量Q1的蛋白在全部样本中的FOT值;
3)从Q1数量的蛋白中选择在至少1/6×S个样本(占总样本数的六分之一)中FOT值大于10-5的数量为Q2的蛋白,并将FOT值小于10-5的值用10-5替换;进一步从Q2数量的蛋白中筛选在至少1/10×S/2个样本人数(占样本人数的十分之一)中T/N大于3倍,或者小于1/3倍的蛋白,共得到数量为Q3个蛋白;
4)获取上述数量为Q3蛋白FOT的T/N的值,并对T/N的值取以10为底的对数变换;
5)采用一致性聚类(consensus clustering)对上述筛选得到的Q3个蛋白在S样本中的log10(T/N)进行分子分型,得到蛋白质组学技术对弥漫型胃癌的3个亚型;
6)根据差异表达的FDR p value和高表达3倍以上的频率与低表达1/3倍以上的频率的差值所占的比例两个指标来定义差异表达蛋白,从初筛数量Q1的蛋白中获取数量为Q4弥漫型胃癌亚型的差异表达蛋白;
7)选取上一步骤中获得的数量Q4的差异表达蛋白和用于一致性聚类的数量Q3的蛋白的交集数量为Q5的蛋白,采用PAM进行分类器的构建,选取分型平均错误率小于某一定值的数量为Q6蛋白作为胃癌分子分型蛋白标志物,建立弥漫型胃癌分子分型的分类器。
在一些实施方式中,步骤7)中所述PAM采用重心邻近法(nearest centroid),蛋白在每个亚型中的重心值=蛋白在每一个亚型中的平均值/蛋白的组内方差。
在一些实施方式中,步骤1)中弥漫型胃癌患者的癌与癌旁组织样本数量为84对,即1/2S样本人数为84,则全部样本S数量为168。
在一些实施方式中,步骤2)中初筛蛋白数量Q1为9186个。
在一些实施方式中,步骤3)中蛋白数量Q2为3619个,1/6×S个样本为28个,1/10×S/2个样本人数为8个;所述蛋白数量Q3为2538个。
在一些实施方式中,步骤5)中确定了弥漫型胃癌的3个亚型,分别命名为PX1、PX2和PX3。
在一些实施方式中,步骤6)中差异表达蛋白数量Q4为2264个。
在一些实施方式中,步骤7)中交集蛋白数量Q5为875个,弥漫型胃癌的3个亚型的分型平均错误率为9.6%,确定的弥漫型胃癌分子分型蛋白标志物数量Q6为199个。
在一些实施方式中,所述3个亚型PX1、PX2和PX3存在以下临床特征差异或分子特征差异:
所述临床特征方面差异包括:(1)三个亚型的肿瘤原发部位差异显著(Fisher’s精确检验P=0.000),一半以上的PX3样本原发部位在胃窦区;(2)三个亚型的癌旁组织性质差异显著(Fisher’s精确检验P=0.001),3/4以上的PX2样本的癌旁只是有浅表性胃炎,不是慢性萎缩性胃炎或肠化生;(3)三个亚型的总生存时间差异显著(Kaplan–Meier分析P=0.038);
所述分子特征方面的差异为:
(A)PX1和PX2中癌和癌旁差异表达的蛋白富集在细胞周期和DNA复制相关通路,而PX3中癌和癌旁差异表达的蛋白未富集在细胞周期或DNA复制相关通路中;
(B)PX1-PX3中癌和癌旁差异表达的蛋白均富集在E2F靶基因和MYC靶基因的基因集中,但PX1在这些通路中的富集程度更高;
(C)PX1中癌和癌旁差异表达的蛋白富集在RNA处理和DNA修复相关通路,而PX2和PX3中癌和癌旁差异表达的蛋白未富集在这些通路中;
(D)PX2和PX3中癌和癌旁差异表达的蛋白富集在上皮-间充质转化过程、补体、细胞外基质和干扰素gamma响应的相关通路中,而PX1中癌和癌旁差异表达的蛋白未富集在这些通路中;
(E)PX3中癌和癌旁差异表达的蛋白富集在NFKB介导的TNF-alpha信号通路、免疫响应和细胞粘附的相关通路中,而PX1和PX2中癌和癌旁差异表达的蛋白未富集在这些通路中。
在一些实施方式中,步骤1)中所述采集数量为1/2S对来自弥漫型胃癌患者的癌与癌旁组织样本的蛋白质组表达谱数据集的操作包括以下步骤:
11)采集数量为1/2S对来自弥漫型胃癌患者的癌与癌旁组织样本,作为1/2S对胃癌样本;
12)胃癌蛋白样品制备;该步骤是将收集的组织样本裂解为全蛋白裂解液,再还原烷基化后用胰酶消化得到的肽即为胃癌蛋白样品;其中所述组织样本为癌组织样本和癌旁组织样本,两者分别制备并且每一患者的两种样本配对;
13)胃癌蛋白样品的质谱检测;
14)胃癌蛋白样品的质谱数据分析;
15)将所有样本对应的胃癌蛋白样品质谱数据归集,得到胃癌患者癌与癌旁的蛋白质组表达谱数据集。
在一些实施方式中,步骤12)中所述裂解的具体操作为:将组织样本置于8M Urea蛋白裂解液和PMSF中,碾碎,随后进行超声裂解,之后置于冰上裂解,期间不时弹动,获得裂解液,其中所述超声裂解的条件为3秒开始,3秒间歇,振幅25%;将裂解液离心后取上清获得全蛋白裂解液;
所述胰酶消化的具体操作为:向完成还原烷基化的裂解液中加入5倍体积的pH8.2Tris-Hcl,终浓度1mM的CaCl2使溶液呈弱碱性;按照胰酶:蛋白的质量比为1:50确定用于消化蛋白的胰酶量,在裂解液中加入3/4胰酶的量,全部转入离心管中,37℃旋转消化过夜;次日添加剩余1/4胰酶,消化4小时,消化完成后加入甲酸终止消化,并离心取上清,获得胃癌蛋白样品。
在一些实施方式中,步骤11)中所述胃癌样本收集为选取术前未接受新辅助放化疗的胃癌患者的癌组织与癌旁组织,要求癌组织样本中肿瘤细胞含量超过50%,癌旁组织样本中不含肿瘤细胞;步骤13)中所述胃癌蛋白样品的质谱检测包括对胃癌蛋白样品的肽进行sRP液相分离为多个组分的肽段,对每个组分的肽段用液相色谱串联质谱系统进行检测,合并多个组分的检测结果得到该胃癌蛋白样品对应的全部蛋白质谱数据;步骤14)中针对步骤13)取得的胃癌蛋白样品的全部蛋白质谱数据在蛋白数据库中搜库、定量确定对应该胃癌蛋白样品中的蛋白种类、数量和各种蛋白的表达量FOT值。
本发明还包括基于84对弥漫型胃癌样本癌与癌旁的蛋白质表达谱数据集筛选用于弥漫性胃癌分子分型的蛋白标志物的方法,筛选得到的199个蛋白标志物,以及蛋白标志物用于分型的方法。
一、84例弥漫型胃癌样本癌与癌旁的组织蛋白质表达谱数据集
1、收集弥漫型胃癌患者信息;采集符合条件的84例弥漫型胃癌患者的癌与癌旁组织;制备肿瘤组织的癌与癌旁的蛋白样品。
2、对癌与癌旁蛋白样品进行蛋白质组质谱测定、数据库搜索、肽段定量及蛋白拼接组装,确定每一个组织蛋白样品中的蛋白种类及各蛋白的定量,共鉴定到11340个基因产物。
二、基于蛋白质表达谱数据集筛选用于弥漫型胃癌分子分型的蛋白标志物
基于蛋白表达谱数据对胃癌进行分型的流程图如图1所示。
1、基于肽段的专属性(unique peptide)和离子分值(ion score)评分筛选蛋白,获得了9186个蛋白在84例样本的癌与癌旁组织中的表达值。
2、基于蛋白的表达水平和在样本中出现的频率和差异表达的频率,筛选得到表达水平较高(FOT大于10-5),出现频率至少28(168/6)个样本,并且在至少8个样本中T/N大于3倍,或者小于1/3倍的2538个蛋白。
3、根据上述筛选得到的2538个蛋白,采用一致性聚类(consensus clustering)对84对样本的T/N的蛋白质组定量数据进行分子分型。最终得到了3个亚型(图2,图3,图4)。三个亚型之间的总生存时间存在差异(图5)。
4、采用SAM软件计算三个亚型中每个亚型T与N的差异表达蛋白,初始数据集采用包含9186个蛋白,84对样本的数据集。我们选取SAM FDR q值<0.01,并且(高表达频率-低表达频率)/总检出频率>0.5的蛋白为显著高表达蛋白,而(高表达频率-低表达频率)/总检出频率<-0.5的蛋白为显著低表达蛋白。第一个亚型中获得了580个显著高表达的蛋白和8个显著低表达的蛋白;第二个亚型中获得了1016个显著高表达的蛋白和572个显著低表达的蛋白;第三个亚型中获得了457个显著高表达的蛋白和551个显著低表达的蛋白。共有2264个差异表达蛋白。
5、采用PAM进行分类器的构建。初始蛋白集合采用上一步骤获得的2264个差异表达蛋白和2538个蛋白的交集共875个蛋白。将84对样本随机分为两个数据集。第一个数据集中包含随机抽取的62对样本(75%);剩下的22对样本(25%)作为第二个数据集。采用PAM对第一个数据集中的样本进行分类,考虑到平均错误率和最终获得的蛋白集合的平衡,我们选取了平均错误率为9.6%的199个蛋白作为分型蛋白。将这199个蛋白的集合用于预测第二个数据集,错误率为4.5%(图7)。
三、应用199个蛋白标志物进行弥漫型胃癌的分子分型
应用199个蛋白标志物进行弥漫型胃癌的分子分型的流程图如图8所示。
1、采集成对的弥漫型胃癌的癌与癌旁的蛋白表达谱数据。
2、检测样本的蛋白质组表达谱。
3、要求定量的蛋白含有至少两条专属肽段(unique peptide),并且离子分值(ionscore)大于等于20。
4、将FOT值小于10-5的值用10-5替换。获取199个蛋白中每个蛋白FOT的T/N的值,并进行以10为底的log变换。
5、将之前用PAM分析62例样本数据获得的模型,即199个蛋白分别在三个弥漫型胃癌亚型中的重心值(表5)用于预测新样本的分型:计算新样本与三个亚型重心的秩相关系数(Spearman’s rank correlation),将这例样本归入秩相关系数最高的亚型。
本发明提供弥漫型胃癌进一步分型获得的三个亚型,提供分子分型的蛋白标志物及其筛选方法和应用,包括84例弥漫型胃癌样本癌与癌旁的组织蛋白表达谱数据,基于蛋白表达谱数据对胃癌进行分型的方法以及用于胃癌分型的199个蛋白标志物。采用这种方法可以对Lauren分型中的弥漫型胃癌进行进一步的分类。可分为三类:一类为细胞周期型(cell cycle subtype)弥漫型胃癌,一类为上皮-间充质转化(Epithelial-mesenchymaltransition,EMT)型弥漫型胃癌,另一类为免疫富集型弥漫型胃癌。本发明为弥漫型胃癌的分类以及设计靶向药物奠定了基础。
附图说明
图1为流程图,显示本发明蛋白标志物的筛选流程。
图2为一致性聚类图,显示84对样本的T/N的蛋白质组定量数据进行分子分型的一致性聚类图;类别1对应浅蓝色,类别2对应深蓝色,类别3对应绿色。
图3为累积分布函数曲线图(CDF),显示类别数为2至6时的累积分布函数曲线图,线条从左向右分别为类别数k=6、k=5、k=4、k=3和k=2。
图4为累计分布函数曲线下面积的折线图,显示图2的累积分布函数曲线(CDF)从k-1到k类增加的曲线下面积。
图5为三个亚型的生存曲线(Kaplan-Meier analysis)图;线条由上向下分别对应PX1、PX2和PX3。
图6为热图,具体为三个亚型中癌与癌旁差异表达的蛋白富集通路的热图。
图7为分型错误率曲线图,显示采用PAM筛选蛋白标志物时采用的阈值及此阈值下的分型错误率。
图8为流程图,显示采用199个蛋白标志物进行分子分型的流程。
具体实施方式
本发明旨在提供一种获取胃癌、具体为弥漫型胃癌分子分型的方法,并进一步提出用于弥漫型胃癌分型的蛋白标志物。通过检测胃癌与癌旁的蛋白表达谱,可以发现癌与癌旁的差异表达蛋白,并进一步发现新的胃癌分子亚型,筛选新的蛋白标志物。
本发明具体通过收集弥漫型胃癌样本,对癌与癌旁组织进行蛋白质组定量,基于蛋白质组表达谱数据对弥漫型样本进行分子分型,筛选弥漫型胃癌亚型中的差异表达蛋白以及筛选用于分型的蛋白标志物。蛋白质组的概念是指细胞内、组织内、体液内或个体内全部种类蛋白的集合。本发明中,蛋白质组特指胃癌或癌旁组织样本中所包括的全部不同种类的蛋白。
为达成以上成果,本发明结合图1所示蛋白标志物的筛选流程和图8所示采用所确定的蛋白标志物进行分子分型的流程就以下几方面内容做出说明:
一、胃癌样本收集
胃癌样本收集是指选取术前未接受新辅助放化疗的胃癌患者的癌组织与癌旁组织,要求癌组织样本中肿瘤细胞含量超过50%,癌旁组织样本中不含肿瘤细胞。以下述描述为实例,其中弥漫型胃癌及其数据仅为示例,不作为限制:
2012年12月至2015年7月于北京肿瘤医院进行全胃切除术或胃大部切除术的胃癌患者共2451例,根据患者资料进行了初筛,排除了术前接受新辅助放化疗患者428例、病理诊断为胃间质瘤患者92例、Lauren分型为肠型胃癌患者751例及混合型胃癌患者644例、保存的新鲜肿瘤组织或癌旁组织太小无法开展研究者390例,2例患者既往具有恶性肿瘤病史。剩余144例弥漫型胃癌进行后续深入筛选。分析上述144例弥漫型胃癌患者的新鲜肿瘤组织及配对癌旁正常组织的苏木精-伊红(Hematoxylin and eosin,H&E)染色结果后,排除了肿瘤组织中肿瘤细胞含量不足50%者56例、癌旁正常组织中含有少量肿瘤细胞者3例、癌旁正常组织中含有大量肌纤维者1例,得到来自84例患者的癌组织与癌旁组织样本84对。
二、胃癌蛋白样品的制备
本发明中,“样本”特指收集的组织样本,用于检测的称为“样品”。
1、称取空EP管重量并做好标记。每一对癌组织(Tumor,T)及配对的癌旁正常组织(Nearby tissue,N)样本保证同时平行处理。组织样本由-80℃取出后,首先在EP管内称重并核对、记录,埋入冰盒。在EP管中加入约5倍体积8M Urea蛋白裂解液,同时加入PMSF,使用一次性研磨棒碾碎,随后进行一分钟超声裂解(3秒开始,3秒间歇,振幅25%),之后置于冰上裂解15分钟,期间不时弹动。
2、裂解液于14000g,15℃条件下离心10分钟,吸取上清至新EP管,记录上清体积,即为全蛋白裂解液。
3、考马斯亮蓝(Bradford)法蛋白浓度测定,记录蛋白浓度。每例样本取相当于100μg全蛋白的裂解液,余下部分存于-80℃冰箱,登记入库。
4、在100μg全蛋白裂解液中,加入10mM DTT,在56℃摇床震荡30分钟,冷却至室温后加入10mM IAA,置于避光条件下30分钟,完成还原烷基化过程。
5、在完成还原烷基化的裂解液中加入5倍体积的pH8.2 Tris-Hcl,终浓度1mM的Cacl2,使溶液呈弱碱性。按照胰酶:蛋白为1:50的比例(质量比),确定用于消化蛋白的胰酶量,在裂解液中加入3/4胰酶的量,全部转入0.5ml离心管中,37℃旋转消化过夜。次日添加剩余1/4胰酶,消化4小时。消化完成后加入1/100体积甲酸终止消化,并离心取上清,于真空泵中抽干,得到消化后肽样品。每一个组织样本对应一个肽样品即胃癌蛋白样品。
三、胃癌蛋白样品的质谱检测
1、sRP液相分离
肽样品首次分离,使用特别制作的用Tip头填充的反向C18分离柱。该反向C18分离柱是在Tip头底部制作C18筛板,将C18填料悬浮于乙腈并填充于Tip头中,制作为分离柱。使用50%乙腈对该分离柱进行两次冲洗,pH=10的缓冲液(10mM碳酸氢铵)平衡分离柱,同时使用pH=10的缓冲液对步骤二中抽干的肽样品复溶,复溶后的肽样品重复四次载入分离柱,之后采用逐渐升高的梯度比例乙腈溶液进行肽样品洗脱,洗脱乙腈浓度依次为6%、9%、12%、15%、18%、21%、25%、30%和35%(体积浓度)。得到的9个浓度洗脱溶液合并(6%和25%合并,9%和30%合并,12%和35%合并)分为6组分的肽段,真空抽干后用于质谱检测。
2.肽样品质谱检测
每个组分的肽段使用Orbitrap Fusion串联Easy-nLC 1000液相色谱串联质谱系统进行检测,合并6个组分的检测结果得到该肽样品对应的全部蛋白的质谱数据。具体操作为:
抽干的肽段使用A液(0.1%甲酸)复溶,上样至trap柱(100μm×2cm,粒子大小,3μm;孔径大小,),后在分析柱上进行分离(150μm×12cm,粒子大小,1.9μm;孔径大小,),梯度为3%-30%的B液(0.1%甲酸溶于乙腈)洗脱,流速600nl/min,合计洗脱75分钟。
质谱检测参数设置如下,母离子使用Orbitrap进行扫描,扫描范围为质荷比(m/z)300-1400,使用120000(200m/z)的分辨率。母离子使用top-speed模式,在Quadrupole进行选择后(1.6m/z分离窗隔离窗口)的母离子使用HCD(higher energy collision induceddissociation,高能碰撞诱导解离)进行碎裂,碎裂能量35%。子离子碎片在线性离子阱中进行检测,质谱参数设置如下:在一级全扫描时,离子阱中允许进入的最大离子数(automatic gain control target,AGC target)为5e5 ions,最大注入时间50ms,MS/MS二级扫描AGC target为5e3,最大注入时间35ms,动态排除设定18秒。液相色谱串联质谱系统使用Xcalibur软件控制进行数据采集。
四、胃癌蛋白样品的质谱数据分析
质谱鉴定得到的图谱,使用搭载Mascot 2.3的Proteome Discover 1.4搜库,数据库使用人源NCBI(human National Center for Biotechnology Information)蛋白数据库(04-07-2013版本,共32015个蛋白种类)。
搜库参数设置如下:
Fusion:母离子质量偏差范围质谱低限(mass tolerances)为20个ppm,子离子为0.5Da。
最大允许两个漏切位点,固定修饰为carbamidomethyl(C),动态修饰为proteinacetyl(protein N-term),oxidation(M)。母离子电量限制在+2,+3,+4,+5和+6。同时使用正反库交叉搜索,肽段水平设定1%错误发现率(false discovery rate,FDR)。
本发明采用无标定量,计算方法为无标记,基于强度的绝对定量方法(intensity-based absolute quantification,iBAQ)。本研究多采用FOT(fraction of total)来表示蛋白定量,即某一蛋白的iBAQ值除以该实验样本整体iBAQ总和的值。
通过该步骤获得每一胃癌蛋白样品(与每一组织样本匹配)对应的蛋白质组数据集,包括蛋白的种类、数量和各种蛋白的定量值。实例中84对样本得到配对的84对蛋白质组数据。
五、胃癌分子分型
对全部样本(实例84对168例样本)对应的蛋白定量FOT值进行了筛查,要求定量的蛋白含有至少两条专属肽段(unique peptide),并且离子分值(ion score)大于等于20。通过筛选,实例中获得了9186个蛋白在84对样本的癌与癌旁组织中的表达值(表1)。
表1:9186个蛋白在84对样本的癌(T)与癌旁(N)组织中的表达值(FOT值,×10-5)
通过分析蛋白表达值的变异系数,发现当蛋白表达值为10-5时,蛋白的变异系数可降到30%左右。9186个蛋白中的6143个蛋白在至少1个组织表达谱中FOT值大于10-5。为了采用有共性的蛋白进行分子分型,进一步筛选出3619个蛋白在至少28个样本(全部样本的六分之一)中FOT值大于10-5,将FOT值小于10-5的值用10-5替换,其中2538个蛋白在至少8个样本(全部样本的十分之一)中T/N大于3倍,或者小于1/3倍。T代表癌组织样本中的蛋白定量值,N代表癌旁组织样本中的蛋白定量值。
根据上述筛选得到的2538个蛋白,对T/N的值取以10为底的log变换(参见表2),采用一致性聚类(consensus clustering)对84对样本的T/N的蛋白质组定量数据进行分子分型。
表2:2538个蛋白在84对样本的癌(T)与癌旁(N)组织中蛋白FOT值比值以10为底的对数(log10(T/N))
根据上述筛选得到的2538个蛋白在84对样本中的log10(T/N)(表2)进行分子分型。分型采用一致性聚类算法(consensus clustering),使用ConsensusClusterPlus的R(https://www.r-project.org/)程序包(Monti S(2003)Machine Learning 52:91-118,Wilkerson M.D.(2010)Bioinformatics 26(12):1572-3)。分型数目的选取如下:如果大部分样本对的一致性指数(consensus index)在0附近(图3),那么CDF曲线下面积就会显著增加,而随着分类数目的增加,CDF曲线下面积增加越来越不明显(图4),部分是由于新的分类中样本数太少,或者样本分类不稳定。因此综合考虑选取分型数目k=3,一致性聚类图参见图2,因此,最终得到了弥漫型胃癌三个亚型PX1、PX2和PX3(如图2中1/2/3所示)。
六、胃癌分型的临床特征和分子特征
根据蛋白质组表达谱数据获得的弥漫型胃癌三个亚型PX1-PX3与现有的TNM分期无关(P=0.266)。三个亚型与肿瘤原发部位有关(P=0.000),一半以上的PX3样本原发部位在胃窦区。此外,三个亚型与癌旁组织性质有关(P=0.001),3/4以上的PX2样本的癌旁只是有浅表性胃炎,不是慢性萎缩性胃炎或肠化生(见表3)。统计检验均采用Fisher’s exacttest。三个亚型的总生存时间存在差异(见图5),统计检验采用Kaplan-Meier analysis。
表3:亚型与临床指标的相关性
三个亚型的临床特征和分子特征均有明显差异。临床特征方面差异:(1)三个亚型的肿瘤原发部位差异显著(Fisher’s精确检验P=0.000),一半以上的PX3样本原发部位在胃窦区;(2)三个亚型的癌旁组织性质差异显著(Fisher’s精确检验P=0.001),3/4以上的PX2样本的癌旁只是有浅表性胃炎,不是慢性萎缩性胃炎或肠化生;(3)三个亚型的总生存时间差异显著(Kaplan–Meier分析P=0.038)。
三个亚型的分子特征方面的差异:(1)PX1和PX2中癌和癌旁差异表达的蛋白富集在细胞周期和DNA复制相关通路,而PX3中癌和癌旁差异表达的蛋白未富集在细胞周期或DNA复制相关通路中;(2)PX1-3中癌和癌旁差异表达的蛋白均富集在E2F靶基因和MYC靶基因的基因集中,但PX1在这些通路中的富集程度更高;(3)PX1中癌和癌旁差异表达的蛋白富集在RNA处理和DNA修复相关通路,而PX2和PX3中癌和癌旁差异表达的蛋白未富集在这些通路中;(4)PX2和PX3中癌和癌旁差异表达的蛋白富集在上皮-间充质转化过程、补体、细胞外基质和干扰素gamma响应的相关通路中,而PX1中癌和癌旁差异表达的蛋白未富集在这些通路中;(5)PX3中癌和癌旁差异表达的蛋白富集在NFKB介导的TNF-alpha信号通路、免疫响应和细胞黏附的相关通路中,而PX1和PX2中癌和癌旁差异表达的蛋白未富集在这些通路中(图6)。通路富集分析均采用Fisher’s精确检验,通路信息来自数据库MsigDB(http://software.broadinstitute.org/gsea/msigdb)。
确定第一亚型为细胞周期型弥漫型胃癌PX1,第二亚型为上皮-间充质转化型弥漫型胃癌PX2,第三亚型为免疫富集型弥漫型胃癌PX3。
七、胃癌分型的差异表达蛋白
根据差异表达的FDR p value和高表达3倍以上的频率与低表达1/3倍以上的频率的差值所占的比例两个指标来定义差异表达蛋白。
首先采用SAM软件计算弥漫型胃癌三个亚型中每个亚型T与N组织表达谱中FDR<0.01的蛋白,其中(高表达3倍以上频率-低表达1/3以上频率)/总检出频率>0.5的蛋白为显著高表达蛋白,而(高表达3倍以上频率-低表达1/3以上频率)/总检出频率<-0.5的蛋白为显著低表达蛋白。结果第一个亚型PX1中获得了580个显著高表达的蛋白和8个显著低表达的蛋白;第二个亚型PX2中获得了1016个显著高表达的蛋白和572个显著低表达的蛋白;第三个亚型PX3中获得了457个显著高表达的蛋白和551个显著低表达的蛋白。三个亚型共3184个显著高表达(上调)的蛋白或显著低表达(下调)的蛋白去冗余后得到2264个蛋白列于表4。
表4:弥漫型胃癌亚型的2264个差异表达蛋白
八、用于胃癌分子分型的标志物
用于胃癌分子分型的分类器的建立,包括发现、验证及测试三个阶段,即将样本FOT的log10(T/N)数据(见表2)随机分成两个独立数据集,在第一个独立数据集中进行训练、验证获得分类器,得到用于弥漫型胃癌分子分型的蛋白标志物,然后采用第二个独立数据集进行测试。实施例中将弥漫型胃癌84对样本随机分为两个数据集,第一个数据集中包含随机抽取的62对样本(75%);剩下的22对样本(25%)作为第二个数据集。
对第一个数据集中的62对样本采用Prediction Analysis of Microarray(PAM)算法(Narashiman and Chu(2002)PNAS 99:6567-6572)进行分类器的构建。PAM采用重心邻近法(nearest centroid),蛋白在每个亚型中的重心值=蛋白在每一个亚型中的平均值/蛋白的组内方差。用于新样本的分类时,先提取新样本的蛋白表达谱,然后计算这一例样本与三个亚型重心的秩相关系数(Spearman’s rank correlation),将这例样本归入秩相关系数最高的亚型。PAM在估计错误率时采用10倍交叉验证法,首先将62对样本随机分为10等分,这10等分中的三个亚型基本也是平均分布的。用90%的样本进行建模,用剩下的10%样本进行测试,如此重复10次,计算测试10次的错误率的平均值。采用2264个差异蛋白(表4)和2538个IFOT大于1的蛋白(表2)的交集共875个蛋白进行计算。软阈值(soft threshold)用于筛选蛋白集合,当阈值为0时,是初始的875个蛋白构成的集合。随着阈值的增大,蛋白集合中包含的蛋白越来越少,错误率逐渐增大(图7,采用PAM筛选蛋白标志物时采用的阈值及此阈值下的平均分型错误率)。考虑到平均错误率和最终获得的蛋白集合的平衡,选取了平均错误率为9.6%(图7中三角形标识,见表6-1)的199个蛋白(表5)建立的模型作为分型模型。
将这199个蛋白的构成的模型用于预测第二个数据集,结果平均错误率为4.5%(见表6-2)。
表5:用于分类器的199个蛋白
表6-1:第一个数据集(62例样本)的分型错误率
表6-2:第二个数据集(22例样本)的分型错误率
因此可确定该199个蛋白作为用于弥漫型胃癌分子分型的标志物,其中:
蛋白RPL23、RPL21、RPS15、RPL38、RPL27、RPS27A、CAMP、RPL36、RPL27A、AZU1、RPL7A、RPL30、MPO、RPL9、PRTN3、RPS20、SBDS、ELANE、MMP9、RPL31、FKBP3、RPL24、RPL10、S100A8、SNRPD2、RPL37A、RPS18、NDUFA6、LTF、RPL22、NDUFS7、BCCIP、S100A9、ATP5C1、HINT2、RBBP9、SRP14、ATP5O、C4BPA、C1R、RPL8、ISOC2可作为PX1亚型的标志物;
蛋白SLC4A2、ATP4B、ATP4A、DBT、CLIC6、VIL1、CKM、GIF、AUH、NNT、IDH2、GLUL、ECI2、BCKDHA、PDIA2、FABP1、C14orf159、ECHDC3、CKMT2、ACSS1、BCKDHB、NDUFB7、CYB5R1、MPC2、MFSD4、NDUFB11、REG3A、APOO、SLC25A4、MLYCD、TXNRD2、ACAT1、SMIM20、TMEM126A、ACADSB、SERPINB5、PCCB、C14orf2、GOT1、CA9、ALDH4A1、CHCHD4、APOOL、FABP3、ADCK3、ATP5F1、ME3、DCXR、Marc2、TPD52L1、NDUFB5、AQP4、CASP1、COA3、FAM162A、LDHD、COX1、NDUFA4、HIP1R、COX4I1、CRAT、PDHX、COQ10B、ATP5L、ADA、ETFDH、BCAT2、MTFP1、NDUFB6、ACADS、TGFBI、GBAS、ND4、FAHD2A、NDUFS7、TMEM214、CTH、COQ9、NDUFA11、ND1、NDUFA3、ALDH6A1、COX17、DLAT、NDUFA9、NUDT19、PPIF、SDHD、PCCA、ALDH5A1、ND5、FHL2、COX6C、ATP5D、PDHB、ALDH1L2、ATPAF2、MCCC2、BDH1、NDUFS2、CLCA1、HINT2、ACO2、FBP2、NDUFS4、OGDH、TPPP3、MTX2、ATP5A1、FGD4、CARKD、PDLIM1、RAP1GAP、COX3、SUOX、REG4、ATP5B、MCCC1、COX7A2L、COX2、NDUFC2、OPLAH、NDUFA7、L2HGDH、AIFM1、PDLIM3、VCAN、OXCT1、PC、CYTB、PYCARD、NADK2、NDUFAB1、C21orf33、SLC25A20、ANPEP、ECHDC2可作为PX2亚型的标志物;
蛋白ATP4A、ATP4B、CLIC6、SLC4A2、GIF、GLUL、CKM、MPO、PRTN3、NNT、AZU1、IDH2、MLYCD、DBT、MMP9、C14orf2、CAMP、C14orf159、PDIA2、NDUFB11、ATP5F1、CKMT2、AUH、MPC2、ALDH4A1、ATP5L、NDUFB5、FABP3、MFSD4、ND4、COX4I1、CYB5R1、ND1、SLC25A4、BCKDHA、SDHD、MMP8、CA9、NDUFB6、S100A12、APOO、BCKDHB、COA3、TXNRD2、TPD52L1、BPI、ME3、HIP1R、ELANE、MTFP1、COX1、ITGAM、PGLYRP1、NDUFA11、TMEM126A、AQP4、VBP1、RPL27A、EPCAM、REG4、SF3A2、RPL22、GSDMB、RBP2、RPL30、SRP14、SNRPE、EIF2AK2、RPL9、RPS27A、SBDS、ACAT2、ME1、TPPP3、DCTPP1、RPL27、UGDH、FKBP3、CNN1、LSM4、OTC、RBM39、REG3A、RBM14、RPL23、HSD17B2、SERPINB5、FABP1、VIL1可作为PX3亚型的标志物。
以下结合实施例说明本发明。
实施例1:
本实施例说明建立胃癌患者癌与癌旁的蛋白质组表达谱数据集的过程。建立数据集的过程包括:
1、胃癌样本收集:选取术前未接受新辅助放化疗的弥漫型胃癌患者的癌与癌旁组织,要求肿瘤细胞含量超过50%并且癌旁组织中不含肿瘤细胞。收集的癌组织与癌旁组织统称胃癌样本。
2、胃癌蛋白样品制备。
(1)称取空EP管重量并做好标记。每一对癌组织(Tumor,T)及配对的癌旁正常组织(Nearby tissue,N)样本保证同时处理。组织样本由-80℃取出后,首先在EP管内称重并核对、记录,埋入冰盒。在EP管中加入约5倍体积8M Urea蛋白裂解液,同时加入PMSF,使用一次性研磨棒碾碎,随后进行一分钟超声裂解(3秒开始,3秒间歇,振幅25%),之后置于冰上裂解15分钟,期间不时弹动。
(2)裂解液于14000g,15℃条件下离心10分钟,吸取上清至新EP管,记录上清体积,即为全蛋白裂解液。
(3)考马斯亮蓝(Bradford)法蛋白浓度测定,记录蛋白浓度。每例样本取相当于100μg全蛋白的裂解液,余下部分存于-80℃冰箱,登记入库。
(4)在100μg全蛋白裂解液中,加入10mM DTT,在56℃摇床震荡30分钟,冷却至室温后加入10mM IAA,置于避光条件下30分钟,完成还原烷基化过程。
(5)在完成还原烷基化的裂解液中加入5倍体积的pH8.2 Tris-Hcl,终浓度1mM的Cacl2,使溶液呈弱碱性。按照胰酶:蛋白为1:50的比例(质量比),确定用于消化蛋白的胰酶量,在裂解液中加入3/4胰酶的量,全部转入0.5ml离心管中,37℃旋转消化过夜。次日添加剩余1/4胰酶,消化4小时。消化完成后加入1/100体积甲酸终止消化,并离心取上清,于真空泵中抽干,得到消化后肽样品。每一组织样本对应一个肽样品。
3、胃癌蛋白样品的质谱检测。
(1)sRP液相分离
肽样品首次分离,使用特别制作的用Tip头填充的反向C18分离柱。该反向C18分离柱是在Tip头底部制作C18筛板,将C18填料悬浮于乙腈并填充于Tip头中,制作为分离柱。使用50%乙腈对该分离柱进行两次冲洗,pH=10的缓冲液(10mM碳酸氢铵)平衡分离柱,同时使用pH=10的缓冲液对步骤二中抽干的肽样品复溶,复溶后的肽样品重复四次载入分离柱,之后采用逐渐升高的梯度比例乙腈溶液进行肽样品洗脱,洗脱乙腈浓度依次为6%、9%、12%、15%、18%、21%、25%、30%和35%(体积浓度)。得到的9个浓度洗脱溶液合并为6组分的肽样品,真空抽干后用于质谱检测。
(2)肽样品质谱检测
每个组分的肽样品使用Orbitrap Fusion串联Easy-nLC 1000液相色谱串联质谱系统进行检测,并得到该肽样品对应的质谱数据。
抽干的肽样品使用A液(0.1%甲酸)复溶,上样至trap柱(100μm×2cm,粒子大小,3μm;孔径大小,),后在分析柱上进行分离(150μm×12cm,粒子大小,1.9μm;孔径大小,),梯度为3%-30%的B液(0.1%甲酸溶于乙腈)洗脱,流速600nl/min,合计洗脱75分钟。
质谱检测参数设置如下,母离子使用Orbitrap进行扫描,扫描范围为质荷比(m/z)300-1400,使用120000(200m/z)的分辨率。母离子使用top-speed模式,在Quadrupole进行选择后(1.6m/z分离窗隔离窗口)的母离子使用HCD(higher energy collision induceddissociation,高能碰撞诱导解离)进行碎裂,碎裂能量35%。子离子碎片在线性离子阱中进行检测,质谱参数设置如下:在一级全扫描时,离子阱中允许进入的最大离子数(automatic gain control target,AGC target)为5e5 ions,最大注入时间50ms,MS/MS二级扫描AGC target为5e3,最大注入时间35ms,动态排除设定18秒。液相色谱串联质谱系统使用Xcalibur软件控制进行数据采集。
4、胃癌蛋白样品的质谱数据分析。
质谱鉴定得到的图谱,使用搭载Mascot 2.3的Proteome Discover 1.4搜库,数据库使用人源NCBI(human National Center for Biotechnology Information)蛋白数据库(04-07-2013版本,共32015个蛋白种类)。
搜库参数设置如下:
Fusion:母离子质量偏差范围质谱低限(mass tolerances)为20个ppm,子离子为0.5Da。
最大允许两个漏切位点,固定修饰为carbamidomethyl(C),动态修饰为proteinacetyl(protein N-term),oxidation(M)。母离子电量限制在+2,+3,+4,+5和+6。同时使用正反库交叉搜索,肽段水平设定1%错误发现率(false discovery rate,FDR)。
本发明采用无标定量,计算方法为无标记,基于强度的绝对定量方法(intensity-based absolute quantification,iBAQ)。本研究多采用FOT(fraction of total)来表示蛋白定量,即某一蛋白的iBAQ值除以该实验样本整体iBAQ总和的值。
通过该步骤获得每一蛋白样品(与每一组织样本匹配)对应的蛋白质组数据集,包括蛋白的种类、数量和各种蛋白的定量值。
实施例1的实施结果参见前述第一至第四部分。
实施例2:
本实施例说明如何采用分类器对弥漫型胃癌患者进行分子分型。
分子分型的过程包括:
1、采集弥漫型胃癌的癌与癌旁的蛋白表达谱数据。
2、对全部样本对应的蛋白定量FOT值进行筛查,要求定量的蛋白含有至少两条专属肽段(unique peptide),并且离子分值(ion score)大于等于20。
3、将FOT值小于10-5的值用10-5替换。
4、获取每个蛋白FOT的log10(T/N)的值。
5、将这199个蛋白的构成的模型,即199个蛋白分别在三个亚型中的重心值(表5)用于预测:计算新样本与三个亚型重心的秩相关系数(Spearman’s rank correlation),将这例样本归入秩相关系数最高的亚型。例如新样本与三个亚型重心的秩相关系数分别为:PX1:0.2,PX2:0.8和PX3:0.1,则新样本归于PX2。
实施例2的实施结果参见前述第五至第八部分。

Claims (6)

1.蛋白标志物在用于制备弥漫型胃癌分子分型的分类器中的应用,其中,所述弥漫型胃癌的分子分型包括弥漫型胃癌PX1亚型、弥漫型胃癌PX2亚型及弥漫型胃癌PX3亚型,
所述弥漫型胃癌PX1亚型的蛋白标志物选自如下任意多种蛋白:RPL23、RPL21、RPS15、RPL38、RPL27、RPS27A、CAMP、RPL36、RPL27A、AZU1、RPL7A、RPL30、MPO、RPL9、PRTN3、RPS20、SBDS、ELANE、MMP9、RPL31、FKBP3、RPL24、RPL10、S100A8、SNRPD2、RPL37A、RPS18、NDUFA6、LTF、RPL22、NDUFS7、BCCIP、S100A9、ATP5C1、HINT2、RBBP9、SRP14、ATP5O、C4BPA、C1R、RPL8、ISOC2;
所述弥漫型胃癌PX2亚型的蛋白标志物选自如下任意多种蛋白:SLC4A2、ATP4B、ATP4A、DBT、CLIC6、VIL1、CKM、GIF、AUH、NNT、IDH2、GLUL、ECI2、BCKDHA、PDIA2、FABP1、C14orf159、ECHDC3、CKMT2、ACSS1、BCKDHB、NDUFB7、CYB5R1、MPC2、MFSD4、NDUFB11、REG3A、APOO、SLC25A4、MLYCD、TXNRD2、ACAT1、SMIM20、TMEM126A、ACADSB、SERPINB5、PCCB、C14orf2、GOT1、CA9、ALDH4A1、CHCHD4、APOOL、FABP3、ADCK3、ATP5F1、ME3、DCXR、Marc2、TPD52L1、NDUFB5、AQP4、CASP1、COA3、FAM162A、LDHD、COX1、NDUFA4、HIP1R、COX4I1、CRAT、PDHX、COQ10B、ATP5L、ADA、ETFDH、BCAT2、MTFP1、NDUFB6、ACADS、TGFBI、GBAS、ND4、FAHD2A、NDUFS7、TMEM214、CTH、COQ9、NDUFA11、ND1、NDUFA3、ALDH6A1、COX17、DLAT、NDUFA9、NUDT19、PPIF、SDHD、PCCA、ALDH5A1、ND5、FHL2、COX6C、ATP5D、PDHB、ALDH1L2、ATPAF2、MCCC2、BDH1、NDUFS2、CLCA1、HINT2、ACO2、FBP2、NDUFS4、OGDH、TPPP3、MTX2、ATP5A1、FGD4、CARKD、PDLIM1、RAP1GAP、COX3、SUOX、REG4、ATP5B、MCCC1、COX7A2L、COX2、NDUFC2、OPLAH、NDUFA7、L2HGDH、AIFM1、PDLIM3、VCAN、OXCT1、PC、CYTB、PYCARD、NADK2、NDUFAB1、C21orf33、SLC25A20、ANPEP、ECHDC2;
所述弥漫型胃癌PX3亚型的蛋白标志物选自如下多种蛋白:ATP4A、ATP4B、CLIC6、SLC4A2、GIF、GLUL、CKM、MPO、PRTN3、NNT、AZU1、IDH2、MLYCD、DBT、MMP9、C14orf2、CAMP、C14orf159、PDIA2、NDUFB11、ATP5F1、CKMT2、AUH、MPC2、ALDH4A1、ATP5L、NDUFB5、FABP3、MFSD4、ND4、COX4I1、CYB5R1、ND1、SLC25A4、BCKDHA、SDHD、MMP8、CA9、NDUFB6、S100A12、APOO、BCKDHB、COA3、TXNRD2、TPD52L1、BPI、ME3、HIP1R、ELANE、MTFP1、COX1、ITGAM、PGLYRP1、NDUFA11、TMEM126A、AQP4、VBP1、RPL27A、EPCAM、REG4、SF3A2、RPL22、GSDMB、RBP2、RPL30、SRP14、SNRPE、EIF2AK2、RPL9、RPS27A、SBDS、ACAT2、ME1、TPPP3、DCTPP1、RPL27、UGDH、FKBP3、CNN1、LSM4、OTC、RBM39、REG3A、RBM14、RPL23、HSD17B2、SERPINB5、FABP1、VIL1;
其中PX1为细胞周期型弥漫型胃癌,PX2为上皮-间充质转化型弥漫型胃癌,PX3为免疫富集型弥漫型胃癌。
2.根据权利要求1所述的应用,其特征在于,所述弥漫型胃癌PX1亚型、所述弥漫型胃癌PX2亚型及所述弥漫型胃癌PX3亚型存在以下临床特征差异或分子特征差异:
所述临床特征差异包括:(1)三个亚型的肿瘤原发部位差异显著,一半以上的PX3样本原发部位在胃窦区;(2)三个亚型的癌旁组织性质差异显著,3/4以上的PX2样本的癌旁只是有浅表性胃炎,不是慢性萎缩性胃炎或肠化生;(3)三个亚型的总生存时间差异显著;
所述分子特征差异为:
(A)PX1和PX2中癌和癌旁差异表达的蛋白富集在细胞周期和DNA复制相关通路,而PX3中癌和癌旁差异表达的蛋白未富集在细胞周期或DNA复制相关通路中;
(B)PX1-PX3中癌和癌旁差异表达的蛋白均富集在E2F靶基因和MYC靶基因的基因集中,但PX1在这些通路中的富集程度更高;
(C)PX1中癌和癌旁差异表达的蛋白富集在RNA处理和DNA修复相关通路,而PX2和PX3中癌和癌旁差异表达的蛋白未富集在这些通路中;
(D)PX2和PX3中癌和癌旁差异表达的蛋白富集在上皮-间充质转化过程、补体、细胞外基质和干扰素gamma响应的相关通路中,而PX1中癌和癌旁差异表达的蛋白未富集在这些通路中;
(E)PX3中癌和癌旁差异表达的蛋白富集在NFKB介导的TNF-alpha信号通路、免疫响应和细胞粘附的相关通路中,而PX1和PX2中癌和癌旁差异表达的蛋白未富集在这些通路中。
3.一种弥漫型胃癌分子分型的分类器的构建方法,其特征在于,所述构建方法包括:
1)采集数量为1/2S对来自弥漫型胃癌患者的癌与癌旁组织样本的蛋白质组表达谱数据集,包括全部样本S归集的蛋白质的数量、种类以及各种蛋白的表达量FOT值;S为包括胃癌组织样本和胃癌旁组织样本全部样本的数量,S/2为患者人数,即样本人数;
2)对全部样本S对应的各种蛋白的表达量FOT值进行筛查,要求定量的蛋白含有至少两条专属肽段(unique peptide),并且离子分值(ion score)大于等于20,获得初筛数量Q1的蛋白在全部样本中的FOT值;
3)从Q1数量的蛋白中选择在至少1/6×S个样本中FOT值大于10-5的数量为Q2的蛋白,并将FOT值小于10-5的值用10-5替换;进一步从Q2数量的蛋白中筛选在至少1/10×S/2个样本人数中T/N大于3倍,或者小于1/3倍的蛋白,共得到数量为Q3个蛋白;
4)获取上述数量为Q3蛋白FOT的T/N的值,并对T/N的值取以10为底的对数变换;
5)采用一致性聚类(consensus clustering)对上述筛选得到的Q3个蛋白在S样本中的log10(T/N)进行分子分型,得到蛋白质组学技术对弥漫型胃癌的3个亚型;
6)根据差异表达的FDR p value和高表达3倍以上的频率与低表达1/3倍以上的频率的差值所占的比例两个指标来定义差异表达蛋白,从初筛数量Q1的蛋白中获取数量为Q4弥漫型胃癌亚型的差异表达蛋白;
7)选取上一步骤中获得的数量Q4的差异表达蛋白和用于一致性聚类的数量Q3的蛋白的交集数量为Q5的蛋白,采用PAM进行分类器的构建,选取分型平均错误率小于某一定值的数量为Q6蛋白作为胃癌分子分型蛋白标志物,建立弥漫型胃癌分子分型的分类器。
4.根据权利要求3所述的构建方法,其中步骤7)中所述PAM采用重心邻近法(nearestcentroid),蛋白在每个亚型中的重心值=蛋白在每一个亚型中的平均值/蛋白的组内方差。
5.根据权利要求3或4所述的构建方法,其中,
步骤1)中弥漫型胃癌患者的癌与癌旁组织样本数量为84对,即1/2S样本人数为84,则全部样本S数量为168;
步骤2)中初筛蛋白数量Q1为9186个;
步骤3)中蛋白数量Q2为3619个,1/6×S个样本为28个,1/10×S/2个样本人数为8个;所述蛋白数量Q3为2538个;
步骤5)中确定了弥漫型胃癌的3个亚型,分别命名为PX1、PX2和PX3;
步骤6)中差异表达蛋白数量Q4为2264个;
步骤7)中交集蛋白数量Q5为875个,弥漫型胃癌的3个亚型的分型平均错误率为9.6%,确定的弥漫型胃癌分子分型蛋白标志物数量Q6为199个。
6.根据权利要求3所述的构建方法,其中所述3个亚型PX1、PX2和PX3存在以下临床特征差异或分子特征差异:
所述临床特征差异包括:(1)三个亚型的肿瘤原发部位差异显著,一半以上的PX3样本原发部位在胃窦区;(2)三个亚型的癌旁组织性质差异显著,3/4以上的PX2样本的癌旁只是有浅表性胃炎,不是慢性萎缩性胃炎或肠化生;(3)三个亚型的总生存时间差异显著;
所述分子特征差异为:
(A)PX1和PX2中癌和癌旁差异表达的蛋白富集在细胞周期和DNA复制相关通路,而PX3中癌和癌旁差异表达的蛋白未富集在细胞周期或DNA复制相关通路中;
(B)PX1-PX3中癌和癌旁差异表达的蛋白均富集在E2F靶基因和MYC靶基因的基因集中,但PX1在这些通路中的富集程度更高;
(C)PX1中癌和癌旁差异表达的蛋白富集在RNA处理和DNA修复相关通路,而PX2和PX3中癌和癌旁差异表达的蛋白未富集在这些通路中;
(D)PX2和PX3中癌和癌旁差异表达的蛋白富集在上皮-间充质转化过程、补体、细胞外基质和干扰素gamma响应的相关通路中,而PX1中癌和癌旁差异表达的蛋白未富集在这些通路中;
(E)PX3中癌和癌旁差异表达的蛋白富集在NFKB介导的TNF-alpha信号通路、免疫响应和细胞粘附的相关通路中,而PX1和PX2中癌和癌旁差异表达的蛋白未富集在这些通路中。
CN202210562342.7A 2017-03-31 2018-02-13 蛋白标志物在制备弥漫型胃癌分子分型的产品中的应用及弥漫型胃癌分子分型的分类器 Active CN115266962B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210562342.7A CN115266962B (zh) 2017-03-31 2018-02-13 蛋白标志物在制备弥漫型胃癌分子分型的产品中的应用及弥漫型胃癌分子分型的分类器

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN2017102059357 2017-03-31
CN201710205935 2017-03-31
CN202210562342.7A CN115266962B (zh) 2017-03-31 2018-02-13 蛋白标志物在制备弥漫型胃癌分子分型的产品中的应用及弥漫型胃癌分子分型的分类器
CN201810149737.8A CN108445097A (zh) 2017-03-31 2018-02-13 弥漫型胃癌的分子分型、用于分型的蛋白标志物及其筛选方法和应用

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN201810149737.8A Division CN108445097A (zh) 2017-03-31 2018-02-13 弥漫型胃癌的分子分型、用于分型的蛋白标志物及其筛选方法和应用

Publications (2)

Publication Number Publication Date
CN115266962A CN115266962A (zh) 2022-11-01
CN115266962B true CN115266962B (zh) 2024-05-31

Family

ID=63192337

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202210562342.7A Active CN115266962B (zh) 2017-03-31 2018-02-13 蛋白标志物在制备弥漫型胃癌分子分型的产品中的应用及弥漫型胃癌分子分型的分类器
CN201810149737.8A Pending CN108445097A (zh) 2017-03-31 2018-02-13 弥漫型胃癌的分子分型、用于分型的蛋白标志物及其筛选方法和应用

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201810149737.8A Pending CN108445097A (zh) 2017-03-31 2018-02-13 弥漫型胃癌的分子分型、用于分型的蛋白标志物及其筛选方法和应用

Country Status (1)

Country Link
CN (2) CN115266962B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109182516A (zh) * 2018-09-07 2019-01-11 河南大学 基于smim20基因的消化道癌辅助诊断试剂盒及其应用
CN109613254B (zh) * 2018-11-06 2022-04-05 上海市公共卫生临床中心 一种用于肿瘤治疗和诊断的靶点标志物pdia2
CN110146636B (zh) * 2019-04-30 2022-05-13 北京谷海天目生物医学科技有限公司 胃癌分型的蛋白标志物的筛选方法、筛选装置及其筛选的蛋白标志物的应用
CN110144005B (zh) * 2019-05-29 2021-07-16 兰州大学 与黄芩素抗癌活性相关的两个氧化磷酸化通路关键蛋白
CN110404074B (zh) * 2019-08-20 2021-07-13 中国医学科学院基础医学研究所 Ogdh抑制物在治疗病毒感染性疾病中的应用
CN111281967B (zh) * 2020-04-07 2023-03-10 复旦大学附属中山医院 Hint2在制备治疗缺血性血管疾病药物中的应用
CN111933211B (zh) * 2020-06-28 2023-10-31 北京谷海天目生物医学科技有限公司 癌症精准化疗分型标志物筛选方法、化疗敏感性的分子分型方法和应用
CN112071363B (zh) * 2020-07-21 2023-11-14 北京谷海天目生物医学科技有限公司 胃黏膜病变蛋白质分子分型、病变进展及胃癌相关蛋白标志物、预测病变进展风险的方法
CN112837815A (zh) * 2021-02-04 2021-05-25 复旦大学附属中山医院 基于Lauren分型和术后残留淋巴结评估胃癌预后的模型及应用
CN113555121B (zh) * 2021-08-25 2024-05-28 南方医科大学南方医院 一种胃癌预后标志物的筛选和分类方法、胃癌预后标志物和检测胃癌预后的试剂及应用
CN113466471B (zh) * 2021-09-06 2021-11-26 中国农业科学院农产品加工研究所 胃粘膜损伤相关生物标志物、筛选方法及用途
CN114694748B (zh) * 2022-02-22 2022-10-28 中国人民解放军军事科学院军事医学研究院 一种基于预后信息与强化学习的蛋白质组学分子分型方法
CN116798520B (zh) * 2023-06-28 2024-08-23 复旦大学附属肿瘤医院 鳞状细胞癌组织起源位点蛋白标志物预测模型的构建方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105039523A (zh) * 2008-11-17 2015-11-11 威拉赛特公司 用于疾病诊断的分子表达谱的方法和组合物

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1908189A (zh) * 2005-08-02 2007-02-07 博奥生物有限公司 体外辅助鉴定肠型胃癌及其分化程度的方法与专用试剂盒
CN101329348A (zh) * 2007-06-18 2008-12-24 许洋 检测胃癌特征蛋白的优化质谱模型及其制备方法和应用
AT505726A2 (de) * 2007-08-30 2009-03-15 Arc Austrian Res Centers Gmbh Set von tumor-markern
CN101685080B (zh) * 2008-06-12 2012-10-31 高尚先 一种生物样品中蛋白质组的分析方法
CN102348979A (zh) * 2009-03-09 2012-02-08 乔治亚大学研究基金公司 胃癌诊断用蛋白标记的鉴定
TWI399541B (zh) * 2009-05-27 2013-06-21 Univ Nat Taiwan 檢測一胃癌預後程度的方法
US20110059452A1 (en) * 2009-08-18 2011-03-10 Vanderbilt University Methods of screening for gastric cancer
CN102435665A (zh) * 2011-09-23 2012-05-02 浙江省新华医院 胰腺癌早期诊断的血清肿瘤标志物、检测方法和诊断模型
US10648035B2 (en) * 2012-11-26 2020-05-12 The Johns Hopkins University Methods and compositions for diagnosing and treating gastric cancer
GB201320061D0 (en) * 2013-11-13 2013-12-25 Electrophoretics Ltd Materials nad methods for diagnosis and prognosis of liver cancer

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105039523A (zh) * 2008-11-17 2015-11-11 威拉赛特公司 用于疾病诊断的分子表达谱的方法和组合物

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
胃癌分子分型研究进展;李加桩;王凯冰;陈丹;郑红艳;隋红;;中国肿瘤(第11期);888-892 *

Also Published As

Publication number Publication date
CN115266962A (zh) 2022-11-01
CN108445097A (zh) 2018-08-24

Similar Documents

Publication Publication Date Title
CN115266962B (zh) 蛋白标志物在制备弥漫型胃癌分子分型的产品中的应用及弥漫型胃癌分子分型的分类器
US20210040562A1 (en) Methods for evaluating lung cancer status
Filella et al. Prostate cancer detection and prognosis: from prostate specific antigen (PSA) to exosomal biomarkers
AU2024200059A1 (en) Methods and systems for analysis of organ transplantation
Yao et al. Identification of gene-expression signatures and protein markers for breast cancer grading and staging
CN103733065A (zh) 用于癌症的分子诊断试验
CN103299188A (zh) 用于癌症的分子诊断试验
US20170073763A1 (en) Methods and Compositions for Assessing Patients with Non-small Cell Lung Cancer
CN105981026A (zh) 生物标志标识方法及用于其的装置和试剂盒
Xiao et al. Identification of five genes as a potential biomarker for predicting progress and prognosis in adrenocortical carcinoma
CN103998622A (zh) 肺癌的多基因预后试验
Ou et al. Novel breast cancer biomarkers identified by integrative proteomic and gene expression mapping
Tian et al. Systematic identification of key genes and pathways in clear cell renal cell carcinoma on bioinformatics analysis
Cantara et al. Molecular signature of indeterminate thyroid lesions: current methods to improve fine needle aspiration cytology (FNAC) diagnosis
CN104818322B (zh) miRNA和Cyfra21‑1联合在检测非小细胞肺癌中的应用
Sacher et al. Prospective validation of rapid plasma genotyping as a sensitive and specific tool for guiding lung cancer care
Dos Santos et al. Transcriptomic analysis reveals a tissue-specific loss of identity during ageing and cancer
WO2021228888A1 (en) Biomarkers for hyperprogressive disease and therapy response after immunotherapy
Wong et al. Regional and bilateral MRI and gene signatures in facioscapulohumeral dystrophy: implications for clinical trial design and mechanisms of disease progression
Wilmott et al. Tumour procurement, DNA extraction, coverage analysis and optimisation of mutation-detection algorithms for human melanoma genomes
Nooh et al. Prediction of blood miRNA-mRNA regulatory network in gastric cancer
CN115762800A (zh) 可以预测黑色素瘤患者预后及免疫治疗应答率的评分系统
CN109207581A (zh) 一种自身免疫性疾病诊断试剂盒和应用
Elzanowska et al. Defining optimal conditions for tumor extracellular vesicle DNA extraction for mutation profiling
CN115678995A (zh) 生物标志物组合及其在预测药物治疗胃癌效果中的应用

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant