Arid1a mutation suppresses TGF-ß signaling and induces cholangiocarcinoma.

Activating KRAS mutations and functional loss of members of the SWI/SNF complex, including ARID1A, are found together in the primary liver tumor cholangiocarcinoma (CC). How these mutations cooperate to promote CC has not been established. Using murine models of hepatocyte and biliary-specific lineage tracing, we show that Kras and Arid1a mutations drive the formation of CC and tumor precursors from the biliary compartment, which are accelerated by liver inflammation. Using cultured cells, we find that Arid1a loss causes cellular proliferation, escape from cell-cycle control, senescence, and widespread changes in chromatin structure. Notably, we show that the biliary proliferative response elicited by Kras/Arid1a cooperation and tissue injury in CC is caused by failed engagement of the TGF-ß-Smad4 tumor suppressor pathway. We thus identify an ARID1A-TGF-ß-Smad4 axis as essential in limiting the biliary epithelial response to oncogenic insults, while its loss leads to biliary pre-neoplasia and CC.

Combined Inhibition of G9a and EZH2 Suppresses Tumor Growth via Synergistic Induction of IL24-Mediated Apoptosis.

G9a and EZH2 are two histone methyltransferases commonly upregulated in several cancer types, yet the precise roles that these enzymes play cooperatively in cancer is unclear. We demonstrate here that frequent concurrent upregulation of both G9a and EZH2 occurs in several human tumors. These methyltransferases cooperatively repressed molecular pathways responsible for tumor cell death. In genetically distinct tumor subtypes, concomitant inhibition of G9a and EZH2 potently induced tumor cell death, highlighting the existence of tumor cell survival dependency at the epigenetic level. G9a and EZH2 synergistically repressed expression of genes involved in the induction of endoplasmic reticulum (ER) stress and the production of reactive oxygen species. IL24 was essential for the induction of tumor cell death and was identified as a common target of G9a and EZH2. Loss of function of G9a and EZH2 activated the IL24-ER stress axis and increased apoptosis in cancer cells while not affecting normal cells. These results indicate that G9a and EZH2 promotes the evasion of ER stress-mediated apoptosis by repressing IL24 transcription, therefore suggesting that their inhibition may represent a potential therapeutic strategy for solid cancers. SIGNIFICANCE: These findings demonstrate a novel role for G9a and EZH2 histone methyltransferases in suppressing apoptosis, which can be targeted with small molecule inhibitors as a potential approach to improve solid cancer treatment.

Subtype-Independent ANP32E Reduction During Breast Cancer Progression in Accordance with Chromatin Relaxation.

BACKGROUND: Chromatin state provides a clear decipherable blueprint for maintenance of transcriptional patterns, exemplifying a mitotically stable form of cellular programming in dividing cells. In this regard, genomic studies of chromatin states within cancerous tissues have the potential to uncover novel aspects of tumor biology and unique mechanisms associated with disease phenotypes and outcomes. The degree to which chromatin state differences occur in accordance with breast cancer features has not been established. METHODS: We applied a series of unsupervised computational methods to identify chromatin and molecular differences associated with discrete physiologies across human breast cancer tumors. RESULTS: Chromatin patterns alone are capable of stratifying tumors in association with cancer subtype and disease progression. Major differences occur at DNA motifs for the transcription factor FOXA1, in hormone receptor-positive tumors, and motifs for SOX9 in Basal-like tumors. We find that one potential driver of this effect, the histone chaperone ANP32E, is inversely correlated with tumor progression and relaxation of chromatin at FOXA1 binding sites. Tumors with high levels of ANP32E exhibit an immune response and proliferative gene expression signature, whereas tumors with low ANP32E levels appear programmed for differentiation. CONCLUSIONS: Our results indicate that ANP32E may function through chromatin state regulation to control breast cancer differentiation and tumor plasticity. This study sets a precedent for future computational studies of chromatin changes in carcinogenesis.

Ten-eleven translocation protein 1 modulates medulloblastoma progression.

