牙鲆(Paralichthys olivaceus)是我国北方沿海地区重要的海水养殖鱼类之一。近年来，牙鲆的大规模工厂化养殖发展迅速，高密度的工厂化养殖和长期的近亲交配导致牙鲆种质退化、病害频发、产量下降，严重妨碍了牙鲆养殖业的可持续发展(王志敏等, 2005; 司飞等, 2017)。

迟缓爱德华氏菌病是由迟缓爱德华氏菌(Edwardsiella tarda)引起的水生动物疾病(陈爱平等, 2011)，俗称腹水病。该病具有感染率高、感染后死亡迅速、患病前期不易被发现等特点，25℃~30℃水温时该病最易发生。迟缓爱德华氏菌主要通过鱼类消化道、表皮受伤破损处、鳃进入鱼体，也可通过接触患病鱼体、摄入携带迟缓爱德华氏菌的食物或生物饵料而相互感染(陈爱平等, 2011)。李宁求等(2012)对爱德华氏菌菌蜕制备条件进行了优化，为鳗源迟缓爱德华氏菌病疫苗开发奠定了基础。陈昌福等(1998)对从日本鳗鲡(Anguilla japonica)分离出的迟缓爱德华氏菌进行致病力测试，注射菌液后的发病症状与自然发病的日本鳗鲡症状相同。

陈松林等(2008)、徐田军等(2010)、王磊等(2013)分别于2007、2009及2012年建立了大量牙鲆家系，通过鳗弧菌(Vibrio anguillarum)感染实验，已经选育出牙鲆抗鳗弧菌家系。陈松林等(2011)自2003年起，对中国牙鲆进行人工感染鳗弧菌，选育出中国本地的抗病群体，然后将中国本地的抗病群体与日本引进的牙鲆群体杂交，获得养殖成活率高、生长较快的群体，将其作为母本，以韩国引进牙鲆的选育群体作为父本进行杂交得到的三杂交后代，育成了牙鲆“鲆优1号”新品种。张英平等(2014)、孙何军等(2015)和郑卫卫等(2016)分别利用迟缓爱德华氏菌对牙鲆家系进行攻毒感染，筛选出了一些抗迟缓爱德华氏菌家系，为培育抗迟缓爱德华氏菌牙鲆新品种奠定基础。齐文山等(2014)研究表明，冷冻保存的异源精子可以诱导建立牙鲆卵裂雌核纯系，有助于构建牙鲆优良纯系。侯吉伦等(2014)以牙鲆两代雌核发育的个体为亲本，再度诱导雌核发育，获得三代雌核发育后代，结果显示，雌核发育后代个体的纯合度、平均遗传相似度进一步提高。

本研究主要通过家系选育方法，以2007、2009、2014年选育出的抗病牙鲆家系为亲本，于2016年培育出28个牙鲆家系，从中筛选出抗病速生牙鲆家系，并对这些家系进行了抗迟缓爱德华氏菌感染能力测试和比较分析。

亲鱼来源包括：日本牙鲆，韩国牙鲆，2009、2010、2012、2013及2014年牙鲆家系。2016年4月开始进行2016年抗迟缓爱德华氏菌牙鲆家系的培育。本次家系建立共使用雌鱼25尾，来自于16个不同家系；雄鱼27尾，来自于10个不同家系。家系命名规则参照王磊等(2013)，即F(family)+家系建立年份+家系编号。

2016年7月15日~9月15日，在山东莱州明波水产有限公司进行牙鲆家系生长性能测定实验，测定方法参照田永胜等(2009)。2016年7月15日，将鱼苗(69~90日龄)分别饲养在水体积为2.5 m³的圆形玻璃钢养殖桶中，每组300尾。实验期间，对鱼苗进行饱食投喂，流水养殖，保持水温22℃，24 h光照。60 d后，从所有牙鲆家系中各随机选取50尾鱼苗进行生长测量。测量指标主要为全长(从吻端到尾鳍末端)和体重。

通过预实验确定半致死浓度，再从每个家系中选取100尾牙鲆鱼苗(169~186日龄)，采用腹腔注射半致死浓度的菌液，注射剂量与预实验相同，并用PBS注射作对照组。攻毒感染实验的水温为(22±1)℃，注射后，每隔6 h检查并取样，及时清理死鱼，记录每条鱼的死亡时间、家系号，测量死鱼的全长、体重。感染后12 d，统计每个家系的存活鱼苗数量，测量其全长和体重。

用SPSS 19.0软件计算各家系的平均全长和体重，统计死亡与存活鱼苗数量，计算各个家系的存活率及各表型的相关性。

