【摘要】【目的】挖掘适应坝上生态条件的高效结瘤大豆种质，鉴定调控大豆-根瘤菌共生结瘤的遗传位点和候选基因，提高大豆共生固氮效率。【方法】以包含260份大豆种质的自然群体为研究对象，在河北坝上室外盆栽条件下接种根瘤菌菌株USDA110。以单株根瘤数量、单株根瘤干重数值作为表型数据，结合大豆260份种质基因型数据，对其进行全基因组关联分析，挖掘大豆共生结瘤相关基因。【结果】共检测到18个与大豆根瘤数量显著关联的SNP，分别位于第2、7、8、13、18和19染色体上。其中，位于第2染色体上的显著关联位点BARC＿2.01＿Chr02＿43161654＿A＿G是控制大豆根瘤数量的主效位点（LOD=3.89），对该位点上下游共200 kb区间进行连锁不平衡分析，在包含BARC＿2.01＿Chr02＿43161654＿A＿G的连锁区间内，共获得10个调控大豆根瘤数量候选基因，通过单倍型分析发现，Glyma.02G243200不同单倍型所对应的材料之间的根瘤数量具有显著差异（P<0.05），在SoyBase数据库中查询该基因的表达模式发现，Glyma.02G243200在根毛中表达。该基因可能是影响大豆根瘤数量的关键基因。共检测到6个与大豆根瘤干重显著关联的SNP，分别位于第6、18和20染色体上。其中，位于第6染色体上的显著关联位点BARC＿2.01＿Chr06＿6069381＿G＿A和BARC＿2.01＿Chr06＿6192925＿T＿C是控制大豆根瘤干重的主效位点（LOD=3.49和LOD=3.35），对BARC＿2.01＿Chr06＿6069381＿G＿A上游100 kb和BARC＿2.01＿Chr06＿6192925＿T＿C下游100 kb区间进行连锁不平衡分析，获得14个调控大豆根瘤干重候选基因，并进行单倍型分析。其中，Glyma.06G079600和Glyma.06G079900不同单倍型所对应的材料之间的根瘤干重具有显著差异（P<0.01,P<0.001），在SoyBase数据库中查询这两个基因的表达模式发现，Glyma.06G079600和Glyma.06G079900在根中表达。这两个基因可能是影响大豆根瘤干重的关键基因。【结论】在第2染色体上发现一个与根瘤数量显著相关的候选基因，在第6染色体上发现2个与根瘤干重显著相关的候选基因。