大豆Glyma04G227700生物信息学分析及与Glyma08G11030互作研究

咖啡酸-O-甲基转移酶(caffeic acid-O-methyltransferase, COMT)是一种催化苯丙素类化合物羟基上氧原子的甲基化酶。前期研究表明,大豆抵御干旱胁迫的F-box基因Glyma08G11030编码蛋白可能与COMT蛋白Glyma04G227700发生互作,为了预测分析大豆COMT基因Glyma04G227700的功能及二者的互作关系,本研究克隆Glyma04G227700基因,并对其进行生物信息学分析,通过实时荧光定量检测方法分析其在大豆不同组织中的表达情况,并利用酵母双杂交系统验证蛋白互作情况。结果显示：Glyma04G227700基因全长1 095 bp,编码366个氨基酸,与野大豆(RZC17879.1 Glycine soja)亲缘关系较近。Glyma04G227700在大豆的根、茎、叶、子叶中都有表达,且在根和茎中表达量较高,在子叶中表达量最低。酵母双杂结果表明Glyma08G11030与Glyma04G227700蛋白不存在互作。     

基于Overview和物理图谱的大豆主茎节数候选基因挖掘

利用与大豆主茎节数性状相关的54个原始QTL位点,应用Overview方法首次以大豆参考基因组物理图谱为背景进行整合和分析,得到11个重演性较好的置信区间,分布在大豆D1b、C2、B1、F、L和I连锁群上,其中L连锁群重演性较好的置信区间较多。对得到的候选区段进行基因注释得到488个候选基因,其中Glyma.11G087300、Glyma.20G014300、Glyma.13G221400、Glyma.06G243500、Glyma.13G052900、Glyma.13G052700参与植物激素信号转到通路(ID:Ko04075),推测这6个基因通过该通路的赤霉素途径和生长素途径参与大豆主茎节数的遗传调控。本研究所挖掘到的与茎生长发育及主茎节数直接相关的通路和候选基因能够为构建理想株型和大豆分子辅助育种提供新思路。     

大豆耐胁迫相关蛋白编码基因GmSAMDC1的克隆及表达分析

S-腺苷甲硫氨酸脱羧酶(S-adenosine methionine decarboxylase,SAMDC)是亚精胺和精胺合成的关键酶,更是多胺生物合成的限速酶之一,在植物耐胁迫反应中发挥重要调控作用。为研究大豆SAMDC编码基因的结构和表达特性,本研究从抗病大豆品种科丰1号中克隆了位于大豆基因组2号染色体上的GmSAMDC1(Glyma. 02G128000)基因,分析结果表明:其完整ORF长度为1 068 bp,编码一个由355个氨基酸组成的包含1个SAM_decarbox结构域的蛋白; GmSAMDC1基因所编码蛋白的理论等电点为4. 86,相对分子质量为38 987. 13 Da,为亲水性蛋白,不含跨膜区;大豆中共有8个GmSAMDC1的同源基因,GmSAMDC1与红车轴草(PNY09439. 1,82. 64%)和拟南芥(At SAMDC3,67. 22%)中的SAMDC蛋白编码基因亲缘关系最近; GmSAMDC1启动子序列包含防卫和胁迫响应元件、植物激素应答元件、光应答元件等许多顺式作用元件; GmSAMDC1在花中的表达量最高,与其大豆同源基因Glyma. 01G071300(序列相似性为94. 6%)、Glyma. 18G278800(66. 8%)和Glyma. 08G25580(66. 5%)的组织特异性表达模式比较相似; Glyma. 02G128000-GFP融合蛋白在细胞膜和细胞质上表达。本研究结果为进一步阐明大豆GmSAMDC1基因在大豆耐胁迫过程中的作用提供了理论依据。     

基于大豆导入系群体结瘤相关性状QTL定位及基因挖掘

