Inference of the genetic basis of fruit texture in highbush blueberries using genome-wide association analyses

Abstract

The global production and consumption of blueberry (Vaccinium spp.), a specialty crop known for its abundant bioactive and antioxidant compounds, has more than doubled over the last decade. To hold this momentum, plant breeders have begun to use quantitative genetics and molecular breeding to guide their decisions and select new cultivars that are improved for fruit quality. In this study, we leveraged our inferences on the genetic basis of fruit texture and chemical components by surveying large breeding populations from northern highbush blueberries (NHBs) and southern highbush blueberries (SHBs), the two dominant cultivated blueberries. After evaluating 1065 NHB genotypes planted at the Oregon State University,and 992 SHB genotypes maintained at the University of Florida for 17 texture-related traits, evaluated over multiple years, our contributions consist of the following: (i) we drew attention to differences between NHB and SHB materials and showed that both blueberry types can be differentiated using texture traits; (ii) we computed genetic parameters and shed light on the genetic architecture of important texture attributes, indicating that most traits had a complex nature with low to moderate heritability; (iii) using molecular breeding, we emphasized that prediction could be performed across populations; and finally (iv) the genomic association analyses pinpointed some genomic regions harboring potential candidate genes for texture that could be used for further validation studies. Altogether, the methods and approaches used here can guide future breeding efforts focused on maximizing texture improvements in blueberries.

Acknowledgements

This work was supported by SCRI-NIFA Award 2019-51181-30015, VacciniumCAP: Leveraging genetic and genomic resources to enable development of blueberry and cranberry cultivars with improved fruit quality attributes and by the University of Florida royalty fund generated by the licensing of blueberry cultivars.

Data availability

The genomic and phenotypic data are available as supplemental material.

Conflict of interest statement

The authors declare no conflict of interest.

Supplementary data

Supplementary data is available at Horticulture Research online.

References

[1] Kalt W, Cassidy A, Howard LR. et al. Recent research on the health benefits of blueberries and their anthocyanins. Adv Nutr. 2020; 11 :224-36

[2] Hancock JF, Lyrene P, Finn CE. et al. Blueberries and cranberries. In: Temperate Fruit Crop Breeding. Springer, 2008,115-50

[3] Retamales JB, Hancock JF. Blueberries. Cabi. 2018.

[4] Gallardo RK, Zhang Q, Dossett M. et al. Breeding trait priorities of the blueberry industry in the United States and Canada. HortScience. 2018; 53 :1021-8

[5] Roth M, Muranty H, Di Guardo M. et al. Genomic prediction of fruit texture and training population optimization towards the application of genomic selection in apple. Hortic Res. 2020; 7 :7

[6] Rivera S, Giongo L, Cappai F. et al. Blueberry firmness-a review of the textural and mechanical properties used in quality evaluations. Postharvest Biol Technol. 2022; 192 :112016

[7] Giongo L, Poncetta P, Loretti P. et al. Texture profiling of blueberries (Vaccinium spp.) during fruit development, ripening and storage. Postharvest Biol Technol. 2013; 76 :34-9

[8] Oh H, Pottorff M, Giongo L. et al. Dissecting the texture characteristics determining the phenotypic variation and postharvest change of fruit texture in blueberry. In: Hortscience. ALEXANDRIA, VA, 2022,S130-0

[9] Ferrão LFV, Amadeu RR, Benevenuto J. et al. Genomic selection in an outcrossing autotetraploid fruit crop: lessons from blueberry breeding. Front. Plant Sci. 2021; 12 :1075

[10] Ferrão LFV, Benevenuto J, Oliveira IB. et al. Insights into the genetic basis of blueberry fruit-related traits using diploid and polyploid models in a GWAS context. Front Ecol Evol. 2018; 6 :107

