Natural variations in MdNAC18 exert major genetic effect on apple fruit harvest date by regulating ethylene biosynthesis genes

Dissecting the genetic control of apple fruit harvest date (AFHD) into multiple Mendelian factors poses a significant challenge in modern genetics. Here, a quantitative trait locus (QTL) for AFHD was fine-mapped to the NAC transcription factor (TF) MdNAC18 within the interval defined by the overlap of QTLs Z03.5/Z03.6 and F03.2/F03.3. One direct target of MdNAC18 is the ethylene biosynthesis gene MdACO1. The single nucleotide polymorphisms (SNPs) SNP517 and SNP958 in the MdNAC18 coding sequence modulated activation of MdACO1 by MdNAC18. SNP1229 in the MdACO1 promoter destroyed the MdNAC18 binding site and thus abolished MdNAC18 binding. SNP517 and SNP958 also affected MdNAC18 activation of the TF gene MdARF5; MdARF5 activates the ethylene biosynthesis gene MdACS1. SNP517 and SNP958 in MdNAC18, SNP1229 and SNP769 (linked to InDel62) in MdACO1, and InDel162 in MdACS1 constituted a genetic variation network. The genetic effect of this network on AFHD was estimated as 60.3 d, accounting for 52.6% of the phenotype variation of the training population. The joint effects of these polymorphisms increased the accuracy of a genomics-assisted prediction (GAP) model for AFHD (r = 0.7125). Together, our results suggest that genetic variation in MdNAC18 affects AFHD by modulating ethylene biosynthesis and provide an optimized GAP model for apple breeding.