wg-blimp: an end-to-end analysis pipeline for whole genome bisulfite sequencing data

Genome-wide DNA methylation changes in human spermatogenesis

Sperm production and function require the correct establishment of DNA methylation patterns in the germline. Here, we examined the genome-wide DNA methylation changes during human spermatogenesis and its alterations in disturbed spermatogenesis. We found that spermatogenesis is associated with remodeling of the methylome, comprising a global decline in DNA methylation in primary spermatocytes followed by selective remethylation, resulting in a spermatids/sperm-specific methylome. Hypomethylated regions in spermatids/sperm were enriched in specific transcription factor binding sites for DMRT and SOX family members and spermatid-specific genes. Intriguingly, while SINEs displayed differential methylation throughout spermatogenesis, LINEs appeared to be protected from changes in DNA methylation. In disturbed spermatogenesis, germ cells exhibited considerable DNA methylation changes, which were significantly enriched at transposable elements and genes involved in spermatogenesis. We detected hypomethylation in SVA and L1HS in disturbed spermatogenesis, suggesting an association between the abnormal programming of these regions and failure of germ cells progressing beyond meiosis.

Circulating Nucleic Acids in Colorectal Cancer: Diagnostic and Prognostic Value

Colorectal cancer (CRC) is the third most prevalent cancer in the world and the fourth leading cause of cancer-related mortality. DNA (cfDNA/ctDNA) and RNA (cfRNA/ctRNA) in the blood are promising noninvasive biomarkers for molecular profiling, screening, diagnosis, treatment management, and prognosis of CRC. Technological advancements that enable precise detection of both genetic and epigenetic abnormalities, even in minute quantities in circulation, can overcome some of these challenges. This review focuses on testing for circulating nucleic acids in the circulation as a noninvasive method for CRC detection, monitoring, detection of minimal residual disease, and patient management. In addition, the benefits and drawbacks of various diagnostic techniques and associated bioinformatics tools have been detailed.

Source cell-type epigenetic memory persists in induced pluripotent cells but is lost in subsequently derived germline cells

Introduction: Retention of source cell-type epigenetic memory may mitigate the potential for induced pluripotent stem cells (iPSCs) to fully achieve transitions in cell fate in vitro. While this may not preclude the use of iPSC-derived somatic cell types for therapeutic applications, it becomes a major concern impacting the potential use of iPSC-derived germline cell types for reproductive applications. The transition from a source somatic cell type to iPSCs and then on to germ-cell like cells (GCLCs) recapitulates two major epigenetic reprogramming events that normally occur during development in vivo-embryonic reprogramming in the epiblast and germline reprogramming in primordial germ cells (PGCs). We examined the extent of epigenetic and transcriptomic memory persisting first during the transition from differentiated source cell types to iPSCs, and then during the transition from iPSCs to PGC-like cells (PGCLCs). Methods: We derived iPSCs from four differentiated mouse cell types including two somatic and two germ cell types and tested the extent to which each resulting iPSC line resembled a) a validated ES cell reference line, and b) their respective source cell types, on the basis of genome-wide gene expression and DNA methylation patterns. We then induced each iPSC line to form PGCLCs, and assessed epigenomic and transcriptomic memory in each compared to endogenous PGCs/M-prospermatogonia. Results: In each iPSC line, we found residual gene expression and epigenetic programming patterns characteristic of the corresponding source differentiated cell type from which each was derived. However, upon deriving PGCLCs, we found very little evidence of lingering epigenetic or transcriptomic memory of the original source cell type. Discussion: This result indicates that derivation of iPSCs and then GCLCs from differentiated source cell types in vitro recapitulates the two-phase epigenetic reprogramming that normally occurs in vivo, and that, to a significant extent, germline cell types derived in vitro from pluripotent cells accurately recapitulate epigenetic programming and gene expression patterns corresponding to equivalent endogenous germ cell types, suggesting that they have the potential to form the basis of in vitro gametogenesis as a useful therapeutic strategy for treatment of infertility.

Approaches for the Analysis and Interpretation of Whole-Genome Bisulfite Sequencing Data

DNA methylation is a covalent modification of DNA that plays important roles in processes such as the regulation of gene expression, transcription factor binding, and suppression of transposable elements. The use of whole-genome bisulfite sequencing (WGBS) enables the genome-wide identification and quantification of DNA methylation patterns at single-base resolution and is the gold standard for the analysis of DNA methylation. However, the computational analysis of WGBS data can be particularly challenging, as many computationally intensive steps are required. Here, we outline step-by-step an approach for the analysis and interpretation of WGBS data. First, sequencing reads must be trimmed, quality-checked, and aligned to the genome. Second, DNA methylation levels are estimated at each cytosine position using the aligned sequence reads of the bisulfite-treated DNA. Third, regions of differential cytosine methylation between samples can be identified. Finally, these data need to be visualized and interpreted in the context of the biological question at hand.

