Next Article in Journal
CD123 and More: How to Target the Cell Surface of Blastic Plasmacytoid Dendritic Cell Neoplasm
Next Article in Special Issue
Targeting the Tumor Microenvironment through mTOR Inhibition and Chemotherapy as Induction Therapy for Locally Advanced Head and Neck Squamous Cell Carcinoma: The CAPRA Study
Previous Article in Journal
New Developments in the Pathogenesis, Therapeutic Targeting, and Treatment of Pediatric Medulloblastoma
Previous Article in Special Issue
The Hidden Link of Exosomes to Head and Neck Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

mRNA Subtype of Cancer-Associated Fibroblasts Significantly Affects Key Characteristics of Head and Neck Cancer Cells

1
Department of Pathological Physiology, Faculty of Medicine, Masaryk University, Kamenice 5, 62500 Brno, Czech Republic
2
Department of Physiology, Faculty of Medicine, Masaryk University, Kamenice 5, 62500 Brno, Czech Republic
3
Department of Chemistry and Biochemistry, Mendel University in Brno, Zemedelska 1, 61300 Brno, Czech Republic
4
Max Planck Institute for the Science of Light, Staudtstraße 2, 91058 Erlangen, Germany
5
Department of Otorhinolaryngology and Head and Neck Surgery, First Faculty of Medicine, University Hospital Motol, Charles University, V Uvalu 84, 15006 Prague, Czech Republic
6
Department of Pathology and Molecular Medicine, 2nd Faculty of Medicine, University Hospital Motol, Charles University, V Uvalu 84, 15006 Prague, Czech Republic
7
BIOCEV, First Faculty of Medicine, Charles University, Prumyslova 595, 25250 Vestec, Czech Republic
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(9), 2286; https://doi.org/10.3390/cancers14092286
Submission received: 25 February 2022 / Revised: 27 April 2022 / Accepted: 27 April 2022 / Published: 3 May 2022
(This article belongs to the Special Issue Tumor Microenvironment of Head and Neck Cancer)

Abstract

:

Simple Summary

The interplay between cancer-associated fibroblasts (CAFs) and head and neck cancer cells and its effect on the metabolic state, mechanobiology, and aggressiveness of cancer cells is described. In vitro coculture techniques were used. Not all CAFs support tumours–the right kind of cooperation between cancer cells and CAFs is needed for tumour growth and progression, and only specific mRNA subtypes of CAFs can support the growth of primary cancer cells or metastases carrying different mutational statuses.

Abstract

Head and neck squamous cell carcinomas (HNSCC) belong among severe and highly complex malignant diseases showing a high level of heterogeneity and consequently also a variance in therapeutic response, regardless of clinical stage. Our study implies that the progression of HNSCC may be supported by cancer-associated fibroblasts (CAFs) in the tumour microenvironment (TME) and the heterogeneity of this disease may lie in the level of cooperation between CAFs and epithelial cancer cells, as communication between CAFs and epithelial cancer cells seems to be a key factor for the sustained growth of the tumour mass. In this study, we investigated how CAFs derived from tumours of different mRNA subtypes influence the proliferation of cancer cells and their metabolic and biomechanical reprogramming. We also investigated the clinicopathological significance of the expression of these metabolism-related genes in tissue samples of HNSCC patients to identify a possible gene signature typical for HNSCC progression. We found that the right kind of cooperation between cancer cells and CAFs is needed for tumour growth and progression, and only specific mRNA subtypes can support the growth of primary cancer cells or metastases. Specifically, during coculture, cancer cell colony supporting effect and effect of CAFs on cell stiffness of cancer cells are driven by the mRNA subtype of the tumour from which the CAFs are derived. The degree of colony-forming support is reflected in cancer cell glycolysis levels and lactate shuttle-related transporters.

1. Introduction

Head and neck squamous cell carcinoma (HNSCC) is the sixth most prevalent cancer worldwide [1] and represents a heterogeneous group of tumours originating from mucosal epithelial cells that line the nasal and oral cavity, oropharynx, hypopharynx, and larynx. In recent years, it has been shown that cells within the tumour stroma drastically modulate the behaviour of tumour cells [2]. The predominant cell populations within the tumour microenvironment (TME) are cancer-associated fibroblasts (CAFs). In HNSCC, the late stages of carcinomas frequently contain about 80% of CAFs [3].
CAFs were originally presumed to represent a homogeneous group of stromal cells whose main role was to provide and maintain a favourable TME to promote the growth, survival, and metastatic spread of tumour cells. However, CAFs display a great heterogeneity regarding their actual tumour-promoting or tumour-inhibiting properties. The wide range and heterogeneity of functions attributed to CAFs raise the question of whether this cell population executes all these functions or whether specific subpopulations of CAFs perform specific tasks, such as supporting proliferation or promoting metastatic potential. This can also implicate possible overlaps of different CAF subpopulations or a possible switch between different functional states. It is crucial to point out that solid tumours display, not only histological and morphological heterogeneity, but also intertumoral, metabolic heterogeneity and metabolic plasticity. This characteristic is not restricted only to cancer cells, but it can also be found in surrounding stromal cells [4,5].
It has been observed that cancer cells and CAFs can establish a metabolic symbiosis in which cancer cells corrupt normal fibroblasts, which in turn produce metabolites for cancer cells to supply their metabolic needs. In addition to glucose, cancer cells can also utilize lactate as an alternative energy source via lactate shuttle, which is regulated by the conversion of lactate into pyruvate, as well as by the transport of lactate across the cancer cell plasma membrane. Interestingly, aerobic glycolysis can also take place in CAFs rather than in epithelial cells, which then, in turn, supplies cancer cells with lactate for OXPHOS. This lactate exchange between cancer cells and CAFs is a process known as the reverse Warburg effect [6]. Nevertheless, metabolic symbiosis may also occur between oxidative CAFs and glycolytic tumour cells, which has also been observed in HNSCC [5].
Apart from metabolic symbiosis, CAFs are known to support cancer cell migration [7], invasiveness, and metastasis formation [8]. Cell mechanic phenotypes largely impact the cell’s ability to migrate and invade due to changes in the cytoskeleton. Accordingly, the epithelial–mesenchymal transition (EMT) is accompanied by significant changes in cell stiffness [9]. Processes enabling migration are energetically demanding, and these demands further increase with increasing cell stiffness [10].
In our study, the classification to the mRNA cluster was based on Lawrence et al. [11] and modified accordingly: Basal phenotype was classified by high TP63 and EGFR and low SOX2 expression. The mesenchymal phenotype was classified by high VIM, DES, and low TP63 expression, the classical phenotype by high expression of SOX2, NFE2L2, and positive smoking history, the atypical phenotype was characterised by high SOX2 expression and p16 positivity. The mean of standardised log2 expression/smoking/p16 status values was calculated, and the patient was attributed to the category with the highest mean value (based on Lawrence et al. [11]). In this study, we investigated how CAFs derived from tumour tissues of different mRNA subtypes influence the proliferation, metabolic, and biomechanical reprogramming of cancer cells. We also investigated the clinicopathological significance of the expression of these metabolism-related genes in tissue samples of HNSCC patients to identify a possible gene signature typical for HNSCC progression. We assume that the right kind of cooperation between cancer cells and CAFs is needed for tumour growth and progression.

2. Material and Methods

2.1. Tumour Samples Collection

The study was approved by the ethical committee of Motol University Hospital, Prague, Czech Republic. All surgical tissue samples were obtained from HNSCC patients after they signed the informed consent documents. Patients were completely clinically examined, and tumour samples were taken from verified HNSCC under general anaesthesia; the inclusion criteria were as follows: histologically confirmed squamous cell carcinoma and no previous oncologic treatment; each tissue specimen underwent pathology quality control. Hematoxylin and eosin-stained sections from each sample were subjected to independent pathology review to confirm squamous cell carcinoma histological consistence. Tumour samples with approximately >60% tumour nuclei and <20% necrosis were considered sufficient for the following analyses. The tissue material harvested during surgery was divided into two pieces, the first was stored in an RNAlater (Ambion, Austin, TX, USA) for later RNA analyses. The second portion of the sample was placed into the cultivation medium (MEM) with the addition of 1% trypsin and 1% antibiotic-antimycotic solution (Santa Cruz Biotechnology, Dallas, TX, USA). The material was maintained at cold temperatures for 16–24 h to enable trypsin diffusion. Primary cell lines were subsequently prepared.

2.2. Ex-Vivo Cell Cultures

The tissue material was first rinsed with 50% ethanol (Merck, Darmstadt, Germany), followed by phosphate buffer (Invitrogen, Waltham, MA, USA). The necrotic and adipose tissues were then removed from the tissue using a sterile scalpel, and the tumour tissue was cut into 3 mm pieces. Next, samples in the trypsin media were kept at 37 °C to activate the trypsin. Subsequently, the tumour pieces were transferred to a sterile tube with phosphate buffer and centrifuged at 4 °C, 2700 rpm for 7 min. After removing the supernatant, the tumour pieces were placed into the culture medium (MEM) with added 10% FBS (Biochrom, Holliston, MA, USA), 1% antibiotic-antimycotic solution, 10 μg/mL gentamicin sulfate, and 10 μg/mL ciprofloxacin (Santa Cruz Biotechnology, Dallas, TX, USA) to prevent contamination. After 7 days the culture medium was replaced with MEM medium with 10% FBS, supplemented with antibiotics (penicillin 100 U/mL and streptomycin 0.1 mg/mL) (Biochrom, Holliston, MA, USA) and 0.4 μg/mL hydrocortisone (Merck, Darmstadt, Germany). Since we were primarily interested in the effect of TME-derived fibroblasts (CAFs), selection using population overgrowth was performed. The passages of CAFs ranged from 1 to 5.

2.3. Cell Lines

Two human head and neck squamous cell carcinoma cell lines were used in this study. FaDu (RRID:CVCL_1218) is an epithelial adherent cell line derived from hypopharyngeal carcinoma. Detroit 562 (RRID:CVCL_1171) is an epithelial adherent cell line derived from the pleural effusion of metastatic pharyngeal carcinoma. Authenticated cell lines were purchased from ATCC (Manassas, VA, USA) within the last three years.

2.4. Cell Lines Cultivation

FaDu and Detroit 562 cells were cultured in MEM medium with 10% FBS, supplemented with antibiotics (penicillin 100 U/mL and streptomycin 0.1 mg/mL) and 0.4 μg/mL hydrocortisone. Cells were maintained at 37 °C in a humidified incubator with 5% CO2 (Sanyo, Moriguchi, Osaka, Japan). The passages of all cell lines performed in our lab ranged from 5 to 10. All experiments were performed with mycoplasma-free cells. Mycoplasma was detected by PCR (primers MYCO_A GGCGAATGGGTGAGTAACACG and MYCO_B CGGATAACGCTTGCGACCTATG).

2.5. Conditioned Media Preparation

Cell-conditioned media were collected from 90% confluent FaDu, Detroit 562 cell lines, and the CAFs were cultivated for 24 h. The collected media were frozen at −20 °C until further use.

2.6. Colony-Forming Assay

The FaDu and Detroit 562 cell lines were cocultured with the CAFs using the Corning® Transwell® cell culture inserts (Merck, Darmstadt, Germany). The FaDu and Detroit 562 cells were seeded in 6-well plates, while the CAFs were seeded in inserts and left to adhere overnight. The seeding density was 100 cancer cells in the bottom chamber and 5000 CAFs in the insert (1:50 ratio) to form colonies in 3 weeks. After cultivation, colonies were fixed with cold methanol, visualised by trypan blue and analysed. For analysis, regions of interest were chosen by registration of each image to the reference image (with the manually labelled area of interest). Next, the colonies were segmented by thresholding of the blue component of the image, where a single fixed threshold was used. Finally, the fraction of areas covered by colonies was computed.

2.7. RNA Isolation and Reverse Transcription

The FaDu and Detroit 562 cell lines were cocultured with the CAFs using Transwell® cell culture inserts. The FaDu and Detroit 562 cells were seeded in 6-well plates at a density of 2 × 104, while the CAFs were seeded at a density of 5 × 103 in inserts and left to adhere overnight. The inserts were then inserted into the wells and cocultures were cultured for 7 days. After 7 days, the medium was collected, and cells were harvested for RNA isolation. TriPure Isolation Reagent (Roche, Bazel, Switzerland) was used for RNA isolation from tissue and harvested cells. The isolated RNA was used for cDNA synthesis. The RNA (1000 ng) was transcribed using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Waltham, MA, USA), which was used according to the manufacturer’s instructions. cDNA (20 μL) prepared from total RNA was diluted with RNase-free water to 100 μL, and the amount of 5 μL was directly analysed.

2.8. Quantitative Real-Time Polymerase Chain Reaction

qRT-PCR was performed using TaqMan gene expression assays with the LightCycler®480 II System (Roche, Bazel, Switzerland) and the amplified cDNA was analysed by the comparative Ct method using PSMB2 as a housekeeping gene. Primer and probe sets for PSMB2 (Hs01009704_m1), SOX2 (Hs01053049_s1), EGFR (Hs01076078_m1), TP63 (Hs00978340_m1), VIM (Hs00958111_m1), DES (Hs00157258_m1), NFE2L2 (Hs00975961_g1), ACTA2 (Hs05005341_m1), MCT1 (Hs00161826_m1), MCT4 (Hs00358829_m1), CD147 (Hs00174305_m1), CAV1 (Hs00971716_m1), were selected from the TaqMan gene expression assays (Thermo Fisher Scientific, Waltham, MA, USA). qRT-PCR was performed under the following amplification conditions: total volume 20 μL, initial incubation at 50 °C/2 min followed by denaturation at 95 °C/10 min, then 45 cycles at 95 °C/15 s and 60 °C/1 min.

2.9. Lactate Assay

L-lactate levels in coculture media were quantitatively measured with a colourimetric L-lactate assay kit (Abcam, Cambridge, UK) according to the manufacturer’s protocol. Briefly, coculture media and reaction mix, consisting of assay buffer, substrate mix, and enzyme mix were added into a 96-well plate. After incubation for 30 min at room temperature, the colour signal was measured at 450 nm.

2.10. Live-Cell Metabolic Assay

The assay was performed according to the manufacturer’s instructions. Briefly, to examine the effect of the HNSCC CAF-derived conditioned media, the FaDu and Detroit 562 cells were seeded at the density of 1 × 104 cells per well in a Seahorse 8-well plate (Agilent Technologies, Santa Clara, CA, USA) and incubated for 24 h at 37 °C in the CAF-derived conditioned media. A day before the analysis, Seahorse cartridge (Agilent Technologies) was hydrated and incubated at 37 °C without CO2. On the day of the experiment, the Seahorse XFp analyzer (Agilent Technologies, Santa Clara, CA, USA) was calibrated using hydrated cartilage. The culture medium was changed to Seahorse XF RPMI medium, pH 7.4 (Agilent Technologies, Santa Clara, CA, USA). After incubating cells for 1 h at 37 °C without CO2, the Seahorse assay was performed by measuring the extracellular acidification rates (ECAR) as a proxy for glycolytic readouts and the oxygen consumption rates (OCR) as a proxy for oxidative phosphorylation readouts. The data were acquired and analysed using Wave Controller 2.4 software (Agilent Technologies, Santa Clara, CA, USA).

2.11. Immunohistochemistry

Representative FFPE tumour tissue blocks were selected for immunohistochemical analysis. Immunohistochemistry was performed on 2 µm sections using specific primary antibodies against p16 (BD Pharmingen™, BD Biosciences, Franklin Lakes, NJ, USA, clone G175-405, dilution 1:100), MCT1 (Abcam, Cambridge, UK, clone SLC16A1, ab90582, 1:100), MCT4 (Merck, Darmstadt, Germany, clone SLC16A3, HPA021451, 1:300), CD147 (Abcam, Cambridge, UK, clone ab666, 1:100), IRF3 (Abcam, Cambridge, UK, clone ab68481, 1:500) and α-SMA (Abcam, Cambridge, UK, clone ab5694, 1:300). Antigen demasking was performed at 99 °C, by pH6 or pH9. As a secondary antibody, a Mouse/Rabbit PolyDetector HRP Label (BioSB, Goleta, CA, USA) was used. All the procedures were performed according to added datasheets.
We assessed the expressions as negative or positive. In positive cases, we graded the intensity into three grades (1+ weak staining; 2+ moderate; 3+ strong positivity). At staining of p16, we evaluated the intensity, percentage of positive cells, and location of expression (nuclear and/or cytoplasmic) in tumour cells. The expressions of the remaining antibodies were assessed in tumour cells, and in peritumoral tissue (inflammatory cells, stromal elements), evaluating intensity, location of expression (cytoplasm, membrane), and coarse estimating the number of positive cells.

2.12. Atomic Force Microscopy

For the analysis of cell stiffness of FaDu and Detroit 562 cells, the cells were cultured in the CAF-derived conditioned media for 24 h. Then, the bioAFM microscope JPK NanoWizard 3 (JPK BioAFM, Berlin, Germany) was placed on the inverted optical microscope Olympus IX-81 (Olympus, Tokyo, Japan) and equipped with a fluorescence and confocal module, thus allowing a combined experiment (AFM-optical combined images). The maximal scanning range of the AFM microscope in the X-Y-Z range was 100-100-15 µm. The typical approach/retract settings were identical with a 15 μm extend/retract length, a setpoint value of 1 nN, a pixel rate of 2048 Hz, and a speed of 30 µm/s. The system operated under closed-loop control. After reaching the selected contact force, the cantilever was retracted. The retraction length of 15 μm was sufficient to overcome any adhesion between the tip and the sample, and to make sure that the cantilever had been completely retracted from the sample surface. A force-distance (FD) curve was recorded at each point of the cantilever approach/retract movement. AFM measurements were obtained at 37 °C with force measurements recorded at a pulling speed of 30 µm/s (extension time 0.5 s).
The Young’s modulus (E) was calculated by fitting the Hertzian–Sneddon model on the FD curves measured as force maps (64 × 64 points) of the region containing either a single cell or multiple cells. JPK data evaluation software was used for the batch processing of measured data. The adjustment of the cantilever position above the sample was carried out under the microscope by controlling the position of the AFM-head by a motorised stage, equipped with a Petri dish heater (JPK) allowing a precise positioning of the sample together with a constant elevated temperature of the sample for the whole period of the experiment. Soft uncoated AFM probes HYDRA-2R-100N (Applied NanoStructures, Mountain View, CA, USA), i.e., silicon nitride cantilevers with silicon tips, are used for stiffness studies because they are maximally gentle to living cells (not causing mechanical stimulation). Moreover, as compared with coated cantilevers, these probes are very stable under elevated temperatures in liquids—thus allowing long-time measurements without nonspecific changes in the measured signal.

2.13. Quantitative Phase Imaging

The migration capacity of the FaDu cells before and after coculture with conditioned media was assessed using Tescan multimodal holographic microscope Q-PHASE. The Q-PHASE is based on the original concept of the coherence-controlled holographic microscope. Measurements were carried out under conditions optimised in our previous study [12]. Briefly, cells were cultivated in flow chambers μ-Slide I Lauer Family (Ibidi GmbH, Gräfelfing, Germany). To image enough cells in one field of view, the Nikon Plan 10/0.30 lens was chosen. Holograms were captured by a CCD camera (XIMEA, Marianka, Slovakia). Complete image reconstruction and image processing were performed in Q-PHASE control software.

2.14. Flow Cytometry

To verify the purity of the CAF cell culture, the cells were stained for the mesenchymal surface marker CD90 and analysed by flow cytometry. The cells were prepared as a single-cell suspension for FACS staining. For CD90 staining, cells were stained with APC-CD90 antibody (Biolegend, San Diego, CA, USA) for 30 min at 4 °C. CD90+ and CD90− cells were identified based on isotype gating. The stained cells were acquired for analysis on a FACS Aria II flow cytometer (BD Biosciences, Franklin Lakes, NJ, USA). Staining for dead cells was performed using SYTOXblue staining. Flow cytometry data were analysed with FlowJo software (BD Biosciences, Franklin Lakes, NJ, USA).

2.15. Real-Time Deformability Cytometry

