Bud sport mutants of apple (Malus domestica Borkh.) trees with a highly blushed colouring pattern are mainly caused by the accumulation of anthocyanins in the fruit skin. Hormones are important factors modulating anthocyanin accumulation. However, a good understanding of the interplay between hormones and anthocyanin synthesis in apples, especially in mutants at the molecular level, remains elusive. Here, physiological and comparative transcriptome approaches were used to reveal the molecular basis of color pigmentation in the skin of ‘Red Delicious’ (G0) and its mutants, including ‘Starking Red’ (G1), ‘Starkrimson’ (G2), ‘Campbell Redchief’ (G3) and ‘Vallee spur’ (G4). Pigmentation in the skin gradually proliferated from G0 to G4. The anthocyanin content was higher in the mutants than in ‘Red Delicious’. The activation of early phenylpropanoid biosynthesis genes, including ASP3, PAL, 4CL, PER, CHS, CYP98A and F3’H, was more responsible for anthocyanin accumulation in mutants at the color break stage. In addition, IAA and ABA had a positive regulatory effect on the synthesis of anthocyanins, while GA had the reverse effect. The down-regulation of AACT1, HMGS, HMGR, MVK, MVD2, IDI1 and FPPS2 involved in terpenoid biosynthesis influences anthocyanin accumulation by positively regulating transcripts of AUX1 and SAUR that contribute to the synthesis of IAA, GID2 to GA, PP2C and SnRK2 to ABA. Furthermore, MYB and bHLH members, which are highly correlated (r=0.882–0.980) with anthocyanin content, modulated anthocyanin accumulation by regulating the transcription of structural genes, including CHS and F3’H, involved in the flavonoid biosynthesis pathway. The present comprehensive transcriptome analyses contribute to the understanding of the the relationship between hormones and anthocyanin synthesis as well as the molecular mechanism involved in apple skin pigmentation.

Heritable DNA methylation is a highly conserved epigenetic mark that is important for many biological processes. In a previous transcriptomic study on the fruit skin pigmentation of apple (Malus domestica Borkh.) cv. ‘Red Delicious’ (G0) and its four continuous-generation bud sport mutants including ‘Starking Red’ (G1), ‘Starkrimson’ (G2), ‘Campbell Redchief’ (G3) and ‘Vallee spur’ (G4), we identified MYB transcription factors (TFs) MdLUX and MdPCL-like involved in regulating anthocyanin synthesis. However, how these TFs ultimately determine the fruit skin colour traits remain elusive. Here, bioinformatics analysis revealed that MdLUX and MdPCL-like contained a well-conserved motif SH[AL]QKY[RF] in their C-terminal region and were located in the nucleus of onion epidermal cells. Overexpression of MdLUX and MdPCL-like in ‘Golden Delicious’ fruits, ‘Gala’ calli and Arabidopsis thaliana promoted the accumulation of anthocyanin, whereas MdLUX and MdPCL-like suppression inhibited anthocyanin accumulation in ‘Red Fuji’ apple fruit skin. Yeast one-hybrid assays revealed that MdLUX and MdPCL-like may bind to the promoter region of the anthocyanin biosynthesis gene MdF3H. Dual-luciferase assays indicated that MdLUX and MdPCL-like activated MdF3H. The whole-genome DNA methylation study revealed that the methylation levels of the ^mCG context at the upstream (i.e., promoter region) of MdLUX and MdPCL-like were inversely correlated with their mRNA levels and anthocyanin accumulation. Hence, the data suggest that MYB_SH[AL]QKY[RF] TFs MdLUX and MdPCL-like promote anthocyanin biosynthesis in apple fruit skins through the DNA hypomethylation of their promoter regions and the activation of the structural flavonoid gene MdF3H.

Q-type C2H2 zinc finger proteins (ZFPs), the largest family of transcription factors, have been extensively studied in plant genomes. However, the genes encoding this transcription factor family have not been explored in grapevine genomes. Therefore, in this study, we conducted a genome-wide identification of ZFP genes in three species of grapevine, namely Vitis vinifera, Vitis riparia, and Vitis amurensis, based on the sequence databases and phylogenetic and their conserved domains. We identified 52, 54, and 55 members of Q-type C2H2 ZFPs in V. vinifera, V. riparia, and V. amurensis, respectively. The physical and chemical properties of VvZFPs, VrZFPs, and VaZFPs were examined. The results showed that these proteins exhibited differences in the physical and chemical properties and that they all were hydrophobic proteins; the instability index showed that the four proteins were stable. The subcellular location of the ZFPs in the grapevine was predicted mainly in the nucleus. The phylogenetic tree analysis of the amino acid sequences of VvZFP, VaZFP, VrZFP, and AtZFP proteins showed that they were closely related and were divided into six subgroups. Chromosome mapping analysis showed that VvZFPs, VrZFPs, and VaZFPs were unevenly distributed on different chromosomes. The clustered gene analysis showed that the motif distribution was similar and the sequence of genes was highly conserved. Exon and intron structure analysis showed that 118 genes of ZFPs were intron deletion types, and the remaining genes had variable numbers of introns, ranging from 2 to 15. Cis-element analysis showed that the promoter of VvZFPs contained multiple cis-elements related to plant hormone response, stress resistance, and growth, among which the stress resistance elements were the predominant elements. Finally, the expression of VvZFP genes was determined using real-time quantitative PCR, which confirmed that the identified genes were involved in response to methyl jasmonate (MeJA), abscisic acid (ABA), salicylic acid (SA), and low-temperature (4 °C) stress. VvZFP10-GFP and VvZFP46-GFP fusion proteins were localized in the nucleus of tobacco cells, and VvZFP10 is the most responsive gene among all VvZFPs with the highest relative expression level to MeJA, ABA, SA and low-temperature (4 °C) stress. The present study provides a theoretical basis for exploring the mechanism of response to exogenous hormones and low-temperature tolerance in grapes and its molecular breeding in the future.

Caffeic acid-O-methyltransferase (COMT) is a key enzyme in lignin synthesis and secondary metabolism in plants, and it participates in the regulation of plant growth and development as well as plants’ stress response. To further investigate the function of COMT in grapevine, a total of 124 COMT family genes were identified from three Vitis species in this study, namely Pinot noir (Vitis vinifera L.), Vitis amurensis, and Vitis riparia. The amino acid sequence encoded by these genes ranged from 55 to 1422 aa, and their molecular mass ranged from 6640.82 to 77,034.43 Da. Subcellular localization prediction inferred that they were mainly located in the plasma membrane and cytoplasm. The prediction of secondary structures showed that α-helix and irregular coiled-coil were primary structural elements. These genes were unevenly distributed across 10 different chromosomes, respectively. Phylogenetic tree analysis of the amino acid sequences of VvCOMT, VaCOMT, VrCOMT, and AtCOMT proteins showed that they were closely related and were divided into four subgroups. The motif distribution was similar among the cluster genes, and the gene sequence was notably conserved. The 124 members of the COMT gene family possessed a variable number of exons, ranging from 2 to 13. The promoter region of all of these COMTs genes contained multiple cis-acting elements related to hormones (e.g., ABA, IAA, MeJA, GA, and SA), growth and development (e.g., endosperm, circadian, meristem, light response), and various stress responses (e.g., drought, low temperature, wounding, anaerobic, defense, and stress). The intraspecies collinearity analysis suggested that there were one pair, three pairs, and six pairs of collinear genes in Va, Pinot noir, and Vr, respectively, and that tandem duplication contributed more to the expansion of these gene family members. In addition, interspecific collinearity revealed that the VvCOMTs had the strongest homology with the VaCOMTs, followed by the VrCOMTs, and the weakest homology with the AtCOMTs. The expression patterns of different tissues and organs at different developmental stages indicated that the VvCOMT genes had obvious tissue expression specificity. The majority of VvCOMT genes were only expressed at higher levels in certain tissues. Furthermore, we screened 13 VvCOMT genes to conduct qRT-PCR verification according to the transcriptome data of VvCOMTs under abiotic stresses (NaCl, PEG, and cold). The results confirmed that these genes were involved in the responses to NaCl, PEG, and cold stress. This study lays a foundation for the exploration of the function of the COMT genes, and is of great importance for the genetic improvement of abiotic stress resistance in grapes.