MethTools--a toolbox to visualize and analyze DNA methylation data

ABSP: an automated R tool to efficiently analyze region-specific CpG methylation from bisulfite sequencing PCR

MOTIVATION

Nowadays, epigenetic gene regulations are studied in each part of the biology, from embryonic development to diseases such as cancers and neurodegenerative disorders. Currently, to quantify and compare CpG methylation levels of a specific region of interest, the most accessible technique is the bisulfite sequencing PCR (BSP). However, no existing user-friendly tool is able to analyze data from all approaches of BSP. Therefore, the most convenient way to process results from the direct sequencing of PCR products (direct-BSP) is to manually analyze the chromatogram traces, which is a repetitive and prone to error task.

RESULTS

Here, we implement a new R-based tool, called ABSP for analysis of bisulfite sequencing PCR, providing a complete analytic process of both direct-BSP and cloning-BSP data. It uses the raw sequencing trace files (.ab1) as input to compute and compare CpG methylation percentages. It is fully automated and includes a user-friendly interface as a built-in R shiny app, quality control steps and generates publication-ready graphics.

AVAILABILITY AND IMPLEMENTATION

The ABSP tool and associated data are available on GitHub at https://github.com/ABSP-methylation-tool/ABSP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Novel role of UHRF1 in the epigenetic repression of the latent HIV-1

BACKGROUND

The multiplicity, heterogeneity, and dynamic nature of human immunodeficiency virus type-1 (HIV-1) latency mechanisms are reflected in the current lack of functional cure for HIV-1. Accordingly, all classes of latency-reversing agents (LRAs) have been reported to present variable ex vivo potencies. Here, we investigated the molecular mechanisms underlying the potency variability of one LRA: the DNA methylation inhibitor 5-aza-2'-deoxycytidine (5-AzadC).

METHODS

We employed epigenetic interrogation methods (electrophoretic mobility shift assays, chromatin immunoprecipitation, Infinium array) in complementary HIV-1 infection models (latently-infected T-cell line models, primary CD4+ T-cell models and ex vivo cultures of PBMCs from HIV+ individuals). Extracellular staining of cell surface receptors and intracellular metabolic activity were measured in drug-treated cells. HIV-1 expression in reactivation studies was explored by combining the measures of capsid p24Gag protein, green fluorescence protein signal, intracellular and extracellular viral RNA and viral DNA.

FINDINGS

We uncovered specific demethylation CpG signatures induced by 5-AzadC in the HIV-1 promoter. By analyzing the binding modalities to these CpG, we revealed the recruitment of the epigenetic integrator Ubiquitin-like with PHD and RING finger domain 1 (UHRF1) to the HIV-1 promoter. We showed that UHRF1 redundantly binds to the HIV-1 promoter with different binding modalities where DNA methylation was either non-essential, essential or enhancing UHRF1 binding. We further demonstrated the role of UHRF1 in the epigenetic repression of the latent viral promoter by a concerted control of DNA and histone methylations.

INTERPRETATION

A better understanding of the molecular mechanisms of HIV-1 latency allows for the development of innovative antiviral strategies. As a proof-of-concept, we showed that pharmacological inhibition of UHRF1 in ex vivo HIV+ patient cell cultures resulted in potent viral reactivation from latency. Together, we identify UHRF1 as a novel actor in HIV-1 epigenetic silencing and highlight that it constitutes a new molecular target for HIV-1 cure strategies.

FUNDING

Funding was provided by the Belgian National Fund for Scientific Research (F.R.S.-FNRS, Belgium), the « Fondation Roi Baudouin », the NEAT (European AIDS Treatment Network) program, the Internationale Brachet Stiftung, ViiV Healthcare, the Télévie, the Walloon Region (« Fonds de Maturation »), « Les Amis des Instituts Pasteur à Bruxelles, asbl », the University of Brussels (Action de Recherche Concertée ULB grant), the Marie Skodowska Curie COFUND action, the European Union's Horizon 2020 research and innovation program under grant agreement No 691119-EU4HIVCURE-H2020-MSCA-RISE-2015, the French Agency for Research on AIDS and Viral Hepatitis (ANRS), the Sidaction and the "Alsace contre le Cancer" Foundation. This work is supported by 1UM1AI164562-01, co-funded by National Heart, Lung and Blood Institute, National Institute of Diabetes and Digestive and Kidney Diseases, National Institute of Neurological Disorders and Stroke, National Institute on Drug Abuse and the National Institute of Allergy and Infectious Diseases.

Transcriptional behavior of the HIV-1 promoter in context of the BACH2 prominent proviral integration gene

Chronic infection with human immunodeficiency virus (HIV)-1 is characterized by accumulation of proviral sequences in the genome of target cells. Integration of viral DNA in patients on long-term antiretroviral therapy selectively persists at preferential loci, suggesting site-specific crosstalk of viral sequences and human genes. This crosstalk likely contributes to chronic HIV disease through modulation of host immune pathways and emergence of clonal infected cell populations. To systematically interrogate such effects, we undertook genome engineering to generate Jurkat cell models that replicate integration of HIV-1 long terminal repeat (LTR) sequences at the BTB and CNC Homolog 2 (BACH2) integration locus. This locus is a prominent HIV-1 integration gene in chronic infection, found in 30 % of long-term treated patients with mapped proviral integrations. Using five clonal models carrying an LTR-driven reporter at different BACH2 intergenic regions, we here show that LTR transcriptional activity is repressed in BACH2 regions associated with proviral-DNA integrations in vivo but not in a control region. Our data indicates that this repression is in part epigenetically regulated, particularly through DNA methylation. Importantly, we demonstrate that transcriptional activity of the LTR is independent of BACH2 gene transcription and vice versa in our models. This suggests no transcriptional interference of endogenous and HIV-1 promoters. Taken together, our study provides first insights into how activity of HIV-1 LTR sequences is regulated at the BACH2 locus as prominent example for a recurrently-detected integration gene in chronic infection. Given the importance of integration-site dependent virus/host crosstalk for chronic HIV disease, our findings for the BACH2 locus have potential implications for future therapeutic strategies.

Quantitative Analysis of DNA Methylation by Bisulfite Sequencing

Changes in deoxyribonucleic acid (DNA) methylation are shown to occur with aging in mammals. Besides changes that seem to be essentially stochastic, methylation levels of certain CpG sites display a strong correlation with age. Collectively, methylation of such CpG sites could be used as "epigenetic clocks" to predict biological age. Numerous versions of the epigenetic clock have been proposed, all of them based on quantitative estimation of the methylation levels of individual CpG sites. Different methods were elaborated for quantitative measurements of DNA methylation, with the most reliable of these based on bisulfite treatment of DNA. We present here a protocol for assessment of the methylation levels of individual CpG sites in target DNA sequences by the direct sequencing of polymerase chain reaction (PCR) amplification products obtained from bisulfate-converted DNA.

Epigenetic Modifications Following Noxious Stimuli in Infants

Purpose:

To recruit healthy full- and preterm infants into genetic research and determine the effectiveness of a noninvasive DNA sampling technique for comparing epigenetic modifications.

Background:

Noxious stimuli during a vulnerable period of infant neuronal plasticity may trigger long-term epigenetic changes affecting neurodevelopment, pain modulation, and reactivity. Recognizing epigenetic pain findings is problematic because parents are reluctant to enroll newborns into genetic research.

Methods:

