Extraction and high-throughput sequencing of oak heartwood DNA: Assessing the feasibility of genome-wide DNA methylation profiling

Tree ring features are affected by environmental factors and therefore are the basis for dendrochronological studies to reconstruct past environmental conditions. Oak wood often provides the data for these studies because of the durability of oak heartwood and hence the availability of samples spanning long time periods of the distant past. Wood formation is regulated in part by epigenetic mechanisms such as DNA methylation. Studies of the methylation state of DNA preserved in oak heartwood thus could identify epigenetic tree ring features informing on past environmental conditions. In this study, we aimed to establish protocols for the extraction of DNA, the high-throughput sequencing of whole-genome DNA libraries (WGS) and the profiling of DNA methylation by whole-genome bisulfite sequencing (WGBS) for oak (Quercus robur) heartwood drill cores taken from the trunks of living standing trees spanning the AD 1776-2014 time period. Heartwood contains little DNA, and large amounts of phenolic compounds known to hinder the preparation of high-throughput sequencing libraries. Whole-genome and DNA methylome library preparation and sequencing consistently failed for oak heartwood samples more than 100 and 50 years of age, respectively. DNA fragmentation increased with sample age and was exacerbated by the additional bisulfite treatment step during methylome library preparation. Relative coverage of the non-repetitive portion of the oak genome was sparse. These results suggest that quantitative methylome studies of oak hardwood will likely be limited to relatively recent samples and will require a high sequencing depth to achieve sufficient genome coverage.

Validation of SFRP1 Promoter Hypermethylation in Plasma as a Prognostic Marker for Survival and Gemcitabine Effectiveness in Patients with Stage IV Pancreatic Adenocarcinoma

No reliable predictive blood-based biomarkers are available for determining survival from pancreatic adenocarcinoma (PDAC). This combined discovery and validation study examines promoter hypermethylation (ph) of secreted frizzled-related protein 1 (SFRP1) in plasma-derived cell-free DNA as an independent prognostic marker for survival and Gemcitabine effectiveness in patients with stage IV PDAC. We conducted methylation-specific polymerase chain reaction analysis of the promoter region of the SFRP1 gene, based on bisulfite treatment. Survival was analyzed with Kaplan–Meier curves, log-rank test, and Cox regression. The discovery cohort included 40 patients, 25 receiving Gem. Gem-treated patients with phSFRP1 had a shorter median overall survival (mOS) (4.4 months) than unmethylated patients (11.6 months). Adjusted Cox-regression yielded a hazard rate (HR) of 3.48 (1.39–8.70). The validation cohort included 58 Gem-treated patients. Patients with phSFRP1 had a shorter mOS (3.2 months) than unmethylated patients (6.3 months). Adjusted Cox regression yielded an HR of 3.53 (1.85–6.74). In both cohorts, phSFRP1 was associated with poorer survival in Gem-treated patients. This may indicate that tumors with phSFRP1 are more aggressive and less sensitive to Gem treatment. This knowledge may facilitate tailored treatment of patients with stage IV PDAC. Further studies are planned to examine phSFRP1 in more intensive chemotherapy regimens.

Detecting cell-of-origin and cancer-specific features of cell-free DNA with Nanopore sequencing

DNA methylation (5mC) is a promising biomarker for detecting circulating tumor DNA (ctDNA), providing information on a cell's genomic regulation, developmental lineage, and molecular age. Sequencing assays for detecting ctDNA methylation involve pre-processing steps such as immunoprecipitation, enzymatic treatment, or the most common method, sodium bisulfite treatment. These steps add complexity and time that pose a challenge for clinical labs, and bisulfite treatment in particular degrades input DNA and can result in loss of informative ctDNA fragmentation patterns. In this feasibility study, we demonstrate that whole genome sequencing of circulating cell-free DNA using conventional Oxford Nanopore Technologies (ONT) sequencing can accurately detect cell-of-origin and cancer-specific 5mC changes while preserving important fragmentomic information. The simplicity of this approach makes it attractive as a liquid biopsy assay for cancer as well as non-cancer applications in emergency medicine.

