小麦穗发芽性状的全基因组关联分析

为了解小麦穗发芽的遗传机制,发掘与小麦穗发芽相关的候选基因,采用整穗发芽法对来自全国的207份小麦品种(系)进行表型鉴定,并结合小麦90K SNP基因芯片,通过TASSLE软件的MLM (Q+K)模型对小麦的整穗发芽率进行全基因组关联分析(genome-wide association study, GWAS)。研究结果表明,不同年份间小麦品种(系)的穗发芽表现出丰富的表型变异,变异系数为0.34和0.25,多态性信息含量(polymorphicinformationcontent,PIC)为0.01～0.38,全基因组LD衰减距离为3 Mb。群体结构分析和主成分分析表明, 207份小麦品种(系)所构成的自然群体结构简单,可分为3个亚群。GWAS检测结果显示,在不同环境间共检测到34个与小麦穗发芽显著关联的SNP标记位点(P≤0.001),分布在小麦3A、3B、4A、4B、5D、6A、6B、6D、7B和7D染色体上,单个位点可解释5.55%～11.63%的表型变异, 16个标记位点在两个及以上环境下均被检测到,其中6B染色体上的标记wsnp_Ex_c14101_22012676在E1、E2及平均环境下被共同检测到,属稳定遗传的位点。通过对表型效应值大且稳定遗传的关联位点进行挖掘,共筛选到13个与小麦穗发芽相关的候选基因。其中TraesCS3A01G589400LC、TraesCS6B01G138600/TraeCS6B01G516700LC/TraesCS6B01G548900LC、TraesCS6D01G103600及TraesCS7B01G200100等基因通过调控植物内源激素-脱落酸(abscisicacid,ABA)的灵敏性进而影响种子的休眠;TraesCS3B01G415900LC、 TraesCS6A01G144700LC及TraesCS6B01G294800等基因编码的F-box蛋白在植物激素的信号转导、光信号转导以及花器官发育等生理过程中起重要作用; TraesCS6A01G108800、TraesCS6B01G138200/TraesCS6B01G293700基因编码Myb转录因子家族蛋白能调控种子中类黄酮的生物合成,对籽粒颜色有重要影响,这些候选基因是与小麦穗发芽相关的重要基因。     

小麦株高相关性状与SNP标记全基因组关联分析

株高是影响小麦产量和控制倒伏的重要因素,研究小麦株高相关性状的遗传机制对高产育种具有指导意义。以205份中国冬麦区小麦品种(系)为材料,利用分布于小麦全基因组的24 355个单核苷酸多态性(SNP)标记对株高相关性状进行关联分析。共发现38个与株高相关性状显著关联(P0.0001)的SNP,分布在1B、2A、2B、3A、3B、3D、4A、4B、5A和6D染色体上。其中,11个位点至少在2个环境中稳定表达,可用于开发CAPS标记。同时,发掘了一批株高性状相关基因的优异等位变异,如降低株高的等位变异Bob White_c48009_52,平均降低株高12.9 cm;控制穗下节间长的等位变异BS00039422_51-C和IAAV1698-A,分别调控穗下节间长5.9 cm和6.6 cm。本研究发掘的控制小麦株高基因位点为在分子水平上研究小麦株高复杂性状提供了有价值的参考。     

调控小麦重要农艺性状基因的染色体定位

为了研究调控小麦重要农艺性状基因的染色体分布,进一步解析控制小麦重要农艺性状基因的遗传基础,本研究以一套"中国春-人工合成小麦"染色体代换系为材料,在大田条件下,对不同株系小麦的株高、穗长和穗数等重要农艺性状进行调查。通过同一性状不同株系间的差异显著性分析,对调控株高、穗长和穗数的基因进行了染色体定位。研究结果表明:1 B、2 B、3 A、4 A、6 A、6 D和7 D染色体上可能携带有调控小麦株高的位点;1 B、1 D、2 B、2 D、5 A、6 A、6 D和7 D染色体上可能携带有调控小麦穗长的基因位点;3 A和6 D染色体上可能携带有调控小麦成穗数的位点。本研究有利于进一步理解调控小麦农艺性状的遗传基础,并对小麦分子育种实践具有一定的指导意义。     

