5-Hydroxymethylcytosine signatures in circulating cell-free DNA as diagnostic and predictive biomarkers for coronary artery disease

An update on the cell-free DNA-derived methylome as a non-invasive biomarker for coronary artery disease

Cardiovascular diseases are the foremost contributor to global mortality, presenting a complex etiology and an expanding array of risk factors. Coronary artery disease characterized by atherosclerotic plaque build-up in the coronary arteries, imposes significant mortality and financial burdens, especially in low- and middle-income nations. The pathogenesis of coronary artery disease involves a multifaceted interplay of genetic, environmental, and epigenetic factors. Epigenetic regulation contributes to the dynamic control of gene expression without altering the underlying DNA sequence. The mounting evidence that highlights the pivotal role of epigenetic regulation in coronary artery disease development and progression, offering potential avenues for the development of novel diagnostic biomarkers and therapeutic targets. Abnormal DNA methylation patterns are linked to the modulation of gene expression involved in crucial processes like lipid metabolism, inflammation, and vascular function in the context of coronary artery disease. Cell-free DNA has become invaluable in tumor biology as a liquid biopsy, while its applications in coronary artery disease are limited, but intriguing. Atherosclerotic plaque rupture causes myocardial infarction, by depriving heart muscles of oxygen, releasing cell-free DNA from dead cardiac cells, and providing a minimally invasive source to explore tissue-specific epigenetic alterations. We discussed the methodologies for studying the global methylome and hydroxy-methylome landscape, their advantages, and limitations. It explores methylome alterations in coronary artery disease, considering risk factors and their relevance in coronary artery disease genesis. The review also details the implications of MI-derived cell-free DNA for developing minimally invasive biomarkers and associated challenges.

Differential gene expression patterns in ST-elevation Myocardial Infarction and Non-ST-elevation Myocardial Infarction

The ST-elevation Myocardial Infarction (STEMI) and Non-ST-elevation Myocardial Infarction (NSTEMI) might occur because of coronary artery stenosis. The gene biomarkers apply to the clinical diagnosis and therapeutic decisions in Myocardial Infarction. The aim of this study was to introduce, enrich and estimate timely the blood gene profiles based on the high-throughput data for the molecular distinction of STEMI and NSTEMI. The text mining data (50 genes) annotated with DisGeNET data (144 genes) were merged with the GEO gene expression data (5 datasets) using R software. Then, the STEMI and NSTEMI networks were primarily created using the STRING server, and improved using the Cytoscape software. The high-score genes were enriched using the KEGG signaling pathways and Gene Ontology (GO). Furthermore, the genes were categorized to determine the NSTEMI and STEMI gene profiles. The time cut-off points were identified statistically by monitoring the gene profiles up to 30 days after Myocardial Infarction (MI). The gene heatmaps were clearly created for the STEMI (high-fold genes 69, low-fold genes 45) and NSTEMI (high-fold genes 68, low-fold genes 36). The STEMI and NSTEMI networks suggested the high-score gene profiles. Furthermore, the gene enrichment suggested the different biological conditions for STEMI and NSTEMI. The time cut-off points for the NSTEMI (4 genes) and STEMI (13 genes) gene profiles were established up to three days after Myocardial Infarction. The study showed the different pathophysiologic conditions for STEMI and NSTEMI. Furthermore, the high-score gene profiles are suggested to measure up to 3 days after MI to distinguish the STEMI and NSTEMI.

Analysis of genome-wide 5-hydroxymethylation of blood samples stored in different anticoagulants: opportunities for the expansion of clinical resources for epigenetic research

BACKGROUND

Elucidating epigenetic mechanisms could provide new biomarkers for disease diagnosis and prognosis. Technological advances allow genome-wide profiling of 5-hydroxymethylcytosines (5hmC) in liquid biopsies. 5hmC-Seal followed by NGS is a highly sensitive technique for 5hmC biomarker discovery in cfDNA. Currently, 5hmC Seal is optimized for EDTA blood collection. We asked whether heparin was compatible with 5hmC Seal as many clinical and biobanked samples are stored in heparin.

