Dynamic changes in the human methylome during differentiation

Analysis of mRNA Expression and DNA Methylation Datasets According to the Genomic Distribution of CpG Sites in Osteoarthritis

Objectives Transcriptional changes in cartilage can impact function by causing degradation such as that which occurs during the development of osteoarthritis (OA). Epigenetic regulation may be a key factor leading to transcriptional changes in OA. In this study, we performed a combined analysis of DNA methylation and gene expression microarray datasets and identified key transcription factors (TFs) central to the regulation of gene expression in OA. Methods A DNA methylation profile dataset (GSE63106) and a gene expression profiling dataset (GSE114007) were extracted from the Gene Expression Omnibus (GEO). We used ChAMP methylation analysis and the Limma package to identify differentially methylation genes (DMGs) and differentially expressed genes (DEGs) from normal and human knee cartilage samples in OA. Function enrichment analysis of DMGs was conducted using the DAVID database. A combined analysis of DEGs and DMGs was conducted to identify key TFs in OA. We then validated the mRNA expression of selected TFs in normal and OA cartilage by RT-qPCR. Primary chondrocytes were cultured and treated with the DNA methylation inhibitor 5-Aza-2-deoxycytidine (5-Aza) for functional validation. Results We identified 2,170 differential methylation sites (DMS) containing 1005 genes and 1986 DEGs between normal human and OA cartilage. Functional analysis of DMGs revealed that focal adhesion, extracellular matrix (ECM)-receptor interactions, the PI3K-Akt signaling pathway, and the FoxO signaling pathway were involved in OA. Integrated analysis showed a subset of 17 TFs. Four TFs (ELF3, SOX11, RARA, and FOXD2) were validated. RT-qPCR results showed the mRNA expression of SOX11, RARA, and FOXD2 were consistent with the results from the mRNA expression data. However, the expression of ELF3 could not be validated. Upon 5-Aza-2'-deoxycytidine (5-Aza) treatment, the mRNA levels of ELF3 and SOX11 were down-regulated, whilst RARA was up-regulated, and FOXD2 showed no significant change in expression level. Conclusions the effect of DNA methylation on the transcriptional regulation is related to the distribution of methylated sites across the genome. Epigenetic studies on the positions of DMS in transcriptional units can inform a better understanding of the function of DNA methylation and its transcription regulation.

Low biological fluctuation of mitochondrial CpG and non-CpG methylation at the single-molecule level

Mammalian cytosine DNA methylation (5mC) is associated with the integrity of the genome and the transcriptional status of nuclear DNA. Due to technical limitations, it has been less clear if mitochondrial DNA (mtDNA) is methylated and whether 5mC has a regulatory role in this context. Here, we used bisulfite-independent single-molecule sequencing of native human and mouse DNA to study mitochondrial 5mC across different biological conditions. We first validated the ability of long-read nanopore sequencing to detect 5mC in CpG (5mCpG) and non-CpG (5mCpH) context in nuclear DNA at expected genomic locations (i.e. promoters, gene bodies, enhancers, and cell type-specific transcription factor binding sites). Next, using high coverage nanopore sequencing we found low levels of mtDNA CpG and CpH methylation (with several exceptions) and little variation across biological processes: differentiation, oxidative stress, and cancer. 5mCpG and 5mCpH were overall higher in tissues compared to cell lines, with small additional variation between cell lines of different origin. Despite general low levels, global and single-base differences were found in cancer tissues compared to their adjacent counterparts, in particular for 5mCpG. In conclusion, nanopore sequencing is a useful tool for the detection of modified DNA bases on mitochondria that avoid the biases introduced by bisulfite and PCR amplification. Enhanced nanopore basecalling models will provide further resolution on the small size effects detected here, as well as rule out the presence of other DNA modifications such as oxidized forms of 5mC.

Single-molecule micromanipulation studies of methylated DNA

Cytosine methylated at the five-carbon position is the most widely studied reversible DNA modification. Prior findings indicate that methylation can alter mechanical properties. However, those findings were qualitative and sometimes contradictory, leaving many aspects unclear. By applying single-molecule magnetic force spectroscopy techniques allowing for direct manipulation and dynamic observation of DNA mechanics and mechanically driven strand separation, we investigated how CpG and non-CpG cytosine methylation affects DNA micromechanical properties. We quantitatively characterized DNA stiffness using persistence length measurements from force-extension curves in the nanoscale length regime and demonstrated that cytosine methylation results in longer contour length and increased DNA flexibility (i.e., decreased persistence length). In addition, we observed the preferential formation of plectonemes over unwound single-stranded "bubbles" of DNA under physiologically relevant stretching forces and supercoiling densities. The flexibility and high structural stability of methylated DNA is likely to have significant consequences on the recruitment of proteins recognizing cytosine methylation and DNA packaging.

Computing the Role of Alternative Splicing in Cancer

Most human genes undergo alternative splicing (AS), and dysregulation of alternative splicing contributes to tumor initiation and progression. Computational analysis of genomic and transcriptomic data enables the systematic characterization of alternative splicing and its functional role in cancer. In this review, we summarize the latest computational approaches to studying alternative splicing in cancer and the current limitations of the most popular tools in this field. Finally, we describe some of the current computational challenges in the characterization of the role of alternative splicing in cancer.

A region-based method for causal mediation analysis of DNA methylation data

Exposure to environmental factors can affect DNA methylation at a 5'-cytosine-phosphate-guanine-3' (CpG) site or a genomic region, which can then affect an outcome. In other words, environmental effects on an outcome could be mediated by DNA methylation. To date, single CpG-site-based mediation analysis has been employed extensively. More recently, however, there has been considerable interest in studying differentially methylated regions (DMRs), both because DMRs are more likely to have functional effects than single CpG sites and because testing DMRs reduces multiple testing. In this report, we propose a novel causal mediation approach under the counterfactual framework to test the significance of total (TE), direct (DE), and indirect effects (IE) of predictors on response variable with a methylated region (MR) as the mediator (denoted as MR-Mediation). Functional linear transformation is used to reduce the possible high dimension of the CpG sites in a predefined MR and to account for their location information. In our simulation studies, MR-Mediation retained the desired Type I error rates for TE, DE, and IE tests. Furthermore, MR-Mediation had better power performance than testing mean methylation level as the mediator in most considered scenarios, especially for IE (i.e., mediated effect) test, which could be more interesting than the other two effect tests. We further illustrate our proposed method by analysing the methylation mediated effect of exposure to gun violence on total immunoglobulin E or atopic asthma among participants in the Epigenetic Variation and Childhood Asthma in Puerto Ricans study.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

DNA Methylation, Mechanisms of FMR1 Inactivation and Therapeutic Perspectives for Fragile X Syndrome

Among the inherited causes of intellectual disability and autism, Fragile X syndrome (FXS) is the most frequent form, for which there is currently no cure. In most FXS patients, the FMR1 gene is epigenetically inactivated following the expansion over 200 triplets of a CGG repeat (FM: full mutation). FMR1 encodes the Fragile X Mental Retardation Protein (FMRP), which binds several mRNAs, mainly in the brain. When the FM becomes methylated at 10-12 weeks of gestation, the FMR1 gene is transcriptionally silent. The molecular mechanisms involved in the epigenetic silencing are not fully elucidated. Among FXS families, there is a rare occurrence of males carrying a FM, which remains active because it is not methylated, thus ensuring enough FMRPs to allow for an intellectual development within normal range. Which mechanisms are responsible for sparing these individuals from being affected by FXS? In order to answer this critical question, which may have possible implications for FXS therapy, several potential epigenetic mechanisms have been described. Here, we focus on current knowledge about the role of DNA methylation and other epigenetic modifications in FMR1 gene silencing.

