Search

Phenotypic plasticity is a cause of glioblastoma therapy failure. We previously showed that suppressing the oligodendrocyte-lineage regulator SOX10 promotes glioblastoma progression. Here, we analyze SOX10-mediated phenotypic plasticity and exploit it for glioblastoma therapy design. We show that low SOX10 expression is linked to neural stem-cell (NSC)-like glioblastoma cell states and is a consequence of temozolomide treatment in animal and cell line models. Single-cell transcriptome profiling of Sox10-KD tumors indicates that Sox10 suppression is sufficient to induce tumor progression to an aggressive NSC/developmental-like phenotype, including a quiescent NSC-like cell population. The quiescent NSC state is induced by temozolomide and Sox10-KD and reduced by Notch pathway inhibition in cell line models. Combination treatment using Notch and HDAC/PI3K inhibitors extends the survival of mice carrying Sox10-KD tumors, validating our experimental therapy approach. In summary, SOX10 suppression mediates glioblastoma progression through NSC/developmental cell-state transition, including the induction of a targetable quiescent NSC state. This work provides a rationale for the design of tumor therapies based on single-cell phenotypic plasticity analysis., (© 2024. The Author(s).)

Abstract: The T-box transcription factor T-bet is known as a master regulator of the T-cell response but its role in malignant B cells has not been sufficiently explored. Here, we conducted single-cell resolved multi-omics analyses of malignant B cells from patients with chronic lymphocytic leukemia (CLL) and studied a CLL mouse model with a genetic knockout of Tbx21. We found that T-bet acts as a tumor suppressor in malignant B cells by decreasing their proliferation rate. NF-κB activity, induced by inflammatory signals provided by the microenvironment, triggered T-bet expression, which affected promoter-proximal and distal chromatin coaccessibility and controlled a specific gene signature by mainly suppressing transcription. Gene set enrichment analysis identified a positive regulation of interferon signaling and negative control of proliferation by T-bet. In line, we showed that T-bet represses cell cycling and is associated with longer overall survival of patients with CLL. Our study uncovered a novel tumor suppressive role of T-bet in malignant B cells via its regulation of inflammatory processes and cell cycling, which has implications for the stratification and therapy of patients with CLL. Linking T-bet activity to inflammation explains the good prognostic role of genetic alterations in the inflammatory signaling pathways in CLL.

In January 2021, the monitoring of circulating variants of SARS-CoV-2 was initiated in Germany under the Corona Surveillance Act, which was discontinued after July 2023. This initiative aimed to enhance pandemic containment, as specific amino acid changes, particularly in the spike protein, were associated with increased transmission and reduced vaccine efficacy. Our group conducted whole genome sequencing using the ARTIC protocol (currently V4) on Illumina's NextSeq 500 platform (and, starting in May 2023, on the MiSeq DX platform) for SARS-CoV-2 positive specimen from patients at Heidelberg University Hospital, associated hospitals, and the public health office in the Rhine-Neckar/Heidelberg region. In total, we sequenced 26,795 SARS-CoV-2-positive samples between January 2021 and July 2023. Valid sequences, meeting the requirements for upload to the German electronic sequencing data hub (DESH) operated by the Robert Koch Institute (RKI), were determined for 24,852 samples, and the lineage/clade could be identified for 25,912 samples. The year 2021 witnessed significant dynamics in the circulating variants in the Rhine-Neckar/Heidelberg region, including A.27.RN, followed by the emergence of B.1.1.7 (Alpha), subsequently displaced by B.1.617.2 (Delta), and the initial occurrences of B.1.1.529 (Omicron). By January 2022, B.1.1.529 had superseded B.1.617.2, dominating with over 90%. The years 2022 and 2023 were then characterized by the dominance of B.1.1.529 and its sublineages, particularly BA.5 and BA.2, and more recently, the emergence of recombinant variants like XBB.1.5. Since the global dominance of B.1.617.2, the identified variant distribution in our local study, apart from a time delay in the spread of new variants, can be considered largely representative of the global distribution. om a time delay in the spread of new variants, can be considered largely representative of the global distribution., Competing Interests: Declaration of competing interest None declared., (Copyright © 2023. Published by Elsevier B.V.)

Recent advances in the genomics of glioblastoma (GBM) led to the introduction of molecular neuropathology but failed to translate into treatment improvement. This is largely attributed to the genetic and phenotypic heterogeneity of GBM, which are considered the major obstacle to GBM therapy. Here, we use advanced human GBM-like organoid (LEGO: Laboratory Engineered Glioblastoma-like Organoid) models and provide an unprecedented comprehensive characterization of LEGO models using single-cell transcriptome, DNA methylome, metabolome, lipidome, proteome, and phospho-proteome analysis. We discovered that genetic heterogeneity dictates functional heterogeneity across molecular layers and demonstrates that NF1 mutation drives mesenchymal signature. Most importantly, we found that glycerol lipid reprogramming is a hallmark of GBM, and several targets and drugs were discovered along this line. We also provide a genotype-based drug reference map using LEGO-based drug screen. This study provides new human GBM models and a research path toward effective GBM therapy., (© 2024. The Author(s).)

Medulloblastomas with extensive nodularity are cerebellar tumors characterized by two distinct compartments and variable disease progression. The mechanisms governing the balance between proliferation and differentiation in MBEN remain poorly understood. Here, we employ a multi-modal single cell transcriptome analysis to dissect this process. In the internodular compartment, we identify proliferating cerebellar granular neuronal precursor-like malignant cells, along with stromal, vascular, and immune cells. In contrast, the nodular compartment comprises postmitotic, neuronally differentiated malignant cells. Both compartments are connected through an intermediate cell stage resembling actively migrating CGNPs. Notably, we also discover astrocytic-like malignant cells, found in proximity to migrating and differentiated cells at the transition zone between the two compartments. Our study sheds light on the spatial tissue organization and its link to the developmental trajectory, resulting in a more benign tumor phenotype. This integrative approach holds promise to explore intercompartmental interactions in other cancers with varying histology., (© 2024. The Author(s).)

Intratumor heterogeneity as a clinical challenge becomes most evident after several treatment lines, when multidrug-resistant subclones accumulate. To address this challenge, the characterization of resistance mechanisms at the subclonal level is key to identify common vulnerabilities. In this study, we integrate whole-genome sequencing, single-cell (sc) transcriptomics (scRNA sequencing), and chromatin accessibility (scATAC sequencing) together with mitochondrial DNA mutations to define subclonal architecture and evolution for longitudinal samples from 15 patients with relapsed or refractory multiple myeloma. We assess transcriptomic and epigenomic changes to resolve the multifactorial nature of therapy resistance and relate it to the parallel occurrence of different mechanisms: (1) preexisting epigenetic profiles of subclones associated with survival advantages, (2) converging phenotypic adaptation of genetically distinct subclones, and (3) subclone-specific interactions of myeloma and bone marrow microenvironment cells. Our study showcases how an integrative multiomics analysis can be applied to track and characterize distinct multidrug-resistant subclones over time for the identification of molecular targets against them., (© 2023 by The American Society of Hematology. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

The location of nucleosomes in the human genome determines the primary chromatin structure and regulates access to regulatory regions. However, genome-wide information on deregulated nucleosome occupancy and its implications in primary cancer cells is scarce. Here, we conducted a genome-wide comparison of high-resolution nucleosome maps in peripheral blood B cells from patients with chronic lymphocytic leukemia (CLL) and healthy individuals at single-base-pair resolution. Our investigation uncovered significant changes of nucleosome positioning in CLL. Globally, the spacing between nucleosomes-the nucleosome repeat length (NRL)-is shortened in CLL. This effect is stronger in the more aggressive IGHV-unmutated CLL subtype than in the IGHV-mutated CLL subtype. Changes in nucleosome occupancy at specific sites are linked to active chromatin remodeling and reduced DNA methylation. Nucleosomes lost or gained in CLL marks differential binding of 3D chromatin organizers such as CTCF as well as immune response-related transcription factors and delineated mechanisms of epigenetic deregulation. The principal component analysis of nucleosome occupancy in cancer-specific regions allowed the classification of samples between cancer subtypes and normal controls. Furthermore, patients could be better assigned to CLL subtypes according to differential nucleosome occupancy than based on DNA methylation or gene expression. Thus, nucleosome positioning constitutes a novel readout to dissect molecular mechanisms of disease progression and to stratify patients. Furthermore, we anticipate that the global nucleosome repositioning detected in our study, such as changes in the NRL, can be exploited for liquid biopsy applications based on cell-free DNA to stratify patients and monitor disease progression., (© 2023 Piroeva et al.; Published by Cold Spring Harbor Laboratory Press.)

Background: Here we present scSNPdemux, a sample demultiplexing pipeline for single-cell RNA sequencing data using natural genetic variations in humans. The pipeline requires alignment files from Cell Ranger (10× Genomics), a population SNP database and genotyped single nucleotide polymorphisms (SNPs) per sample. The tool works on sparse genotyping data in VCF format for sample identification., Results: The pipeline was tested on both single-cell and single-nuclei based RNA sequencing datasets and showed superior demultiplexing performance over the lipid-based CellPlex and Multi-seq sample multiplexing technique which incurs additional single cell library preparation steps. Specifically, our pipeline demonstrated superior sensitivity and specificity in cell-identity assignment over CellPlex, especially on immune cell types with low RNA content., Conclusions: We designed a streamlined pipeline for single-cell sample demultiplexing, aiming to overcome common problems in multiplexing samples using single cell libraries which might affect data quality and can be costly., (© 2023. BioMed Central Ltd., part of Springer Nature.)

In multiple myeloma spatial differences in the subclonal architecture, molecular signatures and composition of the microenvironment remain poorly characterized. To address this shortcoming, we perform multi-region sequencing on paired random bone marrow and focal lesion samples from 17 newly diagnosed patients. Using single-cell RNA- and ATAC-seq we find a median of 6 tumor subclones per patient and unique subclones in focal lesions. Genetically identical subclones display different levels of spatial transcriptional plasticity, including nearly identical profiles and pronounced heterogeneity at different sites, which can include differential expression of immunotherapy targets, such as CD20 and CD38. Macrophages are significantly depleted in the microenvironment of focal lesions. We observe proportional changes in the T-cell repertoire but no site-specific expansion of T-cell clones in intramedullary lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity in multiple myeloma with potential implications for models of cell-cell interactions and disease progression., (© 2023. Springer Nature Limited.)

Single-cell sequencing of RNA (scRNA-seq) has advanced our understanding of cellular heterogeneity and signaling in developmental biology and disease. A large number of complementary assays have been developed to profile transcriptomes of individual cells, also in combination with other readouts, such as chromatin accessibility or antibody-based analysis of protein surface markers. As scRNA-seq technologies are advancing fast, it is challenging to establish robust workflows and up-to-date protocols that are best suited to address the large range of research questions. Here, we review scRNA-seq techniques from mRNA end-counting to total RNA in relation to their specific features and outline the necessary sample preparation steps and quality control measures. Based on our experience in dealing with the continuously growing portfolio from the perspective of a central single-cell facility, we aim to provide guidance on how workflows can be best automatized and share our experience in coping with the continuous expansion of scRNA-seq techniques.

The antiviral response induced by type I interferon (IFN) via the JAK-STAT signaling cascade activates hundreds of IFN-stimulated genes (ISGs) across human and mouse tissues but varies between cell types. However, the links between the underlying epigenetic features and the ISG profile are not well understood. We mapped ISGs, binding sites of the STAT1 and STAT2 transcription factors, chromatin accessibility, and histone H3 lysine modification by acetylation (ac) and mono-/tri-methylation (me1, me3) in mouse embryonic stem cells and fibroblasts before and after IFNβ treatment. A large fraction of ISGs and STAT-binding sites was cell type specific with promoter binding of a STAT1/2 complex being a key driver of ISGs. Furthermore, STAT1/2 binding to putative enhancers induced ISGs as inferred from a chromatin co-accessibility analysis. STAT1/2 binding was dependent on the chromatin context and positively correlated with preexisting H3K4me1 and H3K27ac marks in an open chromatin state, whereas the presence of H3K27me3 had an inhibitory effect. Thus, chromatin features present before stimulation represent an additional regulatory layer for the cell type-specific antiviral response., (© 2023 Muckenhuber et al.)

One goal of regenerative medicine is to rejuvenate tissues and extend lifespan by restoring the function of endogenous aged stem cells. However, evidence that somatic stem cells can be targeted in vivo to extend lifespan is still lacking. Here, we demonstrate that after a short systemic treatment with a specific inhibitor of the small RhoGTPase Cdc42 (CASIN), transplanting aged hematopoietic stem cells (HSCs) from treated mice is sufficient to extend the healthspan and lifespan of aged immunocompromised mice without additional treatment. In detail, we show that systemic CASIN treatment improves strength and endurance of aged mice by increasing the myogenic regenerative potential of aged skeletal muscle stem cells. Further, we show that CASIN modifies niche localization and H4K16ac polarity of HSCs in vivo. Single-cell profiling reveals changes in HSC transcriptome, which underlie enhanced lymphoid and regenerative capacity in serial transplantation assays. Overall, we provide proof-of-concept evidence that a short systemic treatment to decrease Cdc42 activity improves the regenerative capacity of different endogenous aged stem cells in vivo, and that rejuvenated HSCs exert a broad systemic effect sufficient to extend murine health- and lifespan., (© 2022. The Author(s).)

Keratinocyte cancers (KC) are the most prevalent malignancies in fair-skinned populations, posing a significant medical and economic burden to health systems. KC originate in the epidermis and mainly comprise basal cell carcinoma (BCC) and cutaneous squamous cell carcinoma (cSCC). Here, we combined single-cell multi-omics, transcriptomics, and methylomics to investigate the epigenomic dynamics during epidermal differentiation. We identified ~3,800 differentially accessible regions between undifferentiated and differentiated keratinocytes, corresponding to regulatory regions associated with key transcription factors. DNA methylation at these regions defined AK/cSCC subtypes with epidermal stem cell- or keratinocyte-like features. Using cell-type deconvolution tools and integration of bulk and single-cell methylomes, we demonstrate that these subclasses are consistent with distinct cells-of-origin. Further characterization of the phenotypic traits of the subclasses and the study of additional unstratified KC entities uncovered distinct clinical features for the subclasses, linking invasive and metastatic KC cases with undifferentiated cells-of-origin. Our study provides a thorough characterization of the epigenomic dynamics underlying human keratinocyte differentiation and uncovers novel links between KC cells-of-origin and their prognosis., (© 2022 The Authors. Published under the terms of the CC BY 4.0 license.)

Hematopoietic stem cells (HSCs) mediate regeneration of the hematopoietic system following injury, such as following infection or inflammation. These challenges impair HSC function, but whether this functional impairment extends beyond the duration of inflammatory exposure is unknown. Unexpectedly, we observed an irreversible depletion of functional HSCs following challenge with inflammation or bacterial infection, with no evidence of any recovery up to 1 year afterward. HSCs from challenged mice demonstrated multiple cellular and molecular features of accelerated aging and developed clinically relevant blood and bone marrow phenotypes not normally observed in aged laboratory mice but commonly seen in elderly humans. In vivo HSC self-renewal divisions were absent or extremely rare during both challenge and recovery periods. The progressive, irreversible attrition of HSC function demonstrates that temporally discrete inflammatory events elicit a cumulative inhibitory effect on HSCs. This work positions early/mid-life inflammation as a mediator of lifelong defects in tissue maintenance and regeneration., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Upon aging, the function of the intestinal epithelium declines with a concomitant increase in aging-related diseases. ISCs play an important role in this process. It is known that ISC clonal dynamics follow a neutral drift model. However, it is not clear whether the drift model is still valid in aged ISCs. Tracking of clonal dynamics by clonal tracing revealed that aged crypts drift into monoclonality substantially faster than young ones. However, ISC tracing experiments, in vivo and ex vivo, implied a similar clonal expansion ability of both young and aged ISCs. Single-cell RNA sequencing for 1,920 high Lgr5 ISCs from young and aged mice revealed increased heterogeneity among subgroups of aged ISCs. Genes associated with cell adhesion were down-regulated in aged ISCs. ISCs of aged mice indeed show weaker adhesion to the matrix. Simulations applying a single cell-based model of the small intestinal crypt demonstrated an accelerated clonal drift at reduced adhesion strength, implying a central role for reduced adhesion for affecting clonal dynamics upon aging., (© 2022 Hageb et al.)

Virtually all patients with multiple myeloma become unresponsive to treatment over time. Relapsed/refractory multiple myeloma (RRMM) is accompanied by the clonal evolution of myeloma cells with heterogeneous genomic aberrations and profound changes of the bone marrow microenvironment (BME). However, the molecular mechanisms that drive drug resistance remain elusive. Here, we analyze the heterogeneous tumor cell population and its complex interaction network with the BME of 20 RRMM patients by single cell RNA-sequencing before/after treatment. Subclones with chromosome 1q-gain express a specific transcriptomic signature and frequently expand during treatment. Furthermore, RRMM cells shape an immune suppressive BME by upregulation of inflammatory cytokines and close interaction with the myeloid compartment. It is characterized by the accumulation of PD1 + γδ T-cells and tumor-associated macrophages as well as the depletion of hematopoietic progenitors. Thus, our study resolves transcriptional features of subclones in RRMM and mechanisms of microenvironmental reprogramming with implications for clinical decision-making., (© 2021. The Author(s).)

Background & Aims: Primary sclerosing cholangitis (PSC) is characterized by chronic inflammation and progressive fibrosis of the biliary tree. The bile acid receptor TGR5 (GPBAR1) is found on biliary epithelial cells (BECs), where it promotes secretion, proliferation and tight junction integrity. Thus, we speculated that changes in TGR5-expression in BECs may contribute to PSC pathogenesis., Methods: TGR5-expression and -localization were analyzed in PSC livers and liver tissue, isolated bile ducts and BECs from Abcb4 -/- , Abcb4 -/- /Tgr5 Tg and ursodeoxycholic acid (UDCA)- or 24-norursodeoxycholic acid (norUDCA)-fed Abcb4 -/- mice. The effects of IL8/IL8 homologues on TGR5 mRNA and protein levels were studied. BEC gene expression was analyzed by single-cell transcriptomics (scRNA-seq) from distinct mouse models., Results: TGR5 mRNA expression and immunofluorescence staining intensity were reduced in BECs of PSC and Abcb4 -/- livers, in Abcb4 -/- extrahepatic bile ducts, but not in intrahepatic macrophages. No changes in TGR5 BEC fluorescence intensity were detected in liver tissue of other liver diseases, including primary biliary cholangitis. Incubation of BECs with IL8/IL8 homologues, but not with other cytokines, reduced TGR5 mRNA and protein levels. BECs from Abcb4 -/- mice had lower levels of phosphorylated Erk and higher expression levels of Icam1, Vcam1 and Tgfβ2. Overexpression of Tgr5 abolished the activated inflammatory phenotype characteristic of Abcb4 -/- BECs. NorUDCA-feeding restored TGR5-expression levels in BECs in Abcb4 -/- livers., Conclusions: Reduced TGR5 levels in BECs from patients with PSC and Abcb4 -/- mice promote development of a reactive BEC phenotype, aggravate biliary injury and thus contribute to the pathogenesis of sclerosing cholangitis. Restoration of biliary TGR5-expression levels represents a previously unknown mechanism of action of norUDCA., Lay Summary: Primary sclerosing cholangitis (PSC) is a chronic cholestatic liver disease-associated with progressive inflammation of the bile duct, leading to fibrosis and end-stage liver disease. Bile acid (BA) toxicity may contribute to the development and disease progression of PSC. TGR5 is a membrane-bound receptor for BAs, which is found on bile ducts and protects bile ducts from BA toxicity. In this study, we show that TGR5 levels were reduced in bile ducts from PSC livers and in bile ducts from a genetic mouse model of PSC. Our investigations indicate that lower levels of TGR5 in bile ducts may contribute to PSC development and progression. Furthermore, treatment with norUDCA, a drug currently being tested in a phase III trial for PSC, restored TGR5 levels in biliary epithelial cells., Competing Interests: Conflict of interest MR, LS, CK, KF, JS, KD, JH, EL, JV, SFG, GH, SJF, IE, DN, CDF, MGB, GG, SV, JPM, MV, KS, TL, KK, DH, TL, MH have nothing to disclose.VK has served as speaker for Falk Foundation, Albireo and Abbvie. JRH has served on the advisory board for Novartis and Orkla Health. JRH has received research support from Biogen unrelated to the current work. THK has received consulting/speaker reimbursement from Gilead and Intercept unrelated to the current work. CS has served as speaker for Falk Foundation. MT has served as speaker for Falk Foundation, Gilead, Intercept and MSD; he has advised for Albireo, BiomX, Boehringer Ingelheim, Falk Pharma GmbH, Genfit, Gilead, Intercept, Jannsen, MSD, Novartis, Phenex, Regulus and Shire. MT further received travel grants from Abbvie, Falk, Gilead and Intercept and research grants from Albireo, Cymabay, Falk, Gilead, Intercept, MSD and Takeda. MT is also co-inventor of patents on the medical use of NorUDCA filed by the Medical Universities of Graz and Vienna (WO2006119803, WO20099013334). PF received norUDCA, an unrestricted research grant, travel grants and advisory board fees from Dr. Falk Pharma GmbH. PF is listed as co-inventor of patents for the medical use of norUDCA filed by the Medical University of Graz (WO2006119803, WO20099013334). UB received grants for investigator-initiated studies by Dr Falk GmbH, Freiburg, and Intercept, San Diego., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

Background: Cancer evolution depends on epigenetic and genetic diversity. Historically, in multiple myeloma (MM), subclonal diversity and tumor evolution have been investigated mostly from a genetic perspective., Methods: Here, we performed an analysis of 42 MM samples from 21 patients by using enhanced reduced representation bisulfite sequencing (eRRBS). We combined several metrics of epigenetic heterogeneity to analyze DNA methylation heterogeneity in MM patients., Results: We show that MM is characterized by the continuous accumulation of stochastic methylation at the promoters of development-related genes. High combinatorial entropy change is associated with poor outcomes in our pilot study and depends predominantly on partially methylated domains (PMDs). These PMDs, which represent the major source of inter- and intrapatient DNA methylation heterogeneity in MM, are linked to other key epigenetic aberrations, such as CpG island (CGI)/transcription start site (TSS) hypermethylation and H3K27me3 redistribution as well as 3D organization alterations. In addition, transcriptome analysis revealed that intratumor methylation heterogeneity was associated with low-level expression and high variability., Conclusions: We propose that disrupted DNA methylation in MM is responsible for high epigenetic and transcriptomic instability allowing tumor cells to adapt to environmental changes by tapping into a pool of evolutionary trajectories., (© 2021. The Author(s).)

Hepatocellular carcinoma (HCC) can have viral or non-viral causes 1-5 . Non-alcoholic steatohepatitis (NASH) is an important driver of HCC. Immunotherapy has been approved for treating HCC, but biomarker-based stratification of patients for optimal response to therapy is an unmet need 6,7 . Here we report the progressive accumulation of exhausted, unconventionally activated CD8 + PD1 + T cells in NASH-affected livers. In preclinical models of NASH-induced HCC, therapeutic immunotherapy targeted at programmed death-1 (PD1) expanded activated CD8 + PD1 + T cells within tumours but did not lead to tumour regression, which indicates that tumour immune surveillance was impaired. When given prophylactically, anti-PD1 treatment led to an increase in the incidence of NASH-HCC and in the number and size of tumour nodules, which correlated with increased hepatic CD8 + PD1 + CXCR6 + , TOX + , and TNF + T cells. The increase in HCC triggered by anti-PD1 treatment was prevented by depletion of CD8 + T cells or TNF neutralization, suggesting that CD8 + T cells help to induce NASH-HCC, rather than invigorating or executing immune surveillance. We found similar phenotypic and functional profiles in hepatic CD8 + PD1 + T cells from humans with NAFLD or NASH. A meta-analysis of three randomized phase III clinical trials that tested inhibitors of PDL1 (programmed death-ligand 1) or PD1 in more than 1,600 patients with advanced HCC revealed that immune therapy did not improve survival in patients with non-viral HCC. In two additional cohorts, patients with NASH-driven HCC who received anti-PD1 or anti-PDL1 treatment showed reduced overall survival compared to patients with other aetiologies. Collectively, these data show that non-viral HCC, and particularly NASH-HCC, might be less responsive to immunotherapy, probably owing to NASH-related aberrant T cell activation causing tissue damage that leads to impaired immune surveillance. Our data provide a rationale for stratification of patients with HCC according to underlying aetiology in studies of immunotherapy as a primary or adjuvant treatment.

Nonalcoholic steatohepatitis (NASH) is a manifestation of systemic metabolic disease related to obesity, and causes liver disease and cancer 1,2 . The accumulation of metabolites leads to cell stress and inflammation in the liver 3 , but mechanistic understandings of liver damage in NASH are incomplete. Here, using a preclinical mouse model that displays key features of human NASH (hereafter, NASH mice), we found an indispensable role for T cells in liver immunopathology. We detected the hepatic accumulation of CD8 T cells with phenotypes that combined tissue residency (CXCR6) with effector (granzyme) and exhaustion (PD1) characteristics. Liver CXCR6 + CD8 T cells were characterized by low activity of the FOXO1 transcription factor, and were abundant in NASH mice and in patients with NASH. Mechanistically, IL-15 induced FOXO1 downregulation and CXCR6 upregulation, which together rendered liver-resident CXCR6 + CD8 T cells susceptible to metabolic stimuli (including acetate and extracellular ATP) and collectively triggered auto-aggression. CXCR6 + CD8 T cells from the livers of NASH mice or of patients with NASH had similar transcriptional signatures, and showed auto-aggressive killing of cells in an MHC-class-I-independent fashion after signalling through P2X7 purinergic receptors. This killing by auto-aggressive CD8 T cells fundamentally differed from that by antigen-specific cells, which mechanistically distinguishes auto-aggressive and protective T cell immunity.

Tumour heterogeneity encompasses both the malignant cells and their microenvironment. While heterogeneity between individual patients is known to affect the efficacy of cancer therapy, most personalized treatment approaches do not account for intratumour heterogeneity. We addressed this issue by studying the heterogeneity of nodal B-cell lymphomas by single-cell RNA-sequencing and transcriptome-informed flow cytometry. We identified transcriptionally distinct malignant subpopulations and compared their drug-response and genomic profiles. Malignant subpopulations from the same patient responded strikingly differently to anti-cancer drugs ex vivo, which recapitulated subpopulation-specific drug sensitivity during in vivo treatment. Infiltrating T cells represented the majority of non-malignant cells, whose gene-expression signatures were similar across all donors, whereas the frequencies of T-cell subsets varied significantly between the donors. Our data provide insights into the heterogeneity of nodal B-cell lymphomas and highlight the relevance of intratumour heterogeneity for personalized cancer therapy.

The caudal-related homeobox transcription factor CDX2 is expressed in leukemic cells but not during normal blood formation. Retroviral overexpression of Cdx2 induces AML in mice, however the developmental stage at which CDX2 exerts its effect is unknown. We developed a conditionally inducible Cdx2 mouse model to determine the effects of in vivo, inducible Cdx2 expression in hematopoietic stem and progenitor cells (HSPCs). Cdx2-transgenic mice develop myelodysplastic syndrome with progression to acute leukemia associated with acquisition of additional driver mutations. Cdx2-expressing HSPCs demonstrate enrichment of hematopoietic-specific enhancers associated with pro-differentiation transcription factors. Furthermore, treatment of Cdx2 AML with azacitidine decreases leukemic burden. Extended scheduling of low-dose azacitidine shows greater efficacy in comparison to intermittent higher-dose azacitidine, linked to more specific epigenetic modulation. Conditional Cdx2 expression in HSPCs is an inducible model of de novo leukemic transformation and can be used to optimize treatment in high-risk AML.

Fibroblasts are an essential cell population for human skin architecture and function. While fibroblast heterogeneity is well established, this phenomenon has not been analyzed systematically yet. We have used single-cell RNA sequencing to analyze the transcriptomes of more than 5,000 fibroblasts from a sun-protected area in healthy human donors. Our results define four main subpopulations that can be spatially localized and show differential secretory, mesenchymal and pro-inflammatory functional annotations. Importantly, we found that this fibroblast 'priming' becomes reduced with age. We also show that aging causes a substantial reduction in the predicted interactions between dermal fibroblasts and other skin cells, including undifferentiated keratinocytes at the dermal-epidermal junction. Our work thus provides evidence for a functional specialization of human dermal fibroblasts and identifies the partial loss of cellular identity as an important age-related change in the human dermis. These findings have important implications for understanding human skin aging and its associated phenotypes.

Background: In cancer, normal epigenetic patterns are disturbed and contribute to gene expression changes, disease onset, and progression. The cancer epigenome is composed of the epigenetic patterns present in the tumor-initiating cell at the time of transformation, and the tumor-specific epigenetic alterations that are acquired during tumor initiation and progression. The precise dissection of these two components of the tumor epigenome will facilitate a better understanding of the biological mechanisms underlying malignant transformation. Chronic lymphocytic leukemia (CLL) originates from differentiating B cells, which undergo extensive epigenetic programming. This poses the challenge to precisely determine the epigenomic ground state of the cell-of-origin in order to identify CLL-specific epigenetic aberrations., Methods: We developed a linear regression model, methylome-based cell-of-origin modeling (Methyl-COOM), to map the cell-of-origin for individual CLL patients based on the continuum of epigenomic changes during normal B cell differentiation., Results: Methyl-COOM accurately maps the cell-of-origin of CLL and identifies CLL-specific aberrant DNA methylation events that are not confounded by physiologic epigenetic B cell programming. Furthermore, Methyl-COOM unmasks abnormal action of transcription factors, altered super-enhancer activities, and aberrant transcript expression in CLL. Among the aberrantly regulated transcripts were many genes that have previously been implicated in T cell biology. Flow cytometry analysis of these markers confirmed their aberrant expression on malignant B cells at the protein level., Conclusions: Methyl-COOM analysis of CLL identified disease-specific aberrant gene regulation. The aberrantly expressed genes identified in this study might play a role in immune-evasion in CLL and might serve as novel targets for immunotherapy approaches. In summary, we propose a novel framework for in silico modeling of reference DNA methylomes and for the identification of cancer-specific epigenetic changes, a concept that can be broadly applied to other human malignancies.

Tumor-initiating cells are a subpopulation of cells that have self-renewal capacity to regenerate a tumor. Here, we identify stem cell-like chromatin features in human glioblastoma initiating cells (GICs) and link them to a loss of the repressive histone H3 lysine 9 trimethylation (H3K9me3) mark. Increasing H3K9me3 levels by histone demethylase inhibition led to cell death in GICs but not in their differentiated counterparts. The induction of apoptosis was accompanied by a loss of the activating H3 lysine 9 acetylation (H3K9ac) modification and accumulation of DNA damage and downregulation of DNA damage response genes. Upon knockdown of histone demethylases, KDM4C and KDM7A both differentiation and DNA damage were induced. Thus, the H3K9me3-H3K9ac equilibrium is crucial for GIC viability and represents a chromatin feature that can be exploited to specifically target this tumor subpopulation., (© 2019 The Authors. International Journal of Cancer published by John Wiley & Sons Ltd on behalf of UICC.)

Specialized regulatory T (Treg) cells accumulate and perform homeostatic and regenerative functions in nonlymphoid tissues. Whether common precursors for nonlymphoid-tissue Treg cells exist and how they differentiate remain elusive. Using transcription factor nuclear factor, interleukin 3 regulated (Nfil3) reporter mice and single-cell RNA-sequencing (scRNA-seq), we identified two precursor stages of interleukin 33 (IL-33) receptor ST2-expressing nonlymphoid tissue Treg cells, which resided in the spleen and lymph nodes. Global chromatin profiling of nonlymphoid tissue Treg cells and the two precursor stages revealed a stepwise acquisition of chromatin accessibility and reprogramming toward the nonlymphoid-tissue Treg cell phenotype. Mechanistically, we identified and validated the transcription factor Batf as the driver of the molecular tissue program in the precursors. Understanding this tissue development program will help to harness regenerative properties of tissue Treg cells for therapy., Competing Interests: Declaration of Interests The authors declare no competing financial interests., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

The presence of genome-wide DNA hypermethylation is a hallmark of lower grade gliomas (LGG) with isocitrate dehydrogenase (IDH) mutations. Further molecular classification of IDH mutant gliomas is defined by the presence (IDHmut-codel) or absence (IDHmut-noncodel) of hemizygous codeletion of chromosome arms 1p and 19q. Despite the DNA hypermethylation seen in bulk tumors, intra-tumoral heterogeneity at the epigenetic level has not been thoroughly analyzed. To address this question, we performed the first epigenetic profiling of single cells in a cohort of 5 gliomas with IDH1 mutation using single nucleus Assay for Transposase-Accessible Chromatin with high-throughput sequencing (snATAC-seq). Using the Fluidigm HT IFC microfluidics platform, we generated chromatin accessibility maps from 336 individual nuclei, and identified variable promoter accessibility of non-coding RNAs in LGGs. Interestingly, local chromatin structures of several non-coding RNAs are significant factors that contribute to heterogeneity, and show increased promoter accessibility in IDHmut-noncodel samples. As an example for clinical significance of this result, we identify CYTOR as a poor prognosis factor in gliomas with IDH mutation. Open chromatin assay points to differential accessibility of non-coding RNAs as an important source of epigenetic heterogeneity within individual tumors and between molecular subgroups. Rare populations of nuclei that resemble either IDH mutant molecular group co-exist within IDHmut-noncodel and IDHmut-codel groups, and along with non-coding RNAs may be an important issue to consider for future studies, as they may help guide predict treatment response and relapse.A web-based explorer for the data is available at shiny.turcanlab.org.

With ageing, intrinsic haematopoietic stem cell (HSC) activity decreases, resulting in impaired tissue homeostasis, reduced engraftment following transplantation and increased susceptibility to diseases. However, whether ageing also affects the HSC niche, and thereby impairs its capacity to support HSC function, is still widely debated. Here, by using in-vivo long-term label-retention assays we demonstrate that aged label-retaining HSCs, which are, in old mice, the most quiescent HSC subpopulation with the highest regenerative capacity and cellular polarity, reside predominantly in perisinusoidal niches. Furthermore, we demonstrate that sinusoidal niches are uniquely preserved in shape, morphology and number on ageing. Finally, we show that myeloablative chemotherapy can selectively disrupt aged sinusoidal niches in the long term, which is linked to the lack of recovery of endothelial Jag2 at sinusoids. Overall, our data characterize the functional alterations of the aged HSC niche and unveil that perisinusoidal niches are uniquely preserved and thereby protect HSCs from ageing.

Patient-derived 3D cell culture systems are currently advancing cancer research since they potentiate the molecular analysis of tissue-like properties and drug response under well-defined conditions. However, our understanding of the relationship between the heterogeneity of morphological phenotypes and the underlying transcriptome is still limited. To address this issue, we here introduce "pheno-seq" to directly link visual features of 3D cell culture systems with profiling their transcriptome. As prototypic applications breast and colorectal cancer (CRC) spheroids were analyzed by pheno-seq. We identified characteristic gene expression signatures of epithelial-to-mesenchymal transition that are associated with invasive growth behavior of clonal breast cancer spheroids. Furthermore, we linked long-term proliferative capacity in a patient-derived model of CRC to a lowly abundant PROX1-positive cancer stem cell subtype. We anticipate that the ability to integrate transcriptome analysis and morphological patho-phenotypes of cancer cells will provide novel insight on the molecular origins of intratumor heterogeneity.

In chronic lymphocytic leukemia (CLL), a diverse set of genetic mutations is embedded in a deregulated epigenetic landscape that drives cancerogenesis. To elucidate the role of aberrant chromatin features, we mapped DNA methylation, seven histone modifications, nucleosome positions, chromatin accessibility, binding of EBF1 and CTCF, as well as the transcriptome of B cells from CLL patients and healthy donors. A globally increased histone deacetylase activity was detected and half of the genome comprised transcriptionally downregulated partially DNA methylated domains demarcated by CTCF CLL samples displayed a H3K4me3 redistribution and nucleosome gain at promoters as well as changes of enhancer activity and enhancer linkage to target genes. A DNA binding motif analysis identified transcription factors that gained or lost binding in CLL at sites with aberrant chromatin features. These findings were integrated into a gene regulatory enhancer containing network enriched for B-cell receptor signaling pathway components. Our study predicts novel molecular links to targets of CLL therapies and provides a valuable resource for further studies on the epigenetic contribution to the disease., (© 2019 The Authors. Published under the terms of the CC BY 4.0 license.)

The function of somatic stem cells declines with age. Understanding the molecular underpinnings of this decline is key to counteract age-related disease. Here, we report a dramatic drop in the neural stem cells (NSCs) number in the aging murine brain. We find that this smaller stem cell reservoir is protected from full depletion by an increase in quiescence that makes old NSCs more resistant to regenerate the injured brain. Once activated, however, young and old NSCs show similar proliferation and differentiation capacity. Single-cell transcriptomics of NSCs indicate that aging changes NSCs minimally. In the aging brain, niche-derived inflammatory signals and the Wnt antagonist sFRP5 induce quiescence. Indeed, intervention to neutralize them increases activation of old NSCs during homeostasis and following injury. Our study identifies quiescence as a key feature of old NSCs imposed by the niche and uncovers ways to activate NSCs to repair the aging brain., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Biallelic inactivation of SMARCB1, encoding a member of the SWI/SNF chromatin remodeling complex, is the hallmark genetic aberration of atypical teratoid rhabdoid tumors (ATRT). Here, we report how loss of SMARCB1 affects the epigenome in these tumors. Using chromatin immunoprecipitation sequencing (ChIP-seq) on primary tumors for a series of active and repressive histone marks, we identified the chromatin states differentially represented in ATRTs compared with other brain tumors and non-neoplastic brain. Re-expression of SMARCB1 in ATRT cell lines enabled confirmation of our genome-wide findings for the chromatin states. Additional generation of ChIP-seq data for SWI/SNF and Polycomb group proteins and the transcriptional repressor protein REST determined differential dependencies of SWI/SNF and Polycomb complexes in regulation of diverse gene sets in ATRTs., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Processes like cellular senescence are characterized by complex events giving rise to heterogeneous cell populations. However, the early molecular events driving this cascade remain elusive. We hypothesized that senescence entry is triggered by an early disruption of the cells' three-dimensional (3D) genome organization. To test this, we combined Hi-C, single-cell and population transcriptomics, imaging, and in silico modeling of three distinct cells types entering senescence. Genes involved in DNA conformation maintenance are suppressed upon senescence entry across all cell types. We show that nuclear depletion of the abundant HMGB2 protein occurs early on the path to senescence and coincides with the dramatic spatial clustering of CTCF. Knocking down HMGB2 suffices for senescence-induced CTCF clustering and for loop reshuffling, while ectopically expressing HMGB2 rescues these effects. Our data suggest that HMGB2-mediated genomic reorganization constitutes a primer for the ensuing senescent program., (Copyright © 2018 Elsevier Inc. All rights reserved.)

This corrects the article DOI: 10.1038/ni.3799.

This corrects the article DOI: 10.1038/ng.3889.

Regulatory T cells (T reg cells) perform two distinct functions: they maintain self-tolerance, and they support organ homeostasis by differentiating into specialized tissue T reg cells. We found that epigenetic modifications defined the molecular characteristics of tissue T reg cells. Tagmentation-based whole-genome bisulfite sequencing revealed more than 11,000 regions that were methylated differentially in pairwise comparisons of tissue T reg cell populations and lymphoid T cells. Similarities in the epigenetic landscape led to the identification of a common tissue T reg cell population that was present in many organs and was characterized by gain and loss of DNA methylation that included many gene sites associated with the T H 2 subset of helper T cells, such as the gene encoding cytokine IL-33 receptor ST2, as well as the production of tissue-regenerative factors. Furthermore, the ST2-expressing population was dependent on the transcriptional regulator BATF and could be expanded by IL-33. Thus, tissue T reg cells integrate multiple waves of epigenetic reprogramming that define their tissue-restricted specialization.

Several mechanisms of action have been proposed for DNA methyltransferase and histone deacetylase inhibitors (DNMTi and HDACi), primarily based on candidate-gene approaches. However, less is known about their genome-wide transcriptional and epigenomic consequences. By mapping global transcription start site (TSS) and chromatin dynamics, we observed the cryptic transcription of thousands of treatment-induced non-annotated TSSs (TINATs) following DNMTi and HDACi treatment. The resulting transcripts frequently splice into protein-coding exons and encode truncated or chimeric ORFs translated into products with predicted abnormal or immunogenic functions. TINAT transcription after DNMTi treatment coincided with DNA hypomethylation and gain of classical promoter histone marks, while HDACi specifically induced a subset of TINATs in association with H2AK9ac, H3K14ac, and H3K23ac. Despite this mechanistic difference, both inhibitors convergently induced transcription from identical sites, as we found TINATs to be encoded in solitary long terminal repeats of the ERV9/LTR12 family, which are epigenetically repressed in virtually all normal cells.

Cell differentiation is associated with changes in chromatin organization and gene expression. In this study, we examine chromatin structure following differentiation of the human myeloid leukemia cell line (HL-60/S4) into granulocytes with retinoic acid (RA) or into macrophage with phorbol ester (TPA). We performed ChIP-seq of histone H3 and its modifications, analyzing changes in nucleosome occupancy, nucleosome repeat length, eu-/heterochromatin redistribution and properties of epichromatin (surface chromatin adjacent to the nuclear envelope). Nucleosome positions changed genome-wide, exhibiting a specific class of alterations involving nucleosome loss in extended (∼1kb) regions, pronounced in enhancers and promoters. Genes that lost nucleosomes at their promoters showed a tendency to be upregulated. On the other hand, nucleosome gain did not show simple effects on transcript levels. The average genome-wide nucleosome repeat length (NRL) did not change significantly with differentiation. However, we detected an approximate 10 bp NRL decrease around the haematopoietic transcription factor (TF) PU.1 and the architectural protein CTCF, suggesting an effect on NRL proximal to TF binding sites. Nucleosome occupancy changed in regions associated with active promoters in differentiated cells, compared with untreated HL-60/S4 cells. Epichromatin regions revealed an increased GC content and high nucleosome density compared with surrounding chromatin. Epichromatin showed depletion of major histone modifications and revealed enrichment with PML body-associated genes. In general, chromatin changes during HL-60/S4 differentiation appeared to be more localized to regulatory regions, compared with genome-wide changes among diverse cell types studied elsewhere.

Epigenetic modifications and other chromatin features partition the genome on multiple length scales. They define chromatin domains with distinct biological functions that come in sizes ranging from single modified DNA bases to several megabases in the case of heterochromatic histone modifications. Due to chromatin folding, domains that are well separated along the linear nucleosome chain can form long-range interactions in three-dimensional space. It has now become a routine task to map epigenetic marks and chromatin structure by deep sequencing methods. However, assessing and comparing the properties of chromatin domains and their positional relationships across data sets without a priori assumptions remains challenging. Here, we introduce multiscale correlation evaluation (MCORE), which uses the fluctuation spectrum of mapped sequencing reads to quantify and compare chromatin patterns over a broad range of length scales in a model-independent manner. We applied MCORE to map the chromatin landscape in mouse embryonic stem cells and differentiated neural cells. We integrated sequencing data from chromatin immunoprecipitation, RNA expression, DNA methylation, and chromosome conformation capture experiments into network models that reflect the positional relationships among these features on different genomic scales. Furthermore, we used MCORE to compare our experimental data to models for heterochromatin reorganization during differentiation. The application of correlation functions to deep sequencing data complements current evaluation schemes and will support the development of quantitative descriptions of chromatin networks., (Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Epigenetic mechanisms have emerged as links between prenatal environmental exposure and disease risk later in life. Here, we studied epigenetic changes associated with maternal smoking at base pair resolution by mapping DNA methylation, histone modifications, and transcription in expectant mothers and their newborn children. We found extensive global differential methylation and carefully evaluated these changes to separate environment associated from genotype-related DNA methylation changes. Differential methylation is enriched in enhancer elements and targets in particular "commuting" enhancers having multiple, regulatory interactions with distal genes. Longitudinal whole-genome bisulfite sequencing revealed that DNA methylation changes associated with maternal smoking persist over years of life. Particularly in children prenatal environmental exposure leads to chromatin transitions into a hyperactive state. Combined DNA methylation, histone modification, and gene expression analyses indicate that differential methylation in enhancer regions is more often functionally translated than methylation changes in promoters or non-regulatory elements. Finally, we show that epigenetic deregulation of a commuting enhancer targeting c-Jun N-terminal kinase 2 (JNK2) is linked to impaired lung function in early childhood., (© 2016 The Authors. Published under the terms of the CC BY 4.0 license.)

Measuring concentrations of the differentiation-promoting hormone retinoic acid (RA) in glioblastoma tissues would help to understand the reason why RA treatment has been inefficient in clinical trials involving brain tumor patients. Here, we apply a recently established extraction and measurement protocol to screen glioblastoma tissues for the levels of the RA precursor retinol and biologically active RA. Combining this approach with mRNA analyses of 26 tumors and 8 normal brains, we identify a multifaceted disturbance of RA synthesis in glioblastoma, involving multiple aldehyde dehydrogenase 1 family and retinol dehydrogenase enzymes. Through database studies and methylation analyses, we narrow down chromosomal deletions and aberrant promoter hypermethylation as potential mechanisms accounting for these alterations. Employing chromatin immunoprecipitation analyses and cell-culture studies, we further show that chromatin at RA target genes is poised to RA substitution, but most glioblastoma cell cultures are completely resistant to RA treatment. This paradoxical RA response is unrelated to alternative RA signaling through the fatty acid-binding protein 5/peroxisome proliferator-activated receptor delta axis. Our data suggest a multifaceted disturbance of RA synthesis in glioblastoma and contribute to reconsider current RA treatment strategies., (© 2015 Wiley Periodicals, Inc.)

Background: Pluripotent embryonic stem cells (ESCs) have the unique ability to differentiate into every cell type and to self-renew. These characteristics correlate with a distinct nuclear architecture, epigenetic signatures enriched for active chromatin marks and hyperdynamic binding of structural chromatin proteins. Recently, several chromatin-related proteins have been shown to regulate ESC pluripotency and/or differentiation, yet the role of the major heterochromatin proteins in pluripotency is unknown., Results: Here we identify Heterochromatin Protein 1β (HP1β) as an essential protein for proper differentiation, and, unexpectedly, for the maintenance of pluripotency in ESCs. In pluripotent and differentiated cells HP1β is differentially localized and differentially associated with chromatin. Deletion of HP1β, but not HP1α, in ESCs provokes a loss of the morphological and proliferative characteristics of embryonic pluripotent cells, reduces expression of pluripotency factors and causes aberrant differentiation. However, in differentiated cells, loss of HP1β has the opposite effect, perturbing maintenance of the differentiation state and facilitating reprogramming to an induced pluripotent state. Microscopy, biochemical fractionation and chromatin immunoprecipitation reveal a diffuse nucleoplasmic distribution, weak association with chromatin and high expression levels for HP1β in ESCs. The minor fraction of HP1β that is chromatin-bound in ESCs is enriched within exons, unlike the situation in differentiated cells, where it binds heterochromatic satellite repeats and chromocenters., Conclusions: We demonstrate an unexpected duality in the role of HP1β: it is essential in ESCs for maintaining pluripotency, while it is required for proper differentiation in differentiated cells. Thus, HP1β function both depends on, and regulates, the pluripotent state.

Non-coding RNAs can modulate histone modifications that, at the same time, affect transcript expression levels. Here, we dissect such a network in mouse embryonic stem cells (ESCs). It regulates the activity of the reverse transcriptase telomerase, which synthesizes telomeric repeats at the chromosome ends. We find that histone H3 serine 10 phosphorylation set by Aurora kinase B (AURKB) in ESCs during the S phase of the cell cycle at centromeric and (sub)telomeric loci promotes the expression of non-coding minor satellite RNA (cenRNA). Inhibition of AURKB induces silencing of cenRNA transcription and establishment of a repressive chromatin state with histone H3 lysine 9 trimethylation and heterochromatin protein 1 accumulation. This process results in a continuous shortening of telomeres. We further show that AURKB interacts with both telomerase and cenRNA and activates telomerase in trans. Thus, in mouse ESCs, telomere maintenance is regulated via expression of cenRNA in a cell-cycle-dependent manner., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

The global impact of DNA methylation on alternative splicing is largely unknown. Using a genome-wide approach in wild-type and methylation-deficient embryonic stem cells, we found that DNA methylation can either enhance or silence exon recognition and affects the splicing of more than 20% of alternative exons. These exons are characterized by distinct genetic and epigenetic signatures. Alternative splicing regulation of a subset of these exons can be explained by heterochromatin protein 1 (HP1), which silences or enhances exon recognition in a position-dependent manner. We constructed an experimental system using site-specific targeting of a methylated/unmethylated gene and demonstrate a direct causal relationship between DNA methylation and alternative splicing. HP1 regulates this gene's alternative splicing in a methylation-dependent manner by recruiting splicing factors to its methylated form. Our results demonstrate DNA methylation's significant global influence on mRNA splicing and identify a specific mechanism of splicing regulation mediated by HP1., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)