Evaluating Wild Germplasm Introgression into Autotetraploid Blueberry

: Wild germplasm can be classiﬁed as the raw material essential for crop improvement. Introgression of wild germplasm is normally used in breeding to increase crop quality or resilience to evolving biotic and abiotic threats. Here, we explore the potential of introgressing Vaccinium elliottii into commercial blueberry germplasm. Vaccinium elliottii is a wild diploid blueberry species endemic to the southeastern United States that possesses highly desirable and economically important traits for blueberry breeding such as: short bloom to ripe period, adaptation to upland sandy soils, disease resistance, ﬁrmness, and pleasant ﬂavor. To examine the potential of hybridization, we evaluated populations of interspeciﬁc hybrids across multiple stages of breeding (i.e., F 1 , F 2 , and backcrosses) in two crop seasons. We used our extensive pedigree data to generate breeding values for pre-breeding blueberry hybrid populations. Hybrid performance was evaluated considering ﬁtness (i.e., plant vigor and plant height) in addition to evaluating six fruit-quality and marketable-related traits (i.e., size, ﬁrmness, acidity, soluble solids, weight, and yield). Overall, F 2 and backcrosses rapidly achieved market thresholds, presenting values not signiﬁcantly different from commercial blueberry germplasm. Our results conﬁrmed the potential of exploiting the high genetic variability contained in V. elliottii for interspeciﬁc hybridization. Additionally, we developed germplasm resources that can be further evaluated and utilized in the breeding process, advancing selections for fruit quality and environmental adaptation. were generally preferred by consumers. Our results show that no signiﬁcant differences for fruit


Introduction
Global blueberry production has more than tripled since 1985 [1,2]. Thus, the importance of blueberry as a crop has increased significantly [3], requiring continuous improvements to meet market needs. Breeding provides several possibilities to address market requirements [4,5]. One of these possibilities is the identification and development of potential germplasm resources such as native wild germplasm, which is of great importance for crop breeding and pre-breeding as well as for global food security [4,6,7].
A successful example of the use of wild species as genetic resources for blueberry breeding is the development of southern highbush blueberry (SHB) cultivars [3,4,8]. According to Sharpe [9], the three main species that were used to create the SHB germplasm were Northern highbush blueberry (V. corymbosum), Darrow's blueberry (V. darrowii), and Rabbiteye blueberry (V. virgatum). Collectively, the introgression of wild germplasm generated SHB hybrids with low chilling hours, early ripening in berries, increase in fruit quality and productivity, improved tolerance to higher soil pH, and added resistance to multiple diseases [4,10,11]. Current breeding efforts focus on introgressing further disease resistance, abiotic stress adaptation, earlier ripening, improved fruit quality traits such as flavor and texture, and traits conducive to machine harvesting [2,11].
When choosing wild species to introgress into a crop's commercial germplasm, breeding programs should prioritize species with a broad spectrum of desirable traits that will complement and improve the existing germplasm. Elliott's Blueberry (Vaccinium elliottii Chapm.) exhibits several beneficial fruit quality traits such as preferable flavor and fragrance, minimal scarring, small seeds, and juiciness [4,11,12]. Additionally, V. elliottii displays economically important traits for blueberry production, such as tolerance to drought and high pH, as well as presenting early ripening [11]. The species also presents high firmness [12], which is a valuable trait for developing machine harvestable varieties, and a critical breeding focus to decrease overall cost of production [13].
Moreover, as crop pest pressures continue to increase, research into wild species' resistance must continue to meet the escalating pressures that threaten global food production. Notably, wild species can allow for introgression of disease-resistance traits as previously seen in wheat [14], rice [15], and corn [16]. In addition to improved fruit quality and abiotic stress resistance, V. elliottii can also function as a genetic resource for biotic stress resistance. In this sense, V. elliottii introgression could enhance the SHB germplasm with fungal resistance to Phytophthora root rot, blueberry stem canker, and stem blight, as well as entomological resistance to the sharp-nosed leaf hopper [11]. Several attempts to cross V. elliottii wild genotypes into SHB germplasm have been made in the past [12,17,18]. The crosses' success led to the development of highly aromatic, flavorful, and commercially successful cultivars, such as "Carteret," "Snowchaser," and "Kestrel" [12,19]. Norden et al. [12] and Dweikat et al. [17] also obtained vigorous fertile hybrids when crossing V. elliottii into SHB germplasm. These studies generated valuable resources for the introgression of genetic diversity and beneficial traits into breeding populations. Undeniably, wild germplasm has proven to be a valuable historical and modern tool for blueberry breeding programs. More research is required, however, to fully explore the benefits of V. elliottii as a genetic resource for breeding programs, given V. elliottii's intraspecific diversity [8,11,12], and its ability to produce quality hybrids [12,17,18]. Wild introgression from a species such as V. elliottii can allow breeding programs to maintain a dynamic and diverse germplasm to meet evolving environmental and economic demands of the blueberry industry.
Here, we evaluate the introgression of V. elliottii into SHB germplasm. We evaluated progeny's fruit quality and overall vigor traits through multiple crosses between V. elliottii and SHB over two crop seasons. In order to maximize the production of fertile progeny, tetraploid clones of V. elliottii were used. The interspecifically generated populations under investigation comprised different stages of the early breeding process, including the F 1 , F 2 , and backcrosses. The information generated here can help guide genotype selection toward the development of new improved blueberry genotypes.

