Targeting neddylation and sumoylation in chemoresistant triple negative breast cancer.

Triple negative breast cancer (TNBC) accounts for 15-20% of breast cancer cases in the United States. Systemic neoadjuvant chemotherapy (NACT), with or without immunotherapy, is the current standard of care for patients with early-stage TNBC. However, up to 70% of TNBC patients have significant residual disease once NACT is completed, which is associated with a high risk of developing recurrence within two to three years of surgical resection. To identify targetable vulnerabilities in chemoresistant TNBC, we generated longitudinal patient-derived xenograft (PDX) models from TNBC tumors before and after patients received NACT. We then compiled transcriptomes and drug response profiles for all models. Transcriptomic analysis identified the enrichment of aberrant protein homeostasis pathways in models from post-NACT tumors relative to pre-NACT tumors. This observation correlated with increased sensitivity in vitro to inhibitors targeting the proteasome, heat shock proteins, and neddylation pathways. Pevonedistat, a drug annotated as a NEDD8-activating enzyme (NAE) inhibitor, was prioritized for validation in vivo and demonstrated efficacy as a single agent in multiple PDX models of TNBC. Pharmacotranscriptomic analysis identified a pathway-level correlation between pevonedistat activity and post-translational modification (PTM) machinery, particularly involving neddylation and sumoylation targets. Elevated levels of both NEDD8 and SUMO1 were observed in models exhibiting a favorable response to pevonedistat compared to those with a less favorable response in vivo. Moreover, a correlation emerged between the expression of neddylation-regulated pathways and tumor response to pevonedistat, indicating that targeting these PTM pathways may prove effective in combating chemoresistant TNBC.

PTEN in triple-negative breast carcinoma: protein expression and genomic alteration in pretreatment and posttreatment specimens.

Background: Recent advances have been made in targeting the phosphoinositide 3-kinase pathway in breast cancer. Phosphatase and tensin homolog (PTEN) is a key component of that pathway. Objective: To understand the changes in PTEN expression over the course of the disease in patients with triple-negative breast cancer (TNBC) and whether PTEN copy number variation (CNV) by next-generation sequencing (NGS) can serve as an alternative to immunohistochemistry (IHC) to identify PTEN loss. Methods: We compared PTEN expression by IHC between pretreatment tumors and residual tumors in the breast and lymph nodes after neoadjuvant chemotherapy in 96 patients enrolled in a TNBC clinical trial. A correlative analysis between PTEN protein expression and PTEN CNV by NGS was also performed. Results: With a stringent cutoff for PTEN IHC scoring, PTEN expression was discordant between pretreatment and posttreatment primary tumors in 5% of patients (n = 96) and between posttreatment primary tumors and lymph node metastases in 9% (n = 33). A less stringent cutoff yielded similar discordance rates. Intratumoral heterogeneity for PTEN loss was observed in 7% of the patients. Among pretreatment tumors, PTEN copy numbers by whole exome sequencing (n = 72) were significantly higher in the PTEN-positive tumors by IHC compared with the IHC PTEN-loss tumors (p < 0.0001). However, PTEN-positive and PTEN-loss tumors by IHC overlapped in copy numbers: 14 of 60 PTEN-positive samples showed decreased copy numbers in the range of those of the PTEN-loss tumors. Conclusion: Testing various specimens by IHC may generate different PTEN results in a small proportion of patients with TNBC; therefore, the decision of testing one versus multiple specimens in a clinical trial should be defined in the patient inclusion criteria. Although a distinct cutoff by which CNV differentiated PTEN-positive tumors from those with PTEN loss was not identified, higher copy number of PTEN may confer positive PTEN, whereas lower copy number of PTEN would necessitate additional testing by IHC to assess PTEN loss. Trial registration: NCT02276443.

Phylogenetic inference from single-cell RNA-seq data.

Tumors are comprised of subpopulations of cancer cells that harbor distinct genetic profiles and phenotypes that evolve over time and during treatment. By reconstructing the course of cancer evolution, we can understand the acquisition of the malignant properties that drive tumor progression. Unfortunately, recovering the evolutionary relationships of individual cancer cells linked to their phenotypes remains a difficult challenge. To address this need, we have developed PhylinSic, a method that reconstructs the phylogenetic relationships among cells linked to their gene expression profiles from single cell RNA-sequencing (scRNA-Seq) data. This method calls nucleotide bases using a probabilistic smoothing approach and then estimates a phylogenetic tree using a Bayesian modeling algorithm. We showed that PhylinSic identified evolutionary relationships underpinning drug selection and metastasis and was sensitive enough to identify subclones from genetic drift. We found that breast cancer tumors resistant to chemotherapies harbored multiple genetic lineages that independently acquired high K-Ras and ß-catenin, suggesting that therapeutic strategies may need to control multiple lineages to be durable. These results demonstrated that PhylinSic can reconstruct evolution and link the genotypes and phenotypes of cells across monophyletic tumors using scRNA-Seq.