小麦茎秆断裂强度相关性状QTL的连锁和关联分析

小麦茎秆断裂强度与倒伏特性关系密切,并对产量有很大影响。本研究旨在解析茎秆断裂强度的遗传机制,开发与该性状紧密连锁/关联的分子标记。利用山农01-35′藁城9411重组自交系(RIL)群体(含173个F8:9株系)和由205个品种(系)构成的自然群体,借助90 k小麦SNP基因芯片、DArT芯片及传统分子标记技术,在2个环境中对两群体的茎秆断裂强度相关性状进行连锁分析和全基因组关联分析。利用已构建的高密度连锁图谱,在4B染色体的TDURUM_CONTIG63670_287–IACX557和EX_C101685–RAC875_C27536等区段上,检测到9个控制小麦茎秆断裂强度、株高、茎秆第2节间充实度、茎秆第2节壁厚相关性状的加性QTL,可解释表型变异9.40%~36.30%。同时,利用包含24 355个SNP位点的复合遗传图谱,在自然群体中检测到37个与茎秆断裂强度相关性状(P0.0001)的标记,分别位于1A、1B、2B、2D、3A、3B、4A、4B、5A、5B、5D、6B、7A、7B和7D染色体,可解释表型变异7.76%~36.36%。在4B染色体上,以连锁分析检测到控制茎秆断裂强度的RAC875_C27536与关联分析检测到的Tdurum_contig4974_355标记,在复合遗传图谱上的距离为6.7 cM,说明该区段存在控制小麦茎秆断裂强度的重要基因。     

普通小麦氮素利用效率相关性状全基因组关联分析

氮元素在粮食作物生长和发育过程起着不可替代的作用。发掘氮素利用效率相关基因对于提升小麦产量、减少环境污染具有重要意义。植株根系构型(root system architecture,RSA)代表着根系的结构及空间造型,显著受氮素水平影响。本研究在正常供氮和缺氮两种氮素水平下,对160份来自黄淮冬麦区和北部冬麦区普通小麦品系的根系构型相关性状(总根长、总根表面积、总根体积、平均根直径和根尖数)进行统计,并结合小麦660K SNP(single nucleotide polymorphism)芯片基因分型数据对根系相关性状的相对值进行全基因组关联分析,以期发掘氮素利用效率相关位点。本研究检测到34个与氮素利用效率显著关联的SNP位点,可解释6.9%~15.4%的表型变异。关联位点在所有染色体均有分布,主要集中于1A、2B、3B、5B、6A、6B和7A染色体。11个位点与已报道位点重叠或接近,其他23个位点可能为新的位点。另外,在3B染色体上发现一个编码E3泛素连接酶的候选基因。     

基于SNP标记小麦自然群体遗传多样性及复合图谱的构建

本研究利用Illumina i Select 90K(81 587个SNP)基因芯片对中国冬麦区205份小麦育成品种(系)构成的自然群体进行基因分型,分析中国冬麦区小麦的遗传多样性并整合群体的复合遗传图谱。用于检测的81 587个SNP标记位点中,利用32 432个具有品种间多态性的位点进行群体遗传多样性分析,检测到64 864个等位变异,每个位点平均2个等位变异,SNP位点的多态性信息含量(PIC)变化范围为0.05~0.38,平均为0.26;聚类分析将205份小麦材料按亲缘关系划分为4个类群。复合遗传图谱包含24 355个SNP标记位点,覆盖小麦全基因组3 674.16 c M,单个染色体长度为118.91~241.38 c M,标记间平均遗传距离为0.15 c M;小麦的3个染色体组中,B基因组中包含的SNP标记最多(12 321,50.60%),其次是A基因组(9 523,39.10%),D基因组含标记最少(2 511,10.31%)。研究结果为利用该群体进行标记/性状间的关联分析提供遗传信息,并为小麦杂交亲本的选择和新品种培育提供参考。     

