Construction of an Immune Cell Infiltration Score to Evaluate the Prognosis and Therapeutic Efficacy of Ovarian Cancer Patients

Radiogenomic biomarkers for immunotherapy in glioblastoma: A systematic review of magnetic resonance imaging studies

Background

Immunotherapy is an effective "precision medicine" treatment for several cancers. Imaging signatures of the underlying genome (radiogenomics) in glioblastoma patients may serve as preoperative biomarkers of the tumor-host immune apparatus. Validated biomarkers would have the potential to stratify patients during immunotherapy clinical trials, and if trials are beneficial, facilitate personalized neo-adjuvant treatment. The increased use of whole genome sequencing data, and the advances in bioinformatics and machine learning make such developments plausible. We performed a systematic review to determine the extent of development and validation of immune-related radiogenomic biomarkers for glioblastoma.

Methods

A systematic review was performed following PRISMA guidelines using the PubMed, Medline, and Embase databases. Qualitative analysis was performed by incorporating the QUADAS 2 tool and CLAIM checklist. PROSPERO registered: CRD42022340968. Extracted data were insufficiently homogenous to perform a meta-analysis.

Results

Nine studies, all retrospective, were included. Biomarkers extracted from magnetic resonance imaging volumes of interest included apparent diffusion coefficient values, relative cerebral blood volume values, and image-derived features. These biomarkers correlated with genomic markers from tumor cells or immune cells or with patient survival. The majority of studies had a high risk of bias and applicability concerns regarding the index test performed.

Conclusions

Radiogenomic immune biomarkers have the potential to provide early treatment options to patients with glioblastoma. Targeted immunotherapy, stratified by these biomarkers, has the potential to allow individualized neo-adjuvant precision treatment options in clinical trials. However, there are no prospective studies validating these biomarkers, and interpretation is limited due to study bias with little evidence of generalizability.

A modified immune cell infiltration score achieves ideal stratification for CD8+ T cell efficacy and immunotherapy benefit in hepatocellular carcinoma

Immunotherapy, which aims to enhance the function of T cells, has emerged as a novel therapeutic approach for hepatocellular carcinoma (HCC). Nevertheless, the clinical utility of using flow cytometry to assess immune cell infiltration (ICI) is hindered by its cumbersome procedures, prompting the need for more accessible methods. Here, we acquired gene expression profiles and survival data of HCC from TCGA and GSE10186 datasets. The patients were categorized into two clusters of ICI, and a set of 11 characteristic genes responsible for the differentiation performance of these ICI clusters were identified. Subsequently, we successfully developed a modified ICI score (mICIS) by utilizing the expression levels of these genes. The efficacy of our mICIS was confirmed via mass cytometry, flow cytometry, and immunohistochemistry. Our research indicated that the favorable overall survival (OS) rate could be attributed to the improved function of anti-tumor leukocytes rather than their infiltration. Furthermore, we observed that the low score group exhibited lower expression levels of T-cell exhaustion-associated genes, which was confirmed in both HCC tissues from patients and mice, which demonstrated that the benefits of the low scores were due to enhanced active/cytotoxic CD8+ T cells and reduced exhausted CD8+ T cells. Additionally, our mICIS stratified the benefits derived from immunotherapies. Lastly, we observed a misalignment between CD8+ T-cell infiltration and function in HCC. In summary, our mICIS demonstrated proficiency in assessing the OS rate of HCC and offering significant stratified data pertaining to distinct responses to immunotherapy.

Characterization of the immune cell function landscape in head and neck squamous carcinoma to assist in prognosis prediction and immunotherapy

BACKGROUND

The malignant characteristics of cancer depend not only on intrinsic properties of cancer cells but also on the functions of infiltrating immune cells. In this study, we aimed to investigate the functional landscape of immune cells in head and neck squamous cell carcinoma (HNSCC).

METHODS

We employed single-sample gene set enrichment analysis to examine the immunophenotypes of HNSCC based on 29 immune cell functions (ICFs) in TCGA and GSE65858 datasets. We analyzed the clinical features, immune microenvironment, molecular profiles, and biological processes. Additionally, we developed and validated an ICF-based risk score for personalized prognosis prediction. We confirmed the value of the ICF score in our cohort using qRT-PCR and immunohistochemistry. Molecular docking was used to predict potential compounds for immunotherapy.

RESULTS

