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

Whole-genome bisulfite sequencing data analysis learning module on Google Cloud Platform

This study describes the development of a resource module that is part of a learning platform named 'NIGMS Sandbox for Cloud-based Learning' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox at the beginning of this Supplement. This module is designed to facilitate interactive learning of whole-genome bisulfite sequencing (WGBS) data analysis utilizing cloud-based tools in Google Cloud Platform, such as Cloud Storage, Vertex AI notebooks and Google Batch. WGBS is a powerful technique that can provide comprehensive insights into DNA methylation patterns at single cytosine resolution, essential for understanding epigenetic regulation across the genome. The designed learning module first provides step-by-step tutorials that guide learners through two main stages of WGBS data analysis, preprocessing and the identification of differentially methylated regions. And then, it provides a streamlined workflow and demonstrates how to effectively use it for large datasets given the power of cloud infrastructure. The integration of these interconnected submodules progressively deepens the user's understanding of the WGBS analysis process along with the use of cloud resources. Through this module, we can enhance the accessibility and adoption of cloud computing in epigenomic research, speeding up the advancements in the related field and beyond. This manuscript describes the development of a resource module that is part of a learning platform named ``NIGMS Sandbox for Cloud-based Learning'' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox [1] at the beginning of this Supplement. This module delivers learning materials on the analysis of bulk and single-cell ATAC-seq data in an interactive format that uses appropriate cloud resources for data access and analyses.

BSXplorer: analytical framework for exploratory analysis of BS-seq data

BACKGROUND

Bisulfite sequencing detects and quantifies DNA methylation patterns, contributing to our understanding of gene expression regulation, genome stability maintenance, conservation of epigenetic mechanisms across divergent taxa, epigenetic inheritance and, eventually, phenotypic variation. Graphical representation of methylation data is crucial in exploring epigenetic regulation on a genome-wide scale in both plants and animals. This is especially relevant for non-model organisms with poorly annotated genomes and/or organisms where genome sequences are not yet assembled on chromosome level. Despite being a technology of choice to profile DNA methylation for many years now there are surprisingly few lightweight and robust standalone tools available for efficient graphical analysis of data in non-model systems. This significantly limits evolutionary studies and agrigenomics research. BSXplorer is a tool specifically developed to fill this gap and assist researchers in explorative data analysis and in visualising and interpreting bisulfite sequencing data more easily.

RESULTS

BSXplorer provides in-depth graphical analysis of sequencing data encompassing (a) profiling of methylation levels in metagenes or in user-defined regions using line plots and heatmaps, generation of summary statistics charts, (b) enabling comparative analyses of methylation patterns across experimental samples, methylation contexts and species, and (c) identification of modules sharing similar methylation signatures at functional genomic elements. The tool processes methylation data quickly and offers API and CLI capabilities, along with the ability to create high-quality figures suitable for publication.

CONCLUSIONS

BSXplorer facilitates efficient methylation data mining, contrasting and visualization, making it an easy-to-use package that is highly useful for epigenetic research.

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.

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.

Phenotypic maps for precision medicine: a promising systems biology tool for assessing therapy response and resistance at a personalized level

In this perspective we discuss how tumor heterogeneity and therapy resistance necessitate a focus on more personalized approaches, prompting a shift toward precision medicine. At the heart of the shift towards personalized medicine, omics-driven systems biology becomes a driving force as it leverages high-throughput technologies and novel bioinformatics tools. These enable the creation of systems-based maps, providing a comprehensive view of individual tumor's functional plasticity. We highlight the innovative PHENOSTAMP program, which leverages high-dimensional data to construct a visually intuitive and user-friendly map. This map was created to encapsulate complex transitional states in cancer cells, such as Epithelial-Mesenchymal Transition (EMT) and Mesenchymal-Epithelial Transition (MET), offering a visually intuitive way to understand disease progression and therapeutic responses at single-cell resolution in relation to EMT-related single-cell phenotypes. Most importantly, PHENOSTAMP functions as a reference map, which allows researchers and clinicians to assess one clinical specimen at a time in relation to their phenotypic heterogeneity, setting the foundation on constructing phenotypic maps for personalized medicine. This perspective argues that such dynamic predictive maps could also catalyze the development of personalized cancer treatment. They hold the potential to transform our understanding of cancer biology, providing a foundation for a future where therapy is tailored to each patient's unique molecular and cellular tumor profile. As our knowledge of cancer expands, these maps can be continually refined, ensuring they remain a valuable tool in precision oncology.

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.