河南小麦品种穗粒数性状的动态变化及全基因组关联分析

穗粒数是小麦产量构成的重要因子,研究小麦穗粒数及其相关性状的变化规律与遗传机制对于高产育种具有指导意义。以198份河南种植或审定的小麦品种(系)为材料,在4个环境下鉴定穗粒数相关性状表现,选用660K SNP芯片检测材料的基因型并进行全基因组关联分析。发现穗粒数存在着广泛的变异,穗粒数、穗长、总小穗数、可育小穗数和不育小穗数的广义遗传力分别为77.43%、87.54%、85.45%、79.70%和68.08%,2002年以前、2003-2008年、2009-2014年和2015年以来育成品种的平均穗粒数分别为35.71个、36.85个、37.40个和37.98个,随着育种时期的推进总小穗数和可育小穗数也表现出逐渐提高的趋势。共发现41个和穗粒数等性状相关联的位点,分布在除了1D、3A、3D、4B和4D的16条染色体上,单个关联位点的表型变异贡献率(R2)范围为6.19%~20.83%,其中10个位点至少在2个环境中稳定表达。对关联位点进一步分析发掘了多个与穗粒数性状相关的优异等位变异,如7A上AX-109368219(A)能增加穗粒数1.86粒,2A上AX-111674224(C)能提高总小穗数0.96个,2D上AX-111582877(G)能提高可育小穗数0.73个。研究揭示近五十年河南小麦穗粒数性状的动态变化及其规律,获得了多个和穗粒数等性状关联位点及等位变异,为小麦产量性状的遗传改良及分子设计育种提供了技术支持和基因资源。     

野生二粒小麦主要农艺特性融入普通小麦的全基因组关联分析

野生二粒小麦(Triticum turgidum ssp.dicoccoides)是普通小麦的四倍体祖先种,具有广泛的基因型变异,是改良普通小麦的重要种质资源。本研究对不同年份和地点四个环境下的161份野生二粒小麦渗入系材料的株高、分蘖数、小穗数、抽穗期、开花期和千粒重进行表型鉴定,并利用覆盖全基因组的13,116个DArT标记对各农艺性状进行全基因组关联分析,以期发掘性状显著关联标记及相关候选基因。本研究共检测到147个与6个农艺性状相关的稳定标记。其中,一些关联标记与抽穗期和开花期同时相关,并且在2B染色体上聚集成簇。在野生二粒小麦渗入系群体中共推定了21个与农艺性状相关的候选基因,其中位于7A染色体与千粒重显著相关的候选基因与细胞周期蛋白有关。这些标记和候选基因可为克隆优异农艺性状相关基因提供重要信息,从而为野生二粒小麦在普通小麦背景中的遗传改良综合利用提供依据与指导。     

中国西门塔尔牛SCD1基因多态性与肉质性状的相关分析

利用Illumina HD SNP技术对484头中国西门塔尔牛SCD1基因的遗传变异及其与部分肉质性状的相关性进行了分析.结果表明,该群体中存在7个SCD1基因的SNP位点.SCD1基因第4内含子存在A8605G、A9229G与A10050C 3个SNP位点;第5内含子存在1个SNP位点A12304G;3' UTR区存在A13655G、C14790T与A15565G 3个位点.A8605G位点的B等位基因为优势等位基因,频率为0.748,其他位点等位基因频率在0.4 ～0.6之间.x2检验结果显示,SCD1基因的所有SNPs均处于哈代-温伯格平衡状态(P＞0.05).SCD1基因的7个SNPs位点全部为中度多态,杂合度和有效等位基因数均较高.SCD1基因SNP多态性与肉质性状的关联分析表明,SCD1基因对肌内脂肪含量、剪切力、大理石花纹等级和脂肪颜色有一定的影响,可以作为改善牛肉质量进行标记辅助选择的重要候选基因.     

基于SNP分子标记的泰山/泰科麦系列小麦遗传解析

