完整糖蛋白组学揭示肾肿瘤 N-糖基化多层级异质性：基于 Shannon 多样性的位点特异性与功能域关联分析

《Molecular & Cellular Proteomics》：Deep Intact Glycoproteomic Analysis Reveals Multi-Level Glycosylation Heterogeneity in Renal Tumors

　　

在多种恶性肿瘤中可观察到糖基化的广泛变化，但关于被修饰蛋白及其位点的信息仍然有限。研究人员采用完整糖蛋白组学（intact glycoproteomics）对32例手术切除的肾肿瘤及配对的癌旁正常组织中的N-糖基化位点（N-glycosites）进行定位。研究人员从1,940种糖蛋白中鉴定出853,340条总N-连接糖肽（N-linked glycopeptides）和54,354条独特N-连接糖肽，涵盖220种聚糖组成（glycan compositions）。研究人员使用Shannon多样性指数（Shannon diversity index）表征糖基化多样性，并发现肿瘤相关差异富集于映射至细胞黏附（cell adhesion）和免疫（immunity）过程的糖蛋白。尽管在肿瘤中鉴定到更多糖肽，但全糖蛋白组（whole-glycoproteome）的Shannon多样性在组织类型间无显著差异；肿瘤相关差异集中于特定糖蛋白和糖基化位点。在肿瘤中，复合岩藻糖基化（complex fucosylated）和复合唾液酸化（complex sialylated）聚糖更为常见。基于RNA测序（RNA-seq）和uniLASSO的整合分析提名了与位点特异性聚糖类别变化相关的糖基化相关转录本。总蛋白丰度变化与蛋白水平（R² = 0.290）和糖基化位点水平（R² = 0.279）的Shannon多样性变化呈中等相关，但大多数线性变异不能仅由蛋白丰度解释。多样性变化还集中于功能和结构蛋白结构域。上述结果为肾肿瘤中蛋白糖基化变化的独立验证和机制研究提供了优先候选及假设。

研究背景与目的：蛋白质糖基化是丰度高且多样的翻译后修饰，参与蛋白折叠、定位、稳定性、活性及功能相互作用；异常糖基化如唾液酸化、岩藻糖基化增加与多种疾病包括癌症相关。肾细胞癌（RCC）是泌尿系统常见恶性肿瘤，已有研究多聚焦特定蛋白上的特定聚糖变化及其功能后果，但肾肿瘤与癌旁正常组织间糖基化多样性是否显著不同尚不清楚。高通量糖蛋白组学产生的糖基多样性变化难以捕捉和分析，因此研究人员采用Shannon多样性指数（Shannon diversity index）评估全糖蛋白组、单个糖蛋白和特定位点层面的糖基化多样性，以识别肿瘤相关多样性变化集中的蛋白、糖基位点、结构域和通路，并整合总蛋白丰度与RNA测序（RNA-seq）提名候选转录本-糖基位点关联。研究意义在于为肾肿瘤蛋白糖基化变化的独立验证和机制研究提供优先候选与假设。

关键技术方法（≤250字）：样本来自接受根治性肾切除术的患者（N=32），获取32对肾肿瘤与配对癌旁正常冰冻组织；最终病理包括透明细胞肾细胞癌24例、乳头状肾细胞癌1例、嫌色肾细胞癌1例、嗜酸细胞瘤4例、移行细胞癌1例和血管平滑肌脂肪瘤1例。研究采用完整糖肽液相色谱-串联质谱（LC-MS/MS）分析，经强阴离子交换-静电排斥亲水相互作用色谱（SAX-ERLIC）富集糖肽和高pH反相分级，Byonic进行N-糖肽搜索；同时进行总蛋白LC-MS/MS定量和bulk RNA-seq；以Shannon多样性指数计算全糖蛋白组、糖蛋白和糖基位点层级异质性，并用Mann-Whitney U检验及Benjamini-Hochberg校正控制错误发现率（false discovery rate, FDR）；采用uniLASSO（惩罚回归特征选择方法）整合糖基化相关转录本与糖基位点-聚糖类别特征。

