儿童先天性心脏病与神经发育障碍共享遗传通路的多组学研究

《HUMAN MUTATION》：Multiomic Investigation of Shared Genetic Pathways in Paediatric Congenital Heart Disease and Neurodevelopmental Disorders

　　

近期医学进步显著提高了先天性心脏病（CHD）患者的预期寿命；然而，这些儿童仍面临共患神经发育障碍（NDD）的高风险，例如注意缺陷/多动障碍和自闭症谱系障碍。尽管产前环境因素（包括胎盘功能障碍和宫内氧水平改变）以及产后事件（如心脏手术）可能促成这一风险增加，但共享遗传因素也可能同时构成这两种疾病的基础。因此，本研究旨在应用整合多组学方法，探究NDD和/或CHD患者中的基因组、转录组和表观基因组发现。研究人员招募了一个由14个三人核心家系和一个两人核心家系组成的队列：两名先证者同时患有NDD和CHD，两名仅患有CHD，11名仅患有NDD。研究人员对血液样本进行了全基因组测序（WGS）、RNA测序（RNA-seq）和DNA甲基化谱分析。研究人员鉴定出一个大片段缺失（约2.5 Mb）和七个可能致病/致病性（LP/P）变异，其中两个位于与患者表型相关的常染色体显性基因中，五个位于与携带者状态一致的常染色体隐性基因中。这包括染色质重塑基因ARID1B中的一个新生可能致病性剪接破坏变异，该变异通过RNA测序验证。DNA甲基化分析揭示了CHD与NDD患者之间的表观遗传差异，CHD患者显示出升高的生物学年龄加速。与一个包含178名对照的参考队列比较，研究人员识别出四名具有极端甲基化失调的先证者，包括携带ARID1B变异的个体。总体而言，这些发现强调了整合多组学分析在儿科发育疾病队列中的诊断和机制价值。

**论文解读：儿童先天性心脏病与神经发育障碍共享遗传通路的多组学研究**

**研究背景与问题**

先天性心脏病（CHD）是常见的出生缺陷，随着医疗进步，患者生存率显著提高，但共患神经发育障碍（NDD）的风险仍然很高，估计约35%的CHD患者伴有某种形式的NDD，而普通人群仅为3%–11%。NDD包括一系列影响认知、运动、情绪和行为发育的疾病，如注意缺陷/多动障碍（ADHD）和自闭症谱系障碍（ASD）。已知的风险因素包括心血管异常严重程度、新生儿心脏手术需求、住院时间延长及家庭社会经济劣势等。此外，严重CHD胎儿在宫内可能经历脑氧合和营养供应减少，导致脑发育受损；新生儿期围手术期暴露也可能对发育中的大脑施加额外生理应激。然而，这些因素仅能解释约45%的神经发育结局变异，因此研究人员假设遗传因素可能在CHD与NDD的共病中起重要作用。近年研究已开始探索两者间的遗传重叠，发现一些关键基因的变异同时影响心脏和神经发育，但这些分子通路尚不完全清楚。为此，本研究旨在通过整合多组学方法，在NDD和/或CHD先证者中鉴定发育相关基因的遗传变异，并利用RNA测序和DNA甲基化数据描绘与NDD和CHD相关的基因表达及甲基化变化。

**研究内容与结论**

研究人员在美国辛辛那提儿童医院医学中心招募了14个三人核心家系和1个两人核心家系，共15名先证者（2名同时患有NDD和CHD，2名仅患CHD，11名仅患NDD）。对先证者及其生物学父母的外周血样本进行全基因组测序（WGS，≥30×覆盖度）、RNA测序（针对9个三人核心家系和1个两人核心家系，≥2500万对端读数）和DNA甲基化谱分析（使用Illumina Infinium MethylationEPIC BeadChip芯片）。通过GATK最佳实践流程进行变异检测，使用ANNOVAR和Ensembl Variant Effect Predictor注释，并根据ACMG指南分类。对欧洲血统先证者计算多基因风险评分（PRS），使用OUTRIDER鉴定异常表达基因，通过交叉队列差异甲基化分析识别甲基化失调个体。主要结论：共鉴定出1个新生大片段缺失（约2.5 Mb）和7个LP/P变异，包括染色质重塑基因ARID1B的新生剪接位点变异（通过RNA-seq验证致异常剪接），以及MYH6、RERE等基因的变异。DNA甲基化分析显示，CHD患者（特别是左心发育不全综合征HLHS患者）表现出显著的生物学年龄加速（较NDD-only患者平均高9.5岁，较对照平均高7.4岁）。4名先证者存在极端甲基化失调（>10,000个差异甲基化位点），其中携带ARID1B变异的先证者甲基化模式尤为广泛。通路富集分析显示差异甲基化区域富集于肌肉发育、心脏收缩调节和神经元导向等通路。该研究发表在《HUMAN MUTATION》上，强调了整合多组学分析在儿科发育疾病队列中的诊断和机制价值。

**主要关键技术方法**

1. **全基因组测序（WGS）**：对13个三人核心家系和1个两人核心家系进行≥30×覆盖度的Illumina测序，通过GATK最佳实践流程进行变异检测和注释，使用ANNOVAR和Ensembl VEP，并根据ACMG标准分类。2. **RNA测序分析**：对9个三人核心家系和1个两人核心家系进行RNA-seq，使用OUTRIDER鉴定异常表达基因，并通过剪接连接计数验证变异导致的剪接缺陷。3. **DNA甲基化谱分析**：使用Illumina Infinium MethylationEPIC BeadChip芯片对14个三人核心家系和1个两人核心家系进行甲基化检测，通过Hannum和Horvath时钟估计表观遗传年龄，并与两个GEO参考对照数据集（GSE74432, n=53；GSE97362, n=125）整合，使用limma鉴定差异甲基化位点（DMP），以FDR<0.05筛选。4. **多基因风险评分（PRS）分析**：对9名欧洲血统先证者计算ADHD、焦虑、自闭症和抑郁的PRS，与健康老年欧洲对照队列（n=3823）比较。