Three immunophenotypes (Immune-L, Immune-M, and Immune-H) were identified in 769 HNSCC samples. The characteristics of Immune-H were consistent with a "Hot" tumor, Immune-L was similar to a "Cold" tumor, and Immune-M exhibited intermediate features. The ICF risk score was associated with immune checkpoints, infiltrating immune cells, tumor mutation burden, and sensitivities to targeted/chemotherapeutic agents. Gene set variation analysis implicated the involvement of metabolic reprogramming pathways in the high-risk group. The combination of "Tumor Immune Dysfunction and Exclusion" and "Immunophenoscore" algorithms indicated that the low-risk group had a higher likelihood of benefiting from immunotherapy. Finally, we identified Eltrombopag and other compounds that may be beneficial for HNSCC immunotherapy.

CONCLUSION

Our study provides a novel perspective on the tumor microenvironment of HNSCC, aiding in the understanding of HNSCC heterogeneity and the development of personalized/precision medicine.

A methylation- and immune-related lncRNA signature to predict ovarian cancer outcome and uncover mechanisms of chemoresistance

Tumor-associated lncRNAs regulated by epigenetic modification switches mediate immune escape and chemoresistance in ovarian cancer (OC). However, the underlying mechanisms and concrete targets have not been systematically elucidated. Here, we discovered that methylation modifications played a significant role in regulating immune cell infiltration and sensitizing OC to chemotherapy by modulating immune-related lncRNAs (irlncRNAs), which represent tumor immune status. Through deep analysis of the TCGA database, a prognostic risk model incorporating four methylation-related lncRNAs (mrlncRNAs) and irlncRNAs was constructed. Twenty-one mrlncRNA/irlncRNA pairs were identified that were significantly related to the overall survival (OS) of OC patients. Subsequently, we selected four lncRNAs to construct a risk signature predictive of OS and indicative of OC immune infiltration, and verified the robustness of the risk signature in an internal validation set. The risk score was an independent prognostic factor for OC prognosis, which was demonstrated via multifactorial Cox regression analysis and nomogram. Moreover, risk scores were negatively related to the expression of CD274, CTLA4, ICOS, LAG3, PDCD1, and PDCD1LG2 and negatively correlated with CD8+, CD4+, and Treg tumor-infiltrating immune cells. In addition, a high-risk score was associated with a higher IC50 value for cisplatin, which was associated with a significantly worse clinical outcome. Next, a competing endogenous RNA (ceRNA) network and a signaling pathway controlling the infiltration of CD8+ T cells were explored based on the lncRNA model, which suggested a potential therapeutic target for immunotherapy. Overall, this study constructed a prognostic model by pairing mrlncRNAs and irlncRNAs and revealed the critical role of the FTO/RP5-991G20.1/hsa-miR-1976/MEIS1 signaling pathway in regulating immune function and enhancing anticancer therapy.

An immune-related biomarker index for predicting the effectiveness of immunotherapy and prognosis in hepatocellular carcinoma

OBJECTIVE

Currently, there are no recognized biomarkers for predicting the immunotherapy response and prognosis of hepatocellular carcinoma (HCC). This study aimed to establish an immune-related gene prognostic index (IRGPI) for HCC, and to investigate the clinical, immune, molecular, and microenvironmental characteristics of the IRGPI subgroups, as well as their impact on the effectiveness of immune checkpoint inhibitors (ICIs) therapy and patients' prognosis.

METHODS

We analyzed the LIHC dataset (n = 424) from the The Cancer Genome Atlas (TCGA) database and the GSE10140 dataset (n = 84) from the Gene Expression Omnibus (GEO) database using weighted gene co-expression network analysis (WGCNA) and univariate/multivariate Cox regression analysis to identify immune-related hub genes with prognostic significance. Subsequently, The IRGPI was then established with these special genes obtained, and the molecular, immune, and clinicopathological characteristics of the IRGPI subgroups, along with their predictive role in ICIs treatment and HCC prognosis, were investigated.

RESULTS

The IRGPI was composed of nine genes, namely CHGA, GAL, CCR3, MMP7, STC1, UCN, OXT, SOCS2, and GCG. The IRGPI-high group exhibited a worse prognosis in both the TCGA and GEO databases compared to the IRGPI-low group. The IRGPI-high group was primarily associated with adaptive immune response and cell-cell interaction pathways and exhibited a higher frequency of gene mutations (such as TP53 and CTNNB1), higher expression of PD-L1 and CTLA4, a higher proportion of macrophages M0 and follicular helper T cells, and a higher APC_co_inhibition and T_cell_co-inhibition immune score. Furthermore, the IRGPI-high group was associated with worse immune subtypes, clinicopathological characteristics, immunotherapy response, and clinical prognosis.

