棉花产量构成因素性状的全基因组关联分析

【目的】 对铃重、衣分、单株铃数和籽指等棉花产量构成因素性状进行全基因组关联分析（genome-wide association study,GWAS）,发掘与其关联的标记位点、优异等位变异及候选基因,为棉花产量的分子育种提供理论依据。【方法】 以408份陆地棉品种（系）资源为材料,利用Cotton SNP 80K芯片,对6个环境的铃重、衣分、单株铃数和籽指4个产量构成因素性状进行基于混合线性模型（mixed linear model,MLM）的全基因组关联分析,检测与产量构成因素性状显著关联的位点、优异等位变异;进一步依据转录组数据的基因表达量,在显著关联的位点侧翼序列1 Mb区间挖掘可能的候选基因。【结果】 4个产量构成因素性状在不同环境下均表现出广泛的表型变异,其中,单株铃数变异系数最大为16.67%—22.66%,各性状的遗传率为48.4%—92.2%;除铃重与衣分间相关性不显著外,其他性状间均呈显著或极显著相关性;基于6个环境各性状表型数据的最佳线性无偏预测值（best linear unbiased prediction,BLUP）,GWAS共检测到分布于基因组的7个区间内23个与目标性状关联的SNP位点,其中,与铃重关联的位点5个,与衣分关联的位点1个,与单株铃数关联的位点9个,与籽指关联的位点8个,有3个位点（TM21094、TM21102和TM57382）同时与多个目标性状关联;鉴定到7个最优SNP位点的优异等位变异,分别为TM21099（TT）、TM57382（GG）、TM78920（CC）、TM53448（TT）、TM59015（AA）、TM43412（GG）和TM69770（AA）;利用转录组数据分析,在基因组的7个区间筛选到158个与产量形成可能的候选基因,GO富集分析和KEGG代谢途径分析发现,候选基因功能类别多样并参与了多种代谢途径。【结论】 在陆地棉品种（系）群体中共鉴定到23个与产量构成因素性状关联的SNP位点,筛选到158个可能与产量性状相关的候选基因。     

玉米杂交种穗部性状的全基因组关联分析

【目的】玉米穗部性状是产量的重要构成因子,利用全基因组关联分析（genome-wide association study,GWAS）方法解析玉米杂交种穗部性状的遗传基础、挖掘与穗部性状相关的位点,为功能基因克隆和高产玉米品种培育提供参考。【方法】选用115份来源于陕A群和陕B群的优良玉米自交系和4份国内骨干作为亲本,以基于NCⅡ遗传交配设计获得的442份玉米杂交种为材料构建关联群体,调查2个环境中群体材料的穗长、穗粗、穗行数等8个穗部性状;利用tGBS技术检测亲本基因型,推测出F₁杂交种的19 461个高质量SNP,结合杂交种表型和基因型开展基于加性、显性及上位性模型的穗部性状的全基因组关联分析,并利用公共数据库中玉米穗发育相关组织的转录组数据和基因的注释信息预测候选基因。【结果】表型数据分析结果显示,试验群体的8个穗部性状均符合正态分布,表型变异为3.78%—45.25%。方差分析表明,8个穗部性状的环境效应和基因型效应均呈现极显著水平（P<0.001）,广义遗传力为54.15%—68.89%。同时玉米杂交种穗部性状间呈现显著正相关或显著负相关。利用加性和显性模型分别检测到16个和3个显著SNP,上位性模型检测到79个上位性位点。3种模型检测的显著位点累积解释各性状38.21%—60.69%的表型变异,其中,加性模型检测到的显著SNP累积解释的表型变异为0.00—41.26%,上位性模型检测到的位点累积解释的表型变异为15.18%—45.36%。基于加性和显性模型检测的显著SNP的效应分析发现多数位点呈现加性和部分显性效应,仅2个为超显性。进一步分析发现,7个单SNP和5个上位性位点能够解释5%以上的表型变异。根据SNP的位置以及基因的表达信息预测了17个候选基因。【结论】玉米杂交种穗部性状主要受加性、上位性效应影响,显性效应影响较小;加性和显性模型检测的SNP主要表现为加性和部分显性效应,可通过聚合有利等位基因改良目标性状。     

