Aberrations in ion channels interacting with lipid metabolism and epithelial–mesenchymal transition in esophageal squamous cell carcinoma

Abstract

Esophageal squamous cell carcinoma (ESCC) is the most prevalent malignant gastrointestinal tumor. Ion channels contribute to tumor growth and progression through interactions with their neighboring molecules including lipids. The dysregulation of membrane ion channels and lipid metabolism may contribute to the epithelial–mesenchymal transition (EMT), leading to metastatic progression. Herein, transcriptome profiles of patients with ESCC were analyzed by performing differential gene expression and weighted gene co-expression network analysis to identify the altered ion channels, lipid metabolism- and EMT-related genes in ESCC. A total of 1,081 differentially expressed genes, including 113 ion channels, 487 lipid metabolism-related, and 537 EMT-related genes, were identified in patients with ESCC. Thereafter, EMT scores were correlated with altered co-expressed genes. The altered co-expressed genes indicated a correlation with EMT signatures. Interactions among 22 ion channels with 3 hub lipid metabolism- and 13 hub EMT-related proteins were determined using protein–protein interaction networks. A pathway map was generated to depict deregulated signaling pathways including insulin resistance and the estrogen receptor-Ca²⁺ signaling pathway in ESCC. The relationship between potential ion channels and 5-year survival rates in ESCC was determined using Kaplan–Meier plots and Cox proportional hazard regression analysis. Inositol 1,4,5-trisphosphate receptor type 3 (ITPR3) was found to be associated with poor prognosis of patients with ESCC. Additionally, drugs interacting with potential ion channels, including GJA1 and ITPR3, were identified. Understanding alterations in ion channels with lipid metabolism and EMT in ESCC pathophysiology would most likely provide potential targets for the better treatment of patients with ESCC.

1 Introduction

Esophageal squamous cell carcinoma (ESCC) is the most prominent subtype of esophageal cancer worldwide, with a significant increasing incidence rate (Abnet et al., 2018). Several prognostic strategies have been introduced to improve the survival rates of patients with ESCC (Hirano and Kato, 2019). Hitherto, negative results from clinical trials of targeted therapies suggest the need to identify novel biomarkers in regulating significant molecular mechanisms of ESCC progression (Hirano and Kato, 2019). Alterations in ion channels, the membrane proteins that allow the influx and efflux of ions across membranes, have been reported in several tumors (Sharma et al., 2019; Parthasarathi et al., 2022; Aydar et al., 2009; Chioni et al., 2010). Roles of ion channels in apoptosis and their potential in cell regulation and tumor progression were identified (Comes et al., 2015; Fraser et al., 2014; Bustin et al., 2001). The voltage-gated K⁺, Cl⁻, and Ca²⁺, and transient receptor potential (TRP) channels contributed in the progression of several tumors (Shiozaki et al., 2021; Litan and Langhans, 2015). Potassium voltage-gated channels, KCNA1 (Northcott et al., 2012) and KCNA3 (Comes et al., 2013), have been reported to be overexpressed in different malignancies such as medulloblastoma, breast, colon, and prostate cancers. In addition, the growth and migration of glioma cells were observed as regulated by the upregulation of splice variants of large potassium channels (Liu et al., 2002).

Recently, it has been reported that EMT leads to changed plasma membrane ion channel expression (Azimi and Monteith, 2016). Aquaporins, Cl⁻ and Ca²⁺ channels, were found to be involved in the progression of ESCC (Shiozaki et al., 2021). Furthermore, the progression of esophageal cancer was identified to be regulated by different types of ion channels, including voltage-gated K⁺, Cl⁻, and Ca²⁺, and TRP channels expressed by esophageal cancer cells and tissues (Shiozaki et al., 2021).

Alterations in lipid metabolism would result in variations in membrane composition, protein distribution, and cellular functions. This further results in the development and progression of various diseases, including cystic fibrosis, neurological disorders, and cancer (Huang and Freter, 2015). Furthermore, various studies indicate the role of lipids in enhancing migration and proliferation of tumor cells (Wang et al., 2021a; Deng et al., 2019; Kouba et al., 2019). The lipid metabolism-related genes, such as CYP2C9, DGAT1, UGT1A6, INS, HPGDS, and LPL, were identified as prognostically significant in patients with lung adenocarcinoma (Li et al., 2020). Chen and colleagues reported preoperative serum total cholesterol and low-density lipoprotein cholesterol as prognostic factors in ESCC patients who underwent esophagectomy (Chen et al., 2017). Moreover, cellular lipid metabolism involves ion channels and may influence the functions of other ion channels. The conversion of phosphatidylinositol bisphosphate to diacylglycerol or sphingomyelin to ceramide has been reported in both the maintenance of membrane surface potential and opening and closing of channel proteins (Kasimova et al., 2014), (Minke and Cook, 2002), (Huang and Freter, 2015).

The presence of ion channels, lipid metabolism- and EMT-related genes on the membrane enhances the possibility of various interactions between them and their role as a union in several key cellular processes. The alterations in the membrane potential may lead to changes in the lateral spatiotemporal distribution of charged lipids present in the membrane (Downward, 2003; Accardi, 2015). This could further affect the membrane-bound signaling proteins like Ras family GTPases (Downward, 2003; Accardi, 2015). Kouba et al. discussed on devising novel therapeutic strategies by understanding the inter-dependencies between Ca²⁺ and Ca²⁺-activated K⁺ channels with lipid mediators in epithelial ovarian cancer (Kouba et al., 2019). However, currently there is no study to indicate the relationship of ion channels with lipid metabolism-related genes and the EMT program in ESCC.

In this study, transcriptome sequencing (RNA-seq) data were obtained from 12 matched paired samples of patients with ESCC from the Kidwai Memorial Institute of Oncology, Bangalore, India (IOB-KMIO), and the publicly available dataset GSE32424 from the Gene Expression Omnibus (GEO) database was used to identify alterations co-expressions between ion channels, lipid metabolism- and EMT-related genes in ESCC. Thereafter, the deregulated genes were correlated with EMT using EMT scoring techniques. Subsequently, protein–protein interaction (PPI) network and clustering analyses were performed. Furthermore, the expression profiles of ESCC patients from The Cancer Genome Atlas (TCGA) and healthy individuals from Genotype-Tissue Expression (GTEx) data resources were used to compare the identified deregulated ion channels interacting with lipid metabolism- and EMT-related genes in a larger cohort of patients. In total, 1,081 unique genes were found to be differentially expressed in the transcriptome of patients with ESCC that included 113 ion channels, and 487 lipid metabolism- and 537 EMT-related genes. The interactome consisting of the interactions between the three sets of proteins that included 16 significantly altered ion channels was generated. Gap junction alpha-1 (GJA1), transient receptor potential cation channel subfamily V member 3 (TRPV3), and ITPR3 were found interacting with hub proteins AKT1, CCL2, CDH2, FN1, EGFR, and CAV1 and were a part of dense clusters within the PPI network. A pathway map was generated to depict the involvement of ion channels interacting with lipid metabolism- and EMT-related genes in various pathways in ESCC. Deregulation of several pathways, including insulin signaling and estrogen receptor-calcium signaling, was identified in ESCC. ITPR3 was found to be prognostically significant. Additionally, drugs interacting with the potential ion channels were identified. The analyses could identify the importance of interactions between ion channels and lipid metabolism-related genes that might orchestrate the processes involved in EMT. Experimental validation of the potential ion channels is needed to further elucidate their roles as potential targets in ESCC.

2 Materials and methodology

2.1 Institutional Review Board statement

The study was approved by the Institutional Review Board Medical Ethics Committee of the Kidwai Memorial Institute of Oncology, Bangalore, India (KMIO/MEC/021/24 dated 30 November 2016), and conducted according to the guidelines of the Declaration of Helsinki.

2.2 Informed consent statement

