常见的顺式调控变异修改颅面微 somia 中致病性 SHROOM3 变异的外显率

Common cis-regulatory variation modifies the penetrance of pathogenic SHROOM3 variants in craniofacial microsomia

作者信息Hao Zhu, Jiao Zhang, Soumya Rao, Matthew D Durbin, Ying Li, Ruirui Lang, Jiqiang Liu, Baichuan Xiao, Hailin Shan, Ziqiu Meng, Jinmo Wang, Xiaokai Tang, Zhenni Shi, Liza L Cox, Shouqin Zhao, Stephanie M Ware, Tiong Y Tan, Michelle de Silva, Lyndon Gallacher, Ting Liu, Jie Mi, Changqing Zeng, Hou-Feng Zheng, Qingguo Zhang, Stylianos E Antonarakis, Timothy C Cox, Yong-Biao Zhang
PMID40234029
发布时间2025-05
DOI10.1101/gr.280047.124

实验完整度

包含人群遗传关联分析、单倍型分析、双荧光素酶报告基因实验、以及两种 Shroom3 小鼠模型的定量表型分析,证据层级覆盖人类遗传数据和体内功能验证。

主要模型

HEK-293T 细胞 Shroom3 基因陷阱小鼠模型(Shroom3gt/gt 及 Shroom3+/gt) Shroom3 敲除小鼠模型(Shroom3em1(IMPC)Bay/em1(IMPC)Bay)

重点核对

SHROOM3 位点预测致病性变异的单倍型与 eQTL rs10017322 的等位基因组合 eQTL rs344131 位于 SHROOM3 增强子内,其风险等位基因 G 对增强子活性的影响(荧光素酶活性降低) eQTL 组合 rs344131 和 rs61090632 对增强子活性的协同效应 Shroom3 基因剂量(纯合、杂合缺失)对小鼠耳廓、下颌骨和上颌骨长度及不对称的影响 CFM 患者中携带四个风险等位基因的 eQTL 组合富集于患者组

摘要

致病性编码变异已在数千个基因中被鉴定,但携带这些变异的个体中不完全外显的机制仍知之甚少。在本研究中,在一个由 2009 名中国血统颅面微 somia(CFM)患者和 2625 名汉族对照组成的队列中,我们在 CFM 患者和健康个体中均鉴定了 SHROOM3 中的多个预测致病性编码变异。我们发现 CFM 的外显率与包含可能致病性编码 SHROOM3 变异和 CFM 相关 SHROOM3 表达数量性状位点(eQTL)的特定单倍型组合相关。进一步的研究涉及特定的 eQTL 组合,如 rs1001322 或 rs344131,与其他显著的 CFM 相关 eQTL 组合,我们称之为组合 eQTL 表型修饰因子(CePMods)。我们还显示,位于 SHROOM3 调控增强子区域内的 rs344131 对增强子活性具有等位基因特异性效应,从而影响携带任何罕见编码变异的 SHROOM3 等位基因的表达水平。我们的发现还表明,即使在 SHROOM3 中没有罕见的 deleterious 编码变异的情况下,CePMods 也可能作为致病决定因素。这突显了等位基因表达在决定颅面异常(包括小耳畸形和面部不对称)的外显率和严重程度中的关键作用。此外,使用定量表型分析,我们证明小耳畸形和面部不对称均存在于两个独立的 Shroom3 小鼠模型中,其严重程度取决于基因剂量。我们的研究确立了 SHROOM3 作为 CFM 的可能致病基因,并证明 eQTL 是携带可能致病性罕见编码变异的个体中疾病表现外显率修饰的决定因素。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
常见顺式调控变异(eQTL)如何修饰 SHROOM3 致病性编码变异在颅面微 somia(CFM)中的外显率,以及 SHROOM3 表达水平对颅面表型的影响。
核心机制
预测致病性编码变异与高表达等位基因(如 rs10017322 的 G 等位基因)顺式组合时,可能增加致病等位基因的表达水平,从而提高外显率;而低表达等位基因组合则可能降低外显率。此外,eQTL 组合(CePMods),如 rs344131 和 rs61090632 的风险等位基因组合,可通过降低增强子活性来减少 SHROOM3 表达,从而在没有罕见编码变异的情况下也可能导致 CFM 表型。
主要证据
人群遗传关联和单倍型分析:在 CFM 患者中,所有预测致病性变异均与 rs10017322 的高表达 G 等位基因顺式组合,而在健康对照中与低表达 A 等位基因组合;荧光素酶实验证明 rs344131 的 G 等位基因和 rs344131 与 rs61090632 的风险等位基因组合降低增强子活性;两种 Shroom3 小鼠模型(基因陷阱和敲除)显示纯合子和杂合子均出现小耳畸形和面部不对称,且严重程度与基因剂量相关。
研究意义
该研究将 SHROOM3 确立为 CFM 的可能致病基因,并揭示了 eQTL 和 CePMods 作为修饰外显率的因素,为理解不完全外显机制和遗传咨询提供了新见解。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

