Next Article in Journal
A Novel m7G-Related Gene Signature Predicts the Prognosis of Colon Cancer
Next Article in Special Issue
Blood Plasma Metabolome Profiling at Different Stages of Renal Cell Carcinoma
Previous Article in Journal
Lenvatinib as First-Line Treatment for Unresectable Hepatocellular Carcinoma: A Systematic Review and Meta-Analysis
Previous Article in Special Issue
Integrative Analysis Identifies TCIRG1 as a Potential Prognostic and Immunotherapy-Relevant Biomarker Associated with Malignant Cell Migration in Clear Cell Renal Cell Carcinoma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LRRC8A Is a Promising Prognostic Biomarker and Therapeutic Target for Pancreatic Adenocarcinoma

1
Division of Cancer Biology, Laboratory Animal Center, The Fourth Military Medical University, Xi’an 710032, China
2
Department of Geriatrics, Xijing Hospital, The Fourth Military Medical University, Xi’an 710032, China
3
School of Basic Medical Sciences, Medical College of Yan’an University, Yan’an 716000, China
4
Laboratory Animal Center, Guangzhou University of Chinese Medicine, Guangzhou 510405, China
5
Laboratory Animal Center, General Hospital of Northern Theater Command, Shenyang 110016, China
6
College of Basic Medicine, Shaanxi University of Chinese Medicine, Xianyang 712046, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2022, 14(22), 5526; https://doi.org/10.3390/cancers14225526
Submission received: 3 October 2022 / Revised: 18 October 2022 / Accepted: 3 November 2022 / Published: 10 November 2022
(This article belongs to the Collection The Biomarkers for the Diagnosis and Prognosis in Cancer)

Abstract

:

Simple Summary

Pancreatic adenocarcinoma (PAAD) is a highly malignant tumor with a poor prognosis and lack of effective biomarkers. We demonstrated, using bioinformatics analyses, that the leucine-rich repeat-containing 8A (LRRC8A) can affect the prognosis of patients with PAAD. In addition, experiments in vitro showed shRNA of LRRC8A can significantly affect the hallmarks of PAAD, including cell proliferation, cell migration, drug resistance, and immune infiltration. Collectively, our data suggest that LRRC8A could be a promising prognostic biomarker and therapeutic target for PAAD.

Abstract

Pancreatic adenocarcinoma (PAAD) is a highly malignant tumor of the digestive system with increasing morbidity and mortality. The lack of sensitive and reliable biomarkers is one of the main reasons for the poor prognosis. Volume-regulated anion channels (VRAC), which are ubiquitously expressed in the vertebrate cell membrane, are composed of leucine-rich repeat-containing 8A (LRRC8A) and four other homologous family members (LRRC8B–E). VRAC heterogeneous complex is implicated in each of the six “hallmarks of cancer” and represents a novel therapeutic target for cancer. In this study, LRRC8A was speculated to be a promising prognostic biomarker and therapeutic target for PAAD based on a series of bioinformatics analyses. Additional cell experiments and immunohistochemical assays demonstrated that LRRC8A can affect the prognosis of PAAD and is correlated to cell proliferation, cell migration, drug resistance, and immune infiltration. Functional analysis indicated that LRRC8A influences the progression and prognosis of patients with PAAD by the regulation of CD8+ T cells immune infiltration. Taken together, these results can help in the design of new therapeutic drugs for patients with PAAD.

1. Introduction

Pancreatic adenocarcinoma (PAAD) is a deadly malignant tumor, with a 5-year survival rate of only 2–9%. It is the fourth leading cause of cancer-associated mortality [1,2]. Most patients with PAAD remain asymptomatic during the initial stage of cancer development until they progress to advanced metastasis, wherein tumor cells become highly invasive [3]. The lack of sensitive and reliable biomarkers is one of the main reasons for its poor prognosis. Although the application of carbohydrate antigen (CA) 19-9, which is the primary biomarker used at present for PAAD has witnessed incremental progress in the diagnosis of PAAD, it has a limitation that cannot be overlooked, wherein CA19-9 showed false-positive results in patients with chronic pancreatitis, cholangitis, and other diseases [4]. Several novel biomarkers, including HOX antisense intergenic RNA (a type of long non-coding RNA), C4b binding protein alpha, micro-RNA, and exosomes have been identified as potential PAAD biomarkers, but these markers cannot guarantee high reliability, sensitivity, and specificity owing to some inherent limitations, such as inadequate quantities of micro RNA in various biofluids [5,6,7]. Another challenge in the treatment of PAAD is the complex tumor microenvironment, including the presence of tumor-associated macrophages and regulatory T cells within the stroma that create a barrier for conventional therapies and immunotherapies [8]. Although numerous clinical trials targeting the stroma have been conducted, most of them did not provide the expected efficacy based on preclinical studies [9,10]. Therefore, identifying novel biomarkers to predict the prognosis of PAAD and developing effective therapeutic targets are required for the treatment of PAAD.
Volume-regulated anion channels (VRACs), also called volume-sensitive outwardly rectifying anion channels, are ubiquitously expressed in the membrane of vertebrate cells. Swelling-activated chloride currents are activated when VRACs are opened. Leucine-rich repeat-containing 8A (LRRC8A) and four other homologous family members (LRRC8B–E) assemble into VRAC complexes. Since these proteins can exist in a soluble form and as an integral membrane protein, cyst fluid sampling and testing for LRRC8A and its homologous family proteins are feasible using mass spectroscopy. This would provide a potential opportunity for early detection and management [11]. Moreover, as the main chloride channel type, VRAC is implicated in each of the six “hallmarks of cancer” and represents an emerging therapeutic target for cancer [12,13,14]. A dysregulated expression or dysfunction of VRAC can be detected in many cancer types, indicating a correlation between ion channels and tumorigenesis [15]. LRRC8A downregulation can reportedly increase the sensitivity of glioblastoma to temozolomide and carmustine [16]. LRRC8A can regulate the proliferation and metastasis of colorectal cancer, hepatocellular carcinoma, and nasopharyngeal carcinoma [16,17,18,19]. Recent studies also supported the finding that LRRC8D and LRRC8A/LRRC8E play an integral role in the regulation of multidrug resistance and tumor immunotherapy, respectively [20]. Previous studies confirmed the role of calcium-activated chloride channel regulator 1 as a potential suppressor and independent factor in PAAD [21]; however, there are no data regarding the prognostic role of VRAC in PAAD.
In this study, a series of bioinformatics analyses based on several large public databases were conducted, and LRRC8A was identified as a possible potent therapeutic target and prognostic biomarker for PAAD based on cell experiments. Additional functional analysis, immune infiltration analysis, and immunohistochemical assay indicated that LRRC8A influences the progression and prognosis of PAAD by the regulation of CD8+ T cells immune infiltration.