Identification of intergenerational epigenetic inheritance by whole genome DNA methylation analysis in trios

Genome-wide association studies have identified thousands of loci associated with common diseases and traits. However, a large fraction of heritability remains unexplained. Epigenetic modifications, such as the observed in DNA methylation have been proposed as a mechanism of intergenerational inheritance. To investigate the potential contribution of DNA methylation to the missing heritability, we analysed the methylomes of four healthy trios (two parents and one offspring) using whole genome bisulphite sequencing. Of the 1.5 million CpGs (19%) with over 20% variability between parents in at least one family and compatible with a Mendelian inheritance pattern, only 3488 CpGs (0.2%) lacked correlation with any SNP in the genome, marking them as potential sites for intergenerational epigenetic inheritance. These markers were distributed genome-wide, with some preference to be located in promoters. They displayed a bimodal distribution, being either fully methylated or unmethylated, and were often found at the boundaries of genomic regions with high/low GC content. This analysis provides a starting point for future investigations into the missing heritability of simple and complex traits.

Accelerating the alignment processing speed of the comprehensive end-to-end whole-genome bisulfite sequencing pipeline, wg-blimp

Analyzing whole-genome bisulfite and related sequencing datasets is a time-intensive process due to the complexity and size of the input raw sequencing files and lengthy read alignment step requiring correction for conversion of all unmethylated Cs to Ts genome-wide. The objective of this study was to modify the read alignment algorithm associated with the whole-genome bisulfite sequencing methylation analysis pipeline (wg-blimp) to shorten the time required to complete this phase while retaining overall read alignment accuracy. Here, we report an update to the recently published pipeline wg-blimp achieved by replacing the use of the bwa-meth aligner with the faster gemBS aligner. This improvement to the wg-blimp pipeline has led to a more than ×7 acceleration in the processing speed of samples when scaled to larger publicly available FASTQ datasets containing 80-160 million reads while maintaining nearly identical accuracy of properly mapped reads when compared with data from the previous pipeline. The modifications to the wg-blimp pipeline reported here merge the speed and accuracy of the gemBS aligner with the comprehensive analysis and data visualization assets of the wg-blimp pipeline to provide a significantly accelerated workflow that can produce high-quality data much more rapidly without compromising read accuracy at the expense of increasing RAM requirements up to 48 GB.

msPIPE: a pipeline for the analysis and visualization of whole-genome bisulfite sequencing data

BACKGROUND

DNA methylation is an important epigenetic modification that is known to regulate gene expression. Whole-genome bisulfite sequencing (WGBS) is a powerful method for studying cytosine methylation in a whole genome. However, it is difficult to obtain methylation profiles using the WGBS raw reads and is necessary to be proficient in all types of bioinformatic tools for the study of DNA methylation. In addition, recent end-to-end pipelines for DNA methylation analyses are not sufficient for addressing those difficulties.

RESULTS

Here we present msPIPE, a pipeline for DNA methylation analyses with WGBS data seamlessly connecting all the required tasks ranging from data pre-processing to multiple downstream DNA methylation analyses. The msPIPE can generate various methylation profiles to analyze methylation patterns in the given sample, including statistical summaries and methylation levels. Also, the methylation levels in the functional regions of a genome are computed with proper annotation. The results of methylation profiles, hypomethylation, and differential methylation analysis are plotted in publication-quality figures. The msPIPE can be easily and conveniently used with a Docker image, which includes all dependent packages and software related to DNA methylation analyses.

CONCLUSION

msPIPE is a new end-to-end pipeline designed for methylation calling, profiling, and various types of downstream DNA methylation analyses, leading to the creation of publication-quality figures. msPIPE allows researchers to process and analyze the WGBS data in an easy and convenient way. It is available at https://github.com/jkimlab/msPIPE and https://hub.docker.com/r/jkimlab/mspipe .

Analysis and Performance Assessment of the Whole Genome Bisulfite Sequencing Data Workflow: Currently Available Tools and a Practical Guide to Advance DNA Methylation Studies

DNA methylation is associated with transcriptional repression, genomic imprinting, stem cell differentiation, embryonic development, and inflammation. Aberrant DNA methylation can indicate disease states, including cancer and neurological disorders. Therefore, the prevalence and location of 5-methylcytosine in the human genome is a topic of interest. Whole-genome bisulfite sequencing (WGBS) is a high-throughput method for analyzing DNA methylation. This technique involves library preparation, alignment, and quality control. Advancements in epigenetic technology have led to an increase in DNA methylation studies. This review compares the detailed experimental methodology of WGBS using accessible and up-to-date analysis tools. Practical codes for WGBS data processing are included as a general guide to assist progress in DNA methylation studies through a comprehensive case study.

