Nature Medicine | 辅助生殖会增加孩子的新发突变吗？2.4万人WGS把答案拆到了ICSI、促排和胚胎培养每一步

先把尺度看清：一个孩子平均有多少个“父母都没有”的突变？

新发突变（de novo mutations，DNMs）指孩子基因组中新出现、而父母基因组中没有检测到的遗传变异。

在这项研究中，经过质量控制后，研究人员分析了8,328名活产婴儿，最终识别出390,924个新发单核苷酸变异（de novo single-nucleotide variants，dnSNVs）和28,032个新发插入缺失（de novo indels，dnIndels）。按可检测基因组区域校正后，每名孩子平均约有61.61个DNM。

其中，42,811个、约占全部DNM的10.22%，由于变异等位基因频率（variant allele frequency，VAF）低于0.35，被研究人员归类为早期合子后嵌合突变（early post-zygotic mosaic mutations，EPZMs）。这意味着它们可能不是来自父亲的精子或母亲的卵子，而是在受精后的早期胚胎细胞分裂过程中出现。

一个很容易被忽略的数据是：在390,924个dnSNV中，仅0.34%被归为功能性编码变异。绝大多数位于基因附近的非编码区域或基因间区。

因此，后文看到“多了1个”“多了2个”dnSNV时，不能直接理解成“多了1个或2个致病突变”。突变数量、突变位置和临床后果，是三个不同层次的问题。

父亲年龄像持续“加法”，母亲年龄的曲线却出现了弯折

父亲年龄与新发突变的关系过去已经得到大量研究支持，这项大规模中国人群数据再次观察到了非常清晰的效应。

在7,374名单胎儿中，在同时考虑父母双方年龄后，父亲受孕年龄每增加1岁，孩子平均增加1.06个父源DNM，其中约1.00个为dnSNV。

母亲的年龄效应较小：母亲年龄每增加1岁，平均增加0.38个母源DNM，其中0.36个为dnSNV。

但这里出现了一个值得关注的细节。

父亲年龄与dnSNV的关系基本可以用线性模型描述，而母亲的关系并不是一条完全笔直的线。研究人员发现，二次模型对母源dnSNV的拟合优于简单线性模型。在限制性立方样条（restricted cubic spline）分析中，31.81岁以上母亲的年龄效应斜率约为28.92岁以下母亲的1.52倍——0.42对0.28个dnSNV/年。

更有意思的是，这种非线性主要集中在C>G颠换（transversion）。年龄较大的母亲中，T(C>G)T和T(C>G)A等序列背景更加丰富，呈现出与APOBEC相关突变模式相似的特征。研究人员据此提出，卵母细胞老化可能不仅仅是“时间更长”这么简单，也可能伴随DNA损伤耐受或修复过程的改变，例如REV1依赖的跨损伤DNA合成（translesion synthesis，TLS）。

注意，这一部分属于机制推断：数据观察到了突变谱变化，但具体分子机制仍需实验验证。

ART不是一个暴露变量：多出的2.38个dnSNV，来源并不相同

如果只问“ART会不会增加新发突变”，很容易把问题问得过于粗糙。

在调整父母年龄和婴儿性别后，ART受孕的单胎婴儿平均比自然受孕（natural conception，NC）婴儿多2.48个DNM；如果只看dnSNV，则平均多2.38个。

1.02÷2.38约等于42.86%。换句话说，研究中与ART相关的额外dnSNV，接近一半表现为可能发生在受精以后早期胚胎发育阶段的嵌合突变。

这恰恰提示了这项研究最重要的思路转变：

一个ART周期可能包含促排卵、取卵、精子处理、IVF或卵胞浆内单精子注射、胚胎培养、冷冻复苏以及不同发育阶段的胚胎移植。不同操作发生在完全不同的生物学阶段，因此它们留下的突变信号也可能不同。

ICSI指向父源，GnRH-ant指向母源：同是ART，突变路径并不一样

当研究人员进一步拆分具体操作后，出现了两个很有辨识度的信号。

第一个是卵胞浆内单精子注射（intracytoplasmic sperm injection，ICSI）。

在同时考虑其他ART操作、父母年龄和婴儿性别后，与常规体外受精（in vitro fertilization，IVF）相比，ICSI与平均增加1.01个父源dnSNV相关。

如果直接与自然受孕相比，IVF组父源dnSNV增加0.41个，差异没有统计学显著性；ICSI组则增加1.58个，95%置信区间为0.74～2.43。

一个合理质疑是：接受ICSI的男性本来就更多存在少弱畸形精子症（oligoasthenoteratozoospermia，OAT），看到的会不会只是男性不育因素？

研究人员对此进行了分层分析。无论是否诊断OAT，ICSI与较多父源dnSNV之间的关系方向相似，异质性检验P=0.83。因此，男性因素不育本身似乎不足以完全解释这个信号。

但这仍然是观察性证据。是ICSI绕过了受精过程中的某些自然选择，还是显微操作产生的机械或生化应激参与其中，目前无法由这项研究区分。

第二个信号来自促排阶段使用的促性腺激素释放激素拮抗剂（gonadotropin-releasing hormone antagonist，GnRH-ant）。

GnRH-ant使用与平均增加1.69个母源dnSNV相关，而且研究观察到剂量—反应趋势。高剂量组与自然受孕相比平均增加1.65个母源dnSNV，而未使用GnRH-ant的ART组与自然受孕组之间没有显著差异。

这个换算非常直观，但也最容易被误解。它只是利用年龄—突变数量的回归斜率做出的量纲比较，不能理解为“使用GnRH-ant等于母亲老了4.69岁”，更不能直接转化为临床风险。

胚胎培养留下了另一种信号：不是父源，也不是母源

相比精子和卵子来源的突变，胚胎移植阶段出现的是另一幅图景——EPZM。

与自然受孕相比，冷冻囊胚移植（frozen blastocyst-stage embryo transfer）平均多1.23个EPZM，95%置信区间为1.02～1.44。

研究人员进一步查看突变类型，发现最突出的变化集中在C>A颠换。

囊胚移植后代的C>A EPZM比例为13.45%，高于卵裂期胚胎移植的12.28%，也高于自然受孕的11.95%。冷冻胚胎移植中的C>A比例为13.20%，而新鲜胚胎移植为12.04%，后者与自然受孕的11.95%几乎没有差异。

为什么偏偏是C>A？

C>A突变与活性氧（reactive oxygen species，ROS）导致的8-氧代鸟嘌呤（8-oxoguanine，8-oxoG）损伤具有已知的生物学联系。更进一步，囊胚移植与自然受孕之间的差异突变谱，与COSMIC中的SBS18突变特征最为相似，而SBS18同样与ROS介导的氧化DNA损伤有关。

这里仍然只能说“可能”。突变谱的相似性提供了机制线索，并不是对培养环境造成DNA氧化损伤的直接实验检测。

多一个突变，就等于多一分疾病风险吗？答案没有这么简单

研究的价值还在于，它没有停留在“测到了更多突变”，而是把基因组结果与出生和发育结局连接起来。

每增加1个父源dnSNV，与孕周缩短0.006周以及较低出生体重相关，同时早产和低出生体重的比值比（odds ratio，OR）均约为1.02。

单个突变对应的效应量其实很小。这一点很重要：群体水平上稳定的统计学关联，不意味着每增加一个突变就会产生可感知的个体健康变化。

研究还发现，整体父母来源dnSNV数量与1岁时的出生缺陷或神经发育指标没有显著关联。一个例外是位于已知先天性心脏病（congenital heart disease，CHD）基因中的功能性dnSNV，其CHD风险OR达到39，但95%置信区间为6～246，非常宽，反映出这类事件本身极为少见。

在中介分析（mediation analysis）中，父源dnSNV可以统计学解释父母年龄与孕周缩短关系的13.57%。对于ICSI，其增加的父源dnSNV可以解释ICSI与孕周缩短关联的3.99%，以及ICSI与早产风险关联的8.32%。

这同样不能被理解成已经证明了“突变导致早产”。中介分析可以提出可能的路径，但建立因果链仍要求更严格的设计和独立验证。

这意味着，如果这个结果能够被后续研究重复，未来值得研究的也许并不是笼统的“突变越多越差”，而是哪些突变类型、发生在胚胎发育的哪个时间窗口、落在哪些细胞谱系，才与后来的表型有关。

这项研究最容易被误读的，是把“相关”理解为“因果”

这项研究样本规模大，拥有亲子WGS、详细的ART操作记录以及出生和1岁随访数据，这是它的重要优势。但它依然是一项观察性研究。

研究人员明确指出三个主要限制。

首先，EPZM采用VAF<0.35进行定义，因此可能遗漏一部分合子后突变，也无法通过外周血检测捕捉组织特异性嵌合（tissue-specific mosaicism）。

其次，研究对象全部是活产儿。如果某些严重DNM导致胚胎停止发育、流产、死胎或终止妊娠，它们不会进入当前分析。因此，这项研究看到的本身就是经过妊娠和出生筛选后的群体。

第三，即使模型已经调整了大量父母、临床和操作因素，残余混杂（residual confounding）仍无法完全排除。

还有一个时间尺度问题：神经发育只随访到了1岁。对认知、语言、行为等高度动态的表型而言，一岁既重要，又远远不是终点。

比“ART安全吗”更值得关注的，是“哪一步、为什么、能不能优化”

这项研究带来的启发，并不是重新制造一个简单的“自然受孕 versus 辅助生殖”二元对立。

更重要的是，它把ART拆成了不同的生物学阶段：精子阶段出现父源信号，卵母细胞阶段出现母源信号，胚胎培养阶段则出现合子后嵌合信号。

这种分层视角为后续研究提供了更具体的问题。

ICSI相关信号究竟来自受精选择机制还是操作过程？GnRH-ant的剂量关系能否在独立人群中重复？培养时间、氧浓度、培养液成分以及冷冻复苏条件是否会改变C>A型EPZM？这些突变是否存在组织特异性，又是否能在儿童更长期的神经发育随访中重复观察？

回答这些问题，才有可能把“发现一个统计学关联”推进到“找到一个可以优化的生殖医学环节”。