2020年11月2日，国际学术权威刊物《美国科学院院报（PNAS）》在线发表福建农林大学基因组中心张积森课题组完成的“金鱼的进化起源和驯养历史（The evolutionary origin and domestication history of goldfish (Carassius auratus)）” 研究成果，深入探讨了金鱼的起源、演化、驯化和人工选择。

金鱼属于鲤形目、鲤亚科、鲫属，形态优美，身姿奇异，色彩绚丽，有草金、文、蛋、龙金、龙背等五大类，250个品种之多，是世界著名三大观赏鱼类之一，达尔文形容它有着无限多样的颜色和非凡的表型结构变化。金鱼起源于中国，有“国粹”之称，在中国传统文化中有着吉庆、富有的象征，深受老百姓的喜爱和推崇。福州是传统三大金鱼产区之一，有三百多年的发展历史，金鱼养殖产业发展迅速，有“中国金鱼之都”之称。金鱼的研究对进一步促进福州金鱼产业的发展有着重要的意义。

图1 福寿（福州）供图

基于金鱼的重要文化、观赏及科学研究价值，特别是其丰富的变异性在动植物中相当罕见，课题组从以下几个方面做了深入研究。首先，完善并提升高质量基因组。研究人员使用PacBio单分子测序、染色质三维构象捕获 (Hi-C)，结合光学图谱BioNano基因组研究技术，构建了高质量基因图谱，将基因组95.75% 重叠群（contigs）锚定到50条假染色体上，获得单倍型参考基因组，其大小为1.73Gbp，scaffold N50 为 606kb，基因组完整性达94.7%，同时预测了50条染色体中38条的潜在着丝粒区域，更加突出了组装的基因组的完整性。

其次，解析金鱼亚基因组及其在鲤亚科中的演化。多倍化现象在脊椎动物中极为罕见，以金鱼祖先鲫为代表的鲤亚科鱼类在演化历史上可能经历了四轮全基因组加倍事件。针对金鱼的基因组是来自于两个祖先物种的杂交事件(也称异源多倍体)还是来自于同源多倍化的问题，科研人员在前人研究基础上，通过对6个鲤科代表性二倍体（2n=50）物种进行重测序，观察比对率发现，䰾鱼平均80%的读段比对到金鱼的一个染色体簇群上，将金鱼50条染色体明显分为两个亚基因组A（ChrA01-A25）、B（ChrB01-B25），且两个基因组大小相近。研究人员通过结缔组织生长因子基因分析验证了之前关于金鱼为四倍体，是异源杂交染色体加倍的假说，而䰾亚科（Barbinae）为最有可能的二倍体祖先种之一，推测金鱼祖先在1000至1200万年前发生过古老的杂交事件

基于同义替换率Ks计算， A、B亚基因组大约分化于14.84-16.67 百万年前，先于金鱼和鲤鱼的分化（11.72-13.16 百万年前），暗示产生异源多倍体金鱼的祖先种早于金鱼与鲤鱼的分化，同时表明全基因组复制(WGD)事件发生在金鱼和普通鲤鱼的共同血统中，这也支持了他们拥有相同染色体数目的事实。金鱼A、B亚基因组间整体共线性较强，但在6个同源染色体对中存在较大的片段反转，进一步支持了金鱼是异源四倍体的假设。有意思的是金鱼在Chr4/Chr29上存在高度不共线性的染色体区域，功能富集结果显示这些区域在免疫应答、疾病抗性上起着重要的作用，推测对鲤科物种为适应多种病原体而在不同自然条件下栖息和生活具有促进作用。通过结合RNA-seq数据，观察金鱼AB亚基因组在10个不同组织中的差异表达情况，发现同源基因在两个亚基因组之间表现出不对称的表达模式，而全基因组范围内的表达水平总体上更偏好A亚基因组。

图2 金鱼演化示意图

