鱼类抗病育种研究进展

周欣，高风英，卢迈新*

(1.中国水产科学研究院珠江水产研究所 农业农村部热带亚热带水产资源利用与养殖重点实验室，广东 广州 510380; 2.上海海洋大学 水产与生命学院，上海 201306)

近年来，鱼类养殖业发展稳定，2019年中国水产品养殖产量达到5 079.07万t，鱼类占到2 708.61万t，其中，海水鱼类养殖产量160.58万t，淡水鱼类养殖产量2 548.03万t，养殖面积达711万hm2[1]。根据2018年联合国粮食及农业组织发布的报告，随着经济的增长，人们对优质蛋白质的选择趋势也增加，鱼产品作为蛋白质和必需微量营养素的重要来源之一，其发展的健康和稳定十分重要[2]，并且鱼类产品消费的益处远大于危害[3]。然而，鱼类养殖生产活动易受到多种因素的影响，且鱼类生存对水质、水温等生存环境存在着严格的要求，近年来，为了满足水产品日益增长的需求，多数地区的水产养殖模式从粗放型转变为集约型，因此，极易导致养殖动物疾病的暴发[4]，且不同的鱼类物种所面临的病害问题也不同，主要是细菌性和病毒性疾病。如由鳗弧菌Vibrioanguillarum引起的牙鲆Paralichthysolivaceus弧菌病[5]，这对分布在中国、日本和韩国海域的海水养殖业造成了巨大的损失[6-7]。半滑舌鳎Cynoglossussemilaevis是分布在中国沿海地区具有极高经济价值的名贵鱼类，其生长快、肉质鲜美，但近年来随着其养殖规模的扩大，细菌病害威胁越来越严重，其中，哈维氏弧菌Vibrioharveyi是侵入半滑舌鳎的主要细菌性病原之一，能够引起烂鳍等疾病[8-9]。此外，哈维氏弧菌也可诱发大黄鱼Larimichthyscrocea弧菌病[10]，更为严重的是若误食感染哈维氏弧菌的食物或水，会引发人腹泻或败血症、中耳炎等疾病[11]。罗非鱼Oreochromisspp.是联合国粮食与农业组织推荐的21世纪重要淡水养殖品种之一，是全球未来动物性蛋白质的主要来源之一[12]，但罗非鱼集约化养殖模式的大规模推广及养殖环境的恶化，使得病害暴发越来越严重，特别是无乳链球菌Streptococcusagalactiae感染的传播速度和引起的死亡率逐渐增加。2009年以来，在以色列、厄瓜多尔、埃及、泰国和印度等多个国家和地区的罗非鱼中暴发了一种由罗非鱼湖病毒(Tilapia Lake virus, TiLV)引起的病毒性疾病，被世界动物卫生组织列为一种新生病毒，对罗非鱼养殖业造成了严重威胁[13-14]。除此之外，鲫疱疹病毒(Carassiusauratusherpesvirus, CaHV)和草鱼Ctenopharyngodonidellus呼肠孤病毒(Grass carp reovirus, GCRV)等其他多种鱼类病害也频频暴发。因此，为保障鱼类养殖业的绿色健康发展，培育抗病力强的优良品种就显得十分重要。

1 鱼类抗病育种的起源

从病理学研究上来看，人们已经掌握了多种治疗或预防鱼类疾病的方法，如使用抗生素、化疗药物及接种疫苗等。但过度使用抗生素不仅会使细菌产生耐药性，抗生素进入到人体中也会对人造成伤害[11]；化疗药物的使用易造成环境污染及副作用发生；而疫苗接种及益生菌和免疫刺激剂的使用等相关方法在养殖实践中难以大规模实施[15]。因此，国内外专家开始研究鱼类固有的抗病遗传因素，并进行选择育种从而建立抗病品系[16]。与其他物种相比，鱼类的选择育种研究起源较晚。

