Search

Mutations in epigenetic regulators are common in relapsed pediatric acute lymphoblastic leukemia (ALL). Here, we uncovered the mechanism underlying the relapse of ALL driven by an activating mutation of the NSD2 histone methyltransferase (p.E1099K). Using high-throughput drug screening, we found that NSD2-mutant cells were specifically resistant to glucocorticoids. Correction of this mutation restored glucocorticoid sensitivity. The transcriptional response to glucocorticoids was blocked in NSD2-mutant cells due to depressed glucocorticoid receptor (GR) levels and the failure of glucocorticoids to autoactivate GR expression. Although H3K27me3 was globally decreased by NSD2 p.E1099K, H3K27me3 accumulated at the NR3C1 (GR) promoter. Pretreatment of NSD2 p.E1099K cell lines and patient-derived xenograft samples with PRC2 inhibitors reversed glucocorticoid resistance in vitro and in vivo. PRC2 inhibitors restored NR3C1 autoactivation by glucocorticoids, increasing GR levels and allowing GR binding and activation of proapoptotic genes. These findings suggest a new therapeutic approach to relapsed ALL associated with NSD2 mutation. Significance: NSD2 histone methyltransferase mutations observed in relapsed pediatric ALL drove glucocorticoid resistance by repression of the GR and abrogation of GR gene autoactivation due to accumulation of K3K27me3 at its promoter. Pretreatment with PRC2 inhibitors reversed resistance, suggesting a new therapeutic approach to these patients with ALL. This article is highlighted in the In This Issue feature, p. 1

Multiple myeloma (MM) is an incurable disease, whose clinical heterogeneity makes its management challenging, highlighting the need for biological features to guide improved therapies. Deregulation of specific long non-coding RNAs (lncRNAs) has been shown in MM, nevertheless, the complete lncRNA transcriptome has not yet been elucidated. In this work, we identified 40,511 novel lncRNAs in MM samples. lncRNAs accounted for 82% of the MM transcriptome and were more heterogeneously expressed than coding genes. A total of 10,351 overexpressed and 9,535 downregulated lncRNAs were identified in MM patients when compared with normal bone-marrow plasma cells. Transcriptional dynamics study of lncRNAs in the context of normal B-cell maturation revealed 989 lncRNAs with exclusive expression in MM, among which 89 showed de novo epigenomic activation. Knockdown studies on one of these lncRNAs, SMILO (specific myeloma intergenic long non-coding RNA), resulted in reduced proliferation and induction of apoptosis of MM cells, and activation of the interferon pathway. We also showed that the expression of lncRNAs, together with clinical and genetic risk alterations, stratified MM patients into several progression-free survival and overall survival groups. In summary, our global analysis of the lncRNAs transcriptome reveals the presence of specific lncRNAs associated with the biological and clinical behavior of the disease.

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.

To investigate the three-dimensional (3D) genome architecture across normal B cell differentiation and in neoplastic cells from different subtypes of chronic lymphocytic leukemia and mantle cell lymphoma patients, here we integrate in situ Hi-C and nine additional omics layers. Beyond conventional active (A) and inactive (B) compartments, we uncover a highly-dynamic intermediate compartment enriched in poised and polycomb-repressed chromatin. During B cell development, 28% of the compartments change, mostly involving a widespread chromatin activation from naive to germinal center B cells and a reversal to the naive state upon further maturation into memory B cells. B cell neoplasms are characterized by both entity and subtype-specific alterations in 3D genome organization, including large chromatin blocks spanning key disease-specific genes. This study indicates that 3D genome interactions are extensively modulated during normal B cell differentiation and that the genome of B cell neoplasias acquires a tumor-specific 3D genome architecture., The dynamics of genome architecture during human cell differentiation and upon neoplastic transformation remain poorly characterized. Here, the authors integrate in situ Hi-C and nine additional omic layers to characterize the dynamic changes in 3D genome architecture during normal B cell differentiation and in neoplastic cells from chronic lymphocytic leukemia and mantle cell lymphoma patients.

Clinical and genetic risk factors are currently used in multiple myeloma (MM) to stratify patients and to design specific therapies. However, these systems do not capture the heterogeneity of the disease supporting the development of new prognostic factors. In this study, we identified active promoters and alternative active promoters in 6 different B cell subpopulations, including bone-marrow plasma cells, and 32 MM patient samples, using RNA-seq data. We find that expression initiated at both regular and alternative promoters was specific of each B cell subpopulation or MM plasma cells, showing a remarkable level of consistency with chromatin-based promoter definition. Interestingly, using 595 MM patient samples from the CoMMpass dataset, we observed that the expression derived from some alternative promoters was associated with lower progression-free and overall survival in MM patients independently of genetic alterations. Altogether, our results define cancer-specific alternative active promoters as new transcriptomic features that can provide a new avenue for prognostic stratification possibilities in patients with MM.

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, metabolomic and transcriptomic approach to systematically assess metabolic plasticity in diet-induced obese mice after a combined nutritional and exercise intervention. Although most pathological features were successfully reverted, we observed a high degree of metabolic dysfunction irreversibility in visceral white adipose tissue, characterised by abnormal mitochondrial morphology and functionality. Despite two sequential therapeutic interventions and apparent global phenotypic recovery, obesity specifically triggered in visceral adipose a cascade of events progressing from mitochondrial metabolic and proteostatic defects to widespread cellular stress, which compromises its biosynthetic and recycling capacity. Our data indicate that obesity prompts a lasting metabolic fingerprint that leads to a progressive breakdown of metabolic plasticity in white adipose tissue, becoming a significant milestone in disease progression.

