High-resolution spatiotemporal pHe and pO2 imaging in head and neck and oesophageal carcinoma cells
Cancer & Metabolism volume 9, Article number: 21 (2021)
pO2 and pH are physiological parameters relevant for different processes in health and disease, including wound healing and cancer progression. Head and neck squamous cell carcinomas (HNSCC) and oesophageal squamous cell carcinomas (ESCC) have a high rate of local recurrence that is partly related to treatment-resistant residual tumour cells. Hence, novel diagnostic tools are required to visualise potential residual tumour cells and thereby improve treatment outcome for HNSCC and ESCC patients. We developed a device to spatiotemporally measure oxygen consumption rates (OCR) and extracellular acidification rates (ECAR) to distinguish HNSCC and ESCC cells from healthy cells in vitro, exploiting general metabolic differences between cancer cells and healthy cells.
OCR and ECAR were measured via a newly developed device named STO2p-Q (SpatioTemporal O2 and pH Quantification) using the VisiSens technology based on ratiometric fluorescence imaging, facilitating spatiotemporal resolution. Results were confirmed using extracellular flux analyses (Seahorse technology).
STO2p-Q is described and used to measure OCR and ECAR in HNSCC and ESCC cell lines and normal fibroblast and epithelial cells as components of the tumour microenvironment. OCR measurements showed differences amongst HNSCC and ESCC cell lines and between HNSCC/ESCC and normal cells, which on average had lower OCR than HNSCC/ESCC cells. Both OCR and ECAR measurements were independently verified using the Seahorse technology. Additionally, using STO2p-Q, HNSCC/ESCC, and normal cells could be spatially resolved with a resolution in the low millimetre range.
We developed a method to spatiotemporally measure OCR and ECAR of cells, which has many potential in vitro applications and lays the foundation for the development of novel diagnostic tools for the detection of cancerous tissue in HNSCC and ESCC patients in vivo.
Head and neck squamous cell carcinomas (HNSCCs) have a worldwide incidence of over 550,000 cases annually and cause over 380,000 deaths per year . Despite improvements in surgical and adjuvant treatment options including laser- and robotic-based surgery , monoclonal antibody therapies targeting EGFR and PD-1, and molecular subtyping options [3, 4], the 5-year survival rate for HNSCCs remains below 50% . Thus, the clinical outcome of HNSCC patients is persistently poor in comparison to other tumour entities and further therapy improvements are in high demand.
The standard treatment procedure for HNSCC patients consists of complete tumour resection, while preserving functionally essential anatomical structures. Advanced stages of the disease usually require a combination of surgery with irradiation and/or chemotherapeutic drugs such as 5-fluorouracil and cisplatin. For the surgical removal of the tumour, a resection margin of minimum 5 mm is recommended, and incomplete tumour removal is associated with worsened survival rates . Rates of positive surgical margins (PSMs), i.e. incomplete tumour resection, range from 15 to 30% in head and neck cancer patients [6, 7]. Incomplete tumour resection can result from topological constraints and residual tumour cells remaining undetected during surgery. That is why further improving surgical resection techniques is an important clinical aim in addition to better diagnostic tools to support therapy planning by distinguishing cancerous tissue from healthy tissue.
One of the main characteristics of cancer cells is their highly proliferative nature, which requires them to adapt their metabolism to produce the necessary energy and biomass to create new progeny . Otto Warburg observed that cancer cells, even in the presence of oxygen, preferentially metabolise glucose through “aerobic” glycolysis instead of oxidative phosphorylation [9, 10]. This allows cancer cells to survive and proliferate in both aerobic and anaerobic conditions, which is especially relevant in solid tumours, comprised of regions of high and low oxygenation. Additionally, aerobic glycolysis indirectly provides metabolites to produce cellular building blocks, such as fatty acids and non-essential amino acids . However, glycolysis results in the formation of large amounts of lactate, mostly a waste product, which cancer cells extrude at higher rates than normal cells by upregulating acid extruding transporters, e.g. NHE1 and monocarboxylate transporters . Through this mechanism, cancer cells maintain a higher intracellular pH (pHi) than normal cells, which contributes to efficient proliferation, migration, and metabolism . As a further result, the extracellular pH (pHe) of cancer cells in tumours is measurably lower compared to non-malignant tissues with average ranges of pH 6.8–7.0 for tumour vs. 7.4 for normal cells [13, 14].
Both oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) are indicators of cellular metabolic activity. They are usually determined by quantifying decreases in oxygen concentration and pH of the medium surrounding the cells over time, respectively. However, in order to distinguish cancerous from healthy tissue, a diagnostic tool must provide spatial resolution of OCR and ECAR within a tissue. In contrast to electrode-based or soluble sensor-based tools, immobilised pH- or O2-sensitive fluorophores are suitable options for such an application. The VisiSens TD system (PreSens GmbH) in combination with pH- and O2-sensitive sensor foils facilitates spatiotemporal monitoring of pH and O2 concentrations, both in vitro and in vivo [15,16,17,18]. In the present study, we developed a device based on the existing VisiSens TD camera system [19, 20] that allows spatiotemporal pO2 and pH quantification (STO2P-Q). Our results show that the newly developed measurement system can detect differences in OCR and ECAR between HNSCC/ oesophageal squamous cell carcinomas (ESCC) cell lines and normal cells, which lays the foundation for the development of diagnostic tools based on the VisiSens technology.
mCherry was PCR amplified and cloned into the pCAG-3SIP plasmid using EcoRI and NheI restrictions sites. Additionally, the promotor was changed from CMV to EF1.
Cell culture and generation of transgenic cell lines
FaDu, Kyse30, Cal33, and Cal27 were purchased from DSMZ (Braunschweig, Germany), HNEpC from PromoCell (Heidelberg, Germany). HFF-1 cells were initially provided by Prof. Markus Meissner and were subsequently purchased from ATCC (Manassas, USA) (SCRC-1041). FaDu cells were initially isolated from squamous cell carcinoma of the hypopharynx . Kyse30 cells are derived from a well-differentiated invasive oesophageal squamous cell carcinoma . Cal27 were established from a poorly differentiated squamous cell carcinoma of the tongue, and Cal33 from a moderately differentiated squamous cell carcinoma of the tongue . HNEpC are primary nasal epithelial cells isolated from healthy individuals, HFF-1 are fibroblasts isolated from the foreskins of new-born human males. All cell lines were cultured in T25 or T75 cell culture flasks (Sarstedt AG, Nümbrecht, Germany) at 37 °C, 5% CO2, and 100% humidity. Cells were grown to 70–90% confluence and then split 1:3 to 1:10, one to three times a week. FaDu, Cal33, Cal27, and HFF were cultured in DMEM (Gibco, Dublin, Ireland) supplemented with 10% FBS (FBS Superior, Sigma Aldrich, St. Louis, USA) and 1% penicillin/streptomycin (final concentration 100 μg/ml, Gibco, Dublin, Ireland). Kyse30 cells were cultured in RPMI 1640 with l-glutamine (Gibco, Dublin, Ireland) supplemented with 10% FCS and 1% penicillin/streptomycin. HNEpC were cultured in Airway Epithelial Medium containing Supplement Mix according to the manufacturer’s instructions (PromoCell, Heidelberg, Germany).
Kyse30 and FaDu cells were transfected with a pCAG-EF1-mCherry expression plasmid using the MaTra Kit (PromoCell, Heidelberg, Germany) according to manufacturer’s instructions. Two days after transfection, transfected cells were selected by adding 1 μg/ml puromycin to the culture medium. After 2 to 3 weeks, the selected cells were resistant to puromycin and mCherry expression was verified by fluorescence microscopy.
O2 and pH measurements using the newly developed device
Short-term measurements: standard procedure
Each separate batch of sensor foils (PreSens, Regensburg, Germany) was separately calibrated. O2 (SF-RPSu4-NAU-L4/W4-OIW) and pH (SF-HP5R-L4/W4-US) sensor foils were provided by PreSens, Regensburg, Germany. Calibration of O2 sensor foils was performed using a saturated solution of Na2SO3 in H2O for the 0% O2 calibration point and air saturated medium for the 100% O2 calibration point. pH calibration was performed using NaH2PO4/Na2HPO4 buffers with a buffer capacity of 40 mM, substituted with NaCl to a total ionic strength of 140 mM with six different pH values between pH 5 and pH 8.5.
To measure the oxygen consumption and extracellular acidification of individual cell lines, 4.5 × 104 cells were seeded into one well of 2-well silicon inserts (Ibidi, Gräfelfing, Germany) inserted into the wells of a 6-well plate (Sarstedt AG, Nümbrecht, Germany), and cells were allowed to attach overnight. The silicon inserts were removed and 2 ml imaging medium (RPMI 1640 without Phenol Red, with l-glutamine) (Gibco, Dublin, Ireland) supplemented with 10% FCS and 1% penicillin/streptomycin added to each well. The 6-well plate was then transferred to an incubator containing the STO2p-Q device, placed on the frame and the lid replaced by a lid equipped with sensor foil-covered plungers. The sensor foils were equilibrated in PBS for at least 1 h before each measurement. Cells were equilibrated to the incubator’s environment for 1 h before measurement was started. To start measurements, plungers were lowered onto the cells and the VisiSensVS Software (PreSens) programmed to acquire images in 10-s intervals for 30 min (pH measurements) or 20-s intervals for 60 min (pO2 measurements). The exposure time was set between 400 and 600 ms, the gain was set to 3, and LED intensity set to maximum.
After measurement, cells were stained with 1 μg/ml Hoechst 33342 fluorescent DNA staining dye for 10 min. Microscopic images over the entire area of the cell layer were taken in the DAPI channel (excitation 350/50, emission 460/50) using the tile scanning function of the LASX software and a Leica DMi8 microscope (Leica, Wetzlar, Germany) with a × 10 objective (HC PL FLUOTAR 10x/0.32 PH1). Nuclei were counted as described below to normalise OCR and ECAR for cell numbers.
Short-term measurements: inhibitor experiments
Oligomycin (Abcam, ab141829, Cambridge, UK), 2-Deoxy-d-Glucose (2-DG) (Abcam, ab142242, Cambridge, UK), glucose (Merck, Darmstadt, Germany), and Carbonyl cyanide-4-(trifluoromethoxy)phenylhydrazone (FCCP) (Abcam, ab120081, Cambridge, UK) were diluted in DMSO or ddH2O according to manufacturer’s instructions. Cells were prepared according to the standard procedure described above. Drugs were added to the imaging medium and the plate incubated for 20 min on the device’s frame before starting the measurements. pO2 measurements were performed in 20-s intervals for 30 min, and pH measurements in 10-s intervals for 30 min. After the measurements, nuclei were stained and microscopic images acquired as described above.
To create two-dimensional (2D) cell landscapes with defined shapes and areas, 4-well silicon inserts and 4-well silicon μ-inserts (Ibidi, Gräfelfing, Germany), respectively, were inserted into 6-well plates. 7.5 × 104 cells were seeded into each well of the larger silicon inserts, and between 1 × 103 and 8 × 103 cells were seeded into the smaller silicon inserts. Cells were allowed to attach over night before the inserts were removed and 2 ml imaging medium was added to each well of the 6-well plate. Before measurements were started, plates were equilibrated to the incubator environment for 1 h. pH and pO2 were measured in 20-s intervals for 60 min, as described above.
In a second approach 200–800 FaDu or Kyse30 cells were seeded into the wells of a 6-well plate and allowed to grow into visible colonies for 7–14 days. Two millilitres of imaging medium was added to each well before microscopy, pH and pO2 measurements were performed as described above (20-s intervals, 60-min duration).
To create mixed 2D cell landscapes, small drops (5–20 μl) of cell suspensions of both HFF and FaDu (pCAG-EF1-mCherry) or Kyse30 (pCAG-EF1-mCherry) with 1 × 106–2 × 106 cells/ml were pipetted into the wells of a 6-well plate to create different patterns. After a few minutes, the drops were connected using a 10-μl pipet tip. Approximately 0.5–1 ml of medium was carefully added to the edge of the well around the drops to prevent evaporation. Cells were allowed to adhere overnight, and then 2 ml of imaging medium was added to each well. Microscopic images over the entire area covered by cells were taken in phase contrast, and in the TXR channel (excitation 560/40, emission 639/75) with a × 5 or × 10 objective using the tile scanning function of the LASX software and a Leica DMi8 microscope. pH and pO2 measurements were performed as described above (20-s intervals, 60-min duration).
In a different approach 200–800 FaDu (pCAG-EF1-mCherry) or Kyse30 (pCAG-EF1-mCherry) cells were seeded into the wells of a 6-well plate and allowed to grow into visible colonies for 7–14 days. Subsequently, 500,000 HFF were seeded into each well containing colonies. The cells were allowed to adhere overnight and then 2 ml of imaging medium was added to each well. Microscopy, pH and pO2 measurements were performed as described above (20-s intervals, 60-min duration).
In a third approach, 500,000 HFF were seeded into the wells of a 6-well plate, allowed to grow for 3 days and subsequently irradiated with 35 Gy using a CIX2 cabinet irradiator (Xstrahl, Camberley, UK) equipped with a 0.5-mm Cu filter. Two days later 200–800 FaDu (pCAG-EF1-mCherry) or Kyse30 (pCAG-EF1-mCherry) cells were seeded onto the HFF cell layers and allowed to grow into visible colonies for 7–14 days. Two millilitres of imaging medium was added to each well before microscopy, and pH and pO2 measurements were performed as described above (20-s intervals, 60-min duration). For all colony-based experiments colony sizes were measured using Fiji .
Self-adhesive pH sensor spots (PreSens) were inserted into 96-well plates and pre-treated with 100 μl FCS for 1 h at 37 °C before cells were seeded to improve cell adhesion. For O2 measurements, 96-well plates with the entire bottom consisting of O2 sensor foil material were equilibrated with 100 μl PBS for 1 h at 37 °C. PBS and FCS, respectively, were removed and pH sensors washed with PBS once before seeding cells.
To compare the pH and pO2 in FaDu cell cultures depending on cell numbers, 5000, 10,000, 20,000, 40,000, and 60,000 FaDu cells were seeded into 96-well plated equipped with pH sensor spots or with the bottom consisting of O2 sensor foil. For each cell number, six wells were used, i.e. six technical replicates. To compare cell lines, 40,000 cells of FaDu, Kyse30, Cal33, Cal27, HFF, and HNEpC were seeded into six wells each (six technical replicates). O2 and pH calibration buffers were included for each separate batch of sensor foils, as described above. The 96-well plates containing the cells were sealed with a gas permeable foil (Biozym, 600295, Vienna, Austria) to prevent evaporation and left in the laminar flow cabinet for 1 h to improve equal distribution of cells in the wells. Plates were then placed on the frame and pH or O2 measured for 90 h in 1-h intervals.
Extracellular flux analysis (Seahorse XFe96)
Approximately 20 h before measurement, 4 × 104 cells/well were seeded into XF96 microplates in five technical replicates per cell line. One hour before the assay, culture medium was replaced by Seahorse assay medium (XF RPMI 103576-100, Agilent, Santa Clara, USA) containing 2 mM l-glutamine and 5 mM glucose (for OCR measurements). Basal OCR and ECAR, and ECAR measurements after injection of glucose (5 mM final concentration) were performed in three measurement cycles, with a mix duration of 3 min, and measurement duration of 3 min. OCR and ECAR shown here represent the basal OCR and ECAR after injection of glucose, from the third measurement cycle. Immediately after the assay was run, cell numbers were imaged in each well using Hoechst 33342 DNA staining and quantified by a multi-imager (Cytation1, BioTek Instruments, Winooski, USA). Metabolic activities, i.e. OCR and ECAR, were normalised to the respective cell number and calculated per 104 cells.
Software, data analysis, and statistics
OCR and ECAR calculation from short-term measurements
After the application of a noise filter and the sensor calibration using the VisiSensVS software, regions of interest (ROI) were drawn around areas that corresponded to the cell layer, and pO2 concentrations over time exported to a data frame using the “multi-z” tool. For pH measurements, noise filter and sensor calibration were performed using the VisiSensVS software. The images and calibration files were then imported into the VS Live Plugin (Release 0.4.7) (PreSens, Regensburg, Germany). ROI were drawn and pH values over time exported to a data frame. pH and pO2 were plotted using R  and the R package ggplot2 . A linear regression was calculated for the linear range of the pH and O2 curves. The respective slopes of the regression lines represent the decrease of O2 concentration or pH over time (%O2/s and pH/s), thus the OCR and ECAR. To normalise for cell numbers, the slope was divided by the respective cell numbers calculated from fluorescent images of the cells’ nuclei and multiplied by a factor of 107 (for OCR) and 104 (for ECAR) to improve the way of representation of the data, as indicated in the respective figures. The normalised OCR and ECAR were plotted as dotplots using R and ggplot2.
Analysis of long-term measurements
After data acquisition, noise correction and sensor calibration were performed using the VisiSensVS software (PreSens, Regensburg, Germany). Image and calibration files were then loaded into the VDT 96-well Extractor programme (PreSens, Release 0.3.1.2). A one-point adjustment (“OPA”) was performed using an image of a plate with the same pO2 or pH in all wells to correct for uneven illumination of the plate. pO2 and pH values were then automatically exported to a CSV file. The mean and standard deviation were calculated and data were plotted as line graphs using R and ggplot2.
Quantification of cell numbers
Sixteen-bit greyscale images obtained by fluorescence microscopy were loaded into CellProfiler  and nuclei detected using the module “IdentifyPrimaryObjects” with advanced settings that were adjusted according to fluorescence intensity, nuclear size, and density of cells.
Image processing and figure preparation
All statistical analyses were performed using R and the package rstatix. To determine differences between different groups, one-way ANOVA tests with post hoc Tukey HSD were performed. Differences were considered to be statistically significant when p ≤ 0.05.
Fluorescence-based spatiotemporal measurement of O2 consumption and extracellular pH in HNSCC cells
To measure O2 consumption and extracellular acidification in cell cultures, we developed a device for the spatiotemporal pO2 and pH quantification in adherent cell cultures in a 6-well format (STO2p-Q). STO2p-Q enables a temporary reduction of the media volume surrounding the cells, creating a micro-respirator and subsequent measurement of the pO2 or pH in the cellular microenvironment by lowering sensor foil-equipped plungers onto the cells in each well simultaneously. The plungers consist of metal discs (weight 4 g), connected to strings that are routed through the lid of a 6-well plate and attached to a plastic lid that can be lowered and lifted through an attached string routed through a loop and then through the door of a cell culture incubator. This facilitates the operation of STO2p-Q in an incubator without disturbing the environmental conditions, i.e. temperature and CO2 concentration. Within the incubator, the measurement plate is placed on a frame that facilitates the excitation and detection of the fluorescence from below using the VisiSens TD camera readout device. LEDs facilitating excitation were attached to the frame (Fig. 1a).
To test if STO2p-Q enables accurate measurements of OCR and ECAR, effects of known disruptors of mitochondrial respiration and glycolysis, respectively, were measured in a HNSCC and ESCC cell line, respectively. FaDu and Kyse30 cells were seeded into one well of 2-well silicone inserts, inserted into a 6-well plate, and allowed to attach and form a confluent layer. After the addition of the inhibitors, pO2 and pH measurements were performed as described in the methods section. FCCP is an uncoupler of mitochondrial oxidative phosphorylation and is expected to increase cellular oxygen consumption, while oligomycin, an inhibitor of H+ ATP-synthase, inhibits mitochondrial respiration and thus is expected to decrease cellular oxygen consumption. In FaDu cells, addition of 1 μM FCCP increased the cells’ oxygen consumption compared to the DMSO control, while the addition of 1.5 μM oligomycin decreased the OCR significantly compared to the FCCP treatment (Fig. 1b). In Kyse30 cells, the effect of the inhibitors was more pronounced, with FCCP significantly increasing and oligomycin significantly decreasing the OCR compared to the DMSO control (Fig. 1c).
To manipulate the glycolysis rate and thus lactate production and subsequent extracellular acidification, 2-DG and oligomycin in combination with glucose were used. 2-DG is a glucose analogue that competes with glucose for phosphoglucoisomerase but cannot be metabolised, thus inhibiting glycolysis. Oligomycin blocks mitochondrial respiration and in the presence of glucose is expected to increase the ECAR. In FaDu cells the addition of 3 μM oligomycin and 10 mM glucose caused an increase in the ECAR, while 50 mM 2-DG resulted in a lower ECAR compared to the DMSO control, however not statistically significant (Fig. 1d). In Kyse30 cells, the measurements were more variable with the only statistically significant difference observed between cells treated with 2-DG and oligomycin and glucose, with oligomycin and glucose showing a significantly higher ECAR (Fig. 1e).
The results demonstrate that STO2p-Q in combination with the here established data analysis pipeline is able to distinguish differences in OCRs and ECARs between cells treated with different compounds.
Differences in ECAR and OCR between HNSCC cell lines and healthy cells
One of our aims was to develop a method to distinguish healthy tissue from cancerous tissue in patients with HNSCCs. As a first step, we explored potential differences in OCR and ECAR between HNSCC/ESCC cells and normal cells in vitro. In order to identify potential differences, three established HNSCC cell lines (FaDu, Cal33, Cal27), one ESCC cell line (Kyse30), and two normal cell lines (human foreskin fibroblasts, HFF; primary human nasal epithelial cells, HNEpC) were measured using STO2p-Q equipped with O2 or pH sensor foils as described. In Kyse30 cells, the measured OCR was significantly higher than in all other cell lines. FaDu cells displayed the second highest OCR, differing significantly from Cal33, Cal27 and the two normal cell lines, HFF and HNEpC, both of which showed the lowest OCR. The OCR of Cal33 and Cal27 did not differ significantly from HFF and HNEpC (Fig. 2a). To verify our results obtained by STO2p-Q, we performed OCR measurements in all cell lines via the established method of extracellular flux analysis using the Seahorse technology. The results obtained from the Seahorse measurements reflected the results obtained by our method in that Kyse30 displayed the highest basal OCR, followed by FaDu and Cal27, and showing significantly higher OCR than HFF and HNEpC cells. However, in the Seahorse measurements Cal33 showed a significantly lower OCR than HFF and HNEpC (Fig. 2b).
Using STO2p-Q, the pH measurements showed a significantly higher ECAR in HNEpCs compared to all other cell lines except FaDu, while HFFs exhibited the lowest ECAR. Differences between HFF and HNSCC cell lines were not statistically significant; however, a tendency for higher OCR, especially in Kyse30 and FaDu, compared to HFF was detectable (Fig. 2c). The results from the Seahorse measurement confirmed a low ECAR (in the presence of 5 mM glucose) in HFFs, which differed significantly from all other cell lines, while there was no significant difference between the HNSCC/ESCC cell lines (Fig. 2d), comparable to the results obtained with STO2p-Q.
Spatial discrimination of different cell lines based on O2 consumption
Detection of the spatial distribution of tumour cells and healthy cells is particularly relevant for clinical applications. To test the ability of STO2p-Q to identify areas containing tumour cells and to define the detection limit and resolution, different 2D-patterns of cells (“landscapes”) were created. O2 concentrations in the cellular landscapes were measured using the device equipped with O2 sensor foils. First, HNSCC/ESCC cells and normal cells were seeded into 4-well silicon inserts with a growth area of 0.35 cm2/well and with each well containing a different cell type (Fig. 3a). After removal of the silicon inserts, measurements were performed. The wells containing Kyse30 and FaDu cells could be identified after 2 min based on the O2 heat maps and after 12 min the outlines of the areas containing FaDu and Kyse30 cells were clearly detectable (Fig. 3a). Areas containing Kyse30 cells were characterised by the lowest O2 concentration, followed by FaDu, Cal27, Cal33, HNEpC, and HFF, reflecting the previous results from measurements of individual cell lines (see Figs. 2a and 3a).
Next, we were interested in defining an approximate detection limit for HNSCC/ESCC cell densities based on their OCR. Therefore, different cell numbers between 1000 and 8000 cells of FaDu, Kyse30, Cal22, and Cal27 were seeded into 4-well silicon μ-inserts with a growth area of 0.03 cm2/well. O2 concentrations were measured after removal of the silicon insert. Exact cell numbers were determined subsequently by fluorescence microscopy. After 10 min, the wells containing the highest cell numbers could be visualised in all cell lines, with the strongest signal, i.e., lowest O2 concentration, in Kyse30 cells. After 60 min, the two highest cell densities (i.e. from 7800 to 9400 cells and 3300 to 5800 cells, respectively) were detectable in all cell lines, and a weak signal was visible in most cell lines for the second lowest cell density (i.e. 1500 to 2400 cells) (Fig. 3b). This indicates that the detection limit for the HNSCC cells used in this study is at approximately 2000 cells growing in an area of 0.03 cm2, which translates to a cell density of approximately 66,000 cells/cm2.
In an approach to determine the minimum size of a cell colony resolvable by STO2p-Q, FaDu and Kyse30 cells were seeded into 6-well plates at low density and allowed to grow into colonies of different sizes. The cells in Kyse30 colonies were evenly spread out and forming a two-dimensional layer, while FaDu colonies were growing in dense and packed colonies where cells appeared to grow in multiple layers. Kye30 colonies were detectable after 10 min down to a size of 1 mm2 (Additional file 1, A). FaDu colonies could be detected down to a size of 0.6 mm2 within 2 min (Additional file 1, B).
To create a scenario that closely resembles an in vivo situation, HNSCC/ESCC cells and normal cells were seeded together in 2D-landscapes. As Kyse30 and FaDu cells showed the highest OCR in previous measurements, they were used in this experiment, together with HFF, which showed the lowest OCR (see Figs. 2a and 3a). The respective tumour cells stably expressed a cytosolic mCherry fluorescent protein that enabled us to identify the cells by fluorescence microscopy prior to the O2 measurements. After the measurements, the outline of tumour and normal cell areas, obtained from phase contrast and fluorescence micrographs, was overlaid with the O2 heat maps to mark the areas occupied by cells. Areas with Kyse30 cells within a layer of HFF were first detectable after 1 min of measurement. After 3 min, the areas containing the tumour cells were more clearly identifiable based on O2 heat maps. However, areas with an apparently lower cell density or smaller extensions of the tumour cell layer into the HFF layer were not represented in the O2 heat maps (Fig. 3c). In the landscapes with FaDu cells, areas containing tumour cells started being visible in the O2 heat maps after 3 min and were more clearly detectable after 8 min. Similar to the landscapes with Kyse30 cells, areas with lower cell densities and extensions were not resolved (Fig. 3d). The intensity of the signal for FaDu cells was lower than in Kyse30 cells and the time until detection of a signal was longer, which is in accordance with the lower OCR measured in FaDu cells compared to Kyse30 in previous experiments (Fig. 2a).
In an alternative approach, mCherry expressing FaDu or Kyse30 cells were seeded into 6-well plates at low density and allowed to form colonies. Subsequently, HFF cells were seeded into the wells and arranged themselves around the colonies to form a single layer. Areas containing Kyse30 colonies larger than 2 mm2 could be distinguished from HFF cells between 2 and 10 min of measurement (Fig. 3e), while FaDu colonies with a size between 0.5 and 0.8 mm2 were resolvable (Fig. 3f), in accordance with the results from the measurements of individual colonies. Alternatively, HFF cells were seeded first and FaDu and Kyse30 colonies allowed to grow on the layer of HFF cells. While Kyse30 cells formed colonies with low cell densities that could only be detected after approximately 30 to 50 min (Additional file 1, C), FaDu cells formed dense colonies that were detectable between 2 and 10 min (Additional file 1, D).
With these experiments, we demonstrated that it is feasible to spatially distinguish cell types based on their two-dimensional O2 consumption profiles using STO2p-Q and determined the approximate limit of resolution.
Spatial discrimination of different cell lines based on extracellular acidification
To test the ability of STO2p-Q to spatially distinguish different cell types based on extracellular acidification, we repeated the experiments described above, measuring pH instead of O2 concentrations. When seeded into 4-well silicon inserts, areas containing cells were first detectable in the pH heat map after 3 min. After 20 min, areas containing HNEpC cells exhibited the lowest pH (Fig. 4a). This is in accordance with previously measured ECARs for HNEpCs and HFFs, representing the cells with the highest and lowest ECAR, respectively (Fig. 2c).
Next, to determine the detection limit of the pH sensors for tumour cells, different cell numbers were seeded into smaller 4-well silicon inserts. After 10 min, wells containing the highest cell numbers (7900–12,000 cells) could be clearly distinguished in the pH heat map. After 60 min, areas containing the second highest cell numbers (3600–6100 cells) were also clearly identifiable. In Kyse30 and Cal33, a weak signal was identified in wells with the second lowest cell number (1700 and 2500 cells, respectively (Fig. 4b)). Based on these results, the detection limit for pH measurements is in the range of 2000 and 6000 cells per 0.03 cm2 or between 66,000 and 200,000 cells/cm2.
Individual Kyse30 colonies could be detected down to a size of 1.3 mm2 within 10 min (Additional file 2, A), and FaDu colonies down to a size of 0.5 mm2 (Additional file 2, B), similar to the resolution limit of pO2 measurements.
Similar to the experiments with the oxygen sensor foils described above, HFFs and either Kyse30 or FaDu cells were seeded to create cellular landscapes. Areas containing Kyse30 cells started exhibiting a lower pH compared to the surrounding areas after 2 min of measurement. After 20 min those areas were clearly identifiable, however, similar to the oxygen heat maps, smaller extensions of the tumour cell area were not resolved (Fig. 4c). The landscapes with FaDu cells and HFF behaved similarly, with a decrease in pH detectable in the areas containing FaDu cells after 3 min instead of 2 min in Kyse30 cells (Fig. 4d).
Kyse30 colonies over 1.4 mm2 within a layer of HFF cells were detectable between 5 and 15 min of pH measurements (Fig. 4e), while FaDu colonies were detectable at a minimum size of 1.2 mm2 (Fig. 4f). FaDu colonies growing on top of HFF cells could be visualised down to a size of 0.6 mm2 (Additional file 2, D), while Kyse30 colonies had a lower cell density and individual colonies were not clearly detectable based on pH measurements (Additional file 2, C).
These experiments showed that we could spatially distinguish different cell types based on their extracellular acidification using STO2p-Q and defined an approximate resolution limit.
Long-term OCR and ECAR measurements in HNSCC and normal cells
In the experiments described above, short-term measurements with STO2p-Q were used to quantify differences in the cell’s OCR and ECAR. However, we were also interested in investigating how cells consumed oxygen and acidified their environment in long-term culture conditions. For this purpose, 96-well plates with the bottom consisting of oxygen sensor foil material and pH sensor foil spots glued into the bottom of the wells, respectively, were used. The cells were directly seeded onto the sensor foils, which allows for monitoring of the direct cellular microenvironment. Plates were placed onto the frame also used for STO2p-Q, and measurements performed using the same camera system and software. Cells were allowed to grow for 90 h, during which the oxygen or pH were monitored every 60 min.
First, different numbers of FaDu cells were seeded and measurements were started 5 h later. The O2 concentration curves reflected the number of cells seeded, in that the higher the cell number, the faster the oxygen consumption (Fig. 5a).
Next, 40,000 cells of four HNSCC/ESCC cell lines (FaDu, Kyse30, Cal27, Cal33) and two normal cell lines (HFF, HNEpC) were seeded and measurements started after 5 h. Differences between cell lines became most pronounced after 25 h of measurement, with HFF cells showing the lowest oxygen consumption and FaDu and Kyse30 cells consuming oxygen faster than the other cell lines (Fig. 5b).
For pH measurements, different numbers of FaDu cells were seeded and measurements started after 1 h. The course of the pH values over time closely resembled the course of oxygen concentrations for the respective cell numbers, with high cell numbers showing a steep decrease in pH while curves for low cell numbers remained flatter (Fig. 5c).
When comparing different cell lines, pH curves of Kyse30 and HNEpC showed a similar course, representing the cell lines with the fastest acidification. FaDu cells exhibited the strongest acidification overall, while HFF showed the slowest decrease in pH (Fig. 5d).
This measurement setup allowed us to monitor the pH and oxygen concentration in the direct cellular environment over several days without interference. Furthermore, decreases in both O2 concentration and pH correlated with cell numbers and cell type.
While several different techniques exist to measure OCRs and ECARs of cells , they do not support easy spatial measurements. Electrode-based measurement devices rely on either moving the electrode between different localisations  or using multiple electrodes to gain spatial information, which limits both the spatial and temporal resolution. Fluorescence- or luminescence-based oxygen and pH measurement sensors are available in soluble forms or as immobilised sensors, as in the Seahorse system that is commonly used for the determination of OCR and ECAR [31,32,33,34,35]. While the Seahorse system is highly automated and thus convenient, it does not allow for spatially resolved OCR and ECAR measurements.
In this study, we describe the development and application of a newly developed device, which we have termed STO2p-Q referring to spatiotemporal O2 and pH quantification. We have used the STO2p-Q device to measure OCR and ECAR in HNSCC cell lines and normal cells, and found that some of the HNSCC/ESCC cell lines, especially FaDu and Kyse30 cells, consumed more oxygen than the two normal fibroblast and epithelial cell lines, which was confirmed by independent measurements with the Seahorse system. The Seahorse measurements also provided absolute OCR and ECAR for the different cell lines, which STO2p-Q is currently not suitable for as the media volume in the micro-respirator created by lowering the stamps is unknown. The relative differences in OCR enabled the spatial distinction of the two commonly used HNSCC/ESCC cell lines FaDu and Kyse30 from fibroblasts in a co-culture approach based on their oxygen consumption. Similarly, FaDu and Kyse30 cells could be spatially distinguished from HFFs based on their extracellular pH. This lays the foundation for further in vitro and in vivo applications, where a spatial resolution of O2 concentrations and pH in the cellular microenvironment is of importance.
In our experiments, tumour cell colonies down to a size of 0.5 mm2 could be resolved, depending on cell type and density of formed colonies. A limiting factor in gaining higher resolutions might be the rate at which O2 concentrations or pH can be restored in the direct cellular microenvironment through diffusion, thus preventing detection of O2 consumption and acidification by the sensor.
We are currently developing an automated, wireless, and programmable device that will allow for repeated short-term measurements over longer periods and for a more controlled cellular microenvironment. Diffusion could be minimised by decreasing the space between cell surface and sensor foil using stamps with a heavier weight. However, a balance between minimal diffusion and preventing damage to the cells would have to be warranted. By redesigning the stamps to have a small rim, potential horizontal diffusion at the sensor borders could be prevented. This could also facilitate the determination of the media volume in the micro-respirator, thus enabling absolute OCR and ECAR measurements and potentially improving variability between experiments. Furthermore, the resolution was limited by the camera setup used in this study. In future experiments, we will use detection cameras equipped with higher resolution objectives and thus attempt to further improve the detection limit.
The future device will also facilitate tracking changes in cell populations over time, as multiple sequences of measurements can be programmed and a sterile environment can be ensured. Combining STO2p-Q with a suitable excitation source and camera, a simultaneous detection of O2 or pH and fluorescent proteins expressed by certain cell populations is technically feasible. This would facilitate tracking of different cell populations over time and eliminate the requirement for removing the cells from the incubator to perform fluorescence microscopy between measurements.
Another application of the VisiSens technology we described in this study is long-term measurements in a 96-well format. However, these measurements do not allow for a direct determination of OCR and ECAR and comparisons of different cell lines, as cells continue to proliferate over the long measurement periods and would confound the data. However, long-term measurements are valuable for tracking changes within cell cultures in the direct cellular microenvironment over time without disturbances, which is especially important for cell lines sensitive to changes in pH or O2 concentrations.
In a clinical setting, future enhancement of the current technology could help identify residual cancerous tissue in patients in an intraoperative setting. Our in vitro measurements showed that the technology is capable of detecting differences in O2 consumption and extracellular pH between normal cells and some HNSCC/ESCC cell lines, which is a prerequisite for distinguishing healthy from cancerous tissue also in vivo. As local recurrences in HNSCC are very frequent, there is an urgent need for improving both early detection of cancer and for superior techniques to detect cancerous tissue during resection of the primary tumour. A variety of novel approaches has been discussed such as intraoperative molecular imaging , or specimen mapping by near-infrared fluorescence . Especially real time, intraoperative imaging of tumour cells is promising to obtain additional visual information and to avoid residual disease. A modified version of the VisiSens system has been used previously to investigate pO2 and pH gradients in chronic and acute wounds  and in postoperative wounds after radiation therapy . These studies have used transparent flexible sensor foils in combination with a handheld device for excitation of the sensor foils and acquisition of O2 concentrations and pH. In accessible areas of the head and neck, this setup could be used to assess cancerous and surrounding tissues before, during, and after surgery to identify the specific outlines of the tumour. Accordingly, we are currently preparing a study to investigate if cancerous tissue can be detected based on O2 consumption and pH in vivo. One of the main advantages of a setup using sensor foils in combination with a handheld device is that no treatment, e.g. labelling, of the tissue would be necessary and that the measurements can be performed in a timeframe of minutes . The applicability of sensor foils on uneven tissue and in the presence of liquids, e.g. saliva and blood, will have to be assessed in future experiments, as those could be potential obstacles to clear measurements. Similarly, the effect of food, medication, etc. on the variability of pH measurements in the oral cavity might represent a challenge.
Another clinical application, where pH and pO2 are of major importance, is radiation therapy as a common treatment in neck and head cancer. Despite benefits with respect to disease control, irradiation may show interference with wound repair processes. Sufficient oxygenation is mandatory for essential molecular healing processes like cell proliferation and protein synthesis . Chronic wounds in irradiated patients show hypoxic characteristics , qualifying pO2 as a major parameter in monitoring wound healing progresses after surgical interventions. Furthermore, damage of healthy skin through radiotherapy is reflected by an altered pH and pO2 compared to non-irradiated skin .
Therefore, pO2 and pH represent metabolic parameters of clinical and prognostic value and a device that allows for a standardised and easy measurement of these parameters would represent a considerable benefit for patients in numerous different settings.
In this study, we show the development of a device, the STO2p-Q, to spatiotemporally measure OCR and ECAR in cells. Using STO2p-Q, we were able to measure differences in OCR between HNSCC/ESCC cells and normal cells and to spatially resolve the different cell types based on both OCR and ECAR. This is, to the best of our knowledge, the first description of a device capable of spatiotemporally resolving OCR and ECAR at high resolution in cell cultures. The described method will be valuable for different in vitro applications in the investigation of cellular metabolism while also laying the foundation for the development of novel diagnostic tools for HNSCC patients.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request
Extracellular acidification rate(s)
Oesophageal squamous cell carcinoma(s)
Human foreskin fibroblast(s)
Human nasal epithelial cell(s)
Head and neck squamous cell carcinoma(s)
Spatiotemporal O2 and pH quantification
Sodium-hydrogen antiporter 1
Oxygen consumption rate(s)
- pHe :
- pHi :
- pO2 :
Partial oxygen pressure
Positive surgical margin(s)
Region(s) of interest
Vos T, Abajobir AA, Abate KH. Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet. 2017;390:1211–59.
Weiss BG, Ihler F, Anczykowski MZ, Bertlich M, Kitz J, Steiner W, et al. Transoral laser microsurgery for treatment of oropharyngeal cancer in 368 patients. Head Neck. 2019;41(9):3144–58. https://doi.org/10.1002/hed.25806.
Leemans CR, Snijders PJF, Brakenhoff RH. The molecular landscape of head and neck cancer. Nat Rev Cancer. 2018;18(5):269–82. https://doi.org/10.1038/nrc.2018.11.
Mes SW, Leemans CR, Brakenhoff RH. Applications of molecular diagnostics for personalized treatment of head and neck cancer: state of the art. Expert Rev Mol Diagn. 2016;16(2):205–21. https://doi.org/10.1586/14737159.2016.1126512.
Ganly I, Patel S, Shah J. Early stage squamous cell cancer of the oral tongue-clinicopathologic features affecting outcome. Cancer. 2012;118(1):101–11. https://doi.org/10.1002/cncr.26229.
Fakurnejad S, Krishnan G, van Keulen S, Nishio N, Birkeland AC, Baik FM, et al. Intraoperative molecular imaging for ex vivo assessment of peripheral margins in oral squamous cell carcinoma. Front Oncol. 2020;9:1476. https://doi.org/10.3389/fonc.2019.01476.
Chambers KJ, Kraft S, Emerick K. Evaluation of frozen section margins in high-risk cutaneous squamous cell carcinomas of the head and neck. Laryngoscope. 2015;125(3):636–9. https://doi.org/10.1002/lary.24945.
Vander Heiden M, Cantley L, Thompson C. The metabolic requirements of cell proliferation. Science. 2009;324(5930):1029–33. https://doi.org/10.1126/science.1160809.
Warburg O. Über den Stoffwechsel der Carcinomzelle. Naturwissenschaften. 1924;12(50):1131–7. https://doi.org/10.1007/BF01504608.
Warburg O. On the origin of cancer cells. Science. 1956;123(3191):309–14.
Counillon L, Bouret Y, Marchiq I, Pouysségur J. Na+/H+ antiporter (NHE1) and lactate/H+ symporters (MCTs) in pH homeostasis and cancer metabolism. Biochim Biophys Acta - Mol Cell Res. 2016;1863(10):2465–80. https://doi.org/10.1016/j.bbamcr.2016.02.018.
White KA, Grillo-Hill BK, Barber DL. Cancer cell behaviors mediated by dysregulated pH dynamics at a glance. J Cell Sci. 2017;130(4):663–9. https://doi.org/10.1242/jcs.195297.
Ashby BS. pH studies in human malignant tumours. Lancet. 1966;288(7458):312–5. https://doi.org/10.1016/S0140-6736(66)92598-0.
Engin K, Leeper DB, Cater JR, Thistlethwaite AJ, Tupchong L, Mcfarlane JD. Extracellular ph distribution in human tumours. Int J Hyperth. 1995;11(2):211–6. https://doi.org/10.3109/02656739509022457.
Schreml S, Meier RJ, Wolfbeis OS, Landthaler M, Szeimies R-M, Babilas P. 2D luminescence imaging of pH in vivo. Proc Natl Acad Sci. 2011;108(6):2432–7. https://doi.org/10.1073/pnas.1006945108.
Schreml S, Meier RJ, Kirschbaum M, Kong SC, Gehmert S, Felthaus O, et al. Luminescent dual sensors reveal extracellular pH-gradients and hypoxia on chronic wounds that disrupt epidermal repair. Theranostics. 2014;4(7):721–35. https://doi.org/10.7150/thno.9052.
Schreml S, Szeimies RM, Prantl L, Karrer S, Landthaler M, Babilas P. Oxygen in acute and chronic wound healing. Br J Dermatol. 2010;163(2):257–68. https://doi.org/10.1111/j.1365-2133.2010.09804.x.
Zhang J, Nuebel E, Wisidagama DRR, Setoguchi K, Hong JS, Van Horn CM, et al. Measuring energy metabolism in cultured cells, including human pluripotent stem cells and differentiated cells. Nat Protoc. 2012;7(6):1068–85. https://doi.org/10.1038/nprot.2012.048.
Peniche Silva CJ, Liebsch G, Meier RJ, Gutbrod MS, Balmayor ER, van Griensven M. A new non-invasive technique for measuring 3D-oxygen gradients in wells during mammalian cell culture. Front Bioeng Biotechnol. 2020;8. https://doi.org/10.3389/fbioe.2020.00595.
Wolff P, Heimann L, Liebsch G, Meier RJ, Gutbrod M, van Griensven M, et al. Oxygen-distribution within 3-D collagen I hydrogels for bone tissue engineering. Mater Sci Eng C. 2019;95:422–7. https://doi.org/10.1016/j.msec.2018.02.015.
Rangan SRS. A new human cell line (FaDu) from a hypopharyngeal carcinoma. Cancer. 1972;29(1):117–21. https://doi.org/10.1002/1097-0142(197201)29:1<117::AID-CNCR2820290119>3.0.CO;2-R.
Shimada Y, Imamura M. Prognostic significance of cell culture in carcinoma of the oesophagus. Br J Surg. 1993;80(5):605–7. https://doi.org/10.1002/bjs.1800800518.
Gioanni J, Fischel JL, Lambert JC, Demard F, Mazeau C, Zanghellini E, et al. Two new human tumor cell lines derived from squamous cell carcinomas of the tongue: establishment, characterization and response to cytotoxic treatment. Eur J Cancer Clin Oncol. 1988;24(9):1445–55. https://doi.org/10.1016/0277-5379(88)90335-5.
Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nature Methods. 2012;9(7):676–82. https://doi.org/10.1038/nmeth.2019.
R Core Team (2020). R: A language and environment for statistical computing. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2020.
Wickham H. ggplot2. ggplot2; 2009.
Carpenter AE, Jones TR, Lamprecht MR, Clarke C, Kang IH, Friman O, et al. CellProfiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biol. 2006;7(10):R100. https://doi.org/10.1186/gb-2006-7-10-r100.
No Title. p. Inkscape. Available from: http://www.inkscape.org/. Accessed 4 Apr 2019.
Dmitriev RI, Papkovsky DB. Optical probes and techniques for O 2 measurement in live cells and tissue. Cell Mol Life Sci. 2012;69(12):2025–39. https://doi.org/10.1007/s00018-011-0914-0.
Kagawa Y, Matsuura K, Shimizu T, Tsuneda S. Direct measurement of local dissolved oxygen concentration spatial profiles in a cell culture environment. Biotechnol Bioeng. 2015;112(6):1263–74. https://doi.org/10.1002/bit.25531.
Divakaruni AS, Paradyse A, Ferrick DA, Murphy AN, Jastroch M. Analysis and interpretation of microplate-based oxygen consumption and pH data. In: Methods in enzymology; 2014.
Pike Winer LS, Wu M. Rapid analysis of glycolytic and oxidative substrate flux of cancer cells in a microplate. Plos One. 2014;9(10):e109916. https://doi.org/10.1371/journal.pone.0109916.
Wei C, Heitmeier M, Hruz PW, Shanmugam M. Evaluating the efficacy of GLUT inhibitors using a seahorse extracellular flux analyzer. In: Methods in molecular biology; 2018.
Zhang J, Zhang Q. Using seahorse machine to measure OCR and ECAR in cancer cells. In: Methods in molecular biology; 2019.
Martin SD, McGee SL. A systematic flux analysis approach to identify metabolic vulnerabilities in human breast cancer cell lines. Cancer Metab. 2019;7:12. https://doi.org/10.1186/s40170-019-0207-x.
Gao RW, Teraphongphom NT, van den Berg NS, Martin BA, Oberhelman NJ, Divi V, et al. Determination of tumor margins with surgical specimen mapping using near-infrared fluorescence. Cancer Res. 2018
Auerswald S, Schreml S, Meier R, Blancke Soares A, Niyazi M, Marschner S, et al. Wound monitoring of pH and oxygen in patients after radiation therapy. Radiat Oncol. 2019;14(1):199. https://doi.org/10.1186/s13014-019-1413-y.
The authors thank Prof. Markus Meissner (Faculty of Veterinary Medicine, LMU) for providing HFF cells.
The project was funded by the Bayrisches Staatsministerium für Wirtschaft, Energie und Technologie, with the funding programme BayBIO (grant number 1803-003). Open Access funding enabled and organized by Projekt DEAL.
Ethics approval and consent to participate
Consent for publication
Co-authors Robert Meier and Gregor Liebsch are employees of PreSens Precision Sensing GmbH, Regensburg, Germany.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
pO2 measurements in cell colonies, A, B Individual Kyse30 (A) or FaDu (B) cells grown into different size colonies. On the left, brightfield micrographs are shown with the colony size indicated below. On the right, heatmaps represent pO2 after 10, 40 and 60 minutes (A) or 2, 30 and 60 minutes (B). Scale bars: 1 mm C,D Individual Kyse30 (C) or FaDu (D) cells expressing mCherry grown into different size colonies on a layer of HFF cells. On the left, micrographs with the brightfield and mCherry signal overlaid are shown. Kyse30/FaDu colonies are outlined in black or white and were overlaid with the heatmaps representing the pO2 after 30 and 50 minutes (C) or 2 and 10 minutes (D).
pH measurements in cell colonies. A, B Individual Kyse30 (A) or FaDu (B) cells grown into different size colonies. On the left, brightfield micrographs are shown with the colony size indicated below. On the right, heatmaps represent pH after 10, 30 and 60 minutes. Scale bars: 1 mm C,D Individual Kyse30 (C) or FaDu (D) cells expressing mCherry grown into different size colonies on a layer of HFF cells. On the left, micrographs with the brightfield and mCherry signal overlaid are shown. Kyse30/FaDu colonies are outlined in black or white and were overlaid with the heatmaps representing the pH after 15 and 40 minutes (C) or 5 and 15 minutes (D).
About this article
Cite this article
Blancke Soares, A., Meier, R., Liebsch, G. et al. High-resolution spatiotemporal pHe and pO2 imaging in head and neck and oesophageal carcinoma cells. Cancer Metab 9, 21 (2021). https://doi.org/10.1186/s40170-021-00257-6