Informed consent was obtained from all subjects involved in the study.

2.3 Data collection

2.3.1 Gene list collection

The list of ion channels was downloaded from the HGNC database (accessed on 11th March 2022), lipid metabolism-related genes were collected from the literature (Li et al., 2020; Fernandez et al., 2020; Dunn et al., 2019; Moskot et al., 2015; Meng et al., 2021; El-Fasakhany et al., 2001; Russo et al., 2018; Lutz and Cornett, 2013; Astudillo et al., 2012; Sonnweber et al., 2018; Hanna and Hafez, 2018), and EMT-related genes, from dbEMT (v2.0). The total number of genes present in each gene category is provided in Supplementary Table S1.

2.3.2 RNA-seq data collection

Twelve matched paired samples of patients with ESCC were collected from the Kidwai Memorial Institute of Oncology (IOB-KMIO), Bangalore, India. Supplementary File S1 contains the experimental procedure used for transcriptome sequencing. Another RNA-seq dataset, GSE32424, obtained using the Illumina Genome Analyzer IIx platform was downloaded from the GEO database. Furthermore, the gene expression profiles of samples corresponding to ESCC from the TCGA_ESCA cohort were downloaded from the UCSC Xena portal (Goldman et al., 2020) using R software (v4.1.3) UCSCXenaShiny library (Wang et al., 2021b). The healthy samples of the esophagus were downloaded from a GTEx data resource using the UCSCXenaShiny library. The detailed description on all the samples included in these datasets and the criteria used for the selection of healthy samples from GTEx is described in Supplementary Tables S2–S5. The RNA-seq datasets from IOB-KMIO and GSE32424 were used as the discovery datasets, and the TCGA-GTEx dataset was utilized to compare the deregulated genes in a larger cohort of patients with ESCC.

2.4 Data processing

The raw reads from IOB-KMIO and GSE32424 datasets were analyzed using the workflow depicted in Supplementary Figure S1. The raw reads were quality checked using FastQC, quality trimmed using Trimmomatic (v0.39), and clean reads were obtained. Furthermore, clean reads were mapped and aligned to the human genome GRCh38.105 using the STAR RNA-seq aligner (v2.7.10a) (Supplementary Tables S6, S7). Furthermore, the aligned reads were quantified using HTSeq (v1.99.2). R (v4.1) scripts were used to merge the HTSeq output files, and the R Bioconductor package, biomaRt (Durinck et al., 2009), was used to convert the Ensembl IDs to HGNC gene symbols. The ESCC RNA-seq log (read count +1) normalized data from TCGA were transformed to raw read counts using R scripts. Furthermore, the read counts corresponding to ion channels, lipid metabolism- and EMT-related genes were parsed using R scripts. Supplementary Figure S2 depicts the overall workflow used for the identification of ion channels and their interactions with lipid metabolism and EMT.

2.5 Identification of differentially expressed genes

Differentially expressed ion channels, lipid metabolism- and EMT-related genes in the IOB-KMIO and GSE32424 datasets were identified using the R Bioconductor package DESeq2 (v4.1.2) (Love et al., 2014). An adjusted p-value (padj) ≤ 0.05 and log2-fold change (log2FC) > 0.6 were set as the criteria for the selection of upregulated genes, and the genes with padj ≤ 0.05 and log2FC < −0.6 were selected as downregulated genes. A non-redundant list with all DEGs obtained from the two datasets was generated.

2.6 Construction of co-expression networks of DEGs

The expression profiles of ESCC patients from IOB-KMIO and GSE32424 datasets corresponding to DEGs were parsed. Furthermore, these profiles were provided as input to the R Bioconductor package weighted gene co-expression network analysis (WGCNA) (Zhang and Horvath, 2005; Langfelder and Horvath, 2008) to construct the scale-free co-expression networks. The read counts provided as input were DESeq2-normalized. On the basis of a scale-free network model, an adjacency correlation matrix was calculated (Langfelder and Horvath, 2008). A suitable soft threshold power that ensured a scale-free network was selected based on the pickSoftThreshold function. Thereafter, a signed topological overlap matrix (TOM) and the corresponding dissimilarity matrix (1-TOM) were derived from the obtained adjacency matrix. The dissimilarity matrix was used for performing hierarchical clustering with the DynamicTreeCut algorithm (Langfelder et al., 2008) in order to group the genes with similar expression profiles into the same gene modules.

The expression profiles of each obtained modules were epitomized by module eigengenes (MEs) and correlated with binary trait (normal and tumor) values, and a relationship between the module and traits was obtained. The p-value was calculated as a confidence measure for the module–trait relationship. Furthermore, gene significance (GS) and module membership (MM) values were calculated to estimate the association of individual genes in the module with the binary trait and the correlation of gene expression profiles with MEs, respectively. The intramodular gene connectivity, with the binary trait, was estimated on the basis of GS vs. MM plots (Langfelder and Horvath, 2008). Thereafter, the most correlated modules with a significant p-value were selected from each dataset. The obtained weighted gene co-expression networks were exported to Cytoscape with the weight threshold value set to 0.02 for visualization.

2.7 Calculation of EMT scores

In this study, 76-gene EMT signature-based (GS76), multinomial logistic regression-based (MLR), and Kolmogorov–Smirnov test-based (KS) methods developed in the past (Chakraborty et al., 2020; George et al., 2017; Mandal et al., 2021) were used for the calculation of EMT scores of transcriptomic profiles of samples from patients with ESCC. A higher KS and MLR score indicated a more mesenchymal sample, and a higher GS76 score indicated a more epithelial sample. MLR quantifies the extent of EMT on a scale of [0, 2], and KS scores are defined on a scale of [−1, +1] (Chakraborty et al., 2020). The scatter plots between different EMT scores were plotted using the R package “ggplot2” (Wickham, 2016) with correlation coefficients and p-values calculated using the R package “corrplot” (Wei et al., 2017) by Pearson’s method.

2.8 Computation of the correlation of altered ion channels and lipid metabolism-related genes with EMT

The estimated EMT scores of IOB-KMIO samples by the GS76, MLR, and KS methods were correlated with the expression profiles of the module ion channels and lipid metabolism-related genes interacting with EMT-related genes in ESCC.

2.9 Construction of the protein–protein interaction networks of co-expressed gene modules

To infer the interactions between the genes in the modules, protein–protein interaction networks (PPINs) were constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) (v11.0) database (von Mering et al., 2003; Szklarczyk et al., 2021). The interactions with a medium confidence and false discovery rate of 1% were further visualized in Cytoscape (v3.9) (Shannon et al., 2003). The unconnected nodes in the network were discarded. Thereafter, NetworkAnalyser, a Cytoscape plugin, was used to calculate the various properties of the network, including degree centrality and betweenness centrality. The top 15 nodes in the network with the highest degree centrality were selected as the hub nodes in the networks.

2.10 Identification of highly interacting clusters within the protein–protein interaction networks

The Molecular Complex Detection (MCODE) (Bader and Hogue, 2003) tool, a Cytoscape plugin, was used to identify the dense clusters representing the highly interacting nodes within the PPI networks. The parameter maximum depth was set to 100 degree, k-core was maintained as 2, and the node score was set to 0.2. Thereafter, the clusters consisting of ion channels, and lipid metabolism- and EMT-related nodes were selected to identify ion channels interacting with both lipid metabolism- and EMT-related nodes.

2.11 Generation of the pathway map of ion channels interacting with lipid metabolism- and EMT-related genes

An extensive literature search was carried out using PubMed, and reactions like activations, inhibitions, and translocations that describe the interactions between the clustered ion channels, lipid metabolism- and EMT-related genes were annotated. Thereafter, a signaling pathway map consisting of the deregulated ion channels, lipid metabolism- and EMT-related genes was generated using PathVisio (version 3.3.0) in the Graphical Pathway Markup Language (GPML) format.

