Almond diversity and homozygosity define structure, kinship, inbreeding, and linkage disequilibrium in cultivated germplasm, and reveal genomic associations with nut and seed weight

Abstract

Almond [ Prunus dulcis Miller (D.A. Webb)] is the main tree nut species worldwide. Here, genotyping-by-sequencing (GBS) was applied to 149 almond cultivars from the ex situ collections of the Italian Council for Agricultural Research (CREA) and the Spanish National Research Council (CSIC), leading to the detection of 93,119 single-nucleotide polymorphisms (SNPs). The study of population structure outlined four distinct genetic groups and highlighted diversification between the Mediterranean and Californian gene pools. Data on SNP diversity and runs of homozygosity (ROHs) allowed the definition of kinship, inbreeding, and linkage disequilibrium (LD) decay in almond cultivated germplasm. Four-year phenotypic observations, gathered on 98 cultivars of the CREA collection, were used to perform a genome-wide association study (GWAS) and, for the first time in a crop species, homozygosity mapping (HM), resulting in the identification of genomic associations with nut, shell, and seed weight. Both GWAS and HM suggested that loci controlling nut and seed weight are mostly independent. Overall, this study provides insights on the almond cultivation history and delivers information of major interest for almond genetics and breeding. In a broader perspective, our results encourage the use of ROHs in crop science to estimate inbreeding, choose parental combinations minimizing the risk of inbreeding depression, and identify genomic footprints of selection for specific traits.

References

[1] Delplancke, M. et al. Evolutionary history of almond tree domestication in the Mediterranean basin. Mol. Ecol. 22, 1092-1104 (2013).

[2] Acquaah, G. (ed.) Principles of Plant Genetics and Breeding (Wiley-Blackwell, 2012).

[3] De Lorenzis, G. et al. SNP genotyping elucidates the genetic diversity of Magna Graecia grapevine germplasm and its historical origin and dissemination. BMC Plant Biol. 19, 7 (2019).

[4] Pavan, S. et al. Genotyping by sequencing of cultivated lentil (Lens culinaris Medik.) highlights population structure in the Mediterranean gene pool associated with geographic patterns and phenotypic variables . Front. Genet. 10, 872 (2019).

[5] Hamadeh, B., Chalak, L., Coppens d’Eeckenbrugge, G., Benoit, L. & Joly, H. I. Evolution of almond genetic diversity and farmer practices in Lebanon: impacts of the diffusion of a graft-propagated cultivar in a traditional system based on seed-propagation. BMC Plant Biol. 18, 155 (2018).

[6] Rigoldi, M. P., Rapposelli, E., De Giorgio, D., Resta, P. & Porceddu, A. Genetic diversity in two Italian almond collections. Electron. J. Biotechn. 18, 40-45 (2015).

[7] Fernández i Martí, A. et al. Molecular analyses of evolution and population structure in a worldwide almond [ Prunus dulcis (Mill.) D.A. Webb syn. P. amygdalus Batsch] pool assessed by microsatellite markers . Genet. Resour. Crop. Ev. 62, 205-219 (2015).

[8] Font i Forcada, C. et al. Identification of genetic loci associated with quality traits in almond via association mapping. PLoS ONE 10, e0127656 (2015).

[9] Halász, J. et al. Genetic variability is preserved among strongly differentiated and geographically diverse almond germplasm: an assessment by simple sequence repeat markers. Tree Genet. Genomes 15, 12 (2019).

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

[11] Taranto, F., Nicolia, A., Pavan, S., De Vita, P. & D’Agostino, N. Biotechnological and digital revolution for climate-smart plant breeding. Agron 8, 277 (2018).

[12] Pavan, S. et al. Recommendations for choosing the genotyping method and best practices for quality control in crop genome-wide association studies. Front. Genet. https://doi.org/10.3389/fgene.2020.00447 (2020).

[13] Charlesworth, D. & Willis, J. H. The genetics of inbreeding depression. Nat. Rev. Genet. 10, 783-796 (2009).

