Epigenomic analysis detects widespread gene-body DNA hypomethylation in chronic lymphocytic leukemia

DNA methylation and childhood asthma in the inner city

BACKGROUND

Epigenetic marks are heritable, influenced by the environment, direct the maturation of T lymphocytes, and in mice enhance the development of allergic airway disease. Thus it is important to define epigenetic alterations in asthmatic populations.

OBJECTIVE

We hypothesize that epigenetic alterations in circulating PBMCs are associated with allergic asthma.

METHODS

We compared DNA methylation patterns and gene expression in inner-city children with persistent atopic asthma versus healthy control subjects by using DNA and RNA from PBMCs. Results were validated in an independent population of asthmatic patients.

RESULTS

Comparing asthmatic patients (n = 97) with control subjects (n = 97), we identified 81 regions that were differentially methylated. Several immune genes were hypomethylated in asthma, including IL13, RUNX3, and specific genes relevant to T lymphocytes (TIGIT). Among asthmatic patients, 11 differentially methylated regions were associated with higher serum IgE concentrations, and 16 were associated with percent predicted FEV1. Hypomethylated and hypermethylated regions were associated with increased and decreased gene expression, respectively (P < 6 × 10(-12) for asthma and P < .01 for IgE). We further explored the relationship between DNA methylation and gene expression using an integrative analysis and identified additional candidates relevant to asthma (IL4 and ST2). Methylation marks involved in T-cell maturation (RUNX3), TH2 immunity (IL4), and oxidative stress (catalase) were validated in an independent asthmatic cohort of children living in the inner city.

CONCLUSIONS

Our results demonstrate that DNA methylation marks in specific gene loci are associated with asthma and suggest that epigenetic changes might play a role in establishing the immune phenotype associated with asthma.

Epigenetic remodeling in B-cell acute lymphoblastic leukemia occurs in two tracks and employs embryonic stem cell-like signatures

We investigated DNA methylomes of pediatric B-cell acute lymphoblastic leukemias (B-ALLs) using whole-genome bisulfite sequencing and high-definition microarrays, along with RNA expression profiles. Epigenetic alteration of B-ALLs occurred in two tracks: de novo methylation of small functional compartments and demethylation of large inter-compartmental backbones. The deviations were exaggerated in lamina-associated domains, with differences corresponding to methylation clusters and/or cytogenetic groups. Our data also suggested a pivotal role of polycomb and CTBP2 in de novo methylation, which may be traced back to bivalency status of embryonic stem cells. Driven by these potent epigenetic modulations, suppression of polycomb target genes was observed along with disruption of developmental fate and cell cycle and mismatch repair pathways and altered activities of key upstream regulators.

DNA methylation of tumor-suppressor miRNA genes in chronic lymphocytic leukemia

DNA methylation is one of the most important epigenetic modifications of the genome involved in the regulation of numerous cellular processes through gene silencing without altering DNA sequences. miRNAs, a class of single-stranded noncoding RNAs of 19-25 nucleotides in length, function as post-transcriptional regulators of gene expression leading to mRNA cleavage or translational repression of their corresponding target protein-coding genes. Recently, dysregulation of tumor suppressor miRNAs mediated by promoter DNA hypermethylation is implicated in human cancers, including B-cell chronic lymphocytic leukemia (CLL). Moreover, it appears that methylated miRNA genes could be potential biomarkers for CLL diagnosis or therapy. This review will highlight the role of aberrant methylation of miRNA genes in the pathogenesis of CLL.

Whole-epigenome analysis in multiple myeloma reveals DNA hypermethylation of B cell-specific enhancers

While analyzing the DNA methylome of multiple myeloma (MM), a plasma cell neoplasm, by whole-genome bisulfite sequencing and high-density arrays, we observed a highly heterogeneous pattern globally characterized by regional DNA hypermethylation embedded in extensive hypomethylation. In contrast to the widely reported DNA hypermethylation of promoter-associated CpG islands (CGIs) in cancer, hypermethylated sites in MM, as opposed to normal plasma cells, were located outside CpG islands and were unexpectedly associated with intronic enhancer regions defined in normal B cells and plasma cells. Both RNA-seq and in vitro reporter assays indicated that enhancer hypermethylation is globally associated with down-regulation of its host genes. ChIP-seq and DNase-seq further revealed that DNA hypermethylation in these regions is related to enhancer decommissioning. Hypermethylated enhancer regions overlapped with binding sites of B cell-specific transcription factors (TFs) and the degree of enhancer methylation inversely correlated with expression levels of these TFs in MM. Furthermore, hypermethylated regions in MM were methylated in stem cells and gradually became demethylated during normal B-cell differentiation, suggesting that MM cells either reacquire epigenetic features of undifferentiated cells or maintain an epigenetic signature of a putative myeloma stem cell progenitor. Overall, we have identified DNA hypermethylation of developmentally regulated enhancers as a new type of epigenetic modification associated with the pathogenesis of MM.

Molecular pathogenesis of CLL and its evolution

In spite of being the most prevalent adult leukemia in Western countries, the molecular mechanisms driving the establishment and progression of chronic lymphocytic leukemia (CLL) remain largely unknown. In recent years, the use of next-generation sequencing techniques has uncovered new and, in some cases, unexpected driver genes with prognostic and therapeutic value. The mutational landscape of CLL is characterized by high-genetic and epigenetic heterogeneity, low mutation recurrence and a long tail of cases with undefined driver genes. On the other hand, the use of deep sequencing has also revealed high intra-tumor heterogeneity and provided a detailed picture of clonal evolution processes. This phenomenon, in which aberrant DNA methylation can also participate, appears to be tightly associated to poor outcomes and chemo-refractoriness, thus providing a new subject for therapeutic intervention. Hence, and having in mind the limitations derived from the CLL complexity thus described, the application of massively parallel sequencing studies has unveiled a wealth of information that is expected to substantially improve patient staging schemes and CLL clinical management.

Tissue-specific effects of genetic and epigenetic variation on gene regulation and splicing

