Cell Research(IF46)丨华山医院神外科赵曜教授等团队利用多组学技术首次揭示垂体腺瘤蛋白基因组学特征及7大分子分型

2022-12-08 10:36:16, XHL 上海吉凯基因医学科技股份有限公司


 


垂体神经内分泌肿瘤(pituitary neuroendocrine tumor,PitNET) ,又称为垂体腺瘤,是常见的一种颅内肿瘤,发病率大约为70-90人/10万人。垂体腺瘤的治疗受到手术无法完全切除以及化疗和放射治疗的有限性影响,其治疗目前仍面临很大挑战。


整合多组学策略的进展,如蛋白质组学和磷酸化蛋白质组学结合基因组学研究,已经推动了多种不同肿瘤的治疗发展。近期,国家神经疾病医学中心、复旦大学附属华山医院神经外科赵曜教授团队和复旦大学表型组研究员、生命科学学院丁琛教授团队对200例垂体腺瘤组织和7例对照样本进行了多组学(基因组、转录组、蛋白质组和磷酸化蛋白质组)检测和分析,以揭示垂体腺瘤的蛋白基因组学特征,并基于组学结果进行疾病分子分型,为将来的垂体腺瘤的精准治疗提供理论依据。研究成果以“Integrated proteogenomic characterization across major histological types of pituitary neuroendocrine tumors”为题,发表在Cell Research(IF:46.297)上。



多组学样本分组设置及临床特征:

基因组全外显子测序:200例垂体腺瘤组织;

转录组RNA测序:194例垂体腺瘤组织+7例垂体前叶组织;

蛋白质组:200例垂体腺瘤组织+7例垂体前叶组织;

磷酸化蛋白质组:194例垂体腺瘤组织+7例垂体前叶组织;


垂体腺瘤患者的多种TF谱系、病理亚型及其他临床特征如下:

研究结果  

01

垂体腺瘤(PitNET)蛋白基因组图景概览



全外显子测序(WES)检测到了7333个突变基因,包括11092个非沉默点突变和419个小插入或缺失。研究者发现了多个和PitNET功能相关的显著性突变基因(SMG),包括GNAS (11%)、KRT76 (8%)、TCHH (6%)、ZMIZ2 (6%)、DGKZ (6%)、KRTAP9-2 (5%)和USP8 (4%)。对体细胞碱基改变分析显示,与其他五种替代类型相比,PitNET患者携带的C→A转变的比例较高。与TCGA数据库中的其他肿瘤相比,PitNET的肿瘤突变负荷(TMB)在33种癌症类型中处于较低或中等水平。随后,研究者比较了显著性突变基因(SMG)的变异等位基因频率(VAF),发现USP8突变(中位数:0.45)排在第一位,其次是GNAS突变(中位数:0.36)。紧接着,研究者对突变在各PitNET亚型中的富集情况进行分析,结果显示USP8突变在ACTH亚型中富集,而GNAS突变在多激素亚型和GH亚型中富集。体细胞CNA(拷贝数改变)分析发现染色体臂级别的扩增(Chr 5、7、8、9、12q、14q、19、20、21p)和缺失(Chr 1p、2q、11)。



转录组RNA测序检测到了18397个基因(FPKM>1)。蛋白质组检测到了10011个蛋白,其中至少在一种肿瘤亚型的50%样本中被检测到的6160个蛋白进入下游分析。磷酸化蛋白质组检测到了5483个蛋白上的29219个磷酸化位点,其中至少在一种肿瘤亚型的50%样本中被检测到的9905个磷酸化位点进入下游分析。各组学的PCA分析显示,在蛋白质组学上,PIT1谱系(GH、PRL、TSH、silent PIT1和多激素腺瘤亚型)和APG、ACTH、silent TPIT及其他PitNET类型有明显区分。而转录组和磷酸化蛋白质组的PCA结果则显示,PIT1谱系肿瘤与APG组有更高的相似性。



02

基因改变对转录组、蛋白质组和磷酸化蛋白质组的影响



肿瘤组织转录组和蛋白质组的6115个mRNA-蛋白质对中的92.95%正相关。强相关基因富集在和神经系统、EMT、激素代谢过程相关的通路,表明这些通路在PitNET中过激活。总体样本层面上的mRNA-蛋白相关性中位值为0.45,各临床病理亚型样本的相关性从0.42-0.46。



