Dynamic DNA methylation across diverse human cell lines and tissues

An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia

Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

Weaning Induces Stress-Dependent DNA Methylation and Transcriptional Changes in Piglet PBMCs

Changes to the epigenome, including those to DNA methylation, have been proposed as mechanisms by which stress can induce long-term physiological changes in livestock species. Pig weaning is associated with dietary and social stress, both of which elicit an immune response and changes to the hypothalamic–pituitary–adrenal (HPA) axis. While differential methylation following stress has been assessed in model organisms, it remains poorly understood how the pig methylome is altered by stressors in production settings. We quantified changes in CpG methylation and transcript abundance in piglet peripheral blood mononuclear cells (PBMCs) following weaning and also assessed differential patterns in pigs exhibiting high and low stress response as measured by cortisol concentration and lesion scores. Blood was collected from nine gilt piglets 24 h before and after weaning, and whole-genome bisulfite sequencing (WGBS) and RNA-sequencing were performed on six and nine animals, respectively, at both time points. We identified 2,674 differentially methylated regions (DMRs) that were enriched within promoters of genes associated with lymphocyte stimulation and transcriptional regulation. Stress groups displayed unique differential methylation and expression patterns associated with activation and suppression of T cell immunity in low and high stress animals, respectively. Differential methylation was strongly associated with differential expression; specifically, upregulated genes were enriched among hypomethylated genes. We observed post-weaning hypermethylation of the glucocorticoid receptor (NR3C1) promoter and a significant decrease in NR3C1 expression (n = 9, p = 6.1 × 10^–3). Our results indicate that weaning-associated stress elicits genome-wide methylation changes associated with differential gene expression, reduced T cell activation, and an altered HPA axis response.

Preclinical models of pancreatic ductal adenocarcinoma: challenges and opportunities in the era of precision medicine

Pancreatic ductal adenocarcinoma (PDAC) is an extremely lethal malignancy, with an average 5-year survival rate of 9% (Siegel RL, Miller KD, Jemal A. Ca Cancer J Clin. 2019;69(1):7-34). The steady increase in mortality rate indicates limited efficacy of the conventional regimen. The heterogeneity of PDAC calls for personalized treatment in clinical practice, which requires the construction of a preclinical system for generating patient-derived models. Currently, the lack of high-quality preclinical models results in ineffective translation of novel targeted therapeutics. This review summarizes applications of commonly used models, discusses major difficulties in PDAC model construction and provides recommendations for integrating workflows for precision medicine.

Gene-Specific DNA Methylation Linked to Postoperative Cognitive Dysfunction in Apolipoprotein E3 and E4 Mice

BACKGROUND

Geriatric surgical patients are at higher risk of developing postoperative neurocognitive disorders (NCD) than younger patients. The specific mechanisms underlying postoperative NCD remain unknown, but they have been linked to genetic risk factors, such as the presence of APOE4, compared to APOE3, and epigenetic modifications caused by exposure to anesthesia and surgery.

OBJECTIVE

To test the hypothesis that compared to E3 mice, E4 mice exhibit a more pronounced postoperative cognitive impairment associated with differential DNA methylation in brain regions linked to learning and memory.

METHODS

16-month-old humanized apolipoprotein-E targeted replacement mice bearing E3 or E4 were subjected to surgery (laparotomy) under general isoflurane anesthesia or sham. Postoperative behavioral testing and genome-wide DNA methylation were performed.

RESULTS

Exposure to surgery and anesthesia impaired cognition in aged E3, but not E4 mice, likely due to the already lower cognitive performance of E4 prior to surgery. Cognitive impairment in E3 mice was associated with hypermethylation of specific genes, including genes in the Ephrin pathway implicated in synaptic plasticity and learning in adults and has been linked to Alzheimer's disease. Other genes, such as the Scratch Family Transcriptional Repressor 2, were altered after surgery and anesthesia in both the E3 and E4 mice.

CONCLUSION

Our findings suggest that the neurocognitive and behavioral effects of surgery and anesthesia depend on baseline neurocognitive status and are associated with APOE isoform-dependent epigenetic modifications of specific genes and pathways involved in memory and learning.

Methylation Analysis in Monozygotic Twins With Treatment-Resistant Schizophrenia and Discordant Responses to Clozapine

Schizophrenia is a mental illness that involves both genetic and environmental factors. Clozapine, an atypical antipsychotic, is a well-established therapy for treatment-resistant schizophrenia. In this study, we focused on a set of monozygotic twins with treatment-resistant schizophrenia in which one twin effectively responded to clozapine treatment and the other did not. Our previous study generated neurons from induced pluripotent stem (iPS) cells derived from these patients and compared the transcriptome profiles between mock- and clozapine-treated neurons. In this study, we performed genome-wide DNA methylation profiling to investigate the mechanisms underlying gene expression changes. First, we extracted the differentially methylated sites from each twin based on statistical analysis. Then, we combined the DNA methylation profiling with transcriptome profiling from our previous RNA-seq data. Among the genes with altered methylation and expression, we found the different proportions of the genes related to neuronal and synaptic functions between the clozapine responder and non-responder (35.7 and 6.7%, respectively). This trend was observed even when the basal differences between the responder and non-responder was excluded. These results suggest that effective clozapine action may correct the abnormalities of neuronal and synapse functions in schizophrenia via changes in methylation.

DNA methylation and hydroxymethylation have distinct genome-wide profiles related to axonal regeneration

Alterations in environmentally sensitive epigenetic mechanisms (e.g., DNA methylation) influence axonal regeneration in the spinal cord following sharp injury. Conventional DNA methylation detection methods using sodium bisulphite treatment do not distinguish between methylated and hydroxymethylated forms of cytosine, meaning that past studies report a composite of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). To identify the distinct contributions of DNA methylation modifications to axonal regeneration, we collected spinal cord tissue after sharp injury from untreated adult F3 male rats with enhanced regeneration of injured spinal axons or controls, derived from folate- or water-treated F0 lineages, respectively. Genomic DNA was profiled for genome-wide 5hmC levels, revealing 658 differentially hydroxymethylated regions (DhMRs). Genomic profiling with whole genome bisulphite sequencing disclosed regeneration-related alterations in composite 5mC + 5hmC DNA methylation levels at 2,260 differentially methylated regions (DMRs). While pathway analyses revealed that differentially hydroxymethylated and methylated genes are linked to biologically relevant axon developmental pathways, only 22 genes harbour both DhMR and DMRs. Since these differential modifications were more than 60 kilobases on average away from each other, the large majority of differential hydroxymethylated and methylated regions are unique with distinct functions in the axonal regeneration phenotype. These data highlight the importance of distinguishing independent contributions of 5mC and 5hmC levels in the central nervous system, and denote discrete roles for DNA methylation modifications in spinal cord injury and regeneration in the context of transgenerational inheritance.

Integrative Analysis of Glucometabolic Traits, Adipose Tissue DNA Methylation, and Gene Expression Identifies Epigenetic Regulatory Mechanisms of Insulin Resistance and Obesity in African Americans