病例-对照关联分析

鉴定 CFM 患者和健康对照中 SHROOM3 位点的常见风险变异和预测致病性编码变异。

对 2009 名中国 CFM 患者和 2625 名健康汉人进行全基因组测序,进行单变异关联分析,并利用 SKAT 进行基因负担检验。

2

精细定位 eQTL

识别 SHROOM3 区域中可能修饰致病性变异外显率的因果 eQTL。

使用 PLINK 条件分析、CAVIAR、FINEMAP 和 PAINTOR 对关联信号进行精细定位,并结合 GTEx 数据库中的 eQTL 注释。

3

单倍型分析

确定预测致病性编码变异与 eQTL 等位基因在患者和对照中的顺式组合模式。

使用 SHAPEIT5 对 SHROOM3 区域变异进行单倍型构建,并比较 CFM 患者和健康对照中携带预测致病性变异的单倍型上 rs10017322 的等位基因。

4

eQTL 组合分析

评估多个 eQTL 风险等位基因组合(CePMods)与 CFM 表型的关联。

分析 rs10017322 和 rs344131 与其他 eQTL 的两点组合,在 CFM 患者和健康对照中比较风险等位基因组合的频率。

5

荧光素酶报告基因实验

检测个 eQTL 对 SHROOM3 增强子活性的功能性影响。

将含有不同 eQTL 等位基因的 SHROOM3 增强子序列克隆到荧光素酶载体中,转染 HEK-293T 细胞,检测双荧光素酶活性。

6

小鼠模型表型分析

评估 Shroom3 基因剂量对小鼠颅面形态(特别是耳廓、下颌骨和上颌骨)的影响。

对 Shroom3gt/gt、Shroom3+/gt 和野生型胚胎(E14.25-E15.25)进行 OPT 成像,对 Shroom3em1(IMPC)Bay/em1(IMPC)Bay 和野生型 E18.5 胚胎进行 microCT 成像,并基于标志点进行定量测量。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
DMEM 培养基(Gibco)Gibco--
胎牛血清(Gibco)Gibco--
Lipofectamine 3000 转染试剂Thermo FisherL3000008
双荧光素酶检测试剂盒TransGenFR201
CLARIOstar 微孔板读板器BMG Labtech--
Q5 定点突变试剂盒New England BioLabsE0554

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
病例-对照关联分析
队列构成:2009 例中国 CFM 患者和 2625 例健康汉人;测序平台:BGI T7;平均深度约 10×;关联分析参数:MAF>0.001,HWE p>0.001,样本 call rate>98%,SNP call rate>95%
阅读提示:Methods: Study cohorts, High-throughput sequencing and variant calling, Association study
精细定位 eQTL
所使用的软件版本及参数:PLINK 条件分析、CAVIAR v2.2、FINEMAP v1.4.2、PAINTOR v3.0;LD block 划分方法(MAF>0.1)
阅读提示:Methods: Fine mapping, LD block analysis
单倍型分析
单倍型构建软件:SHAPEIT5 v5.1.0;变异筛选标准:MAF>0.005 的常见变异使用 phase_common,MAF<0.005 的罕见变异使用 phase_rare;预测致病性变异定义标准
阅读提示:Methods: Haplotype prediction and combination analysis, Identification of likely pathogenic rare variants
荧光素酶报告基因实验
细胞系:HEK-293T(1101HUM-PUMC000091);增强子克隆载体:GV238;所用 eQTL 位点:rs10017322、rs344131、rs74918804、rs189707263、rs344126、rs61090632;转染试剂:Lipofectamine 3000;检测时间:24小时;检测试剂盒:双荧光素酶检测试剂盒
阅读提示:Methods: Luciferase assay
小鼠模型表型分析
小鼠品系:Shroom3gt(基因陷阱)和 Shroom3em1(IMPC)Bay(敲除);胚胎阶段:E14.25-E15.25(基因陷阱)和 E18.5(敲除);成像方法:OPT(Bioptonics 3001M)和 microCT(IMPC);测量标志点;统计分析:Welch's t-test
阅读提示:Methods: Mouse model, OPT and quantitative assessment of embryonic craniofacial tissue, Quantitative assessment of microCT data