Abstract Background We performed a comparative analysis of the genome-wide DNA methylation profiles from three human embryonic stem cell (HESC) lines. It had previously been shown that HESC lines had significantly higher non-CG methylation than differentiated cells, and we therefore asked whether these sites were conserved across cell lines. Results We find that heavily methylated non-CG sites are strongly conserved, especially when found within the motif TACAG. They are enriched in splice sites and are more methylated than other non-CG sites in genes. We next studied the relationship between allele-specific expression and allele-specific methylation. By combining bisulfite sequencing and whole transcriptome shotgun sequencing (RNA-seq) data we identified 1,020 genes that show allele-specific expression, and 14% of CG sites genome-wide have allele-specific methylation. Finally, we asked whether the methylation state of transcription factor binding sites affects the binding of transcription factors. We identified variations in methylation levels at binding sites and found that for several transcription factors the correlation between the methylation at binding sites and gene expression is generally stronger than in the neighboring sequences. Conclusions These results suggest a possible but as yet unknown functional role for the highly methylated conserved non-CG sites in the regulation of HESCs. We also identified a novel set of genes that are likely transcriptionally regulated by methylation in an allele-specific manner. The analysis of transcription factor binding sites suggests that the methylation state of cis-regulatory elements impacts the ability of factors to bind and regulate transcription.

Abstract Background Chemotaxis is the process by which motile bacteria sense their chemical environment and move towards more favourable conditions. Escherichia coli utilises a single sensory pathway, but little is known about signalling pathways in species with more complex systems. Results To investigate whether chemotaxis pathways in other bacteria follow the E. coli paradigm, we analysed 206 species encoding at least 1 homologue of each of the 5 core chemotaxis proteins (CheA, CheB, CheR, CheW and CheY). 61 species encode more than one of all of these 5 proteins, suggesting they have multiple chemotaxis pathways. Operon information is not available for most bacteria, so we developed a novel statistical approach to cluster che genes into putative operons. Using operon-based models, we reconstructed putative chemotaxis pathways for all 206 species. We show that cheA-cheW and cheR-cheB have strong preferences to occur in the same operon as two-gene blocks, which may reflect a functional requirement for co-transcription. However, other che genes, most notably cheY, are more dispersed on the genome. Comparison of our operons with shuffled equivalents demonstrates that specific patterns of genomic location may be a determining factor for the observed in vivo chemotaxis pathways. We then examined the chemotaxis pathways of Rhodobacter sphaeroides. Here, the PpfA protein is known to be critical for correct partitioning of proteins in the cytoplasmically-localised pathway. We found ppfA in che operons of many species, suggesting that partitioning of cytoplasmic Che protein clusters is common. We also examined the apparently non-typical chemotaxis components, CheA3, CheA4 and CheY6. We found that though variants of CheA proteins are rare, the CheY6 variant may be a common type of CheY, with a significantly disordered C-terminal region which may be functionally significant. Conclusions We find that many bacterial species potentially have multiple chemotaxis pathways, with grouping of che genes into operons likely to be a major factor in keeping signalling pathways distinct. Gene order is highly conserved with cheA-cheW and cheR-cheB blocks, perhaps reflecting functional linkage. CheY behaves differently to other Che proteins, both in its genomic location and its putative protein interactions, which should be considered when modelling chemotaxis pathways.

BACKGROUND: DNA methylation is an important epigenetic modification involved in many biological processes. Bisulfite treatment coupled with high-throughput sequencing provides an effective approach for studying genome-wide DNA methylation at base resolution. Libraries such as whole genome bisulfite sequencing (WGBS) and reduced represented bisulfite sequencing (RRBS) are widely used for generating DNA methylomes, demanding efficient and versatile tools for aligning bisulfite sequencing data. RESULTS: We have developed BS-Seeker2, an updated version of BS Seeker, as a full pipeline for mapping bisulfite sequencing data and generating DNA methylomes. BS-Seeker2 improves mappability over existing aligners by using local alignment. It can also map reads from RRBS library by building special indexes with improved efficiency and accuracy. Moreover, BS-Seeker2 provides additional function for filtering out reads with incomplete bisulfite conversion, which is useful in minimizing the overestimation of DNA methylation levels. We also defined CGmap and ATCGmap file formats for full representations of DNA methylomes, as part of the outputs of BS-Seeker2 pipeline together with BAM and WIG files. CONCLUSIONS: Our evaluations on the performance show that BS-Seeker2 works efficiently and accurately for both WGBS data and RRBS data. BS-Seeker2 is freely available at http://pellegrini.mcdb.ucla.edu/BS_Seeker2/ and the Galaxy server.