Simple combination of multiple somatic variant callers to increase accuracy.

Publications comparing variant caller algorithms present discordant results with contradictory rankings. Caller performances are inconsistent and wide ranging, and dependent upon input data, application, parameter settings, and evaluation metric. With no single variant caller emerging as a superior standard, combinations or ensembles of variant callers have appeared in the literature. In this study, a whole genome somatic reference standard was used to derive principles to guide strategies for combining variant calls. Then, manually annotated variants called from the whole exome sequencing of a tumor were used to corroborate these general principles. Finally, we examined the ability of these principles to reduce noise in targeted sequencing.

Multiomic analysis of homologous recombination-deficient end-stage high-grade serous ovarian cancer.

High-grade serous ovarian cancer (HGSC) is frequently characterized by homologous recombination (HR) DNA repair deficiency and, while most such tumors are sensitive to initial treatment, acquired resistance is common. We undertook a multiomics approach to interrogate molecular diversity in end-stage disease, using multiple autopsy samples collected from 15 women with HR-deficient HGSC. Patients had polyclonal disease, and several resistance mechanisms were identified within most patients, including reversion mutations and HR restoration by other means. We also observed frequent whole-genome duplication and global changes in immune composition with evidence of immune escape. This analysis highlights diverse evolutionary changes within HGSC that evade therapy and ultimately overwhelm individual patients.

Predictors of success in establishing orthotopic patient-derived xenograft models of triple negative breast cancer.

Patient-derived xenograft (PDX) models of breast cancer are an effective discovery platform and tool for preclinical pharmacologic testing and biomarker identification. We established orthotopic PDX models of triple negative breast cancer (TNBC) from the primary breast tumors of patients prior to and following neoadjuvant chemotherapy (NACT) while they were enrolled in the ARTEMIS trial (NCT02276443). Serial biopsies were obtained from patients prior to treatment (pre-NACT), from poorly responsive disease after four cycles of Adriamycin and cyclophosphamide (AC, mid-NACT), and in cases of AC-resistance, after a 3-month course of different experimental therapies and/or additional chemotherapy (post-NACT). Our study cohort includes a total of 269 fine needle aspirates (FNAs) from 217 women, generating a total of 62 PDX models (overall success-rate = 23%). Success of PDX engraftment was generally higher from those cancers that proved to be treatment-resistant, whether poorly responsive to AC as determined by ultrasound measurements mid-NACT (p = 0.063), RCB II/III status after NACT (p = 0.046), or metastatic relapse within 2 years of surgery (p = 0.008). TNBC molecular subtype determined from gene expression microarrays of pre-NACT tumors revealed no significant association with PDX engraftment rate (p = 0.877). Finally, we developed a statistical model predictive of PDX engraftment using percent Ki67 positive cells in the patient's diagnostic biopsy, positive lymph node status at diagnosis, and low volumetric reduction of the patient's tumor following AC treatment. This novel bank of 62 PDX models of TNBC provides a valuable resource for biomarker discovery and preclinical therapeutic trials aimed at improving neoadjuvant response rates for patients with TNBC.

Glycolysis regulates KRAS plasma membrane localization and function through defined glycosphingolipids.

Oncogenic KRAS expression generates a metabolic dependency on aerobic glycolysis, known as the Warburg effect. We report an effect of increased glycolytic flux that feeds into glycosphingolipid biosynthesis and is directly linked to KRAS oncogenic function. High resolution imaging and genetic approaches show that a defined subset of outer leaflet glycosphingolipids, including GM3 and SM4, is required to maintain KRAS plasma membrane localization, with GM3 engaging in cross-bilayer coupling to maintain inner leaflet phosphatidylserine content. Thus, glycolysis is critical for KRAS plasma membrane localization and nanoscale spatial organization. Reciprocally oncogenic KRAS selectively upregulates cellular content of these same glycosphingolipids, whose depletion in turn abrogates KRAS oncogenesis in pancreatic cancer models. Our findings expand the role of the Warburg effect beyond ATP generation and biomass building to high-level regulation of KRAS function. The positive feedforward loop between oncogenic KRAS signaling and glycosphingolipid synthesis represents a vulnerability with therapeutic potential.

Predicting clinical response to everolimus in ER+ breast cancers using machine-learning.

Endocrine therapy remains the primary treatment choice for ER+ breast cancers. However, most advanced ER+ breast cancers ultimately develop resistance to endocrine. This acquired resistance to endocrine therapy is often driven by the activation of the PI3K/AKT/mTOR signaling pathway. Everolimus, a drug that targets and inhibits the mTOR complex has been shown to improve clinical outcomes in metastatic ER+ breast cancers. However, there are no biomarkers currently available to guide the use of everolimus in the clinic for progressive patients, where multiple therapeutic options are available. Here, we utilized gene expression signatures from 9 ER+ breast cancer cell lines and 23 patients treated with everolimus to develop and validate an integrative machine learning biomarker of mTOR inhibitor response. Our results show that the machine learning biomarker can successfully distinguish responders from non-responders and can be applied to identify patients that will most likely benefit from everolimus treatment.

Human NOP2/NSUN1 regulates ribosome biogenesis through non-catalytic complex formation with box C/D snoRNPs.

5-Methylcytosine (m5C) is a base modification broadly found on various RNAs in the human transcriptome. In eukaryotes, m5C is catalyzed by enzymes of the NSUN family composed of seven human members (NSUN1-7). NOP2/NSUN1 has been primarily characterized in budding yeast as an essential ribosome biogenesis factor required for the deposition of m5C on the 25S ribosomal RNA (rRNA). Although human NOP2/NSUN1 has been known to be an oncogene overexpressed in several types of cancer, its functions and substrates remain poorly characterized. Here, we used a miCLIP-seq approach to identify human NOP2/NSUN1 RNA substrates. Our analysis revealed that NOP2/NSUN1 catalyzes the deposition of m5C at position 4447 on the 28S rRNA. We also find that NOP2/NSUN1 binds to the 5'ETS region of the pre-rRNA transcript and regulates pre-rRNA processing through non-catalytic complex formation with box C/D snoRNAs. We provide evidence that NOP2/NSUN1 facilitates the recruitment of U3 and U8 snoRNAs to pre-90S ribosomal particles and their stable assembly into snoRNP complexes. Remarkably, expression of both WT and catalytically inactive NOP2/NSUN1 in knockdown background rescues the rRNA processing defects and the stable assembly of box C/D snoRNP complexes, suggesting that NOP2/NSUN1-mediated deposition of m5C on rRNA is not required for ribosome synthesis.

Tumor Immune Microenvironment Changes by Multiplex Immunofluorescence Staining in a Pilot Study of Neoadjuvant Talazoparib for Early-Stage Breast Cancer Patients with a Hereditary BRCA Mutation.

PURPOSE: The immunological profile of early-stage breast cancer treated with neoadjuvant PARP inhibitors has not been described. The aim of this study was to delineate the changes in the tumor immune microenvironment (TiME) induced by talazoparib. PATIENTS AND METHODS: Patients with operable germline BRCA1/2 pathogenic variant (gBRCA1/2+) breast cancer were enrolled in a feasibility study of neoadjuvant talazoparib. Thirteen patients who received 8 weeks of neoadjuvant talazoparib were available for analysis, including 11 paired pre- and post-talazoparib core biopsies. Treatment-related changes in tumor-infiltrating lymphocytes were examined and immune cell phenotypes and their spatial distribution in the TiME were identified and quantified by multiplex immunofluorescence using a panel of 6 biomarkers (CD3, CD8, CD68, PD-1, PD-L1, and CK). RESULTS: Neoadjuvant talazoparib significantly increased infiltrating intratumoral and stromal T-cell and cytotoxic T-cell density. There was no difference in PD-1 or PD-L1 immune cell phenotypes in the pre- and post-talazoparib specimens and PD-L1 expression in tumor cells was rare in this cohort. Spatial analysis demonstrated that pre-talazoparib interactions between macrophages and T cells may correlate with pathologic complete response. CONCLUSIONS: This is the first study with phenotyping to characterize the immune response to neoadjuvant talazoparib in patients with gBRCA1/2+ breast cancer. These findings support an emerging role for PARP inhibitors in enhancing tumor immunogenicity. Further investigation of combinatorial strategies is warranted with agents that exploit the immunomodulatory effects of PARP inhibitors on the TiME.

Identification of biomarkers of response to preoperative talazoparib monotherapy in treatment naïve gBRCA+ breast cancers.

Germline mutations in BRCA1 or BRCA2 exist in ~2-7% of breast cancer patients, which has led to the approval of PARP inhibitors in the advanced setting. We have previously reported a phase II neoadjuvant trial of single agent talazoparib for patients with germline BRCA pathogenic variants with a pathologic complete response (pCR) rate of 53%. As nearly half of the patients treated did not have pCR, better strategies are needed to overcome treatment resistance. To this end, we conducted multi-omic analysis of 13 treatment naïve breast cancer tumors from patients that went on to receive single-agent neoadjuvant talazoparib. We looked for biomarkers that were predictive of response (assessed by residual cancer burden) after 6 months of therapy. We found that all resistant tumors exhibited either the loss of SHLD2, expression of a hypoxia signature, or expression of a stem cell signature. These results indicate that the deep analysis of pre-treatment tumors can identify biomarkers that are predictive of response to talazoparib and potentially other PARP inhibitors, and provides a framework that will allow for better selection of patients for treatment, as well as a roadmap for the development of novel combination therapies to prevent emergence of resistance.

Molecular Characterization and Prospective Evaluation of Pathologic Response and Outcomes with Neoadjuvant Therapy in Metaplastic Triple-Negative Breast Cancer.

PURPOSE: Metaplastic breast cancer (MpBC) is a rare subtype of breast cancer that is commonly triple-negative and poorly responsive to neoadjuvant therapy in retrospective studies. EXPERIMENTAL DESIGN: To better define clinical outcomes and correlates of response, we analyzed the rate of pathologic complete response (pCR) to neoadjuvant therapy, survival outcomes, and genomic and transcriptomic profiles of the pretreatment tumors in a prospective clinical trial (NCT02276443). A total of 211 patients with triple-negative breast cancer (TNBC), including 39 with MpBC, received doxorubicin-cyclophosphamide-based neoadjuvant therapy. RESULTS: Although not meeting the threshold for statistical significance, patients with MpBCs were less likely to experience a pCR (23% vs. 40%; P = 0.07), had shorter event-free survival (29.4 vs. 32.2 months, P = 0.15), metastasis-free survival (30.3 vs. 32.4 months, P = 0.22); and overall survival (32.6 vs. 34.3 months, P = 0.21). This heterogeneity is mirrored in the molecular profiling. Mutations in PI3KCA (23% vs. 9%, P = 0.07) and its pathway (41% vs. 18%, P = 0.02) were frequently observed and enriched in MpBCs. The gene expression profiles of each histologically defined subtype were distinguishable and characterized by distinctive gene signatures. Among nonmetaplastic (non-Mp) TNBCs, 10% possessed a metaplastic-like gene expression signature and had pCR rates and survival outcomes similar to MpBC. CONCLUSIONS: Further investigations will determine if metaplastic-like tumors should be treated more similarly to MpBC in the clinic. The 23% pCR rate in this study suggests that patients with MpBC should be considered for NAT. To improve this rate, a pathway analysis predicted enrichment of histone deacetylase (HDAC) and RTK/MAPK pathways in MpBC, which may serve as new targetable vulnerabilities.

Chemotherapy Coupled to Macrophage Inhibition Induces T-cell and B-cell Infiltration and Durable Regression in Triple-Negative Breast Cancer.

Immunosuppressive elements within the tumor microenvironment, such as tumor-associated macrophages (TAM), can present a barrier to successful antitumor responses by cytolytic T cells. Here we employed preclinical syngeneic p53 null mouse models of triple-negative breast cancer (TNBC) to develop a treatment regimen that harnessed the immunostimulatory effects of low-dose cyclophosphamide coupled with the pharmacologic inhibition of TAMs using either a small-molecule CSF1R inhibitor or an anti-CSF1R antibody. This therapeutic combination was effective in treating several highly aggressive TNBC murine mammary tumor and lung metastasis models. Single-cell RNA sequencing characterized tumor-infiltrating lymphocytes including Th cells and antigen-presenting B cells that were highly enriched in responders to combination therapy. In one model that exhibited long-term posttreatment tumor regression, high-dimensional imaging techniques identified the close spatial localization of B220+/CD86+-activated B cells and CD4+ T cells in tertiary lymphoid structures that were present up to 6 weeks posttreatment. The transcriptional and metabolic heterogeneity of TAMs was also characterized in two closely related claudin-low/mesenchymal subtype tumor models with differential treatment responses. A murine TAM signature derived from the T12 model was highly conserved in human claudin-low breast cancers, and high expression of the TAM signature correlated with reduced overall survival in patients with breast cancer. This TAM signature may help identify human patients with claudin-low breast cancer that will benefit from the combination of cyclophosphamide and anti-CSF1R therapy. These studies illustrate the complexity of the tumor immune microenvironment and highlight different immune responses that result from rational immunotherapy combinations. SIGNIFICANCE: Immunostimulatory chemotherapy combined with pharmacologic inhibition of TAMs results in durable treatment responses elicited by Th cells and B cells in claudin-low TNBC models.

Prognostic Impact of High Baseline Stromal Tumor-Infiltrating Lymphocytes in the Absence of Pathologic Complete Response in Early-Stage Triple-Negative Breast Cancer.

High stromal tumor-infiltrating lymphocytes (sTILs) are associated with an improved pathologic complete response (pCR) and survival in triple-negative breast cancer (TNBC). We hypothesized that high baseline sTILs would have a favorable prognostic impact in TNBC patients without a pCR after neoadjuvant chemotherapy (NACT). In this prospective NACT study, pretreatment biopsies from 318 patients with early-stage TNBC were evaluated for sTILs. Recursive partitioning analysis (RPA) was applied to search for the sTIL cutoff best associated with a pCR. With ≥20% sTILs identified as the optimal cutoff, 33% patients had high sTILs (pCR rate 64%) and 67% had low sTILs (pCR rate 29%). Patients were stratified according to the sTIL cutoff (low vs. high) and response to NACT (pCR vs. residual disease (RD)). The primary endpoint was event-free survival (EFS), with hazard ratios calculated using the Cox proportional hazards regression model and the 3-year restricted mean survival time (RMST) as primary measures. Within the high-sTIL group, EFS was better in patients with a pCR compared with those with RD (HR 0.05; 95% CI 0.01-0.39; p = 0.004). The difference in the 3-year RMST for EFS between the two groups was 5.6 months (95% CI 2.3-8.8; p = 0.001). However, among patients with RD, EFS was not significantly different between those with high sTILs and those with low sTILs (p = 0.7). RNA-seq analysis predicted more CD8+ T cells in the high-sTIL group with favorable EFS compared with the high-sTIL group with unfavorable EFS. This study did not demonstrate that high baseline sTILs confer a benefit in EFS in the absence of a pCR.

ABCA1 Expression Is Upregulated in an EMT in Breast Cancer Cell Lines via MYC-Mediated De-Repression of Its Proximal Ebox Element.

The ATP-Binding Cassette transporter A1 (ABCA1) reverse cholesterol transport channel has been associated with a number of phenotypes in breast cancer, including reduced proliferation and increased metastatic capacity. It is induced in an epithelial-mesenchymal transition (EMT), but little is known about how this occurs, and whether it is sufficient to promote metastatic phenotypes. To address these questions, we have deciphered the transcriptional regulation of ABCA1 across EMT states and found that it is repressed by MYC via an E-box element in its P1 alternative promoter. De-repression of the promoter by MYC knockdown leads to induction of ABCA1 expression. This indicates that ABCA1 expression is regulated in an EMT, revealing another link between ABCA1 and malignant phenotypes.

CKB inhibits epithelial-mesenchymal transition and prostate cancer progression by sequestering and inhibiting AKT activation.

Epithelial-mesenchymal transition (EMT) contributes to tumor invasion, metastasis and drug resistance. AKT activation is key in a number of cellular processes. While many positive regulators for either EMT or AKT activation have been reported, few negative regulators are established. Through kinase cDNA screen, we identified brain-type creatine kinase (CKB or BCK) as a potent suppressor for both. As a ubiquitously expressed kinase in normal tissues, CKB is significantly downregulated in several solid cancer types. Lower CKB expression is significantly associated with worse prognosis. Phenotypically, CKB overexpression suppresses, while its silencing promotes, EMT and cell migration, xenograft tumor growth and metastasis of prostate cancer cells. AKT activation is one of the most prominent signaling events upon CKB silencing in prostate cancer cells, which is in line with prostate cancer TCGA data. EMT enhanced by CKB silencing is abolished by AKT inhibition. Mechanistically, CKB interacts with AKT and sequestrates it from activation by mTOR. We further elucidated that an 84aa fragment at C-terminus of CKB protein interacts with AKT's PH domain. Ectopic expression of the 84aa CKB fragment inhibits AKT activation, EMT and cell proliferation. Interestingly, molecular dynamics simulation on crystal structures of AKT and CKB independently demonstrates that AKT's PH domain and CKB's 84aa fragment establish their major interaction interface. In summary, we have discovered CKB as a negative regulator of EMT and AKT activation, revealing a new mode of their regulation . We have also demonstrated that CKB downregulation is a poor prognosticator, which is sufficient to promote prostate cancer progression.

Serial single-cell genomics reveals convergent subclonal evolution of resistance as early-stage breast cancer patients progress on endocrine plus CDK4/6 therapy.

Combining cyclin-dependent kinase (CDK) inhibitors with endocrine therapy improves outcomes for metastatic estrogen receptor positive (ER+) breast cancer patients but its value in earlier stage patients is unclear. We examined evolutionary trajectories of early-stage breast cancer tumors, using single cell RNA sequencing (scRNAseq) of serial biopsies from the FELINE clinical trial (#NCT02712723) of endocrine therapy (letrozole) alone or combined with the CDK inhibitor ribociclib. Despite differences in subclonal diversity evolution across patients and treatments, common resistance phenotypes emerged. Resistant tumors treated with combination therapy showed accelerated loss of estrogen signaling with convergent up-regulation of JNK signaling through growth factor receptors. In contrast, cancer cells maintaining estrogen signaling during mono- or combination therapy showed potentiation of CDK4/6 activation and ERK upregulation through ERBB4 signaling. These results indicate that combination therapy in early-stage ER+ breast cancer leads to emergence of resistance through a shift from estrogen to alternative growth signal-mediated proliferation.

Mammary-specific expression of Trim24 establishes a mouse model of human metaplastic breast cancer.

Conditional overexpression of histone reader Tripartite motif containing protein 24 (TRIM24) in mouse mammary epithelia (Trim24COE) drives spontaneous development of mammary carcinosarcoma tumors, lacking ER, PR and HER2. Human carcinosarcomas or metaplastic breast cancers (MpBC) are a rare, chemorefractory subclass of triple-negative breast cancers (TNBC). Comparison of Trim24COE metaplastic carcinosarcoma morphology, TRIM24 protein levels and a derived Trim24COE gene signature reveals strong correlation with human MpBC tumors and MpBC patient-derived xenograft (PDX) models. Global and single-cell tumor profiling reveal Met as a direct oncogenic target of TRIM24, leading to aberrant PI3K/mTOR activation. Here, we find that pharmacological inhibition of these pathways in primary Trim24COE tumor cells and TRIM24-PROTAC treatment of MpBC TNBC PDX tumorspheres decreased cellular viability, suggesting potential in therapeutically targeting TRIM24 and its regulated pathways in TRIM24-expressing TNBC.

A method of sample-wise region-set enrichment analysis for DNA methylomics.

Aim: Gene set analysis has commonly been used to interpret DNA methylome data. However, summarizing the DNA methylation level of a gene is challenging due to variability in the number, density and methylation levels of CpG sites, and the numerous intergenic CpGs. Instead, we propose to use region sets to annotate the DNA methylome. Methods: We developed single sample region-set enrichment analysis for DNA methylome (methyl-ssRSEA) to conduct sample-wise, region-set enrichment analysis. Results: Methyl-ssRSEA can handle both microarray- and sequencing-based platforms and reproducibly recover the known biology from the methylation profiles of peripheral blood cells and breast cancers. The performance was superior to existing tools for region-set analysis in discriminating blood cell types. Conclusion: Methyl-ssRSEA offers a novel way to functionally interpret the DNA methylome in the cell.

Lay abstract Gene set analysis has been a common way to understand the meaning of DNA methylome data. However, organizing the DNA methylation level of a gene is challenging due to variation in the number, density and extent of methylation, of methylation sites, and the substantial number of methylation sites between genes. Instead, we propose to use region sets for the organization. We developed single sample region-set enrichment analysis for DNA methylome (methyl-ssRSEA) to conduct region-set analysis for every sample. Methyl-ssRSEA can handle both microarray- and sequencing-based methods and repeatedly find the known characters from the methylation patterns of peripheral blood cells and breast cancers. The performance was better than existing tools for region-set analysis in differentiating blood cell types. Methyl-ssRSEA offers a novel way to find the features of DNA methylome in the cell.