Understanding how genetic variation affects distinct cellular phenotypes, such as gene expression levels, alternative splicing and DNA methylation levels, is essential for better understanding of complex diseases and traits. Furthermore, how inter-individual variation of DNA methylation is associated to gene expression is just starting to be studied. In this study, we use the GenCord cohort of 204 newborn Europeans’ lymphoblastoid cell lines, T-cells and fibroblasts derived from umbilical cords. The samples were previously genotyped for 2.5 million SNPs, mRNA-sequenced, and assayed for methylation levels in 482,421 CpG sites. We observe that methylation sites associated to expression levels are enriched in enhancers, gene bodies and CpG island shores. We show that while the correlation between DNA methylation and gene expression can be positive or negative, it is very consistent across cell-types. However, this epigenetic association to gene expression appears more tissue-specific than the genetic effects on gene expression or DNA methylation (observed in both sharing estimations based on P-values and effect size correlations between cell-types). This predominance of genetic effects can also be reflected by the observation that allele specific expression differences between individuals dominate over tissue-specific effects. Additionally, we discover genetic effects on alternative splicing and interestingly, a large amount of DNA methylation correlating to alternative splicing, both in a tissue-specific manner. The locations of the SNPs and methylation sites involved in these associations highlight the participation of promoter proximal and distant regulatory regions on alternative splicing. Overall, our results provide high-resolution analyses showing how genome sequence variation has a broad effect on cellular phenotypes across cell-types, whereas epigenetic factors provide a secondary layer of variation that is more tissue-specific. Furthermore, the details of how this tissue-specificity may vary across inter-relations of molecular traits, and where these are occurring, can yield further insights into gene regulation and cellular biology as a whole.

In order to better understand how genetic differences between individuals can cause diseases, it is crucial to understand how genetic variants affect cellular functions in the different tissues that compose the human body. From the umbilical cord of 195 newborn babies, we previously obtained three different cell-types: fibroblasts, T-cells and immortalized B-cells. From every individual in each cell type we measured four features across the genome: 1) genetic differences, 2) DNA methylation, an epigenetic modification of DNA that can affect its functional state, 3) gene expression—the amount of gene activity, 4) alternative splicing—which of the different versions of a gene is manifested. We find thousands of genetic variants of the DNA sequence that affect methylation, gene expression, and splicing. We show that while these genetic effects often affect multiple cell-types, the strength of these effects varies between cell-types. Also epigenetic methylation marks of DNA associate to gene expression and particularly often to splicing. Since abnormalities in gene expression, DNA methylation and alternative splicing are associated to diseases, it is important to continue studying how these traits are inter-related and affected by genetic variation across cell-types.

Higher gene expression variability in the more aggressive subtype of chronic lymphocytic leukemia

Background

Chronic lymphocytic leukemia (CLL) presents two subtypes which have drastically different clinical outcomes, IgVH mutated (M-CLL) and IgVH unmutated (U-CLL). So far, these two subtypes are not associated to clear differences in gene expression profiles. Interestingly, recent results have highlighted important roles for heterogeneity, both at the genetic and at the epigenetic level in CLL progression.

Methods

We analyzed gene expression data of two large cohorts of CLL patients and quantified expression variability across individuals to investigate differences between the two subtypes using different measures and statistical tests. Functional significance was explored by pathway enrichment and network analyses. Furthermore, we implemented a random forest approach based on expression variability to classify patients into disease subtypes.

Results

We found that U-CLL, the more aggressive type of the disease, shows significantly increased variability of gene expression across patients and that, overall, genes that show higher variability in the aggressive subtype are related to cell cycle, development and inter-cellular communication. These functions indicate a potential relation between gene expression variability and the faster progression of this CLL subtype. Finally, a classifier based on gene expression variability was able to correctly predict the disease subtype of CLL patients.

Conclusions

There are strong relations between gene expression variability and disease subtype linking significantly increased expression variability to phenotypes such as aggressiveness and resistance to therapy in CLL.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-014-0125-z) contains supplementary material, which is available to authorized users.

Prenatal nutritional deficiency reprogrammed postnatal gene expression in mammal brains: implications for schizophrenia

BACKGROUND

Epidemiological studies have identified prenatal exposure to famine as a risk factor for schizophrenia, and animal models of prenatal malnutrition display structural and functional brain abnormalities implicated in schizophrenia.

METHODS

The offspring of the RLP50 rat, a recently developed animal model of prenatal famine malnutrition exposure, was used to investigate the changes of gene expression and epigenetic modifications in the brain regions. Microarray gene expression analysis was carried out in the prefrontal cortex and the hippocampus from 8 RLP50 offspring rats and 8 controls. MBD-seq was used to test the changes in DNA methylation in hippocampus depending on prenatal malnutrition exposure.

RESULTS

In the prefrontal cortex, offspring of RLP50 exhibit differences in neurotransmitters and olfactory-associated gene expression. In the hippocampus, the differentially-expressed genes are related to synaptic function and transcription regulation. DNA methylome profiling of the hippocampus also shows widespread but systematic epigenetic changes; in most cases (87%) this involves hypermethylation. Remarkably, genes encoded for the plasma membrane are significantly enriched for changes in both gene expression and DNA methylome profiling screens (p = 2.37×10(-9) and 5.36×10(-9), respectively). Interestingly, Mecp2 and Slc2a1, two genes associated with cognitive impairment, show significant down-regulation, and Slc2a1 is hypermethylated in the hippocampus of the RLP50 offspring.

CONCLUSIONS

Collectively, our results indicate that prenatal exposure to malnutrition leads to the reprogramming of postnatal brain gene expression and that the epigenetic modifications contribute to the reprogramming. The process may impair learning and memory ability and result in higher susceptibility to schizophrenia.

Comparison of DNA methylation profiles in human fetal and adult red blood cell progenitors

Background

DNA methylation is an epigenetic modification that plays an important role during mammalian development. Around birth in humans, the main site of red blood cell production moves from the fetal liver to the bone marrow. DNA methylation changes at the β-globin locus and a switch from fetal to adult hemoglobin production characterize this transition. Understanding this globin switch may improve the treatment of patients with sickle cell disease and β-thalassemia, two of the most common Mendelian diseases in the world. The goal of our study was to describe and compare the genome-wide patterns of DNA methylation in fetal and adult human erythroblasts.

Methods

We used the Illumina HumanMethylation 450 k BeadChip to measure DNA methylation at 402,819 CpGs in ex vivo-differentiated erythroblasts from 12 fetal liver and 12 bone marrow CD34+ donors.

Results

