标题：挑战性染色体进化经典模型 —— 多刺鱼古老同形 ZW 系统的多组学全景解析
期刊：Science Advances
论文：Evolution of a ZW sex determination system in sticklebacks
团队：中国科学院动物研究所郭宝成
标签：进化基因组、性染色体、同形性染色体、刺鱼、趋同演化、生物多样性、新质资源
一、思维导图（脑图）
脊椎动物性染色体演化长期遵循一套教科书式经典模型：常染色体获得性别决定位点→重组抑制逐步扩张→Y/W 积累转座子、大量基因丢失，最终形成形态差异巨大的异形性染色体（如人类 XY、鸟类 ZW）。但大量野外物种不断涌现反例：上亿年演化仍保持同形、无明显分化的性染色体，这类系统的演化逻辑始终悬而未决。
中国科学院动物研究所郭宝成团队以东北亚特有中华九刺鱼（Pungitius sinensis）及其近缘巴氏九刺鱼为核心研究对象，依托 PacBio、Hi-C、家系重测序多组学技术，首次完整解析了刺鱼属古老 ZW 性别决定系统：这套 ZW 系统已存在至少 500 万年，却始终维持 Z、W 高度同形，完全违背 “越古老分化越剧烈” 的经典预测。本研究打破传统演化框架，提出 “弱性拮抗选择 + 无大片段染色体重排” 共同约束性染色体分化的全新机制，同时产出两套高质量鱼类染色体基因组新质资源，为全球脊椎动物性染色体多样性研究提供标杆体系。
三、关键科学问题经典模型核心悖论：按照传统理论，500 万年演化足以让 Z/W 染色体发生大规模序列分化、W 基因丢失，但中华九刺鱼 Z/W 几乎无全局分化，为何古老性染色体能长期维持同形状态？阻碍重组抑制扩张的核心进化力量是什么？
系统演化溯源问题：多刺鱼属内部同时存在 XY（P. pungitius）和 ZW（中华九刺鱼支系）两套完全相反的性别决定系统，ZW 是由 XY 祖先转换而来吗？两套系统分化的时间节点、演化路径如何？
性别决定区（SDR）精细结构未知：此前仅粗略定位 12 号染色体为性染色体，缺乏单倍型分辨率的 SDR 边界、雌雄特异区段（FSR）、高分化区（HDR）精细划分，SDR 内哪些基因调控雌雄分化仍不明确。
重复序列驱动分化机制：转座子 TE 通常是重组抑制、W 退化的核心推手，但本研究中 TE 仅局部富集，全染色体无大范围扩张，TE 爆发的时序与 ZW 起源存在何种关联？
类群普适科学命题：同形性染色体稳定演化是否为脊椎动物普遍模式？刺鱼 ZW 系统能否为鸟类、两栖类古老同形性染色体提供统一解释框架？
高质量染色体基因组组装壁垒中华九刺鱼、巴氏九刺鱼杂合度高、重复序列占比 15% 左右，仅短读长无法搭建连续骨架；必须整合 PacBio HiFi 长读、Nanopore 超长读、Hi-C 染色质互作、二代短读多平台混合组装，流程复杂、测序成本高昂，野外鱼类样本 DNA 降解严重，高完整度核酸获取难度大。
Z/W 单倍型精准拆分技术难点雌雄个体混装基因组会混淆 Z、W 序列，需通过雌性个体长读长单倍型分型，才能分离两套性染色体独立序列，单倍型解析计算量巨大，常规组装流程无法实现。
高分辨率 SDR 定位门槛仅单一野生群体无法区分性别关联变异，必须构建 16 雄 16 雌人工家系，通过全基因组重测序逐窗口计算 FST、雌雄杂合度差异，海量 SNP 变异的统计校正、可视化分析耗时极长。
区分真实分化与测序噪声SDR 内局部覆盖度差异、杂合信号易被误判为性染色体分化，需同步校正测序深度、重复序列干扰、种群遗传背景三层干扰因素。
跨物种比较系统发育约束多刺鱼属存在种间基因渗入、染色体重排，构建可靠物种树、精准估算分化时间，需要数千个单拷贝直系同源基因联合建树，算力消耗巨大。
图1、刺鱼属性染色体演化
整体研究采用基因组组装→家系定位→单倍型精细解析→跨物种溯源→分子机制挖掘五层递进逻辑：
基因组基线搭建：从头组装中华九刺鱼、巴氏九刺鱼两套染色体水平基因组，完成基因、重复元件完整注释，建立多刺鱼属高质量比较基因组基准。
家系群体锁定 SDR：利用人工杂交家系雌雄重测序数据，全基因组滑动窗口计算雌雄 FST、杂合度、测序覆盖度差异，锁定 12 号染色体 2–7 Mb 为性别决定核心区 SDR，2.4–4.8 Mb 为功能核心区间。
Z/W 单倍型分区解析：对雌性个体进行单倍型拆分，对比 Z、W 完整序列，划分雌性特异区 FSR（W 独有基因、LTR 富集）、高分化区 HDR（序列差异显著但无独有基因），全局序列高度共线性。
系统发育溯源 ZW 起源：整合 6 种刺鱼全基因组单拷贝基因构建进化树，分子钟推算 ZW 系统至少 500 万年前在sinensis支系共同祖先中起源；与P. pungitius的 XY 系统独立演化，无同源 SDR 区段，证实刺鱼属发生性染色体转换。
分化停滞机制挖掘：对比 TE 转座子插入爆发时间，发现 TE 扩张与 ZW 起源同步，但未持续扩散；全基因组扫描未发现 SDR 区间大片段倒位；雌雄脑转录组性别偏倚基因数量极低，证明弱性拮抗选择、无结构重排、TE 扩张停滞三重约束抑制重组抑制扩张。
候选基因筛选：SDR 内注释 352 个蛋白编码基因，筛选 dmrt4、cyp26b1 等 14 个性腺分化关键基因，作为性别决定主控候选。
理论整合：基于全部组学证据，提出区别于经典退化模型的同形性染色体长期稳定演化新范式。
图2、家系重测序定位SDR
1. 多平台混合染色体基因组组装（PacBio+Nanopore+Hi-C+BGI）操作：长读搭建骨架，Hi-C 挂载染色体，多轮纠错提升 Contig 连续性，单倍型分型拆分雌性 Z/W 两套序列。优点：获得染色体级、单倍型分辨基因组，可直接对比 Z/W 结构差异，是性染色体精细解析唯一基础手段。缺点：测序、生信分析成本极高，杂合度高物种组装难度翻倍，野外低质量样本产出差。
2. 家系全基因组重测序雌雄分化扫描（FST、滑动窗口杂合度、覆盖度）操作：同一家系雌雄 32 个个体重测序，10kb 窗口统计雌雄遗传分化信号，定位 SDR 区间。优点：排除种群背景干扰，精准区分性别连锁变异，SDR 定位分辨率可达 Mb 级。缺点：需要人工繁育杂交家系，鱼类繁育周期长，样本构建门槛高。
3. 单倍型水平 Z/W 共线性、重复元件时序分析操作：拆分 Z、W 单倍型，注释 TE 类型，利用 Kimura 距离估算转座子插入年代。优点：直观区分 W 特有序列，量化分化时序，解释性染色体分化驱动力。缺点：单倍型分型算法复杂，长读数据依赖度极高。
4. 雌雄性腺 / 脑组织 RNA-seq 差异表达分析操作：雌雄多组织转录组，筛选性别偏倚表达基因，富集性腺发育通路。优点：快速锁定性别分化候选基因，区分性腺与全身性别调控差异。缺点：仅表达关联证据，无法证明基因是性别决定主控因子。
5. 全基因组系统发育 + 分子钟分化时间估算操作：8000 + 单拷贝基因串联建树，化石校正分子钟推算物种、性染色体分化时间。优点：精准溯源 ZW 系统演化节点，区分趋同与同源特征。缺点：依赖高质量直系同源基因集，算力消耗巨大。
七、核心数据集及其独特性1. 原创两套染色体级多刺鱼基因组（核心新质资源）物种：Pungitius sinensis、Pungitius bussei，Contig N50 分别达 10.47 Mb、5.67 Mb，BUSCO 完整度 > 96%，95% 以上序列锚定 21 条染色体；
独特性：全球首批中华九刺鱼支系染色体水平基因组，同时拆分 Z、W 独立单倍型，填补 ZW 刺鱼高质量基因组空白，是脊椎动物同形性染色体研究全新基准资源。
32 尾（16 雌 16 雄）同一家系全基因组 SNP 集，126 万高质量群体变异；独特性：完全消除种群遗传背景噪声，是鱼类 SDR 高分辨定位的稀缺家系资源。
3. Z/W 分拆单倍型基因组与重复元件注释库独立 Z、W 序列，FSR/HDR 分区注释、TE 插入年代定量数据；独特性：首次实现鱼类 ZW 性染色体单倍型分离，可直接对比雌雄染色体结构差异。
4. 多组织雌雄转录组数据库脑、性腺、肌肉、肝脏雌雄转录组，包含数千个性别偏倚基因表达矩阵；独特性：区分性腺特异性别调控与全身性别分化表达，为候选基因筛选提供转录证据。
5. 多刺鱼属跨物种比较基因组集6 种刺鱼全基因组共线性、同源基因集；独特性：覆盖 XY、ZW 两套性别系统近缘物种，直接支撑性染色体转换演化对比。
八、主要结果及其未尽之处（原文图 5）图3、ZW演化时序与分化对比图
（一）五大核心研究结果系统演化结论：中华九刺鱼、巴氏九刺鱼共享同源 ZW 性别决定系统，起源至少 500 万年前；同属近缘P. pungitius为独立 XY 系统，证明多刺鱼属发生性染色体转换，ZW 由 XY 祖先演化形成。
SDR 精细结构：性染色体为 12 号染色体，核心 SDR 区间 2–7 Mb，功能核心 2.4–4.8 Mb；Z/W 分为雌性特异区 FS（W 独有免疫相关基因）、高分化区 HDR，染色体全局共线性极强，整体高度同形。
分化停滞核心机制：500 万年未大规模退化源于三重约束：① 物种雌雄二型弱，性拮抗选择压力极低；② SDR 区间无大片段染色体倒位（倒位是重组抑制扩张核心推手）；③ W 上 LTR 转座子仅 500 万年前短期爆发，后续插入停滞，无法持续扩大非重组区。
基因筛选结果：SDR 注释 352 个蛋白基因，14 个已知性腺 / 性别分化基因（dmrt4、cyp26b1、tesk1），是潜在性别决定主控基因；W 特有 FSR 富集免疫 V 结构域蛋白基因。
演化范式革新：推翻 “性染色体必然随时间逐步退化异形” 经典模型，证明长期同形是稳定演化的正常结果，提出弱选择 + 无结构重排驱动同形染色体新演化框架。
未锁定单一主控性别决定基因：仅筛选候选基因，无基因敲除 / 敲入功能验证，无法确定哪个基因是 ZW 系统性别开关。
野外种群验证缺失：仅单一家系、少量样本，未沿东亚南北梯度采集野生种群，无法验证不同气候种群 SDR 分化程度差异。3 杂交转换机制空白：未开展 ZW×XY 刺鱼杂交实验，不清楚两套性别系统共存时的显隐性、基因渗入如何驱动性染色体转换。
剂量补偿机制未解析：未探究 Z 染色体基因剂量平衡调控，无法对比鸟类 ZW 剂量补偿通路差异。
跨门类普适性不足：仅局限刺鱼类群，未拓展两栖、鸟类古老同形染色体做统一对比。
候选基因体内功能验证建立九刺鱼 CRISPR 基因编辑体系，分别敲除dmrt4、cyp26b1等 SDR 关键基因，观察性腺发育表型，锁定主控性别决定基因。
大范围野外群体演化调查沿中国东北、朝鲜、日本梯度采集中华九刺鱼野生种群，检测 SDR 序列分化、W 特异基因频率，解析气候、种群瓶颈对性染色体分化的影响。
ZW/XY 杂交实验解析转换机制构建P. sinensis（ZW 雌）×P. pungitius（XY 雄）杂交群体，观测后代性别比例、SDR 基因渗入，揭示性染色体转换分子通路。
多门类同形性染色体统一理论构建整合鱼类、鸵鸟、蟾蜍古老同形 ZW/XY 数据，对比 TE 扩张、性拮抗强度、染色体重排，提炼跨脊椎动物通用稳定演化模型。
剂量补偿通路专项研究雌雄全基因组甲基化、转录组对比，挖掘 Z 染色体剂量平衡调控机制，与鸟类 ZW 系统进行演化对比。
刺鱼属完整性染色体演化谱系重建扩增东南亚、东北亚多刺鱼近缘物种基因组，还原 XY→ZW 完整转换时间线，厘清基因渗入、染色体重排在转换中的作用。
应用转化：性别控制育种基于 SDR 雌雄特异分子标记，开发九刺鱼性别快速鉴定试剂盒，为鱼类性别调控育种提供分子工具。
模式物种资源奠基多刺鱼属是脊椎动物性染色体多样性天然 “演化实验室”，兼具 XY、ZW 两套独立系统、频繁染色体转换、近缘可杂交等特征，生物多样性前期分类、野外种群调查为本研究提供核心材料，若无刺鱼多样性基础，无法开展跨系统对比。
野外生态背景支撑不同九刺鱼物种地理隔离、气候分化、繁殖行为（筑巢、雌雄二型）的多样性研究，直接解释 “性拮抗选择弱” 的生态根源，是解析 Z/W 不分化的前提。
比较类群标尺价值刺鱼科包含三刺鱼、四刺鱼、九刺鱼多个独立支系，多样性别系统构成天然对照，支撑 “性染色体转换” 演化推论，验证经典模型缺陷。
全球脊椎动物多样性拓展意义本研究填补硬骨鱼类古老 ZW 系统研究空白，完善脊椎动物性染色体多样性图谱，为理解全球动物性别决定机制演化多样性提供关键一环。
本研究产出多项进化基因组领域原创新质科研资源，具备长期复用价值：
两套高质量单倍型分辨染色体基因组中华九刺鱼、巴氏九刺鱼全新参考基因组，Z/W 染色体独立拆分序列，填补东北亚淡水鱼类基因组空白，可作为硬骨鱼性别演化通用参考数据库。
性别决定区 SDR 分子标记资源W 雌性特异 FSR 基因、特异 SNP、覆盖度差异标记，开发鱼类雌雄快速鉴定分子探针，服务水产育种、野外种群性别调查。
多刺鱼属系统发育分子标记集8395 个单拷贝直系同源基因，可用于多刺鱼属物种界定、种群分化、生物地理研究。
转座子时序注释数据库Z/W 染色体 TE 插入年代、类型、分布全套注释，为鱼类转座子介导演化提供标准化分析数据集。
雌雄多组织转录组表达矩阵性腺、脑组织性别分化基因表达资源，挖掘鱼类雌雄发育调控通路，支撑性别发育分子机制后续研究。
标准化家系重测序分析流程公开可复用的 SDR 定位生信代码，为其他鱼类、两栖类性染色体研究提供通用分析工具链。
长久以来，性染色体演化的教科书模型始终将 “分化、退化、异形化” 视为唯一终极方向，而本研究通过 500 万年古老同形 ZW 刺鱼系统，提供了强有力的反例：性染色体不一定走向衰败，弱性拮抗、缺乏染色体重排可以长期锁住重组抑制，维持 Z/W 高度同源。
依托多组学交叉手段，本研究产出两套鱼类染色体级新质基因组资源，厘清了多刺鱼属 ZW 性别系统的起源与演化轨迹，刷新了脊椎动物性染色体演化认知。但研究仍存在候选基因功能缺失、野外群体证据不足等短板，后续基因编辑、跨种群调查、杂交实验将进一步完善这一全新演化范式。
刺鱼作为全球经典进化模式生物，其多样性别系统为解码动物生殖、物种形成提供天然窗口。这项工作不仅修正演化理论，同时产出的基因组、分子标记资源，也将推动淡水鱼类生物多样性调查、水产性别控制育种等应用方向，让基础进化生物学研究真正兼顾理论突破与产业价值。
转载本文请联系原作者获取授权，同时请注明本文来自朱朝东科学网博客。
链接地址：https://wap.sciencenet.cn/blog-536560-1544254.html?mobile=1
收藏