对不同年份育成的21个小麦品种（系）进行全基因组扫描,通过分析遗传距离和染色体区段/位点,明确其亲缘关系远近和遗传差异。分析可知,获得的2029个SNP基因位点在B基因组拥有较高的遗传多样性,其次是A和D基因组;在7个同源群中,第3和第6同源群呈现出较高的遗传多样性,而第1和第4同源群的遗传多样性较低;21条染色体中,3A、1B、6B染色体的遗传多样性较高,而1A、6A的遗传多样性偏低。对21份供试材料依据审定（育成）年份分析其群体的平均遗传距离,不同年份品种间的平均遗传距离先增大后减小,遗传多样性逐渐降低;21份供试材料间的遗传相似系数在0.69~0.99之间,大致可聚为4个类群,同一年份的品种一般聚在一起,与其系谱关系吻合。构建并分析供试材料的基因型图谱发现,00s、10s和现在育成的小麦品种（系）共有SNP和共有染色体区段分别主要在A、D和B基因组,对应已发表性状同不同年份育种目标吻合。同时发现21份供试材料均含有25个共同SNP位点,分布在1A、5A、6A、7A、2B、3B、6B、1D、2D、3D和7D染色体上,且每条染色体上分布的SNP位点数目均不相同,通过对应已发表性状进一步证实在品种（系）组配与选育过程中注重产量、株高、分蘖数、抽穗期、灌浆速率和抗病等性状的选择。以上研究结果可为今后小麦新品种组配和选育提供参考依据。     

小麦芒长抑制基因B2的精细定位与候选基因分析

芒是小麦重要的穗部器官和形态特征,是小麦长期进化和适应环境的结果,对产量和抗旱性等具有重要影响。目前,对麦芒的遗传与发育还缺乏系统的研究,相关基因克隆或精细定位的研究尚未见报道。本试验利用短芒材料‘六柱头’与长芒材料‘石矮1号’构建的F2群体(SL-F2)对芒的遗传与发育进行了研究。细胞学观察表明,短芒主要是由细胞长度变短引起;遗传分析表明,‘六柱头’的短芒由显性单基因控制;借助Wheat660K SNP芯片的BSA分析和SL-F2群体的精细定位,确定‘六柱头’的芒长抑制基因是前人报道的B2位点,并将其定位到6B染色体4.84Mb的物理区间(471.28～476.12 Mb)内,该区段在中国春与矮抗58间具有良好的共线性。在B2定位区间共有61个基因,其中5个在中国春穗部特异表达,TraesCS6B02G264400在中国春和Azhurnaya幼穗表达差异显著。这些研究结果为B2基因的克隆、小麦芒形成机理的解析及育种中的应用奠定了基础。     

砷胁迫下甘蓝型油菜苗期根、下胚轴和鲜重的全基因组关联分析

油菜是修复土壤重金属污染的理想作物,为筛选甘蓝型油菜耐砷性的显著关联单核苷酸多态性位点及相关候选基因,本研究以140份不同来源的甘蓝型油菜自交系为材料,测定和利用油菜60KSNP芯片对正常和砷胁迫条件下的相对根长(RRL)、相对下胚轴长(RHL)和相对鲜重(RFW)进行了全基因组关联分析。结果表明,与RRL、RHL和RFW显著关联的SNP位点分别为15、20和35个,单个SNP位点表型贡献率分别介于13.31%～24.39%、18.04%～33.82%和20.19%～25.06%之间;其中在A02、A07和C02染色体上同时存在与RRL、RHL和RFW显著关联的LD区间。基于油菜基因组信息在LD区间内共筛选到61个可能与砷胁迫相关的候选基因,其中PHT3;3、PHT1;9、GST、OTC5、NRAMP1和ZIP12等与重金属吸收和转运相关。实时荧光定量PCR分析结果表明, PHT3;3和PHT1;9是与甘蓝型油菜砷离子吸收转运相关的重要候选基因。本研究结果对于甘蓝型油菜耐砷胁迫机理的研究、性状的改良具有重要参考价值。     

国审小麦品种山农20高产稳产、广适性表现及遗传解析

