七氟烷暴露后老龄大鼠海马与人前额叶皮层细胞中铁死亡与神经炎症转录特征有限的整合再分析

《Frontiers in Cell and Developmental Biology》：Integrated re-analysis of bulk and single-cell transcriptomic data reveals limited ferroptosis and neuroinflammation signatures after sevoflurane exposure in the aged rat hippocampus and human prefrontal cortex cells

　　

背景： 实验研究已提示铁死亡（ferroptosis）与七氟烷（sevoflurane）相关的神经损伤有关，但这种信号在转录组系统、物种及神经细胞类型间的一致性仍不明确。因此，研究人员重新评估了公开可用的转录组数据集是否支持七氟烷暴露后存在一致的铁死亡与神经炎症（neuroinflammation）特征。 方法： 研究人员重新分析了两个公共转录组数据集：GSE141242，包含暴露于2.5%七氟烷4小时（n=3）或氧气对照（n=3）的年老大鼠海马组织；以及GSE196239，包含暴露于七氟烷、丙泊酚（propofol）、恢复条件或未处理对照的人胚胎前额叶皮层混合细胞培养物（每个条件一个文库）。质量控制后保留了67,272个单细胞。利用经验贝叶斯线性模型评估了整体差异表达，而针对无重复的单细胞条件，则通过关键基因表达、预定义基因集评分及观察到的细胞类型组成进行描述性总结。此外，研究人员在C6胶质瘤细胞和人星形胶质细胞（HA astrocytes）中进行了定量逆转录聚合酶链反应（qRT-PCR），将细胞分为对照组、3.3%七氟烷6小时组、H2O2（200 μM, 24小时）组或七氟烷+H2O2组，比较了四种条件下铁死亡相关基因的调控情况。 结果： 在年老大鼠海马中，无探针集达到错误发现率（FDR）<0.05（最小FDR=0.071），且两组在主成分空间中未明显分离。在人前额叶皮层培养物中，选定的铁死亡相关基因在单个对照文库与七氟烷文库之间仅显示出微小的描述性差异，预定义的基因集总结也相似。观察到的细胞类型比例差异小于两个百分点，且无广泛细胞类型表现出一致的七氟烷相关铁死亡或炎症模式。由于数据集在生物学和重复性方面存在显著差异，研究人员对其进行了定性比较而非正式的跨物种相关性分析。在四组qRT-PCR实验中，七氟烷和H2O2均降低了C6细胞和HA细胞中的谷胱甘肽过氧化物酶4（GPX4）和溶质载体家族7成员11（SLC7A11）的表达，并增加了酰基辅酶A合成酶长链家族成员4（ACSL4）和前列腺素内过氧化物合酶2（PTGS2）的表达；联合组的变化最大，所有处理组均与对照组存在差异（Dunnett校正后P<0.05）。由于单细胞数据集中每个条件仅有一个文库，所有条件间的单细胞比较均为描述性。 结论： 现有公开转录组数据并未提供强有力的证据表明在所检测的系统中存在由七氟烷引起的协调性铁死亡或神经炎症程序。四组实验显示，七氟烷和H2O2各自调节了铁死亡相关基因，其中联合组变化最大，这提供了两种暴露条件下通路扰动的补充证据。尽管单一药物组明确了每种暴露相关变化的趋势，但仅凭qRT-PCR无法确定铁死亡性细胞死亡或药理学相互作用。仍需开展具有更大生物学样本量和正交铁死亡检测方法的重复性研究，以确定铁死亡是否为七氟烷相关神经毒性的可重复组成部分。

研究背景与科学问题

铁死亡（ferroptosis）是一种铁依赖性的调节性细胞死亡形式，其特征为磷脂过氧化及细胞抗氧化防御系统的失效。GPX4-谷胱甘肽轴是其关键的抑制机制，而ACSL4依赖性多不饱和脂肪酸掺入膜磷脂则会增加铁死亡易感性。铁死亡细胞可能释放损伤相关分子模式（如HMGB1），从而在某些实验背景下将脂质过氧化驱动的细胞损伤与炎症信号联系起来。这些发现推动了铁死亡在神经损伤中的研究，但其存在与否及程度高度依赖于细胞类型、模型和实验条件。

七氟烷（sevoflurane）是广泛使用的吸入麻醉剂，实验工作引发了对长期或重复暴露可能改变神经元和胶质细胞稳态的担忧，尤其是在发育期或衰老神经系统。近期机制研究报告了七氟烷暴露后的铁死亡相关变化，包括H₂O₂积累、GPX4/SLC7A11信号、PRKCD-Hippo信号及PLIN4相关通路的改变。这些研究支持其生物学合理性，但使用了不同的物种、发育阶段、暴露强度、细胞系统和铁死亡检测指标。因此，在独立转录组数据集中能否检测到可重复的铁死亡转录特征仍不确定。

