Search

FOXP3+ regulatory T cells (Treg cells) are key for immune homeostasis. Here, we reveal that nuclear receptor corepressor 1 (NCOR1) controls naïve and effector Treg cell states. Upon NCOR1 deletion in T cells, effector Treg cell frequencies were elevated in mice and in in vitro-generated human Treg cells. NCOR1-deficient Treg cells failed to protect mice from severe weight loss and intestinal inflammation associated with CD4+ T cell transfer colitis, indicating impaired suppressive function. NCOR1 controls the transcriptional integrity of Treg cells, since effector gene signatures were already upregulated in naïve NCOR1-deficient Treg cells while effector NCOR1-deficient Treg cells failed to repress genes associated with naïve Treg cells. Moreover, genes related to cholesterol homeostasis including targets of liver X receptor (LXR) were dysregulated in NCOR1-deficient Treg cells. However, genetic ablation of LXRβ in T cells did not revert the effects of NCOR1 deficiency, indicating that NCOR1 controls naïve and effector Treg cell subset composition independent from its ability to repress LXRβ-induced gene expression. Thus, our study reveals that NCOR1 maintains naïve and effector Treg cell states via regulating their transcriptional integrity. We also reveal a critical role for this epigenetic regulator in supporting the suppressive functions of Treg cells in vivo.

Abstract The significance of transient neuropeptide expression during postnatal brain development is unknown. Here, we show that galanin expression in the ventrobasal thalamus of infant mice coincides with whisker map development and modulates subcortical circuit wiring. Time-resolved neuroanatomy and single-nucleus RNA-seq identified complementary galanin (Gal) and galanin receptor 1 (Galr1) expression in the ventrobasal thalamus and the principal sensory nucleus of the trigeminal nerve (Pr5), respectively. Somatodendritic galanin release from the ventrobasal thalamus was time-locked to the first postnatal week, when Gal1R+ Pr5 afferents form glutamatergic (Slc17a6 +) synapses for the topographical whisker map to emerge. RNAi-mediated silencing of galanin expression disrupted glutamatergic synaptogenesis, which manifested as impaired whisker-dependent exploratory behaviors in infant mice, with behavioral abnormalities enduring into adulthood. Pharmacological probing of receptor selectivity in vivo corroborated that target recognition and synaptogenesis in the thalamus, at least in part, are reliant on agonist-induced Gal1R activation in inbound excitatory axons. Overall, we suggest a neuropeptide-dependent developmental mechanism to contribute to the topographical specification of a fundamental sensory neurocircuit in mice.

ABSTRACTSupplementation of one-carbon (1C) metabolism micronutrients, which include B-vitamins and methionine, is essential for the healthy growth and development of Atlantic salmon (Salmo salar). However, the recent shift towards non-fish meal diets in salmon aquaculture has led to the need for reassessments of recommended micronutrient levels. Despite the importance of 1C metabolism in growth performance and various cellular regulations, the molecular mechanisms affected by these dietary alterations are less understood. To investigate the molecular effect of 1C nutrients, we analysed gene expression and DNA methylation using two types of omics data: RNA sequencing (RNA-seq) and reduced-representation bisulphite sequencing (RRBS). We collected liver samples at the end of a feeding trial that lasted 220 days through the smoltification stage, where fish were fed three different levels of four key 1C nutrients: methionine, vitamin B6, B9, and B12. Our results indicate that the dosage of 1C nutrients significantly impacts genetic and epigenetic regulations in the liver of Atlantic salmon, particularly in biological pathways related to protein synthesis. The interplay between DNA methylation and gene expression in these pathways may play an important role in the mechanisms underlying growth performance affected by 1C metabolism.

Abstract SMNDC1 is a Tudor domain protein that recognizes di-methylated arginines and controls gene expression as an essential splicing factor. Here, we study the specific contributions of the SMNDC1 Tudor domain to protein-protein interactions, subcellular localization, and molecular function. To perturb the protein function in cells, we develop small molecule inhibitors targeting the dimethylarginine binding pocket of the SMNDC1 Tudor domain. We find that SMNDC1 localizes to phase-separated membraneless organelles that partially overlap with nuclear speckles. This condensation behavior is driven by the unstructured C-terminal region of SMNDC1, depends on RNA interaction and can be recapitulated in vitro. Inhibitors of the protein’s Tudor domain drastically alter protein-protein interactions and subcellular localization, causing splicing changes for SMNDC1-dependent genes. These compounds will enable further pharmacological studies on the role of SMNDC1 in the regulation of nuclear condensates, gene regulation and cell identity.

Abstract Metastasis is the major cause of cancer-related deaths. Neuroblastoma (NB), a childhood tumor has been molecularly defined at the primary cancer site, however, the bone marrow (BM) as the metastatic niche of NB is poorly characterized. Here we perform single-cell transcriptomic and epigenomic profiling of BM aspirates from 11 subjects spanning three major NB subtypes and compare these to five age-matched and metastasis-free BM, followed by in-depth single cell analyses of tissue diversity and cell-cell interactions, as well as functional validation. We show that cellular plasticity of NB tumor cells is conserved upon metastasis and tumor cell type composition is NB subtype-dependent. NB cells signal to the BM microenvironment, rewiring via macrophage mgration inhibitory factor and midkine signaling specifically monocytes, which exhibit M1 and M2 features, are marked by activation of pro- and anti-inflammatory programs, and express tumor-promoting factors, reminiscent of tumor-associated macrophages. The interactions and pathways characterized in our study provide the basis for therapeutic approaches that target tumor-to-microenvironment interactions.

Abstract Background Epigenetic alterations are a near-universal feature of human malignancy and have been detected in malignant cells as well as in easily accessible specimens such as blood and urine. These findings offer promising applications in cancer detection, subtyping, and treatment monitoring. However, much of the current evidence is based on findings in retrospective studies and may reflect epigenetic patterns that have already been influenced by the onset of the disease. Methods Studying breast cancer, we established genome-scale DNA methylation profiles of prospectively collected buffy coat samples (n = 702) from a case–control study nested within the EPIC-Heidelberg cohort using reduced representation bisulphite sequencing (RRBS). Results We observed cancer-specific DNA methylation events in buffy coat samples. Increased DNA methylation in genomic regions associated with SURF6 and REXO1/CTB31O20.3 was linked to the length of time to diagnosis in the prospectively collected buffy coat DNA from individuals who subsequently developed breast cancer. Using machine learning methods, we piloted a DNA methylation-based classifier that predicted case–control status in a held-out validation set with 76.5% accuracy, in some cases up to 15 years before clinical diagnosis of the disease. Conclusions Taken together, our findings suggest a model of gradual accumulation of cancer-associated DNA methylation patterns in peripheral blood, which may be detected long before clinical manifestation of cancer. Such changes may provide useful markers for risk stratification and, ultimately, personalized cancer prevention.