2. Materials and Methods

2.1. Pan-Cancer Analysis of LRRC8A-8E Expression by ONCOMINE

The mRNA expression of LRRC8A, LRRC8B, LRRC8C, LRRC8D, and LRRC8E in different human cancer types was analyzed using ONCOMINE (www.oncomine.org (accessed on 11 June 2021)). A p-value of 0.05 and fold change of 2 were set as significance thresholds.

2.2. LRRC8A–8E Expression in PAAD by GEPIA and UALCAN

The transcriptional levels of LRRC8A–8E in PAAD were thoroughly and comprehensively analyzed using GEPIA (http://gepia.cancer-pku.cn/ (accessed on 27 July 2021)) and UALCAN (http://ualcan.path.uab.edu/analysis.html (accessed on 27 July 2021)). The expression data of tumor samples (n = 179) and normal samples (n = 171) in GEPIA were acquired from the Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression. Each box plot of LRRC8A–8E in PAAD was designed using the “Expression DIY” module of GEPIA. The expression levels of LRRC8A–8E based on individual cancer stages in PAAD were analyzed using the TCGA dataset in UALCAN. p < 0.05 was set as a significance threshold.

2.3. Survival Analysis Using Kaplan–Meier Plotter and GEPIA

The correlation between LRRC8A–8E expression and clinical outcomes of PAAD was assessed using a Kaplan–Meier Plotter (https://kmplot.com/analysis/ (accessed on 7 January 2022)) and the module of survival analysis in GEPIA. The Kaplan–Meier plotter can assess the overall survival (OS) and relapse-free survival (RFS) of patients with PAAD (n = 177), whereas the OS and disease-free survival (DFS) analysis of PAAD can be obtained using GEPIA (n = 178).

2.4. Functional Analysis of LRRC8A–8E in Patients with PAAD

The co-expression and network module of LRRC8A–E in patients with PAAD were analyzed using cBioPortal (www.cbioportal.org (accessed on 27 July 2021)), which included 11 studies with 1347 PAAD samples. The Gene Ontology (GO) enrichment heatmap and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the top 50 genes positively associated with LRRC8A were performed using OmicShare tools, which is a free online platform for data analysis (https://www.omicshare.com (accessed on 25 April 2022)). The protein–protein interaction (PP—network of LRRC8A–E was analyzed using the Search Tool for Retrieval of Interacting Genes/Proteins (https://string-db.org/ (accessed on 23 July 2021)) and GeneMANIA (http://www.genemania.org (accessed on 27 July 2021)).

2.5. Correlations between LRRC8A–E Expression and Immune Infiltration Using Tumor Immune Estimation Resource (TIMER)

GEPIA databases and TIMER (http://cistrome.org/TIMER/ (accessed on 28 July 2021)) were used to determine the relationship between LRRC8A–E expression and immune infiltration. The TIMER database contains 10,897 samples across 32 cancer types from TCGA to allow the evaluation of the abundance of immune infiltration.

2.6. Cell Lines

PAAD cell lines PANC1 (CRL-1469™), AsPC-1 (CRL-1682™), and BxPC-3 (CRL-1687) were obtained from the American Type Culture Collection (USA), AsPC-1GEM-R cell line (gemcitabine-resistant AsPC-1) was established in our laboratory, which can grow in the medium containing 18 μg/mL gemcitabine for more than 4 days. All cell lines were cultured in Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum. Mycoplasma was detected and determined to be negative.

2.7. Immunochemistry (IHC) Analysis

PAAD tissues were obtained from the Department of Pathology of Xijing Hospital, and patient consent was reviewed and approved by the Xijing Hospital Clinical Research Ethics Board (certificate number KY2015432). Based on standard protocols [22], IHC analysis was performed using antibodies such as LRRC8A (Biobyt, Cambridge, United Kingdom; orb185053) and CD8 (Proteintech, Rosemont, IL, USA; 66868-1-Ig). The area of positive cells, average optical density values, and integrated optical density values were calculated using Image-Pro Plus analysis software 2020 to evaluate the immunohistochemical results. The R software (version 4.1.0) was used to linearly fit the immunohistochemical quantification data and the function Cook’s distance was used for data visualization.

2.8. Real-Time Polymerase Chain Reaction (PCR) Assays

Real-time PCR was performed using a PrimeScript™ RT Reagent Kit and TB Green® Fast qPCR Mix purchased from Takara (Beijing, China) using a Step One Plus real-time PCR Detection System (Applied Biosystems, Thermo Fisher, Waltham, MA, USA). Relative mRNA expression levels were calculated using the 2−ΔΔCT method, and β-actin mRNA expression was set as the internal reference. The sequences of primers used for PCR are shown in Table 1.

2.9. Western Blotting Analysis

Western blotting analyses were performed using standard protocols [23]. Antibodies against the following proteins were used: LRRC8A (Biobyt, orb185053) and β-actin (Engibody, Dover, DE, USA; AT0001). A grayscale value of western blotting was obtained using the Image J software 2020 and analyzed using the GraphPad Prism 6 software.

2.10. LRRC8A Knockdown

The knockdown of LRRC8A was accomplished by employing short hairpin (shRNA) lentivirus (sense: CCGUCUACUACGUGCACAATT, antisense: UUGUGCACGUAGUAGACGGTT) (Hanbio, Shanghai, China) infection as reported previously [24], and this was verified by real-time PCR and western blotting analysis.

2.11. Cell Viability Assay

Cells were seeded in a 96-well plate and pretreated with 10 μM 4-(2-Butyl-6, 7-dichloro-2-cyclopentyl-indan-1-on-5-yl) oxobutyric acid (DCPIB), a kind of specific VRAC inhibitor. Then, the cells were treated with 100 μM gemcitabine, and cell viability was detected at different times using a Cell Counting Kit-8 (CCK-8) (Solarbio, Beijing, China; CA1210).

2.12. Wound Healing Assay

AsPC-1 and BxPC-3 cells were seeded in a six-well plate and infected with ShRNA lentivirus of LRRC8A. When the cells reached approximately 80% confluence, a pipette tip was used to create a straight line in the middle of the plate. Then, the shed cells were rinsed using phosphate-buffered saline, and a serum-free medium was added. Images of the straight line at 0 h and 48 h were captured, and the relative migration distance was analyzed to elucidate the scribing process of the cells. AsPC-1 and BxPC-3 cells infected with scramble shRNA lentivirus were used as the negative control.

2.13. Transwell Assay

AsPC-1 and BxPC-3 cells with knocked-down LRRC8A were added to the chamber and cultured for 48 h. Then, the cells were treated with 4% paraformaldehyde and stained with crystal violet. Images of cells that passed through the chamber were captured under a microscope.

2.14. Statistical Analysis

All data are presented as means ± standard error of the mean and analyzed using the GraphPad Prism 6 software. The significance was analyzed using Student’s t-tests. p < 0.05 was considered statistically significant.

3. Results

3.1. LRRC8A Is Highly Expressed in PAAD in Pan-Cancer Analysis

The expression of LRRC8A and four other homologous family members (LRRC8B–E) in pan-cancer was retrieved using the ONCOMINE database (Figure 1a). Results revealed that LRRC8A expression was higher in PAAD samples than in normal samples. To further evaluate its expression in PAAD, GEPIA, and UALCAN were used. As shown in Figure 1b, the transcripts per million expression level of LRRC8A in PAAD samples was seven folds higher than that in normal samples (46.52 in tumor samples vs. 6.59 in normal samples). Subsequently, an in-depth analysis by UALCAN demonstrated a significant correlation between LRRC8A expression (p * < 0.05) and the pathological stage of patients with PAAD (Figure 1c). The expression of LRRC8B–8E in different pathological stages of PAAD is shown in Supplementary Figure S1a,b.

3.2. LRRC8A Expression Is Associated with Poor Prognosis in Patients with PAAD

To evaluate the value of differentially expressed LRRC8A–E in the progression of PAAD, its correlation with clinical outcomes was assessed using GEPIA (n = 178) and Kaplan–Meier plotter (n = 177). Using GEPIA, the DFS (LOGRANK test, p = 0.0052; hazards ratio (HR), 1.9 (high); p (HR) = 0.006) and OS rates (logrank test, p = 0.017; HR, 1.7 (high); p (HR) = 0.018) of patients with PAAD with a high LRRC8A expression were significantly worse than those with a low LRRC8A expression (Figure 1d,e). Similarly, the RFS of patients with PAAD using a Kaplan–Meier plotter also indicated that LRRC8A was a detrimental prognostic factor of PAAD (RFS: HR, 4.13, 95% CI, 1.6–10.69; logrank p = 0.0016) (Figure 1f). The correlation of LRRC8B–8E with clinical outcomes of PAAD is shown in Supplementary Figure S1.

3.3. Identification of the Role of LRRC8A in PAAD In Vitro

As mentioned above, uncontrolled proliferation, frequent migration, and multidrug resistance are the fundamental hallmarks of cancer cells. A gemcitabine-resistant AsPC-1GEMR cell line was established to investigate whether LRRC8A was involved in the drug resistance of PAAD (Figure 2a). The CCK-8 assay results revealed that the knockdown of LRRC8A using shRNA or DCPIB, the specific VRAC inhibitor can inhibit the proliferation of PAAD cells and increase their sensitivity to gemcitabine (Figure 2b, Supplementary Figure S1d). Additionally, qRT-PCR and western blot results showed LRRC8A decreased after AsPC-1GEMR cells were pretreated with DCPIB (Supplementary Figure S1d) and the knockdown effect of shLRRC8A in PAAD cell lines was showed in Supplementary Figure S1e.
To evaluate the potential role of LRRC8A in PAAD, qRT-PCR and western blotting were used to detect the expression of LRRC8A in PAAD cell lines AsPC-1, PaNC-1, and BxPC-3. The sequence of cell migration and invasion capacity, angiogenic potential, and tumorigenicity was reportedly BxPC-3 > AsPC-1 > PaNC-1 [25]. A similar mRNA and protein expression gradient of LRRC8A in PAAD cell lines was also observed in Figure 2c,d, respectively. The immunohistochemical analysis of pathological sections obtained from patients with PAAD demonstrated that a high LRRC8A expression was correlated to the tumor grade of PAAD (Figure 2e). To investigate whether LRRC8A affects the migration and invasion abilities of PAAD cells, a wound healing assay with or without LRRC8A knockdown was conducted using AsPC-1GEMR and BxPC-3 cells as showed in Figure 2f, and a transwell invasion assay using AsPC-1 and BxPC-3 cells with or without LRRC8A knockdown were performed in Figure 2g. Results showed a substantial inhibitory effect of LRRC8A knockdown on PAAD cell migration and invasion (Figure 2f,g). Further investigation showed LRRC8A knockdown resulted in the decreased protein expression of EMT markers, such as MMP2 in PAAD cell lines (Figure 2h).
In conclusion, the results suggest that LRRC8A plays a significant role in the tumorigenesis and progression of PAAD.

3.4. Co-Expression, Gene Network, and Enrichment Analysis of LRRC8A in PAAD

The above-mentioned findings suggest that LRRC8A is a promising prognostic biomarker of PAAD. Therefore, it would be crucial to explore the potential relevance and underlying mechanisms of LRRC8A expression in patients with PAAD. Next, the LRRC8A co-expression and its network in patients with PAAD were explored, and the potential co-expression of LRRC8A genes was analyzed using eBioPortal. In PAAD database (TCGA, Firehose Legacy), the top 50 positively correlated genes of LRRC8A were imported to OmicShare to analyze the GO enrichment and KEGG pathways (Figure 3). Figure 3a shows that biological processes, including skin development, hemidesmosome assembly, epidermis development, tissue development, keratinocyte differentiation, epidermal cell differentiation, and cell junction, were associated with the tumorigenesis and progression of PAAD. Cell adhesion molecule binding, cadherin binding, protein binding involved in cell–cell adhesion, protein binding involved in cell adhesion, cadherin binding involved in cell–cell adhesion, and death receptor binding were the most highly enriched molecular functions, whereas adherence junction, anchoring junction, cell junction, membrane-bound vesicles, vesicles, extracellular region, and cell–substrate junction were the most highly enriched cellular components. Based on the KEGG pathway analysis, the phosphoinositide 3-kinase (PI3K)-Akt signaling pathway, focal adhesion, arrhythmogenic right ventricular cardiomyopathy, extracellular membrane–receptor interaction, and mitogen-activated protein kinase signaling pathway were significantly associated with the tumorigenesis and progression of PAAD (Figure 3b). Moreover, given that LRRC8A co-localizes with LRRC8B–E in cells, it is necessary to analyze the PPI network of LRRC8A–E (Figure 3c), the gene–gene interaction network of LRRC8A–E, and the top 50 positively correlated genes of LRRC8A (Figure 3d).

3.5. Correlation of LRRC8A with Immune Infiltration

Interestingly, the co-expression analysis of LRRC8A in PAAD samples showed that purinergic receptor P2Y2 (P2RY2) and ephrin-A receptor 2 (EPHA2) showed a positive correlation with LRRC8A in most database of PAAD, including those of TCGA (PanCancer Atlas and Firehose Legacy), QCMG (Nature 2016) [26], and TIMER (Figure 4). P2RY2 was reported to enhance cancer cell glycolysis by PI3K/AKT-mammalian target of rapamycin (mTOR) signaling, resulting in PAAD progression [27]. However, EPHA2 was reportedly involved in multiple immunity regulation pathways of PAAD [28,29]. The mRNA expression level of EPHA2 decreased considerably when LRRC8A was knocked down in AsPC-1GEM-R and BxPC-3 cells (Figure 5a). Recent studies showed that LRRC8A has a novel role in immune regulation [20,30]; therefore, whether its expression was correlated with immune infiltration levels using TIMER was investigated (Figure 5b). Results showed that LRRC8A expression was significantly correlated with CD8+ T cells (p = 1.83 × 10−3), neutrophils (p = 1.58 × 10−2), and dendritic cells (p = 7.88 × 10−3). Similarly, a correlation analysis of LRRC8A and CD8+ T cells in PAAD (179 samples) using the TISIDB database also showed that LRRC8A expression was significantly correlated with CD8+ T cells (ρ = 0.181, p = 0.0155) (Figure 5c). Further immunohistochemical results of pathological sections obtained from patients with PAAD showed that LRRC8A expression was correlated with CD8 expression (Figure 5d–f), with a correlation coefficient (R2) of 0.83. These findings suggest that LRRC8A affects the patient clinical outcomes by interacting with immune infiltration, especially CD8+ T cells in PAAD. It should be noted that LRRC8C was significantly correlated with B cells (p = 1.67 × 10−7), CD8+ T cells (p = 2.22 × 10−16), CD4+ T cells (p = 8.02 × 10−5), macrophages (p = 7.21 × 10−26), neutrophils (p = 5.30 × 10−20), and dendritic cells (p = 3.14 × 10−27) (Supplementary Figure S2).

4. Discussion

A limitless replicative potential, sustained angiogenesis, evasion of apoptosis, self-sufficiency in growth signals, and insensitivity to anti-growth signals, tissue invasion, and metastasis are classical hallmarks of cancer cells [12,13]. VRAC, one of the most important chloride channels located at the cell membrane and constituted by LRRC8A and four other homologous family members LRRC8B–E, has been implicated in each hallmark of cancer [31]. By using non-specific pharmacological VRAC blockers such as DIDS, DCPIB, and NPPB, VRAC was reported to be involved in the cell cycle regulation in human ovarian cancer, cervical cancer, nasopharyngeal carcinoma, small-cell lung cancer and glioblastoma [14]. Moreover, VRAC activity has also been implicated in the migration and multi-drug resistance in multiple cancers. In particular, VRAC current was record in pancreatic duct adenocarcinoma cell lines [32]. Hence, VRAC represents an emerging therapeutic target of PAAD. As the component of VRAC, LRRC8A–E was reported to be related to cancer cell proliferation, metastasis, and multidrug resistance in some studies, suggesting that LRRC8A–E may be the specific molecular regulated tumor progression [14]. However, the biological function and prognostic value of LRRC8A–E in PAAD have not been well-characterized.
In this study, the expression of LRRC8A–E in pan-cancer was explored using independent datasets in ONCOMINE and TCGA data in GEPIA. LRRC8A expression was higher in PAAD tissues than that in normal tissues, and it increased as PAAD progressed. Cell experiments in vitro also confirmed the high expression of LRRC8A in PAAD cells. The immunohistochemical results of pathological sections obtained from patients with PAAD also demonstrated that the high expression of LRRC8A was associated with a higher tumor grade of PAAD. Moreover, a high LRRC8A expression in patients with PAAD was significantly associated with worse clinical outcomes based on the results obtained using GEPIA and Kaplan–Meier plotter. Fortunately, two research groups almost simultaneously reported that LRRC8 heteromers were essential components of VRAC, and LRRC8A is an indispensable VRAC subunit. Therefore, the critical role and molecular mechanisms of LRRC8A in cancer management need further clarification [33,34].
Because proliferation, migration, invasion, and drug resistance are the main hallmarks of cancer cells, whether LRRC8A knockdown can inhibit the migration, invasion, and drug resistance of PAAD cells in vitro was further investigated. The wound healing assay, transwell invasion assay, and drug sensitivity assay demonstrated the inhibitory effects of LRRC8A knockdown, suggesting its potential as a therapeutic target in PAAD. Then, the functions of LRRC8A were investigated using the co-expression module of eBioPortal and GO enrichment analysis. The KEGG pathway enrichment analysis was conducted to explore the probable mechanisms in PAAD. The neighboring genes of LRRC8A were mainly enriched in the PI3K-AKT pathway and focal adhesion pathway. A previous study found that PI3K/AKT was a downstream signal responsible for the regulatory effects of LRRC8A in cancer cell proliferation [18]. The LRRC8A-PI3K/AKT-mTOR pathways can also regulate the migration of esophageal squamous cell carcinoma [17].
More importantly, EPHA2 showed a positive correlation with LRRC8A in most studies of PAAD, and its knockdown resulted in a decrease in EPHA2 expression, which was demonstrated to be involved in anti-tumor immunity by regulating prostaglandin-endoperoxide synthase. Moreover, LRRC8A is essential for T cell development [35], and it has been recently reported to play an integral role in tumor immunotherapy [20,30]. Heteromeric VRAC channels composed of LRRC8A and either LRRC8C or LRRC8E can transport cGAMP into bystander cells to activate innate immune cGAS–cGAMP–STING pathway in vascular epithelial cells and murine macrophages. The immune infiltration findings in TIMER also suggested that LRRC8A affects patient clinical outcomes by interacting with immune infiltration in PAAD. However, immune infiltration analysis also showed that LRRC8C was significantly correlated with the infiltration levels of various immune cells. It was speculated that heteromeric LRRC8A/C channels mediate cGAMP transport and p53 signaling in T cells [36].

5. Conclusions and Prospects

In conclusion, LRRC8A can affect the prognosis of PAAD and was correlated to cell proliferation, cell migration, drug resistance, and immune infiltration. Specifically, functional analysis indicated that LRRC8A influences the progression and prognosis of PAAD via the regulation of CD8+ T cells immune infiltration. We hope that our results provide novel insights into assisting in the designing of new therapeutic drugs for PAAD.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14225526/s1, Supplementary Figure S1: bioinformatic analysis of LRRC8B–E in PAAD; Supplementary Figure S2: immune infiltration analysis of LRRC8B–8E in PAAD; File S1: the whole western blots.

Author Contributions

Conceptualization, R.X. and C.S.; methodology, Y.H.; validation, Q.X., H.Z. and H.S.; data curation, C.Z.; writing—original draft preparation, R.X.; writing—review and editing, X.W.; visualization, Y.Z.; funding acquisition, R.X. and C.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation Program of China (Grant No. 32070532) and Natural Science Foundation of Shaanxi Province (Grant No. 2022JQ-771; 2021PT-037).

Institutional Review Board Statement

Tumor specimens were obtained with the consent of the patients and all human specimen experiments were approved by the Medical Ethics Committee of Xijing Hospital (No. 2015432).

Informed Consent Statement

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

Data Availability Statement

Authors declared that all and the other data supporting the findings of this study are available within the paper. The raw data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.L.; Miller, K.D.; Fuchs, H.E.; Jemal, A. Cancer Statistics, 2021. CA Cancer J. Clin. 2021, 71, 7–33. [Google Scholar] [CrossRef] [PubMed]
  2. Bazeed, A.Y.; Day, C.M.; Garg, S. Pancreatic Cancer: Challenges and Opportunities in Locoregional Therapies. Cancers 2022, 14, 4257. [Google Scholar] [CrossRef] [PubMed]
  3. Vincent, A.; Herman, J.; Schulick, R.; Hruban, R.H.; Goggins, M. Pancreatic cancer. Lancet 2011, 378, 607–620. [Google Scholar] [CrossRef]
  4. Tian, Y.F.; Ning, C.F.; He, F.; Yin, B.C.; Ye, B.C. Highly sensitive detection of exosomes by SERS using gold nanostar@Raman reporter@nanoshell structures modified with a bivalent cholesterol-labeled DNA anchor. Analyst 2018, 143, 4915–4922. [Google Scholar] [CrossRef]
  5. Xie, Z.; Chen, X.; Li, J.; Guo, Y.; Li, H.; Pan, X.; Jiang, J.; Liu, H.; Wu, B. Salivary HOTAIR and PVT1 as novel biomarkers for early pancreatic cancer. Oncotarget 2016, 7, 25408–25419. [Google Scholar] [CrossRef] [Green Version]
  6. Sogawa, K.; Takano, S.; Iida, F.; Satoh, M.; Tsuchida, S.; Kawashima, Y.; Yoshitomi, H.; Sanda, A.; Kodera, Y.; Takizawa, H.; et al. Identification of a novel serum biomarker for pancreatic cancer, C4b-binding protein α-chain (C4BPA) by quantitative proteomic analysis using tandem mass tags. Br. J. Cancer 2016, 115, 949–956. [Google Scholar] [CrossRef]
  7. Melo, S.A.; Sugimoto, H.; O’Connell, J.T.; Kato, N.; Villanueva, A.; Vidal, A.; Qiu, L.; Vitkin, E.; Perelman, L.T.; Melo, C.A.; et al. Cancer exosomes perform cell-independent microRNA biogenesis and promote tumorigenesis. Cancer Cell 2014, 26, 707–721. [Google Scholar] [CrossRef] [Green Version]
  8. Tan, E.; El-Rayes, B. Pancreatic Cancer and Immunotherapy: Resistance Mechanisms and Proposed Solutions. J. Gastrointest. Cancer 2019, 50, 1–8. [Google Scholar] [CrossRef]
  9. Javle, M.M.; Shroff, R.T.; Xiong, H.; Varadhachary, G.A.; Fogelman, D.; Reddy, S.A.; Davis, D.; Zhang, Y.; Wolff, R.A.; Abbruzzese, J.L. Inhibition of the mammalian target of rapamycin (mTOR) in advanced pancreatic cancer: Results of two phase II studies. BMC Cancer 2010, 10, 368. [Google Scholar] [CrossRef] [Green Version]
  10. Ramanathan, R.K.; McDonough, S.L.; Philip, P.A.; Hingorani, S.R.; Lacy, J.; Kortmansky, J.S.; Thumar, J.; Chiorean, E.G.; Shields, A.F.; Behl, D.; et al. Phase IB/II Randomized Study of FOLFIRINOX Plus Pegylated Recombinant Human Hyaluronidase Versus FOLFIRINOX Alone in Patients With Metastatic Pancreatic Adenocarcinoma: SWOG S1313. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2019, 37, 1062–1069. [Google Scholar] [CrossRef]
  11. Anderson, K.J.; Cormier, R.T.; Scott, P.M. Role of ion channels in gastrointestinal cancer. World J. Gastroenterol. 2019, 25, 5732–5772. [Google Scholar] [CrossRef] [PubMed]
  12. Hanahan, D. Hallmarks of Cancer: New Dimensions. Cancer Discov. 2022, 12, 31–46. [Google Scholar] [CrossRef] [PubMed]
  13. Hanahan, D.; Weinberg, R.A. The hallmarks of cancer. Cell 2000, 100, 57–70. [Google Scholar] [CrossRef] [Green Version]
  14. Xu, R.; Wang, X.; Shi, C. Volume-regulated anion channel as a novel cancer therapeutic target. Int. J. Biol. Macromol. 2020, 159, 570–576. [Google Scholar] [CrossRef] [PubMed]
  15. Prevarskaya, N.; Skryma, R.; Shuba, Y. Ion Channels in Cancer: Are Cancer Hallmarks Oncochannelopathies? Physiol. Rev. 2018, 98, 559–621. [Google Scholar] [CrossRef] [Green Version]
  16. Rubino, S.; Bach, M.D.; Schober, A.L.; Lambert, I.H.; Mongin, A.A. Downregulation of Leucine-Rich Repeat-Containing 8A Limits Proliferation and Increases Sensitivity of Glioblastoma to Temozolomide and Carmustine. Front. Oncol. 2018, 8, 142. [Google Scholar] [CrossRef] [Green Version]
  17. Konishi, T.; Shiozaki, A.; Kosuga, T.; Kudou, M.; Shoda, K.; Arita, T.; Konishi, H.; Komatsu, S.; Kubota, T.; Fujiwara, H.; et al. LRRC8A Expression Influences Growth of Esophageal Squamous Cell Carcinoma. Am. J. Pathol. 2019, 189, 1973–1985. [Google Scholar] [CrossRef]
  18. Lu, P.; Ding, Q.; Li, X.; Ji, X.; Li, L.; Fan, Y.; Xia, Y.; Tian, D.; Liu, M. SWELL1 promotes cell growth and metastasis of hepatocellular carcinoma in vitro and in vivo. EBioMedicine 2019, 48, 100–116. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, H.; Deng, Z.; Zhang, D.; Li, H.; Zhang, L.; Niu, J.; Zuo, W.; Fu, R.; Fan, L.; Ye, J.H.; et al. High expression of leucine-rich repeat-containing 8A is indicative of a worse outcome of colon cancer patients by enhancing cancer cell growth and metastasis. Oncol. Rep. 2018, 40, 1275–1286. [Google Scholar] [CrossRef]
  20. Zhou, C.; Chen, X.; Planells-Cases, R.; Chu, J.; Wang, L.; Cao, L.; Li, Z.; López-Cayuqueo, K.I.; Xie, Y.; Ye, S.; et al. Transfer of cGAMP into Bystander Cells via LRRC8 Volume-Regulated Anion Channels Augments STING-Mediated Interferon Responses and Anti-viral Immunity. Immunity 2020, 52, 767–781.e766. [Google Scholar] [CrossRef]
  21. Hu, D.; Ansari, D.; Zhou, Q.; Sasor, A.; Hilmersson, K.S.; Bauden, M.; Jiang, Y.; Andersson, R. Calcium-activated chloride channel regulator 1 as a prognostic biomarker in pancreatic ductal adenocarcinoma. BMC Cancer 2018, 18, 1096. [Google Scholar] [CrossRef] [PubMed]
  22. Wu, P.; Xu, R.; Chen, X.; Zhao, Y.; Tan, D.; Zhao, Y.; Qin, W.; Zhang, C.; Ge, X.; Shi, C. Establishment and characterization of patient-derived xenografts for hormone-naïve and castrate-resistant prostate cancers to improve treatment modality evaluation. Aging 2020, 12, 3848–3861. [Google Scholar] [CrossRef] [PubMed]
  23. Li, X.; Huo, C.; Xiao, Y.; Xu, R.; Liu, Y.; Jia, X.; Wang, X. Bisdemethoxycurcumin Protection of Cardiomyocyte Mainly Depends on Nrf2/HO-1 Activation Mediated by the PI3K/AKT Pathway. Chem. Res. Toxicol. 2019, 32, 1871–1879. [Google Scholar] [CrossRef] [PubMed]
  24. Huo, C.; Liu, Y.; Li, X.; Xu, R.; Jia, X.; Hou, L.; Wang, X. LRRC8A contributes to angiotensin II-induced cardiac hypertrophy by interacting with NADPH oxidases via the C-terminal leucine-rich repeat domain. Free Radic. Biol. Med. 2021, 165, 191–202. [Google Scholar] [CrossRef]
  25. Deer, E.L.; González-Hernández, J.; Coursen, J.D.; Shea, J.E.; Ngatia, J.; Scaife, C.L.; Firpo, M.A.; Mulvihill, S.J. Phenotype and genotype of pancreatic cancer cell lines. Pancreas 2010, 39, 425–435. [Google Scholar] [CrossRef] [Green Version]
  26. Bailey, P.; Chang, D.K.; Nones, K.; Johns, A.L.; Patch, A.M.; Gingras, M.C.; Miller, D.K.; Christ, A.N.; Bruxner, T.J.; Quinn, M.C.; et al. Genomic analyses identify molecular subtypes of pancreatic cancer. Nature 2016, 531, 47–52. [Google Scholar] [CrossRef]
  27. Hu, L.P.; Zhang, X.X.; Jiang, S.H.; Tao, L.Y.; Li, Q.; Zhu, L.L.; Yang, M.W.; Huo, Y.M.; Jiang, Y.S.; Tian, G.A.; et al. Targeting Purinergic Receptor P2Y2 Prevents the Growth of Pancreatic Ductal Adenocarcinoma by Inhibiting Cancer Cell Glycolysis. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2019, 25, 1318–1330. [Google Scholar] [CrossRef] [Green Version]
  28. Conejo-Garcia, J.R. Breaking barriers for T cells by targeting the EPHA2/TGF-β/COX-2 axis in pancreatic cancer. J. Clin. Investig. 2019, 129, 3521–3523. [Google Scholar] [CrossRef] [Green Version]
  29. Markosyan, N.; Li, J.; Sun, Y.H.; Richman, L.P.; Lin, J.H.; Yan, F.; Quinones, L.; Sela, Y.; Yamazoe, T.; Gordon, N.; et al. Tumor cell-intrinsic EPHA2 suppresses anti-tumor immunity by regulating PTGS2 (COX-2). J. Clin. Investig. 2019, 129, 3594–3609. [Google Scholar] [CrossRef] [Green Version]
  30. Lahey, L.J.; Mardjuki, R.E.; Wen, X.; Hess, G.T.; Ritchie, C.; Carozza, J.A.; Böhnert, V.; Maduke, M.; Bassik, M.C.; Li, L. LRRC8A:C/E Heteromeric Channels Are Ubiquitous Transporters of cGAMP. Mol. Cell 2020, 80, 578–591.e575. [Google Scholar] [CrossRef]
  31. Ramírez, A.; García-Quiroz, J.; Aguilar-Eslava, L.; Sánchez-Pérez, Y.; Camacho, J. Novel Therapeutic Approaches of Ion Channels and Transporters in Cancer. In Targets of Cancer Diagnosis and Treatment; Reviews of Physiology, Biochemistry and Pharmacology Series; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar] [CrossRef]
  32. Zuccolini, P.; Ferrera, L.; Remigante, A.; Picco, C.; Barbieri, R.; Bertelli, S.; Moran, O.; Gavazzo, P.; Pusch, M. The VRAC blocker DCPIB directly gates the BK channels and increases intracellular Ca(2+) in melanoma and pancreatic duct adenocarcinoma cell lines. Br. J. Pharmacol. 2022, 179, 3452–3469. [Google Scholar] [CrossRef] [PubMed]
  33. Qiu, Z.; Dubin, A.E.; Mathur, J.; Tu, B.; Reddy, K.; Miraglia, L.J.; Reinhardt, J.; Orth, A.P.; Patapoutian, A. SWELL1, a plasma membrane protein, is an essential component of volume-regulated anion channel. Cell 2014, 157, 447–458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Voss, F.K.; Ullrich, F.; Münch, J.; Lazarow, K.; Lutter, D.; Mah, N.; Andrade-Navarro, M.A.; von Kries, J.P.; Stauber, T.; Jentsch, T.J. Identification of LRRC8 heteromers as an essential component of the volume-regulated anion channel VRAC. Science 2014, 344, 634–638. [Google Scholar] [CrossRef] [Green Version]
  35. Kumar, L.; Chou, J.; Yee, C.S.; Borzutzky, A.; Vollmann, E.H.; von Andrian, U.H.; Park, S.Y.; Hollander, G.; Manis, J.P.; Poliani, P.L.; et al. Leucine-rich repeat containing 8A (LRRC8A) is essential for T lymphocyte development and function. J. Exp. Med. 2014, 211, 929–942. [Google Scholar] [CrossRef] [PubMed]
  36. Concepcion, A.R.; Wagner, L.E., 2nd; Zhu, J.; Tao, A.Y.; Yang, J.; Khodadadi-Jamayran, A.; Wang, Y.H.; Liu, M.; Rose, R.E.; Jones, D.R.; et al. The volume-regulated anion channel LRRC8C suppresses T cell function by regulating cyclic dinucleotide transport and STING-p53 signaling. Nat. Immunol. 2022, 23, 287–302. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Bioinformatics analyses of LRRC8A in PAAD. (a) Pan-cancer analysis of LRRC8A–E expression by ONCOMINE; (b) the expression of LRRC8A in PAAD was seven-fold higher than that in normal by GEPIA, Tumor (n = 179), Normal (n = 171); (c) LRRC8A expression in PAAD based on individual cancer stages; (d,e) the DFS and OS of patients with PAAD with a high LRRC8A expression were significantly worse than those with a low LRRC8A expression; (f) the RFS of patients with PAAD with a high LRRC8A expression were significantly worse than those with a low LRRC8A expression.
Figure 1. Bioinformatics analyses of LRRC8A in PAAD. (a) Pan-cancer analysis of LRRC8A–E expression by ONCOMINE; (b) the expression of LRRC8A in PAAD was seven-fold higher than that in normal by GEPIA, Tumor (n = 179), Normal (n = 171); (c) LRRC8A expression in PAAD based on individual cancer stages; (d,e) the DFS and OS of patients with PAAD with a high LRRC8A expression were significantly worse than those with a low LRRC8A expression; (f) the RFS of patients with PAAD with a high LRRC8A expression were significantly worse than those with a low LRRC8A expression.
Cancers 14 05526 g001
Figure 2. Roles of LRRC8A in PAAD in vitro. (a) Cell viability of AsPC-1 and AsPC-1GEMR cells on GEM; (b) cell viability of AsPC-1GEMR cells on GEM after pretreated with shLRRC8A or scramble lentivirus (NC); (c) mRNA expression of LRRC8A in PAAD cell lines; (d) LRRC8A protein expression in PAAD cell lines; (e) representative immunohistochemical results of LRRC8A in patients with PAAD and quantitative analysis results; (f) wound healing assay using BxPC-3 cells and AsPC-1GEM-R cells with LRRC8A knockdown and quantitative analysis results; (g) transwell experiments using AsPC-1 and BxPC-3 cells, and quantitative analysis results; (h) expression of LRRC8A and MMP2 after AsPC-1, AsPC-1GEM-R and BxPC-3 cells were treated by shLRRC8A. Error bars indicate SE. * p < 0.05, ** p < 0.01 and *** p < 0.001 by Student’s t-tests. The whole western blots were shown in File S1. ns: not significant.
Figure 2. Roles of LRRC8A in PAAD in vitro. (a) Cell viability of AsPC-1 and AsPC-1GEMR cells on GEM; (b) cell viability of AsPC-1GEMR cells on GEM after pretreated with shLRRC8A or scramble lentivirus (NC); (c) mRNA expression of LRRC8A in PAAD cell lines; (d) LRRC8A protein expression in PAAD cell lines; (e) representative immunohistochemical results of LRRC8A in patients with PAAD and quantitative analysis results; (f) wound healing assay using BxPC-3 cells and AsPC-1GEM-R cells with LRRC8A knockdown and quantitative analysis results; (g) transwell experiments using AsPC-1 and BxPC-3 cells, and quantitative analysis results; (h) expression of LRRC8A and MMP2 after AsPC-1, AsPC-1GEM-R and BxPC-3 cells were treated by shLRRC8A. Error bars indicate SE. * p < 0.05, ** p < 0.01 and *** p < 0.001 by Student’s t-tests. The whole western blots were shown in File S1. ns: not significant.
Cancers 14 05526 g002
Figure 3. Function analysis of LRRC8A in PAAD. (a) GO enrichment analysis (biological process, molecular function, and cellular component) of the top 50 positively correlated genes of LRRC8A; (b) KEGG enrichment analysis of the top 50 positively correlated genes of LRRC8A; (c) protein interaction network of LRRC8A–E; (d) gene–gene interaction network of LRRC8A–E and the top 50 positively correlated genes of LRRC8A.
Figure 3. Function analysis of LRRC8A in PAAD. (a) GO enrichment analysis (biological process, molecular function, and cellular component) of the top 50 positively correlated genes of LRRC8A; (b) KEGG enrichment analysis of the top 50 positively correlated genes of LRRC8A; (c) protein interaction network of LRRC8A–E; (d) gene–gene interaction network of LRRC8A–E and the top 50 positively correlated genes of LRRC8A.
Cancers 14 05526 g003
Figure 4. Correlation analysis of LRRC8A with P2RY2 and EPHA2 in PAAD. (a) LRRC8A showed a positive correlation with P2RY2 and EPHA2 in cBioPortal; (b) LRRC8A showed a positive correlation with P2RY2 and EPHA2 in TIMER.
Figure 4. Correlation analysis of LRRC8A with P2RY2 and EPHA2 in PAAD. (a) LRRC8A showed a positive correlation with P2RY2 and EPHA2 in cBioPortal; (b) LRRC8A showed a positive correlation with P2RY2 and EPHA2 in TIMER.
Cancers 14 05526 g004
Figure 5. Immune infiltration analysis of LRRC8A in PAAD. (a) EPHA2 expression decreased in BxPC3 (p * < 0.05) and AsPC1GEM-R cells (p * < 0.05) after LRRC8A knockdown; (b) immune infiltration of LRRC8A in PAAD; (c) correlation analysis of LRRC8A and CD8+ T cells; (d) representative immunohistochemical results of CD8 in patients with PAAD; (e) correlation analysis of LRRC8A and EPHA2 in patients with PAAD; (f) quantitative analysis results of Figure 5d. ns: not significant.
Figure 5. Immune infiltration analysis of LRRC8A in PAAD. (a) EPHA2 expression decreased in BxPC3 (p * < 0.05) and AsPC1GEM-R cells (p * < 0.05) after LRRC8A knockdown; (b) immune infiltration of LRRC8A in PAAD; (c) correlation analysis of LRRC8A and CD8+ T cells; (d) representative immunohistochemical results of CD8 in patients with PAAD; (e) correlation analysis of LRRC8A and EPHA2 in patients with PAAD; (f) quantitative analysis results of Figure 5d. ns: not significant.
Cancers 14 05526 g005
Table 1. Primer sequences for Real-time PCR.
Table 1. Primer sequences for Real-time PCR.
Gene NameSequences (5′ to 3′)
LRRC8A FORWARDGGAGAGCAGCTACAGCGACA
LRRC8A REVERSETCACTCACCTCCGACAGGAA
EPHA2 FORWARDAGAAGCGCCTGTTCACCAAG
EPHA2 REVERSEGCTCCTCCACGTTCAGCTTC
ACTIN FORWARDCATGTACGTTGCTATCCAGGC
ACTIN REVERSECTCCTTAATGTCACGCACGAT
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, R.; Hu, Y.; Xie, Q.; Zhang, C.; Zhao, Y.; Zhang, H.; Shi, H.; Wang, X.; Shi, C. LRRC8A Is a Promising Prognostic Biomarker and Therapeutic Target for Pancreatic Adenocarcinoma. Cancers 2022, 14, 5526. https://doi.org/10.3390/cancers14225526

AMA Style

Xu R, Hu Y, Xie Q, Zhang C, Zhao Y, Zhang H, Shi H, Wang X, Shi C. LRRC8A Is a Promising Prognostic Biomarker and Therapeutic Target for Pancreatic Adenocarcinoma. Cancers. 2022; 14(22):5526. https://doi.org/10.3390/cancers14225526

Chicago/Turabian Style

Xu, Rong, Yaohua Hu, Qinghua Xie, Caiqin Zhang, Yong Zhao, He Zhang, Hailong Shi, Xiaoming Wang, and Changhong Shi. 2022. "LRRC8A Is a Promising Prognostic Biomarker and Therapeutic Target for Pancreatic Adenocarcinoma" Cancers 14, no. 22: 5526. https://doi.org/10.3390/cancers14225526

APA Style

Xu, R., Hu, Y., Xie, Q., Zhang, C., Zhao, Y., Zhang, H., Shi, H., Wang, X., & Shi, C. (2022). LRRC8A Is a Promising Prognostic Biomarker and Therapeutic Target for Pancreatic Adenocarcinoma. Cancers, 14(22), 5526. https://doi.org/10.3390/cancers14225526

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