Genomic consequences of apple improvement

Abstract

The apple (Malus domestica) is one of the world’s most commercially important perennial crops and its improvement has been the focus of human effort for thousands of years. Here, we genetically characterise over 1000 apple accessions from the United States Department of Agriculture (USDA) germplasm collection using over 30,000 single-nucleotide polymorphisms (SNPs). We confirm the close genetic relationship between modern apple cultivars and their primary progenitor species, Malus sieversii from Central Asia, and find that cider apples derive more of their ancestry from the European crabapple, Malus sylvestris, than do dessert apples. We determine that most of the USDA collection is a large complex pedigree: over half of the collection is interconnected by a series of first-degree relationships. In addition, 15% of the accessions have a first-degree relationship with one of the top 8 cultivars produced in the USA. With the exception of ‘Honeycrisp’, the top 8 cultivars are interconnected to each other via pedigree relationships. The cultivars ‘Golden Delicious’ and ‘Red Delicious’ were found to have over 60 first-degree relatives, consistent with their repeated use by apple breeders. We detected a signature of intense selection for red skin and provide evidence that breeders also selected for increased firmness. Our results suggest that Americans are eating apples largely from a single family tree and that the apple’s future improvement will benefit from increased exploitation of its tremendous natural genetic diversity.

References

[1] Miller, A. J. & Gross, B. L. From forest to field: perennial fruit crop domestication. Am. J. Bot. 98, 1389-414 (2011).

[2] Food and Agriculture Organization of the United Nations. FAOSTAT. http://faostat3.fao.org/browse/rankings/commodities_by_regions/E (2016).

[3] Vavilov, N. I. Pyat’Kontinentov [Five Continents]. (Leningrad: Nauka, 1987).

[4] Velasco, R. et al. The genome of the domesticated apple (Malus x domestica Borkh.) . Nat. Genet. 42, 833-9 (2010).

[5] Duan, N. et al. Genome re-sequencing reveals the history of apple and supports a two-stage model for fruit enlargement. Nat. Commun. 8, 249, https://doi.org/10.1038/s41467-017-00336-7 (2017).

[6] Cornille, A. et al. New insight into the history of domesticated apple: secondary contribution of the European wild apple to the genome of cultivated varieties. PLoS Genet. 8, e1002703 (2012).

[7] Gross, B. L., Henk, A. D., Richards, C. M., Fazio, G. & Volk, G. M. Genetic diversity in Malus xdomestica (Rosaceae) through time in response to domestication. Am. J. Bot. 101, 1770-9 (2014).

[8] World Apple and Pear Association . European Apple and Pear Crop Forecast. http://www.wapa-association.org/asp/page_1.asp?doc_id=447 (2018).

[9] Proctor, J . The ‘McIntosh’ apple. Fruit. Var. J. 44, 50-53 (1990).

[10] U.S. Apple Association . Apple Industry At-a-Glance. http://usapple.org/the-industry/apple-industry-at-a-glance/ (2018).

[11] Elshire, R. J. et al. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 6, e19379 (2011).

[12] Daccord, N. et al. High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat. Genet. https://doi.org/10.1038/ng.3886 (2017).

[13] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-60 (2009).

[14] Glaubitz, J. C. et al. TASSEL-GBS: a high capacity genotyping by sequencing analysis pipeline. PLoS ONE 9, e90346 (2014).

[15] Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156-8 (2011).

[16] Money, D., Migicovsky, Z., Gardner, K. & Myles, S. LinkImputeR: user-guided genotype calling and imputation for non-model organisms. BMC Genomics 18, 523 (2017).

[17] Purcell, S. PLINK v.1.07. http://pngu.mgh.harvard.edu/purcell/plink/ (2009).

[18] Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559-575 (2007).

[19] Larsen, B. et al. Population structure, relatedness and ploidy levels in an apple gene bank revealed through genotyping-by-sequencing. PLoS ONE 13, e0201889 (2018).

[20] Raj, A., Stephens, M. & Pritchard, J. K. fastSTRUCTURE: variational inference of population structure in large SNP data sets. Genetics 197, 573-89 (2014).

[21] Patterson, N., Price, A. L. & Reich, D. Population structure and eigenanalysis. PLoS Genet. 2, e190 (2006).

[22] McVean, G. A genealogical interpretation of principal components analysis. PLoS Genet. 5, e1000686 (2009).

[23] Migicovsky, Z. et al. Genomic ancestry estimation quantifies use of wild species in grape breeding. BMC Genomics 17, 478 (2016).

