KCa1.1, a calcium-activated potassium channel subunit alpha 1, is targeted by miR-17-5p and modulates cell migration in malignant pleural mesothelioma
© Lin et al. 2016
Received: 1 October 2015
Accepted: 20 May 2016
Published: 1 June 2016
Malignant pleural mesothelioma (MPM) is an aggressive, locally invasive, cancer elicited by asbestos exposure and almost invariably a fatal diagnosis. To date, we are one of the leading laboratory that compared microRNA expression profiles in MPM and normal mesothelium samples in order to identify dysregulated microRNAs with functional roles in mesothelioma. We interrogated a significant collection of MPM tumors and normal pleural samples in our biobank in search for novel therapeutic targets.
Utilizing mRNA-microRNA correlations based on differential gene expression using Gene Set Enrichment Analysis (GSEA), we systematically combined publicly available gene expression datasets with our own MPM data in order to identify candidate targets for MPM therapy.
We identified enrichment of target binding sites for the miR-17 and miR-30 families in both MPM tumors and cell lines. RT-qPCR revealed that members of both families were significantly downregulated in MPM tumors and cell lines. Interestingly, lower expression of miR-17-5p (P = 0.022) and miR-20a-5p (P = 0.026) was clearly associated with epithelioid histology. We interrogated the predicted targets of these differentially expressed microRNA families in MPM cell lines, and identified KCa1.1, a calcium-activated potassium channel subunit alpha 1 encoded by the KCNMA1 gene, as a target of miR-17-5p. KCa1.1 was overexpressed in MPM cells compared to the (normal) mesothelial line MeT-5A, and was also upregulated in patient tumor samples compared to normal mesothelium. Transfection of MPM cells with a miR-17-5p mimic or KCNMA1-specific siRNAs reduced mRNA expression of KCa1.1 and inhibited MPM cell migration. Similarly, treatment with paxilline, a small molecule inhibitor of KCa1.1, resulted in suppression of MPM cell migration.
These functional data implicating KCa1.1 in MPM cell migration support our integrative approach using MPM gene expression datasets to identify novel and potentially druggable targets.
KeywordsKCNMA1 miR-17-5p Mesothelioma Therapeutic targets KCa1.1 microRNA Integrative analysis
Malignant pleural mesothelioma (MPM) is an aggressive tumor occurring in the lining of the lungs, induced by exposure to asbestos. MPM has poor prognosis, and palliative chemotherapy is often the only treatment modality that can be offered [1, 2]. Since the adoption of cisplatin and pemetrexed as the standard of care  a decade ago, there has been little therapeutic progress and the identification of new therapeutic targets for MPM is an urgent unmet need.
Genome-wide gene expression profiling studies using microarray and next generation sequencing (NGS) have facilitated identification of disease-specific expression profiles and many are publicly available. While most studies have focused on the identification of a single therapeutic candidate, the greatest challenge remains the interpretation of data within the context of cancer cell biology. Bioinformatic tools such as Gene Set Enrichment Analysis (GSEA)  have improved interpretation of microarray and NGS data for downstream functional validation. In MPM, gene expression profiling studies have identified several novel targets including MMP14 , ALCAM , NME2, CRI1, PDGFC and GSN . However, there was little commonality between these studies, and to date no pharmaceutical approach to targeting these candidates has been developed. At a systems level, pathway analysis has revealed enrichment of genes in MPM belonging to cellular processes such as cellular metabolism, cytoskeletal re-organization, apoptosis, spindle checkpoint and cell cycle progression and regulation [5, 8, 9]. Many of these pathways, however, have not been explored in detail.
Since MPM is characterized by alterations in multiple genes, we hypothesized that a strategy to inhibit and/or restore a single target gene is unlikely to be effective. In comparison, new insights into the involvement of microRNAs in the regulation of MPM growth  have provided an alternative way to inhibit MPM growth with the potential to be successfully translated into a new therapeutic approach for MPM . MicroRNAs are small non-coding RNAs involved in post-transcriptional control of gene expression . They form a complex network where each microRNA regulates multiple mRNAs and each mRNA is regulated by multiple microRNAs. Changes in microRNA expression are associated with proliferation and drug resistance of cancer cells, and microRNAs can act as oncogenes or tumor suppressors [13–16]. Making use of data from our previous studies [17–19], we present here an integrative approach by comparing microRNA and mRNA gene expression datasets to identify enriched biological themes that can be translated into potential druggable targets for MPM, as well as functional data revealing that KCa1.1 is a potential therapeutic target in MPM.
Results and discussion
Identification of target binding site of differentially expressed genes in MPM cell lines and tumors (enriched microRNA binding sites)
Binding sites for miR-17 and miR-30 microRNA families are enriched in all gene expression datasets
The miR-17 family is consistently downregulated in MPM and its predicted targets are associated with multiple cancer-related pathways
We validated microRNA expression using RT-qPCR and found consistent downregulation of multiple miR-17 family members and miR-30e in MPM in a set of 23 normal pleura and 60 formalin-fixed, paraffin-embedded (FFPE) tumor samples from our biobank, as well as a panel of 7 MPM and mesothelial cell lines (Fig. 2b and c). Lower expression of miR-17-5p (P = 0.022) and miR-20a-5p (P = 0.026) was significantly associated with epithelioid histology, as was the co-expressed miR-19b-3p (P = 0.016). Since miR-30e-5p expression did not show any association with clinical characteristics in MPM, we focused further attention on the miR-17 family.
Previous studies addressing the role of miR-17 in cancer biology have shown this to be complex. While miR-17 expression is frequently reported to be upregulated as part of the miR-17 ~ 92 cluster , miR-17-5p downregulation and/or loss of heterozygosity/gene deletion at 13q31 has been reported in various tumors . Thus the role of these microRNAs is likely to be context specific. Pathway enrichment analysis revealed that the predicted targets of miR-17 family are enriched in several key cancer-related signaling pathways previously implicated in MPM biology and treatment, including MAPK signaling , ErbB signaling , Focal adhesion , TNF signaling  and TGF-beta signaling pathways  (Additional file 1: Table S2). Furthermore, the top four enriched pathways share genes linked to MPM cell migration, in particular MAPK1, P = 0.03, −2.2 fold in MPM (MAPK, Axon Guidance, Pancreatic Cancer pathways) and TGFBR2, P = 0.0002, > 5 fold in MPM (Pancreatic Cancer, MAPK, Endocytosis pathways) (Fig. 2d). Specifically, key target genes of the miR-17 family are coordinately regulated in signaling pathways contributing to MPM cell migration (Fig. 2e), where downregulation of miR-17 is associated with increased DUSP2 expression, an inhibitor of MAPK1. Interestingly, DUSP2 has been shown to be downregulated by the miR-17 family member miR-20a . Decreased miR-17 also leads to increased KCNMA1 expression (which in turn interacts with the MAPK signaling pathway , and increased TGFBR2  and altered TGFβ2 signaling. These gene expression changes indicate that all these pathways together contribute to an increase in MPM cell migration (Fig. 2d and e), and collectively these results implicate a broad impact of the miR-17 family in MPM biology.
Reduced miR-17-5p expression is correlated with increased expression of KCNMA1
We investigated the correlation between miR-17-5p expression and mRNA expression of its predicted targets from the same MPM cell lines  and identified that 20 out of 40 target genes demonstrated the expected inverse relationship between miRNA and mRNA levels, P <0.05 (Fig. 2f). Prioritizing these gene targets based firstly on those that showed significant differential gene expression (>5-fold change) and secondly on the availability of small molecule inhibitors/drugs, led to identification of KCNMA1. This gene encodes the calcium-activated potassium channel subunit alpha 1, KCa1.1 previously reported to be upregulated in prostate , breast  and other cancers [39, 40]. RT-qPCR confirmed upregulation of KCNMA1 mRNA in MPM cell lines compared to the normal immortalized mesothelial line MeT-5A (Fig. 2g), which was associated with decreased expression of miR-17-5p (Fig. 2c). We further analyzed KCNMA1 gene expression in fresh-frozen samples of normal pleura and MPM from our biobank. In common with other tumor types, we observed increased expression of KCNMA1 (Fig. 2h) in MPM tumor samples as well as downregulation of miR-17-5p (Fig. 2i). Immunofluorescence studies further indicated that MPM tumor samples with high tumor cell content expressed high levels of KCal.1 (Fig. 2j right); in contrast, no expression was detected in areas without tumor cells (Fig. 2j left).
Potassium ion channels such as KCa1.1 are the subject of increasing interest in cancer research due to their observed effects on cell processes including cell proliferation, cell adhesion, angiogenesis, cell migration and metastasis (reviewed in Pardo ). For example, expression of Kv10.1 channel (EAG1, product of the KCNH1 gene) is normally limited to the brain, but it is overexpressed in many tumor types and knockdown in overexpressing cells was found to reduce viability . Furthermore, there is evidence that altered microRNA expression is involved in the regulation of potassium channels in cancer. In glioblastoma  for example, increased EAG1 protein expression was associated with decreased miR-296-3p levels and in osteosarcoma , increased EAG1 is correlated with reduced miR-34a. Kv11.1 (HERG) is upregulated in pancreatic cancer as a consequence of miR-96 downregulation  and miR-211 expression is inversely correlated with KCa1.1 in melanoma . In each case, ectopic expression of these microRNAs reduced channel expression and inhibited proliferation of the cancer cells.
KCNMA1 gene is a direct target of miR-17-5p and genetic and pharmacological inhibition of its protein, KCa1.1 modulates migration in MPM cells
Pathway analysis shows that KCNMA1 is also involved in a range of cellular processes, including Axon Guidance, Focal Adhesion and Wnt Signaling, so we next investigated the effects of inhibiting KCa1.1 in MPM cells. Silencing KCNMA1 expression with siRNAs or miR-17-5p mimics had limited effects on the proliferation of a panel of MPM cell lines (Additional file 1: Figure S2), and did not lead to changes in the cell cycle (Additional file 1: Figure S3) or induction of apoptosis (Additional file 1: Figure S4) in transfected cells. This is in contrast to the inhibition of breast cancer  and melanoma  cell lines following KCNMA1 knockdown. Both miR-17-5p and KCNMA1-specific siRNAs, however reduced migration in H28 (Fig. 3d) and other MPM cell lines (Additional file 1 Figure S5), and the invasive capacity of cells in a modified agarose spot assay (Additional file 1: Figure S6). Similarly, modulation of KCa1.1 activity with the inhibitor paxilline suppressed migration of MPM cell lines. Paxilline inhibited cell proliferation at concentrations greater than 25 μM (Fig. 3e) and treating MPM cells at a sub-lethal dose (12.5 μM) resulted in significant inhibition of migration (Fig. 3d) and cell colony formation (Fig. 3f). Interestingly, MPM cell lines were less sensitive than breast or prostate cancer lines, in which growth was inhibited by 15 μM paxilline . These data are consistent with previous results implicating a role in cell migration for calcium-activated potassium channels in general  and that KCa1.1 in particular, has been shown to be important in the migration of glioma cells , invasiveness of melanoma cells , and transendothelial migration in metastatic breast cancer [48, 49].
Cell migration and invasion is dependent on calcium flux, which contributes to changes in cellular volume [50, 51]. The blocker-sensitive outward potassium current generated by large conductance potassium channels such as KCa1.1 can functionally couple with other ion channels to control calcium flux, thus facilitating migration and invasion by regulating cell volume . In MPM cells, KCNMA1 knockdown led to an increase in the basal intracellular calcium concentration (Fig. 3g). As KCNMA1 is the major contributor to the resting membrane potential, this is likely due to a partial membrane depolarization and higher activity of Ca2+ channels [53, 54]. Furthermore, KCNMA1 silencing resulted in a more prolonged elevation in submembrane Ca2+ levels in response to high K+ containing buffer (Fig. 3h). Together, these experiments demonstrate a major impact of KCNMA1 silencing on Ca2+ handling by MPM cells. In light of these observations it is interesting that changes in basal Ca2+ levels following KCNMA1 knockdown had negligible effects on cell growth, despite altered Ca2+ handling being linked to apoptosis resistance in MPM cells . Further research is required to understand exactly how the loss of KCNMA1 activity and resultant changes in Ca2+ flux reduces migration and invasion without impacting apoptosis, and is beyond the scope of the current study.
Potassium channels have been implicated in drug resistance of cancer cells. One example is provided by the miR-296-3p-mediated downregulation of EAG1 leading to reversal of anticancer drug resistance in glioma cells . There is conflicting evidence, however, for the role of potassium channels in the toxicity of cisplatin, a platinum-based drug frequently used in chemotherapy for MPM. Increased expression of the potassium channel Kv10.1 (EAG, KCNH1) was linked to the drug resistance of glioblastoma cells , whereas inhibition of the elevated levels of Kv11.1 (HERG, KCNH2) induced by cisplatin selection had no effect on the toxicity of cisplatin in epidermal or liver cancer cells. In contrast, it was found to be essential for cisplatin-induced apoptosis in gastric cancer . We did not, however, observe sensitization of MPM to either cisplatin or gemcitabine in the presence of paxilline (Additional file 1: Figure S7), suggesting that KCa1.1 is not involved in the cytotoxic activity of these drugs in MPM, and is not solely responsible for the potassium flux modulated by amphotericin B and bumetanide previously observed in MPM .
Therapeutic targeting of potassium channels in cancer
Growing evidence suggests that the potassium channels Kv10.1, Kv11.1, KCa1.1 and KCa3.1 have important roles in cancer cell invasion and metastasis [38, 41], and as MPM is a locally invasive cancer with high migratory capacity, that targeting these channels is a potential therapeutic option. Kv11.1 (hERG, encoded by the KCNH2 gene) is overexpressed in many cancers  and our MPM cell lines (1.78-fold upregulation, P = 0.03). First linked to cancer through the finding that resting membrane potential in neuroblastoma cells was linked to proliferation, Kv11.1 has since been found to be expressed at high levels in multiple tumor types . While this potentially makes Kv11.1 a good target, it is also expressed in the heart where it is related to long QT syndrome, making it difficult to envisage a cancer-specific therapeutic window for Kv11.1 inhibition. The Kv10.1 (EAG1, product of the KCNH1 gene) channel is normally only expressed in the brain, but is upregulated in multiple cancer types , and its inhibition or knockdown reduces cell migration and viability . The important role of this channel in cancer was demonstrated by its ability to confer invasive growth characteristics to otherwise non-tumorigenic cell lines in xenograft models . However, its close evolutionary relationship to Kv11.1 reduces the possibility that drugs targeting Kv10.1 will be sufficiently specific.
Other potassium channels may represent preferable targets to control cancer cell migration. For example, inhibition of KCa2.3 reduced bone metastases in a model of breast cancer , and blocking KCa1.1, which was induced by ionizing radiation, inhibited the increase in glioma cell migration . Targeting KCa1.1 in MPM could be achieved using a channel blocker such as paxilline. While an inhibitory dose may be difficult to reach via systemic inhibition, the typical restriction of MPM tumor growth to the pleural cavity provides an alternative loco-regional route for delivery of paxilline or other potassium channel blockers. Alternatively, the post-transcriptional regulation of KCNMA1 by miR-17-5p suggests that KCa1.1 could be an important target of a miR-17-5p mimic-based therapy.
We have shown here that an integrated approach, combining publicly available gene expression datasets, is an effective and practical strategy to identify therapeutic targets in MPM. These targets can be both proteins that are inhibited via traditional pharmacological interventions, as well as microRNAs themselves, as illustrated by a near-complete response in a MPM patient failing on standard chemotherapy, who received a miR-16 mimic packaged in (targeted) nanocells , revealing that microRNA-based treatment concepts are valid . We have demonstrated that KCa1.1 is a direct target of miR-17-5p and inhibition of this combination modulates cell migration in MPM.
Human mesothelioma cell lines H28, H226, H2052, H2452, MSTO and the immortalized normal mesothelial line MeT-5A were obtained from the American Type Cell Culture repository (ATCC, Rockville, USA) and MM05  and REN  were described previously and kindly provided by collaborators. All cells were grown in RPMI with 10 % fetal bovine serum (FBS) at 37 °C with 5 % CO2. REN cells were grown in Ham’s F12 medium. All medium and FBS were obtained from Life Technologies or Sigma.
Formalin-fixed paraffin embedded (FFPE) tumor tissues used in this study were part of a reported series of extrapleural pneumonectomy (EPP) patients collected between 1994 and 2009 from the Royal Prince Alfred Hospital (RPAH) or Strathfield Private Hospital, Sydney . Waiver of consent for the use of these archival samples was granted by the Human Research Ethics Committee (HREC) at Concord Repatriation General Hospital (CRGH), Sydney, Australia (CH62/6/2009/078). Demographics for this population have been described elsewhere . For the fresh-frozen tumor samples and normal pleural tissue samples, all patients gave informed written consent and the project was approved by the HREC at CRGH (HREC/11/CRGH/75) and RPAH (HREC/10/RPAH/599).
Total RNA was extracted from cell lines and fresh-frozen tissue samples using the TRIzol reagent (Life Technologies, Carlsbad, CA) and from FFPE tumors and normal pleura using the RNeasy FFPE kit (Qiagen, Hilden, Germany). RNA was quantified using a Nanophotometer (Implen, Munich, Germany) and quality was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). Samples with RNA integrity numbers (RINs) >8.0 were used for microarray analysis.
Microarray and TaqMan Low Density Array (TLDA) data acquisition
Microarray profiling experiments were performed using NCode Human Non-coding RNA microarrays (Life Technologies) according to MIAME guidelines, as previously described . All expression data is available at the National Centre for Biotechnology Information Gene Expression Omnibus [GSE48174]. TaqMan Low Density Array (TLDA) profiling was performed as previously described , following the profiling protocol without pre-amplification. TLDAs were run on a Viia7 Real-Time machine (Life Technologies). Data was analyzed using the 2-∆∆Cq method , with normalization to the average RNU6B Cq, and calculation of expression values was made relative to MeT-5A. MicroRNAs with Cq values >35 were excluded from analysis.
Quantitative RT-PCR in real time (RT-qPCR)
MicroRNA and mRNA expression levels from microarray analysis were validated using RT-qPCR in the five cell lines assayed using NCode and TLDA arrays. MicroRNA expression levels were also validated in fresh-frozen and FFPE tumor specimens from MPM patients and samples of normal mesothelium using TaqMan microRNA assays (see Additional file 1: Table S3 for assay IDs). Primers were designed for mRNA targets using the Universal Probe Library (UPL) algorithm provided by Roche (Additional file 1: Table S3). To quantify mRNA expression, total RNA (250 ng for cell lines) was reverse transcribed to cDNA using the Superscript III cDNA synthesis kit (Life Technologies) using oligo-dT primers (100 ng/μL). After reverse transcription, cDNA was diluted 1:5 with 2 μL of this product used as template in real-time qPCR. All reactions were run in triplicate on a Viia7 real-time machine (Life Technologies), using KAPA SYBR Fast qPCR master mix (Kapa Biosystems). All reactions had an initial enzyme inactivation step at 95 °C for 10 min followed by 40 cycles of 95 °C for 15 s and 55 °C for 30 s. 18S ribosomal RNA was used as the reference for qPCR data normalization and no template and no-RT samples included as negative controls. For microRNA expression, 100 ng of total RNA was reverse transcribed using pooled microRNA specific primers and the MicroRNA Reverse Transcription Kit (Life Technologies) as previously described . Specific TaqMan assays were used to amplify 10 ng cDNA using the KAPA Probe FAST qPCR master mix. Relative expression levels were calculated using the 2-ΔΔCq method as described  with MeT-5A designated a value of 1 (all fold-change values were calculated relative to this cell line). Genes were deemed technically replicated if the direction of expression was consistent with microarray data, and the magnitude of change was greater than 2-fold.
RNA interference (RNAi), paxilline treatment and growth assays
Knockdown of KCNMA1 was performed using two previously published siRNA sequences  and reintroduction of target microRNAs was performed using microRNA mimics. All siRNAs and microRNA mimics were obtained from Shanghai GenePharma. MPM cells were reverse transfected as described previously . Briefly, 2,500 cells were reverse transfected with varying concentrations of microRNA mimic or siRNA using Lipofectamine RNAiMAX (Life Technologies) at 0.1 μL per well. RNA was extracted 48 h after transfection to confirm knockdown efficiency using RT-qPCR. Cell growth assays were performed over a 5 d period to monitor changes in cell proliferation as described previously . Briefly, at the indicated time points, medium was aspirated from replicate plates, which were then frozen at −80 °C. At the conclusion of the experiments, relative cell numbers were determined by staining with 200 μL/well of SYBR Green I (Life Technologies) 1:4000 in a hypotonic lysis buffer (10 mM Tris HCl (pH 8), 5 mM EDTA, 0.1 % Triton X-100) overnight in the dark at 4 °C and then quantified by fluorimetry, measured using a FLUOstar Optima (BMG LabTech, Ortenberg, Germany) set to 485 nm excitation and 535 nm emission. Fluorescence intensity in siRNA or mimic-transfected cells was normalized to control-transfected cells. Each experiment was performed in triplicate. Growth inhibition using the small molecule inhibitor paxilline (Sigma Aldrich), in the presence or absence of cisplatin or gemcitabine was carried out using the same assay.
AGO2 Immunoprecipitation (AGO2-IP)
Cells were transfected with miR-17-5p or control mimic and AGO2 protein was immunoprecipitated as previously described [46, 47]. Total RNA isolation was carried out using TRIzol reagent, followed by RT-PCR with primers specific for KCNMA1 mRNA. PCR products were run on 2 % agar gel and visualized following ethidium bromide staining. Imaging was carried out using a Gel Logic 2200 Imaging System (Kodak) under non-saturating conditions. Densitometry of band intensity was carried out using the same software.
Immunofluorescence staining of KCa1.1 protein expression
Cells transfected with KCNMA1-specific siRNA, miR-17-5p mimic or control were fixed in paraformaldehyde solution (4 % v/v in PBS, Sigma, St. Louis, MO, USA) for 15 min, washed three times with PBS and permeabilized with 0.2 % Triton X-100 in PBS for 5 min. Fixed cells were blocked with PBS containing 0.1 % Triton and 10 % fetal bovine serum for 1 h at room temperature. For immunostaining, cells were incubated for 2 h at room temperature with 20 μg/mL rabbit anti-KCa1.1 antibody (Sigma) in blocking solution. After three washes with PBS, cells were incubated for 1 h at room temperature with an AlexaFluor 488-labeled goat anti-rabbit antibody (Life Technologies). Nuclear counterstaining was performed with 0.5 μg/mL DAPI. Immunostained cells were imaged with an EVOS fluorescence cell imaging system (Life Technologies).
Colony formation assay
Cells were plated in triplicate in 96-well plates at 2500 cells/well and transferred to 6-well plates 24 h post transfection. After incubation for 10–14 d at 37 °C, cells were fixed with 70 % ethanol and stained with 0.1 % Crystal Violet. The plates were then de-stained with 2 % SDS and absorbance was measured at 562 nM using a FLUOstar Optima.
MPM cell migration was measured using a scratch (or wound-healing) assay. Briefly, transfected cell were plated in 24-well plates and 24 h post-transfection 10 μg/mL Mitomycin C (Sigma) was added to stop cell division; at the same time a scratch was made using a plastic pipette tip. At 24, 48, 72 h post scratch, microscopic imaging was performed with a 20× objective (Leica). Each experiment was performed in triplicate.
Cytosolic and submembrane Ca2+ levels were estimated by overexpressing in MPM cells either soluble GCaMP5 Ca2+ reporter or membrane targeted LCK-GCaMP5 Ca2+ reporter (KCNMA1), respectively . Cells were co-transfected with DNAs coding for reporters and microRNA mimic or siRNA using Lipofectamine 2000 (Life Technologies). Coverslips with transfected cells were placed into glass bottom Petri dishes (MatTek Corporation, Ashland, MA, USA) in 4 mM K+ buffer (150 mM NaCl, 4 mM KCl, 2 mM MgCl2, 10 mM Glucose, 10 mM HEPES, 2 mM CaCl2). Images of cells were captured using an Eclipse TiE fluorescence microscope (Nikon), Plan Apo VC 60× water-immersion objective (Nikon, numerical aperture 1.2), pE-2 CoolLED excitation light source (CooLED, Yorktwon Height, NY, USA) and NIS-elements software (Version 4.0; Nikon). To analyze depolarization-induced Ca2+ influx, 90 mM K+ buffer (64 mM NaCl, 90 mM KCl, 2 mM MgCl2, 10 mM Glucose, 10 mM HEPES, 2 mM CaCl2) was added to cells transfected with LCK-GCaMP5. Images of transfected cells were captured every 0.5 s during the treatment. Fluorescence intensity of the reporter signals was quantified in manually outlined cells using ImageJ (National Institutes of Health) as described .
Differential microarray expression analysis was performed using GeneSpring v12.0 using unpaired t-tests and candidates selected on the basis of statistical significance (P < 0.05) as previously described . Correlation analyses were performed using Pearson-correlation tests and the average expression across all microarray probes for each candidate gene. Group comparisons, correlations and associations were performed using SPSS statistical software and two tailed Mann-Whitney U-tests. A P-value less than 0.05 was considered statistically significant. Pathway enrichment analyses were based on KEGG pathways (Partek Genome Suite v6.4).
ALCAM, activated leukocyte cell adhesion molecule, [GenBank:214, Ensembl:ENSG00000170017; HPRD:03389; MIM:601662; Vega:OTTHUMG00000159192]; CRI1, EP300 interacting inhibitor of differentiation 1, [GenBank:23741, Ensembl:ENSG00000255302; HPRD:06907; MIM:605894; Vega:OTTHUMG00000165911]; FFPE, formalin-fixed paraffin embedded; GSEA, Gene Set Enrichment Analysis; GSN, gelsolin, [Genbank:2934, Ensembl:ENSG00000148180; HPRD:00674; MIM:137350; Vega:OTTHUMG00000020584,]; KCa1.1, calcium-activated potassium channel 1, subunit alpha 1 protein; KCNMA1, calcium-activated potassium channel 1, subunit alpha 1 gene [Gene ID:3778, Ensembl:ENSG00000156113; HPRD:15967; MIM:600150; Vega:OTTHUMG00000018543]; KEGG, Kyoto Encyclopedia of Genes and Genomes; MMP14, matrix metallopeptidase 14 (membrane inserted), [GenBank:4323, Ensembl:ENSG00000157227; HPRD:02856; MIM:600754; Vega:OTTHUMG00000028704]; MPM, Malignant Pleural Mesothelioma; NME2, NME/NM23 nucleoside diphosphate kinase 2, [GenBank: 4831, Ensembl:ENSG00000011052; Ensembl:ENSG00000243678; HPRD:01132; MIM:156491; Vega:OTTHUMG00000154062]; PDGFC, platelet derived growth factor C, [Genbank:56034, HPRD:10529; MIM:608452]
We would like to thank Rebecca Hyland and Kan Chen for curating clinical patient data and collecting patient samples post-surgery respectively. This study was supported by a grant from the Cancer Institute NSW Program Grant (11/TPG/3-06 to NvZ & GR) and a Biaggio Signorelli Fellowship (CW).
YYC, CMW, MBK, MW, KHS carried out experiments and interpreted data. GR, RCYL, YYC and CMW designed experiments and interpreted data. CMW and RCYL did bioinformatics analysis. JJE, MPV, BCM, SK and NvZ provided clinical samples from surgery and clinical data interpretation. SK determined pathology and tumor content of clinical samples. VS and IL conducted calcium flux experiments and generated Fig. 3g and h. CMW, YYC and RCYL generated figures and tables. CMW, GR and RCYL wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kao SC, Reid G, Lee K, Vardy J, Clarke S, van Zandwijk N. Malignant mesothelioma. Intern Med J. 2010;40:742–50.View ArticlePubMedGoogle Scholar
- Tsao AS, Wistuba I, Roth JA, Kindler HL. Malignant pleural mesothelioma. J Clin Oncol. 2009;27:2081–90.View ArticlePubMedGoogle Scholar
- Scherpereel A, Astoul P, Baas P, Berghmans T, Clayson H, de Vuyst P, Dienemann H, Galateau-Salle F, Hennequin C, Hillerdal G, et al. Guidelines of the European Respiratory Society and the European Society of Thoracic Surgeons for the management of malignant pleural mesothelioma. Eur Respir J. 2010;35:479–95.View ArticlePubMedGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Crispi S, Calogero RA, Santini M, Mellone P, Vincenzi B, Citro G, Vicidomini G, Fasano S, Meccariello R, Cobellis G, et al. Global gene expression profiling of human pleural mesotheliomas: identification of matrix metalloproteinase 14 (MMP-14) as potential tumour target. PLoS One. 2009;4:e7016.View ArticlePubMedPubMed CentralGoogle Scholar
- Ishiguro F, Murakami H, Mizuno T, Fujii M, Kondo Y, Usami N, Yokoi K, Osada H, Sekido Y. Activated leukocyte cell-adhesion molecule (ALCAM) promotes malignant phenotypes of malignant mesothelioma. J Thorac Oncol. 2012;7:890–9.View ArticlePubMedGoogle Scholar
- Gordon GJ, Rockwell GN, Jensen RV, Rheinwald JG, Glickman JN, Aronson JP, et al. Identification of novel candidate oncogenes and tumor suppressors in malignant pleural mesothelioma using large-scale transcriptional profiling. Am J Pathol. 2005;166:1827–40.Google Scholar
- Singhal S, Wiewrodt R, Malden LD, Amin KM, Matzie K, Friedberg J, Kucharczuk JC, Litzky LA, Johnson SW, Kaiser LR, Albelda SM. Gene expression profiling of malignant mesothelioma. Clin Cancer Res. 2003;9:3080–97.PubMedGoogle Scholar
- Suraokar MB, Nunez MI, Diao L, Chow CW, Kim D, Behrens C, Lin H, Lee S, Raso G, Moran C, et al. Expression profiling stratifies mesothelioma tumors and signifies deregulation of spindle checkpoint pathway and microtubule network with therapeutic implications. Ann Oncol. 2014;25:1184–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Reid G. MicroRNAs in mesothelioma: from tumour suppressors and biomarkers to therapeutic targets. J Thorac Dis. 2015;7:1031–40.PubMedPubMed CentralGoogle Scholar
- Kao SC, Fulham M, Wong K, Cooper W, Brahmbhatt H, MacDiarmid J, Pattison S, Sagong JO, Huynh Y, Leslie F, et al. A significant metabolic and radiological response after a novel targeted MicroRNA-based treatment approach in malignant pleural mesothelioma. Am J Respir Crit Care Med. 2015;191:1467–9.View ArticlePubMedGoogle Scholar
- Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Garzon R, Calin GA, Croce CM. MicroRNAs in cancer. Annu Rev Med. 2009;60:167–79.View ArticlePubMedGoogle Scholar
- Lee YS, Dutta A. MicroRNAs in cancer. Annu Rev Pathol. 2009;4:199–227.View ArticlePubMedPubMed CentralGoogle Scholar
- Lin SC, Wang CC, Wu MH, Yang SH, Li YH, Tsai SJ. Hypoxia-induced microRNA-20a expression increases ERK phosphorylation and angiogenic gene expression in endometriotic stromal cells. J Clin Endocrinol Metab. 2012;97:E1515–23.View ArticlePubMedGoogle Scholar
- Hossain A, Kuo MT, Saunders GF. Mir-17-5p regulates breast cancer cell proliferation by inhibiting translation of AIB1 mRNA. Mol Cell Biol. 2006;26:8191–201.View ArticlePubMedPubMed CentralGoogle Scholar
- Kirschner MB, Cheng YY, Armstrong NJ, Lin RC, Kao SC, Linton A, Klebe S, McCaughan BC, van Zandwijk N, Reid G. MiR-score: a novel 6-microRNA signature that predicts survival outcomes in patients with malignant pleural mesothelioma. Mol Oncol. 2015;9:715–26.View ArticlePubMedGoogle Scholar
- Reid G, Pel ME, Kirschner MB, Cheng YY, Mugridge N, Weiss J, Williams M, Wright C, Edelman JJ, Vallely MP, et al. Restoring expression of miR-16: a novel approach to therapy for malignant pleural mesothelioma. Ann Oncol. 2013;24:3128–35.View ArticlePubMedGoogle Scholar
- Wright CM, Kirschner MB, Cheng YY, O'Byrne KJ, Gray SG, Schelch K, Hoda MA, Klebe S, McCaughan B, van Zandwijk N, Reid G. Long non coding RNAs (lncRNAs) are dysregulated in Malignant Pleural Mesothelioma (MPM). PLoS One. 2013;8:e70940.View ArticlePubMedPubMed CentralGoogle Scholar
- Williams M, Kirschner MB, Cheng YY, Hanh J, Weiss J, Mugridge N, Wright CM, Linton A, Kao SC, Edelman JJ, et al. miR-193a-3p is a potential tumor suppressor in malignant pleural mesothelioma. Oncotarget. 2015;6(27):23480–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Mogilyansky E, Rigoutsos I. The miR-17/92 cluster: a comprehensive update on its genomics, genetics, functions and increasingly important and numerous roles in health and disease. Cell Death Differ. 2013;20:1603–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang L, Huang J, Yang N, Greshock J, Megraw MS, Giannakakis A, Liang S, Naylor TL, Barchetti A, Ward MR, et al. microRNAs exhibit high frequency genomic alterations in human cancer. Proc Natl Acad Sci U S A. 2006;103:9136–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Grillari J, Hackl M, Grillari-Voglauer R. miR-17-92 cluster: ups and downs in cancer and aging. Biogerontology. 2010;11:501–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Baraniskin A, Birkenkamp-Demtroder K, Maghnouj A, Zollner H, Munding J, Klein-Scory S, Reinacher-Schick A, Schwarte-Waldhoff I, Schmiegel W, Hahn SA. MiR-30a-5p suppresses tumor growth in colon carcinoma by targeting DTL. Carcinogenesis. 2012;33:732–9.View ArticlePubMedGoogle Scholar
- Li X, Zhang Y, Zhang Y, Ding J, Wu K, Fan D. Survival prediction of gastric cancer by a seven-microRNA signature. Gut. 2010;59:579–85.View ArticlePubMedGoogle Scholar
- Kumarswamy R, Mudduluru G, Ceppi P, Muppala S, Kozlowski M, Niklinski J, et al. MicroRNA-30a inhibits epithelial-to-mesenchymal transition by targeting Snai1 and is downregulated in non-small cell lung cancer. Int J Cancer. 2012;130:2044–53.Google Scholar
- Visone R, Pallante P, Vecchione A, Cirombella R, Ferracin M, Ferraro A, Volinia S, Coluzzi S, Leone V, Borbone E, et al. Specific microRNAs are downregulated in human thyroid anaplastic carcinomas. Oncogene. 2007;26:7590–5.View ArticlePubMedGoogle Scholar
- Cheng CW, Wang HW, Chang CW, Chu HW, Chen CY, Yu JC, Chao JI, Liu HF, Ding SL, Shen CY. MicroRNA-30a inhibits cell migration and invasion by downregulating vimentin expression and is a potential prognostic marker in breast cancer. Breast Cancer Res Treat. 2012;134:1081–93.View ArticlePubMedGoogle Scholar
- Wang W, Lin H, Zhou L, Zhu Q, Gao S, Xie H, Liu Z, Xu Z, Wei J, Huang X, Zheng S. MicroRNA-30a-3p inhibits tumor proliferation, invasiveness and metastasis and is downregulated in hepatocellular carcinoma. Eur J Surg Oncol. 2014;40:1586–94.View ArticlePubMedGoogle Scholar
- Zhou S, Liu L, Li H, Eilers G, Kuang Y, Shi S, et al. Multipoint targeting of the PI3K/mTOR pathway in mesothelioma. Br J Cancer. 2014;110:2479–88.Google Scholar
- Kawaguchi K, Murakami H, Taniguchi T, Fujii M, Kawata S, Fukui T, Kondo Y, Osada H, Usami N, Yokoi K, et al. Combined inhibition of MET and EGFR suppresses proliferation of malignant mesothelioma cells. Carcinogenesis. 2009;30:1097–105.View ArticlePubMedGoogle Scholar
- Shapiro IM, Kolev VN, Vidal CM, Kadariya Y, Ring JE, Wright Q, Weaver DT, Menges C, Padval M, McClatchey AI, et al. Merlin deficiency predicts FAK inhibitor sensitivity: a synthetic lethal relationship. Sci Transl Med. 2014;6:237ra268.Google Scholar
- Carbone M, Yang H. Molecular pathways: targeting mechanisms of asbestos and erionite carcinogenesis in mesothelioma. Clin Cancer Res. 2012;18:598–604.View ArticlePubMedPubMed CentralGoogle Scholar
- Stevenson JP, Kindler HL, Papasavvas E, Sun J, Jacobs-Small M, Hull J, Schwed D, Ranganathan A, Newick K, Heitjan DF, et al. Immunological effects of the TGFbeta-blocking antibody GC1008 in malignant pleural mesothelioma patients. Oncoimmunology. 2013;2:e26218.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang J, Yan J. Regulation of BK channels by auxiliary gamma subunits. Front Physiol. 2014;5:401.PubMedPubMed CentralGoogle Scholar
- Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, et al. A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A. 2006;103:2257–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Bloch M, Ousingsawat J, Simon R, Schraml P, Gasser TC, Mihatsch MJ, Kunzelmann K, Bubendorf L. KCNMA1 gene amplification promotes tumor cell proliferation in human prostate cancer. Oncogene. 2007;26:2525–34.View ArticlePubMedGoogle Scholar
- Oeggerli M, Tian Y, Ruiz C, Wijker B, Sauter G, Obermann E, Guth U, Zlobec I, Sausbier M, Kunzelmann K, Bubendorf L. Role of KCNMA1 in breast cancer. PLoS One. 2012;7:e41664.View ArticlePubMedPubMed CentralGoogle Scholar
- Mazar J, DeYoung K, Khaitan D, Meister E, Almodovar A, Goydos J, Ray A, Perera RJ. The regulation of miRNA-211 expression and its role in melanoma cell invasiveness. PLoS One. 2010;5:e13779.View ArticlePubMedPubMed CentralGoogle Scholar
- Rabjerg M, Olivan-Viguera A, Hansen LK, Jensen L, Sevelsted-Moller L, Walter S, Jensen BL, Marcussen N, Kohler R. High expression of KCa3.1 in patients with clear cell renal carcinoma predicts high metastatic risk and poor survival. PLoS One. 2015;10:e0122992.View ArticlePubMedPubMed CentralGoogle Scholar
- Pardo LA, Stuhmer W. The roles of K(+) channels in cancer. Nat Rev Cancer. 2014;14:39–48.View ArticlePubMedGoogle Scholar
- Pardo LA, del Camino D, Sanchez A, Alves F, Bruggemann A, Beckh S, Stuhmer W. Oncogenic potential of EAG K(+) channels. EMBO J. 1999;18:5540–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Bai Y, Liao H, Liu T, Zeng X, Xiao F, Luo L, Guo H, Guo L. MiR-296-3p regulates cell growth and multi-drug resistance of human glioblastoma by targeting ether-a-go-go (EAG1). Eur J Cancer. 2013;49:710–24.View ArticlePubMedGoogle Scholar
- Wu X, Zhong D, Gao Q, Zhai W, Ding Z, Wu J. MicroRNA-34a inhibits human osteosarcoma proliferation by downregulating ether a go-go 1 expression. Int J Med Sci. 2013;10:676–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Feng J, Yu J, Pan X, Li Z, Chen Z, Zhang W, Wang B, Yang L, Xu H, Zhang G, Xu Z. HERG1 functions as an oncogene in pancreatic cancer and is downregulated by miR-96. Oncotarget. 2014;5:5832–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Fan M, Krutilina R, Sun J, Sethuraman A, Yang CH, Wu ZH, Yue J, Pfeffer LM. Comprehensive analysis of microRNA (miRNA) targets in breast cancer cells. J Biol Chem. 2013;288:27480–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Birnie KA, Yip YY, Ng DC, Kirschner MB, Reid G, Prele CM, Musk AW, Lee YC, Thompson PJ, Mutsaers SE, Badrian B. Loss of miR-223 and JNK signaling contribute to elevated Stathmin in malignant pleural mesothelioma. Mol Cancer Res. 2015;13:1106–18.View ArticlePubMedGoogle Scholar
- Khaitan D, Sankpal UT, Weksler B, Meister EA, Romero IA, Couraud PO, Ningaraj NS. Role of KCNMA1 gene in breast cancer invasion and metastasis to brain. BMC Cancer. 2009;9:258.View ArticlePubMedPubMed CentralGoogle Scholar
- Weaver AK, Bomben VC, Sontheimer H. Expression and function of calcium-activated potassium channels in human glioma cells. Glia. 2006;54:223–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Bose T, Cieslar-Pobuda A, Wiechec E. Role of ion channels in regulating Ca2+ homeostasis during the interplay between immune and cancer cells. Cell Death Dis. 2015;6:e1648.View ArticlePubMedPubMed CentralGoogle Scholar
- Deliot N, Constantin B. Plasma membrane calcium channels in cancer: Alterations and consequences for cell proliferation and migration. Biochim Biophys Acta. 1848;2015:2512–22.Google Scholar
- Huang X, Jan LY. Targeting potassium channels in cancer. J Cell Biol. 2014;206:151–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Bentzen BH, Olesen SP, Ronn LC, Grunnet M. BK channel activators and their therapeutic perspectives. Front Physiol. 2014;5:389.View ArticlePubMedPubMed CentralGoogle Scholar
- Cui J, Yang H, Lee US. Molecular mechanisms of BK channel activation. Cell Mol Life Sci. 2009;66:852–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Patergnani S, Giorgi C, Maniero S, Missiroli S, Maniscalco P, Bononi I, Martini F, Cavallesco G, Tognon M, Pinton P. The endoplasmic reticulum mitochondrial calcium cross talk is downregulated in malignant pleural mesothelioma cells and plays a critical role in apoptosis inhibition. Oncotarget. 2015;6:23427–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang R, Tian P, Chi Q, Wang J, Wang Y, Sun L, Liu Y, Tian S, Zhang Q. Human ether-a-go-go-related gene expression is essential for cisplatin to induce apoptosis in human gastric cancer. Oncol Rep. 2012;27:433–40.PubMedGoogle Scholar
- Marklund L, Henriksson R, Grankvist K. Cisplatin-induced apoptosis of mesothelioma cells is affected by potassium ion flux modulator amphotericin B and bumetanide. Int J Cancer. 2001;93:577–83.View ArticlePubMedGoogle Scholar
- Chantome A, Potier-Cartereau M, Clarysse L, Fromont G, Marionneau-Lambot S, Gueguinou M, et al. Pivotal role of the lipid Raft SK3-Orai1 complex in human cancer cell migration and bone metastases. Cancer Res. 2013;73:4852–61.Google Scholar
- Steinle M, Palme D, Misovic M, Rudner J, Dittmann K, Lukowski R, Ruth P, Huber SM. Ionizing radiation induces migration of glioblastoma cells by activating BK K(+) channels. Radiother Oncol. 2011;101:122–6.View ArticlePubMedGoogle Scholar
- Relan V, Morrison L, Parsonson K, Clarke BE, Duhig EE, Windsor MN, Matar KS, Naidoo R, Passmore L, McCaul E, et al. Phenotypes and karyotypes of human malignant mesothelioma cell lines. PLoS One. 2013;8:e58132.View ArticlePubMedPubMed CentralGoogle Scholar
- Smythe WR, Kaiser LR, Hwang HC, Amin KM, Pilewski JM, Eck SJ, Wilson JM, Albelda SM. Successful adenovirus-mediated gene transfer in an in vivo model of human malignant mesothelioma. Ann Thorac Surg. 1994;57:1395–401.View ArticlePubMedGoogle Scholar
- Kao SC, Lee K, Armstrong NJ, Clarke S, Vardy J, van Zandwijk N, Reid G, Burn J, McCaughan BC, Henderson DW, Klebe S. Validation of tissue microarray technology in malignant pleural mesothelioma. Pathology. 2011;43:128–32.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Kirschner MB, Cheng YY, Badrian B, Kao SC, Creaney J, Edelman JJ, Armstrong NJ, Vallely MP, Musk AW, Robinson BW, et al. Increased circulating miR-625-3p: a potential biomarker for patients with malignant pleural mesothelioma. J Thorac Oncol. 2012;7:1184–91.View ArticlePubMedGoogle Scholar
- Reid G, Wallant NC, Patel R, Antonic A, Saxon-Aliifaalogo F, Cao H, Webster G, Watson JD. Potent subunit-specific effects on cell growth and drug sensitivity from optimised siRNA-mediated silencing of ribonucleotide reductase. J RNAi Gene Silenc. 2009;5:321–30.Google Scholar
- Akerboom J, Chen TW, Wardill TJ, Tian L, Marvin JS, Mutlu S, Calderon NC, Esposti F, Borghuis BG, Sun XR, et al. Optimization of a GCaMP calcium indicator for neural activity imaging. J Neurosci. 2012;32:13819–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Sheng L, Leshchyns'ka I, Sytnyk V. Neural cell adhesion molecule 2 promotes the formation of filopodia and neurite branching by inducing submembrane increases in Ca2+ levels. J Neurosci. 2015;35:1739–52.View ArticlePubMedGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4:44–57.View ArticleGoogle Scholar