Molecular and Genetic Bases of Fruit Firmness Variation in Blueberry—A Review

Genetic and epigenetic signatures for improved breeding of cultivated blueberry

Blueberry belongs to the Vaccinium genus and is a highly popular fruit crop with significant economic importance. It was not until the early twentieth century that they began to be domesticated through extensive interspecific hybridization. Here, we collected 220 Vaccinium accessions from various geographical locations, including 154 from the United States, 14 from China, eight from Australia, and 29 from Europe and other countries, comprising 164 Vaccinium corymbosum, 15 Vaccinium ashei, 10 lowbush blueberries, seven half-high blueberries, and others. We present the whole-genome variation map of 220 accessions and reconstructed the hundred-year molecular history of interspecific hybridization of blueberry. We focused on the two major blueberry subgroups, the northern highbush blueberry (NHB) and southern highbush blueberry (SHB) and identified candidate genes that contribute to their distinct traits in climate adaptability and fruit quality. Our analysis unveiled the role of gene introgression from Vaccinium darrowii and V. ashei into SHB in driving the differentiation between SHB and NHB, potentially facilitating SHB's adaptation to subtropical environments. Assisted by genome-wide association studies, our analysis suggested VcTBL44 as a pivotal gene regulator governing fruit firmness in SHB. Additionally, we conducted whole-genome bisulfite sequencing on nine NHB and 12 SHB cultivars, and characterized regions that are differentially methylated between the two subgroups. In particular, we discovered that the β-alanine metabolic pathway genes were enriched for DNA methylation changes. Our study provides high-quality genetic and epigenetic variation maps for blueberry, which offer valuable insights and resources for future blueberry breeding.

Biochemical dynamics during postharvest: Highlighting the interplay of stress during storage and maturation of fresh produce

The lifecycle of fresh produce involves a sequence of biochemical events during their ontology, and these events are particularly significant for climacteric fruits. A high demand during ripening is observed in these plant products, which is reflected in a high rate of respiration and ethylene production. Increased respiratory demand triggers the activation of secondary pathways such as alternate oxidase, which do not experience critical increases in energy consumption in non-climacteric fruit. In addition, biochemical events produced by external factors lead to compensatory responses in fresh produce to counteract the oxidative stress caused by the former. The dynamics of these responses are accompanied by signaling, where reactive oxygen species play a pivotal role in fresh product cell perception. This review aims to describe the protection mechanisms of fresh produce against environmental challenges and how controlled doses of abiotic stressors can be used to improve quality and prolong their shelf-life through the interaction of stress and defense mechanisms.

Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids

Discovery and analysis of genetic variants underlying agriculturally important traits are key to molecular breeding of crops. Reduced representation approaches have provided cost-efficient genotyping using next-generation sequencing. However, accurate genotype calling from next-generation sequencing data is challenging, particularly in polyploid species due to their genome complexity. Recently developed Bayesian statistical methods implemented in available software packages, polyRAD, EBG, and updog, incorporate error rates and population parameters to accurately estimate allelic dosage across any ploidy. We used empirical and simulated data to evaluate the three Bayesian algorithms and demonstrated their impact on the power of genome-wide association study (GWAS) analysis and the accuracy of genomic prediction. We further incorporated uncertainty in allelic dosage estimation by testing continuous genotype calls and comparing their performance to discrete genotypes in GWAS and genomic prediction. We tested the genotype-calling methods using data from two autotetraploid species, Miscanthus sacchariflorus and Vaccinium corymbosum, and performed GWAS and genomic prediction. In the empirical study, the tested Bayesian genotype-calling algorithms differed in their downstream effects on GWAS and genomic prediction, with some showing advantages over others. Through subsequent simulation studies, we observed that at low read depth, polyRAD was advantageous in its effect on GWAS power and limit of false positives. Additionally, we found that continuous genotypes increased the accuracy of genomic prediction, by reducing genotyping error, particularly at low sequencing depth. Our results indicate that by using the Bayesian algorithm implemented in polyRAD and continuous genotypes, we can accurately and cost-efficiently implement GWAS and genomic prediction in polyploid crops.

Metabolomic and transcriptomic analyses of the flavonoid biosynthetic pathway in blueberry (Vaccinium spp.)

