Background Little is known about bridging clonal heterogeneity into the resistance of immune checkpoint inhibitors (ICI). Recent reports showed that excluded tumor-infiltrating lymphocytes (TIL) into stroma assessed by an artificial intelligence (AI)-powered spatial TIL analyzer, Lunit SCOPE IO, was related to loss-of-heterozygosity of HLA genes which would be one of crucial resistance pathways of ICI.¹ In the current study, we hypothesized that Immune-excluded phenotype called by Lunit SCOPE IO would be related to clonal heterogeneity resulted from genome-wide accidents during early carcinogenesis which may cause an improper targeting of TIL for diverse clones with multiple genomic aberrations.

Methods For spatial TIL analysis, we applied Lunit SCOPE IO¹ which automatically detects TIL and segmentizes cancer area and stroma, then it classified Immune phenotype of 1 mm2-sized grid in H&E image. Inflamed score or Immune-excluded score were defined as the proportion of Inflamed phenotype, which is high intra-tumoral TIL density, or Immune-excluded phenotype, which is exclusively high TIL density only in stroma, within a whole-slide image, respectively. We evaluated the correlation of Immune phenotype with APOBEC mutational signature by single-base substitution (SBS) signature 2 and/or SBS13,² whole-genome doubling, and subclonal genome fraction which reflects intra-tumoral heterogeneity,³ and clusters of T cell receptor (TCR) repertoire ⁴ derived from previous reports of The Cancer Genome Atlas (TCGA), consists of 7,467 tumor samples from 22 cancer types.

Results In the TCGA pan-carcinoma database, APOBEC mutational signature was significantly correlated with increased ratio of cancer stroma to cancer epithelium (median 0.866 vs 1.19, fold change +37.4%), and increased TIL density in cancer stroma (median 558 vs 764 / mm2, fold change +36.9%), but it was not correlated with intra-tumoral TIL density (median 63 vs 59 / mm2, fold change -6.3%). Interestingly, Immune-excluded score (IES) called by Lunit SCOPE IO was positively correlated with APOBEC mutational signature as well as expression levels of APOBEC1, APOBEC3A, and APOBEC3B, whole-genome doubling, and subclonal genome fraction, respectively, while Inflamed score (IS) or immune cytolytic activity (GZMA and PRF1 expressions) was negatively or not significantly correlated to those variables (table 1). TCR repertoire was expanded in the tumor samples with high IS (spearman rho = 0.279), but it was not increased in those with high IES (spearman rho = -0.0595).

Conclusions There is a significant correlation between distinct TIL deposition in stroma, or Immune-excluded phenotype, with APOBEC-attributed clonal expansion of cancer, without proper expansion of TCR repertoire.

References

1. Ock CY, Park C, Paeng K, Yoo D, Kim S, Park S, Lee SH, Mok T, Bang YJ. Artificial intelligence-powered spatial analysis of tumor-infiltrating lymphocytes reveals distinct genomic profile of immune excluded phenotype in pan-carcinoma. Cancer Res 2021;81(Supp 13):1908.

2. Alexandrov LB, Kim J, Haradhvala NJ, Huang MN, Tian Ng AW, Wu Y, Boot A, Covington KR, Gordenin DA, Bergstrom EN, Islam SMA, Lopez-Bigas N, Klimczak LJ, McPherson JR, Morganella S, Sabarinathan R, Wheeler DA, Mustonen V, PCAWG Mutational Signatures Working Group, Getz G, Rozen SG, Stratton MR, PCAWG Consortium. The repertoire of mutational signatures in human cancer. Nature 2020;578(7793):94–101.

3. Taylor AM, Shih J, Ha G, Gao GF, Zhang X, Berger AC, Schumacher SE, Wang C, Hu H, Liu J, Lazar AJ, Cancer Genome Atlas Research Network, Cherniack AD, Beroukhim R, Meyerson M. Genomic and functional approaches to understanding cancer aneuploidy. Cancer Cell 2018;33(4):676–689.e3.

4. Zhang H, Liu L, Zhang J, Chen J, Ye J, Shukla S, Qiao J, Zhan X, Chen H, Wu CJ, Fu YX, Li B. Investigation of antigen-specific T-Cell receptor clusters in human cancers. Clin Cancer Res 2020;26(6):1359–1371.