Real-time deformability cytometry (RT-DC) was performed according to Rosendahl et al. [13] using an AcCellerator (Zellmechanik, Dresden, Germany) at 37 °C. After 48 h of indirect coculture of cancer cells with the CAF-derived media, the cells were trypsinised and kept on a roller in the culture medium in an incubator for 30 min to reduce their surface/area ratios. Subsequently, cells were transferred to the measurement buffer CellCarrier (1xPBS with 0.5% methylcellulose and adjusted to a specific viscosity) and measured using RT-DC in 30 μm wide channel microfluidic chips at a total flow of 0.240 μL/s. The acquisition was performed using ShapeIn 2.0.8 software (Zellmechanik, Dresden, Germany). The Young’s modulus was calculated from deformation and cell area using ShapeOut software [14], after filtering with the following setting: cell area 70–2000 μm2 (removal of cell debris/multicellular clusters), porosity 1.00–1.05 (removal of damaged cells). Measurements were performed in duplicates per treatment with a minimum of 1000 filtered cells per measurement.

2.16. Statistical Analysis

Pearson correlations were performed to calculate correlation coefficients between the colony-forming assay and gene expression. A one-sample t-test was used to compare gene expression of cocultures vs. controls. R 4.0.2 [15] was used for analysis with the following packages: ggplot2 [16], gplots 3.0.4 [17], heatmap.plus 1.3 [18], Corrplot 0.84 [19], ggpubr 0.4.0 [20]. Survival analysis was performed with Kaplan–Meier curves and using Cox proportional hazards model with Stepwise model selection (backward and forward) with the packages MASS4 [21], survival 3.2–11 [22] and survminer 0.4.9 [23]. Unless noted otherwise, p-level < 0.05 was considered significant.
Classification to mRNA cluster was based on [11] and modified accordingly: Basal phenotype was classified by high TP63 and EGFR and low SOX2 expression. The mesenchymal phenotype was classified by high VIM, DES, and low TP63 expression, the classical phenotype was characterised by high expression of SOX2 and NFE2L2, and positive smoking history, atypical by high SOX2 expression and p16 positivity. The mean of standardised log2 expression/smoking/p16 status values was calculated and the patient was attributed to the category with the highest mean value.

3. Results

3.1. Clinical Characterization of Patients and Tumours Used for CAF Preparation

A total of 55 patients with head and neck squamous cell cancers were used in the study. The cohort was defined as follows: squamous cell carcinoma histology, curative therapeutic intention using surgery +/− radio or radio/chemotherapy, no distant metastatic dissemination. Tumours with the following locations were included: oropharynx (N = 24), larynx (N = 16), hypopharynx (N = 6), tongue (N = 5), the floor of the mouth (N = 3), and oral cavity (N = 1). The p16 status was analysed in all subjects. TNM8 system was used. The description of patients is shown in Table 1, Tables S1 and S6. Of those patients, a randomly selected subset of 18 patient tumour tissues was used to prepare cancer-associated fibroblast cultures (Tables S2 and S5). Its lineage specificity was confirmed, and these CAF cultures were consequently used for coculture experiments with existing cancer cell lines.

3.2. Lineage Specificity of Patient-Derived CAFs and Model Cancer Cells

Based on the low doubling time of fibroblasts and their rapid growth, population overgrowth was performed. Population overgrowth of patient-derived cell cultures by CAFs was successfully reached after 2 passages (Figure 1a). Our selection method was confirmed by flow cytometry of three randomly selected early passage CAFs. The CD90 (a CAFs marker [24]) positivity of these selected cultures ranged between 79% and 99% (Figure 1b and Figure S1), indicating the successful establishment of the CD90+ CAF population. The CAFs cell cultures showed typical spindle-shaped cell morphology (Figure 1c).
FaDu is an epithelial adherent cell line derived from a primary squamous cell hypopharyngeal tumour that was chosen as the wt-PI3K model. Detroit 562 is an epithelial adherent cell line derived from pleural effusion metastasis of pharyngeal carcinoma with the PI3K activating mutation. Both cell lines are HPV negative. FaDu cells are heterozygous for TP53 missense mutation (p.R248L) and Detroit 562 cells are homozygous for missense mutation p.R175H [25]. First, the CD90/CD44 status of these model cells was verified by flow cytometry (see Figure 1d), as CD44 has been reported to be involved in tumour growth and metastasis in head and neck squamous cell cancer [26]. FaDu showed predominantly CD44/CD90 negativity, while 95.5% of Detroit 562 cells showed CD44+/CD90− status. Furthermore, the mRNA expression profile of FaDu and Detroit cells was performed. Neither FaDu nor Detroit 562 cells showed the mesenchymal phenotype characterised by high expression of EMT markers, such as vimentin (VIM) and desmin (DES). Both cell lines were characterised by high expression of EGFR and NFE2L2. FaDu cells showed low levels of the SOX2 expression relative to the expression of TP63, in Detroit 562 cells the expression of SOX2 was higher relative to the expression of TP63. FaDu had a higher expression of MCT1 and CD147 compared to Detroit 562 cells (see Figure 1e).

3.3. CAF-Cancer Cell Metabolic Symbiosis Is mRNA-Subtype-Specific

Coculture of cancer cells with patient-derived CAFs demonstrated that only some CAFs support the colony-forming of cancer cells and that this phenomenon is cancer cell type-specific (Figure 2a,b). Accordingly, based on a difference to non-cocultured cells, it is possible to create groups of cancer cell colony-supporting and -non-supporting CAFs separately for FaDu and Detroit 562 cells (Figure S2). Interestingly, the degree of cancer cell-colony forming support by CAFs is associated with the tissue-of-origin mRNA subtype (Figure 2c). For the metastatic Detroit 562 cell line, the least supportive CAFs were those derived from the classical (CL) mRNA subtype of HNSCC tumour mass. On the contrary, these CL CAFs were the most supportive of the colony-forming capacity of the primary FaDu cell line. ME CAFs showed the least support for colony-formation of FaDu cells (Figure 2c). It was also obvious, that some mRNA subtypes of CAFs such as AT and BA CAFs can support the growth of both types of tumour cells (primary tumour, metastasis), while CL CAFs supported only the colony-forming capacity of FaDu cells (Figure 2c, Table S3).
Interestingly, a similar mRNA-subtype-related association was also observed for the glycolysis levels (Figure 2d,e). The coculture of cancer cells with CL, AT, and ME CAFs significantly increased glycolysis in FaDu cells. Cancer cell glycolysis levels were in a positive correlation with proliferation (proliferation-promoting) in FaDu cells but were in negative correlation with proliferation in Detroit 562 cells (Figure 2f,g, Table S4). Gene expression of MCT1, MCT4, and CD147 was, therefore, performed to verify whether lactate transport is responsible for this proliferation promotion. These lactate shuttle-related genes were downregulated in cancer cells by coculture with supporting CAFs (Figure 2f and Figure S3). However, the expression of lactate transporters in Detroit 562 derived from metastatic cells was significantly less influenced by coculture. The increased area of colonies was shown to be in a negative correlation with MCTs and CD147 expression only in FaDu primary tumour cells, this effect was not observable in metastatic Detroit 562 cells Figure 2g. Furthermore, the cell division rate assessed by QPI was in negative correlation with MCT4 expression in FaDu cells cocultured with CAFs (Figure 2h,i).
In summary, these data indicate the cancer cell colony-supporting effect is driven by the mRNA subtype of the tumour from which the CAFs are derived, and this degree of colony-forming support is reflected in cancer cell glycolysis level and lactate shuttle-related transporters.

3.4. Cell Stiffness of FaDu Cells Is Associated with Mitochondrial ATP Production and MCT1 Expression

Together with cancer cell migration, changes in the cell mechanical phenotype are important in vitro indicators of cell aggressiveness [27]. As demonstrated previously, distinct cellular mechanic properties support the malignant aggressive potential, including cell migration [28]. Measurement of cell stiffness was therefore performed after coculture. Coculture of tumour cells with CAF-derived conditioned media significantly affects the mechanical properties of cancer cells, as demonstrated by real-time deformability cytometry (Figure S4) and atomic force microscopy (Figure S5). Similarly to the effect of CAFs on colony-forming capacity and lactate shuttle, the coculture effect on cell stiffness is mRNA-subtype specific (Figure 3a,b). Coculture of FaDu cells with the BA subtype of CAFs led to an increase in stiffness of FaDu cells (see Figure 3c). Results obtained from both methods were in positive correlation (see Figure 3d). Furthermore, the value of cell stiffness was negatively associated with mitochondrial ATP production and MCT1 expression in FaDu cells (see Figure 3e–g). Interestingly, CAFs from ME tumours cause the most pronounced decrease of cancer cell stiffness, accompanied by minimal changes in lactate shuttle, and conversely, CAFs from CL tumours cause the most pronounced stiffening accompanied by the highest decrease of MCT1.

3.5. Cancer Cells Can Manipulate Their CAFs

The effect of coculture on the lactate transporters’ gene expression was bidirectional as the coculture of patient-derived CAF cell lines with FaDu and Detroit 562 cells also had a significant effect on the gene expression of MCT4 and CD147 in CAFs (see Figure 4a,b). The CAFs showed an increase in lactate production and MCT4 expression due to coculture with cancer cells. An increase of CD147 expression in the CAFs was specific for the coculture with Detroit 562 cells. Furthermore, MCT1 expression in tumour tissues was significantly higher than the expression in the CAFs (see Figure 4c) but MCT4 expression was not changed between tumour tissue and the CAFs. Therefore, high MCT1 expression in tumour tissue can be attributed to cancer cells.
In conclusion, the data indicate a bi-directional interplay in lactate shuttling between CAFs and cancer cells (Figure 4d), where changes in lactate transporters are accompanied by an increase of extracellular lactate when cocultured.

3.6. Patients with Basal mRNA Subtype of HNSCC and Overexpression of Lactate Transport-Associated Genes Have a Poor Prognosis