[11] Giongo L, Ajelli M, Pottorff M. et al. Comparative multi-parameters approach to dissect texture subcomponents of highbush blueberry cultivars at harvest and postharvest. Postharvest Biol Technol. 2022; 183 :111696

[12] Aleza P, Boutsika A, Bassil NV. et al. Development of a targeted genotyping platform for reproducible results within tetraploid and hexaploid blueberry. Frontiers in Horticulture 2024; 2 :1339310

[13] Mrode RA. Linear Models for the Prediction of Animal Breeding Values. Cabi; 2014:

[14] Cellon C, Amadeu RR, Olmstead JW. et al. Estimation of genetic parameters and prediction of breeding values in an autotetraploid blueberry breeding population with extensive pedigree data. Euphytica. 2018; 214 :1-13

[15] Luby JJ, Finn CE. Quantitative inheritance of plant growth habit in blueberry progenies. J Am Soc Hortic Sci. 1986; 111 :609-11

[16] Oh H, Pottorff M, Giongo L. et al. Exploring shelf-life predictability of appearance traits and fruit texture in blueberry. Postharvest Biol Technol. 2024; 208 :112643

[17] Cappai F, Benevenuto J, Ferrão L. et al. Molecular and genetic bases of fruit firmness variation in blueberry—a review. Agronomy. 2018; 8 :174

[18] Shi Y, Li B, Su G. et al. Transcriptional regulation of fleshy fruit texture. J Integr Plant Biol. 2022; 64 :1649-72

[19] Muñoz-Bertomeu J, Miedes E, Lorences EP. Expression of xyloglucan endotransglucosylase/hydrolase (XTH) genes and XET activity in ethylene treated apple and tomato fruits. J Plant Physiol. 2013; 170 :1194-201

[20] Zhai Z, Feng C, Wang Y. et al. Genome-wide identification of the xyloglucan endotransglucosylase/hydrolase (XTH) and polygalacturonase (PG) genes and characterization of their role in fruit softening of sweet cherry. Int J Mol Sci. 2021; 22 :12331

[21] Wang D, Yeats TH, Uluisik S. et al. Fruit softening: revisiting the role of pectin. Trends Plant Sci. 2018; 23 :302-10

[22] Uluisik S, Chapman NH, Smith R. et al. Genetic improvement of tomato by targeted control of fruit softening. Nat Biotechnol. 2016; 34 :950-2

[23] Phan TD, Bo W, West G. et al. Silencing of the major salt-dependent isoform of pectinesterase in tomato alters fruit softening. Plant Physiol. 2007; 144 :1960-7

[24] Gao Y, Zhang Y-P, Fan Z-Q. et al. Mutagenesis of SlNAC4 by CRISPR/Cas9 alters gene expression and softening of ripening tomato fruit. Textes anciens. 2021; 1 :1-12

[25] Manning K, Tör M, Poole M. et al. A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nat Genet. 2006; 38 :948-52

[26] Li J, Guo T, Liang M. Exogenous BR delayed peach fruit softening by inhibiting pectin degradation enzyme genes. Front Plant Sci. 2023; 14 :1226921

[27] Lai X, Zhu X, Chen H. et al. The MaC2H2-like zinc finger protein is involved in ripening and ripening disorders caused by chilling stress via the regulation of softening-related genes in ‘Fenjiao’banana. Postharvest Biol Technol. 2022; 186 :111817

[28] Saracoglu O, Ozturk B, Yildiz K. et al. Pre-harvest methyl jasmonate treatments delayed ripening and improved quality of sweet cherry fruits. Sci Hortic. 2017; 226 :19-23

[29] Wei J, Wen X, Tang L. Effect of methyl jasmonic acid on peach fruit ripening progress. Sci Hortic (Amsterdam). 2017; 220 :206-13

[30] Zarid M, García-Carpintero V, Esteras C. et al. Transcriptomic analysis of a near-isogenic line of melon with high fruit flesh firmness during ripening. J Sci Food Agric. 2021; 101 :754-77

