CN112273291B - 基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法 - Google Patents
基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法 Download PDFInfo
- Publication number
- CN112273291B CN112273291B CN202011171385.XA CN202011171385A CN112273291B CN 112273291 B CN112273291 B CN 112273291B CN 202011171385 A CN202011171385 A CN 202011171385A CN 112273291 B CN112273291 B CN 112273291B
- Authority
- CN
- China
- Prior art keywords
- breeding
- cryptocaryon irritans
- resistance
- gebv
- large yellow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000009395 breeding Methods 0.000 title claims abstract description 102
- 241001663423 Cryptocaryon irritans Species 0.000 title claims abstract description 55
- 241001596950 Larimichthys crocea Species 0.000 title claims abstract description 41
- 201000010099 disease Diseases 0.000 title claims abstract description 16
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 title claims abstract description 16
- 230000001488 breeding effect Effects 0.000 claims abstract description 95
- 241000251468 Actinopterygii Species 0.000 claims abstract description 62
- 238000000034 method Methods 0.000 claims abstract description 21
- 241000238631 Hexapoda Species 0.000 claims abstract description 14
- 238000003205 genotyping method Methods 0.000 claims abstract description 10
- 230000003042 antagnostic effect Effects 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 238000012795 verification Methods 0.000 claims abstract description 7
- 238000002474 experimental method Methods 0.000 claims description 16
- 238000012163 sequencing technique Methods 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 10
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 claims description 9
- 230000004083 survival effect Effects 0.000 claims description 9
- 239000003550 marker Substances 0.000 claims description 8
- 231100000614 poison Toxicity 0.000 claims description 8
- 239000002574 poison Substances 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 238000012258 culturing Methods 0.000 claims description 6
- 231100000419 toxicity Toxicity 0.000 claims description 6
- 230000001988 toxicity Effects 0.000 claims description 6
- 108090000623 proteins and genes Proteins 0.000 claims description 5
- 238000002790 cross-validation Methods 0.000 claims description 4
- 238000003908 quality control method Methods 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 4
- 238000012098 association analyses Methods 0.000 claims description 3
- 238000004140 cleaning Methods 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 230000012447 hatching Effects 0.000 claims description 3
- 210000003205 muscle Anatomy 0.000 claims description 3
- 238000002205 phenol-chloroform extraction Methods 0.000 claims description 3
- 238000009394 selective breeding Methods 0.000 claims description 3
- 238000012268 genome sequencing Methods 0.000 claims description 2
- 238000011451 sequencing strategy Methods 0.000 claims description 2
- 238000010200 validation analysis Methods 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims 1
- 208000031295 Animal disease Diseases 0.000 abstract description 3
- 238000012360 testing method Methods 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 241000282326 Felis catus Species 0.000 description 8
- WTDRDQBEARUVNC-LURJTMIESA-N L-DOPA Chemical compound OC(=O)[C@@H](N)CC1=CC=C(O)C(O)=C1 WTDRDQBEARUVNC-LURJTMIESA-N 0.000 description 5
- 208000031513 cyst Diseases 0.000 description 5
- 238000012214 genetic breeding Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 241000700605 Viruses Species 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 230000002068 genetic effect Effects 0.000 description 3
- VWDWKYIASSYTQR-UHFFFAOYSA-N sodium nitrate Chemical compound [Na+].[O-][N+]([O-])=O VWDWKYIASSYTQR-UHFFFAOYSA-N 0.000 description 3
- 238000007400 DNA extraction Methods 0.000 description 2
- 241001596954 Larimichthys Species 0.000 description 2
- 101150064138 MAP1 gene Proteins 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 2
- 241000276618 Perciformes Species 0.000 description 2
- 210000000349 chromosome Anatomy 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000011534 incubation Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010079 rubber tapping Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 108700028369 Alleles Proteins 0.000 description 1
- 206010011732 Cyst Diseases 0.000 description 1
- SDKQRNRRDYRQKY-UHFFFAOYSA-N Dioxacarb Chemical compound CNC(=O)OC1=CC=CC=C1C1OCCO1 SDKQRNRRDYRQKY-UHFFFAOYSA-N 0.000 description 1
- 208000035240 Disease Resistance Diseases 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 108010046335 Ferredoxin-NADP Reductase Proteins 0.000 description 1
- 101710168674 Ferredoxin-NADP reductase 1 Proteins 0.000 description 1
- 101710182458 Ferredoxin-NADP reductase, leaf isozyme 1, chloroplastic Proteins 0.000 description 1
- 206010027982 Morphoea Diseases 0.000 description 1
- 235000017517 Tephrosia purpurea Nutrition 0.000 description 1
- 240000004429 Tephrosia purpurea Species 0.000 description 1
- 238000009360 aquaculture Methods 0.000 description 1
- 244000144974 aquaculture Species 0.000 description 1
- 230000003115 biocidal effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 239000003053 toxin Substances 0.000 description 1
- 231100000765 toxin Toxicity 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
- A01K61/00—Culture of aquatic animals
- A01K61/10—Culture of aquatic animals of fish
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
- A01K67/00—Rearing or breeding animals, not otherwise provided for; New or modified breeds of animals
- A01K67/027—New or modified breeds of vertebrates
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
- A01K2207/00—Modified animals
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
- A01K2227/00—Animals characterised by species
- A01K2227/40—Fish
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/124—Animal traits, i.e. production traits, including athletic performance or the like
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/80—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in fisheries management
- Y02A40/81—Aquaculture, e.g. of fish
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Environmental Sciences (AREA)
- Zoology (AREA)
- Organic Chemistry (AREA)
- Analytical Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biodiversity & Conservation Biology (AREA)
- Animal Husbandry (AREA)
- Wood Science & Technology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Genetics & Genomics (AREA)
- Biophysics (AREA)
- Biochemistry (AREA)
- Biotechnology (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Physics & Mathematics (AREA)
- Marine Sciences & Fisheries (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法,属于动物抗病育种领域。方法:建立参考群体;对参考群体群体进行抗性性状测量;建立选育群体;对参考群体和选育群体进行基因分型;建立最佳全基因组选择模型;对选育群体的基因组育种值GEBV进行估算;以一定的选择强度选取选育群体GEBV值排在前面的个体作为亲鱼,并繁殖产生抗刺激隐核虫子一代;对抗性子一代进行刺激隐核虫攻毒验证。抗虫性显著提升。提高了对选育群体基因组育种值的估计准确性,大大缩短育种年限,在一代之内就可获得高抗性的后代,减少了刺激隐核虫对大黄鱼养殖产业造成的经济损失,同时也为其他鱼类的抗病育种提供借鉴和基础,具有广阔的应用前景。
Description
技术领域
本发明属于动物抗病育种领域,尤其是涉及利用全基因组选择技术进行大黄鱼抗刺激隐核虫培育的一种基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法。
背景技术
大黄鱼(Pseudosciaena crocea),隶属于鲈形目(Perciformes),石首鱼科(Sciaenidae),黄鱼属(Larimichthys),主要分布于太平洋西南海域,即中国东南沿海区域。大黄鱼的野生资源由于过度捕捞趋于枯竭,但在科技工作者突破大黄鱼的人工繁殖后,其养殖规模逐年增加,现以成为中国东南沿海各省,尤其是福建省的主要养殖的海水品种,也是中国海水养殖单产量最高的养殖品种。近些年来,由于大黄鱼养殖产业膨胀式的发展,导致沿海内湾网箱养殖区域的水质极度恶化,导致病害频发。在诸多病害中,刺激隐核虫造成的危害最大,俗称白点病,每年的春夏转交之际,由于水温突升,导致其爆发,给大黄鱼养殖产业造成了巨大的经济损失。至今,没有有效的药物可以治疗白点病,而通过遗传选育提升大黄鱼刺激隐核虫抗性便具有极大的实践意义,也避免了药物治疗对水域的进一步污染。
遗传育种的发展经过了漫长了历史,育种技术在不断的革新,使其能够满足不同需求的选育目的。现代育种技术的代表作之一,而且还在蓬勃发展的全基因组选择技术(GS),自从2002年提出以来,经过快20年的发展,已经在动植物遗传选育研究中得到了广泛的应用,使一些遗传力较低,较难测量的性状的选育成为了可能,例如抗病抗逆等性状。在水产动物的抗病遗传选育领域中,很少存在有将GS技术成功应用的案例,而仅仅局限于理论研究。在此技术背景下,本发明聚焦大黄鱼抗刺激隐核虫遗传改良,并将GS技术应用其中,形成一套独特的培育流程,实现了理论成果向实践应用转化。本发明需要大量的测序任务,虽然成本较高,但是随着测序技术的更新换代,测序成本在急速的下降,导致个体分型成本也在下降,这也使本发明在实践过程中,成本可进一步降低。相比于本发明的资金投入,其减少的刺激隐核虫的对养殖产业造成的经济损失是不可估量的。而且,本发明是一个完整的体系,操作清晰,实用性强,可以在水产动物抗病育种工作中进行推广。
发明内容
本发明的目的在于提供有效解决复杂性状选育中亲本无法进行抗性测评等问题的一种基于全基因组选择的大黄鱼抗刺激隐核虫病优良品系的选育方法。
本发明包括以下步骤:
1)建立参考群体;
2)对参考群体进行抗性性状测量;
3)建立选育群体;
4)对参考群体和选育群体进行基因分型;
5)建立最佳全基因组选择模型;
6)对选育群体的基因组育种值(GEBV)进行估算;
7)以一定的选择强度选取选育群体GEBV值排在前面的个体作为亲鱼,并繁殖产生抗刺激隐核虫子一代;
8)对抗性子一代进行刺激隐核虫攻毒验证。
在步骤1)中,所述建立参考群体的具体步骤可为:
①选取来源广泛的大黄鱼2龄健康性成熟亲鱼至少500尾,保证其基因多样性;
②亲鱼不经过刺激隐核虫攻毒挑选,强化培育后,催产并繁殖后代;
③后代即为参考群体,按照商品鱼培育体系进行参考群体的养殖。
在步骤2)中,所述对参考群体进行抗性性状测量的具体步骤可为:
①选取参考群体6月龄幼鱼至少2000尾,将这批鱼进行岸上养殖室内暂养,按照正常的室内养殖规章进行管理,暂养时间为两周;
②培育刺激隐核虫幼虫,并在孵化高峰期计数;
③从暂养的上述大黄鱼个体中,选取至少1000尾进行攻毒实验,分为四组,每组放入3t规格的实验桶,攻毒浓度为每5升水每尾鱼5000个刺激隐核虫幼虫;
④自攻毒起,开始不间断的跟踪每组鱼的攻毒状况;将失去平衡,且10s内不能恢复的个体视为死亡样本,记录死亡时间,并将其捞出实验桶外,剪其背部肌肉保存于纯酒精中;观察时间持续120h,实验结束后将所有存活的个体转移至新的实验桶,进行攻毒后恢复;恢复两周后,仍然存活的个体视为存活样本;如此,大黄鱼抗刺激隐核虫的表型性状采集模式为二分类性状,存活-死亡。
在步骤3)中,所述建立选育群体的具体步骤可为:
①选取即将性成熟且健康的2龄大黄鱼至少500尾作为选育群体;
②对上述选育群体打入电子标签,并剪取鳍条保存于酒精;
③将其转移入新的网箱进行产前强化培养。
在步骤4)中,所述对参考群体和选育群体进行基因分型的具体步骤可为:
①选取保存在酒精中的参考群体死亡样本和存活样本中的鳍条至少各200尾,选育群体样本至少350尾进行DNA提取,方法为酚氯仿法;
②对DNA样品进行基因分型,可以采取任何形式的测序策略,包括简化基因组测序,全基因组重测序,基因芯片等;
③将原始分型数据经过进一步清洗和质控,最终获得可以进行全基因组选择育种分析的可用分型数据集。
在步骤5)中,所述建立最佳全基因组选择模型的具体步骤可为:
①基于参考群体的二分类抗性表型和基因型进行全基因组关联分析,所使用软件为R语言环境下的gaston包,得到所有SNP位点的显著性P值,并将其由小到大排序;
②使用R语言环境下BGLR包中的Bayesian LASSO分析模型建立全基因组选择模型,以10个标记为增量,按照已经排好序的SNP标记,逐渐增加排在前面的标记使用数量,并利用5折交叉验证计算每个标记量下的模型预测准确性;
③判断模型预测准确性最高时,对应的标记数量,相应的此标记密度下的建立的模型即是最佳的全基因组选择模型。
在步骤6)中,所述对选育群体的基因组育种值(GEBV)进行估算的具体步骤可为:
①在最佳全基因组选择模型下,对选育群体的GEBV进行估计,并对其进行排名;
②利用人类中应用十分广泛的多基因风险评分来佐证上述得到的GEBV值的准确性,所使用的软件为PRSice软件。
在步骤7)中,所述以一定的选择强度选取选育群体GEBV值排在前面的个体作为亲鱼,并繁殖产生抗刺激隐核虫子一代的具体步骤可为:
①按照一定的选择强度,选择选育群体GEBV排在前面的个体作为亲鱼繁殖后代,如此,这些亲鱼的后代即为抗性子一代;
②按照正常的商品鱼养殖体系培育抗性子一代。
在步骤8)中,所述对抗性子一代进行刺激隐核虫攻毒验证的具体步骤可为:
①在刺激隐核虫爆发的季节,对抗虫子一代进行攻毒实验,实验流程同参考群体的抗性性状测评的攻毒流程一致;
②与对照组相比,评价其抗虫性。
本发明提供了一种基于全基因组选择的大黄鱼抗刺激隐核虫病优良品系的选育方法,利用本发明提供的培育流程,可直接将其应用到大黄鱼抗刺激隐核虫的选育工作中,也为后续培育稳定的抗虫品系提供最初的原始群体,即抗性子一代。与现有技术相比,本发明至少有以下优点及技术效果:
(1)本发明首次将GS技术应用到大黄鱼抗刺激隐核虫遗传育种中,有效的解决了复杂性状选育中亲本无法进行抗性测评的问题。
(2)本发明在两年之内就可获得抗虫子一代,大大缩短了育种年限。
(3)本发明得到的抗虫子一代的抗虫性能显著高于对照组。
(4)本发明可大大降低刺激隐核虫对大黄鱼养殖产业造成的经济损失。
附图说明
图1为本发明在具体实施过程中的示意图。
图2为本发明在具体实施过程中第⑵个步骤中参考群体抗性性状测量中死亡个体的死亡时间分布图。
图3为本发明在具体实施过程中第⑸个步骤中判断全基因组模型最佳标记的结果图。横坐标表示标记数量,图中上方的图中纵坐标代表模型拟合度(Goodness of Fit),下方的图中纵坐标表示模型预测准确性(predictability)。
图4为本发明在具体实施过程中第⑹个步骤中使用PRSice软件对R语言BGLR包得到的GEBV结果的佐证结果。横坐标表示使用BGLR包得到的选育群体的死亡概率,纵坐标表示使用PRSice软件得到的选育群体多基因死亡风险评分。
图5为本发明在具体实施过程中第⑻个步骤中抗虫子一代的刺激隐核虫攻毒验证结果。横坐标代表死亡时间,纵坐标代表生存概率,CS代表对照组,RS代表抗性子一代。
具体实施方式
下面将结合附图对本发明作进一步的说明和详细描述。
参见图1~5,本发明实施例包括以下步骤:
⑴建立参考群体;
①2017年12月份在福建某海域网箱养殖区,从各鱼排上选取来源广泛的大黄鱼2龄健康性成熟亲鱼500尾,亲鱼来源广泛保证其基因多样性。
②亲鱼不经过刺激隐核虫攻毒挑选,强化培育后,催产并繁殖后代。
③这批亲鱼的后代即为参考群体,按照商品鱼培育体系进行参考群体的养殖。
⑵对参考群体进行抗性性状测量;
①2018年6月中旬选取参考群体6月龄幼鱼4000尾作为攻毒实验材料,鱼的体重规格为6.80±2.52g,体长规格为7.8±1.1cm。将这批鱼进行岸上养殖室内暂养,按照正常的室内养殖规章进行管理,暂养时间为两周。
②培育刺激隐核虫幼虫,具体流程为:收集网箱养殖患病的大黄鱼,将其转移至岸上的养殖桶内;待其包囊脱落并粘附于桶底后,用刷子将包囊与桶底分离,并将其收集并清洗之后转移到1000mL大烧杯中,杯中注满过滤海水,并将其放置于避光的恒温水浴锅中孵化,温度调节为26℃;期间每隔2h观察包囊孵化情况,大约2天就可达到孵化高峰期,此时,利用显微镜计数孵化的幼虫密度。
③从暂养的实验鱼中,选取2000尾进行攻毒实验,分为四组,攻毒实验前一天将其放入3t规格的实验桶,并停止投喂。攻毒时,将计算好的携带有幼虫的培养液倒入实验桶内,并搅拌均匀,攻毒浓度为每5升水每尾鱼5000个刺激隐核虫幼虫。
④自攻毒起,开始不间断的跟踪每组鱼的攻毒状况。将失去平衡,且10s内不能恢复的个体视为死亡样本,记录死亡时间,并将其捞出实验桶外,剪其背部肌肉保存于纯酒精中。观察时间持续120h,实验结束后将所有存活的个体转移至新的实验桶,进行攻毒后恢复。恢复两周后,仍然存活的个体视为存活样本。如此,大黄鱼抗刺激隐核虫的表型性状采集模式为二分类性状,存活-死亡。
⑶建立选育群体;
①2018年12月份中旬,在福建某海域网箱养殖区选取即将性成熟且体表无伤口健康的2龄大黄鱼500尾作为选育群体。
②对上述选育群体打入电子标签,并剪取鳍条保存于酒精。
③进行完上述工作后,将其转移入新的网箱进行产前强化培养,而这些亲鱼即是抗虫子一代的候选亲鱼群体。
⑷对参考群体和选育群体进行基因分型;
①选取保存在酒精中的参考群体死亡样本和存活样本中的鳍条各240尾,还有384尾选育群体样本进行DNA提取,方法为酚氯仿法。
②对DNA样品进行建库测序,建库方式为ddRAD,测序平台为HiSeq2000。
③将原始测序数据经过清洗,获得纯净测序数据,利用stacks2软件将所有混在一起的个体基于不同的barcode分开;之后利用bwa程序将每个个体的测序clean reads与大黄鱼基因组(PRJNA505758)比对,之后利用samtools软件统计测序深度、基因组覆盖率和基因组比对率;之后使用stacks2软件对所有个体进行分型工作,生成vcf文件;之后利用plink软件对分型后的SNP位点进行质控,质控条件为SNP位点检出率>0.9(geno=0.1),个体SNP位点检出率>0.9(mind=0.1),SNP最小等位基因频率>0.05(maf=0.05);之后利用beagle4软件进行缺失基因型填补;最后使用haploview软件搜索tag SNP(高度连锁区域的代表SNP),最终获得可以进行全基因组选择育种分析的可用分型数据集。可用分型数据集包括16515个SNP,以及811个个体,包括451个参考群体个体以及360个选育群体个体。这部分工作在Linux系统下运行,以下是这部分工作的shell执行脚本。
##使用stacks2软件进行分个体工作
pre="path"
n=0
cmd="#!/bin/bash\n"
for i in`ls$pre|grep FKD`
do
f=$(echo$i|awk'{x=substr($1,1,4);y=substr($1,26,1);print x""y}')
cmd=$cmd"process_radtags-P-p$pre/$i/-b$pre/clean.data/adaptor.list/${f}_new_new1-o$pre/clean.data/140bp/--inline_null-c-r-q--renz_1ecoRI-t 130--barcode_dist_2 1-i gzfastq\n"
n=$(($n+1))
if[[$n-ge 16]]
then
echo-e$cmd|qsub-N${f}.tag-o$pre/clean.data/log/rad.tag.log/${f}.log-joe-l nodes=2:ppn=4-q high-
cmd="#!/bin/bash\n"
n=0
fi
done
##利用bwa程序与大黄鱼基因组比对并形成bam文件,利用samtools软件生成.sort和.depth文件
pre="path"
n=0
cmd="#!/bin/bash\n"
for i in`ls$pre/130bp|grep rem-v|grep fq|sed's/\..\.fq\.gz//g'|grepfq-v|sort-k 1,1-u`
do
cmd=$cmd\
bwa aln LG_Chr_genome.fasta$pre/130bp/${i}.1.fq.gz-k 2-t30>$pre/130bp/${i}_read1.sai\nbwa aln LG_Chr_genome.fasta$pre/130bp/${i}.2.fq.gz-k 2-t 30>$pre/130bp/${i}_read2.sai\n bwa sampeLG_Chr_genome.fasta$pre/130bp/${i}_read1.sai$pre/130bp/${i}_read2.sai$pre/130bp/${i}.1.fq.gz$pre/130bp/${i}.2.fq.gz|samtools view-bS->$pre/bam/${i}.bam\n cmd=$cmd"samtools sort-o$pre/sort/${i}.sort$pre/bam/${i}.bam\n/public/tools/sam_utils/samtools-1.8/samtoolsdepth$pre/sort/${i}.sort>$pre/depth/${i}.depth\n"
n=$(($n+1))
if[[$n-ge 49]]
then
echo-e$cmd|qsub-N${i}.process-o$pre/log/process.log/${i}.log-j oe-lnodes=2:ppn=4-q high-
cmd="#!/bin/bash\n"
n=0
fi
done
##使用samtools软件计算基因组比对率、测序深度以及基因组覆盖率
pre="path"
f=""
for i in`ls$pre/bam|sed's/\.bam//g'`
do
f=$f"$i\n"
echo$f
samtools-1.8/samtools flagstat$pre/bam/${i}.bam>>$pre/stat.summary.total/mapping.raw
done
grep"mapped("$pre/stat.summary.total/mapping.raw|awk'{print$5}'|sed's/(//g'>$pre/stat.summary.total/txt
echo-e$f|paste-$pre/stat.summary.total/txt>$pre/stat.summary.total/mappingratio
size=`grep-v">"LG_Chr_genome.fasta|tr"\n""@"|sed's/@//g'|wc|awk'{print$3}'`
for i in`ls$pre/depth|sed's/\.depth//g'`
do
awk-v a=$size-v id=$i'BEGIN{sum=0}{sum+=$3}END{print id"\t"NR/a"\t"sum/NR}'$pre/depth/${i}.depth>>$pre/stat.summary.total/depth
done
##使用stacks2软件进行基因分型工作,并生成vcf文件
pre="path"
stacks2/bin/ref_map.pl-T 6-o$pre/stacks--popmap$pre/popmap--samples$pre/
stacks2/bin/populations-O$pre/population--vcf-M$pre/popmap-P$pre/stacks##利用plink软件基因分型后的数据质控
pre="path"
vcftools--vcf$pre/populations.snps.vcf--plink--out$pre/lc.ci.final
plink--file$pre/lc.ci.final--allow-extra-chr-geno 0.10--recode--out$pre/lc.ci.final
plink--file$pre/lc.ci.final--allow-extra-chr-mind 0.10--recode--out$pre/lc.ci.final
plink--file$pre/lc.ci.final--allow-extra-chr-maf 0.05--recode--out$pre/lc.ci.final
##利用beagle4软件进行缺失基因型填补
pre="path"
plink--file$pre/lc.ci.final--allow-extra-chr--make-bed--out$pre/lc.ci.final plink--bfile$pre/lc.ci.final--allow-extra-chr--recode vcf-iid--out$pre/lc.ci.final
java–jarbeagle/beagle4.jar gt=$pre/lc.ci.final.vcf out=$pre/lc.ci.input.final nthreads=30
gunzip$pre/lc.ci.input.final.vcf.gz
vcftools--vcf$pre/lc.ci.input.final.vcf--plink--out$pre/lc.ci.input.final
##利用Haploview软件找Tag SNP
pre="path"
mkdir$pre/LG
for i in$(seq 1 25)
do
plink--file$pre/lc.ci.input.final--allow-extra-chr--chr LG${i}--recode--out$pre/LG/LG${i}
awk'{print$2"\t"$4}'$pre/LG/LG${i}.map>$pre/LG/LG${i}.map1
mv$pre/LG/LG${i}.map1$pre/LG/LG${i}.map
sed-i's/-9/0/g'$pre/LG/LG${i}.ped
java-jar haploview/Haploview.jar-nogui-pedfile$pre/LG/LG${i}.ped-info$pre/LG/LG${i}.map-pairwiseTagging-log$pre/LG/LG${i}.log-out$pre/LG/LG${i}
cat$pre/LG/LG${i}.TESTS>>$pre/tag.SNP
done
plink--file$pre/lc.ci.input.final--allow-extra-chr--extract$pre/tag.SNP--recode--out$pre/lc.ci.input.tag.final
plink--file$pre/lc.ci.input.tag.final--allow-extra-chr--maf 0.05--recode--out$pre/lc.ci.input.tag.final
⑸建立最佳全基因组选择模型;
①基于参考群体的二分类抗性表型和基因型进行全基因组关联分析,所使用软件为R语言环境下的gaston包,得到所有SNP位点的显著性P值,并将其由小到大排序。
②使用R语言环境下BGLR包中的Bayesian LASSO分析模型建立全基因组选择模型,以10个标记为增量,按照已经排好序的SNP标记,逐渐增加排在前面的标记使用数量,并利用5折交叉验证计算每个标记量下的模型预测准确性。5折交叉验证将参考群体随机分为5份,其中4份为训练集,1份为测试集,即利用训练集建立的模型来预测测试集,预测准确性为测试集的GEBV与真实的观察值之间的相关系数的平方。最终发现,在600个标记下,预测准确性就达到最高。
③判断模型预测准确性最高时,对应的标记数量,相应的此标记密度下的建立的模型即是最佳的全基因组选择模型。
这部分分析工作在linux系统下运行,并配合R语言脚本。以下是具体的执行脚本。
##shell脚本,生成R语言脚本所必须的输入文件
pre="path";file="lc.ci.input.tag.final"
awk'{if(/^[34]/){$6=2;print}else if(/^S/){$6=1;print}else{$6="NA";print}}'$pre/${file}.ped>$pre/test.ped
mv$pre/test.ped$pre/${file}.ped
awk'{if(/HX/){}else{print$1"\t"$2}}'$pre/${file}.ped>$pre/reference.id
plink--file$pre/$file--allow-extra-chr--keep$pre/reference.id--recode--out$pre/${file}.ref#生成只有参考群体的plink格式
plink--file$pre/${file}.ref--allow-extra-chr--maf 0.05--recode--out$pre/${file}.ref
cut-f 2$pre/${file}.ref.map>$pre/test.id
plink--file$pre/$file--allow-extra-chr--extract$pre/test.id--recode--out$pre/$file
rm$pre/test.id
plink--file$pre/$file--allow-extra-chr--make-bed--recode--out$pre/$file
plink--file$pre/$file--allow-extra-chr--freq--out$pre/$file
awk'{printf$1"\t";for(i=7;i<=NF;i=i+2){printf"%s%s\t",$i,$(i+1)};print""}'$pre/${file}.ped|sed's//\//g'|awk'{if(NR==FNR){a[FNR-1]=$3}else{printf$1;for(i=2;i<=NF;i++){split($i,b,"/");printf""(a[i-1]==b[1])+(a[i-1]==b[2])};print""}}'$pre/${file}.frq->$pre/${file}.geno
printf"\nstart generating pheno-data\n"#产生表型文件
awk'{print$1"\t"$6}'$pre/${file}.ped>$pre/${file}.pheno
Rscript$pre/script/binary.gs.best.marker.R$pre/$file
##R语言脚本,binary.gs.best.marker.R。主要分析任务是判断最佳标记数量。
args<-commandArgs(T)
lc.geno<-read.table(paste(args[1],".geno",sep=""),sep="")
lc.map<-read.table(paste(args[1],".map",sep=""),sep="\t")
lc.pheno<-read.table(paste(args[1],".pheno",sep=""),sep="\t")
colnames(lc.pheno)<-c("Taxa","staus")
lc.map<-lc.map[c(2,1,4)]
colnames(lc.map)<-c("SNP","Chromosome","Position")
colnames(lc.geno)<-c("taxa",as.character(lc.map$SNP))
rownames(lc.geno)<-lc.geno$taxa
lc.geno<-lc.geno[-1]
rownames(lc.pheno)<-lc.pheno$Taxa
all.id<-rownames(lc.pheno)
ref.id<-all.id[-grep("HX",all.id)]
sur.id<-ref.id[grep("S",ref.id)];death.id<-ref.id[-grep("S",ref.id)]
lc.pheno[sur.id,2]<-1;lc.pheno[death.id,2]<-0
head(lc.pheno)
lc.geno[1:5,1:5]
library(BGLR);library(sommer);library(gaston)
lc.bim<-read.table(paste(args[1],".bim",sep=""),sep="\t")
lc.fam<-read.table(paste(args[1],".fam",sep=""),sep="")
names(lc.bim)<-c("chr","id","dist","pos","A1","A2")
names(lc.fam)<-c("famid","id","father","mother","sex","pheno")
rownames(lc.fam)<-lc.fam$id
lc.fam[sur.id,6]<-1;lc.fam[death.id,6]<-0
lc.bim$chr<-gsub("LG","",lc.bim$chr)
options(gaston.auto.set.stats=
TRUE,gaston.autosomes=c(1:25),gaston.chr.x=FALSE,gaston.chr.y=FALSE)
lc.geno<-as.matrix(lc.geno)
lc.geno1<-lc.geno[ref.id,]
lc.fam1<-lc.fam[ref.id,]
lc.bed.matrix<-as.bed.matrix(lc.geno1,lc.fam1,lc.bim)
gaston.result<-association.test(lc.bed.matrix,K=GRM(lc.bed.matrix),method="lmm",response="binary")
head(gaston.result)
plink.result<-read.table(paste(args[1],".ref.assoc",sep=""),sep="",header=TRUE)cor(gaston.result$p,plink.result$P)
write.table(gaston.result,paste(args[1],".gaston.result",sep=""),quote=F,row.names=F)
gaston.result<-gaston.result[order(gaston.result$p),]
cat("\n--------OK-------\n")
library(foreach)
library(doParallel)
cl<-makeCluster(40)
registerDoParallel(cores=40)
lc.pheno1<-lc.pheno[ref.id,]
y<-10
x<-floor(nrow(lc.map)/y)
z<-floor(nrow(lc.pheno1)/5)
m<-40
cor<-list()
fit<-list()
for(j in 1:x){
geno<-lc.geno1[,as.character(gaston.result$id[1:(y*j)])]
foreach(i=1:m)%dopar%sample(1:nrow(lc.geno1),z)->tst
foreach(i=1:m)%dopar%lc.pheno1->trait
for(n in 1:m){trait[[n]][tst[[n]],2]<-NA}
foreach(i=1:m)%dopar%BGLR(y=trait[[i]][,2],response_type='ordinal',ETA=list(list(X=geno,model='BL')),nIter=1200,burnIn=200,verbose=F)->ans
foreach(i=1:m,.combine="c")%dopar%cor(ans[[i]]$probs[-tst[[i]],2],lc.pheno1[-tst[[i]],2])->fittness
foreach(i=1:m,.combine="c")%dopar%cor(ans[[i]]$probs[tst[[i]],2],lc.pheno1[tst[[i]],2])->predict
fit[[j]]<-fittness;cor[[j]]<-predict
cat(paste("\nthis is",j,"time\n\n"))
}
stopCluster(cl)
predict<-data.frame(marker=sapply(cor,"[",1))
for(i in 1:m){predict[i]<-sapply(cor,"[",i)}
names(predict)<-paste("rep",1:m,sep="")
predict<-data.frame(number=(1:x)*y,predict)
fitness<-data.frame(marker=sapply(fit,"[",1))
for(i in 1:m){fitness[i]<-sapply(fit,"[",i)}
names(fitness)<-paste("rep",1:m,sep="")
fitness<-data.frame(number=(1:x)*y,fitness)
write.table(predict,paste(args[1],".predict",sep=""),sep="\t",quote=F,row.names=F)
write.table(fitness,paste(args[1],".fitness",sep=""),sep="\t",quote=F,row.names=F)
⑹对选育群体的基因组育种值(GEBV)进行估算;
①在最佳全基因组选择模型下,对选育群体的GEBV进行估计,并对其进行排名。
②利用人类中应用十分广泛的多基因风险评分来佐证上述得到的GEBV值的准确性,所使用的软件为PRSice软件。
这部分工作由以下执行脚本完成。
##R语言脚本,binary.gs.GEBV.R。最佳标记数量下,建立全基因组模型,并计算选育群体GEBV。
args<-commandArgs(T)
lc.geno<-read.table(paste(args[1],".geno",sep=""),sep="")
lc.map<-read.table(paste(args[1],".map",sep=""),sep="\t")
lc.pheno<-read.table(paste(args[1],".pheno",sep=""),sep="\t")
colnames(lc.pheno)<-c("Taxa","staus")
lc.map<-lc.map[c(2,1,4)]
colnames(lc.map)<-c("SNP","Chromosome","Position")
colnames(lc.geno)<-c("taxa",as.character(lc.map$SNP))
rownames(lc.geno)<-lc.geno$taxa
lc.geno<-lc.geno[-1]
rownames(lc.pheno)<-lc.pheno$Taxa
all.id<-rownames(lc.pheno)
slect.id<-all.id[grep("HX",all.id)]
ref.id<-all.id[-grep("HX",all.id)]
sur.id<-ref.id[grep("S",ref.id)];death.id<-ref.id[-grep("S",ref.id)]
lc.pheno[sur.id,2]<-1;lc.pheno[death.id,2]<-0;lc.pheno[slect.id,2]<-NA
print(lc.pheno)
cat("\n--------OK-------\n")
predict<-read.table(paste(args[1],".predict",sep=""),row.names="number",header=T,sep="\t")
predict.mean<-apply(predict,1,mean)
predict.sd<-apply(predict,1,sd)
best.number<-which.max(predict.mean)*10
print(best.number)
cat("\n--------OK-------\n")
gaston.result<-read.table(paste(args[1],".gaston.result",sep=""),header=T)
head(gaston.result)
gaston.result<-gaston.result[order(gaston.result$p),]
head(gaston.result)
lc.geno1<-lc.geno[ref.id,]
lc.pheno1<-lc.pheno[ref.id,]
z<-floor(nrow(lc.pheno1)/5)
library(BGLR)
lc.pheno1->trait
tst<-sample(1:nrow(lc.geno1),z);trait[tst,2]<-NA
geno<-lc.geno1[,as.character(gaston.result$id[1:best.number])]
BGLR(y=trait[,2],response_type='ordinal',
ETA=list(list(X=geno,model='BL')),nIter=1200,burnIn=200,verbose=F)->ans
class(ans$probs);dim(ans$probs)
cv_5<-data.frame(lc.pheno1,ans$probs,stat="train",stringsAsFactors=F)
print(cv_5)
cat("\n--------OK-------\n")
cv_5[tst,"stat"]<-"tst";names(cv_5)[c(3,4)]<-c("death_prob","survival_prob")write.table(cv_5,paste(args[1],".cv_5",sep=""),sep="\t",quote=F,row.names=F)cor(ans$probs[-tst,2],lc.pheno1[-tst,2]);cor(ans$probs[tst,2],lc.pheno1[tst,2])cat("\n--------OK-------\n")
geno<-lc.geno[,as.character(gaston.result$id[1:best.number])]
BGLR(y=lc.pheno[,2],response_type='ordinal',
ETA=list(list(X=geno,model='BL')),nIter=1200,burnIn=200,verbose=F)->ans
GEBV<-data.frame(lc.pheno,ans$probs,stat="train",stringsAsFactors=F)
head(rownames(GEBV));GEBV[slect.id,"stat"]<-"slection"
names(GEBV)[c(3,4)]<-c("death_prob","survival_prob")
write.table(GEBV,paste(args[1],".GEBV",sep=""),sep="\t",quote=F,row.names=F)cat("\n--------OK-------\n")
##shell脚本执行多基因风险评分分析任务
sort-k 8,8g$pre/${file}.gaston.result|cut-f 3-d""|sed 1d|head-670>$pre/effect.snp.id
plink--file$pre/$file--allow-extra-chr--extract$pre/effect.snp.id--make-bed--recode--out$pre/${file}.target
grep-f$pre/effect.snp.id-F-w$pre/${file}.ref.assoc|cut-f 1,2,3,4,7,9,10-d"">$pre/${file}.train.BASE
cut-f 1,2,3,4,7,9,10-d""$pre/${file}.ref.assoc|head-1|cat-$pre/${file}.train.BASE>$pre/txt
mv$pre/txt$pre/${file}.train.BASE
awk'{if(NR>1){$4=1;$5=2;print}else{print}}'$pre/${file}.train.BASE>$pre/txt
mv$pre/txt$pre/${file}.train.BASE
sed-i's/^LG//g'$pre/${file}.train.BASE
awk-v OFS="\t"'{$5=1;$6=2;print}'$pre/${file}.target.bim>$pre/txt
mv$pre/txt$pre/${file}.target.bim
sed-i's/^LG//g'$pre/${file}.target.bim
awk'{$5=1;print}'$pre/${file}.target.fam>$pre/txt
mv$pre/txt$pre/${file}.target.fam
sh$pre/script/binary.gs.PRS.sh$pre/$file
cut-f 4-d""$pre/${file}.prs.best|paste$pre/${file}.GEBV-|cut-f1,2,4,5,6>$pre/${file}.gs.result
awk'{if(NR>1&&$4=="slection"){print}else if(NR==1){print}}'
$pre/${file}.gs.result>$pre/${file}.gs.selection.result
pre=/public/home/fgl/public/softwares/PRSice/
Rscript$pre/PRSice.R--dir$pre\
--prsice$pre/PRSice_linux\
--base${1}.train.BASE\
--target${1}.target\
--thread 1\
--no-clump\
--upper 1\
--scatter-r2\
--stat OR\
--out${1}.prs\
--binary-target T
⑺以一定的选择强度选取选育群体GEBV值排在前面的个体作为亲鱼,并繁殖产生抗刺激隐核虫子一代;
①2019年4月中旬,从360尾有基因型的选育群体中根据电子标签挑选GEBV排名前47的个体作为抗性子一代的亲鱼,如此,选择强度为47/360=13%。挑选完剩余的亲鱼繁殖的后代作为对照组。
②按照正常的商品鱼养殖体系培育抗性子一代和对照组。在苗种培育流程中,抗性子一代已经表现出耐操作刺激等优良表现。
⑻对抗性子一代进行刺激隐核虫攻毒验证。
①在刺激隐核虫爆发的季节,2019年6月中旬,选取抗虫子一代和对照组2月龄幼鱼各150尾进行刺激隐核虫攻毒实验,实验流程同参考群体的抗性性状测评的攻毒流程一致。
②对攻毒结果进行生存分析发现,抗性子一代与对照组存在显著的抗性差异,表现为在96h时,抗虫子一代存活率为59.2%,而对照组为9.9%;在120h,抗虫子一代存活率为40.8%,而对照组已全部死亡;在160h,抗性子一代仍存活34.4%,且已从患病状态中恢复。
以上实施情况,表明本发明提供的基于全基因组选择的大黄鱼抗刺激隐核虫病优良品系的选育方法的可实用性和高效性。
图1为本发明在具体实施过程中的示意图。图2为本发明在具体实施过程中第⑵个步骤中参考群体抗性性状测量中死亡个体的死亡时间分布图。图3为本发明在具体实施过程中第⑸个步骤中判断全基因组模型最佳标记的结果图。结果表明,前600个标记(标记经过P值显著性排名)下,模型预测准确性达到最高。图4为本发明在具体实施过程中第⑹个步骤中使用PRSice软件对R语言BGLR包得到的GEBV结果的佐证结果,结果表明二者具有高度相关性,也佐证了对选育群体GEBV计算的准确性。图5为本发明在具体实施过程中第⑻个步骤中抗虫子一代的刺激隐核虫攻毒验证结果,结果表明,本发明培育出的抗性子一代抗虫性有显著的提升。
本发明将基因组育种技术首次应用到大黄鱼抗刺激隐核虫病的良种选育中。包括:建立参考群体;对参考群体群体进行抗性性状测量;建立选育群体;对参考群体和选育群体进行基因分型;建立最佳全基因组选择模型;对选育群体的基因组育种值GEBV进行估算;以一定的选择强度选取选育群体GEBV值排在前面的个体作为亲鱼,并繁殖产生抗刺激隐核虫子一代;对抗性子一代进行刺激隐核虫攻毒验证。利用该方法培育的抗刺激隐核虫子一代在攻毒验证中,96h存活率较对照组高出49.3%,抗虫性有显著的提升。提高对选育群体基因组育种值的估计准确性,大大缩短了育种年限,在一代之内就可获得高抗性的后代,减少刺激隐核虫对大黄鱼养殖产业造成的经济损失,同时也为其他鱼类的抗病育种提供借鉴和基础,具有广阔的应用前景。
Claims (4)
1.基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法,其特征在于包括以下步骤:
1)建立参考群体:
2)对参考群体进行抗性性状测量:
自攻毒起,开始不间断的跟踪每组鱼的攻毒状况;将失去平衡,且10s内不能恢复的个体视为死亡样本,记录死亡时间,并将其捞出实验桶外,剪其背部肌肉保存于纯酒精中;观察时间持续120 h,实验结束后将所有存活的个体转移至新的实验桶,进行攻毒后恢复;恢复两周后,仍然存活的个体视为存活样本;如此,大黄鱼抗刺激隐核虫的表型性状采集模式为二分类性状,存活-死亡;
3)建立选育群体:
4)对参考群体和选育群体进行基因分型:
5)建立最佳全基因组选择模型:
使用R语言环境下BGLR包中的Bayesian LASSO分析模型建立全基因组选择模型,以10个标记为增量,按照已经排好序的SNP标记,逐渐增加排在前面的标记使用数量,并利用5折交叉验证计算每个标记量下的模型预测准确性;
6)对选育群体的基因组育种值GEBV进行估算;
7)以一定的选择强度选取选育群体GEBV值排在前面的个体作为亲鱼,并繁殖产生抗刺激隐核虫子一代;
8)对抗性子一代进行刺激隐核虫攻毒验证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011171385.XA CN112273291B (zh) | 2020-10-28 | 2020-10-28 | 基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011171385.XA CN112273291B (zh) | 2020-10-28 | 2020-10-28 | 基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112273291A CN112273291A (zh) | 2021-01-29 |
CN112273291B true CN112273291B (zh) | 2021-09-07 |
Family
ID=74372487
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011171385.XA Active CN112273291B (zh) | 2020-10-28 | 2020-10-28 | 基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112273291B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113373245A (zh) * | 2021-07-14 | 2021-09-10 | 广东海洋大学 | 基于全基因组选择的马氏珠母贝金黄壳色性状良种培育方法 |
CN113678767B (zh) * | 2021-08-10 | 2022-08-23 | 中国水产科学研究院黄海水产研究所 | 一种对虾抗病性状的选育方法 |
CN114657262A (zh) * | 2022-03-31 | 2022-06-24 | 集美大学 | 一种抗内脏白点病大黄鱼的分子标记辅助选育方法 |
CN116064846A (zh) * | 2023-01-30 | 2023-05-05 | 中国海洋大学三亚海洋研究院 | 一种评估花鲈生长和抗性性状综合育种值的方法及应用 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106480189A (zh) * | 2016-10-18 | 2017-03-08 | 中国水产科学研究院黄海水产研究所 | 一种基于全基因组选择的鱼类抗病良种培育方法 |
CN107338321A (zh) * | 2017-08-29 | 2017-11-10 | 集美大学 | 一种确定最佳snp数量及其通过筛选标记对大黄鱼生产性能进行基因组选择育种的方法 |
CN110211635A (zh) * | 2019-06-12 | 2019-09-06 | 北京康普森农业科技有限公司 | 用于畜禽基因组选择分析的方法及畜禽育种方法 |
KR102168817B1 (ko) * | 2019-07-17 | 2020-10-22 | 대한민국 | 마이크로새틀라이트 마커를 이용한 향어 선발육종 방법 |
-
2020
- 2020-10-28 CN CN202011171385.XA patent/CN112273291B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106480189A (zh) * | 2016-10-18 | 2017-03-08 | 中国水产科学研究院黄海水产研究所 | 一种基于全基因组选择的鱼类抗病良种培育方法 |
CN107338321A (zh) * | 2017-08-29 | 2017-11-10 | 集美大学 | 一种确定最佳snp数量及其通过筛选标记对大黄鱼生产性能进行基因组选择育种的方法 |
CN110211635A (zh) * | 2019-06-12 | 2019-09-06 | 北京康普森农业科技有限公司 | 用于畜禽基因组选择分析的方法及畜禽育种方法 |
KR102168817B1 (ko) * | 2019-07-17 | 2020-10-22 | 대한민국 | 마이크로새틀라이트 마커를 이용한 향어 선발육종 방법 |
Non-Patent Citations (1)
Title |
---|
大黄鱼高密度连锁图谱构建及刺激隐核虫抗性和体型性状全基因组遗传解析;孔胜楠;《中国优秀硕士学位论文全文数据库 农业科技辑》;20200715(第7期);第9-11页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112273291A (zh) | 2021-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112273291B (zh) | 基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法 | |
Kahilainen et al. | Species introduction promotes hybridization and introgression in Coregonus: is there sign of selection against hybrids? | |
Seamons et al. | Sex-specific patterns of lifetime reproductive success in single and repeat breeding steelhead trout (Oncorhynchus mykiss) | |
Vähä et al. | Temporally stable population‐specific differences in run timing of one‐sea‐winter Atlantic salmon returning to a large river system | |
Limborg et al. | Genetic population structure of European sprat Sprattus sprattus: differentiation across a steep environmental gradient in a small pelagic fish | |
EP2201133A1 (en) | Methods of using genetic markers and related epistatic interactions | |
Dubé et al. | Parentage analyses identify local dispersal events and sibling aggregations in a natural population of Millepora hydrocorals, a free‐spawning marine invertebrate | |
Eldridge et al. | Genetic diversity over multiple generations of supplementation: an example from Chinook salmon using microsatellite and demographic data | |
Barbanti et al. | How many came home? Evaluating ex situ conservation of green turtles in the Cayman Islands | |
Beacham et al. | Variation in migration pattern, broodstock origin, and family productivity of coho salmon hatchery populations in British Columbia, Canada, derived from parentage‐based tagging | |
Ovenden et al. | Genetic population structure of mangrove jack, Lutjanus argentimaculatus (Forsskål) | |
Sunny et al. | Genetic diversity and demography of the critically endangered Roberts’ false brook salamander (Pseudoeurycea robertsi) in Central Mexico | |
Miller et al. | Reproductive dynamics of a native brook trout population following removal of non-native brown trout from a stream in Minnesota, north-central USA | |
Suneetha et al. | Genetic and morphological stock structure of the pearlside, Maurolicus muelleri (Pisces, Sternoptychidae), among Norwegian fjords and offshore area | |
Miller et al. | Introgression of hatchery rainbow trout in naturalized steelhead populations of western Lake Superior | |
Jue | Exploring the structure of genetic variation and the influences of demography on effective population size in the gag grouper Mycteroperca microlepi (Goode & Bean) | |
Smith | The Genomic Basis for Fitness and Ecomorphological Variation in Recovering Populations of Lake Trout (Salvelinus namaycush) in the Great Lakes | |
Sabolić | Genomic patterns of phenotypic adaptation in Italian wall lizard Podarcis siculus (Rafinesque-Schmaltz, 1810) | |
Semenova | Introgressive hybridization in the secondary contact area of the Atlantic herring Clupea harengus and the pacific herring C. pallasii (clupeidae): ecological basis, geographical structure, and temporal variability of the hybridization zone | |
Melland | The genetic mating system and population structure of the green-rumped parrotlet | |
Sanudi | Genetics studies of cichlid Oreochromis shiranus towards sustainable aquaculture and fisheries management in Malawi | |
Azuma et al. | Genetic population structure of Neptune whelk in northern Japan inferred from microsatellite DNA variation | |
Kim | The utility of DNA microsatellite markers in conservation of a Namibian population of the Black Rhinoceros (Diceros bicornis) | |
Thompson | Introgression between wild and selectively-bred Sydney rock oysters (Saccostrea glomerata) | |
Kuo et al. | Variation in the local grey mullet populations (Mugil cephalus) on the western Pacific fringe |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |