Multiple-model GWAS identifies optimal allelic combinations of quantitative trait loci for malic acid in tomato

Malic acid (MA) is an important flavor acid in fruits and acts as a mediator in a series of metabolic pathways. It is important to understand the factors affecting MA metabolism for fruit flavor improvement and to understand MA-mediated biological processes. However, the metabolic accumulation of MA is controlled by complex heredity and environmental factors, making it difficult to predict and regulate the metabolism of MA. In this study, we carried out a genome-wide association study (GWAS) on MA using eight milestone models with two-environment repeats. A series of associated SNP variations were identified from the GWAS, and 15 high-confidence annotated genes were further predicted based on linkage disequilibrium and lead SNPs. The transcriptome data of candidate genes were explored within different tomato organs as well as various fruit tissues, and suggested specific expression patterns in fruit pericarp. Based on the genetic parameters of population differentiation and SNP distribution, tomato MA content has been more influenced by domestication sweeps and less affected by improvement sweeps in the long-term history of tomato breeding. In addition, genotype × environment interaction might contribute to the difference in domestication phenotypic data under different environments. This study provides new genetic insights into how tomato changed its MA content during breeding and makes available function-based markers for breeding by marker-assisted selection.