Investigations on Zinc Isotope Fractionation in Breast Cancer Tissue Using in vitro Cell Culture Uptake-Efflux Experiments.

Zinc (Zn) accumulates in breast cancer tumors compared to adjacent healthy tissue. Clinical samples of breast cancer tissue show light Zn isotopic compositions (δ66Zn) relative to healthy tissue. The underlying mechanisms causing such effects are unknown. To investigate if the isotopic discrimination observed for in vivo breast cancer tissue samples can be reproduced in vitro, we report isotopic data for Zn uptake-efflux experiments using a human breast cancer cell line. MDA-MB-231 cell line was used as a model for triple receptor negative breast cancer. We determined Zn isotope fractionation for Zn cell uptake (Δ66Znuptake) and cell efflux (Δ66Znefflux) using a drip-flow reactor to enable comparison with the in vivo environment. The MDA-MB-231 cell line analyses show Zn isotopic fractionations in an opposite direction to those observed for in vivo breast cancer tissue. Uptake of isotopically heavy Zn (Δ66Znuptake = +0.23 ± 0.05‰) is consistent with transport via Zn transporters (ZIPs), which have histidine-rich binding sites. Zinc excreted during efflux is isotopically lighter than Zn taken up by the cells (Δ66Znefflux = -0.35 ± 0.06‰). The difference in Zn isotope fractionation observed between in vitro MDA-MB-231 cell line experiments and in vivo breast tissues might be due to differences in Zn transporter levels or intercellular Zn storage (endoplasmic reticulum and/or Zn specific vesicles); stromal cells, such as fibroblasts and immune cells. Although, additional experiments using other human breast cancer cell lines (e.g., MCF-7, BT-20) with varying Zn protein characteristics are required, the results highlight differences between in vitro and in vivo Zn isotope fractionation.

Introduction

Over the last decade analysis of natural metal isotopes has emerged as an interdisciplinary field of substantial biomedical potential, including diagnosis and defining disease mechanisms (1–14). Compared to metabolic studies where radioactive and single isotope tracers are used to monitor concentration changes, high precision measurements of natural isotopic fractionations are especially useful because they can inform on mechanisms, e.g., why changes in concentration occur due to altered uptake, secretion or excretion. Two recently published breast cancer studies found that Zn dyshomeostasis linked to carcinogenesis is reflected in enrichment of light Zn isotopes in malignant breast tissue compared to adjacent histologically healthy tissue (3, 5).

Zinc homeostasis in humans is maintained by multiple proteins which tightly regulate intracellular Zn concentrations. Among these proteins are influx-controlling Zn importers (SLC39A; ZIP1-ZIP14), efflux-controlling Zn transporters (SLC30A; ZnT1-ZnT10), and Zn-sequestering proteins (e.g., metallothionein) [e.g., (15–19)]. Histidine-, glutamate-, aspartate- and cysteine-residues serve as binding ligands for Zn in these proteins. In general, histidine-rich loops form the primary Zn-binding sites of ZIPs and ZnTs (20), while Zn-binding on metallothionein occurs on cysteine-rich ligands (21, 22). Cancer-induced Zn dyshomeostasis has been related to up- or down-regulation of Zn proteins (23) and changes in their coordination and ligand chemistry (24). For breast cancer, upregulation of ZIP6, ZIP7, ZIP10, ZnT2 and overexpression of metallothionein have been proposed to implicate higher Zn levels in malignant relative to healthy tissue (25–29).

Isotopic fractionation of Zn isotopes in cells is linked to coordination chemistry. Lower atomic weight is associated with lower bond energy, so heavier isotopes are enriched in the strongest ligand bonds, assuming equilibrium isotope fractionation (30–33). Thus, cysteine (S-ligands) preferentially binds isotopically light Zn, while normally tighter binding histidine (N-ligands) and aspartate (O-ligands) preferentially complex isotopically heavy Zn (30–33). If the normal Zn-ligand binding environments are disturbed by Zn dyshomeostasis, isotopic compositions in various body reservoirs (e.g., blood, urine, tissue) can be perturbed. Thus, Zn isotope analysis of clinical samples could be a method to probe the role of Zn uptake and secretion in cancer development and progression (1–3, 5, 11, 34).

