Next Article in Journal
Clinical Presentation and Risk Factors of Infective Endocarditis in the Elderly: A Systematic Review
Previous Article in Journal
Clinical Applications of Three-Dimensional Printing in Upper Extremity Surgery: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

ANO4 Expression Is a Potential Prognostic Biomarker in Non-Metastasized Clear Cell Renal Cell Carcinoma

1
Faculty of Medicine, Jordan University of Science and Technology, Irbid 22110, Jordan
2
Faculty of Pharmacy, Jordan University of Science and Technology, Irbid 22110, Jordan
3
Faculty of Medicine, Yarmouk University, Irbid 21163, Jordan
4
College of Pharmacy, QU Health, Qatar University, Doha 2713, Qatar
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2023, 13(2), 295; https://doi.org/10.3390/jpm13020295
Submission received: 2 January 2023 / Revised: 28 January 2023 / Accepted: 2 February 2023 / Published: 7 February 2023

Abstract

:
Background: Over the past decade, transcriptome profiling has elucidated many pivotal pathways involved in oncogenesis. However, a detailed comprehensive map of tumorigenesis remains an enigma to solve. Propelled research has been devoted to investigating the molecular drivers of clear cell renal cell carcinoma (ccRCC). To add another piece to the puzzle, we evaluated the role of anoctamin 4 (ANO4) expression as a potential prognostic biomarker in non-metastasized ccRCC. Methods: A total of 422 ccRCC patients with the corresponding ANO4 expression and clinicopathological data were obtained from The Cancer Genome Atlas Program (TCGA). Differential expression across several clinicopathological variables was performed. The Kaplan–Meier method was used to assess the impact of ANO4 expression on the overall survival (OS), progression-free interval (PFI), disease-free interval (DFI), and disease-specific survival (DSS). Univariate and multivariate Cox logistic regression analyses were conducted to identify independent factors modulating the aforementioned outcomes. Gene set enrichment analysis (GSEA) was used to discern a set of molecular mechanisms involved in the prognostic signature. Tumor immune microenvironment was estimated using xCell. Results: ANO4 expression was upregulated in tumor samples compared to normal kidney tissue. Albeit the latter finding, low ANO4 expression is associated with advanced clinicopathological variables such as tumor grade, stage, and pT. In addition, low ANO4 expression is linked to shorter OS, PFI, and DSS. Multivariate Cox logistic regression analysis identified ANO4 expression as an independent prognostic variable in OS (HR: 1.686, 95% CI: 1.120–2.540, p = 0.012), PFI (HR: 1.727, 95% CI: 1.103–2.704, p = 0.017), and DSS (HR: 2.688, 95% CI: 1.465–4.934, p = 0.001). GSEA identified the following pathways to be enriched within the low ANO4 expression group: epithelial–mesenchymal transition, G2-M checkpoint, E2F targets, estrogen response, apical junction, glycolysis, hypoxia, coagulation, KRAS, complement, p53, myogenesis, and TNF-α signaling via NF-κB pathways. ANO4 expression correlates significantly with monocyte (ρ = −0.1429, p = 0.0033) and mast cell (ρ = 0.1598, p = 0.001) infiltration. Conclusions: In the presented work, low ANO4 expression is portrayed as a potential poor prognostic factor in non-metastasized ccRCC. Further experimental studies should be directed to shed new light on the exact molecular mechanisms involved.

1. Introduction

Renal cell carcinoma (RCC) is the most common genitourinary tumor in the United States, with a projected 79,000 new cases and 13,920 deaths in the year 2022 [1]. There is a twofold increase in the likelihood of men being diagnosed with kidney cancer as compared to women [1]. About 60% of cases are detected incidentally and 20% to 30% of patients present with metastatic disease [2]. In most cases, metastatic renal cancer is incurable [2,3]. Rather than being a single disease, RCC refers to a group of heterogeneous tumors that arise from the renal epithelium [4,5]. Each tumor subtype exhibits a unique histology, clinical course, and genetic profile, as well as a different therapeutic response [4,5,6]. Three major histological subtypes of RCC were identified, namely clear cell RCC or kidney renal clear cell carcinoma (ccRCC or KIRC), papillary RCC or kidney renal papillary cell carcinoma (pRCC or KIRP), and chromophobe RCC or kidney chromophobe (chrRCC or KICH) [7,8,9]. The most prevalent and aggressive subtype is the ccRCC, which accounts for more than 70% of all cases [4]. In addition, sarcomatoid and rhabdoid RCC tumors are associated with poor prognosis and resistance to targeted therapies [10].
ccRCC is one of the most inscrutable types of cancer [11]. An array of genetic or acquired factors can contribute to its development [11]. In terms of acquired risk factors, smoking, hypertension, obesity, chronic analgesic use, and diabetes are the most prevalent [11]. Genetically, ccRCC tumors exhibit significant mutation heterogeneity [12]. The vast majority of patients with ccRCC have a deletion in the short arm of chromosome 3 (loss of 3p) [13]. The two most common genetic abnormalities implicated in ccRCC are VHL (von Hippel–Lindau tumor suppressor) and PBRM1 (protein polybromo1) [11,12]. Other genomic alterations involve SETD2 (SET domain containing 2, histone lysine methyltransferase), KDM5C (lysine demethylase 5C), or BAP1 (BRCA1 associated protein 1) [14,15]. The primary treatment for localized ccRCC is surgical resection with no role for conventional chemotherapy and radiotherapy [2,16,17].
Despite the development of numerous drugs targeting genomically prioritized pathways, patients with ccRCC have had limited responses to these treatments [2,18]. These results demonstrate that ccRCC tumorigenesis is a complex process, and hence, a thorough genomic, epigenomic, transcriptomic, and proteogenomic analysis is still further required to fully understand this cancer type in order to discover a curative treatment [13,19]. Therefore, the objective of this study is to evaluate the role of anoctamin 4 (ANO4) expression as a potential prognostic biomarker in non-metastasized ccRCC. ANO4 belongs to the anoctamin (anion channels with 8 transmembrane domains) family [20]. It was formerly known as TMEM16 (transmembrane proteins with 16 domains) [20]. In vertebrates, it includes ten paralogues with high sequence conservation [21]. Despite their close relationship, these proteins perform different functions and show distinct tissue distributions [22]. For example, some members such as ANO1 and ANO2 act as Ca2+-dependent ion channels, whereas others including ANO6 function as Ca2+-dependent scramblases [23]. The roles played by other family members remain unclear and controversial [21]. However, several diseases have been associated with anoctamin proteins, including cancer, muscular dystrophy, asthma, arthritis, and epilepsy [20,24]. ANO4 is primarily expressed in the central nervous system and certain endocrine glands, and it has been linked to a variety of neurological disorders [21]. To date, no study has evaluated ANO4 expression as a potential prognostic biomarker in ccRCC.

2. Materials and Methods

2.1. Clinical and Transcriptomic Data Acquisition and Processing

