BiQ Analyzer HT: locus-specific analysis of DNA methylation by high-throughput bisulfite sequencing

Autoantigen-specific CD4+ T cells acquire an exhausted phenotype and persist in human antigen-specific autoimmune diseases

Pro-inflammatory autoantigen-specific CD4+ T helper (auto-Th) cells are central orchestrators of autoimmune diseases (AIDs). We aimed to characterize these cells in human AIDs with defined autoantigens by combining human leukocyte antigen (HLA)-tetramer-based and activation-based multidimensional ex vivo analyses. In aquaporin4-antibody-positive neuromyelitis optica spectrum disorder (AQP4-NMOSD) patients, auto-Th cells expressed CD154, but proliferative capacity and pro-inflammatory cytokines were strongly reduced. Instead, exhaustion-associated co-inhibitory receptors were expressed together with FOXP3, the canonical regulatory T cell (Treg) transcription factor. Auto-Th cells responded in vitro to checkpoint inhibition and provided potent B cell help. Cells with the same exhaustion-like (ThEx) phenotype were identified in soluble liver antigen (SLA)-antibody-autoimmune hepatitis and BP180-antibody-positive bullous pemphigoid, AIDs of the liver and skin, respectively. While originally described in cancer and chronic infection, our data point to T cell exhaustion as a common mechanism of adaptation to chronic (self-)stimulation across AID types and link exhausted CD4+ T cells to humoral autoimmune responses, with implications for therapeutic targeting.

Targeted DNA Methylation Editing Using an All-in-One System Establishes Paradoxical Activation of EBF3

Cutaneous melanoma is rapidly on the rise globally, surpassing the growth rate of other cancers, with metastasis being the primary cause of death in melanoma patients. Consequently, understanding the mechanisms behind this metastatic process and exploring innovative treatments is of paramount importance. Recent research has shown promise in unravelling the role of epigenetic factors in melanoma progression to metastasis. While DNA hypermethylation at gene promoters typically suppresses gene expression, we have contributed to establishing the newly understood mechanism of paradoxical activation of genes via DNA methylation, where high methylation coincides with increased gene activity. This mechanism challenges the conventional paradigm that promoter methylation solely silences genes, suggesting that, for specific genes, it might actually activate them. Traditionally, altering DNA methylation in vitro has involved using global demethylating agents, which is insufficient for studying the mechanism and testing the direct consequence of gene methylation changes. To investigate promoter hypermethylation and its association with gene activation, we employed a novel approach utilising a CRISPR-SunTag All-in-one system. Here, we focused on editing the DNA methylation of a specific gene promoter segment (EBF3) in melanoma cells using the All-in-one system. Using bisulfite sequencing and qPCR with RNA-Seq, we successfully demonstrated highly effective methylation and demethylation of the EBF3 promoter, with subsequent gene expression changes, to establish and validate the paradoxical role of DNA methylation. Further, our study provides novel insights into the function of the EBF3 gene, which remains largely unknown. Overall, this study challenges the conventional view of methylation as solely a gene-silencing mechanism and demonstrates a potential function of EBF3 in IFN pathway signalling, potentially uncovering new insights into epigenetic drivers of malignancy and metastasis.

Inheritance of perturbed methylation and metabolism caused by uterine malnutrition via oocytes

BACKGROUND

Undernourishment in utero has deleterious effects on the metabolism of offspring, but the mechanism of the transgenerational transmission of metabolic disorders is not well known. In the present study, we found that undernourishment in utero resulted in metabolic disorders of female F1 and F2 in mouse model.

RESULTS

Undernutrition in utero induced metabolic disorders of F1 females, which was transmitted to F2 females. The global methylation in oocytes of F1 exposed to undernutrition in utero was decreased compared with the control. KEGG analysis showed that genes with differential methylation regions (DMRs) in promoters were significantly enriched in metabolic pathways. The altered methylation of some DMRs in F1 oocytes located at the promoters of metabolic-related genes were partially observed in F2 tissues, and the expressions of these genes were also changed. Meanwhile, the abnormal DNA methylation of the validated DMRs in F1 oocytes was also observed in F2 oocytes.

CONCLUSIONS

These results indicate that DNA methylation may mediate the transgenerational inheritance of metabolic disorders induced by undernourishment in utero via female germline.

mir152 hypomethylation as a mechanism for non-syndromic cleft lip and palate

Non-syndromic cleft lip with or without cleft palate (NSCLP), the most common human craniofacial malformation, is a complex disorder given its genetic heterogeneity and multifactorial component revealed by genetic, epidemiological, and epigenetic findings. Epigenetic variations associated with NSCLP have been identified; however, functional investigation has been limited. Here, we combined a reanalysis of NSCLP methylome data with genetic analysis and used both in vitro and in vivo approaches to dissect the functional effects of epigenetic changes. We found a region in mir152 that is frequently hypomethylated in NSCLP cohorts (21-26%), leading to mir152 overexpression. mir152 overexpression in human neural crest cells led to downregulation of spliceosomal, ribosomal, and adherens junction genes. In vivo analysis using zebrafish embryos revealed that mir152 upregulation leads to craniofacial cartilage impairment. Also, we suggest that zebrafish embryonic hypoxia leads to mir152 upregulation combined with mir152 hypomethylation and also analogous palatal alterations. We therefore propose that mir152 hypomethylation, potentially induced by hypoxia in early development, is a novel and frequent predisposing factor to NSCLP.

A transition phase in late mouse oogenesis impacts DNA methylation of the early embryo

A well-orchestrated program of oocyte growth and differentiation results in a developmentally competent oocyte. In late oogenesis, germinal vesicle oocytes (GVOs) undergo chromatin remodeling accompanied by transcriptional silencing from an NSN (non-surrounded nucleolus) to an SN (surrounded nucleolus) chromatin state. By analyzing different cytoplasmic and nuclear characteristics, our results indicate that murine NSN-GVOs transition via an intermediate stage into SN-GVOs in vivo. Interestingly, this transition can also be observed ex vivo, including most characteristics seen in vivo, which allows to analyze this transition process in more detail. The nuclear rearrangements during the transition are accompanied by changes in DNA methylation and Tet enzyme-catalyzed DNA modifications. Early parthenogenetic embryos, derived from NSN-GVOs, show lower DNA methylation levels than SN-derived embryos. Together, our data suggest that a successful NSN-SN transition in oogenesis including proper DNA methylation remodeling is important for the establishment of a developmentally competent oocyte for the beginning of life.