大豆抗倒伏相关性状全基因组关联分析

研究以150份大豆种质为试验材料开展2年3点试验,在收获期测定抗倒伏相关性状,并利用RTM-GWAS方法作关联分析,共检测到99个显著关联SNP位点,其中48个效应显著位点,解释0.28％～3.38％表型变异;89个与环境互作显著位点,解释0.53％～7.04％表型变异;其中,21个位点主效表型变异解释率高于1％;与6个倒伏相关性状显著关联SNP位点中,与茎粗显著关联SNP位点最多,与抗倒指数显著关联SNP位点最少;获得22个有功能注释基因,根据基因表达量和基因功能注释,推测3个基因(Glyma.18G149700、Glyma.04G244100、Glyma.18G011400)可作为与大豆抗倒伏相关候选基因.     

基于GWAS分析的马铃薯主茎数和单株结薯数候选基因挖掘研究

【目的】筛选与马铃薯主茎数和单株结薯数显著关联的单核苷酸多态性(SNP)位点并挖掘候选基因，为高产马铃薯分子育种提供基因资源。【方法】以251份马铃薯核心种质为材料，于2018-2019年在马铃薯主产区4个试验点共8个环境下对其主茎数和单株结薯数进行表型鉴定；提取251份马铃薯种质的DNA，并采用高通量IlluminaHiSeq^TM对DNA进行测序及处理，获得高质量SNP标记，在此基础上，结合表型性状进行全基因组关联（GWAS）分析，对显著关联的SNP位点所在区域的候选基因通过NR、Swiss-port、KEGG、GO等4个数据库进行功能注释。【结果】2018-2019年的检测结果表明，在环境和基因型互作下，251份马铃薯种质间主茎数与单株结薯数有广泛的表型变异。共获得1 209 969个高质量SNP位点，其中与主茎数显著关联的SNP位点有7个，共挖掘到19个候选基因，其中在3个以上环境重叠的候选基因有9个，其可能编码半乳糖醛酸转移酶、转录因子MYB、赤霉素3-β-双加氧酶及细胞色素P450；与单株结薯数显著关联的SNP位点有13个，共挖掘到33个候选基因，去掉2年间重叠的4个候选基因，实际共关联到29个候选基因，其中在3个以上环境下重叠的候选基因有11个，其可能编码果糖-1,6-二磷酸醛缩酶、生长素响应因子IAA₁₃、转录因子WRKY、bHLH113及MADS-box。【结论】在2年8个环境下共鉴定出20个与马铃薯主茎数和单株结薯数显著相关的SNP位点，挖掘到48个候选基因。     

谷子生育期及穗相关性状的QTL定位