Decline in insulin sensitivity due to dysfunction of adipose tissue (AT) is one of the earliest pathogenic events in type 2 diabetes. We hypothesize that differential DNA methylation (DNAm) controls insulin sensitivity and obesity by modulating transcript expression in AT. Integrating AT DNAm profiles with transcript profile data measured in a cohort of 230 African Americans (AAs) from the African American Genetics of Metabolism and Expression cohort, we performed cis-expression quantitative trait methylation (cis-eQTM) analysis to identify epigenetic regulatory loci for glucometabolic trait-associated transcripts. We identified significantly associated cytosine-guanine dinucleotide regions for 82 transcripts (false discovery rate [FDR]-P < 0.05). The strongest eQTM locus was observed for the proopiomelanocortin (POMC; ρ = -0.632, P = 4.70 × 10^-27) gene. Epigenome-wide association studies (EWAS) further identified 155, 46, and 168 cytosine-guanine dinucleotide regions associated (FDR-P < 0.05) with the Matsuda index, S_I, and BMI, respectively. Intersection of EWAS, transcript level to trait association, and eQTM results, followed by causal inference test identified significant eQTM loci for 23 genes that were also associated with Matsuda index, S_I, and/or BMI in EWAS. These associated genes include FERMT3, ITGAM, ITGAX, and POMC In summary, applying an integrative multiomics approach, our study provides evidence for DNAm-mediated regulation of gene expression at both previously identified and novel loci for many key AT transcripts influencing insulin resistance and obesity.

Deconvolution of cellular subsets in human tissue based on targeted DNA methylation analysis at individual CpG sites

Background

The complex composition of different cell types within a tissue can be estimated by deconvolution of bulk gene expression profiles or with various single-cell sequencing approaches. Alternatively, DNA methylation (DNAm) profiles have been used to establish an atlas for multiple human tissues and cell types. DNAm is particularly suitable for deconvolution of cell types because each CG dinucleotide (CpG site) has only two states per DNA strand—methylated or non-methylated—and these epigenetic modifications are very consistent during cellular differentiation. So far, deconvolution of DNAm profiles implies complex signatures of many CpGs that are often measured by genome-wide analysis with Illumina BeadChip microarrays. In this study, we investigated if the characterization of cell types in tissue is also feasible with individual cell type-specific CpG sites, which can be addressed by targeted analysis, such as pyrosequencing.

Results

We compiled and curated 579 Illumina 450k BeadChip DNAm profiles of 14 different non-malignant human cell types. A training and validation strategy was applied to identify and test for cell type-specific CpGs. We initially focused on estimating the relative amount of fibroblasts using two CpGs that were either hypermethylated or hypomethylated in fibroblasts. The combination of these two DNAm levels into a “FibroScore” correlated with the state of fibrosis and was associated with overall survival in various types of cancer. Furthermore, we identified hypomethylated CpGs for leukocytes, endothelial cells, epithelial cells, hepatocytes, glia, neurons, fibroblasts, and induced pluripotent stem cells. The accuracy of this eight CpG signature was tested in additional BeadChip datasets of defined cell mixtures and the results were comparable to previously published signatures based on several thousand CpGs. Finally, we established and validated pyrosequencing assays for the relevant CpGs that can be utilized for classification and deconvolution of cell types.

Conclusion

This proof of concept study demonstrates that DNAm analysis at individual CpGs reflects the cellular composition of cellular mixtures and different tissues. Targeted analysis of these genomic regions facilitates robust methods for application in basic research and clinical settings.

TERT-Regulation and Roles in Cancer Formation

Telomerase reverse transcriptase (TERT) is a catalytic subunit of telomerase. Telomerase complex plays a key role in cancer formation by telomere dependent or independent mechanisms. Telomere maintenance mechanisms include complex TERT changes such as gene amplifications, TERT structural variants, TERT promoter germline and somatic mutations, TERT epigenetic changes, and alternative lengthening of telomere. All of them are cancer specific at tissue histotype and at single cell level. TERT expression is regulated in tumors via multiple genetic and epigenetic alterations which affect telomerase activity. Telomerase activity via TERT expression has an impact on telomere length and can be a useful marker in diagnosis and prognosis of various cancers and a new therapy approach. In this review we want to highlight the main roles of TERT in different mechanisms of cancer development and regulation.

Comparison of DNA Methylation Profiles of Hemostatic Genes between Liver Tissue and Peripheral Blood within Individuals

DNA methylation has become increasingly recognized in the etiology of complex diseases, including thrombotic disorders. Blood is often collected in epidemiological studies for genotyping and has recently also been used to examine DNA methylation in epigenome-wide association studies. DNA methylation patterns are often tissue-specific, thus, peripheral blood may not accurately reflect the methylation pattern in the tissue of relevance. Here, we collected paired liver and blood samples concurrently from 27 individuals undergoing liver surgery. We performed targeted bisulfite sequencing for a set of 35 hemostatic genes primarily expressed in liver to analyze DNA methylation levels of >10,000 cytosine-phosphate-guanine (CpG) dinucleotides. We evaluated whether DNA methylation in blood could serve as a proxy for DNA methylation in liver at individual CpGs. Approximately 30% of CpGs were nonvariable and were predominantly hypo- (<25%) or hypermethylated (>70%) in both tissues. While blood can serve as a proxy for liver at these CpGs, the low variability renders these unlikely to explain phenotypic differences. We therefore focused on CpG sites with variable methylation levels in liver. The level of blood-liver tissue correlation varied widely across these variable CpGs; moderate correlations (0.5 ≤ r < 0.75) were detected for 6% and strong correlations (r ≥ 0.75) for a further 4%. Our findings indicate that it is essential to study the concordance of DNA methylation between blood and liver at individual CpGs. This paired blood-liver dataset is intended as a resource to aid interpretation of blood-based DNA methylation results.

An epigenome-wide association study of Alzheimer's disease blood highlights robust DNA hypermethylation in the HOXB6 gene

A growing number of epigenome-wide association studies have demonstrated a role for DNA methylation in the brain in Alzheimer's disease. With the aim of exploring peripheral biomarker potential, we have examined DNA methylation patterns in whole blood collected from 284 individuals in the AddNeuroMed study, which included 89 nondemented controls, 86 patients with Alzheimer's disease, and 109 individuals with mild cognitive impairment, including 38 individuals who progressed to Alzheimer's disease within 1 year. We identified significant differentially methylated regions, including 12 adjacent hypermethylated probes in the HOXB6 gene in Alzheimer's disease, which we validated using pyrosequencing. Using weighted gene correlation network analysis, we identified comethylated modules of genes that were associated with key variables such as APOE genotype and diagnosis. In summary, this study represents the first large-scale epigenome-wide association study of Alzheimer's disease and mild cognitive impairment using blood. We highlight the differences in various loci and pathways in early disease, suggesting that these patterns relate to cognitive decline at an early stage.

Hypoxia and Oxygen-Sensing Signaling in Gene Regulation and Cancer Progression

Oxygen homeostasis regulation is the most fundamental cellular process for adjusting physiological oxygen variations, and its irregularity leads to various human diseases, including cancer. Hypoxia is closely associated with cancer development, and hypoxia/oxygen-sensing signaling plays critical roles in the modulation of cancer progression. The key molecules of the hypoxia/oxygen-sensing signaling include the transcriptional regulator hypoxia-inducible factor (HIF) which widely controls oxygen responsive genes, the central members of the 2-oxoglutarate (2-OG)-dependent dioxygenases, such as prolyl hydroxylase (PHD or EglN), and an E3 ubiquitin ligase component for HIF degeneration called von Hippel–Lindau (encoding protein pVHL). In this review, we summarize the current knowledge about the canonical hypoxia signaling, HIF transcription factors, and pVHL. In addition, the role of 2-OG-dependent enzymes, such as DNA/RNA-modifying enzymes, JmjC domain-containing enzymes, and prolyl hydroxylases, in gene regulation of cancer progression, is specifically reviewed. We also discuss the therapeutic advancement of targeting hypoxia and oxygen sensing pathways in cancer.

Genome-Wide Assessment of DNA Methylation in Chicken Cardiac Tissue Exposed to Different Incubation Temperatures and CO2 Levels

Temperature and CO₂ concentration during incubation have profound effects on broiler chick development, and numerous studies have identified significant effects on hatch heart weight (HW) as a result of differences in these parameters. Early life environment has also been shown to affect broiler performance later in life; it has thus been suggested that epigenetic mechanisms may mediate long-term physiological changes induced by environmental stimuli. DNA methylation is an epigenetic modification that can confer heritable changes in gene expression. Using reduced-representation bisulfite sequencing (RRBS), we assessed DNA methylation patterns in cardiac tissue of 84 broiler hatchlings incubated at two egg shell temperatures (EST; 37.8°C and 38.9°C) and three CO₂ concentrations (0.1%, 0.4%, and 0.8%) from day 8 of incubation onward. We assessed differential methylation between EST treatments and identified 2,175 differentially methylated (DM) CpGs (1,121 hypermethylated, 1,054 hypomethylated at 38.9° vs. 37.8°) in 269 gene promoters and 949 intragenic regions. DM genes (DMGs) were associated with heart developmental processes, including cardiomyocyte proliferation and differentiation. We identified enriched binding motifs among DM loci, including those for transcription factors associated with cell proliferation and heart development among hypomethylated CpGs that suggest increased binding ability at higher EST. We identified 9,823 DM CpGs between at least two CO₂ treatments, with the greatest difference observed between 0.8 and 0.1% CO₂ that disproportionately impacted genes involved in cardiac muscle development and response to low oxygen levels. Using HW measurements from the same chicks, we performed an epigenome-wide association study (EWAS) for HW, and identified 23 significantly associated CpGs, nine of which were also DM between ESTs. We found corresponding differences in transcript abundance between ESTs in three DMGs (ABLIM2, PITX2, and THRSP). Hypomethylation of an exonic CpG in PITX2 at 38.9°C was associated with increased expression, and suggests increased cell proliferation in broiler hatchlings incubated at higher temperatures. Overall, these results identified numerous epigenetic associations between chick incubation factors and heart development that may manifest in long-term differences in animal performance.

Leveraging locus-specific epigenetic heterogeneity to improve the performance of blood-based DNA methylation biomarkers

BACKGROUND

