Human epigenome data reveal increased CpG methylation in alternatively spliced sites and putative exonic splicing enhancers

Targeting aryl hydrocarbon receptor to prevent cancer in barrier organs

The skin, lung, and gut are important barrier organs that control how the body reacts to environmental stressors such as ultraviolet (UV) radiation, air pollutants, dietary components, and microorganisms. The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that plays an important role in maintaining homeostasis of barrier organs. AhR was initially discovered as a receptor for environmental chemical carcinogens such as polycyclic aromatic hydrocarbons (PAHs). Activation of AhR pathways by PAHs leads to increased DNA damage and mutations which ultimately lead to carcinogenesis. Ongoing evidence reveals an ever-expanding role of AhR. Recently, AhR has been linked to immune systems by the interaction with the development of natural killer (NK) cells, regulatory T (Treg) cells, and T helper 17 (Th17) cells, as well as the production of immunosuppressive cytokines. However, the role of AhR in carcinogenesis is not as straightforward as we initially thought. Although AhR activation has been shown to promote carcinogenesis in some studies, others suggest that it may act as a tumor suppressor. In this review, we aim to explore the role of AhR in the development of cancer that originates from barrier organs. We also examined the preclinical efficacy data of AhR agonists and antagonists on carcinogenesis to determine whether AhR modulation can be a viable option for cancer chemoprevention.

Epigenetic modifier alpha-ketoglutarate modulates aberrant gene body methylation and hydroxymethylation marks in diabetic heart

BACKGROUND

Diabetic cardiomyopathy (DCM) is a leading cause of death in diabetic patients. Hyperglycemic myocardial microenvironment significantly alters chromatin architecture and the transcriptome, resulting in aberrant activation of signaling pathways in a diabetic heart. Epigenetic marks play vital roles in transcriptional reprogramming during the development of DCM. The current study is aimed to profile genome-wide DNA (hydroxy)methylation patterns in the hearts of control and streptozotocin (STZ)-induced diabetic rats and decipher the effect of modulation of DNA methylation by alpha-ketoglutarate (AKG), a TET enzyme cofactor, on the progression of DCM.

METHODS

Diabetes was induced in male adult Wistar rats with an intraperitoneal injection of STZ. Diabetic and vehicle control animals were randomly divided into groups with/without AKG treatment. Cardiac function was monitored by performing cardiac catheterization. Global methylation (5mC) and hydroxymethylation (5hmC) patterns were mapped in the Left ventricular tissue of control and diabetic rats with the help of an enrichment-based (h)MEDIP-sequencing technique by using antibodies specific for 5mC and 5hmC. Sequencing data were validated by performing (h)MEDIP-qPCR analysis at the gene-specific level, and gene expression was analyzed by qPCR. The mRNA and protein expression of enzymes involved in the DNA methylation and demethylation cycle were analyzed by qPCR and western blotting. Global 5mC and 5hmC levels were also assessed in high glucose-treated DNMT3B knockdown H9c2 cells.

RESULTS

We found the increased expression of DNMT3B, MBD2, and MeCP2 with a concomitant accumulation of 5mC and 5hmC, specifically in gene body regions of diabetic rat hearts compared to the control. Calcium signaling was the most significantly affected pathway by cytosine modifications in the diabetic heart. Additionally, hypermethylated gene body regions were associated with Rap1, apelin, and phosphatidyl inositol signaling, while metabolic pathways were most affected by hyperhydroxymethylation. AKG supplementation in diabetic rats reversed aberrant methylation patterns and restored cardiac function. Hyperglycemia also increased 5mC and 5hmC levels in H9c2 cells, which was normalized by DNMT3B knockdown or AKG supplementation.

CONCLUSION

This study demonstrates that reverting hyperglycemic damage to cardiac tissue might be possible by erasing adverse epigenetic signatures by supplementing epigenetic modulators such as AKG along with an existing antidiabetic treatment regimen.

Fine-mapping and replication of EWAS loci harboring putative epigenetic alterations associated with AD neuropathology in a large collection of human brain tissue samples

INTRODUCTION

Our previous epigenome-wide association study (EWAS) of Alzheimer's disease (AD) in human brain identified 71 CpGs associated with AD pathology. However, due to low coverage of the Illumina platform, many important CpGs might have been missed.

METHODS

In a large collection of human brain tissue samples (N = 864), we fine-mapped previous EWAS loci by targeted bisulfite sequencing and examined their associations with AD neuropathology. DNA methylation was also linked to gene expression of the same brain cortex.

RESULTS

Our targeted sequencing captured 130 CpGs (∼1.2 kb), 93 of which are novel. Of the 130 CpGs, 57 sites (only 17 included in previous EWAS) and 12 gene regions (e.g., ANK1, BIN1, RHBDF2, SPG7, PODXL) were significantly associated with amyloid load. DNA methylation in some regions was associated with expression of nearby genes.

DISCUSSION

Targeted methylation sequencing can validate previous EWAS loci and discover novel CpGs associated with AD pathology.

Limb blastema formation: How much do we know at a genetic and epigenetic level?

Regeneration of missing body parts is an incredible ability which is present in a wide number of species. However, this regenerative capability varies among different organisms. Urodeles (salamanders) are able to completely regenerate limbs after amputation through the essential process of blastema formation. The blastema is a collection of relatively undifferentiated progenitor cells that proliferate and repattern to form the internal tissues of a regenerated limb. Understanding blastema formation in salamanders may enable comparative studies with other animals, including mammals, with more limited regenerative abilities and may inspire future therapeutic approaches in humans. This review focuses on the current state of knowledge about how limb blastemas form in salamanders, highlighting both the possible roles of epigenetic controls in this process as well as limitations to scientific understanding that present opportunities for research.

DNA Methylation Signature of Aging: Potential Impact on the Pathogenesis of Parkinson's Disease

Regulation of gene expression by epigenetic modifications means lasting and heritable changes in the function of genes without alterations in the DNA sequence. Of all epigenetic mechanisms identified thus far, DNA methylation has been of particular interest in both aging and age-related disease research over the last decade given the consistency of site-specific DNA methylation changes during aging that can predict future health and lifespan. An increasing line of evidence has implied the dynamic nature of DNA (de)methylation events that occur throughout the lifespan has a role in the pathophysiology of aging and age-associated neurodegenerative conditions, including Parkinson's disease (PD). In this regard, PD methylome shows, to some extent, similar genome-wide changes observed in the methylome of healthy individuals of matching age. In this review, we start by providing a brief overview of studies outlining global patterns of DNA methylation, then its mechanisms and regulation, within the context of aging and PD. Considering diverging lines of evidence from different experimental and animal models of neurodegeneration and how they combine to shape our current understanding of tissue-specific changes in DNA methylome in health and disease, we report a high-level comparison of the genomic methylation landscapes of brain, with an emphasis on dopaminergic neurons in PD and in natural aging. We believe this will be particularly useful for systematically dissecting overlapping genome-wide alterations in DNA methylation during PD and healthy aging, and for improving our knowledge of PD-specific changes in methylation patterns independent of aging process.

Epigenetic adaptations in drug-tolerant tumor cells

Traditional chemotherapy against cancer is often severely hampered by acquired resistance to the drug. Epigenetic alterations and other mechanisms like drug efflux, drug metabolism, and engagement of survival pathways are crucial in evading drug pressure. Herein, growing evidence suggests that a subpopulation of tumor cells can often tolerate drug onslaught by entering a "persister" state with minimal proliferation. The molecular features of these persister cells are gradually unraveling. Notably, the "persisters" act as a cache of cells that can eventually re-populate the tumor post-withdrawal drug pressure and contribute to acquiring stable drug-resistant features. This underlines the clinical significance of the tolerant cells. Accumulating evidence highlights the importance of modulation of the epigenome as a critical adaptive strategy for evading drug pressure. Chromatin remodeling, altered DNA methylation, and de-regulation of non-coding RNA expression and function contribute significantly to this persister state. No wonder targeting adaptive epigenetic modifications is increasingly recognized as an appropriate therapeutic strategy to sensitize them and restore drug sensitivity. Furthermore, manipulating the tumor microenvironment and "drug holiday" is also explored to maneuver the epigenome. However, heterogeneity in adaptive strategies and lack of targeted therapies have significantly hindered the translation of epigenetic therapy to the clinics. In this review, we comprehensively analyze the epigenetic alterations adapted by the drug-tolerant cells, the therapeutic strategies employed to date, and their limitations and future prospects.