【目的】谷子生育期及穗部性状是影响谷子品种适应性及产量的关键因素。通过对相关性状进行QTL定位分析,为探明谷子复杂产量性状的分子遗传机制奠定基础。【方法】以优良品种豫谷18和冀谷19为亲本构建的包含400个家系的RIL群体为试验材料,于2018—2019年分别在4个不同环境下调查谷子抽穗期、抽穗-成熟天数、全生育期及穗长、穗粗和单穗重等穗相关性状的表型值。同时,利用已构建的由1 304个bin标记组成的全长为2 196 cM,标记间平均距离为1.68 cM的高密度遗传连锁图谱。采用复合区间作图法（composite interval mapping,CIM）对生育期及穗部性状进行QTL定位分析,并对所获得的QTL置信区间进行候选基因的预测。【结果】重组自交系群体生育期及穗部性状在4个环境中均表现为连续分布且存在双向超亲分离现象,符合数量性状的遗传特征,适宜进行QTL分析。相关分析表明,谷子抽穗期与全生育期呈极显著正相关,与抽穗-成熟天数呈显著负相关,穗长与穗粗呈显著正相关。共检测到88个与谷子生育期及穗部性状相关的QTL,分布在第1、3、5、6、8和9染色体上。其中45个QTL与抽穗期相关,单个QTL能够解释表型变异的1.61%—27.60%;7个QTL与抽穗-成熟天数相关,单个QTL能够解释表型变异的2.65%—12.14%;20个QTL与全生育期相关,单个QTL能够解释表型变异的1.98%—16.97%;9个QTL与穗长相关,单个QTL能够解释表型变异的3.51%—11.65%;5个QTL与穗粗相关,单个QTL能够解释表型变异的3.74%—8.34%;2个QTL与单穗重相关,单个QTL能够解释表型变异的5.16%—5.20%。本研究共检测到12个主效QTL,其中,qEHD-9-1、qEHD-9-2、qHMD-9-2、qGRP-9-2和qPL-5-1在至少2个环境和BLUP值中被重复检测到。控制生育期的主效QTL（qEHD-9-1、qHMD-9-1、qGRP-9-1）与控制穗长的主效QTL（qPL-9-1）在第9染色体重叠;qEHD-9-2、qHMD-9-3、qGRP-9-2和qPL-9-3也在第9染色体重叠;控制穗长的主效QTL（qPL-5-1）和控制穗粗的QTL（qPD-5-1）在第5染色体重叠。对3个QTL簇的置信区间进行基因注释,筛选出5个与生育期及穗部性状相关的候选基因,其中,2个候选基因在谷子生育期调控和穗部性状发育中均发挥重要作用。【结论】共检测到88个与谷子生育期及穗部性状相关的QTL,12个为主效QTL,其中5个主效QTL在多个环境被重复检测到,且成簇分布。基于基因注释,共筛选了5个与谷子生育期和穗部性状相关的候选基因。     

甘蓝型油菜角果长度全基因组关联分析

【目的】挖掘与油菜角果长度性状显著相关的SNP位点及候选基因,为揭示油菜角果长度性状的遗传基础和分子机制提供理论依据,为油菜产量分子标记辅助选择育种奠定基础。【方法】在江西农业大学试验地和江西省红壤研究所试验地2个环境下考察300份甘蓝型油菜自交系的角果长度性状,利用简化基因组测序技术(specific locus amplified fragment sequencing,SLAF-seq)对300份甘蓝型油菜自交系基因组DNA进行测序并分析,利用获得的均匀分布于甘蓝型油菜基因组上的201 817个群体SNP(single nucleotide polymorphism,SNP)对角果长度性状进行全基因组关联分析(genome-wide association study,GWAS),探测与油菜角果长度显著相关的SNP位点,并基于群体连锁不平衡分析结果搜寻显著SNP位点两侧100 kb范围内的基因,通过BLAST获得关联区域内基因的注释信息,根据注释信息找出与性状相关的候选基因。【结果】农大试验地角果长度表型变异幅度为46.35—107.07 mm;红壤所试验地角果长度表型变异幅度为39.41—101.35 mm,两性状在2个环境下均表现出广泛表型变异。通过一般线性模型(general linear model,GLM)关联分析,农大环境下共检测到121个角果长度显著关联的SNP位点,分布在A04、A06、A08、A09、C02、C03、C06和C09等8条染色体上,其中,A09染色体上分布最多(83个SNP),红壤所环境下检测到22个角果长度显著关联的SNP位点,其中,1个在C09染色体上,其余21个均分布于A09染色体,在两地探测到20个一致性SNP位点;通过混合线性模型(mixed linear model,MLM)分析,农大环境下共检测到5个角果长度显著关联的SNP位点,其中,3个SNP位点与红壤所环境下检测到3个SNP位点一致,所有位点均位于A09染色体上。对MLM关联分析得到的显著SNP位点两侧100 kb区域内基因进行搜寻并进行功能注释,发现多个候选基因参与调节碳水化合物的运输与合成、花器官和种子的发育、信号转导等,它们可能通过上述功能影响油菜角果的生长,导致角果长度的差异。【结论】通过GLM和MLM两种分析方法探测到多个与油菜角果长度性状显著关联的基因位点,并在显著性位点附近搜寻到相关候选基因。     

