Among crop fruit trees, the apricot (Prunus armeniaca) provides an excellent model to study divergence and adaptation processes. Here, we obtain nearly 600 Armeniaca apricot genomes and four high-quality assemblies anchored on genetic maps. Chinese and European apricots form two differentiated gene pools with high genetic diversity, resulting from independent domestication events from distinct wild Central Asian populations, and with subsequent gene flow. A relatively low proportion of the genome is affected by selection. Different genomic regions show footprints of selection in European and Chinese cultivated apricots, despite convergent phenotypic traits, with predicted functions in both groups involved in the perennial life cycle, fruit quality and disease resistance. Selection footprints appear more abundant in European apricots, with a hotspot on chromosome 4, while admixture is more pervasive in Chinese cultivated apricots. Our study provides clues to the biology of selected traits and targets for fruit tree research and breeding.

Malus baccata is one of four wild apple species that can hybridize with the cultivated apple species (Malus domestica). It is widely used in high-latitude apple-producing areas as a rootstock and breeding resource because of its disease resistance, and cold tolerance. A lack of a reference genome has limited the application of M. baccata for apple breeding. We present a draft reference genome for M. baccata. The assembled sequence consisting of 665 Mb, with a scaffold N50 value of 452 kb, included transposable elements (413 Mb) and 46,114 high-quality protein-coding genes. According to a genetic map derived from 390 sibling lines, 72% of the assembly and 85% of the putative genes were anchored to 17 linkage groups. Many of the M. baccata genes under positive selection pressure were associated with plant-pathogen interaction pathways. We identified 2,345 Transcription factor-encoding genes in 58 families in the M. baccata genome. Genes related to disease defense and cold tolerance were also identified. A total of 462 putative nucleotide-binding site (NBS)-leucine-rich-repeat (LRR) genes, 177 Receptor-like kinase (RLK) and 51 receptor-like proteins (RLP) genes were identified in this genome assembly. The M. baccata genome contained 3978 cold-regulated genes, and 50% of these gene promoter containing DREB motif which can be induced by CBF gene. We herein present the first M. baccata genome assembly, which may be useful for exploring genetic variations in diverse apple germplasm, and for facilitating marker-assisted breeding of new apple cultivars exhibiting resistance to disease and cold stress.

Floral induction is an important stage in the apple tree life cycle. In ‘Nagafu No. 2’, which was derived from a ‘Fuji’ bud sport, flower bud formation is associated with serious problems, such as fewer and inferior flower buds, a long juvenile phase, and an alternate bearing phenotype. Moreover, the molecular regulatory mechanisms underlying apple floral induction remain unknown. To characterize these mechanisms, we compared the RNA-sequencing-based transcriptome profiles of buds during floral induction in profusely flowering ‘Qinguan’ and weakly flowering ‘Nagafu No. 2’ apple varieties. Genes differentially expressed between the buds of the two varieties were mainly related to carbohydrate, fatty acid, and lipid pathways. Additionally, the steady up-regulated expression of genes related to the fatty acid and lipid pathways and the down-regulated expression of starch synthesis-related genes in the carbon metabolic pathway of ‘Qinguan’ relative to ‘Nagafu No. 2’ were observed to contribute to the higher flowering rate of ‘Qinguan’. Additionally, global gene expression profiling revealed that genes related to cytokinin, indole-3-acetic acid, and gibberellin synthesis, signalling, and responses (i.e., factors contributing to cell division and differentiation and bud growth) were significantly differentially expressed between the two varieties. The up-regulated expression of genes involved in abscisic acid and salicylic acid biosynthesis via shikimate pathways as well as jasmonic acid production through fatty acid pathways in ‘Qinguan’ buds were also revealed to contribute to the floral induction and relatively high flowering rate of this variety. The differential expression of transcription factor genes (i.e., SPL, bZIP, IDD, and MYB genes) involved in multiple biological processes was also observed to play key roles in floral induction. Finally, important flowering genes (i.e., FT, FD, and AFL) were significantly more highly expressed in ‘Qinguan’ buds than in ‘Nagafu No. 2’ buds during floral induction. A complex genetic network of regulatory mechanisms involving carbohydrate, fatty acid, lipid, and hormone pathways may mediate the induction of apple tree flowering.