As a highly economic small fruit crop, blueberry is enjoyed by most people in terms of color, taste, and rich nutrition. To better understand its coloring mechanism on the process of ripening, an integrative analysis of the metabolome and transcriptome profiles was performed in three blueberry varieties at three developmental stages. In this study, 41 flavonoid metabolites closely related to the coloring in blueberry samples were analyzed. It turned out that the most differential metabolites in the ripening processes were delphinidin-3-O-arabinoside (dpara), peonidin-3-O-glucoside (pnglu), and delphinidin-3-O-galactoside (dpgal), while the most differential metabolites among different varieties were flavonols. Furthermore, to obtain more accurate and comprehensive transcripts of blueberry during the developmental stages, PacBio and Illumina sequencing technology were combined to obtain the transcriptome of the blueberry variety Misty, for the very first time. Finally, by applying the gene coexpression network analysis, the darkviolet and bisque4 modules related to flavonoid synthesis were determined, and the key genes related to two flavonoid 3', 5'-hydroxylase (F3'5'H) genes in the darkviolet module and one bHLH transcription factor in the bisque4 module were predicted. It is believed that our findings could provide valuable information for the future study on the molecular mechanism of flavonoid metabolites and flavonoid synthesis pathways in blueberries.

An Exploration of Pepino (Solanum muricatum) Flavor Compounds Using Machine Learning Combined with Metabolomics and Sensory Evaluation

Flavor is one of the most important characteristics that directly determines the popularity of a food. Moreover, the flavor of fruits is determined by the interaction of multiple metabolic components. Pepino, an emerging horticultural crop, is popular for its unique melon-like flavor. We analyzed metabolomics data from three different pepino growing regions in Haidong, Wuwei, and Jiuquan and counted the status of sweetness, acidity, flavor, and overall liking ratings of pepino fruit in these three regions by sensory panels. The metabolomics and flavor ratings were also integrated and analyzed using statistical and machine learning models, which in turn predicted the sensory panel ratings of consumers based on the chemical composition of the fruit. The results showed that pepino fruit produced in the Jiuquan region received the highest ratings in sweetness, flavor intensity, and liking, and the results with the highest contribution based on sensory evaluation showed that nucleotides and derivatives, phenolic acids, amino acids and derivatives, saccharides, and alcohols were rated in sweetness (74.40%), acidity (51.57%), flavor (56.41%), and likability (33.73%) dominated. We employed 14 machine learning strategies trained on the discovery samples to accurately predict the outcome of sweetness, sourness, flavor, and liking in the replication samples. The Radial Sigma SVM model predicted with better accuracy than the other machine learning models. Then we used the machine learning models to determine which metabolites influenced both pepino flavor and consumer preference. A total of 27 metabolites most important for pepino flavor attributes to distinguish pepino originating from three regions were screened. Substances such as N-acetylhistamine, arginine, and caffeic acid can enhance pepino‘s flavor intensity, and metabolites such as glycerol 3-phosphate, aconitic acid, and sucrose all acted as important variables in explaining the liking preference. While glycolic acid and orthophosphate inhibit sweetness and enhance sourness, sucrose has the opposite effect. Machine learning can identify the types of metabolites that influence fruit flavor by linking metabolomics of fruit with sensory evaluation among consumers, which conduces breeders to incorporate fruit flavor as a trait earlier in the breeding process, making it possible to select and release fruit with more flavor.

GWAS provides new insights into the genetic mechanisms of phytochemicals production and red skin colour in apple

Understanding the genetic architecture of apple phytochemicals, and their interplay with conventional selection traits, is critical for the development of new apple cultivars with enhanced health benefits. Apple accessions (n = 344) used for this genome-wide association study (GWAS) represented the wide diversity of metabolic profiles in the domesticated and wild Malus genepools. Fruit samples were phenotyped for 34 metabolites, including a stable vitamin C glycoside "ascorbic acid 2-β-glucoside" (AA-2βG), and the accessions were genotyped using the Apple 20 K SNP Array. Several fruit quality traits, including red skin over-colour (OCOL), were also assessed. Wild Malus accessions showed at least 2-fold higher average content of several metabolites (e.g. ascorbic acid, chlorogenic acid, phloridzin, and trilobatin) than Malus domestica accessions. Several new genomic regions and potential candidate genes underpinning the genetic diversity of apple phytochemicals were identified. The percentage of phenotypic variance explained by the best SNP ranged between 3% and 21% for the different metabolites. Novel association signals for OCOL in the syntenic regions on chromosomes 13 and 16 suggested that whole genome duplication has played a role in the evolution of apple red skin colour. Genetic correlations between phytochemicals and sensory traits were moderate. This study will assist in the selection of Malus accessions with specific phytochemical profiles to establish innovative genomics-based breeding strategies for the development of apple cultivars with enhanced nutritional value.

Postharvest High-CO2 Treatments on the Quality of Soft Fruit Berries: An Integrated Transcriptomic, Proteomic, and Metabolomic Approach

