CN106755441B - 一种基于多性状的基因组选择进行林木多性状聚合育种的方法 - Google Patents

一种基于多性状的基因组选择进行林木多性状聚合育种的方法 Download PDF

Info

Publication number
CN106755441B
CN106755441B CN201611247628.7A CN201611247628A CN106755441B CN 106755441 B CN106755441 B CN 106755441B CN 201611247628 A CN201611247628 A CN 201611247628A CN 106755441 B CN106755441 B CN 106755441B
Authority
CN
China
Prior art keywords
trait
forest
breeding
clone
snp
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
Application number
CN201611247628.7A
Other languages
English (en)
Other versions
CN106755441A (zh
Inventor
林元震
张卫华
张心菲
程玲
张鑫鑫
周玮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China Agricultural University
Original Assignee
South China Agricultural University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by South China Agricultural University filed Critical South China Agricultural University
Priority to CN201611247628.7A priority Critical patent/CN106755441B/zh
Publication of CN106755441A publication Critical patent/CN106755441A/zh
Application granted granted Critical
Publication of CN106755441B publication Critical patent/CN106755441B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • C12Q1/6895—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for plants, fungi or algae
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • 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/13—Plant traits
    • 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/156—Polymorphic or mutational markers

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biotechnology (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Organic Chemistry (AREA)
  • Genetics & Genomics (AREA)
  • Wood Science & Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Zoology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Theoretical Computer Science (AREA)
  • Mycology (AREA)
  • Botany (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Immunology (AREA)
  • Microbiology (AREA)
  • Medical Informatics (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Biochemistry (AREA)
  • General Engineering & Computer Science (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明公开了一种基于多性状的基因组选择进行林木多性状聚合育种的方法。所述方法依次包括：无性系个体表型数据采集；无性系谱系信息及A矩阵构建，或SNP标记分型数据及G矩阵构建；多性状模型建立和数据分析。本发明的方法可减少林木性状表型测定的工作量，且不受性状与标记关联程度的限制，与常规的传统子代测定技术相比，可实现多性状的定向、精确育种，获得的杂交后代具有可靠的生产力和遗传背景，显著缩短了林木多性状聚合育种周期，且对树种和试验林类型无具体限制，可快速达到林木多性状聚合育种的目的，具有广泛的应用前景。

Description

一种基于多性状的基因组选择进行林木多性状聚合育种的 方法
技术领域
本发明属于林木遗传育种技术领域，具体地涉及一种基于多性状的基因组选择进行林木多性状聚合育种的方法。
背景技术
传统的林木遗传改良需要通过田间对比试验和性状测定，估算各种遗传参数和评价种源、家系、单株或无性系的表现，来培育符合林业生产需求的林木良种。然而林木的生产性状大多属于微效多基因控制，而且林木开花、授粉及子代测定，往往周期比较长，一般在数十年，因而，要进行林木多性状的聚合育种，周期更长，尤其是当性状比较多且性状相关比较复杂时。漫长的世代时间以及大多数复杂性状晚期表达的特性历来都是林木育种所面临的巨大困难和挑战。因此，有必要研发探索一种高效地林木多性状聚合育种的方法。
基因组选择(Genomic Selection，GS)技术是一种基于动植物育种群体，利用芯片技术采集群体个体高密度基因组SNP的基因分型数据，然后构建个体间的基因组关系矩阵，通过GBLUP法获取个体全基因组育种值的新型分子育种技术。它与传统育种技术相比，优点在于可以突破遗传材料生长时间限制，在幼期即可估计个体的基因组育种值，并据此进行早期预测与选择，从而显著缩短育种周期。目前，基因组选择已广泛应用于各国畜禽分子育种中，尤其在奶牛、猪和鸡育种方面。但是目前，还未见有基于多性状的基因组选择进行林木多性状聚合育种的报道。
发明内容
本发明所要解决的技术问题是克服现有技术中林木多性状聚合育种的缺陷和不足，提供一种高效的基于多性状的基因组选择进行林木多性状聚合育种林木多性状聚合育种方法，以解决利用现有的林木育种技术，在进行多性状聚合育种时，不能精准选择亲本以及所需周期长的问题，从而弥补现有技术的不足。该方法是将基因组选择技术与多性状聚合技术相结合形成的一套进行林木多性状聚合育种的新方法体系，无论对于缩短林木良种选育周期，还是对于开展林木多世代育种，都具有非常重要的意义。
本发明的目的是提供一种基于多性状的基因组选择进行林木多性状聚合育种的方法。
本发明的上述目的是通过以下技术方案给予实现的：
一种基于多性状的基因组选择进行林木多性状聚合育种的方法，所述方法具体包括如下步骤：
S1：在无性系测定林中，随机抽样若干个(优选为20～40个，最优选30个)无性系，每个无性系选取若干株(优选为4～7株，最优选5株)，分别测定其表型性状，然后进行性状聚类分析，按正负相关归类，得到性状的归类信息；
S2：根据性状的归类信息，每类各选取容易测定的一个代表性状，对所有参试的无性系进行相应表型测定，得到多性状的表型数据；
S3：根据无性系的来源亲本，构建无性系的谱系数据，构建加性相关矩阵A及其逆矩阵；或对每个无性系进行SNP芯片基因分型，得到无性系的基因分型数据，数据记为0、1、2的SNP标记矩阵，根据无性系分型的SNP标记矩阵，构建无性系的G矩阵；
S4：构建模型，利用REML法对表型数据进行多性状分析，估算环境误差和遗传方差，分析公式为：
y＝Xb+Zu+e
其中，y代表无性系表型测量值，b代表固定效应值，u代表随机效应值，e代表残差效应。
这个模型为通用线性模型，具体要根据数据类型而定。如实施例中给出的两种模型：
模型1，BLUP模型，包含了表型值y，总体均值μ(固定效应值)，无性系Clone(随机效应值)和残差e，其中无性系随机效应的矩阵为加性相关矩阵A。
模型2，GBLUP模型，包含了表型值y，总体均值μ(固定效应值)，无性系Clone(随机效应值)和残差e，其中无性系随机效应的矩阵为基因组关系矩阵G。
S5：根据步骤S4模型获取无性系所有个体(包括亲本)不同性状的全基因组育种值，分别进行排名，然后各自选取全基因组育种值排名靠前的个体作为多性状聚合育种的亲本进行杂交。
本发明所述方法对于试验材料并无树种和试验林类型的具体限制，因此其他树种以及子代测定林和无性系试验林都不受限于此，通过本发明即可快速达到林木多性状聚合育种的目的，具有广泛的应用前景。
优选地，步骤S1所述无性系测定林为桉树无性系试验林。
优选地，步骤S1所述表型性状为生长性状、材性性状或抗性性状中的一个或多个。
更优选地，所述生长性状为树高h、胸径dbh、材积v以及木材吸水率wpro中的一个或多个；所述材性性状为心材比例cpro、木材密度wd、管胞长度tl、管胞宽度tw以及管胞长宽比lrt中的一个或多个。
优选地，步骤S3所述谱系数据采用R程序包kinship2生成。
优选地，步骤S3所述SNP分型数据是采用高通量SNP芯片对无性系多态性的基因分型数据集。
优选地，步骤S3是对每个无性系进行SNP芯片基因分型，得到无性系的基因分型数据，数据记为0、1、2的SNP标记矩阵，根据无性系分型的SNP标记矩阵，构建无性系的G矩阵。实验表明，基因组关系矩阵比加性相关矩阵更能捕获个体间的遗传关系，以此为基础构建的模型2，即利用SNP标记信息构建G矩阵的多性状模型为最优模型。
更优选地，所述SNP分型数据采集的具体步骤为：
提取所有无性系的DNA，采用SNP芯片进行基因型测定。读取SNP芯片数据文件，用Beagle程序进行缺失基因型填补，然后选取SNP最小等位基因频率(MAF)0.01作为质量控制标准，剔除最小等位基因频率(MAF)低于0.01和染色体位置未知的SNP。
再优选地，所述SNP芯片为EucHIP60K.BR(含60728个SNP标记)。
优选地，步骤S4所述构建模型为构建GBLUP模型。
优选地，步骤S4中所述多性状分析为利用SNP标记信息构建G矩阵的多性状GBLUP模型，利用ASReml-R软件对试验数据进行分析，同时结合S3构建的G矩阵，从而估算遗传方差和环境误差，估算遗传方差和环境误差。
同时，上述基于多性状的基因组选择进行林木多性状聚合育种的方法在林木遗传育种上的应用亦在本发明保护范围内。
优选地，为在林木多性状聚合育种上的应用。
与现有技术相比，本发明具有以下有益效果：
本发明提供了一种将基因组选择技术与多性状聚合技术相结合形成的一套进行林木多性状聚合育种的新方法体系，通过林木性状聚类分析、林木杂交子代SNP分型和多性状的基因组选择分析进行林木多性状聚合育种，具有以下优点：
(1)本发明所述方法减少了林木性状表型测定的工作量；
(2)本发明所述方法不受性状与标记关联程度的限制；
(2)本发明所述方法与常规的传统子代测定技术相比，可实现多性状的定向、精确育种；
(4)本发明所述方法获得的杂交后代具有可靠的生产力和遗传背景；
(5)本发明所述方法可显著缩短林木多性状聚合育种周期；
(6)本发明所述方法对于树种和试验林类型没有具体限制，各种树种以及子代测定林和无性系试验林都不受限于此，通过本发明即可快速达到林木多性状聚合育种的目的，具有广泛的应用前景。
附图说明
图1为本实施例1桉树无性系不同性状的聚类图。
图2为本实施例1基于谱系信息的桉树无性系加性相关矩阵图。
图3为本实施例1基于SNP分型信息的桉树无性系基因组关系矩阵图。
具体实施方式
下面结合说明书附图和具体实施例对本发明作出进一步地详细阐述，所述实施例只用于解释本发明，并非用于限定本发明的范围。下述实施例中所使用的试验方法如无特殊说明，均为常规方法；所使用的材料、试剂等，如无特殊说明，为可从商业途径得到的试剂和材料。
本发明的方法可以适用于任何已有的试验设计，利用SNP标记信息构建的基因组关系矩阵G开展林木基因组选择分析，可针对生长、材性、抗性等性状，进行林木多性状聚合育种，尤其适用于杂交子代测定林或无性系试验林的亲本的反向选择。可以满足任意生长阶段的林木多性状聚合育种。
下面对本发明的方法做进一步的描述。
实施例1
本发明的林木多性状聚合育种方法，包括无性系个体表型数据采集、无性系谱系信息及A矩阵构建、SNP标记分型数据及G矩阵构建、多性状模型建立和数据分析步骤；具体如下：
1.无性系个体表型数据采集
(1)试验材料为：桉树无性系试验林，参试无性系共165个，每个无性系大约40株。
抽样30个桉树无性系，每个无性系5株，分别测定树高(h)、胸径(dbh)、材积(v)、心材比例(cpro)、木材密度(wd)、木材吸水率(wpro)、管胞长度(tl)和管胞宽度(tw)及管胞长宽比(lrt)，然后进行性状聚类分析，按正负相关归类，得到性状的两大归类信息；
(2)根据性状的归类信息，每类选取容易测定的一个代表性状，对所有参试的无性系进行相应表型测定，得到多性状的表型数据。
2.无性系谱系信息及A矩阵构建
收集无性系号及其亲本号，采用R程序包kinship2生成谱系，并构建加性相关矩阵A。利用ASReml-R直接构建A的逆矩阵。
3.SNP分型数据及G矩阵构建
SNP分型数据是采用高通量SNP芯片对无性系多态性的基因分型数据集，用于计算无性系间的基因组相关矩阵。
具体构建方法：须提取所有无性系的DNA，采用SNP芯片EucHIP60K.BR(含60728个SNP标记)进行基因型测定。读取SNP芯片数据文件，用Beagle程序进行缺失基因型填补，然后选取SNP最小等位基因频率(MAF)0.01作为质量控制标准，剔除最小等位基因频率(MAF)低于0.01和染色体位置未知的SNP，最终41160个SNP用于分析。G矩阵及其逆矩阵，参照文献Isik等人(2016)的方法(Isik F,Bartholomé J,Farjat A,Chancerel E,Raffin A,Sanchez L,Plomion C,Bouffier L.Genomic selection in maritime pine.PlantScience，2016,242:108–119)。
4.多性状模型建立
为了实现本发明目的，本发明的方案设置了2种分析模型进行了比较，模型方案如下：
模型1，BLUP模型，包含了表型值y，总体均值μ(固定效应)，无性系Clone(随机效应)和残差e，其中无性系随机效应的矩阵为加性相关矩阵A。
模型2，GBLUP模型，包含了表型值y，总体均值μ(固定效应)，无性系Clone(随机效应)和残差e，其中无性系随机效应的矩阵为基因组关系矩阵G。
5.数据分析
(1)由图1可知，抽样测定的9个性状，可分为2类：h、dbh、v和wpro聚为一类，其它5个聚为一类。
第一类性状主要为生长性状，选择胸径dbh；第二类性状为材性性状，选择木材密度wd，作为后续所有无性系的表型测定性状。
(2)图2中，颜色越偏红色，表示个体间的亲缘关系越弱，颜色越偏黄色，表示个体间的亲缘关系越强。因此，基于谱系信息的无性系间加性相关矩阵结果显示，大部分无性系及其亲本之间的亲缘关系比较弱。
(3)图3的颜色图与图2一样，但基于SNP分型信息的无性系间基因组关系矩阵图，明显与谱系信息的加性相关矩阵图不同，前者明显可以看出无性系及其亲本间存在比较强的亲缘关系，表明基因组关系矩阵比加性相关矩阵更能捕获个体间的遗传关系。
(4)基于上述两种模型，利用ASReml-R软件分别估算分差组分，具体结果见表1。
对于模型1(BLUP模型)，性状1的无性系遗传方差Vc为86.83，残差Ve为181.85，性状2的无性系遗传方差Vc为397.34，残差Ve为811.15。
对于模型2(GBLUP模型)，性状1的无性系遗传方差Vc为98.32，残差Ve为181.86，性状2的无性系遗传方差Vc为466.63，残差Ve为811.66。
与模型1比较，两个性状，模型2的残差基本没有变化，但无性系遗传方差都增加了，可见模型2比模型1拟合的结果更好，即基因组选择技术采用的基因组关系矩阵更能捕获个体间的遗传关系。
表1 模型遗传参数估计值
Figure BDA0001197439450000061
注：Vc为无性系遗传方差，Ve为误差。
此外，输出了参数无性系所有个体的两个性状的全基因组育种值(GEBV)和传统育种值(ABV)，只输出GEBV前10％的个体，如表2所示。从表2可知，BLUP模型和GBLUP模型，所获得个体的全基因组育种值GEBV基本上都小于传统育种值ABV，即传统的BLUP模型存在高估育种值的情况。基于上述多性状模型分析的结果可知，GBLUP模型是优选模型，如果每个性状只取最好的10个个体作为亲本，则对于性状1，亲本应选1169、185、1131、1153、1107、1167、1172、1106和p20，对于性状2，亲本应选择163、172、1172、1123、165、151、p20、164、1154和147。值得关注的是，在本轮分析中，原有的12个亲本，仅有亲本P20号保留下来了，表明其遗传基因确实不错。此外，个体1172和亲本P20在两个性状全基因组育种值均排名靠前，而之前的性状聚类分析结果显示，性状1和性状2之间应呈现负相关，因此，个体1172和亲本120应属于相关破坏者，在多性状聚合育种中，对于这类个体，应给予特别关注。
表2 模型遗传参数估计值
Figure BDA0001197439450000071
(5)基于上述分析可知，模型2，即利用SNP标记信息构建G矩阵的多性状模型为最优模型，即本发明的最佳分析模型。通过多性状的最佳模型，可捕获遗传材料间更多的遗传关系，进而精确估算遗传误差，并进一步精确估算全基因组育种值，据此来筛选不同性状全基因组育种值高的个体进行杂交，于是快速达到林木多性状聚合育种的目的。
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。

Claims (10)

1.一种基于多性状的基因组选择进行林木多性状聚合育种的方法，其特征在于，所述方法具体包括如下步骤：
S1：在无性系测定林中，随机抽样若干个无性系，每个无性系选取若干株，分别测定其表型性状，然后进行性状聚类分析，按正负相关归类，得到性状的归类信息；
S2：根据性状的归类信息，每类各选取一个代表性状，对所有参试的无性系进行相应表型测定，得到多性状的表型数据；
S3：根据无性系的来源亲本，构建无性系的谱系数据，构建加性相关矩阵A及其逆矩阵；
或对每个无性系进行SNP芯片基因分型，得到无性系的基因分型数据，数据记为0、1、2的SNP标记矩阵，根据无性系分型的SNP标记矩阵，构建无性系的G矩阵；
S4：构建模型，利用REML法对表型数据进行多性状分析，估算环境误差和遗传方差，分析公式为：
y＝Xb+Zu+e
其中，y代表无性系表型测量值，b代表固定效应值，u代表随机效应值，e代表残差效应，X为固定效应的关联矩阵，Z为随机效应的关联矩阵；
S5：根据步骤S4模型获取无性系所有个体不同性状的全基因组育种值，分别进行排名，然后各自选取全基因组育种值排名靠前的个体作为多性状聚合育种的亲本进行杂交。
2.根据权利要求1所述的基于多性状的基因组选择进行林木多性状聚合育种的方法，其特征在于，步骤S1所述无性系测定林为桉树无性系试验林。
3.根据权利要求1所述的基于多性状的基因组选择进行林木多性状聚合育种的方法，其特征在于，步骤S1所述表型性状为生长性状、材性性状或抗性性状中的一个或多个。
4.根据权利要求3所述的基于多性状的基因组选择进行林木多性状聚合育种的方法，其特征在于，所述生长性状为树高h、胸径dbh、材积v或木材吸水率wpro中的一个或多个；所述材性性状为心材比例cpro、木材密度wd、管胞长度tl、管胞宽度tw或管胞长宽比lrt中的一个或多个。
5.根据权利要求1所述的基于多性状的基因组选择进行林木多性状聚合育种的方法，其特征在于，步骤S3所述谱系数据采用R程序包kinship2生成。
6.根据权利要求1所述的基于多性状的基因组选择进行林木多性状聚合育种的方法，其特征在于，步骤S3所述SNP分型数据是采用高通量SNP芯片对无性系多态性的基因分型数据集。
7.根据权利要求6所述的基于多性状的基因组选择进行林木多性状聚合育种的方法，其特征在于，所述SNP分型数据采集的具体步骤为：
提取所有无性系的DNA，采用SNP芯片进行基因型测定；读取SNP芯片数据文件，用Beagle程序进行缺失基因型填补，然后选取SNP最小等位基因频率0.01作为质量控制标准，剔除最小等位基因频率低于0.01和染色体位置未知的SNP。
8.根据权利要求7所述的基于多性状的基因组选择进行林木多性状聚合育种的方法，其特征在于，所述SNP芯片为EucHIP60K.BR。
9.根据权利要求1所述的基于多性状的基因组选择进行林木多性状聚合育种的方法，其特征在于，其特征在于，步骤S4中所述多性状分析为利用SNP标记信息构建G矩阵的多性状GBLUP模型，估算遗传方差和环境误差。
10.权利要求1～9任一所述的基于多性状的基因组选择进行林木多性状聚合育种的方法在林木遗传育种上的应用。
CN201611247628.7A 2016-12-29 2016-12-29 一种基于多性状的基因组选择进行林木多性状聚合育种的方法 Active CN106755441B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611247628.7A CN106755441B (zh) 2016-12-29 2016-12-29 一种基于多性状的基因组选择进行林木多性状聚合育种的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611247628.7A CN106755441B (zh) 2016-12-29 2016-12-29 一种基于多性状的基因组选择进行林木多性状聚合育种的方法

Publications (2)

Publication Number Publication Date
CN106755441A CN106755441A (zh) 2017-05-31
CN106755441B true CN106755441B (zh) 2020-08-07

Family

ID=58929169

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611247628.7A Active CN106755441B (zh) 2016-12-29 2016-12-29 一种基于多性状的基因组选择进行林木多性状聚合育种的方法

Country Status (1)

Country Link
CN (1) CN106755441B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108157293B (zh) * 2018-02-07 2020-03-20 山东省农业科学院奶牛研究中心 一种基于系谱信息简化选择高生产性能a2a2纯合基因型奶牛的育种方法
CN108371105B (zh) * 2018-03-16 2019-10-25 广东省农业科学院水稻研究所 一种基于核心系谱品种的高密度分子标记辅助聚合育种方法
CN109101786B (zh) * 2018-08-29 2021-02-09 广东省农业科学院动物科学研究所 一种整合显性效应的基因组育种值估计方法
CN109727642B (zh) * 2019-01-22 2021-03-23 隆平农业发展股份有限公司 基于随机森林模型的全基因组预测方法及装置
CN109706265B (zh) * 2019-03-11 2022-04-01 中国林业科学研究院热带林业研究所 一种snp组合及其在桉树无性系鉴定中的应用
CN111524545B (zh) * 2020-04-30 2023-11-10 天津诺禾致源生物信息科技有限公司 全基因组选择育种的方法和装置
CN112331263B (zh) * 2020-10-22 2021-07-23 华南农业大学 一种基于个体遗传竞争与环境空间分析的林木基因组选择方法及其应用
CN116064846A (zh) * 2023-01-30 2023-05-05 中国海洋大学三亚海洋研究院 一种评估花鲈生长和抗性性状综合育种值的方法及应用

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005015989A1 (en) * 2003-08-04 2005-02-24 Monsanto Technology Llc Method for genetic improvement of terminal boars
CN103914631A (zh) * 2014-02-26 2014-07-09 中国农业大学 一种基于snp芯片的综合基因组育种值估计方法及应用
CN103914632A (zh) * 2014-02-26 2014-07-09 中国农业大学 一种快速估计基因组育种值的方法和应用
CN104273023A (zh) * 2014-05-15 2015-01-14 浙江大学 一种水稻聚合育种方法
CN104855311A (zh) * 2015-05-08 2015-08-26 大连海洋大学 壳宽型菲律宾蛤仔新品系制种方法
CN106022005A (zh) * 2016-05-21 2016-10-12 安徽省农业科学院畜牧兽医研究所 一种连续性状和阈性状基因组育种值联合估计的贝叶斯方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005015989A1 (en) * 2003-08-04 2005-02-24 Monsanto Technology Llc Method for genetic improvement of terminal boars
CN103914631A (zh) * 2014-02-26 2014-07-09 中国农业大学 一种基于snp芯片的综合基因组育种值估计方法及应用
CN103914632A (zh) * 2014-02-26 2014-07-09 中国农业大学 一种快速估计基因组育种值的方法和应用
CN104273023A (zh) * 2014-05-15 2015-01-14 浙江大学 一种水稻聚合育种方法
CN104855311A (zh) * 2015-05-08 2015-08-26 大连海洋大学 壳宽型菲律宾蛤仔新品系制种方法
CN106022005A (zh) * 2016-05-21 2016-10-12 安徽省农业科学院畜牧兽医研究所 一种连续性状和阈性状基因组育种值联合估计的贝叶斯方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
"Comparative analysis of the GBLUP,emBayesB, and GWAS algorithms to predict genetic values in large yellow croaker (Larimichthys crocea)";Linsong Dong et al.;《BMC Genomics》;20160614;第17卷;第1-10页 *
"Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa";Keyan Zhao et al.;《NATURE COMMUNICATIONS》;20110913;第2卷;第1-11页 *
"利用紧缩线性模型和贝叶斯模型对猪总产仔数和产活仔数性状的全基因组关联研究";刘小磊等;《遗传 HEREDITAS》;20120910;第34卷(第10期);第1261-1270页 *
"植物关联分析方法的研究进展";冯建英等;《作物学报》;20160511;第42卷(第7期);第945-956页 *
"油松优良家系多性状选择方法研究";刘永红等;《西北农林科技大学学报》;20061231;第34卷(第12期);第115-120页 *

Also Published As

Publication number Publication date
CN106755441A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106755441B (zh) 一种基于多性状的基因组选择进行林木多性状聚合育种的方法
CN104532359B (zh) 玉米DNA指纹库构建及品种分子鉴定SNP核心位点组合-maizeSNP384
CN112695107B (zh) 一种肉用绵羊生长性能snp位点组合及其应用
US20210269887A1 (en) Method and application for rapid and accurate chromosomal location of economic traits in laver
CN111088382A (zh) 一种玉米全基因组snp芯片及其应用
CN105010233A (zh) 一种应用snp辅助选择育种技术培育高繁殖性能种兔的方法
CN108319984A (zh) 基于dna甲基化水平的木本植物叶片表型特征和光合特性预测模型的构建方法及预测方法
CN106755300A (zh) 一种识别猕猴桃杂交亲本对子代基因组贡献比例的方法
Liang et al. Molecular identification and genetic analysis of cherry cultivars using capillary electrophoresis with fluorescence-labeled SSR markers
CN110444253B (zh) 一种适用于混池基因定位的方法及系统
CN101613742B (zh) 一种中华绒螫蟹多元高通量遗传标记系统及遗传分析方法
CN107034306B (zh) 一组snp位点基因分型引物及其在小麦品种鉴定中的应用
CN106498070A (zh) 一种基于基因组LoF位点筛查间接关联猕猴桃种性的方法
CN113755638A (zh) 一种利用srap分子标记快速鉴别芦笋品种的方法
CN117133354B (zh) 一种高效鉴定林木关键育种基因模块的方法
CN110853711B (zh) 一种预测烟草果糖含量的全基因组选择模型及其应用
CN117210597A (zh) 基于snp标记的菊花dna指纹图谱及其构建方法和应用
CN115141893B (zh) 包含7个分子标记的预测猕猴桃果实干物质含量的分子标记组及其应用和试剂盒
CN112331263B (zh) 一种基于个体遗传竞争与环境空间分析的林木基因组选择方法及其应用
CN107535350B (zh) 一种基于ssr标记的玉米组配模式优选方法及系统
CN110853710B (zh) 一种预测烟草淀粉含量的全基因组选择模型及其应用
CN104573409B (zh) 基因定位的多重检验方法
CN116467596B (zh) 水稻粒长预测模型的训练方法、形态预测方法及装置
CN115216526A (zh) 一种基于三体组的基因组特征来评估植物杂种的基因组构成的方法
CN114863991B (zh) 基于二步预测模型建立的提高全基因组预测精度方法

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