Atlas of phenotypic, genotypic and geographical diversity present in the European traditional tomato

Abstract

The Mediterranean basin countries are considered secondary centres of tomato diversification. However, information on phenotypic and allelic variation of local tomato materials is still limited. Here we report on the evaluation of the largest traditional tomato collection, which includes 1499 accessions from Southern Europe. Analyses of 70 traits revealed a broad range of phenotypic variability with different distributions among countries, with the culinary end use within each country being the main driver of tomato diversification. Furthermore, eight main tomato types (phenoclusters) were defined by integrating phenotypic data, country of origin, and end use. Genome-wide association study (GWAS) meta-analyses identified associations in 211 loci, 159 of which were novel. The multidimensional integration of phenoclusters and the GWAS meta-analysis identified the molecular signatures for each traditional tomato type and indicated that signatures originated from differential combinations of loci, which in some cases converged in the same tomato phenotype. Our results provide a roadmap for studying and exploiting this untapped tomato diversity.

References

[1] Gentilcore D . Pomodoro! Columbia University Press; 2010.

[2] López-Terrada M. The History of the Arrival of the Tomato in Europe: An Initial Overview . 2014: 1-17.

[3] Smith AF . The Tomato in America: Early History, Culture, and Cookery . Columbia, S.C: University of South Carolina Press; 1994.

[4] Robbins MD, Sim SC, Yang W et al. Mapping and linkage disequilibrium analysis with a genome-wide collection of SNPs that detect polymorphism in cultivated tomato. J Exp Bot. 2011; 62: 1831-45.

[5] Sim SC, Robbins MD, Deynze A et al. Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.) . 2011; 106: 927-35.

[6] Conesa MÀ, Fullana-Pericàs M, Granell A et al. Mediterranean long shelf-life landraces: an untapped genetic resource for tomato improvement. Front Plant Sci. 2020; 10: 1651.

[7] Esposito S, Cardi T, Campanelli G et al. ddRAD sequencing-based genotyping for population structure analysis in cultivated tomato provides new insights into the genomic diversity of Mediterranean ‘da serbo’ type long shelf-life germplasm. Hortic Res. 2020; 7: 134.

[8] Villand J, Skroch PW, Lai T et al. Genetic variation among tomato accessions from primary and secondary centers of diversity. Crop Sci. 1998; 38: 1339-47.

[9] Rodríguez GR, Muños S, Anderson C et al. Distribution of SUN, OVATE, LC, and FAS in the tomato germplasm and the relationship to fruit shape diversity. Plant Physiol. 2011; 156: 275-85.

[10] Casals J, Pascual L, Cañizares J et al. Genetic basis of long shelf life and variability into Penjar tomato. Genet Resour Crop Evol. 2012; 59: 219-29.

[11] Mazzucato A, Papa R, Bitocchi E et al. Genetic diversity, structure and marker-trait associations in a collection of Italian tomato (Solanum lycopersicum L.) landraces . Theor Appl Genet. 2008; 116: 657-69.

[12] Blanca J, Pons C, Montero-Pau J et al. European vintage tomatoes galore: a result of farmers combinatorial assorting/swapping of a few diversity rich loci. J Exp Botany. 2022: erac072. https://doi.org/10.1093/jxb/erac072.

[13] Ruiz JJ, Garcia-Martinez S, Pico B et al. Genetic variability and relationship of closely related Spanish traditional cultivars of tomato as detected by SRAP and SSR markers. J Am Soc Hortic Sci. 2005; 130: 88-94.

[14] Casals J, Pascual L, Canizares J et al. The risks of success in quality vegetable markets: possible genetic erosion in Marmande tomatoes (Solanum lycopersicum L.) and consumer dissatisfaction . Sci Hortic (Amsterdam). 2011; 130: 78-84.

[15] Parisi M, Aversano R, Graziani G et al. Phenotypic and molecular diversity in a collection of ‘Pomodoro di Sorrento’ Italian tomato landrace. Sci Hortic (Amsterdam). 2016; 203: 143-51.

