引言
随着生育年龄推迟和辅助生殖技术（assisted reproductive technology，ART）的广泛应用，一个越来越难回避的问题是：父母年龄和不同的辅助生殖操作，会不会影响孩子出生时携带的新发突变？这些变化又是否与孩子的健康有关？
近日，《Nature Medicine》的研究报道“Large-scale whole-genome sequencing reveals the landscape and health implications of de novo mutations”，尝试回答这个问题。研究人员对江苏出生队列中的24,030人、7,851个亲子家庭进行了全基因组测序（whole-genome sequencing，WGS），并将基因组数据与出生结局及1岁时的发育随访连接起来。结果比简单的“有影响”或“没影响”复杂得多。
先把尺度看清：一个孩子平均有多少个“父母都没有”的突变？
新发突变（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/年。
需要特别避免的误读：31.81岁并不是一个新的“高龄风险临界值”。它来自这批研究对象年龄分布的模型节点，数据提示的是母源突变累积在大约29～32岁附近开始出现加速趋势，而不是“过了32岁突然发生变化”。
更有意思的是，这种非线性主要集中在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个。
父源dnSNV增加 0.74 个
母源dnSNV增加 0.62 个
EPZM增加 1.02 个
1.02÷2.38约等于42.86%。换句话说，研究中与ART相关的额外dnSNV，接近一半表现为可能发生在受精以后早期胚胎发育阶段的嵌合突变。
这恰恰提示了这项研究最重要的思路转变：
“ART”不是一次单一的暴露。
一个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组与自然受孕组之间没有显著差异。
该研究甚至将这些效应换算成父母自然衰老的突变量级：ICSI的父源效应约相当于1.01年的父亲年龄增长，GnRH-ant的母源效应约相当于4.69年的母亲年龄增长。
这个换算非常直观，但也最容易被误解。它只是利用年龄—突变数量的回归斜率做出的量纲比较，不能理解为“使用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%。
这同样不能被理解成已经证明了“突变导致早产”。中介分析可以提出可能的路径，但建立因果链仍要求更严格的设计和独立验证。
另一个需要关注的结果来自C>A型EPZM：每增加1个C>A EPZM，1岁时认知发育迟缓的风险增加14%，OR为1.14，95%置信区间1.01～1.28。与此同时，EPZM总数量本身与神经发育迟缓并没有显著关系。
这意味着，如果这个结果能够被后续研究重复，未来值得研究的也许并不是笼统的“突变越多越差”，而是哪些突变类型、发生在胚胎发育的哪个时间窗口、落在哪些细胞谱系，才与后来的表型有关。
这项研究最容易被误读的，是把“相关”理解为“因果”
这项研究样本规模大，拥有亲子WGS、详细的ART操作记录以及出生和1岁随访数据，这是它的重要优势。但它依然是一项观察性研究。
研究人员明确指出三个主要限制。
首先，EPZM采用VAF<0.35进行定义，因此可能遗漏一部分合子后突变，也无法通过外周血检测捕捉组织特异性嵌合（tissue-specific mosaicism）。
其次，研究对象全部是活产儿。如果某些严重DNM导致胚胎停止发育、流产、死胎或终止妊娠，它们不会进入当前分析。因此，这项研究看到的本身就是经过妊娠和出生筛选后的群体。
第三，即使模型已经调整了大量父母、临床和操作因素，残余混杂（residual confounding）仍无法完全排除。
还有一个时间尺度问题：神经发育只随访到了1岁。对认知、语言、行为等高度动态的表型而言，一岁既重要，又远远不是终点。
所以，这项研究不能支持“某一种ART操作会导致孩子疾病”这样的结论，更不能仅凭这些数据改变个体患者的治疗方案。
比“ART安全吗”更值得关注的，是“哪一步、为什么、能不能优化”
这项研究带来的启发，并不是重新制造一个简单的“自然受孕 versus 辅助生殖”二元对立。
更重要的是，它把ART拆成了不同的生物学阶段：精子阶段出现父源信号，卵母细胞阶段出现母源信号，胚胎培养阶段则出现合子后嵌合信号。
这种分层视角为后续研究提供了更具体的问题。
ICSI相关信号究竟来自受精选择机制还是操作过程？GnRH-ant的剂量关系能否在独立人群中重复？培养时间、氧浓度、培养液成分以及冷冻复苏条件是否会改变C>A型EPZM？这些突变是否存在组织特异性，又是否能在儿童更长期的神经发育随访中重复观察？
回答这些问题，才有可能把“发现一个统计学关联”推进到“找到一个可以优化的生殖医学环节”。
当辅助生殖进入精细化时代，评价一种技术的目标不应只是“能否获得妊娠和活产”，还应进一步关注，我们能否在维持治疗成功率的同时，把配子和早期胚胎所经历的额外生物学扰动降到更低。
参考文献
Qin N, Dai J, Jiang Y, Xie Y, Lv H, Zhou J, Wang X, Jiang Y, Du J, Chen Y, Xu X, Hu L, Chen C, Xu Z, Zhu J, Liu Y, Wang Y, Song C, Gu Y, Jiang T, Xu B, Han X, Lin Y, Meng Q, Ling X, Ma H, Jin G, Wang C, Shen H, Hu Z. Large-scale whole-genome sequencing reveals the landscape and health implications of de novo mutations. Nat Med. 2026 Aug 7. doi: 10.1038/s41591-026-04585-2. Epub ahead of print. PMID: 42567929.
热门跟贴