Design: Within-subject change over time candidate-gene DNA methylation association study. Setting/ sample: Urban teaching hospital’s neonatal intensive care unit and newborn nursery. Convenience sample of healthy full- (>37 weeks, n = 6) and preterm (<37 weeks, n = 6) infants. Procedure: Parents participated in a genetic presentation prior to informed consent. Infant buccal saliva was collected after admission to the unit and prior to discharge. Analysis: The methylation pattern at the 5′ end of µ-opioid receptor gene ( OPRM1) was examined. DNA was treated with bisulfite to convert all cytosines to uracil residues, leaving methylated cytosines unchanged. Sequencing of untreated and bisulfite-converted DNA was carried out. The sequences of unconverted and bisulfite-converted DNA were aligned with ClustalW, fidelity of the polymerase chain reaction and the sequencing reaction evaluated, and the methylation pattern identified.

Results:

Recruitment and assessment of a noninvasive DNA sampling technique for comparing epigenetic modifications were successful; however, infant stress did not produce a change in OPRM1 methylation expression.

Relevance:

This study established the feasibility of recruiting healthy full-term infants into genetic research and the effectiveness of noninvasive DNA sampling for comparing epigenetic modification in infants.

Increased Expression of Plasma-Induced ABCC1 mRNA in Cystic Fibrosis

The ABCC1 gene is structurally and functionally related to the cystic fibrosis transmembrane conductance regulator gene (CFTR). Upregulation of ABCC1 is thought to improve lung function in patients with cystic fibrosis (CF); the mechanism underlying this effect is unknown. We analyzed the ABCC1 promoter single nucleotide polymorphism (SNP rs504348), plasma-induced ABCC1 mRNA expression levels, and ABCC1 methylation status and their correlation with clinical variables among CF subjects with differing CFTR mutations. We assigned 93 CF subjects into disease severity groups and genotyped SNP rs504348. For 23 CF subjects and 7 healthy controls, donor peripheral blood mononuclear cells (PBMCs) stimulated with plasma underwent gene expression analysis via qRT-PCR. ABCC1 promoter methylation was analyzed in the same 23 CF subjects. No significant correlation was observed between rs504348 genotypes and CF disease severity, but pancreatic insufficient CF subjects showed increased colonization with any form of Pseudomonas aeruginosa (OR = 3.125, 95% CI: 1.192-8.190) and mucoid P. aeruginosa (OR = 5.075, 95% CI: 1.307-28.620) compared to the pancreatic sufficient group. A significantly higher expression of ABCC1 mRNA was induced by CF plasma compared to healthy control plasma (p < 0.001). CF subjects with rs504348 (CC/CG) also had higher mRNA expression compared to those with the ancestral GG genotype (p < 0.005). ABCC1 promoter was completely unmethylated; therefore, we did not detect any association between methylation and CF disease severity. In silico predictions suggested that histone modifications are crucial for regulating ABCC1 expression in PBMCs. Our results suggest that ABCC1 expression has a role in CFTR activity thereby increasing our understanding of the molecular underpinnings of the clinical heterogeneity in CF.

The Biomphalaria glabrata DNA methylation machinery displays spatial tissue expression, is differentially active in distinct snail populations and is modulated by interactions with Schistosoma mansoni

BACKGROUND

The debilitating human disease schistosomiasis is caused by infection with schistosome parasites that maintain a complex lifecycle alternating between definitive (human) and intermediate (snail) hosts. While much is known about how the definitive host responds to schistosome infection, there is comparably less information available describing the snail's response to infection.

METHODOLOGY/PRINCIPLE FINDINGS

Here, using information recently revealed by sequencing of the Biomphalaria glabrata intermediate host genome, we provide evidence that the predicted core snail DNA methylation machinery components are associated with both intra-species reproduction processes and inter-species interactions. Firstly, methyl-CpG binding domain protein (Bgmbd2/3) and DNA methyltransferase 1 (Bgdnmt1) genes are transcriptionally enriched in gonadal compared to somatic tissues with 5-azacytidine (5-AzaC) treatment significantly inhibiting oviposition. Secondly, elevated levels of 5-methyl cytosine (5mC), DNA methyltransferase activity and 5mC binding in pigmented hybrid- compared to inbred (NMRI)- B. glabrata populations indicate a role for the snail's DNA methylation machinery in maintaining hybrid vigour or heterosis. Thirdly, locus-specific detection of 5mC by bisulfite (BS)-PCR revealed 5mC within an exonic region of a housekeeping protein-coding gene (Bg14-3-3), supporting previous in silico predictions and whole genome BS-Seq analysis of this species' genome. Finally, we provide preliminary evidence for parasite-mediated host epigenetic reprogramming in the schistosome/snail system, as demonstrated by the increase in Bgdnmt1 and Bgmbd2/3 transcript abundance following Bge (B. glabrata embryonic cell line) exposure to parasite larval transformation products (LTP).

CONCLUSIONS/SIGNIFICANCE

The presence of a functional DNA methylation machinery in B. glabrata as well as the modulation of these gene products in response to schistosome products, suggests a vital role for DNA methylation during snail development/oviposition and parasite interactions. Further deciphering the role of this epigenetic process during Biomphalaria/Schistosoma co-evolutionary biology may reveal key factors associated with disease transmission and, moreover, enable the discovery of novel lifecycle intervention strategies.

Small-Magnitude Effect Sizes in Epigenetic End Points are Important in Children's Environmental Health Studies: The Children's Environmental Health and Disease Prevention Research Center's Epigenetics Working Group

BACKGROUND

Characterization of the epigenome is a primary interest for children's environmental health researchers studying the environmental influences on human populations, particularly those studying the role of pregnancy and early-life exposures on later-in-life health outcomes.

OBJECTIVES

Our objective was to consider the state of the science in environmental epigenetics research and to focus on DNA methylation and the collective observations of many studies being conducted within the Children's Environmental Health and Disease Prevention Research Centers, as they relate to the Developmental Origins of Health and Disease (DOHaD) hypothesis.

METHODS

We address the current laboratory and statistical tools available for epigenetic analyses, discuss methods for validation and interpretation of findings, particularly when magnitudes of effect are small, question the functional relevance of findings, and discuss the future for environmental epigenetics research.

DISCUSSION

A common finding in environmental epigenetic studies is the small-magnitude epigenetic effect sizes that result from such exposures. Although it is reasonable and necessary that we question the relevance of such small effects, we present examples in which small effects persist and have been replicated across populations and across time. We encourage a critical discourse on the interpretation of such small changes and further research on their functional relevance for children's health.

CONCLUSION

The dynamic nature of the epigenome will require an emphasis on future longitudinal studies in which the epigenome is profiled over time, over changing environmental exposures, and over generations to better understand the multiple ways in which the epigenome may respond to environmental stimuli.

Epigenetic mechanisms of peptidergic regulation of gene expression during aging of human cells

Expression levels of genes encoding specific transcription factors and other functionally important proteins vary upon aging of pancreatic and bronchial epithelium cell cultures. The peptides KEDW and AEDL tissue-specifically affect gene expression in pancreatic and bronchial cell cultures, respectively. It is established in this work that the DNA methylation patterns of the PDX1, PAX6, NGN3, NKX2-1, and SCGB1A1 gene promoter regions change upon aging in pancreatic and bronchial cell cultures in correlation with variations in their expression levels. Thus, stable changes in gene expression upon aging of cell cultures could be caused by changes in their promoter methylation patterns. The methylation patterns of the PAX4 gene in pancreatic cells as well as those of the FOXA1, SCGB3A2, and SFTPA1 genes in bronchial cells do not change upon aging and are unaffected by peptides, whereas their expression levels change in both cases. The promoter region of the FOXA2 gene in pancreatic cells contains a small number of methylated CpG sites, their methylation levels being affected by cell culture aging and KEDW, though without any correlation with gene expression levels. The promoter region of the FOXA2 gene is completely unmethylated in bronchial cells irrespective of cell culture age and AEDL action. Changes in promoter methylation might be the cause of age- and peptide-induced variations in expression levels of the PDX1, PAX6, and NGN3 genes in pancreatic cells and NKX2-1 and SCGB1A1 genes in bronchial cells. Expression levels of the PAX4 and FOXA2 genes in pancreatic cells and FOXA1, FOXA2, SCGB3A2, and SFTPA1 genes in bronchial cells seem to be controlled by some other mechanisms.

Valproic Acid Increases the Hepatic Differentiation Potential of Salivary Gland Cells

The studies of cell plasticity and differentiation abilities are important problems in modern cellular biology. The use of histone deacetylase inhibitor - valproic acid is a promising approach to increasing the differentiation efficiency of various cell types. In this paper we investigate the ability of mouse submandibular salivary gland cells to differentiate into the hepatic direction and the effect of valproic acid on the efficiency of this differentiation. It was shown that the gene expression levels of hepatocyte markers (Aat, Afp, G6p, Pepck, Tat, Cyp3a13) and liver-enriched transcription factors (Hnf-3α, Hnf-3β, Hnf-4α, Hnf-6) were increased after differentiation in salivary gland cells. Valproic acid increases the specificity of hepatic differentiation, reducing the expression levels of the ductal (Krt19, Hhex1, Cyp7a1) and acinar (Ptf1a) markers. After valproic acid exposure, the efficiency of hepatic differentiation also increases, as evidenced by the increase in the gene expression level of Alb and Tdo, and increase in urea production by differentiated cells. No change was found in DNA methylation of the promoter regions of the genes; however, valproic acid treatment and subsequent hepatic differentiation largely affected the histone H3 methylation of liver-enriched genes. Thus, mouse submandibular salivary gland cells are capable of effective differentiation in the hepatic direction. Valproic acid increases the specificity and efficiency of the hepatic differentiation of these cells.

PGC-1α Promoter Methylation in Parkinson's Disease

The etiopathogenesis of sporadic Parkinson's disease (PD) remains elusive although mitochondrial dysfunction has long been implicated. Recent evidence revealed reduced expression of peroxisome proliferator-activated receptor gamma coactivator-1 α (PGC-1α) and downstream regulated nuclear encoded respiratory complex genes in affected brain tissue from PD patients. We sought to determine whether epigenetic modification of the PGC-1α gene could account for diminished expression. In substantia nigra from PD patients but not control subjects, we show significant promoter-proximal non-canonical cytosine methylation of the PGC-1α gene but not an adjacent gene. As neuroinflammation is a prominent feature of PD and a mediator of epigenetic change, we evaluated whether the pro-inflammatory fatty acid, palmitate, would stimulate PGC-1α promoter methylation in different cell types from the CNS. Indeed, in mouse primary cortical neurons, microglia and astrocytes, palmitate causes PGC-1α gene promoter non-canonical cytosine methylation, reduced expression of the gene and reduced mitochondrial content. Moreover, intracerebroventricular (ICV) injection of palmitate to transgenic human α-synuclein mutant mice resulted in increased PGC-1α promoter methylation, decreased PGC-1α expression and reduced mitochondrial content in substantia nigra. Finally we provide evidence that dysregulation of ER stress and inflammatory signaling is associated with PGC-1α promoter methylation. Together, these data strengthen the connection between saturated fatty acids, neuroflammation, ER stress, epigenetic alteration and bioenergetic compromise in PD.

Methylation plotter: a web tool for dynamic visualization of DNA methylation data

Methylation plotter is a Web tool that allows the visualization of methylation data in a user-friendly manner and with publication-ready quality. The user is asked to introduce a file containing the methylation status of a genomic region. This file can contain up to 100 samples and 100 CpGs. Optionally, the user can assign a group for each sample (i.e. whether a sample is a tumoral or normal tissue). After the data upload, the tool produces different graphical representations of the results following the most commonly used styles to display this type of data. They include an interactive plot that summarizes the status of every CpG site and for every sample in lollipop or grid styles. Methylation values ranging from 0 (unmethylated) to 1 (fully methylated) are represented using a gray color gradient. A practical feature of the tool allows the user to choose from different types of arrangement of the samples in the display: for instance, sorting by overall methylation level, by group, by unsupervised clustering or just following the order in which data were entered.

In addition to the detailed plot, Methylation plotter produces a methylation profile plot that summarizes the status of the scrutinized region, a boxplot that sums up the differences between groups (if any) and a dendrogram that classifies the data by unsupervised clustering. Coupled with this analysis, descriptive statistics and testing for differences at both CpG and group levels are provided.

The implementation is based in R/shiny, providing a highly dynamic user interface that generates quality graphics without the need of writing R code. Methylation plotter is freely available at http://gattaca.imppc.org:3838/methylation_plotter/.

SATB1 overexpression promotes malignant T-cell proliferation in cutaneous CD30+ lymphoproliferative disease by repressing p21

Cutaneous CD30(+) lymphoproliferative disease (CD30(+)LPD), characterized by the presence of CD30(+) anaplastic large T cells, comprises the second most common group of cutaneous T-cell lymphoma (CTCL). However, little is known about the pathobiology of the CD30(+) lymphoma cells, as well as the mechanisms of disease progression. Here we report that Special AT-rich region binding protein 1 (SATB1), a thymocyte specific chromatin organizer, is over-expressed in CD30(+) lymphoma cells in most CD30(+)LPDs, and its expression is upregulated during disease progression. Our findings show that SATB1 silencing in CD30(+)LPD cells leads to G1 cell cycle arrest mediated by p21 activation. Using chromatin immunoprecipitation, luciferase assays, and mutational analysis, we demonstrate that SATB1 directly regulates the transcription of p21 in a p53-independent manner. Moreover, DNA demethylation on a specific CpG-rich region of the SATB1 promoter is associated with the upregulation of SATB1 during disease progression. These experiments define a novel SATB1-p21 pathway in malignant CD30(+) T lymphocytes, which provides novel molecular insights into the pathogenesis of CD30(+)LPDs and possibly leads to new therapies.

Progressive 35S promoter methylation increases rapidly during vegetative development in transgenic Nicotiana attenuata plants

BACKGROUND

Genetically modified plants are widely used in agriculture and increasingly in ecological research to enable the selective manipulation of plant traits in the field. Despite their broad usage, many aspects of unwanted transgene silencing throughout plant development are still poorly understood. A transgene can be epigenetically silenced by a process called RNA directed DNA methylation (RdDM), which can be seen as a heritable loss of gene expression. The spontaneous nature of transgene silencing has been widely reported, but patterns of acquirement remain still unclear.

RESULTS

Transgenic wild tobacco plants (Nicotiana attenuata) expressing heterologous genes coding for antimicrobial peptides displayed an erratic and variable occurrence of transgene silencing. We focused on three independently transformed lines (PNA 1.2, PNA 10.1 and ICE 4.4) as they rapidly lost the expression of the resistance marker and down-regulated transgene expression by more than 200 fold after only one plant generation. Bisulfite sequencing indicated hypermethylation within the 35S and NOS promoters of these lines. To shed light on the progress of methylation establishment, we successively sampled leaf tissues from different stages during plant development and found a rapid increase in 35S promoter methylation during vegetative growth (up to 77% absolute increase within 45 days of growth). The levels of de novo methylation were inherited by the offspring without any visible discontinuation. A secondary callus regeneration step could interfere with the establishment of gene silencing and we found successfully restored transgene expression in the offspring of several regenerants.