Regarding implications for a clinical picture, tissue expression of lactate transport-related genes was analysed directly in tumorous tissues of 55 patients used in the study prior to CAF isolation. The clustering based on MCT1, MCT4, CD147, CAV1, and ACTA2 expression resulted in the formation of two clusters (Figure 5a). These clusters demonstrated an association with the mRNA subtype of tumours, with the basal subtype being predominantly located in the second cluster. Interestingly, these two clusters were not associated with any clinically relevant head and neck parameters, including tissue site, p16 status, staging, and smoking. On the other hand, an association with survival was observed: a trend, although insignificant, was observed using a Kaplan–Meier survival analysis (Figure 5b), and, more importantly, using a stepwise Cox model with tissue expression cluster, pN, staging, p16, grade, and smoking as independent predictors, it was demonstrated that “lactate gene cluster” together with tumour stage pN are independent predictors of overall survival (Figure 5c). The cluster of patients with worse overall survival was characterised by high expression of lactate-related transporters as well as high ACTA2 and CAV1 expression in tumorous tissues. Interestingly, when the expression of these genes was dichotomized by median (high/low) none of them were demonstrated to be an independent overall survival predictor, indicating that rather a specific interplay between all those genes is characteristic of tumour aggressiveness and hence patient outcome. Accordingly, CAFs supporting the expression of MCT4 in Detroit 562 cells came significantly more often from patients who later developed recurrence (Table S3).
Regarding tissue MCT1 and MCT4 expression patterns, a partial agreement was observed with mRNA expression data; for MCT1, positivity was detected in the membrane of isolated cancer cells rather than as a homogeneous expression in tumour tissue. In accordance with mRNA MCT1 expression and immunohistochemical analysis, positivities were observed in tumours classified as basal (3/4 samples were positive) as well as atypical (3/4) and classical (3/4). On the other hand, tumours classified as mesenchymal were negative (1/4) for MCT1 expression (Figure 5d). The stromal expression was negative for MCT1. MCT4 followed a similar trend in tumorous tissues and demonstrated isolated positivity in the stroma (rather inflammatory elements) of some tumorous tissues (Figure 5d).

4. Discussion

Cancer-associated fibroblasts (CAFs) play many important roles in the progression of a tumour. They produce proteins of the extracellular matrix and secrete inflammatory and growth factors promoting cancer progression, cancer cell proliferation, and therapy resistance [2]. Many studies indicate that CAFs are necessary for HNSCC progression [5,29]. However, some works indicate that CAFs may also suppress tumour progression in some circumstances [30]. Both the tumour-promoting and tumour-suppressing effects of CAFs were supported by our coculture experiments.
CAFs are not a homogenous cellular type but display significant molecular and metabolic heterogeneity and high levels of metabolic plasticity. In this study, we investigated the interaction between different mRNA subtypes of CAFs and HNSCC cells and its impact on the key characteristics of these cancer cells. The mRNA subtypes of tumour tissues were determined according to Lawrence et al. [11]. Based on our coculture experiments, it was shown that CAFs derived from different mRNA subtypes of HNSCC tumour tissue had different abilities to promote the colony-forming capacity or modulate the mechanical properties of HNSCC cell lines of dissimilar origin, indicating the presence of divergent transcriptional programmes during disease progression. Consequently, fibroblasts in metastases can be activated in a different way than CAFs inhabiting primary tumours [31]. For the metastatic Detroit 562 cell line, the least supportive CAFs were those derived from the CL mRNA subtype of HNSCC tumour mass characteristic by αSMA-low, elevated levels of oxidative stress response genes, including NFE2L2 and SOX2, and a history of heavy smoking. On the other hand, the most supportive were CAFs derived from the ME mRNA subtype of HNSCC tumour mass, characterised by high expression of EMT markers αSMA, vimentin, and desmin and low gene expression of the TP63. The level of support is cancer cell-type-specific; CL CAFs were the most supportive for the colony-forming capacity of the primary FaDu cells, and ME CAFs showed the least support for colony formation of FaDu cells (Figure 2c). It was also obvious that CAFs derived from the AT (αSMA-low, SOX2-high, p16+) or BA (αSMA- and EGFR-high and low SOX2/TP63 ratio) mRNA subtypes could support the growth of both types of tumour cells (primary-tumour ones as well as metastatic-derived ones), while CL and ME CAFs were more specific. The right type of cooperation between CAFs and cancer cells may be needed for tumour progression. For example, the proliferation of primary cancer cells not yet prepared for EMT can be suppressed by the ME subtype of CAFs, as was shown by our results. The CAFs inducing high expression of MCT1 in FaDu cells did not support the proliferation of FaDu cells. MCT1 expression was associated with epithelial to mesenchymal transition (EMT) [32]. In epithelial cells, autocrine ERBB3 activation sustains PI3K signalling. After EMT, downregulation of ERBB3 interrupts autocrine signalling to PI3K, which leads to reduced proliferation after EMT [33]. The increased area of colonies was shown to be in negative correlation with MCT1 expression only in primary tumour cells, this was not true for the metastatic Detroit 562 cells carrying the activation mutation of PI3K. FaDu did not carry PI3K mutations [25] and therefore the colony-forming ability of FaDu cells can be diminished by the EMT process supported by ME CAFs. Furthermore, EMT softens head and neck cancer cells to facilitate migration in 3D environments [34], making changes in the cell mechanical phenotype, an important in vitro indicator of cell aggressiveness [27]. Indeed, coculture of FaDu cells with BA subtype of CAFs increased the cell stiffness of FaDu cells, which implicates EMT inhibition by these fibroblasts.
The differences between FaDu and Detroit 562 cells might not only be related to the primary versus metastatic status of the cells but also (partially) due to the different mutations carried by these cells. In addition to the aforementioned presence of the activating PI3K mutation in Detroit 562 cells, the Detroit 562 cell line has a p53 R175H mutation [35]. Though R175H p53 loses the function of the wild-type p53, p53-R175H aggregates can not only induce coaggregation of wild-type p53 to cause the loss of function, but also induce coaggregation of p63 and p73 to cause a gain of function. p63 and p73 mutations are rare in tumours, but their function can be disrupted by coaggregation with mutant p53, leading to an increase in the oncogenic potential of the affected cells [36]. p63 supports aerobic respiration and p63-depleted cells had a significantly reduced oxygen consumption rate [37]. Depletion of p73 caused altered lysine metabolism and glycolysis, as well as an elevated pentose phosphate pathway and abnormal lipid accumulation. p73 has been shown to regulate metabolic enzymes, such as cytochrome c oxidase subunit IV isoform 1, glucose 6-phosphate dehydrogenase, and glutaminase-2 [38]. In the FaDu cell line, a p53 mutation occurs at codon 248 (p.R248L) [39]. A codon 248 p53 mutation preserves some tumour suppressor function [40]. Homozygous deletion of SMAD4 was also observed in FaDu cells [41,42]. SMAD4 gene depletion can induce autophagy. The protein level of SMAD4 was inversely correlated with autophagy in orthotopic tumour tissue samples [43]. Furthermore, primary tumour cells (FaDu) contained 63.8% cells with CD44-negative status and only 36.1% cells with CD44+ status, 95.5% of metastatic Detroit 562 cells showed CD44+ status. CD44 triggers EMT in many cancers [44,45]. Therefore, EMT stimulating ME subtype of CAFs can probably cooperate better with CD44-positive metastatic Detroit 562 cells. Data indicates that the changes in cancer cell mechanics are linked also with their metabolic reprogramming. In the literature, evidence exists on the level of interplay between stromal stiffness and cancer cell OXPHOS [46] or glycolysis [47]. Nevertheless, the crosstalk between the cell mechanics and metabolic reprogramming is far from being clearly understood. However, coculture experiments performed in this study provide evidence that cancer cell stiffness is in negative correlation with their MCT1 expression and OXPHOS, suggesting that lactate shuttle and CAF-mediated HNSCC metabolic reprogramming are also associated with the support of aggressive potential via cell stiffness modulation [28]. Metabolic reprogramming is most likely used to sustain active processes in the actin cortex of cancer cells [48]. Similarly to colony-forming capacity, cell softening is an mRNA-subtype specific with the most profound softening after ME subtype coculture.
Interestingly, the coculture of cancer cells with CL CAFs significantly increased glycolysis in FaDu cells (Figure 2d and Figure 4d). Glycolysis was rather proliferation-promoting in FaDu cells, but was in negative correlation with proliferation in Detroit 562 cells, Figure 2g. Enhanced glycolysis can provide many important advantages helpful for cancer cells such as the possibility of rapid ATP synthesis, higher flux into biosynthetic pathways, and immune cell evasion [49]. Nevertheless, metastatic cancer cells without functioning oxidative phosphorylation (OXPHOS) show reduced proliferation, low migration and decreased metastatic potential. Mitochondrial ATP synthesis is not required to overcome these obstacles. On the other hand, the respiratory chain and dihydroorotate dehydrogenase (DHODH)-dependent pyrimidine biosynthesis is crucial [50,51].
Consequently, metastatic Detroit 562 cells can prefer higher OXPHOS and CAFs forcing them to glycolysis can be counterproductive for their progression. Furthermore, lactate promotes glutamine uptake and metabolism in oxidative cancer cells and these cancer cells can use lactate as fuel supporting their growth and division [52,53]. All subtypes of CAFs cocultured with cancer cells showed an increase in lactate production and MCT4 expression (see Figure 4b) suggesting metabolic manipulation by cancer cells and the possible establishment of metabolic symbiosis. It was observed that fibroblasts surrounding malignant cells have a low Cav-1 expression, high MCT4 expression and enhanced aerobic glycolysis, which mirrors the simultaneous increase of mitochondrial activity in the adjacent epithelial cancer cells [54]. An increase of the CD147 expression was specific for the coculture of CAFs with Detroit 562 cells. CD147 (extracellular matrix metalloproteinase inducer, EMMPRIN), a transmembrane glycoprotein and the chaperone of MCT1 and MCT4, has been studied as a potential prognostic marker of HNSCC. As its name suggests, CD147 has been thought to act mainly as an inducer of MMPs; however, its role in facilitating the cell membrane localization and functionality of MCTs implicates a dominant role in the regulation of aerobic glycolysis. It has been reported that CD147 plays a role in the transformation of normal fibroblasts into CAFs through the induction of α-SMA expression [55].
Regarding implications for the clinical picture, a group of patients with overexpression of lactate transport-associated genes had poor overall survival rates. This group contained most of the patients with basal mRNA subtype of HNSCC. CAFs supporting the expression of MCT4 in Detroit 562 cells came significantly more often from patients who later developed a recurrence. The BA subtype is characterised by high expression of EGFR, which can signalize CD44-positivity as CD44 upregulates the expression of EGFR, leading to the activation of the PI3K/Akt signalling pathway [45]. The BA subtype was αSMA-high which suggests the presence of myofibroblastic CAFs. Across various cancer types, myofibroblastic CAFs were associated with ECM production [30,56]. In HNSCC, the presence of myofibroblastic CAFs upregulated the expression of galectin-1 [57]. Galectin-1 in ECM modifies the effect of ECM on cells by influencing cell adhesion, cell migration through ECM, and cell death [58]. Targeting galectin-1 production in CAFs inhibited HNSCC metastasis [59]. Furthermore, the BA subtype of CAFs was able to support the growth of both primary tumour cells (FaDu) and metastatic cells (Detroit 562).

