数据挖掘与差异基因筛选
从GEO数据集识别COPD相关的差异表达基因,并与线粒体基因取交集,获得线粒体相关差异表达基因(MitoDEGs)。
分析GSE38974数据集,使用limma包筛选|log2FC|>1和adj.p<0.05的DEGs,与MitoCarta3.0数据库取交集,获得28个MitoDEGs。
Identifying the Key Mitochondria-Related Genes in COPD by Integrating Machine Learning and Bioinformatics Analyses.
包含GEO数据集的机器学习分析、体外细胞CSE刺激实验(qRT-PCR和Western blot)和体内CS诱导小鼠模型验证,并利用IHC和IF验证免疫浸润,涉及多个独立证据层级。
支气管上皮细胞系(Beas-2B和16HBE) CS诱导COPD小鼠模型(C57BL/6J小鼠) GEO数据集(GSE38974和GSE8545)
CSE处理浓度:0%–7% CS暴露时间:6个月,每天10支香烟,5天/周 qRT-PCR引物序列(BAX、DLST、GAPDH) 抗体稀释比例(BAX 1:2000,DLST 1:1000等) ROC分析AUC阈值:>0.7
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
从GEO数据集识别COPD相关的差异表达基因,并与线粒体基因取交集,获得线粒体相关差异表达基因(MitoDEGs)。
分析GSE38974数据集,使用limma包筛选|log2FC|>1和adj.p<0.05的DEGs,与MitoCarta3.0数据库取交集,获得28个MitoDEGs。
探索MitoDEGs在GO和KEGG通路中的富集情况,以揭示其潜在功能。
对28个MitoDEGs进行GO(BP、CC、MF)和KEGG富集分析,使用clusterProfiler包,adj.p<0.05为显著。
利用LASSO和SVM-RFE算法从MitoDEGs中筛选与COPD诊断相关的特征基因。
对28个MitoDEGs应用LASSO回归(10折交叉验证)和SVM-RFE(10折交叉验证),取两种算法的交集,得到5个特征基因:BAX、DLST、FKBP10、FUNDC2、RMDN1。
在验证集GSE8545中验证特征基因的表达差异和诊断效能。
在GSE8545中比较特征基因在COPD和对照组的表达,并绘制ROC曲线,AUC>0.7表示诊断效能。
验证CSE刺激下人支气管上皮细胞中BAX和DLST的表达变化。
用0%–7%浓度的CSE刺激Beas-2B和16HBE细胞,通过qRT-PCR和Western blot检测BAX和DLST的mRNA和蛋白表达。
在CS诱导的COPD小鼠模型中验证BAX和DLST的表达变化及肺组织病理变化。
建立CS诱导的COPD小鼠模型(C57BL/6J小鼠,CS暴露6个月),通过H&E染色评估肺气肿(MLI和DI),Western blot检测BAX和DLST蛋白表达。
通过CIBERSORT分析COPD样本中的免疫细胞浸润情况,并评估其与枢纽基因的相关性。
使用CIBERSORT算法分析GSE38974的免疫细胞浸润,比较22种免疫细胞比例,并分析BAX、DLST与免疫细胞的相关性。
在CS诱导的COPD小鼠肺组织中验证生物信息学发现的免疫细胞浸润和枢纽基因的相关性。
通过IHC检测F4/80、Ly6C、CD8a、CD206的表达,通过IF共染BAX/DLST与F4/80或NKp46/GranzymeB,量化共定位。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 细胞培养 | DMEM培养基 | -- | -- |
| RPMI-1640培养基 | -- | -- | |
| 胎牛血清 | -- | -- | |
| qRT-PCR | SPARKeasy细胞RNA提取试剂盒 | SparkJade Science Co. | -- |
| Evo M-MLV逆转录预混液 | Accurate Biotechnology | -- | |
| SYBR Green Premix Pro Taq HS qPCR试剂盒 | Accurate Biotechnology | -- | |
| Western blot | RIPA裂解液 | Solarbio Technology | -- |
| 蛋白酶抑制剂 | Beyotime | -- | |
| BCA蛋白定量试剂盒 | Beyotime | -- | |
| PVDF膜 | Millipore | -- | |
| BAX抗体 | Abmart | -- | |
| DLST抗体 | Boster | -- | |
| GAPDH抗体 | Proteintech | -- | |
| 抗兔IgG | Proteintech | -- | |
| 抗小鼠IgG | Proteintech | -- | |
| 超敏增强化学发光试剂 | Yeasen Biotechnology | -- | |
| H&E染色 | 伊红 | Servicebio | -- |
| 中性树脂 | Beyotime | -- | |
| IHC | F4/80抗体 | Invitrogen | -- |
| Ly6C抗体 | BioLegend | -- | |
| CD8a抗体 | Abways | -- | |
| CD206抗体 | Abways | -- | |
| IF | 酪胺信号放大试剂盒 | Servicebio | -- |
| NKp46抗体 | Proteintech | -- | |
| 颗粒酶B抗体 | Proteintech | -- | |
| DAPI | -- | -- | |
| 仪器 | BioRadar PCR CFX-Connect实时荧光定量PCR系统 | BioRadar | -- |
| Tanon全自动化学发光图像分析系统 | Tanon | -- | |
| 光学显微镜 | Olympus | -- | |
| 3DHISTECH扫描仪 | 3DHISTECH | -- | |
| 酶标仪 | -- | -- |
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 数据挖掘与分析 | 数据集GSE38974和GSE8545的样本获取、处理平台和批次效应校正方法 阅读提示:参见Methods 2.1数据收集,以及图1和图5、6的图注 |
| 细胞实验 | CSE制备方法(包括香烟品牌、浓度范围0%–7%)、细胞系(Beas-2B和16HBE)的培养条件 阅读提示:参见Methods 2.6 CSE制备和2.7细胞培养 |
| 动物实验 | 小鼠品系、性别、周龄、CS暴露时间和频率 阅读提示:参见Methods 2.9动物模型和Pauwels等参考文献 |
| 分子检测 | qRT-PCR引物序列、退火温度、Western blot抗体稀释比 阅读提示:参见Methods 2.8、2.10和表1 |
| 免疫浸润分析 | CIBERSORT参数(LM22、1000 permutations)和样本筛选标准(p<0.05) 阅读提示:参见Methods 2.5和2.14统计部分 |