脉冲式17β-雌二醇在女性与男性半月板纤维软骨细胞中诱导差异转录组反应

《International Journal of Molecular Sciences》：Biomarkers for Cancer Immunotherapy: Advances, Challenges, and Clinical Applications

　　

骨关节炎（osteoarthritis，OA）在女性中的发病率高于男性，雌激素被认为参与关节稳态的维持，然而其对于半月板纤维软骨细胞（meniscal fibrochondrocytes，MFCs）的影响仍不明确。研究人员对一项公共RNA测序数据集（Gene Expression Omnibus，GEO；GSE199087）进行了重新分析，该数据集包含来自两位年龄匹配供体（一名女性与一名男性，均为47岁）的MFCs。各供体的细胞以三次独立重复的方式处理，分别接受脉冲式17β-雌二醇（17β-estradiol，E2；240 μM，每24小时处理1小时，持续72小时）或溶媒对照。差异表达分析采用limma-voom方法，通路鉴定通过基因集富集分析（Gene Set Enrichment Analysis，GSEA）、Enrichr及iLINCS进行，仅保留受三种方法共同支持的基因本体（Gene Ontology，GO）条目。在女性供体的细胞中，E2上调了DNA复制与细胞周期通路，并降低了MMP1、MMP2、IL12A及MYD88的表达水平。在男性供体的细胞中，E2上调了染色质重塑与发育信号通路，其中NOTCH1为排名居前的前沿基因（leading-edge genes）之一，同时下调抗原呈递通路及线粒体相关基因。在基因水平层面，部分改变为两位供体所共有，包括TIMP1表达下调。鉴于每种性别仅由一名供体代表，上述差异无法确证归因于性别。相关发现具有假设生成价值，并为候选通路在大规模、性别均衡队列中的验证指明了方向。

骨关节炎（osteoarthritis，OA）是一种以关节软骨退变为特征的致残性关节疾病，其在女性中的患病率显著高于男性，尤其在绝经后女性群体中更为突出，这一流行病学特征提示性激素特别是雌激素在OA发病与进展中可能扮演重要角色。已有研究表明，雌激素能够影响组织修复与稳态，17β-雌二醇（17β-estradiol，E2）被认为可通过抑制基质金属蛋白酶（matrix metalloproteinases，MMPs）的表达发挥软骨保护作用，因为MMPs正是驱动细胞外基质（extracellular matrix，ECM）与胶原蛋白降解的主要酶类。然而，半月板在膝关节负荷传递与冲击吸收中发挥关键作用，其细胞中的雌激素作用机制却长期未被阐明。半月板纤维软骨细胞（meniscal fibrochondrocytes，MFCs）具有独特的可塑性，对半月板损伤修复至关重要，但E2对男女MFCs转录组的影响是否存在差异、男女细胞是否通过相同分子通路限制基质分解，均悬而未决。正是针对这一问题，研究人员对公共RNA测序数据集GSE199087进行了重新分析，比较了年龄匹配的一男一女两名47岁供体MFCs对E2处理的响应，并结合三个互补的通路分析平台，最终得出E2在两位供体细胞中激活了不同转录程序的结论：女性供体细胞偏向增殖且MMP降解酶表达降低，男性供体细胞偏向染色质重塑与发育信号，并存在一个包含TIMP1降低的较小共有成分。鉴于每名性别仅由一名供体代表，这些发现属于假设生成性质，但为解释OA患病率的性别差异、开展性别分层的半月板生物学研究指明了候选通路，具有重要意义。

研究人员重新分析了存储于基因表达综合数据库（Gene Expression Omnibus，GEO）的公共RNA测序数据集GSE199087。样本队列来源于经国家疾病研究交换机构（National Disease Research Interchange）获取的遗体捐献者内侧半月板，包含一男一女各一名47岁供体，每组三个生物学重复；数据集中另有16岁与17岁女性供体，因无年龄匹配男性对应者且文库构建与过滤方法不同，未纳入本次分析。关键技术方法包括：对转录本丰度估计进行基因水平汇总后，采用edgeR的TMM法归一化、voom转换，再以线性模型结合经验贝叶斯调节进行差异表达分析；整合基因集富集分析（GSEA，fgsea实现）、Enrichr及iLINCS三种互补平台鉴定通路，仅保留三平台方向一致的基因本体（Gene Ontology，GO）条目；通过前沿基因分析筛选各富集通路的核心贡献基因；最后对两位供体的log₂倍数变化（log2 fold-change，log2FC）进行Spearman相关分析与二项检验。