CONCLUSION

IRGPI is a biomarker with significant potential for predicting the immunotherapy response and prognosis of HCC patients, and is closely related to the immunosuppressive microenvironment and poorer clinicopathological characteristics.

A lipid metabolism-related gene signature reveals dynamic immune infiltration of the colorectal adenoma-carcinoma sequence

BACKGROUND

Lipid metabolism-related genes (LMRGs) have been reported to be correlated with the immune infiltration of colorectal cancer (CRC). This study aimed to investigate the immune infiltration characteristics along the colorectal adenoma-carcinoma sequence (ACS) based on LMRGs.

METHODS

Gene expression data of colorectal adenoma and carcinoma samples were obtained from the public databases. The "limma" package was applied to determine the differentially expressed LMRGs. Unsupervised consensus clustering was used to cluster colorectal samples. The features of the tumor microenvironment were analyzed by the "ESTIMATE", "GSVA", and "TIDE" algorithms.

RESULTS

The expression of 149 differentially expressed LMRGs was defined as the LMRG signature. Based on this signature, the adenoma and carcinoma samples were divided into three clusters. Unexpectedly, these sequential clusters showed a directional relationship and collectively constituted the progressive course of colorectal ACS. Interestingly, the LMRG signature revealed that adenoma progression was accompanied by a progressive loss of immune infiltration and a stepwise establishment of a cold microenvironment, but carcinoma progression was characterized by a progressive gain of immune infiltration and a gradual establishment of a hot microenvironment.

CONCLUSIONS

The LMRG signature reveals dynamic immune infiltration along the colorectal ACS, which substantially changes the understanding of the tumor microenvironment of CRC carcinogenesis and provides novel insight into the role of lipid metabolism in this process.

Machine learning-based integration develops an immune-related risk model for predicting prognosis of high-grade serous ovarian cancer and providing therapeutic strategies

Background

High-grade serous ovarian cancer (HGSOC) is a highly lethal gynecological cancer that requires accurate prognostic models and personalized treatment strategies. The tumor microenvironment (TME) is crucial for disease progression and treatment. Machine learning-based integration is a powerful tool for identifying predictive biomarkers and developing prognostic models. Hence, an immune-related risk model developed using machine learning-based integration could improve prognostic prediction and guide personalized treatment for HGSOC.

Methods

During the bioinformatic study in HGSOC, we performed (i) consensus clustering to identify immune subtypes based on signatures of immune and stromal cells, (ii) differentially expressed genes and univariate Cox regression analysis to derive TME- and prognosis-related genes, (iii) machine learning-based procedures constructed by ten independent machine learning algorithms to screen and construct a TME-related risk score (TMErisk), and (iv) evaluation of the effect of TMErisk on the deconstruction of TME, indication of genomic instability, and guidance of immunotherapy and chemotherapy.

Results

We identified two different immune microenvironment phenotypes and a robust and clinically practicable prognostic scoring system. TMErisk demonstrated superior performance over most clinical features and other published signatures in predicting HGSOC prognosis across cohorts. The low TMErisk group with a notably favorable prognosis was characterized by BRCA1 mutation, activation of immunity, and a better immune response. Conversely, the high TMErisk group was significantly associated with C-X-C motif chemokine ligands deletion and carcinogenic activation pathways. Additionally, low TMErisk group patients were more responsive to eleven candidate agents.

Conclusion

Our study developed a novel immune-related risk model that predicts the prognosis of ovarian cancer patients using machine learning-based integration. Additionally, the study not only depicts the diversity of cell components in the TME of HGSOC but also guides the development of potential therapeutic techniques for addressing tumor immunosuppression and enhancing the response to cancer therapy.

Profiling specific cell populations within the inflammatory tumor microenvironment by oscillating-gradient diffusion-weighted MRI

BACKGROUND

The inflammatory tumor microenvironment (TME) is formed by various immune cells, being closely associated with tumorigenesis. Especially, the interaction between tumor-infiltrating T-cells and macrophages has a crucial impact on tumor progression and metastatic spread. The purpose of this study was to investigate whether oscillating-gradient diffusion-weighted MRI (OGSE-DWI) enables a cell size-based discrimination between different cell populations of the TME.

