Search

Glucocorticoids, which have long served as fundamental therapeutics for diverse inflammatory conditions, are still widely used, despite associated side effects limiting their long-term use. Among their key mediators is glucocorticoid-induced leucine zipper (GILZ), recognized for its anti-inflammatory and immunosuppressive properties. Here, we explore the immunomodulatory effects of GILZ in macrophages through transcriptomic analysis and functional assays. Bulk RNA sequencing of GILZ knockout and GILZ-overexpressing macrophages revealed significant alterations in gene expression profiles, particularly impacting pathways associated with the inflammatory response, phagocytosis, cell death, mitochondrial function, and extracellular structure organization activity. GILZ-overexpression enhances phagocytic and antibacterial activity against Salmonella typhimurium and Escherichia coli, potentially mediated by increased nitric oxide production. In addition, GILZ protects macrophages from pyroptotic cell death, as indicated by a reduced production of reactive oxygen species (ROS) in GILZ transgenic macrophages. In contrast, GILZ KO macrophages produced more ROS, suggesting a regulatory role of GILZ in ROS-dependent pathways. Additionally, GILZ overexpression leads to decreased mitochondrial respiration and heightened matrix metalloproteinase activity, suggesting its involvement in tissue remodeling processes. These findings underscore the multifaceted role of GILZ in modulating macrophage functions and its potential as a therapeutic target for inflammatory disorders, offering insights into the development of novel therapeutic strategies aimed at optimizing the benefits of glucocorticoid therapy while minimizing adverse effects.

Objectives In SLE, deregulation of haematopoiesis is characterised by inflammatory priming and myeloid skewing of haematopoietic stem and progenitor cells (HSPCs). We sought to investigate the role of extramedullary haematopoiesis (EMH) as a key player for tissue injury in systemic autoimmune disorders.Methods Transcriptomic analysis of bone marrow (BM)-derived HSPCs from patients with SLE and NZBW/F1 lupus-prone mice was performed in combination with DNA methylation profile. Trained immunity (TI) was induced through β-glucan administration to the NZBW/F1 lupus-prone model. Disease activity was assessed through lupus nephritis (LN) histological grading. Colony-forming unit assay and adoptive cell transfer were used to assess HSPCs functionalities.Results Transcriptomic analysis shows that splenic HSPCs carry a higher inflammatory potential compared with their BM counterparts. Further induction of TI, through β-glucan administration, exacerbates splenic EMH, accentuates myeloid skewing and worsens LN. Methylomic analysis of BM-derived HSPCs demonstrates myeloid skewing which is in part driven by epigenetic tinkering. Importantly, transcriptomic analysis of human SLE BM-derived HSPCs demonstrates similar findings to those observed in diseased mice.Conclusions These data support a key role of granulocytes derived from primed HSPCs both at medullary and extramedullary sites in the pathogenesis of LN. EMH and TI contribute to SLE by sustaining the systemic inflammatory response and increasing the risk for flare.

Abstract Fatty liver disease or the accumulation of fat in the liver, has been reported to affect the global population. This comes with an increased risk for the development of fibrosis, cirrhosis, and hepatocellular carcinoma. Yet, little is known about the effects of a diet containing high fat and alcohol towards epigenetic aging, with respect to changes in transcriptional and epigenomic profiles. In this study, we took up a multi-omics approach and integrated gene expression, methylation signals, and chromatin signals to study the epigenomic effects of a high-fat and alcohol-containing diet on mouse hepatocytes. We identified four relevant gene network clusters that were associated with relevant pathways that promote steatosis. Using a machine learning approach, we predict specific transcription factors that might be responsible to modulate the functionally relevant clusters. Finally, we discover four additional CpG loci and validate aging-related differential CpG methylation. Differential CpG methylation linked to aging showed minimal overlap with altered methylation in steatosis.

Background Stroke is a leading global cause of human death and disability, with advanced aging associated with elevated incidences of stroke. Despite high mortality and morbidity of stroke, the mechanisms leading to blood‐brain barrier dysfunction and development of stroke with age are poorly understood. In the vasculature of brain, endothelial cells (ECs) constitute the core component of the blood‐brain barrier and provide a physical barrier composed of tight junctions, adherens junctions, and basement membrane. Methods and Results We show, in mice, the incidents of intracerebral bleeding increases with age. After isolating an enriched population of cerebral ECs from murine brains at 2, 6, 12, 18, and 24 months, we studied age‐associated changes in gene expression. The study reveals age‐dependent dysregulation of 1388 genes, including many involved in the maintenance of the blood‐brain barrier and vascular integrity. We also investigated age‐dependent changes on the levels of CpG methylation and accessible chromatin in cerebral ECs. Our study reveals correlations between age‐dependent changes in chromatin structure and gene expression, whereas the dynamics of DNA methylation changes are different. Conclusions We find significant age‐dependent downregulation of the Aplnr gene along with age‐dependent reduction in chromatin accessibility of promoter region of the Aplnr gene in cerebral ECs. Aplnr is associated with positive regulation of vasodilation and is implicated in vascular health. Altogether, our data suggest a potential role of the apelinergic axis involving the ligand apelin and its receptor to be critical in maintenance of the blood‐brain barrier and vascular integrity.

Gonadotropes in the pituitary secrete follicle-stimulating hormone and luteinizing hormone to control gonadal function and fertility, but whether they exert actions on extra-gonadal organs is not fully understood. Here the authors report that gonadotropes regulate liver steatosis independent of the ovaries in mice.

Prenatal exposure to endocrine disrupting chemicals can interfere with development, and has been associated with social-cognitive functioning and adverse health outcomes later in life. Exposure-associated changes of DNA methylation (DNAm) patterns have been suggested as a possible mediator of this relationship. This study investigated whether prenatal low-dose exposure to polychlorinated biphenyls (PCBs) and polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) is associated with altered DNAm patterns across the genome in a Western urban-industrial population. In 142 mother-infant pairs from the Duisburg Birth Cohort Study, PCBs and PCDD/Fs levels were quantified from maternal blood during late pregnancy and associated with DNAm levels in cord blood using the Illumina EPIC beadchip. The epigenome-wide association studies (EWAS) identified 32 significantly differentially methylated positions (DMPs) and eight differentially methylated regions (DMRs) associated with six congeners of PCB and PCDD in females or males (FDRs < 0.05). DMPs and DMRs mapped to genes involved in neurodevelopment, gene regulation, and immune functioning. Weighted gene correlation network analysis (WGCNA) showed 31 co-methylated modules (FDRs < 0.05) associated with one congener of PCDF levels in females. Results of both analytical strategies indicate that prenatal exposure to PCBs and PCDD/Fs is associated with altered DNAm of genes involved in neurodevelopment, gene expression and immune functioning. DNAm and gene expression levels of several of these genes were previously associated with EDC exposure in rodent models. Follow-up studies will clarify whether these epigenetic changes might contribute to the origin for adverse mental and health outcomes.