5. Conclusions

Cancer-associated fibroblasts play many important roles in the progression of HNSCC. Different cellular origins and factors produced by tumour cells shape the phenotype of CAFs and contribute to their heterogeneity and the presence of distinct CAF subtypes. The identity and proportion of these subtypes could differ across normal, premalignant, and malignant states, mutation statuses, and also between metastatic and primary tumours. Based on our results, we assume that the right kind of cooperation between cancer cells and CAFs is needed and premature support of some steps of cancerogenesis, such as EMT, can suppress cancer cells in some circumstances. A better insight into the CAF subpopulations emerging during cancer progression is crucial for successful therapeutic targeting of symbiosis between CAFs and cancer cells. As various CAF subpopulations have distinct roles in tumour progression, targeting them individually may lead to better outcomes. A deep characterization of CAF subtypes will therefore be pivotal for the design of effective combinatorial targeting of CAFs and cancer cells.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14092286/s1, Table S1: Clinical characteristics of the patient cohort; Table S2: Clinical characteristics of a subset of the patient cohort used for CAF generation for coculture experiments; Table S3: Cancer cells cocultured with HNSCC-derived cancer-associated fibroblasts, descriptive statistics and results of ANOVA; Table S4: Cancer cells cocultured with HNSCC-derived cancer-associated fibroblasts, results of Pearson correlation; Table S5: Summary of clinico-pathological characteristics of patients from which CAF cultures were established. All patients were indicated for surgery with a curative intent. The table is supplemented by the mRNA subtype; Table S6: Raw data used for analysis; Figure S1 Lineage specificity of cancer-associated fibroblasts, CD90/SYTOXblue status determined by flow cytometry; Figure S2: Colony area of FaDu and Detroit 562 cells cocultured in Transwell with cancer-associated fibroblasts; Figure S3: MCT1, MCT4, CD147, and CAV1 gene expression in FaDu and Detroit 562 cells cocultured with cancer-associated fibroblasts; Figure S4: FaDu cancer cell stiffness is affected by coculture with head and neck-derived cancer-associated fibroblasts; Figure S5. FaDu and Detroit 562 cancer cell stiffness is affected by coculture with head and neck-derived cancer-associated fibroblasts.

Author Contributions

Conceptualization, M.R., D.K. and M.M.; Data curation, B.P., H.H.P., J.B. (Jan Balvan), J.N., J.G. (Jaromir Gumulec), B.Č. and D.K.; Formal analysis, B.P., H.H.P., M.R., M.K., J.G. (Jochen Guck), M.G., J.P. and M.M.; Funding acquisition, M.M.; Investigation, H.H.P., J.N. and M.G.; Methodology, B.P., J.B. (Jan Balvan), J.N., T.V., B.Č., M.K., J.G. (Jochen Guck) and M.G.; Project administration, M.M.; Resources, J.P., J.B. (Jan Betka) and M.M.; Software, T.V. and J.G. (Jaromir Gumulec); Supervision, M.R., J.P., J.B. (Jan Betka) and M.M.; Validation, T.V., J.G. (Jaromir Gumulec) and J.B. (Jan Betka); Visualization, T.V. and J.G. (Jaromir Gumulec); Writing—original draft, B.P. and M.R.; Writing—review & editing, J.G. (Jaromir Gumulec), J.P., J.B. (Jan Betka) and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Grant Agency of the Czech Republic (GACR—18-03978S), by the Ministry of Health of the Czech Republic (NV18-08-00229), by Technology Agency of the Czech Republic (FW03010186), by the funds from Specific University Research Grant, as provided by the Ministry of Education, Youth and Sports of the Czech Republic in the year 2022 (MUNI/A/1391/2021 and MUNI/A/1133/2021) and by Center for Tumor Ecology—Research of the Cancer Microenvironment Supporting Cancer Growth and Spread (reg. No. CZ.02.1.01/0.0/0.0/16_019/0000785). We acknowledge the core facility CELLIM of CEITEC supported by the Czech-BioImaging large RI project (LM2018129 funded by MEYS CR) for their support in obtaining scientific data presented in this paper. The short-term scientific mission of J. Gumulec was supported by the COST project Action COMULIS CA17121.

Institutional Review Board Statement

The study was approved by the Ethical committee of Motol University Hospital, Prague, Czech Republic, 2412/85, 21 June 2017.

Informed Consent Statement

All surgical tissue samples were obtained from HNSCC patients after they signed the written informed consent.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ACTA2actin alpha 2
AFMatomic force microscopy
ATatypical mRNA subtype
BAbasal mRNA subtype
CAFcancer-associated fibroblast
CAVcaveolin
CD147cluster of differentiation 147
CLclassical mRNA subtype
DHODHdihydroorotate dehydrogenase
ECARextracellular acidification rate
EMTepithelial-mesenchymal transition
FFPEformalin-fixed, paraffin-embedded
HNSCChead and neck squamous cell carcinomas
HPVhuman papilloma virus
MCTmonocarboxylate transporter
MEmesenchymal mRNA subtype
OCRoxygen consumption rate
OXPHOSoxidative phosphorylation
RT-DCreal-time deformability cytometry
TMEtumour microenvironment
VIMvimentin