BACKGROUND: Medulloblastoma (MB) is the most common malignant pediatric brain tumor that originates in the cerebellum and brainstem. Frequent somatic mutations and deregulated expression of epigenetic regulators in MB highlight the substantial role of epigenetic alterations. 5-hydroxymethylcytosine (5hmC) is a highly abundant cytosine modification in the developing cerebellum and is regulated by ten-eleven translocation (TET) enzymes. RESULTS: We investigate the alterations of 5hmC and TET enzymes in MB and their significance to cerebellar cancer formation. We show total abundance of 5hmC is reduced in MB, but identify significant enrichment of MB-specific 5hmC marks at regulatory regions of genes implicated in stem-like properties and Nanog-binding motifs. While TET1 and TET2 levels are high in MBs, only knockout of Tet1 in the smoothened (SmoA1) mouse model attenuates uncontrolled proliferation, leading to a favorable prognosis. The pharmacological Tet1 inhibition reduces cell viability and platelet-derived growth factor signaling pathway-associated genes. CONCLUSIONS: These results together suggest a potential key role of 5hmC and indicate an oncogenic nature for TET1 in MB tumorigenesis, suggesting it as a potential therapeutic target for MBs.

Chromatin Accessibility Identifies Regulatory Elements Predictive of Gene Expression and Disease Outcome in Multiple Myeloma.

PURPOSE: Multiple myeloma is a malignancy of plasma cells. Extensive genetic and transcriptional characterization of myeloma has identified subtypes with prognostic and therapeutic implications. In contrast, relatively little is known about the myeloma epigenome. EXPERIMENTAL DESIGN: CD138+CD38+ myeloma cells were isolated from fresh bone marrow aspirate or the same aspirate after freezing for 1-6 months. Gene expression and chromatin accessibility were compared between fresh and frozen samples by RNA sequencing (RNA-seq) and assay for transpose accessible chromatin sequencing (ATAC-seq). Chromatin accessible regions were used to identify regulatory RNA expression in more than 700 samples from newly diagnosed patients in the Multiple Myeloma Research Foundation CoMMpass trial (NCT01454297). RESULTS: Gene expression and chromatin accessibility of cryopreserved myeloma recapitulated that of freshly isolated samples. ATAC-seq performed on a series of biobanked specimens identified thousands of chromatin accessible regions with hundreds being highly coordinated with gene expression. More than 4,700 of these chromatin accessible regions were transcribed in newly diagnosed myelomas from the CoMMpass trial. Regulatory element activity alone recapitulated myeloma gene expression subtypes, and in particular myeloma subtypes with immunoglobulin heavy chain translocations were defined by transcription of distal regulatory elements. Moreover, enhancer activity predicted oncogene expression implicating gene regulatory mechanisms in aggressive myeloma. CONCLUSIONS: These data demonstrate the feasibility of using biobanked specimens for retrospective studies of the myeloma epigenome and illustrate the unique enhancer landscapes of myeloma subtypes that are coupled to gene expression and disease progression.

Acetylation of KLF5 maintains EMT and tumorigenicity to cause chemoresistant bone metastasis in prostate cancer.

Advanced prostate cancer (PCa) often develops bone metastasis, for which therapies are very limited and the underlying mechanisms are poorly understood. We report that bone-borne TGF-ß induces the acetylation of transcription factor KLF5 in PCa bone metastases, and acetylated KLF5 (Ac-KLF5) causes osteoclastogenesis and bone metastatic lesions by activating CXCR4, which leads to IL-11 secretion, and stimulating SHH/IL-6 paracrine signaling. While essential for maintaining the mesenchymal phenotype and tumorigenicity, Ac-KLF5 also causes resistance to docetaxel in tumors and bone metastases, which is overcome by targeting CXCR4 with FDA-approved plerixafor. Establishing a mechanism for bone metastasis and chemoresistance in PCa, these findings provide a rationale for treating chemoresistant bone metastasis of PCa with inhibitors of Ac-KLF5/CXCR4 signaling.

Aberrant Hypermethylation-Mediated Suppression of PYCARD Is Extremely Frequent in Prostate Cancer with Gleason Score ≥ 7.