2.12 Association of putative ion channels in the survival of patients with ESCC

TCGA-ESCA dataset was used to compare the alterations in the ion channels obtained from IOB-KMIO and GSE32424 datasets. Differentially expressed ion channels, lipid metabolism- and EMT-related genes in the RNA-seq TCGA-GTEx ESCC dataset were identified using the R Bioconductor package DESeq2 (v4.1.2). An adjusted p-value (padj) ≤ 0.05 and log2FC > 0.6 were set as the criteria for the selection of upregulated genes, and the genes with padj ≤ 0.05 and log2FC < −0.6 were selected as downregulated genes.

Furthermore, the clinical data on patients with ESCC were extracted using the R Bioconductor package “RTCGA” for survival analysis. To determine the correlation between ion channel expression and a 5-year survival rate of patients with ESCC, Cox proportional hazard regression analysis was performed. The R Bioconductor package “survival” was utilized to calculate the hazard ratio (HR) with a confidence interval of 95% and the log rank p-values. Kaplan–Meier (KM) survival plots were generated to visualize the differences between the high and low expressions of the putative ion channels using the R Bioconductor package “survminer”. The median value was used as the group cut-off criteria.

2.13 Identification of drugs interacting with ion channels correlating with lipid metabolism- and EMT-related genes

Ion channels are potential clinical targets for restorative intervention due to their localization to the plasma membrane and different regulation mechanisms (Anderson et al., 2019). The list of known drugs interacting with all the human genes was collected from The Drug–Gene Interaction Database (DGIdb) (Freshour et al., 2021). This list was parsed with the list of potential ion channels interacting with lipid metabolism- and EMT-related genes.

3 Results

3.1 DEGs in IOB-KMIO and GSE32424 datasets

A total of 282 genes, including 39 ion channels, 108 lipid metabolism genes, and 135 EMT-related genes, were found as differentially expressed in the IOB-KMIO dataset. In total, 1,013 genes, including 89 ion channels, 433 lipid metabolism genes, and 491 EMT-related genes, were found differentially expressed in the GSE32424 dataset (Table 1 and Supplementary Table S8).

Figure 1(i) represents the overlaps between the upregulated and downregulated genes in IOB-KMIO and GSE32424 datasets. A total of 1,081 DEGs were identified. Figure 1(ii) represents the overlaps between ion channels, and lipid metabolism- and EMT-related genes across IOB-KMIO and GSE32424 datasets.

FIGURE 1

3.2 Co-expression networks of DEGs

WGCNA of the identified DEGs led to the detection of gene modules consisting of genes co-expressed in each of the dataset. For the IOB-KMIO dataset, a soft thresholding power of 14 was chosen and the minimum module size of 50 with a deepSplit parameter set to 0 was used to construct co-expression networks. Similarly, for the GSE32424 dataset, a soft thresholding power of 18 was chosen and the minimum module size of 50 with a deepSplit parameter set to 0 was used to construct co-expression networks. The IOB-KMIO dataset was split into five modules. Based on the p-value and correlation values, the blue and brown modules were selected (Figure 2). GSE32424 was split into two modules—turquoise and blue. Both the modules showed a good correlation among the genes within the modules (Supplementary Figure S3).

FIGURE 2

3.3 Correlation of module genes with EMT

The EMT scores obtained using GS76, MLR, and KS metrics for the samples in IOB-KMIO are provided in Supplementary Table S9. A negative GS76 score with a positive KS score and a higher MLR score indicates the sample to be a mesenchymal sample. Similarly, the sample with a positive GS76 score and a negative KS score with a lower MLR score could be an epithelial sample. The correlation between the scores corresponding to the respective method in the dataset was calculated and is depicted in Figure 3 and Supplementary Table 10. MLR and KS scores correlated positively with each other. GS76 scores correlated negatively with both MLR and KS scores.

FIGURE 3

3.4 Protein–protein interaction networks of co-expressed gene modules

PPI networks obtained through the STRING database for the co-expressed gene modules revealed the predicted PPIs between co-expressed DEGs (Supplementary Figure S4). The lipid metabolism- and EMT-related proteins with multiple interactions were identified as hub nodes in the network (Figure 4; Supplementary Figure S5). Several ion channels were found interacting with the hub nodes in the network. In the network generated from the blue module of the KMIO dataset, the ion channels GJA1, P2RX4, and GJC1 interacted with the hub nodes AKT1, CCL2, CDH2, and FN1. TRPV3 and ITPR3 were found interacting with EGFR in the network generated from the KMIO brown module dataset. In the network from the GSE32424 turquoise module dataset, 13 ion channels, TRPV3, GJB1, GJC1, P2RX4, ITPR3, TRPC6, RYR2, CLIC4, LRRC8A, TRPC1, GJA5, TRPV2, and KCNMA1, interacted with seven hub nodes that included EGFR, CDH2, IL6, AKT1, CAV1, EP300, and IGF1. VDAC2, ITPR2, VDAC1, GJB2, and CFTR were found interacting with HSPA4, MAPK3, and CDH1 in the network from the GSE32424 blue module (Table 2).

FIGURE 4

Furthermore, PPI networks were analyzed using the MCODE algorithm to identify dense networks within the PPI networks. In total, 11 clusters in the IOB-KMIO blue module PPI network, 10 clusters in the IOB-KMIO brown module PPI network, 26 clusters in the GSE32424 turquoise module PPI network, and 20 clusters in the GSE32424 blue module PPI network were identified. These network clusters represented the highly interacting protein clusters within the PPINs. Out of these, five clusters—one from the IOB-KMIO blue module, one from the IOB-KMIO brown module, two from the GSE32424 turquoise module, and one from the GSE32424 blue module—comprised ion channels together with lipid metabolism- and EMT-related proteins (Figures 5A–E) (Supplementary Table S11). Several ion channels present in the cluster including GABRR2, ANO1, and TRPV3 were found to be negatively correlated with GS76 and positively correlated with both MLR and KS scores, indicating that these ion channels have a higher probability toward stabilizing/acquiring a mesenchymal phenotype. Similarly, GABRQ, GABRA3, and TRPM7 negatively correlated with GS76 and positively correlated with MLR and KS scores, indicating that they may have a higher affinity toward maintaining the epithelial phenotype (Supplementary Table S9).

FIGURE 5

3.5 Pathway map of identified altered ion channels, lipid metabolism-related, and EMT-related genes

The significance of the interactions between ion channels, lipid metabolism- and EMT-related genes in various aspects of ESCC tumorigenesis can be laid down by understanding their role in cellular processes. Pathway maps enable the depiction of these cellular processes that may aid in better understanding the disease etiology (Sharma et al., 2020). Figure 6 summarizes the reaction events that are likely to occur in ESCC patients through interactions between deregulated ion channels, lipid metabolism- and EMT-related genes. The pathway map comprises of 77 molecules and 46 reactions. Key cellular processes, including insulin signaling and estrogen receptor-calcium signaling pathways deregulated in ESCC, are depicted in the pathway map.

FIGURE 6

3.6 Prognostically significant ion channels present in the highly interacting clusters

The deregulated ion channels in IOB-KMIO and GSE32424 were compared with the deregulated ion channels in TCGA-ESCC dataset. In total, 240 ion channels were found deregulated in TCGA-ESCC dataset, of which 18 were found deregulated in IOB-KMIO and 49 were present in the deregulated list of ion channels from the GSE32424 dataset. In total, 11 ion channels, GJB2, GRIN2D, TPCN1, KCTD6, RYR2, HTR3B, SCNN1B, CLIC3, KCNE3, CNGA1, and CLCN1, were identified as deregulated in all three datasets (Supplementary Figure S6). Out of 16 ion channels that interacted with EMT- and lipid metabolism-related genes, 14 were found deregulated in at least two datasets. TRPM2 and CLCN1 were found deregulated in all three datasets (Supplementary Table S12).

