We also used a single sample gene set enrichment analysis (ssGSEA) algorithm to establish an ‘ISG15+ Tregs’ signature and found that the ISG15+ Tregs activity scores in BCBrM tissues were significantly higher than that in paired primary breast cancer tissues (Figure 4O). Here, ISG15 is linked to breast cancer.