Genetic architecture and genomic predictive ability of apple quantitative traits across environments

Abstract

Implementation of genomic tools is desirable to increase the efficiency of apple breeding. Recently, the multi-environment apple reference population (apple REFPOP) proved useful for rediscovering loci, estimating genomic predictive ability, and studying genotype by environment interactions (G × E). So far, only two phenological traits were investigated using the apple REFPOP, although the population may be valuable when dissecting genetic architecture and reporting predictive abilities for additional key traits in apple breeding. Here we show contrasting genetic architecture and genomic predictive abilities for 30 quantitative traits across up to six European locations using the apple REFPOP. A total of 59 stable and 277 location-specific associations were found using GWAS, 69.2% of which are novel when compared with 41 reviewed publications. Average genomic predictive abilities of 0.18–0.88 were estimated using main-effect univariate, main-effect multivariate, multi-environment univariate, and multi-environment multivariate models. The G × E accounted for up to 24% of the phenotypic variability. This most comprehensive genomic study in apple in terms of trait- environment combinations provided knowledge of trait biology and prediction models that can be readily applied for marker-assisted or genomic selection, thus facilitating increased breeding efficiency.

References

[1] FAOSTAT (Food and Agriculture Organization of the United Nations, 2019).

[2] Cornille A, Giraud T, Smulders MJM et al. The domestication and evolutionary ecology of apples. Trends Genet. 2014; 30: 57-65. https://doi.org/10.1016/j.tig.2013.10.002.

[3] Way RD, Aldwinckle HS, Lamb RC et al. Apples (Malus) . Acta Hortic. 1991; 3-46. https://doi.org/10.17660/ActaHortic.1991.290.1.

[4] Muranty H, Denancé C, Feugey L et al. Using whole-genome SNP data to reconstruct a large multi-generation pedigree in apple germplasm. BMC Plant Biol. 2020; 20: 2. https://doi.org/10.1186/s12870-019-2171-6.

[5] Migicovsky Z, Gardner KM, Richards C et al. Genomic consequences of apple improvement. Hortic Res. 2021; 8: 9. https://doi.org/10.1038/s41438-020-00441-7.

[6] Urrestarazu J, Denancé C, Ravon E et al. Analysis of the genetic diversity and structure across a wide range of germplasm reveals prominent gene flow in apple at the European level. BMC Plant Biol. 2016; 16: 130. https://doi.org/10.1186/s12870-016-0818-0.

[7] Wannemuehler SD, Luby JJ, Yue C et al. A cost-benefit analysis of DNA informed apple breeding. HortScience. 2019; 54: 1998. https://doi.org/10.21273/hortsci14173-19.

[8] Maliepaard C, Alston FH, van Arkel G et al. Aligning male and female linkage maps of apple (Malus pumila Mill.) using multi-allelic markers . Theor Appl Genet. 1998; 97: 60-73. https://doi.org/10.1007/s001220050867.

[9] Kenis K, Keulemans J, Davey MW . Identification and stability of QTLs for fruit quality traits in apple. Tree Genet Genomes. 2008; 4: 647-61. https://doi.org/10.1007/s11295-008-0140-6.

[10] Jänsch M, GAL B, Weger J et al. Identification of SNPs linked to eight apple disease resistance loci. Mol Breed. 2015; 35: 45. https://doi.org/10.1007/s11032-015-0242-4.

[11] Verma S, Evans K, Guan Y et al. Two large-effect QTLs, Ma and Ma3, determine genetic potential for acidity in apple fruit: breeding insights from a multi-family study . Tree Genet Genomes. 2019; 15: 18. https://doi.org/10.1007/s11295-019-1324-y.

[12] Baumgartner IO, Kellerhals M, Costa F et al. Development of SNP-based assays for disease resistance and fruit quality traits in apple (Malus × domestica Borkh.) and validation in breeding pilot studies. Tree Genet Genomes. 2016; 12: 35. https://doi.org/10.1007/s11295-016-0994-y.

[13] Iezzoni AF, McFerson J, Luby J et al. RosBREED: bridging the chasm between discovery and application to enable DNA-informed breeding in rosaceous crops. Hortic Res. 2020; 7: 177. https://doi.org/10.1038/s41438-020-00398-7.