Epigenetics: New Insights into Mammary Gland Biology

The mammary gland undergoes important anatomical and physiological changes from embryogenesis through puberty, pregnancy, lactation and involution. These steps are under the control of a complex network of molecular factors, in which epigenetic mechanisms play a role that is increasingly well described. Recently, studies investigating epigenetic modifications and their impacts on gene expression in the mammary gland have been performed at different physiological stages and in different mammary cell types. This has led to the establishment of a role for epigenetic marks in milk component biosynthesis. This review aims to summarize the available knowledge regarding the involvement of the four main molecular mechanisms in epigenetics: DNA methylation, histone modifications, polycomb protein activity and non-coding RNA functions.

Combined Effects of Methylated Cytosine and Molecular Crowding on the Thermodynamic Stability of DNA Duplexes

Methylated cytosine within CpG dinucleotides is a key factor for epigenetic gene regulation. It has been revealed that methylated cytosine decreases DNA backbone flexibility and increases the thermal stability of DNA. Although the molecular environment is an important factor for the structure, thermodynamics, and function of biomolecules, there are few reports on the effects of methylated cytosine under a cell-mimicking molecular environment. Here, we systematically investigated the effects of methylated cytosine on the thermodynamics of DNA duplexes under molecular crowding conditions, which is a critical difference between the molecular environment in cells and test tubes. Thermodynamic parameters quantitatively demonstrated that the methylation effect and molecular crowding effect on DNA duplexes are independent and additive, in which the degree of the stabilization is the sum of the methylation effect and molecular crowding effect. Furthermore, the effects of methylation and molecular crowding correlate with the hydration states of DNA duplexes. The stabilization effect of methylation was due to the favorable enthalpic contribution, suggesting that direct interactions of the methyl group with adjacent bases and adjacent methyl groups play a role in determining the flexibility and thermodynamics of DNA duplexes. These results are useful to predict the properties of DNA duplexes with methylation in cell-mimicking conditions.

Hypermethylation of WNT3A gene and non-syndromic cleft lip and/or palate in association with in utero exposure to lead: A mediation analysis

OBJECTIVES

We aim to investigate association between WNT3A methylation and risk of non-syndromic cleft lip and/or palate (NSCL/P), and examine mediating effect of WNT3A methylation on the association of NSCL/P and lead (Pb) exposure in fetuses.

METHODS

DNA methylation of WNT3A in umbilical cord blood was determined among 59 NSCL/P cases and 118 non-malformed controls. Mediation analysis was performed to evaluate the potential mediating effect of WNT3A methylation on association between concentrations of Pb in umbilical cord and risk for NSCL/P. Additionally, an animal experiment in which cleft palates were induced by lead acetate was conducted.

RESULTS

The overall average methylation level of WNT3A was significant higher in NSCL/P cases as compared to controls. The risk for NSCL/P was increased by 1.90-fold with hypermethylation of WNT3A. Significant correlation was observed between concentrations of Pb in umbilical cord and methylation level of WNT3A. The hypermethylation of WNT3A had a mediating effect by 9.32% of total effect of Pb on NSCL/P risk. Gender-specific association between WNT3A methylation and NSCL/P was observed in male fetuses, and the percentage of the mediating effect increased to 14.28%. Animal experiment of mice showed that maternal oral exposure to lead acetate may result in cleft palate in offspring.

CONCLUSION

Hypermethylation of WNT3A was associated with the risk for NSCL/P and may be partly explain the association between exposure to Pb and risk for NSCL/P. The teratogenic and fetotoxic effects of Pb were found in mice.

TET-Mediated Epigenetic Regulation in Immune Cell Development and Disease

TET proteins oxidize 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) and further oxidation products in DNA. The oxidized methylcytosines (oxi-mCs) facilitate DNA demethylation and are also novel epigenetic marks. TET loss-of-function is strongly associated with cancer; TET2 loss-of-function mutations are frequently observed in hematological malignancies that are resistant to conventional therapies. Importantly, TET proteins govern cell fate decisions during development of various cell types by activating a cell-specific gene expression program. In this review, we seek to provide a conceptual framework of the mechanisms that fine tune TET activity. Then, we specifically focus on the multifaceted roles of TET proteins in regulating gene expression in immune cell development, function, and disease.

Distinctive epigenomic alterations in NF1-deficient cutaneous and plexiform neurofibromas drive differential MKK/p38 signaling

Benign peripheral nerve sheath tumors are the clinical hallmark of Neurofibromatosis Type 1. They account for substantial morbidity and mortality in NF1. Cutaneous (CNF) and plexiform neurofibromas (PNF) share nearly identical histology, but maintain different growth rates and risk of malignant conversion. The reasons for this disparate clinical behavior are not well explained by recent genome or transcriptome profiling studies. We hypothesized that CNFs and PNFs are epigenetically distinct tumor types that exhibit differential signaling due to genome-wide and site-specific methylation events. We interrogated the methylation profiles of 45 CNFs and 17 PNFs from NF1 subjects with the Illumina EPIC 850K methylation array. Based on these profiles, we confirm that CNFs and PNFs are epigenetically distinct tumors with broad differences in higher-order chromatin states and specific methylation events altering genes involved in key biological and cellular processes, such as inflammation, RAS/MAPK signaling, actin cytoskeleton rearrangement, and oxytocin signaling. Based on our identification of two separate DMRs associated with alternative leading exons in MAP2K3, we demonstrate differential RAS/MKK3/p38 signaling between CNFs and PNFs. Epigenetic reinforcement of RAS/MKK/p38 was a defining characteristic of CNFs leading to pro-inflammatory signaling and chromatin conformational changes, whereas PNFs signaled predominantly through RAS/MEK. Tumor size also correlated with specific CpG methylation events. Taken together, these findings confirm that NF1 deficiency influences the epigenetic regulation of RAS signaling fates, accounting for observed differences in CNF and PNF clinical behavior. The extension of these findings is that CNFs may respond differently than PNFs to RAS-targeted therapeutics raising the possibility of targeting p38-mediated inflammation for CNF treatment.

Deciphering the multifaceted roles of TET proteins in T-cell lineage specification and malignant transformation

TET proteins are DNA demethylases that can oxidize 5-methylcytosine (5mC) to generate 5-hydroxymethylcytosine (5hmC) and other oxidized mC bases (oxi-mCs). Importantly, TET proteins govern cell fate decisions during development of various cell types by activating a cell-specific gene expression program. In this review, we focus on the role of TET proteins in T-cell lineage specification. We explore the multifaceted roles of TET proteins in regulating gene expression in the contexts of T-cell development, lineage specification, function, and disease. Finally, we discuss the future directions and experimental strategies required to decipher the precise mechanisms employed by TET proteins to fine-tune gene expression and safeguard cell identity.

Modified Forms of Cytosine in Eukaryotes: DNA (De)methylation and Beyond

5-Methylcytosine (5mC) is an epigenetic mark known to contribute to the regulation of gene expression in a wide range of biological systems. Ten Eleven Translocation (TET) dioxygenases oxidize 5mC to 5-hydroxymethylcytosine, 5-formylcytosine, and 5-carboxylcytosine in metazoans and fungi. Moreover, two recent reports imply the existence of other species of modified cytosine in unicellular alga Chlamydomonas reinhardtii and malaria parasite Plasmodium falciparum. Here we provide an overview of the spectrum of cytosine modifications and their roles in demethylation of DNA and regulation of gene expression in different eukaryotic organisms.

DNA methylation changes in response to ocean acidification at the time of larval metamorphosis in the edible oyster, Crassostrea hongkongensis

Unprecedented rate of increased CO₂ level in the ocean and the subsequent changes in carbonate system including decreased pH, known as ocean acidification (OA), is predicted to disrupt not only the calcification process but also several other physiological and developmental processes in a variety of marine organisms, including edible oysters. Nonetheless, not all species are vulnerable to those OA threats, e.g. some species may be able to cope with OA stress using environmentally induced modifications on gene and protein expressions. For example, external environmental stressors including OA can influence the addition and removal of methyl groups through epigenetic modification (e.g. DNA methylation) process to turn gene expression "on or off" as part of a rapid adaptive mechanism to cope with OA. In this study, we tested the above hypothesis through testing the effect of OA, using decreased pH 7.4 as proxy, on DNA methylation pattern of an endemic and a commercially important estuary oyster species, Crassostrea hongkongensis at the time of larval habitat selection and metamorphosis. Larval growth rate did not differ between control pH 8.1 and treatment pH 7.4. The metamorphosis rate of the pediveliger larvae was higher at pH 7.4 than those in control pH 8.1, however over one-third of the larvae raised at pH 7.4 failed to attach on optimal substrate as defined by biofilm presence. During larval development, a total of 130 genes were differentially methylated across the two treatments. The differential methylation in the larval genes may have partially accounted for the higher metamorphosis success rate under decreased pH 7.4 but with poor substratum selection ability. Differentially methylated loci were concentrated in the exon regions and appear to be associated with cytoskeletal and signal transduction, oxidative stress, metabolic processes, and larval metamorphosis, which implies the high potential of C. hongkongensis larvae to acclimate and adapt through non-genetic ways to OA threats within a single generation.

DNA methylation changes in response to ocean acidification at the time of larval metamorphosis in the edible oyster, Crassostrea hongkongensis

Unprecedented rate of increased CO₂ level in the ocean and the subsequent changes in carbonate system including decreased pH, known as ocean acidification (OA), is predicted to disrupt not only the calcification process but also several other physiological and developmental processes in a variety of marine organisms, including edible oysters. Nonetheless, not all species are vulnerable to those OA threats, e.g. some species may be able to cope with OA stress using environmentally induced modifications on gene and protein expressions. For example, external environmental stressors including OA can influence the addition and removal of methyl groups through epigenetic modification (e.g. DNA methylation) process to turn gene expression "on or off" as part of a rapid adaptive mechanism to cope with OA. In this study, we tested the above hypothesis through testing the effect of OA, using decreased pH 7.4 as proxy, on DNA methylation pattern of an endemic and a commercially important estuary oyster species, Crassostrea hongkongensis at the time of larval habitat selection and metamorphosis. Larval growth rate did not differ between control pH 8.1 and treatment pH 7.4. The metamorphosis rate of the pediveliger larvae was higher at pH 7.4 than those in control pH 8.1, however over one-third of the larvae raised at pH 7.4 failed to attach on optimal substrate as defined by biofilm presence. During larval development, a total of 130 genes were differentially methylated across the two treatments. The differential methylation in the larval genes may have partially accounted for the higher metamorphosis success rate under decreased pH 7.4 but with poor substratum selection ability. Differentially methylated loci were concentrated in the exon regions and appear to be associated with cytoskeletal and signal transduction, oxidative stress, metabolic processes, and larval metamorphosis, which implies the high potential of C. hongkongensis larvae to acclimate and adapt through non-genetic ways to OA threats within a single generation.

Cortical cellular diversity and development in schizophrenia

While a definitive understanding of schizophrenia etiology is far from current reality, an increasing body of evidence implicates perturbations in early development that alter the trajectory of brain maturation in this disorder, leading to abnormal function in early childhood and adulthood. This atypical development likely arises from an interaction of many brain cell types that follow distinct developmental paths. Because both cellular identity and development are governed by the transcriptome and epigenome, two levels of gene regulation that have the potential to reflect both genetic and environmental influences, mapping "omic" changes over development in diverse cells is a fruitful avenue for schizophrenia research. In this review, we provide a survey of human brain cellular composition and development, levels of genomic regulation that determine cellular identity and developmental trajectories, and what is known about how genomic regulation is dysregulated in specific cell types in schizophrenia. We also outline technical challenges and solutions to conducting cell type-specific functional genomic studies in human postmortem brain.

Epigenetic interaction of microbes with their mammalian hosts

The interaction of microbiota with its host has the ability to alter the cellular functions of both, through several mechanisms. Recent work, from many laboratories including our own, has shown that epigenetic mechanisms play an important role in the alteration of these cellular functions. Epigenetics broadly refers to change in the phenotype without a corresponding change in the DNA sequence. This change is usually brought by epigenetic modifications of the DNA itself, the histone proteins associated with the DNA in the chromatin, non-coding RNA or the modifications of the transcribed RNA. These modifications, also known as epigenetic code, do not change the DNA sequence but alter the expression level of specific genes. Microorganisms seem to have learned how to modify the host epigenetic code and modulate the host transcriptome in their favour. In this review, we explore the literature that describes the epigenetic interaction of bacteria, fungi and viruses, with their mammalian hosts.

Regulation of ST6GAL1 sialyltransferase expression in cancer cells

The ST6GAL1 sialyltransferase, which adds α2-6 linked sialic acids to N-glycosylated proteins, is overexpressed in a wide range of human malignancies. Recent studies have established the importance of ST6GAL1 in promoting tumor cell behaviors such as invasion, resistance to cell stress and chemoresistance. Furthermore, ST6GAL1 activity has been implicated in imparting cancer stem cell characteristics. However, despite the burgeoning interest in the role of ST6GAL1 in the phenotypic features of tumor cells, insufficient attention has been paid to the molecular mechanisms responsible for ST6GAL1 upregulation during neoplastic transformation. Evidence suggests that these mechanisms are multifactorial, encompassing genetic, epigenetic, transcriptional and posttranslational regulation. The purpose of this review is to summarize current knowledge regarding the molecular events that drive enriched ST6GAL1 expression in cancer cells.

Chromatin Structure and Function in Mosquitoes