Background: Relapse from acute lymphoblastic leukemia (ALL) is one of the most common causes of pediatric cancer-related death. Early relapse of ALL is associated with recurrent mutations of histone methyltransferase NSD2 (nuclear receptor binding SET domain protein 2), which is specific for H3K36me2, suggesting a link to therapy resistance or other mechanisms underlying relapse such as central neural system (CNS) infiltration. NSD2 p.E1099K affects gene expression through disturbing the balance of H3K36me2/H3K27me3. Using CRISPR/Cas9-edited isogenic ALL cell lines and PDX cells, we found that NSD2 p.E1099K drives oncogenic programming, CNS infiltration and glucocorticoid (GC) resistance. However, the molecular mechanisms underlying the relapse of this subtype of ALL are still under investigation. Aim: To elucidate the epigenetic landscape underlying the mechanism of the relapse of pediatric ALL with NSD2 mutation. Methods: We performed in vivo experiments to observe tumor burden, leukemia cell infiltration and survival of the NOD/SCID mice injected with a NSD2 p.E1099K mutation knock-out SEM cell line and knock-in CEM cell line. We determined transcriptome (RNA-Seq), chromatin accessibility (ATAC-Seq) in isogenic RCH-ACV, SEM, RPMI-8402 and CEM cell lines, transcription factor binding and histone modification (ChIP-Seq) and 3D organization (Hi-C) in RCH-ACV cells. Finally, we integrated analysis of RNA-Seq, ATAC-Seq, ChIP-Seq and Hi-C to comprehensively disclose the epigenetic landscape in ALL with NSD2 mutation. Results: NOD/SCID mice xenografted with NSD2 mutant cells developed high tumor burden and infiltration to spleen, liver and brain while the mice injected with WT cells accumulated significant less tumor cells and had extended survival. RNA-Seq analysis showed that reversion of NSD2 mutation to WT caused more genes downregulated while insertion of NSD2 mutation to WT cells led to more genes upregulated (Mutant vs WT: RCH-ACV 838 vs 494, SEM 1567 vs 1158, RPMI-8402 1922 vs 1745, CEM 1809 vs 1031). 50 upregulated genes and 3 downregulated genes were in common in B-ALL and T-ALL with NSD2 mutation. Most of the upregulated genes correlated with neural development and adhesion which might contribute to CNS infiltration (e.g., NCAM1 and NEO1). A few genes were associated with GC resistance such as decreased NR3C1 and increased NR3C2. Accordingly, ATAC-Seq showed that NSD2 mutant cells had more open chromatin peaks than those of WT (RCH-ACV 4853 vs 3212, SEM 10052 vs 7595, RPMI-8402 20392 vs 12133, CEM 10155 vs 6437). ChIP-Seq revealed general large gains of H3K36me2 in intergenic regions in NSD2 mutant cells. Importantly, genes upregulated with NSD2 mutation (e.g., NCAM1 and NEO1) lost H3K27me3 at promoters but gained H3K36me2 at promoters and whole gene bodies, accompanied with increased H3K27ac at enhancers. Conversely, a small subset of genes gained H3K27me3 and lost H3K36me2 in their promoters. Concentrated H3K36me2 in gene bodies diffused and broadened was less prominent and H3K27me3 accumulation became dominant. This for example was associated with repression of NR3C1 to drive GC resistance of NSD2 mutant cells. Genes upregulated in NSD2 mutant cells were enriched for binding sites for lymphoid transcriptional activators such as EBF1 and IRF2. The promoters of the downregulated genes had motifs for transcription factors poorly expressed in lymphoid cells and were enriched for binding sites for the BCL6 transcriptional repressor. Hi-C analysis revealed 430 topologically associated domains (TADs) with increased loop interactions while 136 TADs with decreased interactions were in NSD2 mutant cells compared to WT cells. Overall, 491 regions switched from compartment A to B and 444 regions switched from B to A in NSD2 mutant cells compared to WT cells. Compartment switching from inactive B to active A correlated with upregulated gene expression levels while the reverse was true for switching from A to B. Increased intra-TAD interactions were linked to upregulated genes while decreased intra-TAD interactions were linked to downregulated genes. Conclusions: The NSD2 mutation led to increased tumor burden, CNS infiltration and glucocorticoid resistance due to dysregulation of epigenetic patterns and 3D chromatin architecture, indicating mechanisms underlying the relapse of pediatric ALL and potential therapeutic targets in ALL with NSD2 mutation. Disclosures Licht: Epizyme: Research Funding.

Robust protocols and automation now enable large-scale single-cell RNA and ATAC sequencing experiments and their application on biobank and clinical cohorts. However, technical biases introduced during sample acquisition can hinder solid, reproducible results and a systematic benchmarking is required before entering large-scale data production. Here, we report the existence and extent of gene expression and chromatin accessibility artifacts introduced during sampling and identify experimental and computational solutions for their prevention.

We report a systematic analysis of the biological and clinical implications of DNA methylation variability in five categories of B-cell tumors derived from B cells spanning the entire maturation spectrum. We used 2056 primary samples including training and validation series and show that 88% of the human DNA methylome is dynamically modulated under normal and neoplastic conditions. B-cell tumors display both epigenetic imprints of their cellular origin and de novo, disease-specific epigenetic alterations that in part are related to differential transcription factor binding. These differential methylation patterns were used by a machine-learning approach to create a diagnostic algorithm that accurately classifies 14 B-cell tumor entities and subtypes with different clinical management. Beyond this, we identified extensive patient-specific epigenetic variability targeting constitutively silenced chromatin regions, a phenomenon we could relate to the proliferative history of normal and neoplastic B cells. We observed that, depending on the maturation stage of the tumor cell of origin, mitotic activity leaves different imprints into the DNA methylome. Subsequently, we constructed a novel DNA methylation-based mitotic clock called epiCMIT (epigenetically-determined Cumulative MIToses), whose lapse magnitude represents a strong independent prognostic variable within specific B-cell tumor subtypes and is associated with particular driver genetic alterations. Our findings reveal DNA methylation as a holistic tracker of B-cell tumor developmental history, with implications in the differential diagnosis and prediction of the outcome of the patients.