[14] Chagné D, Vanderzande S, Kirk C et al. Validation of SNP markers for fruit quality and disease resistance loci in apple (Malus × domestica Borkh.) using the OpenArray^® platform. Hortic Res. 2019; 6: 30. https://doi.org/10.1038/s41438-018-0114-2.

[15] Velasco R, Zharkikh A, Affourtit J et al. The genome of the domesticated apple (Malus × domestica Borkh.) . Nat Genet. 2010; 42: 833-9. https://doi.org/10.1038/ng.654.

[16] Daccord N, Celton J-M, Linsmith G et al. High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat Genet. 2017; 49: 1099-106. https://doi.org/10.1038/ng.3886.

[17] Zhang L, Hu J, Han X et al. A high-quality apple genome assembly reveals the association of a retrotransposon and red fruit colour. Nat Commun. 2019; 10: 1494. https://doi.org/10.1038/s41467-019-09518-x.

[18] Sun X, Jiao C, Schwaninger H et al. Phased diploid genome assemblies and pan-genomes provide insights into the genetic history of apple domestication. Nat Genet. 2020; 52: 1423-32. https://doi.org/10.1038/s41588-020-00723-9.

[19] Broggini GAL, Schlathölter I, Russo G et al. Chromosome-scale de novo diploid assembly of the apple cultivar ‘Gala Galaxy’. bioRxiv. 2020; 2020.04.25.058891. https://doi.org/10.1101/2020.04.25.058891.

[20] Bianco L, Cestaro A, Sargent DJ et al. Development and validation of a 20K single nucleotide polymorphism (SNP) whole genome genotyping array for apple (Malus × domestica Borkh) . PLoS One. 2014; 9: e110377. https://doi.org/10.1371/journal.pone.0110377.

[21] Bianco L, Cestaro A, Linsmith G et al. Development and validation of the Axiom^® Apple480K SNP genotyping array . Plant J. 2016; 86: 62-74. https://doi.org/10.1111/tpj.13145.

[22] Migicovsky Z, Gardner KM, Money D et al. Genome to phenome mapping in apple using historical data. The Plant Genome. 2016; 9. https://doi.org/10.3835/plantgenome2015.11.0113.

[23] McClure KA, Gardner KM, Douglas GM et al. A genome-wide association study of apple quality and scab resistance. The Plant Genome. 2018; 11: 170075. https://doi.org/10.3835/plantgenome2017.08.0075.

[24] Hirschhorn JN, Daly MJ . Genome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005; 6: 95-108. https://doi.org/10.1038/nrg1521.

[25] Kumar S, Garrick DJ, Bink MCAM et al. Novel genomic approaches unravel genetic architecture of complex traits in apple. BMC Genomics. 2013; 14: 393. https://doi.org/10.1186/1471-2164-14-393.

[26] Urrestarazu J, Muranty H, Denancé C et al. Genome-wide association mapping of flowering and ripening periods in apple. Front Plant Sci. 2017; 8: 1923. https://doi.org/10.3389/fpls.2017.01923.

[27] Larsen B, Migicovsky Z, Jeppesen AA et al. Genome-wide association studies in apple reveal loci for aroma volatiles, sugar composition, and harvest date. The Plant Genome. 2019; 12: 180104. https://doi.org/10.3835/plantgenome2018.12.0104.

[28] Hu Y, Han Z, Sun Y et al. ERF4 affects fruit firmness through TPL4 by reducing ethylene production. Plant J. 2020; 103: 937-50. https://doi.org/10.1111/tpj.14884.

[29] Minamikawa MF, Kunihisa M, Noshita K et al. Tracing founder haplotypes of Japanese apple varieties: application in genomic prediction and genome-wide association study. Hortic Res. 2021; 8: 49. https://doi.org/10.1038/s41438-021-00485-3.

[30] Meuwissen THE, Hayes BJ, Goddard ME . Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001; 157: 1819.

[31] Meuwissen T. Genomic selection: marker assisted selection on a genome wide scale. J Anim Breed Genet. 2007; 124: 321-2. https://doi.org/10.1111/j.1439-0388.2007.00708.x.

[32] Kumar S, Chagné D, Bink MCAM et al. Genomic selection for fruit quality traits in apple (Malus × domestica Borkh.) . PLoS One. 2012; 7: e36674. https://doi.org/10.1371/journal.pone.0036674.

[33] Kumar S, Molloy C, Muñoz P et al. Genome-enabled estimates of additive and nonadditive genetic variances and prediction of apple phenotypes across environments. G3 Genes|Genomes|Genetics. 2015; 5: 2711-8. https://doi.org/10.1534/g3.115.021105.

[34] Muranty H, Troggio M, Sadok IB et al. Accuracy and responses of genomic selection on key traits in apple breeding. Hortic Res. 2015; 2: 15060. https://doi.org/10.1038/hortres.2015.60.

[35] Roth M, Muranty H, Guardo MD et al. Genomic prediction of fruit texture and training population optimization towards the application of genomic selection in apple. Hortic Res. 2020; 7: 148. https://doi.org/10.1038/s41438-020-00370-5.

[36] Jung M, Roth M, Aranzana MJ et al. The apple REFPOP-a reference population for genomics-assisted breeding in apple. Hortic Res. 2020; 7: 189. https://doi.org/10.1038/s41438-020-00408-8.

[37] García-Ruiz A, Cole JB, VanRaden PM et al. Changes in genetic selection differentials and generation intervals in US Holstein dairy cattle as a result of genomic selection. Proc Natl Acad Sci. 2016; 113: E3995-4004. https://doi.org/10.1073/pnas.1519061113.

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

[39] Howard R, Carriquiry AL, Beavis WD . Parametric and nonparametric statistical methods for genomic selection of traits with additive and epistatic genetic architectures. G3: Genes|Genomes|Genetics. 2014; 4: 1027-46. https://doi.org/10.1534/g3.114.010298.

[40] Cooper M, DeLacy IH . Relationships among analytical methods used to study genotypic variation and genotype-by-environment interaction in plant breeding multi-environment experiments. Theor Appl Genet. 1994; 88: 561-72. https://doi.org/10.1007/BF01240919.

[41] Snape JW, Foulkes MJ, Simmonds J et al. Dissecting gene × environmental effects on wheat yields via QTL and physiological analysis. Euphytica. 2007; 154: 401-8. https://doi.org/10.1007/s10681-006-9208-2.

[42] Lopez-Cruz M, Crossa J, Bonnett D et al. Increased prediction accuracy in wheat breeding trials using a marker × environment interaction genomic selection model. G3: Genes|Genomes|Genetics. 2015; 5: 569-82. https://doi.org/10.1534/g3.114.016097.

[43] Jarquín D, Crossa J, Lacaze X et al. A reaction norm model for genomic selection using high-dimensional genomic and environmental data. Theor Appl Genet. 2014; 127: 595-607. https://doi.org/10.1007/s00122-013-2243-1.

[44] Tsai H-Y, Cericola F, Edriss V et al. Use of multiple traits genomic prediction, genotype by environment interactions and spatial effect to improve prediction accuracy in yield data. PLoS One. 2020; 15: e0232665. https://doi.org/10.1371/journal.pone.0232665.

[45] Lado B, Vásquez D, Quincke M et al. Resource allocation optimization with multi-trait genomic prediction for bread wheat (Triticum aestivum L.) baking quality. Theor Appl Genet. 2018; 131: 2719-31. https://doi.org/10.1007/s00122-018-3186-3.

[46] Gianola D, Fernando RL . A multiple-trait Bayesian LASSO for genome-enabled analysis and prediction of complex traits. Genetics. 2020; 214: 305-31. https://doi.org/10.1534/genetics.119.302934.

[47] Watanabe K, Stringer S, Frei O et al. A global overview of pleiotropy and genetic architecture in complex traits. Nat Genet. 2019; 51: 1339-48. https://doi.org/10.1038/s41588-019-0481-0.