Plant Material
Five breeding stages of the hybridization process between V. elliottii and V. corymbosum were evaluated in this study (i.e., F 1 , F 2 , BC 1 , wild V. elliottii, and SHB). In this paper, pseudo-F 2 and pseudo-backcrossing schemes were designed to circumvent the high inbreeding depression observed for blueberry. In summary, F 1 families were obtained in crosses between tetraploid V. elliottii and selected SHB genotypes, pseudo-F 2 families were the product of crossing two different F 1 hybrid genotypes, and pseudo-BC 1 families were the product of backcrossing F 1 hybrids to a non-related SHB genotype ( Figure 1A). For comparison with parental genotypes, we also evaluated wild V. elliottii families obtained from open-pollinated crosses between wild diploid V. elliottii genotypes and SHB commercial germplasm families. To facilitate understanding, we will omit the expression "pseudo" when talking about each breeding stage.
Two breeding sets were evaluated in this study. The seeds used to compose the first population, hereafter called 2017 nursery, were generated by Norden [12], encompassing two F 2 families and one BC 1 hybrid family. For comparison, data from two southern highbush (SHB) families and one V. elliottii family established in the same nursery were collected. The seeds used to compose the second population, hereafter called 2018 nursery, were generated by Lyrene [20], which included individuals from eight F 1 families, three F 2 families, and two BC 1 families. Additionally, for comparison, data from four SHB families and two V. elliottii families established in the same nursery were collected. Individuals from V. elliottii used in this study were obtained from unreduced gametes and colchicine-treated plants (Supplemental Tables S1 and S2). Tetraploidy of these individuals was confirmed by flow cytometry and visual analysis of the pollen (as described by [12,20]).
As described by Norden [12] and Lyrene [20], these populations were generated through controlled crosses that were performed by emasculating the female recipient before their flowers open, followed by manual application of the pollen to the stigma. Pollen samples were manually collected from selected male parents. Seeds obtained from each cross were germinated in individual 2 L pots. Plants from the 2017 nursery families were transplanted into a high-density blueberry nursery in June 2017, while plants from the 2018 nursery were transplanted into a second high-density blueberry nursery in June 2018. Both nurseries were located at University of Florida Research Station in Citra, Florida (29 • 24 24.18" N, 82 • 8 29.53" W). For both nurseries, each family contained approximately 100 seedlings. Genotypes were spaced at 15 cm between plants and 40 cm between rows. Standard cultivation procedures were applied, including irrigation at a flow of 10-17 L per minute, administered three times a week for 1.5 h with an overhead sprinkler irrigation system as needed, to prevent drought stress. Soil was tilled and amended with pine bark before planting. Slow-release fertilizer was applied approximately monthly from April to June (in mg/L): 15 N, 5 P 2 O 5 , 10 K 2 O (Blueberry Mix, Growers Fertilizer Corporation, Lake Alfred, FL, USA). Manual weed control was performed. At the beginning of the second year, coarse pine bark was applied as a mulch in each nursery to retain moisture and control weeds. The experimental unit of this study consisted of breeding stages (i.e., F 1 , F 2 , BC 1 , V. elliottii, and SHB) sampled within each nursery that were distributed in a complete randomized design.

