Inference of genomic lesions from single-cell RNAseq in myeloma improves functional intra- and inter-clonal analysis.

Smoldering Multiple Myeloma (SMM) is an asymptomatic plasma cell (PC) neoplasm that may evolve with variable frequency into multiple myeloma (MM). SMM is initiated by chromosomal translocations involving the IgH locus or by hyperdiploidy and evolves through acquisition of additional genetic lesions. In this scenario, we aimed at establishing a reliable analysis pipeline to infer genomic lesions from transcriptomic analysis, by combining single-cell RNA sequencing (scRNA-seq) with B-cell receptor sequencing and copy- number abnormality (CNA) analysis to identify clonal PCs at the genetic level along their specific transcriptional landscape. We profiled 20,465 bone marrow (BM) PCs derived from five SMM/MM patients and unbiasedly identified clonal and polyclonal plasma cells. Hyperdiploidy, t(11;14) and t(6;14) were identified at the scRNA level by analysis of chimeric reads. Subclone functional analysis was improved by combining transcriptome with CNA analysis. As examples, we illustrate the different functional properties of a light chain escape subclone in SMM, and of different B-cell and PC subclones in a patient affected by Wäldenstrom Macroglobulinemia and SMM. Overall, our data provide a proof of principle for inference of clinically relevant genotypic data from scRNAseq, which in turn will refine functional annotation of the clonal architecture of PC dyscrasias.

Tumorigenic role of Musashi-2 in aggressive mantle cell lymphoma.

SOX11 overexpression has been associated with aggressive behavior of mantle cell lymphomas (MCL). SOX11 is overexpressed in embryonic and cancer stem cells (CSC) of some tumors. Although CSC have been isolated from primary MCL, their relationship to SOX11 expression and contribution to MCL pathogenesis and clinical evolution remain unknown. Here, we observed enrichment in leukemic and hematopoietic stem cells gene signatures in SOX11+ compared to SOX11- MCL primary cases. Musashi-2 (MSI2) emerged as one of the most significant upregulated stem cell-related genes in SOX11+ MCLs. SOX11 is directly bound to the MSI2 promoter upregulating its expression in vitro. MSI2 intronic enhancers were strongly activated in SOX11+ MCL cell lines and primary cases. MSI2 upregulation was significantly associated with poor overall survival independently of other high-risk features of MCL. MSI2 knockdown decreased the expression of genes related to apoptosis and stem cell features and significantly reduced clonogenic growth, tumor cell survival and chemoresistance in MCL cells. MSI2-knockdown cells had reduced tumorigenic engraftment into mice bone marrow and spleen compared to control cells in xenotransplanted mouse models. Our results suggest that MSI2 might play a key role in sustaining stemness and tumor cell survival, representing a possible novel target for therapeutic interventions in MCL.

Remission of obesity and insulin resistance is not sufficient to restore mitochondrial homeostasis in visceral adipose tissue.

Metabolic plasticity is the ability of a biological system to adapt its metabolic phenotype to different environmental stressors. We used a whole-body and tissue-specific phenotypic, functional, proteomic, metabolomic and transcriptomic approach to systematically assess metabolic plasticity in diet-induced obese mice after a combined nutritional and exercise intervention. Although most obesity and overnutrition-related pathological features were successfully reverted, we observed a high degree of metabolic dysfunction in visceral white adipose tissue, characterized by abnormal mitochondrial morphology and functionality. Despite two sequential therapeutic interventions and an apparent global healthy phenotype, obesity triggered a cascade of events in visceral adipose tissue progressing from mitochondrial metabolic and proteostatic alterations to widespread cellular stress, which compromises its biosynthetic and recycling capacity. In humans, weight loss after bariatric surgery showed a transcriptional signature in visceral adipose tissue similar to our mouse model of obesity reversion. Overall, our data indicate that obesity prompts a lasting metabolic fingerprint that leads to a progressive breakdown of metabolic plasticity in visceral adipose tissue.

Transgene codon usage drives viral fitness and therapeutic efficacy in oncolytic adenoviruses.