We identified 5,937 differentially methylated CpGs that overlap with erythroid enhancers and binding sites for erythropoiesis-related transcription factors. Combining this information with genome-wide association study results, we show that erythroid enhancers define particularly promising genomic regions to identify new genetic variants associated with fetal hemoglobin (HbF) levels in humans. Many differentially methylated CpGs are located near genes with unanticipated roles in red blood cell differentiation and proliferation. For some of these new candidate genes, we confirm the correlation between DNA methylation and gene expression levels in red blood cell progenitors. We also provide evidence that DNA methylation and genetic variation at the β-globin locus independently control globin gene expression in adult erythroblasts.

Conclusions

Our DNA methylome maps confirm the widespread dynamic changes in DNA methylation that occur during human erythropoiesis. These changes tend to happen near erythroid enhancers, further highlighting their importance in erythroid regulation and HbF production. Finally, DNA methylation may act independently of the transcription factor BCL11A to repress fetal hemoglobin production. This provides cues on strategies to more efficiently re-activate HbF production in sickle cell disease and β-thalassemia patients.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-014-0122-2) contains supplementary material, which is available to authorized users.

Molecular basis of chronic lymphocytic leukemia diagnosis and prognosis

BACKGROUNDS

Chronic lymphocytic leukemia (CLL) is the most common type of leukemia in adults and is characterized by a clonal accumulation of mature apoptosis-resistant neoplastic cells. It is also a heterogeneous disease with a variable clinical outcome. Here, we present a review of currently known (epi)genetic alterations that are related to the etiology, progression and chemo-refractoriness of CLL. Relevant literature was identified through a PubMed search (1994-2014) of English-language papers using the terms CLL, signaling pathway, cytogenetic abnormality, somatic mutation, epigenetic alteration and micro-RNA.

RESULTS

CLL is characterized by the presence of gross chromosomal abnormalities, epigenetic alterations, micro-RNA expression alterations, immunoglobulin heavy chain gene mutations and other genetic lesions. The expression of unmutated immunoglobulin heavy chain variable region (IGHV) genes, ZAP-70 and CD38 proteins, the occurrence of chromosomal abnormalities such as 17p and 11q deletions and mutations of the NOTCH1, SF3B1 and BIRC3 genes have been associated with a poor prognosis. In addition, mutations in tumor suppressor genes, such as TP53 and ATM, have been associated with refractoriness to conventional chemotherapeutic agents. Micro-RNA expression alterations and aberrant methylation patterns in genes that are specifically deregulated in CLL, including the BCL-2, TCL1 and ZAP-70 genes, have also been encountered and linked to distinct clinical parameters.

CONCLUSIONS

Specific chromosomal abnormalities and gene mutations may serve as diagnostic and prognostic indicators for disease progression and survival. The identification of these anomalies by state-of-the-art molecular (cyto)genetic techniques such as fluorescence in situ hybridization (FISH), comparative genomic hybridization (CGH), single nucleotide polymorphism (SNP) microarray-based genomic profiling and next-generation sequencing (NGS) can be of paramount help for the clinical management of these patients, including optimal treatment design. The efficacy of novel therapeutics should to be tested according to the presence of these molecular lesions in CLL patients.

DNA methylation: old dog, new tricks?

DNA methylation is an epigenetic modification that is generally associated with repression of transcription initiation at CpG-island promoters. Here we argue that, on the basis of recent high-throughput genomic and proteomic screenings, DNA methylation can also have different outcomes, including activation of transcription. This is evidenced by the fact that transcription factors can interact with methylated DNA sequences. Furthermore, in certain cellular contexts, genes containing methylated promoters are highly transcribed. Interestingly, this uncoupling between methylated DNA and repression of transcription seems to be particularly evident in germ cells and pluripotent cells. Thus, contrary to previous assumptions, DNA methylation is not exclusively associated with repression of transcription initiation.

Epigenetics of idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a complex lung disease of unknown etiology. Development of IPF is influenced by both genetic and environmental factors. Recent work by our and other groups has identified strong genetic predisposition factors for the development of pulmonary fibrosis, and cigarette smoke remains the most strongly associated environmental exposure risk factor. Gene expression profiling studies of IPF lung have taught us quite a bit about the biology of this fatal disease, and those of peripheral blood have provided important biomarkers. However, epigenetic marks may be the missing link that connects the environmental exposure in genetically predisposed individuals to transcriptional changes associated with disease development. Moreover, epigenetic marks represent a promising therapeutic target for IPF. In this review, the disease is introduced, genetic and gene expression studies in IPF are summarized, exposures relevant to IPF and known epigenetic changes associated with cigarette smoke exposure are discussed, and epigenetic studies conducted so far in IPF are summarized. Limitations, challenges, and future opportunities in this field are also discussed.

Epigenetic function of activation-induced cytidine deaminase and its link to lymphomagenesis

Activation-induced cytidine deaminase (AID) is essential for somatic hypermutation and class switch recombination of immunoglobulin (Ig) genes during B cell maturation and immune response. Expression of AID is tightly regulated due to its mutagenic and recombinogenic potential, which is known to target not only Ig genes, but also non-Ig genes, contributing to lymphomagenesis. In recent years, a new epigenetic function of AID and its link to DNA demethylation came to light in several developmental systems. In this review, we summarize existing evidence linking deamination of unmodified and modified cytidine by AID to base-excision repair and mismatch repair machinery resulting in passive or active removal of DNA methylation mark, with the focus on B cell biology. We also discuss potential contribution of AID-dependent DNA hypomethylation to lymphomagenesis.

The genomic landscape of chronic lymphocytic leukaemia: biological and clinical implications

Chronic lymphocytic leukaemia (CLL) remains at the forefront of the genetic analysis of human tumours, principally due its prevalence, protracted natural history and accessibility to suitable material for analysis. With the application of high-throughput genetic technologies, we have an unbridled view of the architecture of the CLL genome, including a comprehensive description of the copy number and mutational landscape of the disease, a detailed picture of clonal evolution during pathogenesis, and the molecular mechanisms that drive genomic instability and therapeutic resistance. This work has nuanced the prognostic importance of established copy number alterations, and identified novel prognostically relevant gene mutations that function within biological pathways that are attractive treatment targets. Herein, an overview of recent genomic discoveries will be reviewed, with associated biological and clinical implications, and a view into how clinical implementation may be facilitated.

Locally disordered methylation forms the basis of intratumor methylome variation in chronic lymphocytic leukemia