Phenotypic Evaluations
During two crop seasons (2017 and 2018 nurseries), information was collected for nine traits relating to fruit quality and plant fitness. The traits evaluated were yield (g), fruit weight (g), fruit diameter (mm), fruit firmness (g mm −1 of compression force), soluble solids content ( • Brix), total titratable acids (TTA), pH, plant height (cm), and plant vigor (rated using a 1-5 scale). The last two traits were only evaluated in the 2018 nursery.
For the evaluation of fruit quality and market-related traits (i.e., weight, diameter, firmness, • Brix, TTA, and pH), 15 fully mature berries from 20 individual plants were sampled from each family. Each plant was flagged and tagged with a sample ID number. Fruit weight (g) was determined using an analytical balance (CP2245, Sartorius Corp., Bohemia, NY, USA). The same 15 berries were then evaluated for firmness (g mm −1 of compression force) and fruit size diameter (mm) using FirmTech II equipment (BioWorks Inc., Wamego, KS, USA). Sensory quality traits were evaluated using the juice obtained from the same 15 berries, thus soluble solids ( • Brix) measures were obtained from pipetting 1 mL of juice onto a digital pocket refractometer (Atago, Inc., Bellevue, WA, USA), and total titratable acids were determined using an automatic titrator (Mettler-Toledo, DL 115, Inc., Columbus, OH, USA). The pH was measured using a glass electrode in the remaining juice (Mettler-Toledo, DL 115, Inc., Columbus, OH, USA). Plant height was measured in the field with a ruler as the distance between the base of the plant to top of the highest branch. Vigor was evaluated using a 1-5 rating scale, where 1 represented low vegetative growth and 5 represented high vegetative growth.

Models
One-step single-trait Bayesian linear mixed models were used to obtain the breeding values for each evaluated individual, as follows where y is a vector of the phenotypic values of the trait being analyzed, µ is the population's overall mean, n is the fixed effect of nursery, b is the random effect of the ith breeding stage, nxb is the random effect of the nursery by breeding stage interaction ∼ N 0, Iσ 2 nxb , g is the random effect of genotype ∼ N 0, Aσ 2 a , where A is the population's pedigree-based relationship matrix, and e is the random residual effect ∼ N 0, Iσ 2 e . Genotype effects were considered nested within nursery, since each nursery was constituted by different breeding populations containing individuals for each breeding stage. For traits measured in a single nursery, the same Equation (1) was used without considering the nursery terms. The variance components for each random variable were the additive σ 2 a , breeding stage σ 2 b , nursery-by-breeding stage interaction σ 2 nxb , and residual σ 2 e . X, Z 1 , Z 2 , and Z 3 were incidence matrices for nursery, breeding stage, nursery-by-breeding stage interaction, and genotype, respectively. A pseudo broad-sense heritability for each trait was obtained considering the ratio of the sum of the additive and breeding stage variances and the total phenotypic variance, and narrow-sense heritability was estimated considering the ratio between the additive variance component and the total phenotypic variance.