本研究于2016年成功建立和培育了28个牙鲆家系，包括F₃代家系4个、F₄代家系23个和对照家系1个(表 1)。对照家系来自山东海阳黄海水产有限公司同时期繁育的普通牙鲆鱼苗。亲鱼培育、苗种繁育参照陈松林等(2008)。

表 1(Tab. 1)

表 1 28个牙鲆家系亲本信息 Table 1 Parental information of 28 flounder families 家系Family 母本Female 父本Male 世代Generation F1601 JS F1406 F3 F1603 F1325 F1411 F4 F1606 F1342 F1422 F3 F1614 F1444 F1411 F4 F1615 F1419 F1421 F4 F1616 F1446 F1446 F4 F1621 JS F1317 F4 F1625 F0927 F1317 F4 F1627 F1013 F1317 F4 F1629 F1432 F1437 F3 F1631 F1263 F1421 F4 F1633 F1263 F1421 F4 F1634 F1413 F1421 F4 F1635 F1432 F1421 F4 F1636 F1432 F1421 F4 F1638 F1413 F1446 F4 F1639 F1430 F1421 F4 F1640 F1423 F1401 F4 F1645 F1002 F1446 F4 F1646 F1002 F1446 F4 F1647 F1260 F1417 F4 F1648 F1260 F1446 F4 F1650 F1260 F1446 F4 F1652 F1421 F1446 F4 F1655 F1432 F1421 F4 F1657 F1432 F1421 F4 F1661 F1430 F1413 F3 F1667 对照组Control group JS：日本群体(Japanese stocks)

表 1 28个牙鲆家系亲本信息 Table 1 Parental information of 28 flounder families

由于前期鱼苗较小，难以测量单个鱼苗的体重，故只进行了鱼苗体长测定。生长性能测定时间为60 d，测定时鱼苗较小，故筛选的速生家系属于早期生长快速的家系，对速生家系的筛选具有参考意义(表 2)。将全长日增长率大于0.160 cm/d的家系定为早期快速生长家系，共5个(F1639、F1621、F1661、F1606、F1648)，其平均全长日增长率为0.168 cm/d，占所有家系的17.9%；将全长日增长率在0.150~0.160 cm/d之间的家系定为早期生长速度一般家系，共9个，平均全长日增长率为0.152 cm/d，占所有家系的32.1%；将全长日增长率小于0.150 cm/d的14个家系定为早期慢速生长家系，平均全长日增长率为0.142 cm/d，占所有家系的50.0%(图 1)。

表 2(Tab. 2)

表 2 28个牙鲆家系的全长日增长率 Table 2 Daily length gain rate of 28 flounder families 家系 Family 初始全长 Primary length (cm) 最终全长 Final length (cm) 最大值 Maximum (cm) 最小值 Minimum (cm) 全长日增长率 Daily length gain rate (cm/d) 1639 2.10±0.34 12.54±1.01 16.0 10.8 0.174 1621 2.95±0.42 13.34±1.34 17.0 11.0 0.173 1661 2.55±0.39 12.48±1.20 15.6 10.2 0.166 1606 2.53±0.42 12.41±1.37 14.9 9.1 0.165 1648 2.03±0.28 11.80±0.85 13.3 9.6 0.163 1647 1.94±0.19 11.38±0.73 13.0 9.6 0.157 1646 2.61±0.27 12.01±0.98 14.0 9.7 0.157 1615 2.73±0.31 12.12±0.64 13.2 10.5 0.157 1638 2.46±0.45 11.76±0.81 13.3 9.7 0.155 1655 2.11±0.25 11.37±0.74 12.9 9.5 0.154 1634 2.33±0.21 11.50±0.87 13.0 8.7 0.153 1635 1.80±0.20 10.92±0.91 12.6 8.9 0.152 1633 2.29±0.42 11.34±1.13 15.6 8.7 0.151 1631 2.40±0.37 11.45±0.72 12.6 9.4 0.154 1640 2.44±0.39 11.37±1.25 14.5 8.6 0.149 1650 2.14±0.27 10.99±0.72 12.7 9.5 0.148 1603 2.51±0.32 11.34±0.94 13.6 9.7 0.147 1616 2.70±0.34 11.49±0.77 13.1 9.6 0.146 1652 2.46±0.28 11.18±0.86 14.8 9.4 0.145 1636 2.63±0.32 11.28±0.63 13.0 9.9 0.144 1601 3.42±0.66 12.05±1.05 15.2 10.3 0.144 1667 2.46±0.33 11.04±0.53 12.2 10.0 0.143 1614 2.43±0.42 10.97±0.81 12.6 9.1 0.143 1627 2.74±0.33 11.26±1.26 14.5 7.2 0.142 1657 2.42±0.38 10.92±0.96 13.2 8.6 0.142 1625 2.87±0.29 11.30±0.74 13.4 10.0 0.140 1645 2.57±0.27 10.91±0.74 12.2 9.0 0.139 1629 2.21±0.22 10.50±1.03 15.7 8.9 0.138 平均值Mean 2.44±0.47 11.54±1.11 17.0 7.2 0.151

