全外显子组测序与结构建模揭示以儿童期起病为主的成人严重肥胖中罕见内源性大麻素系统变异

《Frontiers in Endocrinology》：Rare endocannabinoid system variants identified by whole-exome sequencing and structural modeling in adults with severe obesity predominantly of childhood onset

　　

背景：内源性大麻素系统（endocannabinoid system, ECS）调节能量稳态中的关键过程，但其对始于儿童期或青春期的严重肥胖的贡献仍未完全阐明。本研究旨在识别ECS基因中的罕见变异并评估其潜在结构效应。

方法：从两个专业肥胖项目中招募了93名21岁前起病（BMI ≥ 35 kg/m²）的严重肥胖成人；在有起病信息者中，97.8%报告儿童期起病。研究人员进行了全外显子组测序（whole-exome sequencing, WES），随后对ECS相关基因进行靶向分析。罕见变异（等位基因频率<1%）通过Sanger测序确认，并使用计算机模拟致病性预测和分子建模评估。研究人员还另外使用基于不同方法学原理的三种互补稳定性预测器（FoldX、Pythia和DDMut）以及AlphaMissense，对每个变异位点的全部19种可能替换进行了计算机模拟饱和突变分析。

结果：在12名个体中鉴定出9种罕见变异，其中7种被定义为意义不明变异（variant of uncertain significance, VUS）。在观察到的变异中，FAAH p.(Arg326Pro)显示出最一致的预测影响：它被FoldX列为326位最具去稳定作用的替换（ΔΔG = +3.04 kcal/mol；19种中排名第1），并且是AlphaMissense分类为可能致病（likely pathogenic）的唯一观察变异（0.884）。Arg326距离催化三联体（Lys142–Ser217–Ser241）24.3 ?，与预测的折叠稳定性效应而非活性位点直接破坏一致。相反，在DAGLA第750位，所有19种替换均被FoldX和DDMut预测为中性，并被AlphaMissense预测为可能良性，与该位点的结构耐受性预测一致。经多重检验校正后，携带者与非携带者之间没有任何人体测量或代谢参数显著差异。

结论：在严重肥胖（主要为儿童期起病）成人中鉴定出罕见ECS变异，包括一例此前未报道的FAAH p.(Arg326Pro)替换。这些计算机模拟分析是预测性的，为优先考虑候选变异（尤其是FAAH p.(Arg326Pro)）提供了计算证据；但它们本身并不能确定致病性或与表型的因果关系，这需要功能验证。

背景与目的：肥胖是复杂多因素疾病，全球约14%人口受影响，包括1.24亿儿童。除生活方式外，内源性大麻素系统（endocannabinoid system, ECS）在能量平衡、食欲和脂肪储存中发挥关键调节作用。ECS由大麻素受体CB1/CB2、内源性大麻素及其合成/降解酶组成。常见ECS基因多态性（如FAAH rs324420、CNR1 rs806381）已被报道与肥胖相关，但罕见ECS基因变异在早发严重肥胖中的作用尚不清楚。本研究旨在通过全外显子组测序（whole-exome sequencing, WES）鉴定21岁前起病的严重肥胖成人中ECS基因罕见变异，并评估其结构效应。

研究内容与结论：研究人员从巴西里约热内卢两个专科肥胖项目招募93名严重肥胖成人（体质指数（body mass index, BMI）≥35 kg/m²，起病<21岁；有起病信息者97.8%为儿童期起病），进行WES、靶向ECS基因分析、Sanger验证、计算机致病性预测和分子建模。共在12名个体中发现9个罕见变异，其中7个为意义不明变异（variant of uncertain significance, VUS）。FAAH p.(Arg326Pro)是最一致的候选变异：AlphaMissense评分为0.884（唯一可能致病），FoldX饱和突变将其列为第326位最去稳定替换，且该位点远离催化三联体，提示影响折叠稳定性而非活性位点。DAGLA p.(Ala750Val)在所有预测中均显示耐受。携带者与非携带者的人体测量和代谢参数经多重检验校正后无显著差异。研究意义在于扩展了严重肥胖中ECS罕见变异的谱系，并为功能验证提供了优先候选。

关键技术与方法：研究人员采用Illumina NovaSeq 6000平台进行WES（平均深度52×，99.3% reads比对，98.9%靶区≥20×），通过自定义虚拟面板筛选12个ECS相关基因（ABHD12、ABHD5、ABHD6、CNR1、CNR2、DAGLA、FAAH、GPR55、MGLL、NAAA、NAPEPLD、TRPV1）。罕见变异（等位基因频率<1%）经Sanger双向测序确认。致病性评估使用AlphaMissense、SpliceAI、REVEL和Franklin（ACMG分类）。结构分析使用AlphaFold预测模型或实验结构（PDB），以UCSF ChimeraX建模，并采用FoldX、Pythia、DDMut和AlphaMissense对每个变异位点进行全部19种替换的饱和突变分析。样本来自上述两个项目中符合减重手术条件的成人。

研究结果：
3.1 研究人群表型特征：93名受试者中位BMI 48.9 kg/m²，78%为女性。罕见变异携带者频率为12/93（12.9%；95% CI 7.5–21.2%）。所有携带者均报告儿童期起病。校正后无任何参数组间显著差异。
3.2 计算机模拟分析：共检出9个罕见变异，分布于FAAH（3个，含1个5'UTR）、CNR1、NAAA、DAGLA、TRPV1（2个）和NAPEPLD。7个为VUS；FAAH p.(Arg326Pro)的AlphaMissense评分最高，SpliceAI评分普遍极低，REVEL结果不一。
3.3 结构和饱和分析：FAAH Arg326Pro失去与Leu465/Asp466的氢键和Asp466盐桥，并与Asp468产生空间重叠；Pro546Leu处于高度埋藏位点但属于较保守替换；NAAA Leu325Phe预测去稳定但AlphaMissense低；CNR1 Ile227Val为温和替换；NAPEPLD Arg153His位于二聚体界面，方法间不一致；TRPV1两个变异影响较小；DAGLA Ala750位于低置信度无序区域，所有19种替换均预测中性。各位置耐受性差异显著。
3.4 携带者临床特征：FAAH p.(Arg326Pro)携带者为21岁女性，BMI 45.0 kg/m²，合并代谢综合征和多囊卵巢综合征。

讨论与结论：既往研究集中于ECS常见变异，本研究发现罕见变异并提示FAAH p.(Arg326Pro)可能通过破坏蛋白稳定性影响功能，但未进行酶活测定，不能确定因果关系。讨论指出样本量小、缺少功能验证、起病年龄为分类变量等局限；未来需生化实验、动物模型和更大队列验证。研究结论为：在主要为儿童期起病的严重肥胖成人中鉴定出ECS关键基因的罕见变异；综合结构分析将FAAH p.(Arg326Pro)列为最可能影响蛋白稳定性的候选；其余变异预测不一；DAGLA第750位一致耐受。这些发现是描述性的，计算证据仅用于优先排序，不能单独证明致病性或表型贡献，需进一步功能研究。