METHODS

Sine-shaped OGSE-DWI was combined with the Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion (IMPULSED) approach to measure microscale diffusion distances, here relating to cell sizes. The accuracy of IMPULSED-derived cell radii was evaluated using in vitro spheroid models, consisting of either pure cancer cells, macrophages, or T-cells. Subsequently, in vivo experiments aimed to assess changes within the TME and its specific immune cell composition in syngeneic murine breast cancer models with divergent degrees of malignancy (4T1, 67NR) during tumor progression, clodronate liposome-mediated depletion of macrophages, and immune checkpoint inhibitor (ICI) treatment. Ex vivo analysis of IMPULSED-derived cell radii was conducted by immunohistochemical wheat germ agglutinin staining of cell membranes, while intratumoral immune cell composition was analyzed by CD3 and F4/80 co-staining.

RESULTS

OGSE-DWI detected mean cell radii of 8.8±1.3 µm for 4T1, 8.2±1.4 µm for 67NR, 13.0±1.7 for macrophage, and 3.8±1.8 µm for T-cell spheroids. While T-cell infiltration during progression of 4T1 tumors was observed by decreasing mean cell radii from 9.7±1.0 to 5.0±1.5 µm, increasing amount of intratumoral macrophages during progression of 67NR tumors resulted in increasing mean cell radii from 8.9±1.2 to 12.5±1.1 µm. After macrophage depletion, mean cell radii decreased from 6.3±1.7 to 4.4±0.5 µm. T-cell infiltration after ICI treatment was captured by decreasing mean cell radii in both tumor models, with more pronounced effects in the 67NR tumor model.

CONCLUSIONS

OGSE-DWI provides a versatile tool for non-invasive profiling of the inflammatory TME by assessing the dominating cell type T-cells or macrophages.

Development and verification of a 7-lncRNA prognostic model based on tumor immunity for patients with ovarian cancer

BACKGROUND

Both immune-reaction and lncRNAs play significant roles in the proliferation, invasion, and metastasis of ovarian cancer (OC). In this study, we aimed to construct an immune-related lncRNA risk model for patients with OC.

METHOD

Single sample GSEA (ssGSEA) algorithm was used to analyze the proportion of immune cells in The Cancer Genome Atlas (TCGA) and the hclust algorithm was used to conduct immune typing according to the proportion of immune cells for OC patients. The stromal and immune scores were computed utilizing the ESTIMATE algorithm. Weighted gene co-expression network analysis (WGCNA) and differentially expressed genes (DEGs) analyses were utilized to detect immune cluster-related lncRNAs. The least absolute shrinkage and selection operator (LASSO) regression was conducted for lncRNA selection. The selected lncRNAs were used to construct a prognosis-related risk model, which was then validated in Gene Expression Omnibus (GEO) database and in vitro validation.

RESULTS

We identify two subtypes based on the ssGSEA analysis, high immunity cluster (immunity_H) and low immunity cluster (immunity_L). The proportion of patients in immunity_H cluster was significantly higher than that in immunity_L cluster. The ESTIMATE related scores are relative high in immunity_H group. Through WGCNA and LASSO analyses, we identified 141 immune cluster-related lncRNAs and found that these genes were mainly enriched in autophagy. A signature consisting of 7 lncRNAs, including AL391832.3, LINC00892, LINC02207, LINC02416, PSMB8.AS1, AC078788.1 and AC104971.3, were selected as the basis for classifying patients into high- and low-risk groups. Survival analysis and area under the ROC curve (AUC) of the signature pointed out that this risk model had high accuracy in predicting the prognosis of patients with OC. We also conducted the drug sensitive prediction and found that rapamycin outperformed in patient with high risk score. In vitro experiments also confirmed our prediction.

CONCLUSIONS

We identified 7 immune-related prognostic lncRNAs that effectively predicted survival in OC patients. These findings may offer a valuable indicator for clinical stratification management and personalized therapeutic options for these patients.

A neutrophil extracellular traps-associated lncRNA signature predicts the clinical outcomes in patients with lung adenocarcinoma