Abstract Background Cancer‐associated cachexia (CAC) is a wasting syndrome drastically reducing efficacy of chemotherapy and life expectancy of patients. CAC affects up to 80% of cancer patients, yet the mechanisms underlying the disease are not well understood and no approved disease‐specific medication exists. As a multiorgan disorder, CAC can only be studied on an organismal level. To cover the diverse aetiologies of CAC, researchers rely on the availability of a multifaceted pool of cancer models with varying degrees of cachexia symptoms. So far, no tumour model syngeneic to C57BL/6 mice exists that allows direct comparison between cachexigenic‐ and non‐cachexigenic tumours. Methods MCA207 and CHX207 fibrosarcoma cells were intramuscularly implanted into male or female, 10–11‐week‐old C57BL/6J mice. Tumour tissues were subjected to magnetic resonance imaging, immunohistochemical‐, and transcriptomic analysis. Mice were analysed for tumour growth, body weight and ‐composition, food‐ and water intake, locomotor activity, O2 consumption, CO2 production, circulating blood cells, metabolites, and tumourkines. Mice were sacrificed with same tumour weights in all groups. Adipose tissues were examined using high‐resolution respirometry, lipolysis measurements in vitro and ex vivo, and radioactive tracer studies in vivo. Gene expression was determined in adipose‐ and muscle tissues by quantitative PCR and Western blotting analyses. Muscles and cultured myotubes were analysed histologically and by immunofluorescence microscopy for myofibre cross sectional area and myofibre diameter, respectively. Interleukin‐6 (Il‐6) was deleted from cancer cells using CRISPR/Cas9 mediated gene editing. Results CHX207, but not MCA207‐tumour‐bearing mice exhibited major clinical features of CAC, including systemic inflammation, increased plasma IL‐6 concentrations (190 pg/mL, P ≤ 0.0001), increased energy expenditure (+28%, P ≤ 0.01), adipose tissue loss (−47%, P ≤ 0.0001), skeletal muscle wasting (−18%, P ≤ 0.001), and body weight reduction (−13%, P ≤ 0.01) 13 days after cancer cell inoculation. Adipose tissue loss resulted from reduced lipid uptake and ‐synthesis combined with increased lipolysis but was not associated with elevated beta‐adrenergic signalling or adipose tissue browning. Muscle atrophy was evident by reduced myofibre cross sectional area (−21.8%, P ≤ 0.001), increased catabolic‐ and reduced anabolic signalling. Deletion of IL‐6 from CHX207 cancer cells completely protected CHX207IL6KO‐tumour‐bearing mice from CAC. Conclusions In this study, we present CHX207 fibrosarcoma cells as a novel tool to investigate the mediators and metabolic consequences of CAC in C57BL/6 mice in comparison to non‐cachectic MCA207‐tumour‐bearing mice. IL‐6 represents an essential trigger for CAC development in CHX207‐tumour‐bearing mice.

DNA methylation is involved in regulatory processes throughout the animal kingdom. Here, the authors map DNA methylation in 535 vertebrates and 45 invertebrates, establishing a reference dataset for cross-species analysis and exploring epigenetic variation across vertebrate evolution.

BackgroundEven during long-term combination antiretroviral therapy (cART), people living with HIV (PLHIV) have a dysregulated immune system, characterized by persistent immune activation, accelerated immune ageing and increased risk of non-AIDS comorbidities. A multi-omics approach is applied to a large cohort of PLHIV to understand pathways underlying these dysregulations in order to identify new biomarkers and novel genetically validated therapeutic drugs targets.MethodsThe 2000HIV study is a prospective longitudinal cohort study of PLHIV on cART. In addition, untreated HIV spontaneous controllers were recruited. In-depth multi-omics characterization will be performed, including genomics, epigenomics, transcriptomics, proteomics, metabolomics and metagenomics, functional immunological assays and extensive immunophenotyping. Furthermore, the latent viral reservoir will be assessed through cell associated HIV-1 RNA and DNA, and full-length individual proviral sequencing on a subset. Clinical measurements include an ECG, carotid intima-media thickness and plaque measurement, hepatic steatosis and fibrosis measurement as well as psychological symptoms and recreational drug questionnaires. Additionally, considering the developing pandemic, COVID-19 history and vaccination was recorded. Participants return for a two-year follow-up visit. The 2000HIV study consists of a discovery and validation cohort collected at separate sites to immediately validate any finding in an independent cohort.ResultsOverall, 1895 PLHIV from four sites were included for analysis, 1559 in the discovery and 336 in the validation cohort. The study population was representative of a Western European HIV population, including 288 (15.2%) cis-women, 463 (24.4%) non-whites, and 1360 (71.8%) MSM (Men who have Sex with Men). Extreme phenotypes included 114 spontaneous controllers, 81 rapid progressors and 162 immunological non-responders. According to the Framingham score 321 (16.9%) had a cardiovascular risk of >20% in the next 10 years. COVID-19 infection was documented in 234 (12.3%) participants and 474 (25.0%) individuals had received a COVID-19 vaccine.ConclusionThe 2000HIV study established a cohort of 1895 PLHIV that employs multi-omics to discover new biological pathways and biomarkers to unravel non-AIDS comorbidities, extreme phenotypes and the latent viral reservoir that impact the health of PLHIV. The ultimate goal is to contribute to a more personalized approach to the best standard of care and a potential cure for PLHIV.

