PNAS|金鱼的进化起源和驯化史 |

发布于 2020年11月3日

2020年11月03日，北京百迈客生物科技有限公司合作项目金鱼基因组及群体驯化研究成果以“The evolutionary origin and domestication history of goldfish(Carassius auratus)”为题在国际知名期刊PNAS上在线发表。文中组装获得了高质量的金鱼基因组 (2n=100)的同时，通过比较基因组学研究解析了古代杂交事件产生了两个亚基因组；进一步通过185条代表性金鱼变种和16条野生鲫鱼的重测序分析，揭示了金鱼的起源，并鉴定了与驯化相关的受选择基因区域。研究中全面收集了金鱼品种，获得了与众多解剖特征相关遗传变异信息，其中包括区分传统金鱼进化的特征。除此之外，鉴定了酪氨酸蛋白激酶受体作为透明突变体的候选基因。为金鱼的基因组和多样性数据提供了丰富的资源，使得金鱼成为功能基因组学和驯化研究具前景的模型。

金鱼（Carassius auratus）在古代是由鲫鱼驯化而来的（二者仍然被视为同一物种），是重要的养殖鱼之一，由于金鱼可以生产成千上万的卵，并且在同一个池塘中可以养成数十条这样的小鱼，因此通过强烈的人工选择已经完成了千年的驯化，在描述金鱼时，查尔斯·达尔文（Charles Darwin）曾写道：“经过几乎无限的颜色变化，经历了结构上的修饰。强调了金鱼为调查各种形态特征的遗传学研究提供了丰富的资源。经过1000多年的驯化和育种，金鱼的体型，鳍结构，眼型和颜色都有数百种变异，这使金鱼成为鱼类生理和进化的遗传模型。以此为背景，本研究取得了如下的研究成果：

研究结果1：金鱼基因组组装

研究中主要利用三代测序PacBio RS II，ALLHi-C和光学图谱BioNano基因组研究技术，组装获得一雌性金鱼个体基因组，组装获得金鱼基因组大小约1.73 Gb（流式预估基因组大小1.8 Gb），Scaffold N50= 606 Kb，利用ALLHi-C将95.75%的Contigs锚定到50条假染色体上，进一步利用8,487个markers和50个连锁群的鲤鱼图谱数据，证明scaffolds排序与图谱彼此高度一致。组装获得的高质量金鱼，比近期发表的“Wakin”金鱼基因组草图更完整，基因组完整性达94.7%。共鉴定了56,251个编码基因和10,098个长链非编码transcripts，其中转座元件占基因组33.04 %，预测了50条染色体中38条的潜在着丝粒区域，再次突出了组装的金鱼基因组的接近完整性（图1 A和B）。

图1 金鱼A亚基因组和B亚基因组概况（A）；金鱼Hi-C互作热图（B）

研究结果2：金鱼亚基因组揭秘

在过去的50年里，一直推测金鱼基因组是四倍体。近期研究表明，鲃亚科Barbinae (2n = 50) 中的鲤鱼科可能是金鱼最接近的二倍体祖先，进而形成这样一个假说：鲃亚科（Barbinae）可能是二倍体金鱼谱系的祖先种，研究中进一步进行了6个代表性的二倍体鲃亚科物种（包括条纹小鲃，大鳞高须鱼，长臀鲃，黑鳍袋唇鱼，施氏魮和裂峡鲃）的全基因组鸟枪测序并进行基因组组装。金鱼染色体可清晰分为两个亚基因组，通过reads多序列比对，发现金鱼来自两个祖先世系，其中一个是鲃亚科共有的，进一步将与鲃亚科比例高的reads所在染色体定为A亚基因组（ChrA01~A25)并将其余染色体定为B亚基因组 (ChrB01~B25)（图2）。此外，基于1个核基因的系统进化分析，即结缔组织生长因子like基因在二倍体金鱼中仅有1个拷贝，但在四倍体金鱼中含有2个拷贝，结果显示鲃亚科物种的同源基因与金鱼A亚基因组同源基因聚类到一起（斑马鱼和鲤鱼作外群）。研究结果表明，A亚基因组可能起源于鲃亚科内的祖先物种，而B亚基因组的二倍体祖先可能已经灭绝或起源于未知的鲤科鱼类谱系，推测金鱼的起源可能是由于异源多倍化事件。金鱼的A和B亚基因组分别包含28,133和26,141个基因。同时通过金鱼10个不同组织的RNA-seq测序，探索了A和B两个亚基因组之间的差异表达，发现同源基因在两个亚基因组之间表现出不对称的表达模式，而全基因组范围内的表达水平总体上更偏好A亚基因组。

图2 鲃亚科条纹小鲃Puntius semifasciolatus Reads与金鱼的比对效率

研究结果3：金鱼基因组进化研究

在同源染色体对ChrA01/ChrB01, ChrA04/ChrB04, ChrA05/ChrB05, ChrA09/ChrB09, ChrA11/ChrB11和ChrA15/ChrB15上出现了大片段倒位，尽管整体共线性很强。这些大的片段倒位可能是异源四倍体形成的原因之一。进一步发现，金鱼两个亚基因组（Ks= ~0.17）与鲤鱼（Ks= ~0.21）之间的分歧时间（13.28-16.67 MYA），早于金鱼和鲤鱼之间的分歧时间（7.81〜8.77 MYA）（Ks =〜0.10）（图3 A）。此外，文中假设全基因组复制（WGD）事件发生在金鱼和鲤鱼的共有谱系中，因此它们可能具有相同数量的染色体（2n = 100），染色体基数（2n=50）是二倍体鲤亚科成员的两倍，暗示这两个种都是四倍体。同时，发现斑马鱼于~ 32.8-36.8 MYA，与鲫鱼、鲤鱼和草鱼由共同的祖先（LCA）分化而来（图3 B）。

