Assessment of gene expression of intracellular calcium channels, pumps and exchangers with epidermal growth factor-induced epithelial-mesenchymal transition in a breast cancer cell line

Background Epithelial-mesenchymal transition (EMT) is a process implicated in cancer metastasis that involves the conversion of epithelial cells to a more mesenchymal and invasive cell phenotype. In breast cancer cells EMT is associated with altered store-operated calcium influx and changes in calcium signalling mediated by activation of cell surface purinergic receptors. In this study, we investigated whether MDA-MB-468 breast cancer cells induced to undergo EMT exhibit changes in mRNA levels of calcium channels, pumps and exchangers located on intracellular calcium storing organelles, including the Golgi, mitochondria and endoplasmic reticulum (ER). Methods Epidermal growth factor (EGF) was used to induce EMT in MDA-MB-468 breast cancer cells. Serum-deprived cells were treated with EGF (50 ng/mL) for 12 h and gene expression was assessed using quantitative RT-PCR. Results and conclusions These data reveal no significant alterations in mRNA levels of the Golgi calcium pump secretory pathway calcium ATPases (SPCA1 and SPCA2), or the mitochondrial calcium uniporter (MCU) or Na+/Ca2+ exchanger (NCLX). However, EGF-induced EMT was associated with significant alterations in mRNA levels of specific ER calcium channels and pumps, including (sarco)-endoplasmic reticulum calcium ATPases (SERCAs), and inositol 1,4,5-trisphosphate receptor (IP3R) and ryanodine receptor (RYR) calcium channel isoforms. The most prominent change in gene expression between the epithelial and mesenchymal-like states was RYR2, which was enriched 45-fold in EGF-treated MDA-MB-468 cells. These findings indicate that EGF-induced EMT in breast cancer cells may be associated with major alterations in ER calcium homeostasis.


Background
EMT facilitates cancer cell invasion and metastasis formation, and has also been linked to the acquisition of a stem cell-like phenotype, anchorage-independent growth and chemoresistance in cancer cell lines and clinical samples [1][2][3][4]. The phenotypic changes associated with EGF-induced EMT are well characterised in the human breast cancer cell line MDA-MB-468, and include changes in cell morphology; increased expression of the transcription factor Twist and the intermediate filament protein vimentin; and reduced E-cadherin expression following chronic EGF treatment [5][6][7][8]. In addition, EGF-induced EMT in MDA-MB-468 breast cancer cells has been linked to altered plasma membrane calcium influx [6,7,9].
Increases in cytosolic calcium may occur via calcium entry through calcium permeable plasma membrane channels, including transient receptor potential (TRP), voltage-gated and store-operated calcium channels, or due to the release of calcium from internal stores [10][11][12]. Calcium signals arising via influx or storerelease have distinct spatiotemporal profiles, modes of activation, and cellular responses [13]. These two processes often occur in parallel, for example calcium store release triggering store-operated calcium influx [14], and calcium influx via specific plasma membrane calcium channels potentially causing calcium release from the ER via calcium-induced calcium release (CICR) [15]. Mitochondrial calcium buffering may also alter the nature of cytosolic calcium responses [16][17][18][19]. Thus changes in the expression of calcium channels, pumps and exchangers located on intracellular calcium storing organelles may alter the nature of the cytosolic calcium response to a range of stimuli. However, while altered store-operated [6,9] and non-stimulated [6] plasmalemmal calcium influx have been linked to the acquisition of a more mesenchymal phenotype in breast cancer cells, changes in calcium channels, pumps and exchangers located on intracellular calcium storing organelles have not been assessed during this process.
A phenotypic conversion and cellular heterogeneity comparable to EMT, which is associated with a remodelling of intracellular calcium homeostasis, is observed in vascular smooth muscle cells [20]. In response to cues from the microenvironment, vascular smooth muscle cells can switch from a contractile to a synthetic (proliferative) phenotype [20]. Vascular smooth muscle cell switching is associated with changes in the expression of ER calcium channels and pumps, including reduced expression of RYR3 and SERCA2a [20,21]. Vascular smooth muscle cells are therefore one example of how intracellular calcium channels and pumps may be remodelled to meet the cellular requirements of specific phenotypes.
In this study we assessed whether alterations in mRNA levels of Golgi, mitochondrial and ER calcium channels, pumps and exchangers were associated with EGF-induced EMT in MDA-MB-468 breast cancer cells.

Results and discussion
EGF-induced EMT in MDA-MB-468 cells MDA-MB-468 breast cancer cells treated with EGF (50 ng/mL) for 3, 6 or 12 h showed a significant increase in mRNA levels of Twist ( Figure 1A), a nuclear transcription factor and marker of EMT [22]. This early EMT event was followed by an increase in vimentin mRNA at 12 h post-EGF treatment ( Figure 1B). These findings were consistent with previous studies using this model of EGF-induced EMT [5][6][7][8]. We assessed gene expression of calcium channels, pumps and exchangers located on the Golgi, mitochondria and ER with EMT after treatment with EGF for 12 h.

Golgi: SPCA1 and SPCA2
SPCAs are regulators of Golgi luminal calcium levels [23]. Recently, we showed that SPCA1 is enriched in basal-like breast cancers, a class of aggressive breast cancers with EMT-like features [24]. Given the potential for SPCAs to regulate processes important for carcinogenesis [24,25] and the observed repositioning of the Golgi, in some cell types, to the leading edge during cell migration [26], we hypothesised that altered Golgi calcium homeostasis may be a feature of EGF-induced EMT. However, no significant alterations in the levels of

Mitochondria: MCU and NCLX
Excess mitochondrial calcium accumulation is a trigger for cell death [27]. Downregulation of the MCU in HeLa cervical cancer cells by miR-25 protects cells from calcium-dependent apoptosis, and reduced expression of MCU is a feature of poorly differentiated colonic adenocarcinomas [28]. In contrast, MCU expression appears to be enriched in some breast cancer subtypes associated with a poor clinical prognosis [29]. Although several studies have reported that EMT or expression of EMTassociated transcription factors is associated with apoptosis resistance [30][31][32], the mechanism underpinning this process has not been characterised. We evaluated changes in gene expression of two key mitochondrial calcium accumulation and release proteins-the MCU [33,34] and NCLX [35]-in breast cancer cells induced to undergo EMT with EGF. Our results show no significant changes in the expression of the MCU or NCLX with EGF-induced EMT in MDA-MB-468 cells ( Figure 2C & D). These findings indicate that if EGFinduced EMT in breast cancer cells is associated with apoptosis-resistance, it is not due to altered gene expression of the MCU or NCLX, at least under these study conditions.

Endoplasmic reticulum: SERCAs, IP 3 Rs and RYRs
The ER is a major source of releasable calcium in epithelial cells [36]. Luminal ER calcium levels are maintained by SERCA calcium pumps, which are encoded by three genes (ATP2A/SERCA1-3) [37]. Activation of phospholipase C-coupled receptors (e.g., some G-protein coupled receptors) leads to the production of IP 3 and subsequent activation of the ER calcium-permeable channels IP 3 R1, IP 3 R2 and IP 3 R3 [38]. The nature of the cytosolic calcium signal evoked by IP 3 is highly dependent on the expression profile of IP 3 Rs [39]. For example, in B cells, activation of IP 3 R2 produces regular and sustained calcium oscillations, while IP 3 R3 generates a calcium response that is characterised by a transient cytosolic calcium peak and few oscillatory spikes [39]. IP 3 R subtypes also differ in their sensitivity to IP 3 [39]. To assess possible alterations in ER calcium handling, we assessed both SERCA and IP 3 R isoform mRNA levels in MDA-MB-468 breast cancer cells induced to undergo EMT. No significant change in SERCA1 gene expression was observed ( Figure 3A). However, an increase in SERCA2 and a down-regulation of SERCA3 was associated with EGF-induced EMT ( Figure 3B & C). EGF-induced EMT was also associated with a 2.5-fold increase in IP 3 R1 and a 2.2-fold increase in IP 3 R3 mRNA levels ( Figure 3 D-F). Changes in gene expression of SERCA2, SERCA3, IP 3 R1 and IP 3 R3 suggest that alterations in ER calcium homeostasis may be a characterising feature of EGF-induced EMT in MDA-MB-468 breast cancer cells.
To determine whether alterations in ER calcium signalling with EMT could also be mediated via alterations in the ryanodine receptor family of (sarco)-endoplasmic reticulum calcium release channels, mRNA levels of ryanodine receptor isoforms were assessed during EGFinduced EMT. There are three ryanodine receptor isoforms in humans; RYR1 is predominantly expressed skeletal muscle; RYR2 is the abundant isoform in cardiac muscle; and RYR3 expression occurs in various tissues including the brain [40]. RYR1 and RYR3 were not detected in MDA-MB-468 epithelial breast cancer cells and expression of these isoforms was not induced to detectable levels with EGF-induced EMT (data not shown). Levels of mRNA for the RYR2 subtype were low in MDA-MB-468 breast cancer cells prior to EMT induction with EGF. However, EGF-treatment induced a marked (46-fold) increase in levels of RYR2 mRNA ( Figure 3G). We were unable to detect RYR2 expression in MDA-MB-468 breast cancer cells with immunoblotting using a commercially available antibody (Abcam C3-33, data not shown); however, corresponding to the increase in RYR2 gene expression, MDA-MB-468 breast cancer cells with a more mesenchymal phenotype exhibited a modest but significantly greater response to the ryanodine receptor activator caffeine, as assessed through caffeine-mediated increases in [Ca 2+ ] CYT ( Figure 4A & B).

Conclusions
Alterations in the expression of ER calcium channels and pumps occur in cancer cell lines and clinical cancer samples [12,[41][42][43][44][45]. For example, expression of SERCA3 is significantly reduced in colon carcinomas and colon cancer cell lines compared to normal colonic epithelia [41], and SERCA3 expression is elevated with chemically-induced differentiation of colon, gastric and lung cancer cell lines [41,42]. Therefore, the reduced SERCA3 mRNA levels observed in this study with EGF-induced EMT may reflect a loss of cellular differentiation congruent with the acquisition of a more stem cell-like phenotype. Future studies are needed to assess whether these changes are related to EMTassociated stemness and to assess the various SERCA transcript variants with EGF-induced EMT. Ryanodine receptor expression, using an antibody that detects all ryanodine receptor isoforms, also correlates with breast cancer tumour grade [43]. Furthermore, the induction of RYR2 expression by factors/conditions known to induce EMT in a variety of cell types (including EGF, TGFβ [46] and hypoxia [47]), suggests that RYR2 may be a key calcium channel for the mesenchymal cell phenotype. The association between RYR2 and EMT requires further investigation. We have recently shown that the calcium signal is critically required for EMT induction in breast cancer cells [48] and that, in addition to regulating the conversion to a more mesenchymal cell phenotype, EMT is associated with a remodelling of non-stimulated, agonist-stimulated and store-operated calcium influx [6,7]. This study assessed changes in mRNA levels of intracellular calcium channels, pumps and exchangers with EMT. Our results demonstrate that similar to vascular smooth muscle cell phenotypic switching [20], EMT in breast cancer cells is associated with changes in gene expression of specific ER calcium channels and pumps, the consequence of which may be altered ER calcium storage and signalling. This is in agreement with our previously published data showing altered ER calcium release kinetics with EGF-induced EMT in this cell line [6]. Changes in ER calcium handling may lead to altered responses to phospholipase-C-coupled receptor agonists, and may be important for the invasive properties of mesenchymal-like cancer cells.

Statistical analysis
Statistical analysis was performed using GraphPad Prism (v6.0 for Windows). Statistical tests used in this study are outlined in each figure legend.