Comprehensive analysis of pre-mRNA alternative splicing regulated by m6A methylation in pig oxidative and glycolytic skeletal muscles

BACKGROUND

Different types of skeletal myofibers exhibit distinct physiological and metabolic properties that are associated with meat quality traits in livestock. Alternative splicing (AS) of pre-mRNA can generate multiple transcripts from an individual gene by differential selection of splice sites. N6-methyladenosine (m6A) is the most abundant modification in mRNAs, but its regulation for AS in different muscles remains unknown.  RESULTS: We characterized AS events and m6A methylation pattern in pig oxidative and glycolytic muscles. A tota1 of 1294 differential AS events were identified, and differentially spliced genes were significantly enriched in processes related to different phenotypes between oxidative and glycolytic muscles. We constructed the regulatory network between splicing factors and corresponding differential AS events and identified NOVA1 and KHDRBS2 as key splicing factors. AS event was enriched in m6A-modified genes, and the methylation level was positively correlated with the number of AS events in genes. The dynamic change in m6A enrichment was associated with 115 differentially skipping exon (SE-DAS) events within 92 genes involving in various processes, including muscle contraction and myofibril assembly. We obtained 23.4% SE-DAS events (27/115) regulated by METTL3-meditaed m6A and experimentally validated the aberrant splicing of ZNF280D, PHE4DIP, and NEB. The inhibition of m6A methyltransferase METTL3 could induce the conversion of oxidative fiber to glycolytic fiber in PSCs.

CONCLUSION

Our study suggested that m6A modification could contribute to significant difference in phenotypes between oxidative and glycolytic muscles by mediating the regulation of AS. These findings would provide novel insights into mechanisms underlying muscle fiber conversion.

Aryl hydrocarbon receptor as a DNA methylation reader in the stress response pathway

The aryl hydrocarbon receptor (AhR) mediates various cellular responses upon exposure to exogenous and endogenous stress factors. In these responses, AhR plays a dual role as a stress sensor for detecting various AhR ligands and as a transcription factor that upregulates the expression of downstream effector genes, such as those encoding drug-metabolizing enzymes. As a transcription factor, it selectively binds to the unmethylated form of a specific sequence called the xenobiotic responsive element (XRE). We suggest that AhR is a novel DNA methylation reader, unlike classical methylation readers, such as methyl-CpG-binding protein 2, which binds to methylated sequences. Under physiological conditions of continuous exposure to endogenous AhR ligands, such as kynurenine, methylation states of the individual target XREs must be strictly regulated to select and coordinate the expression of downstream genes responsible for maintaining homeostasis in the body. In contrast, long-term exposure to AhR ligands frequently leads to changes in the methylation patterns around the XRE sequence. These data indicate that AhR may contribute to the adaptive cellular response to various stresses by modulating DNA methylation. Thus, the DNA methylation profile of AhR target genes should be dynamically controlled through a balance between robustness and flexibility under both physiological and stress conditions. AhR is a pivotal player in the regulation of stress response as it shows versatility by functioning as a stress sensor, methylation reader, and putative methylation modulator.

Histone deacetylase in neuropathology

Neuroepigenetics, a new branch of epigenetics, plays an important role in the regulation of gene expression. Neuroepigenetics is associated with holistic neuronal function and helps in formation and maintenance of memory and learning processes. This includes neurodevelopment and neurodegenerative defects in which histone modification enzymes appear to play a crucial role. These modifications, carried out by acetyltransferases and deacetylases, regulate biologic and cellular processes such as apoptosis and autophagy, inflammatory response, mitochondrial dysfunction, cell-cycle progression and oxidative stress. Alterations in acetylation status of histone as well as non-histone substrates lead to transcriptional deregulation. Histone deacetylase decreases acetylation status and causes transcriptional repression of regulatory genes involved in neural plasticity, synaptogenesis, synaptic and neural plasticity, cognition and memory, and neural differentiation. Transcriptional deactivation in the brain results in development of neurodevelopmental and neurodegenerative disorders. Mounting evidence implicates histone deacetylase inhibitors as potential therapeutic targets to combat neurologic disorders. Recent studies have targeted naturally-occurring biomolecules and micro-RNAs to improve cognitive defects and memory. Multi-target drug ligands targeting HDAC have been developed and used in cell-culture and animal-models of neurologic disorders to ameliorate synaptic and cognitive dysfunction. Herein, we focus on the implications of histone deacetylase enzymes in neuropathology, their regulation of brain function and plausible involvement in the pathogenesis of neurologic defects.

The methylation landscape and its role in domestication and gene regulation in the chicken

Domestication is one of the strongest examples of artificial selection and has produced some of the most extreme within-species phenotypic variation known. In the case of the chicken, it has been hypothesized that DNA methylation may play a mechanistic role in the domestication response. By inter-crossing wild-derived red junglefowl with domestic chickens, we mapped quantitative trait loci for hypothalamic methylation (methQTL), gene expression (eQTL) and behaviour. We find large, stable methylation differences, with 6,179 cis and 2,973 trans methQTL identified. Over 46% of the trans effects were genotypically controlled by five loci, mainly associated with increased methylation in the junglefowl genotype. In a third of eQTL, we find that there is a correlation between gene expression and methylation, while statistical causality analysis reveals multiple instances where methylation is driving gene expression, as well as the reverse. We also show that methylation is correlated with some aspects of behavioural variation in the inter-cross. In conclusion, our data suggest a role for methylation in the regulation of gene expression underlying the domesticated phenotype of the chicken.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infection: Triggering a Lethal Fight to Keep Control of the Ten-Eleven Translocase (TET)-Associated DNA Demethylation?

The extended and diverse interference of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in multiple host functions and the diverse associated symptoms implicate its involvement in fundamental cellular regulatory processes. The activity of ten-eleven translocase 2 (TET2) responsible for selective DNA demethylation, has been recently identified as a regulator of endogenous virus inactivation and viral invasion, possibly by proteasomal deregulation of the TET2/TET3 activities. In a recent report, we presented a detailed list of factors that can be affected by TET activity, including recognition of zinc finger protein binding sites and bimodal promoters, by enhancing the flexibility of adjacent sequences. In this review, we summarize the TET-associated processes and factors that could account for SARS-CoV-2 diverse symptoms. Moreover, we provide a correlation for the observed virus-induced symptoms that have been previously associated with TET activities by in vitro and in vitro studies. These include early hypoxia, neuronal regulation, smell and taste development, liver, intestinal, and cardiomyocyte differentiation. Finally, we propose that the high mortality of SARS-CoV-2 among adult patients, the different clinical symptoms of adults compared to children, the higher risk of patients with metabolic deregulation, and the low mortality rates among women can all be accounted for by the complex balance of the three enzymes with TET activity, which is developmentally regulated. This activity is age-dependent, related to telomere homeostasis and integrity, and associated with X chromosome inactivation via (de)regulation of the responsible XIST gene expression.

In silico structural analysis of sequences containing 5-hydroxymethylcytosine reveals its potential as binding regulator for development, ageing and cancer-related transcription factors

The presence of 5-hydroxymethyl cytosine in DNA has been previously associated with ageing. Using in silico analysis of normal liver samples we presently observed that in 5-hydroxymethyl cytosine sequences, DNA methylation is dependent on the co-presence of G-quadruplexes and palindromes. This association exhibits discrete patterns depending on G-quadruplex and palindrome densities. DNase-Seq data show that 5-hydroxymethyl cytosine sequences are common among liver nucleosomes (p < 2.2x10^-16) and threefold more frequent than nucleosome sequences. Nucleosomes lacking palindromes and potential G-quadruplexes are rare in vivo (1%) and nucleosome occupancy potential decreases with increasing G-quadruplexes. Palindrome distribution is similar to that previously reported in nucleosomes. In low and mixed complexity sequences 5-hydroxymethyl cytosine is frequently located next to three elements: G-quadruplexes or imperfect G-quadruplexes with CpGs, or unstable hairpin loops (TCCCAY₆TGGGA) mostly located in antisense strands or finally A-/T-rich segments near these motifs. The high frequencies and selective distribution of pentamer sequences (including TCCCA, TGGGA) probably indicate the positive contribution of 5-hydroxymethyl cytosine to stabilize the formation of structures unstable in the absence of this cytosine modification. Common motifs identified in all total 5-hydroxymethyl cytosine-containing sequences exhibit high homology to recognition sites of several transcription factor families: homeobox, factors involved in growth, mortality/ageing, cancer, neuronal function, vision, and reproduction. We conclude that cytosine hydroxymethylation could play a role in the recognition of sequences with G-quadruplexes/palindromes by forming epigenetically regulated DNA 'springs' and governing expansions or compressions recognized by different transcription factors or stabilizing nucleosomes. The balance of these epigenetic elements is lost in hepatocellular carcinoma.