[16] Koutsika-Sotiriou M, Mylonas I, Tsivelikas A et al. Sci Hortic (Amsterdam) . 2016; 198: 78-85.

[17] Baldina S, Picarella ME, Troise AD et al. Metabolite profiling of Italian tomato landraces with different fruit types. Front Plant Sci. 2016; 7: 1-13.

[18] Andreakis N, Giordano I, Pentangelo A et al. DNA fingerprinting and quality traits of corbarino cherry-like tomato landraces. J Agric Food Chem. 2004; 52: 3366-71.

[19] Mercati F, Longo C, Poma D et al. Genetic variation of an Italian long shelf-life tomato (Solanum lycopersicon L.) collection by using SSR and morphological fruit traits . Genet Resour Crop Evol. 2015; 62: 721-32.

[20] Aflitos S, Schijlen E, De Jong H et al. Exploring genetic variation in the tomato (Solanum section Lycopersicon) clade by whole-genome sequencing. Plant J. 2014; 80: 136-48.

[21] Blanca J, Montero-Pau J, Sauvage C et al. Genomic variation in tomato, from wild ancestors to contemporary breeding accessions. BMC Genomics. 2015; 16: 1-19.

[22] Razifard H, Ramos A, Della Valle AL et al. Genomic evidence for complex domestication history of the cultivated tomato in Latin America. Mol Biol Evol. 2020; 37: 1118-32.

[23] Causse M, Desplat N, Pascual L et al. Whole genome resequencing in tomato reveals variation associated with introgression and breeding events. BMC Genomics. 2013; 14: 791.

[24] Mata-Nicolás E, Montero-Pau J, Gimeno-Paez E et al. Exploiting the diversity of tomato: the development of a phenotypically and genetically detailed germplasm collection. Hortic Res. 2020; 7: 66.

[25] Tieman D, Zhu G, Resende MFR et al. A chemical genetic roadmap to improved tomato flavor. Science. 2017; 355: 391-4.

[26] Lin T, Zhu G, Zhang J et al. Genomic analyses provide insights into the history of tomato breeding. Nat Genet. 2014; 46: 1220-6.

[27] Cortés-Olmos C, Vilanova S, Pascual L et al. SNP markers applied to the characterization of Spanish tomato (Solanum lycopersicum L.) landraces . Sci Hortic (Amsterdam). 2015; 194: 100-10.

[28] Sacco A, Ruggieri V, Parisi M et al. Exploring a tomato landraces collection for fruit-related traits by the aid of a high-throughput genomic platform. PLoS One. 2015; 10: 1-20.

[29] Mazzucato A, Olimpieri I, Siligato F et al. Characterization of genes controlling stamen identity and development in a parthenocarpic tomato mutant indicates a role for the DEFICIENS ortholog in the control of fruit set. Physiol Plant. 2008; 132: 526-37.

[30] Terzopoulos PJ, Bebeli PJ . Phenotypic diversity in Greek tomato (Solanum lycopersicum L.) landraces . Sci Hortic (Amsterdam). 2010; 126: 138-44.

[31] Nankar AN, Tringovska I, Grozeva S et al. Tomato phenotypic diversity determined by combined approaches of conventional and high-throughput tomato analyzer phenotyping. Plan Theory. 2020; 9.

[32] Brigitte E. Analyses Factorielles Simples et Multiples: Objectifs, Méthodes et Interprétation / Brigitte Escofier,... Jérôme Pagès . 3e édition ed. Paris: Dunod; 1998.

[33] Pagès J. Analyse factorielle de données mixtes. Rev Stat Appliquée. 52: 93-111.

[34] Husson F, Josse J, Pages J . Principal component methods-hierarchical clustering-partitional clustering: why would we need to choose for visualizing data. Appl Math Dep. 2010; 1-17.

[35] Devlin B, Roeder K . Genomic control for association studies. Biometrics. 1999; 55: 997-1004.

[36] Yang J, Weedon MN, Purcell S et al. Genomic inflation factors under polygenic inheritance. Eur J Hum Genet. 2011; 19: 807-12.

[37] Evangelou E, Ioannidis JPA . Meta-analysis methods for genome-wide association studies and beyond. Nat Rev Genet. 2013; 14: 379-89.