两个公共数据集为此问题提供了互补机会。GSE141242对七氟烷麻醉后老龄大鼠海马基因表达进行了分析，最初报告了与线粒体、代谢、衰老及神经退行性过程相关的广泛变化。GSE196239通过单细胞RNA测序对人胚胎前额叶皮层混合细胞培养物在七氟烷或丙泊酚暴露后的转录组进行了分析，报告了短暂、细胞类型依赖的转录反应。近期一项对公共整体和单细胞数据集的二次分析也提出了七氟烷神经毒性中SLC7A11相关的铁死亡机制。这些不同的解读使得透明化再分析、明确的表达量尺度控制以及对生物学重复与单个细胞的仔细区分尤为重要。

本研究聚焦于一个保守性问题：这些公共数据集是否显示出七氟烷暴露后与铁死亡和神经炎症一致且可重复的转录模式。研究人员将单细胞比较限制在描述性关键基因表达、预定义基因集总结和观察到的细胞类型组成；两个数据集被并列评估而非作为正式跨物种验证。细胞实验提供了互补的直接暴露背景，在C6胶质瘤细胞和人星形胶质细胞（HA astrocytes）中比较了对照组、3.3%七氟烷6小时组、H₂O₂（200 μM, 24小时）组及联合组，以评估选定铁死亡相关基因变化的方向和幅度。

主要技术方法

1. 1.
   公共转录组数据来源与分析： 重新分析了两个公共数据集。GSE141242使用Affymetrix Rat Transcriptome Array 1.0平台，包含18月龄雄性Sprague-Dawley大鼠海马组织，分别暴露于2.5%七氟烷4小时（n=3）或100%氧气4小时（n=3）。GSE196239使用10x Genomics Chromium v3平台，包含人胚胎前额叶皮层原代混合细胞培养物，分别暴露于3.3%七氟烷6小时、20 μM丙泊酚6小时及相应2小时恢复条件，外加未处理对照（每个条件一个文库）。
2. 2.
   整体微阵列数据处理与差异表达分析： 对GSE141242的RMA归一化log2尺度表达数据，采用经验贝叶斯线性模型（limma包）进行分析，并进行Benjamini-Hochberg多重检验校正，显著性定义为FDR<0.05。
3. 3.
   单细胞RNA测序处理、聚类与注释： 对GSE196239原始计数矩阵进行质量控制，保留67,272个细胞和23,597个基因。数据经归一化、对数转换、主成分分析、Harmony批次校正、UMAP嵌入和Leiden聚类后，通过基于标记基因的评分对九个宽泛细胞状态进行注释。
4. 4.
   描述性单细胞总结与体外验证： 由于每个条件仅有一个文库，未计算条件水平差异表达P值，而是按文库和细胞类型汇总关键铁死亡及炎症基因的平均表达及预定义基因集评分。在C6大鼠胶质瘤细胞和人星形胶质细胞中进行定量逆转录聚合酶链反应（qRT-PCR）验证，检测GPX4、SLC7A11、ACSL4和PTGS2的表达，采用两因素方差分析和Dunnett多重比较检验。

研究结果

单细胞质量控制： 五个单细胞文库经质量控制后，从67,330个输入细胞保留67,272个细胞，各文库保留率为99.9%-100.0%。检测到的基因数、文库大小和线粒体分数在各文库间相当。

前额叶皮层培养物的单细胞景观与注释： 标记基因引导的聚类解析出九种宽泛细胞状态，包括星形胶质细胞（32,819个）、少突胶质细胞（11,310个）、中间神经元（6,396个）、少突胶质前体细胞（5,471个）、增殖细胞（5,265个）、未成熟神经元（3,252个）、小胶质细胞（1,661个）、应激反应胶质细胞（951个）和内皮样细胞（147个）。来源文库在UMAP空间中相互混合。

七氟烷暴露后老龄大鼠海马的转录组反应： 六个样本的主成分分析未显示治疗驱动的清晰分离。探针集效应大小集中在零附近。经验贝叶斯分析在65,535个测试探针集中未识别出任何FDR<0.05的探针集，最小FDR为0.071。

描述性单细胞铁死亡与炎症总结： 在单个对照和七氟烷6小时文库中，选定铁死亡相关基因的平均log1p表达差异很小。例如，GPX4分别为2.034和2.061，SLC7A11为0.273和0.264，ACSL4为0.084和0.083，PTGS2均<0.001，HMGB1为2.001和1.990。铁死亡基因集对比和炎症基因集平均值在各文库间相似，未见七氟烷特异性升高。

