单细胞转录组:解码基因差异之谜
单细胞转录组差异基因分析
单细胞转录组测序(scRNA-seq)能够解析细胞间的异质性,差异基因分析(Differential Gene Expression, DGE)是识别不同细胞群体间关键基因表达变化的核心方法。
数据预处理
原始数据需经过质量控制(QC)、标准化(如CPM、TPM或DESeq2的size factor校正)及批次效应校正(如Harmony或Seurat的CCA)。低质量细胞和基因需过滤,通常保留线粒体基因占比低、检测到足够基因数的细胞。
差异分析方法
常见工具包括MAST、DESeq2(单细胞适配版)、limma和Wilcoxon秩和检验(Seurat默认方法)。MAST适用于零膨胀数据,DESeq2适合高测序深度样本,而Wilcoxon检验对稀疏数据稳健。
关键参数设置
需调整logFC阈值(通常≥0.25)和p值校正方法(如FDR或Bonferroni)。伪批量分析(pseudobulk)可提升统计效力,尤其适用于稀有细胞类型比较。
富集分析
差异基因列表需通过功能富集分析解释其生物学意义,常用方法包括基因本体(GO)、KEGG通路和GSEA(基因集富集分析)。
数据库选择
GO分析涵盖分子功能(MF)、生物过程(BP)和细胞组分(CC)。KEGG提供通路级注释,MSigDB包含更广泛的基因集(如Hallmark集合)。
工具与流程
ClusterProfiler(R包)支持GO/KEGG富集,GSEA需排名基因(按logFC或p值)。结果可视化采用点图、网络图或富集地图(EnrichmentMap)。
注意事项
避免冗余术语(Revigo可去冗余),关注保守通路(如代谢或信号转导)。跨物种分析需基因ID转换(如biomaRt)。
整合分析示例
以Seurat流程为例:
- 差异分析:
markers <- FindMarkers(seurat_obj, ident.1 = "Cluster1", ident.2 = "Cluster2", test.use = "wilcox")
- 富集分析(ClusterProfiler):
go_result <- enrichGO(gene = markers$gene, OrgDb = org.Hs.eg.db, ont = "BP")
dotplot(go_result)
通过差异基因与富集结果的结合,可揭示细胞亚群的特异功能及其在疾病或发育中的潜在作用。
BbS.okapop082.sbs/PoSt/1122_131805.HtM
BbS.okapop083.sbs/PoSt/1122_185166.HtM
BbS.okapop084.sbs/PoSt/1122_900501.HtM
BbS.okapop085.sbs/PoSt/1122_461977.HtM
BbS.okapop086.sbs/PoSt/1122_689423.HtM
BbS.okapop087.sbs/PoSt/1122_768657.HtM
BbS.okapop088.sbs/PoSt/1122_568293.HtM
BbS.okapop090.sbs/PoSt/1122_023300.HtM
BbS.okapop091.sbs/PoSt/1122_511267.HtM
BbS.okapop092.sbs/PoSt/1122_621240.HtM
BbS.okapop082.sbs/PoSt/1122_459033.HtM
BbS.okapop083.sbs/PoSt/1122_499445.HtM
BbS.okapop084.sbs/PoSt/1122_798645.HtM
BbS.okapop085.sbs/PoSt/1122_267821.HtM
BbS.okapop086.sbs/PoSt/1122_765233.HtM
BbS.okapop087.sbs/PoSt/1122_049393.HtM
BbS.okapop088.sbs/PoSt/1122_314276.HtM
BbS.okapop090.sbs/PoSt/1122_972079.HtM
BbS.okapop091.sbs/PoSt/1122_682439.HtM
BbS.okapop092.sbs/PoSt/1122_090890.HtM
BbS.okapop093.sbs/PoSt/1122_380506.HtM
BbS.okapop094.sbs/PoSt/1122_242460.HtM
BbS.okapop095.sbs/PoSt/1122_573563.HtM
BbS.okapop096.sbs/PoSt/1122_931952.HtM
BbS.okapop097.sbs/PoSt/1122_535896.HtM
BbS.okapop098.sbs/PoSt/1122_509103.HtM
BbS.okapop099.sbs/PoSt/1122_534784.HtM
BbS.okapop100.sbs/PoSt/1122_538693.HtM
BbS.okapop101.sbs/PoSt/1122_377341.HtM
BbS.okapop102.sbs/PoSt/1122_465746.HtM
BbS.okapop093.sbs/PoSt/1122_095877.HtM
BbS.okapop094.sbs/PoSt/1122_734683.HtM
BbS.okapop095.sbs/PoSt/1122_261926.HtM
BbS.okapop096.sbs/PoSt/1122_886729.HtM
BbS.okapop097.sbs/PoSt/1122_021469.HtM
BbS.okapop098.sbs/PoSt/1122_498802.HtM
BbS.okapop099.sbs/PoSt/1122_204352.HtM
BbS.okapop100.sbs/PoSt/1122_388393.HtM
BbS.okapop101.sbs/PoSt/1122_335892.HtM
BbS.okapop102.sbs/PoSt/1122_235996.HtM
BbS.okapop093.sbs/PoSt/1122_625000.HtM
BbS.okapop094.sbs/PoSt/1122_940435.HtM
BbS.okapop095.sbs/PoSt/1122_976949.HtM
BbS.okapop096.sbs/PoSt/1122_006067.HtM
BbS.okapop097.sbs/PoSt/1122_611130.HtM
BbS.okapop098.sbs/PoSt/1122_784245.HtM
BbS.okapop099.sbs/PoSt/1122_266208.HtM
BbS.okapop100.sbs/PoSt/1122_595930.HtM
BbS.okapop101.sbs/PoSt/1122_324618.HtM
BbS.okapop102.sbs/PoSt/1122_322643.HtM
BbS.okapop093.sbs/PoSt/1122_903685.HtM
BbS.okapop094.sbs/PoSt/1122_542980.HtM
BbS.okapop095.sbs/PoSt/1122_857754.HtM
BbS.okapop096.sbs/PoSt/1122_049996.HtM
BbS.okapop097.sbs/PoSt/1122_207152.HtM
BbS.okapop098.sbs/PoSt/1122_019525.HtM
BbS.okapop099.sbs/PoSt/1122_409369.HtM
BbS.okapop100.sbs/PoSt/1122_614516.HtM
BbS.okapop101.sbs/PoSt/1122_156942.HtM
BbS.okapop102.sbs/PoSt/1122_267773.HtM
BbS.okapop093.sbs/PoSt/1122_281818.HtM
BbS.okapop094.sbs/PoSt/1122_505076.HtM
BbS.okapop095.sbs/PoSt/1122_398116.HtM
BbS.okapop096.sbs/PoSt/1122_305400.HtM
BbS.okapop097.sbs/PoSt/1122_986656.HtM
BbS.okapop098.sbs/PoSt/1122_717926.HtM
BbS.okapop099.sbs/PoSt/1122_946955.HtM
BbS.okapop100.sbs/PoSt/1122_339440.HtM
BbS.okapop101.sbs/PoSt/1122_173050.HtM
BbS.okapop102.sbs/PoSt/1122_299728.HtM
BbS.okapop093.sbs/PoSt/1122_304728.HtM
BbS.okapop094.sbs/PoSt/1122_660205.HtM
BbS.okapop095.sbs/PoSt/1122_378825.HtM
BbS.okapop096.sbs/PoSt/1122_929483.HtM
BbS.okapop097.sbs/PoSt/1122_473183.HtM
BbS.okapop098.sbs/PoSt/1122_147376.HtM
BbS.okapop099.sbs/PoSt/1122_240135.HtM
BbS.okapop100.sbs/PoSt/1122_365319.HtM
BbS.okapop101.sbs/PoSt/1122_174196.HtM
BbS.okapop102.sbs/PoSt/1122_897205.HtM
顺丰集团工作强度 322人发布