Intratumoral heterogeneity plays a critical role in tumor evolution. To define the contribution of DNA methylation to heterogeneity within tumors, we performed genome-scale bisulfite sequencing of 104 primary chronic lymphocytic leukemias (CLLs). Compared with 26 normal B cell samples, CLLs consistently displayed higher intrasample variability of DNA methylation patterns across the genome, which appears to arise from stochastically disordered methylation in malignant cells. Transcriptome analysis of bulk and single CLL cells revealed that methylation disorder was linked to low-level expression. Disordered methylation was further associated with adverse clinical outcome. We therefore propose that disordered methylation plays a similar role to that of genetic instability, enhancing the ability of cancer cells to search for superior evolutionary trajectories.

DNA hydroxymethylation profiling reveals that WT1 mutations result in loss of TET2 function in acute myeloid leukemia

Somatic mutations in IDH1/IDH2 and TET2 result in impaired TET2-mediated conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). The observation that WT1 inactivating mutations anticorrelate with TET2/IDH1/IDH2 mutations in acute myeloid leukemia (AML) led us to hypothesize that WT1 mutations may impact TET2 function. WT1 mutant AML patients have reduced 5hmC levels similar to TET2/IDH1/IDH2 mutant AML. These mutations are characterized by convergent, site-specific alterations in DNA hydroxymethylation, which drive differential gene expression more than alterations in DNA promoter methylation. WT1 overexpression increases global levels of 5hmC, and WT1 silencing reduced 5hmC levels. WT1 physically interacts with TET2 and TET3, and WT1 loss of function results in a similar hematopoietic differentiation phenotype as observed with TET2 deficiency. These data provide a role for WT1 in regulating DNA hydroxymethylation and suggest that TET2 IDH1/IDH2 and WT1 mutations define an AML subtype defined by dysregulated DNA hydroxymethylation.

Biological and Clinical Implications of Clonal Heterogeneity and Clonal Evolution in Multiple Myeloma

Clonal heterogeneity and clonal evolution have emerged as critical concepts in the field of oncology over the past four decades, largely thanks to the implementation of novel technologies such as comparative genomic hybridization, whole genome/exome sequencing and epigenetic analysis. Along with the identification of cancer stem cells in the majority of neoplasia, the recognition of intertumor and intratumor variability has provided a novel perspective to understand the mechanisms behind tumor evolution and its implication in terms of treatment failure and cancer relapse or recurrence. First hypothesized over two decades ago, clonal heterogeneity and clonal evolution have been confirmed in multiple myeloma (MM), an incurable cancer of plasma cells, almost universally preceded by a pre-malignant conditioned named monoclonal gammopathy of undetermined significance (MGUS). The genetic events and molecular mechanisms underlying such evolution have been difficult to dissect. Moreover, while a role for the bone marrow microenvironment in supporting MM cell survival, proliferation and drug-resistance has been well established, whether it is directly involved in driving evolution from MGUS to MM is at present unclear. We present in this review a historical excursus on the concepts of clonal heterogeneity and clonal evolution in MM with a special emphasis on their role in the progression from MGUS to MM; the contribution of the microenvironment; and the clinical implications in terms of resistance to treatment and disease relapse/recurrence.

Gene body methylation can alter gene expression and is a therapeutic target in cancer

DNA methylation in promoters is well known to silence genes and is the presumed therapeutic target of methylation inhibitors. Gene body methylation is positively correlated with expression, yet its function is unknown. We show that 5-aza-2'-deoxycytidine treatment not only reactivates genes but decreases the overexpression of genes, many of which are involved in metabolic processes regulated by c-MYC. Downregulation is caused by DNA demethylation of the gene bodies and restoration of high levels of expression requires remethylation by DNMT3B. Gene body methylation may, therefore, be an unexpected therapeutic target for DNA methylation inhibitors, resulting in the normalization of gene overexpression induced during carcinogenesis. Our results provide direct evidence for a causal relationship between gene body methylation and transcription.

Divergence of gene body DNA methylation and evolution of plant duplicate genes

It has been shown that gene body DNA methylation is associated with gene expression. However, whether and how deviation of gene body DNA methylation between duplicate genes can influence their divergence remains largely unexplored. Here, we aim to elucidate the potential role of gene body DNA methylation in the fate of duplicate genes. We identified paralogous gene pairs from Arabidopsis and rice (Oryza sativa ssp. japonica) genomes and reprocessed their single-base resolution methylome data. We show that methylation in paralogous genes nonlinearly correlates with several gene properties including exon number/gene length, expression level and mutation rate. Further, we demonstrated that divergence of methylation level and pattern in paralogs indeed positively correlate with their sequence and expression divergences. This result held even after controlling for other confounding factors known to influence the divergence of paralogs. We observed that methylation level divergence might be more relevant to the expression divergence of paralogs than methylation pattern divergence. Finally, we explored the mechanisms that might give rise to the divergence of gene body methylation in paralogs. We found that exonic methylation divergence more closely correlates with expression divergence than intronic methylation divergence. We show that genomic environments (e.g., flanked by transposable elements and repetitive sequences) of paralogs generated by various duplication mechanisms are associated with the methylation divergence of paralogs. Overall, our results suggest that the changes in gene body DNA methylation could provide another avenue for duplicate genes to develop differential expression patterns and undergo different evolutionary fates in plant genomes.

Epigenetic diversity in hematopoietic neoplasms

Tumor cell populations display a remarkable extent of variability in non-genetic characteristics such as DNA methylation, histone modification patterns, and differentiation levels of individual cells. It remains to be elucidated whether non-genetic heterogeneity is simply a byproduct of tumor evolution or instead a manifestation of a higher-order tissue organization that is maintained within the neoplasm to establish a differentiation hierarchy, a favorable microenvironment, or a buffer against changing selection pressures during tumorigenesis. Here, we review recent findings on epigenetic diversity, particularly heterogeneity in DNA methylation patterns in hematologic malignancies. We also address the implications of epigenetic heterogeneity for the clonal evolution of tumors and discuss its effects on gene expression and other genome functions in cancer.

The CD27+ memory B cells display changes in the gene expression pattern in elderly individuals

