首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到17条相似文献,搜索用时 605 毫秒
1.
油菜耐盐相关性状的全基因组关联分析及其候选基因预测   总被引:2,自引:1,他引:1  
【目的】通过对甘蓝型油菜耐盐相关性状进行全基因组关联分析,寻找可能与油菜耐盐性相关的SNP位点,发掘与油菜耐盐性有关的候选基因。【方法】以1.2%NaCl溶液作为培养液,去离子水为对照,对307个不同品系的甘蓝型油菜进行发芽试验。播种后7 d测定幼苗根长、鲜重及发芽率,计算盐胁迫下各性状相对值,并作为评价耐盐的指标。结合油菜60K SNP芯片,利用SPAGeDi v1.4软件对该群体307份甘蓝型油菜进行亲缘关系分析,并计算亲缘关系值的矩阵。利用软件STRUCTURE v2.3.4对关联群体进行了群体结构分析。为有效排除假关联的影响,将群体结构和材料间的亲缘关系考虑进关联分析中,同时进行了PCA+K模型、Q+K模型以及K模型3种混合线性模型分析和比较,根据所有SNP的–lg(P)观察值和期望值,确定每个性状GWAS分析的最优模型。采用TASSEL 5.0软件,在最优模型下对307份材料耐盐各性状的相对值分别与SNP标记进行全基因组关联分析。利用油菜基因组数据库,在显著SNP位点侧翼序列200 kb范围内提取基因。根据拟南芥中已经明确功能的耐盐相关基因,筛选出目标基因组区段内与耐盐相关的油菜同源基因。【结果】全基因组关联分析共检测到164个与根长显著关联的SNP位点,23个与鲜重显著关联的SNP位点,38个与发芽率显著关联的SNP位点。其中与根长、鲜重、发芽率最显著关联的SNP位点分别位于染色体A08、A02和A06,贡献率分别为23.84%、18.59%和31.81%。在这些显著SNP位点侧翼序列200 kb范围内发掘出可能与油菜耐盐性有关的50个候选基因。这些候选基因主要包括转录因子MYB、WRKY、ABI1、b ZIP、ERF1、CZF1、XERICO等以及一些下游受转录因子调控的不同功能基因NHX1、PTR3、CAT1、HKT、CAX1、ACER、STH、STO等。在根长和发芽率2个不同耐盐性状的分析结果中均筛选出位于A03染色体上的耐盐基因BnaA03g14410D。另外,这些耐盐候选基因中包含两组串联重复基因,分别是位于A03染色体上的BnaA03g18900D和BnaA03g18910D,位于C09染色体上的BnaC09g19080D、BnaC09g19090D和BnaC09g19100D。除此之外,耐盐候选基因中还包含2个距离非常近(中间只间隔2个基因)的重复基因BnaC02g39600D和BnaC02g39630D。【结论】共检测到225个与耐盐性状显著关联的SNP位点,筛选出50个可能与油菜耐盐性有关的候选基因。  相似文献   

2.
苯磺隆胁迫下油菜萌发期相关性状的全基因组关联分析   总被引:2,自引:2,他引:0  
目的 研究苯磺隆残留对油菜种子萌发的影响,运用全基因组关联分析(genome-wide association analysis,GWAS)揭示苯磺隆胁迫下油菜萌发期相关性状的遗传因子和候选基因,探究油菜在苯磺隆逆境胁迫下的生理形态所反映的基因调控机制,为耐苯磺隆油菜品种的研究提供参考。方法 以241份甘蓝型油菜品种(系)为材料、25 mg·L -1苯磺隆溶液为处理液、蒸馏水为对照进行发芽试验。发芽7 d测定并计算相对发芽率、相对根长和相对鲜重。结合芸薹属60K SNP芯片分析群体基因型,通过STRUCTURE软件和TASSEL软件分别对该群体进行群体结构分析以及亲缘关系和LD衰减分析。为有效排除假关联的影响,采用一般线性模型(GLM)和混合线性模型(MLM)中的6种模型进行比较,确定每个性状GWAS分析的最优模型。同时,利用TASSEL软件在最优模型下对241份材料的3个性状分别进行全基因组关联分析,根据关联SNP位点的LD区间序列预测候选基因。 结果 241份品种(系)群体可分为P1(94份材料)和P2(147份材料)2个亚群,其中约56.28%的材料之间的亲缘关系值为0。全基因组关联分析(K+PCA模型)共检测到16个与性状显著关联的SNP位点,这些位点可解释9.42%—13.14%的表型变异率。通过分析显著SNP位点的LD区间与甘蓝型油菜对应的区间序列,筛选出25个候选基因可能与油菜耐苯磺隆有关,其中9个为细胞色素P450家族基因,5个参与谷胱甘肽合成或代谢过程,2个为多药耐药相关蛋白基因。同时发现与相对发芽率显著相关的基因ATGSTU19编码谷胱甘肽转移酶,参与毒素分解过程,在各种胁迫反应中起重要作用。在相对根长和相对鲜重共同鉴定到的候选基因BnaC02g27690D功能未知。结论 共检测到16个SNP位点与耐苯磺隆性状显著关联,筛选出25个候选基因可能与油菜耐苯磺隆有关。  相似文献   

