Reconstructing the evolutionary history of multiple myeloma

Timing antigenic escape in multiple myeloma treated with T-cell redirecting immunotherapies

Recent data highlight genomic events driving antigen escape as a recurring cause of chimeric antigen receptor T-cell (CAR-T) and bispecific T-cell engager (TCE) resistance in multiple myeloma (MM). Yet, it remains unclear if these events, leading to clonal dominance at progression, result from acquisition under treatment selection or selection of pre-existing undetectable clones. This differentiation gains importance as these immunotherapies progress to earlier lines of treatment, prompting the need for innovative diagnostic testing to detect these events early on. By reconstructing phylogenetic trees and exploring chemotherapy mutational signatures as temporal barcodes in 11 relapsed refractory MM patients with available whole genome sequencing data before and after CART/TCE treatment, we demonstrated that somatic antigen escape mechanisms for BCMA- and GPRC5D-targeting therapies are acquired post-diagnosis, likely during CART/TCE treatment. Longitudinal tracking of these mutations using digital PCR in 4 patients consistently showed that genomic events promoting antigen escape were not detectable during the initial months of therapy but began to emerge nearly 1 year post therapy initiation. This finding reduces the necessity for a diagnostic panel to identify these events before CART/TCE. Instead, it underscores the importance of surveillance and identifying patients at higher risk of acquiring these events.

G-Quadruplex Forming DNA Sequence Context Is Enriched around Points of Somatic Mutations in a Subset of Multiple Myeloma Patients

Multiple myeloma (MM) is the second most common hematological malignancy, which remains incurable despite recent advances in treatment strategies. Like other forms of cancer, MM is characterized by genomic instability, caused by defects in DNA repair. Along with mutations in DNA repair genes and genotoxic drugs used to treat MM, non-canonical secondary DNA structures (four-stranded G-quadruplex structures) can affect accumulation of somatic mutations and chromosomal abnormalities in the tumor cells of MM patients. Here, we tested the hypothesis that G-quadruplex structures may influence the distribution of somatic mutations in the tumor cells of MM patients. We sequenced exomes of normal and tumor cells of 11 MM patients and analyzed the data for the presence of G4 context around points of somatic mutations. To identify molecular mechanisms that could affect mutational profile of tumors, we also analyzed mutational signatures in tumor cells as well as germline mutations for the presence of specific SNPs in DNA repair genes or in genes regulating G-quadruplex unwinding. In several patients, we found that sites of somatic mutations are frequently located in regions with G4 context. This pattern correlated with specific germline variants found in these patients. We discuss the possible implications of these variants for mutation accumulation and specificity in MM and propose that the extent of G4 context enrichment around somatic mutation sites may be a novel metric characterizing mutational processes in tumors.

The genomic landscape of Vk*MYC myeloma highlights shared pathways of transformation between mice and humans

Multiple myeloma (MM) is a heterogeneous disease characterized by frequent MYC translocations. Sporadic MYC activation in the germinal center of genetically engineered Vk*MYC mice is sufficient to induce plasma cell tumors in which a variety of secondary mutations are spontaneously acquired and selected over time. Analysis of 119 Vk*MYC myeloma reveals recurrent copy number alterations, structural variations, chromothripsis, driver mutations, apolipoprotein B mRNA-editing enzyme, catalytic polypeptide (APOBEC) mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identify frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC. In summary, here we credential the Vk*MYC mouse as a unique resource to explore MM genomic evolution and describe a fully annotated collection of diverse and immortalized murine MM tumors.

Approach to High-Risk Multiple Myeloma

Improving the outcome of high-risk myeloma (HRMM) is a key therapeutic aim for the next decade. To achieve this aim, it is necessary to understand in detail the genetic drivers underlying this clinical behavior and to target its biology therapeutically. Advances have already been made, with a focus on consensus guidance and the application of novel immunotherapeutic approaches. Cases of HRMM are likely to have impaired prognosis even with novel strategies. However, if disease eradication and minimal disease states are achieved, then cure may be possible.

Does the Simultaneous Introduction of Several Pharmaceuticals in the Post-Lenalidomide Era Translate to Better Outcomes in Relapse Refractory Multiple Myeloma? Findings from the Real-World Innovation in Multiple Myeloma (REAL IMM) Study