葡萄果粒质量相关性状全基因组关联分析

【目的】果粒大小是葡萄外观和产量的重要构成因子之一,为受多基因调控的复杂数量性状,挖掘葡萄果粒大小相关性状的关键遗传调控位点和基因,将有助于葡萄产量的提高。【方法】本研究以150份葡萄品种资源为材料,分别于2019年和2020年对葡萄果实单粒重、种子数目和种子质量等进行测定,并结合重测序获得的高密度基因型数据进行全基因组关联分析（genome-wide association study,GWAS）,挖掘调控各性状的遗传位点和基因。【结果】各性状在关联群体中呈现广泛的连续变异,变异系数为39.55%—68.89%;在不同年份均服从正态分布,符合数量性状遗传特征;相关性分析表明葡萄果实单粒重、种子数目和种子质量呈显著正相关。全基因组关联分析共检测到150个与果实单粒重显著关联的SNP,在2019年检测到99个SNP,解释表型变异的14.48%—25.59%;在2020年检测到73个SNP,解释表型变异的16.08%—26.83%;其中24个SNP位点在两个年份均检测到,主要位于1号、5号、11号和16号染色体。相较于果实单粒重,检测到的与种子数目显著关联的SNP较少,2019年检测到1个...     

新疆冬小麦品种资源主要产量性状全基因组关联分析

【目的】高产是小麦育种的永恒主题,利用全基因组关联分析发掘控制小麦产量性状的QTL区段及优异基因,为小麦分子标记辅助选择育种提供理论依据和标记信息。【方法】以新疆本地188个冬小麦品种资源为材料,利用小麦55K SNP芯片进行全基因组扫描,通过对6个不同环境下的株高、穗长、小穗数、可育小穗数、穗粒数、千粒重、粒长、粒宽、籽粒长/宽比9个产量相关性状进行表型鉴定,利用6个环境下各性状数据及最佳线性无偏预测（BLUP）数据,基于混合线性模型（MLM）对表型和基因型进行全基因组关联分析。【结果】经主成分分析,将188个材料分为地方品种和育成品种2个亚群;利用6个环境下各性状数据,9个性状共检测到1 309个显著性SNP标记,其中,每个显著性SNP位点可解释7.259%—70.792%的表型变异。利用BLUP数据,9个性状共检测到66个显著性位点,同时与2个性状关联的共有SNP位点有5个,贡献率波动范围为8.498%—21.877%。将同时与2个性状或2个以上环境关联到的重复位点作为稳定的显著性关联位点,9个性状共检测到38个稳定关联位点,包括株高重复位点5个,穗长重复位点10个,小穗数重复位...     

RTM-GWAS方法应用于大豆RIL群体百粒重QTL检测的功效