Abstract Background and aims Treatment with tumor necrosis factor α (TNFα) antagonists in IBD patients suffers from primary non-response rates of up to 40%. Biomarkers for early prediction of therapy success are missing. We investigated the dynamics of gene expression and DNA methylation in blood samples of IBD patients treated with the TNF antagonist infliximab and analyzed the predictive potential regarding therapy outcome. Methods We performed a longitudinal, blood-based multi-omics study in two prospective IBD patient cohorts receiving first-time infliximab therapy (discovery: 14 patients, replication: 23 patients). Samples were collected at up to 7 time points (from baseline to 14 weeks after therapy induction). RNA-sequencing and genome-wide DNA methylation data were analyzed and correlated with clinical remission at week 14 as a primary endpoint. Results We found no consistent ex ante predictive signature across the two cohorts. Longitudinally upregulated transcripts in the non-remitter group comprised TH2- and eosinophil-related genes including ALOX15, FCER1A, and OLIG2. Network construction identified transcript modules that were coherently expressed at baseline and in non-remitting patients but were disrupted at early time points in remitting patients. These modules reflected processes such as interferon signaling, erythropoiesis, and platelet aggregation. DNA methylation analysis identified remission-specific temporal changes, which partially overlapped with transcriptomic signals. Machine learning approaches identified features from differentially expressed genes cis-linked to DNA methylation changes at week 2 as a robust predictor of therapy outcome at week 14, which was validated in a publicly available dataset of 20 infliximab-treated CD patients. Conclusions Integrative multi-omics analysis reveals early shifts of gene expression and DNA methylation as predictors for efficient response to anti-TNF treatment. Lack of such signatures might be used to identify patients with IBD unlikely to benefit from TNF antagonists at an early time point.

Abstract In the development of orally inhaled drug products preclinical animal models regularly fail to predict pharmacological as well as toxicological responses in humans. Models based on human cells and tissues are potential alternatives to animal experimentation allowing for the isolation of essential processes of human biology and making them accessible in vitro. Here, the generation of a novel monoclonal cell line “Arlo,” derived from the polyclonal human alveolar epithelium lentivirus immortalized cell line hAELVi via single‐cell printing, and its characterization as a model for the human alveolar epithelium as well as a building block for future complex in vitro models is described. “Arlo” is systematically compared in vitro to primary human alveolar epithelial cells (hAEpCs) as well as to the polyclonal hAELVi cell line. “Arlo” cells show enhanced barrier properties with high transepithelial electrical resistance (TEER) of ≈3000 Ω cm2 and a potential difference (PD) of ≈30 mV under air–liquid interface (ALI) conditions, that can be modulated. The cells grow in a polarized monolayer and express genes relevant to barrier integrity as well as homeostasis as is observed in hAEpCs. Successful productive infection with severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) in a proof‐of‐principle study offers an additional, attractive application of “Arlo” beyond biopharmaceutical experimentation.

The composition of the gut microbiota is linked to multiple diseases, including Parkinson’s disease (PD). Abundance of bacteria producing short-chain fatty acids (SCFAs) and fecal SCFA concentrations are reduced in PD. SCFAs exert various beneficial functions in humans. In the interventional, monocentric, open-label clinical trial “Effects of Resistant Starch on Bowel Habits, Short Chain Fatty Acids and Gut Microbiota in Parkinson’s Disease” (RESISTA-PD; ID: NCT02784145), we aimed at altering fecal SCFAs by an 8-week prebiotic intervention with resistant starch (RS). We enrolled 87 subjects in three study-arms: 32 PD patients received RS (PD + RS), 30 control subjects received RS, and 25 PD patients received solely dietary instructions. We performed paired-end 100 bp length metagenomic sequencing of fecal samples using the BGISEQ platform at an average of 9.9 GB. RS was well-tolerated. In the PD + RS group, fecal butyrate concentrations increased significantly, and fecal calprotectin concentrations dropped significantly after 8 weeks of RS intervention. Clinically, we observed a reduction in non-motor symptom load in the PD + RS group. The reference-based analysis of metagenomes highlighted stable alpha-diversity and beta-diversity across the three groups, including bacteria producing SCFAs. Reference-free analysis suggested punctual, yet pronounced differences in the metagenomic signature in the PD + RS group. RESISTA-PD highlights that a prebiotic treatment with RS is safe and well-tolerated in PD. The stable alpha-diversity and beta-diversity alongside altered fecal butyrate and calprotectin concentrations call for long-term studies, also investigating whether RS is able to modify the clinical course of PD.

