Kejun Liu1,2†
Ji Hao
Yang Bu
Bendong Chen- 1Department of Hepatobiliary Surgery, General Hospital of Ningxia Medical University, Yinchuan, China
- 2Ningxia Hepatobiliary and Pancreatic Surgical Diseases Clinical Medical Research Center, Yinchuan, China
- 3School of Clinical Medicine, Ningxia Medical University, Yinchuan, China
- 4Department of Hepatobiliary Surgery, People’s Hospital of Ningxia Hui Autonomous Region, Yinchuan, China
Anillin actin-binding protein (ANLN) is crucially involved in cell proliferation and migration. Moreover, ANLN is significantly in tumor progression in several types of human malignant tumors; however, it remains unclear whether ANLN acts through common molecular pathways within different tumor microenvironments, pathogeneses, prognoses and immunotherapy contexts. Therefore, this study aimed to perform bioinformatics analysis to examine the correlation of ANLN with tumor immune infiltration, immune evasion, tumor progression, immunotherapy, and tumor prognosis. We observed increased ANLN expression in multiple tumors, which could be involved in tumor cell proliferation, migration, infiltration, and prognosis. The level of ANLN methylation and genetic alteration was associated with prognosis in numerous tumors. ANLN facilitates tumor immune evasion through different mechanisms, which involve T-cell exclusion in different cancer types and tumor-infiltrating immune cells in colon adenocarcinoma, kidney renal clear cell carcinoma, liver hepatocellular carcinoma, and prostate adenocarcinoma. Additionally, ANLN is correlated with immune or chemotherapeutic outcomes in malignant cancers. Notably, ANLN expression may be a predictive biomarker for the response to immune checkpoint inhibitors. Taken together, our findings suggest that ANLN can be used as an onco-immunological biomarker and could serve as a hallmark for tumor screening, prognosis, individualized treatment design, and follow-up.
Introduction
Worldwide, cancer remains the leading cause of morbidity and mortality; moreover, it imposes huge health and economic burden on society (Sung et al., 2021). Current therapies mainly include surgery, radiotherapy, and immunotherapy; however, delayed diagnosis or distant metastasis impedes most patients from undergoing radical surgery (Siegel et al., 2021). Recently, immune checkpoint blockade (ICB) has been shown to be efficacious in immunotherapy for various cancer types. Recent advances in ICB therapy, including those targeting programmed cell death protein 1 (PD-1)/PD-ligand 1 (PD-L1) and/or cytotoxic T-lymphocyte-associated protein 4 (CTLA-4), have allowed strong enhancement of anti-tumor immune responses (Hargadon et al., 2018). Although there have been recent positive reports regarding immunotherapy, its effectiveness widely varies among individuals and tumor types. For example, in some patients with cancer, the objective response rate of anti-PD-L1/PD-1 therapy alone is only 10%–30% (Barrueto et al., 2020; Genova et al., 2022). This could be mainly attributed to the fact that current treatments specifically focus on the cancer cells rather than comprehensively addressing the tumor immune microenvironment (TIME).
Tumor cells and the TIME influence each other. Various immune cells, stromal cells, and cytokines inherently present in the TIME can interact with tumor cells. Moreover, regulation of the immune system network has complex interactions with tumor cells and crucially affects tumor development and the response to immunotherapy (Li et al., 2020a). Tumor cells make direct or indirect contact with various TIME components, with continuous and dynamic changes in the surrounding environment that allow them to evade the host immune surveillance (Mao et al., 2021; Liu et al., 2022). Additionally, the TIME is heterogeneous. The composition of tumor-infiltrating immune cells varies according to the cancer type. Therefore, it is important to elucidate the effects of multiple factors in order to accurately predict the overall effect of the TIME on the tumor immune response, which could improve survival (Morad et al., 2021). Moreover, it could inform the identification of predictive biomarkers for the efficacy of ICBs from the TIME perspective. Additionally, elucidating the mechanisms underlying the efficacy of ICBs and tumor-TIME interactions could inform individualized anti-tumor immunotherapy as well as monitoring of treatment efficacy and disease progression.
Anillin actin-binding protein (ANLN), which is a protein-coding gene, is localized on chromosome 7q14.2 and encodes 1,125 amino acid proteins. It is crucially involved in human diseases and multiple cellular functions, including cell growth, cytokinesis, mitosis, and cell migration (Xiao et al., 2020). Through interactions with cytoskeletal components, ANLN is involved in regulating tumor cell growth, invasion, and metastasis (Naydenov et al., 2021). ANLN is overexpressed in different tumor types, including breast cancer (Downs et al., 2019), lung cancer (Suzuki et al., 2005; Deng et al., 2021), pancreatic cancer (Idichi et al., 2017; Wang et al., 2019), liver cancer (Lian et al., 2018), bladder cancer (Zeng et al., 2017), and colorectal cancer (Wang et al., 2016). Moreover, it worsens the prognosis by promoting malignant tumor progression.
However, the pathogenic role of ANLN in various cancers remains unclear. Further, it remains unclear whether ANLN is involved in different tumor immune microenvironments and can influence the therapeutic response through certain common molecular mechanisms. Accordingly, through bioinformatics analysis, we aimed to elucidate the relationship of ANLN with tumor prognosis, immune infiltration, and therapeutic response in order to determine whether ANLN can be used as a biomarker for cancer screening, prognosis, individualized therapy design, and follow-up.
Materials and methods
Data source
We downloaded the clinical information and processed RNA-sequencing data of patients with 33 types of cancer from The Cancer Genome Atlas (TCGA), which is publicly available, and included 9,807 tumor tissues (Blum et al., 2018). Data regarding normal tissue expression were obtained from the Genotype-Tissue Expression (GTEx, gtexportal.org/home/) project (Carithers and Moore, 2015). We used a database of novel junctions identified in TCGA and GTEx to determine the ANLN expression in various types of human cancers and adjacent normal tissues. The UALCAN portal allows the evaluation of epigenetic regulation of gene expression through promoter methylation as well as pan-cancer gene expression analysis (Chandrashekar et al., 2017). We used the “CPTAC Analysis” module provided by the UALCAN portal to compare ANLN expression levels between normal and tumor tissues. Figure 1 shows the study flowchart.
FIGURE 1. The workflow of the study.
The physiological condition and protein interaction of anillin actin-binding protein
The Human Protein Atlas (proteinatlas.org/) web server can map all human proteins using an integration of various omics technologies, including antibody-based imaging, mass spectrometry-based proteomics, transcriptomics, and systems biology (Uhlen et al., 2015; Thul et al., 2017; Uhlen et al., 2017; Karlsson et al., 2021). PROTTER (wlab.ethz.ch/protter/start/) is an online tool for protein-form visualization as well as annotation and prediction of sequence features with experimental proteomics evidence (Omasits et al., 2014). String (cn.string-db.org/) is an online database of known and predicted protein-protein interaction (PPI) networks (confidence score cutoff = 0.7) (Szklarczyk et al., 2021).
Differential anillin actin-binding protein expression in tumor pathological stages and prognostic analysis
Kaplan Meier plotter (kmplot.com/analysis/) is an online tool for analyzing the association of ANLN expression and genetic alterations with therapeutic outcomes in different tumors (Lanczky and Gyorffy, 2021). Based on the median ANLN expression levels, patient samples were allocated to the high and low ANLN expression groups.
The Gene Expression Profiling Interactive Analysis (GEPIA) (gepia.cancer-pku.cn/index.html) is a novel interactive web server that provides customizable functions, including tumor/normal differential expression analysis, profiling according to cancer types or pathological stages, and patient survival analysis. GEPIA uses a standard processing pipeline to analyze the RNA sequencing expression data of 9,736 tumors and 8,587 normal samples from the TCGA and the GTEx projects (Tang et al., 2017). The GEPIA database was used to evaluate differential ANLN expression between tumor tissues and para-tumor tissues, as well as ANLN expression in different tumor pathological stages.
Tumor DNA methylation and gene mutation characteristics
MethSurv (version MethSurv©2017, biit.cs.ut.ee/methsurv/), which is a web tool for survival analysis based on the CpG methylation patterns (Modhukur et al., 2018), was used to analyze DNA methylation sites of ANLN based on the TCGA data. Moreover, we used the Kaplan-Meier test in the MethSurv database to perform overall survival (OS) analysis for each CpG site. Significance was tested using a likelihood ratio test. Density plots were used to visualize cut-off points for the dichotomy of methylation levels (Supplementary Figure S1). The plots were drawn with the β-value as the abscissa and the corresponding sample proportion as the ordinate. The numbers in red indicate the cutoff points used to group patients.
The cBioPortal for Cancer Genomics (cbioportal.org/) is an open-access, open-source resource for interactive exploration of multiple cancer genomics datasets (Cerami et al., 2012). Currently, it stores data regarding the DNA copy number, mRNA and microRNA expression, non-synonymous mutations, protein and phosphoprotein level, DNA methylation, and limited clinical information (Cerami et al., 2012). The ANLN mutation signatures in pan-tumors were investigated on cBioPortal OncoPrint using data obtained from the TCGA Pan-Cancer Atlas (PanCanAtlas) cohort.
Analysis of tumor immune cell infiltration
Tumor IMmune Estimation Resource (TIMER2.0) (timer.- cistrome.org/), which is a web server for analyzing tumor-infiltrating immune cells using TCGA (Li et al., 2020b), was used to analyze the relationship of immune cell infiltration abundance with ANLN expression levels in different cancer types. These infiltrating immune cells include B cells, CD8+ T-cells, CD4+ T-cells, macrophages, neutrophils, dendritic cells (DCs), myeloid-derived suppressor cells (MDSCs), cancer-associated fibroblasts (CAFs), M2-tumor-associated macrophages (TAMs), and T-regulatory lymphocytes (Treg cells) via gene modules (Summerfield et al., 2001). Spearman’s correlation coefficients were used to analyze the correlation between ANLN and immune cell infiltration abundance. Statistical significance was set at p < 0.05. Data were visualized in GraphPad Prism Software (Version 8.0.0 for Windows). A heatmap was used to visualize the correlation of immune cell infiltration abundance with ANLN expression in 33 cancer types in the TCGA.
Analysis of the correlation between the tumor treatment response and anillin actin-binding protein expression
The Tumor Immune Dysfunction and Exclusion (TIDE) database (tide.dfci.harvard.edu) comprises omics data regarding over 33K samples obtained from 188 tumor cohorts in public databases, 998 tumors from 12 clinical studies on ICBs, and eight CRISPR screens. It provides a computational method for modeling the two primary mechanisms of tumor immune evasion (Jiang et al., 2018; Fu et al., 2020). Gene prioritization of ANLN was assessed based on the treatment response to ICBs while knockout phenotypes were identified through CRISPR screening. The z-score in Cox-PH regression was used to assess the effect of gene expression on survival in the ICB-treatment cohort. CRISPR screening was used to assess the role of knockout-mediated lymphocyte-induced tumor death in cancer models. Based on the TIDE biomarker assessment module, we examined the overall predictive utility of ANLN for treatment response outcomes and OS across cancer types; further, we compared ANLN with standardized biomarkers for tumor immune response.
The ROC Plotter (rocplot.org/) is the first online transcriptome-level validation tool for predictive biomarkers. It can link gene expression to treatment response in patients with breast, ovarian, colorectal, and glioblastoma cancer (Fekete and Gyorffy, 2019). We used the ROC plotter to assess the predictive utility of ANLN expression for treatment response in patients with cancer.
Statistical analysis
All expression data were normalized and log-transformed using log2 (data + 1). R software (V3.8.3) was used to compare differential ANLN expression between tumor and normal tissues. Between-group comparisons were performed using t-tests (ns, p ≥ 0.05; *, p < 0.05; **, p < 0.01; ***, p < 0.001). The R package [ggplot2 (version 3.3.3)] was used for visualization. The Wilcoxon rank-sum test was used to examine the relationship of ANLN expression with pathologic/cytogenetic features. For survival analysis, the hazard ratio (HR) and 95% confidence intervals (CI) were derived from Cox proportional hazards models. The p-value was determined using Kaplan-Meier analysis with the log-rank test. Statistical significance was set at p < 0.05.
Results
Topology, localization, expression profile, and protein interaction network of anillin actin-binding protein
Under physiological conditions, the topology of ANLN proteins indicated intracellular membrane localization (Figure 2A). To characterize the intracellular localization of ANLN, we used an indirect immunofluorescence assay to investigate the distribution of ANLN in the endoplasmic reticulum (ER), nucleus, and microtubules of A-431, U-2 OS, and U-251 MG cells. The Human Protein Atlas database showed that ANLN localization was mainly distributed in the nucleus, but not in the ER or microtubules, of these cells (Figure 2B). In addition, ANLN mRNA expressions occur in various normal human tissues, including the immune, visceral, nervous, endocrine, muscle, and reproductive tissues (Figure 2C). Based on the fluorescent ubiquitination-based cell cycle indicator system, single-cell RNA sequencing analysis of U-2 OS cells revealed an association of increased ANLN RNA expression with cell cycle progression (Figure 2D). Moreover, ANLN expression was correlated with the cell cycle progression in the G1, S, and S/G2 phases (Figure 2D). PPI analysis revealed that proteins related to ANLN included RHOA, RACGAP1, KIF23, ECT2, KIF11, KIF20A, CD2AP, DIAPH3, DLGAP5, and BUB1 (Figure 2E).
FIGURE 2. Topology, localization, expression profile, and protein interaction network of ANLN proteins. (A) Topology of ANLN proteins showing intracellular membrane localization. (B) Immunofluorescence staining of the subcellular distribution of ANLN in the nucleus, endoplasmic reticulum, and microtubules of A-431, U-2 OS, and U-251 MG was adopted from the HPA database. (C) Bar graphs of ANLN mRNA expression in various normal human tissues from the TCGA_GTEx database. (D) The plot of single-cell RNA sequencing data from the U-2 OS cell line shows a correlation between ANLN mRNA expression and cell cycle progression. (E) Protein interactions network of ANLN.
Anillin actin-binding protein expression in pancytopenia
Based on the TCGA_GTEx database, there was significantly higher ANLN expression in tumor tissues than in adjacent normal tissues in adrenocortical carcinoma (ACC), bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), cervical squamous cell carcinoma (CESC), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), lymphoid neoplasm diffuse large B-cell lymphoma (DLBC), esophageal carcinoma (ESCA), head and neck squamous cell carcinoma (HNSC), kidney chromophobe (KICH), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), brain lower grade glioma (LGG), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), ovarian serous cystadenocarcinoma (OV), pancreatic adenocarcinoma (PAAD), pheochromocytoma and paraganglioma (PCPG), prostate adenocarcinoma (PRAD), rectum adenocarcinoma (READ), skin cutaneous melanoma (SKCM), stomach adenocarcinoma (STAD), testicular germ cell tumors (TGCT), thyroid carcinoma (THCA), thymoma (THYM), uterine corpus endometrial carcinoma (UCEC), and uterine carcinosarcoma (UCS) (Figure 3A). Further analysis using the “CPTAC Analysis” confirmed ANLN overexpression in BRCA, renal cell carcinoma, COAD, HNSC, LIHC, LUAD, OV, PAAD, and UCEC tumors compared with in normal tissues; however, there was low ANLN expression in the glioblastoma multiforme (Figure 3B).
FIGURE 3. Upregulation of ANLN mRNA expression in pan-cancer. (A) Expression levels of ANLN in TCGA_GTEx database in different cancer types. Log2 (TPM + 1) was used for log scale. (B) Expression levels of ANLN protein in normal and different tumor tissues from the CPTAC database. *p < 0.05, **p < 0.01, and ***p < 0.001.
Correlation of anillin actin-binding protein with pathological stage and prognosis of pan-cancer
To investigate the relationship between ANLN expression and the clinicopathological features of pan-cancer, we evaluated ANLN expression in patients with pathological stages I, II, III, and IV cancers. The TCGA database showed that ANLN expression was significantly upregulated in ACC, BLCA, KIRC, LIHC, and LUAD with advanced pathological stages (Figure 4A). In pan-cancer prognostic analysis, Kaplan-Meier survival curves revealed a significant association of upregulated ANLN expression with poor OS in ACC, BLCA, BRCA, KIRC, KIRP, LIHC, LUAD, mesothelioma (MESO), and PAAD; contrastingly, low ANLN expression was associated with poor OS in THYM (Figure 4B).
FIGURE 4. Correlation of ANLN expression with pathological stage and prognosis of tumors. (A) Correlation of ANLN expression with pathological stages of ACC, BLCA, KIRC, LIHC, and LUAD; (B) Correlation of ANLN expression with a prognosis of patients with ACC, BLCA, BRCA, KIRC, KIRP, LIHC, LUAD, MESO, PAAD, and THYM. *p < 0.05, **p < 0.01, and ***p < 0.001.
DNA methylation and mutational features of anillin actin-binding protein in pan-cancer
DNA methylation directly affects cancer development and progression. We used the MethSurv to investigate DNA methylation of ANLN in terms of the prognostic value of each CpG site in cancer. Prognosis was correlated with the methylation level of four CpG sites, i.e., cg00267323, cg04455682, cg04897631, and cg12284836 (p < 0.05) (Table 1). Patients with low ANLN methylation of CpG site cg12284836 had a worse OS than patients with high ANLN methylation in this site for CESC, COAD, KIRC, LGG, LIHC, and UCEC (Figure 5).
TABLE 1. Effect of hypermethylation level on prognosis in cancer.
FIGURE 5. The Kaplan-Meier survival of the promoter methylation of ANLN.
The cBioPortal (TCGA, Pan-Cancer Atlas) database showed that lung cancer had the highest frequency of ANLN mutation (≈ 23.68%) (Figure 6A). As shown in Figure 6A, ANLN amplification and mutation were the most common alteration type. Moreover, among the patients with tumors, patients with ANLN mutations had worse progression-free survival (PFS); however, there were no significant differences in OS, disease-free survival, and disease-specific survival (DSS) (Figure 6B).
FIGURE 6. The mutational characteristics of ANLN in cancer. (A) The frequency of alterations in different types of mutations was examined using the cBioPortal database. (B) The effect of ANLN mutation status on overall survival (OS), disease-specific survival (DSS), disease-free survival (DFS), and progression-free survival (PFS) of cancer patients was investigated using the cBioPortal database.
Correlation of anillin actin-binding protein with tumor immunity
Based on the TCGA database, ANLN expression in COAD, KIRC, LIHC, and PRAD showed a significant positive correlation with infiltration of six immune cells (B cells, CD8+ T-cells, CD4+ T-cells, macrophages, neutrophils, and DCs) (p < 0.05; Figure 7A). Additionally, we assessed the correlation of ANLN expression with infiltration of four immunosuppressive cells (MDSCs, CAFs, M2-TAMs, and Treg cells), which promote T-cell rejection. ANLN expression was positively correlated with myeloid-derived suppressor cell infiltration in ACC, BLCA, BRCA, BRCA-Basal, BRCA-Her2, BRCA-LumA, BRCA-Lumb, CESC, CHOL, COAD, ESCA, HNSC, HNSC-HPV−, HNSC-HPV+, KICH, KIRC, KIRP, LIHC, LUAD, LUSC, MESO, OV, PAAD, PCPG, PRAD, READ, SARC, SKCM, SKCM-Metastasis, SKCM-Primary, STAD, THYM, UCEC, and uveal melanoma (p < 0.05). ANLN expression was positively correlated with Tregs infiltration in KICH, LIHC, PCPG, and THCA (p < 0.05); ANLN expression was positively correlated with CAF infiltration in ACC, BLCA, BRCA-LumA, ESCA, HNSC, HNSC-HPV-, KICH, KIRC, KIRP, LIHC, LUAD, LUSC, MESO, PAAD, PCPG, PRAD, SKCM, SKCM-Metastasis, THCA, UCEC, and UCS (p < 0.05). ANLN expression was positively correlated with M2-TAM infiltration in KIRC and PRAD (p < 0.05) (Figure 7B).
FIGURE 7. Correlation of ANLN with tumor immune cell infiltration. (A,B) Heat map showing a correlation between ANLN protein and infiltration of six types of immune cells and four types of immunosuppressive cells in multiple tumors; (C) Comparison of ANLN with other biomarkers for efficacy prediction of immunodetectable point-of-care inhibitors.
We compared the predictive utility of ANLN and other standardized biomarkers for the efficacy of ICB. Among the ICB subgroups, 10 had ANLN-AUC > 0.5 (Figure 7C). ANLN showed a higher predictive value than TMB, T.Clonality, and B.Clonality.
Anillin actin-binding protein is strongly correlated with tumor treatment outcomes
We evaluated the effect of ANLN expression on chemotherapy response in a clinical cancer cohort. Patients with breast and colorectal cancer who presented high ANLN expression showed resistance to chemotherapy. Contrastingly, among patients with glioblastoma and ovarian cancer, patients with high ANLN expression showed greater benefits from chemotherapy than those with low ANLN expression (Figure 8A). Additionally, high ANLN expression was associated with clinical benefits derived from ICB therapy (PD-1), and thus longer OS, in patients with bladder cancer and glioblastoma. Contrastingly, in patients with kidney cancer and melanoma, low ANLN expression was associated with clinical benefits derived from ICB therapy (PD-L1 or PD-1), and thus longer OS (Figure 8B). In patients with bladder cancer and melanoma, high ANLN expression was positively correlated with CTL levels (Figure 8B).
FIGURE 8. ANLN expression correlates with treatment response in multiple cancer types. (A) ROC plot of the association between ANLN expression and chemotherapy response in breast, brain, colorectal, and ovarian cancer cohorts. (B) Kaplan-Meier curves (upper panel) as a measure of immunotherapy response between cancer cohorts with high and low ANLN expression. The graph below shows the correlation between ANLN expression and cytotoxic T-cell levels.
Discussion
Recently, there have been increasing studies using Pan-Cancer Analysis of Whole Genomes to reveal genetic mutations, RNA changes, and oncogenes associated with cancer development and progression. Additionally, pan-cancer analysis has revealed common aspects among cancer types regarding key driver gene mutations, tumor-associated signaling pathway activation, immune escape, and microbial features. This suggests that the same therapeutic strategy can be applied for common features in different cancer types (Chen et al., 2021). Therefore, Pan-Cancer analysis may provide novel biomarkers for early diagnosis, prognosis, and treatment efficacy in patients with cancer.
ANLN, which is a key protein in cell division, can combine with GTP-Rho, actin, and cytokinin. ANLN is crucially involved in cell proliferation, differentiation, adhesion, migration, apoptosis, and cycle progression (Naydenov et al., 2021). Abnormal ANLN expression causes dysregulation of cell division and is crucially involved in cancer development, cancer cell proliferation, and metastasis. We found that ANLN mRNA expression was positively correlated with cell cycle progression; moreover, most tumors showed high ANLN expression, which was significantly correlated with poor prognosis and tumor progression. In bladder cancer, there was ANLN overexpression, which could promote tumor cell proliferation, migration, and infiltration; moreover, this was confirmed in vitro experiments. Nuclear ANLN overexpression is significantly correlated with tumor stage, histological grade, and tumor invasion; additionally, it is an independent risk factor for DSS (Liang et al., 2015; Zeng et al., 2017). In gastric cancer, ANLN can promote cell proliferation by activating the Wnt/β-catenin signaling pathway (Pandi et al., 2014). Li et al. (2020c) reported that the CDK1-PLK1-SGOL2-ANLN signaling axis mediating abnormal cell cycle division may be crucially involved in the development of hepatocellular carcinoma (HCC). Moreover, tumor tissues present with high ANLN expression and decreased ANLN methylation in HCC. Further, ANLN silencing may inhibit HCC cell proliferation in the G2/M phase as well as their migration and invasion abilities (Zhang et al., 2021). Contrastingly, the metastatic foci of HCC showed significantly higher cytosolic ANLN expression than the primary tumor tissue; additionally, high cytosolic ANLN expression was an independent risk factor for 5-year OS after hepatic resection. This study suggested that ANLN could be a potential therapeutic target and prognostic biological marker for HCC (Zhang et al., 2021). In gastric cancer, decreased ANLN expression affects downstream DNA damage repair-related pathways, resulting in more irreparable DNA damage in S-phase cells and eventually apoptosis (Wei et al., 2020). Patients with non-small cell lung cancer present with the highest ANLN expression. Inhibiting ANLN protein expression through siRNA could inhibit tumor cell growth and migration as well as induce tumor cell apoptosis (Suzuki et al., 2005). In pancreatic cancer cells, ANLN deficiency promotes miR218-5p expression, which results in apoptosis (Liu et al., 2015). Additionally, ANLN is a key modulator of inducing PI3K/Akt signaling, which suggests that ANLN is indirectly associated with apoptosis induction (Dai et al., 2019). In HNSC, ANLN is overexpressed in tumor tissues, with ANLN knockdown affecting tumor cell proliferation, migration, and invasion. Moreover, exosome-derived ANLN promotes cell differentiation and macrophage polarization via the PTEN/PI3K/Akt signaling pathway, which subsequently stimulates tumor growth. Taken together, ANLN is an independent prognostic factor for patients with HNSC (Guo et al., 2021). Our study used bioinformatics methods to further demonstrate the correlation of ANLN with tumor progression and prognosis, which provides a theoretical basis for the use of ANLN as a potential novel target for biotherapy and a biomarker for tumor surveillance.
The TIME facilitates tumor growth and function by providing a relatively stable living environment for cancer cells to continuously proliferate without interference, which enhances the tumor malignancy potential. Additionally, it contains numerous immune and immunosuppressive cells, including M2-TAMs, CAFs, Tregs, and MDSCs, which can suppress the production, activation, and infiltration of cytotoxic T-cells as well as promote tumor immune evasion, tumor growth, metastasis, and treatment resistance (Liu and Cao, 2016; Monteran and Erez, 2019). The TIME and tumor immune escape are associated with cancer prognosis and treatment (Wu and Dai, 2017). Regarding tumor immune escape mechanisms, tumor immune cell infiltration causes T-cell dysfunction, which promotes tumor evasion of host immune system surveillance as well as facilitates tumor progression, invasion, metastasis, and treatment resistance (Yu et al., 2009; Liu et al., 2022). Another tumor immune escape mechanism is T-cell rejection. Here, CAFs, Tregs, M2-TAMs, and MDSCs in the tumor microenvironment, which act as immunosuppressive cells, block immune cell infiltration (Man et al., 2013; Komohara et al., 2016). Based on the TCGA database, we analyzed the relationship of ANLN with tumor immune cell infiltration. ANLN showed a significant positive correlation with tumor immune cell infiltration only in COAD, KIRC, LIHC, and PRAD, indicating that ANLN could enhance tumor immune infiltration through T-cell dysfunction in these cancer types. Furthermore, ANLN expression was strongly correlated with the levels of immunosuppressive cells in almost all cancer types. This suggests that T-cell rejection is the main mechanism through which ANLN regulates immune cell tumor escape, tumor promotion, and metastasis.
DNA methylation modifications affect a range of biological processes such as eukaryotic cell growth, differentiation, and transformation mechanisms by regulating gene expression. DNA methylation dysregulation causes multiple diseases (Wagner, 2022) and is considered a key event in cancer development and progression, including altered DNA methylation and dysregulation of gene expression in tumor-related genes (Wang J. et al., 2021). The expression of aberrantly methylated genes is associated with the degree of tumor differentiation, TNM stage, and poor prognosis (Skvortsova et al., 2019). Therefore, we investigated the relationship of ANLN methylation levels with prognosis in patients with cancer. We found that prognosis was correlated with the methylation level of four CpG sites; namely, cg00267323, cg04455682, cg04897631, and cg12284836.
Although genetic alterations in tumors are common, premalignant genetic changes are more likely to initiate and promote cancer development. The accumulation of genetic alterations drives the progression of normal cells through hyperplastic and dysplastic stages to invasive cancer and, ultimately, metastatic disease (Garnis et al., 2004). Therefore, analyzing genetic alterations in known oncogenes could further elucidate their role in cancer progression. In our study, lung cancer showed the highest frequency of ANLN alterations (≈ 23.68%); moreover, amplification and mutation were the most common genetic alterations. In addition, patients with altered ANLN genes showed worse PFS. Nonetheless, further experimental studies are warranted to reveal the mechanisms through which ANLN gene alterations are involved in tumor development.
There has been increasing research interest in immune replacement therapy regimens based on PD-1 and PD-L1 inhibitors. PD-1 and PD-L1-based immunotherapy can effectively treat various cancers and improve prognosis (Makuku et al., 2021). We found that high ANLN expression was associated with the clinical benefits of ICB therapy (PD-1) in patients with bladder cancer and glioblastoma, while low ANLN expression was associated with the clinical benefits of ICB therapy (PD-L1 or PD-1) in patients with kidney cancer and melanoma. This suggested that ANLN is crucially involved in tumor-mediated immune escape. Further, regarding the relevance of ANLN in oncologic chemotherapy, in patients with breast and colorectal cancers, ANLN expression was negatively correlated with sensitivity to chemotherapeutic drugs. Patients with breast cancer showed high ANLN expression in tumor tissues; moreover, ANLN downregulation affected the tumor cell cycle and promoted apoptosis while inhibiting tumor cell proliferation, migration, and invasion. Therapy targeting ANLN could provide a similar therapeutic effect as chemotherapy on breast cancer (Wang Z. et al., 2021). However, the mechanism through which ANLN mediates resistance to chemotherapeutic agents should be further investigated.
Conclusion
In summary, our findings showed that ANLN is strongly associated with the prognosis, immune cell infiltration, gene mutations, and tumor treatment in patients with cancer. These findings provide further insight into the role of ANLN in tumorigenesis from an overall perspective. Furthermore, ANLN could be an onco-immunological biomarker and could serve as a hallmark for tumor screening, prognosis, individualized treatment design, and follow-up.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
KJL conceived and designed the experiments, performed the experiments, authored or reviewed drafts of the paper, and approved the final draft. LC and CQL conceived and designed the experiments, performed the experiments, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft. CFT, JH, and YMN analyzed the data, prepared figures and/or tables, and approved the final draft. YB and BDC conceived and designed the experiments, authored or reviewed drafts of the paper, and approved the final draft.
Funding
This study was funded by the Newly recruited master’s training project of the General Hospital of Ningxia Medical University (No. YKDZY2022013) and First-Class Discipline Construction Founded Proiect of NingXia Medical University and the School of Clinical Medicine (No. 2020016).
Acknowledgments
We thank PhD. Yongxue LV for the technological assistance. We would like to thank Editage (www.editage.cn) for English language editing.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.922472/full#supplementary-material
Supplementary Figure S1 | The density plots were used to visualize cut-off points for the dichotomy of methylation levels in patients.
References
Barrueto, L., Caminero, F., Cash, L., Makris, C., Lamichhane, P., and Deshmukh, R. R. (2020). Resistance to checkpoint inhibition in cancer immunotherapy. Transl. Oncol. 13 (3), 100738. doi:10.1016/j.tranon.2019.12.010
Blum, A., Wang, P., and Zenklusen, J. C. (2018). SnapShot: TCGA-analyzed tumors. Cell 173 (2), 530. doi:10.1016/j.cell.2018.03.059
Carithers, L. J., and Moore, H. M. (2015). The Genotype-Tissue expression (GTEx) project. Biopreserv. Biobank. 13 (5), 307–308. doi:10.1089/bio.2015.29031.hmm
Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2 (5), 401–404. doi:10.1158/2159-8290.CD-12-0095
Chandrashekar, D. S., Bashel, B., Balasubramanya, S., Creighton, C. J., Ponce-Rodriguez, I., Chakravarthi, B., et al. (2017). UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia 19 (8), 649–658. doi:10.1016/j.neo.2017.05.002
Chen, F., Wendl, M. C., Wyczalkowski, M. A., Bailey, M. H., Li, Y., and Ding, L. (2021). Moving pan-cancer studies from basic research toward the clinic. Nat. Cancer 2 (9), 879–890. doi:10.1038/s43018-021-00250-4
Dai, X., Chen, X., Hakizimana, O., and Mei, Y. (2019). Genetic interactions between ANLN and KDR are prognostic for breast cancer survival. Oncol. Rep. 42 (6), 2255–2266. doi:10.3892/or.2019.7332
Deng, F., Xu, Z., Zhou, J., Zhang, R., and Gong, X. (2021). ANLN regulated by miR-30a-5p mediates malignant progression of lung adenocarcinoma. Comput. Math. Methods Med. 2021, 9549287. doi:10.1155/2021/9549287
Downs, B., Sherman, S., Cui, J., Kim, Y. C., Snyder, C., Christensen, M., et al. (2019). Common genetic variants contribute to incomplete penetrance: evidence from cancer-free BRCA1 mutation carriers. Eur. J. Cancer 107, 68–78. doi:10.1016/j.ejca.2018.10.022
Fekete, J. T., and Gyorffy, B. (2019). ROCplot.org: validating predictive biomarkers of chemotherapy/hormonal therapy/anti-HER2 therapy using transcriptomic data of 3, 104 breast cancer patients. Int. J. Cancer 145 (11), 3140–3151. doi:10.1002/ijc.32369
Fu, J., Li, K., Zhang, W., Wan, C., Zhang, J., Jiang, P., et al. (2020). Large-scale public data reuse to model immunotherapy response and resistance. Genome Med. 12 (1), 21. doi:10.1186/s13073-020-0721-z
Garnis, C., Buys, T. P., and Lam, W. L. (2004). Genetic alteration and gene expression modulation during cancer progression. Mol. Cancer 3, 9. doi:10.1186/1476-4598-3-9
Genova, C., Dellepiane, C., Carrega, P., Sommariva, S., Ferlazzo, G., Pronzato, P., et al. (2022). Therapeutic implications of tumor microenvironment in lung cancer: focus on immune checkpoint blockade. Front. Immunol. 12, 799455. doi:10.3389/fimmu.2021.799455
Guo, E., Mao, X., Wang, X., Guo, L., An, C., Zhang, C., et al. (2021). Alternatively spliced ANLN isoforms synergistically contribute to the progression of head and neck squamous cell carcinoma. Cell Death Dis. 12 (8), 764. doi:10.1038/s41419-021-04063-2
Hargadon, K. M., Johnson, C. E., and Williams, C. J. (2018). Immune checkpoint blockade therapy for cancer: an overview of FDA-approved immune checkpoint inhibitors. Int. Immunopharmacol. 62, 29–39. doi:10.1016/j.intimp.2018.06.001
Idichi, T., Seki, N., Kurahara, H., Yonemori, K., Osako, Y., Arai, T., et al. (2017). Regulation of actin-binding protein ANLN by antitumor miR-217 inhibits cancer cell aggressiveness in pancreatic ductal adenocarcinoma. Oncotarget 8 (32), 53180–53193. doi:10.18632/oncotarget.18261
Jiang, P., Gu, S., Pan, D., Fu, J., Sahu, A., Hu, X., et al. (2018). Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat. Med. 24 (10), 1550–1558. doi:10.1038/s41591-018-0136-1
Karlsson, M., Zhang, C., Mear, L., Zhong, W., Digre, A., Katona, B., et al. (2021). A single-cell type transcriptomics map of human tissues. Sci. Adv. 7 (31), eabh2169. doi:10.1126/sciadv.abh2169
Komohara, Y., Fujiwara, Y., Ohnishi, K., and Takeya, M. (2016). Tumor-associated macrophages: potential therapeutic targets for anti-cancer therapy. Adv. Drug Deliv. Rev. 99 (Pt B), 180–185. doi:10.1016/j.addr.2015.11.009
Lanczky, A., and Gyorffy, B. (2021). Web-Based survival analysis tool tailored for medical research (KMplot): development and implementation. J. Med. Internet Res. 23 (7), e27633. doi:10.2196/27633
Li, L., Huang, K., Zhao, H., Chen, B., Ye, Q., and Yue, J. (2020c). CDK1-PLK1/SGOL2/ANLN pathway mediating abnormal cell division in cell cycle may be a critical process in hepatocellular carcinoma. Cell Cycle 19 (10), 1236–1252. doi:10.1080/15384101.2020.1749471
Li, T., Fu, J., Zeng, Z., Cohen, D., Li, J., Chen, Q., et al. (2020b). TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 48 (W1), W509–W514. doi:10.1093/nar/gkaa407
Li, Y., Liu, J., Gao, L., Liu, Y., Meng, F., Li, X., et al. (2020a). Targeting the tumor microenvironment to overcome immune checkpoint blockade therapy resistance. Immunol. Lett. 220, 88–96. doi:10.1016/j.imlet.2019.03.006
Lian, Y., Huang, Y., Wang, J., Deng, M., Xia, T., Zeng, M., et al. (2018). Anillin is required for tumor growth and regulated by miR-15a/miR-16-1 in HBV-related hepatocellular carcinoma. Aging 10 (8), 1884–1901. doi:10.18632/aging.101510
Liang, P. I., Chen, W. T., Li, C. F., Li, C. C., Li, W. M., Huang, C. N., et al. (2015). Subcellular localisation of anillin is associated with different survival outcomes in upper urinary tract urothelial carcinoma. J. Clin. Pathol. 68 (12), 1026–1032. doi:10.1136/jclinpath-2015-202958
Liu, X., Hoft, D. F., and Peng, G. (2022). Tumor microenvironment metabolites directing T cell differentiation and function. Trends Immunol. 43 (2), 132–147. doi:10.1016/j.it.2021.12.004
Liu, Y., and Cao, X. (2016). Immunosuppressive cells in tumor immune escape and metastasis. J. Mol. Med. 94 (5), 509–522. doi:10.1007/s00109-015-1376-x
Liu, Z., Xu, Y., Long, J., Guo, K., Ge, C., and Du, R. (2015). MicroRNA-218 suppresses the proliferation, invasion and promotes apoptosis of pancreatic cancer cells by targeting HMGB1. Chin. J. Cancer Res. 27 (3), 247–257. doi:10.3978/j.issn.1000-9604.2015.04.07
Makuku, R., Khalili, N., Razi, S., Keshavarz-Fathi, M., and Rezaei, N. (2021). Current and future perspectives of PD-1/PDL-1 blockade in cancer immunotherapy. J. Immunol. Res. 2021, 6661406. doi:10.1155/2021/6661406
Man, Y. G., Stojadinovic, A., Mason, J., Avital, I., Bilchik, A., Bruecher, B., et al. (2013). Tumor-infiltrating immune cells promoting tumor invasion and metastasis: existing theories. J. Cancer 4 (1), 84–95. doi:10.7150/jca.5482
Mao, X., Xu, J., Wang, W., Liang, C., Hua, J., Liu, J., et al. (2021). Crosstalk between cancer-associated fibroblasts and immune cells in the tumor microenvironment: new findings and future perspectives. Mol. Cancer 20 (1), 131. doi:10.1186/s12943-021-01428-1
Modhukur, V., Iljasenko, T., Metsalu, T., Lokk, K., Laisk-Podar, T., and Vilo, J. (2018). MethSurv: a web tool to perform multivariable survival analysis using DNA methylation data. Epigenomics 10 (3), 277–288. doi:10.2217/epi-2017-0118
Monteran, L., and Erez, N. (2019). The dark side of fibroblasts: cancer-associated fibroblasts as mediators of immunosuppression in the tumor microenvironment. Front. Immunol. 10, 1835. doi:10.3389/fimmu.2019.01835
Morad, G., Helmink, B. A., Sharma, P., and Wargo, J. A. (2021). Hallmarks of response, resistance, and toxicity to immune checkpoint blockade. Cell 184 (21), 5309–5337. doi:10.1016/j.cell.2021.09.020
Naydenov, N. G., Koblinski, J. E., and Ivanov, A. I. (2021). Anillin is an emerging regulator of tumorigenesis, acting as a cortical cytoskeletal scaffold and a nuclear modulator of cancer cell differentiation. Cell. Mol. Life Sci. 78 (2), 621–633. doi:10.1007/s00018-020-03605-9
Omasits, U., Ahrens, C. H., Muller, S., and Wollscheid, B. (2014). Protter: interactive protein feature visualization and integration with experimental proteomic data. Bioinformatics 30 (6), 884–886. doi:10.1093/bioinformatics/btt607
Pandi, N. S., Manimuthu, M., Harunipriya, P., Murugesan, M., Asha, G. V., and Rajendran, S. (2014). In silico analysis of expression pattern of a Wnt/beta-catenin responsive gene ANLN in gastric cancer. Gene 545 (1), 23–29. doi:10.1016/j.gene.2014.05.013
Siegel, R. L., Miller, K. D., Fuchs, H. E., and Jemal, A. (2021). Cancer statistics, 2021. CA. Cancer J. Clin. 71 (1), 7–33. doi:10.3322/caac.21654
Skvortsova, K., Stirzaker, C., and Taberlay, P. (2019). The DNA methylation landscape in cancer. Essays Biochem. 63 (6), 797–811. doi:10.1042/EBC20190037
Summerfield, A., McNeilly, F., Walker, I., Allan, G., Knoetig, S. M., and McCullough, K. C. (2001). Depletion of CD4(+) and CD8(high+) T-cells before the onset of viraemia during classical swine fever. Vet. Immunol. Immunopathol. 78 (1), 3–19. doi:10.1016/s0165-2427(00)00248-8
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA. Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Suzuki, C., Daigo, Y., Ishikawa, N., Kato, T., Hayama, S., Ito, T., et al. (2005). ANLN plays a critical role in human lung carcinogenesis through the activation of RHOA and by involvement in the phosphoinositide 3-kinase/AKT pathway. Cancer Res. 65 (24), 11314–11325. doi:10.1158/0008-5472.CAN-05-1507
Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. (2021). The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 49 (D1), D605–D612. doi:10.1093/nar/gkaa1074
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 45 (W1), W98–W102. doi:10.1093/nar/gkx247
Thul, P. J., Akesson, L., Wiking, M., Mahdessian, D., Geladaki, A., Ait, B. H., et al. (2017). A subcellular map of the human proteome. Science 356 (6340), eaal3321. doi:10.1126/science.aal3321
Uhlen, M., Fagerberg, L., Hallstrom, B. M., Lindskog, C., Oksvold, P., Mardinoglu, A., et al. (2015). Proteomics. Tissue-based map of the human proteome. Science 347 (6220), 1260419. doi:10.1126/science.1260419
Uhlen, M., Zhang, C., Lee, S., Sjostedt, E., Fagerberg, L., Bidkhori, G., et al. (2017). A pathology atlas of the human cancer transcriptome. Science 357 (6352), eaan2507. doi:10.1126/science.aan2507
Wagner, W. (2022). How to translate DNA methylation biomarkers into clinical practice. Front. Cell Dev. Biol. 10, 854797. doi:10.3389/fcell.2022.854797
Wang, A., Dai, H., Gong, Y., Zhang, C., Shu, J., Luo, Y., et al. (2019). ANLN-induced EZH2 upregulation promotes pancreatic cancer progression by mediating miR-218-5p/LASP1 signaling axis. J. Exp. Clin. Cancer Res. 38 (1), 347. doi:10.1186/s13046-019-1340-7
Wang, G., Shen, W., Cui, L., Chen, W., Hu, X., and Fu, J. (2016). Overexpression of Anillin (ANLN) is correlated with colorectal cancer progression and poor prognosis. Cancer Biomark. 16 (3), 459–465. doi:10.3233/CBM-160585
Wang, J., Yang, J., Li, D., and Li, J. (2021). Technologies for targeting DNA methylation modifications: basic mechanism and potential application in cancer. Biochim. Biophys. Acta. Rev. Cancer 1875 (1), 188454. doi:10.1016/j.bbcan.2020.188454
Wang, Z., Hu, S., Li, X., Liu, Z., Han, D., Wang, Y., et al. (2021). MiR-16-5p suppresses breast cancer proliferation by targeting ANLN. BMC Cancer 21 (1), 1188. doi:10.1186/s12885-021-08914-1
Wei, S., Teng, S., Yao, J., Gao, W., Zang, J., Wang, G., et al. (2020). Develop a circular RNA–related regulatory network associated with prognosis of gastric cancer. Cancer Med. 9 (22), 8589–8599. doi:10.1002/cam4.3035
Wu, T., and Dai, Y. (2017). Tumor microenvironment and therapeutic response. Cancer Lett. 387, 61–68. doi:10.1016/j.canlet.2016.01.043
Xiao, J. X., Xu, W., Fei, X., Hao, F., Wang, N., Chen, Y., et al. (2020). Anillin facilitates cell proliferation and induces tumor growth of hepatocellular carcinoma via miR-138/SOX4 axis regulation. Transl. Oncol. 13 (10), 100815. doi:10.1016/j.tranon.2020.100815
Yu, G. P., Chiang, D., Song, S. J., Hoyte, E. G., Huang, J., Vanishsarn, C., et al. (2009). Regulatory T cell dysfunction in subjects with common variable immunodeficiency complicated by autoimmune disease. Clin. Immunol. 131 (2), 240–253. doi:10.1016/j.clim.2008.12.006
Zeng, S., Yu, X., Ma, C., Song, R., Zhang, Z., Zi, X., et al. (2017). Transcriptome sequencing identifies ANLN as a promising prognostic biomarker in bladder urothelial carcinoma. Sci. Rep. 7 (1), 3151. doi:10.1038/s41598-017-02990-9
Keywords: pan-cancer, anillin actin-binding protein (ANLN), immunization, prognosis, bioinformatics analysis
Citation: Liu K, Cui L, Li C, Tang C, Niu Y, Hao J, Bu Y and Chen B (2022) Pan-cancer analysis of the prognostic and immunological role of ANLN: An onco-immunological biomarker. Front. Genet. 13:922472. doi: 10.3389/fgene.2022.922472
Received: 18 April 2022; Accepted: 05 July 2022;
Published: 05 August 2022.
Edited by:
Ehsan Nazemalhosseini-Mojarad, Shahid Beheshti University of Medical Sciences, IranReviewed by:
Mohammad Taheri-Ledari, University of Tehran, IranVahid Chaleshi, Shahid Beheshti University of Medical Sciences, Iran
Copyright © 2022 Liu, Cui, Li, Tang, Niu, Hao, Bu and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yang Bu, boyang1976@163.com; Bendong Chen, bendong8511@hotmail.com
†These authors have contributed equally to this work