Variation in intercellular methylation patterns can complicate the use of methylation biomarkers for clinical diagnostic applications such as blood-based cancer testing. Here, we describe development and validation of a methylation density binary classification method called EpiClass (available for download at https://github.com/Elnitskilab/EpiClass ) that can be used to predict and optimize the performance of methylation biomarkers, particularly in challenging, heterogeneous samples such as liquid biopsies. This approach is based upon leveraging statistical differences in single-molecule sample methylation density distributions to identify ideal thresholds for sample classification.

RESULTS

We developed and tested the classifier using reduced representation bisulfite sequencing (RRBS) data derived from ovarian carcinoma tissue DNA and controls. We used these data to perform in silico simulations using methylation density profiles from individual epiallelic copies of ZNF154, a genomic locus known to be recurrently methylated in numerous cancer types. From these profiles, we predicted the performance of the classifier in liquid biopsies for the detection of epithelial ovarian carcinomas (EOC). In silico analysis indicated that EpiClass could be leveraged to better identify cancer-positive liquid biopsy samples by implementing precise thresholds with respect to methylation density profiles derived from circulating cell-free DNA (cfDNA) analysis. These predictions were confirmed experimentally using DREAMing to perform digital methylation density analysis on a cohort of low volume (1-ml) plasma samples obtained from 26 EOC-positive and 41 cancer-free women. EpiClass performance was then validated in an independent cohort of 24 plasma specimens, derived from a longitudinal study of 8 EOC-positive women, and 12 plasma specimens derived from 12 healthy women, respectively, attaining a sensitivity/specificity of 91.7%/100.0%. Direct comparison of CA-125 measurements with EpiClass demonstrated that EpiClass was able to better identify EOC-positive women than standard CA-125 assessment. Finally, we used independent whole genome bisulfite sequencing (WGBS) datasets to demonstrate that EpiClass can also identify other cancer types as well or better than alternative methylation-based classifiers.

CONCLUSIONS

Our results indicate that assessment of intramolecular methylation density distributions calculated from cfDNA facilitates the use of methylation biomarkers for diagnostic applications. Furthermore, we demonstrated that EpiClass analysis of ZNF154 methylation was able to outperform CA-125 in the detection of etiologically diverse ovarian carcinomas, indicating broad utility of ZNF154 for use as a biomarker of ovarian cancer.

Quantitative analysis questions the role of MeCP2 as a global regulator of alternative splicing

MeCP2 is an abundant protein in mature nerve cells, where it binds to DNA sequences containing methylated cytosine. Mutations in the MECP2 gene cause the severe neurological disorder Rett syndrome (RTT), provoking intensive study of the underlying molecular mechanisms. Multiple functions have been proposed, one of which involves a regulatory role in splicing. Here we leverage the recent availability of high-quality transcriptomic data sets to probe quantitatively the potential influence of MeCP2 on alternative splicing. Using a variety of machine learning approaches that can capture both linear and non-linear associations, we show that widely different levels of MeCP2 have a minimal effect on alternative splicing in three different systems. Alternative splicing was also apparently indifferent to developmental changes in DNA methylation levels. Our results suggest that regulation of splicing is not a major function of MeCP2. They also highlight the importance of multi-variate quantitative analyses in the formulation of biological hypotheses.

pWGBSSimla: a profile-based whole-genome bisulfite sequencing data simulator incorporating methylation QTLs, allele-specific methylations and differentially methylated regions

MOTIVATION

DNA methylation plays an important role in regulating gene expression. DNA methylation is commonly analyzed using bisulfite sequencing (BS-seq)-based designs, such as whole-genome bisulfite sequencing (WGBS), reduced representation bisulfite sequencing (RRBS) and oxidative bisulfite sequencing (oxBS-seq). Furthermore, there has been growing interest in investigating the roles that genetic variants play in changing the methylation levels (i.e. methylation quantitative trait loci or meQTLs), how methylation regulates the imprinting of gene expression (i.e. allele-specific methylation or ASM) and the differentially methylated regions (DMRs) among different cell types. However, none of the current simulation tools can generate different BS-seq data types (e.g. WGBS, RRBS and oxBS-seq) while modeling meQTLs, ASM and DMRs.

RESULTS

We developed profile-based whole-genome bisulfite sequencing data simulator (pWGBSSimla), a profile-based bisulfite sequencing data simulator, which simulates WGBS, RRBS and oxBS-seq data for different cell types based on real data. meQTLs and ASM are modeled based on the block structures of the methylation status at CpGs, whereas the simulation of DMRs is based on observations of methylation rates in real data. We demonstrated that pWGBSSimla adequately simulates data and allows performance comparisons among different methylation analysis methods.

AVAILABILITY AND IMPLEMENTATION

pWGBSSimla is available at https://omicssimla.sourceforge.io.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Unique Epigenomic Landscape Defines Human Erythropoiesis

Mammalian erythropoiesis yields a highly specialized cell type, the mature erythrocyte, evolved to meet the organismal needs of increased oxygen-carrying capacity. To better understand the regulation of erythropoiesis, we performed genome-wide studies of chromatin accessibility, DNA methylation, and transcriptomics using a recently developed strategy to obtain highly purified populations of primary human erythroid cells. The integration of gene expression, DNA methylation, and chromatin state dynamics reveals that stage-specific gene regulation during erythropoiesis is a stepwise and hierarchical process involving many cis-regulatory elements. Erythroid-specific, nonpromoter sites of chromatin accessibility are linked to erythroid cell phenotypic variation and inherited disease. Comparative analyses of stage-specific chromatin accessibility indicate that there is limited early chromatin priming of erythroid genes during hematopoiesis. The epigenome of terminally differentiating erythroid cells defines a distinct subset of highly specialized cells that are vastly dissimilar from other hematopoietic and nonhematopoietic cell types. These epigenomic and transcriptome data are powerful tools to study human erythropoiesis.

Schulz et al. use genome-wide studies of chromatin accessibility, DNA methylation, and transcriptomes in primary human erythroid cells to reveal important characteristics of erythropoiesis. Chromatin accessibility of terminal erythroid differentiation is markedly dissimilar from other hematopoietic cell types. Epigenomic changes are linked to erythroid cell traits and disease genes.

The iNs and Outs of Direct Reprogramming to Induced Neurons

Understanding of cell-type specific transcription factors has promoted progress in methods for cellular reprogramming, such as directly reprogramming somatic cells to induced neurons (iN). Methods for direct reprogramming require neuronal-fate determining gene activation via neuron-specific microRNAs, chemical modulation of key neuronal signaling pathways or overexpression via viral vectors, with some reprogramming strategies requiring a combination of these methods to induce the neuronal-cell fate. These methods have been employed in a multitude of cell types, including fibroblasts, hepatocytes, peripheral blood mononuclear, and T cells. The ability to create iN from skin biopsies and blood samples coupled with recent advancements in artificially inducing age- and disease-associated phenotypes are accelerating the development of disease models for late-onset neurodegenerative disorders. Here, we review how activation of the neuronal transcriptome alters the epigenetic landscape of the donor cell to facilitate reprogramming to neurons. We also discuss the advantages of using DNA binding domains such as CRISPR/dCas9 to overcome epigenetic barriers to induce neuronal-cell fate by activating endogenous neuronal cell-fate determining genes.

Tumor Necrosis Factor α Influences Phenotypic Plasticity and Promotes Epigenetic Changes in Human Basal Forebrain Cholinergic Neuroblasts

TNFα is the main proinflammatory cytokine implicated in the pathogenesis of neurodegenerative disorders, but it also modulates physiological functions in both the developing and adult brain. In this study, we investigated a potential direct role of TNFα in determining phenotypic changes of a recently established cellular model of human basal forebrain cholinergic neuroblasts isolated from the nucleus basalis of Meynert (hfNBMs). Exposing hfNBMs to TNFα reduced the expression of immature markers, such as nestin and β-tubulin III, and inhibited primary cilium formation. On the contrary, TNFα increased the expression of TNFα receptor TNFR2 and the mature neuron marker MAP2, also promoting neurite elongation. Moreover, TNFα affected nerve growth factor receptor expression. We also found that TNFα induced the expression of DNA-methylation enzymes and, accordingly, downregulated genes involved in neuronal development through epigenetic mechanisms, as demonstrated by methylome analysis. In summary, TNFα showed a dual role on hfNBMs phenotypic plasticity, exerting a negative influence on neurogenesis despite a positive effect on differentiation, through mechanisms that remain to be elucidated. Our results help to clarify the complexity of TNFα effects in human neurons and suggest that manipulation of TNFα signaling could provide a potential therapeutic approach against neurodegenerative disorders.

STAT3 and GR Cooperate to Drive Gene Expression and Growth of Basal-Like Triple-Negative Breast Cancer

Breast cancers are divided into subtypes with different prognoses and treatment responses based on global differences in gene expression. Luminal breast cancer gene expression and proliferation are driven by estrogen receptor alpha, and targeting this transcription factor is the most effective therapy for this subtype. By contrast, it remains unclear which transcription factors drive the gene expression signature that defines basal-like triple-negative breast cancer, and there are no targeted therapies approved to treat this aggressive subtype. In this study, we utilized integrated genomic analysis of DNA methylation, chromatin accessibility, transcription factor binding, and gene expression in large collections of breast cancer cell lines and patient tumors to identify transcription factors responsible for the basal-like gene expression program. Glucocorticoid receptor (GR) and STAT3 bind to the same genomic regulatory regions, which were specifically open and unmethylated in basal-like breast cancer. These transcription factors cooperated to regulate expression of hundreds of genes in the basal-like gene expression signature, which were associated with poor prognosis. Combination treatment with small-molecule inhibitors of both transcription factors resulted in synergistic decreases in cell growth in cell lines and patient-derived organoid models. This study demonstrates that GR and STAT3 cooperate to regulate the basal-like breast cancer gene expression program and provides the basis for improved therapy for basal-like triple-negative breast cancer through rational combination of STAT3 and GR inhibitors. SIGNIFICANCE: This study demonstrates that GR and STAT3 cooperate to activate the canonical gene expression signature of basal-like triple-negative breast cancer and that combination treatment with STAT3 and GR inhibitors could provide synergistic therapeutic efficacy.

Correlative study of epigenetic regulation of tumor microenvironment in spindle cell melanomas and cutaneous malignant peripheral nerve sheath tumors

The tumor microenvironment (TME) plays critical roles in tumor growth and progression, however key regulators of gene expression in the TME of cutaneous malignant peripheral nerve sheath tumor (C-MPNST) and spindle cell melanoma (SCM) have not been well elucidated. Herein, we investigate the epigenetic regulation of promoters and gene bodies and their effect on the TME composition of C-MPNSTs and SCMs. A cohort of 30 patients was analyzed using differential gene expression (DGE) and gene set enrichment analysis (GSEA) using the Nanostring platform. Methylation analysis was carried out utilizing an Infinium Methylation EPIC array targeting 866,562 methylation site (CpG) islands. DGE revealed overexpression of genes related to mast cells in the TME of SCMs, and a predominance of exhausted CD8⁺ T cells and macrophages in the TME of C-MPNSTs. Interestingly, we further observed promoter hypermethylation in key overexpressed genes and corresponding gene body hypomethylation. Analysis using ENCODE ChIP-sequencing data identified CTCF as the common transcription factor at the site of the hypomethylated probe. These findings support that the TME composition of C-MPNSTs and SCMs is at least partially independent on promoter methylation status, suggesting a possible relationship between gene body enhancers and expression of key TME genes in both entities.

Discovery of new epigenomics-based biomarkers and the early diagnosis of neurodegenerative diseases

Treatment options for many neurodegenerative diseases are limited due to the lack of early diagnostic procedures that allow timely delivery of therapeutic agents to affected neurons prior to cell death. While notable advances have been made in neurodegenerative disease biomarkers, whether or not the biomarkers discovered to date are useful for early diagnosis remains an open question. Additionally, the reliability of these biomarkers has been disappointing, due in part to the large dissimilarities between the tissues traditionally used to source biomarkers and primarily diseased neurons. In this article, we review the potential viability of atypical epigenetic and/or consequent transcriptional alterations (ETAs) as biomarkers of early-stage neurodegenerative disease, and present our perspectives on the discovery and practical use of such biomarkers in patient-derived neural samples using single-cell level analyses, thereby greatly enhancing the reliability of biomarker application.

Profiling of immune related genes silenced in EBV-positive gastric carcinoma identified novel restriction factors of human gammaherpesviruses

EBV-associated gastric cancer (EBVaGC) is characterized by high frequency of DNA methylation. In this study, we investigated how epigenetic alteration of host genome contributes to pathogenesis of EBVaGC through the analysis of transcriptomic and epigenomic datasets from NIH TCGA (The Cancer Genome Atlas) consortium. We identified that immune related genes (IRGs) is a group of host genes preferentially silenced in EBV-positive gastric cancers through DNA hypermethylation. Further functional characterizations of selected IRGs reveal their novel antiviral activity against not only EBV but also KSHV. In particular, we showed that metallothionein-1 (MT1) and homeobox A (HOXA) gene clusters are down-regulated via EBV-driven DNA hypermethylation. Several MT1 isoforms suppress EBV lytic replication and release of progeny virions as well as KSHV lytic reactivation, suggesting functional redundancy of these genes. In addition, single HOXA10 isoform exerts antiviral activity against both EBV and KSHV. We also confirmed the antiviral effect of other dysregulated IRGs, such as IRAK2 and MAL, in scenario of EBV and KSHV lytic reactivation. Collectively, our results demonstrated that epigenetic silencing of IRGs is a viral strategy to escape immune surveillance and promote viral propagation, which is overall beneficial to viral oncogenesis of human gamma-herpesviruses (EBV and KSHV), considering that these IRGs possess antiviral activities against these oncoviruses.

Pilot study of combined aerobic and resistance exercise on fatigue for patients with head and neck cancer: Inflammatory and epigenetic changes

This pilot study examined whether a combined aerobic resistance exercise program reduced fatigue and the potential inflammatory and epigenetic mechanisms in patients with head and neck cancer (HNC) receiving intensity-modulated radiotherapy. The exercise group (N = 12) received a 3-month supervised aerobic resistance exercise intervention that was initiated before a 6-week radiotherapy regimen; the control group (N = 14) received standard care. Fatigue was measured using Multidimensional Fatigue Inventory-20; physical function measures included a 6-minute walk distance (6MWD), chair stands, bicep curls, and hand grip strength. Inflammatory markers and DNA methylation data were acquired using standardized protocol. Patients were mostly white (93%) and male (81%) with a mean age of 57 years. At the end of the intervention, the exercise group had a marginal decrease in fatigue compared with the control (-5.0 vs. 4.9; P = 0.10). The exercise group had a significantly greater improvement in 6MWD (29.8 vs. -55.5 m; P = 0.04), and a marginally smaller decline in hand grip (-0.3 vs. -5.8 lbs; P = 0.05) at the end of the intervention than the control. No significant difference in inflammatory markers was observed between groups. Lower plasma interleukin (IL) 6, IL1 receptor antagonist, tumor necrosis factor α (TNFα), soluble TNF receptor II and C-reactive protein were significantly associated with increased 6MWD, chair stand, and bicep curl at the end of the intervention (p < 0.05). Among the 1152 differentially methylated sites (DMS) after intervention (p < 0.001), 163 DMS were located in gene promoter regions. Enrichment analysis suggested that the top 10 upstream regulators were associated with tumor (HNF4A, RPP38, HOXA9, SAHM1, CDK7, NDN, RPS15) and inflammation (IRF7, CRKL, ONECUT1). The top 5 diseases or functions annotations of the 62 hypermethylated DMS indicated anti-tumor and anti-inflammatory effects that might be linked to exercise. These findings suggest that exercise may improve physical performance and reduce fatigue, which could be further linked to decreased inflammation, during active radiotherapy for HNC patients. Larger studies are warranted.

Comparative whole genome DNA methylation profiling across cattle tissues reveals global and tissue-specific methylation patterns

BACKGROUND

Efforts to improve animal health, and understand genetic bases for production, may benefit from a comprehensive analysis of animal genomes and epigenomes. Although DNA methylation has been well studied in humans and other model species, its distribution patterns and regulatory impacts in cattle are still largely unknown. Here, we present the largest collection of cattle DNA methylation epigenomic data to date.

RESULTS

Using Holstein cattle, we generated 29 whole genome bisulfite sequencing (WGBS) datasets for 16 tissues, 47 corresponding RNA-seq datasets, and 2 whole genome sequencing datasets. We did read mapping and DNA methylation calling based on two different cattle assemblies, demonstrating the high quality of the long-read-based assembly markedly improved DNA methylation results. We observed large differences across cattle tissues in the methylation patterns of global CpG sites, partially methylated domains (PMDs), hypomethylated regions (HMRs), CG islands (CGIs), and common repeats. We detected that each tissue had a distinct set of PMDs, which showed tissue-specific patterns. Similar to human PMD, cattle PMDs were often linked to a general decrease of gene expression and a decrease in active histone marks and related to long-range chromatin organizations, like topologically associated domains (TADs). We tested a classification of the HMRs based on their distributions relative to transcription start sites (TSSs) and detected tissue-specific TSS-HMRs and genes that showed strong tissue effects. When performing cross-species comparisons of paired genes (two opposite strand genes with their TSS located in the same HMR), we found out they were more consistently co-expressed among human, mouse, sheep, goat, yak, pig, and chicken, but showed lower consistent ratios in more divergent species. We further used these WGBS data to detect 50,023 experimentally supported CGIs across bovine tissues and found that they might function as a guard against C-to-T mutations for TSS-HMRs. Although common repeats were often heavily methylated, some young Bov-A2 repeats were hypomethylated in sperm and could affect the promoter structures by exposing potential transcription factor binding sites.

CONCLUSIONS

This study provides a comprehensive resource for bovine epigenomic research and enables new discoveries about DNA methylation and its role in complex traits.

Absence of MeCP2 binding to non-methylated GT-rich sequences in vivo

MeCP2 is a nuclear protein that binds to sites of cytosine methylation in the genome. While most evidence confirms this epigenetic mark as the primary determinant of DNA binding, MeCP2 is also reported to have an affinity for non-methylated DNA sequences. Here we investigated the molecular basis and in vivo significance of its reported affinity for non-methylated GT-rich sequences. We confirmed this interaction with isolated domains of MeCP2 in vitro and defined a minimal target DNA sequence. Binding depends on pyrimidine 5′ methyl groups provided by thymine and requires adjacent guanines and a correctly orientated A/T-rich flanking sequence. Unexpectedly, full-length MeCP2 protein failed to bind GT-rich sequences in vitro. To test for MeCP2 binding to these motifs in vivo, we analysed human neuronal cells using ChIP-seq and ATAC-seq technologies. While both methods robustly detected DNA methylation-dependent binding of MeCP2 to mCG and mCAC, neither showed evidence of MeCP2 binding to GT-rich motifs. The data suggest that GT binding is an in vitro phenomenon without in vivo relevance. Our findings argue that MeCP2 does not read unadorned DNA sequence and therefore support the notion that its primary role is to interpret epigenetic modifications of DNA.

Interplay between Metabolism, Nutrition and Epigenetics in Shaping Brain DNA Methylation, Neural Function and Behavior

Gene expression in the brain is dramatically regulated by a variety of stimuli. While the role of neural activity has been extensively studied, less is known about the effects of metabolism and nutrition on transcriptional control mechanisms in the brain. Extracellular signals are integrated at the chromatin level through dynamic modifications of epigenetic marks, which in turn fine-tune gene transcription. In the last twenty years, it has become clear that epigenetics plays a crucial role in modulating central nervous system functions and finally behavior. Here, we will focus on the effect of metabolic signals in shaping brain DNA methylation, both during development and adulthood. We will provide an overview of maternal nutrition effects on brain methylation and behavior in offspring. In addition, the impact of different diet challenges on cytosine methylation dynamics in the adult brain will be discussed. Finally, the possible role played by the metabolic status in modulating DNA hydroxymethylation, which is particularly abundant in neural tissue, will be considered.

Matching cell lines with cancer type and subtype of origin via mutational, epigenomic, and transcriptomic patterns

Cell lines are commonly used as cancer models. The tissue of origin provides context for understanding biological mechanisms and predicting therapy response. We therefore systematically examined whether cancer cell lines exhibit features matching the presumed cancer type of origin. Gene expression and DNA methylation classifiers trained on ~9000 tumors identified 35 (of 614 examined) cell lines that better matched a different tissue or cell type than the one originally assigned. Mutational patterns further supported most reassignments. For instance, cell lines identified as originating from the skin often exhibited a UV mutational signature. We cataloged 366 "golden set" cell lines in which transcriptomic and epigenomic profiles strongly resemble the cancer type of origin, further proposing their assignments to subtypes. Accounting for the uncertain tissue of origin in cell line panels can change the interpretation of drug screening and genetic screening data, revealing previously unknown genomic determinants of sensitivity or resistance.

An integrative framework of heterogeneous genomic data for cancer dynamic modules based on matrix decomposition

Cancer progression is dynamic, and tracking dynamic modules is promising for cancer diagnosis and therapy. Accumulated genomic data provide us an opportunity to investigate the underlying mechanisms of cancers. However, as far as we know, no algorithm has been designed for dynamic modules by integrating heterogeneous omics data. To address this issue, we propose an integrative framework for dynamic module detection based on regularized nonnegative matrix factorization method (DrNMF) by integrating the gene expression and protein interaction network. To remove the heterogeneity of genomic data, we divide the samples of expression profiles into groups to construct gene co-expression networks. To characterize the dynamics of modules, the temporal smoothness framework is adopted, in which the gene co-expression network at the previous stage and protein interaction network are incorporated into the objective function of DrNMF via regularization. The experimental results demonstrate that DrNMF is superior to state-of-the-art methods in terms of accuracy. For breast cancer data, the obtained dynamic modules are more enriched by the known pathways, and can be used to predict the stages of cancers and survival time of patients. The proposed model and algorithm provide an effective integrative analysis of heterogeneous genomic data for cancer progression.

Klotho and the Treatment of Human Malignancies

Klotho was first discovered as an anti-ageing protein linked to a number of age-related disease processes, including cardiovascular, renal, musculoskeletal, and neurodegenerative conditions. Emerging research has also demonstrated a potential therapeutic role for Klotho in cancer biology, which is perhaps unsurprising given that cancer and ageing share similar molecular hallmarks. In addition to functioning as a tumour suppressor in numerous solid tumours and haematological malignancies, Klotho represents a candidate therapeutic target for patients with these diseases, the majority of whom have limited treatment options. Here, we examine contemporary evidence evaluating the anti-neoplastic effects of Klotho and describe the modulation of downstream oncogenic signalling pathways, including Wnt/β-catenin, FGF, IGF1, PIK3K/AKT, TGFβ, and the Unfolded Protein Response. We also discuss possible approaches to developing therapeutic Klotho and consider technological advances that may facilitate the delivery of Klotho through gene therapy.

Differential DNA methylation analysis across the promoter regions using methylated DNA immunoprecipitation sequencing profiling of porcine loin muscle

Background and Aim:

Pork leanness and marbling are among the essential traits of consumer preference. To acquire knowledge about universal epigenetic regulations for improving breed selection, a meta-analysis of methylated DNA immunoprecipitation sequencing (MeDIP-seq) profiling data of mixed loin muscle types was performed in this study.

Materials and Methods:

MeDIP-seq profiling datasets of longissimus dorsi muscle and psoas major muscles from male and female pigs of Landrace and Tibetan breeds were preprocessed and aligned to the porcine genome. Analysis of differential methylated DNA regions (DMRs) between the breeds was performed by focusing on transcription start sites (TSSs) of known genes (−20,000-3000 bases from TSS). All associated genes were further reviewed for their functions and predicted for transcription factors (TF) possibly associated with their TSSs.

Results:

When the methylation levels of DMRs in TSS regions of Landrace breed were compared to those of Tibetan breed, 10 DMRs were hypomethylated (Landrace < Tibetan), and 19 DMRs were hypermethylated (Landrace > Tibetan), accordingly (p≤0.001). According to the reviews about gene functions, all associated genes were pieces of evidence for their roles in a variety of muscle and lipid metabolisms. Prediction of the binding TFs revealed the six most abundant binding TFs to such DMRs-associated TSS (p≤0.0001) as follows: ZNF384, Foxd3, IRF1, KLF9, EWSR1-FLI1, HES5, and TFAP2A.

Conclusion:

Common DMRs-associated TSS between the lean-type and the marbled-type loin muscles were identified in this study. Interestingly, the genes associated with such regions were strongly evidenced for their possible roles on the muscle trait characteristics by which further novel research topics could be focused on them in the future.

Weighted elastic net for unsupervised domain adaptation with application to age prediction from DNA methylation data

MOTIVATION

Predictive models are a powerful tool for solving complex problems in computational biology. They are typically designed to predict or classify data coming from the same unknown distribution as the training data. In many real-world settings, however, uncontrolled biological or technical factors can lead to a distribution mismatch between datasets acquired at different times, causing model performance to deteriorate on new data. A common additional obstacle in computational biology is scarce data with many more features than samples. To address these problems, we propose a method for unsupervised domain adaptation that is based on a weighted elastic net. The key idea of our approach is to compare dependencies between inputs in training and test data and to increase the cost of differently behaving features in the elastic net regularization term. In doing so, we encourage the model to assign a higher importance to features that are robust and behave similarly across domains.

RESULTS

We evaluate our method both on simulated data with varying degrees of distribution mismatch and on real data, considering the problem of age prediction based on DNA methylation data across multiple tissues. Compared with a non-adaptive standard model, our approach substantially reduces errors on samples with a mismatched distribution. On real data, we achieve far lower errors on cerebellum samples, a tissue which is not part of the training data and poorly predicted by standard models. Our results demonstrate that unsupervised domain adaptation is possible for applications in computational biology, even with many more features than samples.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/PfeiferLabTue/wenda.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

PRISM: methylation pattern-based, reference-free inference of subclonal makeup

Motivation

Characterizing cancer subclones is crucial for the ultimate conquest of cancer. Thus, a number of bioinformatic tools have been developed to infer heterogeneous tumor populations based on genomic signatures such as mutations and copy number variations. Despite accumulating evidence for the significance of global DNA methylation reprogramming in certain cancer types including myeloid malignancies, none of the bioinformatic tools are designed to exploit subclonally reprogrammed methylation patterns to reveal constituent populations of a tumor. In accordance with the notion of global methylation reprogramming, our preliminary observations on acute myeloid leukemia (AML) samples implied the existence of subclonally occurring focal methylation aberrance throughout the genome.

Results

We present PRISM, a tool for inferring the composition of epigenetically distinct subclones of a tumor solely from methylation patterns obtained by reduced representation bisulfite sequencing. PRISM adopts DNA methyltransferase 1-like hidden Markov model-based in silico proofreading for the correction of erroneous methylation patterns. With error-corrected methylation patterns, PRISM focuses on a short individual genomic region harboring dichotomous patterns that can be split into fully methylated and unmethylated patterns. Frequencies of such two patterns form a sufficient statistic for subclonal abundance. A set of statistics collected from each genomic region is modeled with a beta-binomial mixture. Fitting the mixture with expectation-maximization algorithm finally provides inferred composition of subclones. Applying PRISM for two AML samples, we demonstrate that PRISM could infer the evolutionary history of malignant samples from an epigenetic point of view.

Availability and implementation

PRISM is freely available on GitHub (https://github.com/dohlee/prism).

Supplementary information

Supplementary data are available at Bioinformatics online.

Massively parallel reporter assays of melanoma risk variants identify MX2 as a gene promoting melanoma

Genome-wide association studies (GWAS) have identified ~20 melanoma susceptibility loci, most of which are not functionally characterized. Here we report an approach integrating massively-parallel reporter assays (MPRA) with cell-type-specific epigenome and expression quantitative trait loci (eQTL) to identify susceptibility genes/variants from multiple GWAS loci. From 832 high-LD variants, we identify 39 candidate functional variants from 14 loci displaying allelic transcriptional activity, a subset of which corroborates four colocalizing melanocyte cis-eQTL genes. Among these, we further characterize the locus encompassing the HIV-1 restriction gene, MX2 (Chr21q22.3), and validate a functional intronic variant, rs398206. rs398206 mediates the binding of the transcription factor, YY1, to increase MX2 levels, consistent with the cis-eQTL of MX2 in primary human melanocytes. Melanocyte-specific expression of human MX2 in a zebrafish model demonstrates accelerated melanoma formation in a BRAF^V600E background. Our integrative approach streamlines GWAS follow-up studies and highlights a pleiotropic function of MX2 in melanoma susceptibility.

There are more than 20 known melanoma susceptibility genes. Here, using a massively parallel reporter assay, the authors identify risk-associated variants that alter gene transcription, and demonstrate that expression of one such gene, MX2, leads to the promotion of melanoma in a zebrafish model.

Insight into norcantharidin, a small-molecule synthetic compound with potential multi-target anticancer activities

Norcantharidin (NCTD) is a demethylated derivative of cantharidin, which is an anticancer active ingredient of traditional Chinese medicine, and is currently used clinically as a routine anti-cancer drug in China. Clarifying the anticancer effect and molecular mechanism of NCTD is critical for its clinical application. Here, we summarized the physiological, chemical, pharmacokinetic characteristics and clinical applications of NCTD. Besides, we mainly focus on its potential multi-target anticancer activities and underlying mechanisms, and discuss the problems existing in clinical application and scientific research of NCTD, so as to provide a potential anticancer therapeutic agent for human malignant tumors.

Epigenetic inheritance of circadian period in clonal cells

Circadian oscillations are generated via transcriptional-translational negative feedback loops. However, individual cells from fibroblast cell lines have heterogeneous rhythms, oscillating independently and with different period lengths. Here we showed that heterogeneity in circadian period is heritable and used a multi-omics approach to investigate underlying mechanisms. By examining large-scale phenotype-associated gene expression profiles in hundreds of mouse clonal cell lines, we identified and validated multiple novel candidate genes involved in circadian period determination in the absence of significant genomic variants. We also discovered differentially co-expressed gene networks that were functionally associated with period length. We further demonstrated that global differential DNA methylation bidirectionally regulated these same gene networks. Interestingly, we found that depletion of DNMT1 and DNMT3A had opposite effects on circadian period, suggesting non-redundant roles in circadian gene regulation. Together, our findings identify novel gene candidates involved in periodicity, and reveal DNA methylation as an important regulator of circadian periodicity.

TET1 and 5-Hydroxymethylation Preserve the Stem Cell State of Mouse Trophoblast

The ten-eleven translocation factor TET1 and its conferred epigenetic modification 5-hydroxymethylcytosine (5hmC) have important roles in maintaining the pluripotent state of embryonic stem cells (ESCs). We previously showed that TET1 is also essential to maintain the stem cell state of trophoblast stem cells (TSCs). Here, we establish an integrated panel of absolute 5hmC levels, genome-wide DNA methylation and hydroxymethylation patterns, transcriptomes, and TET1 chromatin occupancy in TSCs and differentiated trophoblast cells. We show that the combined presence of 5-methylcytosine (5mC) and 5hmC correlates with transcriptional activity of associated genes. Hypoxia can slow down the global loss of 5hmC that occurs upon differentiation of TSCs. Notably, unlike in ESCs and epiblast cells, most TET1-bound regions overlap with active chromatin marks and TFAP2C binding sites and demarcate putative trophoblast enhancer regions. These chromatin modification and occupancy patterns are highly informative to identify novel candidate regulators of the TSC state.

•

5hmC to 5mC ratios correlate with gene activity in TS cells

•

TS cell differentiation-associated loss of 5hmC is slowed down in hypoxia

•

TET1 binding in TS cells forms long-range interactions with key trophoblast genes

•

Intergenic TET1 binding sites in TS cells demarcate putative trophoblast enhancers

Integrated data analysis reveals significant associations of KEAP1 mutations with DNA methylation alterations in lung adenocarcinomas

KEAP1 regulates the cytoprotection induced by NRF2 and has been reported to be a candidate tumor suppressor. Recent evidence has shown that mutations in several driver genes cause aberrant DNA methylation patterns, a hallmark of cancer. However, the correlation between KEAP1 mutations and DNA methylation in lung cancer has still not been investigated. In this study, we systematically carried out an integrated multi-omics analysis to explore the correlation between KEAP1 mutations and DNA methylation and its effect on gene expression in lung adenocarcinoma (LUAD). We found that most of the DNA aberrations associated with KEAP1 mutations in LAUD were hypomethylation. Surprisingly, we found several NRF2-regulated genes among the genes that showed differential DNA methylation. Moreover, we identified an 8-gene signature with altered DNA methylation pattern and elevated gene expression levels in LUAD patients with mutated KEAP1, and evaluated the prognostic value of this signature in various clinical datasets. These results establish that KEAP1 mutations are associated with DNA methylation changes capable of shaping regulatory network functions. Combining both epigenomic and transcriptomic changes along with KEAP1 mutations may provide a better understanding of the molecular mechanisms associated with the progression of lung cancer and may help to provide better therapeutic approaches.

Tissue-specific epigenetic inheritance after paternal heat exposure in male wild guinea pigs

External temperature change has been shown to modify epigenetic patterns, such as DNA methylation, which regulates gene expression. DNA methylation is heritable, and as such provides a mechanism to convey environmental information to subsequent generations. Studies on epigenetic response to temperature increase are still scarce in wild mammals, even more so studies that compare tissue-specific epigenetic responses. Here, we aim to address differential epigenetic responses on a gene and gene pathway level in two organs, liver and testis. We chose these organs, because the liver is the main metabolic and thermoregulation organ, and epigenetic modifications in testis are potentially transmitted to the F2 generation. We focused on the transmission of DNA methylation changes to naive male offspring after paternal exposure to an ambient temperature increase of 10 °C, and investigated differential methylated regions of sons sired before and after the paternal exposure using Reduced Representation Bisulfite Sequencing. We detected both a highly tissue-specific epigenetic response, reflected in genes involved in organ-specific metabolic pathways, and a more general regulation of single genes epigenetically modified in both organs. We conclude that genomes are context-specifically differentially epigenetically regulated in response to temperature increase. These findings emphasize the epigenetic relevance in cell differentiation, which is essential for the specific function(s) of complex organs, and is represented in a diverse molecular regulation of genes and gene pathways. The results also emphasize the paternal contribution to adaptive processes.

Genome-Wide DNA Methylation Enhances Stemness in the Mechanical Selection of Tumor-Repopulating Cells

Background

DNA methylation plays essential roles in tumor occurrence and stemness maintenance. Tumor-repopulating cells (TRCs) are cancer stem cell (CSC)-like cells with highly tumorigenic and self-renewing abilities, which were selected from tumor cells in soft three-dimensional (3D) fibrin gels.

Methods

Here, we presented a genome-wide map of methylated cytosines for time-series samples in TRC selection, in a 3D culture using whole-genome bisulfite sequencing (WGBS).

Results

A comparative analysis revealed that the methylation degrees of many differentially methylated genes (DMGs) were increased by the mechanical environment and changed from 2D rigid to 3D soft. DMGs were significantly enriched in stemness-related terms. In 1-day, TRCs had the highest non-CG methylation rate indicating its strong stemness. We found that genes with continuously increasing or decreasing methylation like CREB5/ADAMTS6/LMX1A may also affect the TRC screening process. Furthermore, results showed that stage-specific/common CSCs markers were biased toward changing their methylation in non-CG (CHG and CHH, where H corresponds to A, T, or C) methylation and enriched in gene body region.

Conclusions

WGBS provides DNA methylome in TRC screening. It was confirmed that non-CG DNA methylation plays an important role in TRC selection, which indicates that it is more sensitive to mechanical microenvironments and affects TRCs by regulating the expression of stemness genes in tumor cells.

Methylation in HT22 cells and primary hippocampal neurons with and without isoflurane exposurewhether isoflurane causes

Background

Epigenetic modulation may play a role in anesthesia related phenotypes, such as cognitive impairment or memory loss, especially with exposure to anesthetics in the vulnerable phase of brain development. While isoflurane anesthesia can evoke neuroinflammation and neuroapoptosis in young animals, we investigated in a permanent hippocampal cell line (HT22) and in primary hippocampal neurons in an a priori in vitro analysis, whether isoflurane exposure 1) evokes DNA methylation changes in genes involved in apoptosis and inflammation, and 2) results observed in a permanent hippocampal cell line are comparable to primary hippocampal neurons. In case of methylation changes in specific genes, (3) mRNA analysis was performed to assess possible effects on gene expression.