第三，通过群体研究金鱼家驯化进程。研究人员对185个具有代表性的金鱼品种和16条研究人员对185个具有代表性的金鱼品种和16条野生鲫鱼个体进行重测序，结合33个前人鲤鱼测序数据（三者无生殖隔离），通过SNP位点构建系统发育树显示，金鱼与鲫鱼的亲缘关系高于鲤鱼。主成分分析（PCA）将鲫鱼与金鱼明显分开，进一步佐证了金鱼起源于鲫鱼的假设。主成分分析、系统发育树、群体结构分析、基因流分析等多重证据显示，所有的金鱼品种共享一定数量的血缘，金鱼经历了野生鲫鱼到草金，接着发展到文种，再产生蛋种的家驯化历史，而三分法（草系、蛋系和文系）更符合金鱼的进化历史。

图3 金鱼的种群结构

注：a ,b 不同金鱼种类群体 SNP变异的主成分析；c, 金鱼群体系统发育树；d, Treemix 历史混群检测；e, 群体结构分析

此外，研究人员将遗传变异与许多典型特征（眼睛、背鳍、鳞片、颜色）联系起来，通过遗传分化指数（FST）分析，发现50个受选择区域，选择性清除分析发现393个受选择的基因中，有13个基因与斑马鱼中存在敲除系基因同源，这些敲除系基因与行为(pcdh15a和agrp)、眼部萎缩(ndufs5)、细胞迁移(nav3和zap70)或脑部萎缩（aurka）等相关。推测在整个金鱼群体中这些基因有可能是目标性状在选择过程中的直接靶标，也有可能参与与金鱼主要原始驯化性状有关的表型分化。

第四，金鱼重要典型表型性状的组学研究。利用全基因组关联分析（GWAS）定位到一个与透明鳞相关的基因，为酪氨酸蛋白激酶受体基因，这是一个在中国遗传学研究史上著名的孟德尔性状，为不完全显性(T/ t基因座)，该基因早在1928年被陈桢教授描述过，该基因在斑马鱼中与色素细胞的发育有关，受体酪氨酸激酶家族的成员在人类黑色素瘤的发展中发挥作用。背鳍是区分文种和蛋种分类的主要标志，研究人员通过对文种和蛋种杂交后代F1的背鳍性状分离比进行统计，推断背鳍特征受多个基因位点控制，并具有母体遗传效应。关联分析在13条染色体的8.96 Mb区域上鉴定出了378个与背鳍表型具有相关性的基因，包括背侧纵向血管（itgb1a，rab13和e2f8），背侧异常（kif5ba），背主动脉异常（fev，cxcr4b，gpr183a，tal1，grlx2和uchl5），腹侧（zic2a）和弯曲背鳍（atp1a2a和dhfr），这些基因可能在文种金鱼通过人工选择驯化向蛋种金鱼进化的过程以及背鳍发育过程中发挥了关键作用。本研究中发布的高质量金鱼基因组和多样性数据，为金鱼作为大量候选功能基因的模型以及广泛的家养物种进化案例研究提供了新的知识和资源。

图4 金鱼透明鳞相关性状的全基因组关联分析

注：a，文种和蛋种；b，文种和蛋种的群体主成分分析；c， 正对透明鳞性状的曼哈顿图；d，一般线性模型下GWAS的分位数图；e，位于ChrB15染色体上候选基因位置；f，候选区域的基因型分析；g，候选基因及-UTR调控元件

总的来说，本研究提供了高质量基因组，为金鱼分子育种提供基础，有助于推动品种多样性选择和产业的发展。通过亚基因组的研究，为研究异源多倍体脊椎动物的基因组演化和功能提供了新的思路；通过对金鱼系列表型性状的挖掘与解析，为解析脊椎鱼类多样性的演化机制提供多的研究线索，有助于将金鱼成为脊椎动物生理生化和进化优良遗传模型，来探讨脊椎动物突变的遗传机制。

该研究由福建师范大学陈由强教授、福建农林大学张积森教授联合发起，福建师范大学生命科学学院、福建农林大学基因组与生物技术研究中心、闽江学院、比利时根特大学等单位协作完成。福建师范大学博士生陈多为第一作者，福建农林大学博士生张清、闽江学院唐唯其博士、福建师范大学黄镇讲师、福建农林大学王刚博士为共同第一作者。福建农林大学的张积森教授为本研究通讯作者，福建师范大学陈由强教授、福建农林大学唐海宝教授，比利时根特大学Van de Peer教授、福建师范大学黄镇讲师为共同通讯作者。