In vitro cell lines are powerful simplified models to study processes at the molecular and cellular levels including the study of cancer-induced isotopic fractionation. Various human cell lines have been used to better understand uptake and transport mechanisms causing isotope fractionation of copper (35, 36), iron (37), uranium (38, 39) and Zn (38). Differentiated neuron-like cells (SH-SY5Y) preferentially incorporate and accumulate light isotopes of copper, uranium, and Zn (36, 38). Similar observations have been made for differentiated intestinal (Caco-2) cells showing the preferential uptake of isotopically light iron (37). Under oxidative stress, however, cancerous liver cells incorporate copper that is isotopically heavy (35, 40). These observations contradict the expected correlation between in vitro and in vivo studies on metal ion isotope fractionation.

There is little information about the processes at the molecular and cellular level causing systematic enrichment of light Zn isotopes in malignant breast tissue (3, 5). Our work aimed to improve the mechanistic understanding of the isotopic behavior of Zn in breast cancer, with the long-term objective of biomarker development. We conducted experiments using in vitro breast cancer cell lines (MDA-MB-231) spiked with isotopically-natural Zn and determined the isotopic fractionation induced due to active uptake and efflux. The results highlight differences between in vitro and in vivo Zn isotope processing and will inform future studies on using Zn isotope fractionation as a biomarker (Figure 1).

Figure 1

Schematic of overall approach for biomedical studies to investigate metal stable isotope fractionation.

Methods

Cell Proliferation

The cell proliferation experiments were conducted to investigate if isotopic mass controls cell growth. Isotopically enriched metal powders of 64Zn (>99.00 %) and 68Zn (>99.00 %) (Isoflex, San Francisco, CA, US) were dissolved separately in 6 M HCl. Chloride solutions were diluted with 18 MΩ cm purified H2O (MilliQ) to make 10 mM ZnCl2 solutions which were then added to the exposure medium. The exposure medium consisting of Dulbecco's modified Eagle's medium (DMEM) + 10% fetal calf serum (FBS, Gibco Life, heat inactivated) + 1% penicillin/streptomycin (Sigma-Aldrich, UK) + 10 mM ZnCl2 solution was combined, adjusted to pH 7.4 if necessary and then equilibrated overnight. This exposure medium was added over two passages of MDA-MB-231 cells to enable full isotopic equilibration. At a cell density of 1 × 105, cells were seeded in six well-plates (n = 4) for negative control (0 μM Zn, -ve), positive control (20 μM Zn, +ve), isotopically light (20 μM, 64Zn) and isotopically heavy (20 μM, 68Zn) exposures and grown under normal oxygen conditions. The cells were counted at 48, 72, and 96 h using a Cellometer Auto T4 (Nexcelom Bioscience, MA, US) to monitor the growth rate. To ensure constant cell proliferation the exposure media were refreshed after 48 h for the cells growing for 72 or 96 h.

Uptake-Efflux Experiment

Cell Culturing

In our uptake-efflux experiment, MDA-MB-231 (ATCC HTB-26) cells were first cultured under sterile conditions inside a class II biological safety cabinet at 37°C. Cells were grown in Dulbecco's modified Eagle's medium (DMEM) + 10% fetal calf serum (FBS, Gibco Life, heat inactivated) + 1% penicillin/streptomycin (Sigma-Aldrich, UK) without Zn addition. After 48 h, the cells were transferred to petri dishes containing medium without additional Zn. At this step cells were grown on microscopic slides until they reached 105 cells/cm2. Microscopic slides with MDA-MB-231 were transferred to media spiked with 15 μM isotopically natural Zn solution. A concentration of 15 μM ZnCl2 was selected as it is within the reference range of intracellular and plasma homeostatic Zn (11–24 μM) (41). As “free” Zn is toxic to cells, even at nanomolar levels (42), we added Zn complexed with histidine (Zn-His), a stable octahedral compound (43). MDA-MB-231 cells on microscopic slides were incubated for 24 h at 37°C before being transferred to the drip flow biofilm reactor for the efflux experiment.

Efflux Experiment

Figure 2 shows the experimental setup for the efflux experiment using a drip flow-through biofilm reactor (BioSurface Technologies, Corp., USA). The experiment was conducted in a class II biological safety cabinet. Prior conducting the efflux experiment with a cell line, the setup was tested for Zn blanks using a 20 mM metal-free HEPES input solution. The polysulfone drip flow-through biofilm reactor has a 10° inclination and consists of four chambers. The 20 mM metal-free HEPES input fluid was supplied through four channels imposed by a peristaltic pump and ran through the system for 1 and 2.5 h, respectively, and the effluent was collected separately for each channel.