CONCLUSIONS

The unpredictability of the gene silencing process requires a thorough selection and early detection of unstable plant lines. De novo methylation of the transgenes was acquired solely during vegetative development and did not require a generational change for its establishment or enhancement. A secondary callus regeneration step provides a convenient way to rescue transgene expression without causing undesirable morphological effects, which is essential for experiments that use transformed plants in the analysis of ecologically important traits.

5-methyl-cytosine and 5-hydroxy-methyl-cytosine in the genome of Biomphalaria glabrata, a snail intermediate host of Schistosoma mansoni

Background

Biomphalaria glabrata is the mollusc intermediate host for Schistosoma mansoni, a digenean flatworm parasite that causes human intestinal schistosomiasis. An estimated 200 million people in 74 countries suffer from schistosomiasis, in terms of morbidity this is the most severe tropical disease after malaria. Epigenetic information informs on the status of gene activity that is heritable, for which changes are reversible and that is not based on the DNA sequence. Epigenetic mechanisms generate variability that provides a source for potentially heritable phenotypic variation and therefore could be involved in the adaptation to environmental constraint. Phenotypic variations are particularly important in host-parasite interactions in which both selective pressure and rate of evolution are high. In this context, epigenetic changes are expected to be major drivers of phenotypic plasticity and co-adaptation between host and parasite. Consequently, with characterization of the genomes of invertebrates that are parasite vectors or intermediate hosts, it is also essential to understand how the epigenetic machinery functions to better decipher the interplay between host and parasite.

Methods

The CpGo/e ratios were used as a proxy to investigate the occurrence of CpG methylation in B. glabrata coding regions. The presence of DNA methylation in B. glabrata was also confirmed by several experimental approaches: restriction enzymatic digestion with isoschizomers, bisulfite conversion based techniques and LC-MS/MS analysis.

Results

In this work, we report that DNA methylation, which is one of the carriers of epigenetic information, occurs in B. glabrata; approximately 2% of cytosine nucleotides are methylated. We describe the methylation machinery of B. glabrata. Methylation occurs predominantly at CpG sites, present at high ratios in coding regions of genes associated with housekeeping functions. We also demonstrate by bisulfite treatment that methylation occurs in multiple copies of Nimbus, a transposable element.

Conclusions

This study details DNA methylation for the first time, one of the carriers of epigenetic information in B. glabrata. The general characteristics of DNA methylation that we observed in the B. glabrata genome conform to what epigenetic studies have reported from other invertebrate species.

Weight loss after gastric bypass surgery in human obesity remodels promoter methylation

DNA methylation provides a mechanism by which environmental factors can control insulin sensitivity in obesity. Here, we assessed DNA methylation in skeletal muscle from obese people before and after Roux-en-Y gastric bypass (RYGB). Obesity was associated with altered expression of a subset of genes enriched in metabolic process and mitochondrial function. After weight loss, the expression of the majority of the identified genes was normalized to levels observed in normal-weight, healthy controls. Among the 14 metabolic genes analyzed, promoter methylation of 11 genes was normalized to levels observed in the normal-weight, healthy subjects. Using bisulfite sequencing, we show that promoter methylation of PGC-1α and PDK4 is altered with obesity and restored to nonobese levels after RYGB-induced weight loss. A genome-wide DNA methylation analysis of skeletal muscle revealed that obesity is associated with hypermethylation at CpG shores and exonic regions close to transcription start sites. Our results provide evidence that obesity and RYGB-induced weight loss have a dynamic effect on the epigenome.

A comparison of the whole genome approach of MeDIP-seq to the targeted approach of the Infinium HumanMethylation450 BeadChip(®) for methylome profiling

DNA methylation is one of the most studied epigenetic marks in the human genome, with the result that the desire to map the human methylome has driven the development of several methods to map DNA methylation on a genomic scale. Our study presents the first comparison of two of these techniques - the targeted approach of the Infinium HumanMethylation450 BeadChip® with the immunoprecipitation and sequencing-based method, MeDIP-seq. Both methods were initially validated with respect to bisulfite sequencing as the gold standard and then assessed in terms of coverage, resolution and accuracy. The regions of the methylome that can be assayed by both methods and those that can only be assayed by one method were determined and the discovery of differentially methylated regions (DMRs) by both techniques was examined. Our results show that the Infinium HumanMethylation450 BeadChip® and MeDIP-seq show a good positive correlation (Spearman correlation of 0.68) on a genome-wide scale and can both be used successfully to determine differentially methylated loci in RefSeq genes, CpG islands, shores and shelves. MeDIP-seq however, allows a wider interrogation of methylated regions of the human genome, including thousands of non-RefSeq genes and repetitive elements, all of which may be of importance in disease. In our study MeDIP-seq allowed the detection of 15,709 differentially methylated regions, nearly twice as many as the array-based method (8070), which may result in a more comprehensive study of the methylome.

CpG_MPs: identification of CpG methylation patterns of genomic regions from high-throughput bisulfite sequencing data

High-throughput bisulfite sequencing is widely used to measure cytosine methylation at single-base resolution in eukaryotes. It permits systems-level analysis of genomic methylation patterns associated with gene expression and chromatin structure. However, methods for large-scale identification of methylation patterns from bisulfite sequencing are lacking. We developed a comprehensive tool, CpG_MPs, for identification and analysis of the methylation patterns of genomic regions from bisulfite sequencing data. CpG_MPs first normalizes bisulfite sequencing reads into methylation level of CpGs. Then it identifies unmethylated and methylated regions using the methylation status of neighboring CpGs by hotspot extension algorithm without knowledge of pre-defined regions. Furthermore, the conservatively and differentially methylated regions across paired or multiple samples (cells or tissues) are identified by combining a combinatorial algorithm with Shannon entropy. CpG_MPs identified large amounts of genomic regions with different methylation patterns across five human bisulfite sequencing data during cellular differentiation. Different sequence features and significantly cell-specific methylation patterns were observed. These potentially functional regions form candidate regions for functional analysis of DNA methylation during cellular differentiation. CpG_MPs is the first user-friendly tool for identifying methylation patterns of genomic regions from bisulfite sequencing data, permitting further investigation of the biological functions of genome-scale methylation patterns.

A regulatory potential of the Xist gene promoter in vole M. rossiaemeridionalis

X chromosome inactivation takes place in the early development of female mammals and depends on the Xist gene expression. The mechanisms of Xist expression regulation have not been well understood so far. In this work, we compared Xist promoter region of vole Microtus rossiaemeridionalis and other mammalian species. We observed three conserved regions which were characterized by computational analysis, DNaseI in vitro footprinting, and reporter construct assay. Regulatory factors potentially involved in Xist activation and repression in voles were determined. The role of CpG methylation in vole Xist expression regulation was established. A CTCF binding site was found in the 5' flanking region of the Xist promoter on the active X chromosome in both males and females. We suggest that CTCF acts as an insulator which defines an inactive Xist domain on the active X chromosome in voles.

Genome-wide DNA methylation studies suggest distinct DNA methylation patterns in pediatric embryonal and alveolar rhabdomyosarcomas

