Genome-wide signatures of differential DNA methylation in pediatric acute lymphoblastic leukemia

Background

Although aberrant DNA methylation has been observed previously in acute lymphoblastic leukemia (ALL), the patterns of differential methylation have not been comprehensively determined in all subtypes of ALL on a genome-wide scale. The relationship between DNA methylation, cytogenetic background, drug resistance and relapse in ALL is poorly understood.

Results

We surveyed the DNA methylation levels of 435,941 CpG sites in samples from 764 children at diagnosis of ALL and from 27 children at relapse. This survey uncovered four characteristic methylation signatures. First, compared with control blood cells, the methylomes of ALL cells shared 9,406 predominantly hypermethylated CpG sites, independent of cytogenetic background. Second, each cytogenetic subtype of ALL displayed a unique set of hyper- and hypomethylated CpG sites. The CpG sites that constituted these two signatures differed in their functional genomic enrichment to regions with marks of active or repressed chromatin. Third, we identified subtype-specific differential methylation in promoter and enhancer regions that were strongly correlated with gene expression. Fourth, a set of 6,612 CpG sites was predominantly hypermethylated in ALL cells at relapse, compared with matched samples at diagnosis. Analysis of relapse-free survival identified CpG sites with subtype-specific differential methylation that divided the patients into different risk groups, depending on their methylation status.

Conclusions

Our results suggest an important biological role for DNA methylation in the differences between ALL subtypes and in their clinical outcome after treatment.     

Discovery and expansion of gene modules by seeking isolated groups in a random graph process

Background

A central problem in systems biology research is the identification and extension of biological modules–groups of genes or proteins participating in a common cellular process or physical complex. As a result, there is a persistent need for practical, principled methods to infer the modular organization of genes from genome-scale data.

Results

We introduce a novel approach for the identification of modules based on the persistence of isolated gene groups within an evolving graph process. First, the underlying genomic data is summarized in the form of ranked gene–gene relationships, thereby accommodating studies that quantify the relevant biological relationship directly or indirectly. Then, the observed gene–gene relationship ranks are viewed as the outcome of a random graph process and candidate modules are given by the identifiable subgraphs that arise during this process. An isolation index is computed for each module, which quantifies the statistical significance of its survival time.

Conclusions

The Miso (module isolation) method predicts gene modules from genomic data and the associated isolation index provides a module-specific measure of confidence. Improving on existing alternative, such as graph clustering and the global pruning of dendrograms, this index offers two intuitively appealing features: (1) the score is module-specific; and (2) different choices of threshold correlate logically with the resulting performance, i.e. a stringent cutoff yields high quality predictions, but low sensitivity. Through the analysis of yeast phenotype data, the Miso method is shown to outperform existing alternatives, in terms of the specificity and sensitivity of its predictions.     

An atypical epigenetic mechanism affects uniparental expression of Pol IV-dependent siRNAs

Background

Small RNAs generated by RNA polymerase IV (Pol IV) are the most abundant class of small RNAs in flowering plants. In Arabidopsis thaliana Pol IV-dependent short interfering (p4-si)RNAs are imprinted and accumulate specifically from maternal chromosomes in the developing seeds. Imprinted expression of protein-coding genes is controlled by differential DNA or histone methylation placed in gametes. To identify epigenetic factors required for maternal-specific expression of p4-siRNAs we analyzed the effect of a series of candidate mutations, including those required for genomic imprinting of protein-coding genes, on uniparental expression of a representative p4-siRNA locus.

Results

Paternal alleles of imprinted genes are marked by DNA or histone methylation placed by DNA METHYLTRANSFERASE 1 or the Polycomb Repressive Complex 2. Here we demonstrate that repression of paternal p4-siRNA expression at locus 08002 is not controlled by either of these mechanisms. Similarly, loss of several chromatin modification enzymes, including a histone acetyltransferase, a histone methyltransferase, and two nucleosome remodeling proteins, does not affect maternal expression of locus 08002. Maternal alleles of imprinted genes are hypomethylated by DEMETER DNA glycosylase, yet expression of p4-siRNAs occurs irrespective of demethylation by DEMETER or related glycosylases.

Conclusions

