Abstract
NK cells infiltrating Hepatocellular Carcinoma (HCC) may express residency markers such as Integrin Subunit Alpha 1 (CD49a) that have been associated with nurturing functions in the decidua, and characterized by the production of angiogenic factors as well as loss of cytotoxicity. CIBERSORT, a computational analysis method for quantifying cell fractions from bulk tissue gene expression profiles, was used to estimate the infiltrating immune cell composition of the tumor microenvironment from gene expression profiles of a large cohort of 225 HCCs in the public GEO database. Decidual-like CD49a+ NK cells, in addition to another 22 immune cell populations, were characterized and thoroughly investigated so that HCC cell heterogeneity in a large cohort of 225 HCCs from the public GEO database could be studied. An inverse correlation of the expression of CD49a+ NK-cells and CD8(+) T-cells suggested a negative association with clinical outcomes. This result was confirmed in a further validation cohort of 100 HCC patients from The Cancer Genome Atlas, Liver Hepatocellular Carcinoma (TCGA-LIHC). Cox regression analysis did not identify CD49a+ cells as a variable independently associated with survival. However, a more abundant infiltrate of this subset was present in patients at a more advanced pathological and clinical HCC stage. In conclusion, we found that NK cells, with a decidual-like gene expression profile, are enriched in HCC, and their abundance increases not only in tumor size but also at advanced stages of the disease suggesting that these cells play a role in tumor growth. For this reason, these NK cells may represent a possible new target for immunotherapeutic approaches in HCC.