研究结果部分如下。2.1 差异基因表达的火山图分析：火山图分析显示，女性供体MFCs经E2处理后，显著上调基因涉及细胞周期调控与DNA修复（如MCM8、RAD51、TTK、BRCA1、BUB1B、CDK2），显著下调基因则与炎症及ECM降解相关（如CXCL12、LBP、IL12A、MYD88、MMP1）；男性供体MFCs呈现截然不同的响应，显著上调基因涉及突触组装、染色质重塑与发育信号（如PLXNA3、SRCAP、NOTCH1、SHANK3、LRRC4），显著下调基因则与免疫应答及线粒体功能相关（如HLA-B、HLA-E、NDUFA3、ATP5F1D、OAS1、IFI6）。2.2 GSEA、Enrichr和iLINCS的重叠通路：三种平台共同支持女性供体细胞上调"DNA复制"（GSEA归一化富集分数NES=2.138）与"细胞周期进程正调控"（NES=1.778），下调"肽酶活性"与"成纤维细胞生长因子受体信号通路"；男性供体细胞共同支持上调"突触组装"（NES=1.977）与"染色质重塑"（NES=1.766），下调"MHC I类分子内源性肽抗原加工与呈递"（NES=-2.067）及"肽酶活性调控"（NES=-1.547），其中抗原呈递相关条目仅在男性供体细胞中达到三平台重叠。2.3 各供体的前沿基因：女性供体细胞中，细胞周期通路的核心贡献基因包括RGCC、CDCA8、CDC45与CDC6，DNA复制通路包括DNA2与CDT1，下调的防御应答调控通路以MYD88为前沿基因，下调的肽酶活性通路以MME为前沿基因；男性供体细胞的突触组装通路前沿基因包括PLXNA3、SHANK3、NOTCH1与LRRC4，染色质重塑通路包括ERCC6、BICRAL、BICRA与SUPT16H，下调的MHC I类抗原呈递通路包括HLA-E、HLA-A与HLA-B，肽酶活性调控通路以下调的TIMP1、SPINT2等为代表。2.4 两位供体的基因水平比较：在13,718例共同检出基因中，402例在两位供体中均达显著性阈值，其中344例（85.6%）变化方向一致，全体共享基因的log₂FC呈正相关（Spearman ρ=0.32，p<0.001）；TIMP1在两位供体细胞中均下调，MMP1、MMP2、IL12A、MYD88、MME仅在女性供体细胞中显著变化，HLA-B、NOTCH1、SRCAP、SHANK3、PLXNA3、LRRC4仅在男性供体细胞中显著变化；雌激素受体基因方面，ESR1在男性供体细胞中因E2处理而上调，ESR2仅在女性供体细胞中通过表达过滤，GPER1在两者中均无显著变化。

讨论部分指出，女性供体细胞中E2诱导的增殖与DNA修复基因上调、分解代谢与炎症基因下调的组合，与关节软骨中已报道的雌激素软骨保护作用相吻合，提示雌激素的抗分解代谢作用可能从关节软骨延伸至半月板纤维软骨组织；男性供体细胞中NOTCH1的上调尤其值得关注，鉴于Notch信号在软骨病理生理中的核心地位，该改变可能代表促进软骨细胞肥厚性分化的潜在不良信号。作者强调，由于每种性别仅由一名供体代表，任何单一改变均不能确证归因于性别，最可能代表可复现性别相关响应的候选改变须同时满足"在一个供体中显著而在另一供体中基本不变、为三平台支持的前沿基因、并有独立软骨细胞证据支持"三项条件，其中女性细胞中MMP1的下调完全符合全部条件。讨论同时指出，若观察到的模式确实反映性别差异，则有助于解释女性绝经前后OA发病率的显著上升，并为将性别作为生物学变量纳入雌激素与半月板研究提供依据。研究亦坦承多项局限：单一供体设计严重限制普适性，两位供体报告的种族差异构成潜在混杂因素，所用E2浓度（240 μM）高于生理循环水平，单层培养无法重现体内微环境，结果仅基于转录组且未经定量PCR或蛋白水平验证，基因层面采用未校正多重比较的名义p值，相关结果应视为待验证的候选发现。

研究结论部分翻译如下：本研究的重新分析表明，17β-雌二醇在来自女性与男性供体的人半月板纤维软骨细胞中诱发了不同的转录组反应，同时存在一个包含TIMP1表达下调的较小共有成分。女性供体的细胞激活了增殖、DNA修复与细胞周期通路，同时抑制了炎症及ECM降解基因（如MMP1与MMP2），该程序预计会随绝经时雌激素水平下降而减弱，尽管本研究未直接验证这一推论。男性供体的细胞则启动了染色质重塑与发育信号通路，NOTCH1为主要前沿基因之一，并抑制了抗原呈递通路。由于每个性别仅由一名供体代表，上述差异目前尚不能归因于性别。这些发现提示，关节软骨中报道的雌激素抗分解代谢作用可能延伸至半月板细胞，并确定了需在大规模、性别均衡队列中验证的候选通路。