METHODS

We obtained 60 samples in EDTA matched to 60 samples in heparin from the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial. Samples were comprised of 30 controls and 30 individuals who were later diagnosed with colon cancer. We profiled genome-wide 5hmC in cfDNA using 5hmC-Seal assay followed by NGS. The 5hmC profiling data from samples collected in EDTA were systematically compared to those in heparin across various genomic features.

RESULTS

cfDNA isolation and library construction appeared comparable in heparin vs. EDTA. Typical genomic distribution patterns of 5hmC, including gene bodies and enhancer markers, were comparable in heparin vs. EDTA. 5hmC analysis of cases and controls yielded highly correlated differential features suggesting that both anticoagulants were compatible with 5hmC Seal assay.

CONCLUSIONS

While not currently recommended for the 5hmC-Seal protocol, blood samples stored in heparin were successfully used to generate analysable and biologically relevant genome-wide 5hmC profiling. Our findings are the first to support opportunities to expand the biospecimen resource to heparin samples for 5hmC Seal and perhaps other PCR-based technologies in epigenetic research.

TET (Ten-eleven translocation) family proteins: structure, biological functions and applications

Ten-eleven translocation (TET) family proteins (TETs), specifically, TET1, TET2 and TET3, can modify DNA by oxidizing 5-methylcytosine (5mC) iteratively to yield 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxycytosine (5caC), and then two of these intermediates (5fC and 5caC) can be excised and return to unmethylated cytosines by thymine-DNA glycosylase (TDG)-mediated base excision repair. Because DNA methylation and demethylation play an important role in numerous biological processes, including zygote formation, embryogenesis, spatial learning and immune homeostasis, the regulation of TETs functions is complicated, and dysregulation of their functions is implicated in many diseases such as myeloid malignancies. In addition, recent studies have demonstrated that TET2 is able to catalyze the hydroxymethylation of RNA to perform post-transcriptional regulation. Notably, catalytic-independent functions of TETs in certain biological contexts have been identified, further highlighting their multifunctional roles. Interestingly, by reactivating the expression of selected target genes, accumulated evidences support the potential therapeutic use of TETs-based DNA methylation editing tools in disorders associated with epigenetic silencing. In this review, we summarize recent key findings in TETs functions, activity regulators at various levels, technological advances in the detection of 5hmC, the main TETs oxidative product, and TETs emerging applications in epigenetic editing. Furthermore, we discuss existing challenges and future directions in this field.

PETCH-DB: a Portal for Exploring Tissue-specific and Complex disease-associated 5-Hydroxymethylcytosines

Epigenetic modifications play critical roles in gene regulation and disease pathobiology. Highly sensitive enabling technologies, including microarray- and sequencing-based approaches have allowed genome-wide profiling of cytosine modifications in DNAs in clinical samples to facilitate discovery of epigenetic biomarkers for disease diagnosis and prognosis. Historically, many previous studies, however, did not distinguish the most investigated 5-methylcytosines (5mC) from other modified cytosines, especially the biochemically stable 5-hydroxymethylcytosines (5hmC), which have been shown to have a distinct genomic distribution and regulatory role from 5mC. Notably, during the past several years, the 5hmC-Seal, a highly sensitive chemical labeling technique, has been demonstrated to be a powerful tool for genome-wide profiling of 5hmC in clinically feasible biospecimens (e.g. a few milliliter of plasma or serum). The 5hmC-Seal technique has been utilized by our team in biomarker discovery for human cancers and other complex diseases using circulating cell-free DNA (cfDNA), as well as the characterization of the first 5hmC Human Tissue Map. Convenient access to the accumulating 5hmC-Seal data will allow the research community to validate and re-use these results, potentially providing novel insights into epigenetic contribution to a range of human diseases. Here we introduce the PETCH-DB, an integrated database that was implemented to provide 5hmC-related results generated using the 5hmC-Seal technique. We aim the PETCH-DB to be a central portal, which will be available to the scientific community with regularly updated 5hmC data in clinical samples to reflect current advances in this field. Database URL http://petch-db.org/.