3.
芝麻发芽期耐盐性鉴定方法研究及耐盐候选基因的挖掘   总被引:1,自引:0,他引:1  
【目的】确定芝麻发芽期耐盐性鉴定适宜的NaCl胁迫浓度和评价指标,发掘耐盐种质和耐盐相关重要基因,为芝麻耐盐性大规模鉴定、遗传改良和耐盐机理研究提供方法借鉴和优异基因资源。【方法】以8份耐盐性差异较大的芝麻种质为材料,在不同浓度的NaCl(0、50、100、150、200和250 mmol·L-1)胁迫下发芽,测定其发芽势、成苗率、根长、芽长和鲜重等指标,通过对指标值的方差分析、主成分分析、隶属函数和相关性分析等,筛选芝麻发芽期耐盐性鉴定适宜的NaCl处理浓度和评价指标。以100 mmol·L-1 NaCl溶液为处理浓度,相对成苗率为鉴定指标对71份芝麻核心种质进行发芽期耐盐性鉴定和全基因组关联分析,并对获得的候选基因进行功能注释;通过NaCl胁迫下芝麻幼苗叶片转录组测序和荧光定量PCR分析候选基因的表达模式,筛选耐盐相关候选基因。【结果】对不同浓度NaCl胁迫下8份芝麻材料发芽各指标进行统计和方差分析表明,100 mmol·L-1的NaCl处理下,除发芽势外,发芽指数、活力指数、成苗率、根长、芽长和苗鲜重的标准偏差值均较大,所有指标在0.05水平上均存在显著差异,100 mmol·L-1的 NaCl溶液可以作为芝麻发芽期耐盐性鉴定的适宜胁迫浓度;相对成苗率、相对发芽势、相对发芽指数、相对活力指数、相对芽长、相对根长和相对鲜重7个指标与芝麻发芽期耐盐性关系密切,贡献率较高,可以作为芝麻发芽期耐盐性鉴定的评价指标。71份芝麻核心种质材料的相对成苗率变异比较丰富,符合正态分布;通过对71份芝麻核心种质相对成苗率与SNP标记进行全基因组关联分析,检测到7个与芝麻发芽期耐盐性显著关联的SNP标记(LG5:688003、LG7:9582027、LG10:5274091、LG10:10788493、LG11:11924186、LG14:2128695和LG16:3930301),SNP标记上、下游各100 kb区间内共有基因67个,其中有功能注释的基因34个;通过耐盐材料S04在盐胁迫下的转录组测序和荧光定量PCR对预测的候选基因进行表达模式分析,共有21个候选基因受到盐胁迫诱导显著差异表达。【结论】芝麻发芽期耐盐性鉴定的适宜NaCl胁迫浓度为100 mmol·L-1,相对成苗率等7个指标可以作为适宜的评价指标;检测到与芝麻发芽期耐盐性显著关联的SNP标记7个,并鉴定出耐盐候选基因21个。  相似文献   

4.
甘蓝型油菜角果长度全基因组关联分析   总被引:2,自引:0,他引:2  
【目的】挖掘与油菜角果长度性状显著相关的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两种分析方法探测到多个与油菜角果长度性状显著关联的基因位点,并在显著性位点附近搜寻到相关候选基因。  相似文献   