The neural crest (NC) is an important multipotent embryonic cell population and its impaired specification leads to various developmental defects, often in an anteroposterior (A-P) axial level-specific manner. The mechanisms underlying the correct A-P regionalisation of human NC cells remain elusive. Recent studies have indicated that trunk NC cells, the presumed precursors of childhood tumour neuroblastoma, are derived from neuromesodermal-potent progenitors of the postcranial body. Here we employ human embryonic stem cell differentiation to define how neuromesodermal progenitor (NMP)-derived NC cells acquire a posterior axial identity. We show that TBXT, a pro-mesodermal transcription factor, mediates early posterior NC/spinal cord regionalisation together with WNT signalling effectors. This occurs by TBXT-driven chromatin remodelling via its binding in key enhancers within HOX gene clusters and other posterior regulator-associated loci. This initial posteriorisation event is succeeded by a second phase of trunk HOX gene control that marks the differentiation of NMPs toward their TBXT-negative NC/spinal cord derivatives and relies predominantly on FGF signalling. Our work reveals a previously unknown role of TBXT in influencing posterior NC fate and points to the existence of temporally discrete, cell type-dependent modes of posterior axial identity control.

Abstract Mutations of calreticulin (CALR) are the second most prevalent driver mutations in essential thrombocythemia and primary myelofibrosis. To identify potential targeted therapies for CALR mutated myeloproliferative neoplasms, we searched for small molecules that selectively inhibit the growth of CALR mutated cells using high-throughput drug screening. We investigated 89 172 compounds using isogenic cell lines carrying CALR mutations and identified synthetic lethality with compounds targeting the ATR-CHK1 pathway. The selective inhibitory effect of these compounds was validated in a co-culture assay of CALR mutated and wild-type cells. Of the tested compounds, CHK1 inhibitors potently depleted CALR mutated cells, allowing wild-type cell dominance in the co-culture over time. Neither CALR deficient cells nor JAK2V617F mutated cells showed hypersensitivity to ATR-CHK1 inhibition, thus suggesting specificity for the oncogenic activation by the mutant CALR. CHK1 inhibitors induced replication stress in CALR mutated cells revealed by elevated pan-nuclear staining for γH2AX and hyperphosphorylation of RPA2. This was accompanied by S-phase cell cycle arrest due to incomplete DNA replication. Transcriptomic and phosphoproteomic analyses revealed a replication stress signature caused by oncogenic CALR, suggesting an intrinsic vulnerability to CHK1 perturbation. This study reveals the ATR-CHK1 pathway as a potential therapeutic target in CALR mutated hematopoietic cells.

In most mammals, imprinted genes contain epigenetic marks that differ in each parental genome and control their parent-of-origin-specific expression. Here, the authors map imprinted genes in mouse preimplantation embryos and find that imprinted gene expression in blastocysts is mainly dependent on Polycomb-mediated H3K27me3-associated gene silencing.

The class I histone deacetylases are essential regulators of cell fate decisions in health and disease. While pan- and class-specific HDAC inhibitors are available, these drugs do not allow a comprehensive understanding of individual HDAC function, or the therapeutic potential of isoform-specific targeting. To systematically compare the impact of individual catalytic functions of HDAC1, HDAC2 and HDAC3, we generated human HAP1 cell lines expressing catalytically inactive HDAC enzymes. Using this genetic toolbox we compare the effect of individual HDAC inhibition with the effects of class I specific inhibitors on cell viability, protein acetylation and gene expression. Individual inactivation of HDAC1 or HDAC2 has only mild effects on cell viability, while HDAC3 inactivation or loss results in DNA damage and apoptosis. Inactivation of HDAC1/HDAC2 led to increased acetylation of components of the COREST co-repressor complex, reduced deacetylase activity associated with this complex and derepression of neuronal genes. HDAC3 controls the acetylation of nuclear hormone receptor associated proteins and the expression of nuclear hormone receptor regulated genes. Acetylation of specific histone acetyltransferases and HDACs is sensitive to inactivation of HDAC1/HDAC2. Over a wide range of assays, we determined that in particular HDAC1 or HDAC2 catalytic inactivation mimics class I specific HDAC inhibitors. Importantly, we further demonstrate that catalytic inactivation of HDAC1 or HDAC2 sensitizes cells to specific cancer drugs. In summary, our systematic study revealed isoform-specific roles of HDAC1/2/3 catalytic functions. We suggest that targeted genetic inactivation of particular isoforms effectively mimics pharmacological HDAC inhibition allowing the identification of relevant HDACs as targets for therapeutic intervention.

Liquid biopsies enable minimally invasive applications for diagnosis and treatment monitoring. Here the authors analyse fragmentation patterns of circulating tumour DNA on multiple levels and develop a bioinformatic tool, LIQUORICE, to accurately detect and classify paediatric cancers with low mutational burden.

Chronic rhinosinusitis with nasal polyps is affecting up to 3% of Western populations. About 10% of patients with nasal polyps also suffer from asthma and intolerance to aspirin, a syndrome called aspirin-exacerbated respiratory disease. Although eosinophilic inflammation is predominant in polyps of both diseases, phenotypic differences in the tissue-derived microenvironment, elucidating disease-specific characteristics, have not yet been identified. We sought to obtain detailed information about phenotypic and transcriptional differences in epithelial and immune cells in polyps of aspirin-tolerant and intolerant patients. Cytokine profiles in nasal secretions and serum of patients suffering from aspirin-exacerbated respiratory disease (n = 10) or chronic rhinosinusitis with nasal polyps (n = 9) were assessed using a multiplex mesoscale discovery assay. After enrichment for immune cell subsets by flow cytometry, we performed transcriptomic profiling by employing single-cell RNA sequencing. Aspirin-intolerant patients displayed significantly elevated IL-5 and CCL17 levels in nasal secretions corresponding to a more pronounced eosinophilic type 2 inflammation. Transcriptomic profiling revealed that epithelial and mast cells not only complement one another in terms of gene expression associated with the 15-lipoxygenase pathway but also show a clear type 2-associated inflammatory phenotype as identified by the upregulation of POSTN, CCL26, and IL13 in patients with aspirin-exacerbated respiratory disease. Interestingly, we also observed cellular stress responses indicated by an increase of MTRNR2L12, MTRNR2L8, and NEAT1 across all immune cell subsets in this disease entity. In conclusion, our findings support the hypothesis that epithelial and mast cells act in concert as potential drivers of the pathogenesis of the aspirin-exacerbated respiratory disease.