Genome-wide analysis reveals the association between alternative splicing and DNA methylation across human solid tumors

Background

Dysregulation of alternative splicing (AS) is a critical signature of cancer. However, the regulatory mechanisms of cancer-specific AS events, especially the impact of DNA methylation, are poorly understood.

Methods

By using The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA data for ten solid tumor types, association analysis was performed to characterize the potential link between cancer-specific AS and DNA methylation. Functional and pathway enrichment analyses were performed, and the protein-protein interaction (PPI) network was constructed with the String website. The prognostic analysis was carried out with multivariate Cox regressions models.

Results

15,818 AS events in 3955 annotated genes were identified across ten solid tumor types. The different DNA methylation patterns between tumor and normal tissues at the corresponding alternative spliced exon boundaries were shown, and 51.3% of CpG sites (CpGs) revealed hypomethylated in tumors. Notably, 607 CpGs were found to be highly correlated with 369 cancer-specific AS events after permutation tests. Among them, the hypomethylated CpGs account for 52.7%, and the number of down-regulated exons was 173. Furthermore, we found 38 AS events in 35 genes could serve as new molecular biomarkers to predict patient survival.

Conclusions

Our study described the relationship between DNA methylation and AS events across ten human solid tumor types and provided new insights into intragenic DNA methylation and exon usage during the AS process.

Association of intronic DNA methylation and hydroxymethylation alterations in the epigenetic etiology of dilated cardiomyopathy

In this study, we investigated the role of DNA methylation [5-methylcytosine (5mC)] and 5-hydroxymethylcytosine (5hmC), epigenetic modifications that regulate gene activity, in dilated cardiomyopathy (DCM). A MYBPC3 mutant mouse model of DCM was compared with wild type and used to profile genomic 5mC and 5hmC changes by Chip-seq, and gene expression levels were analyzed by RNA-seq. Both 5mC-altered genes (957) and 5hmC-altered genes (2,022) were identified in DCM hearts. Diverse gene ontology and KEGG pathways were enriched for DCM phenotypes, such as inflammation, tissue fibrosis, cell death, cardiac remodeling, cardiomyocyte growth, and differentiation, as well as sarcomere structure. Hierarchical clustering of mapped genes affected by 5mC and 5hmC clearly differentiated DCM from wild-type phenotype. Based on these data, we propose that genomewide 5mC and 5hmC contents may play a major role in DCM pathogenesis. NEW & NOTEWORTHY Our data demonstrate that development of dilated cardiomyopathy in mice is associated with significant epigenetic changes, specifically in intronic regions, which, when combined with gene expression profiling data, highlight key signaling pathways involved in pathological cardiac remodeling and heart contractile dysfunction.

Tissue-specific epigenetics of atherosclerosis-related ANGPT and ANGPTL genes

Aim:

To understand tissue-specific regulation of angiopoietin/angiopoietin-like (ANGPT/ANGPTL) genes (especially the five genes embedded in introns of host genes) and their association with atherosclerosis.

Methods:

Transcription and epigenomic databases from various normal tissues were examined in the vicinity of ANGPT1, ANGPT2, ANGPTL1, ANGPTL2, ANGPTL3, ANGPTL4 and ANGPTL8.

Results:

We identified tissue-specific enhancer chromatin regions that are likely to regulate transcription of ANGPT/ANGPTL genes and were intragenic, intergenic or host gene-linked. In addition, we found atherosclerosis-linked differentially methylated regions associated with ANGPT2 and with sequences encoding miR-145, a microRNA that targets ANGPT2 mRNA in cancers.

Conclusion:

Our findings implicate enhancers as major contributors to tissue-specific expression of ANGPT/ANGPTL genes, which play critical roles in angiogenesis, atherosclerosis, cancer, and inflammatory and metabolic diseases.

DNMT3B Functions: Novel Insights From Human Disease

DNA methylation plays important roles in gene expression regulation and chromatin structure. Its proper establishment and maintenance are essential for mammalian development and cellular differentiation. DNMT3B is the major de novo DNA methyltransferase expressed and active during the early stage of embryonic development, including implantation. In addition to its well-known role to methylate centromeric, pericentromeric, and subtelomeric repeats, recent observations suggest that DNMT3B acts as the main enzyme methylating intragenic regions of active genes. Although largely studied, much remains unknown regarding how these specific patterns of de novo CpG methylation are established in mammalian cells, and which are the rules governing DNMT3B recruitment and activity. Latest evidence indicates that DNMT3B recruitment is regulated by numerous mechanisms including chromatin modifications, transcription levels, non-coding RNAs, and the presence of DNA-binding factors. DNA methylation abnormalities are a common mark of human diseases involving chromosomal and genomic instabilities, such as inherited disease and cancer. The autosomal recessive Immunodeficiency, Centromeric instability and Facial anomalies syndrome, type I (ICF-1), is associated to hypomorphic mutations in DNMT3B gene, while its altered expression has been correlated with the development of tumors. In both cases, this implies that abnormal DNA hypomethylation and hypermethylation patterns affect gene expression and genomic architecture contributing to the pathological states. We will provide an overview of the most recent research aimed at deciphering the molecular mechanisms by which DNMT3B abnormalities are associated with the onset and progression of these pathologies.

Age-dependent methylation in epigenetic clock CpGs is associated with G-quadruplex, co-transcriptionally formed RNA structures and tentative splice sites

Horvath's epigenetic clock consists of 353 CpGs whose methylation levels can accurately predict the age of individuals. Using bioinformatics analysis, we investigated the conformation, energy characteristics and presence of tentative splice sites of the sequences surrounding the epigenetic clock CpGs, in relation to the median methylation changes in different ages, the presence of CpG islands and their position in genes. Common characteristics in the 100 nt sequences surrounding the epigenetic clock CpGs are G-quadruplexes and/or tentative splice site motifs. Median methylation increases significantly in sequences which adopt less stable structures during transcription. Methylation is higher when CpGs overlap with G-quadruplexes than when they precede them. Median methylation in epigenetic clock CpGs is higher in sequences expressed as single products rather than in multiple products and those containing single donors and multiple acceptors. Age-related methylation variation is significant in sequences without G-quadruplexes, particularly those producing low stability nascent RNA and those with splice sites. CpGs in sequences close to transcription start sites and those which are possibly never expressed (hypothetical proteins) undergo similar extent of age-related median methylation decrease and increase. Preservation of methylation is observed in CpG islands without G-quadruplexes, contrary to CpGs far from CpG islands (open sea). Sequences containing G-quadruplexes and RNA pseudoknots, determining the recognition by H3K27 histone methyltransferase, are hypomethylated. The presented structural DNA and co-transcriptional RNA analysis of epigenetic clock sequences, foreshadows the association of age-related methylation changes with the principle biological processes of DNA and histone methylation, splicing and chromatin silencing.

A survey of transcriptome complexity in Sus scrofa using single-molecule long-read sequencing

Alternative splicing (AS) and fusion transcripts produce a vast expansion of transcriptomes and proteomes diversity. However, the reliability of these events and the extend of epigenetic mechanisms have not been adequately addressed due to its limitation of uncertainties about the complete structure of mRNA. Here we combined single-molecule real-time sequencing, Illumina RNA-seq and DNA methylation data to characterize the landscapes of DNA methylation on AS, fusion isoforms formation and lncRNA feature and further to unveil the transcriptome complexity of pig. Our analysis identified an unprecedented scale of high-quality full-length isoforms with over 28,127 novel isoforms from 26,881 novel genes. More than 92,000 novel AS events were detected and intron retention predominated in AS model, followed by exon skipping. Interestingly, we found that DNA methylation played an important role in generating various AS isoforms by regulating splicing sites, promoter regions and first exons. Furthermore, we identified a large of fusion transcripts and novel lncRNAs, and found that DNA methylation of the promoter and gene body could regulate lncRNA expression. Our results significantly improved existed gene models of pig and unveiled that pig AS and epigenetic modify were more complex than previously thought.