Salivary DNA methylation markers for cancer of oral cavity

PURPOSE: Aberrant DNA methylation plays a crucial role in oral carcinogenesis. Our previous study demonstrated hypermethylation of DAPK1, LRPPRC, RAB6C, and ZNF471 promoters in patients with tongue squamous cell carcinoma compared with normal samples. Methylation profiling using salivary DNA is considered a non-invasive alternative to tissue samples. Hence, the present study tested the DNA methylation status of these four promoters as indicators of oral cancer progression. METHODS: We performed the bisulfite-based targeted next-generation sequencing of four candidate genes in saliva and tissue DNA from normal, premalignant, and squamous cell carcinoma subjects. The clinicopathological association, diagnostic, and prognostic utility of aberrant DNA methylation were evaluated using the TCGA-HNSCC dataset. Using the Xgboost algorithm and logistic regression, CpG sites were prioritized, and Receiver Operating Characteristic was generated. By Log-rank test and Kaplan-Meier (KM) curves, an association between methylation and overall survival (OS), disease-free interval (DFI), and progression-free interval (PFI) were computed. RESULTS: We identified all four genes as significantly hypermethylation in premalignant and malignant samples compared with normal samples. The methylation levels were comparable between saliva and tissue samples with an r-value of 0.6297 to 0.8023 and 0.7823 to 0.9419 between premalignant tissue vs. saliva and OC vs. saliva, respectively. We identified an inverse correlation between DAPK1, LRPPRC, RAB6C, and ZNF471 promoter methylation with its expression. A classifier of 8 differentially methylated CpG sites belonging to DAPK1, RAB6C, and ZNF471 promoters was constructed, showing an AUC of 0.984 to differentiate tumors from normal samples. The differential methylation status of DAPK1, LRPPRC, and ZNF71 promoters was prognostically important. Abnormal expression of all four genes was associated with immune infiltration. CONCLUSIONS: Thus, methylation analysis of these candidate CpG sites from saliva can be helpful as a non-invasive tool for the clinical management of OC.

Heritable Epigenomic Modifications Influence Stress Resilience and Rapid Adaptations in the Brown Planthopper (Nilaparvata lugens)

DNA methylation in insects is integral to cellular differentiation, development, gene regulation, genome integrity, and phenotypic plasticity. However, its evolutionary potential and involvement in facilitating rapid adaptations in insects are enigmatic. Moreover, our understanding of these mechanisms is limited to a few insect species, of which none are pests of crops. Hence, we studied methylation patterns in the brown planthopper (BPH), a major rice pest, under pesticide and nutritional stress, across its life stages. Moreover, as the inheritance of epigenetic changes is fundamentally essential for acclimation, adaptability, and evolution, we determined the heritability and persistence of stress-induced methylation marks in BPH across generations. Our results revealed that DNA methylation pattern(s) in BPH varies/vary with environmental cues and is/are insect life-stage specific. Further, our findings provide novel insights into the heritability of stress-induced methylation marks in BPH. However, it was observed that, though heritable, these marks eventually fade in the absence of the stressors, thereby suggesting the existence of fitness cost(s) associated with the maintenance of the stressed epigenotype. Furthermore, we demonstrate how 5-azacytidine-mediated disruption of BPH methylome influences expression levels of stress-responsive genes and, thereby, highlight demethylation/methylation as a phenomenon underlying stress resilience of BPH.

Next-Generation Bisulfite Sequencing for Targeted DNA Methylation Analysis

Bisulfite sequencing is the "gold-standard" technique for DNA methylation analysis. By combining bisulfite sequencing with high-throughput, next-generation sequencing technology, we can document methylation from many thousands of individual reads (equivalent to alleles or "cells"), for multiple target regions and from many samples simultaneously. Here, we describe a next-generation bisulfite-sequencing assay for targeted DNA methylation analysis which offers scope for the simultaneous interrogation of multiple genomic loci across numerous samples.

Cigarette heavy smoking alters DNA methylation patterns and gene transcription levels in humans spermatozoa

Tobacco smoking is considered the most common reason of death and infertility around the world. This study was designed to assess the impact of tobacco heavy smoking on sperm DNA methylation patterns and to determine whether the transcription level of ALDH3B2, PTGIR, PRICKLE2, and ALS2CR12 genes is different in heavy smokers compared to non-smokers. As a screening study, the 450 K array was used to assess the alteration in DNA methylation patterns between heavy smokers (n = 15) and non-smokers (n = 15). Then, four CpGs that have the highest difference in methylation level (cg16338278, cg08408433, cg05799088, and cg07227024) were selected for validation using deep bisulfite sequencing in an independent cohort of heavy smokers (n = 200) and non-smokers (n = 100). A significant variation was found between heavy smokers and non-smokers in the methylation level at all CpGs within the PRICKLE2 and ALS2CR12 gene amplicon (P < 0.001). Similarly, a significant variation was found in the methylation level at nine out of thirteen CpGs within the ALDH3B2 gene amplicon (P < 0.01). Additionally, eighteen CpGs out of the twenty-six within the PTGIR gene amplicon have a significant difference in the methylation level between heavy smokers and non-smokers (P < 0.01). The study showed a significant difference in sperm global DNA methylation, chromatin non-condensation, and DNA fragmentation (P < 0.001) between heavy smokers and non-smokers. A significant decline was shown in the transcription level of ALDH3B2, PTGIR, PRICKLE2, and ALS2CR12 genes (P < 0.001) in heavy smokers. In conclusion, heavy smoking influences DNA methylation at several CpGs, sperm global DNA methylation, and transcription level of the PRICKLE2, ALS2CR12, ALDH3B2, and PTGIR genes, which affects negatively the semen parameters of heavy smokers.

Bisulfite profiling of the MGMT promoter and comparison with routine testing in glioblastoma diagnostics

Background

Promoter methylation of the DNA repair gene O⁶-methylguanine-DNA methyltransferase (MGMT) is an acknowledged predictive epigenetic marker in glioblastoma multiforme and anaplastic astrocytoma. Patients with methylated CpGs in the MGMT promoter benefit from treatment with alkylating agents, such as temozolomide, and show an improved overall survival and progression-free interval. A precise determination of MGMT promoter methylation is of importance for diagnostic decisions. We experienced that different methods show partially divergent results in a daily routine. For an integrated neuropathological diagnosis of malignant gliomas, we therefore currently apply a combination of methylation-specific PCR assays and pyrosequencing.

Results