国外第一次有文献记载的鱼类选择育种是1920年在美洲红点鲑Salvelinusfontinalis中进行的抗病新品系培育[17]。国内科学系统的鱼类选择育种起源较晚，20世纪70年代出现相关研究[18]。抗病育种相关研究起源于吴维新等[19]将鲤和草鱼杂交一代与草鱼回交获得三倍体草鱼型杂种，并经多年从杂交一代异缘四倍体中选育出可育的二倍体83-2系抗病草鱼[20]。

近年来，随着鱼类遗传学的不断创新发展，以及对鱼类基因组的探索，研究人员建立了一系列的育种技术，包括家系与群体选育、分子标记辅助选育和全基因组选育等。

2 国内鱼类家系及群体选育研究

家系选育是应用较为广泛的育种方式，也为分子育种奠定了基础。其原理是通过逐代提高目的基因的纯合度，建立优良性状家系，获得具有优良性状且稳定的新品种，实际上就是对基因型的选择。近年来，国内外研究人员采用此方法在多种鱼类的抗病育种工作中取得了重要进展，如牙鲆、虹鳟Oncorhynchusmykiss和草鱼等。通过加性遗传方差分析证明，对虹鳟病毒性出血性败血病抗病力进行选育的方法是可行的[21]。通过自然选择和人工感染迟缓爱德华氏菌Edwardsiellatarda，对63个牙鲆家系的生长和抗病性能进行初步测定，从中鉴定得到3个抗病力强的家系[22]；进一步以抗病牙鲆品种为亲本建立了57个家系，选择其中46个家系进行感染试验，从中鉴定出5个抗病家系[23]；从28个优质抗病牙鲆家系中筛选得到7个高抗病力(攻毒感染存活率>66%)的家系[24]；在建立的抗淋巴囊肿(Lymphocystis disease，LD)牙鲆品系中，通过在淋巴囊肿病高发养殖场进行天然攻毒试验，得到G2家系，以此为基础继续选育，有望得到抗LD牙鲆新品种[25]。另外，遗传评定模型分析表明，采用家系选育的方式适合对半滑舌鳎抗哈维氏弧菌性状进行选育[26]。在野生草鱼中已被证实抗GCRV能力较好的个体遗传力较高，为构建草鱼选择育种提供了基础[27]。

除此之外，采用上述类似的方法在其他鱼类选择育种中也取得了一定进展。许氏平鲉Sebastesschlegelii的遗传改良和优质新品系选育就是通过群体选育和家系选育相结合的方法进行的[28]。罗非鱼不同家系的抗链球菌能力存在显著性差异[29]；通过对42个吉富罗非鱼Oreochromisniloticus家系人工腹腔注射感染无乳链球菌，评估筛选出抗病力强且生长速度快的7个家系[30]；吉富罗非鱼经过2个世代的选育，抗病能力便得到了较大提高[31]。罗非鱼不同品种间的抗病力也不同[32-35]，这为新品种的选育提供了基础和材料。此外，镜鲤Cyprinuscarpiovar.specularis抗锦鲤疱疹病家系选育已经进行到了F3代[36]。

以上研究结果均表明，家系选育是培育鱼类抗病优良品种的有效途径，且已在多种鱼类选育中取得成功，证明其是培育抗病新品系的重要手段。

3 鱼类抗病选育的分子标记辅助研究

常规的鱼类育种选育策略是利用细菌或病毒对种群内的个体感染以筛选出抗病能力较强的个体留种选育。但因鱼类的抗病力不仅受遗传影响，而且受环境影响，表型变异不能完全或真实地反映遗传变异[37]，所以这种方法可信度不高，并且选育周期长。随着遗传学研究的发展，尤其是“数量遗传学理论”的提出，动物育种技术进入到快速发展时期。数量遗传学理论的观点是生物的若干性状(如抗病)是由许多微小的基因共同决定的，这些微小基因被称为数量性状位点(quantitative trait locus，QTL)。但若仅研究控制某一性状的多个微小基因，无法准确定位到每一个微小基因并准确估计其单个基因的效应[38]。分子标记技术是传统育种选育技术的辅助手段，即保留了家系选育的优点，也能够对于某些性状的QTL位点进行精确定位，筛选和培育优良品种。其基本原理是利用与目标基因紧密连锁或共分离的分子标记进行区域或全基因组筛选，从而减少连锁冗余，达到提高育种效率的目的。目前，在鱼类抗病育种中应用较多的分子标记辅助选育技术主要有微卫星标记和单核苷酸多态性(SNP)技术两种。