[38] Cong B, Barrero LS, Tanksley SD . Regulatory change in YABBY-like transcription factor led to evolution of extreme fruit size during tomato domestication. Nat Genet. 2008; 40: 800-4.

[39] Xu C, Liberatore KL, Macalister CA et al. A cascade of arabinosyltransferases controls shoot meristem size in tomato. Nat Genet. 2015; 47: 784-92.

[40] Chu YH, Jang JC, Huang Z et al. Tomato locule number and fruit size controlled by natural alleles of lc and fas. Plant Direct. 2019; 3: e00142.

[41] Ballester AR, Molthoff J, de Vos R et al. Biochemical and molecular analysis of pink tomatoes: deregulated expression of the gene encoding transcription factor SlMYB12 leads to pink tomato fruit color. Plant Physiol. 2010; 152: 71 LP- 84.

[42] Barry CS, Giovannoni JJ . Ripening in the tomato Green-ripe mutant is inhibited by ectopic expression of a protein that disrupts ethylene signaling. Proc Natl Acad Sci. 2006; 103: 7923 LP- 7928.

[43] Liu J, Van Eck J, Cong B et al. A new class of regulatory genes underlying the cause of pear-shaped tomato fruit. Proc Natl Acad Sci. 2002; 99: 13302 LP- 13306.

[44] Chakrabarti M, Zhang N, Sauvage C et al. A cytochrome P450 regulates a domestication trait in cultivated tomato. Proc Natl Acad Sci. 2013; 110: 17125 LP- 17130.

[45] Yuste-Lisbona FJ, Fernández-Lozano A, Pineda B et al. ENO regulates tomato fruit size through the floral meristem development network . Proc Natl Acad Sci. 2020; 117: 8187 LP- 8195.

[46] Busch BL, Schmitz G, Rossmann S et al. Shoot branching and leaf dissection in tomato are regulated by homologous gene modules. Plant Cell. 2011; 23: 3595-609.

[47] Vrebalov J, Pan IL, Arroyo AJM et al. Fleshy fruit expansion and ripening are regulated by the tomato SHATTERPROOF gene TAGL1 . Plant Cell. 2009; 21: 3041-62.

[48] van der Knaap E, Sanyal A, Jackson SA et al. High-resolution fine mapping and fluorescence in situ hybridization analysis of sun, a locus controlling tomato fruit shape, reveals a region of the tomato genome prone to DNA rearrangements. Genetics. 2004; 168: 2127 LP- 2140.

[49] Xiao H, Jiang N, Schaffner E et al. A retrotransposon-mediated gene duplication underlies morphological variation of tomato fruit. Science (80-). 2008; 319: 1527 LP- 1530.

[50] Sagi M, Scazzocchio C, Fluhr R . The absence of molybdenum cofactor sulfuration is the primary cause of the flacca phenotype in tomato plants. Plant J. 2002; 31: 305-17.

[51] Nitsch L, Kohlen W, Oplaat C et al. ABA-deficiency results in reduced plant and fruit size in tomato. J Plant Physiol. 2012; 169: 878-83.

[52] Powell ALT, Nguyen CV, Hill T et al. Uniform ripening encodes a Golden 2-like transcription factor regulating tomato fruit chloroplast development. Science (80-). 2012; 336: 1711-5.

[53] Wu S, Zhang B, Keyhaninejad N et al. A common genetic mechanism underlies morphological diversity in fruits and other plant organs. Nat Commun. 2018; 9: 4734.

[54] Schmitz G, Tillmann E, Carriero F et al. The tomato Blind gene encodes a MYB transcription factor that controls the formation of lateral meristems. Proc Natl Acad Sci. 2002; 99: 1064 LP- 1069.

[55] Mu Q, Huang Z, Chakrabarti M et al. Fruit weight is controlled by cell size regulator encoding a novel protein that is expressed in maturing tomato fruits. PLoS Genet. 2017; 13: e1006930.