Pedigree information
A historical and extensive pedigree data set that contained information from 11,866 entries was used to obtain the relationship matrix for the individuals that composed this study. Pedigree information was obtained by combining data from NCGR-Corvallis Vaccinium Catalog [21], the Brooks and Olmo Register of Fruit and Nut Varieties [22,23], and internal pedigree records from the University of Florida blueberry breeding program, including the pedigree of the individuals used in this study. The R package AGHmatrix v. 1.0.2 [24] was used to obtain the pedigree-based relationship matrix (A) considering the autotetraploid model without double reduction as in Kerr et al. [25].

Model implementation
All models were fitted using the package BGLR v. 1.0.5. [26]. Predictions were based on 80,000 Gibbs sampler iterations, considering 40,000 for burn-in, and a thinning of five. Parameters' convergency was evaluated to define the final values used in the analysis. A single-step regression approach was applied to perform all pedigree-BLUP (P-BLUP) analysis. Default hyper-parameters were used, as previously described by Perez and de los Campos [27]. The packages lsmeans [28] and multcomp [29] were used to perform post hoc tests, using the Bonferroni multiple comparisons test to counteract the problem of multiple comparisons considering α = 0.05. Graphical visualizations were obtained using ggplot2 [30]. All analyses were implemented on the R platform [31].

Population Genetics
Large genetic variation was observed within the evaluated populations ( Figure 1). The first two principal components of the pedigree-based matrix (A) used in this study explained approximately 65% of the variance. Clear clusters were observed for F 2 families and for the independent SHB families used as benchmarks in the study. Given that the SHB families evaluated in this study were obtained from multiple biparental crosses, some of those crosses shared a higher relationship than others, therefore, one can see the grouping in the PCA based on the relationship shared among SHB as well as their relationship share with the other breeding stages. The individuals representing V. elliottii, F 1 , BC 1 , and the remaining SHB (i.e., individuals used in the crosses) presented higher genetic similarity and clustered together ( Figure 1). Table 1 summarizes the posterior means of the genetic parameters. All traits displayed genetic variances significantly different than zero, showing that selection can be performed both within and between breeding stages. For firmness, TTA, and plant vigor, the additive variance (σ 2 a ) represented, respectively, 50%, 38%, and 43% of the total variance observed. Conversely, for diameter, weight, and yield, more than 50% of the total variance was explained by breeding stage. Broad sense heritability results (h 2 ) were high (>0.70) for all traits, apart from plant height that presented h 2 = 0.53. Considering all traits, narrow-sense heritability values varied between 0.14 and 0.50, for yield and fruit firmness, respectively. my 2021, 11, x FOR PEER REVIEW 5 of 13 BC1, and the remaining SHB (i.e., individuals used in the crosses) presented higher genetic similarity and clustered together (Figure 1). Table 1 summarizes the posterior means of the genetic parameters. All traits displayed genetic variances significantly different than zero, showing that selection can be performed both within and between breeding stages. For firmness, TTA, and plant vigor, the additive variance ( ) represented, respectively, 50%, 38%, and 43% of the total variance observed. Conversely, for diameter, weight, and yield, more than 50% of the total variance was explained by breeding stage. Broad sense heritability results (h 2 ) were high (> 0.70) for all traits, apart from plant height that presented h 2 = 0.53. Considering all traits, narrow-sense heritability values varied between 0.14 and 0.50, for yield and fruit firmness, respectively.

Hybrid Performance
To assess the performance of hybrid families, nine plant fitness-, fruit quality-, and marketable-related traits were evaluated and compared with the results obtained for commercial germplasm grown in the same conditions as the hybrids (i.e., SHB). Results were also compared with the performance of the wild relative, to confirm improvement towards commercial requirements. Commercial thresholds used by the University of Florida Blueberry breeding program were shown in the plots to facilitate understanding, when existing.
Hybrid performance was above commercial requirements for fruit diameter, fruit weight, and fruit firmness (Figure 2). Backcrossed families did not significantly differ from SHB germplasm regarding fruit diameter (Figure 2A,B) and fruit weight ( Figure 2C). No significant differences were observed for fruit firmness between any of the breeding stages compared (Figure 2C), showing that requirements for size, weight, and firmness can be rapidly achieved after hybridization.