Extraction and high-throughput sequencing of oak heartwood DNA: assessing the feasibility of genome-wide DNA methylation profiling

Tree ring features are affected by environmental factors and therefore are the basis for dendrochronological studies to reconstruct past environmental conditions. Oak wood often provides the data for these studies because of the durability of oak heartwood and hence the availability of samples spanning long time periods of the distant past. Wood formation is regulated in part by epigenetic mechanisms such as DNA methylation. Studies of the methylation state of DNA preserved in oak heartwood thus could identify epigenetic tree ring features informing on past environmental conditions. In this study, we aimed to establish protocols for the extraction of DNA, the high-throughput sequencing of whole-genome DNA libraries (WGS) and the profiling of DNA methylation by whole-genome bisulfite sequencing (WGBS) for oak (\textit{Quercus spp.}) heartwood drill cores taken from the trunks of living standing trees spanning the AD 1776-2014 time period. Heartwood contains little DNA, and large amounts of phenolic compounds known to hinder the preparation of high-throughput sequencing libraries. Whole-genome and DNA methylome library preparation and sequencing consistently failed for oak heartwood samples more than 100 and 50 years of age, respectively. DNA fragmentation increased with sample age and was exacerbated by the additional bisulfite treatment step during methylome library preparation. Relative coverage of the non-repetitive portion of the oak genome was sparse. These results suggest that quantitative methylome studies of oak hardwood will likely be limited to relatively recent samples and will require a high sequencing depth to achieve sufficient genome coverage.

Cardiomyocyte-Specific Circulating Cell-Free Methylated DNA in Esophageal Cancer Patients Treated with Chemoradiation

Thoracic high-dose radiation therapy (RT) for cancer has been associated with early and late cardiac toxicity. To assess altered rates of cardiomyocyte cell death due to RT we monitored changes in cardiomyocyte-specific, cell-free methylated DNA (cfDNA) shed into the circulation. Eleven patients with distal esophageal cancer treated with neoadjuvant chemoradiation to 50.4 Gy (RT) and concurrent carboplatin and paclitaxel were enrolled. Subjects underwent fasting blood draws prior to the initiation and after completion of RT as well as 4–6 months following RT. An island of six unmethylated CpGs in the FAM101A locus was used to identify cardiomyocyte-specific cfDNA in serum. After bisulfite treatment this specific cfDNA was quantified by amplicon sequencing at a depth of >35,000 reads/molecule. Cardiomyocyte-specific cfDNA was detectable before RT in the majority of patient samples and showed some distinct changes during the course of treatment and recovery. We propose that patient-specific cardiac damages in response to the treatment are indicated by these changes although co-morbidities may obscure treatment-specific events.

Profiling a single-stranded DNA region within an rDNA segment that is a loading site for bacterial condensin

Bacterial condensin preferentially loads to single-stranded DNA (ssDNA) in vitro and loads onto rDNA in vivo to support proper chromosome compaction. Thus, the actively transcribing rDNA would provide the ssDNA region for the topological loading of bacterial condensin. We attempted to detect the ssDNA region in the rrnI gene in situ. Non-denaturing sodium bisulfite treatment catalyzed the conversion of cytosines to thymines via uracils (CT-conversion) at locally melted DNA of a bacterial genome. Using next-generation sequencing, we generated an average of 11,000 reads covering each cytosine on the PCR-amplified rDNA segment to obtain the actual CT-conversion rate. In principle, the CT-conversion rate is an accurate guide to detect the formation of the ssDNA segment. We expected that an increment of the CT-conversion rate would reflect a trend toward ssDNA accumulation at a given site within the rDNA. We detected multiple ssDNA segments throughout the rDNA. The deletion mutations of the rDNA that affect the bacterial-condensin loading hindered the ssDNA formation only at the 100-500 bp segment downstream of the promoter. These data support the idea that the ssDNA segment plays a crucial role as the bacterial condensin-loading site and suggest the mechanism of condensin loading onto rDNA.

An effective method to resolve ambiguous bisulfite-treated reads