5.
【目的】近几年随着观光农业的兴起,花色的选育和改良已成为甘蓝型油菜种质资源鉴定和材料创制的重要研究方向。以甘蓝型油菜黄白花分离F2群体为研究对象,通过二代测序技术,对白花性状基因候选区间定位,开发与白花性状连锁的分子标记,为定位白花候选基因和选育白花新材料提供新思路。【方法】以甘蓝型油菜DH纯系黄花Y05和甘蓝型油菜纯系白花W01杂交,观察F1和F2群体的花色分离,分析白花性状遗传模式。在F2群体中选取30株纯白花和30株纯黄花构建DNA叶片子代池和RNA花瓣子代池,对亲本和DNA叶片子代池进行30×重测序,对RNA花瓣子代池进行5×测序。以法国甘蓝型油菜Darmor-bzh、中双11、Darmor、Tapidor为参考序列,重测序QTL-seq分析流程计算2个DNA子代池的SNP-index和delta(SNP-index)。利用R包画出SNP-index和delta(SNP-index)滑窗分析图,鉴定候选区间。转录组MMAPPR分析流程以法国甘蓝型油菜Darmor-bzh为参考序列,计算SNP频率,ED4(Loess fit)检测峰值和鉴定候选区间。利用MISA进行重复序列鉴定,使用Prime3在候选区间进行SSR引物设计,在F2群体中采用聚丙烯酰胺凝胶电泳方法对SSR引物进行筛选。【结果】甘蓝型油菜黄花与白花杂交F2群体中,白花和黄花性状分离比符合3﹕1,暗示白花性状受1对显性主效基因控制。全基因组重测序区间定位结果显示,白花性状基因候选区间在Darmor-bzh C03染色体52—55 Mb。同时以甘蓝型油菜中双11、Darmor、Tapidor分别为参考序列,均鉴定出白花基因候选区间在C03染色体上的一致性和稳定性。转录组测序定位白花性状基因位于Darmor-bzh C03染色体54—55 Mb。转录组测序和重测序定位染色体结果高度一致。在此区间内MISA和Primer3结合设计SSR引物,聚丙烯酰胺凝胶电泳筛选到6个与白花性状紧密连锁共分离的SSR标记。6个SSR标记区间范围在760 kb(52.81—53.57 Mb)。此候选区间与甘蓝、白菜共线性分析,对应白菜A02染色体56.76—57.40 Mb区间,对应甘蓝C03染色体10.99—11.28 Mb区间。【结论】甘蓝型油菜白花性状由1对显性主效基因控制。白花性状基因候选区间在法国甘蓝型油菜Darmor-bzh C03染色体52—55 Mb区间内。此区间760 kb范围内筛选出6个与白花性状基因紧密连锁共分离的SSR标记。  相似文献   

6.
【目的】近几年随着观光农业的兴起,花色的选育和改良已成为甘蓝型油菜种质资源鉴定和材料创制的重要研究方向。以甘蓝型油菜黄白花分离F2群体为研究对象,通过二代测序技术,对白花性状基因候选区间定位,开发与白花性状连锁的分子标记,为定位白花候选基因和选育白花新材料提供新思路。【方法】以甘蓝型油菜DH纯系黄花Y05和甘蓝型油菜纯系白花W01杂交,观察F1和F2群体的花色分离,分析白花性状遗传模式。在F2群体中选取30株纯白花和30株纯黄花构建DNA叶片子代池和RNA花瓣子代池,对亲本和DNA叶片子代池进行30×重测序,对RNA花瓣子代池进行5×测序。以法国甘蓝型油菜Darmor-bzh、中双11、Darmor、Tapidor为参考序列,重测序QTL-seq分析流程计算2个DNA子代池的SNP-index和delta(SNP-index)。利用R包画出SNP-index和delta(SNP-index)滑窗分析图,鉴定候选区间。转录组MMAPPR分析流程以法国甘蓝型油菜Darmor-bzh为参考序列,计算SNP频率,ED 4(Loess fit)检测峰值和鉴定候选区间。利用MISA进行重复序列鉴定,使用Prime3在候选区间进行SSR引物设计,在F2群体中采用聚丙烯酰胺凝胶电泳方法对SSR引物进行筛选。【结果】甘蓝型油菜黄花与白花杂交F2群体中,白花和黄花性状分离比符合3﹕1,暗示白花性状受1对显性主效基因控制。全基因组重测序区间定位结果显示,白花性状基因候选区间在Darmor-bzh C03染色体52—55 Mb。同时以甘蓝型油菜中双11、Darmor、Tapidor分别为参考序列,均鉴定出白花基因候选区间在C03染色体上的一致性和稳定性。转录组测序定位白花性状基因位于Darmor-bzh C03染色体54—55 Mb。转录组测序和重测序定位染色体结果高度一致。在此区间内MISA和Primer3结合设计SSR引物,聚丙烯酰胺凝胶电泳筛选到6个与白花性状紧密连锁共分离的SSR标记。6个SSR标记区间范围在760 kb(52.81—53.57 Mb)。此候选区间与甘蓝、白菜共线性分析,对应白菜A02染色体56.76—57.40 Mb区间,对应甘蓝C03染色体10.99—11.28 Mb区间。【结论】甘蓝型油菜白花性状由1对显性主效基因控制。白花性状基因候选区间在法国甘蓝型油菜Darmor-bzh C03染色体52—55 Mb区间内。此区间760 kb范围内筛选出6个与白花性状基因紧密连锁共分离的SSR标记。  相似文献   

