《Emerging Microbes & Infections》：HLA profiles and susceptibility to Alveolar echinococcosis: evidence from a large cohort with population and household controls
编辑推荐：
泡型棘球蚴病（alveolar echinococcosis, AE）是一种新出现且可能致死的蠕虫病。尽管在地方性流行区感染性多房棘球绦虫（Echinococcus multilocularis, EM）虫卵广泛污染环境，人AE仍罕见，提示宿主遗传易感性具有重
泡型棘球蚴病（alveolar echinococcosis, AE）是一种新出现且可能致死的蠕虫病。尽管在地方性流行区感染性多房棘球绦虫（Echinococcus multilocularis, EM）虫卵广泛污染环境，人AE仍罕见，提示宿主遗传易感性具有重要作用。在这项病例–对照研究中，研究人员通过基于扩增子的二代测序（amplicon-based next-generation sequencing）评估360例AE患者的HLA变异，并与（i）1,440名人群对照进行比较，对照在可行时按年龄、性别和居住地匹配，以及（ii）100名无关家庭对照，以考虑对EM的共享环境暴露。分析包括年龄和性别校正的Firth惩罚logistic回归、考虑连锁不平衡（linkage disequilibrium, LD）的单倍型分析、基于模拟的最大评分检验以及Benjamini–Hochberg错误发现率（false discovery rate, FDR）校正。在双字段分辨率下，C*07:01与AE odds增加相关（OR 2.10；校正后p=0.006），而B*18:01与odds降低相关（OR 0.49；校正后p=0.014）。单倍型分析识别出一个由LD定义的扩展背景，涉及A*29、C*16或两者，常与B*44、DRB1*07或DQB1*02组合。A*29携带对应A01–A24超型类别。A*29–DQB1*02在模型内和全研究校正后仍显著（OR 3.92；校正后p分别为0.005和0.039），而C*16–B*44–DQB1*02在模型内校正后仍显著（OR 2.83；校正后p=0.041），但在全研究校正后不显著。对早发AE（≤35岁；n=63）的探索性分析提示与A*31:01（OR 2.86）和DQB1*05:03（OR 3.44）相关，尽管二者在多重检验校正后均不显著。这些发现支持个体HLA等位基因和扩展单倍型背景在多因素宿主–环境模型中对AE易感性的贡献。
泡型棘球蚴病（alveolar echinococcosis, AE）是由多房棘球绦虫（Echinococcus multilocularis, EM）引起的慢性、潜在致死性蠕虫病，在北半球部分地区包括德国南部流行。成虫寄生于红狐（Vulpes vulpes）肠道，随粪便排出感染性虫卵；人意外摄入后，幼虫阶段可在肝脏建立并呈肿瘤样侵袭性生长。分子流行病学研究显示，地方性流行区厨房花园和食品可被EM DNA污染，但人AE仍罕见（0.1–2.8/10万），提示除暴露外，宿主因素包括遗传易感性参与感染和疾病发生。人类白细胞抗原（human leucocyte antigen, HLA）系统是适应性免疫中心：HLA I类分子向CD8
+ T细胞提呈内源性肽，HLA II类分子向CD4
+ T细胞提呈外源性肽；HLA多态性影响抗原提呈和个体间免疫应答差异，且HLA位点间存在强连锁不平衡（linkage disequilibrium, LD），关联可能涉及遗传等位基因组合而非单个等位基因。既往在瑞士、法国、德国、中国、土耳其、俄罗斯、阿拉斯加和拉脱维亚等的研究结果异质甚至冲突；解释受限因素包括样本量小、分型分辨率和区域HLA背景差异、位点覆盖不全、单倍型背景考虑不足、缺乏暴露匹配对照。因此，研究人员在大型德国队列中评估HLA等位基因组、等位基因和单倍型与AE易感或抵抗的关联，并与人群对照和无关家庭对照联合及分别比较；家庭对照用于更好考虑共享EM环境暴露。还探索HLA纯合性、I类超型和早发AE（≤35岁）。由于欧洲AE平均诊断约54岁，早发可能代表宿主遗传易感性更突出的表型。稳健HLA信号若复制，可作为宿主–暴露风险模型的一部分，为更有针对性监测和精准公共卫生提供信息。
关键技术方法：研究为Ulm大学医院观察性病例对照研究。患者为2024年3月至2025年5月AE门诊成人可能或确诊AE（WHO-IWGE定义）。人群对照来自健康志愿干细胞捐赠者登记，按性别、年龄和诊断前10年居住地匹配，每例患者选4名；家庭对照为与患者无生物学亲缘、诊断前同住≥5年的家庭成员，经两种血清学检测和腹部超声排除AE。关键技术为基于扩增子的二代测序（next-generation sequencing, NGS）对五个HLA位点进行高分辨率分型；统计采用年龄和性别校正的Firth惩罚logistic回归、LD感知单倍型分析（期望最大化算法、加性单倍型回归和单倍型评分/最大评分模拟）以及Benjamini–Hochberg错误发现率（false discovery rate, FDR）校正。
结果部分保留原文小标题并浓缩如下。
参与者特征：比较360例AE患者与1540名对照，包括100名无关家庭对照和1440名人群对照。患者平均年龄63.8岁，家庭对照64.6岁，人群对照59.1岁；患者与家庭对照年龄无显著差异，人群对照显著年轻。性别分布无差异，患者58%女性。
等位基因组和等位基因：单变量显示C*07:01患者更频（35.8% vs 26.8%；OR 1.52，raw p<0.001）；A*29/A*29:02、B*08/B*08:01、B*38/B*38:01、C*16/C*16:01呈名义风险，B*18/B*18:01保护；II类DRB1*14/DRB1*14:01、DQB1*05:03风险，DRB1*01/DRB1*01:01、DQB1*05:01保护。家庭对照比较中DRB1*01:01保护（10.0% vs 19.0%；OR 0.47，raw p=0.014），是唯一在两种对照比较中raw p<0.05者；A*29:02和C*16同方向接近名义。未通过全单变量BH-FDR。多变量一字段模型中B*08、B*38、DRB1*14名义风险，A*66名义保护但AE患者中未观察，稀疏；均未通过模型内BH。二字段模型B*18:01降低约50% odds（OR 0.49，raw p=0.003，BH p=0.014），C*07:01约两倍（OR 2.10，raw p<0.001，BH p=0.006）；A*66:01和B*38:01同方向但未通过BH。B*08:01与C*07:01中度依赖，条件分析提示共享LD信号，不能区分独立因果。人群对照模型相似；家庭对照无显著，C*17和A*02:01名义但未过BH。
纯合性：二字段总体HLA纯合性与AE无关，平均纯合位点0.44 vs 0.47，校正OR 0.95，p=0.476；无位点特异关联。
单倍型分析：成对分析24个名义关联；A*29–DQB1*02最强，频率3.13% vs 1.00%，OR 3.92，raw p<0.001，是唯一在模型内和全研究回归BH后显著（p=0.005和0.039）。最大评分分析10个位点对中5个经BH显著；最强正评分包括A*29–B*44、A*29–C*16、A*29–DRB1*07、A*29–DQB1*02、C*16–DRB1*07；模型级最大评分p 0.0129–0.0227，BH后0.0454。回归中A*29–C*16、A*29–DRB1*07、C*16–DRB1*07 OR 1.98–2.46；A*29–B*44未达名义回归。三字段分析14个名义；仅C*16–B*44–DQB1*02通过模型内BH（3.31% vs 1.36%，OR 2.83，raw p=0.0016，模型内BH p=0.041），全研究BH无。复发风险模式为A*29、C*16或两者与B*44、DRB1*07、DQB1*02组合；7个相关单倍型名义OR 2.29–2.83。评分分析中9/10位点组合最强正评分来自该扩展背景，最大评分BH后p=0.0015–0.0495。DRB1*07呈环境依赖：在含A*29或C*16风险单倍型中出现，与B*57组合名义保护。
HLA I类超型分析：探索性年龄性别校正中，A01–A24交叉反应模式唯一名义关联（OR 1.84，raw p=0.019），BH校正12个后不显著（p=0.233）。该模式包括A*29:01、A*29:02、A*29:03，对应A*29携带组，与等位/单倍型信号一致，但无独立显著超型效应。
早发AE探索性分析：63例≤35岁患者。一字段A*31（OR 3.39）、DRB1*14（OR 2.95）、B*07（OR 0.48）名义，经选定模型内BH显著但未跨全部78个一字段变量。二字段A*31:01（OR 2.86）和DQB1*05:03（OR 3.44）名义，未通过模型内或全家族校正。DQB1*03:02和DRB1*04:01单变量名义，联合时相互影响，提示重叠LD。Leave-one-case-out示DQB1*05:03较A*31:01更稳健。
讨论部分：关键发现与解释：单标记C*07:01最强名义风险；DRB1*01:01唯一在总体、人群和家庭比较中名义保护；A*29:02和C*16人群对照风险，家庭对照同方向但未达名义；均未通过全单变量BH。多变量中B*08、B*38、DRB1*14仍名义，C*07:01约两倍，B*18:01约半；B*08:01/C*07:01相关，不能归因独立因果。单倍型A*29–DQB1*02最稳健，C*16–B*44–DQB1*02唯一三字段通过模型内BH。最大评分支持A*29、C*16或两者与B*44、DRB1*07、DQB1*02组合的扩展风险背景。A*29多变量不显著，不能确定直接或连锁。超型和早发未通过完整校正；HLA纯合性无关。观察关联可能反映HLA依赖的EM来源肽提呈差异，影响T细胞应答，造成AE易感个体差异。
与国际文献整合：既往结果异质，可比性受人群、分辨率、位点、样本量限制。Aydinli et al.报道HLA-A*26在44例土耳其患者更频；本队列A*26患者8.3% vs 对照6.0%，不显著。DRB1*15和DQB1*02、*06、*07在土耳其对照更频，本队列无显著。DRB1*11：Eiermann et al.报告151例欧洲AE患者17.2% vs 453对照27.4%，未校正显著；本队列vs人群对照17.2% vs 22.1%（raw p=0.043），vs家庭19.0%（p=0.680）。Laivacuma et al.报告DRB1*11单倍型风险但仅8例。DRB1*07环境依赖，Li et al.35例中国患者DRB1*07:01较低，Laivacuma et al.风险方向；两项均未含I类位点确定背景。Li et al. DRB1*040x易感；本队列DRB1*04:01患者15.8% vs 12.3%（raw p=0.070），早发23.8% vs 12.3%，条件模型去除DQB1*03:02后名义，C*04–DRB1*04–DQB1*03单倍型名义风险。中国中部诊断平均年龄40岁 vs 欧洲53岁，为潜在平行。Laivacuma et al. DQB1*05:01较高；本队列DRB1*14–DQB1*05单倍型名义风险（OR 1.70，raw p=0.043）。早发DQB1*05:03 7/63（11.1% vs 4.6%，校正OR 3.44，raw p=0.014），leave-one-out名义显著；因两个字段等位基因不同，属等位基因组相似而非直接复制。俄罗斯早期血清学HLA-A1信号重叠有限；A*01未单独关联，但在A*01–DQB1*02名义单倍型和B*08、C*07、DRB1*03频繁背景中。其他报道A11、B35、DRB1*09:01、DRB1*16:01/02未支持。相反，A*29/C*16关联背景未在既往突出，因既往未做正式多座位单倍型分析；人群特异LD、样本量和位点覆盖有限也可解释。
结论与意义（翻译）：在多变量模型中，C*07:01和B*18:01分别显示最强校正风险和保护信号，尽管C*07:01信号出现在涉及B*08:01的相关LD背景中。单倍型分析识别出扩展A*29/C*16相关背景，其中A*29–DQB1*02在全研究校正后仍显著。A01–A24超型模式也反映A*29携带。环境依赖的DRB1*07关联可能有助于解释Li et al.和Laivacuma et al.报道的矛盾结果，强调考虑单倍型背景的重要性。Eiermann et al.和Aydinli et al.报道的DRB1*11保护倾向在人群对照中观察到，但在家庭对照中未观察到，说明考虑环境暴露的价值。探索性早发分析提示与A*31:01和DQB1*05:03的关联，而DRB1*04:01信号与Li et al.报道有潜在平行。这些关联不支持个体风险预测或基于HLA的筛查。相反，AE易感性可能反映遗传背景、环境暴露和免疫状态之间的相互作用。未来研究应优先在独立队列中复制、使用暴露可比对照、进行LD感知单倍型分析和早发疾病研究。
优势与局限：该研究为迄今评估AE中HLA关联的最大队列，超过既往已发表研究AE患者样本量总和。优势包括五个HLA位点高分辨率分型、年龄性别居住地匹配并校正、多重检验校正、考虑LD的互补单倍型方法；无关家庭对照提供共享环境暴露视角。局限包括家庭对照数量少降低统计效能，家庭与人群对照差异可能反映随机抽样或人群构成差异；未正式评估远亲或未识别亲缘；人群对照较年轻，虽中位年龄差小并校正年龄，残余混杂不能排除；超型基于已发表肽结合分类，无实验验证EM肽与特定HLA结合，关联应视为功能研究候选而非等位特异抗原提呈证据；单中心设计限制外推，需国际合作多中心研究。
