# 1. 标题 KPNA异构体特异性关联与胶质母细胞瘤恶性转录程序的联系：一项多队列转录组学分析

《Computational and Systems Oncology》：Isoform-Specific KPNA Associations With Malignant Transcriptional Programs in Glioblastoma: A Multi-Cohort Transcriptomic Analysis

## 3. 摘要译文
胶质母细胞瘤（GBM）表现出显著的转录异质性与可塑性，导致机制解读与治疗靶向十分困难。核输入通路是否与恶性转录程序存在系统性关联，以及在控制主要大体肿瘤混杂因素后，哪些核转运蛋白α（KPNA；importin-α）异构体呈现独立的特异性关联，目前尚不清楚。研究人员开展了一项多队列转录组学分析，以癌症基因组图谱（TCGA）GBM RNA-seq数据集为发现队列（n = 143，配以匹配的肿瘤纯度和拷贝数数据），以非TCGA的胶质瘤纵向分析（GLASS）成人弥漫性胶质瘤队列（n = 191）进行外部验证，并以临床蛋白质组肿瘤分析联盟（CPTAC；n = 99）进行支持性跨平台验证。恶性转录程序评分与KPNA表达之间的关联采用增殖校正的联合模型进行评估，并根据队列特异性的数据可用性额外纳入肿瘤纯度和基因水平拷贝数。KPNA6与KPNA4是优先级最一致的两个家族成员，并在相互关联的转录程序间表现出互补的、程序依赖性的关联。KPNA6优先与黏附与缺氧程序相关联，而KPNA4则与炎症程序（包括肿瘤坏死因子-α（TNFα）–核因子-κB（NF-κB）及干扰素-γ（IFN-γ）信号）相一致。这些模式在协变量校正后依然保持。外部分析具有补充性，但取决于具体队列：GLASS主要支持KPNA6相关程序，而CPTAC支持KPNA4相关程序。尽管拷贝数变异对KPNA表达变异性有贡献，但并不能解释所观察到的转录关联。这些发现确定了KPNA表达与GBM转录程序之间的异构体特异性关联，并为未来对KPNA6和KPNA4的功能研究奠定了基础。

## 4. 论文解读

### 研究背景
胶质母细胞瘤（GBM）即使接受多模式治疗仍为致死性原发脑肿瘤，其生物学特征以显著的转录异质性为主，这使机制推断与治疗靶向都变得复杂。大规模图谱研究确立了可重复的亚型样表达模式及以通路反复改变为标志的基因组图景，但同时表明大体肿瘤常混合多种程序而非分离为单一状态。单细胞与整合框架进一步显示GBM细胞占据一系列动态转录状态，可沿应激相关、侵袭、缺氧与免疫相关程序转变。间质/炎症信号网络及细胞外基质（ECM）重塑程序与侵袭、治疗耐药及免疫髓系区室互动紧密相连。核质运输是调控信号介质与转录因子何时、如何进入基因组的调控层之一。经典核输入依赖核转运蛋白α（KPNA；importin-α）衔接蛋白识别货物蛋白上的核定位信号，并与importin-β协同穿越核孔。7个人KPNA旁系同源基因（KPNA1–KPNA7）功能并非完全可互换，在序列上已分化，可表现出货物偏好及情境依赖性表达，从而介导上游信号与下游转录输出间的异构体特异性耦联。

目前存在的问题是：KPNA家族尚未在GBM中系统性地映射到恶性转录程序并同时处理大体肿瘤转录组固有的主要混杂因素。既有研究多聚焦单一importin-α基因（最突出的是KPNA2），将其与增殖相关表型及致瘤特性相联系；KPNA4亦被牵连于GBM相关过程（包括GBM模型中的间质转化）。然而，当对整个家族进行联合建模时，GBM状态变异究竟最好由更广泛的"核输入上调"现象所捕捉，还是仅反映一小部分KPNA异构体不成比例地追踪不同恶性程序，仍不清楚。此外，大体肿瘤表达受肿瘤细胞度与微环境混杂影响，免疫或炎症特征对此尤为敏感，因此需明确控制纯度及相关协变量；同时拷贝数改变在GBM中普遍存在并可能调节表达剂量，若不加以处理会伪装为调控耦联。为此，研究人员开展回顾性多队列转录组分析，以TCGA为发现队列并整合匹配的肿瘤纯度（ABSOLUTE）与基因水平拷贝数（GISTIC2）以校正混合与剂量效应，以ESTIMATE衍生的纯度/免疫混杂代理作敏感性分析，使用精选基因集与队列内标准化状态评分量化程序活性，在GLASS队列测试可重复性，并以CPTAC GBM RNA测序作为支持性跨平台验证；除家族范围模型外，采用联合建模估计相互校正的异构体特异性关联并优先排序候选基因。