Epigenetic gene silencing by aberrant DNA methylation leads to loss of key cellular pathways in tumorigenesis. In order to analyze the effects of DNA methylation on prostate cancer, we established LNCaP-derived human prostate cancer cells that can pharmacologically induce global reactivation of hypermethylated genes by the methyl-CpG targeted transcriptional activation (MeTA) method. The MeTA suppressed the growth of LNCaP-derived cells and induced apoptosis. Microarray analysis indicated that PYCARD (PYD and CARD domain containing) encoding an apoptosis-inducing factor was upregulated by 65-fold or more after treatment with MeTA. We analyzed DNA methylation statuses using 50 microdissected primary prostate cancer tissues and found an extremely high frequency of tumor-specific promoter hypermethylation of PYCARD (90%, 45/50). Moreover, DNA methylation status was significantly associated with Gleason score (P = 0.0063); the frequency of tumor-specific hypermethylation was 96% (44/46) in tumors with Gleason score ≥ 7, whereas that in tumors with Gleason score 6 was 25% (1/4). Immunohistochemical analyses using these 50 cases indicated that only 8% (4/50) of cancerous tissues expressed PYCARD, whereas 80% (40/50) of corresponding normal prostate epithelial and/or basal cells expressed PYCARD. In addition, there was no relationship between PYCARD immunostaining and the Gleason score in cancerous tissue and surrounding normal tissue. Inducible expression of PYCARD inhibited cell proliferation by induction of apoptosis. These results suggest that aberrant methylation of PYCARD is a distinctive feature of prostate cancers with Gleason score ≥ 7 and may play an important role in escaping from apoptosis in prostatic tumorigenesis.

Epigenetically heterogeneous tumor cells direct collective invasion through filopodia-driven fibronectin micropatterning.

Tumor heterogeneity drives disease progression, treatment resistance, and patient relapse, yet remains largely underexplored in invasion and metastasis. Here, we investigated heterogeneity within collective cancer invasion by integrating DNA methylation and gene expression analysis in rare purified lung cancer leader and follower cells. Our results showed global DNA methylation rewiring in leader cells and revealed the filopodial motor MYO10 as a critical gene at the intersection of epigenetic heterogeneity and three-dimensional (3D) collective invasion. We further identified JAG1 signaling as a previously unknown upstream activator of MYO10 expression in leader cells. Using live-cell imaging, we found that MYO10 drives filopodial persistence necessary for micropatterning extracellular fibronectin into linear tracks at the edge of 3D collective invasion exclusively in leaders. Our data fit a model where epigenetic heterogeneity and JAG1 signaling jointly drive collective cancer invasion through MYO10 up-regulation in epigenetically permissive leader cells, which induces filopodia dynamics necessary for linearized fibronectin micropatterning.

Prognostic significance of an invasive leader cell-derived mutation cluster on chromosome 16q.

BACKGROUND: Intratumoral heterogeneity is defined by subpopulations with varying genotypes and phenotypes. Specialized, highly invasive leader cells and less invasive follower cells are phenotypically distinct subpopulations that cooperate during collective cancer invasion. Because leader cells are a rare subpopulation that would be missed by bulk sequencing, a novel image-guided genomics platform was used to precisely select this subpopulation. This study identified a novel leader cell mutation signature and tested its ability to predict prognosis in non-small cell lung cancer (NSCLC) patient cohorts. METHODS: Spatiotemporal genomic and cellular analysis was used to isolate and perform RNA sequencing on leader and follower populations from the H1299 NSCLC cell line, and it revealed a leader-specific mutation cluster on chromosome 16q. Genomic data from patients with lung squamous cell carcinoma (LUSC; n = 475) and lung adenocarcinoma (LUAD; n = 501) from The Cancer Genome Atlas were stratified by 16q mutation cluster (16qMC) status (16qMC+ vs 16qMC-) and compared for overall survival (OS), progression-free survival (PFS), and gene set enrichment analysis (GSEA). RESULTS: Poorer OS, poorer PFS, or both were found across all stages and among early-stage patients with 16qMC+ tumors within the LUSC and LUAD cohorts. GSEA revealed 16qMC+ tumors to be enriched for the expression of metastasis- and survival-associated gene sets. CONCLUSIONS: This represents the first leader cell mutation signature identified in patients and has the potential to better stratify high-risk NSCLC and ultimately improve patient outcomes.

Genetic heterogeneity within collective invasion packs drives leader and follower cell phenotypes.

Collective invasion, the coordinated movement of cohesive packs of cells, has become recognized as a major mode of metastasis for solid tumors. These packs are phenotypically heterogeneous and include specialized cells that lead the invasive pack and others that follow behind. To better understand how these unique cell types cooperate to facilitate collective invasion, we analyzed transcriptomic sequence variation between leader and follower populations isolated from the H1299 non-small cell lung cancer cell line using an image-guided selection technique. We now identify 14 expressed mutations that are selectively enriched in leader or follower cells, suggesting a novel link between genomic and phenotypic heterogeneity within a collectively invading tumor cell population. Functional characterization of two phenotype-specific candidate mutations showed that ARP3 enhances collective invasion by promoting the leader cell phenotype and that wild-type KDM5B suppresses chain-like cooperative behavior. These results demonstrate an important role for distinct genetic variants in establishing leader and follower phenotypes and highlight the necessity of maintaining a capacity for phenotypic plasticity during collective cancer invasion.