To better rationalize the variation across assays, we compared these standard techniques and assays to deep bisulfite sequencing results in a cohort of 80 malignant astrocytomas. Our deep analysis covers 49 CpG sites of the expanded MGMT promoter, including exon 1, parts of intron 1 and a region upstream of the transcription start site (TSS). We observed that deep sequencing data are in general in agreement with CpG-specific pyrosequencing, while the most widely used MSP assays published by Esteller et al. (N Engl J Med 343(19):1350–1354, 2000. 10.1056/NEJM200011093431901) and Felsberg et al. (Clin Cancer Res 15(21):6683–6693, 2009. 10.1158/1078-0432.CCR-08-2801) resulted in partially discordant results in 22 tumors (27.5%). Local deep bisulfite sequencing (LDBS) revealed that CpGs located in exon 1 are suited best to discriminate methylated from unmethylated samples. Based on LDBS data, we propose an optimized MSP primer pair with 83% and 85% concordance to pyrosequencing and LDBS data. A hitherto neglected region upstream of the TSS, with an overall higher methylation compared to exon 1 and intron 1 of MGMT, is also able to discriminate the methylation status.

Conclusion

Our integrated analysis allows to evaluate and redefine co-methylation domains within the MGMT promoter and to rationalize the practical impact on assays used in daily routine diagnostics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01244-4.

BDNF methylation and mRNA expression in brain and blood of completed suicides in Slovenia

BACKGROUND

Suicide is a major public health problem. Worldwide, around 800000 people die by suicide every year. Suicide is a multifactorial disorder, with numerous environmental and genetic risk factors involved. Among the candidate genes, changes in the BDNF locus at the gene, epigenetic, mRNA, and protein expression levels have been implicated in psychiatric disorders, including suicidal behavior and completed suicides.

AIM

To investigate changes in BDNF methylation and expression of four alternative BDNF transcripts for association with completed suicide.

METHODS

This case-control study included 42 unrelated male Caucasian subjects, where 20 were control subjects who died following acute cardiac arrest, and 22 were suicide victims who died by hanging. DNA and RNA were extracted from brain tissue (Brodmann area 9 and hippocampus) and from blood. DNA methylation and mRNA expression levels were determined by targeted bisulfite next-generation sequencing and reverse-transcription quantitative PCR. Statistical analysis was done by use of two-tailed Student’s t tests for two independent samples, and the Benjamini-Hochberg procedure was implemented for correction for multiple comparisons.

RESULTS

In DNA from brain tissue, there were no significant differences in BDNF methylation between the study groups. However, data showed significantly reduced DNA methylation of the BDNF region upstream of exon I in blood samples of suicide victims compared to the controls (5.67 ± 0.57 vs 6.83 ± 0.64, P_corr = 0.01). In Brodmann area 9 of the brain of the suicide victims but not in their hippocampus, there was higher expression of BDNF transcript I-IX (NM_170731.4) compared to the controls (0.077 ± 0.024 vs 0.05 ± 0.013, P = 0.042). In blood, expression analysis for the BDNF transcripts was not feasible due to extensive RNA degradation.

CONCLUSION

Despite the limitations of the study, the obtained data further support a role for BDNF in suicidality. However, it should be noted that suicidal behavior is a multifactorial disorder with numerous environmental and genetic risk factors involved.

Epigenetic upregulation of TET2 is an independent poor prognostic factor for intrahepatic cholangiocarcinoma

Mutations in IDH1/2 and the epigenetic silencing of TET2 occur in leukaemia or glioma in a mutually exclusive manner. Although intrahepatic cholangiocarcinoma (iCCA) may harbour IDH1/2 mutations, the contribution of TET2 to carcinogenesis remains unknown. In the present study, the expression and promoter methylation of TET2 were investigated in iCCA. The expression of TET2 was assessed in 52 cases of iCCA (small-duct type, n = 33; large-duct type, n = 19) by quantitative PCR, immunohistochemistry (IHC) and a sequencing-based methylation assay, and its relationships with clinicopathological features and alterations in cancer-related genes (e.g., KRAS and IDH1) were investigated. In contrast to non-neoplastic bile ducts, which were negative for TET2 on IHC, 42 cases (81%) of iCCA showed the nuclear overexpression of TET2. Based on IHC scores (area × intensity), these cases were classified as TET2-high (n = 25) and TET2-low (n = 27). The histological type, tumour size, lymph node metastasis and frequency of mutations in cancer-related genes did not significantly differ between the two groups. Overall and recurrence-free survival were significantly worse in patients with TET2-high iCCA than in those with TET2-low iCCA. A multivariate analysis identified the high expression of TET2 as an independent prognostic factor (HR = 2.94; p = 0.007). The degree of methylation at two promoter CpG sites was significantly less in TET2-high iCCA than in TET2-low iCCA or non-cancer tissue. In conclusion, in contrast to other IDH-related neoplasms, TET2 overexpression is common in iCCA of both subtypes, and its high expression, potentially induced by promoter hypomethylation, is an independent poor prognostic factor.

DNA hypomethylation leads to cGAS-induced autoinflammation in the epidermis

DNA methylation is a fundamental epigenetic modification, important across biological processes. The maintenance methyltransferase DNMT1 is essential for lineage differentiation during development, but its functions in tissue homeostasis are incompletely understood. We show that epidermis-specific DNMT1 deletion severely disrupts epidermal structure and homeostasis, initiating a massive innate immune response and infiltration of immune cells. Mechanistically, DNA hypomethylation in keratinocytes triggered transposon derepression, mitotic defects, and formation of micronuclei. DNA release into the cytosol of DNMT1-deficient keratinocytes activated signaling through cGAS and STING, thus triggering inflammation. Our findings show that disruption of a key epigenetic mark directly impacts immune and tissue homeostasis, and potentially impacts our understanding of autoinflammatory diseases and cancer immunotherapy.

Methylation patterns of Tf2 retrotransposons linked to rapid adaptive stress response in the brown planthopper (Nilaparvata lugens)

