Genome-wide association and RNA-seq analyses reveal a potential gene related to linolenic acid in soybean seeds

Linolenic acid (LA) has poor oxidative stability since it is a polyunsaturated fatty acid. Soybean oil has a high LA content and thus has poor oxidative stability. To identify candidate genes that affect the linolenic acid (LA) content in soybean seeds, a genome-wide association study (GWAS) was performed with 1,060 soybean cultivars collected in China between 2019–2021 and which LA content was measured using matrix-assisted laser desorption/ionization time-of-flight imaging mass spectrometry (MALDI-TOF IMS). A candidate gene, GmWRI14, encoding an APETALA2 (AP2)-type transcription factor, was detected by GWAS in cultivars from all three study years. Multiple sequence alignments showed that GmWRI14 belongs to the plant WRI1 family. The fatty acid contents of different soybean lines were evaluated in transgenic lines with a copy of GmWRI14, control lines without GmWRI14, and the gmwri14 mutant. MALDI-TOF IMS revealed that GmWRI14 transgenic soybeans had a lower LA content with a significant effect on seed size and shape, whereas gmwri14 mutants had a higher LA content. compared to control. The RNA-seq results showed that GmWRI14 suppresses GmFAD3s (GmFAD3B and GmFAD3C) and GmbZIP54 expression in soybean seeds, leading to decreased LA content. Based on the RNA-seq data, yeast one-hybrid (Y1H) and qRT-PCR were performed to confirm the transcriptional regulation of FAD3s by GmWRI14. Our results suggest that FAD3 is indirectly regulated by GmWRI14, representing a new molecular mechanism of fatty acid biosynthesis, in which GmWRI14 regulates LA content in soybean seeds.