Multiple myeloma (MM) is a plasma cell neoplasm associated with a broad variety of genetic lesions. In spite of this genetic heterogeneity, MMs share a characteristic malignant phenotype whose underlying molecular basis remains poorly characterized. In the present study, we examined plasma cells from MM using a multi-epigenomics approach and demonstrated that, when compared to normal B cells, malignant plasma cells showed an extensive activation of regulatory elements, in part affecting coregulated adjacent genes. Among target genes up-regulated by this process, we found members of the NOTCH, NF-kB, MTOR signaling, and TP53 signaling pathways. Other activated genes included sets involved in osteoblast differentiation and response to oxidative stress, all of which have been shown to be associated with the MM phenotype and clinical behavior. We functionally characterized MM-specific active distant enhancers controlling the expression of thioredoxin (TXN), a major regulator of cellular redox status and, in addition, identified PRDM5 as a novel essential gene for MM. Collectively, our data indicate that aberrant chromatin activation is a unifying feature underlying the malignant plasma cell phenotype. This research was funded by the European Union's Seventh Framework Programme through the Blueprint Consortium (grant agreement 282510), Fundació La Marató de TV3, Instituto de Salud Carlos III (ISCIII) PI14/01867, PI16/02024, and PI17/00701, TRASCAN (EPICA), MINECO Explora (RTHALMY), Departamento de Salud del Gobierno de Navarra 40/2016, Gilead Fellowship Program (GLD16/00142), Multiple Myeloma Research Foundation Networks of excellence, the International Myeloma Foundation (Brian van Novis) and the Qatar National Research Fund award 7-916-3-237. Furthermore, the authors acknowledge the support of the Generalitat de Catalunya Suport Grups de Recerca AGAUR 2017-SGR-736 and 2017-SGR-1142, TRASCAN-iMMunocelland European Research Council starting grant (MYELOMANEXT), CIBERONC (CB16/12/00489, CB16/12/00369 and CB16/12/00225), co-financed with FEDER funds, the Accelerator award CRUK/AIRC/AECC joint funder-partnership, as well as NCI R01 CA180475 and a MMRF collaborative grant. R.O. was supported by a FPU Fellowship of the Spanish Government, M.K. by an AOI grant of the Spanish Association Against Cancer and N.R. by the Acció instrumental d'incorporació de scientífics i tecnòlegs PERIS 2016 from the Generalitat de Catalunya. This work was partially developed at the Centro Esther Koplowitz (CEK, Barcelona, Spain). We particularly acknowledge the patients for their participation and the Biobank of the University of Navarra for its collaboration.

Background: In multiple myeloma (MM), inactivating mutations and loss of the histone demethylase KDM6A (also named UTX, Ubiquitously transcribed Tetratricopeptide repeat, X chromosome) locus are found in up to 50% of patients and are associated with poorer prognosis. KDM6A belongs to a family of Jumonji-C (Jmj-C)-containing demethylases that work as a scaffold for a multiprotein complex containing H3K4 specific methyltransferases KMT2D and/or KMT2C (MLL2/3), the histone acetyltransferase CBP/p300 and members of the SWI/SNF chromatin-remodeling complex. In a concerted manner this complex activates enhancers by adding activation marks on histones and removing methylation of lysine 27 on histone H3 (H3K27me) associated with gene repression. The H3K27me3 methyltransferase EZH2 that catalyzes the reverse reaction of KDM6A, is often overexpressed in MM and this also correlates with poor prognosis. As well, KMT2C, KMT2D and SWI/SNF components of its complex are all found significantly mutated in MM which emphasize the importance of enhancer deregulation in myeloma. Aim: Our goal is to understand how loss of KDM6A affects gene transcription networks and enhancer function to affect MM pathogenesis or progression. Method: We explored the transcriptional consequences of KDM6A deficiency in patient tumors using the gene expression dataset included within the Multiple Myeloma Research Foundation (MMRF) Researcher Gateway and in MM cell lines using data from the Cancer Cell Line Encyclopedia (CCLE). We further modeled the loss of KDM6A in female MM cell lines using CRISPR-Cas9 ribonucleotide protein (RNP) technology and by re-expressing KDM6A in cells null for KDM6A. Using ChIP-sequencing analysis in isogenic KDM6A knockout or replete cell lines, we identified KDM6A binding site in MM and explored how KDM6A deficiency affects chromatin structure to regulate a transcriptional network. Results: Epithelial-to-mesenchymal-transition (EMT)-related pathways were downregulated in low KDM6A expressing tumor from the MMRF patient cohort. In both male and female patients, low expression of KDM6A was associated with decreased expression of genes involved in immune recognition. In MM cell lines, we observed a significant positive correlation between KDM6A expression level and expression of the transcriptional regulators of Major Histocompatibility complex I and II (MHCII and MCHII), NLRC5 and CIITA, respectively. With ChIP-sequencing analysis we find that KDM6A directly binds NLRC5 and CIITA genes, among many other genes involved in immune surveillance. KDM6A null MM cells displayed upregulation of pathways promoting EMT and downregulation of genes regulating immune function, such as HLA-A/B/C, CCL3/5, HHLA2 and ITGAL. Moreover, re-expression of KDM6A in KDM6A null cells upregulated the expression of these gene along with surface expression of MHCI proteins. Conclusion: Together, our data demonstrates that KDM6A modulates EMT and directly regulates expression of major regulators of immune surveillance. KDM6A deficiency may drive MM towards an EMT program to promote systemic spread and facilitate escape from the immune surveillance. Disclosures No relevant conflicts of interest to declare.

Common Variable Immunodeficiency (CVID) is characterized by impaired antibody production and poor terminal differentiation of the B cell compartment, yet its pathogenesis is still poorly understood. We first reported the occurrence of epigenetic alterations in CVID by high-throughput methylation analysis in CVID-discordant monozygotic twins. Data from a recent whole DNA methylome analysis throughout different stages of normal B cell differentiation allowed us to design a new experimental approach. We selected CpG sites for analysis based on two criteria: one, CpGs with potential association with the transcriptional status of relevant genes for B cell activation and differentiation; and two, CpGs that undergo significant demethylation from naive to memory B cells in healthy individuals. DNA methylation was analyzed by bisulfite pyrosequencing of specific CpG sites in sorted naive and memory B cell subsets from CVID patients and healthy donors. We observed impaired demethylation in two thirds of the selected CpGs in CVID memory B cells, in genes that govern B cell-specific processes or participate in B cell signaling. The degree of demethylation impairment associated with the extent of the memory B cell reduction. The impaired demethylation in such functionally relevant genes as AICDA in switched memory B cells correlated with a lower proliferative rate. Our new results reinforce the hypothesis of altered demethylation during B cell differentiation as a contributing pathogenic mechanism to the impairment of B cell function and maturation in CVID. In particular, deregulated epigenetic control of AICDA could play a role in the defective establishment of a post-germinal center B cell compartment in CVID.