氮元素是植物生长所需的重要营养元素之一,生物固氮是大豆生长所需氮元素的重要来源。本实验室以绥农14为母本,野生豆ZYD00006为父本构建一套覆盖野生大豆全基因组的导入系群体。基于大豆导入系群体对结瘤数目和干重进行QTL定位,定位到3个QTL与根瘤干重相关,分布在N、M和D2这3个连锁群上,定位到5个QTL与根瘤数相关,分布在K、F、J和D2这4个连锁群上,在D2连锁群这两个性状有重叠区段(7.20-7.79Mb)。针对这个重叠区段的63个基因进行基因注释,选择到6个与共生、抗病相关的基因作为候选基因进行下一步验证。qRT-PCR分析表明根瘤菌侵染期间Glyma.17G097000基因表达模式与对照相比差异很大。Glyma.17G097000属于GmHIR基因家族,在大豆基因组中发现了11个家族成员。GmHIR家族基因结构相似性很高,基因表达有组织特异性。大豆HIR蛋白有Stomatins和Prohibitin两个结构域,能够参与离子通道调节等生理过程,与植物抗病和细胞周期有关。GmHIR基因来源于四个祖先,进化过程中是高度保守的。对GmHIR基因家族11个成员qRT-PCR检测,结果显示根瘤菌感染期间Glyma.05G029800、Glyma.09G154400和Glyma.17G097000这三个基因与对照相比表达模式差异较大。结果表明这三个基因可能在大豆共生体系建立中起着重要的作用,参与根瘤菌与大豆共生体系建立过程中的离子通道调节和免疫反应。本研究为大豆-根瘤菌共生机制研究奠定基础并提供有效候选基因。     

大豆胞囊线虫侵染下野生大豆根组织实时定量PCR分析中内参基因的筛选

以高抗和高感HG 0大豆胞囊线虫的野生大豆种质为试验材料,采用q RT-PCR技术检测了ACT11(Glyma.18G290800)、Tubulin-motif(Glyma.19G194800)等24个候选持家基因在不同大豆胞囊线虫侵染时期(接种后9,15和20 d)野生大豆根系中的基因表达情况,除Actin(Glyma.08G182200)出现了非特异扩增外,其它23个持家基因的q RT-PCR扩增效果均理想,扩增特异性高。分析上述23个持家基因的表达丰度,并利用Best Keeper、ge Norm和Normfinder对其表达稳定性进行评价。结果表明:持家基因间的表达丰度不尽相同,而且有些持家基因在不同品种间或不同处理间(接种虫卵和接种水)也存在表达丰度的差异。23个持家基因的表达稳定性也不相同,综合来看,表达最不稳定的基因为Cons9(Glyma.10G152200)、Cons2(Glyma.17G138500)、Cons1(Glyma.15G270900)和Tubulin-motif(Glyma.20G136000),本试验条件下它们不适宜作内参基因;其余持家基因相对稳定,本试验条件下可选择Tubulinmotif(Glyma.15G132200),Cons6/SKIP16(Glyma.12G051100)或Tubulin-motif(Glyma.05G110200)作为内参基因,它们表达量较高,表达最为稳定,其Ct平均值分别为25.7、26.5和25.0,标准偏差分别为0.540、0.575和0.490,ge Norm软件评价其稳定性的M值分别为0.698、0.715和0.727。该研究为野生大豆抗胞囊线虫相关基因的表达分析提供了参考。     

大豆Glyma.18G261700基因生物学分析及功能初探

MYB是参与花青素合成的重要调控基因，数量众多且功能多样。为分析大豆中与花生黑色种皮调控基因AhTc1(MYB家族)同源基因的功能，本研究根据AhTc1基因序列信息，同源克隆大豆MYB家族基因Glyma.18G261700,对其进行生物信息学分析，并检测不同种皮颜色大豆苗期不同部位花青素和基因表达量。结果显示：Glyma.18G261700全长1 572 bp,编码189个氨基酸，属于R2R3型MYB。AhTc1编码的蛋白为疏水性蛋白，而Glyma.18G261700与之不同，编码的蛋白为亲水性蛋白。三级结构中Glyma.18G261700与AhTc1尾部有所区别，在调控种皮颜色中可能具有不同功能。不同种皮颜色大豆幼苗根部花青素含量较低，而Glyma.18G261700基因的表达量较高，表明在大豆幼苗期该基因的大量表达可能抑制根部花青素的合成，但对种皮颜色的调控还需进一步探索。     

不同播期条件下大豆开花期性状的全基因组关联分析

