Skip to main content

Cytotoxic lymphocytes-related gene ITK from a systematic CRISPR screen could predict prognosis of ovarian cancer patients with distant metastasis

Abstract

Background

Ovarian cancer, a highly metastatic malignancy, has benefited tremendously from advances in modern human genomics. However, the genomic variations related to the metastasis remains unclear.

Methods

We filtered various significant genes (n = 6722) associated with metastasis within a large-scale functional genomic CRISPR/Cas9 knock-out library including 122,756 single guide RNAs, and identified ITK (IL2 Inducible T Cell Kinase) as a potential cancer suppressor gene for ovarian cancer metastasis. Downstream bioinformatic analysis was performed for ITK using public databases.

Results

We found that patients in low-ITK group had poor prognosis and more distant metastasis than those in high-ITK group in TCGA and GEO databases. We also demonstrated that ITK combined with the clinical factors could accurately predict prognosis through multiple Cox regression analysis and ROC analysis. Moreover, alterations correlated with distant metastasis emereged with significantly increased expression in SAMRCD1 in low-ITK group, but CD244 and SOCS1 in high-ITK group. Integrated analysis revealed dysregulated molecular processes including predominantly oncogenic signaling pathways in low-ITK group but immune related pathways in high-ITK group, which suggested ITK might inhibit distant metastasis in ovarian cancer. Furtherly, deconvolution of the cellular composition of all samples validated the close correlation between ITK and immune related function especially for cytotoxic lymphocytes.

Conclusions

Together, these data provide insights into the potential role of ITK, with implications for the future development of tansformative ovarian cancer therapeutics.

Introduction

Ovarian cancer, a heterogeneous neoplasm, could not be screened early and 80% of them were typically diagnosed in the late stage. Globally, a total of 21,750 new cases are diagnosed as ovarian cancer each year, with 13,940 cancer-specific deaths [1,2,3]. Notably, peritoneal spread was reported as the primary metastasis in ovarian cancer [4], in which malignant cells circulated in the peritoneal fluid to produce the later spread in abdominal organs [5]. The previous study reported that the cancer stem cell-like properties play a positive role in metastatic spread [6]. In the past decades, various treatments were applied to improve the high mortality but failed, such as the oral contraceptives and hormone replacement therapy [7, 8], genetic molecular alterations [9, 10], immunotherapy, chemotherapy, and the inhibitor based treatment including Platinum-based chemotherapy and poly (ADP-ribose) polymerase inhibitors [11]. Therefore, a new treatment was urgently to develop to inhibit the spread of cancer cells and then reduce the ovarian cancer-related mortality. Furthermore, the cross-talk between intracellular macrophages and disseminated cancer cells are identified as a new target to reduce metastasis and disease recurrence [12].

Previous studies showed that Interleukin-2-inducible T-cell kinase (ITK) is a member of the TEC family of non-receptor tyrosine kinase, which acts an essential mediator of intracellular signal transduction in both T-cells and natural killer (NK) cells [13,14,15,16]. ITK functions as a downstream signaling between T-cell and NK cell surface receptor and regulates multiple aspect of T-cell development and function [17, 18]. It has been demonstrated that the deficiency of ITK is associated with the malfunction of T cell development, and T cell disorderes related human diseases [19]. The numerous evidences of ITK deficiency related human disease have suggested ITK as a promising therapeutic target for various human disease with ITK based molecule inhibitors such as the aminothiazole based ITK inhibitors [20, 21] for suppression of lung inflammation, benzimidazole based ITK inhibitors [22], aminopyrimidine based ITK inhibitor [23], 3-aminopyride-2-ones based ITK inhibitors [24]. In addition to the involvement in inflammatory response and autoimmune diseases, ITK was also involved in oncogenesis[25]. Mutations in ITK directly cause defects in T-cell signaling pathways [26]. A loss of function mutation of ITK led to the occurrence of Hodgkin and non-Hodgkin lymphoma, mononucleosis, lymphoproliferative disease after infections [27, 28]. Generally, ITK plays an important role in the inflammatory processes [23] and oncogenesis [25], which pave the promising way to develop inhibitors of the ovarian cancer.

This study was to investigate the predicitive role of ITK in the prognosis of patients with ovarian cancer. We collected 6,722 significant genes correlated with metastasis within CRISPR/Cas9 library and identified ITK as a key factor to predict clinical outcomes of patients with ovarian cancer in the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. Integraged analysis showed that dysregulated molecular processes including predominantly oncogenic signaling pathways are enriched in low-ITK group but immune related pathways in high-ITK group. Furtherly, we discovered that ITK was positively correlated with CD244 and SOCS1 (a suppressor for cancer) but negatively correlated with SMARCD1. Additionally, deconvolution of the cellular composition of all samples with MCP-counter and CIBERSORT methods validated the immune related analysis and discovered the close correlation between ITK and cytotoxic lymphocytes. Our results provided a new method for ovarian cancer diagnosis and treatement in the clinical practice.

Methods

Patients and samples

All the clinical and genomic data was collected from TCGA and GEO databases. A total of 887 ovarian cancer samples and 12 normal controls were selected into this study, including 379 patients samples from TCGA, 18 patients and 12 normal samples from GSE38666, 380 patients samples from GSE140082, 110 patients samples from GSE17260, respectively. These samples enrolled into this study were collected between March 2010 and November 2019 (Additional file 1: Table S1).

CRISPR/Cas9 knockout in vivo screening and mouse model development

SKOV3 cells were purchased from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). SKOV3 was cultured in McCoy’s 5A Medium Modified (Biological Industries, Israel). Genome-scale CRISPR knock-out (GeCKO) v2.0 pooled library was purchased from Addgene (Watertown, MA, USA) and amplified as described [29]. The library totally was composed of 122,756 sgRNAs, including 1,000 control sgRNAs, targeting 19,050 genes and 1,864 miRNAs. SKOV3 cells were generated and infected with the library carried by Lentivirus, and the multiplicity of infection (MOI) was around 0.5. Then SKOV3 cells expressing sgRNAs were selected with puromycin, 1 × 106 cells were orthotopically injected into NOD-SCID mice. Tumor growth were monitored in every 3–4 days and mice sacrificed after 40 days when tumor burden was evident or general health conditions were failing. Primary tumors and metastatic nodules from different organs in the peritoneal cavity were dissected and isolated in vitro expansion culture for the second round of orthotopic injection, respectively. After 3 rounds selection in an in vivo metastasis model, primary cells and metastatic cells were collected for subsequent high throughput DNA deep sequencing to identify metastasis-related candidate genes. We used the RNAi Gene Enrichment Ranking (RIGER) P-value to analyze sgRNAs that were significantly enriched in the metastatic group (sgRNAMet) or orthotopic group (sgRNAPri) in ovarian cancer mouse model. The screening criteria were the number of enriched sgRNA targeting each gene ≥ 2, P-value < 0.05 and the normalized enrichment score (NES) < − 1.2.

Female NOD-SCID mice with 6–8 weeks old were purchased from SPF Biotechnology (Beijing, China). Human SKOV3 ovarian cancer cells 1 × 106 cells in 10μL sterile PBS orthotopically injected into left ovary of NOD-SCID mice. Ascites volumes, the numbers of metastatic nodules were measured and total cells in ascites were harvested at 40 days post-inoculation for xenograft model. Metastatic behavior was also determined by Caliper Life Science IVIS Lumina II Spectrum imaging system. After 3 rounds of in vivo screening, the behavior of cells were monitored at Day 8, Day 14 and Day 21 from the time of tumor implantation.

Tumor mutation and biological function analysis