International audience; Molecular mechanisms underlying terminal differen­tiation of B cells into plasma cells are major determi­nants of adaptive immunity but remain only partially understood. Here we present the transcriptional and epigenomic landscapes of cell subsets arising from activation of human naive B cells and differenti­ation into plasmablasts. Cell proliferation of acti­vated B cells was linked to a slight decrease in DNA methylation levels, but followed by a committal step in which an S phase-synchronized differentia­tion switch was associated with an extensive DNA demethylation and local acquisition of 5-hydroxyme­thylcytosine at enhancers and genes related to plasma cell identity. Downregulation of both TGF­β1/SMAD3 signaling and p53 pathway supported this final step, allowing the emergence of a CD23­-negative subpopulation in transition from B cells to plasma cells. Remarkably, hydroxymethylation of PRDM1, a gene essential for plasma cell fate, was coupled to progression in S phase, revealing an intri­cate connection among cell cycle, DNA (hydroxy) methylation, and cell fate determination.

Purpose: The classification of medulloblastoma into WNT, SHH, group 3, and group 4 subgroups has become of critical importance for patient risk stratification and subgroup-tailored clinical trials. Here, we aimed to develop a simplified, clinically applicable classification approach that can be implemented in the majority of centers treating patients with medulloblastoma. Experimental Design: We analyzed 1,577 samples comprising previously published DNA methylation microarray data (913 medulloblastomas, 457 non-medulloblastoma tumors, 85 normal tissues), and 122 frozen and formalin-fixed paraffin-embedded medulloblastoma samples. Biomarkers were identified applying stringent selection filters and Linear Discriminant Analysis (LDA) method, and validated using DNA methylation microarray data, bisulfite pyrosequencing, and direct-bisulfite sequencing. Results: Using a LDA-based approach, we developed and validated a prediction method (EpiWNT-SHH classifier) based on six epigenetic biomarkers that allowed for rapid classification of medulloblastoma into the clinically relevant subgroups WNT, SHH, and non-WNT/non-SHH with excellent concordance (>99%) with current gold-standard methods, DNA methylation microarray, and gene signature profiling analysis. The EpiWNT-SHH classifier showed high prediction capacity using both frozen and formalin-fixed material, as well as diverse DNA methylation detection methods. Similarly, we developed a classifier specific for group 3 and group 4 tumors, based on five biomarkers (EpiG3-G4) with good discriminatory capacity, allowing for correct assignment of more than 92% of tumors. EpiWNT-SHH and EpiG3-G4 methylation profiles remained stable across tumor primary, metastasis, and relapse samples. Conclusions: The EpiWNT-SHH and EpiG3-G4 classifiers represent a new simplified approach for accurate, rapid, and cost-effective molecular classification of single medulloblastoma DNA samples, using clinically applicable DNA methylation detection methods. Clin Cancer Res; 24(6); 1355–63. ©2018 AACR.

