Comparative evaluation of reference genomes and cell-type annotation frameworks for single-nucleus transcriptomic analysis in apple

Single-nucleus RNA sequencing (snRNA-seq) provides a powerful tool to profile cell-type-specific transcriptional programs in woody fruit trees. Yet its application in Malus domestica has been limited by heterogeneous reference genomes and the lack of validated marker genes. Here, we first established a comprehensive set of marker genes corresponding to 29 plant cell types and identified their orthologs in apple. Using these resources, we systematically evaluated six publicly available apple genome assemblies as references for snRNA-seq analysis. Our results demonstrate that reference genome selection has a major impact on mapping performance and downstream cell type annotation, with the GDDH13 assembly consistently producing the highest transcript counts per nucleus and the greatest number of identifiable transcriptional clusters (15 clusters representing 11 distinct cell types across 35,557 high-quality nuclei). To achieve robust cell type annotation, we integrated classical marker gene-based annotation with two cross-species computational frameworks — Orthologous Marker Gene Groups (OMGs) and XSpeciesSpanner. These three approaches enabled confident identification of 11 major cell types in apple seedlings and provided cross-validation for ambiguous or previously unannotated clusters. Pseudotime analysis reconstructed the developmental trajectories of both procambial and stomatal lineages, revealing their progressive transcriptional transitions and enabling the delineation of key molecular programs underlying vascular differentiation and early guard-cell formation. Together, this work provides a reference-quality single-nucleus transcriptomic atlas for apple and underlined the critical impact of genome choice and multi-method cell-type annotation on single-cell studies in non-model crops.