Background: Blueberry (Vaccinium spp.) is characterized by the production of berries that are smaller than most common fruits, and the underlying mechanisms of fruit size in blueberry remain elusive. V. corymbosum ‘O’Neal’ and ‘Bluerain’ are commercial southern highbush blueberry cultivars with large- and small-size fruits, respectively, which mature ‘O’Neal’ fruits are 1 ~ 2-fold heavier than those of ‘Bluerain’. In this study, the ontogenetical patterns of ‘O’Neal’ and ‘Bluerain’ hypanthia and fruits were compared, and comparative transcriptomic analysis was performed during early fruit development. Results: V. corymbosum ‘O’Neal’ and ‘Bluerain’ hypanthia and fruits exhibited intricate temporal and spatial cell proliferation and expansion patterns. Cell division before anthesis and cell expansion after fertilization were the major restricting factors, and outer mesocarp was the key tissue affecting fruit size variation among blueberry genotypes. Comparative transcriptomic and annotation analysis of differentially expressed genes revealed that the plant hormone signal transduction pathway was enriched, and that jasmonate-related TIFYs genes might be the key components orchestrating other phytohormones and influencing fruit size during early blueberry fruit development. Conclusions: These results provided detailed ontogenetic evidence for determining blueberry fruit size, and revealed the important roles of phytohormone signal transductions involving in early fruit development. The TIFY genes could be useful as markers for large-size fruit selection in the current breeding programs of blueberry.

Cross-pollination can improve the fruit set and quality of blueberry (Vaccinium spp.) for growers and consumers. However, the xenia effect in southern highbush blueberry remains unclear. Therefore, we selected eight cultivars of southern highbush blueberry (Vaccinium corymbosum L., interspecific hybrids) and applied pollination treatments (i.e., artificial self-pollination, artificial pollination with mixed pollen, or artificial pollination with individual cultivar pollen) to explore the xenia effects on the fruit set and quality of ‘O’Neal’ and ‘Emerald’. Pollen viability tests indicated that all of the cultivars are capable of fertilization. The highest fruit set came from ‘Bluerain’ pollen for ‘O’Neal’, while ‘Gulfcoast’ pollen increased fruit set the most in ‘Emerald’. Principal component analysis revealed that the cross combinations ‘Emerald’ × ‘Gulfcoast’ and ‘O’Neal’ × ‘Gulfcoast’ optimized the external and interior quality of the fruit. SSR was applied to determine which pollen source yielded the most seedlings. Results indicated that ‘Emerald’ × ‘Gulfcoast’ and ‘O’Neal’ × ‘Bluerain’ increased seedling production. Our results demonstrate that the xenia effects of ‘Gulfcoast’ pollen may increase ‘Emerald’ yields and promote fruit quality, while pollen from ‘Bluerain’ or ‘Jewel’ can improve ‘O’Neal’ fruit quality and seed number. Hence, these cross combinations may be utilized in blueberry production to increase fruit set, yield, and quality.

MADS-box is a class of transcriptional regulators that are ubiquitous in plants and plays important roles in the process of plant growth and development. Identification and analysis of blueberry MADS-box genes can lay a foundation for their function investigations. In the present study, 249 putative MADS-box genes were identified in the blueberry genome. Those MADS-box genes were distributed on 47 out of 48 chromosomes. The phylogenetic and evolutionary analyses showed that blueberry MADS-box genes were divided into 131 type I members and 118 type II members. The type I genes contained an average of 1.89 exons and the type II genes contained an average of 7.83 exons. Motif analysis identified 15 conserved motifs, of which 4 were related to the MADS domain and 3 were related to the K-box domain. A variety of cis-acting elements were found in the promoter region of the blueberry MADS-box gene, indicating that the MADS-box gene responded to various hormones and environmental alterations. A total of 243 collinear gene pairs were identified, most of which had a Ka/Ks value of less than 1. Nine genes belonging to SEP, AP3/PI, and AGL6 subfamilies were screened based on transcriptomic data. The expression patterns of those nine genes were also verified using quantitative PCR, suggesting that VcMADS6, VcMADS35, VcMADS44, VcMADS58, VcMADS125, VcMADS188, and VcMADS212 had potential functions in blueberry fruit ripening. The results of this study provide references for an in-depth understanding of the biological function of the blueberry MADS-box genes and the mechanism of blueberry fruit ripening.