开花期是决定大豆生态适应性和产量的重要性状,研究大豆开花期的遗传规律,解析大豆开花调控的分子基础,对大豆品种适应性改良具有重要意义。本研究选用来自北方、黄淮海和南方三大生态区的87份中国品种及31份美国大豆品种,在北京进行春、夏季分期播种,获得10个播期环境下的开花期数据,并利用Illumina BARCSoySNP6K iSelectBeadChip芯片对大豆品种进行SNP分型。通过群体结构分析、聚类分析和PCA分析,将供试大豆品种分为两个亚群。利用TASSEL软件PCA+K模型,结合5 126个SNP标记,对开花期进行全基因关联分析,得到18个显著关联位点,其中14个位点在2个以上的环境中稳定表达,4个位点只在单一环境中被检测到,表明其受光温环境影响较大。通过对显著关联位点进行预测,共发现31个候选基因,其中有2个（Glyma10g36600和Glyma19g37890）已在此前的研究中被证明与大豆开花有关,而Glyma03g33470、Glyma11g15650、Glyma05g01510和Glyma19g07162与拟南芥中相关的开花调控基因同源,推测这4个基因在大豆中也参与开花调控。11个尚未见报道的新位点分别位于3、4、5、9、10、14、16和17号染色体上,可能存在控制大豆开花期的新基因。     

大豆水溶性蛋白质的全基因组关联分析

为进一步解析中国大豆种质水溶性蛋白的遗传基础,为大豆高水溶性蛋白质的分子标记辅助选择育种及品质改良提供理论依据,本研究以224份大豆种质为试验材料,于2017和2018年对大豆水溶性蛋白质含量进行测定,利用1 514个高质量的SNP标记分别对2017、2018年水溶性蛋白质含量及两年均值进行全基因组关联分析,共检测到18个显著关联的SNP标记,这些SNP标记涉及16个位点,有8个位点至少被检测到2次,其余8个位点仅被检测到1次,表明其受环境因素影响较大。16个位点中有7个尚未见报道,分别位于8、11、13、14和15号染色体上,是新发现的控制大豆水溶性蛋白的位点。对表型变异解释率较高且稳定关联的2个位点qWSPC7和qWSPC8-1候选区间内的基因进行预测,共获得25个候选基因,其中有7个基因(Glyma.07g195000、Glyma.08g103100、Glyma.08g108900、Glyma.08g105100、Glyma.08g107800、Glyma.08g107700和Glyma.08G115800)在大豆籽粒、根或根瘤中具有较高的表达水平。这些基因可作为水溶性蛋白质的候选基因,可能具有调控大豆水溶性蛋白质的功能。     

耐低钾大豆品种筛选及低钾胁迫下Lee 68的差异表达基因分析

为了从分子水平上研究大豆幼苗对低钾胁迫的耐性机理,本研究首先依据生物量指标从25份大豆材料中筛选出耐低钾品种Lee 68,并对低钾胁迫下大豆Lee 68幼苗进行转录组测序和分析。通过对转录组Solexa/Illumina高通量测序数据的分析,共得到160 211 759个reads,将获得的数据与大豆Williams 82基因组序列比对,比对率达到92.83%以上。比较组LK_VS_CK中差异表达基因为3 521个,其中下调和上调基因分别为2 393和1 128个。GO功能聚类分析显示LK_VS_CK比较组的差异表达基因主要富集于植物的代谢过程、胁迫响应及信号转导;KEGG pathway分析中,LK_VS_CK比较组有390个差异表达基因显著富集在19个pathway途径中,其中富集最为显著的是代谢途径,包括氨基酸代谢、脂肪酸和类脂代谢以及碳水化合物代谢等;COG分类统计结果表明LK_VS_CK比较组除649个差异表达基因具有一般性功能之外,377个基因被划分到转录因子功能类,343个基因被划分到信号转导机制功能中。结合差异表达基因的功能分析及茉莉酸信号转导调控机制,筛选到了茉莉酸信号途径上可能参与钾离子吸收转运过程的8个关键候选基因,分别是Glyma09g08290、Glyma09g33730、Glyma03g32890、Glyma04g39010、Glyma 08g09720、Glyma12g36310、Glyma16g32821、Glyma19g45260。研究结果对深入研究大豆钾胁迫下与钾离子高效吸收相关基因的调控及克隆研究具有重要的参考价值。     

大豆RING/U-box蛋白Glyma.13G115900的克隆及其对非生物胁迫的应答