Do social insects support Haig's kin theory for the evolution of genomic imprinting?

Although numerous imprinted genes have been described in several lineages, the phenomenon of genomic imprinting presents a peculiar evolutionary problem. Several hypotheses have been proposed to explain gene imprinting, the most supported being Haig's kinship theory. This theory explains the observed pattern of imprinting and the resulting phenotypes as a competition for resources between related individuals, but despite its relevance it has not been independently tested. Haig's theory predicts that gene imprinting should be present in eusocial insects in many social scenarios. These lineages are therefore ideal for testing both the theory's predictions and the mechanism of gene imprinting. Here we review the behavioral evidence of genomic imprinting in eusocial insects, the evidence of a mechanism for genomic imprinting and finally we evaluate recent results showing parent of origin allele specific expression in honeybees in the light of Haig's theory.

Neuroepigenetics of Aging and Age-Related Neurodegenerative Disorders

Neurodegenerative diseases are complex, progressive disorders and affect millions of people worldwide, contributing significantly to the global burden of disease. In recent years, research has begun to investigate epigenetic mechanisms for a potential role in disease etiology. In this chapter, we describe the current state of play for epigenetic research into neurodegenerative disorders including Alzheimer's disease, Parkinson's disease and Huntington's disease. We focus on the recent evidence for a potential role of DNA modifications, histone modifications and non-coding RNA in the etiology of these disorders. Finally, we discuss how new technological and bioinformatics advances in the field of epigenetics could further progress our understanding about the underlying mechanisms of neurodegenerative diseases.

TFAP2E Methylation and Expression Status Does Not Predict Response to 5-FU-based Chemotherapy in Colorectal Cancer

Purpose: A recent study reported that 5-fluorouracil (5-FU)-based chemotherapy is less effective in treating patients with advanced colorectal cancer demonstrating hypermethylation of the TFAP2E gene. The aim of our study was to confirm and validate these findings in large, uniformly treated, well-characterized patient cohorts.Experimental Design: Two cohorts of 783 patients with colorectal cancer: 532 from a population-based, multicenter cohort (EPICOLON I) and 251 patients from a clinic-based trial were used to study the effectiveness of TFAP2E methylation and expression as a predictor of response of colorectal cancer patients to 5-FU-based chemotherapy. DNA methylation status of the TFAP2E gene in patients with colorectal cancer was assessed by quantitative bisulfite pyrosequencing analysis. IHC analysis of the TFAP2E protein expression was also performed.Results: Correlation between TFAP2E methylation status and IHC staining was performed in 607 colorectal cancer samples. Among 357 hypermethylated tumors, only 141 (39.6%) exhibited loss of protein expression. Survival was not affected by TFAP2E hypermethylation in stage IV patients [HR, 1.21; 95% confidence interval (CI), 0.79-1.87; log-rank P = 0.6]. In stage II-III cases, disease-free survival was not influenced by TFAP2E hypermethylation status in 5-FU-treated (HR, 0.91; 95% CI, 0.52-1.59; log-rank P = 0.9) as well as in nontreated patients (HR, 0.88; 95% CI, 0.5-1.54; log-rank P = 0.7).Conclusions:TFAP2E hypermethylation does not correlate with loss of its protein expression. Our large, systematic, and comprehensive study indicates that TFAP2E methylation and expression may not play a major role in predicting response to 5-FU-based chemotherapy in patients with colorectal cancer. Clin Cancer Res; 24(12); 2820-7. ©2018 AACR.

Population-dependent Intron Retention and DNA Methylation in Breast Cancer

Regulation of gene expression by DNA methylation in gene promoter regions is well studied; however, the effects of methylation in the gene body (exons and introns) on gene expression are comparatively understudied. Recently, hypermethylation has been implicated in the inclusion of alternatively spliced exons; moreover, exon recognition can be enhanced by recruiting the methyl-CpG-binding protein (MeCP2) to hypermethylated sites. This study examines whether the methylation status of an intron is correlated with how frequently the intron is retained during splicing using DNA methylation and RNA sequencing data from breast cancer tissue specimens in The Cancer Genome Atlas. Interestingly, hypomethylation of introns is correlated with higher levels of intron expression in mRNA and the methylation level of an intron is inversely correlated with its retention in mRNA from the gene in which it is located. Furthermore, significant population differences were observed in the methylation level of retained introns. In African-American donors, retained introns were not only less methylated compared to European-American donors, but also were more highly expressed. This underscores the need for understanding epigenetic differences in populations and their correlation with breast cancer is an important step toward achieving personalized cancer care.Implications: This research contributes to the understanding of how epigenetic markers in the gene body communicate with the transcriptional machinery to control transcript diversity and differential biological response to changes in methylation status could underlie some of the known, yet unexplained, disparities in certain breast cancer patient populations. Mol Cancer Res; 16(3); 461-9. ©2018 AACR.

ICF-specific DNMT3B dysfunction interferes with intragenic regulation of mRNA transcription and alternative splicing

Hypomorphic mutations in DNA-methyltransferase DNMT3B cause majority of the rare disorder Immunodeficiency, Centromere instability and Facial anomalies syndrome cases (ICF1). By unspecified mechanisms, mutant-DNMT3B interferes with lymphoid-specific pathways resulting in immune response defects. Interestingly, recent findings report that DNMT3B shapes intragenic CpG-methylation of highly-transcribed genes. However, how the DNMT3B-dependent epigenetic network modulates transcription and whether ICF1-specific mutations impair this process remains unknown. We performed a transcriptomic and epigenomic study in patient-derived B-cell lines to investigate the genome-scale effects of DNMT3B dysfunction. We highlighted that altered intragenic CpG-methylation impairs multiple aspects of transcriptional regulation, like alternative TSS usage, antisense transcription and exon splicing. These defects preferentially associate with changes of intragenic H3K4me3 and at lesser extent of H3K27me3 and H3K36me3. In addition, we highlighted a novel DNMT3B activity in modulating the self-regulatory circuit of sense-antisense pairs and the exon skipping during alternative splicing, through interacting with RNA molecules. Strikingly, altered transcription affects disease relevant genes, as for instance the memory-B cell marker CD27 and PTPRC genes, providing us with biological insights into the ICF1-syndrome pathogenesis. Our genome-scale approach sheds light on the mechanisms still poorly understood of the intragenic function of DNMT3B and DNA methylation in gene expression regulation.

A saga of cancer epigenetics: linking epigenetics to alternative splicing

The discovery of an increasing number of alternative splicing events in the human genome highlighted that ∼94% of genes generate alternatively spliced transcripts that may produce different protein isoforms with diverse functions. It is now well known that several diseases are a direct and indirect consequence of aberrant splicing events in humans. In addition to the conventional mode of alternative splicing regulation by 'cis' RNA-binding sites and 'trans' RNA-binding proteins, recent literature provides enormous evidence for epigenetic regulation of alternative splicing. The epigenetic modifications may regulate alternative splicing by either influencing the transcription elongation rate of RNA polymerase II or by recruiting a specific splicing regulator via different chromatin adaptors. The epigenetic alterations and aberrant alternative splicing are known to be associated with various diseases individually, but this review discusses/highlights the latest literature on the role of epigenetic alterations in the regulation of alternative splicing and thereby cancer progression. This review also points out the need for further studies to understand the interplay between epigenetic modifications and aberrant alternative splicing in cancer progression.

High-throughput sequencing offers new insights into 5-hydroxymethylcytosine

