数据收集与预处理
从GEO数据库获取IPF相关微阵列数据集,并合并去批次效应,为后续分析提供基础。
从GEO数据库获取GSE24206、GSE35145、GSE53845、GSE68239作为训练集,GSE70866作为验证集。使用R包inSilicoMerging合并训练集,并采用Johnson等的方法消除批次效应。
Contribution of cuproptosis and immune-related genes to idiopathic pulmonary fibrosis disease.
基于GEO数据集进行差异表达分析、WGCNA、PPI网络和ROC评估,并通过免疫组化及单细胞RNA-seq进行验证,但缺乏动物或细胞功能实验。
GEO微阵列数据集(GSE24206、GSE35145、GSE53845、GSE68239和GSE70866) 人类肺组织标本(5例IPF、3例对照) 单细胞RNA-seq图谱(75,613个细胞)
训练集:GSE24206、GSE35145、GSE53845、GSE68239,99个样本(71例IPF、28例对照) 验证集:GSE70866,196个肺泡灌洗液标本(176例IPF) DEG筛选阈值:调整后p<0.05且|log2FC|>1.5 WGCNA软阈值power=4,黄色模块与IPF显著相关(r=0.63, p<0.001) 免疫组化验证五种蛋白(CFH、STEAP1、HDC、FMO5)的表达
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
从GEO数据库获取IPF相关微阵列数据集,并合并去批次效应,为后续分析提供基础。
从GEO数据库获取GSE24206、GSE35145、GSE53845、GSE68239作为训练集,GSE70866作为验证集。使用R包inSilicoMerging合并训练集,并采用Johnson等的方法消除批次效应。
鉴定IPF与正常对照之间的差异表达基因,并探讨其功能富集。
使用limma包进行差异分析,筛选标准为调整后p<0.05且|log2FC|>1.5,共获得92个DEGs。进行GO和KEGG富集分析。
比较IPF与正常对照之间22种免疫细胞的浸润比例差异。
使用CIBERSORT工具评估22种免疫细胞类型的比例,并比较两组间的差异。
识别与IPF最相关的基因模块,缩小候选基因范围。
使用WGCNA构建共表达网络,软阈值设为4,共获得11个模块,其中黄色模块与IPF显著相关(r=0.63, p<0.001),包含291个基因。
筛选与铜死亡相关的核心基因,并评估其诊断价值。
从GeneCard和Tsvetkov等研究中获取铜死亡相关基因,与DEGs和WGCNA黄色模块基因取交集,鉴定出五个核心基因(NUDT16、FMO5、CFH、HDC、STEAP1)。通过ROC曲线评估诊断价值。
筛选免疫相关核心基因并评估其诊断价值。
从ImmPort数据库获取免疫相关基因,与DEGs取交集获得24个DEIRGs。构建PPI网络,利用CytoNCA识别四个核心基因(CXCL12、CXCR2、CTSG、SPP1)。通过ROC曲线评估诊断价值。
验证核心基因在蛋白水平的表达情况,与转录组结果相互印证。
采用免疫组化检测IPF肺组织和正常肺组织中CFH、STEAP1、HDC、FMO5、CXCL12、CXCR2、SPP1的表达。
在单细胞水平上探究核心基因在不同细胞类型中的表达分布。
对5例IPF和3例对照肺组织进行单细胞RNA测序,分析核心基因在细胞亚群中的表达差异。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 数据预处理 | R | -- | -- |
| inSilicoMerging | -- | -- | |
| 差异表达分析 | limma | -- | 3.40.6 |
| 功能富集分析 | org.Hs.eg.db | -- | 3.1.0 |
| clusterProfiler | -- | 3.14.3 | |
| 免疫浸润分析 | CIBERSORT | -- | -- |
| 蛋白质互作网络 | STRING数据库 | -- | -- |
| 蛋白质互作网络可视化 | Cytoscape | -- | 3.9.1 |
| hub基因筛选 | CytoNCA插件 | -- | -- |
| 免疫组化 | STEAP1抗体 | Bioss | bs-1901R |
| CFH抗体 | Bioss | bs-9525R | |
| HDC抗体 | Bioss | bs-1054R | |
| FMO5抗体 | Bioss | bs13187R | |
| CXCL12抗体 | Bioss | bs-4938R | |
| CXCR2抗体 | Absin | abs-133162 | |
| SPP1抗体 | Affinity | AF0227 | |
| 二抗 | OriGene | PV-6000 | |
| 单细胞解离 | GEXSCOPE®组织解离液 | -- | -- |
| 磷酸盐缓冲液 | HyClone | -- | |
| GEXSCOP®红细胞裂解液 | Singleron | -- | |
| 台盼蓝 | Sigma | -- | |
| 单细胞转录组文库构建 | GEXSCOPE®单细胞RNA文库构建试剂盒 | Singleron Biotechnologies | -- |
| 单细胞测序 | Illumina Novaseq6000测序平台 | Illumina | -- |
| 单细胞测序数据分析 | fastp | -- | V1 |
| STAR | -- | v2.5.3a | |
| featureCounts | -- | v1.6.2 | |
| Seurat包 | -- | 3.0.1 | |
| 统计分析 | SPSS | -- | 22.0 |
| ROC分析 | pROC | -- | 1.17.0.1 |
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 数据预处理 | 训练集数据集ID、样本数、批次效应去除方法 阅读提示:Materials and methods: Collection and data processing of microarray datasets for IPF |
| 差异表达分析 | DEG筛选阈值(调整后p<0.05且|log2FC|>1.5) 阅读提示:Materials and methods: Collection and data processing of microarray datasets for IPF |
| WGCNA | 软阈值power=4、模块合并距离<0.25 阅读提示:Materials and methods: Construction of WGCNA coexpression network |
| 免疫组化 | 抗体品牌、货号、修复条件、孵育时间 阅读提示:Materials and methods: Immunohistochemistry |
| 单细胞RNA-seq | 组织样本数、细胞数、解离方法、文库构建试剂盒、测序平台、聚类分辨率 阅读提示:Materials and methods: Tissue dissociated single cell suspension, Sequencing of single-cell transcriptome libraries, Quality control and analysis of single-cell sequencing data |