大豆单株荚数QTL定位及整合

大豆单株荚数是影响大豆产量的一个重要农艺性状，同时也是育种中的一个重要指标。在大豆产量构成因素(百粒质量、单株荚数、单株粒数等)中，单株荚数与产量相关性较高[1]。同时研究还发现，单株荚数与百粒质量、每荚粒数、主茎节数等重要性状关联度和相关系数较高[2]，单株荚粒数的提高对百粒质量等农艺性状不会有较大的负面影响[3]，说明通过选育较高单株荚数的材料来提高大豆产量是可行的。因此，可以通过选择较多单株荚数的材料进行聚合育种，从而更有效地提高大豆产量。

单株荚数是一个非常复杂的性状，它是由多因素调控的，其中最关键的3个过程分别为花芽分化数、胚珠成功受精比率和成功受精后的胚胎正常发育为荚果的比率[4]。花芽的分化和发育决定子房胚珠数和可育胚珠比率。成功受精胚珠比率是由受精过程所决定的，其中包括花粉育性、花粉与柱头的接触量、花粉粒的萌发、受精条件以及花粉管的导入等过程。受精后的胚珠最终是否能发育为种子由结合子发育过程决定。然而其中任何一个生物学过程都是相当复杂的，因此,将通过研究每一个过程的分子机理和功能，才能最终将单株荚数这个复杂的性状解析开来。

单株荚数是一个典型的数量性状，受多基因控制，同时也易受环境的影响。目前，就国内外已有的关于大豆荚数QTL定位研究的报道可知[5-16]，定位了约50多个QTLs位点，这些位点几乎分布于大豆20条染色体上(除B2和H)；解释表型变异平均值为13.4%，约有69%的QTLs解释表型变异在20%之下，主效荚数QTL也较少。

本研究利用多荚材料C025和少荚材料JD18为亲本，然后利用双亲构建的F2群体定位大豆单株荚数QTL，通过定位的大豆荚数QTL，阐明大豆单株荚数的分子机制；同时整合前人对大豆荚数定位的QTL，解析大豆单株荚数QTL在全基因组的分布以及联系，旨在为选育大豆较高单株荚数提供材料基础并为后续研究大豆单株荚数分子机制奠定基础。

1 材料和方法

1.1 试验材料

本研究使用的亲本材料C025(多荚)和JD18(少荚)来源于山西省农业科学院农作物品种资源研究所，双亲杂交得到182个单株F2群体。2016年5-10月在太原山西省农业科学院东阳基地完成杂交试验得到F1杂交种，2017年5-10月和2018年5-10月(太原2017年、太原2018年)连续2 a在太原山西省农业科学院东阳基地完成亲本和F2群体田间种植和表型鉴定。

1.2 试验方法