Transposable elements (TEs) exhibit vast diversity across insect orders and are one of the major factors driving insect evolution and speciation. Presence of TEs can be both beneficial and deleterious to their host. While it is well-established that TEs impact life-history traits, adaptations and survivability of insects under hostile environments, the influence of the ecological niche on TE-landscape remains unclear. Here, we analysed the dynamics of Tf2 retrotransposons in the brown planthopper (BPH), under environmental fluctuations. BPH, a major pest of rice, is found in almost all rice-growing ecosystems. We believe genome plasticity, attributed to TEs, has allowed BPH to adapt and colonise novel ecological niches. Our study revealed bimodal age-distribution for Tf2 elements in BPH, indicating the occurrence of two major transpositional events in its evolutionary history and their contribution in shaping BPH genome. While TEs can provide genome flexibility and facilitate adaptations, they impose massive load on the genome. Hence, we investigated the involvement of methylation in modulating transposition in BPH. We performed comparative analyses of the methylation patterns of Tf2 elements in BPH feeding on resistant- and susceptible-rice varieties, and also under pesticide stress, across different life-stages. Results confirmed that methylation, particularly in non-CG context, is involved in TE regulation and dynamics under stress. Furthermore, we observed differential methylation for BPH adults and nymphs, emphasising the importance of screening juvenile life-stages in understanding adaptive-stress-responses in insects. Collectively, this study enhances our understanding of the role of transposons in influencing the evolutionary trajectory and survival strategies of BPH across generations.

Influence of tobacco cigarette heavy smoking on DNA methylation patterns and transcription levels of MAPK8IP3, GAA, ANXA2, PRRC2A, and PDE11A genes in human spermatozoa

Background

Tobacco smoking is considered as one of the lifestyles factors that influence the sperm DNA methylation and global sperm DNA methylation and that may affect the sperm phenotype. This study was performed to investigate whether tobacco cigarette heavy smoking influences sperm DNA methylation patterns and semen parameters and to determine whether there is an alteration in the transcription level of MAPK8IP3, GAA, ANXA2, PRRC2A, and PDE11A genes in heavy smokers compared to non-smokers. Thirty samples were subjected to 450K arrays as a screening study to assess the variation in sperm DNA methylation levels between heavy smokers and non-smokers. Five CpG sites have the highest difference in methylation levels (cg07869343, cg05813498, cg09785377, cg06833981, and cg02745784), which are located in the MAPK8IP3, GAA, ANXA2, PRRC2A, and PDE11A genes, respectively, and were selected for further analysis using deep bisulfite sequencing in 280 independent samples (120 proven non-smokers and 160 heavy smokers) with a mean age of 33.8 ± 8.4 years. The global sperm DNA methylation, sperm DNA fragmentation, and chromatin non-condensation were evaluated also.

Results

A significant increase was found in the methylation level at seven, three, and seventeen CpGs within the GAA, ANXA2, and MAPK8IP3 genes amplicon, respectively (P< 0.01) in heavy smokers compared to non-smokers. Additionally, a significant increase was found in the methylation levels at all CpGs within PRRC2A and PDE11A gene amplicon (P< 0.01). A significant increase was found in the level of sperm chromatin non-condensation, DNA fragmentation, and global DNA methylation (P < 0.001) in heavy smokers compared to non-smokers.

Conclusion

These results indicate that tobacco cigarette smoking can alter the DNA methylation level at several CpGs, the status of global DNA methylation, and transcription level of the following genes “MAPK8IP3, GAA, ANXA2, PRRC2A, and PDE11A” in human spermatozoa. These findings may affect negatively semen parameters and men’s fertility.

Locus-Specific DNA Methylation Editing in Melanoma Cell Lines Using a CRISPR-Based System

Simple Summary

DNA methylation is an important modification of the genome that is implicated in the pathogenesis of numerous human diseases, including cancer. DNA methylation changes can alter the expression of critical genes, predisposing to disease progression. Existing techniques that can modify DNA methylation to investigate disease etiology are severely limited with regard to specificity, which means that establishing a causal link between DNA methylation changes and disease progression is difficult. The advent of CRISPR-based technologies has provided a powerful tool for more specific editing of DNA methylation. Here, we describe a comprehensive protocol for the design and application of a CRISPR-dCas9-based tool for editing DNA methylation at a target locus in human melanoma cell lines alongside protocols for downstream techniques used to evaluate subsequent methylation and gene expression changes in methylation-edited cells. Furthermore, we demonstrate highly efficacious methylation and demethylation of the EBF3 promoter across a panel of melanoma cell lines.

Abstract

DNA methylation is a key epigenetic modification implicated in the pathogenesis of numerous human diseases, including cancer development and metastasis. Gene promoter methylation changes are widely associated with transcriptional deregulation and disease progression. The advent of CRISPR-based technologies has provided a powerful toolkit for locus-specific manipulation of the epigenome. Here, we describe a comprehensive global workflow for the design and application of a dCas9-SunTag-based tool for editing the DNA methylation locus in human melanoma cells alongside protocols for downstream techniques used to evaluate subsequent methylation and gene expression changes in methylation-edited cells. Using transient system delivery, we demonstrate both highly efficacious methylation and demethylation of the EBF3 promoter, which is a putative epigenetic driver of melanoma metastasis, achieving up to a 304.00% gain of methylation and 99.99% relative demethylation, respectively. Furthermore, we employ a novel, targeted screening approach to confirm the minimal off-target activity and high on-target specificity of our designed guide RNA within our target locus.

Ten Years of EWAS

Epigenome-wide association study (EWAS) has been applied to analyze DNA methylation variation in complex diseases for a decade, and epigenome as a research target has gradually become a hot topic of current studies. The DNA methylation microarrays, next-generation, and third-generation sequencing technologies have prepared a high-quality platform for EWAS. Here, the progress of EWAS research is reviewed, its contributions to clinical applications, and mainly describe the achievements of four typical diseases. Finally, the challenges encountered by EWAS and make bold predictions for its future development are presented.

Identification of tissue-specific and common methylation quantitative trait loci in healthy individuals using MAGAR

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00415-6.

Epigenomic and transcriptional profiling identifies impaired glyoxylate detoxification in NAFLD as a risk factor for hyperoxaluria

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.

A DNA methylation-based liquid biopsy for triple-negative breast cancer

Here, we present a next-generation sequencing (NGS) methylation-based blood test called methylation DETEction of Circulating Tumour DNA (mDETECT) designed for the optimal detection and monitoring of metastatic triple-negative breast cancer (TNBC). Based on a highly multiplexed targeted sequencing approach, this assay incorporates features that offer superior performance and included 53 amplicons from 47 regions. Analysis of a previously characterised cohort of women with metastatic TNBC with limited quantities of plasma (<2 ml) produced an AUC of 0.92 for detection of a tumour with a sensitivity of 76% for a specificity of 100%. mDETECT_TNBC was quantitative and showed superior performance to an NGS TP53 mutation-based test carried out on the same patients and to the conventional CA15-3 biomarker. mDETECT also functioned well in serum samples from metastatic TNBC patients where it produced an AUC of 0.97 for detection of a tumour with a sensitivity of 93% for a specificity of 100%. An assay for BRCA1 promoter methylation was also incorporated into the mDETECT assay and functioned well but its clinical significance is currently unclear. Clonal Hematopoiesis of Indeterminate Potential was investigated as a source of background in control subjects but was not seen to be significant, though a link to adiposity may be relevant. The mDETECT_TNBC assay is a liquid biopsy able to quantitatively detect all TNBC cancers and has the potential to improve the management of patients with this disease.