Soft fruits are appreciated for their taste qualities and for being a source of health-promoting compounds. However, their postharvest is affected by their high respiratory rates and susceptibility to fungal decay. Our aim here is to provide a perspective on the application of short-term high-CO₂ treatments at a low temperature to maintain the postharvest quality of soft fruits. This work also suggests using a multi-omics approach to better understand the role of the cell wall and phenolic compounds in maintaining quality. Finally, the contribution of high-throughput transcriptomic technologies to understand the mechanisms modulated by the short-term gaseous treatments is also highlighted.

Atypical Climacteric and Functional Ethylene Metabolism and Signaling During Fruit Ripening in Blueberry (Vaccinium sp.)

Climacteric fruits display an increase in respiration and ethylene production during the onset of ripening, while such changes are minimal in non-climacteric fruits. Ethylene is a primary regulator of ripening in climacteric fruits. The ripening behavior and role of ethylene in blueberry (Vaccinium sp.) ripening is controversial. This work aimed to clarify the fruit ripening behavior and the associated role of ethylene in blueberry. Southern highbush (Vaccinium corymbosum hybrids) and rabbiteye (Vaccinium ashei) blueberry displayed an increase in the rate of respiration and ethylene evolution, both reaching a maxima around the Pink and Ripe stages of fruit development, consistent with climacteric fruit ripening behavior. Increase in ethylene evolution was associated with increases in transcript abundance of its biosynthesis genes, AMINOCYCLOPROPANE CARBOXYLATE (ACC) SYNTHASE1 (ACS1) and ACC OXIDASE2 (ACO2), implicating them in developmental ethylene production during ripening. Blueberry fruit did not display autocatalytic system 2 ethylene during ripening as ACS transcript abundance and ACC concentration were not enhanced upon treatment with an ethylene-releasing compound (ethephon). However, ACO transcript abundance was enhanced in response to ethephon, suggesting that ACO was not rate-limiting. Transcript abundance of multiple genes associated with ethylene signal transduction was upregulated concomitant with developmental increase in ethylene evolution, and in response to exogenous ethylene. As these changes require ethylene signal transduction, fruit ripening in blueberry appears to involve functional ethylene signaling. Together, these data indicate that blueberry fruit display atypical climacteric ripening, characterized by a respiratory climacteric, developmentally regulated but non-autocatalytic increase in ethylene evolution, and functional ethylene signaling.

There and back again; historical perspective and future directions for Vaccinium breeding and research studies

The genus Vaccinium L. (Ericaceae) contains a wide diversity of culturally and economically important berry crop species. Consumer demand and scientific research in blueberry (Vaccinium spp.) and cranberry (Vaccinium macrocarpon) have increased worldwide over the crops' relatively short domestication history (~100 years). Other species, including bilberry (Vaccinium myrtillus), lingonberry (Vaccinium vitis-idaea), and ohelo berry (Vaccinium reticulatum) are largely still harvested from the wild but with crop improvement efforts underway. Here, we present a review article on these Vaccinium berry crops on topics that span taxonomy to genetics and genomics to breeding. We highlight the accomplishments made thus far for each of these crops, along their journey from the wild, and propose research areas and questions that will require investments by the community over the coming decades to guide future crop improvement efforts. New tools and resources are needed to underpin the development of superior cultivars that are not only more resilient to various environmental stresses and higher yielding, but also produce fruit that continue to meet a variety of consumer preferences, including fruit quality and health related traits.

De novo sequencing and analysis of the transcriptome of two highbush blueberry (Vaccinium corymbosum L.) cultivars 'Bluecrop' and 'Legacy' at harvest and following post-harvest storage

Fruit firmness and in particular the individual components of texture and moisture loss, are considered the key quality traits when describing blueberry fruit quality, and whilst these traits are genetically regulated, the mechanisms governing their control are not clearly understood. In this investigation, RNAseq was performed on fruits of two blueberry cultivars with very different storage properties, 'Bluecrop' and 'Legacy', at harvest, three weeks storage in a non-modified environment at 4 °C and after three weeks storage at 4 °C followed by three days at 21 °C, with the aim of understanding the transcriptional changes that occur during storage in cultivars with very different post-harvest fruit quality. De novo assemblies of the transcriptomes of the two cultivars were performed separately and a total of 39,335 and 41,896 unigenes for 'Bluecrop' and 'Legacy' respectively were resolved. Differential gene expression analyses were grouped into four cluster profiles based on changes in transcript abundance between harvest and 24 days post-harvest. A total of 290 unigenes were up-regulated in 'Legacy' only, 685 were up-regulated in 'Bluecrop', 252 were up-regulated in both cultivars and 948 were down-regulated in both cultivars between harvest and 24 days post-harvest. Unigenes showing significant differential expression between harvest and following post-harvest cold-storage were grouped into classes of biological processes including stress responses, cell wall metabolism, wax metabolism, calcium metabolism, cellular components, and biological processes. In total 21 differentially expressed unigenes with a putative role in regulating the response to post-harvest cold-storage in the two cultivars were identified from the de novo transcriptome assemblies performed. The results presented provide a stable foundation from which to perform further analyses with which to functionally validate the candidate genes identified, and to begin to understand the genetic mechanisms controlling changes in firmness in blueberry fruits post-harvest.

