整合基因组结构方程模型（Genomic-SEM）与实验验证揭示男性生殖器硬化性苔藓的遗传基础

《Frontiers in Immunology》：Integrating genomic structural equation modeling and experimental validation to unravel the genetic basis of male genital lichen sclerosus

　　

目的：通过整合基因组结构方程模型（Genomic-SEM）与多组学分析和实验验证，研究男性生殖器硬化性苔藓（MGLSc）的遗传结构，并鉴定与疾病发病机制相关的潜在易感位点、候选基因和生物学通路。方法：整合MGLSc相关表型的公开全基因组关联研究（GWAS）汇总统计数据，构建Genomic-SEM框架。使用连锁不平衡分数回归（LDSC）评估遗传相关性并评价模型稳定性。利用FUMA进行功能定位和注释，并通过GWAS减法策略进一步筛选新位点。采用SuSIE和FINEMAP进行精细定位，以优先候选因果变异。应用转录组全关联研究（TWAS）和FOCUS鉴定候选基因。进一步进行基于MAGMA的基因富集分析、分区遗传力分析和多基因风险评分（PRS）分析，以表征相关位点的生物学相关性。最后，进行体外RT-qPCR验证优先基因的表达。结果：Genomic-SEM显示出良好的整体拟合度，并生成了一个包含2,451,318个SNP的MGLSc间接GWAS框架。共208个SNP达到常规全基因组显著性，FUMA注释鉴定出43个风险位点、52个先导SNP和14个候选基因。通过GWAS减法策略进一步鉴定出13个新SNP，包括rs715299和rs10774625。精细定位突出显示了四个高置信度变异，即rs3134608、rs3134952、rs3763307和rs2076524，主要聚集在6号染色体主要组织相容性复合体区域。TWAS将HLA-DPA1鉴定为最显著基因，FOCUS进一步支持其可能的因果作用。MAGMA和富集分析提示，免疫相关和调控区域对MGLSc遗传力贡献显著。PRS分析显示染色体间存在显著异质性，其中6号染色体表现出最强的SNP水平贡献。RT-qPCR证实，与对照组相比，MGLSc样本中HLA-DPA1表达显著降低（P < 0.0001），与生物信息学预测一致。此外，MGLSc样本中多个MAGMA优先基因（包括C4B、DDR1、BBS7、VARS1、PRRT1、PPT2、EGFL8、BTNL2和EME1）下调，而AGPAT1和PBX2上调。结论：本研究通过整合Genomic-SEM、精细定位、转录组优先化和实验验证，系统揭示了MGLSc的遗传基础。研究结果表明，MGLSc受共享的多基因结构影响，该结构在免疫相关位点（尤其是HLA区域）富集。HLA-DPA1成为一个高置信度易感基因，并鉴定出多个新位点和候选基因，为MGLSc的分子机制提供了新见解，并为未来机制研究和转化研究提供了潜在靶点。

**研究背景与问题**
男性生殖器硬化性苔藓（MGLSc）是一种慢性苔藓样炎症性纤维化疾病，主要依据临床表现诊断，包括男性性功能障碍、瘙痒和皮肤病变。其病因尚不明确，既往假说涉及自身免疫机制、免疫失调和感染因素。包皮在MGLSc发展中起重要作用，出生时接受包皮环切术的男性中该病极为罕见。另一假说认为易感上皮细胞暴露于尿液可能导致MGLSc，但具体成分或特征尚未明确。作为一种罕见病，当前研究主要集中于临床病例报告，对遗传和生物学机制的探索几乎空白。为弥补MGLSc机制精确测量的缺失，研究人员设计了针对潜在未测量MGLSc的全基因组关联研究（GWAS），并通过整合基因组结构方程模型（Genomic-SEM）与多组学分析及实验验证，系统揭示其遗传基础。该研究发表在《Frontiers in Immunology》。

**主要关键技术与方法**
（1）基于公开GWAS汇总统计数据，从UK Biobank和FinnGen R12数据库获取MGLSc相关表型（冗余包皮、包茎/嵌顿、尿道狭窄、尿路感染、自身免疫性疾病）的GWAS数据，构建Genomic-SEM框架。
（2）使用连锁不平衡分数回归（LDSC）评估遗传相关性，通过FUMA进行功能注释，并采用GWAS减法策略筛选新位点。
（3）利用SuSIE和FINEMAP进行精细定位，结合转录组全关联研究（TWAS）和FOCUS鉴定候选基因。
（4）进行MAGMA基因富集、分区遗传力分析、多基因风险评分（PRS）分析，以及体外RT-qPCR验证优先基因表达。