Somatic variations can give rise to bud sports with advantageous traits, serving as the foundation for bud sport breeding in perennial plants. Here, we report a fully phased genome assembly of ‘Fuji’ apple, enabling comprehensive identification of somatic variants across 74 clonally propagated ‘Fuji’ varieties. Phylogenetic analysis indicates that spur-type and early-maturation traits in ‘Fuji’ sport varieties arise from multiple independent events. Several putative functional somatic variants have been identified, including a spur-type-specific deletion in the promoter of the TCP transcription factor gene MdTCP11. DNA methylation level of the deletion-associated miniature inverted-repeat transposable element is lower in spur-type varieties compared to standard-type varieties, while the expression of MdTCP11 is significantly higher. Overexpression of MdTCP11 in apple decreases plant height, highlighting its important role in the development of spur-type apple varieties. This study sheds light on the cloning history of ‘Fuji’ and provides valuable resources for apple breeding.

Grafting can improve the agricultural traits of crop plants, especially fruit trees. However, the regulatory networks and differentially expressed microRNAs (miRNAs) related to grafting in apple remain unclear. Herein, we conducted high-throughput sequencing and identified differentially expressed miRNAs among self-rooted Fuji, self-rooted M9, and grafted Fuji/M9. We analyzed the flowering rate, leaf morphology, and nutrient and carbohydrate contents in the three materials. The flowering rate, element and carbohydrate contents, and expression levels of flowering genes were higher in Fuji/M9 than in Fuji. We detected 206 known miRNAs and 976 novel miRNAs in the three materials, and identified those that were up- or downregulated in response to grafting. miR156 was most abundant in Fuji, followed by Fuji/M9, and then self-rooted M9, while miR172 was most abundant in M9, followed by Fuji/M9, and then Fuji. These expression patterns suggest that that these miRNAs were related to grafting. A Gene Ontology (GO) analysis showed that the differentially expressed miRNAs controlled genes involved in various biological processes, including cellular biosynthesis and metabolism. The expression of differentially expressed miRNAs and flowering-related genes was verified by qRT-PCR. Altogether, this comprehensive analysis of miRNAs related to grafting provides valuable information for breeding and grafting of apple and other fruit trees.

Fruit acidity determines the organoleptic quality and nutritive value of most fruits. In litchi, although the organic acid composition of pulps is known, the molecular mechanisms and genes underlying variation in fruit acidity remain elusive. Herein, developing pulps of two contrasting litchi varieties, Huaizhi (HZ, low-acidity) and Boye_No.8 (B8, high-acidity), were subjected to metabolomics and transcriptomics, and the dynamic metabolome and transcriptional changes were determined. Measurements revealed that the dominant acidity-related organic acid in litchi pulps is malate, followed in low levels by citrate and tartrate. Variation in litchi pulps’ acidity is mainly associated with significant differences in malate and citrate metabolisms during fruit development. Malic acid content decreased by 91.43% and 72.28% during fruit ripening in HZ and B8, respectively. The content of citric acid increased significantly in B8, while in HZ it was reduced considerably. Differentially accumulated metabolites and differentially expressed genes analyses unveiled fumarate, succinate, 2-oxoglutarate, GABA (γ-aminobutyric acid), phosphoenolpyruvate, and citrate metabolisms as the key driving pathways of litchi fruits’ acidity variation. The drastic malate and citrate degradation in HZ was linked to higher induction of fumarate and GABA biosynthesis, respectively. Thirty candidate genes, including three key genes (LITCHI026501.m2, fumarase; LITCHI020148.m5, glutamate decarboxylase; and LITCHI003343.m3, glutamate dehydrogenase), were identified for functional studies toward genetic modulation of litchi fruit acidity. Our findings provide insights into the molecular basis of acidity variation in litchi and provide valuable resources for fruit quality improvement.