Using the r package spatstat to assess inhibitory effects of microregional hypoxia on the infiltration of cancers of the head and neck region by cytotoxic t lymphocytes

(1) Background: The immune system has physiological antitumor activity, which is partially mediated by cytotoxic T lymphocytes (CTL). Tumor hypoxia, which is highly prevalent in cancers of the head and neck region, has been hypothesized to inhibit the infiltration of tumors by CTL. In situ data validating this concept have so far been based solely upon the visual assessment of the distribution of CTL. Here, we have established a set of spatial statistical tools to address this problem mathematically and tested their performance. (2) Patients and Methods: We have analyzed regions of interest (ROI) of 22 specimens of cancers of the head and neck region after 4-plex immunofluorescence staining and whole-slide scanning. Single cell-based segmentation was carried out in QuPath. Specimens were analyzed with the endpoints clustering and interactions between CTL, normoxic, and hypoxic tumor areas, both visually and using spatial statistical tools implemented in the R package Spatstat. (3) Results: Visual assessment suggested clustering of CTL in all instances. The visual analysis also suggested an inhibitory effect between hypoxic tumor areas and CTL in a minority of the whole-slide scans (9 of 22, 41%). Conversely, the objective mathematical analysis in Spatstat demonstrated statistically significant inhibitory interactions between hypoxia and CTL accumulation in a substantially higher number of specimens (16 of 22, 73%). It showed a similar trend in all but one of the remaining samples. (4) Conclusion: Our findings provide non-obvious but statistically rigorous evidence of inhibition of CTL infiltration into hypoxic tumor subregions of cancers of the head and neck. Importantly, these shielded sites may be the origin of tumor recurrences. We provide the methodology for the transfer of our statistical approach to similar questions. We discuss why versions of the Kcross and pcf.cross functions may be the methods of choice among the repertoire of statistical tests in Spatstat for this type of analysis.