Cell of Origin and Genetic Alterations in the Pathogenesis of Multiple Myeloma.

B cell activation and differentiation yields plasma cells with high affinity antibodies to a given antigen in a time-frame that allows for host protection. Although the end product is most commonly humoral immunity, the rapid proliferation and somatic mutation of the B cell receptor also results in oncogenic mutations that cause B cell malignancies including plasma cell neoplasms such as multiple myeloma. Myeloma is the second most common hematological malignancy and results in over 100,000 deaths per year worldwide. The genetic alterations that occur in the germinal center, however, are not sufficient to cause myeloma, but rather impart cell proliferation potential on plasma cells, which are normally non-dividing. This pre-malignant state, referred to as monoclonal gammopathy of undetermined significance or MGUS, provides the opportunity for further genetic and epigenetic alterations eventually resulting in a progressive disease that becomes symptomatic. In this review, we will provide a brief history of clonal gammopathies and detail how some of the key discoveries were interwoven with the study of plasma cells. We will also review the genetic and epigenetic alterations discovered over the past 25 years, how these are instrumental to myeloma pathogenesis, and what these events teach us about myeloma and plasma cell biology. These data will be placed in the context of normal B cell development and differentiation and we will discuss how understanding the biology of plasma cells can lead to more effective therapies targeting multiple myeloma.

Multiple myeloma immunoglobulin lambda translocations portend poor prognosis.

Multiple myeloma is a malignancy of antibody-secreting plasma cells. Most patients benefit from current therapies, however, 20% of patients relapse or die within two years and are deemed high risk. Here we analyze structural variants from 795 newly-diagnosed patients as part of the CoMMpass study. We report translocations involving the immunoglobulin lambda (IgL) locus are present in 10% of patients, and indicative of poor prognosis. This is particularly true for IgL-MYC translocations, which coincide with focal amplifications of enhancers at both loci. Importantly, 78% of IgL-MYC translocations co-occur with hyperdiploid disease, a marker of standard risk, suggesting that IgL-MYC-translocated myeloma is being misclassified. Patients with IgL-translocations fail to benefit from IMiDs, which target IKZF1, a transcription factor that binds the IgL enhancer at some of the highest levels in the myeloma epigenome. These data implicate IgL translocation as a driver of poor prognosis which may be due to IMiD resistance.

Structure-Based Engineering of Irreversible Inhibitors against Histone Lysine Demethylase KDM5A.

The active sites of hundreds of human α-ketoglutarate (αKG) and Fe(II)-dependent dioxygenases are exceedingly well preserved, which challenges the design of selective inhibitors. We identified a noncatalytic cysteine (Cys481 in KDM5A) near the active sites of KDM5 histone H3 lysine 4 demethylases, which is absent in other histone demethylase families, that could be explored for interaction with the cysteine-reactive electrophile acrylamide. We synthesized analogs of a thienopyridine-based inhibitor chemotype, namely, 2-((3-aminophenyl)(2-(piperidin-1-yl)ethoxy)methyl)thieno[3,2- b]pyridine-7-carboxylic acid (N70) and a derivative containing a (dimethylamino)but-2-enamido)phenyl moiety (N71) designed to form a covalent interaction with Cys481. We characterized the inhibitory and binding activities against KDM5A and determined the cocrystal structures of the catalytic domain of KDM5A in complex with N70 and N71. Whereas the noncovalent inhibitor N70 displayed αKG-competitive inhibition that could be reversed after dialysis, inhibition by N71 was dependent on enzyme concentration and persisted even after dialysis, consistent with covalent modification.

Insights into the Action of Inhibitor Enantiomers against Histone Lysine Demethylase 5A.

