首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到13条相似文献,搜索用时 437 毫秒
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  
【目的】油菜高产是育种工作的主要研究目标之一。角果密度、主花序有效角果数等性状与产量都有显著或极显著的正相关关系,是油菜高产育种考查的主要性状。为揭示油菜角果密度及其相关性状的遗传机理和分子机制奠定基础。【方法】以不同遗传背景和地理来源的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个与主花序角果密度及其相关性状有关的候选基因。  相似文献   

3.
【目的】解析甘蓝型油菜发芽期根和下胚轴发育及耐盐性的调控位点,筛选油菜耐盐性相关的候选基因,可为油菜耐盐性改良提供依据。【方法】以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个重要候选基因。  相似文献   

4.
【目的】适合机械化收获是当今油菜育种改良和遗传研究的重要目标。该研究以一个自然变异产生的油菜有限花序(denterminate inflorescence 1,di1)突变体为研究对象,通过分析有限花序的遗传模式,开展有限花序性状的基因定位和克隆,以期发掘候选基因,为培育适合机械化收获的油菜新品种提供新思路和新材料,为揭示油菜有限花序遗传机制奠定基础。【方法】以一个稳定遗传的有限花序突变株系FM8与野生型自交系FM7开展正反交,观察F1和F2后代的花序形态,分析有限花序性状的遗传模式。在F2群体中挑选20个有限花序单株和20个野生类型单株构建混合池,对混合池和亲本开展20×和10×覆盖度的全基因组重测序,定位有限花序性状的关联区间。根据关联区间对应到拟南芥基因组的共线性区段和基因注释信息,预测候选基因,并对候选基因进行同源克隆,发掘序列变异,筛选关键基因。【结果】油菜有限花序突变性状表现为初花期主花序和侧枝花序顶部形成一个或若干个顶生花,花序无限生长受阻,导致结角期主枝和侧枝有封顶特征即有限花序。有限花序突变株系与野生型正反交F1均表现为野生型,F2代无限花序与有限花序的分离比符合13﹕3,说明有限花序的遗传受2对隐性基因和1对隐性上位抑制基因互作控制。对混合池及亲本开展全基因组重测序,得到30 123个单核苷酸多态性(SNPs)标记和107 636个插入缺失标记(In Dels)标记,用于有限花序性状的全基因组定位。定位结果共检测获得7个显著关联区间,分布于油菜A08、A09、A10、C08和C09共5条染色体。其中,A10染色体上的关联区间峰值最高,是控制有限花序性状的主效位点。并且,A10染色体关联区间内的14.36—15.07 Mb的区域与C09染色体2个关联区间显示高度同源性。候选基因预测发现位于A08、A09、A10、C08和C09的5个关联区间包含有8个候选基因,包括TERMINAL FLOWER 1(TFL1)、FLOWERING LOCUS C(FLC)、ATBZIP14(FD)、MULTICOPY SUPPRESSOR OF IRA1 4(FVE)和SCHLAFMUTZE(SMZ)。基因序列分析表明di1突变体TFL1、FVE和SMZ的基因编码区存在序列变异,并导致蛋白序列变异。【结论】油菜有限花序突变由2对隐性基因和1对隐性上位抑制基因互作控制。与有限花序性状显著关联的区间有7个,其中,位于染色体A10和C09的关联区间具有高度同源性。TFL1、FVE和SMZ被推断为有限花序性状的候选基因。  相似文献   