The interconnection of ion channels with the survival of patients with ESCC was estimated through KM plots and Cox proportional hazard regression analysis. Ion channels with HR > 1 were selected as major ion channels for the survival of ESCC patients (Figure 7). Six ion channels—GJA1 [HR: 3.24, p (HR): 0.07], GABRQ [HR: 1.68, p (HR): 0.19], GABRA3 [HR: 2.42, p (HR): 0.11], KCNN4 [HR: 3.88, p (HR): 0.04], TRPC1 [HR: 2.72, p (HR): 0.09], and ITPR3 [HR: 4.09, p(HR): 0.04]—exhibited HR > 1. Patients with high expression of ITPR3 (p-value: 0.03) had a better survival rate. Additionally, 782 gene–drug interactions were identified between the 16 ion channels—ANO1, CLCN1, CLCN2, GABRA3, GABRA4, GABRE, GABRP, GABRQ, GABRR2, GJA1, ITPR3, KCNN4, TRPC1, TRPM2, TRPM7, and TRPV3 drugs (Supplementary Table 13).

FIGURE 7

4 Discussion

Alterations in ion channels affect various cellular processes, such as cell volume regulation, migration, and cell proliferation, that act as a form of sustenance in tumorigenesis (Lang and Stournaras, 2014; Comes et al., 2015). Lipid metabolism is a complex process that involves lipid intake, transport, synthesis, and destruction. Diverse signaling pathways control the biosynthesis and breakdown of different lipids under physiological and pathological conditions (Huang and Freter, 2015). These signaling pathways can either be activated or inhibited depending on the needs of the cell and on responses to the changing cellular environment (Huang and Freter, 2015). Furthermore, the EMT-related genes that form the basis of a metastatic transition in tumor cells are associated with the functionalities of both the ion channels and the lipid mediators. Alterations in these molecules may impact the interactions between proteins, resulting in deregulations of the downstream processes. This could activate a cascade of uncontrolled reactions within the cells in a tissue. Elucidating the role of ion channels with lipid metabolism and EMT would aid future studies to come up with better treatment options for the management of ESCC.

The current study aimed to identify potential altered ion channels interacting with lipid metabolism- and EMT-related genes in ESCC. Differentially expressed ion channels, lipid metabolism- and EMT-related genes from the transcriptomic profiles of patients with ESCC were identified. Furthermore, WGCNA led to the identification of significantly co-expressed gene modules that included the three sets of genes. The correlation of the genes in each module with EMT was obtained using EMT scoring techniques that were consistent across the dataset modules. PPI networks further aided in the identification of PPIs among the differentially co-expressed gene modules. Several ion channels were found to interact with the hub proteins in the network.

Three ion channels GJA1, TRPV3, and ITPR3 were found interacting with the hub proteins in the network and clustered together with lipid metabolism- and EMT-related proteins in the highly interacting clusters.

GJA1 interacted with the hub proteins AKT1, CCL2, CDH2, and FN1 and with the other proteins—SOCS3, POMC, MED13L, MED30, VIM, MED24, MED12, THRAP3, MED1, MED27, MED17, aldehyde dehydrogenase 1 family member A1 (ALDH1A1), APOE, mTOR, LIF, LPL, MED22, FABP4, and MED14 (Figures 4A, 5A, D). GJA1 is the key connexin member present in the immune system and is ubiquitously expressed (Gleisner et al., 2017). It interacts with diverse proteins and can affect cancer cell phenotypes, including cell growth and metastasis (Wang et al., 2018; Aasen et al., 2019). The role of GJA1 has been identified in breast cancer bone metastasis (Wang et al., 2018); its interactions with p21-activated protein kinase 1 to activate mitogen-activated protein kinase (MAPK) signaling in HeLa cells (Aasen et al., 2019) and its upregulation leading to micro-metastasis and tumor vasculature in tumor cell–endothelial cell contact areas (Elzarrad et al., 2008) lead to evidence of GJA1 as a potential target.

TRPV3 interacted with the hub protein EGFR and other proteins—ARSB, NPC1, SOAT1, GLA, LPGAT1, LIPA, CSNK2A2, STS, ST3GAL5, CYP27A1, B4GALNT1, NEU3, and ABCC1 (Figures 4B, 5E). TRP channels, such as TRPV1, TRPV3, and TRPV4, can function as ion transport channels, secondary transducers of ionotropic or metabotropic receptors, or main detectors of physical and chemical stimuli (Holzer, 2011). TRPV3 was among the TRP channels found upregulated in the study carried out to estimate the immune microenvironment and prognostic effect of TRP channels in patients with ESCC (Zhao et al., 2022). In non-small cell lung carcinoma, the inhibition of TRPV3 led to the reduction in Ca²⁺ ions that arrested the lung cancer cells at the G1/S stage; reduced the expression of cyclinA, cyclinD1, cyclinE, and p-CaMKII; and increased the P27 level (Li et al., 2016). The association between 43 fatty acid metabolism genes and genetic variability in colorectal cancer indicated TRPV3 with four fatty acid metabolism genes to be associated with a higher risk of colorectal tumor progression (Hoeft et al., 2010), (Holzer, 2011).

ITPR3 was found interacting with EGFR, AKT1, CAV1, ARSB, NPC1, SOAT1, GLA, LIPA, CSNK2A2, LPGAT1, STS, ST3GAL5, CYP27A1, B4GALNT1, NEU3, and ABCC1 (Figures 4B, 5E). ITPR3 dysregulation may interfere with processes leading to endoplasmic reticulum stress (Zheng et al., 2022). The upregulation of ITPR3 was observed in several malignancies, such as cholangiocarcinoma (Ueasilamongkol et al., 2020), colon cancer (Shibao et al., 2010), melanoma (Kuchay et al., 2017), colorectal cancer (Shibao et al., 2010), cervical squamous cell carcinoma cancer (Yang et al., 2017), glioblastoma (Kang et al., 2010), hepatocellular carcinoma (HCC) (Guerra et al., 2019), mesothelioma, and prostate cancer (Bononi et al., 2017).The upregulation of ITPR3 plays an essential role in T-lymphocyte apoptosis (Khan et al., 1996). The interactions of ion channels with lipid mediators could possibly be an indication of their role in the dysregulated cascade of events in tumorigenic environments. Pathways contributing to ESCC tumorigenesis through ion channels that interacted with lipid metabolism- and EMT-related proteins included insulin resistance and estrogen receptor-calcium signaling (Figure 6).

4.1 Insulin resistance

GJA1 was found interacting with a suppressor of cytokine signaling 3 (SOCS3) and other proteins involved in the mechanism leading to tumor progression through insulin resistance-related pathways. Obesity is one of the major risk factors contributing toward the onset of ESCC (Etemadi et al., 2012). Elevated levels of IL-6 observed in obese individuals could lead to an increase in the levels of suppressor of cytokine signaling 1 (SOCS1) and SOCS3 in white adipose tissue and the liver. At the molecular level, this increase could impair the insulin functionality by binding of the SOCS proteins to insulin receptors, IRS1 and IRS2. Another mechanism leading to the development of insulin resistance by overexpression of the SOCS protein was through reduction in the tyrosine phosphorylation of IRS1 and IRS2 required for insulin transduction (Wunderlich et al., 2013). This signal transduction is carried out through phosphoinositide 3-kinases (PI3K) and the MAPK pathway accounting for mitogenesis, transformation, and apoptosis protection in tumor cells (Arcidiacono et al., 2021). Wunderlich and others have shown the involvement of SOCS1 and SOCS3 in insulin resistance through the elevation of interleukin-6 (IL-6), thus leading to the activation of the JAK-STAT pathway.

4.2 JAK/STAT pathway