紧接着,研究者检测了23109个体细胞CNA对同基因座位(cis效应)和其他基因座位(trans效应)上mRNA、蛋白和磷酸化蛋白丰度的影响。Cis效应分析发现,CNA分别影响6545个mRNA表达丰度和684个蛋白质表达丰度,其中有502个显著性cis效应重合。这502个基因富集在蛋白翻译后修饰、NAD+合成和代谢通路。此外,研究者发现,CNA通过cis效应影响15个肿瘤相关基因(CAG)的mRNA和蛋白质表达丰度,CNA 通过cis效应影响6个CAG的mRNA和磷酸化蛋白质表达丰度。



03

PIT1谱系肿瘤中的GNAS突变和GNAS拷贝数增加



研究者的结果中,GNAS是PIT1谱系肿瘤中频率最高的基因突变,包含两个突变热点:R186C/G/L/H和Q212L。GNAS拷贝数增加,一个受到cis效应影响的肿瘤相关基因,对PIT1谱系肿瘤有重要影响。GNAS突变(VAF>5%)或GNAS拷贝数增加均富集在PIT1谱系肿瘤中。


通过整合全外显子数据,研究者将PIT1谱系肿瘤分为三个亚群:野生型(WT)、GNAS拷贝数增加和GNAS突变。与野生型相比,有GNAS突变的样本中,GNAS蛋白水平下调,而与GNAS拷贝数增加的样本中,GNAS mRNA和蛋白水平均上调。与WT相比,在GNAS突变亚群中,参与生长素合成、分泌和作用通路的基因在mRNA、蛋白和磷酸化水平均上调。与WT相比,GNAS拷贝数增加亚群有更高的多基因增殖得分(MGPS)和临床肿瘤体积。细胞周期相关分子中,PRKDC和CDK6是前两个与GNAS拷贝数增加正相关的分子。不管是在PIT1谱系肿瘤还是所有的PitNET中,高CDK6得分的样本占比在高GNAS拷贝数增加样本中比在低GNAS拷贝数增加的样本中更多。该结果提示我们,靶向CDK6的抑制剂也许是对含有GNAS拷贝数增加的PIT1谱系肿瘤有效的治疗策略。



04

基于多组学结果的PitNET分子分型



基于蛋白质组、转录组和磷酸化蛋白质组的数据分别将PitNET分为7个、5个和7个疾病亚型。7个基于蛋白质组的疾病亚型为GHenrich、 EMTPRO、 PRLenrich、TSH/silent PIT1enrich、ACTHenrich、silent TPITenrich和SF1/NULLenrich。Hedgehog信号通路在GHenrich亚型中上调,MYC targets v1在PRLenrich亚型中上调。TSH/silent PIT1enrich中富集了干扰素α响应、抗原提呈通路。SF1/NULLenrich亚型中存在代谢相关通路的上调,如脂肪酸代谢和三羧酸循环。此外,在SF1/NULLenrich亚型中,有更多的男性、更高的患者年龄和更大的肿瘤直径。ACTHenrich亚型富集了USP8突变,ACTHenrich亚型和silent TPITenrich亚型都有更高比例的女性患者。



05

前三大TF谱系垂体腺瘤的分子特征分析



本研究中,垂体腺瘤按TF谱系分类,患者占比前三的分别是PIT1谱系、TPIT谱系和SF1谱系,总计占比超过75%。因此,研究者对这三类垂体腺瘤中的分子特征进行了分析。


(1)在PIT1谱系肿瘤中发现了以EMT(上皮间质转化)为特征的侵袭亚型

与其他基于蛋白质组的疾病亚型相比,在EMTPRO亚型中存在与稳态及EMT相关分子包括GP1BB、FGB、MMP8、FN1和ITGB3的高表达。EMTPRO亚型是高侵袭性的,且肿瘤直径大。鉴于EMT在肿瘤转移中不可忽视的作用,研究者利用转录组和蛋白质组数据分析了基质细胞在其中的贡献。结果表明,不管是RNA还是蛋白质的分析结果都表明,EMTPRO亚型都有更高的基质细胞分值,该结果得到HE染色的验证。上皮marker pan-cytokeratin和间质marker FN1的免疫荧光共染色表明,在EMTPRO亚型中有更高比例的染色面积。以上这些结果表明,EMTPRO亚型中的肿瘤细胞呈现为EMT状态。


EMT-inducing transcription factor(EMT-TF)是EMT表型的重要驱动因子。研究者比较了PIT1谱系肿瘤中基于转录组分型的两个亚型间EMT-TF的表达情况,结果表明5个EMT-TF(SNAI1、SNAI2 、ZEB2 、TWIST1 和TWIST2)在EMTRNA亚型中显著上调。这5个EMT-TF的表达模式在EMTRNA和EMTPRO中高度正相关。