Rhabdomyosarcoma is the most common soft-tissue sarcoma in children. While cytogenetic abnormalities have been well characterized in this disease, aberrant epigenetic events such as DNA hypermethylation have not been described in genome-wide studies. We have analyzed the methylation status of 25,500 promoters in normal skeletal muscle, and in cell lines and tumor samples of embryonal and alveolar rhabdomyosarcoma from pediatric patients. We identified over 1,900 CpG islands that are hypermethylated in rhabdomyosarcomas relative to skeletal muscle. Genes involved in tissue development, differentiation, and oncogenesis such as DNAJA4, HES5, IRX1, BMP8A, GATA4, GATA6, ALX3, and P4HTM were hypermethylated in both RMS cell lines and primary samples, implicating aberrant DNA methylation in the pathogenesis of rhabdomyosarcoma. Furthermore, cluster analysis revealed embryonal and alveolar subtypes had distinct DNA methylation patterns, with the alveolar subtype being enriched in DNA hypermethylation of polycomb target genes. These results suggest that DNA methylation signatures may aid in the diagnosis and risk stratification of pediatric rhabdomyosarcoma and help identify new targets for therapy.

Defects in imprinting and genome-wide DNA methylation are not common in the in vitro fertilization population

OBJECTIVE

To examine an IVF cohort for imprinted and genome-wide DNA methylation abnormalities.

DESIGN

Retrospective study.

SETTING

Research laboratory.

PATIENT(S)

DNA samples from a previously described IVF cohort that comprised 66 IVF-conceived prepubertal children (IVF, n = 34; intracytoplasmic sperm injection, n = 32) and 69 matched naturally conceived controls.

INTERVENTION(S)

DNA methylation was examined at four imprinted gene loci (H19, SNRPN, KCNQ1OT1, and IGF2) and satellite 2 using methylation-sensitive quantitative polymerase chain reaction (MSQ-PCR) followed by bisulfite sequencing at H19, SNRPN, and KCNQ1OT1. Methylated DNA immunoprecipitation (MeDIP) microarray with validation using the Sequenom MassARRAY EpiTYPER(®) platform was also used.

MAIN OUTCOME MEASURE(S)

Percentage of DNA methylation by MSQ-PCR, differential methylation based on microarray signal intensity, and percentage DNA methylation as determined by Sequenom MassARRAY EpiTYPER were compared.

RESULT(S)

No differences in percentage of methylation between the IVF and control group were observed at H19, KCNQ1OT1, SNRPN, or IGF2. Absence of aberrant imprinting was confirmed using bisulfite sequencing. Methylation of satellite 2 repeats (a surrogate for global methylation) showed no difference between the IVF and control groups. MeDIP was used to screen for differences in promoter methylation. Subsequent quantification of methylation of eight candidate genes using the Sequenom MassARRAY EpiTYPER system did not reveal any differential methylation.

CONCLUSION(S)

Low-level imprinting errors are not common in the IVF population.

Evidence for non-CpG methylation in mammals

In mammals, the existence of cytosine methylation on non-CpG sequences is controversial. Here, we adapted a LuminoMetric-based Assay (LUMA) to determine global non-CpG methylation levels in rodent and human tissues. We observed that <1% cytosines in non-CpG motifs were methylated in 3T3-L1 fibroblasts, whereas 7-13% cytosines in non-CpG motifs were methylated in mouse tissues or embryonic fibroblasts. Analysis of cytosine methylation in human, rat, and mouse tissues by bisulfite sequencing revealed non-CpG methylation levels up to 7.5% of all non-CpG cytosines. These levels dropped to 1.5% when a second round of PCR was performed prior to bisulfite sequencing, providing an explanation for the common underestimation of non-CpG methylation levels. Collectively, our results provide evidence that non-CpG methylation exists at substantial levels in mammals.

DNA methylation status predicts cell type-specific enhancer activity

Cell-selective glucocorticoid receptor (GR) binding to distal regulatory elements is associated with cell type-specific regions of locally accessible chromatin. These regions can either pre-exist in chromatin (pre-programmed) or be induced by the receptor (de novo). Mechanisms that create and maintain these sites are not well understood. We observe a global enrichment of CpG density for pre-programmed elements, and implicate their demethylated state in the maintenance of open chromatin in a tissue-specific manner. In contrast, sites that are actively opened by GR (de novo) are characterized by low CpG density, and form a unique class of enhancers devoid of suppressive effect of agglomerated methyl-cytosines. Furthermore, treatment with glucocorticoids induces rapid changes in methylation levels at selected CpGs within de novo sites. Finally, we identify GR-binding elements with CpGs at critical positions, and show that methylation can affect GR-DNA interactions in vitro. The findings present a unique link between tissue-specific chromatin accessibility, DNA methylation and transcription factor binding and show that DNA methylation can be an integral component of gene regulation by nuclear receptors.

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