Hybrid Performance
To assess the performance of hybrid families, nine plant fitness-, fruit quality-, and marketable-related traits were evaluated and compared with the results obtained for commercial germplasm grown in the same conditions as the hybrids (i.e., SHB). Results were also compared with the performance of the wild relative, to confirm improvement towards commercial requirements. Commercial thresholds used by the University of Florida Blueberry breeding program were shown in the plots to facilitate understanding, when existing.
Hybrid performance was above commercial requirements for fruit diameter, fruit weight, and fruit firmness (Figure 2). Backcrossed families did not significantly differ from SHB germplasm regarding fruit diameter (Figure 2A,B) and fruit weight ( Figure 2C). No significant differences were observed for fruit firmness between any of the breeding stages compared ( Figure 2C), showing that requirements for size, weight, and firmness can be rapidly achieved after hybridization. V. elliottii presented a significantly higher amount of sugar (i.e., higher Brix) than both hybrids and SHB germplasm. Nonetheless, there was no significant difference for Brix between most of the hybrids and SHB families ( Figure 3A), showing that commercial V. elliottii presented a significantly higher amount of sugar (i.e., higher Brix) than both hybrids and SHB germplasm. Nonetheless, there was no significant difference for Brix between most of the hybrids and SHB families ( Figure 3A), showing that commercial sugar requirements can be promptly obtained through hybridization. For V. elliottii, however, pH was significantly lower and TTA values were significantly higher when compared with the hybrids and SHB ( Figure 3B,C). Moreover, some hybrid families (F 2 and BC 1 ) performed as well as or slightly better than SHB germplasm when considering pH and TTA values ( Figure 3B,C).
Plant fitness was evaluated by measuring total yield, plant height, and plant vigor. With the exception of F 2 , hybrids presented height values as good as or better than the ones observed for SHB germplasm ( Figure S1). The same pattern was also observed for vigor ( Figure S1), showing that interspecific hybridization between SHB and V. elliottii can produce healthy and robust progeny. However, plant yield can still be improved upon, and promising results are expected, given BC 1 families presented yield measures significantly higher than the other hybrids for nursery 2017 and not significantly different from SHB families for the 2018 nursery ( Figure 3D). Nevertheless, only one round of backcrossing was performed. BC1) performed as well as or slightly better than SHB germplasm when considering pH and TTA values ( Figure 3B,C).
Plant fitness was evaluated by measuring total yield, plant height, and plant vigor. With the exception of F2, hybrids presented height values as good as or better than the ones observed for SHB germplasm ( Figure S1). The same pattern was also observed for vigor ( Figure S1), showing that interspecific hybridization between SHB and V. elliottii can produce healthy and robust progeny. However, plant yield can still be improved upon, and promising results are expected, given BC1 families presented yield measures significantly higher than the other hybrids for nursery 2017 and not significantly different from SHB families for the 2018 nursery ( Figure 3D). Nevertheless, only one round of backcrossing was performed.  Figure 4 summarizes the performance of hybrids, wild germplasm, and SHB germplasm across crop seasons. Our results show that commercial requirements for fruit quality and marketable traits can be quickly achieved, given the BC1 families performed as well as SHB families or slightly better for most of the traits evaluated. Hybrid individuals that presented satisfactory performance were visually evaluated for fruit color, wax coating, scar, and stem detachment. Individual genotypes that reached the quality standards of the blueberry breeding program were flagged and kept in the nurseries for future evaluations, factoring in the possibility of evaluating resistance to biotic and abiotic factors and usage in next stages of breeding.  Figure 4 summarizes the performance of hybrids, wild germplasm, and SHB germplasm across crop seasons. Our results show that commercial requirements for fruit quality and marketable traits can be quickly achieved, given the BC 1 families performed as well as SHB families or slightly better for most of the traits evaluated. Hybrid individuals that presented satisfactory performance were visually evaluated for fruit color, wax coating, scar, and stem detachment. Individual genotypes that reached the quality standards of the blueberry breeding program were flagged and kept in the nurseries for future evaluations, factoring in the possibility of evaluating resistance to biotic and abiotic factors and usage in next stages of breeding.