Abstract Background DNA methylation is an important epigenetic modification involved in many biological processes. Bisulfite treatment coupled with high-throughput sequencing provides an effective approach for studying genome-wide DNA methylation at base resolution. Libraries such as whole genome bisulfite sequencing (WGBS) and reduced represented bisulfite sequencing (RRBS) are widely used for generating DNA methylomes, demanding efficient and versatile tools for aligning bisulfite sequencing data. Results We have developed BS-Seeker2, an updated version of BS Seeker, as a full pipeline for mapping bisulfite sequencing data and generating DNA methylomes. BS-Seeker2 improves mappability over existing aligners by using local alignment. It can also map reads from RRBS library by building special indexes with improved efficiency and accuracy. Moreover, BS-Seeker2 provides additional function for filtering out reads with incomplete bisulfite conversion, which is useful in minimizing the overestimation of DNA methylation levels. We also defined CGmap and ATCGmap file formats for full representations of DNA methylomes, as part of the outputs of BS-Seeker2 pipeline together with BAM and WIG files. Conclusions Our evaluations on the performance show that BS-Seeker2 works efficiently and accurately for both WGBS data and RRBS data. BS-Seeker2 is freely available at http://pellegrini.mcdb.ucla.edu/BS_Seeker2/ and the Galaxy server.

Background

Tuber melanosporum, also known in the gastronomic community as "truffle", features one of the largest fungal genomes (125 Mb) with an exceptionally high transposable element (TE) and repetitive DNA content (>58%). The main purpose of DNA methylation in fungi is TE silencing. As obligate outcrossing organisms, truffles are bound to a sexual mode of propagation, which together with TEs is thought to represent a major force driving the evolution of DNA methylation. Thus, it was of interest to examine if and how T. melanosporum exploits DNA methylation to maintain genome integrity.

Findings

We performed whole-genome DNA bisulfite sequencing and mRNA sequencing on different developmental stages of T. melanosporum; namely, fruitbody ("truffle"), free-living mycelium and ectomycorrhiza. The data revealed a high rate of cytosine methylation (>44%), selectively targeting TEs rather than genes with a strong preference for CpG sites. Whole genome DNA sequencing uncovered multiple TE-enriched, copy number variant regions bearing a significant fraction of hypomethylated and expressed TEs, almost exclusively in free-living mycelium propagated in vitro. Treatment of mycelia with 5-azacytidine partially reduced DNA methylation and increased TE transcription. Our transcriptome assembly also resulted in the identification of a set of novel transcripts from 614 genes.

Conclusions

The datasets presented here provide valuable and comprehensive (epi)genomic information that can be of interest for evolutionary genomics studies of multicellular (filamentous) fungi, in particular Ascomycetes belonging to the subphylum, Pezizomycotina. Evidence derived from comparative methylome and transcriptome analyses indicates that a non-exhaustive and partly reversible methylation process operates in truffles.

Global DNA demethylation in humans is a fundamental process that occurs in pre-implantation embryos and reversion to naive ground state pluripotent stem cells (PSCs). However, the extent of DNA methylation reprogramming in human germline cells is unknown. Here, we performed whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq) of human prenatal germline cells from 53 to 137 days of development. We discovered that the transcriptome and methylome of human germline is distinct from both human PSCs and the inner cell mass (ICM) of human blastocysts. Using this resource to monitor the outcome of global DNA demethylation with reversion of primed PSCs to the naive ground state, we uncovered hotspots of ultralow methylation at transposons that are protected from demethylation in the germline and ICM. Taken together, the human germline serves as a valuable in vivo tool for monitoring the epigenome of cells that have emerged from a global DNA demethylation event.

Background

We investigated how an extremely transposon element (TE)-rich organism such as the plant-symbiotic ascomycete truffle Tuber melanosporum exploits DNA methylation to cope with the more than 45,000 repeated elements that populate its genome.

Results