3.1 鱼类抗病相关微卫星标记的筛选及QTL定位

微卫星标记是均匀分布于真核生物基因组中的简单重复序列，由2～6个核苷酸的串联重复片段构成，由于重复单位的重复次数在个体间呈高度变异性且数量丰富，故该技术应用广泛。采用微卫星标记技术对不同家系间的遗传结构和差异进行分析，构建高密度遗传连锁图谱并对抗病基因进行精准的QTL定位，可为后续的抗病选育奠定基础。

在中国，由微卫星标记构建的鱼类高密度遗传连锁图谱首先在鲤Cyprinuscarpio[39]、半滑舌鳎[40]、牙鲆[41]、鲢Hypophthalmichthysmolitrix[42]和鳙Aristichthysnobilis[43]等几种重要经济鱼类中完成(表1)。目前，将与 LD 抗性相关的遗传连锁图谱连锁群LG15上的主要基因座LD-R作为标记辅助育种的候选微卫星基因座，成功获得牙鲆Paralichthysolivaceus抗病品系[6, 44]。采用基因分型微卫星识别方法已成功对大西洋鲑Salmosalar抗传染性鲑鱼贫血症进行QTL定位[45]。利用微卫星标记技术已成功对4种大菱鲆Scophthalmusmaximus品系的杀鲑气单胞菌Aeromonassalmonicida抗性进行QTL定位[46]。利用微卫星标记技术在吉富品系尼罗罗非鱼中筛选出2个与抗病性状相关的位点[47]。采用混合分离子分析法(bulked segregant analysis，BSA)对半滑舌鳎中169个微卫星标记进行筛选，在感染试验存活率为52.22%的家系中得到了1个可能与抗哈维氏弧菌相关的微卫星标记，并成功对其进行了QTL定位，这为选育抗病新品种提供了基础和依据[9]。目前，在多种鱼类中完成了对抗病基因的精准QTL定位(表2)。

表1 部分鱼类微卫星遗传连锁图谱构建情况

表2 部分鱼类抗病相关QTL定位

3.2 基于候选基因的鱼类抗病相关SNP标记筛选

SNP是第三代分子标记技术，能够建立起更加密集的SNP图谱，遗传稳定性高、位点丰富且分布广泛、易于基因分型，从而定位出与目的性状更加紧密相关的标记或QTL区间[66]。近年来，SNP标记筛选技术在鱼类育种中得到广泛应用(表3)。

表3 部分鱼类抗病相关SNP标记的筛选研究进展

评估SNP位点的多态性并采用Snapshot法进行分型验证，在鲤的11个基因中检测到20个可用于连锁分析的标记[67]。草鱼作为中国主要的淡水养殖鱼类，病害频发严重,缺乏良种，但草鱼的优良性状选育已经取得了重要的进展。目前，利用测序筛查SNP位点，并采用Snapshot进行分型验证的方法，在草鱼抗呼肠孤病毒基因Mx2、TLR3(Toll样受体基因3)、补体C6、RIG-Ⅰ(视黄酸(维甲酸)诱导基因蛋白Ⅰ)、LGP2(RIG-Ⅰ样受体家族成员遗传生理学实验室2,laboratory of genetics and physiology 2)、MAD5(黑色素瘤分化相关基因5， melanoma differentiation-associated gene 5)和IPS-Ⅰ(β-干扰素启动子刺激因子Ⅰ， interferon-b promoter stimulator Ⅰ)中发现多个对草鱼GCRV抗性或易感性性状显著相关的SNP位点及单倍型。这些SNP位点和单倍型可以作为潜在的分子标记辅助筛选草鱼抗出血病品系[68-75]。通过人工感染筛选出异育银鲫Carassiusauratusgibelio易感克隆系和抗性克隆系，进一步鉴定出8个干扰素(interferon, IFN)相关基因与抗CaHV的侵染有关，可为今后异育银鲫抗病育种提供科学依据[76-77]。利用测序筛查SNP位点并采用Snapshot进行分型验证的方法，在罗非鱼IPS-Ⅰ、NOD1(核苷酸结合和寡聚化结构域1,nucleotide binding and oligomerization domain 1)、β2m(β2-微球蛋白，β2-microglobulin)、IKaros(C2H2型锌指蛋白类转录因子)、LBP(脂多糖结合蛋白， lipopolysaccharide binding protein)和MCP-8(肥大细胞蛋白酶-8, mast cell protease-8)基因中筛选到与无乳链球菌抗性/易感相关单倍型或SNP位点[78-83]。这些SNP位点和单倍型可以作为潜在的分子标记辅助筛选罗非鱼抗链球菌病品系。