Chemical modifications of DNA comprise epigenetic mechanisms that contribute to the maintenance of cellular activities and memory. Although the function of 5-methylcytosine (5-mC) has been extensively studied, little is known about the function(s) of relatively rarer and underappreciated cytosine modifications including 5-hydroxymethylcytosine (5-hmC). The discovery that ten-eleven translocation (Tet) proteins mediate conversion of 5-mC to 5-hmC, and other oxidation derivatives, sparked renewed interest to understand the biological role of 5-hmC. Studies examining total 5-hmC levels revealed the highly dynamic yet tissue-specific nature of this modification, implicating a role in epigenetic regulation and development. Intriguingly, 5-hmC levels are highest during early development and in the brain where abnormal patterns of 5-hmC have been observed in disease conditions. Thus, 5-hmC adds to the growing list of epigenetic modifications with potential utility in clinical applications and warrants further investigation. This review discusses the emerging functional roles of 5-hmC in normal and disease states, focusing primarily on insights provided by recent studies exploring the genome-wide distribution of this modification in mammals.

Association of two synonymous splicing-associated CpG single nucleotide polymorphisms in calpain 10 and solute carrier family 2 member 2 with type 2 diabetes

Coding synonymous single nucleotide polymorphisms (SNPs) have attracted little attention until recently. However, such SNPs located in epigenetic, CpG sites modifying exonic splicing enhancers (ESEs) can be informative with regards to the recently verified association of intragenic methylation and splicing. The present study describes the association of type 2 diabetes (T2D) with the exonic, synonymous, epigenetic SNPs, rs3749166 in calpain 10 (CAPN10) glucose transporter (GLUT4) translocator and rs5404 in solute carrier family 2, member 2 (SLC2A2), also termed GLUT2, which, according to prior bioinformatic analysis, strongly modify the splicing potential of glucose transport-associated genes. Previous association studies reveal that only rs5404 exhibits a strong negative T2D association, while data on the CAPN10 polymorphism are contradictory. In the present study DNA from blood samples of 99 Greek non-diabetic control subjects and 71 T2D patients was analyzed. In addition, relevant publicly available cases (40) resulting from examination of 110 Personal Genome Project data files were analyzed. The frequency of the rs3749166 A allele, was similar in the patients and non-diabetic control subjects. However, AG heterozygotes were more frequent among patients (73.24% for Greek patients and 54.55% for corresponding non-diabetic control subjects; P=0.0262; total cases, 52.99 and 75.00%, respectively; P=0.0039). The rs5404 T allele was only observed in CT heterozygotes (Greek non-diabetic control subjects, 39.39% and Greek patients, 22.54%; P=0.0205; total cases, 34.69 and 21.28%, respectively; P=0.0258). Notably, only one genotype, heterozygous AG/CC, was T2D-associated (Greek non-diabetic control subjects, 29.29% and Greek patients, 56.33%; P=0.004; total cases, 32.84 and 56.58%, respectively; P=0.0008). Furthermore, AG/CC was strongly associated with very high (≥8.5%) glycosylated plasma hemoglobin levels among patients (P=0.0002 for all cases). These results reveal the complex heterozygotic SNP association with T2D, and indicate possible synergies of these epigenetic, splicing-regulatory, synonymous SNPs, which modify the splicing potential of two alternative glucose transport-associated genes.

The integration of epigenetics and genetics in nutrition research for CVD risk factors

There is increasing evidence documenting gene-by-environment (G × E) interactions for CVD related traits. However, the underlying mechanisms are still unclear. DNA methylation may represent one of such potential mechanisms. The objective of this review paper is to summarise the current evidence supporting the interplay among DNA methylation, genetic variants, and environmental factors, specifically (1) the association between SNP and DNA methylation; (2) the role that DNA methylation plays in G × E interactions. The current evidence supports the notion that genotype-dependent methylation may account, in part, for the mechanisms underlying observed G × E interactions in loci such asAPOE, IL6and ATP-binding cassette A1. However, these findings should be validated using intervention studies with high level of scientific evidence. The ultimate goal is to apply the knowledge and the technology generated by this research towards genetically based strategies for the development of personalised nutrition and medicine.

Integrative analysis of methylome and transcriptome in human blood identifies extensive sex- and immune cell-specific differentially methylated regions

The relationship between DNA methylation and gene expression is complex and elusive. To further elucidate these relations, we performed an integrative analysis of the methylome and transcriptome of 4 circulating immune cell subsets (B cells, monocytes, CD4(+), and CD8(+) T cells) from healthy females. Additionally, in light of the known sex bias in the prevalence of several immune-mediated diseases, the female datasets were compared with similar public available male data sets. Immune cell-specific differentially methylated regions (DMRs) were found to be highly similar between sexes, with an average correlation coefficient of 0.82; however, numerous sex-specific DMRs, shared by the cell subsets, were identified, mainly on autosomal chromosomes. This provides a list of highly interesting candidate genes to be studied in disorders with sexual dimorphism, such as autoimmune diseases. Immune cell-specific DMRs were mainly located in the gene body and intergenic region, distant from CpG islands but overlapping with enhancer elements, indicating that distal regulatory elements are important in immune cell specificity. In contrast, sex-specific DMRs were overrepresented in CpG islands, suggesting that the epigenetic regulatory mechanisms of sex and immune cell specificity may differ. Both positive and, more frequently, negative correlations between subset-specific expression and methylation were observed, and cell-specific DMRs of both interactions were associated with similar biological pathways, while sex-specific DMRs were linked to networks of early development or estrogen receptor and immune-related molecules. Our findings of immune cell- and sex-specific methylome and transcriptome profiles provide novel insight on their complex regulatory interactions and may particularly contribute to research of immune-mediated diseases.

Responses of bovine early embryos to S-adenosyl methionine supplementation in culture

AIM

There is a growing concern about the potential adverse effects of high dose folic acid (FA) supplementation before and during pregnancy. FA metabolism generates S-adenosyl methionine (SAM) which is an important cofactor of epigenetic programming. We sought to assess the impact of a large dose of SAM on early embryo development.

MATERIALS & METHODS

In vitro cultured bovine embryos were treated with SAM from the eight-cell stage to the blastocyst stage. In addition to the phenotype, the genome-wide epigenetic and transcription profiles were analyzed.

RESULTS

Treatment significantly improved embryo hatching and caused a shift in sex ratio in favor of males. SAM caused genome-wide hypermethylation mainly in exonic regions and in CpG islands. Although differentially expressed genes were associated with response to nutrients and developmental processes, no correspondence was found with the differentially methylated regions, suggesting that cellular responses to SAM treatment during early embryo development may not require DNA methylation-driven changes.

CONCLUSION

Since bovine embryos were not indifferent to SAM, effects of large-dose FA supplements on early embryonic development in humans cannot be ruled out.

In silico analysis of regulatory and structural motifs of the ovine HSP90AA1 gene

Gene promoters are essential regions of DNA where the transcriptional molecular machinery to produce RNA molecules is recruited. In this process, DNA epigenetic modifications can acquire a fundamental role in the regulation of gene expression. Recently, in a previous work of our group, functional features and DNA methylation involved in the ovine HSP90AA1 gene expression regulation have been observed. In this work, we report a combination of methylation analysis by bisulfite sequencing in several tissues and at different developmental stages together with in silico bioinformatic analysis of putative regulating factors in order to identify regulative mechanisms both at the promoter and gene body. Our results show a "hybrid structure" (TATA box + CpG island) of the ovine HSP90AA1 gene promoter both in somatic and non-differentiated germ tissues, revealing the ability of the HSP90AA1 gene to be regulated both in an inducible and constitutive fashion. In addition, in silico analysis showed that several putative alternative spliced regulatory motifs, exonic splicing enhancers (ESEs), and G-quadruplex secondary structures were somehow related to the DNA methylation pattern found. The results obtained here could help explain the differences in cell-type transcripts, tissue expression rate, and transcription silencing mechanisms found in this gene.

Analysis of DNA methylation landscape reveals the roles of DNA methylation in the regulation of drug metabolizing enzymes

Background

Drug metabolizing enzymes (DMEs) exhibit dramatic inter- and intra-individual variability in expression and activity. However, the mechanisms determining this variability have not been fully elucidated. The aim of this study was to evaluate the biological significance of DNA methylation in the regulation of DME genes by genome-wide integrative analysis.

Results

