《Neurobiology of Disease》：Within-state transcriptional remodelling predominates in a 41-gene ageing programme of human prefrontal OPCs
编辑推荐：
（1）一个由41个基因组成的OPC（少突胶质前体细胞）衰老程序在两个外部皮质队列中重现。
（2）状态内变化（within-state change）主导其调整后的细胞平均（cell-mean）年龄对比。
（3）在配对数据中，一种承诺相关状态（commitmen
（1）一个由41个基因组成的OPC（少突胶质前体细胞）衰老程序在两个外部皮质队列中重现。
（2）状态内变化（within-state change）主导其调整后的细胞平均（cell-mean）年龄对比。
（3）在配对数据中，一种承诺相关状态（commitment-associated state）表现出更低的相对占据（relative occupancy）。
（4）在配对供体中，启动子可及性（promoter accessibility）与RNA程序共变。
少突胶质前体细胞（Oligodendrocyte precursor cells, OPCs）在成人人类大脑中持续存在，负责生成少突胶质细胞以进行髓鞘维持和修复。该谱系经历前体、承诺相关及早期髓鞘化状态，且OPCs随年龄增长获得区域与功能异质性。现有实验性衰老模型表明，通过细胞内在、表观遗传及环境机制可导致OPC成熟与再髓鞘化受损。然而，这些分子与细胞效应在人类皮质OPC群体内的组织方式仍不清楚。年龄相关的群体信号可能源于细胞状态比例变化、状态内基因表达变化，或两者兼有。群体平均分析无法区分这些可能性，而衰老图谱研究虽显示分子改变与细胞组成变化可共存，但尚不清楚广义OPC程序反映的是状态混合改变还是跨状态分布的分子变化。因此，确定固定转录衰老程序反映的是状态占据改变、状态内重塑抑或两者兼具，成为解析人类皮质OPC衰老机制的关键科学问题。
为解答上述问题，研究人员基于配对人类前额叶单细胞多组学数据定义了一个固定的41基因双向衰老程序，并在两个独立外部皮质数据集中验证其重现性。研究进一步在349名供体的配对RNA-ATAC（RNA-seq与ATAC-seq）谱中分解细胞平均年龄对比，同时考察承诺相关状态的相对占据及启动子可及性变化。研究结论表明，该41基因程序的调整年龄对比主要由状态内转录重塑主导，而承诺相关状态的相对占据下降构成一个独立的伴随特征。这一发现将生物学关注点从识别单一低丰度细胞状态，转移至解释分布于各OPC状态的分子变化，为理解人类大脑衰老过程中少胶质谱系的调控组织提供了新的理论框架，并提示未来研究需关注状态内分子重塑而非仅仅细胞比例改变。
研究人员首先基于Catching等提供的配对人类前额叶多ome RNA数据（228名供体，年龄20–70岁）及PsychAD衰老队列背外侧前额叶皮质RNA数据（73名供体，年龄24–70岁），通过堆叠供体聚集的CR1模型筛选并固定了包含22个年龄下调和19个年龄上调基因的41基因程序。随后在GSE254569（眶额叶皮质，87名供体）和GSE291605（30名前额叶供体）两个外部数据集中检验该程序的重现性。在Catching配对RNA-ATAC数据（349名供体）中，采用无监督加权K-means聚类（K=6）定义OPC转录状态，并通过层级分解框架将调整后的细胞平均对比拆分为组成、状态内变化及交互项。此外，运用beta-binomial（β-二项）模型与logistic回归分析承诺相关状态的相对占据变化，并基于转录起始位点（transcription start site, TSS）侧翼2 kb窗口评估启动子可及性与RNA程序的共变关系。
3.1 双向OPC衰老程序在两个外部皮质队列中重现
通过在GSE254569和GSE291605中对固定41基因模块进行协变量调整分析，研究人员发现年龄下调与年龄上调基因在眶额叶皮质（orbitofrontal cortex, OFC）和前额叶皮层中均保持预期方向，全部41个基因的关联方向均在外部数据集中得以保留。眶额叶皮质中，年龄下调与上调系数分别为每十年?0.3280和+0.2511，而另一数据集中为?0.3057和+0.2965。所有基因在两个资源中的效应方向一致，证实了固定程序的高度重现性。
3.2 衰老程序在注释OPC区室内持续存在
在作者注释的OPC核中，无论采用何种质控表示，广义衰老程序仍然显著。在Catching配对多ome资源中，去除承诺相关状态后，伪bulk衰老轴系数基本不变，表明该罕见状态对广义年龄相关信号贡献甚微，程序信号并非由单一细胞状态比例变化驱动。
3.3 状态内重塑主导固定程序对比
通过四项分解框架，研究人员将调整后的细胞平均年龄对比分解为组成、状态内变化、交互及重建余项。结果显示，在HBCC和NABEC供体来源中，状态内变化分别占直接对比的99.5%和98.8%，而组成项仅占0.8%和1.2%。该主导性在K=4至12的多种聚类参数及不同表示方法下均持续存在，且移除所选状态后广义伪bulk轴亦基本不变。
3.4 Catching中承诺相关状态显示更低相对占据
研究人员在两个Catching供体来源中均鉴定出一种承诺相关OPC状态，其相对占据随年龄增长而降低。该状态富含OPC-II型相对于OPC-I型的标志物，且存在/缺失模型与beta-binomial模型均支持显著的负向年龄关联，其对数优势系数在HBCC中为?0.300，在NABEC中为?0.337。这一占据变化在数量上可区别于广义转录程序，移除该状态后程序轴基本不变，提示二者为可分离的不同层面。
3.5 启动子可及性与转录衰老程序平行变化
通过基于GENCODE-v32注释的转录起始位点侧翼2 kb启动子窗口分析，研究人员发现下调与上调基因启动子的可及性变化方向与RNA程序一致。配对供体中RNA模块与启动子加性对数比评分呈显著正相关，在HBCC和NABEC中调整协变量后多数关联仍保持显著，提示启动子可及性变化构成了该转录程序的表观遗传调节相关物。
讨论部分强调，41基因程序在未贡献基因选择的外部皮质数据集中重现，其年龄对比在Catching配对数据中由状态内重塑主导，而承诺相关状态相对占据的降低则表现为独立特征。尽管状态组成对该特定评分对比贡献有限，但可能影响谱系输出与修复能力，因而将科学问题从识别单一低占据状态转向解释分布性分子改变。承诺相关占据下降提示可能存在细胞进入减少或转位改变，但由于相对占据不反映绝对丰度，需组织学验证。外部数据集中，眶额叶皮质因具备明确疾病与技术协变量而锚定效应量解释，而GSE291605因年龄-区域不平衡主要支持方向性重现。启动子可及性结果确立了RNA与染色质层面的一致性，但尚不能确立因果关系，需进一步扰动实验验证。此外，部分年龄下调基因较长且在神经元中表达，可能受环境RNA污染影响，尽管经质控校正与神经元代理调整后方向仍保留，但空间转录组学或单分子FISH（fluorescence in situ hybridization）仍是定位转录本所必需的。因此，确定上述变化如何影响少突胶质维持与修复是后续研究的核心任务。
该41基因衰老程序在分析的人类皮质数据集中重现。在Catching配对谱中，其年龄对比由状态内重塑主导，并伴随一个承诺相关状态的独立相对占据变化。确定这些变化如何影响少突胶质维持与修复是下一步的实验步骤。