[56] Su L, Bassa C, Audran C et al. The auxin Sl-IAA17 transcriptional repressor controls fruit size via the regulation of Endoreduplication-related cell expansion. Plant Cell Physiol. 2014; 55: 1969-76.

[57] Zhang Y, Zhang B, Yang T et al. The GAMYB-like gene SlMYB33 mediates flowering and pollen development in tomato. Hortic Res. 2020; 7: 133.

[58] How Kit A, Boureau L, Stammitti-Bert L et al. Functional analysis of SlEZ1 a tomato enhancer of zeste (E(z)) gene demonstrates a role in flower development. Plant Mol Biol. 2010; 74: 201-13.

[59] Yang L, Qi S, Touqeer A et al. SlGT11 controls floral organ patterning and floral determinacy in tomato. BMC Plant Biol. 2020; 20: 562.

[60] Czerednik A, Busscher M, Angenent GC et al. The cell size distribution of tomato fruit can be changed by overexpression of CDKA1. Plant Biotechnol J. 2015; 13: 259-68.

[61] van der Biezen EA, Brandwagt BF, van Leeuwen W et al. Identification and isolation of the FEEBLY gene from tomato by transposon tagging. Mol Gen Genet. 1996; 251: 267-80.

[62] Bannenberg G, Martínez M, Rodríguez MJ et al. Functional analysis of α-DOX2, an active α-dioxygenase critical for Normal development in tomato plants . Plant Physiol. 2009; 151: 1421-32.

[63] Giménez E, Pineda B, Capel J et al. Functional analysis of the Arlequin mutant corroborates the essential role of the ARLEQUIN/TAGL1 gene during reproductive development of tomato. PLoS One. 2010; 5: e14427.

[64] Daminato M, Masiero S, Resentini F et al. Characterization of TM8, a MADS-box gene expressed in tomato flowers. BMC Plant Biol. 2014; 14: 319.

[65] Rodriguez-Leal D, Xu C, Kwon CT et al. Evolution of buffering in a genetic circuit controlling plant stem cell proliferation. Nat Genet. 2019; 51: 786-92.

[66] Molinero-Rosales N, Jamilena M, Zurita S et al. FALSIFLORA, the tomato orthologue of FLORICAULA and LEAFY, controls flowering time and floral meristem identity. Plant J. 1999; 20: 685-93.

[67] Rodrigo MJ, Lado J, Alós E et al. A mutant allele of ζ-carotene isomerase (Z-ISO) is associated with the yellow pigmentation of the “Pinalate” sweet orange mutant and reveals new insights into its role in fruit carotenogenesis. BMC Plant Biol. 2019; 19: 465.

[68] Waseem M, Li N, Su D et al. Overexpression of a basic helix-loop-helix transcription factor gene, SlbHLH22, promotes early flowering and accelerates fruit ripening in tomato (Solanum lycopersicum L.) . Planta. 2019; 250: 173-85.

[69] Kang JH, McRoberts J, Shi F et al. The flavonoid biosynthetic enzyme chalcone isomerase modulates terpenoid production in glandular trichomes of tomato. Plant Physiol. 2014; 164: 1161-74.

[70] Ma N, Feng H, Meng X et al. Overexpression of tomato SlNAC1 transcription factor alters fruit pigmentation and softening. BMC Plant Biol. 2014; 14: 351.

[71] Fantini E, Sulli M, Zhang L et al. Pivotal roles of Cryptochromes 1a and 2 in tomato development and physiology. Plant Physiol. 2019; 179: 732 LP- 748.

[72] Centeno DC, Osorio S, Nunes-Nesi A et al. Malate plays a crucial role in starch metabolism, ripening, and soluble solid content of tomato fruit and affects postharvest softening. Plant Cell. 2011; 23: 162-84.

[73] Sim SC, Van Deynze A, Stoffel K et al. High-density SNP genotyping of tomato (Solanum lycopersicum L.) reveals patterns of genetic variation due to breeding . PLoS One. 2012; 7: e45520.

[74] Anne F, Clint NT, Amy F et al. fw2.2: a quantitative trait locus key to the evolution of tomato fruit size . Science (80-). 2000; 289: 85-8.