Methods

HT22 cells and primary mouse hippocampal neurons were exposed to 3% isoflurane for 4 h and DNA (each 6 single experiments) and RNA (3 single independent experiments) were extracted. Methylation analysis (EpiTect Methyl II PCR Array Systems, Qiagen) included the methylation status of 66 genes involved in apoptosis, cytokine production, inflammatory response, and autoimmunity. Quantitative Real-Time PCR was performed using the Quantitect SYBR Green Kit on a Step One Plus.

Results

Methylation status was markedly different between immortalized HT22 cells and cultured primary hippocampal neurons without isoflurane exposure. Of 66 genes investigated, 29 were methylated to a significantly greater degree in HT22 cells compared to primary hippocampal neurons. In cultured primary hippocampal neurons, in contrast, there was a greater methylation in several genes involved in inflammation, accompanied with significant downregulation of C-X-C motif chemokine 12 with isoflurane exposure (p = 0.023).

Conclusions

We demonstrate marked differences in gene methylation between HT22 cells and cultured primary hippocampal neurons without isoflurane exposure, with a greater methylation of several genes involved in inflammation upon isoflurane exposure and significant downregulation of Cxcl12 mRNA expression in primary hippocampal neurons. Accordingly, further investigations of anesthesia related DNA methylation should be performed with special consideration being given to the choice of cells targeted for such investigations.