[24] Sawler, J. et al. Genomics assisted ancestry deconvolution in grape. PLoS ONE 8, e80791 (2013).

[25] Migicovsky, Z. et al. Genome to phenome mapping in apple using historical data. Plant Genome 9 https://doi.org/10.3835/plantgenome2015.11.0113 (2016).

[26] Patterson, N. et al. Ancient admixture in human history. Genetics 192, 1065-93 (2012).

[27] Weir, B. S. & Cockerham, C. C. Estimating F-statistics for the analysis of population structure. Evolution 38, 1358-1370 (1984).

[28] Scheet, P. & Stephens, M. A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am. J. Hum. Genet. 78, 629-644 (2006).

[29] Szpiech, Z. A. & Hernandez, R. D. selscan: an efficient multithreaded program to perform EHH-based scans for positive selection. Mol. Biol. Evol. 31, 2824-7 (2014).

[30] Jung, S. et al. 15 years of GDR: New data and functionality in the Genome Database for Rosaceae. Nucleic Acids Res. https://doi.org/10.1093/nar/gky1000 (2018).

[31] Alexa, A. & Rahnenfuhrer, J. topGO: Enrichment Analysis for Gene Ontology. (R package version 2.28.0, 2016).

[32] Butts, C. T. network: a package for managing relational data in R. J. Stat. Softw. 24, 1-36 (2008).

[33] Noiton, D. A. & Alspach, P. A. Founding clones, inbreeding, coancestry, and status number of modern apple cultivars. J. Am. Soc. Hortic. Sci. 121, 773-782 (1996).

[34] Hokanson, S., Szewc-McFadden, A., Lamboy, W. & McFerson, J. Microsatellite (SSR) markers reveal genetic identities, genetic diversity and relationships in a Malus× domestica Borkh. core subset collection. Theor. Appl. Genet. 97, 671-683 (1998).

[35] Bradbury, P. J. et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633-2635 (2007).

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

[37] Turner, S. D. qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. J. Open Source Softw. 3, 731, https://doi.org/10.21105/joss.00731 (2018).

[38] Takos, A. M. et al. Light-induced expression of a MYB gene regulates anthocyanin biosynthesis in red apples. Plant Physiol. 142, 1216-32 (2006).

[39] Espley, R. V. et al. Red colouration in apple fruit is due to the activity of the MYB transcription factor, MdMYB10. Plant J. 49, 414-27 (2007).

[40] Ban, Y. et al. Isolation and functional analysis of a MYB transcription factor gene that is a key regulator for the development of red coloration in apple skin. Plant Cell Physiol. 48, 958-70 (2007).

[41] Chagné, D. et al. A functional genetic marker for apple red skin coloration across different environments. Tree Genet. Genomes 12, 67 (2016).

[42] Zhang, L. et al. A high-quality apple genome assembly reveals the association of a retrotransposon and red fruit colour. Nat. Commun. 10, 1494 (2019).

[43] Larsen, B. et al. Genome-wide association studies in apple reveal loci for aroma volatiles, sugar composition, and harvest date. Plant Genome 12 (2019).

[44] Yeats, T. H. et al. Allelic diversity of NAC18.1 is a major determinant of fruit firmness and harvest date in apple (Malus domestica) . https://www.biorxiv.org/content/10.1101/708040v1 (in review).

[45] Kumar, S., Raulier, P., Chagné, D. & Whitworth, C. Molecular-level and trait-level differentiation between the cultivated apple (Malus× domestica Borkh.) and its main progenitor Malus sieversii . Plant Genet. Resour. 12, 330-340 (2014).

[46] Harris, S. A., Robinson, J. P. & Juniper, B. E. Genetic clues to the origin of the apple. Trends Genet. 18, 426-430 (2002).

[47] Cornille, A., Giraud, T., Smulders, M. J. M., Roldán-Ruiz, I. & Gladieux, P. The domestication and evolutionary ecology of apples. Trends Genet. 30, 57-65 (2014).

[48] Volk, G. M. et al. Chloroplast heterogeneity and historical admixture within the genus Malus . Am. J. Bot. 102, 1198-1208 (2015).

[49] Cornille, A. et al. A multifaceted overview of apple tree domestication. Trends Plant Sci. (2019) https://doi.org/10.1016/j.tplants.2019.05.007.