**研究结果**

**3.1 临床信息**：队列包含10名男性先证者（66.7%）和5名女性。4名先证者诊断为左心发育不全综合征（HLHS），其中2名也患有NDD。最常见的NDD类型为ADHD（13/15）。8名先证者至少有一名父母患有NDD，其中6名有更广泛的NDD家族史。

**3.2 WGS变异**：共鉴定出33个稀有潜在有害变异和1个大片段缺失。根据ACMG标准，3个变异/缺失被认为对先证者表型有潜在影响（常染色体显性或新生作用基因），包括：先证者9携带一个新生2.5 Mb缺失（涉及DVL3、ALG3、CLCN2、ECE2）；先证者16携带父源MYH6无义变异（可能与心脏表型相关）和父源RERE无义变异（意义未明）；先证者19携带新生ARID1B剪接位点变异（经RNA-seq验证）。另外5名先证者携带常染色体隐性致病基因的单杂合LP/P变异（携带者状态），包括MYO7A、TWIST2、GPR179、COL18A1和ACADM。全基因组功能丧失分析在CHD9、CNTFR、INTS6、SPTBN1等基因中发现6个新变异。

**3.3 多基因风险评分分析**：队列中ADHD和ASD的PRS呈升高趋势，但未达统计学显著性。个体水平观察支持多基因遗传在NDD风险中的作用，例如先证者4（仅NDD）无单基因LP/P变异，但ADHD PRS在第89百分位；先证者16（NDD+CHD）携带MYH6变异，同时焦虑PRS在第88百分位。

**3.4 基因表达**：3名先证者中14个基因异常表达，其中8个基因（EIF4G1、EIF2B5、DVL3等）在先证者9中下调，与染色体3大片段缺失一致。先证者8携带一个PPID基因的新生意义未明变异，可能解释其异常表达。

**3.5 RNA测序验证致病性剪接缺陷**：对3个预测破坏剪接位点的变异进行验证，其中ARID1B变异（先证者19）导致外显子跳跃（父母中均无），COL18A1变异（先证者15）导致隐蔽剪接供体激活和内含子保留，均引入提前终止密码子。

**3.6 单等位基因表达分析**：整合DNA甲基化和RNA-seq数据，在10个三人核心家系中评估单等位基因表达模式。总体MAE率为0.3%（656个位点），无家族超过1%，未发现广泛等位基因特异性沉默或印记缺陷。

**3.7 表观遗传年龄加速**：Hannum和Horvath时钟均显示HLHS患者生物学年龄加速显著大于NDD-only组（平均高9.5年，p=0.006）和对照组（平均高7.4年，p=0.008）。先证者总体年龄加速高于其父母（先证者平均11.1年，母亲6.8年，父亲6.6年）。EN和Levine时钟也支持CHD先证者表观遗传年龄升高，但未达显著性。

**3.8 差异DNA甲基化分析**：与178名对照比较，4名先证者（13、15、19、20）显示极端甲基化失调（>10,000个DMP）。先证者13有72,355个DMP，以CpG岛高甲基化和开放海低甲基化为特征；先证者19（56,816 DMP）和20（49,618 DMP）甲基化模式相似，均表现为广泛低甲基化；先证者15（10,037 DMP）表现为全基因组高甲基化。

**3.9 甲基化通路富集分析**：差异甲基化区域富集于肌肉组织发育、肌肉细胞分化、心脏收缩调节、轴突导向、突触组织和神经冲动传导等通路，提示心脏和神经发育的共享表观遗传调控。

**讨论与结论总结**

讨论部分指出，该研究通过整合多组学框架，在小型队列中展示了基因组、转录组和表观基因组数据在个体水平提供互补诊断和机制信息的能力。最显著的发现是ARID1B新生剪接变异及其导致的广泛甲基化失调，该基因是SWI/SNF染色质重塑复合物组分，其功能丧失可产生全基因组表观遗传后果。其他极端甲基化失调的先证者显示出不同模式，提示存在异质性分子亚群。HLHS患者的显著生物学年龄加速可能反映慢性缺氧和手术应激，但需更大队列验证。研究局限性包括样本量小、仅两名共病先证者、缺乏匹配对照基因组分析等。结论：整合多组学方法识别出与CHD和NDD相关的分子发现，包括影响染色质重塑的变异和异质性表观遗传失调模式，支持在更大队列中进一步研究共享发育机制。通路富集分析提示肌肉发育和神经元导向之间的可能分子联系，与早期胚胎发生中共享祖细胞群的单细胞转录组研究一致。这些发现强调了整合基因组和表观基因组分析在捕捉疾病复杂性和提高诊断准确性方面的价值。

**研究结论翻译**：研究人员的整合多组学方法识别出与CHD和NDD相关的分子发现，包括影响染色质重塑的变异和异质性表观遗传失调模式。这些观察结果支持在更大队列中进一步研究共享发育机制。通路富集分析，鉴定出肌肉发育和神经元导向的参与，提示心脏和神经发育之间可能存在分子联系。这与最近的单细胞转录组研究一致，该研究指出早期胚胎发生过程中存在一些共享的祖细胞群。总体而言，这些发现强调了整合基因组和表观基因组分析对于捕捉这些疾病的复杂性以及提高受影响儿童诊断精准度的价值。需要更大、表型分层的队列和功能验证研究来阐明因果关系并确定这些观察结果的完整临床意义。