Prospective identification of patients with chronic lymphocytic leukemia (CLL) destined to progress would greatly facilitate their clinical management. Recently, whole-genome DNA methylation analyses identified three clinicobiologic CLL subgroups with an epigenetic signature related to different normal B-cell counterparts. Here, we developed a clinically applicable method to identify these subgroups and to study their clinical relevance. Using a support vector machine approach, we built a prediction model using five epigenetic biomarkers that was able to classify CLL patients accurately into the three subgroups, namely naive B-cell-like, intermediate and memory B-cell-like CLL. DNA methylation was quantified by highly reproducible bisulfite pyrosequencing assays in two independent CLL series. In the initial series (n=211), the three subgroups showed differential levels of IGHV (immunoglobulin heavy-chain locus) mutation (P

Deregulation of long non-coding RNAs (lncRNAs) is emerging as a common feature of different human tumors and their investigation may uncover novel biomarkers and oncogenic mechanisms. Previous studies have suggested that the alteration of some lncRNAs may play an important role in the pathogenesis of multiple myeloma (MM); however, the complete expression landscape of lncRNAs has not been elucidated. In the present work we characterized the lncRNAs transcriptome of MM and determined the potential involvement of lncRNAs in the pathogenesis and clinical behavior of MM. To characterize the MM transcriptome, we performed paired-end strand-specific RNA sequencing (ssRNA-seq) in 38 purified plasma cell (PC) samples from MM patients and in 3 bone marrow PCs (BMPCs) of healthy donors, as well as in distinct normal B-cell populations (Naïve, Centroblasts, Centrocytes, Memory and Tonsilar PCs). We identified 40,511 novel lncRNAs that were expressed, accounting for more than half of MM transcriptome (56%). This group of novel lncRNAs together with previously annotated lncRNAs comprised most (82%) of the MM transcriptome. We studied the transcriptional heterogeneity in MM samples and observed that lncRNAs showed a much more heterogeneous expression than coding genes, suggesting that these elements could contribute to the biological heterogeneity of the disease. Moreover, to determine differentially expressed genes, each MM patient was compared to normal BMPCs, detecting 19,886 lncRNAs deregulated (10,351 overexpressed and 9,535 downregulated) in more than 50% of patients. We then analyzed the transcriptional dynamics of MM considering the different stages of B-cell differentiation and focused on a group of 989 lncRNAs that were upregulated specifically in plasma cells from MM in comparison with the rest of B-cell stages (MM-specific lncRNAs). Next, we aimed to determine whether upregulation of MM-specific lncRNAs in MM was under epigenetic control so we analyzed the distribution of six histone modifications with non-overlapping functions (H3K4me3, H3K4me1, H3K27ac, H3K36me3, H3K27me3, and H3K9me3) of within the lncRNAs of interest by ChIP-seq in MM cases as compared to normal B cell subtypes. We detected 89 lncRNAs with de novo epigenomic activation. These data suggest an epigenetic rewiring in MM where the loci of most MM-specific lncRNAs are in an inactive state in normal cells and become active in MM. We focused on a specific lncRNA, LINC-SMILO, de novo epigenetically active and expressed in MM cells to determine whether upregulation of this lncRNA could play a role in the pathogenesis of the disease. Knockdown of LINC-SMILO in 3 different MM cell lines (MM.1S, MM.1R and KMS-11) using two different shRNAs, resulted in reduced proliferation and induction of apoptosis of myeloma cells. Using low input RNA-seq (MARS-seq), we found that inhibition of LINC-SMILO was associated with activation of ERVs (Endogenous retroviruses) and increase in interferon (IFN) induced genes and activation of IFN pathways, essential for MM cells survival. Finally, we aimed to determine whether the use of specific lncRNAs could improve the current prognostic stratification of MM patients using the IA11 release of CoMMpass data. We analyzed the prognostic value of lncRNAs using COX regression analysis and Backward elimination of Stepwise regression analysis, obtaining that the overexpression of the lncRNA PDLIM1P4 together with 1q amplification and 17p deletion stratified MM patients in three different risk groups (Figure 1). In summary, our study shows that the lncRNA transcriptome is widely altered in MM and suggests that some of the identified lncRNAs have marked prognostic influence and can be used as potential therapeutic targets for MM. Disclosures Paiva: Amgen, Bristol-Myers Squibb, Celgene, Janssen, Merck, Novartis, Roche, and Sanofi; unrestricted grants from Celgene, EngMab, Sanofi, and Takeda; and consultancy for Celgene, Janssen, and Sanofi: Consultancy, Honoraria, Research Funding, Speakers Bureau. San-Miguel:Amgen, Bristol-Myers Squibb, Celgene, Janssen, MSD, Novartis, Roche, Sanofi, and Takeda: Consultancy, Honoraria. Melnick:KDAc Therapeutics: Membership on an entity's Board of Directors or advisory committees; Constellation Pharmaceuticals: Consultancy; Epizyme: Consultancy; Janssenn: Research Funding.

Introduction: Genomic studies of chronic lymphocytic leukemia (CLL) have uncovered >80 potential driver mutations. The vast majority of these mutations affect coding regions, and just two potential drivers have been identified in non-coding elements. Aim: To describe the biological and clinical impact of a recurrent A>C mutation at the third base of the small nuclear RNA U1, the non-coding component of the spliceosome involved in the recognition of the 5' splice site (5'SS). Methods: Whole-genome sequencing (WGS) and RNA-seq from 318 CLL patients were used to identify and characterize a highly recurrent A>C point mutation occurring at position 3 of the U1 snRNA gene (g.3A>C mutation). The U1 wild-type and mutant forms were introduced into three CLL cell lines (JVM3, HG3, MEC1) to validate in vitro the predicted effect of this alteration. We screened two independent cohorts including a total of 1,314 CLL patients for the presence of the mutation using the rhAmp SNP genotyping assay, and integrated the U1 mutational status with well-known driver alterations, IGHV and epigenetic subgroups, and clinical parameters. Results: The U1 mutation was found in 8/78 (10.3%) CLL cases analyzed by WGS. Given its role in 5'SS recognition by base-pairing, we reasoned that this mutation was likely to alter the splicing and expression patterns of CLL. We were able to confirm widespread specific alterations in the transcriptome by comparing RNA-seq data between wild-type and g.3A>C mutated samples. Applying this knowledge to an algorithm aimed to infer the U1 mutational status from expression data, we were able to identify 4 mutated cases among 240 additional cases that had RNA-seq but no WGS. In total, 12/318 (3.8%) CLL patients analyzed by WGS and/or RNA-seq harbored this mutation. This g.3A>C U1 mutation changes the preferential A-U base-pairing between U1 and 5'SS to C-G base-pairing, creating novel splice junctions and altering the splicing pattern of 3,193 introns in 1,519 genes. In addition to altered splicing, 869 genes were differentially expressed between mutated and wild-type cases. We identified specific cancer genes (e.g. MSI2, POLD1, or CD44) and pathways (B-cell receptor signaling, promotion of apoptosis, telomere maintenance, among others) altered by the U1 mutation. To confirm a causal link between this mutation and splicing changes, we introduced exogenous U1 genes with or without the mutation into three cell lines. Subsequent RNA-seq of these cell lines recapitulated the altered splicing and expression patterns observed in CLL patients. We next screened for the presence of the U1 mutation 1,057 patients (cohort 1) using the rhAmp assay and it was found in 30 (2.8%) cases. The distribution of the mutation was similar in Binet stages and CLL vs monoclonal B-cell lymphocytosis. However, the U1 mutation was almost always found in IGHV unmutated CLL (29/30, p=9.0e-11) and within the naïve-like CLL epigenetic subgroup (p=3.7e-7). None of the U1 mutated cases had mutations in the SF3B1 splicing factor. Considering only pre-treatment CLL samples, U1 mutation was associated with a shorter time to first treatment independently of the Binet stage, IGHV mutational status, epigenetic subgroups, and mutations in the well-known CLL drivers SF3B1, NOTCH1, ATMor TP53. In cohort 2 (n=257), this mutation was found in 13 (5.1%) patients, confirming its enrichment in IGHV unmutated cases, naïve-like epigenetic subgroup, and splicing modulation. Despite the relatively small number of pre-treatment samples carrying the U1 mutation (7/178) and short follow-up of the patients (median 2.6 years), the effect of this mutation on time to first treatment in cohort 2 was compatible with the one observed in cohort 1. Finally, we screened for the U1 mutation a cohort of diffuse large B-cell lymphoma (n=108), mantle cell lymphoma (n=101), follicular lymphoma (n=87), splenic marginal zone lymphoma (n=12), acute myeloid leukemia (n=52), and myelodysplastic syndrome (n=67). The mutation was not present in any of the samples analyzed. Conclusions: Here we have reported that the third base of the small nuclear RNA U1 is recurrently mutated in CLL, proved its effect in splicing and gene expression, and shown that this mutation is independently associated with faster disease progression. The g.3A>C U1 mutation represents a novel non-coding driver alteration in CLL with potential clinical and therapeutic implications. Disclosures Ramirez Payer: GILEAD SCIENCES: Research Funding. Terol:Astra Zeneca: Consultancy; Gilead: Research Funding; Abbvie: Consultancy; Janssen: Consultancy, Research Funding; Roche: Consultancy. Lopez-Guillermo:Celgene: Consultancy, Research Funding; Janssen: Research Funding; Roche: Consultancy, Research Funding; Gilead: Consultancy, Research Funding.