[48] Johnston JW, Gunaseelan K, Pidakala P et al. Co-ordination of early and late ripening events in apples is regulated through differential sensitivities to ethylene. J Exp Bot. 2009; 60: 2689-99. https://doi.org/10.1093/jxb/erp122.

[49] Chagné D, Dayatilake D, Diack R et al. Genetic and environmental control of fruit maturation, dry matter and firmness in apple (Malus × domestica Borkh.) . Hortic Res. 2014; 1: 14046. https://doi.org/10.1038/hortres.2014.46.

[50] Costa F, Stella S, Van de Weg WE et al. Role of the genes Md-ACO1 and Md-ACS1 in ethylene production and shelf life of apple (Malus domestica Borkh) . Euphytica. 2005; 141: 181-90. https://doi.org/10.1007/s10681-005-6805-4.

[51] Costa F, Peace CP, Stella S et al. QTL dynamics for fruit firmness and softening around an ethylene-dependent polygalacturonase gene in apple (Malus × domestica Borkh.) . J Exp Bot. 2010; 61: 3029-39. https://doi.org/10.1093/jxb/erq130.

[52] Longhi S, Moretto M, Viola R et al. Comprehensive QTL mapping survey dissects the complex fruit texture physiology in apple (Malus × domestica Borkh.) . J Exp Bot. 2012; 63: 1107-21. https://doi.org/10.1093/jxb/err326.

[53] Longhi S, Hamblin MT, Trainotti L et al. A candidate gene based approach validates Md-PG1 as the main responsible for a QTL impacting fruit texture in apple (Malus × domestica Borkh.) . BMC Plant Biol. 2013; 13: 37. https://doi.org/10.1186/1471-2229-13-37.

[54] Legay S, Guerriero G, Deleruelle A et al. Apple russeting as seen through the RNA-seq lens: strong alterations in the exocarp cell wall. Plant Mol Biol. 2015; 88: 21-40. https://doi.org/10.1007/s11103-015-0303-4.

[55] Johnston JW, Hewett EW, Hertog MLATM . Postharvest softening of apple (Malus domestica) fruit: a review . N Z J Crop Hortic Sci. 2002; 30: 145-60. https://doi.org/10.1080/01140671.2002.9514210.

[56] Rymenants M, van de Weg E, Auwerkerken A et al. Detection of QTL for apple fruit acidity and sweetness using sensorial evaluation in multiple pedigreed full-sib families. Tree Genet Genomes. 2020; 16: 71. https://doi.org/10.1007/s11295-020-01466-8.

[57] van de Weg E, Guardo MD, Jänsch M et al. Epistatic fire blight resistance QTL alleles in the apple cultivar ‘Enterprise’ and selection X-6398 discovered and characterized through pedigree-informed analysis. Mol Breed. 2017; 38: 5. https://doi.org/10.1007/s11032-017-0755-0.

[58] Evans KM, Patocchi A, Rezzonico F et al. Genotyping of pedigreed apple breeding material with a genome-covering set of SSRs: trueness-to-type of cultivars and their parentages. Mol Breed. 2011; 28: 535-47. https://doi.org/10.1007/s11032-010-9502-5.

[59] Heslot N, Yang H-P, Sorrells ME et al. Genomic selection in plant breeding: a comparison of models. Crop Sci. 2012; 52: 146-60. https://doi.org/10.2135/cropsci2011.06.0297.

[60] Strobl C, Boulesteix A-L, Kneib T et al. Conditional variable importance for random forests. BMC Bioinformatics. 2008; 9: 307. https://doi.org/10.1186/1471-2105-9-307.

[61] Hardner CM, Evans K, Brien C et al. Genetic architecture of apple fruit quality traits following storage and implications for genetic improvement. Tree Genet Genomes. 2016; 12: 20. https://doi.org/10.1007/s11295-016-0977-z.

[62] Hardner CM, Hayes BJ, Kumar S et al. Prediction of genetic value for sweet cherry fruit maturity among environments using a 6K SNP array. Hortic Res. 2019; 6: 6. https://doi.org/10.1038/s41438-018-0081-7.