DNA methylation and mRNA expression profiles of human tissues and hepatoma cells were examined by microarrays. The data were combined with GEO datasets of liver tissues, and integrative analysis was performed on selected DME genes. Detailed DNA methylation statuses at individual CpG sites were evaluated by DNA methylation mapping. From analysis of 20 liver tissues, highly variable DNA methylation was observed in 37 DME genes, 7 of which showed significant inverse correlations between DNA methylation and mRNA expression. In hepatoma cells, treatment with a demethylating agent resulted in upregulation of 5 DME genes, which could be explained by DNA methylation status. Interestingly, some DMEs were suggested to act as tumor-suppressor or housekeeper based on their unique DNA methylation features. Moreover, tissue-specific and age-dependent expression of UDP-glucuronosyltransferase 1A splicing variants was associated with DNA methylation status of individual first exons.

Conclusions

Some DME genes were regulated by DNA methylation, potentially resulting in inter- and intra-individual differences in drug metabolism. Analysis of DNA methylation landscape facilitated elucidation of the role of DNA methylation in the regulation of DME genes, such as mediator of inter-individual variability, guide for correct alternative splicing, and potential tumor-suppressor or housekeeper.

Electronic supplementary material

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

Expression of mouse Dab2ip transcript variants and gene methylation during brain development

Dab2ip (DOC-2/DAB2 interacting protein) is a RasGAP protein which shows a growth-inhibitory effect in human prostate cancer cell lines. Recent studies have shown that Dab2ip also plays an important role in regulating dendrite development and neuronal migration during brain development. In this study, we provide a more complete description of the mouse Dab2ip (mDab2ip) gene locus and examined DNA methylation and expression of Dab2ip during cerebellar development. Analysis of cDNA sequences in public databases revealed a total of 20 possible exons for mDab2ip gene, spanning over 172kb. Using Cap Analysis of Gene Expression (CAGE) data available through FANTOM5 project, we deduced five different transcription start sites for mDab2ip. Here, we characterized three different mDab2ip transcript variants beginning with exon 1. These transcripts varied by the presence or absence of exons 3 and 5, which encode a putative nuclear localization signal and the N-terminal region of a PH-domain, respectively. The 5' region of the mDab2ip gene contains three putative CpG islands (CpG131, CpG54, and CpG85). Interestingly, CpG54 and CpG85 are localized on exons 3 and 5. Bisulfate DNA sequencing showed that methylation level of CpG54 remained constant whereas methylation of CpG85 increased during cerebellar development. Real-time PCR analysis showed that the proportion of PH-domain containing mDab2ip transcripts increased during cerebellar development, in correlation with the increase in CpG85 methylation. These data suggest that site-specific methylation of mDab2ip gene during cerebellar development may play a role in inclusion of exon 5, resulting in a Dab2ip transcript variant that encodes a full pleckstrin homology (PH) domain.

Incorporating alternative splicing and mRNA editing into the genetic analysis of complex traits

The nomination of candidate genes underlying complex traits is often focused on genetic variations that alter mRNA abundance or result in non-conservative changes in amino acids. Although inconspicuous in complex trait analysis, genetic variants that affect splicing or RNA editing can also generate proteomic diversity and impact genetic traits. Indeed, it is known that splicing and RNA editing modulate several traits in humans and model organisms. Using high-throughput RNA sequencing (RNA-seq) analysis, it is now possible to integrate the genetics of transcript abundance, alternative splicing (AS) and editing with the analysis of complex traits. We recently demonstrated that both AS and mRNA editing are modulated by genetic and environmental factors, and potentially engender phenotypic diversity in a genetically segregating mouse population. Therefore, the analysis of splicing and RNA editing can expand not only the regulatory landscape of transcriptome and proteome complexity, but also the repertoire of candidate genes for complex traits.

Conserved and divergent patterns of DNA methylation in higher vertebrates

DNA methylation in the genome plays a fundamental role in the regulation of gene expression and is widespread in the genome of eukaryotic species. For example, in higher vertebrates, there is a "global" methylation pattern involving complete methylation of CpG sites genome-wide, except in promoter regions that are typically enriched for CpG dinucleotides, or so called "CpG islands." Here, we comprehensively examined and compared the distribution of CpG sites within ten model eukaryotic species and linked the observed patterns to the role of DNA methylation in controlling gene transcription. The analysis revealed two distinct but conserved methylation patterns for gene promoters in human and mouse genomes, involving genes with distinct distributions of promoter CpGs and gene expression patterns. Comparative analysis with four other higher vertebrates revealed that the primary regulatory role of the DNA methylation system is highly conserved in higher vertebrates.

DNA modifications: function and applications in normal and disease States

Epigenetics refers to a variety of processes that have heritable effects on gene expression programs without changes in DNA sequence. Key players in epigenetic control are chemical modifications to DNA, histone, and non-histone chromosomal proteins, which establish a complex regulatory network that controls genome function. Methylation of DNA at the fifth position of cytosine in CpG dinucleotides (5-methylcytosine, 5mC), which is carried out by DNA methyltransferases, is commonly associated with gene silencing. However, high resolution mapping of DNA methylation has revealed that 5mC is enriched in exonic nucleosomes and at intron-exon junctions, suggesting a role of DNA methylation in the relationship between elongation and RNA splicing. Recent studies have increased our knowledge of another modification of DNA, 5-hydroxymethylcytosine (5hmC), which is a product of the ten-eleven translocation (TET) proteins converting 5mC to 5hmC. In this review, we will highlight current studies on the role of 5mC and 5hmC in regulating gene expression (using some aspects of brain development as examples). Further the roles of these modifications in detection of pathological states (type 2 diabetes, Rett syndrome, fetal alcohol spectrum disorders and teratogen exposure) will be discussed.

Occurrence of HHIP gene CpG island methylation in gastric cancer

The present study aimed to observe the methylation status of the CpG islands at the human hedgehog interacting protein (HHIP) gene in gastric cancer tissues, peritumoral tissues and the AGS cell line, to analyze the association between the methylation status of the CpG islands and the tumorigenesis of gastric cancer. The HHIP mRNA expression in 60 human gastric carcinnoma tissues, peritumoral tissues and the gastric carcinoma AGS cell line were detected by reverse transcription polymerase chain reaction (RT-PCR). The HHIP methylation status of the promoter region in the gastric carcinnoma tissues and peritumoral tissues was detected by methylation-specific PCR (MSP). Prior to and following treatment with methyl transferase inhibitor 5-aza-2'-deoxycitydine (5-aza-dc), the HHIP mRNA expression level, the methylation status of the promoter region and the methylation site loci on the CpG islands in the AGS cells were detected by RT-PCR, MSP and bisulfite sequencing PCR (BSP), respectively. The correlation between the methylation status of the CpG islands at the HHIP promoter region and the HHIP mRNA expression level were analyzed. It was found that the expression level of the HHIP mRNA in the gastric carcinoma tissues was significantly lower than that in the adjacent tissues (0.82±0.38 vs. 1.60±0.26, respectively; P<0.001). No significant correlations were observed between the expression of HHIP mRNA and age, gender, tumor-node-metastasis stage, differentiation degree and presence of lymph node metastasis (P>0.05). The degree of methylation of the HHIP gene promotor in the peritumoral tissues (17.7±3.59%) was significantly lower than that in the gastric cancer tissues (62.9±6.14%) and in the AGS cells (99.7±0.67%) (P<0.05). Compared with prior to 5-aza-dc intervention, the HHIP mRNA expression level in the AGS cells was significantly increased subsequent to intervention (0.21±0.12 vs. 4.68±0.22; P<0.01), while the degree of methylation in the AGS cells was significantly decreased (90.2±0.67 vs. 10.1±0.21%; P<0.01), and the methylation sites in CpG islands were significantly reduced. The degree of HHIP methylation showed a negative correlation with the level of mRNA expression (r=-0.693; P<0.01). It can be hypothesized that a high degree of methylation of the HHIP gene promoter CpG islands in gastric cancer tissues and cells causes a decrease in HHIP mRNA expression, which may be involved in the carcinogenesis of gastric cancer.

New answers to old questions from genome-wide maps of DNA methylation in hematopoietic cells