5.
【目的】甘蓝型油菜波里马细胞质雄性不育(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的功能标记将有助于对该不育系统进行恢复系和保持系的筛选。  相似文献   

6.
基于SNP遗传图谱定位甘蓝型油菜千粒重QTL位点   总被引:1,自引:0,他引:1  
【目的】甘蓝型油菜籽粒重量是构成油菜单株产量的三大因素之一(单株有效角果数、每角果粒数、粒重),是重要的育种目标。通过对5种环境下甘蓝型油菜千粒重进行QTL定位分析,寻找甘蓝型油菜千粒重的QTL及影响本甘蓝型油菜群体千粒重的候选基因。【方法】利用重组自交系群体在德国吉森、重庆北碚5种不同的环境下,测定各株系天然种子千粒重。利用重庆市油菜工程技术研究中心实验室构建的SNP高密度遗传图谱扫描5种环境中的千粒重QTL。该遗传图谱包括2 795个SNP位点,覆盖甘蓝型油菜基因组1 832.9 cM,标记之间的平均距离为0.66 cM。利用Windows QTL Cartographer2.5复合区间作图法对千粒重进行QTL定位。将49个拟南芥粒重相关基因与QTL对应置信区间序列进行同源比较分析(E值<1E–21),找出可能与甘蓝型油菜千粒重关联的候选基因。【结果】5种环境中千粒重变异范围较大,且均呈现正态分布,符合QTL定位要求。在5种环境之间千粒重均表现出正相关,其中,2013北碚与2012北碚、2008年吉森达到极显著水平,相关系数分别为0.248和0.249;2012年北碚与2010年北碚、2011年北碚及2008年吉森达到显著相关,相关系数分别为0.226、0.397和0.190。5种环境中共检测到14个QTL,分布在9条染色体,其中,C03染色体3个,A06、A07和C01各有2个,A03、A05、A08、A10和C02染色体上各有1个,LOD值在2.57-6.05,单个QTL解释的表型变异为4.64%-14.13%。与拟南芥粒重基因进行同源性分析,有16个粒重相关基因落在8个QTL置信区间,匹配E值介于0-2E-21。其中QTL qTSWA07-2区间内筛出7个粒重基因。粒重基因TTG2qTSWA03-1qTSWC02-1 2个QTL区间内均被检测到。AHK3qTSWA07-2qTSWA08-1qTSWC01-1区间内被检测到。【结论】利用该套油菜60K芯片准确定位了5种环境条件千粒重的QTL位点,与拟南芥粒重基因比对出该群体油菜粒重基因,该结果有利于不同材料在使用该套SNP芯片分析及对千粒重QTL位点的比对和候选基因的分析。  相似文献   

7.
【目的】探究不同氮素供应环境下与小麦苗期生物量及氮效率相关性状显著关联的SNP位点,预测相关候选基因,为小麦氮效率的基因克隆及其在育种中的应用提供参考。【方法】以134个小麦品种(系)组成的群体为供试群体,设置低氮、正常氮和高氮3个处理,各处理重复4次,并在2年(2013和2014年)进行了2次完全重复的营养液培养试验。试验对小麦苗期生物量及氮效率相关的14个性状进行了表型鉴定,采用MLM+K+Q混合线性模型,利用90K SNP芯片对小麦生物量及氮效率相关性状进行全基因组关联分析(genome-wide association study,GWAS),获得显著关联的SNP位点。【结果】与正常氮处理相比,低氮处理条件下,根系、地上部及植株氮含量和氮积累量均显著下降,而根生物量和根、植株氮效率均显著增加,高氮处理下,几乎所有鉴定性状均显著增加;14个性状的广义遗传力均在40%以上,其中,植株总干重的遗传力最高(95.73%)。利用9 329个SNP标记进行关联分析,共检测到838个SNP标记位点与供试材料的14个性状存在显著关联(P≤0.001),分布在21条染色体上。有435个(51.91%)SNP标记位点仅在一个关联分析环境中被检测到;有403个位点至少在2个处理环境(包含均值环境)中被检测到与同一性状显著关联(稳定关联标记)。其中8个SNP标记位点至少在3个环境中被检测到。在4个环境下(包括均值环境)均检测到的稳定关联位点有2个:Kukri_c65481_121和tplb0025f09_1052,分别与植株总氮利用效率(total nitrogen use efficiency of plant,TNUE)和根系总氮利用效率(root nitrogen use efficiency,RNUE)显著关联;同时与至少6个性状(生物量及养分效率相关性状)显著关联的SNP标记位点共5个,分别位于1A、1B(3)和2A染色体上;根据小麦基因组注释及LD衰减水平,在同时定位了6个性状的5个SNP位点和2个多环境(4个环境)稳定关联的SNP位点的214 kb的基因组区域中共筛选候选基因84个,根据已知克隆氮效率基因的编码蛋白类型、候选基因功能注释信息及利用植物比较基因组学资源库蛋白序列同源比对分析,筛选到3个候选基因与生物量及氮效率相关。【结论】不同氮素处理显著影响小麦苗期生物量、氮效率相关性状及其相关QTL的表达,大多数SNP位点仅在1个氮素检测环境中被检测到,但也存在环境稳定性较强的位点。生物量及氮效率相关性状之间存在显著相关关系,并在一定程度上受到相同的QTL/基因控制。  相似文献   

8.
【目的】探索甘蓝型油菜不同发育时期种子和角果皮硫苷含量的动态变化规律,通过转录组关联分析解析控制硫苷含量变异的遗传机制。【方法】在由113个油菜品种构成的自然群体中,利用高效液相色谱法检测不同发育时期油菜种子及角果皮中的硫苷含量。提取油菜幼嫩叶片mRNA并进行mRNA-Seq测序,开发了355 536个单核苷酸多态性标记(SNP)和116 098个基因表达量标记(GEM)。将开发的SNPs和GEMs分别导入到混合线性模型(MLM)和回归分析模型中进行关联分析,获得与油菜硫苷表型变异显著相关的遗传位点,进一步通过序列比对和功能预测筛选确定候选基因。【结果】授粉后15 d(15 DAP),油菜种子和角果皮的硫苷含量变异范围分别是1.69-20.45和1.47-25.23 μmol·g-1。25 DAP种子、25 DAP角果皮以及成熟种子的硫苷变异范围分别是2.17-147.21、0.73-130.77和8.87-111.83 μmol·g-1。后两个时期(即25 DAP和成熟期)的硫苷含量表现出较大变异,适合进行转录组关联分析。基于mRNA-Seq测序数据,经质量控制后获得256 397个高质量的SNP标记(MAF>0.01)和53 889个GEM标记(平均基因表达量>0.4)。利用SNP标记对成熟种子、25 DAP种子和25 DAP角果皮硫苷含量进行关联分析,分别检测到167、158和3个显著关联位点(-log10P>6.71)。其中,与成熟种子显著关联的SNP标记形成5个明显的关联峰,分别位于A2、A9、C2、C7和C9染色体上。同时,利用GEM标记进行的关联分析中,分别检测到127、16和24个与成熟种子、25 DAP种子、25 DAP角果皮硫苷含量显著关联的位点(-log10P>6.03),这些位点主要分布在A2、A8、A9、C2和C9染色体上并形成明显的关联峰。其中,成熟种子和25 DAP角果皮在A9、C2和C9染色体上形成共同的关联峰。通过与公共数据库的基因组序列比对和功能注释,共筛选到295个功能基因,其中25个基因涉及硫苷代谢网路途径,直接参与硫苷的调控。其余基因虽然没有直接参与硫苷代谢过程,但是它们大多属于转录因子基因、刺激响应因子基因,涉及细胞过程或是具有催化活性,因而推测在硫苷积累中同样起到重要作用。【结论】油菜角果皮中的硫苷含量较高,而且与种子硫苷含量显著正相关,是硫苷合成或转运的重要器官。共检测到328个SNP显著标记和144个GEM显著标记,筛选到25个与硫苷合成或者转运有关的功能基因以及73个功能未知的新基因。  相似文献   

9.
【目的】近几年随着观光农业的兴起,花色的选育和改良已成为甘蓝型油菜种质资源鉴定和材料创制的重要研究方向。以甘蓝型油菜黄白花分离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标记。  相似文献   

10.
鸡胸腺重和脾脏重性状的全基因组关联   总被引:2,自引:2,他引:0  
【目的】全基因组关联分析(genome-wide association study,GWAS)是复杂性状和疾病相关基因定位的新策略。【方法】试验利用鸡60K SNP芯片对来自50个公鸡家系的728只北京油鸡进行SNP分型检测,采用全基因组关联分析方法对影响100日龄胸腺重和脾脏重的基因进行定位研究。【结果】结果发现24个达5%全基因组水平显著的位点,与100日龄胸腺重和脾脏重显著相关,并在这些位点附近发现Janus kinase 1(JAK1)、 zinc finger DHHC-type containing 8(ZDHHC8)、vav 3 guanine nucleotide exchange factor(VAV3)、SATB homeobox 1(SATB1)等候选基因;84个与这两个性状潜在关联同时达到5%染色体水平显著的位点。【结论】利用GWAS分析策略筛选和鉴定的重要突变位点及候选基因,将为揭示鸡免疫器官发育的分子调控机制和抗病育种分子标记辅助选择提供必要的分子基础。  相似文献   

11.
【目的】解析白菜类作物开花时间的调控位点,定位白菜类作物开花时间相关的候选基因,为白菜类作物抽薹开花时间的遗传改良提供依据。【方法】以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个白菜类作物开花时间相关的候选基因。  相似文献   

12.
【目的】研究白菜型油菜的自交亲和性、角果及种子相关性状的QTL定位,为白菜型油菜自交亲和品种、优良自交系的选育及其产量和品质的遗传改良提供依据。【方法】以自交亲和性较差的菜薹L58和自交亲和性较好的白菜型油菜R-O-18构建的包含117个株系的RIL群体为试验材料,利用已构建的包括372个InDel标记的分子遗传图谱,采用区间作图法(IM)对花期亲和指数、结角率、角果长度、角果喙长度、角果长宽比、角果喙长度/角果长度、每角果种子数、千粒重及种子颜色9个性状分别进行QTL分析。【结果】RIL群体在自交亲和性、角果及种子相关的9个性状上表现为连续变异,并且变异幅度较大,均呈正态分布或偏正态分布,具有典型的数量性状遗传特点。花期亲和指数与结角率、角果长度、角果喙长度、角果长宽比、每角果种子数、千粒重均存在极显著的正相关,其中花期亲和指数与每角果种子数间的相关系数最大,达到0.8487。与其他几个性状的相关性大小依次为:千粒重>结角率>角果喙长度>角果长度>角果长宽比。每角果种子数及千粒重均与结角率、角果长度、角果喙长度、角果长宽比呈现极显著正相关,而每角果种子数与千粒重之间也存在极显著的正相关(0.6477)。结角率、角果长度、角果喙长度与角果长宽比之间呈现显著正相关;而角果喙长/角果长除了与角果长呈负相关,角果喙长度存在极显著正相关性外,与其他性状均不相关。共检测到12个QTL位点,其中6个QTL位于A09连锁群,10个QTL解释大于10%的表型变异。2个控制花期亲和指数的QTL均位于A09连锁群上,分别可解释13.1%和16.7%的表型变异;6个角果相关性状的QTL,分别位于A06、A09和A10连锁群上,单个QTL可解释13.4%-17.7%的表型变异;4个种子相关性状的QTL,分别位于A02、A05、A06和A09连锁群上,可解释7.9%-42.1%的表型变异,位于A09连锁群的种子颜色QTL为主效位点,其加性效应值为1.22。【结论】共检测到12个控制自交亲和性(2个)、角果(6个)及种子相关性状(4个)的QTL,其中,2个控制自交亲和性的QTL,与前人得到的调控自交亲和性主位点(S位点)不同,可能为控制自交亲和性的非主效QTL位点。1个种子颜色QTL为主效位点。  相似文献   

13.
【目的】 对铃重、衣分、单株铃数和籽指等棉花产量构成因素性状进行全基因组关联分析(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号