Obesity and dyslipidemia are associated with partially reversible modifications to DNA hydroxymethylation of apoptosis- and senescence-related genes in swine adipose-derived mesenchymal stem/stromal cells

BACKGROUND

Obesity dysregulates key biological processes underlying the functional homeostasis, fate decisions, and reparative potential of mesenchymal stem/stromal cells (MSCs). Mechanisms directing obesity-induced phenotypic alterations in MSCs remain unclear, but emerging drivers include dynamic modification of epigenetic marks, like 5-hydroxymethylcytosine (5hmC). We hypothesized that obesity and cardiovascular risk factors induce functionally relevant, locus-specific changes in 5hmC of swine adipose-derived MSCs and evaluated their reversibility using an epigenetic modulator, vitamin-C.

METHODS

Female domestic pigs were fed a 16-week Lean or Obese diet (n = 6 each). MSCs were harvested from subcutaneous adipose tissue, and 5hmC profiles were examined through hydroxymethylated DNA immunoprecipitation sequencing (hMeDIP-seq) followed by an integrative (hMeDIP and mRNA sequencing) gene set enrichment analysis. For clinical context, we compared 5hmC profiles of adipose tissue-derived human MSCs harvested from patients with obesity and healthy controls.

RESULTS

hMeDIP-seq revealed 467 hyper- (fold change ≥ 1.4; p-value ≤ 0.05) and 591 hypo- (fold change ≤ 0.7; p-value ≤ 0.05) hydroxymethylated loci in swine Obese- versus Lean-MSCs. Integrative hMeDIP-seq/mRNA-seq analysis identified overlapping dysregulated gene sets and discrete differentially hydroxymethylated loci with functions related to apoptosis, cell proliferation, and senescence. These 5hmC changes were associated with increased senescence in cultured MSCs (p16/CDKN2A immunoreactivity, senescence-associated β-galactosidase [SA-β-Gal] staining), were partly reversed in swine Obese-MSCs treated with vitamin-C, and shared common pathways with 5hmC changes in human Obese-MSCs.

CONCLUSIONS

Obesity and dyslipidemia are associated with dysregulated DNA hydroxymethylation of apoptosis- and senescence-related genes in swine and human MSCs, potentially affecting cell vitality and regenerative functions. Vitamin-C may mediate reprogramming of this altered epigenomic landscape, providing a potential strategy to improve the success of autologous MSC transplantation in obese patients.

Regional gain and global loss of 5-hydroxymethylcytosine coexist in genitourinary cancers and regulate different oncogenic pathways

BACKGROUND

DNA 5-hydroxymethylcytosine (5hmC) is produced by dynamic 5mC oxidation process contributing to tissue specification, and loss of 5hmC has been reported in multiple cancers including genitourinary cancers. However, 5hmC is also cell-type specific, and its variability may exist between differentiated tumor cells and cancer stem cells. Thus, cancer-associated changes in 5hmC may be contributed by distinct sets of tumor cells within the tumor tissues.

RESULTS