Differential DNA methylation and other chromatin modifications associated with epigenetic silencing are not required for maternal-specific expression of p4-siRNAs at locus 08002. These data indicate that there is an as yet unknown epigenetic mechanism causing maternal-specific p4-siRNA expression that is distinct from the well-characterized mechanisms associated with DNA methylation or the Polycomb Repressive Complex 2.     

The diversity of small non-coding RNAs in the diatom Phaeodactylum tricornutum

Background

Marine diatoms constitute a major component of eukaryotic phytoplankton and stand at the crossroads of several evolutionary lineages. These microalgae possess peculiar genomic features and novel combinations of genes acquired from bacterial, animal and plant ancestors. Furthermore, they display both DNA methylation and gene silencing activities. Yet, the biogenesis and regulatory function of small RNAs (sRNAs) remain ill defined in diatoms.

Results

Here we report the first comprehensive characterization of the sRNA landscape and its correlation with genomic and epigenomic information in Phaeodactylum tricornutum. The majority of sRNAs is 25 to 30 nt-long and maps to repetitive and silenced Transposable Elements marked by DNA methylation. A subset of this population also targets DNA methylated protein-coding genes, suggesting that gene body methylation might be sRNA-driven in diatoms. Remarkably, 25-30 nt sRNAs display a well-defined and unprecedented 180 nt-long periodic distribution at several highly methylated regions that awaits characterization. While canonical miRNAs are not detectable, other 21-25 nt sRNAs of unknown origin are highly expressed. Besides, non-coding RNAs with well-described function, namely tRNAs and U2 snRNA, constitute a major source of 21-25 nt sRNAs and likely play important roles under stressful environmental conditions.

Conclusions

P. tricornutum has evolved diversified sRNA pathways, likely implicated in the regulation of largely still uncharacterized genetic and epigenetic processes. These results uncover an unexpected complexity of diatom sRNA population and previously unappreciated features, providing new insights into the diversification of sRNA-based processes in eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-698) contains supplementary material, which is available to authorized users.     

Map of open and closed chromatin domains in Drosophila genome

Background

Chromatin compactness has been considered a major determinant of gene activity and has been associated with specific chromatin modifications in studies on a few individual genetic loci. At the same time, genome-wide patterns of open and closed chromatin have been understudied, and are at present largely predicted from chromatin modification and gene expression data. However the universal applicability of such predictions is not self-evident, and requires experimental verification.

Results

We developed and implemented a high-throughput analysis for general chromatin sensitivity to DNase I which provides a comprehensive epigenomic assessment in a single assay. Contiguous domains of open and closed chromatin were identified by computational analysis of the data, and correlated to other genome annotations including predicted chromatin “states”, individual chromatin modifications, nuclear lamina interactions, and gene expression. While showing that the widely trusted predictions of chromatin structure are correct in the majority of cases, we detected diverse “exceptions” from the conventional rules. We found a profound paucity of chromatin modifications in a major fraction of closed chromatin, and identified a number of loci where chromatin configuration is opposite to that expected from modification and gene expression patterns. Further, we observed that chromatin of large introns tends to be closed even when the genes are expressed, and that a significant proportion of active genes including their promoters are located in closed chromatin.

Conclusions

These findings reveal limitations of the existing predictive models, indicate novel mechanisms of epigenetic regulation, and provide important insights into genome organization and function.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-988) contains supplementary material, which is available to authorized users.     

Role of CpG deserts in the epigenetic transgenerational inheritance of differential DNA methylation regions

Background

Previously a variety of environmental toxicants were found to promote the epigenetic transgenerational inheritance of disease through differential DNA methylation regions (DMRs), termed epimutations, present in sperm. The transgenerational epimutations in sperm and somatic cells identified in a number of previous studies were further investigated.

Results

The epimutations from six different environmental exposures were found to be predominantly exposure specific with negligible overlap. The current report describes a major genomic feature of all the unique epimutations identified (535) as a very low (<10 CpG/100 bp) CpG density in sperm and somatic cells associated with transgenerational disease. The genomic locations of these epimutations were found to contain DMRs with small clusters of CpG within a general region of very low density CpG. The potential role of these epimutations on gene expression is suggested to be important.

Conclusions

Observations suggest a potential regulatory role for lower density CpG regions termed “CpG deserts”. The potential evolutionary origins of these regions is also discussed.