研究结果：
1. 肾肿瘤组织单个糖蛋白显示显著高于癌旁组织的糖基化异质性：通过完整糖肽LC-MS/MS和Shannon多样性指数比较全糖蛋白组、单个糖蛋白和糖基位点层级；全糖蛋白组Shannon多样性在肿瘤与癌旁间无显著差异，但单个糖蛋白水平肿瘤显著更高（P<0.0001）。
2. 肾肿瘤中糖基化多样性增加的糖蛋白参与关键癌症相关生物学过程：对raw p<0.01的100个糖蛋白聚类，86个在肿瘤中多样性增加、14个降低；基因本体生物学过程（Gene Ontology Biological Process, GO-BP）过代表分析显示最显著类别为细胞外结构组织，并涉及免疫反应和细胞黏附。
3. 肾肿瘤中聚糖类别的描述性变化：三个最大糖蛋白簇约占显著蛋白63%；簇#3主要含高甘露糖，除高甘露糖下降外其他亚型增加；簇#6主要含复合唾液酸化+岩藻糖基化；簇#5主要含复合唾液酸化，并见复合岩藻糖基化、岩藻糖基化+唾液酸化及不完整核心增加、复合岩藻糖基化下降；原文强调这些类别变化为描述性。
4. 肾肿瘤与癌旁组织间常见蛋白糖基化差异：对177个在≥40%样本可计算多样性的蛋白进行Mann-Whitney U检验和Benjamini-Hochberg校正；FN1、MAN2B1显示一致且稳健的蛋白表达和糖基化异质性增加，IGHA2、AOC3、COL6A1、POSTN、A2M、CLU呈类似模式，CLU N374岩藻糖基化和唾液酸化增加；MMRN2、PLXNB2、TMED9多样性增加但总蛋白下降；CLPTM1L、UGT1A9、AHSG多样性指数与癌症分期显著相关（Pearson P<0.01）。
5. 蛋白丰度与糖基化多样性相关但未完全解释其变化：总蛋白丰度变化与Shannon多样性变化在蛋白水平（R²=0.290）和糖基位点水平（R²=0.279）呈中等正相关，约71.0%和72.1%的线性变异未由单变量线性模型解释；LAMP1 N249在癌旁更高、N261在肿瘤更高，但仅N249校正后显著，提示位点特异和多因素。
6. 肾肿瘤与癌旁组织间常见糖基化差异具有位点特异性：FN1的8个糖基位点在至少4例患者中检出，N430、N528、N542、N1007在肿瘤中多样性显著增加，其中N430、N528、N542位于胶原结合域，N1007位于中央区近III型纤连蛋白重复边界；CLPTM1L差异主要由N91驱动；UGT1A9主要由N71和N344驱动，复合岩藻糖基化增加并与分期相关；13个蛋白整体蛋白水平无显著差异（p>0.1）但至少一个位点显著变化，包括AZGP1、COL12A1、F2、HSPG2、HYOU1、ITGA1、ITGA3、ITGAV、LAMA5、ORM1、SERPINA3、SERPING1和VCAN。
7. 位点特异性糖基化差异与特定蛋白结构域相关：用PROSITE（蛋白结构域数据库）映射糖基位点；肿瘤中域内较域间糖基化多样性显著增加（p=0.008），正常组织有较弱趋势；域内，肿瘤中复杂交互域位点多样性高于域间，正常中简单交互域位点更高；肿瘤和正常中LAMP域位点多样性均高于域间。
8. 肾肿瘤中糖基化相关酶转录水平发生改变：bulk RNA-seq显示66个糖基转移酶或聚糖加工基因中22个在校正后p<0.01差异表达，除FUT6外均在肿瘤中升高；FUT6显著降低、FUT8显著升高，GFUS显著升高但非岩藻糖基转移酶；FUT11/POFUT4因O-岩藻糖基转移酶属性被排除；ST3GAL2而非ST3GAL1显著升高；原文指出未建立直接酶因果。
9. 整合建模识别候选酶-糖基位点关联：uniLASSO特征选择显示THBS1 N1067（hybrid）和FN1 N542（复合岩藻糖基化+hybrid）在超过5,000个模型中被选择，其他高频位点包括PSAP N332、ACAN N333、HSP90B1 N217、APOD N98、APMAP N160；RNA-seq模型高频基因包括MGAT1、DPAGT1、B3GNT8、ST3GAL3、ST3GAL5；Sankey图仅纳入>3%模型的糖基位点-聚糖类型和>1%模型的酶，定位为假设生成而非因果通路。

讨论部分总结：讨论强调肿瘤相关差异集中于选定蛋白和糖基位点而非全糖蛋白组；全糖蛋白组多样性相似，说明肿瘤变化并非简单反映全局糖基化保真度下降，Shannon多样性指数可捕获位点异质性。差异富集于细胞外基质（extracellular matrix, ECM）蛋白，与既往透明细胞RCC的ECM、细胞黏附和黏着斑富集一致；免疫相关糖蛋白多样性变化及肿瘤中复合唾液酸化增加与唾液酸化糖蛋白作为免疫细胞唾液酸结合免疫球蛋白型凝集素（Siglec）受体配体的免疫检查点作用相关。既往肾癌糖组/糖蛋白组研究观察到高甘露糖、唾液酸化及岩藻糖基化改变，但研究材料和分析层面不同；本研究聚焦蛋白和位点层面的聚糖类别分布及Shannon多样性。讨论指出位点多样性变化与复杂交互域、LAMP域及FN1胶原结合/折叠相关区域相关；蛋白丰度仅部分解释多样性，部分蛋白多样性升高而丰度降低，且同蛋白不同位点可呈相反方向。整合分析提名候选酶-糖基位点关联，但宽泛聚糖类别无法分辨连接方式、核心/外臂岩藻糖基化或异构体，转录升高不必然增加多样性。原文提出结构暴露、酶特异性、降解及功能糖型选择等可能机制，但强调数据未确定机制。局限性包括单队列、无独立生物学重复、DDA采样深度和蛋白丰度影响低丰度糖型检测、聚糖类别结构分辨率有限、RNA-seq关联不建立酶因果；未来需靶向糖肽定量、独立队列和样本再处理、候选酶扰动、定点突变及蛋白特异功能实验。

研究结论（翻译）：综上所述，对肾肿瘤的深度糖蛋白组学分析揭示，特定糖基化位点的聚糖多样性是肾肿瘤与正常组织之间最显著的差异。该多样性可通过Shannon多样性指数（Shannon’s diversity index）捕获，该指标有助于识别糖基化变化的关键区域。最大的变化高度位点特异性，聚集于功能和结构蛋白结构域，并集中于特定功能通路，最突出的是细胞外和免疫蛋白。整合RNA-seq建模识别了候选转录本-糖基位点关联以供未来检验。仍需独立验证和机制实验来确定所提名糖基化变化的功能后果和潜在生物标志物价值。