### 关键技术方法
（1）多队列设计：TCGA-GBM为发现队列（n = 143，含完整纯度与CNV），非TCGA的GLASS成人弥漫性胶质瘤为主要外部验证队列（n = 191，一例一样本，优先初始时间点），CPTAC GBM RNA-seq（n = 99）为支持性跨平台验证；另分析公共GSE131928单细胞RNA测序数据集。数据来自UCSC Xena、PanCanAtlas、Synapse（syn69961520）及CPTAC数据门户等公开数据库。（2）程序评分：基于MSigDB v2025.1 Hallmark与精选（c2.cp）基因集构建EMT、缺氧、炎症、TNFα–NF-κB、IFN-γ反应、黏附及增殖评分，基因z标准化后取均值；剔除KPNA2定义proliferation_noKPNA2以避免循环性。（3）统计建模：OLS回归配HC3稳健标准误，逐步模型M0–M3（纯度、基因CNV、增殖），联合模型同时纳入KPNA1–7，增量R2量化贡献，1000次患者级bootstrap评估选择稳定性，BH控制FDR。（4）验证与分层：GLASS临床敏感性分析（年龄、性别、复发阶段、中心、MGMT）、KPNA6×KPNA4中位数四分位分层（Kruskal–Wallis及BH校正Mann–Whitney U，Cliff's delta效应量）及CNV–表达耦联Spearman分析。

### 研究结果
**TCGA-GBM中KPNA1–7的表达图景**：限制为原发肿瘤并每患者取一样本后，143例具匹配纯度与KPNA4/KPNA6拷贝数数据的肿瘤构成主要发现分析集。中位诊断年龄62岁（IQR 53–70），男性65.5%，中位总生存334天（IQR 147–488），107/142例事件（75.4%）。肿瘤纯度中位0.80（IQR 0.65–0.88）；两基因拷贝数状态以中性为主（KPNA4 75.5%、KPNA6 70.6%为CN = 0）。KPNA7显著低于其他成员（中位log2(TPM + 1) = 0.13），KPNA4与KPNA6跨越较宽动态范围（分别为5.03与5.08）。

**恶性状态评分间的相关结构及其与KPNA轴评分的关系**：状态评分彼此强相关，提示耦联的转录结构；EMT样评分与黏附、缺氧、TNFα–NF-κB强相关（ρ = 0.90、0.80、0.80），炎症反应与TNFα–NF-κB（ρ = 0.86）及IFN-γ反应（ρ = 0.77）强相关，而增殖与免疫/炎症反应弱相关（炎症ρ = ?0.16；IFN-γ ρ = ?0.04）。基因集重叠以EMT与黏附最大（81个共享基因；Jaccard指数0.184）。主成分分析识别出解释89.1%方差的两个成分（PC1 72.7%；PC2 16.4%）；校正联合模型中KPNA6与KPNA4均与PC1正相关（β = 0.359与0.292，均q < 0.001），PC2与KPNA6负相关（β = ?0.270，q < 0.001）而与KPNA4正相关（β = 0.205，q = 0.026）。轴评分筛查显示alpha-1与alpha-2评分与增殖相关性最强（ρ = 0.63与0.65），与黏附、缺氧呈较弱的正相关，与炎症及IFN-γ反应相关性可忽略。