(2)在TPIT谱系肿瘤中发现3种EGFR激活模式

EGFR与多种癌症相关,且被提示可以作为ACTH PitNET的治疗靶点。研究者发现在TPIT谱系中,EGFR的mRNA表达、蛋白丰度和磷酸化修饰水平均比其他类型的肿瘤要高。通过分析EGFR相关通路,将TPIT谱系肿瘤分为3个亚群,分别为USP8突变的ACTH肿瘤、USP8未突变的ACTH肿瘤和silent TPIT肿瘤。在USP8未突变的ACTH肿瘤中,EGFR配体的mRNA显著上调。在silent TPIT肿瘤中,EGFR T693磷酸化显著富集。在TPIT谱系肿瘤中,EGFR表达更高,且EGFR T693磷酸化阳性率在silent TPIT肿瘤中高于ACTH肿瘤。此外,研究者还发现,4种转录因子可能参与了ACTH分泌和调控,或可作为新的治疗靶点。


(3)SF1谱系肿瘤和NULL肿瘤存在VEGF通路和缺氧信号的激活

多组学的疾病分型分析时发现,SF1谱系肿瘤和NULL肿瘤呈现出相似的表达模式,因此合并分型为SF1/NULLenrich亚型。肿瘤相关通路的富集分析显示,在SF1/NULLenrich亚型存在VEGF通路和缺氧信号的激活。进一步分析发现,不管是在mRNA、蛋白质还是磷酸化蛋白质水平,VEGF信号的下游信号,如RAS/RAF/MEK/ERK和PI3K-AKT信号均在SF1/NULLenrich亚型中激活。相比于其他PitNET肿瘤,VEGFR2的IHC染色结果也显示,其表达量在SF1谱系肿瘤和NULL肿瘤高表达。以上结果提示我们,靶向VEGFR2的治疗方案或许是这类病人的有效治疗手段。


06

PitNET的免疫浸润特征



研究者进一步使用xCell分析了肿瘤微环境中各种细胞的占比。一致性聚类分析共得到了4种具有独特免疫和基质特征的肿瘤类型:免疫耗竭型、CD4+T细胞浸润型、内皮细胞型和CD8+T细胞浸润型。CD8+T细胞浸润型包括ACTHenrich亚型和silentTPITenrich亚型,存在EGFR信号和细胞周期通路的上调。内皮细胞型存在CSF1R、CD34和PDL2 mRNA水平的上调及FGFR1蛋白水平的上调。CD4+T细胞浸润型中存在CTLA4的上调,提示我们此类肿瘤可能对免疫检查点相关治疗手段有响应。免疫耗竭型包括TSH/silent PIT1enrich和GHenrich蛋白亚型,存在高评分的Treg细胞和PDL1 mRNA的上调。PDL1 mRNA在TSH/silent PIT1enrich和GHenrich蛋白亚型中上调,提示我们可以尝试使用PDL1疗法。EMTPRO亚型有高的免疫评分和低的PDL1,而高的免疫评分或许与JAK3-STAT6-FOS/JUN调控的细胞因子表达有关,这也为该类肿瘤的治疗提供了潜在靶点。



07

蛋白质组疾病亚型中可用药物靶点的表达情况分析



为了利用这些多组学数据发现更多的潜在PitNET治疗方案,研究者评估了FDA批准的药物靶点的表达水平,包括多巴胺受体2 (DRD2)和生长抑素受体(SSTR2/SSTR5),以及参与DNA修复的O -6-甲基鸟嘌呤-DNA甲基转移酶(MGMT),这些是已知的影响temozolomide(替莫唑胺)疗效的分子。生长抑素激动剂是治疗GH PitNET的主要药物,据报道对ACTH PitNET有效。转录组数据表明,在GH PitNET和TSH/silent PIT1enrich蛋白质组亚型中,SSTR2 和SSTR5表达升高。MGMT mRNA在PRLenrich、TSH/silent PIT1enrich和silent TPITenrich蛋白质组亚型中表达下调。DRD2被报道在PRL和GN PitNET中高表达,与本研究发现的DRD2 mRNA在PRLenrich蛋白质亚型中也显著上调相符。


08

在750个 PitNET组成的独立队列中验证前期发现的结果