【目的】为全面解析大豆重组自交系群体中调控百粒重性状的QTL体系,将限制性两阶段多位点全基因组关联分析方法（RTM-GWAS）和不同定位方法进行比较、优选,为后续候选基因体系探索及分子标记辅助育种设计提供依据。【方法】利用以科丰1号和南农1138-2为亲本衍生的重组自交系群体NJRIKY的427个家系,通过由全基因组39 353个SNP构建的3 683个SNPLDB标记及3个环境下的百粒重表型数据,选用复合区间作图法（CIM）、基于混合线性模型的全基因组关联分析方法（MLM-GWAS）和RTM-GWAS3种方法检测百粒重QTL,通过QTL数目和总的表型变异解释率比较检测功效,挑选最佳定位结果进行NJRIKY群体中的百粒重遗传体系解析。通过候选基因体系的功能注释,挖掘调控大豆百粒重的生物学途径。【结果】科丰1号与南农1138-2的百粒重差异较大,多环境平均数分别为9.0和17.9 g,遗传变异系数为12.4%,遗传率为85.4%,适用于百粒重性状的遗传解析。比较3种方法定位结果表明RTM-GWAS方法表现最佳,检测QTL数目最多（57个）,解释表型变异最多（70.78%）。而CIM仅检测到14个QTL,解释了56.47%的表型变异,MLM-GWAS仅定位到6个QTL,解释了18.47%的表型变异。RTM-GWAS共检测到57个QTL,分布在19条染色体上,表型变异解释率为0.03%—7.57%,其中41个QTL覆盖了已报道的来自30个双亲群体的81个百粒重QTL,16个QTL为新发现位点,包含一个表型变异解释率大于3%的大效应位点Sw-09-2。此外,检测的57个QTL中有20个位点与环境存在互作效应。这57个QTL构成了影响NJRIKY群体百粒重性状的遗传体系。通过SNPLDB标记与预测基因内的SNP进行χ²检验,共筛选到36个候选基因,其中4个候选基因来自大效应QTL,剩余32个候选基因来自小效应QTL。通过GO注释发现这些候选基因功能注释丰富,其中13个候选基因与籽粒发育直接相关,剩余的候选基因功能丰富,包含转运、转录调节因子等,表明不同生物学途径的基因共同调控NJRIKY群体中百粒重性状的表达。【结论】3种定位方法中,高效的RTM-GWAS方法检测到较为全面的NJRIKY群体的百粒重QTL,更适用于双亲RIL群体的QTL定位。不同功能的候选基因共同调控了复杂的百粒重性状的表达。     

RTM-GWAS方法应用于大豆RIL群体百粒重QTL检测的功效

【目的】为全面解析大豆重组自交系群体中调控百粒重性状的QTL体系,将限制性两阶段多位点全基因组关联分析方法(RTM-GWAS)和不同定位方法进行比较、优选,为后续候选基因体系探索及分子标记辅助育种设计提供依据。【方法】利用以科丰1号和南农1138-2为亲本衍生的重组自交系群体NJRIKY的427个家系,通过由全基因组39 353个SNP构建的3 683个SNPLDB标记及3个环境下的百粒重表型数据,选用复合区间作图法(CIM)、基于混合线性模型的全基因组关联分析方法(MLM-GWAS)和RTM-GWAS 3种方法检测百粒重QTL,通过QTL数目和总的表型变异解释率比较检测功效,挑选最佳定位结果进行NJRIKY群体中的百粒重遗传体系解析。通过候选基因体系的功能注释,挖掘调控大豆百粒重的生物学途径。【结果】科丰1号与南农1138-2的百粒重差异较大,多环境平均数分别为9.0和17.9 g,遗传变异系数为12.4%,遗传率为85.4%,适用于百粒重性状的遗传解析。比较3种方法定位结果表明RTM-GWAS方法表现最佳,检测QTL数目最多(57个),解释表型变异最多(70.78%)。而CIM仅检测到14个QTL,解释了56.47%的表型变异,MLM-GWAS仅定位到6个QTL,解释了18.47%的表型变异。RTM-GWAS共检测到57个QTL,分布在19条染色体上,表型变异解释率为0.03%—7.57%,其中41个QTL覆盖了已报道的来自30个双亲群体的81个百粒重QTL,16个QTL为新发现位点,包含一个表型变异解释率大于3%的大效应位点Sw-09-2。此外,检测的57个QTL中有20个位点与环境存在互作效应。这57个QTL构成了影响NJRIKY群体百粒重性状的遗传体系。通过SNPLDB标记与预测基因内的SNP进行χ2检验,共筛选到36个候选基因,其中4个候选基因来自大效应QTL,剩余32个候选基因来自小效应QTL。通过GO注释发现这些候选基因功能注释丰富,其中13个候选基因与籽粒发育直接相关,剩余的候选基因功能丰富,包含转运、转录调节因子等,表明不同生物学途径的基因共同调控NJRIKY群体中百粒重性状的表达。【结论】3种定位方法中,高效的RTM-GWAS方法检测到较为全面的NJRIKY群体的百粒重QTL,更适用于双亲RIL群体的QTL定位。不同功能的候选基因共同调控了复杂的百粒重性状的表达。     