Background The combination of the bisulfite treatment and the next-generation sequencing is an important method for methylation analysis, and aligning the bisulfite-treated reads (BS-reads) is the critical step for the downstream applications. As bisulfite treatment reduces the complexity of the sequences, a large portion of BS-reads might be aligned to multiple locations of the reference genome ambiguously, called multireads. These multireads cannot be employed in the downstream applications since they are likely to introduce artifacts. To identify the best mapping location of each multiread, existing Bayesian-based methods calculate the probability of the read at each position by considering how does it overlap with unique mapped reads. However, $$\sim 25$$ ∼ 25 % of multireads are not overlapped with any unique reads, which are unresolvable for existing method. Results Here we propose a novel method (EM-MUL) that not only rescues multireads overlapped with unique reads, but also uses the overall coverage and accurate base-level alignment to resolve multireads that cannot be handled by current methods. We benchmark our method on both simulated datasets and real datasets. Experimental results show that it is able to align more than 80% of multireads to the best mapping position with very high accuracy. Conclusions EM-MUL is an effective method designed to accurately determine the best mapping position of multireads in BS-reads. For the downstream applications, it is useful to improve the methylation resolution on the repetitive regions of genome. EM-MUL is free available at https://github.com/lmylynn/EM-MUL.

Current Advances in DNA Methylation Analysis Methods

DNA methylation is one of the epigenetic changes, which plays a major role in regulating gene expression and, thus, many biological processes and diseases. There are several methods for determining the methylation of DNA samples. However, selecting the most appropriate method for answering biological questions appears to be a challenging task. The primary methods in DNA methylation focused on identifying the state of methylation of the examined genes and determining the total amount of 5-methyl cytosine. The study of DNA methylation at a large scale of genomic levels became possible following the use of microarray hybridization technology. The new generation of sequencing platforms now allows the preparation of genomic maps of DNA methylation at the single-open level. This review includes the majority of methods available to date, introducing the most widely used methods, the bisulfite treatment, biological identification, and chemical cutting along with their advantages and disadvantages. The techniques are then scrutinized according to their robustness, high throughput capabilities, and cost.

Whole genome 6-methyladenosine sequencing (6-mA-Seq) enables bacterial epigenomics studies

Methylation sequencing using bisulfite treatment has revolutionized the field of molecular biology for eukaryotic systems, unveiling levels of intricacy in regulation of gene expression in response to different environmental conditions. While bacteria also utilize methylation to regulate gene expression, bisulfite sequencing does not work as well in bacteria as in eukaryotes, because bacteria methylate adenosine instead of cytosine. Therefore, global bacterial methylation patterns cannot be studied using common Illumina sequencers. In this work, we demonstrate 6mA-seq, a method that can be used to identify patterns in bacteria that methylate adenosines at GATC sites. Furthermore, this method was used on Escherichia coli cultured on four different carbon sources to demonstrate different methylation patterns due to carbon utilization. 6mA-seq can increase the speed in which epigenetic research is conducted in bacteria.

Analysis of mitochondrial genome methylation using Nanopore single-molecule sequencing

The level and the biological significance of mitochondrial DNA (mtDNA) methylation in human cells is a controversial topic. Using long-read third-generation sequencing technology, mtDNA methylation can be detected directly from the sequencing data, which overcomes previously suggested biases, introduced by bisulfite treatment-dependent methods. We investigated mtDNA from whole blood-derived DNA and established a workflow to detect CpG methylation with Nanopolish. In order to obtain native mtDNA, we adjusted a whole-genome sequencing protocol and performed ligation library preparation and Nanopore sequencing. To validate the workflow, 897bp of methylated and unmethylated synthetic DNA samples at different dilution ratios were sequenced and CpG methylation was detected. Interestingly, we observed that reads with higher methylation in the synthetic DNA did not pass Guppy calling, possibly affecting conclusions about DNA methylation in Nanopore sequencing. We detected in all blood-derived samples overall low-level methylation across the mitochondrial genome, with exceptions at certain CpG sites. Our results suggest that Nanopore sequencing is capable of detecting low-level mtDNA methylation. However, further refinement of the bioinformatical pipelines including Guppy failed reads are recommended.