研究者在一个包含750个PitNET患者,且有超过78%的患者有随访数据的队列(以下简称验证队列)中验证在先前200例患者构成的队列(以下简称发现队列)中发现的结果。验证队列的结果分析发现,高GNAS IHC分值的患者占比在高肿瘤直径组中显著高于低肿瘤直径组,而高IHC分值意味着差的预后。在验证队列中ZEB2和TWIST1在高侵袭组中有更高的IHC分值。此外,在临床病理分型的TSH和silent PIT1亚型中,PDL1都有更高的IHC分值,这与发现队列中发现的PDL1在免疫耗竭型亚群中高表达相符合。VEGFR2的IHC染色结果也表明,在发现队列和验证队列的SF1和NULL肿瘤中均存在VEGFR2的高表达。验证队列中也发现了EGFR在TPIT谱系肿瘤中的独特性表达,与前期结果相符。而EGFR T693磷酸化的IHC染色也表明,高IHC分值的占比在silent TPIT亚型中比在ACTH亚型中更高。总体以上结果表明,发现队列发现的一些诊断标志物或治疗靶点在750例的验证队列中得到验证。


研究总结

本研究对200例垂体腺瘤(PitNET)组织和7个垂体前叶组织进行了多组学(基因组、转录组、蛋白质组和磷酸化蛋白质组)检测和分析。基因组学数据表明,GNAS拷贝数增加可以作为PIT1谱系肿瘤过度增殖的一个可靠的诊断标志。基于蛋白质组的PitNET分子分型共确定了7个亚型。其中,肿瘤过表达上皮细胞转化(EMT)标记聚集成一个更具侵袭性的亚型。进一步的分析确定了针对不同亚型的潜在治疗靶点,包括CDK6、TWIST1、EGFR和VEGFR2。在750名PitNET患者的独立验证队列中进一步证实了这些已识别的分子标记和不同亚群中的变化。这种基于多组学的蛋白基因组分析大大拓展了我们对于PitNET的生理病理学理解,并提出了新的治疗靶点和策略。



吉凯基因凭借多年在靶标筛选及验证服务领域的技术积累,建立的标准化 、工程化 、系统化的GRP平台,为中国研究型医生提供科研服务,加快科研成果转化。其中,多组学平台包含蛋白质组学平台和高通量测序平台


·蛋白质组学平台拥有多台timsTOF Pro、Exploris 480高精度质谱仪,专业领先的Spectronaut Plusar、Mascot等分析软件,提供专业的4D、DIA、TMT、PRM、磷酸化修饰组、olink蛋白质组等检测服务,强大的机器学习算法、IPA分析、蛋白基因组分析服务,系统的生物标志物、分子分型、药物靶点、基因功能研究等解决方案,真正让广大研究型医生的科研工作更省心、更省力、更高效;


·高通量测序平台分为常规测序服务和单细胞测序服务:单细胞测序拥有10x和BD两个平台,提供单细胞RNA-seq、单细胞核测序、单细胞混样RNA-seq、单细胞TCR/BCR、单细胞(RNA+ATAC)、空间转录组测序等服务;常规测序服务提供meRIP-seq(m6A/m1A/m7G/m5C 等RNA甲基化修饰测序)、acRIP-seq(ac4C RNA乙酰化修饰测序)、ATAC-seq、Ribo-seq(翻译组测序) 、mRNA/miRNA/LncRNA/circRNA-seq、全转录组测序(两文库/三文库)、外泌体miRNA/LncRNA-seq、WGS/WES、WGBS、RRBS、BSAS等服务。



1.实验技术干货

2.蛋白质组学研究

3.腺病毒简介及应用

4.临床基础研究思路解析    

5.组织特异性腺相关病毒

6.单细胞测序    

7.慢病毒实验操作指南

8.悬浮细胞专用病毒

9.靶点设计/数据库教程

10.测序技术研究与应用

11.非编码RNA研究技术与应用

12.腺相关病毒选择/应用    

13.表观遗传研究

14.文章解析

15.国自然课题设计思路解析

16.生物信息分析及工具      

17.外泌体研究    

18.肿瘤免疫研究

19.高分文章  

20.吉凯病毒神经方向应用案例 



  • 客服电话: 400-6699-117 转 1000
  • 京ICP备07018254号
  • 电信与信息服务业务经营许可证:京ICP证110310号
  • 京公网安备1101085018
  • 客服电话: 400-6699-117 转 1000
  • 京ICP备07018254号
  • 电信与信息服务业务经营许可证:京ICP证110310号
  • 京公网安备1101085018

Copyright ©2007-2024 ANTPEDIA, All Rights Reserved