The principles and function of chromatin and nuclear architecture have been extensively studied in model organisms, such as Drosophila melanogaster. However, little is known about the role of these epigenetic processes in transcriptional regulation in other insects including mosquitoes, which are major disease vectors and a worldwide threat for human health. Some of these life-threatening diseases are malaria, which is caused by protozoan parasites of the genus Plasmodium and transmitted by Anopheles mosquitoes; dengue fever, which is caused by an arbovirus mainly transmitted by Aedes aegypti; and West Nile fever, which is caused by an arbovirus transmitted by Culex spp. In this contribution, we review what is known about chromatin-associated mechanisms and the 3D genome structure in various mosquito vectors, including Anopheles, Aedes, and Culex spp. We also discuss the similarities between epigenetic mechanisms in mosquitoes and the model organism Drosophila melanogaster, and advocate that the field could benefit from the cross-application of state-of-the-art functional genomic technologies that are well-developed in the fruit fly. Uncovering the mosquito regulatory genome can lead to the discovery of unique regulatory networks associated with the parasitic life-style of these insects. It is also critical to understand the molecular interactions between the vectors and the pathogens that they transmit, which could hold the key to major breakthroughs on the fight against mosquito-borne diseases. Finally, it is clear that epigenetic mechanisms controlling mosquito environmental plasticity and evolvability are also of utmost importance, particularly in the current context of globalization and climate change.

A novel computational strategy for DNA methylation imputation using mixture regression model (MRM)

BACKGROUND

DNA methylation is an important heritable epigenetic mark that plays a crucial role in transcriptional regulation and the pathogenesis of various human disorders. The commonly used DNA methylation measurement approaches, e.g., Illumina Infinium HumanMethylation-27 and -450 BeadChip arrays (27 K and 450 K arrays) and reduced representation bisulfite sequencing (RRBS), only cover a small proportion of the total CpG sites in the human genome, which considerably limited the scope of the DNA methylation analysis in those studies.

RESULTS

We proposed a new computational strategy to impute the methylation value at the unmeasured CpG sites using the mixture of regression model (MRM) of radial basis functions, integrating information of neighboring CpGs and the similarities in local methylation patterns across subjects and across multiple genomic regions. Our method achieved a better imputation accuracy over a set of competing methods on both simulated and empirical data, particularly when the missing rate is high. By applying MRM to an RRBS dataset from subjects with low versus high bone mineral density (BMD), we recovered methylation values of ~ 300 K CpGs in the promoter regions of chromosome 17 and identified some novel differentially methylated CpGs that are significantly associated with BMD.

CONCLUSIONS

Our method is well applicable to the numerous methylation studies. By expanding the coverage of the methylation dataset to unmeasured sites, it can significantly enhance the discovery of novel differential methylation signals and thus reveal the mechanisms underlying various human disorders/traits.

Complex DNA sequence readout mechanisms of the DNMT3B DNA methyltransferase

DNA methyltransferases interact with their CpG target sites in the context of variable flanking sequences. We investigated DNA methylation by the human DNMT3B catalytic domain using substrate pools containing CpX target sites in randomized flanking context and identified combined effects of CpG recognition and flanking sequence interaction together with complex contact networks involved in balancing the interaction with different flanking sites. DNA methylation rates were more affected by flanking sequences at non-CpG than at CpG sites. We show that T775 has an essential dynamic role in the catalytic mechanism of DNMT3B. Moreover, we identify six amino acid residues in the DNA-binding interface of DNMT3B (N652, N656, N658, K777, N779, and R823), which are involved in the equalization of methylation rates of CpG sites in favored and disfavored sequence contexts by forming compensatory interactions to the flanking residues including a CpG specific contact to an A at the +1 flanking site. Non-CpG flanking preferences of DNMT3B are highly correlated with non-CpG methylation patterns in human cells. Comparison of the flanking sequence preferences of human and mouse DNMT3B revealed subtle differences suggesting a co-evolution of flanking sequence preferences and cellular DNMT targets.

HIF3A DNA methylation, obesity and weight gain, and breast cancer risk among Mexican American women

OBJECTIVE

In previous epigenome-wide association studies, Hypoxia inducible Factor 3 Alpha Subunit (HIF3A) DNA methylation has been reported to be associated with body mass index (BMI) and weight change. However, none of these studies have included Mexican Americans.

METHODS

In the current study, we assessed levels of HIF3A methylation in 927 Mexican American women identified from Mano-A-Mano, the Mexican American Cohort study.

RESULTS

Significantly higher methylation levels at three CpG sites (position 46801557, 46801642, and 46801699) were observed in obese women compared to non-obese women (P < 0.05). Furthermore, we found that elevated methylation levels at those three CpG sites were associated with significant weight gain (P < 0.05), defined as an increase in BMI by at least one category between the baseline and the follow-up, with a median follow-up time of 39 months. Then, using pre-diagnostic blood DNA samples, we found increased DNA methylation at CpG 46801642 to be associated with a 1.35-fold increased risk of breast cancer (Hazard Ratio (HR) = 1.35, 95% Confidence Interval (CI): 1.02, 3.01), with a median follow-up time of 127 months. Using the Cancer Genome Atlas (TCGA) data, we further found that levels of HIF3A were significantly higher-methylated and down-regulated in breast tumor than in normal tissues (P < 1 × 10¹² for both).

CONCLUSION

Thus, our results provide evidence to support the role of HIF3A in obesity, weight gain, and the development of breast cancer.

MethHaplo: combining allele-specific DNA methylation and SNPs for haplotype region identification

Background

DNA methylation is an important epigenetic modification that plays a critical role in most eukaryotic organisms. Parental alleles in haploid genomes may exhibit different methylation patterns, which can lead to different phenotypes and even different therapeutic and drug responses to diseases. However, to our knowledge, no software is available for the identification of DNA methylation haplotype regions with combined allele-specific DNA methylation, single nucleotide polymorphisms (SNPs) and high-throughput chromosome conformation capture (Hi-C) data.

Results

In this paper, we developed a new method, MethHaplo, that identify DNA methylation haplotype regions with allele-specific DNA methylation and SNPs from whole-genome bisulfite sequencing (WGBS) data. Our results showed that methylation haplotype regions were ten times longer than haplotypes with SNPs only. When we integrate WGBS and Hi-C data, MethHaplo could call even longer haplotypes.

Conclusions

This study illustrates the usefulness of methylation haplotypes. By constructing methylation haplotypes for various cell lines, we provide a clearer picture of the effect of DNA methylation on gene expression, histone modification and three-dimensional chromosome structure at the haplotype level. Our method could benefit the study of parental inheritance-related disease and hybrid vigor in agriculture.

DNA methylation methods: Global DNA methylation and methylomic analyses

DNA methylation provides a pivotal layer of epigenetic regulation in eukaryotes that has significant involvement for numerous biological processes in health and disease. The function of methylation of cytosine bases in DNA was originally proposed as a "silencing" epigenetic marker and focused on promoter regions of genes for decades. Improved technologies and accumulating studies have been extending our understanding of the roles of DNA methylation to various genomic contexts including gene bodies, repeat sequences and transcriptional start sites. The demand for comprehensively describing DNA methylation patterns spawns a diversity of DNA methylation profiling technologies that target its genomic distribution. These approaches have enabled the measurement of cytosine methylation from specific loci at restricted regions to single-base-pair resolution on a genome-scale level. In this review, we discuss the different DNA methylation analysis technologies primarily based on the initial treatments of DNA samples: bisulfite conversion, endonuclease digestion and affinity enrichment, involving methodology evolution, principles, applications, and their relative merits. This review may offer referable information for the selection of various platforms for genome-wide analysis of DNA methylation.

Targeted bisulfite sequencing: A novel tool for the assessment of DNA methylation with high sensitivity and increased coverage