Figure 2

Schematic of the cell uptake-efflux experiment. A HEPES solution from the Tedlar bags is pumped to the drip flow-through reactor (1 ml min−1, as controlled by a peristaltic pump). The drip flow-through reactor contains 4 microscopic slides with MDA-MB-231 cells pre-grown in medium spiked with 15 μM Zn-histidine (Zn-His). Effluent samples were taken for 0.5 h in metal-free 50 mL centrifuge tubes.

For the cell line efflux experiment, microscopic slides with cells cultured for 24 h in Zn-His spiked medium were washed three times with 20 mM metal-free HEPES solution before placed in the drip flow-through biofilm reactor to remove detached cells. This step ensured that the collected effluent solution of the efflux experiment only contained Zn excreted by the cells. The initial cell density and initial Zn isotope composition of these cells prior to the efflux experiments were determined by harvesting the cells from a microscopic slide with trypsin (Lonza). The harvested cells were counted using the Cellometer Auto T4 cell counter (Nexcelom, Bioscience). The 20 mM metal-free HEPES solution was directly dripped onto the cells with a rate of 1 ml min−1. Due to inclination of the drip flow-through reactor the residence time of the fluid passing the chambers was extremely short (<1 min) resulting in short contact time of the fluid and the cells. The effluent was collected in 30 min intervals in acid-washed metal-free centrifuge tubes (VWR). After the experiment, the cells on the microscopic slides were harvested using trypsin. Harvested cells, effluent and aliquots of the medium, and starting input HEPES solution were stored at −20°C, until required.

Prior to isotope analysis, harvested cells were lysed using a MARS 5 Digestion Microwave System (CEM Corp., UK). The harvested cell suspension was transferred in acid cleaned XP-1500 Plus (PTFE) vessels and 3 ml of quartz sub-boiled distilled nitric acid (15.4 N) and 2 ml hydrogen peroxide (Romil Ltd) were added. A blank sample consisting of only 3 ml of nitric acid and 2 ml hydrogen peroxide was processed for quality control. The samples were pre-digested at room temperature overnight. Microwave-assisted acid digestion was performed by ramping up the temperature stepwise to 210°C and 250 psi over 60 min, and held there for 30 min to ensure complete digestion.

For Zn isotope analysis, all samples (e.g., acid-digested cells, efflux solution) were double-spiked consisting of mixture of 64Zn and 67Zn. The double-spike mixture had a 64Zn/67Zn = 4.2336 and was added to each sample in an ideal sample/spike mixture of 1. To ensure the ideal sample/spike mixture, a small sample aliquot (20 μl) was spiked with an appropriate volume of 64Zn + 67Zn spike and the sample Zn concentration was determined by isotope dilution. Once sample-Zn was known, a larger aliquot of the sample was prepared for Zn isotopic analysis. Zinc was separated and purified using AG-MP1 resin (BioRad, 100–200 mesh) following the method described by (1, 44). The Zn columns with 250 μl AG-MP1 resin were rinsed with 10 mL of each 0.1N HNO3 and double-deionized water, conditioned with 6N HCl, and equilibrated with 4 x 0.5 ml 1N HCl. The spiked sample re-dissolved in 1 ml 1N HCl was loaded on the column and subsequently rinsed with 8 mL of 1N HCl. In the last step, 6 mL of 0.01N HCl was added to elute Zn from the column. Procedural blanks and Zn standard solutions (IRMM-3702, London-Zn) were processed in the same way as samples. The Zn isotope ratios and precise elemental Zn concentrations were measured with Nu Plasma HR-multiple-collector ICP-MS by double-spike and sample-standard bracketing techniques. All Zn isotope values were expressed as delta notation [δ66Zn (‰)] relative to the JMC-Lyon:

The extent of isotope fractionation was determined and described as:

where Δ66Znuptake reflects Zn isotope fractionation for Zn taken up by cells (δ66Zncells) from Zn-His spiked medium (δ66Znmedium) after incubation for 24 h, and Zn isotope fractionation of the efflux (Δ66Znefflux) is described as:

based on the difference between Zn isotopic composition of cells after incubation for 24 h (δ66Zncells) and the Zn isotopic composition of the effluent (δ66Zneffluent) collected for 0.5 h.

