泛素化(Ubiquitination)锚定型特征界定神经内分泌前列腺癌(Neuroendocrine Prostate Cancer, NEPC)：基于公共转录组数据集整合生物信息学分析鉴定枢纽基因(Hub Genes)及单细胞生态系统(Single-cell Ecosystem)洞察

《The Aging Male》：Ubiquitination-anchored signature defines neuroendocrine prostate cancer: hub genes and single-cell ecosystem insights from integrated bioinformatics analysis of public transcriptomic datasets

　　

背景：神经内分泌前列腺癌(Neuroendocrine Prostate Cancer, NEPC)是一种侵袭性强、治疗难治性的表型，常在雄激素受体(Androgen Receptor, AR)通路抑制下出现。研究人员假设泛素化(Ubiquitination)失调是NEPC谱系可塑性(Lineage Plasticity)的基础，整合批量(Bulk)与单细胞(Single-cell)转录组可定义泛素化中心特征及具有诊断与治疗价值的肿瘤微环境(Tumor Microenvironment, TME)回路。方法：研究人员合并公共批量转录组数据集比较NEPC与前列腺腺癌，包括来自GSE32967和GSE104786的发现队列共49例组织样本，并与精选泛素化相关基因全集取交集获得泛素化相关差异表达基因(Ubiquitination-Related Differentially Expressed Genes, URDEGs)；采用共表达网络、功能富集及蛋白质-蛋白质相互作用(Protein–Protein Interaction, PPI)网络拓扑分析鉴定并优先排序候选枢纽基因(Hub Genes)；整合独立scRNA-seq队列进行生物学情境化以映射细胞谱系、枢纽基因定位及细胞间通讯。结果：研究人员鉴定出317个URDEGs可将NEPC与常规前列腺癌明确分离，并聚类为富集细胞周期和有丝分裂程序的NEPC关联模块；网络整合得到11基因枢纽面板，包括AURKA、CCNA2、EZH2、FGFR1及TTK；单细胞数据质控后保留25,325个细胞，UMAP解析为17个簇并注释为8大谱系，单细胞映射进一步揭示谱系偏向性枢纽基因表达及显著的基质-血管以及免疫MIF–CD74/CXCR4信号轴。结论：该整合分析鉴定了公共批量数据集中与NEPC一致相关的泛素化锚定转录组特征及关联枢纽基因网络，并经独立scRNA-seq队列细胞类型解析模式支持；这些发现提名了候选生物标志物与治疗假说，需在前瞻性收集、临床注释队列中进行蛋白水平评估与外部验证。

本文对发表于《The Aging Male》的一项整合生物信息学研究进行解读。该研究针对神经内分泌前列腺癌(Neuroendocrine Prostate Cancer, NEPC)——一种在雄激素剥夺疗法(Androgen Deprivation Therapy, ADT)或雄激素受体(Androgen Receptor, AR)通路抑制剂治疗压力下出现的去势抵抗性前列腺癌(CRPC)罕见但极具侵袭性的亚型——探讨其分子发病机制。现有研究多聚焦个别泛素化(Ubiquitination)通路分子，尚缺乏系统性、多尺度刻画泛素化相关程序区分NEPC与常规前列腺腺癌的全景分析，且NEPC肿瘤微环境(Tumor Microenvironment, TME)中特定细胞群及信号网络仍不清楚。为此，研究人员假设泛素化失调驱动NEPC谱系可塑性(Lineage Plasticity)，通过整合公共批量转录组、加权基因共表达网络分析(Weighted Gene Co-expression Network Analysis, WGCNA)及单细胞RNA测序(Single-cell RNA Sequencing, scRNA-seq)数据，旨在鉴定泛素化相关差异表达基因(Ubiquitination-Related Differentially Expressed Genes, URDEGs)、构建共表达模块、识别枢纽基因(Hub Genes)及其调控网络，并在单细胞分辨率下解析NEPC TME的细胞组成与通讯网络。研究发现鉴定出317个URDEGs及一个11基因枢纽面板（含AURKA、CCNA2、EZH2、FGFR1、TTK等），单细胞图谱揭示8大细胞谱系及MIF–CD74/CXCR₄免疫轴，为NEPC早期诊断生物标志物发掘及靶向治疗策略提供了新线索。

研究人员采用的主要关键技术方法如下：①数据来源与处理——从GEO数据库获取批量转录组发现队列GSE32967（8例NEPC/SCPC与8例前列腺腺癌）和GSE104786（16例原发小细胞癌SCPC与17例高级别腺癌），合并后用sva/ComBat去除批次效应；泛素化相关基因(Ubiquitination-Related Genes, URGs)集通过GeneCards检索（相关性评分＞3）结合文献补充后去重。②差异表达与WGCNA——limma包进行差异分析（|log₂FC|＞0.5，FDR＜0.05），取URGs交集得URDEGs；选取前30%变异基因做WGCNA（软阈值β=6，动态剪切最小模块大小=50），保留|模块-性状相关系数|＞0.3的NEPC关联模块并与URDEGs相交。③功能富集与PPI——clusterProfiler做GO及KEGG富集；STRING建PPI网络，cytoHubba按最大团中心性(Maximal Clique Centrality, MCC)排名取前11基因为Hub Genes。④scRNA-seq分析——独立队列GSE137829经Seurat(v5.0)标准QC、归一化、PCA、UMAP及图聚类（分辨率0.4），SingleR加经典标记注释8大谱系；CellChat推断细胞间通讯(Ligand-Receptor, LR pair)。

Cohort assembly and harmonization（队列组装与调和）

研究人员将GSE32967与GSE104786经平台特异性探针重注释、合并及ComBat批次校正后构建Combined GEO发现队列（共49例，24例NEPC/SCPC，25例腺癌），PCA显示批次有效移除且保留生物学分组；独立scRNA-seq队列GSE137829用于生物学验证而非Hub基因提名。

