《Heredity》：Different sex determination systems in two closely related Eurasian minnow (Phoxinus) species
编辑推荐：
为解决欧亚鱥（Phoxinus）属鱼类性别决定系统不明的问题，研究人员对P. phoxinus和P. csikii进行了全基因组测序分析，结合覆盖度、SNP和k-mer等方法，鉴定出与性别相关的基因组区域。结果揭示这两个近缘物种拥有完全不同的性别决定系统：P. phoxinus为雄性异配（XX/XY），而P. csikii为雌性异配（ZZ/ZW）。这些不同的性别决定机制可能促进了物种间的生殖隔离，为理解鱼类性别决定系统的快速演化与物种形成提供了新见解。
在大自然的演化舞台上，性别决定机制展现出了令人惊叹的多样性，尤其是在硬骨鱼类中，这种多样性达到了极致。与哺乳动物相对固定的XY系统或鸟类相对固定的ZW系统不同，鱼类的性别决定方式可谓“花样百出”，既有遗传机制，也有环境因素（如温度）的影响，甚至在近缘物种间，也常演化出完全不同的性染色体系统。这种快速演化的背后机制是什么？它对物种的生殖隔离和物种形成又扮演着怎样的角色？这些都是演化生物学和遗传学领域令人着迷的课题。
欧亚鱥（Phoxinus），一种广泛分布于欧亚大陆淡水中的小型鱼类，为探究这些问题提供了一个理想的研究模型。这个类群已知存在杂交现象，但其性别决定的遗传基础却一直是个未解之谜。更引人入胜的是，这个类群的物种边界尚不完全清晰，且存在自然杂交区，这为我们研究性别决定系统的差异如何影响甚至塑造物种间的生殖屏障提供了绝佳机会。如果两个能杂交的物种拥有不同的性别决定系统，会产生怎样的遗传后果？这会是阻碍基因交流、促进物种分化的一道关键“壁垒”吗？为了解答这些疑问，研究人员将目光聚焦于两个在系统发育和地理分布上都十分接近的物种：P. phoxinus 和 P. csikii。
为了揭示这两种鱥鱼性别决定的“密码”，研究人员开展了一项整合多种基因组学方法的系统性研究。论文发表在《Heredity》杂志上。研究者们分别采集了26尾P. phoxinus（13雄13雌）和16尾P. csikii（8雄8雌）样本，这些样本分别来自莱茵河、默兹河、多瑙河等不同的流域。研究团队首先对样本进行了全基因组测序，并利用线粒体基因组构建了系统发育树，确认了物种的分类地位。随后，他们采用了一个名为SexFindR的生物信息学分析流程，该流程整合了三种互补的分析策略来鉴定性连锁区域：一是基于基因组覆盖度的分析，用于检测因性染色体分化（如Y或W染色体退化）导致的拷贝数差异；二是基于单核苷酸多态性（SNP）的分析，通过比较两性间的等位基因频率、群体分化指数（FST）、关联分析（GWAS）和核苷酸多样性（π）来寻找与性别高度关联的基因组区域；三是基于k-mer（k核苷酸短序列）的分析，这是一种不依赖于参考基因组的比对，能够独立地鉴定出性别特异的DNA序列。这种多管齐下的方法，大大提高了在非模式生物中鉴定尚未充分分化的同形性染色体的能力和可靠性。
结果
1. 线粒体基因组系统发育
基于12个线粒体蛋白编码基因构建的最大似然系统发育树显示，P. phoxinus 和 P. csikii 各自形成了具有高支持率的单系群，证实了研究中使用的样本分类准确，为后续物种特异性的性别决定分析奠定了基础。
2. P. phoxinus 的性别决定——XX/XY系统的证据
覆盖度分析未发现任何染色体在两性间存在显著的覆盖度差异，这与两者都具有同形性染色体的预期一致。然而，SNP和k-mer分析则揭示了清晰的性别关联信号。
- •
染色体3和12上的性别关联区域：SNP分析鉴定出两个候选的性别关联区域，一个位于3号染色体（47.14–47.16 Mb），另一个位于12号染色体（21.27–21.28 Mb）。k-mer分析进一步支持了3号染色体是潜在的性染色体，所组装的性别特异k-mer序列富集于该染色体上的gpsm1基因内部。
- •
雄性异配（XY）的遗传模式：在这两个区域，基因型分析显示雄性个体普遍呈现杂合状态，而雌性个体则主要为纯合状态。这一模式是典型的XX/XY（雄性异配）系统的特征，即雄性携带不同的X和Y染色体，而雌性携带两条相同的X染色体。
- •
流域特异性模式：有趣的是，这种性别关联信号在不同流域的种群中表现出差异。3号染色体上的信号在莱茵河种群中更强，而12号染色体上的信号则在默兹河种群中更显著。这暗示P. phoxinus 中可能存在种群特异的性别决定位点，或正处于性染色体更替的早期阶段。
- •
功能基因注释：3号染色体关联区域内的连锁区块包含rab14、abl1、exosc2和spr等基因，其中一些与精子发生、细胞信号传导和色素沉着有关。12号染色体区域则包含gabrr2、rragd等基因，与TOR信号通路等相关。
3. P. csikii 的性别决定——ZZ/ZW系统的证据
与P. phoxinus 类似，覆盖度分析也未发现异形性染色体。SNP分析未能鉴定出跨所有指标一致的性别关联区域。然而，不依赖于比对的k-mer分析却取得了关键发现。
- •
k-mer分析指向3号染色体：k-mer分析鉴定出大量性别特异的k-mer序列，其中超过90%为雌性特异。这些序列被组装成133个重叠群，其中许多富集于3号染色体的起始区域（约1.49–1.50 Mb 和 1.95–1.96 Mb），这与P. phoxinus 的性别关联区域位置不同，并非同源。
- •
雌性异配（ZW）的遗传模式：在1.49 Mb附近的区域，研究人员找到了76个SNP。基因型分析显示，雌性个体在这些位点普遍为杂合子，而雄性个体则为纯合子。这正是ZZ/ZW（雌性异配）系统的典型特征，即雌性携带不同的Z和W染色体，雄性携带两条相同的Z染色体。
- •
功能基因注释：该性别关联区域内的基因富集于“卵巢类固醇生成”和“GnRH信号通路”等与雌性生殖相关的通路，为ZZ/ZW系统提供了功能上的佐证。
结论与讨论
本研究通过综合运用全基因组测序和多种生物信息学方法，首次揭示了两种近缘欧亚鱥——P. phoxinus 和 P. csikii——拥有截然不同的性别决定系统。P. phoxinus 表现为雄性异配的XX/XY系统，而P. csikii 则表现为雌性异配的ZZ/ZW系统。这一发现具有多重重要意义：
首先，它为硬骨鱼类性别决定系统的极端可塑性和快速演化提供了又一个鲜活例证。即使在能够杂交的近缘物种间，控制性别这一最基本生物学特征的遗传机制也可能完全不同。这种差异可能是由性染色体周转（sex chromosome turnover）驱动的，即新的性别决定基因在不同染色体上起源并固定，而旧的退化性染色体丢失。
其次，不同的性别决定系统可能构成重要的生殖隔离屏障。当拥有XY系统的雄性P. phoxinus 与拥有ZW系统的雌性P. csikii 杂交时，其后代的性别染色体组成（如XW或YZ）可能导致性别发育异常、性比失衡或杂交后代不育，从而有效阻止基因流动，促进物种维持其独立性。因此，性别决定系统的分化不仅是物种分化的结果，也可能是推动物种形成的驱动力之一。
再者，研究揭示了P. phoxinus 内部可能存在种群特异性的性别决定区域（3号与12号染色体），这暗示其性别决定机制可能比目前认识的更为复杂，或许正处于动态演化过程中，为研究性染色体系统的早期起源和建立提供了宝贵材料。
最后，在方法论上，本研究成功展示了整合覆盖度、SNP和k-mer分析在鉴定同形性染色体方面的强大效力，特别是k-mer分析在参考基因组可能存在偏差或性别分化信号微弱时展现出的独特优势。
总之，这项研究不仅解开了两种欧亚鱥的性别决定之谜，更将其置于演化生物学的大图景中，凸显了性别决定机制在物种分化与生物多样性形成中的关键作用。未来，对这两个物种及其杂交带的进一步研究，将有助于我们更深入地理解性染色体演化、生殖隔离形成以及物种边界维持的遗传与演化机制。