7.
8.
为解析甘蓝型油菜幼苗期下胚轴长、根长与成熟期主要农艺性状及品质性状的关系,以125份甘蓝型油菜品种(系)为材料,进行了室内和田间试验,并对幼苗期下胚轴长、根长与成熟期主要农艺性状及品质性状进行了相关性分析及灰色关联度分析。结果表明,甘蓝型油菜芥酸含量、硫甙含量变异系数较大,易受外界环境影响,而含油量、株高等变异系数较小,遗传力较高。幼苗期下胚轴长与根长呈极显著正相关,且与幼苗期根长关联度最大的也是幼苗期下胚轴长,幼苗期根长与倒伏级呈显著负相关。幼苗期下胚轴长、根长与主要品质性状均呈正相关,其中下胚轴长与芥酸含量呈显著正相关。可见,苗期根系越发达,成熟期抗倒伏性和品质性状越优良,选育抗倒伏和品质优的甘蓝型油菜时可将幼苗期下胚轴长与根长作为重要的参考指标。  相似文献   

9.
【目的】寻找苯磺隆胁迫下油菜种子萌发性状相关的QTL及其耐性基因,为筛选与培育耐苯磺隆油菜种质以及探究油菜种子萌发过程中苯磺隆耐性分子机理奠定基础。【方法】用0.15 mg·kg -1苯磺隆溶液处理由人工合成甘蓝型油菜10D130和甘蓝型油菜常规品种ZS11构建的包含175个株系的高世代重组自交系(RIL)群体,进行种子发芽试验,以蒸馏水为对照,分别测定其相对发芽势、相对发芽率、相对根长和相对干重。然后,利用油菜6K SNP芯片对该RIL群体进行基因分型,通过JoinMap4.0软件构建高密度遗传连锁图谱。基于该遗传图谱,利用MapQTL软件多QTL作图法对4个性状的相对值进行QTL定位,根据各QTL置信区间查找甘蓝型油菜的基因序列,并依次与拟南芥基因组序列进行BLAST,筛选可能与耐苯磺隆胁迫相关的候选基因。【结果】频数分布表明4个相对性状的变异范围较大,且呈连续性分布,符合数量性状表现特征,适宜进行QTL遗传分析。相关分析表明,相对发芽率和相对发芽势呈极显著正相关,相关系数为0.587。构建的遗传图谱包含1 897个多态性SNP标记,覆盖甘蓝型油菜基因组3 214.19 cM,标记之间的平均图距为1.69 cM。利用此图谱共检测到22个相关QTL,表型贡献率变幅为6.4%—12.6%。其中,与相对发芽势、相对发芽率相关的QTL分别有6个和3个,与相对根长和相对干重有关的QTL分别为8个和5个。在A01染色体64.857 cM、55.935 cM和56.645 cM处检测到的相对发芽势与相对发芽率QTL的置信区间完全或者部分重叠。通过分析QTL置信区间上甘蓝型油菜对应的区间序列,筛选到30个可能与油菜耐苯磺隆有关的候选基因,其中包括18个细胞色素P450家族成员、5个糖基转移酶家族基因、1个GSTF相关基因、1个ABC转运蛋白相关基因和1个ALS基因,这些基因均与除草剂抗性机制有关,尤其ALS为磺酰脲类除草剂靶位点酶;另外筛选到1个BHLH和1个JAZ6基因,BHLH与JAZ蛋白可通过相互作用来防御胁迫;检测到1个LSU2蛋白相关基因和1个MATE家族成员,前者参与细胞氧化剂解毒及植物防御反应,后者参与类黄酮、生物碱、金属离子、其他多种代谢物的转运及有毒物质引起的植物胁迫响应。【结论】检测到与相关QTL共22个,筛选出可能与苯磺隆耐性有关的候选基因30个。这些基因通过加速毒性分子的转运与代谢从而响应有毒物质引起的胁迫反应,可能参与植物对苯磺隆的抗性调节与反应机制。  相似文献   