盐胁迫下甘蓝型油菜发芽期下胚轴和根长的全基因组关联分析

【目的】解析甘蓝型油菜发芽期根和下胚轴发育及耐盐性的调控位点,筛选油菜耐盐性相关的候选基因,可为油菜耐盐性改良提供依据。【方法】以317份具有代表性的甘蓝型油菜自交系为材料,在正常生长和盐胁迫条件下进行沙培鉴定,利用芸薹属60K SNP芯片和全基因组关联分析鉴定正常生长与盐胁迫下甘蓝型油菜发芽期根和下胚轴长度显著关联的SNP,并确定其连锁不平衡区间。通过区间内基因的功能注释及盐胁迫下油菜幼苗根和叶片转录组差异表达基因筛选连锁不平衡区间内的重要候选基因,并以实时荧光定量PCR分析候选基因的组织特异性和盐胁迫诱导表达模式,提高候选基因筛选的准确性。【结果】正常生长和盐胁迫下甘蓝型油菜发芽期下胚轴和根长在不同材料间变异较大,频次分布表明目标性状均为数量性状,受多基因调控。全基因组关联分析模型比较表明,MLM+P+K模型为最优模型。以此模型对目标性状进行全基因组关联分析,检测到45个显著关联SNP,其中40个与下胚轴长度显著关联,5个与根长显著关联,单个SNP解释的表型变异分别为9.12%—14.46%和7.67%—8.93%。重复检测的显著相关SNP中,值得注意的是C04染色体的rs8970,同时与4个性状显著关联,表型贡献率为7.67%—12.35%,是唯一在下胚轴长和根长间重复检测到的显著关联SNP。11个重要关联SNP中有6个位于10—442 kb的连锁不平衡区块中。转录组分析表明,11个连锁不平衡区间共包含447个基因,其中15个受盐胁迫诱导表达。转录组和基因功能注释综合分析表明,BnaSRO1、BnaPAGR2、BnaNPH3、BnaMYB124、BnaSAM-Mtase、BnaBIN2、BnaUMAMIT11、BnaEXPA7、BnaRPT3、BnaEF-hand和BnaF3H很可能为各自区间的候选基因。实时荧光定量PCR结果证实除BnaNPH3外,其他基因均在根或下胚轴中受盐胁迫诱导上调表达。组织特异性分析还发现BnaUMAMIT11、BnaPAGR2和BnaEXPA7主要在萌发的根和下胚轴中特异表达,BnaRPT3、BnaBIN2和BnaMYB124虽然呈组成型表达,但在萌发阶段的下胚轴中表达量最高,证实这些基因很可能参与油菜发芽期根和下胚轴生长发育及耐盐性的调节。【结论】全基因组关联分析共鉴定出45个控制油菜发芽期根和下胚轴发育及耐盐性的显著关联SNP。连锁不平衡、转录组和基因功能注释综合分析初步鉴定出11个重要候选基因。     

花生每荚种子数相关性状QTL的定位

