APOE位点的单倍型解析DNA甲基化识别与阿尔茨海默病风险相关的等位基因特异性表观遗传特征

《npj Dementia》：Haplotype-resolved DNA methylation at the APOE locus identifies allele-specific epigenetic signatures relevant to Alzheimer’s disease risk

　　

APOE基因编码一种脂质转运蛋白，是阿尔茨海默病（Alzheimer's disease, AD）发病机制的核心。三种常见等位基因——ε₂（rs7412(C?>?T)）、ε₃（参考序列）和ε₄（rs429358(T?>?C)）——由外显子4中的两个编码变异产生，并赋予不同的AD风险谱，其中ε₄增加风险，ε₂具有保护作用。与ε₃连锁的APOE变异rs769455[T]也与携带APOE ε₄等位基因的非洲裔个体的AD风险增加相关。确定遗传变异如何影响CpG甲基化需要methQTL型分析，但传统的亚硫酸氢盐和基于芯片的方法在区分等位基因特异性效应方面的分辨率有限。在这里，研究人员使用高精度长读长测序（long-read sequencing, LRS），在来自祖先多样性队列的332个死后脑组织样本中生成跨APOE位点的单倍型解析甲基化谱，包括201个欧洲裔个体样本和131个非洲裔及非洲混血个体样本。将每个单倍型作为一个独立观察值，普通最小二乘（ordinary least squares, OLS）回归识别了18个与ε₂、ε₄和rs769455[T]相关的跨APOE位点（TOMM40、APOE、APOC1和APOC4-APOC2基因）的新型差异甲基化CpG位点（differentially methylated CpG sites, DMCs）。这些发现揭示了独特的等位基因特异性甲基化特征，并证明了长读长测序在解析与AD风险相关的表观遗传变异方面的实用性。

**论文解读：APOE位点等位基因特异性甲基化与阿尔茨海默病风险**

**研究背景与问题**

载脂蛋白E（Apolipoprotein E, APOE）基因位于染色体19q13.32，编码一种脂质转运蛋白，对神经元发育、维持和修复至关重要。APOE是晚发型阿尔茨海默病（late-onset Alzheimer's disease, LOAD）最强的常见遗传风险因子。其外显子4中的两个单核苷酸变异（single nucleotide variants, SNVs）——rs7412(C>T)和rs429358(T>C)——定义了三个常见等位基因：ε₂（保护性）、ε₃（参考）和ε₄（风险性）。此外，在非洲裔人群中，与ε₃连锁的错义变异rs769455(C>T)也被发现与AD风险增加相关。除了遗传变异，表观遗传修饰如DNA甲基化（5-甲基胞嘧啶，5-methylcytosine, 5mC）也参与APOE表达调控和AD发病。以往研究主要使用短读长测序（short-read sequencing, SRS）或基于芯片的方法（如Illumina Infinium HumanMethylation450 BeadChip）检测APOE基因簇（包括TOMM40、APOE、APOC1和APOC4-APOC2）中的甲基化，但这些技术存在局限性：亚硫酸氢盐转化导致DNA降解、芯片覆盖CpG位点不全、SRS无法有效分相等，因此无法区分等位基因特异性甲基化。基于此，研究人员利用高精度的牛津纳米孔技术（Oxford Nanopore Technologies, ONT）长读长测序（long-read sequencing, LRS），直接检测单碱基分辨率的核苷酸修饰，并实现单倍型解析的甲基化分析，以解决上述问题。

**研究内容与核心发现**

研究人员利用美国国立卫生研究院阿尔茨海默病及相关痴呆症中心（NIH Center for Alzheimer's and Related Dementias, CARD）长读长测序计划公开的ONT全基因组测序数据，分析了332个死后额叶皮层脑组织样本。这些样本来自两个队列：北美脑表达联盟（North American Brain Expression Consortium, NABEC，n=201，欧洲裔）和人类脑采集核心（Human Brain Collection Core, HBCC，n=131，非洲裔及非洲混血）。通过单倍型分型和CpG甲基化调用，以每个单倍型为独立观察值，采用OLS线性回归模型，鉴定与APOE ε₂、ε₄和rs769455[T]等位基因显著相关的差异甲基化CpG位点。研究共发现18个新型DMCs，分布在APOE基因簇的多个调控区域，揭示了独特的等位基因特异性甲基化特征。这些发现发表在《npj Dementia》上，强调了LRS在解析与AD风险相关的复杂表观遗传变异中的关键作用。

**关键技术方法（≤250字）**