Abstract Background Deep learning has emerged as a versatile approach for predicting complex biological phenomena. However, its utility for biological discovery has so far been limited, given that generic deep neural networks provide little insight into the biological mechanisms that underlie a successful prediction. Here we demonstrate deep learning on biological networks, where every node has a molecular equivalent, such as a protein or gene, and every edge has a mechanistic interpretation, such as a regulatory interaction along a signaling pathway. Results With knowledge-primed neural networks (KPNNs), we exploit the ability of deep learning algorithms to assign meaningful weights in multi-layered networks, resulting in a widely applicable approach for interpretable deep learning. We present a learning method that enhances the interpretability of trained KPNNs by stabilizing node weights in the presence of redundancy, enhancing the quantitative interpretability of node weights, and controlling for uneven connectivity in biological networks. We validate KPNNs on simulated data with known ground truth and demonstrate their practical use and utility in five biological applications with single-cell RNA-seq data for cancer and immune cells. Conclusions We introduce KPNNs as a method that combines the predictive power of deep learning with the interpretability of biological networks. While demonstrated here on single-cell sequencing data, this method is broadly relevant to other research areas where prior domain knowledge can be represented as networks.

Abstract Background 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. Results We find that contamination by cell-free RNA can constitute up to 20% of reads in human primary tissue samples, and we show that the ensuing biases can be removed effectively using a novel bioinformatics algorithm. Applying our method to both human and mouse pancreatic islets treated ex vivo, we obtain an accurate and quantitative assessment of cell-specific drug effects on the transcriptome. We observe that FOXO inhibition induces dedifferentiation of both alpha and beta cells, while artemether treatment upregulates insulin and other beta cell marker genes in a subset of alpha cells. In beta cells, dedifferentiation and insulin repression upon artemether treatment occurs predominantly in mouse but not in human samples. Conclusions This new method for quantitative, error-correcting, scRNA-seq data normalization using spike-in reference cells helps clarify complex cell-specific effects of pharmacological perturbations with single-cell resolution and high quantitative accuracy.

Ibrutinib, a Bruton tyrosine kinase inhibitor, provides effective treatment for chronic lymphocytic leukemia (CLL). Here, the authors describe time-dependent molecular changes to malignant cells and to the immune system in patients undergoing ibrutinib therapy, with can be used for therapy monitoring.

Abstract The balance of activating and inhibitory signals from the low affinity Fc gamma receptors modulates immune responses triggered by IgG antibody‐immune complexes. In homeostasis, this leads to antigen clearance, while in autoimmune diseases to unwanted immune response. Besides the activating receptors FcɣRIIa, FcɣRIIIa, and the inhibitory FcɣRIIb receptor, a third activating receptor, FcɣRIIc, was shown to be expressed on several immune cell types, however, only in the presence of a functional FCGR2C‐ORF allele. FcɣRIIc expression is associated with autoimmune diseases such as idiopathic thrombocytopenic purpura, systemic lupus erythematosus or systemic sclerosis. Thus, the determination of the functional FCGR2C gene resulting in protein expression on immune cells becomes highly relevant, particularly in the context of unwanted immune responses through inadvertent FcɣRIIc activation by molecules targeting stimulation of the inhibitory receptor FcɣRIIb, currently pursued by several pharmaceutical companies. The high degree of homology within the FCGR2/3 gene cluster complicates development of an accurate method for identification of FcɣRIIc expression. Here we describe a comprehensive approach to characterize genetic status of the FCGR2C gene locus consisting of cDNA sequencing, SNaPshot genotyping and low‐coverage next‐generation sequencing. This might enable Mendelian randomization hypothesis testing across autoimmune diseases to personalize therapies and enhance treatment outcomes.

Objective: The loss of forkhead box protein O1 (FoxO1) signaling in response to metabolic stress contributes to the etiology of type II diabetes, causing the dedifferentiation of pancreatic beta cells to a cell type reminiscent of endocrine progenitors. Lack of methods to easily model this process in vitro, however, have hindered progress into the identification of key downstream targets and potential inhibitors. We therefore aimed to establish such an in vitro cellular dedifferentiation model and apply it to identify novel agents involved in the maintenance of beta-cell identity. Methods: The murine beta-cell line, Min6, was used for primary experiments and high-content screening. Screens encompassed a library of small-molecule drugs representing the chemical and target space of all FDA-approved small molecules with an automated immunofluorescence readout. Validation experiments were performed in a murine alpha-cell line as well as in primary murine and human diabetic islets. Developmental effects were studied in zebrafish and C. elegans models, while diabetic db/db mouse models were used to elucidate global glucose metabolism outcomes. Results: We show that short-term pharmacological FoxO1 inhibition can model beta-cell dedifferentiation by downregulating beta-cell-specific transcription factors, resulting in the aberrant expression of progenitor genes and the alpha-cell marker glucagon. From a high-content screen, we identified loperamide as a small molecule that can prevent FoxO inhibitor-induced glucagon expression and further stimulate insulin protein processing and secretion by altering calcium levels, intracellular pH, and FoxO1 localization. Conclusions: Our study provides novel models, molecular targets, and drug candidates for studying and preventing beta-cell dedifferentiation.

Objectives We explore the importance of SARS-CoV-2 sentinel surveillance testing in primary care during a regional COVID-19 outbreak in Austria.Design Prospective cohort study.Setting A single sentinel practice serving 22 829 people in the ski-resort of Schladming-Dachstein.Participants All 73 patients presenting with mild-to-moderate flu-like symptoms between 24 February and 03 April, 2020.Intervention Nasopharyngeal sampling to detect SARS-CoV-2 using real-time reverse transcriptase-quantitative PCR (RT-qPCR).Outcome measures We compared RT-qPCR at presentation with confirmed antibody status. We split the outbreak in two parts, by halving the period from the first to the last case, to characterise three cohorts of patients with confirmed infection: early acute (RT-qPCR reactive) in the first half; and late acute (reactive) and late convalescent (non-reactive) in the second half. For each cohort, we report the number of cases detected, the accuracy of RT-qPCR, the duration and variety of symptoms, and the number of viral clades present.Results Twenty-two patients were diagnosed with COVID-19 (eight early acute, seven late acute and seven late convalescent), 44 patients tested SARS-CoV-2 negative and 7 were excluded. The sensitivity of RT-qPCR was 100% among all acute cases, dropping to 68.1% when including convalescent. Test specificity was 100%. Mean duration of symptoms for each group were 2 days (range 1–4) among early acute, 4.4 days (1–7) among late acute and 8 days (2–12) among late convalescent. Confirmed infection was associated with loss of taste. Acute infection was associated with loss of taste, nausea/vomiting, breathlessness, sore throat and myalgia; but not anosmia, fever or cough. Transmission clusters of three viral clades (G, GR and L) were identified.Conclusions RT-qPCR testing in primary care can rapidly and accurately detect SARS-CoV-2 among people with flu-like illness in a heterogeneous viral outbreak. Targeted testing in primary care can support national sentinel surveillance of COVID-19.