表 2 28个牙鲆家系的全长日增长率 Table 2 Daily length gain rate of 28 flounder families

图 1(Fig. 1)

图 1 28个牙鲆家系的全长日增长率 Figure 1 Daily length gain rate of 28 flounder families

观察结果显示，注射24~72 h后，所有家系鱼摄食量减少，活力下降，出现白色粪便，部分鱼表现腹水症状。第3天开始出现死亡，死亡的个体体色发黑，腹部肿胀，部分个体有肠道外露现象。第4、5天开始大规模死亡，死亡持续到第12天逐渐稳定，总体平均存活率为(44.48±19.72)%(表 3)。其中，F1639、F1601、F1633、F1635、F1657、F1606和F1631 7个家系(占总量的25.0%)的存活率均高于60%，定义为高抗病力家系；11个存活率介于30%~60%之间的家系(占比39.3%)，定义为抗病力一般的家系；10个存活率低于30%的家系(占总量的35.7%)，称为低抗病力的家系(图 2)。

表 3(Tab. 3)

表 3 感染后各牙鲆家系存活率统计 Table 3 Statistic of survival rate of the infected flounder family 家系Family 感染鱼数/尾 Infected fish 存活率 Survival rate (%) F1639 101 77.23 F1601 97 74.23 F1633 42 73.81 F1635 99 70.71 F1657 98 69.39 F1606 83 67.47 F1631 96 66.67 F1634 50 54.00 F1667 97 53.61 F1652 105 53.33 F1638 90 52.22 F1621 96 52.08 F1614 92 51.09 F1640 100 44.00 平均Mean 93 44.48 F1646 74 41.89 F1615 102 39.22 F1655 97 39.18 F1648 94 36.17 F1647 105 27.62 F1625 97 26.80 F1650 100 26.00 F1616 98 25.51 F1661 102 24.51 F1627 96 23.96 F1636 82 21.95 F1603 101 18.81 F1629 102 18.63 F1645 105 15.24

表 3 感染后各牙鲆家系存活率统计 Table 3 Statistic of survival rate of the infected flounder family

图 1(Fig. 1)

图 2 感染后各牙鲆家系存活率统计 Figure 2 Statistic of survival rate of the infected flounder family

使用SPSS 19.0软件对28个牙鲆家系中生长性能与抗病性能进行皮尔森相关性分析(表 4)，P < 0.05说明存活率与生长性状显著相关，其中，全长与存活率的相关系数为0.496，全长日增长率与存活率的相关系数为0.301，体重与存活率的相关系数为0.432，说明存活率与全长、全长日增长率、体重之间呈现中等相关性。

表 4(Tab. 4)

表 4 牙鲆28个家系中各项表型间相关性分析 Table 4 Correlation analysis of phenotypes in 28 flounder families 全长 Body length 全长日增长率 Daily length gain rate 体重 Body weight 存活率 Survival rate 全长Body length 1.000 0.850** 0.974** 0.496** 全长日增长率Daily length gain rate 1.000 0.811** 0.301 体重Body weight 1.000 0.432* 存活率Survival rate 1.000 注：*在0.05水平上显著相关；**在0.01水平上显著相关(双尾检验) Note: * Indicates significant correlation at the 0.05 level, ** indicates significant correlation at the 0.01 level (2-tailed)

表 4 牙鲆28个家系中各项表型间相关性分析 Table 4 Correlation analysis of phenotypes in 28 flounder families

图 1(Fig. 1)

图 3 抗病家系系谱 Figure 3 The family trees of disease resistance family 虚线外框家系为抗病家系，箭头上端为亲本，←为杂交，KS：韩国群体；RS：抗病群体；JS：日本群体；LS：鲈鱼精子 The families in frame with dotted lines are the resistant families, the upper ends of the arrows are the parents of families arrowed, ←represent hybridization, KS: Korean stocks; RS: Disease-resistant stocks; JS: Japanese stocks; LS: Lateolabrax japonicus sperm

由图 1和图 2可见，F1639家系早期生长最快，抗病力也最强，为F₄代中的速生高抗病力家系。F1639的母本来自F1421家系，父本来自F1430家系。F1430家系为F₃代中筛选出的生长快速家系，其母本来自F0917家系，父本来自F0915家系，均为抗病速生家系；F1421家系为F₃代选育出的抗病速生家系，其母本来自F0927家系，父本来自F0915家系。F0927家系为F₂代筛选出的抗病力高家系，是F0750家系个体雌核发育后代。F0750为F₁代选育出的抗病力高家系，亲本为选育出的中国牙鲆抗病选育群体(RS)与日本引进后选育的牙鲆群体(JS)。