DNA methylation is a well-studied epigenetic modification essential for efficient cellular differentiation. Aberrant DNA methylation patterns are a characteristic feature of cancer, including myeloid malignancies such as acute myeloid leukemia. Recurrent mutations in DNA-modifying enzymes were identified in acute myeloid leukemia and linked to distinct DNA methylation signatures. In addition, discovery of Tet enzymes provided new mechanisms for the reversal of DNA methylation. Advances in base-resolution profiling of DNA methylation have enabled a more comprehensive understanding of the methylome landscape in the genome. This review will summarize and discuss the key questions in the function of DNA methylation in the hematopoietic system, including where and how DNA methylation regulates diverse biological processes in the genome as elucidated by recent studies.

DNA methylation dynamics during epigenetic reprogramming in the germline and preimplantation embryos

Methylation of DNA is an essential epigenetic control mechanism in mammals. Messerschmidt et al. review the current understanding of epigenetic dynamics regulating the molecular processes that prepare the mammalian embryo for normal development.

Methylation of DNA is an essential epigenetic control mechanism in mammals. During embryonic development, cells are directed toward their future lineages, and DNA methylation poses a fundamental epigenetic barrier that guides and restricts differentiation and prevents regression into an undifferentiated state. DNA methylation also plays an important role in sex chromosome dosage compensation, the repression of retrotransposons that threaten genome integrity, the maintenance of genome stability, and the coordinated expression of imprinted genes. However, DNA methylation marks must be globally removed to allow for sexual reproduction and the adoption of the specialized, hypomethylated epigenome of the primordial germ cell and the preimplantation embryo. Recent technological advances in genome-wide DNA methylation analysis and the functional description of novel enzymatic DNA demethylation pathways have provided significant insights into the molecular processes that prepare the mammalian embryo for normal development.

Widespread differences in cortex DNA methylation of the "language gene" CNTNAP2 between humans and chimpanzees

CNTNAP2, one of the largest genes in the human genome, has been linked to human-specific language abilities and neurodevelopmental disorders. Our hypothesis is that epigenetic rather than genetic changes have accelerated the evolution of the human brain. To compare the cortex DNA methylation patterns of human and chimpanzee CNTNAP2 at ultra-high resolution, we combined methylated DNA immunoprecipitation (MeDIP) with NimbleGen tiling arrays for the orthologous gene and flanking sequences. Approximately 1.59 Mb of the 2.51 Mb target region could be aligned and analyzed with a customized algorithm in both species. More than one fifth (0.34 Mb) of the analyzed sequence throughout the entire gene displayed significant methylation differences between six human and five chimpanzee cortices. One of the most striking interspecies differences with 28% methylation in human and 59% in chimpanzee cortex (by bisulfite pyrosequencing) lies in a region 300 bp upstream of human SNP rs7794745 which has been associated with autism and parent-of-origin effects. Quantitative real-time RT PCR revealed that the protein-coding splice variant CNTNAP2-201 is 1.6-fold upregulated in human cortex, compared with the chimpanzee. Transcripts CNTNAP2-001, -002, and -003 did not show skewed allelic expression, which argues against CNTNAP2 imprinting, at least in adult human brain. Collectively, our results suggest widespread cortex DNA methylation changes in CNTNAP2 since the human-chimpanzee split, supporting a role for CNTNAP2 fine-regulation in human-specific language and communication traits.

Cross-species analysis of genic GC3 content and DNA methylation patterns

The GC content in the third codon position (GC₃) exhibits a unimodal distribution in many plant and animal genomes. Interestingly, grasses and homeotherm vertebrates exhibit a unique bimodal distribution. High GC₃ was previously found to be associated with variable expression, higher frequency of upstream TATA boxes, and an increase of GC₃ from 5′ to 3′. Moreover, GC₃-rich genes are predominant in certain gene classes and are enriched in CpG dinucleotides that are potential targets for methylation. Based on the GC₃ bimodal distribution we hypothesize that GC₃ has a regulatory role involving methylation and gene expression. To test that hypothesis, we selected diverse taxa (rice, thale cress, bee, and human) that varied in the modality of their GC₃ distribution and tested the association between GC₃, DNA methylation, and gene expression. We examine the relationship between cytosine methylation levels and GC₃, gene expression, genome signature, gene length, and other gene compositional features. We find a strong negative correlation (Pearson’s correlation coefficient r = −0.67, P value < 0.0001) between GC₃ and genic CpG methylation. The comparison between 5′-3′ gradients of CG₃-skew and genic methylation for the taxa in the study suggests interplay between gene-body methylation and transcription-coupled cytosine deamination effect. Compositional features are correlated with methylation levels of genes in rice, thale cress, human, bee, and fruit fly (which acts as an unmethylated control). These patterns allow us to generate evolutionary hypotheses about the relationships between GC₃ and methylation and how these affect expression patterns. Specifically, we propose that the opposite effects of methylation and compositional gradients along coding regions of GC₃-poor and GC₃-rich genes are the products of several competing processes.

DNA methylation is associated with an increased level of conservation at nondegenerate nucleotides in mammals

DNA methylation at CpG dinucleotides can significantly increase the rate of cytosine-to-thymine mutations and the level of sequence divergence. Although the correlations between DNA methylation and genomic sequence evolution have been widely studied, an unaddressed yet fundamental question is how DNA methylation is associated with the conservation of individual nucleotides in different sequence contexts. Here, we demonstrate that in mammalian exons, the correlations between DNA methylation and the conservation of individual nucleotides are dependent on the type of exonic sequence (coding or untranslated), the degeneracy of coding nucleotides, background selection pressure, and the relative position (first or nonfirst exon in the transcript) where the nucleotides are located. For untranslated and nonzero-fold degenerate nucleotides, methylated sites are less conserved than unmethylated sites regardless of background selection pressure and the relative position of the exon. For zero-fold degenerate (or nondegenerate) nucleotides, however, the reverse trend is observed in nonfirst coding exons and first coding exons that are under stringent background selection pressure. Furthermore, cytosine-to-thymine mutations at methylated zero-fold degenerate nucleotides are predicted to be more detrimental than those that occur at unmethylated nucleotides. As zero-fold and nonzero-fold degenerate nucleotides are very close to each other, our results suggest that the "functional resolution" of DNA methylation may be finer than previously recognized. In addition, the positive correlation between CpG methylation and the level of conservation at zero-fold degenerate nucleotides implies that CpG methylation may serve as an "indicator" of functional importance of these nucleotides.

The emerging role of epigenetics in rheumatic diseases

Epigenetics is a key mechanism regulating the expression of genes. There are three main and interrelated mechanisms: DNA methylation, post-translational modification of histone proteins and non-coding RNA. Gene activation is generally associated with lower levels of DNA methylation in promoters and with distinct histone marks such as acetylation of amino acids in histones. Unlike the genetic code, the epigenome is altered by endogenous (e.g. hormonal) and environmental (e.g. diet, exercise) factors and changes with age. Recent evidence implicates epigenetic mechanisms in the pathogenesis of common rheumatic disease, including RA, OA, SLE and scleroderma. Epigenetic drift has been implicated in age-related changes in the immune system that result in the development of a pro-inflammatory status termed inflammageing, potentially increasing the risk of age-related conditions such as polymyalgia rheumatica. Therapeutic targeting of the epigenome has shown promise in animal models of rheumatic diseases. Rapid advances in computational biology and DNA sequencing technology will lead to a more comprehensive understanding of the roles of epigenetics in the pathogenesis of common rheumatic diseases.

Intragenic DNA methylation modulates alternative splicing by recruiting MeCP2 to promote exon recognition

Although the function of DNA methylation in gene promoter regions is well established in transcriptional repression, the function of the evolutionarily conserved widespread distribution of DNA methylation in gene body regions remains incompletely understood. Here, we show that DNA methylation is enriched in included alternatively spliced exons (ASEs), and that inhibition of DNA methylation results in aberrant splicing of ASEs. The methyl-CpG-binding protein MeCP2 is enriched in included ASEs, particularly those that are also highly methylated, and inhibition of DNA methylation disrupts specific targeting of MeCP2 to exons. Interestingly, ablation of MeCP2 results in increased histone acetylation and aberrant ASE-skipping events. We further show that inhibition of histone deacetylase (HDAC) activity leads to exon skipping that shows a highly significant degree of overlap with that caused by MeCP2 knockdown. Together, our data indicate that intragenic DNA methylation operates in exon definition to modulate alternative RNA splicing and can enhance exon recognition via recruitment of the multifunctional protein MeCP2, which thereby maintains local histone hypoacetylation through the subsequent recruitment of HDACs.

