Fig. 5

Analysis of key genes related to anthocyanin accumulation. A Number of DEGs in different comparison groups. B Venn diagram shows DEGs in drought treated fruit vs. control fruit. C Module/trait correlations and corresponding p values. The right panel shows a colour scale for module/trait correlations from − 1 to 1. D Venn diagram shows DEGs in group1 and the blue and red modules. E Heatmap of the expression levels of DEGs in the intersection genes between group1 and the blue and red modules. F qRT-PCR detection of intersection genes between group1 and the blue and red modules, XTH9 (xyloglucan endotransglucosylas), LIP2 (octanoyltransferase), EXL1 (GDSL esterase/lipase), LOX5 (probable linoleate 9s-lipoxygenase), PAL3 (phenylalanine ammonia lyase) and FBP (F-box protein). Error bars represent the standard deviations of the mean in three replicates, and for each stage, * stand for significant differences between seasonal drought treatment and control at p < 0.05. G Correlation analysis of the expression profiles in qRT-PCR and anthocyanin content (*, ** and *** represent significant differences at p < 0.05, p < 0.01 and p < 0.001, respectively between the two sets of data)