Introduction

Cadmium (Cd2+) is a metal pollutant prone to bioaccumulation, mainly in the liver and in the kidney (Chen et al. 2019; Rani et al. 2014). The major sources of Cd2+ exposure for humans are the cigarette smoke and contaminated drinking water and/or food. Cd2+ is considered as a class I human lung carcinogen since 2004 (Waalkes 2003). Moreover, Cd2+ accumulation in tissues is suspected to be linked to breast cancer development (Larsson et al. 2015). In vivo studies showed that a single intraperitoneal dose of Cd2+ (~ 27 nmol/kg) mimics estrogens leading to both endometrial and mammary gland hyperplasia and hypertrophy in ovariectomized female rats (Johnson et al. 2003). Moreover, Cd2+ stimulates MCF7 breast cancer cell proliferation through estrogen receptor (ER) activation (Martinez-Campa et al. 2006). In several breast cancer cell line models, Cd2+ potentiates the interaction between ER and c-Jun leading to the stabilization of this transcription factor complex and its recruitment to the promoters of cyclin D1 and c-myc (Siewit et al. 2010). Thus, Cd2+ has been defined as a metalloestrogen, an endocrine disrupter able to promote the development of estrogen-dependent cancers (Darbre 2006). However, chronic Cd2+ exposure can also induce non-cancer epithelial cell transformation into a highly invasive phenotype independently of ER expression (Benbrahim-Tallaa et al. 2009; Qu et al. 2012), suggesting that Cd2+ acts as a carcinogen through an alternative mechanism that do not involve estrogen receptors. Cd2+ can diffuse across the plasma membrane through several transporters including ion channels (Thevenod 2010). Among them, TRPM7 is a non-selective channel, mainly permeant to Mg2+/Ca2+ ions and also permeant to other divalent cations including Cd2+ (Monteilh-Zoller et al. 2003). It has been recently shown that TRPM7 is the main gatekeeper for Zn2+, Mg2+ and Ca2+ intestinal absorption (Mittermeier et al. 2019). TRPM7 is also involved in Cd2+ uptake, leading to cytotoxic accumulation in human osteoblasts. This latter uptake can possibly explain the occurrence of bone resorption and fracture following Cd2+ poisoning (Martineau et al. 2010). The low selectivity of TRPM7 towards divalent cations may also provide a favorite route for Cd2+ in epithelial cell leading to cell transformation into a neoplastic phenotype. Numerous studies have shown that TRPM7 is involved in cancer progression (Gautier et al. 2016). For instance, TRPM7 is overexpressed in breast cancer and this channel regulates MCF7 cell proliferation (Dhennin-Duthille et al. 2011; Guilbert et al. 2009). Moreover, TRPM7 also regulates breast cancer cell migration and invasion in MDA-MB231 and MDA-MB-435S cell lines (Guilbert et al. 2013; Meng et al. 2013), as well as metastasis formation and dissemination in vivo in a mouse xenograft model (Middelbeek et al. 2012). The present study aims to assess the role of TRPM7 in the acquisition of an aggressive phenotype by non-cancer epithelial cells chronically exposed to cadmium.

Material and methods

Cell culture and cadmium exposure

We used two human non-cancer epithelial cell lines: one from breast (MCF10A), and one from pancreas (hTERT-HPNE). MCF10A were purchased from ATCC® (CRL-10317™) and were exposed with 2.5 µM Cd2+ for 40 weeks according to previously published data (Benbrahim-Tallaa et al. 2009). hTERT-HPNE were purchased by ATCC® (CRL-4023™) and exposed with 1 µM Cd2+ for 30 weeks (Qu et al. 2012). MCF10A were cultured in Dulbecco’s Modified Eagle’s Medium/Nutrient Mixture F-12 (DMEM/F-12) supplemented with 10% Fetal Calf Serum (FCS), hydroxycortisone (50 ng/mL), insulin (0.01 mg/mL), human recombinant EGF (20 ng/mL), and cholera toxin (100 ng/mL). hTERT-HPNE were cultured in base medium: 75% DMEM without glucose (Sigma) and 25% Medium M3 Base (Incell Corp) with 5% FCS. Complete growth medium was supplemented with the following components: human recombinant EGF (10 ng/mL), D-glucose (1 g/L), puromycin (750 ng/mL), l-glutamine (0.3 g/L), and sodium bicarbonate (1.5 g/L). Cells were trypsinized once a week using trypsin–EDTA (Sigma-Aldrich, Inc.) and incubated at 37 °C in a humidified atmosphere containing 5% CO2.

Small interference RNA technique

siRNAs were transfected in non-cancer cells by nucleofection using a Nucleofector™ II device (Lonza). MCF10A and hTERT-HPNE cells (106 cells) were transfected with 2 μg of siRNA according to the optimized protocol recommended by Lonza. Non-cancer cells were nucleofected with siRNA targeting TRPM7 (siTRPM7) or a scrambled siRNA (siControl) as previously described (Rybarczyk et al. 2017). All the experiments were performed 48 h following nucleofection.

TRPM7 overexpression

Zebrafish (ZF) TRPM7 was overexpressed in MCF10A cells by nucleofection (106 cells, program X-005, kit L). MCF10A cells were transfected with 2 µg pcDNA5-TO-ZF-WT-TRPM7 (WT-TRPM7) or 2 µg pcDNA5-TO as previously published (Guilbert et al. 2013). All experiments were performed 72 h after transfection.

Western-blot

Cells were lysed 30 min on ice in radioimmunoprecipitation assay buffer (1% Triton X-100, 1% Na deoxycholate, 150 mM NaCl, 10 mM PO4Na2/K, pH 7.2) supplemented with Sigma P8340 inhibitors cocktail, 2 mM EDTA, and 5 mM orthovanadate. After centrifugation at 13,000 rpm, the proteins in the supernatant were quantified using the BCA method (BioRad). Equal amounts of each protein sample (50 μg) were separated by electrophoresis on sodium dodecyl sulfate (SDS) polyacrylamide gel electrophoresis and blotted onto nitrocellulose membrane (Amersham). Blots were incubated with antibodies raised against TRPM7 (1/1000, Abcam), N-Cadherin (1/1000, Abcam), E-Cadherin (1/1000, Abcam), Vimentin (1/250, ThermoFisher) and Tubulin (1/7000, Sigma-Aldrich). Blots were developed with the enhanced chemiluminescence system using specific peroxidase-conjugated anti-IgG secondary antibodies.

Electrophysiological recordings

To study ion channel currents, we used the conventional patch-clamp technique (Hamill et al. 1981). TRPM7 currents recordings were performed using the whole-cell configuration as previously described (Rybarczyk et al. 2017). Briefly, membrane potential was held at − 40 mV, and currents were elicited by a ramp depolarization from − 100 mV to + 100 mV for 350 ms. Interval between each ramp depolarization was 10 s. After the dialysis of intracellular media by the free Mg intrapipette solution, the Magnesium Inhibited Cation (MIC) current was recorded (Prakriya and Lewis 2002). To measure the MIC current, the difference between the steady-state current activated by the depletion of [Mg2+]i and the basal current recorded few minutes after patch rupture was calculated. Currents were expressed as current densities (in pA.pF−1) by dividing the current intensity (in pA) by the cell capacitance (in pF). All experiments were performed at room temperature.

Manganese quenching and intracellular calcium and magnesium imaging

Non-cancer cells were plated on glass coverslips in 35-mm diameter dishes at a density of 3 × 104 and 5 × 104 cells for hTERT-HPNE and MCF10A, respectively. After 48 h, cells were loaded in cell growth medium at + 37 °C for different times (35 min for hTERT-HPNE and 40 min for MCF-10A) with 3 μM Fura-2 (Sigma-Aldrich) or 3 μM Mag-Fura-2 (Invitrogen) and subsequently washed three times with the extracellular solution (in mM: NaCl 145; KCl 5; HEPES 10; Glucose 5; MgCl2 2; CaCl2 as indicated; MnCl2 as indicated; CdCl2 as indicated; pH adjusted to 7.4 with NaOH). The coverslip was then transferred onto a perfusion chamber on a Zeiss microscope equipped with fluorescence. Manganese (Mn2+) is a cation able to fix Fura-2 fluoroprobe and to quench its fluorescence (Fasolato et al. 1993). To estimate the divalent cation influx, we used the Mn2+ quenching protocol as previously described (Rybarczyk et al. 2017). Briefly, cells were perfused for 1 min with the extracellular solution. Then, Ca2+ was replaced by Mn2+ (2 mM) for 2 min. Cells were excited at 360 nm and fluorescence emission was monitored at 510 nm. After Mn2+ perfusion, the decrease of Fura-2 fluorescence followed a linear decay, and the slope is correlated with the rate of the Mn2+ influx. The calculated slope is obtained by subtracting the slope of basal decreasing Fura-2 fluorescence obtained in basal conditions and after Mn2+ application. Similar experiments were made using Cd2+ (2 mM) instead of Mn2+ to quench Fura-2 fluorescence. Contrarily to Mn2+, Cd2+ entry increases Fura-2 fluorescence. Moreover, MCF10A and hTERT-HPNE were incubated with Fura-2 (as previously indicated) to measure the Ca2+ basal fluorescence ratio. During the perfusion with the extracellular solution (containing 2 mM of Ca2+), cells were excited alternatively at 350 nm and 380 nm. The ratio of Fura-2 fluorescence intensities measured with excitation at 350 and 380 nm (F350/F380) was used as a [Ca2+]i related signal. To evaluate Mg2+ basal ratio, cells were perfused with an extracellular solution containing 1 mM Mg2+, and Mag-Fura-2 was used as a fluoroprobe. Mag-Fura-2 fluorescence intensities (measured following excitation at 330 and 370 nm, F330/F370) were used as a [Mg2+]i related signal.

Cell morphology assessment

After the Cd2+ exposure period, cells were photographed under an inverted microscope (Zeiss) at intermediate confluence. Cells clustering at × 200 magnification was evaluated. Circularity index of cells at × 400 magnification was calculated by Image J software (https://imagej.nih.gov/ij/).

Wound healing assay

The effect of Cd2+ exposure on cellular motility was assessed using the wound healing assay. 3 × 105 cells were cultured in a 35-mm diameter dishes. When the cells have formed a confluent cell monolayer, a linear empty space was generated using sterile plastic pipette tip. The wound closure was estimated by measuring its area for each time (T0, T12h, T16h, T20h, T24h) using the Image J software (https://imagej.nih.gov/ij/).

Cell migration and invasion assays

For invasion tests, cell culture inserts (BD Falcon Cell Culture Inserts, BD Biosciences) were coated with Matrigel (Sigma-Aldrich, Inc.). The upper compartment was seeded with 2 × 104 cells and 4 × 104 of 24-h transfected hTERT-HPNE and MCF10A cells, respectively, in growth medium supplemented with FCS. The lower compartments were also filled with growth medium supplemented with FCS. Thus, migration and invasion assays were performed without addition of chemoattractant. After 24 h of incubation at 37 °C, inserts were washed by phosphate-buffered saline, fixed 10 min by methanol and stained 5 min with hematoxylin. The remaining cells on upper side were removed from the membrane by scrubbing. Cells in 20 contiguous areas at × 400 magnification for each insert were counted under an inverted microscope (Zeiss). For each experiment, the number of migrant or invasive cells per area for each condition was normalized by the mean of non-exposed siControl migrant or invasive cells.

Cell proliferation assays

Cell proliferation was tested by MTT assay. Cells were seeded in 6-well plates in the same density as migration and invasion assays to prove that the variations observed are not due to a seeding difference. Then, the yellow tetrazolium salts of 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide (MTT, Sigma-Aldrich, Inc.) solubilized in culture medium without FCS (0.5 mg/mL) was added to each well and incubated 50 min at + 37 °C in the dark. To dissolve purple formazan crystals, the culture medium was replaced with dimethyl sulfoxide (DMSO, Sigma-Aldrich, Inc.). Absorbance of each well was quantified at 550 nm using an Infinite® 200 Pro reader (Tecan Trading AG).

Statistical analysis

Data are presented as Mean ± S.E.M. and the number of experiments or the number of studied cells is represented by n. Experiments were repeated at least in three different cell passages. Data analysis and figure conception were made using Microcal™ Origin® (Microcal Software, Inc., USA) and Clampfit (Molecular Devices, Inc., USA) software. Statistical analyses were made using Student’s t tests or Mann–Whitney rank sum test depending on sample normality determined by paired Wilcoxon signed-rank test using Sigma-Stat 3.0 (Systat Software, Inc.). To compare the combinate effects of Cd2+ exposure and TRPM7 silencing, a two-way analysis of variance was used followed by Holm-Sidak method using SigmaStat 3.0 (Systat Software, Inc.).

Results

Effect of chronic cadmium exposure on TRPM7 channel expression in non-cancer epithelial cells

Western Blot analyses were performed to evaluate how TRPM7 expression was altered in MCF10A cells exposed to cadmium (Cd2+) when compared to the control cells. Results (displayed in Fig. 1a) show that TRPM7 relative expression increased by almost twofold in Cd2+ cells (n = 5) compared to control cells (n = 6; P < 0.005). Furthermore, TRPM7 knockdown (n = 6) decreased TRPM7 relative expression by almost 75% compared to siControl (n = 6; P < 0.05) in control cells, and by almost 80% in Cd2+ cells (n = 6) compared to siControl (n = 5; P < 0.001). These data show that Cd2+ exposure increased TRPM7 expression in MCF10A cells by twofold. As TRPM7 is also overexpressed in human pancreatic ductal adenocarcinoma (PDAC) (Rybarczyk et al. 2012; Yee et al. 2015), we further assess its expression in the non-cancer pancreatic epithelial cell line, hTERT-HPNE. We confirmed that Cd2+ exposure increased TRPM7 expression in hTERT-HPNE in a similar manner than for MCF10A (Fig. 1b).

Fig. 1
figure 1

TRPM7 expression and electrophysiological activity in human epithelial cells chronically exposed to cadmium. a, b Upper panels. Lysates from MCF10A (a) and hTERT-HPNE (b) exposed or not to cadmium (Cd2+) and transfected with a scrambled siRNA (siControl) or targeting TRPM7 (siTRPM7) were immunoblotted with the anti-TRPM7 and the anti-tubulin antibodies. a, b Lower panels. Quantification of immunoblotting normalized to control cells transfected by a scrambled siRNA in MCF10A (a) and in hTERT-HPNE (b). c–f Recordings of Magnesium Inhibited Cation (MIC) currents in MCF10A cells. c MIC currents were recorded by dialyzing the intracellular media with an intrapipette solution containing cesium (Cs) and EGTA. d Time-course of outward current densities (recorded at + 100 mV) in control (n = 15; black dots) and Cd2+ exposed cells (n = 15; red dots) during the dialyze of the intracellular media and the removal of intracellular Mg2+. e Current–voltage relationships of MIC currents built as steady state for control (n = 15; black trace) and Cd2+ exposed cells (n = 15; red trace). f Current–voltage relationships of MIC currents recorded in MCF10A transfected with a siControl (n = 8; red trace) or with a siTRPM7 (n = 8; violet trace). All results are shown as means ± SEM (effect of Cd2+ exposure: §P < 0.05; effect of siTRPM7: *P < 0.05, ***P < 0.001) (color figure online)

To evaluate if TRPM7 overexpression increased TRPM7 currents at the plasma membrane, MIC currents were recorded. An increase of the outward membrane current densities (recorded at + 100 mV) was observed during the dialysis of the intracellular media by EGTA (Fig. 1c, d). MIC current was calculated as the difference between recorded current after intracellular dialysis and the stable current recorded immediately after the patch rupture. The current–voltage (I–V) relationships of MIC currents were built (Fig. 1e, f). The I–V relationship has some of the typical features of MIC currents: inward current at negative membrane potentials, reversal potential close to 0 mV and strong outward rectification for positive membrane potentials (Prakriya and Lewis 2002). As shown in Fig. 1e, MIC currents recorded in Cd2+ cells (n = 15) were significantly larger than those recorded in control cells (n = 15; P < 0.001). When TRPM7 expression was downregulated, MIC currents were largely reduced in Cd2+ cells (n = 8) when compared to control cells (n = 8; P < 0.001; Fig. 1f). This result shows that the enhanced MIC current recorded in cells exposed to Cd2+ was mainly due to TRPM7 membrane channel activity.

Taken together, these data indicate that Cd2+ exposure increased TRPM7 expression and function leading to enhanced membrane MIC currents.

Effect of cadmium exposure on constitutive cation influxes, and Ca2+ and Mg2+ homeostasis

TRPM7 is a non-selective cation channel that mainly conducts Ca2+ and Mg2+. Moreover, TRPM7 is permeant to other metals including Mn2+ (Monteilh-Zoller et al. 2003). Using the Mn2+-quenching technique, we show that chronic Cd2+ exposure increased the constitutive cation entry in MCF10A cells (n = 164) when compared to control cells (n = 139; P < 0.001; Fig. 2a, b). TRPM7 silencing had no effect on control cells (n = 131) whereas it strongly decreased the constitutive cation entry in Cd2+ cells (n = 94; P < 0.001), indicating that TRPM7 channels are mainly responsible for enhanced constitutive cation entry in MCF10A cells chronically exposed to Cd2+. Next, we studied how Cd2+ exposure alters Ca2+- and Mg2+ intracellular homeostasis. As shown in the Fig. 2c, the cytosolic basal Ca2+ level was larger in MCF10A exposed to Cd2+ (n = 73) when compared to the control ones (n = 119; P = 0.024). Surprisingly, TRPM7 silencing increased Ca2+ level in control cells (n = 85; P = 0.018) whereas it was decreased in Cd2+ ones (n = 71; P < 0.001). The homeostasis of Mg2+ was also studied using the Mg2+-sensitive Mag-Fura-2 fluorescent probe (Fig. 2d). Chronic Cd2+ exposure increased Mg2+ level (n = 72) when compared to control cells (n = 101; P < 0.001). TRPM7 silencing had no effect on control cells (n = 77) whereas it strongly decreased Mg2+ homeostasis in Cd2+ exposed cells (n = 55; P < 0.001).

Fig. 2
figure 2

Role of TRPM7 in constitutive divalent cation entry and in cytosolic calcium and magnesium levels of human epithelial cells chronically exposed to cadmium. a Constitutive divalent cation entry estimated by manganese (Mn2+) quenching of fura-2 fluorescent probe in MCF10A cells. Mn2+-quenching was measured (at 360 nm) in control cells (n = 139; black trace), in control cells transfected with a siRNA targeting TRPM7 (siTRPM7; n = 131; grey trace), in cadmium exposed (Cd2+) cells transfected with a siControl (n = 164; red trace), and in Cd2+ cells transfected with a siTRPM7 (n = 94; violet trace). b Quantification of Mn2+-quench slope in MCF10A cells. c Estimation of cytosolic Ca2+ levels by measurement of fura-2 basal fluorescence in MCF10A cells. Fura-2 fluorescence ratio (F350/380) were measured in control cells (n = 119), in control cells transfected with a siTRPM7 (n = 85), in Cd2+ cells (n = 73), and in Cd2+ cells transfected with a siTRPM7 (n = 71). d Estimation of cytosolic Mg2+ levels by measurement of MagFura-2 basal fluorescence in MCF10A cells. MagFura-2 fluorescence ratio (F330/380) were measured in control cells (n = 101), in control cells transfected with a siTRPM7 (n = 77), in Cd2+ cells (n = 72), and in Cd2+ cells transfected with a siTRPM7 (n = 55). e Constitutive divalent cation entry estimated by Mn2+-quenching of fura-2 fluorescent probe in hTERT-HPNE cells. Mn2+-quenching was measured (at 360 nm) in control cells (n = 70; black trace), in control cells transfected with a siTRPM7 (n = 79; grey trace), in Cd2+ cells transfected with a siControl (n = 111; red trace), and in Cd2+ cells transfected with a siTRPM7 (n = 67; violet trace). f Quantification of Mn2+-quench slope in hTERT-HPNE cells. g Estimation of cytosolic Ca2+ levels by measurement of fura-2 basal fluorescence in hTERT-HPNE cells. Fura-2 fluorescence ratio (F350/380) were measured in control cells (n = 79), in control cells transfected with a siTRPM7 (n = 57), in Cd2+ cells (n = 87), and in Cd2+ cells transfected with a siTRPM7 (n = 47). h Estimation of cytosolic Mg2+ levels by measurement of MagFura-2 basal fluorescence in hTERT-HPNE cells. MagFura-2 fluorescence ratio (F330/380) were measured in control cells (n = 93), in control cells transfected with a siTRPM7 (n = 54), in Cd2+ cells (n = 103), and in Cd2+ cells transfected with a siTRPM7 (n = 63). i Cd2+ entry estimated by Cd2+-induced positive quenching of fura-2 fluorescent probe in MCF10A cells. Cd2+-quenching was measured (at 360 nm) in control cells (n = 14; black trace), in control cells transfected with a siTRPM7 (n = 14; grey trace), in Cd2+ cells transfected with a siControl (n = 23; red trace), and in Cd2+ cells transfected with a siTRPM7 (n = 23; violet trace). j Quantification of Cd2+-quench slope in MCF10A cells. All results are shown as means ± SEM (effect of Cd2+ exposure: §P < 0.05, §§§P < 0.001; effect of siTRPM7: *P < 0.05, ***P < 0.001) (color figure online)

Similar experiments were conducted in hTERT-HPNE cells. Contrarily to what we observed for MCF10A cells, the basal Mn2+-quenching slope was weak in hTERT-HPNE suggesting that there was hardly any constitutive cation entry in non-cancer pancreatic epithelial cells (Fig. 2e). As shown for MCF10A cells, Cd2+ exposure increased the constitutive cation entry (n = 111) when compared to control cells (n = 70; P = 0.004; Fig. 2e, f). TRPM7 silencing had no effect on control cells (n = 79) whereas it almost fully abolished the constitutive cation entry in Cd2+ exposed cells (n = 67; P < 0.05. Chronic Cd2+ exposure further increased Ca2+ homeostasis in hTERT-HPNE cells (n = 87) when compared to the control ones (n = 79; P < 0.05; Fig. 2g). TRPM7 silencing had no effect on control cells (n = 57) whereas it decreased Ca2+ homeostasis in Cd2+ exposed cells (n = 47; P < 0.05). Mg2+ cytosolic level was also larger in hTERT-HPNE exposed to Cd2+ (n = 103) than in control cells (n = 93; P < 0.05, Fig. 2h). In control cells, TRPM7 silencing induced a slight increase of cytosolic Mg2+ (n = 54; P < 0.05), whereas it induced a strong decrease in Cd2+ exposed cells (n = 63; P < 0.05).

TRPM7 is also involved in Cd2+ entry (Martineau et al. 2010; Monteilh-Zoller et al. 2003). Here, we used the Cd2+ quenching of Fura-2 fluorescence technique to assess the entry of Cd2+ in MCF10A cells (Fig. 2i, j). We show that MCF10A cells were permeant to Cd2+ because perfusion with 2 mM Cd2+ induced a positive quench of Fura-2 (n = 14; Fig. 2i, black trace). The quenching slope was larger in Cd2+ exposed cells indicating that these latter were more permeant to external Cd2+ (n = 23; P < 0.05. TRPM7 silencing had no effect in control cells (n = 14) while it decreased Cd2+ entry in Cd2+ exposed cells (n = 23; P < 0.05). These results suggest that TRPM7 channels provide a favorite route for Cd2+ entry in cells chronically exposed to Cd2+. However, TRPM7 silencing had no effect on Cd2+ entry in control cells suggesting that Cd2+ enters through other transporters.

To summarize, chronic Cd2+ exposure induced an increase of the divalent cation constitutive entry in both MCF10A and hTERT-HPNE cells as well as higher cytosolic levels, particularly for Mg2+. These effects were abolished by TRPM7 silencing indicating that TRPM7 overexpression induced by chronic Cd2+ exposure led to enhanced divalent cation constitutive entry and Mg2+ overload.

Effect of cadmium exposure on cell shape and organization

The membrane capacitance recorded using whole-cell patch clamp can be used to estimate membrane surface. We did not observe any difference between Cd2+ exposed (19.9 ± 1.1 pF; n = 36) and control cells (19.5 ± 1.0 pF; n = 90; P > 0.05) indicating that Cd2+ exposure did not change the lipid bilayer surface.

On the other hand, a modification of the cell organization was observed in MCF10A cells following Cd2+ exposure (Fig. 3). Indeed, control cells were organized as clusters with cell–cell junctions whereas Cd2+ exposed cells were more individualized (Fig. 3a). TRPM7 silencing restored cell clustering in Cd2+ exposed cells (Fig. 3b). However, Cd2+ exposure had no effect on hTERT-HPNE cell organization and morphology (data not shown).

Fig. 3
figure 3

Involvement of TRPM7 in MCF10A cell organization and morphology alterations induced by cadmium exposure. a Photography of control MCF10A cells and chronically exposed to cadmium (Cd2+; scale bar = 100 µm). b Photography of control and Cd2+ MCF10A cells transfected with a scrambled siRNA (siC) or with a siRNA targeting TRPM7 (siTRPM7). Yellow arrows show cells with rounded phenotype in Cd2+ cells (scale bar = 50 µm). c Quantification of circularity index in control siC cells (n = 197), in control siTRPM7 cells (n = 193), in Cd2+ siC cells (n = 192), and in Cd2+ siTRPM7 cells (n = 186). All results are shown as means ± SEM (effect of Cd2+ exposure: §P < 0.05; effect of siTRPM7: *P < 0.05) (color figure online)

Moreover, many MCF10A cells exposed to Cd2+ presented a rounded morphology (Fig. 3b, yellow arrows). Circularity index was calculated to confirm the effect of Cd2+ exposure on MCF10A morphology (Fig. 3c). Cd2+ exposed cells displayed indeed a higher circularity index (n = 192) than control cells (n = 197; P < 0.05). While TRPM7 silencing had no effect on control cell morphology (n = 193), it reversed the rounded shape induced by Cd2+ exposure (n = 186; P < 0.05).

Taken together, our data show that Cd2+ exposure induced a change in MCF10A cell phenotype which is mediated by TRPM7 channel expression.

Expression of EMT markers in MCF10A cells exposed to cadmium

The first step of the metastatic cascade is the epithelial-to-mesenchymal transition (EMT), that is characterized by morphological changes (including loss of cell–cell junctions). Here, we assessed the expression of some EMT markers in MCF10A cells following Cd2+ exposure (Fig. 4). Chronic Cd2+ exposure induced a shift of cadherin expression profile with a decrease of E-Cadherin and an increase of N-Cadherin (Fig. 4a). E/N-Cadherin switch was evaluated by E/N Cadherin expression ratio (Fig. 4b). Chronic Cd2+ exposure decreased E/N Cadherin ratio (n = 7; P < 0.05). TRPM7 silencing had no effect in control cells (n = 7; P = 0.099) but restored E/N-Cadherin ratio in Cd2+ exposed cells (n = 7; P < 0.05). Vimentin/E-cadherin ratio (V/E ratio) was also studied as a mesenchymal marker (Murai et al. 2014). V/E ratio was increased in Cd2+ exposed cells (n = 7; P < 0.05). Again, TRPM7 silencing partially restored V/E ratio in Cd2+ exposed cells (n = 7; P < 0.05), but had no effect in control cells (Fig. 4c).

Fig. 4
figure 4

TRPM7 regulates epithelial-to-mesenchymal transition induced by cadmium exposure in MCF10A cells. a Lysates from MCF10A (exposed or not to cadmium (Cd2+) and transfected with a scrambled siRNA (siControl) or targeting TRPM7 (siTRPM7) were immunoblotted with the anti-N-Cadherin, the anti-E-Cadherin, the anti-vimentin and the anti-tubulin antibodies. b E- to N-Cadherin switch quantified by E/N Cadherin expression ratio in control cells (n = 7), in control cells transfected with a siTRPM7 (n = 7), in Cd2+ cells (n = 7), and in Cd2+ cells transfected with a siTRPM7 (n = 7). c Quantification of Vimentin/E-Cadherin (V/E) expression ratio in control cells (n = 7), in control cells transfected with a siTRPM7 (n = 7), in Cd2+ cells (n = 7), and in Cd2+ cells transfected with a siTRPM7 (n = 7). All results are shown as means ± SEM (effect of Cd2+ exposure: §P < 0.05; effect of siTRPM7: *P < 0.05)

Our results show that Cd2+ exposure induced EMT in MCF10A cells. Interestingly, TRPM7 channel inhibition prevented the changes of EMT marker expression induced by chronic Cd2+ exposure.

Role of TRPM7 in cell migration and invasion induced by chronic cadmium exposure

Previous studies have shown that Cd2+ exposure transforms MCF10A cell into a basal-like phenotype, that is characterized by enhanced cell migration and invasiveness (Benbrahim-Tallaa et al. 2009). Here, we confirmed that Cd2+ exposure confers abilities to MCF10A for migration (Fig. 5a–c), since cells chronically exposed to Cd2+ were able to close a scratch wound faster than the control cells (n = 8, P < 0.01). TRPM7 silencing reduced the cell motility in Cd2+ exposed MCF10A (n = 8, P < 0.01) but not in control cells (Fig. 5a, b). We further tested the migration properties using Boyden chamber assays. We confirmed that Cd2+ exposure was able to stimulate MCF10A cell migration (n = 160; P < 0.001). TRPM7 silencing inhibited cell migration in control cells (n = 160; P < 0.001) and this effect was even more pronounced in Cd2+ exposed cells (n = 160; P < 0.001; Fig. 5c). Furthermore, we were able to confirm that Cd2+ exposure enhanced MCF10A cell invasion (n = 120; P < 0.001; Fig. 5d). As for cell migration, TRPM7 silencing inhibited cell invasion in control cells and even more in Cd2+ exposed cells (n = 120; P < 0.001). As Cd2+ exposure is also able to transform epithelial pancreatic cells into a more aggressive phenotype, cell invasion was tested in hTERT-HPNE cells (Fig. 5e). Cd2+ exposure also enhanced cell invasion in hTERT-HPNE (n = 60; P < 0.001), whereas TRPM7 silencing inhibited the cell invasion in Cd2+ exposed cells only (n = 60; P < 0.001).

Fig. 5
figure 5

TRPM7 silencing inhibits human epithelial cell migration and invasion stimulated by chronic cadmium exposure. a Photography of wound-healing assays taken immediately after the scratch (T0) or 24 h later (T24) in MCF10A cells. b Quantification of wound-healing closure in control cells transfected with a scrambled siRNA (siC) (n = 8; black dots), in control cells transfected with a siRNA targeting TRPM7 (siTRPM7; n = 8; grey dots), in cadmium exposed (Cd2+) cells transfected with a siC (n = 8; red dots), and in Cd2+ cells transfected with a siTRPM7 (n = 9; violet dots). Scratch opening measurements were normalized to the control siC condition. c MCF10A cell migration assessment after 24 h in Boyden chambers. Cell migration was normalized to the control siC condition. d MCF10A cell invasion assessment after 24 h in Boyden chambers containing Matrigel. Cell invasion was normalized to the control siC condition. e hTERT-HPNE cell invasion assessment after 24 h in Boyden chambers containing Matrigel. Cell invasion was normalized to the control siC condition. All results are shown as means ± SEM (effect of Cd2+ exposure: §§P < 0.01, §§§P < 0.001; effect of siTRPM7: ***P < 0.001) (color figure online)

TRPM7 overexpression in MCF10A

To prove the promigratory role of TRPM7 in non-cancer epithelial cells, zebrafish (ZF) TRPM7 was overexpressed in MCF10A cells. We observed that TRPM7 overexpression increased cell migration from 1 ± 0.05 to 1.49 ± 0.06 compared to control (n = 120; P < 0.001, Mann–Whitney Rank Sum Test; Fig. 6a). TRPM7 overexpression had no effect on cell viability assessed by MTT assay (n = 15; P > 0.05; Fig. 6b).

Fig. 6
figure 6

Overexpression of TRPM7 in MCF10A cells. a Cell migration in MCF10A transfected with empty vector or wild-type (WT) TRPM7 (n = 120). b Cell viability in MCF10A transfected with empty vector or WT-TRPM7 (n = 15). c Validation of TRPM7 overexpression by RT-PCR in MCF10A cells. All the results are shown as means ± SEM (***P < 0.001, Mann–Whitney Rank Sum Test)

Taken together, our results clearly show that TRPM7 is involved in enhanced cell motility and invasion induced by Cd2+ exposure, in human epithelial cells.

Discussion

Cadmium (Cd2+) is a xenobiotic which is known to induce alteration of gene expression (Chen et al. 2019; Rani et al. 2014). Cd2+ enters into the cells through numerous transporters including TRPM7 (Thevenod et al. 2019). TRPM7 is for instance involved in Cd2+ entry in HEK293 expressing TRPM7 (Monteilh-Zoller et al. 2003) and in human osteoclasts (Martineau et al. 2010). Using Cd2+-quenching of Fura-2, we show that TRPM7 is probably not involved in basal Cd2+ entry in MCF10A cells. This suggests that other transporters are required for cell Cd2+ uptake in human mammary non-cancer cells. In the present study, we show that: (1) Cd2+ exposure increases TRPM7 expression; (2) TRPM7 overexpression is followed by an increase of membrane MIC currents and constitutive divalent cation entry with increased cytosolic Mg2+ levels; (3) Cd2+ exposure induces EMT through TRPM7 expression; (4) migration and invasion are enhanced in Cd2+ exposed cells.

Previous studies have shown that Cd2+ exposure leads to non-cancer epithelial cell transformation into an invasive phenotype, by an unknown mechanism that does not involve estrogen receptor activation (Benbrahim-Tallaa et al. 2009; Qu et al. 2012). Here, we show for the first time that TRPM7 channels participate in this neoplastic transformation. TRPM7 is a non-selective cation channel which is overexpressed in numerous cancers, including breast cancer and pancreatic ductal adenocarcinoma (Dhennin-Duthille et al. 2011; Gautier et al. 2016; Middelbeek et al. 2012; Rybarczyk et al. 2012; Yee et al. 2015). We show that Cd2+ exposure induces TRPM7 overexpression in both human mammary and pancreatic non-cancer epithelial cells, suggesting that TRPM7 could be used as a biomarker for cell transformation. Interestingly, TRPM6 was downregulated following Cd2+ exposure (supplemental data) indicating that TRPM7 is probably the main Ca2+/Mg2+ channel overexpressed in cells exposed to Cd2+. Moreover, this result also confirms that TRPM7 but not TRPM6 is required for the acquisition of an invasive phenotype. To our knowledge, our study provided a unique cellular model of TRPM7 overexpression induced by pollutant exposure.

TRPM7 overexpression induces an increase of both TRPM7 membrane currents and cation constitutive entry, indicating that overexpressed TRPM7 channels are mainly located at the plasma membrane. Consequently, both cytosolic Ca2+ and Mg2+ are increased in Cd2+ exposed cells. These results are quite surprising for mammary cells because we previously showed that TRPM7 regulates Ca2+ homeostasis specifically in non-invasive MCF-7 breast cancer cells (Guilbert et al. 2009). However, TRPM7 silencing had no effect on Ca2+ and Mg2+ homeostasis in both invasive MDA-MB-231 and MDA-MB-435S cancer cells (Guilbert et al. 2013). Furthermore, it has been shown that TRPM7 regulates metastasis formation and dissemination in a model of mouse xenograft using MDA-MB-231 by regulating myosin II-based cellular tension through its kinase domain (Middelbeek et al. 2012). While TRPM7 promotes the formation of Ca2+ sparks and invadosome in neuroblastoma cells, the role of TRPM7 on invadosome formation is independent of Ca2+ influx (Visser et al. 2013). However, the role of Mg2+ has not been assessed in these two latter studies. It may be interesting to study the regulation of TRPM7 kinase domain by Mg2+ in cancer cell migration and invasion.

Regarding the pancreatic cells, we and others previously showed that TRPM7 is overexpressed in pancreatic ductal adenocarcinoma in relation with metastasis and survival (Rybarczyk et al. 2012, 2017; Yee et al. 2015). TRPM7 regulates both pancreatic cancer cell migration and invasion through Mg2+ entry (Rybarczyk et al. 2012, 2017). It has been recently shown that TRPM7 also participates to Mg2+ uptake and cell proliferation in cancer colon cells (Luongo et al. 2018). Here, we show that Cd2+ exposure was able to promote a larger increase of cytosolic Mg2+ than for Ca2+. Taken together with the recent literature, our data strongly suggest that TRPM7 overexpression leads to enhanced cytosolic Mg2+ in epithelial cells exposed to Cd2+.

MCF10A exposed to Cd2+ are not organized in clusters and have a rounded-shaped morphology, unlike control ones. These effects can be reverted by silencing TRPM7. TRPM7 knockdown can disorganize cytoskeleton in mesenchymal MDA-MB-231 which lose their spindle-shaped morphology (Middelbeek et al. 2012). TRPM7 is also known to participate in the epithelial-to-mesenchymal transition (EMT) in breast cancer (Davis et al. 2014; Kuipers et al. 2018; Sun et al. 2018). In the present study, we also demonstrate that Cd2+ exposure increases EMT markers N-Cadherin and Vimentin, while E-Cadherin is decreased in mammary epithelial cells. TRPM7 silencing dramatically reduces EMT marker expression in mammary epithelial cells exposed to Cd2+. EMT marker expression and the acquisition of a rounded-shaped morphology strongly suggest that Cd2+ exposure induces collective-to-amoeboid transition in MCF10A. Interestingly, it has been shown that other factors like hypoxia promote EMT in epithelial cancer cells leading to collective-to-amoeboid migration (Lehmann et al. 2017). However, Cd2+ exposure had no effect on pancreatic cell morphology, nor on EMT marker expression (data not shown). It has been recently shown that pancreatic cancer cells can invade and form metastasis independently of EMT (Chen et al. 2018; Zheng et al. 2015). Moreover, E-Cadherin is required for liver metastasis suggesting that EMT is not necessary for metastasis cascade in pancreatic ductal adenocarcinoma (Reichert et al. 2018).

Enhanced cell invasion in both mammary and pancreatic epithelial cells after chronic Cd2+ exposure has already been published (Benbrahim-Tallaa et al. 2009; Qu et al. 2012). Here, we show that TRPM7 regulates invasiveness in mammary non-cancer epithelial cells. Moreover, TRPM7 silencing almost fully abolishes invasiveness in Cd2+ exposed cells, suggesting that Cd2+-induced cell invasion is mainly regulated by TRPM7 in mammary epithelial cells. In non-cancer pancreatic epithelial cells, TRPM7 does not regulate cell invasion. However, TRPM7 silencing strongly reduces Cd2+-induced cell invasion as already shown for pancreatic cancer cells (Rybarczyk et al. 2017). Thus, it is tempting to speculate that TRPM7 regulated invasion specifically in cancer cell for pancreatic tissue. Taken together these results confirm that TRPM7 is a key player of invasiveness in pancreatic cancer cells. Cd2+ has been identified as a class I human carcinogen for lung cancer since 2004, but there are fewer evidence for a direct relationship between Cd2+ exposure and breast cancer/pancreatic ductal adenocarcinoma occurrence. Nevertheless, our data confirm that chronic Cd2+ exposure transforms human mammary and pancreatic epithelial cells into a more invasive phenotype. In the present study, we show that TRPM7 overexpression induced by Cd2+ exposure is linked to EMT and migration/invasion enhancement in non-cancer epithelial cells. However, there is no evidence showing that Mg2+ and/or Ca2+ entry through TRPM7 was directly involved in the acquisition of this invasive phenotype.

In conclusion, we show that TRPM7 is overexpressed in Cd2+ exposed epithelial cells, leading to Mg2+ homeostasis unbalance and enhanced cell invasiveness (Fig. 7).

Fig. 7
figure 7

Model for TRPM7 overexpression in transformed mammary epithelial cells by chronic cadmium exposure. Chronic exposure to cadmium induces TRPM7 overexpression at the plasma membrane of human epithelial cells leading to increased divalent cation influx, including cadmium influx. In transformed cells, TRPM7 regulates epithelial-to-mesenchymal transition as well as the acquisition of an invasive phenotype