田间试验按照完全随机区组设计，3次重复，每个小区2行，行距(株距为50.0 cm×13.5 cm。大豆成熟时，每个小区选择10个代表单株人工统计荚数。具体调查标准参照《中国大豆品种志》[17]进行。

1.3 QTL定位分析

利用在公共数据库中已发表的大豆遗传连锁图谱(https://soybase.org/)，结合荚数表型鉴定进行QTL定位。QTL扫描采用WinQTL Cartographer 2.5软件(http://statgen.ncsu.edu/qtlcart/WQTLCart.htm)中的复合区间作图法[18]。LOD阈值显著性采用1 000次排布测验分析[19]，基本参数步长、窗口大小、扫描间距和运行速度分别设置为1，10，5，1 cM，P=0.05显著水平检测并确认QTL。采用QTL-元分析方法整合前人研究的不同环境不同群体检测的大豆荚数QTL[20]。

1.4 数据分析

利用SAS V8软件中的PROC ANOVA程序对不同环境单株荚数表型进行方差分析，从而估算遗传力。群体单株荚数频率分布和亲本间单株荚数的显著性分析采用Excel中的函数进行。

2 结果与分析

2.1 父母本及F2群体单株荚数表型变异和遗传力分析

亲本C025和JD18在2个不同环境中的单株荚数均表现出极显著差异(P<0.01)，C025单株总荚数明显高于JD18，是JD18的近2倍，并且在2 a表型考察中较稳定(表1)。F2分离群体单株荚数在不同环境中都表现出广泛的变异(2～150)，并且有些株系有明显的超亲分离，表明控制大豆荚数的基因在2个亲本中都存在。F2分离群体单株荚数在2个环境中都呈现近似正态分布，可以说明单株荚数符合数量性状特征，适合利用QTL定位分析(图1)。通过方差分析评估大豆荚数的遗传力为65%，说明大豆荚数遗传力较大，有利于进行后期的遗传改良。

2.2 大豆单株荚数QTL定位分析

利用1张包括1 015个SSR标记的大豆遗传图谱(https://soybase.org/)，其覆盖大豆基因组2 282.92 cM，长度为 71.31～165.73 cM，相邻标记间的平均距离为2.44 cM，这些SSR标记在 20 条染色体上分布较均匀。利用软件WinQTL Cartographer 2.5对2个环境的单株荚数分别进行全基因组QTL扫描，结果显示，单环境下共检测出33个QTLs，分别位于D1a、N、C1、C2、M、A2、K、O、B1、F、B2、E、J、D2、G、L和I 共17个染色体上，LOD值在2.8～15.1，贡献率在0.2%～56.4%，平均值为14.9%。经过元分析整合后，共定位23个QTLs，分别位于D1a、N、C1、C2、M、A2、K、O、B1、F、B2、E、J、D2、G、L和I 共17个染色体(表2)。其中，有10个QTLs在2个环境中被重复检测到，分别定位于A2、D2、C2、I、B1、C1和J染色体上，LOD值在2.8～15.1，平均值为9.4，解释的平均表型变异0.2%～56.4%，平均值为25.2%，其中有6个QTLs的贡献率高于20%。qPN.C2-4、qPN.C1、qPN.J、qPN.D2-1和qPN.A2-1的贡献率分别为29.5%，35.4%，38.9%，53.9%和56.4%，均高于前人定位的大豆荚数QTLs。通过比较发现，qPN.C2-4和qPN.C1的置信区间与前人的研究重叠，并且本研究的效应值高于前人，同时QTL-qPN.A2-1的贡献率为56.4%，是目前检测到的效应值最大的控制大豆单株荚数的QTL。因此，这3个QTLs可看作控制大豆单株荚数的主效QTL，并对其进一步深入解析。

2.3 大豆荚数QTL整合分析

通过查阅已报道的与大豆荚数定位相关的文献，结果显示，共定位了51个QTLs，在大豆18个染色体上均有分布(在B2和H中没有检测到荚数QTL)，其中，在C2染色体上并且在一个区段内分布较为集中。经过元分析整合前人定位的大豆荚数QTLs，最终确定39个QTLs，LOD值为2.80～6.83，贡献率为6.54%～27.01%，平均值为13.4%(表3)，约有69%的QTLs的贡献率在20%之下，约有5个QTLs可看作主效QTLs。整合的39个QTLs，其中有6个重复检测到的QTLs被整合，分别分布在B1、C1、C2、D1a和L染色体上。

与前人研究的结果相比，本研究新定位了18个大豆单株荚数QTLs，并且首次在B2染色体上定位了1个控制大豆单株荚数的QTL(qPN.B2)，充分补充了大豆单株荚数分子机制研究。同时还发现，本研究定位的23个QTLs中有5个QTLs与前人的研究重叠，分别为qPN.C2-3、qPN.I、qPN.C2-4、qPN.C1和qPN.L，其中，qPN.C2-4和qPN.C1是本研究发现的主效QTL，说明qPN.C2-4和qPN.C1是可多环境多材料重复检测的大豆单株荚数主效QTLs，这可为后续大豆单株荚数分子研究奠定基础。

3 结论与讨论

由于前人对大豆荚数QTL定位研究所用的遗传图谱和标记不一致，以至于目前还未得到一个完整的大豆荚数分子遗传图谱。本研究基于已有研究的大豆遗传图谱，将大豆荚数QTL整合在同一张遗传图谱上。这些QTLs几乎分布于大豆20条染色体上(除H)，其绝大部分的贡献率在20%以下，只有几个达到20%以上[7]，这几个QTLs可为后续大豆单株荚数QTL精细定位和基因克隆奠定基础。

本研究共定位了23个大豆单株荚数QTLs，其中有5个与已报道的QTL重叠，另外18个代表新的大豆单株荚数QTL定位，填补了在B2染色体上没有检测到大豆单株荚数QTL的空白(本研究定位了qPN.B2)。qPN.C2-4和qPN.C1虽然与前人研究重叠[5, 7]，但是本研究的主效QTL贡献率高于前人报道，说明QTL-qPN.C2-4和QTL-qPN.C1是控制大豆单株荚数的主效QTLs位点，可被不同材料和不同环境重复检测，可作为后续精细定位和基因克隆的首选QTLs。同时本研究多环境检测到qPN.A2-1的贡献率为56.4%，是目前检测到的效应值最大的QTL，并且与前人研究没有重叠。因此可以看作是一个新的控制大豆单株荚数的QTL。

整合前人对大豆单株荚数的研究发现，目前仅处于一个初步研究阶段，还未阐明大豆单株荚数分子机理。Wang等[21]在大豆中克隆GmCYP78A10基因，通过分析发现，GmCYP78A10基因的表达影响大豆单株荚数。目前，还未有其他有关基因在大豆中报道，这可能与大豆单株荚数性状复杂并且受环境条件影响较大难以精准鉴定有关。

总之，本研究定位的3个主效QTLs贡献率均较大，完全具备进一步精细定位和基因克隆的基础，从而可为大豆单株荚数遗传改良打下坚实的理论基础和技术服务。

[1] 王秋玲，郭凌云，刘艳，张桂花，刘凤洲. 夏大豆单株产量与主要农艺性状的灰色关联度分析 [J]. 大豆通报，2002(6): 5,12.doi:10.3969/j.issn.1674-3547.2002.06.005.

Wang Q L, Guo L Y, Liu Y, Zhang G H, Liu F Z.Analysis of grey related degree between individual plant yield and main agriculture character of summer soybean [J]. Soybean Bulletin, 2002(6):5,12.

[2] 王芳，李传仁. 大豆单株荚数与主要农艺性状关系的分析 [J]. 种业导刊，2008(11): 24-25,27.doi: 10.3969/j.issn.1003-4749.2008.11.006.

Wang F, Li C R. The analysis of the relationship between pod number per plant and the main agronomic traits in soybean [J]. Journal of Seed Industry Guide, 2008(11): 24-25,27.

[3] 戴瓯和. 夏大豆荚粒性状遗传变异、遗传相关和选择指数研究的初报 [J]. 安徽农业科学， 1982(2): 46-50.doi: 10.13989/j.cnki.0517-6611.1982.02.010.

Dai O H. The report of genetic variation, genetic correlation and selection index research of seed and pod in summer soybean [J]. Journal of Anhui Agricultural Sciences, 1982(2): 46-50.

[4] Qu B, Zhang W, Chen X H, Li N, Cui N, Li T L. Research progress of flower bud differentiation mechanism of plant [J]. Chinese Agricultural Science Bulletin, 2010, 26(24): 109-114.

[5] Zhang D，Cheng H，Wang H，Zhang H Y，Liu C Y，Yu D Y. Identification of genomic regions determining flower and pod numbers development in soybean(Glycine max L.)[J]. Journal of Genetics and Genomics, 2010, 37(8): 545-556. doi:10.1016/S1673-8527(09)60074-6.

[6] Chen Q S, Zhang Z C, Liu C Y, Xin D W, Qiu H M, Shan D P,Shan C Y, Hu G H. QTL analysis of major agronomic traits in soybean [J]. Agricultural Sciences in China, 2007, 6(4): 399-405. doi:10.1016/S1671-2927(07)60062-5.

[7] Sun D S, Li W B, Zhang Z C, Chen Q S, Ning H L, Qiu L J, Sun G L. Quantitative trait loci analysis for the developmental behavior of Soybean(Glycinemax L. Merr.)[J]. Theoretical and Applied Genetics, 2006, 112(4): 665-673.doi:10.1007/s00122-005-0169-y.

[8] Liu W X, Moon Young Kim, Van K, Lee Y H, Li H L, Liu X H, Lee S H. QTL identification of yield-related traits and their association with flowering and maturity in soybean [J]. Journal of Crop Science and Biotechnology, 2011, 14(1): 65-70. doi:10.1007/s12892-010-0115-7.

[9] Palomeque L, Liu L J, Li W B, Hedges B, Cober E R,Rajcan I. QTL in mega-environments:Ⅱ. agronomic trait QTL co-localized with seed yield QTL detected in a population derived from a cross of high-yielding adapted×high-yielding exotic soybean lines[J]. Theoretical and Applied Genetics, 2009, 119(3):429-436. doi:10.1007/s00122-009-1048-8.

[10] Yosuke Kuroda, Akito Kaga, Norihiko Tomooka, Hiroshi Yano, Yoshitake Takada,Shin Kato, Duncan Vaughan. QTL affecting fitness of hybrids between wild and cultivated soybeans in experimental fields [J]. Ecology and Evolution, 2013, 3(7): 2150-2168. doi:10.1002/ece3.606.

[11] Sangrea Shim, Moon Young Kim, Jungmin Ha, Yeong-Ho Lee, Suk-Ha Lee. Identification of QTLs for branching in soybean(Glycine max(L.)Merrill)[J]. Euphytica, 2017, 213(9): 225. doi:10.1007/s10681-017-2016-z.

[12] Bobby R, Bazzelle R, Clark W, Kantartzi S, Meksem K, Akond M, Kassem M A. Genetic analysis of yield components and seed isoflavones contents in the PI 438489B by "Hamilton" recombinant inbred line(RIL)population of soybean(Glycine max(L.)Merr)[J]. Journal of Agricultural Science, 2012, 4(9): 697-702. doi:10.5539/jas.v4n9p98.

[13] 李莹莹，李瑞超，程春光，赵圆圆，刘春燕，齐照明，李灿东，王囡囡，蒋洪蔚，陈庆山. 大豆荚粒数相关QTL的Meta和Overview分析及其候选基因预测 [J]. 农业生物技术学报，2018，26(11): 1821-1833. doi: 10.3969/j.issn.1674-7968.2018.11.001.

Li Y Y, Li R C, Cheng C G, Zhao Y Y, Liu C Y, Qi Z M, Li C D, Wang N N, Jiang H W, Chen Q S. Meta and overview analysis of QTL associated with pod and seed traits and candidate gene mining in Soybean(Glycine max)[J]. Journal of Agricultural Biotechnology, 2018, 26(11): 1821-1833.

[14] 杨胜先，牛远，李梦，魏世平，刘晓芬，吕海燕，章元明. 栽培大豆农艺性状的关联分析及优异等位变异挖掘[J].中国农业科学，2014，47(20):3941-3952. doi:10.3864/j.issn. 0578-1752.2014.20.002.

Yang S X, Niu Y, Li M, Wei S P, Liu X F, L H Y, Zhang Y M. Association mapping of agronomic traits in soybean(Glycine max L. Merr.)and mining of novel alleles[J]. Scientia Agricultura Sinica, 2014, 47(20):3941-3952.

[15] 徐琰，孙晓环，孙霞，王燕平，宗春美，齐玉鑫，白艳凤，任海洋，潘相文，杜维广，孔凡江，刘宝辉.大豆花荚脱落及单株荚数的QTL定位 [J]. 土壤与作物，2015，4(2): 71-76. doi:10.11689/j.issn.2095-2961.2015.02.004.

Xu Y, Sun X H, Sun X, Wang Y P, Zong C M, Qi Y X, Bai Y F, Ren H Y, Pan X W, Du W G, Kong F J, Liu B H. QTL mapping of flower and pod abscission and pod number per plant in soybean [J]. Soil and Crop, 2015, 4(2): 71-76.

[16] Orf J H, Chaseb K, Jarvikb T, Mansurc L M, Cregand P B, Adlerb F R, Lark K G. Genetics of soybean agronomic traits: I. Comparison of three related recombinant inbred populations [J]. Crop Science, 1999, 39(6): 1642-1651.doi: 10.2135/cropsci1999.3961642x.

[17] 吉林省农业科学院大豆研究所. 中国大豆品种志(1978-1992)[M]. 北京: 农业出版社，1993.

Soybean Institute of Jilin Agricultural Sciences. Chinese soybean cultivar catologue(1978-1992)[M]. Beijing: Agriculture Press, 1993.

[18] Zeng Z B. Precision mapping of quantitative trait loci [J]. Genetics, 1994, 136(4): 1457-1468.

[19] Churchill G A，Doerge R W.Empirical threshold values for quantitative trait mapping[J]. Genetics, 1994, 138(3): 963-971.

[20] Goffinet B, Gerber S. Quantitative trait loci: a meta-analysis [J]. Genetics, 2000, 155(1): 463-473.

[21] Wang X B, Li Y H, Zhang H W, Sun G L, Zhang W M, Qiu L J. Evolution and association analysis of GmCYP78A10 gene with seed size/weight and pod number in soybean [J]. Molecular Biology Reports, 2015, 42(2): 489-496. doi:10.1007/s11033-014-3792-3.