《Nature Communications》：Multi-omics integration unravels four molecular subgroups of corticotroph pituitary neuroendocrine tumours with distinct clinicopathological features
编辑推荐：
促肾上腺皮质激素型垂体神经内分泌肿瘤（corticotroph pituitary neuroendocrine tumours, corticotroph PitNETs）/腺瘤是具有异质性的鞍区肿瘤。目前确立的组织病理学分类方法常被认为在充分捕捉此类肿瘤的
促肾上腺皮质激素型垂体神经内分泌肿瘤（corticotroph pituitary neuroendocrine tumours, corticotroph PitNETs）/腺瘤是具有异质性的鞍区肿瘤。目前确立的组织病理学分类方法常被认为在充分捕捉此类肿瘤的临床与生物学复杂性方面存在局限。迄今，促肾上腺皮质激素型PitNETs中尚未建立基于分子的分类。研究人员汇编了270例促肾上腺皮质激素型PitNETs（111例内部、159例外部）的分子数据，涵盖表观基因组、转录组和蛋白质组谱。通过综合整合分析以识别、验证并表征明确的分子亚群。促肾上腺皮质激素型PitNETs可划分为四个稳健且与临床病理特征各异的分子亚群，这些亚群可经SSTR1、GATA3和SSTR5免疫组织化学在显微镜下广泛区分。整合上述分子亚群的分层模型显示出显著的预后效用。研究结果表明支持建立精细化的基于分子的促肾上腺皮质激素型PitNETs分类，其完整临床价值需在前瞻性研究中验证。为促进未来研究，研究人员提供了一个易于使用的表观基因组分类器用于促肾上腺皮质激素型PitNETs。
研究背景与立项依据
促肾上腺皮质激素型垂体神经内分泌肿瘤（corticotroph pituitary neuroendocrine tumours, corticotroph PitNETs）亦称促肾上腺皮质激素型垂体腺瘤，临床表现从“静默性”肿瘤到功能性库欣病（Cushing’s disease, CD）不等，术后病程在生长动力学、复发率及治疗反应上高度异质。现行WHO分类依赖TPIT、SF1、PIT1转录因子免疫表型，并将皮质otroph PitNETs分为稀疏颗粒型（sparsely granulated corticotroph tumours, SGCTs）、致密颗粒型（densely granulated corticotroph tumours, DGCTs）与Crooke细胞型（Crooke cell tumours, CCTs）；该框架被指难以反映真实肿瘤行为，且颗粒判定标准偏描述性、观察者间一致性有限。虽既往发现USP8突变等分子线索并提出三群划分，但缺乏大样本多组学验证与常规可用分类器。为此，研究人员整合内部与外部公开队列共270例，开展表观-转录-蛋白多组学整合，旨在建立可重复、具预后价值的分子分型。
主要技术方法概览
研究人员汇集111例内部（汉堡-埃彭多夫大学医学中心等）与159例外部公开队列（Bujko、Neou、Capper、Salomon、Zhang、Silva-Júnior、Belakhoua等）促肾上腺皮质激素型PitNETs，总样本270例；分子层面对应DNA甲基化（Illumina 450K/EPIC/EPICv2）、RNA-seq、LC-MS/MS蛋白质组。分析流程含gap statistic定簇数、ConsensusClusterPlus共识聚类、UMAP降维、GSEA（clusterProfiler）与单样本GSVA、拷贝数（conumee/GenVisR）、Sanger测序USP8/USP48热点、免疫组化（SSTR1/GATA3/SSTR5/TPIT/SF1/PIT1/ACTH/CAM5.2/Ki-67）双盲评分、有序逻辑回归与rpart决策树，并开发R Shiny版MagicTPiT甲基化分类器（以参考系列代表签名做共识聚类打分）。
研究结果
Corticotroph PitNETs split into four distinct molecular subgroups
基于表观组、转录组、蛋白组独立gap statistic均指向k=4，共识聚类无批次效应，三模态聚类一致性高于随机、以甲基化最稳；208例参考系列按甲基化（缺则RNA）重分为亚群1（n=61）、2（n=48）、3（n=33）、4（n=66）。
Characterisation of genomic and epigenomic features
Heidelberg脑瘤分类器确认多数属“PITAD ACTH”甲基化类但不能分四亚群；累积拷贝数剖面（cumulative copy number profiles, CCNPs）显示亚群1染色体Alterations频繁（5、7-9、12、13q、14q、20增益，2、4、6、10、15q缺失），亚群2-4相对安静。
Characterisation of transcriptome and proteome profiles
RNA-GSEA示亚群1富集分化/发育（stemness），亚群2富集突触囊泡胞吐（synapse assembly），亚群3富集免疫迁移（gonadotrophs相关），亚群4富集跨膜信号与分泌囊泡（secretory vesicle）；蛋白组佐证且亚群4线粒体蛋白偏高。
Characterisation of histomorphology
PAS/ACTH/Crooke（CAM5.2环样）半定量示亚群1与4多DGCT表型、亚群2与3多SGCT；但两观察员组织亚型判定Cohen’s kappa仅0.577，证明WHO亚型与分子亚群映射弱。
SSTR1, GATA3 and SSTR5 expression associates with corticotroph molecular subgroups
SSTR1 mRNA/蛋白在亚群2最高（>50%胞质阳性者80%），GATA3在亚群3最高（核阳性85.7%），SSTR5在亚群4最高（膜阳性76.2%），三者丰度染色互斥且特异，观察员kappa=0.961，优于组织颗粒分型。
Validation of four corticotroph molecular subgroups
62例验证系列（38内部+24 Belakhoua）甲基化再证k=4、亚群1拷贝数紊乱复现、SSTR1/GATA3/SSTR5特异性复现；参考+验证UMAP合并仍稳分四群，治疗史不影响。
Corticotroph molecular subgroups relate to distinct patient and tumour features
合并270例：亚群4最年轻（中位40岁）、女性为主（93.1%）、多为显性CD（89.9%）、微肿瘤比例最高（51.2%）、侵袭最低（23.5%）、Ki-67≥3%达47.7%；亚群2全女、多静默；亚群3静默为主、具gonadotroph签名；亚群1大瘤/静默或功能混合、拷贝数不稳、USP8变异近无；USP8变异74.5%落亚群4，USP48变异在USP8野生中仅见于亚群4（50%）。
Incorporation of corticotroph molecular subgroups into an integrated prognostic stratification model
有序回归（106例随访）示单变量中侵袭（OR=0.033）、亚群1（OR=0.14）为负向预后，女性/微瘤/USP8变异/全切（GTR）为正向；多变量仅侵袭、亚群1、GTR独立。rpart决策树以侵袭→GTR→亚群1分层出高/中/低危三组，优于WHO亚型或单纯Ki-67。
MagicTPiT classifier development
基于参考系列代表β值签名构建R Shiny应用MagicTPiT，输入idat即输出亚群1-4及垂体对照打分；以0.72为匹配阈值，参考+验证190例甲基化样本准确率98.5%、覆盖率98.5%，对照垂体均打1.0，与Heidelberg联用增效。
讨论与结论翻译
研究人员指出四分子亚群具明确表观-转录-蛋白特征：亚群1大体积静默/功能混合、拷贝数不稳、轨迹差；亚群2全女静默、SSTR1+；亚群3对应既往GATA3+gonadotroph签名静默群；亚群4为USP8/USP48变异、年轻女性、微瘤显性CD，即既往三群框架因亚群2欠采样而漏判。WHO三亚型（SGCT/DGCT/CCT）与分子亚群不完全对齐，CCT不独成分子群，颗粒模糊比例高，提示组织学分型局限；SSTR1/GATA3/SSTR5免疫组化可作日常 surrogate。整合模型纳入侵袭、GTR与亚群1优于现行系统，Ki-67≥3%未显独立预后，或需按谱系重校。研究受meta合并异质性与回顾性限制，真实人群占比不可外推，前瞻验证待行。综上，该研究在《Nature Communications》发表，确立促肾上腺皮质激素型PitNETs四级分子分类、提供MagicTPiT开源分类器，将分子分层推入可转化病理实践。