课题组前期对干旱胁迫下大豆转录组的数据分析发现大豆Glyma.13G115900基因编码一个RING/U-box蛋白,其表达水平受干旱胁迫影响显著。本研究以垦丰16大豆为试验材料,克隆Glyma.13G115900基因。氨基酸多重序列比对表明其编码的蛋白与其它物种都具有高度保守的RING/U-box结构域。构建原核表达载体pET-29b-Glyma.13G115900转化到大肠杆菌中,Glyma.13G115900蛋白在大肠杆菌中能够表达。荧光定量PCR分析表明Glyma.13G115900基因的表达量受PEG、NaCl和ABA的影响显著,但基本不受冷胁迫的诱导。经PEG和NaCl处理后,该基因表达量与CK相比呈现出显著下降的趋势,PEG处理的表达量变化比NaCl下调的更明显;在ABA诱导下与CK相比该基因的mRNA丰度呈现出先上升后下降的趋势,在4 h表达量出现峰值,推测该基因可能通过依赖于ABA途径参与非生物胁迫应答。以上结果为进一步深入研究该基因的调控途径奠定基础。     

大豆TOE基因的进化分析及开花调控功能解析

为了进一步解析大豆中重要植物开花和花器官发育调控转录因子AP2的编码基因TOE的进化规律及其对开花功能的调控作用,为大豆TOE基因的功能解析和大豆纬度适应性研究提供基础,本研究利用生物信息学手段对大豆TOE基因进行聚类分析、序列特征分析、染色体区段共线性分析和组织特异性表达分析,预测关键开花基因启动子区段AP2结合位点,并验证不同单倍型大豆开花时间.结果 显示:从PlantTFDB数据库检索到12个大豆TOE基因,GmTOE6b(Glyma.02G087400)为新发现的大豆TOE基因.GmTOE6a和GmTOE6b均只有1个AP2结构域,其余大豆TOE基因均有两个AP2结构域.6个大豆TOE基因与拟南芥TOE1基因聚为一类;2个与拟南芥TOE2基因聚为一类;4个与拟南芥TOE3、AP2聚为一类.12个大豆TOE基因都有且仅有1个miR172靶位点,且该靶位点序列与拟南芥TOE的miR172靶位点序列高度一致.染色体区段共线性分析显示,大豆12个TOE基因按起源方式可以分为3类,6个随大豆基因组复制而产生;4个起源于大豆物种形成之前且与拟南芥TOE基因有共同祖先;2个起源于大豆物种形成之前且与拟南芥TOE基因无共同祖先.大豆开花关键基因GmFT2a和GmFT5a启动子序列中均含有多个AP2结合位点,GmTOE4b和GmTOE5b两个基因均可影响大豆的开花时间.研究结果说明大豆中12个TOE基因极有可能均是miR172的靶基因,虽然其编码蛋白的氨基酸组成非常相似,但它们的进化规律和组织特异性表达规律存在不同,它们在进化过程中可能存在功能分化.GmTOE4b和GmTOE5b可能通过结合GmFT2a和GmFT5a启动子上的顺式元件来调控其基因转录,从而调控开花.     

大豆耐盐相关基因研究进展

大豆（Glycine max（L.）Merill）是植物蛋白质和油脂的重要来源。盐胁迫是造成大豆产量损失的主要非 生物胁迫因素。耐盐基因的挖掘对培育大豆耐盐品种至关重要。本文一方面总结了通过正向遗传学获得的大豆 耐盐相关数量性状位点或基因,如萌发期耐盐性主效基因GmCDF1（Glyma.08g102000）、2个出苗期QTL（分别位于6 号和14号染色体）;苗期耐盐性主效基因GmSALT3（Glyma03g32900）以及位于G连锁群的QTL。随着对大豆耐盐性 研究的不断深入,目前认为大豆萌发期、出苗期、苗期的耐盐性无直接相关性。另一方面总结了通过反向遗传学途 径获得的参与离子运输、转录调控等耐盐性基因,以及通过生物工程技术转入外源基因提高大豆耐盐性的研究进 展,期望为解析大豆耐盐分子机制和耐盐育种提供参考。     

多环境下大豆子粒大小性状的全基因组关联分析