The Molecular Basis of MeCP2 Function in the Brain

MeCP2 is a reader of the DNA methylome that occupies a large proportion of the genome due to its high abundance and the frequency of its target sites. It has been the subject of extensive study because of its link with 'MECP2-related disorders', of which Rett syndrome is the most prevalent. This review integrates evidence from patient mutation data with results of experimental studies using mouse models, cell lines and in vitro systems to critically evaluate our understanding of MeCP2 protein function. Recent evidence challenges the idea that MeCP2 is a multifunctional hub that integrates diverse processes to underpin neuronal function, suggesting instead that its primary role is to recruit the NCoR1/2 co-repressor complex to methylated sites in the genome, leading to dampening of gene expression.

Identification of methylation states of DNA regions for Illumina methylation BeadChip

Background

Methylation of cytosine bases in DNA is a critical epigenetic mark in many eukaryotes and has also been implicated in the development and progression of normal and diseased cells. Therefore, profiling DNA methylation across the genome is vital to understanding the effects of epigenetic. In recent years the Illumina HumanMethylation450 (HM450K) and MethylationEPIC (EPIC) BeadChip have been widely used to profile DNA methylation in human samples. The methods to predict the methylation states of DNA regions based on microarray methylation datasets are critical to enable genome-wide analyses.