Newer methodologies are needed to assess the real-world comparative effectiveness of a "generation" of pharmaceutical innovation versus the prior standard of care. This chart review study aimed to first evaluate the cumulative clinical benefits of pharmaceutical innovation in everyday relapse/refractory multiple myeloma before analyzing findings in the context of respective real-world outcomes from the bortezomib/lenalidomide era. Study endpoints included the 52-week PFS rate in second and third line of therapy (LOT), mPFS-2 across the first and second LOT, the ORR, reasons for discontinuation, and the treatment duration per therapeutic algorithm. Data from 107 patients were collected. The median follow-up was 2.0 years. Of the subjects who met the selection criteria for the second LOT, 72.2% maintained the PFS at 52 weeks. In the third-line setting, the PFS rate at 52 weeks was 63.5%. The mPFS across the first and second, the second, and the third LOTs were 26, 17, and 15 months, respectively. The ORR was 76.1% in the second and 69.7% in the third LOT. After non-response or progression, the main reason for drug discontinuation was treatment intolerability. The second-line ORR and the 52-week PFS rate were similar to previous real-world findings from the bortezomib/lenalidomide era. The cumulative mPFS across the second and third LOTs was higher than the respective mPFS across the first and second LOTs. Despite its limitations, the methodology and findings from this study may be used in future clinical and economic evaluations across all hematological malignancies.

Genomic and immune signatures predict clinical outcome in newly diagnosed multiple myeloma treated with immunotherapy regimens

Despite improving outcomes, 40% of patients with newly diagnosed multiple myeloma treated with regimens containing daratumumab, a CD38-targeted monoclonal antibody, progress prematurely. By integrating tumor whole-genome and microenvironment single-cell RNA sequencing from upfront phase 2 trials using carfilzomib, lenalidomide and dexamethasone with daratumumab ( NCT03290950 ), we show how distinct genomic drivers including high APOBEC mutational activity, IKZF3 and RPL5 deletions and 8q gain affect clinical outcomes. Furthermore, evaluation of paired bone marrow profiles, taken before and after eight cycles of carfilzomib, lenalidomide and dexamethasone with daratumumab, shows that numbers of natural killer cells before treatment, high T cell receptor diversity before treatment, the disappearance of sustained immune activation (that is, B cells and T cells) and monocyte expansion over time are all predictive of sustained minimal residual disease negativity. Overall, this study provides strong evidence of a complex interplay between tumor cells and the immune microenvironment that is predictive of clinical outcome and depth of treatment response in patients with newly diagnosed multiple myeloma treated with highly effective combinations containing anti-CD38 antibodies.

The Vk*MYC Mouse Model recapitulates human multiple myeloma evolution and genomic diversity

Despite advancements in profiling multiple myeloma (MM) and its precursor conditions, there is limited information on mechanisms underlying disease progression. Clincal efforts designed to deconvolute such mechanisms are challenged by the long lead time between monoclonal gammopathy and its transformation to MM. MM mouse models represent an opportunity to overcome this temporal limitation. Here, we profile the genomic landscape of 118 genetically engineered Vk*MYC MM and reveal that it recapitulates the genomic heterogenenity and life history of human MM. We observed recurrent copy number alterations, structural variations, chromothripsis, driver mutations, APOBEC mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identified frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC expression, that drives the progression of monoclonal gammopathy to MM.

Molecular Diagnostics of Plasma Cell Neoplasms

Genetic characterization of myeloma at diagnosis by interphase fluorescence in situ hybridization and next-generation sequencing (NGS) can assist with risk stratification and treatment planning. Measurable residual disease (MRD) status after treatment, as evaluated by next-generation flow cytometry or NGS on bone marrow aspirate material, is one of the most important predictors of prognosis. Less-invasive tools for MRD assessment such as liquid biopsy approaches have also recently emerged as potential alternatives.

LAMP5 may promote MM progression by activating p38

Multiple myeloma (MM) is the second most common tumor of the hematologic system. MM remains incurable at this time. In this study, we used bioinformatics analysis to find key genes in the pathogenesis of MM. We first found that Lysosome associated membrane protein 5 (LAMP5) expression was sequentially increased in healthy donors (HD), monoclonal gammopathy of undetermined significance (MGUS), smoldering multiple myeloma (SMM) and newly diagnosed MM (NDMM), relapsed MM (RMM). We collected bone marrow from patients with NDMM, HD and post-treatment MM (PTMM) and performed qPCR analysis of LAMP5, and found that the expression of LAMP5 is stronger in NDMM than in HD, and decreases after treatment. Western blotting assay also found more expression of LAMP5 in NDMM than in HD. Patients with high LAMP5 expression have a higher DS (Durie-Salmon) stage and worse prognosis. We next verified the expression of LAMP5 in four MM cell lines and silenced LAMP5 expression in RPMI-8226 and AMO-1, and explored the effects of LAMP5 silencing on MM cell apoptosis and cell cycle by flow cytometry and western blotting. Knockdown of LAMP5 promoted apoptosis in MM cells, but had no effect on the cell cycle. Mechanistically, LAMP5 may exert its pro-tumor effects in MM in part through activation of p38 protein. We screened LAMP5 for the first time as a key gene for MM progression and recurrence, and found that LAMP5 may exert its pro-tumor effects in MM through activation of p38 protein.

Predictive Value Analysis of Serum Ig A, Ig G, and TNF-α in Recurrence of Multiple Myeloma

Objective. The study is aimed at analyzing the predictive value of serum Ig A, Ig G, and TNF-α in the recurrence of multiple myeloma (MM). Methods. 136 patients with MM treated in our hospital from January 2010 to January 2017 were followed up for 5 years. Finally, 100 patients who met the inclusion and exclusion criteria and had the complete follow-up visit were selected as the study subjects, with the recurrence of MM as endpoint event, and the observation was taken until the occurrence of endpoint event in patients or the termination of this study. They were divided into the recurrence group (RG) and the nonrecurrence group (NRG) according to whether the endpoint event occurred. The venous blood of patients was collected at the first diagnosis and subsequent visit (at the time of recurrence or termination of the study) to measure the Ig A and Ig G using a full automatic special protein analyzer and the TNF-α level by enzyme-linked immunosorbent assay. The data obtained in this study were analyzed by univariate analysis to choose the factors with difference in statistical significance to draw the ROC curve, and the areas under the curve (AUC) were recorded to analyze the potential mechanism of Ig A, Ig G, and TNF-α in predicting the recurrence of MM. Results. After follow-up visit, there were 62 patients with recurrence (62.0%) and 38 patients without recurrence (38.0%), with no obvious difference in gender, age, body weight, and immune classification between the two groups ( $P>0.05$ ). Compared with the NRG, the levels of soluble interleukin-2 receptor (sIL-2R) and β2-microglobulin (β2-MG) in the RG at the first diagnosis were distinctly higher ( $P<0.001$ ); the levels of Ig A, Ig G, and TNF-α in the RG at the first diagnosis were visibly higher ( $P<0.05$ ); and the levels of Ig A, Ig G, and TNF-α in the RG at the subsequent visit were clearly higher ( $P<0.05$ ). There was a correlation between Ig G, Ig A, and TNF-α and β2-MG at the first diagnosis and the subsequent visit ( $P<0.05$ ); there was a correlation between Ig G and TNF-α, and sIL-2R at the first diagnosis and the subsequent visit ( $P<0.05$ ); and there was a correlation between Ig A and sIL-2R at the subsequent visit ( $P<0.05$ ). The AUC of Ig G, Ig A, and TNF-α in predicting the MM at the first diagnosis were 0.772, 0.776, and 0.778, respectively. Conclusion. The serum Ig A, Ig G, and TNF-α had a predictive value in the recurrence of MM, and TNF-α was correlated with sIL-2R and β2-MG, with the highest AUC and the best predictive value.

Correlation of changes in subclonal architecture with progression in the MMRF CoMMpass study

Multiple myeloma (MM) is a heterogeneous plasma cell proliferative disorder that arises from its premalignant precursor stages through a complex cascade of interactions between clonal mutations and co-evolving microenvironment. The temporo-spatial evolutionary trajectories of MM are established early during myelomatogenesis in precursor stages and retained in MM. Such molecular events impact subsequent disease progression and clinical outcomes. Identification of clonal sweeps of actionable gene targets in MM could reveal potential vulnerabilities that may exist in early stages and thus potentiate prognostication and customization of early therapeutic interventions. We have evaluated clonal evolution at multiple time points in 76 MM patients enrolled in the MMRF CoMMpass study. The major findings of this study are (a) MM progresses predominantly through branching evolution, (b) there is a heterogeneous spectrum of mutational landscapes that include unique actionable gene targets at diagnosis compared to progression, (c) unique clonal gains/ losses of mutant driver genes can be identified in patients with different cytogenetic aberrations, (d) there is a significant correlation between co-occurring oncogenic mutations/ co-occurring subclones e.g., with mutated TP53+SYNE1, NRAS+MAGI3, and anticorrelative dependencies between FAT3+FCGBP gene pairs. Such co-trajectories may synchronize molecular events of drug response, myelomatogenesis and warrant future studies to explore their potential for early prognostication and development of risk stratified personalized therapies in MM.

Kinase Inhibition in Multiple Myeloma: Current Scenario and Clinical Perspectives

Multiple myeloma (MM) is a blood cell neoplasm characterized by excessive production of malignant monoclonal plasma cells (activated B lymphocytes) by the bone marrow, which end up synthesizing antibodies or antibody fragments, called M proteins, in excess. The accumulation of this production, both cells themselves and of the immunoglobulins, causes a series of problems for the patient, of a systemic and local nature, such as blood hyperviscosity, renal failure, anemia, bone lesions, and infections due to compromised immunity. MM is the third most common hematological neoplasm, constituting 1% of all cancer cases, and is a disease that is difficult to treat, still being considered an incurable disease. The treatments currently available cannot cure the patient, but only extend their lifespan, and the main and most effective alternative is autologous hematopoietic stem cell transplantation, but not every patient is eligible, often due to age and pre-existing comorbidities. In this context, the search for new therapies that can bring better results to patients is of utmost importance. Protein tyrosine kinases (PTKs) are involved in several biological processes, such as cell growth regulation and proliferation, thus, mutations that affect their functionality can have a great impact on crucial molecular pathways in the cells, leading to tumorigenesis. In the past couple of decades, the use of small-molecule inhibitors, which include tyrosine kinase inhibitors (TKIs), has been a hallmark in the treatment of hematological malignancies, and MM patients may also benefit from TKI-based treatment strategies. In this review, we seek to understand the applicability of TKIs used in MM clinical trials in the last 10 years.

Genome Instability in Multiple Myeloma: Facts and Factors

Multiple myeloma (MM) is a malignant neoplasm of terminally differentiated immunoglobulin-producing B lymphocytes called plasma cells. MM is the second most common hematologic malignancy, and it poses a heavy economic and social burden because it remains incurable and confers a profound disability to patients. Despite current progress in MM treatment, the disease invariably recurs, even after the transplantation of autologous hematopoietic stem cells (ASCT). Biological processes leading to a pathological myeloma clone and the mechanisms of further evolution of the disease are far from complete understanding. Genetically, MM is a complex disease that demonstrates a high level of heterogeneity. Myeloma genomes carry numerous genetic changes, including structural genome variations and chromosomal gains and losses, and these changes occur in combinations with point mutations affecting various cellular pathways, including genome maintenance. MM genome instability in its extreme is manifested in mutation kataegis and complex genomic rearrangements: chromothripsis, templated insertions, and chromoplexy. Chemotherapeutic agents used to treat MM add another level of complexity because many of them exacerbate genome instability. Genome abnormalities are driver events and deciphering their mechanisms will help understand the causes of MM and play a pivotal role in developing new therapies.

Branching clonal evolution patterns predominate mutational landscape in multiple myeloma

Multiple Myeloma (MM) arises from malignant transformation and deregulated proliferation of clonal plasma cells (PCs) harbouring heterogeneous molecular anomalies. The effect of evolving mutations on clone fitness and their cellular prevalence shapes the progressing myeloma genome and impacts clinical outcomes. Although clonal heterogeneity in MM is well established, which subclonal mutations emerge/persist/perish with progression in MM and which of these can be targeted therapeutically remains an open question. In line with this, we have sequenced pairwise whole exomes of 62 MM patients collected at two time points, i.e., at diagnosis and on progression. Somatic variants were called using a novel ensemble approach where a consensus was deduced from four variant callers (Illumina's Dragen, Strelka2, SomaticSniper and SpeedSeq) and actionable/druggable gene targets were identified. A marked intraclonal heterogeneity was observed. Branching evolution was observed among 72.58% patients, of whom 64.51% had low TMBs (<10) and 61.29% had 2 or more founder clones. The hypermutator patients (with high TMB levels ≥10 to ≤100) showed a significant decrease in their TMBs from diagnosis (median TMB 77.11) to progression (median TMB 31.22). A distinct temporal fall in subclonal driver mutations was identified recurrently across diagnosis to progression e.g., in PABPC1, BRAF, KRAS, CR1, DIS3 and ATM genes in 3 or more patients suggesting such patients could be treated early with target specific drugs like Vemurafenib/Cobimetinib. An analogous rise in driver mutations was observed in KMT2C, FOXD4L1, SP140, NRAS and other genes. A few drivers such as FAT4, IGLL5 and CDKN1A retained consistent distribution patterns at two time points. These findings are clinically relevant and point at consideration of evaluating multi time point subclonal mutational landscapes for designing better risk stratification strategies and tailoring time to time risk adapted combination therapies in future.

Tumor-associated macrophages modulate angiogenesis and tumor growth in a xenograft mouse model of multiple myeloma

Tumor-associated macrophages (TAMs) are closely associated with poor multiple myeloma (MM) prognosis. Therefore, in-depth understanding of the mechanism by which TAM supports MM progression may lead to its effective treatment. We used the MM nude mouse subcutaneous xenograft model to evaluate the efficacy of the macrophage-depleting agent clodronate liposome (Clo) against MM and elucidate the mode of action of this therapy. At the same time, observe whether the elimination of TAM in vivo while silencing the expression of VEGFA has the same effect as in vitro experiments. We also used Clo to eliminate macrophages and reinjected M1 or M2 TAM through mouse tail veins to investigate the effects of various macrophage subtypes on MM xenograft tumor growth. We applied qRT-PCR, immunohistochemistry, and enzyme-linked immunosorbent assay to quantify VEGFA, CD31, and CD163 expression in tumor tissues and sera. Removal of TAMs from the tumor microenvironment impeded tumor growth. The combination of Clo plus VEGFA siRNA had a stronger inhibitory effect on tumor growth than Clo alone, and M2 and M1 macrophages promoted and inhibited tumor growth, respectively. Macrophage depletion combined with cytokine blocking is a promising MM treatment. Targeted M2 macrophage elimination together with cytokine block may be more effective at inhibiting MM growth than either treatment alone. The results of the present study lay an empirical foundation for the development of novel therapeutic strategies for MM.

Genomics of Smoldering Multiple Myeloma: Time for Clinical Translation of Findings?

Smoldering multiple myeloma (SMM) is an asymptomatic disorder of clonal bone marrow (BM) plasma cells (PCs) in between the premalignant condition known as monoclonal gammopathy of undetermined significance and overt multiple myeloma (MM). It is characterized by a deep biological heterogeneity that is reflected in a markedly variable progression risk among patients. Recently proposed risk stratification models mainly rely on indirect markers of disease burden and are unable to identify cases in whom clonal PCs have already undergone the "malignant switch" but major clonal expansion has not occurred yet. In the last years, the application of next-generation sequencing (NGS) techniques has led to profound advances in the understanding of the molecular bases of SMM progression, and in all likelihood, it will contribute to the needed improvement of SMM prognostication. In this Review, we describe the recent advances in characterizing the genomic landscape of SMM and intrinsic determinants of its progression, highlighting their implications in terms of understanding of tumor evolution and prognostication. We also review the main studies investigating the role of the microenvironment in this early disease stage. Finally, we mention the results of the first randomized clinical trials and discuss the potential clinical translability of the genomic insights.

Genome-wide identification of potential biomarkers in multiple myeloma using meta-analysis of mRNA and miRNA expression data

Multiple myeloma (MM) is a plasma cell malignancy with diverse clinical phenotypes and molecular heterogeneity not completely understood. Differentially expressed genes (DEGs) and miRNAs (DEMs) in MM may influence disease pathogenesis, clinical presentation / drug sensitivities. But these signatures overlap meagrely plausibly due to complexity of myeloma genome, diversity in primary cells studied, molecular technologies/ analytical tools utilized. This warrants further investigations since DEGs/DEMs can impact clinical outcomes and guide personalized therapy. We have conducted genome-wide meta-analysis of DEGs/DEMs in MM versus Normal Plasma Cells (NPCs) and derived unified putative signatures for MM. 100 DEMs and 1,362 DEGs were found deranged between MM and NPCs. Signatures of 37 DEMs ('Union 37') and 154 DEGs ('Union 154') were deduced that shared 17 DEMs and 22 DEGs with published prognostic signatures, respectively. Two miRs (miR-16-2-3p, 30d-2-3p) correlated with survival outcomes. PPI analysis identified 5 topmost functionally connected hub genes (UBC, ITGA4, HSP90AB1, VCAM1, VCP). Transcription factor regulatory networks were determined for five seed DEGs with ≥ 4 biomarker applications (CDKN1A, CDKN2A, MMP9, IGF1, MKI67) and three topmost up/ down regulated DEMs (miR-23b, 195, let7b/ miR-20a, 155, 92a). Further studies are warranted to establish and translate prognostic potential of these signatures for MM.

Whole-genome sequencing reveals progressive versus stable myeloma precursor conditions as two distinct entities

Multiple myeloma (MM) is consistently preceded by precursor conditions recognized clinically as monoclonal gammopathy of undetermined significance (MGUS) or smoldering myeloma (SMM). We interrogate the whole genome sequence (WGS) profile of 18 MGUS and compare them with those from 14 SMMs and 80 MMs. We show that cases with a non-progressing, clinically stable myeloma precursor condition (n = 15) are characterized by later initiation in the patient's life and by the absence of myeloma defining genomic events including: chromothripsis, templated insertions, mutations in driver genes, aneuploidy, and canonical APOBEC mutational activity. This data provides evidence that WGS can be used to recognize two biologically and clinically distinct myeloma precursor entities that are either progressive or stable.

Initial Whole-Genome Sequencing of Plasma Cell Neoplasms in First Responders and Recovery Workers Exposed to the World Trade Center Attack of September 11, 2001

PURPOSE

The World Trade Center (WTC) attack of September 11, 2001 created an unprecedented environmental exposure to known and suspected carcinogens. High incidence of multiple myeloma and precursor conditions has been reported among first responders to the WTC disaster. To expand on our prior screening studies, and to characterize the genomic impact of the exposure to known and potential carcinogens in the WTC debris, we were motivated to perform whole-genome sequencing (WGS) of WTC first responders and recovery workers who developed a plasma cell disorder after the attack.

EXPERIMENTAL DESIGN

We performed WGS of nine CD138-positive bone marrow mononuclear samples from patients who were diagnosed with plasma cell disorders after the WTC disaster.

RESULTS

No significant differences were observed in comparing the post-WTC driver and mutational signature landscapes with 110 previously published WGSs from 56 patients with multiple myeloma and the CoMMpass WGS cohort (n = 752). Leveraging constant activity of the single-base substitution mutational signatures 1 and 5 over time, we estimated that tumor-initiating chromosomal gains were windowed to both pre- and post-WTC exposure.

CONCLUSIONS

Although limitations in sample size preclude any definitive conclusions, our findings suggest that the observed increased incidence of plasma cell neoplasms in this population is due to complex and heterogeneous effects of the WTC exposure that may have initiated or contributed to progression of malignancy.

Mouse tumor susceptibility genes identify drug combinations for multiple myeloma

Long-term genetic studies utilizing backcross and congenic strain analyses coupled with positional cloning strategies and functional studies identified Cdkn2a, Mtor, and Mndal as mouse plasmacytoma susceptibility/resistance genes. Tumor incidence data in congenic strains carrying the resistance alleles of Cdkn2a and Mtor led us to hypothesize that drug combinations affecting these pathways are likely to have an additive, if not synergistic effect in inhibiting tumor cell growth. Traditional and novel systems-level genomic approaches were used to assess combination activity, disease specificity, and clinical potential of a drug combination involving rapamycin/everolimus, an Mtor inhibitor, with entinostat, an histone deacetylase inhibitor. The combination synergistically repressed oncogenic MYC and activated the Cdkn2a tumor suppressor. The identification of MYC as a primary upstream regulator led to the identification of small molecule binders of the G-quadruplex structure that forms in the NHEIII region of the MYC promoter. These studies highlight the importance of identifying drug combinations which simultaneously upregulate tumor suppressors and downregulate oncogenes.