图3 金鱼、鲤鱼、草鱼、条纹小鲃和斑马鱼的全基因组复制事件（A）和系统进化分析（B）

研究结果4：金鱼群体进化研究

通过201个个体（16只野生鲫鱼和185只代表性金鱼变种，平均测序深度12.5 x）的系统发育树构建发现普通金鱼与鲫鱼亲缘关系最近（图4 A-C），进一步结合系统发育树和PCA分析确定各金鱼变种的系统分化关系，区分普通金鱼和卵金鱼（无背鳍）和温金鱼（有背鳍），无背鳍金鱼和有背鳍金鱼进一步分为不同的亚群。

图4 金鱼群体结构分析

LD及群体选择分析支持了前期的假设，即金鱼驯化过程中种群遗传瓶颈效应的存在和强烈的人工选择（图5 A）。进一步的遗传多样性分析发现，由于人工选择的结果，从鲫鱼和普通金鱼的驯化过程中遗传多样性的增加，表明了从鲫鱼驯化后，金鱼重要遗传变异的积累（图5 B）；33个金鱼个体和16个野生鲫鱼个体的群体选择分析发现了50个受选择区域（25.2 Mb区间内，共计946个基因）（图5 C），940个候选基因中大部分与形态发生，色素沉着，行为，免疫反应/传染病，能量代谢和对激素的反应有关；进一步通过以上201个个体的选择性清除分析发现393个基因受选择，其中21个代表性候选基因中，有13个与在斑马鱼中存在敲除系的基因同源，这些敲除系显示与行为相关的突变表型（pcdh15a和agrp），眼睛大小减小（ndufs5），细胞迁移（nav3和zap70）或脑大小减小（aurka）。

图5 金鱼群体遗传多样性和群体选择分析

研究结果5：金鱼表型性状GWAS关联分析

背鳍作为金鱼一重要的表型特征，在不同变种中形态不一，背鳍是区分Wen金鱼和Egg金鱼的显着特征，背鳍的缩小是Egg金鱼进化的关键特征。通过Wen金鱼和Egg金鱼杂交后代F1背鳍性状分离比，暗示了了背鳍特征受多个基因位点控制，可能具有母体遗传效应。对96只有背鳍和87只无背鳍金鱼的全基因组关联分析GWAS研究分析发现跨越13条染色体上的8.96 Mb区域上的378个基因与背鳍性状相关（图6），85.2% (322)个基因位于A亚基因组，大部分位于五条染色体上，包括ChrA09，ChrA07，ChrA22，ChrA16和ChrA05，比位于B亚基因组上的背鳍相关基因的数量高得多，这表明背鳍相关基因不对称的亚基因组分布，GO和KEGG分析表明，该性状的候选基因可能参与了“细胞表面受体信号传导”途径，“信号转导”，“跨膜转运”，“骨骼系统发育”和“初级代谢过程和器官氮化合物代谢过程”。进一步通过与斑马鱼直系同源基因基因敲除品系表现出的突变表型分析发现，二氢叶酸还原酶（dhfr）可能在某种程度上与鱼鳍是否存在相关，该基因可能在金鱼人工选择过程中，在有无背鳍这一表型中发挥重要作用。

图6 金鱼背鳍相关性状GWAS关联分析

透明的金鱼突变体首次记录于1579年，是孟德尔遗传学在鱼类甚至脊椎动物中的典型案例，在T / t基因座上的优势表达不完全。本研究中利用48个实验个体和127个对照个体进行全基因组关联分析（图7 A-C），发现在B15染色体臂上的5个SNPs与T / t基因座上的不完全优势遗传模式相吻合。与同源染色体15相比，在该区域发生了倒位。在候选区域ChrB15：28,558−153,105中，鉴定了基因Cau.B15G0000010(ChrB15:48,888−76,549)，其编码酪氨酸蛋白激酶受体UFO-like蛋白，在5′-UTR区的预测转录因子结合位点包含四个SNPs（图7 G）。白细胞酪氨酸激酶与斑马鱼的色素细胞发育有关，并且受体酪氨酸激酶家族的成员在人类黑色素瘤的发育中起着作用，尽管目前该基因家族的作用在金鱼中尚不清楚。在金鱼中，酪氨酸蛋白激酶受体可能调节酪氨酸激酶，从而导致色素沉着变化，在该区域内，编码具有非同义单核苷酸突变的β-分泌酶1-like蛋白（bace1）的Cau.B15G0000020基因与斑马鱼和人的黑素细胞迁移有关，这可能解释在透明鳞片的金鱼上出现黑色污迹或斑点。尽管功能验证仍在进行中，但研究中认为酪氨酸蛋白激酶受体是金鱼T / t基因座的强力候选基因。

图7 金鱼透明鳞片性状GWAS关联分析

研究团队介绍

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