10.
【目的】寻找苯磺隆胁迫下油菜种子萌发性状相关的QTL及其耐性基因,为筛选与培育耐苯磺隆油菜种质以及探究油菜种子萌发过程中苯磺隆耐性分子机理奠定基础。【方法】用0.15 mg·kg~(-1)苯磺隆溶液处理由人工合成甘蓝型油菜10D130和甘蓝型油菜常规品种ZS11构建的包含175个株系的高世代重组自交系(RIL)群体,进行种子发芽试验,以蒸馏水为对照,分别测定其相对发芽势、相对发芽率、相对根长和相对干重。然后,利用油菜6K SNP芯片对该RIL群体进行基因分型,通过JoinMap4.0软件构建高密度遗传连锁图谱。基于该遗传图谱,利用MapQTL软件多QTL作图法对4个性状的相对值进行QTL定位,根据各QTL置信区间查找甘蓝型油菜的基因序列,并依次与拟南芥基因组序列进行BLAST,筛选可能与耐苯磺隆胁迫相关的候选基因。【结果】频数分布表明4个相对性状的变异范围较大,且呈连续性分布,符合数量性状表现特征,适宜进行QTL遗传分析。相关分析表明,相对发芽率和相对发芽势呈极显著正相关,相关系数为0.587。构建的遗传图谱包含1 897个多态性SNP标记,覆盖甘蓝型油菜基因组3 214.19 cM,标记之间的平均图距为1.69 cM。利用此图谱共检测到22个相关QTL,表型贡献率变幅为6.4%—12.6%。其中,与相对发芽势、相对发芽率相关的QTL分别有6个和3个,与相对根长和相对干重有关的QTL分别为8个和5个。在A01染色体64.857 cM、55.935 cM和56.645 cM处检测到的相对发芽势与相对发芽率QTL的置信区间完全或者部分重叠。通过分析QTL置信区间上甘蓝型油菜对应的区间序列,筛选到30个可能与油菜耐苯磺隆有关的候选基因,其中包括18个细胞色素P450家族成员、5个糖基转移酶家族基因、1个GSTF相关基因、1个ABC转运蛋白相关基因和1个ALS基因,这些基因均与除草剂抗性机制有关,尤其ALS为磺酰脲类除草剂靶位点酶;另外筛选到1个BHLH和1个JAZ6基因,BHLH与JAZ蛋白可通过相互作用来防御胁迫;检测到1个LSU2蛋白相关基因和1个MATE家族成员,前者参与细胞氧化剂解毒及植物防御反应,后者参与类黄酮、生物碱、金属离子、其他多种代谢物的转运及有毒物质引起的植物胁迫响应。【结论】检测到与相关QTL共22个,筛选出可能与苯磺隆耐性有关的候选基因30个。这些基因通过加速毒性分子的转运与代谢从而响应有毒物质引起的胁迫反应,可能参与植物对苯磺隆的抗性调节与反应机制。  相似文献   