GJA1 was found interacting with the signal transducer and activator of transcription (STAT), and other proteins involved in this pathway. STAT3 aberrant activation could lead to esophageal cellular transformation, angiogenesis, EMT, and metastasis. Sugase and colleagues have identified the association of the JAK/STAT pathway in the prognosis of patients with ESCC through immunohistochemistry studies (Sugase et al., 2017). The correlation between the phosphorylation levels of STAT3 and JAK1 was significant in predicting the poor prognosis of ESCC patients, indicating their importance in better management of ESCC (Rah et al., 2022).

4.3 PI3K-AKT pathway

ALDH1A1 was found to interact with the ion channels GABRE, GABRR2, GABRQ, GABRA3, ANO1, and GJA1 (Figure 5B). The PI3K-AKT pathway known to be deregulated in ESCC is activated by the overexpression of ALDH1A1. ALDH1A1 increases the expression of AKT serine/threonine kinase 1 (AKT1) and p-AKT. The phosphorylation of AKT activates several substrates that regulate cellular proliferation, survival, and motility (Shi et al., 2019). ALDH1A1 plays a significant role in maintaining stem cell properties in patients with ESCC through the AKT pathway. ALDH1A1 could be a potential target in patients with ESCC (Wang et al., 2020).

4.4 Estrogen receptor-Ca²⁺ signaling pathway

The estrogen receptor-Ca²⁺ signaling pathway was found to be associated with male predominance of ESCC. The estrogen receptors present on the plasma membrane activate multiple intracellular secondary messengers, including Ca²⁺ signaling. An increase in the cytoplasmic Ca²⁺ levels occur due to the release of Ca²⁺ from the endoplasmic reticulum through ion channels, such as ITPR3 (Zhang et al., 2017).

Cox proportional regression analysis with high HR indicated the alterations in the ion channels GJA1, GABRQ, GABRA3, KCNN4, TRPC1, and ITPR3 to possibly be associated with poor survival in patients with ESCC. Patients with low expression of GJA1 (p-value: 0.14), KCNN4 (p-value: 0.56), and TRPC1 (p-value: 0.55) and higher expression of GABRQ (p-value: 0.14), GABRA3 (p-value: 0.57), and ITPR3 (p-value: 0.032) may show better survival. However, the number of patients with ESCC was limited to attain statistical significance to support the expression difference using a log rank p-value <0.05 in survival analysis except ITPR3. Of these, GJA1 and ITPR3 were found interacting with hub proteins in the interactome and with other proteins in the highly interacting network clusters. Several ion channels found deregulated in the present study, including GJA1, and ITPR3, interacted with various known drugs. The presence of ion channels on the plasma membrane makes them accessible targets to therapeutic drugs. Moreover, their putative significance as diagnostic and prognostic markers and the availability of multiple drugs interacting with the potential ion channels make them ideal targets for clinical studies. In vivo mouse model studies showed that SKF96365, a phenylethylimidazole, and an inhibitor of store operated calcium channels (SOCE) resulted in decreased tumor volume by 50% in 4 weeks (Leanza et al., 2016). The inhibition of TRPC6 in ESCC cells leads to the suppression of proliferation and induction of the G2/M phase arrest (Zhang et al., 2018). The inhibition of TRPV2 by tranilast; inhibition of ANO1 by T16Ainh-A01, digallic acid, and tannic acids; and the inhibition of voltage-gated sodium channels by valproic acid were reported as potential targeted therapeutic agents in esophageal squamous cell carcinoma (Shiozaki et al., 2018; Capatina et al., 2022; Bill et al., 2014; Wheler et al., 2014). However, these channels are ubiquitously present, and toxicity toward normal cells on targeting these channels has not been addressed in depth. Further in vivo and in vitro studies could be carried out to elucidate the role of identified potentially altered ion channels interacting with genes of lipid metabolism- and EMT-related genes.

5 Conclusion

The interactome of ion channels with lipid metabolism- and EMT-related genes of patients with ESCC determined potential altered ion channels in ESCC. Furthermore, altered ion channels were found to be involved in several pathways, including insulin resistance, JAK-STAT, PI3K-AKT, and estrogen receptor-Ca²⁺ signaling. Additionally, ITPR3 was associated with poor prognosis of patients with ESCC. Targeting potential ion channel networks with inhibitors could most likely be an alternative therapeutic strategy to improve ESCC patient care. Moreover, experimental studies demonstrating the mechanistic view of these interactions need to be established.

References

• 1

  AasenT.LeitheE.GrahamS. V.KameritschP.MayanM. D.MesnilM.et al (2019). Connexins in cancer: Bridging the gap to the clinic. Oncogene38 (23), 4429–4451. 10.1038/s41388-019-0741-6

  • CrossRef
  • Google Scholar

• 2

  AbnetC. C.ArnoldM.WeiW. Q. (2018). Epidemiology of esophageal squamous cell carcinoma. Gastroenterology154 (2), 360–373. 10.1053/j.gastro.2017.08.023

  • CrossRef
  • Google Scholar

• 3

  AccardiA. (2015). Cell Signaling, Lipids link ion channels and cancer. Science349 (6250), 789–790. 10.1126/science.aad0874

  • CrossRef
  • Google Scholar

• 4

  AndersonK. J.CormierR. T.ScottP. M. (2019). Role of ion channels in gastrointestinal cancer. World J. Gastroenterol.25 (38), 5732–5772. 10.3748/wjg.v25.i38.5732

  • CrossRef
  • Google Scholar

• 5

  ArcidiaconoD.ZaramellaA.FabrisF.Sanchez-RodriguezR.NucciD.FassanM.et al (2021). Insulin/IGF-1 signaling is downregulated in barrett's esophagus patients undergoing a moderate calorie and protein restriction program: A randomized 2-year trial. Nutrients13 (10), 3638. 10.3390/nu13103638

  • CrossRef
  • Google Scholar

• 6

  AstudilloA. M.BalgomaD.BalboaM. A.BalsindeJ. (2012). Dynamics of arachidonic acid mobilization by inflammatory cells. Biochim. Biophys. Acta1821 (2), 249–256. 10.1016/j.bbalip.2011.11.006

  • CrossRef
  • Google Scholar

• 7

  AydarE.YeoS.DjamgozM.PalmerC. (2009). Abnormal expression, localization and interaction of canonical transient receptor potential ion channels in human breast cancer cell lines and tissues: A potential target for breast cancer diagnosis and therapy. Cancer Cell Int.9, 23. 10.1186/1475-2867-9-23

  • CrossRef
  • Google Scholar

• 8

  AzimiI.MonteithG. R. (2016). Plasma membrane ion channels and epithelial to mesenchymal transition in cancer cells. Endocr. Relat. Cancer23 (11), R517–R525. 10.1530/ERC-16-0334

  • CrossRef
  • Google Scholar

• 9

  BaderG. D.HogueC. W. (2003). An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinforma.4, 2. 10.1186/1471-2105-4-2

  • CrossRef
  • Google Scholar

• 10

  BillA.HallM. L.BorawskiJ.HodgsonC.JenkinsJ.PiechonP.et al (2014). Small molecule-facilitated degradation of ANO1 protein: A new targeting approach for anticancer therapeutics. J. Biol. Chem.289 (16), 11029–11041. 10.1074/jbc.M114.549188

  • CrossRef
  • Google Scholar

• 11

  BononiA.GiorgiC.PatergnaniS.LarsonD.VerbruggenK.TanjiM.et al (2017). BAP1 regulates IP3R3-mediated Ca(2+) flux to mitochondria suppressing cell transformation. Nature546 (7659), 549–553. 10.1038/nature22798

  • CrossRef
  • Google Scholar

• 12

  BustinS. A.LiS. R.DorudiS. (2001). Expression of the Ca2+-activated chloride channel genes CLCA1 and CLCA2 is downregulated in human colorectal cancer. DNA Cell Biol.20 (6), 331–338. 10.1089/10445490152122442

  • CrossRef
  • Google Scholar