Abstract Background Understanding the influence of genetic variants on DNA methylation is fundamental for the interpretation of epigenomic data in the context of disease. There is a need for systematic approaches not only for determining methylation quantitative trait loci (methQTL), but also for discriminating general from cell type-specific effects. Results Here, we present a two-step computational framework MAGAR ( https://bioconductor.org/packages/MAGAR ), which fully supports the identification of methQTLs from matched genotyping and DNA methylation data, and additionally allows for illuminating cell type-specific methQTL effects. In a pilot analysis, we apply MAGAR on data in four tissues (ileum, rectum, T cells, B cells) from healthy individuals and demonstrate the discrimination of common from cell type-specific methQTLs. We experimentally validate both types of methQTLs in an independent data set comprising additional cell types and tissues. Finally, we validate selected methQTLs located in the PON1, ZNF155, and NRG2 genes by ultra-deep local sequencing. In line with previous reports, we find cell type-specific methQTLs to be preferentially located in enhancer elements. Conclusions Our analysis demonstrates that a systematic analysis of methQTLs provides important new insights on the influences of genetic variants to cell type-specific epigenomic variation.

Fibroblast-like synoviocytes are important mediators of joint pathology in rheumatoid arthritis (RA). Here the authors show that Lasp1 is epigenetically regulated and highly expressed by these cells in RA and its deletion can limit joint pathology in a mouse model of inflammatory arthritis.

Although epigenome-wide association studies of Alzheimer’s disease have highlighted neuropathology-associated DNA methylation differences, previous studies have been limited in sample size and brain region used. Here, the authors combine data from six DNA methylomic studies of Alzheimer’s disease (N = 1453 unique individuals) to identify differentially methylated loci across cortex.

Abstract Background High-grade serous ovarian cancer (HGSOC) is a major unmet need in oncology. The remaining uncertainty on its originating tissue has hampered the discovery of molecular oncogenic pathways and the development of effective therapies. Methods We used an approach based on the retention in tumors of a DNA methylation trace (OriPrint) that distinguishes the two putative tissues of origin of HGSOC, the fimbrial (FI) and ovarian surface epithelia (OSE), to stratify HGSOC by several clustering methods, both linear and non-linear. The identified tumor subtypes (FI-like and OSE-like HGSOC) were investigated at the RNAseq level to stratify an in-house cohort of macrodissected HGSOC FFPE samples to derive overall and disease-free survival and identify specific transcriptional alterations of the two tumor subtypes, both by classical differential expression and weighted correlation network analysis. We translated our strategy to published datasets and verified the co-occurrence of previously described molecular classification of HGSOC. We performed cytokine analysis coupled to immune phenotyping to verify alterations in the immune compartment associated with HGSOC. We identified genes that are both differentially expressed and methylated in the two tumor subtypes, concentrating on PAX8 as a bona fide marker of FI-like HGSOC. Results We show that: - OriPrint is a robust DNA methylation tracer that exposes the tissue of origin of HGSOC. - The tissue of origin of HGSOC is the main determinant of DNA methylation variance in HGSOC. - The tissue of origin is a prognostic factor for HGSOC patients. - FI-like and OSE-like HGSOC are endowed with specific transcriptional alterations that impact patients’ prognosis. - OSE-like tumors present a more invasive and immunomodulatory phenotype, compatible with its worse prognostic impact. - Among genes that are differentially expressed and regulated in FI-like and OSE-like HGSOC, PAX8 is a bona fide marker of FI-like tumors. Conclusions Through an integrated approach, our work demonstrates that both FI and OSE are possible origins for human HGSOC, whose derived subtypes are both molecularly and clinically distinct. These results will help define a new roadmap towards rational, subtype-specific therapeutic inroads and improved patients’ care.

Mutations in CLRN1 cause Usher syndrome (USH) type III (USH3A), a disease characterized by progressive hearing impairment, retinitis pigmentosa, and vestibular dysfunction. Due to the lack of appropriate disease models, no efficient therapy for retinitis pigmentosa in USH patients exists so far. In addition, given the yet undefined functional role and expression of the different CLRN1 splice isoforms in the retina, non-causative therapies such as gene supplementation are unsuitable at this stage. In this study, we focused on the recently identified deep intronic c.254-649T>G CLRN1 splicing mutation and aimed to establish two causative treatment approaches: CRISPR-Cas9-mediated excision of the mutated intronic region and antisense oligonucleotide (AON)-mediated correction of mRNA splicing. The therapeutic potential of these approaches was validated in different cell types transiently or stably expressing CLRN1 minigenes. Both approaches led to substantial correction of the splice defect. Surprisingly, however, no synergistic effect was detected when combining both methods. Finally, the injection of naked AONs into mice expressing the mutant CLRN1 minigene in the retina also led to a significant splice rescue. We propose that both AONs and CRISPR-Cas9 are suitable strategies to initiate advanced preclinical studies for treatment of USH3A patients.

Adoptive transfer of regulatory T cells (Treg) is a promising new therapeutic option to treat detrimental inflammatory conditions after transplantation and during autoimmune disease. To reach sufficient cell yield for treatment, ex vivo isolated autologous or allogenic Tregs need to be expanded extensively in vitro during manufacturing of the Treg product. However, repetitive cycles of restimulation and prolonged culture have been shown to impact T cell phenotypes, functionality and fitness. It is therefore critical to scrutinize the molecular changes which occur during T cell product generation, and reexamine current manufacturing practices. We performed genome-wide DNA methylation profiling of cells throughout the manufacturing process of a polyclonal Treg product that has proven safety and hints of therapeutic efficacy in kidney transplant patients. We found progressive DNA methylation changes over the duration of culture, which were donor-independent and reproducible between manufacturing runs. Differentially methylated regions (DMRs) in the final products were significantly enriched at promoters and enhancers of genes implicated in T cell activation. Additionally, significant hypomethylation did also occur in promoters of genes implicated in functional exhaustion in conventional T cells, some of which, however, have been reported to strengthen immunosuppressive effector function in Tregs. At the same time, a set of reported Treg-specific demethylated regions increased methylation levels with culture, indicating a possible destabilization of Treg identity during manufacturing, which was independent of the purity of the starting material. Together, our results indicate that the repetitive TCR-mediated stimulation lead to epigenetic changes that might impact functionality of Treg products in multiple ways, by possibly shifting to an effector Treg phenotype with enhanced functional activity or by risking destabilization of Treg identity and impaired TCR activation. Our analyses also illustrate the value of epigenetic profiling for the evaluation of T cell product manufacturing pipelines, which might open new avenues for the improvement of current adoptive Treg therapies with relevance for conventional effector T cell products.

Background: Based on reports on elevated cholesterol levels in cancer cells, strategies to lower cholesterol synthesis have been suggested as an antitumour strategy. However, cholesterol depletion has also been shown to induce tumour-promoting actions in tumour-associated macrophages (TAMs). Methods: We performed lipidomic and transcriptomic analyses of human lung cancer material. To assess whether the TAM phenotype is shaped by secreted factors produced by tumour cells, primary human monocyte-derived macrophages were polarized towards a TAM-like phenotype using tumour cell-conditioned medium. Findings: Lipidomic analysis of lung adenocarcinoma (n=29) and adjacent non-tumour tissues (n=22) revealed a significant accumulation of free cholesterol and cholesteryl esters within the tumour tissue. In contrast, cholesterol levels were reduced in TAMs isolated from lung adenocarcinoma tissues when compared with alveolar macrophages (AMs) obtained from adjacent non-tumour tissues. Bulk-RNA-Seq revealed that genes involved in cholesterol biosynthesis and metabolism were downregulated in TAMs, while cholesterol efflux transporters were upregulated. In vitro polarized TAM-like macrophages showed an attenuated lipogenic gene expression signature and exhibited lower cholesterol levels compared with non-polarized macrophages. A genome-wide comparison by bulk RNA-Seq confirmed a high similarity of ex vivo TAMs and in vitro TAM-like macrophages. Modulation of intracellular cholesterol levels by either starving, cholesterol depletion, or efflux transporter inhibition indicated that cholesterol distinctly shapes macrophage gene expression. Interpretation: Our data show an opposite dysregulation of cholesterol homeostasis in tumour tissue vs. TAMs. Polarization of in vitro differentiated macrophages by tumour cell-conditioned medium recapitulates key features of ex vivo TAMs. Funding: Deutsche Forschungsgemeinschaft (DFG), Landesforschungsf €orderungsprogramm Saarland (LFPP).

Summary: Epigenetic modifications (e.g. DNA methylation) in NAFLD and their contribution to disease progression and extrahepatic complications are poorly explored. Here, we use an integrated epigenome and transcriptome analysis of mouse NAFLD hepatocytes and identify alterations in glyoxylate metabolism, a pathway relevant in kidney damage via oxalate release—a harmful waste product and kidney stone-promoting factor. Downregulation and hypermethylation of alanine-glyoxylate aminotransferase (Agxt), which detoxifies glyoxylate, preventing excessive oxalate accumulation, is accompanied by increased oxalate formation after metabolism of the precursor hydroxyproline. Viral-mediated Agxt transfer or inhibiting hydroxyproline catabolism rescues excessive oxalate release. In human steatotic hepatocytes, AGXT is also downregulated and hypermethylated, and in NAFLD adolescents, steatosis severity correlates with urinary oxalate excretion. Thus, this work identifies a reduced capacity of the steatotic liver to detoxify glyoxylate, triggering elevated oxalate, and provides a mechanistic explanation for the increased risk of kidney stones and chronic kidney disease in NAFLD patients.

Abstract We present the software Condition-specific Regulatory Units Prediction (CRUP) to infer from epigenetic marks a list of regulatory units consisting of dynamically changing enhancers with their target genes. The workflow consists of a novel pre-trained enhancer predictor that can be reliably applied across cell types and species, solely based on histone modification ChIP-seq data. Enhancers are subsequently assigned to different conditions and correlated with gene expression to derive regulatory units. We thoroughly test and then apply CRUP to a rheumatoid arthritis model, identifying enhancer-gene pairs comprising known disease genes as well as new candidate genes.

Monozygotic (MZ) twins are ideal to study the influence of non-genetic factors on complex phenotypes. Here, Souren et al. perform an EWAS in peripheral blood mononuclear cells from 45 MZ twins discordant for multiple sclerosis and identify disease and treatment-associated epigenetic markers.

Abstract As epigenetic studies become more common and lead to new insights into health and disease, the return of individual epigenetic results to research participants, in particular in large-scale epigenomic studies, will be of growing importance. Members of the International Human Epigenome Consortium (IHEC) Bioethics Workgroup considered the potential ethical, legal, and social issues (ELSI) involved in returning epigenetic research results and incidental findings in order to produce a set of ‘Points-to-consider’ (P-t-C) for the epigenetics research community. These P-t-C draw on existing guidance on the return of genetic research results, while also integrating the IHEC Bioethics Workgroup’s ELSI research on and discussion of the issues associated with epigenetic data as well as the experience of a return of results pilot study by the Personal Genome Project UK (PGP-UK). Major challenges include how to determine the clinical validity and actionability of epigenetic results, and considerations related to environmental exposures and epigenetic marks, including circumstances warranting the sharing of results with family members and third parties. Interdisciplinary collaboration and good public communication regarding epigenetic risk will be important to advance the return of results framework for epigenetic science.

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.

Abstract Background Sarcoidosis is a systemic disease of unknown etiology. The disease mechanisms are largely speculative and may include the role microbial patterns that initiate and drive an underlying immune process. The aim of this study was to characterize the microbiota of the lung of patients with sarcoidosis and compare its composition and diversity with the results from patients with other interstitial lung disease (ILD) and historic healthy controls. Methods Patients (sarcoidosis, n = 31; interstitial lung disease, n = 19) were recruited within the PULMOHOM study, a prospective cohort study to characterize inflammatory processes in pulmonary diseases. Bronchoscopy of the middle lobe or the lingula was performed and the recovered fluid was immediately sent for analysis of the pulmonary microbiota by 16sRNA gene sequencing. Subsequent bioinformatic analysis was performed to compare the groups. Results There were no significant differences between patients with sarcoidosis or other ILDs with regard to microbiome composition and diversity. In addition, the abundance of the genera Atopobium, Fusobacterium, Mycobacterium or Propionibacterium were not different between the two groups. There were no gross differences to historical healthy controls. Conclusion The analysis of the pulmonary microbiota based on 16sRNA gene sequencing did not show a significant dysbiosis in patients with sarcoidosis as compared to other ILD patients. These data do not exclude a microbiological component in the pathogenesis of sarcoidosis.

Dimethyl sulfoxide (DMSO) is used to sustain or favor hepatocyte differentiation in vitro. Thus, DMSO is used in the differentiation protocol of the HepaRG cells that present the closest drug-metabolizing enzyme activities to primary human hepatocytes in culture. The aim of our study is to clarify its influence on liver-specific gene expression. For that purpose, we performed a large-scale analysis (gene expression and histone modification) to determine the global role of DMSO exposure during the differentiation process of the HepaRG cells. The addition of DMSO drives the upregulation of genes mainly regulated by PXR and PPARα whereas genes not affected by this addition are regulated by HNF1α, HNF4α, and PPARα. DMSO-differentiated-HepaRG cells show a differential expression for genes regulated by histone acetylation, while differentiated-HepaRG cells without DMSO show gene signatures associated with histone deacetylases. In addition, we observed an interplay between cytoskeleton organization and EMC remodeling with hepatocyte maturation.

The functions of bone marrow plasma cells (BMPC) beyond antibody production are not fully elucidated and distinct subsets of BMPC suggest potential different functions. Phenotypic differences were identified for human BMPC depending on CD19 expression. Since CD19 is a co-stimulatory molecule of the B-cell-receptor (BCR), and IgA+ and IgM+ BMPC express the BCR on their surface, we here studied whether CD19 expression affects cellular responses, such as BCR signaling and the expression of checkpoint molecules. We analyzed 132 BM samples from individuals undergoing routine total hip arthroplasty. We found that both CD19+ and CD19− BMPC expressed BCR signaling molecules. Notably, the BCR-associated kinase spleen tyrosine kinase (SYK) including pSYK was higher expressed in CD19+ BMPC compared to CD19− BMPC. BCR stimulation also resulted in increased kinase phosphorylation downstream of the BCR while expression of CD19 remained stable afterwards. Interestingly, the BCR response was restricted to IgA+ BMPC independently of CD19 expression. With regard to the expression of checkpoint molecules, CD19− BMPC expressed higher levels of co-inhibitory molecule programmed cell death protein-1 (PD-1) than CD19+ BMPC. IgA+ BMPC characteristically upregulated PD-1 upon BCR stimulation in contrast to other PC subsets and inhibition of the kinase SYK abrogated PD-1 upregulation. In contrast, expression of PD-1 ligand, B and T lymphocyte attenuator (BTLA) and CD28 did not change upon BCR activation of IgA+ BMPC. Here, we identify a distinct characteristic of IgA+ BMPC that is independent of the phenotypic heterogeneity of the subsets according to their CD19 expression. The data suggest that IgA+ BMPC underlie different regulatory principles and/or exert distinct regulatory functions.

CD4+ regulatory T cells (Tregs) are key mediators of immunological tolerance and promising effector cells for immuno-suppressive adoptive cellular therapy to fight autoimmunity and chronic inflammation. Their functional stability is critical for their clinical utility and has been correlated to the demethylated state of the TSDR/CNS2 enhancer element in the Treg lineage transcription factor FOXP3. However, proof for a causal contribution of the TSDR de-methylation to FOXP3 stability and Treg induction is so far lacking. We here established a powerful transient-transfection CRISPR-Cas9-based epigenetic editing method for the selective de-methylation of the TSDR within the endogenous chromatin environment of a living cell. The induced de-methylated state was stable over weeks in clonal T cell proliferation cultures even after expression of the editing complex had ceased. Epigenetic editing of the TSDR resulted in FOXP3 expression, even in its physiological isoform distribution, proving a causal role for the de-methylated TSDR in FOXP3 regulation. However, successful FOXP3 induction was not associated with a switch towards a functional Treg phenotype, in contrast to what has been reported from FOXP3 overexpression approaches. Thus, TSDR de-methylation is required, but not sufficient for a stable Treg phenotype induction. Therefore, targeted demethylation of the TSDR may be a critical addition to published in vitro Treg induction protocols which so far lack FOXP3 stability.

Understanding the factors that underlie the epigenetic regulation of genes is crucial to understand the gene regulatory machinery as a whole. Several experimental and computational studies examined the relationship between different factors involved. Here we investigate the relationship between transcription factors (TFs) and histone modifications (HMs), based on ChIP-seq data in cell lines. As it was shown that gene regulation by TFs differs depending on the CpG class of a promoter, we study the impact of the CpG content in promoters on the associations between TFs and HMs. We suggest an approach based on sparse linear regression models to infer associations between TFs and HMs with respect to CpG content. A study of the partial correlation of HMs for the two classes of high and low CpG content reveals possible CpG dependence and potential candidates for confounding factors in our models. We show that the models are accurate, inferred associations reflect known biological relationships, and we give new insight into associations with respect to CpG content. Moreover, analysis of a ChIP-seq dataset in HepG2 cells of the HM H3K122ac, an HM about little is known, reveals novel TF associations and supports a previously established link to active transcription.

Abstract Background Bidirectional promoters (BPs) are prevalent in eukaryotic genomes. However, it is poorly understood how the cell integrates different epigenomic information, such as transcription factor (TF) binding and chromatin marks, to drive gene expression at BPs. Single-cell sequencing technologies are revolutionizing the field of genome biology. Therefore, this study focuses on the integration of single-cell RNA-seq data with bulk ChIP-seq and other epigenetics data, for which single-cell technologies are not yet established, in the context of BPs. Results We performed integrative analyses of novel human single-cell RNA-seq (scRNA-seq) data with bulk ChIP-seq and other epigenetics data. scRNA-seq data revealed distinct transcription states of BPs that were previously not recognized. We find associations between these transcription states to distinct patterns in structural gene features, DNA accessibility, histone modification, DNA methylation and TF binding profiles. Conclusions Our results suggest that a complex interplay of all of these elements is required to achieve BP-specific transcriptional output in this specialized promoter configuration. Further, our study implies that novel statistical methods can be developed to deconvolute masked subpopulations of cells measured with different bulk epigenomic assays using scRNA-seq data.

Spatial mapping of genomic programs in tissue cells is an important step in the understanding of organ function and disease. Here, the authors provide a spatially resolved epigenomic and transcriptomic map of human liver and show porto-central gradients in metabolic and morphogen networks and transcription factor binding sites as a basis to better understand liver regeneration and function.

Abstract Background Partially methylated domains are extended regions in the genome exhibiting a reduced average DNA methylation level. They cover gene-poor and transcriptionally inactive regions and tend to be heterochromatic. We present a comprehensive comparative analysis of partially methylated domains in human and mouse cells, to identify structural and functional features associated with them. Results Partially methylated domains are present in up to 75% of the genome in human and mouse cells irrespective of their tissue or cell origin. Each cell type has a distinct set of partially methylated domains, and genes expressed in such domains show a strong cell type effect. The methylation level varies between cell types with a more pronounced effect in differentiating and replicating cells. The lowest level of methylation is observed in highly proliferating and immortal cancer cell lines. A decrease of DNA methylation within partially methylated domains tends to be linked to an increase in heterochromatic histone marks and a decrease of gene expression. Characteristic combinations of heterochromatic signatures in partially methylated domains are linked to domains of early and middle S-phase and late S-G2 phases of DNA replication. Conclusions Partially methylated domains are prominent signatures of long-range epigenomic organization. Integrative analysis identifies them as important general, lineage- and cell type-specific topological features. Changes in partially methylated domains are hallmarks of cell differentiation, with decreased methylation levels and increased heterochromatic marks being linked to enhanced cell proliferation. In combination with broad histone marks, partially methylated domains demarcate distinct domains of late DNA replication.

Abstract Background Epigenome-wide association studies (EWAS) based on human brain samples allow a deep and direct understanding of epigenetic dysregulation in Alzheimer’s disease (AD). However, strong variation of cell-type proportions across brain tissue samples represents a significant source of data noise. Here, we report the first EWAS based on sorted neuronal and non-neuronal (mostly glia) nuclei from postmortem human brain tissues. Results We show that cell sorting strongly enhances the robust detection of disease-related DNA methylation changes even in a relatively small cohort. We identify numerous genes with cell-type-specific methylation signatures and document differential methylation dynamics associated with aging specifically in neurons such as CLU, SYNJ2 and NCOR2 or in glia RAI1,CXXC5 and INPP5A. Further, we found neuron or glia-specific associations with AD Braak stage progression at genes such as MCF2L, ANK1, MAP2, LRRC8B, STK32C and S100B. A comparison of our study with previous tissue-based EWAS validates multiple AD-associated DNA methylation signals and additionally specifies their origin to neuron, e.g., HOXA3 or glia (ANK1). In a meta-analysis, we reveal two novel previously unrecognized methylation changes at the key AD risk genes APP and ADAM17. Conclusions Our data highlight the complex interplay between disease, age and cell-type-specific methylation changes in AD risk genes thus offering new perspectives for the validation and interpretation of large EWAS results.

Abstract Background The interplay of epigenetic processes and the intestinal microbiota may play an important role in intestinal development and homeostasis. Previous studies have established that the microbiota regulates a large proportion of the intestinal epithelial transcriptome in the adult host, but microbial effects on DNA methylation and gene expression during early postnatal development are still poorly understood. Here, we sought to investigate the microbial effects on DNA methylation and the transcriptome of intestinal epithelial cells (IECs) during postnatal development. Methods We collected IECs from the small intestine of each of five 1-, 4- and 12 to 16-week-old mice representing the infant, juvenile, and adult states, raised either in the presence or absence of a microbiota. The DNA methylation profile was determined using reduced representation bisulfite sequencing (RRBS) and the epithelial transcriptome by RNA sequencing using paired samples from each individual mouse to analyze the link between microbiota, gene expression, and DNA methylation. Results We found that microbiota-dependent and -independent processes act together to shape the postnatal development of the transcriptome and DNA methylation signatures of IECs. The bacterial effect on the transcriptome increased over time, whereas most microbiota-dependent DNA methylation differences were detected already early after birth. Microbiota-responsive transcripts could be attributed to stage-specific cellular programs during postnatal development and regulated gene sets involved primarily immune pathways and metabolic processes. Integrated analysis of the methylome and transcriptome data identified 126 genomic loci at which coupled differential DNA methylation and RNA transcription were associated with the presence of intestinal microbiota. We validated a subset of differentially expressed and methylated genes in an independent mouse cohort, indicating the existence of microbiota-dependent “functional” methylation sites which may impact on long-term gene expression signatures in IECs. Conclusions Our study represents the first genome-wide analysis of microbiota-mediated effects on maturation of DNA methylation signatures and the transcriptional program of IECs after birth. It indicates that the gut microbiota dynamically modulates large portions of the epithelial transcriptome during postnatal development, but targets only a subset of microbially responsive genes through their DNA methylation status.

Abstract Background Genetic predisposition and epigenetic alterations are both considered to contribute to sporadic neurodegenerative diseases (NDDs) such as Parkinson’s disease (PD). Since cell reprogramming and the generation of induced pluripotent stem cells (iPSCs) are themselves associated with major epigenetic remodeling, it remains unclear to what extent iPSC-derived neurons lend themselves to model epigenetic disease-associated changes. A key question to be addressed in this context is whether iPSC-derived neurons exhibit epigenetic signatures typically observed in neurons derived from non-reprogrammed human embryonic stem cells (hESCs). Results Here, we compare mature neurons derived from hESC and isogenic human iPSC generated from hESC-derived neural stem cells. Genome-wide 450 K-based DNA methylation and HT12v4 gene array expression analyses were complemented by a deep analysis of selected genes known to be involved in NDD. Our studies show that DNA methylation and gene expression patterns of isogenic hESC- and iPSC-derived neurons are markedly preserved on a genome-wide and single gene level. Conclusions Overall, iPSC-derived neurons exhibit similar DNA methylation patterns compared to isogenic hESC-derived neurons. Further studies will be required to explore whether the epigenetic patterns observed in iPSC-derived neurons correspond to those detectable in native brain neurons.

Pluripotent cells emanate from the inner cell mass (ICM) of the blastocyst and when cultivated under optimal conditions immortalize as embryonic stem cells (ESCs). The fundamental mechanism underlying ESC derivation has, however, remained elusive. Recently, we have devised a highly efficient approach for establishing ESCs, through inhibition of the MEK and TGF-β pathways. This regimen provides a platform for dissecting the molecular mechanism of ESC derivation. Via temporal gene expression analysis, we reveal key genes involved in the ICM to ESC transition. We found that DNA methyltransferases play a pivotal role in efficient ESC generation. We further observed a tight correlation between ESCs and preimplantation epiblast cell-related genes and noticed that fundamental events such as epithelial-to-mesenchymal transition blockage play a key role in launching the ESC self-renewal program. Our study provides a time course transcriptional resource highlighting the dynamics of the gene regulatory network during the ICM to ESC transition.

The long non-coding RNA (lncRNA) H19 represents a maternally expressed and epigenetically regulated imprinted gene product and is discussed to have either tumor-promoting or tumor-suppressive actions. Recently, H19 was shown to be regulated under inflammatory conditions. Therefore, aim of this study was to determine the function of H19 in hepatocellular carcinoma (HCC), an inflammation-associated type of tumor. In four different human HCC patient cohorts H19 was distinctly downregulated in tumor tissue compared to normal or non-tumorous adjacent tissue. We therefore determined the action of H19 in three different human hepatoma cell lines (HepG2, Plc/Prf5, and Huh7). Clonogenicity and proliferation assays showed that H19 overexpression could suppress tumor cell survival and proliferation after treatment with either sorafenib or doxorubicin, suggesting chemosensitizing actions of H19. Since HCC displays a highly chemoresistant tumor entity, cell lines resistant to doxorubicin or sorafenib were established. In all six chemoresistant cell lines H19 expression was significantly downregulated. The promoter methylation of the H19 gene was significantly different in chemoresistant cell lines compared to their sensitive counterparts. Chemoresistant cells were sensitized after H19 overexpression by either increasing the cytotoxic action of doxorubicin or decreasing cell proliferation upon sorafenib treatment. An H19 knockout mouse model (H19Δ3) showed increased tumor development and tumor cell proliferation after treatment with the carcinogen diethylnitrosamine (DEN) independent of the reciprocally imprinted insulin-like growth factor 2 (IGF2). In conclusion, H19 suppresses hepatocarcinogenesis, hepatoma cell growth, and HCC chemoresistance. Thus, mimicking H19 action might be a potential target to overcome chemoresistance in future HCC therapy.

The insulin-like growth factor 2 (IGF2) mRNA binding proteins (IMPs/IGF2BPs) IMP1 and 3 are regarded as oncofetal proteins, whereas the hepatic IMP2 expression in adults is controversially discussed. The splice variant IMP2-2/p62 promotes steatohepatitis and hepatocellular carcinoma. Aim of this study was to clarify whether IMP2 is expressed in the adult liver and influences progression toward cirrhosis. IMP2 was expressed at higher levels in embryonic compared to adult tissues as quantified in embryonic, newborn, and adult C57BL/6J mouse livers and suggested by analysis of publicly available human data. In an IMP2-2 transgenic mouse model microarray and qPCR analyses revealed increased expression of liver progenitor cell (LPC) markers Bex1, Prom1, Spp1, and Cdh1 indicating a de-differentiated liver cell phenotype. Induction of these LPC markers was confirmed in human cirrhotic tissue datasets. The LPC marker SPP1 has been described to play a major role in fibrogenesis. Thus, DNA methylation was investigated in order to decipher the regulatory mechanism of Spp1 induction. In IMP2-2 transgenic mouse livers single CpG sites were differentially methylated, as quantified by amplicon sequencing, whereas human HCC samples of a human publicly available dataset showed promoter hypomethylation. In order to study the impact of IMP2 on fibrogenesis in the context of steatohepatitis wild-type or IMP2-2 transgenic mice were fed either a methionine-choline deficient (MCD) or a control diet for 2–12 weeks. MCD-fed IMP2-2 transgenic mice showed a higher incidence of ductular reaction (DR), accompanied by hepatic stellate cell activation, extracellular matrix (ECM) deposition, and induction of the LPC markers Spp1, Cdh1, and Afp suggesting the occurrence of de-differentiated cells in transgenic livers. In human cirrhotic samples IMP2 overexpression correlated with LPC marker and ECM component expression. Progression of liver disease was induced by combined MCD and diethylnitrosamine (DEN) treatment. Combined MCD-DEN treatment resulted in shorter survival of IMP2-2 transgenic compared to wild-type mice. Only IMP2-2 transgenic livers progressed to cirrhosis, which was accompanied by strong DR. In conclusion, IMP2 is an oncofetal protein in the liver that promotes DR characterized by de-differentiated cells toward steatohepatitis-associated cirrhosis development with poor survival.

Autoimmune diseases (AID) such as systemic lupus erythematosus (SLE), primary Sjögren's syndrome (pSS), and rheumatoid arthritis (RA) are chronic inflammatory diseases in which abnormalities of B cell function play a central role. Although it is widely accepted that autoimmune B cells are hyperactive in vivo, a full understanding of their functional status in AID has not been delineated. Here, we present a detailed analysis of the functional capabilities of AID B cells and dissect the mechanisms underlying altered B cell function. Upon BCR activation, decreased spleen tyrosine kinase (Syk) and Bruton's tyrosine kinase (Btk) phosphorylation was noted in AID memory B cells combined with constitutive co-localization of CD22 and protein tyrosine phosphatase (PTP) non-receptor type 6 (SHP-1) along with hyporesponsiveness to TLR9 signaling, a Syk-dependent response. Similar BCR hyporesponsiveness was also noted specifically in SLE CD27− B cells together with increased PTP activities and increased transcripts for PTPN2, PTPN11, PTPN22, PTPRC, and PTPRO in SLE B cells. Additional studies revealed that repetitive BCR stimulation of normal B cells can induce BCR hyporesponsiveness and that tissue-resident memory B cells from AID patients also exhibited decreased responsiveness immediately ex vivo, suggesting that the hyporesponsive status can be acquired by repeated exposure to autoantigen(s) in vivo. Functional studies to overcome B cell hyporesponsiveness revealed that CD40 co-stimulation increased BCR signaling, induced proliferation, and downregulated PTP expression (PTPN2, PTPN22, and receptor-type PTPs). The data support the conclusion that hyporesponsiveness of AID and especially SLE B cells results from chronic in vivo stimulation through the BCR without T cell help mediated by CD40–CD154 interaction and is manifested by decreased phosphorylation of BCR-related proximal signaling molecules and increased PTPs. The hyporesponsiveness of AID B cells is similar to a form of functional anergy.

Mouse embryonic stem cells are dynamic and heterogeneous. For example, rare cells cycle through a state characterized by decondensed chromatin and expression of transcripts, including the Zscan4 cluster and MERVL endogenous retrovirus, which are usually restricted to preimplantation embryos. Here, we further characterize the dynamics and consequences of this transient cell state. Single-cell transcriptomics identified the earliest upregulated transcripts as cells enter the MERVL/Zscan4 state. The MERVL/Zscan4 transcriptional network was also upregulated during induced pluripotent stem cell reprogramming. Genome-wide DNA methylation and chromatin analyses revealed global DNA hypomethylation accompanying increased chromatin accessibility. This transient DNA demethylation was driven by a loss of DNA methyltransferase proteins in the cells and occurred genome-wide. While methylation levels were restored once cells exit this state, genomic imprints remained hypomethylated, demonstrating a potential global and enduring influence of endogenous retroviral activation on the epigenome.

Co-localizing chromatin modifications and regulators can exert a combinatorial effect on chromatin structure and function. Here the authors describe reChIP-seq and normR to identify co-localizing proteins in an unbiased genome-wide manner.

Regulatory T cells (Tregs) are an attractive therapeutic tool for several different immune pathologies. Therapeutic Treg application often requires prolonged in vitro culture to generate sufficient Treg numbers or to optimize their functionality, e.g., via genetic engineering of their antigen receptors. However, purity of clinical Treg expansion cultures is highly variable, and currently, it is impossible to identify and separate stable Tregs from contaminating effector T cells, either ex vivo or after prior expansion. This represents a major obstacle for quality assurance of expanded Tregs and raises significant safety concerns. Here, we describe a Treg activation signature that allows identification and sorting of epigenetically imprinted Tregs even after prolonged in vitro culture. We show that short-term reactivation resulted in expression of CD137 but not CD154 on stable FoxP3+ Tregs that displayed a demethylated Treg-specific demethylated region, high suppressive potential, and lack of inflammatory cytokine expression. We also applied this Treg activation signature for rapid testing of chimeric antigen receptor functionality in human Tregs and identified major differences in the signaling requirements regarding CD137 versus CD28 costimulation. Taken together, CD137+CD154− expression emerges as a universal Treg activation signature ex vivo and upon in vitro expansion allowing the identification and isolation of epigenetically stable antigen-activated Tregs and providing a means for their rapid functional testing in vitro.

DNA methylation and demethylation are opposing processes that when in balance create stable patterns of epigenetic memory. The control of DNA methylation pattern formation by replication dependent and independent demethylation processes has been suggested to be influenced by Tet mediated oxidation of 5mC. Several alternative mechanisms have been proposed suggesting that 5hmC influences either replication dependent maintenance of DNA methylation or replication independent processes of active demethylation. Using high resolution hairpin oxidative bisulfite sequencing data, we precisely determine the amount of 5mC and 5hmC and model the contribution of 5hmC to processes of demethylation in mouse ESCs. We develop an extended hidden Markov model capable of accurately describing the regional contribution of 5hmC to demethylation dynamics. Our analysis shows that 5hmC has a strong impact on replication dependent demethylation, mainly by impairing methylation maintenance.

Memory type 1 T helper (T(H)1) cells are characterized by the stable expression of interferon (IFN)-γ as well as by the epigenetic imprinting of the IFNG locus. Among innate cells, NK cells play a crucial role in the defense against cytomegalovirus (CMV) and represent the main source of IFN-γ. Recently, it was shown that memory-like features can be observed in NK cell subsets after CMV infection. However, the molecular mechanisms underlying NK cell adaptive properties have not been completely defined. In the present study, we demonstrated that only NKG2Chi NK cells expanded in human CMV (HCMV) seropositive individuals underwent epigenetic remodeling of the IFNG conserved non-coding sequence (CNS) 1, similar to memory CD8(+) T cells or T(H)1 cells. The accessibility of the CNS1 was required to enhance IFN-γ transcriptional activity in response to NKG2C and 2B4 engagement, which led to consistent IFN-γ production in NKG2C(hi) NK cells. Thus, our data identify epigenetic imprinting of the IFNG locus as selective hallmark and crucial mechanism driving strong and stable IFN-γ expression in HCMV-specific NK cell expansions, providing a molecular basis for the regulation of adaptive features in innate cells.

Aberrant CpG methylation is a universal epigenetic trait of cancer cell genomes. However, human cancer samples or cell lines preclude the investigation of epigenetic changes occurring early during tumour development. Here, we have used MeDIP-seq to analyse the DNA methylome of APC(Min) adenoma as a model for intestinal cancer initiation, and we present a list of more than 13,000 recurring differentially methylated regions (DMRs) characterizing intestinal adenoma of the mouse. We show that Polycomb Repressive Complex (PRC) targets are strongly enriched among hypermethylated DMRs, and several PRC2 components and DNA methyltransferases were up-regulated in adenoma. We further demonstrate by bisulfite pyrosequencing of purified cell populations that the DMR signature arises de novo in adenoma cells rather than by expansion of a pre-existing pattern in intestinal stem cells or undifferentiated crypt cells. We found that epigenetic silencing of tumour suppressors, which occurs frequently in colon cancer, was rare in adenoma. Quite strikingly, we identified a core set of DMRs, which is conserved between mouse adenoma and human colon cancer, thus possibly revealing a global panel of epigenetically modified genes for intestinal tumours. Our data allow a distinction between early conserved epigenetic alterations occurring in intestinal adenoma and late stochastic events promoting colon cancer progression, and may facilitate the selection of more specific clinical epigenetic biomarkers.

The enzymatic control of the setting and maintenance of symmetric and non-symmetric DNA methylation patterns in a particular genome context is not well understood. Here, we describe a comprehensive analysis of DNA methylation patterns generated by high resolution sequencing of hairpin-bisulfite amplicons of selected single copy genes and repetitive elements (LINE1, B1, IAP-LTR-retrotransposons, and major satellites). The analysis unambiguously identifies a substantial amount of regional incomplete methylation maintenance, i.e. hemimethylated CpG positions, with variant degrees among cell types. Moreover, non-CpG cytosine methylation is confined to ESCs and exclusively catalysed by Dnmt3a and Dnmt3b. This sequence position-, cell type-, and region-dependent non-CpG methylation is strongly linked to neighboring CpG methylation and requires the presence of Dnmt3L. The generation of a comprehensive data set of 146,000 CpG dyads was used to apply and develop parameter estimated hidden Markov models (HMM) to calculate the relative contribution of DNA methyltransferases (Dnmts) for de novo and maintenance DNA methylation. The comparative modelling included wild-type ESCs and mutant ESCs deficient for Dnmt1, Dnmt3a, Dnmt3b, or Dnmt3a/3b, respectively. The HMM analysis identifies a considerable de novo methylation activity for Dnmt1 at certain repetitive elements and single copy sequences. Dnmt3a and Dnmt3b contribute de novo function. However, both enzymes are also essential to maintain symmetrical CpG methylation at distinct repetitive and single copy sequences in ESCs.

Differential DNA methylation is an essential epigenetic signal for gene regulation, development, and disease processes. We mapped DNA methylation patterns of 190 gene promoter regions on chromosome 21 using bisulfite conversion and subclone sequencing in five human cell types. A total of 28,626 subclones were sequenced at high accuracy using (long-read) Sanger sequencing resulting in the measurement of the DNA methylation state of 580427 CpG sites. Our results show that average DNA methylation levels are distributed bimodally with enrichment of highly methylated and unmethylated sequences, both for amplicons and individual subclones, which represent single alleles from individual cells. Within CpG-rich sequences, DNA methylation was found to be anti-correlated with CpG dinucleotide density and GC content, and methylated CpGs are more likely to be flanked by AT-rich sequences. We observed over-representation of CpG sites in distances of 9, 18, and 27 bps in highly methylated amplicons. However, DNA sequence alone is not sufficient to predict an amplicon's DNA methylation status, since 43% of all amplicons are differentially methylated between the cell types studied here. DNA methylation in promoter regions is strongly correlated with the absence of gene expression and low levels of activating epigenetic marks like H3K4 methylation and H3K9 and K14 acetylation. Utilizing the single base pair and single allele resolution of our data, we found that i) amplicons from different parts of a CpG island frequently differ in their DNA methylation level, ii) methylation levels of individual cells in one tissue are very similar, and iii) methylation patterns follow a relaxed site-specific distribution. Furthermore, iv) we identified three cases of allele-specific DNA methylation on chromosome 21. Our data shed new light on the nature of methylation patterns in human cells, the sequence dependence of DNA methylation, and its function as epigenetic signal in gene regulation. Further, we illustrate genotype-epigenotype interactions by showing novel examples of allele-specific methylation.

CpG islands were originally identified by epigenetic and functional properties, namely, absence of DNA methylation and frequent promoter association. However, this concept was quickly replaced by simple DNA sequence criteria, which allowed for genome-wide annotation of CpG islands in the absence of large-scale epigenetic datasets. Although widely used, the current CpG island criteria incur significant disadvantages: (1) reliance on arbitrary threshold parameters that bear little biological justification, (2) failure to account for widespread heterogeneity among CpG islands, and (3) apparent lack of specificity when applied to the human genome. This study is driven by the idea that a quantitative score of "CpG island strength" that incorporates epigenetic and functional aspects can help resolve these issues. We construct an epigenome prediction pipeline that links the DNA sequence of CpG islands to their epigenetic states, including DNA methylation, histone modifications, and chromatin accessibility. By training support vector machines on epigenetic data for CpG islands on human Chromosomes 21 and 22, we identify informative DNA attributes that correlate with open versus compact chromatin structures. These DNA attributes are used to predict the epigenetic states of all CpG islands genome-wide. Combining predictions for multiple epigenetic features, we estimate the inherent CpG island strength for each CpG island in the human genome, i.e., its inherent tendency to exhibit an open and transcriptionally competent chromatin structure. We extensively validate our results on independent datasets, showing that the CpG island strength predictions are applicable and informative across different tissues and cell types, and we derive improved maps of predicted "bona fide" CpG islands. The mapping of CpG islands by epigenome prediction is conceptually superior to identifying CpG islands by widely used sequence criteria since it links CpG island detection to their characteristic epigenetic and functional states. And it is superior to purely experimental epigenome mapping for CpG island detection since it abstracts from specific properties that are limited to a single cell type or tissue. In addition, using computational epigenetics methods we could identify high correlation between the epigenome and characteristics of the DNA sequence, a finding which emphasizes the need for a better understanding of the mechanistic links between genome and epigenome.

CpG island methylation plays an important role in epigenetic gene control during mammalian development and is frequently altered in disease situations such as cancer. The majority of CpG islands is normally unmethylated, but a sizeable fraction is prone to become methylated in various cell types and pathological situations. The goal of this study is to show that a computational epigenetics approach can discriminate between CpG islands that are prone to methylation from those that remain unmethylated. We develop a bioinformatics scoring and prediction method on the basis of a set of 1,184 DNA attributes, which refer to sequence, repeats, predicted structure, CpG islands, genes, predicted binding sites, conservation, and single nucleotide polymorphisms. These attributes are scored on 132 CpG islands across the entire human Chromosome 21, whose methylation status was previously established for normal human lymphocytes. Our results show that three groups of DNA attributes, namely certain sequence patterns, specific DNA repeats, and a particular DNA structure, are each highly correlated with CpG island methylation (correlation coefficients of 0.64, 0.66, and 0.49, respectively). We predicted, and subsequently experimentally examined 12 CpG islands from human Chromosome 21 with unknown methylation patterns and found more than 90% of our predictions to be correct. In addition, we applied our prediction method to analyzing Human Epigenome Project methylation data on human Chromosome 6 and again observed high prediction accuracy. In summary, our results suggest that DNA composition of CpG islands (sequence, repeats, and structure) plays a significant role in predisposing CpG islands for DNA methylation. This finding may have a strong impact on our understanding of changes in CpG island methylation in development and disease.