Result

We report a computational approach based on the two layers two-state hidden Markov model (HMM) to identify methylation states of single CpG site and DNA regions in HM450K and EPIC BeadChip. Using this mothed, all CpGs detected by HM450K and EPIC in H1-hESC and GM12878 cell lines are identified as un-methylated, middle-methylated and full-methylated states. A large number of DNA regions are segmented into three methylation states as well. Comparing the identified regions with the result from the whole genome bisulfite sequencing (WGBS) datasets segmented by MethySeekR, our method is verified. Genome-wide maps of chromatin states show that methylation state is inversely correlated with active histone marks. Genes regulated by un-methylated regions are expressed and regulated by full-methylated regions are repressed. Our method is illustrated to be useful and robust.

Conclusion

Our method is valuable for DNA methylation genome-wide analyses. It is focusing on identification of DNA methylation states on microarray methylation datasets. For the features of array datasets, using two layers two-state HMM to identify to methylation states on CpG sites and regions creatively, our method which takes into account the distribution of genome-wide methylation levels is more reasonable than segmentation with a fixed threshold.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6019-0) contains supplementary material, which is available to authorized users.

Methylation Status of the Nanog Promoter Determines the Switch between Cancer Cells and Cancer Stem Cells

Cancer stem cells (CSCs) are the main cause of tumor development, metastasis, and relapse. CSCs are thus considered promising targets for cancer therapy. However, it is hard to eradicate CSCs due to their inherent plasticity and heterogeneity, and the underlying mechanism of the switch between non-CSCs and CSCs remains unclear. Here, it is shown that miR-135a combined with SMYD4 activates Nanog expression and induces the switch of non-CSCs into CSCs. The miR-135a level, once elevated, lowers the methylation level of the CG5 site in the Nanog promoter by directly targeting DNMT1. SMYD4 binds to the unmethylated Nanog promoter to activate Nanog expression in Nanog-negative tumor cells. The in vivo regulation of miR-135a levels could significantly affect both the CSCs proportion and tumor progression. These findings indicate that DNA methylation of the Nanog promoter modulates the switch of non-CSCs into CSCs under the control of the miRNA-135 level. In addition, the related pathways, miR-135a/DNMT1 and SMYD4, involved in these processes are potential targets for CSC-targeted therapy.

Epigenetics, Development, and Psychopathology

Epigenetic mechanisms govern the transcription of the genome. Research with model systems reveals that environmental conditions can directly influence epigenetic mechanisms that are associated with interindividual differences in gene expression in brain and neural function. In this review, we provide a brief overview of epigenetic mechanisms and research with relevant rodent models. We emphasize more recent translational research programs in epigenetics as well as the challenges inherent in the integration of epigenetics into developmental and clinical psychology. Our objectives are to present an update with respect to the translational relevance of epigenetics for the study of psychopathology and to consider the state of current research with respect to its potential importance for clinical research and practice in mental health.