Memory B cells (MBCs) have a very long life-span as compared to naïve B cells (NBCs), remaining viable for years. It could predispose them to suffer misbalances in the gene expression pattern at the long term, which might be involved in the development of age-related B-cell disorders. In order to identify genes whose expression might change during life, we analyzed the gene expression patterns of CD27- NBCs versus CD27+ MBCs in young and old subjects. Using microarray assays we observed that the expression pattern of CD27- NBCs versus CD27+ MBCs is significantly different. Furthermore, in order to evaluate the age effect, we compared the gene expression pattern of young versus aged subjects in both cell populations. Interestingly, we did not find significant differences in the CD27- NBC population between young and aged individuals, whereas we found 925 genes differentially expressed in CD27+ MBCs. Among these genes, 193 were also differentially expressed in CD27+ MBCs as compared to CD27- NBCs, most of them involved in cell survival, cell growth and proliferation, cellular development and gene expression. We conclude that gene expression profiles of CD27- NBCs and CD27+ MBCs are different. Moreover, whereas the gene expression pattern of CD27+ MBCs varies with age, the same does not happen in CD27- NBCs. This suggests that MBCs undergo time-dependent changes which could underlie a higher susceptibility to dysfunction with age. This article is protected by copyright. All rights reserved.

Pan-cancer patterns of DNA methylation

The comparison of DNA methylation patterns across cancer types (pan-cancer methylome analyses) has revealed distinct subgroups of tumors that share similar methylation patterns. Integration of these data with the wealth of information derived from cancer genome profiling studies performed by large international consortia has provided novel insights into the cellular aberrations that contribute to cancer development. There is evidence that genetic mutations in epigenetic regulators (such as DNMT3, IDH1/2 or H3.3) mediate or contribute to these patterns, although a unifying molecular mechanism underlying the global alterations of DNA methylation has largely been elusive. Knowledge gained from pan-cancer methylome analyses will aid the development of diagnostic and prognostic biomarkers, improve patient stratification and the discovery of novel druggable targets for therapy, and will generate hypotheses for innovative clinical trial designs based on methylation subgroups rather than on cancer subtypes. In this review, we discuss recent advances in the global profiling of tumor genomes for aberrant DNA methylation and the integration of these data with cancer genome profiling data, highlight potential mechanisms leading to different methylation subgroups, and show how this information can be used in basic research and for translational applications. A remaining challenge is to experimentally prove the functional link between observed pan-cancer methylation patterns, the associated genetic aberrations, and their relevance for the development of cancer.

Genome-wide DNA methylation profiles in progression to in situ and invasive carcinoma of the breast with impact on gene transcription and prognosis

BACKGROUND

Ductal carcinoma in situ (DCIS) of the breast is a precursor of invasive breast carcinoma. DNA methylation alterations are thought to be an early event in progression of cancer, and may prove valuable as a tool in clinical decision making and for understanding neoplastic development.

RESULTS

We generate genome-wide DNA methylation profiles of 285 breast tissue samples representing progression of cancer, and validate methylation changes between normal and DCIS in an independent dataset of 15 normal and 40 DCIS samples. We also validate a prognostic signature on 583 breast cancer samples from The Cancer Genome Atlas. Our analysis reveals that DNA methylation profiles of DCIS are radically altered compared to normal breast tissue, involving more than 5,000 genes. Changes between DCIS and invasive breast carcinoma involve around 1,000 genes. In tumors, DNA methylation is associated with gene expression of almost 3,000 genes, including both negative and positive correlations. A prognostic signature based on methylation level of 18 CpGs is associated with survival of breast cancer patients with invasive tumors, as well as with survival of patients with DCIS and mixed lesions of DCIS and invasive breast carcinoma.

CONCLUSIONS

This work demonstrates that changes in the epigenome occur early in the neoplastic progression, provides evidence for the possible utilization of DNA methylation-based markers of progression in the clinic, and highlights the importance of epigenetic changes in carcinogenesis.

Genome-wide DNA methylation profiles in progression to in situ and invasive carcinoma of the breast with impact on gene transcription and prognosis

BACKGROUND

Ductal carcinoma in situ (DCIS) of the breast is a precursor of invasive breast carcinoma. DNA methylation alterations are thought to be an early event in progression of cancer, and may prove valuable as a tool in clinical decision making and for understanding neoplastic development.

RESULTS

We generate genome-wide DNA methylation profiles of 285 breast tissue samples representing progression of cancer, and validate methylation changes between normal and DCIS in an independent dataset of 15 normal and 40 DCIS samples. We also validate a prognostic signature on 583 breast cancer samples from The Cancer Genome Atlas. Our analysis reveals that DNA methylation profiles of DCIS are radically altered compared to normal breast tissue, involving more than 5,000 genes. Changes between DCIS and invasive breast carcinoma involve around 1,000 genes. In tumors, DNA methylation is associated with gene expression of almost 3,000 genes, including both negative and positive correlations. A prognostic signature based on methylation level of 18 CpGs is associated with survival of breast cancer patients with invasive tumors, as well as with survival of patients with DCIS and mixed lesions of DCIS and invasive breast carcinoma.

CONCLUSIONS

This work demonstrates that changes in the epigenome occur early in the neoplastic progression, provides evidence for the possible utilization of DNA methylation-based markers of progression in the clinic, and highlights the importance of epigenetic changes in carcinogenesis.

A comprehensive DNA methylation profile of epithelial-to-mesenchymal transition

Epithelial-to-mesenchymal transition (EMT) is a plastic process in which fully differentiated epithelial cells are converted into poorly differentiated, migratory and invasive mesenchymal cells, and it has been related to the metastasis potential of tumors. This is a reversible process and cells can also eventually undergo mesenchymal-to-epithelial transition. The existence of a dynamic EMT process suggests the involvement of epigenetic shifts in the phenotype. Herein, we obtained the DNA methylomes at single-base resolution of Madin-Darby canine kidney cells undergoing EMT and translated the identified differentially methylated regions to human breast cancer cells undergoing a gain of migratory and invasive capabilities associated with the EMT phenotype. We noticed dynamic and reversible changes of DNA methylation, both on promoter sequences and gene-bodies in association with transcription regulation of EMT-related genes. Most importantly, the identified DNA methylation markers of EMT were present in primary mammary tumors in association with the epithelial or the mesenchymal phenotype of the studied breast cancer samples.

p19-INK4d inhibits neuroblastoma cell growth, induces differentiation and is hypermethylated and downregulated in MYCN-amplified neuroblastomas