Identification and validation of DNA methylation changes in pre-eclampsia

INTRODUCTION

Pre-eclampsia (PE) is a dangerous placental condition that can lead to premature labour, seizures and death of mother and infant. Several studies have identified altered placental DNA methylation in PE; however, there is widespread inconsistency between studies and most findings have not been replicated. This study aimed to identify and validate consistent differences in methylation across multiple PE cohorts.

METHODS

Seven publicly available 450K methylation array datasets were analysed to identify consistent differentially methylated positions (DMPs) in PE. DMPs were identified based on methylation difference (≥10%) and significance (p-value ≤ 1 × 10^-7). Targeted deep bisulfite sequencing was then performed to validate a subset of DMPs in an additional independent PE cohort.

RESULTS

Stringent analysis of the seven 450K datasets identified 25 DMPs (associated with 11 genes) in only one dataset. Using more relaxed criteria confirmed 19 of the stringent 25 DMPs in at least four of the remaining six datasets. Targeted deep bisulfite sequencing of eight DMPs (associated with three genes; CMIP, ST3GAL1 and DAPK3) in an independent PE cohort validated two DMPs in the CMIP gene. Seven additional CpG sites in CMIP were found to be significantly differentially methylated in PE.

DISCUSSION

The identification and validation of significant differential methylation in CMIP suggests that the altered DNA methylation of this gene may be associated with the pathogenesis of PE, and may have the potential to serve as diagnostic biomarkers for this dangerous condition of pregnancy.

Bioinformatic tools for DNA methylation and histone modification: A survey

Epigenetic inheritance occurs due to different mechanisms such as chromatin and histone modifications, DNA methylation and processes mediated by non-coding RNAs. It leads to changes in gene expressions and the emergence of new traits in different organisms in many diseases such as cancer. Recent advances in experimental methods led to the identification of epigenetic target sites in various organisms. Computational approaches have enabled us to analyze mass data produced by these methods. Next-generation sequencing (NGS) methods have been broadly used to identify these target sites and their patterns. By using these patterns, the emergence of diseases could be prognosticated. In this study, target site prediction tools for two major epigenetic mechanisms comprising histone modification and DNA methylation are reviewed. Publicly accessible databases are reviewed as well. Some suggestions regarding the state-of-the-art methods and databases have been made, including examining patterns of epigenetic changes that are important in epigenotypes detection.

An aging and p53 related marker: HOXA5 promoter methylation negatively correlates with mRNA and protein expression in old age

The process of aging has been associated with differential patterns of DNA methylation which relate to changes in gene expression. Hence, we aimed to identify genes with significant age-related methylation differences and to study their mRNA and protein expression profile. Genome-wide DNA methylation analysis was performed with the Illumina Infinium Methylation EPIC BeadChip Microarray on bisulfite-converted DNA prepared from monocytes derived from young (average age: 23.8 yrs) and old (average age: 81.5 yrs) volunteers that are separated by at least 50 years of age difference, n=4, respectively. Differentially methylated CpG sites were assigned to the associated genes and validated by deep sequencing analysis (n=20). Demonstrating an age-associated significant increase of methylation in the promoter region (p=1x10^-8), Homeobox A5 (HOXA5), also known to activate p53, emerged as an interesting candidate for further expression analyses by Realtime PCR, ELISA and Western Blot Analysis (n=30, respectively). Consistent with its hypermethylation we observed significant HOXA5 mRNA downregulation (p=0.0053) correlating with significant p53 downregulation (p=0.0431) in the old cohort. Moreover, we observed a significant change in HOXA5 protein expression (p=3x10^-5) negatively correlating with age and promoter methylation thus qualifying HOXA5 for an eligible p53-related aging marker.

Targeted De-Methylation of the FOXP3-TSDR Is Sufficient to Induce Physiological FOXP3 Expression but Not a Functional Treg Phenotype

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.

Micro-RNA-125a mediates the effects of hypomethylating agents in chronic myelomonocytic leukemia

BACKGROUND

Chronic myelomonocytic leukemia (CMML) is an aggressive hematopoietic malignancy that arises from hematopoietic stem and progenitor cells (HSPCs). Patients with CMML are frequently treated with epigenetic therapeutic approaches, in particular the hypomethylating agents (HMAs), azacitidine (Aza) and decitabine (Dec). Although HMAs are believed to mediate their efficacy via re-expression of hypermethylated tumor suppressors, knowledge about relevant HMA targets is scarce. As silencing of tumor-suppressive micro-RNAs (miRs) by promoter hypermethylation is a crucial step in malignant transformation, we asked for a role of miRs in HMA efficacy in CMML.

RESULTS

Initially, we performed genome-wide miR-expression profiling in a Kras^G12D-induced CMML mouse model. Selected candidates with prominently decreased expression were validated by qPCR in CMML mice and human CMML patients. These experiments revealed the consistent decrease in miR-125a, a miR with previously described tumor-suppressive function in myeloid neoplasias. Furthermore, we show that miR-125a downregulation is caused by hypermethylation of its upstream region and can be reversed by HMA treatment. By employing both lentiviral and CRISPR/Cas9-based miR-125a modification, we demonstrate that HMA-induced miR-125a upregulation indeed contributes to mediating the anti-leukemic effects of these drugs. These data were validated in a clinical context, as miR-125a expression increased after HMA treatment in CMML patients, a phenomenon that was particularly pronounced in cases showing clinical response to these drugs.

CONCLUSIONS

Taken together, we report decreased expression of miR-125a in CMML and delineate its relevance as mediator of HMA efficacy within this neoplasia.

Hairpin-bisulfite sequencing of cells exposed to decitabine documents the process of DNA demethylation

Although the mechanism of DNA demethylating drugs has been understood for many years, the direct effect of these drugs on methylation of the complementary strands of DNA has not been formally demonstrated. By using hairpin-bisulphite sequencing, we describe the kinetics and pattern of DNA methylation following treatment of cells by the DNA methyltransferase 1 (DNMT1) inhibitor, decitabine. As expected, we demonstrate complete loss of methylation on the daughter strand following S-phase in selected densely methylated genes in synchronized Jurkat cells. Thereafter, cells showed a heterogeneous pattern of methylation reflecting replication of the unmethylated strand and restoration of methylation.