DNA methylation analysis is increasingly used in stress research. Available methods are expensive, laborious and often limited by either the analysis of short CpG stretches or low assay sensitivity. Here, we present a cost-efficient next generation sequencing-based strategy for the simultaneous investigation of multiple candidate genes in large cohorts. To illustrate the method, we present analysis of four candidate genes commonly assessed in psychoneuroendocrine research: Glucocorticoid receptor (NR3C1), Serotonin transporter (SLC6A4), FKBP Prolyl isomerase 5 (FKBP5), and the Oxytocin receptor (OXTR). DNA methylation standards (100 %; 75 %; 50 %; 25 % and 0 %) and DNA of a female and male donor were bisulfite treated in three independent trials and were used to generate sequencing libraries for 42 CpGs from the NR3C1 1 F promoter region, 84 CpGs of the SLC6A4 5' regulatory region, 5 CpGs located in FKBP5 intron 7, and additional 12 CpGs located in a potential enhancer element in intron 3 of the OXTR. In addition, DNA of 45 patients with borderline personality disorder (BPD) and 45 healthy controls was assayed. Multiplex libraries of all samples were sequenced on a MiSeq system and analyzed for mean methylation values of all CpG sites using amplikyzer2 software. Results indicated excellent accuracy of the assays when investigating replicates generated from the same bisulfite converted DNA, and very high linearity (R² > 0.9) of the assays shown by the analysis of differentially methylated DNA standards. Comparing DNA methylation between BPD and healthy controls revealed no biologically relevant differences. The technical approach as described here facilitates targeted DNA methylation analysis and represents a highly sensitive, cost-efficient and high throughput tool to close the gap between coverage and precision in epigenetic research of stress-associated phenotypes.

DNA Methylation as a Potential Molecular Mechanism in X-linked Dystonia-Parkinsonism

BACKGROUND

X-linked dystonia-parkinsonism is a neurodegenerative movement disorder. The underlying molecular basis has still not been completely elucidated, but likely involves dysregulation of TAF1 expression. In X-linked dystonia-parkinsonism, 3 disease-specific single-nucleotide changes (DSCs) introduce (DSC12) or abolish (DSC2 and DSC3) CpG dinucleotides and consequently sites of putative DNA methylation. Because transcriptional regulation tightly correlates with specific epigenetic marks, we investigated the role of DNA methylation in the pathogenesis of X-linked dystonia-parkinsonism.

METHODS

DNA methylation at DSC12, DSC3, and DSC2 was quantified by bisulfite pyrosequencing in DNA from peripheral blood leukocytes, fibroblasts, induced pluripotent stem cell-derived cortical neurons and brain tissue from X-linked dystonia-parkinsonism patients and age- and sex-matched healthy Filipino controls in a prospective study.

RESULTS

Compared with controls, X-linked dystonia-parkinsonism patients showed striking differences in DNA methylation at the 3 investigated CpG sites. Using methylation-sensitive luciferase reporter gene assays and immunoprecipitation, we demonstrated (1) that lack of DNA methylation because of DSC2 and DSC3 affects gene promoter activity and (2) that methylation at all 3 investigated CpG sites alters DNA-protein interaction. Interestingly, DSC3 decreased promoter activity per se compared with wild type, and promoter activity further decreased when methylation was present. Moreover, we identified specific binding of proteins to the investigated DSCs that are associated with splicing and RNA and DNA binding.

CONCLUSIONS

We identified altered DNA methylation in X-linked dystonia-parkinsonism patients as a possible additional mechanism modulating TAF1 expression and putative novel targets for future therapies using DNA methylation-modifying agents. © 2020 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Gestational Diabetes Mellitus Affects Offspring's Epigenome. Is There a Way to Reduce the Negative Consequences?

Gestational diabetes mellitus (GDM) is the most common pregnancy complication worldwide and may result in short-term and long-term consequences for offspring. The present review highlights evidence of epigenetic programming, mostly from human studies, which occurs in offspring exposed to maternal GDM during different stages of development, paying special attention to the differences in sensitivity of offspring to maternal hyperglycemia as a result of sex-related factors. We also aim to answer the following question: If these epigenetic changes are constant throughout the lifetime of the offspring, how do they present phenotypically?

Comparative DNA methylome analysis of estrus ewes reveals the complex regulatory pathways of sheep fecundity

BACKGROUND/AIMS

Sheep are important livestock with variant ovulation rate and fertility. Dorset sheep is a typical breed with low prolificacy, whereas Small Tail Han sheep with FecB mutation (HanBB) have hyperprolificacy. Our previous studies have revealed the gene expression difference between the ovaries from Dorset and HanBB sheep contributes to the difference of fecundity, however, what leads to these gene expression difference remains unclear. DNA methylation, an important epigenetic process, plays a crucial role in gene expression regulation.

METHODS

In the present study, we constructed a methylated DNA immunoprecipitation combined with high throughput sequencing (MeDIP-seq) strategy to investigate the differentially methylated genes between the Dorset and HanBB ovaries.

RESULTS

Our findings suggest the genes involved in immune response, branched-chain amino acid metabolism, cell growth and cell junction were differentially methylated in or around the gene body regions.

CONCLUSIONS

These findings provide prospective insights on the epigenetic basis of sheep fecundity.

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.

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.

Comprehensive structure-function characterization of DNMT3B and DNMT3A reveals distinctive de novo DNA methylation mechanisms

Mammalian DNA methylation patterns are established by two de novo DNA methyltransferases, DNMT3A and DNMT3B, which exhibit both redundant and distinctive methylation activities. However, the related molecular basis remains undetermined. Through comprehensive structural, enzymology and cellular characterization of DNMT3A and DNMT3B, we here report a multi-layered substrate-recognition mechanism underpinning their divergent genomic methylation activities. A hydrogen bond in the catalytic loop of DNMT3B causes a lower CpG specificity than DNMT3A, while the interplay of target recognition domain and homodimeric interface fine-tunes the distinct target selection between the two enzymes, with Lysine 777 of DNMT3B acting as a unique sensor of the +1 flanking base. The divergent substrate preference between DNMT3A and DNMT3B provides an explanation for site-specific epigenomic alterations seen in ICF syndrome with DNMT3B mutations. Together, this study reveals distinctive substrate-readout mechanisms of the two DNMT3 enzymes, implicative of their differential roles during development and pathogenesis.

DeepH&M: Estimating single-CpG hydroxymethylation and methylation levels from enrichment and restriction enzyme sequencing methods

Increased appreciation of 5-hydroxymethylcytosine (5hmC) as a stable epigenetic mark, which defines cell identity and disease progress, has engendered a need for cost-effective, but high-resolution, 5hmC mapping technology. Current enrichment-based technologies provide cheap but low-resolution and relative enrichment of 5hmC levels, while single-base resolution methods can be prohibitively expensive to scale up to large experiments. To address this problem, we developed a deep learning-based method, "DeepH&M," which integrates enrichment and restriction enzyme sequencing methods to simultaneously estimate absolute hydroxymethylation and methylation levels at single-CpG resolution. Using 7-week-old mouse cerebellum data for training the DeepH&M model, we demonstrated that the 5hmC and 5mC levels predicted by DeepH&M were in high concordance with whole-genome bisulfite-based approaches. The DeepH&M model can be applied to 7-week-old frontal cortex and 79-week-old cerebellum, revealing the robust generalizability of this method to other tissues from various biological time points.