• 13

  CapatinaA. L.LagosD.BrackenburyW. J. (2022). Targeting ion channels for cancer treatment: Current progress and future challenges. Rev. Physiol. Biochem. Pharmacol.183, 1–43. 10.1007/112_2020_46

  • CrossRef
  • Google Scholar

• 14

  ChakrabortyP.GeorgeJ. T.TripathiS.LevineH.JollyM. K. (2020). Comparative study of transcriptomics-based scoring metrics for the epithelial-hybrid-mesenchymal spectrum. Front. Bioeng. Biotechnol.8, 220. 10.3389/fbioe.2020.00220

  • CrossRef
  • Google Scholar

• 15

  ChenP.HanL.WangC.JiaY.SongQ.WangJ.et al (2017). Preoperative serum lipids as prognostic predictors in esophageal squamous cell carcinoma patients with esophagectomy. Oncotarget8 (25), 41605–41619. 10.18632/oncotarget.15651

  • CrossRef
  • Google Scholar

• 16

  ChioniA. M.ShaoD.GroseR.DjamgozM. B. (2010). Protein kinase A and regulation of neonatal Nav1.5 expression in human breast cancer cells: Activity-dependent positive feedback and cellular migration. Int. J. Biochem. Cell Biol.42 (2), 346–358. 10.1016/j.biocel.2009.11.021

  • CrossRef
  • Google Scholar

• 17

  ComesN.BielanskaJ.Vallejo-GraciaA.Serrano-AlbarrasA.MarruecosL.GomezD.et al (2013). The voltage-dependent K(+) channels Kv1.3 and Kv1.5 in human cancer. Front. Physiol.4, 283. 10.3389/fphys.2013.00283

  • CrossRef
  • Google Scholar

• 18

  ComesN.Serrano-AlbarrasA.CaperaJ.Serrano-NovilloC.CondomE.RamonY. C. S.et al (2015). Involvement of potassium channels in the progression of cancer to a more malignant phenotype. Biochim. Biophys. Acta1848 (10), 2477–2492. 10.1016/j.bbamem.2014.12.008

  • CrossRef
  • Google Scholar

• 19

  DengH.ZhouT.MoX.LiuC.YinY. (2019). Low-density lipoprotein promotes lymphatic metastasis of esophageal squamous cell carcinoma and is an adverse prognostic factor. Oncol. Lett.17 (1), 1053–1061. 10.3892/ol.2018.9683

  • CrossRef
  • Google Scholar

• 20

  DownwardJ. (2003). Targeting RAS signalling pathways in cancer therapy. Nat. Rev. Cancer3 (1), 11–22. 10.1038/nrc969

  • CrossRef
  • Google Scholar

• 21

  DunnT. M.TifftC. J.ProiaR. L. (2019). A perilous path: The inborn errors of sphingolipid metabolism. J. Lipid Res.60 (3), 475–483. 10.1194/jlr.S091827

  • CrossRef
  • Google Scholar

• 22

  DurinckS.SpellmanP. T.BirneyE.HuberW. (2009). Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat. Protoc.4 (8), 1184–1191. 10.1038/nprot.2009.97

  • CrossRef
  • Google Scholar

• 23

  El-FasakhanyF. M.UchimuraK.KannagiR.MuramatsuT. (2001). A novel human gal-3-O-sulfotransferase: Molecular cloning, characterization, and its implications in biosynthesis of (SO(4)-3)Galbeta1-4(fucalpha1-3)GlcNAc. J. Biol. Chem.276 (29), 26988–26994. 10.1074/jbc.M100348200

  • CrossRef
  • Google Scholar

• 24

  ElzarradM. K.HaroonA.WilleckeK.DobrowolskiR.GillespieM. N.Al-MehdiA. B. (2008). Connexin-43 upregulation in micrometastases and tumor vasculature and its role in tumor cell attachment to pulmonary endothelium. BMC Med.6, 20. 10.1186/1741-7015-6-20

  • CrossRef
  • Google Scholar

• 25

  EtemadiA.GolozarA.KamangarF.FreedmanN. D.ShakeriR.MatthewsC.et al (2012). Large body size and sedentary lifestyle during childhood and early adulthood and esophageal squamous cell carcinoma in a high-risk population. Ann. Oncol.23 (6), 1593–1600. 10.1093/annonc/mdr494

  • CrossRef
  • Google Scholar

• 26

  FernandezL. P.Gomez de CedronM.Ramirez de MolinaA. (2020). Alterations of lipid metabolism in cancer: Implications in prognosis and treatment. Front. Oncol.10, 577420. 10.3389/fonc.2020.577420

  • CrossRef
  • Google Scholar

• 27

  FraserS. P.Ozerlat-GunduzI.BrackenburyW. J.FitzgeraldE. M.CampbellT. M.CoombesR. C.et al (2014). Regulation of voltage-gated sodium channel expression in cancer: Hormones, growth factors and auto-regulation. Philos. Trans. R. Soc. Lond B Biol. Sci.369 (1638), 20130105. 10.1098/rstb.2013.0105

  • CrossRef
  • Google Scholar

• 28

  FreshourS. L.KiwalaS.CottoK. C.CoffmanA. C.McMichaelJ. F.SongJ. J.et al (2021). Integration of the drug-gene interaction database (DGIdb 4.0) with open crowdsource efforts. Nucleic Acids Res.49 (1), D1144–D1151. 10.1093/nar/gkaa1084

  • CrossRef
  • Google Scholar

• 29

  GeorgeJ. T.JollyM. K.XuS.SomarelliJ. A.LevineH. (2017). Survival outcomes in cancer patients predicted by a partial EMT gene expression scoring metric. Cancer Res.77 (22), 6415–6428. 10.1158/0008-5472.CAN-16-3521

  • CrossRef
  • Google Scholar

• 30

  GleisnerM. A.NavarreteM.HofmannF.Salazar-OnfrayF.TittarelliA. (2017). Mind the gaps in tumor immunity: Impact of connexin-mediated intercellular connections. Front. Immunol.8, 1067. 10.3389/fimmu.2017.01067

  • CrossRef
  • Google Scholar

• 31

  GoldmanM. J.CraftB.HastieM.RepeckaK.McDadeF.KamathA.et al (2020). Visualizing and interpreting cancer genomics data via the Xena platform. Nat. Biotechnol.38 (6), 675–678. 10.1038/s41587-020-0546-8

  • CrossRef
  • Google Scholar

• 32

  GuerraM. T.FlorentinoR. M.FrancaA.Lima FilhoA. C.Dos SantosM. L.FonsecaR. C.et al (2019). Expression of the type 3 InsP(3) receptor is a final common event in the development of hepatocellular carcinoma. Gut68 (9), 1676–1687. 10.1136/gutjnl-2018-317811

  • CrossRef
  • Google Scholar

• 33

  HannaV. S.HafezE. A. A. (2018). Synopsis of arachidonic acid metabolism: A review. J. Adv. Res.11, 23–32. 10.1016/j.jare.2018.03.005

  • CrossRef
  • Google Scholar

• 34

  HiranoH.KatoK. (2019). Systemic treatment of advanced esophageal squamous cell carcinoma: Chemotherapy, molecular-targeting therapy and immunotherapy. Jpn. J. Clin. Oncol.49 (5), 412–420. 10.1093/jjco/hyz034

  • CrossRef
  • Google Scholar

• 35

  HoeftB.LinseisenJ.BeckmannL.Muller-DeckerK.CanzianF.HusingA.et al (2010). Polymorphisms in fatty-acid-metabolism-related genes are associated with colorectal cancer risk. Carcinogenesis31 (3), 466–472. 10.1093/carcin/bgp325

  • CrossRef
  • Google Scholar