4 候选基因MHC多态性与抗病性的相关性

主要组织相容性复合体(major histocompatibility complex，MHC)是编码动物主要组织相容性抗原基因群的统称，分为MHCⅠ类和Ⅱ类分子，由其所编码的细胞表面转膜蛋白可以结合内外源抗原并呈递给T淋巴细胞，从而引发免疫反应。多基因性、高度多态性是与抗病力相关的MHC类基因的两个重要特性，目前公认其是研究与抗病/敏感表型相关联的多态性分子标记最合适的候选基因。Miretti等[87]2005年发表了涵盖人MHC基因区域4.46-Mb高分辨率(1.9 kb)的连锁不平衡图谱及第一代标签SNP，这使得基因分型效率增加了5倍，从而可用于未来与MHC相关的所有疾病研究中。由于鱼类MHC的高度多态性，其在抗病选育中备受关注，如在大西洋鲑、虹鳟和镜鲤等鱼类中发现，MHC位点具有高度多态性且与抗病力有关[88-89]。

4.1 鱼类MHCⅠ 类基因

1990年，研究者第一次扩增出鱼类的MHCⅠA链基因[90]，1993年，在大西洋鲑中测得第一条鱼类MHCⅠA链 cDNA全序列[91]，随后陆续报道了多种鱼类的MHCⅠ类基因cDNA序列。不同的硬骨鱼类中MHC基因的拷贝数不同，但其结构与哺乳动物相似。

4.2 鱼类MHC Ⅱ类基因

1993年，第一条鱼类MHCⅡA基因cDNA序列首次在斑马鱼中被报道[92]。根据目前所得到的鱼类MHCⅡAcDNA分子序列，推测其蛋白质结构与哺乳动物相似。MHCⅡB基因是MHC基因中被研究最多的，其具有丰富的多态性，如牙鲆MHCⅡB基因的多态性[93]。

4.3 MHC 基因多态性与鱼体抗病性的关系

1999年，Van Muiswinkel等[94]提出MHC与鱼类疾病抗性具有重要关系。之后，国内外对多种鱼体抗病性与MHC多态性的关系进行了研究(表4)。如急性病毒感染过程中，虹鳟MHCⅠ基因在病毒感染期间表达显著下调[89]。大西洋鲑MHCⅠ和MHCⅡ类基因的多态性也与抗贫血病毒和疖病有关[95]。

表4 MHC基因多态性与鱼体疾病抗性的关系

