- Primary research
- Open Access
Contribution of voltage-gated sodium channel β-subunits to cervical cancer cells metastatic behavior
Cancer Cell Internationalvolume 19, Article number: 35 (2019)
Voltage-gated sodium (NaV) channels are heteromeric proteins consisting of a single pore forming α-subunit associated with one or two auxiliary β-subunits. These channels are classically known for being responsible of action potential generation and propagation in excitable cells; but lately they have been reported as widely expressed and regulated in several human cancer types. We have previously demonstrated the overexpression of NaV1.6 channel in cervical cancer (CeCa) biopsies and primary cultures, and its contribution to cell migration and invasiveness. Here, we investigated the expression of NaV channels β-subunits (NaVβs) in the CeCa cell lines HeLa, SiHa and CaSki, and determined their contribution to cell proliferation, migration and invasiveness.
We assessed the expression of NaVβs in CeCa cell lines by performing RT-PCR and western blotting experiments. We also evaluated CeCa cell lines proliferation, migration, and invasion by in vitro assays, both in basal conditions and after inducing changes in NaVβs levels by transfecting specific cDNAs or siRNAs. The potential role of NaVβs in modulating the expression of NaV α-subunits in the plasma membrane of CeCa cells was examined by the patch-clamp whole-cell technique. Furthermore, we investigated the role of NaVβ1 on cell cycle in SiHa cells by flow cytometry.
We found that the four NaVβs are expressed in the three CeCa cell lines, even in the absence of functional NaV α-subunit expression in the plasma membrane. Functional in vitro assays showed differential roles for NaVβ1 and NaVβ4, the latter as a cell invasiveness repressor and the former as a migration abolisher in CeCa cells. In silico analysis of NaVβ4 expression in cervical tissues corroborated the downregulation of this protein expression in CeCa vs normal cervix, supporting the evidence of NaVβ4’s role as a cell invasiveness repressor.
Our results contribute to the recent conception about NaVβs as multifunctional proteins involved in cell processes like ion channel regulation, cell adhesion and motility, and even in metastatic cell behaviors. These non-canonical functions of NaVβs are independent of the presence of functional NaV α-subunits in the plasma membrane and might represent a new therapeutic target for the treatment of cervical cancer.
Cervical cancer (CeCa) is the fourth most common cancer in women worldwide, with more than 500,000 new cases each year and an estimated 260,000 deaths, of which approximately 87% occur in developing countries . Despite the implementation of vaccination against Human Papillomavirus infection, it will take decades before the impact on the cervical cancer incidence became clear, giving the natural history of this disease . Therefore, it is necessary to better understand the biology of this type of cancer in order to improve the early detection and treatment.
Voltage-gated sodium (NaV) channels are heteromeric transmembrane proteins that consist of a single pore forming α-subunit associated with one or more auxiliary β-subunits (NaVβs) . These channels are known to be responsible for action potential generation and propagation in excitable cells, but they are also widely expressed and upregulated in a variety of human cancer types, including breast, prostate, colon, ovarian, lung, gastric, melanoma, astrocytoma and uterine cervix, where the NaV α-subunit activity has been associated mainly, but not exclusively, with cell motility and invasiveness [4,5,6,7,8,9,10,11,12,13,14]. So far, four NaVβs (NaVβ1–NaVβ4) have been described. These are multifunctional type I transmembrane proteins that modulate NaV α-subunit gating, localization and transit to the plasma membrane. Additionally, they possess a V-type immunoglobulin repeat in the extracellular domain similar to the family of neural cell adhesion molecules (CAMs), thus, their role as CAMs in a trans-homophilic and trans-heterophilic way, and with other molecules, in presence or absence of the α-subunits has been described in numerous reports and in different expression systems [15,16,17,18,19,20,21,22].
The role of NaVβs in cancer has also been studied, although not as widely as for the α-subunits. In prostate cancer (PCa), it has been found that the overexpression of NaVβs is associated with the increasing metastatic potential of PCa cell lines , and when overexpressing NaVβ2 in LnCaP cells, the cells acquired a fibroblastic-like morphology, became more invasive and had enhanced migration in vitro; in contrast, these cells had a reduced tumor volume in vivo [24, 25]. On the other hand, the overexpression of NaVβ3 in Saos-2 (osteosarcoma) and T98G (glioblastoma) cells, seems to activate an apoptotic pathway mediated by the tumor suppressor protein p53 . However, most of the studies have been done in breast cancer, where it has been demonstrated that NaVβ1 is overexpressed in the low metastatic MCF-7 cell line vs highly metastatic MDA-MB-231 cells . Overexpressing NaVβ1 in these last cells produced a drastic reduction on in vitro transwell migration and an increase in cellular adhesiveness. In contrast, MDA-MB-231 cells stably transfected with NaVβ1 produced larger tumors in an in vivo model compared with those generated by wild type MDA-MB-231 cells. In addition, overexpression of NaVβ1 induced a reduction in caspase activity in these cells . More recently, the potential role of NaVβ4 as a metastasis-suppressor gene has been described in breast cancer cells, associated with the small GTPase RhoA activity and therefore with the cytoskeleton remodeling during migration and invasiveness .
We have previously demonstrated the overexpression of NaV α-subunits in cervical cancer biopsies vs non-cancerous cervix, and furthermore, we reported the functional expression of these proteins in primary cultures derived from human cervical cancer biopsies and their contribution to the cell migration and invasiveness. Specifically, NaV1.6 showed to be the most overexpressed α-subunit and the one responsible of the augmented invasive potential as the specific blockade of this channel was able to diminish the invasive potential with the same effectiveness as the blockade of all the NaV α-subunits by TTX [14, 30]. In the present work, we investigated the electrophysiological activity of NaV channels and the expression of NaVβs in three cervical cancer cell lines (HeLa, SiHa and CaSki) as well as their contribution to three metastatic cell behaviors: in vitro proliferation, transwell migration and invasiveness.
Cell lines and cell culture
Cervical cancer cell lines SiHa, CaSki (both HPV-16 positive) and HeLa (HPV-18 positive) were grown in DMEM supplemented with 10% fetal bovine serum (FBS), 100 U/ml penicillin, and 100 μg/ml streptomycin at 37 °C in a humidified 5% CO2 incubator. All cell culture reagents were purchased from Gibco-Thermo Fisher Scientific (Waltham, MA). The three cell lines were authenticated by short tandem repeat profiling prior to the experiments (data not shown).
Overexpression and downregulation of NaVβs expression
Plasmids containing each of the NaVβs from rat were kindly donated by Dr. L. Isom (University of Michigan). Overexpression of each NaVβ in CeCa cells was achieved by transient transfections using the JetPEI transfection reagent (PolyPlus Transfection; Illkirch, France), according to the manufacturer’s instructions. On the other hand, downregulation of NaVβs expression was performed using two or three different predesigned siRNAs (Sigma-Aldrich; St. Louis, MO) for each NaVβ (100 nM each); in all cases, we chose the one with the best efficiency. Optimal siRNA transfection conditions were stablished previously by using 50 nM of siGLO Green Transfection Indicator (Sigma-Aldrich), a fluorescent siRNA for determining transfection efficiency. Efficiency of both experimental strategies to manipulate the expression levels of each NaVβ was evaluated by RT-PCR and western blot (see Additional file 1).
Detailed methods for whole-cell patch-clamp recordings and protocols have been previously described [14, 30]. Sodium currents of native or transiently transfected CeCa cells were recorded at room temperature (20–23 °C) using an Axopatch 200B amplifier, a Digidata 1320A/D converter, and the pCLAMP 10.0 software (Molecular Devices; Sunnyvale, CA). As a positive control for sodium currents, we used transiently transfected CeCa cells with the NaV1.6 channel.
Standard PCR (RT-PCR)
The NaVβs mRNA expression levels of CeCa cell lines were assessed by RT-PCR with specific primers for each β-subunit and β-actin as control (Table 1) as described before . Briefly, total RNA of each cell line was isolated with Trizol reagent (Invitrogen; Carlsbad, CA) and RT-PCR reactions were performed in a final volume of 25 µl using the Super-Script One-step RT-PCR kit (Invitrogen) with 250 ng of total RNA, 0.5 µl of enzyme mix, 0.2–0.5 µM of each primer, and 12.5 µl of 2X buffer containing 0.4 mM of each dNTP and 2.4 mM MgSO4.
Real-time PCR (qPCR)
Total RNA was extracted using the RNeasy Mini Kit (Qiagen; Hilden, Germany), then reverse-transcribed with the High Capacity cDNA Reverse Transcription kit (Applied Biosystems; Foster City, CA) according to the manufacturer’s instructions using 2 µg of total RNA in a final volume of 20 µl. Real-time PCR was carried out in a Rotor-Gene Q (Qiagen) using Custom TaqMan Gene Expression Assays (Applied Biosystems) as described before . Briefly, 100 ng of cDNA, 0.4 µl of the TaqMan assay (Table 2) and 5 µl of TaqMan Universal PCR Master Mix (Applied Biosystems) were mixed in a final reaction volume of 10 µl for each qPCR reaction. At least three independent experiments were done, and each assay was performed in triplicate. The results were analyzed by the 2−ΔΔCt method  using HPRT1 expression as the normalizing gene control and results are shown as relative expression values of NaVβ1 in HeLa cells.
Total protein from native or transiently transfected CeCa cells was extracted 24, 48, 72 and 96 h post-transfection (with cDNA or siRNAs, for overexpression or inhibition of the NaVβ expression respectively) using RIPA buffer (25 mM Tris–HCl, pH 7.4; 150 mM NaCl; 1% IGEPAL; 1% Sodium deoxycholate, and 1% SDS) supplemented with complete EDTA-free protease inhibitors (Roche, Switzerland), and quantified by Bradford assay. Equal amounts of protein (100 µg) were subjected to SDS-PAGE, transferred onto a polyvinylidene difluoride membrane (Millipore, Billerica MA) and probed overnight with the following primary antibodies: rabbit anti-NaVβ1 (1:3000; LifeSpan BioSciences Inc.; Seattle, WA); rabbit anti-NaVβ2 (1:1000; LifeSpan BioSciences Inc.); rabbit anti-NaVβ3 (1:5000; Abcam; Cambridge, UK), rabbit anti-NaVβ4 (1:3000; Novus Biologicals; Littleton, CO) and a homemade mouse anti-β-actin antibody (1:1000) used as a loading control. Blots were subsequently probed with an anti-rabbit or an anti-mouse (as the case may be) secondary antibody conjugated with horseradish peroxidase (1:10,000; Santa Cruz Biotechnology; Dallas, TX) for 1 h at room temperature and visualized using the SuperSignal West Pico chemiluminescent substrate (Thermo Fisher Scientific). Signal intensity of immunoblots was calculated with ImageJ software.
Cell proliferation was determined using the colorimetric 3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide (MTT) assay. Briefly, cells were plated in 48-well microplates at an initial density of 2 × 103 cells per well. To measure cell proliferation at 24, 48, 72 and 96 h after seeding, 10 µl of MTT solution (5 mg/ml) were added in each well after replacing the culture medium with 200 µl of fresh culture medium. The cells were incubated at 37 °C, and 3 h later, cell culture and MTT solution were removed and 150 µl of dimethyl sulfoxide was added into the plate to dissolve formazan crystals. The absorbance at a wavelength of 562 nm was recorded using an Eppendorf AG 22331 Biophotometer (Eppendorf; Hamburg, Germany). For analyzing the role of NaVβs on proliferation, after overnight incubation, cells were transfected with the respective cDNA to overexpress or the specific siRNA to downregulate each β-subunit expression, and finally, MTT assay was performed 24, 48, 72 and 96 h after transfection.
Cell cycle analysis by flow cytometry
The role of NaVβ1 on cell cycle in SiHa cells was assessed by flow cytometry after staining the cells with propidium iodide (PI; Molecular Probes, Life Technologies; Grand Island, NY). Briefly, 2 × 105 cells were seeded onto 12-well plates, incubated in normal conditions for 24 h, synchronized by serum-reduced conditions (0.2% FBS) for 48 h and then transfected with the specific cDNA or siRNA. After incubating for 48 h, cells were collected and fixed with ice-cold 70% ethanol, stored for at least 24 h at − 20 °C, treated with RNase and stained with PI for 1 h at 37 °C. For each condition, 2 × 104 cells (singlets) were analyzed in an Accuri C6 cytometer (BD Biosciences; San Jose, CA), and the percentages of cells in each cell cycle phase were analyzed with FlowJo software (Tree Star; Ashland, OR) using the Dean-Jett-Fox algorithm by three independent tests for each condition. Appropriate gating was used to select the single cell population and used on all samples.
In vitro migration and invasion assays
Migration and invasion assays were done using culture inserts with 8-µm pore size membranes covered or not with Matrigel (Corning Inc., Corning, NY) for invasion and migration respectively. The upper compartment was seeded with 5 × 104 for migration and 1 × 105 cells for invasion in DMEM supplemented with 5% FBS, and the lower compartment was filled with DMEM supplemented with 15% FBS as a chemoattractant. The number of cells migrating or invading (according to the circumstances) over 48 h at 37 °C was evaluated using the MTT assay described before. Results were compiled as the mean of at least three repeats, each time done by duplicate, and compared with control conditions.
In silico immunohistochemistry (IHC) staining analysis
Expression of NaVβ4 in CeCa and normal cervical tissues was examined using the Human Protein Atlas (http://www.proteinatlas.org/) [32, 33], a Swedish-based online tool for human protein expression analysis. Representative CeCa and normal tissues cores were chosen for illustrative purposes.
All quantitative results are given as the mean ± S.D. Differences in means were tested with an unpaired two-tailed Student’s t test.
Cervical cancer cell lines do not express plasma-membrane voltage-gated sodium currents
Our previous results showed the functional expression of NaV channels in CeCa primary cultures and specifically demonstrated the relevance of the NaV1.6 channel in the invasiveness of cervical cancer [14, 30]. Thus, in order to further investigate the role of NaV channels in the biology of CeCa, we intended to use CeCa cell lines for such experiments. However, we found that SiHa, CaSki and HeLa cell lines do not express functional NaV channels at the plasma-membrane as determined by whole-cell patch-clamp experiments. (n = 30 cells total; see Additional file 2), even though we have previously shown that they do express the NaV1.6 mRNA and at least one isoform of the NaV1.6 protein .
Endogenous expression of NaVβs in cervical cancer cell lines
In previous studies, we described the role of NaV α-subunits in the invasive potential of CeCa primary cultures. In this work, we wanted to analyze the possible role of NaVβs in different CeCa cellular processes. To that end, endogenous expression of NaVβs was first evaluated by RT-PCR using specific primers for each subunit. In most of cases, the mRNA was abundant as evidenced by the generation of the expected amplicon for each β-subunit (Fig. 1a), suggesting the very likely presence of these proteins. Along with the expected bands, additional PCR amplicons were observed when analyzing the expression of NaVβ2 and NaVβ4 by RT-PCR (Fig. 1a). The ≈ 150 bp DNA fragment observed in the NaVβ4 line matches with a splice variant that has been already reported , that would produce a highly truncated protein (if translated) with an unknown functionality. The larger non-specific bands in NaVβ2 and NaVβ4 RT-PCR do not correspond to any reported variants of such isoforms, and we cannot rule out the possibility of a non-specific amplification with these primers. A more detailed study is required to analyze the NaVβs isoforms expressed in the CeCa cells and their potential function and/or contribution to the cancerous phenotype. Nevertheless, to carry out a more quantitative analysis of the NaVβs expression levels in the different CeCa cell lines we performed qPCR experiments. Results exhibited substantial differences in NaVβs mRNA expression, but in all cases, NaVβ1 showed to be highly expressed, whereas NaVβ3 displayed the lowest expression levels, except in SiHa cells (Fig. 1b). In addition, evaluation of protein expression by western blot for each NaVβ indicates that the four NaVβs are expressed in the three CeCa cell lines (Fig. 1c). Interestingly, SiHa was the only cell line that express very similar levels of all four NaVβs (Fig. 1d). Also, SiHa and CaSki cells expressed significantly higher levels of NaVβ1 and NaVβ4 than HeLa cells. While NaVβ3 levels were similar among the CeCa cell lines, and NaVβ2 was the lees expressed except in SiHa cells (Fig. 1d). The western blot results suggest the absence of additional variants of NaVβs isoforms, as we only observed the protein of the expected size.
Proliferation, migration and invasion of CeCa cell lines under basal conditions
In order to evaluate the metastatic potential of each CeCa cell line, we analyzed three key metastatic behaviors: proliferation, migration and invasiveness in vitro. HeLa cells demonstrated to be the most aggressive cell line as these cells proliferated, migrated and invaded the most, followed by SiHa cells, and CaSki cells (Fig. 2).
Increasing or decreasing the expression of NaVβs is not enough to promote the functional activity of α-subunits in the plasma membrane
First, to evaluate the siRNA transfection efficacy, we used the siGLO Green Transfection Indicator, a fluorescent dye attached to a scramble siRNA, and analyzed the cells 24 h later by epifluorescent microscopy. The results indicated that almost 100% of the cells were transfected with the siRNAs, as can be appreciated in Additional file 3. Thus, we followed the same conditions to transfect the specific siRNAs against each of the NaVβs.
It is well known that some of the NaVβs are implicated in the correct trafficking of the NaV α-subunits to the plasma membrane, as has been shown with the heterologous in vitro expression of NaVβ1 and NaVβ2 in Xenopus oocytes [15, 35], and in vivo with NaVβ2 knocked-out mice . These evidences made us wonder if the sole overexpression or even the down regulation of the expression of NaVβs was enough to promote the functional expression of the pore-forming sodium channels in the plasma membrane and therefore the appearance of voltage-gated sodium currents. To address this question, HeLa and SiHa cells were analyzed by the whole-cell patch-clamp technique 24, 48, 72 and 96 h after transfection to increase or decrease the expression of each of the NaVβs. More than 130 cells were successfully investigated with this technique, but we did not find any voltage-gated sodium currents, suggesting that the regulation of the expression of NaVβs is not enough to promote the functional activity of NaV α-subunits in the plasma membrane (see Additional file 2).
Role of NaVβs on proliferation
Next, we investigated whether the expression of NaVβs had an effect on cell proliferation. This was evaluated in CeCa cell lines 24, 48, 72 and 96 h post-transfection with the cDNA or with the specific siRNA for each NaVβ subunit to increase or to reduce the expression respectively. Only NaVβ1 had an effect on proliferation of SiHa cells: the overexpression of NaVβ1 increased cell proliferation, whereas its downregulation decreased cell duplication (Fig. 3a; center panels). This is similar to what was found with the breast cancer cell line MDA-MB-231, where the overexpression of NaVβ1 in these cells led to the development of larger tumors in mice . To corroborate this observation, we decided to perform a cell cycle assay by flow cytometry. First, for standardization purposes, SiHa cells were arrested at the G0/G1 phase by incubation in serum-reduced conditions for 24 and 48 h. To estimate the proportion of cells in different phases of the cell cycle, cells were stained with PI and analyzed by flow cytometry. Results show a consistent arrest of cells when treated for 48 h, as the proportion of the cells in G0/G1 phase increased significantly whereas the percentage of cells in G2/M phase was smaller compared to control (normal medium) conditions (see Additional file 4). Once synchronized cells were transfected with either cDNA or siRNA for NaVβ1 subunit under the experimental conditions indicated above. Results show a significant increase in the proportion of cells in G2/M phase when overexpressing NaVβ1 vs the control condition (JetPEI mock control), as well as a significant decrease in the proportion of cells in G0/G1 phase. In the same sense, the downregulation of the expression of NaVβ1 generated a significantly larger proportion of cells in G0/G1 phase, with no significant changes in the proportion of cells in G2/M phase compared with the N-TER mock control (Fig. 3b). These results suggest a positive role of NaVβ1 in the proliferation of SiHa cells, but not in HeLa or CaSki cells.
Role of NaVβs on migration and invasion
The role of NaVβs on in vitro migration and invasion of CeCa cell lines was evaluated 48 h post-transfection with the cDNA or with the siRNA for each β-subunit to overexpress or to downregulate the expression of NaVβs respectively. Migration results show that overexpression of NaVβs in HeLa cells induces a significant reduction in migration capability (around 50% in all cases), regardless of the type of NaVβ, accompanied by a tendency to increase the invasive capability when the NaVβs are downregulated, being the effect of the treatment with siRNAs against NaVβ1 the most substantial (a 58% increment). In SiHa cells, a similar result was found with NaVβ1 and NaVβ4; a significant reduction on the migration capability (36% and 28% respectively) when they are overexpressed, and an increment (60% and 47% respectively) when they are downregulated, while the manipulation of the expression of NaVβ2 and NaVβ3 did not alter the migration of these cells. Finally, in CaSki cells, only NaVβ1 showed to have a role on cellular migration: the overexpression of this subunit decreased the migration by 14%, whereas the transfection with siRNAs increased it by 41% (Fig. 4a). It is worth of notice that only the effect of NaVβ1 (negative regulation) in migration remains constant and statistically significant in the three cell lines, suggesting a potential role of NaVβ1 in CeCa migration capability.
With regards to the in vitro invasion, results show that the downregulation of NaVβ4 increased the invasiveness of the three cell lines: 64% in HeLa, 43% in SiHa and 63% in CaSki cells. Accordingly, the overexpression of this protein induced a reduction on the invasiveness of SiHa and CaSki cells by 28% and 21% respectively, suggesting a potential role of NaVβ4 as an invasion suppressor. The rest of the NaVβs did not significantly modified the invasion capacity of the CeCa cells (Fig. 4b).
Given all these results, we decided to corroborate whereas the diminished expression of NaVβ4 was associated with the presence of CeCa vs normal cervix tissues. For that purpose, we analyzed in silico the expression NaVβ4 in CeCa biopsies and normal cervical tissues from the Human Protein Atlas. Although the database only has a small number of IHC normal cervix and CeCa biopsies stained for NaVβ4 (n = 3 for normal cervix, and n = 12 for cervical cancer), the relationship between the expression of the protein and the presence or absence of cervical cancer is evident: the staining results indicate that three out of three normal cervix tissues have moderate NaVβ4 expression, while in the CeCa biopsies no signal of NaVβ4 was detected in twelve out of twelve cases (see Additional file 5). Taken together, these results strongly suggest that NaVβ4 expression is downregulated in CeCa, leading to an increase in the invasive capacity of the cells, and thus, the metastatic potential.
For many years, NaVβs were considered only as “auxiliary” NaV subunits commissioned of regulating the α-subunit kinetics and cellular localization. Nowadays plenty of evidence shows the multifunctionality of these proteins, participating in a wide variety of cellular processes like cell adhesion, cytoskeleton remodeling, filopodia and invadopodia formation, transcriptional regulation, cell cycle regulation and even apoptosis [21, 29].
For the present work, we chose three different CeCa cell lines in order to cover a wide range of CeCa cases, and to analyze the potential role of NaVβs in the malignant behavior of this disease. SiHa and CaSki cells are HPV-16 positive (the most frequently detected viral type in CeCa, found in approximately 57% of cases) whereas HeLa cells arise from a HPV-18 positive tumor (the second most frequent virus; 16% of cases) . Furthermore, HeLa and SiHa cells derived from a primary tumor, while CaSki cells were taken from a metastatic tumor. In addition, the patients’ ethnicity and the cellular origins of the three cell lines are different [38,39,40].
Here we demonstrated the presence (expression) of NaVβs in the three different CeCa cell lines (Fig. 1), even in the absence of functional NaV α-subunit expression in the plasma membrane (evidenced by whole-cell patch-clamp recordings, see Additional file 2), suggesting that their role in cancer cell biology is independent of the pore-forming NaV α-subunit.
When we evaluated the aggressiveness of each cell line, by analyzing the basal proliferation, migration and invasiveness in vitro, HeLa cells demonstrated to be the most aggressive cell line, as these cells proliferated, migrated and invaded the most, followed by SiHa cells. CaSki cells were the less proliferative, and they migrated and invaded the less (Fig. 2). As far as we are concerned, this is the first report about the differences in the metastatic behavior among these CeCa cell lines.
We have previously reported the role of NaV1.6 (the most overexpressed NaV α-subunit in CeCa biopsies) in the invasiveness potential of CeCa primary cultures . In agreement with our recent findings , in the present work we did not find voltage-gated sodium currents in any of the CeCa cell lines studied with the patch-clamp technique. This could be due to at least two possibilities: either the sodium channels are being correctly synthetized but are not reaching the plasma membrane for some unknown reason; or these cells are expressing a non-functional splicing variant. In this regard, it has been reported that NaVβs participate in the appropriate traffic and localization of the NaV α-subunits in the plasma membrane [15, 35, 36]. Moreover, it has been shown that NaVβ2 can be cleaved by secretases, generating small intracellular peptides capable of reaching the cell nucleus and promoting the transcriptional upregulation of genes including SCN1A, which encodes for NaV1.1 [41, 42]. Hence, we thought that maybe by modifying the expression of these proteins we could promote the functional localization of conducting NaV channels in the plasma membrane, however our results suggest this is not the case, as the manipulation of the expression of the NaVβs in HeLa and SiHa cells was not enough to promote the appearance of voltage-dependent sodium currents (see Additional file 2). Yet, there are many other molecules such as hormones and growth factors that could be involved in the correct trafficking of the channel. In fact, we have examined the possibility that NaVs functional expression in CeCa cell lines could be regulated by β-estradiol or EGF, but we did not find sodium currents under any of the experimental conditions tested (data not shown). Another likely explanation could be related to the oxygen conditions in which the cell lines are cultured. It has been reported that hypoxia can change both the expression and activity of NaVs. Tumor cells usually experience severe hypoxia, which has been correlated with a more robust expression and activity of NaVs [43, 44]. On the contrary, CeCa cell lines have been cultured in normoxic tissue culture conditions, which might have modified the NaVs expression and/or trafficking.
It is also possible that CeCa cell lines are expressing a non-conducting channel with a different and unknown function and/or localization in the cell. It remains to be investigated whether NaV α-subunits are expressed in intracellular compartments, as well as their possible role in invasion. The existence of functional intracellular sodium channels has been demonstrated in macrophages where the activity of this sodium channel contributes to cellular invasion through a mechanism involving sodium ions release from cationic intracellular stores, followed by a subsequent mitochondrial calcium release mediated by the Na+/Ca2+ exchanger, which in turn facilitates cytoskeletal remodeling and invadopodia formation . In fact, recent data from our group suggests the presence of intracellular NaV1.6 channels in CeCa cell lines . However, more experiments are needed to fully elucidate this discrepancy between NaVs functional expression in CeCa primary cultures vs cell lines.
Regarding to the effect of NaVβs on cell proliferation, we found that only NaVβ1 had an effect on this process and only in SiHa cells: the overexpression of NaVβ1 increased cell proliferation, whereas the downregulation of the expression of this protein decreased it (Fig. 3a). This is in agreement with observations reported by Nelson et al.  demonstrating that NaVβ1 enhances breast tumor growth and metastasis in vivo, by increasing cell proliferation and reducing apoptosis. The fact that we only observed this effect in one out of three CeCa cell lines shows that it is a particularity of the cell line (SiHa) rather than a general role of NaVβ1 on proliferation in cervical cancer.
On the other hand, it has been suggested that NaVβ3 mediates a p53-dependent apoptotic pathway . In breast cancer cells, for example, this subunit is totally absent, although it is normally expressed in non-cancerous breast cells (Sanchez-Sandoval et al., unpublished data). However, the overexpression or down regulation of this subunit did not alter the proliferation rate of the CeCa cell lines analyzed here, suggesting a lack of the pro-apoptotic activity of NaVβ3 in these cells. This might be related to the p53 protein status in CeCa cell lines. It has been reported that the association between the E6 protein (one of the two main HPV oncogenes) with p53 leads to the specific ubiquitination and degradation of p53 protein , therefore inactivating any pro-apoptotic effect due to the NaVβ3 expression in basal conditions. This is one of the several ways the HPV genome secures the growing of the tumor. However, the intermediate steps that NaVβ3 must undergo to promote apoptosis is still unknown, as well as why NaVβ3 had no pro-apoptotic effect in CeCa cells when overexpressed. A possible explanation could be the existence of a positive (and necessary) feedback regulation between p53 (or one of its multiple target genes) and NaVβ3, leading to a lack of pro-apoptotic activity of NaVβ3 in the absence of functional p53.
On the contrary, our migration experiments show a clear participation (negative regulation) of NaVβ1 on the ability of CeCa cells to perform transwell migration across of an 8 µm-pore membrane. When we induced the overexpression of this subunit, the migration rate decreased in the three CeCa cell lines, being HeLa cells the most affected (almost 50% inhibition in cell migration). Correspondingly, the inhibition of NaVβ1 expression by siRNAs led to significant increments in cell migration of HeLa, SiHa and CaSki cells (Fig. 4a). Thus, NaVβ1 could be acting as a cell adhesion molecule, promoting the cell–cell adhesion and therefore making it difficult to move across the membrane. It is also worth of notice that, among the three cell lines, HeLa have the lowest NaVβ1 expression compared to the other two CeCa cell lines (Fig. 1), and they also have the greatest migration capacity. In the same sense, CaSki cells are the ones with the bigger expression of NaVβ1, and the lowest migration behavior (Fig. 2b).
A study made on breast cancer cell lines revealed that the overexpression of NaVβ1 in highly metastatic MDA-MB-231 cells increased cell–cell adhesion and decreased in vitro migration, consistent with the proposed role of NaVβ1 as a cell adhesion molecule . However, the trans-homophilic interactions through NaVβ1 can also mediate the process outgrowth in vivo, as in neurons, generating an elongated morphology and therefore potentiating metastasis, as has been reported for the same breast cancer cells but using mouse models instead of transwell experiments . Something similar could be happening with CeCa cell lines. Thus, the persistent contribution of NaVβ1 in the three CeCa cell lines migration is indicative of its role as a negative regulator of the in vitro migration, regardless the HPV type present in the CeCa cell line. Although it remains to be investigated whether this behavior also occurs in vivo.
With respect to the contribution of NaVβs in invasiveness, our results show that downregulation of NaVβ4 expression resulted in a significant increase in the in vitro invasiveness of the three cell lines studied here (Fig. 4b), suggesting the role of this subunit as a metastasis suppressor protein. These observations agree with the in silico analysis of NaVβ4 expression performed with the Human Protein Atlas data: normal cervix biopsies show moderated NaVβ4 expression, while in CeCa biopsies it seems to be absent (see Additional file 5). Taken together, these results strongly suggest that NaVβ4 expression is downregulated in CeCa, which in turn increases the invasive capacity of the cells. In accordance with this interpretation, out of the three CeCa cell lines studied in the present work, the HeLa cell line showed the lowest expression of NaVβ4 and the highest index of cell invasion.
In breast cancer cell lines, it has been found that increasing NaVβ1 expression promotes the overexpression of NaV1.5, the main NaV α-subunit involved in the invasiveness of this type of cancer . However, according to our results, this is not the case for CeCa cell lines, as the overexpression of NaVβ1 or NaVβ4 did not promote the expression of conducting NaV1.6 channels in the plasma membrane, the principal NaV α-subunit involved in the invasive potential of these cells. Thus, we conclude that the pathways used by these proteins to induce cell migration and invasion, are independent of the activity of NaV1.6 or other NaV α-subunits in the plasma membrane of CeCa cell lines.
Other groups have also reported the downregulation of NaVβ4 expression in cancer vs normal tissue. This is the case of prostate cancer , papillary thyroid cancer , and breast cancer . The latter study, performed by the group of Dr. Sebastien Roger at the University of Tours, reported solid evidence about the possible mechanistic pathway through which NaVβ4 prevents metastasis in cancer cells. They showed an increased activity of the small GTPase RhoA (crucial for cytoskeleton remodeling during migration and invasion) when NaVβ4 was silenced in breast cancer MDA-MB-231 cells. This effect was independent of NaV1.5 channel activity, the main NaV α-subunit implicated in BCa invasiveness. A similar mechanism might be taking place in CeCa cells: a downregulation of the NaVβ4 subunit that leads to an increased activity of RhoA and therefore, a more robust invasive potential, which is independent of the NaV1.6 expression in the plasma membrane. However, we speculate that additional cellular mechanisms might be involved in this phenomenon, as the effects of NaVβ4 in cellular migration were not significant in our study, meaning that the sole change on cellular motility cannot explain the dramatic increase on cellular invasiveness when NaVβ4 is downregulated. A more recent study showed that NaVβ4 is downregulated in papillary thyroid cancer (PTC) compared with normal thyroid tissues. More interestingly, the authors also found that preserved NaVβ4 expression might independently predict favorable recurrence-free survival in classical PTC . Altogether, these data reinforce the notion of using NaVβ4 as a biomarker for cancer metastasis and a potential new therapeutic target for the treatment of cervical cancer.
In conclusion, our results show that even though cervical cancer cell lines do not generate NaV currents, they do express NaVβ subunits, which could play a role in proliferation (NaVβ1 only in SiHa cells), migration (NaVβ1 in the three cell lines), and invasion (NaVβ4 in the three cell lines studied). The results support the recent knowledge about NaVβs as multifunctional proteins involved in cell processes like ion channel regulation, cell adhesion and motility, and even in metastatic cell behaviors. Finally, the evidence presented here contributed to propose NaVβ1 and NaVβ4 as prognostic markers for cervical cancer tumor progression.
cell adhesion molecules
Dulbecco’s Modified Eagle’s Medium
fetal bovine serum
- NaV :
voltage-gated sodium channels
voltage-gated sodium channels β-subunits
voltage-gated sodium channel coded by the SCN8A α-subunit
papillary thyroid cancer
real-time polymerase chain reaction
reverse transcription polymerase chain reaction
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.
Rana MM, Huhtala H, Apter D, Eriksson T, Luostarinen T, Natunen K, et al. Understanding long-term protection of human papillomavirus vaccination against cervical carcinoma: cancer registry-based follow-up. Int J Cancer J Int du Cancer. 2013;132(12):2833–8.
Catterall WA. Voltage-gated sodium channels at 60: structure, function and pathophysiology. J Physiol. 2012;590(Pt 11):2577–89.
Fraser SP, Diss JK, Chioni AM, Mycielska ME, Pan H, Yamaci RF, et al. Voltage-gated sodium channel expression and potentiation of human breast cancer metastasis. Clin Cancer Res. 2005;11(15):5381–9.
Brisson L, Gillet L, Calaghan S, Besson P, Le Guennec JY, Roger S, et al. Na(V)1.5 enhances breast cancer cell invasiveness by increasing NHE1-dependent H(+) efflux in caveolae. Oncogene. 2011;30(17):2070–6.
Laniado ME, Lalani EN, Fraser SP, Grimes JA, Bhangal G, Djamgoz MB, et al. Expression and functional analysis of voltage-activated Na+ channels in human prostate cancer cell lines and their contribution to invasion in vitro. Am J Pathol. 1997;150(4):1213–21.
Diss JK, Archer SN, Hirano J, Fraser SP, Djamgoz MB. Expression profiles of voltage-gated Na(+) channel alpha-subunit genes in rat and human prostate cancer cell lines. Prostate. 2001;48(3):165–78.
House CD, Vaske CJ, Schwartz AM, Obias V, Frank B, Luu T, et al. Voltage-gated Na+ channel SCN5A is a key regulator of a gene transcriptional network that controls colon cancer invasion. Cancer Res. 2010;70(17):6957–67.
Gao R, Shen Y, Cai J, Lei M, Wang Z. Expression of voltage-gated sodium channel alpha subunit in human ovarian cancer. Oncol Rep. 2010;23(5):1293–9.
Roger S, Rollin J, Barascu A, Besson P, Raynal PI, Iochmann S, et al. Voltage-gated sodium channels potentiate the invasive capacities of human non-small-cell lung cancer cell lines. Int J Biochemis Cell Biol. 2007;39(4):774–86.
Xia J, Huang N, Huang H, Sun L, Dong S, Su J, et al. Voltage-gated sodium channel Nav 1.7 promotes gastric cancer progression through MACC1-mediated upregulation of NHE1. Int J Cancer. 2016;139(11):2553–69.
Carrithers MD, Chatterjee G, Carrithers LM, Offoha R, Iheagwara U, Rahner C, et al. Regulation of podosome formation in macrophages by a splice variant of the sodium channel SCN8A. J Biol Chem. 2009;284(12):8114–26.
Xing D, Wang J, Ou S, Wang Y, Qiu B, Ding D, et al. Expression of neonatal Nav1.5 in human brain astrocytoma and its effect on proliferation, invasion and apoptosis of astrocytoma cells. Oncol Rep. 2014;31(6):2692–700.
Hernandez-Plata E, Ortiz CS, Marquina-Castillo B, Medina-Martinez I, Alfaro A, Berumen J, et al. Overexpression of NaV 1.6 channels is associated with the invasion capacity of human cervical cancer. Int J Cancer J Int du Cancer. 2012;130(9):2013–23.
Isom LL, Ragsdale DS, De Jongh KS, Westenbroek RE, Reber BF, Scheuer T, et al. Structure and function of the beta 2 subunit of brain sodium channels, a transmembrane glycoprotein with a CAM motif. Cell. 1995;83(3):433–42.
Malhotra JD, Kazen-Gillespie K, Hortsch M, Isom LL. Sodium channel beta subunits mediate homophilic cell adhesion and recruit ankyrin to points of cell-cell contact. J Biol Chem. 2000;275(15):11383–8.
Ratcliffe CF, Westenbroek RE, Curtis R, Catterall WA. Sodium channel beta1 and beta3 subunits associate with neurofascin through their extracellular immunoglobulin-like domain. J Cell Biol. 2001;154(2):427–34.
Malhotra JD, Koopmann MC, Kazen-Gillespie KA, Fettman N, Hortsch M, Isom LL. Structural requirements for interaction of sodium channel beta 1 subunits with ankyrin. J Biol Chem. 2002;277(29):26681–8.
Yu FH, Westenbroek RE, Silos-Santiago I, McCormick KA, Lawson D, Ge P, et al. Sodium channel beta4, a new disulfide-linked auxiliary subunit with similarity to beta2. J Neurosci. 2003;23(20):7577–85.
McEwen DP, Isom LL. Heterophilic interactions of sodium channel beta1 subunits with axonal and glial cell adhesion molecules. J Biol Chem. 2004;279(50):52744–52.
Brackenbury WJ, Isom LL. Na channel beta subunits: overachievers of the ion channel family. Front Pharmacol. 2011;2:53.
Yereddi NR, Cusdin FS, Namadurai S, Packman LC, Monie TP, Slavny P, et al. The immunoglobulin domain of the sodium channel beta3 subunit contains a surface-localized disulfide bond that is required for homophilic binding. FASEB J. 2013;27(2):568–80.
Diss JK, Fraser SP, Walker MM, Patel A, Latchman DS, Djamgoz MB. Beta-subunits of voltage-gated sodium channels in human prostate cancer: quantitative in vitro and in vivo analyses of mRNA expression. Prostate Cancer Prostatic Dis. 2008;11(4):325–33.
Jansson KH, Lynch JE, Lepori-Bui N, Czymmek KJ, Duncan RL, Sikes RA. Overexpression of the VSSC-associated CAM, beta-2, enhances LNCaP cell metastasis associated behavior. Prostate. 2012;72(10):1080–92.
Jansson KH, Castillo DG, Morris JW, Boggs ME, Czymmek KJ, Adams EL, et al. Identification of beta-2 as a key cell adhesion molecule in PCa cell neurotropic behavior: a novel ex vivo and biophysical approach. PLoS ONE. 2014;9(6):e98408.
Adachi K, Toyota M, Sasaki Y, Yamashita T, Ishida S, Ohe-Toyota M, et al. Identification of SCN3B as a novel p53-inducible proapoptotic gene. Oncogene. 2004;23(47):7791–8.
Chioni AM, Brackenbury WJ, Calhoun JD, Isom LL, Djamgoz MB. A novel adhesion molecule in human breast cancer cells: voltage-gated Na + channel beta1 subunit. Int J Biochemis Cell Biol. 2009;41(5):1216–27.
Nelson M, Millican-Slater R, Forrest LC, Brackenbury WJ. The sodium channel beta1 subunit mediates outgrowth of neurite-like processes on breast cancer cells and promotes tumour growth and metastasis. Int J Cancer J Int du Cancer. 2014;135(10):2338–51.
Bon E, Driffort V, Gradek F, Martinez-Caceres C, Anchelin M, Pelegrin P, et al. SCN4B acts as a metastasis-suppressor gene preventing hyperactivation of cell migration in breast cancer. Nat Commun. 2016;7:13648.
Diaz D, Delgadillo DM, Hernandez-Gallegos E, Ramirez-Dominguez ME, Hinojosa LM, Ortiz CS, et al. Functional expression of voltage-gated sodium channels in primary cultures of human cervical cancer. J Cell Physiol. 2007;210(2):469–78.
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(4):402–8.
Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419.
Uhlen M, Oksvold P, Fagerberg L, Lundberg E, Jonasson K, Forsberg M, et al. Towards a knowledge-based human protein atlas. Nat Biotechnol. 2010;28(12):1248–50.
Lopez-Charcas O, Espinosa AM, Alfaro A, Herrera-Carrillo Z, Ramirez-Cordero BE, Cortes-Reynosa P, et al. The invasiveness of human cervical cancer associated to the function of NaV1.6 channels is mediated by MMP-2 activity. Sci Rep. 2018;8(1):12995.
Isom LL, De Jongh KS, Patton DE, Reber BF, Offord J, Charbonneau H, et al. Primary structure and functional expression of the beta 1 subunit of the rat brain sodium channel. Science. 1992;256(5058):839–42.
Chen C, Bharucha V, Chen Y, Westenbroek RE, Brown A, Malhotra JD, et al. Reduced sodium channel density, altered voltage dependence of inactivation, and increased susceptibility to seizures in mice lacking sodium channel beta 2-subunits. Proc Natl Acad Sci USA. 2002;99(26):17072–7.
Crosbie EJ, Einstein MH, Franceschi S, Kitchener HC. Human papillomavirus and cervical cancer. Lancet. 2013;382(9895):889–99.
Gey GO, Coffman WD, Kubicek MT. Tissue culture studies of the proliferative capacity of cervical carcinoma and normal epithelium. Can Res. 1952;12(4):264–5.
Friedl F, Kimura I, Osato T, Ito Y. Studies on a new human cell line (SiHa) derived from carcinoma of uterus. I. Its establishment and morphology. Proc Soc Exp Biol Med. 1970;135(2):543–5.
Pattillo RA, Hussa RO, Story MT, Ruckert AC, Shalaby MR, Mattingly RF. Tumor antigen and human chorionic gonadotropin in CaSki cells: a new epidermoid cervical cancer cell line. Science. 1977;196(4297):1456–8.
Kim DY, Carey BW, Wang HB, Ingano LAM, Binshtok AM, Wertz MH, et al. BACE1 regulates voltage-gated sodium channels and neuronal activity. Nat Cell Biol. 2007;9(7):755–64.
Kim DY, Ingano LAM, Carey BW, Pettingell WH, Kovacs DM. Presenilin/gamma-secretase-mediated cleavage of the voltage-gated sodium channel beta 2-subunit regulates cell adhesion and migration. J Biol Chem. 2005;280(24):23251–61.
Ju YK, Saint DA, Gage PW. Hypoxia increases persistent sodium current in rat ventricular myocytes. J Physiol-London. 1996;497(2):337–47.
Fearon IM, Brown ST. Acute and chronic hypoxic regulation of recombinant hNa(v)1.5 alpha subunits. Biochem Bioph Res Co. 2004;324(4):1289–95.
Scheffner M, Huibregtse JM, Vierstra RD, Howley PM. The HPV-16 E6 and E6-AP complex functions as a ubiquitin-protein ligase in the ubiquitination of p53. Cell. 1993;75(3):495–505.
Gong Y, Yang J, Wu W, Liu F, Su A, Li Z, et al. Preserved SCN4B expression is an independent indicator of favorable recurrence-free survival in classical papillary thyroid cancer. PLoS ONE. 2018;13(5):e0197007.
ALSS conducted all the experiments, analyzed the results, and wrote the original version of the manuscript. JCG conceived the idea for the project, analyzed results and wrote the manuscript with ALSS. Both authors read and approved the final manuscript.
We are very grateful with Dr. Lori L. Isom for the donation of the four NaVβs vectors; and Dr. Laura Ongay, from Unidad de Biología Molecular at Instituto de Fisiología Celular, UNAM, for expert technical support.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its additional information files].
Consent for publication
Ethics approval and consent to participate
This work was supported by grants from CONACYT-México (167790-B) and PAPIIT-DGAPA-UNAM (IN206917) to JCG. Ana Laura Sanchez-Sandoval is a doctoral student from Programa de Doctorado en Ciencias Bioquímicas, Universidad Nacional Autónoma de México (UNAM) and received fellowship from CONACYT. The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.