Genome-wide association study identifies candidate genes related to the linoleic acid content in soybean seeds

Soybean (Glycine max (L.) Merrill) oil is a complex mixture of five fatty acids (palmitic, stearic, oleic, linoleic, and linolenic). The high content of linoleic acid (LA) contributes to the oil having poor oxidative stability. Therefore, soybean seed with a lower LA content is desirable. To investigate the genetic architecture of LA, we performed a genome-wide association study (GWAS) using 510 soybean cultivars collected from China. The phenotypic identification results showed that the content of LA varied from 36.22% to 72.18%. The GWAS analysis showed that there were 37 genes related to oleic acid content, with a contribution rate of 7%. The candidate gene Glyma.04G116500.1 (GmWRI14) on chromosome 4 was detected in three consecutive years. The GmWRI14 showed a negative correlation with the LA content and the correlation coefficient was −0.912. To test whether GmWRI14 can lead to a lower LA content in soybean, we introduced GmWRI14 into the soybean genome. Matrix-assisted laser desorption/ionization time-of-flight imaging mass spectrometry (MALDI-TOF IMS) showed that the overexpression of GmWRI14 leads to a lower LA content in soybean seeds. Meanwhile, RNA-seq verified that GmWRI14-overexpressed soybean lines showed a lower accumulation of GmFAD2-1A and GmFAD2-1B than control lines. Our results indicate that the down-regulation of the FAD2 gene triggered by the transcription factor GmWRI14 is the underlying mechanism reducing the LA level of seed. Our results provide novel insights into the genetic architecture of LA and pinpoint potential candidate genes for further in-depth studies.