m6A调节因子介导的RNA甲基化修饰模式调节骨关节炎中的免疫微环境

m6A Regulator-Mediated RNA Methylation Modification Patterns Regulate the Immune Microenvironment in Osteoarthritis.

作者信息Yang Duan, Cheng Yu, Meiping Yan, Yuzhen Ouyang, Songjia Ni
PMID35812736
发布时间2022-06-23
DOI10.3389/fgene.2022.921256

实验完整度

包含生信分析、RT-qPCR和Western blot验证,有体外细胞实验(IL-1β诱导软骨细胞)验证表达,但缺乏功能验证和体内实验。

主要模型

人软骨组织样本(正常和OA) IL-1β诱导的软骨细胞体外模型

重点核对

GSE114007数据集,38个样本(18对照,20 OA) IL-1β浓度:40 ng/ml,处理48小时 软骨细胞系:CP-H107(Procell) RT-qPCR和Western blot验证基因表达

摘要

表观遗传调控,尤其是RNA N6-甲基腺苷(m6A)修饰,在免疫应答中发挥重要作用。然而,m6A在骨关节炎(OA)免疫微环境中的调节作用仍不清楚。因此,我们系统研究了38个样本中23个m6A调节因子介导的RNA修饰模式,并探讨了m6A修饰的免疫微环境特征。接下来,我们构建了一个新的OA m6A列线图、一个m6A-转录因子-miRNA网络和一个药物网络。健康与OA样本表现出不同的m6A调节因子表达模式。YTHDF3在OA样本中表达上调,并与II型辅助细胞和TGFb家族成员受体呈正相关。此外,23个m6A调节因子介导了三种不同的RNA修饰模式;在模式3中,YTHDF3、II型T辅助细胞和TGFb家族成员受体的表达水平上调。内质网氧化应激和线粒体自噬相关通路与模式3和23个m6A调节因子相关。通过RT-qPCR,我们验证了SREBF2和EGR1作为YTHDF3和IGF2BP3的转录因子与OA的发展密切相关,hsa-miR-340作为YTHDF3和IGF2BP3的miRNA参与OA的发展,我们还通过western blotting检测了IGF2BP3、YTHDF3、EGR1和SREBF2的蛋白表达水平,结果与PCR一致。总之,构建的列线图有助于预测OA风险。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
m6A调节因子介导的RNA甲基化修饰模式是否以及如何调节骨关节炎的免疫微环境?
核心机制
YTHDF3和IGF2BP3在OA中高表达,与免疫细胞浸润和TGFb信号相关,可能通过内质网应激和线粒体自噬通路影响OA进展。
主要证据
从GEO数据集(GSE114007)分析23个m6A调节因子表达,鉴定三种m6A修饰模式;通过RT-qPCR和Western blot验证IGF2BP3、YTHDF3、EGR1、SREBF2在IL-1β诱导的软骨细胞中的表达变化。
研究意义
揭示了m6A修饰在OA免疫微环境中的潜在调控机制,为OA的临床诊断提供了参考列线图。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

获取公开数据并进行差异表达分析

识别OA与正常样本中m6A调节因子的表达差异。

从GEO数据库下载GSE114007数据集,包含38个样本,使用DEseq2包分析m6A基因差异表达,并绘制火山图、热图等。

2

构建m6A预测模型

评估m6A调节因子对OA的诊断能力,并构建风险预测模型。

使用随机森林筛选关键m6A调节因子,再用LASSO回归构建预测模型和列线图。

3

免疫微环境关联分析

探索m6A调节因子与免疫细胞浸润及免疫应答的关系。

采用ssGSEA估计免疫细胞浸润和免疫反应活性,并计算与m6A调节因子的相关性。

4

m6A修饰模式识别

根据m6A调节因子表达鉴定不同的m6A修饰模式。

使用ConsensusClusterPlus对OA样本进行无监督聚类,确定最佳k值,并进行PCA验证。

5

通路富集分析

揭示不同m6A修饰模式下的生物学通路差异。

使用GSVA对HALLMARK和GO基因集进行富集分析,比较不同模式间的通路活化差异。

6

体外细胞实验验证

验证生信分析发现的m6A调节因子及转录因子、miRNA的表达变化。

使用IL-1β诱导软骨细胞模拟OA,通过RT-qPCR和Western blot检测相关基因和蛋白的表达。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
TRIzol试剂--15,596–026;
PrimeScript RT试剂盒(含gDNA Eraser)TakaraRR047A
SYBR Green qPCR混合液BeyotimeD7260
全细胞裂解液KeyGEN Biotech--
PVDF膜Millipore--
ECL化学发光液KeyGEN Biotech--
DMEM培养基----
0.25%胰蛋白酶----
磷酸盐缓冲液----
软骨细胞(CP-H107)ProcellCP-H107
白细胞介素-1βPEPROTECH--
IGF2BP3抗体Thermo FisherMA5-27484;
YTHDF3抗体Thermo FisherPA5-107309;
EGR1抗体Thermo FisherMA5-15008;
SREBF2抗体Thermo Fisher1:2000
GAPDH抗体Beyotime Biotechnology1:1000

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
数据获取与差异分析
数据集GSE114007,包含38个样本(18对照,20 OA);差异分析使用DEseq2
阅读提示:Materials and Methods: Data Acquisition and Difference Analysis
m6A预测模型构建
随机森林和LASSO回归的具体参数;训练和验证集划分比例(70/30)
阅读提示:Materials and Methods: Analysis of m6A Regulatory Factors in OA; Results: Construction of an m6A Gene Prediction Model
免疫浸润分析
ssGSEA分析的基因集来源(ImmPort数据库)
阅读提示:Materials and Methods: Analysis of the Correlation Between m6A Regulators and Immune Characteristics
m6A修饰模式聚类
聚类数k=3选择依据,ConsensusClusterPlus参数
阅读提示:Materials and Methods: Identification of m6A Modification Pattern; Results: Modification Patterns of m6A RNA Methylation
细胞培养与处理
软骨细胞来源(Procell, CP-H107);IL-1β浓度40 ng/ml,处理48小时
阅读提示:Materials and Methods: Cell Culture
基因表达验证
RT-qPCR引物序列(补充表S1);Western blot抗体信息及稀释比例
阅读提示:Materials and Methods: RT-qPCR and Western Blotting