The ability to adapt to environmental stress, including therapeutic insult, contributes to tumor evolution and drug resistance. In suboptimal conditions, the integrated stress response (ISR) promotes survival by dampening cytosolic translation. We show that ISR-dependent survival also relies on a concomitant up-regulation of mitochondrial protein synthesis, a vulnerability that can be exploited using mitoribosome-targeting antibiotics. Accordingly, such agents sensitized to MAPK inhibition, thus preventing the development of resistance in BRAFV600E melanoma models. Additionally, this treatment compromised the growth of melanomas that exhibited elevated ISR activity and resistance to both immunotherapy and targeted therapy. In keeping with this, pharmacological inactivation of ISR, or silencing of ATF4, rescued the antitumoral response to the tetracyclines. Moreover, a melanoma patient exposed to doxycycline experienced complete and long-lasting response of a treatment-resistant lesion. Our study indicates that the repurposing of mitoribosome-targeting antibiotics offers a rational salvage strategy for targeted therapy in BRAF mutant melanoma and a therapeutic option for NRAS-driven and immunotherapy-resistant tumors.
Introduction
Overcoming resistance to cancer therapy remains a major clinical challenge. Despite recent breakthroughs in targeted therapy and immune checkpoint blockade (ICB; Sosman et al., 2012; Larkin et al., 2015), the clinical response to these therapies is often incomplete and/or transient. These findings indicate that combining various therapeutic agents and modalities will be needed to overcome treatment resistance (Holohan et al., 2013). Given that the development of new drugs is very slow (Pushpakom et al., 2019) and is affected by high attrition rates (Waring et al., 2015), the repurposing of existing de-risked compounds offers a fast track for the clinical implementation of novel and high-order drug combinations. However, a deeper understanding of the mechanisms underlying drug resistance is essential before novel combinations can be rationally designed.
The most commonly accepted explanation for the inexorable development of therapy resistance invokes specific genetic alterations that are acquired by chance before or during treatment (Holohan et al., 2013). However, recent findings indicate that a subset of cancer cells are capable of surviving therapeutic insult by engaging specific adaptation programs that confer them with drug-tolerant phenotypes (Sharma et al., 2010; Roesch et al., 2013; Ravindran Menon et al., 2015; Su et al., 2017; Shen et al., 2019). These drug-tolerant persister (DTP) cells provide a pool, commonly referred to as minimal residual disease (MRD), from which stable resistance is established. These findings indicate that eradicating MRD before stable resistance is acquired may offer new promising therapeutic avenues (Boumahdi and de Sauvage, 2020; Rambow et al., 2019). However, characterization of the cellular composition of MRD using single-cell approaches has recently highlighted that implementing such a strategy will come with its own challenges. Indeed, coemergence (within the same MRD lesion) of four very distinct drug-tolerant states was observed following exposure of proto-oncogene B-Raf (BRAF) mutant melanoma to MAPK-targeted therapy (Rambow et al., 2018). These included the starved melanoma cell (SMC) state sharing transcriptomic features of nutrient-deprived cells (Rambow et al., 2018), a neural crest stem-like cell (NCSC) state, an invasive or mesenchymal-like state that was recently renamed undifferentiated state (Tsoi et al., 2018; Rambow et al., 2018), and a hyperdifferentiated state. The two de-differentiated states, NCSC and undifferentiated/mesenchymal, which harbor cancer stem cell features, are considered particularly important drivers of tumor recurrence (Boshuizen et al., 2018; Rambow et al., 2018). Unfortunately, there are currently no clinically compatible approaches known to efficiently cotarget these two distinct cell populations. Whether these (and the other additional DTP subpopulations), despite harboring distinct gene expression signatures, exhibit common and actionable vulnerabilities has therefore become a key question.
Although the metabolic profile of cancer cells varies across patients, tumor types, and subclones within a tumor, there is emerging evidence that mitochondrial bioenergetics, biosynthesis, and signaling are required for tumorigenesis. Accordingly, several recent studies have demonstrated that mitochondrial biology is potentially a promising new area for cancer therapy (Weinberg and Chandel, 2015; Jagust et al., 2019). Critically, as our understanding of the biology of MRD across multiple tumor types increases, it is becoming clear that DTP cells exhibit a strong dependence on mitochondrial biology, even higher than that of their drug-naive counterparts (Jagust et al., 2019; Sharon et al., 2019; Davis et al., 2020). For instance, relapse-initiating cells in B-progenitor acute lymphoblastic leukemia are characterized by elevated levels of mitochondrial metabolism (Dobson et al., 2020). Similarly, a slow-cycling population of melanoma cells that emerge in cultures exposed to vemurafenib or cisplatin exhibits elevated oxidative phosphorylation, and targeting mitochondrial respiration blocks their emergence and delays drug resistance (Roesch et al., 2013). It remains unclear, however, whether metabolic reprogramming is the only underlying cause of this increased mitochondrial dependence and whether (all) DTP cells that coemerge following BRAF and MAPK kinase (MEK) coinhibition are equally sensitive to mitochondria-targeting agents.
DTP cells from multiple cancer types often exhibit elevated activation of the integrated stress response (ISR; Deng and Haynes, 2017) following drug exposure (Almanza et al., 2019; McConkey, 2017). As the activating transcription factor 4 (ATF4) is a key mediator of the ISR, its role in drug tolerance is increasingly recognized (Rzymski et al., 2009; Jewer et al., 2020; Ojha et al., 2019). ISR is an adaptive translation program, triggered by several intracellular and extracellular stressors, that converges on the phosphorylation and activation of the eukaryotic translation initiation factor 2 (eIF2) complex. Activation of this pathway eventually results in reduced global translation and entry into a quiescent state essential to survive the environmental stress (Min and Spencer, 2019). However, when the stress becomes overwhelming, the same pathway triggers apoptosis (Verheyden et al., 2019). Emerging evidence indicates that in mammals, the ISR is essential to convey mitochondrial disfunctions to the nucleus through a process known as retrograde signaling (Quirós et al., 2017). Conversely, the effect of ISR activation on mitochondrial activities is less clear. The observation that the rate of mitochondrial translation adapts to the influx of nuclear-encoded mitochondrial proteins (Richter-Dennerlein et al., 2016) raised the possibility that DTP cells may also reduce their mitochondrial translational rate. In contrast to this prediction, however, we show herein that activation of ISR in DTP cells promotes selective translation of a subset of mRNAs encoding for mitochondrial proteins, and thereby generates an anterograde signaling from the cytosol-boosting mitochondrial translation. This mechanism may explain the exquisite sensitivity of DTP cells to mitochondria-targeting agents such as uncouplers, and importantly identifies mitochondrial translation as one critical sensitive node.
Interestingly, specific antibiotics can be repurposed to inhibit mitochondrial protein synthesis. Like the bacteria from which they originate, the mitochondrial translational machinery of eukaryotic cells, in particular mitoribosomes, are inhibited by tetracyclines (Kalghatgi et al., 2013; Moullan et al., 2015). The antitumor efficacy of these antibiotics has been demonstrated in vitro (Ahler et al., 2013) and in vivo in hematologic malignancies (Škrtić et al., 2011; D’Andrea et al., 2016; Zhang et al., 2017; Ravà et al., 2018), with several clinical trials in acute myeloid leukemia and double-hit lymphomas currently ongoing (Reed et al., 2016). Critically, we show herein that targeting mitochondrial protein synthesis with antibiotics of the tetracycline family prevented emergence of most (three out of four) drug-tolerant subpopulations and delayed and even prevented the development of resistance to MAPK inhibition in BRAF mutant preclinical patient-derived xenograft (PDX) melanoma models (Patton et al., 2021) as well as in one melanoma patient. Finally, we also show efficacy of this treatment on models derived from patients with limited therapeutic options such as uveal melanomas (UMs), BRAF WT cutaneous melanomas, or melanomas with intrinsic or acquired resistance to targeted therapies and ICB.
Together, these findings indicate that targeting mitochondrial translation with antibiotics of the tetracycline family should be exploited to rationally design anticancer therapeutic regimens. Importantly, given the widespread clinical use of such agents, these approaches could be rapidly implemented into the clinic.
Results
Activation of the ISR increases mitochondrial translation
Phenotype switching into an undifferentiated drug-tolerant state can be induced in vitro by activating the ISR, leading to an ATF4-dependent down-regulation of the microphthalmia-associated transcription factor (MITF; Falletta et al., 2017). Accordingly, exposure of drug-naive melanoma cells to salubrinal, an ISR agonist (Boyce et al., 2005), increased levels of ATF4 and caused a concomitant down-regulation of MITF (Fig. 1 A).
Considering that mitochondrial translation adapts to the influx of nuclear-encoded mitochondrial proteins, the observation that drug-tolerant cells down-regulate cytosolic protein synthesis suggests that these cells may also reduce the activity of their mitochondrial translation machinery. Surprisingly, however, puromycin incorporation assay followed by mitoplast isolation upon salubrinal treatment demonstrated that ISR activation caused instead a dramatic increase in mitochondrial translation, despite the expected overall decrease in cytosolic translation (Fig. 1, B and C). To further investigate the underlying mechanism, we identified translationally regulated mRNAs upon salubrinal treatment by performing polysome profiling analyses followed by RNA sequencing. We identified 382 transcripts whose association with ribosomes significantly (adjusted P value <0.05) changed in response to ISR activation and thus in response to phenotype switching and acquisition of therapy resistance. As expected, the vast majority of the transcripts (90%) were depleted from the ribosomal fractions upon ISR activation, while only 10% showed enrichment (Fig. S1 A). Among those, 2.3% were mitochondrial mRNAs and, accordingly, Ingenuity Pathway Analysis showed enrichment for mitochondria-related terms (Fig. S1 B). These findings were further validated in both BRAF mutant (Fig. 1 D) and neuroblastoma RAS (NRAS) mutant (Fig. S1 C) cell lines by performing polysome profiling followed by quantitative RT-PCR (RT-qPCR) upon salubrinal treatment. Since the influx of nuclear-encoded mitochondrial proteins from the cytosol directly regulates mitochondrial translation (Richter-Dennerlein et al., 2016), the latter observation provides a likely explanation for the observed increase in mitochondrial protein synthesis.
Strikingly, in silico analysis highlighted that 69.3% (122 of 176) of all ribosome-bound mRNAs in salubrinal-treated cells contain G-quadruplexes (rG4s; P = 2.4 × 1005 Benjamini–Hochberg corrected value; Fig. 1, E and F). In keeping with this, validated and/or predicted rG4s were identified in 88% (15 of 17) of the ribosome-associated mRNAs encoding mitochondrial proteins. rG4s, which are noncanonical four-stranded structures formed by guanine-rich sequences, have been implicated in all steps of RNA metabolism (Dumas et al., 2021). The above observations suggest that rG4s may contribute to the specific recruitment of selected transcripts to the translation machinery in cells undergoing ISR (Fig. 1, E and F).
Overall, these data identify mitochondrial translation as a putative Achilles’ heel of cells engaging the ISR pathway, such as drug-tolerant cells. Of note, ISR signaling is activated not only in response to drug treatment: western blot analysis detected activation of the ISR, as demonstrated by the detection of P-eIF2α, ATF4, and C/EBP homologous protein (CHOP), in non–drug-exposed human melanoma lesions, including lesions that were resistant to both targeted and immune therapies (Fig. 1 G). These data therefore indicate that targeting mitochondrial translation may offer a promising therapeutic approach in a substantial number of patients/cases.
Tigecycline overcomes acquired resistance to MAPK inhibitors in vivo and significantly increases overall survival (OS) and progression-free survival (PFS)
We reasoned that the increased recruitment of transcripts encoding mitochondrial proteins in cells harboring an activated ISR may sensitize them to antibiotics that target mitochondrial translation. To test this hypothesis, we exposed several human melanoma cell lines (Table 1) that harbor transcriptional profiles reminiscent of two critical drug-tolerant states, namely invasive/undifferentiated and NCSCs, to increasing concentrations of both tigecycline and doxycycline. Both antibiotics promoted a significant and dose-dependent decrease in cell growth in all cell lines tested (Fig. 2, A and B; and Fig. S2, A–C). Importantly, western blot analysis for ISR regulators confirmed that the ISR pathway is activated in these cell lines before treatment. Moreover, exposure to the antibiotics alone or in combination with MAPK inhibitors (MAPKi) further exacerbated proteotoxic stress, as illustrated by an increase in CHOP activation (Fig. 2 C). Conversely, treatment with antibiotics of MM034, which do not display expression of ATF4 (Fig. 2 C), did induce a certain degree of growth inhibition; however, the decrease did not reach statistical significance (Fig. 2 A). These observations indicated that tetracyclines may offer a therapeutic option for ablation of the residual tumor cells emerging following exposure to MAPKi. We therefore tested whether addition of tigecycline to dabrafenib-trametinib (DT), a standard-of-care treatment for patients with BRAF mutant melanoma, delayed or prevented the onset of resistance in two different BRAFV600E melanoma PDXs (Mel-006 and Mel-015). As expected, all mice from both cohorts initially responded to DT but eventually developed resistance (Fig. 3 and Fig. S2, D and E). In the Mel-015 cohort, the addition of antibiotics, whether added from the start of the treatment (dabrafenib-trametinib-tigecycline [DTT]) or after lesions reached MRD (DT + T [tigecycline]), significant delayed the development of resistance (Fig. 3 A) and increased PFS and OS (Figs. 3 B and S2 D). Strikingly, in the Mel-006 cohort, DTT treatment resulted in complete remission in 11 of 14 mice (Fig. 3, C and D; and Fig. S2 E). Importantly, the curative nature of this treatment was confirmed by interrupting the DTT treatment in another cohort of mice after tumors had regressed (Fig. 3, C and D; and Fig. S2 E). Notably, dabrafenib plus tigecycline (D + Tige) was as effective as DT treatment at inducing tumor regression in this particular PDX model. This treatment even significantly increased PFS and OS compared with DT (Fig. 3, C and D; and Fig. S3 B). Thus, addition of tigecycline successfully delays and/or prevents acquisition of resistance to targeted therapy in BRAF mutant skin melanoma preclinical models.
Tigecycline suppresses known drug-resistant cell populations
We have now provided evidence that tetracyclines affect the growth of both mesenchymal and NCSC de-differentiated drug-tolerant cells in vitro. To assess the sensitivity of these and other drug-tolerant subpopulations in vivo, we treated Mel-006 with either D + Tige or DT until reaching MRD. Compared with DT treatment, the D + Tige combination successfully eradicated most of the NCSCs (NGFR+AQP1+), undifferentiated/mesenchymal-like cells (AXL+), and SMCs (CD36+) drug-tolerant cells (Fig. 4, A–C). Similarly, in two different BRAF mutant PDX cohorts (Mel-006 and Mel-015) treated with the triple DTT combination, antibiotic treatment led to the eradication of NCSCs (NGFR+AQP1+), undifferentiated/mesenchymal-like cells (AXL+) and pseudo-starved cell population SMCs (CD36+). The hyperpigmented cell population (MITF+MLANA+) was the only drug-tolerant population to resist tigecycline combinatorial treatments (Fig. 4, B–E). Note that, consistent with tigecycline targeting mitochondrial translation and not metabolism, spatial metabolomics indicated that the antibiotic-containing treatment (DTT) did not induce any significant metabolic shift compared with DT alone (Fig. S3, A–N). These data indicate that a combination of antibiotics and MAPKis affects known drug-resistant populations by targeting mitochondrial translation in both de-differentiated (NCSCs) and undifferentiated/invasive states in the absence of significant metabolic rewiring.
Tetracyclines restrain the growth of therapy-resistant melanoma
Data in Fig. 2 C indicate that cells intrinsically resistant to MAPK inhibition (see MM011) and/or immunotherapy (see IGR37) may also be sensitive to mitochondrial translation inhibition. Therefore, we exposed BRAF WT melanoma cell lines, including a UM line and the immunotherapy-resistant IGR37 and YUMM 1.7 cell lines (together with the sensitive counterpart YUMMER 1.7; see Table 1) to increasing concentrations of tigecycline or doxycycline (Fig. 5, A and B; and Fig. S4, A and B). Both antibiotics promoted a significant and dose-dependent decrease in cell growth in all cell lines. We next tested the efficacy of these antibiotics on the in vivo growth of several therapy-resistant PDX models (Fig. 5, C–E). Daily treatment with tigecycline significantly delayed the growth of NRAS mutant PDXs (Mel-020 and Mel-083), including one model derived from a patient that had progressed on immunotherapy (Mel-083; Fig. 5, C and D; and Fig. S4 E). Milder, yet significant, results were also obtained in the drug-naive BRAFV600E mutant cutaneous melanoma models (Mel-006 and Mel-015), which exhibit sensitivity to MAPKi (Fig. S4, C and D). Likewise, C57BL/6 mice engrafted with the immunotherapy-resistant YUMM 1.7 cells (Fig. S4 F) displayed sensitivity to daily monotreatment of tigecycline (Fig. S4 G). Furthermore, tigecycline treatment of the UM PDX model Mel-077 and NRASQ61R mutant cutaneous melanoma Mel-083, both derived from patients that had progressed on immune checkpoint inhibitors (ICI, pembrolizumab + temozolomide for Mel-077 and ipilimumab + nivolumab for Mel-083), was also sufficient to significantly delay tumor growth (Fig. 5, D and E; and Fig. S4 H). Lastly, to investigate whether the antibiotic treatment negatively impact the response to ICIs, the Mel-006 model was engrafted in hu-NOG-EXL mice humanized with CD34+ cells (Fig. S5 G). Mice were treated with anti–PD-1 alone or in combination with tigecycline (50 mg/kg) once tumors reached 100 mm3 (Fig. 5 F). Treatment with tigecycline did not overtly ameliorate, nor impair, the response to anti–PD-1 (Fig. 5 F). Similar results were obtained with YUMMER allografts (data not shown). Together, these data indicate that although tetracycline does not sensitize to ICI, monotreatment with this antibiotic is sufficient to significantly delay progression of melanoma lesions that are insensitive to targeted and immune therapies in vivo.
Doxycycline sensitizes refractory lesions to MAPKi in preclinical settings and in one patient
We also assessed whether tetracyclines can overcome intrinsic resistance to targeted therapy. To this end, we derived a cutaneous melanoma PDX model (Mel-007) from a lesion that exhibited intrinsic resistance to a combination of BRAF-MEK inhibitors. Exposure of a cohort of these mice to DT versus DTT showed that the addition of tigecycline promoted a more robust antitumor effect and even promoted tumor regression, in contrast to DT treatment (Fig. 6 A).
Doxycycline is a broad-spectrum antibiotic that also belongs to the tetracycline family. It exerts its antibacterial action by binding to the 30S ribosomal subunit and thereby blocking ribosome biogenesis. It can be administered orally for extended periods of time with only minor adverse effects in patients (Tan et al., 2011).
To assess the relevance of our in vitro and preclinical findings, we followed the clinical course of one 73-yr-old female patient at our clinic. The patient was diagnosed with stage III melanoma in 2011 and with a relapse of her disease in 2016 after an episode of jaundice, which was related to a tumor mass in her gallbladder and bile duct. Histology and immunohistochemistry of the bile duct biopsy were performed (Fig. 6 B). H&E (left) shows a highly anaplastic tumor with large epithelioid cells and pleomorphic nuclei; the presence of brown pigment associated with these cells was highly suggestive for a melanoma metastasis. The immunostaining for S100 (right), one of the most sensitive markers to identify melanoma, was diffusely positive at nuclear and cytoplasmic levels, therefore confirming the diagnosis of melanoma metastasis. Baseline positron-emission tomography/computed tomography (PET-CT) showed metastatic tumor at the aforementioned locations as well as in the liver and in the gallbladder (Fig. 6 C). Subsequently BRAF-MEK inhibitor therapy was initiated. Response assessment after four cycles showed that while the liver metastasis had responded to the treatment, the gallbladder did not. Due to a skin infection (an acne-like rash) the patient received doxycycline during her fifth BRAF-MEK inhibitor cycle for a total of 12 d. Strikingly, the gallbladder metastasis started to regress after six cycles, as measured by CT scan, and the PET-CT scan after eight cycles showed a complete response. The patient was under BRAF-MEK inhibitor therapy and had a persistent response for >36 mo (Fig. 6 C). This particular case study highlights the potential clinical benefit of using tetracyclines to sensitize intrinsically resistant lesions to MAPK inhibitors.
Activation of the ISR predicts durable responses to an antibiotic treatment
We postulated that the differential efficacy of the antibiotic combinatorial treatment in preventing the development of resistance between Mel-006 and Mel-015 BRAFV600E mutant melanoma PDXs may result from a difference in the kinetic and/or extent at which the ISR pathway is activated before applying the therapy. Consistently, ATF4 activation was readily detectable by western blotting, before the start of the therapy (T0), in the Mel-006 PDX model, an observation that may explain the high sensitivity to the antibiotic treatment and the long-lasting response (Figs. 1 G and 7 A). In contrast, ATF4 activation could only be detected in lesions from Mel-015 following therapy exposure (Figs. 1 G and 7 B). Importantly, basal ATF4 levels were also detectable in drug-naive lesions from the tigecycline-sensitive models Mel-077 UM model and the Mel-020 and Mel-083 NRASQ61R mutant models (Fig. S5, A–C). This observation may explain why the OS and PFS were comparable in cohorts exposed to the antibiotics at MRD or from the beginning of the treatment. These data suggest that the extent of the response to mitoribosome-targeting antibiotics is dependent on the level of chronic ISR activation. In keeping with this, cotreatment with the ISR inhibitor ISRIB, which blocks signaling downstream of the eIF2α kinase and thereby dampens ATF4 activation and ISR engagement (Fig. S5 D), reduced CHOP induction (Fig. S5 D) and abolished the apoptotic response observed in a series of melanoma cells lines exposed to doxycycline (Fig. 7 C). Consistently, silencing of ATF4 was also sufficient to overcome doxycycline-induced growth inhibition (Fig. S5, E and F). Altogether, these data indicate that dependency on mitochondrial translation is a consequence of ATF4-mediated ISR activation. Importantly, these findings also identify ATF4 as a predictive biomarker for the antitumor efficacy of the antibiotic treatment.
Discussion
Cancer cells have the ability to survive a wide range of insults—including therapy exposure—through activation of an ATF4-ISR–dependent translation remodeling. We show herein that this pathway, in addition to promoting an overall decrease in cytosolic translation, enhances mitochondrial protein synthesis and selective translation of mRNAs encoding mitochondrial proteins in the cytosol and thereby establishes an unexpected vulnerability to inhibitors of mitoribosome assembly or function. Remarkably, this vulnerability can be targeted using the FDA-approved and widely used mitoribosome-targeting antibiotics. Accordingly, we show that tigecycline and doxycycline can decrease the growth and viability of a wide range of melanoma cells, exhibiting elevated ISR activity, including de-differentiated DTP cells or cells that have acquired resistance to targeted and immune therapies. We show that ATF4 inhibition or pharmacological inactivation of ISR rescues the antitumor effects of such antibiotics, thus confirming that this dependence relies on elevated ISR activity.
It has been proposed that internal ribosome entry site–dependent translation is one mechanism allowing translation of transcripts encoding survival proteins in cells exposed to stress and engaging an ISR. In fact, the translation of ATF4 itself depends on the internal ribosome entry site (Chan et al., 2013). Here we show that transcripts that are recruited to the ribosomes are enriched in rG4 sequences, raising the possibility that these elements also contribute to the selective translation of a series of transcripts, including those encoding mitochondrial proteins, in stressed cells.
An emerging theme in cancer biology and drug resistance is the reliance of different lesions and DTP cells, across multiple cancer types, on mitochondrial integrity and functions (Chen, 2012; Haq et al., 2013; Jagust et al., 2019; Faubert et al., 2020; Davis et al., 2020). There is even evidence that primary resistance to immunotherapy may also be driven by the up-regulation of mitochondrial translation (Jerby-Arnon et al., 2018; Poźniak et al., 2019). For instance, it has been shown that DTP cells exhibit elevated levels of oxidative phosphorylation, and that targeting mitochondrial respiration can delay resistance to targeted therapy (Roesch et al., 2013). Unfortunately, the use of mitochondrial respiration uncouplers to treat obesity and other metabolic diseases led to serious systemic side effects in patients (Childress et al., 2018). These observations have dampened the excitement toward these molecules and delayed their exploitation in the clinic. Importantly, we show herein that metabolic reprogramming is not the only underlying cause of the increased mitochondrial dependency of DTP cells. Our data therefore provide an alternative approach for the targeting of mitochondrial biology, through the repurposing of mitoribosome-targeting antibiotics. We argue that this approach may even be less susceptible to the development of resistance through metabolic adaptation. Accordingly, we show that tetracyclines addition to the standard-of-care targeted therapy (i.e., DT) extended PFS and often prevented relapse in a BRAF mutant preclinical PDX model (Mel-006). Mechanistically, we showed that antibiotic treatment compromises the emergence and/or survival of multiple drug–tolerant subpopulations, including NCSCs, mesenchymal/undifferentiated, and pseudo-starved subpopulations at MRD. This is an exciting observation given that these de-differentiated melanoma cells were shown to contribute to the development of resistance to various anti-melanoma therapies, including T cell transfer therapy (Landsberg et al., 2012; Lee et al., 2020; Mehta et al., 2018; Rambow et al., 2018; Boshuizen et al., 2020).
Importantly, although providing a significant increase in survival, the triple DTT combination did not prevent relapse in another BRAF mutant preclinical PDX model (Mel-015). The relatively lower benefit from the tigecycline treatment in this model correlated with undetectable levels of ISR activation (as measured by lack of ATF4 expression) before treatment. ATF4 levels also predicted responses of therapy-resistant melanomas, such as UM and BRAF WT cutaneous melanomas, to the tetracycline treatment. These observations offer a new avenue for the treatment of these diseases with no, or only limited, therapeutic options and indicate that ATF4 levels are predictive of responses to tetracycline exposure, either as a single agent or in combination with standard of care.
Recent reports highlighted that brain metastasis across different cancers (Chen et al., 2007; Chen, 2012; Cheng et al., 2019), including melanoma (Fischer et al., 2019; Sundstrøm et al., 2019), display increased oxygen consumption and are therefore more dependent on OXPHOS. Considering that tetracycline family antibiotics can cross the blood–brain barrier (Nau et al., 2010), this approach may also be a viable option for patients with this life-threatening and frequently untreatable condition.
Note that, in addition to their ability to inhibit mitochondrial translation, it remains possible that the antitumor effect of tetracyclines is multifactorial. Recent findings highlighted that mitochondrial translation deficiencies impair lysosome acidification and, thereby, autophagy (Yagi et al., 2021). Given the importance of autophagy in the biology of drug tolerance (Kinsey et al., 2019), this observation may also partly explain the efficacy of tetracyclines in delaying/preventing the development of resistance to MAPKi. Moreover, both dabrafenib and trametinib are partially metabolized by CYP3A4 (Lawrence et al., 2014), which is inhibited by tetracycline. Its coadministration may therefore raise the active concentrations of dabrafenib and trametinib and thereby enhance the antitumor efficacy of DT (Bassi et al., 2004). We also cannot rule out that, in addition to cancer cell intrinsic effects, the dramatic responses observed in the preclinical and clinical settings are a consequence of tetracycline-mediated remodeling of the tumor microenvironment, such as the (intratumor) microbiome. This may be important, as the efficacy of anti–PD-1 therapy has been shown to be affected by the microbiome (Matson et al., 2018; McQuade et al., 2019). However, most common symbiotic bacteria are resistant to doxycycline, which is a broad-spectrum antibiotic that can be administered orally for long periods of time with no or minor toxic effects (Tan et al., 2011). Its administration is therefore unlikely to affect significantly the patients’ microbiome.
In conclusion, our data indicate that the repurposing of mitoribosome-targeting antibiotics offers a new promising therapeutic avenue for the treatment of a large spectrum of melanoma patients, including those with limited therapeutic options. Moreover, we provide evidence that patient stratification should be guided by ATF4 levels, which could be used as an accompanying biomarker to predict efficacy. Importantly, this drug-repositioning strategy can be easily implemented in the clinic.
Materials and methods
Study design
The objectives of this study were (a) to identify novel metastatic melanoma vulnerabilities, (b) to prevent or delay acquisition of therapy resistance, (c) to sensitize therapy-resistant lesions to therapy, and (d) to exploit this knowledge for therapeutic benefit in a preclinical setting. To test this objective, we made use of several clinically relevant mouse models of metastatic cutaneous and UM with different mutational backgrounds and therapy sensitivities (i.e., with intrinsic or acquired therapy resistance profiles). Written informed consent was obtained from patients, and all procedures involving human samples were approved by the Universitair Ziekenhuis (UZ) Leuven/Katholieke Universiteit (KU) Leuven Medical Ethical Committee (S63799) and performed in accordance with the principles of the Declaration of Helsinki and with General Data Protection Regulations. The mouse experiments were approved by the KU Leuven animal ethical committee under Ethical Committee for Animal Testing (Ethische Commissie Dierproven, ECD) P038-2015 and P164-2019 and performed in accordance with the internal, national, and European guidelines of animal care and use. Mice were evaluated for OS, PFS, and tumor growth daily. Cell lines were monitored for cell growth and cell death. The melanoma patient was monitored for tumor load and tumor response to the treatments. Each in vitro experiment was repeated at least three times (as specified in the figure legend) to ensure reproducibility. For in vivo experiments, sample size was calculated using the software G power v3.1.9.2 with a power of 80–95% and an α-error probability of 0.05%, assuming an effect size of 0.25 (based on previous experiments and pilot studies). This number was not altered during the experiment. This study included several cell cultures, two animal models, and one melanoma patient. Investigators were not blinded.
Cell lines
MM383 and WM852 cell lines (a kind gift from Göran Jönsson), SK-MEL-28 cells, and the IGR37 cell line (a kind gift from Irwin Davidson) were grown in RPMI 1640 (Gibco BRL Invitrogen) supplemented with 10% FBS (Gibco BRL Invitrogen). The patient-derived low-passage MM cell lines (a gift from G.-E. Ghanem) were grown in F-10 (Gibco BRL Invitrogen), supplemented with 10% FBS (Gibco BRL Invitrogen). The YUMM 1.7 and YUMMER 1.7 cell lines (a kind gift from Marcus Bosenberg; Meeth et al., 2016; Wang et al., 2017) were grown in DMEM/F12 (50:50; Gibco BRL Invitrogen) medium supplemented with 10% FBS (Gibco BRL Invitrogen). The UM92.1 cell line (a kind gift from Aart Jochemsen) were grown in a 1:1 mixture of RPMI 1640 (Gibco BRL Invitrogen) and DMEM/F12 (Gibco BRL Invitrogen) medium supplemented with 10% FBS (Gibco BRL Invitrogen).
All cell lines used are of human origin, save for the YUMM 1.7 and YUMMER 1.7 cell lines, derived from male C57BL/6J (Mus musculus). Gender of the patients from whom the cell cultures were derived is as follows (female: F; male: M): SK-MEL-28, M; MM034 and MM011, F; MM099, M; MM165, M; IGR37, M; YUMM 1.7, M; UM92.1, F. All cell lines were confirmed negative for mycoplasma before the study using the MycoAlert Mycoplasma Detection Kit (Lonza) according to the manufacturer’s specifications.
Cell viability assays
For colony formation assays, cells were plated in 6-well plates at the appropriate density (1.5 × 104/well in the case of the MM011, MM034, MM099, MM165, MM383, WM852, IGR37, YUMM 1.7, and YUMMER 1.7 cells and 8 × 104/well in the case of the UM92.1 cells). Cells were then treated with increasing amounts of tigecycline and cultured for 3 d (YUMM 1.7, YUMMER 1.7, and UM 92.1), 5 d (WM852 and IGR37), or 7 d (MM011, MM034, MM165, and MM383). Cells were then washed twice with PBS, fixed, and stained for 15 min with 1% crystal violet in 35% methanol solution. Surface occupancy was measured using ImageJ.
For siRNA-mediated knockdown, cells were plated in duplicate in 6-well plates at the confluency described above and transfected with Lipofectamine 2000 (Thermo Fisher Scientific) according to manufacturer instructions with 25 nM of nontargeting siRNA (control, SMART-pool; Dharmacon) or with 25 nM siATF4 (SMART-pool; Dharmacon). 72 h after transfection, cells were collected (for protein extraction) or stained with crystal violet and quantified as described above.
For IncuCyte proliferation assays, cells were plated in 96-well plates (TPP) at the appropriate density (MM011, MM034, MM165, MM099, IGR37, and MM383, 1.5 × 104/well; WM852, 5 × 103/well; YUMM 1.7, YUMMER 1.7, and UM 92.1, 2.5 × 103/well). Cells were treated with increasing amounts of tigecycline or doxycycline and cultured for 72 h. Apoptotic cells were labeled with IncuCyte Caspase 3/7 Green Apoptosis Assay Reagent (Essen BioScience). Four images per well were taken at 2-h intervals using an IncuCyte ZOOM system (Essen BioScience). The percentage of cell confluency and fluorescent green counts indicating apoptotic cells were measured and analyzed with IncuCyte ZOOM software.
Polysome profiling
SK-MEL-28 cells (5 × 15-cm dishes per each condition) were plated to have 70% confluency after 72 h. The following day, cells were treated with 20 µM salubrinal (Sigma-Aldrich). 72 h after the start of the treatment, cells were treated with 100 µg/ml of cycloheximide (Sigma-Aldrich) for 12 min at 37°C, collected, and resuspended in lysis buffer (30 mM Tris-HCl, 150 mM KCl, 10 mM MgCl2 supplemented with 1 mM DTT [Sigma-Aldrich], 100 µg/ml cycloheximide, 20 U/µl SUPERase-IN RNase Inhibitor [Invitrogen, Thermo Fisher Scientific], and 1× Halt Protease and Phosphatase Inhibitor Single-Use Cocktail [Life Technologies] before the start of the experiment). Lysates were then incubated agitating at 4°C for 35 min, and then centrifuged at 17,000 rcf for 15 min at 4°C. Lysates were loaded on a sucrose gradient (the linear sucrose gradient 5–20% was generated from two different solutions, sucrose 5 and 20%, made with buffer G (20 mM Tris-HCl, 100 mM KCl, 10 mM MgCl2 supplemented with 1 mM DTT and 100 µg/ml cycloheximide before the start of the experiment). Samples were then centrifugated in an SW41Ti rotor (Beckman Coulter) at 37,000 rpm for 170 min at 4°C. The fractions were obtained with a Biological LP System (Bio-Rad). 14 fractions were collected from each sample, with each fraction having a final volume of 600 µl. From the initial 14 fractions, 4 final samples were obtained (by pulling together some of them): 40S, 60S, 80S, and polysomes.
RNA extraction
For RNA extraction from the fractions, each fraction was digested at 37°C for 90 min in a digestion mix of proteinase K (final concentration 100 µg/ml; Sigma-Aldrich) and 1% SDS. Phenol acid chloroform (5:1; Sigma-Aldrich) and 10 mM NaCl were then added. Samples were centrifuged at 16,000 rcf for 5 min at 4°C. The upper aqueous phase was transferred in a new tube, and 1 ml of isopropanol was added. Samples were stored at −80°C overnight to precipitate the RNA. The following day, the samples were centrifuged at 16,000 rcf for 40 min at 4°C. The supernatant was discarded, and the pellet was washed with 500 µl of 70% EtOH, then centrifuged again at 16,000 rcf for 5 min at 4°C. Pellet were air dried and resuspended in diethylpyrocarbonate water. No DNase treatment was performed (not needed). Input RNAs were extracted with TRIzol (Invitrogen) according to the manufacturer’s instructions.
RT-qPCR
RNA was reverse transcribed using the High-Capacity complementary DNA Reverse Transcription Kit (Thermo Fisher Scientific) on a Veriti 96-well thermal cycler (Thermo Fisher Scientific). Gene expression was measured by qPCR on a QuantStudio 5 (Thermo Fisher Scientific) and normalized in qbase + 3.0 (Biogazelle) using 28S and 18S as reference genes (for polysome profiling experiments) or the average of HPRT, TBP, and UBC. Sequences of the primers are indicated in Table S1.
RNA sequencing and data analysis
Samples were prepared for sequencing with TruSeq Stranded Total RNA kit (Illumina) according to the manufacturer’s instructions. Libraries were sequenced on an Illumina NextSeq according to the manufacturer’s instructions, generating PE high output configuration cycles: (R1: 100) − (I1: 6) − (I2: 0) − (R2: 100). Quality control of raw reads was performed with FastQC v0.11.5. Adapters were filtered with eautils v1.2.2.18. Splice-aware alignment was performed with STAR against the human hg19 human genome assembly. The number of allowed mismatches was 2. Reads that mapped to more than one site of the reference genome were discarded. The minimal score of alignment quality to be included in count analysis was 10. Resulting Sequence Alignment/Map (SAM) and binary SAM alignment files were handled with Samtools v0.1.19.24. Quantification of reads per gene was performed with HT-Seq count v0.5.3p3. Count-based differential expression analysis was done with R-based (R Foundation for Statistical Computing) Bioconductor package DESeq. Reported P values were adjusted for multiple testing with the Benjamini–Hochberg procedure, which controls false discovery rate. Differentially expressed genes were selected at false discovery rate 0.05. For the comparison of transcriptional and ribosome association changes in Fig. S1 A, the R package ggplot2 was used. The data generated in this study are accessible through GEO with accession number GSE177075.
rG4 element analysis
rG4 elements were predicted within the mRNA of the selected gene lists by using the QGRS Mapper tool (https://bioinformatics.ramapo.edu/QGRS/index.php; Kikin et al., 2006) with default parameters, and transcript sequences from the GENCODE V32 release, hg38 human genome assembly (https://www.gencodegenes.org). Being the rG4 score distribution bimodal with two well-separated peaks, the high-confidence score threshold was selected to be 21, i.e., the value corresponding to the point at which the first peak fell to 0 predicted rG4s before the start of the second, higher-scoring peak. Experimentally determined in vitro rG4s were obtained from the Kwok et al. (2016) and Guo and Bartel (2016) datasets. Enrichment of rG4s in the selected gene lists was computed by Fisher’s exact test, using the whole transcriptome as a background. Finally, the enrichment of regulatory factor binding sites (RNA-binding proteins and microRNAs) and cis-elements sites in the 5′ and 3′ UTRs of rG4-containing mRNAs was computed by the Regulatory Enrichment function of the AURA2 database (Dassi et al., 2014), using a 0.05 threshold on the BH-adjusted enrichment P value.
Puromycin incorporation assay (SUNsET)
SUNsET was performed as described in Schmidt et al. (2009). Briefly, ≈80% confluent adherent cells were washed twice in 1× PBS and subsequently pulsed with puromycin-containing medium (InvivoGen, 10 µg ml−1) for 10 min. The cells were then supplemented with normal medium for 60 min before downstream applications (chase). Due to the fact that puromycin is a structural analogue of aminoacyl transfer RNAs, it gets incorporated into the nascent polypeptide chain and prevents elongation. Therefore, when used for reduced amounts of time, puromycin incorporation in neosynthesized proteins directly reflects the rate of mRNA translation in vitro. Puromycin incorporation was measured by western blotting using an antibody that recognizes puromycin.
Cellular fractionation and mitoplast isolation
Briefly, mitochondria were purified from 4–6 × 107 cells using a mitochondria isolation kit for cultured cells (Thermo Fisher Scientific) according to manufacturer instructions; all buffers were supplemented with 60 U ml−1 SUPERase-In (Ambion) and 1× Halt Protease and Phosphatase Inhibitor Single-Use Cocktail (Life Technologies). Mitoplasts were obtained by incubating purified mitochondria in RNase A–containing hypotonic buffer (Hepes, pH 7.2, supplemented with 1× Halt Protease and Phosphatase Inhibitor Single-Use Cocktail and 10 µg/ml RNase A [Roche]) for 20 min on ice and subsequently incubated for 10 additional min at room temperature in order to remove all possible cytosolic RNA contaminants. The purified mitoplasts were then washed three times with mitoplast isolation buffer (250 nM mannitol, 5 mM Hepes, pH 7.2, 0.5 mM EGTA, 1 mg/ml BSA supplemented with 60 U ml−1 Superase-In [Ambion] and 1× Halt Protease and Phosphatase Inhibitor Single-Use Cocktail).
Ascorbic acid peroxidase (APEX) assay
SK-MEL-28 cells were seeded at 70% confluency in 15-cm dishes (one per condition) and transiently transfected with 10 µg mito-V5-APEX2 plasmid (#72480; Addgene; Lam et al., 2015) using Lipofectamine 2000 (Thermo Fisher Scientific). 24 h after transfection, the medium was removed, and the cells were subsequently treated with either 1× PBS (control) or 30 µM tigecycline. 24 h after treatment, which corresponds to 48 h after transfection, cells were labeled by a 30-min incubation with 500 µM biotin-phenol at 5% CO2, 37°C (Iris Biotech), followed by a 1-min incubation with 1 mM H2O2 at room temperature. Medium was quickly removed, and the biotin-phenol reaction was quenched by three washes with quencher solution (10 mM sodium ascorbate, 5 mM Trolox, and 10 mM sodium azide diluted in PBS). The cells were then also washed in 1× PBS, scraped in 4 ml cold 1× PBS, and collected for lysis. Each sample was lysed in 50 µl polysome lysis buffer supplemented with 1× Halt Protease and Phosphatase Inhibitor Single-Use Cocktail.
Western blotting
Western blotting experiments were performed using the following primary antibodies: vinculin (V9131, clone hVIN-1; Sigma-Aldrich, 1:5,000), histone 3 (#4499, clone D1H2; Cell Signaling Technology, 1:1,000), MITF (ab12039; Abcam, 1:1,000), ATF4 (#11815, clone D4B8; Cell Signaling Technology, 1:1,000), ATF5 (SAB4500895; Sigma-Aldrich, 1:500), eIF2α-Tot (#5324, clone D7D3; Cell Signaling Technology, 1:1,000), eIF2α-Phospho-S51 (#3398, clone D9G8; Cell Signaling Technology, 1:1,000), CHOP (#2895, clone L63F7, 1:1000), puromycin (MABE343, clone 12D10; Merck-Millipore, 1:10,000), streptavidin (#N1000; ThermoFisher Scientific, 1:5,000), β-actin (#4970, clone 13E5; Cell Signaling Technology, 1:1,000), and Tom20 (sc-17764, clone F-10; Santa Cruz Biotechnology, 1:1,000). The following HRP-linked secondary antibodies were used: anti-mouse IgG (NA931-1ML; Sigma-Aldrich, 1:10,000) and anti-rabbit IgG (NA934-1ML; Sigma-Aldrich, 1:10,000). Relative protein levels were measured using ImageJ.
PDX experiments
The cutaneous melanoma PDX models are part of the Trace collection (https://www.uzleuven-kuleuven.be/lki/trace/trace-leuven-pdx-platform) and were established using metastatic melanoma lesions derived from patients undergoing surgery as part of standard treatment at UZ Leuven. The UM model derives from a lymph node metastasis and is a kind gift of M. Herlyn (The Wistar Institute, Philadelphia, PA). Written informed consent was obtained from patients and all procedures involving human samples were approved by the UZ Leuven/KU Leuven Medical Ethical Committee (S63799) and performed in accordance with the principles of the Declaration of Helsinki and with GDPR regulations. The experiments were approved by the KU Leuven animal ethical committee under ECDs P038-2015 and P164-2019 and performed in accordance with the internal, national, and European guidelines of animal care and use. Tumor pieces were implanted subcutaneously in the interscapular fat pad of female NMRI nude BomTac:NMRI-Foxn1nu, 4-wk-old females (Taconic Biosciences). Mice were maintained in a semi-specific pathogen–free facility under standard housing conditions with continuous access to food and water. The health and welfare of the animals was supervised by a designated veterinarian. The KU Leuven animal facilities comply with all appropriate standards (cages, space per animal, temperature [22°C], light, humidity, food, and water), and all cages are enriched with materials that allow the animals to exert their natural behavior. Mice used in the study were maintained on a diurnal 12-h light/dark cycle. PDX models Mel-006, Mel-015, and Mel-020 were derived from a female, male, and female drug-naive patients, respectively. The UM Mel-077 sample was derived from a male patient progressing on pembrolizumab and temozolomide. Mel-018, Mel-021, and Mel-078 were derived from male, female, and male patients, respectively.
Once tumors reached 250 mm3 for Mel-077 and Mel-020, 500 mm3 for Mel-015, or 1,000 mm3 for Mel-006, the mice were enrolled into treatment cohorts. Mice from distinct batches were randomly assigned to the different experimental groups. Mice were treated daily by oral gavage with a capped dose of 600 dabrafenib µg–6 trametinib µg (DT), in 250 µl total volume, and tigecycline (50 mg/kg) was administered daily by i.p. injection.
Humanized huNOG-EXL (NOD.Cg-Prkdcscid Il2rgtm1Sug Tg(SV40/HTLV-IL3,CSF2)10-7Jic/JicTac), engrafted at the age of 6 wk with CD34+ hematopoietic stem cells from two different donors, were purchased from Taconic Biosciences. The percentage of humanization was 20–80% as estimated by FACS for CD45+ cells 10 wk after engraftment. The mice were engrafted with the Mel-006 human melanoma model described before at the age of 23 wk. When the tumor reached the size of 100 mm3, treatment with tigecycline (50 mg/kg i.p.) and/or nivolumab (10 mg/kg twice a week i.p.) was initiated. This experiment was approved by the KU Leuven animal ethical committee under ECD P210/2018.
No specific randomization method was used. According to animal welfare guidelines, mice were sacrificed when tumors reached a volume of 2,500 mm3 or when body weight decreased >20% from the initial weight. Mice used in this study never reached or exceeded these limits.
Allografts
5 × 104 YUMMER 1.7 cells were injected subcutaneously in the interscapular fat pad of 4-wk-old C57BL/6 males (males were chosen since females frequently rejected the implantation and had sporadic tumor ulceration). Mice were maintained under standard housing conditions with continuous access to food and water. The health and welfare of the animals was supervised by a designated veterinarian. The KU Leuven animal facilities comply with all appropriate standards (cages, space per animal, temperature (22°C), light, humidity, food, water), and all cages are enriched with materials that allow the animals to exert their natural behavior. Mice used in the study were maintained on a diurnal 12 h light/dark cycle. The experiments were approved by the KU Leuven animal ethical committee under ECD P049-2019 and performed in accordance with the internal, national, and European guidelines of animal care and use. Once the tumors reached 50 mm3 (for YUMMER 1.7), mice were enrolled in the treatment cohort. Mice from distinct batches were randomly assigned to the different experimental groups. Mice were treated with tigecycline (Bio Connect) i.p., 50 mg/kg daily, with anti–PD-1 (Ultra-LEAF Purified anti-mouse CD279 antibody, clone RMP1-14) i.p., 10 mg/kg twice per week, for 3 wk, for a total of six injections.
No specific randomization method was used. According to animal welfare guidelines, mice were sacrificed when tumors reach a volume of 2,500 mm3 or when body weight decreased >20% from the initial weight. Mice used in this paper never reached nor exceeded these limits.
Immunofluorescence on PDX biopsies
Fluorescent staining was performed using Opal staining reagents, which use individual tyramide signal amplification–conjugated fluorophores to detect various targets within an immunofluorescence assay. In brief, samples were fixed with 4% paraformaldehyde and embedded in paraffin. Serially cut sections of 5 µm were stained with H&E for routine light microscopy and used for immunohistochemistry.
Depending on the antibody, antigen retrieval (AR) was performed in citrate buffer at pH 6.0 or EDTA buffer at pH 9. Deparaffinized sections were then incubated overnight with primary antibodies against AQP1 (#AB2219; Millipore, 1:3,000), NGFR (8238; Cell Signaling Technology, 1:1,000), MITF (#HPA003259; Sigma-Aldrich, 1:100), MLANA (#HPA048662; Sigma-Aldrich, 1:200), CD36 (#HPA002018; Sigma-Aldrich, 1:200), AXL (#AF154; R&D Systems, 1:50), and S100 (#Z0311; Dako, 1:100). Subsequently, the slides were washed in PBS, pH 7.2, and incubated for 10 min at room temperature with Opal Polymer HRP Mouse Plus Rabbit secondaries (PerkinElmer). After another wash in PBS, the slides were incubated at room temperature for 10 min with one of the following Alexa Fluorescent tyramides (PerkinElmer) included in the Opal 4 color kit (NEL810001KT; Akoya Biosciences) to detect antibody staining, prepared according to the manufacturer’s instructions: Opal 520, Opal 570, and Opal 690 (dilution 1:50). Stripping of primary and secondary antibodies was performed by placing the slides in a plastic container filled with AR; microwave technology was used to bring the liquid to 100°C (2 min), and the sections were then microwaved for an additional 15 min at 75°C. Slides were allowed to cool in the AR buffer for 15 min at room temperature and were then rinsed with deionized water and 1× Tris-buffered saline with Tween 20. After three additional washes in deionized water, the slides were counterstained with DAPI for 5 min and mounted with ProLong Gold Antifade Mountant (P10144; Thermo Fisher Scientific). Slides were scanned for image acquisition using Zeiss AxioScan Z.1 and ZEN2 software. Images were quantified using QuPath. In brief, three different regions per section were quantified. Every dot in the graph represents the percentage of each quantified population. For MITF/NGFR/AQP1 Opal staining, melanocytic/proliferative cells were quantified based on MITFhigh expression, and only NGFR+/AQP1+ cells were considered to be NCSCs. In AXL/CD36/MLANA Opal staining, mesenchymal cells were quantified based on AXL positivity, and to identify SMCs, we quantified the presence of CD36+ cells.
Matrix-assisted laser desorption/ionization (MALDI)
Frozen cryosections were removed from the freezer, quickly placed into cryomolds, covered with a 10% solution of gelatin, and frozen in a 2-methylbutane/liquid nitrogen bath. Samples were then stored at −80°C until further processing. For cryosectioning, samples were taken from −80°C and placed into a cryostat (Leica CM1950) chamber for 1 h before cutting. Sections of 10-µm thickness were thaw-mounted onto indium tin oxide (ITO) glass (Ossila), washed subsequently with hexane and 2-propanol, and coated with a solution composed of 0.05% poly-L-lysine (Sigma-Aldrich) and 0.1% NP-40 detergent substitute. Sections were dried in a desiccator for 15 min, sealed in foil inside a plastic container, evacuated, and stored at −80°C until further processing. Samples were then transferred to room temperature, left to temper for 15 min, and dried in a desiccator for 15 min. Samples were then sprayed with a matrix solution of 9-aminoacridine (5 mg/ml in 70% ethanol [vol/vol]). A TM-Sprayer 3 (HTX-Technologies) was used to cover the tissue with matrix solution with the following settings: nozzle temperature 95°C, 16 passes, flow 0.06 ml/min, nozzle velocity 1,200 mm/min, track spacing 3 mm, CC pattern, gas pressure 10 psi, gas flow rate 2 liter/min, drying time 5 s, nozzle height 40 mm, and propelled with 50% methanol (vol/vol). After spraying, samples were dried in a desiccator for an additional 15 min.
Spectral images were taken with a rapifleX MALDI–time of flight (TOF)/TOF spectrometer (Bruker) in negative mode using the 355-nm smartbeam 3D laser with a spatial resolution of 50 × 50 µm within the range of m/z 20–1,100, with constant laser fluence of 44% and laser frequency of 5 kHz. 200 shots were accumulated from every position. Samples were measured in random order. Calibration was done externally using red phosphorus (Sigma-Aldrich), achieving precision ≤5 ppm.
For histology, tissues on ITO glass with matrix were washed three times with 70% ethanol, each time for 2 min. After washing in water, samples were stained with H&E, dehydrated through an ethanol series, washed in xylene (Leica ST5020), and mounted in Pertex (Leica CV5030). Samples were then scanned with Axio Scan.Z1 (Zeiss) operated by Zen 2 software (blue edition; Zeiss). Scanned images were processed with white balance corrected to background and exported as JPG files using ZEN 2.5 lite software (Zeiss). For image processing, single raw data files were processed at once using SCiLS Lab software (v2020a, SCiLS; Bruker). Raw data were smoothed with a convolution algorithm (width 20), and all subsequent computations and image rendering were done with root mean square normalization, with the exception of ATP/ADP ratio and lactate/pyruvate ratio, for which ADP or pyruvate raw peak values, respectively, were used to calculate the ratio. Regions of interest were corrected with H&E scans. Spectra were exported as CSV files and processed in mMass software (v5.5.0, Martin Strohalm). Peak picking was performed, and matrix-related peaks were removed. The resulting peak list was imported back to SCiLS Lab and used for computations.
Receiver operating characteristic analysis was performed to find discriminative values for two selected regions. All tissues from the respective groups (control, DT, or DTT) were used for computations. The operation was performed with all individual spectra, with a random subset of spectra set up properly according to the number of spectra in respective regions. Peaks with an area under the curve value >0.9 or <0.1 were considered significant. Based on morphological assessment, five region types were selected and marked by a trained pathologist: tumor tissue, fibrotic tissue, necrosis, white adipose tissue, and muscle tissue. Classification Pipeline was performed in SCiLS software to assign tissue type automatically. Resulting areas of tumor tissue were selected, and receiver operating characteristic analysis was performed to find differences in metabolite expression. All tissues from the respective groups (control, DT, or DTT) were used for computations. The operation was performed with all individual spectra with a random subset of spectra set up properly accordingly to the number of spectra in respective regions. Peaks with an area under the curve value above 0.9 or less than 0.1, respectively, were considered as significant.
Identification was done using modified HMDB v4 database (using compounds from serum, cerebrospinal fluid, urine, sweat, and saliva, excluding lipids, drugs, and metabolites; https://hmdb.ca/) considering [M-H]−, [M-H-H2O]−, [M + Cl]−, [M + Na-2H]−, [M + Na-2H-H2O]−, [M + K-2H]−, and [M + K-2H-H2O]− adducts within a tolerance of 30 ppm.
Case report
The patient was treated with BRAF-MEK inhibitors for stage IV malignant melanoma according to standard of care. PET-CT scans were done at baseline and for confirmation of remission. Every 2 mo, the patient underwent radiographical assessment with either PET-CT or regular CT scan.
Statistical analyses
All data are graphed as mean ± SEM. In animal experiments, n indicates the number of animals used. When comparing two treatment cohorts over time, statistical significance was calculated by two-way ANOVA. When comparing three or more groups, statistical significance was calculated by one-way ANOVA with Dunnett’s multiple comparisons test. When comparing Kaplan–Meier plots, statistical significance was calculated by log-rank (Mantel–Cox) test. Statistical tests used are specified in the corresponding figure legends. Throughout all figures, statistical significance was considered at P < 0.05 (*, P < 0.05; **, P < 0.01; ***, P < 0.001; and ****, P < 0.0001). All statistical analyses were performed with GraphPad Prism v8.4.2 (464), April 8, 2020, for Mac OS Catalina.
Online supplemental material
Fig. S1 shows the increase in mitochondrial translation upon activation of the ISR as seen through sequencing of a BRAF mutant cell line and confirmed by qPCR in a NRAS mutant cell line. Fig. S2 shows that tetracyclines affect multiple drug–tolerant states and thus increase OS in PDX models. In Fig. S3, a spatial metabolomic experiment shows that addition of tigecycline to targeted therapy does not significantly rewire melanoma metabolism in vivo. Fig. S4 shows that tigecycline significantly affects the growth of immune and targeted therapy–resistant cells and lesions. Fig. S5 shows that the effect of tetracyclines can be rescued by knocking down ATF4. Table S1 lists primers used in the study.
Acknowledgments
The authors wish to thank Ghanem E. Ghanem (Jules Bordet Institute, Brussels, Belgium), Marcus Bosenberg (Yale University, New Haven, CT), Irwin Davidson (Institute of Genetic and Molecular and Cellular Biology, Strasbourg, France), Aart Jochemsen (Leiden University, Leiden, Netherlands), and Göran Jönsson (Lund University, Lund, Sweden) for kindly making gifts of some of the melanoma cell lines used in this study. We also thank Amanda Lund for comments and suggestions. We are grateful to Petr Makovicky and Ashkan Zareie (Czech Centre for Phenogenomics) for their help with the MALDI-TOF analysis.
The present study was funded by Kom op tegen Kanker (Stand Up To Cancer), the Flemish Cancer Society, and by a KU Leuven C1 grant. E. Leucci is funded by the Melanoma Research Alliance Amanda and Jonathan Eilian young investigator award (2018). Trace staff are supported by Stichting Tegen Kanker grant 2016-054. V. Katopodi is a recipient of a Research Foundation – Flanders PhD fellowship (1S47519N). S. Cinque is a recipient of a Research Foundation – Flanders PhD fellowship (1SD1620N). The Czech Centre for Phenogenomics and Trace are members of the EuroPDX Consortium. L. Kucer and J. Rozman are supported by the Czech Academy of Sciences (RVO 68378050, LM2018126), the Ministry of Education, Youth and Sports (project CZ.02.1.01/0.0/0.0/16_013/0001789), and European Structural and Investment Funds.
Author contributions: R. Vendramin designed and performed in vitro and allograft experiments and supervised the work of A. Konnova and Z. Knezevic. A. Konnova performed cell treatments (tigecycline, doxycycline, and ISRIB) followed by colony assays and Incucyte measurements and collected the samples for western blots in Figs. 2 C and S5 D. Z. Knezevic performed the drug treatments in combination with ATF4 silencing and collected the samples for western blots in Fig. S5 F. R. Vendramin, A. Konnova, and Z. Knezevic performed allograft experiments. S. Cinque and V. Katopodi performed the in vitro experiments in Fig. 1 and ran western blots. S. Adnane performed the UM PDX experiment, contributed to the in vitro experiments, and extracted proteins from tumors. Y. Verheyden performed bioinformatic analysis. P. Karras and F.M. Bosisio performed and interpreted multiplex IHC. E. Demesmaeker and L. Rizzotto performed and designed cutaneous melanoma PDX experiments. I. Gladwyn-Ng assisted with the design of the humanized PDX experiments. L. Kucera and J. Rozman performed and analyzed spatial metabolomics. E. Dassi and S. Millevoi designed and analyzed rG4-related experiments. O. Bechter provided and analyzed the clinical case. E. Leucci and R. Vendramin analyzed and interpreted the data. E. Leucci acquired funding. J-C. Marine and E. Leucci designed the research study and wrote the manuscript with input from all the authors.
References
Competing Interests
Disclosures: R. Vendramin and E. Leucci reported a patent to PCT/EP2020/059416 pending "UVEAL MELANOMA TREATMENT." E. Leucci reported non-financial support from Taconic Biosciences during the conduct of the study. I. Gladwyn-Ng is an employee of Taconic Biosciences. No other disclosures were reported.
Author notes
J-C. Marine and E. Leucci contributed equally to this paper.