Comparison of sequencing-based methods to profile DNA methylation and identification of monoallelic epigenetic modifications.

Analysis of DNA methylation patterns relies increasingly on sequencing-based profiling methods. The four most frequently used sequencing-based technologies are the bisulfite-based methods MethylC-seq and reduced representation bisulfite sequencing (RRBS), and the enrichment-based techniques methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA binding domain sequencing (MBD-seq). We applied all four methods to biological replicates of human embryonic stem cells to assess their genome-wide CpG coverage, resolution, cost, concordance and the influence of CpG density and genomic context. The methylation levels assessed by the two bisulfite methods were concordant (their difference did not exceed a given threshold) for 82% for CpGs and 99% of the non-CpG cytosines. Using binary methylation calls, the two enrichment methods were 99% concordant and regions assessed by all four methods were 97% concordant. We combined MeDIP-seq with methylation-sensitive restriction enzyme (MRE-seq) sequencing for comprehensive methylome coverage at lower cost. This, along with RNA-seq and ChIP-seq of the ES cells enabled us to detect regions with allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression

Features of mammalian microRNA promoters emerge from polymerase II chromatin immunoprecipitation data

Background: MicroRNAs (miRNAs) are short, non-coding RNA regulators of protein coding genes. miRNAs play a very important role in diverse biological processes and various diseases. Many algorithms are able to predict miRNA genes and their targets, but their transcription regulation is still under investigation. It is generally believed that intragenic miRNAs (located in introns or exons of protein coding genes) are co-transcribed with their host genes and most intergenic miRNAs transcribed from their own RNA polymerase II (Pol II) promoter. However, the length of the primary transcripts and promoter organization is currently unknown. Methodology: We performed Pol II chromatin immunoprecipitation (ChIP)-chip using a custom array surrounding regions of known miRNA genes. To identify the true core transcription start sites of the miRNA genes we developed a new tool (CPPP). We showed that miRNA genes can be transcribed from promoters located several kilobases away and that their promoters share the same general features as those of protein coding genes. Finally, we found evidence that as many as 26% of the intragenic miRNAs may be transcribed from their own unique promoters. Conclusion: miRNA promoters have similar features to those of protein coding genes, but miRNA transcript organization is more complex. © 2009 Corcoran et al

Identification of Type 1 Diabetes-Associated DNA Methylation Variable Positions That Precede Disease Diagnosis

Monozygotic (MZ) twin pair discordance for childhood-onset Type 1 Diabetes (T1D) is similar to 50%, implicating roles for genetic and non-genetic factors in the aetiology of this complex autoimmune disease. Although significant progress has been made in elucidating the genetics of T1D in recent years, the non-genetic component has remained poorly defined. We hypothesized that epigenetic variation could underlie some of the non-genetic component of T1D aetiology and, thus, performed an epigenome-wide association study (EWAS) for this disease. We generated genome-wide DNA methylation profiles of purified CD14(+) monocytes (an immune effector cell type relevant to T1D pathogenesis) from 15 T1D-discordant MZ twin pairs. This identified 132 different CpG sites at which the direction of the intra-MZ pair DNA methylation difference significantly correlated with the diabetic state, i.e. T1D-associated methylation variable positions (T1D-MVPs). We confirmed these T1D-MVPs display statistically significant intra-MZ pair DNA methylation differences in the expected direction in an independent set of T1D-discordant MZ pairs (P = 0.035). Then, to establish the temporal origins of the T1D-MVPs, we generated two further genome-wide datasets and established that, when compared with controls, T1D-MVPs are enriched in singletons both before (P = 0.001) and at (P = 0.015) disease diagnosis, and also in singletons positive for diabetes-associated autoantibodies but disease-free even after 12 years follow-up (P = 0.0023). Combined, these results suggest that T1D-MVPs arise very early in the etiological process that leads to overt T1D. Our EWAS of T1D represents an important contribution toward understanding the etiological role of epigenetic variation in type 1 diabetes, and it is also the first systematic analysis of the temporal origins of disease-associated epigenetic variation for any human complex disease

DNA Methylation of the First Exon Is Tightly Linked to Transcriptional Silencing

Tissue specific patterns of methylated cytosine residues vary with age, can be altered by environmental factors, and are often abnormal in human disease yet the cellular consequences of DNA methylation are incompletely understood. Although the bodies of highly expressed genes are often extensively methylated in plants, the relationship between intragenic methylation and expression is less clear in mammalian cells. We performed genome-wide analyses of DNA methylation and gene expression to determine how the pattern of intragenic methylation correlates with transcription and to assess the relationship between methylation of exonic and intronic portions of the gene body. We found that dense exonic methylation is far more common than previously recognized or expected statistically, yet first exons are relatively spared compared to more downstream exons and introns. Dense methylation surrounding the transcription start site (TSS) is uncoupled from methylation within more downstream regions suggesting that there are at least two classes of intragenic methylation. Whereas methylation surrounding the TSS is tightly linked to transcriptional silencing, methylation of more downstream regions is unassociated with the magnitude of gene expression. Notably, we found that DNA methylation downstream of the TSS, in the region of the first exon, is much more tightly linked to transcriptional silencing than is methylation in the upstream promoter region. These data provide direct evidence that DNA methylation is interpreted dissimilarly in different regions of the gene body and suggest that first exon methylation blocks transcript initiation, or vice versa. Our data also show that once initiated, downstream methylation is not a significant impediment to polymerase extension. Thus, the consequences of most intragenic DNA methylation must extend beyond the modulation of transcription magnitude