Here, we applied a sensitive immunoprecipitation-based method (hMeDIP-seq) to analyze 5hmC changes during genitourinary carcinogenesis (including prostate, urothelial and kidney). We confirmed the tissue-specific distribution of 5hmC in genitourinary tissues and identified regional gain and global loss of 5hmC coexisting in genitourinary cancers. The genes with gain of 5hmC during tumorigenesis were functionally enriched in regulating stemness and hypoxia, whereas were associated with poor clinical prognosis irrespective of their differences in tumor type. We identified that gain of 5hmC occurred in soft fibrin gel-induced 3D tumor spheres with a tumor-repopulating phenotype in two prostate cancer cell lines, 22RV1 and PC3, compared with conventional two-dimensional (2D) rigid dishes. Then, we defined a malignant signature derived from the differentially hydroxymethylated regions affected genes of cancer stem-like cells, which could predict a worse clinical outcome and identified phenotypically malignant populations of cells from prostate cancer tumors. Notably, an oxidation-resistant vitamin C derivative, ascorbyl phosphate magnesium, restored 5hmC and killed the cancer stem cell-like cells leading to apoptosis in prostate cancer cell lines.

CONCLUSIONS

Collectively, our study dissects the regional gain of 5hmC in maintaining cancer stem-like cells and related to poor prognosis, which provides proof of concept for an epigenetic differentiation therapy with vitamin C by 5hmC reprogramming.

Association Between Circulating Cell-Free DNA Level at Admission and the Risk of Heart Failure Incidence in Acute Myocardial Infarction Patients

Plasma cell-free DNA (cfDNA) was elevated in patients with acute myocardial infarction (AMI) or heart failure (HF). However, whether cfDNA could serve as a predictor for risk of HF after AMI remains unknown. In this study, we conducted a pilot prospective cohort study in which 98 AMI patients were enrolled from a single center to assess the association between cfDNA levels at admission and risk of HF in an AMI population. Patients with cfDNA above the median level (14.39 ng/mL) showed higher low-density lipoprotein cholesterol, cardiac troponin I (cTnI), and soluble suppression of tumorigenicity 2 (sST2) levels compared with patients below the median. cfDNA was positively correlated with cTnI (r = 0.377, p < 0.001) and sST2 (r = 0.443, p < 0.001). Within a median follow-up of about 345 days, 46 patients (52.6%) developed HF. Multivariate Cox analysis showed that a higher cfDNA (above the cutoff value: 9.227 ng/mL) was an effective risk predictor (C-index = 0.74, 95% confidence interval [CI]: 0.733-0.748) for HF incidence after AMI (adjusted hazard ratio [HR]: 2.805; 95% CI: 1.087-7.242; p = 0.033). Moreover, a linear association was observed between cfDNA and risk of HF incidence adjusted for by age, gender, and history of chronic kidney disease (p for linear trend = 0.044). Taken together, the cfDNA levels at admission are associated with the incidence of HF in AMI patients. A positive correlation between cfDNA and the fibrotic factor sST2 was proved, but the underlying mechanisms require further study.

5-Hydroxymethylcytosine Signatures in Circulating Cell-Free DNA as Early Warning Biomarkers for COVID-19 Progression and Myocardial Injury

Background: The symptoms of coronavirus disease 2019 (COVID-19) range from moderate to critical conditions, leading to death in some patients, and the early warning indicators of the COVID-19 progression and the occurrence of its serious complications such as myocardial injury are limited.

Methods: We carried out a multi-center, prospective cohort study in three hospitals in Wuhan. Genome-wide 5-hydroxymethylcytosine (5hmC) profiles in plasma cell-free DNA (cfDNA) was used to identify risk factors for COVID-19 pneumonia and develop a machine learning model using samples from 53 healthy volunteers, 66 patients with moderate COVID-19, 99 patients with severe COVID-19, and 38 patients with critical COVID-19.

Results: Our warning model demonstrated that an area under the curve (AUC) for 5hmC warning moderate patients developed into severe status was 0.81 (95% CI 0.77–0.85) and for severe patients developed into critical status was 0.92 (95% CI 0.89–0.96). We further built a warning model on patients with and without myocardial injury with the AUC of 0.89 (95% CI 0.84–0.95).

Conclusion: This is the first study showing the utility of 5hmC as an accurate early warning marker for disease progression and myocardial injury in patients with COVID-19. Our results show that phosphodiesterase 4D and ten-eleven translocation 2 may be important markers in the progression of COVID-19 disease.