**TCGA-GBM中KPNA4与KPNA6的家族范围与聚焦关联**：同时纳入KPNA1–7并校正纯度、基因水平CN及proliferation_noKPNA2的多变量模型解释了相当比例的方差（各程序R2范围0.362–0.536），KPNA6与KPNA4呈现最一致的程序特异性关联。KPNA6对黏附（β = 0.185，q = 9.54 × 10^?7）、EMT（β = 0.178，q = 5.32 × 10^?4）、缺氧（β = 0.161，q = 2.83 × 10^?3）及TNFα–NF-κB（β = 0.172，q = 2.33 × 10^?3）呈FDR显著正系数；KPNA4对TNFα–NF-κB（β = 0.213，q = 1.92 × 10^?4）、EMT（β = 0.185，q = 5.32 × 10^?4）、缺氧（β = 0.105，q = 3.49 × 10^?2）、IFN-γ反应（β = 0.210，q = 8.29 × 10^?3）及炎症反应（β = 0.134，q = 1.42 × 10^?2）呈FDR显著正关联。增量解释方差显示互补的程序依赖模式：KPNA6在黏附（ΔR2 = 0.117）与缺氧（0.083）最大，KPNA4在EMT（0.077）、炎症（0.046）、TNFα–NF-κB（0.097）及IFN-γ反应（0.078）居首。联合模型中KPNA6仍与黏附（β = 0.175，95% CI 0.124–0.227，q = 1.53 × 10^?10）与缺氧（β = 0.148，q = 6.41 × 10^?5）显著关联，并独立关联EMT、炎症反应与TNFα–NF-κB，但与IFN-γ反应在相互校正后无关联（q = 0.491）；KPNA4则保留与IFN-γ反应（β = 0.158，q = 8.74 × 10^?4）、炎症反应（β = 0.113，q = 2.97 × 10^?4）、TNFα–NF-κB（β = 0.145，q = 1.25 × 10^?4）、EMT（β = 0.0958，q = 4.20 × 10^?3）及缺氧（β = 0.0839，q = 2.37 × 10^?3）的关联，而与黏附无关（q = 0.255）。

**CNV–表达耦联及协变量校正关系支持TCGA-GBM中两种不同的KPNA6与KPNA4相关程序**：阈值化CN状态与两基因表达相关，KPNA6（ρ = 0.34，p = 4.0 × 10^?5）强于KPNA4（ρ = 0.17，p = 0.048），支持发现模型中显式校正CNV，并提示联合模型中的KPNA–程序关联不太可能仅由拷贝数剂量解释。协变量校正后的程序评分在四分位间呈现与联合模型一致的模式：较高KPNA6表达与升高的缺氧和黏附评分相关，较高KPNA4表达与炎症和IFN-γ相关评分增加相关。

**GLASS成人弥漫性胶质瘤队列中KPNA6与KPNA4关联的外部验证**：该队列（n = 191）包括154例初始时间点与37例复发样本；IDH野生型135例、突变型46例、不可用10例；IV级155例、III级15例、II级20例。全队列中KPNA6在增殖校正联合模型下与全部六个转录程序显著正相关，而KPNA4关联较弱。在128例具临床注释的IDH野生型IV级GBM敏感性分析中，KPNA6关联保留而KPNA4关联减弱；KPNA6与全部六个程序的关联在校正年龄与性别后、以及在102例无既往治疗记录的初始时间点病例中仍FDR显著；在128例队列中额外校正复发阶段与贡献中心后，黏附与缺氧仍显著（q = 0.010与0.035）；在中心校正的102例分析中黏附仍显著（q = 0.025）而缺氧减弱（q = 0.065）；在60例具MGMT甲基化状态者中黏附、缺氧与炎症仍显著；KPNA4在这些分析中无FDR显著关联。

**CPTAC与单细胞分析**：CPTAC GBM RNA-seq（n = 99）的增殖校正联合模型中KPNA4与全部六个程序评分正相关，而KPNA6无FDR显著关联；额外校正匹配的KPNA4与KPNA6拷贝数（n = 98）后，KPNA4六个关联中的五个仍显著，黏附关联减弱，KPNA6仍不显著；ESTIMATE衍生的免疫与基质成分敏感性分析支持TCGA联合模型关联的稳健性；程序间相关结构在CPTAC中大体保留。GSE131928单细胞分析显示KPNA4与KPNA6均在恶性与非恶性区室中表达；KPNA4在恶性细胞与任何匹配非恶性区室间无显著差异，KPNA6在恶性细胞中高于匹配巨噬细胞（n = 9对样本，q = 0.023），与少突胶质细胞和T细胞的比较不显著，两基因均未局限于恶性区室。