[14] Lansari, A., Iezzoni, A. F. & Kester, D. E. Morphological variation within collections of Moroccan almond clones and Mediterranean and North American cultivars. Euphytica 78, 27-41 (1994).

[15] Alonso Segura, J. M. & Socias i Company, R. Negative inbreeding effects in tree fruit breeding: self-compatibility transmission in almond. Euphytica 115, 151-158 (2007).

[16] Martínez-García, P. J., Dicenta, F. & Ortega, E. Anomalous embryo sac development and fruit abortion caused by inbreeding depression in almond (Prunus dulcis) . Sci. Hortic. 133, 23-30 (2012).

[17] 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).

[18] Conomos, M. P., Reiner, A. P., Weir, B. S. & Thornton, T. A. Model-free estimation of recent genetic relatedness. Am. J. Hum. Genet. 98, 127-148 (2016).

[19] Kardos, M., Luikart, G. & Allendorf, F. W. Measuring individual inbreeding in the age of genomics: marker-based measures are better than pedigrees. Heredity 115, 63-72 (2015).

[20] Kardos, M., Taylor, H. R., Ellegren, H., Luikart, G. & Allendorf, F. W. Genomics advances the study of inbreeding depression in the wild. Evol. Appl. 9, 1205-1218 (2016).

[21] Goszczynski, D. et al. Runs of homozygosity in a selected cattle population with extremely inbred bulls: descriptive and functional analyses revealed highly variable patterns. PLoS ONE 13, e0200069 (2018).

[22] Yengo, L., Wray, N. R. & Visscher, P. M. Extreme inbreeding in a European ancestry sample from the contemporary UK population. Nat. Commun. 10, 3719 (2019).

[23] Wu, G. A. et al. Genomics of the origin and evolution of Citrus . Nat 554, 311-316 (2018).

[24] Fernández i Martí, A., Font i Forcada, C. & Socias i Company, R. Genetic analysis for physical nut traits in almond. Tree Genet. Genomes 9, 455-465 (2013).

[25] Sorkheh, K. et al. Correlations between quantitative tree and fruit almond traits and their implications for breeding. Sci. Hortic. 125, 323-331 (2010).

[26] Spiegel-Roy, P. & Kochba, J. Inheritance of nut and kernel traits in almond (Prunus amygdalus Batsch) . Euphytica 30, 167-174 (1981).

[27] Dicenta, F., Garcia, J. E. & Carbonell, E. A. Heritability of fruit characters in almond. J. Hort. Sci. 68, 121-126 (1993).

[28] Sánchez-Pérez, R. et al. Mutation of a bHLH transcription factor allowed almond domestication. Sci 364, 1095-1098 (2019).

[29] Alioto, T. et al. Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence. Plant J. 101, 455-472 (2020).

[30] Bazakos, C., Hanemian, M., Trontin, C., Jiménez-Gómez, J. M. & Loudet, O. New strategies and tools in quantitative genetics: how to go from the phenotype to the genotype. Annu. Rev. Plant Biol. 68, 435-455 (2017).

[31] Lander, E. S. & Botstein, D. Homozygosity mapping: a way to map human recessive traits with the DNA of inbred children. Science 236, 1567-1570 (1987).

[32] Purfield, D. C., McParland, S., Wall, E. & Berry, D. P. The distribution of runs of homozygosity and selection signatures in six commercial meat sheep breeds. PLoS ONE 12, e0176780 (2017).

[33] Keller, M. C. et al. Runs of homozygosity implicate autozygosity as a schizophrenia risk factor. PLoS Genet. 8, e1002656 (2012).

[34] Kim, E.-S. et al. Effect of artificial selection on runs of homozygosity in U.S. Holstein cattle. PLoS ONE 8, e80813 (2013).

[35] Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655-1664 (2009).

[36] Pickrell, J. K. & Pritchard, J. K. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 8, e1002967 (2012).

[37] Socias i Company, R. & Gradziel, T. M. (eds) Almonds: Botany, Production and Uses (CABI Publishing, 2017).