Mutation data was downloaded from cBioPortal for Cancer Genomics (https://www.cbioportal.org/), and then submitted the query regarding “ITK” in the input box of MutationMapper. The lollipop plot was generated according to the unique paitent ID.

To retrieve the biological function of ITK, the GSEA anlaysis (https://www.gsea-msigdb.org/gsea/index.jsp) was performed in ovarian cancer by using TCGA dataset. Interested genes (n = 4,469) were adopted to process by GSEA against the curated gene sets (Hallmark, KEGG and GO) from the molecular signature database (MSigDB) to investigate the significant enrichment of signaling pathways (FDR < 0.05 and minimum gene set size ≥ 3). Then, we used R package tidyverse to visualize these results. Meanwile, we used Metascape tool (https://metascape.org/gp/index.html#/main/step1) to complete the functional network of these significant genes.

We applied two deconvolution methods—CIBERSORT algorithm [30] to estimate the relative cellular fraction of 22 immune cell types, and MCP-counter [31] software to compute the absolute abundance scores of 8 major immune cell types (CD8+ T cells, CD3+ T cells, Natural killer (NK) cells, cytotoxic lymphocytes, B lymphocytes, myeloid dendritic cells, monocytic lineage cells and neutrophils), fibroblasts and endothelial cells. We input log2-transformed genes expression value for these two methods, and used the LM22 leukocyte gene signature as the input gene signature for CIBERSORT. The deconvolution profiles were then clustered hierarchically and compared between low and high-ITK groups.

Statistical analysis

The raw RNA data was downloaded from TCGA data set and then processed by limma package to identify DEGs (|Fold Change|≥ 2 and FDR q-value < 0.05) between regional and distant metastasis subgroup. We defined a mean as the cutoff value in the expression of ITK in each set. We applied Pearson’s chi-square to identify the statistical significance in the basic characteristics. The t test was used to distinguish the distributive difference of ITK expression between low and high groups. The correlation analysis between ITK and the other 4,469 genes was evaluated with Pearson correlation analysis in the TCGA cohort. We applied ROC curves (R package: timeROC) analysis to identify the sensitivity and specificity of ITK and other clinical factors. We quantified the performance of these factors with the area under ROC curve. All the factors were evaluated based on their values from the TCGA and GEO data sets, separately. In the overall survival (OS) analysis, the Kaplan–Meier method was performed to explore the relationship between risk factors and overall survival, and compare the survival curves with log-rank test. Then we used multiple Cox proportional hazards regression analysis to evaluate the independency of this gene. Hazard rations (HRs) and 95% confidence intervals (CIs) were produced in each cohort. We used R software (version 4.0.2) to perform all the statistical analysis, and defined p < 0.05 as the significant threshold.

Results

Screen the significant genes associated with metastasis by using CRISPR/Cas9 library

To investigate the significant genes correlated with metstasis, we established the orthotopic transplant tumor model in NOD-SCID mice with relatively low metastatic capacity human EOC cell SKOV3. We first transfected the human CRISPR knockout library (GeCKO v2.0) into SKOV3 cells, the CRISPR knockout library contains 19,050 sgRNAs that specifically targeting protein encoding genes, along with 1864 miRNAs and 1000 non-targeted sgRNAs [32]. Subsequently, we collected orthotopic tumor tissues and metastatic tissues from the surface of various organs including ascites in the peritoneal cavity, isolated and expanded the culture in vitro for the next round of orthotopic transplant (Fig. 1A). In vivo imaging showed the growth and metastasis of tumors in CRISPR Vec and CRISPR Lib on days 8, 14 and 21 at the third round, it was observed that on days 14 and 21, NOD-SCID mice in CRISPR Lib developed significant metastasis on the ventral side, whereas intraperitoneal tumors were confined in orthotopic and did not metastasize (Fig. 1B). Furthermore, we observed the significant changes of Photon Flux in CRISPR Lib and CRISPR Vec at day 14 (p < 0.01) and day 21 (p < 0.001; Fig. 1C), which indicated that SKOV3 cells transduced with the human CRISPR knockout library had strong ability of tumor metastasis compared with CRISPR Vec. The numbers of metastatic nodules on the surfaces of intestines were shown in Additional file 2: Figure S1A. Bioluminescence indicated the presence of tumors from the primary site in CRISPR Vec and presence of tumors from the primary site and intestines in CRISPR Lib (Additional file 2: Figure S1B). After 3 rounds of in vivo screening, highly metastatic SKOV3 cells were performed to high-throughput sgRNA ligbrary sequencing. The RIGER P analysis was used to analyze the data of sgRNA library sequencing, As a result, a total of 91 genes were identified as inhibited or promoted factors for the ovarian tumor metastasis in the established model (Fig. 1D). Notebably, ITK was one of the the top three genes to inhibit metastasis in the mouse model of ovarianS cancer. These results demonstrated that ITK might be an important tumor supressor gene for ovarian cancer metastasis.

Fig. 1
figure1

Establishment of orthotopic transplant tumor model in NOD-SCID mice and high-throughput screening of genes associated with metastasis using CRISPR/Cas9 library, A Flow chart of ovarian cancer orthotopic mouse model construction. B In vivo imaging of CRISPR-Lib and CRISPR-Vec at different stage, n = 18. C Photo Flux change of CRISPR-Lib and CRISPR-Vec,** p < 0.01, *** p < 0.001. D Volcano plot of metastasis-related candidate genes. Red dots represent inhibit genes in ovarian tumor metastasis model

ITK expression was significantly correlated with distant metastasis and prognosis of patients with ovarian cancer

To explore the correlation between the distant metastasis and the prognosis of ovarian cancer patients, we compared 91 significant genes from DNA-seq filtration with 217 tumor specific genes with significance between regional and distant metastaisis subgroup from TCGA cohort and obtained one overlapped gene ITK (Fig. 2A and Additional file 3: Table S2). And then we mapped mutations on a linear protein and its domanis, and found 2 truncating mutations (X238 splice and X505 splice) and 6 missense mutations, of which 4 shallow deletion were produced following with these mutations (Fig. 2B). Meanwhile, we discovered patients with these serious DNA mutation were almost low-ITK expression, un-regional metastasis and dead with tumor (Fig. 2B), which could partially explain the low levels of ITK expression in tumor tissues and indicate its prognostic role. Then we discovered that ITK was with significantly lower expression (p = 0.036) in tumor (n = 9) than that in normal (n = 7) samples. Similar result was presented between regional (n = 23) and distant (n = 284) metastasis groups (p = 0.0006; Fig. 2C). This result suggested that ITK might inhibit distant metastasis of ovarian cancer. However, these discoveries should be validated by biological experiments.

Fig. 2
figure2

ITK expression was significantly correlated with distant metastasis and prognosis of patients with ovarian cancer. A The correlation between DNA sequences and TCGA data set. B ITK mutation analysis in ovarian cancer patients. C The ITK expression level analysis in GSE38666 and TCGA data sets. D Kaplan–Meier Survival Analysis of ITK expression level in TCGA data set, p = 0.036. E The hazard ratio evaluation with different clinical factor. (Cox regression model, *p < 0.05, ***p < 0.001)

To further explore the prognostic role of ITK in ovarial cancer, we divided these samples into two groups (low and high-ITK group) based on the mean value as a cut off threshold. Then we discovered that patients with low value expression of ITK (median survial months: 25.5 (95%CI: 18.4–28.7)) had higher risk than those with high values (median survial months: 18.4 (95%CI: 16.3–22.2); p = 0.036; Fig. 2D). Also, multivariate analysis with Cox regression including another 5 prognostic factors revealed that ITK expression was significantly associated with OS as a continuous variable in TCGA cohort (HRs = 1.61 (95%CI: 1.02–2.50); p = 0.041; Fig. 2E). These results confirmed the ability of ITK in predicting survival as an independent factor.

Valiation of ITK for survival prediction in GEO data sets

To validate our discoveries, we selected another two GEO data sets (GSE140082 and GSE17260) to evaluate the prognostic power of ITK. Similarly, we classified patients of each cohort into two subgroups (low and high value) based on the mean value as the cutoff point. In GSE140082, patients with low expression values had significantly poor clinical outcomes (median survial months: NA for OS and 28.9 (95%CI:20.0-NA) for PFS) than those with high expression values (median survial months: 39.8 (95%CI: 35.3-NA) for OS and 18.2 (95%CI:16.4–19.8) for PFS; p < 0.05; Fig. 3A, C). Then, multivariate analysis with Cox regression including another 5 prognostic factors reveled that ITK expression was significantly associated with OS or PFS as a continuous variable in this cohort (HRs > 1.80; p < 0.05; Fig. 3B, D). Finally, we generated the similar results as above in GSE17260 (Fig. 3E, F). Integrately, ITK could predict prognosis of patients with ovarian cancer independently as a favorable biomarker.

Fig. 3
figure3

Validation of ITK for survival prediction in GEO data sets. Kaplan–Meier Survival Analysis of ITK expression level in GSE140082 (A), GSE140082 (C), GSE17260 (E). The hazard ratio evaluation with different clinical factor (B, D and F, Cox regression model, *p < 0.05, **P < 0.01, ***p < 0.001)

ITK had strong diagnositic power in the prognostic prediction

To confirm its predictive power in prognosis correlated with metastasis, we performed ROC analysis and computed the area under curve (AUC) value for the significant clinical factors and ITK expression. Although the AUC value of ITK was lower than the clinical factors (CF), ITK combined CF could reach significantly higher AUC values than any single AUC value from CF (AUC > 0.84; p < 0.05; Fig. 4A–C). These data indicated that ITK could add diagnostic power to the clinicophological prognostic features.

Fig. 4
figure4

ROC curves compare the prognostic power of ITK with clinical risk factors in TCGA and GEO data sets. A TCGA. B GSE140082. C GSE17260. P-values show the AUC value of ITK vs. the ones of another features. ROC receiver operator characteristics, AUC area under curves, CI confidence interval

Dysregulated signaling pathways of ITK between two subgroups with different prognosis

Since ITK expression affects the prognosis of ovarian cancer patients with regional and distant metastasis, we explored the correlation between ITK and 4,469 differentially expressed genes, and we discovered that 359 out of 4469 genes were positively correlated with ITK (Pearson r > 0.3) but 81 genes were negatively correlated with ITK (Pearson r < − 0.3; Fig. 5A and Additional file 4: Table S3). Moreover, alterations correlated with distant metastasis emereged with significantly increased expression in SAMRCD1 in low-ITK group, but CD244 and SOCS1 in high-ITK group (p < 0.001; Fig. 5B and Additional file 5: Table S4), which indicated as oncogene or tumor suppressors in numerous cancers and potential therapeutic targets. Integrated analysis revealed dysregulated molecular processes including predominantly oncogenic signaling pathways in low-ITK group but immune related pathways in high-ITK group, which suggested ITK might inhibit distant metastasis in ovarian cancer (Fig. 5C, D). Furtherly, we identified the function and signaling pathways of the three key genes (CD244, SOCS1 and SAMRCD1). As the favorable biomarkers, CD244 and SOCS1, were mainly expressed in the immune pathway. Oppositely, the oncogene, SMARCD1 was mainly expressed in apoptosis pathways (Fig. 5E). These results indicated that patients with poor prognosis predominantly had more oncogenic signaling pathways correlated with metastasis but immune-enriched pathways in low-ITK group.

Fig. 5
figure5

Gene Set Enrichment Analysis (GSEA) to identify key function of ITK. A The correlativity analysis between ITK expression and the TCGA gene data set (Pearso’s correlation test, r > 0.3 or r < 0.3 as threshold value), B The positive correlation genes (CD244 p = 9.1e−10, SOCS1 p = 4.8e−09) and negative correlation (SMARCD1 p = 0.0045) (t test). C The pathways analysis of the high correlated genes. D, E the cumulative distribution function analysis with high ITK level and low ITK level respectively

Deconvolution of the celluar compositon of ovarian cancer specimens revealed strong correlation between ITK and cytotoxic lymphocytes

To identify the immune-related role of ITK in the tumor microenvironment (TME), we applied MCP-counter and CIBERSORT methods to estimate absolute abundance scores of 8 major immune cell types, endothelial cells and fibroblasts (Fig. 6A) and the relative cellular fraction of 22 immune cell types (Fig. 6B). The abundance scores from MCP-counter showed significantly lower abundance of cytotoxic lymphoctyes (q = 1.11E−3.5), monocytic lineage (q = 8.89E−29), CD8 T cells (q = 2.61E−25), NK cells (q = 1.11E−19) and myeloid dendritic cells (q = 3.79E−19) in low-ITK group compared with high-ITK group (Fig. 6C). And based on CIBERSORT method, macrophages.M1(q = 1.22E−14), CD8 T cells (q = 5.88E−7) and tregs (q = 0.019) were also significantly lower in low-ITK group (Fig. 6D). Furtherly, we analyzed all the markers in these two alogrithms corresponding to different immune cells between low and high-ITK groups, and we found that CD8A, EOMES, KLRC4 and KLRD1, expressed in cytotoxic lymphocytes were strong correlation with ITK in the high-ITK group (Fig. 6E, F). These analysis suggested that ITK, as a favorable factor of TME, played an important role in the immune function of cytoxic lymphocytes.

Fig. 6
figure6

Specific immune landscape. A, B The yellow-black heatmaps showing quantification of ITK expression level and immune cells computed by MCP-counter and CIBERSORT-counter method. C The immune cells in high MCP counter scores (cytotoxic lymphocytes q = 1.11E−3.5, monocytic lineage q = 8.89E−29, CD8 T cells q = 2.61E−25, NK cells q = 1.11E−19 and myeloid dendritic cells q = 3.79E−19), D The immune cells in high CIBERSORT counter scores (M1 q = 1.22E−14, CD8 T cells q = 5.88E−7 and tregs q = 0.019), E, F The correlativity analysis between ITK expression and the marker of different immune cells (Pearso’s correlation test)

Discussion

Ovarian cancer has been extensively studied so far [33], yet the role of ITK in ovarian cancer has not been elucidated in many publications. However, an increasing number of ovarian cancer patients with metastasis from the public data sets had a poor prognosis[34, 35]. Therefore, we provided a landscape of ITK in ovarian cancer from the mouse model to patients, which could help us understand the correlation between ITK expression and the prognosis of patients with ovarian cancer.

The mouse model of ovarian cancer metastasis was successfully constructed by the CRISPR genome wide knock-out library, and the genes were sequenced through RIGER method. Our data showed that ITK expression was significantly inhibited in the metastatic mouse model, which suggested its essential role in the occurrence and development of ovarian cancer metastasis. Then we explored the correlation between low ITK expression and its DNA mutation. We observed two truncating mutations (X238 and X505) and four shallow deletions (D202A, L254V, A261P, and X505) involved in the mutated ITK. The multi-locus and multi-modal mutations may lead to being a low expression of ITK in ovarian cancer [36]. Next, we investigated whether the alteration of ITK expression in the ovarian cancer metastasis mouse model was consistent with that in the ovarian cancer patients. We analyzed and validated these results in TCGA as well as GSE38666 data set. The expression of ITK was significantly lower in patients with ovarian cancer than normal samples, especially in the patients with un-regional metastasis. Moreover, these patients with metastasis had poor prognosis in the low-ITK group than those in high-ITK group. All the data suggested that ITK was a tumor suppressor gene, especially a gene that inhibited cancer metastasis both in mouse models and patients. There was a strong correlation between the ITK expression and the prognosis of patients. With the increase of gene expression value, the overall survival rate of patients is gradually increased. Then the multivariable Cox regression model was used to evaluate the independency of ITK in predicting prognosis in TCGA and GEO data sets. The similar results in both data sets showed that the HRs of patients with low ITK expression was more significant than those with high ITK expression, which suggested that the expression of ITK can be used as an independent prognostic factor to predict the prognosis of patients with ovarian cancer. Furtherly, we evaluated the diagnostic power of ITK expression by using ROC curve analysis. And we found that ITK could add the diagnostic power of the clinical factors, although the AUC value of ITK was not very high in each data set.

To more fully assess the function of ITK in ovarian cancer, we furtherly performed correlation analysis and found that the ITK expression was positively correlated with 359 genes (including cancer suppressor gene CD244 and SOCS1, r > 0.3) and negatively associated with 81 genes (including oncogene SMARCD1, r < − 0.3). In a previous study, the receptor encoded by the CD244 gene was thought to modulate NK-cell cytolytic activity and predominantly displayed inhibitory signaling in tumor-associated immune cells. The receptors' activity and expression and their ligands were correlated with tumor progression, prognosis, and inflammatory responses [37]. Also, it has been reported that the expression level of SOCS1 was significantly increased in human ovarian cancer and might function as a diagnostic biomarker [38]. Our results showed that ITK expression was positively correlated with cancer suppressor gene CD244 (p = 9.1e−10) and SOCS1 (p = 4.8e−09) among patients with a favorable prognosis, suggesting that the function of ITK may be similar with CD244 and SOCS1. Moreover, oncogene SMARCD1 was upregulated in ovarian cells exposed to ascites, allowing a non-stimulatory effect on cancer cell migration [39]. The ITK expression was negatively correlated with SMARCD1(p = 0.0045) in patients with ovarian metastases. All these data showed that ITK could act as a cancer suppressor gene in patients with ovarian cancer and had a similar function to CD244 and SOCS1.

We integrated the pathways enrichment analysis between ITK and another related genes. The results revealed that these genes in high-ITK group were positively correlated with the immune pathway, and genes in low-ITK group were associated with apoptosis and oncogenic pathway. This indicated that patients with low ITK expression were with poor prognosis, mainly because of the factors involved in the cancer development. In contrast, patients with high ITK expression were with favorable prognosis, mainly due to have more factors involved in immune pathways. Finally, we continued to verify the ITK expression in immunity. The results showed that high ITK expression was mainly enriched in different immune cells, such as cytotoxic lymphocytes, monocytic lineage, NK cells, and myeloid dendritic cells. Moreover, ITK expression was positively correlated with markers of these immune cells in patients with favorable prognosis. Generally, it suggested that ITK was indeed involved in the immune function and produced sound prognostic effects.

The biological validation is in some way weak in this study. We understand that it is better to reveal the potential ITK transcriptional mechanism by detecting the RNA and protein expression in corresponding cancer cells and tumor tissues. However, we predominantly focused on exploring the clinical analysis with bioinformatics methods to reveal the mechanism of ITK with low expression in the ovarian cancer patients. And the further efforts will be paid to validate the discoveries about the expression and function of ITK with biological experiments in the next study.

In summary, we identified a strong correlation between differential ITK expression and ovarian cancer development and prognosis, which was verified by TCGA and GEO data sets analysis. Moreover, our enrichment analysis results highlighted the prognostic role of ITK in ovarian cancer were involved in immunology function. Novel prognostic markers are urgently required. This research will provide a framework for future study and a factor for the diagnosis of ovarian cancer, which allows us to define the prognosis status of ovarian patients.

Availability of data and materials

All data generated or analyzed during this are included in this article.

References

  1. 1.

    Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.

    Article  Google Scholar 

  2. 2.

    Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.

    Article  Google Scholar 

  3. 3.

    Sant M, Chirlaque-Lopez MD, Agresti R, Sánchez-Pérez MJ, Holleczek B, Bielska-Lasota M, Dimitrova N, Innos K, Katalinic A, Langseth H, Larrañaga N, Rossi S, Siesling S, Minicozzi P. Survival of women with cancers of breast and genital organs in Europe 1999–2007: Results of the EUROCARE-5 study. Eur J Cancer. 2015;51(15):2191–205.

    Article  Google Scholar 

  4. 4.

    Parashar D, Nair B, Geethadevi A, George J, Nair A, Tsaih SW, Kadamberi IP, Gopinadhan Nair GK, Lu Y, Ramchandran R, Uyar DS, Rader JS, Ram PT, Mills GB, Pradeep S, Chaluvally-Raghavan P. Peritoneal spread of ovarian cancer harbors therapeutic vulnerabilities regulated by FOXM1 and EGFR/ERBB2 signaling. Cancer Res. 2020;80(24):5554–68.

    CAS  Article  Google Scholar 

  5. 5.

    McPherson A, Roth A, Laks E, Masud T, Bashashati A, Zhang AW, Ha G, Biele J, Yap D, Wan A, Prentice LM, Khattra J, Smith MA, Nielsen CB, Mullaly SC, Kalloger S, Karnezis A, Shumansky K, Siu C, Rosner J, Chan HL, Ho J, Melnyk N, Senz J, Yang W, Moore R, Mungall AJ, Marra MA, Bouchard-Côté A, Gilks CB, Huntsman DG, McAlpine JN, Aparicio S, Shah SP. Divergent modes of clonal spread and intraperitoneal mixing in high-grade serous ovarian cancer. Nat Genet. 2016;48(7):758–67.

    CAS  Article  Google Scholar 

  6. 6.

    Bapat SA, Mali AM, Koppikar CB, Kurrey NK. Stem and progenitor-like cells contribute to the aggressive behavior of human epithelial ovarian cancer. Cancer Res. 2005;65(8):3025–9.

    CAS  Article  Google Scholar 

  7. 7.

    Havrilesky LJ, Moorman PG, Lowery WJ, Gierisch JM, Coeytaux RR, Urrutia RP, Dinan M, McBroom AJ, Hasselblad V, Sanders GD, Myers ER. Oral contraceptive pills as primary prevention for ovarian cancer: a systematic review and meta-analysis. Obstet Gynecol. 2013;122(1):139–47.

    CAS  Article  Google Scholar 

  8. 8.

    Wentzensen N, Poole EM, Trabert B, White E, Arslan AA, Patel AV, Setiawan VW, Visvanathan K, Weiderpass E, Adami HO, Black A, Bernstein L, Brinton LA, Buring J, Butler LM, Chamosa S, Clendenen TV, Dossus L, Fortner R, Gapstur SM, Gaudet MM, Gram IT, Hartge P, Hoffman-Bolton J, Idahl A, Jones M, Kaaks R, Kirsh V, Koh WP, Lacey JV, Lee IM, Lundin E, Merritt MA, Onland-Moret NC, Peters U, Poynter JN, Rinaldi S, Robien K, Rohan T, Sandler DP, Schairer C, Schouten LJ, Sjöholm LK, Sieri S, Swerdlow A, Tjonneland A, Travis R, Trichopoulou A, van den Brandt PA, Wilkens L, Wolk A, Yang HP, Zeleniuch-Jacquotte A, Tworoger SS. Ovarian Cancer Risk Factors by Histologic Subtype: An Analysis From the Ovarian Cancer Cohort Consortium. J Clin Oncol. 2016;34(24):2888–98.

    CAS  Article  Google Scholar 

  9. 9.

    Network CGAR. Integrated genomic analyses of ovarian carcinoma. Nature. 2011;474(7353):609–15.

    Article  Google Scholar 

  10. 10.

    Patch AM, Christie EL, Etemadmoghadam D, Garsed DW, George J, Fereday S, Nones K, et al. Whole-genome characterization of chemoresistant ovarian cancer. Nature. 2015;521(7553):489–94.

    CAS  Article  Google Scholar 

  11. 11.

    Fong PC, Yap TA, Boss DS, Carden CP, Mergui-Roelvink M, Gourley C, De Greve J, Lubinski J, Shanley S, Messiou C, A’Hern R, Tutt A, Ashworth A, Stone J, Carmichael J, Schellens JH, de Bono JS, Kaye SB. Poly(ADP)-ribose polymerase inhibition: frequent durable responses in BRCA carrier ovarian cancer correlating with platinum-free interval. J Clin Oncol. 2010;28(15):2512–9.

    CAS  Article  Google Scholar 

  12. 12.

    Etzerodt A, Moulin M, Doktor TK, Delfini M, Mossadegh-Keller N, Bajenoff M, Sieweke MH, Moestrup SK, Auphan-Anezin N, Lawrence T. Tissue-resident macrophages in omentum promote metastatic spread of ovarian cancer. J Exp Med. 2020;217:4.

    Article  Google Scholar 

  13. 13.

    Zhong Y, Dong S, Strattan E, Ren L, Butchar JP, Thornton K, Mishra A, Porcu P, Bradshaw JM, Bisconte A, Owens TD, Verner E, Brameld KA, Funk JO, Hill RJ, Johnson AJ, Dubovsky JA. Targeting interleukin-2-inducible T-cell kinase (ITK) and resting lymphocyte kinase (RLK) using a novel covalent inhibitor PRN694. J Biol Chem. 2015;290(10):5960–78.

    CAS  Article  Google Scholar 

  14. 14.

    Felices M, Berg LJ. The Tec kinases Itk and Rlk regulate NKT cell maturation, cytokine production, and survival. J Immunol. 2008;180(5):3007–18.

    CAS  Article  Google Scholar 

  15. 15.

    Khurana D, Arneson LN, Schoon RA, Dick CJ, Leibson PJ. Differential regulation of human NK cell-mediated cytotoxicity by the tyrosine kinase Itk. J Immunol. 2007;178(6):3575–82.

    CAS  Article  Google Scholar 

  16. 16.

    Schaeffer EM, Debnath J, Yap G, McVicar D, Liao XC, Littman DR, Sher A, Varmus HE, Lenardo MJ, Schwartzberg PL. Requirement for Tec kinases Rlk and Itk in T cell receptor signaling and immunity. Science. 1999;284(5414):638–41.

    CAS  Article  Google Scholar 

  17. 17.

    Andreotti AH, Schwartzberg PL, Joseph RE, Berg LJ. T-cell signaling regulated by the Tec family kinase, Itk. Cold Spring Harb Perspect Biol. 2010;2(7):a002287.

    Article  Google Scholar 

  18. 18.

    Huang L, Ye K, McGee MC, Nidetz NF, Elmore JP, Limper CB, Southard TL, Russell DG, August A, Huang W. Interleukin-2-inducible T-cell kinase deficiency impairs early pulmonary protection against. Front Immunol. 2019;10:3103–15.

    CAS  Article  Google Scholar 

  19. 19.

    Ghosh S, Drexler I, Bhatia S, Adler H, Gennery AR, Borkhardt A. Corrigendum: interleukin-2-inducible T-cell kinase deficiency-new patients, new insight? Front Immunol. 2018;9:2197–204.

    Article  Google Scholar 

  20. 20.

    Das J, Furch JA, Liu C, Moquin RV, Lin J, Spergel SH, McIntyre KW, Shuster DJ, O’Day KD, Penhallow B, Hung CY, Doweyko AM, Kamath A, Zhang H, Marathe P, Kanner SB, Lin TA, Dodd JH, Barrish JC, Wityak J. Discovery and SAR of 2-amino-5-(thioaryl)thiazoles as potent and selective Itk inhibitors. Bioorg Med Chem Lett. 2006;16(14):3706–12.

    CAS  Article  Google Scholar 

  21. 21.

    Lin TA, McIntyre KW, Das J, Liu C, O’Day KD, Penhallow B, Hung CY, Whitney GS, Shuster DJ, Yang X, Townsend R, Postelnek J, Spergel SH, Lin J, Moquin RV, Furch JA, Kamath AV, Zhang H, Marathe PH, Perez-Villar JJ, Doweyko A, Killar L, Dodd JH, Barrish JC, Wityak J, Kanner SB. Selective Itk inhibitors block T-cell activation and murine lung inflammation. Biochemistry. 2004;43(34):11056–62.

    CAS  Article  Google Scholar 

  22. 22.

    Riether D, Zindell R, Kowalski JA, Cook BN, Bentzien J, Lombaert SD, Thomson D, Kugler SZ, Skow D, Martin LS, Raymond EL, Khine HH, O’Shea K, Woska JR, Jeanfavre D, Sellati R, Ralph KL, Ahlberg J, Labissiere G, Kashem MA, Pullen SS, Takahashi H. 5-Aminomethylbenzimidazoles as potent ITK antagonists. Bioorg Med Chem Lett. 2009;19(6):1588–91.

    CAS  Article  Google Scholar 

  23. 23.

    von Bonin A, Rausch A, Mengel A, Hitchcock M, Krüger M, von Ahsen O, Merz C, Röse L, Stock C, Martin SF, Leder G, Döcke WD, Asadullah K, Zügel U. Inhibition of the IL-2-inducible tyrosine kinase (Itk) activity: a new concept for the therapy of inflammatory skin diseases. Exp Dermatol. 2011;20(1):41–7.

    Article  Google Scholar 

  24. 24.

    Charrier JD, Miller A, Kay DP, Brenchley G, Twin HC, Collier PN, Ramaya S, Keily SB, Durrant SJ, Knegtel RM, Tanner AJ, Brown K, Curnock AP, Jimenez JM. Discovery and structure-activity relationship of 3-aminopyrid-2-ones as potent and selective interleukin-2 inducible T-cell kinase (Itk) inhibitors. J Med Chem. 2011;54(7):2341–50.

    CAS  Article  Google Scholar 

  25. 25.

    Pechloff K, Holch J, Ferch U, Schweneker M, Brunner K, Kremer M, Sparwasser T, Quintanilla-Martinez L, Zimber-Strobl U, Streubel B, Gewies A, Peschel C, Ruland J. The fusion kinase ITK-SYK mimics a T cell receptor signal and drives oncogenesis in conditional mouse models of peripheral T cell lymphoma. J Exp Med. 2010;207(5):1031–44.

    CAS  Article  Google Scholar 

  26. 26.

    Singleton KL, Gosh M, Dandekar RD, Au-Yeung BB, Ksionda O, Tybulewicz VL, Altman A, Fowell DJ, Wülfing C. Itk controls the spatiotemporal organization of T cell activation. Sci Signal. 2011;4(193):ra66.

    Article  Google Scholar 

  27. 27.

    Parvaneh N, Filipovich AH, Borkhardt A. Primary immunodeficiencies predisposed to Epstein–Barr virus-driven haematological diseases. Br J Haematol. 2013;162(5):573–86.

    CAS  Article  Google Scholar 

  28. 28.

    Tangye SG, Palendira U, Edwards ES. Human immunity against EBV-lessons from the clinic. J Exp Med. 2017;214(2):269–83.

    Article  Google Scholar 

  29. 29.

    Joung J, Konermann S, Gootenberg JS, Abudayyeh OO, Platt RJ, Brigham MD, Sanjana NE, Zhang F. Genome-scale CRISPR-Cas9 knockout and transcriptional activation screening. Nat Protoc. 2017;12(4):828–63.

    CAS  Article  Google Scholar 

  30. 30.

    Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling Tumor Infiltrating Immune Cells with CIBERSORT. Methods Mol Biol. 2018;1711:243–59.

    CAS  Article  Google Scholar 

  31. 31.

    Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, Selves J, Laurent-Puig P, Sautès-Fridman C, Fridman WH, de Reyniès A. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016;17(1):218–38.

    Article  Google Scholar 

  32. 32.

    Shalem O, Sanjana NE, Hartenian E, Shi X, Scott DA, Mikkelson T, Heckl D, Ebert BL, Root DE, Doench JG, Zhang F. Genome-scale CRISPR-Cas9 knockout screening in human cells. Science. 2014;343(6166):84–7.

    CAS  Article  Google Scholar 

  33. 33.

    Nash Z, Menon U. Ovarian cancer screening: Current status and future directions. Best Pract Res Clin Obstet Gynaecol. 2020;65:32–45.

    Article  Google Scholar 

  34. 34.

    Lili LN, Matyunina LV, Walker LD, Benigno BB, McDonald JF. Molecular profiling predicts the existence of two functionally distinct classes of ovarian cancer stroma. Biomed Res Int. 2013;2013:846387.

    Article  Google Scholar 

  35. 35.

    Yoshihara K, Tajima A, Yahata T, Kodama S, Fujiwara H, Suzuki M, Onishi Y, Hatae M, Sueyoshi K, Kudo Y, Kotera K, Masuzaki H, Tashiro H, Katabuchi H, Inoue I, Tanaka K. Gene expression profile for predicting survival in advanced-stage serous ovarian cancer across two independent datasets. PLoS One. 2010;5(3):e9615.

    Article  Google Scholar 

  36. 36.

    Linka RM, Risse SL, Bienemann K, Werner M, Linka Y, Krux F, Synaeve C, Deenen R, Ginzel S, Dvorsky R, Gombert M, Halenius A, Hartig R, Helminen M, Fischer A, Stepensky P, Vettenranta K, Köhrer K, Ahmadian MR, Laws HJ, Fleckenstein B, Jumaa H, Latour S, Schraven B, Borkhardt A. Loss-of-function mutations within the IL-2 inducible kinase ITK in patients with EBV-associated lymphoproliferative diseases. Leukemia. 2012;26(5):963–71.

    CAS  Article  Google Scholar 

  37. 37.

    Buller CW, Mathew PA, Mathew SO. Roles of NK Cell Receptors 2B4 (CD244), CS1 (CD319), and LLT1 (CLEC2D) in Cancer. Cancers (Basel). 2020;12(7):1755–70.

    CAS  Article  Google Scholar 

  38. 38.

    Li X, Kong C, Fan Y, Liu J, Lu W, Meng C, Li A, Zhai A, Yan B, Song W, Han X. Demethylation of SOCS1 mediates its abnormally high expression in ovarian cancer. Oncol Lett. 2019;18(2):1330–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Meunier L, Puiffe ML, Le Page C, Filali-Mouhim A, Chevrette M, Tonin PN, Provencher DM, Mes-Masson AM. Effect of ovarian cancer ascites on cell migration and gene expression in an epithelial ovarian cancer in vitro model. Transl Oncol. 2010;3(4):230–8.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (61802209), the Open Fund of Tianjin Central Hospital of Gynecology Obstetrics/Tianjin Key Laboratory of Human Development and Reproductive Regulation (2020XHY03), a project funded by the China Postdoctoral Science Foundation (BS000098), and the Fundamental Research Funds of the Central Universities, Nankai University (63181329, 63191422).

Author information

Affiliations

Authors

Contributions

MX, YP, and SZ were responsible for the conception and design of the study. MX, JC, SH, and WX contributed to the data analysis. MX and RX conducted the biological validation. MX, YP, and SZ drafted the manuscript. All authors read, and approved the final manuscript.

Corresponding authors

Correspondence to Yongjun Piao or Shuangtao Zhao.

Ethics declarations

Ethics approval and consent to participate

All the experiments involving mice were conducted according to the guidelines established by the Nankai University Animal Care and Use Committee (NUACUC) by skilled experimenters under an approved protocol, which was in accordance with the principles and procedures outlined in the NIH Guide for the Care and Use of Laboratory Animals.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Basical characteristics.

Additional file 2: Figure S1.

A), The numbers of metastatic nodules on the surfaces of intestines in CRISPR-Vec and CRISPR-Lib in the third round screening (White arrow represents tumor in primary site; blue arrow represents metastatic nodules), B) Bioluminescence.

Additional file 3: Table S2.

All genes ttest result between Regional and Unregional groups.

Additional file 4: Table S3.

Correlation analysis among the 4751 genes.

Additional file 5: Table S4.

Key genes correlated with ITK.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Xu, M., Huang, S., Chen, J. et al. Cytotoxic lymphocytes-related gene ITK from a systematic CRISPR screen could predict prognosis of ovarian cancer patients with distant metastasis. J Transl Med 19, 447 (2021). https://doi.org/10.1186/s12967-021-03119-3

Download citation

Keywords

  • Ovarian cancer
  • Metastasis
  • CRISPR screen
  • Bioinformatic analysis