Despite improvements in chemotherapy that have increased the 5-year survival rates of pediatric ALL to close to 90%, 15-20% of patients may relapse with a very poor prognosis. Pediatric ALL patients, particularly those in relapse can harbor a specific point mutation (E1099K) in NSD2 (nuclear receptor binding SET domain protein 2) gene, also known as MMSET or WHSC1, which encodes a histone methyl transferase specific for H3K36me2. To understand the biology of mutant NSD2, we used CRISPR-Cas9 gene editing to disrupt the NSD2E1099K mutant allele in B-ALL cell lines (RCH-ACV and SEM) and T-ALL cell line (RPMI-8402) or insert the E1099K mutation into the NSD2WT T-ALL cell line (CEM) and B-ALL cell line (697). Cell lines in which the NSD2E1099K mutant allele is present display increased global levels of H3K36me2 and decreased H3K27me3. NSD2E1099Kcells demonstrate enhanced cell growth, colony formation and migration. NSD2E1099K mutant cell lines assayed by RNA-Seq exhibit an aberrant gene signature, mostly representing gene activation, with activation of signaling pathways, genes implicated in the epithelial mesenchymal transition and prominent expression of neural genes not generally found in hematopoietic tissues. Accordingly, NSD2E1099K cell lines showed prominent tropism to the central neural system in xenografts. To understand why this NSD2 mutations are identified prominently in children who relapse early from therapy for ALL, we performed high-throughput screening in our isogenic cell lines with the National Center for Advancing Translation Science (NCATS) Pharmaceutical Collection and other annotated chemical libraries and found that NSD2E1099K cells are resistant to glucocorticoids (GC) but not to other chemotherapeutic agents used to treat ALL such as vincristine, doxorubicin, cyclophosphamide, methotrexate, and 6-mercaptopurine. Accordingly, patient-derived-xenograft ALL cells with NSD2E1099K mutation were resistant to GC treatment. Reversion of NSD2E1099K mutation to NSD2WT restored GC sensitivity to both B- and T-ALL cell lines, which was accompanied by cell cycle arrest in G1 and induced-apoptosis. Furthermore, knock-in of the NSD2E1099K mutation conferred GC resistance to ALL cell lines by triggering cell cycle progression, proliferation and anti-apoptotic processes. Mice with NSD2E1099K xenografts were completely resistant to GC treatment while treatment of mice injected with isogenic NSD2WT cells led to significant tumor reduction and survival benefit. To illustrate these biological phenotypes and understand the molecular mechanism of GC resistance driven by NSD2E1099Kmutation, we investigated the GC-induced transcriptome, GC receptor (GR) binding sites and related epigenetic changes in isogenic ALL cell lines in response to GC treatment. RNA-Seq showed that GC transcriptional response was almost completely blocked in NSD2E1099K cells, especially in T-ALL cell lines, correlating with their lack of biological response. GC treatment activated apoptotic pathways and downregulated cell cycle and DNA repair pathways only in NSD2WT cells. The critical pro-apoptotic regulators BIM and BMF failed to be activated by GC in NSD2E1099K cells but were prominently activated when the NSD2 mutation was removed. Chromatin immunoprecipitation sequencing (ChIP-Seq) showed that, the NSD2E1099K mutation blocked the ability of GR and CTCF to bind most GC response elements (GREs) such as those within BIM and BMF. While GR binding in NSD2WT cells was accompanied by increased H3K27 acetylation and gene expression, this failed to occur in NSD2 mutant cells. Furthermore, we found that GR RNA and protein levels were repressed in ALL cells expressing NSD2E1099K and GC failed to induce GR expression in these cells. Paradoxically, while H3K27me3 levels were generally decreased in NSD2E1099K cells, we saw increased levels of H3K27me3 at the GRE within the GR gene body where GR itself and CTCF normally bind, suggesting a novel role for the polycomb repressive complex 2 and EZH2 inhibitors for this form of GC resistance. In conclusion, these studies demonstrate that NSD2E1099K mutation may play an important role in treatment failure of pediatric ALL relapse by interfering with the GR expression and its ability to bind and activate key target genes. Gene editing screens are being performed to understand how to overcome this resistance. Disclosures No relevant conflicts of interest to declare.

Multiple myeloma (MM) is a plasma cell neoplasm associated with a broad variety of genetic lesions. In spite of this genetic heterogeneity, MMs share a characteristic malignant phenotype whose underlying molecular basis remains poorly characterized. In the present study, we examined plasma cells from MM using a multi-epigenomics approach and demonstrated that when compared to normal B cells, malignant plasma cells showed an extensive activation of regulatory elements, in part affecting co-regulated adjacent genes. Among target genes upregulated by this process, we found members of the NOTCH, NFkB, mTOR1 signaling and p53 signaling pathways. Other activated genes included sets involved in osteoblast differentiation and response to oxidative stress, all of which have been shown to be associated with the MM phenotype and clinical behavior. We functionally characterized MM specific active distant enhancers controlling the expression of thioredoxin (TXN), a major regulator of cellular redox status, and in addition identified PRDM5 as a novel essential gene for MM. Collectively our data indicates that aberrant chromatin activation is a unifying feature underlying the malignant plasma cell phenotype.

Array-based DNA methylation analysis in classical Hodgkin lymphoma reveals new insights into the mechanisms underlying silencing of B cell-specific genes

Acute lymphoblastic leukemia (ALL) is the most common diagnosed pediatric cancer. Despite improvements in chemotherapy that have increased the 5-year survival rate to close to 90%, 15-20% of these patients may relapse with the majority of such children succumbing to this disease. Pediatric ALL patients, particularly those in relapse can harbor a specific point mutation (E1099K) in NSD2 (nuclear receptor binding SET domain protein 2) gene, also known as MMSET or WHSC1, which encodes a histone methyl transferase specific for H3K36me2. To understand the biological processes mediated by mutant NSD2, we used CRISPR-Cas9 gene editing to disrupt the NSD2E1099K mutant allele in two B-ALL cell lines (RCH-ACV and SEM) and one T-ALL cell line (RPMI-8402) and inserted the E1099K mutation into three ALL cell lines (697, CEM, MOLT4). Cell lines in which the NSD2E1099K mutant allele is present display increased global levels of H3K36me2 and decreased H3K27me3. NSD2E1099Kcells compared to cells in which the mutation is removed demonstrate enhanced cell growth, colony formation and migration. NSD2 mutant cell lines assayed by RNA-Seq exhibit an aberrant gene signature, mostly representing gene activation, with activation of signaling pathways, genes implicated in the epithelial mesenchymal transition and prominent expression of neural genes not generally found in hematopoietic tissues. Accordingly, NSD2E1099K cell lines showed prominent tropism to the central neural system (CNS) in xenografts. The NSD2 mutation is found prominently in children who relapse early from therapy for ALL, and NSD2E1099K cells are particularly resistant to glucocorticoids (GC). Reversion of NSD2E1099K mutation to wild type NSD2 conferred glucocorticoid sensitivity to both B and T cell lines. GC response upon disruption of mutant NSD2 was accompanied by cell cycle arrest and apoptosis. Mice xenografted with NSD2E1099K cells were completely resistant to GC treatment while treatment of mice injected with isogenic NSD2 wild-type cells led to significant tumor reduction and survival extension. RNA-Seq analysis showed that GC transcriptional response was almost completely blocked in NSD2E1099K cells, correlating with their lack of biological response. GC treatment activated apoptotic pathways and downregulated cell cycle and DNA repair pathways only in NSD2 wild-type cells. Furthermore, in NSD2 mutant cells, there was lower basal expression level of glucocorticoid receptor (GR) and GR levels were not significantly induced by GC. Accordingly, after treatment with GC, there was significantly less DNA-binding activity of the GR in NSD2E1099K cells than that of NSD2 wild-type cells. The key pro-apoptotic regulators Bim and BMF failed to be activated by GC in NSD2E1099K cells but were prominently activated when the NSD2 mutation was removed. In conclusion, these studies demonstrate that the NSD2E1099K mutation may play an important role in treatment failure of pediatric ALL relapse by causing GC resistance. Future studies will determine how NSD2 which generally activates genes paradoxically blocks the ability of GC and the GR to induce critical pro-death genes. Disclosures Licht: Celgene: Research Funding.