Queuine links translational control in eukaryotes to a micronutrient from bacteria

In eukaryotes, the wobble position of tRNA with a GUN anticodon is modified to the 7-deaza-guanosine derivative queuosine (Q34), but the original source of Q is bacterial, since Q is synthesized by eubacteria and salvaged by eukaryotes for incorporation into tRNA. Q34 modification stimulates Dnmt2/Pmt1-dependent C38 methylation (m⁵C38) in the tRNA^Asp anticodon loop in Schizosaccharomyces pombe. Here, we show by ribosome profiling in S. pombe that Q modification enhances the translational speed of the C-ending codons for aspartate (GAC) and histidine (CAC) and reduces that of U-ending codons for asparagine (AAU) and tyrosine (UAU), thus equilibrating the genome-wide translation of synonymous Q codons. Furthermore, Q prevents translation errors by suppressing second-position misreading of the glycine codon GGC, but not of wobble misreading. The absence of Q causes reduced translation of mRNAs involved in mitochondrial functions, and accordingly, lack of Q modification causes a mitochondrial defect in S. pombe. We also show that Q-dependent stimulation of Dnmt2 is conserved in mice. Our findings reveal a direct mechanism for the regulation of translational speed and fidelity in eukaryotes by a nutrient originating from bacteria.

IGF2 mRNA Binding Protein 2 Transgenic Mice Are More Prone to Develop a Ductular Reaction and to Progress Toward Cirrhosis

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.

Two are better than one: HPoxBS - hairpin oxidative bisulfite sequencing

The controlled and stepwise oxidation of 5mC to 5hmC, 5fC and 5caC by Tet enzymes is influencing the chemical and biological properties of cytosine. Besides direct effects on gene regulation, oxidised forms influence the dynamics of demethylation and re-methylation processes. So far, no combined methods exist which allow to precisely determine the strand specific localisation of cytosine modifications along with their CpG symmetric distribution. Here we describe a comprehensive protocol combining conventional hairpin bisulfite with oxidative bisulfite sequencing (HPoxBS) to determine the strand specific distribution of 5mC and 5hmC at base resolution. We apply this method to analyse the contribution of local oxidative effects on DNA demethylation in mouse ES cells. Our method includes the HPoxBS workflow and subsequent data analysis using our developed software tools. Besides a precise estimation and display of strand specific 5mC and 5hmC levels at base resolution we apply the data to predict region specific activities of Dnmt and Tet enzymes. Our experimental and computational workflow provides a precise double strand display of 5mC and 5hmC modifications at single base resolution. Based on our data we predict region specific Tet and Dnmt enzyme efficiencies shaping the distinct locus levels and patterns of 5hmC and 5mC.

Locus-Specific DNA Methylation Analysis by Targeted Deep Bisulfite Sequencing

DNA methylation, i.e., the methylation of cytosine at carbon atom C5, has an important role in the regulation of gene expression. The methylation status of each cytosine in a specific genomic region can be determined by targeted deep bisulfite sequencing at single-molecule resolution. Here we describe the design of PCR primers that are used to amplify specific sequences from bisulfite-converted DNA, the preparation of sequencing libraries, the sequencing of these libraries on the MiSeq system, as well as the analysis of the sequence reads. Using appropriate software tools such as amplikyzer2, it is easy to analyze complex multiplexed samples with several regions of interest, to determine the mean methylation values of all CpG dinucleotides in a region or of each CpG dinucleotide across all or selected reads, and to compare these values between different samples and between different alleles within a sample.

Division of labour: tRNA methylation by the NSun2 tRNA methyltransferases Trm4a and Trm4b in fission yeast

Enzymes of the cytosine-5 RNA methyltransferase Trm4/NSun2 family methylate tRNAs at C48 and C49 in multiple tRNAs, as well as C34 and C40 in selected tRNAs. In contrast to most other organisms, fission yeast Schizosaccharomyces pombe carries two Trm4/NSun2 homologs, Trm4a (SPAC17D4.04) and Trm4b (SPAC23C4.17). Here, we have employed tRNA methylome analysis to determine the dependence of cytosine-5 methylation (m⁵C) tRNA methylation in vivo on the two enzymes. Remarkably, Trm4a is responsible for all C48 methylation, which lies in the tRNA variable loop, as well as for C34 in tRNA^Leu_CAA and tRNA^Pro_CGG, which are at the anticodon wobble position. Conversely, Trm4b methylates C49 and C50, which both lie in the TΨC-stem. Thus, S. pombe show an unusual separation of activities of the NSun2/Trm4 enzymes that are united in a single enzyme in other eukaryotes like humans, mice and Saccharomyces cerevisiae. Furthermore, in vitro activity assays showed that Trm4a displays intron-dependent methylation of C34, whereas Trm4b activity is independent of the intron. The absence of Trm4a, but not Trm4b, causes a mild resistance of S. pombe to calcium chloride.

DNA methylation of the promoter region of bnip3 and bnip3l genes induced by metabolic programming

BACKGROUND

Environmental changes of biotic or abiotic nature during critical periods of early development may exert a profound influence on physiological functions later in life. This process, named developmental programming can also be driven through parental nutrition. At molecular level, epigenetic modifications are the most likely candidate for persistent modulation of genes expression in later life.

RESULTS

In order to investigate epigenetic modifications induced by programming in rainbow trout, we focused on bnip3 and bnip3l paralogous genes known to be sensitive to environmental changes but also regulated by epigenetic modifications. Two specific stimuli were used: (i) early acute hypoxia applied at embryo stage and (ii) broodstock and fry methionine deficient diet, considering methionine as one of the main methyl-group donor needed for DNA methylation. We observed a programming effect of hypoxia with an increase of bnip3a and the four paralogs of bnip3l expression level in fry. In addition, parental methionine nutrition was correlated to bnip3a and bnip3lb1 expression showing evidence for early fry programming. We highlighted that both stimuli modified DNA methylation levels at some specific loci of bnip3a and bnip3lb1.

CONCLUSION

Overall, these data demonstrate that methionine level and hypoxia stimulus can be of critical importance in metabolic programming. Both stimuli affected DNA methylation of specific loci, among them, an interesting CpG site have been identified, namely - 884 bp site of bnip3a, and may be positively related with mRNA levels.