4.3.1MHCⅡA基因多态性与抗病性 荷包红鲤CyprinuscarpioLinnaeusMHCⅡA类基因具有较高的多态性，等位基因Cyca-DXA4和Cyca-DXA5与对嗜水链球菌的敏感性显著相关，而Cyca-DXA17、Cyca-DXA18和Cyca-DXA19与其对嗜水链球菌病抗性相关[96]。团头鲂MHCⅡA基因的多态性与抗细菌性败血症性状也具有显著关联[84]。在牙鲆60个家系的MHCⅡA和MHCⅡB基因中，筛选到与弧菌病抗性密切相关的抗病/敏感家系MHC等位基因，其可作为分子标记进行抗弧菌病牙鲆品系的选育[97-98]。通过新加坡石斑鱼虹彩病毒(Singapore grouper iridovirus，SGIV)侵染试验研究了MHCⅡA基因多态性与斜带石斑鱼Epinepheluscoioides抗病/敏感性的关系，结果得到了37个高多态性等位基因，经过验证分析发现，与SGIV敏感性和SGIV抗性相关的等位基因各一个[99]。在黄河鲤中发现多个维氏气单胞菌Aeromonasveronii易感性相关的等位基因，经验证MHCⅡ类基因在黄河鲤抵抗病原入侵的过程中有着不可或缺的重要作用[100]。大西洋鲑MHCⅡ基因多态性与对鲑疮痂鱼虱Lepeophtheirussalmonis的敏感性间具有相关性[101]。在虹鳟中采用聚合酶链反应-SSCP分析技术等方法检测MHCⅡA类DAA基因(尤其是外显子2)的多态性，并研究其与对传染性造血器官坏死病病毒(IHNV)抵抗力的关系，结果发现与抗性有关的两个等位基因Onmy-DAA*1301和Onmy-DAA*0304[102]。从鮸鱼Miichthysmiiuy中发现共有48个(Mimi-DAB)功能基因，并预测有14个MHCⅡB类位点与抗原结合位点相关；还鉴定出40个MHCⅡA类(Mimi-DAA)功能基因，预测有5个位点与肽抗原的结合有关[103-104]。尼罗罗非鱼个体中MHCⅡA类基因多态性与无乳链球菌感染的易感程度间的关联分析表明，等位基因orni-daa*1101与无乳链球菌的易感性有关；而MHCⅡA等位基因第2外显子可用于探索疾病易感性或抗性与MHCⅡA的多样性间的联系[105]。Wynne等[106]利用MHCⅡ类基因作为分子标记成功地筛选出了抗阿米巴鳃病(Amoebic gill disease, AGD)的大西洋鲑品系。

4.3.2MHCⅡB基因多态性与抗病性 在关于MHCⅡB基因的研究中，研究人员通过分析抗病和不抗病牙鲆的MHCⅡB基因多态性，筛选得到与抗病相关联的MHCⅡB等位基因[93]。采用RACE法成功克隆出赤点石斑鱼EpinephelusakaaraMHCⅡB基因的cDNA全序列，并利用单链构象多态性(SSCP)分析技术证明了赤点石斑鱼MHCⅡB的表达多态性[107]。对斜带石斑鱼MHCⅡB基因多态性也进行了研究，关联分析结果证明，等位基因EPCO-DBB*1001与SGIV抗性显著相关[108]。在对半滑舌鳎MHCⅡB基因多态性与抗病相关性的研究中，得到9个与鳗弧菌病抗性/敏感性相关的位点[109]。在虹鳟中鉴定出45个不同的MHCⅡB类等位基因，其中Onmy-DAB*0201、Onmy-DAB*0904和Onmy-DAB*1102与IHNV易感性相关[110]。在大菱鲆抗病性与MHCⅡB基因多态性的关系研究中发现,与抗病/敏感相关的MHCⅡB基因型各1个[111]。在“全红”体色瓯江彩鲤Cyprinuscarpiovar.color中也发现了与抗病/易感病群体相关的MHCⅡB基因型[112]。卵形鲳鯵TrachinotusovatusMHCⅡB 基因第2外显子多态性与美人鱼发光杆菌Photobacteriumdamsela抗性/敏感性有关，关联分析结果表明，1个等位基因与其抗性有关，3个等位基因与其敏感性有关[113]。根据尼罗罗非鱼F1代基因型和F3代链球菌抗性/易感性状进行关联分析，结果表明，MHCⅡB等位基因Orni-DAB*0107、Orni-DAB*0201和Orni-DAB*0302与链球菌抗性显著相关, 等位基因Orni-DAB*0701与链球菌易感显著相关[114]。在该研究之前，对尼罗罗非鱼的研究也表明，MHCⅡB等位基因Orni-DAB*0201与链球菌抗性显著相关[115]。因此，MHCⅡB等位基因Orni-DAB*0201可以作为尼罗罗非鱼抗链球菌病选育的分子标记，用于将来的抗病辅助育种。