MicroRNAs (miRNAs) are short non-coding RNA molecules that regulate post-transcriptional gene expression. They influence a wide range of physiological functions, including neuronal processes, and are regulated by various mechanisms, such as DNA methylation. This epigenetic mark is recognized by transcriptional regulators such as the methyl CpG binding protein Mecp2. Rett syndrome is a complex neurological disorder that has been associated with mutations in the gene coding for Mecp2. Thus, we examined the possible miRNA misregulation caused by Mecp2 absence in a mouse model of Rett syndrome. Using miRNA expression microarrays, we observed that the brain of Rett syndrome mice undergoes a disruption of the expression profiles of miRNAs. Among the significantly altered miRNAs (26%, 65 of 245), overall downregulation of these transcripts was the most common feature (71%), whilst the remaining 30% were upregulated. Further validation by quantitative RT-PCR demonstrated that the most commonly disrupted miRNAs were miR-146a, miR-146b, miR-130, miR-122a, miR-342 and miR-409 (downregulated), and miR-29b, miR329, miR-199b, miR-382, miR-296, miR-221 and miR-92 (upregulated). Most importantly, transfection of miR-146a in a neuroblastoma cell line caused the downregulation of IL-1 receptor-associated kinase 1 (Irak1) levels, suggesting that the identified defect of miR-146a in Rett syndrome mice brains might be responsible for the observed upregulation of Irak1 in this model of the human disease. Overall, we provide another level of molecular deregulation occurring in Rett syndrome that might be useful for understanding the disease and for designing targeted therapies.

International audience; We analyzed the DNA methylome of ten subpopulations spanning the entire B cell differentiation program by whole-genome bisulfite sequencing and high-density microarrays. We observed that non-CpG methylation disappeared upon B cell commitment, whereas CpG methylation changed extensively during B cell maturation, showing an accumulative pattern and affecting around 30% of all measured CpG sites. Early differentiation stages mainly displayed enhancer demethylation, which was associated with upregulation of key B cell transcription factors and affected multiple genes involved in B cell biology. Late differentiation stages, in contrast, showed extensive demethylation of heterochromatin and methylation gain at Polycomb-repressed areas, and genes with apparent functional impact in B cells were not affected. This signature, which has previously been linked to aging and cancer, was particularly widespread in mature cells with an extended lifespan. Comparing B cell neoplasms with their normal counterparts, we determined that they frequently acquire methylation changes in regions already undergoing dynamic methylation during normal B cell differentiation.

Although Burkitt lymphomas and follicular lymphomas both have features of germinal center B cells, they are biologically and clinically quite distinct. Here we performed whole-genome bisulfite, genome and transcriptome sequencing in 13 IG-MYC translocation-positive Burkitt lymphoma, nine BCL2 translocation-positive follicular lymphoma and four normal germinal center B cell samples. Comparison of Burkitt and follicular lymphoma samples showed differential methylation of intragenic regions that strongly correlated with expression of associated genes, for example, genes active in germinal center dark-zone and light-zone B cells. Integrative pathway analyses of regions differentially methylated in Burkitt and follicular lymphomas implicated DNA methylation as cooperating with somatic mutation of sphingosine phosphate signaling, as well as the TCF3-ID3 and SWI/SNF complexes, in a large fraction of Burkitt lymphomas. Taken together, our results demonstrate a tight connection between somatic mutation, DNA methylation and transcriptional control in key B cell pathways deregulated differentially in Burkitt lymphoma and other germinal center B cell lymphomas.

This work was funded by the Spanish Ministry of Economy and Competitivity (MINECO) through the Instituto de Salud Carlos III (ISCIII), Red Tematica de Investigacion del Cancer (RTICC) and Instituto Nacional de Bioinformatica (INB) from ISCIII, and Consolider CSD2007-00050. C.L.-O. is an investigator of the Botın Foundation and E.C. of the ICREA-Academia program. N.L.-B., D.T., and A.G.-P. acknowledge funding from the Spanish Ministry of Science and Technology (grant number SAF2009-06954). S.E. is supported by a fellowship from La Caixa., Ferreira, P.G., Jares, P., Rico, D., Gómez-López, G., Martínez-Trillos, A., Villamor, N., Ecker, S., González-Pérez, A., Knowles, D.G., Monlong, J., Johnson, R., Quesada, V., Djebali, S., Papasaikas, P., López-Guerra, M., Colomer, D., Royo, C., Cazorla, M., Pinyol, M., Clot, G., Aymerich, M., Rozman, M., Kulis, M., Tamborero, D., Gouin, A., Blanc, J., Gut, M., Gut, I., Puente, X.S., Pisano, D.G., Martin-Subero, J.I., López-Bigas, N., López-Guillermo, A., Valencia, A., López-Otín, C., Campo, E., Guigó, R.