Bisulfite sequencing is a widely used method for measuring DNA methylation in eukaryotic genomes. The assay provides single-base pair resolution and, given sufficient sequencing depth, its quantitative accuracy is excellent. High-throughput sequencing of bisulfite-converted DNA can be applied either genome wide or targeted to a defined set of genomic loci (e.g. using locus-specific PCR primers or DNA capture probes). Here, we describe BiQ Analyzer HT (http://biq-analyzer-ht.bioinf.mpi-inf.mpg.de/), a user-friendly software tool that supports locus-specific analysis and visualization of high-throughput bisulfite sequencing data. The software facilitates the shift from time-consuming clonal bisulfite sequencing to the more quantitative and cost-efficient use of high-throughput sequencing for studying locus-specific DNA methylation patterns. In addition, it is useful for locus-specific visualization of genome-wide bisulfite sequencing data.

Heterochromatic genes undergo epigenetic changes and escape silencing in immunodeficiency, centromeric instability, facial anomalies (ICF) syndrome

Immunodeficiency, Centromeric Instability, Facial Anomalies (ICF) syndrome is a rare autosomal recessive disorder that is characterized by a marked immunodeficiency, severe hypomethylation of the classical satellites 2 and 3 associated with disruption of constitutive heterochromatin, and facial anomalies. Sixty percent of ICF patients have mutations in the DNMT3B (DNA methyltransferase 3B) gene, encoding a de novo DNA methyltransferase. In the present study, we have shown that, in ICF lymphoblasts and peripheral blood, juxtacentromeric heterochromatic genes undergo dramatic changes in DNA methylation, indicating that they are bona fide targets of the DNMT3B protein. DNA methylation in heterochromatic genes dropped from about 80% in normal cells to approximately 30% in ICF cells. Hypomethylation was observed in five ICF patients and was associated with activation of these silent genes. Although DNA hypomethylation occurred in all the analyzed heterochromatic genes and in all the ICF patients, gene expression was restricted to some genes, every patient having his own group of activated genes. Histone modifications were preserved in ICF patients. Heterochromatic genes were associated with histone modifications that are typical of inactive chromatin: they had low acetylation on H3 and H4 histones and were slightly enriched in H3K9Me(3), both in ICF and controls. This was also the case for those heterochromatic genes that escaped silencing. This finding suggests that gene activation was not generalized to all the cells, but rather was restricted to a clonal cell population that may contribute to the phenotypic variability observed in ICF syndrome. A slight increase in H3K27 monomethylation was observed both in heterochromatin and active euchromatin in ICF patients; however, no correlation between this modification and activation of heterochromatic genes was found.

MethVisual - visualization and exploratory statistical analysis of DNA methylation profiles from bisulfite sequencing

Background

Exploration of DNA methylation and its impact on various regulatory mechanisms has become a very active field of research. Simultaneously there is an arising need for tools to process and analyse the data together with statistical investigation and visualisation.

Findings

MethVisual is a new application that enables exploratory analysis and intuitive visualization of DNA methylation data as is typically generated by bisulfite sequencing. The package allows the import of DNA methylation sequences, aligns them and performs quality control comparison. It comprises basic analysis steps as lollipop visualization, co-occurrence display of methylation of neighbouring and distant CpG sites, summary statistics on methylation status, clustering and correspondence analysis. The package has been developed for methylation data but can be also used for other data types for which binary coding can be inferred. The application of the package, as well as a comparison to existing DNA methylation analysis tools and its workflow based on two datasets is presented in this paper.

Conclusions

The R package MethVisual offers various analysis procedures for data that can be binarized, in particular for bisulfite sequenced methylation data. R/Bioconductor has become one of the most important environments for statistical analysis of various types of biological and medical data. Therefore, any data analysis within R that allows the integration of various data types as provided from different technological platforms is convenient. It is the first and so far the only specific package for DNA methylation analysis, in particular for bisulfite sequenced data available in R/Bioconductor enviroment. The package is available for free at http://methvisual.molgen.mpg.de/ and from the Bioconductor Consortium http://www.bioconductor.org.

Chronic high-fat diet in fathers programs β-cell dysfunction in female rat offspring

The global prevalence of obesity is increasing across most ages in both sexes. This is contributing to the early emergence of type 2 diabetes and its related epidemic. Having either parent obese is an independent risk factor for childhood obesity. Although the detrimental impacts of diet-induced maternal obesity on adiposity and metabolism in offspring are well established, the extent of any contribution of obese fathers is unclear, particularly the role of non-genetic factors in the causal pathway. Here we show that paternal high-fat-diet (HFD) exposure programs β-cell 'dysfunction' in rat F(1) female offspring. Chronic HFD consumption in Sprague-Dawley fathers induced increased body weight, adiposity, impaired glucose tolerance and insulin sensitivity. Relative to controls, their female offspring had an early onset of impaired insulin secretion and glucose tolerance that worsened with time, and normal adiposity. Paternal HFD altered the expression of 642 pancreatic islet genes in adult female offspring (P < 0.01); genes belonged to 13 functional clusters, including cation and ATP binding, cytoskeleton and intracellular transport. Broader pathway analysis of 2,492 genes differentially expressed (P < 0.05) demonstrated involvement of calcium-, MAPK- and Wnt-signalling pathways, apoptosis and the cell cycle. Hypomethylation of the Il13ra2 gene, which showed the highest fold difference in expression (1.76-fold increase), was demonstrated. This is the first report in mammals of non-genetic, intergenerational transmission of metabolic sequelae of a HFD from father to offspring.

MethylViewer: computational analysis and editing for bisulfite sequencing and methyltransferase accessibility protocol for individual templates (MAPit) projects

Bisulfite sequencing is a widely-used technique for examining cytosine DNA methylation at nucleotide resolution along single DNA strands. Probing with cytosine DNA methyltransferases followed by bisulfite sequencing (MAPit) is an effective technique for mapping protein-DNA interactions. Here, MAPit methylation footprinting with M.CviPI, a GC methyltransferase we previously cloned and characterized, was used to probe hMLH1 chromatin in HCT116 and RKO colorectal cancer cells. Because M.CviPI-probed samples contain both CG and GC methylation, we developed a versatile, visually-intuitive program, called MethylViewer, for evaluating the bisulfite sequencing results. Uniquely, MethylViewer can simultaneously query cytosine methylation status in bisulfite-converted sequences at as many as four different user-defined motifs, e.g. CG, GC, etc., including motifs with degenerate bases. Data can also be exported for statistical analysis and as publication-quality images. Analysis of hMLH1 MAPit data with MethylViewer showed that endogenous CG methylation and accessible GC sites were both mapped on single molecules at high resolution. Disruption of positioned nucleosomes on single molecules of the PHO5 promoter was detected in budding yeast using M.CviPII, increasing the number of enzymes available for probing protein-DNA interactions. MethylViewer provides an integrated solution for primer design and rapid, accurate and detailed analysis of bisulfite sequencing or MAPit datasets from virtually any biological or biochemical system.

Tissue culture-induced novel epialleles of a Myb transcription factor encoded by pericarp color1 in maize

Plants regenerated from tissue culture often display somaclonal variation, that is, somatic and often meiotically heritable phenotypic variation that can result from both genetic and epigenetic modifications. To better understand the molecular basis of somaclonal variation, we have characterized four unique tissue culture-derived epialleles of the pericarp color1 (p1) gene of maize (Zea mays L.). The progenitor p1 allele, P1-wr, is composed of multiple head-to-tail tandemly arranged copies of the complete gene unit and specifies brick-red phlobaphene pigmentation in the cob glumes. The novel epialleles identified in progeny plants regenerated from tissue culture showed partial to complete loss of p1 function indicated by pink or colorless cob glumes. Loss of pigmentation was correlated with nearly complete loss of p1 steady-state transcripts. DNA gel-blot analysis and genomic bisulfite sequencing showed that silencing of the epialleles was associated with hypermethylation of a region in the second intron of P1-wr. Presence of Unstable factor for orange1 (Ufo1), an unlinked epigenetic modifier of p1, restored the cob glume pigmentation in the silenced alleles, and such reactivation was accompanied by hypomethylation of the p1 sequence. This observation confirmed that silencing of the epialleles is indeed due to epigenetic modifications and that the p1 epialleles were capable of functioning in the presence of the correct trans-acting factors. While the low-copy regions of the genome generally undergo hypomethylation during tissue culture, our study shows that the tandemly repeated genes are also prone to hypermethylation and epigenetic silencing.

Analysis of bisulfite sequencing data from plant DNA using CyMATE

Amplifying and sequencing DNA after bisulfite treatment of genomic DNA reveals the methylation state of cytosine residues at the highest resolution possible. However, a thorough analysis is required for statistical evaluation of methylation at all sites in each genomic region. Several software tools were developed to assist in quantitative evaluation of bisulfite sequencing data from complex methylation patterns occurring in plants. This chapter describes the application of Cytosine Methylation Analysis Tool for Everyone (CyMATE). From aligned sequences, CyMATE quantifies and illustrates general and pattern-specific methylation at CG, CHG, and CHH (H = A, C, or T) sites, both per sequence and per position. CyMATE is also able to perform a quality control of sequences and to detect redundancy among individual clones. The software is able to reveal methylation patterns on complementary strands by handling data from hairpin bisulfite sequencing. The tool is freely available for non-commercial use at http://www.cymate.org .

Analysis of DNA methylation in plants by bisulfite sequencing

Methylation of cytosines is a very important epigenetic modification of genomic DNA in many different eukaryotes, and it is frequently involved in transcriptional regulation of genes. In plants, DNA methylation is regulated by a complex interplay between several methylating and demethylating enzymes. Analysis of the resulting cytosine methylation patterns with the highest resolution is achieved after sodium bisulfite treatment, deaminating nonmethylated cytosines to uracil. Subsequent PCR and sequence analysis of individual amplicons displays the degree, position, and sequence context of methylation of every cytosine residue in individual genomic sequences. We describe the application of bisulfite sequencing for the analysis of DNA methylation at defined individual sequences of plant genomic DNA.

Non-CpG methylation of the PGC-1alpha promoter through DNMT3B controls mitochondrial density

Epigenetic modification through DNA methylation is implicated in metabolic disease. Using whole-genome promoter methylation analysis of skeletal muscle from normal glucose-tolerant and type 2 diabetic subjects, we identified cytosine hypermethylation of peroxisome proliferator-activated receptor gamma (PPARgamma) coactivator-1 alpha (PGC-1alpha) in diabetic subjects. Methylation levels were negatively correlated with PGC-1alpha mRNA and mitochondrial DNA (mtDNA). Bisulfite sequencing revealed that the highest proportion of cytosine methylation within PGC-1alpha was found within non-CpG nucleotides. Non-CpG methylation was acutely increased in human myotubes by exposure to tumor necrosis factor-alpha (TNF-alpha) or free fatty acids, but not insulin or glucose. Selective silencing of the DNA methyltransferase 3B (DNMT3B), but not DNMT1 or DNMT3A, prevented palmitate-induced non-CpG methylation of PGC-1alpha and decreased mtDNA and PGC-1alpha mRNA. We provide evidence for PGC-1alpha hypermethylation, concomitant with reduced mitochondrial content in type 2 diabetic patients, and link DNMT3B to the acute fatty-acid-induced non-CpG methylation of PGC-1alpha promoter.

Trans-generation inheritance of methylation patterns in a tobacco transgene following a post-transcriptional silencing event

We have studied the inheritance of the epigenetic state of tobacco transgenes whose expression was post-transcriptionally silenced by an invertedly repeated silencer locus. We show that, in hybrids, the coding region of the target neomycin phosphotransferase (nptII) gene was almost exclusively methylated at CG configurations, and dense non-CG methylation occurred in the 3' untranslated region. Homologous sequences in the silencer locus were heavily methylated at both CG and non-CG motifs. After segregation of the silencer locus, the CG methylation but not the non-CG methylation of the target genes was transmitted to the progeny. In the segregants, we observed an overall increase of CG methylation in the target genes, associated with a re-distribution from the 3' end of the coding region towards the middle. This pattern was inherited with some fluctuation for at least two additional generations in the absence of a detectable T-DNA-derived small RNA fraction. Thus CG methylation is not cleared during meiosis and may be inherited over generations without RNA signals being present. These epi-allelic variants re-expressed the reporter gene immediately after segregation of the trigger, showing that relatively dense CG methylation (approximately 60-80%) imprinted on most of the coding region (>500 bp) did not reduce expression compared with the parental non-methylated locus. We propose that the genic CG methylation seen in euchromatic regions of the genome may originate from ancient post-transcriptional gene silencing events as a result of adventitiously produced methylation-directing RNA molecules.

QUMA: quantification tool for methylation analysis

Bisulfite sequencing, a standard method for DNA methylation profile analysis, is widely used in basic and clinical studies. This method is limited, however, by the time-consuming data analysis processes required to obtain accurate DNA methylation profiles from the raw sequence output of the DNA sequencer, and by the fact that quality checking of the results can be influenced by a researcher's bias. We have developed an interactive and easy-to-use web-based tool, QUMA (quantification tool for methylation analysis), for the bisulfite sequencing analysis of CpG methylation. QUMA includes most of the data-processing functions necessary for the analysis of bisulfite sequences. It also provides a platform for consistent quality control of the analysis. The QUMA web server is available at http://quma.cdb.riken.jp/.

Genomic imprinting of IGF2 in marsupials is methylation dependent

Background-

Parent-specific methylation of specific CpG residues is critical to imprinting in eutherian mammals, but its importance to imprinting in marsupials and, thus, the evolutionary origins of the imprinting mechanism have been the subject of controversy. This has been particularly true for the imprinted Insulin-like Growth Factor II (IGF2), a key regulator of embryonic growth in vertebrates and a focal point of the selective forces leading to genomic imprinting. The presence of the essential imprinting effector, DNMT3L, in marsupial genomes and the demonstration of a differentially methylated region (DMR) in the retrotransposon-derived imprinted gene, PEG10, in tammar wallaby argue for a role for methylation in imprinting, but several studies have found no evidence of parent-specific methylation at other imprinted loci in marsupials.

Results-

We performed the most extensive search to date for allele-specific patterns of CpG methylation within CpG isochores or CpG enriched segments across a 22 kilobase region surrounding the IGF2 gene in the South American opossum Monodelphis domestica. We identified a previously unknown 5'-untranslated exon for opossum IGF2, which is flanked by sequences defining a putative neonatal promoter, a DMR and an active Matrix Attachment Region (MAR). Demethylation of this DMR in opossum neonatal fibroblasts results in abherrant biallelic expression of IGF2.

Conclusion-

The demonstration of a DMR and an active MAR in the 5' flank of opossum IGF2 mirrors the regulatory features of the 5' flank of Igf2 in mice. However, demethylation induced activation of the maternal allele of IGF2 in opossum differs from the demethylation induced repression of the paternal Igf2 allele in mice. While it can now be concluded that parent-specific DNA methylation is an epigentic mark common to Marsupialia and Eutheria, the molecular mechanisms of transcriptional silencing at imprinted loci have clearly evolved along independent trajectories.

CpG PatternFinder: a Windows-based utility program for easy and rapid identification of the CpG methylation status of DNA

The bisulfite genomic sequencing technique is one of the most widely used techniques to study sequence-specific DNA methylation because of its unambiguous ability to reveal DNA methylation status to the order of a single nucleotide. One characteristic feature of the bisulfite genomic sequencing technique is that a number of sample sequence files will be produced from a single DNA sample. The PCR products of bisulfite-treated DNA samples cannot be sequenced directly because they are heterogeneous in nature; therefore they should be cloned into suitable plasmids and then sequenced. This procedure generates an enormous number of sample DNA sequence files as well as adding extra bases belonging to the plasmids to the sequence, which will cause problems in the final sequence comparison. Finding the methylation status for each CpG in each sample sequence is not an easy job. As a result CpG PatternFinder was developed for this purpose. The main functions of the CpG PatternFinder are: (i) to analyze the reference sequence to obtain CpG and non-CpG-C residue position information. (ii) To tailor sample sequence files (delete insertions and mark deletions from the sample sequence files) based on a configuration of ClustalW multiple alignment. (iii) To align sample sequence files with a reference file to obtain bisulfite conversion efficiency and CpG methylation status. And, (iv) to produce graphics, highlighted aligned sequence text and a summary report which can be easily exported to Microsoft Office suite. CpG PatternFinder is designed to operate cooperatively with BioEdit, a freeware on the internet. It can handle up to 100 files of sample DNA sequences simultaneously, and the total CpG pattern analysis process can be finished in minutes. CpG PatternFinder is an ideal software tool for DNA methylation studies to determine the differential methylation pattern in a large number of individuals in a population. Previously we developed the CpG Analyzer program; CpG PatternFinder is our further effort to create software tools for DNA methylation studies.

CyMATE: a new tool for methylation analysis of plant genomic DNA after bisulphite sequencing

Cytosine methylation is a hallmark of epigenetic information in the DNA of many fungi, vertebrates and plants. The technique of bisulphite genomic sequencing reveals the methylation state of every individual cytosine in a sequence, and thereby provides high-resolution data on epigenetic diversity; however, the manual evaluation and documentation of large amounts of data is laborious and error-prone. While some software is available for facilitating the analysis of mammalian DNA methylation, which is found nearly exclusively at CG sites, there is no software optimally suited for data from DNA with significant non-CG methylation. We describe CyMATE (Cytosine Methylation Analysis Tool for Everyone) for in silico analysis of DNA sequences after bisulphite conversion of plant DNA, in which methylation is more divergent with respect to sequence context and biological relevance. From aligned sequences, CyMATE includes and distinguishes methylation at CG, CHG and CHH (where H = A, C or T), and can extract both quantitative and qualitative data regarding general and pattern-specific methylation per sequence and per position, i.e. data for individual sites in a sequence and the epigenetic divergence within a sample. In addition, it can provide graphical output from alignments in either an overview or a 'zoom-in' view as pdf files. Detailed information, including a quality control of the sequencing data, is provided in text format. We applied CyMATE to the analysis of DNA methylation at transcriptionally silenced promoters in diploid and polyploid Arabidopsis and found significant hypermethylation, high stability of the methylated state independent of chromosome number, and non-redundant patterns of mC distribution. CyMATE is freely available for non-commercial use at http://www.gmi.oeaw.ac.at/CyMATE.

Sequence analysis and editing for bisulphite genomic sequencing projects

Bisulphite genomic sequencing is a widely used technique for detailed analysis of the methylation status of a region of DNA. It relies upon the selective deamination of unmethylated cytosine to uracil after treatment with sodium bisulphite, usually followed by PCR amplification of the chosen target region. Since this two-step procedure replaces all unmethylated cytosine bases with thymine, PCR products derived from unmethylated templates contain only three types of nucleotide, in unequal proportions. This can create a number of technical difficulties (e.g. for some base-calling methods) and impedes manual analysis of sequencing results (since the long runs of T or A residues are difficult to align visually with the parent sequence). To facilitate the detailed analysis of bisulphite PCR products (particularly using multiple cloned templates), we have developed a visually intuitive program that identifies the methylation status of CpG dinucleotides by analysis of raw sequence data files produced by MegaBace or ABI sequencers as well as Staden SCF trace files and plain text files. The program then also collates and presents data derived from independent templates (e.g. separate clones). This results in a considerable reduction in the time required for completion of a detailed genomic methylation project.

DNA methylation immediately adjacent to active histone marking does not silence transcription

Active promoters generally contain histone H3/H4 hyperacetylation and tri-methylation at H3 lysine 4, whereas repressed promoters are associated with DNA methylation. Here we show that the repressed erythroid-specific carbonic anhydrase II (CAII) promoter has active histone modifications localized around the transcription start, while high levels of CpG methylation are present directly upstream from these active marks. Despite the presence of active histone modifications, the repressed promoter requires hormone-induced activation for efficient preinitiation complex assembly. Transient and positional changes in histone H3/H4 acetylation and local changes in nucleosome density are evident during activation, but the bipartite epigenetic code is stably maintained. Our results suggest that active histone modifications may prevent spreading of CpG methylation towards the promoter and show that repressive DNA methylation immediately adjacent to a promoter does not necessarily repress transcription.

Detailed methylation analysis of CpG islands on human chromosome region 9p21

Deletion of 9p21 is the most commonly reported chromosomal abnormality in pediatric acute lymphoblastic leukemia, and published data suggest that the maternal chromosome is preferentially deleted. Preferential maternal deletion of 9p21 and reports of a differentially methylated region (DMR) and of parental effects in mice with lymphoma suggest there may be an unrecognized imprinted locus in this region. To screen for DMRs, we used the mcrBC/HpaII screening method and peripheral-blood DNA. Of 36 CpG islands within an 8.5-Mb region of 9p21, seven were identified as putative DMRs and were further analyzed by bisulfite sequencing. Neither any of the CpG islands nor a previously published putative DMR nearby showed evidence of differential parental methylation; however, the published DMR did demonstrate sequence-dependent differential methylation. Our data, which showed heterogeneous and low-level methylation of CpG islands, have obvious implications for methylation studies.

The trans-silencing capacity of invertedly repeated transgenes depends on their epigenetic state in tobacco

We studied the in trans-silencing capacities of a transgene locus that carried the neomycin phosphotransferase II reporter gene linked to the 35S promoter in an inverted repeat (IR). This transgene locus was originally posttranscriptionally silenced but switched to a transcriptionally silenced epiallele after in vitro tissue culture. Here, we show that both epialleles were strongly methylated in the coding region and IR center. However, by genomic sequencing, we found that the 1.0 kb region around the transcription start site was heavily methylated in symmetrical and non-symmetrical contexts in transcriptionally but not in posttranscriptionally silenced epilallele. Also, the posttranscriptionally silenced epiallele could trans-silence and trans-methylate homologous transgene loci irrespective of their genomic organization. We demonstrate that this in trans-silencing was accompanied by the production of small RNA molecules. On the other hand, the transcriptionally silenced variant could neither trans-silence nor trans-methylate homologous sequences, even after being in the same genetic background for generations and meiotic cycles. Interestingly, 5-aza-2-deoxy-cytidine-induced hypomethylation could partially restore signaling from the transcriptionally silenced epiallele. These results are consistent with the hypothesis that non-transcribed highly methylated IRs are poor silencers of homologous loci at non-allelic positions even across two generations and that transcription of the inverted sequences is essential for their trans-silencing potential.

Stat1 expression is not sufficient to regulate the interferon signaling pathway in cellular immortalization

DNA hypermethylation in gene promoters is an epigenetic mechanism regulating gene expression in cellular immortalization, an important step in carcinogenesis. Previously, we studied the genes dysregulated during immortalization using spontaneously immortalized fibroblasts from patients with Li-Fraumeni syndrome (LFS), who carry a germline mutation in the tumor suppressor gene p53. We found that multiple interferon (IFN) signaling pathway genes were regulated by epigenetic silencing. In this study we focused on a key regulator of that pathway, the signal transducer and transcription activator 1 (Stat1) gene. Although Stat1 is downregulated after cellular immortalization and upregulated in immortal MDAH041 cells after 5-aza-2'-deoxycytidine (5-aza-dC) treatment, we detected no methylation of the Stat1 promoter region in these cells before or after immortalization. To analyze the function of Stat1 in immortalization, we expressed Stat1 in immortal MDAH041 cells by stable infection, expecting to induce IFN-regulated genes or cellular senescence or both. However, the overexpression of Stat1 alone was not sufficient to repress the proliferation rate of immortal MDAH041 cells or induce senescence in immortal MDAH041 cells. We concluded that factor(s) additional to Stat1 (whether IFN dependent or not) are required for the immortalization of LFS fibroblasts.

DNA methylation in the CTCF-binding site I and the expression pattern of the H19 gene: does positive expression predict poor prognosis in early stage head and neck carcinomas?

Loss of allele-specific expression by the imprinted genes IGF2 and H19 has been correlated with a differentially methylated region (DMR) upstream to the H19 gene. The H19-DMR contains seven potential CCCTC-binding factor (CTCF) binding sites. CTCF is a chromatin insulator and a multifunctional transcription factor whose binding to the H19-DMR is suppressed by DNA methylation. Our study included a group of 41 head and neck squamous cell carcinoma (HNSCC) samples. The imprinting status of the H19 gene was analyzed in 11 out of 35 positive cases for H19 gene expression, and only 1 of them showed loss of imprinting. We detected a significant correlation (P = 0.041, Fisher's exact test) between H19 expression and tumor recurrence. Among H19 positive cases, six were T2, in which five developed recurrence and/or metastasis. Inversely, in the group of tumors that showed no H19 gene expression, 5 out of 24 were T2 and only 1 presented regional recurrence. These data support the hypothesis that H19 expression could be used as a prognostic marker to indicate recurrence in early stage tumors. We also examined the methylation of the CTCF binding site 1 in a subgroup of these samples. The H19 gene silencing and loss of imprinting were not correlated with the methylation pattern of the CTCF binding site 1. However, the significant correlation between H19 expression and tumor recurrence suggest that this transcript could be a marker for the progression of HNSCC.