山农20是山东农业大学利用常规育种与分子标记辅助选择育成的多抗、广适、高产且稳产小麦品种，先后通过国家黄淮南片、北片、新疆和甘肃审定。山农20含有长穗偃麦草和美国种质血缘，表现抗病、大穗和大粒等优异农艺性状。通过基因功能标记或与基因紧密连锁的微卫星标记的检测，分析了山农20含有的控制重要农艺性状的关键基因，得出山农20含有Pm2、Pm24和Pm30白粉病抗性基因，在田间表现出较好的白粉病抗性；含有Yr15和Yr26抗条锈基因、Lr21抗叶锈基因和Ses1抗纹枯病基因，在田间表现出较好的条锈病、叶锈病及纹枯病抗性；聚合了控制小麦产量三因素的基因位点，即粒重位点QGw6A-29和QTaGW4B.4、穗粒数位点QGns2B-2、穗长位点QS12D-3、分蘖成穗位点QMtw5D-1和QMtw6A，这与其高产、稳产和广适性是一致的。以上信息为深入认识山农20重要农艺性状分子机理提供了线索，对未来小麦遗传改良中高效利用该品种的重要基因具有实用价值。     

小麦转录因子基因TaNAC67参与调控穗长和每穗小穗数

NAC转录因子是植物特有的一类转录因子,在植物生长发育和逆境胁迫应答反应中发挥着重要作用。前期研究表明, TaNAC67参与对多种逆境胁迫的应答,过量表达能增强拟南芥的抗逆性。为进一步揭示其在调控小麦主要农艺性状发育方面的作用,本研究以36份普通小麦组成的高多态性群体为材料,测序分析了TaNAC67-6A、TaNAC67-6B、TaNAC67-6D序列多态性,发现TaNAC67-6A启动子区–1516 nt有1个A/G转换SNP,在–873～–748 nt处有1个126 bp的InDel; TaNAC67-6B启动子区–2014和–1916 nt处各有1个C/T转换SNP; TaNAC67-6D启动子区–1795 nt有1个T/G颠换SNP,编码区357 nt处有1个C/T转换SNP。根据多态性分别开发了功能分子标记,扫描由282份普通小麦构成的自然群体,并将基因型和表型性状检测结果进行关联分析,TaNAC67-6A、TaNAC67-6B的标记与表型性状无显著关联,而TaNAC67-6D的2个标记SNP-D-1和SNP-D-2分别与小麦穗长和每穗小穗数显著相关。单倍型分析发现,自然群体中存在3种主要单倍型,其中Hap-6D-3是增加穗长和每穗小穗数的最优单倍型,在我国小麦育种历史中受到了正向选择。转基因水稻表型分析发现,TaNAC67过表达能显著增加水稻主穗穗长、穗分枝和穗粒数,验证了小麦关联分析结果。因此,TaNAC67-6D可用于改良农作物穗部性状,其分子标记可用于小麦分子标记辅助选择育种。     

不同氮水平下水稻钾吸收及全基因组关联分析

以业已完成深度重测序的134份水稻地方种质资源为材料,在大田栽培条件下按不施氮(N0),施96 kg hm~(–2)纯氮(N1)和施192 kg hm~(–2)纯氮(N2) 3种氮肥水平,检测了分蘖盛期植株的钾含量、植株干重和钾积累。结果表明,水稻钾含量、植株干重和钾积累在N0、N1、N2三种氮肥处理下均呈正态分布,表型变异丰富。植株干重和钾积累与施氮量呈极显著正相关,钾含量与施氮量相关性不显著;钾含量与植株干重呈负相关,钾含量与钾积累相关性不显著,而植株干重与钾积累呈极显著正相关。在3个施氮水平下,籼稻的钾含量极显著低于粳稻,籼稻的干重和钾积累都极显著高于粳稻。全基因组关联分析发现,在3个氮肥水平下检测到12个显著相关位点,其中钾积累、钾含量和植株干重的显著相关位点分别有2、5和5个。在N1水平下,位于第6染色体上与钾含量相关的SNP (Chr6_1,524,776)的显著性峰候选区包含与钾离子转运蛋白互作的RUPO基因。根据钾含量的差异,鉴定出3个钾含量与低氮响应有关的SNP位点, 1个位点与高氮响应有关,而位于第10染色体上的显著性位点Chr10_2,822,026对低氮和高氮均有响应,该区域的4个候选基因的表达在不同氮水平间存在差异。     