Interspecific Hybridization in Blueberry
Interspecific hybridization is undeniably beneficial for modern agricultural advancements [7]. The use of wild species as raw material for breeding, both through classical and biotechnological techniques, can provide the genetic variability required to respond to environmental and demographical changes [6]. This process has contributed to the creation of commercial crops, such as sugarcane [32], banana [33], cotton [34], and bread wheat [35], as well as improving resistance and tolerance to biotic and abiotic factors. Prominent examples include bacterial blight resistance in rice [36], verticillium wilt resistance in canola [37], resistance to Phytophthora infestans and Globodera pallida in potatoes [38], and drought tolerance in wheat [39]. Introgression of traits from wild relatives has also improved fruit quality and marketable traits, such as fruit flavor, fruit size, and fruit texture in strawberry, tomatoes, and peppers [40,41].
Improvements obtained from interspecific hybridization in blueberry are not limited to the incorporation of novel traits, but also relate to the geographic expansion of the crop [3,19]. In 1948, Ralph H. Sharpe at the University of Florida developed the ground-breaking Southern Highbush Blueberry (SHB) that allowed for this expansion of blueberries into warmer climates [2]. The SHB blueberry encompasses cultivated hybrids that exhibit low cold hour requirements (i.e., chilling hours and hours of cold below 7 °C, e.g., ~300 h). This achievement resulted from the interspecific hybridization between HB, LB, Rabbiteye, and wild native diploid species [48,49], resulting in the development of hybrids with high vigor, resistance to diseases and abiotic stressors (i.e., heat and humidity), early maturation, vertical architecture, better fruit quality (e.g., firmness, color, and flavor), and traits consistent with commercial requirements [8,50]. Nevertheless, the global market is

Interspecific Hybridization in Blueberry
Interspecific hybridization is undeniably beneficial for modern agricultural advancements [7]. The use of wild species as raw material for breeding, both through classical and biotechnological techniques, can provide the genetic variability required to respond to environmental and demographical changes [6]. This process has contributed to the creation of commercial crops, such as sugarcane [32], banana [33], cotton [34], and bread wheat [35], as well as improving resistance and tolerance to biotic and abiotic factors. Prominent examples include bacterial blight resistance in rice [36], verticillium wilt resistance in canola [37], resistance to Phytophthora infestans and Globodera pallida in potatoes [38], and drought tolerance in wheat [39]. Introgression of traits from wild relatives has also improved fruit quality and marketable traits, such as fruit flavor, fruit size, and fruit texture in strawberry, tomatoes, and peppers [40,41].
Improvements obtained from interspecific hybridization in blueberry are not limited to the incorporation of novel traits, but also relate to the geographic expansion of the crop [3,19]. In 1948, Ralph H. Sharpe at the University of Florida developed the groundbreaking Southern Highbush Blueberry (SHB) that allowed for this expansion of blueberries into warmer climates [2]. The SHB blueberry encompasses cultivated hybrids that exhibit low cold hour requirements (i.e., chilling hours and hours of cold below 7 • C, e.g.,~300 h). This achievement resulted from the interspecific hybridization between HB, LB, Rabbiteye, and wild native diploid species [48,49], resulting in the development of hybrids with high vigor, resistance to diseases and abiotic stressors (i.e., heat and humidity), early maturation, vertical architecture, better fruit quality (e.g., firmness, color, and flavor), and traits consistent with commercial requirements [8,50]. Nevertheless, the global market is constantly requiring improvements and innovation in fruit quality, yield, and resources to meet the escalating pressures that threaten production such as global warming and biotic and abiotic resistance/tolerance. Breeding provides several perspectives to address these demands. Most notably, introgression of germplasm from wild species has demonstrated the ability to provide the genetic variability for valuable trait improvements in blueberries.