为探究多环境下大豆子粒大小性状的分子遗传基础,挖掘与子粒大小性状相关的SNP位点和候选基因,利用150份大豆种质资源在2019年和2020年6个环境条件下对大豆子粒粒长、粒宽、粒厚和百粒重性状进行表型测定,并进行全基因组关联分析。结果表明：在CMLM（压缩混合线性）模型下,在6个环境条件下检测到896个与子粒大小性状显著关联的SNP位点,分布于20条染色体。不同性状检测到72个重叠的SNP位点。检测到39个稳定遗传的SNP位点,贡献率为10.68%~24.93%。通过稳定性与重叠性分析,获得35个稳定表达的SNP位点,贡献率为10.92%~23.16%。在粒宽、粒厚及百粒重性状中同时检测到显著关联的SNP位点最多,位点rs16533609的贡献率最高（16.51%）。根据稳定表达的SNP筛选候选基因,推测Glyma.03G006600、Glyma.04G077100、Glyma.08G203600、Glyma.12G195400、Glyma.17G039800、Glyma.18G202100和Glyma.20G215700等7个基因对大豆子粒大小性状有调控作用。     

大豆Glyma.05G222700.2基因生物信息学与逆境表达分析

为研究Glyma.05G222700.2基因编码的丝氨酸/苏氨酸蛋白激酶在抗非生物胁迫过程的功能和原理,促进大豆抗逆候选基因的开发利用,本研究通过生物信息学方法对大豆Glyma.05G222700.2基因进行同源序列、蛋白结构、进化树和转录组分析,通过qRT-PCR分析盐胁迫下大豆不同组织中该基因的表达情况。结果表明:该基因编码区长2 040 bp,编码697个氨基酸,预测分子量为656.54 kD,pI8.026。多序列比对发现Glyma.05G222700.2蛋白包含1个Pkinase结构域。进化树分析表明该蛋白与野大豆、刺毛黧豆、赤豆一致性较高。转录组数据表明Glyma.05G222700.2基因在大豆各组织中均有表达,其中在种子中表达量最高,在根中表达量最低。qRT-PCR结果发现Glyma.05G222700.2基因在毛状根、茎、叶中均有表达,在茎中表达量最高,在叶中表达量最低;在毛状根中盐胁迫12 h表达量达到极值,盐胁迫24 h表达量下降;在茎中表达量呈上升趋势,在24 h表达量达到极值;在叶中表达量不稳定,盐胁迫2 h该基因不表达,盐胁迫6 h该基因表达量达到极值并高于对照,盐胁迫12和24 h表达量低于对照。推断该基因可能在大豆抵抗盐胁迫过程中起到重要作用。     

大豆GmIDD及其同源基因的生物信息学分析及功能预测

利用NCBI和Phytozome数据库进行GmIDD基因的同源基因检索,获取该基因在大豆不同染色体中5个拷贝的序列信息。通过对功能结构域的分析结果显示GmIDD基因5个拷贝均含有2个C2H2型(72~92 aa、114~142 aa)、2个C2HC型(149~169 aa、176~195 aa)锌指结构域。GmIDD基因cDNA序列全长为1 227 bp,编码408个氨基酸。系统进化树分析结果显示GmIDD基因5个拷贝中除了Glyma.14g10940和Glyma.17G228000外,其它IDD蛋白间进化距离较远,表明不同大豆IDD基因在功能上可能存在一定差异。长短日照处理下GmIDD表达量的分析结果显示GmIDD基因受短日照诱导表达,因此推测GmIDD可能参与大豆开花的调控过程。     

大豆病程相关蛋白PR-5及其同源蛋白TPLs的生物信息学分析

大豆病程相关蛋白PR-5又称为类甜蛋白TLPs,广泛分布于高等植物体内,各种生物和非生物胁迫诱导(病原微生物、渗透胁迫,机械损伤和植物激素等)均可诱导其表达。本研究以PR-5编码基因Glyma01g42670.1为目的基因,对该基因及其同源蛋白的结构和功能进行了分析;通过生物信息学的方法,共筛选获得42条大豆TLPs蛋白质序列,并用Mega 4.1软件构建其系统进化树;同时对Soy Base中的大豆TLPs相关的表达数据进行分析。结果表明:大豆病程相关蛋白PR-5的编码基因Glyma01g42670.1编码240个氨基酸的多肽链,p I 6.26,是一种酸性蛋白,属于疏水蛋白,该蛋白质具有典型的TLP家族的保守结构域,具有指导PR-5蛋白的跨膜转移(定位)的N端信号肽序列,PR-5蛋白不含有跨膜螺旋结构,主要存在于质外体中,属于分泌蛋白,推测该蛋白具有内切葡聚糖酶的催化活性,进化分析表明:该编码基因与拟南芥的NP_192902.1基因属于直系同源基因,它们是具有相同功能和共同起源的基因,组织表达特异性分析表明该基因主要在根部和花中表达。     