5 全基因组抗病选择育种研究

2001年Meuwissen等[117]提出了全基因组选择的概念。用覆盖全基因组的分子标记进行辅助选育的手段，是对传统标记辅助选育的创新和改进，这已经成为动物遗传改良的研究热点及运用到育种生产中的重要手段。其主要方法是通过全基因组中大量的高密度遗传标记估计出不同染色体片段或单个标记效应值，然后将个体全基因组范围内片段或标记效应值累加，从而获得基因组估计育种值。目前，利用大黄鱼基因组育种值筛选抗性亲本取得了较好的效果[118]。所以全基因组测序及高密度遗传图谱的建立和遗传育种值的预测模型在全基因组选择技术中十分重要。

5.1 全基因组测序

国外对鱼类全基因组测序的研究起步较早，始于20世纪90年代中期。2002年研究者完成第一个全基因组测序的鱼类——红鳍东方鲀Takifugurubripes[119]。随着测序技术的发展，越来越多的鱼类全基因组被测序，如青鳉Oryziaslatipes[120]、斑马鱼Daniorerio[121]、虹鳟[122]和大西洋鲑[123]等28种鱼类。

国内对鱼类全基因组测序的研究起步较晚，2014年，中国完成了半滑舌鳎的全基因组测序和精细图谱的绘制[124]；同年，鲤的全基因组测序完成[125]。目前，研究人员已完成了包括4种滩涂鱼Scartelaoshistophorus、Boleophthalmuspectinirostris、Periophthalmodonschlosseri、Periophthalmusmagnuspinnatus[126]，以及大黄鱼[127]、草鱼[128]、牙鲆[129]、鞍带石斑鱼Epinepheluslanceolatus[130]等共23种鱼类的全基因组测序[131]。

5.2 基于SNP标记的精细遗传图谱构建

目前，多种鱼类的高密度遗传连锁图谱已构建完成，如在牙鲆中构建了12 712个SNP标记平均间隔为0.47 cM的高密度遗传连锁图谱，并定位得到9个与鳗弧菌抗性相关的QTL位点[85]。除此之外，团头鲂Megalobramaamblycephala[132]、中华鲟Acipensersinensis[133]、鳙[134]、鲤[135]等其他鱼类的高密度遗传连锁图谱也构建完成[131]。

5.3 全基因组抗病育种研究

5.3.1 基于高通量测序技术的全基因组SNP位点的分型 随着高通量SNP标记检测成本的不断降低，以及基因分型和测序技术的发展，全基因组选择技术开始被成功应用于育种实践当中。目前，基于2b-RAD高通量测序的基因型分型技术已建立并被广泛应用[136]。采用2b-RAD方法进行基于SNPs的基因分型，可以有效地捕获抗病力的遗传变异[137-138]。根据草鱼基因组SNP电子芯片对感染GCRV后的全同胞群体草鱼中存活个体和死亡个体的基因分型进行基因组重测序，从筛选得到的具有显著差异的1 822个基因型位点中获得9个抗病相关基因和10个显性等位基因位点，进一步采用高通量SNP分型技术筛选野生群体中携带抗病基因的草鱼个体，以此为亲本人工繁育出抗病群体[128]。降低基因分型成本的策略：1)选择极端表型(EP)或预选择SNP以构建候选基因的低密度标记，而且单标记分析是预选择SNP中节省候选基因分型成本的最可行方法；2)在基因组选择中，利用极端表型选择和单标记分析相结合的方法，显示出了较好的结果[139]。Wang等[140]也建立了一种高通量、低成本的基因型分析技术。

