用共存网络解析拟南芥叶片EXO70亚基空间表达图谱

原文：The expression pattern of EXO70 subunits in single-cell stereo-seq of Arabidopsis thaliana leaves using coexistence network analysis
期刊：Scientific Reports | DOI: 10.1038/s41598-026-44501-9
作者：Xin Yu; Jingmin Shang
机构：School of Life Sciences, Qilu Normal University, Jinan 250200, China；School of Chemistry and Chemical Engineering, Qilu Normal University, Jinan 250200, China

1. 研究内容：解决低丰度基因家族在单细胞空间组学中的解析难题

在拟南芥（Arabidopsis thaliana）中，exocyst 复合体的 70 kDa 组分——EXO70 家族至少包含 23 个亚型（isoforms）。这些蛋白负责把囊泡定向运输到质膜，是胞吐与极性分泌过程中的关键节点。不同 EXO70 亚型在叶片不同细胞类型（表皮、叶肉、维管等）中可能承担差异化功能，但其在叶片组织内的表达异质性一直缺乏系统刻画。

问题在于：单细空间转录组（如 stereo-seq）虽然能在细胞分辨率下绘制基因表达图谱，但单个细胞中检测到的表达基因比例仅约 2.5%，属于高度稀疏数据。对于 EXO70 这类“低表达但功能关键”的基因家族，传统共表达网络（如 WGCNA、ARACNE 或 Spearman 相关矩阵）往往难以捕捉有意义的关联。因此，作者提出用“共存网络（coexistence network）”框架，把稀疏矩阵转为二元存在/缺失矩阵，从而分析低表达基因家族在复杂组织中的空间共现模式。本研究的通讯作者来自齐鲁师范学院（Qilu Normal University）生命科学学院。作者从 StomicsDB 获取拟南芥叶片单细胞 stereo-seq 数据（ accession ID STDS000010421），覆盖上表皮、下表皮、海绵组织叶肉细胞、栅栏组织叶肉细胞、维管细胞与保卫细胞等六类细胞。每个切片平均约 630 个细胞，检测基因约 18,310 个；但高表达基因（均值>1）平均仅约 105 个，不足 1%。

在 EXO70 相关基因中，表达水平整体很低：保卫细胞中仅 SEC3A 与 EXO70D2 有表达；除 EXO70G2 外，其余基因在表皮与叶肉细胞中均有表达；维管与上表皮细胞表达更低。

作者先用 STRING 数据库构建以 EXO70 家族为种子、扩展至 VETH 与 VAP 等囊泡相关蛋白的小型互作网络，用于后续分析。

关键结果包括：

• 在上表皮细胞中，传统 Spearman 共表达矩阵仅发现 EXO70H7 与 SEC3A 显著相关（p<0.05）；而共存矩阵能揭示更多弱但存在的关联。
• 通过置换检验（permutation test，10,000 次），在上表皮细胞中识别出两对显著共存基因：EXO70A1 与 EXO70H7、EXO70F1 与 EXO70B1（p<0.05）。
• 在海绵组织叶肉细胞中，EXO70A1、EXO70H7 与 SEC3A 表现出显著共现，提示 exocyst 亚基在特定叶肉细胞中的功能协同。
• SEC 家族基因在叶肉共存网络中富集，可能参与叶绿体发育或类囊体膜运输相关过程。

2. 方法与技术亮点

本研究的方法核心是“共存网络”流程：将稀疏表达矩阵二值化（非零为 1），通过矩阵转置与稀疏矩阵乘法计算共存矩阵，再除以“总存在数减去共存数”的分母矩阵，得到共存比率 P(CE)。随后用置换检验比较观测共存概率与随机期望 P(G_A|G_B)，以 p<0.05 判定显著共存。

相比依赖表达量大小的传统相关分析，共存分析只关注“是否在同一细胞中出现”，因此更适合单细胞空间组学中低检出率、低表达基因家族的互作挖掘。分析代码已开源在 GitHub（setereo_seq_gcoa）。

3. 研究影响与展望

该工作一方面揭示了 EXO70 亚型在拟南芥叶片不同细胞类型中的功能特化：例如上表皮中 EXO70H7 与 EXO70A1 的共现可能关联胁迫响应与细胞壁修饰；海绵叶肉中 EXO70A1、EXO70H7 与 SEC3A 的协同则指向囊泡运输的空间组织。另一方面，它验证了共存网络可作为分析低丰度基因家族（如转录因子、信号激酶）在复杂组织中互作的可扩展框架。

作者也指出局限：多数 EXO70 亚基受转录后调控，RNA 丰度不一定等于蛋白水平（如 EXO70H4 启动子广泛活跃但蛋白仅限定于毛状体）。因此未来需结合原生启动子驱动的荧光融合蛋白进行验证。总体而言，该研究为植物空间组学中“低表达但关键”的基因网络提供了方法论范例。

本文图片基于论文 PDF 生成，版权归原论文所有，仅供学习参考，请以原文为准。