Genome-wide association study identifies QTL for eight fruit traits in cultivated tomato (Solanum lycopersicum L.)

Abstract

Genome-wide association study (GWAS) is effective in identifying favorable alleles for traits of interest with high mapping resolution in crop species. In this study, we conducted GWAS to explore quantitative trait loci (QTL) for eight fruit traits using 162 tomato accessions with diverse genetic backgrounds. The eight traits included fruit weight, fruit width, fruit height, fruit shape index, pericarp thickness, locule number, fruit firmness, and brix. Phenotypic variations of these traits in the tomato collection were evaluated with three replicates in field trials over three years. We filtered 34,550 confident SNPs from the 51 K Axiom^® tomato array based on < 10% of missing data and > 5% of minor allele frequency for association analysis. The 162 tomato accessions were divided into seven clusters and their membership coefficients were used to account for population structure along with a kinship matrix. To identify marker-trait associations (MTAs), four phenotypic data sets representing each of three years and combined were independently analyzed in the multilocus mixed model (MLMM). A total of 30 significant MTAs was detected over data sets for eight fruit traits at P < 0.0005. The number of MTA per trait ranged from one (brix) to seven (fruit weight and fruit width). Two SNP markers on chromosomes 1 and 2 were significantly associated with multiple traits, suggesting pleiotropic effects of QTL. Furthermore, 16 of 30 MTAs suggest potential novel QTL for eight fruit traits. These results facilitate genetic dissection of tomato fruit traits and provide a useful resource to develop molecular tools for improving fruit traits via marker-assisted selection and genomic selection in tomato breeding programs.

References

[1] Food & Agriculture Organization of the United Nation. FAOSTAT statistical database. https://search.library.wisc.edu/catalog/999890171702121 (2020).

[2] Celik, I., Gurbuz, N., Uncu, A. T., Frary, A. & Doganlar, S. Genome-wide SNP discovery and QTL mapping for fruit quality traits in inbred backcross lines (IBLs) of solanum pimpinellifolium using genotyping by sequencing. BMC Genomics. 18, 1 (2017).

[3] Lippman, Z. & Tanksley, S. D. Dissecting the genetic pathway to extreme fruit size in tomato using a cross between the small-fruited wild species Lycopersicon pimpinellifolium and L. esculentum var. Giant Heirloom. Genetics. 158, 413-422 (2001).

[4] Munos, S. et al. Increase in tomato locule number is controlled by two single-nucleotide polymorphisms located near WUSCHEL. Plant Physiol. 156, 2244-2254 (2011).

[5] Ranc, N. et al. Genome-wide association mapping in tomato (Solanum lycopersicum) is possible using genome admixture of Solanum lycopersicum var cerasiforme. G3 (Bethesda). 2, 853-864 (2012).

[6] Rodriguez, G. R., Kim, H. J. & van der Knaap, E. Mapping of two suppressors of OVATE (sov) loci in tomato. Heredity. 111, 256-264 (2013).

[7] Xu, J. et al. Phenotypic diversity and association mapping for fruit quality traits in cultivated tomato and related species. Theor. Appl Genet. 126, 567-581 (2013).

[8] Holland, J. B. Genetic architecture of complex traits in plants. Curr. Opin. Plant Biol. 10, 156-161 (2007).

[9] Perez-de-Castro, A. M. et al. Application of genomic tools in plant breeding. Curr. Genomics. 13, 179-195 (2012).

[10] Yu, J. & Buckler, E. S. Genetic association mapping and genome organization of maize. Curr. Opin. Biotechnol. 17, 155-160 (2006).

[11] Gupta, P. K., Rustgi, S. & Kulwal, P. L. Linkage disequilibrium and association studies in higher plants: present status and future prospects. Plant Mol. Biol. 57, 461-485 (2005).

[12] Zhu, C., Gore, M., Buckler, E. S. & Yu, J. Status and prospects of association mapping in plants. Plant Genome-Us. 1, 5-20 (2008).