Decreased DNA methylation of a CpG site in the HBAP1 gene in plasma DNA from pregnant women

Objective

The objective of this study is to identify potential CpG site(s) or DNA methylation pattern(s) in the pseudo α-globin 1 gene (HBAP1 gene), the gene which locates in α-thalassemia-1 deletion mutation, to differentiate plasma DNA between pregnant and non-pregnant women.

Method

DNA methylation profiles of placenta and peripheral blood from the MethBase database were compared to screen differentially methylated regions. This region was confirmed the differential by methylation-sensitive high resolution melt (MS-HRM) analysis. The differential region was used to compare DNA methylation profile of plasma DNA between pregnant and non-pregnant women by bisulfite amplicon sequencing in three levels: overall, individual CpG sites and individual molecules (DNA methylation patterns).

Result

Using MethBase data, four consecutive CpG sites in the HBAP1 gene were identified as regions of differential DNA methylation between placenta and peripheral blood. The confirmation by MS-HRM showed the differential DNA methylation profile between the placenta and plasma from non-pregnant women. The comparison of DNA methylation profiles between the plasma of pregnant and non-pregnant women showed that, in the overall levels of the four CpG sites, DNA methylation of pregnant women was detected at lower levels than non-pregnant women. In the individual CpG site level, only the second CpG site showed differential DNA methylation between the groups. In the DNA methylation pattern level, there was no strongly significant differences in DNA methylation patterns between the pregnant and non-pregnant groups.

Conclusion

Our result demonstrated that, in the plasma from pregnant women, only one of the four CpG sites displays a decrease in DNA methylation compared with non-pregnant women. It indicates that this CpG site might be useful for determining the presence or absence of fetal wild-type α-globin gene cluster allele in maternal plasma.

CD137+CD154- Expression As a Regulatory T Cell (Treg)-Specific Activation Signature for Identification and Sorting of Stable Human Tregs from In Vitro Expansion Cultures

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 alterations in iPSC- and hESC-derived neurons: potential implications for neurological disease modeling

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.

Impact of cigarette-smoking on sperm DNA methylation and its effect on sperm parameters

DNA methylation is an epigenetic modification of the genome. The purpose of this study was to determine the influence of cigarette-smoking on sperm DNA methylation from a genomewide survey of sperm samples and to ascertain its effect on sperm parameters. Twenty-eight sperm DNA samples (from 14 fertile smokers as a case study and 14 proven fertile nonsmokers as controls) were subjected to Infinium 450K BeadChip arrays to identify the changes in the DNA methylation level between the two groups. Then, deep bisulphite sequencing was used to validate five CpGs on 78 samples. The results from the Infinium 450K found that only 11 CpGs showed a significant difference in DNA methylation between the case and the control groups. Five CpGs of the eleven (cg00648582, cg0932376, cg19169023, cg23841288 and cg27391564) underwent deep bisulphite sequencing where cg00648582, related to the PGAM5 gene, and the cg23841288 CpGs, related to the PTPRN2 gene amplicons, showed a significant increase in their DNA methylation level in more than one CpG in the case group. In contrast, a significant decrease was found at cg19169023 and at its various neighbouring CpGs in the TYRO3 gene-related amplicons. Furthermore, this study demonstrated a significant correlation between the variation in sperm DNA methylation level and standard sperm parameters in the case group.

Hairpin Bisulfite Sequencing: Synchronous Methylation Analysis on Complementary DNA Strands of Individual Chromosomes

The accurate and quantitative detection of 5-methylcytosine is of great importance in the field of epigenetics. The method of choice is usually bisulfite sequencing because of the high resolution and the possibility to combine it with next generation sequencing. Nevertheless, also this method has its limitations. Following the bisulfite treatment DNA strands are no longer complementary such that in a subsequent PCR amplification the DNA methylation patterns information of only one of the two DNA strand is preserved. Several years ago Hairpin Bisulfite sequencing was developed as a method to obtain the pattern information on complementary DNA strands. The method requires fragmentation (usually by enzymatic cleavage) of genomic DNA followed by a covalent linking of both DNA strands through ligation of a short DNA hairpin oligonucleotide to both strands. The ligated covalently linked dsDNA products are then subjected to a conventional bisulfite treatment during which all unmodified cytosines are converted to uracils. During the treatment the DNA is denatured forming noncomplementary ssDNA circles. These circles serve as a template for a locus specific PCR to amplify chromosomal patterns of the region of interest. As a result one ends up with a linearized product, which contains the methylation information of both complementary DNA strands.

Nucleosome Occupancy and Methylome Sequencing (NOMe-seq)

Various methodologies are available to interrogate specific components of epigenetic mechanisms such as DNA methylation or nucleosome occupancy at both the locus-specific and the genome-wide level. It has become increasingly clear, however, that comprehension of the functional interactions between epigenetic mechanisms is critical for understanding how cellular transcription programs are regulated or deregulated during normal and disease development. The Nucleosome Occupancy and Methylome sequencing (NOMe-seq) assay allows us to directly measure the relationship between DNA methylation and nucleosome occupancy by taking advantage of the methyltransferase M.CviPI, which methylates unprotected GpC dinucleotides to create a footprint of chromatin accessibility. This assay generates dual nucleosome occupancy and DNA methylation information at a single-DNA molecule resolution using as little as 200,000 cells and in as short as 15 min reaction time. DNA methylation levels and nucleosome occupancy status of genomic regions of interest can be subsequently interrogated by cloning PCR-amplified bisulfite DNA and sequencing individual clones. Alternatively, NOMe-seq can be combined with next-generation sequencing in order to generate an integrated global map of DNA methylation and nucleosome occupancy, which allows for comprehensive examination as to how these epigenetic components correlate with each other.

Epigenetic Basis of Circadian Rhythm Disruption in Cancer

Self-sustained and synchronized to environmental stimuli, circadian clocks are under genetic and epigenetic regulation. Recent findings have greatly increased our understanding of epigenetic plasticity governed by circadian clock. Thus, the link between circadian clock and epigenetic machinery is reciprocal. Circadian clock can affect epigenetic features including genomic DNA methylation, noncoding RNA, mainly miRNA expression, and histone modifications resulted in their 24-h rhythms. Concomitantly, these epigenetic events can directly modulate cyclic system of transcription and translation of core circadian genes and indirectly clock output genes. Significant findings interlocking circadian clock, epigenetics, and cancer have been revealed, particularly in breast, colorectal, and blood cancers. Aberrant methylation of circadian gene promoter regions and miRNA expression affected circadian gene expression, together with 24-h expression oscillation pace have been frequently observed.