宁麦9号与扬麦158株高及其构成因素的遗传解析

宁麦9号与扬麦158是我国长江中下游麦区的主栽品种和骨干亲本,长江中下游麦区近3年来审定品种中80%都是其衍生后代,研究其性状的遗传具重要意义。以宁麦9号与扬麦158为亲本构建的包含282个家系的重组自交系群体为材料,利用Illumina 90k芯片对群体进行基因型分析,建立高密度遗传图谱。连续3个生长季对株高及节间长度、穗长等株高构成因素进行测定,结合遗传图谱对株高及相关性状进行QTL定位,获得14个控制株高及其构成因素的稳定表达位点。通过进一步位置比对,聚焦到6个染色体区段,初步明确了各节间对株高的遗传调控机制。同时,将6个染色体区段中同源性较低的连锁标记转化为适用于高通量筛选的KASP标记,利用101份区域试验材料进行标记效应验证,结果显示聚合Qph-2D与Qph-5A.1两个位点具有较高的选择效率,继续聚合Q2A后,中选材料显著减少,可能降低选择效率;对Q2A与Q5A两个一因多效位点的选择建议以降低株高的等位变异为主;Qd1-5D可作为穗下节间(D1)的选择标记对株高展开优化选择。期望以上结果能为长江中下游麦区的小麦株高遗传改良提供帮助。     

玉米隐花色素CRY1b和CRY2基因转录丰度对不同光质处理的响应

玉米株高、开花期、产量、品质等性状与环境中的光密切相关。隐花色素是一类蓝光和近紫外光的受体, 主要参与植物的光形态建成及动、植物的生物钟调控。通过研究玉米隐花色素基因对不同光处理的表达模式, 可为进一步研究其对玉米光形态建成的作用奠定基础。本研究采用RT-PCR技术克隆了玉米ZmCRY1b和ZmCRY2基因; 利用生物信息学相关网站和软件对其编码蛋白的结构域及氨基酸进行了系统发育分析; 利用qRT-PCR分析了玉米自交系B73中ZmCRY1b和ZmCRY2基因在不同组织、以及响应不同光质及长日照和短日照处理的转录丰度。研究发现, 玉米与拟南芥、水稻和小麦的CRY蛋白有相同的结构域及较高的氨基酸序列的一致性, 表明它们具有相似的功能。ZmCRY1b和ZmCRY2基因主要在玉米的叶片中表达; 二者能迅速响应各种持续光质、黑暗到不同光质转换及长日照和短日照处理, 且ZmCRY1b在各种处理下的转录丰度均高于ZmCRY2, 可能暗示ZmCRY1b在玉米中功能更强。以上研究结果表明, ZmCRY1b和ZmCRY2基因均能有效地响应各种光质和光周期处理, 并在玉米的光形态建成中发挥重要作用。本研究为进一步探明ZmCRY1b和ZmCRY2基因的功能及其在玉米品种改良中的应用提供了研究基础。     

水稻高秆染色体片段代换系Z1377的鉴定及重要农艺性状QTL定位