The Use of Vaccinium elliottii as a Genetical Resource
Vaccinium elliottii is continually recognized as an important source of genetic variability to improve favorable traits for blueberry. Its ability to endure extreme environmental conditions, produce highly aromatic berries with a short bloom to ripe period, and maintain an upright vegetative architecture make this wild species a strong candidate for interspecific hybridization with SHB [12,44]. Successful released varieties such as 'Snowchaser', 'Kestrel', and 'Carteret' were produced by introgressing V. elliottii germplasm into SHB [12,19]. However, given the high variability shown within the species [12], more studies are required to verify V. elliottii's potential as a genetic resource for blueberry breeding. Here, using a large breeding population composed of multiple hybrid stages, we show that promising results are expected when introgressing V. elliottii germplasm into SHB. Specifically, our study corroborated results from previous analysis, showing that market thresholds can be rapidly obtained in the hybrids [12], in addition to generating genetic variability using a wild species as genetic resource for breeding. Different from previous studies, such as Norden et al. [12] and Lyrene [20], we implemented a pedigree-based analysis using cross records and historic pedigrees for each hybrid, thus improving predictions for the families analyzed. To our knowledge, this is the first study characterizing such a large population of hybrids involving V. elliottii, as well as the first time that extensive pedigree data has been applied to generate breeding values to pre-breeding hybrid populations for blueberry.
We show that high genetic variability can be obtained when using V. elliottii as breeding resource. This is an important development for blueberry breeding, given previous studies illustrated declining heterozygosity of cultivated blueberry resulting from F. Coville's selective breeding founding event [19]. Both global food security and production are intrinsically connected to the wise use, conservation, and identification of biodiversity and wild germplasm genetic resources. These genetic resources can be classified as the raw material containing the genetic diversity necessary to build crops' resilience to evolving environmental and anthropogenic threats [6,7]. A broad-based foundation population, such as the one generated here, could be a new genetic resource for breeding [44]. This development broadens not only the genetic base of breeding populations but can also assist in developing successful commercial cultivars that can face escalating pressures threatening crop production. Given tetraploid intersectional hybrids between SHB and V. elliottii can be easily intercrossed with SHB genotypes [12], the segregating populations generated here can be further evaluated and used in different stages of breeding to capture positive and valuable traits.
One of the primary goals of breeding programs is the identification of superior genotypes when considering fruit quality. Current breeding efforts for blueberry focus on improving fruit firmness, fruit size, and flavor (i.e., sugar content, acidity, and volatiles) [2]. Fruit firmness is a key trait for blueberry breeding as it influences production/marketability from harvesting to transportability, shelf life, and consumer preference. Hand harvesting is normally associated with specialty crops, such as blueberry, which are destined for fresh market given the fragility of the fruit. Nevertheless, high costs and shortage of labor are directing the market to favor machine harvestable varieties [51]. Fruit firmness is a primary factor that allows for machine harvesting of blueberries, as firmer fruits better resist physical impacts [2,13,52]. Regarding postharvest behavior, firm fruits can withstand harvest handling and subsequent transportation better than soft fruits [53]. Even small differences in fruit firmness can affect shelf life [54]. High berry firmness is also associated with consumer liking. Gilbert et al. [55] has shown that firm fruits and crispness were generally preferred by consumers. Our results show that no significant differences for fruit firmness were observed between hybrids and SHB genotypes. Given that V. elliottii presents high firmness, commercial thresholds were quickly obtained through hybridization, with F 1 individuals even displaying high values for firmness ( Figure 2). Our outcomes corroborate results obtained in previous studies, showing that firmness values over 200 g mm −1 of compression force can be obtained in the first stages of hybridization [12].
Specific disadvantages can be observed when crossing SHB and V. elliottii. Among the major challenges is increasing fruit size and achieving commercial yield requirements [18]. For consumers' first purchases, appearance dictates preference. Gilbert et al. [55] illustrated that blueberry fruit size can additionally affect consumer preferences. The same is observed for other crop species, such as kiwi [56]. Lyrene and Sherman [18] and Lyrene [8] reported that berry size of BC 1 genotypes have the tendency to average below the normal commercial standards. Given the small berry size for V. elliottii, fruit size may be a difficult trait to recover when using the species for interspecific hybridization. Here, we show promising results in recovering fruit size when interbreeding V. elliottii and SHB. In our study, with only one round of backcrossing, genotypes presented diameter similar to the ones observed for SHB families across crop seasons ( Figure 4). Our results probably differed from previous studies due to the selected V. elliottii stock. Given V. elliottii's high diversity, one can expect resulting hybrids to display high phenotypic and genetic variability depending on the stock chosen. Similar results were observed, as expected, for fruit weight ( Figures 2C and 4), since fruit size and fruit weight are positively and significantly correlated [5,57]. The effect of fruit weight and size as components of total yield for blueberry [58] can also be seen in the promising yield values observed for BC 1 and F 2 stages across crop season (Figure 4), considering that only one round of backcrossing was performed.
Beyond fruit weight, size, and yield, flavor is a fundamental component of fruit quality that affects consumer preference and marketability. Consumers are willing to pay more for better tasting varieties, justifying breeding for traits involved with flavor perception [59]. Flavor is a complex trait that is expensive to routinely evaluate in breeding programs [60]. However, sugar content and acidity greatly contribute to overall consumer liking and are easier and less expensive to phenotype than volatile profiles. Ferrão et al. [60] verified that overall liking of blueberry is positively correlated with high sugar content and negatively correlated with high acidity. Here, we show that we can achieve commercial requirements for total soluble solids and TTA in the first stage of hybridization (i.e., F 1 ). Our results show that some BC 1 and F 2 families presented sugar content and acidity equivalent to or better than the values presented by SHB individuals. Additionally, V. elliottii presented a significantly higher quantity of total sugars than SHB and all hybrid stages evaluated. However, V. elliottii also presented a significantly higher acidity and a lower pH than SHB and the hybrids. Similar results were also obtained by Norden et al. [12]. These results show that even though V. elliottii fruits presented more sugars, the perception of sweetness can be altered by the higher quantity of acids. V. elliottii germplasm could be further examined to screen for genotypes low in acidity (i.e., higher pH and lower TTA) with greater positive volatile and sugar contents, generating information for better stocks. Meeting these requirements could facilitate generation of hybrids with a better flavor perception and could further improve values for these traits.

Conclusions
By analyzing a large population comprising five different breeding stages across two seasons, here, we illustrate that the introgression of wild germplasm for blueberry breeding can generate a large amount of diversity and can generate promising and quick results towards fruit quality improvement. The use of V. elliottii through interspecific hybridization resulted in hybrids that quickly attained commercial benchmarks for most of the traits analyzed (e.g., fruit size, fruit weight, sugars, and acidity), and displayed height and vigor comparable to SHB germplasm. The segregating populations generated here can be used in different stages of breeding to introgress positive and valuable traits into commercial SHB germplasm. These findings also present encouraging implications for the blueberry breeding industry, as future research can be performed to evaluate hybrid environmental adaptations and disease resistance, with the potential to corroborate the success of introgression of these traits from V. elliottii into SHB material.