[63] Crossa J, Pérez-Rodríguez P, Cuevas J et al. Genomic selection in plant breeding: methods, models, and perspectives. Trends Plant Sci. 2017; 22: 961-75. https://doi.org/10.1016/j.tplants.2017.08.011.

[64] Cuevas J, Crossa J, Montesinos-López OA et al. Bayesian genomic prediction with genotype × environment interaction kernel models. G3 Genes|Genomes|Genetics. 2017; 7: 41-53. https://doi.org/10.1534/g3.116.035584.

[65] Cirilli M, Micali S, Aranzana MJ et al. The multisite PeachRefPop collection: a true cultural heritage and international scientific tool for fruit trees. Plant Physiol. 2020; 184: 632-46. https://doi.org/10.1104/pp.19.01412.

[66] Laurens F, Aranzana MJ, Arus P et al. An integrated approach for increasing breeding efficiency in apple and peach in Europe. Hortic Res. 2018; 5: 11. https://doi.org/10.1038/s41438-018-0016-3.

[67] Browning SR, Browning BL . Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007; 81: 1084-97. https://doi.org/10.1086/521987.

[68] Rodríguez-Álvarez MX, Boer MP, van Eeuwijk FA et al. Correcting for spatial heterogeneity in plant breeding experiments with P-splines. Spatial Statistics. 2018; 23: 52-71. https://doi.org/10.1016/j.spasta.2017.10.003.

[69] Bates D, Mächler M, Bolker B et al. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015; 67: 1-48.

[70] Gabriel KR . The biplot graphic display of matrices with application to principal component analysis. Biometrika. 1971; 58: 453-67. https://doi.org/10.1093/biomet/58.3.453.

[71] Tibbs Cortes L, Zhang Z, Yu J . Status and prospects of genome-wide association studies in plants. The Plant Genome. 2021; 14: e20077. https://doi.org/10.1002/tpg2.20077.

[72] Huang M, Liu X, Zhou Y et al. BLINK: a package for the next level of genome-wide association studies with both individuals and markers in the millions. GigaScience. 2018; 8. https://doi.org/10.1093/gigascience/giy154.

[73] Tang Y, Liu X, Wang J et al. GAPIT version 2: an enhanced integrated tool for genomic association and prediction. The Plant Genome. 2016; 9. https://doi.org/10.3835/plantgenome2015.11.0120.

[74] Ward JH . Hierarchical grouping to optimize an objective function. J Am Stat Assoc. 1963; 58: 236-44. https://doi.org/10.1080/01621459.1963.10500845.

[75] Breiman L. Random forests. Mach Learn. 2001; 45: 5-32. https://doi.org/10.1023/A:1010933404324.

[76] Habier D, Fernando RL, Kizilkaya K et al. Extension of the bayesian alphabet for genomic selection. BMC Bioinformatics. 2011; 12: 186. https://doi.org/10.1186/1471-2105-12-186.

[77] de los Campos G, Gianola D, Rosa GJM et al. Semi-parametric genomic-enabled prediction of genetic values using reproducing kernel Hilbert spaces methods. Genet Res. 2010; 92: 295-308. https://doi.org/10.1017/S0016672310000285.

[78] VanRaden PM . Efficient methods to compute genomic predictions. J Dairy Sci. 2008; 91: 4414-23. https://doi.org/10.3168/jds.2007-0980.

[79] Wright MN, Ziegler A . Ranger: a fast implementation of random forests for high dimensional data in C++ and R. J Stat Softw. 2017; 77: 1-17. https://doi.org/10.18637/jss.v077.i01.

[80] Pérez P, de los Campos G . Genome-wide regression and prediction with the BGLR statistical package. Genetics. 2014; 198: 483-95. https://doi.org/10.1534/genetics.114.164442.

[81] de los Campos G, Sorensen D, Gianola D . Genomic heritability: what is it. PLoS Genet. 2015; 11: e1005048. https://doi.org/10.1371/journal.pgen.1005048.

[82] Lehermeier C, de Los Campos G, Wimmer V et al. Genomic variance estimates: with or without disequilibrium covariances? J Anim Breed Genet. 2017; 134: 232-41. https://doi.org/10.1111/jbg.12268.