Isomers of chiral drugs can exhibit marked differences in biological activities. We studied the binding and inhibitory activities of 12 compounds against KDM5A. Among them are two pairs of enantiomers representing two distinct inhibitor chemotypes, namely, ( R)- and ( S)-2-((2-chlorophenyl)(2-(piperidin-1-yl)ethoxy)methyl)-1 H-pyrrolo[3,2- b]pyridine-7-carboxylic acid (compounds N51 and N52) and ( R) - and ( S) -N-(1-(3-isopropyl-1 H-pyrazole-5-carbonyl)pyrrolidin-3-yl)cyclopropanecarboxamide (compounds N54 and N55). In vitro, the S enantiomer of the N51/N52 pair (N52) and the R enantiomer of the N54/N55 pair (N54) exhibited about 4- to 5-fold greater binding affinity. The more potent enzyme inhibition of KDM5A by the R-isoform for the cell-permeable N54/N55 pair translated to differences in growth inhibitory activity. We determined structures of the KDM5A catalytic domain in complex with all 12 inhibitors, which revealed the interactions (or lack thereof) responsible for the differences in binding affinity. These results provide insights to guide improvements in binding potency and avenues for development of cell permeable inhibitors of the KDM5 family.

Orphan CpG islands define a novel class of highly active enhancers.

CpG islands (CGI) are critical genomic regulatory elements that support transcriptional initiation and are associated with the promoters of most human genes. CGI are distinguished from the bulk genome by their high CpG density, lack of DNA methylation, and euchromatic features. While CGI are canonically known as strong promoters, thousands of 'orphan' CGI lie far from any known transcript, leaving their function an open question. We undertook a comprehensive analysis of the epigenetic state of orphan CGI across over 100 cell types. Here we show that most orphan CGI display the chromatin features of active enhancers (H3K4me1, H3K27Ac) in at least one cell type. Relative to classical enhancers, these enhancer CGI (ECGI) are stronger, as gauged by chromatin state and in functional assays, are more broadly expressed, and are more highly conserved. Likewise, ECGI engage in more genomic contacts and are enriched for transcription factor binding relative to classical enhancers. In human cancers, these epigenetic differences between ECGI vs. classical enhancers manifest in distinct alterations in DNA methylation. Thus, ECGI define a class of highly active enhancers, strengthened by the broad transcriptional activity, CpG density, hypomethylation, and chromatin features they share with promoter CGI. In addition to indicating a role for thousands of orphan CGI, these findings suggests that enhancer activity may be an intrinsic function of CGI in general and provides new insights into the evolution of enhancers and their epigenetic regulation during development and tumorigenesis.

MAX is an epigenetic sensor of 5-carboxylcytosine and is altered in multiple myeloma.

The oncogenic transcription factor MYC and its binding partner MAX regulate gene expression by binding to DNA at enhancer-box (E-box) elements 5΄-CACGTG-3΄. In mammalian genomes, the central E-box CpG has the potential to be methylated at the 5-position of cytosine (5mC), or to undergo further oxidation to the 5-hydroxymethyl (5hmC), 5-formyl (5fC), or 5-carboxyl (5caC) forms. We find that MAX exhibits the greatest affinity for a 5caC or unmodified C-containing E-box, and much reduced affinities for the corresponding 5mC, 5hmC or 5fC forms. Crystallization of MAX with a 5caC modified E-box oligonucleotide revealed that MAX Arg36 recognizes 5caC using a 5caC-Arg-Guanine triad, with the next nearest residue to the carboxylate group being Arg60. In an analysis of >800 primary multiple myelomas, MAX alterations occurred at a frequency of ∼3%, more than half of which were single nucleotide substitutions affecting a basic clamp-like interface important for DNA interaction. Among these, arginines 35, 36 and 60 were the most frequently altered. In vitro binding studies showed that whereas mutation of Arg36 (R36W) or Arg35 (R35H/L) completely abolished DNA binding, mutation of Arg60 (R60Q) significantly reduced DNA binding, but retained a preference for the 5caC modified E-box. Interestingly, MAX alterations define a subset of myeloma patients with lower MYC expression and a better overall prognosis. Together these data indicate that MAX can act as a direct epigenetic sensor of E-box cytosine modification states and that local CpG modification and MAX variants converge to modulate the MAX-MYC transcriptional network.

Structural Basis for KDM5A Histone Lysine Demethylase Inhibition by Diverse Compounds.