High Resolution Detection and Analysis of CpG Dinucleotides Methylation Using MBD-Seq Technology

Methyl-CpG binding domain protein sequencing (MBD-seq) is widely used to survey DNA methylation patterns. However, the optimal experimental parameters for MBD-seq remain unclear and the data analysis remains challenging. In this study, we generated high depth MBD-seq data in MCF-7 cell and developed a bi-asymmetric-Laplace model (BALM) to perform data analysis. We found that optimal efficiency of MBD-seq experiments was achieved by sequencing ∼100 million unique mapped tags from a combination of 500 mM and 1000 mM salt concentration elution in MCF-7 cells. Clonal bisulfite sequencing results showed that the methylation status of each CpG dinucleotides in the tested regions was accurately detected with high resolution using the proposed model. These results demonstrated the combination of MBD-seq and BALM could serve as a useful tool to investigate DNA methylome due to its low cost, high specificity, efficiency and resolution

High Resolution Methylome Map of Rat Indicates Role of Intragenic DNA Methylation in Identification of Coding Region

DNA methylation is crucial for gene regulation and maintenance of genomic stability. Rat has been a key model system in understanding mammalian systemic physiology, however detailed rat methylome remains uncharacterized till date. Here, we present the first high resolution methylome of rat liver generated using Methylated DNA immunoprecipitation and high throughput sequencing (MeDIP-Seq) approach. We observed that within the DNA/RNA repeat elements, simple repeats harbor the highest degree of methylation. Promoter hypomethylation and exon hypermethylation were common features in both RefSeq genes and expressed genes (as evaluated by proteomic approach). We also found that although CpG islands were generally hypomethylated, about 6% of them were methylated and a large proportion (37%) of methylated islands fell within the exons. Notably, we obeserved significant differences in methylation of terminal exons (UTRs); methylation being more pronounced in coding/partially coding exons compared to the non-coding exons. Further, events like alternate exon splicing (cassette exon) and intron retentions were marked by DNA methylation and these regions are retained in the final transcript. Thus, we suggest that DNA methylation could play a crucial role in marking coding regions thereby regulating alternative splicing. Apart from generating the first high resolution methylome map of rat liver tissue, the present study provides several critical insights into methylome organization and extends our understanding of interplay between epigenome, gene expression and genome stability

Chromosome-wide DNA methylation analysis predicts human tissue-specific X inactivation

X-chromosome inactivation (XCI) results in the differential marking of the active and inactive X with epigenetic modifications including DNA methylation. Consistent with the previous studies showing that CpG island-containing promoters of genes subject to XCI are approximately 50% methylated in females and unmethylated in males while genes which escape XCI are unmethylated in both sexes; our chromosome-wide (Methylated DNA ImmunoPrecipitation) and promoter-targeted methylation analyses (Illumina Infinium HumanMethylation27 array) showed the largest methylation difference (D = 0.12, p < 2.2 E−16) between male and female blood at X-linked CpG islands promoters. We used the methylation differences between males and females to predict XCI statuses in blood and found that 81% had the same XCI status as previously determined using expression data. Most genes (83%) showed the same XCI status across tissues (blood, fetal: muscle, kidney and nerual); however, the methylation of a subset of genes predicted different XCI statuses in different tissues. Using previously published expression data the effect of transcription on gene-body methylation was investigated and while X-linked introns of highly expressed genes were more methylated than the introns of lowly expressed genes, exonic methylation did not differ based on expression level. We conclude that the XCI status predicted using methylation of X-linked promoters with CpG islands was usually the same as determined by expression analysis and that 12% of X-linked genes examined show tissue-specific XCI whereby a gene has a different XCI status in at least one of the four tissues examined

Notes for genera: basal clades of Fungi (including Aphelidiomycota, Basidiobolomycota, Blastocladiomycota, Calcarisporiellomycota, Caulochytriomycota, Chytridiomycota, Entomophthoromycota, Glomeromycota, Kickxellomycota, Monoblepharomycota, Mortierellomycota, Mucoromycota, Neocallimastigomycota, Olpidiomycota, Rozellomycota and Zoopagomycota)

Compared to the higher fungi (Dikarya), taxonomic and evolutionary studies on the basal clades of fungi are fewer in number. Thus, the generic boundaries and higher ranks in the basal clades of fungi are poorly known. Recent DNA based taxonomic studies have provided reliable and accurate information. It is therefore necessary to compile all available information since basal clades genera lack updated checklists or outlines. Recently, Tedersoo et al. (MycoKeys 13:1--20, 2016) accepted Aphelidiomycota and Rozellomycota in Fungal clade. Thus, we regard both these phyla as members in Kingdom Fungi. We accept 16 phyla in basal clades viz. Aphelidiomycota, Basidiobolomycota, Blastocladiomycota, Calcarisporiellomycota, Caulochytriomycota, Chytridiomycota, Entomophthoromycota, Glomeromycota, Kickxellomycota, Monoblepharomycota, Mortierellomycota, Mucoromycota, Neocallimastigomycota, Olpidiomycota, Rozellomycota and Zoopagomycota. Thus, 611 genera in 153 families, 43 orders and 18 classes are provided with details of classification, synonyms, life modes, distribution, recent literature and genomic data. Moreover, Catenariaceae Couch is proposed to be conserved, Cladochytriales Mozl.-Standr. is emended and the family Nephridiophagaceae is introduced