Uncontrolled cell cycle entry, resulting from deregulated CDK-RB1-E2F pathway activity, is a crucial determinant of neuroblastoma cell malignancy. Here we identify neuroblastoma-suppressive functions of the p19-INK4d CDK inhibitor and uncover mechanisms of its repression in high-risk neuroblastomas. Reduced p19-INK4d expression was associated with poor event-free and overall survival and neuroblastoma risk factors including amplified MYCN in a set of 478 primary neuroblastomas. High MYCN expression repressed p19-INK4d mRNA and protein levels in different neuroblastoma cell models with conditional MYCN expression. MassARRAY and 450K methylation analyses of 105 primary neuroblastomas uncovered a differentially methylated region within p19-INK4d. Hypermethylation of this region was associated with reduced p19-INK4d expression. In accordance, p19-INK4d expression was activated upon treatment with the demethylating agent, 2'-deoxy-5-azacytidine, in neuroblastoma cell lines. Ectopic p19-INK4d expression decreased viability, clonogenicity and the capacity for anchorage-independent growth of neuroblastoma cells, and shifted the cell cycle towards the G1/0 phase. p19-INK4d also induced neurite-like processes and markers of neuronal differentiation. Moreover, neuroblastoma cell differentiation, induced by all-trans retinoic acid or NGF-NTRK1-signaling, activated p19-INK4d expression. Our findings pinpoint p19-INK4d as a neuroblastoma suppressor and provide evidence for MYCN-mediated repression and for epigenetic silencing of p19-INK4d by DNA hypermethylation in high-risk neuroblastomas.

Chronic lymphocytic leukaemia: could immunological tolerance mechanisms be the origin of lymphoid neoplasms?

Immunological tolerance theory in chronic lymphocytic leukaemia (CLL): we suggest that B cells that express B-cell receptors (BCR) that recognize their own BCR epitopes are viewed by immune system as 'dangerous cells'. BCR autonomous signalling may induce constant receptor editing and mistakes in allelic exclusion. The fact that whole BCR recognizes a self-antigen or foreing antigen may be irrelevant in early B cell development. In early B cells, autonomous signalling induced by recognition of the BCR's own epitopes simulates an antigen-antibody engagement. In the bone marrow this interaction is viewed as recognition of self-molecules and induces receptor editing. In mature B cells autonomous signalling by the BCR may promote 'reversible anergy' and also may correct self-reactivity induced by the somatic hypermutation mechanisms in mutated CLL B cells. However, in unmutated CLL B cells, BCR autonomous signalling in addition to self-antigen recognition augments B cell activation, proliferation and genomic instability. We suggest that CLL originates from a coordinated normal immunologic tolerance mechanism to destroy self-reactive B cells. Additional genetic damage induced by tolerance mechanisms may immortalize self-reactive B cells and transform them into a leukemia.

The tissue is the issue: improved methylome analysis from paraffin-embedded tissues by application of the HOPE technique

Alterations in the DNA methylome are characteristic for numerous diseases and a typical hallmark of cancer. Therefore, DNA methylation is currently under investigation in research labs and has also entered diagnostics. Recently, protocols like the BeadChip technology have become commercially available to study DNA methylation in an array format and semiquantitative fashion. However, it is known that fixation of the sample material with formalin prior to BeadChip analysis can affect the results. In this study we compared the influence of fixation on the outcome of BeadChip analysis. From six patients each a lung cancer tissue sample and a corresponding tumor-free lung tissue sample were collected. The samples were separated into three pieces. One piece of each sample was fixed with formalin, another one by the non-cross-linking HOPE technique (Hepes-glutamic acid buffer mediated Organic solvent Protection Effect). Subsequently, both became paraffin embedded. As a reference, the remaining third piece was cryopreserved. In addition we used three adenocarcinoma cell lines (H838, A549, and H1650) to validate the results from patient tissues. We show that using the HOPE technique instead of formalin largely prevents the introduction of formalin-fixation related artifacts. An ANOVA analysis significantly separated HOPE- and cryopreserved from formalin-fixed samples (FDR<0.05), while differences in the methylation data obtained from HOPE-fixed and cryopreserved material were minor. Consequently, HOPE fixation is superior to formalin fixation if a subsequent BeadChip analysis of paraffin-embedded sample material is intended.

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.

Evolving understanding of the CLL genome

Over the past few years, massively parallel sequencing technologies have revealed with high resolution the tremendous genetic and epigenetic heterogeneity in chronic lymphocytic leukemia (CLL). We have learned how the molecular architecture differs not only between affected individuals but also within samples and over time. These insights have catalyzed our understanding of the pathobiology of CLL and point to critical signaling pathways in the development and progression of the disease. Several key driver alterations have been identified, which serve to refine prognostic schemata but also to inspire the development of new therapeutic strategies. Ongoing advances in technology promise to further elucidate the molecular basis of CLL, and this knowledge is anticipated to aid us in understanding and addressing the clinical challenge presented by the vast variability in the clinical course of patients with CLL.

Expression of myeloperoxidase in acute myeloid leukemia blasts mirrors the distinct DNA methylation pattern involving the downregulation of DNA methyltransferase DNMT3B

Myeloperoxidase (MPO) has been associated with both a myeloid lineage commitment and favorable prognosis in patients with acute myeloid leukemia (AML). DNA methyltransferase inhibitors (decitabine and zeburaline) induced MPO gene promoter demethylation and MPO gene transcription in AML cells with low MPO activity. Therefore, MPO gene transcription was directly and indirectly regulated by DNA methylation. A DNA methylation microarray subsequently revealed a distinct methylation pattern in 33 genes, including DNA methyltransferase 3 beta (DNMT3B), in CD34-positive cells obtained from AML patients with a high percentage of MPO-positive blasts. Based on the inverse relationship between the methylation status of DNMT3B and MPO, we found an inverse relationship between DNMT3B and MPO transcription levels in CD34-positive AML cells (P=0.0283). In addition, a distinct methylation pattern was observed in five genes related to myeloid differentiation or therapeutic sensitivity in CD34-positive cells from AML patients with a high percentage of MPO-positive blasts. Taken together, the results of the present study indicate that MPO may serve as an informative marker for identifying a distinct and crucial DNA methylation profile in CD34-positive AML cells.

Epigenetically regulated microRNAs and their prospect in cancer diagnosis

Epigenetic alterations have been reported to deregulate the expression of many transcripts, including noncoding RNAs that have no apparent protein-coding capacity. Recently, as the result of numerous studies focused on miRNAs, novel sequencing technologies have made available the transcription profile of the entire human genome. miRNAs as drivers of tumor-suppressive and oncogenic functions have been found to be dysregulated in numerous cancer types. However, the functions of epigenetically regulated genetic elements other than protein-coding genes are still a matter of debate. In this review, the authors focus mainly on describing the epigenetic regulation of miRNAs in cancer. They also discuss the role of miRNAs as potential diagnostic and/or prognostic biomarkers.