• 36

  HolzerP. (2011). Transient receptor potential (TRP) channels as drug targets for diseases of the digestive system. Pharmacol. Ther.131 (1), 142–170. 10.1016/j.pharmthera.2011.03.006

  • CrossRef
  • Google Scholar

• 37

  HuangC.FreterC. (2015). Lipid metabolism, apoptosis and cancer therapy. Int. J. Mol. Sci.16 (1), 924–949. 10.3390/ijms16010924

  • CrossRef
  • Google Scholar

• 38

  KangS. S.HanK. S.KuB. M.LeeY. K.HongJ.ShinH. Y.et al (2010). Caffeine-mediated inhibition of calcium release channel inositol 1,4,5-trisphosphate receptor subtype 3 blocks glioblastoma invasion and extends survival. Cancer Res.70 (3), 1173–1183. 10.1158/0008-5472.CAN-09-2886

  • CrossRef
  • Google Scholar

• 39

  KasimovaM. A.TarekM.ShaytanA. K.ShaitanK. V.DelemotteL. (2014). Voltage-gated ion channel modulation by lipids: Insights from molecular dynamics simulations. Biochim. Biophys. Acta1838 (5), 1322–1331. 10.1016/j.bbamem.2014.01.024

  • CrossRef
  • Google Scholar

• 40

  KhanA. A.SoloskiM. J.SharpA. H.SchillingG.SabatiniD. M.LiS. H.et al (1996). Lymphocyte apoptosis: Mediation by increased type 3 inositol 1,4,5-trisphosphate receptor. Science273 (5274), 503–507. 10.1126/science.273.5274.503

  • CrossRef
  • Google Scholar

• 41

  KoubaS.OuldamerL.GarciaC.FontaineD.ChantomeA.VandierC.et al (2019). Lipid metabolism and Calcium signaling in epithelial ovarian cancer. Cell Calcium81, 38–50. 10.1016/j.ceca.2019.06.002

  • CrossRef
  • Google Scholar

• 42

  KuchayS.GiorgiC.SimoneschiD.PaganJ.MissiroliS.SarafA.et al (2017). PTEN counteracts FBXL2 to promote IP3R3-and Ca(2+)-mediated apoptosis limiting tumour growth. Nature546 (7659), 554–558. 10.1038/nature22965

  • CrossRef
  • Google Scholar

• 43

  LangF.StournarasC. (2014). Ion channels in cancer: Future perspectives and clinical potential. Philos. Trans. R. Soc. Lond B Biol. Sci.369 (1638), 20130108. 10.1098/rstb.2013.0108

  • CrossRef
  • Google Scholar

• 44

  LangfelderP.HorvathS. (2008). Wgcna: an R package for weighted correlation network analysis. BMC Bioinforma.9, 559. 10.1186/1471-2105-9-559

  • CrossRef
  • Google Scholar

• 45

  LangfelderP.ZhangB.HorvathS. (2008). Defining clusters from a hierarchical cluster tree: The dynamic tree cut package for R. Bioinformatics24 (5), 719–720. 10.1093/bioinformatics/btm563

  • CrossRef
  • Google Scholar

• 46

  LeanzaL.ManagoA.ZorattiM.GulbinsE.SzaboI. (2016). Pharmacological targeting of ion channels for cancer therapy: In vivo evidences. Biochim. Biophys. Acta1863 (6), 1385–1397. 10.1016/j.bbamcr.2015.11.032

  • CrossRef
  • Google Scholar

• 47

  LiJ.LiQ.SuZ.SunQ.ZhaoY.FengT.et al (2020). Lipid metabolism gene-wide profile and survival signature of lung adenocarcinoma. Lipids Health Dis.19 (1), 222. 10.1186/s12944-020-01390-9

  • CrossRef
  • Google Scholar

• 48

  LiX.ZhangQ.FanK.LiB.LiH.QiH.et al (2016). Overexpression of TRPV3 correlates with tumor progression in non-small cell lung cancer. Int. J. Mol. Sci.17 (4), 437. 10.3390/ijms17040437

  • CrossRef
  • Google Scholar

• 49

  LitanA.LanghansS. A. (2015). Cancer as a channelopathy: Ion channels and pumps in tumor development and progression. Front. Cell Neurosci.9, 86. 10.3389/fncel.2015.00086

  • CrossRef
  • Google Scholar

• 50

  LiuX.ChangY.ReinhartP. H.SontheimerH.ChangY. (2002). Cloning and characterization of glioma BK, a novel BK channel isoform highly expressed in human glioma cells. J. Neurosci.22 (5), 1840–1849. 10.1523/JNEUROSCI.22-05-01840.2002

  • CrossRef
  • Google Scholar

• 51

  LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.15 (12), 550. 10.1186/s13059-014-0550-8

  • CrossRef
  • Google Scholar

• 52

  LutzC. S.CornettA. L. (2013). Regulation of genes in the arachidonic acid metabolic pathway by RNA processing and RNA-mediated mechanisms. Wiley Interdiscip. Rev. RNA4 (5), 593–605. 10.1002/wrna.1183

  • CrossRef
  • Google Scholar

• 53

  MandalS.TejaswiT.JanivaraR.SrikrishnanS.ThakurP.SahooS.et al (2021). Transcriptomic-based quantification of the epithelial-hybrid-mesenchymal spectrum across biological contexts. Biomolecules12 (1), 29. 10.3390/biom12010029

  • CrossRef
  • Google Scholar

• 54

  MengQ.HuX.ZhaoX.KongX.MengY. M.ChenY.et al (2021). A circular network of coregulated sphingolipids dictates lung cancer growth and progression. EBioMedicine66, 103301. 10.1016/j.ebiom.2021.103301

  • CrossRef
  • Google Scholar

• 55

  MinkeB.CookB. (2002). TRP channel proteins and signal transduction. Physiol. Rev.82 (2), 429–472. 10.1152/physrev.00001.2002

  • CrossRef
  • Google Scholar

• 56

  MoskotM.Jakobkiewicz-BaneckaJ.SmolinskaE.BaneckiB.WegrzynG.Gabig-CiminskaM. (2015). Activities of genes controlling sphingolipid metabolism in human fibroblasts treated with flavonoids. Metab. Brain Dis.30 (5), 1257–1267. 10.1007/s11011-015-9705-x

  • CrossRef
  • Google Scholar

• 57

  NorthcottP. A.DubucA. M.PfisterS.TaylorM. D. (2012). Molecular subgroups of medulloblastoma. Expert Rev. Neurother.12 (7), 871–884. 10.1586/ern.12.66

  • CrossRef
  • Google Scholar

• 58

  ParthasarathiK. T. S.MandalS.SinghS.GundimedaS.JollyM. K.PandeyA.et al (2022). In silico analysis of ion channels and their correlation with epithelial to mesenchymal transition in breast cancer. Cancers (Basel)14 (6), 1444. 10.3390/cancers14061444

  • CrossRef
  • Google Scholar

• 59

  RahB.RatherR. A.BhatG. R.BabaA. B.MushtaqI.FarooqM.et al (2022). JAK/STAT signaling: Molecular targets, therapeutic opportunities, and limitations of targeted inhibitions in solid malignancies. Front. Pharmacol.13, 821344. 10.3389/fphar.2022.821344

  • CrossRef
  • Google Scholar

• 60

  RussoA.Biselli-ChicoteP. M.Kawasaki-OyamaR. S.Castanhole-NunesM. M. U.ManigliaJ. V.de Santi NetoD.et al (2018). Differential expression of prostaglandin I2 synthase associated with arachidonic acid pathway in the oral squamous cell carcinoma. J. Oncol.2018, 6301980. 10.1155/2018/6301980

  • CrossRef
  • Google Scholar

• 61

  ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res.13 (11), 2498–2504. 10.1101/gr.1239303

  • CrossRef
  • Google Scholar

