Interestingly, we also observed a parabolic expression in genes involved in immune exhaustion (e.g. TCF7, LAG3, CTLA4 and PDCD1) (Fig. 4e), which is consistent with the peak of immune landscape remodelling at S2 tumours as shown by our data above. Here, PDCD1 is linked to neoplasm.