[75] Sierra-Orozco E, Shekasteband R, Illa-Berenguer E et al. Identification and characterization of GLOBE, a major gene controlling fruit shape and impacting fruit size and marketability in tomato. Hortic Res. 2021; 8: 138.

[76] Sacco A, Di Matteo A, Lombardi N et al. Quantitative trait loci pyramiding for fruit quality traits in tomato. Mol Breed. 2013; 31: 217-22.

[77] Alonge M, Wang X, Benoit M et al. Major impacts of widespread structural variation on gene expression and crop improvement in tomato. Cell. 2020; 182: 145-161.e23.

[78] Domínguez M, Dugas E, Benchouaia M et al. The impact of transposable elements on tomato diversity. Nat Commun. 2020; 11: 4058.

[79] Soyk S, Lemmon ZH, Sedlazeck FJ et al. Duplication of a domestication locus neutralized a cryptic variant that caused a breeding barrier in tomato. Nat Plants. 2019; 5: 471-9.

[80] Racine JS . RStudio: a platform-independent IDE for R and Sweave. J Appl Econ. 2012; 27: 167-72.

[81] Josse J, Husson F . missMDA: a package for handling missing values in multivariate data analysis. J Stat Software. 2016; 1.

[82] Comtois D. Summary Tools: Tools to Quickly and Neatly Summarize Data R Package Version 0.9.4. 2019.

[83] Frankfort-Nachmias C. Social Statistics for a Diverse Society / Chava Frankfort-Nachmias, Anna Leon-Guerrero . Second ed. Thousand Oaks, Calif.: Pine Forge Press; 2000.

[84] Metsalu T, Vilo J . ClustVis: a web tool for visualizing clustering of multivariate data using principal component analysis and heatmap. Nucleic Acids Res. 2015; 43: W566-70.

[85] Wickham H, Averick M, Bryan J et al. Welcome to the Tidyverse. J Open Source Softw. 2019; 4: 1686.

[86] Villanueva RAM, Chen ZJ . ggplot2: Elegant graphics for data analysis (2nd ed.). Meas Interdiscip Res Perspect 2019; 17: 160-7.

[87] Friendly M . Mosaic displays for multi-way contingency tables. J Am Stat Assoc. 1994; 89: 190-200.

[88] Friendly M, Meyer D, Zeileis A. Discrete Data Analysis with R: Visualization and Modeling Techniques for Categorical and Count Data. 2015.

[89] Simko TW and V. R package “corrplot”: Visualization of a Correlation Matrix. (Version 084). 2017.

[90] Lê S, Josse J, Husson F . FactoMineR: an R package for multivariate analysis. J Stat Software. 2008;Vol 1.

[91] Kassambara A. Practical Guide To Principal Component Methods in R: PCA, M(CA), FAMD, MFA, HCPC, Factoextra . Sthda.com, 2017.

[92] Kassambara A. Multivariate analysis II: practical guide to principal component methods in R. Sthda. 2017; 170.

[93] Kelso NV, Patterson T . Natural earth vector. Cartogr Perspect. 2009; 0: 45-50.

[94] Money D, Gardner K, Migicovsky Z et al. LinkImpute: fast and accurate genotype imputation for nonmodel organisms. G3 (Bethesda). 2015; 5: 2383-90.

[95] Bradbury PJ, Zhang Z, Kroon DE et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007; 23: 2633-5.

[96] Breseghello F, Sorrells ME . Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) . Cultivars Genetics. 2006; 172: 1165-77.

[97] Tang Y, Liu X, Wang J et al. GAPIT version 2: an enhanced integrated tool for genomic association and prediction. Plant Genome. 2016; 9.

[98] Willer CJ, Li Y, Abecasis GR . METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010; 26: 2190-1.

[99] Szklarczyk D, Gable AL, Lyon D et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019; 47: D607-13.

[100] Quinlan AR, Hall IM . BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26: 841-2.

[101] Kanai M, Akiyama M, Takahashi A et al. Genetic analysis of quantitative traits in the Japanese population links cell types to complex human diseases. Nat Genet. 2018; 50: 390-400.