Cell Wall Calcium and Hemicellulose Have a Role in the Fruit Firmness during Storage of Blueberry (Vaccinium spp.)

The firmness of blueberry is one of its most significant quality attributes. Modifications in the composition of the cell wall have been associated with changes in the fruit firmness. In this work, cell wall components and calcium concentration in two blueberry cultivars with contrasting firmness phenotypes were evaluated at harvest and 30 days cold storage (0 °C). High performance anion-exchange chromatography with pulse amperometric detector (HPAEC-PAD) analysis was performed using the "Emerald" (firmer) and "Jewel" (softer) blueberry cultivars, showing increased glucose in the firmer cultivar after cold storage. Moreover, the LM15 antibody, which recognizes xyloglucan domains, displayed an increased signal in the Emerald cultivar after 30 d cold storage. Additionally, the antibody 2F4, recognizing a homogalacturonan calcium-binding domain, showed a greater signal in the firmer Emerald blueberries, which correlates with a higher calcium concentration in the cell wall. These findings suggest that xyloglucan metabolism and a higher concentration of cell wall calcium influenced the firmness of the blueberry fruit. These results open new perspectives regarding the role of cell wall components as xyloglucans and calcium in blueberry firmness.

Genomic insight into the developmental history of southern highbush blueberry populations

Interspecific hybridization is a common breeding approach for introducing novel traits and genetic diversity to breeding populations. Southern highbush blueberry (SHB) is a blueberry cultivar group that has been intensively bred over the last 60 years. Specifically, it was developed by multiple interspecific crosses between northern highbush blueberry [NHB, Vaccinium corymbosum L. (2n = 4x = 48)] and low-chill Vaccinium species to expand the geographic limits of highbush blueberry production. In this study, we genotyped polyploid blueberries, including 105 SHB, 17 NHB, and 10 rabbiteye blueberry (RE) (Vaccinium virgatum Aiton), from the accessions planted at Poplarville, Mississippi, and accessions distributed in Japan, based on the double-digest restriction site-associated DNA sequencing. The genome-wide SNP data clearly indicated that RE cultivars were genetically distinct from SHB and NHB cultivars, whereas NHB and SHB were genetically indistinguishable. The population structure results appeared to reflect the differences in the allele selection strategies that breeders used for developing germplasm adapted to local climates. The genotype data implied that there are no or very few genomic segments that were commonly introgressed from low-chill Vaccinium species to the SHB genome. Principal component analysis-based outlier detection analysis found a few loci associated with a variable that could partially differentiate NHB and SHB. These SNP loci were detected in Mb-scale haplotype blocks and may be close to the functional genes related to SHB development. Collectively, the data generated in this study suggest a polygenic adaptation of SHB to the southern climate, and may be relevant for future population-scale genome-wide analyses of blueberry.

High-Resolution Linkage Map and QTL Analyses of Fruit Firmness in Autotetraploid Blueberry

Blueberry (Vaccinium corymbosum and hybrids) is an autotetraploid crop whose commercial relevance has been growing steadily during the last 20 years. However, the ever-increasing cost of labor for hand-picking blueberry is one main constraint in competitive marketing of the fruit. Machine harvestability is, therefore, a key trait for the blueberry industry. Understanding the genetic architecture of traits related to machine harvestability through Quantitative Trait Loci (QTL) mapping is the first step toward implementation of molecular breeding for faster genetic gains. Despite recent advances in software development for autotetraploid genetic mapping, a high-resolution map is still not available for blueberry. In this study, we crafted a map for autotetraploid low-chill highbush blueberry containing 11,292 SNP markers and a total size of 1,953.97 cM (average density of 5.78 markers/cM). This map was subsequently used to perform QTL analyses in 2-year field trials for a trait crucial to machine harvesting: fruit firmness. Preliminary insights were also sought for single evaluations of firmness retention after cold storage, and fruit detachment force traits. Significant QTL peaks were identified for all the traits and overlapping QTL intervals were detected for firmness across the years. We found low-to-moderate QTL effects explaining the phenotypic variance, which suggest a quantitative nature of these traits. The QTL intervals were further speculated for putative gene repertoire. Altogether, our findings provide the basis for future fine-mapping and molecular breeding efforts for machine harvesting in blueberry.