Backgrounds: Neutrophil extracellular traps (NETs) play an important role in the occurrence, metastasis, and immune escape of cancers. We aim to investigate Long non-coding RNAs (lncRNAs) that are correlated to NETs to find some potentially useful biomarkers for lung adenocarcinoma (LUAD), and to explore their correlations with immunotherapy and chemotherapy, as well as the tumor microenvironment. Methods: Based on the The Cancer Genome Atlas (TCGA) database, we identified the prognosis-related lncRNAs which are associated with NETs using cox regression. The patients were then separated into two clusters based on the expression of NETs-associated lncRNAs to perform tumor microenvironment analysis and immune-checkpoint analysis. Least absolute shrinkage and selection operator (LASSO) regression was then performed to establish a prognostic signature. Furthermore, nomogram analysis, tumor mutation burden analysis, immune infiltration analysis, as well as drug sensitivity analysis were performed to test the signature. Results: Using univariate cox regression, we found 10 NETs-associated lncRNAs that are associated with the outcomes of LUAD patients. Also, further analysis which separated the patients into 2 clusters showed that the 10 lncRNAs had significant correlations with the tumor microenvironment. Using LASSO regression, we finally constructed a signature to predict the outcomes of the patients based on 4 NETs-associated lncRNAs. The 4 NETs-associated lncRNAs were namely SIRLNT, AL365181.3, FAM83A-AS1, and AJ003147.2. Using Kaplan-Meier (K-M) analysis, we found that the risk model was strongly associated with the survival outcomes of the patients both in the training group and in the validation group 1 and 2 (p < 0.001, p = 0.026, and p < 0.01). Using receiver operating characteristic (ROC) curve, we tested the sensitivity combined with the specificity of the model and found that the risk model had a satisfactory level of 1-year, 3-year, and 5-year concordance index (C-index) (C = 0.661 in the training group, C = 0.679 in validation group 1, C = 0.692 in validation group 2). We also explored the immune microenvironment and immune checkpoint correlation of the risk model and found some significant results. Conclusion: We constructed a NETs-associated lncRNA signature to predict the outcome of patients with LUAD, which is associated with immunephenoscores and immune checkpoint-gene expression.

Identification of Tumor Microenvironment Scoring Scheme Based on Bioinformatics Analysis of Immune Cell Infiltration Pattern of Ovarian Cancer

Background

Tumor microenvironment (TME) is the crucial mediator of tumor progression, and the TME model based on immune cell infiltration to characterize ovarian cancer is considered to be a promising strategy.

Methods

Sample data of three ovarian cancer cohorts were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The scores of 22 kinds of immune cells were calculated based on CIBERSORT, and the TME clusters (TMECs) of ovarian cancer was determined by ConsensusClusterPlus. Genomic subtype was identified by non-negative matrix factorization (NMF). A TME scoring scheme was constructed using k-means algorithm and principal component analysis (PCA) to quantify the TME infiltration pattern of individuals.

Results

Four TME subtypes of ovarian cancer samples were defined: TMEC1, TMEC2, TMEC3, and TMEC4. There were also significant differences in overall survival (OS) among the four TMEC, and the OS of TMEC3 was the longest. The difference analysis of TMEC3 and the other three TMECs respectively identified the DEGs and took the intersection, and 585 DEGs were obtained. Two genomic subtypes were identified by NMF based on the expression of 585 genes, which were called GeneC1 and GeneC2. The TME scoring scheme constructed by k-means and PCA algorithm was used to calculate the TME score of ovarian cancer in TCGA. High-TME score was significantly correlated with shorter survival time, older age, lower immunoactivated molecules, and immune checkpoint gene expression.

Conclusions

This study highlighted the complexity and diversity of TME infiltration patterns in ovarian cancer and constructed a set of TME scoring scheme to reveal TME infiltration patterns and provided new insights into the landscape of TME.

Exploration of the Immunotyping Landscape and Immune Infiltration-Related Prognostic Markers in Ovarian Cancer Patients

Background

Increasing evidence indicates that immune cell infiltration (ICI) affects the prognosis of multiple cancers. This study aims to explore the immunotypes and ICI-related biomarkers in ovarian cancer.

Methods