细胞类型特异性描述性铁死亡总结： 对照和七氟烷6小时文库中，各宽泛细胞类型的铁死亡相关基因平均值大致相似。驱动基因和抑制基因平均值在不同宽泛细胞类型间的差异大于在两个观察文库之间的差异。七氟烷减去对照的铁死亡对比差值在各宽泛细胞类型间范围为-0.089至-0.001。

细胞类型特异性描述性炎症总结： 对照和七氟烷6小时文库中，各宽泛细胞类型的炎症相关基因平均值大致相似，观察到的差异小而方向混杂。小胶质细胞炎症基因集平均值最高，但在七氟烷6小时文库中并未升高。低表达炎症基因的检测频率差异通常低于一个百分点。

细胞类型组成与表达状态动态： 各文库间观察到的细胞类型组成变化不大，对照与七氟烷6小时比较显示小的谱系特异性差异。七氟烷6小时相对于对照的最大变化仍低于两个百分点。

整合复制与选定基因审计： 整体样本间Pearson相关系数为0.981-0.995，无治疗驱动的分组。选定基因的整体效应置信区间跨越零，经FDR校正后均不显著。单细胞文库中相同选定基因的平均表达在对照和七氟烷6小时之间相似。低表达炎症基因在两个文库中均少于1%的细胞中被检测到。

七氟烷和H₂O₂暴露后铁死亡相关基因失调的体外验证： 在C6胶质瘤细胞和人星形胶质细胞中，所有三个处理组（Sevo、H₂O₂、Sevo+H₂O₂）的GPX4、SLC7A11、ACSL4和PTGS2表达均与对照组存在显著差异（Dunnett校正后P<0.05）。保护性标志物GPX4/SLC7A11下调，促铁死亡标志物ACSL4/PTGS2上调，联合组变化最大。

讨论与结论

本再分析得出比七氟烷统一激活铁死亡及下游神经炎症模型更为谨慎的解读。在年老大鼠海马整体数据集中，无探针集达到FDR<0.05。在人前额叶皮层单细胞数据集中，单文库比较仅为描述性；经典铁死亡基因仅显示微小差异，预定义基因集总结相似，观察到的组成变化低于两个百分点。总体而言，公共转录组证据有限且非确认性，不支持稳健的跨系统铁死亡或神经炎症特征。

这一阴性或弱铁死亡信号并不意味着七氟烷暴露后不可能发生铁死亡。原始研究及直接实验研究已报告了七氟烷相关的铁死亡，差异可能源于暴露强度和时间、发育阶段、神经区室、实验平台、统计阈值及定义铁死亡的特定分子标准。本研究结果缩小了主张范围：在这些特定公共数据集中，在本研究应用的分析框架下，未能独立重现强烈的协调性铁死亡转录程序。

两个转录组系统不应被视为正式跨物种验证。它们同时在物种、脑区、发育状态、组织结构及暴露背景上存在差异。单细胞统计设计存在额外限制：GSE196239每个实验条件只有一个文库，单个细胞不能被视为条件水平推断的独立生物学重复。因此，本研究未使用细胞水平调整P值、细胞类型差异表达基因计数或富集结果来支持条件水平主张。

体外实验采用了四组暴露设计。两种暴露均以相同方向改变了铁死亡相关基因表达，联合组幅度最大。然而，仅凭基因表达变化不能确立铁死亡性细胞死亡或药理学相互作用。qRT-PCR结果为七氟烷和H₂O₂暴露后铁死亡相关通路扰动提供了补充证据。

细胞身份仍是重要的描述性变异来源。铁死亡和炎症基因集值在不同宽泛细胞类型间的差异大于在观察到的对照与七氟烷文库之间的差异。然而，由于条件比较缺乏生物学重复文库，当前数据未能识别出独特的铁死亡易感群体。

研究结论翻译： 本整合再分析表明，七氟烷暴露在所检测系统中与转录变化相关，但可用公共数据集并未显示出稳健且可重复的铁死亡或神经炎症特征。年老大鼠海马数据集未显示FDR显著的探针集，而人单细胞数据集在每个条件仅有一个文库的情况下，在经典基因、预定义基因集总结及观察到的组成方面仅显示出微小描述性变化。相比之下，七氟烷和H₂O₂均降低了C6胶质瘤细胞和人星形胶质细胞中GPX4/SLC7A11的表达并增加了ACSL4/PTGS2的表达，联合组变化最大，这提供了两种暴露条件下铁死亡相关通路扰动的补充证据。尽管定向基因模式与铁死亡相关调控兼容，但这些发现不应被解释为铁死亡性细胞死亡或协同作用的证据。总体而言，结果支持一种谨慎的、产生假说的解读，并强调了开展具有更大生物学样本量、样本水平统计学方法和正交铁死亡检测方法的重复性直接暴露研究的必要性。