Arming oncolytic adenoviruses with therapeutic transgenes is a well-established strategy for multimodal tumour attack. However, this strategy sometimes leads to unexpected attenuated viral replication and a loss of oncolytic effects, preventing these viruses from reaching the clinic. Previous work has shown that altering codon usage in viral genes can hamper viral fitness. Here, we have analysed how transgene codon usage impacts viral replication and oncolytic activity. We observe that, although transgenes with optimized codons show high expression levels at the first round of infection, they impair viral fitness and are therefore not expressed in a sustained manner. Conversely, transgenes encoded by suboptimal codons do not compromise viral replication and are thus stably expressed over time, allowing a greater oncolytic activity both in vitro and in vivo. Altogether, our work shows that fine-tuning codon usage leads to a concerted optimization of transgene expression and viral replication paving the way for the rational design of more efficacious oncolytic therapies.

Investigating Neuroimaging Correlates of Early Frailty in Patients With Behavioral Variant Frontotemporal Dementia: A MRI and FDG-PET Study.

Frailty is a dynamic clinical condition characterized by the reduction of interconnections among different psychobiological domains, which leads to a homeostatic vulnerability. The association between physical frailty and cognitive dysfunctions is a possible predictor of poor prognosis in patients with neurodegenerative disorders. However, this construct has not been fully analyzed by a multidimensional neuropsychogeriatric assessment matched with multimodal neuroimaging methods in patients with behavioral variant frontotemporal dementia (bvFTD). We have investigated cognitive dysfunctions and frailty status, assessed by both a neuropsychological evaluation and the Multidimensional Prognostic Index (MPI), in a sample of 18 bvFTD patients and compared to matched healthy controls. Gray matter (GM) volume (as assessed by voxel-based morphometry) and metabolism (on 18fluorodeoxyglucose positron emission tomography) were first separately compared between groups, then voxelwise compared and correlated to each other within patients. Linear regression of the MPI was performed on those voxels presenting a significant correlation between altered GM volume and metabolism. The neuropsychological assessment reflected the diagnoses and the functional-anatomical alterations documented by neuroimaging analyses. In particular, the majority of patients presented significant executive dysfunction and mood changes in terms of apathy, depression, and anxiety. In the overall MPI score, the patients fell in the lower range (indicating an early frailty status). On imaging, they exhibited a bilateral decrease of GM density and hypometabolism involving the frontal pole, the anterior opercular region, and the anterior cingulate cortex. Greater atrophy than hypometabolism was observed in the bilateral orbitofrontal cortex, the triangular part of the inferior frontal gyrus, and the ventral striatum, whereas the contrary was detected in the bilateral dorsal anterior cingulate cortex and pre-supplementary motor area. MPI scores significantly correlated only with the co-occurrence of a decrease of GM density and hypometabolism in the right anterior insular cortex, but not with the separate pathological phenomena. Our results show a correlation between a specific pattern of co-occurring GM atrophy and hypometabolism with early frailty in bvFTD patients. These aspects, combined with executive dysfunction and mood changes, may lead to an increased risk of poor prognosis, highlighting a potentially critical and precocious role of the insula in the pathogenesis of frailty.

Predictive factors of a worse response to radioactive Iodine-I131 treatment in hyperthyroidism: outcome analysis in 424 patients. A single centre experience.

PURPOSE: Aim of our study was to search for variables associated with worse outcomes in patients treated with radioactive iodine (RAI) for hyperthyroidism by a dosimetric-based approach. METHODS: Four hundred twenty-four patients with hyperthyroidism related to Toxic Multinodular Goiter (TMG; n = 213), Grave's disease (GD; n = 150) and toxic adenoma (TA; n = 61) treated with RAI between 2000 and 2018 and with at least 12 months follow-up were retrospectively evaluated. Association between outcomes (response vs. no response) at 6 and 12 months and baseline TSH values, anti-thyroid drugs (ATD) duration and posology, RAI absorbed dose and dimensional reduction of target mass at ultrasound was evaluated by Mann-Whitney test. Risk factors for response vs. no-response were analysed by binary logistic regression model. RESULTS: Overall response rate was 78.7 and 83% at 6 and 12 months, respectively. Both at 6 and 12 months higher TSH baseline values (p < 0.001), lower ATD duration (p = 0.004 and p = 0.043), lower ATD posology (p = 0.014 and p = 0.005), and lower dose to target (DT) (327 vs. 373 Gy, p = 0.003) were associated to response. Longer ATD duration and higher ATD posology were independent risk factors for no response at 6 and 12 months in GD and TMG, with no response at 6 months in TA subgroups. CONCLUSIONS: Low TSH levels, longer duration and higher posology of ATD were associated with worse response to RAI. These data confirm that RAI therapy should be considered earlier in patients' management to allow better outcome and avoid ATD toxicity.

