细胞类型特异性eQTL揭示复杂性状的遗传架构

《Nature》：Cell-type-specific eQTLs underlie the genetic architecture of complex traits

　　

复杂性状的遗传效应主要通过调控基因表达实现；然而该过程尚不清楚。对基因表达遗传效应（表达数量性状位点（eQTL））的研究可揭示遗传变异与复杂性状之间的基因调控层。既往研究未能有效捕获细胞类型特异性eQTL，而后者可能对复杂性状十分重要。研究人员通过将方差组分模型应用于群体规模单细胞RNA测序（scRNA-seq）数据，无偏地表征了细胞类型特异性eQTL。利用OneK1K队列的外周血单核细胞，研究人员证实细胞类型特异性eQTL富集复杂性状遗传力，而在细胞类型共享eQTL中未观察到此现象。研究人员还发现eQTL特异性与具有更强选择约束、增强子复杂性更高及基因网络连通性更强的基因相关，这三种特征在复杂性状中相对于已知eQTL富集。转录组范围内，反式（trans）eQTL多为细胞类型特异性（60%特异），而顺式（cis）eQTL多为共享（30%特异）。研究人员使用第二个scRNA-seq数据集复现发现，并证明细胞类型共享与特异性eQTL在不同祖源间一致。结果表明eQTL细胞类型特异性是基因调控的关键特征，并部分解释为何已知eQTL在复杂性状的基因调控效应中匮乏。

研究背景方面，连接复杂性状遗传效应与其潜在基因调控机制的挑战推动了利用批量组织RNA测序（RNA-seq）开展群体规模基因表达研究（如GTEx）。尽管已鉴定数千个贡献转录丰度变异的eQTL，已知eQTL仅解释复杂性状遗传效应的约10%–30%，提示基因调控效应可能依赖批量组织未显现的上下文（如细胞类型、环境或发育阶段）。既往研究聚焦于批量组织中达统计显著的eQTL，因效能有限仅识别少量疑似eQTL，且偏向跨细胞类型共享的超大效应，这些效应在复杂性状遗传架构与相关基因特征中匮乏。近期scRNA-seq的eQTL研究提升了识别细胞类型特异性eQTL的效能，但仍聚焦显著eQTL，整体角色不清。研究人员提出基于scRNA-seq数据变异分解的替代方案，无偏量化细胞类型共享与特异性eQTL的总体贡献，确立其在复杂性状遗传架构中的重要性，并部分解释已知eQTL在功能相关基因特征中匮乏的原因。该论文发表于《Nature》。

主要关键技术方法方面，研究人员开发细胞类型知情遗传混合模型分析（CIGMA），将GREML（基于基因组亲缘关系的约束最大似然）思路扩展至scRNA-seq伪批量数据，联合拟合细胞类型共享随机效应α_l与细胞类型c特异性随机效应γ_lc，并预扣除细胞内变异δ_ic；利用OneK1K队列（928人、7种PBMC细胞类型、平均约1190细胞/人、10288基因、cis SNP取基因体±500 kb且MAF>5%）估计σ_g²与v_c及特异性\bar{v}/(σ_g²+ \bar{v})；以基因型置换、细胞置换及伪批量模拟校准无偏性与效能；辅以CLUES（含SLE患者与欧洲祖源对照的PBMC百万级细胞）与ImmVar作跨祖源复现；下游用LDSC（连锁不平衡分数回归）检验顶级cs-eGene对10个复杂性状遗传力富集，并用抽象中介模型量化介导比例。

研究结果部分，CIGMA量化细胞类型特异性eQTL小节指出，模拟显示CIGMA在无eQTL时校准良好，细胞置换保留共享eQTL且特异性估计为空，随个体数N增加无偏且精度提升，检测cs-eGene效能从N=1000约40%升至N=2000近100%；相较忽略δ的CIGMA、GREML（CTP/OP）、GxEMM、BOLT-REML、OTD，CIGMA因建模细胞间变异与共享/特异双路径，在真实scRNA-seq下无偏。OneK1K中的细胞类型特异性eQTL小节指出，实测鉴定193个显著cs-eGene（如CTLA4、SMDT1），cis eQTL特异性31.1%、trans eQTL特异性59.4%，特异性随cis窗口距TSS增大而升（5 kb时30%至1 Mb时41%），OP混合使cis特异性降至5%。eQTL特异性与基因特征小节指出，约束更强（LOEUF更低）基因共享与特异cis遗传方差更低但特异性更高，高特异性伴更多增强子与网络连通性（共享eQTL富集于最不连通基因十分位数），高特异eGene GO富集于免疫细胞项，低特异富集于细胞质翻译。特异性eQTL构成复杂性状基础小节指出，顶级200 cs-eGene在7个血相关性状中4个富集（LDSC P≤0.05），共享eGene与批量eGene无富集；抽象中介显示cs-eGene介导血相关性状遗传力9.0%（五倍于其他基因），CD4_ET富集克罗恩与溃疡性结肠炎，CD8_ET与NK富集SLE。CLUES与ImmVar中的eQTL特异性小节指出，CLUES中cis中位遗传力4%、特异性63%（高于OneK1K 34%），LOEUF/增强子/连通性关联方向一致，OneK1K cs-eGene在CLUES中P值显著富集（36/183 P<0.05），σ_g²与\bar{v}跨队列相关（均P<1×10^-16，均值r=0.27）。

讨论部分总结，已有理论称已知eQTL偏共享大效应而复杂性状由弱细胞类型特异性eQTL驱动；本研究无偏量化支持该假说：细胞类型特异性eQTL富集关键基因特征与复杂性状遗传力，而共享eQTL不富集，提示复杂性状“缺失的基因调控”主要由需超大样本才可单个检出的eQTL解释。CIGMA优势在消除检出效能偏差，局限在估计更噪、不能共定位GWAS、未计非高斯与复杂细胞关系、伪批量忽略胞内异质；更细细胞状态分区或进一步提升特异性富集。结论翻译：一种解释已知eQTL与复杂性状遗传效应有限重叠的主导理论认为，已知eQTL主要是跨细胞类型共享的大效应，而复杂性状主要由细胞类型特异性的弱eQTL驱动。研究人员利用scRNA-seq对eQTL的无偏量化支持该假说：细胞类型特异性eQTL富集于关键基因特征与复杂性状遗传力，不同于细胞类型共享eQTL。这表明复杂性状下缺失的基因调控（如同复杂性状缺失遗传力）将主要由无海量样本无法单个检出的eQTL解释。研究人员还确立eQTL细胞类型特异性是稳健基因特征，trans相比cis约两倍富集，且与远端cis调控、更强选择约束、更多增强子及更高基因网络连通性相关。