Deregulation of transcription factors is frequently observed in cancer. Here, the authors develop a blood-based approach to map transcription factor binding, revealing both patient-specific and tumour-specific aberrations

The regulation of tumor inflammation is incompletely understood and the role of B cells is unclear. Here, the authors show that a specific subtype of B cells is induced in melanoma and required to recruit T lymphocytes and elicit inflammation.

Abstract DNA methylation is a widely investigated epigenetic mark with important roles in development and disease. High-throughput assays enable genome-scale DNA methylation analysis in large numbers of samples. Here, we describe a new version of our RnBeads software - an R/Bioconductor package that implements start-to-finish analysis workflows for Infinium microarrays and various types of bisulfite sequencing. RnBeads 2.0 (https://rnbeads.org/) provides additional data types and analysis methods, new functionality for interpreting DNA methylation differences, improved usability with a novel graphical user interface, and better use of computational resources. We demonstrate RnBeads 2.0 in four re-runnable use cases focusing on cell differentiation and cancer.

During the past few years, our understanding of molecular mechanisms and cellular interactions relevant to malignant blood cell disorders has improved substantially. New insights include a detailed knowledge about disease-initiating exogenous factors, endogenous (genetic, somatic, epigenetic) elicitors or facilitators of disease evolution, and drug actions and interactions that underlie efficacy and adverse event profiles in defined cohorts of patients. As a result, precision medicine and personalized medicine are rapidly growing new disciplines that support the clinician in making the correct diagnosis, in predicting outcomes, and in optimally selecting patients for interventional therapies. In addition, precision medicine tools are greatly facilitating the development of new drugs, therapeutic approaches, and new multiparametric prognostic scoring models. However, although the emerging roles of precision medicine and personalized medicine in hematology and oncology are clearly visible, several questions remain. For example, it remains unknown how precision medicine tools can be implemented in healthcare systems and whether all possible approaches are also affordable. In addition, there is a need to define terminologies and to relate these to specific and context-related tools and strategies in basic and applied science. To discuss these issues, a working conference was organized in September 2019. The outcomes of this conference are summarized herein and include a proposal for definitions, terminologies, and applications of precision and personalized medicine concepts and tools in hematologic neoplasms. We also provide proposals aimed at reducing costs, thereby making these applications affordable in daily practice.

The BTB zinc finger transcription factor MAZR (also known as PATZ1) controls, partially in synergy with the transcription factor Runx3, the development of CD8 lineage T cells. Here we explored the role of MAZR as well as combined activities of MAZR/Runx3 during cytotoxic T lymphocyte (CTL) and memory CD8+ T cell differentiation. In contrast to the essential role of Runx3 for CTL effector function, the deletion of MAZR had a mild effect on the generation of CTLs in vitro. However, a transcriptome analysis demonstrated that the combined deletion of MAZR and Runx3 resulted in much more widespread downregulation of CTL signature genes compared to single Runx3 deletion, indicating that MAZR partially compensates for loss of Runx3 in CTLs. Moreover, in line with the findings made in vitro, the analysis of CTL responses to LCMV infection revealed that MAZR and Runx3 cooperatively regulate the expression of CD8α, Granzyme B and perforin in vivo. Interestingly, while memory T cell differentiation is severely impaired in Runx3-deficient mice, the deletion of MAZR leads to an enlargement of the long-lived memory subset and also partially restored the differentiation defect caused by loss of Runx3. This indicates distinct functions of MAZR and Runx3 in the generation of memory T cell subsets, which is in contrast to their cooperative roles in CTLs. Together, our study demonstrates complex interplay between MAZR and Runx3 during CTL and memory T cell differentiation, and provides further insight into the molecular mechanisms underlying the establishment of CTL and memory T cell pools.

Atopic dermatitis (AD) typically starts in infancy or early childhood, showing spontaneous remission in a subset of patients, while others develop lifelong disease. Despite an increased understanding of AD, factors guiding its natural course are only insufficiently elucidated. We thus performed suction blistering in skin of adult patients with stable, spontaneous remission from previous moderate-to-severe AD during childhood. Samples were compared to healthy controls without personal or familial history of atopy, and to chronic, active AD lesions. Skin cells and tissue fluid obtained were used for single-cell RNA sequencing and proteomic multiplex assays, respectively. We found overall cell composition and proteomic profiles of spontaneously healed AD to be comparable to healthy control skin, without upregulation of typical AD activity markers (e.g., IL13, S100As, and KRT16). Among all cell types in spontaneously healed AD, melanocytes harbored the largest numbers of differentially expressed genes in comparison to healthy controls, with upregulation of potentially anti-inflammatory markers such as PLA2G7. Conventional T-cells also showed increases in regulatory markers, and a general skewing toward a more Th1-like phenotype. By contrast, gene expression of regulatory T-cells and keratinocytes was essentially indistinguishable from healthy skin. Melanocytes and conventional T-cells might thus contribute a specific regulatory milieu in spontaneously healed AD skin.

The differentiation of naïve CD4+ T cells into T helper (Th) subsets is key for a functional immune response and has to be tightly controlled by transcriptional and epigenetic processes. However, the function of cofactors that connect gene-specific transcription factors with repressive chromatin-modifying enzymes in Th cells is yet unknown. Here we demonstrate an essential role for nuclear receptor corepressor 1 (NCOR1) in regulating naïve CD4+ T cell and Th1/Th17 effector transcriptomes. Moreover, NCOR1 binds to a conserved cis-regulatory element within the Ifng locus and controls the extent of IFNγ expression in Th1 cells. Further, NCOR1 controls the survival of activated CD4+ T cells and Th1 cells in vitro, while Th17 cell survival was not affected in the absence of NCOR1. In vivo, effector functions were compromised since adoptive transfer of NCOR1-deficient CD4+ T cells resulted in attenuated colitis due to lower frequencies of IFNγ+ and IFNγ+IL-17A+ Th cells and overall reduced CD4+ T cell numbers. Collectively, our data demonstrate that the coregulator NCOR1 shapes transcriptional landscapes in CD4+ T cells and controls Th1/Th17 effector functions.

Abstract The ATP-binding cassette transporter TAPL translocates polypeptides from the cytosol into the lysosomal lumen. TAPL can be divided into two functional units: coreTAPL, active in ATP-dependent peptide translocation, and the N-terminal membrane spanning domain, TMD0, responsible for cellular localization and interaction with the lysosomal associated membrane proteins LAMP-1 and LAMP-2. Although the structure and function of ABC transporters were intensively analyzed in the past, the knowledge about accessory membrane embedded domains is limited. Therefore, we expressed the TMD0 of TAPL via a cell-free expression system and confirmed its correct folding by NMR and interaction studies. In cell as well as cell-free expressed TMD0 forms oligomers, which were assigned as dimers by PELDOR spectroscopy and static light scattering. By NMR spectroscopy of uniformly and selectively isotope labeled TMD0 we performed a complete backbone and partial side chain assignment. Accordingly, TMD0 has a four transmembrane helix topology with a short helical segment in a lysosomal loop. The topology of TMD0 was confirmed by paramagnetic relaxation enhancement with paramagnetic stearic acid as well as by nuclear Overhauser effects with c6-DHPC and cross-peaks with water.

Summary: Forkhead box protein P3+ (FOXP3+) regulatory T cells (Treg cells) play a key role in maintaining tolerance and immune homeostasis. Here, we report that a T cell-specific deletion of the transcription factor MAZR (also known as PATZ1) leads to an increased frequency of Treg cells, while enforced MAZR expression impairs Treg cell differentiation. Further, MAZR expression levels are progressively downregulated during thymic Treg cell development and during in-vitro-induced human Treg cell differentiation, suggesting that MAZR protein levels are critical for controlling Treg cell development. However, MAZR-deficient Treg cells show only minor transcriptional changes ex vivo, indicating that MAZR is not essential for establishing the transcriptional program of peripheral Treg cells. Finally, the loss of MAZR reduces the clinical score in dextran-sodium sulfate (DSS)-induced colitis, suggesting that MAZR activity in T cells controls the extent of intestinal inflammation. Together, these data indicate that MAZR is part of a Treg cell-intrinsic transcriptional network that modulates Treg cell development. : FOXP3+ Treg cells are essential for maintaining tolerance and immune homeostasis. Andersen et al. reveal that MAZR is an important factor in regulating the delicate balance of Treg cell generation and report that MAZR expression levels play a key role in controlling Treg cell development and differentiation. Keywords: MAZR, PATZ1, FOXP3, regulatory T cells, Treg, DSS-induced colitis

Abstract Background Mammalian species exhibit a wide range of lifespans. To date, a robust and dynamic molecular readout of these lifespan differences has not yet been identified. Recent studies have established the existence of ageing-associated differentially methylated positions (aDMPs) in human and mouse. These are CpG sites at which DNA methylation dynamics show significant correlations with age. We hypothesise that aDMPs are pan-mammalian and are a dynamic molecular readout of lifespan variation among different mammalian species. Results A large-scale integrated analysis of aDMPs in six different mammals reveals a strong negative relationship between rate of change of methylation levels at aDMPs and lifespan. This relationship also holds when comparing two different dog breeds with known differences in lifespans. In an ageing cohort of aneuploid mice carrying a complete copy of human chromosome 21, aDMPs accumulate far more rapidly than is seen in human tissues, revealing that DNA methylation at aDMP sites is largely shaped by the nuclear trans-environment and represents a robust molecular readout of the ageing cellular milieu. Conclusions Overall, we define the first dynamic molecular readout of lifespan differences among mammalian species and propose that aDMPs will be an invaluable molecular tool for future evolutionary and mechanistic studies aimed at understanding the biological factors that determine lifespan in mammals.

Cancer evolution is a step-wise non-linear process that may start early in life or later in adulthood, and includes pre-malignant (indolent) and malignant phases. Early somatic changes may not be detectable or are found by chance in apparently healthy individuals. The same lesions may be detected in pre-malignant clonal conditions. In some patients, these lesions may never become relevant clinically whereas in others, they act together with additional pro-oncogenic hits and thereby contribute to the formation of an overt malignancy. Although some pre-malignant stages of a malignancy have been characterized, no global system to define and to classify these conditions is available. To discuss open issues related to pre-malignant phases of neoplastic disorders, a working conference was organized in Vienna in August 2015. The outcomes of this conference are summarized herein and include a basic proposal for a nomenclature and classification of pre-malignant conditions. This proposal should assist in the communication among patients, physicians and scientists, which is critical as genome-sequencing will soon be offered widely for early cancer-detection.

Fanconi anemia is a complex disease affecting multiple DNA repair proteins that resolve DNA crosslinks which can block vital processes. Here the authors use parallel genome-wide screens that identify the BLM helicase complex as a suppressor of Fanconi anemia phenotypes.

Summary: Regulatory B cells (Bregs) play a critical role in the control of autoimmunity and inflammation. IL-10 production is the hallmark for the identification of Bregs. However, the molecular determinants that regulate the transcription of IL-10 and control the Breg developmental program remain unknown. Here, we demonstrate that aryl hydrocarbon receptor (AhR) regulates the differentiation and function of IL-10-producing CD19+CD21hiCD24hiBregs and limits their differentiation into B cells that contribute to inflammation. Chromatin profiling and transcriptome analyses show that loss of AhR in B cells reduces expression of IL-10 by skewing the differentiation of CD19+CD21hiCD24hiB cells into a pro-inflammatory program, under Breg-inducing conditions. B cell AhR-deficient mice develop exacerbated arthritis, show significant reductions in IL-10-producing Bregs and regulatory T cells, and show an increase in T helper (Th) 1 and Th17 cells compared with B cell AhR-sufficient mice. Thus, we identify AhR as a relevant contributor to the transcriptional regulation of Breg differentiation. : The transcriptional control of interleukin-10 (IL-10) in regulatory B cells (Bregs) remains undefined. Piper et al. identify the aryl hydrocarbon receptor (AhR) as an important transcription factor involved in Breg differentiation and show a direct role of AhR in the regulation of IL-10 transcription.

Abstract Reverse genetic screens have driven gene annotation and target discovery in model organisms. However, many disease‐relevant genotypes and phenotypes cannot be studied in lower organisms. It is therefore essential to overcome technical hurdles associated with large‐scale reverse genetics in human cells. Here, we establish a reverse genetic approach based on highly robust and sensitive multiplexed RNA sequencing of mutant human cells. We conduct 10 parallel screens using a collection of engineered haploid isogenic cell lines with knockouts covering tyrosine kinases and identify known and unexpected effects on signaling pathways. Our study provides proof of concept for a scalable approach to link genotype to phenotype in human cells, which has broad applications. In particular, it clears the way for systematic phenotyping of still poorly characterized human genes and for systematic study of uncharacterized genomic features associated with human disease.

Chronic lymphocytic leukemia (CLL) is characterized by substantial clinical heterogeneity. Here, the authors report the genome-wide chromatin accessibility maps for 88 CLL samples from 55 patients using ATAC-seq, and 10 matched RNA-seq datasets, providing a resource for studying epigenome deregulation in CLL.

The cellular response to replication stress requires the DNA-damage-responsive kinase ATM and its cofactor ATMIN; however, the roles of this signaling pathway following replication stress are unclear. To identify the functions of ATM and ATMIN in response to replication stress, we utilized both transcriptomics and quantitative mass-spectrometry-based phosphoproteomics. We found that replication stress induced by aphidicolin triggered widespread changes in both gene expression and protein phosphorylation patterns. These changes gave rise to distinct early and late replication stress responses. Furthermore, our analysis revealed previously unknown targets of ATM and ATMIN downstream of replication stress. We demonstrate ATMIN-dependent phosphorylation of H2AX and of CRMP2, a protein previously implicated in Alzheimer’s disease but not in the DNA damage response. Overall, our dataset provides a comprehensive resource for discovering the cellular responses to replication stress and, potentially, associated pathologies.

Gestützt auf 16 Empfehlungen sollten Anstrengungen gemacht werden, um folgendes Ziel zu erreichen: Bis 2025 ist die gesamte wissenschaftliche Publikationstätigkeit in Österreich auf Open Access umgestellt. Das bedeutet, dass alle wissenschaftlichen Publikationen, die aus Unterstützungen mit öffentlichen Mitteln hervorgegangen sind, ohne Zeitverzögerung und in der finalen Version im Internet frei zugänglich sind (Gold Open Access). Die notwendigen Mittel werden den AutorInnen zur Verfügung gestellt oder die Kosten der Publikationsorgane werden direkt von den Wissenschaftsorganisationen getragen.Gestützt auf 16 Empfehlungen sollten Anstrengungen gemacht werden, um folgendes Ziel zu erreichen: Bis 2025 ist die gesamte wissenschaftliche Publikationstätigkeit in Österreich auf Open Access umgestellt. Das bedeutet, dass alle wissenschaftlichen Publikationen, die aus Unterstützungen mit öffentlichen Mitteln hervorgegangen sind, ohne Zeitverzögerung und in der finalen Version im Internet frei zugänglich sind (Gold Open Access). Die notwendigen Mittel werden den AutorInnen zur Verfügung gestellt oder die Kosten der Publikationsorgane werden direkt von den Wissenschaftsorganisationen getragen.

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.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Transcription factor fusion proteins can transform cells by inducing global changes of the transcriptome, often creating a state of oncogene addiction. Here, we investigate the role of epigenetic mechanisms in this process, focusing on Ewing sarcoma cells that are dependent on the EWS-FLI1 fusion protein. We established reference epigenome maps comprising DNA methylation, seven histone marks, open chromatin states, and RNA levels, and we analyzed the epigenome dynamics upon downregulation of the driving oncogene. Reduced EWS-FLI1 expression led to widespread epigenetic changes in promoters, enhancers, and super-enhancers, and we identified histone H3K27 acetylation as the most strongly affected mark. Clustering of epigenetic promoter signatures defined classes of EWS-FLI1-regulated genes that responded differently to low-dose treatment with histone deacetylase inhibitors. Furthermore, we observed strong and opposing enrichment patterns for E2F and AP-1 among EWS-FLI1-correlated and anticorrelated genes. Our data describe extensive genome-wide rewiring of epigenetic cell states driven by an oncogenic fusion protein.

Gestützt auf 16 Empfehlungen sollten Anstrengungen gemacht werden, um folgendes Ziel zu erreichen: Bis 2025 ist die gesamte wissenschaftliche Publikationstätigkeit in Österreich auf Open Access umgestellt. Das bedeutet, dass alle wissenschaftlichen Publikationen, die aus Unterstützungen mit öffentlichen Mitteln hervorgegangen sind, ohne Zeitverzögerung und in der finalen Version im Internet frei zugänglich sind (Gold Open Access). Die notwendigen Mittel werden den AutorInnen zur Verfügung gestellt oder die Kosten der Publikationsorgane werden direkt von den Wissenschaftsorganisationen getragen.Gestützt auf 16 Empfehlungen sollten Anstrengungen gemacht werden, um folgendes Ziel zu erreichen: Bis 2025 ist die gesamte wissenschaftliche Publikationstätigkeit in Österreich auf Open Access umgestellt. Das bedeutet, dass alle wissenschaftlichen Publikationen, die aus Unterstützungen mit öffentlichen Mitteln hervorgegangen sind, ohne Zeitverzögerung und in der finalen Version im Internet frei zugänglich sind (Gold Open Access). Die notwendigen Mittel werden den AutorInnen zur Verfügung gestellt oder die Kosten der Publikationsorgane werden direkt von den Wissenschaftsorganisationen getragen.

To determine the dynamics of allelic-specific expression during mouse development, we analyzed RNA-seq data from 23 F1 tissues from different developmental stages, including 19 female tissues allowing X chromosome inactivation (XCI) escapers to also be detected. We demonstrate that allelic expression arising from genetic or epigenetic differences is highly tissue-specific. We find that tissue-specific strain-biased gene expression may be regulated by tissue-specific enhancers or by post-transcriptional differences in stability between the alleles. We also find that escape from X-inactivation is tissue-specific, with leg muscle showing an unexpectedly high rate of XCI escapers. By surveying a range of tissues during development, and performing extensive validation, we are able to provide a high confidence list of mouse imprinted genes including 18 novel genes. This shows that cluster size varies dynamically during development and can be substantially larger than previously thought, with the Igf2r cluster extending over 10 Mb in placenta.

DNA methylation is a central epigenetic modification that is established by de novo DNA methyltransferases. The mechanisms underlying the generation of genomic methylation patterns are still poorly understood. Using mass spectrometry and a phosphospecific Dnmt3a antibody, we demonstrate that CK2 phosphorylates endogenous Dnmt3a at two key residues located near its PWWP domain, thereby downregulating the ability of Dnmt3a to methylate DNA. Genome-wide DNA methylation analysis shows that CK2 primarily modulates CpG methylation of several repeats, most notably of Alu SINEs. This modulation can be directly attributed to CK2-mediated phosphorylation of Dnmt3a. We also find that CK2-mediated phosphorylation is required for localization of Dnmt3a to heterochromatin. By revealing phosphorylation as a mode of regulation of de novo DNA methyltransferase function and by uncovering a mechanism for the regulation of methylation at repetitive elements, our results shed light on the origin of DNA methylation patterns.

The conversion of somatic cells into pluripotent stem cells via overexpression of reprogramming factors involves epigenetic remodeling. DNA methylation at a significant proportion of CpG sites in induced pluripotent stem cells (iPSCs) differs from that of embryonic stem cells (ESCs). Whether different sets of reprogramming factors influence the type and extent of aberrant DNA methylation in iPSCs differently remains unknown. In order to help resolve this critical question, we generated human iPSCs from a common fibroblast cell source using either the Yamanaka factors (OCT4, SOX2, KLF4 and cMYC) or the Thomson factors (OCT4, SOX2, NANOG and LIN28), and determined their genome-wide DNA methylation profiles. In addition to shared DNA methylation aberrations present in all our iPSCs, we identified Yamanaka-iPSC (Y-iPSC)-specific and Thomson-iPSC (T-iPSC)-specific recurrent aberrations. Strikingly, not only were the genomic locations of the aberrations different but also their types: reprogramming with Yamanaka factors mainly resulted in failure to demethylate CpGs, whereas reprogramming with Thomson factors mainly resulted in failure to methylate CpGs. Differences in the level of transcripts encoding DNMT3b and TET3 between Y-iPSCs and T-iPSCs may contribute partially to the distinct types of aberrations. Finally, de novo aberrantly methylated genes in Y-iPSCs were enriched for NANOG targets that are also aberrantly methylated in some cancers. Our study thus reveals that the choice of reprogramming factors influences the amount, location, and class of DNA methylation aberrations in iPSCs. These findings may provide clues into how to produce human iPSCs with fewer DNA methylation abnormalities.

Exome sequencing of primary tumors identifies complex somatic mutation patterns. Assignment of relevance of individual somatic mutations is difficult and poses the next challenge for interpretation of next generation sequencing data. Here we present an approach how exome sequencing in combination with SNP microarray data may identify targets of chromosomal aberrations in myeloid malignancies. The rationale of this approach is that hotspots of chromosomal aberrations might also harbor point mutations in the target genes of deletions, gains or uniparental disomies (UPDs). Chromosome 11 is a frequent target of lesions in myeloid malignancies. Therefore, we studied chromosome 11 in a total of 813 samples from 773 individual patients with different myeloid malignancies by SNP microarrays and complemented the data with exome sequencing in selected cases exhibiting chromosome 11 defects. We found gains, losses and UPDs of chromosome 11 in 52 of the 813 samples (6.4%). Chromosome 11q UPDs frequently associated with mutations of CBL. In one patient the 11qUPD amplified somatic mutations in both CBL and the DNA repair gene DDB1. A duplication within MLL exon 3 was detected in another patient with 11qUPD. We identified several common deleted regions (CDR) on chromosome 11. One of the CDRs associated with de novo acute myeloid leukemia (P=0.013). One patient with a deletion at the LMO2 locus harbored an additional point mutation on the other allele indicating that LMO2 might be a tumor suppressor frequently targeted by 11p deletions. Our chromosome-centered analysis indicates that chromosome 11 contains a number of tumor suppressor genes and that the role of this chromosome in myeloid malignancies is more complex than previously recognized.

H1 linker histones facilitate higher-order chromatin folding and are essential for mammalian development. To achieve high-resolution mapping of H1 variants H1d and H1c in embryonic stem cells (ESCs), we have established a knock-in system and shown that the N-terminally tagged H1 proteins are functionally interchangeable to their endogenous counterparts in vivo. H1d and H1c are depleted from GC- and gene-rich regions and active promoters, inversely correlated with H3K4me3, but positively correlated with H3K9me3 and associated with characteristic sequence features. Surprisingly, both H1d and H1c are significantly enriched at major satellites, which display increased nucleosome spacing compared with bulk chromatin. While also depleted at active promoters and enriched at major satellites, overexpressed H1(0) displays differential binding patterns in specific repetitive sequences compared with H1d and H1c. Depletion of H1c, H1d, and H1e causes pericentric chromocenter clustering and de-repression of major satellites. These results integrate the localization of an understudied type of chromatin proteins, namely the H1 variants, into the epigenome map of mouse ESCs, and we identify significant changes at pericentric heterochromatin upon depletion of this epigenetic mark.

DNA methylation plays an important role in development and disease. The primary sites of DNA methylation in vertebrates are cytosines in the CpG dinucleotide context, which account for roughly three quarters of the total DNA methylation content in human and mouse cells. While the genomic distribution, inter-individual stability, and functional role of CpG methylation are reasonably well understood, little is known about DNA methylation targeting CpA, CpT, and CpC (non-CpG) dinucleotides. Here we report a comprehensive analysis of non-CpG methylation in 76 genome-scale DNA methylation maps across pluripotent and differentiated human cell types. We confirm non-CpG methylation to be predominantly present in pluripotent cell types and observe a decrease upon differentiation and near complete absence in various somatic cell types. Although no function has been assigned to it in pluripotency, our data highlight that non-CpG methylation patterns reappear upon iPS cell reprogramming. Intriguingly, the patterns are highly variable and show little conservation between different pluripotent cell lines. We find a strong correlation of non-CpG methylation and DNMT3 expression levels while showing statistical independence of non-CpG methylation from pluripotency associated gene expression. In line with these findings, we show that knockdown of DNMTA and DNMT3B in hESCs results in a global reduction of non-CpG methylation. Finally, non-CpG methylation appears to be spatially correlated with CpG methylation. In summary these results contribute further to our understanding of cytosine methylation patterns in human cells using a large representative sample set.