A Cyclin D1-Dependent Transcriptional Program Predicts Clinical Outcome in Mantle Cell Lymphoma.

PURPOSE: Mantle cell lymphoma (MCL) is characterized by the t(11;14)(q13;q32) translocation leading to cyclin D1 overexpression. Cyclin D1 is a major cell-cycle regulator and also regulates transcription, but the impact of cyclin D1-mediated transcriptional dysregulation on MCL pathogenesis remains poorly understood. The aim of this study was to define a cyclin D1-dependent gene expression program and analyze its prognostic value. EXPERIMENTAL DESIGN: We integrated genome-wide expression analysis of cyclin D1-silenced and overexpressing cells with cyclin D1 chromatin-binding profiles to identify a cyclin D1-dependent transcriptional program in MCL cells. We analyzed this gene program in two MCL series of peripheral blood samples (n = 53) and lymphoid tissues (n = 106) to determine its biological and clinical relevance. We then obtained a simplified signature of this program and evaluated a third series of peripheral blood MCL samples (n = 81) by NanoString gene expression profiling to validate our findings. RESULTS: We identified a cyclin D1-dependent transcriptional program composed of 295 genes that were mainly involved in cell-cycle control. The cyclin D1-dependent gene program was overexpressed in MCL tumors directly proportional to cyclin D1 levels. High expression of this program conferred an adverse prognosis with significant shorter overall survival of the patients. These observations were validated in an independent cohort of patients using a simplified 37-gene cyclin D1 signature. The cyclin D1-dependent transcriptional program was also present in multiple myeloma and breast tumors with cyclin D1 overexpression. CONCLUSIONS: We identified a cyclin D1-dependent transcriptional program that is overexpressed in MCL and predicts clinical outcome.

Cell Plasticity-Related Phenotypes and Taxanes Resistance in Castration-Resistant Prostate Cancer.

The prostatic tumor cells plasticity is involved in resistance to hormone-therapy, allowing these cells to survive despite androgen receptor inhibition. However, its role in taxanes resistance has not been fully established. Gene expression of plasticity-related phenotypes such as epithelial-mesenchymal transition (EMT), stem cell-like and neuroendocrine (NE) phenotypes was studied in vitro, in silico, in circulating tumor cells (CTCs) (N=22) and in tumor samples (N=117) from taxanes-treated metastatic castration-resistant prostate cancer (mCRPC) patients. Docetaxel (D)-resistant cells presented a more pronounced EMT phenotype than cabazitaxel (CZ)-resistant cells. In silico analysis revealed ESRP1 down-regulation in taxane-exposed mCRPC samples. Cell plasticity-related changes occurred in CTCs after taxanes treatment. Tumor EMT phenotype was associated with lower PSA progression-free survival (PFS) to D (P<0.001), and better to CZ (P=0.002). High ESRP1 expression was independently associated with longer PSA-PFS (P<0.001) and radiologic-PFS (P=0.001) in D and shorter PSA-PFS in the CZ cohort (P=0.041). High SYP expression was independently associated with lower PSA-PFS in D (P=0.003) and overall survival (OS) in CZ (P=0.002), and high EZH2 expression was associated with adverse OS in D-treated patients (P=0.013). In conclusion, EMT profile in primary tumor is differentially associated with D or CZ benefit and NE dedifferentiation correlates with adverse taxanes clinical outcome.

Genomic and epigenomic insights into the origin, pathogenesis, and clinical behavior of mantle cell lymphoma subtypes.

Mantle cell lymphoma (MCL) is a mature B-cell neoplasm initially driven by CCND1 rearrangement with 2 molecular subtypes, conventional MCL (cMCL) and leukemic non-nodal MCL (nnMCL), that differ in their clinicobiological behavior. To identify the genetic and epigenetic alterations determining this diversity, we used whole-genome (n = 61) and exome (n = 21) sequencing (74% cMCL, 26% nnMCL) combined with transcriptome and DNA methylation profiles in the context of 5 MCL reference epigenomes. We identified that open and active chromatin at the major translocation cluster locus might facilitate the t(11;14)(q13;32), which modifies the 3-dimensional structure of the involved regions. This translocation is mainly acquired in precursor B cells mediated by recombination-activating genes in both MCL subtypes, whereas in 8% of cases the translocation occurs in mature B cells mediated by activation-induced cytidine deaminase. We identified novel recurrent MCL drivers, including CDKN1B, SAMHD1, BCOR, SYNE1, HNRNPH1, SMARCB1, and DAZAP1. Complex structural alterations emerge as a relevant early oncogenic mechanism in MCL, targeting key driver genes. Breakage-fusion-bridge cycles and translocations activated oncogenes (BMI1, MIR17HG, TERT, MYC, and MYCN), generating gene amplifications and remodeling regulatory regions. cMCL carried significant higher numbers of structural variants, copy number alterations, and driver changes than nnMCL, with exclusive alterations of ATM in cMCL, whereas TP53 and TERT alterations were slightly enriched in nnMCL. Several drivers had prognostic impact, but only TP53 and MYC aberrations added value independently of genomic complexity. An increasing genomic complexity, together with the presence of breakage-fusion-bridge cycles and high DNA methylation changes related to the proliferative cell history, defines patients with different clinical evolution.

Molecular classification and therapeutic targets in extrahepatic cholangiocarcinoma.

BACKGROUND & AIMS: Cholangiocarcinoma (CCA), a deadly malignancy of the bile ducts, can be classified based on its anatomical location into either intrahepatic (iCCA) or extrahepatic (eCCA), each with different pathogenesis and clinical management. There is limited understanding of the molecular landscape of eCCA and no targeted therapy with clinical efficacy has been approved. We aimed to provide a molecular classification of eCCA and identify potential targets for molecular therapies. METHODS: An integrative genomic analysis of an international multicenter cohort of 189 eCCA cases was conducted. Genomic analysis included whole-genome expression, targeted DNA-sequencing and immunohistochemistry. Molecular findings were validated in an external set of 181 biliary tract tumors from the ICGC. RESULTS: KRAS (36.7%), TP53 (34.7%), ARID1A (14%) and SMAD4 (10.7%) were the most prevalent mutations, with ∼25% of tumors having a putative actionable genomic alteration according to OncoKB. Transcriptome-based unsupervised clustering helped us define 4 molecular classes of eCCA. Tumors classified within the Metabolic class (19%) showed a hepatocyte-like phenotype with activation of the transcription factor HNF4A and enrichment in gene signatures related to bile acid metabolism. The Proliferation class (23%), more common in patients with distal CCA, was characterized by enrichment of MYC targets, ERBB2 mutations/amplifications and activation of mTOR signaling. The Mesenchymal class (47%) was defined by signatures of epithelial-mesenchymal transition, aberrant TGFß signaling and poor overall survival. Finally, tumors in the Immune class (11%) had a higher lymphocyte infiltration, overexpression of PD-1/PD-L1 and molecular features associated with a better response to immune checkpoint inhibitors. CONCLUSION: An integrative molecular characterization identified distinct subclasses of eCCA. Genomic traits of each class provide the rationale for exploring patient stratification and novel therapeutic approaches. LAY SUMMARY: Targeted therapies have not been approved for the treatment of extrahepatic cholangiocarcinoma. We performed a multi-platform molecular characterization of this tumor in a cohort of 189 patients. These analyses revealed 4 novel transcriptome-based molecular classes of extrahepatic cholangiocarcinoma and identified ∼25% of tumors with actionable genomic alterations, which has potential prognostic and therapeutic implications.

Androgen Receptor and Its Splicing Variant 7 Expression in Peripheral Blood Mononuclear Cells and in Circulating Tumor Cells in Metastatic Castration-Resistant Prostate Cancer.

Androgen receptor (AR) signaling remains crucial in castration-resistant prostate cancer (CRPC). Since it is also essential in immune cells, we studied whether the expression of AR full-length (ARFL) and its splicing variant ARV7 in peripheral blood mononuclear cells (PBMC) predicts systemic treatment response in mCRPC in comparison with circulating-tumor cells (CTC). We measured ARFL and ARV7 mRNA in PBMC and CTC from patients prior to receiving abiraterone (AA), enzalutamide (E), or taxanes by a pre-amplification plus quantitative reverse-transcription PCR. They were also tested in LNCaP-ARV7-transfected and in 22RV1 docetaxel-resistant (22RV1DR) cells. We studied 171 PBMC from 136 patients and from 24 non-cancer controls, and 47 CTC from 22 patients. High PBMC ARV7 levels correlated with worse AA/E and better taxane response. In taxane-treated patients high PBMC ARFL also correlated with longer progression-free survival (PFS). High ARV7 and ARFL expression were independently associated with better biochemical-PFS. Conversely, high CTC ARV7 and ARFL correlated with shorter radiological-PFS and overall survival, respectively. High ARV7 in 22RV1DR and LNCaP-ARV7 cells correlated with taxane resistance. In conclusion, ARFL and ARV7 at PBMC or CTC have a different predictive role in the taxane response, suggesting a potential influence of the AR pathway from PBMC in such response modulation.

