Next Article in Journal
Where Artificial Intelligence Can Take Us in the Management and Understanding of Cancerization Fields
Next Article in Special Issue
Does Genotype-Specific Phenotype in Patients with Multiple Endocrine Neoplasia Type 2 Occur as Current Guidelines Predict?
Previous Article in Journal
Special Issue Editorial: Emerging Therapies in Renal Cell Carcinoma: The Road to a Cure?
Previous Article in Special Issue
Nasopharynx Battlefield: Cellular Immune Responses Mediated by Midkine in Nasopharyngeal Carcinoma and COVID-19
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Targeting Translation and the Cell Cycle Inversely Affects CTC Metabolism but Not Metastasis

1
Division of Molecular Medicine, Department of Internal Medicine, University of New Mexico Health Sciences Center, Albuquerque, NM 87131, USA
2
Animal Models Shared Resource, University of New Mexico Comprehensive Cancer Center, Albuquerque, NM 87120, USA
3
Division of Hematology and Oncology, Department of Internal Medicine, University of New Mexico Comprehensive Cancer Center, Albuquerque, NM 87120, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2023, 15(21), 5263; https://doi.org/10.3390/cancers15215263
Submission received: 27 September 2023 / Revised: 26 October 2023 / Accepted: 30 October 2023 / Published: 2 November 2023
(This article belongs to the Special Issue Recent Advances in Rare Cancers: From Bench to Bedside and Back)

Abstract

:

Simple Summary

This study identifies specific ribosomal proteins that drive MBM and extracranial metastasis via CTCs. Dual targeting of cellular translation and proliferation prevents metastasis of aggressive CTC subsets with high RPL/RPS expression. We report the first-ever real-time metabolic flux analysis of patient-derived melanoma CTCs and altered carbohydrate metabolism during impaired translation in a melanoma-CTC-derived clone.

Abstract

Melanoma brain metastasis (MBM) is significantly associated with poor prognosis and is diagnosed in 80% of patients at autopsy. Circulating tumor cells (CTCs) are “seeds” of metastasis and the smallest functional units of cancer. Our multilevel approach has previously identified a CTC RPL/RPS gene signature directly linked to MBM onset. We hypothesized that targeting ribogenesis prevents MBM/metastasis in CTC-derived xenografts. We treated parallel cohorts of MBM mice with FDA-approved protein translation inhibitor omacetaxine with or without CDK4/CDK6 inhibitor palbociclib, and monitored metastatic development and cell proliferation. Necropsies and IVIS imaging showed decreased MBM/extracranial metastasis in drug-treated mice, and RNA-Seq on mouse-blood-derived CTCs revealed downregulation of four RPL/RPS genes. However, mitochondrial stress tests and RT-qPCR showed that omacetaxine and palbociclib inversely affected glycolytic metabolism, demonstrating that dual targeting of cell translation/proliferation is critical to suppress plasticity in metastasis-competent CTCs. Equally relevant, we provide the first-ever functional metabolic characterization of patient-derived circulating neoplastic cells/CTCs.

Graphical Abstract

1. Introduction

Melanoma is the most aggressive form of skin cancer with poor therapeutic options and clinical outcomes [1,2,3]. Metastatic dissemination to the brain occurs in 60% of advanced melanoma patients and in up to 80% of patients at autopsy [1,4]. Patients with melanoma brain metastasis (MBM) have a life expectancy of only 4–6 months, associated with a rapid decrease in life quality [1,2,5,6]. Localized therapies include surgical removal and radiation, while metastatic melanoma is effectively treated with immunotherapy and targeted agents [1,7,8]. Although innovative therapeutic breakthroughs are constantly developed, substantial advances remain elusive and most MBM patients die due to MBM progression, hemorrhage, and intra- and extracranial metastases [9,10].
Metastatic dissemination to the brain and other organ sites occurs due to continuous shedding of circulating tumor cells (CTCs) from primary and/or metastatic tumors into the bloodstream [2,11,12,13]. CTCs are “seeds” of fatal metastasis and the smallest functional units of cancer. Hypoxia is one of the major factors that triggers CTC shedding and enhanced metastasis [14,15]. Due to the detrimental effects of shear stress in blood, CTC clusters are dissociated and the vast majority of CTCs die prior to successful dissemination to distant organs [14,16,17,18,19]. CTCs that do disseminate can remain dormant for years, e.g., in the bone marrow as depot, or initiate clinically undetectable micrometastases [20,21,22]. CTCs employ diverse invasive behaviors, such as amoeboid cellular invasion, CTC cluster streaming, or collective cell migration and invasion [14]. The latter is the major CTC behavioral method for establishment and progression of a solid tumor [14]. Significantly, high numbers of single CTCs and/or CTC clusters are directly associated with cancer severity and clinical outcomes [12,23,24,25,26]. CTC clusters are composed of leader and follower cells, with leaders preparing a low-resistance path to shield followers from biomechanical stressors in the blood such as shear stress, tension, and compression [14,27]. By traveling with companions, including immune and stromal cells, CTCs can increase their survival chances in the harsh microenvironment of blood [14]. We have previously reported that growing MBM induces exponential shedding of CTCs into the bloodstream [28].
We employed an unbiased multilevel approach to identify a unique 21-member RPL/RPS signature of melanoma CTCs that drives MBM onset and progression (“RPL” means proteins of the large ribosomal 60S subunit, while “RPS” refers to proteins of the small ribosomal 40S subunit [28]). Translational reprogramming of cancer cells occurs early on in tumorigenesis to initiate tumor growth, dissemination of CTCs into the bloodstream, and metastasis [29]. Cancer cells can rewire ribogenesis to preferentially translate oncogenic transcripts and promote metastatic spread, neoplastic growth, and survival, and this deregulation has been linked to increased metastatic spread of cancer [29,30,31,32]. For example, upregulation of RPL15 gene expression in breast cancer promotes massive metastasis and triggers enhanced translation of other ribosomal proteins [29,30], while inhibition of RPL27A hampers dissemination and invasive potential of breast cancer cells [33].
In synchrony with the above studies, we have identified 21 RPL/RPS genes significantly upregulated in CTCs that drive MBM onset [28]. Because there is highly regulated cooperation between ribosome assembly and cell-cycle progression, and since inhibition of translation blocks cell proliferation [31,34], we hypothesized that targeting ribosome biogenesis and/or cell proliferation prevents MBM development in CTC-derived xenografts (CDXs). We evaluated this by inhibiting these two vital processes in MBM CDXs using the translational inhibitor omacetaxine and/or the cell-cycle inhibitor palbociclib. Omacetaxine inhibits peptide elongation by directly binding to the A-site cleft of the ribosome and, while the exact mechanism of ribosomal inhibition remains to be elucidated [35], it is FDA-approved for treating chronic myelogenous leukemia [36,37]. Conversely, palbociclib is a CDK4/6 inhibitor that has been widely used in many cancer therapies [38]. Here, we demonstrate that targeting either the cell translational machinery and/or cell proliferation results in a significant decrease in MBM/metastasis. Additionally, we singled out 4 key players from the 21-member CTC RPL/RPS gene signature as significantly downregulated following drug treatment and revealed carbohydrate metabolism to be dysregulated during impaired ribogenesis. Metabolic rewiring is fundamental to metastatic progression and metastasis-competent CTCs are hypothesized to be capable of flexibly deploying both oxidative and anabolic metabolism to meet energy demands, overcome oxidative stress, and expand available metastatic niches [39,40,41,42,43,44]. We report the first-ever metabolic flux analysis of patient-derived neoplastic cell fractions/CTCs in real-time, alongside omacetaxine and/or palbociclib treatment of a melanoma CTC-derived clone to elucidate melanoma CTC metabolism during inhibition of ribogenesis.

2. Materials and Methods

2.1. Patients’ Blood Collection

Metastatic melanoma patient samples were accrued according to protocols approved by the Institutional Review Board at UNM Health Sciences Center (UNM-HSC) in Albuquerque, New Mexico. Patients’ clinical parameters are shown in Table S4. Patient samples were collected after completed informed written consent was obtained. Under strict aseptic conditions, peripheral blood (12–18 mL) from each patient was collected at the middle of vein puncture in a sodium-etheylenediamine tetraacetic acid (EDTA) tube as part of their doctor’s visit. Blood samples were provided immediately to the laboratory after collection for CTC isolation and examination, and processed immediately or fixed in a CellSaveTM tube (Menarini Silicon Biosystems, Inc., Bologna, Italy).

2.2. CTC-Derived Xenografts (CDXs)

All in vivo studies were approved by Institutional Animal Care and Use Committee (IACUC) protocol. The generation of CTC-derived xenografts (CDXs) was executed using 6 to 8 week old NBSGW (NOD.Cg-KitW-41J Tyr + Prkdcscid Il2rgtm1Wjl/ThomJ) immunodeficient mice. Thirty mice received 50 µL (4 mg/mL) of low-molecular weight heparin (NDC, Cat# 63323-533-01, Lake Zurich, IL, USA) retro-orbitally 10 min before intracardiac injection of CTC-derived CTC-Luc2 clone to decrease thromboembolism-related mortality in mice [45]. Following anesthesia with isoflurane (2.5%, 1 L/min O2 flow), animals were injected under aseptic conditions into the left ventricle with 5.0 × 105 cells in 50 µL of PBS using a sterile 0.5 mL U-100 insulin syringe with a 29G×/2″ needle (Beckton Dickinson, Cat# 58324702, Franklin Lakes, NJ, USA). Successful intracardiac injection was confirmed by blood backflow into the syringe prior to cell injection. MBM development was confirmed as early as 24 h post-injection using the Xenogen IVIS Spectrum animal imager (PerkinElmer, Waltham, MA, USA) following intraperitoneal injection of D-luciferin (Sigma, Cat# 2591-17-5; 150 mg/kg, St. Louis, MO, USA). IVIS imaging was performed on weekly basis to monitor tumor development. One week post-injection, 20 MBM mice were equally divided between four cohorts: (1) placebo mice which received vehicle (n = 5); (2) omacetaxine mice which received IP injection of omacetaxine (Chem Scene, Cat# CS-2872; 0.5 mg/kg, Monmouth Junction, NJ, USA), n = 5; (3) palbociclib mice which received oral gavage of the palbociclib (Selleckchem, PD-0332991; 25 mg/kg, Houston, TX, USA), n = 5; (4) mice which received combined treatment, n = 5. Drug treatment was performed every weekday, as previously described, over a period of three weeks [38,46]. After daily monitoring of animals over the period of four weeks for changes in overall health status (sudden weight loss, rough coat, distress, difficulty with ambulation, seizures, or difficulty in obtaining food or water), all the animals were sacrificed for subsequent necropsies. Approximately 800 to 900 µL of mouse blood was collected into an EDTA-containing MiniCollect tube (Greiner Bio-One, Cat# K3E K3EDTA, Monroe, NC, USA). Brain, liver, lungs, spleen, sternum, and skull-cap tissues were fixed in 10% neutral buffered formalin for further pathological assessment. Since most melanoma cells produce melanin, melanoma macrometastases were visually identified as brown-to-black pigmented areas [47].

2.3. CTC/CTC Cluster Capture

Mouse blood (100–150 µL) was collected retro-orbitally using EDTA-coated glass Pasteur pipette into a MiniCollect tube (Greiner Bio-One, Cat# K3E K3EDTA, Monroe, NC, USA) and loaded onto the CTC ParsortixTM microfluidic chip (8 µm) within 1 h of mouse necropsies. CTC and CTC clusters were captured by the filtration/microfluidic Parsortix PR1 platform (Angle Europe Ltd., Surrey, UK) using 6.5 µM cartridges (Angle, PLC). For immunofluorescence staining, CTCs were captured by Parsortix based upon their size and deformability, and subjected to on-cassette staining, as previously described [25,28]. CTCs were designated as cells positive for human Mel-A (Alexa Fluor 594-tagged, Santa Cruz Biotechnology, Cat#sc-20032, Dallas, TX, USA) and human DAPI (Thermo Fisher, Cat# D3571, Waltham, MA, USA) staining, but negative for human CD45 (FITC-tagged, BioLegend, Cat#103108, San Diego, CA, USA). CTC visualization and quantification were performed using Zeiss LSM800 microscope (10–40× magnification) and ZEN Blue 2.1 software (Carl Zeiss Microscopy; Oberkochen, Germany). For RNA-Seq analyses, CTCs were captured and harvested using Parsortix, according to manufacturer’s guidelines.

2.4. RNA Sequencing

RNA was isolated from enriched and flow-through fractions (10–20 × 103 cells) after ParsortixTM CTC harvesting (ANGLE, LLC, Plymouth Meeting, PA, USA). Enriched fraction contained human CTCs, while flow-through fraction consisted of immune cells. RNA analysis, cDNA synthesis, and library preparation were performed using the microarray platform SMARTer Universal Low Input RNA kit for RNA sequencing (Clontech, Cat#634946, San Jose, CA, USA). Fragmented RNA was processed using the Ion Plus Fragment Library kit (Thermo Fisher, Cat#4471252, Waltham, MA, USA), as previously reported [48,49]. The samples were sequenced using the Ion Proton S5/XL platform (Thermo Fisher, Waltham, MA, USA) at the Analytical and Translational Genomics (ATG) Shared Resource Core at the University of New Mexico Comprehensive Cancer Center (UNM-CCC).

2.5. RNA-Seq Analyses and Bioinformatics

Sequence alignment was performed using tmap (v5.10.11) to a BED file carrying non-overlapping exon regions from the UCSC genome browser (GRCh38/hg38). Exons were quantified using the HTSeq (v0.11.1) Python platform [50,51]. Generation of gene counts was performed by summing counts across exons. Library size was normalized and differential analysis was executed using edgeR [50,52]. An adjusted p-value of cutoff of 0.05 was used for analyzing edgeR and DESeq values. Heatmap3 was used to generate heatmaps and execute cluster analyses. Pathway enrichment analyses were conducted using clusterProfiler, GAGE, Pathview, and topGO software programs, R Package Version 2.42.0 [50,52]. The Reactome pathway database was used to generate pathway discrimination analyses [53]. The crystal structure of the ribosome and its components was generated using the PyMOL molecular visualization system, with the specific structure code being 4UG0 [54]. The relevant proteins were colored using PyMOL and the individually referenced proteins within the protein databank [55].

2.6. Cell Culture

Early-passage brain-metastatic melanoma CTC-derived clonal cell line (CTC clone; generated in Dr. Marchetti’s laboratory) [28,56] or human lung-metastatic melanoma MeWo cell line (ATCC, Cat# HTB-65, Manassas, VA, USA) were cultured in DMEM/F12 (Thermo Fisher, Cat# 11320082), supplemented with 10% fetal bovine serum (FBS) (Gibco, Cat# A4766801, Billings, MT, USA) and 1X GlutaMAX (Gibco, Cat# 35050-061), termed “standard media” (Figure S2B). Cells were maintained at 37 °C in a humidified 5% CO2 incubator and passaged using 0.05% Trypsin-EDTA (Gibco, Cat# 25300054) before reaching confluency. Routine testing for Mycoplasma was performed using the Mycoplasma Detection Assay (MycoAlert, Lonza, Morristown, NJ, USA) every 20 passages. Cell-line authentication was conducted using PCR-based assay routinely. Luciferase-tagged CTC clone cells were generated as previously reported [57]. Prior to injection into mice, cells were checked for morphological and phenotypic changes using microscopy.

2.7. Blood Separation and Lineage Sorting

Whole blood obtained as above was subjected to magnetic labeling for erythrocyte and granulocyte depletion per the manufacturer’s instructions using the Custom Whole Blood PBMC Iso Kit (Miltenyi Cat# 130-126-359, Auburn, CA, USA), autoMACS Running Buffer (Miltenyi, Cat# 130-091-221), and Custom Sedimentation Kit2 (Miltenyi, Cat# 130-126-357). Lineage sorting was performed using LS Positive Selection Columns (Miltenyi, Cat# 130-126-357), anti-CD45 microbeads (Miltenyi, Cat# 130-045-801), and the QuadroMACs Magnetic Separator (Miltenyi, Cat# 130-091-051). After granulocyte and erythrocyte depletion, cells were washed in autoMACS buffer, and re-suspended in 300 µL of autoMACS and 100 µL of Fc Block (BioLegend, Cat# 422302, San Diego, CA, USA), each per 5 × 107 cells, and incubated at 4 °C for 10 min. An amount of 100 µL of CD45 microbeads per 5 × 107 cells were added prior to 30 min incubation at 4 °C. At 15 min into the 30 min incubation, 100 µL of Dead Cell Removal Beads (Miltenyi, Cat# 130-090-101) per 1 × 107 cells was added. Fractions were then subject to magnetic separation in LS Positive Selection Columns and washed with PBS, and the CD45− and CD45+ (Lin− and Lin+, respectively) fractions were then used for metabolic flux analysis.

2.8. Metabolic Flux Analysis on Patient-Derived PBMCs and CTC-Derived Clone

Mitochondrial stress tests (MSTs) (Agilent, Cat# 103015-100, Santa Clara, CA, USA) and Glycolytic Rate Assays (GRAs) (Agilent, Cat# 103344-100) were performed according to the manufacturer’s instructions. Assay medium for all experiments was XF DMEM, pH 7.4 (Agilent, Cat# 103575-100), supplemented with 10 mM glucose (Agilent, Cat# 103577-100), 1 mM pyruvate (Agilent, Cat# 103578-100), and 2 mM glutamine (Agilent, Cat# 103579-100). MST port drugs were used at final well concentrations as follows: oligomycin (3 µM), FCCP (3 µM) [58,59], and rotenone/antimycin A (Rot/AA, 1 µM) [60]. GRA port drugs were 2-deoxyglucose (2-DG, 5 µM) and Rot/AA (500 mM). Post-assay normalization on suspension cells was performed by including Hoechst Dye (Life Technologies, Cat# H3570) with the final drug injection of the SeahorseTM assay (Agilent, Inc.), yielding a final well concentration of 16 µM. Plates were then analyzed using a BioTek Cytation 1 Cell Imaging Multimode Reader running Gen5 3.11 Software programmed to output MeanTsf values, which were multiplied by input viability to yield the final normalization unit for PBMC assays. Normalization was performed using Agilent Wave 2.6.1.56 Software. Prior to use, Seahorse Cell Culture Plates (Agilent, Cat# 102416-100) were coated overnight in a Biosafety Cabinet with 100 µL per well of 100 µg/mL PureCol-S Solution (ECM Protein, Cat# CC300), and the next day washed with 100 µL PBS per well and allowed to dry for at least 2 h prior to suspension cell seeding. Patient-derived whole blood was magnetically depleted of erythrocytes and granulocytes as above prior to labeling with CD45 magnetic beads and sorting into CD45− (Lin−) and CD45+ (Lin+) populations. After washing with PBS, cells were re-suspended in assay medium at a density of 1 × 105 cells/100 µL per well into the collagen-coated plates and centrifuged for 1 min at 200× g with 0 braking. An amount of 80 µL of assay medium was then added to bring the final well volume to 180 µL. Cells were placed in a non-CO2 incubator at 37C for 1 h prior to running flux analysis.
The CTC-derived clone and MeWo control cells were plated at a density of 1 × 104/100 µL in 96-well Seahorse Cell Culture Plates, and allowed to settle for one hour at room temperature before 100 µL of drug or placebo was added to yield final well concentrations of 1.5 µM palbociclib (LabShops, Cat# C1250-100mg, Claymont, DE, USA), 0.5 µM omacetaxine (ChemScene, Cat# 26833-87-4, Monmouth Junction, NJ, USA), both, or neither. Cells were then incubated in standard conditions for 3 days prior to flux analysis. For substrate oxidation experiments, 10 µM BPTES (Agilent, Cat# 103674-100, Santa Clara, CA, USA) or 10 µM UK5099 (Agilent, Cat# 103873-100) was added 24 h before flux analysis. For glutamine deprivation during drug treatment, DMEM without Phenol Red, D-glucose, L-glutamine, or sodium pyruvate (Gibco, Cat# A14430-01, Billings, MT, USA) was manually supplemented with 10 mM glucose (Agilent, Cat# 103577-100), 1 mM pyruvate (Agilent, Cat# 103578-100), 10% FBS, and 5 µM Hepes (Cytiva, Cat# SH30237.01, Marloboro MA, USA), with or without 2 mM glutamine (Agilent, Cat# 103579-100). We confirmed that this custom-supplemented media recapitulated the effects of our drug treatments compared to standard media, such that effects of glutamine deprivation were not due to the change in media formulation (Figure S2B). Adherent cells were washed twice with 200 µL of supplemented assay medium and brought to a final volume of 180 µL per well prior to non-CO2 incubation. For adherent cells after 3 day drug treatment, flux data were normalized to MeanTsf values.

2.9. mRNA Isolation and RT-qPCR

Total mRNA was isolated using the RNeasy Micro Kit (Qiagen, Cat# 74004, Hilden, Germany) per manufacturer’s instructions. A maximum of 5 × 105 cells were input for lysis and RNA extraction, and stored at −80 °C until ready for analysis. Amounts of 2 ng of total RNA per well were loaded into 0.1 mL MicroAMP fast 96-well Reaction Plates (Applied Biosystems, Cat# 4346907, Waltham, MA, USA) in 20 µL/well reactions and sealed with Optical Adhesive Covers (Applied Biosystems, Cat# 4360954). Each 20 µL reaction contained 5 µL of Taqman Fast Virus 1-step MasterMix (Thermo Fisher, Cat# 4444434), 1 µL of normalized RNA, 1 L of 20× Taqman Gene Expression Assays, and 13 µL of RNAse-free H2O supplied in the RNeasy Microkits. Taqman Gene Expression Assays were β-Actin (Hs01060665_g1), PGC-1α (Hs00173304_m1), PFKFB4 (Hs00894603_m1), GLUT1 (Hs00892681_m1), GLS (Hs01014020_m1), LDHA (Hs01378790_g1), LDHB (Hs00929956_m1), HIF1α (H Hs00153153_m1s), and PC (Hs01085875_g1). RT-qPCR was performed in duplicate with 2 ng/µL total RNA per well and following the cycling parameters specified for Fast Virus 1-step MasterMix for 40 cycles in a QuantStudio 6 Pro RT-PCR Analyzer, and final well ΔCT values were output by QuantStudio Version 6 software. Reported mRNA fold changes reflect 2−ΔΔCT values. * p = 0.05, ** p = 0.01.

3. Results

3.1. The Translation Machinery and/or Cell Proliferation Are Critical for Successful CTC Micro/Macrometastases

First, 30 NBSGW mice were injected with an MBM CTC-derived clone. This CTC-derived melanoma clonal cell line was generated in Dr. Marchetti’s laboratory. It resulted from sequential selections (three) of cells injected intracardiacally in NSG mice, with the original clone being isolated from MBM. The CTC clonal cell line was shown to possess high-metastatic capabilities, found to consistently promote metastasis within two weeks from cell injection into animals, and to generate MBM as early as 24 h [28]. MBM mice were identified via IVIS imaging as early as 24 h post-injection. After one week, 20 MBM mice were equally divided into four cohorts (five mice in each cohort): placebo mice, mice treated with RPL inhibitor omacetaxine, mice treated with CDK4/CDK6 inhibitor palbociclib, and mice with combined drug treatment. Tumor progression was observed using IVIS which showed a higher tumor burden in placebo mice over time when compared to the drug-treated mouse groups (Figure 1A). Tridimensional reconstruction of tumor dissemination in all four animal groups was assessed using 3D IVIS tomography, showing downregulation of spatial metastasis in the drug-treated animals (Figure 1A). Notably, the MBM 3D signal was absent in the omacetaxine group with or without palbociclib. Total flux of intra- and extracranial metastases was quantified using IVIS which revealed a statistically significant decrease in signal in those mice that received individual or dual drug treatment (Figure 1B,C, with statistical analyses presented in Tables S1 and S2, respectively).
Whole body IVIS quantification of the total flux in the mouse revealed similar metastatic trends (Figure S1B). Parallel animal necropsies confirmed a higher number of MBM mice and more extensive tumor burden in placebo mice (Figure 2A). A higher number of placebo mice developed extracranial macrometastases to distant organs, such as liver and bone marrow, when compared to drug-treated cohorts (Figure 2B). Micrometastasis quantification of tumor cells in mouse brains was consistent with macrometastatic analyses—placebo brains exhibited higher infiltration of tumor cells, when compared to brains from other groups (Figure 2C). Mouse sections of placebo brains had a 2–3-fold increase in tumor cells when compared to drug-treated brain sections (Figure 2D).

3.2. Inhibition of Translation and/or Cell Proliferation Decreased CTC Numbers and Altered Gene Signatures

To examine CTCs and changes in gene expression following drug treatment, we captured and analyzed CTCs from all animal groups. In each group, equal volumes of blood from three MBM CDXs were combined following retro-orbital blood draw and analyzed using the CTC Parsortix platform at the end of the study. Human CTCs were captured from mouse blood based on their size and deformability, and were then immunostained for human Mel-A Alexa Fluor 594, human FITC-CD45, and DAPI nuclear marker. CTCs were defined as Mel-A+/CD-45-/DAPI+ cells, visualized and quantified (Figure 3A and Table 1). Importantly, drug-treated CDXs exhibited a lower number of CTCs in their blood when compared to the positive control (Figure 3 and Table 1). Second, these observations were especially true for mice treated with omacetaxine with or without palbociclib. With omacetaxine treatment, the number of single CTCs decreased 4–6-fold in comparison with control (Figure 3 and Table 1). Third, blood from drug-treated mice did not contain any small or large CTC clusters. This finding is indicative of CTC decreased metastatic capacity and better therapy response in the drug-treated groups [61] (Figure 3A and Table 1). Additionally, CTCs from metastatic melanoma patients were captured and visualized using Parsortix (Figure 3B). As a negative control, blood of healthy donors was used to confirm the absence of CTCs (Figure 3C).
In parallel to ParsortixTM staining, two blood samples from each animal group were processed by Parsortix to harvest CTCs for single-cell RNA-Seq analysis. RNA-Seq examination was performed to compare gene expression levels in different experimental groups, with libraries being aligned to the human genome. Hierarchical clustering revealed differential clustering of CTCs from placebo vs. drug-treated mouse groups (Figure 4, left panel). The displayed volcano plots show global gene expression of CTC populations with significant log2 fold change in the placebo samples versus the drug-treated samples (Figure 4, middle panel). Next, the differential gene expression between samples was generated using the ReactomeTM pathway database and resulted in the list of statistically significant pathways (Figure 4, right panel). Importantly, Reactome analyses of placebo versus omacetaxine groups generated an “Infectious disease” pathway that contains 41 RPL/RPS genes (Figure 4A, right panel—highlighted in yellow). Significantly, this pathway was present in the list of molecular pathways previously reported to be associated with the CTC RPL/RPS gene signature of MBM [28]. Another consequential pathway was “Metabolism of carbohydrates” which links RPL inhibition with metabolic alterations (Figure 4A, right panel—highlighted in red).
The second Reactome analysis against hallmark genes differentially expressed in placebo versus palbociclib generated a list of significant molecular pathways involved in the senescence and inhibition of cell proliferation (Figure 4B, right panel—highlighted in green). This result was expected since palbociclib inhibits CDK4/6 and therefore cell proliferation [62]. The third Reactome analysis between placebo and dual drug treatment produced a list of 35 upregulated significant pathways (Figure 4C, right panel). There were several important findings regarding this analysis. First, 34 out of 35 pathways had RPL/RPS genes and were involved in translational processes. Second, 33 out of 35 pathways were the exact signaling pathways related to the RPL/RPS gene signature of MBM [28]. Third, the only pathway that did not express any RPL/RPS genes was “Metabolism of carbohydrates” (highlighted in red). Fourth, the “Metabolism of carbohydrates” molecular pathway was found exclusively in omacetaxine-treated samples, with or without palbociclib (Figure 4A,C, right panels).
Furthermore, the mean value of each RPL/RPS gene in the MBM gene signature was evaluated. Notably, 4 of 21 RPL/RPS genes (RPL23, RPL35A, RPL6, and RPS18) were significantly downregulated in drug-treated groups when compared to placebo (Table 2). The fold-change difference in relation to placebo was quantified for each of the four members (Figure 5A). The crystal structure of the ribosome shows the protein configuration and highlights specific locations of RPL23, RPL35A, RPL6, and RPS18 proteins (Figure 5B). Strikingly, two of these four ribosomal proteins, RPL35A and RPL6, directly interact via hydrogen bonding, while the other two are found in opposite parts in the ribosome (Figure 5B). This could potentially indicate the important roles of these RPL/RPS proteins in ribosome assembly by being master regulators and recruiting other RPL/RPS members. Additionally, gene expression of proliferation markers was evaluated to determine palbociclib targets of proliferation (Table S5). MYC, MYCL, MYCN, PCNA, E2F3, E2F4, E2F7, and E2F8 proliferation biomarkers were downregulated in CTCs from drug-treated mice, when compared to placebo. Notably, the proto-oncogene c-MYC has been reported to promote translational potential by increasing ribosome biogenesis [31,63].

3.3. Metabolic Profiling of CTCs

It has been previously reported that metabolic alterations can be associated with ribosome defects and thus can promote the oncogenic potential of ribosomal proteins [64,65]. Given that our Reactome analyses showed an altered carbohydrate metabolism after omacetaxine treatment (Figure 4A,C), alongside substantial decreases in CTC numbers in CDX models relative to placebo and palbociclib (Figure 3 and Table 1), we undertook in vitro real-time metabolic profiling of patient-derived neoplastic blood fractions. These were complemented with omacetaxine and/or palbociclib treatment of CTC-derived clones prior to metabolic flux analysis.
Real-time metabolic flux analysis via mitochondrial stress tests on Lin− and Lin+ cell fractions from patient-derived blood revealed substantial heterogeneity between patients. In the case of SeahorseTM Patient 1, heterogeneity within the same patient with two months between blood collections was detected (longitudinal CTC monitoring) (Figure 6; patient parameters are listed in Table S4). Two-way ANOVA of healthy donors and patient Lin− and Lin+ blood fractions showed all comparisons as non-significant (Figure 6B). However, the two-way unpaired Mann Whitney test revealed significantly elevated respiration in the patient Lin− fraction relative to the healthy donor Lin− fraction (Figure S3A), an effect which persisted even after removal of Seahorse Patient 3′s exceptionally large and metabolically active Lin− fraction (Figure S3B). Seahorse Patient 3′s Lin− fraction exhibited significantly elevated basal respiration, basal ECAR, and ATP-coupled respiration relative to the patient’s own Lin+ fraction and one of the healthy donor. Though the healthy donor’s Lin− fraction yielded only enough cells for two replicate wells which precludes statistical analysis, it was unambiguously distinct from the patient’s Lin− fraction (Figure 6C), suggesting an exceptional property of the circulating neoplastic cell fraction at the time of collection.
Three-day drug treatment of CTC clones with omacetaxine and/or palbociclib prior to metabolic flux analysis revealed strikingly similar results to those seen in CDX models: omacetaxine treatment resulted in substantial decreases in cell numbers, relative to both palbociclib- and placebo-treated cells (Figure 3, Table 1 and Figure 7A right panel). Real-time metabolic flux analysis via mitochondrial stress test (MST) and Glycolytic Rate Assay (GRA) revealed that CTC clones have lower baseline mitochondrial respiration and higher basal glycolysis than weakly metastatic MeWos (Figure S4A,B). Additionally, CTC clones proved highly resistant to metabolic insult during flux analysis. They mounted far smaller compensatory responses to mitochondrial proton gradient decoupler FCCP during MSTs, and a smaller response to electron transport chain inhibition via Rot/AA during both GRAs and MSTs compared to weakly metastatic MeWos (Figure S4). However, the CTC-derived clone did respond robustly to the loss of glycolytic capacity in response to 2-DG injection (Figure 7C, left panel, and Figure S4, higher concentrations of mitochondrial proton gradient decoupler FCCP inhibitor not shown). Strikingly, cells retained a higher level of metabolic activity after palbociclib treatment than with omacetaxine, evidenced by higher per-cell OCR and ECAR after palbociclib treatment, and a decrease in both parameters after omacetaxine treatment (Figure 7 and Figure S4), bolstering Reactome analyses to show altered carbohydrate metabolism in this condition (Figure 4A,C).
Omacetaxine and/or palbociclib treatment alongside mitochondrial pyruvate transport inhibitor UK5099 revealed a greater degree of metabolic flexibility and competence remaining following palbociclib treatment. Only palbociclib-treated cells had excess pyruvate to direct toward aerobic glycolysis when mitochondrial transport was blocked by UK5099, suggesting that omacetaxine treatment left cells without spare fermentative capacity. While no treatment difference was significant in the presence of one-day treatment with glutaminase inhibitor BPTES, glutamine deprivation experiments alongside 3-day BPTES and inhibitor treatments revealed decreased cell counts in the absence of glutamine in both placebo- and palbociclib-treated cells, while omacetaxine treatment left cells largely without proliferative capacity in the presence or absence of glutamine (Figure S2).
From Reactome analyses of omacetaxine-treated cells, we selected key markers of carbohydrate metabolism that showed alteration relative to placebo-treated CDX mice and evaluated their expression at the mRNA level in CTC clones after inhibitor treatment by RT-qPCR. In omacetaxine-treated cells, two-way ANOVA revealed significant upregulation of GLS, GLUT1, and HIF1α, and significant downregulation of LDHB, PC, PFKFB4, and PGC1α relative to placebo-treated cells (Figure 7). Both omacetaxine and palbociclib showed substantially reduced oxidative capacity during glutamine deprivation, especially in the presence of 3-day treatment with BPTES (Figure S2A,C), in line with prior research demonstrating glutamine’s flexible role in maintaining oxidative metabolism to sustain proliferation [66].

4. Discussion

This work provides first-time evidence that examining the melanoma CTC RPL/RPS gene signature and metabolic changes in patients’ CTCs is a critical prognostic factor of melanoma metastasis in general and MBM in particular. These findings were achieved by performing inhibitory drug studies in vivo and in vitro, in addition to metabolic profiling of a drug-treated CTC clone and the first-ever report of metabolic flux analysis on patient neoplastic cell fractions in real-time. Our previous work has employed a multilevel analysis to determine a common CTC signature for RPL/RPS genes which is linked to MBM onset and progression [28]. In this study, we targeted ribosomal biogenesis by using the translational inhibitor omacetaxine with or without the CDK4/CDK6 inhibitor palbociclib. Our results indicate that dysregulated ribogenesis inhibits MBM and progression of secondary metastasis. The identification of 4 out of 21 RPL/RPS genes downregulated following treatment further refines our approach to evaluate and prevent melanoma metastasis. Deregulation of ribosomal protein expression can accompany metabolic alterations in promoting cancer metastasis and resistance to therapy [31]. We expanded on these concepts by evaluating them in the context of metabolic flux and carbohydrate metabolism. Downregulation of the four RP subunits following omacetaxine treatment co-occurred with substantially decreased cell viability and proliferation as well as a quiescent metabolic state in vitro. Conversely, downregulation of the same four units in responses to palbociclib treatment reduced cell proliferation but did not result in reduced metabolic competency at the bulk flux level.
Each drug and the combination of drugs in this study were carefully selected based on specific mechanisms of action, and the important interplays between protein synthesis and cell-cycle progression [30,31,34]. Omacetaxine prevents the initial elongation step of translation by binding to the A-site cleft in the large ribosomal subunit [35]. However, the exact mechanism of inhibition is unknown. Palbociclib selectively inhibits cycLin−dependent kinases CDK 4/6 [38]. To our knowledge, this is the first reported evidence that four specific RPL/RPS genes (RPL23, RPL35A, RPL6, and RPS18) are downregulated in response to translation and proliferation inhibition (Table 2). Moreover, the 3D reconstruction of the ribosome showed the exact localization of these four RPL/RPS proteins (Figure 5B). Interestingly, two out of four RPL/RPS proteins, specifically RPL6 and RPL35A, are in direct proximity to each other and physically interact (Figure 5B). It is intriguing to speculate that these four RPL/RPS members act as master regulators of ribosome biosynthesis and they recruit other RPL/RPS proteins to finalize ribosome assembly. This work sets the stage for further mechanistic studies to examine the function of each RPL/RPS member and their binding partners.
Numerous studies have reported multiple cellular functions of the RPL/RPS proteins other than ribosomal assembly [31,64,67,68]. RPL/RPS extra-ribosomal functions include cell differentiation, cellular redox homeostasis, tumor suppression, cancer progression, immune signaling, p53 activation, oncogenic transformation, abnormal maturation of tumor cells, and DNA repair [64,68]. Changes in RPL/RPS expression can alter transcriptional-specific translational output, considering that upregulation of ribosomal proteins has been observed in cancer [69]. We have previously reported that the 21 RPL/RPS genes of the CTC signature are associated with melanoma brain metastasis development and progression [28].
Ribosome composition has been reported to be highly heterogeneous and dynamic [67], and different cell types have different ribosome composition due to having different binding partners [67,68,70]. Notably, deregulation of individual ribosomal proteins have specific patterns across cancer subtypes that can predict clinical outcomes [71,72] and specific somatic mutations in the RPL/RPS genes have been observed in various cancer types, including melanoma [64]. Importantly, oncogene-dependent cancer types have been shown to have enhanced activation of protein synthesis associated with metabolic rewiring of cells [73,74], and high ribosome biosynthesis is vital for proliferating and metastasizing cancer cells [75]. Disruption of ribosome biosynthesis results in ribosomopathies [64,76]. For example, deregulation of ribosomal protein synthesis via RPL15, and in turn protein translation, have been reported to be associated with breast cancer metastasis [30]. Alterations of RPL/RPS composition in cancer cells must be tightly regulated, since the imbalanced expression of RPL/RPS proteins can initiate the p53 response, leading to cell cycle arrest and apoptosis [64,68,77].
The diversity of RPL/RPS members and their specificity raise three intriguing possibilities for specific RPL/RPS contributions reducing the onset of MBM/extracranial metastasis: (1) the inhibition of the four RPL/RPS members (RPL22, RPL35A, RPL6, and RPS18) can alter the composition of ribosomes which may become specialized for the translation of specific oncogenic transcripts, similar to what has been observed in embryonic stem cells [78,79]; (2) the presence and action of extra-ribosomal RPL/RPS functions, e.g., the activation of p53-dependent/independent pathways in response to stress, resulting in cell cycle arrest and/or apoptosis [64]; (3) the involvement of RPL/RPS members in the translational control of circulating immune cells, affecting the development of metastasis. Examples are CTCs interacting with T cells where ribosomes lacking some RPs influence CD8+ T cell immunosurveillance [80], and neutrophils or myeloid-derived suppressor cells which are known to generate clusters with CTCs and whose translational capacity regulates the CTC metastatic propensities [25,81]. Collectively, these findings underscore the diversity of roles of RPL/RPS members in a diverse range of settings, connecting their canonical (“Ribogenesis”) and non-canonical functions (“Ribosomopathy”) with cancer.
Multiple oncogenes also known to modulate metabolic function, in particular MYC, have been linked to increased ribosomal synthesis and translation [29,63,82]. Our results demonstrate that targeting either cell ribogenesis machinery or cell proliferation results in a significant decrease in MYC and MYCL gene expression, along with gene expression of other proliferative biomarkers, e.g., PCNA, E2F3, E2F4, and E2F8 (Table S5). MYC controls transcription of many biomarkers, including RPL/RPS proteins [29,83], and c-MYC has been shown to promote expression of translational machinery, and ribosome and tRNA biogenesis directly [84]. Interestingly, alongside loss of c-MYC activity after inhibition of ribogenesis in the presence of both omacetaxine and palbociclib, only omacetaxine-treated CTCs showed upregulation of HIF1α at the mRNA level (Figure 7B). HIF1α and c-MYC exert overlapping actions to facilitate metabolic flexibility during metastatic progression, both complementary and antagonistic [84]. While both c-MYC and HIFα signaling can promote anabolic function directly and indirectly, HIFα does so by preventing glycolytic end-products from entering the mitochondria [84]. This elevation of HIF1α mRNA expression may represent a failed compensatory response to the loss of c-MYC’s anabolic functionality and overall anabolic incompetence during impaired ribogenesis, but not cell-cycle inhibition, as evidenced by the decreased basal glycolysis in omacetaxine but not palbociclib-treated cells. It is possible that upregulated mRNA expression of GLS and GLUT1 represent part of this same failed compensatory response (Figure 7B).
Our results build on published studies since we show the link between proliferative inhibition and downregulation of RPL/RPS gene expression (Table S5). Additionally, we show that the deregulation of carbohydrate metabolism that occurs alongside impaired ribogenesis represents a possible therapeutic vulnerability window of CTCs. Metabolic flux analysis on the brain metastatic CTC-derived clone showed that omacetaxine treatment substantially reduced its proliferation and viability, and pushed these cells toward a quiescent metabolic state (Figure 7A and Figure S2C). While palbociclib also reduced metastasis in our CDX models, it caused elevation of both basal OCR and basal ECAR, and resulted in a smaller decrease in total CTC numbers in vitro and in vivo than did omacetaxine (Figure 3, Table 1, Figure 7A right panel, and Figure S2C). It is possible that the anti-metastatic effects of palbociclib in vivo are independent of bulk metabolic flux alteration, or that altered patterns of interaction with CDX host cells and/or tissue microenvironments are required for this effect.
Given our results showing that omacetaxine treatment co-occurred with decreased metastasis and impaired carbohydrate metabolism in the CTC-derived clone, we suggest a high level of anabolic activity may occur in the presence or absence of elevated oxidative activity, and may represent a trait of metastasis-competent CTCs during dissemination. It is important to note that, due to the heterogeneity between tumors, CTCs derived from them, and the properties of derived metastatic tumors, the relatively lower level of oxidative activity in the MBM-competent CTC clone is independent of the High-OXPHOS and Low-OXPHOS MBM paradigm [40]. Taken together, metabolic flux and RT-qPCR analyses show that omacetaxine treatment of CTCs results in a decrease in metabolic flexibility relative to palbociclib, specifically the inability to re-direct pyruvate as a substrate to compensate for metabolic insult. Inhibitor treatment alongside mitochondrial pyruvate carrier (MPC) inhibitor UK5099 revealed a lack of excess pyruvate after omacetaxine treatment to re-direct from mitochondrial transport to glycolytic flux (Figure 7C). Downregulation of PC (pyruvate carboxylase), which yields oxaloacetate from pyruvate, suggests a loss of anaplerotic function after omacetaxine treatment.
Additionally, we found downregulation of LDHB but not LDHA after omacetaxine treatment (Figure 7B). LDHA preferentially converts pyruvate to lactate, while LDHB preferentially performs the reverse reaction, suggesting further loss of flexibility in using pyruvate as a substrate during impaired ribogenesis. PFKFB4, which confers anabolic flexibility by preferentially re-directing glucose to the pentose phosphate pathway [85,86] and enhances melanoma cell migration through RAS/AKT-dependent signaling independent of glycolysis [87], was also downregulated at the mRNA level after omacetaxine but not palbociclib treatment (Figure 7B).
Notably, we report the first-ever metabolic flux analysis in patient-derived CTCs in real-time. While there was considerable variation between and within patients, we saw overall elevated oxidative activity in patient Lin− compared to healthy donor Lin− samples (Figure S3A), which persisted even after excluding the exceptional Lin− fraction of Seahorse Patient 3 (Figure S3B). Seahorse Patient 3′s neoplastic fraction showed substantial elevation of both OCR and ECAR relative to comparable fractions (Figure 6C), and the extraordinary character of this patient’s Lin− fraction may reveal a rare snapshot of a particularly malignant CTC-shedding event.
The size of obtainable patient Lin− fractions may often preclude direct drug testing on patient-derived neoplastic cells during metabolic flux analysis, so examining direct effects of targeted inhibition on patient CTCs remains elusive. However, our findings of elevated oxidative metabolism in patient Lin− fractions and substantial elevation of metabolic activity (Seahorse Patient 3) highlight the interface of oxidative and anabolic metabolism as a potential therapeutic vulnerability during CTC dissemination. PC is a key mediator of anaplerosis in cancer cells and a central node for metabolic plasticity during cancer progression as it can promote both catabolic and anabolic metabolism [88]. PC is more highly expressed in cancer cells than in healthy cells; its overexpression is linked with the increased invasiveness of cancer cells and targeting PC with small molecule inhibitors has shown promise in suppressing cancer progression [89]. Additionally, inhibition of PC, along with mitochondrial pyruvate carriers (MPCs), has shown promise in sensitizing Tyrosine Kinase Inhibitor (TKI)-resistant leukemic stem cells to treatment [90]. We found that omacetaxine-treated CDXs show decreased CTC numbers and reduced MBM/extracranial metastasis, likely without the metabolic flexibility conferred by PC expression. Altogether, our findings suggest that disrupting PC function may be critical in clinical interventions to suppress and/or prevent metastasis.
Our results are altogether in line with previous findings showing heterogeneity in the oxidative capacity of neoplastic cells and the flexible use of carbohydrate metabolism to expand available metastatic niches [39,40,41,42,43,44]. We expand the current knowledge of CTC biology by demonstrating elevated oxidative metabolism in patient-derived Lin−/CTC populations relative to that of healthy donor blood. Employing the MBM CTC-derived clone, we demonstrate that impaired anabolic competency during dysregulated ribogenesis reveals a vulnerability for MBM-competent CTCs during dissemination. Targeting pyruvate metabolism can represent a potent and broadly viable therapeutic avenue for impairing CTC metastatic capacity, particularly in drug-resistant cancers. Future research characterizing PC expression in patient Lin−/CTCs at the mRNA/protein levels, in combination with PC expression assessment during glutamine deprivation, MPC inhibition, and metabolic flux analysis, will further elucidate how to leverage this CTC vulnerability for therapeutic efficacy.
We postulate that translationally activated, metastasis-competent CTCs could be novel therapeutic targets to combat cancer. Nevertheless, this study has some limitations. First, a limited number of melanoma patients was used in these analyses; therefore, we cannot claim that these observations are applicable to all melanoma patients. Heterogeneity between and within patients adds an additional level of complexity in the way of reaching definite conclusions, and further studies are warranted. Second, the drug inhibitory investigations used a limited number of CDX mice, which, although statistically significant, may introduce inherent sampling bias. Third, because we did not execute an inhibitory study in humanized mice, we cannot exclude the possibility that the drugs could affect interactions of CTCs with circulatory immune cells that could alter the immune response. Fourth, while our CTC enrichment method removes RBCs, PBMCs, granulocytes, and platelets upstream of real-time metabolic flux analysis, our current knowledge of CTC biology suggests that it is unlikely that the entirety of collected Lin−/CTC fractions are all CTCs per se. Nevertheless, marker-agnostic approaches to neoplastic cell characterization are necessary to fully capture the entire spectrum of metastasis-competent CTCs [43,91]. Our analyses indicated a trend toward greater size of patient Lin−/CTC fractions (Figure S3C), as well as discrete metabolic characteristics of sorted Lin−/CTC fractions from patients relative to those from healthy donor blood (Figure 6C, Figure S3A,B). We therefore suggest that our analyses reveal unique metabolic properties of melanoma patient Lin−/CTC fractions. Fifth, bulk flux analysis of CTCs after drug treatment may miss subtler shifts in metabolic state. While our targeted, Reactome-driven examination of the critical dimensions of carbohydrate metabolism demonstrated the power of this approach for rigorously characterizing drug responses, future functional assays alongside proteomic and metabolomic approaches will more fully characterize CTC metabolic vulnerabilities during targeted therapeutic disruption of the onco-ribosome, as well as revealing further heterogeneity of CTC metabolism [42,43,91,92]. Palbociclib’s anti-metastatic effects in vivo show that the loss of metabolic flexibility accompanying impaired carbohydrate metabolism downstream of disrupted ribogenesis represents only one possible pathway to preventing metastasis. However, omacetaxine yielded greater cell death, and decreased metabolic competency and flexibility in vitro, as well as fewer CTCs in vivo, than did treatment with palbociclib, which may be critical for complete elimination of rare subpopulations of metastasis-competent CTC before and/or during dissemination.

5. Conclusions

This study represents pioneering work that identifies the relevance of targeting the CTC RPL/RPS signature to suppress melanoma metastasis in general and MBM in particular. The systemic examination of patients’ RPL/RPS gene expression in CTCs could be pivotal to prescribing more targeted treatment based on the needs of the patient. Future experimental and analytical innovations built on our findings will aid in novel therapeutic interventions and improvement of clinical outcomes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15215263/s1, Figure S1: Assessment of mouse metastatic patterns by IVIS imaging; Figure S2: Glutamine deprivation experiments in the presence/absence of glutamine and glutaminase inhibitor BPTES (10 µM); Figure S3: Mitochondrial respiration; Figure S4: Baseline metabolic flux results in the melanoma CTC-derived clone; Table S1: Statistical analyses for IVIS flux quantification of melanoma brain metastasis; Table S2: Statistical analyses for IVIS flux quantification of extracranial metastasis; Table S3: Top 12 downregulated genes in MBM; Table S4: Clinico-pathological characteristics of melanoma patients used in this study; Table S5: Individual and mean values (in cpm) of proliferation biomarkers.

Author Contributions

Conceptualization, T.Y.B., S.D.M. and D.M.; methodology, B.T.; performance of experiments, T.Y.B., S.D.M., I.V.L., M.C.O. and M.L.; validation, T.Y.B., S.D.M., B.T. and D.M; investigation, T.Y.B., S.D.M. and D.M.; B.T. provided blood samples, patients’ parameters, and clinical expertise; writing—original draft preparation, T.Y.B., S.D.M. and D.M.; writing—review and editing, T.Y.B., S.D.M. and D.M.; visualization, T.Y.B., S.D.M., M.C.O. and D.M.; supervision, D.M.; project administration, T.Y.B., S.D.M. and D.M.; funding acquisition, D.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NIH grant 1R01 CA216991, the UNM-CCC Translational Initiative (TSI), the Oxnard Foundation (to D.M.), by the NCI P30CA118100-16 grant supporting the UNM Comprehensive Cancer Center.

Institutional Review Board Statement

The Institutional Review Board of the University of New Mexico Health Sciences Center approved this study, per protocols (HRRC #15-507 and HRCC #19-513).

Informed Consent Statement

Informed consent was obtained from all melanoma patients involved in the study, per IRB-approved protocols (HRRC #15-507 and HRCC #19-513).

Data Availability Statement

Not applicable.

Acknowledgments

We are deeply appreciative of all melanoma patients who consented to provide blood samples for this study. We gratefully acknowledge Mara Steinkamp and Lillian Fitzpat-rick (UNM-CCC Animal Models Shared Resources) for their extensive help with in vivo work throughout the study. We extend our gratitude to Sharina Desai and the AIM Core for assistance with metabolic studies, to Hannah Taylor-Noren for PyMOL assistance, and to Isaac Lira for consenting melanoma patients to donate blood for these studies. In addition, we thank Olufunmilola Oyebamiji for providing bioinformatics assistance, and appreciate the assistance provided by ATG, HTR, Biostatistics, and Fluorescent Microscopy Shared Resources of UNM-CCC.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. In, G.K.; Poorman, K.; Saul, M.; O’Day, S.; Farma, J.M.; Olszanski, A.J.; Gordon, M.S.; Thomas, J.S.; Eisenberg, B.; Flaherty, L. Molecular profiling of melanoma brain metastases compared to primary cutaneous melanoma and to extracranial metastases. Oncotarget 2020, 11, 3118. [Google Scholar] [CrossRef]
  2. Eroglu, Z.; Holmen, S.L.; Chen, Q.; Khushalani, N.I.; Amaravadi, R.; Thomas, R.; Ahmed, K.A.; Tawbi, H.; Chandra, S.; Markowitz, J. Melanoma central nervous system metastases: An update to approaches, challenges, and opportunities. Pigment Cell Melanoma Res. 2019, 32, 458–469. [Google Scholar] [CrossRef]
  3. Cohen, J.V.; Tawbi, H.; Margolin, K.A.; Amravadi, R.; Bosenberg, M.; Brastianos, P.K.; Chiang, V.L.; de Groot, J.; Glitza, I.C.; Herlyn, M.J.P.c. Melanoma central nervous system metastases: Current approaches, challenges, and opportunities. Pigment Cell Melanoma Res. 2016, 29, 627–642. [Google Scholar] [CrossRef]
  4. Zakrzewski, J.; Geraghty, L.N.; Rose, A.E.; Christos, P.J.; Mazumdar, M.; Polsky, D.; Shapiro, R.; Berman, R.; Darvishian, F.; Hernando, E. Clinical variables and primary tumor characteristics predictive of the development of melanoma brain metastases and post-brain metastases survival. Cancer 2011, 117, 1711–1720. [Google Scholar] [CrossRef] [PubMed]
  5. Davies, M.A.; Liu, P.; McIntyre, S.; Kim, K.B.; Papadopoulos, N.; Hwu, W.J.; Hwu, P.; Bedikian, A. Prognostic factors for survival in melanoma patients with brain metastases. Cancer 2011, 117, 1687–1696. [Google Scholar] [CrossRef] [PubMed]
  6. Sampson, J.H.; Carter, J.H.; Friedman, A.H.; Seigler, H.F. Demographics, prognosis, and therapy in 702 patients with brain metastases from malignant melanoma. J. Neurosurg. 1998, 88, 11–20. [Google Scholar] [CrossRef] [PubMed]
  7. Luke, J.J.; Flaherty, K.T.; Ribas, A.; Long, G.V. Targeted agents and immunotherapies: Optimizing outcomes in melanoma. Nat. Rev. Clin. Oncol. 2017, 14, 463–482. [Google Scholar] [CrossRef]
  8. Fife, K.; Colman, M.; Stevens, G.; Firth, I.; Moon, D.; Shannon, K.; Harman, R.; Petersen-Schaefer, K.; Zacest, A.; Besser, M. Determinants of outcome in melanoma patients with cerebral metastases. J. Clin. Oncol. 2004, 22, 1293–1300. [Google Scholar] [CrossRef]
  9. Berghoff, A.S.; Schur, S.; Füreder, L.M.; Gatterbauer, B.; Dieckmann, K.; Widhalm, G.; Hainfellner, J.; Zielinski, C.C.; Birner, P.; Bartsch, R. Descriptive statistical analysis of a real life cohort of 2419 patients with brain metastases of solid cancers. ESMO Open 2016, 1, e000024. [Google Scholar] [CrossRef]
  10. Zoga, E.; Wolff, R.; Ackermann, H.; Meissner, M.; Rödel, C.; Tselis, N.; Chatzikonstantinou, G. Factors associated with hemorrhage of melanoma brain metastases after stereotactic radiosurgery in the era of targeted/immune checkpoint inhibitor therapies. Cancers 2022, 14, 2391. [Google Scholar] [CrossRef]
  11. Gupta, G.P.; Massagué, J. Cancer metastasis: Building a framework. Cell 2006, 127, 679–695. [Google Scholar] [CrossRef]
  12. Giuliano, M.; Shaikh, A.; Lo, H.C.; Arpino, G.; De Placido, S.; Zhang, X.H.; Cristofanilli, M.; Schiff, R.; Trivedi, M.V. Perspective on circulating tumor cell clusters: Why it takes a village to metastasize. Cancer Res. 2018, 78, 845–852. [Google Scholar] [CrossRef]
  13. Vishnoi, M.; Boral, D.; Marchetti, D. Transcriptomic Analysis of CTCs. In Circulating Tumor Cells: Advances in Liquid Biopsy Technologies; Springer: Berlin/Heidelberg, Germany, 2023; pp. 277–299. [Google Scholar]
  14. Ildiz, E.S.; Gvozdenovic, A.; Kovacs, W.J.; Aceto, N. Travelling under pressure—Hypoxia and shear stress in the metastatic journey. Clin. Exp. Metastasis 2023, 40, 375–394. [Google Scholar] [CrossRef] [PubMed]
  15. Godet, I.; Shin, Y.J.; Ju, J.A.; Ye, I.C.; Wang, G.; Gilkes, D.M. Fate-mapping post-hypoxic tumor cells reveals a ROS-resistant phenotype that promotes metastasis. Nat. Commun. 2019, 10, 4862. [Google Scholar] [CrossRef] [PubMed]
  16. Werner-Klein, M.; Scheitler, S.; Hoffmann, M.; Hodak, I.; Dietz, K.; Lehnert, P.; Naimer, V.; Polzer, B.; Treitschke, S.; Werno, C. Genetic alterations driving metastatic colony formation are acquired outside of the primary tumour in melanoma. Nat. Commun. 2018, 9, 1–17. [Google Scholar] [CrossRef] [PubMed]
  17. Klein, C.A. Selection and adaptation during metastatic cancer progression. Nature 2013, 501, 365–372. [Google Scholar] [CrossRef]
  18. Gomis, R.R.; Gawrzak, S. Tumor cell dormancy. Mol. Oncol. 2017, 11, 62–78. [Google Scholar] [CrossRef]
  19. Marrella, A.; Fedi, A.; Varani, G.; Vaccari, I.; Fato, M.; Firpo, G.; Guida, P.; Aceto, N.; Scaglione, S. High blood flow shear stress values are associated with circulating tumor cells cluster disaggregation in a multi-channel microfluidic device. PLoS ONE 2021, 16, e0245536. [Google Scholar] [CrossRef]
  20. Dianat-Moghadam, H.; Azizi, M.; Eslami-S, Z.; Cortés-Hernández, L.E.; Heidarifard, M.; Nouri, M.; Alix-Panabières, C. The role of circulating tumor cells in the metastatic cascade: Biology, technical challenges, and clinical relevance. Cancers 2020, 12, 867. [Google Scholar] [CrossRef]
  21. Luzzi, K.J.; MacDonald, I.C.; Schmidt, E.E.; Kerkvliet, N.; Morris, V.L.; Chambers, A.F.; Groom, A.C. Multistep nature of metastatic inefficiency: Dormancy of solitary cells after successful extravasation and limited survival of early micrometastases. Am. J. Pathol. 1998, 153, 865–873. [Google Scholar] [CrossRef]
  22. Pantel, K.; Alix-Panabières, C.; Riethdorf, S. Cancer micrometastases. Nat. Rev. Clin. Oncol. 2009, 6, 339–351. [Google Scholar] [CrossRef] [PubMed]
  23. Lucci, A.; Hall, C.S.; Patel, S.P.; Narendran, B.; Bauldry, J.B.; Royal, R.E.; Karhade, M.; Upshaw, J.R.; Wargo, J.A.; Glitza, I.C.J.C.C.R. Circulating tumor cells and early relapse in node-positive melanoma. Clin. Cancer Res. 2020, 26, 1886–1895. [Google Scholar] [CrossRef] [PubMed]
  24. Aceto, N.; Bardia, A.; Miyamoto, D.T.; Donaldson, M.C.; Wittner, B.S.; Spencer, J.A.; Yu, M.; Pely, A.; Engstrom, A.; Zhu, H. Circulating tumor cell clusters are oligoclonal precursors of breast cancer metastasis. Cell 2014, 158, 1110–1122. [Google Scholar] [CrossRef]
  25. Sprouse, M.L.; Welte, T.; Boral, D.; Liu, H.N.; Yin, W.; Vishnoi, M.; Goswami-Sewell, D.; Li, L.; Pei, G.; Jia, P. PMN-MDSCs enhance CTC metastatic properties through reciprocal interactions via ROS/Notch/Nodal signaling. Int. J. Mol. Sci. 2019, 20, 1916. [Google Scholar] [CrossRef]
  26. Donato, C.; Kunz, L.; Castro-Giner, F.; Paasinen-Sohns, A.; Strittmatter, K.; Szczerba, B.M.; Scherrer, R.; Di Maggio, N.; Heusermann, W.; Biehlmaier, O. Hypoxia triggers the intravasation of clustered circulating tumor cells. Cell Rep. 2020, 32, 108105. [Google Scholar] [CrossRef]
  27. Michor, F.; Liphardt, J.; Ferrari, M.; Widom, J. What does physics have to do with cancer? Nat. Rev. Cancer 2011, 11, 657–670. [Google Scholar] [CrossRef] [PubMed]
  28. Bowley, T.Y.; Lagutina, I.V.; Francis, C.; Sivakumar, S.; Selwyn, R.G.; Taylor, E.; Guo, Y.; Fahy, B.N.; Tawfik, B.; Marchetti, D. The RPL/RPS gene signature of melanoma CTCs associates with brain metastasis. Cancer Res. Commun. 2022, 2, 1436–1448. [Google Scholar] [CrossRef] [PubMed]
  29. Micalizzi, D.S.; Ebright, R.Y.; Haber, D.A.; Maheswaran, S. Translational regulation of cancer metastasis. Cancer Res. 2021, 81, 517–524. [Google Scholar] [CrossRef]
  30. Ebright, R.Y.; Lee, S.; Wittner, B.S.; Niederhoffer, K.L.; Nicholson, B.T.; Bardia, A.; Truesdell, S.; Wiley, D.F.; Wesley, B.; Li, S. Deregulation of ribosomal protein expression and translation promotes breast cancer metastasis. Science 2020, 367, 1468–1473. [Google Scholar] [CrossRef]
  31. Elhamamsy, A.R.; Metge, B.J.; Alsheikh, H.A.; Shevde, L.A.; Samant, R.S. Ribosome biogenesis: A central player in cancer metastasis and therapeutic resistance. Cancer Res. 2022, 82, 2344–2353. [Google Scholar] [CrossRef]
  32. Robichaud, N.; Sonenberg, N.; Ruggero, D.; Schneider, R. Translational control in cancer. Cold Spring Harb. Perspect. Biol. 2019, 11, a032896. [Google Scholar] [CrossRef] [PubMed]
  33. Zhao, W.; Li, X.; Nian, W.; Wang, J.; Wang, X.; Sun, L.; Zhu, Y.; Tong, Z. Ribosome proteins represented by RPL27A mark the development and metastasis of triple-negative breast cancer in mouse and human. Front. Cell Dev. Biol. 2021, 9, 716730. [Google Scholar] [CrossRef]
  34. Volarević, S.a.; Stewart, M.J.; Ledermann, B.; Zilberman, F.; Terracciano, L.; Montini, E.; Grompe, M.; Kozma, S.C.; Thomas, G. Proliferation, but not growth, blocked by conditional deletion of 40 S ribosomal protein S6. Science 2000, 288, 2045–2047. [Google Scholar] [CrossRef]
  35. Gandhi, V.; Plunkett, W.; Cortes, J.E. Omacetaxine: A protein translation inhibitor for treatment of chronic myelogenous leukemia. Clin. Cancer Res. 2014, 20, 1735–1740. [Google Scholar] [CrossRef] [PubMed]
  36. Wetzler, M.; Segal, D. Omacetaxine as an anticancer therapeutic: What is old is new again. Curr. Pharm. Des. 2011, 17, 59–64. [Google Scholar] [CrossRef] [PubMed]
  37. Rosshandler, Y.; Shen, A.Q.; Cortes, J.; Khoury, H.J. Omacetaxine mepesuccinate for chronic myeloid leukemia. Expert Rev. Hematol. 2016, 9, 419–424. [Google Scholar] [CrossRef]
  38. Turner, N.C.; Ro, J.; André, F.; Loi, S.; Verma, S.; Iwata, H.; Harbeck, N.; Loibl, S.; Huang Bartlett, C.; Zhang, K. Palbociclib in hormone-receptor–positive advanced breast cancer. N. Engl. J. Med. 2015, 373, 209–219. [Google Scholar] [CrossRef]
  39. Scott, D.A.; Richardson, A.D.; Filipp, F.V.; Knutzen, C.A.; Chiang, G.G.; Ze’ev, A.R.; Osterman, A.L.; Smith, J.W. Comparative metabolic flux profiling of melanoma cell lines: Beyond the Warburg effect. J. Biol. Chem. 2011, 286, 42626–42634. [Google Scholar] [CrossRef]
  40. Fischer, G.M.; Guerrieri, R.A.; Hu, Q.; Joon, A.Y.; Kumar, S.; Haydu, L.E.; McQuade, J.L.; Vashisht Gopal, Y.; Knighton, B.; Deng, W. Clinical, molecular, metabolic, and immune features associated with oxidative phosphorylation in melanoma brain metastases. Neuro-Oncol. Adv. 2021, 3, vdaa177. [Google Scholar] [CrossRef]
  41. Bulfoni, M.; Turetta, M.; Del Ben, F.; Di Loreto, C.; Beltrami, A.P.; Cesselli, D. Dissecting the heterogeneity of circulating tumor cells in metastatic breast cancer: Going far beyond the needle in the haystack. Int. J. Mol. Sci. 2016, 17, 1775. [Google Scholar] [CrossRef]
  42. Alix-Panabieres, C.; Cayrefourcq, L.; Mazard, T.; Maudelonde, T.; Assenat, E.; Assou, S. Molecular portrait of metastasis-competent circulating tumor cells in colon cancer reveals the crucial role of genes regulating energy metabolism and DNA repair. Clin. Chem. 2017, 63, 700–713. [Google Scholar] [CrossRef] [PubMed]
  43. Micalizzi, D.S.; Maheswaran, S.; Haber, D.A. A conduit to metastasis: Circulating tumor cell biology. Genes Dev. 2017, 31, 1827–1840. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, W.; Xu, F.; Yao, J.; Mao, C.; Zhu, M.; Qian, M.; Hu, J.; Zhong, H.; Zhou, J.; Shi, X.J. Single-cell metabolic fingerprints discover a cluster of circulating tumor cells with distinct metastatic potential. Nat. Commun. 2023, 14, 2485. [Google Scholar] [CrossRef] [PubMed]
  45. Stocking, K.L.; Jones, J.C.; Everds, N.E.; Buetow, B.S.; Roudier, M.P.; Miller, R.E. Use of Low-Molecular–Weight Heparin to Decrease Mortality in Mice after Intracardiac Injection of Tumor Cells. Comp. Med. 2009, 59, 37–45. [Google Scholar]
  46. Devlin, J.R.; Hannan, K.M.; Hein, N.; Cullinane, C.; Kusnadi, E.; Ng, P.Y.; George, A.J.; Shortt, J.; Bywater, M.J.; Poortinga, G. Combination therapy targeting ribosome biogenesis and mRNA translation synergistically extends survival in MYC-driven lymphoma. Cancer Discov. 2016, 6, 59–70. [Google Scholar] [CrossRef]
  47. Slominski, A.; Tobin, D.J.; Shibahara, S.; Wortsman, J. Melanin pigmentation in mammalian skin and its hormonal regulation. Physiol. Rev. 2004, 84, 1155–1228. [Google Scholar] [CrossRef]
  48. Frerich, C.A.; Brayer, K.J.; Painter, B.M.; Kang, H.; Mitani, Y.; El-Naggar, A.K.; Ness, S.A. Transcriptomes define distinct subgroups of salivary gland adenoid cystic carcinoma with different driver mutations and outcomes. Oncotarget 2018, 9, 7341. [Google Scholar] [CrossRef]
  49. Brown, R.B.; Madrid, N.J.; Suzuki, H.; Ness, S.A. Optimized approach for Ion Proton RNA sequencing reveals details of RNA splicing and editing features of the transcriptome. PLoS ONE 2017, 12, e0176675. [Google Scholar] [CrossRef]
  50. Pauken, C.M.; Kenney, S.R.; Brayer, K.J.; Guo, Y.; Brown-Glaberman, U.A.; Marchetti, D. Heterogeneity of Circulating Tumor Cell Neoplastic Subpopulations Outlined by Single-Cell Transcriptomics. Cancers 2021, 13, 4885. [Google Scholar] [CrossRef]
  51. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef]
  52. Alexa, A.; Rahnenfuhrer, J. topGO: Enrichment Analysis for Gene Ontology. R Package Version 2.28.0. Cranio 2016. Package. Available online: https://bioconductor.org/packages/release/bioc/html/topGO.html (accessed on 5 March 2023).
  53. Croft, D.; O’Kelly, G.; Wu, G.; Haw, R.; Gillespie, M.; Matthews, L.; Caudy, M.; Garapati, P.; Gopinath, G.; Jassal, B. Reactome: A database of reactions, pathways, and biological processes. Nucleic Acids Res. 2011, 39, D691–D697. [Google Scholar] [CrossRef]
  54. Lill, M.A.; Danielson, M.L. Computer-aided drug design platform using PyMOL. J. Comput.-Aided Mol. Des. 2011, 25, 13–19. [Google Scholar] [CrossRef]
  55. Khatter, H.; Myasnikov, A.G.; Natchiar, S.K.; Klaholz, B.P. Structure of the human 80S ribosome. Nature 2015, 520, 640–645. [Google Scholar] [CrossRef] [PubMed]
  56. Vishnoi, M.; Boral, D.; Liu, H.; Sprouse, M.L.; Yin, W.; Goswami-Sewell, D.; Tetzlaff, M.T.; Davies, M.A.; Oliva, I.C.G.; Marchetti, D. Targeting USP7 identifies a metastasis-competent state within bone marrow–resident melanoma CTCs. Cancer Res. 2018, 78, 5349–5362. [Google Scholar] [CrossRef]
  57. Lee, A.S.; Wu, J.C. Imaging of embryonic stem cell migration in vivo. In Stem Cell Migration; Springer: Berlin/Heidelberg, Germany, 2011; pp. 101–114. [Google Scholar]
  58. Tessema, B.; Riemer, J.; Sack, U.; König, B. Cellular Stress Assay in Peripheral Blood Mononuclear Cells: Factors Influencing Its Results. Int. J. Mol. Sci. 2022, 23, 13118. [Google Scholar] [CrossRef]
  59. König, B.; Lahodny, J. Ozone high dose therapy (OHT) improves mitochondrial bioenergetics in peripheral blood mononuclear cells. Transl. Med. Commun. 2022, 7, 1–12. [Google Scholar] [CrossRef] [PubMed]
  60. De Vitis, C.; Capalbo, C.; Torsello, A.; Napoli, C.; Salvati, V.; Loffredo, C.; Blandino, G.; Piaggio, G.; Auciello, F.R.; Pelliccia, F. Opposite Effect of Thyroid Hormones on Oxidative Stress and on Mitochondrial Respiration in COVID-19 Patients. Antioxidants 2022, 11, 1998. [Google Scholar] [CrossRef] [PubMed]
  61. Amintas, S.; Bedel, A.; Moreau-Gaudry, F.; Boutin, J.; Buscail, L.; Merlio, J.-P.; Vendrely, V.; Dabernat, S.; Buscail, E. Circulating tumor cell clusters: United we stand divided we fall. Int. J. Mol. Sci. 2020, 21, 2653. [Google Scholar] [CrossRef]
  62. Hamilton, E.; Infante, J.R. Targeting CDK4/6 in patients with cancer. Cancer Treat. Rev. 2016, 45, 129–138. [Google Scholar] [CrossRef]
  63. Van Riggelen, J.; Yetil, A.; Felsher, D.W. MYC as a regulator of ribosome biogenesis and protein synthesis. Nat. Rev. Cancer 2010, 10, 301–309. [Google Scholar] [CrossRef]
  64. Kang, J.; Brajanovski, N.; Chan, K.T.; Xuan, J.; Pearson, R.B.; Sanij, E. Ribosomal proteins and human diseases: Molecular mechanisms and targeted therapy. Signal Transduct. Target. Ther. 2021, 6, 323. [Google Scholar] [CrossRef] [PubMed]
  65. Sulima, S.O.; Kampen, K.R.; De Keersmaecker, K. Cancer Biogenesis in Ribosomopathies. Cells 2019, 8, 229. [Google Scholar] [CrossRef] [PubMed]
  66. Filipp, F.V.; Ratnikov, B.; De Ingeniis, J.; Smith, J.W.; Osterman, A.L.; Scott, D.A. Glutamine-fueled mitochondrial metabolism is decoupled from glycolysis in melanoma. Pigment Cell Melanoma Res. 2012, 25, 732–739. [Google Scholar] [CrossRef] [PubMed]
  67. Genuth, N.R.; Barna, M. The discovery of ribosome heterogeneity and its implications for gene regulation and organismal life. Mol. Cell 2018, 71, 364–374. [Google Scholar] [CrossRef]
  68. Guimaraes, J.C.; Zavolan, M. Patterns of ribosomal protein expression specify normal and malignant human cells. Genome Biol. 2016, 17, 1–13. [Google Scholar] [CrossRef]
  69. de Las Heras-Rubio, A.; Perucho, L.; Paciucci, R.; Vilardell, J.; ME, L.L. Ribosomal proteins as novel players in tumorigenesis. Cancer Metastasis Rev. 2014, 33, 115–141. [Google Scholar] [CrossRef]
  70. Bortoluzzi, S.; d’Alessi, F.; Romualdi, C.; Danieli, G.A. Differential expression of genes coding for ribosomal proteins in different human tissues. Bioinformatics 2001, 17, 1152–1157. [Google Scholar] [CrossRef]
  71. Henry, J.L.; Coggin, D.L.; King, C.R. High-Level Expression of the Ribosomal Protein L19 in Human Breast Tumors That Overexpress erb B-2. Cancer Res. 1993, 53, 1403–1408. [Google Scholar]
  72. Takagi, M.; Absalon, M.J.; McLure, K.G.; Kastan, M.B. Regulation of p53 translation and induction after DNA damage by ribosomal protein L26 and nucleolin. Cell 2005, 123, 49–63. [Google Scholar] [CrossRef]
  73. Pelletier, J.; Thomas, G.; Volarević, S. Ribosome biogenesis in cancer: New players and therapeutic avenues. Nat. Rev. Cancer 2018, 18, 51–63. [Google Scholar] [CrossRef]
  74. Penzo, M.; Montanaro, L.; Treré, D.; Derenzini, M. The Ribosome Biogenesis-Cancer Connection. Cells 2019, 8, 55. [Google Scholar] [CrossRef] [PubMed]
  75. Thomas, G. An encore for ribosome biogenesis in the control of cell proliferation. Nat. Cell Biol. 2000, 2, E71–E72. [Google Scholar] [CrossRef] [PubMed]
  76. Li, D.; Wang, J. Ribosome heterogeneity in stem cells and development. J. Cell Biol. 2020, 219, e202001108. [Google Scholar] [CrossRef]
  77. Zhang, Y.; Wolf, G.W.; Bhat, K.; Jin, A.; Allio, T.; Burkhart, W.A.; Xiong, Y. Ribosomal protein L11 negatively regulates oncoprotein MDM2 and mediates a p53-dependent ribosomal-stress checkpoint pathway. Mol. Cell Biol. 2003, 23, 8902–8912. [Google Scholar] [CrossRef]
  78. Shi, Z.; Fujii, K.; Kovary, K.M.; Genuth, N.R.; Röst, H.L.; Teruel, M.N.; Barna, M. Heterogeneous ribosomes preferentially translate distinct subpools of mRNAs genome-wide. Mol. Cell 2017, 67, 71–83.e77. [Google Scholar] [CrossRef] [PubMed]
  79. Jiao, L.; Liu, Y.; Yu, X.Y.; Pan, X.; Zhang, Y.; Tu, J.; Song, Y.H.; Li, Y. Ribosome biogenesis in disease: New players and therapeutic targets. Signal Transduct. Target. Ther. 2023, 8, 15. [Google Scholar] [CrossRef]
  80. Wei, J.; Kishton, R.J.; Angel, M.; Conn, C.S.; Dalla-Venezia, N.; Marcel, V.; Vincent, A.; Catez, F.; Ferré, S.; Ayadi, L.; et al. Ribosomal Proteins Regulate MHC Class I Peptide Generation for Immunosurveillance. Mol. Cell 2019, 73, 1162–1173.e1165. [Google Scholar] [CrossRef]
  81. Szczerba, B.M.; Castro-Giner, F.; Vetter, M.; Krol, I.; Gkountela, S.; Landin, J.; Scheidmann, M.C.; Donato, C.; Scherrer, R.; Singer, J.; et al. Neutrophils escort circulating tumour cells to enable cell cycle progression. Nature 2019, 566, 553–557. [Google Scholar] [CrossRef]
  82. Dang, C.V. MYC, metabolism, cell growth, and tumorigenesis. Cold Spring Harb. Perspect. Med. 2013, 3, a014217. [Google Scholar] [CrossRef]
  83. Ruggero, D. The role of Myc-induced protein synthesis in cancer. Cancer Res. 2009, 69, 8839–8843. [Google Scholar] [CrossRef]
  84. Gordan, J.D.; Thompson, C.B.; Simon, M.C. HIF and c-Myc: Sibling rivals for control of cancer cell metabolism and proliferation. Cancer Cell 2007, 12, 108–113. [Google Scholar] [CrossRef] [PubMed]
  85. Yi, M.; Ban, Y.; Tan, Y.; Xiong, W.; Li, G.; Xiang, B. 6-Phosphofructo-2-kinase/fructose-2,6-biphosphatase 3 and 4: A pair of valves for fine-tuning of glucose metabolism in human cancer. Mol. Metab. 2019, 20, 1–13. [Google Scholar] [CrossRef] [PubMed]
  86. Jaworska, M.; Szczudło, J.; Pietrzyk, A.; Shah, J.; Trojan, S.E.; Ostrowska, B.; Kocemba-Pilarczyk, K.A. The Warburg effect: A score for many instruments in the concert of cancer and cancer niche cells. Pharmacol. Rep. 2023, 75, 876–890. [Google Scholar] [CrossRef] [PubMed]
  87. Sittewelle, M.; Kappès, V.; Zhou, C.; Lécuyer, D.; Monsoro-Burq, A.H. PFKFB4 interacts with ICMT and activates RAS/AKT signaling-dependent cell migration in melanoma. Life Sci. Alliance 2022, 5, 12. [Google Scholar] [CrossRef] [PubMed]
  88. Kiesel, V.A.; Sheeley, M.P.; Coleman, M.F.; Cotul, E.K.; Donkin, S.S.; Hursting, S.D.; Wendt, M.K.; Teegarden, D. Pyruvate carboxylase and cancer progression. Cancer Metab. 2021, 9, 1–13. [Google Scholar] [CrossRef]
  89. Lin, Q.; He, Y.; Wang, X.; Zhang, Y.; Hu, M.; Guo, W.; He, Y.; Zhang, T.; Lai, L.; Sun, Z.; et al. Targeting Pyruvate Carboxylase by a Small Molecule Suppresses Breast Cancer Progression. Adv. Sci. 2020, 7, 1903483. [Google Scholar] [CrossRef]
  90. Rattigan, K.M.; Brabcova, Z.; Sarnello, D.; Zarou, M.M.; Roy, K.; Kwan, R.; de Beauchamp, L.; Dawson, A.; Ianniciello, A.; Khalaf, A.; et al. Pyruvate anaplerosis is a targetable vulnerability in persistent leukaemic stem cells. Nat. Commun. 2023, 14, 4634. [Google Scholar] [CrossRef]
  91. Zhang, L.; Ridgway, L.D.; Wetzel, M.D.; Ngo, J.; Yin, W.; Kumar, D.; Goodman, J.C.; Groves, M.D.; Marchetti, D. The identification and characterization of breast cancer CTCs competent for brain metastasis. Sci. Transl. Med. 2013, 5, ra148–ra180. [Google Scholar] [CrossRef]
  92. Alix-Panabières, C.; Pantel, K. Challenges in circulating tumour cell research. Nat. Rev. Cancer 2014, 14, 623–631. [Google Scholar] [CrossRef]
Figure 1. IVIS imaging and analysis after treating CTC-derived clone-injected MBM mice with the RPL inhibitor (omacetaxine), the CDK4/6 inhibitor (palbociclib), or the combination of drugs. Placebo group was used as a positive control. (A) tumor burden of each mouse was monitored via IVIS imaging weekly. Three representative mice in each group are shown (5 mice in each cohort). IVIS quantification of intra- (B) or extracranial (C) mouse regions shows decrease in tumor burden in drug-treated groups, compared to placebo. See “Materials and Methods” for experimental details, and Tables S1 and S2 for statistical analyses. * p = 0.05, ** p = 0.01.
Figure 1. IVIS imaging and analysis after treating CTC-derived clone-injected MBM mice with the RPL inhibitor (omacetaxine), the CDK4/6 inhibitor (palbociclib), or the combination of drugs. Placebo group was used as a positive control. (A) tumor burden of each mouse was monitored via IVIS imaging weekly. Three representative mice in each group are shown (5 mice in each cohort). IVIS quantification of intra- (B) or extracranial (C) mouse regions shows decrease in tumor burden in drug-treated groups, compared to placebo. See “Materials and Methods” for experimental details, and Tables S1 and S2 for statistical analyses. * p = 0.05, ** p = 0.01.
Cancers 15 05263 g001
Figure 2. CTC-driven micro- and macrometastasis quantification. (A) representative mouse brains from four animal groups are shown. (B) Bar graphs show the percentage of macrometastases from distant sites. Incidence of metastasis was found to be statistically significant between placebo- and drug-treated groups. (C) Melanoma cells in the brain sections from four animal groups using H&E staining (yellow arrows). (D) Quantification of tumor cells in the brain sections. See “Materials and Methods” for experimental details. ns = non-significant, * p = 0.05, ** p = 0.01.
Figure 2. CTC-driven micro- and macrometastasis quantification. (A) representative mouse brains from four animal groups are shown. (B) Bar graphs show the percentage of macrometastases from distant sites. Incidence of metastasis was found to be statistically significant between placebo- and drug-treated groups. (C) Melanoma cells in the brain sections from four animal groups using H&E staining (yellow arrows). (D) Quantification of tumor cells in the brain sections. See “Materials and Methods” for experimental details. ns = non-significant, * p = 0.05, ** p = 0.01.
Cancers 15 05263 g002
Figure 3. Enumeration and examination of CTCs captured by the Parsortix platform. Red box shows enumeration of large CTC clusters in each group. Visualization of CTCs and CTC clusters in mouse blood (A), melanoma patients’ blood (B), or healthy donor blood (C). See “Materials and Methods” for experimental details.
Figure 3. Enumeration and examination of CTCs captured by the Parsortix platform. Red box shows enumeration of large CTC clusters in each group. Visualization of CTCs and CTC clusters in mouse blood (A), melanoma patients’ blood (B), or healthy donor blood (C). See “Materials and Methods” for experimental details.
Cancers 15 05263 g003
Figure 4. RNA-Seq analyses of Parsortix-harvested CTC-enriched cell populations in four different mouse groups. Mouse blood from each experimental group was collected, pooled, and analyzed using Parsortix. Cellular fractions enriched with CTCs were harvested, followed by RNA isolation and transcriptomic profiling. Hierarchical clustering of gene expression profiling showed differential clustering of CTCs from placebo mice versus CTCs from omacetaxine-treated mice (A), CTCs from placebo mice versus CTCs from palbociclib-treated mice (B), or CTCs from placebo mice versus CTCs from mice receiving combinatory drug treatment (C). Corresponding heat maps are shown (left panels, red = upregulated, blue = downregulated). Scatter plots display global gene ex-pression of CTC populations with significant log2 fold change in the placebo group (green dots) versus the drug-treated groups (red dots; middle panels). Right panels show statistically significant pathways resulting from the hierarchical clustering using Reactome analysis. Pathways highlighted in red are the molecular pathways related to metabolic profiling. Highlighted in yellow are the pathways containing the CTC RPL/RPS gene signature of MBM [28]. Molecular pathways associated with proliferation are highlighted in green. See “Materials and Methods” for experimental details.
Figure 4. RNA-Seq analyses of Parsortix-harvested CTC-enriched cell populations in four different mouse groups. Mouse blood from each experimental group was collected, pooled, and analyzed using Parsortix. Cellular fractions enriched with CTCs were harvested, followed by RNA isolation and transcriptomic profiling. Hierarchical clustering of gene expression profiling showed differential clustering of CTCs from placebo mice versus CTCs from omacetaxine-treated mice (A), CTCs from placebo mice versus CTCs from palbociclib-treated mice (B), or CTCs from placebo mice versus CTCs from mice receiving combinatory drug treatment (C). Corresponding heat maps are shown (left panels, red = upregulated, blue = downregulated). Scatter plots display global gene ex-pression of CTC populations with significant log2 fold change in the placebo group (green dots) versus the drug-treated groups (red dots; middle panels). Right panels show statistically significant pathways resulting from the hierarchical clustering using Reactome analysis. Pathways highlighted in red are the molecular pathways related to metabolic profiling. Highlighted in yellow are the pathways containing the CTC RPL/RPS gene signature of MBM [28]. Molecular pathways associated with proliferation are highlighted in green. See “Materials and Methods” for experimental details.
Cancers 15 05263 g004
Figure 5. Four RPL/RPS genes (RPL23, RPL35A, RPL6, and RPS18) were inhibited by drug treatment. (A) representation of fold-change difference between drug-treated groups and placebo. (B) The 3D crystal structure of a ribosome generated by the PyMOL software shows specific localization of RPL23, RPL35A, RPL6, and RPS18 (black arrows). The left two panels show direct physical interaction between RPL6 (highlighted in purple) and RPL35A (highlighted in blue) via hydrogen bonding.
Figure 5. Four RPL/RPS genes (RPL23, RPL35A, RPL6, and RPS18) were inhibited by drug treatment. (A) representation of fold-change difference between drug-treated groups and placebo. (B) The 3D crystal structure of a ribosome generated by the PyMOL software shows specific localization of RPL23, RPL35A, RPL6, and RPS18 (black arrows). The left two panels show direct physical interaction between RPL6 (highlighted in purple) and RPL35A (highlighted in blue) via hydrogen bonding.
Cancers 15 05263 g005
Figure 6. Metabolic flux in patient-derived Lin−/CTC fractions. (A) Representative Mito Stress Test kinetics after lineage separation of two patients and a healthy donor sample. (B) Pooled metabolic parameters for the entire patient dataset. (C) OCR kinetics and metabolic parameters of Patient 1, showing basal respiration, basal ECAR, and ATP-coupled respiration to be elevated in the patient Lin−/CTC fraction. Violin plots all represent normalized data. All analysis was performed using two-way ANOVA with Kruskal–Wallis multiple comparisons test. OCR parameters at or below background threshold that showed appropriate well kinetics were treated as “0” values. Basal ECAR from Mito Stress Tests represent per-well averages of the three ECAR reads preceding the first drug injection. Analyzed patient samples all displayed robust flux kinetics and all experiments included control healthy donor cells on the same Seahorse plate. ns = non-significant, **** p < 0.0001.
Figure 6. Metabolic flux in patient-derived Lin−/CTC fractions. (A) Representative Mito Stress Test kinetics after lineage separation of two patients and a healthy donor sample. (B) Pooled metabolic parameters for the entire patient dataset. (C) OCR kinetics and metabolic parameters of Patient 1, showing basal respiration, basal ECAR, and ATP-coupled respiration to be elevated in the patient Lin−/CTC fraction. Violin plots all represent normalized data. All analysis was performed using two-way ANOVA with Kruskal–Wallis multiple comparisons test. OCR parameters at or below background threshold that showed appropriate well kinetics were treated as “0” values. Basal ECAR from Mito Stress Tests represent per-well averages of the three ECAR reads preceding the first drug injection. Analyzed patient samples all displayed robust flux kinetics and all experiments included control healthy donor cells on the same Seahorse plate. ns = non-significant, **** p < 0.0001.
Cancers 15 05263 g006
Figure 7. Cell numbers and carbohydrate metabolism in CTC-derived clone after inhibitor treatment. (A) Representative images showing cell density and morphology after treatment, as well as cell counts (top right panel) using a Biotek Cytation 1 with Gen 5 3.11 software after metabolic flux analysis. (B) Fold-change RT-qPCR analyses, and analysis using two-way ANOVA with Kruskal–Wallis multiple comparison test. (C) Glycolytic Rate Assay kinetics during 3-day inhibitor treatment in the presence or absence of 10 µM BPTES (top panel) or 10 µM UK5099 (bottom panel) for the last 24 h of inhibitor treatment. ns = non-significant, * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.
Figure 7. Cell numbers and carbohydrate metabolism in CTC-derived clone after inhibitor treatment. (A) Representative images showing cell density and morphology after treatment, as well as cell counts (top right panel) using a Biotek Cytation 1 with Gen 5 3.11 software after metabolic flux analysis. (B) Fold-change RT-qPCR analyses, and analysis using two-way ANOVA with Kruskal–Wallis multiple comparison test. (C) Glycolytic Rate Assay kinetics during 3-day inhibitor treatment in the presence or absence of 10 µM BPTES (top panel) or 10 µM UK5099 (bottom panel) for the last 24 h of inhibitor treatment. ns = non-significant, * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.
Cancers 15 05263 g007
Table 1. Parsortix CTC quantification from CDXs.
Table 1. Parsortix CTC quantification from CDXs.
CTC ClustersPlaceboOmacetaxinePalbociclibPal/Oma
CTCs/mL
Single-cell60153510
2-cell200100
3-cell10000
4-cell10000
5-cell or greater5000
Table 2. Inhibitory effect on RPL/RPS gene expression (cpm value).
Table 2. Inhibitory effect on RPL/RPS gene expression (cpm value).
PlaceboOmacetaxinePalbociclibBoth DrugsCTC Clone
RPL122.141.714.817.9511.4
RPL135.13.145.553.2177.59
RPL18A12.431.771.7316.78
RPL1981.1108.64154.07143.3343.91
RPL263.291.713.274.356.5
RPL3715.3515.9822.966.7273.96
RPL381.913.855.62.3927.73
RPL73.292.4211.7315.61
RPL7A8.725.271146.48
RPS188.963.851.755.8222.42
RPL63.05111.657.07
RPL35A4.621.715.573.7869.09
RPL2313.012.434.817.3735.69
RPS1260.1748.0633.6381.217.56
RPS15A111120.35
RPS24111.754.3524.82
RPS2611113.51
RPS28111112.91
RPS55.112.549.4215.37
RPS71.911.712.53.125.03
RPSA12.394.567.0920.0528.92
Average11.210.1712.5614.8732.51
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bowley, T.Y.; Merkley, S.D.; Lagutina, I.V.; Ortiz, M.C.; Lee, M.; Tawfik, B.; Marchetti, D. Targeting Translation and the Cell Cycle Inversely Affects CTC Metabolism but Not Metastasis. Cancers 2023, 15, 5263. https://doi.org/10.3390/cancers15215263

AMA Style

Bowley TY, Merkley SD, Lagutina IV, Ortiz MC, Lee M, Tawfik B, Marchetti D. Targeting Translation and the Cell Cycle Inversely Affects CTC Metabolism but Not Metastasis. Cancers. 2023; 15(21):5263. https://doi.org/10.3390/cancers15215263

Chicago/Turabian Style

Bowley, Tetiana Y., Seth D. Merkley, Irina V. Lagutina, Mireya C. Ortiz, Margaret Lee, Bernard Tawfik, and Dario Marchetti. 2023. "Targeting Translation and the Cell Cycle Inversely Affects CTC Metabolism but Not Metastasis" Cancers 15, no. 21: 5263. https://doi.org/10.3390/cancers15215263

APA Style

Bowley, T. Y., Merkley, S. D., Lagutina, I. V., Ortiz, M. C., Lee, M., Tawfik, B., & Marchetti, D. (2023). Targeting Translation and the Cell Cycle Inversely Affects CTC Metabolism but Not Metastasis. Cancers, 15(21), 5263. https://doi.org/10.3390/cancers15215263

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