11.
油菜主花序角果密度及其相关性状的全基因组关联分析   总被引:2,自引:2,他引:0  
【目的】油菜高产是育种工作的主要研究目标之一。角果密度、主花序有效角果数等性状与产量都有显著或极显著的正相关关系,是油菜高产育种考查的主要性状。为揭示油菜角果密度及其相关性状的遗传机理和分子机制奠定基础。【方法】以不同遗传背景和地理来源的213份甘蓝型油菜品种(系)构成的自然群体为研究对象,利用芸薹属60K Illumina Infinium SNP芯片对该群体进行基因型分型。分别于2015年和2016年在成熟期调查该群体主花序有效长和主花序有效角果数,计算主花序角果密度。利用Structure 2.3.4软件对该群体进行群体结构分析,Tassel 5.1.0软件分析亲缘关系和染色体连锁不平衡的衰减;然后基于最优模型对主花序角果密度及其相关性状进行全基因组关联分析(genome-wide association analysis,GWAS),依据关联SNP位点的LD区间序列,预测与性状相关的重要关联候选基因。【结果】群体结构分析显示,213份甘蓝型油菜分为P1和P2亚群,P1亚群包含50份材料(23.5%),P2亚群包含163份材料(76.5%),基本上和油菜的地理栽培属性一致;亲缘关系发现约89.74%材料之间的亲缘关系值小于0.2,其中约有59.91%材料的亲缘关系值为0。总体来看,整个自然群体材料之间的亲缘关系比较远。对A、C基因组进行连锁不平衡分析发现,A和C基因组的r2随着遗传距离的增加而下降,A基因组的衰减距离整体比C基因组的衰减距离小。GWAS分析两年数据共检测到17个SNP位点与主花序角果密度及其相关性状关联。其中与主花序角果密度和主花序有效长相关的SNP标记分别有7个和9个,并分别解释11.34%—15.96%和9.67%—13.10%的表型变异;与主花序有效角果数相关联的标记有1个,解释11.56%的表型变异。通过分析关联SNP位点的LD区间与甘蓝型油菜对应的区间序列,找到22个与主花序角果密度及其相关性状有关的候选基因,其中BnaA01g16940D、BnaC01g38800D和BnaA04g09170D等主要通过调控赤霉素和生长素等内源激素的合成和信号转导来控制主花序角果密度及其相关性状;BnaA01g16970D、BnaA03g29180D、BnaA03g29810D、BnaC01g39680D和 BnaC03g32770D通过对分生组织的调控来改变表型;BnaC09g18690D和 BnaC09g09210D等主要通过控制细胞分裂生长等过程改变表型。【结论】检测到17个SNP标记与油菜主花序角果密度、主花序有效长和主花序有效角果数关联,筛选出22个与主花序角果密度及其相关性状有关的候选基因。  相似文献   

12.
【目的】根系作为植株吸收水分和养分的重要器官,对玉米生长及产量的形成至关重要。研究玉米根系结构的遗传机制指导玉米高产育种实践。【方法】以111份玉米优异自交系为材料,于2017年在北京、陕西永寿、山西定襄和河南原阳4个环境下对玉米地下节根层数(RLN)、地下节根总条数(TRN)、地下节根角度(RA)、地下节根面积(RS)、地下节根体积(RV)和地下节根干重(RDW)等6个玉米根系相关性状进行调查。取4个环境的平均值作为6个根系相关性状的表型数据,对6个相关性状进行统计分析和相关性分析,对不同年代、不同类群自交系的地下节根相关性状进行差异分析。基于该群体全基因组152 352个高质量SNP标记,利用FarmCPU模型进行全基因组关联分析获得显著关联SNP位点,并在LD衰减距离范围内查找候选基因,对候选基因的功能进行富集分析。【结果】表型分析表明,6个地下节根性状均呈现正态分布,且均显示出较高的遗传力;相关性分析结果表明,地下节根层数和总条数均与地下节根角度和面积呈负相关,地下节根的角度、面积、体积和干重等4个性状之间相互呈现显著正相关关系;不同年代的玉米地下根系结构存在差异,地下节根层数和总条数在年代的更替间表现出下降的趋势,地下节根角度和面积在年代更替间表现出上升的趋势,根干重和根体积在各年代间无显著差异;玉米地下根系结构在类群间也存在差异,旅大红骨类群的6个地下节根性状值均高于其余类群。全基因组关联分析共检测到26个SNP位点与地下节根层数、总条数、体积和干重性状显著关联(P<0.00001),其中11个显著关联位点定位于前人报道的根系QTL区间内,2个显著关联SNP在地下节根层数和总条数中均被检测到。基于显著关联SNP位点共挖掘到177个候选基因,其中135个具有功能注释,Zm00001d037368可能为控制地下节根层数和总条数的一因多效候选基因。候选基因功能的富集分析结果显示,候选基因的功能主要涉及植物体内的代谢调节、应激反应、运输活性、催化活性、结合蛋白及细胞成分等。【结论】玉米自交系的根系结构在不同年代间和不同类群间存在不同程度的差异,采用全基因组关联分析策略挖掘控制玉米根系结构的相关遗传位点及候选基因,共检测到26个显著关联的SNP位点。  相似文献   