Whole-genome bisulfite sequencing performed on different developmental stages reveals a high fraction of methylated cytosines with a strong preference for CpG sites. The methylation pattern is highly similar among samples and selectively targets TEs rather than genes. A marked trend toward hypomethylation is observed for TEs located within a 1 kb distance from expressed genes, rather than segregated in TE-rich regions of the genome. Approximately 300 hypomethylated or unmethylated TEs are transcriptionally active, with higher expression levels in free-living mycelium compared to fruitbody. Indeed, multiple TE-enriched, copy number variant regions bearing a significant fraction of hypomethylated and expressed TEs are found almost exclusively in free-living mycelium. A reduction of DNA methylation, restricted to non-CpG sites and accompanied by an increase in TE expression, is observed upon treatment of free-living mycelia with 5-azacytidine.

Conclusions

Evidence derived from analysis of the T. melanosporum methylome indicates that a non-exhaustive, partly reversible, methylation process operates in truffles. This allows for the existence of hypomethylated, transcriptionally active TEs that are associated with copy number variant regions of the genome. Non-exhaustive TE methylation may reflect a role of active TEs in promoting genome plasticity and the ability to adapt to sudden environmental changes.

Transition from primed to naive pluripotency is associated with dynamic changes in transposable element (TE) expression and demethylation of imprinting control regions (ICRs). In mouse, ICR methylation and TE expression are each regulated by TRIM28; however, the role of TRIM28 in humans is less clear. Here, we show that a null mutation in TRIM28 causes significant alterations in TE expression in both the naive and primed states of human pluripotency, and phenotypically this has limited effects on self-renewal, instead causing a loss of germline competency. Furthermore, we discovered that TRIM28 regulates paternal ICR methylation and chromatin accessibility in the primed state, with no effects on maternal ICRs. Taken together, our study shows that abnormal TE expression is tolerated by self-renewing human pluripotent cells, whereas germline competency is not.

Remodeling DNA methylation in mammalian genomes can be global, as seen in preimplantation embryos and primordial germ cells (PGCs), or locus specific, which can regulate neighboring gene expression. In PGCs, global and locus-specific DNA demethylation occur in sequential stages, with an initial global decrease in methylated cytosines (stage I) followed by a Tet methylcytosine dioxygenase (Tet)-dependent decrease in methylated cytosines that act at imprinting control regions (ICRs) and meiotic genes (stage II). The purpose of the two-stage mechanism is unclear. Here we show that Dnmt1 preserves DNA methylation through stage I at ICRs and meiotic gene promoters and is required for the pericentromeric enrichment of 5hmC. We discovered that the functional consequence of abrogating two-stage DNA demethylation in PGCs was precocious germline differentiation leading to hypogonadism and infertility. Therefore, bypassing stage-specific DNA demethylation has significant consequences for progenitor germ cell differentiation and the ability to transmit DNA from parent to offspring.

Inherent genetic programming and environmental factors affect fetal growth in utero. Epidemiologic data in growth-altered fetuses, either intrauterine growth restricted (IUGR) or large for gestational age (LGA), demonstrate that these newborns are at increased risk of cardiometabolic disease in adulthood. There is growing evidence that the in utero environment leads to epigenetic modification, contributing to eventual risk of developing heart disease or diabetes. In this study, we used reduced representation bisulfite sequencing to examine genome-wide DNA methylation variation in placental samples from offspring born IUGR, LGA, and appropriate for gestational age (AGA) and to identify differential methylation of genes important for conferring risk of cardiometabolic disease. We found that there were distinct methylation signatures for IUGR, LGA, and AGA groups and identified over 500 differentially methylated genes (DMGs) among these group comparisons. Functional and gene network analyses revealed expected relationships of DMGs to placental physiology and transport, but also identified novel pathways with biologic plausibility and potential clinical importance to cardiometabolic disease. Specific loci for DMGs of interest had methylation patterns that were strongly associated with anthropometric presentations. We further validated altered gene expression of these specific DMGs contributing to vascular and metabolic diseases (SLC36A1, PTPRN2, CASZ1, IL10), thereby establishing transcriptional effects toward assigning functional significance. Our results suggest that the gene expression and methylation state of the human placenta are related and sensitive to the intrauterine environment, as it affects fetal growth patterns. We speculate that these observed changes may affect risk for offspring in developing adult cardiometabolic disease.