Genome-wide association study on resistance of cultivated soybean to Fusarium oxysporum root rot in Northeast China

Abstract

Abstract Background Fusarium oxysporum is a prevalent fungal pathogen that diminishes soybean yield through seedling disease and root rot. Preventing Fusarium oxysporum root rot (FORR) damage entails on the identification of resistance genes and developing resistant cultivars. Therefore, conducting fine mapping and marker development for FORR resistance genes is of great significance for fostering the cultivation of resistant varieties. In this study, 350 soybean germplasm accessions, mainly from Northeast China, underwent genotyping using the SoySNP50K Illumina BeadChip, which includes 52,041 single nucleotide polymorphisms (SNPs). Their resistance to FORR was assessed in a greenhouse. Genome-wide association studies utilizing the general linear model, mixed linear model, compressed mixed linear model, and settlement of MLM under progressively exclusive relationship models were conducted to identify marker-trait associations while effectively controlling for population structure. Results The results demonstrated that these models effectively managed population structure. Eight SNP loci significantly associated with FORR resistance in soybean were detected, primarily located on Chromosome 6. Notably, there was a strong linkage disequilibrium between the large-effect SNPs ss715595462 and ss715595463, contributing substantially to phenotypic variation. Within the genetic interval encompassing these loci, 28 genes were present, with one gene Glyma.06G088400 encoding a protein kinase family protein containing a leucine-rich repeat domain identified as a potential candidate gene in the reference genome of Williams82. Additionally, quantitative real-time reverse transcription polymerase chain reaction analysis evaluated the gene expression levels between highly resistant and susceptible accessions, focusing on primary root tissues collected at different time points after F. oxysporum inoculation. Among the examined genes, only this gene emerged as the strongest candidate associated with FORR resistance. Conclusions The identification of this candidate gene Glyma.06G088400 improves our understanding of soybean resistance to FORR and the markers strongly linked to resistance can be beneficial for molecular marker-assisted selection in breeding resistant soybean accessions against F. oxysporum.