The in-house London Zn was used as standard for sample-standard bracketing and the precision of the unprocessed London Zn solution was ± 0.04‰ (n = 26) over two analytical days. The precision of δ66Zn for processed pure Zn IRMM-3702 reference material was −0.13 ± 0.08‰ (2SD, n = 6). If not stated otherwise, the precision of 2SD for δ66Zn (‰) refers to repeated analysis of London Zn bracketing standards.

Results

Zinc Isotope Enrichment Proliferation Tests

Figure 3 shows the average cell growth of MDA-MB-231 in a 64Zn and 68Zn isotope enriched medium as well as the control without Zn (-ve) and natural Zn (+ve). In all experiments, regardless of the isotopic composition of the Zn in solution, the cell density increased by 5-fold within 96 h reaching an average density of 5.3 × 105 ± 0.5 × 105 cells. If cell growth is driven by the mass of an isotope, we would expect to see faster growth for the lighter isotopes compared t the heavier isotopes of Zn (Figure 3 and Table 1).

Figure 3

Zinc isotope enrichment proliferation after 96 h under different medium Zn conditions: -ve control = negative control without Zn; +ve control = positive control with 20 μM “natural” Zn; 64-Zinc = 20 μM 64Zn and 68-Zinc = 20μM 68Zn only.

Table 1

Data summary for Zn isotope enrichment proliferation tests.

Experimental condition	Time (hrs)	Average (x105 cells)	Stdev (x105 cells)
-ve control	0	1.00	–
-ve control	48	1.82	0.66
-ve control	72	3.30	0.93
-ve control	96	5.94	1.49
+ve control	0	1.00	–
+ve control	48	1.42	0.40
+ve control	72	3.39	1.74
+ve control	96	5.20	0.79
64Zn	0	1.00	–
64Zn	48	1.86	0.05
64Zn	72	3.21	0.99
64Zn	96	5.35	1.61
68Zn	0	1.00	–
68Zn	48	0.93	0.54
68Zn	72	2.98	1.62
68Zn	96	4.60	1.68
-ve control	96	5.94	1.49
+ve control	96	5.20	0.79
64Zn	96	5.35	1.61
68Zn	96	4.60	1.68

Cell Efflux

Control of blank levels. We assessed blank contributions because Zn usually has a high background level, which can significantly alter the measured isotopic composition. The microwave control blank was 0.9 ng Zn ml−1 and procedural blanks for the anion exchange chromatography were 0.24 ± 0.20 ng Zn ml−1 (n = 4). The Zn blank for the experimental setup ranged between 0.5 and 2.7 ng Zn ml−1 after 1 h effluent collection and decreased to 0.6–1.4 ng Zn ml−1 after 2.5 h (Table 2). This corresponds to 0.7–4.8% of Zn (mean: 1.9 ± 2.0% Zn) collected in the effluent of the cell line efflux experiment. The sum of all blanks (microwave, procedural and efflux) corresponds to 0.008 and 0.81% of the initial Zn in medium and harvested cells before the efflux experiment. This observation demonstrates that precise and reliable Zn isotope data can be achieved for cell line experiments using the drip flow-through biofilm reactor.

Table 2

Data summary for uptake-efflux experiment.

Experiment	Sample type	Time (h)	Flow reactor channel	Zn concentration (ng/ml)	% Zn rel. efflux	Number of analysis
Blank test	Fluid	1	1	NA		–
	Fluid	1	2	0.48		1
	Fluid	1	3	2.09		1
	Fluid	1	4	2.66		1
Blank test	Fluid	2.5	1	0.68		1
	Fluid	2.5	2	1.38		1
	Fluid	2.5	3	0.55		1
	Fluid	2.5	4	NA		–
Experiment	Sample type	Time (h)	Flow reactor channel	Zn concentration (ng/ml)	δ66ZnIRMM3702 (‰)	Number of analysis
Efflux	Fluid	0.5	1	105	−0.19 (0.05)	1
	Fluid	0.5	2	40	−0.07 (0.04)	1
	Fluid	0.5	3	70	−0.11 (0.05)	1
	Fluid	0.5	4	55	NA	1
Initial Zn medium	Fluid	0			−0.03 (0.07)	2
Cells before efflux	Solid	–		302	+0.2 (0.07)	1
IRMM-3702 (processed)					−0.13 (0.08)	6

Mass balance. Our current experimental design does not allow for a complete isotopic mass balance. Cellular concentrations and isotopic changes of Zn can only be monitored assuming a similar distribution and behavior of cells on each microscopic slide. Based on the Zn concentration of harvested cells from two microscopic slides, we calculated that about 0.15% of total 15 μM Zn from the medium was incorporated into MDA-MB-231 cells. With an initial MDA-MB-231 cell density of 2.6 × 106 ± 0.62 × 106 cells (n = 2), the Zn uptake yielded in 8.9 × 10−6 nM Zn/cell (= accounting for 5.4 × 109 Zn atoms/cell). The initial medium spiked with isotopic naturally-distributed 15 μM Zn-His had a δ66Zn of −0.03 ± 0.02‰ (2 SD, n = 2). MDA-MB-231 cells treated with 15 μM Zn-His had a △66Zn of +0.20‰ prior to the efflux experiment. The intracellular Zn of MDA-MB-231 is isotopically heavier compared to the initial medium (Δ66Znuptake = +0.23 ± 0.1‰; 2SD, Figure 4). Zinc in the effluent collected from four channels corresponds to 4.03 ± 1.66 × 10−7 nM Zn/cell (= accounting for 2.4 × 108 Zn atoms/cell). The δ66Zn of the efflux solutions collected from the four channels range between −0.07 and −0.19‰ with an average of −0.12 ± 0.12‰ (2 SD).

Figure 4

Possible pathways that may lead to Zn isotope fractionation (Δ66Zn) based on the uptake-efflux experiment with the MDA-MB-231 cell line. Zn was provided as Zn complexed with histidine (Zn-His). Zn isotope fractionation by cell uptake (Δ66Znuptake) can be caused by (1) active import of Zn-His via Zn importer proteins (ZIP) or (2) non-quantitative diffusion of Zn via unspecified channels. Zn isotope fractionation via efflux (Δ66Znefflux) of Zn out of the cells can be caused by (1) active export of Zn-His via Zn transporter proteins (ZnT) or (2) non-quantitative diffusion of Zn via unspecified channels.

Discussion

The results show nearly identical proliferation of MDA-MB-231 under single isotope and natural Zn exposure conditions; these observations demonstrate that isotopic mass does not control cell growth. The addition of 20 μM Zn did not enhance or suppress the proliferation, as the cell density for light (64Zn) or heavy (68Zn) isotopically spiked Zn conditions are similar to those for the no-Zn control and Zn with natural isotope composition (Figure 3).

Metabolism in cancerous breast tissue is generally thought to be enhanced in order to provide sufficient energy and anabolic substrate for proliferation (45, 46). Adaptation in Zn-protein expression in neoplasms is not driven by preferential uptake of any Zn isotopologue. This means that the isotopic shift reported in breast tumors (3, 5) likely reflects metabolic changes for Zn in cancerous breast cells, and not an energetic advantage related to the isotopic composition possessed by the neoplasm. As such, observed isotopic shifts in breast cancer tissue can be used purely as a monitor of neoplastic driven metabolic changes. Furthermore, single stable isotope tracer studies are established as providing better insight into the processes under investigation without inducing or exacerbating pathological processes (47).

Notably, our MDA-MB-231 cell line analyses show isotopic fractionations in the opposite direction to those determined for in vivo breast cancer tissue samples. Malignant breast tumors have significantly lower δ66Zn than healthy tissues which has been interpreted as an indicator for increased metallothionein expression (3, 5). The difference in Zn isotope fractionation observed between in vitro MDA-MB-231 cell line experiments and in vivo breast tissue might be due to differences in Zn transporter levels or intercellular Zn storage mechanisms (e.g., metallothionein expression).

A positive direction of Zn isotope fractionation for the cell uptake (Δ66Znuptake = +0.23 ± 0.1‰, 2SD, n = 2) can be mechanistically explained by Zn adsorption, diffusion through non-selective channels into the cell or an active cell uptake via ZIP proteins. Adsorption studies on bacterial cells and biofilms show preferential binding of heavy Zn isotopes with equilibrium isotopic fractionation between +0.50 and +1.30‰ (48, 49). However, adsorption may be less likely in our studies as the 2:1 Zn (His)2 complexes are uncharged (43) and cells were rinsed with metal-free HEPES to remove any loosely adsorbed Zn. Instead, the preferential uptake of 66Zn by MDA-MB-231 may reflect previously reported upregulation of histidine-rich ZIPs in breast cancer cells (25–29). ZIP proteins (ZIP1-ZIP14) control intracellular Zn levels by transporting Zn into the cells across the cell membrane. Upregulation of ZIPs has been observed in most cancers to meet the demand of increased rate of proliferation and metabolism (46). An indicator for ZIP upregulation is one-order of magnitude higher Zn uptake of 5.4 × 109 Zn atoms/cells by MDA-MB-231 in our experiments compared to Zn required by mammalian cells (108 Zn atoms/cell) (50).

Zinc isotope fractionation for Zn uptake by cells can be associated with Zn speciation, Zn-binding domains of proteins or exchange reactions. A study with neuron-like cells showed a preferential uptake of light 64Zn with an isotopic fractionation of 1.14‰ (39). However, uptake and efflux of Zn-His should not fractionate Zn isotopes due to a lack in coordination change for complexed Zn species (51). Caldelas and Weiss (51), however, argued that non-quantitative uptake of complexed Zn can result in an enrichment of 66Zn in plant cells. We can assume a similar isotopic effects of non-quantitative uptake for Zn-His because only 0.15% of Zn from the medium has been taken up by MDA-MB-231 cells. However, our data cannot directly resolve whether diffusion through non-selective channels or uptake via ZIPs, or other proteins, is the predominant pathway for Zn uptake by MDA-MB-231 cells.

Negative Δ66Znefflux suggests diffusion via non-selective channels. During diffusion light 64Zn diffuses faster out of the cell than 66Zn. The efflux experiment shows that MDA-MB-231 cells preferentially 'excrete' light 64Zn with Δ66Znefflux of −0.42 ± 0.12‰ (Figure 4). An existing concentration gradient with high Zn in cells and no Zn in the efflux solution may allow diffusion of Zn-His via non-selective channels. Coutaud et al. (49) reported preferential excretion of lighter Zn isotopes with Δ66Zn of −0.50 ± 0.20‰ due to desorption from the biofilm. Also, lighter Zn isotopes can diffuse faster through biological membranes or due to a gradient within the boundary layer (52). This leaves the cell culture enriched in heavy 66Zn over time resulting in opposite isotopic signal of what has been observed for in vivo malignant breast tissue. These in vivo studies found that Zn isotope compositions are about 0.17‰ lower in malignant breast than healthy tissues (3, 5).

The ZnTs family (ZnT1-ZnT10) regulates intracellular Zn levels by transporting Zn from the cytoplasm into the extracellular space. To decrease the toxic effects of high Zn in breast cancer cells, ZnT1 is overexpressed to prevent apoptotic cell death and increase Zn secretion. Wang et al. (53) observed that the mRNA expression level of ZnT-1 increases significantly in MDA-MB-231 cells when they are exposed to Zn sulfate. In vivo, ZnTs transport only free Zn2+ across cell membranes (54) and with that most likely isotopically heavy Zn2+ out of the cell due to their histidine-rich structure. In our experiments, however, Zn-His might be able to pass cell membranes via different channels and the diffusion out of the cell is similar to Zn-His uptake. This explains the preferential efflux of light 64Zn due to faster diffusion of light 64Zn. Obviously conducting an uptake-efflux experiment with free Zn2+ could be a better approach, but even if Zn is added to the medium as free Zn2+, culture medium contains millimolar levels of ligands (e.g., histidine, cysteine and phosphate) which rapidly complex free Zn2+ (42). Moreover, conducting an experiment with “free” Zn can result in apoptosis because “free” Zn is a cell toxin at even nanomolar levels (42).

In vivo, tumors are heterogenous, have a three-dimensional shape and are often mixed with connective tissue, immune cells, and stromal tissues (55–57). More specifically, breast cancer varies in its origin and genetic lesions, leading to distinct phenotypes (58, 59) which is also reflected in different ZIP and ZnT expression between various breast cancer cell lines. Therefore, we tested our findings from MDA-MB-231 in a large panel of breast cancer cell lines (60) to further our understanding in cell lines representing patient subtypes (59) (Table 3). Using a well-defined list of human ZnT and ZIP proteins (23), we identified ten Zn proteins which are differentially expressed (FDR-adjusted P < 0.1) in MDA-MB-231 compared to other basal-like breast cancer cell lines (e.g., Her2, Luminal A and Luminal B) (Table 3). It shows that more ZIPs and ZnTs are up- than downregulated which might explain the high uptake and efflux. Other studies on the hormone-dependent breast cancer cell lines T-47D and MCF-7 showed that ZnT2 and metallothionein are overexpressed, providing protection from Zn hyperaccumulation and preventing apoptosis by either removing Zn from the cell or redistributing it among cellular compartments (27, 29). Consequently, whether breast cancer cell lines and the experimental approach with initially ligand complexed Zn (Zn-His) reflect that of in vivo carcinoma remains an important issue to resolve before drawing any reliable conclusion on Zn isotope fractionation in clinical samples.