[50] Sun, X. et al. Phased diploid genome assemblies and pan-genomes provide insights into the genetic history of apple domestication. Nat. Genet. 52, 1423-1432, https://doi.org/10.1038/s41588-020-00723-9 (2020).

[51] Myles, S. et al. Genetic structure and domestication history of the grape. Proc. Natl Acad. Sci. 108, 3530-3535 (2011).

[52] Dobrenel, T. et al. Regulation of plant growth and metabolism by the TOR kinase. Biochem. Soc. Trans. 39, 477-481 (2011).

[53] Kelsey, D. & Brown, S. `McIntosh Wijcik': a columnar mutation of `McIntosh' apple proving useful in physiology and breeding research. Fruit Var. J. 46, 83-87 (1992).

[54] Luby, J. J. & Bedford, D. S. Cultivars as consumer brands: trends in protecting and commercializing apple cultivars via intellectual property rights. Crop Sci. 55, 2504-2510 (2015).

[55] Komatsu, H. Red Fuji in Japan: choosing the best strain for your business strategy. Compact Fruit Tree 31, 44-45 (1998).

[56] Roach, M. J. et al. Population sequencing reveals clonal diversity and ancestral inbreeding in the grapevine cultivar Chardonnay. PLoS Genet. 14, e1007807 (2018).

[57] Ballard, J. K. Granny Smith: an Important Apple for the Pacific Northwest (Pullman, 1981).

[58] Muranty, H. et al. Using whole-genome SNP data to reconstruct a large multi-generation pedigree in apple germplasm. BMC Plant Biol. 20, 2 (2020).

[59] Warschefsky, E. J. & von Wettberg, E. J. B. Population genomic analysis of mango (Mangifera indica) suggests a complex history of domestication . N. Phytol. 222, 2023-2037 (2019).

[60] Cabe, P. R., Baumgarten, A., Onan, K., Luby, J. J. & Bedford, D. S. Using microsatellite analysis to verify breeding records: a study of ‘honeycrisp’ and other cold-hardy apple cultivars. HortScience 40, 15-17 (2005).

[61] Howard, N. P. et al. Elucidation of the ‘Honeycrisp’ pedigree through haplotype analysis with a multi-family integrated SNP linkage map and a large apple (Malus x domestica) pedigree-connected SNP data set . Hortic. Res. 4, 17003 (2017).

[62] Migicovsky, Z. & Myles, S. Exploiting wild relatives for genomics-assisted breeding of perennial crops. Front. Plant Sci. 460 (2017).

[63] Valois, S., Merwin, I. A. & Padilla-Zakour, O. I. Characterization of fermented cider apple varieties grown in upstate New York. J. Am. Pomol. Soc. 60, 113-128 (2006).

[64] Sanoner, P., Guyot, S., Marnet, N., Molle, D. & Drilleau, J.-F. Polyphenol profiles of French cider apple varieties (Malus domestica sp.) . J. Agric. Food Chem. 47, 4847-4853 (1999).

[65] Hamadziripi, E. T., Theron, K. I., Muller, M. & Steyn, W. J. Apple compositional and peel color differences resulting from canopy microclimate affect consumer preference for eating quality and appearance. HortScience 49, 384-392 (2014).

[66] Hazzouri, K. M. et al. Genome-wide association mapping of date palm fruit traits. Nat. Commun. 10, 4680 (2019).

[67] Huang, D. et al. Subfunctionalization of the Ruby2-Ruby1 gene cluster during the domestication of citrus. Nat. Plants 4, 930-941 (2018).

[68] Legun, K. Ever-redder apples: How aesthetics shape the biology of markets. In Biological Economies: Experimentation and the politics of agri-food frontiers (eds Heron, R. L., Campbell, H., Lewis, N. & Carolan, M.) 127-140 (Routledge, 2016).

[69] Hampson, C. R. & Kemp, H. in Apples: Botany, Production and Uses (eds. Ferree, D. C. & Warrington, I. J.) 61-89 (CABI Publishing, 2003).

[70] Laurens, F. et al. An integrated approach for increasing breeding efficiency in apple and peach in Europe. Hortic. Res. 5, 11 (2018).

[71] Wu, J. et al. Diversification and independent domestication of Asian and European pears. Genome Biol. 19, 77 (2018).

[72] Urrestarazu, J. et al. Genome-wide association mapping of flowering and ripening periods in apple. Front. Plant Sci. 8, 1923 (2017).

[73] McClure, K. A. et al. A genome-wide association study of apple quality and scab resistance. Plant Genome 11, 170075 (2018).