• 62

  SharmaJ.BalakrishnanL.KaushikS.KashyapM. K. (2020). Editorial: Multi-Omics approaches to study signaling pathways. Front. Bioeng. Biotechnol.8, 829. 10.3389/fbioe.2020.00829

  • CrossRef
  • Google Scholar

• 63

  SharmaJ.DebB.GeorgeI. A.KapilS.CoralK.KakkarN.et al (2019). Somatic mutations profile of a young patient with metastatic urothelial carcinoma reveals mutations in genes involved in ion channels. Front. Oncol.9, 435. 10.3389/fonc.2019.00435

  • CrossRef
  • Google Scholar

• 64

  ShiN.YuH.ChenT. (2019). Inhibition of esophageal cancer growth through the suppression of PI3K/AKT/mTOR signaling pathway. Onco Targets Ther.12, 7637–7647. 10.2147/OTT.S205457

  • CrossRef
  • Google Scholar

• 65

  ShibaoK.FiedlerM. J.NagataJ.MinagawaN.HirataK.NakayamaY.et al (2010). The type III inositol 1,4,5-trisphosphate receptor is associated with aggressiveness of colorectal carcinoma. Cell Calcium48 (6), 315–323. 10.1016/j.ceca.2010.09.005

  • CrossRef
  • Google Scholar

• 66

  ShiozakiA.KudouM.IchikawaD.FujiwaraH.ShimizuH.IshimotoT.et al (2018). Esophageal cancer stem cells are suppressed by tranilast, a TRPV2 channel inhibitor. J. Gastroenterol.53 (2), 197–207. 10.1007/s00535-017-1338-x

  • CrossRef
  • Google Scholar

• 67

  ShiozakiA.MarunakaY.OtsujiE. (2021). Roles of ion and water channels in the cell death and survival of upper gastrointestinal tract cancers. Front. Cell Dev. Biol.9, 616933. 10.3389/fcell.2021.616933

  • CrossRef
  • Google Scholar

• 68

  SonnweberT.PizziniA.NairzM.WeissG.TancevskiI. (2018). Arachidonic acid metabolites in cardiovascular and metabolic diseases. Int. J. Mol. Sci.19 (11), 3285. 10.3390/ijms19113285

  • CrossRef
  • Google Scholar

• 69

  SugaseT. T.SeradaS.NakatsukaR.FujimotoM.OhkawaraT.HaraH.et al (2017). Suppressor of cytokine signaling-1 gene therapy induces potent antitumor effect in patient-derived esophageal squamous cell carcinoma xenograft mice. Int. J. cancer140, 2608–2621. 10.1002/ijc.30666

  • CrossRef
  • Google Scholar

• 70

  SzklarczykD.GableA. L.NastouK. C.LyonD.KirschR.PyysaloS.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 (1), D605–D612. 10.1093/nar/gkaa1074

  • CrossRef
  • Google Scholar

• 71

  UeasilamongkolP.KhamphayaT.GuerraM. T.RodriguesM. A.GomesD. A.KongY.et al (2020). Type 3 inositol 1,4,5-trisphosphate receptor is increased and enhances malignant properties in cholangiocarcinoma. Hepatology71 (2), 583–599. 10.1002/hep.30839

  • CrossRef
  • Google Scholar

• 72

  von MeringC.HuynenM.JaeggiD.SchmidtS.BorkP.SnelB. (2003). String: A database of predicted functional associations between proteins. Nucleic Acids Res.31 (1), 258–261. 10.1093/nar/gkg034

  • CrossRef
  • Google Scholar

• 73

  WangH.TianL.LiuJ.GoldsteinA.BadoI.ZhangW.et al (2018). The osteogenic niche is a calcium reservoir of bone micrometastases and confers unexpected therapeutic vulnerability. Cancer Cell34 (5), 823–839.e7. 10.1016/j.ccell.2018.10.002

  • CrossRef
  • Google Scholar

• 74

  WangM. M.GuoC. H.LiF. L.XuR. P.LiuZ.PanY. Q.et al (2021a). Family history of esophageal cancer modifies the association of serum lipids and malignant esophageal lesions: A nested case-control study from the "endoscopic screening for esophageal cancer in China" trial. Chin. Med. J. Engl.134 (9), 1079–1086. 10.1097/CM9.0000000000001432

  • CrossRef
  • Google Scholar

• 75

  WangM. M.XiongY.ZhaoL.GuK.LiY.ZhaoF.et al (2021b). UCSCXenaShiny: An R/CRAN package for interactive analysis of UCSC Xena data. Bioinformatics38, 527–529. 10.1093/bioinformatics/btab561

  • CrossRef
  • Google Scholar

• 76

  WangW.HeS.ZhangR.PengJ.GuoD.ZhangJ.et al (2020). ALDH1A1 maintains the cancer stem-like cells properties of esophageal squamous cell carcinoma by activating the AKT signal pathway and interacting with beta-catenin. Biomed. Pharmacother.125, 109940. 10.1016/j.biopha.2020.109940

  • CrossRef
  • Google Scholar

• 77

  WeiT.SimkoV.LevyM.XieY.JinY.ZemlaJ. (2017). Package ‘corrplot’. Statistician. 56, 316–324.

  • Google Scholar

• 78

  WhelerJ. J.JankuF.FalchookG. S.JacksonT. L.FuS.NaingA.et al (2014). Phase I study of anti-VEGF monoclonal antibody bevacizumab and histone deacetylase inhibitor valproic acid in patients with advanced cancers. Cancer Chemother. Pharmacol.73 (3), 495–501. 10.1007/s00280-014-2384-1

  • CrossRef
  • Google Scholar

• 79

  WickhamH. (2016). ggplot2: Elegant graphics for data analysis. Berlin, Germany: Springer.

  • Google Scholar

• 80

  WunderlichC. M.HovelmeyerN.WunderlichF. T. (2013). Mechanisms of chronic JAK-STAT3-SOCS3 signaling in obesity. JAKSTAT2 (2), e23878. 10.4161/jkst.23878

  • CrossRef
  • Google Scholar

• 81

  YangY. C.ChangT. Y.ChenT. C.LinW. S.ChangS. C.LeeY. J. (2017). ITPR3 gene haplotype is associated with cervical squamous cell carcinoma risk in Taiwanese women. Oncotarget8 (6), 10085–10090. 10.18632/oncotarget.14341

  • CrossRef
  • Google Scholar

• 82

  ZhangB.HorvathS. (2005). A general framework for weighted gene co-expression network analysis. Stat. Appl. Genet. Mol. Biol.4, 17. 10.2202/1544-6115.1128

  • CrossRef
  • Google Scholar

• 83

  ZhangG.WangX.XueQ. (2018). Combined targeted ion channel therapy: Can it be an alternative choice for esophageal cancer patients?Med. Hypotheses117, 59–62. 10.1016/j.mehy.2018.06.010

  • CrossRef
  • Google Scholar

• 84

  ZhangZ.HeQ.FuS.ZhengZ. (2017). Estrogen receptors in regulating cell proliferation of esophageal squamous cell carcinoma: Involvement of intracellular Ca(2+) signaling. Pathol. Oncol. Res.23 (2), 329–334. 10.1007/s12253-016-0105-2

  • CrossRef
  • Google Scholar

• 85

  ZhaoF.GaoS.QinX.NiuR.LiZ.WangC.et al (2022). Comprehensive analysis of TRP channel-related genes for estimating the immune microenvironment, prognosis, and therapeutic effect in patients with esophageal squamous cell carcinoma. Front. Cell Dev. Biol.10, 820870. 10.3389/fcell.2022.820870

  • CrossRef
  • Google Scholar

• 86

  ZhengW.BaiX.ZhouY.YuL.JiD.ZhengY.et al (2022). Transcriptional ITPR3 as potential targets and biomarkers for human pancreatic cancer. Aging (Albany NY)14 (10), 4425–4444. 10.18632/aging.204080

  • CrossRef
  • Google Scholar