The KDM5/JARID1 family of Fe(II)- and α-ketoglutarate-dependent demethylases removes methyl groups from methylated lysine 4 of histone H3. Accumulating evidence supports a role for KDM5 family members as oncogenic drivers. We compare the in vitro inhibitory properties and binding affinity of ten diverse compounds with all four family members, and present the crystal structures of the KDM5A-linked Jumonji domain in complex with eight of these inhibitors in the presence of Mn(II). All eight inhibitors structurally examined occupy the binding site of α-ketoglutarate, but differ in their specific binding interactions, including the number of ligands involved in metal coordination. We also observed inhibitor-induced conformational changes in KDM5A, particularly those residues involved in the binding of α-ketoglutarate, the anticipated peptide substrate, and intramolecular interactions. We discuss how particular chemical moieties contribute to inhibitor potency and suggest strategies that might be utilized in the successful design of selective and potent epigenetic inhibitors.

A hot L1 retrotransposon evades somatic repression and initiates human colorectal cancer.

Although human LINE-1 (L1) elements are actively mobilized in many cancers, a role for somatic L1 retrotransposition in tumor initiation has not been conclusively demonstrated. Here, we identify a novel somatic L1 insertion in the APC tumor suppressor gene that provided us with a unique opportunity to determine whether such insertions can actually initiate colorectal cancer (CRC), and if so, how this might occur. Our data support a model whereby a hot L1 source element on Chromosome 17 of the patient's genome evaded somatic repression in normal colon tissues and thereby initiated CRC by mutating the APC gene. This insertion worked together with a point mutation in the second APC allele to initiate tumorigenesis through the classic two-hit CRC pathway. We also show that L1 source profiles vary considerably depending on the ancestry of an individual, and that population-specific hot L1 elements represent a novel form of cancer risk.

Factors affecting the persistence of drug-induced reprogramming of the cancer methylome.

Aberrant DNA methylation is a critical feature of cancer. Epigenetic therapy seeks to reverse these changes to restore normal gene expression. DNA demethylating agents, including 5-aza-2'-deoxycytidine (DAC), are currently used to treat certain leukemias, and can sensitize solid tumors to chemotherapy and immunotherapy. However, it has been difficult to pin the clinical efficacy of these agents to specific demethylation events, and the factors that contribute to the durability of response remain largely unknown. Here we examined the genome-wide kinetics of DAC-induced DNA demethylation and subsequent remethylation after drug withdrawal in breast cancer cells. We find that CpGs differ in both their susceptibility to demethylation and propensity for remethylation after drug removal. DAC-induced demethylation was most apparent at CpGs with higher initial methylation levels and further from CpG islands. Once demethylated, such sites exhibited varied remethylation potentials. The most rapidly remethylating CpGs regained >75% of their starting methylation within a month of drug withdrawal. These sites had higher pretreatment methylation levels, were enriched in gene bodies, marked by H3K36me3, and tended to be methylated in normal breast cells. In contrast, a more resistant class of CpG sites failed to regain even 20% of their initial methylation after 3 months. These sites had lower pretreatment methylation levels, were within or near CpG islands, marked by H3K79me2 or H3K4me2/3, and were overrepresented in sites that become aberrantly hypermethylated in breast cancers. Thus, whereas DAC-induced demethylation affects both endogenous and aberrantly methylated sites, tumor-specific hypermethylation is more slowly regained, even as normal methylation promptly recovers. Taken together, these data suggest that the durability of DAC response is linked to its selective ability to stably reset at least a portion of the cancer methylome.

Gene integrated set profile analysis: a context-based approach for inferring biological endpoints.

The identification of genes with specific patterns of change (e.g. down-regulated and methylated) as phenotype drivers or samples with similar profiles for a given gene set as drivers of clinical outcome, requires the integration of several genomic data types for which an 'integrate by intersection' (IBI) approach is often applied. In this approach, results from separate analyses of each data type are intersected, which has the limitation of a smaller intersection with more data types. We introduce a new method, GISPA (Gene Integrated Set Profile Analysis) for integrated genomic analysis and its variation, SISPA (Sample Integrated Set Profile Analysis) for defining respective genes and samples with the context of similar, a priori specified molecular profiles. With GISPA, the user defines a molecular profile that is compared among several classes and obtains ranked gene sets that satisfy the profile as drivers of each class. With SISPA, the user defines a gene set that satisfies a profile and obtains sample groups of profile activity. Our results from applying GISPA to human multiple myeloma (MM) cell lines contained genes of known profiles and importance, along with several novel targets, and their further SISPA application to MM coMMpass trial data showed clinical relevance.