Apple (Malus × domestica) trees are vulnerable to freezing temperatures. However, there has been only limited success in developing cold-hardy cultivars. This lack of progress is due at least partly to lack of understanding of the molecular mechanisms of freezing tolerance in apple. In this study, we evaluated the potential roles for two R2R3 MYB transcription factors (TFs), MYB88 and the paralogous FLP (MYB124), in cold stress in apple and Arabidopsis. We found that MYB88 and MYB124 positively regulate freezing tolerance and cold-responsive gene expression in both apple and Arabidopsis. Chromatin-Immunoprecipitation-qPCR and electrophoretic mobility shift assays showed that MdMYB88/MdMYB124 act as direct regulators of the COLD SHOCK DOMAIN PROTEIN 3 (MdCSP3) and CIRCADIAN CLOCK ASSOCIATED 1 (MdCCA1) genes. Dual luciferase reporter assay indicated that MdCCA1 but not MdCSP3 activated the expression of MdCBF3 under cold stress. Moreover, MdMYB88 and MdMYB124 promoted anthocyanin accumulation and H O detoxification in response to cold. Taken together, our results suggest that MdMYB88 and MdMYB124 positively regulate cold hardiness and cold-responsive gene expression under cold stress by C-REPEAT BINDING FACTOR (CBF)-dependent and CBF-independent pathways.

Water deficit is one of the main challenges for apple (Malus × domestica) growth and productivity. Breeding drought tolerant cultivars depends on a thorough understanding of the drought responses of apple trees. Here, we identified the zinc-finger protein B-BOX 7/CONSTANS-LIKE 9 (MdBBX7/MdCOL9), which plays a positive role in apple drought tolerance. Overexpression of MdBBX7 enhanced drought tolerance, whereas knocking down MdBBX7 expression reduced it. Chromatin immunoprecipitation (ChIP)-seq analysis identified one cis-element of MdBBX7, CCTTG, as well as its known binding motif, the T/G box. ChIP-seq and RNA-seq identified 1197 direct targets of MdBBX7, including ETHYLENE RESPONSE FACTOR (ERF1), EARLY RESPONSIVE TO DEHYDRATION 15 (ERD15), and GOLDEN2-LIKE 1 (GLK1) and these were further verified by ChIP-qPCR and electronic mobility shift assays. Yeast two-hybrid screen identified an interacting protein of MdBBX7, RING-type E3 ligase MYB30-INTERACTING E3 LIGASE 1 (MIEL1). Further examination revealed that MdMIEL1 could mediate the ubiquitination and degradation of MdBBX7 by the 26S proteasome pathway. Genetic interaction analysis suggested that MdMIEL1 acts as an upstream factor of MdBBX7. In addition, MdMIEL1 was a negative regulator of the apple drought stress response. Taken together, our results illustrate the molecular mechanisms by which the MdMIEL1-MdBBX7 module influences the response of apple to drought stress.

DNA-binding one zinc-finger (Dof) proteins constitute a family of transcription factors with a highly conserved Dof domain that contains a C2C2 zinc-finger motif. Although several studies have demonstrated that Dof proteins are involved in multiple plant processes, including development and stress resistance, the functions of these proteins in drought stress resistance are largely unknown. Here, we report the identification of the MdDof54 gene from apple and document its positive roles in apple drought resistance. After long-term drought stress, compared with nontransgenic plants, MdDof54 RNAi plants had significantly shorter heights and weaker root systems; the transgenic plants also had lower shoot and root hydraulic conductivity, as well as lower photosynthesis rates. By contrast, compared with nontransgenic plants, MdDof54-overexpressing plants had higher photosynthesis rates and shoot hydraulic conductivity under long-term drought stress. Moreover, compared with nontransgenic plants, MdDof54-overexpressing plants had higher survival percentages under short-term drought stress, whereas MdDof54 RNAi plants had lower survival percentages. MdDof54 RNAi plants showed significant downregulation of 99 genes and significant upregulation of 992 genes in response to drought, and 366 of these genes were responsive to drought. We used DAP-seq and ChIP-seq analyses to demonstrate that MdDof54 recognizes cis-elements that contain an AAAG motif. Taken together, our results provide new information on the functions of MdDof54 in plant drought stress resistance as well as resources for apple breeding aimed at the improvement of drought resistance.

The evolutionary history of the Malus genus has not been well studied. In the current study, we presented genetic evidence on the origin of the Malus genus based on genome sequencing of 297 Malus accessions, revealing the genetic relationship between wild species and cultivated apples. Our results demonstrated that North American and East Asian wild species are closer to the outgroup (pear) than Central Asian species, and hybrid species including natural (separated before the Pleistocene, about 2.5 Mya) and artificial hybrids (including ornamental trees and rootstocks) are between East and Central Asian wild species. Introgressions from M. sylvestris in cultivated apples appeared to be more extensive than those from M. sieversii, whose genetic background flowed westward across Eurasia and eastward to wild species including M. prunifolia, M. × asiatica, M. × micromalus, and M. × robust. Our results suggested that the loss of ancestral gene flow from M. sieversii in cultivated apples accompanied the movement of European traders around the world since the Age of Discovery. Natural SNP variations showed that cultivated apples had higher nucleotide diversity than wild species and more unique SNPs than other apple groups. An apple ERECTA‐like gene that underwent selection during domestication on 15th chromosome was identified as a likely major determinant of fruit length and diameter, and an NB‐ARC domain‐containing gene was found to strongly affect anthocyanin accumulation using a genome‐wide association approach. Our results provide new insights into the origin and domestication of apples and will be useful in new breeding programmes and efforts to increase fruit crop productivity.

SUMMARY Apple (Malus domestica) trees are vulnerable to freezing temperatures. Cold resistance in woody perennial plants can be improved through biotechnological approaches. However, genetic engineering requires a thorough understanding of the molecular mechanisms of the tree's response to cold. In this study, we demonstrated that the Mdm‐miR160–MdARF17–MdWRKY33 module is crucial for apple freezing tolerance. Mdm‐miR160 plays a negative role in apple freezing tolerance, whereas MdARF17, one of the targets of Mdm‐miR160, is a positive regulator of apple freezing tolerance. RNA sequencing analysis revealed that in apple, MdARF17 mediates the cold response by influencing the expression of cold‐responsive genes. EMSA and ChIP‐qPCR assays demonstrated that MdARF17 can bind to the promoter of MdWRKY33 and promotes its expression. Overexpression of MdWRKY33 enhanced the cold tolerance of the apple calli. In addition, we found that the Mdm‐miR160–MdARF17–MdWRKY33 module regulates cold tolerance in apple by regulating reactive oxygen species (ROS) scavenging, as revealed by (i) increased H2O2 levels and decreased peroxidase (POD) and catalase (CAT) activities in Mdm‐miR160e OE plants and MdARF17 RNAi plants and (ii) decreased H2O2 levels and increased POD and CAT activities in MdmARF17 OE plants and MdWRKY33 OE calli. Taken together, our study uncovered the molecular roles of the Mdm‐miR160–MdARF17–MdWRKY33 module in freezing tolerance in apple, thus providing support for breeding of cold‐tolerant apple cultivars.