Increasing amount of evidence indicates that the deregulation of non-coding elements is a common feature of cancer and therefore, its investigation may uncover new molecular oncogenic mechanisms. In multiple myeloma (MM), the altered expression of a small number of long non-coding RNAs (lncRNAs) has been associated with progression and decreased survival, suggesting that these elements may play a more important role in this disease than previously expected. Nevertheless, an extensive high-throughput analysis that characterizes the deregulation of lncRNAs in MM has not yet been performed. To characterize the transcriptome, including all genomic types of lncRNAs, of MM we performed a paired end strand-specific RNA sequencing (ssRNA-seq) in 38 purified plasma cell (PC) samples from MM patients, as well as PC samples from tonsils (TPCs, n=5) and bone marrow (BMPCs, n=3) of healthy donors as controls. Principal component analysis (PCA) demonstrated that normal PC samples from tonsil and bone marrow cluster separately, suggesting that in spite of being the same cell type, their coding and non coding transcriptomes are very different. Therefore, we selected BMPCs as the normal counterparts for comparison with BM of MM samples. PCA analysis also demonstrated that the well known heterogeneity of MM patients rely not only on the coding transcriptome but also on the lncRNA expression profile. Comparison of MM to BMPCs samples showed 70 previously annotated lncRNAs that were deregulated in MM patients, with 3 lncRNAs showing higher and 67 lower expression than normal BMPCs. Moreover, we identified 40.552 novel MM-specific lncRNAs that were present in at least 3 of the 38 patients, highlighting the magnitude of the deregulation of these non coding elements in MM. To determine the functional role of altered lncRNAs in the biology of MM plasma cells we focused on the study of LINC-MSL1 (Myeloma-Specific LncRNA 1). Analysis of the expression of this lncRNA at different stages of B-cell differentiation (Naïve, Germinal Center, Memory and PC) indicated that it is not expressed at any stage, except for a modest expression in BMPCs. Interestingly, its overexpression was detected in 40% of MM specimens when compared to normal BMPCs which was validated by qPCR in an independent cohort of MM patients. To determine whether the expression of this lncRNA is regulated by epigenetic mechanisms, we studied the DNA methylation state of this gene. DNA methylation analysis in MM demonstrated that the CpGs located upstream of LINC-MSL1 were differentially methylated in comparison with normal counterpart BMPC. These CpGs showed 70% DNA methylation in control samples, about 40% in MGUS, whereas the average of MM was about 20%, showing a remarkable hypomethylation. We validated these results by pyrosequencing, which showed a significantly lower DNA methylation at the promoter region in comparison with B cell populations from tonsil, normal BMPCs and cell lines that do not overexpress LINC-MSL1. We also have observed a gain of active chromatin states analyzed by ChiP-seq in the promoter region of LINC-MSL1 in MM patient samples. These data suggest that epigenetic mechanisms, namely the progressive hypomethylation and the gain of active histone modifications, are the cause of the overexpression of LINC-MSL1 in MM. To analyze the role of the overexpression of LINC-MSL1 in MM, we engineered two MM cell lines that show high levels of LINC-MSL1, MM.1S and MM.1R, to express shRNAs against this lncRNA. Knockdown of LINC-MSL1 by two different shRNAs resulted in a reduced proliferation of the cell lines over time. This effect was not associated with a cell cycle arrest but with a marked increased in the percentage of Annexin V-positive apoptotic cells, indicating that the overexpression of LINC-MSL1 is necessary for the survival of MM cells. All together, these data demonstrate that the alteration of lncRNAs is an important an unexplored feature that contributes to MM pathogenesis. The overexpression of LINC-MSL1 is essential for MM survival and is very specific of MM BMPCs, suggesting it could be a relevant therapeutic target. Disclosures Paiva: Celgene: Honoraria, Research Funding; Janssen: Honoraria; Takeda: Honoraria, Research Funding; Sanofi: Consultancy, Research Funding; EngMab: Research Funding; Amgen: Honoraria; Binding Site: Research Funding. Melnick:Janssen: Research Funding.

Neuroblastoma (NB) pathogenesis has been reported to be closely associated with numerous genetic alterations. However, underlying DNA methylation patterns have not been extensively studied in this developmental malignancy. Here, we generated microarray-based DNA methylation profiles of primary neuroblastic tumors. Stringent supervised differential methylation analyses allowed us to identify epigenetic changes characteristic for NB tumors as well as for clinical and biological subtypes of NB. We observed that gene-specific loss of DNA methylation is more prevalent than promoter hypermethylation. Remarkably, such hypomethylation affected cancer-related biological functions and genes relevant to NB pathogenesis such as CCND1, SPRR3, BTC, EGF and FGF6. In particular, differential methylation in CCND1 affected mostly an evolutionary conserved functionally relevant 3′ untranslated region, suggesting that hypomethylation outside promoter regions may play a role in NB pathogenesis. Hypermethylation targeted genes involved in cell development and proliferation such as RASSF1A, POU2F2 or HOXD3, among others. The results derived from this study provide new candidate epigenetic biomarkers associated with NB as well as insights into the molecular pathogenesis of this tumor, which involves a marked gene-specific hypomethylation.

Decapping is an essential step in multiple pathways of mRNA degradation. Previously, we synthesized mRNAs containing caps that were resistant to decapping, both to dissect the various pathways for mRNA degradation and to stabilize mRNA for more sustained protein expression. mRNAs containing an α-β CH2 group are resistant to in vitro cleavage by the decapping enzyme hDcp2 but poorly translated. mRNAs containing an S substitution at the β-phosphate are well translated but only partially resistant to hDcp2. We now describe seven new cap analogs substituted at the β-phosphate with BH3 or Se, or substituted at either the α-β or β-γ O with NH. The analogs differ in affinity for eIF4E and efficiency of in vitro incorporation into mRNA by T7 RNA polymerase. Luciferase mRNAs capped with these analogs differ in resistance to hDcp2 hydrolysis in vitro, translational efficiency in rabbit reticulocyte lysate and in HeLa cells, and stability in HeLa cells. Whereas mRNAs capped with m27,2′-OGppSpG were previously found to have the most favorable properties of translational efficiency and stability in mammalian cells, mRNAs capped with m7GppBH3pm7G are translated with the same efficiency but are more stable. Interestingly, some mRNAs exhibit a lag of up to 60 min before undergoing first-order decay (t1/2 ≅ 25 min). Only mRNAs that are efficiently capped, resistant to decapping in vitro, and actively translated have long lag phases.