Glutathione S-transferases (GSTs) are proteases with multiple physiological functions and play an important role in plant responses to abiotic stresses. Nevertheless, there is a paucity of systematic research on GST genes in Prunus genus. Here, 330 GST genes in four Prunus species were identified for the first time and classified into eight subgroups based on protein sequence and conserved structure, among which Tau subfamily genes had the largest number. The amino acid lengths of GST-encoded proteins in the four species ranged from 66 to 1152 aa, most of which were soluble proteins and located in the cytoplasm and chloroplasts. The GST family was propelled by tandem duplications, yet robust purifying selection constrained its divergence. Conserved motif and domain analysis revealed that the majority of PmGSTs exhibited a highly conserved GST-N structure. The expression pattern of PmGSTs exhibited tissue specificity and spatiotemporal specificity. qRT-PCR validated the transcriptome results and 11 genes were differentially expressed in varieties with different flower and stem colors. In addition, we discovered an anthocyanin-related gene PmGSTF2, which can effectively restore the anthocyanin and proanthocyanidin deficiency-related phenotypes of the Arabidopsis tt19 mutant. Recombinant PmGSTF2 enhanced the water solubility of cyanidin and cyanidin-3-O-glucoside in vitro. Moreover, PmMYBa1 could directly bind to the promoter of PmGSTF2 and activate its expression. The findings revealed that GSTs were preserved in Prunus species and that PmGSTF2 was critical in regulating anthocyanin accumulation.