碱基日读 每天读完真实文献,只讲有用的那一段

孟德尔随机化联合单细胞转录组,在干燥综合征中筛出CD7和LAPTM5

2026-10-03·碱基·生信方法

要点 CD7 遗传保护效应与全血表达方向一致,但 MR 表型宽泛,小样本验证尚不充分

原发性干燥综合征(primary Sjögren’s disease, pSD)是以外分泌腺功能受损为核心的慢性自身免疫病。多年来有不少 GWAS 信号,但从染色体位点跳到「这个基因在疾病里做什么」,中间的鸿沟一直没填好。这篇发在 Clinical Rheumatology 的分析尝试用整合路线来填——孟德尔随机化(Mendelian randomization, MR)确定方向,单细胞 RNA 测序(scRNA-seq)和全血表达数据分别做独立佐证(PMID: 42821046)。

MR 跑出 10 个关联,CD7 和 LAPTM5 两处有共定位支撑

MR 的核心逻辑是借遗传变异做工具变量,估计「某基因表达升高,患病风险如何变化」。基因型出生时就确定,不受后天疾病状态影响,这条推断在方向性上比直接比较患者和对照的表达量要硬一些。

这项分析用的 eQTL 数据来自全血,疾病表型定义是电子病历里的西卡综合征/干燥综合征谱(PheCode 709.2)。跑出了 10 个名义显著的基因-表型关联,CD7 的 OR 是 0.50(P = 0.0033),意味着遗传预测 CD7 表达升高对应风险降低;LAPTM5 的 OR 是 1.29(P = 0.0088),方向相反(PMID: 42821046)。

共定位(colocalization)分析对两处都给出了支持,说明 eQTL 信号和疾病信号共享同一遗传变异,而不是相邻位点的连锁不平衡在作怪。反向 MR——把表型当暴露、基因表达当结果——结果阴性,逆向因果的解释进一步减弱。

两处局限不能绕开。第一,eQTL 来自全血,pSD 的主要靶组织是唾液腺和外分泌腺,两个组织的调控模式不一定通用。第二,PheCode 709.2 覆盖的是更宽泛的西卡谱系,不等于标准化诊断的 pSD 专病队列。10 个名义关联里真正能在专病样本里复现几个,现在还是开放的问题。

GSE51092 队列中 CD7 表达下降

MR 预测 CD7 高表达保护,那患者里 CD7 应该偏低。全血转录组数据集 GSE51092 有 190 名患者和 32 名对照,CD7 在患者中的折叠变化是 0.74,多重比较校正后 q = 2.8 × 10⁻⁴,方向和遗传预测一致(PMID: 42821046)。

这一步的价值在于:两条独立来源指向同一侧。一条是遗传预测,一条是实际观测,都说「患者 CD7 偏低」。

同一队列里还出现了一个 12 基因巨核细胞/血小板(megakaryocyte/platelet)特征,标准化均数差(standardized mean difference, SMD)为 −0.31,95% CI −0.56 到 −0.06,Welch q = 0.045。这条线不是预设假设,是分析中涌现出来的。190 对 32 的样本量不算宽裕,全血也不是 pSD 的靶部位,能不能外推到腺体里还是未知数。

单细胞数据里的巨核/血小板样群体

scRNA-seq 数据来自 GSE157278,包含 10 名捐献者的外周血单核细胞(peripheral blood mononuclear cells, PBMCs)。分析在全部 10 名捐献者里都识别出了巨核细胞/血小板样群体,用标志基因面板和差异表达分析做了表征,与此前流式细胞术文献里「系统性自身免疫病患者循环巨核细胞扩增」的报告方向一致(PMID: 42821046)。在这个群体里,CD7 和 LAPTM5 的表达变化与乳酸化评分之间的关联在六种细胞类型中最为突出,但作者自己把这放在探索层面,没有做强推断。

n = 10,全外周血,不是靶组织。这个群体在 pSD 里是否有功能性意义,目前只是关联。

乳酸化评分:五种算法,一个失败

乳酸化(lactylation)近年被认为可能参与自身免疫病的发病,是这篇文章专门加进来的一条探索线。作者从 332 篇文献里整理出 330 个乳酸化相关基因,用五种算法在 scRNA-seq 数据里为每个细胞打分,检验是否在特定细胞类型里富集(PMID: 42821046)。

五种算法给出的分数相互分歧,文库复杂度校正之后不再有稳健的细胞类型特异性信号。 样本层面,CD14⁺ 单核细胞的评分在患者中名义偏高(P = 0.016,q = 0.095),NK 细胞同样(P = 0.032,q = 0.095),但两者 q 值均未过 0.05 阈值。

方法依赖性这么强,说明测量本身还不稳定。作者的态度是克制的:乳酸化评分是转录组层面的代理指标,建立实质性联系需要在疾病相关组织里做直接的生化实验。读成阴性结果是合理的。

RT-qPCR:10 对 10,结果非显著

最后在 PBMCs 里做了 RT-qPCR,样本量是 10 名患者和 10 名对照。LAPTM5 中位折叠变化 2.9,CD7 中位折叠变化 1.7,两者均未达到统计显著性(PMID: 42821046)。

10 对 10 的把握度极低,几乎任何效应量都会淹没在噪声里。只能当方向性参考,构不成独立验证。


简讯

另两篇值得扫一眼。stFormer(PMID: 42815472)是专为空间转录组(spatial transcriptomics, ST)数据设计的 Transformer 基础模型,用交叉注意力(cross-attention)模块把空间配体基因信息引入细胞表征,并提出有偏交叉注意力(biased cross-attention)统一单细胞分辨率与低分辨率 Visium 数据之间的技术差异。预训练语料约 410 万份公开人类 Visium 数据,在细胞聚类、批次效应校正、细胞类型预测和基因功能预测上超过了已有单细胞基础模型 scFoundation,还支持通过 in silico 扰动推断细胞间配体-受体信号通路。

第三篇(PMID: 42815804)在真实实验室污染场景下验证了臭氧消毒降低 mNGS 微生物核酸背景的效果,是这类条件下的首份系统报告。对有假阳性率控制需求的宏基因组测序实验室,这是一个可以直接参考的操作数据点。

孟德尔随机化干燥综合征单细胞转录组eQTL多组学整合

同主题文章

相关工具

多重检验校正在线计算:BH (FDR)、BY 与 Bonferroni

粘贴一列 p 值,一次得到 BH 的 q 值、BY 与 Bonferroni 校正值,三者并列对照。

log2FC 与倍数变化 (Fold Change) 在线互转

log2FC ↔ 倍数、上调下调百分比,附常用阈值对照表。

比尔-朗伯定律计算器:吸光度、消光系数与浓度互算

输入吸光度、摩尔消光系数、浓度与光程中的任意三个,立刻求出第四个——酶活测定、NAD(H)动力学与小分子定量的日常换算。

DNA / RNA 拷贝数在线计算(质粒标准品、qPCR、ddPCR)

由浓度和长度算出每微升拷贝数、每管拷贝数与单拷贝质量,并生成标准曲线稀释方案。

没有广告,不收集数据,每个公式都验算过。维护它需要持续投入——如果帮你省了时间,请杯咖啡。
请杯咖啡 了解更多 →