13.
黄瓜幼苗下胚轴长度GWAS分析及候选基因挖掘   总被引:1,自引:1,他引:0  
【目的】挖掘与黄瓜幼苗下胚轴长度显著相关的SNP位点及候选基因,为揭示下胚轴长度的遗传基础和分子机制提供理论依据,为短下胚轴分子标记辅助选择育种奠定基础。【方法】以95份黄瓜核心种质为试验材料,分别于2016年春季、2017年春季、2017年秋季和2018年春季在中国农业科学院南口试验基地塑料大棚进行种植,在两叶一心期调查黄瓜幼苗的下胚轴长度;利用Structure 2.3.4软件分析群体结构,Haploview 软件分析连锁不平衡的衰减;基于最优模型对下胚轴长度进行全基因组关联分析(GWAS),依据关联SNP位点的LD区间序列,预测与下胚轴长度相关的重要关联候选基因,并利用荧光定量PCR对预测基因进行表达模式分析。【结果】共检测到8个显著关联的位点(Hl1.1、Hl1.2、Hl2.1、Hl3.1、Hl3.2、Hl4.1、Hl5.1、Hl6.1),分别位于1、2、3、4、5、6号染色体,其中,Hl2.1、Hl3.1、Hl3.2、Hl5.1、Hl6.1等5个位点被重复检测到两次以上。通过分析关联SNP位点的LD区间序列,获得Csa1G074930、Csa1G475980、Csa2G381650、Csa3G141820、Csa4G051570、Csa3G627150、Csa5G174640、Csa6G362970 8个与黄瓜下胚轴长度有关的候选基因,其中既有光形态建成、泛素化、激素信号通路等调控基因,也有调控网络下游参与细胞生长发育,调节细胞大小,直接调控黄瓜下胚轴长度的基因。多基因在不同黄瓜材料中的有机分布,形成了具有不同下胚轴长度的黄瓜种质。基因表达分析显示Csa1G074930、Csa1G475980、Csa2G381650、Csa4G051570、Csa5G174640在短下胚轴材料中高表达。Csa3G141820、Csa3G627150在长下胚轴材料中高表达。【结论】检测到Hl1.1、Hl1.2、Hl2.1、Hl3.1、Hl3.2、Hl4.1、Hl5.1、Hl6.1等8个与黄瓜下胚轴长度密切关联的SNP位点,挖掘到Csa1G074930、Csa1G475980、Csa2G381650、Csa3G141820、Csa4G051570、Csa3G627150、Csa5G174640、Csa6G362970等8个调控下胚轴长度的候选基因。  相似文献   

14.
【目的】解析白菜类作物开花时间的调控位点,定位白菜类作物开花时间相关的候选基因,为白菜类作物抽薹开花时间的遗传改良提供依据。【方法】以116份白菜类作物组成的自然群体作为研究材料,分别种植在温室与露地2个独立的环境中,进行开花时间调查。同时,提取试验材料的DNA样品进行深度为1.2x的重测序,对测序数据用Pooled Mapping法进行过滤、与参考基因组比对,获得全基因组高密度SNP集合。经过条件过滤后,对高质量的SNP集合进行生物信息分析,包括试验材料的群体结构分析和全基因组连锁不平衡分析。从高质量的SNP集合中,随机挑选出2 000个变异位点,用Phy ML软件以最大似然法对116份试验材料进行系统发育树分析。用全部的高质量SNP集合位点通过软件Haploview进行全基因组连锁不平衡分析。最后,将高质量的SNP集合与开花时间数据结合,通过TASSEL和GAPIT软件包以及R程序语言进行全基因组关联分析。根据强关联峰值信号点位置和连锁不平衡区间定位开花时间候选位点,再通过白菜与同源物种拟南芥的基因共线性关系以及基因功能注释分析来预测白菜类作物开花时间相关的候选基因。【结果】不同种植条件下、不同类型的白菜类作物在开花时间上存在广泛差异。试验材料在露地环境下的开花时间高峰期明显早于温室环境下的材料;试验材料在露地环境下的开花时间总体表现出偏正态分布,而在温室环境下,开花时间各个阶段呈现出较为均衡的分布。温室与露地环境下的开花时间呈显著正相关。通过生物信息学分析最终得到的高质量SNP位点共103万个。试验材料的群体结构分析表明在系统发育树上各亚群内部分布较为集中,不同亚群之间的分布与材料的地理起源密切相关。全基因组衰减平均LD为2.3 kb,表明在116份白菜类作物构建的群体内存在较为频繁的重组和突变。对不同条件下的开花时间进行全基因组关联分析,用复合模型检测到54个(P4)强关联峰值信号点,一般模型检测到87个(P5)。通过进一步分析强关联信号点的连锁不平衡(linkage disequilibrium,LD)区段,得到存在强连锁关系(r20.33)的峰值信号点共33个(温室环境下27个,露地环境下19个)。其中,在温室与露地环境下的共定位位点13个。根据33个关联候选位点,再通过白菜与同源物种拟南芥的基因共线性关系以及基因功能注释分析筛选出白菜类作物开花时间相关的候选基因14个,其中温室与露地环境下共定位候选基因3个(FUL、PHYB和FPF1)。在露地条件下定位到开花关键基因FT1。【结论】不同条件下开花时间的相关性分析表明,遗传效应在开花早晚中起着决定性作用。全基因组关联分析共鉴定出33个与开花时间相关的显著关联信号。通过连锁不平衡分析、白菜与同源物种拟南芥的基因共线性关系以及基因功能注释分析初步鉴定出14个白菜类作物开花时间相关的候选基因。  相似文献   

