《Frontiers in Psychiatry》：Increased frequency of rare damaging variants in PLIN4 in patients with anorexia nervosa
编辑推荐：
背景：神经性厌食症（anorexia nervosa, AN）是一种严重精神疾病，特征为对体重增加的病理性恐惧、持续限制能量摄入以及广泛的身体意象障碍。尽管 AN 具有显著遗传性，但促成易感性的罕见编码变异仍未完全明确。基于家系的全外显子组测序（whole-e
背景：神经性厌食症（anorexia nervosa, AN）是一种严重精神疾病，特征为对体重增加的病理性恐惧、持续限制能量摄入以及广泛的身体意象障碍。尽管 AN 具有显著遗传性，但促成易感性的罕见编码变异仍未完全明确。基于家系的全外显子组测序（whole-exome sequencing, WES）研究有助于识别罕见、高效应蛋白编码变异，从而阐明疾病机制。方法：对 AN 高密度家系进行 WES 以识别候选基因。在154名欧洲血统的无关先证者中，评估32个候选基因中罕见（次要等位基因频率<1%）、推测有害[组合注释依赖耗竭（combined annotation-dependent depletion, CADD）>15，包括预测功能丧失（loss-of-function, LOF）]变异的频率。携带者计数与 gnomAD v4 非芬兰欧洲人（non-Finnish European, NFE）参考进行比较，使用 TRAPD 风格覆盖度校正分母（每个基因合格位点最大等位基因数/2）。通过单侧 Fisher 精确检验评估显著性，在32个候选基因上作 Bonferroni 校正。结果：PLIN4 是唯一在 Bonferroni 校正显著性下 AN 病例中富集的基因：154例患者中17例（11.0%）携带罕见损伤性变异，而 gnomAD 非芬兰欧洲人为3.2%[比值比（odds ratio, OR）3.86，95%置信区间（confidence interval, CI）2.35至6.35；单侧 Fisher p = 9.98e?06；Bonferroni 校正 p = 3.19e?04]。移除索引先证者后结果保持（n = 153；OR 3.64，95% CI 2.18至6.07；校正 p = 1.13e?03）。154例患者中6例携带预测 LOF 变异；6例中4例至少有一名 AN 受累的一级或二级亲属携带相同变异，提供共分离证据。仅 LOF 负荷单独未通过 Bonferroni 校正。结论：AN 患者携带 PLIN4 中罕见、推测有害变异的统计学显著过量，该阈值在32个候选基因的 Bonferroni 校正后仍存在。这些发现提示脂滴生物学是 AN 病理生理学的候选贡献因素。
研究人员围绕神经性厌食症（anorexia nervosa, AN）的罕见编码变异易感性开展研究。AN 是严重精神疾病，特征为病理性恐惧体重增加、持续限制能量摄入和广泛身体意象障碍，可进展为危及生命的消瘦。部分心理治疗可改善青少年短期体重恢复，但相当比例患者从现有心理治疗或药物治疗中未获持久获益。AN 风险来自环境与遗传因素交互。全基因组关联研究（genome-wide association study, GWAS）识别常见、效应温和变异，但多位于非编码调控区，限制机制推断；家系分离研究分析多代家系，识别罕见、高效应变异，常位于蛋白编码区，功能后果更易表征。此前家系研究在 AN 多发病例家系中识别 ESRRA 和 HDAC4 罕见编码变异，并推动单基因小鼠模型。Walker et al. 报告迄今最大拷贝数变异（copy number variant, CNV）负荷分析，在7,414例 AN 中未发现罕见 CNV 显著富集，使罕见点突变负荷研究成为 AN 罕见变异生物学的主要发现范式。类似策略涉及脂代谢基因、脑源性神经营养因子（brain-derived neurotrophic factor, BDNF）和肽类神经递质基因。因此，全面理解 AN 遗传结构需要表征常见和罕见变异贡献。研究人员报告，PLIN4 中预测 LOF 和推测有害（CADD > 15）变异在 AN 患者中相对人群对照显著富集。PLIN4 编码 perilipin 4，是脂滴（lipid droplet, LD）包被蛋白 perilipin 家族成员。其他 perilipin 基因未见类似富集，提示 PLIN4 特异性贡献。
主要关键技术方法：研究队列来自 Precision Psychiatry (Plano, TX) 专科门诊转诊的154名无关 AN 先证者，主要为出生性别女性（148/154，96.1%）和非芬兰欧洲人（NFE）血统（150/154，97.4%）；PLIN4 队列为先前 BBOX1 队列（n = 183）完整子集，联合队列183例。对 AN 高密度家系行 WES，Novogene 100×，比对 GRCh38/hg38，GATK 联合基因分型，ANNOVAR 注释，CADD v1.6 致病性先验。家系发现44个共享罕见编码变异，32个 CADD > 15 候选。对32候选基因以 gnomAD v4 NFE 外显子为参考，采用 TRAPD 风格覆盖校正分母，单侧 Fisher 精确检验，Bonferroni 校正32个基因；并做仅 LOF 敏感性分析。
研究结果：
Index family pedigree：索引家系先证者患治疗抵抗 AN，母亲和两个姐妹均有 AN。选择该家系因一级和二级 AN 亲属密度最高。分析罕见（次要等位基因频率<0.01）编码变异在受累成员中共享，发现44个，32个 CADD > 15，其中2个预测 LOF。最高 CADD 变异映射到 SYNC、CDK15、MAS1、CSMD3、GRIK4、HAL、EAPP 和 PLIN4。虽然这些基因先前未与 AN 风险或摄食行为调节相关，但 PLIN4 的 rs8887 曾与膳食 omega-3 多不饱和脂肪酸敏感性增强相关，PLIN4 蛋白在细胞脂质储存和全身脂肪分布中有已知作用，故选择 PLIN4 随访。
PLIN4 variants in the unrelated AN cohort：筛查154名先证者的预测 LOF。rs780147097 为多等位；无其他先证者携带相同26碱基对插入，但另两名先证者在同染色体位置携带32碱基对插入。另有三名无关先证者各携带 PLIN4 不同预测 LOF；每例均有女性一级或二级亲属（分别为父系姑母、母亲、女儿）携带相同 LOF 等位基因并有确诊 AN。总计6/154（3.9%）无关 AN 先证者和6名女性亲属携带预测 PLIN4 LOF。随后筛查罕见（次要等位基因频率<0.01）、推测有害（CADD > 15）PLIN4 变异，另有11名先证者携带7个不同推测有害单核苷酸变异。合并频率17/154（11.0%）对比 gnomAD v4 NFE 3.2%，显著升高（OR 3.86，95% CI 2.35–6.35；单侧 Fisher p = 9.98e
?06）。其他 perilipin 家族成员未见显著富集。
Burden analysis of 32 genes with multiple-testing correction：为处理 PLIN4 事后选择，对32候选基因与 gnomAD v4 NFE 对照进行单侧 Fisher 负担测试，Bonferroni 校正32。仅 PLIN4 在富集方向达校正显著。n = 154：17例携带，OR 3.86，95% CI 2.35–6.35，p = 9.98e
?06，校正 p = 3.19e
?04。n = 153（移除索引先证者）：16例携带，OR 3.64，95% CI 2.18–6.07，p = 3.52e
?05，校正 p = 1.13e
?03。KIAA0754 是唯一另一达 Bonferroni 显著基因，但方向为缺失（n = 153：13/153，OR 0.43，95% CI 0.24–0.74，双侧 Fisher p = 7.88e
?04，校正 p = 0.025；n = 154：OR 0.42，双侧 p = 7.93e
?04，校正 p = 0.025），不符合 AN 风险升高，可能反映高频良性变异扩大对照分母。无其他基因达显著。
LOF-only sensitivity analysis：仅 LOF 负担单独报告。n = 154：6例携带，OR 4.56（95% CI 2.08–10.01），单侧 p = 3.8e
?03，校正 p = 0.12；n = 153：5例携带，OR 3.86（95% CI 1.65–9.04），单侧 p = 0.016，校正 p = 0.51。两种框定均未通过 Bonferroni。通过多重校正的统计信号是 CADD > 15 + LOF 联合负担，而非仅 LOF。然而 LOF 携带者提供共分离证据：n = 154 的6名无关 LOF 携带者中4名至少有一名女性亲属有 AN 或回避/限制性食物摄入障碍（avoidant/restrictive food intake disorder, ARFID）史并携带相同 PLIN4 LOF。四个共分离家系包括索引先证者母亲和两个姐妹（一级）携带 rs780147097；另一先证者父系姑母（二级）；一先证者母亲（一级）；一先证者女儿（一级）。家系水平和队列水平 LOF 计数不同，不互换。
EAPP co-carriage disclosure：索引先证者同时携带 PLIN4 和 EAPP LOF（rs370344686；p.Glu19Ter）。EAPP 负担在32候选 Bonferroni 校正后不显著（n = 154：1/154，OR 0.98；n = 153：0/153，OR 0.33），不解释为独立 AN 风险位点。PLIN4 发现在排除先证者后仍 Bonferroni 显著。
讨论部分，研究人员首先讨论显著性阈值和多重检验背景：负担 p 值未达全基因组 5e
?08，但本研究为家系锚定候选基因负担测试而非无假设全基因组扫描，适当标准是32候选基因 Bonferroni 校正，结果通过。LOF 携带者提供生物共分离证据，联合负担提供跨队列多重检验存活统计信号，两者互补，均必要但单独不足。生物学背景：研究人员此前报道 AN 患者 BBOX1 罕见推测有害变异富集；PLIN4 队列为 BBOX1 队列子集，联合183例中28/183（15.3%）携带 PLIN4 或 BBOX1 损伤变异（16仅 PLIN4，11仅 BBOX1，1两者），若独立证实提示脂肪酸储存和线粒体输入受损为共享机制。PLIN4 LOF 与体重、体脂百分比、腰臀比差异相关，携带者选择性减少 gynoid 脂肪并增加躯干脂肪；与 AN 成人恢复期躯干脂肪增加相关。LD 功能受损可解释 AN 体重恢复后高甘油三酯血症；预测 PLIN4 LOF 与血浆甘油三酯升高相关。本队列6名无关 LOF 携带者中5名有高甘油三酯血症或高脂血症临床史，余1名有淋巴细胞性结肠炎，可能减弱高脂血症表型。全部携带者中8/17（47.1%）有高脂血症临床诊断，非携带者11.7%（16/137），OR 6.72，95% CI 2.09–21.55，Fisher 双侧 p = 0.0011。携带者发病年龄较早（p = 0.047）和最低体重指数（body mass index, BMI）较低（p = 0.051）呈趋势，AN 亚型分布无差异；这些比较为描述性，未 Bonferroni 校正。PLIN4 主要表达于脂肪组织和骨骼肌，成人脑表达极低，但特定病理背景下有神经元表达报道；神经元 LD 来源脂肪酸可作为支持突触传递的代谢燃料储备，是否参与需直接实验。罕见变异图景：Walker et al. 未发现罕见 CNV 富集，Verebi et al. 发现 HCRTR1 剪接供体变异，与 BBOX1、PLIN4 等独立信号汇聚于食欲素信号、肉碱生物合成和 LD 生物学。局限性：单中心临床转诊、治疗抵抗 AN、以 NFE 为主，非人群样本；病例来自专科门诊，对照来自 gnomAD v4 NFE，存在确定偏倚，未达同队列病例/对照金标准；p 值未达全基因组显著但候选基因 Bonferroni 合适；其他推测有害变异仅基于 in silico 预测有争议；既往 GWAS 报告 PLIN4 与 AN 关联未达全基因组显著（p = 0.0012）；PLIN4 第3外显子含 99 碱基对基序的可变数目串联重复，需长读长测序；未来需确定变异如何影响蛋白功能。重复计划：优先在外部确定 AN 队列独立重复，已联系 Broad Institute 的 Genetics of Eating Disorders Study，计划招募10,000例进食障碍患者和10,000例对照，使用混合外显子基因组测序。临床意义：目前不推荐 AN 患者常规临床基因检测 PLIN4；发现需独立重复，单个 PLIN4 LOF 绝对风险增加未确定。若儿童偶然发现 PLIN4 LOF，当前最佳实践为家系级联检测、青春期高风险窗口行为和体重轨迹监测、遗传咨询强调风险增加但非确定性；不推荐无症状个体预测性检测。
结论翻译：AN 患者携带 PLIN4 中罕见、推测有害变异的统计学显著过量。联合 CADD > 15 + LOF 负担在索引家系筛选的所有32个候选基因 Bonferroni 校正后仍显著；仅 LOF 携带者在多个家系中提供独立共分离证据。这些发现提示脂肪组织和骨骼肌中的 LD 生物学是 AN 易感性的候选贡献因素，值得在外部确定队列中独立重复。
