Miniature inverted-repeat transposable elements (MITEs) are widely distributed in the plant genome and can be methylated. However, whether DNA methylation of MITEs is associated with induced allelic expression and drought tolerance is unclear. Here, we identified the drought-inducible MdRFNR1 (root-type ferredoxin-NADP⁺ oxidoreductase) gene in apple (Malus domestica). MdRFNR1 plays a positive role in drought tolerance by regulating the redox system, including increasing NADP⁺ accumulation and catalase and peroxidase activities and decreasing NADPH levels. Sequence analysis identified a MITE insertion (MITE-MdRF1) in the promoter of MdRFNR1-1 but not the MdRFNR1-2 allele. MdRFNR1-1 but not MdRFNR1-2 expression was significantly induced by drought stress, which was positively associated with the MITE-MdRF1 insertion and its DNA methylation. The methylated MITE-MdRF1 is recognized by the transcriptional anti-silencing factors MdSUVH1 and MdSUVH3, which recruit the DNAJ domain-containing proteins MdDNAJ1, MdDNAJ2, and MdDNAJ5, thereby activating MdRFNR1-1 expression under drought stress. Finally, we showed that MdSUVH1 and MdDNAJ1 are positive regulators of drought tolerance. These findings illustrate the molecular roles of methylated MITE-MdRF1 (which is recognized by the MdSUVH–MdDNAJ complex) in induced MdRFNR1-1 expression as well as the drought response of apple and shed light on the molecular mechanisms of natural variation in perennial trees.

The apple is a favorite fruit for human diet and is one of the most important commercial fruit crops around the world. Investigating metabolic variations during fruit development can provide a better understanding on the formation of fruit quality. The present study applied a widely targeted LC-MS-based metabolomics approach with large-scale detection, identification and quantification to investigate the widespread metabolic changes during “Pinova” apple development and ripening. A total of 462 primary and secondary metabolites were simultaneously detected, and their changes along with the four fruit-development stages were further investigated. The results indicated that most of the sugars presented increasing accumulation levels while organic acid, including Tricarboxylic acid cycle (TCA) intermediates, showed a distinct decreasing trend across the four fruit-development stages. A total of 207 secondary metabolites consisted of 104 flavonoids and 103 other secondary metabolites. Many flavonoids maintained relatively high levels in the early fruit stage and then rapidly decreased their levels at the following developmental stages. Further correlation analyses of each metabolite–metabolite pair highlighted the cross talk between the primary and secondary metabolisms across fruit development and ripening, indicating the significant negative correlations between sugars and secondary metabolites. Moreover, transcriptome analysis provided the molecular basis for metabolic variations during fruit development. The results showed that most differentially expressed genes (DEGs) involved in the TCA cycle were upregulated from the early fruit stage to the preripening stage. The extensive downregulation of controlling genes involved in the flavonoid pathway is probably responsible for the rapid decrease of flavonoid content at the early fruit stage. These data provide a global view of the apple metabolome and a comprehensive analysis on metabolomic variations during fruit development, providing a broader and better understanding on the molecular and metabolic basis of important fruit quality traits in commercial apples.

Drought stress is an important environmental factor limiting apple yield and fruit quality. Previously, we identified GRETCHEN HAGEN3.6 (GH3.6) as a negative regulator of drought stress in apple trees. Using transgenic MdGH3 RNAi (knocking down MdGH3.6 and its five homologs) plants as rootstock can increase drought tolerance, water use efficiency, flowering, and fruit quality of the Fuji scion after drought stress. However, the molecular mechanism behind this phenomenon is still unknown. Here, we performed transcriptome sequencing of the grafted plants (Fuji/GL-3 where Fuji was used as the scion and non-transgenic GL-3 was used as the rootstock, and Fuji/MdGH3 RNAi where MdGH3 RNAi was used as the rootstock) under control and drought conditions. Under control conditions, 667 up-regulated genes and 176 down-regulated genes were identified in the scion of Fuji/MdGH3 RNAi, as compared to the scion of Fuji/GL-3. Moreover, 941 up-regulated genes and 2226 down-regulated genes were identified in the rootstock of MdGH3 RNAi plants relative to GL-3. GO terms of these differentially expressed genes (DEGs) in scion and rootstock showed associations with plant growth, fruit development, and stress responses. After drought stress, 220 up-regulated and 452 down-regulated genes were identified in MdGH3 RNAi rootstock, as compared to GL-3. Significantly enriched GO terms included response to abiotic stimulus, cell division, microtubule-based process, metabolic and biosynthetic process of flavonoid, pigment, and lignin. The comparison between the scion of Fuji/MdGH3 RNAi and Fuji/GL-3 yielded a smaller number of DEGs; however, all of them were significantly enriched in stress-related GO terms. Furthermore, 365 and 300 mRNAs could potentially move from MdGH3 RNAi rootstock to scion under control and drought conditions, respectively, including FIDDLEHEAD (FDH), RESPONSIVE TO DESICCATION 26 (RD26), ARS-binding factor 2 (ABF2), WRKY75, and ferritin (FER). Overall, our work demonstrates the effects of rootstock on scion at the transcriptional level after drought stress and provides theoretical support for further understanding and utilization of MdGH3 RNAi plants.