15.
【目的】 研究不同基因型海岛棉萌发期的耐盐特性,分析盐胁迫下海岛棉种子活力指标的变化规律,为培育耐盐海岛棉品种提供支持。【方法】 以35份海岛棉种子作为材料,以NaCl溶液(150 mmol/L)模拟盐胁迫条件,统计分析发芽势、发芽率、根长、下胚轴长、苗鲜重、苗干重6项指标。【结果】 (1)盐胁迫条件下,海岛棉种子各指标受影响程度不同,其中下胚轴受影响程度最大,苗干重受影响最小;(2)与对照样本比较,其中发芽势、发芽率、根长、下胚轴长、苗鲜重5项指标最大值、最小值、平均数均低于对照,苗干重的最大值、最小值、平均值则高于对照。【结论】 35份海岛棉种子分为3类:第Ⅰ类属于高耐盐品种,第Ⅱ类属于中等耐盐品种,第Ⅲ类属于低耐盐品种。  相似文献   

16.
【目的】甘蓝型油菜波里马细胞质雄性不育(pol CMS)在中国已被广泛应用于杂交种育种,其育性恢复程度表现出受1对主效基因的控制,并受微效修饰基因的影响。通过全基因组关联分析方法挖掘育性恢复位点,并对候选基因进行比较分析。【方法】通过芸薹属60K SNP芯片对308份甘蓝型油菜自然群体进行基因型分型,并用pol CMS系301A作母本,与上述材料分别进行杂交得到308份F1,每份F-_1分别于2013年和2014年进行种植,每年2次重复,于始花期根据花粉育性和花蕊发育情况调查F1植株的育性等级,同时对测交父本自然群体进行群体结构分析和亲缘关系评估,并结合测交父本的基因型分型结果和F_1的育性等级进行全基因组关联分析(GWAS)。从GWAS分析中显著的SNP左右100 kb区间或与显著SNP处于同一单体型块(R~20.5)的区间内预测候选基因,并对候选基因进行QTL比较分析和单体型或等位基因的效应分析。【结果】方差分析结果显示,两年F_1的育性等级存在显著差异(P0.01),但相关分析发现,两年的育性等级存在显著的正相关(r=0.52,P0.001)。群体结构分析显示,所有测交父本被分为3个亚群(冬性、春性和半冬性),亲缘关系分析发现,任何2个材料之间平均亲缘关系值为0.072,73%的任意材料间亲缘关系值小于0.1,其中,约53%的材料亲缘关系值为0。GWAS分析共检测到13个与育性恢复程度显著关联的SNP,构成了6个候选区间,分别位于A01、A09、C03、C06和C08 5条染色体上,单个SNP解释的表型变异介于2.53%—9.96%。从中共预测到6个与育性恢复位点相关的候选基因,其中4个编码的蛋白含有恢复基因特有的PPR保守基序。共线性分析发现,4个候选基因中的2个(Bna A09g46700D和Bna C08g40710D)位于A09和C08染色体部分同源区间,且与已克隆的pol CMS育性恢复位点ORF2同源。另外2个新鉴定到的候选基因(Bna C03g45840D和Bna C06g13000D)连锁的SNP等位基因或单体型变化都与育性等级显著相关(P0.001)。【结论】通过GWAS分析鉴定到多个与油菜育性恢复有关的候选基因,开发基于与这些基因连锁位点或SNP的功能标记将有助于对该不育系统进行恢复系和保持系的筛选。  相似文献   

17.
【目的】 对铃重、衣分、单株铃数和籽指等棉花产量构成因素性状进行全基因组关联分析(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个可能与产量性状相关的候选基因。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号