基于共线性分析的大豆种子硬实性QTL定位及候选基因预测

大豆种子硬实度关乎食品与饲料加工。为提高大豆种子硬实度分子标记辅助选择的效率，发现相关候选基因，利用染色体片段代换系（chromosome segment substitution lines, CSSL）对大豆种子硬实性进行了两年的QTL定位研究，结合前人得到的25个种子硬实性QTL，利用MCScanX对整个大豆基因组进行分析，生成共线性区组，评估大豆种子硬实性QTL之间的共线性，确定了位于Gm02片段的中心QTL。利用多个数据库分析hub-QTL区段的基因，锚定了8个与大豆种子硬实性相关的候选基因。从CSSL群体中选择种子硬实性不同的两个品系和轮回亲本，用于随后的实时荧光定量PCR分析，发现候选基因Glyma.02G269400和Glyma.02G269500在CSSL群体中硬实性不同的2个品系及轮回亲本绥农14中的表达差异显著。Glyma.02G262600在绥农14中的相对表达量约是CSSL-136的5倍，而在CSSL-200中表达量中等，推测该基因抑制大豆种皮的形成。     

大豆橙花叔醇合酶基因克隆及功能鉴定

通过同源检索得到大豆橙花叔醇合酶(nerolidol synthase,Gm NES)基因Glyma13g32380,采用RT-PCR方法从经水杨酸(SA)处理的大豆材料中,克隆到c DNA片段,其包括完整的ORF(1 605bp)。共编码534个氨基酸,蛋白分子量约62.338 k D。该蛋白序列与番茄等植物的NES序列同源性较高,包含DDXXD保守结构域。将该ORF构建入p ET28a表达载体,在大肠杆菌BL21(DE3)中表达为可溶性蛋白。通过与法尼基焦磷酸合酶基因在大肠杆菌中进行微生物代谢工程共表达,气质联用检测到橙花叔醇的生成。同时还检测到SA模拟生物胁迫处理诱导大豆叶片产生橙花叔醇,并且在1.0 mmol/L SA处理下,Gm NES基因表达量明显上调。本研究克隆了大豆橙花叔醇合酶基因,并鉴定了其生化功能及其在植物体内的生理相关性,为研究大豆萜类代谢及对病虫害的防御反应和抗性品种选育提供了理论依据。     

利用ISSR标记进行菜用大豆和普通大豆的分类

ISSR标记技术因具有很高的多肽性而被广泛应用于遗传多样性研究.进行菜用大豆品种多样性研究将有助于育种家进行杂交亲本的合理选配.采用ISSR标记技术对37个菜用大豆品种(系)和13个普通大豆品种(系)进行了分类研究.结果表明:50个ISSR引物中有11个引物没有多肽性,其他39个引物共产生132条谱带,其中81条多肽性谱带,占61.4%.引物扩增的谱带数在l～7条之间,大小为830～3530 bp.每个引物平均扩增的谱带数和多肽性带数分别为3.3和2.1.聚类分析结果表明,可将供试的50个品种划分为3个类组,扁茎大豆Taikadai-zu单独一组,7个普通大豆(多数为中国选育)划归为B组,C组由日本菜用大豆和4个普通大豆组成.C组还可以划分成若干亚组,同一亚组的品种具有相似的性状,比如有一亚组的品种均是褐脐、茶色种皮、白花、不易炸荚.系选的品种和原始品种在分类图中也显示出其密切的亲缘关系.基于ISSR标记的分类可以反映菜用大豆的遗传关系,菜用大豆和普通大豆之间存在遗传多样性,合理保存菜用大豆种质将有助于未来大豆品质改良计划.     

大豆“冀黄13”突变体筛选及突变体库的建立

利用甲基磺酸乙酯(ethyl methane sulfonate,EMS)、叠氮化钠(sodium azide,NaN3)和N离子束分别诱变处理大豆品种"冀黄13"的种子。经M2选择,M3、M4代验证,共筛选出茎器官突变体15份,叶器官突变体19份,花器官突变体31份,种子器官突变体3份以及成熟期突变体18份,其中有9份突变体属于复合性状的突变,最终获得77份稳定遗传的突变体。通过F2群体的表型分析,确定1份雄性不育突变体和1份黑种皮色突变体的表型均由单隐性基因控制。本研究所获得的突变体和所构建突变体库将为大豆功能基因组学研究和遗传改良奠定基础。