High-density linkage map construction and identification of loci regulating fruit quality traits in blueberry

Abstract

Fruit quality traits play a significant role in consumer preferences and consumption in blueberry (Vaccinium corymbosum L). The objectives of this study were to construct a high-density linkage map and to identify the underlying genetic basis of fruit quality traits in blueberry. A total of 287 F₁ individuals derived from a cross between two southern highbush blueberry cultivars, ‘Reveille’ and ‘Arlen’, were phenotyped over three years (2016–2018) for fruit quality-related traits, including titratable acidity, pH, total soluble solids, and fruit weight. A high-density linkage map was constructed using 17k single nucleotide polymorphisms markers. The linkage map spanned a total of 1397 cM with an average inter-loci distance of 0.08 cM. The quantitative trait loci interval mapping based on the hidden Markov model identified 18 loci for fruit quality traits, including seven loci for fruit weight, three loci for titratable acidity, five loci for pH, and three loci for total soluble solids. Ten of these loci were detected in more than one year. These loci explained phenotypic variance ranging from 7 to 28% for titratable acidity and total soluble solid, and 8–13% for pH. However, the loci identified for fruit weight did not explain more than 10% of the phenotypic variance. We also reported the association between fruit quality traits and metabolites detected by Proton nuclear magnetic resonance analysis directly responsible for these fruit quality traits. Organic acids, citric acid, and quinic acid were significantly (P < 0.05) and positively correlated with titratable acidity. Sugar molecules showed a strong and positive correlation with total soluble solids. Overall, the study dissected the genetic basis of fruit quality traits and established an association between these fruit quality traits and metabolites.

References

[1] Mengist, M. F. et al. Diversity in metabolites and fruit quality traits in blueberry enables ploidy and species differentiation and establishes a strategy for future genetic studies. Front. Plant Sci. 11, 370 (2020).

[2] Mengist, M. F. et al. Development of a genetic framework to improve the efficiency of bioactive delivery from blueberry. Sci. Rep. 10, 1-13 (2020).

[3] Farneti, B. et al. Development of a novel phenotypic roadmap to improve blueberry quality and storability. Front. Plant Sci. 11, 1140 (2020).

[4] Gilbert, J. L. et al. Identifying breeding priorities for blueberry flavor using biochemical, sensory, and genotype by environment analyses. PLoS ONE 10, 1-21 (2015).

[5] Evans E. A. & Ballen F. H. An overview of US blueberry production, trade, and consumption, with special reference to Florida. University of Florida, IFAS, Extension, FE952 Available online: https://edis.ifas.ufl.edu/publication/fe952 2014.

[6] Yang, B. & Kortesniemi, M. Clinical evidence on potential health benefits of berries. Curr. Opin. Food Sci. 2, 36-42 (2015).

[7] Gallardo, R. K. et al. Breeding trait priorities of the blueberry industry in the United States and Canada. HortScience 53, 1021-1028 (2018).

[8] Brazelton C., Kayla Y. & Bauer N. 2016 Global Blueberry Statistics and Intelligence Report. 2017. https://www.internationalblueberry.org.

[9] Rawandoozi, Z. J. et al. Identification and characterization of QTLs for fruit quality traits in peach through a multi-family approach. BMC genomics 21, 1-18 (2020).

[10] Calle, A. & Wünsch, A. Multiple-population QTL mapping of maturity and fruit-quality traits reveals LG4 region as a breeding target in sweet cherry (Prunus avium L.). Horticulture Res. 7, 1-13 (2020).

[11] Benevenuto, J., Ferrão, L. F. V., Amadeu, R. R. & Munoz, P . How can a high-quality genome assembly help plant breeders? Gigascience 8, giz068 (2019).

[12] Zeng, Q. et al. High altitude is beneficial for antioxidant components and sweetness accumulation of Rabbiteye Blueberry. Front. Plant Sci. 11, 1492 (2020).

[13] Zhang, J. et al. Evaluation of sugar and organic acid composition and their levels in highbush blueberries from two regions of China. J. Integr. Agriculture 19, 2352-2361 (2020).

[14] Colantonio, V. et al. Metabolomic selection for enhanced fruit flavor. bioRxiv 9, 30280, https://doi.org/10.1101/2020.09.17.302802 (2020).

[15] Rowland, L. J. et al. Generating genomic tools for blueberry improvement. Int. J. Fruit. Sci. 12, 276-287 (2012).

[16] Lobos, G. A. & Hancock, J. F. Breeding blueberries for a changing global environment: a review. Front. Plant Sci. 6, 1-14 (2015).

[17] Retamales, J. B. & Hancock, J. F. Blueberries. (CABI, Wallington, UK, 2018).

[18] Rowland, L. J. & Levi, A. RAPD-based genetic linkage map of blueberry derived from a cross between diploid species (Vaccinium darrowi and V. elliottii). Theor. Appl. Genet.: Int. J. Plant Breed. Res. 87, 863-868 (1994).

[19] Hancock, J. F., Lyrene, P., Finn, C. E., Vorsa, N. & Lobos, G. A. Blueberries and cranberries. Temp. Fruit. Crop Breed.: Germplasm Genomics 290, 115-149 (2008).

[20] Rowland, L. J., Ogden, E. L. & Vinyard, B. T. Phenotypic evaluation of a hybrid diploid blueberry population for plant development and fruit quality traits. Agronomy 10, 1067 (2020).

[21] Cappai, F. et al. High-resolution linkage map and QTL analyses of fruit firmness in autotetraploid blueberry. Front. Plant Sci. 11, 1767 (2020).

