CN113678767A - 一种对虾抗病性状的选育方法 - Google Patents
一种对虾抗病性状的选育方法 Download PDFInfo
- Publication number
- CN113678767A CN113678767A CN202110914524.1A CN202110914524A CN113678767A CN 113678767 A CN113678767 A CN 113678767A CN 202110914524 A CN202110914524 A CN 202110914524A CN 113678767 A CN113678767 A CN 113678767A
- Authority
- CN
- China
- Prior art keywords
- population
- individual
- snp
- breeding
- reference population
- 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.)
- Granted
Links
- 238000009395 breeding Methods 0.000 title claims abstract description 126
- 241000238557 Decapoda Species 0.000 title claims abstract description 78
- 208000035240 Disease Resistance Diseases 0.000 title claims abstract description 35
- 230000001488 breeding effect Effects 0.000 claims abstract description 110
- 238000012360 testing method Methods 0.000 claims abstract description 36
- 230000004083 survival effect Effects 0.000 claims abstract description 31
- 208000015181 infectious disease Diseases 0.000 claims abstract description 18
- 230000001717 pathogenic effect Effects 0.000 claims abstract description 14
- 230000013011 mating Effects 0.000 claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims description 37
- 238000000034 method Methods 0.000 claims description 35
- 238000005516 engineering process Methods 0.000 claims description 15
- 238000003908 quality control method Methods 0.000 claims description 14
- 201000010099 disease Diseases 0.000 claims description 9
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 9
- 238000012268 genome sequencing Methods 0.000 claims description 6
- 108090000623 proteins and genes Proteins 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 4
- 238000009396 hybridization Methods 0.000 claims description 4
- 230000000717 retained effect Effects 0.000 claims description 2
- 244000052769 pathogen Species 0.000 description 9
- 230000000694 effects Effects 0.000 description 7
- 239000003550 marker Substances 0.000 description 7
- 210000003205 muscle Anatomy 0.000 description 6
- 238000003205 genotyping method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000009394 selective breeding Methods 0.000 description 5
- 230000000996 additive effect Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000007614 genetic variation Effects 0.000 description 4
- 241000238553 Litopenaeus vannamei Species 0.000 description 3
- 239000000654 additive Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 108700028369 Alleles Proteins 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002790 cross-validation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 238000004321 preservation Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000002791 soaking Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000007400 DNA extraction Methods 0.000 description 1
- 241000607272 Vibrio parahaemolyticus Species 0.000 description 1
- 241000700605 Viruses Species 0.000 description 1
- 238000000540 analysis of variance Methods 0.000 description 1
- 238000003975 animal breeding Methods 0.000 description 1
- 238000012098 association analyses Methods 0.000 description 1
- 230000001580 bacterial effect Effects 0.000 description 1
- 244000052616 bacterial pathogen Species 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000012248 genetic selection Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000008774 maternal effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009385 viral infection Effects 0.000 description 1
- 230000000007 visual effect Effects 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/50—Culture of aquatic animals of shellfish
- A01K61/59—Culture of aquatic animals of shellfish of crustaceans, e.g. lobsters or shrimps
-
- 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)
- Environmental Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Zoology (AREA)
- Animal Husbandry (AREA)
- Biodiversity & Conservation Biology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种对虾抗病性状的选育方法，所述方法包括：构建对虾家系的测试群体和保种群体；对测试群体进行病原感染，获取每个死亡个体的个体存活时间；基于测试群体构建参考群体；获得参考群体的SNP分型数据；估计参考群体的个体基因组育种值；获得保种群体的SNP分型数据；基于参考群体的所述个体存活时间、参考群体共有SNP位点的分型数据和保种群体共有SNP位点的分型数据，估计保种群体的个体基因组育种值；根据保种群体的个体基因组育种值选择保种个体进行留种和配种。应用本发明，可以提高对虾抗病性状的选择准确性。
Description
技术领域
本发明属于水产动物育种技术领域，具体地说，是涉及对虾育种技术，更具体地说，是涉及一种对虾抗病性状的选育方法。
背景技术
对虾是水产种业最具代表性的物种之一，对高质量苗种需求尤为迫切。以凡纳滨对虾为例，每年需亲虾150万对以上，苗种1.5万亿尾以上。以规模化家系为基础的选择育种技术是当前支撑对虾种业发展的核心技术。
对虾疾病一直困扰着其养殖业的发展，在对虾抗病良种选育过程中经常需要对养殖对虾进行人工感染实验，以研究对虾对疾病的抵抗能力。对于抗病性状来讲，如果使用感染实验中存活下来的对虾作为留种亲虾，有携带并传播病原的风险。因此，对于抗病性状的育种，通常采用保种群体与测试群体相分离的模式。具体做法是，核心群体中的每个全同胞家系个体分为两部分，一部分个体利用可视嵌入型荧光标记进行标记，混合后在特定的条件下进行性状测试；另外一部分留在单独的保种池，养至亲虾规格。依据性状测试群体的遗传评估结果，从保种群体中选择优秀的家系制定配种方案生产下一代家系。该模式虽然可以确保核心保种群体始终处于高等级的生物安保水平，但对于核心保种群体，选择仅在家系水平实施，只利用了一半的加性遗传变异，降低了育种目标性状的选择准确性和选择反应，严重影响了选育效率。
基因组选择育种最早由Meuwissen于2001年提出，其利用覆盖全基因组的SNP标记对参考群体和候选群体进行分型，根据参考群体的表型和基因型信息，获得每个SNP标记的效应，利用SNP效应对候选群体进行预测和评估。近年来，基因组选择育种方法在水产动物选育中进行了应用，例如，公开号为CN 110867208A的中国专利申请，公开了一种提高水产动物全基因组选择育种效率的方法，通过低密度标记实现准确的全基因组选择，充分利用了效应SNP标记，提高了全基因组选择准确率。但是，其仅是对育种群体进行全基因组关联分析，获得每个SNP位点的效应，筛选最优的用于全基因组选择的标记组合，作为效应SNP标记，利用效应SNP标记进行选择育种分析，达到降低分型成本的目的。在对虾抗病选育技术领域，还未见有基于基因组选择的方法。因此，如何将基因组选择育种技术应用到对虾抗病性状选育中，以提高选育效率，是亟待解决的技术问题。
发明内容
本发明的目的在于提供一种对虾抗病性状的选育方法，利用基因组选择方法估计保种群体的个体基因组育种值，利用个体基因组育种值进行个体选择，提高对虾抗病性状的选育准确性。
为实现上述发明目的，本发明采用下述技术方案予以实现：
一种对虾抗病性状的选育方法，其特征在于，所述方法包括：
对已构建的对虾家系，从每个家系中选择多个个体，分别构成测试群体和保种群体；
对所述测试群体进行病原感染，获取每个死亡个体的个体存活时间；
从所述测试群体内每个家系的死亡个体中选取多个个体，构建参考群体；
获取所述参考群体中每个个体的DNA，采用高通量SNP分型技术获得参考群体的SNP分型数据；
基于所述参考群体的所述个体存活时间和所述参考群体的SNP分型数据，估计参考群体的个体基因组育种值；
根据所述参考群体的所述个体基因组育种值计算家系平均育种值，按照所述家系平均育种值从高到低对所有家系进行排名；
从所有家系中选择排名在前的指定数量的家系，作为分析家系，从所述保种群体的分析家系中选择多个保种个体，采用高通量SNP分型技术获得保种群体的SNP分型数据；
将所述参考群体的SNP分型数据和所述保种群体的SNP分型数据进行合并分析，获得共有SNP位点，确定参考群体共有SNP位点的分型数据和保种群体共有SNP位点的分型数据；
基于所述参考群体的所述个体存活时间、所述参考群体共有SNP位点的分型数据和所述保种群体共有SNP位点的分型数据，估计保种群体的个体基因组育种值；
根据所述保种群体的个体基因组育种值选择保种个体进行留种和配种。
在其中一个优选实施例中，基于所述参考群体的所述个体存活时间和所述参考群体的SNP分型数据，估计参考群体的个体基因组育种值，具体包括：
基于所述参考群体的所述个体存活时间和所述参考群体的SNP分型数据，利用SSBR-A模型估计参考群体的个体基因组育种值；
基于所述参考群体的个体存活时间、所述参考群体共有SNP位点的分型数据和所述保种群体共有SNP位点的分型数据，估计保种群体的个体基因组育种值，具体包括：
基于所述参考群体的个体存活时间、所述参考群体共有SNP位点的分型数据和所述保种群体共有SNP位点的分型数据，利用SSBR-A模型估计保种群体的个体基因组育种值。
在其中一个优选实施例中，获取所述参考群体中每个个体的DNA，采用高通量SNP分型获得参考群体的SNP分型数据，具体包括：
获取所述参考群体中每个个体的DNA，利用简化基因组测序方法或基因芯片的方法进行高通量SNP分型，并对获得的SNP进行质量控制，质量控制后保留的SNP分型数据确定为参考群体的SNP分型数据。
在其中一个优选实施例中，从所述保种群体的分析家系中选择多个保种个体，采用高通量SNP分型技术获得保种群体的SNP分型数据，具体包括：
获取所述分析家系中选择的每个保种个体的DNA，利用简化基因组测序方法或基因芯片的方法进行高通量SNP分型，并对获得的SNP进行质量控制，质量控制后保留的SNP分型数据确定为保种群体的SNP分型数据。
在其中一个优选实施例中，获取所述分析家系中选择的每个保种个体的DNA，具体包括：
获取所述分析家系中选择的每个保种个体的触须，从触须中提取保种个体的DNA。
在其中一个优选实施例中，从所有家系中选择排名在前的指定数量的家系，作为分析家系，具体包括：
将所有家系中排名前N%的家系确定为所述分析家系；所述N满足：10＜N＜50。
在其中一个优选实施例中，所述对虾家系采用下述方法构建：
获取对虾种质群体，通过双列杂交方式构建育种基础群体，从所述基础群体中选取亲虾，采用巢式交配设计构建所述对虾家系。
与现有技术相比，本发明的优点和积极效果是：本发明提供的对虾抗病性状的选育方法中，构建测试群体和保种群体，测试群体用于病原感染，再基于参与病原感染的测试群体建立参考群体，基于参考群体中的个体存活时间、参考群体的SNP分型数据及保种群体的SNP分型数据估计保种群体的个体基因组育种值，最后根据保种群体的个体基因组育种值选择保种个体进行留种、配种，从而，将对虾抗病性状的选育从现有技术中的家系水平提高到保种群体中的个体水平，充分利用了对虾加性遗传变异，提高了育种抗病性状的选择准确性和选择反应；而且，通过选择测试群体中的部分个体构建参考群体，仅获取参考群体的SNP分型数据，再基于参考群体的SNP分型数据估计参考群体的个体基因组育种值，进而根据个体基因组育种值获取家系平均育种值，根据家系平均育种值选择排名靠前的家系作为分析家系，基于分析家系获得保种群体的SNP分型数据，从而能够减少测试群体及保种群体的分型数量，降低了分型成本，同时提升选育效率。
结合附图阅读本发明的具体实施方式后，本发明的其他特点和优点将变得更加清楚。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明对虾抗病性状的选育方法一个实施例的流程图；
图2为采用本发明对虾抗病性状的选育方法一个具体实施例中保种群体的个体基因组育种值分布示意图；
图3为采用不同模型估计个体基因组育种值预测准确性的结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下将结合附图和实施例，对本发明作进一步详细说明。
请参见图1，该图所示为本发明对虾抗病性状的选育方法一个实施例的流程图。
如图1所示，该实施例采用下述过程实现对虾抗病性状的选育。
步骤101：对已构建的对虾家系，从每个家系中选择多个个体，分别构成测试群体和保种群体。
对虾家系的构建，可以采用现有技术中构建家系的方法。
在一些优选实施例中，采用下述方法构建对虾家系：
获取对虾种质群体，通过双列杂交的方式构建育种基础群体，从基础群体中选取亲虾，采用巢式交配设计构建对虾家系。通过双列杂交方式构建育种基础群体，能够包含丰富的遗传变异；通过巢式交配设计构建家系，有助于在遗传评估时从加性效应中剖分出共同环境效应，使育种值估计更加准确。
在家系中的对虾个体长至一定长度后，从每个家系中选择部分个体，构成测试群体，用于抗病测试；选择另一部分个体，构成保种群体。一般的，每个家系中构成测试群体的个体数量小于构成保种群体的个体数量。对于保种群体，在生物安全环境下进行养殖。
步骤102：对测试群体进行病原感染，获取每个死亡个体的个体存活时间。
对于测试群体，通过可视性荧光标记组合进行不同家系的区分，然后进行病原感染。病原感染时，针对病毒病原，采用单尾口饲的方式进行人工病毒感染，具体方法可采用现有技术来实现；针对细菌病原，以家系为单位采用细菌终浓度为107cfu/mL的菌液浸泡2小时以上进行对虾人工病原感染。优选的，都选择一个较高的病原浓度进行感染，保证90%的个体在一定时间（例如，小于30天）内能够全部死亡。
感染结束后，参与抗病测试的对虾全部转入同一养殖池中进行混合养殖，正常投饵、换水，每隔2小时观察对虾死亡情况，并收集死亡个体，取其肌肉组织置于-80℃保存，以便于后续使用。从第一尾对虾死亡开始记录个体死亡时间和所属家系，并给死亡个体编号。死亡高峰时，每隔1小时收集一次死虾；待90%以上的对虾死亡时，结束测试。死亡定义为：对虾失去平衡，倒向一侧，并对外界刺激失去反应。从感染结束开始计算，根据个体死亡时间，确定每个死亡个体的个体存活时间。
步骤103：从测试群体内每个家系的死亡个体中选取多个个体，构建参考群体。
为便于后续基因分型，从测试群体内每个家系的死亡个体中选取多个个体，所选取的个体要求具有保存较好的肌肉组织，以便于提取DNA。每个家系选取的多个死亡个体作为样本，构建参考群体。
步骤104：获取参考群体中每个个体的DNA，采用高通量SNP分型技术获得参考群体的SNP分型数据。
从参考群体中每个个体所保存的肌肉组织中提取其DNA，然后，采用高通量SNP分型技术获得参考群体的SNP分型数据。提取DNA、获得SNP分型数据的过程，可采用现有技术实现。
在其他一些优选实施例中，在获取参考群体中每个个体的DNA之后，利用简化基因组测序方法或基因芯片的方法进行高通量SNP分型，并对获得的SNP进行质量控制，质量控制后保留的SNP分型数据确定为参考群体的SNP分型数据。通过对SNP进行质量控制，能够剔除数据质量较差的SNP和个体，提高后续分析的准确性。
步骤105：基于参考群体的个体存活时间和参考群体的SNP分型数据，估计参考群体的个体基因组育种值。
基于个体存活时间和SNP分型数据估计个体基因组育种值，可以采用现有技术中的多种模型来实现。例如，采用GBLUP（genomic best linear unbiased prediction，基因组最佳线性无偏预测）模型、ssGBLUP（single-step genomic best linear unbiasedprediction，一步法基因组最佳线性无偏预测）模型、Bayes A（贝叶斯A）模型等。
针对对虾抗病性状的选育，经过分析和验证，优选采用SSBR-A（single-stepBayesian regression A，一步法贝叶斯回归A）模型估计参考群体的个体基因组育种值，能够获得准确性更高的个体基因组育种值。
步骤106：根据参考群体的个体基因组育种值计算家系平均育种值，按照家系平均育种值从高到低对所有家系进行排名。
根据个体基因组育种值计算家系平均育种值的方法，采用现有技术来实现。
步骤107：从所有家系中选择排名在前的指定数量的家系作为分析家系，从保种群体的分析家系中选择多个保种个体，采用高通量SNP分型技术获得保种群体的SNP分型数据。
为降低后续基因分型成本以及提高选育效率，从所有家系中选择排名在前的指定数量的家系作为分析家系，从保种群体的分析家系中选择多个保种个体，作为进行基因分型的保种个体，采用高通量SNP分型技术，获得保种群体的SNP分型数据。
在一些优选实施例中，从所有家系中选择排名在前的指定数量的家系，作为分析家系，具体包括：将所有家系中排名前N%的家系确定为分析家系。也即，将排名为前N%的所有家系的数量确定为指定数量，相应的家系作为分析家系。其中，N满足：10＜N＜50，既能减少分型的样本数量，节省分型成本，还能提高选择强度以及选育效率。
在获得保种群体的SNP分型数据时，首先需要获得分析家系中的保种个体的DNA。优选的，获取每个保种个体的触须提取其DNA。例如，获取保种个体断落的触须，然后提取DNA。然后，采用高通量SNP分型技术获得保种群体的SNP分型数据。提取DNA、获得SNP分型数据的具体过程，可采用现有技术实现。
在其他一些优选实施例中，在获取保种群体中每个个体的DNA之后，利用简化基因组测序方法或基因芯片的方法进行高通量SNP分型，并对获得的SNP进行质量控制，质量控制后保留的SNP分型数据确定为保种群体的SNP分型数据。通过对SNP进行质量控制，能够剔除数据质量较差的SNP和个体，提高后续分析的准确性。
步骤108：将参考群体的SNP分型数据和保种群体的SNP分型数据进行合并分析，获得共有SNP位点，确定参考群体共有SNP位点的分型数据和保种群体共有SNP位点的分型数据。
步骤109：基于参考群体的个体存活时间、参考群体共有SNP位点的分型数据和保种群体共有SNP位点的分型数据，估计保种群体的个体基因组育种值。
同样的，基于个体存活时间和SNP分型数据估计个体基因组育种值，可以采用现有技术中的多种模型来实现。例如，采用GBLUP模型、ssGBLUP模型、Bayes A模型等。
而针对对虾抗病性状的选育，经过分析和验证，优选采用SSBR-A模型估计保种群体的个体基因组育种值，能够获得准确性更高的个体基因组育种值。
步骤110：根据保种群体的个体基因组育种值选择保种个体进行留种和配种。
具体的，选择个体基因组育种值排名靠前的保种个体进行留种和配种，实现抗病性能优良的保种个体的选育。
采用上述实施例的对虾抗病性状的选育方法中，构建测试群体和保种群体，实现测试群体与保种群体的分离。测试群体用于病原感染，再基于参与病原感染的测试群体建立参考群体，基于参考群体中的个体存活时间、参考群体的SNP分型数据及保种群体的SNP分型数据估计保种群体的个体基因组育种值，最后根据保种群体的个体基因组育种值选择保种个体进行留种、配种，从而，将对虾抗病性状的选育从现有技术中的家系水平提高到保种群体中的个体水平，充分利用了对虾加性遗传变异，提高了抗病性状育种的选择准确性和选择反应。而且，通过选择测试群体中的部分个体构建参考群体，仅获取参考群体的SNP分型数据，再基于参考群体的SNP分型数据估计参考群体的个体基因组育种值，进而根据个体基因组育种值获取家系平均育种值，根据家系平均育种值选择排名靠前的家系作为分析家系，基于分析家系获得保种群体的SNP分型数据，从而能够减少测试群体及保种群体的分型样本数量，同时提升选择强度以及选育效率。
下面以一个具体实例进一步描述基于本发明提供的方法进行对虾抗病性状的选育的过程。
2020年5月份，在某养殖基地构建凡纳滨对虾家系40个，其中包含6个母系半同胞家系。于7月份从每个家系中选取20尾4～5cm的个体，作为抗病测试群体，通过在个体尾部肌肉中注射不同的可视性荧光标记组合以区分不同家系，随后进行副溶血弧菌（Vibrioparahaemolyticus, VPAHPND）感染。另外再从每个家系随机选取100尾个体作为保种群体在遗传育种中心小水体生物安全环境下进行独立保种。
感染时，以家系为单位放入单独的塑料容器中，采用终浓度为107cfu/mL的VPAHPND菌液浸泡3小时。感染结束后，将所有待测试个体全部转入同一大型混养池中，水温调整26-28℃，每日按照所有对虾体重的5％分4次投喂人工配合饲料。每隔2小时观察对虾死亡情况，收集死亡个体，取肌肉组织置于-80℃保存。从第一尾对虾死亡开始记录个体死亡时间和所属家系，并给死亡个体编号，将存活时间记作个体抗病性表型；死亡高峰时，每个小时收集一次死虾；待90%以上的对虾死亡时，结束测试。
抗病测试于2020年7月25日18:00开始，至2020年9月8日16:00结束，剩余21尾存活个体，记录有存活时间表型数据的个体数量为762尾。第一尾死亡个体与最后一尾死亡个相差1074个小时，家系平均存活时间从138.3小时至428.5小时。方差分析显示，家系间的平均存活时间差异显著（P < 0.001），表明具有较大的选择潜力和较好的选择效果。
从有存活时间表型记录的个体中挑选了447个肌肉组织保存较好的样本构建参考群体。采用SNP芯片对447尾个体进行SNP分型，获得45815个SNP位点。使用PLINK1.9软件进行质控，剔除最小等位基因频率小于0.05与基因型检出率小于0.9的位点，剔除个体检出率小于0.8的个体，最终保留444尾个体和35927个位点。
基于参考群体的个体存活时间和SNP分型数据，利用SSBR-A模型估计个体抗病表型的基因组育种值。
使用SSBR-A模型估计的参考群体的个体基因组育种值计算家系平均育种值，按家系平均育种值从高到低选择前10个家系（排名前25%），作为分析家系。从每个分析家系的保种个体中挑选健康的雌虾和雄虾各5尾，将每个保种个体的其中一个眼柄套上带有唯一编号的环状眼标作为个体标识，并取其一根触须提取基因组DNA。
共计取样100尾个体（50尾雌虾与50尾雄虾），利用SNP芯片对100个样本进行高通量SNP分型，获得42300个SNP位点。使用PLINK1.9软件进行质控，剔除最小等位基因频率小于0.05与基因型检出率小于0.9的位点，剔除个体检出率小于0.9的个体，最终保留97尾个体和31212个位点。
将参考群体的SNP分型数据和保种群体的SNP分型数据进行合并分析，获得共有SNP位点，确定参考群体共有SNP位点的分型数据和保种群体共有SNP位点的分型数据。然后，使用SSBR-A模型估计100尾保种个体的个体基因组育种值。保种群体中100尾保种个体的个体基因组育种值分布示意图如图2所示。图2中，横坐标为家系序号，1-10共10个家系；纵坐标为个体基因组育种值。从图2可以看出，家系内和家系间的个体基因组育种值都存在明显差异。最后，从保种群体中选择个体基因组育种值最高的25尾雌虾和25尾雄虾用于留种和后续配种，从而可以进一步提高遗传选择进展。
图3所示为采用不同模型估计个体基因组育种值预测准确性的结果示意图。
具体的，基于上述具体实例中获取的参考群体的个体存活时间和SNP分型数据，分别利用GBLUP模型、ssGBLUP 模型、Bayes A模型、Bayes B（贝叶斯B）模型、Bayes C（贝叶斯C）模型、SSBR-A模型、SSBR-B（single-step Bayesian regression B，一步法贝叶斯回归B）模型和SSBR-C（single-step Bayesian regression C，一步法贝叶斯回归C）模型，估计个体抗病表型的基因组育种值。然后，采用五倍交叉验证的方法评估不同模型对基因组育种值的预测准确性。预测准确性定义为个体估计育种值与表型值间的皮尔森相关系数。交叉验证时，首先将参考群体所有个体随机均匀地分成5组，取其中1组数据作为候选群，表型设为缺失，其余4组数据合并在一起作为参考群用以计算候选群个体的育种值，然后估计育种值预测准确性。共计重复5次，使每组数据有且仅有一次机会作为候选群进行育种值估算。最后取5次的平均值作为每种模型的个体基因组育种值预测准确性。图3示出了采用各种模型估计的预测准确性的结果示意图，图中，横坐标为采用的不同模型，纵坐标为预测准确性。从图3可以看出，基于SSBR-A模型的预测准确性为0.294，明显高于其他模型的预测准确性（0.17-0.29），表明采用SSBR-A模型进行对虾抗病性状的选育，能够获得准确性更高的个体基因组育种值。
以上实施例仅用以说明本发明的技术方案，而非对其进行限制；尽管参照前述实施例对本发明进行了详细的说明，对于本领域的普通技术人员来说，依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或替换，并不使相应技术方案的本质脱离本发明所要求保护的技术方案的精神和范围。
Claims (7)
1.一种对虾抗病性状的选育方法，其特征在于，所述方法包括：
对已构建的对虾家系，从每个家系中选择多个个体，分别构成测试群体和保种群体；
对所述测试群体进行病原感染，获取每个死亡个体的个体存活时间；
从所述测试群体内每个家系的死亡个体中选取多个个体，构建参考群体；
获取所述参考群体中每个个体的DNA，采用高通量SNP分型技术获得参考群体的SNP分型数据；
基于所述参考群体的所述个体存活时间和所述参考群体的SNP分型数据，估计参考群体的个体基因组育种值；
根据所述参考群体的所述个体基因组育种值计算家系平均育种值，按照所述家系平均育种值从高到低对所有家系进行排名；
从所有家系中选择排名在前的指定数量的家系，作为分析家系，从所述保种群体的分析家系中选择多个保种个体，采用高通量SNP分型技术获得保种群体的SNP分型数据；
将所述参考群体的SNP分型数据和所述保种群体的SNP分型数据进行合并分析，获得共有SNP位点，确定参考群体共有SNP位点的分型数据和保种群体共有SNP位点的分型数据；
基于所述参考群体的所述个体存活时间、所述参考群体共有SNP位点的分型数据和所述保种群体共有SNP位点的分型数据，估计保种群体的个体基因组育种值；
根据所述保种群体的个体基因组育种值选择保种个体进行留种和配种。
2.根据权利要求1所述的对虾抗病性状的选育方法，其特征在于，基于所述参考群体的所述个体存活时间和所述参考群体的SNP分型数据，估计参考群体的个体基因组育种值，具体包括：
基于所述参考群体的所述个体存活时间和所述参考群体的SNP分型数据，利用SSBR-A模型估计参考群体的个体基因组育种值；
基于所述参考群体的个体存活时间、所述参考群体共有SNP位点的分型数据和所述保种群体共有SNP位点的分型数据，估计保种群体的个体基因组育种值，具体包括：
基于所述参考群体的个体存活时间、所述参考群体共有SNP位点的分型数据和所述保种群体共有SNP位点的分型数据，利用SSBR-A模型估计保种群体的个体基因组育种值。
3.根据权利要求1所述的对虾抗病性状的选育方法，其特征在于，获取所述参考群体中每个个体的DNA，采用高通量SNP分型技术获得参考群体的SNP分型数据，具体包括：
获取所述参考群体中每个个体的DNA，利用简化基因组测序方法或基因芯片的方法进行高通量SNP分型，并对获得的SNP进行质量控制，质量控制后保留的SNP分型数据确定为参考群体的SNP分型数据。
4.根据权利要求1所述的对虾抗病性状的选育方法，其特征在于，从所述保种群体的分析家系中选择多个保种个体，采用高通量SNP分型技术获得保种群体的SNP分型数据，具体包括：
获取所述分析家系中选择的每个保种个体的DNA，利用简化基因组测序方法或基因芯片的方法进行高通量SNP分型，并对获得的SNP进行质量控制，质量控制后保留的SNP分型数据确定为保种群体的SNP分型数据。
5.根据权利要求4所述的对虾抗病性状的选育方法，其特征在于，获取所述分析家系中选择的每个保种个体的DNA，具体包括：
获取所述分析家系中选择的每个保种个体的触须，从触须中提取保种个体的DNA。
6.根据权利要求1所述的对虾抗病性状的选育方法，其特征在于，从所有家系中选择排名在前的指定数量的家系，作为分析家系，具体包括：
将所有家系中排名前N%的家系确定为所述分析家系；所述N满足：10＜N＜50。
7.根据权利要求1至6中任一项所述的对虾抗病性状的选育方法，其特征在于，所述对虾家系采用下述方法构建：
获取对虾种质群体，通过双列杂交方式构建育种基础群体，从所述基础群体中选取亲虾，采用巢式交配设计构建所述对虾家系。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110914524.1A CN113678767B (zh) | 2021-08-10 | 2021-08-10 | 一种对虾抗病性状的选育方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110914524.1A CN113678767B (zh) | 2021-08-10 | 2021-08-10 | 一种对虾抗病性状的选育方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113678767A true CN113678767A (zh) | 2021-11-23 |
| CN113678767B CN113678767B (zh) | 2022-08-23 |
Family
ID=78579283
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110914524.1A Active CN113678767B (zh) | 2021-08-10 | 2021-08-10 | 一种对虾抗病性状的选育方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113678767B (zh) |
Citations (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060150273A1 (en) * | 2003-06-03 | 2006-07-06 | Tapper Brian A | Grass endophytes |
| AU2007200633A1 (en) * | 2000-07-27 | 2007-03-08 | The Johns Hopkins University | Growth differentiation factor receptors, agonists and antagonists thereof, and methods of using same |
| CN101011037A (zh) * | 2007-01-27 | 2007-08-08 | 中国水产科学研究院黄海水产研究所 | 鲈鱼冷冻精子诱导半滑舌鳎卵减数雌核发育方法 |
| CN101016519A (zh) * | 2006-09-29 | 2007-08-15 | 武汉大学 | 对虾抗病毒促生长双功能工程菌株及构建方法和应用 |
| CN101120660A (zh) * | 2007-08-28 | 2008-02-13 | 中国水产科学研究院黄海水产研究所 | 鱼类异源冷冻精子诱导雌核发育的方法 |
| CN101185426A (zh) * | 2007-11-20 | 2008-05-28 | 中国水产科学研究院黄海水产研究所 | 无特定病原中国对虾零换水苗种生产方法 |
| CN101215565A (zh) * | 2008-01-11 | 2008-07-09 | 华中农业大学 | 一种作为猪标记辅助选择的与免疫性状相关的分子标记的克隆及应用 |
| CN101292635A (zh) * | 2008-06-11 | 2008-10-29 | 大连水产学院 | 菲律宾蛤仔杂交系的制种方法 |
| CN101317554A (zh) * | 2008-06-11 | 2008-12-10 | 大连水产学院 | 以放射肋为标记的菲律宾蛤仔制种工艺 |
| CN101386857A (zh) * | 2008-10-21 | 2009-03-18 | 华中农业大学 | 一种与番茄抗病性状相关的分子标记克隆及应用 |
| CN101433188A (zh) * | 2008-12-10 | 2009-05-20 | 中国水产科学研究院黄海水产研究所 | 罗非鱼选择育种个体标记方法 |
| US20100037330A1 (en) * | 2008-06-11 | 2010-02-11 | Kannika Siripattarapravat | Efficient Somatic Cell Nuclear Transfer In Fish |
| WO2010068999A1 (en) * | 2008-12-19 | 2010-06-24 | Genetysis Pty Ltd | Method for identifying an animal suitable for breeding |
| CN101854797A (zh) * | 2007-06-22 | 2010-10-06 | 孟山都技术公司 | 针对性状表现和表达选择基因座的方法和组合物 |
| US20100313819A1 (en) * | 2003-09-26 | 2010-12-16 | Matsumoto Institute Of Microorganisms Co., Ltd. | Amphidiploid aquatic animal and method of breeding the same |
| CN101979663A (zh) * | 2010-11-26 | 2011-02-23 | 厦门大学 | 日本囊对虾2种形态变异体的识别方法 |
| JP2014213836A (ja) * | 2013-04-30 | 2014-11-17 | 株式会社ブリヂストン | 乗用車用空気入りラジアルタイヤ |
| CN104686430A (zh) * | 2015-04-07 | 2015-06-10 | 雷州市海威水产养殖有限公司 | 一种对虾耐低溶氧家系选育方法 |
| RU2014110321A (ru) * | 2014-03-18 | 2015-09-27 | Государственное научное учреждение Всероссийский научно-исследовательский институт генетики и разведения сельскохозяйственных животных Российской академии сельскохозяйственных наук | Способ геномной селекции крупного рогатого скота |
| CN104998277A (zh) * | 2015-05-21 | 2015-10-28 | 华南农业大学 | 一种抗弓形虫组方药物及其筛选方法 |
| CN105918173A (zh) * | 2016-05-03 | 2016-09-07 | 中国水产科学研究院黄海水产研究所 | 一种三疣梭子蟹家系高效选育方法 |
| CN106480189A (zh) * | 2016-10-18 | 2017-03-08 | 中国水产科学研究院黄海水产研究所 | 一种基于全基因组选择的鱼类抗病良种培育方法 |
| US20190161780A1 (en) * | 2016-08-10 | 2019-05-30 | Ajinomoto Co., Inc. | Production Method for L-Amino Acid |
| CA2998211A1 (en) * | 2018-03-16 | 2019-09-16 | Commonwealth Scientific And Industrial Research Organisation | Plants producing modified levels of medium chain fatty acids |
| CN111128306A (zh) * | 2020-01-06 | 2020-05-08 | 中国水产科学研究院黄海水产研究所 | 一种罗非鱼基因组选择育种方法 |
| AU2020100975A4 (en) * | 2020-06-10 | 2020-08-27 | Beijing Academy Of Agriculture And Forestry Sciences | Snp combination for germplasm resource diversity analysis and molecular breeding of chinese cabbage, and use thereof |
| AU2020102926A4 (en) * | 2020-10-21 | 2020-12-17 | Yellow Sea Fisheries Research Institute, Chinese Academy Of Fishery Sciences | A method for culturing Panulirus ornatus in cages of seawater pond |
| CN112164424A (zh) * | 2020-08-03 | 2021-01-01 | 南京派森诺基因科技有限公司 | 一种基于无参考基因组的群体进化分析方法 |
| CN112273291A (zh) * | 2020-10-28 | 2021-01-29 | 厦门大学 | 基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法 |
-
2021
- 2021-08-10 CN CN202110914524.1A patent/CN113678767B/zh active Active
Patent Citations (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2007200633A1 (en) * | 2000-07-27 | 2007-03-08 | The Johns Hopkins University | Growth differentiation factor receptors, agonists and antagonists thereof, and methods of using same |
| US20060150273A1 (en) * | 2003-06-03 | 2006-07-06 | Tapper Brian A | Grass endophytes |
| US20100313819A1 (en) * | 2003-09-26 | 2010-12-16 | Matsumoto Institute Of Microorganisms Co., Ltd. | Amphidiploid aquatic animal and method of breeding the same |
| CN101016519A (zh) * | 2006-09-29 | 2007-08-15 | 武汉大学 | 对虾抗病毒促生长双功能工程菌株及构建方法和应用 |
| CN101011037A (zh) * | 2007-01-27 | 2007-08-08 | 中国水产科学研究院黄海水产研究所 | 鲈鱼冷冻精子诱导半滑舌鳎卵减数雌核发育方法 |
| CN101854797A (zh) * | 2007-06-22 | 2010-10-06 | 孟山都技术公司 | 针对性状表现和表达选择基因座的方法和组合物 |
| CN101120660A (zh) * | 2007-08-28 | 2008-02-13 | 中国水产科学研究院黄海水产研究所 | 鱼类异源冷冻精子诱导雌核发育的方法 |
| CN101185426A (zh) * | 2007-11-20 | 2008-05-28 | 中国水产科学研究院黄海水产研究所 | 无特定病原中国对虾零换水苗种生产方法 |
| CN101215565A (zh) * | 2008-01-11 | 2008-07-09 | 华中农业大学 | 一种作为猪标记辅助选择的与免疫性状相关的分子标记的克隆及应用 |
| CN101292635A (zh) * | 2008-06-11 | 2008-10-29 | 大连水产学院 | 菲律宾蛤仔杂交系的制种方法 |
| US20100037330A1 (en) * | 2008-06-11 | 2010-02-11 | Kannika Siripattarapravat | Efficient Somatic Cell Nuclear Transfer In Fish |
| CN101317554A (zh) * | 2008-06-11 | 2008-12-10 | 大连水产学院 | 以放射肋为标记的菲律宾蛤仔制种工艺 |
| CN101386857A (zh) * | 2008-10-21 | 2009-03-18 | 华中农业大学 | 一种与番茄抗病性状相关的分子标记克隆及应用 |
| CN101433188A (zh) * | 2008-12-10 | 2009-05-20 | 中国水产科学研究院黄海水产研究所 | 罗非鱼选择育种个体标记方法 |
| WO2010068999A1 (en) * | 2008-12-19 | 2010-06-24 | Genetysis Pty Ltd | Method for identifying an animal suitable for breeding |
| CN101979663A (zh) * | 2010-11-26 | 2011-02-23 | 厦门大学 | 日本囊对虾2种形态变异体的识别方法 |
| JP2014213836A (ja) * | 2013-04-30 | 2014-11-17 | 株式会社ブリヂストン | 乗用車用空気入りラジアルタイヤ |
| RU2014110321A (ru) * | 2014-03-18 | 2015-09-27 | Государственное научное учреждение Всероссийский научно-исследовательский институт генетики и разведения сельскохозяйственных животных Российской академии сельскохозяйственных наук | Способ геномной селекции крупного рогатого скота |
| CN104686430A (zh) * | 2015-04-07 | 2015-06-10 | 雷州市海威水产养殖有限公司 | 一种对虾耐低溶氧家系选育方法 |
| CN104998277A (zh) * | 2015-05-21 | 2015-10-28 | 华南农业大学 | 一种抗弓形虫组方药物及其筛选方法 |
| CN105918173A (zh) * | 2016-05-03 | 2016-09-07 | 中国水产科学研究院黄海水产研究所 | 一种三疣梭子蟹家系高效选育方法 |
| US20190161780A1 (en) * | 2016-08-10 | 2019-05-30 | Ajinomoto Co., Inc. | Production Method for L-Amino Acid |
| CN106480189A (zh) * | 2016-10-18 | 2017-03-08 | 中国水产科学研究院黄海水产研究所 | 一种基于全基因组选择的鱼类抗病良种培育方法 |
| CA2998211A1 (en) * | 2018-03-16 | 2019-09-16 | Commonwealth Scientific And Industrial Research Organisation | Plants producing modified levels of medium chain fatty acids |
| CN111128306A (zh) * | 2020-01-06 | 2020-05-08 | 中国水产科学研究院黄海水产研究所 | 一种罗非鱼基因组选择育种方法 |
| AU2020100975A4 (en) * | 2020-06-10 | 2020-08-27 | Beijing Academy Of Agriculture And Forestry Sciences | Snp combination for germplasm resource diversity analysis and molecular breeding of chinese cabbage, and use thereof |
| CN112164424A (zh) * | 2020-08-03 | 2021-01-01 | 南京派森诺基因科技有限公司 | 一种基于无参考基因组的群体进化分析方法 |
| AU2020102926A4 (en) * | 2020-10-21 | 2020-12-17 | Yellow Sea Fisheries Research Institute, Chinese Academy Of Fishery Sciences | A method for culturing Panulirus ornatus in cages of seawater pond |
| CN112273291A (zh) * | 2020-10-28 | 2021-01-29 | 厦门大学 | 基于全基因组选择的大黄鱼抗刺激隐核虫病的选育方法 |
Non-Patent Citations (6)
| Title |
|---|
| TRYGVE GJEDREM: "Disease Resistant Fish and Shellfish Are Within Reach : A Review", 《MARINE SCIENCE AND ENGINEERING》 * |
| 孔杰: "中国对象抗病品种的选育", 《青岛 全国水产养殖研讨会论文集》 * |
| 宁岳等: "牡蛎育种研究进展", 《厦门大学学报（自然科学版）》 * |
| 李健: "《中国对虾和三疣梭子蟹遗传育种》", 30 June 2016, 中国海洋大学出版社 * |
| 蔡春波等: "运用SNP芯片评估马身猪保种群体的遗传结构", 《畜牧兽医学报》 * |
| 蔡生力: "《水产养殖学概论》", 30 November 2015, 海洋出版社 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113678767B (zh) | 2022-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gutierrez et al. | Genomic selection for growth traits in Pacific oyster (Crassostrea gigas): potential of low-density marker panels for breeding value prediction | |
| Gjøen et al. | Past, present, and future of genetic improvement in salmon aquaculture | |
| CN113699250B (zh) | 与肉鸡饲料转化效率性状相关的分子标记及其应用 | |
| US11690361B2 (en) | Transition cow index | |
| EP3882360B1 (en) | Genomic selection method of huaxi cattle | |
| Bigi et al. | Genetic characterization of four native Italian shepherd dog breeds and analysis of their relationship to cosmopolitan dog breeds using microsatellite markers | |
| CN112002371A (zh) | 一种白羽肉鸡剩余采食量的基因组选择方法 | |
| Brown | Genetic management and selective breeding in farmed populations of gilthead seabream (Sparus aurata) | |
| KR101135443B1 (ko) | 성장이 빠른 육종 넙치의 육종 방법 | |
| CN102382878B (zh) | 一种鉴别建鲤不同家系的分子标记方法 | |
| CN107711685A (zh) | 优质高抗快繁新品鸡基因筛选及组合方法 | |
| Jerry et al. | Application of DNA parentage analyses for determining relative growth rates of Penaeus japonicus families reared in commercial ponds | |
| Phavaphutanon et al. | Evaluation of quantitative trait loci for hip dysplasia in Labrador Retrievers | |
| CN113678767B (zh) | 一种对虾抗病性状的选育方法 | |
| Sukhavachana et al. | Improving genomic prediction accuracy for harvest traits in Asian seabass (Lates calcarifer, Bloch 1790) via marker selection | |
| CN113549699B (zh) | 一种白羽肉鸡产蛋数的基因组选择方法 | |
| CN110771547B (zh) | 一种基于Gompertz模型的卵形鲳鲹体质量育种方法 | |
| CN112021224B (zh) | 一种基于扩繁群体的虾类选择育种方法 | |
| US20220307092A1 (en) | Method for predicting survival following streptococcus iniae infection | |
| CN105002293A (zh) | 一种利用微卫星标记鉴定混养松浦镜鲤家系的方法 | |
| Mengistu et al. | interaction between aerated and non-aerated ponds and the impact of aeration on genetic parameters in Nile tilapia (Oreochromis niloticus) | |
| Bosse | The hybrid nature of pig genomes: unraveling the mosaic haplotype structure in wild and commercial Sus scrofa populations | |
| Metallinos | Canine molecular genetic testing | |
| Le Cordeur | Genetic and phenotypic characterisation of commercial dusky kob (Argyrosomus japonicus) cohorts | |
| Kristjánsson et al. | Using GBLUP selection for increased Piscirickettsia salmonis (SRS) resistance in Atlantic salmon (Salmo salar) |
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 |