5.3.2 全基因组育种值的预测方法 在全基因组抗病育种技术研究过程中，研究人员开发了多种全基因组育种值预测方法(表5)。对虹鳟的研究发现，所有全基因组育种值的预测方法精度都高于系谱的BLUP(P-BLUP)方法，其中，贝叶斯C(BayesC)的精度相对提高最多[141-142]。在大黄鱼Larimichthyscrocea育种中，采用基于测序的基因分型构建NGS测序文库，得到约30 000个SNPs，根据拟合预测精度曲线估计基因组预测中特定性状的算法达到理想预测精度所需的训练规模，并成功评估了预选择SNPs的可行性[143]。目前，Dong等[144]开发了两种新的大黄鱼基因组育种值的计算策略：一种是针对非基于马尔科夫链-蒙特卡洛算法(MCMC)的快速混合正态分布(FMixP)策略，另一种是针对基于MCMC的MixP(MMixP)策略。该研究中使用4种方法BayesA、BayesCπ、MMixP和FMixP来比较预测结果，结果表明，对于模拟数据，BayesCπ、MMixP和FMixP的性能优于BayesA，FMixP在计算上比基于MCMC的方法快得多，但是所有方法对于大黄鱼基因组育种值的预测能力都非常相似，该研究结果可能被应用于基因组育种值预测中；另外，Dong等[145]还提出了一种计算速度明显快于BayesC的快速贝叶斯C(fBayesC)方法，在基因组育种值预测上是可行的；该研究者还研究了对SNP基因型进行标准化与否条件下先验超参数的两种策略，并分别计算了7种预测模型(岭回归BLUP、BayesA、BayesB、BayesCπ、快速BayesB(fBayesB)、FMixP和MMixP)在两种策略下的基因组育种值，提出了先进行SNP基因型标准化，再设定先验分布参数值的建议[146]。“鲆优2号”牙鲆新品种培育过程中，采用BayesCπ和GBLUP两种算法进行独立个体育种值估算，从而建立了抗病牙鲆品系的基因组选择育种技术，成功培育出“鲆优2号”牙鲆新品种[147]。对半滑舌鳎良种选育研究表明，通过BayesCπ模型进行全基因组选择育种会取得较好效果，为通过芯片分型对半滑舌鳎全基因组选择育种奠定了基础[148]。

表5 全基因组选择计算方法

5.3.3 基于全基因组关联分析筛选抗病相关SNP位点 全基因组关联分析(GWAS)是筛选抗病相关SNP位点的重要手段。在鲈抗病研究中，采用一种简化的高通量测序技术——限制性位点相关DNA测序方法得到了病毒性神经坏死鲈的全基因组SNP数据，并对其进行GWAS，检测了个体SNP与VNN抗性间的关联[150]。对大西洋鲑的育种研究表明，基因组预测精度比传统预测模型获得的基因组预测精度高18%[151]。对大西洋鲑进行GWAS，获得5个与抗AGD显著相关的SNP位点[152]。在半滑舌鳎抗病性研究中，通过GWAS、Fst和核苷酸多样性筛选结合的方法识别对抗病性重要的位点，在基因组重测序中鉴定的 1 016 774 个SNPs位点中检测到 33 个与抗病性显著关联的SNPs位点[153]。

5.3.4 基于GWAS的QTL定位 采用GWAS法进行QTL定位，可以为鱼类抗病提供有价值的候选基因。如采用GWAS法发现了4个与斑点叉尾鮰Ictaluruspunctatus的柱状病抗性相关的QTL[149]，并发现了9个与斑点叉尾鮰爱德华氏菌Edwardsiellaictalurid免疫性相关的QTL[154]。在虹鳟中进行GWAS,检测到10个与抗病相关的QTL[155]。在感染锦鲤疱疹病毒(Koi herpesvirus,KHV)镜鲤的重要QTL或SNP周围发现了许多与免疫相关的基因，包括tnfa、hif1a、galectin-8、rootletin和palladin，并利用公开的KHV感染鲤的转录组数据，获得了候选基因的表达模式，这为鲤对KHV抗性的遗传结构研究提供了许多候选基因[114]。在鲤的LG 44上鉴定出一个与KHV抗性相关的QTL[156]。对虹鳟30 060个SNP进行GWAS，检测到几个QTL可能与抗冷水细菌病有关[157]。对大西洋鲑进行GWAS，在第二染色体的QTL区域发现3个可能与抗AGD有关的候选基因[152]。在大西洋鲑Ssa27和Ssa12染色体上发现两个与心肌炎病毒(piscine myocarditis virus，PMCV)有关的QTL，其中包括多个免疫相关基因，如MHCⅡ类基因等[158]。

6 存在问题及展望

在世界范围内，水产养殖业的发展速度快于所有其他食用动物生产部门，但却面临着巨大的疾病感染风险[4]。鱼类病害种类多样、病因复杂且传播性强，解决病害问题并减少病害带来的损失对整个水产养殖行业具有重要意义。近年来，随着育种技术的不断创新发展，育种水平不断提高，传统选育与分子标记辅助选育技术的结合，以及全基因组选择模型的不断优化，动植物育种工作都取得了很大的进展，同时在鱼类育种研究中也取得了一定进步。

6.1 抗病育种面临的问题

1)鱼类基因组资源开发发展迅速，但是后基因组时代，鱼类基因功能研究包括抗病相关基因功能研究相对滞后。影响鱼类基因功能研究的主要因素是相对于哺乳类，鱼类工具细胞系的建立相对滞后，目前建立的细胞系普遍存在转染效率低下的现象，严重影响了鱼类基因功能的研究、验证。另外，尚未建立以CRISPR/Cas9等为代表的基因组编辑技术和平台，严重影响了重要抗病性状关键基因的发掘及其在抗病育种中的应用。

2)鱼类抗病育种相关的基础研究及抗病分子育种的研究相对滞后。目前，对鱼类免疫系统的免疫机制研究还处于初级阶段。对鱼类免疫系统及其调控机制进行深入系统的研究，是推动鱼类抗病育种的基础，而抗病分子育种的研究是抗病育种的重要手段。

3)鱼类抗病育种的基因组选择研究尚处于初级阶段。价格低廉的分子标记筛选高通量分型平台尚未建立，鱼类抗病性状的界定手段还无统一的标准，全基因组育种值的统计方法还处于摸索阶段。

6.2 未来重点研究方向

1)鱼类免疫系统的解析。由于鱼类免疫知识的缺乏，限制了免疫系统进化研究、疫苗研发及抗病品种(系)的选育，因此，将来需要加强鱼类免疫系统的解析研究，包括鱼类先天免疫和获得性免疫系统的研究，同时也包括了免疫负调控因子的挖掘及其在先天免疫中作用机制的解析。

2)鱼类全基因结构解析及基因资源的深度挖掘。包括全基因组测序和精细图谱构建，鱼类基因组的拼接、组装和生物信息学分析技术，阐释不同鱼类基因组的结构特征、进化规律和功能，进行基因资源的深度挖掘，包括参考基因组序列及其注释(编码和非编码调节元件)、全基因组多态性标记、高效基因分型平台、高密度和高分辨率连锁图及转录组资源(包括非编码转录本)。

3)加强鱼类抗病性状的遗传解析及分子育种的基础研究。借助基因组编辑平台和鱼类工具细胞系培养技术，筛选抗病性状关键基因，研究和验证基因的具体功能和调控机制，阐明基因对抗病性状的决定机制，解析抗病性状基因的调控网络，为分子设计育种奠定理论基础；深入开展抗病性状的GWAS和基因组选择及分子设计的基础研究，包括全基因组育种值的统计方法研究、分子标记筛选的高通量分型平台技术的建立(如抗病育种芯片)、鱼类抗病性状界定方法标准的建立；拓展抗病性状相关分子标记和关键基因的筛选与鉴定研究，建立基因与抗病性状的关联性；定位抗病性状基因的QTL，明确基因型与抗病表型的关系。

4)确定抗病性状遗传和表观遗传调控。通过基因组甲基化分析、组蛋白修饰及小核糖核酸(RNA)分析等技术手段，解析环境因素介导的抗病性状形成的表观遗传调控机制；建立甲基化修饰因子在水产动物中的检测技术体系，筛选与抗病性状相关的表观遗传标记。