The genome-wide molecular signature of transcription factors in leukemia

Transcription factors control expression of genes essential for the normal functioning of the hematopoietic system and regulate development of distinct blood cell types. During leukemogenesis, aberrant regulation of transcription factors such as RUNX1, CBFβ, MLL, C/EBPα, SPI1, GATA, and TAL1 is central to the disease. Here, we will discuss the mechanisms of transcription factor deregulation in leukemia and how in recent years next-generation sequencing approaches have helped to elucidate the molecular role of many of these aberrantly expressed transcription factors. We will focus on the complexes in which these factors reside, the role of posttranslational modification of these factors, their involvement in setting up higher order chromatin structures, and their influence on the local epigenetic environment. We suggest that only comprehensive knowledge on all these aspects will increase our understanding of aberrant gene expression in leukemia as well as open new entry points for therapeutic intervention.

A symbiotic liaison between the genetic and epigenetic code

With rapid advances in sequencing technologies, we are undergoing a paradigm shift from hypothesis- to data-driven research. Genome-wide profiling efforts have given informative insights into biological processes; however, considering the wealth of variation, the major challenge still remains in their meaningful interpretation. In particular sequence variation in non-coding contexts is often challenging to interpret. Here, data integration approaches for the identification of functional genetic variability represent a possible solution. Exemplary, functional linkage analysis integrating genotype and expression data determined regulatory quantitative trait loci and proposed causal relationships. In addition to gene expression, epigenetic regulation and specifically DNA methylation was established as highly valuable surrogate mark for functional variance of the genetic code. Epigenetic modification has served as powerful mediator trait to elucidate mechanisms forming phenotypes in health and disease. Particularly, integrative studies of genetic and DNA methylation data have been able to guide interpretation strategies of risk genotypes, but also proved their value for physiological traits, such as natural human variation and aging. This Review seeks to illustrate the power of data integration in the genomic era exemplified by DNA methylation quantitative trait loci. However, the model is further extendable to virtually all traceable molecular traits.

Recurrent epimutations activate gene body promoters in primary glioblastoma

Aberrant DNA hypomethylation may play an important role in the growth rate of glioblastoma (GBM), but the functional impact on transcription remains poorly understood. We assayed the GBM methylome with MeDIP-seq and MRE-seq, adjusting for copy number differences, in a small set of non-glioma CpG island methylator phenotype (non-G-CIMP) primary tumors. Recurrent hypomethylated loci were enriched within a region of chromosome 5p15 that is specified as a cancer amplicon and also encompasses TERT, encoding telomerase reverse transcriptase, which plays a critical role in tumorigenesis. Overall, 76 gene body promoters were recurrently hypomethylated, including TERT and the oncogenes GLI3 and TP73. Recurring hypomethylation also affected previously unannotated alternative promoters, and luciferase reporter assays for three of four of these promoters confirmed strong promoter activity in GBM cells. Histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq on tissue from the GBMs uncovered peaks that coincide precisely with tumor-specific decrease of DNA methylation at 200 loci, 133 of which are in gene bodies. Detailed investigation of TP73 and TERT gene body hypomethylation demonstrated increased expression of corresponding alternate transcripts, which in TP73 encodes a truncated p73 protein with oncogenic function and in TERT encodes a putative reverse transcriptase-null protein. Our findings suggest that recurring gene body promoter hypomethylation events, along with histone H3K4 trimethylation, alter the transcriptional landscape of GBM through the activation of a limited number of normally silenced promoters within gene bodies, in at least one case leading to expression of an oncogenic protein.

The prognostic impact of minimal residual disease in patients with chronic lymphocytic leukemia requiring first-line therapy

A proportion of patients with chronic lymphocytic leukemia achieve a minimal residual disease negative status after therapy. We retrospectively evaluated the impact of minimal residual disease on the outcome of 255 consecutive patients receiving any front-line therapy in the context of a detailed prognostic evaluation, including assessment of IGHV, TP53, NOTCH1 and SF3B1 mutations. The median follow-up was 73 months (range, 2-202) from disease evaluation. The median treatment-free survival durations for patients achieving a complete response without or with minimal residual disease, a partial response and no response were 76, 40, 11 and 11 months, respectively (P<0.001). Multivariate analysis revealed that three variables had a significant impact on treatment-free survival: minimal residual disease (P<0.001), IGHV status (P<0.001) and β2-microglobulin levels (P=0.012). With regards to overall survival, factors predictive of an unfavorable outcome were minimal residual disease positivity (P=0.014), together with advanced age (P<0.001), unmutated IGHV status (P=0.001), TP53 mutations (P<0.001) and elevated levels of β2-microglobulin (P=0.003). In conclusion, for patients requiring front-line therapy, achievement of minimal residual disease negativity is associated with significantly prolonged treatment-free and overall survival irrespective of other prognostic markers or treatment administered.

Somatic mutation as a mechanism of Wnt/β-catenin pathway activation in CLL

One major goal of cancer genome sequencing is to identify key genes and pathways that drive tumor pathogenesis. Although many studies have identified candidate driver genes based on recurrence of mutations in individual genes, subsets of genes with nonrecurrent mutations may also be defined as putative drivers if they affect a single biological pathway. In this fashion, we previously identified Wnt signaling as significantly mutated through large-scale massively parallel DNA sequencing of chronic lymphocytic leukemia (CLL). Here, we use a novel method of biomolecule delivery, vertical silicon nanowires, to efficiently introduce small interfering RNAs into CLL cells, and interrogate the effects of 8 of 15 mutated Wnt pathway members identified across 91 CLLs. In HEK293T cells, mutations in 2 genes did not generate functional changes, 3 led to dysregulated pathway activation, and 3 led to further activation or loss of repression of pathway activation. Silencing 4 of 8 mutated genes in CLL samples harboring the mutated alleles resulted in reduced viability compared with leukemia samples with wild-type alleles. We demonstrate that somatic mutations in CLL can generate dependence on this pathway for survival. These findings support the notion that nonrecurrent mutations at different nodes of the Wnt pathway can contribute to leukemogenesis.

Pathway bridge based multiobjective optimization approach for lurking pathway prediction