研究人员使用了ONT PromethION平台的LRS数据（NABEC样本采用R9.4.1化学，HBCC样本采用R10.4.1化学），从332个死后脑组织样本（来源：NABEC队列201例欧洲裔，HBCC队列131例非洲裔/非洲混血）中提取DNA，经片段化后测序。利用PEPPER-Margin-DeepVariant工具进行单倍型分型，通过modkit进行CpG甲基化位点检测和频率计算。针对每个CpG位点，采用OLS线性回归模型，以等位基因携带状态（0/1）为自变量，甲基化频率为因变量，校正年龄、性别、死后间隔、脑库来源、遗传主成分和甲基化主成分。同时，使用未分相的甲基化数据进行了等位基因平均分析，并利用样本对应的单核RNA测序数据估算细胞类型比例，评估其对甲基化的影响。

**研究结果**

**1. Comprehensive resolution of haplotypes and allele-specific methylation through long-read sequencing**

通过比较ONT LRS与Illumina 450k芯片在APOE基因簇中的CpG检测能力，发现LRS在NABEC样本中检测到1556个CpG位点，而450k芯片仅检测到46个，LRS检测数量约为芯片的33倍。两种平台共有的46个位点（对应于850k芯片的76个位点）的甲基化频率高度相关（Pearson r=0.98），验证了LRS数据的可靠性。LRS的覆盖范围扩展到内含子、基因间区及调控区域，突破了芯片仅覆盖启动子和外显子的限制。

**2. Haplotype-resolved long-read sequencing refines phylogenetic structure**

利用FastTree基于近似最大似然法构建系统发育树。未分相数据仅显示基因型水平的分离，等位基因聚类不连续且部分重叠；而单倍型解析数据则形成了对应于ε₂、ε₃、ε₄和rs769455[T]四个清晰分支的聚类，表明单倍型分相显著提升了APOE位点遗传结构的解析能力。此外，使用同一份样本分别用R9和R10化学测序的数据验证，确认聚类差异来自遗传序列而非测序化学批次效应。

**3. Allele-specific differentially methylated CpG sites in the APOE gene cluster region**

在NABEC队列中，鉴定出8个显著的DMCs（BH-FDR校正P<0.05）。其中，3个与ε₂相关（如cpg_chr19_44914329位于未表征单外显子基因ENSG00000280087.1），6个与ε₄相关（如cpg_chr19_44917997位于APOC1内含子3的短散在核元件SINE中，其甲基化降低与ε₄高度相关）。最显著位点cpg_chr19_44917997与APOC1内含子变异rs12721046(G>A)连锁，该变异已知在非西班牙裔白人中增加AD风险。在HBCC队列中，鉴定出11个DMCs，其中1个与ε₄相关（cpg_chr19_44896082位于APOE内含子2附近的CpG岛），10个与rs769455[T]相关（如cpg_chr19_44904817显示低甲基化，位于APOE与TOMM40之间）。所有18个DMCs均为新发现。与之对比，未分相等位基因平均分析仅在NABEC中检出3个DMCs且显著性降低，在HBCC中未检出任何DMC，证明单倍型分相是检测等位基因特异性甲基化的关键。

**讨论与结论**

该研究是迄今利用LRS对人类脑组织中APOE位点进行等位基因特异性甲基化分析最全面的一项。18个DMCs富集于启动子、CpG岛、增强子、内含子和SINE等潜在调控区域，为AD风险的表观遗传机制提供了精细图谱。值得注意的是，之前唯一与芯片探针重叠的位点cpg_chr19_44914329在血液和脑中呈现相反的甲基化方向（血液中低甲基化、脑中高甲基化），提示组织特异性差异。ε₂特异性甲基化模式在人类脑中首次报道，为理解其保护作用提供了表观遗传线索。研究局限性包括：不同测序化学可能引入技术变异；非洲裔样本量较小，且未检测到跨队列共享DMCs，可能与遗传背景和统计效力有关；表达关联分析（eQTL和eQTM）未发现显著结果，可能受限于短读长RNA-seq的分辨率；细胞类型比例分析未显示显著影响，但单核数据来自邻近组织而非精确同一区域。未来需要结合长读长RNA-seq、单细胞分辨率和功能实验（如等位基因特异性报告实验、CRISPR表观基因组编辑）来验证甲基化变化的调控后果。

**研究结论翻译：** 总之，本研究强调了长读长测序在解析人类脑组织中APOE位点复杂单倍型结构和表观遗传景观方面的优势。这些发现为APOE遗传变异如何塑造甲基化模式提供了新见解，并为剖析连接APOE单倍型与AD相关细胞功能障碍的分子机制奠定了基础。