[22] Mengist, M. F. et al. Genetic mapping of quantitative trait loci for tuber-cadmium and zinc concentration in potato reveals associations with maturity and both overlapping and independent components of genetic control. Theor. Appl. Genet. 131, 929-945 (2018).

[23] Hackett, C. A., Boskamp, B., Vogogias, A., Preedy, K. F. & Milne, I. TetraploidSNPMap: Software for linkage analysis and QTL mapping in autotetraploid populations using SNP dosage data. J. Heredity 108, 438-442 (2017).

[24] Bourke, P. M. et al. polymapR-linkage analysis and genetic map construction from F1 populations of outcrossing polyploids. Bioinformatics 34, 3496-3502 (2018).

[25] Colle, M. et al. Haplotype-phased genome and evolution of phytonutrient pathways of tetraploid blueberry. GigaScience 8, giz012 (2019).

[26] McCallum, S. et al. Construction of a SNP and SSR linkage map in autotetraploid blueberry using genotyping by sequencing. Mol. Breed. 36, 1-24 (2016).

[27] Zhou, Q. et al. Haplotype-resolved genome analyses of a heterozygous diploid potato. Nat. Genet. 52, 1018-1023 (2020).

[28] Ferrão, L. F. V. et al. Insights into the genetic basis of blueberry fruit-related traits using diploid and polyploid models in a GWAS context. Front. Ecol. Evolution 6, 107 (2018).

[29] Wu, B. et al. Genome-wide identification, expression patterns, and functional analysis of UDP glycosyltransferase family in peach (Prunus persica L. Batsch). Front. Plant Sci. 8, 389 (2017).

[30] Hong, Z., Zhang, Z., Olson, J. M. & Verma, D. P. S. A novel UDP-glucose transferase is part of the callose synthase complex and interacts with phragmoplastin at the forming cell plate. Plant Cell 13, 769-779 (2001).

[31] Moremen, K. W. & Haltiwanger, R. S. Emerging structural insights into glycosyltransferase-mediated synthesis of glycans. Nat. Chem. Biol. 15, 853-864 (2019).

[32] Shi, C.-Y. et al. Citrus PH5-like H+-ATPase genes: identification and transcript analysis to investigate their possible relationship with citrate accumulation in fruits. Front. Plant Sci. 6, 135 (2015).

[33] Cohen, S. et al. The PH gene determines fruit acidity and contributes to the evolution of sweet melons. Nat. Commun. 5, 1-9 (2014).

[34] Verma, S. 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 Genetics and Genomes 2019; 15 https://doi.org/10.1007/s11295-019-1324-y.

[35] Fong, S. K. et al. A low malic acid trait in cranberry fruit: genetics, molecular mapping, and interaction with a citric acid locus. Tree Genetics and Genomes 2021; 17 https://doi.org/10.1007/s11295-020-01482-8.

[36] Diaz-Garcia, L. et al. Massive phenotyping of multiple cranberry populations reveals novel QTLs for fruit anthocyanin content and other important chemical traits. Mol. Genet. Genomics 293, 1379-1392 (2018).

[37] Cline, B. 2007 Blueberry Workshop Agent Training Sponsored by: NC Blueberry Council. 2007 https://smallfruits.org/files/2019/06/SRFC-training-Jun-19-21-2007.pdf.

[38] Diaz-Garcia, L. et al. GiNA, an Efficient and high-throughput software for horticultural phenotyping. PLoS ONE 11, e0160439 (2016).

[39] Ulrich, E. L. et al. BioMagResBank. Nucleic acids Res. 36, D402-D408 (2007).

[40] Sumner, L. W. et al. Proposed minimum reporting standards for chemical analysis. Metabolomics 3, 211-221 (2007).

[41] Wei, T. et al. Package ‘corrplot’. Statistician 56, 316-324 (2017).

[42] Panta, G. R., Rowland, L. J., Arora, R., Ogden, E. L. & Lim, C. C. Inheritance of cold hardiness and dehydrin genes in diploid mapping populations of blueberry. J. Crop Improv. 10, 37-52 (2004).

[43] Bolger, A. et al. The genome of the stress-tolerant wild tomato species. Nat. Genet. 46, 1034-1039 (2014).

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

[45] Garrison E., Marth G. Haplotype-based variant detection from short-read sequencing. arXiv preprint arXiv:12073907 2012.

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

[47] Gerard, D., Ferrão, L. F. V., Garcia, A. A. F. & Stephens, M. Genotyping polyploids from messy sequencing data. Genetics 210, 789-807 (2018).

[48] Preedy, K. F. & Hackett, C. A. A rapid marker ordering approach for high-density genetic linkage maps in experimental autotetraploid populations using multidimensional scaling. Theor. Appl. Genet. 129, 2117-2132 (2016).

[49] Voorrips, R. E. MapChart: software for the graphical presentation of linkage maps and QTLs. J. Heredity 93, 77-78 (2002).

[50] Wickham H. ggplot2: elegant graphics for data analysis. springer, 2016.

[51] Rezvoy, C., Charif, D., Guéguen, L. & Marais, G. A. B. MareyMap: an R-based tool with graphical interface for estimating recombination rates. Bioinformatics 23, 2188-2189 (2007).

[52] Hackett, C. A., Bradshaw, J. E. & Bryan, G. J. QTL mapping in autotetraploids using SNP dosage information. Theor. Appl. Genet. 127, 1885-1904 (2014).

[53] Schwarz, G. Estimating the dimension of a model. Ann. Stat. 6, 461-464 (1978).