The ICI levels were quantified with the CIBERSORT and ESTIMATE algorithms. The unsupervised consensus clustering method determined immunotypes based on the ICI profiles. Characteristic genes were identified with the Boruta algorithm. Then, the ICI score, a novel prognostic marker, was generated with the principal component analysis of the characteristic genes. The relationships between the ICI scores and clinical features were revealed. Further, an ICI signature was integrated after the univariate Cox, lasso, and stepwise regression analyses. The accuracy and robustness of the model were tested by three independent cohorts. The roles of the model in the immunophenoscores (IPS), tumor immune dysfunction and exclusion (TIDE) scores, and immunotherapy responses were also explored. Finally, risk genes (GBP1P1, TGFBI, PLA2G2D) and immune cell marker genes (CD11B, NOS2, CD206, CD8A) were tested by qRT-PCR in clinical tissues.

Results

Three immunotypes were identified, and ICI scores were generated based on the 75 characteristic genes. CD8 TCR pathways, chemokine-related pathways, and lymphocyte activation were critical to immunophenotyping. Higher ICI scores contributed to better prognoses. An independent prognostic factor, a three-gene signature, was integrated to calculate patients’ risk scores. Higher TIDE scores, lower ICI scores, lower IPS, lower immunotherapy responses, and worse prognoses were revealed in high-risk patients. Macrophage polarization and CD8 T cell infiltration were indicated to play potentially important roles in the development of ovarian cancer in the clinical validation cohort.

Conclusions

Our study characterized the immunotyping landscape and provided novel immune infiltration-related prognostic markers in ovarian cancer.

Identification of an Immune Gene-Based Cisplatin Response Model and CD27 as a Therapeutic Target against Cisplatin Resistance for Ovarian Cancer

Objective. Evidence demonstrates that the immune microenvironment is extensively associated with chemotherapy response of ovarian cancer (OV). Herein, this study is aimed at establishing a cisplatin response prediction model for OV on the basis of immune genes. Methods. The expression profiles of cisplatin-sensitive and cisplatin-resistant OV specimens were integrated from multiple public datasets. The abundance scores of 22 immune cells were estimated with CIBERSORT algorithm. Differentially expressed immune genes (DEGs) were determined between cisplatin-sensitive and cisplatin-resistant groups. Thereafter, a cisplatin response model was constructed based on prognostic DEGs with logistic regression analysis. The prediction performance was validated in independent cohorts. The possible relationships between the model and immunotherapy were then assessed. Results. Treg scores were significantly decreased in cisplatin-resistant than cisplatin-sensitive OV specimens, with the opposite results for naïve B cells and activated dendritic cells. Fourteen prognostic DEGs were identified and used to develop a cisplatin-response model. The response scores, estimated by the model, showed favorable performance in discriminating cisplatin-response and nonresponse samples. The response scores also presented significantly negative correlations with three well-known cisplatin-resistant pathways and a positive correlation with the expression of CD274 (PD-L1). Moreover, the decreased CD27 expression was observed in cisplatin-resistant groups, and OV specimens with higher CD27 expressions were more sensitive to cisplatin treatment. Conclusion. Altogether, our findings proposed a cisplatin response prediction model and identified CD27 that might be involved in cisplatin resistance. Further investigations suggested that CD27 could be a promising immunotherapeutic target for cisplatin-resistant subset of OV.

Natural Killer Cells: the Missing Link in Effective Treatment for High-Grade Serous Ovarian Carcinoma

OPINION STATEMENT

Ovarian cancer (OC), especially high-grade serous cancer (HGSC), is a highly heterogeneous malignancy with limited options for curative treatment and a high frequency of relapse. Interactions between OC and the immune system may permit immunoediting and immune escape, and current standard of care therapies can influence immune cell infiltration and function within the tumor microenvironment. Natural killer (NK) cells are involved in cancer immunosurveillance and immunoediting and can be activated by therapy, but deliberate approaches to maximize NK cell reactivity for treatment of HGSC are in their infancy. NK cells may be the ideal target for immunotherapy of HGSC. The diverse functions of NK cells, and their established roles in immunosurveillance, make them attractive candidates for more precise and effective HGSC treatment. NK cells' functional capabilities differ because of variation in receptor expression and genetics, with meaningful impacts on their anticancer activity. Studying HGSC:NK cell interactions will define the features that predict the best outcomes for patients with the disease, but the highly diverse nature of HGSC will likely require combination therapies or approaches to simultaneously target multiple, co-existing features of the tumor to avoid tumor escape and relapse. We expect that the ideal therapy will enable NK cell infiltration and activity, reverse immunosuppression within the tumor microenvironment, and enable effector functions against the diverse subpopulations that comprise HGSC.