Gut memories do not fade: epigenetic regulation of lasting gut homing receptor expression in CD4+ memory T cells

The concept of a "topographical memory" in lymphocytes implies a stable expression of homing receptors mediating trafficking of lymphocytes back to the tissue of initial activation. However, a significant plasticity of the gut-homing receptor α₄β₇ was found in CD8⁺ T cells, questioning the concept. We now demonstrate that α₄β₇ expression in murine CD4⁺ memory T cells is, in contrast, imprinted and remains stable in the absence of the inducing factor retinoic acid (RA) or other stimuli from mucosal environments. Repetitive rounds of RA treatment enhanced the stability of de novo induced α₄β₇. A novel enhancer element in the murine Itga4 locus was identified that showed, correlating to stability, selective DNA demethylation in mucosa-seeking memory cells and methylation-dependent transcriptional activity in a reporter gene assay. This implies that epigenetic mechanisms contribute to the stabilization of α₄β₇ expression. Analogous DNA methylation patterns could be observed in the human ITGA4 locus, suggesting that its epigenetic regulation is conserved between mice and men. These data prove that mucosa-specific homing mediated by α₄β₇ is imprinted in CD4⁺ memory T cells, reinstating the validity of the concept of "topographical memory" for mucosal tissues, and imply a critical role of epigenetic mechanisms.

Aberrations in sperm DNA methylation patterns of males suffering from reduced fecundity

The purpose of this study was to evaluate the aberrations in sperm DNA methylation patterns of males suffering from reduced fecundity. A total of 108 males (65 males suffering from reduced fecundity as cases and 43 proven fertile males as a control) were included in the study. Thirty samples were subjected to 450K arrays as a screening phase, and then, three CpG sites located in the following genes: TYRO3, CGβ and FAM189A1 were selected to validate on 78 samples using deep bisulphite sequencing. A significant difference in the methylation level was found between cases and controls at all CpGs in TYRO3 gene-related amplicon (CpG1, p ≤ .003, CpG2, p ≤ .0001, CpG3, p ≤ .003 and CpG4, p ≤ .030) and CpG1 in CGβ gene-related amplicon (p ≤ .0001). Besides, a significant difference was found at two CpGs (CpG1, p ≤ .004 and CpG2, p ≤ .002) tested in the FAM189A1 gene-related amplicon. A significant correlation was found between the methylation level at CpG1 in the FAM189A1 gene and the different types of sperm motility. In conclusion, an alteration in the methylation levels of sperm DNA from males with reduced fecundity was showed. In addition, a relationship between variations in the methylation level of these CpGs and sperm motility has been observed.

CpG location and methylation level are crucial factors for the early detection of oral squamous cell carcinoma in brushing samples using bisulfite sequencing of a 13-gene panel

Background

Oral squamous cell carcinoma (OSCC) is usually diagnosed at an advanced stage and is commonly preceded by oral premalignant lesions. The mortality rates have remained unchanged (50% within 5 years after diagnosis), and it is related to tobacco smoking and alcohol intake. Novel molecular markers for early diagnosis are urgently needed. The purpose of this study was to evaluate the diagnostic value of methylation level in a set of 18 genes by bisulfite next-generation sequencing.

Methods

With minimally invasive oral brushing, 28 consecutive OSCC, one squamous cell carcinoma with sarcomatoid features, six high-grade squamous intraepithelial lesions (HGSIL), 30 normal contralateral mucosa from the same patients, and 65 healthy donors were evaluated for DNA methylation analyzing 18 target genes by quantitative bisulfite next-generation sequencing. We further evaluated an independent cohort (validation dataset) made of 20 normal donors, one oral fibroma, 14 oral lichen planus (OLP), three proliferative verrucous leukoplakia (PVL), and two OSCC.

Results

Comparing OSCC with normal healthy donors and contralateral mucosa in 355 CpGs, we identified the following epigenetically altered genes: ZAP70, ITGA4, KIF1A, PARP15, EPHX3, NTM, LRRTM1, FLI1, MIR193, LINC00599, PAX1, and MIR137HG showing hypermethylation and MIR296, TERT, and GP1BB showing hypomethylation. The behavior of ZAP70, GP1BB, H19, EPHX3, and MIR193 fluctuated among different interrogated CpGs. The gap between normal and OSCC samples remained mostly the same (Kruskal-Wallis P values < 0.05), but the absolute values changed conspicuously. ROC curve analysis identified the most informative CpGs, and we correctly stratified OSCC and HGSIL from normal donors using a multiclass linear discriminant analysis in a 13-gene panel (AUC 0.981). Only the OSCC with sarcomatoid features was negative. Three contralateral mucosa were positive, a sign of a possible field cancerization. Among imprinted genes, only MIR296 showed loss of imprinting. DNMT1, TERC, and H19 together with the global methylation of long interspersed element 1 were unchanged. In the validation dataset, values over the threshold were detected in 2/2 OSCC, in 3/3 PVL, and in 2/14 OLP.

Conclusions

Our data highlight the importance of CpG location and correct estimation of DNA methylation level for highly accurate early diagnosis of OSCC.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0386-7) contains supplementary material, which is available to authorized users.

Differential methylation is associated with non-syndromic cleft lip and palate and contributes to penetrance effects

Non-syndromic cleft lip and/or palate (NSCLP) is a common congenital malformation with a multifactorial model of inheritance. Although several at-risk alleles have been identified, they do not completely explain the high heritability. We postulate that epigenetic factors as DNA methylation might contribute to this missing heritability. Using a Methylome-wide association study in a Brazilian cohort (67 NSCLP, 59 controls), we found 578 methylation variable positions (MVPs) that were significantly associated with NSCLP. MVPs were enriched in regulatory and active regions of the genome and in pathways already implicated in craniofacial development. In an independent UK cohort (171 NSCLP, 177 controls), we replicated 4 out of 11 tested MVPs. We demonstrated a significant positive correlation between blood and lip tissue DNA methylation, indicating blood as a suitable tissue for NSCLP methylation studies. Next, we quantified CDH1 promoter methylation levels in CDH1 mutation-positive families, including penetrants, non-penetrants or non-carriers for NSCLP. We found methylation levels to be significantly higher in the penetrant individuals. Taken together, our results demonstrated the association of methylation at specific genomic locations as contributing factors to both non-familial and familial NSCLP and altered DNA methylation may be a second hit contributing to penetrance.