Intragenic DNA methylation in transcriptional regulation, normal differentiation and cancer

Ever since the discovery of DNA methylation at cytosine residues, the role of this so called fifth base has been extensively studied and debated. Until recently, the majority of DNA methylation studies focused on the analysis of CpG islands associated to promoter regions. However, with the upcoming possibilities to study DNA methylation in a genome-wide context, this epigenetic mark can now be studied in an unbiased manner. As a result, recent studies have shown that not only promoters but also intragenic and intergenic regions are widely modulated during physiological processes and disease. In particular, it is becoming increasingly clear that DNA methylation in the gene body is not just a passive witness of gene transcription but it seems to be actively involved in multiple gene regulation processes. In this review we discuss the potential role of intragenic DNA methylation in alternative promoter usage, regulation of short and long non-coding RNAs, alternative RNA processing, as well as enhancer activity. Furthermore, we summarize how the intragenic DNA methylome is modified both during normal cell differentiation and neoplastic transformation.

Integrative analysis of tissue-specific methylation and alternative splicing identifies conserved transcription factor binding motifs

The exact role of intragenic DNA methylation in regulating tissue-specific gene regulation is unclear. Recently, the DNA-binding protein CTCF has been shown to participate in the regulation of alternative splicing in a DNA methylation-dependent manner. To globally evaluate the relationship between DNA methylation and tissue-specific alternative splicing, we performed genome-wide DNA methylation profiling of mouse retina and brain. In protein-coding genes, tissue-specific differentially methylated regions (T-DMRs) were preferentially located in exons and introns. Gene ontology and evolutionary conservation analysis suggest that these T-DMRs are likely to be biologically relevant. More than 14% of alternatively spliced genes were associated with a T-DMR. T-DMR-associated genes were enriched for developmental genes, suggesting that a specific set of alternatively spliced genes may be regulated through DNA methylation. Novel DNA sequences motifs overrepresented in T-DMRs were identified as being associated with positive and/or negative regulation of alternative splicing in a position-dependent context. The majority of these evolutionarily conserved motifs contain a CpG dinucleotide. Some transcription factors, which recognize these motifs, are known to be involved in splicing. Our results suggest that DNA methylation-dependent alternative splicing is widespread and lay the foundation for further mechanistic studies of the role of DNA methylation in tissue-specific splicing regulation.

Detecting splicing variants in idiopathic pulmonary fibrosis from non-differentially expressed genes

Idiopathic pulmonary fibrosis (IPF) is an interstitial lung disease of unknown cause that lacks a proven therapy for altering its high mortality rate. Microarrays have been employed to investigate the pathogenesis of IPF, but are presented mostly at the gene-expression level due to technologic limitations. In as much as, alternative RNA splicing isoforms are increasingly identified as potential regulators of human diseases, including IPF, we propose a new approach with the capacity to detect splicing variants using RNA-seq data. We conducted a joint analysis of differential expression and differential splicing on annotated human genes and isoforms, and identified 122 non-differentially expressed genes with a high degree of "switch" between major and minor isoforms. Three cases with variant mechanisms for alternative splicing were validated using qRT-PCR, among the group of genes in which expression was not significantly changed at the gene level. We also identified 35 novel transcripts that were unique to the fibrotic lungs using exon-exon junction evidence, and selected a representative for qRT-PCR validation. The results of our study are likely to provide new insight into the pathogenesis of pulmonary fibrosis and may eventuate in new treatment targets.

Genetic and epigenetic variants contributing to clofarabine cytotoxicity

2-chloro-2-fluoro-deoxy-9-D-arabinofuranosyladenine (Clofarabine), a purine nucleoside analog, is used in the treatment of hematologic malignancies and as induction therapy for stem cell transplantation. The discovery of pharmacogenomic markers associated with chemotherapeutic efficacy and toxicity would greatly benefit the utility of this drug. Our objective was to identify genetic and epigenetic variants associated with clofarabine toxicity using an unbiased, whole genome approach. To this end, we employed International HapMap lymphoblastoid cell lines (190 LCLs) of European (CEU) or African (YRI) ancestry with known genetic information to evaluate cellular sensitivity to clofarabine. We measured modified cytosine levels to ascertain the contribution of genetic and epigenetic factors influencing clofarabine-mediated cytotoxicity. Association studies revealed 182 single nucleotide polymorphisms (SNPs) and 143 modified cytosines associated with cytotoxicity in both populations at the threshold P ≤ 0.0001. Correlation between cytotoxicity and baseline gene expression revealed 234 genes at P ≤ 3.98 × 10(-6). Six genes were implicated as: (i) their expression was directly correlated to cytotoxicity, (ii) they had a targeting SNP associated with cytotoxicity, and (iii) they had local modified cytosines associated with gene expression and cytotoxicity. We identified a set of three SNPs and three CpG sites targeting these six genes explaining 43.1% of the observed variation in phenotype. siRNA knockdown of the top three genes (SETBP1, BAG3, KLHL6) in LCLs revealed altered susceptibility to clofarabine, confirming relevance. As clofarabine's toxicity profile includes acute kidney injury, we examined the effect of siRNA knockdown in HEK293 cells. siSETBP1 led to a significant change in HEK293 cell susceptibility to clofarabine.

Extensive alternative splicing of the repressor element silencing transcription factor linked to cancer

The repressor element silencing transcription factor (REST) is a coordinate transcriptional and epigenetic regulator which functions as a tumor suppressor or an oncogene depending on cellular context, and a truncated splice variant REST4 has been linked to various types of cancer. We performed a comprehensive analysis of alternative splicing (AS) of REST by rapid amplification of cDNA ends and PCR amplification of cDNAs from various tissues and cell lines with specific primers. We identified 8 novel alternative exons including an alternate last exon which doubles the REST gene boundary, along with numerous 5'/3' splice sites and ends in the constitutive exons. With the combination of various splicing patterns (e.g. exon skipping and alternative usage of the first and last exons) that are predictive of altered REST activity, at least 45 alternatively spliced variants of coding and non-coding mRNA were expressed in a species- and cell-type/tissue-specific manner with individual differences. By examining the repertoire of REST pre-mRNA splicing in 27 patients with kidney, liver and lung cancer, we found that all patients without exception showed differential expression of various REST splice variants between paired tumor and adjacent normal tissues, with striking cell-type/tissue and individual differences. Moreover, we revealed that exon 3 skipping, which causes no frame shift but loss of a domain essential for nuclear translocation, was affected by pioglitazone, a highly selective activator of the peroxisome proliferator-activated receptor gamma (PPARγ) which contributes to cell differentiation and tumorigenesis besides its metabolic actions. Accordingly, this study demonstrates an extensive AS of REST pre-mRNA which redefines REST gene boundary and structure, along with a general but differential link between REST pre-mRNA splicing and various types of cancer. These findings advance our understanding of the complex, context-dependent regulation of REST gene expression and function, and provide potential biomarkers and therapeutic targets for cancer.

DNA-methylation effect on cotranscriptional splicing is dependent on GC architecture of the exon-intron structure

DNA methylation is known to regulate transcription and was recently found to be involved in exon recognition via cotranscriptional splicing. We recently observed that exon–intron architectures can be grouped into two classes: one with higher GC content in exons compared to the flanking introns, and the other with similar GC content in exons and introns. The first group has higher nucleosome occupancy on exons than introns, whereas the second group exhibits weak nucleosome marking of exons, suggesting another type of epigenetic marker distinguishes exons from introns when GC content is similar. We find different and specific patterns of DNA methylation in each of the GC architectures; yet in both groups, DNA methylation clearly marks the exons. Exons of the leveled GC architecture exhibit a significantly stronger DNA methylation signal in relation to their flanking introns compared to exons of the differential GC architecture. This is accentuated by a reduction of the DNA methylation level in the intronic sequences in proximity to the splice sites and shows that different epigenetic modifications mark the location of exons already at the DNA level. Also, lower levels of methylated CpGs on alternative exons can successfully distinguish alternative exons from constitutive ones. Three positions at the splice sites show high CpG abundance and accompany elevated nucleosome occupancy in a leveled GC architecture. Overall, these results suggest that DNA methylation affects exon recognition and is influenced by the GC architecture of the exon and flanking introns.