Methylation-Based Therapies for Colorectal Cancer

Colorectal carcinogenesis (CRC) is caused by the gradual long-term accumulation of both genetic and epigenetic changes. Recently, epigenetic alterations have been included in the classification of the CRC molecular subtype, and this points out their prognostic impact. As epigenetic modifications are reversible, they may represent relevant therapeutic targets. DNA methylation, catalyzed by DNA methyltransferases (DNMTs), regulates gene expression. For many years, the deregulation of DNA methylation has been considered to play a substantial part in CRC etiology and evolution. Despite considerable advances in CRC treatment, patient therapy response persists as limited, and their profit from systemic therapies are often hampered by the introduction of chemoresistance. In addition, inter-individual changes in therapy response in CRC patients can arise from their specific (epi)genetic compositions. In this review article, we summarize the options of CRC treatment based on DNA methylation status for their predictive value. This review also includes the therapy outcomes based on the patient's methylation status in CRC patients. In addition, the current challenge of research is to develop therapeutic inhibitors of DNMT. Based on the essential role of DNA methylation in CRC development, the application of DNMT inhibitors was recently proposed for the treatment of CRC patients, especially in patients with DNA hypermethylation.

Profiling the DNA methylation patterns of imprinted genes in abnormal semen samples by next-generation bisulfite sequencing

PURPOSE

Changes in DNA methylation modifications have been associated with male infertility. With the development of assisted reproductive technologies (ARTs), abnormal DNA methylation in sperm, especially in imprinted genes, may impact the health of offspring and requires an in-depth study.

METHODS

In this study, we collected abnormal human semen samples, including asthenospermic, oligospermic, oligoasthenospermic and deformed sperm, and investigated the methylation of imprinted genes by reduced representation bisulfite sequencing (RRBS) and bisulfite amplicon sequencing on the Illumina platform.

RESULTS

The differentially methylated regions (DMRs) of imprinted genes, including H19, GNAS, MEG8 and SNRPN, were different in the abnormal semen groups. MEG8 DMR methylation in the asthenospermic group was significantly increased. Furthermore, higher methylation levels of MEG8, GNAS and SNRPN DMR in the oligospermic and oligoasthenospermic groups and a decrease in the H19 DMR methylation level in the oligospermic group were observed. However, the methylation levels of these regions varied greatly among the different semen samples and among individual sperm within the same semen sample. The SNP rs2525883 genotype in the H19 DMR affected DNA methylation. Moreover, DNA methylation levels differed in the abnormal semen groups in the non-imprinted genomic regions, including repetitive sequence DNA transposons and long/short interspersed nuclear elements (LINEs and SINEs).

CONCLUSION

Our study established that imprinted gene DMRs, such as H19, GNAS, SNRPN and MEG8, were differentially methylated in the abnormal semen groups with obvious inter- and intra-sample heterogeneities. These results suggest that special attention needs to be paid to possible epigenetic risks during reproduction.

A cost-effective approach to DNA methylation detection by Methyl Sensitive DArT sequencing

Several studies suggest the relation of DNA methylation to diseases in humans and important phenotypes in plants drawing attention to this epigenetic mark as an important source of variability. In the last decades, several methodologies were developed to assess the methylation state of a genome. However, there is still a lack of affordable and precise methods for genome wide analysis in large sample size studies. Methyl sensitive double digestion MS-DArT sequencing method emerges as a promising alternative for methylation profiling. We developed a computational pipeline for the identification of DNA methylation using MS-DArT-seq data and carried out a pilot study using the Eucalyptus grandis tree sequenced for the species reference genome. Using a statistic framework as in differential expression analysis, 72,515 genomic sites were investigated and 5,846 methylated sites identified, several tissue specific, distributed along the species 11 chromosomes. We highlight a bias towards identification of DNA methylation in genic regions and the identification of 2,783 genes and 842 transposons containing methylated sites. Comparison with WGBS, DNA sequencing after treatment with bisulfite, data demonstrated a precision rate higher than 95% for our approach. The availability of a reference genome is useful for determining the genomic context of methylated sites but not imperative, making this approach suitable for any species. Our approach provides a cost effective, broad and reliable examination of DNA methylation profile on MspI/HpaII restriction sites, is fully reproducible and the source code is available on GitHub (https://github.com/wendelljpereira/ms-dart-seq).

A cost-effective approach to DNA methylation detection by Methyl Sensitive DArT sequencing

Several studies suggest the relation of DNA methylation to diseases in humans and important phenotypes in plants drawing attention to this epigenetic mark as an important source of variability. In the last decades, several methodologies were developed to assess the methylation state of a genome. However, there is still a lack of affordable and precise methods for genome wide analysis in large sample size studies. Methyl sensitive double digestion MS-DArT sequencing method emerges as a promising alternative for methylation profiling. We developed a computational pipeline for the identification of DNA methylation using MS-DArT-seq data and carried out a pilot study using the Eucalyptus grandis tree sequenced for the species reference genome. Using a statistic framework as in differential expression analysis, 72,515 genomic sites were investigated and 5,846 methylated sites identified, several tissue specific, distributed along the species 11 chromosomes. We highlight a bias towards identification of DNA methylation in genic regions and the identification of 2,783 genes and 842 transposons containing methylated sites. Comparison with WGBS, DNA sequencing after treatment with bisulfite, data demonstrated a precision rate higher than 95% for our approach. The availability of a reference genome is useful for determining the genomic context of methylated sites but not imperative, making this approach suitable for any species. Our approach provides a cost effective, broad and reliable examination of DNA methylation profile on MspI/HpaII restriction sites, is fully reproducible and the source code is available on GitHub (https://github.com/wendelljpereira/ms-dart-seq).

The Epigenetic Progenitor Origin of Cancer Reassessed: DNA Methylation Brings Balance to the Stem Force

Cancer initiation and progression toward malignant stages occur as the results of accumulating genetic alterations and epigenetic dysregulation. During the last decade, the development of next generation sequencing (NGS) technologies and the increasing pan-genomic knowledge have revolutionized how we consider the evolving epigenetic landscapes during homeostasis and tumor progression. DNA methylation represents the best studied mark and is considered as a common mechanism of epigenetic regulation in normal homeostasis and cancer. A remarkable amount of work has recently started clarifying the central role played by DNA methylation dynamics on the maintenance of cell identity and on cell fate decisions during the different steps of normal development and tumor evolution. Importantly, a growing number of studies show that DNA methylation is key in the maintenance of adult stemness and in orchestrating commitment in multiple ways. Perturbations of the normal DNA methylation patterns impair the homeostatic balance and can lead to tumor initiation. Therefore, DNA methylation represents an interesting therapeutic target to recover homeostasis in tumor stem cells.

A Conformational Switch in the Zinc Finger Protein Kaiso Mediates Differential Readout of Specific and Methylated DNA Sequences

Recognition of the epigenetic mark 5-methylcytosine (mC) at CpG sites in DNA has emerged as a novel function of many eukaryotic transcription factors (TFs). It remains unclear why the sequence specificity of these TFs differs for CpG-methylated motifs and consensus motifs. Here, we dissect the structural and dynamic basis for this differential DNA binding specificity in the human zinc finger TF Kaiso, which exhibits high affinity for two consecutive mCpG sites in variable contexts and also for a longer, sequence-specific Kaiso binding site (KBS). By integrating structural analysis and DNA binding studies with targeted protein mutagenesis and nucleotide substitutions, we identify distinct mechanisms for readout of methylated and KBS motifs by Kaiso. We show that a key glutamate residue (E535), critical for mCpG site recognition, adopts different conformations in complexes with specific and methylated DNA. These conformational differences, together with intrinsic variations in DNA flexibility and/or solvation at TpG versus mCpG sites, contribute to the different DNA affinity and sequence specificity. With methylated DNA, multiple direct contacts between E535 and the 5' mCpG site dominate the binding affinity, allowing for tolerance of different flanking DNA sequences. With KBS, Kaiso employs E535 as part of an indirect screen of the 5' flanking sequence, relying on key tyrosine-DNA interactions to stabilize an optimal DNA conformation and select against noncognate sites. These findings demonstrate how TFs use conformational adaptation and exploit variations in DNA flexibility to achieve distinct DNA readout outcomes and target a greater variety of regulatory and epigenetic sites than previously appreciated.

Mitochondrial DNA alterations may influence the cisplatin responsiveness of oral squamous cell carcinoma

Cisplatin is the first-line chemotherapeutic agent for the treatment of oral squamous cell carcinoma (OSCC). However, the intrinsic or acquired resistance against cisplatin remains a major obstacle to treatment efficacy in OSCC. Recently, mitochondrial DNA (mtDNA) alterations have been reported in a variety of cancers. However, the role of mtDNA alterations in OSCC has not been comprehensively studied. In this study, we evaluated the correlation between mtDNA alterations (mtDNA content, point mutations, large-scale deletions, and methylation status) and cisplatin sensitivity using two OSCC cell lines, namely SAS and H103, and stem cell-like tumour spheres derived from SAS. By microarray analysis, we found that the tumour spheres profited from aberrant lipid and glucose metabolism and became resistant to cisplatin. By qPCR analysis, we found that the cells with less mtDNA were less responsive to cisplatin (H103 and the tumour spheres). Based on the findings, we theorised that the metabolic changes in the tumour spheres probably resulted in mtDNA depletion, as the cells suppressed mitochondrial respiration and switched to an alternative mode of energy production, i.e. glycolysis. Then, to ascertain the origin of the variation in mtDNA content, we used MinION, a nanopore sequencer, to sequence the mitochondrial genomes of H103, SAS, and the tumour spheres. We found that the lower cisplatin sensitivity of H103 could have been caused by a constellation of genetic and epigenetic changes in its mitochondrial genome. Future work may look into how changes in mtDNA translate into an impact on cell function and therefore cisplatin response.

Prediction model construction of mouse stem cell pluripotency using CpG and non-CpG DNA methylation markers

Background

Genome-wide studies of DNA methylation across the epigenetic landscape provide insights into the heterogeneity of pluripotent embryonic stem cells (ESCs). Differentiating into embryonic somatic and germ cells, ESCs exhibit varying degrees of pluripotency, and epigenetic changes occurring in this process have emerged as important factors explaining stem cell pluripotency.

Results

Here, using paired scBS-seq and scRNA-seq data of mice, we constructed a machine learning model that predicts degrees of pluripotency for mouse ESCs. Since the biological activities of non-CpG markers have yet to be clarified, we tested the predictive power of CpG and non-CpG markers, as well as a combination thereof, in the model. Through rigorous performance evaluation with both internal and external validation, we discovered that a model using both CpG and non-CpG markers predicted the pluripotency of ESCs with the highest prediction performance (0.956 AUC, external test). The prediction model consisted of 16 CpG and 33 non-CpG markers. The CpG and most of the non-CpG markers targeted depletions of methylation and were indicative of cell pluripotency, whereas only a few non-CpG markers reflected accumulations of methylation. Additionally, we confirmed that there exists the differing pluripotency between individual developmental stages, such as E3.5 and E6.5, as well as between induced mouse pluripotent stem cell (iPSC) and somatic cell.

Conclusions

In this study, we investigated CpG and non-CpG methylation in relation to mouse stem cell pluripotency and developed a model thereon that successfully predicts the pluripotency of mouse ESCs.

Towards a comprehensive catalogue of validated and target-linked human enhancers

The human gene catalogue is essentially complete, but we lack an equivalently vetted inventory of bona fide human enhancers. Hundreds of thousands of candidate enhancers have been nominated via biochemical annotations; however, only a handful of these have been validated and confidently linked to their target genes. Here we review emerging technologies for discovering, characterizing and validating human enhancers at scale. We furthermore propose a new framework for operationally defining enhancers that accommodates the heterogeneous and complementary results that are emerging from reporter assays, biochemical measurements and CRISPR screens.

Comparative Analysis of CpG Sites and Islands Distributed in Mitochondrial DNA of Model Organisms

Simple Summary

In recent years, the existence of methylation of mammalian mitochondrial DNA (mtDNA) has been discussed. The current state of knowledge indicates that mtDNA is poorly methylated; in fact, it only accounts for 2–8% methylated sites and its pattern is unknown. The lack of comprehensive information on the mtDNA methylation pattern prompted us to investigate the distribution of guanine-cytosine-rich sequences (CpG) in different animal species. The aim of the study was to determine the localization of CpG sites and islands in mtDNA of model organisms. The CpG sites and islands found in vertebrates and invertebrates indicate a diversified pattern of CpG distribution. Generally, the number of observed CpG sites of the mitochondrial genome was higher in the analysed vertebrates than in the invertebrates. However, there was no relationship between the frequency of the CpG sites in the mitochondrial genome and the complexity of the analysed organism. The distribution of the CpG sites for transfer RNA (tRNA) coding genes was usually cumulated in a larger CpG region in the vertebrates.

Abstract

The information about mtDNA methylation is still limited, thus epigenetic modification remains unclear. The lack of comprehensive information on the comparative epigenomics of mtDNA prompts comprehensive investigations of the epigenomic modification of mtDNA in different species. This is the first study in which the theoretical CpG localization in the mtDNA reference sequences from various species (12) was compared. The aim of the study was to determine the localization of CpG sites and islands in mtDNA of model organisms and to compare their distribution. The results are suitable for further investigations of mtDNA methylation. The analysis involved both strands of mtDNA sequences of animal model organisms representing different taxonomic groups of invertebrates and vertebrates. For each sequence, such parameters as the number, length, and localization of CpG islands were determined with the use of EMBOSS (European Molecular Biology Open Software Suite) software. The number of CpG sites for each sequence was indicated using the newcpgseek algorithm. The results showed that methylation of mtDNA in the analysed species involved mitochondrial gene expression. Our analyses showed that the CpG sites were commonly present in genomic regions including the D-loop, CYTB, ND6, ND5, ND4, ND3, ND2, ND1, COX3, COX2, COX1, ATP6, 16s rRNA, and 12s rRNA. The CpG distribution in animals from different species was diversified. Generally, the number of observed CpG sites of the mitochondrial genome was higher in the vertebrates than in the invertebrates. However, there was no relationship between the frequency of the CpG sites in the mitochondrial genome and the complexity of the analysed organisms. Interestingly, the distribution of the CpG sites for tRNA coding genes was usually cumulated in a larger CpG region in vertebrates. This paper may be a starting point for further research, since the collected information indicates possible methylation regions localized in mtDNA among different species including invertebrates and vertebrates.

Precise genomic mapping of 5-hydroxymethylcytosine via covalent tether-directed sequencing

5-hydroxymethylcytosine (5hmC) is the most prevalent intermediate on the oxidative DNA demethylation pathway and is implicated in regulation of embryogenesis, neurological processes, and cancerogenesis. Profiling of this relatively scarce genomic modification in clinical samples requires cost-effective high-resolution techniques that avoid harsh chemical treatment. Here, we present a bisulfite-free approach for 5hmC profiling at single-nucleotide resolution, named hmTOP-seq (5hmC-specific tethered oligonucleotide–primed sequencing), which is based on direct sequence readout primed at covalently labeled 5hmC sites from an in situ tethered DNA oligonucleotide. Examination of distinct conjugation chemistries suggested a structural model for the tether-directed nonhomologous polymerase priming enabling theoretical evaluation of suitable tethers at the design stage. The hmTOP-seq procedure was optimized and validated on a small model genome and mouse embryonic stem cells, which allowed construction of single-nucleotide 5hmC maps reflecting subtle differences in strand-specific CG hydroxymethylation. Collectively, hmTOP-seq provides a new valuable tool for cost-effective and precise identification of 5hmC in characterizing its biological role and epigenetic changes associated with human disease.

This study describes hmTOP-seq, a bisulfite-free approach for profiling of the epigenetic mark 5-hydroxymethylcytosine (5hmC) at single-nucleotide resolution, based on direct sequence readout primed at an in situ tethered DNA oligonucleotide.

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.

A concise review of human brain methylome during aging and neurodegenerative diseases

DNA methylation at CpG sites is an essential epigenetic mark that regulates gene expression during mammalian development and diseases. Methylome refers to the entire set of methylation modifications present in the whole genome. Over the last several years, an increasing number of reports on brain DNA methylome reported the association between aberrant methylation and the abnormalities in the expression of critical genes known to have critical roles during aging and neurodegenerative diseases. Consequently, the role of methylation in understanding neurodegenerative diseases has been under focus. This review outlines the current knowledge of the human brain DNA methylomes during aging and neurodegenerative diseases. We describe the differentially methylated genes from fetal stage to old age and their biological functions. Additionally, we summarize the key aspects and methylated genes identified from brain methylome studies on neurodegenerative diseases. The brain methylome studies could provide a basis for studying the functional aspects of neurodegenerative diseases.

Promoter DNA Hypermethylation and Paradoxical Gene Activation

DNA methylation is a stable epigenetic modification that contributes to the spatiotemporal regulation of gene expression. The manner in which DNA methylation contributes to transcriptional control is dependent on the biological context, including physiological state and the properties of the DNA itself. Classically, dense promoter DNA methylation is associated with transcriptional repression. However, growing evidence suggests that this association may not always hold true, and promoter hypermethylation now also appears to be associated with high transcriptional activity. Furthermore, in a selection of contexts, increasing levels of promoter methylation correlate directly with increased gene expression. These findings postulate a context-dependent model whereby epigenetic contributions to transcriptional regulation occur in a more complex and dynamic manner. We present current evidence documenting promoter hypermethylation and high levels of gene expression, offer insights into the possible mechanisms by which this occurs, and discuss the potential implications for both research and clinical applications.

Increased WISP1 expression in human osteoarthritic articular cartilage is epigenetically regulated and decreases cartilage matrix production

OBJECTIVES

Previously, we have shown the involvement of Wnt-activated protein Wnt-1-induced signaling protein 1 (WISP1) in the development of OA in mice. Here, we aimed to characterize the relation between WISP1 expression and human OA and its regulatory epigenetic determinants.

METHODS

Preserved and lesioned articular cartilage from end-stage OA patients and non-OA-diagnosed individuals was collected. WISP1 expression was determined using immunohistochemistry and damage was classified using Mankin scoring. RNA expression and DNA methylation were assessed in silico from genome-wide datasets (microarray analysis and RNA sequencing, and 450 k-methylationarrays, respectively). Effects of WISP1 were tested in pellet cultures of primary human chondrocytes.

RESULTS

WISP1 expression in cartilage of OA patients was increased compared with non-OA-diagnosed controls and, within OA patients, WISP1 was even higher in lesioned compared with preserved regions, with expression strongly correlating with Mankin score. In early symptomatic OA patients with disease progression, higher synovial WISP1 expression was observed as compared with non-progressors. Notably, increased WISP1 expression was inversely correlated with methylation levels of a positional CpG-dinucleotide (cg10191240), with lesioned areas showing strong hypomethylation for this CpG as compared with preserved cartilage. Additionally, we observed that methylation levels were allele-dependent for an intronic single-nucleotide polymorphism nearby cg10191240. Finally, addition of recombinant WISP1 to pellets of primary chondrocytes strongly inhibited deposition of extracellular matrix as reflected by decreased pellet circumference, proteoglycan content and decreased expression of matrix components.

CONCLUSION

Increased WISP1 expression is found in lesioned human articular cartilage, and appears epigenetically regulated via DNA methylation. In vitro assays suggest that increased WISP1 is detrimental for cartilage integrity.

Novel Quinoline Compounds Active in Cancer Cells through Coupled DNA Methyltransferase Inhibition and Degradation

DNA methyltransferases (DNMTs) play a relevant role in epigenetic control of cancer cell survival and proliferation. Since only two DNMT inhibitors (azacitidine and decitabine) have been approved to date for the treatment of hematological malignancies, the development of novel potent and specific inhibitors is urgent. Here we describe the design, synthesis, and biological evaluation of a new series of compounds acting at the same time as DNMTs (mainly DNMT3A) inhibitors and degraders. Tested against leukemic and solid cancer cell lines, 2a-c and 4a-c (the last only for leukemias) displayed up to submicromolar antiproliferative activities. In HCT116 cells, such compounds induced EGFP gene expression in a promoter demethylation assay, confirming their demethylating activity in cells. In the same cell line, 2b and 4c chosen as representative samples induced DNMT1 and -3A protein degradation, suggesting for these compounds a double mechanism of DNMT3A inhibition and DNMT protein degradation.

DNA hypomethylation during MSC chondrogenesis occurs predominantly at enhancer regions

Regulation of transcription occurs in a cell type specific manner orchestrated by epigenetic mechanisms including DNA methylation. Methylation changes may also play a key role in lineage specification during stem cell differentiation. To further our understanding of epigenetic regulation in chondrocytes we characterised the DNA methylation changes during chondrogenesis of mesenchymal stem cells (MSCs) by Infinium 450 K methylation array. Significant DNA hypomethylation was identified during chondrogenic differentiation including changes at many key cartilage gene loci. Integration with chondrogenesis gene expression data revealed an enrichment of significant CpGs in upregulated genes, while characterisation of significant CpG loci indicated their predominant localisation to enhancer regions. Comparison with methylation profiles of other tissues, including healthy and diseased adult cartilage, identified chondrocyte-specific regions of hypomethylation and the overlap with differentially methylated CpGs in osteoarthritis. Taken together we have associated DNA methylation levels with the chondrocyte phenotype. The consequences of which has potential to improve cartilage generation for tissue engineering purposes and also to provide context for observed methylation changes in cartilage diseases such as osteoarthritis.