### 讨论与结论
讨论指出三项值得注意的发现：其一，KPNA6与KPNA4总体上是最一致被优先排序的家族成员，其选择稳定性互补且程序依赖（bootstrap显示两者在六个程序中的平均前二选择频率最高，分别为69.5%与71.8%）；其二，完全校正的相互校正联合模型中两者呈现部分不同的关联谱——KPNA6与黏附/ECM及缺氧程序关联特别强，KPNA4与免疫/炎症程序联系更强，并在EMT与TNFα/NF-κB信号上有共享关联；其三，关键关联在外部非TCGA队列GLASS中得到确认，支持可重复性并揭示优先捕捉IFN-γ反应变异的基因存在队列特异性差异。核输入生物学为此提供可能的组织框架：KPNA作为经典核输入衔接蛋白，异构体在功能上并非可互换，可呈现优先货物互作与情境依赖表达；家族范围模型提示有限亚群（以KPNA6与KPNA4最突出）在与其他成员联合评估时解释多个恶性程序的不成比例方差，支持（但不证明）特定importin-α异构体可能作为GBM转录环路中状态依赖节点的概念，而非统一"核输入上调"现象。KPNA6（importin α7）与ECM/黏附及缺氧程序相关，而缺氧与ECM重塑是GBM病理生理标志特征；KPNA4（importin α3）与免疫/炎症程序关联最强，与KPNA4在其他细胞情境中介导NF-κB与STAT3核输入以及NF-κB/STAT3网络与GBM间质转化及炎症信号的关联相呼应，但分析为相关性且未鉴定具体货物。两者关联模式呈梯度而非互斥，与共享的间质样轨迹上应激、侵袭与炎症程序共激活一致。讨论还分析TCGA与GLASS差异的可能原因（GLASS样本跨中心与时间点、包含复发、验证模型未纳入匹配纯度与拷贝数协变量、免疫/基质混杂对炎症与IFN特征影响大、单细胞分析显示两基因不局限于恶性区室等）。CN整合提供了基因组剂量的额外控制，显式校正后关联仍存在。探索性生存分析未显示二分类KPNA6或KPNA4表达能可靠分层总生存。局限性包括：分析均为大体转录组关联，不能确立因果或鉴定货物蛋白；大体表达仍为恶性与非恶性细胞群的复合；外部验证改变队列结构与协变量可得性；基因集评分依赖精选特征与队列内标准化；四分位分层仅提供组合结构可视化，正式组间检验（n = 142）显示黏附、缺氧与炎症存在显著四分位间差异而IFN-γ未达显著，缺氧与炎症在6-low/4-low与6-high/4-high间对比最大。研究优势包括多组学发现队列配匹配纯度与CNV协变量、系统比较多变量模型中全部KPNA家族成员以考虑基因间依赖、以增量解释方差量化相对重要性、采用相互校正的KPNA6 + KPNA4模型检验非冗余性并在外部队列验证。未来工作应超越相关性，开展单细胞与空间转录组学解析恶性状态特异性与免疫生态位，通过CRISPR介导的KPNA6与KPNA4敲低或过表达测试异构体丰度改变是否影响候选货物转录因子的核定位并转变下游程序，以及importin互作组蛋白质组学鉴定货物。

**结论译文**：在以TCGA为核心、以GLASS数据进行外部验证的多队列分析中，KPNA6与KPNA4被确定为GBM恶性转录程序最具信息量的KPNA家族相关因子。相互校正模型支持一种加权化的关联模式划分：KPNA6优先与黏附/ECM及缺氧程序相一致，而KPNA4优先与炎症信号状态相一致。这些发现确定了KPNA表达与GBM转录状态之间的异构体特异性关联，并推动功能研究以判定这些关联反映的是因果作用还是对肿瘤状态动态的适应性反应。