MethylScore, a pipeline for accurate and context-aware identification of differentially methylated regions from population-scale plant whole-genome bisulfite sequencing data

Whole-genome bisulfite sequencing (WGBS) is the standard method for profiling DNA methylation at single-nucleotide resolution. Different tools have been developed to extract differentially methylated regions (DMRs), often built upon assumptions from mammalian data. Here, we present MethylScore, a pipeline to analyse WGBS data and to account for the substantially more complex and variable nature of plant DNA methylation. MethylScore uses an unsupervised machine learning approach to segment the genome by classification into states of high and low methylation. It processes data from genomic alignments to DMR output and is designed to be usable by novice and expert users alike. We show how MethylScore can identify DMRs from hundreds of samples and how its data-driven approach can stratify associated samples without prior information. We identify DMRs in the A. thaliana 1,001 Genomes dataset to unveil known and unknown genotype-epigenotype associations .

Sequences to Differences in Gene Expression: Analysis of RNA-Seq Data

RNA-Seq is now a routinely employed assay to measure gene expression. As the technique matured over the last decade, so have dedicated analytic tools. In this chapter, we first describe the mainstream as well as the most up-to-date protocols and their implications on downstream analysis. We then detail the steps entailing RNA-Seq analysis in three main stages: (i) preprocessing and data preparation, (ii) upstream processing, and (iii) high-level analyses. We review the most recent and relevant tools as one workflow following a stepwise order. The chapter further encompasses in-depth features of these tools. Details of the required code are made available throughout the chapter, as well as of the underlying statistics. We illustrate these steps with analysis of publicly available RNA-Seq data.

Comprehensive methylation profile of CSF cfDNA revealed pathogenesis and diagnostic markers for early-onset Parkinson's disease

Background: Early-onset Parkinson's disease (EOPD) is one uncommon Parkinson's disease subtype with characteristic clinicopathological features. The full epigenomic profile of EOPD is largely unknown. Methods: We performed the first study to investigate the EOPD full methylation profile of cerebrospinal fluid (CSF) cell-free DNA (cfDNA) from 26 EOPD patients and 10 control patients. Results: 2220 differentially methylated genes were identified in EOPD. Hypermethylation far outweighed hypomethylation in gene numbers. Clustering and enrichment analyses identified aberrant neuronal function and immune response. Weighted correlation network analysis demonstrated significant correlation between methylation signatures and clock drawing test (CDT), mini-mental state examination (MMSE), education, working status, alcohol drinking history and Hamilton anxiety scale (HAMA). Several key networking genes in EOPD aberrant methylation were also identified. Conclusions: The methylation profile and signatures of CSF cfDNA were revealed for the first time in EOPD. Aberrant methylation signatures were correlated with education, working status, alcohol drinking history, CDT, MMSE and HAMA.

Whole-genome methylation analysis of testicular germ cells from cryptozoospermic men points to recurrent and functionally relevant DNA methylation changes

BACKGROUND

Several studies have reported an association between male infertility and aberrant sperm DNA methylation patterns, in particular in imprinted genes. In a recent investigation based on whole methylome and deep bisulfite sequencing, we have not found any evidence for such an association, but have demonstrated that somatic DNA contamination and genetic variation confound methylation studies in sperm of severely oligozoospermic men. To find out whether testicular germ cells (TGCs) of such patients might carry aberrant DNA methylation, we compared the TGC methylomes of four men with cryptozoospermia (CZ) and four men with obstructive azoospermia, who had normal spermatogenesis and served as controls (CTR).

RESULTS

There was no difference in DNA methylation at the whole genome level or at imprinted regions between CZ and CTR samples. However, using stringent filters to identify group-specific methylation differences, we detected 271 differentially methylated regions (DMRs), 238 of which were hypermethylated in CZ (binominal test, p < 2.2 × 10-16). The DMRs were enriched for distal regulatory elements (p = 1.0 × 10-6) and associated with 132 genes, 61 of which are differentially expressed at various stages of spermatogenesis. Almost all of the 67 DMRs associated with the 61 genes (94%) are hypermethylated in CZ (63/67, p = 1.107 × 10-14). As judged by single-cell RNA sequencing, 13 DMR-associated genes, which are mainly expressed during meiosis and spermiogenesis, show a significantly different pattern of expression in CZ patients. In four of these genes, the promoter is hypermethylated in CZ men, which correlates with a lower expression level in these patients. In the other nine genes, eight of which downregulated in CZ, germ cell-specific enhancers may be affected.

CONCLUSIONS

We found that impaired spermatogenesis is associated with DNA methylation changes in testicular germ cells at functionally relevant regions of the genome. We hypothesize that the described DNA methylation changes may reflect or contribute to premature abortion of spermatogenesis and therefore not appear in the mature, motile sperm.