【目的】花生是重要的油料作物和经济作物,高产一直是花生育种的主要目标,决定产量的因素是单位面积的种子数和仁重。单位面积种子数是种植密度、每株荚数和每荚种子数的乘积。因此,对花生每荚果种子数相关性状进行QTL分析,有助于发掘该性状相关基因/位点,为花生产量相关性状分子育种提供重要的理论依据。【方法】以四粒红×冀农黑3号构建的RIL群体为研究材料,于2018年（E1）和2020年（E2）在河北省保定市河北农业大学清苑试验站（115°30′E,38°40′N）种植鉴定,收获时调查统计单仁果数、双仁果数以及多仁果数表型值,利用河北农业大学花生创新团队实验室构建的高密度遗传图谱,采用QTL Icimapping V4.2中的完备区间作图法对2个环境下的每荚种子数相关性状进行QTL定位与分析。【结果】单仁果率与双仁果率均呈正态分布,多仁果率呈偏正态分布。3个性状的QTL定位分析结果表明,共检测到11个QTL,可解释4.66%—22.34%的表型变异,加性效应为-9.35—9.42。其中,定位到5个多仁果率QTL,可解释3.19%—22.34%的表型变异,有1个QTL的加性效应为负值（-4.77）,来自冀农黑3号,其余4个QTL的加性效应为正值（3.59—9.42）,均来自母本四粒红;定位到2个单仁果率QTL,可解释4.97%—6.43%的表型变异,加性效应均为负值（-4.45和-4.54）,均来自冀农黑3号;定位到4个双仁果率QTL,可解释3.46%—20.87%的表型变异,加性效应均为负值（-9.35—-3.84）,均来自冀农黑3号。这些QTL中,6个为主效QTL,其中,qRMSPA05被重复检测到,且可遗传表型变异为16.58%—17.34%,加性效应为7.69—8.12。【结论】定位6个主效QTL和1个主效稳定的多仁果率QTL,有助于改良花生产量性状,可以作为遗传改良的重要候选区段,用于分子标记辅助选择与精细定位研究。     

玉米株型相关性状的全基因组关联分析

【目的】玉米株型性状与植株产量、光合效率、抗倒性等密切相关,是理想株型设计育种的基础,通过对玉米多个株型相关性状进行全基因组关联分析,构建玉米理想株型,为抗倒伏、耐密性的玉米新品种选育提供理论基础。【方法】以284份温带、亚热带和热带材料组成的关联群体为研究对象,在郑州与浚县2个环境下对玉米总叶片数（LN）、穗上叶片数（LNAN）、穗上叶片数与总叶片数比值（LNAN/LN）、株高（PH）、穗位高（EH）、穗位高与株高的比值（EH/PH）等6个玉米株型相关性状进行调查,借助覆盖玉米全基因组约56万个SNP标记,进行全基因组关联分析,以期为玉米新品种的选育提供理论依据。【结果】2个环境下,6个性状均表现为正态分布;大多数性状之间均存在高度正相关或负相关;方差分析表明这6个株型相关性状的基因型与环境以及基因型×环境的互作均达到显著水平。在选择最优关联分析模型时,发现Q模型假阳性较高,Q+K模型对假阳性的控制过于严格,而K模型的表现最好;在2个环境下,全基因组关联分析（K模型）共检测到56个显著SNP-性状关联（P≤3.99E-6）,这些SNPs共涉及5个性状的17个位点,每个位点解释的表型变异从7.97%-10.56%不等。同时发现有4个位点能够在2个环境中同时被检测到,表明这4个位点受环境影响较小,在不同环境下可以稳定存在。通过对显著关联的SNP上下游各50 kb范围内候选基因进行搜索,共发现80个候选基因,其中42个具有功能注释。例如,与株高和穗位高显著相关的基因GRMZM2G161293编码一个具有乙酰葡糖转移酶活性的蛋白,该酶催化UDP-N-乙酰氨基葡萄糖生成糖过程中的N-乙酰葡糖氨基残基的转移,可能通过影响玉米籽粒可溶性糖含量进而影响产量,推测其为最可能的候选基因。【结论】K模型对假阳性的控制效果最好,基于K模型的GWAS结果,一共检测到17个株型性状相关的位点。     