[13] Hamilton, J. P. et al. Single nucleotide polymorphism discovery in cultivated tomato via sequencing by synthesis. Plant Genome-Us. 5, 17-29 (2012).

[14] Sim, S. C. et al. High-density SNP genotyping of tomato (Solanum lycopersicum L.) reveals patterns of genetic variation due to breeding. PLoS ONE. 7, e45520 (2012).

[15] Consortium, T. G. S. et al. Exploring genetic variation in the tomato (Solanum section Lycopersicon) clade by whole-genome sequencing. Plant J. 80, 136-148 (2014).

[16] Liu, H. et al. An evaluation of genotyping by sequencing (GBS) to map the Breviaristatum-e (ari-e) locus in cultivated barley. BMC Genomics. 15, 104 (2014).

[17] Yamamoto, E. et al. A simulation-based breeding design that uses whole-genome prediction in tomato. Sci. Rep. 6, 19454 (2016).

[18] Segura, V. et al. An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations. Nat. Genet. 44, 825-830 (2012).

[19] Zhang, Z. et al. Mixed linear model approach adapted for genome-wide association studies. Nat. Genet. 42, 355-360 (2010).

[20] Ruggieri, V. et al. An association mapping approach to identify favourable alleles for tomato fruit quality breeding. BMC Plant Biol. 14, 337 (2014).

[21] Sauvage, C. et al. Genome-wide association in tomato reveals 44 candidate loci for fruit metabolic traits. Plant Physiol. 165, 1120-1132 (2014).

[22] Burgos, E. et al. Validated MAGIC and GWAS population mapping reveals the link between vitamin E content and natural variation in chorismate metabolism in tomato. Plant J.: Cell Mol. Biol. 105, 907-923 (2021).

[23] Tieman, D. et al. A chemical genetic roadmap to improved tomato flavor. Science. 355, 391-394 (2017).

[24] Zhang, J. et al. Genome-wide association mapping for tomato volatiles positively contributing to tomato flavor. Front Plant Sci. 6, 1042 (2015).

[25] Zhao, J. et al. Meta-analysis of genome-wide association studies provides insights into genetic control of tomato flavor. Nat. Commun. 10, 1534 (2019).

[26] Lin, T. et al. Genomic analyses provide insights into the history of tomato breeding. Nat. Genet. 46, 1220-1226 (2014).

[27] Sacco, A. et al. Exploring a tomato landraces collection for fruit-related traits by the aid of a high-throughput genomic platform. PLoS ONE. 10, e0137139 (2015).

[28] Phan, N. T. et al. Identification of loci associated with fruit traits using genome-wide single nucleotide polymorphisms in a core collection of tomato (Solanum lycopersicum L.). Sci. Hortic.-Amst. 243, 567-574 (2019).

[29] Mata-Nicolás, E. et al. Exploiting the diversity of tomato: the development of a phenotypically and genetically detailed germplasm collection. Hortic. Res. 7, 1-14 (2020).

[30] Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611-2620 (2005).

[31] Li, M. X., Yeung, J. M., Cherny, S. S. & Sham, P. C. Evaluating the effective numbers of independent tests and significant p-value thresholds in commercial genotyping arrays and public imputation reference datasets. Hum. Genet 131, 747-756 (2012).

[32] Adhikari, P., McNellie, J. & Panthee, D. R. Detection of quantitative trait loci (QTL) associated with the fruit morphology of tomato. Genes (Basel) 11, 1117 (2020).

[33] Brekke, T. D., Stroud, J. A., Shaw, D. S., Crawford, S. & Steele, K. A. QTL mapping in salad tomatoes. Euphytica 215, 115 (2019).

[34] Liu, X., Huang, M., Fan, B., Buckler, E. S. & Zhang, Z. Iterative usage of fixed and random effect models for powerful and efficient genome-wide association studies. PLoS Genet 12, e1005767 (2016).

[35] Safaei, M. et al. Four genetic loci control compact plant size with yellow pear-shaped fruit in ornamental tomato (Solanum lycopersicum L.). Plant Genome-Us 13, e20017 (2020).