Table 3

Differential gene expression analysis for ZnTs and ZIPs between the basal-like (including MDA-MB-231) and other breast cancer cell lines (Her2, Luminal A and Luminal B).

Protein coding gene	Proteins	logFC	P-value	Q-value	Direction of expression
SLC39A6	ZIP6	−1.172	0.001	0.007	Down
SLC39A9	ZIP9	−0.521	0.000	0.001	Down
SLC39A11	ZIP11	−0.899	0.000	0.003	Down
SLC30A9	ZnT9	−0.318	0.036	0.095	Down
SLC39A4	ZIP4	0.899	0.012	0.041	Up
SLC39A8	ZIP8	0.952	0.010	0.040	Up
SLC39A10	ZIP10	0.430	0.040	0.095	Up
SLC39A14	ZIP14	0.938	0.010	0.040	Up
SLC30A3	ZnT3	0.769	0.006	0.033	Up
SLC30A6	ZnT6	0.243	0.036	0.095	Up
SLC30A5	ZnT5	−0.226	0.068	0.148	Normal
SLC30A10	ZnT10	0.113	0.079	0.158	Normal
SLC39A12	ZIP12	0.021	0.137	0.254	Normal
SLC39A2	ZIP2	0.373	0.153	0.262	Normal
SLC30A4	ZnT4	−0.207	0.320	0.510	Normal
SLC39A7	ZIP7	−0.164	0.356	0.510	Normal
SLC39A1	ZIP1	−0.162	0.361	0.510	Normal
SLC39A13	ZIP13	−0.194	0.419	0.558	Normal
SLC30A8	ZnT8	0.146	0.522	0.659	Normal
SLC30A2	ZnT2	0.095	0.555	0.667	Normal
SLC39A5	ZIP5	−0.027	0.653	0.706	Normal
SLC30A7	ZnT7	0.062	0.687	0.706	Normal
SLC39A3	ZIP3	0.061	0.701	0.706	Normal
SLC30A1	ZnT1	0.081	0.706	0.706	Normal

FC, fold change; For each gene, log2FC was calculated as the difference in means of mRNA abundance between the basal-like and other cell lines. Welch's t-test was used to quantify statistical significance in the difference in mRNA abundance of the two groups. P-values were adjusted for multiple comparisons (Q-values) using the Benjamini-Hochberg method. P-values < 0.05 are considered being significant. Blue: Downregulated ZnTs and ZIPs in MDA-MB-23 relative to other basal-like breast cancer cell lines. Green: Upregulated ZnTs and ZIPs in MDA-MB-23 relative to other basal-like breast cancer cell lines.

Future Directions

Our results support the relevance of in vitro cell culture models to identify and understand factors driving changes in Zn isotopic composition caused by cancer development and progression. However, with respect to in vivo relevance, the results should be regarded as preliminary, and additional experiments are required using other human breast cancer cell lines, e.g., MCF-7, BT-20, with different expressions of ZnT and ZIP (Table 1). This is important, because, for instance, upregulation of ZIP6 in estrogen treated MCF-7 (61) (Table 1) might lead to more pronounced uptake of heavy 66Zn. As almost every cancer is unique, it may be advantageous to determine average Zn isotope fractionation for cancer cells based on Zn biology (e.g., MT expression, ZnT and ZIP regulation, cell growth rate) moreover the use of 3D cell cultures as spheroids may be a representative approach and enable other stromal cell types to be incorporated into the workstream.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

KS and FL designed the experiments. KS performed the data analysis, interpretation of the data, and drafting of the manuscript. ALH provided lab space at MRC Weatherall Institute of Molecular Medicine to conduct the uptake-efflux experiment. The isotopic analyses were conducted in the Oxford Earth Sciences isotope facility of ANH. SH provided the bioinformatic data and wrote the bioinformatic parts of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was funded by Cancer Research UK (No. C5255/A18085), through the CRUK Oxford Centre and John Fell Oxford University Press Research Fund (No. 153/035).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.