The proliferative history shapes the DNA methylome of B-cell tumors and predicts clinical outcome.

We report a systematic analysis of the DNA methylation variability in 1,595 samples of normal cell subpopulations and 14 tumor subtypes spanning the entire human B-cell lineage. Differential methylation among tumor entities relates to differences in cellular origin and to de novo epigenetic alterations, which allowed us to build an accurate machine learning-based diagnostic algorithm. We identify extensive patient-specific methylation variability in silenced chromatin associated with the proliferative history of normal and neoplastic B cells. Mitotic activity generally leaves both hyper- and hypomethylation imprints, but some B-cell neoplasms preferentially gain or lose DNA methylation. Subsequently, we construct a DNA methylation-based mitotic clock called epiCMIT, whose lapse magnitude represents a strong independent prognostic variable in B-cell tumors and is associated with particular driver genetic alterations. Our findings reveal DNA methylation as a holistic tracer of B-cell tumor developmental history, with implications in the differential diagnosis and prediction of clinical outcome.

Whole-genome DNA hyper-methylation in iPSC-derived dopaminergic neurons from Parkinson's disease patients.

BACKGROUND: Parkinson's disease (PD) is characterized by the loss of midbrain dopaminergic neurons (DAn). Previously, we described the presence of DNA hyper- and hypo-methylation alterations in induced pluripotent stem cells (iPSC)-derived DAn from PD patients using the Illumina 450K array which prominently covers gene regulatory regions. METHODS: To expand and contextualize previous findings, we performed the first whole-genome DNA bisulfite sequencing (WGBS) using iPSC-derived DAn from representative PD subjects: one sporadic PD (sPD) patient, one monogenic LRRK2-associated PD patient (L2PD), and one control. RESULTS: At the whole-genome level, we detected global DNA hyper-methylation in the PD which was similarly spread across the genome in both sPD and L2PD and mostly affected intergenic regions. CONCLUSION: This study implements previous epigenetic knowledge in PD at a whole genome level providing the first comprehensive and unbiased CpG DNA methylation data using iPSC-derived DAn from PD patients. Our results indicate that DAn from monogenic or sporadic PD exhibit global DNA hyper-methylation changes. Findings from this exploratory study are to be validated in further studies analyzing other PD cell models and patient tissues.

Parkinsonian traits in amyotrophic lateral sclerosis (ALS): a prospective population-based study.

BACKGROUND: Amyotrophic lateral sclerosis (ALS) is characterized by a spectrum of phenotypes, but only a few studies have addressed the presence of parkinsonian (PK) symptoms. The aim of our study was to investigate the occurrence of PK features in a prospective population-based cohort of ALS patients, determining their demographic, clinical, neuropsychological and genetic characteristics, and identifying their morphological and functional imaging correlates. METHODS: A consecutive series of ALS patients were enrolled and prospectively followed for 2 years. Patients were classified according to the presence (ALS-PK) or absence (ALS) of PK signs, and they underwent neuropsychological testing, genetic analysis for the main ALS and PD genes, brain MRI and 18F-FDG-PET. ALS-PK patients underwent 123I-ioflupane SPECT. RESULTS: Out of 114 eligible patients, 101 (64 men; mean age at onset 65.1 years) were recruited. Thirty-one patients (30.7%) were classified as ALS-PK. Compared to ALS patients, ALS-PK patients were more frequently male, but did not differ for any other clinical, demographic or neuropsychological factors. 123I-ioflupane SPECT was normal in all but two ALS-PK patients. At 18F-FDG-PET, ALS-PK patients showed a relative hypometabolism in left cerebellum and a relatively more preserved metabolism in right insula and frontal regions; MRI fractional anisotropy was reduced in the sagittal stratum and increased in the retrolenticular part of the internal capsule. CONCLUSIONS: In our study, about 30% of ALS patients showed PK signs. Neuroimaging data indicate that PK signs are due to the involvement of brain circuitries other than classical nigrostriatal ones, strengthening the hypothesis of ALS as a complex multisystem disease.

Differential expression of long non-coding RNAs are related to proliferation and histological diversity in follicular lymphomas.

Long non-coding RNAs (lncRNAs) comprise a family of non-coding transcripts that are emerging as relevant gene expression regulators of different processes, including tumour development. To determine the possible contribution of lncRNA to the pathogenesis of follicular lymphoma (FL) we performed RNA-sequencing at high depth sequencing in primary FL samples ranging from grade 1-3A to aggressive grade 3B variants using unpurified (n = 16) and purified (n = 12) tumour cell suspensions from nodal samples. FL grade 3B had a significantly higher number of differentially expressed lncRNAs (dif-lncRNAs) with potential target coding genes related to cell cycle regulation. Nine out of the 18 selected dif-lncRNAs were validated by quantitative real time polymerase chain reaction in an independent series (n = 43) of FL. RP4-694A7.2 was identified as the top deregulated lncRNA potentially involved in cell proliferation. RP4-694A7.2 silencing in the WSU-FSCCL FL cell line reduced cell proliferation due to a block in the G1/S phase. The relationship between RP4-694A7.2 and proliferation was confirmed in primary samples as its expression levels positively related to the Ki-67 proliferation index. In summary, lncRNAs are differentially expressed across the clinico-biological spectrum of FL and a subset of them, related to cell cycle, may participate in cell proliferation regulation in these tumours.

Cyclin D1 overexpression induces global transcriptional downregulation in lymphoid neoplasms.

Cyclin D1 is an oncogene frequently overexpressed in human cancers that has a dual function as cell cycle and transcriptional regulator, although the latter is widely unexplored. Here, we investigated the transcriptional role of cyclin D1 in lymphoid tumor cells with cyclin D1 oncogenic overexpression. Cyclin D1 showed widespread binding to the promoters of most actively transcribed genes, and the promoter occupancy positively correlated with the transcriptional output of targeted genes. Despite this association, the overexpression of cyclin D1 in lymphoid cells led to a global transcriptional downmodulation that was proportional to cyclin D1 levels. This cyclin D1-dependent global transcriptional downregulation was associated with a reduced nascent transcription and an accumulation of promoter-proximal paused RNA polymerase II (Pol II) that colocalized with cyclin D1. Concordantly, cyclin D1 overexpression promoted an increase in the Poll II pausing index. This transcriptional impairment seems to be mediated by the interaction of cyclin D1 with the transcription machinery. In addition, cyclin D1 overexpression sensitized cells to transcription inhibitors, revealing a synthetic lethality interaction that was also observed in primary mantle cell lymphoma cases. This finding of global transcriptional dysregulation expands the known functions of oncogenic cyclin D1 and suggests the therapeutic potential of targeting the transcriptional machinery in cyclin D1-overexpressing tumors.

The reference epigenome and regulatory chromatin landscape of chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is a frequent hematological neoplasm in which underlying epigenetic alterations are only partially understood. Here, we analyze the reference epigenome of seven primary CLLs and the regulatory chromatin landscape of 107 primary cases in the context of normal B cell differentiation. We identify that the CLL chromatin landscape is largely influenced by distinct dynamics during normal B cell maturation. Beyond this, we define extensive catalogues of regulatory elements de novo reprogrammed in CLL as a whole and in its major clinico-biological subtypes classified by IGHV somatic hypermutation levels. We uncover that IGHV-unmutated CLLs harbor more active and open chromatin than IGHV-mutated cases. Furthermore, we show that de novo active regions in CLL are enriched for NFAT, FOX and TCF/LEF transcription factor family binding sites. Although most genetic alterations are not associated with consistent epigenetic profiles, CLLs with MYD88 mutations and trisomy 12 show distinct chromatin configurations. Furthermore, we observe that non-coding mutations in IGHV-mutated CLLs are enriched in H3K27ac-associated regulatory elements outside accessible chromatin. Overall, this study provides an integrative portrait of the CLL epigenome, identifies extensive networks of altered regulatory elements and sheds light on the relationship between the genetic and epigenetic architecture of the disease.