Ubiquitin-related differentially expressed genes（泛素相关差异表达基因）

在NEPC vs 前列腺腺癌差异分析中获1,316个差异表达基因(Differentially Expressed Genes, DEGs)，与2,718个URGs取交集得到317个URDEGs；热图显示URDEGs可清晰区分NEPC与前列腺腺癌样本。

Weighted gene Co-expression network analysis（加权基因共表达网络分析）

WGCNA软阈值β=6达到无标度拓扑拟合R²≈0.90，生成5个模块（MEgrey、MEbrown、MEblue、MEturquoise、MEyellow），其中MEyellow、MEbrown、MEblue、MEturquoise与NEPC表型|r|＞0.30被保留为疾病关联模块；四模块基因与317个URDEGs相交获得PC–NEPC模块基因集。

Functional enrichment of PC–NEPC module genes（PC–NEPC模块基因的功能富集）

对317个URDEGs行GO（BP/CC/MF）及KEGG富集，显著富集于有丝分裂细胞周期时相转换、染色体分离、纺锤体组装等生物学过程及细胞周期(KEGG: hsa04110)通路，分子功能富集泛素样蛋白转移酶活性。

PPI network and regulatory architecture of hub genes（枢纽基因的PPI网络及调控架构）

STRING构建模块基因PPI网络，cytoHubba MCC算法筛选前11基因为枢纽基因：AURKA、CCNA2、EZH2、FGFR1、TTK、CCNB1、CDC20、BUB1、CDKN1A、UBE2C、PLK1；GeneMANIA显示其与20个功能相似蛋白共关联；ChIPBase和StarBase分别预测到25个转录因子(Transcription Factors, TFs)及50个microRNAs(miRNAs)调控上述Hub Genes。

Single-cell quality control and cell-type mapping（单细胞质量控制与细胞类型映射）

GSE137829经QC过滤后保留25,325个细胞，UMAP聚为17簇并注释为8大谱系——T细胞、B细胞、髓系细胞、肥大细胞、成纤维细胞、肌成纤维细胞、内皮细胞、上皮细胞；Hub基因呈谱系偏向性表达：AURKA富集于上皮细胞/肌成纤维细胞/肥大细胞，CCNA2于上皮及髓系细胞，EZH2广泛表达于内皮/上皮/B/髓系/肥大/成纤维细胞，FGFR1于肥大/成纤维-肌成纤维/内皮/上皮细胞，TTK主要在上皮细胞。

Single-cell differential expression and AUCell activity scores（单细胞差异表达及AUCell活性评分）

FindAllMarkers鉴定各簇标记基因与注释吻合；AUCell对各簇Top10标记基因集打分显示肥大细胞、B细胞及髓系细胞签名活性最高。

GO and KEGG enrichment of lineage-specific single-cell DEGs（谱系特异性单细胞DEGs的GO与KEGG富集）

分别对肥大细胞、B细胞、髓系细胞Top10 scDEGs富集分析：肥大细胞富集胞外基质相关、肾素-血管紧张素系统及5-羟色胺摄取调节；B细胞富集B细胞受体(B Cell Receptor, BCR)信号及免疫球蛋白复合体；髓系细胞富集补体激活及吞噬突触修剪。

Single-cell cell communication landscape（单细胞细胞通讯全景）

CellChat推断八谱系间配体-受体互作，全局互作频次与强度突出成纤维细胞-内皮细胞及内皮细胞内部通讯（基质-血管轴）；免疫区室中鉴定显著MIF–(CD74/CXCR₄)回路，最强互作见于T细胞→肥大细胞及肥大细胞→B细胞。

讨论(DISCUSSION)部分总结

研究人员指出317个URDEGs富集细胞周期和有丝分裂程序符合NEPC高增殖及ADT诱导可塑性特征，泛素化修饰可能是NEPC状态的关键门控环节。WGCNA捕获四个NEPC关联共表达模块，交叠URDEGs得精炼基因集；11基因枢纽面板涵盖有丝分裂调控(AURKA、TTK、CCNA2、CCNB1、CDC20、BUB1、PLK1)、染色质重塑(EZH2)、受体酪氨酸激酶信号(FGFR1)、细胞周期抑制(CDKN1A)及泛素偶联酶(UBE2C)，多层次TF/miRNA调控网络支持其汇聚节点地位。单细胞图谱确认8谱系并揭示Hub基因在上皮细胞主导表达及跨基质免疫细胞散布，提示NEPC表型受许可性微环境(Niche)支持；CellChat突显基质-血管轴及MIF–CD74/CXCR₄免疫轴可作为肿瘤外潜在可成药靶点。该泛素化锚定特征旨在补充而非替代传统神经内分泌标记物（如CHGA、SYP），提供正交分子信息。局限性包括纯回顾性公共数据二次分析未做湿实验验证、批次效应残留可能、临床注释不全限制亚组分析、scRNA-seq样本量有限等，需前瞻性标注队列及蛋白水平验证。

结论(Conclusions)部分翻译

本研整合批量转录组与scRNA-seq数据定义了公共数据集中区分NEPC/SCPC与前列腺腺癌的泛素化中心分子特征。研究人员提名了一个含AURKA、EZH2及FGFR1在内的11基因枢纽面板作为候选生物标志物，并映射出TME中显著的MIF–CD74/CXCR₄免疫通讯轴。这些发现为NEPC分子基础提供了新线索；然而所有结论源自计算机(In Silico)分析，枢纽基因及通路的诊断效用与治疗潜力需在独立前瞻性队列中通过免疫组织化学及功能实验进行严格验证。