Computational challenges in detection of cancer using cell-free DNA methylation

Cell-free DNA(cfDNA) methylation profiling is considered promising and potentially reliable for liquid biopsy to study progress of diseases and develop reliable and consistent diagnostic and prognostic biomarkers. There are several different mechanisms responsible for the release of cfDNA in blood plasma, and henceforth it can provide information regarding dynamic changes in the human body. Due to the fragmented nature, low concentration of cfDNA, and high background noise, there are several challenges in its analysis for regular use in diagnosis of cancer. Such challenges in the analysis of the methylation profile of cfDNA are further aggravated due to heterogeneity, biomarker sensitivity, platform biases, and batch effects. This review delineates the origin of cfDNA methylation, its profiling, and associated computational problems in analysis for diagnosis. Here we also contemplate upon the multi-marker approach to handle the scenario of cancer heterogeneity and explore the utility of markers for 5hmC based cfDNA methylation pattern. Further, we provide a critical overview of deconvolution and machine learning methods for cfDNA methylation analysis. Our review of current methods reveals the potential for further improvement in analysis strategies for detecting early cancer using cfDNA methylation.

Research progress in predicting DNA methylation modifications and the relation with human diseases

DNA methylation is an important mode of regulation in epigenetic mechanisms, and it is one of the research foci in the field of epigenetics. DNA methylation modification affects a series of biological processes, such as eukaryotic cell growth, differentiation and transformation mechanisms, by regulating gene expression. In this review, we systematically summarized the DNA methylation databases, prediction tools for DNA methylation modification, machine learning algorithms for predicting DNA methylation modification, and the relationship between DNA methylation modification and diseases such as hypertension, Alzheimer's disease, diabetic nephropathy, and cancer. An in-depth understanding of DNA methylation mechanisms can promote accurate prediction of DNA methylation modifications and the treatment and diagnosis of related diseases.

Whole-Genome Differentially Hydroxymethylated DNA Regions among Twins Discordant for Cardiovascular Death

Epigenetics is a mechanism underlying cardiovascular disease. It is unknown whether DNA hydroxymethylation is prospectively associated with the risk for cardiovascular death independent of germline and common environment. Male twin pairs middle-aged in 1969–1973 and discordant for cardiovascular death through December 31, 2014, were included. Hydroxymethylation was quantified in buffy coat DNA collected in 1986–1987. The 1893 differentially hydroxymethylated regions (DhMRs) were identified after controlling for blood leukocyte subtypes and age among 12 monozygotic (MZ) pairs (Benjamini–Hochberg False Discovery Rate < 0.01), of which the 102 DhMRs were confirmed with directionally consistent log₂-fold changes and p < 0.01 among additional 7 MZ pairs. These signature 102 DhMRs, independent of the germline, were located on all chromosomes except for chromosome 21 and the Y chromosome, mainly within/overlapped with intergenic regions and introns, and predominantly hyper-hydroxymethylated. A binary linear classifier predicting cardiovascular death among 19 dizygotic pairs was identified and equivalent to that generated from MZ via the 2D transformation. Computational bioinformatics discovered pathways, phenotypes, and DNA motifs for these DhMRs or their subtypes, suggesting that hydroxymethylation was a pathophysiological mechanism underlying cardiovascular death that might be influenced by genetic factors and warranted further investigations of mechanisms of these signature regions in vivo and in vitro.

Cancer Biomarkers Discovery of Methylation Modification With Direct High-Throughput Nanopore Sequencing

