1、基于基因组微卫星的团头鲂养殖群体遗传多样性分析和指纹图谱构建张 芹，王延晖，王冰柯，屈长义，刘英杰（河南省水产科学研究院，河南郑州；商丘市水产技术推广站，河南商丘）摘要 利用基因组微卫星检测 个团头鲂（）养殖群体的遗传多样性，个位点共检测到 个等位基因，大小在 。各位点观测等位基因数（）在 个，平均为 个；有效等位基因数（）在，平均为。各位点检测到的香农信息指数（）在，平均为。观测杂合度（）在，平均值为；期望杂合度（）在，平均值为。多态信息含量指数（）在，平均值为。群体间的遗传分化指数（）为，群体间的基因交流（）为，说明 和 个群体之间有一定程度基因流。个位点检测到的近交系数（），其他 个位点
2、近交系数均。研究表明，和 两个人工养殖群体中均有较高的遗传多样性，群体间有一定程度的基因交流，指纹图谱的构建可以用来进行团头鲂种质资源品种间的鉴定以及品种内的亲缘关系研究，尤其是可以区分没有亲缘关系数据背景的群体，可以指导没有家系构建能力的苗种繁育场进行繁殖策略的升级。关键词 团头鲂；基因组；多态性中图分类号 文献标识码 文章编号（）：开放科学（资源服务）标识码（）：，（，），；基金项目 河南省基本科研业务费项目（）；国家淡水水产种质资源库（：）；河南省农业产业技术体系（）；现代农业产业技术体系建设专项（）。作者简介 张芹（），女，河南新乡人，高级水产师，从事鱼类遗传育种研究。通信作者，高级水
3、产师，从事生态渔业研究。鸣 谢 感谢华中农业大学王卫民教授在团头鲂样品采集过程中给予的帮助！收稿日期 团头鲂（）为我国特有的淡水经济鱼类，因其生长快、成活率高、易捕捞、肉质好等优点，在全国各地均有养殖。团头鲂“华海 号”水产新品种经第五届全国水产原种和良种审定委员会审定通过，具有遗传性状稳定、生长快、成活率高等优良特征。“伊河鲂”在分类上属于鲂属团头鲂，是近年来开发黄河流域土著鱼类品种时，选育出来的一个优良新品系。微卫星标记（）由于重复性好，多态性高，广泛应用于水生动物的种质资源鉴定以及遗传多样性研究，尤其是利用 的多态性对水生生物进行亲子鉴定的较多。李宝玉等利用微卫星标记对团头鲂二倍体和三倍
4、体群体进行遗传特性分析，等利用微卫星技术对团头鲂群体不同家系的肌间刺数量进行研究，而有关团头鲂基因组微卫星研究鲜见报道。笔者对团头鲂基因组 序列进行分析，筛选得到若干多态性微卫星 标记，笔者利用这些标记，对 个团头鲂养殖群体进行遗传多样性研究和指纹图谱的构建，旨在为团头鲂种质资源品种间的鉴定以及品种内的亲缘关系研究提供科学依据。材料与方法 材料 在河南省嵩县陆浑水库伊河鲂鱼种厂采集团头鲂伊河群体 尾，以下简称 群体；从湖北省鄂州国家级团头鲂原种场采集团头鲂“华海 号”群体 尾，以下简称 群体；在湖北省鄂州市梁子湖采集野生团头鲂个体 尾，以下简称 个体。所有个体测量体长、体重，剪尾鳍用无水乙醇保
5、存备用。方法 团头鲂 的提取。团头鲂 用磁珠法基因组 提取试剂盒（天根生化）提取，经 琼脂糖凝胶电泳检测，用超微量分光光度计（）检测 浓度和纯度，置于 低温保存。安徽农业科学，（）：扩增及 检测。挑选出 对微卫星引物，使用 通用接头序列（）加到每对引物的 引物的 方向，并合成带不同荧光基团的 接头序列。反应体系（）：，基因组 模板 （），。扩增条件：预变性 （变性 ，退火 ，延伸 ）共 个循环，终延伸 ，保存。荧光毛细管电泳检测。取 荧光 产物进行琼脂糖凝胶电泳鉴定，检测 条带是否单一、片段大小是否与预期一致。条带单一且大小相符的，对照 的浓度进行定量，将所有产物稀释至相同的浓度范围，取 稀释
6、产物加 甲酰胺（含 内标）变性后上 测序仪 进行毛细管荧光电泳检测。结果整理。将毛细管电泳检测结果导入 分析软件中对原始数据进行条带分型，每对引物导出 峰图文件，结果汇总统计后进行分析。各位点等位基因数（）、有效等位基因数（）、观测杂合度（）、期望杂合度（）等遗传多样性参数使用 软件计算，各位点的多态信息含量指数（）由 软件计算，群体 平衡偏离用固定指数 衡量，使用 软件进行群体的哈代温伯格平衡检验。群体间的聚类分析运用 软件对其进行 方法构建系统发生树。群体间遗传分化系数（）、地方群体内近交系数（）、基因差异分化系数（）、分子方差、分析由 软件进行运算，热点图为 语言绘制。连锁不平衡使用 包
7、 的 方法进行分析，图用 语言绘制。结果与分析 各位点遗传多样性参数 由表 可知，采用的基因组微卫星位点均为多态性位点，个位点共检测到 个等位基因，大小在 。各位点观测等位基因数目（）在 个，平均数为 个；有效等位基因数（）在 ，平均值为。各位点检测到的香农信息指数（）在，平均值为。观测杂合度（）在，平均值为；期望杂合度（）在，平均值为。多态信息含量指数（）在，平均值为。各位点检测到的遗传分化指数（）值，平均值为。在 和 两个位点检测到的近交系数（），其他位点近交系数均。基因流（）为，平均值为。群体在、和 位点上显著偏离哈代温伯格平衡，群体在、位点上显著偏离哈代温伯格平衡，其中、位点在 个群体
8、中都显著偏离哈代温伯格平衡。表 个位点遗传多样性参数 位点等位基因数 有效等位基因数 香农信息指数观测杂合度 期望杂合度 固定指数 多态信息含量指数 遗传分化系数地方群体内近交系数 基因差异分化系数 基因流哈代温伯格平衡 平均值 注：表示无显著差异（）；表示显著差异（）；表示极显著差异（）；表示极显著差异（）。：（）；（），（），（）安徽农业科学 年 个群体各位点的遗传多样性参数由表 可知，群体检测到的观测等位基因数目平均值为，有效等位基因数平均值为，香农信息指数平均值为，观测杂合度平均值为，期望杂合度平均值为；群体检测到的观测等位基因数目平均值为，有效等位基因数平均值为，香农信息指数平均值为
9、，观测杂合度平均值为，期望杂合度平均值为。表 群体遗传多样性指数 群体 等位基因数 有效等位基因数 香农信息指数 观测杂合度 期望杂合度 群体遗传结构 和 个群体之间的 为，说明二者间存在中等程度的分化。个群体间的基因流（）为，说明二者间有一定程度基因交流。团头鲂群体的分子方差分析结果（表）显示，群体间的差异占差异来源的，个体间的差异较小，差异的主要来源在个体内部，占总差异的。个体内部的差异是指由杂合等位基因引起的遗传差异，大小与个体杂合位点数相关，即个体的遗传多样性，由此可知 个养殖群体的个体内部存在着丰富的遗传多样性。根据遗传距离绘制的群体间系统发生树（图），群体和 群体首先聚合在一起，然
10、后再与 群体聚合在一起，说明 群体和 群体之间的遗传距离更近，而 个养殖群体离 群体之间的遗传距离更远。表 分子变异方差分析 （）差异来源 自由度 估算的差异值差异值的百分比 群体间 个体间 个体内 合计 图 群体系统发生树 （）用来研究样本群体组成的相似性或相异性（图），群体的个体几乎集中分布在同一个区域，群体内的个体也是集中分布在另一个区域，而 的个体离 个集中分布的区域均较远，单独分布在一个区域。说明这 个群体之间存在一定程度的分化，个体之间有较明显的遗传差异。图 团头鲂群体的 分析 在连锁不平衡中，一般通过 值平方来表征连锁不平衡程度，个位点两两之间 值平方在（图），说明位点之间存在一
11、定程度的连锁不平衡。利用在线工具 进行群体结构分析，参数设置如下：混合模型（）；值为 ，值为；值设置为，每个 值重复 卷 期 张 芹等 基于基因组微卫星的团头鲂养殖群体遗传多样性分析和指纹图谱构建图 连锁不平衡程度示意 次。软件对每个 值模拟的结果，都会对应产生最大似然值（），最大似然值是取自然对数后输出的数值（）。越大，说明 值越接近于真实情况。一般随着 值升高，值也会不断升高，但会慢慢进入平台期，选择最优 值的目标是要找到拐点。该研究最佳 值为（图），根据最佳 值估计群体数为（图），说明试验所检测的样品最佳可划分为 个群体。利用 指纹图谱鉴别 个团头鲂群体 根据 个养殖群体在 个微卫星位点
12、上的扩增结果，按照表 中的编号顺序，利用 作图软件构建这 个群体的 指纹模式图（图）。从指纹图谱库中可以筛选出 个特异的微卫星 标记（、和）用以鉴别这 个群体，每个标记至少在其中一个群体中可以扩增出独有的条带，从而将这 个群体区分开。而将这 个标记配合使用，则可以将群体中的个体区分开来，这 个特异微卫星标记产生的多态位点可作为团头鲂育种过程中种质鉴定的依据。图 团头鲂群体遗传结构 值曲线图 （）图 值估计群体数 安徽农业科学 年图 微卫星 指纹模式图 讨论 群体遗传多样性分析 该研究观测等位基因数和有效等位基因数差别较大，则等位基因在 个群体中分布不均匀。团头鲂养殖群体检测到的平均等位基因数目
13、高于团头鲂二倍体和三倍体群体、鳜鱼、鲂鲌杂交及回交 群体、三角鲂翘嘴鲌、团头鲂翘嘴鲌 种杂交后代、团头鲂耐低氧新品系雌核发育群体和海蜇群体；观测等位基因数和有效等位基因数低于大口黑鲈和刺参群体，高于花鲈群体。等位基因情况反映出团头鲂 个养殖群体与其他水产养殖品种相比处于较高水平，一方面说明 个养殖群体保持了较高的遗传多样性，另一方面说明选择得到的基因组微卫星位点具有较高的遗传多样性，能够反映出团头鲂群体的遗传背景。该研究团头鲂养殖群体检测到的观测杂合度和期望杂合度均低于团头鲂群体和鳜鱼群体，高于三角鲂翘嘴鲌、团头鲂翘嘴鲌 种杂交后代、海蜇群体和花鲈群体。观测杂合度高于刺参群体，低于团头鲂耐低氧
14、新品系雌核发育群体；期望杂合度低于刺参群体，高于大口黑鲈群体和团头鲂耐低氧新品系雌核发育群体。观测杂合度和期望杂合度与鲂鲌杂交及回交 群体相比处于中等水平。综合来看，该研究团头鲂养殖群体的观测杂合度和期望杂合度与其他水产养殖品种相比处于中等水平。通过固定指数可以衡量种群中基因型的实际频率是否偏离 平衡，刺参和团头鲂耐低氧新品系中大量位点表现为杂合子过剩。群体在 个位点上显著偏离哈代温伯格平衡，群体在 个位点上显著偏离哈代温伯格平衡，说明人工选择育种过程给养殖群体造成了一定的选择压力。其中，个位点在 个群体中都显著偏离哈代温伯格平衡，可见这 个位点在不同的选育策略中都承受了较大的选择压力，可在今
15、后的选育过程中予以特别关注，以期进一步了解与位点相关的性状以及选育过程对其带来的影响。该研究中微卫星位点检测到的多态信息含量指数（）在，平均值为，高于鳜鱼群体、三角鲂翘嘴鲌、团头鲂翘嘴鲌 种杂交后代、团头鲂耐低氧新品系和花鲈群体，低于团头鲂群体和刺参群体，与鲂鲌杂交及回交群体接近。个团头鲂养殖群体的多态性较高，与等位基因衡量的结果一致。个微卫星位点中有 个位点的 值均大于，其中大于 的位点有 个，说明试验中采用的微卫星位点是高度多态的，能够为群体遗传信息监测提供丰富、高质量的数据。群体遗传结构分析 遗传分化指数（）是反映群体间遗传差异的重要参数，各位点检测到的遗传分化指数（）值，平均为；在检测
16、的 个位点中，个位点属于轻度分化位点（），个位点属于中度分化位点（），个位点属于高度分化位点（）。群体间的 为，说明 个群体属于中度分化群体，群体间遗传变异不显著。和 个群体间的基因流（）为，说明 个养殖群体之间有一定程度基因交流。基因差异分化系数（）可以不必考虑地方群体形成历史、迁移方式、群体中等位基因数和基因型频率，只需将整个群体中杂合体频率分解为地方群体间变异和地方群体内变异。各位点检测到的 范围在，群体间有一定的遗传分化。地方群体内近交系数（）对地方群体内的繁育系统具有指示意义，在 和 个微卫星位点检测到的近交系数，其他 个位点近交系数，个群体近交现象很少，个群体在 个位点的 值的范围在，接近于，说明群体比较接近于随机交配群体。各位点检测到的近交系数说明 个养殖群体在选育种过程中参与繁殖的亲本较多，群体保持较好的遗传多样性。基因流（）均大于，就能有效抑制由遗传漂变而引起的遗传分化。群体间基因流为，群体间存在一定程度基因流，个位点检测到的基因流（）大于，其中 个位点基因流大于。推测原因可能是伊河鲂群体选育时间较短，河南地区养殖的团头鲂苗种较多来源于湖北地区，养殖群体逃逸及放流的个