F1639家系的感染存活率为77.23%，比2016年构建的28个家系的平均存活率(44.48%)高32.75%。F1421家系为F₃代选育出的抗病力高家系，存活率为47.06%，比F₃代平均存活率(20.29%)高26.77%。F0927家系为F₂代筛选出的抗病力高家系，感染存活率为64.97%，比F₂代平均存活率(28.24%)高36.73%。F0750为F₁代选育出的抗病力高家系，感染存活率为64.05%，比F₁代平均存活率(32.60%)高31.45%。

牙鲆等多种水产养殖动物感染迟缓爱德华氏菌后会出现摄食量和活力明显下降，腹腔膨大，腹膜软化出血，内脏肿大甚至外漏，肝、肾等组织局灶性坏死等临床症状(张晓君等, 2005; 刘云等, 2000)，这些症状与本研究所观察到的结果基本一致。另外，比较人工腹腔注射感染与自然感染，人工注射感染的发病症状表现突出，发病情况较为普遍，不同家系之间的抗病力差异区分明显。

2016年本研究共建立了28个牙鲆家系，包括F₃代家系4个、F₄代家系23个和对照家系1个。通过迟缓爱德华氏菌攻毒感染实验，成功筛选出抗迟缓爱德华氏菌F₃家系2个(F1601和F1606)，F₄代家系5个(F1633、F1635、F1639、F1657和F631)。分析表明，2016年抗病力强的家系的父本来自2014年抗迟缓爱德华氏菌病家系(F1406、F1421和F1422)，其祖先或祖先之一67%来自抗鳗弧菌家系。由此可见，近年来本实验室不断选育的牙鲆抗病家系效果明显。在抗病力强的7个家系中，F1639家系的早期阶段的全长日增长率(0.174 cm/d)最高，为早期快速生长家系；F1606为早期快速生长家系，F1633和F1635为生长速度一般家系。综合比较这些家系的抗病和生长性能可以看出，F1639在抗病和生长性能方面都优于其他家系，可作为以后研究的F₄代重点家系。马爱军等(2014)从2010年开始大菱鲆抗鳗弧菌病选育研究成功选育出5个抗病力较强的优良家系，可作为抗鳗弧菌选育的核心育种群体。朱佳杰等(2012)构建了吉富罗非鱼(GIFT strain Oreochromis niloticus)家系，并进行无乳链球菌(Streptococcus agalactiae)攻毒感染实验，利用家系选育方法对吉富罗非鱼进行了抗病品种选育的研究。

为了说明抗病选育对各世代抗病力的提高情况，比较本实验室所培育的牙鲆F₁~F₄世代家系攻毒感染后的存活率。其中，F₁代为2007年陈松林等(2008)建立的抗鳗弧菌病家系，F₂代家系为2009年王磊等(2013)建立的抗鳗弧菌病家系，F₃代家系为2014年孙何军等(2015)建立的抗迟缓爱德华氏菌病家系。本实验建立的大多数为F₄代抗迟缓爱德华氏菌病家系。

F₁代与F₂代进行的是抗鳗弧菌选育，F₂代家系对鳗弧菌攻击的抗病力显著增强，抗病力高家系的攻毒感染存活率由F₁代的51.86%提高至F₂代的60.68%。F₃代与F₄代进行的是抗迟缓爱德华氏菌选育，抗病力高家系的攻毒感染存活率由F₃代的40.59%提高至F₄代的72.14%。在家系选育中，虽然家系培育占用场地较多、牙鲆亲鱼性成熟时间较长，但家系选育能够增加有效群体数量，控制近交系数的增量，延长世代间隔，减少基因流失(黄永春等, 2013)，因此，家系选育是进行鱼类抗病良种培育的重要技术手段。

经过F₁代、F₂代选育出的牙鲆家系具有生长迅速、养殖成活率较高、抗鳗弧菌感染的优点。本研究通过家系选育，建立了牙鲆抗迟缓爱德华氏菌病家系，成功筛选到抗病力强、生长快速的F₄代家系(F1639)，为进一步培育抗病高产优质的牙鲆优良品种奠定了重要基础。本研究表明，经过多代的家系选育和攻毒感染实验，牙鲆的抗病力得到逐代提升，实现了抗病力性状向后代的稳定遗传。本研究结果对于其他鱼类抗病良种培育，具有重大的应用价值和广阔的推广前景。