Cancer is a complex disease, driven by a combination of genetic and epigenetic alterations. DNA and RNA methylation modifications are the most common epigenetic events that play critical roles in cancer development and progression. Bisulfite converted sequencing is a widely used technique to detect base modifications in DNA methylation, but its main drawbacks lie in DNA degradation, lack of specificity, or short reads with low sequence diversity. The nanopore sequencing technology can directly detect base modifications in native DNA as well as RNA without harsh chemical treatment, compared to bisulfite sequencing. Furthermore, CRISPR/Cas9-targeted enrichment nanopore sequencing techniques are straightforward and cost-effective when targeting genomic regions are of interest. In this review, we mainly focus on DNA and RNA methylation modification detection in cancer with the current nanopore sequencing approaches. We also present the respective strengths, weaknesses of nanopore sequencing techniques, and their future translational applications in identification of epigenetic biomarkers for cancer detection and prognosis.

A human tissue map of 5-hydroxymethylcytosines exhibits tissue specificity through gene and enhancer modulation

DNA 5-hydroxymethylcytosine (5hmC) modification is known to be associated with gene transcription and frequently used as a mark to investigate dynamic DNA methylation conversion during mammalian development and in human diseases. However, the lack of genome-wide 5hmC profiles in different human tissue types impedes drawing generalized conclusions about how 5hmC is implicated in transcription activity and tissue specificity. To meet this need, we describe the development of a 5hmC tissue map by characterizing the genomic distributions of 5hmC in 19 human tissues derived from ten organ systems. Subsequent sequencing results enabled the identification of genome-wide 5hmC distributions that uniquely separates samples by tissue type. Further comparison of the 5hmC profiles with transcriptomes and histone modifications revealed that 5hmC is preferentially enriched on tissue-specific gene bodies and enhancers. Taken together, the results provide an extensive 5hmC map across diverse human tissue types that suggests a potential role of 5hmC in tissue-specific development; as well as a resource to facilitate future studies of DNA demethylation in pathogenesis and the development of 5hmC as biomarkers.

5-Hydroxymethylcytosine signature in circulating cell-free DNA as a potential diagnostic factor for early-stage colorectal cancer and precancerous adenoma

Approximately 85% colorectal cancers (CRCs) are thought to evolve through the adenoma‐to‐carcinoma sequence associated with specific molecular alterations, including the 5‐hydroxymethylcytosine (5hmC) signature in circulating cell‐free DNA (cfDNA). To explore colorectal disease progression and evaluate the use of cfDNA as a potential diagnostic factor for CRC screening, here, we performed genome‐wide 5hmC profiling in plasma cfDNA and tissue genomic DNA (gDNA) acquired from 101 samples (63 plasma and 38 tissues), collected from 21 early‐stage CRC patients, 21 AD patients, and 21 healthy controls (HC). The gDNA and cfDNA 5hmC signatures identified in gene bodies and promoter regions in CRC and AD groups were compared with those in HC group. All the differential 5hmC‐modified regions (DhMRs) were gathered into four clusters: Disease‐enriched, AD‐enriched, Disease‐lost, and AD‐lost, with no overlap. AD‐related clusters, AD‐enriched and AD‐lost, displayed the unique 5hmC signals in AD patients. Disease‐enriched and Disease‐lost clusters indicated the general 5hmC changes when colorectal lesions occurred. Cancer patients with a confirmable adenoma history segmentally gathered in AD‐enriched clusters. KEGG functional enrichment and GO analyses determined distinct differential 5hmC‐modified profiles in cfDNA of HC individuals, AD, and CRC patients. All patients had comprehensive 5hmC signatures where Disease‐enriched and Disease‐lost DhMR clusters demonstrated similar epigenetic modifications, while AD‐enriched and AD‐lost DhMR clusters indicated complicated subpopulations in adenoma. Analysis of CRC patients with adenoma history showed exclusive 5hmC‐gain characteristics, consistent with the ‘parallel’ evolution hypothesis in adenoma, either developed through the adenoma‐to‐carcinoma sequence or not. These findings deepen our understanding of colorectal disease and suggest that the 5hmC modifications of different pathological subtypes (cancer patients with or without adenoma history) could be used to screen early‐stage CRC and assess adenoma malignancy with large‐scale follow‐up studies in the future.