[31] Zhao Y, Christensen SK, Fankhauser C. et al. A role for flavin monooxygenase-like enzymes in auxin biosynthesis. Science. 2001; 291 :306-9

[32] Tatsuki M, Soeno K, Shimada Y. et al. Insertion of a transposon-like sequence in the 5 ′-flanking region of the YUCCA gene causes the stony hard phenotype. Plant J. 2018; 96 :815-27

[33] Deng H, Chen Y, Liu Z. et al. SlERF. F12 modulates the transition to ripening in tomato fruit by recruiting the co-repressor TOPLESS and histone deacetylases to repress key ripening genes. Plant Cell. 2022; 34 :1250-72

[34] Guo J-E, Hu Z, Li F. et al. Silencing of histone deacetylase SlHDT3 delays fruit ripening and suppresses carotenoid accumulation in tomato. Plant Sci. 2017; 265 :29-38

[35] Lin H, Ma L, Guo Q. et al. Berry texture QTL and candidate gene analysis in grape (Vitis vinifera L.). Hortic Res. 2023; 10 :uhad226

[36] Ma L, Sun L, Guo Y. et al. Transcriptome analysis of table grapes (Vitis vinifera L.) identified a gene network module associated with berry firmness. PLoS One. 2020; 15 :e0237526

[37] Rothkegel K, Espinoza A, Sanhueza D. et al. Identification of DNA methylation and transcriptomic profiles associated with fruit Mealiness in Prunus persica (L.) Batsch. Front Plant Sci. 2021; 12 :684130

[38] Mengist MF, Bostan H, De Paola D. et al. Autopolyploid inheritance and a heterozygous reciprocal translocation shape chromosome genetic behavior in tetraploid blueberry (Vaccinium corymbosum ). New Phytol. 2023; 237 :1024-39

[39] Lee W-P, Stromberg MP, Ward A. et al. MOSAIK: a hash-based algorithm for accurate next-generation sequencing short-read mapping. PLoS One. 2014; 9 :e90581

[40] Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. arXiv Prepr arXiv12073907. 2012;

[41] Danecek P, Auton A, Abecasis G. et al. The variant call format and VCFtools. Bioinformatics. 2011; 27 :2156-8

[42] Gerard D, Ferrão LFV, Garcia AAF. et al. Genotyping polyploids from messy sequencing data. Genetics. 2018; 210 :789-807

[43] de Bem OI, Resende MFR, Ferrão LFV. et al. Genomic prediction of autotetraploids; influence of relationship matrices, allele dosage, and continuous genotyping calls in phenotype prediction. G3 (Bethesda). 2019; 9 :1189-98

[44] Wright S. The genetical structure of populations. Ann Eugenics. 1949; 15 :323-54

[45] Amadeu RR, Garcia AAF, Munoz PR. et al. AGHmatrix: genetic relationship matrices in R. Bioinformatics. 2023; 39 :btad445

[46] Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000; 155 :945-59

[47] Earl DA, VonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012; 4 :359-61

[48] Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008; 25 :1-18

[49] Breiman L. Bagging predictors. Mach Learn. 1996; 24 :123-40

[50] Butler DG, Cullis BR, Gilmour AR. et al. { ASReml } -R reference manual. 2009

[51] Hernandez CO, Wyatt LE, Mazourek MR. Genomic prediction and selection for fruit traits in winter squash. G three. 2020; 10 :3601-10

[52] Yang J, Zaitlen NA, Goddard ME. et al. Advantages and pitfalls in the application of mixed-model association methods. Nat Genet. 2014; 46 :100-6

[53] Törönen P, Medlar A, Holm L. PANNZER2: a rapid functional annotation web server. Nucleic Acids Res. 2018; 46 :W84-8

[54] Cantalapiedra CP, Hernández-Plaza A, Letunic I. et al. eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol Biol Evol. 2021; 38 :5825-9