Ovarian carcinoma immunoreactive antigen-like protein 2 (OCIAD2) is a protein with unknown function. Frequently methylated or downregulated, OCIAD2 has been observed in kinds of tumors, and TGFβ signaling has been proved to induce the expression of OCIAD2. However, current pathway analysis tools do not cover the genes without reported interactions like OCIAD2 and also miss some significant genes with relatively lower expression. To investigate potential biological milieu of OCIAD2, especially in cancer microenvironment, a nova approach pbMOO was created to find the potential pathways from TGFβ to OCIAD2 by searching on the pathway bridge, which consisted of cancer enriched looping patterns from the complicated entire protein interactions network. The pbMOO approach was further applied to study the modulator of ligand TGFβ1, receptor TGFβR1, intermediate transfer proteins, transcription factor, and signature OCIAD2. Verified by literature and public database, the pathway TGFβ1-TGFβR1-SMAD2/3-SMAD4/AR-OCIAD2 was detected, which concealed the androgen receptor (AR) which was the possible transcription factor of OCIAD2 in TGFβsignal, and it well explained the mechanism of TGFβ induced OCIAD2 expression in cancer microenvironment, therefore providing an important clue for the future functional analysis of OCIAD2 in tumor pathogenesis.

Next-generation sequencing in chronic lymphocytic leukemia

The use of next-generation sequencing (NGS) has expanded our knowledge of the genomic alterations in chronic lymphocytic leukemia (CLL) and provides new tools for analyzing leukemic clonal architecture. Recent studies have demonstrated substantial differences in genomic alterations between mutated and unmutated IGHV subgroups, which reflect distinct molecular pathways and mutagenic mechanisms in the pathogenesis of the disease. The mutational profile of CLL can be characterized by a relatively low number of somatic mutations per case, few recurrent mutations at moderate frequency (5%-15%) and a long tail of recurrent lower frequency somatic mutations. Functional and clinical studies of novel mutations have uncovered new mechanisms involved in the pathogenesis of the disease, revealing new insights into CLL molecular evolution that could ultimately translate into improvements in the management of patients. The clonal architecture of CLL shows striking heterogeneity between patients, which could have important clinical implications. In summary, NGS studies of CLL are expanding our fundamental knowledge on the molecular mechanisms involved in the pathogenesis of the disease and offering new perspectives for the clinical management of the patients.

TET proteins and 5-methylcytosine oxidation in the immune system

DNA methylation in the form of 5-methylcytosine (5mC) is essential for normal development in mammals and influences a variety of biological processes, including transcriptional regulation, imprinting, and the maintenance of genomic stability. The recent discovery of TET proteins, which oxidize 5mC to 5-hydroxymethylcytosine, 5-formylcytosine, and 5-carboxylcytosine, has changed our understanding of the process of DNA demethylation. Here, we summarize our current knowledge of the roles of DNA methylation and TET proteins in cell differentiation and function. The intensive research on this subject has so far focused primarily on embryonic stem (ES) cells and neurons. In addition, we summarize what is known about DNA methylation in T-cell function.

Evolution of DNA methylation is linked to genetic aberrations in chronic lymphocytic leukemia

Although clonal selection by genetic driver aberrations in cancer is well documented, the ability of epigenetic alterations to promote tumor evolution is undefined. We used 450k arrays and next-generation sequencing to evaluate intratumor heterogeneity and evolution of DNA methylation and genetic aberrations in chronic lymphocytic leukemia (CLL). CLL cases exhibit vast interpatient differences in intratumor methylation heterogeneity, with genetically clonal cases maintaining low methylation heterogeneity and up to 10% of total CpGs in a monoallelically methylated state. Increasing methylation heterogeneity correlates with advanced genetic subclonal complexity. Selection of novel DNA methylation patterns is observed only in cases that undergo genetic evolution, and independent genetic evolution is uncommon and is restricted to low-risk alterations. These results reveal that although evolution of DNA methylation occurs in high-risk, clinically progressive cases, positive selection of novel methylation patterns entails coevolution of genetic alteration(s) in CLL.

The relationship between DNA methylation, genetic and expression inter-individual variation in untransformed human fibroblasts

BACKGROUND

DNA methylation plays an essential role in the regulation of gene expression. While its presence near the transcription start site of a gene has been associated with reduced expression, the variation in methylation levels across individuals, its environmental or genetic causes, and its association with gene expression remain poorly understood.

RESULTS

We report the joint analysis of sequence variants, gene expression and DNA methylation in primary fibroblast samples derived from a set of 62 unrelated individuals. Approximately 2% of the most variable CpG sites are mappable in cis to sequence variation, usually within 5 kb. Via eQTL analysis with microarray data combined with mapping of allelic expression regions, we obtained a set of 2,770 regions mappable in cis to sequence variation. In 9.5% of these expressed regions, an associated SNP was also a methylation QTL. Methylation and gene expression are often correlated without direct discernible involvement of sequence variation, but not always in the expected direction of negative for promoter CpGs and positive for gene-body CpGs. Population-level correlation between methylation and expression is strongest in a subset of developmentally significant genes, including all four HOX clusters. The presence and sign of this correlation are best predicted using specific chromatin marks rather than position of the CpG site with respect to the gene.

CONCLUSIONS

Our results indicate a wide variety of relationships between gene expression, DNA methylation and sequence variation in untransformed adult human fibroblasts, with considerable involvement of chromatin features and some discernible involvement of sequence variation.

Bisulfighter: accurate detection of methylated cytosines and differentially methylated regions

Analysis of bisulfite sequencing data usually requires two tasks: to call methylated cytosines (mCs) in a sample, and to detect differentially methylated regions (DMRs) between paired samples. Although numerous tools have been proposed for mC calling, methods for DMR detection have been largely limited. Here, we present Bisulfighter, a new software package for detecting mCs and DMRs from bisulfite sequencing data. Bisulfighter combines the LAST alignment tool for mC calling, and a novel framework for DMR detection based on hidden Markov models (HMMs). Unlike previous attempts that depend on empirical parameters, Bisulfighter can use the expectation-maximization algorithm for HMMs to adjust parameters for each data set. We conduct extensive experiments in which accuracy of mC calling and DMR detection is evaluated on simulated data with various mC contexts, read qualities, sequencing depths and DMR lengths, as well as on real data from a wide range of biological processes. We demonstrate that Bisulfighter consistently achieves better accuracy than other published tools, providing greater sensitivity for mCs with fewer false positives, more precise estimates of mC levels, more exact locations of DMRs and better agreement of DMRs with gene expression and DNase I hypersensitivity. The source code is available at http://epigenome.cbrc.jp/bisulfighter.