**研究结果**
3.1 Genomic-SEM construction：通过LDSC分析显示自身免疫病遗传贡献最高，模型拟合优度良好（CFI=1.0，SRMR=0.144），生成含2,451,318个SNP的间接GWAS框架。
3.2 Hierarchical evaluation of genomic-SEM：在P<5×10^-8阈值下鉴定出208个显著SNP，P<5×10^-12和P<5×10^-16阈值下分别有150个和113个SNP，表明Genomic-SEM增强检测能力。
3.3 Assessment the stability of genomic-SEM via LDSC regression：通过参数筛选保留583,468个有效SNP，总观测尺度遗传力为1.9088×10^-9，证明观察到的膨胀源于多基因遗传信号而非混杂。
3.4 Evaluation of genomic-SEM using FUMA software：FUMA注释鉴定出43个风险遗传位点、52个先导SNP、14个候选基因；GWAS减法策略进一步识别13个新SNP，如rs715299和rs10774625，多数位于基因间区或内含子区，提示非编码调控元件的重要作用。
3.5 Finemap：精细定位确定四个高置信度变异（rs3134608、rs3134952、rs3763307、rs2076524），均位于6号染色体主要组织相容性复合体（MHC）区域，支持免疫相关位点主导MGLSc遗传关联。
3.6 Integrative transcriptome-wide association study and colocalization analysis：TWAS鉴定HLA-DPA1为最显著基因，FOCUS进一步支持其因果作用，且HLA-DPA1表达与MGLSc呈负相关（Z<0）。
3.7 Enrichment analysis of pathways, cell types, and Mendelian disease genes：MAGMA基因富集分析显示显著富集于溃疡性结肠炎、特发性膜性肾病、哮喘、自身免疫性甲状腺疾病等免疫相关疾病，MendelVar分析提示涉及内体组织、线粒体分裂、转录偶联核苷酸切除修复等生物学过程。
3.8 Heritability contributions across genomic regions：LDSC分析显示遗传力贡献主要富集于编码区和表达调控区，其中编码区和增强子区域效应最显著。
3.9 Polygenic risk scores results：PRS分析显示染色体间遗传贡献异质性显著，6号染色体SNP贡献最强，且与FINEMAP、MAGMA结果一致。
3.10 Sensitivity analysis results excluding autoimmune traits：排除自身免疫表型后模型拟合度下降（CFI=0.94，SRMR=0.269），但MHC信号效应大小和方向一致，HLA-DPA1关联仍可检测，提示局部泌尿生殖因素与系统性自身免疫背景共同塑造MGLSc遗传结构。
3.11 Validation of prognostic gene expression：RT-qPCR验证显示MGLSc样本中HLA-DPA1表达显著低于对照（P<0.0001），与生物信息学预测一致；MAGMA优先基因中C4B、DDR1、BBS7、VARS1、PRRT1、PPT2、EGFL8、BTNL2、EME1下调，AGPAT1和PBX2上调，支持免疫失调在MGLSc发病中的作用。

**讨论与结论**
讨论部分指出，MGLSc作为一种罕见慢性炎症性皮肤病，严重影响男性生活质量，且与阴茎鳞状细胞癌相关。通过Genomic-SEM分析揭示了四个单变量表型（冗余包皮、包茎/嵌顿、尿道狭窄、尿路感染、自身免疫病）之间的共享遗传因素，其中尿路感染和自身免疫病在MGLSc发病中起关键作用。新发现的SNP（如rs715299、rs10774625）表现出多效性，与免疫炎症、代谢、血液/泌尿系统等广泛表型相关，提示MGLSc遗传基础涉及免疫失调、自身免疫易感性和代谢紊乱的复杂网络。精细定位结果集中于MHC区域，进一步支持抗原呈递和免疫失调是核心机制。TWAS和FOCUS定位HLA-DPA1为高置信度易感基因，RT-qPCR验证其下调表达，同时C4B、BTNL2、DDR1等基因的差异表达从补体调控、T细胞免疫抑制和炎症纤维化重塑角度支持MGLSc的免疫介导性质。研究局限性包括样本主要来自欧洲人群、缺乏其他种族验证、男性特异性效应可能被稀释、实验验证样本量小等。结论：本研究通过整合Genomic-SEM、精细定位、转录组优先化和实验验证，系统揭示了MGLSc的遗传基础，表明MGLSc受共享的多基因结构影响，特别是免疫相关位点（HLA区域）富集。HLA-DPA1被确定为高置信度易感基因，多个新位点和候选基因的发现为MGLSc分子机制提供了新见解，并为未来机制研究和转化研究提供了潜在靶点。