籼稻骨干亲本产量相关性状遗传效应

【目的】研究水稻产量相关性状的遗传效应、遗传率和遗传相关,旨在加深理解水稻产量相关性状遗传体系。【方法】以几个重要的恢复系和不育系为材料,采用加性-显性-母体效应遗传模型(ADM模型)和统计分析方法研究水稻产量相关性状遗传效应、遗传率和遗传相关。【结果】研究表明水稻产量相关性状同时受加性效应、显性效应和母体效应的控制,穗总粒数、单株总粒数、千粒重和单株重以加性效应占主导,株高和单株有效穗则由显性效应和母体效应共同控制;狭义遗传率以千粒重的表现最高,单株重次之;不育系中粤泰A能增加单株有效穗,Y58S和中9A能增加单株总粒数;恢复系中绵恢725能增加穗总粒数、单株粒数、单株重和千粒重;此外,遗传相关研究表明通过增加穗总粒数和单株总粒数可以达到增加单株重的目的,而要选到穗总粒数较多的材料有效穗不能太多。【结论】在水稻产量相关性状育种实践中,应根据水稻产量相关性状的遗传方式不同采用不同的选育方法。穗总粒数、单株总粒数、千粒重、单株重和播始历期可在常规聚合育种上通过世代综合选择,使加性效应得以稳定遗传,而株高和单株有效穗在杂种优势利用上有一定的潜力。千粒重和单株重早期世代选择效果较好,株高和单株有效穗高代选择比较适宜。     

小麦苗期根系性状的全基因组关联分析与优异位点挖掘

【目的】植物根系对水分及营养的获取、作物的生长发育和产量的形成至关重要。挖掘小麦苗期根系性状显著关联的SNP位点,预测相关候选基因,为解析小麦根系建成遗传机制及选育具有优良根系构型的小麦品种奠定基础。【方法】以189份小麦品种组成的自然群体为供试材料,调查2种培养条件（霍格兰营养液和去离子水）下培育21 d的苗期根系总长度（TRL）、根系总表面积（TRA）、根系总体积（TRV）、根系平均直径（ARD）及根系干重（RDW）等5个根系性状,试验进行2次重复,同时结合小麦660K SNP芯片的分型结果进行全基因组关联分析（genome-wide association study,GWAS）。此外,通过序列比对、结构域分析和注释信息预测候选基因,并采用竞争性等位基因特异性PCR(kompetitive allele specific PCR,KASP)技术开发根系性状的分子标记。【结果】霍格兰营养液培养条件下的根系性状变异范围较大,根系整体粗短;而去离子水条件下的根系细长、侧根较多。选用贝叶斯信息与连锁不平衡迭代嵌套式模型（BLINK）、压缩式混合线性模型（CMLM）、固定随机循环概率模型（...     

陆海BC4F2和BC4F3代换系的评价及纤维产量与品质相关QTL的检测

 【目的】利用SSR标记对陆海BC₄F₂和BC₄F₃代换系进行评价并检测纤维产量与品质相关的QTL,为筛选棉花染色体单片段代换系、精细定位纤维品质QTL、实现分子聚合育种奠定基础。【方法】利用GGT32（graphical genotyping）软件分析每个代换系的基因型组成,采用SAS PROC GLM的单向方差分析方法检测影响各性状的QTL。【结果】检测到50个单片段代换系,其中9株含有纯合的海岛棉片段,并筛选出12个代换片段少、纤维品质优良的代换系。共检测到15个控制产量性状和19个控制纤维品质的QTL,集中分布在12个连锁群中,解释的表型变异率在2.80%—14.13%。【结论】4个上半部平均长度QTL在2个世代中稳定遗传,1个上半部平均长度QTL在前人研究论文中检测到,部分标记位点同时控制几个不同的性状,并发现增效基因不全来自高值亲本。