This work aims to investigate the prognostic utility of ANO4 mRNA expression in non-metastasized ccRCC (KIRC). KIRC clinical and transcriptomic data from The Cancer Genome Atlas Program (TCGA) were accessed using the California Santa Cruz Cancer Genomics Browser (UCSC Xena, http://xena.ucsc.edu, accessed on the 1 October 2022); a web-based platform for visualizing and analyzing public genomic data resources [25]. Experimental genotypic profiling was performed using Illumina HiSeq 2000 RNA sequencing platform to obtain level 3 data. ANO4 expression from tumor samples with normal adjacent tissues was downloaded in an RNA-Seq by expectation maximization (RSEM) normalized count transformed as log 2   ( x + 1 ) . Clinicopathological variables incluedg age, gender, American Joint Committee on Cancer (AJCC) stage, International Society of Urologic Pathologists (ISUP) grade alongside the TNM scoring system which comprises tumor size, lymph node involvement, and metastasis status. Primary end points were the overall survival (OS), progression-free interval (PFI), disease-free interval (DFI), and disease-specific survival (DSS). The KIRC cohort was curated by omitting patients with an OS time of 0, metastasis status of M1 and MX, and patients without ANO4 expression data. KIRC cohort was also divided into two subsets based on ANO4 expression (high vs low expression), with a cut-point determined by X-tile software [26].

2.2. Statistical Analysis

IBM SPSS statistical package for Windows v.26 (Armonk, NY, USA) and GraphPad prism v.9.3.1 (San Diego, CA, USA) were utilized for statistical analysis and graph generation. Nominal data were presented as frequency (percentage). Mean ± standard deviation of the mean (SD) or standard error of the mean (SEM) were used to present normally distributed continuous variables, while non-normally distributed data were presented as median (interquartile range (IQR)). Normality was evaluated using the Kolmogorov–Smirnov test and the Shapiro–Wilk test aided with quantile–quantile (Q–Q) plots. Comparison between ANO4 expression status against clinicopathological variables was performed as follows: Chi-square test or Fisher’s exact test for categorical variables, paired t-test for normally distributed paired samples, unpaired t-test and Welch’s corrected unpaired t-test for normally distributed non-paired samples according to variance equality, and finally, Wilcoxon matched pairs test and Mann–Whitney U-test for non-normally distributed data.
Kaplan–Meier survival methods were utilized to evaluate the impact of ANO4 expression status in relation to OS, PFI, DFI, and DSS. Statistical difference across survival curves was detected using a log-rank test reporting the p-value, 95% confidence intervals (95% CI), and hazard ratios (HR). Univariate and multivariate Cox logistic regression models were applied to test clinicopathological variables and ANO4 expression as independent prognostic indicators. Variables were dichotomized for Cox logistic regression analysis as follows: age (ref. ≤ 53 years), gender (ref. Female), pT stage (ref. T1 + T2), pN stage (ref. N0), AJCC stage (ref. Stage 1 + 2), ISUP grade (ref. Grade 1 + 2), and ANO4 expression (ref. Low). The age cut-off value (53 years) was defined previously based on age-related differentially expressed genes optimized by the TCGA and Surveillance Epidemiology and End Results (SEER) database [27]. All statistical tests were two-sided and a p ≤ 0.05 was considered statistically significant.

2.3. Gene Set Enrichment Analysis (GSEA)

GSEA was performed to discern a set of molecular mechanisms involved in the prognostic ANO4 signature. GSEA performs genome-wide transcriptional profiling across two expressional groups (high vs low ANO4 expression) against a set of genes representing pivotal processes involved in oncogenesis/tumorigenesis [28]. Three molecular signatures databases were involved in the analysis, namely: hallmark gene sets, the C2 positional gene sets (Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in cancer), and C5 ontology gene sets (BP (biological process), CC (cellular component), and MF (molecular function)). The number of permutations was set to 1000 with the “gene set” permutation type. The chip platform was Human_UniProt_IDs. Gene sets were considered significantly enriched with an adjusted p value ≤ 0.05 and a false discovery rate (FDR) < 0.25.

2.4. Protein–Protein Interaction (PPI) Network Construction

Identification and retrieval of ANO4-related interacting genes were browsed using STRING v.11 (https://string-db.org, accessed on 1 October 2022), an online PPI networking resource based on functional interactions annotated using gene-enrichment analysis, GO/KEGG classification systems, and high-throughput text-mining [29]. Significant interactions were labeled if a combined score ≥ 0.4 was observed. The maximum number of interactions was limited to 50.

2.5. Immune Infiltration Analysis

Immune cell infiltration was enumerated from transcriptomes using xCell, an online platform with pre-calculated TCGA immune infiltration estimates that deploys a curve-fitting approach with a novel spillover compensation method [30]. Spearman’s correlation test was used to report Spearman’s rank correlation coefficient (rho, ρ) and the correlation significance.

3. Results

The presented work studies the usefulness of ANO4 expression as a potential prognostic biomarker in non-metastasized ccRCC using a modified TCGA-KIRC cohort. It consisted of 422 patients with a median age of 61 (IQR: 51–71) and a male predominance (275, 65.2%). The majority of patients with non-metastasized ccRCC were found to have tumor grade 2 (204, 48.3%), followed by grade 3 (162, 38.4%) and grade 4 (39, 9.2%). In regard to tumor staging, the cohort was mainly distributed across two stages including Stage 1 (243, 57.6%) and stage 3 (121, 28.7%), which was definitely in concordance with the pathological T and N scoring distribution. Table 1 represents the baseline demographical and clinicopathological characteristics of the modified TCGA-KIRC cohort. X-tile software divided the patients based on ANO4 expression harboring two groups (low expression, n = 126 and high expression, n =296).

3.1. Low ANO4 Expression Is Correlated with Poor Clinicopathological Features in Non-Metastasized ccRCC

In comparison with normal kidney tissue, the expression of ANO4 was significantly upregulated (p < 0.0001) in the non-metastasized KIRC (Figure 1A). The latter finding was confirmed by analyzing the expression of ANO4 using age and gender-matched participants (p < 0.0001), as illustrated in Figure 1B. Although tumor tissues exhibit a high ANO4 expression compared to normal tissue, intriguingly, low ANO4 expression within tumor samples is associated with advanced demographics and clinicopathological features. ANO4 expression did not exhibit a tractable significant difference across age subgroups (p = 0.1989), as seen in Figure 1C. On the other hand, samples obtained from male participants have a significant low ANO4 expression (p < 0.0001) in contrast to female samples (Figure 1D). The advanced tumor grade group, including grades 3 and 4, has significant reduced ANO4 expression (p = 0.0497) in comparison to grades 1 and 2 (Figure 1E). The same pattern tends to be observed also in the tumor stage (p = 0.004, Figure 1F) and pT (p < 0.0048, Figure 1G). Despite the previous associations, a non-significant difference across groups with and without lymph involvement (p = 0.524) was observed (Figure 1H). Table 2 summarizes the differences between low and high-expression groups based on the dichotomized demographics and clinicopathological features.

3.2. Low ANO4 Expression Is Associated with Poor OS, PFI, and DSS

Kaplan–Meier survival analysis was used to explore the impact of ANO4 expression (low vs high expression) on the OS, PFI, DFI, and DSS. The low ANO4 expression group has a poor OS (HR = 2.093, 95% CI: 1.368–3.203, p = 0.0001) compared to the high expression group (Figure 1I). Of note, the same finding was observed regarding the PFI (HR = 2.485, 95% CI: 1.544–4.001, p < 0.0001, Figure 1J) and the DSS (HR = 3.879, 95% CI: 2.068–7.278, p < 0.0001, Figure 1K). However, the DFI did not remarkably differ across test groups (HR = 1.422, 95% CI: 0.4486–4.505, p = 0.513) as seen in Figure 1L. Further analysis involving curves stratification based on the dichotomized demographics and clinicopathological features was performed. Such analysis focuses on the effect of ANO4 expression on the OS in specific sub-groups. In patients below the age of 53, ANO4 expression did not impact the OS (HR = 1.715, 95% CI: 0.5709–5.153, p = 0.268, Figure 2A). On the contrary, the low ANO4 expression group has a poor OS in the subgroup with an age above 53 (HR = 1.977, 95% CI: 1.257–3.109, p = 0.001, Figure 2B). Furthermore, gender-related OS analysis revealed a significant difference among the expression groups with the same effect of low ANO4 expressing (Male: HR = 1.758, 95% CI: 1.055–2.931, p = 0.0203) vs (Female: HR = 3.246, 95% CI: 1.421–7.416, p < 0.0001) as presented in Figure 2C,D, respectively. In both low- and high-grade groups, low ANO4 expression still exhibits the same survival effect (G1 + G2: HR = 1.946, 95% CI: 0.9375–4.038, p = 0.0378) vs (G3 + G4: HR = 1.972, 95% CI: 1.180–3.296, p = 0.0043) as presented in Figure 2E,F, respectively. While the low-stage group follows the same low ANO4 expression effect (Stage 1 + 2: HR = 2.025, 95% CI: 1.086–3.775, p = 0.0085, Figure 2G), a non-significant survival difference was observed in the high-stage group (Stage 3 + 4: HR = 1.712, 95% CI: 0.9721–3.016, p = 0.0518, Figure 2H). The same effect was also in harmony regarding the pT, that low ANO4 expression is still a poor survival indicator in both low pT (T1 + T2: HR = 2.035, 95% CI: 1.105–3.747, p = 0.0072, Figure 2I) and high pT (T3 + T4: HR = 1.752, 95% CI: 0.9844–3.117, p = 0.0444, Figure 2J). Finally, survival curves generated in sub-sets with and without lymph node involvement were impacted in the same trend by low ANO4 expression (N0: HR = 2.045, 95% CI: 1.155–3.623, p = 0.0059, Figure 2K) vs (N1 + NX: HR = 2.138, 95% CI: 1.131–4.042, p = 0.006, Figure 2L).

3.3. The Independent Prognostic Value of ANO4 Expression Evaluated by Univariate and Multivariate Cox Logistic Regression Analysis

Univariate Cox logistic regression analysis of the OS identified five independent prognostic indicators, including age (HR: 0.418, 95% CI: 0.225–0.686, p = 0.001), grade (HR: 0.559, 95% CI: 0.378–0.827, p = 0.004), stage (HR: 0.388, 95% CI: 0.266–0.566, p < 0.001), pT (HR: 0.386, 95% CI: 0.265–0.564, p < 0.001), and ANO4 expression (HR: 2.099, 95% CI: 1.434–3.073, p < 0.001), while the gender (p = 0.596) and pN (p = 0.349) were not prognostic. However, multivariate Cox logistic regression highlighted the age (HR: 0.510, 95% CI: 0.305–0.852, p = 0.010) and ANO4 expression (HR: 1.686, 95% CI: 1.120–2.540, p = 0.012) as significant predictors, while the grade (p = 0.082), stage (p = 0.617), and pT (p = 0.912) were not significant. In regard to PFI, univariate Cox logistic regression analysis illustrated a prognostic utility of all indicators except for the pN (p = 0.617) which moves in concordance with the OS results. The significant PFI predictors are age (HR: 0.588, 95% CI: 0.360–0.960, p = 0.034), gender (HR: 0.621, 95% CI: 0.388–0.992, p = 0.046), grade (HR: 0.444, 95% CI: 0.284–0.693, p < 0.001), stage (HR: 0.249, 95% CI: 0.163–0.379, p < 0.001), pT (HR: 0.264, 95% CI: 0.174–0.402, p < 0.001), and ANO4 expression (HR: 2.521, 95% CI: 1.654–3.840, p < 0.001). Upon pooling proposed predictors in PFI multivariate Cox logistic regression, stage (HR: 0.129, 95% CI: 0.092–0.579, p = 0.008) and ANO4 expression (HR: 1.727, 95% CI: 1.103–2.704, p = 0.017) were the only significant predictors. A total of six predictors were identified in the univariate Cox logistic regression analysis of the DSS. They include age (HR: 0.097, 95% CI: 0.248–0.995, p = 0.048), grade (HR: 0.328, 95% CI: 0.174–0.619, p = 0.001), stage (HR: 0.201, 95% CI: 0.112–0.359, p < 0.001), pT (HR: 0.210, 95% CI: 0.118–0.374, p < 0.001), and ANO4 expression (HR: 3.894, 95% CI: 2.201–6.889, p < 0.001). On the other hand, the DSS multivariate Cox logistic regression portrayed ANO4 expression as the only predictor (HR: 2.688, 95% CI: 1.465–4.934, p = 0.001). Tumor stage was the only predictor in the DFI-related univariate (HR: 0.311, 95% CI: 0.112–0.865, p = 0.025) and multivariate (HR: 0.258, 95% CI: 0.079–0.848, p = 0.026) Cox logistic regression analysis. Table 3 summarizes the Cox regression results discussed in this section.

3.4. ANO4-Related Signaling Pathways Based on GSEA

The first GSEA involved the hallmark gene sets in which a total of 34 out of 50 gene sets were enriched in the low ANO4 expression group, with 21 gene sets having FDR < 25%. In addition, only 14 gene sets reached the level of statistical significance (p ≤ 0.05). The enrichment results include the following sets: epithelial–mesenchymal transition (ES: 0.53, p < 0.001, FDR: < 0.001), G2-M checkpoint (ES: 0.47, p < 0.001, FDR < 0.001), estrogen response-late (ES: 0.44, p < 0.001, FDR = 0.003), E2F targets (ES: 0.45, p < 0.001, FDR = 0.003), apical junction (ES: 0.40, p < 0.001, FDR = 0.023), glycolysis (ES: 0.40, p = 0.001, FDR = 0.023), hypoxia (ES: 0.39, p = 0.001, FDR = 0.021), coagulation (ES: 0.40, p = 0.004, FDR = 0.024), KRAS signaling pathway (ES: 0.38, p = 0.003, FDR = 0.035), complement (ES: 0.37, p = 0.005, FDR = 0.044), p53 pathway (ES: 0.37, p = 0.007, FDR = 0.042), estrogen response-early (ES: 0.36, p = 0.008, FDR < 0.044), myogenesis (ES: 0.36, p = 0.008, FDR = 0.054), and TNF-α signaling pathway via NF-κB (ES: 0.35, p = 0.009, FDR = 0.061). The high ANO4 expression group enrichment analysis revealed a total of 16 enriched gene sets out of 50, with only three genes with FDR < 25% and a p ≤ 0.05. The enriched sets include fatty acid metabolism (ES: 0.39, p < 0.001, FDR = 0.043), pancreatic beta cells (ES: 0.48, p = 0.014, FDR = 0.024), and bile acid metabolism (ES: 0.39, p = 0.004, FDR = 0.037). Figure 3 displays the statistically significant enrichment plots in both the low and high ANO4 expression groups.
Figure 4A represents the top 15 KEGG-enriched pathways within low and high ANO4 expression groups. The poor prognostic effect of ANO4 downregulation is hypothesized to be mediated through the following pathways: ribosome (ES: 0.57, p < 0.001), α-linolenic acid metabolism (ES: 0.71, p = 0.00177), glycosaminoglycan biosynthesis–keratan sulfate (ES: 0.68, p = 0.010695), glycosaminoglycan biosynthesis–chondroitin sulfate (ES: 0.61, p = 0.006689), complement and coagulation cascades (ES: 0.49, p < 0.001), arachidonic acid metabolism (ES: 0.48, p = 0.018377), NOD-like receptor signaling pathway (ES: 0.47, p = 0.010671), proteasome (ES: 0.49, p = 0.024961), basal cell carcinoma (ES: 0.46, p = 0.01506), base excision repair (ES: 0.50, p = 0.024958), cell cycle (ES: 0.41, p = 0.004104), homologous recombination (ES: 0.54, p = 0.033333), p53 pathway (ES: 0.44, p = 0.023495), taste transduction (ES: 0.46, p = 0.018576), and Wnt signaling pathway (ES: 0.38, p = 0.007022). On the other hand, the high ANO4 expression group was associated with significant enrichment of 11 pathways out of the top 15 pathways, with an FDR < 0.25. The pathways are propanoate metabolism (ES: −0.65, p < 0.001), valine, leucine, and isoleucine degradation (ES: −0.60, p < 0.001), PPAR signaling pathway (ES: −0.54, p < 0.001), fatty acid metabolism (ES: −0.54, p = 0.00489), peroxisome (ES: 0.56, p = 0.005277), renin–angiotensin system pathway (ES: −0.55, p = 0.020725), pyruvate metabolism (ES: −0.54, p = 0.020725), proximal tubule bicarbonate reclamation (ES: −0.50, p = 0.028369), TCA cycle (ES: −0.41, p = 0.008772), nitrogen metabolism (ES: −0.48, p = 0.01355), and histidine metabolism (ES: −0.54, p = 0.026667). Figure 4B represents the GO results highlighting the top five enriched pathways in each subcategory, including MF, CC, and BP.

3.5. The PPI Network of ANO4 Based on String Analysis

Figure 5A represents the ANO4-related PPI. The generated PPI contains 11 nodes and 12 edges with an average node degree of 2.18, an average local clustering coefficient of 0.883, and PPI enrichment p-value of 0.329. The nodes include the following PLSCR5, CCDC181, ACOT7, SLC17A8, ASCL1, ASH2L, TSPAN16, DYSF, SPIC, and GAS2L3.

3.6. The Correlation between ANO4 Expression and Tumor Immune Infiltrate

Pre-calculated TCGA immune infiltrate estimates were obtained from xCell. ANO4 expression correlates significantly with monocyte (ρ = −0.1429, p = 0.0033) and mast cell (ρ = 0.1598, p = 0.001) infiltration. On the other hand, non-significant correlations were observed regarding neutrophils (ρ = −0.04231, p = 0.3859), macrophages (ρ = −0.0744, p = 0.1270), NK cells (ρ = 0.07838, p = 0.1079), basophils (ρ = −0.09535, p = 0.0503), B-cells (ρ = −0.02063, p = 0.6726), plasma cells (ρ = 0.06586, p = 0.1769), CD4+ T-cells (ρ = 0.02585, p = 0.5964), and CD8+ T-cells (ρ = 0.04374, p =.3700). Figure 4B–K represents the correlation graphs of the tumor microenvironment.

4. Discussion

Anoctamins (ANO) are a group of anion channels with eight transmembrane domains, and numerous cellular functions [20]. The coding gene family, formerly known as transmembrane proteins with 16 domains (TMEM16), was only just discovered in 2014 thanks to bioinformatic analysis [31]. The new name (anoctamins), proposed by Yang et al., has replaced TMEM16 in GenBank despite objections, and it has been given the HUGO nomenclature seal of approval [31,32]. This family of transmembrane proteins constitutes ten paralogues (ANO1 through ANO10/ TMEM16A to TMEM16K [excluding the J alphabet from the naming]) dispersed across various human tissues [20]. All anoctamins feature eight hydrophobic helices, which have been previously hypothesized to be transmembrane domains according to hydropathy analysis, however, these findings are still debatable [31]. Duran et al. reported that, unlike ANO1 and ANO2 which have a clear Ca+2 activated C-l channel (CaCC) functionality, other members did not pursue such function, since they could not produce a C-l current through Ca+2 activation [33]. Furthermore, ANOs 3 through 7 were determined to be intracellular proteins probably residing in the endoplasmic reticulum without being trafficked to the cell membrane [33]. Despite being scarce, this family has been the subject of an expanding body of literature due to its association with numerous pathologies. To illustrate, ANO1 has been linked to a variety of malignancies. Along with ANO5, which has been connected to specific types of muscular dystrophy, ANO6 and ANO10 have also been linked to Scott syndrome and autosomal recessive spinocerebellar ataxia, respectively [31].
Anoctamins have been involved in a wide variety of cellular functions, including epithelial cell secretion, neuronal activation, smooth muscle contractions, skeletal muscle membrane repair, sensory transduction, and carcinogenesis [20,31,33]. On the molecular level, they behave as CaCCs, along with a scramblase activity, in which they express phospholipids across the membrane from the cytoplasmic side to the extracellular side while requiring Ca+2 [20,31,33]. ANO1 has been the most extensively studied ANO protein and has been shown to mediate Cl−1 secretion in secretory epithelia of the respiratory, gastrointestinal and renal systems, and sweat glands [34]. Under the activation of noxious heat, its role in heat sensation via somatosensory neurons is obvious [35]. Furthermore, ANO1 regulates vascular and bronchial smooth muscle tone. Because of this, ANO1 may play a role in the pathophysiological mechanisms underlying hypertension and asthma, respectively [34,36]. According to Sun et al., the interaction between RANK and ANO1 in osteoclasts causes increased bone resorption and decreased bone mass. In individuals with osteoporosis, this makes ANO1 a viable therapeutic target [37]. ANO1 antagonist can have a synergistic effect with commercially used tocolytics on human uterine smooth muscles [38]. Finally, and perhaps most importantly, it is believed to escalate multiple carcinogenic processes including cellular proliferation, migration, and metastasis [39]. ANO2, a member of the ANO family that is expressed in olfactory sensory cells of the olfactory epithelium may have a function in the olfaction process [40]. It has been found that ANO3 is expressed in the dorsal root ganglia, which controls nociception [34]. ANO5 is most abundantly expressed in the musculoskeletal system including bones, chondrocytes, cardiac, and skeletal muscles [41]. The role of ANO6 as a procoagulant has been confirmed through human and animal studies resulting in Cl−1 influx and swelling of platelets [42]. ANO6 has a well-known scramblase activity, which is also shared with ANO3, ANO4, and ANO7 [43]. Overall, except for ANO1, which has been thoroughly investigated, other members of the ANO family need even more in-depth research, and we are only at the beginning.
ANO1, which is found on locus 11q13, is amplified in multiple tumors, particularly in gastrointestinal stromal tumors (GIST) and head and neck squamous cell carcinoma (HNSCC), in addition to numerous others, such as lung adenocarcinoma, chondroblastoma, esophageal squamous cell carcinoma, salivary gland tumors, oral squamous cell carcinoma, leiomyosarcoma of the uterus, glioma, breast, colorectal and prostate cancer [39,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61]. ANO1 was a cancer biomarker prior to it is identification as a chloride channel in GIST known as DOG1, and it has received several other names including ORAOV2, and TAOS-2 among oncologists [31,39]. However, its overexpression has, unfortunately, been a poor prognostic factor, leading to greater mortality rates, tumor growth, and invasiveness, as well as its positive correlation with distant metastasis, migration, and tumor grading [45,46,47,50,61]. The carcinogenesis orchestrated by ANO1 has implicated numerous cell signaling pathways. ANO1-mediated proliferation was shown by Duvvuri et al. to be accompanied by extracellular signal-regulated kinase (ERK)1/2 activation, cyclin D1 evocation, and mitogen-activated protein kinase (MAPK) activation [45]. ANO1 upregulates calcium/calmodulin-dependent protein kinase II (CaMKII) and epidermal growth factor receptor (EGFR) expression, the latter of which regulates the MAPK or PI3K-AKT pathway. Therefore, ANO1 plays a role in governing cell variability through EGFR-AKT/SRC/MAPK and CaMKII signaling pathways [34]. ezrin-radixin-moesin proteins form a cross-link between the cell membrane and actin filaments of the cytoskeleton, hence its involvement in cell migration. Together with its physical association with ANO1, this might provide a clue to the role of ANO1 in EGF-driven migratory and invasive properties [34]. In addition to its connection to the cytoskeleton, its ability to control cell size aids in the development of shrunken cells that can move via diapedesis through inter-endothelial gaps [39,62]. Another intriguing association lies between the sonic hedgehog signaling pathway and ANO1 which is known to coordinate cellular growth and differentiation [39,63].
Anaplastic thyroid carcinoma appears to acquire its undesirable traits and aggressiveness through ANO1 overexpression; ANO1 knockdown greatly reduces the tumor’s aggressive behavior [64]. In gastric cancer, ANO5 has been linked to a negative prognostic role; when it is knocked down, it causes apoptosis, reduces cell proliferation, and arrests the cell cycle at the G1/S transition [65]. According to research by Pan et al., osteosarcomas express ANO5, a pro-tumorigenic factor that increases tumor size, grade, and metastasis. The instability and destruction of nel-like proteins 1 and 2 enable such activity [66]. ANO5 expression is scarce in healthy pancreatic tissue, but it is elevated in pancreatic cancer, where it contributes to the disease’s proliferation and migration [67]. ANO7 is only expressed in prostate cells, has an unknown function, and is negatively correlated with the prognosis of prostate cancer. Its low expression has been linked to positive surgical margin, lymph node metastasis, high classical and quantitative Gleason grades, advanced tumor stage, high Ki67 labelling index, and early biochemical recurrence [68].
We focused our bioinformatic research on ANO4 because it is one of the ANO family members that has not received enough attention. It is mostly expressed in the cervix, ovaries, prostate, adrenal glands, and central nervous system [69]. Given this, it has been suggested that it contributes to a variety of neurological diseases, such as schizophrenia, multiple sclerosis, Alzheimer’s disease, and anxiety disorders [70,71,72,73,74]. Reichhart et al. reported that ANO4 acts as a Ca2+-dependent phospholipid scramblase and monovalent nonselective ion channel [69]. Leitzke et al. reported ANO4 to modulate disintegrin-like metalloproteases ADAM 10 and 17 sheddase activity, evident through the increasing ADAM10 and 17 substrates: transforming growth factor alpha (TGF-α), amphiregulin (AREG), and betacellulin. They also demonstrated that the effects brought upon by the overexpression of ANO4 are due to it is scramblase activity, consequently resulting in diminishing AREG and increased cellular proliferation [75]. Maniero et al. recognized ANO4 as a significant gene expressed in the zona glomerulosa cells despite its contradictory effects on aldosterone secretion [76,77]. To the best of our knowledge, this is the first article to relate the expression of ANO4 with tumor involvement.
In conclusion, ANO4 expression was upregulated in tumor samples compared to normal kidney tissue. Albeit the latter finding, low ANO4 expression is associated with advanced clinicopathological variables such as tumor grade, stage, and pT. In addition, low ANO4 expression is linked to shorter OS, PFI, and DSS. Multivariate Cox logistic regression analysis identified ANO4 expression as an independent prognostic variable in OS (HR: 1.686, 95% CI: 1.120–2.540, p = 0.012), PFI (HR: 1.727, 95% CI: 1.103–2.704, p = 0.017), and DSS (HR: 2.688, 95% CI: 1.465–4.934, p = 0.001). GSEA identified the following pathways to be enriched within the low ANO4 expression group: epithelial–mesenchymal transition, G2-M checkpoint, E2F targets, estrogen response, apical junction, glycolysis, hypoxia, coagulation, KRAS, complement, p53, myogenesis, and TNF-α signaling via NF-κB pathways. ANO4 expression correlates significantly with monocyte (ρ = −0.1429, p = 0.0033) and mast cell (ρ = 0.1598, p = 0.001) infiltration.

Author Contributions

Conceptualization, T.E.-E. and B.M.A.M.; Methodology, A.H.A.S., Y.O.A.Z., K.A.-K. and I.H.I.; Software, A.H.A.S., K.A.-K. and A.A.A.M.; Validation, Y.O.A.Z., A.A.L., I.H.I., S.A.S. and F.Q.A.; Formal analysis, B.M.A.M.; Investigation, A.H.A.S. and A.A.L.; Resources, A.A.L., A.A.A.M. and F.Q.A.; Data curation, A.H.A.S. and Y.O.A.Z.; Writing—original draft, A.H.A.S., Y.O.A.Z., K.A.-K. and S.A.S.; Writing—review & editing, T.E.-E., I.H.I., B.M.A.M. and F.Q.A.; Visualization, T.E.-E.; Project administration, F.Q.A.; Funding acquisition, F.Q.A. All authors have read and agreed to the published version of the manuscript.

Funding

The article processing charges were funded jointly by Qatar National Library and Qatar University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The leading author would thank the 6th year medical school surgery final exam for giving the time and inspiration to learn the basis of bioinformatics which was the core nucleus of this work. In addition, the authors thank Qatar National Library for providing the open access support of this article.

Conflicts of Interest

All the authors have no conflict of interest to declare.

References

  1. Cancer Facts & Figures 2022. Am. Cancer Soc. 2022, 80. Available online: https://www.cancer.org/research/cancer-facts-statistics/all-cancer-facts-figures/cancer-facts-figures-2022.html (accessed on 1 October 2022).
  2. Petejova, N.; Martínek, A. Renal cell carcinoma: Review of etiology, pathophysiology and risk factors. Biomed. Pap. 2016, 160, 183–194. [Google Scholar] [CrossRef] [PubMed]
  3. Riveros, C.; Ranganathan, S.; Xu, J.; Chang, C.; Kaushik, D.; Morgan, M.; Miles, B.J.; Muhammad, T.; Anis, M.; Aghazadeh, M.; et al. Comparative real-world survival outcomes of metastatic papillary and clear cell renal cell carcinoma treated with immunotherapy, targeted therapy, and combination therapy. Urol. Oncol. Semin. Orig. Investig. 2023. [Google Scholar] [CrossRef] [PubMed]
  4. Baldewijns, M.M.L.; van Vlodrop, I.J.H.; Schouten, L.J.; Soetekouw, P.M.M.B.; de Bruïne, A.P.; van Engeland, M. Genetics and epigenetics of renal cell cancer. Biochim. Biophys. Acta (BBA)—Rev. Cancer 2008, 1785, 133–155. [Google Scholar] [CrossRef] [PubMed]
  5. Kase, A.M.; George, D.J.; Ramalingam, S. Clear Cell Renal Cell Carcinoma: From Biology to Treatment. Cancers 2023, 15, 665. [Google Scholar] [CrossRef]
  6. Moch, H.; Cubilla, A.L.; Humphrey, P.A.; Reuter, V.E.; Ulbright, T.M. The 2016 WHO Classification of Tumours of the Urinary System and Male Genital Organs—Part A: Renal, Penile, and Testicular Tumours. Eur. Urol. 2016, 70, 93–105. [Google Scholar] [CrossRef]
  7. Kovacs, G.; Akhtar, M.; Beckwith, B.J.; Bugert, P.; Cooper, C.S.; Delahunt, B.; Eble, J.N.; Fleming, S.; Ljungberg, B.; Medeiros, L.J.; et al. The Heidelberg classification of renal cell tumours. J. Pathol. 1997, 183, 131–133. [Google Scholar] [CrossRef]
  8. Shuch, B.; Amin, A.; Armstrong, A.J.; Eble, J.N.; Ficarra, V.; Lopez-Beltran, A.; Martignoni, G.; Rini, B.I.; Kutikov, A. Understanding Pathologic Variants of Renal Cell Carcinoma: Distilling Therapeutic Opportunities from Biologic Complexity. Eur. Urol. 2015, 67, 85–97. [Google Scholar] [CrossRef]
  9. John, A.; Spain, L.; Hamid, A.A. Navigating the Current Landscape of Non-Clear Cell Renal Cell Carcinoma: A Review of the Literature. Curr. Oncol. 2023, 30, 923–937. [Google Scholar] [CrossRef]
  10. Hahn, A.W.; Lebenthal, J.; Genovese, G.; Sircar, K.; Tannir, N.M.; Msaouel, P. The significance of sarcomatoid and rhabdoid dedifferentiation in renal cell carcinoma. Cancer Treat. Res. Commun. 2022, 33, 100640. [Google Scholar] [CrossRef]
  11. Nabi, S.; Kessler, E.R.; Bernard, B.; Flaig, T.W.; Lam, E.T. Renal cell carcinoma: A review of biology and pathophysiology. F1000Research 2018, 7, 307. [Google Scholar] [CrossRef] [Green Version]
  12. Brugarolas, J. Molecular Genetics of Clear-Cell Renal Cell Carcinoma. J. Clin. Oncol. 2014, 32, 1968–1976. [Google Scholar] [CrossRef] [PubMed]
  13. Clark, D.J.; Dhanasekaran, S.M.; Petralia, F.; Pan, J.; Song, X.; Hu, Y.; da Veiga Leprevost, F.; Reva, B.; Lih, T.-S.M.; Chang, H.-Y.; et al. Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma. Cell 2019, 179, 964–983.e31. [Google Scholar] [CrossRef] [PubMed]
  14. Hakimi, A.A.; Pham, C.G.; Hsieh, J.J. A clear picture of renal cell carcinoma. Nat. Genet. 2013, 45, 849–850. [Google Scholar] [CrossRef] [PubMed]
  15. Kapur, P.; Peña-Llopis, S.; Christie, A.; Zhrebker, L.; Pavía-Jiménez, A.; Rathmell, W.K.; Xie, X.-J.; Brugarolas, J. Effects on survival of BAP1 and PBRM1 mutations in sporadic clear-cell renal-cell carcinoma: A retrospective analysis with independent validation. Lancet Oncol. 2013, 14, 159–167. [Google Scholar] [CrossRef]
  16. Blanco, A.I.; Teh, B.S.; Amato, R.J. Role of Radiation Therapy in the Management of Renal Cell Cancer. Cancers 2011, 3, 4010–4023. [Google Scholar] [CrossRef]
  17. Diamond, E.; Molina, A.; Carbonaro, M.; Akhtar, N.; Giannakakou, P.; Tagawa, S.; Nanus, D. Cytotoxic chemotherapy in the treatment of advanced renal cell carcinoma in the era of targeted therapy. Crit. Rev. Oncol. 2015, 96, 518–526. [Google Scholar] [CrossRef]
  18. Hsieh, J.J.; Le, V.H.; Oyama, T.; Ricketts, C.J.; Ho, T.H.; Cheng, E.H. Chromosome 3p Loss-Orchestrated VHL, HIF, and Epigenetic Deregulation in Clear Cell Renal Cell Carcinoma. J. Clin. Oncol. 2018, 36, 3533–3539. [Google Scholar] [CrossRef]
  19. Linehan, W.M.; Ricketts, C.J. The Cancer Genome Atlas of renal cell carcinoma: Findings and clinical implications. Nat. Rev. Urol. 2019, 16, 539–552. [Google Scholar] [CrossRef]
  20. Benarroch, E.E. Anoctamins (TMEM16 proteins). Neurology 2017, 89, 722–729. [Google Scholar] [CrossRef]
  21. Picollo, A.; Malvezzi, M.; Accardi, A. TMEM16 Proteins: Unknown Structure and Confusing Functions. J. Mol. Biol. 2015, 427, 94–105. [Google Scholar] [CrossRef] [Green Version]
  22. Brunner, J.D.; Lim, N.K.; Schenck, S.; Duerst, A.; Dutzler, R. X-ray structure of a calcium-activated TMEM16 lipid scramblase. Nature 2014, 516, 207–212. [Google Scholar] [CrossRef] [PubMed]
  23. Pedemonte, N.; Galietta, L.J. Structure and Function of TMEM16 Proteins (Anoctamins). Physiol. Rev. 2014, 94, 419–459. [Google Scholar] [CrossRef] [PubMed]
  24. Le, S.C.; Liang, P.; Lowry, A.J.; Yang, H. Gating and Regulatory Mechanisms of TMEM16 Ion Channels and Scramblases. Front. Physiol. 2021, 12, 787773. [Google Scholar] [CrossRef] [PubMed]
  25. Goldman, M.J.; Craft, B.; Hastie, M.; Repečka, K.; McDade, F.; Kamath, A.; Banerjee, A.; Luo, Y.; Rogers, D.; Brooks, A.N.; et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat. Biotechnol. 2020, 38, 675–678. [Google Scholar] [CrossRef]
  26. Camp, R.L.; Dolled-Filhart, M.; Rimm, D.L. X-Tile: A New Bio-Informatics Tool for Biomarker Assessment and Outcome-Based Cut-Point Optimization. Clin. Cancer Res. 2004, 10, 7252–7259. [Google Scholar] [CrossRef]
  27. Tang, F.; Lu, Z.; He, C.; Zhang, H.; Wu, W.; He, Z. 53 years old is a reasonable cut-off value to define young and old patients in clear cell renal cell carcinoma: A study based on TCGA and SEER database. BMC Cancer 2021, 21, 638. [Google Scholar] [CrossRef]
  28. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef]
  29. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2018, 47, D607–D613. [Google Scholar] [CrossRef]
  30. Aran, D.; Hu, Z.; Butte, A.J. xCell: Digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017, 18, 220. [Google Scholar] [CrossRef]
  31. Duran, C.; Hartzell, H.C. Physiological roles and diseases of tmem16/anoctamin proteins: Are they all chloride channels? Acta Pharmacol. Sin. 2011, 32, 685–692. [Google Scholar] [CrossRef] [Green Version]
  32. Yang, Y.D.; Cho, H.; Koo, J.Y.; Tak, M.H.; Cho, Y.; Shim, W.-S.; Park, S.P.; Lee, J.; Lee, B.; Kim, B.-M.; et al. TMEM16A confers receptor-activated calcium-dependent chloride conductance. Nature 2008, 455, 1210–1215. [Google Scholar] [CrossRef] [PubMed]
  33. Duran, C.; Qu, Z.; Osunkoya, A.O.; Cui, Y.; Hartzell, H.C. ANOs 3–7 in the anoctamin/Tmem16 ClChannel family are intracellular proteins. Am. J. Physiol.-Cell Physiol. 2011, 302, C482–C493. [Google Scholar] [CrossRef] [PubMed]
  34. Oh, U.; Jung, J. Cellular functions of TMEM16/anoctamin. Pflügers Arch.—Eur. J. Physiol. 2016, 468, 443–453. [Google Scholar] [CrossRef]
  35. Cho, H.; Oh, U. Anoctamin 1 Mediates Thermal Pain as a Heat Sensor. Curr. Neuropharmacol. 2013, 11, 641–651. [Google Scholar] [CrossRef] [PubMed]
  36. Huang, F.; Zhang, H.; Wu, M.; Yang, H.; Kudo, M.; Peters, C.J.; Woodruff, P.G.; Solberg, O.D.; Donne, M.L.; Huang, X.; et al. Calcium-activated chloride channel TMEM16A modulates mucin secretion and airway smooth muscle contraction. Proc. Natl. Acad. Sci. USA 2012, 109, 16354–16359. [Google Scholar] [CrossRef]
  37. Sun, W.; Guo, S.; Li, Y.; Li, J.; Liu, C.; Chen, Y.; Wang, X.; Tan, Y.; Tian, H.; Wang, C.; et al. Anoctamin 1 controls bone resorption by coupling ClChannel activation with RANKL-RANK signaling transduction. Nat. Commun. 2022, 13, 2899. [Google Scholar] [CrossRef]
  38. Hyuga, S.; Parry, R.C.; Danielsson, J.; Vink, J.; Fu, X.W.; Wu, A.; Dan, W.; Yim, P.D.; Gallos, G. Anoctamin 1 antagonism potentiates conventional tocolytic-mediated relaxation of pregnant human uterine smooth muscle. J. Physiol. Sci. 2021, 71, 7. [Google Scholar] [CrossRef] [PubMed]
  39. Wanitchakool, P.; Wolf, L.; Koehl, G.E.; Sirianant, L.; Schreiber, R.; Kulkarni, S.; Duvvuri, U.; Kunzelmann, K. Role of anoctamins in cancer and apoptosis. Philos. Trans. R. Soc. B: Biol. Sci. 2014, 369, 20130096. [Google Scholar] [CrossRef]
  40. Stephan, A.B.; Shum, E.Y.; Hirsh, S.; Cygnar, K.D.; Reisert, J.; Zhao, H. ANO2 is the cilial calcium-activated chloride channel that may mediate olfactory amplification. Proc. Natl. Acad. Sci. USA 2009, 106, 11776–11781. [Google Scholar] [CrossRef]
  41. Mizuta, K.; Tsutsumi, S.; Inoue, H.; Sakamoto, Y.; Miyatake, K.; Miyawaki, K.; Noji, S.; Kamata, N.; Itakura, M. Molecular characterization of GDD1/TMEM16E, the gene product responsible for autosomal dominant gnathodiaphyseal dysplasia. Biochem. Biophys. Res. Commun. 2007, 357, 126–132. [Google Scholar] [CrossRef]
  42. Fernández, D.I.; Kuijpers, M.J.E.; Heemskerk, J.W.M. Platelet calcium signaling by G-protein coupled and ITAM-linked receptors regulating anoctamin-6 and procoagulant activity. Platelets 2020, 32, 863–871. [Google Scholar] [CrossRef] [PubMed]
  43. Suzuki, J.; Fujii, T.; Imao, T.; Ishihara, K.; Kuba, H.; Nagata, S. Calcium-dependent Phospholipid Scramblase Activity of TMEM16 Protein Family Members *. J. Biol. Chem. 2013, 288, 13305–13316. [Google Scholar] [CrossRef] [PubMed]
  44. Dixit, R.; Kemp, C.; Kulich, S.; Seethala, R.; Chiosea, S.; Ling, S.; Ha, P.K.; Duvvuri, U. TMEM16A/ANO1 is differentially expressed in HPV-negative versus HPV-positive head and neck squamous cell carcinoma through promoter methylation. Sci. Rep. 2015, 5, 16657. [Google Scholar] [CrossRef] [PubMed]
  45. Duvvuri, U.; Shiwarski, D.J.; Xiao, D.; Bertrand, C.; Huang, X.; Edinger, R.S.; Rock, J.R.; Harfe, B.D.; Henson, B.J.; Kunzelmann, K.; et al. TMEM16A Induces MAPK and Contributes Directly to Tumorigenesis and Cancer Progression. Cancer Res. 2012, 72, 3270–3281. [Google Scholar] [CrossRef]
  46. Ayoub, C.; Wasylyk, C.; Li, Y.; Thomas, E.; Marisa, L.; Robé, A.; Roux, M.; Abecassis, J.; de Reyniès, A.; Wasylyk, B. ANO1 amplification and expression in HNSCC with a high propensity for future distant metastasis and its functions in HNSCC cell lines. Br. J. Cancer 2010, 103, 715–726. [Google Scholar] [CrossRef]
  47. Ruiz, C.; Martins, J.R.; Rudin, F.; Schneider, S.; Dietsche, T.; Fischer, C.A.; Tornillo, L.; Terracciano, L.M.; Schreiber, R.; Bubendorf, L.; et al. Enhanced Expression of ANO1 in Head and Neck Squamous Cell Carcinoma Causes Cell Migration and Correlates with Poor Prognosis. PLoS ONE 2012, 7, e43265. [Google Scholar] [CrossRef]
  48. Jia, L.; Liu, W.; Guan, L.; Lu, M.; Wang, K. Inhibition of Calcium-Activated Chloride Channel ANO1/TMEM16A Suppresses Tumor Growth and Invasion in Human Lung Cancer. PLoS ONE 2015, 10, e0136584. [Google Scholar] [CrossRef]
  49. Wu, H.; Guan, S.; Sun, M.; Yu, Z.; Zhao, L.; He, M.; Zhao, H.; Yao, W.; Wang, E.; Jin, F.; et al. Ano1/TMEM16A Overexpression Is Associated with Good Prognosis in PR-Positive or HER2-Negative Breast Cancer Patients following Tamoxifen Treatment. PLoS ONE 2015, 10, e0126128. [Google Scholar] [CrossRef]
  50. Britschgi, A.; Bill, A.; Brinkhaus, H.; Rothwell, C.; Clay, I.; Duss, S.; Rebhan, M.; Raman, P.; Guy, C.T.; Wetzel, K.; et al. Calcium-activated chloride channel ANO1 promotes breast cancer progression by activating EGFR and CAMK signaling. Proc. Natl. Acad. Sci. USA 2013, 110, E1026–E1034. [Google Scholar] [CrossRef]
  51. Sui, Y.; Sun, M.; Wu, F.; Yang, L.; Di, W.; Zhang, G.; Zhong, L.; Ma, Z.; Zheng, J.; Fang, X.; et al. Inhibition of TMEM16A Expression Suppresses Growth and Invasion in Human Colorectal Cancer Cells. PLoS ONE 2014, 9, e115443. [Google Scholar] [CrossRef]
  52. Sauter, D.R.P.; Novak, I.; Pedersen, S.F.; Larsen, E.H.; Hoffmann, E.K. ANO1 (TMEM16A) in pancreatic ductal adenocarcinoma (PDAC). Pflügers Arch.-Eur. J. Physiol. 2015, 467, 1495–1508. [Google Scholar] [CrossRef] [PubMed]
  53. Akpalo, H.; Lange, C.; Zustin, J. Discovered on gastrointestinal stromal tumour 1 (DOG1): A useful immunohistochemical marker for diagnosing chondroblastoma. Histopathology 2012, 60, 1099–1106. [Google Scholar] [CrossRef] [PubMed]
  54. Chênevert, J.; Duvvuri, U.; Chiosea, S.; Dacic, S.; Cieply, K.; Kim, J.; Shiwarski, D.; Seethala, R.R. DOG1: A novel marker of salivary acinar and intercalated duct differentiation. Mod. Pathol. 2012, 25, 919–929. [Google Scholar] [CrossRef] [PubMed]
  55. Kashyap, M.; Marimuthu, A.; Kishore, C.J.H.; Peri, S.; Keerthikumar, S.; Prasad, T.S.K.; Mahmood, R.; Rao, S.; Ranganathan, P.; Sanjeeviah, R.C.; et al. Genomewide mRNA profiling of esophageal squamous cell carcinoma for identification of cancer biomarkers. Cancer Biol. Ther. 2009, 8, 36–46. [Google Scholar] [CrossRef]
  56. Shi, Z.-Z.; Shang, L.; Jiang, Y.-Y.; Hao, J.-J.; Zhang, Y.; Zhang, T.-T.; Lin, D.-C.; Liu, S.-G.; Wang, B.-S.; Gong, T.; et al. Consistent and Differential Genetic Aberrations between Esophageal Dysplasia and Squamous Cell Carcinoma Detected By Array Comparative Genomic Hybridization. Clin. Cancer Res. 2013, 19, 5867–5878. [Google Scholar] [CrossRef]
  57. Huang, X.; Gollin, S.M.; Raja, S.; Godfrey, T.E. High-resolution mapping of the 11q13 amplicon and identification of a gene, TAOS1, that is amplified and overexpressed in oral cancer cells. Proc. Natl. Acad. Sci. USA 2002, 99, 11369–11374. [Google Scholar] [CrossRef]
  58. Aloi, M.; Lionetti, P.; Barabino, A.; Guariso, G.; Costa, S.; Fontana, M.; Romano, C.; Lombardi, G.; Miele, E.; Alvisi, P.; et al. Phenotype and Disease Course of Early-onset Pediatric Inflammatory Bowel Disease. Inflamm. Bowel Dis. 2014, 20, 597–605. [Google Scholar] [CrossRef]
  59. Sah, S.P.; McCluggage, W.G. DOG1 immunoreactivity in uterine leiomyosarcomas. J. Clin. Pathol. 2012, 66, 40–43. [Google Scholar] [CrossRef]
  60. Liu, J.; Liu, Y.; Ren, Y.; Kang, L.; Zhang, L. Transmembrane protein with unknown function 16A overexpression promotes glioma formation through the nuclear factor-κB signaling pathway. Mol. Med. Rep. 2014, 9, 1068–1074. [Google Scholar] [CrossRef]
  61. Liu, W.; Lu, M.; Liu, B.; Huang, Y.; Wang, K. Inhibition of Ca2+-activated Cl channel ANO1/TMEM16A expression suppresses tumor growth and invasiveness in human prostate carcinoma. Cancer Lett. 2012, 326, 41–51. [Google Scholar] [CrossRef]
  62. Cuddapah, V.A.; Sontheimer, H. Ion channels and tranporters in cancer. 2. Ion channels and the control of cancer cell migration. Am. J. Physiol. Physiol. 2011, 301, C541–C549. [Google Scholar] [CrossRef] [PubMed]
  63. Rock, J.R.; Lopez, M.C.; Baker, H.; Harfe, B.D. Identification of genes expressed in the mouse limb using a novel ZPA microarray approach. Gene Expr. Patterns 2007, 8, 19–26. [Google Scholar] [CrossRef] [PubMed]
  64. Kim, J.-Y.; Youn, H.Y.; Choi, J.; Baek, S.K.; Kwon, S.Y.; Eun, B.K.; Park, J.-Y.; Oh, K.H. Anoctamin-1 affects the migration and invasion of anaplastic thyroid carcinoma cells. Anim. Cells Syst. 2019, 23, 294–301. [Google Scholar] [CrossRef] [PubMed]
  65. Fukami, T.; Shiozaki, A.; Kosuga, T.; Kudou, M.; Shimizu, H.; Ohashi, T.; Arita, T.; Konishi, H.; Komatsu, S.; Kubota, T.; et al. Anoctamin 5 regulates the cell cycle and affects prognosis in gastric cancer. World J. Gastroenterol. 2022, 28, 4649–4667. [Google Scholar] [CrossRef] [PubMed]
  66. Pan, R.; Lu, Q.; Ren, C.; Li, H.; Zeng, F.; Tian, X.; Chen, H. Anoctamin 5 promotes osteosarcoma development by increasing degradation of Nel-like proteins 1 and 2. Aging 2021, 13, 17316–17327. [Google Scholar] [CrossRef]
  67. Song, H.-Y.; Zhou, L.; Hou, X.-F.; Lian, H. Anoctamin 5 regulates cell proliferation and migration in pancreatic cancer. Int. J. Clin. Exp. Pathol. 2019, 12, 4263–4270. [Google Scholar]
  68. Marx, A.; Koopmann, L.; Höflmayer, D.; Büscheck, F.; Hube-Magg, C.; Steurer, S.; Eichenauer, T.; Clauditz, T.S.; Wilczak, W.; Simon, R.; et al. Reduced anoctamin 7 (ANO7) expression is a strong and independent predictor of poor prognosis in prostate cancer. Cancer Biol. Med. 2021, 18, 245–255. [Google Scholar] [CrossRef]
  69. Reichhart, N.; Schöberl, S.; Keckeis, S.; Alfaar, A.S.; Roubeix, C.; Cordes, M.; Crespo-Garcia, S.; Haeckel, A.; Kociok, N.; Föckler, R.; et al. Anoctamin-4 is a bona fide Ca2+-dependent non-selective cation channel. Sci. Rep. 2019, 9, 2257. [Google Scholar] [CrossRef]
  70. Hendrickx, D.A.E.; van Scheppingen, J.; van der Poel, M.; Bossers, K.; Schuurman, K.G.; van Eden, C.G.; Hol, E.M.; Hamann, J.; Huitinga, I. Gene Expression Profiling of Multiple Sclerosis Pathology Identifies Early Patterns of Demyelination Surrounding Chronic Active Lesions. Front. Immunol. 2017, 8, 1810. [Google Scholar] [CrossRef]
  71. Sherva, R.; Tripodis, Y.; Bennett, D.A.; Chibnik, L.B.; Crane, P.K.; de Jager, P.L.; Farrer, L.A.; Saykin, A.J.; Shulman, J.M.; Naj, A. Genome-wide association study of the rate of cognitive decline in Alzheimer’s disease. Alzheimers Dement 2014, 10, 45–52. [Google Scholar] [CrossRef]
  72. Webb, B.T.; Guo, A.-Y.; Maher, B.S.; Zhao, Z.; van den Oord, E.J.; Kendler, K.S.; Riley, B.P.; Gillespie, N.A.; Prescott, C.A.; Middeldorp, C.M.; et al. Meta-analyses of genome-wide linkage scans of anxiety-related phenotypes. Eur. J. Hum. Genet. 2012, 20, 1078–1084. [Google Scholar] [CrossRef] [PubMed]
  73. Terracciano, A.; Sanna, S.; Uda, M.; Deiana, B.; Usala, G.; Busonero, F.; Maschio, A.; Scally, M.; Patriciu, N.; Chen, W.-M.; et al. Genome-wide association scan for five major dimensions of personality. Mol. Psychiatry 2008, 15, 647–656. [Google Scholar] [CrossRef] [PubMed]
  74. Athanasiu, L.; Mattingsdal, M.; Kähler, A.K.; Brown, A.; Gustafsson, O.; Agartz, I.; Giegling, I.; Muglia, P.; Cichon, S.; Rietschel, M.; et al. Gene variants associated with schizophrenia in a Norwegian genome-wide study are replicated in a large European cohort. J. Psychiatr. Res. 2010, 44, 748–753. [Google Scholar] [CrossRef] [PubMed]
  75. Leitzke, S.; Seidel, J.; Ahrens, B.; Schreiber, R.; Kunzelmann, K.; Sperrhacke, M.; Bhakdi, S.; Reiss, K. Influence of Anoctamin-4 and-9 on ADAM10 and ADAM17 Sheddase Function. Membranes 2022, 12, 123. [Google Scholar] [CrossRef] [PubMed]
  76. Maniero, C.; Scudieri, P.; Shaikh, L.H.; Zhao, W.; Gurnell, M.; Galietta, L.J.; Brown, M.J. ANO4 (Anoctamin 4) Is a Novel Marker of Zona Glomerulosa That Regulates Stimulated Aldosterone Secretion. Hypertension 2019, 74, 1152–1159. [Google Scholar] [CrossRef]
  77. Maniero, C.; Zhou, J.; Shaikh, L.H.; Azizan, E.A.B.; McFarlane, I.; Neogi, S.; Scudieri, P.; Galietta, L.J.; Brown, M.J. Role of ANO4 in regulation of aldosterone secretion in the zona glomerulosa of the human adrenal gland. Lancet 2015, 385, S62. [Google Scholar] [CrossRef]
Figure 1. Differential expression of ANO4 between normal kidney tissues and tumor samples (A) and its age and gender-matched comparison (B). The correlation between ANO4 expression and clinicopathological variables including age (C), gender (D), grade (E), stage (F), pT (G), and pN (H). Kaplan–Meier curves demonstrating the impact of ANO4 expression on the OS (I), PFI (J), DSS (K), and DFI (L). Data are presented as mean ± SEM.
Figure 1. Differential expression of ANO4 between normal kidney tissues and tumor samples (A) and its age and gender-matched comparison (B). The correlation between ANO4 expression and clinicopathological variables including age (C), gender (D), grade (E), stage (F), pT (G), and pN (H). Kaplan–Meier curves demonstrating the impact of ANO4 expression on the OS (I), PFI (J), DSS (K), and DFI (L). Data are presented as mean ± SEM.
Jpm 13 00295 g001
Figure 2. Kaplan–Meier curves demonstrating the impact of ANO4 expression on the OS within subgroups of dichotomized demographics and clinicopathological variables including age (A,B), gender (C,D), grade (E,F), stage (G,H), pT (I,J), and pN (K,L).
Figure 2. Kaplan–Meier curves demonstrating the impact of ANO4 expression on the OS within subgroups of dichotomized demographics and clinicopathological variables including age (A,B), gender (C,D), grade (E,F), stage (G,H), pT (I,J), and pN (K,L).
Jpm 13 00295 g002
Figure 3. GSEA plots illustrating significantly enriched pathways in low (AM) and high (NP) ANO4 expression groups.
Figure 3. GSEA plots illustrating significantly enriched pathways in low (AM) and high (NP) ANO4 expression groups.
Jpm 13 00295 g003
Figure 4. Bubble plot illustrating the top 15 KEGG pathways enriched within the low and high ANO4 expression groups (A). Top five biological processes (BP), cellular components (CC), and molecular functions (MF) enriched within the low and high ANO4 expression groups (B). * Denotes for an FDR > 25%.
Figure 4. Bubble plot illustrating the top 15 KEGG pathways enriched within the low and high ANO4 expression groups (A). Top five biological processes (BP), cellular components (CC), and molecular functions (MF) enriched within the low and high ANO4 expression groups (B). * Denotes for an FDR > 25%.
Jpm 13 00295 g004
Figure 5. PPI network of ANO4 based on string analysis (A). The correlation between ANO4 expression and immune microenvironment (BK).
Figure 5. PPI network of ANO4 based on string analysis (A). The correlation between ANO4 expression and immune microenvironment (BK).
Jpm 13 00295 g005
Table 1. Baseline clinicopathological characteristics of the TCGA non-metastasized ccRCC cohort (n = 422).
Table 1. Baseline clinicopathological characteristics of the TCGA non-metastasized ccRCC cohort (n = 422).
CharacteristicData
Age 61 (51–71)
Gender
  Male275 (65.2)
  Female147 (34.8)
Stage
  1243 (57.6)
  254 (12.8)
  3121 (28.7)
  44 (0.9)
Grade
  110 (2.4)
  2204 (48.3)
  3162 (38.4)
  439 (9.2)
  Unknown 7 (1.7)
pT
  T1245 (58.1)
  T255 (13.0)
  T3119 (28.3)
  T43 (0.7)
pN
  N0200 (47.4)
  N110 (2.4)
  NX212 (50.2)
Data are presented as n (%) or median (IQR).
Table 2. Baseline clinicopathological characteristics of the TCGA non-metastasized ccRCC cohort according to ANO4 expression.
Table 2. Baseline clinicopathological characteristics of the TCGA non-metastasized ccRCC cohort according to ANO4 expression.
CharacteristicsLow Expression (n = 126)High Expression (n = 296)χ2p-Value
Age
  ≤ 5330 (7.1)100 (23.7)4.1250.042
  > 5396 (22.7)196 (46.4)
Gender
  Male95 (22.5)180 (42.7)8.2830.004
  Female31 (7.3)116 (27.5)
Grade
  G1 + G255 (13.3)159 (38.3)3.6820.055
  G3 + G469 (16.6)132 (31.8)
Stage
  1 + 277 (18.2)220 (52.1)7.4020.007
  3 + 449 (11.6)76 (18.0)
pT stage
  T1 + T279 (18.7)221 (52.4)6.1550.013
  T3 + T447 (11.1)75 (17.8)
pN stage
  N065 (15.4)135 (32.0)1.2670.260
  N1 + NX61 (14.5)161 (38.2)
Data are presented as n (%).
Table 3. Univariate and multivariate Cox logistic regression analysis evaluating the utility of ANO4 expression and relevant clinicopathological covariates in predicting the overall survival.
Table 3. Univariate and multivariate Cox logistic regression analysis evaluating the utility of ANO4 expression and relevant clinicopathological covariates in predicting the overall survival.
Covariates *Univariate Cox Logistic Regression Multivariate Cox Logistic Regression
HR95% CIp-ValueHR95% CIp-Value
OS
  Age (≤53 vs. >53)0.4180.225–0.6860.0010.5100.305–0.8520.010
  Gender (Female vs. Male)1.1110.754–1.6360.5961.1160.738–1.6880.602
  Grade (G1 + G2 vs. G3 + G4)0.5590.378–0.8270.0040.6890.453–1.0480.082
  Stage (1 + 2 vs. 3 + 4)0.3880.266–0.566<0.0010.5960.078–4.5290.617
  pT (T1 + T2 vs. T3 + T4)0.3860.265–0.564<0.0010.8910.116–6.8700.912
  pN (N0 vs. N1 + NX)1.1980.821–1.7490.3491.0900.738–1.6100.665
  ANO4 expression (Low vs. High)2.0991.434–3.073<0.0011.6861.120–2.5400.012
PFI
  Age (≤53 vs. >53)0.5880.360–0.9600.0340.6550.388–1.1070.114
  Gender (Female vs. Male)0.6210.388–0.9920.0460.6270.381–1.0300.065
  Grade (G1 + G2 vs. G3 + G4)0.4440.284–0.693<0.0010.6280.392–1.0070.054
  Stage (1 + 2 vs. 3 + 4)0.2490.163–0.379<0.0010.1290.092–0.5790.008
  pT (T1 + T2 vs. T3 + T4)0.2640.174–0.402<0.0012.4720.545–11.2150.214
  pN (N0 vs. N1 + NX)1.1120.733–1.6870.6171.1170.727–1.7180.613
  ANO4 expression (Low vs. High)2.5211.654–3.840<0.0011.7271.103–2.7040.017
DSS
  Age (≤53 vs. >53)0.0970.248–0.9950.0480.6810.325–1.4230.307
  Gender (Female vs. Male)0.7290.392–1.3570.3190.8390.436–1.6150.599
  Grade (G1 + G2 vs. G3 + G4)0.3280.174–0.6190.0010.5190.260–1.0350.063
  Stage (1 + 2 vs. 3 + 4)0.2010.112–0.359<0.0010.2170.026–1.8270.160
  pT (T1 + T2 vs. T3 + T4)0.2100.118–0.374<0.0011.4150.166–12.1010.751
  pN (N0 vs. N1 + NX)1.2970.738–2.2780.3661.1460.639–2.0050.648
  ANO4 expression (Low vs. High)3.8942.201–6.889<0.0012.6881.465–4.9340.001
DFI
  Age (≤53 vs. >53)1.4940.534–4.1770.4441.9330.619–6.0240.257
  Gender (Female vs. Male)0.4850.154–1.5260.2160.6500.183–2.3070.505
  Grade (G1 + G2 vs. G3 + G4)0.5520.182–1.6740.2930.9840.277–3.4960.980
  Stage (1 + 2 vs. 3 + 4)0.3110.112–0.8650.0250.2580.079–0.8480.026
  pT (T1 + T2 vs. T3 + T4)0.3110.112–0.8650.025--
  pN (N0 vs. N1 + NX)0.7680.271–2.1770.6201.1020.346–3.5060.869
  ANO4 expression (Low vs. High)1.4380.482–4.2930.5151.5430.471–5.0590.474
* Initial sub-covariate was used a reference.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Al Sharie, A.H.; Al Zu’bi, Y.O.; El-Elimat, T.; Al-Kammash, K.; Abu Lil, A.; Isawi, I.H.; Al Sharie, S.; Abu Mousa, B.M.; Al Malkawi, A.A.; Alali, F.Q. ANO4 Expression Is a Potential Prognostic Biomarker in Non-Metastasized Clear Cell Renal Cell Carcinoma. J. Pers. Med. 2023, 13, 295. https://doi.org/10.3390/jpm13020295

AMA Style

Al Sharie AH, Al Zu’bi YO, El-Elimat T, Al-Kammash K, Abu Lil A, Isawi IH, Al Sharie S, Abu Mousa BM, Al Malkawi AA, Alali FQ. ANO4 Expression Is a Potential Prognostic Biomarker in Non-Metastasized Clear Cell Renal Cell Carcinoma. Journal of Personalized Medicine. 2023; 13(2):295. https://doi.org/10.3390/jpm13020295

Chicago/Turabian Style

Al Sharie, Ahmed H., Yazan O. Al Zu’bi, Tamam El-Elimat, Kinda Al-Kammash, Alma Abu Lil, Israa H. Isawi, Sarah Al Sharie, Balqis M. Abu Mousa, Abubaker A. Al Malkawi, and Feras Q. Alali. 2023. "ANO4 Expression Is a Potential Prognostic Biomarker in Non-Metastasized Clear Cell Renal Cell Carcinoma" Journal of Personalized Medicine 13, no. 2: 295. https://doi.org/10.3390/jpm13020295

APA Style

Al Sharie, A. H., Al Zu’bi, Y. O., El-Elimat, T., Al-Kammash, K., Abu Lil, A., Isawi, I. H., Al Sharie, S., Abu Mousa, B. M., Al Malkawi, A. A., & Alali, F. Q. (2023). ANO4 Expression Is a Potential Prognostic Biomarker in Non-Metastasized Clear Cell Renal Cell Carcinoma. Journal of Personalized Medicine, 13(2), 295. https://doi.org/10.3390/jpm13020295

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop