Chronic lymphocytic leukaemia (CLL) is characterized by substantial clinical heterogeneity, despite relatively few genetic alterations. To provide a basis for studying epigenome deregulation in CLL, here we present genome-wide chromatin accessibility maps for 88 CLL samples from 55 patients measured by the ATAC-seq assay. We also performed ChIPmentation and RNA-seq profiling for ten representative samples. Based on the resulting data set, we devised and applied a bioinformatic method that links chromatin profiles to clinical annotations. Our analysis identified sample-specific variation on top of a shared core of CLL regulatory regions. IGHV mutation status-which distinguishes the two major subtypes of CLL-was accurately predicted by the chromatin profiles and gene regulatory networks inferred for IGHV-mutated versus IGHV-unmutated samples identified characteristic differences between these two disease subtypes. In summary, we discovered widespread heterogeneity in the chromatin landscape of CLL, established a community resource for studying epigenome deregulation in leukaemia and demonstrated the feasibility of large-scale chromatin accessibility mapping in cancer cohorts and clinical research.

Genome-wide DNA methylation mapping uncovers epigenetic changes associated with animal development, environmental adaptation, and species evolution. To address the lack of high-throughput methods for DNA methylation analysis in non-model organisms, we developed an integrated approach for studying DNA methylation differences independent of a reference genome. Experimentally, our method relies on an optimized 96-well protocol for reduced representation bisulfite sequencing (RRBS), which we have validated in nine species (human, mouse, rat, cow, dog, chicken, carp, sea bass, and zebrafish). Bioinformatically, we developed the RefFreeDMA software to deduce ad hoc genomes directly from RRBS reads and to pinpoint differentially methylated regions between samples or groups of individuals (http://RefFreeDMA.computational-epigenetics.org). The identified regions are interpreted using motif enrichment analysis and/or cross-mapping to annotated genomes. We validated our method by reference-free analysis of cell-type-specific DNA methylation in the blood of human, cow, and carp. In summary, we present a cost-effective method for epigenome analysis in ecology and evolution, which enables epigenome-wide association studies in natural populations and species without a reference genome.

A large subset of mammalian imprinted genes show extra-embryonic lineage (EXEL) specific imprinted expression that is restricted to placental trophectoderm lineages and to visceral yolk sac endoderm (ysE). Isolated ysE provides a homogenous in vivo model of a mid-gestation extra-embryonic tissue to examine the mechanism of EXEL-specific imprinted gene silencing, but an in vitro model of ysE to facilitate more rapid and cost-effective experiments is not available. Reports indicate that ES cells differentiated into cystic embryoid bodies (EBs) contain ysE, so here we investigate if cystic EBs model ysE imprinted expression. The imprinted expression pattern of cystic EBs is shown to resemble fetal liver and not ysE. To investigate the reason for this we characterized the methylome and transcriptome of cystic EBs in comparison to fetal liver and ysE, by whole genome bisulphite sequencing and RNA-seq. Cystic EBs show a fetal liver pattern of global hypermethylation and low expression of repeats, while ysE shows global hypomethylation and high expression of IAPEz retroviral repeats, as reported for placenta. Transcriptome analysis confirmed that cystic EBs are more similar to fetal liver than ysE and express markers of early embryonic endoderm. Genome-wide analysis shows that ysE shares epigenetic and repeat expression features with placenta. Contrary to previous reports, we show that cystic EBs do not contain ysE, but are more similar to the embryonic endoderm of fetal liver. This explains why cystic EBs reproduce the imprinted expression seen in the embryo but not that seen in the ysE.

The pluripotent ground state is defined as a basal state free of epigenetic restrictions, which influence lineage specification. While naive embryonic stem cells (ESCs) can be maintained in a hypomethylated state with open chromatin when grown using two small-molecule inhibitors (2i)/leukemia inhibitory factor (LIF), in contrast to serum/LIF-grown ESCs that resemble early post-implantation embryos, broader features of the ground-state pluripotent epigenome are not well understood. We identified epigenetic features of mouse ESCs cultured using 2i/LIF or serum/LIF by proteomic profiling of chromatin-associated complexes and histone modifications. Polycomb-repressive complex 2 (PRC2) and its product H3K27me3 are highly abundant in 2i/LIF ESCs, and H3K27me3 is distributed genome-wide in a CpG-dependent fashion. Consistently, PRC2-deficient ESCs showed increased DNA methylation at sites normally occupied by H3K27me3 and increased H4 acetylation. Inhibiting DNA methylation in PRC2-deficient ESCs did not affect their viability or transcriptome. Our findings suggest a unique H3K27me3 configuration protects naive ESCs from lineage priming, and they reveal widespread epigenetic crosstalk in ground-state pluripotency.

Single-cell RNA-seq (scRNA-seq) is emerging as a powerful tool to dissect cell-specific effects of drug treatment in complex tissues. This application requires high levels of precision, robustness, and quantitative accuracy-beyond those achievable with existing methods for mainly qualitative single-cell analysis. Here, we establish the use of standardized reference cells as spike-in controls for accurate and robust dissection of single-cell drug responses.

It was recently demonstrated that embryonic glucagon-producing cells in the pancreas can regenerate and convert into insulin-producing β-like cells through the constitutive/ectopic expression of the Pax4 gene. However, whether α cells in adult mice display the same plasticity is unknown. Similarly, the mechanisms underlying such reprogramming remain unclear. We now demonstrate that the misexpression of Pax4 in glucagon(+) cells age-independently induces their conversion into β-like cells and their glucagon shortage-mediated replacement, resulting in islet hypertrophy and in an unexpected islet neogenesis. Combining several lineage-tracing approaches, we show that, upon Pax4-mediated α-to-β-like cell conversion, pancreatic duct-lining precursor cells are continuously mobilized, re-express the developmental gene Ngn3, and successively adopt a glucagon(+) and a β-like cell identity through a mechanism involving the reawakening of the epithelial-to-mesenchymal transition. Importantly, these processes can repeatedly regenerate the whole β cell mass and thereby reverse several rounds of toxin-induced diabetes, providing perspectives to design therapeutic regenerative strategies.

The histone acetyl reader bromodomain-containing protein 4 (BRD4) is an important regulator of chromatin structure and transcription, yet factors modulating its activity have remained elusive. Here we describe two complementary screens for genetic and physical interactors of BRD4, which converge on the folate pathway enzyme MTHFD1 (methylenetetrahydrofolate dehydrogenase, cyclohydrolase and formyltetrahydrofolate synthetase 1). We show that a fraction of MTHFD1 resides in the nucleus, where it is recruited to distinct genomic loci by direct interaction with BRD4. Inhibition of either BRD4 or MTHFD1 results in similar changes in nuclear metabolite composition and gene expression; pharmacological inhibitors of the two pathways synergize to impair cancer cell viability in vitro and in vivo. Our finding that MTHFD1 and other metabolic enzymes are chromatin associated suggests a direct role for nuclear metabolism in the control of gene expression.

People living with HIV (PLWH) are at increased risk for developing skin and mucosal malignancies despite systemic reconstitution of CD4+ T cells upon antiretroviral therapy (ART). The underlying mechanism of chronic tissue-related immunodeficiency in HIV is unclear. We found that skin CD4+ tissue-resident memory T (Trm) cells were depleted after HIV infection and replenished only upon early ART initiation. TCR clonal analysis following early ART suggested a systemic origin for reconstituting CD4+ Trm cells. Single-cell RNA sequencing in PLWH that received late ART treatment revealed a loss of CXCR3+ Trm cells and a tolerogenic skin immune environment. Human papilloma virus-induced precancerous lesion biopsies showed reduced CXCR3+ Trm cell frequencies in the mucosa in PLWH versus HIV- individuals. These results reveal an irreversible loss of CXCR3+ Trm cells confined to skin and mucosa in PLWH who received late ART treatment, which may be a precipitating factor in the development of HPV-related cancer.

In the context of type 1 diabetes research and the development of insulin-producing β-cell replacement strategies, whether pancreatic ductal cells retain their developmental capability to adopt an endocrine cell identity remains debated, most likely due to the diversity of models employed to induce pancreatic regeneration. In this work, rather than injuring the pancreas, we developed a mouse model allowing the inducible misexpression of the proendocrine gene Neurog3 in ductal cells in vivo. These animals developed a progressive islet hypertrophy attributed to a proportional increase in all endocrine cell populations. Lineage tracing experiments indicated a continuous neo-generation of endocrine cells exhibiting a ductal ontogeny. Interestingly, the resulting supplementary β-like cells were found to be functional. Based on these findings, we suggest that ductal cells could represent a renewable source of new β-like cells and that strategies aiming at controlling the expression of Neurog3, or of its molecular targets/co-factors, may pave new avenues for the improved treatments of diabetes.

Rhabdomeric opsins (r-opsins) are light sensors in cephalic eye photoreceptors, but also function in additional sensory organs. This has prompted questions on the evolutionary relationship of these cell types, and if ancient r-opsins were non-photosensory. A molecular profiling approach in the marine bristleworm Platynereis dumerilii revealed shared and distinct features of cephalic and non-cephalic r-opsin1-expressing cells. Non-cephalic cells possess a full set of phototransduction components, but also a mechanosensory signature. Prompted by the latter, we investigated Platynereis putative mechanotransducer and found that nompc and pkd2.1 co-expressed with r-opsin1 in TRE cells by HCR RNA-FISH. To further assess the role of r-Opsin1 in these cells, we studied its signaling properties and unraveled that r-Opsin1 is a Gαq-coupled blue light receptor. Profiling of cells from r-opsin1 mutants versus wild-types, and a comparison under different light conditions reveals that in the non-cephalic cells light - mediated by r-Opsin1 - adjusts the expression level of a calcium transporter relevant for auditory mechanosensation in vertebrates. We establish a deep-learning-based quantitative behavioral analysis for animal trunk movements and identify a light- and r-Opsin-1-dependent fine-tuning of the worm's undulatory movements in headless trunks, which are known to require mechanosensory feedback. Our results provide new data on peripheral cell types of likely light sensory/mechanosensory nature. These results point towards a concept in which such a multisensory cell type evolved to allow for fine-tuning of mechanosensation by light. This implies that light-independent mechanosensory roles of r-opsins may have evolved secondarily.

Atopic dermatitis (AD) is a chronic inflammatory skin disease characterized by type 2 cytokine-driven skin inflammation and epithelial barrier dysfunction. The latter is believed to allow the increased penetration of chemicals, toxins, and allergens into the skin. House dust mite allergens, particularly Der p 2, are important triggers in sensitized individuals with AD; the precise actions of these allergens in epithelial biology remain, however, incompletely understood. In this study, we compared the effects of the protein allergen Der p 2 and a mix of non-IgE-reactive Der p 2 peptides on skin cells using patch tests in AD patients and healthy participants. We then analyzed mRNA expression profiles of keratinocytes by single-cell RNA-sequencing. We report that existing barrier deficiencies in the non-lesional skin of AD patients allow deep penetration of Der p 2 and its peptides, leading to local microinflammation. Der p 2 protein specifically upregulated genes involved in the innate immune system, stress, and danger signals in suprabasal KC. Der p 2 peptides further downregulated skin barrier genes, in particular the expression of genes involved in cell-matrix and cell-cell adhesion. Peptides also induced genes involved in hyperproliferation and caused disturbances in keratinocyte differentiation. Furthermore, inflammasome-relevant genes and IL18 were overexpressed, while KRT1 was downregulated. Our data suggest that Der p 2 peptides contribute to AD initiation and exacerbation by augmenting hallmark features of AD, such as skin inflammation, barrier disruption, and hyperplasia of keratinocytes.

Neutrophils are evolutionarily conserved innate immune cells playing pivotal roles in host defense. Zebrafish models have contributed substantially to our understanding of neutrophil functions but similarities to human neutrophil maturation have not been systematically characterized, which limits their applicability to studying human disease. Here we show, by generating and analysing transgenic zebrafish strains representing distinct neutrophil differentiation stages, a high-resolution transcriptional profile of neutrophil maturation. We link gene expression at each stage to characteristic transcription factors, including C/ebp-β, which is important for late neutrophil maturation. Cross-species comparison of zebrafish, mouse, and human samples confirms high molecular similarity of immature stages and discriminates zebrafish-specific from pan-species gene signatures. Applying the pan-species neutrophil maturation signature to RNA-sequencing data from human neuroblastoma patients reveals association between metastatic tumor cell infiltration in the bone marrow and an overall increase in mature neutrophils. Our detailed neutrophil maturation atlas thus provides a valuable resource for studying neutrophil function at different stages across species in health and disease.

The transcriptional response to infection with the bacterium Listeria monocytogenes (Lm) requires cooperative signals of the type I interferon (IFN-I)-stimulated JAK-STAT and proinflammatory NF-κB pathways. Using ChIP-seq analysis, we define genes induced in Lm-infected macrophages through synergistic transcriptional activation by NF-κB and the IFN-I-activated transcription factor ISGF3. Using the Nos2 and IL6 genes as prime examples of this group, we show that NF-κB functions to recruit enzymes that establish histone marks of transcriptionally active genes. In addition, NF-κB regulates transcriptional elongation by employing the mediator kinase module for the recruitment of the pTEFb complex. ISGF3 has a major role in associating the core mediator with the transcription start as a prerequisite for TFIID and RNA polymerase II (Pol II) binding. Our data suggest that the functional cooperation between two major antimicrobial pathways is based on promoter priming by NF-κB and the engagement of the core mediator for Pol II binding by ISGF3.

Methods for single-cell genome and transcriptome sequencing have contributed to our understanding of cellular heterogeneity, whereas methods for single-cell epigenomics are much less established. Here, we describe a whole-genome bisulfite sequencing (WGBS) assay that enables DNA methylation mapping in very small cell populations (μWGBS) and single cells (scWGBS). Our assay is optimized for profiling many samples at low coverage, and we describe a bioinformatic method that analyzes collections of single-cell methylomes to infer cell-state dynamics. Using these technological advances, we studied epigenomic cell-state dynamics in three in vitro models of cellular differentiation and pluripotency, where we observed characteristic patterns of epigenome remodeling and cell-to-cell heterogeneity. The described method enables single-cell analysis of DNA methylation in a broad range of biological systems, including embryonic development, stem cell differentiation, and cancer. It can also be used to establish composite methylomes that account for cell-to-cell heterogeneity in complex tissue samples.

The recent discovery that genetically modified α cells can regenerate and convert into β-like cells in vivo holds great promise for diabetes research. However, to eventually translate these findings to human, it is crucial to discover compounds with similar activities. Herein, we report the identification of GABA as an inducer of α-to-β-like cell conversion in vivo. This conversion induces α cell replacement mechanisms through the mobilization of duct-lining precursor cells that adopt an α cell identity prior to being converted into β-like cells, solely upon sustained GABA exposure. Importantly, these neo-generated β-like cells are functional and can repeatedly reverse chemically induced diabetes in vivo. Similarly, the treatment of transplanted human islets with GABA results in a loss of α cells and a concomitant increase in β-like cell counts, suggestive of α-to-β-like cell conversion processes also in humans. This newly discovered GABA-induced α cell-mediated β-like cell neogenesis could therefore represent an unprecedented hope toward improved therapies for diabetes.

Type 1 diabetes is characterized by the destruction of pancreatic β cells, and generating new insulin-producing cells from other cell types is a major aim of regenerative medicine. One promising approach is transdifferentiation of developmentally related pancreatic cell types, including glucagon-producing α cells. In a genetic model, loss of the master regulatory transcription factor Arx is sufficient to induce the conversion of α cells to functional β-like cells. Here, we identify artemisinins as small molecules that functionally repress Arx by causing its translocation to the cytoplasm. We show that the protein gephyrin is the mammalian target of these antimalarial drugs and that the mechanism of action of these molecules depends on the enhancement of GABAA receptor signaling. Our results in zebrafish, rodents, and primary human pancreatic islets identify gephyrin as a druggable target for the regeneration of pancreatic β cell mass from α cells.

Pancreatic islets of Langerhans contain several specialized endocrine cell types, which are commonly identified by the expression of single marker genes. However, the established marker genes cannot capture the complete spectrum of cellular heterogeneity in human pancreatic islets, and existing bulk transcriptome datasets provide averages across several cell populations. To dissect the cellular composition of the human pancreatic islet and to establish transcriptomes for all major cell types, we performed single-cell RNA sequencing on 70 cells sorted from human primary tissue. We used this dataset to validate previously described marker genes at the single-cell level and to identify specifically expressed transcription factors for all islet cell subtypes. All data are available for browsing and download, thus establishing a useful resource of single-cell expression profiles for endocrine cells in human pancreatic islets.

Hematopoietic stem cells (HSCs) sustain hematopoiesis throughout life. HSCs exit dormancy to restore hemostasis in response to stressful events, such as acute blood loss, and must return to a quiescent state to prevent their exhaustion and resulting bone marrow failure. HSC activation is driven in part through the phosphatidylinositol 3-kinase (PI3K)/AKT/mTORC1 signaling pathway, but less is known about the cell-intrinsic pathways that control HSC dormancy. Here, we delineate an ERK-dependent, rate-limiting feedback mechanism that controls HSC fitness and their re-entry into quiescence. We show that the MEK/ERK and PI3K pathways are synchronously activated in HSCs during emergency hematopoiesis and that feedback phosphorylation of MEK1 by activated ERK counterbalances AKT/mTORC1 activation. Genetic or chemical ablation of this feedback loop tilts the balance between HSC dormancy and activation, increasing differentiated cell output and accelerating HSC exhaustion. These results suggest that MEK inhibitors developed for cancer therapy may find additional utility in controlling HSC activation.

In early-stage mycosis fungoides (MF), the most common primary cutaneous T-cell lymphoma, limited skin involvement with patches and plaques is associated with a favorable prognosis. Nevertheless, approximately 20-30% of cases progress to tumors or erythroderma, resulting in poor outcome. At present, factors contributing to this switch from indolent to aggressive disease are only insufficiently understood.

Methylation of cytosines is a prototypic epigenetic modification of the DNA. It has been implicated in various regulatory mechanisms across the animal kingdom and particularly in vertebrates. We mapped DNA methylation in 580 animal species (535 vertebrates, 45 invertebrates), resulting in 2443 genome-scale DNA methylation profiles of multiple organs. Bioinformatic analysis of this large dataset quantified the association of DNA methylation with the underlying genomic DNA sequence throughout vertebrate evolution. We observed a broadly conserved link with two major transitions-once in the first vertebrates and again with the emergence of reptiles. Cross-species comparisons focusing on individual organs supported a deeply conserved association of DNA methylation with tissue type, and cross-mapping analysis of DNA methylation at gene promoters revealed evolutionary changes for orthologous genes. In summary, this study establishes a large resource of vertebrate and invertebrate DNA methylomes, it showcases the power of reference-free epigenome analysis in species for which no reference genomes are available, and it contributes an epigenetic perspective to the study of vertebrate evolution.