株高是水稻重要的农艺性状, 往往与产量相关性状密切关联, 在水稻育种中有重要利用价值。本研究以日本晴为受体、缙恢35为供体亲本, 经表型和分子标记双重选择, 鉴定了一个水稻高秆染色体片段代换系Z1377。Z1377共含有18个代换片段, 平均代换长度为2.95 Mb。与日本晴相比, Z1377的株高、倒一节间至倒四节间长、穗长、一次枝梗数、二次枝梗数、粒长、每穗实粒数、总粒数显著增加; 粒宽显著变细, 有效穗数、结实率显著减少, 但仍达86.75%。用日本晴与Z1377杂交构建的次级F₂群体共检测到16个相关QTL, 分布于第2、第3、第4、第5、第6、第7和第9染色体。其中有8个可能与已克隆基因等位, 如GW2、EUI1、ZFP185等, 另8个如qPH3等尚未见报道。Z1377的株高由一个主效QTL (qPH3)和一个微效QTL (qPH5)控制, 其中qPH3的贡献率达28.59%。而且, 在F₂群体中, 高秆和矮秆基本呈现双峰分布, 经卡平方测验, 符合3∶1分离比, 表明高秆对矮秆显性, 并主要由qPH3负责。这将为该主效基因的精细定位和克隆奠定基础, 同时为进一步选育含2~3个代换片段的中高优良染色体片段代换系并应用于育种奠定基础。     

小麦骨干亲本繁6产量相关性状关键基因组区段的遗传效应

繁6是中国小麦育种最重要的骨干亲本之一, 明确其优良特性的遗传机制对小麦育种具有重要意义。本研究鉴定了39个繁6衍生品种的7个产量相关性状, 并利用全基因组SSR标记分析了繁6中控制这些性状重要遗传区段和基因位点在子代中的遗传效应。表型鉴定结果表明, 繁6产量相关性状在不同世代衍生品种中表现无显著差异, 表明这些性状在衍生品种选育过程中受到选择并稳定遗传。利用已获得的控制小麦产量相关“一致性” QTL区段的417个SSR标记进行分子扫描, 发现11个繁6特异SSR标记在其衍生后代中被高频率遗传。性状-标记关联分析表明, 21个来自繁6的特异SSR标记与产量相关性状极显著关联(P < 0.01)。同时鉴定出分别位于2A和5A染色体的Xgdm93.3-Xgwm526.2和Xbarc100-Xgwm156.1区段, 前者控制株高和小穗数, 后者控制千粒重。本研究证实, 上述两个与小麦产量相关性状关联的位点或区段在小麦产量育种进程中受到强烈的人为定向连续选择, 并在四川乃至西南麦区小麦产量育种中发挥了重要作用。     

用全基因组关联作图和共表达网络分析鉴定油菜种子硫苷含量的候选基因

油菜籽饼粕是畜禽养殖中重要的蛋白原料, 但饼粕中的硫苷是一种抗营养物质, 食用过多会对禽畜产生毒害, 因此挖掘油菜籽粒硫苷含量的候选基因对油菜种子低硫苷育种具有重要现实意义。本研究连续4年种植1个含157份材料的油菜自然群体, 结合重测序数据对种子硫苷含量进行全基因组关联分析(GWAS), 并对15份低硫苷和15份高硫苷材料进行种子发育早期的转录组测序, 通过权重基因共表达网络分析(WGCNA)鉴定种子硫苷含量的候选基因。用GWAS共检测到45个与种子硫苷含量显著相关的SNP, 单个位点解释的表型变异为13.5%~23.3%, 主要分布在A09、C02和C09染色体的3个区间中, 覆盖5个已知的硫苷代谢基因。用WGCNA分析发现高、低硫苷材料之间的2275个差异表达基因, 可分为12个基因模块, 其中1个模块的基因显著富集在已知的硫苷生物合成途径, 对该模块内163个基因的权重分析得到13个候选基因。经检测, GWAS和WGCNA共得到的18个候选基因中, 有14个候选基因的表达量与种子硫苷含量显著相关(r = 0.376~0.638, P<0.05)。用两种方法鉴定到1个共同的候选基因BnaC02g41790D (基因名MAM1), 与该基因连锁的5个SNP构成5种单体型, 等位基因效应分析发现, 自然群体中63%的材料(99/157)为Hap 5, 平均硫苷含量为50.79 μmol g^-1, 与另外4种单体型(95.04~110.28 μmol g^-1)存在极显著差异(P<0.01)。本研究结合GWAS和WGCNA两种方法鉴定了油菜种子硫苷含量的候选基因, 可为复杂性状候选基因的筛选提供参考。