References

  1. Stewart, B.W.; Wild, C.P. World Cancer Report 2014; WHO Press: Geneva, Switzerland, 2014. [Google Scholar]
  2. Peltanova, B.; Raudenska, M.; Masarik, M. Effect of tumor microenvironment on pathogenesis of the head and neck squamous cell carcinoma: A systematic review. Mol. Cancer 2019, 18, 63. [Google Scholar] [CrossRef] [PubMed]
  3. Ansems, M.; Span, P.N. The tumor microenvironment and radiotherapy response; a central role for cancer-associated fibroblasts. Clin. Transl. Radiat. Oncol. 2020, 22, 90–97. [Google Scholar] [CrossRef] [PubMed]
  4. Gouirand, V.; Guillaumond, F.; Vasseur, S. Influence of the Tumor Microenvironment on Cancer Cells Metabolic Reprogramming. Front. Oncol. 2018, 8, 117. [Google Scholar] [CrossRef] [PubMed]
  5. Kumar, D.; New, J.; Vishwakarma, V.; Joshi, R.; Enders, J.; Lin, F.; Dasari, S.; Gutierrez, W.R.; Leef, G.; Ponnurangam, S.; et al. Cancer-Associated Fibroblasts Drive Glycolysis in a Targetable Signaling Loop Implicated in Head and Neck Squamous Cell Carcinoma Progression. Cancer Res. 2018, 78, 3769–3782. [Google Scholar] [CrossRef] [Green Version]
  6. Bonuccelli, G.; Tsirigos, A.; Whitaker-Menezes, D.; Pavlides, S.; Pestell, R.G.; Chiavarina, B.; Frank, P.G.; Flomenberg, N.; Howell, A.; Martinez-Outschoorn, U.E.; et al. Ketones and lactate “fuel” tumor growth and metastasis. Cell Cycle 2010, 9, 3506–3514. [Google Scholar] [CrossRef]
  7. Erdogan, B.; Ao, M.; White, L.M.; Means, A.L.; Brewer, B.M.; Yang, L.; Washington, M.K.; Shi, C.; Franco, O.E.; Weaver, A.M.; et al. Cancer-associated fibroblasts promote directional cancer cell migration by aligning fibronectin. J. Cell Biol. 2017, 216, 3799–3816. [Google Scholar] [CrossRef] [Green Version]
  8. Wen, S.; Hou, Y.; Fu, L.; Xi, L.; Yang, D.; Zhao, M.; Qin, Y.; Sun, K.; Teng, Y.; Liu, M. Cancer-associated fibroblast (CAF)-derived IL32 promotes breast cancer cell invasion and metastasis via integrin β3–p38 MAPK signalling. Cancer Lett. 2018, 442, 320–332. [Google Scholar] [CrossRef]
  9. Mierke, C.T. Mechanical Cues Affect Migration and Invasion of Cells From Three Different Directions. Front. Cell Dev. Biol. 2020, 8, 946. [Google Scholar] [CrossRef]
  10. Zanotelli, M.R.; Rahman-Zaman, A.; VanderBurgh, J.A.; Taufalele, P.V.; Jain, A.; Erickson, D.; Bordeleau, F.; Reinhart-King, C.A. Energetic costs regulated by cell mechanics and confinement are predictive of migration path during decision-making. Nat. Commun. 2019, 10, 4185. [Google Scholar] [CrossRef] [Green Version]
  11. Lawrence, M.S.; Sougnez, C.; Lichtenstein, L.; Cibulskis, K.; Lander, E.; Gabriel, S.B.; Getz, G.; Ally, A.; Balasundaram, M.; Birol, I.; et al. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature 2015, 517, 576–582. [Google Scholar]
  12. Vicar, T.; Raudenska, M.; Gumulec, J.; Balvan, J. The Quantitative-Phase Dynamics of Apoptosis and Lytic Cell Death. Sci. Rep. 2020, 10, 1566. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Rosendahl, P.; Plak, K.; Jacobi, A.; Kraeter, M.; Toepfner, N.; Otto, O.; Herold, C.; Winzi, M.; Herbig, M.; Ge, Y. Real-time fluorescence and deformability cytometry. Nat. Methods 2018, 15, 355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Müller, P.; O’Connell, E.; Schloögel, M. Shape-Out Version 2.6.4: Analysis Software for Real-Time Deformability Cytometry [Software]; Max Planck Institute: Dresden, Germany, 2021; Volume 2021, Available online: https://github.com/ZELLMECHANIK-DRESDEN/ShapeOut2 (accessed on 28 April 2020).
  15. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  16. Wickham, H. Ggplot2: Elegant Graphics for Data Analysised; Springer: New York, NY, USA, 2016. [Google Scholar]
  17. Warnes, G.R.; Bolker, B.; Bonebakker, L.; Gentleman, R.; Huber, W.; Liaw, A.; Lumley, T.; Maechler, M.; Magnusson, A.; Moeller, S.; et al. Gplots: Various R Programming Tools for Plotting Data; ScienceOpen: Berlin, Germany, 2020. [Google Scholar]
  18. Day, A. Heatmap.Plus: Heatmap with More Sensible Behavior. 2012. Available online: https://rdrr.io/cran/heatmap.plus/ (accessed on 25 March 2022).
  19. Wei, T.; Simko, V. R Package “Corrplot”: Visualization of a Correlation Matrix. 2017. Available online: https://scirp.org/reference/referencespapers.aspx?referenceid=3067218 (accessed on 25 March 2022).
  20. Kassambara, A. Ggpubr: ‘ggplot2’ Based Publication Ready Plots. 2020. Available online: https://github.com/kassambara/ggpubr (accessed on 25 March 2022).
  21. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S, 4th ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
  22. Therneau, T.M.; Grambsch, P.M. Modeling Survival Data: Extending the Cox Modeled; Springer: New York, NY, USA, 2000. [Google Scholar]
  23. Kassambara, A.; Kosinski, M.; Biecek, P. Survminer: Drawing Survival Curves Using ‘ggplot2’. R Package Version 0.4.9. 2021, Volume 2021. Available online: https://rpkgs.datanovia.com/survminer/ (accessed on 25 March 2022).
  24. Schliekelman, M.J.; Creighton, C.J.; Baird, B.N.; Chen, Y.; Banerjee, P.; Bota-Rabassedas, N.; Ahn, Y.H.; Roybal, J.D.; Chen, F.; Zhang, Y.; et al. Thy-1(+) Cancer-associated Fibroblasts Adversely Impact Lung Cancer Prognosis. Sci. Rep. 2017, 7, 6478. [Google Scholar] [CrossRef]
  25. Wirtz, E.D.; Hoshino, D.; Maldonado, A.T.; Tyson, D.R.; Weaver, A.M. Response of head and neck squamous cell carcinoma cells carrying PIK3CA mutations to selected targeted therapies. JAMA Otolaryngol. Head Neck Surg. 2015, 141, 543–549. [Google Scholar] [CrossRef] [Green Version]
  26. Chen, J.; Zhou, J.; Lu, J.; Xiong, H.; Shi, X.; Gong, L. Significance of CD44 expression in head and neck cancer: A systemic review and meta-analysis. BMC Cancer 2014, 14, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Mierke, C.T. The Role of the Optical Stretcher Is Crucial in the Investigation of Cell Mechanics Regulating Cell Adhesion and Motility. Front. Cell Dev. Biol. 2019, 7, 184. [Google Scholar] [CrossRef] [Green Version]
  28. Guck, J.; Schinkinger, S.; Lincoln, B.; Wottawah, F.; Ebert, S.; Romeyke, M.; Lenz, D.; Erickson, H.M.; Ananthakrishnan, R.; Mitchell, D.; et al. Optical deformability as an inherent cell marker for testing malignant transformation and metastatic competence. Biophys. J. 2005, 88, 3689–3698. [Google Scholar] [CrossRef] [Green Version]
  29. Magan, M.; Wiechec, E.; Roberg, K. CAFs affect the proliferation and treatment response of head and neck cancer spheroids during co-culturing in a unique in vitro model. Cancer Cell Int. 2020, 20, 599. [Google Scholar] [CrossRef]
  30. Biffi, G.; Tuveson, D.A. Diversity and Biology of Cancer-Associated Fibroblasts. Physiol. Rev. 2021, 101, 147–176. [Google Scholar] [CrossRef]
  31. Nielsen, S.R.; Quaranta, V.; Linford, A.; Emeagi, P.; Rainer, C.; Santos, A.; Ireland, L.; Sakai, T.; Sakai, K.; Kim, Y.S.; et al. Macrophage-secreted granulin supports pancreatic cancer metastasis by inducing liver fibrosis. Nat. Cell Biol. 2016, 18, 549–560. [Google Scholar] [CrossRef] [Green Version]
  32. Zhang, G.; Zhang, Y.; Dong, D.; Wang, F.; Ma, X.; Guan, F.; Sun, L. MCT1 regulates aggressive and metabolic phenotypes in bladder cancer. J. Cancer 2018, 9, 2492–2501. [Google Scholar] [CrossRef] [PubMed]
  33. Salt, M.B.; Bandyopadhyay, S.; McCormick, F. Epithelial-to-mesenchymal transition rewires the molecular path to PI3K-dependent proliferation. Cancer Discov. 2014, 4, 186–199. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Chen, Y.-Q.; Lan, H.-Y.; Wu, Y.-C.; Yang, W.-H.; Chiou, A.; Yang, M.-H. Epithelial-mesenchymal transition softens head and neck cancer cells to facilitate migration in 3D environments. J. Cell. Mol. Med. 2018, 22, 3837–3846. [Google Scholar] [CrossRef] [PubMed]
  35. De Smet, F.; Saiz Rubio, M.; Hompes, D.; Naus, E.; De Baets, G.; Langenberg, T.; Hipp, M.S.; Houben, B.; Claes, F.; Charbonneau, S.; et al. Nuclear inclusion bodies of mutant and wild-type p53 in cancer: A hallmark of p53 inactivation and proteostasis remodelling by p53 aggregation. J. Pathol. 2017, 242, 24–38. [Google Scholar] [CrossRef]
  36. Xu, J.; Reumers, J.; Couceiro, J.R.; De Smet, F.; Gallardo, R.; Rudyak, S.; Cornelis, A.; Rozenski, J.; Zwolinska, A.; Marine, J.-C.; et al. Gain of function of mutant p53 by coaggregation with multiple tumor suppressors. Nat. Chem. Biol. 2011, 7, 285–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Viticchiè, G.; Agostini, M.; Lena, A.M.; Mancini, M.; Zhou, H.; Zolla, L.; Dinsdale, D.; Saintigny, G.; Melino, G.; Candi, E. p63 supports aerobic respiration through hexokinase II. Proc. Natl. Acad. Sci. USA 2015, 112, 11577–11582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. He, Z.; Agostini, M.; Liu, H.; Melino, G.; Simon, H.U. P73 regulates basal and starvation-induced liver metabolism in vivo. Oncotarget 2015, 6, 33178–33190. [Google Scholar] [CrossRef]
  39. Kim, M.S.; Li, S.L.; Bertolami, C.N.; Cherrick, H.M.; Park, N.H. State of p53, Rb and DCC tumor suppressor genes in human oral cancer cell lines. Anticancer Res. 1993, 13, 1405–1413. [Google Scholar]
  40. Mukhopadhyay, T.; Roth, J.A. A codon 248 p53 mutation retains tumor suppressor function as shown by enhancement of tumor growth by antisense p53. Cancer Res. 1993, 53, 4362–4366. [Google Scholar]
  41. Qiu, W.; Schönleben, F.; Li, X.; Su, G.H. Disruption of transforming growth factor beta-Smad signaling pathway in head and neck squamous cell carcinoma as evidenced by mutations of SMAD2 and SMAD4. Cancer Lett. 2007, 245, 163–170. [Google Scholar] [CrossRef] [Green Version]
  42. Lin, L.-H.; Chang, K.-W.; Cheng, H.-W.; Liu, C.-J. SMAD4 Somatic Mutations in Head and Neck Carcinoma Are Associated with Tumor Progression. Front. Oncol. 2019, 9, 1379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wang, F.; Xia, X.; Yang, C.; Shen, J.; Mai, J.; Kim, H.-C.; Kirui, D.; Kang, Y.; Fleming, J.B.; Koay, E.J.; et al. SMAD4 Gene Mutation Renders Pancreatic Cancer Resistance to Radiotherapy through Promotion of Autophagy. Clin. Cancer Res. 2018, 24, 3176–3185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Xu, H.; Tian, Y.; Yuan, X.; Wu, H.; Liu, Q.; Pestell, R.G.; Wu, K. The role of CD44 in epithelial-mesenchymal transition and cancer development. OncoTargets Ther. 2015, 8, 3783–3792. [Google Scholar]
  45. Cho, S.H.; Park, Y.S.; Kim, H.J.; Kim, C.H.; Lim, S.W.; Huh, J.W.; Lee, J.H.; Kim, H.R. CD44 enhances the epithelial-mesenchymal transition in association with colon cancer invasion. Int. J. Oncol. 2012, 41, 211–218. [Google Scholar]
  46. Mah, E.J.; Lefebvre, A.E.; McGahey, G.E.; Yee, A.F.; Digman, M.A. Collagen density modulates triple-negative breast cancer cell metabolism through adhesion-mediated contractility. Sci. Rep. 2018, 8, 17094. [Google Scholar] [CrossRef] [Green Version]
  47. Park, J.S.; Burckhardt, C.J.; Lazcano, R.; Solis, L.M.; Isogai, T.; Li, L.; Chen, C.S.; Gao, B.; Minna, J.D.; Bachoo, R.; et al. Mechanical regulation of glycolysis via cytoskeleton architecture. Nature 2020, 578, 621–626. [Google Scholar] [CrossRef]
  48. Yubero, M.L.; Kosaka, P.M.; San Paulo, A.; Malumbres, M.; Calleja, M.; Tamayo, J. Effects of energy metabolism on the mechanical properties of breast cancer cells. Commun. Biol. 2020, 3, 1–9. [Google Scholar] [CrossRef]
  49. Liberti, M.V.; Locasale, J.W. The Warburg Effect: How Does it Benefit Cancer Cells? Trends Biochem. Sci. 2016, 41, 211–218. [Google Scholar] [CrossRef] [Green Version]
  50. Tan, A.S.; Baty, J.W.; Dong, L.F.; Bezawork-Geleta, A.; Endaya, B.; Goodwin, J.; Bajzikova, M.; Kovarova, J.; Peterka, M.; Yan, B.; et al. Mitochondrial genome acquisition restores respiratory function and tumorigenic potential of cancer cells without mitochondrial DNA. Cell Metab. 2015, 21, 81–94. [Google Scholar] [CrossRef] [Green Version]
  51. Bajzikova, M.; Kovarova, J.; Coelho, A.R.; Boukalova, S.; Oh, S.; Rohlenova, K.; Svec, D.; Hubackova, S.; Endaya, B.; Judasova, K.; et al. Reactivation of Dihydroorotate Dehydrogenase-Driven Pyrimidine Biosynthesis Restores Tumor Growth of Respiration-Deficient Cancer Cells. Cell Metab. 2019, 29, 399–416. [Google Scholar] [CrossRef] [Green Version]
  52. Pérez-Escuredo, J.; Dadhich, R.K.; Dhup, S.; Cacace, A.; Van Hée, V.; De Saedeleer, C.J.; Sboarina, M.; Rodriguez, F.; Fontenille, M.-J.; Brisson, L.; et al. Lactate promotes glutamine uptake and metabolism in oxidative cancer cells. Cell Cycle 2016, 15, 72–83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Hui, S.; Ghergurovich, J.M.; Morscher, R.J.; Jang, C.; Teng, X.; Lu, W.; Esparza, L.A.; Reya, T.; Zhan, L.; Guo, J.Y.; et al. Glucose feeds the TCA cycle via circulating lactate. Nature 2017, 551, 115–118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Dimmer, K.S.; Friedrich, B.; Lang, F.; Deitmer, J.W.; Bröer, S. The low-affinity monocarboxylate transporter MCT4 is adapted to the export of lactate in highly glycolytic cells. Biochem. J. 2000, 350, 219–227. [Google Scholar] [CrossRef] [PubMed]
  55. Xu, J.; Lu, Y.; Qiu, S.; Chen, Z.-N.; Fan, Z. A novel role of EMMPRIN/CD147 in transformation of quiescent fibroblasts to cancer-associated fibroblasts by breast cancer cells. Cancer Lett. 2013, 335, 380–386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Öhlund, D.; Handly-Santana, A.; Biffi, G.; Elyada, E.; Almeida, A.S.; Ponz-Sarvise, M.; Corbo, V.; Oni, T.E.; Hearn, S.A.; Lee, E.J.; et al. Distinct populations of inflammatory fibroblasts and myofibroblasts in pancreatic cancer. J. Exp. Med. 2017, 214, 579–596. [Google Scholar] [CrossRef]
  57. Valach, J.; Fík, Z.; Strnad, H.; Chovanec, M.; Plzák, J.; Čada, Z.; Szabo, P.; Šáchová, J.; Hroudová, M.; Urbanová, M.; et al. Smooth muscle actin-expressing stromal fibroblasts in head and neck squamous cell carcinoma: Increased expression of galectin-1 and induction of poor prognosis factors. Int. J. Cancer 2012, 131, 2499–2508. [Google Scholar] [CrossRef]
  58. He, J.; Baum, L.G. Galectin interactions with extracellular matrix and effects on cellular function. Methods Enzym. 2006, 417, 247–256. [Google Scholar]
  59. Wu, M.H.; Hong, H.C.; Hong, T.M.; Chiang, W.F.; Jin, Y.T.; Chen, Y.L. Targeting galectin-1 in carcinoma-associated fibroblasts inhibits oral squamous cell carcinoma metastasis by downregulating MCP-1/CCL2 expression. Clin. Cancer Res. 2011, 17, 1306–1316. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Lineage specificity of cancer-associated fibroblast (CAF) culture and cancer cells. (a) schematic of stromal cell isolation from tumour tissue fragments. (b) CD90 lineage specificity of CAF culture, positivity shown for viable = SYTOX™ Blue-negative cell subpopulation (SYTOX™ Blue stains dead cells with compromised plasma membranes but will not cross intact cell membranes), for more cases see Figure S1. (c) morphology of CAFs, scalebar indicate 10 μm. (d) CD44/CD90 status of FaDu and Detroit 562 head and neck cancer cell lines. (e) mRNA subtype-relevant genes, mean expression in FaDu and Detroit 562 cancer cell lines.