[38] Marchese, A., Bošković, R. I., Martínez‐García, P. J. & Tobutt, K. R. The origin of the self‐compatible almond ‘Supernova’. Plant Breed. 127, 105-107 (2008).

[39] Lawson, D. J., van Dorp, L. & Falush, D. A tutorial on how not to over-interpret STRUCTURE and ADMIXTURE bar plots. Nat. Commun. 9, 3258 (2018).

[40] Wickson, E. J. (ed.) The California Fruits and How to Grow Them (Dewey & Co., 1889).

[41] Wood, M. N. (ed.) Almond Varieties in the United States (U.S. Department of Agriculture, 1925).

[42] Kester, D. E. Almond cultivar and breeding programs in California. Acta Hortic. 373, 13-28 (1994).

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

[44] Kester, D. E., Gradziel, T. M. & Grasselly, C. ALMOND (PRUNUS). Acta Hortic. 290, 699-758 (1991).

[45] Jules, J. & Moore, J. N. (eds) Fruit Breeding, Volume III: Nuts (John Wiley & Sons, Inc., 1996).

[46] Richter, A. A. Almond. Trudy gosudarstvennogo Nikitskogo botanicheskogo sada [Proc. of the State Nikita Botanical Garden] 47 [in Russian] (1972).

[47] Denisov, V. P. Almond genetic resources in the USSR and their use in production and breeding. Acta Hortic. 224, 299-306 (1988).

[48] Xu, L. et al. Genomic patterns of homozygosity in Chinese local cattle. Sci. Rep. 9, 16977 (2019).

[49] Song, W. et al. Rice PCR1 influences grain weight and Zn accumulation in grains. Plant Cell Environ. 38, 2327-2339 (2015).

[50] Collins, C., Dewitte, W. & Murray, J. A. H. D-type cyclins control cell division and developmental rate during Arabidopsis seed development. J. Exp. Bot. 63, 3571-3586 (2012).

[51] Liu, S. et al. Tomato AUXIN RESPONSE FACTOR 5 regulates fruit set and development via the mediation of auxin and gibberellin signaling. Sci. Rep. 8, 2971 (2018).

[52] Rodriguez, C. E. et al. Peach fruit development: a comparative proteomic study between endocarp and mesocarp at very early stages underpins the main differential biochemical processes between these tissues. Front. Plant Sci. 10, 715 (2019).

[53] Lombardo, V. A. et al. Metabolic profiling during peach fruit development and ripening reveals the metabolic networks that underpin each developmental stage. Plant Physiol. 157, 1696-1710 (2011).

[54] Liu, J. J., Sturrock, R. & Ekramoddoullah, A. K. The superfamily of thaumatin-like proteins: its origin, evolution, and expression towards biological function. Plant Cell Rep. 29, 419-436 (2010).

[55] Anžlovar, S. & Dermastia, M. The comparative analysis of osmotins and osmotin-like PR-5 proteins. Plant Biol. 5, 116-124 (2003).

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

[57] Kolde, R. Pheatmap: pretty heatmaps. R package version 1.0.12. http://CRAN.R-project.org/package=pheatmap (2019).

[58] Wickham, H. ggplot2: Elegant Graphics for Data Analysis. https://ggplot2.tidyverse.org (2016).

[59] Chang, C. C. et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience 4, 1-16 (2015).

[60] Hill, W. G. & Weir, B. S. Variances and covariances of squared linkage disequilibria in finite populations. Theor. Popul. Biol. 33, 54-78 (1988).

[61] Bates, D., Maechler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1-48 (2015).

[62] Zhang, L. et al. cgaTOH: extended approach for identifying tracts of homozygosity. PLoS ONE 8, e57772 (2013).

[63] Jung, S. et al. GDR (Genome Database for Rosaceae): integrated web-database for Rosaceae genomics and genetics data. Nucleic Acids Res. 36, D1034-D1040 (2008).

[64] Kang, H. M. et al. Variance component model to account for sample structure in genome-wide association studies. Nat. Genet. 42, 348-354 (2010).