[36] Gautier, H. et al. How does tomato quality (sugar, acid, and nutritional quality) vary with ripening stage, temperature, and irradiance? J. Agric. Food Chem. 56, 1241-1250 (2008).

[37] Frary, A. et al. fw2.2: a quantitative trait locus key to the evolution of tomato fruit size. Science 289, 85-88 (2000).

[38] Liu, J., Van Eck, J., Cong, B. & Tanksley, S. D. A new class of regulatory genes underlying the cause of pear-shaped tomato fruit. Proc. Natl Acad. Sci. USA. 99, 13302-13306 (2002).

[39] Liu, X., Geng, X., Zhang, H., Shen, H. & Yang, W. Association and genetic identification of loci for four fruit traits in tomato using InDel markers. Front Plant Sci. 8, 1269 (2017).

[40] Grandillo, S., Ku, H. M. & Tanksley, S. D. Identifying the loci responsible for natural variation in fruit size and shape in tomato. Theor. Appl. Genet. 99, 978-987 (1999).

[41] Hernández-Bautista, A. et al. Fruit size QTLs affect in a major proportion the yield in tomato. Chil. J. Agric. Res. 75, 402-409 (2015).

[42] Bernardo, R. Bandwagons I, too, have known. Theor. Appl. Genet. 129, 2323-2332 (2016).

[43] Heffner, E. L., Sorrells, M. E. & Jannink, J. L. Genomic selection for crop improvement. Crop Sci. 49, 1-12 (2009).

[44] Crossa, J. et al. Genomic selection in plant breeding: methods, models, and perspectives. Trends Plant Sci. 22, 961-975 (2017).

[45] Liabeuf, D., Sim, S. C. & Francis, D. M. Comparison of marker-based genomic estimated breeding values and phenotypic evaluation for selection of bacterial spot resistance in tomato. Phytopathology. 108, 392-401 (2018).

[46] Spindel, J. E. et al. Genome-wide prediction models that incorporate de novo GWAS are a powerful new tool for tropical rice improvement. Heredity. 116, 395-408 (2016).

[47] Zhang, H., Yin, L., Wang, M., Yuan, X. & Liu, X. Factors affecting the accuracy of genomic selection for agricultural economic traits in maize, cattle, and pig populations. Front Genet. 10, 189 (2019).

[48] Kabelka, E., Franchino, B. & Francis, D. M. Two loci from Lycopersicon hirsutum LA407 confer resistance to strains of Clavibacter michiganensis subsp. michiganensis. Phytopathology. 92, 504-510 (2002).

[49] Browning, B. L., Zhou, Y. & Browning, S. R. A One-Penny Imputed Genome from Next-Generation Reference Panels. Am. J. Hum. Genet. 103, 338-348 (2018).

[50] Brewer, M. T. et al. Development of a controlled vocabulary and software application to analyze fruit shape variation in tomato and other plant species. Plant Physiol. 141, 15-25 (2006).

[51] Rohlf, F. J. & Sokal, R. R. Biometry: the principles and practice of statistics in biological research. (Freeman New York, 1981).

[52] Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945-959 (2000).

[53] Nei, M . Estimation of average heterozygosity and genetic distacne from a small number of individuals. Genetics 89, 583-590 (1978).

[54] Kamvar, Z. N., Brooks, J. C. & Grunwald, N. J. Novel R tools for analysis of genome-wide population genetic data with emphasis on clonality. Front Genet 6, 208 (2015).

[55] Lipka, A. E. et al. GAPIT: genome association and prediction integrated tool. Bioinformatics 28, 2397-2399 (2012).

[56] VanRaden, P. M. Efficient methods to compute genomic predictions. J. Dairy Sci. 91, 4414-4423 (2008).

[57] Knoch, D. et al. Strong temporal dynamics of QTL action on plant growth progression revealed through high-throughput phenotyping in canola. Plant Biotechnol. J. 18, 68-82 (2020).