Figure 1. Lineage specificity of cancer-associated fibroblast (CAF) culture and cancer cells. (a) schematic of stromal cell isolation from tumour tissue fragments. (b) CD90 lineage specificity of CAF culture, positivity shown for viable = SYTOX™ Blue-negative cell subpopulation (SYTOX™ Blue stains dead cells with compromised plasma membranes but will not cross intact cell membranes), for more cases see Figure S1. (c) morphology of CAFs, scalebar indicate 10 μm. (d) CD44/CD90 status of FaDu and Detroit 562 head and neck cancer cell lines. (e) mRNA subtype-relevant genes, mean expression in FaDu and Detroit 562 cancer cell lines.
Cancers 14 02286 g001
Figure 2. Cancer cell colony-forming capacity and lactate shuttling are reprogrammed by CAFs through coculture. (a) schematic of Transwell coculture and indirect coculture with conditioned media. (b) colony-forming assay, cancer cells Transwell-cultured with CAFs; colony-supporting indicate p < 0.05 for colony size vs. non-cocultured control (see Figure S2). Measurements were performed in triplicate. (c) cancer cell colony area affected by CAFs derived from tumour tissues of varying mRNA subtypes. CL, classical, AT, atypical, BA, basal, ME, mesenchymal. CAFs derived from tumour tissues of different mRNA subtypes affect FaDu (d) and Detroit 562 (e) glycolytic ATP production differently. Performed in triplicate; PER = Proton efflux rate, rate of protons extruded into the extracellular medium during glycolysis. (f) expression of lactate-shuttle-relevant genes in cancer cell lines is associated with tumour tissue-of-origin mRNA subtype and colony-supporting status after CAF coculture. (g) correlation heatmap of the colony-forming assay with lactate-shuttle-relevant genes, lactate in medium, cells’ Young’s modulus and ATP production from glycolysis. An asterisk indicates a statistically significant correlation at p < 0.05 according to the Fischer-exact test for cluster vs. colony-supporting status and mRNA subtype. (h) Division rate of FaDu cells cocultured with CAFs relative to non-cocultured cells (0 indicates non-cocultured). N indicate the number of divisions detected in timelapse in 74 fields of view. (i) cell division rate correlates negatively with MCT4 in FaDu. OM = oligomycin; R = rotenone; AA = antimycin A; mRNA sub. = mRNA subtype; ERT-DC = Young’s modulus, real-time deformability cytometry, E (AFM) = Young’s modulus, atomic force microscopy.
Figure 2. Cancer cell colony-forming capacity and lactate shuttling are reprogrammed by CAFs through coculture. (a) schematic of Transwell coculture and indirect coculture with conditioned media. (b) colony-forming assay, cancer cells Transwell-cultured with CAFs; colony-supporting indicate p < 0.05 for colony size vs. non-cocultured control (see Figure S2). Measurements were performed in triplicate. (c) cancer cell colony area affected by CAFs derived from tumour tissues of varying mRNA subtypes. CL, classical, AT, atypical, BA, basal, ME, mesenchymal. CAFs derived from tumour tissues of different mRNA subtypes affect FaDu (d) and Detroit 562 (e) glycolytic ATP production differently. Performed in triplicate; PER = Proton efflux rate, rate of protons extruded into the extracellular medium during glycolysis. (f) expression of lactate-shuttle-relevant genes in cancer cell lines is associated with tumour tissue-of-origin mRNA subtype and colony-supporting status after CAF coculture. (g) correlation heatmap of the colony-forming assay with lactate-shuttle-relevant genes, lactate in medium, cells’ Young’s modulus and ATP production from glycolysis. An asterisk indicates a statistically significant correlation at p < 0.05 according to the Fischer-exact test for cluster vs. colony-supporting status and mRNA subtype. (h) Division rate of FaDu cells cocultured with CAFs relative to non-cocultured cells (0 indicates non-cocultured). N indicate the number of divisions detected in timelapse in 74 fields of view. (i) cell division rate correlates negatively with MCT4 in FaDu. OM = oligomycin; R = rotenone; AA = antimycin A; mRNA sub. = mRNA subtype; ERT-DC = Young’s modulus, real-time deformability cytometry, E (AFM) = Young’s modulus, atomic force microscopy.
Cancers 14 02286 g002
Figure 3. Cancer cell stiffness is affected by CAF coculture. (a) schematic of indirect (conditioned media-based) coculture experimental scheme used for cancer cell mechanophenotypisation. (b) Real-time deformability cytometry (RT-DC) scatterplots show the dependence of cell deformation on cell size. One dot represents a single cell. Representative scatterplots of FaDu cells cocultured with a CAF-conditioned medium (left) and contour plots of these scatterplots with representative cells (image width 25 µm). CAFs were derived from four tumour tissue mRNA subtypes: BA, basal; AT, atypical; CL, classical; ME, mesenchymal. For all patients tested see Figure S4. (c) Young’s modulus determined by RT-DC (ERT-DC) for four major mRNA subtypes. Performed in duplicates per CAF coculture treatment, N indicates the number of cocultures, asterisk indicates statistical significance at p < 0.05. (d) Young’s modulus determined by RT-DC and atomic force microscopy weakly correlates; for AFM measurements see Figure S4. (e) The oxygen consumption rate (OCR) of FaDu cells is affected by the CAF mRNA subtype and is in negative correlation with Young’s modulus (g), and MCT1 FaDu expression (f) OM = oligomycin; R = rotenone; AA = antimycin A; ERT-DC = Young modulus, real-time deformability cytometry, EAFM=Young modulus, atomic force microscopy.
Figure 3. Cancer cell stiffness is affected by CAF coculture. (a) schematic of indirect (conditioned media-based) coculture experimental scheme used for cancer cell mechanophenotypisation. (b) Real-time deformability cytometry (RT-DC) scatterplots show the dependence of cell deformation on cell size. One dot represents a single cell. Representative scatterplots of FaDu cells cocultured with a CAF-conditioned medium (left) and contour plots of these scatterplots with representative cells (image width 25 µm). CAFs were derived from four tumour tissue mRNA subtypes: BA, basal; AT, atypical; CL, classical; ME, mesenchymal. For all patients tested see Figure S4. (c) Young’s modulus determined by RT-DC (ERT-DC) for four major mRNA subtypes. Performed in duplicates per CAF coculture treatment, N indicates the number of cocultures, asterisk indicates statistical significance at p < 0.05. (d) Young’s modulus determined by RT-DC and atomic force microscopy weakly correlates; for AFM measurements see Figure S4. (e) The oxygen consumption rate (OCR) of FaDu cells is affected by the CAF mRNA subtype and is in negative correlation with Young’s modulus (g), and MCT1 FaDu expression (f) OM = oligomycin; R = rotenone; AA = antimycin A; ERT-DC = Young modulus, real-time deformability cytometry, EAFM=Young modulus, atomic force microscopy.
Cancers 14 02286 g003
Figure 4. Coculture of the CAFs with cancer cell lines increases lactate in medium and lactate shuttle-relevant gene expression in CAFs. (a) scheme of measurement, (b) levels of lactate in medium and lactate-shuttle-relevant genes expression level in CAFs, asterisk indicate p < 0.05 in paired t-test. (c) MCT1 and MCT4 gene expression in tumour tissue and corresponding tumour tissue-derived CAFs. p indicates the result of paired t-test, individual dots represent individual patient-derived CAFs. (d) summary of the coculture data –shown as an integrative heatmap of measured data on cocultured cells, shown as a mean-per-mRNA subtype (CL, classical, AT, atypical, BA, basal, ME, mesenchymal).
Figure 4. Coculture of the CAFs with cancer cell lines increases lactate in medium and lactate shuttle-relevant gene expression in CAFs. (a) scheme of measurement, (b) levels of lactate in medium and lactate-shuttle-relevant genes expression level in CAFs, asterisk indicate p < 0.05 in paired t-test. (c) MCT1 and MCT4 gene expression in tumour tissue and corresponding tumour tissue-derived CAFs. p indicates the result of paired t-test, individual dots represent individual patient-derived CAFs. (d) summary of the coculture data –shown as an integrative heatmap of measured data on cocultured cells, shown as a mean-per-mRNA subtype (CL, classical, AT, atypical, BA, basal, ME, mesenchymal).
Cancers 14 02286 g004
Figure 5. Lactate shuttle-relevant gene expression clusters in tumour tissues are associated with patient survival. (a) heatmap of tumour tissue lactate shuttle-relevant gene expression in 55 primary head and neck tumours. p values on the right side indicate Fischer’s exact test for p16 status tissue site, node status, cancer stage, smoking, and mRNA subtype. Clustering results in two clusters, which demonstrated an overall survival trend in the Kaplan-Meier chart (b). (c) Cox regression, stepwise model. Gene expression cluster is an overall survival predictor independent of tumour node status (pN), p16 status, and tumour stage, which were used as input in the Cox model. Asterisk highlights displayed p-values lower than 0.05. (d) MCT1 and MCT4 expression in the stroma and tumorous tissues, representative basal, mesenchymal, and atypical mRNA subtype tumours. MCT1 cytoplasmic positivity and isolated membrane positivity in tumour cells, of AT tumours (first column), MCT4 weak positivity in tumour tissue and strong positivity of inflammatory elements and sparse stromal cells in peritumoral stroma (S) on a periphery of tumour nests (T, indicated by the arrow, second column); MCT1 positivity of isolated tumour cells in a basal layer of the nests, BA tumours (indicated by the arrow, third column); MCT1 negativity in tumorous tissues (last column). 40×, scalebar indicates 100 μm, detail width 80 μm. (e) MCT1/4 tumour tissue positivity and mRNA subtype in immunohistochemical staining.
Figure 5. Lactate shuttle-relevant gene expression clusters in tumour tissues are associated with patient survival. (a) heatmap of tumour tissue lactate shuttle-relevant gene expression in 55 primary head and neck tumours. p values on the right side indicate Fischer’s exact test for p16 status tissue site, node status, cancer stage, smoking, and mRNA subtype. Clustering results in two clusters, which demonstrated an overall survival trend in the Kaplan-Meier chart (b). (c) Cox regression, stepwise model. Gene expression cluster is an overall survival predictor independent of tumour node status (pN), p16 status, and tumour stage, which were used as input in the Cox model. Asterisk highlights displayed p-values lower than 0.05. (d) MCT1 and MCT4 expression in the stroma and tumorous tissues, representative basal, mesenchymal, and atypical mRNA subtype tumours. MCT1 cytoplasmic positivity and isolated membrane positivity in tumour cells, of AT tumours (first column), MCT4 weak positivity in tumour tissue and strong positivity of inflammatory elements and sparse stromal cells in peritumoral stroma (S) on a periphery of tumour nests (T, indicated by the arrow, second column); MCT1 positivity of isolated tumour cells in a basal layer of the nests, BA tumours (indicated by the arrow, third column); MCT1 negativity in tumorous tissues (last column). 40×, scalebar indicates 100 μm, detail width 80 μm. (e) MCT1/4 tumour tissue positivity and mRNA subtype in immunohistochemical staining.
Cancers 14 02286 g005
Table 1. Clinico-pathological characteristics of 55 patient cohort. Shown as number (% per group). p-values for Fischer exact test for factor vs. mRNA subtype and for a lactate shuttle-relevant gene expression clusters (see further in Figure 5), respectively. AT, atypical, BA, basal, CL, classical, ME, mesenchymal. N, number of patients.
Table 1. Clinico-pathological characteristics of 55 patient cohort. Shown as number (% per group). p-values for Fischer exact test for factor vs. mRNA subtype and for a lactate shuttle-relevant gene expression clusters (see further in Figure 5), respectively. AT, atypical, BA, basal, CL, classical, ME, mesenchymal. N, number of patients.
Factor/LevelNmRNA SubtypeLactate Shuttle Gene Cluster
ATBACLME12
Gender, p-value = 0.974/0.519
F144 (28.6%)3 (21.4%)4 (28.6%)3 (21.4%)9 (64.3%)5 (35.7%)
M4114 (34.1%)8 (19.5%)10 (24.4%)9 (22.0%)30 (73.2%)11 (26.8%)
Tumor location, p-value < 0.001/0.286
floor of the mouth30 (0.0%)1 (33.3%)1 (33.3%)1 (33.3%)2 (66.7%)1 (33.3%)
hypopharynx61 (16.7%)0 (0.0%)3 (50.0%)2 (33.3%)4 (66.7%)2 (33.3%)
larynx160 (0.0%)5 (31.2%)9 (56.2%)2 (12.5%)12 (75.0%)4 (25.0%)
oral cavity10 (0.0%)1 (100.0%)0 (0.0%)0 (0.0%)0 (0.0%)1 (100.0%)
oropharynx2417 (70.8%)2 (8.3%)1 (4.2%)4 (16.7%)19 (79.2%)5 (20.8%)
tongue50 (0.0%)2 (40.0%)0 (0.0%)3 (60.0%)2 (40.0%)3 (60.0%)
p16 status, p < 0.001/0.3594
p16 neg.350 (0.0%)10 (28.6%)14 (40.0%)11 (31.4%)23 (65.7%)12 (34.3%)
p16 pos.2018 (90.0%)1 (5.0%)0 (0.0%)1 (5.0%)16 (80.0%)4 (20.0%)
pN, p-value = 0.03059/1.000
<23712 (32.4%)9 (24.3%)12 (32.4%)4 (10.8%)26 (70.3%)11 (29.7%)
>2186 (33.3%)2 (11.1%)2 (11.1%)8 (44.4%)13 (72.2%)5 (27.8%)
Stage, p value < 0.001/1.000
I–II2816 (57.1%)5 (17.9%)3 (10.7%)4 (14.3%)20 (71.4%)8 (28.6%)
III–IV272 (7.4%)6 (22.2%)11 (40.7%)8 (29.6%)19 (70.4%)8 (29.6%)
Smoking status, p < 0.001/0.536
01814 (77.8%)0 (0.0%)0 (0.0%)4 (22.2%)14 (77.8%)4 (22.2%)
1374 (10.8%)11 (29.7%)14 (37.8%)8 (21.6%)25 (67.6%)12 (32.4%)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peltanová, B.; Holcová Polanská, H.; Raudenská, M.; Balvan, J.; Navrátil, J.; Vičar, T.; Gumulec, J.; Čechová, B.; Kräter, M.; Guck, J.; et al. mRNA Subtype of Cancer-Associated Fibroblasts Significantly Affects Key Characteristics of Head and Neck Cancer Cells. Cancers 2022, 14, 2286. https://doi.org/10.3390/cancers14092286

AMA Style

Peltanová B, Holcová Polanská H, Raudenská M, Balvan J, Navrátil J, Vičar T, Gumulec J, Čechová B, Kräter M, Guck J, et al. mRNA Subtype of Cancer-Associated Fibroblasts Significantly Affects Key Characteristics of Head and Neck Cancer Cells. Cancers. 2022; 14(9):2286. https://doi.org/10.3390/cancers14092286

Chicago/Turabian Style

Peltanová, Barbora, Hana Holcová Polanská, Martina Raudenská, Jan Balvan, Jiří Navrátil, Tomáš Vičar, Jaromír Gumulec, Barbora Čechová, Martin Kräter, Jochen Guck, and et al. 2022. "mRNA Subtype of Cancer-Associated Fibroblasts Significantly Affects Key Characteristics of Head and Neck Cancer Cells" Cancers 14, no. 9: 2286. https://doi.org/10.3390/cancers14092286

APA Style

Peltanová, B., Holcová Polanská, H., Raudenská, M., Balvan, J., Navrátil, J., Vičar, T., Gumulec, J., Čechová, B., Kräter, M., Guck, J., Kalfeřt, D., Grega, M., Plzák, J., Betka, J., & Masařík, M. (2022). mRNA Subtype of Cancer-Associated Fibroblasts Significantly Affects Key Characteristics of Head and Neck Cancer Cells. Cancers, 14(9), 2286. https://doi.org/10.3390/cancers14092286

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop